Spatial associations of long-term exposure to diesel particulate matter with seasonal and annual mortality due to COVID-19 in the contiguous United States
BMC Public Health volume 23, Article number: 423 (2023)
People with certain underlying respiratory and cardiovascular conditions might be at an increased risk for severe illness from COVID-19. Diesel Particulate Matter (DPM) exposure may affect the pulmonary and cardiovascular systems. The study aims to assess if DPM was spatially associated with COVID-19 mortality rates across three waves of the disease and throughout 2020.
We tested an ordinary least squares (OLS) model, then two global models, a spatial lag model (SLM) and a spatial error model (SEM) designed to explore spatial dependence, and a geographically weighted regression (GWR) model designed to explore local associations between COVID-19 mortality rates and DPM exposure, using data from the 2018 AirToxScreen database.
The GWR model found that associations between COVID-19 mortality rate and DPM concentrations may increase up to 77 deaths per 100,000 people in some US counties for every interquartile range (0.21 μg/m3) increase in DPM concentration. Significant positive associations between mortality rate and DPM were observed in New York, New Jersey, eastern Pennsylvania, and western Connecticut for the wave from January to May, and in southern Florida and southern Texas for June to September. The period from October to December exhibited a negative association in most parts of the US, which seems to have influenced the year-long relationship due to the large number of deaths during that wave of the disease.
Our models provided a picture in which long-term DPM exposure may have influenced COVID-19 mortality during the early stages of the disease. That influence appears to have waned over time as transmission patterns evolved.
In 2020, more than 20 million cases of coronavirus disease 2019 (COVID-19) were identified in the United States (U.S.), and more than 350,000 people died [5, 9, 19]. In addition to age, socioeconomic status, access to healthcare, physical environment, and education have been identified as social determinants of COVID-19 hospitalization and mortality [30, 36]. Several studies have observed a disproportionate share of COVID-19 incidence and mortality among predominantly Black U.S. communities, which may be partly attributable to social and economic inequalities and preexisting comorbidities [7, 8, 33 36, 40, 55], as well as to disproportionately high exposures to air pollution .
The impact of particulate matter exposures on COVID-19 outcomes have also been evaluated, with some studies centered on diesel particulate matter (DPM). In an investigation of the role of long-term exposure (2000-2016) to air pollution during the first months of the pandemic, Wu et al.  found that an increase of 1 μg/m3 in particulate matter with a nominal diameter of 2.5 μm (PM2.5) was associated with an 11% increase in the COVID-19 death rate for January 1-June 18, 2020. Bozack et al.  performed a similar analysis to test associations of COVID-19 intensive care unit (ICU) admission and mortality with long-term concentrations of PM2.5, nitrogen dioxide, and black carbon for the period March 8-August 30, 2020 in New York City. They noted an association of ICU admission and mortality with long-term PM2.5 concentrations (collected December, 2018-December, 2019). Petroni et al.  investigated the association of COVID-19 mortality with respiratory hazard index calculated across 3223 U.S. counties using emissions data for 2014 and COVID-19 data through May 13, 2020. They observed a 9% increase in COVID-19 mortality per unit increase in respiratory hazard index, which includes DPM. Their analyses with only DPM demonstrated an increased effect of 182% in the mortality rate ratio with a 0.5 μg/m3 increase in DPM concentration. Hendryx and Luo  studied the association of long-term exposure to ozone (obtained from 2016), PM2.5 (obtained from 2016), and DPM (obtained from 2014) with COVID-19 prevalence and mortality through May 31, 2020. They showed an increase of 14.3 deaths per 100,000 U.S. residents for each DPM concentration increase of 1 μg/m3 in a single-pollutant model adjusted for demographic, health, smoking, and COVID-19 testing covariates. These findings collectively suggest that long-term PM exposure may predispose an individual to COVID-19 mortality. However, association of COVID-19 mortality with long-term DPM may change over time with the evolution of the coronavirus and changes in policies and personal behaviors. Our understanding of the effect of long-term DPM exposure on COVID-19 mortality during different waves of the disease and over the locations impacted by those waves remain unknown, hampering anticipation of disease hotspots.
DPM is composed of a complex mixture of black carbon and organic carbon. Studies have shown that 80-90% of particles emitted by diesel engines are smaller than 2.0 μm [12, 24], small enough to penetrate the alveoli . Long-term DPM exposure has been associated with adverse respiratory and cardiovascular effects [12, 37, 41]. Diesel engines power school buses, heavy-duty trucks, a variety of off-road heavy equipment, shipping, and commercial boating [12, 24]. DPM emissions are higher in urban areas, where most of the global population lives [12, 41]. Likewise, greater DPM concentrations have been observed in socioeconomically disadvantaged communities [8, 12].
Our study explores spatial associations between long term average concentrations of DPM, as a metric for past air pollution exposure, and COVID-19 mortality across each pandemic wave and throughout 2020 in the U.S. The objectives of the study are 1) to assess if living near DPM sources increased the risk of death from COVID-19, 2) to estimate how associations between mortality and long-term exposure to DPM (using the U.S. Environmental Protection Agency  broad definition of long-term exposure measured over “months to years”) may have changed over time with changes in the Coronavirus and in the population’s behavior, and 3) to test if models accounting for spatial autocorrelation improve model estimates. Data for air pollution, health, demographic, and social determinants of health were merged for this analysis, and global and local models were both applied to examine these relationships.
Two measures of mortality were considered for our study: mortality rate (defined as the number of dead per 100,000 people in a defined geographic area) and case fatality rate (CFR, defined as the percent dead among confirmed COVID-19 cases in a defined geographic area). County-level number of COVID-19 deaths and CFR were obtained from the publicly-available Johns Hopkins Coronavirus Resource Center  for the period January 1- December 31, 2020. Total population shapefile data were obtained from the U.S. Census Bureau  for the mortality rate calculation. Use of the mortality rate in our model may point to an impact of long-term DPM exposure on COVID-19 mortality in the general population. Use of the CFR in our models may indicate an effect of long-term DPM exposure on mortality among those who are already infected with COVID-19. An advantage of the CFR is that positive associations may indicate that long-term DPM exposure causes death due to COVID-19. However, a disadvantage of the CFR is that two measures (COVID-19 cases and deaths) are estimated, so it is susceptible to errors due to undercounting both the COVID-19 mortality count and the COVID-19 case count. The mortality rate is only susceptible to errors in the death count.
Data for potential confounders associated with the measures of COVID-19 mortality and DPM, including access to health care, education, poverty, demographics, transportation, and occupation were obtained from the American Community Survey (ACS ;) and the County Health Rankings (CHR ;) (Table 1). The variables tested as potential confounders are similar to those used in other studies investigating factors associated with COVID-19 that also tested for confounders and observed associations with variables relating to socioeconomic status, demographics, and healthcare availability that could potentially be correlated with air pollution (e.g., [30, 44]).
Long-term average DPM concentrations were obtained from the 2018 AirToxScreen database, the most recently modeled concentrations of hazardous air pollutants and select other pollutants . EPA used a hybrid model that coupled a Community Multiscale Air Quality (CMAQ) chemical transport model to the American Meteorological Society/Environmental Protection Agency Regulatory Model (AERMOD), a dispersion model, to estimate AirToxScreen air pollutant concentrations at the census tract level through a multi-step process. CMAQ v5.2 was first run over a 12 km × 12 km grid based on DPM emissions inputs from the National Emissions Inventory . Next, the AERMOD dispersion model was run for each source using the same inputs but with receptors distributed over census tract centroids. Finally, concentrations estimated by AERMOD along the census tract centroids were scaled by the ratio of the CMAQ concentration to the average of the AERMOD concentrations over that same grid cell. This formulation allows for more accurate representation of the chemistry and physics of the DPM than the AERMOD dispersion model can provide alone, while maintaining the finer census tract level spatial resolution of the dispersion model. Because the concentrations are calculated from the annual emissions, the concentrations are annual averages.
We tested the associations of COVID-19 mortality rate and CFR with long-term DPM concentrations across the contiguous United States for time periods coinciding with each COVID-19 wave in 2020: January 1-May 31, 2020, June 1-September 30, 2020, and October 1-December 31, 2020. We also ran the models for the entire year: January 1-December 31, 2020.
We used regression analysis to examine spatial non-stationarity in the relationship between the measures of COVID-19 mortality and DPM while accounting for potentially confounding effects. This work is similar to spatial modeling approaches used by Sun et al.  and Rahman et al. . Sun et al.  investigated different spatial regression models and compared them with an ordinary least squares (OLS) regression model to explain the transmission pattern of COVID-19. County-level race/ethnicity and socio-economic covariates were included in their models. We adapted their approach by focusing on associations of COVID-19 mortality with DPM and by investigating different waves of the disease. Three global models, OLS, spatial lag model (SLM), and spatial error model (SEM), were run to produce a nationwide effect estimate. One local model, geographically weighted regression (GWR), produced effect estimates at the county scale. The R Statistical Software version 4.0.5 was used to run all code. We performed spatial regression modeling with the following libraries: spdep, spgwr, and spatialreg.
OLS models are designed to minimize the sum of squared differences between the true data and the prediction across the dataset . Mollalo et al.  studied county-level variations of COVID-19 incidence in the U.S. From a list of 35 demographic, socio-economic, topographic and environmental variables, they used a stepwise forward selection procedure and then checked for multicollinearity to determine the most significant predictors of COVID-19. Then, using the same selected explanatory variables, they tested their model using OLS and several spatial models including SEM, SLM, and GWR (described below). Accounting for spatial autocorrelation in their model improved performance over OLS. Karaye and Horney  also compared OLS to spatial regression models to analyze the impact of social vulnerability on COVID-19 cases. Spatial autocorrelation of the residuals may compromise the validity of the OLS model and produce biased estimators [25, 28]. The model assumptions of zero mean, independence, heteroscedasticity, and normal distribution are met for the case where OLS is a complete and correct model in which the variables capture all of the spatial variation without specifying spatial positions [10, 43]. Spatial autocorrelation in residuals may occur due to an omitted variable. Heteroscedasticity, or dependence of the residuals on the fitted values, may result in part from spatial autocorrelation . This was evaluated in the OLS using the Breusch-Pagan test for heteroscedasticity of the residuals. The SLM and SEM employ generalized frameworks that apply a transformation to the data to improve heteroscedasticity of the data using appropriate control of the error term and calculate efficient maximum-likelihood estimates .
SLMs estimate an autocorrelation parameter (“spatial lag”) using a weighted average of the response variable across neighboring areas, testing if neighboring observations affect one another [16, 26, 47]. As the autocorrelation parameter approaches zero, the SLM approaches the OLS . In SEMs, errors across neighboring areas are autocorrelated (“spatial error”) [16, 23]. SEMs estimate the relationship between the residuals in a spatial region and those in adjacent regions . The spatial structure is in the residuals, meaning that some important predictors are omitted in the model .
SLM and SEM have only one spatial dependence parameter. The single-valued characteristic makes it impossible for global spatial models to reveal local spatial patterns [6, 14]. Another limitation of global spatial models is that the model is dependent on the spatial weighting matrix 6. In contrast, GWR allows for local models to be fit to each observation using spatial distance as a weighting factor for the influence of all other points . To determine local associations between COVID-19 cases in the U.S. and demographic, socio-economic, topographic and environmental parameters, Mollalo et al.  examined two local models including GWR. The variables incorporated in the model are the same set used for OLS, SLM, and SEM. Similarly, Karaye and Horney  compared GWR to OLS to understand the spatially varying effect in the relationship between social vulnerability and COVID-19 case counts. The main advantage of GWR as a local model is the ability to test for spatial variability among the effects of different variables in the model [6, 14, 25]. Another strength is that GWR has the same model structure as the OLS, which facilitates comparison between the two models .
For our spatial autoregressive models, we estimated spatial relationships between regions based on contiguous boundaries shared between 2 or more counties, assuming that COVID-19 spread in a county is influenced by adjacent counties. For GWR, a cross-validation function minimizes the root mean square prediction error that defines the weight matrix. We evaluated spatial autocorrelation among contiguous cells in the model residuals using Moran’s I . Statistically significant Moran’s I indicates either correlation or anticorrelation among neighboring units. Additionally, we used Lagrange multiplier test statistics to understand whether the spatial lag or spatial error pattern is more important for interpreting the local results.
The level of urgency of the COVID-19 outbreak contributed to uncertain policy decisions and interventions in health in compressed timeframes coupled with the complex social, economic and political events of 2020 . Effects related to pandemic waves could have influenced the importance of specific variables during these different times of the year. Therefore, a set of different covariates have been integrated into the model for each time period. To determine which covariates to include in the regression models of COVID-19 mortality, we applied a stepwise selection algorithm for each season (Table 1). Then, the same covariates were incorporated in the best model for OLS, SLM, SEM, and GWR for each specific wave (Table 2), based on the following framework:
The confounder selection procedure was based on minimizing the Akaike information criterion (AIC) after controlling for multicollinearity. We used this same process for each of the three waves and throughout 2020 to find the most significant models for determining the nationwide and local associations between COVID-19 mortality and DPM concentration.
County-level annual average DPM concentration varied from 0.000202 to 1.72 μg/m3 with a nationwide median of 0.204 μg/m3. Elevated DPM concentration could be observed at specific points corresponding to cities (Fig. 1). Across New York metropolitan area counties, which was greatly impacted during the first wave of the pandemic, the average DPM concentration was 0.425 μg/m3.
During the January-to-May wave, the highest cumulative numbers of COVID-19 deaths were found in roughly the same regions as elevated DPM (Fig. 2a). As 2020 progressed, most counties experienced a higher mortality rate. The New York region exhibited lower cumulative mortality rate during the October-to-December wave of our study (Fig. 2c), with a mean of 98 deaths per 100,000 compared with the January-to-May wave, which had a mean of 280 deaths per 100,000 (Fig. 2a). As shown in Fig. 2a and b, cumulative mortality rate increased substantially from the first wave to the second wave in the Southeast region. In the West region, New Mexico, Arizona and California displayed the same pattern as the Southeast, with a large increase during the second wave. For the September-to-December wave, COVID-19 mortality rate increased across almost all of the US, exhibiting nearly the same pattern as for the all-year distribution (Fig. 2c and d).
At a global level, all models demonstrated a statistically significant association between long-term average DPM concentration and COVID-19 mortality rate for the first 9 months of 2020, as represented by the January-to-May and June-to-September waves (Table 3). SLM and SEM produced slightly higher associations for the June-to-September wave. For the October-to-December wave, none of the global models were found to produce positive associations or to be statistically significant. For the entire year, both the OLS and SLM produced positive associations, while the SEM produced a negative association.
OLS did not seem to be the most appropriate model to study spatial associations between COVID-19 mortality rate and DPM. Smaller associations for the spatial autoregression models compared with OLS suggested that the OLS covariates were positively biased due to spatial autocorrelation. Moran’s I and visual inspection of the residuals maps (Supplemental Fig. S1) indicated spatial clusters of high values and of low values. The Breusch-Pagan test provided support for heteroscedasticity of the residuals (p < < 10− 6) in the OLS models, which may have been partially attributed to spatial autocorrelation . The SLM and SEM models provided modest improvements in model fit, as indicated by slightly higher values of coefficient of determination (R2). Model fit testing indicates that the SLM and SEM provided comparable fits, based on the Lagrange multiplier test.
The local spatial differences estimated using the GWR model are presented as a range of values (Table 3). The mean COVID-19 mortality rate – DPM association for the GWR is identical to that of the OLS, but overall R2 for the GWR indicates improved performance over all global models. Spatial distribution of the DPM coefficients indicates changing conditions across the country during the three parts of the year (Fig. 3). During the January-to-May wave, associations were mostly positive across the U.S. (Fig. 3a), up to an increase of 76.94 deaths per 100,000 for every interquartile range (IQR) increase in DPM concentration. During the June-to-September wave, about half of the contiguous US presented a positive association (Fig. 3b), while associations were more negative for the October-to-December wave (Fig. 3c). Year-round COVID-19 mortality rate associations with DPM were similar to those for the October-to-December wave, likely due to the large number of cases during that timeframe. Local variations in R2 across the waves showed high (> 70%) values in the Northeast and Southwest during the January-to-May and June-to-September waves and in the year-long model. High R2 persisted into the October-to-December wave for the Southwest, albeit with a smaller area (Fig. 4). Low values of R2 (< 40%) were observed in the areas with greatest decrease in mortality with increasing DPM concentration, suggesting much greater uncertainty in those associations than in the positive ones seen in the New York area during the first wave. Moreover, COVID-19 mortality rate was statistically significantly associated with DPM concentration during the January-to-May and June-to-September waves but not during the October-to-December wave.
Among all potentially confounding covariates incorporated in the global models, fraction Black race and fraction American Indian ethnicity were statistically significantly positive in all global models. In addition to these two covariates, inactivity was significant in the June-to-September and October-to-December waves and in the year-long model, and the confounders Hispanic, Mining or Agriculture, Public Transportation, Time to Work, Income Inequality, and Population Density were significant at different time periods of the model. A negative relationship was found for smoking for the January-to-May wave, while a strong positive association was obtained for the October-to-December wave. In evaluations of the effects of smoking on COVID-19 incidence or mortality, inconsistent results have been found in the literature [1, 27, 31]. Benowitz et al.  highlighted the need for further investigations to better understand the mechanisms and effect of smoking on COVID-19 related outcomes. Correlation analysis suggests that DPM was moderately correlated with COVID-19 mortality rate in the January-to-May wave (Fig. 5a) and for the year-long (Fig. 5d) model, but that correlation was reduced in the June-to-September (Fig. 5b) and October-to-December (Fig. 5c) waves. Population density was correlated with COVID-19 mortality rate (ρ = 0.74) and DPM (ρ = 0.59) in the January-to-May wave. Fraction using public transportation was moderately correlated with COVID-19 mortality rate (ρ = 0.48) and DPM (ρ = 0.62) in the January-to-May wave. Therefore, population density and public transportation usage had the potential to act as confounders in a model testing the association between COVID-19 mortality and DPM concentration for the January-to-May wave. A model including just DPM in the SLM and SEM produced effect estimates of 7.728 and 9.713 deaths due to COVID-19 per 100,000 people for an IQR change in DPM (with R2 = 0.2 for both models), respectively. Inclusion of the covariates in the model produced effect estimates of 1.710 and 1.845 deaths due to COVID-19 per 100,000 people for an IQR change in DPM (with R2 = 0.42 and 0.45), respectively. These differences suggest that the final models controlled for those confounders. In the model that only included the covariates, the effect estimates for population density and use of public transportation were slightly lower than in the full model, while R2 for the SLM and SEM were the same as for the models including DPM.
Detailed results from the models examining associations between CFR and DPM are included in the Supplemental Material (see Supplemental text, Supplemental Table S1, and Supplemental Figs. S2-S5). The relationships between CFR and DPM were similar to associations between mortality rate and DPM at each wave and throughout 2020. During the January-to-May wave, associations were positive and strongest in the Northeast. Associations were visible through parts of the Midwest and the Pacific Northwest. Associations persisted in the Northeast for the June-to-September and October-to-December waves, but the magnitude of the associations was lower than for January-to-May. Negative associations were observed across the Southern and Mountain states for the June-to-September and October-to-December waves. Model fit (R2) was consistently lower for the CFR models across model type and wave compared with the models testing associations between mortality rate and DPM.
Our study analyzed the spatial correlation of COVID-19 mortality rate and case fatality rate with long-term DPM concentration as a surrogate for exposure across the continental United States during three waves of the COVID-19 pandemic during 2020. Our results suggested that long-term exposure to DPM may have been an important factor in COVID-19 mortality during the first two waves of the disease and that long-term DPM exposure may have been more highly influential during the January-May wave. Sidell et al.  examined associations between air pollution exposure and COVID-19 incidence for monthly and annual averages of PM2.5, nitrogen dioxide (NO2), and ozone (O3) over four waves corresponding to those in our study plus January-February, 2021 for a Southern California cohort. They similarly observed that PM2.5 had a larger effect during the first wave and that the effect diminished over time. A spatial autocorrelation term was controlled for in these models, but Sidell et al.  did not incorporate local methods. Differences in the outcome variable and the specific exposure also necessitates a further examination of spatial and temporal patterns.
Our results indicate that the OLS model does not account for the spatial associations of COVID-19 mortality rate or CFR with DPM concentrations. These results are similar to those of Sidell et al.  and Mollalo et al. , although their studies considered COVID-19 incidence rate rather than mortality. Mollalo et al.  used OLS, SLM, SEM, and two versions of the GWR to model COVID-19 incidence and mortality for the time period of January 22-April 9, 2020 and found notable spatial associations of both COVID-19 incidence and mortality with several predictors. The study of Hendryx and Luo , covering the January-to-May wave, revealed strong associations of COVID-19 prevalence and mortality with long-term DPM and PM2.5 concentrations. Their study estimated a coefficient of 14.3-18.7 deaths per 100,000 U.S. residents for each increase of 1 μg/m3 in DPM concentration. Inflation of the DPM effect shown in their results is possibly due to the correlation between covariates and their mixed linear multiple regression model that does not account for spatial correlation. Stakhovych and Bijmolt  emphasized that correlated spatial errors lead to bias and uncertainty in the OLS results. Moreover, LeSage and Fischer  noted that spatial correlation in the OLS error terms is a sufficient motivation to employ spatial autoregression models for discerning spatial relationships between dependent and independent variables.
The spatial global models outperformed the OLS model fit for all models for both mortality metrics. This improved performance may be related to spatial autocorrelation. A difference in coefficients and R2 among the OLS, SLM, and SEM models was not observed for mortality rate during the June-to-September wave. Kim  reported an inflated effect of spatial autocorrelation on OLS predictor coefficients, suggesting less spatial autocorrelation during the June-to-September wave consistent with Bini et al.  and Smith and Lee .
Among the modeling techniques analyzed for our study, GWR provided the best model fit, based on estimated global R2. Our results revealed where and when local long-term exposure to DPM may have been associated with COVID-19 mortality, consistent with results from both Karaye and Horney  and Mollalo et al.  regarding patterns of local prevalence and local mortality of the disease based on local R2. Some areas in the Northeast and West regions presenting a high R2 in our mortality rate model align with Mollalo et al.  estimates for incidence rate. As noted by Fotheringham et al. , our GWR results illustrate the need to account for local phenomena.
Socio-economic disparity could explain the non-stationary effect of DPM exposure on COVID-19 mortality rate, due to drastic differences between contiguous areas. Socially vulnerable communities, including minoritized racial groups, have seen spatially associated COVID-19 incidence . This is consistent with the strong association we observed for the fraction Black confounder in the mortality rate model (Table 3). Moreover, Paolella et al.  pointed out spatial associations among fine particulate matter concentration, health effects, and minoritized groups and found out that finer spatial resolution revealed substantially higher fine particulate matter concentrations in Black and Hispanic communities.
The differences among associations of COVID-19 mortality rate and DPM concentrations found by the SLM and SEM for the year-long time period, when SLM was demonstrated to be more significant by a Lagrange test, helped to illustrate that neighboring effects were more relevant in modeling the spatial relationship with COVID-19 deaths than unobserved latent variables contained in the error term. Counties near other counties with high COVID-19 incidence are likely to have higher incidence. Nonetheless, since the weighting matrix chosen for our study was based on spatial adjacency, the county size differences between the Eastern and Western U.S. may have affected the parameter estimates creating more uncertainty in the larger counties . Some variability in the association between COVID-19 and DPM exposure within counties might not have been captured, although DPM sources are more likely to be found in urban areas. However, since the SLM and SEM for the year-long time period were not statistically significant, other models should be considered when data are combined across multiple waves.
Several limitations of this study need to be acknowledged with respect to the input data. It is possible that, with more data and/or more time, the associations would disappear. Exposure measurement error could bias the results . Our spatial modeling approach is intended to account for spatial exposure measurement errors. However, errors from applying cross-sectional analyses persist. Although we studied different waves of the disease, our models were not truly longitudinal. Long-term exposure to DPM was estimated using concentrations from 2018. The dataset likely includes higher DPM concentrations than for 2020 given reduced driving patterns during 2020 and, to a lesser extent, fleet turnover. This suggests that the magnitude of the effects of DPM calculated by our study and these other studies were underestimated. Widely reported undercounting of cases and deaths during the January-to-May wave would further contribute to this underestimation .
The set of potential confounders employed in our models was chosen to evaluate the influence of factors other than DPM potentially associated with COVID-19 outcomes . However, it was impossible to represent all influential factors in the relationship between each wave of COVID-19 mortality and long-term DPM concentrations, so uncertainty in the potential for confounding existed [18, 54]. Furthermore, the study was designed at county level. Spatial variation within counties was not captured, and the difference in county size could have caused uncertainty since the weighting matrix defined for our SLM, SEM and GWR accounted for spatial adjacency. Therefore, associations at scales finer than county-level, including individual- and neighborhood-level associations, could not be inferred . Despite these limitations, our study included a rigorous analysis of spatial relationships for different time periods and tested a variety of potential confounders to minimize these limitations.
Our study built on previous findings by exploring associations of COVID-19 mortality rate with long-term DPM concentrations across the first three waves of the pandemic. In doing so, our models provided a picture in which long-term DPM exposure may have influenced COVID-19 mortality during the early stages of the disease, as observed specifically for the periods of January-to-May and June-to-September, 2020. Waning influence of DPM during the October-to-December wave suggested that person-to-person disease transmission regardless of past DPM exposures may have become more influential in the spread of COVID-19 and in mortality rates once the Coronavirus became widespread throughout the U.S. Further investigation might focus on factors associated with COVID-19 mortality rate during the October-to-December wave. Although COVID-19 data were available beyond this period, the introduction of vaccines during 2021 were likely to have been so influential that combination of the 2 years of data may have produced misleading conclusions.
Availability of data and materials
The datasets supporting the conclusions of this article are in the following repositories:
COVID-19 mortality data can be found in the Johns Hopkins University Global Coronavirus (COVID-19) Database, https://data.world/covid-19-data-resource-hub/covid-19-case-counts/workspace/file?filename=COVID-19+Cases.csv.
Diesel particulate matter data can be found in the U.S. Environmental Protection Agency 2014 National Air Toxics Assessment, https://www.epa.gov/national-air-toxics-assessment/2014-nata-assessment-results#self.
Community characteristics data can be found in the U.S. Census Bureau 2014-2018 American Community Survey, https://data.census.gov/cedsci/table?d=ACS%205-Year%20Estimates%20Data%20Profiles&g=0400000US22.050000&tid=ACSDP5Y2018.DP05&hidePreview=false&tp=true&moe=true&vintage=2018
American Community Survey
American Meteorological Society/Environmental Protection Agency Regulatory Model
Akaike information criterion
County Health Rankings
Coronavirus disease 2019
Diesel particulate matter
Environmental Protection Agency
Geographically weighted regression
Intensive care unit
National Air Toxics Assessment
- NO2 :
Ordinary least squares
- O3 :
- PM2.5 :
Particulate matter with a nominal diameter of 2.5 μm
- R2 :
Coefficient of determination
Spatial error model
Spatial lag model
Benowitz NL, Goniewicz ML, Halpern-Felsher B, Krishnan-Sarin S, Ling PM, O'Connor RJ, et al. Tobacco product use and the risks of SARS-CoV-2 infection and COVID-19: current understanding and recommendations for future research. Lancet Respir Med. 2022;10(9):900–15. https://doi.org/10.1016/S2213-2600(22)00182-5.
Bini LM, Diniz-Filho JAF, Rangel TF, Akre TS, Albaladejo RG, Albuquerque FS, et al. Coefficient shifts in geographical ecology: an empirical evaluation of spatial and non-spatial regression. Ecography. 2009;32(2):193–204. https://doi.org/10.1111/j.1600-0587.2009.05717.x.
Bozack A, Pierre S, DeFelice N, Colicino E, Jack D, Chillrud SN, et al. Long-term air pollution exposure and COVID-19 mortality. Am J Respir Crit Care Med. 2022;205(6):651–62. https://doi.org/10.1164/rccm.202104-0845OC.
Carroll RJ, Ruppert D. Transformation and Weighting in Regression, 1st Edition. Boca Raton: Chapman Hall/CRC; 1988.
Centers for Disease Control (CDC). COVID-19; 2020. https://www.cdc.gov/coronavirus/2019-ncov/index.html. Accessed 30 Aug 2021
Chi G, Zhu J. Spatial Regression Models for the Social Sciences. Thousand Oaks: SAGE Publications, Inc.; 2019.
Clark A, Jit M, Warren-Gash C, Guthrie B, Wang HHX, Mercer SW, et al. Global regional and national estimates of the population at increased risk of severe COVID-19 due to underlying health conditions in 2020: a modelling study. Lancet Glob Health. 2020;8(8):e1003–e1017. https://doi.org/10.1016/S2214-109X(20)30264-3.
Clougherty JE, Shmool JL, Kubzansky LD. The role of non-chemical stressors in mediating socioeconomic susceptibility to environmental chemicals. Current Environmental Health Reports. 2014;1:302–13. https://doi.org/10.1007/s40572-014-0031-y.
Davis JA, Meng Q, Sacks JD, Dutton SJ, Wilson WE, Pinto JP. Regional variations in particulate matter composition and the ability of monitoring data to represent population exposures. Sci Total Environ. 2011;409:5129–35. https://doi.org/10.1016/j.scitotenv.2011.08.013.
DeAngelis DL, Yurek S. Spatially explicit modeling in ecology: A review. Ecosystems. 2017;20:284–300. https://doi.org/10.1007/s10021-016-0066-z.
Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. 2020;20(5):533-534. https://doi.org/10.1016/S1473-3099(20)30120-1.
Douglas JA, Archer RS, Alexander SE. Ecological determinants of respiratory health: Examining associations between asthma emergency department visits, diesel particulate matter, and public parks and open space in Los Angeles, California. Prev Med Rep. 2019;14:100855. https://doi.org/10.1016/j.pmedr.2019.100855.
Dubrow JK. Local data and upstream reporting as sources of error in the administrative data undercount of Covid 19. Int J Soc Res Method. 2021. https://doi.org/10.1080/13645579.2021.1909337.
Fotheringham AS, Brunsdon C, Charlton M. Geographically Weighted Regression: The Analysis of Spatially Varying Relationships. New York: Wiley; 2003.
Fotheringham AS, Charlton ME, Brunsdon C. Geographically weighted regression: a natural evolution of the expansion method for spatial data analysis. Environ Plan A. 1998;30:1905–27. https://doi.org/10.1068/a301905.
Gao C, Feng Y, Tong X, Lei Z, Chen S, Zhai S. Modeling urban growth using spatially heterogeneous cellular automata models: Comparison of spatial lag, spatial error and GWR. Comput Environ Urban Syst. 2020;81:101459. https://doi.org/10.1016/j.compenvurbsys.2020.101459.
Goldberger AS. Classical Linear Regression. In: Econometric Theory. New York: Wiley; 1964.
Hendryx M, Luo J. COVID-19 prevalence and fatality rates in association with air pollution emission concentrations and emission sources. Environ Pollut. 2020;265:115126. https://doi.org/10.1016/j.envpol.2020.115126.
Johns Hopkins University. Global Coronavirus (COVID-19) Data; 2020. Accessed 26 Apr 2020. https://data.world/covid-19-data-resource-hub/covid-19-case-counts/workspace/file?filename=COVID-19+Cases.csv
Karaye IM, Horney JA. The impact of social vulnerability on COVID-19 in the US: an analysis of spatially varying relationships. Am J Prev Med. 2020;59:317–25. https://doi.org/10.1016/j.amepre.2020.06.006.
Kim D. Predicting the magnitude of residual spatial autocorrelation in geographical ecology. Ecography. 2021;44:1121–30. https://doi.org/10.1111/ecog.05403.
Lancaster K, Rhodes T, Rosengarten M. Making evidence and policy in public health emergencies: lessons from COVID-19 for adaptive evidence-making and intervention. Evidence Policy. 2020. https://doi.org/10.1332/174426420X15913559981103.
Le Gallo J, Baumont C, Dall’Erba S, Ertur C. On the property of diffusion in the spatial error model. Appl Econ Lett. 2005;12:533–6. https://doi.org/10.1080/13504850500120722.
Lee K-H, Jung H-J, Park D-U, Ryu S-H, Kim B, Ha K-C, et al. Occupational Exposure to Diesel Particulate Matter in Municipal Household Waste Workers. PLoS One. 2015;10(8):e0135229. https://doi.org/10.1371/journal.pone.0135229.
LeSage J, Pace RK. Introduction to Spatial Econometrics. London: Chapman and Hall/CRC; 2009.
LeSage JP, Fischer MM. Spatial growth regressions: model specification, estimation and interpretation. Spat Econ Anal. 2008;3:275–304. https://doi.org/10.1080/17421770802353758.
Liu K, He M, Zhuang Z, He D, Li H. Unexpected positive correlation between human development index and risk of infections and deaths of COVID-19 in Italy. One Health. 2020;10:100174. https://doi.org/10.1016/j.onehlt.2020.100174.
Loonis V, de Bellefon MP. Handbook of Spatial Analysis: Theory and practical application with R: Insee Méthodes; 2018. p. 131.
Mikati I, Benson A, Luben TL, Sacks JD, Richmond-Bryant J. Disparities in distribution of particulate matter emission sources by race and poverty status. Am J Pub Health. 2018;108:480–5. https://doi.org/10.2105/AJPH.2017.304297.
Mollalo A, Vahedi B, Rivera KM. GIS-based spatial modeling of COVID-19 incidence rate in the continental United States. Sci Tot Environ. 2020;138884. https://doi.org/10.1016/j.scitotenv.2020.138884.
Moran PAP. Notes on Continuous Stochastic Phenomena. Biometrika. 1950;37:17–23. https://doi.org/10.2307/2332142.
Paolella DA, Tessum CW, Adams PJ, Apte JS, Chambliss S, Hill J, et al. Effect of model spatial resolution on estimates of fine particulate matter exposure and exposure disparities in the United States. Environ Sci Technol. 2018;5:436–41. https://doi.org/10.1021/acs.estlett.8b00279.
Peek ME, Simons RA, Parker WF, Ansell DA, Rogers SO, Edmonds BT. COVID-19 among African Americans: an action plan for mitigating disparities. Am J Pub Health. 2021;111:286–92. https://doi.org/10.2105/AJPH.2020.305990.
Petrilli CM, Jones SA, Yang J, Rajagopalan H, O’Donnell L, Chernyak Y, et al. Factors associated with hospital admission and critical illness among 5279 people with coronavirus disease 2019 in New York City: prospective cohort study. BMJ. 2020;369. https://doi.org/10.1136/bmj.m1966.
Petroni M, Hill D, Younes L, Barkman L, Howard S, Howell IB, et al. Hazardous air pollutant exposure as a contributing factor to COVID-19 mortality in the United States. Environ Res Lett. 2020;15:0940a9. https://doi.org/10.1088/1748-9326/abaf86/meta.
Phillips N, Park IW, Robinson JR, Jones HP. The Perfect Storm: COVID-19 Health Disparities in US Blacks. J Racial Ethn Health Disparities. 2020;1-8. https://doi.org/10.1007/s40615-020-00871-y.
Pronk A, Coble J, Stewart PA. Occupational exposure to diesel engine exhaust: a literature review. J Expo Sci Environ Epidemiol. 2009;19:443–57. https://doi.org/10.1038/jes.2009.21.
Qian H, Rahbek C, Rodríguez MÁ, Rueda M, Ruggiero A, Sackmann P, et al. Coefficient shifts in geographical ecology: an empirical evaluation of spatial and non-spatial regression. Ecography. 2009;32:193–204. https://doi.org/10.1111/j.1600-0587.2009.05717.x.
Rahman MH, Zafri NM, Ashik FR, Waliullah M. GIS-based spatial modeling to identify factors affecting COVID-19 incidence rates in Bangladesh. medRxiv. 2020. https://doi.org/10.1101/2020.08.16.20175976.
Reyes MV. The disproportional impact of COVID-19 on African Americans. Health and Hum Rights. 2020;22:299.
Ristovski ZD, Miljevic B, Surawski NC, Morawska L, Fong KM, Goh F, et al. Respiratory health effects of diesel particulate matter. Respirology. 2012;17(2):201–12. https://doi.org/10.1111/j.1440-1843.2011.02109.x.
Robert Wood Johnson Foundation. (2020) County Health Rankings. Accessed 25 Aug 2020. http://www.countyhealthrankings.org.
Schabenberger O, Gotway CA. Statistical Methods for Spatial Data Analysis. London: CRC Press; 2017.
Sidell MA, Chen Z, Huang BZ, Chow T, Eckel SP, Martinez MP, et al. Ambient air pollution and COVID-19 incidence during four 2020-2021 case surges. Environ Res. 2022;208:112758. https://doi.org/10.1016/j.envres.2022.112758.
Smith TE, Lee KL. The effects of spatial autoregressive dependencies on inference in ordinary least squares: a geometric approach. J Geogr Syst. 2012;14(1):91–124. https://doi.org/10.1007/s10109-011-0152-x.
Stakhovych S, Bijmolt TH. Specification of spatial models: A simulation study on weights matrices. Papers Regional Sci. 2009;88:389–408. https://doi.org/10.1111/j.1435-5957.2008.00213.x.
Sun F, Matthews SA, Yang TC, Hu MH. A spatial analysis of the COVID-19 period prevalence in US counties through June 28, 2020: Where geography matters? Ann Epidemiol. 2020. https://doi.org/10.1016/j.annepidem.2020.07.014.
U.S. Census Bureau. (2020) 2014-2018 American Community Survey. Accessed 24 Apr 2020. https://data.census.gov/cedsci/table?d=ACS%205-Year%20Estimates%20Data%20Profiles&g=0400000US22.050000&tid=ACSDP5Y2018.DP05&hidePreview=false&tp=true&moe=true&vintage=2018.
U.S. Census Bureau. (2021) Decennial Census (2010) https://www.census.gov/data/developers/data-sets/decennial-census.2010.html#list-tab-40E33DZYH5NNN4DKDV. Accessed 16 July 2021.
U.S. Environmental Protection Agency. Integrated Science Assessment for Particulate Matter; 2019. U.S. Environmental Protection Agency, Research Triangle Park, NC. EPA/600/R-19/188
U.S. Environmental Protection Agency. Technical Support Document. EPA’s Air Toxics Screening Assessment; 2022. 2017 Air ToxScreen TSD. U.S. Environmental Protection Agency, Research Triangle Park. https://www.epa.gov/system/files/documents/2022-03/airtoxscreen_2017tsd.pdf
Villeneuve PJ, Goldberg MS. Methodological considerations For epidemiological studies of air pollution and the SARS and COVID-19 coronavirus outbreaks. Environ Health Perspect. 2020;128(9):095001. https://doi.org/10.1289/EHP7411.
Williamson EJ, Walker AJ, Bhasharan K, Bacon S, Bates C, Morton CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584:430–6. https://doi.org/10.1038/s41586-020-2521-4.
Wu X, Nethery RC, Sabath MB, Braun D, Dominici F. Air pollution and COVID-19 mortality in the United States: Strengths and limitations of an ecological regression analysis. Sci Adv. 2020;6. https://doi.org/10.1126/sciadv.abd4049.
Yancy CW. COVID-19 and African Americans. JAMA. 2020;323(19):1891–2. https://doi.org/10.1001/jama.2020.6548.
We would like to thank anonymous reviewers for taking the time and effort to review the manuscript. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
MEM is grateful for the support of the North Carolina State University Center for Geospatial Analytics. JRB and MEM are grateful for the support of the National Institute for Environmental Health Sciences (P42 ES013638).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Deaths per 100 confirmed cases per change in independent variable. Figure S1. Spatial distribution of residuals for OLS models for mortality rate. Figure S2. Spatial distribution of COVID-19 case fatality rate for (a, top left) January-May, (b, top right) June-September, (c, bottom left) October-December, and (d, bottom right) all of 2020. Figure S3. Map of associations between COVID-19 case fatality rate and long-term DPM concentration for U.S. counties. Figure S4. Spatial distribution of local R2 for the GWR model for case fatality rate. Figure S5. Pearson correlation matrix for case fatality rate for (a, top left) January-May, (b, top right) June-September, (c, bottom left) October-December, and (d, bottom right) all of 2020.
About this article
Cite this article
Mathieu, M.E., Gray, J. & Richmond-Bryant, J. Spatial associations of long-term exposure to diesel particulate matter with seasonal and annual mortality due to COVID-19 in the contiguous United States. BMC Public Health 23, 423 (2023). https://doi.org/10.1186/s12889-023-15064-5