Next Article in Journal
Predictive Utility of Changes in Optic Nerve Sheath Diameter after Cardiac Arrest for Neurologic Outcomes
Next Article in Special Issue
Fostering Holistic Development with a Designed Multisport Intervention in Physical Education: A Class-Randomized Cross-Over Trial
Previous Article in Journal
Developmental Hip Dysplasia: An Epidemiological Nationwide Study in Italy from 2001 to 2016
Previous Article in Special Issue
Children’s Active School Travel: Examining the Combined Perceived and Objective Built-Environment Factors from Space Syntax
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 18
Issue 12
10.3390/ijerph18126588
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Urban Green Spaces, Greenness Exposure and Species Richness in Residential Environments and Relations with Physical Activity and BMI in Portuguese Adolescents

by
Juliana Melo
1,*,
Ana Isabel Ribeiro
2,3,
Susana Aznar
4,5,
Andreia Pizarro
1,3 and
Maria Paula Santos
1,2
1
Research Center in Physical Activity, Health and Leisure (CIAFEL), Faculty of Sports, University of Porto (FADEUP), Rua Dr. Placido Costa, 91, 4200-450 Porto, Portugal
2
EPIUnit, Instituto de Saúde Pública da Universidade do Porto, Rua das Taipas 135, 4050-600 Porto, Portugal
3
Laboratory for Integrative and Translational Research in Population Health (ITR), Rua das Taipas 135, 4050-600 Porto, Portugal
4
PAFS Research Group, Faculty of Sports Sciences, University of Castilla-La Mancha (UCLM), Avda. Carlos III s/n, 45071 Toledo, Spain
5
Biomedical Research Networking Center on Frailty and Healthy Aging (CIBERFES), 28029 Madrid, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(12), 6588; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126588
Submission received: 1 April 2021 / Revised: 8 June 2021 / Accepted: 16 June 2021 / Published: 18 June 2021
(This article belongs to the Special Issue Psychosocial and Environmental Correlates of Physical Activity in Youth and Childhood)
Versions Notes

Abstract

:
Environmental factors play an important role in obesity-related behaviors. Evidence indicates significant associations between weight and urban green spaces in adults, but it is not clear whether this relationship applies to adolescents. Therefore, our aim was to determine the associations between urban green spaces, greenness exposure and species richness in residential environments with physical activity and body mass index. Sixty-two adolescents between 12 and 18 years of age answered a self-administered questionnaire, providing information on height, weight, age, sex and home address. Data on socioeconomic deprivation were obtained from the European Index of Deprivation for Small Portuguese Areas. Physical activity levels were assessed using accelerometers. Urban green space counts and the normalized difference vegetation index values were measured using buffers along the roads with distances of 300, 500, 1000 and 1500 m from each participant’s residence. To quantify the species richness, the species richness index was used. Linear regression models were fitted to analyze whether urban green spaces, exposure to green spaces and species richness counts for each distance were associated with physical activity and self-reported body mass index. We did not find significant associations between the independent variables and the probability of overweight or obesity. The relationship between environmental variables, adolescents’ physical activity and body weight seems to be complex and further studies may contribute to better understanding of the topic.

1. Introduction

In recent decades, obesity and overweight have shown rapid increases and high prevalence in children and adolescents. This scenario has amplified the concern for the health of young people worldwide and has motivated researches to clarify causal factors of obesity for these age groups [1].
A high body mass index (BMI) is strongly associated with an increased risk of developing chronic diseases, such as cardiovascular diseases [2], diabetes mellitus [3], cancer [4] and some musculoskeletal disorders [5]. In addition to a greater risk of developing chronic non-communicable diseases, overweight or obesity during childhood tends to persist into adulthood [6,7].
Several studies have shown that children and adolescents who engage in more physical activity (PA) have a lower BMI than those who are less active [8,9]. In this context, physical inactivity and sedentary behaviors are increasingly being identified as determining factors for pediatric overweight or obesity [10,11].
In this context, socioecological models suggest that individual, interpersonal, environmental and political factors influence health and PA behaviors [12]. The main theory of socioecological models is that influences at various levels can interact and exert synergistic effects on behaviors, and in this sense interventions at various levels may be more effective for changing behaviors [13]. Environmental characteristics, such as access to recreation facilities, aesthetics and infrastructure for pedestrians, seem to consistently predict general and recreational PA, while neighborhood mobility and street connectivity seem to be consistently related to PA for transportation [14].
Recent studies have demonstrated that certain characteristics of the built environment, such as street connectivity and transport systems, are capable of promoting PA; therefore, environmental interventions seem to be a sustainable strategy to encourage young populations to engage in or increase PA [15]. Moreover, the influence of urban environment seems to be a good strategy when designing future community-based interventions to promote more active lifestyles among adolescents, particularly when adolescents’ self-efficacy and moderate-to-vigorous physical activity (MVPA) are low [16].
As the existing studies did not simultaneously consider the impacts of all the characteristics of the built environment, it is still difficult to indicate the characteristics that most affect the levels of PA of adolescents [17]. Regarding the perceived barriers, it is known that a lack of motivation [18], lack of time [19], lack of knowledge of the existing PA facilities or program, lack of transportation, concerns with neighborhood security and lack of access to recreational facilities are often reported [20]. Regarding studies that have used measures of the built environment, most of them have reported inconsistent associations with street connectivity and good walkability with PA but consistent relationships between the availability of recreational facilities and mixed land use [21].
Among recreational facilities, urban green spaces (UGS) in particular have been attracting attention [22]. This type of space is an easily modifiable element for the promotion of PA and the reduction of sedentary behaviors.
The relationship between PA and BMI is well documented and the presence of UGS in the residential environment has been increasingly described as one of the determinants highly associated with overweight and obesity [8,9].
The negative relationship between the number of UGS in the residential area BMI has been consistently reported among child populations [23,24]; however, few studies have investigated this connection among adolescents and there is insufficient evidence to make consistent associations [25].
Although many studies have demonstrated positive impacts of UGS on human health and well-being [26,27], it is not clear whether solutions based on an increase in urban nature can really contribute positively to increase PA levels and consequently reduce in overweight and obesity levels [28,29].
Given this scenario, green space measures, such as exposure to green spaces and biodiversity, can be assessed to help understand and comprehensively clarify the possible relationships between UGS and health [30].
We aimed to explore whether the number of UGS, the level of green exposure and the amount of species richness in the residential environment were associated with PA and BMI in Portuguese adolescents, and whether these relationships change with adjustment to the socioeconomic deprivation index.

2. Materials and Methods

2.1. Study Area

Oporto is the second city and the fourth most populous municipality in Portugal. The municipality, with an area of 41.42 km2 [31], has a population of around 237,000 inhabitants within its administrative limits, being subdivided into seven parishes [32].

2.2. Study Design

Data were collected in 4 of the 7 public high schools in Oporto city. The schools were selected for convenience, according to the interest of teachers and school managers in participating in the study. Participation was requested from all adolescents enrolled in physical education classes and whose teachers agreed to participate in the research.

2.3. Sample

The calculations for determining the sample size with a confidence level of 95% and a margin of error of 5% resulted in a total of 382 subjects [33].
We obtained a total of 311 participants residing in the metropolitan area of Porto; however, since data referring to the species richness index were only available for the city of Porto, 249 subjects were excluded because they lived outside the limits of this municipality. The final sample size comprised 62 adolescents (32 girls, 30 boys) between 12 and 18 years of age. Prior to data collection, both participants and their parents received a complete explanation of the purpose of the investigation and signed a written consent form to participate in the study.
Ethical approval for this study was obtained from the Ethic Committee of the Faculty of Sports of the University of Porto (CEFADE 22.2013). Consent was also collected from the Education Authority and the School Board.

2.4. Individual Anthropometric and Sociodemographic Data

Information on height, weight, age, gender and residential address were obtained through a self-administered questionnaire.
The measure commonly used to define the weight status is the BMI and associations between BMI and health outcomes within and between populations are often used to help determine possible causes of illness [34].
The BMI values were calculated as weight in kilograms divided by height in meters squared (kg/m2), while the BMI z-scores were categorized into normal weight (1 standard deviation (SD)), overweight (>+1SD) and obesity (>+2SD) [34]. For analysis purposes, the overweight and obesity categories were grouped into just one category.

2.5. Socioeconomic Deprivation

Data regarding the socioeconomic deprivation of the sample subjects were obtained from the European Index of Deprivation for Small Portuguese Areas (EDI-PT) [35].
The EDI-PT is a sensitive measure to capture health inequalities, as it is associated with the development and health of the population [35]. The present variable was categorized into terciles (1 for least private and 3 for more private), which were then classified as “low”, “medium” and “high” deprivation, respectively. This strategy aimed to maximize the heterogeneity of results [36].

2.6. Physical Activity

Physical activity (PA) levels were assessed using the Actigraph model GT3X+ accelerometers (Actigraph, Pensacola, FL, USA). Participants used an accelerometer attached to an elastic belt positioned on the right side of the pelvic girdle throughout waking hours for seven consecutive days, except during bathing and water activities.
The epoch was set to 30 s, and for data analysis purposes only participants with a minimum of 600 min of daily use for 4 weekdays and 1 weekend day (or with 3000 min of use over 4 days) were included [37].
The cut-off points proposed by Evenson [38] were used to categorize moderate-to-vigorous PA (MVPA).

2.7. Assessment of the Geographical Accessibility of Urban Green Spaces

A total of 226 public green spaces were considered, with no restrictions on size, location or characteristics, covering the universe of public green spaces available in the study area that could be freely (free of charge) used by the population. Polygons of green spaces and entrance places were obtained from digital city maps, as described in previous studies [39].
Using the location of each child’s residence as a starting point, the following measures based on a network of geographic accessibility to green spaces were calculated: number of available green spaces within 300, 500, 1000 and 1500 m from the residence and school (score). When the green spaces were delimited (that is, with fences), we used the distance to the entrance; otherwise, the distance to the limit was used.
For the calculations, we used ArcGIS software version 10.5 and the Network Analyst extension, using an updated street network dataset provided by the Environmental Systems Research Institute.

2.8. Normalized Difference Vegetation Index

The normalized difference vegetation index (NDVI) was used to analyze the greenness exposure [40]. To capture the surrounding greenery, we calculated the average NDVI using buffers along the road with distances at 300, 500, 1000 and 1500 m from each participant’s residence. The NDVI is calculated based on the ground surface reflectance of the visible red (VISR) and near-infrared (NIR) wavelengths using the formula presented in Equation (1):
NDVI = (NIR − VISR)/(NIR + VIRSR)
The underlying principle employed in the NDVI calculation is that chlorophyll in healthy vegetation absorbs radiation in the VIRS region (630–690 nm) of the electromagnetic spectrum and reflects radiation in the NIR region (760–900 nm). It is a unitless index ranging from −1 to 1, with higher values corresponding to a higher density of healthy vegetation. For this study, images with 5% or less of Landsat 8 cloud cover (spatial resolution: 30 m) were used during the spring–summer (vegetation peak) period of 2016–2017 [41].
ArcMap 10.5 was used to process satellite images and QGIS 3.8.2 was used to extract the average NDVI for each participant.

2.9. Species Richness Assessment

In order to quantify the species richness, we used the species richness index (SRI). This index was obtained from a species richness report previously compiled by Porto City Council and the Research Center for Biodiversity and Genetic Resources (CIBIO) [42,43]. The number, type and area of habitats in Porto were retrieved, followed by an estimation of the number of species living in those habitats. Included species consisted of four vertebrate groups (amphibians, birds, reptiles and small mammals), totaling 89 different species.
To express the concentrations of species in the habitats that surrounded the residences, we summed the numbers of different species present in the habitats located within different network-based buffers (300, 500, 1000 and 1500 m) multiplied by the total area of each habitat, then subsequently we divided the resulting values by the buffer area in square meters.
Such an approach has been used in previous studies, such as the EXALAR XXI project (POCI-01-0145-FEDER-030193; https://www.exalar21.com/home|accessed on 13 November 2019) [44,45].

2.10. Statistical Analysis

Adjusted linear regression models were used to analyze whether UGS counts, exposure to green spaces (NDVI) and SRI for road distances of 300, 500, 1000 and 1500 m were associated with PA and BMI. The models were applied to the entire sample and by gender. The EDI-PT of the residential neighborhood was included in the model as an adjustment variable.
The likelihood ratio tests were performed to compare the power of a predictive model. Collinearity diagnostics were run on the variables included in each model. Odds ratios were calculated for 95% confidence intervals.
The intraclass correlation coefficient (ICC) was also calculated to provide a better understanding of the proportion of variance due to group differences [46].
All statistical analyses were performed using Statistical Package for the Social Sciences (SPSS) version 25 (IBM Corp, Armonk, NY, USA).

3. Results

The sample characteristics are presented in Table 1. About 51% of the students were girls. The sample deprivation index values were predominantly (71.1%) medium and high. The BMI average for girls was 21.7 ± 3.6 and for boys was 22.2 ± 2.9 kg/m2. About 83% of the adolescents did not achieve daily PA recommendations of 60 min of MVPA, with the average MVPA being about 41 min.
Descriptive statistics for each variable by gender and by BMI categories are shown in Table 2. About 23% (n = 14) of the sample was overweight or obese. The average BMI for the normal weight group was approximately 20 kg/m2 for both sexes.
In relation to the UGS average at 300 m, the residential environment for boys with normal weight contained fewer green spaces when compared with the environment for boys with overweight / obesity (p = 0.04).
In this study, the ICC value (ICC = 0.08) suggested that clustering between neighborhoods and schools was non-existent; therefore, a multilevel analysis was not used [46].
The crude regression models performed individually with each predictor (Tables S1–S3, S6–S8) did not show significant results. Models including all predictors simultaneously did not show significant results (see Table 3). In addition, we tested the model stratified by gender (see Tables S4, S5, S9 and S10) and found that the socioeconomic deprivation index was inversely correlated to BMI in boys (r = −0.4), but the models result did not show statistical significance.
Since socioeconomic deprivation seems to play an important role in the relationship with overweight and obesity, we developed a model with adjustment for the socioeconomic deprivation index (see Table 4); however, the results were the same, showing no significant statistics.
Likelihood ratio tests showed that the addition of variables to each subsequent model did not significantly improve the model fit. Collinearity diagnostics did not reveal major concerns, with no variance inflation factor scores above three.

4. Discussion

An increasing number of studies have sought to assess the relationship between urban nature and human health, using many different research approaches and obtaining heterogeneous results.
In our study, the socioeconomic deprivation index was inversely related to BMI in boys and no significant results were found for the other variables. Potential confounding variables were included but the size of the effect did not change. Regardless, our findings raise some interesting points for discussion.
There is a general consensus that green spaces close to residential areas support human health and that their distance to a subject’s home is important [47].
The European Environment Agency (EEA) [48] recommends a walk of up to 15 min to the nearest UGS. Nielsen and Hansen [49] claim that a distance of 300 to 400 m is seen as a limit value, after which the frequency of use begins to decrease. However, Kaczynski, Besenyi [50] demonstrated that the existing recommendations do not appear to be relevant for adolescent populations. Ding and Gebel [51] corroborated this hypothesis when reporting that the proximity of the UGS was not associated with PA in adolescents; in addition, they concluded that the use of objective measures increased the consistency of the results.
In this context, McCrorie et al. [52] stated that part of the recent growth in the PA environmental determinants can be attributed to the advancement of technology. The use of objective approaches enables accurate measurements, allowing a more detailed representation of the environment.
In particular, the objective measures of PA, accessibility to the UGS and exposure to green spaces seem to be very relevant with regard to the identification of environmental determinants about the use of UGS and the relationship with the health of populations [52,53,54].
Except for the weight and height obtained by self-report, in our study we mainly used information from objective measures, with the latter method seeming to be more advantageous compared to the subjective measures; thus, we consider that our results are likely to be quite reliable, although we emphasize that the environmental determinants related to UGS in favor of health seem complex and interrelated [55].
In addition, despite the unanimity of the benefits of UGS for health populations, certain environmental determinants of the use and provision of UGS may not be universal, since each city has a specific structure and a peculiar group of inhabitants [56].
Although UGS have generally been treated as a type of homogeneous environment, some studies have separated the spaces into relatively broad typologies, suggesting that some types affect health to varying degrees; that is, not all UGS are equal in their health benefits [57].
The composition and configuration of urban green spaces are usually measured through metrics of the landscape [11]. Among the dozens of existing methods, the NDVI has been gaining prominence in studies that analyze green spaces and their relationship with health [58].
Although some studies have postulated that more green closer to home may be important for health [59], in our study we did not find any significant associations between NDVI, PA and BMI in adolescents.
Recently Ekkel and de Vries [60] conducted a review study on green spaces and health, showing that associations vary within and between studies, mainly due to the different buffer sizes, a factor that probably ends up leading to an inconsistent direction of the effects.
Eventually, if similar approaches are replicated elsewhere, they may be able to provide more information about the mechanisms by which exposure to green spaces may or may not influence PA and BMI in juvenile populations.
In this context, the relationship between UGS and health is becoming prominent in the urban planning agenda and in development policies, which already refer to the positive effects of the presence, use and design of UGS [29,56].
Despite the UGS being integrated into policies and urban planning development, the role of these spaces in supporting biodiversity and the links between biodiversity and human health have received insufficient attention, with the few existing studies on these subjects presenting inconclusive results [61].
This study was one of the first to explore the influence of species richness on adolescents’ PA and BMI, for which we found no significant results.
Interestingly, the results from the study by Schwartz et al. [62] highlighted that the loss of urban biodiversity can result in compromised ecosystem functions, which in turn can negatively influence human health.
In contrast, the study by Leemans and De Groot [63] concluded that a decline in ecosystems was accompanied by steady gains in human well-being at the global level. Using this logic, Hough [64] explained that this discovery seems intuitive, since economic or industrial development tends to threaten species and at the same time increases human life expectancy.
However, Mills et al. [65] found that industrialized urban habitats are poor in biodiversity, which ends up discouraging contact with beneficial environmental microbiotics for the human body. When associated with inadequate diets, excessive intake of antibiotics, physical inactivity and other factors seem to be directly linked to the epidemic of non-communicable diseases in societies with these characteristics.
In addition to the interaction between UGS and biodiversity, the interaction of SES with the environmental context has also gained some relevance. Our results showed the existence of an inverse relationship between the socioeconomic deprivation index and the boys’ BMI values. Some studies conducted with children and adolescents have highlighted the fact that areas with greater material deprivation or lower levels of wealth have higher prevalence of obesity or overweight, regardless of the socioeconomic position of the household [66,67].
However, Gordon-Larsen et al. [68] disagreed when reporting that the amount of physical activity can be linked to economic factors, whereby lower socioeconomic areas in many cases are built with fewer and worse leisure areas and also involve greater distances to reach the spaces that offer PA opportunities. Furthermore, neighborhoods can be perceived as more or less safe, and consequently more or less suitable for leisure activities for young people.
According to Noonan [69], there are complex paths that link socioeconomic disadvantages with PA and overweight or childhood obesity. The disparity between the methods used in the analysis of socioeconomic deprivation could be a relevant factor when it comes to obtaining harmonious conclusions. Evidence on this topic suggests that these socioeconomic disparities in health are the result of differentiated access to built environmental and economic resources [70]; however, until now, the analysis of socioeconomic differences in terms of amounts of PA and obesity trends in adolescents has been limited [71]. Additionally, it does not seem there is a general consensus on what constitutes socioeconomic deprivation [72].
The method used for analyzing socioeconomic deprivation could be an important factor in terms of inconsistencies, since there are several existing measures, including the education of parents or guardians and household income.
Although this study does not report significant results related to PA and overweight or obesity, further studies may be important to improve knowledge in the search for an understanding of the mechanisms by which access to green spaces, greenery and the presence of species richness in urban environments influence PA and weight in adolescent populations.

5. Study Limitations

Studies of the determinants of urban health are complex because cities are constantly changing, resulting in different living conditions within and between cities. In this sense, the availability of green spaces can vary considerably between different urban areas. Currently, there are no standards that accurately establish the quantity or ideal characteristics of green spaces associated with benefits for human health.
Although this study tried to fill a knowledge gap, certain limitations need to be recognized. First, since it was a cross-sectional study, it was not possible to make an inference about causal directions. Second, although self-reported anthropometrics are valuable sources of data, the self-reported bias can have important consequences for the accuracy of screening for overweight and obesity [73]. Previous studies have shown that self-reported weights are significantly lower than the measured weights [74,75]; however, the use of BMI resulting from self-reported data shows good performance, moderate sensitivity and high specificity. With that said, this appears to be an alternative when direct BMI measurement is not available [76].
Third, we had no information about possible confounding variables, such as concerns about safety, traffic or related to overall aesthetics environments’, which made impossible testing their effects. In addition, it is likely that there are other attributes of green spaces (e.g., quality, size) that we did not capture that could influence the levels of adolescent PA, and which consequently could impact BMI values.
Fourth, despite being considered an a priori estimate of the ideal sample size for this research, the smaller geographical extent of the species richness dataset forced us to conduct the study using a smaller sample. This may have influenced the results, especially with regard to the probability of detecting statistical significant associations between variables [77]. Furthermore, determining the sample size when applying multilevel models is important, as the statistical power depends on the total sample size for each level [78]. In this context, considering that statistical power and p values depend on both the effect size and the sample size, the sample used probably did not provide enough power to find statistical differences [77].
Despite these weaknesses and the lack of evidence on the relationships between the amount of UGS, green exposure and the amount of species richness inserted in the residential environment with PA and BMI for Portuguese adolescents, the results of this study are a relevant contribution to the literature as they provide are groundbreaking.

6. Conclusions

In this study, we did not find relations between the number of urban green spaces, green exposure and species richness in residence areas with PA and incidence rates of overweight and obesity among a Portuguese adolescent population.
The relationships between the amount of UGS, exposure to green spaces and the richness of species with PA and BMI seem to be complex, and further studies are needed for a better understanding of this topic.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijerph18126588/s1.

Author Contributions

Conceptualization J.M., A.I.R., A.P. and M.P.S.; methodology, J.M., A.I.R. and A.P.; data collection, J.M. and A.P.; analysis, J.M., A.I.R. and A.P.; writing—original draft preparation, J.M., A.P. and M.P.S.; writing—review and editing, J.M., A.P., S.A. and M.P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FEDER through the Operational Program for Competitiveness and Internationalization and received national funding from the Foundation for Science and Technology—FCT (Portuguese Ministry of Science, Technology and Higher Education), under the Research Center in Physical Activity, Health and Leisure (CIAFEL) (UIDB/00617/2020) and also under the Epidemiology Research Unit (EPIUnit)—Institute of Public Health of the University of Porto (UIDB/04750/2020). Ana Isabel Ribeiro was supported by national funds through FCT, under the ‘Stimulus of Scientific Employment—Individual Support’ program within the contract CEECIND/02386/2018.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethic Committee of the Faculty of Sports of the University of Porto (process code CEFADE 22.2013, approved on 21 January 2014).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data supporting the findings of this study are available within the article and its supplementary materials.

Acknowledgments

The authors would like to thank the research volunteers, as well as the Editorial Office of IJERPH and all anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jastreboff, A.M.; Kotz, C.M.; Kahan, S.; Kelly, A.S.; Heymsfield, S.B. Obesity as a Disease: The Obesity Society 2018 Position Statement. Obesity 2019, 27, 7–9. [Google Scholar] [CrossRef] [Green Version]
  2. Ortega, F.B.; Lavie, C.J.; Blair, S.N. Obesity and cardiovascular disease. Circ. Res. 2016, 118, 1752–1770. [Google Scholar] [CrossRef] [Green Version]
  3. Singh, G.K.; Siahpush, M.; Kogan, M.D. Neighborhood socioeconomic conditions, built environments, and childhood obesity. Health Aff. 2010, 29, 503–512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lauby-Secretan, B.; Scoccianti, C.; Loomis, D.; Grosse, Y.; Bianchini, F.; Straif, K.; International Agency for research on cancer handbook working group. Body fatness and cancer—Viewpoint of the IARC working group. N. Engl. J. Med. 2016, 375, 794–798. [Google Scholar] [CrossRef] [Green Version]
  5. Molina-Garcia, P.; Migueles, J.H.; Cadenas-Sanchez, C.; Esteban-Cornejo, I.; Mora-Gonzalez, J.; Rodriguez-Ayllon, M.; Plaza-Florido, A.; Vanrenterghem, J.; Ortega, F.B. A systematic review on biomechanical characteristics of walking in children and adolescents with overweight/obesity: Possible implications for the development of musculoskeletal disorders. Obes. Rev. 2019, 20, 1033–1044. [Google Scholar] [CrossRef] [PubMed]
  6. Ryder, J.R.; Jacobs, D.R.; Sinaiko, A.R.; Kornblum, A.P.; Steinberger, J. Longitudinal changes in weight status from childhood and adolescence to adulthood. J. Pediatr. 2019, 214, 187–192. [Google Scholar] [CrossRef] [PubMed]
  7. World Health Organization. Diet, Nutrition, and the Prevention of Chronic Diseases: Report of a Joint WHO/FAO Expert Consultation; World Health Organization: Geneva, Switzerland, 2003; Volume 916. [Google Scholar]
  8. James, P.; Hart, J.E.; Arcaya, M.C.; Feskanich, D.; Laden, F.; Subramanian, S. Neighborhood Self-Selection: The Role of Pre-Move Health Factors on the Built and Socioeconomic Environment. Int. J. Environ. Res. Public Health 2015, 12, 12489–12504. [Google Scholar] [CrossRef] [Green Version]
  9. Smith, L.; Foley, L.; Panter, J. Activity spaces in studies of the environment and physical activity: A review and synthesis of implications for causality. Health Place 2019, 58, 102113. [Google Scholar] [CrossRef] [PubMed]
  10. Kenney, E.L.; Gortmaker, S.L. United States adolescents’ television, computer, videogame, smartphone, and tablet use: Associations with sugary drinks, sleep, physical activity, and obesity. J. Pediatr. 2017, 182, 144–149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Myers, A.; Gibbons, C.; Finlayson, G.; Blundell, J. Associations among sedentary and active behaviours, body fat and appetite dysregulation: Investigating the myth of physical inactivity and obesity. Br. J. Sports Med. 2017, 51, 1540–1544. [Google Scholar] [CrossRef]
  12. Sallis, J.F.; Cervero, R.B.; Ascher, W.L.; Henderson, K.A.; Kraft, M.K.; Kerr, J. An ecological approach to creating active living communities. Annu. Rev. Public Health 2006, 27, 297–322. [Google Scholar] [CrossRef] [Green Version]
  13. Sallis, J.F.; Owen, N.; Fisher, E. Ecological models of health behavior. Health Behav. Theory Res. Pract. 2015, 5, 43–64. [Google Scholar]
  14. Bauman, A.E.; Reis, R.S.; Sallis, J.F.; Wells, J.; Loos, R.; Martin, B.W. Correlates of physical activity: Why are some people physically active and others not? Lancet 2012, 380, 258–271. [Google Scholar] [CrossRef]
  15. Andersen, H.B.; Christiansen, L.B.; Klinker, C.D.; Ersbøll, A.K.; Troelsen, J.; Kerr, J.; Schipperijn, J. Increases in use and activity due to urban renewal: Effect of a natural experiment. Am. J. Prev. Med. 2017, 53, e81–e87. [Google Scholar] [CrossRef]
  16. Aznar, S.; Queralt, A.; García-Massó, X.; Villarrasa-Sapiña, I.; Molina-García, J. Multifactorial combinations predicting active vs inactive stages of change for physical activity in adolescents considering built environment and psychosocial factors: A classification tree approach. Health Place 2018, 53, 150–154. [Google Scholar] [CrossRef] [PubMed]
  17. Ding, D.; Sallis, J.F.; Kerr, J.; Lee, S.; Rosenberg, D.E. Neighborhood environment and physical activity among youth: A review. Am. J. Prev. Med. 2011, 41, 442–455. [Google Scholar] [CrossRef] [PubMed]
  18. Kubik, M.Y.; Lytle, L.; Fulkerson, J.A. Fruits, vegetables, and football: Findings from focus groups with alternative high school students regarding eating and physical activity. J. Adolesc. Health 2005, 36, 494–500. [Google Scholar] [CrossRef] [PubMed]
  19. Hannay, J.; Dudley, R.; Milan, S.; Leibovitz, P.K. Combining Photovoice and focus groups: Engaging latina teens in community assessment. Am. J. Prev. Med. 2013, 44, S215–S224. [Google Scholar] [CrossRef] [Green Version]
  20. Finkelstein, D.M.; Petersen, D.M.; Schottenfeld, L.S. Promoting children’s physical activity in low-income communities in colorado: What are the barriers and opportunities? Prev. Chronic Dis. 2017, 14, E134. [Google Scholar] [CrossRef]
  21. Carlson, J.A.; Saelens, B.E.; Kerr, J.; Schipperijn, J.; Conway, T.L.; Frank, L.D.; Chapman, J.E.; Glanz, K.; Cain, K.L.; Sallis, J.F. Association between neighborhood walkability and GPS-measured walking, bicycling and vehicle time in adolescents. Health Place 2015, 32, 1–7. [Google Scholar] [CrossRef] [Green Version]
  22. An, R.; Shen, J.; Yang, Q.; Yang, Y. Impact of built environment on physical activity and obesity among children and adolescents in China: A narrative systematic review. J. Sport Health Sci. 2019, 8, 153–169. [Google Scholar] [CrossRef]
  23. Sanders, T.; Feng, X.; Fahey, P.P.; Lonsdale, C.; Astell-Burt, T. Greener neighbourhoods, slimmer children? Evidence from 4423 participants aged 6 to 13 years in the Longitudinal Study of Australian children. Int. J. Obes. 2015, 39, 1224–1229. [Google Scholar] [CrossRef] [PubMed]
  24. Peirson, L.; Fitzpatrick-Lewis, D.; Morrison, K.M.; Warren, R.; Ali, M.; Raina, P. Treatment of overweight and obesity in children and youth: A systematic review and meta-analysis. CMAJ Open 2015, 3, E35–E46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Lachowycz, K.; Jones, A.P. Greenspace and obesity: A systematic review of the evidence. Obes. Rev. 2011, 12, e183–e189. [Google Scholar] [CrossRef] [PubMed]
  26. Kabisch, N.; van den Bosch, M.; Lafortezza, R. The health benefits of nature-based solutions to urbanization challenges for children and the elderly—A systematic review. Environ. Res. 2017, 159, 362–373. [Google Scholar] [CrossRef] [PubMed]
  27. Sandifer, P.A.; Sutton-Grier, A.E.; Ward, B.P. Exploring connections among nature, biodiversity, ecosystem services, and human health and well-being: Opportunities to enhance health and biodiversity conservation. Ecosyst. Serv. 2015, 12, 1–15. [Google Scholar] [CrossRef] [Green Version]
  28. Dallimer, M.; Irvine, K.N.; Skinner, A.M.J.; Davies, Z.G.; Rouquette, J.R.; Maltby, L.; Warren, P.H.; Armsworth, P.R.; Gaston, K.J. Biodiversity and the Feel-good factor: Understanding associations between self-reported human well-being and species richness. BioScience 2012, 62, 47–55. [Google Scholar] [CrossRef] [Green Version]
  29. Douglas, O.; Lennon, M.; Scott, M. Green space benefits for health and well-being: A life-course approach for urban planning, design and management. Cities 2017, 66, 53–62. [Google Scholar] [CrossRef] [Green Version]
  30. Lai, H.; Flies, E.J.; Weinstein, P.; Woodward, A. The impact of green space and biodiversity on health. Front. Ecol. Environ. 2019, 17, 383–390. [Google Scholar] [CrossRef]
  31. Instituto Português de Geografia. Carta Administrativa Oficial de Portugal (CAOP); Instituto Português de Geografia: Lisbon, Portugal, 2013.
  32. Instituto Nacional de Estatística. Censos 2011—Quadros de Apuramento por Freguesia; Instituto Nacional de Estatística: Lisbon, Portugal, 2012.
  33. Taherdoost, H. Determining sample size; how to calculate survey sample size. Int. J. Econ. Manag. Syst. 2017, 2, 3. [Google Scholar]
  34. World Health Organization. Physical status: The Use of and Interpretation of Anthropometry, Report of a WHO Expert Committee; WHO: Geneva, Switzerland, 1995. [Google Scholar]
  35. Ribeiro, A.I.; Launay, L.; Guillaume, E.; Launoy, G.; Barros, H. The Portuguese version of the European Deprivation Index: Development and association with all-cause mortality. PLoS ONE 2018, 13, e0208320. [Google Scholar] [CrossRef] [Green Version]
  36. Pigeyre, M.; Rousseaux, J.; Trouiller, P.; Dumont, J.; Goumidi, L.; Bonte, D.; Dumont, M.-P.; Chmielewski, A.; Duhamel, A.; Amouyel, P.; et al. How obesity relates to socio-economic status: Identification of eating behavior mediators. Int. J. Obes. 2016, 40, 1794–1801. [Google Scholar] [CrossRef]
  37. Migueles, J.H.; Cadenas-Sanchez, C.; Ekelund, U.; Nyström, C.D.; Mora-Gonzalez, J.; Löf, M.; Labayen, I.; Ruiz, J.R.; Ortega, F.B. Accelerometer data collection and processing criteria to assess physical activity and other outcomes: A systematic review and practical considerations. Sports Med. 2017, 47, 1821–1845. [Google Scholar] [CrossRef] [PubMed]
  38. Evenson, K.R.; Catellier, D.J.; Gill, K.; Ondrak, K.S.; McMurray, R.G. Calibration of two objective measures of physical activity for children. J. Sports Sci. 2008, 26, 1557–1565. [Google Scholar] [CrossRef] [PubMed]
  39. Hoffimann, E.; Barros, H.; Ribeiro, A.I. Socioeconomic Inequalities in Green Space Quality and Accessibility—Evidence from a Southern European City. Int. J. Environ. Res. Public Health 2017, 14, 916. [Google Scholar] [CrossRef] [PubMed]
  40. Rufo, J.C.; Paciência, I.; Ribeiro, A.I. Green environments and allergic diseases in children: A scoping review. Curr. Epidemiol. Rep. 2019, 6, 442–448. [Google Scholar] [CrossRef]
  41. Paciência, I.; Moreira, A.; Moreira, C.; Rufo, J.C.; Sokhatska, O.; Rama, T.; Hoffimann, E.; Santos, A.C.; Barros, H.; Ribeiro, A.I. Neighbourhood green and blue spaces and allergic sensitization in children: A longitudinal study based on repeated measures from the Generation XXI cohort. Sci. Total Environ. 2021, 772, 145394. [Google Scholar] [CrossRef]
  42. Gomes, V.; Ribeiro, R.; Carretero, M. Effects of urban habitat fragmentation on common small mammals: Species versus communities. Biodivers. Conserv. 2011, 20, 3577–3590. [Google Scholar] [CrossRef]
  43. Guilherme, F.; Carretero, M.A.; Farinha-Marques, P. Vertebrates and invertebrates of European cities: Selected non-avian fauna. In Vertebrates and Invertebrates of European Cities; Kelcey, J.G., Ed.; Springer: New York, NY, USA, 2015. [Google Scholar]
  44. Rufo, J.C.; Paciência, I.; Hoffimann, E.; Moreira, A.; Barros, H.; Ribeiro, A.I. The neighbourhood natural environment is associated with asthma in children: A birth cohort study. Allergy 2021, 76, 348–358. [Google Scholar] [CrossRef]
  45. Rufo, J.C.; Ribeiro, A.I.; Paciência, I.R.; Delgado, L.; Moreira, A.M.A. The influence of species richness in primary school surroundings on children lung function and allergic disease development. Pediatr. Allergy Immunol. 2020, 31, 358–363. [Google Scholar] [CrossRef]
  46. Merlo, J.; Chaix, B.; Ohlsson, H.; Beckman, A.; Johnell, K.; Hjerpe, P.; Råstam, L.; Larsen, K. A brief conceptual tutorial of multilevel analysis in social epidemiology: Using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J. Epidemiology Community Health 2006, 60, 290–297. [Google Scholar] [CrossRef] [Green Version]
  47. Sugiyama, T.; Cerin, E.; Owen, N.; Oyeyemi, A.L.; Conway, T.L.; Van Dyck, D.; Schipperijn, J.; Macfarlane, D.J.; Salvo, D.; Reis, R.S.; et al. Perceived neighbourhood environmental attributes associated with adults’ recreational walking: IPEN Adult study in 12 countries. Health Place 2014, 28, 22–30. [Google Scholar] [CrossRef] [Green Version]
  48. Stanners, D.; Bourdeau, P. Europe’s Environment: The Dobris Assessment; European Environment Agency: Copenhagen, Danmark, 1995.
  49. Nielsen, T.S.; Hansen, K.B. Do green areas affect health? Results from a Danish survey on the use of green areas and health indicators. Health Place 2007, 13, 839–850. [Google Scholar] [CrossRef]
  50. Kaczynski, A.T.; Besenyi, G.M.; Stanis, S.A.W.; Koohsari, M.J.; Oestman, K.B.; Bergstrom, R.D.; Potwarka, L.R.; Reis, R.S. Are park proximity and park features related to park use and park-based physical activity among adults? Variations by multiple socio-demographic characteristics. Int. J. Behav. Nutr. Phys. Act. 2014, 11, 146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Ding, D.; Gebel, K. Built environment, physical activity, and obesity: What have we learned from reviewing the literature? Health Place 2012, 18, 100–105. [Google Scholar] [CrossRef] [Green Version]
  52. McCrorie, P.R.; Fenton, C.; Ellaway, A. Combining GPS, GIS, and accelerometry to explore the physical activity and environment relationship in children and young people—A review. Int. J. Behav. Nutr. Phys. Act. 2014, 11, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Klompmaker, J.O.; Hoek, G.; Bloemsma, L.D.; Gehring, U.; Strak, M.; Wijga, A.H.; Brink, C.V.D.; Brunekreef, B.; Lebret, E.; Janssen, N.A. Green space definition affects associations of green space with overweight and physical activity. Environ. Res. 2018, 160, 531–540. [Google Scholar] [CrossRef] [PubMed]
  54. Hunter, R.F.; Christian, H.; Veitch, J.; Astell-Burt, T.; Hipp, J.; Schipperijn, J. The impact of interventions to promote physical activity in urban green space: A systematic review and recommendations for future research. Soc. Sci. Med. 2015, 124, 246–256. [Google Scholar] [CrossRef] [Green Version]
  55. Markevych, I.; Schoierer, J.; Hartig, T.; Chudnovsky, A.; Hystad, P.; Dzhambov, A.; de Vries, S.; Triguero-Mas, M.; Brauer, M.; Nieuwenhuijsen, M.J.; et al. Exploring pathways linking greenspace to health: Theoretical and methodological guidance. Environ. Res. 2017, 158, 301–317. [Google Scholar] [CrossRef] [PubMed]
  56. Schipperijn, J.; Stigsdotter, U.K.; Randrup, T.B.; Troelsen, J. Influences on the use of urban green space—A case study in Odense, Denmark. Urban For. Urban Green. 2010, 9, 25–32. [Google Scholar] [CrossRef]
  57. James, P.; Banay, R.F.; Hart, J.E.; Laden, F. A Review of the health benefits of greenness. Curr. Epidemiol. Rep. 2015, 2, 131–142. [Google Scholar] [CrossRef] [Green Version]
  58. Mueller, W.; Steinle, S.; Pärkkä, J.; Parmes, E.; Liedes, H.; Kuijpers, E.; Pronk, A.; Sarigiannis, D.; Karakitsios, S.; Chapizanis, D.; et al. Urban greenspace and the indoor environment: Pathways to health via indoor particulate matter, noise, and road noise annoyance. Environ. Res. 2020, 180, 108850. [Google Scholar] [CrossRef]
  59. Annerstedt, M.; Östergren, P.-O.; Björk, J.; Grahn, P.; Skärbäck, E.; Währborg, P. Green qualities in the neighbourhood and mental health—results from a longitudinal cohort study in Southern Sweden. BMC Public Health 2012, 12, 337. [Google Scholar] [CrossRef] [Green Version]
  60. Ekkel, E.D.; de Vries, S. Nearby green space and human health: Evaluating accessibility metrics. J. Landsc. Urban Plan. 2017, 157, 214–220. [Google Scholar] [CrossRef]
  61. Hansen, A.Y.; Meyer, M.R.U.; Lenardson, J.D.; Hartley, D. Built environments and active living in rural and remote areas: A review of the literature. Curr. Obes. Rep. 2015, 4, 484–493. [Google Scholar] [CrossRef] [PubMed]
  62. Schwartz, M.W.; Brigham, C.A.; Hoeksema, J.D.; Lyons, K.G.; Mills, M.H.; Van Mantgem, P. Linking biodiversity to ecosystem function: Implications for conservation ecology. Oecologia 2000, 122, 297–305. [Google Scholar] [CrossRef] [PubMed]
  63. Leemans, R.; de Groot, R. Millennium Ecosystem Assessment: Ecosystems and Human Well-Being: A Framework for Assessment; Island Press: Washington, DC, USA, 2003. [Google Scholar]
  64. Hough, R.L. Biodiversity and human health: Evidence for causality? Biodivers. Conserv. 2014, 23, 267–288. [Google Scholar] [CrossRef]
  65. Mills, J.G.; Brookes, J.D.; Gellie, N.J.; Liddicoat, C.; Lowe, A.J.; Sydnor, H.R.; Thomas, T.; Weinstein, P.; Weyrich, L.S.; Breed, M.F. Relating urban biodiversity to human health with the ‘holobiont’concept. J. Front. Microbiol. 2019, 10, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Janssen, I.; Boyce, W.F.; Simpson, K.; Pickett, W. Influence of individual- and area-level measures of socioeconomic status on obesity, unhealthy eating, and physical inactivity in Canadian adolescents. Am. J. Clin. Nutr. 2006, 83, 139–145. [Google Scholar] [CrossRef] [PubMed]
  67. Villanueva, R.; Albaladejo, R.; Astasio, P.; Ortega, P.; Santos, J.; Regidor, E. Socio-economic environment, area facilities and obesity and physical inactivity among children. Eur. J. Public Health 2015, 26, 267–271. [Google Scholar] [CrossRef]
  68. Gordon-Larsen, P.; Nelson, M.C.; Page, P.; Popkin, B. Inequality in the built environment underlies key health disparities in physical activity and obesity. Pediatrics 2006, 117, 417–424. [Google Scholar] [CrossRef] [Green Version]
  69. Noonan, R.J. Prevalence of Childhood Overweight and Obesity in Liverpool between 2006 and 2012: Evidence of Widening Socioeconomic Inequalities. Int. J. Environ. Res. Public Health 2018, 15, 2612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Wu, J.; He, Q.; Chen, Y.; Lin, J.; Wang, S. Dismantling the fence for social justice? Evidence based on the inequity of urban green space accessibility in the central urban area of Beijing. Environ. Plan B Urban Anal. City Sci. 2020, 47, 626–644. [Google Scholar] [CrossRef]
  71. Chung, A.; Backholer, K.; Wong, E.; Palermo, C.; Keating, C.; Peeters, A. Trends in child and adolescent obesity prevalence in economically advanced countries according to socioeconomic position: A systematic review. Obes. Rev. 2016, 17, 276–295. [Google Scholar] [CrossRef] [PubMed]
  72. Denny, S.; Lewycka, S.; Utter, J.; Fleming, T.; Peiris-John, R.; Sheridan, J.; Rossen, F.; Wynd, D.; Teevale, T.; Bullen, P.; et al. The association between socioeconomic deprivation and secondary school students’ health: Findings from a latent class analysis of a national adolescent health survey. Int. J. Equity Health 2016, 15, 109. [Google Scholar] [CrossRef] [Green Version]
  73. Sherry, B.; Jefferds, M.E.; Grummer-Strawn, L.M. Accuracy of adolescent self-report of height and weight in assessing overweight status: A literature review. Arch. Pediatr. Adolesc. Med. 2007, 161, 1154–1161. [Google Scholar] [CrossRef]
  74. Karchynskaya, V.; Kopcakova, J.; Klein, D.; Gába, A.; Madarasova-Geckova, A.; Van Dijk, J.P.; De Winter, A.F.; Reijneveld, S.A. Is bmi a valid indicator of overweight and obesity for adolescents? Int. J. Environ. Res. Public Health 2020, 17, 4815. [Google Scholar] [CrossRef]
  75. Gosse, M. How accurate is self-reported BMI? Nutr. Bull. 2014, 39, 105–114. [Google Scholar] [CrossRef]
  76. He, J.; Cai, Z.; Fan, X. Accuracy of using self-reported data to screen children and adolescents for overweight and obesity status: A diagnostic meta-analysis. Obes. Res. Clin. Pract. 2017, 11, 257–267. [Google Scholar] [CrossRef]
  77. Hackshaw, A. Small studies: Strengths and limitations. Eur. Respir. J. 2008, 32, 1141–1143. [Google Scholar] [CrossRef]
  78. Snijders, T.A.B. Power and sample size in multilevel linear models. Encycl. Stat. Behav. Sci. 2005. [Google Scholar] [CrossRef]
Table
Table 1. Sample characteristics.
Table 1. Sample characteristics.
Normal WeightOverweight and Obesity
Total Sample
(n = 62)		Total Sample
(n = 48)		Girls
(n = 25)		Boys
(n = 23)		Total Sample
(n = 14)		Girls
(n = 7)		Boys
(n = 7)
Age (years), M ± SD16.4 ± 1.516.3 ± 1.416.0 ± 1.516.6 ± 1.316.1 ± 1.316.2 ± 1.616.0 ± 1.1
Weight (kg), M ± SD61.4 ± 11.158.3 ± 8.854.9 ± 8.661.9 ± 7.874.0 ± 10.270.3 ± 8.377.6 ± 11.2
Height (m), M ± SD1.67 ± 0.101.67 ± 0.11.62 ± 0.91.72 ± 0.71.65 ± 0.11.6 ± 0.11.7 ± 0.1
BMI (kg/m2), M ± SD21.9 ± 3.220.7 ± 2.020.6 ± 2.420.8 ± 1.626.9 ± 1.827.1 ± 2.226.6 ± 1.4
MVPA Average Daily (minutes), M ± SD41.2 ± 20.841.5 ± 22.238.4 ± 18444.1 ± 25.742.2 ± 16.642.2 ± 13.942.2 ± 20.1
EDI-PT, M ± SD2.1 ± 0.82.19 ± 0.82.8 ± 0.82.3 ± 0.81.7 ± 0.82.1 ± 0.61.4 ± 0.7
BMI: body mass index; MVPA: moderate-to-vigorous physical activity; EDI-PT: Index of Deprivation for Small Portuguese Areas.
Table
Table 2. Descriptive statistics for each variable by gender and by body mass index categories.
Table 2. Descriptive statistics for each variable by gender and by body mass index categories.
Girls (n = 32)Boys (n = 30)
Normal Weight
(n = 25)		Overweight and Obesity
(n = 7)		p
(within Groups)		Normal Weight
(n = 23)		Overweight and Obesity
(n = 7)		p
(within Groups)
BMI (kg/m2), M ± SD20.7 ± 2.527.2 ± 2.20.7120.8 ± 1.726.6 ±1.40.62
MVPA Average Daily (minutes), M ± SD38.4 ± 18.542.2 ± 13.90.2544.1 ± 25.842.3 ± 20.10.63
EDI-PT (terciles), M ± SD2.1 ± 0.82.2 ± 0.70.372.3 ± 0.81.4 ± 0.80.48
UGS Counts [300 m], M ± SD0.4 ± 0.50.5 ± 0.70.270.2 ± 0.40.7 ± 0.70.04 *
UGS Counts [500 m], M ± SD1.0 ± 0.91.0 ± 1.10.530.8 ± 0.71.4 ± 1.00.53
UGS Counts [1000 m], M ± SD3.4 ± 2.13.4 ± 3.20.053.0 ± 2.44.8 ± 2.60.80
UGS Counts [1500 m], M ± SD6.8 ± 3.67.4 ± 3.90.645.8 ± 3.78.2 ± 3.30.70
NDVI [300 m], M ± SD0.17 ± 0.040.15 ± 0.030.080.18 ± 0.050.16 ± 0.030.16
NDVI [500 m], M ± SD0.18 ± 0.040.15 ± 0.030.710.18 ± 0.050.16 ± 0.020.36
NDVI [1000 m], M ± SD0.18 ± 0.030.17 ± 0.030.590.19 ± 0.030.17 ± 0.010.18
NDVI [1500 m], M ± SD018 ± 0.030.17 ± 0.030.410.19 ± 0.020.17 ± 0.020.22
SRI [300 m], M ± SD3.4 ± 3.92.4 ± 2.40.403.3 ± 3.33.1 ± 2.80.33
SRI [500 m], M ± SD3.7± 2.62.6 ± 2.20.473.6 ± 2.73.4± 2.50.85
SRI [1000 m], M ± SD3.7 ± 1.53.1 ± 1.40.594.1 ± 1.43.2 ± 0.90.33
SRI [1500 m], M ± SD3.7 ± 1.13.4 ± 1.30.394.0 ± 1.13.2 ± 0.50.06
Note: * significant difference (p < 0.05); BMI: body mass index; MVPA: moderate-to-vigorous physical activity; EDI-PT: Index of Deprivation for Small Portuguese Areas; UGS: urban green space; SRI: species richness index; NDVI: normalized difference vegetation index.
Table
Table 3. The results of the multiple linear regression model for body mass index considering the urban green space, normalized difference vegetation index and species richness index values for the sample without adjustment.
Table 3. The results of the multiple linear regression model for body mass index considering the urban green space, normalized difference vegetation index and species richness index values for the sample without adjustment.
ModelVariables 95% Confidence Interval for B
BSE BβLower BoundUpper Boundtp
1UGS Counts [300 m]0.5080.7580.091−1.0102.0770.6700.505
SRI [300 m]−0.2180.191−0.224−0.6000.165−1.1380.260
NDVI [300 m]6.74215.6080.086−24.50237.9860.4320.667
Constant21.4452.442 16.55726.338.782<0.001 **
Model 0.582
2UGS Counts [500 m]−0.3380.486−0.094−1.3100.635−0.6940.490
SRI [500 m]−0.0450.268−0.036−0.5830.492−0.1690.866
NDVI [500 m]−7.59818.662−0.088−44.95428.759−0.4070.685
Constant23.9762.833 18.30529.6478.463<0.001 **
Model 0.789
3UGS [1000 m]−0.0380.226−0.028−0.4910.416−0.1670.868
SRI [1000 m]−0.2980.563−0.124−1.4260.831−0.5280.599
NDVI [1000 m]10.09430.7920.086−51.56671.7540.3280.774
Constant21.4434.817 11.79831.0884.452<0.001 **
Model 0.956
4UGS [1500 m]−0.0230.203−0.025−0.4290.384−0.1120.911
SRI [1500 m]−0.3140.711−0.105−1.7271.119−0.4280.670
NDVI [1500 m]11.94037.4030.091−62.9319.8110.3190.751
Constant21.2166.725 7.75534.5773.1550.033 *
Model 0.974
Note: ** extremely significant difference (p < 0.001); * significant difference (p < 0.05); UGS: urban green space; SRI: species richness index; NDVI: normalized difference vegetation index.
Table
Table 4. The results of the multiple linear regression model for body mass index considering urban green space, normalized difference vegetation index and species richness index values with adjustment.
Table 4. The results of the multiple linear regression model for body mass index considering urban green space, normalized difference vegetation index and species richness index values with adjustment.
ModelVariables 95% Confidence Interval for B
BSE BβLower BoundUpper Boundtp
1UGS Counts [300 m]0.1450.7940.026−1.4451.7340.1820.856
SRI [300 m]−0.2140.190−0.221−0.5940.165−1.1320.263
NDVI [300 m]9.46715.5880.121−21.74940.6820.6070.546
EDI-PT−0.7950.557−0.199−1.9110.321−1.4270.159
Constant22.7882.597 17.58827.9898.774<0.001 **
Model 0.410
2UGS Counts [500 m]−0.5750.494−0.160−1.5630.414−1.1640.249
SRI [500 m]−0.1020.265−0.080−0.6330.428−0.3850.701
NDVI [500 m]−2.54118.502−0.029−39.59234.509−0.1370.891
EDI-PT 0.9880.540−0.248−2.0680.093−1.8300.072
Constant25.5822.913 19.74931.4158.783<0.001 **
Model 0.360
3UGS Counts [1000 m]−0.1290.226−0.096−0.5810.324−0.5690.572
SRI [1000 m]−0.6150.574−0256−1.7660.535−1.0710.289
NDVI [1000 m]20.23730.5310.172−40.92481.3980.6630.510
EDI-PT −1.0530.544−0.267−2.1430.038−1.9340.058
Constant23.2684.799 13.65532.8814.849<0.001 **
Model 0.405
4UGS Counts [1500 m]−0.0640.201−0.071−0.4660.339−0.3170.753
SRI [1500 m]−0.4860.707−0.167−1.9030.930−0.6870.495
NDVI [1500 m]14.70536.8340.113−59.05388.4630.3990.691
EDI-PT−0.8940.593−0.224−1.9400.153−1.7100.093
Constant23.5276.753 10.00537.0493.484<0.001 **
Model 0.538
Note: ** extremely significant difference (p < 0.001); * significant difference (p < 0.05); UGS: urban green space; SRI: species richness index; NDVI: normalized difference vegetation index; MVPA: moderate-to-vigorous physical activity; EDI-PT: European Index of Deprivation for Small Portuguese Areas.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Melo, J.; Ribeiro, A.I.; Aznar, S.; Pizarro, A.; Santos, M.P. Urban Green Spaces, Greenness Exposure and Species Richness in Residential Environments and Relations with Physical Activity and BMI in Portuguese Adolescents. Int. J. Environ. Res. Public Health 2021, 18, 6588. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126588

AMA Style

Melo J, Ribeiro AI, Aznar S, Pizarro A, Santos MP. Urban Green Spaces, Greenness Exposure and Species Richness in Residential Environments and Relations with Physical Activity and BMI in Portuguese Adolescents. International Journal of Environmental Research and Public Health. 2021; 18(12):6588. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126588

Chicago/Turabian Style

Melo, Juliana, Ana Isabel Ribeiro, Susana Aznar, Andreia Pizarro, and Maria Paula Santos. 2021. "Urban Green Spaces, Greenness Exposure and Species Richness in Residential Environments and Relations with Physical Activity and BMI in Portuguese Adolescents" International Journal of Environmental Research and Public Health 18, no. 12: 6588. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126588

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop