key: cord-0800963-3op1peh3 authors: Björk, Jonas; Mattisson, Kristoffer; Ahlbom, Anders title: Impact of winter holiday and government responses on mortality in Europe during the first wave of the COVID-19 pandemic date: 2021-02-24 journal: Eur J Public Health DOI: 10.1093/eurpub/ckab017 sha: 423dbe5832612d27354ffde28d74a3ee769bdc60 doc_id: 800963 cord_uid: 3op1peh3 BACKGROUND: This aggregated population study investigated the impact of the seemingly quasi-randomly assigned school winter holiday in week 6–10 (February to early March) on excess mortality in 219 European regions (11 countries) during the COVID-19 pandemic in the spring 2020. A secondary aim was to evaluate the impact of government responses to the early inflow of infected cases. METHODS: Data on government responses week 8–14 were obtained from the Oxford COVID-19 Government Response Tracker. Regional data on total all-cause mortality during week 14–23 in 2020 were retrieved from Eurostat and national statistical agencies and compared with the average mortality during same period 2015–2019. Variance-weighted least square regression was used with mortality difference as dependent variable with adjustment for country, population density and age distribution. RESULTS: Being a region with winter holiday exclusively in week 9 was in the adjusted analysis associated with 16 weekly excess deaths (95% confidence interval 13 to 20) per million inhabitants during week 14–23, which corresponds to 38% of the excess mortality in these regions. A more stringent response implemented in week 11, corresponding to 10 additional units on the 0–100 ordinal scale, was associated with 20 fewer weekly deaths (95% confidence interval 18 to 22) per million inhabitants. CONCLUSIONS: Winter holiday in week 9 was an amplifying event that contributed importantly to the excess mortality observed in the study regions during the spring 2020. Timely government responses to the resulting early inflow of cases reduced the excess in mortality. Early cases of COVID-19 known to be infected in Europe (21 cases by 21 February, week 8, 2020) were linked to two different clusters, one in Bavaria, Germany, and one in Haute-Savoie, France. 1 On 22 February the Italian authorities reported clusters of cases in Lombardy and additional cases from three other regions in northern Italy. 2 These initial clusters were all located in close proximity to the Alps, the most extensive mountain range in Europe. Late detection of the index cases generally delayed isolation of further local cases in the region. 1 Additionally, the large number of tourists at ski resorts in the Alps during the school winter holidays in February and March most likely contributed importantly to the further spread of the virus across Europe. This assertion is supported by genetic analyses of SARS-CoV-2 in for example Denmark and Iceland, which have associated specific haplotypes with travel to ski resort areas in Austria (haplotype A2a2a) and Italy (A2a1). 3, 4 Travelling during public holidays has further been found to be important determinants of the COVID-19 transmissions out from the virus epicentres in Wuhan, China, and Qom, Iran. [5] [6] [7] [8] In the present study we focus on 11 countries in Europe that all have school winter holidays and were not part of the initial epicentre: Germany, France, UK, Belgium, Netherlands, Luxemburg, Sweden, Denmark, Norway, Finland and Iceland. Striking regional differences in excess mortality during March -June 2020 have been observed for these countries ( Figure 1A ). 9 It is not known, however, to what extent the exact week of the school winter holiday, which varied from week 6 (early February) to week 10 (early March) across regions, 10 (Figure 1B) , has contributed to the variation in mortality. We hypothesized that the exact week of the winter holiday was not only important for the initial spread of the disease, but also essentially quasi-random with respect to structural determinants of the further development of the pandemic and its effect on the mortality across regions. http://mc.manuscriptcentral.com/ejph The governments in these 11 European countries reacted differently and with different timing in the implementation of non-pharmaceutical interventions at the population level as a response to the COVID-19 pandemic during the spring 2020. 11 As example of the differences in response, France and Germany, reacted stringently and early, UK reacted stringent but with a delay, while Sweden maintained a less stringent reaction throughout this period ( Figure S1 ). Previous studies have suggested that a timely response in relation to when initial infections occurred had critical and long-term effects. 12, 13 However, so far no systematic effort has been made to quantify this effect on mortality while taking the influence of amplifying events or other population-level factors into account. The overall aim of the present ecological study of 219 regions in 11 European countries was to exploit the effect of the seemingly quasi-randomly assigned school winter holiday week on the excess mortality during spring 2020. As a secondary aim, the impact of stringency and timing of the initial government responses was evaluated. The present study includes population data 2015 -2020 for 219 regions in 11 European countries (B, D, DK, F, FIN, ISL, LUX, N, NL, S and UK) that generally all have school winter holidays for at least one week. An underlying assumption when defining the study area was that a significant part of the population having winter holiday in the included countries travel to the Alps for skiing, 14 including the areas in Austria and northern Italy affected early during the COVID-19 pandemic. Regional data for the present study were originally extracted from European Data Journalism Network, 9 and subsequently verified and updated by data from Eurostat and national statistical agencies. Most data were gathered on the NUTS (Nomenclature des Unités Territoriales Statistiques) level 3 (169 regions), but Germany and UK only had mortality data Data on school winter holiday weeks for each study region were directly accessible from a Swiss website 10 , and verified for a subsample. Nine regions in Germany either did not have a full week winter holiday (n = 7), or had the holiday already week 6 (n = 2). All other 210 regions had their winter holiday between week 7 and 10 ( Figure 1B and Table S1 ). France had two-week winter holidays (week 7 + 8, 8 + 9 or 9 + 10). The proportion of the study population exposed to winter holidays the different weeks in each country is presented in Table S2 . Scotland (UK NUTS 1) was not included in the study as its winter holiday week varied on finer geographical levels than available mortality data. Rest of UK had their winter holiday solely in week 8 and Belgium solely in week 9. Daily country-level data on initial non-pharmaceutical interventions at the population level were obtained from the Oxford COVID-19 Government Response Tracker. 11 Specifically, we used their stringency index (0 -100) which is an average of 9 ordinal items, school closing, work place closing, cancelation of public events, restrictions on gatherings, closing of public transports, stay at home requirements, restrictions on internal movement, international travel controls and public information campaigns. We calculated weekly averages for each country during week 8 (17 -23 Feb) to week 14 (30 Mar -5 Apr) based on data for Monday through Friday ( Figure S1 ). No specific regional data on government response were available, which implies that all regions within a country were assigned the same value for the stringency index a given week. We considered this appropriate given that the initial responses were mainly acting on national rather than on regional levels. Statistical analyses were conducted in Stata SE 14.2 (Stata Corp.) and IBM SPSS Statistics 25 (SPSS Corp.). We calculated the mortality difference (MD) for week 14 to 23 for each region , The weighted analysis implies that observations for regions with larger population sizes, and thus more precise estimates of MD, are more influential than those for smaller regions with fewer inhabitants. Holiday week was used as categorical in the analysis with regions with no winter holiday or holiday in week 6 (all in Germany) as a collapsed reference category due to small numbers. The two-week holidays in France were classified by placing one half of the deaths and person-weeks in the first week and the other half in the second week. The regression analysis for the association between winter holiday week and MD was presented both unadjusted and with adjustment for country (categorical without reference category), the population density of the region on the log-scale, the proportions of individuals below 20 (most of them exposed to school holidays) and the proportion above 70 years of age (as an indicator of the vulnerability of the population). In separate linear regression analyses, we investigated the association between countrylevel stringency of the government response (index 0 -100) during week 8 -14, one week at a time, and the regional MD for the entire follow up. The dependent variable, holiday-, density and age-adjusted MD, and its variance-weights were obtained from the previous analytical step. We also estimated the percentage of the observed country-level differences in mortality that could be attributed to regional differences in winter holiday week, population density and government response stringency, respectively. A total of 586 977 deaths occurred in the study area during week 14 -23 Elevated MDs were observed in non-French regions having their winter holidays in week 8 or 9, and for the French regions irrespectively of the timing of winter holiday ( Figure 2 and Table S2 ). There were differences in population density, age distribution and country affiliation between regions that had their winter holiday early (week 6 -7), in the middle (week 8 -9) or late in the season (week 10; Table S1 ). In particular, regions having winter holiday week 9 or 10 were generally less dense than regions having holiday earlier or not at all. Thus, adjustments for population factors led to marked changes in the estimate effects associated with winter holiday week ( Model 4) . Population density of the region was clearly associated with mortality. A tenfold higher population density (e.g. 1 000 rather than 100 inhabitants per km 2 ) was associated with 12.5 and 20.2 excess deaths per million habitants and week in the non-French and French regions, respectively. Age distribution of the population was essentially unrelated to the regional excess mortality. Overall, an estimated 6.1% (95% CI 4.8 to 7.5%) of the excess deaths in the study area was attributed to having winter holiday in week 9, but this proportion varied substantially across different countries. Large number of excess deaths was attributed to winter holiday week 9 in the regions of Germany (Bavaria), Belgium (all 11 regions), Netherlands (4 regions) and Sweden (Stockholm and 5 other regions; Table 2 ). By contrast, the winter holiday week did not seem to explain any of the excess mortality in UK or France. As a comparison, 19% of the excess mortality (95% CI 14 to 24) was attributed to population density above 100 persons per km 2 . Population density explained large proportions of the excess mortality in Germany (66%) and France (27%) but only a minor proportion in for example Sweden (4%; Table 2 ). Timing and stringency of the government response in relation to excess mortality was illustrated for all 219 regions and the subset of 38 regions with winter holiday in week 9 in Figure S2 . Response data were analysed further on the country-level after adjustment for winter holiday week, population density and age distribution of the regions (Table S4) Overall, 20% of the excess mortality (95% CI 17 to 22%) in the study area was attributed to response stringency in week 11 that was below average (Table 2) . For Sweden, we estimated that 55% of the excess mortality (95% CI 49 to 61%) was associated with a less stringent response. A salient finding of the present study was that the timing of winter holiday contributed importantly to the excess mortality observed in the investigated 219 European regions during the spring 2020. The winter holiday was held in week 9 in 38 regions constituting one fifth of the study population, and was in these regions estimated to explain 38% of the excess mortality. Another important finding was that a timely government response to the resulting early inflow of cases was associated with lower excess mortality. By contrast, the government responses early in the pandemic (week 8) or late when the transmission was already widespread (week 13) did not seem to explain any of the observed differences in excess mortality. No effect of the winter holiday on the mortality was discerned for the French regions despite that the majority (75 out of 96 regions) had school holiday that included week 9. A possible explanation could be that the French population do not go abroad on winter holiday to the same extent as other European countries. 15 The proportion of trips in Q1 (first quarter) with at least one night stay with an international destination, was only 10% in France, compared to for example 73% in Belgium, 29% in Denmark and 27% in Germany and Sweden according to the most recent European quarterly travel statistics in 2011. 15 . Not only the relative but also the absolute number of persons travelling abroad in relation to population size was low in France compared with the other countries. 15 UK citizens travel abroad during the winter to a large extent (36% of all trips in 2011 Q1) but had school holiday already in week 8, i.e. one week earlier than the peak inflow of cases returning from the Alps. The Paris is large also during the winter, 14 which thus seems as a more likely explanation for the for the initial spread of the pandemic in France and UK than returning travellers from the Alps. Having winter holiday in week 10 was not associated with excess mortality, which was unexpected as the initial epicentre in the Alps had developed further by then. One possible explanation could be that the week 10 regions were mostly rural and located in northern part of the study region with fewer travellers to the Alps. Comparisons between countries and regions are essential in order to improve understanding on why some have markedly higher burdens of COVID-19 mortality than others, 16 our study was that we were able to estimate the effect of government response while controlling for the regional differences in initial spread caused by the timing of the winter holiday, and for differences in population density and age structure. The results confirmed the importance of a timely response to the COVID-19 pandemic suggested by previous real-world studies. 12, 13, 19 Belgium, Sweden and UK stood out as the three countries in our study area Cross-national differences in socioeconomic conditions has been put forward as a potential explanation for some of the observed differences in excess mortality during the COVID-19 pandemic. 20 A population-based study from Sweden with follow up until early May, 2020, found similar socioeconomic gradients in COVID-19 mortality as for mortality from all other causes. 21 Additionally, an elevated mortality was observed among immigrants from low-and middle-income countries also after adjusting for socioeconomic characteristics. There are differences in this socioeconomic gradient in health across the Nordic countries, but neither the socioeconomic structure nor the strength of the gradient appear to be particularly disadvantageous in Sweden. 22 Thus, higher population vulnerability is not likely to be an important explanation for the excess mortality observed in Sweden compared to its neighbours. Our study had several limitations. The use of aggregated data means that results can be subject to ecological inference fallacy, i.e. that the observed associations are not necessarily reflecting true or correctly estimated associations on the individual level. 16 The ecological analysis is generally sensitive to choice of adjustment factors. We decided to include adjustment for country, which effectively implies that associations between winter holiday week and excess mortality in each country are pooled together in the regression analysis. But it also means that only countries where the winter holiday week differed across regions contribute to the estimated exposure effect. We used excess mortality as outcome measure to avoid bias due to differences in registration of covid-19 deaths in different countries, but this also means that we cannot rule out that some of the excess was unrelated to the pandemic. The calculation of the expected mortality in a region was only adjusted for changes in population size, and did not take time trends, yearly fluctuations (e.g. mean temperature a given year or severity of the flu season) or other population changes into account. 13 However, the quasi-experimental nature of the study makes it unlikely that socioeconomic conditions, yearly fluctuations or changes in population structure were systematically associated with the timing of the winter holiday. A further limitation was that we that we only estimated the effect of the winter holiday for each region separately and did not consider spill-over effects to neighbouring regions from for example commuting or regional travelling. In this respect, it is therefore likely that our estimates represent the lower bound of the winter holiday effect. It should also be noted that government response was only evaluated in the initial phase of the pandemic. Lockdown measures and other restrictions that last for longer time periods may have substantial adverse health, psychosocial and financial consequences, 17, 23-25 which were not assessed in our study. In conclusion, the timing of the winter holiday contributed importantly to variation in the excess mortality observed in the study area during the spring 2020. It is imperative that authorities respond quickly when such amplifying events are identified in order to limit their public health impact. JB conceived the study with a design that was further developed together with AA and KM. JB did the data acquisition and analysis with assistance regarding the mapping from KM. JB drafted the first version of the manuscript, where after all authors revised it critically. All authors have read and approved the final manuscript and agree to be accountable for all aspects of the work. The present study was based on secondary, aggregated data that are hosted and maintained by Outbreak of novel coronavirus disease 2019 (COVID19): situation in Italy. Stockholm: European Centre for Disease Prevention and Control SARS-CoV-2 transmission chains from genetic data: a Danish case study Spread of SARS-CoV-2 in the Icelandic Population Travel rush during Chinese Spring Festival and the 2019-nCoV Prediction of Epidemic Spread of the 2019 Novel Coronavirus Driven by Spring Festival Transportation in China: A Population-Based Study Nowruz travelers and the COVID-19 pandemic in Iran The association between the incidence of COVID-19 and the distance from the virus epicenter in Iran A fraction of European regions account for a majority of COVID-19 deaths: European Data Journalism Network Variation in government responses to covid-19 COVID-19 pandemic-related lockdown: response time is more important than its strictness Magnitude, demographics and dynamics of the effect of the first wave of the COVID-19 pandemic on all-cause mortality in 21 industrialized countries Analysing spatiotemporal patterns of tourism in Europe at high-resolution with conventional and big data sources Number of trips -quarterly data: European Union Nussbaumer-Streit B, Mayr V, Dobrescu AI et al. Quarantine alone or in combination with other public health measures to control COVID-19: a rapid review Impact assessment of non-pharmaceutical interventions against coronavirus disease 2019 and influenza in Hong Kong: an observational study Cross-Country Comparisons of Covid-19: Policy, Politics and the Price of Life A population-based cohort study of socio-demographic risk factors for COVID-19 deaths in Sweden Income-related health inequalities in the Nordic countries: examining the role of education, occupational class, and age Mitigating the wider health effects of covid-19 pandemic response The psychological impact of quarantine and how to reduce it: rapid review of the evidence Lockdown, quarantine measures, and social distancing: Associations with depression, anxiety and distress at the beginning of the COVID-19 pandemic among adults from Germany Log10(Density) centered at 100 persons per km 2 , proportion below 20 and above 70 at 20% and 15% All-cause mortality difference (MD) per week and million inhabitants in 219 European regions, contrasting 2020 week 14 -23 with the same period in 2015 -2019 (2016 -2019 for Germany and Netherlands), and stratified by the timing of the school winter holiday in each region (week 6 to 10). Regions with no winter holiday were collapsed with week 6. The French regions with two weeks winter holiday were placed in between the first and second week We want to acknowledge the European Data Journalism Network who collated and published regional data from Eurostat and national statistical agencies that formed the basis for our study. The authors have declared no conflicts of interest.