Research Design
The Green Schoolyards Project consisted of two separate studies: a) a serial cross-sectional study focused on physical sites within parks, and b) a prospective cohort study focused on individual students. The cross-sectional study was designed to examine the associations between heat index, children’s physical activity, and interaction with green features for multiple sites (e.g., playground, track, and soccer field) per park on multiple times per day. The cohort study followed a selection of students from each school affiliated with these parks to assess the impact of green features—the amount of which differs per park—on student’s physical activity levels during recess for different heat index conditions, along with their connection to nature, social-emotional learning skills, misconduct, and standardized test scores.
Conceptual Model
Research design of the Green Schoolyards Project was based on the social-ecological model and attention restoration theory. The social-ecological model states multiple levels of influence—individual, social, environmental, and policy—impact health behaviors (e.g., physical activity), and these influences interact across levels to impact behavior (22). Effective health interventions focus on behavior-specific influences, and intervene at multiple levels of influence. In general, attention restoration theory proposes exposure to natural environments replenishes the ability to concentrate, a cognitive resource that can be depleted (23, 24). Natural settings are considered restorative because these environments provide a) respite from everyday stressors, b) space that feels extensive, c) objects that fascinate and require little concentration, and d) intrinsic compatibility with humans (24).
From these theoretical underpinnings, we developed a conceptual model of the Green Schoolyards Project that illustrates the concepts of interest and relations between these concepts (Figure 1). The model provides our logic for testing green features in school parks as a strategy for improving health, psychosocial, and academic outcomes of children. We posited that trees and gardens will decrease outdoor heat and that nature trails, which surround gardens but are not directly vegetated, will not significantly impact outdoor heat. When subjected to outdoor heat, children will engage in less physical activity due to thermal discomfort. Trees, gardens, and nature trails will inherently provide more opportunities for children to physically interact with nature, yet outdoor heat will have a mixed impact on these physical interactions: children will seek natural settings providing protection from outdoor heat (e.g., park sites under tree shade) while not seeking and potentially avoiding natural settings without heat protection (e.g., unshaded park sites). The relation between children’s physical interactions with nature and their physical activity levels will depend on the interaction with nature (e.g., climbing trees versus resting under tree shade) and whether the setting is intended for physical activity (e.g., playground under tree canopy versus picnic tables under tree canopy). Lastly, children’s physical interactions with nature will offer opportunities for beneficial outcomes, such as mentally connecting with nature and restoring ability to concentrate. Improved concentration will increase children’s social-emotional learning skills and standardized tests scores, and decrease their misconduct at school.
Figure 1. Conceptual model of Green Schoolyards Project.
Project Sites and Green Features
Three elementary school parks within a school district in Central Texas, US, were used for this study. The project was a comparative analysis between similar parks with different levels of green features. Three schools met our initial selection criteria: serving populations greater than 85% economically disadvantaged Latinx; located in ZIP codes with low Nature Factor scores; joint-use agreements between the school district and the city Parks and Recreation Department permitting the surrounding community to use parks after school hours; and equivalent park features (e.g., playgrounds, soccer fields, running tracks, and basketball courts).
The selection criterion of Nature Factor is defined as the sum of Nature Factor Ratings of all parks within a ZIP code (25). Nature Factor Rating is the sum of four park-level ratings: park acreage rating, Trust for the Public Land land use rating, National Recreation and Parks Association park status rating, and tree canopy rating (26). High values for Nature Factor Ratings (e.g., high park acreage, designed lands, open park status, and high levels of tree canopy) correspond to high Nature Factor scores (i.e., higher levels of nature present in that ZIP code). The three selected schools are in ZIP codes with Nature Factor scores of 121, 118, and 198, respectively, which are relatively low compared to those of other ZIP codes (n = 53; min. = 0; max. = 712; mean = 150; standard deviation = 143).
The three school parks were characterized by different profiles of green features: the “intervention park” had added green features (i.e., trees, wildflower meadow, and nature trail); the “low-green park” had relatively low amounts of historical green features (i.e., trees); and the “high-green park” had relatively high amounts of historical green features (i.e., trees, wildlife habitat garden, and nature trail).
Landscape vegetation—including trees and gardens—was planted in each park to provide visual enhancement; learning, play, and recreation opportunities; and a calming and welcoming environment. The school district is responsible for maintaining trees, and individual schools and their communities are responsible for maintaining gardens and nature trails. The wildflower meadow at the intervention park is a 1,383m2 pollinator-supporting grassland consisting of dense, tall wildflowers and grasses native to Texas (27). The meadow includes signage and a weaving, 20 meter-long, 1.5 meter-wide, nature trail—a secondary or tertiary path that branches from a primary track, trail, or pathway and composed of materials such as decomposed granite, compacted dirt, and stepping stones. To encourage use of the meadow, a nature themed-story in English- and Spanish-language was installed on temporary signage placed along the trail. The wildlife habitat garden at the high-green park is a 247m2 dense cover of native grasses and shrubs that serves as a haven for local and migratory species by providing food, cover, and places to raise young (28). This garden includes a nature trail that is nearly identical to the trail in the wildflower meadow at the intervention park. For the Green Schoolyards Project, both wildflower meadow and wildlife habitat garden are categorized as gardens because of similar characteristics. The low-green park had no gardens or nature trails.
For each school park, we calculated tree canopy cover using i-Tree Canopy, a publicly available tool that uses random point sampling to estimate the percentage of tree canopy cover for a predefined area (29). Although trees at the intervention park were more abundant and evenly distributed than trees at parks at the other schools, tree canopy cover was only 8.5% (standard error = 1.97) because trees planted were saplings. The low-green park had 11.5% (standard error = 2.26) tree canopy cover, with most trees clustered in the far northwest corner. The high-green park had 22.5% (standard error = 2.95) tree canopy cover, with relatively large trees on the park’s periphery.
Project Period
We designed data collection to take place on 18 days over the fall semester in 2019, which will be duplicated in 2020 for a total of 36 study days. Each year, study days consist of two September weeks (i.e., five weekdays and one weekend day per week) and one November week (i.e., five weekdays and one weekend day). We selected September and November because these months have historically high and moderate temperature conditions, respectively: the weather station at the city’s major airport recorded monthly mean average air temperatures of 26.4°C in September and 15.3°C in November from 2009 to 2018 (30). Prior to undertaking any project activities, we received approval of project protocols by the institutional review board at The University of Texas Health Science Center at Houston (HSC-SPH-19-0502) and the school district. We also received informed consent from participants’ parents and written assent from study participants.
Measurement of Heat Index
We measured heat index—the combination term for air temperature and relative humidity that captures what the temperature feels like (31)—by semi-permanently installing 10 HOBO MX2302A external air temperature/relative humidity sensor data loggers (Onset Computer Corporation, MA) at each park. Previous studies have used comparable networks of in situ sensors to monitor microclimatic conditions of a given area (32-34). Measurement of near-surface air temperatures is advantageous over the use of land surface temperatures as a proxy for air temperatures, as research has shown land surface temperatures are not directly comparable to air temperatures (35, 36). In situ measurement of air temperatures has been found to be more useful for estimating short-term, actual temperature exposures than using land surface temperatures or the percentage of impervious surface (37). Designed for outdoor use, the MX2302A model collects air temperature (±0.2°C from 0 to 70°C) and relative humidity data (±2.5% from 10% to 90%), and is configured to wirelessly link with the free HOBOmobile app on a cell phone or tablet (38), permitting efficient collection of air temperature and relative humidity data by project members.
Before deploying HOBO sensors, we encased each sensor in an RS3-B solar radiation shield (Onset Computer Corporation, MA), which results in improved temperature measurement accuracy by protecting the sensor from absorption of incoming solar radiation and resultant heat gain. In another attempt to improve temperature measurement accuracy, we attached the radiation shield (encasing the external sensor) and data logger to a 2x2” piece of weather-treated lumber, which serves as a physical buffer minimizing heat transfer between the sensor and the installation surface in the park, such as a metal swing set pole (Figure 2). We programmed the sensors to record air temperature and relative humidity every five minutes, consistent with previous studies (34, 39, 40).
Figure 2. HOBO MX2302A external air temperature/relative humidity sensor data logger on swing set at low-green park.
For each park, we selected 10 sites for HOBO sensors based on land cover (e.g., grass, pavement, and mulch); land use (e.g., soccer field, basketball court, and playground); comparability across parks; and even spatial coverage (Figure 3). To include a highly impervious area for comparison within park sites, we installed one of the 10 sensors at each park’s parking lot, just outside park boundaries. To capture air temperature and relative humidity experienced by humans, we installed sensors at two meters above ground level, similar to previous studies (33, 34, 39). To promote community awareness of our project and deter vandalism, we attached a small laminated tag with a description of the sensor and our contact information to each sensor. Sensors were installed the day before a study week, and removed at the end of each study week.
Measurement of Green Features
We identified the location, type, and quantity of green features using four-band, 60cm orthoimagery taken in November 2018 by the US Department of Agriculture’s National Agriculture Imagery Program (41). Within a geographic information system (ArcGIS 10.6.1, ESRI, Redlands, CA, USA), we digitized polygons of trees, gardens, and nature trails, an established technique deemed appropriate for the relatively small park areas (i.e., intervention = 21,448 m2; low-green = 27,923 m2; high-green = 16,187 m2) (42).
Direct Observation of Parks
The cross-sectional study utilized the System for Observing Play and Recreation in Communities (SOPARC), a validated direct observation tool for assessing the conditions and users of park sites (43). Following SOPARC protocol, we divided each park into target areas intended for physical activity, such as basketball courts and soccer fields (Figure 3). On study days during school (i.e., 7:00 and 12:00) and after school (i.e., 16:00 and 18:00), study staff—in pairs for interrater reliability—administered SOPARC by walking from target area to target area and recording what they observed. We adapted SOPARC to measure physical activity levels of children aged 1–12 years old and these children’s interactions with green features. Although previous research has used direct observation to examine the influence of nature on children’s play (44), no research employs SOPARC to quantify the number of children’s interactions with different green features at multiple park sites.
On study days for each target area, trained staff used SOPARC to record the date and time of observation and target area conditions (i.e., whether area was accessible, usable, equipped, supervised, dark, empty, and organized); then, if applicable, to record the primary physical activity (e.g., playing basketball) of female and male children; scan for the physical activity levels of female and male children; and lastly scan for the number of female and male children interacting with green features (i.e., no interaction, under tree canopy or touching tree, interacting with garden, and on nature trail). To employ SOPARC, staff used the iSOPARC application on an electronic tablet for scan counts and input data on a data collection form (see Additional File 1).
Figure 3. HOBO sensors and SOPARC target areas at (A1–A2) intervention, (B1–B2) low-green, and (C1–C2) high-green parks.
Cohort Study Sample
For the cohort study, we recruited 40 3rd and 40 4th grade students per school over two years, to achieve a final sample size of 30 students per grade, after accounting for attrition. From mid-August 2019 through early September 2019, we recruited participants by convenience sampling. Participant incentives were a total of $35 US dollars/year worth of supermarket gift cards (i.e., $10 for each September study week and $15 for each November study week).
Measurement of Geographic Location and Physical Activity
On study days during recess (i.e., 30-minute period of unstructured play under teacher supervision), the cohort sample wore elastic belts around their waist equipped with a Qstarz BT-Q1000XT Global Positioning System (GPS) device (Qstarz Intl Co., Taipei, Taiwan) and an Actigraph wGT3X-BT accelerometer (ActiGraph LLC, FL) to measure geographic location and physical activity levels, over time (45, 46). We set sampling rates of 15 seconds for GPS devices and accelerometers (47, 48). For the separate recess periods per grade, belt distribution began five minutes before recess start, and belt collection occurred once teachers signaled recess end.
Time-Matching of Geographic Location, Physical Activity, and Heat Index
Data from GPS devices, accelerometers, and HOBO sensors will be time-matched, allowing us to know a student’s location, student’s physical activity intensity level, and heat index at that location at 15-second intervals throughout recess. The location of green features will be joined to the time-matched device data, within GIS. Although previous studies have matched children’s geographic location and physical activity levels over time (45-48), this study enables assessment of a child’s experienced heat index and physical activity level at any particular location.
Measurement of Connection to Nature, Psychosocial and Academic Outcomes, and School Policies
For baseline data collection in November 2019, we collected data on the cohort sample and school policies from three sources. First, we administered aloud a written survey—in both English and Spanish language—to the cohort sample, asking them about their connection to nature using two adapted instruments: Inclusion of Nature with Self and Connection to Nature Index (see Additional File 2) (49-51). Second, the school district provided individual data for each student on sociodemographic characteristics, social-emotional learning skills (from a student climate survey), number of disciplinary actions for misconduct, and standardized test scores. Lastly, we distributed an annual survey—adapted from a previous study (52)—to ask school principals about policies impacting park access, greening at school parks, and student physical activity.
Statistical Analysis of Baseline Data
For the selection of baseline findings shared within, we performed analysis over several steps. To develop thermal profiles for each park, we first calculated heat index—from air temperature and relative humidity data recorded by HOBO sensors—using validated equations utilized by US National Weather Service (53, 54). Heat index is measured in degrees Fahrenheit (55). Within GIS, we used inverse distance weighting to create spatially continuous thermal profiles for each park. A common interpolation method in urban heat island measurement (56), inverse distance weighting permits estimation of unsampled heat index values between HOBO sensors by averaging sampled heat index values from sensors surrounding each prediction location. We used SOPARC data to understand how children interact with trees during time periods with different temperature conditions, summing observed counts of children under tree canopy or touching trees by sex of child, park, and study period.