Abstract
Background Limited evidence is available on the health effects of particulate matter (i.e. PM2.5, particulate matter with an aerodynamic diameter < 2.5μm; PM10, < 10μm; PM2.5-10, 2.5-10μm) during the pandemic of COVID-19 in Italy.
Objectives To examine the associations between all-cause mortality and daily PM2.5, PM2.5-10, and PM10 in the pandemic period, and compare them to the normal periods (2015-2019) in Italy.
Methods We collected daily data regarding all-cause (stratified by age and gender), and PM2.5, PM2.5-10, and PM10 for 107 Italian provinces from 1, January 2015 to 31, May 2020. A time-stratified case-cross design with the distributed lag non-linear model was used to examine the association between PM and all-cause mortality during the first three months of the COVID-19 outbreak (March to May in 2020) and the same months in 2015-2019. We also compared the counts and fractions of death attributable to PM in two periods.
Results Overall, Italy saw an increase in daily death counts while slight decreases in PM concentrations in 2020 pandemic period compared to same months of 2015-2019. Mortality effects were significant in lag 0-3 days for PM2.5, lag 0-2 for PM10, and lag 0-1 for PM2.5-10. Each 10 µg/m3 increase in PM was associated much higher increase in daily all-cause mortality during 2020 pandemic period compared to the same months during 2015-2019 [increased mortality rate: 7.24 % (95%CI: 4.84%, 9.70%) versus 1.69% (95%CI: 1.12%, 2.25%) for PM2.5; 3.45 % (95%C: 2.58%, 4.34%) versus 1.11% (95%CI: 0.79%, 1.42%) for PM10, 4.25% (95%CI: 2.99%, 5.52%) versus 1.76% (95%CI: 1.14%, 2.38%) for PM2.5-10]. The counts and fractions of deaths attributable to PM were higher in 2020 than the normal periods for PM2.5 (attributable death counts: 20,062 in 2020 versus 3,927 per year in 2015-2019; attributable fractions: 10.2% versus 2.4%), PM10 (15,112 versus 3,999; 7.7% versus 2.5%), and PM2.5-10 (7,193 versus 2303; 3.7% versus 1.4%).
Conclusions COVID-19 pandemic increased the vulnerability and excess cases of all-cause mortality associated with short-term exposure to PM2.5, PM2.5-10 and PM10 in Italy, despite a decline in air pollution level. This suggests using historical PM-mortality association to calculate health benefits associated with reduction in PMs has big uncertainties.
1. Introduction
Coronavirus Disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus1 has substantially affected human society, including healthcare, economic structure, social relationships. The measures and responses to control virus transmission can protect human health, but it also results in unprecedented side-effects. Even though the severe health impacts of the COVID-19 emergence remain the top priority, it is still unknown how the pandemic may affect the association between environmental exposure and health, notably the health impacts of air pollution, particularly for ambient particulate matter (PM) air pollutant which is an important risk factor for cardiovascular and respiratory health outcomes.2,3
The COVID-19 pandemic has observed a noteworthy decline in anthropogenic PM in many countries, such as China,4,5 Morocco,6 Malaysia,7 India,8 Brazil,9 United States,10 and Spain11 due to the reduced social activities and vehicle exhaust emission.12,13 To control the COVID-19 outbreak, Italian government imposed a quarantine, restricting the migrations of the population and additional lockdown restrictions mandating the temporary closure of non-essential shops and businesses first in ten municipalities of the province of Lodi in Lombardy and neighboring municipalities in the northern region at around February 21, 2020, and then in early-March the quarantine measures were expanded to the entire country.14 Alicandro et al15 suggested that Italy’s first wave of the COVID-19 pandemic has ended in May because no excess mortality but COVID-19 deaths were probably over-registered. During the first wave of COVID-19 pandemic (March to May) in Italy, the severe limitation of people movements determined a significant reduction of PM10, PM2.5 pollutants concentration mainly due to vehicular traffic, especially over hard-hit northern Italy (i.e. Milan16 and Roma17). It is expected that reductions in PM could reduce burden of the air pollution-related diseases, if the magnitude of mortality/morbidity risks associated with PM would not change. However, it is still unclear whether this hypothesis is correct, because the exposure-response association could be changed due to the alteration of people’s behaviors and reorganized medical resources during the COVID-19 pandemic.
Previous assessments have focused on the effects of PM on the death of COVID-19.18-21 However, the number of COVID-19 deaths in official reports may be significantly lower than the actual death toll due to the insufficient testing availability, especially during rapid growth phases of the epidemic.22 In addition, indirect deaths from the pandemic may increase because of the limitation of health care resources, especially for the patients who were scheduled for surgery or treated but suspended due to the pandemic.23,24 On the other hand, some susceptible population are concerned about being infected when going to a hospital.25
The overlap of vulnerable populations related to COVID-19 and unrelated makes the measurement of the pandemic impacts complexed. Unlike the deaths attributed to COVID-19, all-cause mortality could be an accurate measurement to estimate the direct and indirect effects of the pandemic on deaths.26 Therefore, understanding whether the COVID-19 pandemic could switch the association of PM with total all-cause mortality or not would provide important clues regarding health promotion and air pollution control.
In this study, we aimed to quantify the association between short-term exposure to PM and all-cause mortality in Italy during the first wave of COVID-19 pandemic and compare it with that of the same months in 2015-2019 to identify the changes in PM-mortality association.
2. Methods
2.1 Study setting
Italy is a country consisting of a peninsula delimited by the Alps and surrounded by several islands, covering a total area of 301,340 km2. There are currently 107 provinces (second level constituent entities) in Italy, within 20 regions (first-level constituent entities). Its northern regions (Lombardia, Veneto, Piemonte, Emilia Romagna) together host 39% of the national population, and approximately one-half of the Italian GDP is produced there. Such a spatial concentration of economic activities involves the industrial manufacturing sectors to the largest extent, and the consequent high level of emissions is at least in part responsible for heavy pollution in the region.27 During the COVID-19 pandemic, Italy has been one of the worst countries affected by the spreading of coronavirus, especially the northern regions.28 Here, we performed analyses with daily air quality data from national monitoring stations and daily all-cause mortality data at province level.
2.2 Mortality data
We collected daily all-cause mortality from Italian National Institute of Statistics (ISTAT) from January 1, 2015 to May 31, 2020. The mortality data covers 7,357 municipalities in 107 provinces, representing 95.0% of Italy population since January 1, 2015. Daily counts of all-cause mortality were aggregated at province level and stratified by gender and age-specific groups (< 65 years and ≥ 65 years).
2.3 Environmental exposure data
Air pollution data were downloaded from European Environment Agency (EEA) air quality database. The database includes hourly PM2.5, PM10, ozone (O3), nitrogen dioxide (NO2), sulfur dioxide (SO2), and carbon monoxide (CO). The 24-h average concentrations of PM2.5, PM10, NO2, SO2, and CO were calculated as daily concentration, while we used the maximum 8-h average concentrations of O3 as its daily value. Daily PM2.5-10 was calculated as the difference between 24-h average PM10 and 24-h average PM2.5.29 Daily concentrations of air pollutants at station level in the same province were averaged one daily times series. If a province only had one monitoring station, data from this station was used to represent the exposure level of this province.
To allow adjustment for the meteorological factors, we collected the ERA5 hourly surface (at 2 meters above the land surface) ambient temperature and ambient dew point temperature at 0.1°×0.1° spatial resolution from the ERA5-Land hourly data. Hourly data were averaged into daily values. We calibrated the collected temperature data with the observed meteorological data through random forest models (see Supplementary materials for detail), and then we linked the data to the centroid of each municipality based on longitude and latitude. We then calculated daily mean relative humidity from the calibrated ERA5 daily mean temperature and ERA5 daily mean dew point temperature, using the algorithm provided by the “humidity” R package.30 Weather data at municipality level were aggregated into province level by averaging observations of all municipalities within the province.
2.4 Statistical analysis
A time-stratified case-crossover design was used to examine the association between PM air pollution and all-cause mortality. The design compares the air pollution exposure in the case period when events occurred with air pollution exposures in control periods to compare the differences in exposure which might explain the differences in the daily number of cases. In this study, the cases and controls were matched by day of the week in the same calendar month and the same province, to control for the effects of day of the week, seasonality, long-term trend, and spatial variation on deaths or pollution. The Quasi-Poisson regression31 allowing for over-dispersion was applied to perform time-stratified case-crossover design. To determine an appropriate lag time (i.e., the number of days between exposure and the estimated effect) for PM to be used in the main analyses, we compared a variety of lag days and choose all significant lags as the maximum lag day (see Supplementary Figure S1). A linear function was used for PM concentrations while a 3 degrees of freedom natural cubic spline was used for lags. Our initial analyses showed that significant mortality effects were observed in lag 0, 1, 2, and 3 days for PM2.5, lag 0, 1, and 2 days for PM10, and lag 0, 1 days for PM2.5-10. Therefore, we used cumulative effects along lag 0-3 days for PM2.5, lag 0-2 for PM10, and lag 0-1 for PM2.5-10 for subsequent analyses. We have controlled for potential nonlinear and lagged confounding effects of weather conditions, with 3 degrees of freedom natural cubic spline for 21-day moving averages32,33 of daily mean temperature and daily mean relative humidity, respectively.
To compare the associations between PM air pollution and all-cause mortality during COVID pandemic and pre-outbreak periods, we performed above analyses for 2020 COVID pandemic period (from March 1 to May 31, 2020) and the same months during 2015-2019, respectively. The analyses for PM2.5, PM10 and PM2.5-10 were performed separately, to avoid their high collinearity. Fixed effect meta-regression was used to compare the magnitude of the mortality risks associated with PM air pollution in different sub-groups.
To estimate the burden of mortality attributable to PM, the attributable number deaths (AD) caused by PM were calculated every day and total AD was generated by summing the AD during the study period.34 The corresponding attributable fractions (AF) of mortality were calculated by dividing the total AD by the death toll.
2.5 Sensitivity analyses
Sensitivity analyses were performed to examine the robustness of the results. We tested the variation of the PM pollution-mortality association in normal period by replacing the study period of 2015-2019 with every single year. To evaluate potential impacts of gaseous pollutants on the associations between PM air pollution and mortality, we also performed multi-pollutant models through adjusting for different combinations of NO2, CO, O3, and SO2.
All the analyses were performed by the R software (v. 3.6.1). The ‘dlnm’ package was used to perform the distributed lag non-linear models for PM, describing simultaneously the linear relationship along air pollution and non-linear along lags; the ‘gnm’ package was used to perform conditional Poisson regression.35 The “mvmeta” package was used to perform meta-regression.36 The relative risks (RRs) with 95% confidence intervals (CIs) per 10 µg/m3 change in PM concentration were reported. For all statistical tests, a p value of 0.05 (two-tailed) was considered statistically significant.
3. Results
There were 1,000,459 (51.9% females; 89.4% aged ≥ 65 years) all-cause deaths in the years from March to May during 2015 and 2020. The death counts in the pandemic period were significantly higher than the same months during 2015-2019. For instance, the daily average mortality counts have risen by 24.61% from 837 ± 62 in 2015-2019 to 1043 ± 325 in 2020 for men and by 19.98% for women from 911 ± 80 to 1093 ± 271 (Table 1). There were no significant trends in mean temperature and relative humidity from 2015 to 2020, but the pandemic period witnessed a slight decrease in relative humidity versus an increase in mean temperature compared with the same months during 2015-2019 (Table 1).
There were slight reductions in PM10, PM2.5 and PM2.5-10 concentrations over 2015-2020 (Table 1). The average (± SD) PM2.5 concentration during March and May reduced from 15.60 ± 10.23 μg/m3 in 2015 to 12.52 ± 7.57 μg/m3 in 2020, PM10 from 23.05 ± 12.25 μg/m3 to 20.54 ± 13.95 μg/m3, and PM2.5-10 from 8.86 ± 5.83 μg/m3 to 8.14 ± 7.43 μg/m3,with all the differences being statistically significant (p < 0.001).
Figure 1 shows the spatial variation in all-cause deaths, PM2.5, PM10 and PM2.5-10 during March and May in 2020 and 2015-2019, respectively. The figure presented the difference by subtracting the average of 2015-2019 from the daily value in 2020. The northern region had increased number of deaths in 2020. The PM concentrations declined slightly from 2015-2019 to 2020.
Cumulative RRs along lag 0-3 days for PM2.5, lag 0-2 for PM10, and lag 0-1 for PM2.5-10 for all-cause mortality and group-specific mortality are shown in Figure 2. Each 10 µg/m3 increase in PM was associated much higher increase in daily all-cause mortality during 2020 pandemic period compared to the same months during 2015-2019 [increased mortality risk: 7.24 % (95%CI: 4.84, 9.70) versus 1.69% (95%CI: 1.12, 2.25) for PM2.5; 3.45 % (95%C: 2.58, 4.34) versus 1.11% (95%CI: 0.79, 1.42) for PM10, 4.25 % (95%CI: 2.99, 5.52) versus 1.76% (95%CI: 1.14, 2.38) for PM2.5-10]. Such disparity in the PM-mortality associations were consistent among different gender and age groups.
Table 2 shows the attributable mortality fractions and attributable deaths associated with PM2.5, PM10 and PM2.5-10 during March and May in 2020 and average values in 2015-2019. AFs and ADs were higher in 2020 than 2015-2019. We estimated that 10.21% (95% CI: 7.13, 13.31) of deaths were attributable to PM2.5 in first three months of pandemic in 2020, whereas the average AF was only 2.44 (95% CI: 1.63, 3.23) in 2015-2019, such disparity was consistent across all sex and age groups and was similar for PM2.5-10 and PM10. During the March to May, 20,062, 15,112, and 7,193 all-cause deaths were estimated to be attributable to PM2.5, PM10 and PM2.5-10, in 2020, approximately 5 times higher than the average values in 2015-2019.
Sensitivity results in Figure 3 consistently showed a stronger PM-mortality association during 2020 pandemic period compared to the associations during the same months in each year of 2015-2019, despite the variations after adjusting for different gaseous pollutants. Likewise, when we used different lag days of PM (lag 0-4 and lag 0-5 for PM2.5, lag 0-3 and lag 0-4 for PM10, lag 0-2 and lag 0-3 for PM2.5-10), the effects remained higher during 2020 pandemic period than the same period in 2015-2019 (Figure S2).
4. Discussion
To our best knowledge, this the first study to investigate the relationship between PM air pollution and daily all-cause mortality during the COVID-19 pandemic period in the world. In this study, we examined the effects of PM (PM2.5, PM10, and PM2.5-10) on daily all-cause mortality in 107 Italian provinces and compared the mortality risks and mortality burdens associated with PM before and during the COVID-19 pandemic. The mortality risks (vulnerability) and burden associated with PM2.5, PM10 and PM2.5-10 in the pandemic were significantly higher than risks estimated in 2015-2019. People aged ≥ 65 years were consistently at higher risk than younger people in both pandemic period and normal periods.
Exposure to PM air pollution has been identified as the risk factor for excess mortality.37,38 Our results are in line with a previous multi-country epidemiological study,39 in which it observed a 0.65% (95% CI: 0.26%, 1.04%) increase in all-cause mortality risk per 10 μg/m3 increase at lag 0-1 of PM10 in 18 Italian cities during 2006-2015. Although the effect estimate is slightly lower than ours for PM10 at lag 0-2 [increased risk:1.11% (95% CI: 0.79%, 1.42%)], such difference in the PM10-mortality associations could be ascribed to the heterogeneity in period and season, as PM-mortality association might be stronger during cold months that we chose.40
The plausible pathways by which ambient PM air pollution may impact increased mortality risks have been well documented. Biologically, particles less than 10 μm in diameter41 have been acknowledged to be responsible for oxidative stress and chronic inflammation status that causes the hyper-activation of the immune system and the life-threatening respiratory disorders.42,43 The rapid and direct passage of smaller sized particles such as PM2.5 even escape into the blood stream and directly affect other organs far away from the lungs through systemic circulation.44 In addition, another plausible pathway for these adverse impacts is that PM may act as both a carrier and substrate of other toxic substance (e.g. carbonaceous elements,45 hazardous metals,46 etc.) and thus influence their fate and transport in the environment and reaching susceptible receptors.
In this study, we found that PM-attributed deaths during COVID-19 pandemic in 2020 were 5 times higher than the same months during 2015-2019, despite a lower PM level. One potential explanation is that PM may contribute to the COVID-19 related deaths. Most current studies have reported a positive association between ambient PM2.5 or PM10 and COVID-19 deaths,27,47-53 although they were limited by the inaccurate official reports of COVID-19 deaths. A review study highlighted the potential role of PM in the spread of COVID-19, focusing on Italian cities in which correspondence between poor air quality and COVID-19 induced mortality was the starkest yet.54 First, COVID-19 could have an air transmission55,56 and atmospheric PM could create a suitable environment for transporting the virus at greater distances than those considered for close contact.57 Second, PM has been shown to induce inflammation in lung cells58 and exposure to PM could increase the susceptibility and severity of the COVID-19 patient symptoms.59
However, the increase of PM-related mortality during pandemic period could also be from non-COVID-19 causes, particularly for those vulnerable populations (e.g., elderly people, people with chronic diseases). For example, during the early days of the COVID-19 outbreak, the regional disparities in health-care resource availability and accessibility60 could play an important role in the change of the vulnerability to PM related deaths. In March, though Italy’s health system has 3.2 hospital beds per 1000 people (as compared with 2.8 in the United States), it was over occupied to meet the needs of rapid surge of COVID-19 cases.61 The health needs created by the coronavirus pandemic went well beyond the capacity of national health system, diagnostic, therapeutic, and preventive interventions were scarce and rationed.62,63 In the worst scenario, patients with PM-related diseases would die while waiting for needed resources (e.g. ICU, ventilator and acute care).64 This hypothesis is supported by observations showing a significant decrease in hospitalization rates for acute coronary syndrome (13.3 admissions per day versus 18.0)65 and acute myocardial infarction (a reduction of 52.1% in North Italy, 59.3% in Central Italy, and 52.1% in South Italy)66 compared with the equivalent time in 2019 in Italy.
We also investigated the different effects between gender and age groups in the pandemic and normal periods. With regards to the age-stratified analysis, older people generally have a higher risk of PM-related death risk in both pandemic period and pre-pandemic periods, which was consistent with most previous studies.67,68 Though the detailed reasons are still unclear, pre-existing cardiovascular and respiratory diseases are more prevalent in elderly people, and it may enable the elderly to be more susceptible to ambient PM air pollution.
There are some limitations to this study. First, like most case-crossover or time-series studies, we used province-level air pollution to represent the individual-level exposure, which is likely to cause random exposure assessment error and thus underestimate the PM-mortality associations. However, since the same design was applied to the 2020 and 2015-2019 period, this error is not likely to affect our main findings. Second, we were not able to assess the actual indoor exposures to PM, which might be important because people would increase the indoor time during the lockdown period. Finally, we cannot exclude the COVID-19 deaths from our analyses because the daily data were not available at province level. Therefore, we were uncertain about whether the increase in PM-mortality association was due to COVID-19 deaths or deaths due to other causes.
In conclusion, with a large nationwide data set covering 107 Italian provinces, we observed significantly increased impacts of PM on all-cause mortality during the pandemic period compared to pre-COVID-19 periods. This suggests the historical exposure-response relationship between PM and mortality may underestimate the health impacts of PM during the COVID-19 pandemic, although air pollution concentrations declined. To date, the second wave of pandemic is ongoing in many Europe countries and no sign shows the end of the COVID-19 pandemic worldwide. Therefore, we would expect our findings to be helpful to guide future public health policy and environmental protection strategies such as making early warning systems for highly polluted regions, and giving close attention to vulnerable groups.
Data Availability
Original data are open sources.
Contributors
TY, RX, WY, YG, and SL designed the study. TY, RX, WY, and ZC collected data. TY did statistical analyses. TY and WY wrote the manuscript. RX, YG and SL provided scientific comments on the manuscript. All authors provided critical conceptual input, analyzed and interpreted data, and critically revised the report.
Funding
No specific funding for this study. TY, and RX were supported by China Scholarship Council funds (number 201906320051 for TY and 201806010405 for RX); SL by an Early Career Fellowship of the Australian National Health and Medical Research Council (number APP1109193); and YG by Career Development Fellowships of the Australian National Health and Medical Research Council (numbers APP1107107 and APP1163693)
Competing interests
The authors declare that they have no competing interests.