{'fraction_non_alphanumeric': 0.02408676741817381, 'fraction_numerical': 0.013406006551600737, 'mean_word_length': 5.1075991930060525, 'pattern_counts': {'":': 0, '<': 9, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
COVID-19 Pandemic with Human Mobility Across Countries

This study develops a holistic view of the novel coronavirus (COVID-19) spread worldwide through a spatial–temporal model with network dynamics. By using a unique human mobility dataset containing 547 166 flights with a total capacity of 101 455 913 passengers from January 22 to April 24, 2020, we analyze the epidemic correlations across 22 countries in six continents and particularly the changes in such correlations before and after implementing the international travel restriction policies targeting different countries. Results show that policymakers should move away from the previous practices that focus only on restricting hotspot areas with high infection rates. Instead, they should develop a new holistic view of global human mobility to impose the international movement restriction. The study further highlights potential correlations between international human mobility and focal countries’ epidemic situations in the global network of COVID-19 pandemic. Electronic supplementary material The online version of this article (10.1007/s40305-020-00317-6) contains supplementary material, which is available to authorized users.


Introduction
Although countries worldwide have gradually restricted nonessential international travel, particularly those from high-infection areas, it does not seem to have had the desired effect of stopping the spread of the novel coronavirus  in the past few weeks. The epidemic continues to be severe in the European and North American regions. It has even gradually spread to more areas, such as Africa and South America, leading to controversy over the pathway of the epidemic transmission [1,2]. The urgent situation needs a better assessment of the COVID-19 spread under global human mobility [3,4].
Besides non-pharmaceutical interventions within each country [5], restrictions on nonessential international travel from epidemic areas are proposed as a critical strategy to slow the spread [6]. It intends to cut off people's outbound movement from a country when it becomes a hot spot of the epidemic. This approach may not be as effective as policymakers expect if they ignore the current reality of global human mobility. Another important decision-making factor that has been overlooked is how to impose internal movement control policies, such as curfew or other forms of travel restriction, in the infected areas. These issues are particularly important in designing a global strategy to respond to the COVID-19 dynamics and to recover social and economic activities. Figure 1 summarizes three different views on understanding the role of human mobility in the pandemic. The first view is a focal view that addresses the role of internal mobility and movement controls in tackling epidemic in a single country while ignoring the influence of inbound and outbound mobility [7][8][9][10]. The second one is a dyadic view that explores the role of international mobility and travel restriction between a centric country (usually a hotspot area in the epidemic) and some other countries [11][12][13]. This approach usually does not take into account the impact of population movements other than those of the centric country, and in addition, only the one-way impact of the centric country on other countries. The third one is a holistic view that simultaneously analyzes the multiple paths of international mobility and travel restrictions in a global network and thus allows for the correlation of each path in the network. In this way, the networked approach allows the study of simultaneous and asymmetric effects between multiple countries through different mobility pathways.
Global network-based analysis is now more urgently needed than the previous focal view approach [7][8][9][10] or the dyadic view approach [11][12][13], that is, investigating epidemic developments between countries in the global network where each country is a node and the inter-country human mobility between two countries is represented by an edge. In this way, scholars and policymakers can gain the holistic insights on the virus' spread across the globe and develop effective and worldwide-coordinated measures to contain the pandemic [14].
This study responds to this urgent call by developing a spatial-temporal model with network dynamics [15,16] to understand the correlation between the COVID-19 epidemics in different countries accounting for the inter-country human mobility and international travel restrictions targeting different countries. Our analysis of the COVID-19 pandemic starts from January 22, 2020, when the epidemic was officially reported by the World Health Organization (WHO), to April 24, 2020, when the number Three views on the role of human mobility in the pandemic. Note: The solid and dashed lines indicate unrestricted human mobility before public interventions and restricted mobility after public interventions, respectively. A (solid or dotted) line with the arrow from country A to country B stands for the correlation between the epidemics in these two countries accounting for (restricted or unrestricted) international mobility. An arrow connector on a country itself represents the correlation between the internal spread and (controlled or unrestricted) internal movement of COVID-19 infections exceeded 2.6 million worldwide. Based on the timing of the outbreak and the scale in terms of the accumulated confirmed infection numbers in the COVID-19 pandemic, we selected 22 countries from six continents, which accounted for 86.70% of total infection amount worldwide as of April 24, 2020. Detailed information of the countries is provided in Table S1 in supplementary materials. We collected the daily number of confirmed COVID-19 infections in every country from the WHO. We also use a unique dataset containing 547 166 flights with a total capacity of 101 455 913 passengers among the 22 countries from January 22 to April 24, 2020. In the follow-up analyses, the information is aggregated to the countryday level, resulting in 23 885 daily pairs of international travel movement across six continents. Figure S1 illustrates a clear decreasing trend since early March regarding the capacity and number of international flights among these 22 countries. We also compile data from multiple sources containing international travel restriction and internal movement control policies (such as curfew and other forms of domestic travel restriction), as summarized in Table S1.
By treating each country as a node in the global network and the international human mobility as the connections between countries, we model the dynamic process of the epidemic. The unit of analysis in this study is the number of confirmed infections in country i on day t. At the beginning of the epidemic, every country had zero infection and human mobility was unrestricted between countries. Likewise, there were neither non-pharmaceutical interventions within countries nor international mobility restriction between countries. Therefore, countries around the world had formed a fully connected human mobility network.
During the pandemic period, however, travel restrictions were introduced by countries gradually and applied to travelers from different countries at different time points, which led to a dynamically changing global human mobility network. The restriction reduces the human mobility between countries. If the action, saying the travel restriction from a hotspot country A to another country B, was introduced in time and did take effect, we would expect to see a weaker association between the epidemic in country A and the daily new cases in country B in the controlled international travel period than in the pre-controlled period, which is equivalent to a positive effect of such correlation before the introduction of international travel restriction with reference in the restricted period. This would effectively imply a decreasing association of epidemic situations from country A to country B with the introduction of international travel restriction, largely due to the fewer number of travelers from country A to country B. However, if the international travel restriction is not effectively implemented by country B regarding timing and targeting countries, we might even observe a stronger association between the daily new cases in country B and the epidemic in country A after country B bans the entry of international travelers from country A, which is equivalent to a negative association between the daily increment in country B and the epidemic in country A in the pre-controlled period for international travel with reference in the controlled period.
In addition, reducing the spread of the epidemic does not solely depend on the international travel restriction and would also rely on the effective internal movement control within each country. Specifically, even if country B reduces international movement from the hotspot country A, the virus may still have the chance to spread within country B (e.g., via internal spread or from the third country), so that the impact of the epidemic in country B cannot be reduced without effective internal movement control policies. Therefore, if country B effectively implements internal movement control, we would also expect the control to work with international mobility restriction and observe a decreasing epidemic trend in country B, after a certain lagged period.
This idea is applied to all the countries in the network, taking into account the multiple paths of international mobility and travel restrictions between any two of them simultaneously. Through such a global network perspective, we aim to examine the correlations of COVID-19 epidemics in respective countries and the extent to which the human mobility and the introduction of international travel restriction targeting hotspot countries influence the pandemic.

Data
We collected multi-sourced datasets in this study. The daily numbers of confirmed infections of all countries were collected from the WHO COVID-19 Dashboard. 1 Based on the timing of the outbreak occurred and the scale in terms of their accumulated confirmed infection numbers during the outbreak, we selected 22 countries from six continents: 6 countries from Asia (China, Iran, India, Japan, South Korea, and Turkey), 10 countries from Europa (Belgium, France, Germany, Italy, Netherlands, Portugal, Russian Federation, UK, Spain, and Switzerland), 2 countries from North America (Canada and USA), 2 countries from South America (Brazil, Ecuador), 1 country from Africa (South Africa), and 1 country from Oceania (Australia). These 22 countries accounted for 86.70% of the total infection amount worldwide as of April 24, 2020. In addition to the WHO data, we also collected data from and the COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University [17] to replicate the analysis.
We further use a unique human mobility dataset containing daily global commercial flights among the selected 22 countries between January 22, 2020, and April 24, 2020, from a leading data consulting company in the civil aviation industry, VariFlight. 2 The dataset contains information on the origin country, destination country, date, the number of flights, and total capacity in terms of the maximum available seats. In summary, this dataset provides information on 547 166 flights between 23 855 pairs of origin and destination countries that covers a total capacity of 101 455 913 passengers during the period.
Moreover, we collected data on the international travel restriction between countries and the internal movement control within every country from Oxford COVID-19 Government Response Tracker 3 and GardaWorld Crisis24 Global Portal, 4 which timely document epidemic prevention policies in countries around the world. In a global network, the path connection from country A to country B is set by value 1 at day 0, representing the availability of international movement from A to B. Such connection is disabled on day t, due to either the internal movement control in country A or the entry ban for travelers from country A into B. Thus, on day t onward, the connection from A to B is set to be 0 to represent such disconnection. Table S1 summarizes the timing and scale of the outbreak, as well as travel restriction policies by the 22 countries.

Dynamic Spatial Network Model
The unit of our analysis is the number of confirmed infections in country i on day t. Given the spatial nature of this research, we develop a dynamic network model based on the spatial-temporal features [15,16] to examine the extent to which the number of newly confirmed infections in each country is correlated with the cumulative number of infections in each of the remaining 21 countries, by taking into account (1) the mobility volume between countries, (2) the introduction of inter-country travel restriction policies targeting different countries at different time point, and (3) internal movement restriction within each country. Specifically, the number of newly confirmed infections in the country i (i 1 to 22) on day t (t 1 to 94, starting at January 22 till April 24), NCI i,t , is given by where ε it is the error term, the coefficient α 1i measures the τ -day lagged effect of new infections in the same country (NCI i,t−τ ), and α 2i captures the ω-day lagged effect related to internal movement restriction in the same country (PreIR i,t−ω ), such as social distancing and home-stay orders, where Specifically, the significantly positive value of α 2i represents the positive association between fully allowed internal movement (with reference to the restriction) and the new infections in the same country i after ω days. This effectively means that there was the less severe epidemic in the controlled period (PreIR i,t−ω 0) than in the pre-controlled period for internal movement (PreIR i,t−ω 1) and thus implies the effectiveness of internal movement control policy in country i. On contrary, the significantly negative value of α 2i would suggest the more severe epidemic in the controlled period (PreIR i,t−ω 0) than in the pre-controlled period (PreIR i,t−ω 1), while the insignificant value of α 2i means there is no statistical difference on the epidemic in country i before and after introducing the internal movement control policy.
Furthermore, C I jt is the number of cumulative confirmed infections in each of the rest of countries other than country i by day t ( j i) and FCap i jt stands for the human mobility from country j to i, such as the number of direct flights or maximum flight seats in the direct flights from country j to i on date t, which might be further restricted by the international travel control policies introduced by country i targeting different countries on different days. To capture such dynamic international travel controls, W t is the time-dependent spatial weight matrix of the international mobility network between all considered countries on day t, with its element w i jt defined as Specifically, the travel from country j to i on day t might be restricted due to the respective policies in the origin or destination countries. In the origin country j, its internal movement restriction on day t might make regular travelers difficult or even impossible to leave country j. The destination country i may also ban the entry of travelers from country j. Therefore, α 3i j (i.e., the coefficient of the term FCap i j,t−θ W t−θ C I j,t−θ ) measures the θ -day lagged correlation between the epidemic in country j and the daily new cases in country i, accounting for international mobility capacity and the international travel restrictions, where the θ -day lag can be largely attributed to the incubation period. The significantly positive value of α 3i j means that the correlation between the cumulative epidemic in country j and the θday lagged new infections in country i is stronger in the pre-controlled period for the international mobility (w i jt 1) than in the international mobility controlled period (w i jt 0). In other words, the positive value of α 3i j suggests the weaker correlation effect in the international mobility controlled period and thus implies the effectiveness of introducing international mobility control policies on containing the COVID-19 epidemic in country i. Conversely, the significantly negative effects of α 3i j means that the stronger correlation between the cumulative epidemic in country j and the θ -day lag new infections in country i in the international mobility controlled period (w i jt 0) than in the pre-controlled period (w i jt 0), which might be due to the delayed implementation of international travel restriction or the myopia of only banning international travelers from the existing hotspot but ignoring the other countries with emerging epidemics. In addition, the insignificant value of α 3i j means the intercountry correlation of epidemics remains largely unchanged with the introduction of international travel restrictions from country j to i.

Estimating the Model
The econometric model we develop above has a system of 22 simultaneous equations with the time-dependent spatial weight matrix W t . This unique feature makes our model differs from the conventional spatial models, including the dynamic spatial panel model, which typically use the time-invariant spatial weight matrix based on the geographical characteristics [16]. Therefore, we follow the process below to estimate the model.

2.
Organize the data into a time series format, with each row representing the date and each column for the different data variables on that date, such as the coded elements of the spatial weight matrix, the number of infection cases in every country and the numbers of flights and maximum available seats in every pair of international mobility route. 3. Create time-lagged variables where necessary. In the empirical analysis, we take τ 1 to one-day lagged effect of new infections, ω 14 to check the effect of internal movement restriction within each country, and θ 14 to capture the influence of 14-day incubation period of the COVID-19 on its international spread. 4. Construct an interdependent system of linear regression equations, and use seemingly unrelated regression (SUR) [18][19][20] to estimate the parameters as specified in Eq. (1). The SUR is also known as joint generalized least squares (JGLS), and it is a generalization of OLS for multi-equation systems by allowing for the correlation of the error term of every equation. In an m-equation system for T-period observations, The variance-covariance matrix for the error term vector e can be written as where ⊗ is the Kronecker multiplication operator, Σ is an m* m positive definite symmetric matrix with σ 11 as the variance of the error in the first equation and σ i j as the covariance between the errors of the ith and jth equations, and I is the identity matrix. The estimation process of SUR is as follows: a. First apply OLS to every equation and obtain the residual of every equation, e i (i 1, 2, · · · , m). b. Since Σ is typically unknown, the elements of Σ can then be estimated usinĝ c. Flexible generalized least square (FGLS) estimators are used to estimate the coefficientsβ where X V −1 X −1 is the variance-covariance matrix of the estimated coefficients.
The model is implemented using the SUR estimator in the SYSLIN procedure in SAS. The SAS code is available in supplementary materials. Table 1 presents the estimation results, where the red and blue colors highlight the positive and negative effects, respectively, both significant at 5% level. Overall, the 14-day lagged internal movement control policy is found to have significant negative effects (i.e., α 2i < 0) in 18 countries, which suggests the epidemics became more severe in these countries in spite of the introduction of movement control policies such as stay-at-home order.

Results
The correlations (and the changes of correlations) between the epidemics across different countries are presented in the 22 * 22 matrix under the coefficients α 3i j in Table 1. First, the USA, South Africa, Switzerland, the Netherlands, and Brazil are positively correlated (i.e., in red cells) with the epidemics in three or more countries, which means such correlations (with 14-day lag) have been weaker after introducing international travel restriction across these countries. At the same time, the daily new cases in each of these five countries are also negatively correlated (i.e., in blue cells) with the epidemics in at least three other countries, which implies the 14-day lagged epidemic became more severe in these five countries even after introducing the international travel restriction in respective countries. Overall, there were 30 significantly positive inter-country correlations and 29 significantly negative inter-country correlations as presented in Table 1, which indicates that in a global network the effect of international travel restriction on containing COVID-19 pandemic is not as straightforward as normally expected. The full estimation results are summarized in Table  S2.
As shown in Table 1, the international travel restriction for the travelers from China is found to have insignificant effect on the changes in the association between the epidemic in China and the outbreaks in the remaining 21 countries. That is, all the associations between China's cumulative infection number and the daily new infection numbers in other 19 countries with direct flights are insignificant (note that China has no direct flights to Brazil and Ecuador). The main reason may be due to China's early and immediate restriction not only on internal movement but also on nonessential international travel outbound from China in late January 2020. What is more, the number of new infections in China is not significantly correlated with the cumulative infections in other countries except Italy (α ITA→CHN 3 0.539, p < 0.001), which means that the correlation between the epidemic in Italy and daily new cases in China has been weaker after restricting foreigners entering China and such correlation between other countries (except Italy) and China has not changed significantly. Other similar cases are South Korea and Iran, in which the inter-country epidemic correlations have not changed significantly before and after the international travel restriction, given the insignificant correlations (1) between their cumulative infection numbers and Table 1 Estimation results The abbreviations for country names are provided in Table S1 other countries' lagged daily new cases and (2) between other countries cumulative infections and their own lagged daily new cases.
However, other Asian countries such as Japan and India experienced a different situation as they have more changes in their connections with other countries. The correlation between the cumulative cases in Japan and other countries' outbreaks has not changed significantly, but we find negative effects on the correlations between the lagged number of daily new cases in Japan and the cumulative infection numbers in the USA ( α USA→JPN 3 −0.123, p < 0.05) and Russia (α RUS→JPN 3 −0.112, p < 0.05). This means the correlation between these two countries' epidemic and the daily increment in Japan has been stronger, even after the international travel restrictions imposed on American and Russian travelers to visit Japan. Although the outbreak of COVID-19 in India was quite late, its epidemic size (i.e., cumulative number of infections) has a significant association with the numbers of daily new cases in European countries, such as Italy in the pre-controlled period for international travel (α IND→ITA 0.269, p < 0.01), which suggests the inter-country correlations between these countries epidemic and the daily new cases in the USA have been weaker after introducing international travel restrictions to the USA in March 2020. In comparison, even after the USA introduced the restriction for foreigners to enter, its lagged number of new infections is found to have stronger correlation with the cumulative infection sizes in several other countries than in the pre-controlled period, as indicated by the negative correlation coefficients with Spain (α ESP→USA the UK (α GBR→ESP 3 0.228, p < 0.001) in the pre-controlled period, and Italy (α ITA→ESP 3 −0.142, p < 0.05) in the controlled international travel period. Furthermore, the negative correlations between the epidemic in Spain and the numbers of new cases in the Netherlands and the USA indicate such association becomes stronger even after the international travel has been restricted from Spain to these two countries. Similarly, the number of new cases in Germany is significantly associated with the epidemics in the UK and the Netherlands as well as the USA and South Africa. Meanwhile, the epidemic in Germany also has strong correlations with the daily increment in Belgium, Switzerland, the USA, Brazil, and South Africa. Importantly, the correlations related to Germany are mostly positive, except with the Netherlands and South Africa, which implies most associations became weaker after imposing international travel restriction to/from Germany.
Our analysis also reveals the third-wave epidemics differ in West Asia (e.g., Turkey), East Europe (e.g., Russia), South America (e.g., Brazil and Ecuador), and Africa (e.g., South Africa). Specifically, the correlation between the lagged number of daily new cases in Turkey and the Canadian cumulative epidemic size became stronger in the controlled international mobility period than in the pre-controlled period (α CAN→TUR 3 −0.092, p < 0.01); the Turkish correlation between its cumulative infection number and the daily increment of South Africa was weaker with the international travel restriction, but its correlation with the daily increments in France, Switzerland, USA, and Brazil is stronger in the controlled international travel period than in the pre-controlled period. The cumulative epidemic in Russia is significantly associated with the recent upsurge of the COVID-19 new cases in Japan (α RUS→JPN 3 −0.112, p < 0.05). In comparison, the number of daily new cases in Brazil is strongly correlated with the cumulative epidemics in the USA, Germany, and South Africa in the pre-controlled international travel period, and with the Netherlands, France, and Turkey in the controlled period. In addition, the cumulative infection in South Africa has stronger correlation with the daily increment in the USA in the precontrolled period, but the daily new cases in South Africa are more correlated with the epidemic size in the USA in the controlled period. In addition, the daily increment in South Africa is more correlated with the cumulative infection numbers in Germany and Brazil in the controlled period, and with that of UK, the Netherlands, Switzerland, and Turkey in the pre-controlled period.
In summary, our analysis provides insightful findings on the changes in correlations of the COVID-19 pandemic across countries with the introduction of international travel restrictions, by taking into account the international human mobility. As illustrated in Fig. 2, the size of each node represents the total infection numbers in the corresponding country by April 24, 2020. A red (blue) line represents the positive (negative) correlation between each pair of countries with the strength of the line representing the absolute value of the standardized correlation coefficient.
To check the robustness of our results, we repeat our analysis by using the number of international flights instead of flight capacity for international mobility volume (i.e., FCap i jt ) in the model. The results are consistent with that in the main analysis and summarized in Table S3. We further consider the situation without flight information available but keeping the time-dependent spatial weight matrix (W t ) to indicate whether a country implemented international travel restriction with another country  Table S4 show that the correlation patterns we observe are largely held unchanged and may simulate the situation of international travels with transitions in third countries. Finally, we replicate the analysis by using another COVID-19 dataset maintained by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University [17] which contains more sources of information from CDC and online media sources in many countries. Results show largely consistent pattern of the correlation on the COVID-19 epidemics across countries and continents, although country-specific impact coefficients vary in the global network. For instance, China is still found to be fully independent from the epidemics in all other countries. Also, the USA is found to be consistently correlated with the European countries regarding their respective epidemics. The results based on the data from Johns Hopkins University are presented in Tables S5 to S7, corresponding to three models with the international flight capacity, the number of international flights and without flight information. The significant associations presented in Table  S5 are also illustrated in Figure S2.
Considering potential collinearities between the estimated coefficients, we conduct a stepwise approach that adds countries' variables continent by continent. Using the USA as an example, the stepwise analysis in Table S8 shows the consistent effects, that is, the correlations between other countries' epidemics and the daily new infections in the USA, remain largely same and robust regarding the effect sizes and significance levels.

Discussions
By using the global human mobility data, this study provides the first comprehensive insights into the correlation between the epidemics across a number of countries with the introduction of international travel restrictions during the first 3 months of the COVID-19 pandemic. Specifically, although the epidemic started in Asia first, the study finds that the epidemics in Asian countries have diverse degrees of associations with countries in other continents. China, for example, was one of the first countries to face an outbreak and severely affected, but introducing international travel restriction targeting travelers from China is found to have little change on the correlation between the early epidemic in China and the outbreak later in other countries. In contrast, the lagged daily infection in Japan shows more correlations with other counties' epidemics even with the international travel restriction. The study further reveals intensive correlations between the epidemics across North American and European countries either before or after the banning of trans-Atlantic travel. Among many European countries, we also find a close circle of mutual correlations in different periods of COVID-19 epidemics.
The first important takeaway from the study is the importance of a holistic view of the COVID-19 transmission. Restrictions on nonessential international human mobility between countries are one critical strategy to fight against epidemic outbreaks [6]; however, little was known before regarding the proper ways to implement the strategy. As such, policymakers might focus only on the hotspot areas with a high transmission rate (which is not wrong) but overlook the pathway effect beyond the epidemic areas in the global network. Unfortunately, over the past 3 months, we have observed such a narrow perspective repeatedly occurring in many countries, focusing only on cutting off hotspot outbreak areas. New York Governor Andrew Cuomo recently start to question where the coronavirus that hit New York state came from: "…We closed the front door with the China travel ban, which was right…but we left the back door open because the virus had left China by the time we did the China travel ban" [21]. Through this study, we show that the establishment of a new holistic and networked view of pandemic transmission on a global scale is imperative, with the following important implications.
First of all, it cannot be simply assumed that the international travel ban from highrisk areas will be sufficient to control the epidemic once and for all while ignoring the dynamic spread of the epidemic in the global network. Policymakers in every country should keep abreast of population movements between countries and, through scientifically rigorous analysis, and foresee the pattern of epidemic transmission in the network. Accordingly, they should dynamically adjust the corresponding international restriction strategy promptly.
Second, the global perspective should also take into account the domestic movement restriction already implemented and its effects. Despite the early outbreak in Asian countries, their impact on other countries is not the same. Because of strict national and international restrictions imposed by China, the epidemic in China was developing in isolation: It did not affect other countries or was affected by other countries' epidemics. In contrast, the lack of strict and timely domestic restriction in some other countries, combined with their lack of global view on timely international restriction, contributed to the pandemic, which means these countries were subsequently affected by the epidemic in the second and third waves.
By the same token, as epidemics are gradually brought under control in the coming months, the effects and extent of national and international movement restriction should both be taken into account when countries decide to resume international economic and social activities. At the same time, it's also important to remain vigilant about the third-wave outbreak that may appear in Latin America and Africa. The holistic view should be applied promptly to guide the internal and international movement restriction policies in the areas.
Finally, our analysis also shows that there may exist factors outside the epidemic that are influencing countries' decisions on international travel restrictions. For instance, when the USA started to ban non-US citizens or permanent residents who had been in China in the past 14 days to enter the USA on February 2, 2020, there were 16 640 infected patients in China with the infection rate of 11.95 per million people. In comparison, when the USA banned the entry of European people on March 13, Italy itself had 17 660 patients with an infection rate of 292.29 per million people. Therefore, this also adds another dimension of restriction timing that future research should be aware of.
As one of the first few studies to focus on a holistic view of the epidemic, the paper cannot avoid its limitations. First, our estimation model captures how potential infections spread as people move out of their destinations by multiplying the flight numbers and the number of cumulative infections in the destination countries. Since there is no reported infection number of China, either by WHO or JHU before January 22, the potential spread values from China to any other countries remain zero. It means that the study cannot estimate the spread situation before January 22 and the findings of the study are limited to the situation after January 22. Since the flight dataset shows continuous outbound flights from China to other countries after 22 January (albeit with a decrease in number), potential infection transmission from China to other countries can still be well captured by the mobility data.
Second, all the restrictions were still in effect at the time of writing; therefore, we mainly focused on the short-term effects rather than a longitudinal examination. We recommend that scholars exercise caution when extrapolating our conclusions to longer periods. Third, some reliability concerns still exist related to the confirmed number reported and policy implementation; however, we could analyze only the available data. Furthermore, examining the restriction implementation process and efforts of the actions in a more detailed manner would increase the rigor and power of the initial analysis presented in this paper. Models that can predict the dynamic future is also the focus of the next steps.