Evidence and magnitude of the effects of meteorological changes on SARS-CoV-2 transmission

Importance Intensity and duration of the COVID-19 pandemic, and planning required to balance concerns of saving lives and avoiding economic collapse, could depend significantly on whether SARS-CoV-2 transmission is sensitive to seasonal changes. Objective Hypothesis is that increasing temperature results in reduced SARS CoV-2 transmission and may help slow the increase of cases over time. Setting Fifty representative Northern Hemisphere countries meeting specific criteria had sufficient COVID-19 case and meteorological data for analysis. Methods Regression was used to find the relationship between the log of number of COVID-19 cases and temperature over time in 50 representative countries. To summarize the day-day variability, and reduce dimensionality, we selected a robust measure, Coefficient of Time (CT), for each location. The resulting regression coefficients were then used in a multivariable regression against meteorological, country-level and demographic covariates. Results Median minimum daily temperature showed the strongest correlation with the reciprocal of CT (which can be considered as a rate associated with doubling time) for confirmed cases (adjusted R2 = 0.610, p = 1.45E-06). A similar correlation was found using median daily dewpoint, which was highly colinear with temperature, and therefore was not used in the analysis. The correlation between minimum median temperature and the rate of increase of the log of confirmed cases was 47% and 45% greater than for cases of death and recovered cases of COVID-19, respectively. This suggests the primary influence of temperature is on SARS-CoV-2 transmission more than COVID-19 morbidity. Based on the correlation between temperature and the rate of increase in COVID-19, it can be estimated that, between the range of 30 to 100 degrees Fahrenheit, a one degree increase is associated with a 1% decrease—and a one degree decrease could be associated with a 3.7% increase—in the rate of increase of the log of daily confirmed cases. This model of the effect of decreasing temperatures can only be verified over time as the pandemic proceeds through colder months. Conclusions The results suggest that boreal summer months are associated with slower rates of COVID-19 transmission, consistent with the behavior of a seasonal respiratory virus. Knowledge of COVID-19 seasonality could prove useful in local planning for phased reductions social interventions and help to prepare for the timing of possible pandemic resurgence during cooler months.


Introduction
Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), is a newly-identified enveloped, non-segmented, positive sense RNA virus. It is responsible for Coronavirus Disease 2019 (COVID-2019), which was designated a pandemic by the World Health Organization (WHO) on March 11, 2020 [1]. SARS-CoV-2 belongs to a large family of human corona viruses such as MERS-CoV and SARS-CoV-1, which are respiratory pathogens associated with a range of respiratory and non-respiratory outcomes. Like other respiratory infectious agents, most coronaviruses display marked seasonality [2]. Seasonality in respiratory viruses is characterized by increased transmission in cooler, less humid months and decreased transmission in warmer and more humid months. It has been demonstrated in the laboratory that for most respiratory viruses, including Influenza, SARS-CoV-1 and SARS-CoV-2, the stability of the envelope surrounding and protecting these viruses is temperature sensitive-degrading more rapidly at higher temperatures and showing increased stability at cooler temperatures [3]. In addition to temperature, the transmission of respiratory viruses demonstrating seasonality are also impacted by humidity, where lower humidity leads to smaller droplets that are more volatile and capable of lodging deeper in the lung [3]. Usually, rather than humidity, dewpoint is substituted instead in meteorological reports. Dew point is the temperature, at constant pressure, to which the air needs to be cooled in order to achieve a 100% relative humidity [4]. At that point the air cannot hold more water in the gas form. SARS-CoV-1, which shows seasonal changes with decreased rates of transmission with elevated temperatures, is genetically 80% identical to SARS-CoV-2 [4,5].
Because SAR-CoV-2 is a recently identified human pathogen, seasonal variations in its transmission have not become evident, though there is much speculation that seasonal change with the boreal summer might decrease infection rates of SARS-CoV-2 and help flatten the epidemic curve of COVID-19 [6]. Several studies, varying in their geographical sampling,

Funding:
The authors received no specific funding for this work. MAAR is an employee of Pupa but contributed to this work independently of this affiliation. Pupa had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests:
The authors of the study have read the journal's policy, and have the following competing interests to declare: MAAR is an employee of Pupa but contributed to this work independently of this affiliation. AK is currently employed by MyMD Pharmaceuticals; however, this affiliation was not held at the time this study was conducted. This does not alter our adherence to PLOS ONE policies on sharing data and materials. There are no patents, products in development or marketed products associated with this research to declare. The authors have declared that no other competing interests exist.
methodology and findings, have endeavored to shed light on this subject. With respect to geography, four out of seven of the peer-reviewed studies that attempted to establish whether SARS-CoV-2 has seasonal changes in its transmissibility examined rates in China [7][8][9][10], two others restricted themselves to either Japan or France [11,12], and two examined a collection of cities or countries around the world [7,13]. Methodologies included the use of generalized additive model, poisson regression analysis, and univariate and multiple regression modeling [7,[9][10][11]14]. The findings concerning whether SARS-CoV-2 is likely to display seasonality in its rates of transmission varied widely, including: a) two studies reported an asymmetric impact of temperature (one being biphasic at 50˚F, and another demonstrating effects only below 37.4˚F [8,15], b), one showed no influence of temperature on COVID-19 rates [9], c) two showed inconclusive results with respect to the likelihood of seasonality of SARS-CoV-2 [8,13], d) one showed distribution of substantial community outbreaks of COVID-19 along restricted latitude, temperature, and humidity measurement [13] and e) one showed a sensitivity to temperature that was of such large magnitude (i.e. a 50% decrease in rate of transmission with every 1 degree increase in average temperature) that it seemed incongruous with observed variations across the globe [14]. A better understanding of the seasonal dependence of SARS-CoV-2 infections, and the identification of environmental conditions that regulate its spread, would be helpful in assessing the differential impact of the virus in different geographic locations under various meteorological conditions. Such knowledge would also help in predicting the timing and degree of potential slowing of transmission in the summer and resurgence of the virus during winter [16].
To understand the seasonal dependence of SARS-CoV-2 infections in the absence of yearlong incidence data, we examined the effect of atmospheric temperature, the main variable with season [16], on SARS CoV-2 transmissibility using currently available data on the rates of SARS-CoV-2 infection during winter and spring from fifty representative countries in the Northern Hemisphere.

Geographical sampling criteria
Geographical changes in cumulative COVID-19 confirmed deaths, and recovered cases from 1/22/20 through 4/6/20 were estimated from databases made available by the Johns Hopkins University [17]. The start date of January 22, 2020 was chosen because that is when the database started tracking cases. For analysis, we selected only the Northern Hemisphere countries because even though the number of cases have been increasing on both sides of the equator, only the Northern Hemisphere experienced an increase in temperature during the study period. The opposite was true in the Southern Hemisphere. As a result, the impact of temperature on the rate of new cases cannot be investigated simultaneously on pooled data from both hemispheres. We excluded tropical regions because of the minimal change in temperature throughout the year. Italy was excluded because of its "massive underreporting and undertesting of COVID-19," with an estimate of only 2% of total cases being reported based on calculations involving crude case fatality risks [18].
We also excluded from analysis China [19] and Russia [20]. We excluded China and other countries with multiple data sites (for both COVID-19 cases and NOAA meteorological data) within its borders. Each country had differing responses to the increasing rates of COVID-19 within its borders (e.g. the timing and availability of testing, government policy decisions about when and how robustly to respond with non-pharmaceutical interventions, etc). Because each country that was included contributed data representing a single response to COVID-19, the variability in the response would be randomized (some countries more and some less in various aspects of their response). The inclusion of a single country with numerous measures within its borders would risk creating a systematic bias for that country in terms of a potential effect on the rate of transmission of COVID-19. China had 33 provinces for which cases were collected from JHU COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (https://github.com/ CSSEGISandData/COVID-19). Thus, to preserve methodological conformity, countries like China with separate reports of cases from multiple regions were excluded.
There is also an abundance evidence to support the dramatic under-reporting of cases from China before the time the JHU COVID-19 Data Repository database begins on 1/22/20, thus suggesting that trying to base results on non-JHU sources prior to that date would be fraught with large inaccuracies. For example, a recent publication by Krantz et al [21] calculated the ratio of reported COVID-19 cases and contrasted them with model-based predictions of COVID-19 for 8 major countries. The ratio of reported to infected cases for China was between 1:149 to 1:1104. In comparison, France was 1:5 and Germany was 1:4. Based on cremation related information from China, He et al [22] were able to calculate COVID-19 case estimates projected on February 7, 2020 in Wuhan ranged from 305,000 to 1,272,000 for infections that were at least 10 times the official figure of 13,603. The same source of information implied starting time of the outbreak is October 2019, rather than the government reported date of 12/31/19 [22].
Including Russia presented a problem of size. A country of Russia's magnitude had only a single JHU CSSE COVID-19 data entry with a single latitude and longitude, which could clearly not accurately represent the appropriate temperature corresponding to the rate of increase of cases for the entire country. Russia has the largest land area of any country in the world, with vastly differing temperatures across its vast expanse, being comprised of 17 million square kilometers. The largest country included in the 50 analyzed in our manuscript was Algeria, that has a land area of 2.4 million square kilometers, which is 7-fold smaller than Russia. One temperature for such a vast country with a huge range of local climates would not yield accurate data. Thus, Russia was excluded.
Other countries were excluded due to limited access to, or the lack of, reliable data on COVID-19 cases or meteorological data. To enable analysis of cases during consecutive days in a log linear fashion, we established the following criteria to ensure sufficient numbers of confirmed cases: to be included, during time period 1/22/20-4/6/20, countries needed to have a median number of cumulative confirmed cases >1 or the mean number of cumulative confirmed cases >75.
The resulting 50 representative countries (Table 1) had a mean latitude of 38 +/-15 SD, range 1-65 degrees, that covers 70% of the total range (1-90 degrees), with each country being separated from the next by a mean of 1.3 +/-1.7 SD degrees. Countries had a mean longitude of 34+/-44 SD, range -103 to 138, that covered 70% of the total range (-180 to 180 degrees), with each country separated from the next by a mean of 6 +/-13 SD degrees. The resulting wide range of temperatures across the 50 countries was as follows: 1) mean average temperature 51.1 +/-16.9˚F, median 44.9˚F, range 22.6-85.7˚F, and 2) mean minimal temperature 42.0 +/-16.3˚F, median 36.6˚F, range 15.5-77.4˚F.

Meteorological variables
Daily weather data for each of the 50 countries were obtained from National Oceanic and Atmospheric Administration (NOAA) [23]. Daily minimum temperatures (Tmin) were used because they appear to be more important for virus transmissibility [24]. We used Dew Point (DP) as a measure of humidity because it is an absolute assessment of water vapor content that does not depend upon temperature, and it is the only country meteorological data related to humidity that is tracked by NOAA. The daily Tmin and DP for each country were obtained from NOAA station locations closest to the specific latitude and longitude that the JHU database listed for each country. The medians of the daily Tmin and daily DP were derived for each country for whatever portion of 1/22/20 through 4/6/20 that they had cases of confirmed COVID-19 reported (Table 1). Median was chosen because it is less affected by outliers and skewed data than the mean.

Software packages used in analyses
Statistical analysis was performed with StatPlus software for the Mac [25]. The graphing software employed was GraphPad Prism 8 for the Mac [26].

Corrections applied to account for multiple comparisons
The unit variable was country. At the unit level, dependent and independent variables were evaluated. There were three correlated outcome measures: confirmed, deaths, or recovered cases. For each of these three outcomes, an aggregate measure CT or CT -1 was calculated and then associated with multiple independent variables. Due to evaluating three outcomes in the same population, a multiple adjustment was required. Two simple and conservative approaches that have been used include the following: 1. to compare observed p values (p) with adjusted alpha (= alpha/3) or 2. inflate the observed p values by 3 (p � 3) and compare at fixed alpha level [27,28] We adopted the latter approach. Also note that because CT and CT -1 are one-to-one transformations, there was no need to adjust for using CT or CT -1 as an dependent variable.

Relationship of CT to doubling time (Td)
To clarify the relative magnitude of our findings, we explored the relationship between CT and doubling time (T d ), a standard measure of rate of viral transmission [29]. Doubling times are the intervals at which the cumulative incidence of cases doubles. Under conditions of exponential, with a constant growth rate r, the doubling time equals (ln 2)/r [30]. Thus, an increase in the doubling time represents a decreased rate of transmission.
To calculate T d from our data obtained in our first regression (see Results section for details), we took the following steps: to begin, we used our first regression of log of confirmed cases vs time to interpolate on which days the number of cases doubled for each country, utilizing a centered second order polynomial model to fit the data (mean R 2 = 0.958, Standard Deviation = +/-0.025). Because we were using log of confirmed cases, we converted doubling of cases to log form (e. g. 2, 4, 8, 16, 32, 64 corresponds to 0.301, 0.602, 0.903, 1.204, 1.505, 1.806). Since the period of time we measured cases and temperature was 75 days long (1/22/20-4/6/ 20), we included doubling times up to day 64. We next calculated the average T d (i.e. number of days) between each interval corresponding to a doubling of confirmed cases. We then performed a linear regression on T d and CT, and T d and Tmin Med, with the following results:

Analysis of the impact of changing temperature on the rate of rise of COVID-19 cases
To investigate the relationship between Tmin and COVID-19 transmissibility, we devised a two-step process. In the first step, starting with the date of the first case for each country we examined linear regressions of the log of the daily cumulative COVID-19 cases versus time in sequential days (Fig 1). From this analysis we calculated the resulting regression coefficients (which we termed coefficients of time (CT) with units of log cases/day) for each country (Table 1). To determine consecutive days, we converted calendar dates to sequential serial numbers that can be used in calculations-e.g. 1/22/2020 was converted to 43852 because it is 43,851 days after 1/1/1900. The CT of log confirmed cases vs sequential days for all fifty countries yielded a mean overall R 2 of 0.919 +/-0.049 SD, range 0.792-0.986, and p-value mean 4.20 x10 -14 ± 2.61 x10 -13 , range 1.00 x10 -15 -1.84 x 10 −12 (Table 1).
In the second step, we investigated whether CT was associated with the median daily T min (Tmin_Med) in each country for the same time period that cases had been sampled. Because of the approximate proportionality of fitted CT with the standardized residuals when performing linear regression with Tmin_Med, we used a reciprocal transformation of the dependent variable [31] Thus, for the second regression for all countries we used the reciprocal of CT (CT -1 ) (Fig 2).
To account for susceptibility variation resulting from age differences that could affect both the comorbidity levels as well as social gathering behavior, we used each country's population median age (MA) as another covariate [32]. Because the 50 countries varied widely in size, and the density of the population could impact transmission, we used land area per capita (LAPC) as an additional covariate [33]. Finally, each country differed with respect to when it reported its identified cases, so that between the target dates of 1/22/20-4/6/20, there were different durations of time over which the course of transmission occurred. To account for this, the total number of days for each country during which they had cases of COVID-19 between 1/ 22/20 and 4/6/20 (Days of Cases) was used as a final covariate. In multivariable regression analysis of CT -1 versus the median of Tmin_Med (Table 2), we found that the rate of increase of the log of confirmed cases was significantly associated with the covariates Tmin_Med, LAPC, MA and Days of Cases (adjusted R 2 = 0.610, p = 1.45 x 10 −6 ; see  Table 2).

Sensitivity of the correlation between CT -1 and Tmin_Med
In our two-step model, first we aggregated the outcome measure CT or CT -1 at the country level and applied a standard linear regression; second, using multivariable regression analysis we associated CT or CT -1 with local meteorological and demographic independent variables. Since the initial aggregation might have played a role in determining response measures, we performed a sensitivity analysis by reducing sampling systematically (by approximately 10-20%) in the first stage [34]. This permitted us to test the robustness of the findings of our model.

Analysis of dewpoint and its influence on COVID-19 rate of rise of cases
In addition to regression of CT -1 with Tmin_Med, we also examined the median daily DP (DP_Med) for each country. We found that the rate of increase of the log of confirmed cases per unit time using covariates DP_Med, LAPC, MA and Days of Cases, when compared to  Table 1) for each country were converted to their reciprocal values (CT -1 ). This was done because of an approximate proportionality of fitted CT with the standardized residuals when performing linear regression with Tmin MEDIAN [31]. A second regression was then performed using the CT -1 for each country as the dependent variable and independent variables consisting of median minimal temperature [Tmin_MEDIAN] with or without regional demographic variables (i.e. land area and population median age). The total number of days for each country during which they had cases of COVID-19 between 1/22/20 and 4/6/20 (Days of Cases) was also used as a covariate. Shown are the resulting estimated regression coefficients and the R 2 for each regression. To control for multiple comparisons, since the same CT -1 and Tmin was used for Confirmed, Deaths and Recovered cases, P values were multiplied by 3. P values are displayed with progressively more stars for each 10-fold decrease in value according to this pattern: 0.05, 0.005, 0.0005, 0.00005, etc = � , �� , ��� , ���� , etc. NS stands for Not Significant (P>0.05). To control for multiple comparisons-since the same CT -1 and Tmin was used for Confirmed, Deaths and Recovered cases-p values were multiplied by 3 to account for the repeated measures of cases (See S1 and S2 Tables for further details). All countries were included in the Confirmed Case Analysis. For Death Cases there were 4 countries not included for the following three reasons: 1) Insufficient data for Kuwait (3 days of 0 log cases) and Vietnam (no days of deaths), 2) Inconsistent and limited data for Slovakia (4 days 0 log cases, 10 days no reported cases, then 5 days 0 log cases followed by a single 0.3 log cases), 3) Outlier: Iceland. Recovered Cases had one exclusion: Outlier: Slovakia. Outliers were determined using the ROUT method of Graph Pad's Prism software. Briefly, ROUT first fits a model to the data using a robust method where outliers have little impact.
Then it uses a new outlier detection method, based on the false discovery rate, to decide which points are far enough from the prediction of the model to be called outliers [26]. https://doi.org/10.1371/journal.pone.0246167.t002

PLOS ONE
Weather effects on SARS-CoV-2 transmission Tmin_Med, showed a comparably robust and statistically significant association (adjusted R 2 = 0.599, p = 2.68 x10 -6 ). Because of the collinear relationship between Tmin_Med and DP_Med, which in linear regression resulted in an adjusted R 2 of 0.76, (p = 1.4 x 10 −16 ), we assumed that, since Tmin_Med had the higher correlation with cases, it was likely the primary factor underlying the association between weather and COVID-19 transmission.

Magnitude of effect of temperature change on rate of transmission of COVID-19
From the second regression described above, the following equation was derived: CT ¼ 0:116 À 0:000911 x Tmin Median It was therefore possible to obtain CT for temperatures ranging from 30˚F to 100˚F-i.e. 0.0885 (log cases/day) to 0.0248 (log cases/day), respectively. This analysis indicated that an increase in temperature from 30 to 100˚F was associated with a 72% decrease in the rate of number of confirmed cases per day (CT 0.0885 to 0.0248), whereas a decrease in temperature from 100˚to 30˚F represents a 257% increase in the rate of change of confirmed cases (CT 0.0248 to 0.0885). This corresponds to 1.03% decrease in confirmed cases for each degree increase in temperature, and a 3.68% increase for each degree decrease in temperature (e.g. -72/(100-30) and 368/(100-30)).

Correlation with rates of death and recovered cases of COVID-19
The correlation between the reciprocal of the rise in the log of the number of deaths and recovered cases of COVID-19 (CT -1 ) vs Tmin_Med was 47% and 45% less than that of confirmed cases, respectively (see Table 2). Moreover, the covariates of LAPC and MA for death and recovered cases had regression p values of 0.5 or greater, suggesting that the impact of these covariates on confirmed cases was specific.

Relationship of CT to doubling time (Td)
As described in the method section, we explored the relationship between CT and doubling time (T d ), a standard measure of rate of viral transmission [29]. Doubling times are the intervals at which the cumulative incidence of cases doubles. Under conditions of exponential growth, with a constant growth rate r, the doubling time equals (ln 2)/r [30]. Thus, an increase in the doubling time represents a decreased rate of transmission. We calculated T d from our data obtained in our first regression, and we took the following steps: we used our first regression of log of confirmed cases vs time to interpolate on which days the number of cases doubled for each country utilizing a centered second order polynomial model to fit the data (mean R 2 = 0.958, Standard Deviation = +/-0.025). Because we were using log of confirmed cases we converted doubling of cases to log form (e.g. 2, 4, 8, 16,

Details of the model formulations (including equations in appropriate context)
We began by performing a two-step regression to derive the relationship between temperature and rate of increase in COVID-19 cases.
First regression involved the rate of rise of confirmed cases of COVID-19 in 50 countries, and the associated local meteorological and demographic variables: A second regression was then performed using the inverse of CT for each country as the dependent variable and the independent variables consisted of associated meteorological variables (i.e. median minimal temperature [Tmin_Med] or median dewpoint [DP]) with or without regional demographic variables (i.e. land area per capita and population median age):

Discussion
The major finding of this study is that the rates of transmission of SARS-CoV-2 infections are robustly associated with the ambient atmospheric temperature. The results of our analysis of longitudinal data from 50 representative countries in the Northern Hemisphere suggest that changes in the minimum daily temperature are associated with alterations in SARS CoV-2 transmissibility and viability. In addition, we found that the median daily dewpoint had a significant effect on the cumulative confirmed COVID-19 cases over time. However, because of the collinearity between dewpoint and atmospheric temperature over time, it was not possible to deduce whether the two atmospheric parameters independently affect SAR CoV-2 transmission. Since dewpoint is less predictable, and temperature more frequently monitored and used as a meteorological variable, we focused on changes in temperature for its potential predictive value.
We found that in 50 representative countries spread out over 70% of the possible latitude and longitude across the Northern Hemisphere, a slower rise in the log of daily cumulative number of confirmed cases was associated with higher daily minimum temperatures and dewpoints over the period from 1/22/20 to 4/6/20. This finding was robust, yielding an adjusted R 2 = 0.610, p = 1.45 x 10 −6 (using covariates Tmin_Med, median age, land area per capita, and days of cases). The correlation between the reciprocal of the rise in the log of the daily number of death and recovered cases (CT -1 ) vs Tmin_Med was almost 50% less than that of confirmed cases. This is not surprising, given that the impact of temperature on seasonal respiratory viruses has been found to have a consistent influence on the stability and transmission of infection, and not on the morbidity of infected individuals [2]. In view of these findings, atmospheric temperature would be expected to influence confirmed cases, and not deaths or the number of recovered subjects (e.g. once an infection occurred, temperature exerts a small or absent effect on the course of COVID-19 morbidity).
This study was able to achieve a clear and robust association between meteorological changes (i.e. temperature and dewpoint) and transmissibility of COVD-19, and to quantitate the magnitude of these correlations, because of the following steps taken: 1. we chose systematic criteria to select 50 representative countries that were selected to increase the probability that the data gathered (e.g. confirmed, deaths, recovered) and the rate of spread of COVID-19 would be comparable or at least randomly biased; 2. we took great care to ensure that we used weather stations located as close to the longitude and latitude corresponding to the cases in each country as presented by the Johns Hopkins real time COVID-19 tracking database; 3. our two-step regression model permitted us to relate the rate of rise of cases to changes in meteorological variables; 4. we focusing our analysis on the initial rise in cases across the various countries. During this early phase there was a delayed response on the part of the country (usually both the government and the populace) to institute and practice non-pharmaceutical interventions to slow the spread of the virus. We were therefore able to capture a period where weather and not social interventions played the predominant role in impacting transmission. And 5. after investigating numerous independent variables, we were able to select those that mediated the greatest association between meteorological variables and rate of transmission in our multivariable regression model.
Although our study is a systematic and quantitative analysis of the dependence of SAR-CoV-2 infection rates on temperature dewpoint that has demonstrated a robust and significant correlation (Table 2), it has several notable limitations. This is an ecological study and is therefore potentially subject to ecological fallacy [35]. Moreover, our data are not direct measures of individual-to-individual transmission, so cause and effect relationships cannot be established. However, the biological plausibility for our hypothesis is supported by previous work showing that the transmission rates of SARS-CoV-2 and similar viruses depend upon atmospheric temperature [2,36].
In temperate regions, the annual occurrence of respiratory viral diseases during the winter season-from the common cold to influenza-has been appreciated for several thousand years [2]. A similar seasonal pattern of infections has been reported for SARS-CoV, which was prevalent mostly during winter months [2]. Despite this, there is sparse evidence on the seasonal behavior of the novel SARS-CoV-2 [6], and there are conflicting reports on how its transmission is affected by meteorological conditions. In this context, our study provides the most comprehensive and up-to-date evidence for a robust and significant impact of temperature and dewpoint on SARS-CoV-2 transmissibility.
Further research is needed to clarify whether the association between temperature and SARS-CoV-2 transmissibility described here has a biological underpinning, and whether our model is accurate for seasonal predictions, surveillance and preparedness related to the spread of the virus over different geographical areas. Nonetheless, should the association between temperature and the rate of transmission hold true (either because of an underlying causal relationship or because of the empirical relationship described here), it may be possible to approximate the impact of weather on the rates of transmission of SARS-CoV-2. Because atmospheric parameters are measurable, and routinely tracked on a daily basis across the world, determining the contribution of temperature to the rate of transmission should make it possible to quantify the contribution of non-weather factors (e.g. social and non-pharmaceutical interventions (NPI))-since observed rate of transmission is equal to the contribution of meteorological variables plus the contribution of NPI. Once calibrated to the contributions of meteorological variables and social interventions and other NPI, such estimates may be able to inform and guide systematic academic, business and governmental decisions about when and how to impose or relax shelter-in-place or social distancing guidelines. Such evidence-based policies may succeed in minimizing social and economic disruptions due to the pandemic while at the same time optimizing public health and safety. Finally, our analysis predicts that, between the range of 30 to 100˚F, the decrease in the rate of COVID-19 transmission with increasing temperature (1% per degree F) is smaller than the increase rate of transmission with decreasing temperatures (3.7% per degree F) in the 50 representative countries in the Northern Hemisphere countries examined. If this model is correct, then it implies that effort invested in containing, minimizing, and ideally eliminating the spread of COVID-19 during Spring and Summer months could pay off significantly in the Fall and Winter, due to potentially disproportionate effects of decreasing temperatures compared to increasing temperatures on the current pandemic.
Supporting information S1 Table. The regression coefficients (abbreviated CT for coefficients of time) from the linear regression of the log of cumulative COVID-19 cases vs time (see Table 1