Skip to main content

The association between first and second wave COVID-19 mortality in Italy



The relation between the magnitude of successive waves of the COVID-19 outbreak within the same communities could be useful in predicting the scope of new outbreaks.


We investigated the extent to which COVID-19 mortality in Italy during the second wave was related to first wave mortality within the same provinces. We compared data on province-specific COVID-19 2020 mortality in two time periods, corresponding to the first wave (February 24–June 30, 2020) and to the second wave (September 1–December 31, 2020), using cubic spline regression.


For provinces with the lowest crude mortality rate in the first wave (February–June), i.e. < 22 cases/100,000/month, mortality in the second wave (September–December) was positively associated with mortality during the first wave. In provinces with mortality greater than 22/100,000/month during the first wave, higher mortality in the first wave was associated with a lower second wave mortality. Results were similar when the analysis was censored at October 2020, before the implementation of region-specific measures against the outbreak. Neither vaccination nor variant spread had any role during the study period.


These findings indicate that provinces with the most severe initial COVID-19 outbreaks, as assessed through mortality data, faced milder second waves.

Peer Review reports


The COVID-19 pandemic in many countries has been characterized by waves of infection. These recurrent local exacerbations of the pandemic present an opportunity to study the spread of the virus within a population. Key features of the pandemic are still not well understood, such as the susceptibility of the population to subsequent waves after the first outbreak, the threshold for herd immunity, the role of superspreaders [1,2,3,4,5,6,7] as well as of meteorological and environmental factors [8,9,10,11].

We compared the two COVID-19 waves within Italy, where the geographical distribution of the SARS-CoV-2 infection spread and the COVID-19 incidence was uneven during the first wave of the pandemic [12]. From an analysis within provinces of the official case counts during the first and second wave of SARS-CoV-2 infection, we found evidence of an association in the occurrence between the two periods. When incidence in the first wave was low (< 500 cases/100,000/day), the second wave incidence tended to be positively correlated, whereas a high first-wave incidence was strongly and inversely correlated with the second wave incidence. This observation raised the possibility that greater spread at the beginning of the pandemic could have induced some protection at the population level, resulting in a milder second wave, despite low levels (10–20%) of anti-SARS-CoV-2 antibody seroprevalence even in the areas most severely hit.

One problem with our earlier finding is that the incidence of SARS-CoV-2 infection is difficult to measure. It depends on the severity of the clinical course, the availability of molecular tests and the policy underlying their use at the population level. In Italy, for instance, the implementation of nasopharyngeal swabs and molecular tests to detect SARS-CoV-2 infection in the first wave was restricted primarily to symptomatic patients with suspected COVID-19. In the summer of 2020, however, the public-health policy was modified to extend testing to contacts of cases and to individuals without clinical symptoms, nearly tripling the number of daily molecular tests from the first to the second wave time periods [13]. An additional change of policy, starting on November 6, 2020, during the second wave, was the mandatory adoption of differential area-specific public-health measures [14, 15].

The area-specific relation between the two waves may be better assessed by a more objective measures of COVID-19 spread such as mortality rate, instead of incidence rates that are subject to swab testing policy and availability. Recently, COVID-19 death data for Italy in 2020 have become available with greater geographic detail, with figures available for the 107 provinces rather than the 20 larger regions within Italy. We used these mortality data to revisit the relation of the dynamics of the first and second waves of the outbreak.


Study population and outcome

We used data made freely available from public sources. We accessed 2020 monthly province-specific COVID-19 deaths as reported by the Italian National Institute of Statistics – ISTAT [16]. Provinces are the administrative entities that are intermediate between smaller municipalities and larger regions within Italy. Within 2020, we examined two time frames: we considered the first wave to be the period from the pandemic onset (February 24) to June 30, when the number of new infection cases had substantially dropped after the March–April peak. From July to the end of August the number of new cases remained low, mostly related to screening of subjects returning from vacation with swab testing. Cases increased again beginning in early September. Therefore, we considered the second wave to be the period from September 1 to December 31, 2020. During these periods, neither circulation of virus variants nor implementation of SARS-CoV-2 vaccination at the population level had yet started in Italy. Only a handful of Italians had been vaccinated by the end of December, but large scale administration started in January 2021 [17, 18].

Based on these data and the Italian population at January 1, 2020, and January 1, 2021 available at the ISTAT website [19], we computed wave-specific COVID-19 mortality during the first and second waves. For the former period, we used as reference population data the provincial population size as of January 1, 2020. For the second period and for the overall 2020 mortality we averaged the population size at January 1, 2020 and January 1, 2021 for each province. We retrieved the seroprevalence data for each Italian region and province determined through a national survey made available by ISTAT and the National Ministry of Health with reference to the period May 25, 2020 through July 15, 2020 [20].

Data analysis

We examined the relation between first and second wave province-specific crude mortality rates from COVID-19 by modeling the mean mortality rate using a restricted cubic spline regression [21], a model that fits a curvilinear pattern to the data as previously reported [22,23,24]. In particular, we fitted a restricted cubic spline model that weighted provincial data by population size, performing both a crude analysis and adjusting for potential confounders. Variables included in the model were aging index (i.e. a ratio between resident population aged ≥65 years and those aged ≤14 years), percentage commuting outside the municipality of residence on a daily basis, and percentage of dwellings occupied by only one resident (available at a provincial level from the National Institute of Statistics [19]). In the regression model, we used three knots at fixed percentiles (10th, 50th and 90th) of the first wave distribution, and computed a pointwise 95% confidence interval (CI) [21, 25]. We used ‘mkspline’, ‘regress’, and ‘xbrcsplinei’ routines of Stata software (v17.0, College Station-TX, 2021) for all analyses.


Table 1 and Fig. 1 lists province-specific mortality rates in the investigated periods. National cumulative mortality (cumulative incidence of death) during the first wave averaged 58.2/100,000 persons, ranging from 6.5/100,000 persons in the Southern Italy, Basilicata region, to 165/100,000 persons in Northern Italy, Lombardy region. Corresponding figures for the second wave period were 67.6/100,000 persons for the national average, ranging from 19.4/100,000 in the Calabria region to 192/100,000 persons in the Aosta Valley region. In particular, seven provinces, five of the in Lombardy region (Bergamo, Brescia, Cremona, Lodi, and Pavia) and two in Emilia-Romagna region (Parma and Piacenza) experienced the highest mortality rates during the first wave. Second wave mortality was lower than first wave, always ≤100/100,000 persons and generally lower than the regional or national averages. The provinces with the highest first wave mortality, namely Piacenza (with the overall highest mortality in Italy with 333/100,000 persons), Cremona, Lodi and Bergamo, experienced the highest absolute decrease in mortality rates.

Table 1 Number of SARS-CoV-2 cases, COVID-19 deaths and COVID-19 mortality rates (deaths/100,000/wave timeframe) in the 1st and 2nd waves in 2020 divided by province
Fig. 1
figure 1

COVID-19 mortality in the Italian provinces in the first two pandemic waves and total 2020 mortality

In the spline analysis, the relation between mortality in the first and second wave was U-shaped (Fig. 2), with a direct association between the estimates of the two waves up to 88 deaths /100,000, and above that an inverse pattern, with low second wave figures for those provinces with the highest rates during the first wave (88 through 192 deaths/100,000). The latter corresponded to a mortality rate of 22/100,000/month, as only 5 deaths were recorded in February 2020 following the first disease diagnosis in Italian residents on February 20. Crude analysis, ignoring control of potential confounders such as proportion elderly, living alone, and degree of mobility, gave similar results. When we repeated the main analysis by limiting the second wave to September–October 2020, i.e. by removing the November–December period when region-specific public health and social distancing policies were first allowed and implemented, results were not substantially changed, though mortality in this truncated period was much lower (Fig. 2).

Fig. 2
figure 2

Population-weighted cubic spline regression analysis of the relation between first (from February 24 to June 30, 2020) and (A) second wave (from September 1 to December 31, 2020) mortality rates of SARS-CoV-2 infection in Italian provinces or (B) second wave before the beginning of mobility restrictions (from September 1 to October 31, 2020) mortality rates of SARS-CoV-2 infection in Italian provinces adjusted for aging index, percentage commuting outside the municipality of residence on a daily basis, and percentage of dwellings occupied by only one resident. Shadow area is the 95% confidence interval of the predicted mean mortality rate (solid line); circles are province-specific values weighted for population size


Italy was the first country where the SARS-CoV-2 infection swept out of control, following the initial outbreak in China, where it was kept under control and geographically confined. The early arrival of the pandemic in Italy allowed a longer time frame to monitor and study the behavior of the epidemic. Italy did not allow the adoption of area-specific policies to curtail the outbreak until November, using instead a nationwide approach imposed by the Italian Government. It has been suggested recently that severity of COVID-19 pneumonia in hospitalized subjects was lower during the second wave [26], possibly due to SARS-CoV-2 genomic variation [27]. However, during the time frame of this study there was no evidence of spread of SARS-CoV-2 variants [17], thus avoiding any confounding effect from differential geographic spread of viral lineages with different transmission or virulence features. In addition, recent assessment of fatality rates between first and second waves in Italy showed that they were comparable [28]. Nonetheless, we cannot rule out that improvement in COVID-19 treatment may have decreased mortality during the pandemic, especially in the second wave, despite the lack of fully effective treatments [29,30,31]. Finally, vaccination had no role in curbing either of the two waves analyzed here, since almost no vaccine doses were administered in Italy in 2020. Indeed, the vaccination campaign had not begun yet in the investigated period, apart from a handful of ‘demonstration’ vaccinations administered to health professionals on December 27, 2020 at the national hub of the Spallanzani Hospital in Rome. Large scale administration started in January 2021 [32].

In contrast to our earlier analysis that relied on SARS-CoV-2 infection incidence, when area-specific mortality was still unavailable, we used mortality from COVID-19 to gauge the severity of the waves.

Mortality is a far more reliable indicator than incidence for documenting the actual trends of the pandemic, because the detection and notification of newly-infected cases from the general population depends on the availability and implementation of molecular testing, which in turn depends on health authority guidelines and policies [33]. Therefore, the number and type of individuals to be subject to SARS-CoV-2 testing may differ considerably across geographic areas or time periods independently of real differences in incidence, due to changing referral guidelines (such as the presence of COVID-19 symptoms) and availability of molecular tests. In Italy, first wave health policy allowed COVID-19 molecular testing only for COVID-19 symptomatic individuals, while in the second wave much more testing occurred, with systematic extension to asymptomatic or weakly symptomatic individuals and therefore a great increase in number of notifications. As a consequence, the ‘official’ SARS-CoV-2 infection incidence strongly increased in the second wave. This increase, however, could have been an artifact due to the aforementioned changes in testing policy and availability, especially in light of the roughly comparable numbers of deaths in the two waves. However, incidence is less reliable compared to mortality in documenting the real trend of the pandemic, also due to its inherent differences compared with SARS-CoV-2 infection incidence and with the hybrid pattern of the related assessments at the population level, due to changes in molecular testing availability and policy.

We were able to confirm the previously reported incidence patterns for SARS-CoV-2 in Italy during 2020 [12]. The key findings were a positive association between the two waves when the epidemic was not particularly strong in the first part of the year, and a negative association when the first wave was severe. This finding was similar when using a different time frame for the second wave to account for the possible influence of local policies. Local policies were permitted beginning as of November 6, 2020. In addition, the pattern we investigated is based on ‘real’ data, compared to studies presenting predictive modeling [34, 35], thus strengthening our findings.

Provinces that were relatively lightly affected by the first wave had a similar experience in the second wave, with rates closely mirroring those experienced during the first phase of the outbreak. In contrast, the Italian provinces severely hit by the first wave, such as Bergamo, Lodi, Cremona and Piacenza, experienced a marked reduction of mortality in the second wave, and presumably in the spread of COVID-19. The reasons for such an inverse pattern, however, were not directly identifiable through this study and may only be hypothesized. One possibility is a higher likelihood of immunity for superspreaders and individuals with high mobility and propensity to transmission during the first wave [36]. Another possible explanation is the development of an additional cellular immunity against SARS-CoV-2 in the most severely hit population during the first wave [37, 38]. The provinces most severely affected by the first wave showed a population-based seroprevalence in the 10–25% range, far below the estimated threshold for herd immunity of 50–70% [39]. It is possible that the prevalence of immunity, despite being below the herd immunity threshold, was high enough in provinces with a severe first wave to interact with behavioral factors and lessen the intensity of the second wave.

Availability of data and materials

We used data made freely available from public sources. All data generated or analysed during this study are included in this published article.


  1. Kault D. Superspreaders, asymptomatics and COVID-19 elimination. Med J Aust. 2020;213(10):447–8 e1.

    Article  PubMed  Google Scholar 

  2. Xia Y, Zhong L, Tan J, Zhang Z, Lyu J, Chen Y, et al. How to understand "herd immunity" in COVID-19 pandemic. Front Cell Dev Biol. 2020;8:547314.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Okell LC, Verity R, Watson OJ, Mishra S, Walker P, Whittaker C, et al. Have deaths from COVID-19 in Europe plateaued due to herd immunity? Lancet. 2020;395(10241):e110–e1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Randolph HE, Barreiro LB. Herd immunity: understanding COVID-19. Immunity. 2020;52(5):737–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Kwok KO, Lai F, Wei WI, Wong SYS, Tang JWT. Herd immunity - estimating the level required to halt the COVID-19 epidemics in affected countries. J Inf Secur. 2020;80(6):e32–e3.

    Article  CAS  Google Scholar 

  6. Omer SB, Yildirim I, Forman HP. Herd immunity and implications for SARS-CoV-2 control. JAMA. 2020;324(20):2095–6.

    Article  CAS  PubMed  Google Scholar 

  7. Papachristodoulou E, Kakoullis L, Parperis K, Panos G. Long-term and herd immunity against SARS-CoV-2: implications from current and past knowledge. Pathog Dis. 2020;78:ftaa025.

    Article  CAS  PubMed  Google Scholar 

  8. Marques M, Rovira J, Nadal M, Domingo JL. Effects of air pollution on the potential transmission and mortality of COVID-19: a preliminary case-study in Tarragona Province (Catalonia, Spain). Environ Res. 2021;192:110315.

    Article  CAS  PubMed  Google Scholar 

  9. Filippini T, Rothman KJ, Cocchio S, Narne E, Mantoan D, Saia M, et al. Associations between mortality from COVID-19 in two Italian regions and outdoor air pollution as assessed through tropospheric nitrogen dioxide. Sci Total Environ. 2021;760:143355.

    Article  CAS  PubMed  Google Scholar 

  10. Cherrie M, Clemens T, Colandrea C, Feng Z, Webb DJ, Weller RB, et al. Ultraviolet a radiation and COVID-19 deaths in the USA with replication studies in England and Italy. Br J Dermatol. 2021;185(2):363–70.

    Article  CAS  PubMed  Google Scholar 

  11. Filippini T, Rothman KJ, Goffi A, Ferrari F, Maffeis G, Orsini N, et al. Satellite-detected tropospheric nitrogen dioxide and spread of SARS-CoV-2 infection in northern Italy. Sci Total Environ. 2020;739:140278.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Vinceti M, Filippini T, Rothman KJ, Di Federico S, Orsini N. SARS-CoV-2 infection incidence during the first and second COVID-19 waves in Italy. Environ Res. 2021;197:111097.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Romagnani P, Gnone G, Guzzi F, Negrini S, Guastalla A, Annunziato F, et al. The COVID-19 infection: lessons from the Italian experience. J Public Health Policy. 2020;41(3):238–44.

    Article  PubMed  Google Scholar 

  14. Presidente del Consiglio dei Ministri. DPCM 6 Novembre 2020. 2020 [Accessed August 30, 2021].

  15. Vinceti M, Filippini T, Rothman KJ, Ferrari F, Goffi A, Maffeis G, et al. Lockdown timing and efficacy in controlling COVID-19 using mobile phone tracking. EClinicalMedicine. 2020;25:100457.

    Article  PubMed  PubMed Central  Google Scholar 

  16. ISTAT. COVID-19 deaths in Italy. Report March 5, 2021. [Accessed August 30, 2021].

  17. ELIXIR-IT. SARS-CoV-2 lineages of potential epidemiological interest in Italy. 2021 [Accessed August 30, 2021].

  18. Ministry of Health. Anti-COVID vaccination campaign. 2021 [Accessed August 30, 2021].

  19. ISTAT. Populations and households. Italian National Institute of Statistics. 2020 [Accessed August 31, 2021].

  20. ISTAT. First results of the seroprevalence investigation on SARS-COV-2. Italian National Institute of Statistics. 2020 [Accessed August 30, 2021].

  21. Harrell FE. Regression modeling strategies: with applications to linear models, logistic regression, and survival analysis. New York: Springer; 2001.

    Book  Google Scholar 

  22. Filippini T, Adani G, Malavolti M, Garuti C, Cilloni S, Vinceti G, et al. Dietary habits and risk of early-onset dementia in an Italian case-control study. Nutrients. 2020;12(12):3682.

    Article  PubMed Central  Google Scholar 

  23. Filippini T, Mandrioli J, Malagoli C, Costanzini S, Cherubini A, Maffeis G, et al. Risk of amyotrophic lateral sclerosis and exposure to particulate matter from vehicular traffic: a case-control study. Int J Environ Res Public Health. 2021;18(3):973.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Vinceti M, Rothman KJ, Crespi CM, Sterni A, Cherubini A, Guerra L, et al. Leukemia risk in children exposed to benzene and PM10 from vehicular traffic: a case-control study in an Italian population. Eur J Epidemiol. 2012;27(10):781–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Orsini N. XBRCSPLINE: Stata module to tabulate differences in predicted responses after restricted cubic spline models: Boston College Department of economics; 2009.

    Google Scholar 

  26. Borghesi A, Golemi S, Carapella N, Zigliani A, Farina D, Maroldi R. Lombardy, northern Italy: COVID-19 second wave less severe and deadly than the first? A preliminary investigation. Infect Dis (Lond). 2021;53(5):370–5.

    Article  Google Scholar 

  27. Toyoshima Y, Nemoto K, Matsumoto S, Nakamura Y, Kiyotani K. SARS-CoV-2 genomic variations associated with mortality rate of COVID-19. J Hum Genet. 2020;65(12):1075–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Filippini T, Zagnoli F, Bosi M, Giannone ME, Marchesi C, Vinceti M. An assessment of case-fatality and infection-fatality rates of first and second COVID-19 waves in Italy: COVID-19 fatality rate in Italy. Acta Biomed. 2021;92:e2021420.

    Article  PubMed  Google Scholar 

  29. COVID-19 RISK and Treatments (CORIST) Collaboration. Use of hydroxychloroquine in hospitalised COVID-19 patients is associated with reduced mortality: findings from the observational multicentre Italian CORIST study. Eur J Intern Med. 2020;82:38–47.

  30. Docea AO, Tsatsakis A, Albulescu D, Cristea O, Zlatian O, Vinceti M, et al. A new threat from an old enemy: reemergence of coronavirus (review). Int J Mol Med. 2020;45:1631–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Di Castelnuovo A, Costanzo S, Antinori A, Berselli N, Blandi L, Bonaccio M, et al. Lopinavir/ritonavir and darunavir/cobicistat in hospitalized COVID-19 patients: findings from the multicenter Italian CORIST study. Front Med (Lausanne). 2021;8:639970.

    Article  Google Scholar 

  32. CPD - Italian Civil Protection Department. Anti COVID-19 Vaccine Report. 2021 [Accessed August 30, 2021].

  33. De Matteis A, Turkmen Ceylan FB, Adrianopoli M, Ertek S. The contribution of testing in the fight against COVID-19: evidence from Italy. Eur J Pub Health. 2021;31(4):913–7.

    Article  Google Scholar 

  34. Qian Y, Xie W, Zhao J, Xue M, Liu S, Wang L, et al. Investigating the effectiveness of re-opening policies before vaccination during a pandemic: SD modelling research based on COVID-19 in Wuhan. BMC Public Health. 2021;21(1):1638.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Ahmed DA, Ansari AR, Imran M, Dingle K, Bonsall MB. Mechanistic modelling of COVID-19 and the impact of lockdowns on a short-time scale. PLoS One. 2021;16(10):e0258084.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Signorelli C, Zucchi A, Tersalvi CA, Ciampichini R, Beato E, Balzarini F, et al. High seroprevalence of SARS-CoV-2 in Bergamo: evidence for herd immunity or reason to be cautious? Int J Public Health. 2020;65(9):1815–7.

    Article  PubMed  Google Scholar 

  37. Braun J, Loyal L, Frentsch M, Wendisch D, Georg P, Kurth F, et al. SARS-CoV-2-reactive T cells in healthy donors and patients with COVID-19. Nature. 2020;587(7833):270–4.

    Article  CAS  PubMed  Google Scholar 

  38. Mateus J, Grifoni A, Tarke A, Sidney J, Ramirez SI, Dan JM, et al. Selective and cross-reactive SARS-CoV-2 T cell epitopes in unexposed humans. Science. 2020;370(6512):89–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Fontanet A, Cauchemez S. COVID-19 herd immunity: where are we? Nat Rev Immunol. 2020;20(10):583–4.

    Article  CAS  PubMed  Google Scholar 

Download references


not applicable.


This study was supported by the grants ‘UNIMORE FAR 2019 and 2020 Interdisciplinare’ of the Fondazione di Modena’ to Dr. Vinceti and Dr. Filippini and by the FISR 2020-COVID19 grant by the Italian Ministry of the University and the Research to Dr. Vinceti.

Author information

Authors and Affiliations



MV and TF designed the original study, and with KJR analyzed and interpreted the data, and wrote the manuscript. NO designed and carried out data analysis with TF and SDF. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Marco Vinceti.

Ethics declarations

Ethics approval and consent to participate

not applicable.

Consent for publication

not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vinceti, M., Filippini, T., Rothman, K.J. et al. The association between first and second wave COVID-19 mortality in Italy. BMC Public Health 21, 2069 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: