#### 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 long-term 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 self-organized 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 scale-free 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 distance-decay law compose three basic laws of urban geography. Zipf’s law indicates the rank-size 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 rank-size 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 rank-size 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 rank-size pattern emerges from the edge of chaos [28]. A new discovery is that Chinese city-size distribution can be described by the three-parameter Zipf’s model instead of the two-parameter Zipf’s model [29]. The three-parameter 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 rank-size 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 self-similar 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 city-size distribution. Suppose that the cities in a geographical region are grouped into* M* classes in a top-down 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 top-down 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 top-level cities, and are the mean population size and urban area of the top-level 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 Beckmann-Davis 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 power-law relations including the three-parameter Zipf’s law can be derived from the above exponential laws [21, 29]. The power-law models are as below:where , /, , , , and . Equation (4) is termed the hierarchical* size-number scaling relation* of cities, which is equivalent to the Pareto law of city-size distribution, and* D* is the fractal dimension of urban hierarchies measured with urban population. Equation (5) is termed the hierarchical* area-number scaling relation* of cities, which is equivalent to the Pareto law of city-area 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 three-parameter Zipf’s model. The Zipf distribution is a signature of self-similar hierarchies associated with fractal patterns and self-organized processes.

#### 3. Empirical Analysis

##### 3.1. Study Area and Methods

The models of hierarchical allometry can be applied to the rank-size 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 Urban-Rural Development of China (MOHURD).

A number of algorithms can be employed to estimate the power exponent values, including the least squares method (LSM), maximum-likelihood 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 city-size 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 rank-size 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 self-similar hierarchy. The two-parameter Zipf models of urban population and area can be expressed aswhere denotes the rank of a city, and refer to the population size and built-up 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) Built-up area**

Then, we investigate the allometric scaling relation based on the rank-size distribution. This is a kind of transversal allometry based on cross-sectional 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 built-up area of Chinese cities yields models of cross-sectional 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 power-law 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 cross-sectional 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 self-similar hierarchy can be constructed as follows. Suppose that we divide all the cities into* M* classes/levels in a top-down 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].

**(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 self-organized process of urbanization.* From 1991 to 1992 year, the urban population-area relation follows the power-law 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 power-law 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 medium-sized 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 rank-size 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 rank-population distribution was not consistent with the rank-area distribution. As a result, the allometric scaling relation of Chinese cities degenerated from power-law relation to quasilinear relation. This lends further support to the inference that urban area-population 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 land-use 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 scale-free range in which the allometric relation is more significant*. On log-log 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 rank-size distributions and the corresponding size-number 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 rank-size distribution or size-number 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.

**(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 self-affine 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 Maximum-Likelihood 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 maximum-likelihood estimation (MLE) is developed by Clauset et al. [34] to address the power-law distribution in empirical data. Combining the methods of maximum-likelihood fitting with the goodness-of-fit tests based on the Kolmogorov-Smirnov statistic and likelihood ratios, the approach is employed by many scholars to identify various power-law distributions in both natural and social sciences. Unfortunately, the MLE-based approach cannot be directly applied to our datasets in this work. The reasons are as follows. First, the MLE-based method is developed for power-law frequency distributions, while this study is devoted to allometric scaling relations. Second, the MLE-based method is suitable for binned data, while this study is based on pairs of cascade sequences.

However, the MLE-based 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 rank-size allometry [29, 46]. The MLE-based 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 area-population relation follows the rank-size allometric scaling law, it will follow the hierarchical allometric scaling law [18]. Thus, we can apply the MLE-based 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 maximum-likelihood fitting. For a power-law 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 MLE-based approach and the OLS-based 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.

The MLE-based 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 MLE-based allometric scaling exponent values are abnormal, ranging from 0.47 to 2. Accordingly, the OLS-based results are also abnormal, ranging from 0.5 to 0.74. There are significant difference between the MLE-based results and the OLS-based results (Table 5). From 2003 to 2014, the MLE-based results and the OLS-based 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 MLE-based method developed by Clauset et al. [34] is rather more sensitive to the power-law relation than the OLS-based method. Second, the similarity and differences between the MLE-based results and the OLS-based results are helpful for us to examine the evolution of allometric scaling. Third, as a whole, the MLE-based 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 population-area relation follow the power law, but this local allometric exponent exceeded the proper upper limit (). The second stage is the allometry-breaking 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 population-area 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 built-up 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 top-down force coming from command economy and government intervention, and the other is the bottom-up 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 rank-size 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 two-parameter Zipf law, but Chinese cities follow the three-parameter Zipf’s law. Using hierarchical scaling to replace the rank-size scaling, we can catch the characters of the three-parameter Zipf distributions [29]. Secondly, the early studies are based on the rank-size distribution instead of hierarchical structure of urban systems. This paper is based on datasets of cascade sequences abstracting from self-similar 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 rank-size 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 rank-size 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 self-similar hierarchies. One is the OLS-based approach, and the other is the MLE-based 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 MLE-based 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 OLS-based approach is better than the MLE-based 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 two-parameter Zipf distribution was replaced by the three-parameter 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 cross-sectional 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 number-based 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 area-based 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 rank-size 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 rank-size 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 power-law relation degenerated to a linear relation. From 2003 on, the rank-size 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 rank-size scaling are evolutional laws.* Second, the allometric scaling exponent reflects the characters of the rank-size distribution and the man-land 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* land-use 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 rank-size 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 power-law relations. From 2003 onwards, both the urban population and area conform to Zipf’s law, and thus the rank-size allometry and hierarchical allometry emerged from the process of urbanization. Both the OLS-based results and MLE-based 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.