Neighbourhood and path-based greenspace in three European countries: associations with objective physical activity
BMC Public Health volume 21, Article number: 282 (2021)
Greenspace has been associated with health benefits in many contexts. An important pathway may be through outdoor physical activity. We use a novel approach to examine the link between greenspace microenvironments and outdoor physical activity levels in the HEALS study conducted in Edinburgh (UK), the Netherlands, and Athens and Thessaloniki (Greece).
Using physical activity tracker recordings, 118 HEALS participants with young children were classified with regard to daily minutes of moderate to vigorous physical activity (MVPA); 60 were classified with regard to the metabolic equivalent task (MET)-minutes for each of the 1014 active trips they made. Greenspace indicators were generated for Normalised Difference Vegetation Index (NDVI), tree cover density (TCD), and green land use (GLU). We employed linear mixed-effects models to analyse (1) daily MVPA in relation to greenspace within 300 m and 1000 m of residential addresses and (2) trip MET-minutes in relation to average greenspace within a 50 m buffer of walking/cycling routes. Models were adjusted for activity, walkability, bluespace, age, sex, car ownership, dog ownership, season, weekday/weekend day, and local meteorology.
There was no clear association between MVPA-minutes and any residential greenspace measure. For example, in fully adjusted models, a 10 percentage point increase in NDVI within 300 m of home was associated with a daily increase of 1.14 (95% CI − 0.41 to 2.70) minutes of MVPA. However, we did find evidence to indicate greenspace markers were positively linked to intensity and duration of activity: in fully adjusted models, 10 percentage point increases in trip NDVI, TCD, and GLU were associated with increases of 10.4 (95% CI: 4.43 to 16.4), 10.6 (95% CI: 4.96 to 16.3), and 3.36 (95% CI: 0.00 to 6.72) MET-minutes, respectively. The magnitude of associations with greenspace tended to be greater for cycling.
More strenuous or longer walking and cycling trips occurred in environments with more greenspace, but levels of residential greenspace did not have a clear link with outdoor MVPA. To build on our research, we suggest future work examine larger, more diverse populations and investigate the influence of greenspace for trip purpose and route preference.
Increased residential greenspace (e.g., parks) or greenness (e.g., street trees) has shown to be associated with beneficial health, such as better self-reported health and reduced all-cause and cardiovascular mortality . Research has now progressed to explore potential causal mechanisms. As strong links have been made between physical activity (PA) and numerous health outcomes, particularly for cardiovascular outcomes , an important pathway to health may be access to areas in which to engage in PA. Moreover, though still an active research area, exercise specifically undertaken in green areas may enhance the proven benefits of PA .
Nevertheless, research on the importance of greenspace for exercise has produced mixed results. Cross-sectional studies relying on self-reported data to assess the relationship between residential greenspace and PA identified positive associations in populations in Australia , Canada , and the US , while other work in Denmark , Netherlands , and Scotland  found no such links. With the emergence of low-cost GPS-equipped sensors and devices , researchers can now better track objective measures of PA and actual greenspace use, though these studies too have found equivocal results: the amount of residential greenspace was related to higher levels of overall moderate to vigorous PA (MVPA) , but in another study, associations were found only with PA when undertaken within green areas (i.e., not overall PA) .
Recommendations from agencies, including the World Health Organization (WHO), prescribe a minimum weekly dose of 150 min of moderate intensity or 75 min of vigorous PA, yet a recent global survey found over a quarter of individuals were not achieving these salubrious levels . Though greenspace may help promote active travel and facilitate outdoor PA, for example, through appealing tree-lined streets or accessible parks, other neighbourhood attributes, such as overall walkability (e.g., street connectivity, population density, mixed use development) and access to services, have been found to be more important [14, 22]. Even if a positive link with greenspace is established, a further complicating factor is that self-selection may bias findings if healthier individuals choose to live in greener areas with more options for outdoor exercise ; if present, this bias would result in exaggerated health benefits of greenspace.
Our study explored two distinct research questions to advance our understanding of the association of greenspace and PA within the built environment: 1) whether the availability of residential greenspace is associated with increased MVPA and 2) whether individuals choose routes with on average higher greenspace levels for longer/more active journeys. In addition, for the second question, we also assessed the greenspace associations separately for walking and cycling trips.
Study design and population
Data were obtained from the EU-funded study, Health and Environment-wide Associations based on Large population Surveys (HEALS; http://www.heals-eu.eu), which employed indoor and personal sensors to characterise the environments of families with young children. The study included a sample of households concentrated in Edinburgh, UK; Utrecht and elsewhere in the Netherlands; and Thessaloniki and Athens, Greece. Individuals aged 18 years or older with a young child (< 3 years of age) were eligible to participate in the HEALS study (n = 131) and were recruited through advertising via universities, childcare groups, and word of mouth. Informed written consent was provided by all participants. Personal monitoring periods lasted approximately 1 week during 2015 and 2016 and entailed indoor monitoring of air pollutants and noise and the participant wearing a physical activity tracker device. Questionnaires were developed in the HEALS study to gather household data, including socioeconomic position (SEP) (see supplementary material).
We assigned three indicators of urban greenspace: the Normalised Difference Vegetation Index (NDVI), tree cover density (TCD), and green land use (GLU), similar to a previous analysis using the HEALS dataset published by the authors . Each indicator provides potentially overlapping, but distinct, perspectives of greenspace: NDVI (− 1 to + 1) represents the overall greenness of a given area, TCD provides the percentage (0–100%) of an area covered by the canopy of trees as visible from satellites, and GLU indicates areas used for specific types of green land (parks, forests, sports pitches, etc.) (see Fig. 1).
For each study area, NDVI values were calculated using Sentinel-2 satellite images available from the Copernicus Open Access Hub at 10-m spatial and five-day temporal resolutions. NDVI raster data with values of <− 0.1 represent water or ice and were excluded from greenness calculations . Images from summer with cloud coverage of < 10% were selected to maximise spatial contrasts of greenness. Images produced within 1 year of the personal monitoring periods were retrieved, except for those in and around Edinburgh, due to cloud coverage (See Table S1 for exact image dates). Average annual TCD based on Sentinel-2 and Landsat 8 satellite images (20 m spatial resolution) for Europe in 2015 was also obtained from the Copernicus Hub. Coastal waters were excluded in the calculation of TCD values. GLU was based on CORINE land use data (2012), which has been refined subsequently through data fusion with other spatial datasets (e.g., Urban Atlas, OpenStreet Map) and is publicly available as a 100 m raster dataset . Unlike the original CORINE dataset, this enhanced version distinguishes between green and non-green sport and leisure facilities. The following categories were combined to create a GLU map: green urban areas, green sport and leisure facilities, broad-leaved forest, coniferous forest, mixed forest, natural grasslands, moors and heathland, sclerophyllous vegetation, and transitional woodland-shrub. Mean values of NDVI and TCD, and the proportion of GLU, were calculated in 300 m and 1000 m radial buffers around home addresses. These sizes were selected to represent a reasonable walking distance to greenspace (300 m; ) and to reflect a larger, neighbourhood scale (1000 m; ). Additional details of the methods for each indicator can be found in Mueller et al. .
During the personal monitoring periods, study participants wore a Fitbit flex device (original version) on their wrist (Fitbit Inc., San Francisco, CA, USA)  and installed the ‘Moves’ app (moves-app.com)  on their mobile phones; participants were asked to keep their Fitbit and phone with them whenever possible. Fitbits recorded the total number of steps completed each minute and the Moves app recorded GPS locations and the duration, distance, and activity (i.e., walking, running, cycling, vehicle transport) based on its algorithm to identify discrete trips. The Fitbit flex has been found to reliably record steps compared to gold standards (Optogait system and ActivePAL device) , and the Moves app can correctly record the location and type of separate trips [4, 47].
To take advantage of both the physical activity sensor and mobile phone app deployed in the HEALS study, we derived two PA metrics that made use of the particular data provided by each sensor: daily minutes of MVPA steps (Fitbit) and Metabolic Equivalent Task minutes (MET-minutes) (Moves app); METs represent the energy cost of an activity relative to a resting state . Daily steps were calculated by summing minutes with ≥100 steps as recorded by the Fitbit flex (equivalent to ≥3 METs)  across the monitoring period. These daily values were then divided by the number of days with at least 12 h of data (i.e., 75% complete data, assuming 8 h of sleep), where at least four such days had been recorded during the monitoring period. Out of 133 individuals who were provided Fitbits (some households had multiple participants), 124 (93%) provided sufficient data for analysis.
MET-minutes were calculated by assigning a specific MET to those trips identified by the Moves app to be ‘walking’, ‘running’, or ‘cycling’, depending on the activity; average speed (based on distance and duration, as recorded by Moves); and overall grade change (steepness) during each trip using values set out in Ainsworth et al. . To account for steepness in the calculation of METs, topographical GIS maps (30 m resolution) were acquired from the Japan Aerospace Exploration Agency, based on the Advanced Land Observing Satellite (ALOS-2; ). Where no METs were specified by Ainsworth et al.  for a given combination of activity/speed/grade, values were interpolated or extrapolated (n = 3) (see Table S2 for a complete list of METs used in analysis). METs were multiplied by the duration of each trip to calculate MET-minutes. GPS points were converted to lines in QGIS v.3.10.1  and visual inspection was used to remove trips either with straight lines that did not appear to follow road networks or that traversed bodies of water (n = 16). Only six trips were assigned as ‘running,’ which were subsequently excluded from analysis. Values above five standard deviations (SD) in excess of the mean were excluded for MET-minutes (n = 7) and duration (n = 2). To select trips that occurred outdoors, those of < 3 min in duration or < 100 m in distance were excluded from analysis. As with the daily steps calculation, Moves data were used only from individuals with at least 18 h (i.e., 75%) of complete data on four or more days during the monitoring period. Out of 123 individuals who downloaded Moves onto their phones, 69 (56%) provided sufficient data for analysis. Since few (n = 4) participants in Thessaloniki generated sufficient Moves data, this study centre was excluded from the trip-based analysis.
As certain features of the built environment may be more likely to encourage physical activity , we calculated walk scores to capture the degree of walkability of residential and travel environments. Similar to previous studies (e.g., [19, 22, 57]), walk scores were calculated based on GIS data using three factors: population density, intersection counts, and land use mix. As well as walking, these same built environment factors may also encourage cycling . Population density was based on global 1 × 1 km gridded estimates for 2015 . Intersection counts were calculated using QGIS via road networks from OpenStreetMap shapefiles downloaded during March–April 2019 from Geofabrik (https://download.geofabrik.de/). Auto-oriented (i.e., non-pedestrian accessible) roads were removed by deleting feature classes for ‘motorway’, ‘service’, or ‘trunk’, and the processing tool in QGIS, ‘v.clean’, was employed to identify intersections of two or more distinct roads. Land use mix was based on the refined CORINE dataset, including ‘commercial/service facilities’, ‘public facilities’, and ‘sport and leisure green/built-up’. Z-scores of each walk score (i.e., mean population density, total intersection counts, and presence of specific land uses) were calculated across all home addresses for the 300 m and 1000 m buffers and were summed to create a walk score. Walk scores calculated separately within and across study areas were highly correlated for both 300 m (r = 0.88) and 1000 m (r = 0.91) buffers; the latter metric was used for analysis.
To examine the association of greenspace and MET-minutes between different trips taken by the same individual, linear buffers of 50 m were generated for each trip for which mean values of NDVI and TCD, as well as the proportion of GLU, were calculated; a smaller buffer size has been shown to be most strongly associated with MVPA . To account for different trip distances, the number of intersections within each trip buffer was divided by the total distance, which was then used to calculate walk scores in a similar fashion as described above.
As well as walkability, we adjusted for bluespace, daily meteorology, and season as other environmental factors. We accounted for bluespace by identifying any bodies of water in residential and trip buffers, as bluespace has been shown to be positively correlated with physical activity, especially walking [16, 38]. We included in our definition of bluespace the following CORINE land cover types: ‘water courses’, ‘water bodies’, ‘coastal lagoons’, ‘estuaries’, and ‘sea and ocean’. We obtained for the dates of the personal monitoring periods weather data, including daily maximum temperature and wind speed, and total precipitation from the US National Centers for Environmental Information  from the following stations (latitude, longitude): Edinburgh Royal Botanic Garden (55.967, − 3.210); Schiphol, Netherlands (52.316, 4.790); and Hellinikon, Greece (37.900, 23.750). Season was assigned to each monitoring period based on the majority of dates that occurred in a given season. As noted above, during the monitoring periods, participants also completed questionnaires on SEP and other information, including employment status (e.g., working, in school, caring for family), highest education completed, car ownership, and household pets.
We used mixed regression methods to examine associations between greenspace and physical activity metrics. Each greenspace metric (mean NDVI score, mean TCD, and proportion of GLU,) was rescaled such that regression coefficients represented the change in outcome for a 10 percentage point increase in the relevant parameter, an approach adopted by Mueller et al. .
Models were developed to assess:
the between-individuals association of MVPA with residential greenspace (seeking to answer the question of whether people living in greener areas have higher levels of MVPA),
the association, within individuals, of MET-minutes with trip-based greenspace (seeking to answer the question of whether longer/more active journeys are undertaken in areas with more greenspace compared with shorter/less active journeys).
For (i), with daily MVPA-minutes as the outcome, regression models with a random intercept for study centre were separately developed for residential greenspace metric at 300 m and 1000 m buffers around the home. Model results are presented with various levels of pre-specified confounder adjustment: (1) an unadjusted model, (2) a model adjusted for age using cubic splines with three knots, sex, season, and bluespace (any), and (3) a model with additional adjustment for car ownership, dog ownership, walk score, education, and employment.
For (ii), regression models with random intercepts for both study centre and individual and robust standard errors were separately developed for each of the three greenspace metrics: NDVI, TCD, and GLU. Results are again presented with adjustment for different sets of pre-specified confounders: (1) an unadjusted model, (2) a model with adjustment for age, sex, season, and bluespace (any), and (3) a model with additional adjustments for education, employment status, walk score, day of week, weather conditions on the day of activity, mean residential greenspace (1000 m buffer), car ownership, and dog ownership. Effect modification by activity (i.e., walking and cycling) was examined by including in the models an interaction term between greenspace metric and activity. Cubic splines were included into the model for age and temperature. Geospatial analysis was performed using QGIS and statistical analysis was undertaken using Stata v15 .
A total of 131 households enrolled in the HEALS study across the four study centres, with personal monitoring periods spanning from March 2015 to June 2016. There were 118 and 60 individuals who provided sufficient data and for whom covariate data were available in the neighbourhood and trip-based greenspace analyses, respectively. Descriptive characteristics pertaining to those individuals are presented in Table 1. The mean duration of MVPA-minutes was just under 12 min per day, with a maximum of nearly 40 min. The number of trips recorded for each participant ranged from one to 96, with a mean of 30.3 (SD = 23.8); the mean trip duration was just over 9 min. There was a total of 1014 trips, of which 676 (66.7%) were walking and 338 (33.3%) cycling; 89.9% (n = 304) of the cycling trips were in the Netherlands. The mean METs for each trip was 3.8; when accounting for duration, mean MET-minutes equated to 37.0.
Mean residential greenspace values were slightly higher for the 1000 m compared to the 300 m buffer (Table 1). The average trip-based NDVI was 0.27, with minimum and maximum values of − 0.04 and 0.83, respectively. Trip-based TCD levels ranged from 0 to 73.5%, with 85.5% (n = 864) of trips containing tree cover. The percentage of trips with any GLU was 31.2% (n = 316), with three (0.3%) trips occurring entirely in places of GLU. The greenspace metrics were weakly to moderately correlated, with NDVI and TCD consistently having the strongest associations. Greenspace metrics were mostly negatively correlated with walk score. There was little apparent correlation between residential greenspace metrics and daily MVPA-minutes. By contrast, trip-based greenspace was moderately correlated with MET-minutes, with coefficient values ranging from 0.44 (GLU) to 0.59 (TCD) (Table 2).
The analysis of residential greenspace and MVPA-minutes did not provide clear evidence of associations with greenspace at either the 300 m or 1000 m buffers (Table 3). Coefficients of the increase in MVPA were generally small, and confidence intervals included 0 in fully adjusted models for all greenspace metrics (Table 3). Of the covariates, only walk score in the NDVI model (300 m buffer) showed a clear positive trend (1.13 MVPA-minutes [95% CI: 0.03 to 2.23]) per 1-unit increase in walk scores in fully adjusted models (data not shown).
All average trip-based greenspace coefficients were positively associated with MET-minutes in the unadjusted and adjusted models. NDVI and TCD were most strongly related to MET-minutes, compared to GLU, with very similar coefficient values (10.41 [95% CI: 4.43 to 16.39] and 10.63 [95% CI: 4.96 to 16.30] additional MET-minutes per 10 percentage point increase, respectively). Although less precise, estimates of the absolute increase in MET-minutes for cycling trips were consistently higher than those for walking (Table 4). Select environmental covariates also were positively linked with MET-minutes across the greenspace models, particularly walk score and the presence of bluespace (data not shown).
Proximity to greenspace, typically in a residential setting, has been associated with a host of positive health outcomes. In this study, we used objective indicators to explore greenspace and outdoor PA as a potential underlying mechanism for health. We found no evidence to suggest individuals who lived in greener neighbourhoods engaged in greater levels of MVPA than those residing in less green areas. On the other hand, we found strong support that individuals choose greener settings for physically active travel of higher intensity and/or longer duration.
We found no clear evidence that the amount of greenspace around the home was associated with overall MVPA. A similar finding has been reported in some studies [53, 54] but not in others [23, 43], with some of the earlier work examining comparable residential greenspace metrics and objective PA, the majority of which examined GLU as the exposure of interest. The number of parks within a 1 km residential buffer, but not the residential distance to the nearest park, was associated with objective MVPA in a group of US adults . Likewise, the number of parks within 500 m and 1 km buffers was also found to be the strongest indicator for MVPA minutes in an eight-country study; park area within those same buffer sizes (a metric similar to the GLU metric in the current study) did not indicate a correlation with PA . Sallis et al.  also found parks within 500 m of residential addresses to be positively associated with objective MVPA, after adjusting for walkability features (also significant), in a large sample of individuals from 10 countries. A study examining GLU (i.e., parks and other green land uses) and objective MVPA in Dutch adults aged 45–65 years found positive results, but only with smaller buffers (25–400 m) . Triguero-Mas et al.  found overall MVPA activity was not associated with GLU situated within 300 m of home addresses in European adults, but was associated with contact and exercise specifically in natural outdoor environments; researchers did not account for walkability. We identified only one previous study that examined residential NDVI, which found no statistical links with overall objectively measured MVPA, and an inverse relationship with MVPA within a 1 km home buffer, in a sample of adult trail users in the US, . We are unaware of any previous studies that compare the amount of residential tree canopy to objective measures of PA.
While some studies have found positive correlations between residential greenspace and objective MVPA, albeit mainly with the number of nearby parks, the existing evidence is neither consistent nor comprehensive. Our study found a mix of positive and negative greenspace effects, which may have achieved statistical significance (in either direction) with a larger sample size. Sample size notwithstanding, there are several reasons that may explain the lack of stronger findings: walkability indicators have typically been shown to be as or more important than nearby greenspace  (identified in the current study), the physical environment may be less important to influence exercise in parents of young children , PA in nearby parks has been found to constitute a small proportion of overall PA , and perhaps most pertinent is that MVPA may have occurred outside the 300 m and 1000 m buffers employed in the present study. Most participants in our study owned a car; Hillsdon et al.  found that car owners engaged in more than 60% of outdoor PA outside of the neighbourhood, as defined by an 800 m residential buffer. Thus, the amount of greenspace within a residential area may not be as important for people with access to a vehicle.
In our analysis of trip-specific data, we found positive links between the amount of vegetation (NDVI) and tree coverage, and to a lesser degree GLU, with longer and more active journeys. Few studies have used a GPS approach to combine greenspace exposure with objective PA in adults, but all have found some indication of a positive trend with PA. James et al.  assessed momentary exposure to NDVI, as opposed to trip-level averages as analysed in the current study, in female nurses in the US and found a positive relationship with accelerometer counts per minute, particularly when walkability was low. A study of a similar design to that of James et al. recruited trail users in the US and found NDVI to be positively associated with a higher likelihood of MVPA . Houston  used a land cover map (including greenspace as tree canopy, irrigated grass cover, or non-irrigated grass cover/bare soil) and identified significant positive associations with the likelihood of adults engaging in MVPA. The amount of GLU at trip origin and end was associated with a higher probability of walking in a study in France, which found that trip-level characteristics outweighed those of the residential environment . A study of adults in Barcelona that also used the Moves app found both the proportion of large parks and tree density along routes to be positively associated with walking minutes .
We found higher effects of greenspace on cycling compared to walking, though the former had a wider range of possible effects. Few previous studies have examined greenspace with objective adult physical activity measures of both walking and cycling. Le et al.  quantified the built environment surrounding bicycle and pedestrian counters in 20 US cities and found a greater positive effect on cycling than walking (though greenspace and bluespace were combined in their analysis). Our results with objective measures support studies of self-reported cycling. Commuters in Barcelona were more likely to be cyclists with higher greenness in the study/work environment; interestingly, the greenness of the route was not significant, though commuting journeys were estimated by shortest distance rather than those actually travelled . Questionnaire respondents in Stockholm reported greenery to be one of the most important factors to stimulate cycle commuting . Although we looked at all active trips (i.e., not just those for commuting), our results build on this earlier research to suggest that greenness, through both overall vegetation and trees, might enhance and encourage all active transport by providing a more pleasant route.
We examined both residential and active transport environments, which provided an opportunity to compare and contrast these exposures using the same dataset. We found no evidence to support the residential environment being associated with objective MVPA, though our analysis was based on steps and therefore would have only pertained to walking or running. This analysis also only related to the availability of greenspace, not necessarily its use. We also examined greenspace levels of the entire route for those trips involving walking or cycling. Whereas contemporaneous momentary designs (i.e., matching exposure and PA at points in time) are more likely to reveal typical behaviours in certain settings (e.g., less PA in commercial areas and more PA in natural areas, such as greenspaces) , our analysis took into account average characteristics of the entire route. Therefore, our approach was more equipped to answer the question: given an individual has decided to undertake PA, how is greenspace associated with the intensity and duration of activity? In other words, how does the presence of greenspace factor in the selection of environments through which individuals choose to travel or exercise? We found clear evidence indicating both NDVI and TCD as greenspace markers were positively linked to intensity and duration of activity, while adjusting for other characteristics of the built environment. Certain such characteristics, namely walk score, were consistently related to higher levels of PA; nevertheless, the different scales of greenspace markers and walk score render it difficult to identify which is the more influential factor for PA. We also found positive links to MVPA with the proportion of GLU along a route, but not specifically for cycling trips. The use of a particular greenspace for a specific activity, namely cycling in this case, may be more dependent on certain features, including size, cycling routes, and wooded areas, which were not quantified explicitly in the overall area-based GLU metric employed in our study . In addition, the GLU map we used was based on a lower spatial resolution (100 m) than the NDVI (10 m) or TCD (20 m) metrics. Therefore, the use of this coarser resolution, with greater aggregation of features and potential exclusion of smaller parks, might help explain the weaker associations we observed between GLU and PA indicators .
Strengths and limitations
Our study had several key strengths. We assessed the importance of both the residential and active route settings, thus developing dynamic and multi-contextual environmental exposures , with two objective MVPA indicators. We also used three different objective indicators to help characterise greenspace features of the built environment, with two different residential buffer sizes to help address the modifiable areal unit problem . These advantages notwithstanding, there were some limitations to our research. Although we did not explicitly address reasons for choosing residential locations, we attempted to control for self-selection in the trip-based analysis by including residential greenspace levels and found our results to be unchanged. Several greenspace and PA studies have attempted to account for self-selection by including reasons for choosing to live in their neighbourhood (e.g., access to places that support PA, access to local services). Associations with PA have persisted after adjustment for such factors [24, 34]. Therefore, it is not likely that residential self-selection would have strongly biased our results. However, it would have been beneficial for our analysis, and understanding of the importance and role of greenspace, to know the purpose(s) of each trip.
While the Moves app has been shown to accurately provide location, speed, and duration, the software has had challenges identifying multi-modal trips, which may have been included as discrete events in our analysis . In addition, there was a lower proportion of participants with complete Moves data than that provided by the Fitbit; this might be due to phones running out of batteries or being switched off. Our sample size was quite modest, and our study demographic was limited to parents of young children, which could restrict the generalisability of key findings. Although Candelaria et al.  found little difference in the amount of objective MVPA recorded between parents of young and older children and non-parents, the mean MVPA-minutes in our sample was much lower than Candelaria et al. and in studies with other demographics  (~ 12 vs > 30 mins/day). If MVPA steps were underestimated in our study, any association with residential greenspace levels might have been hindered. The majority of our study sample was university educated and owned a car, indicative of a higher SEP; lower SEP individuals might experience different relationships between greenspace and MVPA . As noted above, the environments of study/work may be important, but we did not have this information for all study participants. We also were not able to distinguish whether study subjects were currently working or on maternity/paternity leave. We characterised surrounding streets and intersections using maps from 2019, though personal monitoring took place over 2015–2016; therefore, some misclassification of walkability may have been introduced by any road network changes occurring in the intervening years, but it is expected that any impact on our results would have been minimal. Each subject participated in only one personal monitoring period in the HEALS study; repeating data collection with participants during different times of the year may provide insights into the role of temporal/seasonal factors of greenspace and PA.
We examined PA as a potential explanatory pathway for observed associations between health and greenspace, assessing both residential and trip-specific environments. We found little evidence to suggest residential greenspace was associated with higher levels of MVPA, regardless of where that may take place. On the other hand, we found clear, positive associations between intensity and duration of activities with the average amount of greenness and tree coverage along a route, which was true for both walking and even more so for cycling. We suggest future research to build on this proposed model of specific pathways by examining larger, more diverse populations, while also investigating the influence of greenspace for trip purpose and route preference.
Availability of data and materials
Access to the personal data used in this study is governed by the HEALS Data Management Plan. The environmental datasets obtained in this study are publicly available from the specific references herein. The datasets generated and/or analysed during the current study are not publicly available due to reasons of ensuring anonymity for study subjects. Specific data may be made available from the corresponding author on reasonable request.
Green land use
Health and Environment-wide Associations based on Large population Surveys
Metabolic task equivalent
Moderate to vigorous physical activity
Normalised Difference Vegetation Index
Tree cover density
World Health Organization
Ainsworth BE, Haskell WL, Herrmann SD, Meckes N, Bassett DR Jr, Tudor-Locke C, Greer JL, Vezina J, Whitt-Glover MC, Leon AS. 2011 compendium of physical activities: a second update of codes and MET values. Med Sci Sports Exerc. 2011;43(8):1575–81.
Astell-Burt T, Feng X, Kolt GS. Green space is associated with walking and moderate-to-vigorous physical activity (MVPA) in middle-to-older-aged adults: findings from 203 883 Australians in the 45 and up study. Br J Sports Med. 2014;48(5):404–6.
Browning M, Lee K. Within what distance does “greenness” best predict physical health? A systematic review of articles with GIS buffer analyses across the lifespan. Int J Environ Res Public Health. 2017;14(7):675.
Bucher D, Cellina F, Mangili F, Raubal M, Rudel R, Rizzoli AE, Elabed O. Exploiting fitness apps for sustainable mobility-challenges deploying the GoEco! App. In: ICT for sustainability 2016. Amsterdam: Atlantis Press; 2016.
Candelaria JI, Sallis JF, Conway TL, Saelens BE, Frank LD, Slymen DJ. Differences in physical activity among adults in households with and without children. J Phys Act Health. 2012;9(7):985–95.
Carson V, Rosu A, Janssen I. A cross-sectional study of the environment, physical activity, and screen time among young children and their parents. BMC Public Health. 2014;14(1):61.
Center for International Earth Science Information Network - CIESIN - Columbia University. Gridded Population of the World, version 4 (GPWv4): population density, revision 11. Palisades: NASA Socioeconomic Data and Applications Center (SEDAC); 2018. https://doi.org/10.7927/H49C6VHW.
Chaix B, Kestens Y, Duncan DT, Brondeel R, Méline J, El Aarbaoui T, Pannier B, Merlo J. A GPS-based methodology to analyze environment-health associations at the trip level: case-crossover analyses of built environments and walking. Am J Epidemiol. 2016;184(8):579–89.
Chaix B, Meline J, Duncan S, Merrien C, Karusisi N, Perchoux C, Lewin A, Labadi K, Kestens Y. GPS tracking in neighborhood and health studies: a step forward for environmental exposure assessment, a step backward for causal inference? Health Place. 2013;21:46–51.
Cohen-Cline H, Turkheimer E, Duncan GE. Access to green space, physical activity and mental health: a twin study. J Epidemiol Community Health. 2015;69(6):523–9.
Cole-Hunter T, Donaire-Gonzalez D, Curto A, Ambros A, Valentin A, Garcia-Aymerich J, Martínez D, Braun LM, Mendez M, Jerrett M, Rodriguez D. Objective correlates and determinants of bicycle commuting propensity in an urban environment. Transp Res Part D Transp Environ. 2015;40:132–43.
Diaz KM, Krupka DJ, Chang MJ, Peacock J, Ma Y, Goldsmith J, Schwartz JE, Davidson KW. Fitbit®: an accurate and reliable device for wireless physical activity tracking. Int J Cardiol. 2015;185:138–40.
Evenson KR, Furberg RD. Moves app: a digital diary to track physical activity and location. Br J Sports Med. 2017;51(15):1169–70.
Gascon M, Götschi T, de Nazelle A, Gracia E, Ambròs A, Márquez S, Marquet O, Avila-Palencia I, Brand C, Iacorossi F, Raser E. Correlates of walking for travel in seven European cities: the PASTA project. Environ Health Perspect. 2019;127(9):097003.
Gascon M, Sánchez-Benavides G, Dadvand P, Martínez D, Gramunt N, Gotsens X, Cirach M, Vert C, Molinuevo JL, Crous-Bou M, Nieuwenhuijsen M. Long-term exposure to residential green and blue spaces and anxiety and depression in adults: a cross-sectional study. Environ Res. 2018;162:231–9.
Gascon M, Zijlema W, Vert C, White MP, Nieuwenhuijsen MJ. Outdoor blue spaces, human health and well-being: a systematic review of quantitative studies. Int J Hyg Environ Health. 2017;220(8):1207–21.
Giles-Corti B, Donovan RJ. Socioeconomic status differences in recreational physical activity levels and real and perceived access to a supportive physical environment. Prev Med. 2002;35(6):601–11.
Guthold R, Stevens GA, Riley LM, Bull FC. Worldwide trends in insufficient physical activity from 2001 to 2016: a pooled analysis of 358 population-based surveys with 1· 9 million participants. Lancet Glob Health. 2018;6(10):e1077–86.
Hajna S, Ross NA, Joseph L, Harper S, Dasgupta K. Neighbourhood walkability and daily steps in adults with type 2 diabetes. PLoS One. 2016;11(3):e0151544.
Hillsdon M, Coombes E, Griew P, Jones A. An assessment of the relevance of the home neighbourhood for understanding environmental influences on physical activity: how far from home do people roam? Int J Behav Nutr Phys Act. 2015;12(1):100.
Houston D. Implications of the modifiable areal unit problem for assessing built environment correlates of moderate and vigorous physical activity. Appl Geogr. 2014;50:40–7.
James P, Hart JE, Hipp JA, Mitchell JA, Kerr J, Hurvitz PM, Glanz K, Laden F. GPS-based exposure to greenness and walkability and accelerometry-based physical activity. Cancer Epidemiol Biomark Prev. 2017;26(4):525–32.
Jansen M, Kamphuis CB, Pierik FH, Ettema DF, Dijst MJ. Neighborhood-based PA and its environmental correlates: a GIS-and GPS based cross-sectional study in the Netherlands. BMC Public Health. 2018;18(1):233.
Kaczynski AT, Mowen AJ. Does self-selection influence the relationship between park availability and physical activity? Prev Med. 2011;52(1):23–5.
Kantomaa MT, Tikanmäki M, Kankaanpää A, Vääräsmäki M, Sipola-Leppänen M, Ekelund U, Hakonen H, Järvelin MR, Kajantie E, Tammelin TH. Accelerometer-measured physical activity and sedentary time differ according to education level in young adults. PLoS One. 2016;11(7):e0158902.
Kärmeniemi M, Lankila T, Ikäheimo T, Puhakka S, Niemelä M, Jämsä T, Koivumaa-Honkanen H, Korpelainen R. Residential relocation trajectories and neighborhood density, mixed land use and access networks as predictors of walking and bicycling in the northern Finland birth cohort 1966. Int J Behav Nutr Phys Act. 2019;16(1):88.
Klein Tank AMG, Wijngaard JB, Können GP, Böhm R, Demarée G, Gocheva A, Mileta M, et al. Daily dataset of 20th-century surface air temperature and precipitation series for the European Climate Assessment. Int J Climatol. 2002;22(12):1441–53.
Kooiman TJ, Dontje ML, Sprenger SR, Krijnen WP, van der Schans CP, de Groot M. Reliability and validity of ten consumer activity trackers. BMC Sports Sci Med Rehabil. 2015;7(1):24.
Kwan MP. The limits of the neighborhood effect: contextual uncertainties in geographic, environmental health, and social science research. Ann Am Assoc Geogr. 2018;108(6):1482–90.
Labib SM, Lindley S, Huck JJ. Spatial dimensions of the influence of urban green-blue spaces on human health: a systematic review. Environ Res. 2019;180:108869.
Le HT, Buehler R, Hankey S. Correlates of the built environment and active travel: evidence from 20 US metropolitan areas. Environ Health Perspect. 2018;126(07):077011.
Loh M, Sarigiannis D, Gotti A, Karakitsios S, Pronk A, Kuijpers E, Annesi-Maesano I, Baiz N, Madureira J, Oliveira Fernandes E, Jerrett M. How sensors might help define the external exposome. Int J Environ Res Public Health. 2017;14(4):434.
Maas J, Verheij RA, Spreeuwenberg P, Groenewegen PP. Physical activity as a possible mechanism behind the relationship between green space and health: a multilevel analysis. BMC Public Health. 2008;8(1):206.
McCormack GR. Neighbourhood built environment characteristics associated with different types of physical activity in Canadian adults. Health Promot Chronic Dis Prev Can. 2017;37(6):175–85.
McMorris O, Villeneuve PJ, Su J, Jerrett M. Urban greenness and physical activity in a national survey of Canadians. Environ Res. 2015;137:94–100.
Mueller W, Steinle S, Pärkkä J, Parmes E, Liedes H, Kuijpers E, Pronk A, Sarigiannis D, Karakitsios S, Chapizanis D, Maggos T. Urban greenspace and the indoor environment: pathways to health via indoor particulate matter, noise, and road noise annoyance. Environ Res. 2020;180:108850.
Ord K, Mitchell R, Pearce J. Is level of neighbourhood green space associated with physical activity in green space? Int J Behav Nutr Phys Act. 2013;10:127.
Pasanen TP, White MP, Wheeler BW, Garrett JK, Elliott LR. Neighbourhood blue space, health and wellbeing: the mediating role of different types of physical activity. Environ Int. 2019;131:105016.
QGIS Development Team. QGIS geographic information system: Open Source Geospatial Foundation Project; 2014. https://qgis.org
Rosina K, Batista e Silva F, Vizcaino P, Herrera MM, Freire S, Schiavina M. Increasing the detail of European land use/cover data by combining heterogeneous data sets. Int J Digit Earth. 2020;13(5):602–26. https://doi.org/10.1080/17538947.2018.1550119.
Rowe DA, Welk GJ, Heil DP, Mahar MT, Kemble CD, Calabró MA, Camenisch K. Stride rate recommendations for moderate-intensity walking. Med Sci Sports Exerc. 2011;43(2):312–8.
Saelens BE, Sallis JF, Frank LD, Cain KL, Conway TL, Chapman JE, Slymen DJ, Kerr J. Neighborhood environment and psychosocial correlates of adults’ physical activity. Med Sci Sports Exerc. 2012;44(4):637–46.
Sallis JF, Cerin E, Conway TL, Adams MA, Frank LD, Pratt M, Salvo D, Schipperijn J, Smith G, Cain KL, Davey R. Physical activity in relation to urban environments in 14 cities worldwide: a cross-sectional study. Lancet. 2016;387(10034):2207–17.
Schipperijn J, Bentsen P, Troelsen J, Toftager M, Stigsdotter UK. Associations between physical activity and characteristics of urban green space. Urban For Urban Green. 2013;12(1):109–16.
Schipperijn J, Cerin E, Adams MA, Reis R, Smith G, Cain K, Christiansen LB, Van Dyck D, Gidlow C, Frank LD, Mitáš J. Access to parks and physical activity: an eight country comparison. Urban For Urban Green. 2017;27:253–63.
Shanahan DF, Franco L, Lin BB, Gaston KJ, Fuller RA. The benefits of natural environments for physical activity. Sports Med. 2016;46(7):989–95.
Stamatelopoulou A, Chapizanis D, Karakitsios S, Kontoroupis P, Asimakopoulos DN, Maggos T, Sarigiannis D. Assessing and enhancing the utility of low-cost activity and location sensors for exposure studies. Environ Monit Assess. 2018;190(3):155.
StataCorp. Stata statistical software: release 15. College Station: StataCorp LLC; 2017.
Stewart OT, Moudon AV, Littman AJ, Seto E, Saelens BE. Why neighborhood park proximity is not associated with total physical activity. Health Place. 2018;52:163–9.
Tadono T, Nagai H, Ishida H, Oda F, Naito S, Minakawa K, Iwamoto H. Generation of the 30 M-mesh global digital surface model by Alos prism: ISPRS-International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences; 2016. p. 157–62.
Tamura K, Wilson JS, Goldfeld K, Puett RC, Klenosky DB, Harper WA, Troped PJ. Accelerometer and GPS data to analyze built environments and physical activity. Res Q Exerc Sport. 2019;90(3):395–402.
Thornton CM, Kerr J, Conway TL, Saelens BE, Sallis JF, Ahn DK, Frank LD, Cain KL, King AC. Physical activity in older adults: an ecological approach. Ann Behav Med. 2017;51(2):159–69.
Triguero-Mas M, Donaire-Gonzalez D, Seto E, Valentín A, Smith G, Martínez D, Carrasco-Turigas G, Masterson D, Van den Berg M, Ambròs A, Martínez-Íñiguez T. Living close to natural outdoor environments in four European cities: adults’ contact with the environments and physical activity. Int J Environ Res Public Health. 2017;14(10):1162.
Troped PJ, Wilson JS, Matthews CE, Cromley EK, Melly SJ. The built environment and location-based physical activity. Am J Prev Med. 2010;38(4):429–38.
Twohig-Bennett C, Jones A. The health benefits of the great outdoors: a systematic review and meta-analysis of greenspace exposure and health outcomes. Environ Res. 2018;166:628–37.
van den Bosch MA, Mudu P, Uscila V, Barrdahl M, Kulinkina A, Staatsen B, Swart W, Kruize H, Zurlyte I, Egorov AI. Development of an urban green space indicator and the public health rationale. Scand J Public Health. 2016;44(2):159–67.
van Dyck D, Cardon G, Deforche B, Sallis JF, Owen N, De Bourdeaudhuij I. Neighborhood SES and walkability are related to physical activity behavior in Belgian adults. Prev Med. 2010;50:S74–9.
Vich G, Marquet O, Miralles-Guasch C. Green streetscape and walking: exploring active mobility patterns in dense and compact cities. J Transp Health. 2019;12:50–9.
Wahid A, Manek N, Nichols M, Kelly P, Foster C, Webster P, Kaur A, Friedemann Smith C, Wilkins E, Rayner M, Roberts N. Quantifying the association between physical activity and cardiovascular disease and diabetes: a systematic review and meta-analysis. J Am Heart Assoc. 2016;5(9):e002495.
Wahlgren L, Schantz P. Exploring bikeability in a metropolitan setting: stimulating and hindering factors in commuting route environments. BMC Public Health. 2012;12(1):168.
We thank all of the participants in the HEALS study. Dr. Paul Kelly is thanked for his valuable input in the analysis.
The HEALS project was funded from the European Union’s Seventh Programme for research, technological development and demonstration under grant agreement No 603946. The funder had no role in the design of the study; collection, analysis, and interpretation of data; or in writing the manuscript.
Ethics approval and consent to participate
Ethical approval was received for each study centre (UK: Heriot Watt University Ethics Review Board 2015–07, London School of Hygiene and Tropical Medicine Observational / Interventions Research Ethics Committee 17028; Netherlands: Medisch Ethische Toetsingscommissie Brabant NW2015–07; Athens: National Centre for Scientific Research ‘Demokritos’ Ethics Review Board 2015–04: 260/2015–1671; Thessaloniki: Aristotle University Ethics Committee 140540/2018). Informed written consent was provided by all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Mueller, W., Wilkinson, P., Milner, J. et al. Neighbourhood and path-based greenspace in three European countries: associations with objective physical activity. BMC Public Health 21, 282 (2021). https://doi.org/10.1186/s12889-021-10259-0
- Physical activity