COVID-19 superspreading in cities versus the countryside ======================================================== * Andreas Eilersen * Kim Sneppen ## Abstract So far, the COVID-19 pandemic has been characterised by an initial rapid rise in new cases followed by a peak and a more erratic behaviour that varies between regions. This is not easy to reproduce with traditional SIR models, which predict a more symmetric epidemic. Here, we argue that superspreaders and population heterogeneity are the core factors explaining this discrepancy. We do so through an agent-based lattice model of a disease spreading in a heterogeneous population. We predict that an epidemic driven by superspreaders will spread rapidly in cities, but not in the countryside where the sparse population limits the maximal number of secondary infections. This suggests that mitigation strategies should include restrictions on venues where people meet a large number of strangers. Furthermore, mitigating the epidemic in cities and in the countryside may require different levels of restrictions. ## INTRODUCTION At its onset the COVID-19 pandemic shocked the world, with the number of new cases and deaths growing more than 20 % per day in the main hotspots [1]. With a growth rate this high, the disease was expected to spread through the population in less than six months without mitigation, and to reach a peak after three months, at which point 30 % of the population would have had the disease [2]. This, however, was not how the initial wave of the epidemic played out [3]. While most of the epidemic undoubtedly was halted due to mitigation efforts, it is striking that even in countries that have implemented a very light lockdown, such as Sweden, the epidemic has peaked long before herd immunity was achieved. Furthermore, even societies that slowly reopened businesses and public life did not experience an immediate explosive resurgence of the epidemic expected given the low levels of immunity and the speed with which the disease spread initially [3]. Here, we will propose an agent-based lattice model of an infectious disease that spreads in a geographically heterogeneous population. The model is a simplified depiction of the dynamics of COVID-19 with its characteristic parameters. We will examine the effect of the heterogeneous infection pattern that is so characteristic of this disease, using a gamma distributed infectiousness with dispersion factor *k* = 0.1 from [4]. Infection heterogeneity is a feature of several epidemic diseases [5], and plays a particularly important role for COVID-19 [4, 6–8]. In individual events, a single person has caused dozens of infections [9]. At the same time, the attack rate within households has been reported to be very low, at less than 20 % [10], despite prolonged close contact. This suggests that the majority of COVID-19 patients infect very little. In [11] an agent-based model was used to demonstrate that the Achilles heel of an epidemic driven by superspreaders was public social contacts, while the repeated contacts to smaller family and work groups were less dangerous. Looking at COVID-19 data from the United States in the analysis by [12], it is clear that a rapidly spreading epidemic occurs primarily in densely populated areas. In less densely populated areas, the epidemic onset is delayed, and in rural areas the epidemic never really starts, with most cases appearing to be spillover from the cities. The peak daily per capita mortality in [12] varies by a factor of ten between the most and least densely populated areas in the USA. The density dependence of COVID epidemics could in principle be explained by the higher chance of meeting infected people in dense areas. However, only 1% of COVID-19 spread was outdoors in China [13], suggesting that only indoor meetings counts. Furthermore, only few infections happens within households [10], suggesting that it is social visits and meetings in confined areas that facilitate infection [14]. We will base our model on these observations. ## MODEL Our model plays out on a lattice of side *L* with periodic boundary conditions. Each lattice site may either be empty or contain one agent. The agents can be in one of three states, susceptible, infectious, or recovered. Individuals interact with their neighbours with a frequency *f**meet* which is fixed along with the mean infection probability per meeting to give the desired number of secondary infections, *R*. The distance that agents travel we draw from the distribution ![Graphic][1] with a mean of *r* = 10 sites unless stated otherwise. The real probability distribution of travel distances for cars has in Italy been determined to be an exponential function below 20 kilometers, and then a powerlaw with a steeper cutoff around 500 kilometers [15]. The interaction radius *r* is not by itself a meaningful parameter. Rather, it is the number of neighbours within this radius, given by ![Graphic][2] (where *ρ* is the population density) which determines the behaviour of the system, as will be demonstrated below. COVID-19 is known to be transmitted before the onset of symptoms [16]. For simplicity, rather than using the SEIR modeling framework, we combined the exposed, presymptomatic state and the overlapping infectious period (both estimated at about five days) into one single infectious period of 1*/γ* = 10 days [17, 18]. Infected agents randomly leave the infectious state with a probability of *γ* per day, meaning that the duration of the disease is exponentially distributed. When an infectious agent *i* interacts with a susceptible agent, the susceptible agent will become infected with a probability *p**i* that is specific to the infecting agent. We draw these probabilities from a gamma distribution with dispersion factor *k* = 0.1 [4, 5], within the range observed for COVID-19 [4, 6, 8]. The distribution is normalised to give an average reproduction number *R* of 3 at a population density of 1. The geographical heterogeneity of the population is modeled by placing a square “city” of side *L/*5 on a lattice with periodic boundary conditions and lattice size *L*. The city has the population density 1, i.e. all sites are occupied, and the city population is thus *L*2/25. The city is surrounded by “countryside” with a population density *ρ* and total population of 24*ρL*2/25. Importantly we assume that the rate of interactions per agent is kept the same in both city and countryside, meaning that we assume that people are equally social. If an agent in the countryside attempts to interact with an empty site, the attempt is counted as failed, and new attempts are made until the number of contacts is the same as in the city, where all sites are occupied. Thus, people in the countryside interact with a smaller set of people while still spending the same amount of time on social activities. In a wider perspective this proposes that density dependence of disease spreading is more due to difference in diversity of contacts than due to differences in time spent around other people. Thereby our model assumes an infection rate that depends on density, but not in a simple linear fashion as sometimes assumed [19]. We seed the disease within the city. Even if we were to seed it randomly, the city would usually be hit early on provided that the epidemic catches on. ## RESULTS In fig. 2 we consider a disease with heterogeneous infection rates and study how simulated epidemic trajectories depend on population density. One sees that the dynamics resembles that of a SIR model at *ρ* values much lower and much higher than the critical density. In the low regime, the epidemic only spreads in the city which is relatively well-mixed, whereas in the high regime the countryside begins to resemble the city more and more. In both extremes, the fraction of infected individuals rise and fall symmetrically. The figure illustrates how the epidemic is “stretched out” in the intermediate density range. The epidemic has the longest lifetime when *ρ* is just above the percolation threshhold *ρ**crit*, as the disease still spreads in the countryside, but is nonetheless slowed down by te lower population density. ![FIG. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F1.medium.gif) [FIG. 1.](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F1) FIG. 1. Model: A superspreader in a city interacts a little with a lot of people and will infect some fraction of them. On the other hand, a superspreader outside the city will interact a lot with each of a smaller set of people. The superspreader then infects practically all of them, but there is a lower cap on the number of secondary infections. ![FIG. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F2.medium.gif) [FIG. 2.](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F2) FIG. 2. Epidemic trajectories: Infected fraction of the population over time changes with the countryside population density *ρ*. In the low-*ρ* regime, increasing the population density stretches the curve, as the epidemic spreads further from the city. When *ρ* is above ∼0.06, the epidemic again approaches the behaviour of a SIR model, as the epidemic now spreads unhindered across the whole system. Around *ρ**crit*, there is a large variation in the duration of the epidemic. The parameters used are *γ* = 0.1, *r* = 10, *f**meet* = 10 and dispersion parameter *k* = 0.1. In fig. 3 (a, b) we measure the attack rate of the epidemic in a homogeneously distributed population in order to find the percolation threshhold, below which the epidemic will stop propagating. Panel (a) identifies *ρ**crit* ≈0.01 for a homogeneously spreading disease. In contrast, a disease with an overdispersion of *k* = 0.1 has a much higher critical density *ρ**crit* ≈0.04, as seen in (b). The same figure demonstrates that it is not the density alone that determines the ability of the disease to percolate, but rather the number of neighbours, proportional to ![Graphic][3]. As expected, the disease with homogeneous infectivity can spread already when each person has one neighbour on average, whereas the overdispersed disease requires almost four neighbours to percolate. The overdispersed simulation and the simulation with homogeneous infection were done with same average disease transmission rate and the factor ∼4 difference in critical density comes about because a disease with *k* = 0.1 has 10 % of the infected being responsible for 80 % of the infections. Thus most people do not transmit the disease, and it is therefore the density of the few people who do spread the disease that sets the critical threshhold. ![FIG. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F3.medium.gif) [FIG. 3.](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F3) FIG. 3. Comparison of models without and with superspreaders. Panel (a) shows the attack rate as a function of the number of neighbours within the radius of interaction ![Graphic][4] in a population where everyone infects with the same rate. (b) shows attack rate with heterogeneous infection rates, using a gamma distribution with dispersion factor *k* = 0.1. The two overlaid curves demonstrate that the parameter *r* does not affect the physics of the system, and what really determines the ability of the disease to percolate is the number of neighbours, proportional to ![Graphic][5]. (c) Epidemic trajectory when superspreaders dominate (blue) and when infectiousness is evenly distributed (red) for equal countryside population density and radius of interaction ![Graphic][6]. When superspreaders are the main drivers of the epidemic, it is strongly impeded once the city has reached herd immunity. When everyone infects equally, the epidemic simply spreads radially out from the city, leading to a “second wave” in the countryside. Parameters are as in fig. 2. Since our model analysis centers on superspreaders as a main driver of the epidemic, fig. 3(c) compares epidemics with and without superspreaders. It can be seen that with no superspreaders, the epidemic will spread un-hindered in the countryside, albeit more slowly since the countryside is geographically larger than in the city. This leads to a graph similar to two superimposed SIRlike models. If superspreaders are present, however, the epidemic may spread both slowly and erratically in the countryside and continue long after herd immunity is achieved in the city. The dependence of the percolation threshhold on the dispersion parameter *k* is further investigated in fig. 4, which shows the attack rate as a function of both the number of neighbours and *k*. The figure shows that the percolation threshhold increases drastically at low *k*, meaning that a superspreader-driven epidemic requires more social contacts per person in order to spread. Once the epidemic becomes sufficiently overdispersed (*k <* 0.05), it is no longer viable. ![FIG. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F4.medium.gif) [FIG. 4.](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F4) FIG. 4. Dependence of attack rate on density and *k*. Since the variable determining percolation is not the absolute density, but the number of neighbours, we plot ![Graphic][7] on the x-axis rather than *ρ*. It is seen that the disease percolates much more easily at a higher *k*, which implies a more homogeneous infectivity. The more overdispersed the disease (corresponding to lower *k*), the more neighbours are required for percolation. In fig. 5, we try to replicate the data compiled by [12] and see that local disease incidence in a model with heterogeneous infectivity is indeed much more population density dependent than a model assuming homogeneous infectivity. This fits well with the cited data, which suggest a strong dependence of COVID-19 incidence on population density. It has already been known for years that the spread of epidemics is population density dependent [19]. Here, we show that this dependence is enhanced by heterogeneous infectivity. Importantly, as opposed to traditional disease models, we assume that everyone is equally social, but that the set of available contacts is smaller in sparsely populated regions. The significance of this will be discussed further below. ![FIG. 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F5.medium.gif) [FIG. 5.](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F5) FIG. 5. Epidemic trajectory in regions of varying density. When infectivity is homogeneous (top), the epidemic is a lot less sensitive to a lower population density than when the epidemic is driven by superspreaders (bottom). Here, the epidemic is nearly absent in the low-density regions, and appears to be driven by spillover. Inset illustrates the layout of the lattice. ![FIG. 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F6.medium.gif) [FIG. 6.](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F6) FIG. 6. Simulation with multiple cities: Infected individuals are shown in red, susceptibles are green, empty sites are white and recovered are black. City size distribution mimics Zipf’s law [23], such that there is one city of with 40000 in-habitants, 4 with 10000, and so forth. The graphs show the fraction of the population currently infected as a function of time for three example runs of the simulation. Here, *ρ* = 0.03 and the other parameters are the same as the above figures. The delayed onset and erratic behaviour of the countryside epidemic obviously depends on the density and other characteristics of the countryside and the assumed travelling pattern of individuals. Also, real countryside contains a diverse pattern of smaller and large settlements. Therefore we considered a system with several cities distributed on a lattice with side length *L* = 1000. The system mimics the observed city size distribution which is fairly close to the Zipf law [20]. The random distribution of smaller cities is not entirely naturalistic, since evidence suggests that real cities are organised in a fractal pattern [21, 22]. A figure using a city distribution closer to a fractal can be found in the supplement. For illustrative purposes we chose a density in the countryside that is close to but below the percolation threshhold for the disease (*ρ* = 0.03). With this below-critical spreading we observe cities that are spared and cities that are nearly completely infected. In reality assuming near critical spreading in the countryside would not be not necessary for global spread, since people occasionally travel long distances [15], facilitating rare direct transmissions between distant cities. ## DISCUSSION Based on the above, we suggest that the lopsided appearance of COVID-19 epidemic curves can be explained by heterogeneous infection ability combined with a geographically heterogeneous population. However, our model makes a number of assumptions and breaks from reality, whose importance must be discussed. First and foremost, the distribution of full and empty lattice sites does not represent the geographical distribution of people as such. Rather, it represents the density of contacts. The density ratio between people in the countryside compared to an inner city is much lower than used in our model. On the other hand, persons living outside the cities may be more mobile than city-dwellers. The density contrast of agents on our lattice represents the combined effect of these factors. A further complication is the distribution of disease duration and incubation periods. We here assume a simple exponential distribution of infectious period duration, whereas the real mechanics of COVID-19 includes presymptomatic transmission and a broad gamma distribution of incubation periods [17]. A different infectious period distribution might complicate our findings. Therefore, we examine the effect of a gamma distributed infectious period duration in the supplement. We find that, while changing this distribution has an effect on the percolation threshhold, it does not change our fundamental conclusions. It has long been discussed how disease transmission rate depends on population density [19, 24, 25]. We here present a new way of looking at this problem. A classical SIR model with pseudo-mass action transmission would assume a simple linear dependence of *R* on population density, implicitly assuming that people become more social by living in a densely populated area. Our model makes a different assumption: People will be equally social regardless of population density, but when population density is lower, the groups that they spend time with will be less diverse. The number of possible unique contacts for each individual declines linearly with population density, but people in sparsely populated areas are likely to have multiple encounters with the same persons, as their communities are smaller. If, however, the epidemic is dominated by superspreaders who only need one or a few encounters to transmit the disease, the duration of each encounter becomes less important as even a rather brief contact to a superspreader is likely to lead to infection. Instead, what limits the action of superspreaders is their number of unique contacts. Superspreaders that interact with only a small, tight-knit group will inherently be highly limited in how many secondary infections they can generate. If superspreading makes a disease vulnerable to variations in population density, we should conversely expect to see that diseases with a homogeneous infectivity, i.e. a dispersion factor close to or above 1, exhibit little variation with population density. One example of a disease with a homogeneous infectivity is influenza, with an estimated dispersion factor of *k* = 0.94 for the 1918 pandemic flu [26]. When examining the incidence of seasonal influenza, which we assume to have a similar dispersion factor, [27] find no consistent variation with population density. This is a point in favour of the link between superspreading and population density dependence. Finally, the effect of lockdowns and changes of social behaviour is important. A previous paper [11] suggests that even moderate mitigation may limit the action of superspreaders by reducing the maximal number of people any person can be in contact with. If this hypothesis is true, bans on gatherings and a reduction in public social life would lead to early peaks in the number of new cases. Our study compounds this finding, and suggests that a change in behaviour is not strictly necessary to cause an epidemic peak well before herd immunity has been achieved. Mitigation strategies that primarily target cities may well be sufficiently effective in bringing down the epidemic. However, large events like funerals, weddings, or festivals are not included in our model and will of course facilitate spreading in any location. Despite some caveats, our model reproduces the main aspects of geographical heterogeneity and suggests a new view on density dependence of epidemic dynamics. An epidemic with a large heterogeneity in infection rates is predicted to be most intense in large cities while it slowly tapers off in the countryside. This is consistent with what we see in data from the COVID-19 pandemic [12, 28]. Our results thus favour the hypothesis that the COVID-19 pandemic is driven by superspreaders, and that the observed quick exponential growth phase, early peak, and slow recovery phase are consequences of combining heterogeneity of infectivity with a heterogeneous population density. ## Data Availability An implementation of the model described in this article can be found on figshare under the DOI https://doi.org/10.6084/m9.figshare.12919709.v1 [https://doi.org/10.6084/m9.figshare.12919709.v1](https://doi.org/10.6084/m9.figshare.12919709.v1) ## supplement ### Different distribution of infectious periods In order to take into account that the duration of illness caused by COVID-19 does not actually follow a neat exponential distribution, we have here considered the case where the duration of infection is gamma distributed as well, though with a dispersion factor greater than one (*k* ≈4). We repeat the central figures 3 and 5 from the main article in order to see if a different distribution of infectious period durations changes our results. The gamma distribution should give fewer people with a very short disease duration, though both distributions have a relatively long tail. The mean duration of infectiousness is the same as in the main article (10 days). An important complication, that individuals with a long period of illness are likely to be hospitalised or otherwise isolated for part of it and thus less infectious, has not been accounted for here. We find that, while the different distribution of infectious periods does lower the percolation thresh-hold slightly and also diminishes the population density dependence of the epidemic, it does not fundamentally change our conclusions. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F7.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F7) Figure 1: (a) and (b) show attack rates as a function of ![Graphic][8], a measure of the number of neighbours within a radius of *r*. (a) shows the case with no infectivity overdispersion (*k* = 1) and (b) shows the case where infectivity is overdispersed with *k* = 0.1. We see that the system behaves similar to the main model, but with a slightly lower percolation threshhold in both cases. (c) shows the probability distribution function for disease duration used here. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F8.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F8) Figure 2: A comparison of the effect of population density on an epidemic that is not overdispersed (top) and an overdispersed epidemic (bottom). As above, the duration of illness is gamma distributed. The effect of density is somewhat reduced in the overdispersed case, but still noticeably greater than the non-overdispersed case. The smaller number of people with very short disease durations using this distribution seemingly also helps the epidemic spread more effectively. ### Fractal city distribution The random distribution of cities used in fig. 6 of the main article is not entirely realistic. Rather, there is evidence suggesting that real geographical distribution of cities follows a fractal pattern [1,2]. We therefore repeat the figure here with a (partially) fractal distribution of cities. The results of this simulation are fairly similar to the random city distribution. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.09.04.20188359/F9.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/10/14/2020.09.04.20188359/F9) Figure 3: A simulation where, instead of a random distribution, the cities are laid out in a fractal pattern. We see that this makes little difference from the random distribution shown in the main article. The lower panel shows epidemic trajectories from three test runs. ## ACKNOWLEDGMENTS We thank Lone Simonsen, Viggo Andreasen, and Bjarke Frost Nielsen for enlightening discussions. Our research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No. [740704]. ## Footnotes * Added figure exploring the connection between overdispersion factor k and population density * Received September 4, 2020. * Revision received October 13, 2020. * Accepted October 14, 2020. * © 2020, 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].Remuzzi and G. Remuzzi, Covid-19 and italy: what next?, The Lancet (2020). 2. [2]. R. M. Anderson, H. Heesterbeek, D. Klinkenberg, and. D. Hollingsworth, How will country-based mitigation measures influence the course of the covid-19 epidemic?, The Lancet 395, 931 (2020). 3. [3]. E. Dong, H. Du, and L. Gardner, An interactive web-based dashboard to track covid-19 in real time, The Lancet infectious diseases 20, 533 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30120-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.09.04.20188359.atom) 4. [4]. A. Endo, S. Abbott, A. J. Kucharski, S. Funk, et al., Estimating the overdispersion in covid-19 transmission using outbreak sizes outside china, Wellcome Open Research 5, 67 (2020). 5. [5]. J. O. Lloyd-Smith, S. J. Schreiber, P. E. Kopp, and W. M. Getz, Superspreading and the effect of individual variation on disease emergence, Nature 438, 355 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature04153&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16292310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.09.04.20188359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233300200048&link_type=ISI) 6. [6]. Q. Bi, Y. Wu, S. Mei, C. Ye, X. Zou, Z. Zhang, X. Liu, L. Wei, S. A. Truelove, T. Zhang, et al., Epidemiology and transmission of covid-19 in 391 cases and 1286 of their close contacts in shenzhen, china: a retrospective cohort study, The Lancet Infectious Diseases (2020). 7. [7]. D. Adam, P. Wu, J. Wong, E. Lau, T. Tsang, S. Cauchemez, G. Leung, and B. Cowling, Clustering and superspreading potential of severe acute respiratory syndrome coronavirus 2 (sars-cov-2) infections in hong kong, Research Square - preprint (2020). 8. [8]. D. Miller, M. A. Martin, N. Harel, T. Kustin, O. Tirosh, M. Meir, N. Sorek, S. Gefen-Halevi, S. Amit, O. Vorontsov, et al., Full genome viral sequences in-form patterns of sars-cov-2 spread into and within israel, medRxiv (2020). 9. [9]. Y. Shen, C. Li, H. Dong, Z. Wang, L. Martinez, Z. Sun, A. Handel, Z. Chen, E. Chen, M. H. Ebell, F. Wang, Yi, H. Wang, X. Wang, A. Wang, B. Chen, Y. Qi, L. Liang, Y. Li, F. Ling, J. Chen, and G. Xu, Community outbreak investigation of sars-cov-2 transmissionamong bus riders in eastern china, JAMA Internal Medicine (2020). 10. [10]. S. Y. Park, Y.-M. Kim, S. Yi, S. Lee, B.-J. Na, C. B. Kim, J.-i. Kim, H. S. Kim, Y. B. Kim, Y. Park, et al., Early release-coronavirus disease outbreak in call center, south korea, Emerging Infectious Diseases (2020). 11. [11]. K. Sneppen, R. J. Taylor, and L. Simonsen, Impact of superspreaders on dissemination and mitigation of covid-19, medRxiv (2020). 12. [12]. N. Medhi, Covid-19 cases in the u.s. by population density, [https://www.endcoronavirus.org/projects-1/covid-19-cases-in-the-us-by-population-density](https://www.endcoronavirus.org/projects-1/covid-19-cases-in-the-us-by-population-density) (2020), accessed: 20-08-2020. 13. [13]. C. J. Carlson, J. D. Chipperfield, B. M. Benito, R. J. Telford, and R. B. O’Hara, Species distribution models are inappropriate for covid-19, Nature Ecology & Evolution 4, 770 (2020). 14. [14]. H. Nishiura, H. Oshitani, T. Kobayashi, T. Saito, T. Sunagawa, T. Matsui, T. Wakita, M. Covid, and M. Suzuki, Closed environments facilitate secondary transmission of coronavirus disease 2019 (covid-19), medRxiv (2020). 15. [15]. L. Pappalardo, S. Rinzivillo, Z. Qu, D. Pedreschi, and F. Giannotti, Understanding the patterns of car travel, The European Physical Journal Special Topics 215, 61 (2013). 16. [16]. L. C. Tindale, J. E. Stockdale, M. Coombe, E. S. Garlock, W. Y. V. Lau, M. Saraswat, L. Zhang, D. Chen, J. Wallinga, and C. Colijn, Evidence for transmission of covid-19 prior to symptom onset, Elife 9, e57149 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.57149&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32568070&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.09.04.20188359.atom) 17. [17]. H. Nishiura, N. M. Linton, and A. R. Akhmetzhanov, Serial interval of novel coronavirus (covid-19) infections, International journal of infectious diseases (2020). 18. [18]. Q. Li, X. Guan, P. Wu, X. Wang, L. Zhou, Y. Tong, R. Ren, K. S. Leung, E. H. Lau, J. Y. Wong, et al., Early transmission dynamics in wuhan, china, of novel coronavirus–infected pneumonia, New England Journal of Medicine (2020). 19. [19]. R. M. Anderson, B. Anderson, and R. M. May, Infectious diseases of humans: dynamics and control (Oxford university press, 1992). 20. [20]. X. Gabaix, Zipf’s law for cities: an explanation, The Quarterly journal of economics 114, 739 (1999). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1162/003355399556133&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000082063300002&link_type=ISI) 21. [21]. S.-H. Yook, H. Jeong, and A.-L. Barabási, Modeling the internet’s large-scale topology, Proceedings of the National Academy of Sciences 99, 13382 (2002). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiOTkvMjEvMTMzODIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMC8xNC8yMDIwLjA5LjA0LjIwMTg4MzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 22. [22]. T. Mori, T. E. Smith, and W.-T. Hsu, Common power laws for cities and spatial fractal structures, Proceedings of the National Academy of Sciences 117, 6469 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE3LzEyLzY0NjkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMC8xNC8yMDIwLjA5LjA0LjIwMTg4MzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 23. [23]. K. T. Soo, Zipf’s law for cities: a cross-country investigation, Regional science and urban Economics 35, 239 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.regsciurbeco.2004.04.004&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228611400002&link_type=ISI) 24. [24]. B. Finkenst ä dt, M. Keeling, and B. Grenfell, Patterns of density dependence in measles dynamics, Proceedings of the Royal Society of London. Series B: Biological Sciences 265, 753 (1998). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1998.0357&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9628034&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.09.04.20188359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000073580200003&link_type=ISI) 25. [25]. B. Finkenst ä dt and B. Grenfell, Empirical determinants of measles metapopulation dynamics in england and wales, Proceedings of the Royal Society of London. Series B: Biological Sciences 265, 211 (1998). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1998.0284&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9493407&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.09.04.20188359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071993600007&link_type=ISI) 26. [26]. C. Fraser, D. A. Cummings, D. Klinkenberg, D. S. Burke, and N. M. Ferguson, Influenza transmission in households during the 1918 pandemic, American journal of epidemiology 174, 505 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwr122&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21749971&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.09.04.20188359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000294356800001&link_type=ISI) 27. [27]. E. Bonabeau, L. Toubiana, and A. Flahault, The geographical spread of influenza, Proceedings of the Royal Society of London. Series B: Biological Sciences 265, 2421 (1998). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1998.0593&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9921681&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.09.04.20188359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000077911400012&link_type=ISI) 28. [28]. M. Smith, K. Yourish, S. Almukhtar, K. Collins, D. Ivory, and A. Harmon, Covid-19 data repository, [https://github.com/nytimes/covid-19-data](https://github.com/nytimes/covid-19-data) (2020), accessed: 03-09-2020. ## References 1. [1].1. Soon-Hyung Yook, 2. Hawoong Jeong, and 3. Albert-László Barabási Modeling the Internet’s large-scale topology, Soon-Hyung Yook, Hawoong Jeong, and Albert-László Barabási, PNAS 99 (21) pp. 13382–13386 (2002) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiOTkvMjEvMTMzODIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMC8xNC8yMDIwLjA5LjA0LjIwMTg4MzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 2. [2].1. Tomoya Mori, 2. Tony E. Smith, and 3. Wen-Tai Hsu Common power laws for cities and spatial fractal structures, Tomoya Mori, Tony E. Smith, and Wen-Tai Hsu, PNAS 117 (12) pp. 6469–6475 (2020) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE3LzEyLzY0NjkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMC8xNC8yMDIwLjA5LjA0LjIwMTg4MzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: F3/embed/inline-graphic-4.gif [5]: F3/embed/inline-graphic-5.gif [6]: F3/embed/inline-graphic-6.gif [7]: F4/embed/inline-graphic-7.gif [8]: F7/embed/inline-graphic-8.gif