key: cord-0986660-dounwt47 authors: Zhang, Yahua; Zhang, Anming; Wang, Jiaoe title: Exploring the roles of high-speed train, air and coach services in the spread of COVID-19 in China date: 2020-05-26 journal: Transp Policy (Oxf) DOI: 10.1016/j.tranpol.2020.05.012 sha: 514762bd6ecd023c793b266919cfa6728ca4f066 doc_id: 986660 cord_uid: dounwt47 To understand the roles of different transport modes in the spread of COVID-19 pandemic across Chinese cities, this paper looks at the factors influencing the number of imported cases from Wuhan and the spread speed and pattern of the pandemic. We find that frequencies of air flights and high-speed train (HST) services out of Wuhan are significantly associated with the number of COVID-19 cases in the destination cities. The presence of an airport or HST station at a city is significantly related to the speed of the pandemic spread, but its link with the total number of confirmed cases is weak. The farther the distance from Wuhan, the lower number of cases in a city and the slower the dissemination of the pandemic. The longitude and latitude coordinates do not have a significant relationship with the number of total cases but can increase the speed of the COVID-19 spread. Specifically, cities in the higher longitudinal region tended to record a COVID-19 case earlier than their counterparties in the west. Cities in the north were more likely to report the first case later than those in the south. The pandemic may emerge in large cities earlier than in small cities as GDP is a factor positively associated with the spread speed. Patients with pneumonia of unknown causes were first detected in Wuhan, the capital city of Year migration. To halt the spread of the virus, an emergent city lockdown measure was imposed in Wuhan, the epicentre of COVID-19, on 23 January. The lockdown was later extended to the whole province. The shutdown affected a population of about 60 million, including 11 million in Wuhan. Although this radical move was questioned by many people due to the high economic (and other) costs associated, this approach appeared to be vindicated as the curve for the spread of the disease outside Hubei flattened in mid-February 2020, and the same pattern emerged in Wuhan in late February 2020. On 14 March 2020, only Wuhan was classified as a high-risk area and the rest of the province was medium-or low-risk. In the next few days, the lockdown gradually loosened throughout Hubei Province. On 8 April 2020, the city of Wuhan reopened after more than 10 weeks in lockdown. Air and train services have resumed and the city's recovery is being closely watched by the world. Wuhan is one of China's most important transport hubs, and is therefore very well connected to other areas in China (Zhu et al., 2019) . The present research aims to examine the role that each transport mode played in diffusing the cases of COVID-19. Various studies have considered the link between the movement of people and the imported cases of COVID-19. Zhao et al. (2020) found that the number of air passengers from Wuhan and local population can be used to explain the number of cases in the infected cities. Ai et al. (2020) reported a significant and positive relationship between population movement and the number of COVID-19 cases. They argued that some cases could be avoided and prevented if the city closure was implemented earlier. However, most of these studies only consider one transport mode or the total movement of people regardless of the transport means. More specifically, this paper investigates the factors influencing the number of imported cases from Wuhan and the spread speed and pattern of the pandemic. The gravity model is used with a consideration of the factors of the frequencies of high-speed train (HST), coach and air services (flights) between Wuhan and the other domestic cities, the cities' GDP and distance from Wuhan as well as longitude and latitude, and the presence of an airport/HST Station. HST, coach (inter-city bus) and air services were the most important means of transporting the five million people out of Wuhan before the lockdown. We believe that exploring the relationship between transport means and the speed and spread pattern of infectious diseases such as COVID-19 is important to both the general public, policy makers and health professionals. Such relationship can also be used for future prediction and control or mitigation when similar incidences take place. We find significant links between flight and HST frequencies out of Wuhan and the number of COVID-19 cases in the destination cities. The presence of an airport or HST station at a city is significantly related to the speed of the pandemic spread, but its link with the number of total confirmed cases is weak. The farther the distance from Wuhan, the lower number of cases in a city and the slower the dissemination of the pandemic. The longitude and latitude values do not have a significant relationship with the number of total cases, but can increase the speed of the COVID-19 spread. Cities in the south were more likely to report the first case earlier than those in the north. The pandemic may emerge in large cities earlier than in small cities. The paper is organised as follows. Section 2 describes background and reviews the relevant studies, and Section 3 discusses our empirical methodology. The results are reported and discussed in Section 4. Finally, Section 5 contains the concluding remarks. Wuhan is one of China's most important transport hubs. It has three railway stations, one international airport, and a significant inland port that lies at the confluence of the Yangtze and Hanjiang rivers. China's first long-haul HST line, the Wuhan-Guangzhou HST, was launched in 2009 . This HST line was later extended to Beijing and Hong Kong, forming the longest north-south trunk line of China's "4+4" HST network (four northsouth and four east-west trunk lines; see, e.g., Wang et al., 2017) . Wuhan is also the midpoint of the Shanghai-Wuhan-Chengdu HST line, which is one of the four east-west trunk lines. Therefore, China's most important cities such as Beijing, Shanghai, Guangzhou, Chengdu, Shenzhen and Hong Kong can be reached within five hours from Wuhan by HST. Inside Hubei Province, the Wuhan-Shiyan HST commenced services at the end of December 2019, connecting Wuhan to Hubei's key industrial cities such as Xiangyang, and Shiyan. With all these connections, Wuhan has been called the heart of China's HST network and it is ranked the fourth place in China by the weighted degree and closeness centrality, and the second place by weighted betweenness centrality (Jiao et al., 2017) . It has been ranked the fourth most important rail-air combined node in China's city network (after Shanghai, Beijing and Guangzhou; see Zhu et al., 2018) , and fourth in terms of the total rail passengers handled and first in terms of the number of transfer passengers in recent years. 1 Furthermore, Wuhan Tianhe International Airport is the largest airport in Central China, ranking the top-15 airports in China by all three centrality indexes (Wang et al., 2011) . It handled 27 million passengers in 2019, an increase of 10.8% from 2018. Besides, Wuhan plays an important role in China's intercity coach network . The weekly frequencies of flights, HST and coaches out of Wuhan are illustrated in Figures 1, 2 and 3, respectively. 1 Wuhan is also ranked among the top in criticality and importance. For example, its closure may reduce the performance of HST network by 10.3% ; and it is ranked the second in rail-air network by criticality (Li et al., 2019) . provinces were Henan, Hunan, Anhui, Chongqing, and Jiangxi. It is expected that the population that flowed out of Wuhan may have significantly impacted on the outbreaks in other cities. In fact, since the outbreak of COVID-19, tracking passengers from Wuhan and putting them into quarantine were the top priority of almost every city from late January to February 2020, as most of the confirmed cases were imported cases or close contacts of the imported cases. A few centuries ago, infectious diseases could spread only as fast and as far as people could walk, horses could gallop and ships could sail (Tatem et al., 2006) . The emergence of modern transport modes has accelerated the spread of new diseases in an unprecedented speed and put more people at risk. In this century, human beings have been hit by a series epidemic diseases including SARS, H1N1 pandemic (in 2009), MERS, and most recently, COVID-19 at a global level, owing largely to the growing availability of affordable air travel. The risk of transmission of respiratory infections on airplanes was a major concern to the public and airline industry. Some studies suggested that air transport is a mode that contributes to the accelerating and amplifying influenza propagation, as such transmission can occur in the flight or at the airports (see a comprehensive survey by Browne et al., 2016) . However, other studies on the ventilation systems and patient outcomes showed that the dissemination of pathogens during the flight occurs rarely (see a good survey by Leder and Newman, 2005) . For example, quite a few studies suggested that in-flight transmission of SARS was not common. Wilder-Smith et al. (2003) reported that in-flight transmission occurred in one of the three flights with SARS patients on board. The authors note that the risk of in-flight transmission is lower than that reported for influenza, but may increase if super-spreaders are on board. However, while the transmission risk in the flight may be low, flights can carry people with the virus to new places. Bowen Jr and Laroe (2006) attempted to examine the link between air transport accessibility and the speed of SARS diffusion. The accessibility to China was measured by flight frequency and directness of scheduled airlines services. They found that airline network accessibility was a determinant of the speed with which SARS arrived in infected countries, but its influence diminished in later weeks of the outbreak due in part to the public health measures such as health screening at the airports. Research into the role of ground transport in the propagation of infectious diseases is limited. Among the few studies, Troko et al. (2011) found that the use of public buses or trams is a significant factor for the acquisition of acute respiratory infection. The risk can be greater among occasional bus or tram users. Therefore, it is important to exercise good respiratory hygiene and refrain from making unnecessary travels by public transport during the pandemic outbreak season. Cui (20011) Using the laboratory-confirmed cases of H1N1 for the period from May 2009 to April 2010, Cai (2019) evaluated the effects of airports and railway stations on arrival days and peak days of the virus. The authors found that airports and railways stations in Chinese prefectures significantly advanced the arrival days, but they were not significantly associated with the peak days of the pandemic. Du et al. (2020) is one of the few studies that estimated the probability of transportation of COVID-19 from Wuhan to 369 other cities in China before the quarantine. They expected that in 130 cities, the COVID-19 risk was greater than 50% and this risk was greater than 99% in the four largest metropolitan areas. However, as with many other studies on COVID-19, they only attributed the spread of COVID-19 to the travel or movement of people without looking at how different transport modes have specifically contributed to the dissemination of the virus. The present research aims to address this issue. The gravity model is used to identify the factors, particularly different transport modes, associated with the number of COVID-19 cases in Chinese cities. Our first dependent variable is the cumulative confirmed cases of COVID-19 reported on 15 February 2020, most of which were imported cases with a travel history to Wuhan. 2 This date was chosen because of the following several reasons. First, there was a shortage of testing kits for COVID-19 before mid-February and patients needed to meet strict criteria to become eligible for a test. Second, patients have an incubation period before various symptoms develop. The WHO estimates that the incubation period ranges from 1 to 14 days, with a median of 5 to 6 days. 3 China reported a case with 24 days of incubation. In addition, the quality of the testing kits from the early days may not have been satisfactory, as they reported a high proportion of false-negative results. The definition of a confirmed case set by China's National Health Commission changed over time. From January to March 2020, seven versions of the Novel Coronavirus Pneumonia Diagnosis and Treatment Plan were issued. The fifth edition even recommended that Hubei Province use CT scan results as confirmation of the infections in order to have a quick identification in the presence of high false-negative risk. Although by 15 February, there had been community-transmission cases, the percentage was relatively low and most of these cases were associated with close contact with those who had travel history to Wuhan. It would not distort the results significantly by including them in the imported cases. As a robustness check, we also use the number of cases on 1 February 2020 as the dependent variable, but it should be noted that the number of cases on this date may not accurately reflect the true cases for the abovementioned reasons. The imported cases are analogous to the trade flows in the gravity model, which has been widely used in cross-country empirical analyses of international trade flows. The gravity model has also been applied to air transport to identify the determinants of bilateral air passenger or air cargo flows (Zhang and Lu, 2013; Zhang and Zhang, 2016) . It has been reported that a large portion of the international trade matrix consists of zero trade, either because data are missing or some economies simply do not trade (Helpman et al., 2008) . This is also the case in our study where about one tenth of the cities report zero cases during our study period. The problem of zero or missing trade flows when using gravity models has long been ignored in the international trade literature . Many researchers simply discard the zero flows from the sample, which results in the loss of information. Some add a small constant value to the zero values to allow the estimation of log-linear equations (van Bergeijk and Brakman, 2010) , but this method does not have a theoretical basis. A more acceptable approach was introduced in Santos Silva and Tenreyro (2006), where Poisson pseudo-maximum likelihood estimator (PPML) technique is proposed. This approach can deal with the zero problem in the dependent variable and is also capable of coping with the heteroscedasticity problem that is common in the gravity data, in which case, the parameters of log-linearised gravity models estimated by ordinary least squares (OLS) will be highly misleading. Other nonlinear transformations such as the nonlinear least squares method and Tobit regressions are also not working if the errors are heteroscedastic. Therefore, a gravity type model with the PPML estimation approach is used in this study. The empirical gravity model employed is expressed as follows: The variables are explained as follows: • CASE215 is the number of cumulative confirmed COVID-19 cases on 15 February 2020 in each city. We will also replace it with the data of 1 February 2020 (CASE201). The data are from China National Health Commission's daily report. • GDP is the gross domestic product of each city. The data are obtained from China City Statistics Book 2018 and the 2017 data are used. • AIR denotes weekly flight frequency from Wuhan to each city in 2019. The data source is Official Aviation Guide (OAG). Bowen Jr. and Laroe (2006) argue that it is important to establish the relationship between the diffusion of the disease and the schedule data as in a potential pandemic crisis, the schedule data are easily obtained while the actual passenger data are only available much later. 4 • COACH denotes weekly coach (bus) services from Wuhan to each city. The information is extracted from Xinxin Travel (www.cncn.com) using 2018 data. • HST denotes weekly HST frequency from Wuhan to each city in late 2019. The data are taken from China Railway's booking website (www.12306.cn). • HUB is a dummy taking the value of one if a city has an airport or HST station, and zero otherwise. • lnLAT is the latitude of each city's administration centre in logarithm. lnLONG is the longitude of each city's administration centre in logarithm. These two variables are included because there have been reports suggesting that the virus is most active in a temperature range between 5 to 11 degrees Celsius (Sajadi et al., 2020) . Sajadi et al. (2020) We are also interested in the link between the transport modes and the emergence of the first case of COVID-19 in each city, which is to be referred to as the speed of the transmission. Before we run the regressions, we conduct a simple t-test to see whether cities with and without a particular transport hub showed a significant difference in the number of cases using the 15 February data. A city is called a transport hub in this study if it has an airport or an HST station. The Levene test is used to check the equality of variance of the data. The Levene test is not significant at the 5% level but significant at the 10% level. Therefore, the homogeneity of variance assumption can be assumed at 5%. We find that cities without airports or HST stations recorded an average of 30.7 cases, while those with airports or HST stations reported an average of 87.7. The difference is not statistically significant at a significance level of 5% if equal variances are assumed. Interestingly, if unequal variances are assumed, the difference is statistically significant. This is also the case when the 1 February 2020 data are used. Therefore, it seems that the link between the presence of an airport or HST station and the number of cumulative confirmed cases is not very strong based on the t-test. We should therefore resort to the multiple regression results and include more transport-related variables. For the arrival day (ARRDAY) variable, we find that cities with airports or HST stations recorded its first case of COVID-19 within a significantly shorter time (17.7 days) calculated from 10 January 2020, than those without such facilities (24.8 days) at the level of 1%. It seems that a transport infrastructure such as airport or HST station can speed up the spread of the virus. Specifications (3) and (4) include the HUB dummy while specifications (1) and (2) do not. The distance variable is included in specifications (1), (3) and (5) while in specifications (2) and (4), this variable is replaced by travel time (TIME). As can be seen, all the specifications have the expected signs, and the levels of significance are largely consistent. However, the presence of an airport or an HST station does not have a significant impact on the number of cases after other variables are controlled for. It is understood that the effect of the HUB dummy can be diluted as it might be correlated with other variables. For example, cities with airport or HST stations tend to be large or at least medium-sized cities, and thus have higher GDP. The correlation matrix of all the independent variables suggest that moderate correlations exist between some variables, but none of them are greater than 0.6. The problem of multicollinearity may arise, but the consequence is not serious particularly because one of the purposes of building the gravity is for future prediction purpose. In fact, multicollinearity does not result in biased estimations. Achen (1982) pointed out that the only effect of multicollinearity is that it makes it hard to get coefficient estimates with small standard error and thus can lead to insignificance of coefficients. Flight frequency is positively and significantly associated with the number of confirmed cases at the level of 1%. If the weekly frequency increases by one flight, we would expect that the number of the cases increases by 1% according to specification (3) in Table 2 . The frequency of HST also has an impact on the number of cases at the 10% level for specifications (1) to (4). An increase in HST frequency by one unit is associated with a much smaller increase in the number of cases compared with air travel. Coach services seem not to have any significant links with the imported cases, which is a bit surprising. One possible reason is that the vast majority of coach services are for short trips. Unlike air travel that requires passengers to report to the airport at least one hour before the departure, and go through several formalities before boarding, Coach passengers normally arrive at the station a short time before the departure and during the trip, they tend to stay in their seats and do not move around as they would do on the HST train. This implies that the risk of coach passengers' exposure to the COVID-19 virus is relatively low. The distance and travel time are negatively and significantly related to the confirmed cases, which is consistent with the COVID-19 distribution pattern. Considering that Wuhan is The distance variable is positively related to the arrival day variable and statistically significant, implying that the farther the distance from Wuhan, the later the emergence of the confirmed COVID-19 case. GDP has a significantly negative impact on the speed of the spread and its impact is larger for higher quantiles. The air transport frequency could have facilitated the spread of COVID-19 at lower quartiles, but it is not significantly associated with the arrival day variable at the 0.5 and 0.75 quantiles. The frequencies of HST and coach services are not statistically significant. The HUB variable, or the presence of an airport or an HST station can increase the speed of the COVID-19 transmission as this variable is consistently negative and significant at medium and higher quantiles. This is consistent with Cai et al. (2019) who reported that the presence of airports or high-ranking railway stations in Chinese prefectures significantly advanced arrival days of the H1N1 pandemic. Note that the interpretation of these coefficients in quantile regressions is the same as that for the OLS estimates. For example, at the 0.25 quantile, every one unit increase in flight frequency will lead to a reduction in the arrival days by 0.013. The effects of the variables of GDP, flight frequency, distance, hub, latitude, and longitude coordinates can be more clearly seen from the graphs in Figure 5 . The graphs show the change in quantile coefficients along with the confidence interval of 95% (the shaded area). The quantile process plots allow us to readily identify which predictors are associated with different parts of the response distribution. The black dashed lines are the OLS estimates and their confidence intervals. The coefficient of the log-transformed latitude does not vary much and is no different from its OLS estimate for most of the quantiles less than 0.8. The quantile process plot for the log-transformed longitude variable is quite interesting. The coefficient decreases first but starts to increase after the 0.8 quantile. Its coefficient in the lower quantiles (less than 0.5) well exceeds the OLS estimates. In the flight frequency panel, we can see that the coefficient does not differ much from the OLS estimates and is relatively stable across quantiles. The positive effects of the distance become stronger as the quantile increases. In contrast, the negative effect of the transport hub gradually becomes significant and stronger as the quantile rises. To understand the roles of different transport modes in the spread of COVID-19 pandemic imported cases from other countries became a major concern. The government therefore decided to limit the number of international flights to China. Chinese airlines were only allowed to operate one route to any specific country with one flight each week. Each foreign airline was only allowed to operate one route to China with one flight per week. With these measures, China has quickly brought the imported cases down to a manageable level. We also found that the magnitude of the HST frequency coefficient is much smaller than that of the flight frequency in the gravity model. This may have implications to the measures used to contain the pandemic. Many countries have adopted a combination of containment and mitigation measures to delay the surges of patients and flatten the curve to avoid overwhelming the health care system. One extreme measure is to lock down a city and cut off all the transport links with other cities. This normally involves the shutdown and halt of most businesses and outdoor activities, which comes at a huge economic cost. Another extreme is to achieve herd immunity. This option allows the infections to rise in a controlled way until 60% of the population gets infected, after which it would be harder for the virus to spread (Chang, 2020) . Between the two methods, a discretionary control policy can be used. For the transport links, these measures can be considered: reduce the air transport frequency without cancelling all the flights; cut air services first while keeping most of the HST services, considering the lesser impact of HST in spreading the disease; tighten the restrictions on residents' movements when the risk is high and loosen a bit when the pandemic is under control. This research reports some interesting spatial distribution patterns of the COVID-19 cases. The farther the distance from Wuhan, the lower number of cases in a city and the slower speed for the pandemic to be disseminated. The longitude and latitude coordinates do not have an impact on the number of cases, but are significantly associated with the speed of the COVID-19 spread. Specifically, cities in the east tended to record a COVID-19 case earlier than their west counterparts. Cities in South China may detect the first case earlier than those in North China. The pandemic might spread in large cities first before they arrive to small cities, as GDP is a factor positively associated with the spread speed. All these results suggest that the outbreak of a pandemic in a large city lying in the centre of a country such as Wuhan can have devastating impact on the whole nation and that the extreme measure to lock down the city might have been a correct move to reduce such impact on other parts of China. This is particularly important to Hubei's neighbouring provinces and the well-developed areas in East and South China such as Shanghai, and Guangdong. Interpreting and using regression Population movement, city closure and spatial transmission of the 2019-nCoV infection in China Airline networks and the international diffusion of severe acute respiratory syndrome (SARS) The roles of transportation and transportation hubs in the propagation of influenza and coronaviruses: a systematic review The hidden geometry of complex, Network-Driven Contagion Phenomena Roles of Different Transport Modes in the Spatial Spread of the 2009 Influenza A (H1N1) Pandemic in Mainland China The decision Australia needs to make now on how we end the coronavirus epidemic Risk for Transportation of Coronavirus Disease from Wuhan to Other Cities in China Estimating trade flows: Trading partners and trading volumes. The quarterly journal of economics Impacts of high-speed rail lines on the city network in China Quantile Regression Regression quantiles Respiratory infections during air travel A comprehensive method for the robustness assessment of highspeed rail network with operation data: A case in China Vulnerability analysis and critical area identification of public transport system: A case of high-speed rail and air transport coupling system in China How Wenzhou, 900 km from Wuhan, went into total lockdown. Sixth Tone Temperature and latitude analysis to predict potential spread and seasonality for COVID-19. Available at SSRN 3550308 Global transport networks and infectious disease spread Five million people left Wuhan: An escape or normal travel? Is public transport a risk factor for acute respiratory infection The Gravity Model in International Trade: Advances and Applications Should China further expand its high-speed rail network? Consider the low-cost carrier factor Exploring the network structure and nodal centrality of China's air transport network-A complex network approach Inter-city connections in China: High-speed train vs. inter-city coach Low risk of transmissio n of severe acute respiratory syndrome on airplanes: the Singapore experience Impacts of high-speed rail on airlines, airports and regional economies: A survey of recent research Gravity models in air transport research: A survey and an application Low cost carriers in China and its contribution to passenger traffic flow Determinants of air passenger flows in China and gravity model: deregulation, LCCs, and high-speed rail Tracking the spread of novel coronavirus (2019-nCoV) based on big data Connectivity of intercity passenger transportation in China: A multi-modal and network approach Measuring multi-modal connections and connectivity radiations of transport infrastructure in China The authors wish to thank Haoran Yang and Delin Du for collecting some of the research data, and Tao Li and Jessica Zhang for helpful comments. This work is financially supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant XDA19040402).