Understanding Chinese Urban Form: The Universal Fractal Pattern of Street Networks over 298 Cities

: Urban form can be reﬂected by many city elements, such as streets. A street network serves as the backbone of a city and reﬂects a city’s physical structure. A street network’s topological measures and statistical distributions have been widely investigated in recent years, but previous studies have seldom characterized the heavy-tailed distribution of street connectivities from a fractal perspective. The long-tail distribution of street connectivities can be fractal under the new, third deﬁnition: a set or pattern is fractal if the scaling of far more small things than large ones recurs at least twice. The number of recurred scaling patterns of far more less-connected streets than well-connected ones greatly helps in measuring the scaling hierarchy of a street network. Moreover, it enables us to examine the potential fractality of urban street networks at the national scale. In this connection, the present study aims to contribute to urban morphology in China through the investigation of the ubiquity of fractal cities from the lens of street networks. To do this, we generate hundreds of thousands of natural streets from about 4.5 million street segments over 298 Chinese cities and adopted power-law detection as well as three fractal metrics that emerged from the third deﬁnition of fractal. The results show that almost all cities bear a fractal structure in terms of street connectivities. Furthermore, our multiple regression analysis suggests that the fractality of street networks is positively correlated with urban socioeconomic status and negatively correlated with energy consumption. Therefore, the fractal metrics can be a useful supplement to traditional street-network conﬁguration measures such as street lengths.


Introduction
It is widely recognized that geographic features such as mountains and rivers are neither smooth nor regular, so they cannot be precisely described on the basis of Euclidean geometry. Instead, fractal geometry has the capacity to fill this gap. The very word "fractal" originates from Latin and literally means "broken" or "irregular". The development of fractal geometry went through three definitions, all of which center on the term "self-similarity": a part of an object or pattern is similar to the whole [1][2][3]. For the first and second definitions, fractals refer to shapes that are either strictly or statistically self-similar. A classic example of strictly self-similar shapes is a Koch curve [4], whose parts are exactly the scaled copies of the original curve. Mandelbrot [5] relaxed the strictness by adding some randomness to those "scaled copies" to make the constructed curve capable of describing real-world features, such as a coastline. Strict and statistical self-similarities both maintain the same power law Section 2.2) for the identification of the fractal or scaling structure. We have successfully processed and derived the natural streets for almost all prefectural-level cities in China, which could be a valuable data source for the research community. Second, we designed an innovative methodology in combining different fractal metrics under the theoretical framework of the third definition. Starting from the head/tail breaks-induced ht-index [6], we can assess the fractality of the street network for each individual city and the overall complexity of all street networks at a country scale. Next, due to the lack of sensitivity of the ht-index, we employed two alternatives: the cumulative rate of growth (CRG, [30]) and the ratio of areas in a rank-size plot (RA, [31]) (see more details in Section 2.3). We found that using them collectively can differentiate cities within the same scaling hierarchy, making it possible to formulate a more comprehensive and systematic understanding of Chinese urban forms. Lastly, we sought relationships between the fractality of street networks and different kinds of urban metrics and evaluated the explanatory power of urban fractals on urban socioeconomic status and energy consumption at the national scale.
The remainder of this paper proceeds as follows. Section 2 introduces the data sets and the proposed methodological framework. Section 3 presents the visualization and statistical results regarding the universal fractal pattern over nearly 300 cities in China, as well as the correlations between each fractal metric and urban quantity. Section 4 further discusses the fractal structure of street networks, before Section 5 draws conclusions and points to future research directions.

Data and Data Processing
Three data sets were used in this study: (1) 298 city boundaries, (2) the national street network, and (3) China grid data sets of GDP, population, and CO 2 emissions. As illustrated in Figure 1, we selected 298 prefectural-level cities that were spread all over China as the study areas. Each city boundary is further adopted as the unit for data processing. The national street network was downloaded from a Germany free Geodata platform GEOFABRIK [32]. We used ArcGIS Interoperability Toolbox to process the large collection of polylines with all road types into 4,419,603 segments to ensure that each segment's end vertex is the street junction where at least three polylines intersect, and then we clipped the entire set of segments using each urban boundary to obtain the city-level streets. The national GDP and population grid data were sourced from the National Resources and Environment Database of the Chinese Academy of Sciences [33]. The data were collected in 2015 and spatialized by remote sensing technology with a ground resolution of 1 km (for more details of this approach, see [34]). City-level GDP and population statistics were calculated successfully by aggregating the cell values within each city boundary. The city-level CO 2 emissions were derived directly from the previous work [35], which provides related statistics based on CHRED (China high-resolution emission database; see more details in [36]). Note that the urban-scale CO 2 emissions in Cai et al.'s [35] work are only available for 283 cities.

Transformation from Segment-Segment Topology to Street-Street Topology
Segment-segment topology is widely adopted in current GIS systems. For example, the ArcGIS Network Analyst extension builds up the network analysis based on the relationship between segments that are separated at a street junction node. The segment-segment topology, or segment-based analysis in general, is remarkable because it offers many urban applications, such as navigation. However, the segment-segment topology is essentially a geometric perspective because the modeling unit of the street network is based on geometric primitives, such as points, polylines, and polygons. The geometric primitives represent space mechanically, as their lengths or connectivities tend to follow a normal-like distribution (Figure 2a,b). ISPRS Int. J. Geo-Inf. 2020, 9,192 4 of 17

Transformation from Segment-Segment Topology to Street-Street Topology
Segment-segment topology is widely adopted in current GIS systems. For example, the ArcGIS Network Analyst extension builds up the network analysis based on the relationship between segments that are separated at a street junction node. The segment-segment topology, or segmentbased analysis in general, is remarkable because it offers many urban applications, such as navigation. However, the segment-segment topology is essentially a geometric perspective because the modeling unit of the street network is based on geometric primitives, such as points, polylines, and polygons. The geometric primitives represent space mechanically, as their lengths or connectivities tend to follow a normal-like distribution (Figure 2a, b).
In order to view the space more organically, one can consider transforming the segmentsegment topology into the street-street topology. A street refers to the notion of a "natural street" or "named street", which is made up of a set of neighboring segments with good continuity or shared names [37]. The good continuity can be understood as small deflection angles (that is, those smaller than 45 degrees) across adjacent segment pairs. Compared to segments, streets such as the modeling unit are coherent and meaningful in our perception [22]. Moreover, the street-street topology allows us to investigate the street network from a complex-network perspective [38], as the street lengths and connectivities under such a setting tend to be scale-free or heterogeneous, indicating far more short/less-connected streets than long/well-connected ones (Figure 2d). Although it is possible that street length correlates with street connectivity, we have only used street connectivity for the following experiments because geometric lengths are not a major factor for understanding the urban In order to view the space more organically, one can consider transforming the segment-segment topology into the street-street topology. A street refers to the notion of a "natural street" or "named street", which is made up of a set of neighboring segments with good continuity or shared names [37]. The good continuity can be understood as small deflection angles (that is, those smaller than 45 degrees) across adjacent segment pairs. Compared to segments, streets such as the modeling unit are coherent and meaningful in our perception [22]. Moreover, the street-street topology allows us to investigate the street network from a complex-network perspective [38], as the street lengths and connectivities under such a setting tend to be scale-free or heterogeneous, indicating far more short/less-connected streets than long/well-connected ones (Figure 2d). Although it is possible that street length correlates with street connectivity, we have only used street connectivity for the following experiments because geometric lengths are not a major factor for understanding the urban space. This can be easily understood by the fact that a well-connected street is commonly found in a city, whereas a long street is often not a well-connected street and is usually found in the countryside. ISPRS Int. J. Geo-Inf. 2020, 9, 192 5 of 17 space. This can be easily understood by the fact that a well-connected street is commonly found in a city, whereas a long street is often not a well-connected street and is usually found in the countryside. (Note: From the network perspective of street segments and streets, segment-segment topology (Panels a and b) is a mechanistic structure that lacks the scaling pattern of far more smalls than larges, whereas the street-street topology (Panels c and d) is organic, as it holds such a pattern clearly).

Fractal Analysis of Street Connectivities
The fractal structure of urban street networks can be characterized by power law or other heavytailed distribution models regarding the street connectivities. As the power law model may be too strict to identify a fractal, we derive the scaling hierarchies of streets using head/tail breaks [39] and compute three metrics-the ht-index, CRG, and RA-to measure the potential fractality of the street network. (1) The simplest way to identify a power law distribution is by ranking the data values from largest to smallest and then creating a log-log plot of the sorted data. The distribution will be a straight line if the data exactly follows a power-law distribution. However, it is very difficult for this method to deal with the messy up-and-downs in the tail for many real-world data sets [40]. To address this problem, a robust method using maximum likelihood estimation, proposed by Clauset et al. [41], can help detect accurately whether data is a power law or not. The method firstly examines the power law exponent α, which is denoted by Equation (2): where is the smallest value from where the data is power-law distributed. The acceptable range for exponent α is from 1 to 3. To test how good the data can fit with a power law, Clauset et al. [41] suggest doing a Kolmogorov-Smirnoff test. This test compares the data with the samples based on an ideal power law distribution, using the derived x-min and alpha values. The goodness-of-fit index (Note: From the network perspective of street segments and streets, segment-segment topology (Panels a and b) is a mechanistic structure that lacks the scaling pattern of far more smalls than larges, whereas the street-street topology (Panels c and d) is organic, as it holds such a pattern clearly).

Fractal Analysis of Street Connectivities
The fractal structure of urban street networks can be characterized by power law or other heavy-tailed distribution models regarding the street connectivities. As the power law model may be too strict to identify a fractal, we derive the scaling hierarchies of streets using head/tail breaks [39] and compute three metrics-the ht-index, CRG, and RA-to measure the potential fractality of the street network.

Power Law Detection
A power law model can be expressed by Equation (1): The simplest way to identify a power law distribution is by ranking the data values from largest to smallest and then creating a log-log plot of the sorted data. The distribution will be a straight line if the data exactly follows a power-law distribution. However, it is very difficult for this method to deal with the messy up-and-downs in the tail for many real-world data sets [40]. To address this problem, a robust method using maximum likelihood estimation, proposed by Clauset et al. [41], can help detect accurately whether data is a power law or not. The method firstly examines the power law exponent α, which is denoted by Equation (2): where x min is the smallest value from where the data is power-law distributed. The acceptable range for exponent α is from 1 to 3. To test how good the data can fit with a power law, Clauset et al. [41] suggest doing a Kolmogorov-Smirnoff test. This test compares the data with the samples based on an ideal power law distribution, using the derived x-min and alpha values. The goodness-of-fit index p-value, ranging from 0 to 1, is then generated and used to determine the extent of fitness for the data to a power law model. If p-value ≥ 0.01, we accept a set of data being power-law distributed.

Integration of Ht-Index, CRG, and RA for Fractal Measurement
Jiang and Yin [6] proposed the third definition of fractal as "the recurrence of far more smalls than larges" and induced a new metric ht-index, coming from the new classification scheme head/tail breaks [39]. To elaborate, given data with a heavy-tailed distribution, we can always split the data into a great many small values (that is, the tail) and a few large values (the head). Note that there is an imbalance between head percentage and tail percentage (for example, 40/60). The splitting operation is called head/tail division. We can apply this division iteratively to the head part of the data until the updated head and tail percentages reach a balance, such as 50/50. The entire process demonstrates how the head/tail breaks method [39] works. Since the process is recursive, we would normally achieve a series of mean values. The number of obtained mean values i plus one is the so-called ht-index, as expressed by Equation (3): Figure 3a shows the first three iterations of a Koch curve [4]. The first iteration contains one segment with a scale of 1; the second contains four segments with a scale of 1/3; the third contains 16 segments with a scale of 1/9. The fractal dimension, i.e., the power law exponent between its details and scales, equals to the slope of log-log plot as shown in Figure 3b. To apply the head/tail breaks method, we rank all scales of the Koch curve from the largest to the smallest (a rank-size plot, Figure 3c). The induced ht-index value of the Koch curve is 3, meaning that the scaling pattern of far more short segments than long ones appears twice. Two mean values were computed: In this way, we can derive the curve's three scaling hierarchical levels: 16 segments whose scales are below m 2 , four segments whose scales are between m 1 and m 2 , and one segment whose scale is above m 1. The proposal of the ht-index is significant in many aspects of computing the fractality of a geographic feature. The ht-index extends the usage of other heavy-tailed distribution mathematic models to quantify a fractal. Furthermore, as the fractal dimension does not focus on the change of scales among different iteration stages, the ht-index supplements the fractal dimension in terms of measuring the growth of a fractal as the iteration goes. However, the ht-index suffers from distinguishing two fractals with slightly different patterns but the same ht-index value, as the ht-index value is formatted as an integer. To solve the insensitivity of the ht-index, several fractal metrics have been developed, such as CRG and RA [30,31].
CRG is calculated by adding the ratios between consecutive mean values, which is expressed by Equation (4): where i denotes the i th mean values during the head/tail breaks process. Taking the Koch curve in Figure 3a again as a working example, its CRG equals m 2 /m 1 = 0.46/0.2 = 2.38, which supports decimal values, so that the CRG can capture the growth of a fractal more sensitively. RA is an index that examines the extent to which the numeric values are long-tail distributed. To compute RA, we can obtain the ratio between the areas of two polygons A and B based on the rank-size distribution (Figure 3d), as denoted by Equation (5).
where S is the area of a polygon formed by the distribution line.  RA is an index that examines the extent to which the numeric values are long-tail distributed. To compute RA, we can obtain the ratio between the areas of two polygons A and B based on the rank-size distribution (Figure 3d), as denoted by Equation (5).
where is the area of a polygon formed by the distribution line.
In sum, it can be said that the ht-index is the primary metric responsible for the fractal characterization. The CRG and RA indices were, to a large extent, the alternatives to the ht-index, and their major advantage is capturing the sensitivity that the ht-index cannot. More specifically, the CRG sums the differences between every two consecutive hierarchical levels, while RA summarizes the overall distance of the whole population from a heavy-tailed distribution to a normal distribution. Here, we add two more iterations of a Koch curve for clearer illustration. As Table 1 and Figure 4 show, CRG and RA changes nonlinearly as the iteration goes, while the ht-index values appear to be linear. It is also interesting to find that CRG and RA are complementary to each other. Since CRG is an upward trend whose range is [0, ∞ (the higher the CRG, the more scaling hierarchical levels the data), while RA is downward from 1 to 0 (the lower the RA, the more heavy-tailed distributed the data). Together with CRG and RA, we can differentiate fractals with the same scaling hierarchy or the same ht-index value. However, at present, these three metrics are rarely combined for a more comprehensive understanding of a fractal. In this study, we integrate the ht-index, CRG, and RA to implement the fractal analysis of street connectivities for each urban street network. In sum, it can be said that the ht-index is the primary metric responsible for the fractal characterization. The CRG and RA indices were, to a large extent, the alternatives to the ht-index, and their major advantage is capturing the sensitivity that the ht-index cannot. More specifically, the CRG sums the differences between every two consecutive hierarchical levels, while RA summarizes the overall distance of the whole population from a heavy-tailed distribution to a normal distribution. Here, we add two more iterations of a Koch curve for clearer illustration. As Table 1 and Figure 4 show, CRG and RA changes nonlinearly as the iteration goes, while the ht-index values appear to be linear. It is also interesting to find that CRG and RA are complementary to each other. Since CRG is an upward trend whose range is [0, ∞) (the higher the CRG, the more scaling hierarchical levels the data), while RA is downward from 1 to 0 (the lower the RA, the more heavy-tailed distributed the data). Together with CRG and RA, we can differentiate fractals with the same scaling hierarchy or the same ht-index value. However, at present, these three metrics are rarely combined for a more comprehensive understanding of a fractal. In this study, we integrate the ht-index, CRG, and RA to implement the fractal analysis of street connectivities for each urban street network.

Regression Analysis
Prior studies have assessed the growth of the Chinese prefectural-level city system using urban quantities; for example, GDP and population [42]. This study inspects the linkage between Chinese urban quantities and street structure with the identified fractal pattern. To do this, we conduct a multiple regression analysis to unpack relationships between fractal metrics and collected urban indicators including socioeconomic factors (GDP and population) and carbon emissions. Apart from three fractal metrics, we follow the previous study [43] by introducing the city area, number of streets (#Streets), and total length of streets ( ∀ ) as variables to the regression model to estimate comprehensively the explanatory power of street-network configuration measures on cities' socioeconomic status. We take the logarithm of some metrics to scale all variables in a similar magnitude. Therefore, the multiple regression model takes the following form as Equation (6) denotes: We also run two multiple regression models respectively with log(population) and log(CO2 emissions) as the dependent variable. The correlation analysis was conducted based on the Regression tool in Microsoft Excel software.

Regression Analysis
Prior studies have assessed the growth of the Chinese prefectural-level city system using urban quantities; for example, GDP and population [42]. This study inspects the linkage between Chinese urban quantities and street structure with the identified fractal pattern. To do this, we conduct a multiple regression analysis to unpack relationships between fractal metrics and collected urban indicators including socioeconomic factors (GDP and population) and carbon emissions. Apart from three fractal metrics, we follow the previous study [43] by introducing the city area, number of streets (#Streets), and total length of streets (l ∀Streets ) as variables to the regression model to estimate comprehensively the explanatory power of street-network configuration measures on cities' socioeconomic status. We take the logarithm of some metrics to scale all variables in a similar magnitude. Therefore, the multiple regression model takes the following form as Equation (6) denotes: log(GDP) = β 0 + β 1 * log(Area) + β 2 * log(#Streets) + β 3 * log(l ∀Streets ) + β 4 * ht + β 5 * CRG + β 6 * RA.
We also run two multiple regression models respectively with log(population) and log(CO 2 emissions) as the dependent variable. The correlation analysis was conducted based on the Regression tool in Microsoft Excel software.

The Universal Fractal Pattern of 298 Urban Street Networks
The natural streets were generated for all urban street networks; the number of segments and the number of streets ranged from 173,360 and 34,924 (Shanghai) to 175 and 49 (Wuzhishan). The number of segments and the number of streets for each city both obey the power-law distribution ( Figure 5). The number of resulting streets for each city on average is 4.4 times smaller than the number of segments. We applied the head/tail breaks method on the number of streets for all cities and produced an ht-index of 4. From Table 2, we noted that each head percentage was very low (around 30%). This indicates that only a minority of cities possess more than 10,000 streets and most cities have 4000 or fewer streets, given the data sets provided by OSM.
Next, we computed the ht-indexes in terms of street connectivities for each city. Almost all cities had an ht-index value greater than or equal to 3, except for one city (Wuzhishan) whose ht-index value was only 1. The top four cities in terms of ht-index are Shanghai, Dongguan, Hangzhou, and Jinhua, all of which are cities located in highly developed areas of China (for example, around the Yangtze River Delta). Most of the important cities, such as municipalities directly under the central government and provincial capital cities, possess ht-indexes above 5. We can see from Table 3 that more than 80% of cities have ht-indexes above 4, indicating that in each of those cities, the scaling pattern of far more less-connected streets than well-connected ones recurs at least four times. These underlying scaling hierarchical levels at the city scale help us identify the ubiquity of fractal patterns of Chinese urban street networks. We also conducted the power-law detection for the degree of connectivity of each urban street network. In total, there are 199 cities whose street connectivities pass the power-law test. We took a closer look at cities with the same street hierarchical level (that is, the same ht-index value) and found an obvious trend that street networks with a higher ht-index value tend to pass the strict power law examination, as the percentage of "power-law-distributed cities" generally grows when the ht-index increases ( Table 3).

The Universal Fractal Pattern of 298 Urban Street Networks
The natural streets were generated for all urban street networks; the number of segments and the number of streets ranged from 173,360 and 34,924 (Shanghai) to 175 and 49 (Wuzhishan). The number of segments and the number of streets for each city both obey the power-law distribution ( Figure 5). The number of resulting streets for each city on average is 4.4 times smaller than the number of segments. We applied the head/tail breaks method on the number of streets for all cities and produced an ht-index of 4. From Table 2, we noted that each head percentage was very low (around 30%). This indicates that only a minority of cities possess more than 10,000 streets and most cities have 4000 or fewer streets, given the data sets provided by OSM.  Next, we computed the ht-indexes in terms of street connectivities for each city. Almost all cities had an ht-index value greater than or equal to 3, except for one city (Wuzhishan) whose ht-index value was only 1. The top four cities in terms of ht-index are Shanghai, Dongguan, Hangzhou, and Jinhua, all of which are cities located in highly developed areas of China (for example, around the Yangtze River Delta). Most of the important cities, such as municipalities directly under the central government and provincial capital cities, possess ht-indexes above 5. We can see from Table 3 that more than 80% of cities have ht-indexes above 4, indicating that in each of those cities, the scaling pattern of far more less-connected streets than well-connected ones recurs at least four times. These underlying scaling hierarchical levels at the city scale help us identify the ubiquity of fractal patterns of Chinese urban street networks. We also conducted the power-law detection for the degree of connectivity of each urban street network. In total, there are 199 cities whose street connectivities pass the power-law test. We took a closer look at cities with the same street hierarchical level (that is, the same ht-index value) and found an obvious trend that street networks with a higher ht-index value tend to pass the strict power law examination, as the percentage of "power-law-distributed cities" generally grows when the ht-index increases (Table 3). Table 3. Ht-indexes of street connectivities for 298 cities. (Note: the number in brackets indicates the number of cities whose street connectivities pass the power-law test).   To investigate the urban street networks more precisely, we further calculated the CRG and RA indexes for each city. For visualization purposes, we normalized both CRG and RA within the range [0, 1], and used 1 -RA to show the results in order to keep the CRG and RA values in the same direction; that is, the more upper right, the more fractal the city ( Figure 5). Specifically, we can say that the x-axis represents the increment of summed differences of every two consecutive hierarchical levels, while the y-axis shows the closeness of a power-law distribution.
The plot can be explained from both global and local perspectives (Figures 6 and 7, respectively). Globally, the distribution of cities using CRG and RA largely follows the distribution using the ht-index. However, because of the great increase of sensitivity, the differences between cities within the same or across different ht-indexes can be effectively detected. For example, Shanghai, which has both the largest CRG and 1 -RA, was the most fractal city among the top four cities with an ht-index of 9; Beijing, with an ht-index of 7 but a perfect power-law (alpha = 0.58, p-value = 0.26, and 1 -RA = 1), resided uppermost and, more interestingly, in the middle of cities whose ht-indexes were 8.
The plot can be explained from both global and local perspectives (Figures 6 and 7, respectively). Globally, the distribution of cities using CRG and RA largely follows the distribution using the htindex. However, because of the great increase of sensitivity, the differences between cities within the same or across different ht-indexes can be effectively detected. For example, Shanghai, which has both the largest CRG and 1 -RA, was the most fractal city among the top four cities with an ht-index of 9; Beijing, with an ht-index of 7 but a perfect power-law (alpha = 0.58, p-value = 0.26, and 1 -RA = 1), resided uppermost and, more interestingly, in the middle of cities whose ht-indexes were 8. Locally, we can spot that CRG and RA correlated very well at each ht-index, as the cities with the same ht-index value clustered and distributed around the trend line ( Figure 6). For instance, the correlation between CRG and RA can be best depicted using a second-order polynomial function when the ht-index = 5, 6, and 7 (Figure 7b-d), while the CRG linearly correlates with RA when the ht-index = 8 (Figure 7a). The good correlation indicates the positive relationship between the summed differences among hierarchical levels and the nearness of being a power law. As the CRG and RA are not 100% correlated, the trend line, in this connection, works as a "ruler" to discriminate what matters more in a fractal pattern: being more power law (above the trend line) or having more growth between its hierarchical levels (below the trend line). From Figure 7, we found one thing in common for each correlation plot: in most cases, cities with higher economic (based on derived urban GDP) or political status (if a provincial capital city or city under central government) are located at the more upper-right part of the figure area or above the trend line. This implies that the degree of connectivity in street networks in well-developed cities moves toward a power-law distribution, but less so in the less developed cities. Let us zoom into the top 20 cities in China (ranked by GDP) to see the appearance of fractal patterns in the representative cities. As Table 4 shows, the fractal pattern for each of these welldeveloped cities was striking, which was represented by the large ht-index and power law statistics. Note that only Nanjing's street connectivities were non-power-law-distributed, and this is why Locally, we can spot that CRG and RA correlated very well at each ht-index, as the cities with the same ht-index value clustered and distributed around the trend line ( Figure 6). For instance, the correlation between CRG and RA can be best depicted using a second-order polynomial function when the ht-index = 5, 6, and 7 (Figure 7b-d), while the CRG linearly correlates with RA when the ht-index = 8 (Figure 7a). The good correlation indicates the positive relationship between the summed differences among hierarchical levels and the nearness of being a power law. As the CRG and RA are not 100% correlated, the trend line, in this connection, works as a "ruler" to discriminate what matters more in a fractal pattern: being more power law (above the trend line) or having more growth between its hierarchical levels (below the trend line). From Figure 7, we found one thing in common for each correlation plot: in most cases, cities with higher economic (based on derived urban GDP) or political status (if a provincial capital city or city under central government) are located at the more upper-right part of the figure area or above the trend line. This implies that the degree of connectivity in street networks in well-developed cities moves toward a power-law distribution, but less so in the less developed cities.
Let us zoom into the top 20 cities in China (ranked by GDP) to see the appearance of fractal patterns in the representative cities. As Table 4 shows, the fractal pattern for each of these well-developed cities was striking, which was represented by the large ht-index and power law statistics. Note that only Nanjing's street connectivities were non-power-law-distributed, and this is why Nanjing was located below the trend line (Figure 6c). The underlying hierarchies of urban street networks in terms of street connectivities can be effectively visualized by a spectrum of colors from blue to red in Figure 8. Despite the various shapes of different city boundaries, the ubiquitous fractal pattern of far more less-connected streets than well-connected ones between these cities was remarkable. The scarcity of the red lines (the most connected streets) versus the plenty of blue lines (the poorest connected streets), plus other streets in between, formed the urban structure in a fractal fashion.  A spectrum of color is adopted to visualize the city's street structure, with red and blue representing the best-connected and poorest-connected streets, respectively; the other colors (in cyan, yellow, and orange) represent the streets in between).

Correlations between Fractal Metrics and Urban Quantities
The investigation of how the underlying fractal pattern is related to urban quantities was conducted based on the multiple regression analysis, with the results shown in Table 5. The adjust R 2 values are the goodness of fit for the models. Together with traditional street measures, R 2 values were the largest for urban GDP (0.57), followed, in order, by population (0.459) and CO2 emissions (0.285). We also calculated the significance values using a t-test and ensured that the relationships between the independent variables and dependent variable are significant at the 1% significance level. Overall, it can be seen that the ht-index underperforms either CRG or RA with respect to the significance level. We suggest that CRG or RA being more significant than the ht-index is probably due to its increase of sensitivity.
Tables 5 summarizes the effects of different types of street-configuration measures on urban socioeconomic status (log (GDP) and log (Population)) and energy consumption (log (CO2)). The fractality of street networks is positively associated with both GDP and population, while it has a (Note: A spectrum of color is adopted to visualize the city's street structure, with red and blue representing the best-connected and poorest-connected streets, respectively; the other colors (in cyan, yellow, and orange) represent the streets in between).

Correlations between Fractal Metrics and Urban Quantities
The investigation of how the underlying fractal pattern is related to urban quantities was conducted based on the multiple regression analysis, with the results shown in Table 5. The adjust R 2 values are the goodness of fit for the models. Together with traditional street measures, R 2 values were the largest for urban GDP (0.57), followed, in order, by population (0.459) and CO 2 emissions (0.285). We also calculated the significance values using a t-test and ensured that the relationships between the independent variables and dependent variable are significant at the 1% significance level. Overall, it can be seen that the ht-index underperforms either CRG or RA with respect to the significance level. We suggest that CRG or RA being more significant than the ht-index is probably due to its increase of sensitivity.  Table 5 summarizes the effects of different types of street-configuration measures on urban socioeconomic status (log (GDP) and log (Population)) and energy consumption (log (CO 2 )). The fractality of street networks is positively associated with both GDP and population, while it has a slight negative correlation with CO 2 emissions. It is interesting to note that street fractal metrics have a weaker effect than the traditional street measures, such as l ∀Streets and #Streets, while they have a stronger effect than city size. More specifically, a one-unit increase in l ∀Streets is associated with a 1.66-unit increase of urban GDP, which is about 4 times bigger than that linked with a one-unit increase in CRG. Similarly, the increase in population with a one-unit increase in l ∀Streets is five-fold as much as that with one-unit decrease in RA, ceteris paribus. Note that the decrease of RA is consistent with the increase of the ht-index and CRG, as smaller RA values means more striking fractals. For CO 2 emissions, a one-unit increase in street total length and number of streets is correlated with 1.02-unit and 1.57-unit increases, respectively, whereas a one-unit increase in RA (less fractal) is linked with a 0.035-unit increase.

Further Discussion of this Study
The augmented openness and attainability of geospatial big data helps us to analyze and grasp urban forms at both intracity and intercity levels. Street networks largely reflect an overall picture of an urban physical environment. Based on the derived natural streets over 298 cities, we can identify the fractal structure and scaling statistics of China urban street networks and further find the universal fractal pattern among a big variety of urban spaces. In this section, we further discuss the developed methodology and the obtained results and how they contribute to our knowledge in current Chinese urban studies.
Geographic features are fractal in essence, as is the urban space that involves a large collection of geographic features. The invention of the third definition of fractal [6] enables us to perceive the universal fractal pattern at various scales of geographic space. It should be stressed that given the right perspective, almost all geographic features can be seen as fractals. For example, previous work explores the underlying fractality from an individual smooth curve [1] from the very perspective of recursive bends rather than its contained vertices or segments. Similarly, for the collection of linear features-that is, the street network-one must take the perspective of street-street topology instead of segment-segment or junction-junction topology for seeing the fractal. In this study, we employed such perspective and demonstrated successfully the universal fractal pattern of street networks among over 298 cities.
The detected universal fractal pattern of urban street networks can be captured quantitatively by the power-law distribution of street connectivities and the underlying scaling hierarchies. Power-law distribution is widely observed in natural and societal phenomena, and it is usually regarded as an indicator for sustainability [29]. The power-law detection for each city was conducted using the robust maximum likelihood estimation method and helps us identify that about 67% of Chinese urban street networks hold a power law nature. For the remaining 33% of cities, the scaling hierarchical levels derived from head/tail breaks show, as well, the striking fractal patterns of far more less-connected streets than well-connected ones. Thus, the ht-index is superior to power-law metrics with respect to the identification of the universal fractal pattern of Chinese street networks. Besides, the ht-index is good for discriminating urban street networks across different scaling hierarchies; its two alternatives, CRG and RA, help us to differentiate the cities within the same scaling hierarchy because of the increment of sensitivity. In this regard, the integration of three fractal measures can provide us with a more comprehensive and also a more precise picture of fractal cities in China.
The study also attempted to estimate the explanatory power of the fractality of street networks on urban socioeconomic status and fuel consumption. As the urban indicators such as GDP are normally associated with multiple factors of a street network, the study adopted the multiple regression approach by involving other traditional street-network geometric measures in order to make the model closer to the reality and reduce the probability of confounding effects. As a result, the correlation coefficient R 2 value of the multiple regression model, as in the case with GDP (R 2 = 0.57), is higher than that of a single regression model based on one single street geometric property or fractal metric (e.g., R 2 = 0.41 with l ∀Streets only or R 2 = 0.23 with CRG only). Therefore, it can be said that multiple street-network measures are suitable for an exploratory analysis of each urban quantity than that of a single measure.
Our multiple regression analysis further suggests that cities with more strikingly fractal street patterns (i.e., more scaling hierarchical levels of their street networks) tend to have better economic status. The fractal analysis results in Table 4 also confirm that those most fractal cities are top-tier metropolitans in China such as Beijing and Shanghai. If we take a closer look at the coefficient β, it should be noted that the positive effect of l ∀Streets to urban socioeconomic status is stronger than any of three fractal metrics. This indicates that a city's economic development is more directly represented by the geometric details of its street network than by the underlying fractal pattern. However, it does not mean that the fractal metrics are less meaningful. As shown by the regression analysis results with CO 2 emissions, the fractal structure of a street network can contribute to the energy efficiency of a city. The opposite signs of β between fractal metrics and l ∀Streets indicate that street fractalities can perform differently from street geometric properties. In this regard, the fractality of street networks can be a useful supplement to traditional street-network configuration measures.

Conclusions
A city is not only comprised of traffic zones or buildings, but also the collection of individual coherent streets. This article explores the universal fractal pattern of 298 cities in China from the perspective of street networks. The paper applies an integration of three fractal metrics-ht-index, CRG, and RA-for fractal or scaling analysis of street connectivities. The results show that almost all cities bear a fractal structure, as the number of scaling hierarchical levels of the street network for each city is above or equal to 3. We have proven that CRG and RA can greatly compensate for the insensitivity of the ht-index. Besides, the multiple regression analysis results suggest that the fractal structure of the urban street network positively correlates with urban socioeconomic status and negatively correlates with urban energy consumption. Further research will focus on the application of the universality of fractal pattern under the framework of the third definition to the urban planning or policy-decision making by involving other geometric (such as azimuth) and structural properties (such as betweenness) of street networks and more detailed urban indicators.