How congestion shapes cities: from mobility patterns to scaling

How congestion shapes cities: from mobility patterns to scaling

Rémi Louf Institut de Physique Théorique, CEA, CNRS-URA 2306, F-91191, Gif-sur-Yvette, France
   Marc Barthelemy Institut de Physique Théorique, CEA, CNRS-URA 2306, F-91191, Gif-sur-Yvette, France

The recent availability of data for cities has allowed scientists to exhibit scalings which present themselves in the form of a power-law dependence on population of various socio-economical and structural indicators. We propose here a stochastic theory of urban growth which accounts for some of the observed scalings and we confirm these predictions on US and OECD empirical data. In particular, we show that the dependence on population size of the total number of miles driven daily, the total length of the road network, the total traffic delay, the total consumption of gasoline, the quantity of emitted and the relation between area and population of cities, are all governed by a single parameter which characterizes the sensitivity to congestion. Our results suggest that diseconomies associated with congestion scale superlinearly with population size, implying that –despite polycentrism– cities whose transportation infrastructure rely heavily on traffic sensitive modes are unsustainable.

I Introduction

The recent availability of an unprecedented amount of data has made possible quantitative studies of cities Fujita:1999 (); Batty:2007 (); Marshall:2004 (), opening the way to a Science of Cities. In particular, the discovery of allometric scaling relationships in cities has driven the quantitative research on urban systems in the past years. Indeed, there is a great amount of evidence to show that different socio-economic indicators in cities, such as the GDP, the crime rate, the number of patents as well as different structural indicators such as the total length of the road network, the urbanized land area, etc., exhibit robust scaling relationships with respect to population Newman:1989 (); Makse:1995 (); Pumain:2006 (); Bettencourt:2007 (); Samaniego:2008 (); Rozenfeld:2008 (); Pan:2013 (). The existence of these simple scaling relationship hints at the existence of universal processes shared by urban systems, and thus at the possibility of modeling cities.

A common trait shared by all complex systems –including cities– is the existence of a large variety of processes occuring over a wide range of time and spatial scales. The main obstacle to the understanding of these systems therefore resides in uncovering the hierarchy of processes and in singling out the few ones which govern their dynamics. Albeit difficult, the hierarchisation of processes is of prime importance. A failure to do so leads to models which are either too complex to give any real insight into the phenomenon, or too simple and abstract to have any resemblance with reality. As a matter of fact, despite numerous attempts Fujita:1982 (); Makse:1995 (); Batty:2008 (); Frasco:2013 (); Bettencourt:2010 (); Bettencourt:2013 (), a theoretical understanding of many observed empirical regularities in cities is still missing.

In the present study, we show that the spatial structure of the mobility pattern controls the behaviour of many quantities in urban systems. Indeed, cities are not only defined by the spatial organisation of places fulfilling different functions –shops, places of residence, workplaces, etc.– but also by the way indivduals move among them. Understanding where people live, where and how they travel within the city thus appears as a necessary step towards a scientific theory of cities.

Although an increasing amount of data about mobility is now available Gonzalez:2008 (), we still lack a simple model explaining the dominant mechanisms governing the formation and evolution of mobility patterns. Many factors such as geographical constraints, facilities location and available transportation –to name a few– can impact the mobility and it thus appears as an intricate issue. Here, we tackle the problem of mobility by making simplifying –yet not simple– assumptions, trying to grasp the most important parameters which define the problem. We thus build upon a simple out-of-equilibrium model previously developed Louf:2013 (). This model, among other things, accounts for the polycentric transition of cities and gives a prediction for the number of centers as a function of population. We show that this framework allows us to predict the behaviour of many quantities related to mobility and the structure of cities: the scaling with population of the total time wasted in congestion, transport related emissions, total travelled distance, total lane miles and surface area.

Our results allow us give a quantitative insight into two important debates around urban systems. First, we are able to discuss the benefits of polycentricity and quantify some of its aspects. Then, maybe more importantly, we are able to put into perspective the sustainability of urban systems.

Ii Results

ii.1 Naive scalings

We start by presenting some naive arguments to estimate the scaling exponents for the area , the total daily distance driven and the total lane miles . Although these predictions turn out to be wrong, naive scalings are useful as a first approach to the problem as they allow us understand how the different quantities relate to one another.

ii.1.1 Surface area

First, we would like to estimate the dependence of the area of a city on its population –a long standing problem in the field Makse:1995 (). A first crude approach would be to assume that cities evolve in such a way that their population density remains constant. This assumption straighforwardly implies that the area should scale linearly with population


where is the average surface occupied by each individual (the assumption of a constant density is then equivalent to the one of a constant average surface per capita).

ii.1.2 Total length of roads

We would now like to estimate the total length of all the roads within a city. If we consider that the network formed by streets is such that all the nodes (intersections) are connected to their closest neighbour, the typical length of a road segment is given by


where is the number of intersections. Previous studies of road networks in different regions, and over extended time periods Strano:2012 (); Barthelemy:2013 (), have shown that the number of intersections is proportional to the population size. Therefore, the typical length of a road segment (between two intersections) varies with the population size as


and the total length of the network should then scale as


Using the naive scaling for the dependence of on population size given previously in Eq. 1 we finally get


ii.1.3 Total daily commuting distance

Individual constraint

We would also like to estimate the total commuting distance . The first constraint on this distance comes from individuals’s limitations and behaviour. We make here the simple assumption that individuals choose their residence and work place such that their total commuting distance is fixed (or at least smaller than a certain value) and equal on average to . In that case, we would simply have


(by constant, we mean independent from the population size of the city).

The city structure constraint

An additional contraint on is given by the structure of the city Samaniego:2008 (). Indeed, the individual commuting distance is also related to the total suface area of the city and the location of activity centers.

If we first assume that the city is monocentric, individuals are all commuting to the same center and the typical commuting distance is controlled by the typical size of the city of order


On the other hand, if we assume that the city is completely decentralized, the typical commuting distance is of order the nearest neighbour distance , and we obtain


ii.2 Comparison of naive scalings with empirical results

The comparison of the naive exponents with the exponents measured on US data is shown in Table 1 (see the Methods section for details about the data). There are important discrepancies, which we discuss in the following.

First, we note that the naive scaling for the surface area predicts a value of the exponent that is quantitatively –and worse, qualitatively– different from that observed. Indeed, we find that for real cities


with . While the naive argument implies a linear dependence of the surface area with population, we find a sublinear scaling in the data, which is a qualitatively different behavior (Table 1). This disagreement on this basic quantity will naturally impact the scaling of the other quantities.

The data also show that can be considered reasonably independent from (with a value of approximately miles for the US, see Fig. 1), in agreement with the individual constraint assumption (Eq. 6). This finding is also in agreement with the results drawn from census data in Germany by Wilkerson:2013 (). Although this assumption of a constant distance is simple and verified on the US data, we think that it deserves to be systematically tested on other datasets for other countries and cities.

Figure 1: Constant daily driven distance per capita. (a) daily total driven distance per capita as a function of population for 441 urbanised area in the US in 2010. The data shown in the plot are compatible with a population-independent behaviour. (b) Histogram of the daily total driven distance per capita for the same cities. The average daily driven distance for these cities is miles, and the standard deviation miles.

Finally, the scaling of given in the extreme cases of a monocentric city structure and a totally decentralized city structure disagree with the value measured on data (see Table 1). This suggests that most cities have a structure that is neither completely centralized, nor totally decentralized. In particular, this result cast some doubts about the study Bettencourt:2013 () which assumes implicitely that cities are always monocentric. Any situation between the two previous extreme cases would give a scaling of the form where . One can easily see that this expression is consistent with that of and if


which is indeed what we observe empirically (up to error bars). This preliminary analysis thus leads us to the conclusion that, in order to compute the various exponents, we need to better describe the structure of commuting patterns. In other words, we need to find a description of cities that goes beyond the naive monocentric or totally decentralized views, and which accounts for the observed sub-linear scaling of the surface area .

ii.3 Beyond naive scaling: modeling mobility patterns

We begin with the assumption that mobility patterns are mostly driven by the daily commuting and we would like to understand how an individual, given his household location, will choose his job location. We assume that this choice will be determined by two dominant factors: the expected wage at a given job, and the commuting time to this job’s location. Indeed, places with high average salaries are attractive, but having to spend a sensible amount of time commuting every day is less desirable. We assume there are potential activity centers in the city, each characterized by an average wage at location . This wage is endogenously determined and depends a priori on many factors such as agglomeration effects Glaeser:2001 (), the type of industry, etc. Although it is in principle possible to write down equations to determine the wage (as attempted in Fujita:1982 () for instance), not only is it impossible to solve them, but also not necessarily useful. A similar situation arises in physics when one studies the behaviour of atoms made of a large number of electrons. Physicists found out Dyson:1962 () that, in fact, a statistical description of these systems relying on random matrices could lead to predictions which agree with experimental results. We would like to import this idea of replacing a complex quantity such as wages –which depends on so many factors and interactions– by a random one in spatial economics. So, we treat the wage as if it was exogenous and random Louf:2013 (), that is we write where represents the typical income in this city and is a random number chosen uniformly in . Furthermore, we assume that the commuting time does not only depend on the distance between the two places, but also on the traffic between those two locations. An individual living at will thus commute to the center which corresponds to the best trade-off between income and commuting time, thus to the center such that the quantity


is maximum Louf:2013 (). The quantity is the euclidean distance between and (both supposed to be scattered randomly across the city), the total incoming traffic at , the capacity of the underlying transportation network, and is an exponent describing the sensitivity of the network to congestion. The quantity is the maximum distance that people can financially travel daily, defined as the ratio between the typical individual income and the transportation costs per unit of distance.

This simple model displays a surprisingly rich behaviour Louf:2013 (). In particular, it accounts for the monocentric to polycentric transition observed in most cities. It has been a well-known fact for quite some time that as cities grow, they evolve from a monocentric organisation where all the activities are concentrated in the same geographical area –usually the central business district– to a more distributed, polycentric organisation Fujita:1982 (); McMillen:2003 (). Several theories in spatial economics exist Fujita:1999 (), but are not satisfactory for many reasons. Among other things, they do not take congestion into account and have no predictive, testable content Bouchaud:2008 (). Within this framework, congestion is actually responsible for the transition, and the number of activity centers in a city of population size is on average given by




Using data of employment per Zip Code Area in the US Louf:2013 (), we showed that


where we measure ( confidence interval (CI)). In other words, the number of centers scales sublinearly with population size.

Figure 2: Mobility and city structure and their impact on agglomeration economies and diseconomies. (a) Variation of the daily total driven distance with the population for 441 urbanized areas in the US in 2010. The dashed line shows the power-law fit with exponent (). (b) Variation of the land area with population for 3540 urbanised areas in the US in 2010. The fit assuming a power-law dependence gives an exponent (). Both exponents are smaller than 1, as predicted by our theory. (c) Variation of the total lane miles with population for 363 urbanised areas in the US. A power law fit (dashed line) gives (). The sublinear behaviour –which agrees with our prediction– means that larger cities need to spend less in infrastructure per capita than smaller ones. (d) Variation of the total delay due to congestion with population for 97 urbanised areas in the US. A power law fit gives an exponent (). The superlinear behaviour agrees with the prediction given by our model and challenges the claims of sustainability of cities.

ii.4 Computing the exponents

ii.4.1 Area

At this stage, the number of centers is a function of population and the area


and we need an additional equation in order to get a closed system. Here we focus on the area and its evolution with the population size, which reflects the growth process of the city. In the following, we will investigate two different approaches. It is worth noting that both approaches give results in qualitative agreement, showing that some stylized facts —such as super- or sublinearity— are very robust.

Fitting procedure.

In the absence of knowledge of the processes responsible for urban sprawl, we can assume that the area behaves as


where is the exponent to be determined, through fits on data. The empirical value for the exponent for the US data is . Once this exponent is given we can then compute the various exponent for the quantities of interest (see the following and table 2). We get for the number of centers


which is sublinear as long as , in agreement with the empirical results for US cities. As we will see, this approach yields the same qualitative behaviours as those predicted with the method of the next section. In other words, even if the main mechanism behind urban sprawl is not congestion, the conclusions of this paper are not affected as long as the area scales sublinearly with population.

Coherent growth.

Let us now assume that the scaling of with population is determined by the number of activity centers and the constant commuting length of individuals. This means that the growth of the area is controlled by the appearance of new activity centers. if we assume that a city is organized around activity centers and that the attraction basin of each of these centers are spatially separated Louf:2013 (), we then have where is the area of each subcenter’s attraction basin. This area is related to the average individual commuting distance by , and we obtain


This leads to expression for the number of centers


which is always smaller than , also in agreement with the empirical results for US cities. We can now also compute the scaling of the surface area


We further assume that is a fraction of the longest possible journey individuals can afford, that is to say


It is important to note that if is independent from , the quantitative predictions of our model would still hold. The final expression for the area is then here given by


where . The exponent is smaller than whatever , which implies that the density of cities increases sublinearly with population. In other words, the density of cities increases with population. We verify this prediction in Table 2, with data about land area of urbanized areas in the US (Figure 2). We find which is not too far from the theoretical value , equal to in this case.

Because the area of an urban system results from centuries of evolution, we do not a priori expect our model –where individual vehicles are assumed to be the only vector of mobility– to give a prediction valid for all countries and all times. Nevertheless, these results give us reasons to believe that the spatial structure of the journey-to-work commuting should still be the dominant factor in the dependence of land area on population.

ii.4.2 Total commuting distance

Using Eq. 6 and Eq. 22 we are now able to compute


We plot for urbanized areas in the US on Figure 2, and one can check in Table 2 that the exponent predicted from the previously measured value of agrees well with the exponent measured on the data. In the fitting case, the exponent would simply be given by .

ii.4.3 Total length of roads

If we use the previously derived expression for the area , we find


The quantity is less than , which implies that scales sublinearly with the city’s population size. In other words, larger cities need less roads per capita than smaller ones: we recover the fact that agglomeration of people in urban centers involves economies of scale for infrastructures. Within the fitting assumption (Eq.16), we would obtain .

ii.4.4 Total delay due to congestion

Unfortunately, agglomeration in cities does not only generate economies. Congestion, for instance, is a major diseconomy associated with the concentration of people in a given area. A simple way to quantify the impairement caused by traffic congestion is through the total delay it generates. If we make the first order approximation that the average free-flow speed is the same for everyone, the total delay due to congestion is given –according to our model– by


If we assume that all the centers share the same number of commuters –a reasonable assumption within our model Louf:2013 ()– we obtain


which, using the expressions for and given in Eq. 23 and Eq. 22 respectively, gives


The total commuting time corresponding to the same distance but without congestion scales as and thus less rapidly than the total delay which scales super-linearly with population (even when polycentricity is taken into account). This means that, for the largest cities, delays due to congestion actually dominate the time spent in traffic, and that economical losses per capita due to the time lost in congestion –and the corresponding strain on people’s life– increase with the size of the city.

In the fitting assumption Eq. 16, and using the same arguments for the calculation of , we easily obtain for the exponent the value .

ii.4.5 Transport related emission. Gasoline consumption

Another diseconomy associated with congestion is the quantity of emitted by cars and the gasoline consumed by motor vehicles. This amount not only depends on the distance that has been driven, but also on the traffic during the journey. It indeed turns out that for the same length driven, a car burns more oil when the traffic is heavy than when the road is clear. Within our model, the presence of traffic is seen in the time spent to cover a given distance, and we write that the quantity of emitted by a vehicle is proportional to the total time spent in traffic, leading to


where is the average quantity of produced per unit time. In the polycentric case with subcenters, the typical trip length is given by and we obtain


The first term in brackets is a constant, and the quantity of is thus dominated by congestion effects at large populations


and the total daily transport-related emission per capita thus scales as


The quantity of emitted per capita in cities thus increases with the size of the city, a consequence of congestion. This prediction agrees with the exponent we measure (Figure 3) on data gathered for US and OECD cities (Data about the area and population of urbanised areas can be found on the Census Bureau website DataUSA1 (), data about congestions in urban areas can be found in the Urban Mobility Report DataUSA2 (), and data about the total lane miles and the daily total miles driven in urbanized areas can be found on on the Federal Highway administration website DataUSA3 ()). We are aware that the scaling of with population size is controversial, with results varying from one study to another. Although a systematic meta-analysis of these results is beyond the scope of this paper, we note that the authors of Fragkias:2013 () are concerned with the total emissions of , while this paper is only concerned with emissions due to transportations. Moreover, our prediction agrees well with the exponent of measured by the authors of Oliveira:2014 () on the same dataset, but with a different definition of cities. Finally, our prediction also agrees with measurements made in Rybski () for developing countries.

Figure 3: Variation of emissions due to transport with city size. In blue, excess (in tons) due to congestion, as given by the Urban Mobility Report (2010) for 101 metropolitan areas in the US DataUSA2 (). In green, we show the estimated emissions (in tons) due to transports, as given by the OECD for 268 metropolitan areas in 28 different countries (Data about the total emissions due to transportation in major metropolitan area in the OECD can be found online OECD ()). The dashed yellow lines represent the least-square fit assuming a power-law dependency with multiplicative noise, which gives respectively for the US data and for the OECD data.

Another important related quantity is the the consumption of gasoline which in principle is proportional to the emission of and the time spent driving. The total daily gasoline consumption is then given by


where is the average quantity of gasoline needed per unit time. From this expression, we see that the total daily gasoline consumption per capita scales as


and is therefore not a simple function of the city density, in contrast with what was suggested by the seminal paper of Newman and Kenworthy Newman:1989 (). At this stage however, more data about gasoline consumption is needed to test this prediction and draw definitive conclusions.

Iii Discussion

iii.1 Monocentric versus polycentric

Although polycentricity emerges naturally from our model as a result of congestion, many circumstances can prevent or foster the appearance of new activity centers in a city. There are many debates as to whether policies should favour polycentric or monocentric developpement of cities. Most of them are based on ideologies and opinions about how cities should be, very few are based on a quantitative understanding of the city as a complex system. Although this only represents a small part of the debate, our model allows to quantify the effect of polycentricity on the total delay due to congestion.

We can indeed compute the total delay due to congestion in the case of a monocentric configuration. In this situation, all the population commutes to a single destination and we have


It follows, using the expression given above for


From the fact that , we indeed find that the total delay due to congestion is worse for monocentric cities than it is for polycentric cities with the same population, which agrees with the usual intuition. More precisely the ratio of delays is given by


where the exponent is of order . Therefore, even though diseconomies associated with polycentric cities scale superlinearly with population, it would be even worse if we did not let cities evolve from the monocentric case. The same reasoning applies to the consumption of gasoline and the emissions. This suggests that, everything else being equal, polycentricity should be favoured for quality of life and environmental reasons.

Figure 4: Scaling down. We consider a population and see how indicators change when we compare it with a system with many cities and the same total population. (a) Variation of the yearly delay per capita due to congestion with the number of cities (normalised by the value corresponding to the single city case). (b) Evolution of the infrastructure length with the number of cities (normalised by the value corresponding to the single city case). Relative gains in terms of commuting time per person decrease faster than infrastructure costs increase, suggesting that life in cities could be improved at a relatively low cost by decentralisation.

iii.2 Megapolis versus urban villages

Also, given the superlinear behaviour of the diseconomies associated with living in cities, it is clear that we would be better off living in several smaller cities rather than a single huge city. However, due to the economies of scale realised in large cities, we can wonder whether this is also economically reasonable. If we assume that the total cost of a city of population is the sum of its infrastructure cost and the economical losses due to congestion we have


where is the average cost of a kilometer of roads, the average hourly wage and the planning horizon in years (this expression is not exhaustive, as the costs dues to emissions and gasoline consumptions are not included). The infrastructure needs maintenance, and its cost depends on the planning horizon as well and can be written where is the construction cost in and the maintenance cost in .

We assume that the population is distributed among cities of the same size (see Figure 4). The total lane miles for the cities reads where is the total lane for one city. The total congestion delay for cities is and we thus obtain the total cost for cities


The number of cities which minimises the total cost is obtained when , leading to (for )


(the actual number of cities is of course an integer, and can be taken as the nearest integer from for instance). It is then economically advantageous to divide the population in several cities if . To illustrate this point, we compute the number of cities which would minimise the cost for a world population . The World Bank estimates the maintenance cost of roads to be of the order of , and the average hourly wage to be of the order of , the value of is taken from the measures on US data, , and . We then obtain


which gives an average city size of . This result is to put in perspective with the fact that the world hosts or so cities with over inhabitants and that this number is still increasing.

iii.3 The most economical population distribution

The previous results assume that we split a large city into many cities of the same size. The cities are however organized in various sizes distributed according to something that can be approximated by a Pareto distribution, as known since Zipf’s work Zipf:1949 (). It is still unclear why we observe such a convergence Batty:2006 (); Cristelli:2012 (). We propose here a new perspective to this debate by asking: Assuming cities are distributed according to a Pareto distribution, what value of the exponent minimises the overall cost? Indeed from above the total cost for a population size is given at large times by


We assume that the population is distributed according to


with and a cut-off population (which is at most equal to the world’s population). The average cost is then given by


leading to


The only consistent solution is obtained for . The dominant term for is given by


The optimal power law distribution minimizes the average cost and is such that . We obtain the following equation


and in the limit we obtain the optimal value for


Numerically, and , leading to . It is interesting to note that this value would lead to a rank-plot exponent () not far from those measured on different countries around the world Soo:2005 (). Although we do not pretend that the above reasoning provides a definitive answer to the Zipf puzzle, it nevertheless suggests that the broad diversity of population might derive from economical considerations, and that there may be a connection between the Zipf law exponent and optimality considerations.

iii.4 Outlook

The superlinear increase of congestion delay with population, and thereby of gasoline consumption and of emissions, has terrible consequences on the economy, the environment, health and well-being. The outlook is nothing short of grim in our ever-urbanising world. As the proportion of human beings living in cities dramatically increases –the UN expects the world population to be urban in 2050 UN:2011 ()– wages are likely to increase Bettencourt:2007 () but not enough to compensate for the negative effects of congestion. As a result, if the individual car stays the dominant transportation mode, cities will put more strain on people’s life, while acting as catalysts for the production of greenhouse gas, responsible for an overall increase of the planet’s temperature Oreskes:2004 (). It is currently believed that advantages associated with living in a large city outweigh the costs. Our results reveal however the existence of very rapidly growing problems such as congestion and emissions, which inevitably begs the question of the sustainability of large cities. It might be time to cut down considerably the use of individual vehicles, or to consider the possibility of living in smaller or medium sized cities: the infrastructure costs () may be larger, but the impact on the environment ( emissions) and on the well-being of people (delays in congestion) would be beneficial (see Figure 3).

The most striking fact about the above results is that despite the apparence of complexity that is conveyed by cities, most of their structure can be explained by the very simple and universal desire for the best achievable balance between income and commuting costs. Our model unifies mobility patterns, spatial structure of cities and allometric scalings in a framework that can be built upon. More work is needed in order to integrate information about firm locations, the influence of public transportation on mobility patterns Roth:2011 (), the effect of the integration of cities into urban systems Rozenblat:2007 (), to understand the fluctuations around the average trends, and to test the validity of the model on different sets of data. We believe however that the results presented here represent a crucial step towards a scientific understanding of cities.

Iv Methods

iv.1 Data

As recently stressed in Arcaute:2013 (), when trying to identify patterns accross cities, one must be careful and consistent in the definition of city boundaries. These authors indeed found out that the scaling exponents measures for several quantities are usually sensitive to the definition chosen for the city. In order to make our results reproducible, we detail in the following the data source and the corresponding city definition.

Total distance driven and lane miles The daily commuting vehicle-miles as well as the total lane miles data were obtained for the year 2011 from the Federal Highway Administration for 441 Urban Areas (as defined by the Census Bureau) in the US.

Area The surface area data were obtained for the year 2010 from the Census Bureau for 3540 Urban Areas (as defined by the Census Bureau) in the US. It is interesting to note that the dependence of the surface area of Metropolitan Statistical Areas with population is a lot less clear-cut, implying that, with respect to surface area, the definition of UAs delineates more coherent systems than the definition of MSAs.

Values of In order to compute a value for , we use for the average wage at the county level, provided by the Bureau of Labor Statistics. For , the transportation cost per unit distance, we use the average gas price per state as given by the U.S. Energy Information Administration, and assume that all vehicles burn the same quantity of gas per unit distance on average. Interestingly, while we have assumed a constant throughout this paper, we have noticed that its effect on the different scalings was not negligible (Compare the results for and between Table 1 and Table 2 for instance), implying that has a small, yet non-zero dependence on the population. This probably comes from the dependence of the average wage on population Bettencourt:2007 (). We leave the investigation of this dependence for further studies.

Total delay and emissions The excess and the total delay due to traffic congestion were obtained for the year 2012 from the Urban Mobility Report for 97 Urban Areas in the US. Also, the quantity of emissions due to transportation was obtained from the OECD for 268 metropolitan areas accross 28 countries for the year 2008. It is worth noting here that the US definition of Urban Area and the OECD definition of Metropolitan Area are qualitatively different, added to the fact that OECD data cover many different countries. Yet, the measured values of the exponent are compatible with each other.

As far as the United States are concerned, we present results for Urban Areas only. Indeed, when data were available for both MSA and Urban Area, we found out that the MSA data did not exhibit as clear-cut regularities as the Urban Area data did. We believe that this effect is due to the lack of a unique, quantitative definition of a city which makes In this work, we assumed that Urban Areas designate areas which are coherent with respect to the quantities we are measuring, and leave the crucial issue of city definition for further studies.

V Acknowledgments

We thank Giulia Carra, Riccardo Gallotti and Thomas Louail for stimulating discussions. MB acknowledges funding from the EU commission through project EUNOIA (FP7-DG.Connect-318367).

Vi Author Contributions Statement

RL, MB designed, performed research and wrote the paper.

Vii Competing financial interests

The author(s) declare no competing financial interests.


  • (1) Fujita, M., Krugman, P. R., & Venables, A. J. The spatial economy: cities, regions and international trade (The MIT press, Cambridge, USA, 1999).
  • (2) Batty, M. Cities and complexity: understanding cities with cellular automata, agent-based models, and fractals (The MIT press, Cambridge, USA 2007)
  • (3) Marshall, S. Streets and patterns (Routledge, Abingdon, 2004)
  • (4) Newman, P. W. & Kenworthy, J. R. Gasoline consumption and cities: a comparison of US cities with a global survey. Journal of the American Planning Association 55, 24–37 (1989).
  • (5) Makse, H.A., Havlin, S., Stanley, H.E. Modelling urban growth. Nature 377, 608-612 (1995).
  • (6) Pumain, D., Paulus, F., Vacchiani-Marcuzzo, C., & Lobo, J. An evolutionary theory for interpreting urban scaling laws. Cybergeo: European Journal of Geography. (2006)
  • (7) Bettencourt, L.M.A., Lobo, J., Helbing, D., Kühnert, C. & West, G.B. Growth, innovation, scaling, and the pace of life in cities. Proc. Natl. Acad. Sci. U.S.A. 104, 7301-7306 (2007).
  • (8) Samaniego, H., & Moses, M.E., Cities as organisms: Allometric scaling of urban road networks Journal of Transport and Land Use 1, 21-39 (2008).
  • (9) Rozenfeld, H. D. et al. Laws of population growth. Proc. Natl Acad. Sci. USA 105, 18702–18707 (2008).
  • (10) Pan, W., Ghoshal, G., Krumme, C., Cebrian, M., & Pentland, A. Urban characteristics attributable to density-driven tie formation. Nature communications, 4:1961 (2013).
  • (11) Fujita, M. & Ogawa, H. Multiple equilibria and structural transition of non-monocentric urban configurations. Regional science and urban economics 12, 161-196 (1982).
  • (12) Batty, M. The Size, Scale, and Shape of Cities. Science 319, 769–771 (2008).
  • (13) Frasco, G.F., Sun, J., Rozenfeld, H.D. & ben-Avraham, D. Spatially distributed social complex networks ArXiv:1306.0257 (2013)
  • (14) Bettencourt, L. & West, G. A unified theory of urban living. Nature 467, 912–913 (2010).
  • (15) Bettencourt, L. M. A. The Origins of Scaling in Cities. Science 340, 1438–1441 (2013).
  • (16) Gonzalez, M., Hidalgo, C.A., Barabasi, A.-L.. Understanding individual human mobility patterns. Nature 453, 779-782 (2008).
  • (17) Louf, R. & Barthelemy, M. Modeling the polycentric transition of cities. Physical Review Letters 111, 198702 (2013)
  • (18) Strano, E., Nicosia, V., Latora, V., Porta, S. & Barthelemy, M. Elementary processes governing the evolution of road networks Scientific Reports 2:296 (2012).
  • (19) Barthelemy, M., Bordin, P., Berestycki, H., & Gribaudi, M. Self-organization versus top-down planning in the evolution of a city. Scientific Reports 3:2153 (2013).
  • (20) Wilkerson, G., Khalili, R. & Schmid, S. Urban Mobility Scaling: Lessons from Little Data. arXiv:1401.0207 (2013).
  • (21) Glaeser, E. L. & David C. M. Cities and skills. Journal of Labor Economics 19, 316–342 (2001).
  • (22) F.J. Dyson, Journal of Mathematical Physics 3, 140-156 (1962).
  • (23) McMillen, D. P. & Smith, S. C. The number of subcenters in large urban areas. Journal of Urban Economics 53, 321–338 (2003).
  • (24) Bouchaud, J.-P. Economics needs a scientific revolution. Nature 455, 1181-1181 (2008).
  • (25) Barthelemy, M. Spatial networks. Physics Reports 499, 1–101 (2011).
  • (26) Oliveira, E. A., Andrade Jr., J. S. & Makse, H. A. Large cities are less green. Scientific Reports 4 4235 (2014)
  • (27) Fragkias, M., Lobo, J., Strumsky, D. & Seto, K. C. Does Size Matter? Scaling of CO2 Emissions and U.S. Urban Areas. PLoS ONE 8, e64727 (2013).
  • (28) Rybski, D., Sterzel, T., Reusser, D.E., Winz, A.-L., Fichtner, C., Kropp, J.P. Cities as nuclei of sustainability? arXiv:1304.4406 (2013).
  • (29) Zipf, G.K. Human behavior and the principle of least effort (Addison-Wesley Press, Cambridge, 1949).
  • (30) Batty, M. Rank clocks. Nature 444, 592–596 (2006).
  • (31) Cristelli, M., Batty, M. & Pietronero, L. There is More than a Power Law in Zipf. Scientific Reports 2, 812 (2012).
  • (32) Soo, K. T. Zipf’s Law for cities: a cross-country investigation. Regional science and urban Economics 35, 239–263 (2005).
  • (33) United Nations, The world urban propects: The 2011 revision. (2011)
  • (34) Oreskes, N. The scientific consensus on climate change. Science 306, 1686 (2004)
  • (35) Roth, C., Kang, S. M., Batty, M. & Barthelemy, M. Structure of Urban Movements: Polycentric Activity and Entangled Hierarchical Flows. PLoS ONE 6, e15923 (2011).
  • (36) Rozenblat, C., & Pumain, D. “Firm linkages, innovation and the evolution of urban systems” in Cities in globalization: practices, policies and theories, [130-156] (Routledge, London, 2007).
  • (37) Arcaute, E. et al. City boundaries and the universality of scaling laws. arXiv:1301.1674 (2013)
  • (38) Census bureau: (date of access:01/04/2014)
  • (39) Urban mobility report: (date of access: 01/04/2014)
  • (40) Federal Highway administration: (date of access: 01/04/2014).
  • (41) (Date of access: 01/04/2014).

Appendix A Tables

Quantity Naive exponent Measured value
Table 1: This table displays the value of the exponent governing the behavior with the population obtained by naive arguments and the value obtained from empirical data. The discrepancies reveal the failure of the naive scaling arguments and the necessity to go further and model mobility patterns. The data used for this table can be found in DataUSA1 (); DataUSA2 (); DataUSA3 ().
Quantity Theoretical dependence on Predicted value of the exponents Measured value
in self-consistant case self-consistant case fitting case
 DataUSA1 (); DataUSA2 (); DataUSA3 ()
 DataUSA1 (); DataUSA2 (); DataUSA3 ()
 DataUSA1 (); DataUSA2 (); DataUSA3 ()
 DataUSA1 (); DataUSA2 (); DataUSA3 ()
 OECD ()
 Oliveira:2014 ()
 DataUSA1 (); DataUSA2 (); DataUSA3 ()
 DataUSA1 (); DataUSA2 (); DataUSA3 ()
Table 2: This table displays the predicted theoretical behavior and the empirical observations versus the population size for different quantities: is the daily total driven distance, is the area of the city, is the total length of the road network, is the daily total delay due to congestion, is the yearly total consumption of gasoline and is the total emissions emitted yearly due to transportation. In the third column, we show the predicted values of the exponent of using the value of measured on US employment data Louf:2013 (), and in the fourth column, the value of the exponents directly measured on data about US and OECD cities. The measured values are in good agreement with the prediction. In particular, the exponents for and are consistent with our prediction that their difference should be .
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description