Early evaluation of Wuhan City travel restrictions in response to the 2019 novel coronavirus outbreak ===================================================================================================== * Huaiyu Tian * Yidan Li * Yonghong Liu * Moritz Kraemer * Bin Chen * Jun Cai * Bingying Li * Bo Xu * Qiqi Yang * Peng Yang * Yujun Cui * Yimeng Song * Pai Zheng * Quanyi Wang * Ottar Bjornstad * Ruifu Yang * Oliver G. Pybus * Bryan Grenfell * Christopher Dye ## Abstract An ongoing outbreak of a novel coronavirus (2019-nCoV) was first reported in China in December 2019 and has spread to other countries. On January 23rd 2020 China shut down transit in and out of Wuhan city, a major transport hub and conurbation of 11 million inhabitants, to contain the outbreak. By combining epidemiological and human mobility data we find that the travel ban slowed the dispersal of 2019-nCoV from Wuhan to other cities in China by 2.91 days (95% CI: 2.54-3.29). This delay provided extra time to establish and reinforce other control measures that are essential to halt the epidemic. The ongoing diffusion of 2019-nCoV provides an opportunity to examine how travel restrictions impede the spatial dispersal of an emerging infectious disease. **One Sentence Summary** The Wuhan city travel shutdown delayed the dispersal of 2019-nCoV infection to other cities in China Key words * 2019 novel coronavirus * Wuhan * spatiotemporal transmission * China ## Main text At end of December 2019, less than a month before Chinese New Year (Spring Festival), a cluster of pneumonia cases caused by an unknown pathogen were reported in Wuhan City, the largest transport hub in Central China. A novel coronavirus provisionally named 2019-nCoV 1 has since been identified as the etiological agent. Human-to-human transmission of 2019-nCoV has been confirmed 2,3 and by January 28th 2020 all provinces in China except Tibet had reported 2019-nCoV infections, and the population at risk of infection sin the country stands at >1.2 billion people. The increasing movement of people for Chinese New Year is expected to spread the virus further throughout cities in China and elsewhere. To prevent further diffusion of 2019-nCoV from its source, Wuhan prohibited all transport in and out of the city as of 10:00 on January 23rd. To our knowledge, this is the largest attempted movement restriction or quarantine in human history to prevent infectious disease spread. By January 25th, 30 provinces in China had raised their public health response level to the highest state of emergency (level-1). Here we present an analysis and quantification of the consequences and importance of the Wuhan travel prohibition on the ongoing spread of 2019-nCoV across China. Assessment and measurement of the effects of large-scale interventions are crucial for the design of efficient responses against this and future epidemics. To this end, we collected all case reports for 262 cities across China and noted the onset time of the first case of 2019-nCoV (arrival time) and the timing of the local response. We analyse these dates together with high-resolution data on human movement among cities in China, obtained from a large dataset of geolocated mobile phone records, spanning 2018, from the Tencent network. This dataset describes patterns of human mobility across China before 2019-nCoV was discovered. In order to quantify the effect of the Wuhan travel shutdown on 2019-nCoV epidemic spread we analyzed the arrival time of 2019-nCoV from Wuhan to each city as a function of geographic distance (between city centers) and of human movement by air, train, and road (as recorded by Tencent’s location-based services database). Spatial spread of 2019-nCoV (Fig. 1A) was rapid, with 262 cities reporting cases within only 28 days (for comparison, the 2009-H1N1pdm took 132 days to reach the current extent of 219-nCoV in China). Most cities with early arrival dates were in southeast China among which there is greater mobility and higher population density. The rate at which cities first reported 2019-nCoV peaked on January 23rd (the day of the Wuhan travel ban), with 60 reports, after which the spatial diffusion of 2019-nCoV slowed. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/11/2020.01.30.20019844/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/02/11/2020.01.30.20019844/F1) Figure 1. Spatial diffusion of 2019-nCoV in China. (A) Cumulative number of cities reporting disease. Arrival days, defined as the time interval (days) from the date of the first case in the first infected city (Wuhan) to the date of the first case in each newly infected city (a total of 262 cities), to characterize the inter-city transmission rate of 2019-nCoV and 2009-H1N1pdm, respectively. Dashed line shows the date of Wuhan shutdown. (B) Before (blue) and after (red) the intervention. The blue line and points show the fitted regression of arrival times up to the shutdown on day 24 (January 23rd, vertical dashed line). Grey points show the expected arrival times after day 24, without the shutdown. The red line and points show the fitted regression of delayed arrival times after the shutdown on day 24. Each observation (point) represents one city. Error bars give ±2 standard deviations. (C) Map of arrival time delayed by the shutdown of Wuhan. Colors represent the change in arrival time (days) after January 23rd. The arrival time is estimated using the data before the shutdown of Wuhan. (D) Human movement outflows from Wuhan city to other cities in 2018. The warmer and thick lines denote higher volume of outflows (high-connectivity) while the cool and thin lines denote a lower volume of outflows (low-connectivity). Table 1 shows that the Wuhan travel intervention slowed disease spread. As expected, the time of arrival of 2019-nCoV in each city increased with distance from Wuhan, and decreased with passenger flow from Wuhan. Thus the epidemic arrived sooner in those cities that had larger population and had more travelers from Wuhan. On average, the Wuhan shutdown delayed the arrival time of 2019-nCoV in other cities by 2.91 days (95%CI: 2.54-3.29 days), added to the expected arrival time of 22.3 days (95%CI: 20.6-24.0 days) without the shutdown (Fig. 1B). More than 130 cities, covering more than half geographic area and population of China, benefited from this intervention (Fig. 1C). View this table: [Table 1.](http://medrxiv.org/content/early/2020/02/11/2020.01.30.20019844/T1) Table 1. Estimating the impact of the Wuhan travel ban on 2019-nCoV diffusion to other cities in China. Our analysis evaluates a unique intervention against an emerging infectious disease – the cessation of travel from a large, well-connected city in an industrialized country (Fig. 1D). We find that this intervention was effective in slowing 2019-nCoV invasion of new locations. However, other measures need to be reinforced to halt the ongoing epidemic. The Wuhan travel ban provided extra time to make progress on these responses. In addition to governmental responses, public awareness is critical for controlling the spread of this novel coronavirus. Early detection of 2019-nCoV cases in new locations in China is needed in the coming weeks and months to prevent other cities becoming major exporters of 2019-nCoV. China is better equipped than in the past to meet this challenge and a direct reporting system for notifiable epidemic diseases was established in the country after the 2003 SARS outbreak. However, the experience of 2019-nCoV in China suggests that urbanization and the development of modern transport systems increase the urgency of control measures against emerging infectious diseases, as indicated by the faster spatial spread of 2019-nCoV than H1N1pdm in China. While emerging data from mainland China on 2019-nCoV epidemiology and virus genomics have generated important insights into the origin, transmissibility, and severity of this unfolding epidemic, 1,4 many uncertainties remain.5 For example, the role of mobility network structure on disease diffusion is unknown, and the degree to which the delay in spatial spread (due to the Wuhan travel ban) will impact on the final size of the epidemic is unclear. Our analysis is a preliminary step and further analyses and models are required to evaluate the impact of quarantine and other interventions. ## Materials and Methods ### Epidemiological and demographic data We collected raw data from the daily official reports of the health commission of 34 provincial-level administrative units and 341 city-level units. Only laboratory-confirmed cases were used. We constructed a real-time database recording the date of the first reported case in all newly-infected cities with daily updates from December 31st 2019 to January 28th 2020. Population sizes for each city were collected from the China City Statistical Yearbook ([http://olap.epsnet.com.cn/](http://olap.epsnet.com.cn/)). The population sizes recorded for 2018 were used. We calculated the distance between Wuhan and each city reporting 2019-nCoV cases. The location of each city is geocoded by the latitude and longitude coordinates of the city center. ### Human mobility data Human movement can be observed directly from mobile phone data, through the location-based services (LBS) employed by popular Tencent applications, such as WeChat and QQ. Average movement outflows from Wuhan City to other cities, by air, train, and road, were calculated from the migration flows database ([https://heat.qq.com/](https://heat.qq.com/)) over the entire 2018. ### Statistical model The association between distance, human movement, interventions and epidemic timing of 2019-nCoV was assessed with a regression analysis using a General Linear Model framework (GLM). The best model that emerged from this analysis is: ![Formula][1] where *TotalFlow**j* represents the average passenger volume from Wuhan City to city *j* by airplane, train, and road. *Pop**j* is the population in city *j. Lat**j* and *Lon* *j* represent *latitude* and *longitude* of city *j*. The dummy variable *shutdown**j* is used to identify whether the arrival time of the newly infected city *j* is influenced by the shutdown of Wuhan, where 0 represents no-intervention and 1 represents intervention. The dependent variable *Y**j* is the arrival time of the epidemic in city *j*, which measures the spatial spread of 2019-nCoV. The *β**i* are the regression coefficients. ## Data Availability Population sizes for each city were collected from the China City Statistical Yearbook (http://olap.epsnet.com.cn/).Human movement can be observed directly from mobile phone data, through the location-based services (LBS) employed by popular Tencent applications, such as WeChat and QQ. Average movement outflows from Wuhan City to other cities, by air, train, and road, were calculated from the migration flows database (https://heat.qq.com/) over the entire 2018. ## Acknowledgements We thank the many thousands of CDC staff and local health workers in China who collected data and continue to strive to contain 2019-nCoV in China and elsewhere. Funding for this study was provided by the Beijing Natural Science Foundation (JQ18025); Beijing Advanced Innovation Program for Land Surface Science; National Natural Science Foundation of China (81673234); Young Elite Scientist Sponsorship Program by CAST (YESS)(2018QNRC001); HT and OGP acknowledge support from the Oxford Martin School. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. * Received January 30, 2020. * Revision received February 11, 2020. * Accepted February 11, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Zhu N, Zhang D, Wang W, et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. New Eng J Med 2020:DOI: 10.1056/NEJMoa2001017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31978945&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F11%2F2020.01.30.20019844.atom) 2. 2.Cai J, Xu B, Chan KKY, et al. Roles of Different Transport Modes in the Spatial Spread of the 2009 Influenza A (H1N1) Pandemic in Mainland China. Int J Environ Res Public Health 2019;16:222. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph16020222&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30646629&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F11%2F2020.01.30.20019844.atom) 3. 3.Wang C, Horby PW, Hayden FG, Gao GF. A novel coronavirus outbreak of global health concern. Lancet 2020:DOI:[https://doi.org/10.1016/S0140-6736(20)30185-9](https://doi.org/10.1016/S0140-6736(20)30185-9). 4. 4.Chan JFW, Yuan S, Kok KH, et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet 2020:DOI:[https://doi.org/10.1016/S0140-6736(20)30154-9](https://doi.org/10.1016/S0140-6736(20)30154-9). 5. 5.Heymann DL. Data sharing and outbreaks: best practice exemplified. Lancet 2020:DOI:[https://doi.org/10.1016/S0140-6736(20)30184-7](https://doi.org/10.1016/S0140-6736(20)30184-7). [1]: /embed/graphic-3.gif