Monitoring the Covid-19 epidemics in Italy from mortality data

The monitoring of the coronavirus pandemic (Covid-19) is usually based on indicators, such as the number of new positive cases per day, which heavily depend on the number of tests administered to the population. Since the number of tests increases as the disease spreads, the resulting epidemic curves are unable to fully represent the actual evolution of the Covid-19 pandemic in a given region [1]. In particular, in Italy the number of tests made has increased by more than a factor 4 during March 2020, as a combined result of changing the guidelines and the capacity to perform them. All subsequent epidemiological indicators, number of infected, recovered and dead, have been all affected by biases which are difficult to be estimated and corrected in absence of further information about the tests.

Overall mortality data can be used instead as an alternative source to monitor the status of Covid-19. The potential of mortality data is well known and discussed in several papers [2,3,4,5]. The dataset under study is provided by Istituto Nazionale di Statistica (ISTAT) [6] and Sistema di Sorveglianza Mortalità Giornaliera (SiSMG) [7]. It provides the number of deaths per day divided in age categories and gender, up to 4th April 2020, for ISTAT data, and up to 25th April for SiSMG data, i.e. after about two months from the first official direct transmission Covid-19 case reported in Italy. The dataset also contains the historical data, for the same months, for the years 2015-2019. A selection is applied to the municipalities included in the ISTAT dataset: only towns presenting an excess of 20% with respect to the historical average are present, inducing a potential bias which decreases with the dataset coverage. The combined dataset (ISTAT+SiSMG) refers to about 50% of the infected population, taking into account the regional spread of the epidemics. Big cities, like Milan, Turin and Rome are present. Lombardy, which accounts for about one half of the Covid19 deaths in Italy, has the largest coverage, around 72%. Mortality data are also being made available in other countries [8]: these data are not publicly available in a form that allows a complete statistical analysis.

The time-series data for the mortality excess, for different cities and regions, have been obtained by subtracting an average historical model, taking into account the seasonal mortality trend. This is obtained from the average of 2015-2019, normalised to the number of deaths observed in the first 3 weeks of February. Such a model is designed to best extrapolate the observed mortality rate at the beginning of February 2020 into March and April. As an example, the time-series before (left) and after (right) the historical model subtraction is shown in Fig. 1 for the city of Brescia. The same procedure can be applied to any data category, in particular to other cities and regions. The mortality excess is compared with the official Covid-19 deaths provided daily by the Italian Protezione Civile [9], as shown in the same figure (right). Confirmed deaths for Covid-19 are only available at regional granularity: values for cities have been obtained from the corresponding region, taking into account the fraction of active Covid-19 cases in the province as a function of time, and then rescaling for the population. This calculation assumes a uniform distribution of deaths within the same province; an uncertainty of 20% is assigned, obtained from a comparison between this estimate and the real value for a few cities where mortality data was available.

• Download figure
• Open in new tab

Fig. 1:

Left. distribution of deaths in Brescia as a function of time. Right. baseline-subtracted distribution of deaths in Brescia compared to the official Covid-19 distribution. The ratio of the two distributions is displayed below.

• Download figure
• Open in new tab

Fig. 2:

Baseline-subtracted distribution of the age at death compared to the official Covid-19 distribution.

The mortality excess time-series seems to follow the laboratory-confirmed cases, with the start of the raise at the end of February, when the first Italian case of Covid-19 has been reported. The excess of mortality is higher than the confirmed Covid-19 cases at early days of the pandemic, indicating that at the beginning of the outbreak a large fraction of cases hasn’t been identified by the official monitoring. Later, the two values tend to be closer, suggesting that the effective sampling of tested people has increased with time. Similar behaviour is observed for other cities and regions. A large spread of the ratio between the overall mortality excess and confirmed Covid-19 deaths is present: it is above two for many cities and, in some cases, turns out to be very large, as for the city of Genoa.

Tab. 1 summarizes these values. We also report the aggregated result for the full Italian sample.

The effective coverage of the ISTAT+SiSMG dataset is taken into account to estimate the total number of official deaths. The uncertainty assigned reflects the partial coverage of the data.

Age and gender of the mortality excess differ significantly from the laboratory-confirmed Covid-19 deaths. The following Fig.2 shows the comparison between the distributions subtracted of historical data for cities in the Northern part of Italy (Brescia, Milan, Genoa, and Bologna), where the impact of the virus has been stronger with a mortality of about 0.2%, and for cities of Center-South (Rome, Perugia, Bari, Potenza, Palermo, Messina), where the mortality is about a factor 10 less, with the official Covid-19 deaths. Both distributions look quite different from the official data and also between themselves. This is an indication that an important fraction of the excess in deaths is due to older people, dying at home or at the hospice, without even accessing the healthcare system. The difference between North and Center-South may also suggest that in cities with a smaller Covid-19 mortality the impact of the deaths not directly due to the virus but connected to lockdown could be much more significant.

These data indicate that mortality information can be used to provide a less biased time profile of the pandemic. In general, there is a large excess of deaths, more pronounced at the beginning of the pandemic. The additional information regarding age, gender, and regional differences with respect to the official data may also provide a powerful handle to disentangle the actual Covid-19-related contribution from the deaths due to the stress of the healthcare system or to the reduced access to emergency rooms. The impact of this contribution should be minor, given that the difference between the excess of mortality and the official Covid-19 deaths becomes smaller, as time passes. If made available, the information about the cause of the death and the type of hospitalization, which are usually stored in a mortality database, will disentangle the different components, helping to identify the pure Covid-19 contribution. Similar comparisons can be also done between different countries, shedding light on the mechanisms of the spread of the infection and on its real impact.