1095384 research-article2022 EABXXX10.1177/00139165221095384Environment and BehaviorEzpeleta et al. Article Environment and Behavior Greenspace Exposure and 2022, Vol. 54(5) 893­–916 © The Author(s) 2022 Article reuse guidelines: Obsessive-Compulsive sagepub.com/journals-permissions DOI: 10.1177/00139165221095384 https://doi.org/10.1177/00139165221095384 Behaviors in journals.sagepub.com/home/eab Schoolchildren Lourdes Ezpeleta1 , José Blas Navarro2, Lucía Alonso3, Núria de la Osa1, Albert Ambrós3, Mònica Ubalde3,4, Eva Penelo2, and Payam Dadvand3,4 Abstract Green environments are associated with improved child brain development and mental health. We study cross-sectionally the association of the availability of greenspace at home and school with obsessive-compulsive behaviors (OCB) in primary schoolchildren. Greenspace and tree cover surrounding home and school of 378 children aged 9 to 10 in Barcelona (Spain) were characterized using satellite-based indices [Normalized Difference Vegetation Index (NDVI), Modified Soil Adjusted Vegetation Index (MSAVI), Vegetation Continuous Field (VCF)] across buffers of 100, 300, and 500 m, and distance to the nearest green space. OCB was assessed with the Spence Children’s Anxiety Scale-Parent version. Linear and mixed effects models showed that greenspace at school, but not at home, was significantly related to a reduction in OCB across buffers, with benefits for girls and also children with graduate parents. Higher greenspace around the school might be associated with less obsessive-compulsive behavior in primary schoolchildren, especially in girls and those with higher socioeconomic status. 1 Departament de Psicologia Clínica i de la Salut, Universitat Autònoma de Barcelona, Spain 2 Departament de Psicobiologia i Metodologia de les Ciències de la Salut, Universitat Autònoma de Barcelona, Spain 3 Barcelona Institute for Global Health (ISGlobal), Spain 4 Universitat Pompeu Fabra (UPF), Barcelona. CIBER Epidemiología y Salud Pública (CIBERESP), Spain Corresponding Author: Lourdes Ezpeleta, Departament de Psicologia Clínica i de la Salut. Edifici B, Universitat Autònoma de Barcelona, Bellaterra (Barcelona) 08193, Spain. Email: lourdes.ezpeleta@uab.cat 894 Environment and Behavior 54(5) Keywords greenspace, pediatric psychology, natural environment, behavioral development The ongoing urbanization worldwide has substantially reduced the contact of the individuals with nature. A growing body of evidence has associated con- tact with nature, including greenspace, with health benefits such as improv- ing mental health (J. Zhang, Yu, et al., 2020). In adults, exposure to greenspace has been associated with significant benefits for psychological wellbeing and mental health (Bratman et al., 2019; Dadvand et al., 2016). Similarly, in chil- dren, a number of studies have reported that exposure to residential and school greenspace is associated with improved neurodevelopmental out- comes such as higher IQ levels (Bijnens et al., 2020; Lee et al., 2021), improvement in academic performance (Donovan et al., 2020), more proso- cial behavior (Amoly et al., 2014; Putra et al., 2021), improved early psycho- motor development (Kabisch et al., 2019; Liao et al., 2019), enhanced cognitive development (Dadvand et al., 2015, 2017), improved self-satisfac- tion and well-being (Dadvand et al., 2019; Feng & Astell-Burt, 2017; McCormick, 2017; McEachan et al., 2018), and is also associated with ben- eficial anatomical changes in developing brain (Dadvand et al., 2018). Specifically, exposure to greenspace has been related with lower risk of behavioral and emotional problems such as externalizing behavior problems (Bijnens et al., 2020; Lee et al., 2019), anxiety and depression (Larson et al., 2018; Liao et al., 2020), attention deficit/hyperactivity disorder (Amoly et al., 2014; Lee et al., 2019; Liao et al., 2020; Markevych et al., 2014; Vanaken & Danckaerts, 2018; Yang et al., 2019), schizophrenia (Engemann et al., 2018, 2020), and poorer general mental health (Andrusaityte et al., 2020). Obsessive-compulsive disorder (OCD) refers to the presence of unwanted recurrent thoughts, urges, or images (obsessions) and/or unnecessary repeti- tive behaviors or mental acts performed in response to obsession (compul- sions) (American Psychiatric Association, 2013). OCD thoughts and rituals are identified as nonsensical. In about half of the cases the disorder starts in childhood, with prevalence until adolescence ranging from 0.1% to 4% (Heyman et al., 2001; Rapoport et al., 2000). The severe form of this disorder could be associated with serious distress and impairment in the daily life of not only the affected child, but also the family, which is usually involved in the rituals (Monzani et al., 2020). According to a recent systematic review and meta-analysis of the available evidence, persistence rates for full pediat- ric OCD diagnosis are about 40% at least 1 year later, with shorter duration being a predictor of remission and hence a better prognosis (Liu, Cui et al., Ezpeleta et al. 895 2021). These findings highlight the need for identifying the measures that may help to early hinder the progression and continuity of symptoms to full OCD. The proposed mechanisms through which greenspace could be beneficial on health include: (a) reducing harm such as mitigating exposure to noise, heat or air pollution, (b) restoring capacities, mainly through reducing physi- ological stress and recovering attention, and (c) building capacities such as improving physical fitness or promoting social contact and cohesion (Markevych et al., 2017). Although through some of these mechanisms, greenspace could also have an impact on OCD, the available evidence on such an impact is still very scarce. Therefore, the aim of our study was to evaluate the association of exposure to greenspace at home and school with obsessive-compulsive behaviors (OCB) in primary schoolchildren. Materials and Methods Study Setting The study was conducted in the Barcelona metropolitan area (636 km2) located in the northeast coast of Spain with a population of 5.6 million inhab- itants, which is the second most populated province in Spain (Instituto Nacional de Estadística, 2020). The area has a Mediterranean climate, with hot and dry summers and mild winters (average annual temperature is 21.2°C during the day and 15.1°C at night) and much of the precipitation occurring during spring and autumn months (on average 78 rainy days a year). Study Population This study was part of a well-established ongoing cohort that aims to identify the risk factors of psychological problems in children (Ezpeleta et al., 2014). Following a two-phase sample design an initial sample of 2,283 children from 54 schools (25.9% semi-public and 74.1% public) were randomly selected from the census of kindergarten schools across our study area in 2009 to 2010 (Figure 1). A total of 1,341 families (58.7%) agreed to partici- pate. In the second stage, children were screened for psychological problems using the parent-administered Strengths and Difficulties Questionnaire (SDQ) (Goodman, 1997). All the children who screened positively for behav- ioral problems and an additional 30% of the children with scores below the cut-off screening scores were enrolled in the longitudinal study. Of those included, 622 children (89.4 %) and their parents agreed to participate in the longitudinal study, with no statistically significant differences in sex (p = .82) 896 Environment and Behavior 54(5) Figure 1. Study sample and follow-ups. or type of school (p = .85) found between participants and drop-outs; how- ever, children who did agree to participate were of higher socioeconomic status (SES) (86.2% of high SES for participants vs. 73.6% of high SES for drop-outs; p = .01). Children were followed annually (i.e., six follow-ups until the assessment of OCB) for the evaluation of their neurodevelopment and mental health status and the corresponding risk factors. Further details regarding this cohort has been published previously (Ezpeleta et al., 2014). The final sample for the present study included 378 children at age 9 years (i.e., sixth follow-up – cross-sectional analysis) when the evaluation of the OCB was conducted (Figure 1), 244 children (64.6%) from the screening positive initial group. There was no bias in sex (p = .62) neither in type of Ezpeleta et al. 897 school (p = .18) due to attrition. Nevertheless, children who remained in the study at age 9 had a higher SES (p < .01) than those at the first follow-up. The project was approved by the Ethics Committee on Animal and Human Experimentation of the Universitat Autònoma de Barcelona (Ethics approval CEEAH 4324). Families were recruited at the schools and gave written con- sent for the assessment. The families who agreed were contacted by tele- phone and interviewed at the school for each annual assessment. Exposure Assessment Assessment of greenspace exposure included two different aspects: namely surrounding greenspace and access to greenspace, and was conducted, sepa- rately, for home as well as school. Surrounding greenspace. We applied the following three satellite-based indi- ces to identify greenspace, separately for all vegetation types together as well as trees, surrounding home and school for each participant. (1) The Normalized Difference Vegetation Index (NDVI) which is an indicator of photosynthetically active landcover, calculated based on the land surface reflectance of red and near infrared parts of the spec- trum. Its values range from −1 to 1 with higher values indicating the higher greenspace (NCLB, 2000). To englove the study population there were used two images from 2015 (Landsat 8 OLI/TIRS) from the greenest season, April and May (Landsat images for 2015 U.S. Geological Survey, 2015). (2) The Vegetation Continuous Field (VCF) is an indicator of tree canopy cover, and displays the percentage of land covered by the woody veg- etation with a height greater than 5 m (Sexton et al., 2013). The higher VCF values indicate the higher percentage of tree cover. VCF is derived from Landsat images and at periods of 5 years from 2000. Here were used 2015, and also were used two images to spatially cover the study population (Land Cover Land Use Change (NASA) (2015). The address of residence for each participant were obtained at baseline and during each follow-up. We defined home surrounding greenspace as the average of NDVI, VCF, and MSAVI across buffers of 100, 300, and 500 m in a circular buffer area (representing immedi- ate, intermediate, and neighborhood areas, respectively) around all the home addresses where children had lived between age 3 years (baseline visit) and 9 years (assessment of OCD), weighted by the time that the child had resided in each address. Similarly, we defined 898 Environment and Behavior 54(5) school surrounding greenspace as average of NDVI, VCF, and MSAVI across the aforementioned buffer sizes around all the schools where children had attended between age 3 and 9 years, weighted by the time that the child had spent at each school. Access to green spaces. The distance to the closest green space (DCG) was applied as an indicator of access to green spaces We applied the Urban Atlas of the Copernicus Land Monitoring Service for 2012 (European Union, 2020) to identify green spaces in our study area. We then calculated the Euclidean distance separately from home and school to their nearest green space. As such higher values corresponded to less exposure to greenspace. Assessment of Obsessive-Compulsive Behaviors Pediatric OCD has a typical onset in later childhood or adolescence (Geller, 2006). Accordingly we chose age 9 to 10 to characterize this condition in our cohort, using the Spence Children’s Anxiety Scale-Parent version (SCAS-P) (Spence, 1998) questionnaire. SCAS-P assesses the severity of anxiety symp- toms following the anxiety disorder categories proposed by the DSM-IV (American Psychiatric Association, 2000). It contains 38 items (0: never to 3: always) organized in six domains of anxiety one of which being OCD. The obsessive-compulsive scale includes six items about checking, repeating behaviors, ritualistic behavior, intrusive recurrent thoughts, special thoughts, and bad thoughts. Whiteside et al. (2012) studied the convergent validity through examining correlations of the subscale with ratings from standardized clinical diagnostic interviews, such as the Children’s Yale-Brown Obsessive Compulsive Scale (CY-BOCS) (Scahill et al., 1997) and Anxiety Disorders Interview Schedule for DSM-IV: Child Version (ADIS-C) (Silverman et al., 1994). SCAS obses- sive-compulsive scale parent-reports were positively correlated with the total score and subscale scores on the CY-BOCS and the number of obsessions and compulsions and distress associated with OCD symptoms as reported on the ADIS:C. SCAS obsessive-compulsive scale scores also showed criterion validity given that mean scores differed significantly between patients with an OCD diagnosis, patients without an OCD diagnosis, and children from the community, and the mean from the OCD group was significantly greater, with large effect sizes, than the means of the other two groups. In the clinical sample pre-treatment scores of SCAS obsessive-compulsive scale were sig- nificantly greater than post-treatment scores for parent report with large effect sizes (d = 1.28). Ezpeleta et al. 899 The raw data was transformed to T-scores according to the test manual (Spence, n.d.), which were used as the indicator of OCB. Ordinal alpha for SCAS-P obsessive-compulsive scale was .91 in our study sample. Covariates Data on sex, ethnicity, and parental education were obtained through inter- view as a part of the routine demographic data collection in our cohort. Neighborhood sociodemographic status characteristics such as the percentage of unemployed people, and the percentage of people with no or primary edu- cation at the census tract where home and school were located were also abstracted. The Life Events Checklist (Johnson & McCutcheon, 1980) was also applied to obtain data on child’s stressful life events during the year previ- ous to the assessment (moving home, change school, parent’s divorce, etc.) from the mother or caregiver. Statistical Analysis Main analyses. Given the multistage sampling process, the analyses were weighted by the inverse probability of the selection in the second phase of sampling (Lohr, 2019), as we did in our previous studies based on the data from this cohort (Ezpeleta et al., 2014). For evaluating the association between the indicators of exposure to greenspace at home and OCB, we developed linear mixed models including OCB as the outcome, each measure of greenspace exposure at each buffer (one at a time) as predictor, and school as the random effect. In the analysis of greenspace at school all the children going to the same school had the same greenspace value, with no within- school variability. Therefore, the random effect of school was not included in analysis of the association of school greenspace with OCB. Accordingly, for estimating the association of greenspace exposure at school with OCB, we applied linear regression models. All models were adjusted for a priori selected set of covariates including (i) sex, (ii) parental education (with vs. without university education), (iii) ethnicity (White vs. other), (iv) two indi- cators of neighborhood sociodemographic status (both for the home address and the school’s address), and (v) number of stressful life events the child has lived during the year preceding the assessment of the OCD. Exposures were scaled by their interquartile range (IQR) in order to provide comparable results for different exposure indicators. Because the low percentage of miss- ing, data pairwise analysis was done using for each model all the available data. All the statistical analyses were carried out using Stata 16 (StataCorp, 2019). The type I error was set at 0.05 value. 900 Environment and Behavior 54(5) Table 1. Description of the Sociodemographic Characteristics and Obsessive- Compulsive Symptoms in Study Participants (N = 378). Characteristics Description Age (years)a 9.7 (0.5) Sex; % Female 49.2 Parent’s education: university studies; % Yes 68.8 Ethnicity; % White 93.1 Number of stressful life events 1.0 (2.0) % of unemployment at neighborhood 18.7 (10.6) % of Primary or no education at neighborhood 7.2 (7.3) SCAS-P Obsessive-Compulsive 47.0 (2.0) Note. SCAS-P OC: Spence Children’s Anxiety Scale-Parent version. a Continuous variables described by their median (IQR). Sensitivity and stratified analyses. As a first sensitivity analysis, we used Modi- fied Soil-Adjusted Vegetation Index (MSAVI) to assess home and school sur- rounding greenspace and repeated the analyses for this alternative set of exposure. The MSAVI improves NDVI measures by minimizing the effect of the soil background and as such could provide a more precise estimation of the vegetation where there are areas with few chlorophyll or high proportions of bare soil (Qi et al., 1994). Similar to NDVI, MSAVI values range from −1 to 1 with higher values indicating more greenspace (NCLB, 2000). Addition- ally, as a second sensitivity analysis and to explore influence of ethnicity on our main results, we repeated the analyses by limiting the participants to White children. Given the small percentage of non-white children (6.9%), it was not possible to conduct a separate analysis for this sub-group. Finally, as a third sensitivity analysis we re-estimated models by further adjusting for the best job occupation from both parents (unemployed, manual/unskilled worker, non-manual/skilled worker). To explore the potential variation of our evaluated associations across strata of sex, SES, and presence of behavior problems, we stratified our anal- yses at home and school surrounding greenspace (at 300 m buffer) and DCG, by sex, by parental education (with or without university education), and by having or not behavior problems at recruitment at age 3. Results The sociodemographic characteristics of study participants are described in Table 1. Of 378 participants included in our study, 49.2% were girls, 68.8% had at least one parent with university studies, 93.1% where white and had a Ezpeleta et al. 901 Table 2. Spearman’s Correlation Coefficients Between Indicators of Exposure to Greenspace at Home (Above Diagonal) and School (Below Diagonal). 1 2 3 4 5 6 7 1. NDVI 100 m. — .86 .80 .85 .81 .78 −.54 2. NDVI 300 m. .90 — .96 .70 .91 .92 −.65 3. NDVI 500 m. .88 .98 — .61 .84 .94 −.62 4. VCF 100 m. .91 .81 .79 — .77 .65 −.46 5. VCF 300 m. .90 .97 .93 .81 — .91 −.61 6. VCF 500 m. .87 .96 .97 .79 .96 — −.60 7. DCG m. −.69 −.76 −.75 −.60 −.76 −.78 — Note. NDVI = normalized difference vegetation index; VCF = vegetation continuous field (% of tree covering); DCG = distance to closest green space; italic: above diagonal correlation values at home. median of 1 stressful life event during the year preceding the OCD evaluation. Table 2 shows the Spearman’s correlation matrix between the different greenspace measures and buffer sizes both at residence and at school. The pattern of relationships was similar in both locations, with slightly higher values at home than at school. There were strong correlations between NDVI and VCF measures when compared with each other and moderate to strong correlations when they were compared with the DCG. As expected, the stron- gest correlations between NDVI and VCF were found for the same buffer size. Table 3 presents the description of greenspace measures and its associa- tion with OCB, separately for home and school greenspace measures. Significant associations were found for school NDVI and VCF, with negative coefficients reflecting that a greater exposure to greenspace was associated with lower OCB (i.e., less OCD symptoms). The magnitude of the associa- tion was greater for smaller buffers (100 m.) than for the largest (300 and 500 m). An increment of one IQR in NDVI measured within 100 and 500 m buffer was associated with a significant decrease of between −1.24 (95% confidence intervals (CIs: −2.36, −0.13) and −1.70 (95% CIs: −3.35, −0.04) points in OCB scale. Similarly, an increment of one IQR in VCF was associ- ated with a significant decrease of between −0.57 (95% CIs: −0.92, −0.22) and −1.22 (95% CIs: −2.25, −0.19) points in OCB. For the school DCG, and also for the home surrounding greenspace and DCG, although we observed beneficial associations with OCB, none of these associations attained statisti- cal significance. 902 Environment and Behavior 54(5) Table 3. Adjusteda Regression Coefficients for One Interquartile (IQR) Increase in Each Indicator of Greenspace Exposure in Association with Obsessive- Compulsive Score. Homeb Schoolc Greenspace Association Greenspace Association a Median (IQR) B (CI 95%) Median (IQR) Ba (CI 95%) NDVI 100 m buffer 0.176 (0.095) −0.51 (−1.48, 0.46) 0.180 (0.156) −1.70 (−3.35, −0.04)** 300 m buffer 0.194 (0.116) −0.59 (−1.95, 0.77) 0.238 (0.167) −1.24 (−2.36, −0.13)** 500 m buffer 0.197 (0.128) −0.49 (−1.99, 1.02) 0.236 (0.170) −1.27 (−2.51, −0.03)** VCF 100 m buffer 2.498 (2.062) −0.14 (−0.48, 0.20) 2.811 (3.002) −1.22 (−2.25, −0.19)** 300 m buffer 2.904 (2.095) −0.11 (−0.38, 0.17) 3.635 (3.148) −0.57 (−0.92, −0.22)** 500 m buffer 3.043 (2.329) −0.13 (−0.43, 0.17) 3.571 (4.278) −0.85 (−1.41, −0.30)** DCG (m) 185.3 (244.3) 1.09 (−0.24, 2.42) 206.6 (288.0) 0.73 (−0.76,2.22) Note. NDVI = Normalized Difference Vegetation Index; VCF = Vegetation Continuous Field (% of tree covering); DCG = Distance to closest green space. a Models were adjusted for sex, ethnicity, parent’s education level (university vs. other), percentage of unemployment, percentage of illiterate and total number of stressful events. CI: Confidence interval. IQR: Interquartile range. b Based on mixed effects models with school as the random effect. c Based on linear regression models. **p < .05. Sensitivity and Stratified Analyses Supplemental Table S1 presents the results of the analyses using MSAVI to assess home and school surrounding greenspace. The associations for both locations were non-significant but showed the same direction than those observed for NDVI and VCF in main analyses, with coefficients of greater magnitude for school than for home. Otherwise, after limiting analyses to only white children, the observed associations were in line with those of main analyses with regards to the direction and statistical significance (Supplemental Table S2). Also, when adding parents’ occupation as adjusting term, regres- sion coefficients were similar to those without, although statistical signifi- cance for NDVI at 100 and 500 m was lost (Supplemental Table S3). Figure 2 presents the association between greenspace exposure at home/ school and OCB stratified by sex. For home surrounding greenspace confi- dence intervals included the zero-value reflecting non-significant associa- tions. At school, a protective association of greenspace exposure for girls when measured with NDVI or VCF was obtained. Ezpeleta et al. 903 Figure 2. Adjusteda Associations (95% Confidence Intervals) Between One Interquartile Increase in Greenspace Exposure at Home/School and Obsessive- Compulsive Score Stratified by Sex. a Models were adjusted for ethnicity, parent’s education level (university vs. other), percentage of unemployment, percentage of illiterate, and total number of stressful events. Figure 3 shows the association between greenspace exposure at home/ school and OCB stratified by parental education (at least one parent with university education versus both parents without university education). For NDVI and VCF, we did not find a notable difference between the two strata in terms of the direction of the associations: in both strata we observed a protective association between these greenspace exposure indicators and OCB, which was statistically significant only for those with at least one par- ent having university education when greenspace exposure was obtained at school. For the DCG, the association varied between locations for those chil- dren with parents without any university education, showing statistically sig- nificant protective association against OCB at home but null at school. Figure 4 presents the association between greenspace exposure at home/ school and OCB stratified by having or not behavior problems at age 3. 904 Environment and Behavior 54(5) Figure 3. Adjusteda Associations (95% Confidence Intervals) Between One Interquartile Increase in Greenspace Exposure at Home/School and Obsessive- Compulsive Score Stratified by Parental Education (Without University Studies/ With University Studies). a Models were adjusted for sex, ethnicity, percentage of unemployment, percentage of illiterate, and total number of stressful events. Results exhibited a very similar pattern for both groups with two exceptions: DCG measured at home, showing a protective significant association with OCB only for children without behavior problems, and NDVI measured at school, which associated with a reduction of OCB in the group with behavior problems. Discussion To our knowledge this is the first study that evaluated the association between greenspace at both home and school and OCB in community children. Our Ezpeleta et al. 905 Figure 4. Adjusteda Associations (95% Confidence Intervals) Between One Interquartile Increase in Greenspace Exposure at Home/School and Obsessive- Compulsive Score Stratified by Behavior Problems at age 3 (Without/With). a Models were adjusted for sex, ethnicity, parent’s education level (university vs. other), percentage of unemployment, percentage of illiterate, and total number of stressful events. study benefited from data from a well-established cohort of children with prospectively collected data on important covariates (e.g., sociodemographic factors and stressful life events) and characterization of OCB using a widely- used tool together with a comprehensive characterization of greenspace exposure at both home and school. We found that greenspace surrounding the school had a potentially beneficial association with obsessive-compulsive behavior, reflected as lower OCB scores. For the greenspace surrounding home, and also home and school DCG, we also observed suggestions for a possible protective association; however, they were not statistically signifi- cant. Furthermore, we observed some suggestions for a potentially stronger association for girls, and perhaps those children with higher SES. 906 Environment and Behavior 54(5) We are aware of only one previous study reporting on the association of greenspace exposure and OCD. This study conducted in Denmark reported higher greenspace surrounding home (measured as NDVI average) was asso- ciated with lower risk of OCD (Engemann et al., 2019). This finding was in line with our observed association between home surrounding greenspace and OCB; however, our observed associations were not statistically signifi- cant. In contrast with this previous study, our study also characterized greens- pace at school, for which we observed stronger associations with OCB compared to those for greenspace at home. Our observed stronger association for greenspace at school is in line with other studies that also found that school but not home surrounding greenspace indicators were significantly associated with cognitive function such as working memory and inattentive- ness (Dadvand et al., 2015), attention deficit/hyperactivity disorder (Yang et al., 2019), and biological markers such as blood pressure or cholesterol (Ribeiro et al., 2019). Children spend large amounts of time at the school, and it is at the schoolyards where they carry out most of the outside activities such as sports, socializing with other peers, and other recreational leisure activities during the week. Greenspace at school may be restorative (Kelz et al., 2015; Y. Zhang, Mavoa, et al., 2020) and can promote reduction of stress, relax- ation, and distraction from repetitive thoughts and behaviors. The associa- tions found at school are especially important for public health because: (a) focused interventions at schools are more feasible and affect more children compared to interventions at neighborhoods and across cities; (b) children are individuals in development with great plasticity to environmental influ- ences such as greenspace at school and their associated benefits; (c) increas- ing greenspace at school in some way is feasible in most cases; and (d) increasing greenspace at school could contribute to improve children’s men- tal health. We observed stronger associations with OCB for the indicator of general greenspace (NDVI) compared to the indicator of the tree cover (VCF) at school. These results are in line with findings of an Australian study in adults, which also showed a slightly stronger association of their indicator of general greenspace with psychological distress, depression or anxiety, and general health compared with their indicator of tree cover (Astell-Burt & Feng, 2019). Similarly, another study showed that higher home surrounding greenspace, measured as NDVI average, was associated with lower cortisol levels (an indi- cator of stress) in the cord blood, while for tree cover surrounding the home address, measured using VCF, no association was found (Boll et al., 2020). The same mechanisms by which greenspace could affect mental health and wellbeing, which are yet to be fully established, could also help to explain the associations found with OCD. Genetic and environmental factors could Ezpeleta et al. 907 influence biological and psychological characteristics that ultimately lead to OCD (Rapoport & Shaw, 2015). Heritability rates of OCD in childhood range between 45% and 65% (van Grootheest et al., 2005). Nonetheless, environ- mental influences, such as perinatal injury, childhood infection, or life stress- ful events have also been suggested to be relevant (Grisham et al., 2008). These environmental influences may result in diverse types of stress that, in turn, may alter different biological systems, such as the immune system (Bottaccioli et al., 2019). OCD is a neuropsychiatric disorder that involves the cortico-striato-thalamo-cortical (CSTC) loop, implicated in action regula- tion, whose abnormal hyperactivation could lead to repetitive thoughts and actions (obsessions and compulsions) (Dougherty et al., 2018). A controlled experiment in adults has shown that being exposed to nature could reduce rumination and the activity in the brain area related with ruminations (sub- genual prefrontal cortex) (Bratman et al., 2015), which is connected to the CSTC brain loop. Rumination is a cognitive process which consists of repeti- tive and unproductive thoughts that are frequent in OCD. Therefore, contact with nature, through distracting from the ruminations, could reduce OCD- related behaviors and levels. Since early childhood the exposure to greenspace seems to favor neurode- velopment, as shown by cognitive and psychomotor execution (Dadvand et al., 2015; Liao et al., 2019). Children with OCD show greater executive dysfunction in flexibility (shift), working memory, planning, and inhibition (Garcia-Delgar et al., 2018; Negreiros et al., 2020; Ornstein et al., 2010). As such, the ability of greenspace exposure to enhance the development of cog- nitive and executive functions could reduce the risk of OCD. In line with these findings, our previous study in a sample of schoolchildren in Barcelona (Dadvand et al., 2018) identified enhancements in some brain regions (pre- frontal cortex, premotor cortex, and cerebellar hemispheres) in urban chil- dren raised in greener neighborhoods. These brain regions were also associated with OCD (Apergis-Schoute et al., 2018) and with better scores in working memory and attention, suggesting the potential contribution of exposure to greenspace to brain development (Dadvand et al., 2018). Finally, another alternative path from lack of greenspaces to OCB might be through the stress produced by the environmental influences (i.e., air pollution, noise, and heat), which may affect the immune system (Bottaccioli et al., 2019) and in turn the basal ganglia involved in the control of voluntary movements (Boileau, 2011). Along this line, proximity of the school to greenspaces is an environmental factor that has been associated with a lower allostatic load in children (Ribeiro et al., 2019). The allostatic load refers to how chronic envi- ronmental stress favors the dysregulation of multiple biological systems (hor- mones, metabolism, immune system, cardiovascular, etc.) of the individuals, 908 Environment and Behavior 54(5) making them more vulnerable to disease (McEwen & Stellar, 1993). The mechanism by which allostatic load improves might be through the enhance- ment of cardiovascular, metabolic and other biologic systems because activ- ity increases when greenspace is available, the ameliorative effects of greenspace on quality of air, or the relaxation effects associated with expo- sure to greenspace (Ribeiro et al., 2019). We found suggestions for a potentially stronger association of surrounding greenspace with OCB in girls; however, for the DCG, we did not observe any notable difference between girls and boys. In which sex the influence of greenspace is greater is a matter of debate and, until now the available evi- dence is heterogeneous. Studies in adults have suggested that psychological and physiological characteristics might influence how each sex is affected by greenspace. For instance, whereas in men the objective amount of greenspace could better associate with health, in women the associations could be stron- ger with more subjective indicators such as quality of greenspaces or their restorative value (Bolte et al., 2019). However, our findings suggest that greenspace might be contributing to the slight preponderance of OCD in pre- pubertal boys versus girls (Rapoport & Shaw, 2015). With regards to SES, while for surrounding greenspace we did not observe a considerable differ- ence across the strata of parental education, for DCG, we observed some suggestions for a potentially stronger association for those children from lower SES. Although there is not full agreement about the direction of SES in the association with greenspace, a recent systematic review has reported that the benefits for physical health of greenspace might be more in low SES than in high SES groups (Rigolon et al., 2021). The results of our study should be interpreted in light of some limitations. First, our study, given its cross-sectional design had a limited capacity to establish a causal relationship. Moreover, given our modest sample size, our stratified analyses, particularly based on parental education, could have been under-powered. Also, because children in the age when we assessed their OCB (i.e., 9 years) could have questionable reliability and validity in answer- ing questionnaires (Schniering et al., 2000), only parents reported on OCB, which could have resulted in the misclassification of our outcome variable. Additionally, our results should be generalized with caution to families of low SES given that these families were under-represented in the analyzed sample. Furthermore, our assessment of exposure to greenspace did not include quality aspects such as safety, esthetic, biodiversity, facilities, ameni- ties, and organized activities, which could have influenced our evaluated associations. The ongoing urbanization worldwide is leading an increasing number of children to live in urban areas, where residents might be at higher risk of increasing obsessive-compulsive behaviors (Rigolon et al., 2021) and often Ezpeleta et al. 909 have less access to greenspace. Obsessive-compulsive behaviors often have an early onset and a large proportion of cases progress to full OCD. Therefore, childhood could be a window of opportunity for primary and secondary pre- vention of this condition. We observed that higher greenspace at school might reduce OCB in primary schoolchildren. Our findings, if confirmed by future studies, could provide an evidence base for policymakers to implement inter- ventions (e.g., greening the schools and its surroundings) to reduce the inci- dence of this condition and its accompanying burden on the family and society in our rapidly urbanizing world. Further longitudinal studies based on large sample sizes with repeated characterization of OCB are required to con- firm our findings in other settings with different climates and cultures. Declaration of Conflicting Interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Spanish Ministry of Science, Innovation and Universities [Grant PGC2018-095239-B-I00 (MICIU/FEDER)] Ethical Approval The project was approved by the Ethics Committee on Animal and Human Experimentation of the Universitat Autònoma de Barcelona. Families gave written consent for the assessment. ORCID iD Lourdes Ezpeleta https://orcid.org/0000-0002-8957-083X Supplemental Material Supplemental material for this article is available online. References American Psychiatric Association. (2000). DSM-IV Diagnostic and statistical manual of mental disorders (4th text rev. ed.). American Psychiatric Press. American Psychiatric Association. (2013). Diagnostic and statistical manual of men- tal disorders (5th ed.). American Psychiatric Association. Amoly, E., Dadvand, P., Forns, J., López-Vicente, M., Basagaña, X., Julvez, J., Alvarez-Pedrerol, M., Nieuwenhuijsen, M. J., & Sunyer, J. (2014). Green and blue spaces and behavioral development in Barcelona schoolchildren: The 910 Environment and Behavior 54(5) BREATHE project. Environmental Health Perspectives, 122(12), 1351–1358. https://doi.org/10.1289/ehp.1408215 Andrusaityte, S., Grazuleviciene, R., Dedele, A., & Balseviciene, B. (2020). The effect of residential greenness and city park visiting habits on preschool children’s mental and general health in Lithuania: A cross-sectional study. International Journal of Hygiene and Environmental Health, 223(1), 142–150. https://doi. org/10.1016/j.ijheh.2019.09.009 Apergis-Schoute, A. M., Bijleveld, B., Gillan, C. M., Fineberg, N. A., Sahakian, B. J., & Robbins, T. W. (2018). Hyperconnectivity of the ventromedial prefrontal cortex in obsessive-compulsive disorder. Brain and Neuroscience Advances, 2, 1–10. https://doi.org/10.1177/2398212818808710 Astell-Burt, T., & Feng, X. (2019). Association of urban green space with mental health and general health among adults in Australia. Jama Network Open, 2(7), e198209. https://doi.org/10.1001/jamanetworkopen.2019.8209. Bijnens, E. M., Derom, C., Thiery, E., Weyers, S., & Nawrot, T. S. (2020). Residential green space and child intelligence and behavior across urban, suburban, and rural areas in Belgium: A longitudinal birth cohort study of twins. PLoS Medicine, 17(8), e1003213. https://doi.org/10.1371/journal.pmed.1003213 Boileau, B. (2011). A review of obsessive-compulsive disorder in children and ado- lescents. Dialogues in Clinical Neuroscience, 13(4), 401–411. Boll, L. M., Khamirchi, R., Alonso, L., Llurba, E., Pozo, Ó. J., Miri, M., & Dadvand, P. (2020). Prenatal greenspace exposure and cord blood cortisol levels: A cross- sectional study in a middle-income country. Environment International, 144. Article 106047. https://doi.org/10.1016/j.envint.2020.106047 Bolte, G., Nanninga, S., & Dandolo, L. (2019). Sex/gender differences in the associa- tion between residential green space and self-rated health-a sex/gender-focused systematic review. International Journal of Environmental Research and Public Health, 16(23), 4818. https://doi.org/10.3390/ijerph16234818 Bottaccioli, A. G., Bottaccioli, F., & Minelli, A. (2019). Stress and the psyche-brain- immune network in psychiatric diseases based on psychoneuroendocrineimmu- nology: A concise review. Annals of the New York Academy of Sciences, 1437(1), 31–42. https://doi.org/10.1111/nyas.13728 Bratman, G. N., Anderson, C. B., Berman, M. G., Cochran, B., de Vries, S., Flanders, J., Folke, C., Frumkin, H., Gross, J. J., Hartig, T., Kahn, Ph, Jr, Kuo, M, Lawler, J. J, Levin, P. S, Lindahl, T, Meyer-Lindenberg, A, Mitchell, R, Ouyang, Z, Roe, J, . . . Daily, G. C. (2019). Nature and mental health: An ecosystem service perspective. Science Advances, 5(7), eaax0903. https://doi.org/10.1126/sciadv. aax0903 Bratman, G. N., Hamilton, J. P., Hahn, K. S., Daily, G. C., & Gross, J. J. (2015). Nature experience reduces rumination and subgenual prefrontal cortex activa- tion. Proceedings of the National Academy of Sciences of the United States of America, 112(28), 8567–8572. https://doi.org/10.1073/pnas.1510459112 Dadvand, P., Bartoll, X., Basagaña, X., Dalmau-Bueno, A., Martinez, D., Ambros, A., Cirach, M., Triguero-Mas, M., Gascon, M., Borrell, C., & Nieuwenhuijsen, M. J. (2016). Green spaces and general health: Roles of mental health status, Ezpeleta et al. 911 social support, and physical activity. Environment International, 91, 161–167. https://doi.org/10.1016/j.envint.2016.02.029 Dadvand, P., Hariri, S., Abbasi, B., Heshmat, R., Qorbani, M., Motlagh, M. E., Basagaña, X., & Kelishadi, R. (2019). Use of green spaces, self-satisfaction and social con- tacts in adolescents: A population-based CASPIAN-V study. Environmental Research, 168, 171–177. https://doi.org/10.1016/j.envres.2018.09.033 Dadvand, P., Nieuwenhuijsen, M. J., Esnaola, M., Forns, J., Basagaña, X., Alvarez- Pedrerol, M., Rivas, I., López-Vicente, M., De Castro Pascual, M., Su, J., Jerrett, M., Querol, X., & Sunyer, J. (2015). Green spaces and cognitive development in primary schoolchildren. Proceedings of the National Academy of Sciences of the United States of America, 112(26), 7937–7942. https://doi.org/10.1073/ pnas.1503402112 Dadvand, P., Pujol, J., Macià, D., Martínez-Vilavella, G., Blanco-Hinojo, L., Mortamais, M., Alvarez-Pedrerol, M., Fenoll, R., Esnaola, M., Dalmau-Bueno, A., López-Vicente, M., Basagaña, X., Jerrett, M., Nieuwenhuijsen, M. J., & Sunyer, J. (2018). The association between lifelong greenspace exposure and 3-dimensional brain magnetic resonance imaging in Barcelona schoolchildren. Environmental Health Perspectives, 126(2), 027012. https://doi.org/10.1289/ EHP1876 Dadvand, P., Tischer, C., Estarlich, M., Llop, S., Dalmau-Bueno, A., López- Vicente, M., Valentín, A., de Keijzer, C., Fernández-Somoano, A., Lertxundi, N., Rodriguez-Dehli, C., Gascon, M., Guxens, M., Zugna, D., Basagaña, X., Nieuwenhuijsen, M. J., Ibarluzea, J., Ballester, F., & Sunyer, J. (2017). Lifelong residential exposure to green space and attention: A population-based prospec- tive study. Environmental Health Perspectives, 125(9), 097016. https://doi. org/10.1289/EHP694 Donovan, G. H., Michael, Y. L., Gatziolis, D., & Hoyer, R. W. (2020). The rela- tionship between the natural environment and individual-level academic perfor- mance in Portland, Oregon. Environment and Behavior, 52(2), 164–186. https:// doi.org/10.1177/0013916518796885 Dougherty, D. D., Brennan, B. P., Stewart, S. E., Wilhelm, S., Widge, A. S., & Rauch, S. L. (2018). Neuroscientifically informed formulation and treatment planning for patients with obsessive-compulsive disorder. A review. JAMA Psychiatry, 75(10), 1081–1087. https://doi.org/10.1001/jamapsychiatry.2018.0930 Engemann, K., Pedersen, C. B., Agerbo, E., Arge, L., Børglum, A. D., Erikstrup, C., Hertel, O., Hougaard, D. M., McGrath, J. J., Mors, O., Mortensen, P. B., Nordentoft, M., Sabel, C. E., Sigsgaard, T., Tsirogiannis, C., Vilhjálmsson, B. J., Werge, T., Svenning, J.-C., & Horsdal, H. T. (2020). Association between childhood green space, genetic liability, and the incidence of schizophrenia. Schizophrenia Bulletin, 46(6), 1629–1637. https://doi.org/10.1093/schbul/sbaa058 Engemann, K., Pedersen, C. B., Arge, L., Tsirogiannis, C., Mortensen, P. B., & Svenning, J. C. (2018). Childhood exposure to green space – A novel risk- decreasing mechanism for schizophrenia? Schizophrenia Research, 199, 142– 148. https://doi.org/10.1016/j.schres.2018.03.026 912 Environment and Behavior 54(5) Engemann, K., Pedersen, C. B., Arge, L., Tsirogiannis, C., Mortensen, P. B., & Svenning, J.-C. (2019). Residential green space in childhood is associated with lower risk of psychiatric disorders from adolescence into adulthood. Proceedings of the National Academy of Sciences of the United States of America, 116(11), 5188–5193. https://doi.org/10.1073/pnas.1807504116 European Union. (2020). Mapping guide v6.1 for a European Urban Atlas. https:// land.copernicus.eu/user-corner/technical-library/urban_atlas_2012_2018_map- ping_guide_v6-1.pdf Ezpeleta, L., de la Osa, N., & Doménech, J. M. (2014). Prevalence of DSM-IV disorders, comorbidity and impairment in 3-year-old Spanish preschoolers. Social Psychiatry and Psychiatric Epidemiology, 49(1), 145–155. https://doi. org/10.1007/s00127-013-0683-1 Feng, X., & Astell-Burt, T. (2017). Residential green space quantity and quality and child well-being: A longitudinal study. American Journal of Preventive Medicine, 53(5), 616–624. https://doi.org/10.1016/j.amepre.2017.06.035 Garcia-Delgar, B., de la Serna, E., Lera, S., Varela, E., Andrés-Perpiñá, S., & Lázaro, L. (2018). Neuropsychological functioning in pediatric obsessive-compulsive disorder: One size does not fit all. Journal of Obsessive-Compulsive and Related Disorders, 16, 36–42. https://doi.org/10.1016/j.jocrd.2017.12.005 Geller, D. A. (2006). Obsessive-compulsive and spectrum disorders in children and adolescents. Psychiatric Clinics of North America, 29(2), 353–370. https://doi. org/10.1016/j.psc.2006.02.012 Goodman, R. (1997). The strengths and Difficulties Questionnaire: A research note. Journal of Child Psychology and Psychiatry, 38, 581–586. https://doi. org/10.1111/j.1469-7610.1997.tb01545.x Grisham, J. R., Anderson, T. M., & Sachdev, P. S. (2008). Genetic and environmental influences on obsessive-compulsive disorder. European Archives of Psychiatry and Clinical Neuroscience, 258(2), 107–116. https://doi.org/10.1007/s00406- 007-0789-0 Heyman, I., Fombonne, E., Simmons, H., Ford, T., Meltzer, H., & Goodman, R. (2001). Prevalence of obsessive–compulsive disorder in the British nationwide survey of child mental health. The British Journal of Psychiatry, 179, 324–329. https://doi.org/10.1192/bjp.179.4.324 Instituto Nacional de Estadística. (2020). Estadística del padrón continuo a 1 de enero de 2020. INE. https://www.ine.es/dyngs/INEbase/es/operacion.htm?c=Es tadistica_C&cid=1254736177012&menu=resultados&idp=1254734710990 Johnson, J. H., & McCutcheon, S. M. (1980). Assessing life stress in older chil- dren and adolescents: preliminary findings with the Life Events Checklist. In I. Saranson & C. D. Spielberger (Eds.), Stress and anxiety (pp. 111–125). Hemisphere Publishing Corporation. Kabisch, N., Alonso, L., Dadvand, P., & van Den Bosch, M. (2019). Urban natu- ral environments and motor development in early life. Environmental Research, 179, 108774. https://doi.org/10.1016/j.envres.2019.108774 Kelz, C., Evans, G. W., & Röderer, K. (2015). The restorative effects of redesign- ing the schoolyard: A multi-methodological, quasi-experimental study in rural Ezpeleta et al. 913 Austrian middle schools. Environment and Behavior, 47(2), 119–139. https://doi. org/10.1177/0013916513510528 Land Cover Land Use Change (NASA). (2015). P197r031_TC_2015, p198r031_ TC_2015. 30m pixel. https://lcluc.umd.edu/metadata/global-30m-landsat-tree- canopy-version-4 Landsat images for 2015 U.S. Geological Survey. (2015). LC08_L1TP_197031_20 150508_20170409_01_T1, LC08_L1TP_198031_20150413_20170409_01_T1. [Satellite Image: LANDSAT 8 OLI/TIRS.] 30m pixel. Earth Explorer Viewer. https://www.usgs.gov/landsat-missions/landsat-headlines-2015 Larson, L. R., Barger, B., Ogletree, S., Torquati, J., Rosenberg, S., Gaither, C. J., Bartz, J. M., Gardner, A., Moody, E., & Schutte, A. (2018). Gray space and green space proximity associated with higher anxiety in youth with autism. Health & Place, 53, 94–102. https://doi.org/10.1016/j.healthplace.2018.07.006 Lee, K. S., Kim, B. N., Cho, J., Jang, Y. Y., Choi, Y. J., Lee, W. S., Han, C., Bae, H. J., Lim, Y. H., Kim, J. I., Shin, C. H., Lee, Y. A., & Hong, Y. C. (2021). Associations between surrounding residential greenness and intelligence quo- tient in 6-year-old children. The Science of the Total Environment, 759, 143561– 143561. https://doi.org/10.1016/j.scitotenv.2020.143561 Lee, M., Kim, S., & Ha, M. (2019). Community greenness and neurobehavioral health in children and adolescents. The Science of the Total Environment, 672, 381–388. https://doi.org/10.1016/j.scitotenv.2019.03.454 Liao, J., Yang, S., Xia, W., Peng, A., Zhao, J., Li, Y., Zhang, Y., Qian, Z., Vaughn, M. G., Schootman, M., Zhang, B., & Xu, S. (2020). Associations of exposure to green space with problem behaviours in preschool-aged children. International Journal of Epidemiology, 49(3), 944–953. https://doi.org/10.1093/ije/dyz243 Liao, J., Zhang, B., Xia, W., Cao, Z., Zhang, Y., Liang, S., Hu, K., Xu, S., & Li, Y. (2019). Residential exposure to green space and early childhood neurode- velopment. Environment International, 128, 70–76. https://doi.org/10.1016/j. envint.2019.03.070 Liu, J., Cui, Y., Yu, L., Wen, F., Wang, F., Yan, J., Yan, C., & Li, Y. (2021). Long- term outcome of pediatric obsessive-compulsive disorder: A meta-analysis. Journal of Child and Adolescent Psychopharmacology, 31, 95–101. https://doi. org/10.1089/cap.2020.0051 Lohr, S. L. (2019). Sampling design and analysis (2nd ed.). Chapman and Hall/CRC. Markevych, I., Schoierer, J., Hartig, T., Chudnovsky, A., Hystad, P., Dzhambov, A. M., de Vries, S., Triguero-Mas, M., Brauer, M., Nieuwenhuijsen, M. J., Lupp, G., Richardson, E. A., Astell-Burt, T., Dimitrova, D., Feng, X., Sadeh, M., Standl, M., Heinrich, J., & Fuertes, E. (2017). Exploring pathways linking greenspace to health: Theoretical and methodological guidance. Environmental Research, 158, 301–317. https://doi.org/10.1016/j.envres.2017.06.028 Markevych, I., Tiesler, C. M., Fuertes, E., Romanos, M., Dadvand, P., Nieuwenhuijsen, M. J., Berdel, D., Koletzko, S., & Heinrich, J. (2014). Access to urban green spaces and behavioural problems in children: Results from the giniplus and lisaplus studies. Environment International, 71, 29–35. https://doi.org/10.1016/j. envint.2014.06.002 914 Environment and Behavior 54(5) McCormick, R. (2017). Does access to green space impact the mental well-being of children: A systematic review. Journal of Pediatric Nursing, 37, 3–7. https://doi. org/10.1016/j.pedn.2017.08.027 McEachan, R. R. C., Yang, T. C., Roberts, H., Pickett, K. E., Arseneau-Powell, D., Gidlow, C. J., Wright, J., & Nieuwenhuijsen, M. (2018). Availability, use of, and satisfaction with green space, and children’s mental wellbeing at age 4 years in a multicultural, deprived, urban area: Results from the Born in Bradford cohort study. Lancet Planetary Health, 2(6), E244–E254. https://doi.org/10.1016/ s2542-5196(18)30119-0 McEwen, B. S., & Stellar, E. (1993). Stress and the individual. Mechanisms lead- ing to disease. Archive of Internal Medicine, 153(18), 2093–2101. https://doi. org/10.1001/archinte.153.18.2093 Monzani, B., Vidal-Ribas, P., Turner, C., Krebs, G., Stokes, C., Heyman, I., Mataix- Cols, D., & Stringaris, A. (2020). The role of paternal accommodation of paedi- atric OCD symptoms: Patterns and implications for treatment outcomes. Journal of Abnormal Child Psychology, 48(10), 1313–1323. https://doi.org/10.1007/ s10802-020-00678-9 NASA earth observatory. (2000). Measuring vegetation (NDVI & EVI). https:// earthobservatory.nasa.gov/features/MeasuringVegetation/measuring_vegeta- tion_2.php Negreiros, J., Best, J. R., Yamin, D. F., Belschner, L., Lin, S., & Stewart, S. E. (2020). Test-based versus parent ratings of executive function in pediatric obsessive- compulsive disorder. Journal of Obsessive-Compulsive and Related Disorders, 24. Article 100495. https://doi.org/10.1016/j.jocrd.2019.100495 Ornstein, T. J., Arnold, P., Manassis, K., Mendlowitz, S., & Schachar, R. (2010). Neuropsychological performance in childhood OCD: A preliminary study. Depression and Anxiety, 27(4), 372–380. https://doi.org/10.1002/da.20638 Putra, I. G. N. E., Astell-Burt, T., Cliff, D. P., Vella, S. A., & Feng, X. (2021). Association between green space quality and prosocial behaviour: A 10-year multilevel longitudinal analysis of Australian children. Environmental Research, 196, 110334–110334. https://doi.org/10.1016/j.envres.2020.110334 Qi, J., Chehbouni, A., Huete, A. R., Kerr, Y. H., & Sorooshian, S. (1994). A modified soil adjusted vegetation index. Remote Sensing of Environment, 48(2), 119–126. https://doi.org/10.1016/0034-4257(94)90134-1 Rapoport, J. L., Inoff-Germain, G., Weissman, M. M., Greenwald, S., Narrow, W. E., Jensen, P. S., Lahey, B. B., & Canino, G. (2000). Childhood obsessive-com- pulsive disorder in the NIMH MECA study: Parent versus child identification of cases. Methods for the epidemiology of child and Adolescent Mental Disorders. Journal of Anxiety Disorders, 14(6), 535–548. https://doi.org/10.1016/s0887- 6185(00)00048-7 Rapoport, J. L., & Shaw, P. (2015). Obsessive compulsive disorder. In A. Thapar, D. S. Pine, J. F. Leckman, S. Scott, M. J. Snowling, & E. Taylor (Eds.), Rutter’s child and adolescent psychiatry (6th ed., pp. 841–857). Wiley. Ribeiro, A. I., Tavares, C., Guttentag, A., & Barros, H. (2019). Association between neighbourhood green space and biological markers in school-aged children. Ezpeleta et al. 915 Findings from the Generation XXI birth cohort. Environment International, 132, 105070. https://doi.org/10.1016/j.envint.2019.105070 Rigolon, A., Browning, M. H. E. M., McAnirlin, O., & Yoon, H. V. (2021). Green space and health equity: A systematic review on the potential of green space to reduce health disparities. International Journal of Environmental Research and Public Health, 18(5), 2563. https://doi.org/10.3390/ijerph18052563 Scahill, L., Riddle, M. A., McSwiggin-Hardin, M., Ort, S. I., King, R. A., Goodman, W. K., Cicchetti, D., & Leckman, J. F. (1997). Children’s Yale-Brown Obsessive Compulsive Scale: Reliability and validity. Journal of the American Academy of Child and Adolescent Psychiatry, 36(6), 844–852. https://doi. org/10.1097/00004583-199706000-00023 Schniering, C. A., Hudson, J. L., & Rapee, R. M. (2000). Issues in the diagnosis and assessment of anxiety disorders in children and adolescents. Clinical Psychology Review, 20(4), 453–478. https://doi.org/10.1016/s0272-7358(99)00037-9 Sexton, J. O., Song, X.-P., Feng, M., Noojipady, P., Anand, A., Huang, C., Kim, D.-H., Collins, K. M., Channan, S., DiMiceli, C., & Townshend, J. R. (2013). Global, 30-m resolution continuous fields of tree cover: Landsat-based rescal- ing of MODIS vegetation continuous fields with lidar-based estimates of error. International Journal of Digital Earth, 6(5), 427–448. https://doi.org/10.1080/1 7538947.2013.786146 Silverman, W. K., Albano, A. M., & Barlow, D. H. (1994). Anxiety Disorders inter- view schedule for DSM-IV: Child and Parent versions. Graywind Publications Inc. Spence, S. H. (1998). A measure of anxiety symptoms among children. Behaviour Research and Therapy, 36(5), 545–566. Spence, S. H. (n.d). Spence Children’s Anxiety Scale. T-Scores - Using the SCAS- Parent T-Scores. SCAS PARENT: T-Scores - The Spence Children’s Anxiety Scale (scaswebsite.com) StataCorp. (2019). Stata. Vanaken, G.-J., & Danckaerts, M. (2018). Impact of green space exposure on chil- dren’s and adolescents’ mental health: A systematic review. International Journal of Environmental Research and Public Health, 15(12), 2668. https://doi. org/10.3390/ijerph15122668 van Grootheest, D. S., Cath, D. C., Beekman, A. T., & Boomsma, D. I. (2005). Twin studies on obsessive-compulsive disorder: A review. Twin Research and Human Genetics, 8(5), 450–458. https://doi.org/10.1375/183242705774310060 Whiteside, S. P. H., Gryczkowski, M. R., Biggs, B. K., Fagen, R., & Owusu, D. (2012). Validation of the Spence Children’s Anxiety Scale’s obsessive compul- sive subscale in a clinical and community sample. Journal of Anxiety Disorders, 26(1), 111–116. https://doi.org/10.1016/j.janxdis.2011.10.002 Yang, B. Y., Zeng, X.-W., Markevych, I., Bloom, M. S., Heinrich, J., Knibbs, L. D., Dharmage, S. C., Lin, S., Jalava, P., Guo, Y., Jalaludin, B., Morawska, L., Zhou, Y., Hu, L.-W., Yu, H.-Y., Yu, Y., & Dong, G.-H. (2019). Association between greenness surrounding schools and kindergartens and attention-deficit/hyper- activity disorder in children in China. Jama Network Open, 2(12), e1917862. https://doi.org/10.1001/jamanetworkopen.2019.17862 916 Environment and Behavior 54(5) Zhang, J., Yu, Z., Zhao, B., Sun, R., & Vejre, H. (2020). Links between green space and public health: A bibliometric review of global research trends and future prospects from 1901 to 2019. Environmental Research Letters, 15(6), 063001. https://doi.org/10.1088/1748-9326/ab7f64 Zhang, Y., Mavoa, S., Zhao, J., Raphael, D., & Smith, M. (2020). The association between green space and Adolescents’ mental well-being: A systematic review. International Journal of Environmental Research and Public Health, 17(18), 6640. https://doi.org/10.3390/ijerph17186640 Author Biographies Lourdes Ezpeleta, Ph.D., is a Professor at the Universitat Autònoma de Barcelona. Her major research interest are developmental psychopathology, and risk and protec- tive factors in child and adolescent psychopathology. José Blas Navarro, Ph.D., is an associate Professor of Methods in Health Science at the Universitat Autònoma de Barcelona. His research is about multilevel modeling applied to psychopathology. Lucía Alonso is a Biostatistician at Barcelona Institute for Global Health (ISGlobal). Núria de la Osa, Ph.D., is an associate Professor at the Universitat Autònoma de Barcelona. Her major research interest is psychological assessment. Albert Ambrós is a Geographic Information Systems (GIS) technichian at the ISGlobal. His research is to assess the impact of the built environment on people’s wellbeing and health. Mònica Ubalde is a biologist by training and a Public Health epidemiologist. She currently develops her research work at the Urban Planning, Environment and Health Initiative, within the Air pollution and Urban Environment Program at ISGlobal. Eva Penelo, Ph.D., is a Professor at the Universitat Autònoma de Barcelona. Her major research interests are applied psychometrics and growth curve modeling Payam Dadvand, MD, PhD, is an environmental epidemiologist, currently working as an associate research professor at the Barcelona Institute for Global Health (ISGlobal). His research interests include, among others, studying the effects of natu- ral environment on child mental and physical health and development.