City-level greenness exposure is associated with COVID-19 incidence in China

Accumulating studies have suggested an important role of environmental factors (e.g. air pollutants) on the occurrence and development of coronavirus disease 2019 (COVID-19). Evidence concerning the relationship of greenness on COVID-19 is still limited. This study aimed to assess the association between greenness and COVID-19 incidence in 266 Chinese cities. A total of 12,377 confirmed COVID-19 cases were identified through February 29th, 2020. We used the average normalized difference vegetation index (NDVI) during January and February 2020 from MOD13A2 product, to represent the city-level greenness exposure. A generalized linear mixed-effects model was used to estimate the association between NDVI exposure and COVID-19 incidence using COVID-19 cases as the outcome. We evaluated whether the association was modified by population density, GDP per capita, and urbanization rate, and was mediated by air pollutants. We also performed a series of sensitivity analyses to discuss the robustness of our results. Per 0.1 unit increment in NDVI was negatively associated with COVID-19 incidence (IRR: 0.921, 95% CI: 0.898, 0.944) after adjustment for confounders. Associations with COVID-19 incidence were stronger in cities with lower population density, lower GDP per capita, and lower urbanization rate. We failed to detect any mediation effect of air pollutants on the association between NDVI and COVID-19 incidence. Sensitivity analyses also indicated consistent estimates. In conclusion, our study suggested a beneficial association between city-level greenness and COVID-19 incidence. We could not establish which mechanisms may explain this relationship. •The study is the first to discuss the association of NDVI with COVID-19 incidence in 266 Chinese cities.•NDVI exposure was negatively associated with COVID-19 incidence.•Stronger protective effect was observed in cities with lower population density.•Stronger protective effect was observed in cities with lower GDP per capita.•Stronger protective effect was observed in cities with lower urbanization rate.