We assessed the relation between COVID-19 waves in Italy, which was severely affected during the pandemic. We evaluated the hypothesis that a larger impact from the first wave (February-May 2020) predicts a smaller peak during the second wave (September-October 2020), in the absence of local changes in public health interventions and area-specific differences in time trends of environmental parameters. Based on publicly available data on province-specific SARS-CoV-2 infections and both crude and multivariable cubic spline regression models, we found that for provinces with the lowest incidence rates in the first wave, the incidence in the second wave increased roughly in proportion with the incidence in the first wave until an incidence of about 500-600 cases/100,000 in the first wave. Above that value, provinces with higher incidences in the first wave experienced lower incidences in the second wave. It appears that a comparatively high cumulative incidence of infection, even if far below theoretical thresholds required for herd immunity, may provide noticeable protection during the second wave. We speculate that, if real, the mechanism for this pattern could be depletion of most susceptible individuals and of superspreaders in the first wave. A population learning effect regarding cautious behavior could have also contributed. Since no area-specific variation of the national policy against the SARS-CoV-2 outbreak was allowed until early November 2020, neither individual behaviors nor established or purported environmental risk factors of COVID-19, such as air pollution and meteorological factors, are likely to have confounded the inverse trends we observed in infection incidence over time.
