key: cord-342268-azi9i2a8 authors: Zhao, Shi; Zhuang, Zian; Cao, Peihua; Ran, Jinjun; Gao, Daozhou; Lou, Yijun; Yang, Lin; Cai, Yongli; Wang, Weiming; He, Daihai; Wang, Maggie H title: Quantifying the association between domestic travel and the exportation of novel coronavirus (2019-nCoV) cases from Wuhan, China in 2020: a correlational analysis date: 2020-02-20 journal: J Travel Med DOI: 10.1093/jtm/taaa022 sha: doc_id: 342268 cord_uid: azi9i2a8 nan In the end of 2019, a novel coronavirus (2019-nCoV) emerged in Wuhan, China and was causing a serious outbreak of acute respiratory illness. 1 Wuhan locates in the centre of mainland China with a population of 14 million and is very conveniently connected to other parts of China through airlines and highspeed rails. 2 As of 31 January 2020 (5:00 p.m., GMT + 8), there were 9809 confirmed 2019-nCoV cases in mainland China, including 213 deaths and 180 discharges. 3, 4 Cases infected in Wuhan were also detected in many foreign countries or regions including Thailand, Japan, Republic of Korea, the United States, Canada and some European countries. 4 The World Health Organization (WHO) has declared that the novel coronavirus outbreak is a public health emergency of international concern. Official reports on the newly confirmed cases are released very rapidly (several times a day) after January 16, 3-5 as the official diagnosis protocol was released by the WHO on January 17. 6 Recent studies indicated the likelihood of travel-related risks of 2019-nCoV spreading both domestically and internationally. 7 , 8 Many major cities in mainland China reported the finding of imported cases, including Beijing, Shanghai, Guangzhou and Shenzhen. The outbreak is still on-going with an increasing trend in daily new cases. 3, 4 Before the Wuhan lockdown (official travel restriction) on January 23, virtually all cases found in other major cities were exported cases from Wuhan. Population flow data between major cities in mainland China are available online due to the rapid development of internet in recent decades, see https://qianxi.baidu.com/ (in Chinese). In this work, we quantified the association between the domestic travel load and the number of cases exported from Wuhan to other cityclusters in mainland China. Our city-clusters are the pool of top five cities in the top 10 provinces (in number of cumulative cases). Thus, we included 10 city-clusters in the analysis, and the details of the selected city-clusters can be found in Supplementary Data S1. We examined the association between the load of domestic passengers departed from Wuhan and the number of confirmed cases to the 10 city-clusters (including the three municipalities, Beijing, Shanghai and Chongqing). Data sets of the daily numbers of domestic passengers were obtained from the location-based services database of Baidu company from January 1 to 20. We selected the top 10 provincial regions (except Hubei) with the largest cumulative number of cases, which accounts for 68% of all cases reported outside Hubei, before the implementation of the city lockdown, on 23 January 2020. Cases from other provinces are scattered and showed no clear pattern. The daily numbers of passengers from Wuhan to the city-clusters of each province were adopted to measure the load of domestic passengers departed from Wuhan to that cluster. From now on, province means the city-clusters in that province. The daily number of cases time series are obtained via the online outbreak situation reports. 3 , 4 Daily cases for each cluster are scaled from daily province total. The detail is given in Supplementary Data. The association was formulated as follows: Here, c i,t represents the daily number of new cases in the ith provincial region on day t. The function E(·) is the expectation; the 'province i ' denotes the dummy variable for the ith provincial region in accounting for the heterogenicity among different provinces, and thus, α i is a locality-varying interception term. Term ξ i is the daily number of passengers from Wuhan to the ith province. The time index was denoted by t, and the term τ modelled the delay from exposure to being detected. The ε represented the daily number of new cases in Wuhan. Hence, the product term (ξ ·ε) was proportional to the rate of cases exported from Wuhan to other places, and this setting was consistent with the framework in Imai et al., 9 which estimated the outbreak size through exported 2019-nCoV cases overseas, as well as in more complex frameworks. 10 The β is the regression coefficient to quantify the association between the load of passengers multiplied by the local infectivity in Wuhan and the number of cases reported outside Wuhan. Hence, the change rate, denoted by ∆ = [exp(β × 100) − 1] × 100%, may be interpreted as the expected percentage change in the daily number of cases offsite found associated with per 100 increase in the daily number of passengers departed from Wuhan where there was one new case daily. We estimated β for three means of transportation, and the P-value less than 0.05 was considered as statistical significance. We also considered a baseline version of the model by replacing term 'β·ξ i,t−τ ·ε t−τ ' in Equation (1) with 'β 0 ·ξ i,t−τ ', which ignored the variation in the force of infection in Wuhan, thus a constant, as in Equation (2) log The likelihood ratio (LR) test was adopted to justify the model structure in Equation (1) against the baseline form in Equation (2) . For the delay term (τ ), it was expected to be equivalent to the incubation period, which was reported to be 5.2 days (95% CI: 4.1-7.0) in. 1 Finally, we examined the association with τ varying at 3, 4, 5, 6 and 7 days in sensitivity analysis. The LR test yielded statistically significant outcomes, which suggests that the model in Equation (1) is more reasonable than the baseline form in Equation (2) . Through the goodness-of-fit in terms of the McFadden's pseudo-R-squared, we found that τ = 5 days attains the best fitting performance to explain the patterns of the cases offsite detected. This well matched the estimate of the mean incubation period of the infection at 5.2 in. 1 We found a statistically significant positive association between the load of passengers multiplied by the local infectivity in Wuhan and the number of cases reported outside Wuhan, see Table 1 . We estimated that per 100 cases increase in the daily number of newly reported cases in Wuhan together with per 100 persons increase in the daily number of passengers departed Wuhan were likely to cause a 16.25% (95% CI: 14.86-17.66%) increase in the daily number of cases offsite detected on average. The sensitivity analysis by varying τ suggested this fundamental relationship holds, and the details can be found in Table 1 and the Supplementary Data S2. This analysis has limitations. The data used in this study are in the early phase of the outbreak and may not represent subsequent waves. We aim at quantifying the association between domestic travel and 2019-nCoV exportation in China. Although the differences in case ascertainment in different cities clusters are considered in the model by the dummy term 'province i ', the association could also bear heterogeneity. Temporal and spatial correlations were not addressed in this simple modelling analysis due to lack of data. The 2019-nCoV surveillance data were too scattered and short at this early stage to consider temporal and spatial correlation. The correlation between population flow, number of cases offsite detected, and source infection prevalence was addressed in this work. Our modelling framework would be extended to a more complex and realistic form for exploring the potential spatial correlations, and benefit from more detailed disease surveillance and travel population flow data. Supplementary data are available at JTM online. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia The association between domestic train transportation and novel coronavirus outbreak in China, from 2019 to 2020: a data-driven correlational report Situation report of the pneumonia cases caused by the novel coronavirus', released by the National Health Commission of the People's Republic of China Novel Coronavirus (2019-nCoV) situation reports, released by the World Health Organization (WHO) News press and situation reports of the pneumonia caused by novel coronavirus Laboratory testing for 2019 novel coronavirus (2019-nCoV) in suspected human cases: Interim guidance, released on January 17, 2020 by the World Health Organization (WHO) Pneumonia of unknown Etiology in Wuhan, China: potential for international spread via commercial air travel Nowcasting and forecasting the Wuhan 2019-nCoV outbreak Preprint published by the School of Public Health of the University of Hong Kong Estimating the potential total number of novel Coronavirus (2019-nCoV) cases in Wuhan City, China Preprint published by the Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study The authors would like to acknowledge colleagues for helpful comments. The authors declared no competing interests. All authors conceived the study, carried out the analysis, discussed the results, drafted the first manuscript, critically read and revised the manuscript and gave final approval for publication. The funding agencies had no role in the design and conduct of the study; collection, management, analysis and interpretation of the data; preparation, review or approval of the manuscript; or decision to submit the manuscript for publication. The ethical approval or individual consent was not applicable. All data and materials used in this work were publicly available and also available based on request. Not applicable.