Effect modification of greenness on the association between heat and mortality: A multi-city multi-country study

Identifying how greenspace impacts the temperature-mortality relationship in urban environments is crucial, especially given climate change and rapid urbanization. However, the effect modification of greenspace on heat-related mortality has been typically focused on a localized area or single country. This study examined the heat-mortality relationship among different greenspace levels in a global setting.

We collected daily ambient temperature and mortality data for 452 locations in 24 countries and used Enhanced Vegetation Index (EVI) as the greenspace measurement. We used distributed lag non-linear model to estimate the heat-mortality relationship in each city and the estimates were pooled adjusting for city-specific average temperature, city-specific temperature range, city-specific population density, and gross domestic product (GDP). The effect modification of greenspace was evaluated by comparing the heat-related mortality risk for different greenspace groups (low, medium, and high), which were divided into terciles among 452 locations.

Cities with high greenspace value had the lowest heat-mortality relative risk of 1·19 (95% CI: 1·13, 1·25), while the heat-related relative risk was 1·46 (95% CI: 1·31, 1·62) for cities with low greenspace when comparing the 99^th temperature and the minimum mortality temperature. A 20% increase of greenspace is associated with a 9·02% (95% CI: 8·88, 9·16) decrease in the heat-related attributable fraction, and if this association is causal (which is not within the scope of this study to assess), such a reduction could save approximately 933 excess deaths per year in 24 countries.

Our findings can inform communities on the potential health benefits of greenspaces in the urban environment and mitigation measures regarding the impacts of climate change.

This publication was developed under Assistance Agreement No. RD83587101 awarded by the U.S. Environmental Protection Agency to Yale University. It has not been formally reviewed by EPA. The views expressed in this document are solely those of the authors and do not necessarily reflect those of the Agency. EPA does not endorse any products or commercial services mentioned in this publication. Research reported in this publication was also supported by the National Institute on Minority Health and Health Disparities of the National Institutes of Health under Award Number R01MD012769. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Also, this work has been supported by the National Research Foundation of Korea (2021R1A6A3A03038675), Medical Research Council-UK (MR/V034162/1 and MR/R013349/1), Natural Environment Research Council UK (Grant ID: NE/R009384/1), Academy of Finland (Grant ID: 310372), European Union's Horizon 2020 Project Exhaustion (Grant ID: 820655 and 874990), Czech Science Foundation (22-24920S), Emory University's NIEHS-funded HERCULES Center (Grant ID: P30ES019776), and Grant CEX2018-000794-S funded by MCIN/AEI/ 10.13039/501100011033 The funders had no role in the design, data collection, analysis, interpretation of results, manuscript writing, or decision to publication.