Next Article in Journal
Racial and Sex Differences between Urinary Phthalates and Metabolic Syndrome among U.S. Adults: NHANES 2005–2014
Previous Article in Journal
Precautionary Behaviors during the Second and Third Phases of the COVID-19 Pandemic: Comparative Study in the Latin American Population
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Country-Level Excess All-Cause Mortality and the Impacts of Air Pollution and Human Activity during the COVID-19 Epidemic

1
Department of Land Surveying and Geo-Informatics, The Hong Kong Polytechnic University, Hong Kong
2
Research Institute for Sustainable Urban Development, The Hong Kong Polytechnic University, Hong Kong
3
School of Geography, South China Normal University, Guangzhou 510000, China
4
College of Geography and Environment, Shandong Normal University, Jinan 250000, China
5
Department of Geography and Resource Management, The Chinese University of Hong Kong, Hong Kong
6
Institute of Space and Earth Information Science, The Chinese University of Hong Kong, Hong Kong
7
Department of Human Geography and Spatial Planning, Utrecht University, 3584 CB Utrecht, The Netherlands
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(13), 6883; https://doi.org/10.3390/ijerph18136883
Submission received: 17 May 2021 / Revised: 22 June 2021 / Accepted: 23 June 2021 / Published: 26 June 2021

Abstract

:
The impact of Coronavirus Disease 2019 (COVID-19) on cause-specific mortality has been investigated on a global scale. However, less is known about the excess all-cause mortality and air pollution-human activity responses. This study estimated the weekly excess all-cause mortality during COVID-19 and evaluated the impacts of air pollution and human activities on mortality variations during the 10th to 52nd weeks of 2020 among sixteen countries. A SARIMA model was adopted to estimate the mortality benchmark based on short-term mortality during 2015–2019 and calculate excess mortality. A quasi-likelihood Poisson-based GAM model was further applied for air pollution/human activity response evaluation, namely ground-level NO2 and PM2.5 and the visit frequencies of parks and workplaces. The findings showed that, compared with COVID-19 mortality (i.e., cause-specific mortality), excess all-cause mortality changed from −26.52% to 373.60% during the 10th to 52nd weeks across the sixteen countries examined, revealing higher excess all-cause mortality than COVID-19 mortality in most countries. For the impact of air pollution and human activities, the average country-level relative risk showed that one unit increase in weekly NO2, PM2.5, park visits and workplace visits was associated with approximately 1.54% increase and 0.19%, 0.23%, and 0.23% decrease in excess all-cause mortality, respectively. Moreover, compared with the impact on COVID-19 mortality, the relative risks of weekly NO2 and PM2.5 were lower, and the relative risks of weekly park and workplace visits were higher for excess all-cause mortality. These results suggest that the estimation based on excess all-cause mortality reduced the potential impact of air pollution and enhanced the influence of human activities compared with the estimation based on COVID-19 mortality.

1. Introduction

Coronavirus Disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was first identified in December 2019, and has caused 2,695,004 deaths worldwide as of 18 March 2021. Countries around the world proposed many policies to mitigate the spread of the disease, including lockdowns and the restriction of human mobility [1,2]. However, the lockdown policies themselves may have triggered other health issues. Research revealed that the restriction of social contact can cause severe mental health issues [3,4,5] and physical problems such as obesity [6]. In fact, there were death cases around the world from people who were at high risk of other diseases (e.g., cancer) who suffered from COVID-19. The data may have led to biases in COVID-19 death rate evaluations—which have been defined herein as excess all-cause mortality [7]. The reason for the present work’s focus on the excess all-cause mortality was to assess the variation between the number of deaths that occurred and the number of expected deaths in the absence of COVID-19 [8,9]. When compared with the death counts of COVID-19 (i.e., cause-specific mortality), excess all-cause mortality should depict the influence and impact of the COVID-19 epidemic with lower biases [10].
Excess all-cause mortality during COVID-19 was elucidated by researchers in many countries. For instance, Modi, et al. [11] indicated that excess mortality was approximately two times higher than COVID-19 deaths up until April 2020. Similar findings were found in the United States during the first three months of COVID-19, suggesting overstated excess mortality, compared with the official COVID-19 mortality [12]. In Germany, the excess mortality rate varied temporally in different weeks [13]. Moreover, a significant association was found between excess mortality and older-age groups [14,15]. The differences in excess mortality were also correlated by gender, with more men than women dying from COVID-19 [16,17]. Considering its spatial and temporal variation in different countries, it was deemed necessary to collate the excess all-cause mortality with fine temporal scale—and on a global scale.
Current research estimating the excess all-cause mortality associated with COVID-19 can be divided into two approaches. The first approach consists of calculating the average death rate over the past several years [8,18]. The second approach involves estimating excess mortality based on statistical models. For instance, Gibertoni et al. [15] estimated the number of observed deaths using time as an independent variable based on a linear regression model. General additional models using weeks of the year as indicator functions have been adopted to establish a semiparametric model for excess mortality estimation [12]. Moreover, spatiotemporal trends have been considered as a method to form a simulation-based approach for mortality count prediction [19]. Researchers have also used various transmission suppression levels and different relative risks to estimate excess mortality based on different COVID-19 incidence scenarios [20]. However, differences in estimation accuracy among these approaches remain unclear. Furthermore, to what extent the excess all-cause mortality was influenced by the COVID-19 outbreak on a global scale is also unknown.
Research has suggested that rising COVID-19 cases vary not only by demographic characteristics but also by environmental and human activity-related conditions [21,22,23,24]. Air pollutants, including nitrogen dioxide (NO2) and fine particulate matter with a diameter less than 2.5 μm (PM2.5), have been shown to be ambient anthropogenic emissions associated with acute respiratory infections and asthma cases [25,26]. Researchers evaluated the potential impact of decreasing NO2 and PM2.5 levels—which were associated with declining COVID-19 mortality trends [27,28]. The results revealed that chronic exposure to air pollutants delays the recovery period after COVID-19 infection and leads to more severe conditions [29,30]. Meanwhile, regular human activity patterns, including work and leisure, were significantly changed due to restrictions on social contact. Those restrictions, in turn, influenced the transmission rates of COVID-19 [31,32]. Research also suggested that green spaces facilitate social distancing, and could potentially mitigate the spread of COVID-19 [33]. In addition, growing trends related to working from home could be assessed in view of the incidence of COVID-19 [34]. Despite the fact that significant research has been done on the environmental and socioeconomic responses to COVID-19 incidence and mortality, their potential impact on the excess all-cause mortality required deeper investigation.
Faced with these challenges, this study estimated the excess all-cause mortality during COVID-19 and further investigated the potential impact of air pollution and human activity in 16 countries. Specifically, weekly excess deaths were estimated and compared by proposing three-year average mortality and Seasonal Autoregressive Integrated Moving Average (SARIMA) model-based mortality estimations. These were further evaluated by investigating the potential influence of exposure to NO2 and PM2.5 and the frequency with which people visited parks and workplaces. This study could help policymakers and stakeholders to better understand the COVID-19 pandemic under high-risk conditions.

2. Materials and Methods

2.1. Study Area and Datasets

This study investigated the excess all-cause mortality and the influence of air pollution and human activity based on the country level, and selected sixteen countries as its study area, including: Belgium, Chile, Croatia, the Czech Republic, France, Germany, Hungary, Lithuania, Luxembourg, Norway, Poland, Spain, Switzerland, and the United States.
We collected weekly country-wide all-cause mortality from the Short-Term Mortality Fluctuations data (STMF) in the Human Mortality Database (HMD) from the years 2015 to 2020. While weekly mortality rates from the 10th through the 52nd week during COVID-19 in 2020 were considered observed data, the mortality data from 2015 to 2019 were utilized to estimate the benchmark of weekly mortality without COVID-19 in 2020. On this basis, the excess all-cause mortality was considered to be the difference between the observed and benchmark mortality rates in 2020. Accordingly, the following quantification of the COVID-19 mortality, air pollution and human activities were all within the same time period: from the 10th through 52nd week, based on one-week intervals.
The COVID-19 mortality data were collected from the Johns Hopkins University Center for Systems Science and Engineering (JHU CSSE). Daily mortality data were summarized into weekly data to meet the time consistency of the weekly all-cause mortality.
The air pollution—including NO2 and PM2.5—was quantified using ground station data obtained from the OpenAQ Platform. The independent ground station data were aggregated into country-level data based on the population-weighted means method, in which global population data in 2019 from LandScan were utilized within the 10km buffers of air stations. The obtained air pollutants were first aggregated into daily data by filling missing values based on KalmanSmoother [35] and were subsequently summed into weekly data.
We focused on two aspects of human activity changes: business and leisure. In this study, we adapted data on the frequencies of daily visits to parks and workplaces from the Google community mobility report to depict human activity patterns. The quantification of the visit frequency values was presented as positive and negative percentages compared with the baseline day, calculated with the median value from a 5-week period during January and February in 2020. Daily human activities were aggregated into weekly data by summing the frequencies of daily visits to parks and workplaces.
In addition, we adapted meteorological data from the Global Surface Summary of the Day (GSOD) from the Integrated Surface Hourly (ISH) dataset, including daily mean temperature and total precipitation, as controlling data. Weekly data were aggregated based on the average daily temperature and the total daily precipitation in each week.

2.2. Excess Mortality Estimation

Excess mortality was defined as the difference between observed and benchmark death counts. While the observed mortality was the actual death count during a time period, the benchmark mortality was defined as the estimated death count according to historical trends, excluding the influence of unexpected disease (i.e., COVID-19).
To estimate the excess all-cause mortality during COVID-19 in each country, this study adopted two approaches to quantifying the benchmark mortality trends from the 10th to the 52nd week of 2020. The first approach was to estimate the benchmark death count based on three-year average mortality rates. Specifically, the average weekly death count within the predefined time period was calculated, and this was considered the benchmark based on the counterfactual condition of COVID-19.
The second approach was to predict the benchmark death count using a Seasonal Autoregressive Integrated Moving Average (SARIMA) model based on historical mortality trends from 2015 to 2019. The SARIMA model is presented as follows:
S A R I M A ( p , d , q ) ( P , D , Q ) t  
where p , d , q ,   P , D , Q indicate, respectively: autoregressive order, difference order, and moving average order, seasonal autoregressive order, difference order, and moving average order. The proposed SARIMA models—with varied parameters—were estimated using the Kwiatkowski–Phillips–Schmidt–Shin (KPSS) tests, Bayesian information criterion (BIC) and the Ljung–Box test—the last of which was used to test the randomness of the temporal data. BIC was applied to estimate parameters (including p , q ,   P , Q ) for better model performance. KPSS was performed to estimate the parameters d and D. It should be noted that the estimated values based on the SARIMA model were converted into death counts, with negative values converted into 0.
To evaluate the performance of the above two approaches among 16 selected countries, two indicators, including average error and the percentage of average errors lower than 0.1, were proposed. The average error in the ith country is calculated as follows:
P r e d i c t e d   e r r o r i , t = A c t u a l   d e a t h i , t P r e d i c t e d   d e a t h i , t A c t u a l   d e a t h i , t
A v e r a g e   e r r o r i = | t = 1 n P r e d i c t e d   e r r o r i , t / n |  
in which P r e d i c t e d   e r r o r i , t refers to the estimated error of death count in the ith country during the tth week. A c t u a l   d e a t h i , t and P r e d i c t e d   d e a t h i , t represent the actual and estimated death count in the ith country during the tth week. n is the total number of weeks. The percentage of average errors lower than 0.1 is calculated as:
P e r c e n t a g e i = n ( | P r e d i c t e d   e r r o r i , t | < 0.1 ) / n .
where n indicates the number of weeks with predicted errors lower than 0.1. Higher percentages are associated with higher accuracy.
Considering the COVID-19 pandemic in 2020, the performances of both the three-year average mortality model and the SARIMA model were evaluated based on weekly mortality from 2015–2019, with data in 2015–2018 used for modelling and data in 2019 used for model assessment. Consequently, the predicted death count for 2020 was considered the benchmark and used to estimate differences with actual mortality in the same year.

2.3. Association Analysis

The impacts of air pollution and human activity on excess all-cause mortality in each country were estimated based on the quasi-likelihood Poisson-based GAM model. In particular, air pollution and human activities were quantified by weekly ground-based NO2 and PM2.5 levels and the frequency of visits to parks and workplaces. These were additionally controlled by weekly precipitation and temperature. Lag effects of 0, 1 and 2 weeks were used in the regression model to reduce the biases caused by time lagging. On this basis, interactions between weekly mortality and air pollution/human activity impacts were proposed to distinguish the differences between excess all-cause mortality and COVID-19 mortality in each country, modeled as follows:
log E ( Y i ) = α + β 1 n o 2 i , t + β 2 p m 2.5 i , t + β 3 p a r k i , t + β 4 w o r k i , t + s ( p r e c i , t ) + s ( t e m p i , t )
where E ( Y i ) refers to the expected weekly death counts in the ith country. Variables n o 2 i , t , p m 2.5 i , t , p a r k i , t and w o r k i , t represent weekly ground-based NO2 and PM2.5 levels and the frequencies of visits to parks and workplaces in the ith country during the tth week, respectively. α indicates the intercept and β 1 , β 2 , β 3 and β 4  represent the coefficients of the corresponding variables. s() is the smoother, based on the penalized smoothing spline, with the number of knots being determined by generalized cross-validation (GCV).

3. Results

3.1. Excess All-Cause Mortality during COVID-19

The performances of benchmark mortality based on the three-year average estimation and the SARIMA-based estimation are displayed in Figure 1. Average errors in most countries (based on the three-year average estimation) were below 0.02, while the average errors in SARIMA-based estimations were around 0.01. The percentages of average errors lower than 0.1 were approximately 0.9 for the three-year average estimation, while the percentages for SARIMA-based estimations were higher in most countries. Thus, we concluded that SARIMA-based estimation performed better than three-year average estimation in benchmark mortality estimation. The excess all-cause mortality during COVID-19 was further calculated based on the difference between the observed all-cause mortality and the benchmark all-cause mortality estimations based on the SARIMA model.
Figure 2 shows the observed and benchmark all-cause mortality rates over the 10th–52nd weeks of 2020 in 16 studied countries. Significant peaks of increasing death count in observed mortality rates were noted in comparison with the benchmark mortality in many countries. For instance, drastically increased mortality rates were seen in the differences between observed and benchmark data during the 20th to 30th weeks in Chile. Meanwhile, several peaks were observed during the 10th to 20th weeks, 30th to 35th weeks, and 40th to 52nd weeks in Belgium. On the other hand, continuously-increasing observed mortality rates (compared with the benchmark) were observed through the 10th–52nd weeks in the United States, indicating the large number of deaths caused by the COVID-19 outbreaks there.
Figure 3 shows temporal weekly mortality trends during the 10th–52nd weeks in 16 countries, in order to further distinguish the excess all-cause mortality. Diverse temporal mortality patterns were observed in different countries. In 12 countries—Chile, Croatia, the Czech Republic, France, Hungary, Lithuania, Luxembourg, Netherlands, Poland, Spain, Switzerland and the United States—the variations observed in excess all-cause mortality were consistent with those of COVID-19 mortality. Essentially, the death count of excess all-cause mortality was significantly larger than COVID-19 mortality. This indicated that the outbreaks of COVID-19 had a drastic influence on increasing all-cause mortality. On the other hand, additional peaks of increasing death count for excess all-cause mortality were noted in four countries: Belgium, Germany, Norway and Portugal. The inconsistency between all-cause and COVID-19 mortality can be explained by the potential impact of COVID-19. For instance, studies have reported an increased risk of death in hospitalized COVID-19 patients with cancers, fever, or acute respiratory distress syndrome [36,37,38,39].
Table 1 summarizes the COVID-19 and excess all-cause death counts and the percentages of changes for all 16 countries. Chile and Luxembourg showed 12.94% and 26.52% decreases from COVID-19 to excess mortalities, respectively. In all remaining countries, increasing death counts were observed, ranging from 1.81% to 373.60%. This indicated the significant discrepancy between COVID-19 and excess mortalities, suggesting the need to investigate excess all-cause mortality during COVID-19 and the potential impacts of environmental and socioeconomic conditions in facilitating control of the virus.

3.2. Air Pollution/Human Activity Impacts Analysis

Table 2 summarizes the minimum, maximum, mean, and standard deviation values of air pollution data, human activity data, meteorological data, and excess all-cause mortality. On this basis, the impacts of air pollution and human activity were estimated on weekly excess all-cause mortality and COVID-19 mortality using the quasi-likelihood Poisson-based GAM model. The performances of GAM with lag effects of 0, 1, and 2 weeks were evaluated based on R square values. As shown in Figure 4, the performances of models based on COVID-19 mortality were relatively higher than those based on excess all-cause mortality. Moreover, no significances were observed among the models with different lag effects, with most values within 0.5–0.9. Model performances also varied among different countries; Croatia obtained relatively higher R square values than others. Considering that the R square values of most countries were higher than 0.5, the performances of models based on both excess and COVID-19 mortality were ensured.
To distinguish the impacts of air pollution and human activities on both excess all-cause and COVID-19 mortalities, relative risks were calculated for each factor with a 95% confidence interval. Figure 5 shows the relative risks of weekly NO2 on the death count variation with lag effects of 0,1, and 2 weeks. For the estimation based on the excess all-cause mortality, relative risks varied from 0.99 to 1.05 in most countries. Lag effects showed less impact on the relative risk estimation of NO2 in most countries. However, the lag effects (1 or 2 weeks) in Belgium indicated a decreasing trend of NO2 on the increasing excess mortality, while less impact was observed with no lag effect. For the estimation based on COVID-19 mortality, the relative risks of NO2 were approximately within 1 to 1.08, suggesting that increasing NO2 levels corresponded to increasing COVID-19 mortality in most countries. Similar to the excess mortality-based estimation, the relative risks in Belgium also revealed negative patterns with lag effects. Compared with the COVID-19 mortality estimation, the relative risks of NO2 estimated by excess mortality were lower, suggesting less impact of NO2 when estimating potential environmental factors. On the other hand, despite the lockdown policies—which significantly reduced air pollutant emissions—the death counts of both excess and COVID-19 mortality were still positively related to ground-level NO2 emissions.
Figure 6 shows the relative risks of PM2.5 based on excess all-cause and COVID-19 mortalities. While the ranges of relative risks were approximately within 0.96 to 1.04 (based on excess mortality), the relative risks based on COVID-19 narrowed down to 0.97–1.02. This indicated that the impact of PM2.5 on COVID-19 mortality has been overestimated compared with excess mortality. For instance, in Lithuania, Luxembourg, Portugal and the United States, no significant relative risks were observed with a 2-week lag effect based on excess mortality, whereas positive and negative relative risks from 0.99 to 1.01 were estimated in those countries based on COVID-19 mortality with a 2-week lag effect. Overall, lower levels of PM2.5 corresponded to increasing death counts in most countries. Additionally, 1- or 2-week lag effects increased these potential risks.
Figure 7 displays the influence of the frequencies of park visits on mortality variation. In particular, the relative risks in Chile showed significant differences from those in other countries. While relative risks—ranging from 0.995 to 1.01, based on COVID-19 mortality with no lag effect—revealed an overall positive trend between park visits and death counts, lower frequencies of park visits were related to increasing death counts when estimated based on excess mortality with lag effects of 0,1, and 2 weeks and when estimated based on COVID-19 mortality with lag effects of 1 and 2 weeks. A negative association was also revealed in the relationship between park visits and COVID-19 mortality in Lithuania. For other countries, higher relative risks of park visits were related to excess all-cause mortality with negative association—compared with those based on COVID-19 mortality with lag effects of 0,1, and 2 weeks. Specifically, a 1-unit decrease in park visits corresponded to approximately 1% increase in death count for excess mortality.
Figure 8 reveals the relative risks of workplace visitation to the excess all-cause and COVID-19 mortalities. Similar to park visits, the relative risks based on excess all-cause mortality were higher than those based on COVID-19 mortality. In terms of the excess all-cause mortality estimation, the relative risks in Chile showed a positive association between workplace visits and death count, ranging from 1 to 1.03. On the other hand, relative risks in other countries ranged from approximately 0.99 to 1 with lag effects of 0,1, and 2 weeks. This indicated a negative association between workplace visits and death counts in most countries. For the estimation based on COVID-19 mortality, a lesser influence was noted for workplace visits (approximately 0.995 to 1.005). The relative risks based on COVID-19 mortality also varied among studied countries, changing between negative and positive trends with lag effects in several countries, including Poland, Spain, and Luxembourg. On the other hand, consistent positive trends of relative risks were shown among different lag effects in some countries, including Chile and Portugal.
Table 3 displays the average country-level relative risks (without 95% confidence interval) of weekly NO2, PM2.5, park visits, and workplace visits on COVID-19 and excess mortality, with lag effects of 0, 1, and 2 weeks. Generally, positive relationships were revealed between weekly NO2 and COVID-19/excess mortality, while negative associations were shown between weekly PM2.5, park visits, workplace visits and COVID-19/excess mortality with lag effects of 0, 1, and 2 weeks. Mean relative risks of four variables (calculated based on the average relative risk with lag effects of 0, 1, and 2 weeks) showed values of 1.0171, 0.998, 0.9981, and 0.9979, associated with COVID-19 mortality, and values of 1.0154, 0.9981, 0.9977, and 0.9977, associated with excess all-cause mortality. Regardless of the positive or negative directions, the data suggested that weekly NO2 and PM2.5 levels were associated with lower relative risks for excess mortality than COVID-19 mortality—and that weekly park and workplace visits were associated with higher relative risks for excess mortality than COVID-19 mortality.

4. Discussion

This study estimated the excess all-cause mortality during COVID-19 and found differences in the impacts of air pollution and human activity. Because of lockdowns and social distancing policies put in place to control COVID-19 transmission, anthropogenic emissions and regular human activities in many countries were significantly influenced. In particular, levels of air pollutants—including NO2 and PM2.5, which are considered major sources of anthropogenic emissions in public health issues—were altered. Meanwhile, researchers have discussed the potential impacts of working locations and parks due to movement restrictions [40,41]. This study adapted data on exposure to NO2 and PM2.5 and the frequencies of visits to parks and workplaces to depict variances in air pollution and human activity patterns. The findings revealed that the analysis based on excess all-cause mortality reduced the potential impact of NO2 and PM2.5 emissions and enhanced the influence of visits to parks and workplaces, as compared with the estimation based on COVID-19 mortality.
Although movement restrictions significantly decreased NO2 emissions, a positive association was still observed in NO2 related to excess mortality. Opposite trends were revealed for the association between PM2.5 and excess mortality in most countries. It should be noted that the findings of this study were not consistent with many other studies, which suggested an inverse association between exposure to NO2 and mortality and a positive association between PM2.5 and mortality [42,43]. The negative association between PM2.5 and excess mortality could be attributed to lockdown protocols and restricted social distancing policies, which significantly changed air pollutant emissions. In addition, the association discrepancies may be due to the heterogeneity of environmental and socioeconomic conditions among the studied countries. Controlling factors—such as meteorological and demographic patterns—could also change the accuracy of impact estimation. The use of different statistical methods in different studies is another factor that could cause estimation bias [44]. The inverse association between visits to parks and workplaces and excess all-cause mortality is in line with governmental response policies.
Limitations exist in this study. First, the inconsistency of COVID-19 mortality data and all-cause mortality data among countries should be considered. The statistics of COVID-19 cases and all-cause mortality cases varied widely among the studied countries, which could lead to biased estimations of excess all-cause mortality in specific countries. Incomplete COVID-19 cases in many countries could also cause this issue. Moreover, the number of deaths classified as COVID-19 deaths which may have been attributable to other high-risk conditions is still uncertain. While many countries, faced with this issue, have tended to attribute deaths to COVID-19, there are also countries that consider other underlying causes of death in addition to COVID-19 [10]. Second, due to the limited availability of weekly all-cause mortality data, this study focused on country-level analyses and ignored the heterogeneity of air pollution and human activities within countries. Thus, this study mainly investigated variations in air pollution, human activity, and excess mortality responses among countries, and compared the differences between excess all-cause mortality and COVID-19 mortality. Results suggested that it would be necessary to calculate excess all-cause mortality on a finer scale in order to support the estimation of potential air pollution and human activity impacts on regional, national and global scales.

5. Conclusions

This study estimated excess all-cause mortality during COVID-19, based on a SARIMA model, and further investigated the impacts of air pollution and human activity on mortality variations in sixteen countries using a quasi-likelihood Poisson-based GAM model. The results showed that COVID-19 dramatically increased death counts compared with benchmark mortality trends, leading to high rates of excess all-cause mortality. Moreover, the excess all-cause mortality changed from −26.52% to 373.60% when compared with COVID-19 mortality during the 10th–52nd weeks of 2020 in sixteen countries. For the impacts of air pollution and human activity, an increase of 1.54% and decreases of 0.19%, 0.23%, and 0.23% in excess all-cause mortality rates were revealed with a 1-unit increase of weekly NO2, PM2.5, park visits, and workplace visits, respectively. In addition, compared with the estimated relative risks of COVID-19 mortality, lower relative risks for weekly NO2 and PM2.5, and higher relative risks for weekly park and workplace visits were shown for excess all-cause mortality. This suggests that, compared with excess all-cause mortality, COVID-19 mortality overestimated the impact of air pollution and underestimated the influence of human activity. These findings provide an alternative perspective for policymakers and stakeholders to better assess the COVID-19 pandemic.

Author Contributions

Conceptualization, Y.M. and M.S.W.; data curation, methodology and validation, Y.M.; formal analysis, Y.M. and R.Z.; writing—original draft preparation, Y.M.; writing—review and editing, Y.M., M.S.W., H.X., M.-P.K. and R.Z; project administration and funding acquisition, M.S.W.; All authors have read and agreed to the published version of the manuscript.

Funding

Man Sing Wong thanks the funding support from a grant by the General Research Fund (Grant no. 15603920), the Collaborative Research Fund (Grant no. C7064-18GF), and the Research Institute for Sustainable Urban Development (Grant no. 1-BBWD), the Hong Kong Polytechnic University. Hanfa Xing expresses thanks for the funding support from a grant by the National Natural Science Foundation of China (Grant no. 41971406). Mei-Po Kwan was supported by grants from the Hong Kong Research Grants Council (General Research Fund Grant no. 14605920; Collaborative Research Fund Grant no. C4023-20GF) and a grant from the Research Committee on Research Sustainability of Major Research Grants Council Funding Schemes of the Chinese University of Hong Kong.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schlosser, F.; Maier, B.F.; Jack, O.; Hinrichs, D.; Zachariae, A.; Brockmann, D. COVID-19 lockdown induces disease-mitigating structural changes in mobility networks. Proc. Natl. Acad. Sci USA 2020, 117, 32883–32890. [Google Scholar] [CrossRef] [PubMed]
  2. Fisher, D.; Wilder-Smith, A. The global community needs to swiftly ramp up the response to contain COVID-19. Lancet 2020, 395, 1109–1110. [Google Scholar] [CrossRef]
  3. Brooks, S.K.; Webster, R.K.; Smith, L.E.; Woodland, L.; Wessely, S.; Greenberg, N.; Rubin, G.J. The psychological impact of quarantine and how to reduce it: Rapid review of the evidence. Lancet 2020, 395, 912–920. [Google Scholar] [CrossRef] [Green Version]
  4. Liu, J.J.; Bao, Y.; Huang, X.; Shi, J.; Lu, L. Mental health considerations for children quarantined because of COVID-19. Lancet Child. Adolesc. Health 2020, 4, 347–349. [Google Scholar] [CrossRef] [Green Version]
  5. Xiang, Y.-T.; Yang, Y.; Li, W.; Zhang, L.; Zhang, Q.; Cheung, T.; Ng, C.H. Timely mental health care for the 2019 novel coronavirus outbreak is urgently needed. Lancet Psychiatry 2020, 7, 228–229. [Google Scholar] [CrossRef] [Green Version]
  6. Pietrobelli, A.; Pecoraro, L.; Ferruzzi, A.; Heo, M.; Faith, M.; Zoller, T.; Antoniazzi, F.; Piacentini, G.; Fearnbach, S.N.; Heymsfield, S.B. Effects of COVID-19 lockdown on lifestyle behaviors in children with obesity living in Verona, Italy: A longitudinal study. Obesity 2020, 28, 1382–1385. [Google Scholar] [CrossRef] [PubMed]
  7. Lai, A.G.; Pasea, L.; Banerjee, A.; Denaxas, S.; Katsoulis, M.; Chang, W.H.; Williams, B.; Pillay, D.; Noursadeghi, M.; Linch, D. Estimating excess mortality in people with cancer and multimorbidity in the COVID-19 emergency. MedRxiv 2020. [Google Scholar] [CrossRef]
  8. Bilinski, A.; Emanuel, E.J. COVID-19 and excess all-cause mortality in the US and 18 comparison countries. JAMA 2020, 324, 2100–2102. [Google Scholar] [CrossRef]
  9. Arolas, H.P.; Acosta, E.; López-Casasnovas, G.; Lo, A.; Nicodemo, C.; Riffe, T.; Myrskylä, M. Years of life lost to COVID-19 in 81 countries. Sci. Rep. 2021, 11, 3504. [Google Scholar] [CrossRef]
  10. Leon, D.A.; Shkolnikov, V.M.; Smeeth, L.; Magnus, P.; Pechholdová, M.; Jarvis, C.I. COVID-19: A need for real-time monitoring of weekly excess deaths. Lancet 2020, 395, e81. [Google Scholar] [CrossRef]
  11. Modi, C.; Boehm, V.; Ferraro, S.; Stein, G.; Seljak, U. Total covid-19 mortality in Italy: Excess mortality and age dependence through time-series analysis. medRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  12. Rivera, R.; Rosenbaum, J.E.; Quispe, W. Excess mortality in the United States during the first three months of the COVID-19 pandemic. Epidemiol. Infect. 2020, 148, e264. [Google Scholar] [CrossRef]
  13. Stang, A.; Standl, F.; Kowall, B.; Brune, B.; Böttcher, J.; Brinkmann, M.; Dittmer, U.; Jöckel, K.-H. Excess mortality due to COVID-19 in Germany. J. Infect. 2020, 81, 797–801. [Google Scholar] [CrossRef]
  14. Nogueira, P.J.; de Araújo Nobre, M.; Nicola, P.J.; Furtado, C.; Carneiro, A.V. Excess mortality estimation during the COVID-19 pandemic: Preliminary data from Portugal. Acta Med. Port. 2020, 33, 376–383. [Google Scholar] [CrossRef]
  15. Gibertoni, D.; Adja, K.Y.C.; Golinelli, D.; Reno, C.; Regazzi, L.; Lenzi, J.; Sanmarchi, F.; Fantini, M.P. Patterns of COVID-19 related excess mortality in the municipalities of Northern Italy during the first wave of the pandemic. Health Place 2021, 67, 102508. [Google Scholar] [CrossRef]
  16. Krieger, N.; Chen, J.T.; Waterman, P.D. Excess mortality in men and women in Massachusetts during the COVID-19 pandemic. Lancet 2020, 395, 1829. [Google Scholar] [CrossRef]
  17. Lancet, T. The gendered dimensions of COVID-19. Lancet 2020, 395, 1168. [Google Scholar] [CrossRef]
  18. Sinnathamby, M.A.; Whitaker, H.; Coughlan, L.; Bernal, J.L.; Ramsay, M.; Andrews, N. All-cause excess mortality observed by age group and regions in the first wave of the COVID-19 pandemic in England. Eurosurveillance 2020, 25, 2001239. [Google Scholar] [CrossRef] [PubMed]
  19. Blangiardo, M.; Cameletti, M.; Pirani, M.; Corsetti, G.; Battaglini, M.; Baio, G. Estimating weekly excess mortality at sub-national level in Italy during the COVID-19 pandemic. PLoS ONE 2020, 15, e0240286. [Google Scholar] [CrossRef] [PubMed]
  20. Banerjee, A.; Pasea, L.; Harris, S.; Gonzalez-Izquierdo, A.; Torralbo, A.; Shallcross, L.; Noursadeghi, M.; Pillay, D.; Sebire, N.; Holmes, C. Estimating excess 1-year mortality associated with the COVID-19 pandemic according to underlying conditions and age: A population-based cohort study. Lancet 2020, 395, 1715–1725. [Google Scholar] [CrossRef]
  21. Dong, Z.; Fan, X.; Wang, J.; Mao, Y.; Luo, Y.; Tang, S. Data-related and methodological obstacles to determining associations between temperature and COVID-19 transmission. Environ. Res. Lett. 2021, 16, 034016. [Google Scholar] [CrossRef]
  22. Chang, S.; Pierson, E.; Koh, P.W.; Gerardin, J.; Redbird, B.; Grusky, D.; Leskovec, J. Mobility network models of COVID-19 explain inequities and inform reopening. Nature 2021, 589, 82–87. [Google Scholar] [CrossRef] [PubMed]
  23. Pullano, G.; Valdano, E.; Scarpa, N.; Rubrichi, S.; Colizza, V. Evaluating the effect of demographic factors, socioeconomic factors, and risk aversion on mobility during the COVID-19 epidemic in France under lockdown: A population-based study. Lancet Digit. Health 2020, 2, e638–e649. [Google Scholar] [CrossRef]
  24. Grantz, K.H.; Meredith, H.R.; Cummings, D.A.; Metcalf, C.J.E.; Grenfell, B.T.; Giles, J.R.; Mehta, S.; Solomon, S.; Labrique, A.; Kishore, N. The use of mobile phone data to inform analysis of COVID-19 pandemic epidemiology. Nat. Commun. 2020, 11, 4961. [Google Scholar] [CrossRef] [PubMed]
  25. Burnett, R.; Chen, H.; Szyszkowicz, M.; Fann, N.; Hubbell, B.; Pope, C.A.; Apte, J.S.; Brauer, M.; Cohen, A.; Weichenthal, S. Global estimates of mortality associated with long-term exposure to outdoor fine particulate matter. Proc. Natl. Acad. Sci USA 2018, 115, 9592–9597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Achakulwisut, P.; Brauer, M.; Hystad, P.; Anenberg, S.C. Global, national, and urban burdens of paediatric asthma incidence attributable to ambient NO2 pollution: Estimates from global datasets. Lancet Planet. Health 2019, 3, e166–e178. [Google Scholar] [CrossRef] [Green Version]
  27. Venter, Z.S.; Aunan, K.; Chowdhury, S.; Lelieveld, J. Air pollution declines during COVID-19 lockdowns mitigate the global health burden. Environ. Res. 2021, 192, 110403. [Google Scholar] [CrossRef]
  28. Barnett-Itzhaki, Z.; Levi, A. Effects of chronic exposure to ambient air pollutants on COVID-19 morbidity and mortality-A lesson from OECD countries. Environ. Res. 2021, 195, 110723. [Google Scholar] [CrossRef]
  29. Domingo, J.L.; Rovira, J. Effects of air pollutants on the transmission and severity of respiratory viral infections. Environ. Res. 2020, 109650. [Google Scholar] [CrossRef]
  30. Wu, X.; Nethery, R.C.; Sabath, B.M.; Braun, D.; Dominici, F. Exposure to air pollution and COVID-19 mortality in the United States. MedRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, L. Emerging study on the transmission of the Novel Coronavirus (COVID-19) from urban perspective: Evidence from China. Cities 2020, 103, 102759. [Google Scholar] [CrossRef]
  32. Ugolini, F.; Massetti, L.; Calaza-Martínez, P.; Cariñanos, P.; Dobbs, C.; Ostoić, S.K.; Marin, A.M.; Pearlmutter, D.; Saaroni, H.; Šaulienė, I. Effects of the COVID-19 pandemic on the use and perceptions of urban green space: An international exploratory study. Urban. For. Urban. Green. 2020, 56, 126888. [Google Scholar] [CrossRef]
  33. Venter, Z.S.; Barton, D.N.; Gundersen, V.; Figari, H.; Nowell, M. Urban nature in a time of crisis: Recreational use of green space increases during the COVID-19 outbreak in Oslo, Norway. Environ. Res. Lett. 2020, 15, 104075. [Google Scholar] [CrossRef]
  34. Brynjolfsson, E.; Horton, J.J.; Ozimek, A.; Rock, D.; Sharma, G.; TuYe, H.-Y. COVID-19 and Remote Work: An Early Look at US Data; 0898-2937; National Bureau of Economic Research: Cambridge, MA, USA, 2020. [Google Scholar]
  35. Bishop, G.; Welch, G. An introduction to the kalman filter. Proc. SiggraphCourse 2001, 8, 41. [Google Scholar]
  36. Rugge, M.; Zorzi, M.; Guzzinati, S. SARS-CoV-2 infection in the Italian Veneto region: Adverse outcomes in patients with cancer. Nat. Cancer 2020, 1, 784–788. [Google Scholar] [CrossRef]
  37. Albiges, L.; Foulon, S.; Bayle, A.; Gachot, B.; Pommeret, F.; Willekens, C.; Stoclin, A.; Merad, M.; Griscelli, F.; Lacroix, L. Determinants of the outcomes of patients with cancer infected with SARS-CoV-2: Results from the Gustave Roussy cohort. Nat. Cancer 2020, 1, 965–975. [Google Scholar] [CrossRef]
  38. Subbiah, V. A global effort to understand the riddles of COVID-19 and cancer. Nat. Cancer 2020, 1, 943–945. [Google Scholar] [CrossRef]
  39. Iftimie, S.; López-Azcona, A.F.; Vicente-Miralles, M.; Descarrega-Reina, R.; Hernández-Aguilera, A.; Riu, F.; Simó, J.M.; Garrido, P.; Joven, J.; Camps, J. Risk factors associated with mortality in hospitalized patients with SARS-CoV-2 infection. A prospective, longitudinal, unicenter study in Reus, Spain. PLoS ONE 2020, 15, e0234452. [Google Scholar] [CrossRef] [PubMed]
  40. Geng, D.C.; Innes, J.; Wu, W.; Wang, G. Impacts of COVID-19 pandemic on urban park visitation: A global analysis. J. For. Res. 2021, 32, 553–567. [Google Scholar] [CrossRef]
  41. Kwok, C.Y.T.; Wong, M.S.; Chan, K.L.; Kwan, M.-P.; Nichol, J.E.; Liu, C.H.; Wong, J.Y.H.; Wai, A.K.C.; Chan, L.W.C.; Xu, Y. Spatial analysis of the impact of urban geometry and socio-demographic characteristics on COVID-19, a study in Hong Kong. Sci. Total Environ. 2021, 764, 144455. [Google Scholar] [CrossRef]
  42. De Angelis, E.; Renzetti, S.; Volta, M.; Donato, F.; Calza, S.; Placidi, D.; Lucchini, R.G.; Rota, M. COVID-19 incidence and mortality in Lombardy, Italy: An ecological study on the role of air pollution, meteorological factors, demographic and socioeconomic variables. Environ. Res. 2021, 195, 110777. [Google Scholar] [CrossRef] [PubMed]
  43. Zoran, M.A.; Savastru, R.S.; Savastru, D.M.; Tautan, M.N. Assessing the relationship between ground levels of ozone (O3) and nitrogen dioxide (NO2) with coronavirus (COVID-19) in Milan, Italy. Sci. Total Environ. 2020, 740, 140005. [Google Scholar] [CrossRef] [PubMed]
  44. Villeneuve, P.J.; Goldberg, M.S. Methodological considerations for epidemiological studies of air pollution and the SARS and COVID-19 coronavirus outbreaks. Environ. Health Perspect. 2020, 128, 095001. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The performance of the three-year average mortality and SARIMA-based mortality estimations. (a) Three-year average mortality; (b) SARIMA-based mortality.
Figure 1. The performance of the three-year average mortality and SARIMA-based mortality estimations. (a) Three-year average mortality; (b) SARIMA-based mortality.
Ijerph 18 06883 g001
Figure 2. The all-cause observed and benchmark mortality during the 10th–52nd weeks in 2020 in 16 countries (The benchmark mortality is estimated by SARIMA model).
Figure 2. The all-cause observed and benchmark mortality during the 10th–52nd weeks in 2020 in 16 countries (The benchmark mortality is estimated by SARIMA model).
Ijerph 18 06883 g002
Figure 3. The excess and COVID-19 mortalities during the 10th–52nd weeks in 2020 in 16 countries. Excess mortality is calculated based on the benchmark estimated by the SARIMA model.
Figure 3. The excess and COVID-19 mortalities during the 10th–52nd weeks in 2020 in 16 countries. Excess mortality is calculated based on the benchmark estimated by the SARIMA model.
Ijerph 18 06883 g003
Figure 4. Comparison of R square values of GAM, based on excess mortality and COVID-19 mortality with lag effects of 0,1, and 2 weeks.
Figure 4. Comparison of R square values of GAM, based on excess mortality and COVID-19 mortality with lag effects of 0,1, and 2 weeks.
Ijerph 18 06883 g004
Figure 5. The relative risks of NO2 related to excess mortality and COVID-19 mortality, with lag effects of 0,1, and 2 weeks.
Figure 5. The relative risks of NO2 related to excess mortality and COVID-19 mortality, with lag effects of 0,1, and 2 weeks.
Ijerph 18 06883 g005
Figure 6. The relative risks of PM2.5 related to excess mortality and COVID-19 mortality with lag effects of 0,1 and 2 weeks.
Figure 6. The relative risks of PM2.5 related to excess mortality and COVID-19 mortality with lag effects of 0,1 and 2 weeks.
Ijerph 18 06883 g006
Figure 7. The relative risks of park visit frequency, related to excess mortality and COVID-19 mortality with lag effects of 0,1, and 2 weeks.
Figure 7. The relative risks of park visit frequency, related to excess mortality and COVID-19 mortality with lag effects of 0,1, and 2 weeks.
Ijerph 18 06883 g007
Figure 8. The relative risks of workplace visitation frequency, related to excess mortality and COVID-19 mortality, with lag effects of 0,1, and 2 weeks.
Figure 8. The relative risks of workplace visitation frequency, related to excess mortality and COVID-19 mortality, with lag effects of 0,1, and 2 weeks.
Ijerph 18 06883 g008
Table 1. COVID-19 and excess all-cause mortality during the 10th–52nd weeks in 2020.
Table 1. COVID-19 and excess all-cause mortality during the 10th–52nd weeks in 2020.
CountryCOVID-19 MortalityExcess All-Cause MortalityVariation
Belgium19,31719,9563.31%
Chile16,44314,315−12.94%
Croatia3671564053.64%
Czech Republic11,05817,22455.76%
France62,63463,7671.81%
Germany30,36650,38565.93%
Hungary904713,90553.70%
Lithuania12545939373.60%
Luxembourg460338−26.52%
Netherland11,09517,53858.07%
Norway42182996.91%
Poland27,11871,237162.69%
Portugal661913,352101.72%
Spain51,74482,35559.16%
Switzerland7316910624.47%
United States333,278445,24233.59%
Table 2. Summary of ground-level air pollution data, human activity data, meteorological data, and excess all-cause mortality.
Table 2. Summary of ground-level air pollution data, human activity data, meteorological data, and excess all-cause mortality.
Min.Max.MeanSD
Weekly NO2 (μg/m3)0.000227.039117.32735.680
Weekly PM2.5 (μg/m3)0.000204.57366.25632.085
Weekly park visits−582.0002859.000250.076483.875
Weekly workplace visits−526.00086.000−191.499102.288
Daily precipitation (inch)0.0002.6160.5170.492
Daily temperature ( )27.31380.23155.11111.324
Weekly excess death count01.29 × 1083.61 × 1061.25 × 107
Table 3. Average country-level relative risks of weekly NO2, PM2.5, park visits, and workplace visits.
Table 3. Average country-level relative risks of weekly NO2, PM2.5, park visits, and workplace visits.
Weekly NO2Weekly PM2.5Weekly Park VisitsWeekly Workplace Visits
Lag 0COVID-19 mortality1.01660.99790.99870.9968
Excess mortality1.01550.99830.9980.9968
Lag 1COVID-19 mortality1.01620.99850.99820.9979
Excess mortality1.01450.99830.99750.9983
Lag 2COVID-19 mortality1.01850.99770.99760.999
Excess mortality1.01620.99770.99760.9979
MeanCOVID-19 mortality1.01710.9980.99810.9979
Excess mortality1.01540.99810.99770.9977
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Meng, Y.; Wong, M.S.; Xing, H.; Kwan, M.-P.; Zhu, R. Assessing the Country-Level Excess All-Cause Mortality and the Impacts of Air Pollution and Human Activity during the COVID-19 Epidemic. Int. J. Environ. Res. Public Health 2021, 18, 6883. https://doi.org/10.3390/ijerph18136883

AMA Style

Meng Y, Wong MS, Xing H, Kwan M-P, Zhu R. Assessing the Country-Level Excess All-Cause Mortality and the Impacts of Air Pollution and Human Activity during the COVID-19 Epidemic. International Journal of Environmental Research and Public Health. 2021; 18(13):6883. https://doi.org/10.3390/ijerph18136883

Chicago/Turabian Style

Meng, Yuan, Man Sing Wong, Hanfa Xing, Mei-Po Kwan, and Rui Zhu. 2021. "Assessing the Country-Level Excess All-Cause Mortality and the Impacts of Air Pollution and Human Activity during the COVID-19 Epidemic" International Journal of Environmental Research and Public Health 18, no. 13: 6883. https://doi.org/10.3390/ijerph18136883

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop