Skip to main content

Transmissibility of COVID-19 in 11 major cities in China and its association with temperature and humidity in Beijing, Shanghai, Guangzhou, and Chengdu

Abstract

Background

The new coronavirus disease COVID-19 began in December 2019 and has spread rapidly by human-to-human transmission. This study evaluated the transmissibility of the infectious disease and analyzed its association with temperature and humidity to study the propagation pattern of COVID-19.

Methods

In this study, we revised the reported data in Wuhan based on several assumptions to estimate the actual number of confirmed cases considering that perhaps not all cases could be detected and reported in the complex situation there. Then we used the equation derived from the Susceptible-Exposed-Infectious-Recovered (SEIR) model to calculate R0 from January 24, 2020 to February 13, 2020 in 11 major cities in China for comparison. With the calculation results, we conducted correlation analysis and regression analysis between R0 and temperature and humidity for four major cities in China to see the association between the transmissibility of COVID-19 and the weather variables.

Results

It was estimated that the cumulative number of confirmed cases had exceeded 45 000 by February 13, 2020 in Wuhan. The average R0 in Wuhan was 2.7, significantly higher than those in other cities ranging from 1.8 to 2.4. The inflection points in the cities outside Hubei Province were between January 30, 2020 and February 3, 2020, while there had not been an obvious downward trend of R0 in Wuhan. R0 negatively correlated with both temperature and humidity, which was significant at the 0.01 level.

Conclusions

The transmissibility of COVID-19 was strong and importance should be attached to the intervention of its transmission especially in Wuhan. According to the correlation between R0 and weather, the spread of disease will be suppressed as the weather warms.

Background

On December 8, 2019, the first case of unexplained pneumonia was officially reported in Wuhan, the capital of Hubei Province in China [1]. There have been reports of the new coronavirus disease (coronavirus disease 2019, COVID-19, named by the World Health Organization on February 11, 2020) since December 2019 [1, 2]. As was reported by the National Health Commission of the People’s Republic of China, the number of confirmed cases had reached 63 851 by February 13, 2020 in China, including 1380 deaths. On the same day, Hubei Province alone totally had 51 986 confirmed cases including 1318 deaths, accounting for 81.4% and 95.5% of the whole country respectively. Among them there were 35 991 confirmed cases and 1016 deaths in Wuhan, accounting for 69.2% and 77.1% of the number in Hubei Province respectively [3]. The cumulative number of confirmed cases keeps rising, indicating the strong transmissibility of COVID-19, especially in Wuhan, Hubei Province. Therefore, it is of great importance to adopt reasonable indicators to assess the transmission ability of the disease, based on which effective intervention and control measures could be put forward [4, 5].

The basic reproduction number (R0) refers to the expected number of cases generated from a single case when all people are susceptible to infection [6]. It is widely used to evaluate the transmission ability of an emerging infectious disease and determine what degree of control measures should be taken to eradicate the disease [710]. When R0>1, the disease starts to spread; and when R0<1, the disease is effectively controlled [11]. R0 is influenced by many other factors except for the characteristics of the disease itself, such as conditions of the environment, policies of the government, people’s awareness of infectious diseases, and social behavior. Therefore, we can use R0 to measure the transmissibility of COVID-19 and analyze its influencing factors, which provides data support for suggestion-proposing and decision-making.

Research on transmissible diseases like influenza [12], severe acute respiratory syndrome (SARS) [13] and Middle East respiratory syndrome (MERS) [14] has found that disease transmission is associated with temperature and humidity of the environment [1520]. In terms of biological methods, influenza virus spread was found to be promoted by cold temperature and low relative humidity with the guinea pig as a model host [21]; besides, an experiment on the SARS coronavirus indicated that high temperature and high humidity suppressed the spread of the virus [22]; similarly, MERS coronavirus was more stable when temperature or humidity was lower [23]. In terms of statistical methods, case studies of SARS in four major cities in China suggested that the transmissibility had a close relationship with temperature and its variation [24]; and a regression equation was derived to show how temperature, relative humidity, and wind velocity affected the transmission of SARS [25]. Thus we wonder if the spread of COVID-19 follows a similar pattern. Considering that R0 is useful for measuring the transmission ability of infectious diseases, we conducted association analyses between R0 and temperature, relative humidity, and absolute humidity respectively. Statistical methods such as correlation and regression were adopted for the analysis.

This paper measured the transmissibility of COVID-19 with R0 and analyzed its correlation with temperature and humidity. First, we revised the epidemiological data in Wuhan to make R0 more accurate. Second, we calculated R0 and compared the average value and developing trend of R0 in 11 cities including Wuhan. Third, we conducted correlation and regression analysis between R0 and temperature and humidity to see the association between R0 and weather.

Methods

Data acquisition and preprocessing

The daily accumulative number of confirmed cases and new additions are reported by the National Health Commission of the People’s Republic of China as well as the health commission of each province on the official website. An R package has been developed to access the epidemiological data directly [26]. The R package was used by us to acquire the number of total cases and new additions from January 18, 2020 to February 13, 2020 in Wuhan, Hubei Province considering that the situation there was complex and needed much attention. Besides, we also collected the daily-reported accumulative number of confirmed cases from January 24, 2020 to February 13, 2020 in 10 Chinese major cities outside Hubei Province including Beijing, Chengdu, Chongqing, Guangzhou, Hangzhou, Hefei, Nanjing, Shanghai, Shenzhen, and Zhengzhou (listed by initials) for further calculation, estimation, and analysis. The reasons for selecting those 10 cities were that they were first-tier cities or capital cities in China with the top number of cases. Certainly, Wuhan also met the criteria. Those cities could well represent the process status of the disease based on which disposal measures could be put forward.

As for Wuhan, it was estimated by Imperial College London, UK that the total number of confirmed diagnoses had reached 4000 by January 18, 2020 [27], which was much higher than the officially reported number. So we attempted to revise the data in Wuhan to infer the actual transmissibility of the new coronavirus. With the substantial enhancement of case detection and reporting, the differences between the official numbers and the estimates are predicted to be fewer and fewer. There are several assumptions for the data-preprocessing procedure:

  1. 1)

    The first case appeared on December 8, 2019 in Wuhan and transmission started from that day on [1, 28].

  2. 2)

    The cumulative number of cases Y(t) by day t since the first single case followed the exponential function Y(t)=eλt in early development [29].

  3. 3)

    The cumulative number of cases on January 18, 2020 was 4000, that was, Y(41)=4000 [27].

  4. 4)

    From February 13, 2020 on, all cases in Wuhan can be confirmed and the number of daily new cases is correct, given that the number of newly confirmed diagnoses on February 12, 2020 in Wuhan increased significantly, exceeding 10 000.

Based on those assumptions, the data-revising procedure in Wuhan is as follows:

  1. 1)

    According to assumption 2 and 3, the exponential growth rate is estimated as λ=ln[Y(41)]/41.

  2. 2)

    According to assumption 2 and 3, the number of new additions on January 18, 2020 equals Y(41)−Y(40)=4000−eλ40=733.

  3. 3)

    According to assumption 4, the number of new additions on February 13, 2020 is 2997, which is consistent with the officially reported number.

  4. 4)

    According to assumption 2, the daily number of new additions y(t) can be calculated by

    $$\begin{array}{@{}rcl@{}} y(t)&=&Y(t)-Y(t-1)\\ &=&e^{\lambda t}-e^{\lambda (t-1)}\\ &=&e^{\lambda t}\left(1-e^{-\lambda}\right). \end{array} $$
    (1)

    Thus

    $$\begin{array}{@{}rcl@{}} ln[y(t)]=ln\left(1-e^{-\lambda}\right)+\lambda t. \end{array} $$
    (2)

    So the relationship between ln[y(t)] and t is linear. Replace ln(1−eλ) with a and λ with b in Eq. (2), and the coefficients a and b of the linear equation can be determined by substituting y(41)=733 and y(67)=2997 into the equation respectively.

  5. 5)

    The number of new additions each day from January 19, 2020 to February 12, 2020 can be calculated through the equation y(t)=ea+bt, where a and b are the coefficients obtained in procedure 4.

  6. 6)

    With the daily number of new additions known, the daily cumulative number of cases from January 19, 2020 to February 13, 2020 can be calculated by Y(t)=Y(t−1)+y(t),t=42,43,...,67.

As for other cities outside Hubei Province, it is assumed that the officially reported data is accurate. Based on the relationship ln[Y(t)]=λt, we performed logarithmic fitting between the cumulative number of diagnoses and time and inferred that transmission started on December 27, 2019 outside Hubei Province.

Calculation of the basic reproduction number

The basic reproduction number indicates the average number of people infected by a patient during the infectious period in the absence of control interventions [6]. It is also denoted R0, which measures the transmissibility of infectious diseases. There are several ways to estimate R0, including formula derivation [30, 31] and model fitting [3234].

We describe the transmission pattern of COVID-19 with the Susceptible-Exposed-Infectious-Recovered (SEIR) model. In the exposed stage, an individual infection is not able to infect others. The duration of the exposed stage TE is also called the latent period. While in the infectious stage with a duration of TI, an infected person does infect susceptible people. Assuming that the cumulative number of confirmed diagnoses increases exponentially in the early stages of an epidemic, the relationship between the basic reproduction number R0 and the exponential growth rate λ can be written as [35].

$$\begin{array}{@{}rcl@{}} R_{0} = \left(1+\lambda T_{E}\right)\left(1+\lambda T_{I}\right). \end{array} $$
(3)

The serial interval Tg is the sum of TE and TI. Let f=TE/Tg be the ratio of the latent period to the serial interval, and then the basic reproduction number can be expressed as [29].

$$\begin{array}{@{}rcl@{}} R_{0} &=& 1+\lambda\left(T_{E}+T_{I}\right)+\lambda^{2}T_{E}T_{I}\\ &=& 1+\lambda T_{g}+\lambda^{2}T_{E}\left(T_{g}-T_{E}\right)\\ &=& 1+\lambda T_{g}+\lambda^{2}fT_{g}\left(T_{g}-fT_{g}\right)\\ &=& 1+\lambda T_{g}+f(1-f)\left(\lambda T_{g}\right)^{2}. \end{array} $$
(4)

The exponential growth rate is λ=ln[Y(t)]/t, where t is the number of days required to generate the cumulative number of Y(t) cases from the first case. According to the research on the first 425 patients with confirmed COVID-19, the mean latent period TE = 5.2 (days) and the mean serial interval Tg = 7.5 (days) [36]. Adopting these values, we can calculate the ratio of the latent period to the serial interval by f=TE/Tg=5.2/7.5=0.69.

Correlation and regression analysis between R0 and weather

Correlation analysis is a commonly used statistical method to study the relationship between variables [37]. Regression analysis determines the quantitative relationship between two variables in statistics [38]. Among all kinds of regression methods, linear regression establishes the relationship between the dependent variable Y and the independent variable X with a linear equation Y=a+bX [39]. There are two coefficients in the equation, a as the intercept and b as the slope. We performed correlation analysis and linear regression between R0 and weather variables with the statistical analysis software IBM SPSS Statistics 25. The procedure is listed below.

  1. 1)

    We collected the data of the daily average temperature and relative humidity from January 24, 2020 to February 13, 2020 in four Chinese major cities which were Beijing (the capital of China), Shanghai (the municipality of China), Guangzhou (the capital of Guangdong Province) and Chengdu (the capital of Sichuan Province). We calculated absolute humidity from the temperature and relative humidity.

  2. 2)

    We imported the data of temperature, relative humidity, and absolute humidity together with R0 into the SPSS software and added cities as the classification label.

  3. 3)

    Through correlation analysis, the Pearson correlation coefficients between R0 and temperature, relative humidity, and absolute humidity were calculated respectively.

  4. 4)

    Through regression analysis, the intercept a and the slope b of the linear equation were estimated with R0 as the dependent variable Y and temperature, relative humidity or absolute humidity as the independent variable X.

  5. 5)

    We split the data by the city label and repeated procedure 3 and 4 for each city separately.

Sensitivity analysis of R0

To analyze the sensitivity of R0 to the three key parameters in Eq. (4): R0=1+λTg+f(1−f)(λTg)2, we differentiated R0 to λ, Tg and f respectively:

$$\begin{array}{@{}rcl@{}} \frac{\partial R_{0}}{\partial \lambda} &=& T_{g} + 2f\left(1-f\right)\lambda T_{g}^{2}, \end{array} $$
(5)
$$\begin{array}{@{}rcl@{}} \frac{\partial R_{0}}{\partial T_{g}}&=& \lambda + 2f\left(1-f\right)\lambda^{2} T_{g},\\ \frac{\partial R_{0}}{\partial f}&=&\left(1-2f\right)\left(\lambda T_{g}\right)^{2}. \end{array} $$
(6) (7)

The sensitivity of the basic reproduction number R0 to the exponential growth rate λ, the serial interval Tg, and the latent period ratio f can be estimated according to the range of variables and the scale of partial derivatives.

Results

Comparisons of transmission among different cities

The comparison between officially reported data and revised data in Wuhan is presented in Fig. 1 with important points marked on it. The estimated number of cumulative cases was higher than the official number every day, and it had reached 46 933 by February 13, 2020, which was 1.3 times that of the official number 35 991. The unusual high peak of new cases on February 12, 2020 was smoothed by revision.

Fig. 1
figure 1

Comparisons between reported and revised data in Wuhan

The calculation results of the basic reproduction number R0 from January 24, 2020 to February 13, 2020 in 11 Chinese major cities are shown in Fig. 2. The values with the label “Wuhan” were calculated using the officially reported number of cases, while those with “Wuhan (revised)” were calculated using the revised number of cases. In this way, the broken line of “Wuhan” reflects the changing trend of R0, and the one of “Wuhan (revised)” reflects the value size of R0. It is assumed that the cumulative number of confirmed cases reported officially in cities outside Hubei Province is accurate, so the broken lines of the other 10 cities represent not only trends but also actual values.

Fig. 2
figure 2

Calculation results of the basic reproduction number

As can be seen from Fig. 2, R0 in Wuhan is significantly higher than those in cities outside Hubei Province. Besides, R0 in cities outside Hubei Province has begun to decrease, while R0 in Wuhan does not show a significant downward trend.

For a more detailed analysis, the average basic reproduction number of the 21 days in each city and the date of the inflection point are presented in Table 1. The cities are listed by the average R0 from high to low. The inflection point refers to the day after which R0 shows a downward trend.

Table 1 The average R0 and the inflection point of each city (listed by the average R0)

It can be seen from Table 1 that the average R0 in Wuhan far exceeds those in other cities, which is 0.3 higher than that in Chongqing, the city which ranks second. It should be noted that the average R0 in Wuhan is calculated with the revised data to better fit the real value. In fact, the average basic reproduction number calculated with the officially reported data is also much higher than those in other cities, which is 2.4.

The inflection points of cities outside Hubei Province range from January 30 to February 3, while the inflection point of Wuhan had not appeared because the number of confirmed cases had kept increasing rapidly by February 13, 2020. Although R0 in Wuhan reaches a peak on February 12, it cannot be determined that February 12 is the inflection point. Because since that day, Hubei Province has included the number of clinically diagnosed cases into the number of confirmed cases. The modification of the diagnostic criteria leads to a sudden increase of newly confirmed patients, which explains why R0 is particularly high on February 12.

Correlation between R0 and temperature

The Pearson correlation coefficients and significance between R0 and temperature are shown in Table 2. The row of “Summary” suggests that calculated as a whole, the correlation between R0 and temperature is statistically significant at the 0.01 level. The correlation coefficient is -0.459, so R0 and temperature have a negative correlation, which means that R0 decreases as the temperature increases. The higher the temperature, the lower the transmission capability. As for the analysis of each city, R0 negatively correlates with temperature in Shanghai and Chengdu, correlation significant at the 0.01 level. Correlation is not significant in Beijing and Guangzhou. Over the study period, the average R0 in Beijing, Shanghai, Guangzhou, and Chengdu are 2.3, 2.2, 2.2, and 2.0 respectively and the average temperatures are -1.0 C, 7.9 C, 14.9 C, and 9.9 C respectively. There is not a significant relationship between the average R0 in a city versus its average temperature (r=−0.486, P>0.5).

Table 2 Correlation analysis between R0 and temperature

Linear regression was performed on the data for all cities combined as well as the data in Shanghai and Chengdu which showed a significant correlation. Table 3 presents the linear regression results. Replace a and b in the equation R0=a+bT (where T is temperature) with the corresponding actual values in Table 3, and correlation between R0 and temperature can be expressed more precisely. For example, the linear regression equation of Shanghai is written as R0=2.424−0.026T. It can be inferred from b<0 that R0 negatively correlates with temperature in Shanghai, which is consistent with the correlation analysis result above.

Table 3 Linear regression analysis of temperature to R0

We plotted every pair of temperature and R0 in a city or the whole data on the scatter figure to make correlation more intuitive, which was presented in Fig. 3. The regression lines followed the corresponding linear regression equations.

Fig. 3
figure 3

Scatter plot of temperature and basic reproduction number

Correlation between R0 and relative humidity

The Pearson correlation coefficients and significance between R0 and relative humidity are presented in Table 4. According to the first row, the correlation between R0 and relative humidity is statistically significant at the 0.01 level in general. The correlation coefficient is -0.391, indicating that R0 decreases as the relative humidity increases. As for the analysis of each city, R0 negatively correlates with relative humidity in Beijing and Shanghai, which is significant at the 0.01 level. While the correlation is significantly positive in Chengdu at the 0.01 level, which implies that the transmission ability and relative humidity have consistent trends there. Correlation is not significant in Guangzhou.

Table 4 Correlation analysis between R0 and relative humidity

The correlation was significant in Beijing, Shanghai, and Chengdu, and thus we conducted linear regression on the data of the three cities as well as the summary of all cities. The linear regression results are presented in Table 5. Replace a and b in the equation R0=a+bRH (where RH is relative humidity) with the corresponding actual values in Table 5, and the correlation between R0 and relative humidity can be expressed with a quantitative method.

Table 5 Linear regression analysis of relative humidity to R0

The scatterplots and corresponding regression lines of relative humidity and R0 summarized across all cities and by individual cities are presented in Fig. 4.

Correlation between R0 and absolute humidity

The Pearson correlation coefficients and significance between R0 and absolute humidity are presented in Table 6. The negative correlation between R0 and absolute humidity is significant in general as well as in Beijing, Shanghai and Guangzhou and the absolute values of the Pearson correlation coefficients for absolute humidity are larger than those for relative humidity, indicating that the relationship is stronger for absolute humidity than relative humidity. The correlation is not significant in Chengdu.

Fig. 4
figure 4

Scatter plot of relative humidity and basic reproduction number

Table 6 Correlation analysis between R0 and absolute humidity

We conducted linear regression on the data of Beijing, Shanghai, Guangzhou as well as the summary of all cities. The linear regression results are presented in Table 7. Replace a and b in the equation R0=a+bAH (where AH is absolute humidity) with the corresponding actual values in Table 7, and the correlation between R0 and absolute humidity can be expressed with a quantitative method.

Table 7 Linear regression analysis of absolute humidity to R0

The scatterplots and corresponding regression lines of absolute humidity and R0 summarized across all cities and by individual cities are presented in Fig. 5.

Fig. 5
figure 5

Scatter plot of absolute humidity and basic reproduction number

Sensitivity of R0 to parameters

Substitute the variables in Eqs. (47) with λ=0.1372 (which is the average λ from January 24 to February 13 in Beijing), Tg=7.5 and f=0.69, and the specific values can be calculated:

$$\begin{array}{@{}rcl@{}} R_{0}&=&2.3, \end{array} $$
(8)
$$\begin{array}{@{}rcl@{}} \frac{\partial R_{0}}{\partial \lambda} &=&10.8,\\ \frac{\partial R_{0}}{\partial T_{g}}&=&0.197,\\ \frac{\partial R_{0}}{\partial f}&=&-0.41. \end{array} $$
(9) (10) (11)

When the variables fluctuate within a small range around the given value, R0 increases as λ or Tg increases and decreases as f increases. λ, Tg and f range at the scales of 10−2, 100 and 10−1 respectively. And the scales of their partial derivatives are 101, 10−1 and 10−1. Thus the fluctuation scales of R0 are 10−1, 10−1 and 10−2 corresponding to λ, Tg and f, which implies that R0 is more sensitive to λ and Tg than f. The accuracy of parameters or variables is important for the estimation of the basic reproduction number. As the research on COVID-19 progresses, we can get more precise data and better describe the transmission pattern of the new coronavirus. But the calculation in this paper still makes sense, considering that we focus on relative values instead of absolute values of R0 in comparison and correlation analysis. Results are reasonable as long as we use the consistent equation and parameters to calculate R0. By comparison, we can see that the control of COVID-19 is especially urgent in Wuhan and people in other cities should also attach importance to inhibiting the spread of the disease. The vigilance cannot languish until R0 drops below 1.

Discussion

Differences between correlation and causation

In this paper, we discovered the negative correlation between the transmissibility of COVID-19 and temperature and humidity. However, it should be emphasized that correlation is different from causation. According to the Oxford Dictionary, correlation is a connection between two things in which one thing changes as the other does, while causation is the process of one event causing or producing another event. We are not able to infer the causal relationship between two variables solely based on the correlation between them. Correlation is the necessary and insufficient condition of causation. Our results indicated that the transmissibility of COVID-19 was likely to decrease as the temperature and humidity increased. But it did not mean that the increase of temperature or humidity was the cause of the decrease of the transmissibility. We were not able to control other variables in the observation, such as population migration and interventions, which might also affect the transmissibility of COVID-19. So perhaps future work is needed to find out if the changes in temperature or humidity cause the changes in the transmissibility. For example, biological experiments can be conducted by setting the temperature or humidity as the independent variable and the transmissibility of the coronavirus as the dependent variable and controlling other irrelevant variables with the elimination method, constant method, matching method or randomization. Nevertheless, this paper makes sense in terms of confirming that the transmissibility of COVID-19 has a correlation with temperature and humidity and that there is probably a causation relationship between them which deserves further research.

Effects of temperature and humidity on the transmission of COVID-19

A recent study indicated that temperature and relative humidity held no significant associations with the transmissibility of COVID-19 [40]. It is a very comprehensive and well-conducted research, but we took a step further to take the time series into account by using everyday temperature and humidity. The results show that the overall correlation between R0 and temperature or humidity is significantly negative, which is consistent with the results of the biological and statistical research on other infectious diseases. It could be explained in several aspects. First, in terms of biological characteristics, a lot of research has confirmed that viruses decay more quickly at high temperature and high humidity [19, 41, 42]. Second, in terms of the transmission media, viruses spread as droplets or aerosols, which maintain large particle sizes at high humidity and thus can settle rapidly or be blocked by masks, nasal cavity, etc [19]. Third, in terms of human immunity, high temperature and high humidity protect the immune organs and benefit people’s health. To sum up, the spread of COVID-19 is likely to weaken at relatively high temperature and humidity and special attention should be paid to the prevention and control of COVID-19 in the coming winter.

As for the correlation in each city, R0 negatively correlates with both temperature and humidity in Shanghai; R0 negatively correlates with humidity in Beijing, while the correlation with temperature is not significant; R0 negatively correlates with absolute humidity in Guangzhou, while the correlation with temperature and relative humidity is not significant; R0 negatively correlates with temperature in Chengdu, while the correlation with relative humidity is positive and the correlation with absolute humidity is not significant. The deviation of the results may be due to several factors.

First, considering that COVID-19 began in winter, people’s activity and virus transmission mainly occur indoors. In China, the cities north of the Qinling Mountains-Huaihe River Line have central heating indoors in winter. Beijing is north of the Qinling Mountains-Huaihe River Line and Shanghai, Guangzhou, and Chengdu are south of the line. Therefore, the indoor temperature is probably much higher than the outdoor temperature in Beijing, while the indoor temperature may follow a similar pattern as the outdoor temperature in the other three cities. The indoor temperature is probably higher in Beijing than that in Shanghai, Guangzhou, and Chengdu. Although the indoor temperature and the outdoor temperature may have some association, it would be better if we could measure the indoor temperature directly. As for humidity, it has been found that outdoor absolute humidity can be more reliably used as a proxy for indoor exposure compared with relative humidity [43, 44]. Therefore, the correlation between R0 and absolute humidity may better reveal the situation indoors than relative humidity. Actually, the Pearson correlation coefficients between R0 and absolute humidity are larger than those between R0 and relative humidity, proving that the relationship is stronger for absolute humidity than relative humidity.

Second, although Beijing, Shanghai, Guangzhou, and Chengdu are all first-tier cities in China with many similarities like buildings, there are some differences between Chengdu and the other cities that may help explain the positive correlation with relative humidity. Chengdu is located in the southwest of China, the west of Sichuan Basin and the hinterland of Chengdu Plain with a subtropical monsoon humid climate, different from Beijing which has a warm temperate semi-humid continental monsoon climate. The air is more humid in Chengdu than that in Beijing. The climate in Chengdu is similar to the subtropical monsoon climate in Shanghai and Guangzhou, but Chengdu is an inland city while Shanghai and Guangzhou are coastal cities.

Third, the effect of weather on COVID-19 is complicated. The joint distribution between weather and potential confounders should be taken into account. For example, population movement might trigger the transmission of COVID-19 [45]. As for the effects of interventions, we have plotted the time series of temperatures from January 24, 2020 to February 13, 2020 in Beijing, Shanghai, Guangzhou, and Chengdu in Additional file 1: Figure S1. It could be seen from the figure that the temperature kept fluctuating during this period. Considering that the strength of interventions was relatively steady without big fluctuations, which was different from the trends of temperature, perhaps the effects of interventions could be separated from the trends in temperature.

Conclusions

In this paper, we calculated and compared the basic reproduction number of COVID-19 in 11 major cities in China and analyzed its association with temperature and humidity in Beijing, Shanghai, Guangzhou, and Chengdu to find out the transmissibility of COVID-19 in different cities and its changing trend with the weather. We conclude that the spread of COVID-19 is most violent in Wuhan, Hubei Province and R0 negatively correlates with temperature, relative humidity, and absolute humidity. Therefore, effective action should be taken to control the transmission of COVID-19 especially in Hubei Province and the transmissibility is predicted to be reduced as the weather warms.

Availability of data and materials

Not applicable.

Abbreviations

COVID-19:

Corona virus disease 2019

R 0 :

Basic reproduction number

SARS:

Severe acute respiratory syndrome

MERS:

Middle east respiratory syndrome

SEIR:

Susceptible-exposed-infectious-recovered

References

  1. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. N Engl J Med. 2020; 382(13):1199–207. https://doi.org/10.1056/NEJMoa2001316.

    Article  CAS  Google Scholar 

  2. Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. Lancet. 2020; 395(10225):689–97. https://doi.org/10.1016/S0140-67362030260-9.

    Article  CAS  Google Scholar 

  3. NHC. National Health Commission of the People’s Republic of China-Epidemic notification. 2020. http://www.nhc.gov.cn/xcs/yqtb/list_gzbd.shtml. Accessed 14 Feb 2020.

  4. Wiwanitkit V, Shi B, Xia S, Yang GJ, Zhou XN, Liu J. Research priorities in modeling the transmission risks of H7N9 bird flu. Infect Dis Poverty. 2013; 2(1):17. https://doi.org/10.1186/2049-9957-2-17.

    Article  Google Scholar 

  5. Lei L, Richards JS, Li ZH, Gong YF, Zhang SZ, Xiao N. A framework for assessing local transmission risk of imported malaria cases. Infect Dis Poverty. 2019; 8(1):43. https://doi.org/10.1186/s40249-019-0552-6.

    Article  Google Scholar 

  6. Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, Griffin J, Baggaley RF, Jenkins HE, Lyons EJ, et al. Pandemic potential of a strain of influenza A (H1N1): early findings. Science. 2009; 324(5934):1557–61. https://doi.org/10.1126/science.1176062.

    Article  CAS  Google Scholar 

  7. Fine P, Eames K, Heymann DL. “Herd immunity”: a rough guide. Clin Infect Dis. 2011; 52(7):911–6. https://doi.org/10.1093/cid/cir007.

    Article  Google Scholar 

  8. Cao Z, Zhang Q, Lu X, Pfeiffer D, Jia Z, Song H, et al. Estimating the effective reproduction number of the 2019-nCoV in China. medRxiv. 2020. https://doi.org/10.1101/2020.01.27.20018952.

  9. Liu T, Hu J, Xiao J, He G, Kang M, Rong Z, et al. Time-varying transmission dynamics of Novel Coronavirus Pneumonia in China. bioRxiv. 2020. https://doi.org/10.1101/2020.01.25.919787.

  10. Kretzschmar M, Teunis PFM, Pebody RG. Incidence and reproduction numbers of pertussis: estimates from serological and social contact data in five European countries. PLoS Med. 2010; 7(6):e1000291. https://doi.org/10.1371/journal.pmed.1000291.

    Article  Google Scholar 

  11. Anderson RM, Anderson B, May RM. Infectious diseases of humans: dynamics and control. Oxford: Oxford university press; 1992.

    Google Scholar 

  12. Mills CE, Robins JM, Lipsitch M. Transmissibility of 1918 pandemic influenza. Nature. 2004; 432(7019):904–6. https://doi.org/10.1038/nature03063.

    Article  CAS  Google Scholar 

  13. Ksiazek TG, Erdman D, Goldsmith CS, Zaki SR, Peret T, Emery S, Tong S, Urbani C, Comer JA, Lim W, et al. A novel coronavirus associated with severe acute respiratory syndrome. N Engl J Med. 2003; 348(20):1953–66. https://doi.org/10.1056/NEJMoa030781.

    Article  CAS  Google Scholar 

  14. Kucharski AJ, Althaus CL. The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission. Euro Surveill. 2015; 20(25):14–8. https://doi.org/10.2807/1560-7917.es2015.20.25.21167.

    Article  CAS  Google Scholar 

  15. Hemmes JH, Winkler KC, Kool SM. Virus survival as a seasonal factor in influenza and polimyelitis. Nature. 1960; 188(4748):430–1. https://doi.org/10.1038/188430a0.

    Article  CAS  Google Scholar 

  16. Shama J, Kohn M. Absolute humidity modulates influenza survival, transmission, and seasonality. Proc Natl Acad Sci U S A. 2009; 106(9):3243–8. https://doi.org/10.1073/pnas.0806852106.

    Article  Google Scholar 

  17. Steel J, Palese P, Lowen AC. Transmission of a 2009 pandemic influenza virus shows a sensitivity to temperature and humidity similar to that of an H3N2 seasonal strain. J Virol. 2011; 85(3):1400–2. https://doi.org/10.1128/JVI.02186-10.

    Article  CAS  Google Scholar 

  18. Dalziel BD, Kissler S, Gog JR, Viboud C, Bjørnstad ON, Metcalf CJE, et al. Urbanization and humidity shape the intensity of influenza epidemics in US cities. Science. 2018; 362(6410):75–9. https://doi.org/10.1126/science.aat6030.

    Article  CAS  Google Scholar 

  19. Marr LC, Tang JW, Van Mullekom J, Lakdawala SS. Mechanistic insights into the effect of humidity on airborne influenza virus survival, transmission and incidence. J R Soc Interface. 2019; 16(150):20180298. https://doi.org/10.1098/rsif.2018.0298.

    Article  CAS  Google Scholar 

  20. Ehelepola NDB, Ariyaratne K, Buddhadasa WMNP, Ratnayake S, Wickramasinghe M. A study of the correlation between dengue and weather in Kandy City, Sri Lanka (2003-2012) and lessons learned. Infect Dis Poverty. 2015; 4(1):42. https://doi.org/10.1186/s40249-015-0075-8.

    Article  CAS  Google Scholar 

  21. Lowen AC, Mubareka S, Steel J, Palese P. Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathog. 2007; 3(10):1470–6. https://doi.org/10.1371/journal.ppat.0030151.

    Article  CAS  Google Scholar 

  22. Chan KH, Peiris JSM, Lam SY, Poon LLM, Yuen KY, Seto WH. The effects of temperature and relative humidity on the viability of the SARS coronavirus. Adv Virol. 2011; 2011(1–7). https://doi.org/10.1155/2011/734690.

  23. van Doremalen N, Bushmaker T, Munster VJ. Stability of Middle East respiratory syndrome coronavirus (MERS-CoV) under different environmental conditions. Euro Surveill. 2013; 18(38):20590. https://doi.org/10.2807/1560-7917.es2013.18.38.20590.

    Article  Google Scholar 

  24. Tan J, Mu L, Huang J, Yu S, Chen B, Yin J. An initial investigation of the association between the SARS outbreak and weather: with the view of the environmental temperature and its variation. J Epidemiol Community Health. 2005; 59(3):186–92. https://doi.org/10.1136/jech.2004.020180.

    Article  Google Scholar 

  25. Yuan J, Yun H, Lan W, Wang W, Sullivan SG, Jia S, Bittles AH. A climatologic investigation of the SARS-CoV outbreak in Beijing, China. Am J Infect Control. 2006; 34(4):234–6. https://doi.org/10.1016/j.ajic.2005.12.006.

    Article  Google Scholar 

  26. Wu T, Ge X, Yu G. An R package and a website with real-time data on the COVID-19 coronavirus outbreak. medRxiv. 2020. https://doi.org/10.1101/2020.02.25.20027433.

  27. Imai N, Dorigatti I, Cori A, Donnelly C, Riley S, Ferguson NM. Report 2 - Estimating the potential total number of novel Coronavirus (2019-nCoV) cases in Wuhan City, China. Imperial College London. 2020. https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-2-update-case-estimates-covid-19/. Accessed 14 Feb 2020.

  28. Zhou T, Liu Q, Yang Z, Liao J, Yang K, Bai W, Lu X, Zhang W. Preliminary prediction of the basic reproduction number of the Wuhan novel coronavirus 2019-nCoV. J Evid Based Med. 2020; 13(1):3–7. https://doi.org/10.1111/jebm.12376.

    Article  Google Scholar 

  29. Lipsitch M, Cohen T, Cooper B, Robins JM, Ma S, James L, et al. Transmission dynamics and control of severe acute respiratory syndrome. Science. 2003; 300(5627):1966–70. https://doi.org/10.1126/science.1086616.

    Article  CAS  Google Scholar 

  30. Khan A, Naveed M, Dur-E-Ahmad M, Imran M. Estimating the basic reproductive ratio for the Ebola outbreak in Liberia and Sierra Leone. Infect Dis Poverty. 2015; 4(1):13. https://doi.org/10.1186/s40249-015-0043-3.

    Article  CAS  Google Scholar 

  31. Chowell G, Castillo-Chavez C, Fenimore PW, Kribs-Zaleta CM, Arriola L, Hyman JM. Model parameters and outbreak control for SARS. Emerg Infect Dis. 2004; 10(7):1258–63. https://doi.org/10.3201/eid1007.030647.

    Article  Google Scholar 

  32. Khan A, Hassan M, Imran M. Estimating the basic reproduction number for single-strain dengue fever epidemics. Infect Dis Poverty. 2014; 3(1):12. https://doi.org/10.1186/2049-9957-3-12.

    Article  Google Scholar 

  33. Althaus CL. Estimating the reproduction number of Ebola virus (EBOV) during the 2014 outbreak in West Africa. PLoS Curr. 2014; 6:1–5. https://doi.org/10.1371/currents.outbreaks.91afb5e0f279e7f29e7056095255b288.

    Google Scholar 

  34. Chen TM, Rui J, Wang QP, Zhao ZY, Cui JA, Yin L. A mathematical model for simulating the phase-based transmissibility of a novel coronavirus. Infect Dis Poverty. 2020; 9(1):24. https://doi.org/10.1186/s40249-020-00640-3.

    Article  CAS  Google Scholar 

  35. Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc Biol Sci. 2007; 274(1609):599–604. https://doi.org/10.1098/rspb.2006.3754.

    Article  CAS  Google Scholar 

  36. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. N Engl J Med. 2020; 382(13):1199–207. https://doi.org/10.1056/NEJMoa2001316.

    Article  CAS  Google Scholar 

  37. Zhang T, Shen S, Cheng C, Song C, Ye S. Long-range correlation analysis of soil temperature and moisture on A’rou hillsides, Babao River Basin. J Geophys Res-Atmos. 2018; 123(22):12606–20. https://doi.org/10.1029/2018JD029094.

    Google Scholar 

  38. Jike M, Itani O, Watanabe N, Buysse DJ, Kaneita Y. Long sleep duration and health outcomes: A systematic review, meta-analysis and meta-regression. Sleep Med Rev. 2018; 39:25–36. https://doi.org/10.1016/j.smrv.2017.06.011.

    Article  Google Scholar 

  39. Dumuid D, Stanford TE, Martin-Fernandez JA, Pedisic Z, Maher CA, Lewis LK, et al. Compositional data analysis for physical activity, sedentary time and sleep research. Stat Methods Med Res. 2018; 27(12):3726–38. https://doi.org/10.1177/0962280217710835.

    Article  Google Scholar 

  40. Yao Y, Pan J, Liu Z, Meng X, Wang W, Kan H, Wang W. No association of COVID-19 transmission with temperature or UV radiation in Chinese cities. Eur Respir J. 2020; 55(5):2000517. https://doi.org/10.1183/1399300300517-2020.

    Article  Google Scholar 

  41. Yang W, Marr LC. Mechanisms by which ambient humidity may affect viruses in aerosols. Appl Environ Microbiol. 2012; 78(19):6781–8. https://doi.org/10.1128/AEM.01658-12.

    Article  CAS  Google Scholar 

  42. Wolkoff P. Indoor air humidity, air quality, and health - An overview. Int J Hyg Environ Health. 2018; 221(3):376–90. https://doi.org/10.1016/j.ijheh.2018.01.015.

    Article  Google Scholar 

  43. Nguyen JL, Schwartz J, Dockery DW. The relationship between indoor and outdoor temperature, apparent temperature, relative humidity, and absolute humidity. Indoor Air. 2014; 24(1):103–12. https://doi.org/10.1111/ina.12052.

    Article  CAS  Google Scholar 

  44. Nguyen JL, Dockery DW. Daily indoor-to-outdoor temperature and humidity relationships: a sample across seasons and diverse climatic regions. Int J Biometeorol. 2016; 60(2):221–9. https://doi.org/10.1007/s00484-015-1019-5.

    Article  Google Scholar 

  45. Liu K, Ai S, Song S, Zhu G, Tian F, Li H, et al. Population movement, city closure in Wuhan and geographical expansion of the 2019-nCoV pneumonia infection in China in January 2020. Clin Infect Dis. 2020. https://doi.org/10.1093/cid/ciaa422.

Download references

Acknowledgements

Not applicable.

Funding

This work is supported by National Key R&D Program of China (No. 2017YFC0803300), National Science Foundation of China (Grant Nos. 7204100828, 91646201, U1633203), High-tech Discipline Construction Fundings for Universities in Beijing (Safety Science and Engineering) and Beijing Key Laboratory of City Integrated Emergency Response Science.

Author information

Authors and Affiliations

Authors

Contributions

HZ conceived the study. XG and YZ collected the data. XG analyzed the data and drafted the manuscript. XG and HZ interpreted the results. XG, HZ and YZ revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hui Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1

Figure S1. The time series of temperature in Beijing, Shanghai, Guangzhou and Chengdu.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guo, XJ., Zhang, H. & Zeng, YP. Transmissibility of COVID-19 in 11 major cities in China and its association with temperature and humidity in Beijing, Shanghai, Guangzhou, and Chengdu. Infect Dis Poverty 9, 87 (2020). https://doi.org/10.1186/s40249-020-00708-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40249-020-00708-0

Keywords