key: cord-0716806-frlc6694 authors: Sanchez-Lorenzo, A.; Vaquero-Martinez, J.; Calbo, J.; Wild, M.; Santurtun, A.; Lopez-Bustins, J.-A.; Vaquero, J.-M.; Folini, D.; Anton, M. title: Anomalous atmospheric circulation favored the spread of COVID-19 in Europe date: 2020-05-01 journal: nan DOI: 10.1101/2020.04.25.20079590 sha: df4368d6df1c13661c83da7618e649cabc2a4428 doc_id: 716806 cord_uid: frlc6694 The current pandemic caused by the coronavirus SARS-CoV-2 is having negative health, social and economic consequences worldwide. In Europe, the pandemic started to develop strongly at the end of February and beginning of March 2020. It has subsequently spread over the continent, with special virulence in northern Italy and inland Spain. In this study we show that an unusual persistent anticyclonic situation prevailing in southwestern Europe during February 2020 (i.e. anomalously strong positive phase of the North Atlantic and Arctic Oscillations) could have resulted in favorable conditions, in terms of air temperature and humidity, in Italy and Spain for a quicker spread of the virus compared with the rest of the European countries. It seems plausible that the strong atmospheric stability and associated dry conditions that dominated in these regions may have favored the virus's propagation, by short-range droplet transmission as well as likely by long-range aerosol (airborne) transmission. The world is currently undergoing a pandemic associated with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is a new coronavirus first noticed in late 2019 in the Hubei province, China 1,2 . The virus has a probable bat 3, 4 Respiratory virus infections can be transmitted via direct and indirect contact, or by means of particles (droplets or aerosols) emitted after a cough or sneeze or during conversation by an infected person. The large particles (>5 μm diameter) are referred to as respiratory droplets and tend to settle down quickly on the ground, usually within one meter of distance. The small particles (<5 μm in diameter) are referred to as droplet nuclei and are related to an airborne transmission. These particles can remain suspended in the air for longer periods of time and can reach a longer distance from the origin 6 . These small aerosol particles are inhalable and can penetrate all the way down to the alveolar space in the lungs 7 , where cell receptors for some infectious respiratory viruses are located, including . 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 preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint the angiotensin converting enzyme II (ACE2) used by SARS-CoV-2 to infect the individual 3 . Airborne transmission has been suggested to play a key role in some diseases like tuberculosis or measles, and even in coronaviruses [8] [9] [10] . A recent study has described that the SARS-CoV-2 virus can remain viable at least up to 3 hours in airborne conditions 11 . Respiratory droplets and aerosols loaded with pathogens can reach distances up to 7-8 meters under some specific conditions such as a turbulence gas cloud emitted after a cough of an infected person 12 . A study performed in Wuhan, the capital of the Hubei province, has shown that the SARS-CoV-2 virus could be found in several health care institutions, as well as in some crowded public areas of the city. It also highlights a potential resuspension of the infectious aerosols from the floors or other hard surfaces with the walking and movement of people 13 . Another study has also shown evidence of potential airborne transmission in a health care institution 14 . Recent studies have pointed out a main role of temperature and humidity in the spread of COVID-19. Warm conditions and wet atmospheres tend to reduce the transmission of the disease [15] [16] [17] [18] [19] [20] [21] [22] . For example, it has also been pointed out that the main first outbreaks worldwide occurred during periods with temperatures around 5-11ºC, never falling below 0ºC, and specific humidity of 3-6 g/kg aproximately 18 . The first major outbreak in Europe was reported in Northern Italy in late February 2020. Following that, several major cases have been reported in Spain, Switzerland and France in early March, with a subsequent spread over many parts of Europe. At present (28 th March 2020) Italy and Spain are still the two main contributors of cases and deaths in the . 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) The copyright holder for this preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint continent, with major health, political and socio-economic implications. The main hypothesis of this work is that the atmospheric circulation pattern in February 2020 has helped to shape the spatial pattern of the outbreak of the disease in Europe. The main atmospheric circulation pattern during February 2020 was characterized by an anomalous anticyclonic system over the western Mediterranean basin, centered between Spain and Italy, and lower pressures over Northern Europe centered over the Northern Sea and Iceland ( Figure 1 , Figure S1 ). This spatial configuration represents the well-known North Atlantic Oscillation (NAO) 23,24 in its positive phase, which is the teleconnection pattern linked to dry conditions in southern Europe whereas the opposite occurs in northern Europe 25 . Figure S2 show maps for February 2020 for several meteorological fields that provide clear evidence of the stable atmospheric circulation in southern Europe, with a tendency towards very dry (i.e., lack of precipitation) and calm conditions. As suggested in an earlier analysis 18 , the SARS-CoV-2 virus seems to be transmitted most effectively in dry conditions with daily mean air temperatures between around 5ºC and 11ºC, which are the conditions shown in Figure 2 for the major part of Italy and Spain. By contrast, northern Europe has experienced mainly wet and windy conditions due to an anomalous strong westerly circulation that is linked to rainy conditions. These spatial patterns fit with the well-known climate features associated over Europe during positive phases of the NAO 26 . The Arctic Oscillation (AO), which is a teleconnection pattern very much linked to NAO, . 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) The copyright holder for this preprint this version posted May 1, 2020. We argue that this spatial configuration of the atmospheric circulation might have played a key role in the modulation of the early spread of the COVID-19 outbreaks over Europe. It is known that some cases were reported already in mid-January in France, with subsequent cases in Germany and other countries 29 . Thus, the SARS-CoV-2 virus was already in Europe in early 2020, but it may only have started to extend rapidly when suitable atmospheric conditions for its spread were reached. It is possible that these proper conditions were met in February, mainly in Italy and Spain, due to the anticyclonic conditions previously mentioned. The link between the COVID-19 spread and atmospheric circulation has been tested as follows. We have extracted the monthly anomalies of sea level pressure (SLP) and 500 hPa geopotential height for February 2020 over each grid point of the 15 capitals of the European countries ( Figure S4 ) with the highest number of COVID-19 cases reported so far (see Data and Methods). Figure 3 (top) shows that there is a statistically significant . 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 preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint (R 2 =0.481, p<0.05) second order polynomial fit between the anomalies of the 500 hPa and the total cases per population. Italy, Spain, and Switzerland, which are the only countries with more than 1,000 cases/million inhabitants in our dataset, clustered together in regions with very large positive anomalies of 500 hPa geopotential heights. For the total number of deaths the fit is also statistically significant for a second order polynomial regression (R 2 =0.50, p<0.05), and it shows clearly how Italy and Spain are out of scale compared to the rest of the European countries. Similar results are obtained using SLP fields (not shown). These results evidence that it seems plausible that the positive phase of the NAO, and the atmospheric conditions associated with it, provided optimal conditions for the spread of the COVID-19 in southern countries like Spain and Italy, where both the start and the most severe impacts of the outbreak in Europe were located. To test this hypothesis further we have also analyzed the COVID-19 and meteorological data within Spain (see Data and Methods, Figure S5 ). The results show that mean temperature and specific humidity variables have the strong relation with COVID and fit with an exponential function ( Figure 4 ). They indicate that lower mean temperatures (i.e., average of around 8-11ºC) and lower specific humidity (e.g., <6 g/Kg) conditions are related to a higher number of cases and deaths in Spain. Nevertheless, it is worth mentioning that both meteorological variables are highly correlated (R 2 =0.838, p<0.05) and are not independent of each other. The temperatures as low as 8-10ºC are only reached in a few regions such as Madrid, Navarra, La Rioja, Aragon, Castille and Leon and Castilla-La Mancha. These areas are mainly located in inland Spain where drier conditions were reported the weeks before the outbreak. The rest of Spain experienced higher temperatures and consequently were out of the 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. The copyright holder for this preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint of higher potential for the spread of the virus, as reported so far in the literature [15] [16] [17] [18] [19] [20] [21] [22] . In addition, higher levels of humidity also limit the impact of the disease, and therefore the coastal areas seem to benefit from lower rates of infection. Thus, in the southern regions of Spain (all of them with more than 13ºC and higher levels of specific humidity) we found lower rates of infection and deceases. This is in line with the spatial pattern in Italy, with the most (least) affected regions by COVID-19 mainly located in the North (South). In contrast, when the whole of Europe is considered on a country by country basis (see above and Figure 3 ), we find the opposite, a clear gradient with more severity from North to South as commented previously. The spatial pattern of COVID-19 described above has some intriguing resemblances with the 1918 influenza pandemic, which is the latest deadly pandemic in modern history of Europe. The excess-mortality rates across Europe in the 1918 flu also showed a clear northsouth gradient, with a higher mortality in southern European countries (i.e., Portugal, Spain or Italy) as compared to northern regions, an aspect that is not explained by socio-economic or health factors 30 . In Spain, a south-north gradient is also reported in the 1918 flu after controlling for demographic factors 31 . The central and northern regions of Spain experienced higher rates of mortality, and this has been suggested to be linked to more favorable climate conditions for influenza transmission as compared to the southern regions 31 . Interestingly, the SLP anomalies of the months before the major wave of this pandemic (which occurred in October-November 1918) shows a clear south-north dipole with positive anomalies in southern Europe centered over the Mediterranean, and negative ones in northern Europe ( Figure S6 ). In other words, the NAO was also in its positive phase just before the major outbreak of the 1918 influenza pandemic. This resembles the spatial . 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 preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint patterns described above for the current COVID-19 outbreak, both in terms of the spatial distribution of the mortality of the pandemic over Europe as well as in prevailing atmospheric circulation conditions before the major outbreak. These intriguing coincidences need further research in order to better understand the spatial and temporal distribution of large respiratory-origin pandemics over Europe. Taking into account these results, we claim that the major initial outbreaks of COVID-19 in is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2020. Although the outbreak of a pandemic is controlled by a high number of biological, health, political, social, economic and environmental factors, with complex and non-linear interrelationships between them, the results of this study indicate that an anomalous atmospheric circulation may explain why the COVID-19 outbreak in Europe developed more easily (or faster) in the south-west (mainly north of Italy and inland of Spain). Specifically, the extreme positive phase of the AO and NAO during February 2020 could have modulated the beginning of the major outbreaks of COVID-19 in Europe. This detected anomalous atmospheric pattern, which produces dry conditions over southwestern Europe, may have provided optimal meteorological conditions for the virus propagation. In the context of anthropogenic climate change, it has been shown that in future emissions scenarios a poleward expansion of the Hadley cell is expected 41 , which in turn is in line with a tendency to increase the frequency of positive phases of the NAO 42 ( Figure S8 ). This should be taken into account for planning against future epidemics and pandemics that arise from respiratory viruses. Interestingly, the conditions during the last major pandemic experienced in Europe (the Spanish flu in 1918), seem to resemble the current spatial pattern of affectation with more cases in the South of Europe as compared to the North. Equally, the dominant atmospheric situation was strongly affected by anticyclonic (cyclonic) conditions in the South (North) of . 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 preprint this version posted May 1, 2020. . 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 preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint . 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 preprint this version posted May 1, 2020. . 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 preprint this version posted May 1, 2020. . 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 preprint this version posted May 1, 2020. NCEP/NCAR 1 , ERA5 2 and ERA20C 3 atmospheric data are used in this manuscript. The maps and data have been retrieved by using the tools and websites referenced in the main text, and more details about the spatial and temporal resolution, vertical levels, assimilation schemes, etc. can be consulted in their references. In brief, an atmospheric reanalysis like those used here is a climate data assimilation project which aims to assimilate historical atmospheric observational data spanning an extended period, using a single consistent assimilation scheme throughout, with the aim of providing continuous gridded data for the whole globe. The Artic Oscillation (AO) index has been extracted from the Climate Prediction Center of the National Oceanic and Atmospheric Administration (NOAA). The AO index is constructed by projecting the daily 1000 hPa height anomalies poleward of 20°N onto the loading pattern of the AO, this latter being defined as the leading mode of Empirical Orthogonal Function (EOF) analysis of monthly mean 1000 hPa height. More details and data: https://www.cpc.ncep.noaa.gov/products/precip/CWlink/daily_ao_index/ao.shtml. COVID-19 data on country basis were obtained on March 26 th , 2020 from the website https://www.worldometers.info/coronavirus/, which it is mainly based on the data provided by the Coronavirus COVID-19 Global Cases by the Center for Systems Science and Engineering (CSSE) at the Johns Hopkins University. Data from Spain on regional scale . 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 preprint this version posted May 1, 2020. we have extracted the monthly anomalies of sea level pressure (SLP) and 500 hPa geopotential height for February 2020 over each grid point of the 15 capitals of the European countries. We have selected the SLP and 500 hPa fields in order to summarize the meteorological conditions over each location, as it is known that several meteorological variables can be involved in the transmission of respiratory viruses 4, 5 . With this approach we also avoid the lack of properly updated data for all potential meteorological variables involved in the COVID-19 spread, which needs further research as soon as the pandemic ends and a more reliable and complete database of both COVID-19 impact and meteorological data can be compiled 6 . For Spain, several meteorological variables with high-quality records were obtained from the Spanish Meteorology Agency (AEMET) based on surface observations for each of the capital cities of the provinces inside each autonomous region. Specifically, monthly averages for February 2020 of 2-m temperature, 2-m maximum temperature, 2-m minimum temperature (°C), air pressure (hPa), wind speed (km h −1 ), specific humidity (g Kg -1 ), relative humidity (%), total precipitation (mm) and days of more than 1 mm of precipitation. An arithmetic average has been calculated for the autonomous regions with more than one province. We have checked Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations for two future scenarios (RCP4.5 and RCP8.5) at the end of the 21th century. . 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) The copyright holder for this preprint this version posted May 1, 2020. . 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) The copyright holder for this preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint Clinical features of patients infected with 2019 novel coronavirus in Wuhan , China A pneumonia outbreak associated with a new coronavirus of probable bat origin Identifying SARS-CoV-2 related coronaviruses in Malayan pangolins The role of particle size in aerosolised pathogen transmission : A review COVID-19): how COVID-19 spreads Recognition of aerosol transmission of infectious agents : a commentary Evidence of Airborne Transmission of the Severe Acute Respiratory Syndrome Virus Aerosol and Surface Stability of SARS-CoV-2 as Compared with SARS-CoV-1 Turbulent Gas Clouds and Respiratory Pathogen Emissions: Potential Implications for Reducing Transmission of COVID-19 Aerodynamic Characteristics and RNA Concentration of SARS-CoV-2 Aerosol in Wuhan Hospitals during COVID-19 Outbreak Transmission Potential of SARS-CoV-2 in Viral Shedding Observed at the University of The role of absolute humidity on transmission rates of the COVID-19 outbreak Role of temperature and humidity in the modulation of the doubling time of COVID-19 cases. under Rev The NCEP / NCAR 40-Year Reanalysis Project ERA5: Fifth generation of ECMWF atmospheric reanalyses of the global climate Influenza Virus Transmission Is Dependent on Relative Humidity and Temperature The effects of weather and climate on the seasonality of influenza: What we know and what we need to know . 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 preprint this version posted May 1, 2020. . 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)The copyright holder for this preprint this version posted May 1, 2020. . https://doi.org/10.1101/2020.04.25.20079590 doi: medRxiv preprint