Abstract
Evidence of the benefits of greenspaces or greenness to human wellbeing in the context of COVID-19 is fragmented and sometimes contradictory. This calls for a meta-analysis of existing studies to clarify the matter. Here, we identified 621 studies across the world, which were then filtered down to 13 relevant studies covering Africa, Asia, Europe, and USA. These studies were meta-analysed, with the impacts of greenspaces on COVID-19 infection rate quantified using regression estimates whereas impacts on mortality was measured using mortality rate ratios. We found evidence of significant negative correlations between greenness and both COVID-19 infection and mortality rates. We further found that the impacts on COVID-19 infection and mortality are moderated by year of publication, greenness metrics, sample size, health and political covariates. This clarification has far-reaching implications on policy development towards the establishment and management of green infrastructure for the benefits of human wellbeing.
Introduction
Global human population is changing rapidly following an exponential growth path, putting tremendous pressures on natural resources. Currently, it is approximately 7.9 billion people1 and is predicted to reach above 9 billion by 2050 or 11 billion by 21002. In response, nature fights back in various ways to bring down global population to a sustainable level. One of these ways is through global pandemics, e.g., COVID-19. Indeed, the world has been witnessing COVID-19 pandemic since 2020, with over half a million of infection cases and over 20,000 deaths in 20203. In 2022, these figures grew tremendously, reaching over 600 million cumulative cases with over 6 million cumulative deaths4. Subsequently, various studies, using different metrics of greenness (the total amount of vegetation in an area), were conducted across the globe to investigate whether greenness act as buffer infrastructure against the spread of COVID-19 infection rates and severity.
However, the findings reported in these studies are mixed5,6. For example, ref.7 found that a 0.1 increase in NDVI is linked to 4.1% reduction in COVID-19 incidence rate ratio in the USA. A similar pattern was observed using street-level indicators of greenness8. The mitigating effects of greenness have also been reported elsewhere: in China and India, an increase in greenness shows a strong negative association with the spread of COVID-19 infections and mortalities9,10. These negative effects may be interpreted as follows: activities of the Natural Killer (NK) cells in human body are boosted with frequent exposure to vegetation3,11 – NK cells, as part of the immune system, attack to eliminate virus-infected cells12. Also, by safeguarding against air pollution, vegetation contributes to lower health risks that may aggravate the severity of COVID-19 infection13,14. Additionally, green infrastructure often provides spacious environment for physical exercise, recreation, and social events with reduced chances of person-to-person contact15,16. As opposed to these negative correlations between greenness and COVID-19 infection rates, reports of positive correlations are also documented. For example, ref.6 found that urban greenspaces were associated with an increase in the spread of COVID-19 infections (see also ref.15,16).
These mixed findings could be linked to the differences in how COVID-19 severity was measured, e.g., as hospitalization rates, mortality rates, admission rate to ICU, etc. Additional sources of differences in findings may be linked to differences in sample size, type and number of covariates considered, and choice of statistical tests17,18. Furthermore, the mixed findings may be linked to differences in how greenness was measured in different studies. Indeed, greenness was variously measured as street trees, botanical gardens, natural forests and grasslands, and residential gardens or as amount of greenness captured in NDVI or EVI or as quality of green infrastructure7,8,19,20. For example, ref.15 measured greenness as ‘green space density’ which is the proportion of specific vegetation types in a given spatial unit which they correlated with COVI-19 infection risk measured as ‘venue density’ (number of buildings visited by confirmed COVID-19 positive cases). Since greenspaces are attraction sites, they attract increasing number of visitors, thus increasing the infection risks, and leading to a positive correlation between greenness and infection rate15. Furthermore, the mixed findings may be linked to the use of various confounding factors in the model of COVID-19 infection and mortality rates. These factors may be age21, ethnicity22, and poverty level23, among others.
The emergence of conflicting findings presents a challenge with regards to the generalization of the benefits of greenness to human wellbeing in the context of COVID-19 pandemic. In such context, a meta-analysis of existing evidence presents an opportunity to integrate the conflicting reported effects of greenness on COVID-19 infection rates and severity to investigate whether generalization is possible. Scientifically rigorous methodologies are increasingly adopted in various studies to improve the validity of findings and lower between-study heterogeneity. These include the use of larger sample sizes, use of multiple predictors, choice of relevant statistical tests and covariates, and use of fine spatial scales23-26. Regardless of these advances, consolidation of measured effect sizes and determination of between-study heterogeneity is still needed. To date, several studies have investigated the relationships between the provision and quantity of greenness and its effects on the spread and severity of COVID-195,7,10,20. However, in the context of conflictual findings reported, a meta-analysis imposes itself or becomes an obligation if we are to clarify how greenness or green infrastructure relates to COVID-19. In the present study, our main objective is to provide such clarifications.
Results
Characteristics of studies included in the meta-analysis
A total of 621 studies across the world (Figure 1A) were identified through the search of Scopus, PubMed, and Google scholar platforms. After removing irrelevant and duplicate studies, 25 studies remained, covering Africa, Asia, Europe, and USA (Figure 1B). A review of the 25 full-text articles resulted in a removal of 12 studies that were either review/commentary in nature or did not report the required statistical parameters for meta-analysis.
The study characteristics are summarized in Table S1. Nine studies that tested the relationships between greenness and COVID-19 infections and four studies that investigated the relationships between greenness and COVID-19 mortality rates were included in the final synthesis. Most of the studies (nine out of 13) were conducted in the United States of America (USA) whereas China, England, India, and South Africa each had one study (Figure 1B). A total of 7 out of 13 studies used more than one predictor of COVID-19 impact in each study with normalised difference vegetation index (NDVI) and abundance of greenness as the mostly used measures of greenness (Figure 2A). Because multiple predictors are used in a single study, a total of 45 different correlations between infection rates and greenness were produced in all 13 studies and 14 correlations between mortality rates and greenness were produced from four studies of COVID-19. We classified covariates into five broad groups: climatic, demographic, economic, health, and political. All 13 studies considered at least one demographic covariate in their analyses, and only four studies included climatic, demographic, economic, health, and political covariates (Figure 2B).
Greenness and COVID-19 infections
We found a statistically significant negative effect of greenness on COVID-19 infections (β = −0.08, 95% CI: −0.1396 – −0.0252; t=-2.90; p=0.006) with a prediction interval of [−0.3601 – 0.1954] (95% CI) (Figure 3). Between-study heterogeneity variance was estimated at τ2= 0.0184 (95% CI: 0.0185 – −0.0813), with an I2 value of 94.1% (95% CI: 92.9% – 95.1%). Subgroup analyses reveal that between-study heterogeneity can be attributed to year of publication (X2=8.24; p=0.02), choice of predictors (X2=129.68; p<0.01), and use of political covariates (X2=8.27; p<0.01) (see Table 1 and Figures S1 – S6).
Greenness and COVID-19 mortalities
We found that an increase in greenness was strongly linked to lower mortality rate ratio (MRR= 0.9272; 95% CI: 0.8788 – 0.9783; t=-3.05; p=0.009) with a prediction interval of [0.7683 – 1.1189] (95% CI) (Figure 4). Furthermore, an estimated 0.0069 between-study heterogeneity variance (95% CI: 0.0032 – 0.0228) was observed with an I2 value of 92% (95% CI: 88.3% – 94.5%). We also found that year of publication (X2=19.10; p<0.01), sample size (X2=7.92; p<0.01), choice of predictors (X2=14.92; p<0.01), and use of health (X2=7.92; p<0.01) and political (X2=22.75; p<0.01) covariates strongly impact the degree of heterogeneity (see Table 2 & Figures S7 – S13).
Publication bias
Existence of publication bias was investigated using the Funnel approach and Orwin fail-safe number. The presence of funnel plot symmetry (Figure 5A) indicated a lack of publication bias for studies that investigate the effect of greenness on COVID-19 infections (Fail-safe N: 45). Publication bias was, however, observed for studies that test the relationship between greenness and COVID-19 mortalities (Figure 5B; Fail- safe N: 14).
Discussion
Our meta-analysis provides evidence that an increase in abundance or exposure to greenness is associated with a significant reduction in COVID-19 infection rates and death cases10,19,27. However, we found high heterogeneity between the studies that were included in the meta-analysis. Subgroup analyses revealed that heterogeneity in studies on COVID-19 infections and mortality is strongly predicted by the studies’ years of publication, choices of predictors (metrics of greenness), and inclusion of political covariates. Additionally, sample size and consideration of health covariates strongly affect heterogeneity of studies on COVID-19 mortalities.
The sensitivity of effect size to year of publication can be attributed to availability of data to adequately model the impact of COVID-19. The spread of COVID-19 and increased global testing for COVID-19 infection accelerated overtime, thus allowing successive studies to have an increasingly larger data pool28,29. This may also impact sample sizes that are adopted in each study. As more regions produce more data on COVID-19 infections and mortality, their eligibility to be included in studies investigating the correlations between COVID-19 and greenness may enhance study designs. In our subgroup analysis, we found that studies that used smaller sample sizes (n<2000) are likely to report larger effect sizes compared to studies with larger sample sizes. Given the importance of selecting an appropriate sample size30, the need to define an appropriate sample size for investigating the health benefits of green infrastructure remains critical.
The diversity of greenness metrics, ranging from street trees to large forests, presents a unique challenge while measuring their impacts. Commonly, studies that cover large study areas use vegetation indices such as NDVI or EVI which are retrieved from satellite imagery31-33. Since health benefits of greenness are usually felt closer to the greenness34,35, several studies consider local greenness such as household gardens36, street trees37,38, and local parks39,40 in their analysis. However, this approach is only feasible when focusing on smaller areas. In some cases, subjective measures of greenness were used41,42. Our findings in the present study suggest that the choice of greenness metrics adopted in different studies affects the its effect size. The use of NDVI, EVI or vegetation canopy size produces large effects of greenness against COVID-19 infections and mortalities. In contrast, studies that use proximity or visitation patterns are likely to report marginal effects.
All studies in our meta-analysis have included demographic covariates, and 92% of studies included economic covariates. While modelling the effects of greenness, the inclusion of demographic variables such as population density and age structure, as well as economic indicators such as gross domestic product (GDP) and household income level as covariates have been largely adopted7,9,27,43. Furthermore, the use of health covariates featured in several studies5,19,26. However, consideration of political covariates in the modelling of greenness benefits to human wellbeing in the context of COVID-19 is only starting to emerge5,6. Political factors such as promulgation of mobility restrictions15,44 and face-masks mandates45 have shown to be significant predictors of COVID-19 impacts, although their inclusion in studies linking greenness to COVID-19 infection and severity remains limited. We found that the use of political covariates significantly affects the effect size. Inclusion of political covariables resulted in a greater effect size in studies of COVID-19 mortality and in a smaller effect size in studies of COVID-19 infections. This may suggest that existing policies are more effective in reducing COVID-19 fatalities than curbing the spread of infections.
Overall, meta-analysing studies from Africa, Asia, Europe, and USA, we found strong support for beneficial effects of greenness to human in the face of COVID-19 infection and severity, suggesting that positive correlations reported in some studies between greenness vs. infection and mortality rates15,16 might simply imply that the greenness metrics used in those studies (e.g., green space density or accessibility to greenspaces) do not fully capture important facets of greenness. This calls for a need to homogenize greenness metrics in studies to come. There is also a need for homogenization of COVID-19 severity metrics since we could not include hospitalization rate in the present study as a measure of COVID-19 severity because very limited studies have investigated hospitalization rate. Lastly, our results showed high degree of between-study heterogeneity which can be explained by year of publication, sample size, and choice of predictor variables and covariates. However, evidence from existing studies show that green infrastructure moderates the impact of COVID-19 by reducing prevalence of infections and associated mortalities.
Nevertheless, our findings have some far-reaching implications for the establishment and management of green infrastructure: greenspaces must be acknowledged as critical infrastructure that has substantial broader public health values, and as such, deserve enough fundings from governments worldwide, especially in the developing world.
Methods
Study selection
The Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRIMSA) guidelines46 were followed to search for literature that focus on green infrastructure and its impact on COVID-19. All search results were reviewed for relevance based on their title and abstract to be considered for meta-analysis (Figure 6). Furthermore, reference lists of all included articles were reviewed to identify studies that meet the inclusion criteria.
Search strategy
Literature search was limited to PubMed, Scopus, and Google Scholar. The following search string was used to search for literature on the 17th of April 2023: (“Greenspace“ or “green space“ or “greenery“ or “greenness“ or “vegetation“ or “trees“ or “forest” or “grass” or “grassland”) and (“COVID-19“ or “SARS-CoV-2“ or “coronavirus“ or “COVID”). We did not apply any restrictions on publication date in the search.
Eligibility criteria
Inclusion criteria for this study were as follow: (a) original research that investigates effects of green infrastructure on COVID-19 infections and related mortalities; (b) full- text is available; (c) publication is in English; (d) required statistical parameters for meta-analysis are reported in the main article or supplementary files (i.e., regression estimates for predicting COVID-19 infections, and mortality rate ratios for predicting COVID-19 mortalities). Exclusion criteria were review or commentary articles, articles without required parameters, and articles not in English (see Figure 6).
Data extraction and analysis
A predetermined template was used to collect study characteristics which are surname of first author, year of publication, country of study, measure of green infrastructure, temporal extent of study, sample size, measure of COVID-19, effect type, effect size, standard error or confidence interval, and list of covariates. All data analysed in this study are available as Supplemental Information (Appendix 1).
All analyses were conducted in R version 4.2.347 (see R script in Supplemental Information). Two separate meta-analyses were conducted, focusing on impacts of greenness on COVID-19 infections (meta-analysis 1) and COVID-19 mortalities (meta-analysis 2). Regression estimates were used as pre-calculated effect type when analysing COVID-19 infections, and mortality rate ratios (MRR) were used as pre-calculated effect type when analysing COVID-19 mortalities. Subsequently, subgroup analyses were applied to the same data to test the effects of predictor variables, sample size, and selection of covariates. Random models were selected in each analysis using metagen function found in the “Metafor” R library48.
Outcomes are reported as pooled regression estimates for COVID-19 infections and as pooled MRR for COVID-19 deaths. Furthermore, in each case, a 95% confidence interval (CI), t-value, and p-values are reported with p<0.005 considered as an indicator of statistical significance. Between-studies heterogeneity was quantified using Higgins & Thompson’s I2 statistic49 with the I2 value of less than 25%, 50% and 75% indicating low, moderate, or high heterogeneity, respectively. Heterogeneity variance and prediction interval were also reported to measure the extent of between-study heterogeneity. Publication bias was tested using the Funnel approach50 (Sterne & Egger 2001) and the Orwin’s fail-safe number51.
Data Availability
All data produced in the present work are contained in the manuscript.
Author Contributions
KY conceived the project, BP collected the data, BP analysed the data, KY and BP wrote the paper.
Data Availability
All data generated or analysed in this study are included in this article (and its Supplementary Information files).
Conflict of Interest Statement
The authors declare no conflict of interests.
Code availability
Code to replicate all results in this paper is available as Supplemental Information.
SUPPLEMENTAL INFORMATION
TABLES
Table S1. Characteristics of all studies included in the present meta-analysis.
Figures S1-S6. Subgroup analyses of between-study heterogeneity can be attributed to year of publication (X2=8.24; p=0.02), choice of predictors (X2=129.68; p<0.01), and use of political covariates (X2=8.27; p<0.01).
Figures S7-S13. Between-study heterogeneity of variance showing that year of publication (X2=19.10; p<0.01), sample size (X2=7.92; p<0.01), choice of predictors (X2=14.92; p<0.01), and use of health (X2=7.92; p<0.01) and political (X2=22.75; p<0.01) covariates strongly impact the degree of heterogeneity.
Appendices
Appendix 1. Data collected and analysed in this study.
Appendix 2. R script used to reproduce the present study.
Acknowledgements
KY acknowledges the South Africa’s National Research Foundation (grant # SRUG22051210107).