Research Article  Open Access
A Hierarchical Allometric Scaling Analysis of Chinese Cities: 1991–2014
Abstract
The law of allometric scaling based on Zipf distributions can be employed to research hierarchies of cities in a geographical region. However, the allometric patterns are easily influenced by random disturbance from the noises in observational data. In theory, both the allometric growth law and Zipf’s law are related to the hierarchical scaling laws associated with fractal structure. In this paper, the scaling laws of hierarchies with cascade structure are used to study Chinese cities, and the method of analysis is applied to analyzing the change trend of the allometric scaling exponents. The results show that the hierarchical scaling relations of Chinese cities became clearer and clearer from 1991 to 2014 year; the global allometric scaling exponent values fluctuated around 0.85, and the local scaling exponent approached 0.85. The Hurst exponent of the allometric parameter change is greater than 0.5, indicating persistence and a longterm memory of urban evolution. The main conclusions can be reached as follows: the allometric scaling law of cities represents an evolutionary order rather than an invariable rule, which emerges from selforganized process of urbanization, and the ideas from allometry and fractals can be combined to optimize spatial and hierarchical structure of urban systems in future city planning.
1. Introduction
Cities as systems (individuals) and systems of cities (groups) are scalefree complex systems, which cannot be effectively described by the traditional mathematical methods based on characteristic scales in many respects. The ideas from scaling can be used to research urban systems (see, e.g., [1–7]). Two correlated scaling laws are often employed to analyze a hierarchy of cities: one is Zipf’s law, and the other is the law of allometric growth. Zipf’s law, allometric growth law, and distancedecay law compose three basic laws of urban geography. Zipf’s law indicates the ranksize pattern of cities in a geographical region [8–12], and the allometric growth law describes the relationship between size and shape in the growth of human settlements [13–17]. In fact, Zipf’s law reflects urban growth [1, 9], and the allometric scaling law can be derived from dual Zipf’s models of ranksize distributions of urban population and area [18]. This suggests that Zipf’s law and the allometric scaling law represent different sides of the same coin. Zipf’s law proved to be equivalent to a hierarchical scaling law, and the ranksize allometric scaling can be replaced by the hierarchical allometric scaling [19]. Hierarchy mirrors a universal structure in natural and social systems [20]. Zipf’s law is a signature of hierarchical structure. Based on hierarchical scaling, we can develop a new approach to studying urban systems.
China bears a large set of cities with a long history. Studies on Chinese cities will help us understand the hidden order of complex systems. There is a dispute about whether the size distribution of Chinese cities follows Zipf’s law [22–27]. In fact, Zipf’s law is a rule of evolution rather than that of existence. The ranksize pattern emerges from the edge of chaos [28]. A new discovery is that Chinese citysize distribution can be described by the threeparameter Zipf’s model instead of the twoparameter Zipf’s model [29]. The threeparameter Zipf’s law suggests an incomplete hierarchy with cascade structure. This implies that the hierarchical scaling law can be used to research the allometry and ranksize pattern of Chinese cities.
This paper is devoted to making a hierarchical allometric analysis of Chinese systems of cities. The aim of this study is as follows. First, we try to reveal the evolutional process and characteristics of the allometric scaling in Chinese cities. Second, we attempt to bring to light the causality behind the hierarchical allometry of Chinese cities. Third, we will sum up a general framework of hierarchical allometric analysis of cities. The trait of this case study rests with large samples, continuous time series, and new angle of view. By this work, we can obtain useful geographical information about spatiotemporal evolution of Chinese cities and new knowledge about scaling in cities. Moreover, the study lends further support to the suggestions that the geographical laws are evolutional laws and there are inherent relationships between Zipf’s law and the law of allometric growth of cities. The remaining parts of the paper are organized as below. In Section 2, the basic mathematical models of hierarchical structure are presented and explained; in Section 3, empirical analyses of hierarchical allometric scaling in Chinese cities are made by means of two algorithms, to show the evolutional regularity of Chinese cities; in Section 4, several questions are discussed, and a general process of allometric analysis is proposed for urban studies. Finally, in Section 5, the article will be concluded by summarizing the mains of this study.
2. Models
2.1. Hierarchical Scaling Law of Cascade Structure
The mathematical models of selfsimilar hierarchies of cities can be expressed as a set of exponential functions and power functions. Using these models, we can make allometric scaling analysis based on citysize distribution. Suppose that the cities in a geographical region are grouped into M classes in a topdown order according to the generalized 2^{n} principle [19, 30, 31]. The cascade structure of the urban system can be modeled by three exponential equations as follows:where m refers to the topdown ordinal number of city level (), denotes the number of cities of order ; correspondingly, and represent the mean population size and urban area at the th level. The meaning of the parameters is as below: refers to the number of the toplevel cities, and are the mean population size and urban area of the toplevel cities, / is the interclass number ratio of cities, is the population size ratio, and is the urban area ratio. Equations (1), (2), and (3) compose the mathematical expressions of the generalized 2^{n} rule [19, 21], which is based on BeckmannDavis models [30, 32]. In theory, if as given, then it will follow that and vice versa. Here the arrow denotes “approach” or “be close to.” If , the generalized 2^{n} rule will return to the normal 2^{n} rule presented [31].
2.2. Hierarchical Allometric Rescaling
The cascade structure of a hierarchy of cities suggests allometry, fractal, and scaling. A set of powerlaw relations including the threeparameter Zipf’s law can be derived from the above exponential laws [21, 29]. The powerlaw models are as below:where , /, , , , and . Equation (4) is termed the hierarchical sizenumber scaling relation of cities, which is equivalent to the Pareto law of citysize distribution, and D is the fractal dimension of urban hierarchies measured with urban population. Equation (5) is termed the hierarchical areanumber scaling relation of cities, which is equivalent to the Pareto law of cityarea distribution, and d can be treated as the fractal dimension of urban hierarchies measured with urban area. According to Chen [19], both D and d are actually paradimension rather than real fractal dimension. Equation (6) is termed the hierarchical allometric scaling relation between urban area and population, and b is the allometric scaling exponent of an urban hierarchy. The inverse functions of (4) and (5) are equivalent to Zipf’s models of urban population size and urban area distributions. The allometric scaling exponent is actually the ratio of the fractal dimension of urban area size distribution to the fractal dimension of urban population size distribution; that is, . Generally speaking, , and is close to 1. Thus, value comes between 0 and 1.
The exponential laws and the power laws reflect two relations of a hierarchy, respectively: longitudinal relations and latitudinal relations. The longitudinal relations are the associations across different classes, while the latitudinal relations are the correspondences between different measures such as city population size and urbanized area [21]. Formally, these hierarchical and measurement relations can be illustrated by dual hierarchies (Figure 1). If cities in a region follow Zipf’s law, they can be organized into a hierarchy with cascade structure. On the other hand, if an urban system bears cascade structure, the cities in the system follow Zipf’s law or can be described with a threeparameter Zipf’s model. The Zipf distribution is a signature of selfsimilar hierarchies associated with fractal patterns and selforganized processes.
3. Empirical Analysis
3.1. Study Area and Methods
The models of hierarchical allometry can be applied to the ranksize distribution of Chinese cities. The study area contains the whole mainland of China, in which there are about 660 officially approved cities. The city number is different in different years. The city size is measured by the urban population within urbanized area, while the urban area is represented by the area of built district. A number of datasets of city sizes and urban area are available for this research, including the statistical data from 1991 to 2014 (24 years). All the observational data of city population size and urban area come from the Ministry of Housing and UrbanRural Development of China (MOHURD).
A number of algorithms can be employed to estimate the power exponent values, including the least squares method (LSM), maximumlikelihood method (MLM), and major axis method (MAM). Because of the advocacy of Newman [33] and Clauset et al. [34], the MLM is treated as the only standard approach to fitting power law to observational data. In fact, the precondition of effectively making use of the MLM is that the observational data meet the joint normal distribution. However, for many social and economic systems such as cities, the observational variables do not satisfy the joint normal distribution. In this case, the LSM has its merits. The analytical process is as follows. Step 1, preliminary analysis: using Zipf’s law and the allometric growth law, we can examine the primary observational datasets based on urban population and area measurement. Step 2, hierarchical reconstruction: if the citysize or urban area distribution follows Zipf’s law and the relationships between urban area and population size follow the allometric scaling law, we can reconstruct the hierarchy by organizing the cities according to the model of cascade structure, (1). Step 3, cascade analysis: using the exponential laws, (2) and (3), we can investigate the cascade structure of system of cities. Step 4, allometric scaling analysis: using the power laws, (4), (5), and (6), we can research the hierarchical allometry of cities by urban size and area. Among various equations, (6) plays the most important role in our case study.
3.2. Results and Findings
Now, let us make an allometric scaling analysis of hierarchy of Chinese cities step by step. First of all, we should examine the ranksize distribution of cities which can be described by Zipf’s formulation. Zipf’s law is universal rule followed by many natural and social systems (see, e.g., [9, 28, 35–42]). This suggests that urban evolution follows the general laws of nature under certain condition. Zipf’s distribution proved to be a signature of hierarchy with cascade structure [21]. If Chinese cities follow Zipf’s law, they can be organized into a selfsimilar hierarchy. The twoparameter Zipf models of urban population and area can be expressed aswhere denotes the rank of a city, and refer to the population size and builtup area of the th city, and are the population size and urban area of the largest city, and and are the Zipf scaling exponents. Using the OLS calculation, we can fit (7) to the observational data of urban population and area of Chinese cities. For example, where population size is concerned, in 2000, there were 665 officially approved cities, but only the largest 250 cities conform to Zipf’s law; in 2010, however, there were 656 officially approved cities, but the largest 550 cities comply with Zipf’s law approximately (Figure 2).
(a) Population size
(b) Builtup area
Then, we investigate the allometric scaling relation based on the ranksize distribution. This is a kind of transversal allometry based on crosssectional data [43, 44]. In fact, from (7) it follows thatin which and , where and [18]. Fitting (8) to the observational data of population size and builtup area of Chinese cities yields models of crosssectional allometry (Figure 3). In 2010, the relationship between urban population and area follow the law of allometric scaling. However, in 2000, the urban population and area take on a linear relationship rather than a powerlaw pattern. In other words, in earlier years, Chinese cities departed from the allometric scaling law to some extent. This indicates that the scaling law of crosssectional allometry is indeed associated with Zipf’s law. Since Chinese cities failed to comply with Zipf’s law in early years, the allometric scaling relation could not form. What is more, this lends further support to the suggestion that the laws of human geography are evolutional laws instead of existent laws [4, 45].
(a) 2000 (665 cities)
(b) 2010 (656 cities)
In the third step, we can check the cascade structure of urban hierarchy. The selfsimilar hierarchy can be constructed as follows. Suppose that we divide all the cities into M classes/levels in a topdown order according to the generalized 2^{n} rule [4, 19, 21, 31]. The result is one city in the first class (rank 1), two cities in the second class (ranks 2 and 3), four cities in the third class (ranks 4, 5, 6, and 7), and so on (Figure 1). This cascade structure is based on city population rather than urban area. Generally speaking, urban area development lags behind population size growth of a city. Then, we can calculate the average population size and average urban area of the cities at different levels (Table 1). If the average population size and average urban area follow the exponential laws, we can describe the size and area decay using (2) and (3). The results show that the average population size and area take on exponent decay patterns (Figure 4). This suggests that the relationships between average urban population and area in different classes follow the power laws and can be described by (4), (5), and (6). The results show that the first and the last classes are sometimes exceptional points, and the other points form a scaling range and comply with the power laws (Figure 5). The scaling exponent of (4) gives the fractal dimension of population size distribution (), while the scaling exponent of (5) yields the fractal dimension of urban area distribution () [19].
 
Note. These results are similar to the results based on urban census data, which were displayed in Chen [29], but there are subtle differences. Unfortunately, we have no datasets for urbanized area of Chinese cities. 
(a) Population size
(b) Urban area
(a) Population size
(b) Urban area
In the fourth step, we should research the hierarchical allometric scaling relationships between urban population and area. This step is the focus of this study. By statistical analysis and double logarithm plots, we can discover a number of new phenomena and trends (Figure 6). First, the allometric scaling relation emerges from the selforganized process of urbanization. From 1991 to 1992 year, the urban populationarea relation follows the powerlaw indicative allometric scaling, but the (local) allometric exponent is greater than 1. The proper allometric scaling exponent should come between 2/3 and 1 (Table 2). From 1993 to 2002 year, the local allometric relation based on scaling range degenerated to a linear relation. In particular, during the period from 1998 to 2000, both the global allometry and local allometry degenerated to linear relations (Table 3). In other words, it is the linear function rather than the power function that can be best fitted to the observational data of 1998, 1999, and 2000. This indicates that the allometric scaling relation between urban population and area has not yet emerged before 2003. The powerlaw relation can be constrainedly fitted to the observational data from 1993 to 2002 year, but the points do not well match the trend lines. In fact, at the early stage, the mediumsized cities with population from 50 to 150 thousand did not develop where urbanized area is concerned. From 2003 on, the allometric scaling relation emerged, and this relation became significant in 2004. Once again, this lends further support to the suggestion that geographical laws are evolutionary laws instead of everlasting laws. What is more, the allometric scaling exponent ranged from 2/3 to 1. However, from 2002 on, the land of the outsize cities with population greater than 400 thousand was overused because of real estate bubble economy. Second, the allometric scaling law depends on Zipf’s law. Generally speaking, population growth precedes the development of urban area indicative of land use, but the ranksize distribution of urban area is more stable than that of population. However, for Chinese cities, the thing is abnormal. For example, from 1998 to 2000, the 500 largest Chinese cities followed Zipf’s law of urban area, but only about 250 largest cities complied with the common Zipf law of city population. In other words, the rankpopulation distribution was not consistent with the rankarea distribution. As a result, the allometric scaling relation of Chinese cities degenerated from powerlaw relation to quasilinear relation. This lends further support to the inference that urban areapopulation allometry proceeds from the scaling relation between Zipf’s distribution of city population and that of urban area [18]. Third, on the whole, the allometric scaling exponent value changed around constant. In theory, the scaling exponent of allometric growth () is the ratio of the fractal dimension of urban landuse form () to that of urban population distribution () [18]. Thus, the allometric exponent is supposed to fluctuate between 2/3 and 1 and approaches 0.85 because and [46]. However, in practice, the computation yields a fraction sequence such as , where refers to the set of natural numbers [47]. Fourth, there is scalefree range in which the allometric relation is more significant. On loglog plots, two data points often take on outliers: one is first class, and the other is the last class. On the one hand, number one is always a special one in theory [19]. The top class indicating the largest city is usually different from other classes and departs from the trend line. On the other, the bottom class also manifests an exceptional value in a plot due to undergrowth of small cities or absence of partial data [29]. The exceptional values are clear in the plots of ranksize distributions and the corresponding sizenumber scaling relation (see Figure 5). The other points coming between the first and last classes form a straight line and represent a scaling range. However, the outliers of ranksize distribution or sizenumber scaling often become insignificant in the hierarchical allometric relation (Figure 6). Despite this, if we fit the model of allometric scaling to the data points within the scaling range, the allometric exponent values will rise, and the mean is about 0.95 (Table 2). Please note that the local allometry is based on the data points within the scaling range in this context. Accordingly, the global allometry is based on all the data points, including the data points inside and outside the scaling range.
 
Note. () The notations are as follows: a—proportionality coefficient, —allometric scaling exponent, — goodness of fit, —size dimension (the fractal dimension of population size distribution), —area dimension (the fractal dimension of urban area distribution), and —fractal dimension quotient or dimension ratio. () Removing the first class and the last class yields the scaling range. 
 
Note. The global analysis is based on all the classes, which can be represented by all the data points on the loglog plots. The local analysis can be made by removing the first class and the last class. For the powerlaw relation, the local analysis is based on the scaling range. 
(a) 1991
(b) 1995
(c) 2000
(d) 2005
(e) 2010
(f) 2014
3.3. Analysis of Allometric Scaling Exponents
The changing trend of the allometric parameter values reflects the trend of urban evolution. As shown above, four sets of parameters have been evaluated. The first is the allometric scaling exponent based on the global analysis, ; the second is the estimated value of by the ratio of the fractal dimension of city population size distribution (size dimension) to that of urban area distribution (area dimension) based on scaling ranges, ; the third is the allometric scaling exponent based on the local analysis, ; the fourth is the estimated value of by the ratio of the size dimension to area dimension based on scaling ranges, . The global allometric exponent b and its estimated result D/d fluctuated around the mean value 0.85, but the local allometric exponent and its estimated result descended and take on a trend of damped vibration—being gradually decayed and approaching 0.85 (the empirical value is about 6/7).
The method of the rescaled range analysis can be employed to predict the future direction of parameter change, which indicates the trend of urban evolution. The rescaled range analysis is also termed R/S analysis, which is an approach of nonlinear time series analysis put forward by Hurst et al. [48]. The basic parameter of the R/S analysis is Hurst exponent, which is associated with the selfaffine fractal dimension and autocorrelation coefficient [49, 50]. The Hurst exponent value comes between 0 and 1. If the Hurst exponent is equal to 1/2, the time series bears no long memory and urban change has no autocorrelation; if the Hurst exponent is greater than 1/2 significantly, the time series bears persistence, and urban change has positive autocorrelation; if the Hurst exponent is less than 1/2 significantly, the time series bears antipersistence and urban change has negative autocorrelation. The calculations indicate that the method is suitable for time series analysis of the allometric scaling exponents (Figure 7). The results show that all the Hurst exponent values are greater than 1/2 (Table 4). This suggests that global allometric exponent will change around the mean value as usual, but the local allometric exponent will continue to decay before it arrives at the expected value, 0.85. An inference is that the global allometric exponent and the local allometric exponent will reach the same goal by different routes.

(a) value
(b) value
(c) value
(d) value
3.4. The Results Based on MaximumLikelihood Fitting
The algorithm adopted in this study is in fact the ordinary least squares (OLS) method, which can be easily applied to double logarithmic linear regression. In fact, OLS is a conventional approach to estimating the power exponents. However, the OLS algorithm is not always the best one for fitting power laws to observational data. A new method based on the maximumlikelihood estimation (MLE) is developed by Clauset et al. [34] to address the powerlaw distribution in empirical data. Combining the methods of maximumlikelihood fitting with the goodnessoffit tests based on the KolmogorovSmirnov statistic and likelihood ratios, the approach is employed by many scholars to identify various powerlaw distributions in both natural and social sciences. Unfortunately, the MLEbased approach cannot be directly applied to our datasets in this work. The reasons are as follows. First, the MLEbased method is developed for powerlaw frequency distributions, while this study is devoted to allometric scaling relations. Second, the MLEbased method is suitable for binned data, while this study is based on pairs of cascade sequences.
However, the MLEbased approach can be indirectly applied to the hierarchical allometric scaling analysis. As indicated above, an allometric scaling relation can be derived from two Zipf’s distributions, and the hierarchical allometry is equivalent to the ranksize allometry [29, 46]. The MLEbased method developed by Clauset and his coworkers can be used to research Zipf’s distributions. If urban population size distribution and area size distributions follow Zipf’s law simultaneously, the relation between urban area and population will follow the allometric scaling law [4]; if the urban areapopulation relation follows the ranksize allometric scaling law, it will follow the hierarchical allometric scaling law [18]. Thus, we can apply the MLEbased approach to Zipf’s distributions of urban population and area in different years, respectively. Zipf’s distribution is mathematically equivalent to Pareto distribution [25]. After the Pareto exponents are estimated, we can estimate the allometric scaling exponent based on maximumlikelihood fitting. For a powerlaw distribution function , the corresponding density distribution can be obtained by differential as follows:where is cumulative distribution function (CDF), is density distribution function (DDF), denotes the Pareto exponent; that is, the cumulative scaling exponent, represents the density scaling exponent, and refers to differential operator and to the proportional relation. Replacing the general unknown quantity by urban population () and area () yields two density distribution functions as below:in which and represent the density scaling exponents of urban population and area distributions and the other symbols are in essence the same as in (4) and (5). Thus the allometric scaling exponent can be estimated by the formula , as shown above. For the purpose of comparison, both the MLEbased approach and the OLSbased approach are applied to the datasets of Zipf’s distributions of Chinese urban population and area from 1991 to 2014. The estimated values of the allometric scaling exponents are tabulated as Table 5.
 
Note. The MLEbased allometric scaling exponent and the OLSbased allometric scaling exponent values are both based on Zipf’s distributions, and are estimated by the ranksize allometric scaling relation. values are different from the results shown in Table 2, which is based on hierarchical allometric scaling. 
The MLEbased results lend further support to the findings shown in Section 3.2 (Results and Findings). In theory, the reasonable density scaling exponent comes between 1.5 and 3, while the logical allometric scaling exponent varies from 2/3 to 1. However, from 1991 to 2002, the MLEbased allometric scaling exponent values are abnormal, ranging from 0.47 to 2. Accordingly, the OLSbased results are also abnormal, ranging from 0.5 to 0.74. There are significant difference between the MLEbased results and the OLSbased results (Table 5). From 2003 to 2014, the MLEbased results and the OLSbased results seem to reach the same goal by different routes (Figure 8). As discussed above, the year 2003 is a turning point. Before 2003, the allometric scaling was not significant. From 2003 on, the allometric scaling relation between urban population and area emerged. The findings are as follows. First, the MLEbased method developed by Clauset et al. [34] is rather more sensitive to the powerlaw relation than the OLSbased method. Second, the similarity and differences between the MLEbased results and the OLSbased results are helpful for us to examine the evolution of allometric scaling. Third, as a whole, the MLEbased method is not suitable for this study, but the two methods can supplement each other.
4. Discussion
Hierarchical scaling suggests a universal law of complex systems such as cities. The scaling law is followed by central place network [21], river composition [51–54], earthquake energy distribution [55, 56], animals’ blood vessels [57–59], and so on. Based on the case studies of hierarchy of Chinese cities, a complete process of allometric scaling analysis is developed. The traditional forms of Zipf’s law and allometric growth law can be integrated into the new analytical framework (Figure 9). Using the hierarchical scaling method, we can research city development and urban evolution from a new angle of view. By the empirical analysis, we obtain an insight into Chinese urban system. First, Zipf’s law and the allometric growth law are not “iron laws” for Chinese cities. The development of allometric scaling in Chinese cities can be divided into three stages. The first stage is the abnormal allometry phase (before 1992), during which the urban populationarea relation follow the power law, but this local allometric exponent exceeded the proper upper limit (). The second stage is the allometrybreaking phase (1993–2002), during which the global allometric scaling or even the local allometric scaling degenerated to linear relations. The third stage is normal allometry phase (after 2003), during which the urban populationarea relation follows the power law, and the allometric exponent comes between 2/3 and 1 (). Second, the scaling relations are significantly influenced by fast urbanization and real estate industry of China. Fast urbanization results in rapid growth of urban population, and realty industry leads to builtup area change hastily. The relationships between urbanization and real estate industry were imprinted on Chinese city development. During the period of fast urbanization, the real estate bubble led to many rounds of “creates city” campaigns in China. Third, two types of forces impact on urban evolution of China. One is the topdown force coming from command economy and government intervention, and the other is the bottomup force coming from market economy and individual actions. The former suggests the Invisible Hand, and the latter indicates the Visible Hand. In this case, the allometric scaling and ranksize patterns of Chinese cities seem to struggle between order and chaos. Despite all these, the allometric scaling exponent is very stable in the main. In fact, a scaling exponent is often close to a ratio such as 1/2, 2/3, 3/4, 4/5, and 5/6 [18, 47]. The allometric exponent values of Chinese cities seemed to change around 6/7 from 1991 to 2014.
Compared with previous studies on the allometry and scaling of Chinese cities, this work bears significant characteristics. Firstly, the early studies are based on the twoparameter Zipf law, but Chinese cities follow the threeparameter Zipf’s law. Using hierarchical scaling to replace the ranksize scaling, we can catch the characters of the threeparameter Zipf distributions [29]. Secondly, the early studies are based on the ranksize distribution instead of hierarchical structure of urban systems. This paper is based on datasets of cascade sequences abstracting from selfsimilar hierarchy of cities. In this instance, the allometric scaling can be revealed obviously. By reconstructing a hierarchy with cascade structure, we can reduce the random disturbance of noises in observational data. In particular, the ranksize patterns can be brought to light. The allometric scaling relation comes from Zipf’s distributions of urban population and area. Allometric growth law can provide circumstantial evidence of Zipf’s law. Maybe some persons prefer the common ranksize scaling to the hierarchical scaling analysis, because the traditional method is simple and the sample of cities seems to be bigger. As a matter of fact, the effect of statistical analysis depends on the association of degree of freedom with level of significance rather than sample size. As indicated above, the common Zipf’s law and allometric growth law are contained in the hierarchical rescaling process. Third, two approaches can be employed to estimate the allometric scaling exponents based on selfsimilar hierarchies. One is the OLSbased approach, and the other is the MLEbased approach. The former can be directly applied to this allometric scaling analysis, while the latter can only be indirectly applied to the hierarchical series data. The MLEbased method proposed by Clauset et al. [34] is designed for binned data rather than cascade sequences. If the allometric scaling relation of cities is well developed, the two approaches will lead to similar computational results. However, if the allometric scaling is not well developed, the OLSbased approach is better than the MLEbased approach. This suggests that the direct approach is better than the indirect approach.
The deficiencies of this study are as follows. First, the sample comprises about 660 officially approved cities. In China, there are thousands of human settlements that can be treated as cities. However, only the observational data of the officially approved cities are available for quantitative studies. Perhaps just because of this, the one or twoparameter Zipf distribution was replaced by the threeparameter Zipf distribution. Second, the data are based on statistical investigation rather than census. On the one hand, the census data are not continuous in time; on the other, we cannot find the land survey data for urbanized area. What is more, city population is not just within the limit of urban area. In fact, the datasets of natural data have better quality [40, 60]. These data have been employed to make allometric scaling analysis and the results are satisfying [61]. From the datasets of natural cities of America, England, France, and Germany, we can derive the same hierarchical allometric scaling as that of China. Unfortunately, we have no long sample path data of urban area and population from the time series of natural cities. Thus, only the crosssectional analysis rather than the dynamic analysis can be made for the time being. Third, the hierarchy is constructed by city number rather than population size. There are two approaches to reconstructing a hierarchy of human settlements [31]: one is by size presented by Davis [30], and the other is by number developed by Chen [4]. For simplicity, this study is only focused on numberbased hierarchy. Fourth, the hierarchical allometric rescaling is based on population size rather than urban area. Urban population is not always consistent with urban area; that is, greater than does not necessarily indicate that is greater than . We can rank cities either by population size or by urban area (Figure 10). There is subtle difference between the two results of hierarchies. Due to limited space, the areabased hierarchy is not discussed in this article. Despite all these shortcomings, the allometric scaling and its evolution patterns of Chinese cities are brought to light.
(a) Ranked by area
(b) Ranked by population
5. Conclusions
A new analytical framework is presented by means of the case study of Chinese cities. In this framework, Zipf’s law and allometric growth law are integrated into hierarchical scaling analysis of urban systems. Using the analytical process, we brought to light the spatiotemporal evolution properties of Chinese cities. By this study, we can get insight into the ranksize patterns and allometric scaling of Chinese system of cities. The main points of this work are as follows. First, Chinese cities follow allometric scaling law that is based on Zipf’s law. However, the ranksize pattern and allometric scaling emerged from the nonlinear dynamics of urbanization. On the whole, the hierarchy of Chinese cities follows the allometric scaling law. However, from 1998 to 2000, the powerlaw relation degenerated to a linear relation. From 2003 on, the ranksize scaling and allometric scaling become clearer and more significant over time. This lends further support to the suggestion that both the allometric scaling and the ranksize scaling are evolutional laws. Second, the allometric scaling exponent reflects the characters of the ranksize distribution and the manland relation of Chinese cities. Where the global level based on all the cities is concerned, the allometric scaling exponent values were relatively stable and approached the theoretically expected value 0.85. However, where the local level is based on the cities within the scaling range, the allometric scaling exponent values departed from 0.85 and approached 0.93. Because of fast urbanization and real estate industry, the allometric exponent values took on a rising trend for a time. This implies that the per capita landuse area in Chinese cities does not fall but rises due to the increase of city sizes. In fact, the real estate bubble economy results in the waste of land resource utilization in China. After 2009, the local allometric exponent trended down to 0.85. Third, the allometric scaling of Chinese cities is associated with Zipf’s distributions of urban population and area. If and only if the urban population size and urban area follow Zipf’s law, the allometric scaling relation between urban area and population emerges. The ranksize allometry can be converted into hierarchical allometry of urban systems. Before 2003, the allometric scaling relation degenerated to a linear relation on one occasion because the city population size failed to comply with Zipf’s law. In this study, the linear relationships are technically treated as powerlaw relations. From 2003 onwards, both the urban population and area conform to Zipf’s law, and thus the ranksize allometry and hierarchical allometry emerged from the process of urbanization. Both the OLSbased results and MLEbased results support this conclusion meanwhile.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was sponsored by the National Natural Science Foundations of China (Grant nos. 41590843 and 41671167). The support is gratefully acknowledged.
References
 M. Batty and P. A. Longley, Fractal Cities: A Geometry of form And Function, Academic Press, London, UK, 1994.
 L. Bettencourt, “The origins of scaling in cities,” American Association for the Advancement of Science. Science, vol. 340, no. 6139, pp. 1438–1441, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 L. M. A. Bettencourt, J. Lobo, D. Helbing, C. Kühnert, and G. B. West, “Growth, innovation, scaling, and the pace of life in cities,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 17, pp. 7301–7306, 2007. View at: Publisher Site  Google Scholar
 Y. G. Chen, Fractal Urban Systems: Scaling, Symmetry, and Spatial Complexity, Science Press, Beijing, China, 2008 (Chinese).
 P. Frankhauser, La Fractalité des Structures Urbaines (The Fractal Aspects of Urban Structures), Economica, Paris, France, 1994.
 J. Lobo, L. M. A. Bettencourt, D. Strumsky, and G. B. West, “Urban scaling and the production function for cities,” PLoS ONE, vol. 8, no. 3, Article ID e58407, 2013. View at: Publisher Site  Google Scholar
 D. Rybski, S. V. Buldyrev, S. Havlin, F. Liljeros, and H. A. Makse, “Scaling laws of human interaction activity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 31, pp. 12640–12645, 2009. View at: Publisher Site  Google Scholar
 C. Carroll, “National citysize distributions: what do we know after 67 years of research?” Progress in Human Geography, vol. 6, no. 1, pp. 1–43, 1982. View at: Google Scholar
 X. Gabaix, “Zipf's law for cities: an explanation,” Quarterly Journal of Economics, vol. 114, no. 3, pp. 739–767, 1999. View at: Publisher Site  Google Scholar
 X. Gabaix and Y. M. Ioannides, “The evolution of city size distributions,” in Handbook of Urban and Regional Economics, J. V. Henderson and J. F. Thisse, Eds., vol. 4, pp. 2341–2378, NorthHolland Publishing Company, Amsterdam, Netherlands, 2004. View at: Publisher Site  Google Scholar
 P. Krugman, “Confronting the mystery of urban hierarchy,” Journal of the Japanese and International Economies, vol. 10, no. 4, pp. 399–418, 1996. View at: Publisher Site  Google Scholar
 G. K. Zipf, Human Behavior and the Principle of Least Effort, AddisonWesley, Reading, Mass, USA, 1949.
 M. Batty, “The size, scale, and shape of cities,” Science, vol. 319, no. 5864, pp. 769–771, 2008. View at: Publisher Site  Google Scholar
 G. Dutton, “Foreword: size and shape in the growth of human communities,” Ekistics, vol. 36, pp. 142–243, 1973. View at: Google Scholar
 S. Nordbeck, “Urban allometric growth,” Geografiska Annaler B, vol. 53, no. 1, pp. 54–67, 1971. View at: Publisher Site  Google Scholar
 Y. Lee, “An allometric analysis of the US urban system: 1960–1680,” Environment and Planning A, vol. 21, no. 4, pp. 463–476, 1989. View at: Publisher Site  Google Scholar
 C. P. Lo and R. Welch, “Chinese urban population estimates,” Annals of the Association of American Geographers, vol. 67, no. 2, pp. 246–253, 1977. View at: Publisher Site  Google Scholar
 Y. Chen, “The spatial meaning of Pareto's scaling exponent of citysize distributions,” Fractals, vol. 22, no. 12, Article ID 14500017, 2014. View at: Publisher Site  Google Scholar
 Y. Chen, “The mathematical relationship between Zipf's law and the hierarchical scaling law,” Physica A: Statistical Mechanics and its Applications, vol. 391, no. 11, pp. 3285–3299, 2012. View at: Publisher Site  Google Scholar
 D. Pumain, Hierarchy in Natural and Social Sciences, vol. 3, SpringerVerlag, Dordrecht, Netherlands, 2006. View at: Publisher Site
 Y. Chen, “Zipf's law, hierarchical structure, and cardsshuffling model for urban development,” Discrete Dynamics in Nature and Society, vol. 2012, Article ID 480196, 21 pages, 2012. View at: Publisher Site  Google Scholar
 G. Anderson and Y. Ge, “The size distribution of Chinese cities,” Regional Science and Urban Economics, vol. 35, no. 6, pp. 756–776, 2005. View at: Publisher Site  Google Scholar
 L. Benguigui and E. BlumenfeldLieberthal, “A dynamic model for city size distribution beyond Zipf 's law,” Physica A: Statistical Mechanics and its Applications, vol. 384, no. 2, pp. 613–627, 2007. View at: Publisher Site  Google Scholar
 L. Benguigui and E. BlumenfeldLieberthal, “Beyond the power law—a new approach to analyze city size distributions,” Computers, Environment and Urban Systems, vol. 31, no. 6, pp. 648–666, 2007. View at: Publisher Site  Google Scholar
 Y. Chen, R. Chen, N. S. Ai, and H. Q. Li, “On the fractal property of citysize distributions,” Economical Geography, vol. 13, no. 3, pp. 48–53, 1993 (Chinese). View at: Google Scholar
 K. Gangopadhyay and B. Basu, “City size distributions for India and China,” Physica A: Statistical Mechanics and Its Applications, vol. 388, no. 13, pp. 2682–2688, 2009. View at: Publisher Site  Google Scholar
 X. Y. Ye and Y. C. Xie, “Reexamination of Zipf’s law and urban dynamics in China: a regional approach,” The Annals of Regional Science, vol. 49, no. 1, pp. 135–156, 2012. View at: Google Scholar
 P. Bak, How Nature Works: The Science of SelfOrganized Criticality, Springer, New York, NY, USA, 1997.
 Y. Chen, “The evolution of Zipf's law indicative of city development,” Physica A: Statistical Mechanics and Its Applications, vol. 443, pp. 555–567, 2016. View at: Publisher Site  Google Scholar
 K. Davis, “World urbanization: 1950–1970,” in Systems of Cities, I. S. Bourne and J. W. Simons, Eds., pp. 92–100, Oxford University Press, New York, NY, USA, 1978. View at: Google Scholar
 B. Jiang and X. Yao, Geospatial Analysis and Modeling of Urban Structure and Dynamics, Springer, Berlin, Germany, 2010.
 M. J. Beckmann, “City hierarchies and distribution of city sizes,” Economic Development and Cultural Change, vol. 6, no. 3, pp. 243–248, 1958. View at: Google Scholar
 M. E. J. Newman, “Power laws, Pareto distributions and Zipf's law,” Contemporary Physics, vol. 46, no. 5, pp. 323–351, 2005. View at: Publisher Site  Google Scholar
 A. Clauset, C. R. Shalizi, and M. E. Newman, “Powerlaw distributions in empirical data,” SIAM Review, vol. 51, no. 4, pp. 661–703, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 R. L. Axtell, “Zipf distribution of U.S. firm sizes,” Science, vol. 293, no. 5536, pp. 1818–1820, 2001. View at: Publisher Site  Google Scholar
 M. Batty, “Rank clocks,” Nature, vol. 444, no. 7119, pp. 592–596, 2006. View at: Publisher Site  Google Scholar
 X. Gabaix, “Zipf's law and the growth of cities,” American Economic Review, vol. 89, no. 2, pp. 129–132, 1999. View at: Publisher Site  Google Scholar
 A. Hernando, D. Puigdomènech, D. Villuendas, C. Vesperinas, and A. Plastino, “Zipf's law from a Fisher variationalprinciple,” Physics Letters, Section A: General, Atomic and Solid State Physics, vol. 374, no. 1, pp. 18–21, 2009. View at: Publisher Site  Google Scholar
 B. H. Hong, K. E. Lee, and J. W. Lee, “Power law in firms bankruptcy,” Physics Letters, Section A: General, Atomic and Solid State Physics, vol. 361, no. 12, pp. 6–8, 2007. View at: Publisher Site  Google Scholar
 B. Jiang and T. Jia, “Zipf's law for all the natural cities in the United States: a geospatial perspective,” International Journal of Geographical Information Science, vol. 25, no. 8, pp. 1269–1281, 2011. View at: Publisher Site  Google Scholar
 J. Shao, P. C. Ivanov, B. Uroević, H. E. Stanley, and B. Podobnik, “Zipf rank approach and crosscountry convergence of incomes,” Europhysics Letters, vol. 94, no. 4, Article ID 48001, 2011. View at: Publisher Site  Google Scholar
 M. H. R. Stanley, S. V. Buldyrev, S. Havlin, R. N. Mantegna, M. A. Salinger, and H. E. Eugene, “Zipf plots and the size distribution of firms,” Economics Letters, vol. 49, no. 4, pp. 453–457, 1995. View at: Publisher Site  Google Scholar
 D. Pumain and F. MoriconiEbrard, “City size distributions and metropolisation,” GeoJournal, vol. 43, no. 4, pp. 307–314, 1997. View at: Publisher Site  Google Scholar
 M. J. Woldenberg, “An allometric analysis of urban land use in the United States,” Ekistics, vol. 36, pp. 282–290, 1973. View at: Google Scholar
 J. Feng and Y. G. Chen, “Spatiotemporal evolution of urban form and landuse structure in Hangzhou, China: evidence from fractals,” Environment and Planning B: Planning and Design, vol. 37, no. 5, pp. 838–856, 2010. View at: Publisher Site  Google Scholar
 Y. G. Chen, “Characterizing growth and form of fractal cities with allometric scaling exponents,” Discrete Dynamics in Nature and Society, vol. 2010, Article ID 194715, 22 pages, 2010. View at: Publisher Site  Google Scholar
 Y. G. Chen, “An allometric scaling relation based on logistic growth of cities,” Chaos, Solitons and Fractals, vol. 65, pp. 65–77, 2014. View at: Google Scholar
 H. E. Hurst, R. P. Black, and Y. M. Simaika, Longterm Storage: An Experimental Study, Constable, London, UK, 1965.
 J. Feder, Fractals, Plenum Press, New York, NY, USA, 1988. View at: Publisher Site  MathSciNet
 B. B. Mandelbrot, The Fractal Geometry of Nature, W. H. Freeman and Company, New York, NY, USA, 1982.
 R. E. Horton, “Erosional development of streams and their drainage basins; hydrophysical approach to quantitative morphology,” Bulletin of the Geological Society of America, vol. 56, no. 3, pp. 275–370, 1945. View at: Publisher Site  Google Scholar
 I. RodriguezIturbe and A. Rinaldo, Fractal River Basins: Chance and SelfOrganization, Cambridge University Press, Cambridge, UK, 2001.
 S. A. Schumm, “Evolution of drainage systems and slopes in badlands at Perth Amboy, New Jersey,” Bulletin of the Geological Society of America, vol. 67, no. 5, pp. 597–646, 1956. View at: Publisher Site  Google Scholar
 A. N. Strahler, “Hypsometric (areaaltitude) analysis of erosional topography,” Bulletin of the Geological Society of America, vol. 63, no. 11, pp. 1117–1142, 1952. View at: Publisher Site  Google Scholar
 B. Gutenberg and C. F. Richter, Seismicity of the Earth and Associated Phenomenon, Princeton University Press, Princeton, NJ, USA, 2nd edition, 1954.
 D. L. Turcotte, Fractals and Chaos in Geology and Geophysics, Cambridge University Press, 1997. View at: Publisher Site  MathSciNet
 Y. G. Chen, “Fractals and fractal dimension of systems of blood vessels: an analogy between artery trees, river networks, and urban hierarchies,” Fractal Geometry and Nonlinear Analysis in Medicine and Biology, vol. 1, pp. 26–32, 2015. View at: Google Scholar
 Z. L. Jiang and G. C. He, “Geometrical morphology of the human coronary arteries,” Journal of Third Military Medical University, vol. 11, no. 2, pp. 85–91, 1989 (Chinese). View at: Google Scholar
 Z. L. Jiang and G. C. He, “Geometrical morphology of coronary arteries in dog,” Chinese Journal of Anatomy, vol. 13, no. 3, pp. 236–241, 1990 (Chinese). View at: Google Scholar
 B. Jiang and X. Liu, “Scaling of geographic space from the perspective of city and field blocks and using volunteered geographic information,” International Journal of Geographical Information Science, vol. 26, no. 2, pp. 215–229, 2011. View at: Publisher Site  Google Scholar
 Y. G. Chen and B. Jiang, “Hierarchical scaling in systems of natural cities,” https://arxiv.org/abs/1608.05770. View at: Google Scholar
Copyright
Copyright © 2017 Yanguang Chen and Jian Feng. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.