Mounting scientific evidence supporting the health benefits of nature exposure is likely to increase the supply of high-quality green spaces in urban environments and promote their active use. This might be one way to improve health and welfare in cities.

In this cross-sectional study, higher frequency of green space visits was associated with lower frequency of psychotropic, antihypertensive and asthma medication use, and the association was not dependent on socioeconomic status. Amounts of residential green and blue spaces or green and blue views from home were not associated with psychotropic, antihypertensive and asthma medication use.

Medication use indicates an illness that is serious enough to require the use of a prescribed drug. However, medications have rarely been used as indicators of health status in studies on the health effects of nature exposure. This study aimed to investigate whether the amounts of residential green and blue space, frequency of green space visits and green or blue views from home are separately associated with self-reported use of psychotropic (anxiolytics, hypnotics and antidepressants), antihypertensive and asthma medication in Finland. These medicine groups were chosen because mental health problems, insomnia, hypertension and asthma are major disease groups from a public health perspective, and information on the relationship between nature exposure and asthma is lacking. Another aim was to assess potential effect modification by socioeconomic status (SES) because there is some evidence that the health benefits of nature exposure might be most evident among lowest socioeconomic groups—for example, due to the lack of access to other health-promoting resources. 23 24

Despite inconsistent evidence, there are plausible pathways that can explain potential beneficial effects of nature exposure on physical and mental health and well-being, such as increased physical activity, reduced stress, social cohesion and beneficial immunological reactions, 4 11 18 20–22 but these mechanisms and the interlinkages between them are not yet fully understood.

Inconsistent evidence may stem from differences, for example, in study designs; characteristics of green space exposure, such as definitions, quality and accessibility 15 ; cultural and individual characteristics, such as environmental preferences 4 16 and intrinsic motivation to nature visits. 17 Therefore, it is possible that such associations do not exist in some populations. A majority of evidence on the association between green space and health is derived from objectively measured residential green space, and information on the actual active or passive use of green space is often lacking. 5 18 However, tentative evidence suggests that visiting green spaces or even having views over green spaces may be more strongly associated with health and well-being than the mere amount of green space in the residential area. 19

Nature exposure (ie, exposure to green and blue spaces) is thought to be beneficial, especially for mental, cardiovascular and respiratory health. 1–4 For example, many reviews have reported an association between green space exposure and mental 1 5 and cardiovascular 1 6–8 health, mainly based on cross-sectional studies. However, reviews of longitudinal and experimental studies have evaluated the evidence as inconsistent or suggestive. 9–12 In particular, findings for asthma have been inconsistent. 11 In addition to some protective associations, exposure to different types of green spaces has been linked to adverse asthma outcomes in some, 13 but not all, 14 studies. Regarding blue space exposure, a review of mainly cross-sectional studies suggests an association between exposure to blue space and mental well-being, while evidence regarding cardiovascular health is more limited. 2

To assess the consistency and robustness of the findings, various alternative and additional variables were introduced individually into the model 3 for the sensitivity analyses. First, we tested alternative exposure variables: the amounts of green and blue space within a 300 m buffer (instead of a 1 km buffer), all green spaces (including fields and pastures that were not included in the original variable) within 1 km of the respondent’s home and all natural spaces (the sum of all green and blue spaces) within 1 km of the respondent’s home. Second, self-reported doctor-diagnosed depression, hypertension and asthma were used as the alternative outcomes. Third, the following alternative covariates were used: area-level unemployment rate (%, instead of area-level annual mean income), categorised area-level annual mean income (instead of continuous) and categorised NO 2 from road traffic (instead of continuous). Fourth, additional covariates were included in model 3 individually: BMI (as a spline), leisure-time physical activity, pet ownership, living in a detached/single family house and city. Finally, model 3 analyses were restricted to those who had lived longer than 1 year in their current apartment (93% of the respondents).

To explore the potential effect modification by SES, the interaction term between each of the exposure variables and SES indicators was included in model 3 for all the four exposures and three outcomes. For the interaction analyses, education and annual household income were compressed into three categories instead of the original four, and continuous area-level annual mean income and area-level unemployment rate were categorised into tertiles. If the p value for the interaction term was <0.10, a stratified analysis was performed. To produce an overall p value for the interaction term, a logistic regression model (the LOGISTIC procedure in SAS) was used instead of a logistic generalised additive model; thus, age was categorised (25–34, 35–44, 45–54, 55–64, 65–74 and ≥75 years) instead of treating it as a spline.

Respondents with missing values for any of the variables used in model 3 were excluded from all the analyses to produce comparable results from models 1–3. To test the potential effect of exclusions, model 1 was run without excluding all exposure–outcome pairs. The ORs and their CIs remained essentially similar to the results from model 1 with exclusions.

In the logistic regression analyses, the crude model was adjusted for age and sex (model 1), and the second model (model 2) was additionally adjusted for socioeconomic indicators (marital status, employment status, education, annual household income and area-level annual mean income), health behaviour (alcohol use, smoking and physical activity at work) and recreational property during the warm season. Finally, other environmental exposures (traffic noise and NO 2 from road traffic) were added (model 3). The results are presented as OR estimates and their 95% confidence intervals (95% CIs).

Multivariate associations between green and blue space indicators and medication use were modelled using a logistic generalised additive model (GAMPL procedure in SAS). Thin-plate regression splines were used to explore whether the exposure–response relationships were linear. Based on the non-linear curves in these analyses, age and BMI (in the sensitivity analyses) were added to the models as splines.

Residential outdoor levels of NO 2 from road traffic (μg/m 3 ) were based on emission and dispersion modelling. 29 The dispersion model used was the ‘Contaminants in the Air from a Road-Finnish Meteorological Institute’ (CAR-FMI), 30 which considers the chemical reactions causing nitrogen oxides to become oxidised. The highest modelled annual average concentration in 2010–2014 at the nearest outdoor receptor point was used as a proxy for the home outdoor concentration. This was the only indicator available, but it was not expected to have a significant effect on the spatial variability of the estimated exposures. The grid size for NO 2 varied from 25 m, close to the major roads, to 500 m in rural areas.

Residential exposure to traffic noise (L den dB) was estimated based on façade noise maps from the Sito Consulting Company. Calculations for 2017 were carried out in accordance with the Environmental Noise Directive 2002/49/EC using the Common Noise Assessment Methods in Europe (CNOSSOS-EU) and data from 2016.

Variables likely to be associated with exposure and outcome were selected a priori based on published literature 5 6 9 or rational judgement and included as covariates ( online supplemental figure 3 ). Details of age, sex, marital status, education, employment status, annual household income, smoking, alcohol use, physical activity at work and use of recreational property (ie, spending time at holiday houses, secondary residences and summer cottages) during the warm season (weeks per year) were obtained from the questionnaire. The zip area-level annual mean income (€) and unemployment rate (%, used in sensitivity analyses) were acquired from Statistics Finland. Variables that might be mediators of the examined associations (body mass index (BMI) and leisure-time physical activity, 28 and pet ownership) or have less often been examined as potential confounders (living in a detached/single family house and living longer than 1 year in the current apartment) were introduced into model 3 in the sensitivity analyses.

Self-reported use of psychotropic medication (anxiolytics, hypnotics and antidepressants), antihypertensives and asthma medications was used as outcomes. Other medicine groups enquired in the questionnaire were analgesics for different purposes, hypnotics, hay fever medication, acetosalicylic acid for the prevention of myocardial infarction, anticoagulants and antibiotics. The respondents were asked about the last time they had used different medicine groups, and the response options were: (1) within the last week; (2) 1–4 weeks ago; (3) 1–12 months ago; (4) more than a year ago and (5) never. For the analyses, medication use was dichotomised by combining the first three options (yes) and the last two options (no). Self-reports on doctor-diagnosed depression, hypertension and asthma were used as alternative outcomes in sensitivity analyses. A vast majority (95%) of the respondents had congruent responses between antihypertensive medication use and doctor-diagnosed hypertension and asthma medication use and doctor-diagnosed asthma. There were no comparable disease questions in the questionnaire about the use of psychotropic medication.

Information about green and blue views from home was obtained by asking the respondents whether they had a view of green or blue space from any of the windows of their dwelling. In the questionnaire, green areas were defined as forests, parks, fields, meadows, boglands and rocks (and any playgrounds and play fields in them), and blue areas as sea, lakes and rivers. The response options were: (1) no; (2) yes, I look at it seldom; (3) yes, I look at it occasionally and (4) yes, I look at it often. Options 2 and 3 were combined for the analyses.

To assess the frequency of green space visits, respondents were asked how often they spent time or exercised outdoors in green spaces (parks, forests, meadows, etc) during the warm season (May–September). The response options were: (1) never; (2) <1 time per week; (3) 1–2 times per week; (4) 3–4 times per week and (5) ≥5 times per week. The first two options were combined for the analyses.

In the Urban Atlas 2012, class 50 000 (water) includes the total area of lakes and rivers (minimum mapping unit of 1 ha). For sea areas, only coastal areas are included in class 50 000, whereas open sea areas are unclassified. Thus, open sea areas were taken from the National Land Survey of Finland Topographic Database 12/2015. If water data from the National Land Survey of Finland Topographic Database had been used instead of the Urban Atlas, 95 respondents out of 5987 would have been classified as blue space ‘>10%’ instead of ‘>0–10%’. This difference is negligible (1.6% of the respondents) ( online supplemental figure 2 ). More than half (54%) of the respondents had no blue space around their homes; thus, the amount of blue space within 1 km radius around the respondents’ home (%) was categorised (0%, >0%–10%, >10%). Cut-off points for the categories were selected to ensure a sufficient number of observations in each category.

The amounts of residential green and blue spaces were calculated as proportions of green space (%) and blue space (%) in a buffer around the centre point coordinates of the respondent’s residential building, based on high-resolution land use and land cover data from the Urban Atlas 2012 26 ( online supplemental figure S1 ). The amount of green space within a 1 km radius around the respondent’s home (%) was chosen as the main green space exposure variable because 1 km is a widely used approximation of a 10 min walking distance. 27 Green space consisted of the following Urban Atlas land use classes (class code): green urban areas, meaning public green areas for predominantly recreational use, such as gardens, zoos, parks, castle parks and cemeteries (14 100); forests (31 000); herbaceous vegetation associations, including natural grassland and moors (32 000); and wetlands (40 000). Arable land—for example, fields (21 000) and pastures (23 000), was excluded because it was not thought to be accessible. To help the interpretation of results of the logistic regression analyses, results on the amount of residential green space (used as continuous variable) were scaled to represent 10% increase in exposure. In the sensitivity analyses, the amount of green space within a 300 m radius around the respondent’s home was used as an alternative exposure variable instead of 1 km because 300 m is a commonly used threshold for accessibility. 27

In comparison with the general population of Helsinki in 2015 and Espoo and Vantaa in 2016, based on Statistics Finland’s free-of-charge statistical databases, the proportion of women was higher among the questionnaire respondents (57% among the respondents vs 52% among the general population). Most respondents belonged to the age groups 45–54, 55–64 and 65–74 years. The proportion of respondents in the age group 45–54 years was similar to (19%), and slightly higher in age groups 55–64 years (20% vs 16%) and 65–74 years (20% vs 13%) than, that in the general population of Helsinki, Espoo and Vantaa. Additionally, the prevalence of medication use among the respondents was approximately similar to the prevalence of the corresponding diseases in Finland, based on the Duodecim online Health Library.

The Finnish Institute for Health and Welfare conducted the Helsinki Capital Region Environmental Health Survey in Helsinki, Espoo and Vantaa between 2015 and 2016 using a mailed self-administered questionnaire. 25 These three cities comprise the capital region and largest urban area in Finland. The survey gathered information on how the urban population experiences exposure to various environmental factors in their living environment, related health risks and personal health status. The former Population Register Centre (currently Digital and Population Data Services Agency) provided three random samples from native Finnish-speaking residents aged 25 years or older in Helsinki (n=8000, data collected between May and August in 2015), Espoo (n=4000, June−August 2016) and Vantaa (n=4000, June−August 2016). The total eligible sample size was 16 000, and a response rate of 46% resulted in 7321 participants. To maximise the number of observations, each exposure–endpoint combination was allowed to have its own total number of observations (n=5881–6031).

The only indication of effect modification based on the p value and consistency of effect estimates across categories was observed between the frequency of green space visits and annual household income in the model for psychotropic medication use. The effects of visiting green spaces were stronger among those reporting lowest annual household income (<€30 000), whereas these associations were non-existent among those reporting highest annual household income (>€90 000) ( online supplemental table 8 ).

The frequency of green space visits was associated with decreased odds of medication use. In the main model, when the reference was <1 green space visit per week, the odds for psychotropic medication use were 0.67 (95% CI 0.55 to 0.82) for those visiting green spaces 3–4 times per week and 0.78 (0.63 to 0.96) for those visiting green spaces at least five times per week ( table 3 ). The respective odds for antihypertensive medication use were 0.64 (0.52 to 0.78) and 0.59 (0.48 to 0.74) ( table 4 ), and for asthma medication use, 0.74 (0.58 to 0.94) and 0.76 (0.59 to 0.99) ( table 5 ). Overall, the effect estimates remained essentially similar for the different models with different sets of confounders.

Amounts of residential green and blue spaces were not associated with self-reported use of psychotropic ( table 3 ), antihypertensive ( table 4 ), or asthma ( table 5 ) medications, although OR estimates for the amount of residential green space and asthma medication use (OR=1.06, 95% CI 0.99 to 1.14) and amount of residential blue space >0–10% and psychotropic medication (1.10, 0.95 to 1.28) were >1. Similarly, green or blue views from home were not associated with medication use ( tables 3–5 ), but the OR estimates were >1 for looking at the green view seldom or occasionally and both antihypertensive (1.07, 0.88 to 1.30) and asthma (1.11, 0.88 to 1.39) medications.

The characteristics of the study population are presented in tables 1 and 2 , respectively. The Pearson correlation coefficients between continuous environmental exposures and SES indicators, means of continuous exposures and SES indicators by categorical exposures and SES indicators, and frequencies of categorical exposures by SES indicators are presented in online supplemental tables S1−S7 . For example, those who reported visiting green spaces often had a slightly higher amount of residential green spaces around their home ( online supplemental table 2 ), but the differences between the categories were negligible.

Discussion

In this cross-sectional study in the Helsinki capital region of Finland, a higher frequency of green space visits was associated with a lower frequency of psychotropic, antihypertensive and asthma medication use, and the association was not dependent on SES. The amounts of residential green and blue spaces or green and blue views from home were not associated with medication use, although there were some increased estimates, and due to limited sample size, associations cannot be ruled out.

Previous evidence on the effects of nature exposure on mental, cardiovascular and respiratory health is somewhat mixed.1 2 5–12 In line with our results for the association between nature exposure and mental health, in a survey across 18 countries (n=16 703), no clear associations were observed between green or blue space exposure indicators and antidepressive or anxiety medication use,31 and in a cross-sectional study of ageing cohorts in four European studies (n=16 189), no association between distance to the nearest green space and self-reported depressed affect was observed.32 In contrast to our results, a recent longitudinal study from Finland (n=19 851) observed an association between high residential greenness and low risk of doctor-diagnosed depression33; however, no information on visits to green space were used. Contrary to our results, a Spanish cross-sectional study (n=958) found that surrounding greenness was associated with decreased odds of self-reported benzodiazepine use.34

For hypertension, a cross-sectional study using data from Barcelona and Brussels (n=5735) reported an association between distance to the nearest green space and self-reported hypertension in Barcelona but not in Brussels.35 Another cross-sectional study from Brazil (n=3418) found a beneficial association between increased number of street trees and decreased distance from urban parks and decreased prevalence of hypertension diagnosis.36

The associations between nature exposure and respiratory diseases, such as asthma, are thought to be bidirectional. Vegetation can improve air quality and is a source of beneficial microbial agents; however, it is also a source of allergens, which may exacerbate allergic responses and interact with air pollution.37 Additionally, beneficial effects have been observed mainly in children and adolescents.38 When studying asthma, the results can be affected by geographical study area, types of vegetation and other environmental exposures, as well as seasonality.11 A review of mainly cross-sectional studies reported a suggestive association between increased green space exposure and decreased incidence of asthma.8 In the present study, more frequent green space visits were associated with less frequent asthma medication use, but this association was attenuated after adjusting for BMI. This might be explained by the fact that obesity is a known risk factor for asthma.39 In contrast, in the present study the OR estimate for asthma medication use increased with increasing amounts of residential green space. Similar increased estimates have been previously observed—for example, in a Spanish cross-sectional study among children (n=3178), regarding proximity to green space and self-reported current asthma.14

In the present study, the frequency of green space visits was the only type of nature exposure that showed an inverse association with medication use. This finding is in line with tentative evidence emphasising the importance of actual use of green space in relation to mental health,19 31 and it suggests that the same holds true for other health conditions, such as asthma and hypertension. The absence of an association between objectively evaluated nature exposure and health in this and some earlier studies can be explained, for instance, by the lack of information on the quality of green areas,40 or the fact that a mere amount of accessible green space does not necessarily indicate an actual interaction between the individual and surrounding environment.10 41 It is also possible that in certain populations, an association between nature exposure and health cannot be observed. It should also be noted that Finland has a high forest coverage, high frequency of nature visits and the majority of outdoor recreation occurs in municipality-owned areas.42 43 In addition, Finnish cities are relatively green44 45; thus, it is possible that the amount of green space is not a limiting factor in this study, but that those willing to use green space can access it with minimal effort.

For blue space exposure, there is suggestive evidence of an association with better mental well-being.2 For example, a recent cross-sectional study conducted in the UK (n=21 097) observed an inverse association between residential freshwater coverage and mental health,28 whereas a Spanish cross-sectional study did not.34

Evidence regarding green or blue views from homes is scarce. Contrary to our results, a small Spanish cross-sectional study (n=465) suggests that enjoying a view of green space from home could decrease the risk of self-reported depression.21

In the present study, many of the observed associations between residential green spaces and medications were attenuated when BMI was included in the models, indicating that the potential association is likely mediated, at least partially, by BMI. This was also reflected in the observation that those who reported visiting green spaces often had a slightly lower BMI than those who visited green spaces less often.

No consistent effect modification by SES was observed in this study, and there was practically no variation in the frequency of green space visits or the amount of residential green space according to SES. However, in line with previous findings, it was observed that the inverse association between the frequency of green space visits and psychotropic medication use was most evident among those reporting the lowest annual household income. In addition to differences in, for example, study populations and designs, the lack of a consistent effect modification could be explained by the fact that in Europe, public green spaces, such as parks, have better quality in lower-SES neighbourhoods compared with, for example, those in North America.24 Similar inconsistent effect modification by SES was observed in a recent Finnish longitudinal study on residential greenness and doctor-diagnosed depression.33