# The effect of human mobility and control measures on the COVID-19 epidemic in China

Moritz U G Kraemer, Chia-Hung Yang, Bernardo Gutierrez, Chieh-Hsi Wu, Brennan Klein, David M Pigott, Louis Du Plessis, Nuno Faria, Ruoran Li, William P Hanage, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Moritz U G Kraemer, Chia-Hung Yang, Bernardo Gutierrez, Chieh-Hsi Wu, Brennan Klein, et al. The effect of human mobility and control measures on the COVID-19 epidemic in China. Science, 2020, 368 (6490), pp.493-497. 10.1126/science.abb4218. hal-02614055

HAL Id: hal-02614055

https://hal.sorbonne-universite.fr/hal-02614055v1

Submitted on 18 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
Cite as: M. U. G. Kraemer et al., *Science*
10.1126/science.abb4218 (2020).

# The effect of human mobility and control measures on the COVID-19 epidemic in China

Moritz U. G. Kraemer<sup>1,2,3*</sup>, Chia-Hung Yang<sup>4</sup>, Bernardo Gutierrez<sup>1,5</sup>, Chieh-Hsi Wu<sup>6</sup>, Brennan Klein<sup>4</sup>, David M. Pigott<sup>7</sup>, Open COVID-19 Data Working Group<sup>†</sup>, Louis du Plessis<sup>1</sup>, Nuno R. Faria<sup>1</sup>, Ruoran Li<sup>8</sup>, William P. Hanage<sup>8</sup>, John S. Brownstein<sup>2,3</sup>, Maylis Layan<sup>9,10</sup>, Alessandro Vespignani<sup>4,11</sup>, Huaiyu Tian<sup>12</sup>, Christopher Dye<sup>1</sup>, Oliver G. Pybus<sup>1,13*</sup>, Samuel V. Scarpino<sup>4*</sup>

The ongoing COVID-19 outbreak expanded rapidly throughout China. Major behavioral, clinical, and state interventions have been undertaken to mitigate the epidemic and prevent the persistence of the virus in human populations in China and worldwide. It remains unclear how these unprecedented interventions, including travel restrictions, affected COVID-19 spread in China. We use real-time mobility data from Wuhan and detailed case data including travel history to elucidate the role of case importation on transmission in cities across China and ascertain the impact of control measures. Early on, the spatial distribution of COVID-19 cases in China was explained well by human mobility data. Following the implementation of control measures, this correlation dropped and growth rates became negative in most locations, although shifts in the demographics of reported cases were still indicative of local chains of transmission outside Wuhan. This study shows that the drastic control measures implemented in China substantially mitigated the spread of COVID-19.

The outbreak of COVID-19 spread rapidly from its origin in Wuhan, Hubei Province, China (1). A range of interventions have been implemented following the detection in late December 2019 of a cluster of pneumonia cases of unknown etiology, and identification of the causative virus SARS-CoV-2 in early January 2020 (2). Interventions include improved rates of diagnostic testing, clinical management, rapid isolation of suspected cases, confirmed cases and contacts and, most notably, restrictions on mobility (hereafter called cordon sanitaire) imposed on Wuhan city on 23rd January. Travel restrictions were subsequently imposed on 14 other cities across Hubei Province and partial movement restrictions were enacted in many cities across China. Initial analysis suggests that the Wuhan cordon sanitaire resulted in an average delay of COVID-19 spread to other cities of 3 days (3), but the full extent of the effect of the mobility restrictions and other types of interventions on transmission has not been examined quantitatively (4–6). Questions remain over how these interventions affected the spread of SARS-CoV-2 to locations outside of Wuhan. We here use real-time mobility data, crowdsourced line-list data of cases with reported travel

history, and timelines of reporting changes to identify early shifts in the epidemiological dynamics of the COVID-19 epidemic in China, from an epidemic driven by frequent importations to local transmission.

## Human mobility predicts the spread and size of epidemics in China

As of 1st March 2020, 79,986 cases of COVID-19 were confirmed in China (Fig. 1a) (7). Reports of cases in China were mostly restricted to Hubei until 23rd January 2020 (81% of all cases), after which most provinces reported rapid increases in cases (Fig. 1a). We built a line list dataset from reported cases in China with information on travel history and demographic characteristics (8). We note that the majority of early cases (before 23rd January 2020, Materials and Methods) reported outside of Wuhan had known travel history to Wuhan (57%) and were distributed across China (Fig. 1b), highlighting the importance of Wuhan as a major source of early cases. However, initial testing was focused mainly on travelers from Wuhan, potentially biasing estimates of travel related infections upwards (Materials and Methods). Among cases known

<footnote n="1"><sup>1</sup>Department of Zoology, University of Oxford, Oxford, UK. <sup>2</sup>Harvard Medical School, Harvard University, Boston, MA, USA. <sup>3</sup>Boston Children's Hospital, Boston, MA, USA.</footnote>

<footnote n="4"><sup>4</sup>Network Science Institute, Northeastern University, Boston, MA, USA. <sup>5</sup>School of Biological and Environmental Sciences, Universidad San Francisco de Quito USFQ, Quito, Ecuador. <sup>6</sup>Mathematical Sciences, University of Southampton, Southampton, UK. <sup>7</sup>Institute for Health Metrics and Evaluation, Department of Health Metrics, University of Washington, Seattle, WA, USA. <sup>8</sup>Harvard T.H. Chan School of Public Health, Boston, MA, USA. <sup>9</sup>Mathematical Modelling of Infectious Diseases Unit, Institut Pasteur, UMR2000, CNRS, Paris, France. <sup>10</sup>Sorbonne Universite, Paris, France. <sup>11</sup>ISI Foundation, Turin, Italy. <sup>12</sup>State Key Laboratory of Remote Sensing Science, College of Global Change and Earth System Science, Beijing Normal University, Beijing, China. <sup>13</sup>Department of Pathobiology and Population Sciences, The Royal Veterinary College, London, UK.</footnote>

<footnote>*Corresponding author. Email: s.scarpino@northeastern.edu (S.V.S.); oliver.pybus@zoo.ox.ac.uk (O.G.P.); moritz.kraemer@zoo.ox.ac.uk (M.U.G.K.)</footnote>

<footnote>†Members of the Open COVID-19 Data Working Group are listed in the supplementary materials.</footnote>
to have traveled from Wuhan before 23rd January 2020, the time from symptom onset to confirmation was 6.5 days (SD: 4.2; fig. S2), providing opportunity for onward transmission at the destination. More active surveillance reduced this interval to 4.8 days (SD: 3.03; fig. S2) for those who travelled after 23rd January 2020.

To identify accurately a timeframe for evaluating early shifts in SARS-CoV-2 transmission in China, we first estimated from case data the average incubation period of COVID-19 infection (i.e. the duration between time of infection and symptom onset (9, 10)). Since infection events are typically not directly observed, we estimate incubation period from the span of exposure during which infection likely occurred. Using detailed information on 38 cases for whom both the dates of entry to and exit from Wuhan are known, we estimate the mean incubation period to be 5.1 days (std. dev. = 3.0 days; fig. S1), similar to previous estimates from other data (11, 12). In subsequent analyses we add an upper estimate of one incubation period (mean + 1 standard deviation = 8 days) to the date of Wuhan shutdown, in order to delineate the date before which cases recorded in other provinces might represent infections acquired in Hubei (i.e., 1st February 2020; Fig. 1a).

In order to understand whether the volume of travel within China could predict the epidemic outside of Wuhan, we analyzed real-time human mobility data from Baidu Inc., together with epidemiological data from each province (Materials and Methods). We investigated spatio-temporal disease spread to elucidate the relative contribution of Wuhan to transmission elsewhere and evaluate how the cordon sanitaire may have impacted it.

Among cases reported outside Hubei province in our dataset, we observe 515 cases with known travel history to Wuhan and a symptom onset date before 31st January 2020, compared with only 39 after 31st January, 2020, illustrating the effect of travel restrictions (Figs. 1b and 2a and fig. S3). We confirm the expected decline of importation with real-time human mobility data from Baidu Inc. Movements of individuals out of Wuhan increased in the days before the Lunar New Year and the establishment of the cordon sanitaire, before rapidly decreasing to almost no movement (Fig. 2, a and b). The travel ban appears to have prevented travel in and out of Wuhan around the time of the Lunar New Year celebration (Fig. 2a) and likely reduced further dissemination of SARS-CoV-2 from Wuhan.

To test the contribution of the epidemic in Wuhan to seeding epidemics elsewhere in China we build a naïve COVID-19 GLM (13) model of daily case counts (Materials and Methods). We estimate the epidemic doubling time outside Hubei to be 4.0 days (range across provinces of 3.6 - 5.0 days) and estimate the epidemic doubling time within Hubei to be 7.2 days, consistent with previous reports (5, 12, 14, 15). Our

model predicts daily case counts across all provinces with relatively high accuracy (as measured with a pseudo-R² from a negative binomial GLM) throughout early February 2020, and when accounting for human mobility (Fig. 2c and tables S1 and S2), consistent with an exploratory analysis (6).

We find that the magnitude of the early epidemic (total number of cases until February 10, 2020) outside of Wuhan is remarkably well predicted by the volume of human movement out of Wuhan alone (R² = 0.89 from a log-linear regression using cumulative cases; fig. S8). Therefore cases exported from Wuhan prior to the cordon sanitaire appear to have contributed to initiating local chains of transmission, both in neighboring provinces (e.g., Henan) and in more distant provinces, (e.g., Guangdong and Zhejiang; Figs. 1a and 2b). Further, the frequency of introductions from Wuhan are also predictive of the size of the early epidemic in other provinces (controlling for population size) and thus the probability of large outbreaks (fig. S8).

After 1st February 2020 (corresponding to one mean + one SD incubation period after the cordon sanitaire and other interventions were implemented), the correlation of daily case counts and human mobility from Wuhan decreased (Fig. 2c), indicating that variability among locations in daily case counts was better explained by factors unrelated to human mobility, such as local public health response. This suggests that while travel restrictions may have reduced the flow of case importations from Wuhan, other local mitigation strategies aimed at halting local transmission increased in importance later.

We estimate also the growth rates of the epidemic in all other provinces (Materials and Methods). Interestingly, we find that all provinces outside Hubei experienced faster growth rates between January 9th – January 22nd, 2020 (Fig. 3, a and b, and fig. S4b) which was the time before travel restrictions and substantial control measures were implemented (Fig. 3c and fig. S6); this is also apparent from the case counts by province (fig. S6). In the same period, variation in the growth rates are almost entirely explained by human movements from Wuhan (Fig. 3c and fig. S9), consistent with theory of infectious disease spread in highly coupled metapopulations (16, 17). Following the implementation of drastic control measures across the country, growth rates become negative (Fig. 3b), indicating that transmission was successfully mitigated. The correlation of growth rates and human mobility from Wuhan becomes negative, i.e., provinces with larger mobility from Wuhan prior to the cordon sanitaire (but also larger number of cases overall) have more rapidly declining growth rates of daily case counts. This could be due partly to travel restrictions but also to the fact that control measures may have been more drastic in locations with larger outbreaks driven by local transmission (see more detail in section “Current role of imported cases in Chinese provinces”).

The travel ban coincides with increased testing capacity across provinces in China. An alternative hypothesis is that
the observed epidemiological patterns outside Wuhan are the result of increased testing capacity. We test this hypothesis by including differences in testing capacity before and after the rollout of large-scale testing in China on 20th of January 2020 (the date that COVID-19 became a Class B notifiable disease (18, 19)) and test the impact of this binary variable on the predictability of daily cases (Materials and Methods). We plot the relative improvement in the prediction of our model (based on normalized residual error) of (i) a model that includes daily mobility from Wuhan and (ii) a model that includes testing availability (see more details in Materials and Methods). Overall, the inclusion of mobility data from Wuhan produces a significant improvement in the model's prediction (delta-BIC > 250, (20)) over a naive model that considers only autochthonous transmission with a doubling time of 2-8 days (Fig. 3b). Of the 27 provinces in China reporting cases through February 6th, 2020, we find that in 12 provinces the largest improvements in prediction can be achieved using mobility only (fig. S5). In 10 provinces, both testing and mobility improve the model's prediction, and in only one province (Hunan) is testing the most important factor improving model prediction (fig. S5). We conclude that laboratory testing during the early phase of the epidemic is critical, however, mobility out of Wuhan remains the main driver of spread prior to the cordon sanitaire. Large scale molecular and serological data will be important to investigate further the exact magnitude of the impact of human mobility compared to other factors.

## Current role of imported cases in Chinese provinces

Since case counts outside Wuhan have decreased (Fig. 3b), we can further investigate the current contribution of imported cases to local epidemics outside Wuhan by investigating case characteristics. Age and sex distributions can reflect heterogeneities in the risk of infection within affected populations. To investigate meaningful shifts in the epidemiology of the COVID-19 outbreak through time, we examined age and sex data for cases from different periods of the outbreak, and from individuals with and without travel from Wuhan. However, details of travel history exist for only a fraction of confirmed cases and this information is particularly scant for some provinces (e.g., Zhejiang and Guangdong). Consequently, we grouped confirmed cases into four categories: (I) early cases with travel history (early = reported before 1st Feb), (II) early cases without travel history, (III) later cases with travel history (later = reported between 1st – 10th Feb), (IV) later cases without travel history.

Using crowdsourced case data, we found that cases with travel history (categories I and III) had similar median ages and sex ratios in both the early and later phases of the outbreak (41 vs 42 years old, 50% interquartile interval: 32.75 vs 30.75 and 54.25 vs 53.5 respectively, p-value > 0.1; 1.47 vs. 1.45

males per female, respectively; Fig. 4d and fig. S7). Early cases with no information on travel history (category II) had a similar median age and sex ratio to those with known travel history (42 years old (50% interquartile interval: 30.5 – 49.5, p-value > 0.1) and 1.80 males per female; Fig. 4d). However, the sex ratio of later cases without reported travel history (category IV) shifted to approximately 1:1 (57 male vs. 62 female, X² test, p-value < 0.01), as expected under a null hypothesis of equal transmission risk (Fig. 4, a, b, and d; see also (21, 22) and the materials and methods) and the median age in this group increased to 46 (50% interquartile interval: 34.25 – 58, t test: p-value < 0.01) (Fig. 4 a, b, and c, and fig. S7). We hypothesize that many of the cases with no known travel history in the early phase were indeed travelers that contributed to disseminating SARS-CoV-2 outside of Wuhan. The shift toward more equal sex ratios and older ages in non-travellers after 31st January 2020 confirm the finding that epidemics outside Wuhan were then driven by local transmission dynamics. The case definition changed to include cases without travel history to Wuhan after 23rd January 2020 (Materials and Methods).

## Discussion

Containment of respiratory infections is particularly difficult if they are characterized by relatively mild symptoms or transmission before the onset of disease (23, 24). Intensive control measures, including travel restrictions, have been implemented to limit the spread of COVID-19 in China. Here, we show that travel restrictions are particularly useful in the early stage of an outbreak when it is confined to a certain area that acts as a major source. However, travel restrictions may be less effective once the outbreak is more widespread. The combination of interventions implemented in China were clearly successful in mitigating spread and reducing local transmission of COVID-19, although in this work it was not possible to definitively determine the impact of each intervention. Much further work is required to determine how to balance optimally the expected positive effect on public health with the negative impact on freedom of movement, the economy, and society at large.

## REFERENCES AND NOTES

1. S. Chen, J. Yang, W. Yang, C. Wang, T. Bärnighausen, COVID-19 control in China during mass population movements at New Year. *Lancet* **395**, 764–766 (2020). doi:10.1016/S0140-6736(20)30421-9 Medline

2. N. Zhu, D. Zhang, W. Wang, X. Li, B. Yang, J. Song, X. Zhao, B. Huang, W. Shi, R. Lu, P. Niu, F. Zhan, X. Ma, D. Wang, W. Xu, G. Wu, G. F. Gao, 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/NEJMoa2001017 Medline

3. H. Tian, Y. Liu, Y. Li, C.-H. Wu, B. Chen, M. U. G. 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. Grenfell, O. Pybus, C. Dye, The impact of transmission control measures during the first 50 days of the COVID-19 epidemic in China. *medRxiv* 2020.01.30.20019844 [Preprint]. 10 March 2020; https://doi.org/10.1101/2020.01.30.20019844
4. Z. Du, L. Wang, S. Cauchemez, X. Xu, X. Wang, B. J. Cowling, L. A. Meyers, Risk for transportation of 2019 novel coronavirus disease from Wuhan to other cities in China. *Emerg. Infect. Dis.* **26** (2020). doi:10.3201/eid2605.200146 *Medline*

5. J. T. Wu, K. Leung, G. M. Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: A modelling study. *Lancet* **395**, 689–697 (2020). doi:10.1016/S0140-6736(20)30260-9 *Medline*

6. S. Zhao, Z. Zhuang, P. Cao, J. Ran, D. Gao, Y. Lou, L. Yang, Y. Cai, W. Wang, D. He, M. H. Wang, Quantifying the association between domestic travel and the exportation of novel coronavirus (2019-nCoV) cases from Wuhan, China in 2020: A correlational analysis. *J. Travel Med.* 10.1093/jtm/taaa022 (2020). doi:10.1093/jtm/taaa022 *Medline*

7. World Health Organization (WHO), *Coronavirus Disease 2019 (COVID-19) Situation Report* – 27 (2020): https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200216-sitrep-27-covid-19.pdf?sfvrsn=78c0eb78_2

8. B. Xu, B. Gutierrez, S. Mekaru, K. Sewalk, L. Goodwin, A. Loskill, E. L. Cohn, Y. Hswen, S. C. Hill, M. M. Cobo, A. E. Zarebski, S. Li, C.-H. Wu, E. Hulland, J. D. Morgan, L. Wang, K. O'Brien, S. V. Scarpino, J. S. Brownstein, O. G. Pybus, D. M. Pigott, M. U. G. Kraemer, Epidemiological data from the COVID-19 outbreak, real-time case information. *Sci. Data* **7**, 106 (2020). doi:10.1038/s41597-020-0448-0

9. S. Cauchemez, C. Fraser, M. D. Van Kerkhove, C. A. Donnelly, S. Riley, A. Rambaut, V. Enouf, S. van der Werf, N. M. Ferguson, Middle East respiratory syndrome coronavirus: Quantification of the extent of the epidemic, surveillance biases, and transmissibility. *Lancet Infect. Dis.* **14**, 50–56 (2014). doi:10.1016/S1473-3099(13)70304-9 *Medline*

10. J. Lessler, N. G. Reich, R. Brookmeyer, T. M. Perl, K. E. Nelson, D. A. T. Cummings, Incubation periods of acute respiratory viral infections: A systematic review. *Lancet Infect. Dis.* **9**, 291–300 (2009). doi:10.1016/S1473-3099(09)70069-6 *Medline*

11. J. A. Backer, D. Klinkenberg, J. Wallinga, Incubation period of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China, 20-28 January 2020. *Euro Surveill.* **25**, 20-28 (2020). doi:10.2807/1560-7917.ES.2020.25.5.2000062 *Medline*

12. Q. Li, X. Guan, P. Wu, X. Wang, L. Zhou, Y. Tong, R. Ren, K. S. M. Leung, E. H. Y. Lau, J. Y. Wong, X. Xing, N. Xiang, Y. Wu, C. Li, Q. Chen, D. Li, T. Liu, J. Zhao, M. Li, W. Tu, C. Chen, L. Jin, R. Yang, Q. Wang, S. Zhou, R. Wang, H. Liu, Y. Luo, Y. Liu, G. Shao, H. Li, Z. Tao, Y. Yang, Z. Deng, B. Liu, Z. Ma, Y. Zhang, G. Shi, T. T. Y. Lam, J. T. K. Wu, G. F. Gao, B. J. Cowling, B. Yang, G. M. Leung, Z. Feng, Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. *N. Engl. J. Med.* NEJMoa2001316 (2020). doi:10.1056/NEJMoa2001316 *Medline*

13. T. J. Hastie, D. Pregibon, "Generalized linear models" in *Statistical Models in S*, J. M. Chambers, T. J. Hastie, Eds. (Wadsworth & Brooks/Cole, 1992), pp. 195–246.

14. J. Riou, C. L. Althaus, Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020. *Euro Surveill.* **25**, 1–5 (2020). doi:10.2807/1560-7917.ES.2020.25.4.2000058 *Medline*

15. A. R. Tuite, D. N. Fisman, Reporting, epidemic growth, and reproduction numbers for the 2019 novel coronavirus (2019-nCoV) epidemic. *Ann. Intern. Med.* (2020). doi:10.7326/M20-0358 *Medline*

16. M. J. Keeling, O. N. Bjørnstad, B. T. Grenfell, "Metapopulation dynamics of infectious diseases" in *Ecology, Genetics and Evolution of Metapopulations*, I. Hanski, O. E. Gaggiotti, Eds. (Elsevier, 2004), pp. 415–445.

17. D. J. Watts, R. Muhamad, D. C. Medina, P. S. Dodds, Multiscale, resurgent epidemics in a hierarchical metapopulation model. *Proc. Natl. Acad. Sci. U.S.A.* **102**, 11157–11162 (2005). doi:10.1073/pnas.0501226102 *Medline*

18. World Health Organization, *Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19) (WHO, 2020): https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf*.

19. Z. Wu, J. M. McGoogan, Characteristics of and important lessons from the coronavirus disease, 2019 (covid-19) outbreak in china: summary of a Report of 72,314 cases from the Chinese Center for Disease Control and Prevention. *JAMA* **219**, 17–20 (2020). doi:10.1001/jama.2020.2648 *Medline*

20. K. P. Burnham, D. R. Anderson, Multimodel inference. *Sociol. Methods Res.* **33**, 261–304 (2004). doi:10.1177/0049124104268644

21. Novel Coronavirus Pneumonia Emergency Response Epidemiology Team, The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China [in Chinese]. *Zhonghua Liu Xing Bing Xue Za Zhi* **41**, 145–151 (2020). *Medline*

22. E. Goldstein, V. E. Pitzer, J. J. O'Hagan, M. Lipsitch, Temporally varying relative risks for infectious diseases: Implications for infectious disease control. *Epidemiology* **28**, 136–144 (2017). doi:10.1097/EDE.0000000000000571 *Medline*

23. C. Fraser, S. Riley, R. M. Anderson, N. M. Ferguson, Factors that make an infectious disease outbreak controllable. *Proc. Natl. Acad. Sci. U.S.A.* **101**, 6146–6151 (2004). doi:10.1073/pnas.0307506101 *Medline*

24. R. Li, S. Pei, B. Chen, Y. Song, T. Zhang, W. Yang, J. Shaman, Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (COVID-19). *medRxiv* 2020.02.14.20023127 [Preprint]. 17 February 2020; https://doi.org/10.1101/2020.02.14.20023127

25. M. U. G. Kraemer, C.-H. Yang, B. Gutierrez, C.-H. Wu, B. Klein, D. Pigott, Open COVID-19 Data Working Group, L. du Plessis, N. R. Faria, R. Li, W. P. Hanage, J. S. Brownstein, M. Layan, A. Vespignani, H. Tian, C. Dye, O. G. Pybus, S. V. Scarpino, Code for: The effect of human mobility and control measures on the COVID-19 epidemic in China. *Zenodo* (2020); https://doi.org/10.5281/zenodo.3714914

26. B. Xu, M. U. G. Kraemer: Open COVID-19 Data Curation Group, Open access epidemiological data from the COVID-19 outbreak. *Lancet Infect. Dis.* **309**, 30119 (2020). doi:10.1016/S1473-3099(20)30119-5 *Medline*

27. R. E. Ramshaw, I. D. Letourneau, A. Y. Hong, J. Hon, J. D. Morgan, J. C. P. Osborne, S. Shirude, M. D. Van Kerkhove, S. I. Hay, D. M. Pigott, A database of geopositioned Middle East Respiratory Syndrome Coronavirus occurrences. *Sci. Data* **6**, 318 (2019). doi:10.1038/s41597-019-0330-0 *Medline*

28. M. L. McHugh, The chi-square test of independence. *Biochem. Med. (Zagreb)* **23**, 143–149 (2013). doi:10.11613/ BM.2013.018 *Medline*

29. J. H. McDonald, *Handbook of Biological Statistics* (Sparky House, ed. 3, 2014).

30. N. R. Faria, M. U. G. Kraemer, S. C. Hill, J. Goes de Jesus, R. S. Aguiar, F. C. M. Iani, J. Xavier, J. Quick, L. du Plessis, S. Dellicour, J. Thézé, R. D. O. Carvalho, G. Baele, C.-H. Wu, P. P. Silveira, M. B. Arruda, M. A. Pereira, G. C. Pereira, J. Lourenço, U. Obolski, L. Abade, T. I. Vasylyeva, M. Giovanetti, D. Yi, D. J. Weiss, G. R. W. Wint, F. M. Shearer, S. Funk, B. Nikolay, V. Fonseca, T. E. R. Adelino, M. A. A. Oliveira, M. V. F. Silva, L. Sacchetto, P. O. Figueiredo, I. M. Rezende, E. M. Mello, R. F. C. Said, D. A. Santos, M. L. Ferraz, M. G. Brito, L. F. Santana, M. T. Menezes, R. M. Brindeiro, A. Tanuri, F. C. P. Dos Santos, M. S. Cunha, J. S. Nogueira, I. M. Rocco, A. C. da Costa, S. C. V. Komninakis, V. Azevedo, A. O. Chieppe, E. S. M. Araujo, M. C. L. Mendonça, C. C. Dos Santos, C. D. Dos Santos, A. M. Mares-Guia, R. M. R. Nogueira, P. C. Sequeira, R. G. Abreu, M. H. O. Garcia, A. L. Abreu, O. Okumoto, E. G. Kroon, C. F. C. de Albuquerque, K. Lewandowski, S. T. Pullan, M. Carroll, T. de Oliveira, E. C. Sabino, R. P. Souza, M. A. Suchard, P. Lemey, G. S. Trindade, B. P. Drumond, A. M. B. Filippis, N. J. Loman, S. Cauchemez, L. C. J. Alcantara, O. G. Pybus, Genomic and epidemiological monitoring of yellow fever virus transmission potential. *Science* **361**, 894–899 (2018). doi:10.1126/science.aat7115 *Medline*

31. M. U. G. Kraemer, R. C. Reiner Jr., O. J. Brady, J. P. Messina, M. Gilbert, D. M. Pigott, D. Yi, K. Johnson, L. Earl, L. B. Marczak, S. Shirude, N. Davis Weaver, D. Bisanzio, T. A. Perkins, S. Lai, X. Lu, P. Jones, G. E. Coelho, R. G. Carvalho, W. Van Bortel, C. Marsboom, G. Hendrickx, F. Schaffner, C. G. Moore, H. H. Nax, L. Bengtsson, E. Wetter, A. J. Tatem, J. S. Brownstein, D. L. Smith, L. Lambrechts, S. Cauchemez, C. Linard, N. R. Faria, O. G. Pybus, T. W. Scott, Q. Liu, H. Yu, G. R. W. Wint, S. I. Hay, N. Golding, Past and future spread of the arbovirus vectors Aedes aegypti and Aedes albopictus. *Nat. Microbiol.* **4**, 854–863 (2019). doi:10.1038/s41564-019-0376-y *Medline*

32. G. Jin, J. Yu, L. Han, S. Duan, The impact of traffic isolation in Wuhan on the spread of 2019-nCoV. *medRxiv* 2020.02.04.20020438 [Preprint]. 5 February 2020; https://doi.org/10.1101/2020.02.04.20020438

33. S. Ai, G. Zhu, F. Tian, H. Li, Y. Gao, Y. Wu, Q. Liu, H. Lin, Population movement, city closure and spatial transmission of the 2019-nCoV infection in China. *medRxiv* 2020.02.04.20020339 [Preprint]. 5 February 2020; https://doi.org/10.1101/2020.02.04.20020339

34. Aurora Big Data, 2017 Mobile Map App Research Report: Which of the Highest, the Baidu, and Tencent Is Strong? [in Chinese] (2017); https://baijiahao.baidu.com/s?id=1590386747028939917&wfr=spider&for=pc.
35. D. Bates, M. Mächler, B. Bolker, S. Walker, Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 201–210 (2015). doi:10.18637/jss.v067.i01

36. J. Friedman, T. Hastie, R. Tibshirani, Regularization paths for generalized linear models via coordinate descent. *J. Stat. Softw.* **33**, 1–22 (2010). doi:10.18637/jss.v033.i01Medline

37. B. Ripley, B. Venables, D. M. Bates, K. Hornik, A. Gebhardt, D. Firth, MASS (R Package Version 7.3-51.5, 2019); https://cran.r-project.org/package=MASS

38. M. U. G. Kraemer, C.-H. Yang, B. Gutierrez, C.-H. Wu, B. Klein, D. M. Pigott, Open COVID-19 Data Working Group, L. du Plessis, N. R. Faria, R. Li, W. P. Hanage, J. S. Brownstein, M. Layan, A. Vespignani, H. Tian, C. Dye, O. G. Pybus, S. V. Scarpino, Code for: The effect of human mobility and control measures on the COVID-19 epidemic in China, Zenodo (2020); https://doi.org/10.5281/zenodo.3714914

39. H. Heinzl, M. Mittlböck, Pseudo R-squared measures for Poisson regression models with over- or underdispersion. *Comput. Stat. Data Anal.* **44**, 253–271 (2003). doi:10.1016/S0167-9473(03)00062-8

40. A. Zeileis, C. Kleiber, S. Jackman, Regression models for count data in R. *J. Stat. Softw.* **27**, 1–25 (2008). doi:10.18637/jss.v027.i08

## ACKNOWLEDGMENTS

We would like to thank all individuals who are collecting epidemiological data of the COVID-19 outbreak around the world. **Funding:** HT, OGP and MUGK acknowledge support from the Oxford Martin School. MUGK is supported by a Branco Weiss Fellowship. BG is supported by a Universities of Academic Excellence Scholarship Program of the Secretariat for Higher Education, Science, Technology, and Innovation of the Republic of Ecuador (ARSEQ-BEC-003163-2017). NRF is supported by a Sir Henry Dale Fellowship. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. WPH was supported by the National Institute of General Medical Sciences (#U54GM088558). **Authors contributions:** MUGK, OGP, SVS developed the idea and research. MUGK and SVS wrote the first draft of the manuscript and all other authors discussed results and edited the manuscript. MUGK, BG, SVS, DMP and the open COVID-19 data working group collected and validated epidemiological data. RL and MUGK collected intervention data. C-HY, BK, and SVS collected and processed human mobility data. **Competing interests:** SVS is on the advisory board for BioFire Diagnostics Trend Surveillance, which includes paid consulting. AV reports past grants and personal fees from Metabiota Inc., outside the submitted work. All other authors declare no competing interests. **Data and materials availability:** Code and data are available on the following GitHub repository: https://github.com/Emergent-Epidemics/covid19_cordon and permanently on Zenodo (25). This work is licensed under a Creative Commons Attribution 4.0 International (CC BY 4.0) license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. This license does not apply to figures/photos/artwork or other content included in the article that is credited to a third party; obtain authorization from the rights holder before using such material.

## SUPPLEMENTARY MATERIALS

science.sciencemag.org/cgi/content/full/science.abb4218/DC1

Materials and Methods

Supplementary Text

Figs. S1 to S9

Tables S1 and S2

List of members of the Open COVID-19 Data Working Group

References (26–40)

3 March 2020; accepted 23 March 2020

Published online 25 March 2020

10.1126/science.abb4218
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Number of cases and key dates during the epidemic. (a) The epidemic curve of the COVID-19 outbreak in provinces in China. Vertical lines and boxes indicate key dates such as implementation of the cordon sanitaire of Wuhan (grey) and the end of the first incubation period after the travel restrictions (red). The thin grey line represents the closure of Wuhan seafood market on 1st January 2020. The width of each horizontal tube represents the number of reported cases in that province. (b) Map of COVID-19 confirmed cases (n = 554) that had reported travel history from Wuhan before travel restrictions were implemented on January 23, 2020. Colors of the arrows indicate date of travel relative to the date of travel restrictions.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>Fig. 2. Human mobility, spread and synchrony of COVID-19 outbreak in China. (a) Human mobility data extracted in real time from Baidu. Travel restrictions from Wuhan and large scale control measures started on January 23, 2020. Dark and red lines represent fluxes of human movements for 2019 and 2020, respectively. (b) Relative movements from Wuhan to other provinces in China. (c) Timeline of the correlation between daily incidence in Wuhan and incidence in all other provinces, weighted by human mobility.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>Fig. 3. Human mobility explains early epidemic growth rate in China. (a) Daily counts of cases in China. (b) Time series of province-level growth rates of the COVID-19 epidemic in China. Estimates of the growth rate were obtained by performing a time-series analysis using mixed-effect model of lagged, log linear daily case counts in each province (Materials and Methods). Above the red line are positive growth rates and below are growth rates are negative rates. Blue indicates dates before the implementation of the cordon sanitaire and green after. (c) Relationship between growth rate and human mobility at different times of the epidemic. Blue indicates before the implementation of the cordon sanitaire and green after.</figcaption></figure>
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 4. Shifting age and sex distributions through time. (a) Age and sex distributions of confirmed cases with known travel history to Wuhan. (b) Age and sex distributions of confirmed cases that had no travel history. (c) Median age for cases reported early (before 1st Feb) and those reported later (between 1st – 10th February. Full distributions are shown in fig. S7. (d) Change through time in the sex ratio of (i) all reported cases in China with no reported travel history, (ii) cases reported in Beijing without travel history, and (iii) cases known to have travelled from Wuhan.</figcaption></figure>