An investigation of transmission control measures during the first 50 days of the COVID-19 epidemic in China

Abstract

Responding to an outbreak of a novel coronavirus (agent of COVID-19) in December 2019, China banned travel to and from Wuhan city on 23 January and implemented a national emergency response.

We investigated the spread and control of COVID-19 using a unique data set including case reports, human movement and public health interventions.

The Wuhan shutdown was associated with the delayed arrival of COVID-19 in other cities by 2.91 days (95%CI: 2.54-3.29). Cities that implemented control measures pre-emptively reported fewer cases, on average, in the first week of their outbreaks (13.0; 7.1-18.8) compared with cities that started control later (20.6; 14.5-26.8).

Suspending intra-city public transport, closing entertainment venues and banning public gatherings were associated with reductions in case incidence.

The national emergency response appears to have delayed the growth and limited the size of the COVID-19 epidemic in China, averting hundreds of thousands of cases by 19 February (day 50).

On 31 December 2019, less than a month before the 2020 Spring Festival holiday, including the Chinese Lunar New Year, a cluster of pneumonia cases caused by an unknown pathogen was reported in Wuhan, a city of 11 million inhabitants and the largest transport hub in Central China. A novel coronavirus (1,2) was identified as the etiological agent (3,4) and human-to-human transmission of the virus (SARS-CoV-2) has been since confirmed (5,6). Further spatial spread of this disease was of great concern in view of the upcoming Spring Festival (“chunyun”) during which there are typically three billion travel movements over the 40-day holiday period, which runs from 15 days before the Spring Festival (Chinese Lunar New Year) to 25 days afterwards (7,8).

As there is currently neither a vaccine nor a specific drug treatment for COVID-19, a range of public health (non-pharmaceutical) interventions has been used to control the epidemic. In an attempt to prevent further dispersal of COVID-19 from its source, all transport was prohibited in and out of Wuhan city from 10:00h on 23 January 2020, followed by the whole of Hubei Province a day later. In terms of the population covered, this appears to be the largest attempted cordon sanitaire in human history.

On 23 January, China also raised its national public health response to the highest state of emergency—Level 1 of 4 levels of severity in the Chinese Emergency System, defined as an “extremely serious incident” (9). As part of the national emergency response, and in addition to the Wuhan city travel ban, suspected and confirmed cases have been isolated, public transport by bus and subway rail suspended, schools and entertainment venues have been closed, public gatherings banned, health checks carried out on migrants (“floating population”), travel prohibited in and out of cities, and information widely disseminated. Despite all these measures, COVID-19 remains a danger in China. Control measures taken in China potentially hold lessos for other countries around the world.

Although the spatial spread of infectious diseases has been intensively studied (10–15), including explicit studies of the role of human movement (16,17), the effectiveness of travel restrictions and social distancing measures in preventing the spread of infection is uncertain. For COVID-19, coronavirus transmission patterns and the impact of interventions are still poorly understood (6,7). We therefore carried out a quantitative analysis to investigate the role of travel restrictions and transmission control measures during the first 50 days of the COVID-19 epidemic in China, from 31 December 2019 to 19 February 2020 (Fig. 1). This period encompassed the 40 days of the Spring Festival holiday, 15 days before the Chinese Lunar New Year on 25 January and 25 days afterwards. The analysis is based on a unique geocoded repository of data on COVID-19 epidemiology, human movement, and public health (non-pharmaceutical) interventions. These data include the numbers of COVID-19 cases reported each day in each city of China, information on 4.3 million human movements from Wuhan city, and data on the timing and type of transmission control measures implemented across cities of China.

We first investigated the role of the Wuhan city travel ban, comparing travel in 2020 with that in previous years and exploring how holiday travel links to the dispersal of infection across China. During Spring Festival travel in 2017 and 2018, there was an average outflow of 5.2 million people from Wuhan city during the 15 days before the Chinese Lunar New Year. In 2020, this travel was interrupted by the Wuhan city shutdown, but 4.3 million people travelled out of the city between 11 January and the implementation of the ban on 23 January (Fig. 2A) (7). In 2017 and 2018, travel out of the city during the 25 days after the Chinese Lunar New Year averaged 6.7 million people each year. In 2020, the travel ban prevented almost all of that movement and markedly reduced the number of exportations of COVID-19 from Wuhan (7,8).

The dispersal of COVID-19 from Wuhan was rapid (Fig. 3A). A total of 262 cities reported cases within 28 days. For comparison, the 2009 influenza H1N1 pandemic took 132 days to reach the same number of cities in China (see methods in Supplementary Materials). The number of cities providing first reports of COVID-19 peaked at 59 per day on 23 January, the date of the Wuhan travel ban.

The total number of cases reported from each province by 30 January, one week after the Wuhan shutdown, was strongly associated with the total number of travellers from Wuhan (r= 0.98,P< 0.01; excluding Hubei,r= 0.69,P< 0.01) (Fig. 2, B and C). COVID-19 arrived sooner in those cities that had larger populations and had more travellers from Wuhan (Table 1 and table S1). However, the Wuhan travel ban was associated with a delayed arrival time of COVID-19 in other cities by an estimated 2.91 days (95%CI: 2.54 to 3.29 days) on average (Fig. 3B and Table 1).

This delay provided extra time to prepare for the arrival of COVID-19 in more than 130 cities across China but would not have curbed transmission after infection had been exported to new locations from Wuhan.Figure 1 shows the timing and implementation of emergency control measures in 342 cities across China (see also figs. S2 and S4). School closure, the isolation of suspected and confirmed patients, plus the disclosure of information was implemented in all cities. Public gatherings were banned and entertainment venues closed in 220 cities (64.3%). Intra-city public transport was suspended in 136 cities (39.7%) and inter-city travel was prohibited by 219 cities (64.0%). All three measures were applied in 136 cities (table S2).

Cities that implemented a Level 1 response (any combination of control measures) (figs. S2 and S4) pre-emptively, before discovering any COVID-19 cases, reported 33.3% (95%CI: 11.1-44.4%) fewer laboratory-confirmed cases during the first week of their outbreaks (13.0, 95%CI: 7.1-18.8,n= 125) compared with cities that started control later (20.6 cases, 95%CI: 14.5-26.8,n= 171), with a statistically significant difference between the two groups (Mann-WhitneyU= 8197,z= –3.4,P< 0.01). A separate analysis using regression models shows that, among specific control measures, cities that suspended intra-city public transport and/or closed entertainment venues and banned public gatherings, and did so sooner, had fewer cases during the first week of their outbreaks (Table 2 and table S3). This analysis provided no evidence that the prohibition of travel between cities, which was implemented after and in addition to the Wuhan shutdown on 23 January, reduced the number of cases in other cities across China. These results are robust to the choice of statistical regression model (table S3).

The reported daily incidence of confirmed cases peaked in Hubei province (including Wuhan) on 4 February (3156 laboratory-confirmed cases, 5.33/100,000 population in Hubei), and in all other provinces on 31 January (875 cases, 0.07/100,000 population) (fig. S1). The low level of peak incidence per capita, the early timing of the peak, and the subsequent decline in daily case reports, suggest that transmission control measures were not only associated with a delay in the growth of the epidemic, but also with a marked reduction in the number of cases. By fitting an epidemic model to the time series of cases reported in each province (fig. S3), we estimate that the (basic) case reproduction number (R0) was 3.15 prior to the implementation of the emergency response on 23 January (Table 3). As control was scaled-up from 23 January onwards (stage 1), the case reproduction number declined to 0.97, 2.01 and 3.05 (estimated asC1R0) in three groups of provinces, depending on the rate of implementation in each group (Table 3 and table S4). Once the implementation of interventions was 95% complete everywhere (stage 2), the case reproduction number had fallen to 0.04 on average (C2R0), far below the replacement rate (<< 1) and consistent with the rapid decline in incidence (Fig. 4A,Table 3, fig. S5, and table S4).

Based on the fit of the model to daily case reports from each province, and on the preceding statistical analyses, we investigated the possible effects of control measures on the trajectory of the epidemic outside Wuhan city (Fig. 4B). Our model suggests that, without the Wuhan travel ban or the national emergency response, there would have been 744,000 (± 156,000) confirmed COVID-19 cases outside Wuhan by 19 February, day 50 of the epidemic. With the Wuhan travel ban alone, this number would have decreased to 202,000 (± 10,000) cases. With the national emergency response alone (without the Wuhan travel ban), the number of cases would have decreased to 199,000 (± 8500). Thus, neither of these interventions would, on their own, have reversed the rise in incidence by 19 February (Fig. 4B). But together and interactively, these control measures offer an explanation of why the rise in incidence was halted and reversed, limiting the number of confirmed cases reported to 29,839 (fitted model estimate 28,000 ± 1400 cases), 96% fewer than expected in the absence of interventions.

In summary, this analysis shows that transmission control (non-pharmaceutical) measures initiated during Chinese Spring Festival holiday, including the unprecedented Wuhan city travel ban and the Level 1 national emergency response, were strongly associated with, though not necessarily the cause of, a delay in epidemic growth and a reduction in case numbers during the first 50 days of the COVID-19 epidemic in China.

The number of people who have developed COVID-19 during this epidemic, and therefore the number of people who were protected by control measures, is not known precisely, given that cases were almost certainly under-reported. However, in view of the small fraction of people known to have been infected by 19 February (75,532 cases, 5.41 per 100,000 population), it is unlikely that the spread of infection was halted and epidemic growth reversed because the supply of susceptible people had been exhausted. This implies that a large fraction of the Chinese population remains at risk of COVID-19; control measures may need to be reinstated, in some form, if there is a resurgence of transmission. Further investigations are needed to verify that proposition, and population surveys of infection are needed to reveal the true number of people who have been exposed to this novel coronavirus.

We could not investigate the impact of all elements of the national emergency response because many were introduced simultaneously across China. However, this analysis shows that suspending intra-city public transport, closing entertainment venues and banning public gatherings, which were introduced at different times in different places, were associated with the overall containment of the epidemic. However, other factors are likely to have contributed to control, especially the isolation of suspected and confirmed patients and their contact, and it is not yet clear which parts of the national emergency response were most effective. We did not find evidence for prohibiting travel between cities or provinces reduced the numbers of COVID-19 cases outside Wuhan and Hubei, perhaps because such travel bans were implemented as a response to, rather than in anticipation of, the arrival of COVID-19.

This study has drawn inferences, not from controlled experiments, but from statistical and mathematical analyses of the temporal and spatial variation in case reports, human mobility and transmission control measures. With that caveat, control measures were strongly associated with the containment of COVID-19, potentially averting hundreds of thousands of cases by 19 February, day 50 of the epidemic. Whether the means and the outcomes of control can be replicated outside China, and which of the interventions are most effective, are now under intense investigation as the virus continues to spread worldwide.

Supplementary materials

Materials and Methods

Figs S1 to S7

Tables S1 to S4

References (19–28)

References and Notes

1. ↵
   1. N.Zhu,
   2. D.Zhang,
   3. W.Wang,
   4. X.Li,
   5. B.Yang,
   6. J.Song,
   7. X.Zhao,
   8. B.Huang,
   9. W.Shi,
   10. R.Lu,
   11. P.Niu,
   12. F.Zhan,
   13. X.Ma,
   14. D.Wang,
   15. W.Xu,
   16. G.Wu,
   17. G. F.Gao,
   18. W.Tan;China Novel Coronavirus Investigating and Research Team
   ,A novel coronavirus from patients with pneumonia in China, 2019.N. Engl. J. Med.382,727–733(2020).doi:10.1056/NEJMoa2001017pmid:31978945 
2. ↵
   1. R.Lu,
   2. X.Zhao,
   3. J.Li,
   4. P.Niu,
   5. B.Yang,
   6. H.Wu,
   7. W.Wang,
   8. H.Song,
   9. B.Huang,
   10. N.Zhu,
   11. Y.Bi,
   12. X.Ma,
   13. F.Zhan,
   14. L.Wang,
   15. T.Hu,
   16. H.Zhou,
   17. Z.Hu,
   18. W.Zhou,
   19. L.Zhao,
   20. J.Chen,
   21. Y.Meng,
   22. J.Wang,
   23. Y.Lin,
   24. J.Yuan,
   25. Z.Xie,
   26. J.Ma,
   27. W. J.Liu,
   28. D.Wang,
   29. W.Xu,
   30. E. C.Holmes,
   31. G. F.Gao,
   32. G.Wu,
   33. W.Chen,
   34. W.Shi,
   35. W.Tan
   ,Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding.Lancet395,565–574(2020).doi:10.1016/S0140-6736(20)30251-8pmid:32007145 
3. ↵
   1. F.Wu,
   2. S.Zhao,
   3. B.Yu,
   4. Y.-M.Chen,
   5. W.Wang,
   6. Z.-G.Song,
   7. Y.Hu,
   8. Z.-W.Tao,
   9. J.-H.Tian,
   10. Y.-Y.Pei,
   11. M.-L.Yuan,
   12. Y.-L.Zhang,
   13. F.-H.Dai,
   14. Y.Liu,
   15. Q.-M.Wang,
   16. J.-J.Zheng,
   17. L.Xu,
   18. E. C.Holmes,
   19. Y.-Z.Zhang
   ,A new coronavirus associated with human respiratory disease in China.Nature579,265–269(2020).doi:10.1038/s41586-020-2008-3pmid:32015508 
4. ↵
   1. P.Zhou,
   2. X.-L.Yang,
   3. X.-G.Wang,
   4. B.Hu,
   5. L.Zhang,
   6. W.Zhang,
   7. H.-R.Si,
   8. Y.Zhu,
   9. B.Li,
   10. C.-L.Huang,
   11. H.-D.Chen,
   12. J.Chen,
   13. Y.Luo,
   14. H.Guo,
   15. R.-D.Jiang,
   16. M.-Q.Liu,
   17. Y.Chen,
   18. X.-R.Shen,
   19. X.Wang,
   20. X.-S.Zheng,
   21. K.Zhao,
   22. Q.-J.Chen,
   23. F.Deng,
   24. L.-L.Liu,
   25. B.Yan,
   26. F.-X.Zhan,
   27. Y.-Y.Wang,
   28. G.-F.Xiao,
   29. Z.-L.Shi
   ,A pneumonia outbreak associated with a new coronavirus of probable bat origin.Nature579,270–273(2020).doi:10.1038/s41586-020-2012-7pmid:32015507 
5. ↵
   1. J.Cai,
   2. B.Xu,
   3. K. K. Y.Chan,
   4. X.Zhang,
   5. B.Zhang,
   6. Z.Chen,
   7. B.Xu
   ,Roles of Different Transport Modes in the Spatial Spread of the 2009 Influenza A(H1N1) Pandemic in Mainland China.Int. J. Environ. Res. Public Health16,222(2019).doi:10.3390/ijerph16020222pmid:30646629 
6. ↵
   1. C.Wang,
   2. P. W.Horby,
   3. F. G.Hayden,
   4. G. F.Gao
   ,A novel coronavirus outbreak of global health concern.Lancet395,470–473(2020).doi:10.1016/S0140-6736(20)30185-9pmid:31986257 
7. ↵
   1. S.Chen,
   2. J.Yang,
   3. W.Yang,
   4. C.Wang,
   5. T.Bärnighausen
   ,COVID-19 control in China during mass population movements at New Year.Lancet395,764–766(2020).doi:10.1016/S0140-6736(20)30421-9pmid:32105609 
8. ↵
   1. M. U. G.Kraemer,
   2. C.-H.Yang,
   3. B.Gutierrez,
   4. C.-H.Wu,
   5. B.Klein,
   6. D. M.Pigott,
   7. L.du Plessis,
   8. N. R.Faria,
   9. R.Li,
   10. W. P.Hanage,
   11. J. S.Brownstein,
   12. M.Layan,
   13. A.Vespignani,
   14. H.Tian,
   15. C.Dye,
   16. O. G.Pybus,
   17. S. V.Scarpino
   ,The effect of human mobility and control measures on the COVID-19 epidemic in China.Scienceeabb4218(2020).doi:10.1126/science.abb4218 
9. ↵Chinadaily, “Tibet activates highest-level public health alert” (Jan 30, 2020).https://www.chinadaily.com.cn/a/202001/29/WS5e318a36a3101282172739c1.html. 
10. ↵
    1. B. T.Grenfell,
    2. O. N.Bjørnstad,
    3. J.Kappey
    ,Travelling waves and spatial hierarchies in measles epidemics.Nature414,716–723(2001).doi:10.1038/414716apmid:11742391 
11. 1. D.Brockmann,
    2. D.Helbing
    ,The hidden geometry of complex, network-driven contagion phenomena.Science342,1337–1342(2013).doi:10.1126/science.1245200pmid:24337289 
12. 1. A.Wesolowski,
    2. N.Eagle,
    3. A. J.Tatem,
    4. D. L.Smith,
    5. A. M.Noor,
    6. R. W.Snow,
    7. C. O.Buckee
    ,Quantifying the impact of human mobility on malaria.Science338,267–270(2012).doi:10.1126/science.1223467pmid:23066082 
13. 1. N. M.Ferguson,
    2. D. A.Cummings,
    3. S.Cauchemez,
    4. C.Fraser,
    5. S.Riley,
    6. A.Meeyai,
    7. S.Iamsirithaworn,
    8. D. S.Burke
    ,Strategies for containing an emerging influenza pandemic in Southeast Asia.Nature437,209–214(2005).doi:10.1038/nature04017pmid:16079797 
14. 1. K. E.Jones,
    2. N. G.Patel,
    3. M. A.Levy,
    4. A.Storeygard,
    5. D.Balk,
    6. J. L.Gittleman,
    7. P.Daszak
    ,Global trends in emerging infectious diseases.Nature451,990–993(2008).doi:10.1038/nature06536pmid:18288193 
15. ↵
    1. D. M.Morens,
    2. G. K.Folkers,
    3. A. S.Fauci
    ,The challenge of emerging and re-emerging infectious diseases.Nature430,242–249(2004).doi:10.1038/nature02759pmid:15241422 
16. ↵
    1. C.Viboud,
    2. O. N.Bjørnstad,
    3. D. L.Smith,
    4. L.Simonsen,
    5. M. A.Miller,
    6. B. T.Grenfell
    ,Synchrony, waves, and spatial hierarchies in the spread of influenza.Science312,447–451(2006).doi:10.1126/science.1125237pmid:16574822 
17. ↵
    1. A.Wesolowski,
    2. T.Qureshi,
    3. M. F.Boni,
    4. P. R.Sundsøy,
    5. M. A.Johansson,
    6. S. B.Rasheed,
    7. K.Engø-Monsen,
    8. C. O.Buckee
    ,Impact of human mobility on the emergence of dengue epidemics in Pakistan.Proc. Natl. Acad. Sci. U.S.A.112,11887–11892(2015).doi:10.1073/pnas.1504964112pmid:26351662 
18. ↵H. Tian, Y. Liu, Y. Li, C.-H. Wu, B. Chen, M. U. Kraemer, B. Li, J. Cai, B. Xu, Q. Yang, B. Wang, P. Yang, Y. Cui, Y. Song, P. Zheng, Q. Wang, O. N. Bjornstad, R. Yang, B. T. Grenfell, O. G. Pybus, C. Dye, Code for: An investigation of transmission control measures during the first 50 days of the COVID-19 epidemic in China.Zenodo(2020);doi:10.5281/zenodo.3727336 
19. ↵
    1. J. T.Wu,
    2. K.Leung,
    3. G. M.Leung
    ,Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: A modelling study.Lancet395,689–697(2020).doi:10.1016/S0140-6736(20)30260-9pmid:32014114 
20. B. Ripley, B. Venables, D. M. Bates, K. Hornik, A. Gebhardt, D. Firth, M. B. Ripley, Package ‘MASS’.CRAN Repository, See http://cran.r-projectorg/web/packages/MASS/MASS.pdf (2013). 
21. B. R. Kirkwood, J. A. Sterne,Essential medical statistics. (John Wiley & Sons, 2010). 
22. D. E. Bailey,Probability and Statistics. (John Wiley & Sons, 1971). 
23. A. Canty, B. Ripley, boot: Bootstrap R (S-Plus) functions. R package version 1.3-24. (2019). 
24. A. C. Davison, D. V. Hinkley,Bootstrap methods and their application. (Cambridge university press, Cambridge, 1997). 
25. R. M. Anderson, R. M. May,Infectious Diseases of Humans: Dynamics and Control. (Oxford Univ Press, Oxford, 1992). 
26. 1. A.Morton,
    2. B. F.Finkenstädt
    ,Discrete time modelling of disease incidence time series by using Markov chain Monte Carlo methods.J. R. Stat. Soc.54,575–594(2005).doi:10.1111/j.1467-9876.2005.05366.x 
27. 1. Q.Li,
    2. X.Guan,
    3. P.Wu,
    4. X.Wang,
    5. L.Zhou,
    6. Y.Tong,
    7. R.Ren,
    8. K. S. M.Leung,
    9. E. H. Y.Lau,
    10. J. Y.Wong,
    11. X.Xing,
    12. N.Xiang,
    13. Y.Wu,
    14. C.Li,
    15. Q.Chen,
    16. D.Li,
    17. T.Liu,
    18. J.Zhao,
    19. M.Liu,
    20. W.Tu,
    21. C.Chen,
    22. L.Jin,
    23. R.Yang,
    24. Q.Wang,
    25. S.Zhou,
    26. R.Wang,
    27. H.Liu,
    28. Y.Luo,
    29. Y.Liu,
    30. G.Shao,
    31. H.Li,
    32. Z.Tao,
    33. Y.Yang,
    34. Z.Deng,
    35. B.Liu,
    36. Z.Ma,
    37. Y.Zhang,
    38. G.Shi,
    39. T. T. Y.Lam,
    40. J. T.Wu,
    41. G. F.Gao,
    42. B. J.Cowling,
    43. B.Yang,
    44. G. M.Leung,
    45. Z.Feng
    ,Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia.N. Engl. J. Med.382,1199–1207(2020).doi:10.1056/NEJMoa2001316pmid:31995857 
28. ↵
    1. O. N.Bjørnstad,
    2. R. A.Ims,
    3. X.Lambin
    ,Spatial population dynamics: Analyzing patterns and processes of population synchrony.Trends Ecol. Evol.14,427–432(1999).doi:10.1016/S0169-5347(99)01677-8pmid:10511718