Air pollution and COVID-19 mortality in the United States: strengths and limitations of an ecological regression analysis (Updated April 24, 2020)
Notice: In the revision on April 24, 2020, we have updated our analysis using data up to April 22, and importantly in which we have adjusted for additional confounding factors that also reflect the timing of the epidemic's spread, the timing of the social distancing policies and the population age distribution. Consequently, we have revised our finding as that an increase of 1 μg/m3 in PM2.5 is associated with an 8% increase in the COVID-19 death rate (95% confidence interval [CI]: 2%, 15%).
Figure: COVID-19 mortality rate ratios (MRR) per 1 μg/m3 increase in PM2.5 and 95% CI using daily cumulative COVID-19 death counts from April 18, 2020 to September 7, 2020.
Xiao Wu MS, Rachel C. Nethery PhD, M. Benjamin Sabath MA, Danielle Braun PhD, Francesca Dominici PhD
All authors are part of the Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, 02115, USA
Lead authors: Xiao Wu and Rachel C. Nethery
Corresponding and senior author: Francesca Dominici, PhD
Background: United States government scientists estimate that COVID-19 may kill tens of thousands of Americans. Many of the pre-existing conditions that increase the risk of death in those with COVID-19 are the same diseases that are affected by long-term exposure to air pollution. We investigated whether long-term average exposure to fine particulate matter (PM2.5) is associated with an increased risk of COVID-19 death in the United States.
Design: A nationwide, cross-sectional study using county-level data.
Data sources: COVID-19 death counts were collected for more than 3,000 counties in the United States (representing 98% of the population) up to April 22, 2020 from Johns Hopkins University, Center for Systems Science and Engineering Coronavirus Resource Center.
Methods: We fit negative binomial mixed models using county-level COVID-19 deaths as the outcome and county-level long-term average of PM2.5 as the exposure. In the main analysis, we adjusted by 20 potential confounding factors including population size, age distribution, population density, time since the beginning of the outbreak, time since state’s issuance of stay-at-home order, hospital beds, number of individuals tested, weather, and socioeconomic and behavioral variables such as obesity and smoking. We included a random intercept by state to account for potential correlation in counties within the same state. We conducted more than 68 additional sensitivity analyses.
Results: We found that an increase of only 1 μg/m3 in PM2.5 is associated with an 8% increase in the COVID-19 death rate (95% confidence interval [CI]: 2%, 15%). The results were statistically significant and robust to secondary and sensitivity analyses.
Conclusions: A small increase in long-term exposure to PM2.5 leads to a large increase in the COVID-19 death rate. Despite inherent limitations of the ecological study design, our results underscore the importance of continuing to enforce existing air pollution regulations to protect human health both during and after the COVID-19 crisis. The data and code are publicly available so our analyses can be updated routinely.
Data and Code:
Our data and code is available on github here. (Updated Sep 18, 2020)
Manuscript and Supplemental Material
- By using the contents on this website and the Github repo, you agree to cite:
Wu, X., Nethery, R.C., Sabath, B.M., Braun, D. and Dominici, F., 2020. Air pollution and COVID-19 mortality in the United States: strengths and limitations of an ecological regression analysis. Science Advances (in press).
A pre-print version can be found at: Exposure to air pollution and COVID-19 mortality in the United States. Xiao Wu, Rachel C. Nethery, Benjamin M. Sabath, Danielle Braun, Francesca Dominici. medRxiv 2020.04.05.20054502; doi: https://doi.org/10.1101/2020.04.05.20054502
We appreciate the work of Aaron Van Donkelaar, Randall Martin, and his team for providing us with access to their estimates of PM2.5 exposure. Their data (V4.NA.02.MAPLE) can be found on Randall Martin's website here: https://sites.wustl.edu/acag/datasets/surface-pm2-5/
The data was produced as part of the following paper:
van Donkelaar, A., R. V. Martin, C. Li, R. T. Burnett, Regional Estimates of Chemical Composition of Fine Particulate Matter using a Combined Geoscience-Statistical Method with Information from Satellites, Models, and Monitors, Environ. Sci. Technol., doi: 10.1021/acs.est.8b06392, 2019.
We would like to thank Lena Goodwin and Stacey Tobin for editorial assistance in the preparation of this manuscript.