Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# The universal visitation law of human mobility

## Abstract

Human mobility impacts many aspects of a city, from its spatial structure1,2,3 to its response to an epidemic4,5,6,7. It is also ultimately key to social interactions8, innovation9,10 and productivity11. However, our quantitative understanding of the aggregate movements of individuals remains incomplete. Existing models—such as the gravity law12,13 or the radiation model14—concentrate on the purely spatial dependence of mobility flows and do not capture the varying frequencies of recurrent visits to the same locations. Here we reveal a simple and robust scaling law that captures the temporal and spatial spectrum of population movement on the basis of large-scale mobility data from diverse cities around the globe. According to this law, the number of visitors to any location decreases as the inverse square of the product of their visiting frequency and travel distance. We further show that the spatio-temporal flows to different locations give rise to prominent spatial clusters with an area distribution that follows Zipf’s law15. Finally, we build an individual mobility model based on exploration and preferential return to provide a mechanistic explanation for the discovered scaling law and the emerging spatial structure. Our findings corroborate long-standing conjectures in human geography (such as central place theory16 and Weber’s theory of emergent optimality10) and allow for predictions of recurrent flows, providing a basis for applications in urban planning, traffic engineering and the mitigation of epidemic diseases.

## Access options

Rent or Buy article

Get time limited or full article access on ReadCube.

from\$8.99

All prices are NET prices.

## Data availability

Raw mobility data are not publicly available to preserve privacy. Grid-cell-level data to reproduce the findings of this study can be requested from the corresponding author.

## Code availability

The code to replicate this research can be requested from the corresponding author.

## References

1. 1.

Batty, M. The New Science of Cities (MIT Press, 2013).

2. 2.

Barthelemy, M. The Structure and Dynamics of Cities (Cambridge Univ. Press, 2016).

3. 3.

Louail, T. et al. Uncovering the spatial structure of mobility networks. Nat. Commun. 6, 6007 (2015).

4. 4.

Anderson, R. M. & May, R. M. Infectious Diseases of Humans: Dynamics and Control (Oxford Univ. Press, 1991).

5. 5.

Eubank, S. et al. Modelling disease outbreaks in realistic urban social networks. Nature 429, 180–184 (2004).

6. 6.

Wesolowski, A. et al. Quantifying the impact of human mobility on malaria. Science 338, 267–270 (2012).

7. 7.

Jia, J. S. et al. Population flow drives spatio-temporal distribution of COVID-19 in China. Nature 582, 389–394 (2020).

8. 8.

Bettencourt, L. M. A. The origins of scaling in cities. Science 340, 1438–1441 (2013).

9. 9.

Glaeser, E. L., Kallal, H. D., Scheinkman, J. A. & Shleifer, A. Growth in cities. J. Polit. Econ. 100, 1126–1152 (1992).

10. 10.

Fujita, M., Krugman, P. R. & Venables, A. J. The Spatial Economy: Cities, Regions, and International Trade (MIT Press, 1999).

11. 11.

Sveikauskas, L. The productivity of cities. Q. J. Econ. 89, 393–413 (1975).

12. 12.

Zipf, G. K. The P1P2/D hypothesis: on the intercity movement of persons. Am. Sociol. Rev. 11, 677–686 (1946).

13. 13.

Erlander, S. & Stewart, N. F. The Gravity Model in Transportation Analysis: Theory and Extensions (CRC Press, 1990).

14. 14.

Simini, F., González, M. C., Maritan, A. & Barabási, A.-L. A universal model for mobility and migration patterns. Nature 484, 96–100 (2012).

15. 15.

Zipf, G. K. Human Behavior and the Principle of Least Effort (Addison-Wesley, 1949).

16. 16.

Christaller, W. Die zentralen Orte in Süddeutschland (Gustav Fischer, 1933).

17. 17.

Pan, W., Ghoshal, G., Krumme, C., Cebrian, M. & Pentland, A. Urban characteristics attributable to density-driven tie formation. Nat. Commun. 4, 1961 (2013).

18. 18.

Schläpfer, M. et al. The scaling of human interactions with city size. J. R. Soc. Interface 11, 20130789 (2014).

19. 19.

Li, R. et al. Simple spatial scaling rules behind complex cities. Nat. Commun. 8, 1841 (2017).

20. 20.

Anas, A., Arnott, R. & Small, K. A. Urban spatial structure. J. Econ. Lit. 36, 1426–1464 (1998).

21. 21.

Henderson, V. & Thisse, J.-F. Handbook of Regional and Urban Economics: Cities and Geography (Elsevier, 2004).

22. 22.

Ewing, R. & Hamidi, S. Compactness versus sprawl: a review of recent evidence from the United States. J. Plann. Lit. 30, 413–432 (2015).

23. 23.

Chang, S. et al. Mobility network models of COVID-19 explain inequities and inform reopening. Nature 589, 82–87 (2021).

24. 24.

Axhausen, K. W. & Gärling, T. Activity-based approaches to travel analysis: conceptual frameworks, models, and research problems. Transp. Rev. 12, 323–341 (1992).

25. 25.

Barbosa, H. et al. Human mobility: models and applications. Phys. Rep. 734, 1–74 (2018).

26. 26.

Stouffer, S. A. Intervening opportunities: a theory relating mobility and distance. Am. Sociol. Rev. 5, 845–867 (1940).

27. 27.

Noulas, A., Scellato, S., Lambiotte, R., Pontil, M. & Mascolo, C. A tale of many cities: universal patterns in human urban mobility. PLoS ONE 7, e37027 (2012); correction 7, https://doi.org/10.1371/annotation/ca85bf7a-7922-47d5-8bfb-bcdf25af8c72 (2012).

28. 28.

Yan, X.-Y., Wang, W.-X., Gao, Z.-Y. & Lai, Y.-C. Universal model of individual and population mobility on diverse spatial scales. Nat. Commun. 8, 1639 (2017).

29. 29.

Mazzoli, M. et al. Field theory for recurrent mobility. Nat. Commun. 10, 3895 (2019).

30. 30.

Brockmann, D., Hufnagel, L. & Geisel, T. The scaling laws of human travel. Nature 439, 462–465 (2006).

31. 31.

Gonzalez, M. C., Hidalgo, C. A. & Barabási, A.-L. Understanding individual human mobility patterns. Nature 453, 779–782 (2008); addendum 453, 779–782 (2008).

32. 32.

Song, C., Koren, T., Wang, P. & Barabási, A.-L. Modelling the scaling properties of human mobility. Nat. Phys. 6, 818–823 (2010).

33. 33.

Alessandretti, L., Aslak, U. & Lehmann, S. The scales of human mobility. Nature 587, 402–407 (2020).

34. 34.

Bertaud, A. Order Without Design: How Markets Shape Cities (MIT Press, 2018).

35. 35.

Louail, T. et al. From mobile phone data to the spatial structure of cities. Sci. Rep. 4, 5276 (2014).

36. 36.

Zhong, C. et al. Revealing centrality in the spatial structure of cities from human activity patterns. Urban Stud. 54, 437–455 (2017).

37. 37.

Rozenfeld, H. D. et al. Laws of population growth. Proc. Natl Acad. Sci. USA 105, 18702–18707 (2008).

38. 38.

Cao, W., Dong, L., Wu, L. & Liu, Y. Quantifying urban areas with multi-source data based on percolation theory. Remote Sens. Environ. 241, 111730 (2020).

39. 39.

Batty, M. The size, scale, and shape of cities. Science 319, 769–771 (2008).

40. 40.

Rozenfeld, H. D., Rybski, D., Gabaix, X. & Makse, H. A. The area and population of cities: new insights from a different perspective on cities. Am. Econ. Rev. 101, 2205–2225 (2011).

41. 41.

Pappalardo, L. et al. Returners and explorers dichotomy in human mobility. Nat. Commun. 6, 8166 (2015).

42. 42.

Weber, A. & Friedrich, C. J. Alfred Weber’s Theory of the Location of Industries (Univ. Chicago Press, 1929).

43. 43.

Kreps, D. M. Game Theory and Economic Modelling (Oxford Univ. Press, 1990).

44. 44.

Myerson, R. B. Game Theory (Harvard Univ. Press, 2013).

45. 45.

Blondel, V. D., Decuyper, A. & Krings, G. A survey of results on mobile phone datasets analysis. EPJ Data Sci. 4, 10 (2015).

46. 46.

Calabrese, F., Diao, M., Di Lorenzo, G., Ferreira, J. Jr & Ratti, C. Understanding individual mobility patterns from urban sensing data: a mobile phone trace example. Transp. Res. Part C 26, 301–313 (2013).

47. 47.

de Montjoye, Y.-A., Smoreda, Z., Trinquart, R., Ziemlicki, C. & Blondel, V. D. D4D-Senegal: the second mobile phone data for development challenge. Preprint at https://arxiv.org/abs/1407.4885 (2014).

48. 48.

Blondel, V. D. et al. Data for development: the D4D challenge on mobile phone data. Preprint at https://arxiv.org/abs/1210.0137 (2012).

49. 49.

Jiang, S. et al. The TimeGeo modeling framework for urban mobility without travel surveys. Proc. Natl Acad. Sci. USA 113, E5370–E5378 (2016); correction 113, E7137 (2016).

50. 50.

Song, C., Qu, Z., Blumm, N. & Barabási, A.-L. Limits of predictability in human mobility. Science 327, 1018–1021 (2010).

51. 51.

Wilson, A. G. The use of entropy maximising models in the theory of trip distribution, mode split and route split. J. Transp. Econ. Policy 3, 108–126 (1969).

52. 52.

Masucci, A. P., Serras, J., Johansson, A. & Batty, M. Gravity versus radiation models: on the importance of scale and heterogeneity in commuting flows. Phys. Rev. E 88, 022812 (2013).

Download references

## Acknowledgements

We thank L. M. A. Bettencourt, R. Sinatra, P. Herthogs, G. Du, Y. Qin and Y. Liu for helpful discussions; W. Cao for assistance with performing the CCA analysis; and S. Grauwin for providing the MATLAB code for the radiation model. We acknowledge Airsage, ORANGE/SONATEL and Singtel for providing the data. This work was supported by the National Science Foundation (grant number PHY1838420), the AT&T Foundation, the MIT SMART programme, the MIT CCES programme, Audi Volkswagen, BBVA, Ericsson, Ferrovial, GE, the MIT Senseable City Lab Consortium, the John Templeton Foundation (grant number 15705), the Eugene and Clare Thaw Charitable Trust, Toby Shannan, the Charities Aid Foundation of Canada, the US Army Research Office Minerva Programme (grant number W911NF-12-1-0097), the Singapore National Research Foundation (FI 370074016) and the National Natural Science Foundation of China (grant number 41801299).

## Author information

Authors

### Contributions

G.B.W., C.R., M. Schläpfer, L.D., K.O. and P.S. designed the research. L.D., M. Schläpfer, K.O. and M.V. processed and analysed the data. G.B.W. and M. Schläpfer identified the data collapse to a single, universal curve. K.O., L.D., S.A. and P.S. developed the PEPR model. H.S., M. Schläpfer and G.B.W. developed the theoretical argument. M.V. and P.S. tested the Fermat–Toricelli–Weber metric. L.D., M. Schläpfer, K.O., M. Szell, P.S., G.B.W. and C.R. wrote the paper. C.R. and G.B.W. contributed equally as senior authors. All authors discussed the results and reviewed the manuscript.

### Corresponding author

Correspondence to Lei Dong.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

Peer review information Nature thanks Laura Alessandretti, Marc Barthelemy and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Extended data figures and tables

### Extended Data Fig. 1 The spatio-temporal structure of movement in cities.

Panels show visitor influx maps for Greater Boston for different parameters (rf). The colour of each grid cell (500 m × 500 m) indicates the value of the spectral flow ρ. Remarkably, visitor influx maps for the same quantity v = rf are nearly identical, as is clear from viewing along the diagonals indicated by the coloured arrows in the figure. Hence, doubling the visitation frequency f (from top row to bottom row) results in the same quantitative decrease of the influx as doubling the travel distance r (from left column to right column).

### Extended Data Fig. 2 Empirical power-law exponents of the distance–frequency distribution.

a, Histogram of the exponents for all locations in the Greater Boston area. The values were determined using ordinary least squares minimization to a linear relation of the logarithmically transformed variables. The red line shows η = 2, consistent with our theoretical argument. b, Corresponding histogram of the R2 values.

### Extended Data Fig. 3 Universality of the scaling relation ρ ∝ (rf)−2 across Greater Boston.

The panels depict the data for individual locations (500 m × 500 m grid cells), ranked according to the total number of visitors from neighbouring cells. Shown are locations of rank 1–30 (from top left to bottom right). The geographic coordinates of each location (latitude and longitude of the centre point of the grid cell) are indicated. The straight lines denote the inverse square of rf (slope = −2), consistent with our theoretical argument.

### Extended Data Fig. 4 Universality of the scaling relation ρ ∝ (rf)−2 across Portugal.

The panels depict the data for individual locations (1 km × 1 km grid cells), ranked according to the total number of visitors from neighbouring cells. Shown are locations of rank 1–30 (from top left to bottom right). The geographic coordinates of each location (latitude and longitude of the centre point of the grid cell) are indicated. The straight lines denote the inverse square of rf (slope = −2), consistent with our theoretical argument.

### Extended Data Fig. 5 Universality of the scaling relation ρ ∝ (rf)−2 across Dakar.

The panels depict the data for individual locations (1 km × 1 km grid cells), ranked according to the total number of visitors from neighbouring cells. Shown are locations of rank 1–30 (from top left to bottom right). The geographic coordinates of each location (latitude and longitude of the centre point of the grid cell) are indicated. The straight lines denote the inverse square of rf (slope = −2), consistent with our theoretical argument.

### Extended Data Fig. 6 Universality of the scaling relation ρ ∝ (rf)−2 across Singapore.

The panels depict the data for individual locations (500 m × 500 m grid cells), ranked according to the total number of visitors from neighbouring cells. Shown are locations of rank 1–30 (from top left to bottom right). The geographic coordinates of each location (latitude and longitude of the centre point of the grid cell) are indicated. The straight lines denote the inverse square of rf (slope = −2), consistent with our theoretical argument.

### Extended Data Fig. 7 Simulation results of the EPR model.

a, b, Generated number of visits (a) and attractiveness values μi (b). c, d, The EPR model generates the rf scaling of the population flows with a scaling exponent that is in remarkable agreement with the data. The generated visitor counts, Ni(rf), are shown in c, and the resulting rf scaling of the spectral flows, ρi(rf), is shown in d. The generated attractiveness values μi are rather homogeneous and uniform across space, which is in contrast to the empirical data (b). Model parameters are taken from Song et al.32 (Methods).

### Extended Data Fig. 8 Estimation of the magnitude of flows from population density ρpop.

The schematic shows a zoom-in on the immediate vicinity of a destination location j (small values of r), where it is reasonable to assume that ρpop(j) ≈ constant. Hence, the local population density imposes an upper bound on the influx, ∫ρjdfρpop(j). A simple boundary condition of the continuous model then dictates that the minimum visiting frequency of all individuals living directly on the boundary (each being assigned to a point at r = rj) assumes the minimum frequency with which the individuals living inside the attracting location return home, fmin ≈ fhome. The minimum distance rmin for locations from which individuals visit with minimum frequency fmin < fhome increases with decreasing value of fmin.

### Extended Data Fig. 9 CPT and radius of attraction.

a, Schematics of CPT, showing the spatial arrangement of three tiers of centres (see Supplementary Information for details). This hierarchical arrangement of central places results in the most efficient transport network. bf, Average travel distance per visit rf to perform activities with fixed visiting frequency f across all locations in Greater Boston (b), Singapore (c), Dakar (d), Abidjan (e) and Lisbon (f). We find a clear inverse relation, rf 1/f. The quantity rf can be interpreted as the characteristic distance associated with the level of specialization of the functions provided by the locations.

### Extended Data Fig. 10 Fermat–Torricelli–Weber (FTW) efficiency of collective human movements.

a, The schematic shows how the FTW efficiency is computed (see Supplementary Information). The effective distance travelled by the visitors of a specific location (cell) can be minimized by moving it on the grid. The efficiency is $$\Delta {\mathcal{D}}/{\mathcal{D}}$$, which is the ratio between the reduction of the effective travel distance of all visitors when moving the cell from its actual location to the optimum FTW point ($$\Delta {\mathcal{D}}$$) and the actual effective travel distance of all visitors to that cell ($${\mathcal{D}}$$). b, Density plots representing the number of cells with a given number of visits and FTW efficiency for the Greater Boston area (for the month of August 2009). The FTW efficiency is computed for each cell based on visits made by visitors who live at distances larger than a given threshold value rthr. For rthr = 0 (top left), the density of locations is particularly high where the FTW efficiency is very high. As the number of visits is increased, the distribution becomes narrower and the FTW efficiency increases. This pattern is generally also valid for larger values of rthr but becomes weaker.

## Supplementary information

### Supplementary Information

This file contains supplementary text, supplementary figures 1 – 28, supplementary equations 1 – 10 and supplementary references.

## Rights and permissions

Reprints and Permissions

## About this article

### Cite this article

Schläpfer, M., Dong, L., O’Keeffe, K. et al. The universal visitation law of human mobility. Nature 593, 522–527 (2021). https://doi.org/10.1038/s41586-021-03480-9

Download citation

• Received:

• Accepted:

• Published:

• Issue Date:

## Further reading

• ### Trip frequency is key ingredient in new law of human travel

• Laura Alessandretti
• Sune Lehmann

Nature (2021)

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.