TY - JOUR T1 - Changing transmission dynamics of COVID-19 in China: a nationwide population-based piecewise mathematical modelling study JF - medRxiv DO - 10.1101/2020.03.27.20045757 SP - 2020.03.27.20045757 AU - Jiawen Hou AU - Jie Hong AU - Boyun Ji AU - Bowen Dong AU - Yue Chen AU - Michael P. Ward AU - Wei Tu AU - Zhen Jin AU - Jian Hu AU - Qing Su AU - Wenge Wang AU - Zheng Zhao AU - Shuang Xiao AU - Jiaqi Huang AU - Wei Lin AU - Zhijie Zhang Y1 - 2020/01/01 UR - http://medrxiv.org/content/early/2020/03/30/2020.03.27.20045757.abstract N2 - Background The first case of COVID-19 atypical pneumonia was reported in Wuhan, China on December 1, 2019. Since then, at least 33 other countries have been affected and there is a possibility of a global outbreak. A tremendous amount of effort has been made to understand its transmission dynamics; however, the temporal and spatial transmission heterogeneity and changing epidemiology have been mostly ignored. The epidemic mechanism of COVID-19 remains largely unclear.Methods Epidemiological data on COVID-19 in China and daily population movement data from Wuhan to other cities were obtained and analyzed. To describe the transmission dynamics of COVID-19 at different spatio-temporal scales, we used a three-stage continuous-time Susceptible-Exposed-Infectious-Recovered (SEIR) meta-population model based on the characteristics and transmission dynamics of each stage: 1) local epidemic from December 1, 2019 to January 9, 2020; 2) long-distance spread due to the Spring Festival travel rush from January 10 to 22, 2020; and 3) intra-provincial transmission from January 23, 2020 when travel restrictions were imposed. Together with the basic reproduction number (R0) for mathematical modelling, we also considered the variation in infectivity and introduced the controlled reproduction number (Rc) by assuming that exposed individuals to be infectious; we then simulated the future spread of COVID across Wuhan and all the provinces in mainland China. In addition, we built a novel source tracing algorithm to infer the initial exposed number of individuals in Wuhan on January 10, 2020, to estimate the number of infections early during this epidemic.Findings The spatial patterns of disease spread were heterogeneous. The estimated controlled reproduction number (Rc) in the neighboring provinces of Hubei province were relatively large, and the nationwide reproduction number ‐ except for Hubei ‐ ranged from 0.98 to 2.74 with an average of 1.79 (95% CI 1.77‐1.80). Infectivity was significantly greater for exposed than infectious individuals, and exposed individuals were predicted to have become the major source of infection after January 23. For the epidemic process, most provinces reached their epidemic peak before February 10, 2020. It is expected that the maximum number of infections will be approached by the end of March. The final infectious size is estimated to be about 58,000 for Wuhan, 20,800 for the rest of Hubei province, and 17,000 for the other provinces in mainland China. Moreover, the estimated number of the exposed individuals is much greater than the officially reported number of infectious individuals in Wuhan on January 10, 2020.Interpretation The transmission dynamics of COVID-19 have been changing over time and were heterogeneous across regions. There was a substantial underestimation of the number of exposed individuals in Wuhan early in the epidemic, and the Spring Festival travel rush played an important role in enhancing and accelerating the spread of COVID-19. However, China’s unprecedented large-scale travel restrictions quickly reduced Rc. The next challenge for the control of COVID-19 will be the second great population movement brought by removing these travel restrictions.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThis work was supported by the National Natural Science Foundation of China from ZJZ(81673239,81973102), WL(11925103 & 61773125), the National Key R&D Program of China (Grant no. 2018YFC0116600), and the STCSM (Grant no. 18DZ1201000).Author DeclarationsAll relevant ethical guidelines have been followed; any necessary IRB and/or ethics committee approvals have been obtained and details of the IRB/oversight body are included in the manuscript.YesAll necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesThe daily cumulative numbers of confirmed and suspected COVID-19 cases for the period from December 1, 2019 to February 8, 2020 were obtained from the National, Hubei Provincial and Wuhan Municipal Health and Family Planning Commissions. The daily population movement data into and out of Wuhan by types of transportation (i.e., long-distance bus, train and air) for the period January 10-23, 2020 were retrieved from an app-derived real-time location database provided by Baidu Huiyan Company. ER -