Big Geodata Reveals Spatial Patterns of Built Environment Stocks Across and Within Cities in China
Zhou Huang
,
Yi Bao
,
Ruichang Mao
,
Han Wang
,
Ganmin Yin
,
Lin Wan
,
Houji Qi
,
Qiaoxuan Li
,
Hongzhao Tang
,
Qiance Liu
,
Linna Li
,
Bailang Yu
,
Qinghua Guo
,
Yu Liu
,
Huadong Guo
,
Gang Liu
The patterns of material accumulation in buildings and infrastructure accompanied by rapid urbanization offer an important, yet hitherto largely missing stock perspective for facilitating urban system engineering and informing urban resources, waste, and climate strategies. However, our existing knowledge on the patterns of built environment stocks across and particularly within cities is limited, largely owing to the lack of sufficient high spatial resolution data. This study leveraged multi-source big geodata, machine learning, and bottom-up stock accounting to characterize the built environment stocks of 50 cities in China at 500 m fine-grained levels. The per capita built environment stock of many cities (261 tonnes per capita on average) is close to that in western cities, despite considerable disparities across cities owing to their varying socioeconomic, geomorphology, and urban form characteristics. This is mainly owing to the construction boom and the building and infrastructure-driven economy of China in the past decades. China’s urban expansion tends to be more “vertical” (with high-rise buildings) than “horizontal” (with expanded road networks). It trades skylines for space, and reflects a concentration-dispersion-concentration pathway for spatialized built environment stocks development within cities in China. These results shed light on future urbanization in developing cities, inform spatial planning, and support circular and low-carbon transitions in cities.
Zhou Huang, Yi Bao, Ruichang Mao, Han Wang, Ganmin Yin, Lin Wan, Houji Qi, Qiaoxuan Li, Hongzhao Tang, Qiance Liu, Linna Li, Bailang Yu, Qinghua Guo, Yu Liu, Huadong Guo, Gang Liu.
Big Geodata Reveals Spatial Patterns of Built Environment Stocks Across and Within Cities in China.
Engineering, 2024, 34 (3) : 143-153 DOI:10.1016/j.eng.2023.05.015
Urbanization is one of the most important global megatrends of the past century [1], [2]. The next three decades will see another 2.5 billion rural residents moving into urban areas—90% of them being from Asia and Africa—and approximately 68% of the global population will live in cities by 2050 [3]. Urbanization represents a process of population concentration [4], expansion of construction land [5], together with the accumulation of materials in buildings [6] and infrastructure [7] (built environment) that defines the physical space of urban activities and provides key services such as shelter and mobility [8]. The construction, maintenance, and demolition of urban built environment stocks result in major sustainability challenges for cities [9], [10], [11], such as resource demand [12], energy use [13], greenhouse gas (GHG) emissions [14], [15], and construction and demolition waste generation [16]. Therefore, understanding the patterns of urban built environment stock development is important to facilitate urban system engineering, inform the circular and low-carbon transition of existing cities, and shed light on future urbanization in the Global South [17].
Previous studies on the patterns of built environment stocks were mostly focused at the regional or national scales, particularly for the temporal dynamics of key construction materials, such as steel [18], [19] and cement [20], and sectors, including buildings [21] and subways [22]. A handful of efforts were made at the urban scale [23], [24], [25]; however, our existing knowledge of the spatial patterns of built environment stocks across and particularly within cities is limited, largely owing to the lack of high spatial resolution data—which is highly data- and labor-intensive [26]. Emerging new types of urban big data, such as point-of-interest data [27], and technologies, such as remote sensing and deep learning [28], [29], offer an opportunity to address such gaps. However, this is not fully captured in the literature.
China is a living laboratory for global urbanization over the past four decades [30]. It has experienced a boom in the development of urban built environment stocks. This benefits economic growth and the well-being of urban residents in China. It also results in significant environmental challenges, including construction and demolition waste generation [16], [31] and GHG emissions [32], [33]. A thorough benchmarking and understanding of built environment stocks across and within cities at different levels of development is essential for China’s endeavor to improve the quality of new urbanization, build zero-waste and eco-cities, and achieve its climate ambition of “peaking before 2030 and neutrality before 2060” [33]. However, existing characterization of built environment stocks for Chinese cities either focus on specific construction materials [24] or sectors [34], [35], [36] without spatial resolution, or are limited to a few cities, such as Beijing [37] and Shanghai [38], or urban areas, such as Tiexi District of Shenyang [39], that cannot support cross-city comparison.
Here, we aimed to address these knowledge gaps by leveraging various urban systems engineering methods involving big geodata, machine learning, and bottom-up stock accounting. We quantified the built environment stocks of 50 Chinese cities and explored their spatial patterns across and within them. Our findings help to inform waste management, urban mining, climate change mitigation, and spatial planning, and the circular and low-carbon transition of Chinese cities, shedding further light on the sustainable urban development of other cities worldwide.
2. Materials and methods
The overall workflow for characterizing the spatialized built environment stocks of the 50 cities in China is shown in Fig. 1. Briefly, multi-source geodata was initially collected and the material composition intensity (MCI) database was established, followed by leveraging bottom-up stock accounting and machine learning approaches to calculate gridded building and infrastructure material stocks of 50 Chinese cities on 500 m fine-grained levels. Finally, a spatial analysis was conducted to reveal the pathway for spatialized urban built environment stock development across and within cities in China.
2.1. Scope and data sources
The built environment stocks considered in this analysis included different types of buildings, such as agricultural, commercial, educational, historical, industrial, mixed, municipal, parking, public, residential, sports, and storage, and transport infrastructure, including roads, railways, and subways. Other types of infrastructure including ports and pipelines contribute very little to the total urban built environment stock and are considered negligible. The year of reference for the calculation was 2018; it was largely based on the availability of multi-source big geodata, and 50 Chinese cities were selected for analysis. They cover all provincial capitals and cities of economic, cultural, and location importance. Together, they account for 50% of the gross domestic product (GDP), 32% of the population, and 15% of the built-up land area of all cities in China, and are deemed representative and sufficient for our comparison.
The important datasets of these 50 selected cities include building footprints—from Baidu, the largest online map portal in China, building age—mainly from real-estate company websites, land use—on a 500 m resolution for five land use categories from the Essential Urban Land Use Categories in China database [40], except for the three sample cities (Beijing, Guangzhou, and Shenzhen) that are based on a 30 m fine-grained level for 12 land use categories, points of interest (POIs)—from Amap, the largest mobile online map platform in China, transport infrastructure—mainly from OpenStreetMap, gridded population—from WorldPop the mainland of China dataset [41], socioeconomic development—mainly from the municipal statistical yearbook, and MCI data collected from various sources (Supplementary material Sections S1.1-1.6).
In particular, a China-specific building MCI database was compiled from various sources, including the bills of quantities, expert interviews, and literature, covering over 2000 sample buildings constructed between 1963 and 2017. They were classified into 12 building typologies (agricultural, commercial, educational, historical, industrial, mixed, municipal, parking, public, residential, sport, and storage). The road and subway MCIs were collected from construction bills provided by several construction companies in China (Supplementary material Section S1.7).
2.2. Bottom-up and spatially refined building stocks of three sample cities
Three cities, Beijing, Guangzhou, and Shenzhen, were selected as the training samples based on data availability. In particular, the 12 building types from 30 m fine-grained land use data were collected in the three cities from the corresponding municipal planning administrations (Supplementary Fig. S2). A bottom-up and spatially refined building stock accounting method was used for these three cities, based on our previous study [37] and shown in Eq. (1). Ten types of construction material were considered: cement, steel, timber, brick, gravel, sand, asphalt, lime, glass, and ceramic.
where MSm,i represents the building stock of material m present in building type i, BFi (measured in m2) is the area of the one-floor building footprint, NF represents the number of building floors, and MCIm,i (kg·m−2) is the composition intensity of material m of type i.
2.3. Machine learning for estimating building stocks for the other 47 cities
Accurate building data with attributes of function, year of construction, and 30 m fine-grained land use data were unavailable for the other 47 cities. Therefore, we leveraged machine learning models to estimate building stocks using the gridded stock values of the three training sample cities aggregated at a 500 m resolution. We combined building attributes and POI attributes to encode each grid with a vector and utilized the random forest model to build the mapping from the grid vector to its material stocks. The model was trained using 80% of the data from Beijing, Guangzhou, and Shenzhen, validated using the remaining 20%, and eventually applied to estimate the building material stock values for each grid of the other 47 cities.
2.4. Transportation material stock calculation
The material stock value was computed for the urban transportation systems in all 50 cities based on the lengths and MCIs of railways, subways, and roads, as shown in Eq. (2). Road MCIs cover five levels: expressways, first-, second-, third-, and fourth-class roads. Railway lines, subway lines, and subway stations were considered for the railway and subway stock estimation.
where MSm,j is the transportation stock of material m in the transportation construction sector j (road, railway, and subway), TLj is the length of transportation (measured in m) in sector j and MCIm,j is the composition intensity of material m in sector j (measured in kg·m−1). S is the subway station, and MCIm,S is the composition intensity of material m in subway station S.
2.5. Spatial statistics for pattern identification
Material stock and population values were used from 500 m × 500 m grids for pattern identification in spatial statistics (Supplementary material Section S3). The grids with the top 1% stock values are regarded as building stock centers. A density-based clustering algorithm (DBSCAN) [42] was used to aggregate the high-stock grids together (Section S3.1). The dispersion index was calculated for each DBSCAN cluster to quantify the spatial dispersion of grids (Section S3.2). The city-level building stock per capita was calculated to understand the role of economic development. The unevenness of the grid-level building stock was determined using the Gini coefficient and the Lorentz curve (Section S3.3). Fitting the grid stock distribution of all 50 cities showed that they conformed to a two-parameter exponential distribution regardless of size, location, and economic development levels (Section S3.4).
3. Results and discussion
3.1. Patterns of urban built environment stocks across cities
The urban built environment stocks of 50 Chinese cities increased to 110 Gt in 2018, which is larger than the total global resource extraction in 2017 (92 Gt) [43]. Buildings (64.3%) and roads (33.4%) dominated the overall construction material stock, whereas other infrastructure only contributed a small share (2.3%). Nonmetallic minerals represented by gravel (51 Gt), cement (26 Gt), sand (17 Gt), and brick (12 Gt) were responsible for 96% of the total types of materials. Steel (1.9 Gt), timber (0.6 Gt), lime (0.5 Gt), and other materials (totaling 1.0 Gt) were used in relatively low quantities (Supplementary Figure Fig. S10).
The total urban built environment stock ranged from 350 Mt in Lhasa to 6771 Mt in Beijing when compared across cities. Stock quantities in 49 out of 50 Chinese cities (2202 Mt on average; except for Lhasa with 350 Mt) and stock densities (4.97 t·m−2 on average) in all 50 Chinese cities were substantially higher than in many western cities—for example, 67 Mt and 0.22 t·m−2 in Odense [12], and 380 Mt and 0.96 t·m−2 in Vienna [44], while per capita stocks (261 tonnes per capita (t·cap−1) in China) were at approximately the same level—329 t·cap−1 in Odense [12], 247 t·cap−1 in Wakayama [12], 210 t·cap−1 in Vienna [44], 209 t·cap−1 in Padua [45], and 272 t·cap−1 in the United Kingdom [46] (Supplementary Table S19). These differences could be explained by the large size and population, but limited built-up area in most Chinese cities [47], together with construction-driven urbanization and real estate-based economic development in the past decades [48].
The urban built environment stocks were unevenly distributed across Chinese cities with varying levels of socioeconomic development (Fig. 2 and Supplementary Figs. S11-14). Cities with large stocks were mostly distributed in the east (35.88%), north (18.58%), and southwest (11.66%), whereas relatively lower amounts were found in cities in the south (9.74%), northeast (9.56%), central (9.11%), and northwest (5.46%) (Supplementary Fig. S15). In particular, the top 10 cities with the largest stocks were all distributed in China’s major urban agglomerations. This includes 3.3 Gt in Hangzhou and 6.8 Gt in Beijing, represented by 25.9 Gt in the Yangtze River Delta and 17.9 Gt in Jing-Jin-Ji Metropolitan Region, which account for 62% of the total. However, cities in the northeast (323 t·cap−1, 3.8 kg·CNY−1) and northwest (346 t·cap−1, 4.0 kg·CNY−1) show significantly higher values than other regions (Supplementary Fig. S16) on a per capita level and per GDP level. This was mainly related to the shrinking population in the northeast [49] and the low population density in the northwest [50], suggesting that material occupancy does not translate into economic growth in these areas [51]. A consideration of the built-up areas shows that cities in the north have the densest stocks (6.1 t·m−2), while the southwest and northwest have lower stock density (both approximately 4.0 t·m−2 on average, Fig. S16). This reflects the geomorphological and socioeconomic characteristics of western China, which has more abundant land and less dense populations than the east [52].
Statistically linear trends between the urban built environment stocks and socioeconomic factors confirm that cities with larger populations and areas and richer cities (in terms of GDP) tend to accumulate more construction materials than smaller and poorer cities (Figs. S12-14). This is also reflected by the fact that the urban built environment stocks tended to increase with an increasing tier rank and urbanization rate after categorizing the 50 cities into six tiers (first, new-first, second, third, fourth, and fifth) using an official classification system based on urban development factors including commercial vitality, transportation convenience, resident activity, lifestyle diversity, and future adaptability (Fig. S11). In this context, megacities with lower stocks per capita (such as 174 t·cap−1 in Chongqing), per square meters built-up area (such as 2.8 t·m−2 in Guangzhou), and per GDP (such as 1.3 kg·CNY−1 in Shenzhen) may have various sustainable paths of stock accumulation and socioeconomic growth that deserves more in-depth analysis to identify leapfrogging opportunities for other yet-to-be developed cities in China and beyond.
The urban form reflecting the physical layouts, structures, and functions of a city is an important driver of the varying levels of material stocks in the urban built environment [53]. Most construction materials were stocked in residential (43%) and industrial (22%) areas, followed by commercial (18%), public (14%) and infrastructure (4%) areas on average across the 50 cities (Fig. 3(a)). However, these proportions vary by city according to their socioeconomic characteristics. For example, Beijing is the capital of China and has the largest share (31%) of stock in public areas (such as education, culture, and healthcare). Meanwhile, Quanzhou and Foshan are two important manufacturing cities in south China that have the largest share of stocks in industrial areas at 52% and 47%, respectively. Furthermore, the variations between material stocks and land areas in different land-use categories clearly reveal the role of urban forms in determining built environment stocks. For example, an average commercial area (often dense and high-rise) accounts for only 3% of land use, but contributes 18% of the total stock, whereas public areas are often sparse and low, account for 49% of land use, but only contribute 14% of the total stock (Figs. 3(a) and (b)).
The building-to-road (BtR) stock ratio in China (5.47 on average for all 50 cities) is notably higher than that in European cities (3.45 in Salford Quays in Manchester [25], 3.13 in Odense city center [12], and 2.94 in Gothenburg [54]) and industrialized counties (1.65 in Japan [55], 1.12 in Germany [56], and 0.91 in Austria [57]). This result of lower road network densities is in line with earlier findings on the city [58] and national [59] levels in China and suggests that China’s urban expansion tends to be more “vertical” than “horizontal.” Further road and infrastructure development, particularly in residential and commercial areas, through better spatial planning or smart design and integration of buildings and roads has become an urgent need to optimize urban services and residents’ well-being (Fig. 3(c)) [60]. Spatially, larger BtR ratio grids were mostly located in city centers, while lower value grids were found in city outskirts; the BtR ratio at the grid level was identified following a log-normal distribution across cities (Supplementary Fig. S21).
3.2. Patterns of spatially refined urban built environment stocks within cities
The gridded building material stocks were clustered in groups of patches using a DBSCAN when presented at a high spatial resolution on the 500 m × 500 m grid level [42] that detects categories based on the closeness of spatial distribution. In contrast, infrastructure material stocks generally follow the distribution of road lines and are spread throughout the city. Furthermore, the spatial patterns of building material stocks in Chinese cities suggest that there are three major phases of urban development: monocentric concentration, multicentric dispersion, and multicentric concentration.
Table 1 and Supplementary Table S23 present the spatial characteristics (number of clusters, dispersion index (DI) of clusters, and Gini index of building material stocks) of the three phases of building material stock growth as the average per capita GDP of cities increases from 70 649, 105 964, and 120 403 CNY, respectively. Fig. 4 shows such spatial patterns in representative case cities (Nanyang, Chongqing, Zhengzhou, and Beijing).
• A city in its relatively early development stage forms very few clusters (exemplified by Nanyang out of 11 cities experiencing a monocentric concentration phase), and a limited number of top building stock grids are compactly distributed in this city with a relatively low dispersion index (DI = 9, Fig. 4(a)).
• An increasing number of grids with large building stocks emerge on the outskirts of urban areas as cities continue to develop and accumulate materials in current grids (as new urban districts and satellite towns). Such multicentric-dispersion phases are featured by spreading clusters (increasing from 2 to 28) and growing dispersion index (increasing from 20 to 221) and can be observed in 22 cities (exemplified by Chongqing and Zhengzhou in Figs. 4(b) and (c), respectively).
• The growing urban built environment stocks gradually help ease communication and mobility in cities by upgrading transportation and telecommunications, attracting more population and businesses, and boosting the urban economy [61]. Accordingly, construction activities and building stock clusters emerge in subsidiary centers that link the central and outskirt clusters. Seventeen cities were identified in this multicentric concentration phase with a closer distance between building stock clusters (DI under 20), including Beijing (Fig. 4(d)).
The 11 cities at the monocentric-concentration phase are Guiyang, Haikou, Handan, Lanzhou, Lhasa, Linyi, Nanchang, Nanyang, Urumqi, Xining, and Zhoukou; the 22 cities at the multicentric-dispersion phase are Baoding, Changsha, Chongqing, Dalian, Foshan, Fuzhou, Harbin, Hangzhou, Luoyang, Nanning, Nantong, Ningbo, Qingdao, Quanzhou, Shijiazhuang, Suzhou, Tangshan, Weifang, Wenzhou, Wuxi, Yinchuan, and Zhengzhou; and the 17 cities at the multicentric-concentration phase are Beijing, Changchun, Chengdu, Guangzhou, Hefei, Hohhot, Jinan, Kunming, Nanjing, Shanghai, Shenyang, Shenzhen, Taiyuan, Tianjin, Wuhan, Xiamen, and Xi’an. The average values are shown for the number of clusters, dispersion index, and Gini index, with the ranges shown in parentheses.
Chinese cities demonstrate an “equilibrium-disequilibrium-equilibrium” pathway of building material stock development corresponding to the “concentration-dispersion-concentration” pattern observed above. This was initially shown in the changes in the Gini index values (average from 0.58, 0.65, and 0.58 for the three phases; Table 1 and Figs. 4(e)-(h)). Moreover, a two-parameter exponential distribution pattern was revealed for the building material stock growth of the 50 cities (exemplified by the probability density functions of Beijing, Suzhou, and Linyi in Fig. 5(a) and detailed for other cities in Supplementary materials Section S3.4 and Appendix A.6).
The distribution of the 50 cities in the four quadrants defined by the location parameter (horizontally with an increasing proportion of low-stock grids) and the scale parameter (vertically with an increasing evenness of stock distribution) of their respective two-parameter exponential distributions are shown in Fig. 5(b). The cluster centers of cities in the monocentric-concentration phase (green star) and in the multicentric-concentration phase (blue star) appear in the second quadrant. This indicated a relatively uniform distribution of building material stocks. In contrast, the cluster center of cities in the multicentric-dispersion phase (orange star) is located in the fourth quadrant. This represents a relatively uneven spatial distribution of building material stocks. The high number of cities with this disequilibrium status (22 out of 50) suggests an urgency for more optimized planning of urban built environment stocks and more coordinated development of urban and rural areas [62], [63]. In this context, cities with higher, but more equalized built environment stocks (such as Changchun with eight clusters, 13 for DI, and 0.46 for GI) can shed light on stock accumulation pathways for other Chinese cities.
The spatially refined building material stocks and gridded population showed a linear correlation with a breakpoint for smaller population grids (R2 = 0.98, on average) and larger ones (R2 = 0.94, on average) in all 50 cities at 500 m resolution. These breakpoints are determined from the continuous piecewise linear function algorithm [64] and mostly vary between 1000 and 3000 for gridded populations—exemplified by Nanyang, Chongqing, Zhengzhou, and Beijing in Figs. 4(i)-(l) and detailed for other cities in Supplementary material Section S3.5). That is, the stock growth rate in grids with a smaller population (Klow = 526) was significantly higher than that in grids with larger populations (Khigh = 98). This implies that more materials were required in the initial stage of urbanization.
3.3. Discussion and implications
Our results reveal that the built environment stocks of many cities in China are close to or higher than those of mature cities in industrialized countries at the per capita level or per area level. This is in line with earlier findings on China’s stock patterns of major construction materials such as cement [20], [65] and aggregates [66] at the national level. Such patterns reflect the construction boom and real estate- and infrastructure-driven urbanization in the past decades in China. Many cities in China are building high-rise residential and non-residential buildings owing to their large population and increasingly limited land area, thus trading skylines for space. This suggests that understanding urban development from a physical stock perspective provides an important and complementary angle for characterizing and informing urbanization that is largely missing in the current literature on urbanization that focuses mostly on population growth [67] and land use change [68].
The spatially refined patterns of urban built environment stocks across and within cities clearly show the role of urban socioeconomic development, such as population and GDP, geomorphology, such as location and land area, and urban form, such as BtR ratio and land use structure, in determining the total volume and the sectoral and spatial distribution of stocks. Therefore, the pace of built environment construction, coordination of buildings and infrastructure development, and tailored approaches for spatial planning and urban resource management should consider the varying stages of urban development in different cities [69], [70]. These spatiotemporally explicit patterns could shed light on future urbanization in western China and other cities in the world to bypass the disequilibrium stage and avoid spatial lock-ins, and provide the public, government, or industry stakeholders with insights into optimized urban spatial planning and urban system engineering towards smart resources, waste, and climate strategies and circular and low-carbon transitions of cities.
Such implications for urban system engineering primarily apply to resource and waste management perspectives. For example, high-resolution mapping of urban built environment stocks allows for an in-depth understanding of urban resource efficiency and forecasting of the quantity, composition, location, and value of future construction and demolition waste generation. Currently, construction and demolition waste in China is mostly dumped or landfilled with only 5% recycled [71]. This challenge will escalate considering the continued urbanization and construction boom in the foreseeable future in China. Understanding built environment stocks with a high spatiotemporal resolution provides a characterization of the urban resource cadaster [12] and enables the circular transition of cities [72]. This includes waste management and urban mining based on spatial and logistics optimization to minimize economic costs and maximize reuse and recycling.
Moreover, understanding the embodied climate impacts of urban built environment stocks facilitates discussions on accounting and mitigating GHG emission during the construction and operation of a city. The carbon replacement value (CRV) [73], [74] concept was adopted to approximate the emissions of constructing a city that would be generated if the existing stock of a city was replaced using current technologies and materials. The overall CRV emissions of the 50 selected cities were estimated to be 32 Gt. This equals 60% of the global GHG emissions in 2019 [75] or 90% of the global CO2 emissions in 2021 [76]. The CRV emissions in Beijing (2.29 Gt or 1.61 t·m−2), Shanghai (2.12 Gt or 2.12 t·m−2), Chengdu (1.33 Gt or 2.57 t·m−2), Suzhou (1.27 Gt or 2.75 t·m−2), and Tianjin (1.26 Gt or 1.25 t·m−2) are among the top five corresponding to the highest urban built environment stocks. This is much larger than those in European (such as 11 Mt in Odense [73]) and Australian (such as 24 Mt in Melbourne [77]) cities. Urban built environment stocks are essential to provide residents with basic services. The CRV emissions in these 50 cities can be used as a benchmark for the climate quota of the other 287 prefecture-level cities in China to reach the same level of services. Relatively low operational emissions and high CRV emissions were observed in two cities that dominate high technology and service-based economies: Shenzhen and Chengdu. In this context, cities with developed economies, upgraded industry structures, low operational emissions, low built environment stocks, and low CRV emissions (for example, Changsha, Quanzhou, Nantong, and Xiamen in the third quadrant of Supplementary material Fig. S22) may be regarded as a model for the low-carbon transition of other small- and medium-sized cities in China and beyond. At least 251 Gt of construction materials equaling 71 Gt of CRV emissions (or approximately seven times the current annual carbon emissions of China [78]) is needed for China’s further urban expansion, assuming that the correlations between urban built environment stocks and socioeconomic parameters (Figs. S12-14) are applicable to the other 287 cities. Therefore, mitigation strategies ranging from technological innovations for emission-intensive construction materials (particularly steel [79] and cement [80]) to improving material efficiency [81], prolonging the building lifetime [31], and optimizing spatial plannin [82], [83] are urgently required to reduce the emissions of constructing a city. All these should be included in future model development and policy formulation from an urban system engineering perspective for the circular and low-carbon transition of cities.
3.4. Method validation and limitations
The proposed machine learning method offers a swift and effective approach for approximating urban material stocks when data incompleteness hinders traditional methods of building stock calculations. Four cross-validation approaches were designed to validate the proposed model and explore the consequent uncertainties. The prediction errors of the model were 0.23%, 3.79%, 0.119%, and 0.02% in the four different validation sets. This indicated that this grid-based method performs well in predicting building material stocks regardless of the city. The composition of building materials differs to some extent in northern and southern China owing to the different demands for winter heating. However, the large number of grid samples, together with our three sampling cities covering the north (Beijing) and south (Guangzhou and Shenzhen) allow the machine learning model to map building attributes to stock values in the sampling cities and make accurate predictions in other cities.
A comparison of our aggregated city-level building stock results for 2018 with a bottom-up accounting study [84] showed that the average differences were below 20% (2.59 Gt vs 3.46 Gt for Chongqing, 4.19 Gt vs 4.75 Gt for Shanghai, and 3.46 Gt vs 2.78 Gt for Tianjin). We believe that the building archetype classification and reference year may contribute to these differences. Furthermore, the building stock results in our study are significantly higher than those of two other studies in a few Chinese cities [27], [38]. This gap is mainly related to the different years of estimation (2018 in this study vs 2010 in previous literature) and the underestimation of nonresidential building stocks in these two studies (see Table S19). Overall, our study provides the first quantification of urban built environment stocks across and within cities on a large scale (50 cities).
The proposed machine learning method helped approximate urban material stocks across and within cities with overall good modelling performance. However, some limitations should be acknowledged. First, the absolute stock results bear unavoidable uncertainties and limitations owing to data gaps. Ideally, MCI data should be specific to each building or transportation infrastructure. We collected over 2000 building samples constructed from 1963-2017 and used the mean value to derive the MCIs for different building typologies. This per square meter indicator is very convenient to scale up, but does not linearly scale with the floor area. For example, a lower material intensity was identified in large-area buildings compared with smaller buildings [85]. More material intensities are required to withstand wind and earthquake loads with increasing building height [86]. Thus, the absolute values of building stock can have an uncertainty of up to ±30%. Second, the data collection process may introduce errors since building coverage may be incomplete, and the vector maps of buildings and infrastructure may differ from reality, leading to uncertainties in the results. The performance of the machine learning model can be compromised when handling materials such as ceramics and glass that have a lower percentage composition. Nevertheless, their effects were relatively negligible owing to their insignificant contributions to the overall quantity. Finally, we used the process-based emission factors from Carbon Emission Accounts and Datasets (CEADs) and the Chinese Life Cycle Database (CLCD) for the CRV calculation using eBalance software. This does not cover the environmental impacts associated with supply chains. Therefore, our CRVs may be underestimated owing to truncation errors [87] in emission accounting. Overall, this study provides the first large-scale quantification of urban built environment stocks, although these limitations should be considered when interpreting the results.
4. Conclusions
Multi-source big geodata, machine learning, and bottom-up stock accounting were leveraged to characterize the built environment stocks of 50 selected cities in China at 500 m fine-grained levels. This large-scale empirical analysis helped to reveal considerable disparities in stocks across Chinese cities owing to their varying socioeconomic, geomorphological, and urban form characteristics. In particular, building material stock development in Chinese cities appeared to follow a two-parameter exponential distribution pattern and a concentration-dispersion-concentration pathway. Our results offer an important, yet hitherto largely missing stock perspective for characterizing and informing urbanization. This informs urban planners and policy makers on spatial planning, and facilitates urban system engineering towards the circular and low-carbon transition of cities. The modeling framework could be extended and validated using more cities to shed more light on future urbanization in China and in other countries.
Acknowledgments
This work is financially supported by the National Natural Science Foundation of China (71991484, 42271471, 72088101, and 41830645), Danish Agency for Higher Education and Science (International Network Project, 0192-00056B), and the Fundamental Research Funds for the Central Universities (Peking University). We acknowledge invaluable comments from Professor Michael F. Goodchild of University of California, Santa Barbara, on an earlier draft of this manuscript.
Compliance with ethics guidelines
Zhou Huang, Yi Bao, Ruichang Mao, Han Wang, Ganmin Yin, Lin Wan, Houji Qi, Qiaoxuan Li, Hongzhao Tang, Qiance Liu, Linna Li, Bailang Yu, Qinghua Guo, Yu Liu, Huadong Guo, and Gang Liu declare that they have no conflict of interest or financial conflicts to disclose.
E.Kalnay, M.Cai. Impact of urbanization and land-use change on climate. Nature, 423 (6939) (2003), pp. 528-531.
[2]
M.R.Montgomery. The urban transformation of the developing world. Science, 319 (5864) (2008), pp. 761-764.
[3]
United Nations Department of Economic and Social Affairs. 68% of the world population projected to live in urban areas by 2050, says UN [Internet]. New York City:United Nations Department of Economic and Social Affairs; 2018 May 16 [cited 2023 May 17]. Available from: https://www.un.org/development/desa/en/news/population/2018-revision-of-world-urbanization-prospects.html
[4]
T.Liu, Y.Qi, G.Cao, H.Liu. Spatial patterns, driving forces, and urbanization effects of China’s internal migration: county-level analysis based on the 2000 and 2010 censuses. J Geogr Sci, 25 (2) (2015), pp. 236-256.
[5]
T.Liu, H.Liu, Y.Qi. Construction land expansion and cultivated land protection in urbanizing China: insights from national land surveys, 1996-2006. Habitat Int, 46 (2015), pp. 13-22.
[6]
X.Zhong, M.Hu, S.Deetman, B.Steubing, H.X.Lin, G.A.Hernandez, et al. Global greenhouse gas emissions from residential and commercial building materials and mitigation strategies to 2060. Nat Commun, 12 (1) (2021), Article 6126.
[7]
M.Jiang, P.Behrens, T.Wang, Z.Tang, Y.Yu, D.Chen, et al. Provincial and sector-level material footprints in China. Proc Natl Acad Sci USA, 116 (52) (2019), pp. 26484-26490.
[8]
F.Creutzig, L.Niamir, X.Bai, M.Callaghan, J.Cullen, J.Díaz-José, et al. Demand-side solutions to climate change mitigation consistent with high levels of well-being. Nat Clim Chang, 12 (1) (2022), pp. 36-46.
[9]
S.Chen, B.Chen, K.Feng, Z.Liu, N.Fromer, X.Tan, et al. Physical and virtual carbon metabolism of global cities. Nat Commun, 11 (1) (2020), Article 182.
[10]
X.Zhong, S.Deetman, A.Tukker, P.Behrens. Increasing material efficiencies of buildings to address the global sand crisis. Nat Sustain, 5 (5) (2022), pp. 389-392.
[11]
Intergovernmental Panel on Climate Change (IPCC). Climate change 2014—mitigation of climate change. Working group III contribiution to the fifth assessment report of the Intergovernmental Panel on Climate Change. New York City: Cambridge University Press; 2014.
[12]
M.Lanau, G.Liu. Developing an urban resource cadaster for circular economy: a case of Odense. Denmark Environ Sci Technol, 54 (7) (2020), pp. 4675-4685.
[13]
C.A.Kennedy, I.Stewart, A.Facchini, I.Cersosimo, R.Mele, B.Chen, et al. Energy and material flows of megacities. Proc Natl Acad Sci USA, 112 (19) (2015), pp. 5985-5990.
[14]
World Green Building Council. New report:the building and construction sector can reach net zero carbon emissions by 2050. Report. London: World Green Building Council; 2019 Sep 23.
H.Nagendra, X.Bai, E.S.Brondizio, S.Lwasa. The urban south and the predicament of global sustainability. Nat Sustain, 1 (7) (2018), pp. 341-349.
[18]
T.Wang, D.B.Müller, T.E.Graedel. Forging the anthropogenic iron cycle. Environ Sci Technol, 41 (14) (2007), pp. 5120-5129.
[19]
D.B.Müller, T.Wang, B.Duval, T.E.Graedel. Exploring the engine of anthropogenic iron cycles. Proc Natl Acad Sci USA, 103 (44) (2006), pp. 16111-16116.
[20]
Z.Cao, L.Shen, A.N.Løvik, D.B.Müller, G.Liu. Elaborating the history of our cementing societies: an in-use stock perspective. Environ Sci Technol, 51 (19) (2017), pp. 11468-11475.
[21]
N.Heeren, S.Hellweg. Tracking construction material over space and time: prospective and geo-referenced modeling of building stocks and construction material flows. J Ind Ecol, 23 (1) (2019), pp. 253-267.
[22]
R.Mao, Y.Bao, H.Duan, G.Liu. Global urban subway development, construction material stocks, and embodied carbon emissions. Humanit Soc Sci Commun, 8 (1) (2021), p. 83.
[23]
Q.Liu, Z.Cao, X.Liu, L.Liu, T.Dai, J.Han, et al. Product and metal stocks accumulation of China’s megacities: patterns, drivers, and implications. Environ Sci Technol, 53 (8) (2019), pp. 4128-4139.
[24]
C.Huang, J.Han, W.Q.Chen. Changing patterns and determinants of infrastructures’ material stocks in Chinese cities. Resour Conserv Recycling, 123 (2017), pp. 47-53.
[25]
H.Tanikawa, S.Hashimoto. Urban stock over time: spatial material stock analysis using 4D-GIS. Build Res Inform, 37 (5-6) (2009), pp. 483-502.
[26]
M.Lanau, G.Liu, U.Kral, D.Wiedenhofer, E.Keijzer, C.Yu, et al. Taking stock of built environment stock studies: progress and prospects. Environ Sci Technol, 53 (15) (2019), pp. 8499-8515.
[27]
J.Guo, A.Miatto, F.Shi, H.Tanikawa. Spatially explicit material stock analysis of buildings in eastern China metropoles. Resour Conserv Recycling, 146 (2019), pp. 45-54.
Y.Bao, Z.Huang, H.Wang, G.Yin, X.Zhou, Y.Gao. High‐resolution quantification of building stock using multi‐source remote sensing imagery and deep learning. J Ind Ecol, 27 (1) (2023), pp. 350-361.
Z.Cao, G.Liu, H.Duan, F.Xi, G.Liu, W.Yang. Unravelling the mystery of Chinese building lifetime: a calibration and verification based on dynamic material flow analysis. Appl Energy, 238 (2019), pp. 442-452.
[32]
L.Hong, N.Zhou, W.Feng, N.Khanna, D.Fridley, Y.Zhao, et al. Building stock dynamics and its impacts on materials and energy demand in China. Energy Policy, 94 (2016), pp. 47-55.
[33]
H.Wang, X.Lu, Y.Deng, Y.Sun, C.P.Nielsen, Y.Liu, et al. China’s CO2 peak before 2030 implied from characteristics and growth of cities. Nat Sustain, 2 (8) (2019), pp. 748-754.
[34]
R.Mao, H.Duan, D.Dong, J.Zuo, Q.Song, G.Liu, et al. Quantification of carbon footprint of urban roads via life cycle assessment: case study of a megacity—Shenzhen. China J Clean Prod, 166 (2017), pp. 40-48.
[35]
T.Wang, J.Zhou, Y.Yue, J.Yang, S.Hashimoto. Weight under steel wheels: material stock and flow analysis of high-speed rail in China. J Ind Ecol, 20 (6) (2016), pp. 1349-1359.
[36]
Z.Guo, D.Hu, F.Zhang, G.Huang, Q.Xiao. An integrated material metabolism model for stocks of urban road system in Beijing. China Sci Total Environ, 470-471 (2014), pp. 883-894.
[37]
R.Mao, Y.Bao, Z.Huang, Q.Liu, G.Liu. High-resolution mapping of the urban built environment stocks in Beijing. Environ Sci Technol, 54 (9) (2020), pp. 5345-5355.
[38]
J.Han, W.Q.Chen, L.Zhang, G.Liu. Uncovering the spatiotemporal dynamics of urban infrastructure development: a high spatial resolution material stock and flow analysis. Environ Sci Technol, 52 (21) (2018), pp. 12122-12132.
[39]
J.Guo, T.Fishman, Y.Wang, A.Miatto, W.Wuyts, L.Zheng, et al. Urban development and sustainability challenges chronicled by a century of construction material flows and stocks in Tiexi. China J Ind Ecol, 25 (1) (2021), pp. 162-175.
[40]
P.Gong, B.Chen, X.Li, H.Liu, J.Wang, Y.Bai, et al. Mapping essential urban land use categories in China (EULUC-China): preliminary results for 2018. Sci Bull, 65 (3) (2020), pp. 182-187.
[41]
WorldPop. Global high resolution population denominators project—funded by the bill and melinda gates foundation [Internet]. Otara: WorldPop; 2020 Feb 1 [cited 2023 May 17]. Available from: https://www.worldpop.org/geodata/summary?id=24924
[42]
M.Ester, H.P.Kriegel, J.Sander, X.Xu. A density-based algorithm for discovering clusters in large spatial databases with noise. Proceedings of the Second International Conference on Knowledge Discovery and Data Mining (KDD-96); 1996 Aug 2 -4, AAAI Press, Oregon, Portland. Cambridge (1996), pp. 226-231.
[43]
United Nations Environment Programme (UNEP). Global resources outlook 2019:natural resources for the future we want. Report. Paris: International Resource Panel (IRP); 2019.
[44]
F.Kleemann, J.Lederer, H.Rechberger, J.Fellner. GIS-based analysis of Vienna’s material stock in buildings. J Ind Ecol, 21 (2) (2017), pp. 368-380.
[45]
A.Miatto, H.Schandl, L.Forlin, F.Ronzani, P.Borin, A.Giordano, et al. A spatial analysis of material stock accumulation and demolition waste potential of buildings: a case study of Padua. Resour Conserv Recycling, 142 (2019), pp. 245-256.
[46]
J.Streeck, D.Wiedenhofer, F.Krausmann, H.Haberl. Stock-flow relations in the socio-economic metabolism of the United Kingdom 1800-2017. Resour Conserv Recycling, 161 (2020), Article 104960.
[47]
X.Bai, J.Chen, P.Shi. Landscape urbanization and economic growth in China: positive feedbacks and sustainability dilemmas. Environ Sci Technol, 46 (1) (2012), pp. 132-139.
[48]
Z.Liu, Z.Deng, G.He, H.Wang, X.Zhang, J.Lin, et al. Challenges and opportunities for carbon neutrality in China. Nat Rev Earth Environ, 3 (2) (2021), pp. 141-155.
[49]
X.Meng, Y.Long. Shrinking cities in China: evidence from the latest two population censuses 2010-2020. Environ Plann A, 54 (3) (2022), pp. 449-453.
[50]
W.Qi, S.Liu, M.Zhao, Z.Liu. China’s different spatial patterns of population growth based on the “Hu Line”. J Geogr Sci, 26 (11) (2016), pp. 1611-1625.
[51]
J.Tan, K.Lo, F.Qiu, W.Liu, J.Li, P.Zhang. Regional economic resilience: resistance and recoverability of resource-based cities during economic crises in northeast China. Sustainability, 9 (12) (2017), p. 2136.
[52]
B.Gao, Q.Huang, C.He, Z.Sun, D.Zhang. How does sprawl differ across cities in China? A multi-scale investigation using nighttime light and census data. Landsc Urban Plan, 148 (2016), pp. 89-98.
[53]
J.R.VandeWeghe, C.Kennedy. A spatial analysis of residential greenhouse gas emissions in the Toronto Census Metropolitan Area. J Ind Ecol, 11 (2) (2007), pp. 133-144.
[54]
P.Gontia, L.Thuvander, B.Ebrahimi, V.Vinas, L.Rosado, H.Wallbaum. Spatial analysis of urban material stock with clustering algorithms: a northern European case study. J Ind Ecol, 23 (6) (2019), pp. 1328-1343.
[55]
H.Tanikawa, T.Fishman, K.Okuoka, K.Sugimoto. The weight of society over time and space: a comprehensive account of the construction material stock of Japan, 1945-2010. J Ind Ecol, 19 (5) (2015), pp. 778-791.
[56]
G.Schiller, F.Müller, R.Ortlepp. Mapping the anthropogenic stock in Germany: metabolic evidence for a circular economy. Resour Conserv Recycling, 123 (2017), pp. 93-107.
[57]
DaxbeckH, BuschmannH, NeumayerS, BrandtB. Methodology for mapping of physical stocks. Sixth framework programme priority. Report. Austria: Resource Management Agency (RMA); 2009 Oct 6. Contract No.: 044409.
[58]
G.Zhao, X.Zheng, Z.Yuan, L.Zhang. Spatial and temporal characteristics of road networks and urban expansion. Land, 6 (2) (2017), p. 30.
[59]
J.Hong, Z.Chu, Q.Wang. Transport infrastructure and regional economic growth: evidence from China. Transportation, 38 (5) (2011), pp. 737-752.
[60]
R.F.M.Ameen, M.Mourshed, H.Li. A critical review of environmental assessment tools for sustainable urban design. Environ Impact Assess Rev, 55 (2015), pp. 110-125.
[61]
S.Thacker, D.Adshead, M.Fay, S.Hallegatte, M.Harvey, H.Meller, et al. Infrastructure for sustainable development. Nat Sustain, 2 (4) (2019), pp. 324-331.
[62]
Y.Song. Rising Chinese regional income inequality: the role of fiscal decentralization. China Econ Rev, 27 (2013), pp. 294-309.
[63]
LiKQ. Report on the work of the government [Internet]. Beiing:State Council of the People’s Republic of China; 2019 May 16 [cited 2023 May 17]. Available from: https://english.www.gov.cn/premier/speeches/2019/03/16/content_281476565265580.htm
[64]
JekelCF, VenterG. pwlf: a Python library for fitting 1D continuous piecewise linear functions [Internet]. Online: GetHub, Inc.;2019 Feb 6 [cited 2023 May 17]. Available from: https://github.com/cjekel/piecewise_linear_fit_py
[65]
Z.Cao, R.J.Myers, R.C.Lupton, H.Duan, R.Sacchi, N.Zhou, et al. The sponge effect and carbon emission mitigation potentials of the global cement cycle. Nat Commun, 11 (1) (2020), p. 3777.
[66]
Z.Ren, M.Jiang, D.Chen, Y.Yu, F.Li, M.Xu, et al. Stocks and flows of sand, gravel, and crushed stone in China (1978-2018): evidence of the peaking and structural transformation of supply and demand. Resour Conserv Recycling, 180 (2022), Article 106173.
[67]
H.Buhaug, H.Urdal. An urbanization bomb? Population growth and social disorder in cities. Glob Environ Change, 23 (1) (2013), pp. 1-10.
[68]
L.Lai, X.Huang, H.Yang, X.Chuai, M.Zhang, T.Zhong, et al. Carbon emissions from land-use change and management in China between 1990 and 2010. Sci Adv, 2 (11) (2016), Article e1601063.
[69]
T.Elmqvist, E.Andersson, N.Frantzeskaki, T.McPhearson, P.Olsson, O.Gaffney, et al. Sustainability and resilience for transformation in the urban century. Nat Sustain, 2 (4) (2019), pp. 267-273.
[70]
R.York, E.A.Rosa, T.Dietz. STIRPAT, IPAT and ImPACT: analytic tools for unpacking the driving forces of environmental impacts. Ecol Econ, 46 (3) (2003), pp. 351-365.
[71]
B.Huang, X.Wang, H.Kua, Y.Geng, R.Bleischwitz, J.Ren. Construction and demolition waste management in China through the 3R principle. Resour Conserv Recycling, 129 (2018), pp. 36-44.
[72]
A.Ajayebi, P.Hopkinson, K.Zhou, D.Lam, H.M.Chen, Y.Wang. Spatiotemporal model to quantify stocks of building structural products for a prospective circular economy. Resour Conserv Recycling, 162 (2020), 105026.
[73]
M.Lanau, L.Herbert, G.Liu. Extending urban stocks and flows analysis to urban greenhouse gas emission accounting: a case of Odense. Denmark J Ind Ecol, 25 (4) (2021), pp. 961-978.
[74]
D.B.Müller, G.Liu, A.N.Løvik, R.Modaresi, S.Pauliuk, F.S.Steinhoff, et al. Carbon emissions of infrastructure development. Environ Sci Technol, 47 (20) (2013), pp. 11739-11746.
[75]
UNEP Copenhagen Climate Centre. Emissions gap report 2020. Report. Copenhagen: UNEP Copenhagen Climate Centre; 2020 Dec 9.
[76]
Z.Liu, Z.Deng, S.J.Davis, C.Giron, P.Ciais. Monitoring global carbon emissions in 2021. Nat Rev Earth Environ, 3 (4) (2022), pp. 217-219.
[77]
A.Stephan, A.Athanassiadis. Quantifying and mapping embodied environmental requirements of urban building stocks. Build Environ, 114 (2017), pp. 187-202.
[78]
Z.Mi, Y.M.Wei, B.Wang, J.Meng, Z.Liu, Y.Shan, et al. Socioeconomic impact assessment of China’s CO2 emissions peak prior to 2030. J Clean Prod, 142 (2017), pp. 2227-2236.
[79]
B.Yu, X.Li, Y.Qiao, L.Shi. Low-carbon transition of iron and steel industry in China: carbon intensity, economic growth and policy intervention. J Environ Sci, 28 (2015), pp. 137-147.
[80]
T.Gao, L.Shen, M.Shen, L.Liu, F.Chen, L.Gao. Evolution and projection of CO2 emissions for China’s cement industry from 1980 to 2020. Renew Sustain Energy Rev, 74 (2017), pp. 522-537.
[81]
C.F.Dunant, M.P.Drewniok, M.Sansom, S.Corbey, J.M.Cullen, J.M.Allwood. Options to make steel reuse profitable: an analysis of cost and risk distribution across the UK construction value chain. J Clean Prod, 183 (2018), pp. 102-111.
[82]
S.H.Wang, S.L.Huang, P.J.Huang. Can spatial planning really mitigate carbon dioxide emissions in urban areas? A case study in Taipei. Taiwan Landsc Urban Plan, 169 (2018), pp. 22-36.
[83]
Y.Bao, Z.Huang, L.Li, H.Wang, J.Lin, G.Liu. Evaluating the human use efficiency of urban built environment and their coordinated development in a spatially refined manner. Resour Conserv Recycling, 189 (2023), p. 106723.
[84]
L.Song, J.Han, N.Li, Y.Huang, M.Hao, M.Dai, et al. China material stocks and flows account for 1978-2018. Sci Data, 8 (1) (2021), p. 303.
[85]
A.Stephan, R.H.Crawford. The relationship between house size and life cycle energy demand: implications for energy efficiency regulations for buildings. Energy, 116 (2016), pp. 1158-1171.
[86]
J.Helal, A.Stephan, R.H.Crawford. The influence of structural design methods on the embodied greenhouse gas emissions of structural systems for tall buildings. Structures, 24 (2020), pp. 650-665.
[87]
R.H.Crawford, A.Stephan, F.Prideaux. The EPiC database: hybrid embodied environmental flow coefficients for construction materials. Resour Conserv Recycling, 180 (2022), 106058.