key: cord-0864188-o3gcmt0e authors: Buonanno, P.; Galletta, S.; Puca, M. title: Spatial dynamics of SARS-Cov-2 and reduced risk of contagion: evidence from the second Italian epidemic wave date: 2020-11-10 journal: nan DOI: 10.1101/2020.11.08.20227934 sha: 176923be7fd5be556d0d4e0cdd1175eddc182f79 doc_id: 864188 cord_uid: o3gcmt0e We highlight a negative association between the severity of the first-wave of SARS-Cov-2 and the spread of the virus during the second-wave. Analyzing data of a sample of municipalities from the Italian region of Lombardy, we find that a one standard deviation increase in excess of mortality during the first-wave is associated with a reduction of approximately 30% in the number of detected infected individuals in the initial phase of the second-wave. Our findings may reflect a behavioral response in more severely hit areas as well as a cross-protection between successive waves. The COVID-19 pandemic has already produced a broad negative effect around the world. As of the end of October, more than 50 million SARS-CoV-2 infections have been detected worldwide, while the number of reported victims passes the 1 million. The surge of new cases during the fall of 2020 in different parts of the globe -in particular many European countries -seems to confirm what most experts were predicting, the arriving of a second-wave of infections (Xu and Li, 2020) . Recent contributions (Cacciapaglia et al., 2020; Fan, Yang, Lin, Zhao, Yang, and He, Fan et al.; Grech and Cuschieri, 2020; Leung et al., 2020; López and Rodó, 2020) forecast and projected the time evolution of COVID-19 pandemic across regions to provide quantitative evidence in support of policymaker and government decision to tackle the spread of the pandemic. Like other countries across the world, Italy, and in particular Lombardy, is experiencing a second wave. Figure 1 shows the first wave in Lombardy in March-May 2020 (spring wave) followed by a second and higher peak in October-November 2020 (fall wave). The COVID-19 pandemic evolution presents some similarities with the 1918 Pandemic Influenza (Spanish Flu) that occurred in three waves from the spring of 1918 to the winter of 1919. As widely discussed in Barry (2004) , the 1918 Pandemic Influenza 2 . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint has been characterized by a mild first wave in spring 1918, an extremely lethal second wave in fall 1918, followed by a less severe third wave in winter and spring 1919. The exposure to influenza in the spring of 1918 provided some immunization during the following waves. As noted in Barry et al. (2008) for the USA: "the intensity of the first wave may have differed across US cities and countries and may partly explain geographical variation in pandemic mortality rates in the fall". In this paper, we aim to study the relationship between the severity of the first-wave of SARS-Cov-2 occurred in Spring 2020 and the spread of the virus during the second-wave in Fall 2020. We focus on municipalities of the Italian region of Lombardy that has been severely hit during the first-wave and sadly known as the European epicenter. Exploiting both a traditional regression framework and a spatial econometric model to account for possible spatial effects and cross-municipalities interactions, we find that a one standard deviation increase in the excess of mortality during the first-wave is associated with a reduction of approximately 30% in the number of detected infected individuals. This could be due both to a behavioral response of individuals in the more severely hit areas and a higher immunization in the first wave that could contribute to a slower spread of the virus during the second wave of the epidemic. Our results suggest that policymakers and health authorities should collaborate to design containment measures tailored only on small geographic entities and that should consider immune protection acquired during previous waves, despite herd immunity is far to be reached (Britton et al., 2020; Brett and Rohani, 2020) . The rest of the paper is organized as follows. We describe data and the empirical methodology in section 2. Section 3 presents the estimation results, and section 5 concludes. Data. We combine data from several sources. Specifically, we create a measure of excess mortality, which compares the number of deaths observed in Lombardy during the first COVID-19 wave (between January 1 and May 31, 2020), to the average number of deaths observed during the same time-span for all the years between 2015 and 2019. We gather mortality data released on October 22, 2020 from the Italian National Statistical Institute (ISTAT). 1 The dataset includes the total mortality of 1,501 out of the 1,506 municipalities of Lombardy, virtually covering the entire regional population. 1 See https://www.istat.it/it/archivio/240401 for more information 3 . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 10, 2020. ; https://doi.org/10. 1101 We merge this information with data on the number of cases observed during the more recent (and not yet ended) second COVID-19 wave, from September 1 to November 1, 2020. 2 The dataset contains, for each Lombardy's municipality, the number of new cases. We use these variables in section 3 to compare the number of COVID-19 cases detected during the second-wave with the number of the excess of deaths observed in the first-wave. To clean our estimations from possible confounding factors, we also collect information on: (i) the age structure of the population and compute the share of individuals over 60 years (Share of over 60 ) from ISTAT; (ii) an indicator of the air quality (Pollution) based on the prevalence of fine particulate matters PM 2.5 over the period 2017 to 2019 (Pollution) collected from ARPA Lombardia (Regional Environmental Protection Agency); and (iii) a proxy for municipality civic capital, that is the share of recycling of urban solid waste (Civic capital ) collected from ISTAT, which should capture the propensity of individuals to contribute to a public good and to comply with legal and social norms. Table 1 reports basic summary statistics for the variables used in our analysis. Methods. In section 3 we analyze our dataset using both ordinary least-squares (OLS) and generalised spatial two stage least squares (GS2SLS) estimator of Kelejian and Prucha (1998) . Moreover, in our OLS regression we account for spatial correlation across nearby units using Conley standard errors (Conley, 1999) . Specifically, we estimate both 2 We obtained data from the Lombardy Region, which are available at https://datawrapper.dwcdn. net/iMArO/10/ and https://www.datawrapper.de/_/567DW/ 4 . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 10, 2020. ; https://doi.org/10.1101/2020.11.08.20227934 doi: medRxiv preprint the following model where y is the second-wave number of detected cases per 100k inhabitants, x is the excess mortality per 100k inhabitants observed during the first epidemic wave, α is the constant, Z is a set of control variables, and ε is a disturbance term. We are now ready to illustrate the results of our estimations. Basics. We start our investigation analyzing the geographic distributions of the firstwave excess mortality and the second-wave number of cases. Figure 2 reports for each municipality of Lombardy, the first-wave excess of mortality per 100k inhabitants (lefthand-side panel) and the second-wave cumulative number of positive cases per 100k inhabitants (right-hand-side panel). Even a rough inspection of the maps reveals that the number of second-wave cases is lower in the municipalities that were more severely hit during the first epidemic wave, suggesting that the risk of contagion might be lower in those areas. . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 10, 2020. ; https://doi.org/10.1101/2020.11.08.20227934 doi: medRxiv preprint (Pollution) and (iii) a proxy of a municipality civic capital, that is the share of recycling of urban solid waste (Civic capital ). Robust standard errors are presented in parentheses. *, ** and *** denote rejection of the null hypothesis of the coefficient being equal to 0 at 10%, 5% and 1% significance level, respectively. In column 6, standard errors are computed using Conley standard errors (Conley, 1999) . To confirm our intuition, we perform a more in-depth empirical investigation with regression analysis. Table 2 reports ordinary least-squares regressions with the number of COVID-19 cases per 100k inhabitants detected during the second wave (Secondwave cases) as the dependent variable and the excess of mortality per 100k inhabitants observed during the first wave (First-wave excess mortality) as the main explanatory variable. The excess-mortality coefficient has a negative sign, and it is statistically significantly different from zero at 5% level. This relationship is robust to the inclusion of a minimal set of demographic controls such as population density (Population density), the share of individuals older than 60 (Share of over 60 ), a proxy of air pollution (Pollution), and a proxy of civic capital measured by the share of solid-waste recycling (Civic capital ). We report in columns (1)-(5) robust standard errors and, in column (6) suggests that an increase of one standard deviation of the first-wave excess mortality corresponds approximately to a 30% reduction of detected second-wave cases. Finally, we illustrate this negative relationship in Figure 3 , which is the visual representation of the model (5) of Table 2 , i.e., the binned scatter plot of a regression model including the whole set of control variables. 6 . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this this version posted November 10, 2020. ; Spatial dynamics. As an additional check we control for possible spatial effects and crossmunicipalities interactions (e.g. Drukker et al., 2013) . There is no reason to believe that infections and contagions follow the administrative boundaries of municipalities. Omitting to take them into account may reduce the efficiency of our estimates and bias them. To address this issue, we estimate a spatial model using the generalized spatial two-stage least squares (GS2SLS) estimator of Kelejian and Prucha (1998) . We report in Table 3 the results of spatial estimations of the correlation between first-wave excess mortality and second-wave cases, which reproduces the same specification of column 5 of 2. We employ a contiguity matrix, and we implement a spatial autoregressive model (column (1)), a spatial error model (column (2)), and a model that combines the two by considering both a spatial lag and a spatial error structure (column (3)). Allowing for a spatial structure in our data does not alter our baseline estimates: the coefficient of First-wave excess mortality is negative and statistically significant at 1%. The role of civic capital. Resent research has highlighted the role of civic capital in slowing down the spread of COVID-19 (Durante et al., 2020) . To further explore such 7 . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 10, 2020. ; https://doi.org/10.1101/2020.11.08.20227934 doi: medRxiv preprint Notes: This table presents the results of a spatial model estimated using the generalized spatial two-stage least squares (GS2SLS) estimator of Kelejian and Prucha (1998) . Included controls are the same as in the specification of column 5 of Table 2 . We employ a contiguity matrix. A spatial lag model, a spatial error model, and a model that combines the two by considering both a spatial lag and a spatial error structure are respectively presented in columns (1), (2), and (3). λ is the spatial lag term, while ρ is the spatial error. The dependent variable is Second-wave cases, the COVID-19 cases per 100k inhabitants detected during the second wave. The main explanatory variable is First-wave excess mortality the excess of mortality observed during the first wave per 100k inhabitants. Controls include: (i) the share of individuals over 60 years (Share of over 60 ); (ii) an indicator of the air quality (Pollution) based on the prevalence of fine particulate matters PM 2.5 over the period 2017 to 2019 (Pollution) and (iii) a proxy of a municipality civic capital, that is the share of recycling of urban solid waste (Civic capital ). Robust standard errors are presented in parentheses. *, ** and *** denote rejection of the null hypothesis of the coefficient being equal to 0 at 10%, 5% and 1% significance level, respectively. a phenomenon, we decompose the effect of First-wave excess mortality on Second-wave cases by interacting the former with the quartiles of Civic capital. While we omit the estimation results (available upon request) for brevity, we observe that the the reduction in the number of cases is significantly stronger for municipalities with higher levels of civic capital. Specifically, the coefficients of the interactions between Second-wave cases and the quartiles of Civic capital are all negative, and statistically significantly different from zero at 1% level. These results are consistent with the idea that civic-minded individuals internalize more than others the effect of their behavior on the spread of an infectious disease. 8 . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 10, 2020. ; Two potentially simultaneous effects could explain the observed negative association between the first-wave severity and the current second-wave reduced risk of contagion. On the one hand, people living in areas that were severely affected during the first-wave are likely to have changed their behavior in a way that would make less likely the spread of the virus (e.g., more accurate use of protecting devices, more frequent sanitizing procedures, etc.). This might not be true in areas where, instead, many individuals did not have actual direct or indirect exposure to the consequences of the virus (e.g., a dead friend or relative). However, since the inclusion of a proxy of civic capital as a control variable did not affect our estimates, the observed heterogeneity could be explained by alternative mechanisms. On the other hand, severe exposure to the virus during the first-wave may have immunized most of the population, contributing to a much slower spread of the virus during the second wave of the epidemic. As reported in previous work (Buonanno et al., 2020; Perico et al., 2020) certain provinces in Lombardy reported a very high level of predicted infections (considering a hypothetical fatality rate of 1%, the share of the population infected is estimated to be: 40% in the province of Bergamo, 36% in Cremona, 29% in Lodi and 19% in Brescia). Our results show that policymakers interested in trading off the benefit of containment measures with the cost of shutting down economic activities should therefore tailor their policies to relatively small geographical areas. Despite a robust and significant negative correlation between the first-wave excess mortality and second-wave contagions, we cannot convincingly disentangle whether the effect is due to a sort of herd immunity or to a behavioral response of individuals. Most likely, the two effects are simultaneous and complement each other. We analyzed the correlation between the exposure to the first-wave COVID-19 pandemic in Lombardy and the current second-wave risk of contagion. We report that municipalities more severely hit by the epidemic during the first wave are experiencing fewer cases during the second wave. Although the results should be interpreted cautiously in light of the assumptions and limitations inherent to our approach, our results 9 . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 10, 2020. ; https://doi.org/10. 1101 suggest that policymakers and health authorities should collaborate to design containment measures that are tailored only on small geographic entities. . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 10, 2020. ; https://doi.org/10.1101/2020.11.08.20227934 doi: medRxiv preprint The great influenza: the epic story of the deadliest plague in history Cross-Protection between Successive Waves of the 1918-1919 Influenza Pandemic: Epidemiological Evidence from US Army Camps and from Britain Transmission dynamics reveal the impracticality of COVID-19 herd immunity strategies A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2 Estimating the severity of COVID-19: Evidence from the italian epicenter Second wave covid-19 pandemics in europe: a temporal playbook GMM estimation with cross sectional dependence On two-step estimation of a spatial autoregressive model with autoregressive disturbances and endogenous regressors Asocial capital: Civic culture and social distancing during covid-19 Decreased case fatality rate of covid-19 in the second wave: A study in 53 countries or regions COVID-19: A global and continental overview of the second wave and its (relatively) attenuated case fatality ratio A generalized spatial two-stage least squares procedure for estimating a spatial autoregressive model with autoregressive disturbances First-wave covid-19 transmissibility and severity in china outside hubei after control measures, and second-wave scenario planning: a modelling impact assessment The end of social confinement and covid-19 re-emergence risk SARS-CoV-2 and Lombardy: Testing the impact of the first wave of the pandemic Beware of the second wave of COVID-19. The