Arxiv:2101.03073V1 [Physics.Soc-Ph] 8 Jan 2021 Physical Measurements Such As Weight Or Size of Objects Are Always Confined to Specific Scales

Home , Abiriba, Abu Tig, Africans in Guangzhou, Cheboksary, Dashiqiao, Denov

Association between population distribution and urban GDP scaling

Haroldo V. Ribeiro1,*, Milena Oehlers2, Ana I. Moreno-Monroy3, J¨urgenP. Kropp2,4, Diego Rybski2,5,†,

1 Departamento de F´ısica,Universidade Estadual de Maringá– Maringá,PR 87020-900, Brazil 2 Potsdam Institute for Climate Impact Research – PIK, Member of Leibniz Association, P.O. Box 601203, 14412 Potsdam, Germany 3 OECD, 1 Rue AndréPascal, 75016, Paris, France; Department of Geography and Planning, University of Liverpool, Chatham St, Liverpool L69 7ZT, United Kingdom 4 Institute for Environmental Science and Geography, University of Potsdam, 14476 Potsdam, Germany 5 Department of Environmental Science Policy and Management, University of California Berkeley, 130 Mulford Hall #3114, Berkeley, CA 94720, USA

* [email protected], † [email protected]

Abstract

Urban scaling and Zipf’s law are two fundamental paradigms for the science of cities. These laws have mostly been investigated independently and are often perceived as disassociated matters. Here we present a large scale investigation about the connection between these two laws using population and GDP data from almost five thousand consistently-defined cities in 96 countries. We empirically demonstrate that both laws are tied to each other and derive an expression relating the urban scaling and Zipf exponents. This expression captures the average tendency of the empirical relation between both exponents, and simulations yield very similar results to the real data after accounting for random variations. We find that while the vast majority of countries exhibit increasing returns to scale of urban GDP, this effect is less pronounced in countries with fewer small cities and more metropolises (small Zipf exponent) than in countries with a more uneven number of small and large cities (large Zipf exponent). Our research puts forward the idea that urban scaling does not solely emerge from intra-city processes, as population distribution and scaling of urban GDP are correlated to each other.

Introduction

arXiv:2101.03073v1 [physics.soc-ph] 8 Jan 2021 Physical measurements such as weight or size of objects are always confined to specific scales. However, the outcomes of several natural phenomena and socio-economic processes can extend across multiple orders of magnitude [1]. Examples of such systems are as diverse as earthquakes [2], stock market fluctuations [3], casualties in human insurgencies [4], and fracture of materials [5]. These systems usually share non-trivial statistical regularities manifested in the form of power-law distributions or power-law relations, the so-called scaling laws. Cities are among these systems as they occur in sizes from thousands to tens of million inhabitants. Urban systems are also well known

January 11, 2021 1/17 to follow scaling relations in time and space, as summarized by Batty [6] in the “seven laws of urban scaling”. Zipf’s law [7, 8] and Bettencourt-West law [9] are two of the best-known scaling laws emerging in urban systems. Initially observed by Auerbach [8] and then popularized by Zipf [7], Zipf’s law for cities states that the distribution of city populations, for a given region or country, is approximated by an inverse power-law function (or the rank-size rule), implying that there are plenty of small cities and very few metropolises. Bettencourt-West law [9], better known as “urban scaling”, establishes a power-law relation between urban indicators and city population. Urban scaling is well illustrated by the concept of agglomeration economies in cities [10, 11], in which the power-law association between urban wealth and population implies that urban wealth increases more than proportionally with city population. Fostered by a recent deluge of highly-detailed city data, researchers from diverse disciplines such as geography, economics, or physics have devoted ongoing efforts in identifying and understanding fundamental principles and regularities underlying urban systems [1,6, 12–14]. While most of these works are either concerned with Zipf’s law [15–25] or urban scaling [26–34], very few have tackled the relationship between both. Zipf’s law and urban scaling have mostly been studied independently because it is commonly assumed that both laws are independent descriptions of urban systems across countries. The few investigations on possible connections between Zipf’s law and urban scaling have a more local character and explore this association within countries [35, 36]. The work of Gomez-Lievano et al. [35] is seminal in this regard and has shown that (under certain conditions) urban scaling can be related to Zipf’s law when urban indicators are also power-law distributed. This connection is obtained from probability distributions and shows that the urban scaling exponent has an upper limit determined by both the Zipf exponent and the power-law exponent that characterizes the urban indicator distribution. This result however does not imply that these three exponents are directly associated with each other. In fact, random permutations of population and urban indicator values only have the effect of removing a possible correlation between the two variables, but do not have the effect of changing their distribution. To date, there has been no attempt to empirically investigate the association between urban scaling and Zipf’s law across a large number of countries. The paucity of such studies certainly reflects the lack of consistent and comparable data across countries. A convincing comparison between these two laws demands a unified city definition across the globe as well as measures of population and urban indicators based on this same definition. And it was not until very recently that satisfying this requirement has become possible. Here we use a new harmonized city definition for investigating Zipf’s law and urban scaling over almost five thousand cities in 96 countries. By analyzing population and gross domestic product (GDP) of these cities per country, we estimate the Zipf and the urban scaling exponents to probe for a possible relation between these two scaling laws. Our empirical results show that both exponents are indeed related to each other and that a functional form of this association can be exactly derived from scaling relations emerging at the country level. We demonstrate that Zipf’s law and urban scaling imply a power-law relation between total urban population and total urban GDP of countries, where the country scaling exponent is dependent on the Zipf and the urban scaling exponents. Because country-aggregated values of urban population and GDP are fixed, there is only one country scaling exponent for total GDP, which in turn associates the urban scaling exponent to the Zipf exponent. We verify the integrity of our model by estimating the country-scaling exponents from the empirical relationship between the Zipf and the urban scaling exponents, and also by showing that numerical simulations yield results very similar to those obtained from real-world data. The connection between both

January 11, 2021 2/17 exponents shows that urban scaling does not only emerge from processes occurring within the city boundaries; instead, it suggests that the population distribution of an urban system does affect urban scaling or vice versa. For the particular case of urban GDP, while almost all countries exhibit increasing returns to scale (urban scaling exponent greater than one), our findings indicate this effect is smaller in countries with a more balanced number of small and large cities (small Zipf exponent) than in countries with a more unbalanced number of small and large cities (large Zipf exponent).

Results Empirical connection between both scaling laws We start by revisiting Zipf’s law and the urban scaling law. Zipf’s law for cities [7] establishes that the rank r and the city population s of an urban system are related through a power-law function r ∼ s−α , (1) where α > 0 is the Zipf exponent. The estimated α varies across countries and epochs, but α ≈ 1 is typically found in empirical studies [18, 23]. Urban scaling laws [9] commonly refer to power-law associations between an urban indicator y and the city population s within an urban system, that is,

y = c sβ , (2)

where c > 0 is a constant and β is the urban scaling exponent. The value of β depends on the type of urban indicator, but increasing returns to scale (β > 1) are usually reported for socio-economic indicators, decreasing returns to scale (β < 1) for infrastructure indicators, and constant returns to scale (β = 1) for indicators related to individual needs [9]. The values of α and β are also susceptible to different city definitions [23, 37], and we thus need a unified definition across different urban systems to test a possible association between their values. To do so, we use a generalized definition of functional urban areas (FUA) recently proposed by Moreno-Monroy et al. [38, 39]. The concept of FUA was initially developed for countries of the OECD and Europe as a unified definition of metropolitan areas, consisting of high-density urban cores and their surrounding areas of influence or commuting areas. The generalized definition we use (so-called eFUA or GHSL-FUA) represents an extension of this concept to countries of the entire globe. By considering the areas of influence of urban cores, eFUAs give a less fragmented representation of the city size distribution because dense clusters proximate to urban cores are rightly not considered as independent cities. While eFUAs delineate world-wide comparable city boundaries, the majority of urban indicators are available at the level of local administrative units, that in addition to not being centralized into a global data set, may change from country to country. We avoid this problem by considering a global gridded data set for GDP [40] (see Materials and Methods). Thus, by combining these two data sources (see Materials and Methods and Fig. 1 in S1 Appendix), we create a consistent data set comprising the population (s) and the GDP (y) of 4,571 cities from 96 countries. Having this harmonized data set, we estimate the Zipf exponent α and urban scaling exponent β after grouping the values of s and y by country (see Materials and Methods). The scatter plot in Fig. 1A depicts the values of β versus α, where the insets show Zipf’s law and urban scaling of city GDP for three countries. As these three examples illustrate, the data follows Zipf’s law and urban scaling with deviations comparable with other studies about these two laws (see S1 Appendix, Sec. 4 for all countries). The

January 11, 2021 3/17 A High-GDP Mid-GDP Low-GDP B

2.0 =

1.5 Zipf exponent,

1.0 1.2 1.4 1.6

1.0 = 1 C

Urban scaling exponent, = 1.06±0.01 = 0.26±0.05 = 0.79±0.03 0.5 1.22

0.5 1.0 1.5 2.0 2.5 3.0 Zipf exponent, Urban scaling exponent,

1.0 1.5 2.0 2.5

Australia Turkey Russia 2 11.0 1.0 11 1 11 1 10.5 0.5 10 = 0.72±0.17 = 1.00±0.02 = 1.06±0.12 = 1.17±0.06 0 = 1.60±0.24 = 1.27±0.08 10 10.0 0 9 r log City rank, r log City rank, r log City rank, y log City GDP, y log City GDP, 0.0 y log City GDP, 5.5 6.0 6.5 5.5 6.0 6.5 5 6 7 5 6 7 6.0 6.5 7.0 6.0 6.5 7.0 City population, logs City population, logs City population, logs City population, logs City population, logs City population, logs Fig 1. Association between the urban scaling exponent β and the Zipf exponent α. (A) Values of β versus α for each country in our data set. The three different markers distinguish countries according to the tercile values of the total urban GDP distribution (for instance, high-GDP countries have highest ≈ 33% GDP values). The horizontal dashed line shows the β = 1 while the inclined dashed line represents the β = α relationship. As we shall discuss in the next section, the continuous line shows the model of Eq. (9) adjusted to data and the gray shaded region stands for the 95% confidence band. The colored background represents the different intervals of α defined in Eq. (9). The insets (indicated by blue arrows) illustrate Zipf’s law for city population (left) and the urban scaling relationship of city GDP (right) for three countries (see S1 Appendix, Sec. 4 for all considered countries). We have verified that the model of Eq. (9) provides a significantly better description of data when compared with a null model where β is independent of α (intercept-only model: β = constant, see Fig. 4 in S1 Appendix), and that statistical correlations between β and α are significantly enhanced in the (orange shaded) region where our model predicts a linear correspondence between the exponents (Fig. 5 in S1 Appendix). Maps in panels (B) and (C) show the color-coded values of α and β for each country (light gray indicates missing values).

world maps in Figs. 1B and C indicate the regional distribution of the values of α and β. In line with the meta-analysis of Refs. [18, 23], we find Zipf exponents α roughly distributed around 1 with average value and standard deviation equal to 1.24 and 0.38, respectively (Fig. 2A in S1 Appendix). In turn, ≈ 94% of the countries exhibit urban scaling exponents with GDP (β) larger than 1 (Fig. 2B in S1 Appendix), and the average value and standard deviation of β are 1.16 and 0.17, respectively. This result agrees with the idea that economic indicators display increasing returns to scale [9, 41]. The results shown in Fig. 1 suggest a positive association between the values of β and α, that is, an increase of the Zipf exponent α appears to come along with a rise in the urban scaling exponent β (or vice versa). This behavior is also perceptible in the world maps, where we note that regional differences in the values of α (Fig. 1B) are similar to the ones observed for β (Fig. 1C). However, these visual similarities are mainly induced by the largest countries in land area, and a more careful analysis of these maps reveals important differences, especially in the African continent and

January 11, 2021 4/17 Central Asia. These differences are more evident in the scatter plot of Fig. 1A, where we observe a considerable spread in the association between both exponents. We have further quantified the association between the values of β and α by estimating the Spearman rank correlation and the Pearson correlation within a sliding window of size ∆α0 centered in α0. For different values of ∆α0, we find that the correlations peak around α0 ≈ 1.1 and start to decrease and become not statistically significant as the sliding window moves away from this value of α0 (Fig. 3 in S1 Appendix). This analysis thus suggests that the overall association between both exponents is non-linear and that β may approach constant values for large and small values of α.

Country scaling and the association between both exponents To better describe the empirical connection between the urban scaling exponent and the Zipf exponent, we hypothesize that this association can be derived from the relationship between country-aggregated values of urban population (S) and urban GDP (Y ). As we shall show, the combination of Zipf’s and the urban scaling laws implies a power-law relation between total urban population S and total urban GDP Y . This country scaling law is characterized by an exponent γ that is a function of β and α, which in turn yields a mathematical expression for the relation between β and α. To do so, we start by noticing that Zipf’s law implies a power-law dependence for the absolute −(α+1) number of cities with population s, that is, p(s) = ks for smin < s < smax, where smin and smax represent a lower and an upper cutoff associated to the smallest and largest city population of a particular urban system, and k is a normalization constant (see S1 Appendix, Sec. 1 for details). By combining this frequency function p(s) with urban scaling (Eq. 2), we can write the total urban GDP of a country as

s Xmax Z smax Y = csβp(s) ≈ kc sβ−α−1ds , (3) s s=smin min

where the normalization constant k is determined by

s Xmax Z smax S = sp(s) ≈ k s−αds . (4) s s=smin min

To solve these equations, we consider that smin and smax are power-law functions of S

δ smin = aS , (5) θ smax = bS , (6)

where a > 0 and b > 0 are constants, and δ > 0 and θ > 0 are country scaling exponents. Thus, by plugging Eqs. (5) and (6) into Eq. (3), we find

c(α − 1)S aβbαSδβ+θα − aαbβSθβ+δα Y = , (7) α − β abαSδ+θα − aαbSθ+δα

that in the limit of S 1 yields γ Y = Y0S , (8)

where Y0 is a constant and γ = γ(α, β, δ, θ) is the country scaling exponent for total urban GDP. As detailed in S1 Appendix, Sec. 1, the exact form of γ depends on conditions imposed on the exponents α, β, δ, and θ (cases A.1-A.8 in S1 Appendix) which in turn emerge from determining the dominant term of Eq. (7) in the limit of S 1. Therefore, the combination of Zipf’s law (Eq. 1) with the urban scaling (Eq. 2)

January 11, 2021 5/17 and the country scaling relations of smin and smax (Eqs. 5 and 6) leads us to the country scaling of total urban GDP (Eq. 8), where the exponent γ depends on α, β, δ and θ. The country scaling relations of GDP and smax have previously been empirically observed in Refs. [42–44]. We have also verified that Eqs. (8), (5), and (6) hold well for the country-aggregated values in our data set with γ > 1 and δ < θ (Fig. 6 in S1 Appendix). More importantly, while α and β are intra-country exponents having different values for each country, γ, δ and θ are inter-country exponents and there is only one value for each across all countries. Thus, we can solve γ = γ(α, β, δ, θ) for β (S1 Appendix, Sec. 1 for details) to find  1 + γ−1 0 < α ≤ 1  θ γ+δ−1 δ γ−1 β = θ + 1 − θ α 1 < α < 1 + δ , (9)  γ−1 γ−1 1 + δ α ≥ 1 + δ for γ > 1 and δ < θ. This piecewise relationship implies that β is a constant up to α = 1, where it starts to increase linearly until reaching the diagonal line given by β = α, and then continues as another constant. It is worth mentioning that by writing β as a function of α, we are not assuming any causal direction for the association between both exponents. Indeed, we could also solve γ = γ(α, β, δ, θ) for α but this yields a non-functional dependency, that is, the horizontal plateaus defined by Eq. (9) become vertical lines when writing α as function of β. We have adjusted Eq. (9) to the empirical relation between α and β (see Materials and Methods) and the best fitting parameters are γ = 1.06 ± 0.01, δ = 0.26 ± 0.05, and θ = 0.79 ± 0.03. The solid line in Fig. 1A represents the best fit of Eq. (9) and the colored background indicates the different ranges of α defined in Eq. (9). Despite the large variations in the data, our model captures the average tendency of the empirical relation between β and α. We have verified that the statistical correlations between the values of β and α are only statistically significant in the mid-range of α values, while there are no significant correlations within the lower and higher ranges of α values (Fig. 5 in S1 Appendix). Furthermore, the Akaike and Bayesian information criteria indicate that it is at least 100 times more likely that the empirical data come from our model (Eq. 9) than from a null-model assuming no relationship between both exponents (intercept-only model: β = constant, see Fig. 4 in S1 Appendix). We also find that the estimated parameters from Eq. (9) are quite robust against thresholds for the total GDP; the adjusted values of γ, δ, and δ barely change, even if we only consider the countries with top 50% GDP values (Fig. 7 in S1 Appendix). In addition to the previous model validation, the adjustment of Eq. (9) allows us to verify the consistency of our modeling approach through the country scaling relationships. Specifically, the country scaling exponents estimated from Eq. (9) should describe well the empirical country scaling relationships. To do so, we have adjusted the country scaling relationships of Eqs. (5), (6), and (8) by considering only the prefactors (a, b, and Y0) as fitting parameters and fixing country scaling exponents (δ, θ and γ) to the values estimated from Eq. (9). Figures 2A-C show that the three country scaling relations describe quite well the empirical data. Furthermore, the country scaling exponents obtained by fitting Eq. (9) to the empirical relationship between β and α agree well with the estimates directly obtained from the country scaling relations, that is, by fitting Eqs. (5), (6), and (8) to the country-aggregated data, as shown in Fig. 2D. The values of δ and θ estimated from both approaches are not significantly different, while the values of γ are both above one but fitting the country scaling directly from data yields a slightly larger value. We believe this result provides support for our model, especially when considering the level of observed variation in the relationship between β and α as well as that in the country scaling relations.

January 11, 2021 6/17 n x i

A B a C D Estimated from the empirical m m 13 = 1.06 s 6.0 = 0.26 8.0 = 0.79 relation between and s

log Y0 = 3.53±0.06 g log a = 3.18±0.03 log b = 1.02±0.02 g 1.25

o Estimated from the country l o Y

l scaling relations , 7.5 , n g 12 n o o

i 1.00 o l t i 5.5 t , a l a P l 7.0 u u D 11 p 0.75 p o G o

p n

y 6.5 a t

y 0.50 i

5.0 t b i c r

10 c t

U t s s

e 6.0 0.25 l e l g a Country scaling exponents 9 Pearson's R = 0.78 Pearson's R = 0.40 r Pearson's R = 0.91 a m 4.5 0.00 L S 6.5 7.0 7.5 8.0 8.5 6.5 7.0 7.5 8.0 8.5 6.5 7.0 7.5 8.0 8.5 Country's urban population, log S Country's urban population, log S Country's urban population, log S Fig 2. Country scaling relationships. (A) Scaling relation between total urban GDP (Y ) and total urban population (S) for all countries in our data set. (B) Scaling law between the smallest city population of a country (smin) and the country’s urban population (S). (C) Scaling law between the largest city population of a country (smax) and the country’s urban population (S). In these three panels, markers represent the values for each country and the dashed lines are the country scaling relationships of Eqs. (8), (5), and (6), where the exponents γ, δ, and θ are obtained from fitting the model of Eq. (9) to the empirical relation between β and α. Only the constants Y0, a and b are adjusted to data and their best fitting values are shown in the panels (± standard errors). (D) Comparison between estimates of the country scaling exponent obtained by fitting Eq. (9) to the empirical association between β and α (bars in dark colors) and by directly fitting the values of the country scaling relationships (bars in light colors, see Fig. 6 in S1 Appendix for the adjusted scaling laws). Error bars represent 95% bootstrap confidence intervals of the parameters. We notice that both approaches yield similar estimates, which are statistically indistinguishable in the cases of δ and θ.

Simulating the association between both exponents The level of variation in the relation between the exponents β and α (as well as in the country scaling relationships) is significantly high and hampers a more visual comparison with our model. In addition to indicating that our description of the association between both exponents is far from perfect, these variations also reflect possible estimation errors in the population and GDP values, errors related to the definition of city boundaries, and mainly errors associated with estimating the Zipf and urban scaling exponents (as data does not perfectly follow the Zipf’s and urban scaling laws). As these errors emerge from different sources and do not seem to affect the association between β and α systematically, we have treated them as random variations and investigated their role via an in silico experiment. As summarized in the Materials and Methods, we have generated artificial data at the city level (city population s and city GDP y values) by considering Zipf’s law (Eq. 1), urban scaling (Eq. 2), country scaling relations (Eqs. 5 and 6), and our model (Eq. 9). These simulations take as inputs the real values of total urban population S and the exponent α of each country to generate replicas of these urban systems (set of values of s and y), from which we can investigate the relationship between β and α under different levels of variation in the urban and country scaling relations. To do so, we have fixed the variation intensity (that is, the standard deviation of Gaussian random variation around the scaling laws) in the simulated country scaling relations to a level comparable with the empirical data, and varied the variation intensity in the simulated urban scaling σy in percentages of the value observed in the real data (σy = 100% means the standard deviation of the simulation is equal to the standard deviation in the empirical values of β). The first four panels of Fig. 3A show examples of simulated relationships (red circles) in comparison with empirical data (blue dots) for different values of σy. As expected, the simulated relationship perfectly agrees with our model (Eq. 9) when there is no random variation in the urban scaling. More importantly, we note that the scattering of simulated data becomes visually very similar to the empirical data as the intensity of the random variation increases up to σy = 100% (Fig. 8 in S1 Appendix). We can also use the simulated data to corroborate our numerical experiment by verifying the

January 11, 2021 7/17 country scaling relations. Figures 3B-D show the simulated country scaling relations in comparison with the behavior of Eqs. (5), (6), and (8) with parameters estimated from real data. We note that simulated scaling relations for total urban GDP (Y ) and smallest city population (smin) follow very closely the adjusted behavior of the empirical data. However, the simulated results for the largest city population (smax) systematically underestimate the trend observed in the empirical data. This happens because, in our simulations, we have used a random number generator associated with a truncated power-law distribution (between smin and smax) for mimicking population values according to Zipf’s law. Since large populations are rare and the number of simulated cities is finite, the simulated values for the largest city population do not get close enough to the upper bound imposed by the truncated power-law distribution. Consequently, the simulated values of smax underestimate the empirical ones. We solve this issue by replacing the truncated power-law behavior by a power-law distribution with exponential cutoff (Materials and Methods). This modification does not alter the country scaling relations of total urban GDP and smallest city population, that is, the results of Fig. 3B and 3C are not affected by the introduction of the exponential cutoff. Similarly, this modification does not affect the relationship between β and α and the simulated associations with exponential cutoff are very similar to those obtained with the truncated power-law (as indicated by the right-most plot in Fig. 3A). Indeed, the inclusion of this exponential cutoff only modifies the country scaling relation of the largest city population by increasing the simulated values of smax. Figure 3E shows an example of simulated results after after replacing the truncated cut-off with a exponential one for the country scaling relation between smax and S. We observe that the simulated values of smax obtained with the exponential cutoff are closer to the empirical data behavior than those obtained with the truncated power-law distribution (Fig. 3D).

Discussion

We have shown that the combination of Zipf’s law and urban scaling implies a country scaling relationship, where the exponent is a function of the Zipf and the urban scaling exponents. While the Zipf and the urban scaling exponents vary from country to country, there is only one country scaling exponent for a given indicator, which in turn implies a direct association between the urban scaling and the Zipf exponents. In qualitative terms, our results agree with the more holistic idea that urban scaling exponents do not solely emerge from processes occurring within the city boundaries; instead, cities do not represent a closed and non-interacting system and what happens in the entire system (such as flow of people and goods) may affect urban scaling. Similarly to what happens in transportation theory where the product of two cities’ populations is usually assumed to be proportional to the commuting flow between them (such as in gravity models), the population of cities is likely to represent an important factor for the interactions among cities. Under this assumption, the distribution of city population (summarized by the Zipf exponent) may thus represent an indirect proxy for interactions between cities, and the association observed between β and α summarizes how the population distribution affects urban scaling and vice versa. Theoretically, the connection between the exponents β and α would imply that instead of unrelated, Zipf’s law and urban scaling are indeed the two sides of the same coin. However, the non-negligible variability observed in the empirical relationship do not corroborate with such a simple conclusion but suggest that other factors (such as level of socio-economic development and the particular history of an urban system) beyond population distribution may also have a significant effect on the urban scaling of

January 11, 2021 8/17 Empirical data 2.0 A = = 1.06 = 0.26 1.5 = 0.79

= 1 1.0

y = 100% Urban scaling exponent, y = 0% y = 50% y = 75% y = 100% (exponential cutoff) 0.5 0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 Zipf exponent, Zipf exponent, Zipf exponent, Zipf exponent, Zipf exponent, min 13 B C max 8.0 D max 8.5 E 8.0 5.4 7.5 12 7.5 7.0 11 5.2 7.0 6.5 6.5 Pearson'sR = 0.86 Pearson'sR = 0.76 Pearson'sR = 0.85 Pearson'sR = 0.92 10 5.0 6.0 y = 100% y = 100% y = 100% y = 100% Y log Urban GDP, 6.0 (exponential cutoff) = 1.06 = 0.26 5.5 = 0.79 = 0.79 9 logY 0 = 3.53 4.8 loga= 3.18 logb= 1.02 5.5 logb= 1.02 s log Largest city population, s log Largest city population, s log Smallest city population, 6.0 6.5 7.0 7.5 8.0 8.5 6.0 6.5 7.0 7.5 8.0 8.5 6.0 6.5 7.0 7.5 8.0 8.5 6.0 6.5 7.0 7.5 8.0 8.5 Country's urban population, logS Country's urban population, logS Country's urban population, logS Country's urban population, logS Fig 3. Simulating the connection between α and β, and the country scaling relationships. (A) Simulated relationships between α and β under different levels of random variation in the urban scaling law (Eq. 2). Here σy is the percentage of the standard deviation in the empirical values of β. We have a perfect agreement between simulations and the model of Eq. (9) when σy = 0%, and the results become very similar to the empirical data (small blue dots) as the intensity of the random variation increases. (B)-(D) Simulated country scaling laws for σy = 100%. The dashed lines represent the scaling relationships of Eqs. (8), (5), and (6), with parameters estimated from empirical data. We notice that the simulated scaling laws of total urban GDP (Y ) and smallest city population (smin) follow well the empirically adjusted relationships, while the simulated values for the largest city population (smax) underestimate the empirical values. This occurs because large populations are rare and do not get close enough to the imposed maximum smax. (E) Scaling relationship between smax and S when considering that city population values are drawn from a power-law distribution with exponential cutoff. We notice this change makes the simulated results very similar to the empirical ones. The right-most plot in panel A shows a simulated relation between β and α when considering the exponential cutoff.

GDP. Understanding the relative importance of population distribution on urban scaling (and vice versa) for different indicators is an important future contribution that is currently limited by the availability of other world-wide comparable city indicators. Even when considering the unexplained variation in the data, the connection between the two scaling laws uncovered by our work indicates the existence of universal processes governing both laws; however, finding out this commonality for arbitrary urban indicators still represents a challenging task. In the context of urban GDP, our results show that urban systems with small values of the Zipf exponent also tend to present lower increasing returns to scale of GDP (low values of β). An urban system described by a small Zipf exponent has a more balanced population distribution, and consequently, fewer small cities and more large cities when compared with urban systems described by larger Zipf exponents. Thus, countries with proportionally more metropolises tend to have less pronounced increasing returns to scale than those having a small number of large cities. We hypothesize that urban systems with a large number of metropolises may also have a more integrated market whereby these large cities cooperate and develop specialized economic activities. As a result, urban systems with more metropolises would have a smaller degree of agglomeration of economic activities in large cities and so weaker increasing returns to scale for city GDP. In contrast, countries with a small number of large cities have to concentrate almost all complex economic activities in relatively fewer metropolises, which in turn intensifies the increasing returns to scale of urban GDP (high values of β).

January 11, 2021 9/17 It is also interesting to note that high-GDP countries present relatively smaller values of Zipf and urban scaling exponents than mid-GDP and low-GDP countries (Fig. 1A and Fig. 9 in S1 Appendix), that is, developed countries tend to have more metropolises and less pronounced increasing returns to scale of urban GDP. This latter point agrees with the more local observation that “rich cities” of the European Union (West cities) also exhibit smaller scaling exponents for GDP than their “poor” counterparts (East cities) [27]. Large values of β may express an urban system with high economic performance, but because β alone does not define the total urban GDP, large values of β also indicate a significant imbalance between the economic productivity of small and large cities. The association between β and α also suggests that part of this economic inequality may reflect the unbalanced distribution of population. It is also worth mentioning the possibility of the existence of different urban planning regimes [45] that may prevent sharp population agglomerations in developed countries and thus also partially explain the negative association between β and total urban GDP. Our data do not allow a dynamic analysis nor the identification of the causal direction of the association between the exponents β and α. Still, a possible explanation for the observed differences among countries with different levels of development is that the economic policies in less developed countries have focused on large cities, fostering this unbalanced situation and creating cities larger than an hypothetical economically optimal. From a labor-market perspective, these large cities may attract inhabitants from smaller cities, changing both the urban scaling and the Zipf exponents. However, these new inhabitants may mostly find low-paying jobs or even become unemployed, which in turn might partially explain the poorer overall economic performance of less efficient urban systems. The association between both exponents is even more crucial because the world urban population may increase up to 90% by 2100 [46]. This urbanization process is likely to be even more intense in developing countries and has the potential to further undermine their economic performance. Thus, it is not only important to discuss which scaling is desirable but also the population distribution within urban systems.

Materials and Methods Data The data used in this work is a product of two data sets. First, we use the recently released GHSL-FUA or eFUA [38, 39] definition of functional urban areas (FUAs). The eFUAs uses gridded population data from Global Human Settlement Layer (GHSL) [47, 48] and an automated classification approach for producing 9,031 urban boundaries (and population counts) over the entire globe (188 countries) for the year 2015. The eFUAs comprise high-density urban centers and their surrounding commuting zones and aim to capture the functional extent of cities. Second, we use the gridded GDP data provided by Kummu et al. [40]. This data set combines sub-national and national GDP data from different sources with population gridded data (from the GHSL and the HYDE 3.2 [49]) to define three gridded global datasets: Gross Domestic Production per capita (5 arc-min resolution), Human Development Index (5 arc-min resolution), and Gross Domestic Production (30 arc-sec-min resolution). We have used only this latter file representing gridded values of total GDP with a resolution of 30 arc-sec (1 km at equator) for the year 2015 (the same information is also available for the years 1990 and 2000). To define GDP consistently at the grid level across countries, Kummu et al. first calculate the GDP across gridcells within each subnational unit of a given country as the corresponding subnational GDP per capita value multiplied by the gridcell population. Next, the authors sum over these gridcell GDP values and divide by

January 11, 2021 10/17 the sum gridcell populations in a country to define the population-weighted national GDP per capita. They then calculate the ratio between this population-weighted national GDP per capita and the subnational GDP per capita. Finally, they multiply this ratio by the national GDP per capita to obtain the final subnational GDP per capita values (that vary across subnational units in each country), which they then multiply by the population in each gridcell within each subnational unit to obtain the final GDP for every gridcell in a given country. This method ensures that the sum over GDP per capita values at the gridcell level always coincides with officially reported GDP per capita values for each country, and that there is global consistency because the method relies on secondary sources of reported subnational GDP per capita and internationally consistent population grids. The gridded GDP values are reported in 2011 international US dollars using purchasing power parity rates (total GDP-PPP). We overlay the gridded GDP data with eFUA boundary polygons and aggregate the GDP cell values within each polygon for associating a GDP value to each eFUA. We illustrate this procedure in Fig. 1 in S1 Appendix. Next, we group the resulting data by country and select the countries having at least 10 eFUAs. We also removed from our analysis 46 eFUA with null GDP (16 from India; 9 from Ethiopia; 5 from Pakistan; 3 from Sudan; 2 from Niger, Congo, and Chad; and 1 from Argentina, Benin, Egypt, Indonesia, Myanmar, Senegal, and Uganda). These criteria lead us to a data set comprising population s and GDP y of 8,650 functional urban areas from 96 countries.

Estimating the Zipf exponent Zipf’s law (Eq. 1) implies that the complementary cumulative distribution function (CDF) of city population is a power law, F (s) ∼ s−ζ , with ζ ≈ α. We use this connection to estimate the values of α from the data. Specifically, we applied the approach of Clauset-Shalizi-Newman [50] to obtain the exponent α via its maximum Pn likelihood estimate α = n/ ( i=1 ln si/s˜min) , wheres ˜min is the lower bound of the power-law regime, si is the population of the i-th city for a given country such that si ≥ s˜min, and n is the number of city populations in the power-law regime. The value ofs ˜min is also estimated from data by minimizing the Kolmogorov-Smirnov “distance” between the empirical CDF√ of city populations and F (s). The standard error in the Zipf exponent SEα = α/ n can be obtained from the width of the likelihood maximum [50]. We have used the Clauset-Shalizi-Newman method as implemented in the Python module powerlaw [51]. In addition to being a quite popular approach for fitting power-law distributions, Bhattacharya et al. [52] have recently proven that the Clauset-Shalizi-Newman approach yields an unbiased and consistent estimator, that is, as data increases indefinitely the estimated parameters converge (in distribution) to the true values. We show the CDF and Zipf’s law adjusted to each country in S1 Appendix, Sec. 4, where a good agreement is observed in the vast majority of cases. After estimating s˜min, we filter out all cities with population smaller than this threshold in all other analyses, leading us to 4,571 functional urban areas from 96 countries. Thus, the urban scaling laws involve only cities belonging to the power-laws regime, and the countries’ urban GDP (Y ) and countries’ urban population (S) are the aggregated values of urban GDP (yi) and urban population (si) of cities belonging to the power-law regime for each country.

Estimating the urban scaling and the country scaling exponents Urban scaling and country scaling laws are generically represented by a power-law relation between a dependent variable z and an independent variable x

z = g xν , (10)

January 11, 2021 11/17 where g is a prefactor and ν is the power-law exponent. Equation (10) is linearized via a logarithmic transformation log z = log g + ν log x , (11) where log z and log x now represent the dependent and independent variables, log g is the intercept and ν the slope, both being regression coefficients of a corresponding linear model. We have estimated the values of log g and ν by adjusting Eq. (11) to the log-transformed data via robust linear regression with the Huber loss function [53], as implemented in the Python module statsmodels [54]. We further estimate standard errors and confidence intervals for the parameters log g and ν via bootstrapping [55]. We show the urban scaling law adjusted to each country in S1 Appendix, Sec. 4, where a good agreement is observed in the vast majority of cases. The adjusted country scaling laws are shown in Fig. 6 in S1 Appendix. In the case of Figs. 2A-C, we have fixed the power-law exponents (regression coefficients) γ, δ and θ to values obtained from fitting Eq. (9) to the relation β versus α, and only the prefactors (intercepts of the linear model) log Y0, log a, and log b have been considered as free parameters in the regression model.

Fitting our model to the relationship between β and α Our model for the relationship between β and α is completely defined in S1 Appendix, Sec. 2. Depending on whether γ > 1 or γ < 1 and also on whether δ < θ or δ > θ, we have an expression of β as function of α (Eqs. S56-S59, in S1 Appendix) and Eq. (9) is a particular case when γ > 1 and δ < θ. We have adjusted the complete model (that is, without assuming anything about the parameters γ, δ, and θ) to the empirical relation between α and β via the L-BFGS-B nonlinear optimization algorithm [56], as implemented in the Python module lmfit [57] and without any constraint. The standard errors and the confidence intervals on the parameters γ, δ and θ are estimated via bootstrapping [55]. The best fitting parameters (± standard errors) are γ = 1.06 ± 0.01, δ = 0.26 ± 0.05, and θ = 0.79 ± 0.03. This leads to Eq. (9) because the best fitting parameters yield γ > 1 and δ < θ. Figures 10, and 11 in S1 Appendix depict different versions of Fig. 1A where we label all countries and show the standard errors in β and α.

Simulating relations between β and α We simulate the relationship between β and α by generating data at the city level. For a given country population P with Zipf exponent α, we start by generating a list of m city populations S = {s1, . . . , si, . . . , sm} until satisfying the constraint Pm N (0,σγ ) i=1 si ≈ S 10 , where N (0, σγ ) is a Gaussian random variable with zero mean and standard deviation σγ . Each si is drawn from a power-law distribution −(α+1) p(s) ∼ s (compatible with Zipf’s law) within the interval smin < s < smax, where smin and smax are obtained from the country scaling relations (Eqs. 5 and 6) with δ N (0,σδ ) θ N (0,σθ ) multiplicative random variations, that is, smin ∼ S 10 and smin ∼ S 10 , where N (0, σδ) and N (0, σθ) are Gaussian random variables with zero mean and standard deviations σδ and σθ, respectively. We next generate a list of urban indicators β N (0,σy ) Y = {y1, . . . , yi, . . . , ym}, where yi = c si 10 and N (0, σy) is a Gaussian random variable with zero mean and standard deviation σy. In the expression for yi, the value of β is obtained from our model (Eq. 9) while the value of c is chosen to satisfy the Pm condition i=1 yi ≈ Y , where Y is the total urban GDP (S1 Appendix, Sec. 3). The random variation controlled by the parameters σγ , σδ, and σθ mimics the variability observed in the empirical country scaling relationships, and we set their values equal to the standard deviation of the bootstrap estimates of the country scaling exponents (γ, δ,

January 11, 2021 12/17 and θ). On the other hand, the random variation controlled by σy mimics the variability in the urban scaling relationships, and we set its value as a fraction of the standard deviation of the empirical values of β. We have thus applied this procedure by using all empirical values of α and Y to obtain the simulated values of β, Y , smin, smax from the lists S and Y under different values of σy (Fig. 3A). We have also considered a modification of this procedure where the population values were drawn from a power-law distribution with exponential cutoff [50], that is, −(α+1) p(s) ∼ s exp(−s/s0)(s > smin), where s0 is an additional parameter. This modification is necessary for reproducing the empirical behavior of the country scaling between smax and S. Because large city populations are very rare, the simulated values of si obtained from the upper-truncated power-law distribution do not get close enough to the imposed maximum value (smax). This results in the underestimation of smax, as shown in Fig. 3D. After replacing the upper-truncated behavior by the exponential cutoff, we note that the simulated country scaling of smax becomes very similar to the 7 empirical relation (Fig. 3E). In this simulation, we have chosen the value s0 = 3 × 10 to make the simulated values of smax closer to the scaling law adjusted from the empirical data. It is worth mentioning that this change has no effect on the relationship between β and α nor on the other country scaling relations (see the two right-most plots of Fig. 3A). S1 Appendix, Sec. 3 shows more details on how we have implemented this simulation.

Acknowledgments

We thank K. Schmidheiny, J. Suedekum, and S. Thies for useful discussions. This work emerged from ideas discussed at the symposium Cities as Complex Systems (Hanover, July 13th-15th, 2016) which was generously funded by VolkswagenFoundation. HVR acknowledges the support of CNPq (Grants Nos. 440650/2014-3 and 303642/2014-9). DR thanks the Alexander von Humboldt Foundation for financial support under the Feodor Lynen Fellowship. DR is also grateful to the Leibniz Association (project IMPETUS) for financially supporting this research.

Supporting information

S1 Appendix. Supplementary Figures (1-11) and Supplementary Text (Sections 1-4) supporting the results discussed discussed in the main text.

References

1. West GB. Scale: The Universal Laws of Growth, Innovation, Sustainability, and the Pace of Life in Organisms, Cities, Economies, and Companies. New York: Penguin; 2017. 2. Gutenberg B, Richter CF. Frequency of earthquakes in California. Bull Seismol Soc Am. 1944;34(4):185–188. 3. Bouchaud JP. Power laws in economics and finance: some ideas from physics. Quant Finance. 2001;1(1):105–112. doi:10.1080/713665538. 4. Bohorquez JC, Gourley S, Dixon AR, Spagat M, Johnson NF. Common ecology quantifies human insurgency. Nature. 2009;462:911–914. doi:10.1038/nature08631.

January 11, 2021 13/17 5. Ribeiro HV, Costa LS, Alves LGA, Santoro PA, Picoli S, Lenzi EK, et al. Analogies Between the Cracking Noise of Ethanol-Dampened Charcoal and Earthquakes. Phys Rev Lett. 2015;115:025503. doi:10.1103/PhysRevLett.115.025503. 6. Batty M. The New Science of Cities. Cambridge, MA: MIT Press; 2013. 7. Zipf GK. Human Behavior and the Principle of Least Effort: An Introduction to Human Ecology. Manfield Centre, CT: Martino Publishing; 2012.

8. Auerbach F. Das Gesetz der Bev¨olkerungskonzentration. Petermanns Geogr Mitteilungen. 1913;59(74):73–76. 9. Bettencourt LMA, Lobo J, Helbing D, K¨uhnertC, West GB. Growth, innovation, scaling, and the pace of life in cities. Proc Nat Acad Sci USA. 2007;104(17):7301–7306. doi:10.1073/pnas.0610172104.

10. Marshall A. Principles of Economics. London and New York: MacMillan & Co; 1890. 11. Sveikauskas L. The productivity of cities. Q J Econ. 1975;89(3):393–413. doi:10.2307/1885259.

12. Barthelemy M. The structure and dynamics of cities. New York: Cambridge University Press; 2016. 13. Barthelemy M. The statistical physics of cities. Nat Rev Phys. 2019;1(6):406–415. doi:10.1038/s42254-019-0054-2. 14. Lobo J, Alberti M, Allen-Dumas M, Arcaute E, Barthelemy M, Bojorquez Tapia LA, et al. Urban science: Integrated theory from the first cities to sustainable metropolises. Report submitted to the NSF on the Present State and Future of Urban Science. 2020;doi:10.2139/ssrn.3526940. 15. Gabaix X. Zipf’s Law and the Growth of Cities. Am Econ Rev. 1999;89(2).

16. Gabaix X, Ioannides YM. The evolution of city size distributions. In: Handbook of Regional and Urban Economics. vol. 4. Elsevier; 2004. p. 2341–2378. 17. Nitsch V. Zipf zipped. J Urban Econ. 2005;57(1):86–100. doi:10.1016/j.jue.2004.09.002. 18. Soo KT. Zipf’s Law for cities: a cross-country investigation. Reg Sci Urban Econ. 2005;35(3):239–263. doi:10.1016/j.regsciurbeco.2004.04.004. 19. Batty M. Rank clocks. Nature. 2006;444(7119):592–596. doi:10.1038/nature05302. 20. Berry BJL, Okulicz-Kozaryn A. The city size distribution debate: Resolution for US urban regions and megalopolitan areas. Cities. 2012;29(SI1):S17–S23. doi:10.1016/j.cities.2011.11.007.

21. Schmidheiny K, Suedekum J. The pan-European population distribution across consistently defined functional urban areas. Econ Lett. 2015;133:10–13. doi:10.1016/j.econlet.2015.05.013. 22. Li X, Wang X, Zhang J, Wu L. Allometric scaling, size distribution and pattern formation of natural cities. Palgrave Commun. 2015;1(1):15017. doi:10.1057/palcomms.2015.17.

January 11, 2021 14/17 23. Cottineau C. MetaZipf. A dynamic meta-analysis of city size distributions. PLoS One. 2017;12(8):e0183919. doi:10.1371/journal.pone.0183919.

24. Zhang L, Zeng G, Li D, Huang HJ, Stanley HE, Havlin S. Scale-free resilience of real traffic jams. Proc Nat Acad Sci USA. 2019;116(18):8673–8678. doi:10.1073/pnas.1814982116. 25. Mori T, Smith TE, Hsu WT. Common power laws for cities and spatial fractal structures. Proc Nat Acad Sci USA. 2020;117(12):6469–6475. doi:10.1073/pnas.1913014117. 26. Bettencourt LM. The origins of scaling in cities. Science. 2013;340(6139):1438–1441. doi:10.1126/science.1235823. 27. Strano E, Sood V. Rich and poor cities in Europe. An urban scaling approach to mapping the European economic transition. PLoS One. 2016;11(8):e0159465. doi:10.1371/journal.pone.0159465. 28. Gomez-Lievano A, Patterson-Lomba O, Hausmann R. Explaining the prevalence, scaling and variance of urban phenomena. Nat Hum Behav. 2017;1:0012. doi:10.1038/s41562-016-0012. 29. Depersin J, Barthelemy M. From global scaling to the dynamics of individual cities. Proc Nat Acad Sci USA. 2018;115(10):2317–2322. doi:10.1073/pnas.1718690115. 30. Ribeiro FL, Meirelles J, Ferreira FF, Neto CR. A model of urban scaling laws based on distance dependent interactions. Royal Soc Open Sci. 2017;4(3):160926. doi:10.1098/rsos.160926.

31. Keuschnigg M. Scaling trajectories of cities. Proc Nat Acad Sci USA. 2019;116(28):13759–13761. doi:10.1073/pnas.1906258116. 32. Pumain D, Rozenblat C. Two metropolisation gradients in the European system of cities revealed by scaling laws. Environ Plan B Urban Anal City Sci. 2019;46(9):1645–1662. doi:10.1177/2399808318785633.

33. Balland PA, Jara-Figueroa C, Petralia SG, Steijn MP, Rigby DL, Hidalgo CA. Complex economic activities concentrate in large cities. Nat Hum Behav. 2020;4:248–254. doi:10.1038/s41562-019-0803-3. 34. Bettencourt LMA, Yang VC, Lobo J, Kempes CP, Rybski D, Hamilton MJ. The interpretation of urban scaling analysis in time. J R Soc Interface. 2020;17(163):20190846. doi:10.1098/rsif.2019.0846. 35. Gomez-Lievano A, Youn H, Bettencourt LMA. The statistics of urban scaling and their connection to Zipf’s law. PLoS One. 2012;7(7):e40393. doi:10.1371/journal.pone.0040393.

36. Alves LGA, Ribeiro HV, Lenzi EK, Mendes RS. Empirical analysis on the connection between power-law distributions and allometries for urban indicators. Physica A. 2014;409:175–182. doi:10.1016/j.physa.2014.04.046. 37. Arcaute E, Hatna E, Ferguson P, Youn H, Johansson A, Batty M. Constructing cities, deconstructing scaling laws. J R Soc Interface. 2014;12(102):20140745. doi:10.1098/rsif.2014.0745.

January 11, 2021 15/17 38. Schiavina M, Moreno-Monroy A, Maffenini L, Veneri P. GHSL-OECD Functional Urban Areas. JRC Technical Report. 2019;.

39. Moreno-Monroy AI, Schiavina M, Veneri P. Metropolitan areas in the world. Delineation and population trends. J Urban Econ. 2020; p. 103242. doi:10.1016/j.jue.2020.103242. 40. Kummu M, Taka M, Guillaume JH. Gridded global datasets for gross domestic product and Human Development Index over 1990–2015. Sci Data. 2018;5:180004.

41. Bettencourt LM, Lobo J, Strumsky D, West GB. Urban scaling and its deviations: Revealing the structure of wealth, innovation and crime across cities. PLoS One. 2010;5(11):e13541. doi:10.1371/journal.pone.0013541. 42. Pumain D, Moriconi-Ebrard F. City size distributions and metropolisation. Geojournal. 1997;43(4):307–314. doi:10.1023/A:1006833025206.

43. Zhang J, Yu T. Allometric scaling of countries. Physica A. 2010;389(21):4887–4896. doi:10.1016/j.physa.2010.06.059. 44. Rybski D, Reusser DE, Winz AL, Fichtner C, Sterzel T, Kropp JP. Cities as nuclei of sustainability? Environ Plan B. 2017;44(3):425–440. doi:10.1177/0265813516638340.

45. Hsieh CT, Moretti E. Housing constraints and spatial misallocation. American Economic Journal: Macroeconomics. 2019;11(2):1–39. doi:https://doi.org/10.1257/mac.20170388. 46. Jiang L, O’Neill BC. Global urbanization projections for the Shared Socioeconomic Pathways. Global Environ Chang. 2017;42:193–199. doi:10.1016/j.gloenvcha.2015.03.008. 47. Freire S, Pesaresi M. GHS population grid, derived from GPW4, multitemporal (1975, 1990, 2000, 2015). European Commission, Joint Research Centre (JRC). 2015;.

48. Pesaresi M, Ehrlich D, Florczyk AJ, Freire S, Julea A, Kemper T, et al. The global human settlement layer from landsat imagery. In: 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS); 2016. p. 7276–7279. 49. Klein Goldewijk K, Beusen A, Janssen P. Long-term dynamic modeling of global population and built-up area in a spatially explicit way: HYDE 3.1. The Holocene. 2010;20(4):565–573. doi:10.1177/0959683609356587. 50. Clauset A, Shalizi CR, Newman MEJ. Power-Law Distributions in Empirical Data. SIAM Rev. 2009;51(4):661–703. doi:10.1137/070710111. 51. Alstott J, Bullmore DP. powerlaw: a Python package for analysis of heavy-tailed distributions. PLoS One. 2014;9(1). doi:10.1371/journal.pone.0085777.

52. Bhattacharya A, Chen B, van der Hofstad R, Zwart B. Consistency of the PLFit estimator for power-law data. arXiv:200206870. 2020;. 53. Huber PJ, Ronchetti EM. Robust Statistics. 2nd ed. New Jersey: Wiley & Sons; 2009.

54. Seabold S, Perktold J. statsmodels: Econometric and statistical modeling with Python. In: 9th Python in Science Conference; 2010.

January 11, 2021 16/17 55. Efron B, Tibshirani RJ. An Introduction to the Bootstrap. Abingdon: Chapman & Hall/CRC; 1994.

56. Byrd RH, Lu P, Nocedal J, Zhu C. A limited memory algorithm for bound constrained optimization. SIAM J Sci Comput. 1995;16(5):1190–1208. doi:10.1137/0916069. 57. Newville M, Stensitzki T, Allen DB, Ingargiola A. LMFIT: Non-Linear Least-Square Minimization and Curve-Fitting for Python. Zenodo. 2014;doi:10.5281/zenodo.11813.

January 11, 2021 17/17 S1 Appendix for

Association between population distribution and urban GDP scaling

H. V. Ribeiro , M. Oehlers, A. I. Moreno-Monroy, J. P. Kropp, D. Rybski Corresponding∗ authors. E-mail: hvr@dﬁ.uem.br or [email protected] ∗ ∗

1 Supplementary Figures

Gridded Log(GDP) 0.81 - 3.39 3.39 - 3.83 3.83 - 4.26 4.26 - 5.30 5.30 - 7.25

FIG. 1. Illustration of the procedure used for obtaining population s and GDP y of cities over the entire globe. The world map shows the boundaries of GHSL-OECD Func- tional Urban Areas (GHS-FUA or eFUA). This definition of urban areas has been recently released by the European Commission’s science and knowledge service, with the aim of providing an internationally-comparable and harmonized definition of urban areas. The eFUA definition does not rely on local administrative boundary units and commuting flows data; instead, it uses gridded data from the Global Human Settlement Layer (GHSL) and an automated classification approach for producing urban boundaries [JRC Tech. Rep. (2019) and J. Urban. Econ., 103242 (2020)]. The shape file is provided for the year of 2015 together with population values for each eFUA. We combine eFUA boundary data with gridded GDP data provided by Kummu et al. [Sci. Data 5, 180004 (2018)]. In more detail, we use the total GDP-PPP (purchasing power parity, in constant 2011 international USD) data with 30 arc-sec resolution ( 1 km at the equator, file name: GDP PPP 30arcsec v2.nc) for the year 2015 to overlay with eFUA boundaries. Next, we aggregate ≈ all grid cell values within each eFUA boundary to associate a GDP value to each urban area of the world. We illustrate this idea by highlighting a small area in the southern part of Brazil, where we also plot the GDP grid cells.

2 A B 25 25

20 20

15 15

10 10 Frequency Frequency

5 5

0 0 0.5 1.0 1.5 2.0 2.5 0.8 1.0 1.2 1.4 1.6 Zipf exponent, Urban scaling exponent,

FIG. 2. Probability distributions of the Zipf exponent α and the urban scaling exponent β. (A) The bars represent a histogram for the values of α and the continuous line is a kernel density estimation (not normalized) for the same data. The vertical dashed line indicates the value α 1. The average value and the standard deviation of α are 1.24 and 0.38, respectively. (B) The bars = represent a normalized histogram for the values of β and the continuous line is a kernel density estimation (not normalized) for the same data. The vertical dashed line indicates the value β 1. The average value and the standard deviation of β are 1.16 and 0.17, respectively. =

3 r

A B , 0.4 0.4 n o

i , t n a l o 0.2 0.2 i e t r a r l o e c r

r 0.0

k 0.0 o n c

a r n

n 0.2

0.2 s a r a m e r a P 0.4

e 0.4 p S

0.8 1.0 1.2 1.4 0.8 1.0 1.2 1.4

Center of Zipf exponent window, 0 Center of Zipf exponent window, 0

FIG. 3. Correlation between the urban scaling exponent β and the Zipf exponent α over sliding windows of values of α. (A) Spearman rank correlation (ρr) of the relationship between β and α estimated within a sliding window centered in α . (B) Pearson correlation (ρ) of the relationship between β and α estimated within a sliding window′ centered in α . For both correlation measures, the size of sliding windows is ∆α 0.4 and the step is 0.02. In addition,′ the red star markers indicate the correlation values that statistically′ significant at the 95% confidence = level (p-values smaller than 0.05), while the gray circles represent non-significant correlation values (p-values larger than 0.05). Similar results are obtained for ∆α between 0.15 and 0.55. In both panels, the continuous black lines around zero represent average values′ of the correlation measures estimated after shuffling the values of β (1000 realizations), and the gray shaded region stands for the 95% bootstrap confidence intervals. We notice that the correlation measures are larger and statistically significant for values α around 1.1 and that the significant values cannot be explained ′ by the distribution of β values. This range containing the largest and significant values of ρ and ρr mostly correspond with mid-range of α values in Eq. (9) [Eq. S56 in Supplementary Text], where our model predicts a linear association between β and α (see main text for details).

4 1333 more likelly 15 6 > C I B

C d 103 more likelly I n a

B 10 l l e 6 e d d >

o o C I d m m Strong or very strong C A C n I evidence for our model I

a A B

- -

C t t

I s s n n

A 5 o o c c IC IC A B = = IC IC A B 0

FIG. 4. Comparison between the model of Eq. (9) [Eq. S56 in Supplementary Text] and a constant (intercept-only) model (β const). The left bar shows the diﬀerence between the

AIC (Akaike information criterion) values (∆AIC) estimated for an intercept-only model (AICconst) = and for the model of Eq. (9) [Eq. S56 in Supplementary Text] (AICmodel). The model with the lowest AIC value gives a better description of the data. Thus, ∆AIC 0 indicates that the model of Eq. (9) [Eq. S56 in Supplementary Text] is a better description for data. As a rule of thumb, it is > common to consider 6 ∆AIC 10 as a strong evidence and ∆AIC 10 as a very strong evidence (our case, since ∆AIC 14.4) in favor of the model of Eq. (9) [Eq. S56 in Supplementary Text]. < < ≥ In addition, because the AIC value is proportional to the negative of log-maximum-likelihood of ≈ the model, the quantity e∆AIC approximates how more likely the model of Eq. (9) [Eq. S56 in Supplementary Text] is to describe the data when compared with an intercept-only model. Thus, we notice that Eq. (S56) is 1333 times more likely to describe the relationship between β and α than the intercept-only model. The right bar shows the diﬀerence between the BIC (Bayesian ≈ information criterion) values (∆BIC) estimated for an intercept-only model (BICconst) and for the model of Eq. (9) [Eq. S56 in Supplementary Text] (BICmodel). The interpretation of BIC and AIC values are analogous. The results indicate a strong evidence (∆BIC 9.3) in favor of Eq. (9) [Eq. S56 in Supplementary Text] and our model is 103 times more likely to describe the ≈ empirical data than the intercept-only model. It is worth remembering that both AIC and BIC ≈ values account for the fact that Eq. (9) [Eq. S56 in Supplementary Text] has more free parameters than the intercept-only model (3 versus 1 parameter); moreover, BIC values penalize the number of parameters more heavily.

5 p-value = 0.001 * 1.22 p-value = 0.003 * 1.22 r 0.6 A 0.5 B 0.5 p-value = 0.009 0.4 p-value = 0.032 0.4 0.3 p-value = 0.125 0.3 p-value = 0.078 0.2 0.2 0.1 0.1 Pearson correlation,

Spearman rank correlation, 0.0 0.0

* ) * )

* ) * ) Linear region Linear region 1) ( 1) ( (1 < < ( (1 < < ( Low plateau region High plateau region Low plateau region High plateau region

FIG. 5. Correlations in each range of values of α defined in Eq. (9) [Eq. S56 in Supple- mentary Text]. (A) Spearman rank correlation ρr between β and α in the low plateau (α 1), linear regime (1 α α ), and high plateau (α α ) regions of Fig. 1 in the main text. Panel (B) ≤ shows the same for the Pearson∗ correlation coefficient∗ ρ. In both panels, we indicate above the bars < < ≥ the p-values testing the significance of the correlation coefficient (permutation test). We notice that

ρr and ρ are significantly higher in the linear region. We further observe that the correlations are γ 1 statistically significant at 99% confidence level only in the linear region. Here α 1 δ 1.22 [Eq. (9) of the main text or Eq. S56 in Supplementary Text, and Fig. 15]. ∗ − = + ≈

6 n x i

A B a C m m = 1.31±0.11 s 6.00 = 0.20±0.05 = 0.73±0.05

s 8.0

log Y0 = 1.73±0.79 g log a = 3.65±0.36 log b = 1.47±0.37 g

13 o l o Y

, 5.75

, 7.5 n g n o o i o l t i 12 t , a 5.50 l a P l

u 7.0 u D p p o G o

p 5.25

p n

y 6.5 a t y i t b i c r

t 5.00 10 U t s s

e 6.0 l e l g a

4.75 r 9 a m L S 6.5 7.0 7.5 8.0 8.5 6.5 7.0 7.5 8.0 8.5 6.5 7.0 7.5 8.0 8.5 Country population, log S Country population, log S Country population, log S

FIG. 6. Fitting country scaling relationships directly to data. (A) Scaling relation between total urban GDP of countries (Y ) and country’s urban population values (S). (B) Scaling law between the smallest city population of a country (smin) and country’s urban population (S). (C) Scaling law between the largest city population of a country (smax) and country’s urban population (S). In the three panels markers represent the values for each country and the dashed lines are the country scaling relationships of Eqs. (8), (5), and (6) [Eqs. (S6), (S7), and (S8) in Supplementary

Text], where the exponents γ, δ, and θ (and the prefactors log Y0, log a, and log b) are obtained via robust linear regression of the log-transformed data. The best fitting parameters are shown in the panels ( standard errors) and gray shaded regions stand for the 95% confidence band of the models. As discussed in the main text, the exponents estimated directly from the scaling laws ± are very similar to the values obtained by fitting Eq. (9) [Eq. S56 in Supplementary Text] to the relationship between β and α (see Fig. 2D of main text for a comparison).

7 A B C l l

e e 15

d 15 d 1.0 o o m m C C I I

0.8 A B 10 - -

10 t t

s AIC > 6 s BIC > 6 0.6 n n Strong or very strong

o Strong or very strong o

c evidence for our model c evidence for our model C C 0.4 I 5 I 5 A B

= =

0.2 C C I I

A 0 B 0 Country scaling exponents 0 10 20 30 40 50 0 10 20 30 40 50 0 10 20 30 40 50 GDP percentile threshold, %GDP * GDP percentile threshold, %GDP * GDP percentile threshold, %GDP *

FIG. 7. Robustness of the model of Eq. (9) [Eq. S56 in Supplementary Text] against different GDP threshold. (A) The colored markers show the values of γ, δ, and θ estimated by fitting (via the L-BFGS algorithm) Eq. (9) [Eq. S56 in Supplementary Text] to the relationship between β and α after selecting countries with total urban GDP higher than the percentile GDP threshold (%GDP ). The shaded regions are 95% bootstratp confidence intervals and the dashed lines represent the∗ values when considering all data. We notice the exponent estimates are quite stable. (B) Difference between the AIC (Akaike information criterion) values (∆AIC) estimated for an intercept-only model (AICconst) and for the model of Eq. (9) [Eq. S56 in Supplementary

Text] (AICmodel) as a function of the percentile GDP threshold (%GDP ). (C) Diﬀerence between the BIC (Bayesian information criterion) values (∆BIC) estimated for∗ an intercept-only model

(BICconst) and for the model of Eq. (9) [Eq. S56 in Supplementary Text] (BICmodel) as a function of the percentile GDP threshold (%GDP ). Results show strong evidence (∆AIC 6 and ∆BIC 6) for the model of Eq. (9) [Eq. S56 in Supplementary∗ Text] at practically all threshold values. > >

8 2.5

2.0

1.5

1.0

0.5 Empirical

Probability distribution Simulations with x=100%

Simulations with x=100% (exponential cutoff) 0.0 0.4 0.2 0.0 0.2 0.4 Fluctuations around the model

FIG. 8. Probability distributions of the random variations around the model of Eq. (9) [Eq. S56 in Supplementary Text]. The diﬀerent curves show kernel density estimations of the random variations around the model of Eq. (9) when considering the empirical data and the simulations obtained for σy 100%. We observe that the distributions obtained from the simulations are very similar to the one emerging from the empirical data, especially when including = an exponential cutoﬀ in the power-law distribution associated with the simulated city sizes.

9 A B slope = -0.11±0.04 slope = -0.07±0.02

, 2.5 Pearson = -0.32 (p-value = 0.002) t 1.6 Pearson = -0.49 (p-value = 0.000) n e n

, o t p n 2.0 1.4 x e e n

o g p n i

x 1.5 l 1.2 e a

f c p s i

Z 1.0 n 1.0 a b r

0.5 U 0.8 9 10 11 12 13 9 10 11 12 13 Urban GDP, logY Urban GDP, logY

FIG. 9. Association between the exponents α and β and total urban GDP Y . (A) Scatter plot of Zipf exponent α versus the total urban GDP Y . (B) Scatter plot of urban scaling exponent β versus the total urban GDP Y . The dashed line in each panel represents a linear regression with slope indicated within the plots (where we also shoe the Pearson correlation coeﬃcient and the corresponding p-values). We observe that both exponents also signiﬁcantly associated with urban GDP Y , such that an increase in Y tends to be followed by a decrease in α and β.

10 1.4

KEN

THA 1.3 BENPER SRB CZE

YEM CODNGA VEN NER CHL SDN BGD HTI LKA HND 1.2 AGO TUR TZA OMN COL TGO SEN UZB GHACIVINDMEX ECU CHE MAR GIN PHL VNM ESP SWE FRA BEL MLI UGA 1.1 GRC IRN MMR BOL NLD LBYCMR IDN SYR JPN POL ZWE AZE PAKDEUDOMCHN 2.0 = SAU ZAFUSA CAN ARG TUN ITA TJK 1.0 KOR BRABLR GBR

SSD SOM BDI MOZ ETH 0.8 1.0 1.2 1.4 1.5 AFG GTM UKRROU TCD KEN NIC PER THAHUN SRB BEN CZEPRK RUSBFA YEM ZMB NER CODCHLNGA VENSDN MYS LKAAGO HND BGD HTI COL TURTZA OMNTGO NPL SENUZB GHACIVINDMEX ECU PHLMAR VNMCHEGINESP BEL GRC MLIIRN SWE UGAFRAMMR BOL NLD EGY IDNPOLSYR CMRLBYZWEJPN AZE = 1 SAUZAFUSA PAKDEUDOMCHN DZA IRQ TWN AUSKORCAN BRAARGBLRTUNITA GBR TJK CUB 1.0 KAZ

TKM = 1.06±0.01 = 0.26±0.05

Urban scaling exponent, = 0.79±0.03 0.5 0.5 1.0 1.5 2.0 2.5 3.0 Zipf exponent,

FIG. 10. Association between the urban scaling exponent β and the Zipf exponent α. The markers indicate the values of β versus α for each country in our data set. The three marker colors group countries according to the tercile values of the total urban GDP distribution (for instance, High-GDP countries have highest 33% GDP values, as the main text). The horizontal dashed line shows the β 1, while the inclined dashed line represents the β α relationship. The ≈ continuous line shows the model of Eq. (9) [Eq. S56 in Supplementary Text] adjusted to the data = = and the gray shaded region stands for the 95% bootstrap confidence band. The colored shaded regions represent the different intervals of α defined in Eq. (9) [Eq. S56 in Supplementary Text]. All countries are labeled according to their ISO codes (the three letters over each marker). The blue frame indicates a region that has been magnified on the right.

11 = 1.06±0.01 = 0.26±0.05 2.0 = 0.79±0.03

SSDSSD

BDI BDIMOZ SOMSOM MOZ

ETHETH

1.5 AFGAFG

ROU

GTM ROUUKR GTM UKR TCDTCD

NIC KENKEN NIC

HUN THA HUN PER THA PERBEN PRK SRB BEN RUS BFA SRB CZECZE PRK RUSBFA ZMBZMB YEM

YEM COD VEN NGA NER CODNGACHL VEN MYS NER CHL SDN BGD MYS SDNHTI LKA HND BGD HTI LKA AGO HND TUR AGO TURTZATZA OMNOMN TGO COL NPL COL TGO NPL UZB SEN GHA IND MEX ECU SENUZB GHAINDMEXCIV ECU CHE CIV MAR GIN PHL CHE MAR VNM PHL GINESP VNM ESPSWE FRA MLI UGA BEL SWE FRA IRN GRC UGA BEL GRC MLIIRN MMR BOLBOL MMR NLD

LBY CMR IDN NLD SYR CMRLBY JPN AZE EGY IDN POL SYR PAK ZWE CHN JPN DEU AZE EGY POL PAKCHNZWEDOM DZA IRQ USA DEU ZAF DZA IRQ SAU USA DOM CAN SAUZAFARG TUN TJK ITA TWN AUS BRA BLR KOR ARG CAN TUN GBR TJK 1.0 TWN AUS KOR BRABLR ITA GBR CUBCUB KAZKAZ

TKMTKM Urban scaling exponent,

0.5

0.5 1.0 1.5 2.0 2.5 3.0 Zipf exponent,

FIG. 11. Association between the urban scaling exponent β and the Zipf exponent α. The same as Fig. (10) but with error bars indicating standard errors. The three colors group countries according to the tercile values of the total urban GDP distribution (for instance, High- GDP countries have highest 33% GDP values). The horizontal dashed line shows the β 1, while the inclined dashed line represents the β α relationship. The continuous line shows the model of ≈ = Eq. (9) [Eq. S56 in Supplementary Text] adjusted to data and the gray shaded region stands for = the 95% bootstrap confidence band. The colored shaded regions represent the different intervals of α defined in Eq. (9) [Eq. S56 in Supplementary Text]. All countries are labeled according to their ISO codes (the three letters over each marker).

12 Supplementary Text

1. CONNECTION BETWEEN URBAN SCALING AND ZIPF EXPONENTS

Let us start by considering an urban system (a country) composed of m cities with populations s1, s2, . . . , sm . Suppose now that we sort these population values s in descending order smax, . . . , s, . . . , smin (where smax and { } smin are respectively the largest and smallest city population in the system) { } and associate a rank variable r deﬁned by 1, 2,...,r,...,m to each city population. Zipf’s law for cities states that rank r and city population s are { } related via r s α , (S1) − where α 0 is a power-law exponent∼ or the Zipf exponent. Zipf’s law implies that the complementary cumulative distribution of city population F s can > be written as F s s ζ, with ζ α, and, consequently, the probability ( ) density distribution P S− of city population is ( ) ∼ ≈ α 1 P (s ) s for smin s smax . (S2) −( + ) From Eq. (S2) we can( ) also∼ write the frequency< < function or the absolute number of cities with population s as

α 1 p s ks for smin s smax , (S3) −( + ) where k is a normalization( ) = constant determined< by<

smax smax S sp s k s αds (S4) s s smin min − and S is the total population= =Q of the( ) urban≈ S system. On the other hand, the urban scaling hypothesis establishes a power-law relationship between an urban indicator y (e.g. the urban GDP) and the population s of a city, that is,

y csβ , (S5) where c is a positive constant, and β= 0 is the urban scaling exponent. Our work has shown that the values of β and α are related to each other, and > 13 we propose that the relationship between these two exponents emerge from country scaling relationships. To derive the relationship between β and α, we consider S as the total population of an urban system and Y as the total quantity of the city indicator y for the whole urban system (e.g. total urban GDP). It is worth noticing that the values of S and Y represent the empirical values of total population and total GDP of all cities of a country within the power-law regime (that is, all cities with population larger than the lower cutoff obtained after fitting the Zipf’s law to the data, see Material and Methods in the main text for details). In addition to the urban scaling that considers relations between population and indicators of cities within one country, we further assume three different country scaling relations between population and indicators across countries. The first one is a generalization of the urban scaling hypothesis at the country level

γ Y Y0S , (S6)

= where Y0 is a constant, and γ is the country scaling exponent associated with the indicator Y . The other two relationships relate the population of the smallest (smin) and largest (smax) cities in an urban system to the total population S, that is,

δ smin aS (S7) s bSθ , (S8) max = = where a and b are positive constants, and δ 0 and θ 0 are two additional country scaling exponents. We have empirically veriﬁed these three country > > scaling relationships in our data (Fig.6 and Fig. 2 of the main text); however and as we shall prove, the relationship of Eq. (S6) can also be derived from Zipf’s law (Eq. S3), urban scaling (Eq. S5) and the two country scalings for smin and smax (Eqs. S7 and S8). Qualitatively speaking, we consider that the existence of country scaling relationships constrains the population and urban indicator values in a country which in turn ties the values of β and α. To ﬁnd this connection, we notice that the total indicator Y can be ex-

14 pressed by using Zipf’s law and urban scaling, that is,

smax smax smax Y yp s csβp s kc sβ α 1ds . (S9) s s smin s smin min − − By determining= the=Q constant( ) = k=Qfrom Eq.( ) ≈ (S4)S and solving the integral in Eq. (S9), we should obtain the country scaling of Eq. (S6) in terms of the exponents α and β. Thus, if we assume γ, δ, and θ to be ﬁxed, the changes in α and β must be constrained in order to yield the particular value of γ imposed by country scaling in Eq. (S6).

Case A: α β 1

To start our≠ calculations,≠ we first consider the general case where α β 1. Under this assumption, Eq. (S4) yields α α ≠ ≠ α 1 Ssmaxsmin k α α , (S10) smaxsmin smaxs ( − ) min and then by replacing smin with= Eq. (S7) and smax with Eq. (S8), we find − α 1 ab αSS δ θ α k . (S11) abαSδ θα aαbS(θ+δα) ( − )( ) By solving the integral in Eq.= (S9),+ we find + kc − Y sβ α sβ α , (S12) β α max min − − and after replacing k by Eq.= (S11), smin by− Eq. (S7), and smax by Eq. (S8), − we have c α 1 S aβbαSδβ θα aαbβSθβ δα Y α β abαSδ+θα aαbSθ δα+ ( − ) − = c α 1 S aβbαSλ1 + aαbβSλ2 + − − (S13) α β abαSλ3 aαbSλ4 ( − ) − = c α 1 S Λ1 Λ2 − −, α β Λ3 Λ4 ( − ) − where = − − λ1 δβ θα

λ2 θβ δα = + (S14) λ δ θα 3 = + λ θ δα . 4 = + =15+ and

β α λ1 Λ1 a b S

α β λ2 Λ2 a b S = (S15) Λ abαSλ3 3 = Λ aαbSλ4 . 4 = = Because S represents the total population of an urban system or a country, we take the limit of large S in Eq. (S13). To do so, we need to know the conditions under which Λ1 dominates over Λ2 in the numerator of Eq. (S13), that is, when

λ1 λ2 δβ θα θβ δα (S16) > δ θ β α 0 , + > + ( − )( − ) > and also the conditions under which Λ3 dominates over Λ4 in the denominator of Eq. (S13), that is, when

λ3 λ4 δ θα θ δα (S17) > δ θ 1 α 0 . + > + ( − )( − ) > By investigating the inequality in Eq. (S16), we end up with the decision tree of Figure 12. Similarly, we ﬁnd the decision tree of Figure 13 for the inequality in Eq. (S17). Having the results of Figures 12 and 13, we can ﬁnd the behavior of Y for large S in Eq. (S13) for each of the following conditions.

16 α > β yes no

δ > θ δ < θ

yes yes no no

False and True and False and True and Λ2 dominates Λ1 dominates Λ2 dominates Λ1 dominates

FIG. 12. Decision tree associated with the inequality λ1 λ2 (Eq. S16) that deﬁnes whether Λ1 dominates over Λ2 (or vice-versa) in Eq. (S13). > α > 1 yes no

δ > θ δ < θ

yes yes no no

False and True and False and True and Λ4 dominates Λ3 dominates Λ4 dominates Λ3 dominates

FIG. 13. Decision tree associated with the inequality λ3 λ4 (Eq. S17) that deﬁnes whether Λ3 dominates over Λ4 (or vice-versa) in Eq. (S13). >

Case A.1: α β, α 1, and δ θ

Under these> conditions,> we ﬁnd< that Λ1 and Λ3 dominate in Eq. (S13), yielding

c α 1 aβ 1 Y S1 λ1 λ3 α β − + − (S18) c(α − 1)aβ 1 ≈ S1 δ βδ . α − β − ( − ) − + ≈ By comparing the result of Eq. (S18−) with the country scaling in Eq. (S6), we ﬁnd γ 1 β 1 for α β & α 1 & δ θ . (S19) δ − It is worth noticing= that+ the previous{ conditions> > depend< on} β, but by com-

17 bining α β & α 1 and eliminating β, we can rewrite Eq. (S19) as γ 1 γ 1 > β 1 > for α max 1, 1 & δ θ . (S20) δ δ − − Equation (S20=) thus+ predicts that β> is a constant, + that is, not< dependent on γ 1 α when α max 1, 1 δ & δ θ . − Case A.2:{ >α β,α +1, and δ θ< }

For this case,> we ﬁnd< that Λ1 we have< < } γ θ 1 θ γ 1 β 1 α for γ 1 & 1 α 1 & δ θ . (S23) δ δ θ + − − Equation= (S23+) thus− predicts a linear < decreasing+ correspondence< < < between β γ 1 and α when γ 1 & 1 θ α 1 & δ θ . − Case A.3: α{ <β, α 1,+ and<δ <θ < }

Under these< assumptions,> Λ2 < } that after eliminating β from the condition leads to γ δ 1 δ γ 1 β 1 α for γ 1 & 1 α 1 & δ θ . (S26) θ θ δ + − − Similarly= to Eq.+ (S23− ), the result of> Eq. (S26< )< predicts+ a linear< increasing γ 1 correspondence between β and α when γ 1 & 1 α 1 δ & δ θ . It is worth noticing that the conditions underlying Eqs. (S23−) and (S26) { > < < + < } are mutually exclusive, that is, only one of the linear correspondences exists depending on whether γ 1 (super-linear country scaling, Eq. S26) or γ 1 (sub-linear country scaling, Eq. S23). > < Case A.4: α β, α 1, and δ θ

For this case,< we ﬁnd< that Λ2 (Eq. S20). < 19 Case A.5: α β, α 1, and δ θ

Under these> conditions,> Λ2 and> Λ4 dominate the behavior of Eq. (S13) for large S, yielding c α 1 bβ 1 Y S1 λ2 λ4 α β − + − (S30) c(α − 1)bβ 1 ≈ S1 θ βθ . α − β − ( − ) − + By comparing the previous result≈ with Eq. (S6), we ﬁnd − γ 1 β 1 for α β & α 1 & δ θ , (S31) θ − and after eliminating= +β from the condition,{ > we have> > } γ 1 γ 1 β 1 for α max 1, 1 & δ θ . (S32) θ θ − − = + > + > γ 1 Equation (S32) thus predicts that β is a constant when α max 1, 1 θ & δ θ . − { > + > Case} A.6: α β, α 1, and δ θ

For this case,> Λ2 and< Λ3 dominate> the behavior of Eq. (S13) for large S and produce c 1 α aα 1bβ α Y S1 λ2 λ3 α β− − + − (S33) c(1 − α)aα 1bβ α ≈ S1 δ βθ α δ θ . α − β− − ( − ) − + + ( − ) By comparing the previous≈ result with Eq. (S6), we ﬁnd − γ δ 1 δ β 1 α for α β & α 1 & δ θ , (S34) θ θ + − and after= eliminating+ β from− the condition,{ > we have< > } γ δ 1 δ γ 1 β 1 α for γ 1 & 1 α 1 & δ θ . (S35) θ θ δ + − − Equation= (S35+) predicts − a linear decreasing < correspondence+ < < between> β and α γ 1 when γ 1 & 1 δ α 1 & δ θ . − { < + < < > 20} Case A.7: α β, α 1, and δ θ

For this case,< Λ1 and> Λ4 dominate> the behavior of Eq. (S13), yielding c α 1 aβ αbα 1 Y S1 λ1 λ4 β α− − + − (S36) c(α − 1)aβ αbα 1 ≈ S1 θ βδ α θ δ . β − α− − ( − ) − + + ( − ) By comparing the previous≈ result with Eq. (S6), we ﬁnd − γ θ 1 θ β 1 α for α β & α 1 & δ θ , (S37) δ δ + − and after= eliminating+ β from− the condition,{ < we have> > } γ θ 1 θ γ 1 β 1 α for γ 1 & 1 α 1 & δ θ . (S38) δ δ θ + − − Equation= (S38+) predicts − a linear increasing > correspondence< < + between> β and α γ 1 when γ 1 & 1 α 1 θ & δ θ . Similarly to cases A.2 and A.3, the conditions underlying Eqs.− (S35) and (S38) are mutually exclusive, that { > < < + > } is, only one of the linear correspondences exists depending on whether γ 1 (super-linear country scaling, Eq. S38) or γ 1 (sub-linear country scaling, > Eq. S35). < Case A.8: α β, α 1, and δ θ

Under these< conditions,< Λ1 and> Λ3 dominate the behavior of Eq. (S13), yielding c 1 α aβ 1 Y S1 λ1 λ3 β α − + − (S39) c(1 − α)aβ 1 ≈ S1 δ βδ . β − α − ( − ) − + By comparing the previous result≈ with Eq. (S6), we ﬁnd − γ 1 β 1 for α β & α 1 & δ θ , (S40) δ − and after eliminating= +β from the condition,{ < we have< > } γ 1 γ 1 β 1 for α min 1, 1 & δ θ . (S41) δ δ − − = + <21 + > γ 1 Equation (S41) predicts that β is a constant when α min 1, 1 δ & δ θ . − { < + > Similarly to cases A.1-A.4, the combined behavior of Eqs. (S32), (S35), } (S38) and (S41) produces a functional dependence of β on α composed by three continuous line segments: an initial horizontal plateau (Eq. S41) followed by an increasing (Eq. S38, when γ 1) or decreasing (Eq. S35, when γ 1) linear function followed by another horizontal plateau (Eq. S32). > Case< B: α 1 and β 1

We now consider= the≠ particular case α 1 and β 1 in order to ensure the continuous behavior of β as function of α at α 1. Under these assumptions, = ≠ the normalization constant in Eq. (S4) is = S k , (S42) ln smax smin = and by replacing smin and smax with Eqs. (S7) and (S8), we ﬁnd S k . (S43) b θ δ ln aS = − The solution of the integral in Eq. (S9 ) for α 1 is kc Y sβ 1 sβ=1 , (S44) β 1 max min − − = − and after replacing k by Eq. (S43−), smin by Eq. (S7), and smax by Eq. (S8), we ﬁnd cS bβ 1S θ βθ aβ 1S δ βδ Y β 1 − −ln+ b Sθ δ− − + ⎛ a− ⎞ (S45) = cS Ω Ω − − ⎝ 1 2 , ⎠ β 1 ln b Sθ δ ⎛ a− ⎞ = − where − ⎝ ⎠ β 1 θ βθ Ω1 b S . (S46) Ω aβ− 1S− δ+βδ 2 = − − + = 22 It is worth noticing that Eq. (S45) does not produce a “pure” power-law behavior for large S due to the logarithmic function in its denominator. How- ever, the logarithmic function changes much slower than the power-law functions in the numerator, allowing us to approximate the behavior of Eq. (S45) by a power-law function for large S.

β > 1 yes no

δ > θ δ < θ

yes yes no no

False and True and False and True and Ω2 dominates Ω1 dominates Ω2 dominates Ω1 dominates

FIG. 14. Decision tree associated with the inequality in Eq. (S47) that deﬁnes whether Ω1 dominates over Ω2 (or vice-versa) in Eq. (S45).

To estimate the behavior of Eq. (S45) for large values of S, we need to determine the conditions under which Ω1 dominates over Ω2, that is, when θ βθ δ βδ (S47) δ θ 1 β 0 . − + > − + The inequality in Eq. (S47() leads− )( us− to) > the decision tree of Figure 14, and from these results, we can ﬁnd the behavior of Y for large S in Eq. (S45) for each of the following conditions.

Case B.1: β 1 and δ θ

Under these> conditions,< Ω1 dominates the behavior of Eq. (S45) for large S, leading us to cbβ 1 S1 θ βθ Y . (S48) β −1 ln b−S+θ δ ⎛ a ⎞ By comparing the previous= result with Eq.− (S6) and considering that γ − ⎝ ⎠ 1 θ βθ, we ﬁnd = γ 1 − + β 1 for α 1 & β 1 & δ θ . (S49) θ − = + {23 = > < } This result is the same as obtained for case A.4 and ensures the continuity of the behavior of β as a function of α at the point α 1.

Case B.2: β 1 and δ θ =

For this case,< Ω2 dominates< the behavior of Eq. (S45) for large S, leading to caβ 1 S1 δ βδ Y . (S50) 1 −β ln b−S+θ δ ⎛ a ⎞ By comparing the previous= result with Eq.− (S6) and considering that γ − ⎝ ⎠ 1 δ βδ, we ﬁnd γ 1 = β 1 for α 1 & β 1 & δ θ . (S51) − + δ − This result is the= same+ as obtained{ for= case A.1< and ensures< } the continuity of the behavior of β as a function of α at the point α 1.

Case B.3: β 1 and δ θ =

For this case,> Ω2 dominates> the behavior of Eq. (S45) for large S, leading us to caβ 1 S1 δ βδ Y . (S52) β 1− ln b−S+θ δ − ⎛ a ⎞ It is worth noticing that Eq.= (S52) yields positive− values for Y since the − ⎝ ⎠ logarithmic function is negative when δ θ and for large S. By comparing the previous result with Eq. (S6) and considering that γ 1 δ βδ, we ﬁnd > γ 1 β 1 for α 1 & β 1 & δ =θ .− + (S53) δ − This result is the= same+ as obtained{ for= case A.8> and ensures> } the continuity of the behavior of β as a function of α at the point α 1.

Case B.4: β 1 and δ θ =

For this case,< Ω1 dominates> the behavior of Eq. (S45) for large S, leading us to cbβ 1 S1 θ βθ Y . (S54) 1 β− ln b−S+θ δ − ⎛ a ⎞ = − − 24⎝ ⎠ It is worth noticing that Eq. (S54) yields positive values for Y since the logarithmic function is negative when δ θ and for large S. By comparing the previous result with Eq. (S6) and considering that γ 1 θ βθ, we ﬁnd > γ 1 β 1 for α 1 & β 1 & δ =θ −. + (S55) θ − This result is the= same+ as obtained{ for= case A.5< and ensures> } the continuity of the behavior of β as a function of α at the point α 1.

25 2. FUNCTIONAL DEPENDENCE BETWEEN URBAN SCALING AND ZIPF EXPONENTS

At this point, we can combine the results of cases A1-A8 and write the functional dependence of β on α for each of the following conditions.

For γ 1 and δ θ: By combining the cases A.1, A.3 and A.4, we ﬁnd

γ 1 > < 1 θ 0 α 1 β ⎧γ δ 1− 1 δ α 1 α 1 γ 1 . (S56) ⎪ +θ θ < ≤ δ ⎪ + − − ⎪1 γ 1 α 1 γ 1 = ⎨ δ+ − < < δ+ ⎪ − − Equation (S56) represents⎪ + three continuous line≥ + segments: an initial horizon- ⎩ γ 1 tal plateau (α 1) followed by an increasing linear function (1 α 1 δ ) followed by another horizontal plateau (α 1 γ 1) higher than the ﬁrst one.− ≤ δ < < + Figure 15 illustrates the typical behavior of Eq.− (S56). ≥ + > 1, < =

, t n

e )

n 1

o 1 ( + 1 p + = 1 + x

e 1 = = 1 + g n i l a c s = 1 n a b r U

1 = 1 = 1 +

Zipf exponent,

FIG. 15. Illustration of the relation between β and α when γ 1 and δ θ predicted by Eq. (S56).

> <

26 For γ 1 and δ θ: By combining the cases A.5, A.7 and A.8, we ﬁnd

γ 1 > > 1 δ 0 α 1 β ⎧γ θ 1− 1 θ α 1 α 1 γ 1 . (S57) ⎪ +δ δ < ≤ θ ⎪ + − − ⎪1 γ 1 α 1 γ 1 = ⎨ θ+ − < < θ+ ⎪ − − Equation (S57) represents⎪ + three continuous line≥ + segments: an initial horizon- ⎩ γ 1 tal plateau (α 1) followed by an increasing linear function (1 α 1 θ ) followed by another horizontal plateau (α 1 γ 1) higher than the ﬁrst− ≤ θ < < + one. It is worth noticing that we can pass back and− forth from Eq. (S56) to ≥ + Eq. (S57) by replacing δ by θ.

> 1, > =

, t n

e )

n 1

o 1 ( + 1 p + = 1 + x

e 1 = = 1 + g n i l a c s = 1 n a b r U

1 = 1 = 1 +

Zipf exponent,

FIG. 16. Illustration of the relation between β and α when γ 1 and δ θ predicted by Eq. (S57).

> >

27 For γ 1 and δ θ: By combining the cases A.1, A.2 and A.4, we ﬁnd

γ 1 γ 1 < < 1 θ 0 α 1 θ β ⎧γ θ 1− 1 θ α 1 γ 1 α −1 . (S58) ⎪ +δ δ < θ ≤ + ⎪ + − − ⎪1 γ 1 α 1 = ⎨ δ+ − + < < ⎪ − Equation (S58) represents⎪ + three continuous line≥ segments: an initial horizon- γ 1 ⎩ γ 1 tal plateau (α 1 θ ) followed by a decreasing linear function (1 θ α 1) followed by another− horizontal plateau (α 1) lower than the ﬁrst− one. ≤ + + < < < 1, < ≥ =

= 1

, t

n 1

e = 1 + n = o + p

x 1 e

g ( 1 n i l ) a c

s 1 = 1 + n a b r U

1 = 1 + = 1

Zipf exponent,

FIG. 17. Illustration of the relation between β and α when γ 1 and δ θ predicted by Eq. (S58).

< <

28 For γ 1 and δ θ: By combining the cases A.5, A.6 and A.8, we ﬁnd

γ 1 γ 1 < > 1 δ 0 α 1 δ β ⎧γ δ 1− 1 δ α 1 γ 1 α −1 . (S59) ⎪ +θ θ < δ ≤ + ⎪ + − − ⎪1 γ 1 α 1 = ⎨ θ+ − + < < ⎪ − Equation (S59) represents⎪ + three continuous line≥ segments: an initial horizon- γ ⎩1 γ 1 tal plateau (α 1 δ ) followed by a decreasing linear function (1 δ α 1) followed by another− horizontal plateau (α 1) lower than the− ﬁrst ≤ + + < one. It is worth noticing that we can pass back and forth from Eq. (S58) to < ≥ Eq. (S59) by replacing δ by θ.

< 1, > =

= 1

, t

n 1

e = 1 + n = o + p

x 1 e

g ( 1 n i l ) a c

s 1 = 1 + n a b r U

1 = 1 + = 1

Zipf exponent,

FIG. 18. Illustration of the relation between β and α when γ 1 and δ θ predicted by Eq. (S59).

< >

29 For δ θ ϕ: By combining the results of Eqs. (S56), (S57), (S58), and (S59), we ﬁnd = = γ 1 β 1 for α 0 , (S60) ϕ − regardless of γ 1 or γ 1.= Thus,+ we have that> β is constant for all values of α when δ θ. It is also interesting to note that β 1 when γ 1, that is, > < an isometric country scaling implies in isometric urban scaling relationships. = = = For γ 1: By combining the results of Eqs. (S56), (S57), (S58), and (S59), we ﬁnd = β 1 (S61) regardless of γ 1 or γ 1. Thus, constant= returns to scale at the country level imply in constant returns to scale at the city level for all values of α. > <

30 3. NUMERICAL SIMULATIONS OF THE CONNECTION BETWEEN URBAN SCALING AND ZIPF EXPONENTS

In order to understand how variations in the country scaling and urban scaling relationships aﬀect the exact expressions of Sec.2, we have designed a numerical experiment to simulate the empirical relation between β and α. We begin by generating data at the city level, that is, given the total population S, total indicator Y , and Zipf exponent α for a country, we generate a list of m population values

s1, s2, . . . , si, . . . , sm , (S62) where each si is a randomS number= { drawn from a truncated} power-law distribution within the interval smax, smin with exponent α 1 (Eq. S2), and a list of corresponding urban indicators ( ) + y1, y2, . . . , yi, . . . , ym , (S63) where each indicator yi followsY = { the urban scaling} law (Eq. S5) with randon variations

β α,γ,δ,θ 0,σy yi c s 10 i (S64) log y log c( β )α, γ,N δ,( θ )log s 0, σ , i = i y where c is a constant, =β α, γ,+ δ, θ( represents) the+ functionalN ( ) relationship between the urban scaling exponent and the Zipf exponent (deﬁned by Eqs. S56- ( ) S59) for given values of γ, δ and θ), and 0, σy is Gaussian random variable 2 with zero mean and variance σy. The term 0, σy introduces random vari- N ( ) ations into the urban scaling relationship, but has a small impact on the N ( ) country scaling relationships. Thus, in order to account for random variations in the country scaling relationships, we have included analogous terms into Eqs. (S7) and (S8), that is,

δ 0,σδ smin a S 10 (S65) log s log a δNlog( S) 0, σ min = δ and = + + N ( ) θ 0,σθ smax a S 10 , (S66) log s log a θNlog( S) 0, σ max = θ = +31 + N ( ) where 0, σδ and 0, σθ are Gaussian random variables with zero mean, 2 2 and variances σγ and σ , respectively. Finally, we have further considered N ( ) N ( δ ) random variations in the country scaling relationship between Y and S by adding a random term to the total country population, that is,

S S 10 0,σγ , (S67) N ( ) where 0, σγ is a Gaussian random→ variable with zero mean and variance 2 σγ. N ( ) In our simulations the values si and yi should satisfy the constraints m m i 1 si S and i 1 yi Y . We satisfy the first constraint by iteratively m generating values of si and appending to until i 1 si S. To fulfill the = ≈ = ≈ ∑second constraint,∑ we define the values of c by comparing Eq. (S6) with S = ≤ results of Case A (Eqs. S18, S21, S24, S27, S30, S33∑, S36, and S39), that is,

α β α 1 aβ 1 for α β & α 1 & δ θ − ⎧ α β− for α β & α 1 & δ θ ⎪(1 −α)aβ αbα 1 { > > < } ⎪ β−α− − ⎪ α 1 β α for α β & α 1 & δ θ ⎪(α− 1)a b { > < < } ⎪ − ⎪ β α − − ⎪(1 −α)bβ 1 for α β & α 1 & δ θ c α, β, δ, θ, a, b, Y0 Y0 ⎪ { < > < } , ⎪ − − ⎪ α β for α β & α 1 & δ θ ⎪(α− 1)bβ 1 { < < < } ⎪ − ( ) = × ⎨ α −β for α β & α 1 & δ θ ⎪(1 −α)aα 1bβ α { > > > } ⎪ β−α− − ⎪( − ) β α α 1 for α β & α 1 & δ θ ⎪ α 1 a b { > < > } ⎪ − − − ⎪ β α for α β & α 1 & δ θ ⎪(1 −α)aβ 1 { < > > } ⎪ − ⎪ − (S68) ⎪( − ) { < < > } where the value of Y0 is obtained⎩ from the ﬁt of Eq. (S6) to the empirical relationship between Y and S (Fig. 2 of main text). Algorithm1 shows the numerical procedures used in our simulations. This algorithm deﬁnes the function generate country data that takes the country population (S), the power-law exponents (α, γ, δ, and θ), some other constants (Y0, a, and b), and the parameters related to the intensity of the random variations (σy, σγ, σδ, and σθ) as arguments and returns the simulated lists of populations and urban indicators . Having these lists, we estimate the simulated exponent β˜ in the same way we have proceeded with S Y the empirical data, that is, via robust linear regression of the relationship

32 Algorithm 1 Algorithm for generating a list of city populations ( ) and a list of urban indicators ( ) for a country with total population S. S ( , ) = generate country data(S,α,γ,δ,θ,Y0, a, b, σy, σγ, σδ, σθ) S S 10 0,σγ Y µ, σ is Gaussian random variate with mean µ and variance σ2 S Y Nδ( ) 0,σδ smin a S 10 = ▷ N ( ) θ N (0,σθ) smax b S 10 = S˜ 0 N ( ) = = S = [] β β α, γ, δ, θ Deﬁned from Eqs. (S56), (S57), (S58), and (S59) Y = [] c c α, β, δ, θ, a, b, Y0 Deﬁned from Eq. (S68) = ( ) ▷ while S˜ S do = ( ) ▷ si α 1, smin, smax ν, xmin, xmax is power-law random variate < with exponent ν in the interval xmin, xmax = P( + ) ▷ P( ) β 0,σy yi c s 10 i ▷ ( ) N ( ) S˜ S˜ si = si Appends si to the list of city populations = + yi Appends yi to the list of city indicators S ←Ð ▷ S end while Y ←Ð ▷ Y

log yi versus log si. We also obtain the simulated values for total population S˜ and total indicator Y˜ by summing up the values of si and yi. Similarly, we estimate the simulated values ofs ˜min ands ˜max by taking the maximum and minimum values within the list . For our simulations, we use the total population (of all cities within the S power-law regime) of each country in our data set as values of S. For each of these values of S, we also use the corresponding estimated Zipf exponents α. In addition, we set up the values γ, δ, and θ equal to the best fitting values of the empirical relationship between β and α. The constants Y0, a, and b do not affect the relationship betweenα ˜ and β˜, and have been chosen to match the fitted values obtained from the country scaling relationships (with fixed power-laws exponents). We also set the parameters σγ, σδ and σθ equal to the standard deviations of the bootstrap estimates of the empirical exponents γ,

δ, and θ, respectively. Finally, we have varied the parameter σy and obtained the simulated relationship between β˜ andα ˜. The same approach is used when considering that si follows a power-law distribution with exponential

33 cutoﬀ, that is, α 1 p s s exp s sc (S69) −( + ) where sc 0 is an additional( ) ∼ parameter.(− The~ ) only change necessary to implement this case is to replace the truncated power-law random number > α 1, smin, smax in Algorithm1 by an a random number generator associated with Eq. (S69), that is, a α 1, smin, sc . P( + ) P( + )

34 4. ZIPF AND URBAN SCALING PLOTS FOR EACH COUNTRY

In this section, we show Zipf’s law and urban scaling law adjusted for every country in our data set. In what follows, each row of panels shows the results for a given country (name is indicated within the plots). The ﬁrst column shows the complementary cumulative distribution F s of city population s (gray continuous line) and the adjusted power-law distribution ( ) (dashed line). The second column shows the rank plot, where each city is named within the plots. The dashed lines show the adjusted Zipf’s law. The power-law exponent α ( standard error) and the lower cutoﬀ population size ssmin ( standard error) are estimated from the approach of Clauset- ± Shalizi-Newman [SIAM Review 51, 661 (2009)]. The third column shows ± urban scaling between urban GDP and city population, where each city is also named within the plots. The dashed line represents the adjusted scaling law. The urban scaling exponent β ( standard error) is estimated via robust linear regression on the log-transformed quantities. ±

1.5 Qala i Naw Chaman 0.0 Pul-i-Alam Bagram Afghanistan (AFG) Maymana Afghanistan (AFG) Afghanistan (AFG) KhanabadKhost Kabul Mihtarlam 10 Samangan Farah Girishk = 1.14±0.22 Sheberghan Shindand = 1.14±0.22 = 1.48±0.27 ) s = 78159 Ghazni Sar-e Pol s = 78159 s min Guzarah Charikar min ( 1.0 Taloqan Lashkar Gah y r

F 0.5 Kandahar KunduzPuli Khumri Jalalabad g g 9

g Mazar i sharif TaloqanKunduz o o Jalalabad Khost o l l Mazar i sharif Sheberghan

l Herat

Puli Khumri

, , Khanabad ,

k 0.5 Girishk P n Kandahar Lashkar Gah Pul-i-Alam Sar-e Pol 1.0 n

o Herat Mihtarlam D i a 8 Farah t Bagram r G

u Samangan ShindandCharikar

y Guzarah b y t

i Ghazni t i r i

t 0.0 Kabul C

1.5 C s i 7 Maymana d

Chaman

Complementary cumulative 0.5 2.0 Qala i Naw

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s 11 Balombo Caxito Bocoio Chinguar 0.0 Luau Vila Nova do Seles Luanda Porto Amboim UNNAMED Angola (AGO) Catabola Cafunfo Angola (AGO) Angola (AGO) Cazombo Katchiungo Quipungo Waku Kungo 1.5 Soyo Gabela Caconda = 0.92±0.14 Camacupa M'banza Congo = 0.92±0.14 Dundo = 1.18±0.06

) N'dalatando Sumbe Calilongue Ganda s smin = 62880 Andulo smin = 62880 Bailundo Matala ( Cubal Namibe y r F 0.5 Menongue Ondjiva 10 g g Luena Saurimo g 1.0 Huambo Lubango o o Kuito Uíge

o Benguela l l

l Lobito

Kuito MalanjeCabinda Lobito Luena Uíge , , Malanje Bailundo , Menongue Saurimo

k Cubal P

n Namibe Benguela Andulo n Cabinda N'dalatando

o Sumbe D

i 1.0 Dundo Matala a 9 Calilongue t

r Camacupa UNNAMEDCafunfo M'banza Congo G

u 0.5 Huambo Porto Amboim Caxito Gabela

y Vila NovaKatchiungo do Seles Ondjiva

b Waku Kungo

y Soyo t

i Lubango Luau t i Cazombo

r Caconda i Chinguar t Balombo C C s i 1.5 Quipungo Ganda d

8 0.0 Luanda Catabola

Complementary cumulative Bocoio

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

35 Concepción Buenos Aires Rio Grande Dassa-ZouméSan Francisco VenadoTindji Tuerto Tartagal EldoradoNecochea 0.0 OlavarríaViedma SavalouJunín GualeguaychuAllada Goya ArgentinaBenin (BEN) (ARG) Pergamino San Ramón de la Nueva Orán ArgentinaBenin (BEN) (ARG) ArgentinaBenin (BEN) (ARG) Oberá Puerto Madryn TanguietaRafaela VillaBanikoara María 11.5 Villa Carlos Paz General Roca Trelew Luján Cotonou Presidencia Roque Sáenz Peña RíoNikki Gallegos SantaLokossa Rosa San Carlos de Bariloche ==0.19.32±1±0.01.227 Reconquista Tandil ==0.19.32±1±0.01.227 9.5 Villa Mercedes Río Cuarto = 1.0218±0.0163 1.5 San MartínKétou SanPobè Nicolás de los Arroyos ) ) Concordia Zárate 1.0 San Rafael San Luis smsimn i=n =6557241627 NatitingouLa Rioja Malanville smsimn i=n =6557241627 s s Comodoro Rivadavia San Fernando del Valle de Catamarca Paraná ( ( Formosa 0.5 WhydahBahía Blanca San Salvador de Jujuy 11.0 Porto-Novo y y r r Kandi F F Posadas Corrientes 0.5 Neuquén Resistencia g g g g Djougou 9.0 Parakou g g Santiago del Estero Santa FeBohicon o o o o 1.0 Salta o o San Juan l l l l Rosario Cordoba l l

10.5

Mar del Plata , , , , La Plata , , 0.5 Parakou Mendoza k 1.0 k P P n n San Miguel de Tucumán Mendoza Porto-Novo n n Mar del Plata La Plata o o 8.5 San Miguel de Tucumán D D i i a 1.0 a DjougouKandi Santa Fe t t Bohicon r r 0.5 10.0 Banikoara Neuquén G Rosario G

u u Río Gallegos Malanville

Comodoro Rivadavia San FernandoSan del JuanValle de Catamarca WhydahBahía Blanca Salta y y Cordoba LokossaKétou Pobè b b y y Nikki t t i i t t i i Trelew San RafaelSantiagoSan Luisdel Estero r r i 1.5 i 0.0 Cotonou PuertoNatitingou Madryn PosadasSanCorrientes SalvadorResistencia de Jujuy t t Paraná Villa Carlos Paz Río Cuarto C C Reconquista Zárate 8.0 VillaSan Mercedes Martín San Nicolás de los Arroyos C C Santa Rosa s s 9.5 La Rioja Rio Grande Tandil i i Rafaela Concordia Formosa General Roca Villa María San Francisco d d 1.5 0.0 Buenos Aires Pergamino Luján AlladaViedma

Necochea OlavarríaTanguietaVenadoSavalouJunín Tuerto Gualeguaychu Presidencia Roque SáenzOberá Peña Dassa-Zoumé 2.0 9.0 Tartagal San Ramón de la Nueva Orán Goya Complementary cumulative Complementary cumulative 7.5 Eldorado 0.5 San Carlos de Bariloche 0.5 ConcepciónTindji 5.0 5.0 5.5 6.05.5 6.5 6.0 7.0 5.0 5.0 5.5 6.05.5 6.5 6.0 7.0 5.0 5.0 5.5 6.05.5 6.5 6.0 7.0 City population, log s City population, log s City population, log s

BoulsaRockingham Gaoua Melbourne Kombissiri Ziniaré 0.0 1.5 ToowoombaSolenzo KongoussiCairns Ouagadougou BurkinaAustralia Faso (AUS) (BFA) Leo Niangoloko BurkinaAustralia Faso (AUS) (BFA) 11.25 BurkinaAustralia Faso (AUS) (BFA) Sydney Tougan Nouna 1.2 AitkenvaleKoupéla Manni Djibo BitouBelconnen Yako Gourcy 9.5 ==10.6.762±±00.3.107 SidéradougouHobart Woden ==1.06.67±2±0.03.017 ==11.2.070±±00.1.02 Bogandé Pouytenga ) ) Houndé Dédougou ssmminin==51605702856 Campbelltown MaroochydoreBanfora smsminin==56150078256 s s 0.50.2 1.0 Dori 11.00 BrisbanePerth ( ( 1.0 Kaya TenkodogoGeelong y y r r Wollongong F F Bobo-Dioulasso Ouahigouya Koudougou 9.0 g g g g g g Newcastle Adelaide Gold Coast

o 10.75 o o 0.4 o 0.5 Fada N'Gourma o o 1.0 l

0.8 l l l l l Bobo-Dioulasso

, , , , , , Adelaide k k 8.5 P P n n Perth n n 10.50 Dori o o 0.6 D 0.6 D i i 1.5 0.0 Tenkodogo Koudougou a a Ouagadougou Fada N'Gourma t t OuahigouyaBanforaNewcastle Gold Coast r r Houndé G G

u u

Brisbane Pouytenga y y b b 8.0 y y t t Kaya i i Gaoua t t i i 10.25 NiangolokoNouna Dédougou r r i 0.4 i 0.8 Manni t t 2.0 0.5 Gourcy Sydney Leo C C KoupélaBoulsaYako Woden Wollongong C C s s Solenzo KongoussiBitou i i Ziniaré GeelongBelconnen Campbelltown Maroochydore Kombissiri d d 7.5

10.00 1.0 0.2 Djibo Sidéradougou 2.5 1.0 AitkenvaleRockinghamHobart Tougan Complementary cumulative Complementary cumulative Cairns Toowoomba Bogandé 1.2 0.0 Melbourne 9.757.0 5.0 5.0 5.5 5.5 6.0 6.0 6.5 6.5 5.0 5.5 5.5 6.0 6.0 6.5 6.5 5.0 5.5 5.5 6.0 6.0 6.5 6.5 CCityit yp pooppuulalatiotionn, ,l ologgss CCitiyty p pooppuulalatitoionn, ,l ologgss CCiittyy ppooppuullaattiioonn,, llooggss BakuDhaka 0.0 ShamakhiRajshahi KurdamirMymensingh BangladeshAzerbaijan (BGD)(AZE) SabirabadLakshmipurBardaSylhet BangladeshAzerbaijan (AZE)(BGD) 11.00 AzerbaijanBangladesh (AZE) (BGD) 1.0 Masall Sheki 10.5 ==1.12.73±4±0.03.236 Khulna Tangail ==11.2.374±±00.3.326 = 1.04±0.14 1.0 C lilabad irvan = 1.20±0.08 ) ) Dhaka Dhaka smsminin==10635664958 ssmminin==160336546598 s s NakhchivanBogura 10.75 ( ( MingachevirBrahmanbaria

0.5 y y r r F F 0.5 Lankaran Sirajganj g g g g g g 0.5 10.0 Dhaka Ganja 10.50 o o o o Chattogram o o l l l l l l

Ganja

, , , , , , 0.5 Chattogram Comilla k 1.0 k P P n n n n Comilla 10.25 o o 0.0 D Baku D i i a a 9.5 t t Dhaka r r Nakhchivan G G

Sabirabad u u

1.0 irvan y y b b y y t t 10.00 Sirajganj i i Dhaka t t i i Mingachevir r r i 1.5 i Kurdamir Brahmanbaria Shamakhi Lankaran t t C C 0.0 Sheki BoguraDhaka C Dhaka C s s 0.5 i i 9.0 Khulna d d 9.75

Masall Barda

Rajshahi Sylhet 1.52.0 Tangail Complementary cumulative Complementary cumulative 1.0 9.50 LakshmipurC lilabadMymensingh 0.5 8.5 6.00 5.06.25 6.505.5 6.75 6.07.00 7.256.5 6.00 5.06.25 6.505.5 6.75 6.07.00 7.256.5 6.00 5.06.25 6.505.5 6.75 6.07.00 7.256.5 Ciity popullatiion, log s CCiittyy ppooppuullaattiioonn,, llooggss CCitiyty p pooppuulalatitoionn, ,l ologgss

0.000.0 1.0 NyanzaZhodzina Lac 9.0 Minsk BurundiBelarus (BDI) (BLR) KayogoroOrsha Makamba BurundiBelarus (BDI) (BLR) BelarusBurundi (BLR) (BDI) 1.0 Barysaw 10.5 Bujumbura Kirundo Pinsk Rumonge ==1.606.9±50±.500.26 Baranovichi ==10.6.965±±00.5.206 ==11.0.606±±00.0.337 ) ) smsin ==8196094509 Kayanza sms in ==8910649059

s min min s Babruysk Muyinga

( 0.25 ( 0.8 Brest 8.5 y r y r F F 0.5 0.5 Hrodna g g g g Ngozi 10.0 g g Homel Mahilyow o o Gitega o o o o HrodnaVitebsk l l Gitega l l l

l 0.6

0.50

Mahilyow

, ,

, Brest ,

, Ngozi , 8.0 Kayanza k k P P n Muyinga n Vitebsk Babruysk

n Kirundo n o o D D i i Baranovichi

a 0.4 a t t r 0.751.0 r Homel 9.5 PinskBarysaw G

0.0 G

Bujumbura u u

y y b b y

y 7.5 t t i i t i t i r i r 0.2 i Zhodzina t t C C C C s s i i 1.00 NyanzaMakamba Lac Rumonge d d

9.0 1.5 0.50.0 Minsk 7.0 Complementary cumulative Complementary cumulative 1.25 Kayogoro 0.2 Orsha 5.05.0 5.2 5.55.4 5.6 6.05.8 6.0 5.005.0 5.255.2 5.505.4 5.755.6 6.005.8 6.256.0 5.00 5.05.255.25.505.4 5.755.6 6.005.8 6.256.0 CCitiyty p pooppuulalatitoionn, ,l ologgss CCitityy ppooppuulalattiioonn,, llooggss City population, log s 1.0 Riberalta Santa CruzBrussels de la Sierra 0.0 Bolivia (BOL) 1.0 Bruges Bolivia (BOL) Bolivia (BOL) 0.0 Leuven El Alto [La Paz] Belgium (BEL) Trinidad Belgium (BEL) 10.00 Belgium (BEL) = 0.80±0.24 Montero = 0.80±0.24 = 1.07±0.08 = 0.78±0.25 Mons = 0.78±0.25 11.0 = 1.09±0.10 Cochabamba ) ) smin = 84248 Potosí Kortrijksmin = 84248 s s 0.2 smin = 179108 0.8 Tarija smin = 179108 Antwerp ( ( 0.2 0.8 9.75 y r y r F F Ghent 10.8 g g Sucre g g g g o o Charleroi o o o Oruro l l o l l l 0.4

0.6 l 9.50

0.6 , , 0.4 , , , ,

k 10.6 Liège k P Tarija n Oruro P

n Ghent n

n Liège o Cochabamba D o i D a i a t t 0.6 r 9.25 r G

Charleroi

0.4 G u

0.4

u 0.6 10.4 Kortrijk y y b y t b Santa Cruz deAntwerp la Sierra y i t t i i t i r i r i

t Potosí

t Leuven C

C 9.00 C s C

s 0.8 i i 10.2 0.20.2 Montero d 0.8 d

8.75 Bruges Trinidad Complementary cumulative 1.0 10.0 Sucre Complementary cumulative Mons 1.0 0.00.0 El Alto [La Paz] Brussels Riberalta 5.255.0 5.505.5 5.75 6.0 6.00 6.25 5.0 5.6 5.55.8 6.06.0 6.2 5.0 5.6 5.55.8 6.06.0 6.2 City poCpituyl aptoiopnu,l alotigons, log s City pCopituyl aptoiopnu,l alotiogns, log s City pCopityu lpaotiopnu,la ltoiogns, log s

36 Santa Cruz do Sul São Paulo CatanduvaVarginha Sinop ItapetiningaBotucatu Santa Cruz do CapibaribeJaú AtibaiaIlhéus AngraPatos dosde Minas Reis JequiéBento Gonçalves Vitória de SantoAlagoinhas Antão ResendeBragança Paulista 0.0 PousoParanaguá Alegre GuarapuavaPoços de Caldas Teixeira Arapongasde Freitas Araruama Brazil (BRA) BarreirasLages TeresópolisAraguaína Brazil (BRA) Brazil (BRA) Parnaíba Nova Friburgo SobradinhoSobral PlanaltinaAraçatuba CastanhalPalmas CachoeiroÁguas Lindas de Itapemirim de Goiás PassoRio GrandeFundo DouradosCamaçari RioChapecó Verde JaraguáRondonópolis do Sul 2.0 ParauapebasItabuna SantarémMogi Guaçu RioMaraba Claro IndaiatubaArapiraca Sete LagoasMarília MacaéRio das Ostras = 0.92±0.08 PresidenteDivinópolis Prudente GuaratinguetáGovernador Valadares = 0.92±0.08 São Carlos Santa Maria = 1.00±0.06 Ciudad del Este Imperatriz Rio de Janeiro AraraquaraMossoró Cabo Frio Vitória da Conquista UberabaCaruaru ) Cascavel Pelotas CriciúmaLimeira Ponta Grossa s = 127971 Boa Vista AnápolisRio Branco s = 127971 min FrancaBauru Juazeiro do Norte min s Montes Claros Campos dos Goytacazes Piracicaba Petrolina Porto Velho Campina Grande Belo Horizonte

( Volta Redonda 11 Blumenau Ipatinga São José do Rio Preto Taubaté Juiz de Fora Campinas Caxias do Sul Macapá y 0.5 r Maringá Feira de Santana F Joinville Gama Porto Alegre Itajaí Londrina SaoSantos Goncalo Goiânia Uberlândia Novo Hamburgo Ribeirao Preto Ceilândia g g 1.5 Cuiabá Campo Grande Fortaleza g SorocabaFlorianópolis Vila Velha Recife Aracaju Sao Jose dos Campos Sorocaba o o João Pessoa Sao Jose dos BelémCamposManaus o Teresina Ribeirao Preto l l Natal

l Maceió Florianópolis

Vila Velha Itajaí Novo Hamburgo São Luís TaubatéJoinville

, São José do Rio Preto , Gama PiracicabaCaxias do Sul , Santos Ceilândia Uberlândia FrancaBauruBlumenau VoltaCampo Redonda GrandeAracaju Natal

k Sao Goncalo Manaus Limeira Campos dosJuiz Goytacazes de Fora CuiabáMaceióSão Luís

P 10 Anápolis São Carlos n 1.0 Criciúma Ipatinga João Pessoa Presidente PrudenteAraraquaraMarília Montes Claros Goiânia Rio Claro Cabo Frio n Belém IndaiatubaGuaratinguetáPelotas Curitiba

o Mogi Guaçu RioSanta dasPorto Ostras Maria Velho D Araçatuba Macaé Uberaba i 1.0 RioChapecó Verde Jaraguá do Sul Macapá Campinas Porto Alegre SobradinhoPalmas a PassoRio GrandeFundo LagesBragançaSeteNovaDivinópolis Paulista CachoeiroLagoas ÁguasFriburgo Lindas deBoaGovernador Itapemirim de Vista Goiás Valadares t Jaú AraruamaAraguaínaPlanaltina Rio BrancoCampina Grande r ItapetiningaAtibaia Teresópolis CaruaruJuazeiro do Norte Angra dos ReisCatanduvaBotucatuBentoResende Gonçalves G Mossoró Curitiba Santa Cruz do Sul

u Salvador Salvador Pouso Alegre ParauapebasPoços de CaldasDourados Patos de MinasVarginha RondonópolisSantarém Petrolina y Maraba

b Castanhal y Imperatriz t Vitória de Santo Antão Arapiraca i Fortaleza t i Santa Cruz do CapibaribeSinop Recife Sobral Londrina r i Maringá Teresina

t 1.5 Ciudad del Este

C 9 C s Cascavel Ponta Grossa i 0.5 Belo Horizonte Feira de Santana

d Rio de Janeiro

Arapongas Guarapuava Paranaguá Vitória da Conquista Itabuna Camaçari Barreiras Teixeira de Freitas 2.0 AlagoinhasIlhéus Complementary cumulative ParnaíbaJequié 0.0 São Paulo 8 5.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s Toronto Thunder Bay Peterborough 0.00 1.50 Nanaimo Brantford 11.5 Canada (CAN) Red Deer Kingston Canada (CAN) Canada (CAN) Moncton Saint-Jean-sur-Richelieu Saint-Jérôme Guelph = 0.82±0.14 Trois-Rivières Abbotsford = 0.82±0.14 = 1.01±0.02 St. John's Sherbrooke ) smin = 103243 1.25 Regina Barrie smin = 103243 s 0.25 Montreal Saskatoon St. Catharines ( y r DetroitVictoria Vancouver F 11.0 Halifax London g g

g Calgary 0.50 1.00 Edmonton o o Kitchener Winnipeg o l l Ottawa l

, , Quebec , Edmonton k P n Quebec

0.75 n 0.75 o

Ottawa D Winnipeg i 10.5 a

t Calgary Kitchener r G

London u

y b

y Detroit t i 1.00 0.50 t i Vancouver St. Catharines r i Halifax Saskatoon Victoria t C C s Montreal i Regina Barrie 10.0 St. John's

d Red Deer Abbotsford Trois-Rivières Sherbrooke 1.25 0.25 Guelph Kingston Saint-JérômeSaint-Jean-sur-Richelieu Nanaimo Brantford Moncton Peterborough Complementary cumulative Thunder Bay 1.50 0.00 Toronto 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s Zurich 0.00 Switzerland (CHE) Thun Switzerland (CHE) Switzerland (CHE) 1.00 Montreux Biel/Bienne = 1.21±0.34 Lugano Olten = 1.21±0.34 11.00 = 1.13±0.05 ) smin = 132808 smin = 132808 s 0.25 Sankt Gallen Lucerne ( y r

F 0.75 Lausanne 10.75 g g

g Bern o 0.50 o o l l l

Geneva

, , , 0.50 Basel k P n 10.50 Basel n o

0.75 Geneva D i

a Lausanne t r G

u 0.25 Bern

y b y t

i 10.25 t i r 1.00 i Lucerne t C C s Sankt Gallen i 0.00 Zurich Olten d

10.00 1.25 LuganoBiel/Bienne MontreuxThun

Complementary cumulative 0.25 1.50 9.75 5.25 5.50 5.75 6.00 6.25 5.25 5.50 5.75 6.00 6.25 5.25 5.50 5.75 6.00 6.25 City population, log s City population, log s City population, log s Santiago 0.0 Curicó Chile (CHL) 1.0 Chillán Antofagasta Chile (CHL) Chile (CHL) = 1.17±0.32 Puerto Montt Iquique = 1.17±0.32 = 1.21±0.06 Talca ) Rancagua smin = 175987 smin = 175987 11.0 s

( Temuco La Serena y r F 0.5 g g g 0.5 Concepción o o o l l l

Viña del Mar [Valparaíso] 10.5

, , , k P n n o D

i 1.0 Viña del Mar [Valparaíso] a t r G 0.0 u Santiago

y Concepción b y t 10.0 i t i r i

t La Serena C C

s Rancagua i 1.5 d 0.5 Talca Temuco Puerto Montt 9.5 Chillán Antofagasta Complementary cumulative Curicó Iquique

5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s

HuozhouHunchun Shanghai LaotingQinyangQuyangBeipiaoHuaijiCaijiapoJJiawangWangchengnchéngji ng 3.0 XiaochangQingShe FengqiuHeishanXiayiHengxianLintaoWuchangLuyiDingyuan FengduWuyangLiuzhi FengchengYunxiaoDong'e TianchangPuchengJiulongbaLuannanZhijiangShangliLingbiJiutai Gong'anNanchuanJiaTongchengWuheDonggangXianJuDiaobingshanZhong JingbianZhangshuDongxianDeqingJintangLinShuiXuyiYu YuanyangHengyangLiulinWulianXilinhotLongxiYanqingDengfengAnping 12.0 YuanpingRongxianJinchangGaoyiWen KuqaLinglingZhuangheXinxingShizuishanMengcheng GuangshuiGongnongGaochengYichengJinxiangTaojiangQiyangXinminShaxiBaidi YanglingQianxiNeihuangDehuiTongrenYuanjiangDehua DanjiangkouDanchengShangluoDengfengYuchengHong'anYiyangBiyangLujiang LianyuanShayuanTangyinShifangQionghaiTaikangLijiangYingshanWugang LongshanZhaoyuanYunyangLianjiangXingguoNingduYong'anPu'er LaibinNingguoLushanXincaiLongchangTongchuanXinfengYongji ChongmingLongchuanDangyangMianxianFudingGanguYanlingZijin YingdeJiaziKaizhouNanxianSongziGuantaoWenshangPulandian FangchengEmeishanYongchunHanchengGuoyangXinzhouChangqiGu'anDazhu HuairenGejiuNanbuShuchengJieshouSuiQianjiang 0.0 RongchengYingchengChizhouXiangheTiantaiJinxian XipingZhongweiLantianLangzhongZhongxiangFu'anHuarong JianshuiWujiXinle YunyangMengziXuwenDalad County ShengfangLiangmataiQixianSanheBaodeFunanJuye RenshouDongpingLixinJinzhouHuangpiLingyuanTongliangYangchun XingchengLvsigangFushunChangliYan'anRunanYuganAnlu YushuRudongZhechengYichuanShuangchengPoyangXiangshanZaoyang China (CHN) RongchangShenmuFeixianTongluDafengJieshiChibi MengzhouRushanYiyuanQuxianLaolingBoxingDuyun China (CHN) China (CHN) JuanchengJiayuguanFenyangLiuyangFula'erjiYima BaiyinDianjiangJiangshanSanchaheShulanZhangjiajieHuazhou GongzhulingXiangxiangShenxianGaochunDunhuaZhaoCi LinzhouXinyiChipingHejianAnda GuannanMachengJizhouNingjinWeishiTang BazhongAnxiYushanGuipingFengshun ShaodongDanzhouNanpingJishouEnpingXinshiYa'an XinqiaoZunhuaQionglaiQingyangQian'anWuxueBazhou HuangchuanWenshanSheyangLuchengJinghaiJiexiuLili UlanQijiangDangshanHamiWeiBobai Hot ZhongmuGuanyunHulunbuirPingyaoAnningLishui YongzhouLhasaLuanzhouLingbaoChangshouJiuquan Beijing Guang'anLianshuiPingtanAnji YunfuWuhuaGaopingHuangshan YunmengZhangyeZhenpingQingheGaoyouQidongYintai NingdeRuijinQingtianGaoanLin'anEnshiFengxian JianyangHanchuanKaramayYangguSihongYutianHotan HezhouTianmenBotouSanmingLingshan PingliangLuodingXiangyinXiaoyiBijieDayi AchengMiyunBinyangMancheng FenghuaSuihuaJurongLüliangTaiheXinji RuzhouBaishanDonghaiHainingXuanchengHongtong MeihekouXuanweiLinquanGanyuHuyi XinghuaLinxiaNanxunXinyiKaiyuanBinhai YizhenNankangFuningGushiCenxi ZhaodongWendengLishuiQinzhouXianningHebi BaichengNingxiangChengdeJiangyouZiyangKaili XingyiPingguLanxiLilingBaoying WafangdianXinzhouHoumaTaishanXushuiAksu UlanqabLianjiangLongkouNinghaiBayannur GaozhouLaizhouSiyangDongtai HuairouShuangyashanZhuozhouLufengCaoxian HuaiyangChangyuanLuchengPinghuHaimenXinjin JianhuTangheShangcaiLeping 1.07 0.04 ShanchengLeizhouChaohuRugaoWu'an ShuozhouQidongChangxingShanxianHaianTieling 1.07 0.04 = ± Pujiang = ± GuangfengWuchuanChanggeLinqing HangyuluWuweiDengzhouXichang = 1.03±0.01 FengfengFengrunYingtanPingdu SuiningConghuaHechuanJi'an GuchengXuanhuaBozhouJuxianJiyuan RenqiuRanghuluShanweiSuizhou TianshanHuaxianAnshunLouxing XinhuaDingzhouLinhaiDali DashiqiaoChuzhouHeyuanJiangjinLeiyang ShiheziWeinanGuigangJingmenYining/Qulja ) XiaoxiangJintan GaobeidianYuncheng Changji/SanjiYongchengXinzhengBaoshanLinying LiyangWuzhongBinzhouHengshuiJingdezhen YongchuanLiupanshuiMeishanXianyouSuining PizhouXingningGaomiFuqing 2.5 ChenzhouZhaotongLeshanYuxi SanmenxiaLinquSanyaXiantaoXinxing Dujiangyan DongshengLiaoyuanLongyanTonghuaHaiyan ChaoyangFuyangNeijiangXintai 11.5 WuzhouYichun Shuyang ShouguangHuaihuaJinzhongDaqing DanyangXitaoHuludaoZhoushan s = 166393 ZhangqiuXinyangBeibeiYanji YuhuanDazhouSipingXinyu s = 166393 Chengdu AnkangQuzhouLixian ZhumadianShengzhouJixi s min Shangrao Yangzhong min HaichengSuzhouDeyang KaipingKashgarMeizhou TianjinSuzhou JinchengYingkouQujingLu'an TonglingZuchengXinmi Hangzhou SongyuanWenshuiQingzhouHegang FuzhouPukouYiyangXiangcheng PanzhihuaYangjiangTianshuiTongliao SuqianLaiwuHanzhongZigong XiaoganXingpingTaixingZhuji HuzhouWeihaiYibin

( 0.5 Anqing Zhoukou YuzhouBeihai YongnianRizhaoEzhou ZouchengLaiyangShiyanYixing DandongBayuquanLiaocheng LiaoyangChifeng YangquanJiamusiZhaoqing Nanjing Guangzhou ZhangjiakouShaoguan ChangzhiQingyuanJinhua LangfangWanzhouDezhou PanjinMudanjiangZhangzhou

Linfen y r Gongyi JiujiangNanchong LianyungangFuxinHeze TengzhouBenxiChangde F Ma'anshanCangzhou PuyangKaifeng ZhenjiangJingzhouLuzhou PingxiangTai'anPingdingshan YueyangShaoyangQiqihar XuchangMaomingQinhuangdao Shenyang GanzhouYichangJining HuangdaoYongkangXingtai ChangshuZunyi TaizhouJinzhou MianyangFuyang BaojiLuoheZhanjiang WenzhouWuhan Jieyang ShangqiuYanchengZhuzhou XiangtanNanyang ChongqingQuanzhou g g BengbuHuaibei GuilinYueqing Xiamen JiangyinYulin ZhongshanFancheng g 2.0 HengyangWuhu Jiaozuo Ningbo Guangzhou 11.0 Zhangjiagang HuangshiXining Anyang QingdaoXi'an FushunPutian Tianjia an Cixi Changzhou YangzhouHuai'an JilinHandan Dalian Jinan Guangzhou o o Zhuhai Yiwu YinchuanBaoding LiuzhouAnshan TaizhouZhengzhouFuzhou o Datong HaikouXinxiang Binhai Zibo l l Weifang Jiaxing Changchun GuangzhouBinhai HuizhouHohhot Guangzhou l Yantai Baotou ZhuhaiNantong Shijiazhuang

1.0 Guangzhou Luoyang Xuzhou Nantong Xuzhou Harbin TangshanLinyi Lanzhou Linyi GuangzhouHefei

Taizhou , , JiaxingNanning Guiyang , Changzhou Ürümqi Taiyuan Cixi Kunming YangzhouHuai'an YantaiYiwu NanchangChangsha GuangzhouÜrümqi Zibo Guangzhou

k Taiyuan Dalian Changsha WeifangZhangjiagangAnshan Ningbo P Jiangyin Tangshan Lanzhou n Changchun Fuzhou Yancheng FushunYueqing Baotou 1.5 Nanchang Jinan 10.5 Putian HohhotNanning Changshu Taizhou n Shijiazhuang Harbin Luoyang Guangzhou YongkangYinchuanGuangzhouJilin Kunming o Qingdao ZhenjiangMianyang JinzhouBaodingHandan Hefei Wenzhou D Lianyungang GuangzhouDatong Xinxiang i Jining Yixing Huizhou a Xiamen Jinhua HuangdaoHuangshi Chongqing HezeLuzhou Tai'an Haikou Guiyang t Zhengzhou Tengzhou XiningFanchengAnyang TaixingZhuji WanzhouSuqianFuxin NanchongBenxi

r Quanzhou Panjin Liuzhou LangfangHuzhouLiaocheng Hengyang 1.5 Shenyang ShuyangLiaoyangDezhouBayuquanZhangzhouGuangzhou JiaozuoTianjia an YangzhongLaiyang Baoji G Quzhou Zoucheng Rizhao Zhuzhou Nanyang Guangzhou ZigongJingzhou Xingtai PizhouPanzhihuaFuyangZhoushanShengzhouPukouYibin YichangShangqiuQinhuangdaoWuhuLuoheXiangtanZhongshan u Suining YingkouDanyang Yulin Xi'an ZhangqiuQingzhouYuhuan Laiwu HuaibeiZhanjiang Huairou Zucheng HaiyanHaichengBeibeiDeyangXintai ChifengWeihaiCangzhouDandongQiqihar PingdingshanGuilin Tianjin Jintan YueyangShaoyangFuyangBengbuXuchang Pinggu MeishanShouguangLiyang Nanjing YongchuanQidongLeshanLinhai DujiangyanNeijiangHuludao Gongyi Kaifeng

y Maoming Haian BinzhouGaomiChaoyangZhangjiakouLinfenGanzhouMudanjiangPuyang Jiangjin SongyuanLinquXingping ChangzhiChangde b Jinghai Miyun Longyan TongliaoDazhouShiyan EzhouYangquanJiamusi Guangzhou y Funing BaoyingJianhuPujiangSuiningDongshengYanji Siping Ma'anshan Pingxiang t 10.0 SiyangHaimen DashiqiaoChangxingXianyou Fuqing HanzhongYongnian 1.0 Wuhan DongtaiPinghuXinjin Hechuan TianshuiXiaogan Jiujiang i FengxianXinyiRugaoJuxian HegangZhoukou Jizhou Pingdu LiaoyuanXichangTonghua Shaoguan Qingyuan

t Ninghai i SihongFenghuaGanyu DonghaiBinhaiLanxi YuzhouXinmi Zhaoqing Yizhen Xinghua Tieling JinzhongWenshuiXiantao Chongming Sheyang Jurong Nanxun Caoxian HuaihuaAnkangJinchengJixi YiyangXiangcheng r i Guannan LaizhouHainingLishui Suzhou Guanyun Jiangyou ShanxianWeinanDaqing Beihai Zunyi Hangzhou LianshuiGaoyou WafangdianZiyang Tianshan JingmenGaobeidianXinyangWuzhongLixianYangjiangLu'anZhumadian Yanqing Lili Xinqiao KaiyuanFengrunLongkou Shihezi Kashgar Anqing RudongLishuiAnji Qingtian WuweiDingzhou XitaoMeizhouFuzhou t Lin'an Xinxing Lvsigang LingyuanJianyangYanggu WendengLinqing RenqiuSuizhouHengshui Dafeng Dayi FengfengXinzheng ShangraoSanmenxiaKaiping Gaochun Changshou UlanqabJiyuanXiaoxiangLeiyangRanghuluYichunSuzhou Tongling Tonglu Guang'anJiangshanChiping GuchengXuanhua WuzhouYuncheng Xinyu C Qionglai Chenzhou Jingdezhen FeixianJuye XiangshanYa'anLaolingPingtanQijiang Sanming Wu'anBozhouZhuozhouYongchengLouxingLinyingHangyuluDengzhouYining/QuljaSanya Baicheng BayannurChangji/SanjiHuaxianShangcai Xingning YuchengTiantai Zhenping NingdeHouma ChanggeChuzhou Xinhua C DengfengRongchangFushun Meihekou Huaiyang Heyuan s Chengdu Shaxi ShenxianQuxian Shuozhou Jieyang YunyangDongping Huyi DazhuJuanchengKaizhouNanbu Zhangye ShanweiGuigang Xuyi LongchangYiyuanBoxingDianjiangAnxi ChengdeTianmenXinzhouXushuiChangyuanShanchengLinxiaXianning Conghua Qujing Deqing EmeishanXingchengGongzhulingLangzhongTongliangRushan Xinji ShuangyashanTanghe i YongchunShifang NanpingPingliangUlan HotBotou ZhaodongYingtan Ji'an 2.0 LinShuiJiawangZhuangheYan'anPulandianWeishiHanchuanLingbao Guangfeng JiulongbaJinxiangRongchengXianJu RenshouHulunbuirYintai AchengManchengBaishanChaohuWuchuanLiling YunmengShulanYutianZunhua GaozhouTaishan Yuxi FengchengYunxiaoDong'eGu'an AnluWenshang BazhongSuihuaWeiJiuquanGaopingNingxiangGushi LianjiangZhaotong XinminLianjiangShengfangDehuaJiayuguanFu'an QingheQingyangLuanzhouAksuHongtongLeizhouLufengLepingBaoshan JintangZhaoyuanHeishanYong'anFuding Lucheng Linquan Hebi DiaobingshanZhongChangli ZhongmuXiaoyiBazhouNankangRuzhou GongnongNanchuanDonggangWulian TangZaoyangBaiyinHejianQian'anHamiLüliang Qinzhou Beipiao ChibiDunhuaNingjinZhaoYushuJiexiuQidong Dali Fengdu LiangmataiSanheHuangchuanZhongweiLantianJishouCiMengzhouPingyaoXiangyinHotan EnshiLingshan d ShenmuSongziGuantaoMachengYimaQianjiang AndaLuodingWuxueTaihe WenHong'anYanlingXiangheWujiFenyangJinzhouXinshi DangshanWuhuaRuijinHezhouCenxiBinyangXuancheng Liupanshui BaidiGanguShizuishanTongchuanZhongxiangDalad Yongzhou 9.5 Jiutai YingchengDehuiQixianXinleRunanDanzhouHuazhouLinzhouYunfu Anshun Yunyang LiuyangZhecheng Lucheng 0.5 LuannanDanjiangkouDangyangFangchengGuoyang NanxianShaodongEnping Guangzhou BaodeHuairenJieshouHuangpiShuangcheng CaijiapoMianxianGong'anLongxi Fula'erjiSuiHuarongYichuan ZhijiangYichengYiyangXiayiXinzhouAnpingMengchengXiangxiang ZhangjiajieXinyi BobaiHuangshanGaoan GaochengDanchengTaojiangLintao TaikangLushanJieshiLixin ZhangshuJinchangGaoyiHanchengJiaChangqiYanglingXiping KaramayGuipingFengshunXuanwei QingGuangshuiShangluoBiyangXilinhotWugangXincai Yushan Beijing XiaochangPuchengYuanpingHunchunFengqiuZijinNeihuangLianyuanFunan Kaili Qinyang YangchunPoyang WuyangJingbianSheLaoting ChizhouJinxianLaibinYugan Xingyi Huozhou LinglingTangyinYingshanJiazi QiyangYuWuchangHengyangDengfengYuanjiangShayuanXinfengYongjiShuchengXuwen County LongshanQuyangWangchengLujiangLongchuanXingguoLiulinQionghai Anning ShangliLingbi Luyi YingdeWenshan Bijie TianchangRongxianHengxianTongchengNingguo Lhasa DongxianDingyuanYuanyang Sanchahe XinxingJianshui Wuhe Huaiji Gejiu J nchéngji ng Mengzi Ningdu Pu'er Lijiang Duyun

Qianxi Kuqa Complementary cumulative Tongren 2.5 Liuzhi 0.0 Shanghai 9.0 5.5 6.0 6.5 7.0 7.5 5.5 6.0 6.5 7.0 7.5 5.5 6.0 6.5 7.0 7.5 City population, log s City population, log s City population, log s

37 Ferkéssédougou 0.0 Méagui Sinfra Abidjan Cote d Ivoire (CIV) Abengourou Soubré Cote d Ivoire (CIV) Cote d Ivoire (CIV) Duekoue Bouaflé = 1.27±0.31 1.0 Divo Man = 1.27±0.31 = 1.14±0.04

) 10.0 smin = 89903 Gagnoa San-Pédro smin = 89903 s

( Yamoussoukro y 0.5 r Korhogo F g g g 0.5 Daloa o o o l l

l Bouaké 9.5

, , , k

1.0 P n

n Bouaké o D i a t 0.0 r Abidjan G

Korhogo u

y 9.0 b y

t Daloa i t i r 1.5 i Yamoussoukro

t San-Pédro

C Gagnoa C s

i Man 0.5 Divo d

8.5 Duekoue 2.0 MéaguiAbengourou Soubré Complementary cumulative FerkéssédougouSinfra 1.0 Bouaflé 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Mamfe Yaounde Dschang Nkongsamba Ngaoundal 0.0 Bafia 10.0 Cameroon (CMR) EdéaBanki Cameroon (CMR) Cameroon (CMR) Mbalmayo Gamboru Douala 1.5 Batouri Ndop Meiganga Wum = 1.04±0.17 Tcholliré Foumbot = 1.04±0.17 Banyo Foumban = 1.05±0.13 Kumbo ) Limbé smin = 75908 TikoMbouda smin = 75908 9.5 s Mora Ebolowa

( Garoua Kribi Bertoua y r

F Kumba Kaele 0.5 Bamenda g g Yagoua Guider Bafoussam g Maroua 1.0 9.0 Ngaoundéré o o Ngaoundéré Bafoussam o l l N'Djamena

l Kumba

Wum Gamboru

, Bertoua , N'Djamena Kumbo

, Bamenda Yagoua

k FoumbanMbouda

P Limbé n Meiganga Tiko

n 8.5 Mokolo MamfeNgaoundal

o Ndop

D Foumbot

i Maroua

a BafiaDschang Kaele t NkongsambaMbalmayo Kribi 1.0 r Ebolowa Guider Mokolo G

u Batouri Mora

y 0.5 Garoua Edéa b y t

i 8.0 t i r Douala i Banki t Tcholliré C C s i

d 7.5 1.5 0.0 Yaounde Complementary cumulative 7.0 Banyo

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Dungu Kinzao ButusandeLubao Kamanyola Kinshasa Punia Kamituga 0.0 Kabondo-Dianda Kasaji Tshela Isiro Kenge Elila Matadi Kasumbalesa Iga Barriere Mulongo Demba Malemba Nkulu Mukubu Bandundu Kipushi Lisala Basankusu Sake Buta Aru Baraka Kambove Lubumbashi Mbuji-Mayi Kabinda Kindu ) Boma Bumba 1.5 Mwene-Ditu Kolwezi Kongolo Boma s ButondoAriwara Kilwa Rubaya ( 9 Manono Gemena Bukavu y r Mbandaka Butembo Beni F 0.5 Likasi Matadi Kananga Kamina Manono Bunia Kalemie Kisangani

Uvira g g Tshikapa g KikwitMongbwalu Kikwit Goma o o Kilwa Likasi o Mahagi l

l Tshikapa Kipushi l

Mbandaka Butembo 1.0 Kalemie Mukubu , , Bunia Goma KolweziKamina Mongbwalu , Ariwara Malemba NkuluKindu KasumbalesaMwene-DituButa k Kabinda Gemena

P Lisala n Kisangani Bukavu Bandundu Kabondo-Dianda Kongolo Mahagi n 8 Kambove o 1.0 Aru KamanyolaIsiroBaraka D Kasaji i

a Kananga LubaoDunguDemba t Iga Barriere Butondo r Beni Tshela Sake Bumba

G Basankusu

Mulongo Uvira u Rubaya

y Kenge b y t

i 0.5 t i Lubumbashi Kinzao r i

t Mbuji-Mayi C KamitugaElila C s Butusande

i 1.5

d Democratic Republicof the Congo (COD) Democratic Republicof the Congo (COD) 7

1.17 0.15 1.17 0.15 Democratic Republicof the Congo (COD) = ± 0.0 = ± Kinshasa Complementary cumulative smin = 123312 smin = 123312 Punia = 1.22±0.10 2.0 5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s 2.0 Orito Bogota Planeta Rica El Carmen de Bolívar ArjonaTurbaco PuertoPuerto Boyacá Asís ManuelaMorroa 0.0 ChiquinquiráChigorodó SantaSan Andrés Rosa de Cabal Calarcá Cereté Colombia (COL) FundaciónEspinal Sabanalarga Colombia (COL) Colombia (COL) San José del Guaviare Arauca San Vicente del Caguán SantanderLa Dorada de Quilichao MontelíbanoTiquisio Lorica Garzón MaganguéCaucasia AguachicaOcaña Sogamoso Tierralta Manaure Uribia Duitama Buga Ipiales 11.0 = 0.97±0.11 Pitalito Apartadó = 0.97±0.11 = 1.16±0.05 Ciénaga Fusagasugá ) Piedecuesta Cartago Tumaco Quibdó ZipaquiráFacatativá Girardot s s = 57330 Rionegro s = 57330 min Jamundí Florencia min Maicao Tunja

( 1.5 Yopal Tuluá Medellín Barrancabermeja Palmira Cali y 0.5 r Sahagún F Riohacha Sincelejo ArmeniaPopayán 10.5 El esfuerzo Bucaramanga g g Neiva g Valledupar Montería Barranquilla o o Villavicencio

o Pasto l l l

Manizales Santa Marta

, , , 1.0 Ibagué Pereira 10.0 Villavicencio Cartagena k P n 1.0 Cúcuta Cartagena n o

D Ibagué Pereira

i Valledupar Cúcuta

a El esfuerzo

t Bucaramanga Manizales r Barranquilla Barrancabermeja G

Neiva u

Santa Marta 9.5 Yopal Montería

y TunjaArmenia

b Tuluá y Riohacha t Palmira Pasto i 0.5 Cali t i Popayán AraucaBuga r i Zipaquirá 1.5 MaicaoCartagoFacatativá t Medellín PiedecuestaApartadóSahagún C Sincelejo San Andrés Rionegro C s GirardotJamundíDuitamaFusagasugá

i San JoséChigorodó del GuaviareAguachica Sogamoso 9.0 Orito MaganguéPitalito Ocaña

d ChiquinquiráManuela Florencia La Dorada Ipiales Calarcá Ciénaga SabanalargaUribia Quibdó El Carmen de Bolívar Garzón Tumaco MontelíbanoArjonaTiquisioCereté 0.0 Bogota EspinalTurbacoSanta Rosa de Cabal Fundación TierraltaCaucasia Puerto Asís 2.0 San Vicente delPlaneta Caguán Rica Lorica Santander de Quilichao Complementary cumulative 8.5 Puerto Boyacá Morroa Manaure

5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

Bayamo Havana 0.0 Cuba (CUB) Cuba (CUB) Cuba (CUB) Santa Clara 10.6 = 1.62±0.61 Camagüey = 1.62±0.61 = 0.99±0.02 ) smin = 203170 smin = 203170 s 0.5 Holguín ( y r F Santiago de Cuba 10.4 g 0.5 g g o o o l l l

, , ,

k 10.2 P n 0.0

n Havana o D i a t r G

u 1.0

y b y t 10.0 i t i Santiago de Cuba r i t C C s Holguín i 0.5 d 9.8 Camagüey 1.5 Santa Clara Complementary cumulative Bayamo

5.4 5.6 5.8 6.0 6.2 5.4 5.6 5.8 6.0 6.2 5.4 5.6 5.8 6.0 6.2 City population, log s City population, log s City population, log s

38 Prague 0.00 1.00 Hradec Králové 11.00 Czech Republic (CZE) Pardubice Haví ov Czech Republic (CZE) Czech Republic (CZE)

= 1.42±0.43 Liberec = 1.42±0.43 Olomouc = 1.27±0.06 ) 0.25 smin = 156917 smin = 156917 s 0.75 10.75 Ústí nad Labem (

Pilsen y r F g g g 0.50 0.50 o o Ostrava 10.50 o l l l

, ,

, Brno k P n 0.75 n 0.25 o Brno D 10.25 i a

t Ostrava r G

y b y 1.00 t i t i 0.00 Prague r i 10.00 t C C s Pilsen i Ústí nad Labem d 1.25 0.25 HavíOlomoucov 9.75 Liberec Pardubice Complementary cumulative 1.50 0.50 Hradec Králové 9.50 5.2 5.4 5.6 5.8 6.0 6.2 5.2 5.4 5.6 5.8 6.0 6.2 5.2 5.4 5.6 5.8 6.0 6.2 City population, log s City population, log s City population, log s 11.50 Oldenburg Erfurt Rostock Dortmund 0.0 Aschaffenburg Gießen Germany (DEU) Würzburg Ulm Germany (DEU) Germany (DEU) Regensburg Lübeck Halle (Saale) HeilbronnMagdeburg Reutlingen Darmstadt 1.5 Osnabrück Kiel = 1.02±0.15 Kassel Freiburg im Breisgau = 1.02±0.15 = 1.03±0.02

) Chemnitz 11.25 Hamburg Münster Koblenz Brunswick CologneFrankfurtBerlin am Main

s s = 238770 s = 238770 min HeerlenAugsburg min

( Mönchengladbach Saarbruecken Munich y r Leipzig Bielefeld Stuttgart F Karlsruhe Wuppertal 11.00 g g Dresden g 0.5 Bremen o o Dusseldorf Hanover Mannheim o l l

l Nuremberg

1.0 Nuremberg , , Mannheim ,

k 10.75 P n Stuttgart Munich BremenDusseldorf n

o Hanover D i Karlsruhe a t Hamburg Wuppertal r

G Bielefeld Frankfurt am Main Dresden u 1.0 10.50

y AugsburgLeipzig b y

t Heerlen i SaarbrueckenMönchengladbach t i Kassel r 0.5 i Cologne Freiburg im Breisgau t Münster

C DarmstadtKoblenz Brunswick C s Reutlingen

i Berlin 10.25 Heilbronn Chemnitz Osnabrück Regensburg

d Würzburg Gießen Ulm Kiel

Aschaffenburg 1.5 Halle (Saale) LübeckMagdeburg 10.00 Oldenburg Complementary cumulative Erfurt Rostock 0.0 Dortmund 5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s Santo Domingo 0.0 Dominican Republic (DOM) San Juan de la Maguana Barahona Dominican Republic (DOM) Dominican Republic (DOM) 1.0 Baní Bonao = 1.05±0.28 Moca Puerto Plata = 1.05±0.28 = 1.03±0.02 ) smin = 76888 La Vega San Francisco de Macoris smin = 76888 10.5 s

( San Pedro de Macorís y r F 0.5 Higüey g g g o o 0.5

o La Romana l l Santiago de los Caballeros l

, , Santiago de los Caballeros 10.0 , k P n n o D i a

t 1.0 r G

y 0.0 Santo Domingo La Romana b y t i t i San Pedro de Macorís Higüey r i 9.5

t La Vega C C s

i San Francisco de Macoris Moca Puerto Plata

d 1.5

Baní Bonao

0.5 San Juan de la Maguana

Complementary cumulative 9.0 Barahona

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Mascara Bou Saada Tizi Ouzou Algiers Khenchela Touggourt 0.0 Ain Turk Relizane Algeria (DZA) Aflou Laghouat Algeria (DZA) Algeria (DZA) 1.5 Medea Saida Guelma Bashar M'Sila El Eulma Souk Ahras Tiaret Jijel Bordj Bou Arreridjj = 1.63±0.26 Bejaia = 1.63±0.26 = 1.03±0.02 Tebessa Ouargla ) Biskra El Oued smin = 143042 Sidi Bel Abbes Skikda smin = 143042 s Mostaganem Tlemcen ( 0.5 10.5 1.0 Ash Shalif y r Batna F Setif Djelfa g g g Annaba o o Blida o

l Oran l l

1.0 , ,

, 0.5 Constantine Constantine k P

n Oran n

o 10.0 D Blida i a

t Annaba r G

u 1.5

y 0.0

b Algiers y

t Djelfa i Setif t i Ash Shalif r i Batna Tlemcen t Mostaganem C C s

i Sidi Bel Abbes Skikda BiskraOuarglaEl Oued 2.0 Guelma Tebessa d 9.5 BejaiaEl Eulma Jijel TiaretBordj Bou Arreridjj

0.5 Souk Ahras AflouM'Sila SaidaBashar KhenchelaMedea Ain TurkTouggourtLaghouat MascaraRelizane

Complementary cumulative Tizi Ouzou Bou Saada 2.5 1.0 5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s

0.0 Ecuador (ECU) Milagro Esmeraldas Ecuador (ECU) Ecuador (ECU) Guayaquil Quevedo City of Loja 10.50 Quito = 1.42±0.35 1.0 Ibarra Riobamba = 1.42±0.35 = 1.14±0.07 Santa Elena Machala ) smin = 162496 smin = 162496 s Portoviejo Cantón Santo Domingo de los Tsachilas ( 10.25 y r F Manta 0.5 Ambato g g g o o o l l

l 10.00

0.5 Cuenca

, , , k

Quito P n Cuenca n o D i a

t 9.75

1.0 r Cantón Santo Domingo de los Tsachilas G

y b y t i t i 0.0 Ambato r Guayaquil i t 9.50 C Manta C s Santa ElenaPortoviejoMachala i d 1.5 Milagro 9.25 Esmeraldas QuevedoIbarra Complementary cumulative 0.5 Riobamba City of Loja

5.25 5.50 5.75 6.00 6.25 6.50 5.25 5.50 5.75 6.00 6.25 6.50 5.25 5.50 5.75 6.00 6.25 6.50 City population, log s City population, log s City population, log s

39 MutubisAtmeeda Cairo Beni Adi NazlatBilqas Dahrut Al BalakusQotur Dayr Abu Maqrufah 0.0 Al Inab al Saghirah MinshatAl Zayyatin al Maghaliqah al Qibliyah Egypt (EGY) GeziretHurghada Nikla Ezbet al Sibil Qiblli Egypt (EGY) Egypt (EGY) Hawsh Isa Sadat City Ihnasya al Madinah Delingat Qaryat Salah al Din Kafr al Hajj Shirbini Shubra Khit Ezbet Ahmad al Ballasi Shamma Kafr Dinshawai Zayyan Al Arish Nazlet Abu Haseiba 9th Idku = 0.66±0.08 Fashn Faris = 0.66±0.08 11 = 1.04±0.03 Kafr al Hawi Isna Abu Huminus Abu Manna Gharb Alexandria

) Itay el-Barud Marsa Matruh Maghagheh 6th of October Mansoura s = 55976 Abu Tig Nazlet al Mughrebi s = 55976 min Fazarah bil Qaryah min s Beni Mazar Qena Tala Rosetta ( 1.5 Kafr Al-Sheikh Shimm al Basal al Bahari Samalut y

r Aswan Girga Mallawi Luxor Asyut F Damanhur An Nazlah Faiyum 0.5 SuezDesouk Al-Minya Sohag Al-Manzala Tanta g g Damietta El Mahalla El Kubra g Al Ismailiya Al-Zaqaziq Port Said Port SaidBanha Beni Suef o o Banha Beni Suef 10 o Damanhur Al Ismailiya l l Suez l Al-Zaqaziq DamiettaAl-Manzala

El Mahalla El Kubra MallawiDesouk

, ,

, Al-MinyaTanta AswanKafr Al-SheikhSamalut k Beni Mazar P n 1.0 An Nazlah Faiyum Shimm al Basal al Bahari Abu Huminus Tala n FashnMaghaghehAbu TigQenaNazletRosetta al Mughrebi o 6th of October D ItayFazarah el-Barud bil Qaryah i 1.0 Girga BilqasShamma a Sohag

t Al BalakusNazlet Abu Haseiba Ezbet al Sibil QiblliIdku Marsa Matruh r Kafr9th al Hawi Isna Nazlat Dahrut Abu Manna Gharb

G Shubra Khit Ihnasya al MadinahZayyanDelingat Faris u 9 Qotur Ezbet Ahmad al Ballasi Al Zayyatin al Qibliyah Minshat al Maghaliqah MutubisAtmeeda SadatKafr alCity Hajj Shirbini Luxor Geziret Nikla y Al Inab al SaghirahHawsh Isa KafrAl Arish Dinshawai

b Asyut y t i t i

r Dayr Abu Maqrufah i Qaryat Salah al Din t

C 0.5

Alexandria C s Beni Adi i 1.5 Mansoura d 8 Complementary cumulative Hurghada 0.0 Cairo 5.0 5.5 6.0 6.5 7.0 7.5 5.0 5.5 6.0 6.5 7.0 7.5 5.0 5.5 6.0 6.5 7.0 7.5 City population, log s City population, log s City population, log s

Torrevieja Madrid 1.5 Vitoria-Gasteiz Almeria 0.0 Jerez Cádiz Spain (ESP) Elx / Elche Gijón Spain (ESP) Spain (ESP) CórdobaSantander Oviedo Pamplona = 1.26±0.23 Castellón de la Plana Valladolid = 1.26±0.23 = 1.11±0.05 Barcelona A Coruña Donostia / San Sebastián ) smin = 229333 Tarragona Santa Cruz de Tenerife smin = 229333 s Vigo Granada ( Alacant / Alicante Palma de Mallorca y r

F 11.0 0.5 1.0 Las Palmas de Gran Canaria Murcia g g g Zaragoza o o Bilbao Valencia o l l l

, ,

, Málaga

k Seville P n Bilbao n o

D Seville i a 0.5 t 1.0 Valencia 10.5 Zaragoza r Málaga G

u Barcelona y b

y Palma de Mallorca t

i Murcia t i Las Palmas de Gran Canaria r i Donostia / San Sebastián

t Tarragona Pamplona C Alacant / AlicanteVigo C s Valladolid i Castellón de laA Plana Coruña Vitoria-Gasteiz SantaGranada Cruz de Tenerife d 1.5 0.0 Madrid 10.0 OviedoSantander Elx / Elche Gijón Jerez CórdobaCádiz Complementary cumulative Almeria Torrevieja

5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s

GojoAykelRobe Town Debre Tabor KaramileKebri Dehar DanglaFik' BekojiGelemso ArbeDambidollo Gona Digna 0 AdigratShone BilaMeki ChagniKersa AmboNefas Meewcha Ethiopia (ETH) DodollaBoditi ChuckoJimma Ethiopia (ETH) 10 Ethiopia (ETH) Ebbenat Nekemte DebarkAsella Gut'OAsassa Addis Ababa Debre BirhanAdet SekaHarar GundoDebre Meskel Markos Areka Kobo MotaJinka = 2.75±0.34 MorochoKefole Chiro = 2.75±0.34 Aletawendo Kombolcha = 1.59±0.23 Wereta / Bule Hora ) Bonsa Bota Gonder smin = 169827 Merawi Alaba Kulito smin = 169827

s Yirgalem Arsi Negelle Quiha ( 1 Bishoftu Jijiga y r DessieHawassa F 1 Shashamane Dire Dawa 9 Adama g g

g Shashamane Dire Dawa o o Adama Bishoftu HawassaBahir Dar

o Bahir Dar Bonsa Bota l l Kefole / BuleJijiga Hora l Yirgalem

DessieQuiha Asella Seka GonderAlaba Kulito , , ArbeDangla Gona Jimma , DodollaMeki Arsi Negelle Debre BirhanNekemte

k Debre Markos KersaRobeAmbo P n 0 Gojo Town ArekaChiroKombolcha Addis Ababa Dambidollo Digna

n Asassa Debre TaborGundoAykel MeskelMorochoChucko Harar

o Adigrat D

i 2 Fik'Shone Merawi a Bila

t 8 Wereta Boditi r G

Karamile

u GelemsoAletawendo

Chagni

y Ebbenat b

y Kobo t Bekoji i t i Gut'O r

i Mota Debark t

C Kebri Dehar C

s 1 i 3 d 7

Nefas Meewcha Complementary cumulative 2 Adet Jinka 5.25 5.50 5.75 6.00 6.25 6.50 5.25 5.50 5.75 6.00 6.25 6.50 5.25 5.50 5.75 6.00 6.25 6.50 City population, log s City population, log s City population, log s

Reims Paris Nimes Valence Brest 0.0 Le Mans Le Havre France (FRA) 1.5 Dijon Angers France (FRA) France (FRA) Metz Valenciennes Orléans Perpignan Caen Nancy Clermont-Ferrand = 1.40±0.24 Avignon Tours = 1.40±0.24 = 1.10±0.02 Saint-ÉtienneMulhouse ) Rennes Grenoble smin = 285722 Montpellier smin = 285722 s Toulon 11.5 ( 0.5 Rouen Strasbourg y r 1.0 Lens F Nantes g g Bordeaux Lille g o o

o Nice l l

l Toulouse

, , , 1.0 11.0 k

0.5 P n Marseille

n Lyon

o Lyon D i a t Marseille r G

Nice

y Toulouse

b Bordeaux y

t Lille i t 1.5 i Nantes r 0.0 Paris i 10.5

t Lens C Rouen Strasbourg C Toulon s Grenoble i Rennes Saint-Étienne Montpellier d

Avignon Mulhouse Tours Nancy 2.0 OrléansCaenAngersClermont-Ferrand ValenceLe MansDijonMetzLeValenciennes Havre 0.5 Perpignan Brest

Complementary cumulative 10.0 Reims Nimes

5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

York London Lincoln Bedford Exeter Blackpool Chelmsford 0.0 Larbert Northampton Dundee Corby United Kingdom (GBR) Aberdeen Milton Keynes United Kingdom (GBR) United Kingdom (GBR) Chester Ipswich Colchester Preston Luton Swindon Plymouth Reading Peterborough Cambridge = 1.19±0.17 1.5 Norwich Oxford = 1.19±0.17 = 1.00±0.02 Blackburn Hull ) Gloucester Swansea Stoke-on-Trent smin = 264763 Coventry smin = 264763 11.5 s Blackwater Southend-on-Sea Middlesbrough ( Bournemouth Edinburgh y r Brighton F 0.5 Leicester Belfast Cardiff g g

g Sheffield Bristol o o 1.0 o Portsmouth Nottingham l l l

Birmingham

, , Liverpool Manchester , Newcastle upon Tyne 11.0 Leeds k P n PortsmouthGlasgow 1.0 n Glasgow

o Leeds D

i Liverpool a t

r Nottingham

G Newcastle upon Tyne

0.5 Birmingham Brighton Bristol

y Sheffield b y t i Blackwater t i Manchester Belfast Cardiff r i 10.5 Bournemouth t Oxford C LeicesterEdinburgh C s 1.5 Southend-on-Sea

i Reading Gloucester Stoke-on-TrentCoventryMiddlesbrough d MiltonNorwich Keynes SwanseaHull PeterboroughPlymouthBlackburnCambridge 0.0 Luton Swindon London Colchester Preston DundeeAberdeen Ipswich LarbertChester Bedford Northampton ChelmsfordCorby

Complementary cumulative Blackpool 10.0 LincolnExeter 2.0 York 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

40 Dunkwa-on-OfinNalerigu Wenchi Karaga Accra Yeji Konongo 0.0 Lomé Nkoranza Dambai Salaga Ghana (GHA) Assin Fosu Akatsi Ghana (GHA) Ghana (GHA) UNNAMED Nkawkaw Atebubu New Tafo Akim Suhum Asamankese 1.5 Dormaa Ahenkro Wa Kumasi Winneba Navrongo 10 = 1.25±0.18 Walewale Kintampo = 1.25±0.18 = 1.14±0.05 Bimbila Gushiegu Akim Oda ) Garu Savelugu smin = 56239 Berekum Bolgatanga smin = 56239 s Somanya Agona Swedru

( Takoradi [Sekondi-Takoradi] 0.5 Sunyani Techiman y r Obuasi

F Hohoe 1.0 Bawku Ho Tamale g g

g Cape Coast KoforiduaYendi Koforidua o o 9 o l l Ho

l Hohoe

Cape Coast Sunyani Obuasi Tamale BolgatangaTechimanYendi , 1.0 , Akim OdaSomanyaWaGaru BawkuAgona Swedru , WinnebaBerekumSavelugu k WalewaleNavrongo

P Kintampo n DormaaSuhum AhenkroAkatsi Asamankese 0.5 Takoradi [Sekondi-Takoradi] Assin Fosu NkawkawNew Tafo Akim n Atebubu Gushiegu Nkoranza o Dunkwa-on-Ofin D

i Kumasi Yeji Salaga a

t WenchiLomé r Konongo G

u 8 Dambai 1.5 Karaga y b y t i t i r 0.0 Accra i t UNNAMED C C s Bimbila i d 2.0 7 0.5 Complementary cumulative

Nalerigu 2.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Banankoro Conakry 0.0 Guinea (GIN) BokéGuéckédou Guinea (GIN) 9.5 Guinea (GIN) Macenta Sinko = 1.24±0.30 Sangarédi Mamou = 1.24±0.30 1.0 Faranah Kissidougou = 1.12±0.07 ) smin = 61481 smin = 61481 s Kamsar Labé ( y r

F Siguiri 0.5 Kindia 9.0 g g g o o

o 0.5 Nzérékoré l l l

, , Kankan , k P n

1.0 n Nzérékoré

o 8.5 D i a

t Kankan r

G Siguiri

u 0.0 Conakry Kindia y b y

t Labé i t i r i Kissidougou t Faranah C C s 1.5 Macenta i 8.0 Sangarédi

d Boké Mamou 0.5 Kamsar Sinko Guéckédou Complementary cumulative 2.0 7.5 Banankoro 5.0 5.5 6.0 5.0 5.5 6.0 5.0 5.5 6.0 City population, log s City population, log s City population, log s 1.00 Athens 0.00 Chalkida Chania Greece (GRC) Greece (GRC) 11.0 Greece (GRC) Ioannina = 0.91±0.29 Volos = 0.91±0.29 = 1.09±0.04

) 0.75 smin = 76841 Larissa smin = 76841 s 0.25

( Heraklion y r F g g g 0.50 0.50 Patras 10.5 o o o l l l

Thessaloniki , ,

, Thessaloniki k P n 0.75 0.25 n o D i a t r G

10.0 y b

0.00 y 1.00 t Athens i t i r i t C C s i PatrasHeraklion d 1.25 0.25

9.5 Chalkida Larissa Volos Chania Complementary cumulative 1.50 0.50 Ioannina

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Fray Bartolomé de las Casas Guatemala City El Estor Poptún 0.0 San Cristobal Verapaz Joyabaj 1.5 Cubulco Santa María Nebaj Guatemala (GTM) Aguacatán Playa Grande Ixcán Guatemala (GTM) Guatemala (GTM) Morales Puerto BarriosChiquimula Chichicastenango Panajachel [Sololá] = 1.26±0.22 Ciudad Vieja [Antigua Guatemala] Retalhuleu = 1.26±0.22 = 1.40±0.20 MalacatánSanta Lucía Cotzumalguapa

) Jutiapa Coatepeque smin = 56714 Santa Cruz del Quiché Flores smin = 56714 s Mazatenango ( 0.5 Santa Cruz Barillas 10 y r 1.0 Jalapa Huehuetenango F Totonicapán Escuintla g g g o o

o San Pedro Sacatepéquez [San Marcos] l l Chimaltenango

l Quetzaltenango

, , , 1.0 k 0.5 Quetzaltenango P n Mazatenango Chimaltenango n Cobán RetalhuleuSanta LucíaEscuintla Cotzumalguapa o Ciudad Vieja [Antigua Guatemala]

D Cobán i

a 9 Puerto Barrios t Morales Coatepeque r G

u San Pedro SacatepéquezEl Estor [San Marcos] MalacatánChiquimula Flores y Totonicapán b y t

i Santa Cruz del Quiché t 1.5 i 0.0 r Guatemala City i Panajachel [Sololá] Jalapa t Jutiapa

C Chichicastenango Huehuetenango C s Poptún

i Santa María Nebaj Aguacatán San Cristobal Verapaz

d Joyabaj 8 Santa Cruz Barillas 2.0 0.5 Cubulco Playa Grande Ixcán Complementary cumulative Fray Bartolomé de las Casas

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Santa Rosa Tegucigalpa 0.00 Honduras (HND) Danli Honduras (HND) Honduras (HND) 1.00 Siguatepeque San Pedro Sula = 1.12±0.31 Barrio El Centro Puerto Cortés = 1.12±0.31 = 1.19±0.05 ) smin = 60541 Comayagua smin = 60541 s 0.25 El Progreso ( 9.5 y r 0.75 F Villanueva g g Choluteca g o 0.50 o o l l l

, , 0.50

, La Ceiba k P n

n La Ceiba o San Pedro Sula 0.75 D i

a 9.0 t r G

Villanueva

u 0.25

y b y t i t i El Progreso r 1.00 i Comayagua t Puerto Cortés C C s

i Choluteca 0.00 Tegucigalpa Barrio El CentroSiguatepeque d

8.5 1.25 Danli

Complementary cumulative 0.25 Santa Rosa 1.50 5.0 5.5 6.0 5.0 5.5 6.0 5.0 5.5 6.0 City population, log s City population, log s City population, log s

41 Petite-Rivière-de-l'Artibonite 0.0 Mirebalais Jérémie Port-au-Prince Haiti (HTI) La Croix Dessalines Haiti (HTI) Haiti (HTI) Saint Michel de l Attalaye Gros Morne Hinche Ouanaminthe = 1.37±0.30 Les Cayes Anse-à-Galets = 1.37±0.30 9.5 = 1.20±0.12

) 1.0 smin = 57858 Saint-Louis-du-NordJacmel smin = 57858 s Port-de-Paix Léogâne ( 0.5 y r F Petit-GoâveSaint-Marc 9.0 g g Ouanaminthe g o o

o 0.5 Cap-Haitien Cap-Haitien l l l

Gonaïves , , Gonaïves , 1.0

k Port-de-Paix P

n 8.5 n o

D Léogâne

i Saint-Louis-du-NordPetit-Goâve

a Saint-Marc

t Les Cayes r 0.0 Port-au-Prince Gros Morne Jacmel G

Anse-à-Galets

y Hinche b y 8.0 1.5 t Mirebalais i t i Jérémie r i t C

C Petite-Rivière-de-l'Artibonite s

i Dessalines 0.5 La Croix d

7.5 2.0

Complementary cumulative 1.0 7.0Saint Michel de l Attalaye 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s Budapest 0.0 Hungary (HUN) Nyíregyháza Gy r Hungary (HUN) 11.0 Hungary (HUN) Pécs = 1.56±0.55 Szeged = 1.56±0.55 = 1.30±0.10 ) smin = 152840 smin = 152840 s 0.5 Miskolc ( y r

F 0.5 Debrecen g g g 10.5 o o o l l l

, , , 0.0 Budapest k P n

1.0 n o D i a t r G

u 10.0 y b y t i t i r 0.5 i t

C Debrecen C s 1.5 i Gy r

d Miskolc Pécs 9.5 Szeged

Complementary cumulative 1.0 Nyíregyháza 2.0 5.25 5.50 5.75 6.00 6.25 6.50 5.25 5.50 5.75 6.00 6.25 6.50 5.25 5.50 5.75 6.00 6.25 6.50 City population, log s City population, log s City population, log s

Sarolangun Jakarta ArgaMannaNunukanTeluk Makmur Dalam PasangkayuBajawa MeulabohBaturajaGunungParigi Tua PangkalanTomohonTondano BrandanBuol PangkalanSidorejoPagar Alam Kerinci PangkajeneKotarajaMajene BontangUteun Geulinggang KualaSerui TungkalUnit Kota 1 MarabahanMuaraTebing Enim Suluh MuarateweUjungBesuki Batu SidikalangPerawangBatu Licin 0.0 SelatpanjangSidogedeBarabai SolokBantaengLepar Samura Tobelo Ampana TanjungPagojenganTarutung Pandan PinrangTolitoliNgabang Indonesia (IDN) LhokseumaweMasohi NgawiSentaniTanjung Balai Karimun Indonesia (IDN) Indonesia (IDN) Palangka Raya KandanganRembangBangko TanjungBaganSemurup Siapi-api Redeb AtambuaKalabahiKalianda MekarjayaGunungsitoliSigli TanjungNangaToboali Pinoh Enim PosoLamonganBalangnipa AsembagusBiakDompu Kota PayakumbuhPalopoPrabumulih CiherangKotabaruBima BaubauPanyabunganBireuen SengkangNegaraNabire ManokwariPurworejo KolakaLahat LabuanPare-PareKogimage Bajo KotamobaguTembilahanAmuntai PangkajeneBanjar KotabumiLangsa RahaCurupSumuradem MaumereBabatEnde LubuklinggauJayapura = 0.84±0.06 RantepaoWates MekarjayaWatampone = 0.84±0.06 Kuala Kapuas Merauke = 1.05±0.04 WonosariTernate TakatidungSumbawa Besar BulukumbaLuwuk KisaranSitubondo PemangkatPandan Singkawang Rantau PrapatBlora TambolakaMuara Bungo ) Bitung DuriMamuju WaingapuPringsewu MuncarTuban Bandung TarakanSorong SangattaSintang 11 Surabaya 2.0 Genteng Tanjung Balai s = 78887 Banyuwangi s = 78887 min Tebing TinggiDumai PalabuhanratuPacitan min s Gunung RutengSugih PadangSampit Sidempuan Samarinda SragenTimika BondowosoPalu Medan PangkalpinangKetapang SampangKel Rijali

( Pandeglang Purworejo 0.5 BukittinggiBengkulu PurwodadiRangkasbitung Yogyakarta Bojonegoro Banda Aceh Jombang PematangBangkalan Siantar Kendari JakartaSemarang y r Trenggalek SingarajaKupang Subang Balikpapan Batam City F Manado Indramayu Pekanbaru Metro SumenepPamekasan PalembangTasikmalaya Tegal GorontaloKebumen Lumajang Madiun CirebonJakarta Balikpapan Probolinggo Jember Cilacap Kuningan MakassarDenpasarPurwokerto g g Cilegon Selong Garut Pati g Tulungagung Jambi Mojokerto Pontianak Batam City Temanggung Magelang Banjarmasin Tarakan Tulungagung Pekalongan o o Padang Pekanbaru TemanggungPadang o Mojokerto Mataram Jambi

1.5 l l 10 Dumai Bandar Lampung Samarinda Bandar Lampung Duri JakartaJakarta l Banjarmasin Pontianak Magelang

Jakarta CilegonPamekasanLumajangProbolinggo

Pekalongan Bangkalan Jakarta Bontang Metro ,

, Trenggalek Madiun BojonegoroSampangGorontaloIndramayu Mataram , 1.0 Pangkalpinang Jember Perawang Kendari Kuningan Palembang Selatpanjang Tembilahan SingarajaManado Tanjung Balai Karimun TubanPematangPadang Siantar Sidempuan k Pati BondowosoBandaKebumenCilacap Aceh Garut Tanjung Enim Kotabumi TebingWatampone TinggiSangatta Purwodadi P PacitanBukittinggi n Denpasar SidogedeSidikalangBagan Siapi-api MuncarTimika UjungBaturaja Batu Tanjung RedebBabatBulukumbaJayapuraSitubondo Purworejo Tasikmalaya Ciherang WonosariPurworejo Bengkulu Balangnipa GentengMuara Bungo Palu SubangSumenep n Tarutung Sorong Tanjung Balai Cirebon Nunukan PrabumulihPandanLubuklinggauMamujuKetapang Selong Makassar Pangkalan Kerinci ManokwariPringsewuBanyuwangiSragen Sampit o Biak Kota Palabuhanratu Ngawi LangsaPanyabungan D Unit 1 SengkangGunungsitoliRantau PrapatBitungRaha

i PangkajeneSumurademPandeglangKupang Purwokerto Tanjung PandanAsembagusBangko NegaraAmuntaiLamonganLahat BanjarBloraMerauke Rangkasbitung a Jakarta Rembang Kotamobagu Kolaka t Palangka Raya SemurupPosoPayakumbuhGunung Sugih Pagojengan Pinrang PalopoBaubauPare-Pare Sintang r Pangkalan Brandan TebingLepar Suluh BimaSamura CurupKisaranSingkawang 1.0 Jakarta Pasangkayu Gunung TuaBantaeng KualaBireuenWates Kapuas Sumbawa Besar Tegal BarabaiSidorejoSolok G 9 Lhokseumawe Sentani PangkajeneSeruiKotaraja KotaPagar AlamDompuSigli

u Luwuk Besuki Waingapu Kel Rijali BuolTelukKandangan DalamMarabahan Takatidung MaumereEnde ParigiToboaliBatu Licin y Semarang Jombang Uteun GeulinggangTolitoli Kogimage

b 1.5 Tondano Pemangkat y Ampana t Tomohon Ternate MannaMeulaboh Ruteng i Muara EnimMekarjaya t i Nanga PinohKotabaru Kuala TungkalArga Makmur Tambolaka r

i Rantepao Medan Kalianda Mekarjaya

t Yogyakarta MajeneSarolangun Kalabahi

C Tobelo Ngabang C s Nabire i 0.5 Bandung Atambua 8 Bajawa d 2.0 Surabaya

Muaratewe Labuan Bajo Complementary cumulative Masohi 0.0 Jakarta 7 5.0 5.5 6.0 6.5 7.0 7.5 5.0 5.5 6.0 6.5 7.0 7.5 5.0 5.5 6.0 6.5 7.0 7.5 City population, log s City population, log s City population, log s

RasraRanaghat Delhi [New Delhi] GudiyathamSujangarhNajibabadSawai Madhopur Bhadrak HabraVidishaUdupiBijnorMadhubani NanjanaguduKishanganjDungarpurHazaribaghFaridpurSalon PuriBhusawalPollachiMahuvaKalpiBandikui DoddaballapuraAmbikapurKrishnagiriBhongaonWardhaAuraiya NawabganjJhumriDehriNawadaKarauliNeyveli Tilaiya PorbandarGauriganjSingrauliMotihariRaichurManjha KasganjLaharpurBilariContaiAizawl BhadravathiChittorgarhEtmadpurAlirajpurBoisar AliganjRudauliSataraRajsamandJoura NotunKannaujSivani AmethiDighaPali MasaurhiChhindwaraSiwanMusiriBasirhat MandiMuhammadabad GobindgarhAnandSangamnerMorbi Nagar RafiganjBuxarKishangarhBurhanpurHassan 0.0 NarkatiaganjShivpuriKundaBatalaGuna ValsadVizianagaramViluppuramPanaji Dhaulpur GangapurJunagadhYavatmalBeawarPatan DurgapurOraiBansiKarimnagarChitradurga Mumbai SurendranagarNautanwaDatia KadapaBankuraBhagalpur India (IND) KurukshetraChhatarpurBisalpurMirganjKarur SasaramMandyaShahganjDewas India (IND) India (IND) MedininagarTilharSirsa LalganjTinsukiaTezpurDausa CuddalorePalanpurJolarpetBaberuSidhi KhammamMaholiMuhammadabadTirunelveli ChhibramauMohmadiBalotraRatlamKatni NavsariDobhiBarmerKolar Gold Fields KaimganjRudrapurDamohBaheri DhubriEluruTandaBathindaPhalodi GopiganjJalalpurBharuchBidar KhandwaMaharajganjHindaunChandrapur KanchipuramNawabIchalkaranjiKarnal Ganj TiruvannamalaiUtraulaSidhauliThanjavur 2.5 RajapalaiyamBharatpurGodhraBeed BhilwaraKorbaKhurjaRewa, Madhya Pradesh Parbhani UNNAMEDNizamabadBhiwaniJalna NagapattinamColonelganjSiuriAnantapur KumbakonamMahesanaPuruliyaRoorkee BanswaraThoothukudiHajipurAkbarpur ShikohabadDharmapuriChandausiLalitpur ShivamoggaKhagaDahodObra ChitrakootNaugarhMorenaNadiad Dham UNNAMEDHathrasShahabadThakurdwara CoochHaldwaniDeoghar BeharArrah KurnoolBhindKakinadaKhalilabad Chennai DhupguriShillongSaharsaPurnia BhingaSikarTumakuru Kolkata Palakkad GandhidhamSamastipur = 1.27±0.06 Mansa NagarSatnaLaturHisar BettiahRajamahendravaramUjjain = 1.27±0.06 SansarpurDimapurDeesa AmrohaBrahmapurVijayapura 1.14 0.03 Ahmedabad Jangipur = ± WaidhanMunger Dindigul 11 DarjeelingKatiharRohtak SandilaGhazipur ChhapraTulsipurEtah BulandshahrLauiyahDavanagere Nandangarh Bengaluru GiridihUnnaoAjmer BhiwadiDhulianRampur Pune AhmadnagarTirupatiBallari SilcharNanparaBihar Sharif Surat ) KharupetiaMainpuri JaynagarBalrampur Majilpur Hyderabad SitamarhiBandaPilibhit PiprakothiPanipatAkola AzamgarhBhadohiDhule BalliaCanningRaiganj BhavnagarEtawahJhansi JalgaonMauSiliguri Kochi s = 294424 NeeleshwaramYamunanagarPatiala RourkelaNanded-WaghalaAmravati s = 294424 Deoria Ambala min JamnagarKharagpurBudaun AnandKalaburagiNellore min s SambhalSagar Nagercoil ThiruvananthapuramKozhikode JalpaiguriBasti DarbhangaAgartalaMalegaon ChandigarhNagpur RaeAnantnagAlwar Bareli WarangalUdaipur Guruvayur 0.5 BikanerSitapur KushinagarFirozabadJaunpur Ludhiana PratapgarhDharwad Bokaro Nashik Kollam Varanasi ( Cuttack Aurangabad AlappuzhaThrissurBelagavi BahraichKrishnanagar VadodaraAsansolCoimbatore MathuraGuntur MirzapurFatehpur Thalassery Salem Jaipur Nibua RaiganjDhing PuducherrySangli Rajkot Lucknow MuzaffarnagarLakhimpur GayaKolhapur Kannur

Bilaspur y r Gonda Solapur VijayawadaAmritsarVasai-Virar Madurai Kanpur 2.0 ShahjahanpurFarrukhabad SilvassaKota Hisar VelloreMangaluru VisakhapatnamPrayagraj Indore Haridwar BerhamporeTiruppur TamlukAgra F Imphal JammuDehradun KolhapurJalandharEttumanoorTiruchirappalliErodeBhopalMeerut BegusaraiSolapur Jalandhar Patiala AlappuzhaAmbala Sangli KolkataMysuru BhubaneshwarMalappuram DurgSaharanpur Rohtak Thrissur Silvassa Gorakhpur Patna Hardoi Ettumanoor PanajiKarnal Yamunanagar BelagaviGuntur Puducherry Faizabad Sultanpur AmritsarTiruchirappalli 10 Sirsa PanipatNanded-WaghalaMalappuramSrinagar GwaliorAligarh Srinagar Ahmadnagar Malegaon Jodhpur Vellore Mysuru Kurukshetra Bhiwani Dhule AkolaAmravatiAnand Warangal Jabalpur Guwahati Vasai-Virar CoochBathinda Behar BhiwadiJalgaonBhubaneshwar Moradabad g

g Rajkot Tiruppur Anantnag Guwahati Thalassery Jodhpur Neyveli PalakkadDeesaLatur TirupatiJamnagarMuzaffarnagarKalaburagiNelloreNagercoilKrishnanagarKotaDurg RaipurBareilly Kolkata Mandi GobindgarhBatala JalnaNadiad NeeleshwaramBhavnagarDindigulKharagpurJalpaiguriNibuaBikaner RaiganjImphalSultanpurJammu Ranchi Dhanbad g Vijayawada Erode BallariDeoriaAlwarDharwadAgartalaCuttack KuchaiSaharanpur Kot Mangaluru Raipur Yavatmal ParbhaniBeedDharmapuriMahesanaNagapattinamRajamahendravaramMathuraBilaspurHaridwarHardoi Muzaffarpur Berhampore Kannur Sangamner Rajapalaiyam KurnoolBhadohi RaiganjSiliguriShahjahanpurLakhimpurDhingUdaipurJaunpur Jabalpur JunagadhJolarpetIchalkaranjiBharuchNizamabad KakinadaVijayapuraSagarBastiRourkelaFarrukhabadFatehpurGondaDehradun Jamshedpur Gwalior Bareilly PatanPalanpur Kumbakonam DimapurAnantapurThoothukudiShivamoggaDahodGandhidhamEtahAjmer DavanagereSambhalRaeCanning PratapgarhBareliSitapur FirozabadBahraich Kuchai Kot Jamshedpur WardhaSurendranagarBhusawalSatara GodhraTirunelveliNavsari DhupguriShillongSiuri Tumakuru Jaynagar MajilpurBegusaraiMirzapur o Chandrapur Sikar o Bidar Eluru Aligarh KrishnagiriBoisarMorbiPollachiContaiKanchipuramViluppuram TiruvannamalaiWaidhan JangipurJhansiBulandshahrBudaun Ranchi Moradabad Bhadravathi CuddaloreBharatpur AzamgarhBrahmapurPilibhitDhulianSilchar Darbhanga Gaya o Visakhapatnam Karur UNNAMEDKhammamPuruliyaRoorkeeMorenaThanjavurSansarpurObra AmrohaBandaGhazipurMau Salem HassanKarimnagarChitradurgaKadapaMandya BhilwaraSatna NanparaBallia l l Vizianagaram KhalilabadUjjain Aurangabad Muzaffarpur GudiyathamBaberuKishangarhChitrakootKolarChandausi NaugarhGold DhamKhurja FieldsKatiharBhindMainpuriUnnaoSitamarhiEtawahRampur DoddaballapuraMuhammadabadMahuvaBeawar Rewa, ThakurdwaraMadhya Pradesh Bokaro

l Madurai 1.0 Nashik Bhadrak NajibabadUdupi Nawab GanjDausaMuhammadabadKorba UNNAMEDHathrasShahabad Piprakothi ChittorgarhRaichurChhatarpurPali Bankura Sidhauli Balrampur Kushinagar Khandwa Faizabad Tamluk Bijnor BaheriValsadDhaulpurShahganjTandaHindaun Khaga RanaghatKarauliBilariBisalpurTilharMusiriBasirhatDamohGopiganjKatni AkbarpurMungerChhapraBettiah Bihar Sharif Ludhiana Dhanbad SawaiGangapur MadhopurDatiaRatlamDurgapurBansiDewasShikohabadDhubriLalitpurArrahPurniaColonelganj Porbandar GunaRajsamandKaimganjRudrapurTinsukiaTezpurBarmerHaldwani SamastipurSandila AmbikapurAuraiyaEtmadpurBandikuiMohmadiChhindwaraOrai Deoghar Kharupetia , Burhanpur Mansa Nagar , Gorakhpur Vadodara SujangarhVidisha NarkatiaganjKannaujSivaniHabraPuriKundaAizawlAliganjRudauliJalalpurMaholi BanswaraHajipur Lauiyah Nandangarh Bhopal NanjanaguduAmethiShivpuriMadhubani Saharsa Giridih , Chandigarh Notun DighaNautanwaBalotra 1.5 GauriganjRasraMedininagarNawabganjChhibramauBhagalpur Meerut Motihari KasganjRafiganj Coimbatore SalonManjha MirganjBuxarSasaramLalganjDobhi Kozhikode KishanganjHazaribagh KalpiDehriSiwan

k Kollam Masaurhi Maharajganj 9 DungarpurFaridpur Agra Nawada Guruvayur P Joura Phalodi n Thiruvananthapuram Utraula Nagpur Alirajpur Tulsipur Bhinga

n Patna Indore BhongaonSingrauli

o Prayagraj Jaipur D LaharpurSidhi i Asansol Anand Nagar

a Jhumri Tilaiya

t Kochi Kanpur 1.5 r Lucknow G

Varanasi u 1.0 Surat Pune y 8 b y t i Hyderabad

Ahmedabad t i r i t

C Chennai Bengaluru C s

i 2.0

d 0.5 Mumbai 7 Kolkata

Complementary cumulative 2.5 Darjeeling 0.0 Delhi [New Delhi] 6 5.5 6.0 6.5 7.0 7.5 5.5 6.0 6.5 7.0 7.5 5.5 6.0 6.5 7.0 7.5 City population, log s City population, log s City population, log s Tehran NazarabadDogonbadan Nahavand TaleshEsfarayen BonabMianeh NurabadShadegan ZarrinshahrAlvand Bandar-eKhash Mahshahr 0.0 2.0 MarandDorud BehbahanFasa Ahar Shahin Shahr Iran (IRN) Bandar-eTorbat-e Anzali JamAligudarzShahreza Iran (IRN) Iran (IRN) Eslamabad-e GharbDarab SemnanKashmar SalmasBaneh ShirvanLahijan Shahr-eBam Kord MarivanKuhdasht IranshahrQuchan Ardakan Marvdasht JahromMahabad ShahrudKazerun Ilam Saveh ZabolRafsanjan = 1.08±0.11 SaqqezIzeh Shoushtar = 1.08±0.11 = 1.09±0.07 Sabzevar Maragheh Malayer PakdashtBagh-e Zahra Miandoab Khoy ) 11 Chalus Yasuj Najafabad Birjand s = 98035 Sirjan Qaem Shahr s = 98035 min Borujerd Amol min s Chabahar Gonbad-e Qabus Isfahan Pardis Kashan Ahwaz

( Neyshabur Bojnurd 0.5 1.5 Babol Sari Jiroft Gorgan Mashhad y r Abadan F Qazvin Khorramabad Zanjan Sanandaj Shiraz Tabriz Dezful Dezful

Qarchak g g Arak AbadanQarchak Qom g Bandar Abbas Ardabil Yazd Yasuj o o Pardis Gonbad-eGorgan QabusArak Kermanshah o Zahedan Babol Kerman l l 10 Sari Kerman Bagh-e Zahra Hamedan Rasht l Ardabil YazdUrmia

Rasht Izeh QaemAmol Shahr QazvinZanjan Hamedan DogonbadanShadeganBandar-e Mahshahr Pakdasht Kashan , , Behbahan ChalusShoushtar , 1.0 Najafabad NeyshaburKhorramabad 1.0 Kermanshah JiroftBojnurd Urmia Ilam BandarSanandaj Abbas k Torbat-e Jam Saveh Sirjan P

n Sabzevar Shahrud BorujerdMaragheh Qom Ahwaz Quchan Jahrom Zahedan n NazarabadShahin Shahr Khoy Zarrinshahr MarvdashtFasaShahrezaSemnan o KashmarMalayerArdakan D Marand Kazerun i Eslamabad-eShahr-e Gharb Kord LahijanMahabad a MianehAlvand ShirvanSaqqezKuhdashtRafsanjan t Shiraz Bam Bandar-e AnzaliAhar Miandoab r Baneh Tabriz Bonab Marivan DorudEsfarayen G Aligudarz

Talesh Zabol u

9 Birjand Nahavand

y Nurabad Salmas b y 1.5 t 0.5 Isfahan i t i r Mashhad i

t Darab Chabahar C C s i

d Iranshahr 8 2.0 0.0 Tehran Complementary cumulative

Khash

5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

42 Baghdad 0.0 Iraq (IRQ) 1.0 Karbala Kirkuk Iraq (IRQ) Iraq (IRQ) = 2.24±0.71 Hindiya Ad Diwaniyah = 2.24±0.71 = 1.03±0.04

) Najaf smin = 830298 Mosul smin = 830298 10.8 s ( y r F 0.5 0.5 Erbil g g

g Basra o o o l l

l 10.6

, , , k P n n

o 1.0 0.0 Baghdad D i a t r G

u 10.4

y b y t i

t Basra i r i Erbil t C C s

i 1.5 0.5

d 10.2

HindiyaNajaf Kirkuk Mosul Ad Diwaniyah

Complementary cumulative Karbala 2.0 1.0 10.0 6.0 6.2 6.4 6.6 6.8 6.0 6.2 6.4 6.6 6.8 6.0 6.2 6.4 6.6 6.8 City population, log s City population, log s City population, log s

Anzio Milan Cosenza 0.00 Lecco Pordenone Italy (ITA) Trento Italy (ITA) Italy (ITA) Massa Reggio Calabria 1.50 Reggio nell'Emilia Triest Parma Udine Rome = 1.03±0.17 MessinaTaranto = 1.03±0.17 1.00 0.06 Treviso Rimini = ±

) Como Modena 0.25 smin = 165011 Vicenza Cagliari smin = 165011 s SalernoMestre 11.0 ( 1.25 Verona y

r Pescara Naples F Brescia Bari g

g Turin g 0.50 Bologna Padua o o

o 1.00 l l Genoa l Catania

Florence

, , Bergamo , Palermo k

Bergamo P n 0.75 10.5 BolognaGenoa n Padua o

0.75 D i Brescia a Florence Palermo t

r Turin G Verona u

b 1.00 Mestre y

t Catania i Como t i r 0.50 i Rome VicenzaRiminiModenaPescara Bari t

C Treviso C s Reggio nell'EmiliaParma i Naples 10.0 Udine Cagliari 1.25 Trento Triest

d Lecco Salerno Massa 0.25 Pordenone Anzio Messina Taranto

Complementary cumulative 1.50 Reggio Calabria 0.00 Milan Cosenza 5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Matsuyama Tokyo 1.5 Kagoshima Kinokawa [Wakayama] 0.0 Fukuyama Tokushima Japan (JPN) Kanazawa Takamatsu Japan (JPN) Japan (JPN) Toyohashi Moriyama [Otsu] Utsunomiya Toyama = 1.17±0.22 Niigata Mito = 1.17±0.22 12.0 = 1.04±0.02 Mishima Kumamoto ) smin = 584761 Hamamatsu Shizuoka smin = 584761 s Naha Okayama

( Osaka [Kyoto] 1.0 Sendai Kitakyushu y 0.5 r F Hiroshima Sapporo Nagoya g g g o o Takasaki [Maebashi] 11.5 o l l Fukuoka l

, , ,

k 0.5

Nagoya P n 1.0 n

o Osaka [Kyoto] D i a t r G

11.0 Fukuoka y

b Takasaki [Maebashi] y t i t i Sapporo r 0.0 i Hiroshima t 1.5 Tokyo C Sendai Kitakyushu C s Okayama i Hamamatsu Mishima Shizuoka d Toyohashi Utsunomiya Naha NiigataMoriyamaKumamoto [Otsu] 10.5 Kanazawa ToyamaMito Fukuyama Takamatsu 0.5 Tokushima Complementary cumulative 2.0 Kinokawa [Wakayama] MatsuyamaKagoshima

6.0 6.5 7.0 7.5 6.0 6.5 7.0 7.5 6.0 6.5 7.0 7.5 City population, log s City population, log s City population, log s

Turkestan Almaty 0.0 Kazakhstan (KAZ) Atyrau Kyzylorda Kazakhstan (KAZ) 10.75 Kazakhstan (KAZ) 1.0 Oral Aktau = 1.93±0.50 TarazOskemen = 1.93±0.50 = 0.95±0.29 ) smin = 225331 Pavlodar Semey smin = 225331 s Nur-Sultan ( 10.50 Atyrau y r Karaganda F 0.5 Aktobe g g g o o

o 0.5 Nur-Sultan l l l

10.25

, ,

, Shymkent k P n Kyzylorda Aktobe n o D i

1.0 a Oral t r 10.00 G

u OskemenSemey y 0.0 Shymkent

b Almaty y t i t i r i t

C KaragandaAktau C

s Pavlodar

i 9.75

d 1.5

0.5

Complementary cumulative 9.50 TurkestanTaraz

5.4 5.6 5.8 6.0 6.2 5.4 5.6 5.8 6.0 6.2 5.4 5.6 5.8 6.0 6.2 City population, log s City population, log s City population, log s

Kimilili Nairobi Kilifi Kakuma 0.0 1.5 Ukunda Mtwapa Kenya (KEN) Nyahururu Busia Kenya (KEN) Kenya (KEN) Naivasha Machakos Embu DadaabMeru Bungoma Nyeri = 1.21±0.21 Hagadera Refugee Camp Wajir = 1.21±0.21 = 1.34±0.11 Kitengela El Wak ) smin = 60180 Kakamega Mandera smin = 60180 s Malindi Ongata Rongai 10 ( 0.5 1.0 Thika Kitale y r F Garissa Kisii Mombasa g g g Kisumu o o Eldoret o l l l

, , , 1.0 0.5 Nakuru ThikaGarissa k 9

P Nakuru n Mombasa

n Kisumu

o Eldoret D i a t r Malindi G

u Kisii

UkundaKilifi KakamegaNyeri

y Kitengela OngataKitale Rongai b y

t Wajir

i 1.5 0.0 Nairobi Nyahururu t i Hagadera RefugeeBusia Camp r i Naivasha MachakosMeru Mandera t El Wak C 8 EmbuMtwapa C s Kakuma i Bungoma Kimilili d 2.0 0.5 Complementary cumulative 7 Dadaab 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

43 Jeju-si 12.0 Seoul 0.0 Geoje-si Iksan-si South Korea (KOR) Gunsan-si Chuncheon-si South Korea (KOR) South Korea (KOR) Suncheon-si Wonju-si = 0.81±0.17 Dongducheon-si Mokpo-si = 0.81±0.17 = 1.00±0.04 Jinju-si Gumi-si ) smin = 185039 Pohang-si Jeonju smin = 185039 s

( 1.0 11.5 Changwon Cheonan-si y r F Cheongju-si g 0.5 g Ulsan g o o o l l Gwangju Daegu Busan l

Daejeon , ,

, 11.0

k Daejeon P n Gwangju

n 0.5 Daegu o D i a t

r Busan Ulsan

G Cheonan-si

1.0 u Cheongju-si

Changwon y b

y 10.5 t i t i r i Pohang-si

t Gumi-si Jeonju C C s i 0.0 Seoul Jinju-si d

Geoje-si Suncheon-siGunsan-siJeju-si 1.5 10.0 Wonju-si Iksan-siChuncheon-si Dongducheon-si Complementary cumulative Mokpo-si

5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

Derna Tripoli 0.00 Libya (LBY) Sirte Libya (LBY) Libya (LBY) 1.00 Marj 10.25 Bayda = 1.04±0.29 Sabha = 1.04±0.29 = 1.05±0.04 ) smin = 74873 Tobruk smin = 74873 s 0.25 Zliten ( 10.00 y r F 0.75 Benghazi Qabilat alturki g g g Al Zawiyah o o o l 0.50 l

l 9.75 Misrata

, , , 0.50 Qabilat alturki Al Zawiyah

k Misrata P n n o D i

a 9.50 t 0.75 Benghazi r G

y 0.25

b Tobruk y t i t i r i Sabha

t 9.25 Bayda Zliten 1.00 C C s

i Marj

d 0.00 Tripoli

9.00 Derna 1.25 Complementary cumulative Sirte 0.25 8.75 5.00 5.25 5.50 5.75 6.00 5.00 5.25 5.50 5.75 6.00 5.00 5.25 5.50 5.75 6.00 City population, log s City population, log s City population, log s 0.00 1.00 Badulla 11.0 Colombo Sri Lanka (LKA) Ambalangoda Trincomalee Sri Lanka (LKA) Sri Lanka (LKA)

Batticaloa = 0.87±0.26 Kalmunai = 0.87±0.26 = 1.19±0.05 ) 0.25 smin = 94688 smin = 94688

s 0.75 Jaffna ( Galle 10.5 y r F g g g 0.50 0.50

Matara o o o l l l

, ,

, Kandy k P n 0.75 10.0 Kandy n 0.25 o D i a t r G

u Matara y b y

1.00 t i t i 0.00 Colombo r i 9.5 Galle t C C s Jaffna i

d 1.25

0.25 Kalmunai Batticaloa Ambalangoda 9.0 Complementary cumulative 1.50 Trincomalee 0.50 Badulla

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s Casablanca Ain El Aouda Bejaad Al Hoceima Ain Taoujdate Sidi Yahya du Gharb Ouled-Teima Azrou 0.0 Midelt Zagora Morocco (MAR) Chefchaouen Kasba Tadla Morocco (MAR) Morocco (MAR) Taroudant Souk Sebt Oulad Nemma Essaouira Melilla Ceuta Deroua 10.5 Ouazzane Azemmour Tiznit Ben Guerir Youssoufia Sidi Kacem = 0.98±0.13 Tiflet Sidi Bennour = 0.98±0.13 Souk El Arba Du Gharb Ouarzazate = 1.12±0.05 Sefrou Oued Zem

) Errachidia 1.5 Taourirt Fkih Ben Salah smin = 53133 El Kelaa des SraghnaGuercif smin = 53133 Rabat s Guelmim Larache Khenifra

( Sidi Slimane BerrechidBerkane y r Khemisset

F Settat 10.0 0.5 Taza AgadirMarrakesh Ksar el-Kebir Tangier g g Khouribga Beni Mellal Fez g El Jadida Nador o o

o Meknes l l Safi Tétouan l KenitraTétouan

, , 1.0 , Kenitra Oujda 9.5 Oujda k P

n Safi Meknes Berrechid

n Tangier Khouribga o Beni Mellal D i 1.0 a El Jadida t Settat Nador r Agadir Ksar el-Kebir G

Marrakesh Melilla u Khemisset OuedFkih BenZem Salah y

b 9.0 y Berkane t Souk El Arba Du Gharb Larache i t i El Kelaa desSidi Sraghna Slimane r 0.5 Fez i Ceuta DerouaTaourirtSidi Kacem Taza t Youssoufia AzemmourSefrouGuelmim C Sidi Bennour Kasba TadlaTiflet

C Ben Guerir s Rabat ErrachidiaOuarzazateKhenifra

i Souk Sebt Oulad Nemma Tiznit 1.5 MideltSidi Yahya du Gharb Guercif d Bejaad Azrou Taroudant 8.5 Ouled-TeimaChefchaouenAin El Aouda

Ain TaoujdateOuazzaneAl Hoceima

Essaouira Complementary cumulative 0.0 Casablanca Zagora 8.0 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Aguascalientes 11.75 Mexico City 0.0 Mexico (MEX) Mérida Cuernavaca Mexico (MEX) Mexico (MEX) 1.0San Luis Potosí Querétaro = 1.29±0.33 Torreon Ciudad Juárez = 1.29±0.33 11.50 = 1.14±0.11 )

s smin = 939745 León Tijuana smin = 939745 ( y r

F Toluca 0.5 Puebla 11.25 g g

g Monterrey o o o l l

l 0.5 Monterrey

, , , 11.00

k Guadalajara P n Guadalajara n o D i a t

1.0 r G

u 10.75

y b y t

i 0.0 t i Mexico City r i Tijuana t

C Toluca Puebla C s 10.50 Torreon i Ciudad Juárez d QuerétaroLeón 1.5 10.25San Luis PotosíAguascalientes Mérida

Complementary cumulative 0.5 Cuernavaca

6.00 6.25 6.50 6.75 7.00 7.25 6.00 6.25 6.50 6.75 7.00 7.25 6.00 6.25 6.50 6.75 7.00 7.25 City population, log s City population, log s City population, log s

44 10.0 Bamako 0.0 Mali (MLI) Mopti Bougouni Mali (MLI) Mali (MLI) 1.0 Kita Timbuktu = 1.08±0.29 San Niono = 1.08±0.29 = 1.10±0.13 ) s = 54841 GaoSocoura s = 54841 s min min 9.5

( Kayes y r Koutiala F 0.5 g g g 0.5 Ségou o o o l l

l 9.0

Sikasso

, , , k P n KayesSégouKoutialaSikasso n o

1.0 D

i Gao a t r

G 8.5 0.0 u Bamako Mopti y b

y Kita t i t i r i San t C C s Bougouni i 1.5 8.0 Timbuktu Socoura d 0.5 Complementary cumulative 7.5 Niono 2.0 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

2.0 Ngathaingchaung Yangon TaungdwingyiDedaye MadayaTatkon HopongShwegyin UNNAMEDBhamoNyaungNattalin Shwe 0.0 PayathonzuHlegu YamethinMyan Aung Kyaukpadaung Kyaukme Myanmar (MMR) KawthaungMuse KyaiktoThongwa Myanmar (MMR) Myanmar (MMR) TwanteMudon NyaunglebinMogok AungpanHsipaw KayanPyu ShweboDaik-U Taungup Pyawbwe LoikawMinbya MyawaddyWundwin Tan Soe Tachileik MaubinDawei = 1.48±0.17 TangyanAunglan Pyapon = 1.48±0.17 Za Lun Pakokku = 1.08±0.12 Danu Phyu Sittwe ) UNNAMED Pantanaw 10 Mandalay 1.5 Kengtung Myingyan s = 73994 Thaton Magway s = 73994 s min Wakema Labutta min 0.5 Kyaiklat Pyay

( Taungoo Mawlamyinegyunn Pyin Oo Lwin Meiktila y r Hpa-An F Monywa Myaungmya Hinthada Mawlamyine g g Myitkyina

g Pathein 1.0 Lashio NaypyidawLashioMyeikTaunggyi Monywa o o Taunggyi

o 1.0 Naypyidaw l l Hpa-AnMeiktilaMawlamyineBago

l Pyin Oo Lwin

Shwebo Magway

Myeik , , Bago 9 KengtungTaungoo , MuseMyawaddyKyaukme Pakokku Kyaukpadaung PyawbweTangyanTan Soe Myingyan Pathein

k Aungpan Payathonzu Wundwin LoikawDawei Pyay P n Tatkon Hsipaw 0.5 Bhamo Myitkyina Mudon Sittwe n HleguNyaung ShweThongwa Kyaiklat o Mandalay D Twante TachileikMinbya

i 1.5 Yamethin Pyapon a UNNAMED Thaton t AunglanKayan TaungdwingyiHopong Hinthada r NyaunglebinPyu Myan Aung G

Daik-U

u NattalinKyaikto Maubin Madaya y Kawthaung MogokDanu Phyu

b Dedaye y t

i 0.0 t i Yangon r 2.0 i 8 t Shwegyin Taungup Myaungmya C C s Labutta i UNNAMEDZa Lun

d Mawlamyinegyunn

0.5 Pantanaw 2.5 Wakema Complementary cumulative 1.0 7 Ngathaingchaung 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

1.0 Maputo 0.0 Mozambique (MOZ) Tete Mozambique (MOZ) Mozambique (MOZ) Quelimane Nacala = 1.66±0.55 = 1.66±0.55 = 1.66±0.32

) Chimoio s = 221384 Beira s = 221384 s min min 9.5 ( y r

F 0.5 0.5 Milange g g g Nampula o o o l l l

, , , 9.0 k P n n o D

i 0.0

a Maputo t

1.0 r Nampula G

Beira u

y b y t i t i r i 8.5 Chimoio t C C s i

d 0.5 Quelimane 1.5 Tete

Complementary cumulative 8.0 Milange Nacala

5.4 5.6 5.8 6.0 6.2 6.4 5.4 5.6 5.8 6.0 6.2 6.4 5.4 5.6 5.8 6.0 6.2 6.4 City population, log s City population, log s City population, log s 1.0 11.5 Kuala Lumpur 0.0 Malaysia (MYS) Kota Bharu Melaka City Malaysia (MYS) Malaysia (MYS) Kuching = 1.50±0.47 George Town = 1.50±0.47 = 1.21±0.26

) Kota Kinabalu smin = 495448 smin = 495448 s Ipoh ( y r

F 11.0 0.5 0.5 Bukit Mertajam g g g Singapore o o o l l l

, ,

, Singapore

k Bukit Mertajam P n n

o George Town

D 10.5 i a t 0.0 Kuala Lumpur 1.0 r G

u Kuching y b y t i t i Melaka City r i

t Ipoh

C Kota Kinabalu C s i 10.0 d

1.5 0.5

Complementary cumulative Kota Bharu

5.75 6.00 6.25 6.50 6.75 5.75 6.00 6.25 6.50 6.75 5.75 6.00 6.25 6.50 6.75 City population, log s City population, log s City population, log s

Aguié Niamey 0.00 Niger (NER) Keïta Tchadoua Niger (NER) 9.0 Niger (NER) Gazaoua Filingué 1.25 Baléyara Agadez = 0.99±0.21 = 0.99±0.21 Tahoua Galmi Magaria = 1.22±0.17 Zinder ) smin = 56927 Diffa Tanout smin = 56927 Maradi s 0.25 ( MadaouaTessaoua y r 8.5 Dosso F 1.00 Dogondoutchi Tillabéri g g g o o

o Gaya l 0.50 l

l Birni N'Konni

Agadez Tessaoua

, , Gaya , 0.75 Dogondoutchi Birni N'Konni k P n Dosso 8.0 Galmi n

o Maradi D i

a Diffa t 0.75 Magaria r G

Aguié u Tchadoua Tillabéri

y 0.50

b Zinder y

t Filingué i t i r i

t 7.5 1.00 C Tahoua Baléyara C s i Gazaoua

d 0.25 Madaoua Keïta 1.25

Complementary cumulative 7.0 0.00 Niamey Tanout 4.75 5.00 5.25 5.50 5.75 6.00 4.75 5.00 5.25 5.50 5.75 6.00 4.75 5.00 5.25 5.50 5.75 6.00 City population, log s City population, log s City population, log s

45 2.5 MaigatariJakuskoDambanUba IdembiaRano UNNAMEDDundewaBagudoGulak OgwashiBunzaMafaraAbiriba Uku ObinombaDamaskImorunKoi Lagos Isa Ilawe-Ekiti TudunJemmaKabba WadaIkoloare IgbetiOwodeIdo Ekiti BokkosAmabaKamba BununuIgboMayaGaya Ora Belwa Dass DongaAdikpoIdanreIlaro JahunEffonNuman Alaiye 0.0 EkogboroDukkuAwgu ShendamAgaieTudunMonguno Wada Kpam Egbe GamawaBiantubuShellemGiade BuniIkotAlesa AkpanataYadi EruwaDamaturu AwalaTitun SongGwarzoNshebori 11 Nigeria (NGA) Mashi Shema Nigeria (NGA) Nigeria (NGA) WurnoGummiDoma ArgunguElu BirninAkwangaGombi Kudu GumelKokoNew Bussa NassarawanDoka EggonAshafaBichi TarzohoEpeAramoko-Ekiti LangtangKatsiroOmema DengiDambartaOnne UduagbaIshiaguEnwan IgbohoGembuDamagum Enugu-EzikeEjigboIdah BagaOtunTakum Ijero-EkitiOmokoUromi UmundiBoriBajoga MonguDikwaIkole Arna MalumfashiEkpomaUgbokolo EbuteTugaAyangba NamaIrele ZinnaSamunakar UmuKaltungo Lawlaw KishiMisauIrri OkitipupaIkomNingi MichikaOzizza = 1.18±0.08 KontagoraKeffi EdieneOmenamaKaura-Namoda = 1.18±0.08 IleGoniri Oluji IwereIjagbo Offa = 1.22±0.04 OkukuOyigbo Lepuolekazaure UgepIkot Ekpene NguruIkerre AgborWudil Onitsha ) Lawmu UmunzeOre Omuo-ObadoreSangasumi GashuaMasuri HadejiaAnkpa JaudiriZuru Gwoza GeidamGubio AuchiEket Amaigbo Ibadan 2.0 VandeSagamu Ikya BiuLokoja Biliri Damboa s s = 87687 Funtua s = 87687 Katsina Ala min ZakiUghelli Biam IseWukari min BugamaDutse IwoJega AzareIkire BonnyOwo ( 0.5 Gamboru Lafia Kano Birnin Kebbi Obudu OtukpoJalingo KrakramaMogho Benin CityAbuja IseyinBida Iragbiji

Mubi y r Port Harcourt AbakalikiOron Ijebu-Igbo Kaduna F Potiskum GusauSaki 10 IkoroduEnuguAba Gwagwalada Makurdi KafanchanGombe Sapele Abeokuta OsogboWarri Oboama Nguru Ijebu-Ode IfeUmuahia Ilorin Ogbomosho Gboko ZariaUyo g

g Daura Nsukka Minna Akure g Ilesa Ado Ekiti IlesaOndoYenagoa Oyo SokotoOkene Ondo Igbe Ogbomosho KatsinaIgbe Jos Yenagoa Oyo Ijebu-Ode AdoCalabar Ekiti o o Bauchi Akure Iseyin Nsukka Maiduguri

o 1.5 Katsina Jimeta

l Saki l Ijebu-IgboSapele Calabar Sokoto Zaki Biam Gwagwalada l Abeokuta OwoKafanchanIragbijiDauraMakurdiUmuahiaGbokoMinnaBauchi

Ife 1.0 Okene IkerreVandeKatsina IkyaUghelli AlaIkire Jimeta Obudu Ikorodu Iwere IseIwo

Osogbo , AbakalikiBida , Kishi Zaria Ijero-EkitiEbute IreleIle OlujiSagamuOreAgborAuchiOtukpo Gusau , Warri Uduagba Igboho Bugama Gombe Uyo Idanre EjigboOkitipupaEkpomaZuru OmokoUromi Ugbokolo Ankpa Lafia k Maiduguri Koi Umundi OmenamaBirninOmuo-Obadore KebbiFuntua Jos Enwan Otun LokojaWukariPotiskum P Aramoko-EkitiKontagora n IlaroOwode Ikom JalingoJega Ilorin 9 EkogboroEnugu-EzikeIshiaguIdahElu Ijagbo OffaEket OronMogho Lawmukazaure UmunzeIrri BiliriAzare Enugu AwalaIgbeti Ayangba Ediene

n Ilawe-EkitiUmu Lawlaw Hadejia Eruwa TitunIdo Ekiti NsheboriEpe Keffi OzizzaUgepDutse Oboama Nguru Ikole ZinnaGeidam o Ogwashi Uku Okuku Lepuole Wudil Aba NassarawanAdikpoBiantubuMashi Eggon KokoNewKaltungo Bussa Bonny D LangtangShemaDengiMalumfashiSamunakarNguruKaura-NamodaBiu i KpamEffonNuman AlaiyeTarzoho MafaraKabbaIsaDamaturuObinomba Dambarta a Abiriba Takum Sangasumi Mubi 1.0 Benin City Tudun WadaKambaIkolo GummiKatsiroEgbeBichiTuga OnneNama Gashua Shendam t Port Harcourt DukkuWurno BagaBajogaMisauGubioIkot Ekpene DokaIgbo Ashafa Ora Gumel Bori Damboa Jemma are GombiAgaieMongu Arna Gamboru r AwguBirnin Kudu Argungu Masuri JakuskoDambanUNNAMED Gwarzo G DamaskBuni Yadi

Gwoza 1.5 AmabaDundewaAkwanga Ningi Abuja Idembia u Kaduna BokkosMaigatari DamagumOyigboJaudiri BununuDomaIkot Dass Akpanata GamawaShellem Bunza ImorunGayaTudun Wada Michika Krakrama y Rano Maya Belwa Gembu b

y Bagudo OmemaAlesa

t Jahun Amaigbo Uba i Ibadan GulakGiade t i Monguno r i Song

t 8

C Donga Dikwa

0.5 C

s Onitsha i 2.0 Kano d

0.0 7

Complementary cumulative Lagos 2.5 Goniri

5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

Jinotega 10.00 0.0 Managua Nicaragua (NIC) Jinotepe Esteli Nicaragua (NIC) Nicaragua (NIC)

= 1.73±0.58 Matagalpa = 1.73±0.58 = 1.35±0.28

) Masaya 9.75 s = 102401 s = 102401 s min min ( 0.5 Chinandega y r F 0.5 León g g 9.50 g o o o l l l

, , ,

k 9.25 P n 0.0 n Managua o D

i 1.0 a t

r Masaya G

u 9.00 y b y t

i Jinotepe t i r i t C C s 0.5 i Chinandega León 1.5 8.75 MatagalpaEsteli d

Complementary cumulative 8.50 1.0 Jinotega 5.00 5.25 5.50 5.75 6.00 5.2 5.4 5.6 5.8 6.0 5.2 5.4 5.6 5.8 6.0 City population, log s City population, log s City population, log s

Roosendaal Rotterdam [The Hague] 0.0 Apeldoorn 11.25 Netherlands (NLD) Venlo Hilversum Netherlands (NLD) Netherlands (NLD) Almere Heerlen = 1.38±0.29 Veenendaal Alkmaar = 1.38±0.29 = 1.06±0.04 Zwolle Maastricht ) smin = 211662 Breda Groningen smin = 211662 s 11.00 Amsterdam

( 1.0 Tilburg Amersfoort y r F 's-Hertogenbosch g 0.5 g Enschede g o o o l l Eindhoven 10.75 l

Utrecht , , Utrecht ,

k Eindhoven

P Arnhem n

n 0.5 Arnhem o D i a

t 10.50 r Amsterdam G

u 1.0 Groningen y Enschede b y

t 's-Hertogenbosch i t i r i Tilburg Amersfoort

t Breda C

C 10.25 s Alkmaar i Veenendaal 0.0 Zwolle Maastricht

d Rotterdam [The Hague]

Hilversum 1.5 Heerlen 10.00 VenloAlmereRoosendaal Complementary cumulative Apeldoorn

5.50 5.75 6.00 6.25 6.50 5.50 5.75 6.00 6.25 6.50 5.50 5.75 6.00 6.25 6.50 City population, log s City population, log s City population, log s

Birtamod Kathmandu 0.00 Nepal (NPL) 1.00 Nepal (NPL) Nepal (NPL) Tilottama Damak 10.0

= 0.74±0.22 Hetauda = 0.74±0.22 = 1.16±0.11

) Nepalgunj smin = 50827 smin = 50827 s 0.25 0.75 ( Bharatpur y r F Pokhara 9.5 g g g Pokhariya o o o

0.50 l l 0.50 l Pokhariya Itahari

, , , k P n Itahari 9.0 n o D i

0.75 a t 0.25 r Pokhara G

Bharatpur y b y t i t i r i

t 8.5

1.00 C C s 0.00 Kathmandu i Damak Birtamod d

Nepalgunj Hetauda 1.25 0.25 8.0 Complementary cumulative Tilottama

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s 1.0 0.0 Al Amerat Oman (OMN) Al Ain Oman (OMN) Al GhubraSeeb = 1.36±0.43 Balad Bani Bu Ali = 1.36±0.43 Nizwa 10.2 ) s = 75594 0.8 s = 75594 Al Ain s 0.2 min min ( Ibri y r F Sohar 10.0 Sohar g g g 0.6 Salalah o 0.4 o o l l l

, , Salalah , 9.8 k P n 0.4 n

o 0.6 D i Seeb a

t Ibri r G

u 9.6 Balad Bani Bu Ali y b y t i t 0.8 i 0.2 r i t

C Nizwa C s 9.4 i d 1.0 0.0 Al Ghubra Oman (OMN) 9.2

Complementary cumulative Al Amerat = 1.17±0.27 1.2 0.2 5.0 5.2 5.4 5.6 5.0 5.2 5.4 5.6 5.0 5.2 5.4 5.6 City population, log s City population, log s City population, log s

46 Mithankot Karachi AttockBalakot 11 Kahror Pakka Toba Tek Singh DeraAhmadpur Murad Jamali East PirDir Murad DinaMoro NawabshahUsta Muhammad 0.0 HavelianKahuta KabirwalaThul Hasilpur Layyah Pakistan (PAK) KhanpurTopi MangowalAlipur Pakistan (PAK) Pakistan (PAK) 2.0 PindiArifwala Gheb MundaPasrur Qala Fateh Jang Battagram Lahore ShujaabadKhanewal KandhkotSowari Tando MuhammadLodhran Khan BindiDera Ismail Khan MundaTarakai Dadu Khwazakhela GojraPakpattan ShahdadpurSangla Hill Turbat Hangu KhuzdarKhairpur = 1.02±0.10 Rana KhemrajDaharki UNNAMED = 1.02±0.10 = 1.03±0.06 Peshawar Kotli ChunianChiniot Rawalpindi [Islamabad] ShahdadkotPanjgur Khatako Shah Faisalabad ) Tando Adam Sadiqabad Shikarpur Burewala Timergara Mandi Bahauddin Multan 121657 Hafizabad Muzaffargarh 121657 s = Jacobabad Takht-e-Bhai s = Hyderabad s min min Jampur Narowal Gujranwala Muzaffarabad Shakar Garh Jhang Chakwal

( Sialkot Rahimyar Khan Swabi Badin Tando Allayhar 10 y r Kohat Kharian Mardan F 0.5 Pattoki Chakdara Bahawalpur Abbottabad 1.5 Dera Ghazi Khan Mirpur Khas Sargodha Kamra Gujrat Sahiwal Sukkur Jhang SahiwalBannu g g Jhelum Ahmadpur East g Larkana Kamra Okara Swat City Bahawalpur Haripur Kasur Wah Cantonment Okara Khanpur Quetta o

o Kharian Kasur Nawabshah Sukkur o MuzaffargarhLarkana Haripur l l Wah Cantonment Kohat Bannu Dera Ismail KhanChiniot Pattoki Jhelum l Khanewal ShikarpurHafizabadMuzaffarabad

Quetta Tando Adam SadiqabadBurewalaDeraMandi Ghazi BahauddinChakwal Khan Sargodha Hasilpur Chakdara

, Daharki , TobaPirShahdadpur Murad Tek Singh JacobabadChunian NarowalSwabi

, Swat City ShujaabadLodhran Kotli Mirpur Khas Gujrat GojraPanjgur Tando Muhammad AlipurKhanPasrur Rahimyar Takht-e-BhaiKhan

k Pakpattan Kahror PakkaKahuta TarakaiMunda QalaKhairpurKhatako Shah

Abbottabad P Arifwala Shakar Garh

n Mardan 9 Dina Jampur 1.0 Mithankot BalakotUstaKabirwala Muhammad Timergara n 1.0 TopiKhwazakhelaLayyah Hyderabad Sialkot Sangla Hill o Havelian Mangowal UNNAMED D i BattagramSowari a Dera Murad Jamali t Multan Bindi r Gujranwala ThulMundaShahdadkot Fateh Jang G Turbat Rana Khemraj u Moro Tando Allayhar

Hangu y Rawalpindi [Islamabad] b

y Badin t Faisalabad i Kandhkot t i r i Attock t 1.5

C 8

0.5 C s Peshawar Dir i Lahore Dadu d

Khuzdar Complementary cumulative 2.0 Pindi Gheb 0.0 Karachi 7

5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s Lima 0.0 Peru (PER) Juliaca Tacna Peru (PER) Peru (PER) 1.0 Chimbote Pucallpa

= 1.21±0.35 Huancayo Cusco = 1.21±0.35 = 1.29±0.05

) 11.0 smin = 258428 smin = 258428 s Piura Chiclayo ( y r F 0.5 0.5 g g Arequipa g o o

o Trujillo 10.5 l l l

, , , k P n Arequipa n o

1.0 D i a

t 0.0

r Lima G

u 10.0 Trujillo y b y t i t i r i Chiclayo t C C s HuancayoPiura i 1.5 Cusco d 0.5 Tacna 9.5 Chimbote Pucallpa

Complementary cumulative Juliaca 2.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

Calapan Quezon City [Manila] IrigaCalbayog 11.5 Candelaria Isabela City BayombongSorsogon City 0.0 ValenciaPuray Laoag Daet Philippines (PHL) SurigaoSan CityJose Panabo Philippines (PHL) Philippines (PHL) Jolo Ormoc Bayambang Digos Santiago Kalibo Vigan Lemery Koronadal Roxas City Tabaco Ozamiz Dipolog Tuguegarao = 0.91±0.12 Tagbilaran Urdaneta = 0.91±0.12 Pagadian Gapan = 1.12±0.04 Puerto Princesa Dumaguete

) 1.5 San Fernando Balanga 11.0 Butuan Legazpi smin = 77409 Marawi smin = 77409 s Tagum Cotabato Lucena

( Tacloban 0.5 OlongapoCabanatuan y r Tarlac City F Iligan Naga General Santos g g Zamboanga City 10.5

g Batangas City o o 1.0 Baguio Bacolod o l l

l Angeles [San Fernando] Cebu City

Iloilo City Cagayan de Oro , , , 1.0 k

Dagupan P 10.0 n Davao City

n Dagupan o

D Davao City

i Iloilo City a t 0.5 Bacolod r Angeles [San Fernando] Batangas City Cagayan de Oro G

Cebu City Baguio 9.5 Olongapo y b y t CabanatuanTarlac CityGeneralZamboanga Santos City i 1.5

t Lucena i r i MarawiIligan Balanga t TaclobanTagum C Puerto PrincesaButuanGapanDumaguete

C TagbilaranSan FernandoCotabato s Lemery Naga i 0.0 Quezon City [Manila] 9.0 KaliboPagadianRoxasTuguegaraoUrdaneta City San Jose Dipolog Legazpi d Candelaria Ozamiz

SantiagoJolo BayambangKoronadalVigan Surigao City LaoagDigos 2.0 Calapan OrmocTabaco IsabelaValencia City Panabo Bayombong Calbayog Puray Daet Complementary cumulative 8.5 0.5 IrigaSorsogon City

5.0 5.5 6.0 6.5 7.0 7.5 5.0 5.5 6.0 6.5 7.0 7.5 5.0 5.5 6.0 6.5 7.0 7.5 City population, log s City population, log s City population, log s

Pi a Tczew Leszno Gronowo Górne 0.00 Tomaszów Mazowiecki Warsaw Gniezno Poland (POL) Lubin Ostrów Wielkopolski Poland (POL) Poland (POL) widnica 11.0 Wejherowo Ostrowiec wi tokrzyski Katowice 1.5 Legnica Inowroc aw = 0.88±0.14 KoszalinJelenia Góra = 0.88±0.14 Gorzów Wielkopolski Kalisz = 1.04±0.04 Olsztyn ) P ock smin = 87293 Opole Nowy S cz smin = 87293 s 0.25 Wa brzych Radom ( Toru Kielce y r F Bia ystokRzeszów 10.5 WroclawKrakow Cz stochowa Szczecin Gdansk g g Poznan g Bydgoszcz Bielsko-Bia a o o o l 0.50 l Rybnik Lodz

l 1.0

RybnikLublin Bielsko-Bia a , ,

, Radom Cz stochowaBydgoszcz Lublin k Wroclaw P n Poznan Wa brzych

n 10.0 o

D P ock Bia ystok i a

t 0.75 Gdansk Toru Rzeszów r Kielce

Lodz G

u Opole

y Nowy S cz

b Olsztyn

y Kalisz t

i Legnica t i Gorzów Wielkopolski Jelenia Góra r 0.5 i Krakow Lubin t widnica

C 9.5 Inowroc aw Leszno GnieznoOstrów Wielkopolski C

s 1.00 Ostrowiec wi tokrzyski Wejherowo i Warsaw Pi a Tczew d Tomaszów Mazowiecki Gronowo Górne

1.25 Szczecin Complementary cumulative 9.0 Koszalin 0.0 Katowice 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Hoechang Chonnae Kosan 10.0 Kowon T Chongnamktaedonggok 0.0 TaeanJi'an Yonan Hoeryong North Korea (PRK) UNNAMEDUnchon Kangdong North Korea (PRK) Tongrim Jonchon Sinpho Jaeryong Pyongyang Kumya Anak Kujang Sukchon Hamju Kusong Phyongwon Koksan Jongphyong Rason = 1.51±0.20 1.5 Musan Sinpyong = 1.51±0.20 Sonchon Hongwon Dandong Omong-ri ) Ryongchon Munchon Tokchon smin = 76818 Sinchon Jongju smin = 76818 9.5 s Kaechon Sillyeonpo Huichon ( 0.5 Hyesan Kimch'aek y r Songnim F Kaesong Ch ngjin Haeju Kanggye Hamhung g g

g 1.0 Anju Nampo o o 9.0 o Dandong

P'y ngs ng l l

l Hyesan

Anju Sariw n , 1.0 , Musan Haeju , Wonsan Sariw n Kaesong k KanggyeNampo

P P'y ngs ng n WonsanSongnim

n Jongju Hoeryong o Kimch'aek 0.5 Ch ngjin D i

a 8.5 T ktaedonggok t Hamhung Tokchon r Munchon SillyeonpoKaechon

G Rason

u Huichon

Sinchon Taean SinphoKujang SinpyongHongwon 1.5 y Jaeryong Ji'an Yonan Kusong b y t

i Kowon ChongnamOmong-ri t i Tongrim Jonchon r i HoechangSonchon Ryongchon t Jongphyong C 0.0 Pyongyang 8.0 C s Chonnae PhyongwonKumya i Hamju Unchon KangdongAnak d

2.0 Kosan Sukchon North Korea (PRK) 0.5 7.5 UNNAMED Complementary cumulative = 1.28±0.11 Koksan 2.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

47 Buz u Bucharest 0.0 Olteni a Baia Mare 11.0 Romania (ROU) Sibiu Br ila Romania (ROU) Romania (ROU) AradBac u = 1.62±0.35 Târgu Mure Pite ti = 1.62±0.35 = 1.41±0.08 Oradea Gala i ) smin = 136262 1.0 Craiova Ploie ti smin = 136262 s

( Brasov Ia i y r F 0.5 10.5 g g Constan a g Timi oara o o o l l l

, , 0.5 Cluj-Napoca , k P n C l ra i n

o 1.0 D Cluj-NapocaConstan a Timi oara i 10.0

a Brasov t r G

u Ploie ti y b y

t Craiova i Ia i t i 0.0 Bucharest Pite ti C l ra i r i Oradea

t Sibiu

C Arad C s

i 1.5 Târgu Mure Gala i 9.5 Br ila

d Bac u

Buz u 0.5 Baia Mare Complementary cumulative Olteni a 2.0 5.25 5.50 5.75 6.00 6.25 5.25 5.50 5.75 6.00 6.25 5.25 5.50 5.75 6.00 6.25 City population, log s City population, log s City population, log s

Ulan-Ude Moscow IvanovoBryansk BelgorodKirov Kaliningrad 0.0 Cheboksary Stavropol Novokuznetsk Lipetsk Russia (RUS) Ryazan Tula Russia (RUS) Russia (RUS) Kemerovo Orenburg Naberezhnye Chelny Yaroslavl 1.5 Khabarovsk Astrakhan Penza Ulyanovsk Izhevsk Tomsk = 1.60±0.24 Irkutsk Barnaul = 1.60±0.24 = 1.27±0.08 Tyumen Tolyatti

) Vladivostok Makhachkala Krasnodar 11.5 smin = 442432 Saratov smin = 442432 s Krasnoyarsk UfaVoronezh ( 0.5 Omsk Kazan y r F 1.0 Chelyabinsk Samara Tyumen g g

g Rostov-on-Don Nizhny Novgorod Saint Petersburg o o

o 11.0 l l Volgograd

l Yekaterinburg

Barnaul

1.0 , , ,

k 0.5

Novosibirsk P n

n Saint Petersburg o D i a Krasnoyarsk t

r 10.5 G

NovosibirskYekaterinburg u 1.5

Kazan Samara y

b Volgograd

y Omsk Nizhny Novgorod t

i 0.0 Chelyabinsk t i Moscow Rostov-on-DonUfa r i VladivostokTomsk Krasnodar t NaberezhnyeKhabarovsk ChelnyIrkutsk Tolyatti C SaratovVoronezh C

s Kemerovo Novokuznetsk i Belgorod IzhevskOrenburg Lipetsk 2.0 10.0 Yaroslavl

d Kaliningrad

Ryazan Ulyanovsk 0.5 TulaMakhachkalaAstrakhan Penza Cheboksary Stavropol KirovUlan-Ude Bryansk Complementary cumulative 2.5 9.5 Ivanovo

6.0 6.5 7.0 6.0 6.5 7.0 6.0 6.5 7.0 City population, log s City population, log s City population, log s Riyadh 0.00 Al Qurayyat Mad nat Yanbu a in yah 11.5 Saudi Arabia (SAU) MuhayilAbu Arish Saudi Arabia (SAU) Saudi Arabia (SAU) Sabya Sakaka = 0.87±0.17 Arar Abha = 0.87±0.17 = 1.02±0.03 1.25 Al Hawiyah Najran Jeddah ) smin = 142452 Al-Kharj Hafar Al-Batin smin = 142452 s 0.25 Hail Unaizah Dammam ( y r F 1.00 Buraydah Al Jubayl Industrial City 11.0 g g Mecca g Tabuk Khamis Mushait o 0.50 o o MedinaHofuf l l

l At Taif

Hofuf , , , k P n 0.75 At Taif Medina Khamis Mushait n o

0.75 D i Mecca 10.5 a Tabuk t r G Buraydah u Muhayil Al Jubayl Industrial City

y Unaizah

b 0.50 y

t Dammam i t i Hail r i Najran 1.00 Al-Kharj Hafar Al-Batin t Mad nat Yanbu a in yah C Jeddah Arar C s

i Al Hawiyah 10.0 Sakaka d 0.25

Al Qurayyat 1.25 Abu Arish Abha Complementary cumulative 0.00 Riyadh Sabya 1.50 5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Muglad UNNAMED UNNAMED UNNAMED New Halfa Ed Dueim Khartoum 0.0 Dilling Al Hilaliyah Sudan (SDN) Al Quwaysi Mellit Sudan (SDN) Sudan (SDN) ad-Damer Abu Jibeha Tandalti Jebel Aulia 1.5 Sennar UNNAMED Ghubaysh Shendi = 1.33±0.20 Atbara Al Rahad = 1.33±0.20 Rabak Geneina 10 = 1.21±0.16 Barah ) Giad Industrial Complex Saraf `Umrah smin = 133795 Ad-Damazin En Nahud smin = 133795 s Al-Qadarif Kassala ( 0.5 Kas Umm Ruwaba y r Al Hasahisa Al-Ubayyid F Kosti Ad Da'einAl-Fashir 1.0 Al ManagilAd Da'ein WadKosti Madan Port Sudan g g ad-Damer Nyala g Al-Qadarif Kassala Wad Madan Port Sudan Atbara o Geneina o En Nahud Umm Ruwaba o l l Ad-Damazin

l 9 Rabak Kas

UNNAMED

Buram Buram , ,

, Al-Fashir NewGiad Halfa Industrial Complex MugladGhubaysh Al Hasahisa

1.0 k Mellit Jebel Aulia P n Dilling Al Rahad UNNAMED Al Managil n 0.5 Barah o Al-Ubayyid Al QuwaysiTandalti D Sennar Saraf `Umrah i

a UNNAMED t Nyala r G

b Ed Dueim y Abu Jibeha t

i 8 t 1.5 i r i Shendi t

C 0.0 Al Hilaliyah

Khartoum C s i d

2.0 7 Complementary cumulative 0.5 UNNAMED

5.5 6.0 6.5 5.5 6.0 6.5 5.5 6.0 6.5 City population, log s City population, log s City population, log s 1.5 GuinguinéoGoudomp 10.0 Dakar 0.0 Dagana Tivaouane Senegal (SEN) Gossas Bignona Senegal (SEN) Senegal (SEN) Richard-Toll Kedougou Madina Gounass Sédhiou = 0.95±0.18 Fatick Bambey = 0.95±0.18 = 1.14±0.11 ) s = 55105 Vélingara Dahra Djoloff s = 55105 9.5 s min Koungheul Louga min

( Touba Saint-Louis Kaffrine y r F 1.0 Ziguinchor Kolda MBour g 0.5 g g 9.0 Kaolack Thiès Tambacounda Tambacounda o o Diourbel o Diourbel l l

l Ziguinchor

Richard-Toll Saint-LouisLouga Kolda

, ,

, Kaolack Dahra Djoloff k Thiès P n 8.5 Bambey n o D

i Madina Gounass a Tivaouane t 0.5 MBour r Vélingara G

1.0 u

Fatick Bignona y Touba Gossas b

y Sédhiou

t 8.0 Kedougou Guinguinéo i t i r i Kaffrine Goudomp t C C s i

d 7.5 1.5 0.0 Dakar Koungheul

Complementary cumulative 7.0 Dagana

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

48 Mogadishu 0.00 Somalia (SOM) Berbera Garoowe Somalia (SOM) Somalia (SOM) 1.25 Dara SalaamDhuusamareeb = 1.15±0.26 Las Anod Beledweyne = 1.15±0.26 Qardho City Erigavo = 1.65±0.24 Hargeisa ) 0.25 smin = 72726 Merca Buurhakaba smin = 72726 s 8 Burao

( 1.00 Galkayo Xawo Tako Baidoa y r

F Galkayo Bosaso Kismayo g g Bosaso Kismayo g 0.50 Qardho City Merca o o

o 0.75 l l Beledweyne l Borama Borama

, , Burao Las Anod , 7 Xawo Tako

k Berbera Buurhakaba P n Erigavo

0.75 n Garoowe o

0.50 D Dhuusamareeb i Baidoa a t r G

u Hargeisa y b y t

i 1.00 t i r 0.25 i 6 t C C s i

d 1.25 0.00 Mogadishu 5 Complementary cumulative 1.50 0.25 Dara Salaam 5.00 5.25 5.50 5.75 6.00 5.00 5.25 5.50 5.75 6.00 5.00 5.25 5.50 5.75 6.00 City population, log s City population, log s City population, log s

0.8 10.5 Belgrade 0.0 Serbia (SRB) Subotica Serbia (SRB) Serbia (SRB) = 1.09±0.44 Kragujevac = 1.09±0.44 = 1.27±0.33

) 0.6 0.2 97325 97325 s smin = smin =

( Ni y r F 0.4 10.0 g g g 0.4 Novi Sad o o o l l l

, , Ni , 0.2 k P n 0.6 n o D i

a 9.5 t r G

u 0.0 Belgrade Novi Sad

0.8 y b y t i t i Subotica r i t C C s i 1.0 0.2 d 9.0

0.4

Complementary cumulative 1.2 Kragujevac

5.00 5.25 5.50 5.75 6.00 5.00 5.25 5.50 5.75 6.00 5.00 5.25 5.50 5.75 6.00 City population, log s City population, log s City population, log s

0.0 Aweil Juba South Sudan (SSD) Malakal Yambio South Sudan (SSD) 1.0 Torit Renk = 2.48±0.69 Rumbek = 2.48±0.69

) 9.0 Wau 0.2 smin = 153875 smin = 153875 s Kuajok Malakal ( Wau Yambio Bor 0.8 Yei y r F Rumbek Nimule g g Bor g 0.4 o o Yei 8.5 o l l

l Renk

0.6

, , Aweil , k P n 0.6 Torit n o D i a

t 0.4

r 8.0 G

u Nimule

0.8 y b y t i t i r i

t 0.2 C C s i 1.0 7.5 d

0.0 Juba South Sudan (SSD)

Complementary cumulative 1.2 Kuajok = 1.70±0.66 0.2 7.0 5.2 5.3 5.4 5.5 5.6 5.7 5.2 5.3 5.4 5.5 5.6 5.7 5.2 5.3 5.4 5.5 5.6 5.7 City population, log s City population, log s City population, log s 1.00 Stockholm 0.00 Sweden (SWE) Linköping Västerås Sweden (SWE) Sweden (SWE) Örebro 11.00 = 1.23±0.39 Lund = 1.23±0.39 = 1.10±0.02 ) s = 128003 0.75 s = 128003 s 0.25 min Uppsala min

( Helsingborg y r F 10.75 g g g 0.50 Malmö

0.50 o o o l l

l Gothenburg

, , Gothenburg , 10.50 k P n

n 0.25

o 0.75 D i a t r G

u Malmö

y 10.25 b y t i t 1.00 i 0.00 Stockholm r i t C C s i Helsingborg

d 10.00 1.25 0.25 UppsalaLund

Örebro

Complementary cumulative Linköping 1.50 0.50 9.75 Västerås 5.25 5.50 5.75 6.00 6.25 5.25 5.50 5.75 6.00 6.25 5.25 5.50 5.75 6.00 6.25 City population, log s City population, log s City population, log s

Abu Kamal Damascus 0.0 Kobani`Uthman Syria (SYR) As Safirah Al-Tabqah Syria (SYR) Syria (SYR) Baniyas As Suwayda 10.0 Aleppo = 0.92±0.19 Jablah Idlib = 0.92±0.19 = 1.04±0.03 Bza`a Al Qurayya ) smin = 57365 smin = 57365 s Tartus Manbij ( 1.0 Deir Ez Zor y r Ar Raqqah F Homs g g Qamishli

g 0.5 Al-Hasakah 9.5 o o o l l l

Hama

, , Latakia , k P

n Latakia

n Hama o

0.5 D i Homs a t

r 9.0 G

u Ar Raqqah 1.0 Aleppo Qamishli y b y t

i Tartus Al-Hasakah t i r i Deir Ez Zor t C C s Bza`a

i Manbij As SuwaydaAl Qurayya d Jablah

0.0 8.5 Idlib Damascus Baniyas`Uthman Al-Tabqah

1.5 As Safirah Kobani Complementary cumulative Abu Kamal

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

49 0.0 Mao Doba 9.5 N'Djamena Chad (TCD) Laï Umm Hajar Chad (TCD) Chad (TCD) Bongor Mongo = 2.36±0.59 1.0 Ati Koumra = 2.36±0.59 = 1.39±0.28 ) smin = 184757 Fianga Adré smin = 184757 s

( Sarh y r F 0.5 Pala 9.0 g g g 0.5 Moundou Sarh o o

o Abéché l l l

Abéché Moundou

, ,

, Doba Pala k P n 1.0 Bongor Adré n Ati

o Laï D i a

t 8.5

r Mao

0.0 G

N'Djamena u

y Koumra b y t i t i Umm Hajar r i t

1.5 C C s i

d 8.0 Mongo

0.5

Complementary cumulative 2.0 Fianga

5.4 5.6 5.8 6.0 5.4 5.6 5.8 6.0 5.4 5.6 5.8 6.0 City population, log s City population, log s City population, log s Lomé 0.0 Notsé 9.5 Togo (TGO) 1.0 Tchamba Tsévié Togo (TGO) Togo (TGO) = 1.46±0.40 Blitta-Gare Anié = 1.46±0.40 Kpalimé Sansanné-Mango = 1.16±0.21 ) smin = 76038 smin = 76038 s Dapaong ( Atakpamé 9.0 y 0.5 r F 0.5 Sokodé g g g Kara o o o l l l

, , Dapaong , 8.5 Sokodé k P n 1.0 n o

0.0 D Kara i Lomé Atakpamé a t

r Kpalimé

G Tchamba

y b y t i t i 8.0 r i Notsé t 1.5 C C s 0.5 i Anié

d Tsévié Blitta-Gare 7.5

Complementary cumulative 2.0 1.0 Sansanné-Mango

5.0 5.5 6.0 5.0 5.5 6.0 5.0 5.5 6.0 City population, log s City population, log s City population, log s

Ubon Ratchathani Phra Nakhon Si Ayutthaya 12.0 Lampang Udon Thani Bangkok 0.0 Lopburi Thailand (THA) Surat Thani Nakhon Si Thammarat Thailand (THA) Thailand (THA) SongkhlaPhitsanulok Laem Chabang Khon Kaen = 1.48±0.32 1.0 Pattaya Nakhon Ratchasima = 1.48±0.32 = 1.30±0.20 ) Rayong 11.5 smin = 168877 Nakhon Pathom smin = 168877 s 0.5

( Hat Yai Phuket y r F 0.5 Chonburi g g

g Chiang Mai 11.0 o 1.0 o o l l l

, , , Rayong k 0.0 Bangkok P n 10.5 Chonburi n

o 1.5 D i a t r Laem ChabangPattaya Phuket G

u 0.5

y Phra Nakhon Si Ayutthaya b

y 10.0 t i t i r 2.0 i

t Nakhon Pathom Chiang Mai

C Hat Yai C s i 1.0 Songkhla d 9.5 Surat Thani

Khon Kaen 2.5 PhitsanulokNakhon Ratchasima LampangLopburi Nakhon Si Thammarat

Complementary cumulative 1.5 Udon Thani 9.0 Ubon Ratchathani 3.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s Dushanbe 0.0 Tajikistan (TJK) Qabodiyon Tajikistan (TJK) Tajikistan (TJK) 1.0 Khistevarz Isfara 9.5

= 1.42±0.43 Proletar UNNAMED = 1.42±0.43 = 1.01±0.10 ) smin = 62746 smin = 62746 s Vahdat Hisor ( y r F 0.5 0.5 g g Kulob g o o 9.0

o Khujand l l l

Khujand

, , , k

P Vahdat n n o D i

1.0 a t Kulob r 0.0 Dushanbe G

u Khistevarz Hisor y b y t

i 8.5 t i Proletar r i t C C s

i Isfara

d 1.5

0.5

Qabodiyon

Complementary cumulative 8.0 UNNAMED

4.75 5.00 5.25 5.50 5.75 6.00 4.75 5.00 5.25 5.50 5.75 6.00 4.75 5.00 5.25 5.50 5.75 6.00 City population, log s City population, log s City population, log s 0.75 0.00 Garada aýak Turkmenistan (TKM) Turkmenabat Turkmenistan (TKM) Turkmenistan (TKM) = 1.87±0.76 = 1.87±0.76 10.2 = 0.83±0.36 Ashgabat ) 0.25 s = 130106 0.50 Dashoguz s = 130106

s min min ( Mary y r F g g 0.25 g 0.50 o o Mary o 10.0 l l l

, , , k P n 0.75 0.00 Ashgabat n o D i a t r G

u 9.8 Turkmenabat 1.00 y 0.25

b Dashoguz y t i t i r i t C C s i 1.25 0.50 d 9.6

Complementary cumulative 1.50 0.75 Garada aýak

5.2 5.4 5.6 5.8 6.0 5.2 5.4 5.6 5.8 5.2 5.4 5.6 5.8 City population, log s City population, log s City population, log s

50 10.50 0.00 1.00 Tunis Tunisia (TUN) Medenine Qafsa Tunisia (TUN) Tunisia (TUN)

= 0.98±0.31 Ksar Hellal = 0.98±0.31 = 1.01±0.02

) Gabes 10.25

s 0.25 smin = 95794 0.75 smin = 95794 ( Bizerte y r

F Hammamet g g 10.00 g o o

o 0.50 0.50 l l Sousse l

Sousse , , , Sfax k

Sfax P n 9.75 n o D i

0.75 a 0.25 t r G

y b

y 9.50 t i t i Hammamet r i Bizerte

t 1.00 0.00 Tunis C C s i Ksar Hellal Gabes d 9.25

1.25 0.25 Qafsa Complementary cumulative 9.00 Medenine

5.0 5.5 6.0 5.0 5.5 6.0 5.0 5.5 6.0 City population, log s City population, log s City population, log s

2.0 Istanbul BingölElbistan EregliCizre SamandaNi de ManavgatViran ehir TekirdaKadirli 0.0 KastamonuMu CeyhanÇanakkale Keskin Alanya Turkey (TUR) ÇorluSilopi Erci Turkey (TUR) Turkey (TUR) Silivri KaramanRize ABolur Ere li Yumurtal k ErzurumTokat 11.5 FethiyeSiirt Yalova Turgutlu GiresunSiverek AksarayU ak Karabük Kütahya Zonguldak = 1.06±0.12 Arnavutköy negöl = 1.06±0.12 Günyurdu Çorum = 1.17±0.06 Afyonkarahisar Ordu ) Tarsus Düzce Adiyaman Aydin Osmaniye Bal kesir Isparta s smin = 105953 skenderun smin = 105953 1.5 Kapakl Sivas Manisa Batman Izmir ( Elazig Trabzon 11.0 0.5 Sakarya [Adapazar ]

Antakya y r Van Ankara F anl urfaMalatya Kahramanmara Eskisehir Bursa g g Denizli Samsun g Diyarbakir Adana Sincan o o

o zmit Gaziantep l l Kayseri 10.5 zmit l Antalya

Mersin 1.0 Mersin Konya , , Kayseri , 1.0 k Antalya Sincan P n DenizliEskisehirAdanaKonya n

o Diyarbakir Gaziantep D 10.0 SakaryaDüzce [Adapazar ] Samsun i Bursa a ManisaElazig t KapaklKahramanmara Antakya r negöl Malatya G anl urfa u Aydin skenderun Günyurdu Bal kesir Trabzon 0.5 Bolu Yalova Izmir Tarsus y TurgutluFethiyeArnavutköyOsmaniye SivasVan b SilivriAksaray y

t ZonguldakÇorumIsparta i 1.5 Ankara Batman t i Çorlu U ak Karabük Ni de Kütahya r i 9.5 Afyonkarahisar Manavgat Çanakkale TekirdaRizeEre li t

C KaramanGiresun Kastamonu Kadirli C Ordu s Keskin Siverek Eregli SamandaCeyhan i Tokat YumurtalAlanyak Erzurum Elbistan

d Viran ehir Adiyaman

0.0 ErciSiirt Istanbul 9.0 A r 2.0 Silopi

Cizre Mu Complementary cumulative 8.5 Bingöl 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s New Taipei [Taipei] 0.0 Taiwan (TWN) 1.0 Xinying District Taiwan (TWN) Taiwan (TWN) Miaoli 11.5 = 0.50±0.15 Hualien = 0.50±0.15 = 1.00±0.01 ) smin = 85307 Luodong smin = 85307 s 0.2 Chiayi Taichung ( 0.8 Kaohsiung y r F g g Hsinchu g 11.0 o o o l l Tainan l 0.4

, , 0.6 Tainan , Hsinchu k P n n o

Kaohsiung D i a t

0.6 r 10.5 Chiayi G

0.4 y b y t i Taichung t i r i t C C s Luodong i 0.8 0.2 d 10.0 Miaoli Hualien

Complementary cumulative 1.0 0.0 New Taipei [Taipei] Xinying District 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

Singida Dar es Salaam 1.5 Ifakara Bukala 10.5 0.0 Iwambi Mpanda Tanzania (TZA) Mtwara Murubona Tanzania (TZA) Tanzania (TZA) Tunduma Ndiuka KatoroGeita = 1.11±0.21 Kahama Shinyanga = 1.11±0.21 = 1.17±0.07 Bukoba Sumbawanga ) smin = 79532 Musoma Tabora smin = 79532 s Songea Moshi 10.0 ( 1.0 Tanga Dodoma y r

F 0.5 Kilima Hewa Morogoro g g g o o

o Mbeya l l Arusha l 9.5

, , , Zanzibar City k P n 1.0 0.5 Zanzibar City ArushaMwanza

n Mbeya o

Mwanza D i a t

r 9.0 G

u Morogoro

y b

y Songea t i t i Kilima Hewa r i TangaMoshi t 0.0 Dar es Salaam Katoro 1.5 C Sumbawanga Dodoma C

s Tunduma MusomaNdiuka i 8.5 Tabora Mpanda IwambiKahamaGeita d Mtwara Bukoba Bukala

Singida MurubonaShinyanga

Complementary cumulative 2.0 0.5 8.0 Ifakara

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Busesa 0.0 Hoima Opuyo Kampala Uganda (UGA) Mubende Fort Portal Uganda (UGA) Uganda (UGA) BusiaLugazi = 1.34±0.31 Kasese Koboko = 1.34±0.31 = 1.10±0.09 1.0 Mbale Masaka ) smin = 63987 smin = 63987 9.5

s Lira Mbarara

( 0.5 Gulu y r Iganga F g g g 0.5 Arua o o o l l Jinja 9.0 l

, 1.0 , , k P n n

o GuluArua D Jinja i 0.0 a Kampala Iganga

t Lira r

G 8.5

1.5 y b y t

i Mbale t i Opuyo r i Kasese Koboko t Masaka C

0.5 C s Mubende Mbarara

i Fort Portal 8.0 Busia d Lugazi

2.0 Hoima

Complementary cumulative 1.0 Busesa 2.5 7.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

51 Ternopil Kyiv 1.5 Horlivka Lutsk 0.0 Kremenchuk Kamianske Ukraine (UKR) Khmelnytskyi Sumy Ukraine (UKR) Ukraine (UKR) SievierodonetskChernihiv Cherkasy Zhytomyr = 1.59±0.29 Sevastopol Rivne = 1.59±0.29 = 1.40±0.10 Poltava Kherson ) smin = 256042 Ivano-Frankivsk Chernivtsi smin = 256042 s Vinnytsia Simferopol 10.5 ( Mariupol Luhansk y r F 0.5 1.0 Mykolaiv Kryvyi Rih g g g Donetsk o o Zaporizhzhia Lviv o l l

l Dnipro Kharkiv

, ,

, Dnipro

k Odesa 10.0 Odesa P n

n Kryvyi Rih

o Zaporizhzhia D i a

t 1.0 0.5 Donetsk r G

u Mariupol Kharkiv Lviv y b y t i t i Mykolaiv r i Poltava Luhansk

t Kamianske

C 9.5 Horlivka C

s Kremenchuk Simferopol i Vinnytsia Sevastopol d SievierodonetskCherkasy 1.5 0.0 Kyiv Ivano-Frankivsk RivneKherson SumyChernihivZhytomyr KhmelnytskyiLutsk Chernivtsi Complementary cumulative 9.0 Ternopil

5.50 5.75 6.00 6.25 6.50 5.50 5.75 6.00 6.25 6.50 5.50 5.75 6.00 6.25 6.50 City population, log s City population, log s City population, log s

Lowell New York SantaNuevo Clarita Laredo Hemet AnchorageColumbusRacine PortlandSaginawAppleton DaytonaLubbock Beach TallahasseeEvansville WestMontgomeryEugene Chester ManchesterBrownsvilleLincoln VisaliaKalamazoo 0.0 FraminghamAntioch PortHesperia Arthur Ann ArborPeoria KilleenSavannah United States (USA) SantaShreveport Rosa YorkNaples United States (USA) United States (USA) Davenport Springfield CorpusWaterbury Christi Huntsville BrocktonPalmdale LawrencePensacola FayettevilleStockton RockfordMobile LafayetteOxnard SalemLancaster ModestoLansing BoiseReno 12.0 TemeculaFlint PortFort St. Wayne Lucie Winston-SalemCanton LittleYoungstown Rock ReadingIndio LakelandWichita 2.0 WorcesterSpokane SouthDes Moines Bend Los Angeles = 0.92±0.07 LexingtonMadison KnoxvilleWaukegan = 0.92±0.07 = 1.02±0.01 Syracuse Ogden PalmProvo BayColumbiaToledo CharlestonAkron ColoradoBakersfield Springs ) HarrisburgCape Coral GreensboroGrand Rapids SpringfieldAlbany Tulsa s = 261876 Baton Rouge RochesterBradenton s = 261876 min BirminghamFresno Concord min s McAllen Albuquerque Chicago Omaha NewDayton Orleans Ciudad HonoluluJuárez Allentown ( Tucson Richmond Hartford Bridgeport Buffalo Jacksonville Nashville Oklahoma City y r Salt Lake City Louisville F Memphis Raleigh 0.5 Providence Virginia Beach HoustonWashingtonPhiladelphia D.C. Milwaukee Columbus Dallas ClevelandIndianapolis Miami g g Austin Kansas City 11.5 g 1.5 Cincinnati Charlotte AtlantaSan Jose Pittsburgh San Antonio Sacramento

o Seattle o Baltimore Phoenix

o Las Vegas l l Portland Orlando l BostonTijuana

Tampa St. Louis Detroit Minneapolis [Saint Paul] , , Boston Denver , Denver Tijuana Minneapolis [Saint Paul] Baltimore k

P Tampa

n St. Louis Seattle Detroit Sacramento LasPortland Vegas Orlando n Phoenix San Jose Pittsburgh

o San Antonio

D Austin i 1.0 ClevelandCincinnati Kansas City a 1.0 Providence Virginia Beach AtlantaWashington D.C. ColumbusCharlotte t 11.0 Milwaukee

r Hartford Indianapolis Bridgeport G Buffalo Raleigh u Miami Philadelphia Richmond y b

y Memphis t Honolulu AllentownJacksonvilleLouisville i SaltNashvilleRochester Lake City Oklahoma City t i Houston Springfield FresnoOmahaTucsonConcord r Dallas i Albany McAllen Dayton Ciudad JuárezBradentonNew Orleans t

C Birmingham Albuquerque SyracuseBaton Rouge Tulsa C WorcesterHarrisburg Bakersfield s Cape CoralColorado Springs

i 0.5 PalmAkron Bay GreensboroGrand Rapids Chicago Waukegan SpokaneIndioMadison Toledo 1.5 CharlestonWichita 10.5 Reading Des Moines WaterburyBrockton d ModestoTemeculaLawrence LakelandKnoxvilleOgdenColumbia Los Angeles LexingtonProvoYoungstownSouth Bend Oxnard Stockton CantonLancasterPort St. Lucie Palmdale Rockford Framingham Winston-SalemRenoFort Wayne FlintSalem Little Rock LafayetteHesperiaYorkLansing AntiochSantaDavenportPeoria Rosa Pensacola CorpusFayetteville ChristiVisaliaNaplesMobileBoise ManchesterShreveport Springfield Lowell PortKilleen ArthurHuntsville LincolnAnchorageWest ChesterSavannah BrownsvilleHemetAnn Arbor TallahasseeAppleton RacineEugeneKalamazoo Complementary cumulative SantaNuevoDaytonaLubbock Clarita MontgomeryLaredo Beach PortlandSaginawEvansville 0.0 Columbus New York 10.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

Yozyovon Tashkent Poytugh Khujayli Sherobod 0.0 Gallaorol Buloqboshi 10.5 Uzbekistan (UZB) Chust Jarqurgon Uzbekistan (UZB) Uzbekistan (UZB) Pakhtaobod Urgut Xonqa Yaypan Koson Guzar Gulistan 1.5 Uzun Chirchiq = 1.00±0.15 Beruniy Hazorasp = 1.00±0.15 = 1.14±0.08 AltiarikShahrixon ) Rishtan Kosonsoy smin = 70003 Yangikurgan smin = 70003

s Navoiy Quva Khiva ( Bukhara Denov y r 10.0 F Gijduvan Termez 0.5 Katta-Kurgan

Urgench g g g Nukus Turtkul Shahrisabz o o Qarshi o l l Samarkand

l 1.0 Namangan

Shahrisabz Kokand Gijduvan Fergana

, , Jizzax , Bukhara Andijan

k Jizzax Termez

P 9.5 n Andijan n o D

i Navoiy Urgench a Fergana t Kokand r Qarshi Denov

1.0 G

u Uzun Turtkul y b

y Quva t Katta-Kurgan i t i 0.5 Namangan Hazorasp Nukus r i 9.0 t Gulistan Yangikurgan

C Chirchiq Beruniy Shahrixon C Khiva s Samarkand Rishtan i Koson

d Yozyovon Altiarik Kosonsoy

Guzar XonqaYaypan 1.5 Khujayli Urgut 8.5 Chust Buloqboshi Poytugh Gallaorol Complementary cumulative Jarqurgon 0.0 Tashkent PakhtaobodSherobod

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s Caracas Biruaca 11.0 0.00 Ciudad Ojeda Venezuela (VEN) Valera Coro Venezuela (VEN) Venezuela (VEN) Punto Fijo Santa Teresa del Tuy Maracaibo AcariguaCabimas = 1.31±0.26 1.25 Charallave Mérida = 1.31±0.26 = 1.22±0.27

) Catia La Mar 0.25 smin = 241935 Barinas smin = 241935 s Cumana Ciudad Bolívar 10.5 ( Porlamar Guarenas Ciudad Guayana y r

F Maracay 1.00 Barcelona MaturinSan Cristóbal g g g 0.50 Punto Fijo Coro CabimasMérida Maturin Barquisimeto o o 10.0 Ciudad Bolívar o Ciudad Guayana l l Barcelona Ciudad Ojeda l

, ,

, 0.75

k Barquisimeto P n 0.75 Maracay BarinasPorlamar n San Cristóbal Valencia o

D Catia La Mar i AcariguaCumana

a 9.5 Guarenas t

r Charallave G

u 0.50

Valencia Santa Teresa del Tuy y b y

1.00 t i t i r Maracaibo i t

C 9.0 C s 0.25 i Biruaca

d 1.25

0.00 Caracas 8.5 Complementary cumulative 1.50 Valera

5.50 5.75 6.00 6.25 6.50 5.50 5.75 6.00 6.25 6.50 5.50 5.75 6.00 6.25 6.50 City population, log s City population, log s City population, log s

Trà Vinh City Ho Chi Minh City Hi p Ph c Cao Lanh Ph Lý 0.0 Tuy Hòa Ninh Bình Vietnam (VNM) Dalat Ca Mau Vietnam (VNM) Vietnam (VNM) Sóc Tr ng i Ngh a 1.5 Pleiku Phan Rang - Tháp Chàm 11.0 B c Liêu Buôn Ma Thu t Phan Thiet Qu ng Ngãi = 1.14±0.18 R ch Giá = 1.14±0.18 Cai L y Ha Long = 1.11±0.05 H i D ng

) Phúc Th smin = 188898 Quy Nh n Tây Ninh smin = 188898 Hanoi s Thái Nguyên Vinh

( Hue Thanh Hóa y 0.5 r Vung Tau F Nha Trang 10.5 1.0 Sa Dec City My Tho g g g Can Tho o o à N ng o l l l

, , Nam nh , Haiphong k P n 1.0 Long Xuyen

n 10.0 Haiphong o

0.5 Long Xuyen D Nam nh i a t

r Hanoi Can Tho G

u Vung Tau à N ng y b y t Nha Trang i t i r i Sa Dec City My Tho t 1.5 Phúc Th Ha Long C 9.5 C s 0.0 R ch Giá Hue Thanh Hóa i Ho Chi Minh City Cai L y Tây Ninh Sóc Tr Hing p Ph c d QuyThái Nh Nguyênn Vinh i Ngh a H i D ng Ca Mau Ph Lý Buôn Ma Thu t B c LiêuNinh Bình Qu ng Ngãi Cao Lanh DalatPleiku Phan Rang - Tháp Chàm Complementary cumulative 2.0 TràPhan Vinh ThietCity 0.5 9.0 Tuy Hòa

5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

52 Ad Dala Sana'a 0.00 Al Qatn Hajja 10.0 Yemen (YEM) Al Milhi Rada'a Yemen (YEM) Yemen (YEM) 1.25 Al Beidha Shihr = 0.93±0.19 Lawdar Harad = 0.93±0.19 = 1.23±0.17 ) 0.25 smin = 59316 Seiyun Tarim smin = 59316 9.5 Al Mansurah Ta'izz s AmranAbs ( 1.00 Hodeidah y r Sa'dah F Al Mutadarririn Ibb g g

g Yarim 0.50 Dhamar 9.0 o o o l l Dhamar l

0.75

, , Ibb Al Mutadarririn , Hodeidah Seiyun Tarim Abs k P n 8.5 0.75 Al Qatn n

o Harad

D Rada'a i Shihr a

t 0.50

r Al Mansurah Hajja G

Amran u Al Milhi

y Al Beidha b

Ta'izz y 8.0 1.00 t i t i r i Sa'dahYarim t 0.25 C C s i 7.5 Ad Dala d 1.25

0.00 Sana'a

Complementary cumulative 1.50 7.0 Lawdar 0.25 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s

Kroonstad Johannesburg Ermelo Naas Uitkomsdal Mankweng 0.0 Stellenbosch Sikhalasenkosi South Africa (ZAF) Worcester Mokopane South Africa (ZAF) South Africa (ZAF) Mpumalanga Sasolburg 11.0 Salubindza eMbalenhle Vosman 1.5 Jouberton Potchefstroom = 0.92±0.14 KwaMhlanga Bronville = 0.92±0.14 = 1.02±0.02 George Mafikeng Middelburg ) Botshabelo smin = 87465 Emalahleni Somerset West smin = 87465 s Kimberley Paarl Cape Town ( eGobhoza Phuthaditjhaba Durban y r VanderbijlparkOsizweni F 0.5 Polokwane Rustenburg 10.5 g g g East LondonBloemfontein Pretoria Klipgat o o

o 1.0

Edendale l l Tsakane l

, , , Evaton Port Elizabeth Evaton Port Elizabeth k P n

n 10.0

o Klipgat D i Pretoria Tsakane a

t 1.0 r East London G

u Edendale

Bloemfontein Vanderbijlpark y 0.5 Polokwane

b Durban y t Rustenburg i t i eGobhoza r Cape Town i t Emalahleni C 9.5 PhuthaditjhabaKimberley

C Osizweni s MafikengMiddelburg PotchefstroomPaarl i KwaMhlanga Somerset West Botshabelo eMbalenhle George

d Jouberton Vosman 1.5 Salubindza Bronville Mokopane Uitkomsdal StellenboschErmelo MankwengSasolburg 0.0 MpumalangaNaas Johannesburg Sikhalasenkosi

Complementary cumulative 9.0 Kroonstad Worcester

5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 5.0 5.5 6.0 6.5 7.0 City population, log s City population, log s City population, log s

Chililabombwe Chinsali Mkushi Mazabuka Lusaka 0.0 1.5 Isoka Samfya Zambia (ZMB) UNNAMED Mumbwa Zambia (ZMB) 10 Zambia (ZMB) Mwinilunga Kalomo Mpulungu Luanshya Kafue Lundazi Monze = 1.54±0.27 Nchelenge Mongu = 1.54±0.27 = 1.25±0.10 SolweziMbala ) smin = 77386 TundumaMufulira smin = 77386 Ndola s Choma Livingstone Kitwe ( 0.5 1.0 Kapiri Mposhi Chipata KabweKasama y r 9 Kapiri Mposhi F Chipata Mansa Kabwe Chingola Mbala LivingstoneChingola Lundazi ChomaMufulira g g Tunduma

g UNNAMED MwinilungaMpulunguLuanshyaSolwezi Kasama NchelengeKalomoMonze o o Mansa Mumbwa o ChililabombweKafue l l

l Mongu

Isoka

, ,

, 1.0 0.5 MazabukaSamfya Ndola 8 Chinsali k P n Kitwe n o D i a t r G

b 1.5 0.0 y t Lusaka i 7 t i r i t C C s i d 2.0 0.5 6 Complementary cumulative Mkushi

5.0 5.5 6.0 5.0 5.5 6.0 5.0 5.5 6.0 City population, log s City population, log s City population, log s Harare 0.0 MaronderaChegutu Zimbabwe (ZWE) Triangle Chinhoyi Zimbabwe (ZWE) Zimbabwe (ZWE) = 1.07±0.27 Norton Kadoma = 1.07±0.27 9.5 1.0 King MineMasvingo = 1.03±0.07 ) smin = 54609 smin = 54609 s Bindura Glenclova ( y r Bulawayo F Mutare 0.5 Kwekwe 9.0 g g g o o o l l 0.5

l Gweru

, , , Bulawayo Gweru k P n 8.5 Mutare n o D i

1.0 a

t Kwekwe r Masvingo G Triangle Bindura u Norton

y b

0.0 y t Harare

i Kadoma t i Marondera r i 8.0 Chegutu t

C Chinhoyi C s i

d 1.5

King Mine 0.5 7.5 Complementary cumulative Glenclova

5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 5.0 5.5 6.0 6.5 City population, log s City population, log s City population, log s