Access the full text.
Sign up today, get DeepDyve free for 14 days.
Background: Research suggests that exposure to the natural environment can improve mood, however, current reviews are limited in scope and there is little understanding of moderators. Objective: To conduct a comprehensive systematic review and meta-analysis of the evidence for the effect of short-term exposure to the natural environment on depressive mood. Methods: Five databases were systematically searched for relevant studies published up to March 2018. Risk of bias was evaluated using the Cochrane Risk of Bias (ROB) tool 1.0 and the Risk of Bias in Non- Randomised Studies of Interventions (ROBINS-I) tool where appropriate. The Grades of Recommendation, Assessment, Development, and Evaluation (GRADE) approach was used to assess the quality of evidence overall. A random-effects meta-analysis was performed. 20 potential moderators of the effect size were coded and the machine learning-based MetaForest algorithm was used to identify relevant moderators. These were then entered into a meta-regression. Results: 33 studies met the inclusion criteria. Effect sizes ranged from −2.30 to 0.84, with a pooled effect size of γ = -0.30 95% CI [-0.50 to -0.10]. However, there was significant residual heterogeneity between studies and risk of bias was high. Type of natural environment, type of built environment, gender mix of the sample, and region of study origin, among others, were identified as relevant moderators but were not significant when entered in a meta-regression. The quality of evidence was rated very low to low. An assessment of publication bias was inconclusive. Conclusions: A small effect was found for reduction in depressive mood following exposure to the natural environment. However, the high risk of bias and low quality of studies limits confidence in the results. The variation in effect size also remains largely unexplained. It is recommended that future studies make use of reporting guidelines and aim to reduce the potential for bias where possible. Keywords: systematic review, meta-analysis, natural environment, built environment, green space, depression Highlights: • A systematic review and meta-analysis is conducted of depressive mood following natural environment exposure • Studies published up until March 2018 were searched; 33 studies met the eligibility criteria • A small effect size for reduction in depressive mood following exposure is found • Studies are found to be highly biased and of low or very low quality • No significant moderators of the effect size are identified 1 Introduction Depression is understood to have a lifetime prevalence of 10.8% among the global population (Lim et al., 2018). It is a leading contributor to the global disease burden (Ferrari et al., 2013), and at its worst can lead to suicide (Hawton, Casañas I Comabella, Haw, & Saunders, 2013). Characterised by depressed mood, loss of interest or enjoyment, and lack of energy, depression is estimated to affect over 300 million people worldwide (World Health Organization, 2017). The natural environment is increasingly recognised as a potential buffer to poor mental health (Gascon et al., 2015; Houlden, Weich, de Albuquerque, Jarvis, & Rees, 2018; Rautio, Filatova, Lehtiniemi, & Miettunen, 2017; Van Den Berg et al., 2015). There are a number of existing theories linking green space and health (Hartig, Mitchell, de Vries, & Frumkin, 2014; Nieuwenhuijsen, Khreis, Triguero-Mas, Gascon, & Dadvand, 2017); two focus on the effects on mental health specifically: attention restoration theory (Kaplan, 1995) and stress reduction theory (Ulrich et al., 1991). Attention restoration theory proposes that the natural environment provides a ‘soft fascination’ whereby a person can pay attention without effort. Stress reduction theory suggests that the presence of nature brings about a psycho- evolutionary response related to safety and survival, and therefore produces positive emotions. These pathways have received support in the literature: a number of cross-sectional studies (Beyer, Kaltenbach, Szabo, Bogar, & Nieto, 2014; Gascon et al., 2018; McEachan et al., 2015; Reklaitiene et al., 2014) and a few longitudinal studies (Alcock, White, Wheeler, Fleming, & Depledge, 2014; Astell- Burt, Mitchell, & Hartig, 2014) have found a relationship between increased surrounding green space and reduced risk of depression, and suicide also (Helbich et al., 2018; Min et al., 2017). Moreover, a recent systematic review of 28 studies found limited evidence of a causal relationship between quantity of and access to surrounding residential green space and mental health in adults (Gascon et al., 2015). However, it is unclear in these studies whether the association between green space and mental health is the result of use of green space or via another mechanism (Van den Bosch & Meyer-Lindenberg, 2019). Indeed, viewing green space from an indoor environment has been shown to have beneficial mental health effects (Ulrich, 1984; Ulrich et al., 1991). There is also a risk of self-selection bias, whereby associations might be attributed to those with fewer mental health problems moving into greener neighbourhoods. Nevertheless, studies examining the relationship between direct exposure to the natural environment and mood have shown improvements after only a short period of time (Barton & Pretty, 2010; Shanahan et al., 2016; Shanahan, Fuller, Bush, Lin, & Gaston, 2015). Barton and Pretty (2010) assessed the effect on mood from exercising in nature and reported that the greatest benefit to mood came following 5 minutes of exercise, with diminishing returns thereafter. Moreover, Shanahan et al. (2016) found that the odds of reporting depression were significantly lower when respondents visited green space for an average of 30 minutes or more. In this way, short-term visits to the natural environment might represent a cheap and feasible intervention to improve mood. A clear understanding of the evidence base is then necessary in order to develop effective interventions. Previous reviews of this topic are restricted in scope. For example, Lee et al. (2017) and Oh et al. (2017) examine the effect on depressive symptoms and health and well-being respectively, but both consider exposure to a specific type of natural environment only – forests. This disregards current research that considers the impact of other natural environments, such as parks (Song et al., 2013, 2014), agricultural land (Lee, Park, Ohira, Kagawa, & Miyazaki, 2015) and streetscape greenery (Helbich et al., 2019). Hansen, Jones, and Tocchini (2017) also review the effects of forest therapy on physical and mental health, but only include studies from Japan or China. While the concept of forest-bathing originated in Japan, resulting in a large proportion of this research coming from there, the number of studies from Europe and North America that also examine the effects of exposure to the natural environment is growing (Stigsdotter, Corazon, Sidenius, Kristiansen, & Grahn, 2017; Triguero-Mas et al., 2017). Lastly, McMahan and Estes (2015) include all types of natural environment and investigate its effect on positive and negative mood, but exclude all studies that are not randomised controlled studies. In sum, previous reviews are narrow in focus: an overall understanding of how short-term exposure to the natural environment affects depressive mood is missing in the literature. In this review the findings of previous reviews are built on and the scope extended by including studies of varying designs and place of origin, and also adopting a wider understanding of ‘natural environment’. All types of green space are considered – any open land with natural vegetation, planned or otherwise. A more complete picture of the effect of the natural environment on depressive mood is therefore presented. Blue space is not included in this review as a systematic review that considers blue space and mental health outcomes has recently been published (Gascon, Zijlema, Vert, White, & Nieuwenhuijsen, 2017). Furthermore, a meta-analysis is conducted and an exploratory approach to moderator analysis is used. The machine learning-based MetaForest algorithm is applied to identify relevant moderators of the effect size (van Lissa, 2017). Moderators entered include age, gender mix of sample, present health condition, type of natural environment and region of study origin. This technique has not previously been applied in the green space-health literature, and therefore provides a novel contribution to a rapidly expanding field of research. This systematic review and meta-analysis aims to summarise and critically appraise the evidence on the effect of short-term exposure to the natural environment on depressive mood. A secondary aim was to identify any potential moderators of this relationship. The final aim was to evaluate the quality of the evidence available. 2 Methods The protocol was registered on PROSPERO (Roberts, Hagedoorn, Kellar, & Helbich, 2018). The review followed PRISMA guidelines (Moher, Liberati, Tetzlaff, Altman, & Group, 2009) (for checklist see Table S1 in the supplementary materials). 2.1 Eligibility criteria Only original peer-reviewed research was eligible; abstracts, conference proceedings and grey literature were excluded. All geographical areas were eligible, but only references written in English were included. 2.2 PECO statement A PECO was developed in order to inform the development of search strategies and guide the screening of relevant studies. Population: Any human adult population Exposure: Exposure to the natural environment Comparator: Exposure to the built environment Outcome: Depressive mood In terms of populations considered, any adult population was eligible, regardless of physical or mental health status. Exposure was defined by placement of participants in direct contact with the exposure environment, within the context of a (randomised or non-randomised) trial (e.g. crossover, parallel group, factorial). Exposure duration was not limited, but based on a previous systematic review of the effects of short-term, direct exposure to the natural environment on health and well-being, it was expected that a single exposure would last approximately one hour with exceptions for those that applied repeated exposures (Bowler, Buyung-Ali, Knight, & Pullin, 2010). Representations of an environment using virtual reality, pictures or video were excluded. Environments were deemed as ‘natural’ if they were defined by a high level of greenery and had not been extensively modified by human activity. In contrast, the built environment was defined as a predominantly man-made environment with a low level of greenery. Studies where participants were exposed to more than two environments but included the natural and built environment were eligible, however, only data from the natural and built environment were included in the meta-analysis. The primary outcome was depressive mood. This should be measured quantitatively either by the researcher(s) or the participant. The measure must refer to the current emotional state, rather than depressive mood over a longer preceding period. It might be measured independently, or as part of a wider mood or affect assessment. Measurements could be recorded pre- and post-exposure, or post- exposure only. Studies that measured well-being or quality of life were excluded as they were understood to be concepts distinct from depression. 2.3 Search strategy A literature search was conducted on five databases: Medline, PsychINFO, Embase, Scopus and Web of Science from inception to March 2018. Search terms were related to the natural environment (such as ‘natural environment’, ‘green space’, ‘open space’ or ‘park’) and to depression or depressive mood (such as ‘depression’, ‘depressive symptoms’, ’mood disorder’, and ‘mental health’). The full search strategies are available in the supplementary material. 2.4 Study selection Records from each database were downloaded and merged in Endnote. Duplicates were removed. The titles and where available, abstracts, of the remaining references were screened independently by two researchers according to the PECO statement and eligibility criteria. Percent agreement was 99%. All eligible references were then evaluated at the full-text level. Full papers were screened independently by two authors using the same criteria. Percent agreement was 98%. Reference lists of included studies were also hand-searched for relevant studies. 2.5 Data extraction Relevant data were extracted by the first author using a standardised form. This included information on the study design, participant information, intervention and control environments, procedural descriptions, outcomes and outcome measures, results and conclusions. All data extraction forms were checked by a second researcher, and any disagreement was resolved through discussion until consensus was reached. 2.6 Quality assessment The Cochrane Risk of Bias (ROB) Tool 1.0 (Higgins et al., 2011) was used to assess included randomised studies. The tool gives an overall risk of bias for randomised trials by scoring them across seven domains: random sequence generation, allocation concealment, blinding of participants and personnel, blinding of outcome assessment, incomplete outcome data, selective reporting, and any other sources of bias. The Risk of Bias In Non-Randomised Studies of Interventions (ROBINS-I) (Sterne et al., 2016) tool was used for non-randomised studies. Studies are scored on seven domains: confounding, participant selection, classification of interventions, deviation from intended intervention, missing data, measurement of outcomes, and reporting bias. All included studies were independently evaluated by two researchers. Appraisals were discussed between the researchers until consensus was reached. Quality of evidence was assessed using the Grading of Recommendations Assessment, Development, and Evaluation (GRADE) guidelines (GRADE Working Group, 2004; Guyatt et al., 2008). These guidelines apply a set of predetermined domains that either increase or decrease the level of confidence in the evidence. Domains that reduce confidence in the evidence are: risk of bias, inconsistency of results, indirectness of evidence, imprecision, and publication bias. On the other hand, a large magnitude of effect, confounding that increases effect magnitude, and a dose-response gradient can increase confidence. Two researchers discussed the domains for each outcome until consensus was reached. 2.7 Data synthesis Narrative summary Studies are first narratively discussed in the context of the type of engagement performed by participants in the exposure environment: active engagement (e.g. walk, run), passive engagement (e.g. sit, stand), or a combination of both. Meta-analysis All studies were considered for inclusion in a meta-analysis. Two types of effect size were calculated, where appropriate: standardized mean difference (SMD) (Hedges' G) for randomized controlled experiments (Hedges, 1981), and standardized mean change rates (SMCR) for pretest- posttest designs (Morris & DeShon, 2002). Hedges’ G is recommended when sample sizes are below 20, and expresses the difference of the means in units of the pooled standard deviation. Furthermore, it can be interpreted similarly to Cohen’s d (e.g. 0.2 refers to a ‘small’ effect) (Fritz, Morris, & Richler, 2012). Mean (or mean change pre- and post-exposure), standard deviation (or standard error) and cell count (n) for all depression outcomes in each included study were extracted. In the first instance, data was extracted directly from the studies. When the data was not available, authors were contacted for further information. Where data could not be provided or contact could not be made, data was extracted from figures using an online ruler (A Ruler for Windows). Two studies did not report the standard deviation (Hartig, Mang, & Evans, 1991; Hartig, Nyberg, Nilsson, & Garling, 1999), therefore an estimate was derived by taking the average from similar papers that used the same outcome measure. One study released a corrigendum after the literature search had been conducted, therefore the author was contacted and the correct data used in the meta-analysis (Lee et al., 2019). It was not possible to retrieve data for two studies (Li et al., 2016; Perkins, Searight, & Ratwik, 2011). Only data from the natural environment and built environment were considered. Data from other settings e.g. blue space, passive lab setting were not included in meta-analysis. Two studies had multiple ‘green’ environments: Sonntag-Öström et al. (2014) examined differences across a spruce forest, a forest with a lake, and a forest with a rocky outcrop, and Tyrväinen et al. (2014) considered both a park and a forest. The first environment listed was selected for both studies. Moreover, it was felt that the lake and ‘rocky outcrop’ present in the alternative environments of Sonntag-Öström et al. (2014) may act as confounders. Only data from timepoints closest to the start and end time of the exposure were extracted; baseline or follow-up measures, or measurements taken during exposure, were not included in analysis. For cases of multiple exposure to the same environment, data was extracted from immediately before the first exposure, and immediately after the final exposure. Moderator analysis For each study a number of potential moderators were coded, for example study design, region of study origin, and mean age of the sample. For the full list of moderators, see Table 1. For studies that had multiple exposures to the same environment, exposure time was measured cumulatively. Table 1. Moderators coded for meta-analysis Moderator Potential codes Study design Crossover design Parallel groups Factorial design Region of study origin Asia Europe US Mean age of sample - Student sample Yes No Gender mix of sample Male Female Mixed Female (%) - Health condition of sample Healthy Poor mental health High blood pressure Chronic heart failure Chronic obstructive pulmonary disease Total exposure time (minutes) - Time between crossover to other environment (if Same day appropriate) Next day Next week Longer Type of natural environment Forest Park Biodiverse area Agricultural area Type of built environment Downtown Residential Other Number of natural environments - Number of built environments - Baseline measurement Yes No Measurement conducted at environment Post-exposure only Pre and post-exposure (either side multiple exposures) Pre and post-exposure for each exposure Measurement conducted during exposure Yes No Follow-up measurement Yes No Activity category Passive Active Mixed Primary depression measure - Secondary depression measure (if appropriate) - Missingness was very limited; three variables (gender mix, proportion of female, time between environments) had some missing values, ranging from 5-32%. Since complete data was required for analysis, single imputation was applied using a non-parametric missing value imputation by means of a random forest algorithm (Stekhoven & Bühlmann, 2012). The number of moderators coded (n=20) was large relative to the sample size. Consequently, including all moderators in a meta-regression risks overfitting the model. Therefore we apply the R package ‘metaforest’ (van Lissa, 2017); an exploratory approach to identify potentially relevant moderators in meta-analysis. The approach is based on the machine-learning algorithm ‘random forests’, which are robust to overfitting. First, the approach ranks moderators in terms of their importance in predicting the effect size. Second, partial dependence plots are produced which visualize the association of each moderator with the effect size, while accounting for the average effect of all other moderators. Lastly, a measure of ‘predictive performance’, or the ‘out-of-bag’ 𝑅 , is calculated for each moderator. In other words, an estimate is given of how much variance the moderators would explain if a new sample of data were provided. Moderators that consistently displayed negative variable importance (i.e., that showed a reduction in predictive performance) were dropped. Moderators that improved predictive performance were then entered into a linear meta-regression in order to understand their association with the effect size. For categorical variables, contrast coding is applied, such that the levels of one variable are compared with the mean of the subsequent levels. For ordinal variables, orthogonal polynomial coding is applied, and the linear, quadratic and cubic trends considered. Publication bias was first assessed by visual examination of funnel plots. Standard error was used as the measure of study size, plotted on the vertical axis, with effect estimates plotted on the horizontal axis (Sterne & Egger, 2001). A symmetrical, inverted funnel indicates absence of bias. In addition, funnel plot asymmetry was tested using Begg’s test, which examines the association between the effect estimates and their variances. Lastly, file drawer analysis was completed (Rosenthal, 1979). This 𝑜𝑜𝑏 calculates the number of studies averaging null results that would have to be added to nullify the summary effect (i.e. reduce the combined significance level (p-value) to a target alpha level (e.g. 0.05)). The meta-analysis was completed in R (R Core Team, 2018) using the metafor (Viechtbauer, 2010) and MetaForest packages (van Lissa, 2017). The full reproducible code is available online at: OSF.IO/S2JV4. 3 Results Key characteristics of each included study are shown in Table 2. The initial database search yielded 8,958 results, of which 2,336 were removed as duplicates. 6,622 titles and abstracts were screened and those deemed potentially relevant were retrieved as full texts. 76 studies were identified for full text screening. A further nine were retrieved from checking reference lists. In total, 33 studies met the inclusion criteria. For the flow diagram of this process, see Figure S1. Study Characteristics Studies came from 10 different countries, with most originating from Japan (n=13) or the US (n=5). Nine were published in Europe. The majority were randomised crossover studies (n=16) or non- randomised crossover studies (n=5). Seven studies used parallel groups, three had a factorial design and two were single-group crossover studies. Table 2. Study characteristics Reference Country Sample Design Intervention Setting(s) Depression Outcome Measure(s) Active engagement during exposure Berman et USA n=20 (8 men) Randomised 2.8 mile, 50-55 Natural environment: PANAS - Decreases in negative affect observed after both al., 2013 Mean age: 26 crossover min walk. Ann Arbor Arboretum the nature walk and the urban walk. MDD diagnosis Built environment: - No significant effect of location on negative downtown Ann Arbor. affect, but significant main effect of time. Bodin & Sweden n=12 (6 men) Randomised 60 min running Natural environment: NMS - The declines in anxiety/depression and anger Hartig, 2003 Male mean age: crossover route (max nature reserve. (anxiety/depression, subscales from pre-test to post-run were 39.7, SD 6.1 (4-period 2- 14km). Built environment: route anger subscales) significant. Female mean age: treatment) through Uppsala city - Change in NMS in the park environment was not 37.0, SD 7.0 significantly different to change in the urban environment. Gidlow et England n=38 (23 men) Randomised 30 min walk. Natural environment: BRUMS - Mood improved from baseline in all al., 2016 (UK) Mean age: 40.9, crossover country park within the environments, with a significant main effect of SD 17.6. (3-period, 3- city time on TMD, post-walk and 30 minutes after Healthy treatment) Built environment: quiet leaving the environment. residential streets - No significant main effect for environment and no Blue site: footpath significant environment*time interaction effect. besides a canal Hartig et al., USA n=34 (17 men) Randomised 40 min walk Natural environment: OHS; ZIPERS - The mean sadness score on ZIPERS was not 1991 Mean age: 20 parallel (sitting in passive Santiago Oaks Regional significantly different for the natural environment (Study 2) College students group site). Park, California group compared to the other groups. Built environment: Santa Ana, California Passive site: University of California campus Jia et al., China n=18 Randomised 90 min walk Natural environment: 65-item POMS - Depression was significantly reduced in the forest 2016 COPD patients parallel Repeated in the White Horse Mountain group between pre and post-exposure.. Mean age: 70.1 group morning and National Forest Park. - No significant difference seen in the urban group, afternoon for 3 Built environment: and score not significantly different between days. Hangzhou city groups. Johansson, Sweden n=20 2x2 factorial 40 min walk Natural environment: NMS - Statistically significant main effect of time for all Hartig, & Mean age: 23.3 design With/without Municipal park affect measures. Staats, 2011 Students friend as within- - Negative affect was not significantly modified by subject factor. environment or social context. Built environment: Street walk in mixed land use area Lee et al., Japan n=48 men Randomised 12-15 min walk. Natural environment: 30-item POMS - No significant change before and after forest and 2014 Mean age: 21.1 crossover forest urban walks on the POMS depression subscale. Mean BMI: 21.3 Built environment: urbanised area near forest Li et al., Japan n=19 men Single 2.6km, 80 min Natural environment: 65-item POMS - Significant decrease in the depression subscale 2016 Mean age: 51.2, group walk. Repeated in forest park after walking in the forest in the morning compared SD 8.8 crossover morning and Built environment: urban to before walking. High-normal or afternoon. Visited area in Nagano - No information on afternoon or Built hypertension other site one prefecture. environment; assumed results were not significant. week later. Mao, Cao, China n=34 Randomised 90min walk with Natural environment: 65-item POMS - The forest group had a significantly lower et al., 2012 Age range for parallel a 20 min break. White Horse Mountain depression subscale score between baseline and inclusion: 60-75 group Repeated in the National Forest Park. post-intervention. years. morning and Built environment: - No significant change in the urban group. Patients with afternoon for 7 location in Hangzhou city essential days. hypertension. Mao, Lan, et China n=20 men Randomised 90 min walk with Natural environment: 65-item POMS - Depression subscale score was significantly lower al., 2012 Mean age: 20.79, parallel a 10 min break. evergreen forest in than that of the urban group post-intervention. SD 0.54. group Repeated in the Hangzhou, China Students afternoon. Built environment: nearby urban area Mao et al., China n=33 Randomised 90 min walk, Natural environment: 65-item POMS - Significant decrease for the forest group in the 2017 Mean age: 72.2 parallel twice a day for forest site in Pan’an depression subscale compared with baseline score. CHF patients group four days. county - Post-intervention score for this subscale also Built environment: significantly lower than the post-intervention score downtown area of for the urban group. Hangzhou Perkins et USA n=26 (7 men) Parallel 20 min walk Natural environment: 65-item POMS - For all settings, change in depression score pre al., 2011 Age range 19-24 group wooded trail and post-intervention not significantly different. Students (3 groups) Built environments: mixed residential/business neighbourhood; parking lot. Roe & Scotland n=24 (11 in good Non- 60 min walk. Natural environment: MACL (hedonic - Significant positive change in mood following the Aspinall, (UK) mental health, 13 randomised Repeated in other Plean Country Park, tone, energy, stress). rural walk; no significant change following the 2011 in poor mental crossover setting one week Stirlingshire urban walk for the good health group. (Study 2) health) later. Built environment: - Significant positive change in mood for both the Poor health: Stirling town centre urban and rural walks for the poor health group. clinically diagnosed mental health problem Shin, Shin, Korea n=60 (35 men, 25 Randomised 4.5km, 50-55 min Natural environment: 65-item POMS - All POMS subscales, including depression, Yeoun, & women) crossover walk. Natural forested park significantly improved following exposure to the Kim, 2011 Mean age 23.27 Repeated in other Built environment: forest. Students setting one week downtown Cheongj - Depression subscale score increased following later. the urban exposure, but this was not significant. Song et al., Japan n=13 men Non- 15 min walk. Natural environment: 30-item POMS - No significant difference between settings was 2013 Mean age 22.5, SD randomised Rested for 20 Kashiwanoha Park in observed in the POMS depression subscale. 3.1 crossover mins then Chiba, Japan Students repeated in other Built environment: city setting. area around the park Song et al., Japan n=17 men Non- 15 min walk. Natural environment: 30-item POMS - No significant difference between settings was 2014 Mean age: 1.2, SD randomised Repeated in other Kashiwanoha park in observed in the POMS depression subscale. 1.7 crossover setting one hour Kashiwa City Students later. Built environment: city area around the park Song, Ikei, Japan n=19 men Randomised 17 min walk. Natural environment: 30-item POMS - Depression subscale of POMS was significantly Kobayashi, Mean age: 58, SD crossover Repeated in other Akasawa natural lower after walking in forest than walking in the et al., 2015 10.6 setting the next recreation forest urban area. High-normal or day. Built environment: Ina hypertensive blood city pressure Song, Ikei, Japan n=23 men Non- 15 min walk. Natural environment: 30-item POMS - No significant difference between settings was Igarashi, Mean age 22.3, SD randomised Rested for 20 min Kashiwa-no-ha Park observed in the POMS depression subscale. Takagaki, & 1.2 crossover then repeated in Built environment: City Miyazaki, Students other setting. area around the park Stigsdotter, Denmark n=51 women Non- 15 min walk. Natural environment: 65-item POMS - Depression scores reduced after both the forest Corazon, Age range 20-36 randomised Other setting Octavia health forest and urban walk, however the change was Sidenius, University students crossover visited within 2 Built environment: significant only for the urban walk. Kristiansen, weeks. historic downtown area - The depression score was significantly higher in & Grahn, of Copenhagen the urban environment pre-exposure compared to 2017 the forest. Triguero- Catalonia n=26 Randomised Participants asked Natural environment: 29-item POMS - TMD was significantly lower in the green site Mas, (Spain) Eligible if Mental case- to “spend time” in Collserola Natural Park compared with the built environment. Gidlow, et Health Inventory crossover environment. Built environment: - TMD was significantly lower in the blue site al., 2017 score in the lower 30 min + 180 min. Eixample neighbourhood compared to the green site. th 50 percentile in Barcelona Blue site: Castelldefels beach Passive engagement during exposure Bielinis, Poland n=62 (36 men) Randomised 15 min standing Natural environment: PANAS, 65-item - Depression score on POMS after exposure to Takayama, Mean age: 21.45, parallel in environment. urban forest (deciduous, POMS forest environment significantly lower to score Boiko, SD 0.18 group Repeated in other broad-leaved) after exposure to the urban area. Omelan, & Students setting in the Built environment: - Negative affect was significantly higher Bielinis, afternoon. Olsztyn city, Poland following exposure to the urban area compared to 2017 after exposure to the forest environment. Hartig et al., USA n=101 2x2 factorial Participants had 5 Natural environment: PANAS; ZIPERS - Those in the natural environment reported lower 1999 Mean age: 20.6 design mins to draw Botanical gardens on sadness in ZIPERS and negative affect in PANAS (Study 3) Students environment, then Berkeley campus than those in the urban environment. 10 mins sitting. Urban: Busy traffic - The difference between groups was not intersection significant. Igarashi et Japan n=17 women Randomised 10 mins sitting. Natural environment: 30-item POMS - Score on depression subscale on POMS al., 2015 Mean age: 46.1 crossover Repeated in other kiwifruit orchard significantly lower following orchard visit Mean BMI: 21.4 environment Built environment: compared to the building site visit. immediately after building site first Joung et al., Korea n=8 Single 15 mins sitting Natural environment: POMS - POMS depression subscale was not significantly 2015 Mean age: 22.0. group local forest different between environments after exposure. Students crossover Built environment: Daejeon city (conducted on rooftop) Lee et al., Japan n=12 men Randomised 15 mins sitting Natural environment: 30-item POMS - Depression subscale score was higher in the urban 2011 Mean age: 21 crossover Repeated in other forest in Tsurui village, area than the forest following exposure. No Mean BMI: 22.5 setting the next Hokkaido significant differences were observed. Students day Built environment: commercial area of Kushiro town, Hokkaido Lee et al., Japan n=12 men Randomised 15 mins sitting in Natural environment: POMS - Depression subscale score was significantly lower 2015 Mean age: 22.3 crossover environment paddy field in Ukiha city for the rural environment compared to the urban Students Repeated in other in southern Japan environment post-exposure. setting the next Built environment: day Hakata railway station Sonntag- Sweden n=20 women Randomised 40 mins sitting in Multiple natural Adapted mood - Significant differences between different Öström et Mean age: 41.6, crossover environment environments: Forest by questionnaire based environments found for all mood scales, except al., 2014 SD 7.3 (4-period, 4- lake, rocky outcrop, on POMS and exhausted-alert. Mean level of treatment) spruce forest ZIPERS - Participants rated higher on the scales relaxed, burnout: 5.7 (7- Built environment: Umea happy, harmonious, peaceful and clearheaded in all point Burnout city forest environments compared to the city. Questionnaire) Tsunetsugu Japan n=46 males Non- 15 mins sitting in Natural environment: 65-item POMS - No significant change in the depression subscale et al., 2013 Mean age: 21.1, randomised environment four forests in central and score of POMS. SD 1.1 crossover Repeated next day Western Japan Students in other Built environment: four environment. urban areas in central and Western Japan Combination of active and passive engagement during exposure Hartig, USA n=112 (56 men) 2x2 factorial Cognitive task Natural environment: OHS; ZIPERS - Those that walked in the nature reserve Evans, Students design prior to walk as Audubon Society’s Starr experienced more positive emotion than those Jamner, Mean age 20.8 SD within-subject Ranch Sanctuary, walking in the urban environment, within the no- Davis, & 3.7 factor. California task condition. The main effect of environment Gärling, Healthy 10 mins sitting, Built environment: among the task subjects was not significant. 2003 then 50 mins Orange city, California walk. Park, Japan n=280 Randomised 14 mins sitting, 16 12 forest sites and 12 POMS - The POMS depression subscale score Tsunetsugu, Mean age: 21.7, crossover mins walk. urban areas across Japan significantly decreased following the viewing of Kasetani, SD 1.5 Repeated in other the city area and improved in the forest area. Kagawa, & Students environment the - When walking, the change in the average POMS Miyazaki, next day. depression subscale score was also significantly 2010 different between the forest and city areas. Park et al., Japan n=168 men Randomised 15 mins sitting in 14 forest sites and 14 POMS - After performing both activities, TMD scores 2011 Mean age: 20.4, crossover the morning, 15 urban areas across Japan were significantly lower for the forest areas than SD 4.1 mins walk in the for the urban areas. However, no significant Students afternoon. differences were observed for the depression Repeated in other subscale. environment the next day Takayama Japan n=45 men Randomised 15 mins walk in 4 forest sites and 4 urban PANAS; POMS - Depression score was significantly lower after the et al., 2014 crossover the morning, 15 areas across Japan. experiment in the forest environment; no mins viewing in significant change to depression in the urban the afternoon. environment. Repeated in other - No significant difference between environments environment the for the depression subscale post-exposure next day - No statistical difference between before the experiment and after viewing either in negative affect or positive affect. Tyrväinen Finland n=77 (6 men) Randomised 15 minutes Multiple green sites: PANAS - People had fewer negative emotions in the forest et al., 2014 Mean age: 47.6 crossover viewing, 30 mins Alppipuisto (urban park), compared to the city walk (approx. and Keskuspiusto (large - Interaction between environment and time was 2km) urban woodland) not significant. At least one week Built environment: between each visit Helsinki city centre Note. MACL: Mood Adjective Checklist; NMS: Negative Mood Scale; OHS: Overall Happiness Scale; PANAS: Positive and Negative Affect Schedule; POMS: Profile of Mood States; TMD: Total Mood Disturbance, ZIPERS: Zuckerman Inventory of Personal Reactions Participants Sample sizes ranged from 8 participants (Joung et al., 2015) to 280 (Park et al., 2010). On the whole samples were small with 76% of studies (n=25) including less than 50 participants. Participants were typically young, with just over half of studies (n=18) recruiting college or university students. Some studies specified a clinical population. This included persons with: major depressive disorder (Berman et al., 2013), high-normal or hypertension (Li et al., 2016; Mao et al., 2012a; Song et al., 2015), congestive heart failure (CHF) (Mao et al., 2017), chronic obstructive pulmonary disease (COPD) (Jia et al., 2016), a mental health problem (Roe & Aspinall, 2011), a high level of burnout (Sonntag-Ostrom et al., 2014), and a poor Mental Health Inventory score (Triguero-Mas et al., 2017). Intervention Exposure time ranged from 10 minutes to 90 minutes, with 15 minutes being the most common (n=11). Some studies had multiple exposures within one day (Li et al., 2016; Park et al., 2011; Takayama et al., 2014; Triguero-Mas et al., 2017). Three studies from the same researcher had a considerably longer exposure time whereby participants completed a walk in the morning and afternoon for a period of 2 days (Mao et al., 2012b), 4 days (Mao et al., 2017) and 7 days (Mao et al., 2012a). For crossover studies, most often the second environment would be visited the following day (n=8), Nine studies specified a length of time ranging from “at least five days apart” to “within two weeks”. One study only indicated that visits were undertaken within the same season (Sonntag-Ostrom et al., 2014). In contrast, in five studies participants visited both environments on the same day, with one study giving participants only one minute to turnaround to face the other environment and three minutes to rest before measurements began again (Igarashi et al., 2015). Sixteen studies had participants actively engage with the environment – most asked participants to walk, and one asked participants to complete a run (Bodin and Hartig, 2003). One study allowed participants to choose what to do, only asking them to “spend time” in the environment (Triguero-Mas et al., 2017). This was coded as active engagement as it was assumed participants would move somewhat within the exposure area. A second group of studies had participants passively engage: seven asked participants to sit and view the environment, and one had participants stand due to the cold weather (Bielenis et al., 2018). Five used a combined approach whereby participants walked in and then viewed the environment or vice versa (Hartig et al., 2003; Park et al., 2010; Park et al., 2011; Takayama et al., 2014; Tyrvainen et al., 2014). Outcome Measures The most frequently used mood measure was the Profile of Mood States (POMS) (n=22). Also used was the Positive and Negative Affect Scale (PANAS) (n=5) and Zuckerman Inventory of Personal Reactions (ZIPERS) (n=3). The Overall Happiness Scale (OHS) and Negative Mood Scale (NMS) were used twice each, and the Mood Adjective Checklist (MACL) and Brunel Mood Scale (BRUMS) were all used once. One study used a bespoke questionnaire that was based on the POMS and ZIPERS (Sonntag-Ostrom et al., 2014). Most studies (n=24) took mood measurements pre- and post-exposure (including studies that had more than one exposure per environment). Nine studies measured mood at post-exposure only. In addition to pre-post measurement, five studies took a baseline measurement before traveling to the exposure environment (Gidlow et al., 2016; Lee et al., 2011, Lee et al., 2015; Park et al., 2010; Triguero-Mas et al., 2017). Furthermore, two studies took a second post-exposure measurement: Gidlow et al. (2016) measured mood immediately after the exposure had finished, and then again 30 minutes later, while Triguero-Mas et al. (2017) completed a final measurement upon return to the lab. Two studies captured mood during the exposure period (Hartig et al. 2003; Li et al., 2016). Setting Most studies used forests as their natural environment (n=16), followed by urban or country parks (n=11). Four used natural environments characterised by their biodiversity (nature reserve, botanical garden) and two used more agricultural settings (kiwifruit orchard: Igarashi et al., 2015, paddy field: Lee et al., 2015). For the comparative built environment, most studies described a location in a downtown, urban area (n=27). One study indicated the location was in an urban area, but participants viewed the area from a rooftop (Joung et al. 2015). Two studies used a residential street (Gidlow et al., 2016; Triguero-Mas et al., 2017), one used a building site (Igarashi et al., 2015), and one a railway station (Lee et al., 2015). Some studies had additional environments that were not explored in this review: a canal path (Gidlow et al., 2016), a beach (Triguero-Mas et al., 2017), forest with rocky outcrop, forest by a lake (Sonntag- Ostrom et al., 2014) and a lab setting (Hartig et al., 1991). Risk of Bias 25 randomised studies were evaluated using the Cochrane Risk of Bias 1.0 tool (Higgins et al., 2011); 8 non-randomised studies were evaluated using the ROBINS-I (Sterne et al., 2016) (see Figure S2 and Table S2 respectively). Concerning the randomised studies, two studies described their method of randomisation and therefore were assigned a low risk of bias in this domain (Gidlow et al., 2016; Sonntag-Ostrom et al., 2014). One study was given high risk of bias because participants were assigned to an exposure group based on participant availability (Triguero-Mas et al., 2017). Remaining studies did not describe their method of randomisation, and so were rated as unclear. A number of important confounders were identified in the non-randomised studies, meaning six (of 8 total) were marked as having moderate or serious risk of bias. Many confounders were possible, but in completing the assessment particular attention was paid to: the weather, food, alcohol and caffeine consumption; social interaction with other participants or researchers; the environment participants were exposed to immediately before measurements started; and the length of time between the experimental and control environment exposures (if applicable). For example, Li et al. (2016) prohibited alcohol, caffeine and smoking during the study period, and participants were not allowed to speak to each other during their walk in the exposure environment. However, they state that the weather was sunny for the built environment exposure, and rainy and cloudy for the natural environment exposure. It was also not clear how the participants travelled to the exposure environments, therefore the type of environment they were exposed to prior to measurement and possible social interactions were not known. For these reasons, the study was marked as having serious risk of bias. Two studies gave too little information on the confounders listed to make an informed decision and were marked as ‘no information’ (Joung et al., 2015; Tsunetsugu et al., 2013). All non-randomised studies were judged as low risk for selection bias because the selection of participants was not related to the intervention or the outcome. Due to the nature of the interventions it was judged that there was also little risk of misclassification of intervention and control sites. No study was deemed to have ‘deviated from intended intervention’. In terms of blinding, all studies were judged as highly biased: blinding is impossible due to the nature of the studies. Some studies attempted to minimise bias by not informing participants which environment would be visited first (Gidlow et al., 2016; Sonntag-Öström et al., 2014). Nevertheless, the outcome was subjective, and participants were likely to be aware of the hypothesis being tested. It was judged likely that attrition was related to the outcome unless otherwise stated. For this reason, 12 of the randomised studies had a high risk of attrition bias. Studies that did not explain different numbers of participants reported in the methods and results were assigned an unclear rating (n=5). The remaining randomised studies (n=8) reported no drop outs, and therefore had a low risk of bias. Within the non-randomised studies, Roe and Aspinall (2011) noted that attrition was concentrated in their ‘poor health’ group only, therefore was rated as having serious risk of bias. All other studies had complete data or the proportion missing was limited. Two non-randomised studies did not report full results, resulting in a serious risk of reporting bias. All other studies received an unclear (using ROB) or moderate (using ROBINS-I) risk of bias since full data was reported but did not have associated study protocols. Narrative data synthesis Active engagement interventions (n=20) Eleven studies reported a significant decrease in depression pre and post-exposure to the natural environment. For example, Mao et al. (2012a) and (2017) had participants walk 90 minutes twice a day for 7 and 4 days respectively, both reporting a significant decrease in depressive mood in the forest environment compared to the pre-exposure score. Shin et al. (2011) had participants walk in a forest for 50-55 minutes, and in a city the following week. All POMS subscales, including depression, were found to significantly improve following the forest exposure. However, six of the eleven studies were not able to demonstrate that the change in mood was significantly different to that observed in the built environment (Berman et al. 2013; Bodin and Hartig, 2003; Gidlow et al., 2016; Jia et al., 2016; Johansson et al., 2011; Roe and Aspinall, 2011). Four studies showed no significant change in mood pre and post exposure to a natural environment (Stigsdotter et al. 2017; Hartig et al., 1991; Lee et al., 2014; Perkins et al., 2011). Li et al. (2016) found a significant decrease in the depression subscale after walking in the forest compared to baseline, but no results from the built environment are presented. Five studies assessed mood at post-exposure only (Mao et al., 2012b; Song et al., 2013; Song et al., 2014; Song et al., 2015a; Song et al., 2012b). Two reported that the POMS depression subscale score was significantly lower following the forest visit than the built environment visit (Mao et al., 2012b; Song et al., 2015a); the remaining three studies found no significant difference between environments. Passive engagement interventions (n=8) Four studies compared depressive mood pre- and post-exposure, with two finding a significant reduction. Bielinis et al. (2017) examined change in mood following winter forest bathing in young students. They were asked to stand in a forest and built environment for 15 minutes. The depression score was significantly lower after exposure to the forest compared to the built environment. Lee et al. (2015) also had an exposure period of 15 minutes. They report that depression was significantly lower in the rural environment post-exposure. There was no significant change in mood in the other two studies (Lee et al., 2011; Tsunetsugu et al., 2013). Four studies measured post-exposure only, also with two reporting significant results. Igarashi et al. (2015) asked women to sit for ten minutes in an orchard, and then in a building site. Depression was rated significantly lower after sitting in the orchard than the building site. Sonntag-Ostrom et al. (2014) sent 20 female patients with exhaustion disorder to four different environments in Sweden, testing for mood, attention and physiological response. Environments were three forest environments (spruce forest, forest with a lake, forest with rocky outcrop) and a built environment. Patients reported feeling significantly more happy, relaxed, harmonious, peaceful and clearheaded in all forest environments compared to the city environment. The remaining two studies (Hartig et al., 1999; Joung et al., 2015) did not find a significant difference in mood following exposure. Joung et al. (2015) note this may be explained by the fact that participants sat on a rooftop to observe an urban area, therefore preventing full immersion of the participant in the environment. Combination of active and passive engagement interventions (n=5) Park et al. (2010) and Takayama et al. (2014) both reported a significant reduction in the POMS depression subscale after a 15 minute walk and 15 minute viewing session in a forest. Park et al. (2011) followed the same procedure of 15 minutes walking and viewing in 14 forest and built environment sites across Japan, but no significant change in depressive mood was observed. Tyrvainen et al. (2014) compared results across three environments: an urban park, urban woodland, and city centre. People experienced fewer negative emotions in the woodland compared to the park and city centre, but there was no interaction between place and time. Similar results were found in Johansson et al. (2011), who compared participants walking in a park and down a street, and with and without a friend accompanying them. There was a significant main effect of time, but change in negative affect was not modified by environment or social context. Finally, Hartig et al. (2003) asked half of the participants to complete a cognitively demanding task. Those who did not complete the task experienced more positive emotion following the natural environment walk than the built environment walk, however for those who completed the task, there was no significant main effect of environment. Meta-analysis Descriptive statistics Observed effect sizes ranged from -2.30 to 0.84. The unweighted mean effect size was 𝑀 = −0.29, 𝑆𝐷 = 0.60, which can be interpreted as a small effect. Six studies reported two effect sizes (i.e. two outcomes). Because of this, a three-level meta-analysis was first used to estimate the amount of within- study- and between-studies variance (Van den Noortgate, López-López, Marín-Martínez, & Sánchez- Meca, 2014). As indicated in Table 3, the within-study variance component did not differ significantly from zero, 𝜎 < 0.01, 95% CI [0, 0.28]. The between-studies variance component, on the other hand, was significant, 𝜎 = 0.31, 95% CI [0.09, 0.55]. Thus, the variation in observed effect sizes was primarily due to differences between studies. As the within-study variance component was near-zero, and the Akaike information criterion (AIC) and Bayesian Information Criterion (BIC) for a model with within- studies variance constrained to zero were lowest out of all models compared, there was no advantage to the multilevel approach. Therefore a random-effects meta-analysis, which only includes a between- studies variance component, was conducted. Table 3. Comparing the fit of different multi-level models df AIC BIC ll LRT p Full three-level model 3 79.86 84.83 -36.92 Between-studies variance constrained 2 93.43 96.76 -44.72 15.59 0.000 Within-studies variance constrained 2 77.84 81.17 -36.92 0.00 1.000 Both variance components constrained 1 256.27 257.94 -127.14 180.43 0.000 Summary effect size The summary effect from random-effect meta-analysis was significantly different from zero, 𝛾 = −0.30,𝑝 < 0.01, 95% CI[−0.50,−0.10]. The random effect was also significant, indicating that there was residual heterogeneity between studies, 𝜏 = 0.31, 𝑆𝐸 = 0.09, 𝑄 (39) = 277.97, 𝑝 < 0.01. 𝑑𝑖𝑠𝑒𝑟 This is reflective of the diversity of studies included in the meta-analysis. See Figure 1 for a forest plot of the included studies. Figure 1. Forest plot of study effect sizes Publication bias Visual inspection of the funnel plot in Figure 2 was inconclusive with regard to publication bias. There was a lack of studies in both lower corners of the funnel plot, indicating that most small-sample studies reported effects with values close to the average weighted effect size. There were also some relatively high-powered studies with large, negative effects. However, Beggs’ test of funnel asymmetry was non- significant (𝑧 = 0.93,𝑝 = 0.35). File drawer analysis indicated that 846 unpublished or unretrieved studies averaging null results would have to be added to render the average unweighted effect size non- significant. Thus, the extent of publication bias was hard to ascertain. Figure 2. Funnel plot to assess potential publication bias Moderation analysis To investigate the source of heterogeneity, a random-effects MetaForest analysis was conducted with 10,000 iterations and replicated 100 times to ensure the reliability of findings. The replicated variable importance metrics can be seen in Figure S4. All variables that reduced the predictive performance of the model were dropped and the remaining eight carried forward to optimize the model. The estimated 2 2 predictive performance in new data was positive; cross-validated 𝑅 = 0.42, out-of-bag 𝑅 = 𝑐𝑣 0.05.The relative variable importance of the moderators in the final model is shown in Figure S5. The model identifies the proportion of females in the sample, the type of natural and built environment, the type of effect size, the time between natural and built environment visits, the country of study origin, the gender mix of the sample, and whether or not a baseline measurement was taken to be the most important moderators of the effect size from the 20 that were entered. Partial dependency plots (Figure S6) were produced to examine the influence of each moderator on the effect size, while averaging over all other moderators. The model predicts that for a sample with a lower proportion of women, the effect size is larger. The effect size was also larger for agricultural, biodiverse and forest environments, compared to a park environment. Categories within other moderators showed similar relationships with the effect size. The eight most important moderators were entered into a meta- regression (Table 4), however, none were significant in a linear model. 𝑜𝑜𝑏 Table 4. Meta regression model with most important moderators Variable Estimate SE Z p CI Intercept -0.09 0.32 -0.28 0.78 [-0.72, 0.54] Female 0.01 0.01 0.60 0.55 [-0.02, 0.04] Natural environment: ABF vs Park -0.14 0.20 -0.70 0.48 [-0.53, 0.25] Natural environment: A vs BF -0.35 0.49 -0.71 0.48 [-1.31, 0.61] Natural environment: B vs F 0.07 0.31 0.22 0.83 [-0.54, 0.68] Type of ES: SMCR vs SMD 0.15 0.18 0.87 0.38 [-0.19, 0.50] Time between environments: Linear 0.97 0.62 1.57 0.12 [-0.24, 2.19] Time between environments: Quad. 0.58 0.50 1.17 0.24 [-0.39, 1.56] Time between environments: Cubic -1.02 0.60 -1.71 0.09 [-2.19, 0.15] Country: Asia vs USEurope 0.14 0.16 0.90 0.37 [-0.17, 0.45] Country: US vs Europe -0.08 0.27 -0.29 0.77 [-0.62, 0.46] Sex: Mixed vs FemaleMale -0.79 1.11 -0.72 0.47 [-2.96, 1.37] Sex: Male vs Female -2.16 2.11 -1.03 0.30 [-6.30, 1.97] Baseline measurement: No vs Yes 0.13 0.22 0.61 0.54 [-0.29, 0.55] Built environment: DowntownOther vs 0.20 0.27 0.72 0.47 [-0.34, 0.73] Residential Built environment: Downtown vs Other -0.40 0.51 -0.78 0.43 [-1.40, 0.60] Quality of evidence A summary of findings table is presented in Figure S3. For randomised studies (n=25), initial confidence is high. However, the studies were downgraded due to serious risk of bias, inconsistency between studies, and plausible confounding. It was judged that because all studies had received a serious risk of bias in the individual study assessments, the overall body evidence would equally be deemed to have a serious risk of bias. Concerning the inconsistency domain, during the meta-analysis it was found that there was significant residual heterogeneity, therefore this was also marked with serious concerns. Lastly, a number of confounding variables were identified during the study-level bias assessments, and so this judgement was also applied to the body of evidence as a whole. In line with GRADE guidelines, the non-randomised studies (n=8) started as low quality due to residual confounding. In addition to the aforementioned judgements, these studies were further downgraded for imprecision due to small sample sizes and wide confidence intervals. Publication bias, an overall large effect and a dose response gradient were not identified. Overall, the randomised studies were deemed to be of low quality and the non-randomised studies of very low quality. 4 Discussion In this review and meta-analysis, 33 studies that investigated the effect of direct, short-term exposure to the natural environment on depressive mood were synthesised. Effect sizes ranged from -2.30 to 0.84, with an unweighted mean effect size of -0.29. However, risk of bias and quality assessments determined the current evidence is highly biased and of poor quality. Confidence in our conclusions is therefore limited, and the summary effect must be interpreted with caution. The meta-analysis also revealed significant residual heterogeneity between studies, which remains largely unexplained following moderator analysis. Nevertheless, this review is in line with a previous meta-analyses of the effect of natural environment exposure on mood (r=-0.12, McMahan & Estes, 2015). The results also complement previous systematic reviews on the mental health benefits of, for example, forest bathing, horticultural activities, and green exercise (Bowler, Buyung-Ali, Knight, & Pullin, 2010; Hansen et al., 2017; Lee et al., 2017; Oh et al., 2017). In particular the findings converge with previous research into the presence of a dose- response relationship, which has demonstrated a boost to mood following a short period of exposure to nature. The most common exposure time identified in this review was 15 minutes - the experiments of Barton and Pretty 2010; and Shanahan et al. 2015, 2016) showed changes in mood following 5 minutes and 30 minutes of exposure time respectively. Overall, the review contributes to a growing evidence base concerning the mental health benefits of exposure to nature. During moderator analysis using the MetaForest approach, it was found that proportion of female participants, type of natural environment and built environment, time between environments, baseline measurement, region of study origin, gender mix of the sample and type of effect size were important moderators of the effect size. This analysis draws some similarities with the results of a previous meta- analysis that found that type of emotion assessment, type of exposure to nature, location of study, and mean age significantly moderated the effect of nature on positive mood (McMahan & Estes, 2015). In addition the current analysis finds several between-study moderators to be relevant, which is reflective of the diversity of included studies. On the other hand, none were significant when entered in a meta- regression. The lack of significance might be explained by the potential for bias in the included studies. A number of confounders were identified which may have influenced the results. It is not known to what extent carryover effects, whereby the effect of one environment might be ‘carried over’ to the next, might contribute to results: a wide range in the duration of time between environments was found. Another issue more generally relates to the issue of blinding participants and outcome assessors. It is essentially impossible to blind persons involved in interventions of this kind, since awareness of the environment is necessary. Van den Berg (2017) explains that these issues represent a key challenge in encouraging ‘green prescriptions’ (greening a person’s environment, or taking them to a green environment, in order to promote health). Health professionals are inclined toward the results of randomised controlled trials, and often this approach is not appropriate for a nature-based intervention. There were attempts to reduce bias in some studies. For example, two studies did not give prior warning of the order in which environments would be visited (Gidlow et al., 2016; Sonntag-Öström et al., 2014). These two studies also reported a clear process of randomisation. No study protocols could be found for studies included in this review, however a recent study of a park prescription program has done this (Razani et al., 2018). It is recommended that future research in the area take steps to reduce bias and improve quality where possible, in order to build a strong clinical evidence base. This will work to persuade policymakers and health professionals of the mental health benefits of exposure to nature. Strengths and limitations This review provides the most comprehensive and up-to-date findings on the effect of short-term exposure to the natural environment on depressive mood. Its key strengths are its broad range of included studies, and a fully reproducible and transparent meta-analysis. However, this review also had some limitations. First, it was limited to English articles only. This prevented articles written in other languages from being included, however, a previous review that included relevant articles written in Korean did not find dissimilar results to this review (Lee et al., 2017). Second, it was not possible to retrieve data for two studies to enter into the meta-analysis. Both studies reported a reduction in depressive mood following exposure to the natural environment (Li et al., 2016; Perkins et al., 2011). Lastly, this review was concerned with short-term exposure only and does not address long-term effects of repeated exposure. It is assumed that repeated exposures would be cumulatively beneficial, and indeed a recent review found that long-term exposure to increased green and blue space in the residential environment is associated with improved mental health (Gascon et al., 2015). Future research Three suggestions for further research are made. First, future meta-analysis would benefit from improved descriptions and reporting of studies. For example, studies should provide an objective description of the experimental and control environments. This might be achieved by measuring the Normalized Difference Vegetation Index of the area, making use of street view imagery, or calculating percent tree canopy. Next, a detailed procedural description is required to fully understand the environmental context within which participants are placed. Moreover, it is recommended that appropriate guidelines are followed during reporting, for example, CONSORT (Schulz, Altman, & Moher, 2010) for randomised trials and TREND (Des Jarlais, Lyles, & Crepaz, 2004) for non- randomised studies. This ensures studies are fully described in a standardised manner. Second, the MetaForest analysis revealed eight moderators that were associated with the effect size. In particular, the type of natural environment and proportion of females in the sample were the two most important moderators. The partial dependent plots showed that a larger effect size was associated with a lower proportion of women, and also in agricultural, biodiverse and forest environments, compared to the park environment. Further, the majority of studies also had young, usually male, university students as their participants. This reduces generalisability to other populations. It is therefore suggested that future research continues to explore the potential moderating role of type of environment and type of population group. This is important to understand in order to develop effective interventions to promote mood. Lastly, increasing research is applying technology such as Global Positioning System, wearables, and ecological momentary assessment to investigate mental state over time and space (Bakolis et al., 2018; Birenboim, Dijst, Scheepers, Poelman, & Helbich, 2019; Chaix, 2018; Helbich, 2018). This represents the next step in this field of research whereby pre- and post-measures can be reformulated into a more dynamic approach. This removes the need for experimental procedure as participants can be followed in their daily life, and the effects of varying exposure duration and potential accumulation effects and long-term mental health benefits might be considered. 5 Conclusions This review and meta-analysis finds a reduction in depressive mood following short-term exposure to the natural environment, however, studies were highly biased and of low quality. It is therefore unclear whether these findings would be replicated in higher quality studies. No significant moderators of the effect size were identified. More rigorous studies are required to improve our understanding of the relationship between the natural environment and mood. Acknowledgements We wish to thank all authors that provided additional data and/or clarification to complete the meta- analysis, namely: Christopher Gidlow, Terry Hartig, Gemma Hurst, Juyoung Lee, Genxiang Mao, Yoshifumi Miyazaki, Jenny Roe, Russ Searight, Chorong Song, Norimasa Takayama, Margarita Triguero-Mas. Funding This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 714993). Declaration of interest None. Conflict of interest The authors have no conflict of interest to declare. References Alcock, I., White, M. P., Wheeler, B. W., Fleming, L. E., & Depledge, M. H. (2014). Longitudinal effects on mental health of moving to greener and less green urban areas. Environmental Science and Technology, 48(2), 1247–1255. https://doi.org/10.1021/es403688w Astell-Burt, T., Mitchell, R., & Hartig, T. (2014). The association between green space and mental health varies across the lifecourse. A longitudinal study. Journal of Epidemiology and Community Health, 68(6), 578–583. https://doi.org/10.1136/jech-2013-203767 Bakolis, I., Hammoud, R., Smythe, M., Gibbons, J., Davidson, N., Tognin, S., & Mechelli, A. (2018). Urban Mind: Using Smartphone Technologies to Investigate the Impact of Nature on Mental Well-Being in Real Time. BioScience, 68(2), 134–145. https://doi.org/10.1093/biosci/bix149 Barton, J., & Pretty, J. (2010). What is the best dose of nature and green exercise for improving mental health- A multi-study analysis. Environmental Science and Technology, 44(10), 3947– 3955. https://doi.org/10.1021/es903183r Berman, M. G., Kross, E., Krpan, K. M., Askren, M. K., Burson, A., Deldin, P. J., … Jonides, J. (2013). Interacting with Nature Improves Cognition and Affect for Individuals with Depression. Journal of Affective Disorders, 140(3), 300–305. https://doi.org/10.1016/j.jad.2012.03.012.Interacting Beyer, K. M. M., Kaltenbach, A., Szabo, A., Bogar, S., & Nieto, F. J. (2014). Exposure to Neighborhood Green Space and Mental Health : Evidence from the Survey of the Health of Wisconsin. International Journal of Environmental Research & Public Health, 11, 3453–3472. https://doi.org/10.3390/ijerph110303453 Bielinis, E., Takayama, N., Boiko, S., Omelan, A., & Bielinis, L. (2017). The effect of winter forest bathing on psychological relaxation of young Polish adults. Urban Forestry and Urban Greening, 29, 276–283. https://doi.org/10.1016/j.ufug.2017.12.006 Birenboim, A., Dijst, M., Scheepers, F., Poelman, M., & Helbich, M. (2019). Wearables and location tracking technologies for mental-state sensing in outdoor environments. Professional Geographer. Bodin, M., & Hartig, T. (2003). Does the outdoor environment matter for psychological restoration gained through running? Psychology of Sport and Exercise, 4, 141–153. https://doi.org/10.1016/S1469-0292(01)00038-3 Bowler, D. E., Buyung-Ali, L. M., Knight, T. M., & Pullin, A. S. (2010). A systematic review of evidence for the added benefits to health of exposure to natural environments. BMC Public Health. https://doi.org/10.1186/1471-2458-10-456 Chaix, B. (2018). Mobile Sensing in Environmental Health and Neighborhood Research. Ssrn. https://doi.org/10.1146/annurev-publhealth-040617-013731 Des Jarlais, D. C., Lyles, C., & Crepaz, N. (2004). Improving the Reporting Quality of Nonrandomized Evaluations of Behavioral and Public Health Interventions: The TREND Statement. American Journal of Public Health, 94(3), 361–366. https://doi.org/10.2105/AJPH.94.3.361 Ferrari, A. J., Charlson, F. J., Norman, R. E., Patten, S. B., Freedman, G., Murray, C. J. L., … Whiteford, H. A. (2013). Burden of Depressive Disorders by Country, Sex, Age, and Year: Findings from the Global Burden of Disease Study 2010. PLoS Medicine, 10(11). https://doi.org/10.1371/journal.pmed.1001547 Fritz, C. O., Morris, P. E., & Richler, J. J. (2012). Effect size estimates: Current use, calculations, and interpretation. Journal of Experimental Psychology: General, 141(1), 2–18. https://doi.org/10.1037/a0024338 Gascon, M., Sánchez-Benavides, G., Dadvand, P., Martínez, D., Gramunt, N., Gotsens, X., … Nieuwenhuijsen, M. (2018). Long-term exposure to residential green and blue spaces and anxiety and depression in adults: A cross-sectional study. Environmental Research, 162(January), 231–239. https://doi.org/10.1016/j.envres.2018.01.012 Gascon, M., Triguero-Mas, M., Martínez, D., Dadvand, P., Forns, J., Plasència, A., & Nieuwenhuijsen, M. (2015). Mental Health Benefits of Long-Term Exposure to Residential Green and Blue Spaces: A Systematic Review. International Journal of Environmental Research and Public Health, 12(4), 4354–4379. https://doi.org/10.3390/ijerph120404354 Gascon, M., Zijlema, W., Vert, C., White, M. P., & Nieuwenhuijsen, M. J. (2017). Outdoor blue spaces, human health and well-being: A systematic review of quantitative studies. International Journal of Hygiene and Environmental Health, 220(8), 1207–1221. https://doi.org/10.1016/j.ijheh.2017.08.004 Gidlow, C. J., Jones, M. V., Hurst, G., Masterson, D., Clark-Carter, D., Tarvainen, M. P., … Nieuwenhuijsen, M. (2016). Where to put your best foot forward: Psycho-physiological responses to walking in natural and urban environments. Journal of Environmental Psychology, 45, 22–29. https://doi.org/10.1016/j.jenvp.2015.11.003 GRADE Working Group. (2004). Grading quality of evidence and strength of recommendations. BMJ (Clinical Research Ed.), 328(7454), 1490. https://doi.org/10.1136/bmj.328.7454.1490 Guyatt, G. H., Oxman, A., Vist, G., Kunz, R., Falck-Ytter, Y., Alonso-Coello, P., & Schünemann, H. J. (2008). GRADE an emerging consensus on rating quality of evidence and stren. British Medical Journal, 336(7650), 924–926. https://doi.org/10.1136/bmj.39489.470347.AD Hansen, M. M., Jones, R., & Tocchini, K. (2017). Shinrin-yoku (Forest bathing) and nature therapy: A state-of-the-art review. International Journal of Environmental Research and Public Health, 14(8). https://doi.org/10.3390/ijerph14080851 Hartig, T., Evans, G. W., Jamner, L. D., Davis, D. S., & Gärling, T. (2003). Tracking restoration in natural and urban field settings. Journal of Environmental Psychology, 23(2), 109–123. https://doi.org/10.1016/S0272-4944(02)00109-3 Hartig, T., Mang, M., & Evans, G. W. (1991). Restorative Effects of Natural Environment Experiences. Environment and Behavior, 23(1), 3–26. https://doi.org/10.1177/0013916591231001 Hartig, T., Mitchell, R., de Vries, S., & Frumkin, H. (2014). Nature and Health. Annual Review of Public Health, 35(1), 207–228. https://doi.org/10.1146/annurev-publhealth-032013-182443 Hartig, T., Nyberg, L., Nilsson, L.-G., & Garling, T. (1999). Testing for Mood Congruent Recall With Environmentally Induced Mood. Journal of Environmental Psychology, 19, 353–367. https://doi.org/10.1006/jevp.1999.0142 Hawton, K., Casañas I Comabella, C., Haw, C., & Saunders, K. (2013). Risk factors for suicide in individuals with depression: A systematic review. Journal of Affective Disorders. https://doi.org/10.1016/j.jad.2013.01.004 Hedges, L. V. (1981). Distribution Theory for Glass’s Estimator of Effect Size and Related Estimators. Journal of Educational Statistics, 6(2), 107–128. Helbich, M. (2018). Toward dynamic urban environmental exposure assessments in mental health research. Environmental Research, 161(October 2017), 129–135. https://doi.org/10.1016/j.envres.2017.11.006 Helbich, M., de Beurs, D., Kwan, M. P., O’Connor, R. C., & Groenewegen, P. P. (2018). Natural environments and suicide mortality in the Netherlands: a cross-sectional, ecological study. The Lancet Planetary Health, 2(3), e134–e139. https://doi.org/10.1016/S2542-5196(18)30033-0 Helbich, M., Yao, Y., Liu, Y., Zhang, J., Liu, P., & Wang, R. (2019). Using deep learning to examine street view green and blue spaces and their associations with geriatric depression in Beijing, China. Environment International, 126(January), in press. https://doi.org/10.1016/J.ENVINT.2019.02.013 Higgins, J. P. T., Altman, D. G., Gøtzsche, P. C., Jüni, P., Moher, D., Oxman, A. D., … Cochrane Statistical Methods Group. (2011). The Cochrane Collaboration’s tool for assessing risk of bias in randomised trials. BMJ, 343(7303), d5928. https://doi.org/10.1136/bmj.d5928 Houlden, V., Weich, S., de Albuquerque, J. P., Jarvis, S., & Rees, K. (2018). The relationship between greenspace and the mental wellbeing of adults: A systematic review. PLoS ONE (Vol. 13). https://doi.org/10.1371/journal.pone.0203000 Igarashi, M., Miwa, M., Ikei, H., Song, C., Takagaki, M., & Miyazaki, Y. (2015). Physiological and psychological effects of viewing a kiwifruit (Actinidia deliciosa ‘Hayward’) orchard landscape in summer in Japan. International Journal of Environmental Research and Public Health, 12, 6657–6668. https://doi.org/10.3390/ijerph120606657 Jia, B. B., Yang, Z. X., Mao, G. X., Lyu, Y. D., Wen, X. L., Xu, W. H., … Wang, G. F. (2016). Health Effect of Forest Bathing Trip on Elderly Patients with Chronic Obstructive Pulmonary Disease. Biomedical and Environmental Sciences, 29(3), 212–218. https://doi.org/10.3967/bes2016.026 Johansson, M., Hartig, T., & Staats, H. (2011). Psychological benefits of walking: Moderation by company and outdoor environment. Applied Psychology: Health and Well-Being, 3(3), 261–280. https://doi.org/10.1111/j.1758-0854.2011.01051.x Joung, D., Kim, G., Choi, Y., Lim, H., Park, S., Woo, J. M., & Park, B. J. (2015). The prefrontal cortex activity and psychological effects of viewing forest landscapes in Autumn season. International Journal of Environmental Research and Public Health, 12, 7235–7243. https://doi.org/10.3390/ijerph120707235 Kaplan, S. (1995). The restorative benefits of nature: Toward an integrative framework. Journal of Environmental Psychology, 15(3), 169–182. https://doi.org/10.1016/0272-4944(95)90001-2 Lee, I., Choi, H., Bang, K.-S., Kim, S., Song, M., & Lee, B. (2017). Effects of Forest Therapy on Depressive Symptoms among Adults: A Systematic Review. International Journal of Environmental Research and Public Health, 14(3), 321. https://doi.org/10.3390/ijerph14030321 Lee, J., Park, B.-J., Tsunetsugu, Y., Ohira, T., Kagawa, T., & Miyazaki, Y. (2011). Effect of forest bathing on physiological and psychological responses in young Japanese male subjects. Public Health, 125, 93–100. https://doi.org/10.1016/j.puhe.2010.09.005 Lee, J., Park, B. J., Ohira, T., Kagawa, T., & Miyazaki, Y. (2015). Acute Effects of Exposure to a Traditional Rural Environment on Urban Dwellers: A Crossover Field Study in Terraced Farmland. International Journal of Environmental Research and Public Health, 12, 1874–1893. https://doi.org/10.3390/ijerph120201874 Lee, J., Park, B. J., Tsunetsugu, Y., Ohira, T., Kagawa, T., & Miyazaki, Y. (2019). Corrigendum to “Effect of forest bathing on physiological and psychological responses in young Japanese male subjects” [Public Health 125 (2) (February 2011), 93– 100](S0033350610003203)(10.1016/j.puhe.2010.09.005). Public Health, 169, 201. https://doi.org/10.1016/j.puhe.2019.03.002 Lee, J., Tsunetsugu, Y., Takayama, N., Park, B. J., Li, Q., Song, C., … Miyazaki, Y. (2014). Influence of forest therapy on cardiovascular relaxation in young adults. Evidence-Based Complementary and Alternative Medicine, 2014. https://doi.org/10.1155/2014/834360 Li, Q., Otsuka, T., Kobayashi, M., Wakayama, Y., Inagaki, H., Katsumata, M., … Kagawa, T. (2016). Effects of forest environments on cardiovascular and metabolic parameters. Evidence-Based Complementary and Alternative Medicine. https://doi.org/http://dx.doi.org/10.1155/2016/2587381 Lim, G. Y., Tam, W. W., Lu, Y., Ho, C. S., Zhang, M. W., & Ho, R. C. (2018). Prevalence of Depression in the Community from 30 Countries between 1994 and 2014 /692/699/476/1414 /692/499 article. Scientific Reports, 8(1), 1–10. https://doi.org/10.1038/s41598-018-21243-x Mao, G. X., Cao, Y. B., Lan, X. G., He, Z. H., Chen, Z. M., Wang, Y. Z., … Yan, J. (2012). Therapeutic effect of forest bathing on human hypertension in the elderly. Journal of Cardiology, 60, 495–502. https://doi.org/10.1016/j.jjcc.2012.08.003 Mao, G. X., Cao, Y., Wang, B., Wang, S., Chen, Z., Wang, J., … Yan, J. (2017). The salutary influence of forest bathing on elderly patients with chronic heart failure. International Journal of Environmental Research and Public Health, 14(368). https://doi.org/10.3390/ijerph14040368 Mao, G. X., Lan, X. G., Cao, Y. B., Chen, Z. M., He, Z. H., Lv, Y. D., … Yan, J. (2012). Effects of short-term forest bathing on human health in a broad-leaved evergreen forest in Zhejiang Province, China. Biomedical and Environmental Sciences Journal, 25(3), 317–324. https://doi.org/10.3967/0895 Mceachan, R. R. C., Prady, S. L., Smith, G., Fairley, L., Cabieses, B., Gidlow, C., … Royal, B. (2015). The association between green space and depressive symptoms in pregnant women : moderating roles of socioeconomic status and physical activity, 1–7. https://doi.org/10.1136/jech-2015-205954 McMahan, E. A., & Estes, D. (2015). The effect of contact with natural environments on positive and negative affect: A meta-analysis. Journal of Positive Psychology, 10(6), 507–519. https://doi.org/10.1080/17439760.2014.994224 Min, K. bok, Kim, H. J., Kim, H. J., & Min, J. young. (2017). Parks and green areas and the risk for depression and suicidal indicators. International Journal of Public Health, 62(6), 647–656. https://doi.org/10.1007/s00038-017-0958-5 Moher, D., Liberati, A., Tetzlaff, J., Altman, D. G., & Group, T. P. (2009). Preferred Reporting Items for Systematic Reviews and Meta-Analyses : The PRISMA Statement. PLoS Medicine, 6(7), e1000097. https://doi.org/10.1371/journal.pmed.1000097 Morris, S. B., & DeShon, R. P. (2002). Combining effect size estimates in meta-analysis with repeated measures and independent-groups designs. Psychological Methods, 7(1), 105–125. https://doi.org/10.1037/1082-989X.7.1.105 Nieuwenhuijsen, M. J., Khreis, H., Triguero-Mas, M., Gascon, M., & Dadvand, P. (2017). Fifty Shades of Green: Pathway to Healthy Urban Living. Epidemiology, 28(1), 63–71. https://doi.org/10.1097/EDE.0000000000000549 Oh, B., Lee, K. J., Zaslawski, C., Yeung, A., Rosenthal, D., Larkey, L., & Back, M. (2017). Health and well-being benefits of spending time in forests: Systematic review. Environmental Health and Preventive Medicine, 22(1), 1–11. https://doi.org/10.1186/s12199-017-0677-9 Park, B. J., Furuya, K., Kasetani, T., Takayama, N., Kagawa, T., & Miyazaki, Y. (2011). Relationship between psychological responses and physical environments in forest settings. Landscape and Urban Planning, 102, 24–32. https://doi.org/10.1016/j.landurbplan.2011.03.005 Park, B. J., Tsunetsugu, Y., Kasetani, T., Kagawa, T., & Miyazaki, Y. (2010). The physiological effects of Shinrin-yoku (taking in the forest atmosphere or forest bathing): Evidence from field experiments in 24 forests across Japan. Environmental Health and Preventive Medicine, 15, 18– 26. https://doi.org/10.1007/s12199-009-0086-9 Perkins, S., Searight, H. R., & Ratwik, S. (2011). Walking in a Natural Winter Setting to Relieve Attention Fatigue: A Pilot Study. Psychology, 2(8), 777–780. https://doi.org/10.4236/psych.2011.28119 R Core Team. (2018). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Rautio, N., Filatova, S., Lehtiniemi, H., & Miettunen, J. (2017). Living environment and its relationship to depressive mood: A systematic review. International Journal of Social Psychiatry, 20764017744582. https://doi.org/https://dx.doi.org/10.1177/0020764017744582 Razani, N., Morshed, S., Kohn, M. A., Wells, N. M., Thompson, D., Alqassari, M., … Rutherford, G. W. (2018). Effect of park prescriptions with and without group visits to parks on stress reduction in low-income parents: SHINE randomized trial. PLoS ONE, 13(2), 1–17. https://doi.org/10.1371/journal.pone.0192921 Reklaitiene, R., Grazuleviciene, R., Dedele, A., Virviciute, D., Vensloviene, J., Tamosiunas, A., … Nieuwenhuijsen, M. J. (2014). The relationship of green space, depressive symptoms and perceived general health in urban population. Scandinavian Journal of Public Health, 42(7), 669–676. https://doi.org/10.1177/1403494814544494 Roberts, H., Hagedoorn, P., Kellar, I., & Helbich, M. (2018). The effect of short-term exposure to natural environments on depression: a systematic review. Retrieved from http://www.crd.york.ac.uk/PROSPERO/display_record.php?ID=CRD42018102612 Roe, J., & Aspinall, P. (2011). The restorative benefits of walking in urban and rural settings in adults with good and poor mental health. Health and Place, 17, 103–113. https://doi.org/10.1016/j.healthplace.2010.09.003 Rosenthal, R. (1979). The “File drawer problem” and Tolerance for Null Results. Psychological Bulletin, 86(3), 638–641. Schulz, K. F., Altman, D. G., & Moher, D. (2010). CONSORT 2010 Statement: updated guidelines for reporting parallel group randomised trials. BMJ, 340(mar23 1), c332–c332. https://doi.org/10.1136/bmj.c332 Shanahan, D. F., Bush, R., Gaston, K. J., Lin, B. B., Dean, J., Barber, E., & Fuller, R. A. (2016). Health Benefits from Nature Experiences Depend on Dose. Scientific Reports, 6(June), 1–10. https://doi.org/10.1038/srep28551 Shanahan, D. F., Fuller, R. A., Bush, R., Lin, B. B., & Gaston, K. J. (2015). The Health Benefits of Urban Nature: How Much Do We Need? BioScience, 65(5), 476–485. https://doi.org/10.1093/biosci/biv032 Shin, W. S., Shin, C. S., Yeoun, P. S., & Kim, J. J. (2011). The influence of interaction with forest on cognitive function. Scandinavian Journal of Forest Research, 26, 595–598. https://doi.org/10.1080/02827581.2011.585996 Song, C., Ikei, H., Igarashi, M., Miwa, M., Takagaki, M., & Miyazaki, Y. (2014). Physiological and psychological responses of young males during spring-time walks in urban parks. Journal of Physiological Anthropology, 33(8), 1–7. https://doi.org/10.1186/1880-6805-33-8 Song, C., Ikei, H., Igarashi, M., Takagaki, M., & Miyazaki, Y. (2015). Physiological and psychological effects of a walk in Urban parks in fall. International Journal of Environmental Research and Public Health, 12, 14216–14228. https://doi.org/10.3390/ijerph121114216 Song, C., Ikei, H., Kobayashi, M., Miura, T., Taue, M., Kagawa, T., … Miyazaki, Y. (2015). Effect of forest walking on autonomic nervous system activity in middle-aged hypertensive individuals: A pilot study. International Journal of Environmental Research and Public Health, 12, 2687– 2699. https://doi.org/10.3390/ijerph120302687 Song, C., Joung, D., Ikei, H., Igarashi, M., Aga, M., Park, B.-J., … Miyazaki, Y. (2013). Physiological and psychological effects of walking on young males in urban parks in winter. Journal of Physiological Anthropology, 32(18), 1–5. https://doi.org/10.1186/1880-6805-32-18 Sonntag-Öström, E., Nordin, M., Lundell, Y., Dolling, A., Wiklund, U., Karlsson, M., … Slunga Järvholm, L. (2014). Restorative effects of visits to urban and forest environments in patients with exhaustion disorder. Urban Forestry and Urban Greening, 13, 344–354. https://doi.org/10.1016/j.ufug.2013.12.007 Stekhoven, D. J., & Bühlmann, P. (2012). Missforest-Non-parametric missing value imputation for mixed-type data. Bioinformatics, 28(1), 112–118. https://doi.org/10.1093/bioinformatics/btr597 Sterne, J. A. C., Hernán, M. A., Reeves, B. C., Savović, J., Berkman, N. D., Viswanathan, M., … Higgins, J. P. (2016). ROBINS-I: A tool for assessing risk of bias in non-randomised studies of interventions. BMJ (Online), 355, 4–10. https://doi.org/10.1136/bmj.i4919 Sterne, J. A., & Egger, M. (2001). Funnel plots for detecting bias in meta-analysis: guidelines on choice of axis. Journal of Clinical Epidemiology, 54(10), 1046–1055. Stigsdotter, U. K., Corazon, S. S., Sidenius, U., Kristiansen, J., & Grahn, P. (2017). It is not all bad for the grey city – A crossover study on physiological and psychological restoration in a forest and an urban environment. Health and Place, 46, 145–154. https://doi.org/10.1016/j.healthplace.2017.05.007 Takayama, N., Korpela, K., Lee, J., Morikawa, T., Tsunetsugu, Y., Park, B. J., … Kagawa, T. (2014). Emotional, restorative and vitalizing effects of forest and urban environments at four sites in Japan. International Journal of Environmental Research and Public Health, 11, 7207–7230. https://doi.org/10.3390/ijerph110707207 Triguero-Mas, M., Gidlow, C. J., Martínez, D., De Bont, J., Carrasco-Turigas, G., Martínez-Íñiguez, T., … Nieuwenhuijsen, M. J. (2017). The effect of randomised exposure to different types of natural outdoor environments compared to exposure to an urban environment on people with indications of psychological distress in Catalonia. PLoS ONE, 12(3), 1–17. https://doi.org/10.1371/journal.pone.0172200 Tsunetsugu, Y., Lee, J., Park, B. J., Tyrväinen, L., Kagawa, T., & Miyazaki, Y. (2013). Physiological and psychological effects of viewing urban forest landscapes assessed by multiple measurements. Landscape and Urban Planning, 113, 90–93. https://doi.org/10.1016/j.landurbplan.2013.01.014 Tyrväinen, L., Ojala, A., Korpela, K., Lanki, T., Tsunetsugu, Y., & Kagawa, T. (2014). The influence of urban green environments on stress relief measures: A field experiment. Journal of Environmental Psychology, 38, 1–9. https://doi.org/10.1016/j.jenvp.2013.12.005 Ulrich, R. (1984). View through a window may influence recovery from surgery. Science, 224(4647), 420–421. https://doi.org/10.1126/science.6143402 Ulrich, R., Simonst, R. F., Lositot, B. D., Fioritot, E., Milest, M. A., & Zelsont, M. (1991). Stress Recovery During Exposure To Natural and Urban Environments. Journal of Environmental Psychology, 11, 201–230. Van den Berg, A. E. (2017). From Green Space to Green Prescriptions: Challenges and Opportunities for Research and Practice. Frontiers in Psychology, 8(February), 8–11. https://doi.org/10.3389/fpsyg.2017.00268 Van Den Berg, M., Wendel-vos, W., Poppel, M. Van, Kemper, H., Mechelen, W. Van, & Maas, J. (2015). Health benefits of green spaces in the living environment: A systematic review of epidemiological studies. Urban Forestry & Urban Greening, 14(4), 806–816. https://doi.org/10.1016/j.ufug.2015.07.008 van den Bosch, M., & Meyer-Lindenberg, A. (2019). Environmental Exposures and Depression: Biological Mechanisms and Epidemiological Evidence. Annual Review of Public Health, 40(1), 1–21. https://doi.org/10.1146/annurev-publhealth-040218-044106 Van den Noortgate, W., López-López, J. A., Marín-Martínez, F., & Sánchez-Meca, J. (2014). Meta- analysis of multiple outcomes: a multilevel approach. Behavior Research Methods, 47(4), 1274– 1294. https://doi.org/10.3758/s13428-014-0527-2 van Lissa, C. J. (2017). MetaForest: Exploring heterogeneity in meta-analysis using random forests. Open Science Framework. https://doi.org/10.17605/OSF.IO/KHJGB Viechtbauer, W. (2010). Conducting Meta-Analyses in R with the metafor Package. Journal of Statistical Software, 36(3). https://doi.org/10.18637/jss.v036.i03 World Health Organization. (2017). Depression and Other Common Mental Disorders: Global Health Estimates. SUPPLEMENTARY MATERIALS Table S1. PRISMA checklist Reported Section/topic # Checklist item on page # TITLE Title 1 Identify the report as a systematic review, meta-analysis, or both. 1 ABSTRACT Structured summary 2 Provide a structured summary including, as applicable: background; objectives; data sources; study eligibility criteria, 2 participants, and interventions; study appraisal and synthesis methods; results; limitations; conclusions and implications of key findings; systematic review registration number. INTRODUCTION Rationale 3 Describe the rationale for the review in the context of what is already known. 3-4 Objectives 4 Provide an explicit statement of questions being addressed with reference to participants, interventions, comparisons, 5 outcomes, and study design (PICOS). METHODS Protocol and registration 5 Indicate if a review protocol exists, if and where it can be accessed (e.g., Web address), and, if available, provide 5 registration information including registration number. Eligibility criteria 6 Specify study characteristics (e.g., PICOS, length of follow-up) and report characteristics (e.g., years considered, 5 language, publication status) used as criteria for eligibility, giving rationale. Information sources 7 Describe all information sources (e.g., databases with dates of coverage, contact with study authors to identify 6 additional studies) in the search and date last searched. Search 8 Present full electronic search strategy for at least one database, including any limits used, such that it could be Supp repeated. material Study selection 9 State the process for selecting studies (i.e., screening, eligibility, included in systematic review, and, if applicable, 6 included in the meta-analysis). Data collection process 10 Describe method of data extraction from reports (e.g., piloted forms, independently, in duplicate) and any processes 6 for obtaining and confirming data from investigators. Data items 11 List and define all variables for which data were sought (e.g., PICOS, funding sources) and any assumptions and 6 simplifications made. Risk of bias in individual 12 Describe methods used for assessing risk of bias of individual studies (including specification of whether this was 6-7 studies done at the study or outcome level), and how this information is to be used in any data synthesis. Summary measures 13 State the principal summary measures (e.g., risk ratio, difference in means). 7 Reported Section/topic # Checklist item on page # Synthesis of results 14 Describe the methods of handling data and combining results of studies, if done, including measures of consistency 7-8 (e.g., I ) for each meta-analysis. Risk of bias across studies 15 Specify any assessment of risk of bias that may affect the cumulative evidence (e.g., publication bias, selective 6-7,9 reporting within studies). Additional analyses 16 Describe methods of additional analyses (e.g., sensitivity or subgroup analyses, meta-regression), if done, indicating 7-8 which were pre-specified. RESULTS Study selection 17 Give numbers of studies screened, assessed for eligibility, and included in the review, with reasons for exclusions at 10 + Supp each stage, ideally with a flow diagram. material Study characteristics 18 For each study, present characteristics for which data were extracted (e.g., study size, PICOS, follow-up period) and 10-18 provide the citations. Risk of bias within studies 19 Present data on risk of bias of each study and, if available, any outcome level assessment (see item 12). 18-19 Results of individual studies 20 For all outcomes considered (benefits or harms), present, for each study: (a) simple summary data for each 21-23 intervention group (b) effect estimates and confidence intervals, ideally with a forest plot. Synthesis of results 21 Present results of each meta-analysis done, including confidence intervals and measures of consistency. 21-23 Risk of bias across studies 22 Present results of any assessment of risk of bias across studies (see Item 15). 23+26 Additional analysis 23 Give results of additional analyses, if done (e.g., sensitivity or subgroup analyses, meta-regression [see Item 16]). 24- 26+Supp Material DISCUSSION Summary of evidence 24 Summarize the main findings including the strength of evidence for each main outcome; consider their relevance to 26-27 key groups (e.g., healthcare providers, users, and policy makers). Limitations 25 Discuss limitations at study and outcome level (e.g., risk of bias), and at review-level (e.g., incomplete retrieval of 27-28 identified research, reporting bias). Conclusions 26 Provide a general interpretation of the results in the context of other evidence, and implications for future research. 26-29 FUNDING Funding 27 Describe sources of funding for the systematic review and other support (e.g., supply of data); role of funders for the 29 systematic review. Table S2. ROBINS-I assessments Joung et al. 2015 No Low Low Low Low Serious Moderate Serious information Li et al., (2016) Serious Low Low Low Low Serious Serious Serious Roe & Aspinall Moderate Low Low Low Serious Serious Moderate Serious (2011) Song et al., Serious Low Low Low Low Serious Serious Serious (2013) Song et al., Moderate Low Low Low Low Serious Moderate Serious (2014) Song et al., Moderate Low Low Low Low Serious Moderate Serious (2015b) Stigsdotter et al., Serious Low Low Low Low Serious Moderate Serious (2017) Tsunetsugu et al. No Low Low Low Low Serious Moderate Serious (2013) information Confounding Participant selection Classification of interventions Deviation from intended intervention Missing data Measurement of outcomes Reporting bias Overall bias Figure S1. Flow diagram of study selection PRISMA 2009 Flow Diagram Records identified through Additional records identified database searching through other sources (n= 8,958) (n=9) Records after duplicates removed (n=6,622) Records screened Records excluded (n=6,622) (n=6,546) Full-text articles assessed Full-text articles excluded, for eligibility with reasons (n=76) (n=52) Review papers: 17 Does not measure depression: 6 No comparison with urban Studies included in environment: 29 quantitative synthesis (n=33) Figure S2. Risk of bias graph and table 1.3 Included 1.2 Eligibility 1.4 Screening 1.1 Identification Figure S3. Summary of Findings table № of patients Effect Certainty assessment Certainty Importance № of Other the natural the urban Relative Absolute Study design Risk of bias Inconsistency Indirectness Imprecision studies considerations environment environment (95% CI) (95% CI) Depressive mood 25 randomised very serious serious not serious not serious all plausible residual 912 902 - SMD 0.05 ⨁⨁◯◯ IMPORTANT trials confounding would SD higher LOW suggest spurious (0.4 lower effect, while no effect to 0.5 was observed higher) Depressive mood b c 8 observational very serious serious not serious serious all plausible residual 179 175 - SMD 0.38 ⨁◯◯◯ IMPORTANT studies confounding would SD lower VERY LOW suggest spurious (0.59 lower effect, while no effect to 0.16 was observed lower) CI: Confidence interval; SMD: Standardised mean difference Explanations a. Issues with randomisation method, allocation concealment, blinding, incomplete outcome data and selective reporting. b. Random effects meta-analysis found residual heterogeneity between studies. c. Small sample size and wide confidence intervals. Figure S4. Replicated MetaForest for variable preselection Figure S5. Variable importance for final model Note. The MetaForest analysis was tuned using the eight moderators above. The optimal model used uniform weights, 2 candidate variables at each split, and a minimal terminal node size of 3. Figure S6. Partial dependence plots Search strategies for all databases searched MEDLINE ((((people[MeSH Major Topic] OR "people" OR adults[MeSH Major Topic] OR "adults" OR volunteers[MeSH Major Topic] OR "volunteers" OR participants OR subjects OR students OR respondents)) AND (intervention[MeSH Major Topic] OR "intervention" OR experiment OR randomized[MeSH Major Topic] OR "randomized" OR crossover OR "case-crossover" OR "pre-post" OR comparison OR "non-randomized" OR "quasi-experiment" OR controlled OR control OR "control group")) AND ("natural environment" OR "natural outdoor environment" OR outdoors OR outside OR "green space" OR forest[MeSH Major Topic] OR "forest" OR "woodland" OR garden OR allotment OR countryside OR "open space" OR landscape OR parkland OR park NOT parkin*)) AND ("mental health" OR depression[MeSH Major Topic] OR "depression" OR mood[MeSH Major Topic] OR "mood" OR "mood disorder" OR "mood change" OR "major depressive disorder" OR "depressive symptoms") Sort by: Author Filters: Humans; English PsychINFO via OvidSP people.mp. adults.mp. volunteers.mp. experimental subjects/ participants.mp. subjects.mp. students.mp. respondents.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] 1 or 2 or 3 or 4 or 5 or 6 or 7 or 8 intervention/ intervention.mp. experiment.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] randomized.mp. crossover.mp. "case-crossover".mp. "pre-post".mp. comparison.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "non-randomized".mp. controlled.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] control.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "control group".mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] 10 or 11 or 12 or 13 or 14 or 15 or 16 or 17 or 18 or 19 or 20 or 21 forest.mp. woodland.mp. "natural environment".mp. "natural outdoor environment".mp. outdoors.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] outside.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "green space".mp. garden.mp. allotment.mp. countryside.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] (park not parkin*).mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] parkland.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "open space".mp. landscape.mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "urban area".mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "urban environment".mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] urban environments/ exp Built Environment/ 23 or 24 or 25 or 26 or 27 or 28 or 29 or 30 or 31 or 32 or 33 or 34 or 35 or 36 or 37 or 38 or 39 or 40 "mental health".mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "depression (emotion)"/ exp MAJOR DEPRESSION/ depression.mp. mood.mp. "mood change".mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "depressive symptoms".mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] "mood disorder".mp. [mp=title, abstract, heading word, table of contents, key concepts, original title, tests & measures] 42 or 43 or 44 or 45 or 46 or 47 or 48 or 49 9 and 22 and 41 and 50 EMBASE #1 participants OR 'students'/exp OR students OR subjects OR 'volunteer'/exp OR 'volunteer' OR people OR respondents #2 'control group' OR control OR 'non-randomized' OR 'controlled clinical trial (topic)' OR 'pre-post' OR comparison OR 'case-crossover' OR 'crossover procedure' OR 'randomized' OR 'randomized controlled trial' OR 'experiment' OR 'intervention study' #3 ('open space' OR 'allotment' OR 'garden' OR 'green space' OR 'natural outdoor environment' OR 'natural environment' OR outdoors OR outside OR 'forest' OR woodland OR countryside OR landscape OR 'recreational park' OR parkland OR 'urban area' OR 'urban environment') NOT parkin* #4 'mental health' OR 'mood change' OR 'depressive symptoms' OR 'mood' OR 'mood disorder' OR 'major depression' OR 'depression' #1 AND #2 AND #3 AND #4 AND [humans]/lim AND [english]/lim Scopus TITLE-ABS-KEY ( people OR volunteers OR subjects OR students OR participants OR respondents ) AND ( TITLE-ABS-KEY ( intervention OR experiment OR "randomized controlled trial" OR randomized OR crossover OR "case-crossover" OR pre-post OR comparison OR non-randomized OR exposure OR controlled OR control OR "control group" ) ) AND ( TITLE-ABS-KEY ( "green space" OR "natural outdoor environment" OR "natural environment" OR outdoors OR outside OR "open space" OR countryside OR allotment OR garden OR forest OR woodland OR landscape OR parkland OR park AND NOT parkin* ) ) AND ( TITLE-ABS-KEY ( "mental health" OR depression OR depressive OR "depressive symptoms" OR mood OR "mood disorder" ) ) AND ( LIMIT-TO ( LANGUAGE , "English " ) ) Web of Science (participants OR students OR subjects OR volunteers OR people) AND ('control group' OR 'non-randomized' OR 'pre-post' OR 'case-crossover' OR 'crossover procedure' OR randomized OR 'randomized controlled trial' OR intervention) AND (park OR 'open space' OR allotment OR garden OR 'green space' OR 'natural outdoor environment' OR 'natural environment' OR forest OR woodland) AND ('mood change' OR 'depressive symptoms' OR mood OR 'mood disorder' OR 'major depression' OR depression)
Statistics – arXiv (Cornell University)
Published: Jul 23, 2019
Access the full text.
Sign up today, get DeepDyve free for 14 days.