Estimating infection-related human mobility networks based on time series data of COVID-19 infection in Japan ============================================================================================================= * Tetsuya Yamada * Shoi Shi ## Abstract Comprehensive and evidence-based countermeasures against emerging infectious diseases have become increas-ingly important in recent years. COVID-19 and many other infectious diseases are spread by human movement and contact, but complex transportation networks in 21 century make it difficult to predict disease spread in rapidly changing situations. It is especially challenging to estimate the network of infection transmission in the countries that the traffic and human movement data infrastructure is not yet developed. In this study, we devised a method to estimate the network of transmission of COVID-19 from the time series data of its infection and applied it to determine its spread across areas in Japan. We incorporated the effects of soft lockdowns, such as the declaration of a state of emergency, and changes in the infection network due to government-sponsored travel promotion, and predicted the spread of infection using the Tokyo Olympics as a model. The models used in this study are available online, and our data-driven infection network models are scalable, whether it be at the level of a city, town, country, or continent, and applicable anywhere in the world, as long as the time-series data of infections per region is available. These estimations of effective distance and the depiction of infectious disease networks based on actual infection data are expected to be useful in devising data-driven countermeasures against emerging infectious diseases worldwide. ## Introduction **S**ince the identification of the coronavirus disease 2019 (COVID-19) in December 2019 [1, 2], over 180 million laboratory-confirmed infections and 4 million deaths have been observed worldwide as of 30 June 2021 [3]. The development of vaccines has certainly reduced the number of new infections worldwide [4, 5], allowing life before COVID-19 to gradually return to normal. However, in countries where vaccination has been delayed, there is a risk that the infection will continue to spread, and the emergence of new mutant virus strains has been observed [6, 7, 8], thus, proper measures against COVID-19 are still required. In the past year and a half, effective countermeasures against COVID-19 have been established, such as social distancing and wearing of face masks. The effectiveness of these non-pharmaceutical interventions (NPIs) has also been quantitatively evaluated [9, 10, 11, 12], together with the need for appropriate application where necessary depending on the active status of COVID-19. Otherwise, governments have to take the painful and severe measures again, the lockdowns, to prevent the spread of infection. Lockdowns are very effective infection control measures [13, 14, 10]; however, the trade-off between the impact on the economy and on the mental health of individuals cannot be ignored [15, 16, 17, 18]; thus, large-scale lockdown should be avoided as much as possible. As an alternative, selective lockdowns, which prevent the spread of infection by locking down a city when the number of infected individuals increases in that city, or other specific interventions have also been applied with some success in the United Kingdom and in other countries [19, 20, 21, 22, 23]. Since the spread of infection between cities is caused by the movement of people [24, 25, 26], infection does not always spread to cities that neighbor geographically, but spreads between cities where there is a continuous flow of individuals [27, 28]. Therefore, it is useful to create a network for predicting the spread of infection by calculating effective distances based on traffic and mobility of individuals [29, 30]. In fact, this method has been successful in predicting the spread of infection between countries and cities [31]. Conversely, effective distances are easier to estimate in countries where the human movement and traffic data is easy to be obtained and are well developed, although, this is limited to a few countries. Even in the countries or regions where the human movement and traffic data infrastructure is not yet developed, it is required to share a scheme for estimating infection networks based on the minimum necessary information and taking appropriate countermeasures worldwide. In addition, the power of effective distances is that they vary in a time-dependent and environment-dependent manner, as they are influenced by the number of infected people in the mobile population and the number of individuals taking preventive measures, which is difficult to be captured from movement-data-based network. In this context, it is useful to calculate effective distances and estimate the network structure of infection transmission based on the actual number of infected individuals. In this context, the dynamics of infected individuals in Japan may be useful to estimate effective distances of viral transmission. In the past year and a half, although there has been no explosion of infections similar to the tens of thousands of new infections and thousands of deaths per day observed in the United States and Europe, there has been a constant rate of infections observed in Japan. The inability to impose legally binding restrictions on the movement of people and the opening of commercial facilities has made it difficult to completely control the number of infected people; although, the cooperative nature of the Japanese might prevent the explosion of infections [32, 33, 34]. In other words, in Japan, it is possible to observe a discrepancy between effective distances based on movement data and on the number of infected people, which is difficult in countries or regions that have experienced a rapid increase in infection and repeated lockdowns that completely restrict movement. In addition, Japan has had three peaks of infection, two declarations of the state of emergency (i.e., with moderate restriction on the movement of people), and two stages of travel campaigns, during which the government encouraged domestic travels by subsidizing travel expenses, between 18 March 2020 and 12 March 2021, which have provided abundant data for describing changes in effective distances (Fig. 1). Thus, the unique dynamics of SARS-CoV-2 infection in Japan is a suitable model for calculating the effective distance of the number of infected people, and provides a new lockdown strategy based on an effective real-life-based network. ![Fig. 1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F1.medium.gif) [Fig. 1](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F1) Fig. 1 Timeline of the pandemic of coronavirus disease 2019 (COVID-19) in Japan. The period from 18 March 2020 to 12 March 2021 was considered. The top row indicates the period during which a state of emergency (i.e., “soft lockdown”, or mild restrictions on movement, alcohol-serving businesses, and large events) was declared. The second row indicates the period during which travel campaign, or “Go To Travel”, during which Japanese government encouraged the citizens to travel domestically by subsidizing travel expenses. The travel campaign had two stages: (1) the first from 22 July 2020 to 30 September 2020, during which Tokyo was excluded from the campaign, and (2) the second from 1 October 2020 to 27 December 2020, during which Tokyo was included and the subsidy was increased from stage 1. The main panel exhibits the 7-day backward moving average of the number of newly confirmed cases in each prefecture, which are highlighted in blue for Tokyo. Depending on the phase of the pandemic, the period was divided into six different periods. In this study, by combining the “Susceptible-Exposed-Infectious-Recovered” SEIR model and the diffusion equation [35, 36, 37], we simulated the infection propagation from Tokyo, the capital of Japan, to other pre-fectures and calculated the effective distance separating them. We then tested the effective distance model envisioning three different declarations of the state of emergency and a travel campaign. A map of Japan based on effective distance was drawn, and this map was used to predict the dynamics of infection associated with the Tokyo 2020 Olympics. ## Results ### Diffusion of infected population from Tokyo recapitulated the propagation of COVID-19 pandemic in Japan To model the propagation of the COVID-19 pandemic in Japan as a diffusion process of the infected population from Tokyo and to estimate effective distances, we first recapitulated a time series using the SEIR model in which the number of infected people in Tokyo depended on the distinct phases of the pandemic, or periods P1, P2, and P3 (Fig. 1). We calibrated the parameters for the number of potentially newly infected individuals in Tokyo during each period using the Markov chain Monte Carlo (MCMC) algorithm, and the summary of the estimated parameters is listed in Suplementary Table 1, which were generally consistent with previous studies on COVID-19 [35, 37]. Models using the estimated parameters accurately recapitulated the temporal dynamics and amplitude of the outbreaks and surges in infections (Supplementary Figure 1a-c). Leveraging these models as sources, we subsequently set up diffusion equations that were calculated as one-dimensional partial differential equations (PDEs). Thus, human mobility outcomes differed and depended on whether the pandemic was on a spreading or shrinking trend; hence, we solved the diffusion equations dividing each period into periods before and after a peak of infection in Tokyo, which generated six periods (i.e., P1a, P1b, P2a, P2b, P3a, and P3b) (Fig. 1). After numerically solving the PDEs, we fit the number of newly confirmed cases in each of the 47 prefectures of Japan to the solutions and estimated the effective distance from Tokyo (Fig. 2, Supplementary Figure 2). The calculated results at the estimated effective distance, particularly during the first part of each period (i.e., P1a, P2a, and P3a), recapitulated the temporal evolution of the number of newly confirmed cases in most prefectures (Fig. 2). Consequently, by considering the propagation of the pandemic as a diffusion process, we estimated the effective distance from Tokyo based on only the time series of newly confirmed cases of COVID-19 in each prefecture. ![Fig. 2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F2.medium.gif) [Fig. 2](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F2) Fig. 2 Effective distance from Tokyo is estimated from diffusion process. **a,b**, Solutions of diffusion equations for each period (a: before peaks of infection, b: after peaks of infection). Solutions for *x* ∈ {2*n* | 0 ≤ *n* ≤ 50, *n* ∈ Z} over time are shown. The source of the diffusion process (*x* = 0) is the solution of the SEIR model, whose parameters were calibrated based on the MCMC algorithm (Supplementary Figure 1a-c). ### The effective distance estimated based on the diffusion process from Tokyo revealed the mobility of individuals during the pandemic To characterize the estimated effective distance based on the diffusion process from Tokyo, we applied linear regression in which the geographical distances from Tokyo or effective distances based on inter-regional passenger traffic are explanatory variables [31]. As for the effective distance derived from passenger traffic data, we used inter-prefecture passenger volume data obtained in Japan in 2019, built a weighted directed network of mobility probability between prefectures, and estimated it as the shortest path from Tokyo in the graph. The effective distance based on passenger traffic network generally predicted the effective distance based on the diffusion process, whereas the geographical distance did not, indicating that human mobility from Tokyo to other prefectures largely explained the time series of infections (Fig. 3a, Supplementary Figure 3a-c). Notably, in the Kanto region of Japan, which consists of seven prefectures centered around Tokyo, both the distance metrics predicted the effective distance derived from the diffusion equations with high accuracy, implying that, around Tokyo, the dynamics of infections could be largely explained by inter-prefecture travel between Tokyo and other prefectures, which correlated with the geographical distance on a local scale (Fig. 3b, Supplementary Figure 3a). Conversely, in prefectures far away from Tokyo, there would be some factors other than the diffusion of infected individuals from Tokyo that contributed to dynamics of infections. To examine the factors that contributed to the deviations of the effective distance based on the diffusion process from the predicted values in Fig. 3a, we investigated residuals from the predicted values. The residuals of some regions exhibited unequal distributions around the baseline Supplementary Figure 4a,b. For example, residuals of the Kinki regions and the Tohoku regions mostly distributed above and below the baseline, respectively Supplementary Figure 4a,b. These results imply that the deviations from the predicted values attributed to the local, or regional factors other than travels from Tokyo, such as within-prefecture spread or containment of the coronavirus cases or travels from nearby prefectures. Therefore, our findings revealed that the effective distance based on the changes in evolution of infection over time in each prefecture putatively reflected both human mobility and other local effects from nearby prefectures during the pandemic. ![Fig. 3](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F3.medium.gif) [Fig. 3](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F3) Fig. 3 Effective distance derived from diffusion equations contains both passenger traffic and other local information that contribute to spread of the pandemic. **a,b**,Relationship between effective distance based on diffusion process from Tokyo and passenger traffic data in 2019 for all prefectures and (a) and the prefectures in Kanto region (b). A point represents each prefecture, and the prefectures whose effective distance is larger than 100 were excluded. Linear regression by least mean square and correlation coefficients are shown for each panel. ### Effective distance changes dynamically depended on stages of the pandemic Another advantage of the effective distance derived from time series data of infections over that derived from passenger volume data is that human mobility can be estimated dynamically without the need for extensive surveys. Such dynamic estimations based on different stages of the pandemic are particularly important to evaluate effects of specific policies (e.g., “soft lockdown” and travel campaigns) and to define a better strategy. Thus, to investigate changes in human mobility over different stages of the pandemic, we first examined the distribution of the effective distance over the respective periods. Histograms of effective distances revealed that the shapes of the distributions shifted depending on different periods, and particularly during spreading periods P1a, P2a, and P3a (Fig. 4a, Supplementary Figure 6a). Regarding the distribution of the effective distance, P2, during which the number of infected people declined without the need for any strict movement restrictions, our model showed that the largest mean and effective distance of each prefecture distributed widely around the mean (Fig. 4a). Conversely, P1 exhibited the smallest mean and deviation from the mean, whereas P3 exhibited an intermediate mean, but the distribution was not unimodal, indicating the existence of distinct trends in human mobility between prefectures (Fig. 4a). Considering that the Japanese government declared a state of emergency during P1 and P3 to contain the pandemic, the distribution of the effective distance of P2 would be considered desirable to limit the diffusion of coronavirus pandemic while maintaining passenger mobility. Furthermore, during P3, the effective distance in some prefectures became smaller while others did not, presumably reflecting the increase in travelers from Tokyo to other prefectures due to the second stage travel campaign in which travels to and from Tokyo became subject to the subsidization. These results demonstrated that infection-related human mobility could be captured during a pandemic. Subsequently, we investigated the changes in the effective distance for each prefecture by region. The general trend of the effective distance across the entire country was conserved over time (Fig. 4b and Supplementary Figure 6b), which indicated that, other than diffusion from Tokyo, there was a local connectivity between nearby prefectures that contributed to the spreading of the virus and altering their distance from Tokyo. Furthermore, considering the combination with other plots (Fig. 4a, Supplementary Figure 5a,b), it can be implied that closer effective distances, in particular like those found in major cities in each region, chiefly contributes to the rapid increase in coronavirus cases against which the government needed to impose restrictions. Consequently, the effective distance estimated from the time series of the newly confirmed cases evolved temporally, which might provide useful information for determining the appropriate mobility level during the pandemic. ![Fig. 4](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F4.medium.gif) [Fig. 4](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F4) Fig. 4 Effective distance changes dynamically depending on stages of the pandemic. **a**,Distribution of effective distances in periods P1a, P2a, and P3a. The vertical line in magenta indicates the mean value of the distribution. Prefectures whose effective distance is larger than 100 were excluded. **b**,Effective distance of each prefecture in periods P1a, P2a, and P3a, ordered by “Prefecture Code” defined International Organization for Standardization. Regions in which prefectures belong to are colored differently. ### A distorted map of Japan, based on effective distances and local interactions, re-vealed the non-uniform spread of the pandemic across Japan To visualize how the pandemic spread across Japan across different periods, we combined information on effective distances and local interactions between prefectures, which we defined as the geographical distance between adjacent prefectures that had large passenger traffic in the 2019 survey, to create a distorted map of Japan. The position (i.e., longitude and latitude) of a capital city in each prefecture was determined as the solution of a mathematical optimization equation using the non-linear least squares calculation [38]. In general, regions to the west of Tokyo, largely the Kinki and Kyusyu regions, reduced in size, whereas prefectures around Tokyo were expanded (for example, the Chubu regions were stretched in the north-west direction) in the distorted maps compared to the original one (Fig. 5a-d). These deformations indicated that the coronavirus tended to rapidly spread to the western regions of Japan, while it required additional time towards spread to the north and north-western directions, particularly towards the Tohoku and Chubu regions. Thus, distorting maps based on effective distance and local interactions would be useful to characterize the spread of the coronavirus geographically and to develop specific policies aimed at preventing a large spike in the infectious population. ![Fig. 5](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F5.medium.gif) [Fig. 5](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F5) Fig. 5 Distorted maps of Japan based on effective distance and local interactions reveal non-uniform spreading of the virus across Japan. **a**,Original map of Japan. Each point indicates the capital of each prefecture, colored by region. **b-c**,Distorted maps of Japan based on effective distance and local connectivity in periods P1a, P2a, and P3a (Supplementary Figure 7). ### Quantitative estimation of the scale of the pandemic in each prefecture in different sets of effective distances Although effective distance allows us to estimate how the virus spreads across Japan, it is still unclear how much the effective distance contributes to determine the scale of the pandemic quantitatively, which would be a useful parameter when developing policies that balance both economic activity and medical resources in each prefecture. To evaluate the effects of the effective distance, we conducted simulations in which only the effective distance was changed from P3. If the effective distance was the same as P1, the total number of infected people would increase to more than twice in peripheral or rural areas, but would not increase to a similar degree in populated prefectures (Fig. 6a). Conversely, if the effective distance was the same as P2, there would be no clear differences between the total number of infected individuals between the central or peripheral prefectures in Japan (Fig. 6a). These results indicated that the effective distance could largely affect the scale of the pandemic, and these effects were heterogeneous between prefectures. Finally, we considered the increase in human mobility that could result from a very large public event such as the 2020 Olympic and the Paralympic Games, which were to be held in Tokyo and several other prefectures, and estimated the impact on the pandemic (Supplementary Figure 9). Our findings showed that prefectures that held several matches, particularly prefectures that usually did not share strong connectivity with Tokyo, would experience a large increase in the number of patients (Fig. 6b). For example, the number of infected cases in Fukushima prefecture would increase by 80%, while other host prefectures would increase by 40 to 60%. This indicated that the impact of an increase in human mobility due to huge public events were more severe in prefectures that usually had a larger effective distance from Tokyo. Based on these quantitative estimations, we recommend that policymakers should correctly control for the effective distance or distribute medical resources across the whole country. Therefore, effective distance enables us to estimate the scale of the pandemic in a single prefecture based on changes human mobility and would be helpful in devising countermeasures against COVID-19 and other infectious diseases over changing situations. ![Fig. 6](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F6.medium.gif) [Fig. 6](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F6) Fig. 6 Simulations quantitatively estimates the effects of effective distance on a scale of the pandemic in each prefecture. **a**,Estimated increase or decrease of total infection during an equivalent period as P3 when the effective distance from Tokyo is set to the same as P1 (left) or P2 (right) and source (Tokyo) is the same as P3. **b**,Estimated increase in infections of infection during an equivalent period as P3 when every prefecture approaches Tokyo because of a large-scale event such as the Olympic. Some prefectures where some games are held are closer to Tokyo than others. The results of the simulated shifts are shown in Supplementary Figure 9. ## Discussion In this study, we calculated the effective distance of the spread of COVID-19 in Japan by applying the SEIR model to Tokyo, the capital city of Japan. Despite the presence of social factors that altered human mobility, such as two declarations of the state of emergency and travel campaigns, the spread of COVID-19 was demon-strated mainly by assuming the spread from Tokyo during all phases of the pandemic (Supplementary Figure 1). More specifically, our SEIR model-based predictions revealed that the effective distance from Tokyo to other prefectures changed over each phase, in which the effective distance between Tokyo and Hokkaido was close for P1 and P3, but far for P2, for example (Supplementary Figure 1). Comparing the fitting parameter and effective distance in P1a, P2a and P3a, we can see the effect of government-sponsored travel promotion. In par-ticular, the effective distance between Tokyo and prefectures with many tourist attractions has increased. These findings provided an explanation for the spread of the two outbreaks in Hokkaido and many other prefectures, which could not be explained merely by passenger traffic data alone. In Fig. 6, we estimated the impact of an huge public event such as the Tokyo Olympics by creating an infection model that took into account the increase in the mobility of people between prefectures having event venues. Since the mobility of people between Tokyo and each prefecture was predicted to be a critical factor for pandemic control [24, 25, 26], the first priority for policy makers should be to prevent the infectious situation in Tokyo from being transferred to local areas, and in this sense, the decision of the Japanese government to hold the Olympic games without spectators is commendable. Conversely, although, spectators have been excluded from participating in the events, athletes and media personnel remained a group of individuals moving around the event venues, so it is important to adequately consider the control strategy. Estimating the effective distance from actual infection data is also a useful method for constructing an infection network for the spread of infection, constructing all infection networks based on the SEIR model, and simulating the spread of infection when an outbreak occurs in a city is computationally infeasible. However, once the network structure is drawn, a feasible simulation becomes available. For example, in the case of this study, by assuming the spread of infection from Tokyo, we were able to explain the infection situation in other prefectures with high accuracy (Supplementary Figure 3), suggesting that the infection network in Japan was a simple radial structure centered on Tokyo. Conversely, if the same method is applied to the United States or European countries, where multiple cities or countries, respectively, experienced surges of infection simultaneously [39, 40], it may be possible to obtain a more complex structural organization of the infection network. By combining the effective distance estimated from infection data with the effective distance based on a transportation network and movement data, we believe that it is possible to establish more effective preventive policies to avoid the spread of infection. It is difficult to quantitatively evaluate the contribution of various modes of transportation, such as airplanes, cars, and trains, to infectious diseases. By comparing the two effective distances, it is possible to determine the means of transportation that contribute to the transmission of the disease and to take selective and effective preventive measures. However, the proposed method has some limitations. Importantly, the assumption of the diffusion process from Tokyo is qualitatively different between periods before and after peaks of infection and does not allow to quantitatively compare effective distances between these periods. Thus, improvements in methodology to the estimate effective distance based on the uniform assumption remain to be addressed in future studies. ## Methods ### Dataset We downloaded the data of the daily number of individuals tested positive in 47 prefectures from the COVID-19 Situation Report in Japan by the Toyo Keizai Online COVID-19 Task Team [41]. The data were originally provided by the Ministry of Health, Labor, and Welfare. Returnees on government charter flights from Wuhan, airport quarantine and passengers, and cases on the *Diamond Princess* cruise ship were excluded. Due to the data availability from all of the prefectures, we took into consideration the data from 18 March 2020 to 13 March 2021, and divided it into three periods based on the three waves of the pandemic: (P1) from 18 March 2020 to 16 May 2020, (P2) from 17 May 2020 to 3 October 2020, and (P3) from 4 October 2020 to 12 March 2021. Each of these periods was further subdivided into two periods (a and b) depending on the time of the peak of infection. The timeline of the study and the major announcements from the Japanese government are illustrated in Fig. 1. The survey data on passenger traffic between prefectures in 2019 were from the Ministry of Land, Infrastruc-ture, Transport, and Tourism. The data included the passenger volumes on railways, cars, ships, and airlines. The population data for each prefecture, the distance between prefectures, which is defined by the distance between prefectural capitals, and the GeoPackage of for Japan were downloaded from the Statistics Bureau of Japan, the Geospatial Information Authority in Japan, and Database of Global Administrative Areas (GADM), respectively. ### SEIR model of COVID-19 pandemic in Tokyo To describe the course of the COVID-19 pandemic in Tokyo, which is the capital and most populated city and likely the epicenter of the pandemic in Japan, we constructed a basic SEIR model, that included the following compartments: Susceptible (S), Exposed (E), Infectious (I), and Recovered (R). The choice of our model was motivated by several previously published studies on the COVID-19 pandemic [30, 35, 37]. The local dynamics of transmission in the model was given by ![Formula][1] In the model, the susceptible population (*S*) becomes exposed to the viral agent upon contact with the infectious population (*I*) and thus, becomes the exposed population (*E*). The transmission rate of the virus is expressed as *β*, and *βI* represents the force of infection. The exposed population becomes an infectious population at the rate *E* (i.e., calculated as the inverse of the latent period), and the infectious population becomes the recovered population (*R*) at the rate *ρ*, the or recovery rate. ### Parameter estimation using the MCMC algorithm Model parameters were estimated using the Bayesian framework by sampling the posterior parameter distri-bution via an affine-invariant MCMC algorithm using the emcee v3 toolkit [42, 43]. For data reliability, the data used to calibrate the model consisted of the 7-day backward moving average of newly confirmed cases in Tokyo. Considering the rapid increase in the spread of the coronavirus pandemic in Tokyo, we assumed that the count data of newly confirmed cases followed a Gaussian distribution with a mean given by *ξI* and variance given by 1*/τ*. *E* and *I*, the number of exposed and infectious individuals on the first day during the period, respectively, were also estimated. For simplicity, the number of recovered individuals, or *R*, was set at zero; thus, the number of susceptible individuals, or *S*, was calculated as *S* = *N**t* *E* *I* where *N**t* represented the total population of Tokyo. Consequently, we estimated seven parameters for each period: *β, E, ρ, ξ, E*, *I*, and *τ*. We used the same uniform distributions for all periods as the prior distributions for each parameter, as described in Suplementary Table 1. The statistics of the posterior distributions are also listed in Suplementary Table 1. The convergence of the MCMC chains are judged based on the integrated autocorrelation time in accordance with a previous study [42]. When the number of samples satisfied the inequality given by, ![Formula][2] where ![Graphic][3] (*M*) was the estimated integrated autocorrelation time, the relative error in the estimation of the posterior distribution was considered sufficiently small. ### Model of the spread of the infected population from Tokyo and estimation of ef-fective distance The spread of the COVID-19 pandemic in prefectures other than Tokyo was modeled as the diffusion process of infected individuals from Tokyo. We denoted the number of newly confirmed cases as *u*(*x, t*), as follows: ![Formula][4] where *D* was the diffusion coefficient. The parameters used in the simulation are summarized in [TABLE]. For simplicity, the initial conditions for the PDE were expressed as: ![Formula][5] which meant that the infectious population was located only in Tokyo in the initial condition. For the boundary conditions, we applied the Dirichlet boundary condition for *x* = 0 and the Neumann boundary condition for *x* = *L*, or ![Formula][6] The parameters for the PDE are summarized in Suplementary Table 2. The PDE was replaced by a series of ordinary differential equations (ODEs) by finite difference approximations, which were solved using the SciPy software. The solutions of the PDE, or the diffusion process of the infectious population from Tokyo depending on *x* and *t* values, are shown in Fig. 2a,b. To estimate the effective distance from Tokyo at different stages of the pandemic, we considered the diffusion process for six different periods depending on the period P1 to P3 and the peaks of the number of new cases in Tokyo during the periods: P1a, P1b, P2a, P2b, P3a, and P3b. The fitting of the observed data to the simulated data was also conducted for each period separately, and the *x*, or the effective distance of each prefecture from Tokyo based on the diffusion process, was determined by choosing *x* with a minimum root mean square error (RMSE) value. The estimated effective distance for all periods are summarized in Suplementary Table 3. If the effective distance of Tokyo, *x**T*, estimated from fitting was not equal to zero, the effective distance of all prefectures was adjusted by subtracting *x**T* for subsequent analyses. ### Effective distance based on the inter-prefecture network of traffic and mobility data of individuals Since the contagion process is predominantly caused by the movement of people, we can also estimate the effective distance between prefectures based on passenger traffic data. The effective distance derived by this method was used to test the validity of the effective distance derived from the diffusion equation and for its characterization (Fig. 3a,b, Supplementary Figure 3a). We defined the connectivity matrix P as follows: ![Formula][7] where *F**nm* represented the inter-prefecture passenger traffic on railways, cars, ships and airlines. Thus, the matrix P quantified the fraction of the passenger flux from prefecture *n* to *m*. Given this flux-fraction matrix P, the effective distance *d**nm* from prefecture *n* to *m* is defined as ![Formula][8] We then constructed a weighted directed graph network consisting of 47 nodes, representing the prefectures in Japan, and edges from prefecture *n* to *m* with weights of *d**nm* (i.e., 47*C*2 = 1081 edges in total). Consequently, the effective distance from Tokyo to prefecture *n*, or *D**n*, was estimated as the length of the shortest path from Tokyo to the prefecture *n* in the graph network. ### Linear regression of the effective distance based on the diffusion process by other distance metrics The effective distance derived from the diffusion equations was predicted by a linear model without intercept, or *y* = *ax*, as Tokyo is at the origin. Models were fit to minimize the residual sum of squares. Coefficient of determination (*R*2) was determined as ![Graphic][9], where *u* is the residual sum of squares and *v* is the total sum of squares. ### Transformation of a map of Japan based on the effective distance from Tokyo and local connectivity We created a transformed map of Japan according to the Mercator projection (i.e., the longitude and latitude corresponding to the *x* and *y* coordinates in the Euclidean plane, respectively) depending on the effective dis-tance from Tokyo, which was estimated from the diffusion equation and local interactions between prefectures (Supplementary Figure 7). Local interactions between prefectures were defined as the geographical distance between adjacent prefectures, including those connected by roads. Furthermore, for simplicity, the local con-nectivity was only considered for adjacent prefectures that had substantial passenger traffic, or 100,000 people per year in at least one direction. The unit of effective distance was converted to that of geographical distance assuming a linear relationship between them locally, as shown in Supplementary Figure 3c,d. In total, 111 links between prefectures (i.e., 46 links that expressed effective distance from Tokyo to other prefectures and 65 links of local interactions) were used to transform a map of Japan to represent the spread of the pandemic. Distorting a map based on a given proximity can be formulated as a nonlinear least squares problem as follows: ![Formula][10] where *ij* denoted the links between points *i* and *j*, and *x* and *y* represented the longitude and latitude of a prefecture, respectively. *L* denoted the set of 111 links between prefectures, and *t**ij* denoted the proximity (i.e., effective distance or local geographical distance). To solve this problem, we adopted the algorithm suggested by Shimizu and Inoue in 2009 [38]. In short, this algorithm considers the norm minimization of bearing changes instead of the norm minimization of variation of coordinates, which is solved by the Levenberg-Marquardt method, which uses bearings between points as the initial condition, and is a natural premise when transforming a map. Next, the nonlinear least squares problem was rewritten as ![Formula][11] where ![Graphic][12] was the approximate bearing of link *ij* in the transformed map. When *α* = 1, we obtained: ![Formula][13] This meant that we could approximate the solution for the initial nonlinear least squares problem by an iteration of solving two independent linear squares problems. When solving these linear squares problems, we assigned different weights to equations derived from the effective distance and those from local interactions (10:1). The weight was determined so that the distorted maps represented effective distance from Tokyo as well as geographical relationships between nearby prefectures (Supplementary Figure 8a,b). These weights enabled the effective distance from Tokyo to be expressed while retaining geographical information![Graphic][14]. When retaining geographical information, the solution was regarded as reaching convergence. ### Simulations to estimate effects of the effective distance on scales of the pandemic Period P3 was selected as the base period because we considered that the period reflected most recent human mobility, thus was suitable for estimating an impact of Tokyo 2020 Olympic and Paralympic on the pandemic. The rate of change was calculated as the ratio of estimated total number of infections in simulated distances to the total number of infections in the original distances (i.e., the effective distance during P3). The prefectures in which 2020 Olympic and Paralympic Games were to held (i.e., Tokyo, Kanagawa, Saitama, Chiba, Ibaraki, Shizuoka, Fukushima, Miyagi, and Hokkaido) were set three times closer than other prefectures, which modeled increased human mobility due to the games (Supplementary Figure 9). ## Data Availability All raw data and source codes used in this manuscript are available at [https://github.com/tyamadat/seir\_pde](https://github.com/tyamadat/seir_pde). [https://github.com/tyamadat/seir\_pde](https://github.com/tyamadat/seir_pde) ## Data availability Raw data used in this manuscript are available at [https://github.com/tyamadat/seir\_pde](http://github.com/tyamadat/seir_pde). ## Code availability All source codes used in this manuscript are open and available at [https://github.com/tyamadat/seir\_pde](http://github.com/tyamadat/seir_pde). ## Author contributions T.Y. and S.S. designed the study, conducted the simulations, analyzed the data, and wrote the manuscript. ![Supplementary Figure 1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F7.medium.gif) [Supplementary Figure 1](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F7) Supplementary Figure 1 SEIR model recapitulates time evolution of the number of newly confirmed cases in Tokyo for each period. **a-c**,Time evolution of infection derived from the SEIR model whose parameters are fit by the MCMC algorithm. ![Supplementary Figure 2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F8.medium.gif) [Supplementary Figure 2](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F8) Supplementary Figure 2 Diffusion process from Tokyo recapitulates time series data of infections in other prefectures during period P1a. Cumulative number of newly confirmed cases that were observed (blue) and estimated by the diffusion process from Tokyo (orange) in each prefecture during P1a. ![Supplementary Figure 3](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F9.medium.gif) [Supplementary Figure 3](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F9) Supplementary Figure 3 Relationship bewteen effective distance based on diffusion and other distances. **a**,Relationships between effective distance based on diffusion process from Tokyo and that on passenger traffic data in 2019 for the whole country (blue + orange) and Kanto region (orange). **b**,**c**,Relationships between effective distance based on diffusion process from Tokyo and Geographical distance from Tokyo for the all prefectures and (a) and the prefectures in Kanto region (b). A point indicates each prefecture. Linear regression by least mean square and correlation coefficient are shown for each panel. ![Supplementary Figure 4](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F10.medium.gif) [Supplementary Figure 4](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F10) Supplementary Figure 4 Residuals of the effective distance based on the diffusion process. Residuals of the effective distance derived from the diffusion process compared to the predicted values based on linear regressions based on effective distance derived from the passenger traffic data (Fig. 3a, Supplementary Figure 5a). ![Supplementary Figure 5](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F11.medium.gif) [Supplementary Figure 5](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F11) Supplementary Figure 5 Distribution of effective distance in different periods. **a-b**,Effective distance in each prefecture for each period, sorted in ascending order. ![Supplementary Figure 6](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F12.medium.gif) [Supplementary Figure 6](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F12) Supplementary Figure 6 Effective distance changes dynamically depending on stages of the pandemic. **a**,Distributions of effective distance in periods P1b, P2b, and P3b. The vertical line in magenta indicates the mean value of the distribution. Prefectures whose effective distance was larger than 100 were excluded. **b**,Effective distance of each prefecture in P1b, P2b, and P3b, ordered by “Prefecture Code” defined International Organization for Standardization. Regions in which prefectures belong to are colored differently. ![Supplementary Figure 7](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F13.medium.gif) [Supplementary Figure 7](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F13) Supplementary Figure 7 Local interactions with large passenger volumes. Links of adjacent prefectures in which more than 100,000 individuals traveled across the border in at least one direction in 2019 are shown. ![Supplementary Figure 8](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F14.medium.gif) [Supplementary Figure 8](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F14) Supplementary Figure 8 Distorted maps of Japan based on effective distance in different weight ratio between effective distance and local interactions. A map of Japan was distorted to represent the effective distance in each period. Distorted maps were determined as solutions for weighted non-linear least squares problems. Various weight ratio between the effective distance and the geographical distances between nearby prefectures which had high passenger volumes was tested (a, b, Fig. 5). ![Supplementary Figure 9](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/04/2021.08.02.21261486/F15.medium.gif) [Supplementary Figure 9](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/F15) Supplementary Figure 9 Putative map of Japan during the Tokyo 2020 Olympic used for the simulation in Fig. 6b. Distorted map based on effective distance during period P3a (left) and modified period P3a (P3’a: Olympic) (right). In period P3’a, the prefecture in which some games will be held, highlighted in the left panel, approached Tokyo more than others. View this table: [Supplementary Table 1.](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/T1) Supplementary Table 1. Parameters in the SEIR model and the estimated values. Seven parameters were estimated using the MCMC algorithm so that the SEIR model represented the time series data of infections in each period. The median of the sampler of each parameter was chosen as a representative parameter. View this table: [Supplementary Table 2.](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/T2) Supplementary Table 2. Parameters in the diffusion equations. The parameters for the diffusion equations were chosen so that the solutions represented the time course of the number of newly confirmed cases in each prefecture. View this table: [Supplementary Table 3.](http://medrxiv.org/content/early/2021/08/04/2021.08.02.21261486/T3) Supplementary Table 3. Effective distance in each prefecture in the respective periods. The effective distance of each prefecture from Tokyo was estimated in each period based on the diffusion process. ## Acknowledgements We would like to thank M. Katori for helpful discussions. This work was supported by JSPS KAKENHI Grant – Grant-in-Aid for Early-Career Scientists (19K16487 and 21K15136, S.S.), Transformative Research Areas A (20H05894 and 20H05903, S.S.), and Moonshot R and D – MILLENNIA Program Grant (JPMJMS2023-25, S.S.). * Received August 2, 2021. * Revision received August 2, 2021. * Accepted August 4, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1]. Fan Wu, Su Zhao, Bin Yu, Yan Mei Chen, Wen Wang, Zhi Gang Song, Yi Hu, Zhao Wu Tao, Jun Hua Tian, Yuan Yuan Pei, Ming Li Yuan, Yu Ling Zhang, Fa Hui Dai, Yi Liu, Qi Min Wang, Jiao Jiao Zheng, Lin Xu, Edward C. Holmes, and Yong Zhen Zhang. A new coronavirus associated with human respiratory disease in China. Nature, 579(7798):265–269, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2008-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 2. [2]. Chaolin Huang, Yeming Wang, Xingwang Li, Lili Ren, Jianping Zhao, Yi Hu, Li Zhang, Guohui Fan, Jiuyang Xu, Xiaoying Gu, Zhenshun Cheng, Ting Yu, Jiaan Xia, Yuan Wei, Wenjuan Wu, Xuelei Xie, Wen Yin, Hui Li, Min Liu, Yan Xiao, Hong Gao, Li Guo, Jungang Xie, Guangfa Wang, Rongmeng Jiang, Zhancheng Gao, Qi Jin, Jianwei Wang, and Bin Cao. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The Lancet, 395(10223):497–506, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0146736(20)30183-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 3. [3].World Health Organization. WHO coronavirus disease (COVID-19) dashboard, 2021. 4. [4]. Fernando P. Polack, Stephen J. Thomas, Nicholas Kitchin, Judith Absalon, Alejandra Gurtman, Stephen Lockhart, John L. Perez, Gonzalo Pérez Marc, Edson D. Moreira, Cristiano Zerbini, Ruth Bailey, Kena A. Swanson, Satrajit Roychoudhury, Kenneth Koury, Ping Li, Warren V. Kalina, David Cooper, Robert W. Frenck, Laura L. Hammitt, Özlem Türeci, Haylene Nell, Axel Schaefer, Serhat Ünal, Dina B. Tresnan, Susan Mather, Philip R. Dormitzer, Uğur şahin, Kathrin U. Jansen, and William C. Gruber. Safety and efficacy of the BNT162b2 mRNA Covid-19 Vaccine. New England Journal of Medicine, 383(27):2603–2615, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2034577&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 5. [5]. Lindsey R. Baden, Hana M. El Sahly, Brandon Essink, Karen Kotloff, Sharon Frey, Rick Novak, David Diemert, Stephen A. Spector, Nadine Rouphael, C. Buddy Creech, John McGettigan, Shishir Khetan, Nathan Segall, Joel Solis, Adam Brosz, Carlos Fierro, Howard Schwartz, Kathleen Neuzil, Lawrence Corey, Peter Gilbert, Holly Janes, Dean Follmann, Mary Marovich, John Mascola, Laura Polakowski, Julie Ledger-wood, Barney S. Graham, Hamilton Bennett, Rolando Pajon, Conor Knightly, Brett Leav, Weiping Deng, Honghong Zhou, Shu Han, Melanie Ivarsson, Jacqueline Miller, and Tal Zaks. Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine. New England Journal of Medicine, 384(5):403–416, 2021. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2035389&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 6. [6]. Adam S. Lauring and Emma B. Hodcroft. Genetic variants of SARS-CoV-2 - What do they mean? JAMA, 325(6):529–531, 2021. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 7. [7]. Salim S. Abdool Karim and T de Oliveira. New SARS-CoV-2 variants — Clinical, public health, and vaccine implications. New England Journal of Medicine, 384(19):1866–1868, 2021. 8. [8]. Pango Lineages. Global report investigating novel coronavirus haplotypes, 2021. 9. [9]. Derek K. Chu, Elie A. Akl, Stephanie Duda, Karla Solo, Sally Yaacoub, Holger J. Schünemann, Amena El-harakeh, Antonio Bognanni, Tamara Lotfi, Mark Loeb, Anisa Hajizadeh, Anna Bak, Ariel Izcovich, Carlos A. Cuello-Garcia, Chen Chen, David J. Harris, Ewa Borowiack, Fatimah Chamseddine, Finn Schünemann, Gian Paolo Morgano, Giovanna E.U. Muti Schünemann, Guang Chen, Hong Zhao, Ignacio Neumann, Jeffrey Chan, Joanne Khabsa, Layal Hneiny, Leila Harrison, Maureen Smith, Nesrine Rizk, Paolo Giorgi Rossi, Pierre AbiHanna, Rayane El-khoury, Rosa Stalteri, Tejan Baldeh, Thomas Piggott, Yuan Zhang, Zahra Saad, Assem Khamis, and Marge Reinap. Physical distancing, face masks, and eye protection to prevent person-to-person transmission of SARS-CoV-2 and COVID-19: a systematic review and meta-analysis. The Lancet, 395(10242):1973–1987, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)31142-9&link_type=DOI) 10. [10]. Seth Flaxman, Swapnil Mishra, Axel Gandy, H. Juliette T. Unwin, Thomas A. Mellan, Helen Coupland, Charles Whittaker, Harrison Zhu, Tresnia Berah, Jeffrey W. Eaton, Mélodie Monod, Pablo N. Perez-Guzman, Nora Schmit, Lucia Cilloni, Kylie E.C. Ainslie, Marc Baguelin, Adhiratha Boonyasiri, Olivia Boyd, Lorenzo Cattarino, Laura V. Cooper, Zulma Cucunubá, Gina Cuomo-Dannenburg, Amy Dighe, Bimandra Djaafara, Ilaria Dorigatti, Sabine L. van Elsland, Richard G. FitzJohn, Katy A.M. Gaythorpe, Lily Geidelberg, Nicholas C. Grassly, William D. Green, Timothy Hallett, Arran Hamlet, Wes Hinsley, Ben Jeffrey, Edward Knock, Daniel J. Laydon, Gemma Nedjati-Gilani, Pierre Nouvellet, Kris V. Parag, Igor Siveroni, Hayley A. Thompson, Robert Verity, Erik Volz, Caroline E. Walters, Haowei Wang, Yuanrong Wang, Oliver J. Watson, Peter Winskill, Xiaoyue Xi, Patrick G.T. Walker, Azra C. Ghani, Christl A. Donnelly, Steven Riley, Michaela A.C. Vollmer, Neil M. Ferguson, Lucy C. Okell, and Samir Bhatt. Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature, 584(7820):257–261, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 11. [11]. Laura Matrajt and Tiffany Leung. Effectiveness of social distancing for COVID-19. CDC, 26(8), 2020. 12. [12]. Colin J. Worby and Hsiao Han Chang. Face mask use in the general population and optimal resource allocation during the COVID-19 pandemic. Nature Communications, 11(1):1–9, 2020. 13. [13]. Neil M. Ferguson, Daniel Laydon, Gemma Nedjati-Gilani, Natsuko Imai, Kylie Ainslie, Marc Baguelin, Sangeeta Bhatia, Adhiratha Boonyasiri, Zulma Cucunubá, Gina Cuomo-Dannenburg, Amy Dighe, Ilaria Dorigatti, Han Fu, Katy Gaythorpe, Will Green, Arran Hamlet, Wes Hinsley, Lucy C. Okell, Sabine van Elsland, Hayley Thompson, Robert Verity, Erik Volz, Haowei Wang, Yuanrong Wang, Patric GT Walker, Caroline Walters, Peter Winskill, Charles Whittaker, Christl A Donnelly, Steven Riley, and Azra C. Ghani. Report 9: Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand. Imperial College London, 2020. 14. [14]. Shiyi Cao, Yong Gan, Chao Wang, Max Bachmann, Shanbo Wei, Jie Gong, Yuchai Huang, Tiantian Wang, Liqing Li, Kai Lu, Heng Jiang, Yanhong Gong, Hongbin Xu, Xin Shen, Qingfeng Tian, Chuanzhu Lv, Fujian Song, Xiaoxv Yin, and Zuxun Lu. Post-lockdown SARS-CoV-2 nucleic acid screening in nearly ten million residents of Wuhan, China. Nature Communications, 11(1):1–7, 2020. 15. [15]. Abdulkadir Atalan. Is the lockdown important to prevent the COVID-9 pandemic? Effects on psychology, environment and economy-perspective. Annals of Medicine and Surgery, 56(June):38–42, 2020. 16. [16]. Abel Brodeur, Andrew E. Clark, Sarah Fleche, and Nattavudh Powdthavee. COVID-19, lockdowns and well-being: Evidence from Google Trends. Journal of Public Economics, 193:104346, 2021. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 17. [17]. Shweta Singh, Deblina Roy, Krittika Sinha, Sheeba Parveen, Ginni Sharma, and Gunjan Joshi. Impact of COVID-19 and lockdown on mental health of children and adolescents: A narrative review with recommendations. Psychiatry Research, 293(May), 2020. 18. [18]. Jasper Verschuur, Elco E. Koks, and Jim W. Hall. Global economic impacts of COVID-19 lockdown measures stand out in highfrequency shipping data. PLoS ONE, 16(4th April):1–16, 2021. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0249644&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 19. [19]. Nils Haug, Lukas Geyrhofer, Alessandro Londei, Elma Dervic, Amélie Desvars-Larrive, Vittorio Loreto, Beate Pinior, Stefan Thurner, and Peter Klimek. Ranking the effectiveness of worldwide COVID-19 government interventions. Nature Human Behaviour, 4(12):1303–1312, 2020. 20. [20]. Akihiro Nishi, George Dewey, Akira Endo, Sophia Neman, Sage K. Iwamoto, Michael Y. Ni, Yusuke Tsugawa, Georgios Iosifidis, Justin D. Smith, and Sean D. Young. Network interventions for managing the COVID-19 pandemic and sustaining economy. Proceedings of the National Academy of Sciences of the United States of America, 117(48):30285–30294, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzQ4LzMwMjg1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMDQvMjAyMS4wOC4wMi4yMTI2MTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 21. [21]. Sheryl L. Chang, Nathan Harding, Cameron Zachreson, Oliver M. Cliff, and Mikhail Prokopenko. Modelling transmission and control of the COVID-19 pandemic in Australia. Nature Communications, 11(1):1–13, 2020. 22. [22]. Jan M. Brauner, Sören Mindermann, Mrinank Sharma, David Johnston, John Salvatier, Tomáš Gavenčiak, Anna B. Stephenson, Gavin Leech, George Altman, Vladimir Mikulik, Alexander John Norman, Joshua Teperowski Monrad, Tamay Besiroglu, Hong Ge, Meghan A. Hartwick, Yee Whye Teh, Leonid Chindelevitch, Yarin Gal, and Jan Kulveit. Inferring the effectiveness of government interventions against COVID-19. Science, 371(6531), 2021. 23. [23]. Matan Yechezkel, Amit Weiss, Idan Rejwan, Edan Shahmoon, Shachaf Ben-Gal, and Dan Yamin. Human mobility and poverty as key drivers of COVID-19 transmission and control. BMC Public Health, 21(1):1–13, 2021. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33451323&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 24. [24]. Moritz U G Kraemer, Chia-Hung Yang, Bernardo Gutierrez, Chieh-Hsi Wu, Brennan Klein, David M Pigott, Open Covid-19 Data, Working Group, Louis Du Plessis, Nuno R Faria, Ruoran Li, William P Hanage, John S Brownstein, Maylis Layan, Alessandro Vespignani, Huaiyu Tian, Christopher Dye, Oliver G Pybus, and Samuel V Scarpino. The effect of human mobility and control measures on the COVID-19 epidemic in China. Science, 368(6490):493–497, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80OTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8wNC8yMDIxLjA4LjAyLjIxMjYxNDg2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 25. [25]. Gregory A. Wellenius, Swapnil Vispute, Valeria Espinosa, Alex Fabrikant, Thomas C. Tsai, Jonathan Hennessy, Andrew Dai, Brian Williams, Krishna Gadepalli, Adam Boulanger, Adam Pearce, Chaitanya Kamath, Arran Schlosberg, Catherine Bendebury, Chinmoy Mandayam, Charlotte Stanton, Shailesh Bavadekar, Christopher Pluntke, Damien Desfontaines, Benjamin H. Jacobson, Zan Armstrong, Bryant Gipson, Royce Wilson, Andrew Widdowson, Katherine Chou, Andrew Oplinger, Tomer Shekel, Ashish K. Jha, and Evgeniy Gabrilovich. Impacts of social distancing policies on mobility and COVID-19 case growth in the US. Nature Communications, 12(1):6–12, 2021. 26. [26]. Pierre Nouvellet, Sangeeta Bhatia, Anne Cori, Kylie E.C. Ainslie, Marc Baguelin, Samir Bhatt, Adhiratha Boonyasiri, Nicholas F. Brazeau, Lorenzo Cattarino, Laura V. Cooper, Helen Coupland, Zulma M. Cucunuba, Gina Cuomo-Dannenburg, Amy Dighe, Bimandra A. Djaafara, Ilaria Dorigatti, Oliver D. Eales, Sabine L. van Elsland, Fabricia F. Nascimento, Richard G. FitzJohn, Katy A.M. Gaythorpe, Lily Geidelberg, William D. Green, Arran Hamlet, Katharina Hauck, Wes Hinsley, Natsuko Imai, Benjamin Jeffrey, Edward Knock, Daniel J. Laydon, John A. Lees, Tara Mangal, Thomas A. Mellan, Gemma Nedjati-Gilani, Kris V. Parag, Margarita Pons-Salort, Manon Ragonnet-Cronin, Steven Riley, H. Juliette T. Unwin, Robert Verity, Michaela A.C. Vollmer, Erik Volz, Patrick G.T. Walker, Caroline E. Walters, Haowei Wang, Oliver J. Watson, Charles Whittaker, Lilith K. Whittles, Xiaoyue Xi, Neil M. Ferguson, and Christl A. Donnelly. Reduction in mobility and COVID-19 transmission. Nature Communications, 12(1):1–9, 2021. 27. [27]. Pavel Skums, Alexander Kirpich, Pelin Icer Baykal, Alex Zelikovsky, and Gerardo Chowell. Global transmission network of SARS-CoV-2: From outbreak to pandemic. medRxiv, (March), 2020. 28. [28]. Shoi Shi, Shiori Tanaka, Ryo Ueno, Stuart Gilmour, Yuta Tanoue, Takayuki Kawashima, Shuhei Nomura, Akifumi Eguchi, Hiroaki Miyata, and Daisuke Yoneoka. Travel restrictions and sars-cov-2 transmission: An effective distance approach to estimate impact. Bulletin of the World Health Organization, 98(8):518–529, 2020. 29. [29]. Marino Gatto, Enrico Bertuzzo, Lorenzo Mari, Stefano Miccoli, Luca Carraro, Renato Casagrandi, and Andrea Rinaldo. Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures. Proceedings of the National Academy of Sciences of the United States of America, 117(19):10484–10491, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzE5LzEwNDg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMDQvMjAyMS4wOC4wMi4yMTI2MTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 30. [30]. Serina Chang, Emma Pierson, Pang Wei Koh, Jaline Gerardin, Beth Redbird, David Grusky, and Jure Leskovec. Mobility network models of COVID-19 explain inequities and inform reopening. Nature, 589(7840):82–87, 2021. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 31. [31]. Dirk Brockmann and Dirk Helbing. The hidden geometry of complex, network-driven contagion phenomena. Science, 342(6164):1337–1342, 2013. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNDIvNjE2NC8xMzM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMDQvMjAyMS4wOC4wMi4yMTI2MTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 32. [32]. Akiko Iwasaki and Nathan D Grubaugh. Why does Japan have so few cases of COVID-19? EMBO Molecular Medicine, 12(5):10–12, 2020. 33. [33]. Takahiro Yabe, Kota Tsubouchi, Naoya Fujiwara, Takayuki Wada, Yoshihide Sekimoto, and Satish V. Ukkusuri. Non-compulsory measures sufficiently reduced human mobility in Tokyo during the COVID-19 epidemic. Scientific Reports, 10(1):1–9, 2020. 34. [34]. Bo Yan, Xiaomin Zhang, Long Wu, Heng Zhu, and Bin Chen. Why do countries respond differently to COVID-19? A comparative study of Sweden, China, France, and Japan. American Review of Public Administration, 50(6-7):762–769, 2020. 35. [35]. Weston C. Roda, Marie B. Varughese, Donglin Han, and Michael Y. Li. Why is it difficult to accurately predict the COVID-19 epidemic? Infectious Disease Modelling, 5:271–281, 2020. 36. [36]. Alex Viguerie, Guillermo Lorenzo, Ferdinando Auricchio, Davide Baroli, Thomas J.R. Hughes, Alessia Patton, Alessandro Reali, Thomas E. Yankeelov, and Alessandro Veneziani. Simulating the spread of COVID-19 via a spatially-resolved susceptible–exposed–infected–recovered–deceased (SEIRD) model with heterogeneous diffusion. Applied Mathematics Letters, 111:106617, 2021. 37. [37]. Leonardo López and Xavier Rodó. A modified SEIR model to predict the COVID-19 outbreak in Spain and Italy: Simulating control scenarios and multi-scale epidemics. Results in Physics, 21, 2021. 38. [38]. Eihan Shimizu and Ryo Inoue. A new algorithm for distance cartogram construction. International Journal of Geographical Information Science, 23(11):1453–1470, 2009. 39. [39]. Jing Yuan, Minghui Li, Gang Lv, and Z. Kevin Lu. Monitoring transmissibility and mortality of COVID-19 in Europe. International Journal of Infectious Diseases, 95:311–315, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F04%2F2021.08.02.21261486.atom) 40. [40]. Anthony R. Ives and Claudio Bozzuto. Estimating and explaining the spread of COVID-19 at the county level in the USA. Communications Biology, 4(1), 2021. 41. [41]. Toyo Keizai Inc. Coronavirus disease (COVID-19) situation report in Japan, 2021. 42. [42]. J Goodman and J Weare. Ensemble samplers with affine invariance. Communications in Applied Mathematics and Computational Science, 5(1):65–80, 2010. 43. [43]. Daniel Foreman-Mackey, Will M. Farr, Manodeep Sinha, Anne M. Archibald, David W. Hogg, Jeremy S. Sanders, Joe Zuntz, Peter K. G. Williams, Andrew R. J. Nelson, Miguel de Val-Borro, Tobias Erhardt, Ilya Pashchenko, and Oriol Abril Pla. emcee v3: A Python ensemble sampling toolkit for affine-invariant MCMC. Journal of Open Source Software, 43(4):1864, 11 2019. [1]: /embed/graphic-7.gif [2]: /embed/graphic-8.gif [3]: /embed/inline-graphic-1.gif [4]: /embed/graphic-9.gif [5]: /embed/graphic-10.gif [6]: /embed/graphic-11.gif [7]: /embed/graphic-12.gif [8]: /embed/graphic-13.gif [9]: /embed/inline-graphic-2.gif [10]: /embed/graphic-14.gif [11]: /embed/graphic-15.gif [12]: /embed/inline-graphic-3.gif [13]: /embed/graphic-16.gif [14]: /embed/inline-graphic-4.gif