Abstract
Background Italy was the first country outside China to experience the impact of the COVID-19 pandemic, which resulted in a significant health burden. This study presents an analysis of the excess mortality across the 107 Italian provinces, stratified by sex, age group, and period of the outbreak.
Methods The analysis was performed using a two-stage interrupted time series design using daily mortality data for the period January 2015 – May 2020. In the first stage, we performed province-level quasi-Poisson regression models, with smooth functions to define a baseline risk while accounting for trends and weather conditions and to flexibly estimate the variation in excess risk during the outbreak. Estimates were pooled in the second stage using a mixed-effects multivariate meta-analysis.
Results In the period 15 February – 15 May 2020, we estimated an excess of 47,490 (95% empirical confidence intervals: 43,984 to 50,362) deaths in Italy, corresponding to an increase of 29.5% (95%eCI: 26.8 to 31.9%) from the expected mortality. The analysis indicates a strong geographical pattern, with the majority of excess deaths occurring in northern regions, where few provinces experienced up to 800% increase during the peak in late March. There were differences by sex, age, and area both in the overall impact and in its temporal distribution.
Conclusions This study offers a detailed picture of excess mortality during the first months of the COVID-19 pandemic in Italy. The strong geographical and temporal patterns can be related to implementation of lockdown policies and multiple direct and indirect pathways in mortality risk.
Key Messages
This study evaluated mortality trends in Italy during the COVID-19 pandemic, reporting an excess of 47,490 (95% empirical confidence intervals: 43,984 to 50,362) deaths in the period 15 February – 15 May 2020, corresponding to an increase of 29.5% (95%eCI: 26.8 to 31.9%) from the expected mortality.
There is a strong geographical pattern, with 71.0% of the estimated excess deaths occurring in just three northern regions (Lombardy, Veneto, and Emilia-Romagna), and few provinces showing increases in mortality up to 800% during the peak of the pandemic.
The impact was slightly higher is men compared to women, with 24,655 and 23,125 excess deaths respectively, and varied by age, with higher mortality in the group 70-79 years old and evidence of a lower but measurable risk even in people less than 60.
The analysis by week suggests differential trends, with more delayed impacts in women and elderly, and the risk limited to the early period in Central and Southern Italy, likely related to the implementation of lockdown policies and contributions from direct and indirect risk pathways.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This work was supported by the Medical Research Council-UK (Grant ID: MR/M022625/1) and the European Union's Horizon 2020 Project Exhaustion (Grant ID: 820655).
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
NA
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
The data, software and code for replicating the analysis and full set of results are made available. The analysis was performed in the R software environment. The scripts for downloading the data and for fully replicating the analysis and results are available in a GitHub repository (https://github.com/gasparrini/ItalyCOVIDdeath). The data (release on 18/06/2020 with updated mortality up to 15th of May) are extracted from the ISTAT website (https://www.istat.it/it/archivio/240401) and stored in the repository. A web application is made available at (https://mscortichini.shinyapps.io/app20200703/) to visualize maps of the excess risk and extract the full set of results, which are also available in the Online Appendix of this article.
https://www.istat.it/it/archivio/240401