Research on the Geographical Pattern, Evolution Model, and Driving Mechanism of Carbon Emission Density from Urban Industrial Land in the Yangtze River Economic Belt of China

: To achieve the goals of “carbon peaking and carbon neutrality”, this paper puts forward the connotation and measurement method for the carbon emission intensity of urban industrial land and conducts an empirical study with the Yangtze River Economic Belt (YREB) as an example. We defined the carbon intensity of urban industrial land as the industrial carbon emissions per unit area of land, which is a spatial mapping of urban industrial economic development and carbon spillover and a key indicator for urban and territorial spatial planning oriented towards the “dual carbon” goal. Findings: The carbon emission density of industrial land in the YREB varied greatly between cities and exhibited significant positive spatial autocorrelation. In addition, the geographical pattern and spatio-temporal evolution model of the urban industrial land carbon emission density had a very complex driving mechanism, and different factors had significant synergistic effects. Therefore, it is suggested that while striving towards the goal of “dual carbon”, the government should incorporate the carbon emission density indicator of urban industrial land into the urban and territorial spatial planning system, and based on the threshold of the medium suitable density, they should design differentiated management policies according to concrete urban policies and encourage co-operation among cities to jointly promote carbon emission management of urban industrial land. In policy design, emphasis should also be placed on highlighting the interactive effects of foreign direct investment, fiscal expenditure, and the number of patent authorizations as well as constructing a combination of policies centered around them to better leverage the impacts of globalization, government intervention, and innovation.


Research Background
More and more governments, scholars, and businesses are now beginning to take notice of global warming caused by greenhouse gas emissions.Land-use change, as a principal element influencing global climate change and the carbon cycle, alters natural carbon processes at the surface and also affects the rate of the carbon cycle by changing the intensity of human energy consumption.The industrial sector is one of the most important sources of carbon emissions, and there is a close link between the spatial layout of industrial land and the types of industries and carbon emission density [1].Therefore, it is a growing concern that strict control and the adjustment of land-use structures, especially strengthening the management of carbon emissions from industrial land use, can effectively mitigate global climate change.The World Cities Report 2022: Envisaging the Future of Cities, published by UN Habitat, stresses managing urban density in the context of multiple urban challenges as a key challenge for sustainable urban construction in the future.From the history of urban construction in China and abroad, it has been found that density has not only been applied to the measurement of urban form, but has also gradually evolved into an important starting point for dealing with urban problems, controlling urban form, and regulating the relationship between man and land in different periods [2,3].Therefore, density research has outstanding theoretical and practical significance for urban construction management, and it is also the core proposition related to the sustainable construction of cities in the future.
Urban planning and land spatial planning, as the government's public intervention on the built environment, aim to ensure space supply, protect public interests, and overcome market failure.For urban and land spatial planning, density is an important indicator for measuring the coordination of land resource allocation, urbanization level, and spatial development quality [4,5].Carbon emission density refers to the carbon emissions per unit land.As a new key indicator in urban and land spatial planning, the carbon intensity reflects the spatial evolution of cities and their regions in terms of land-use intensity, the level of economic development, and urban construction, serving as an index that represents how much of the city's "dual carbon" goal is achieved.The change in the carbon emission density of urban industrial land, in essence, is the result of the continuous inflow and outflow of urban industrial resources and enterprise clusters, employed population and human capital, industry and innovation capital, and other factors.It is accompanied by changes in the population density, employment density, economic density, and building density in urban spaces.Therefore, an empirical study on the geographical pattern, evolutionary pattern, and driving mechanism of carbon density of industrial land in urban agglomerations based on the spatial measurement model will provide key evidence and information for urban planning and territorial spatial planning to serve the goal of "carbon peaking and carbon neutrality", with a great theoretical and practical value.

Literature Review
The carbon emission density of industrial land is an emerging tool to measure the relationship between the scale of industrial carbon emissions and the area of industrial land occupied.Therefore, the literature review needs to take into account the research results on carbon emissions of land use and industrial development and use the ideas and methodological innovations from urban and territorial spatial densities, such as population density and building density.

Land Use and Carbon Emissions
Carbon emissions come from a variety of sources, with carbon emissions from land use ranking second after fossil fuel consumption.Economic development leads to changes in land use and has a significant impact on the carbon cycle in nature.Therefore, an indepth analysis of the relationship between land-use carbon emissions and economic growth is of great significance to the scientific development of energy conservation and emission reduction programs, and it is also the key to promoting regional carbon neutrality as well as the sustainable development of the society [6].Since the 1970s, an increasing number of scholars have been involved in research on land-use carbon emissions, summarized into the following areas and features: research related to land-use carbon sources and sinks from the perspective of the carbon cycle in terrestrial ecosystems, mainly dealing with the accounting of land-use carbon stocks and fluxes and their influencing factors [7,8]; correlation analysis of the relationship between carbon emissions and economic growth or factors of production based on land carbon emissions using econometric methods such as the gray correlation model [9], the coupled coordination degree model [10], and the decoupling model [11]; and research related to low-carbon development policies for land use under the "dual carbon" goal, that is, to portray the evolution models of land-use carbon emissions through methods such as spatial autocorrelation and cold-hot spot analysis and to further propose their policy implications.For example, Sidorczuk-Pietraszko [12] studied the spatial differences in household carbon intensities between different regions in Poland by means of exponential decomposition analysis (IDA); Rong [13] also found that there is an obvious spatial correlation between land-use carbon emissions in Chinese provinces, and thus suggested that China should strengthen the interprovincial cooperation in order to achieve an overall emission reduction across the country.

Carbon Emissions from Industries
The industrial sector, as an important source of carbon emissions from human activities, shoulders the important responsibility of energy conservation and emission reduction, and it faces the important and urgent task of how to control and reduce carbon emissions while pursuing the high-quality development of industry and promote the mutual harmonization of economic development and resources and the environment.As a result, the study of industrial carbon emissions has become a hot topic.At present, the research on industry-related carbon emissions mainly focuses on the following three areas: The first is the research on carbon emission measurement and energy conservation and emission reduction in specific industrial sectors or industries, including the communications and transportation industry [14,15], tourism [16], agriculture [17], services [18], and construction [19], with topics involving carbon emission measurement and prediction methodology, influencing factors, spatio-temporal analysis, and the link between carbon emissions and economic development.The second is the research on the influencing factors of industrial carbon emissions, including the energy structure, industrial structure, economic output, population size, energy intensity, urbanization rate, technological progress, and foreign trade and is mostly based on index decomposition analysis, structural decomposition analysis, IPAT or STIRPAT models, and other econometrics and regression analysis methods [20,21].The third is the research on industrial carbon transfer, carbon leakage, implied carbon emissions, and the differential responsibility for carbon reduction.With the evolution of regional economic development patterns, inter-regional industrial transfers have become more frequent, producing different environmental effects on the transferring places.In this context, some scholars have carried out research on the scale calculation, influencing factors, and mechanism of industrial carbon transfer and analyzed the relevance of industrial carbon emissions and its transmission network characteristics along the industrial chain [22].Research on the differential responsibility for carbon emission reductions has mainly centered on the "polluter pays" principle, from which the principles of "responsibility within national territories" and "producer burden of pollution" have been derived [23].Depending on the principles of "producer burden of pollution" and "consumer burden of pollution", scholars have put forward new principles such as "common sharing between importing and exporting countries" [12], "common sharing between producers and consumers" [13], and "common sharing between upstream and downstream of the industrial chain" according to different research objects [24,25].

Urban and Territorial Spatial Density
In human civilizations, cities are spatial carriers where factors such as population, resources, technology, and information are highly concentrated, and density, as a visual representation of the degree of spatial concentration of these factors, has become a key indicator for understanding the efficiency of the distribution of the various factors in cities and their spatial formations in a complex mega-system.Since the emergence of urban planning and territorial spatial planning, attention has been paid to the study of urban density, mainly on the following areas: The first is the research on the urban density measurement indicators and evaluation mechanism.Urban density research originated from population density analysis [26] and then expanded to building density [27], economic density [28], environmental density [29], employment density [30], land density [31], and road network density [32].The second is the research on the spatio-temporal differentiation of urban density.Leveraging the exploratory spatial data analysis, Theil coefficient, coefficient of variation, Markov chain, kernel density estimation, and other models, it explores and empirically analyzes the laws of spatial differentiation, dynamic evolution, and spatio-temporal evolution of urban density at different scales [33,34].The third is research on the interaction between urban density and other factors.It analyzes the impact of urban density on the economic efficiency, technological innovation, and land-use patterns and intensities.For example, Hong [35] found that when the urban population size is less than 1 million, increasing the urban density promotes carbon emission reduction; on the contrary, it will lead to more carbon emissions.Lee [36] and Li [37] analyzed the impact of urban morphology-particularly population density and building density-on household carbon emissions and household electricity consumption.It also explores the driving factors and mechanisms of the dynamic evolution and spatio-temporal evolution of urban density, revealing the key influencing factors and their pathways of action.For example, using geographically weighted regression models and Baidu heat maps, Zhang [38] found, in a case study of Wuhan, that residential workplaces, education levels, and the quality of urban landscapes are key indicators of real-time urban human activity density.Of note, the carbon emission density shares a similar concept with the population density, characterized by spatial clustering and mobility.Therefore, we can define the carbon intensity of industrial land as the ratio of the carbon emissions from the production of industrial enterprises, as well as the carbon emissions from the liquidity of the business and product transportation, to the area of industrial land.The measurement of carbon emissions and industrial land use is not limited to narrowly defined production processes and carriers but includes supporting links such as transportation and productive services [39].

Research Gap
In general, scholars have conducted in-depth analyses on land use and carbon emissions, industrial carbon emissions, and urban density, providing valuable references for this study.However, there are still some vulnerable spots, which are mainly as follows: First, industrial development, land use, and carbon emissions are closely related, but existing studies generally separate the three to specifically study carbon emissions from industrial development or land use without incorporating them into a unified research framework.Simply calculating the carbon emissions of the industry without taking other factors into account often fails to provide a comprehensive understanding of the industry's carbon emissions, leading to the unsatisfactory implementation of carbon emission management policies for industrial land.With the factor of "land area occupied by industrial activity space" taken, the implementation of industrial carbon emissions to the land, focusing on the industrial carbon emission scale to the industrial carbon emission density, so as to organically integrate land use, industrial development and carbon emission management will be of great value to the policy design of energy conservation and emission reduction, industrial development, and land use under the "dual carbon" goal.
Second, most of the existing studies regard the study area as a separate entity and rarely discuss the spatial effect of carbon emission density from the perspective of geography and spatial interaction.No spatial econometric model has been selected for the regression analysis, which affects the accuracy of the results and conclusions.Cities and regions constitute an open system where information and resources flow and are exchanged constantly between cities or regions, and spatial elements tend to overflow.The first law of geography states that things that are closer in space are more dependent on each other; that is, the spatial autocorrelation of the attribute values is more significant.Wu [40] and Dong [41] have conducted exploratory research on carbon emission density, but it should be noted there are weaknesses in their studies.The former focused on the carbon intensity of agricultural land rather than industrial land, and its analytical results and conclusions are more applicable to rural planning than to urban and land spatial planning.The latter used the extended Kaya constant equation and the log-mean Divisia index model rather than a spatial econometric model and failed to capture the impact of spatial effects.

Research Question
The objective of this paper is to discuss the following issues: First, in the dimensions of urban planning and territorial spatial planning disciplines, what is the connotation and value of the concept of carbon emission density of industrial land, and how can it be measured quantitatively?Second, what are the regular characteristics of the geographical pattern and spatio-temporal evolution model of carbon emission density of industrial land?Third, what are the factors and how do they affect the carbon emission intensity of industrial land?This article takes the YREB as an example to empirically analyze the above issues through spatial econometric models, aiming to provide key information and a basis for dual-carbon-goal-oriented urban industrial planning and national spatial planning.

Study Area
The YREB is home to 11 provincial administrative regions, including Jiangsu, Zhejiang, Anhui, Jiangxi, Hubei, Hunan, Sichuan, Yunnan, Guizhou, Shanghai, and Chongqing, and the administrative units of this study area include 110 cities.In the study area, 16 autonomous prefectural regions and 3 provincial counties were excluded due to a significant amount of missing data or poor data quality, the former being Enshi in Hubei, Xiangxi in Hunan, Aba in Sichuan, Ganzi in Sichuan, Liangshan in Sichuan, Qianxinan in Guizhou, Qiandongnan in Guizhou, Qiannan in Guizhou, Chuxiong in Yunnan, Honghe in Yunnan, Wenshan in Yunnan, Xishuangbanna in Yunnan, Dali in Yunnan, Dehong in Yunnan, Nujiang in Yunnan, and Diqing in Yunnan and the latter being Xiantao, Qianjiang, and Tianmen in Hubei (Figure 1).The YREB is a major strategic space for regional development at the national level, alongside the Beijing-Tianjin-Hebei region, the Guangdong-Hong Kong-Macao Greater Bay Area, and the Yellow River Basin, with a good industrial base and urban construction.The YREB plays a crucial role in the construction of China's new plan and is highly representative in both global industrial clusters and urban agglomerations.The Chinese central government put forward the Outline of the YREB Development Plan in 2016, as well as the Implementation Plan for the Development of the YREB under the 14th Five-Year Plan and the Opinions on Several Policy Measures to Further Promote the High-Quality Development of the YREB in 2021 and 2023.All these documents call for industry transformation and green development in the YREB, the synergistic promotion of ecological environmental protection, and sustainable economic development, so as to build a beautiful China model, where human beings and nature coexist in harmony.The YREB accounts for up to 40% of the country's population size and economic output and more than one-third of its carbon dioxide emissions.Therefore, strengthening the carbon emission management of industrial land for the high-quality green development of the YREB has become a key issue that needs to be solved urgently.

Carbon Emission Density of Industrial Land (CEDIL)
The concept of density has its origins in physics and is a measure of mass in a given volume, measured by dividing the mass of an object by its volume.Sociology, economics, geography, ecology, and urban and territorial spatial planning have drawn on the physical concept of density to introduce the concepts of population density, economic density, species density, and urban density, all using area instead of volume in the measurement of these emerging concepts.Drawing on the measurements of emerging concepts such as population density and economic density, this study proposes the concept of carbon intensity, which is measured by dividing the total amount of carbon emissions in a given area by its area.In particular, the carbon intensity of industrial land use is equal to the total industrial carbon emissions of each city divided by the total area of industrial land use, in tons per square kilometer.With   to represent the total industrial land area (square kilometer) of the th city in the study area and   to represent the total industrial carbon emissions of the th city, the carbon intensity (CEDIL  ) of the industrial land is calculated as follows: CEDIL  stands for carbon emissions per unit of industrial land area, which is a relative concept compared to the absolute concept of total carbon emissions and represents the average carbon emissions from industrial land within a city.Carbon intensity and carbon emissions per capita are also common relative concepts, with the former referring to the total carbon emissions of a region divided by its GDP and the latter referring to the total carbon emissions of a region divided by its population.Total carbon emissions, carbon intensity, and per capita carbon emissions are commonly used indicators in urban and territorial spatial planning; however, they have a major drawback: weak correlation with land.Urban and territorial spatial planning manages the land and the elements it carries (population, buildings, businesses, etc.), and density control, especially density zoning, is a key tool for planning.Notably, the carbon intensity is further mapped by integrating both land and carbon emissions into a spatial zoning scheme for carbon intensity.Therefore, CEDIL  is more advantageous than the indicators carbon emission scale, intensity, and per capita emissions and is more applicable to the preparation and management of urban and territorial spatial planning.

Boston Consulting Group (BCG) Matrix
The Boston Matrix originated in the field of corporate development strategy analysis and is a good method for quantitatively measuring the development trends of a company's business or product market.By integrating two indicators, such as the relative share and growth rate of a business or products, it is possible to select the most appropriate company plan from among four types of alternative development strategies (star, gazelle, cow, and dog) [42].With the average growth rate (growth rate,   ) from 2010 to 2020 to represent the trend of change in the time dimension, and the relative share (relative share,   ) in 2020 to quantitatively measure the regional ecological niche of each city from a spatial perspective, the Boston Matrix enables the quantitative portrayal of the spatio-temporal evolutionary pattern of CEDIL in the YREB.CEDIL − and CEDIL − represent the carbon emission density of industrial land in the th city in the end and base periods, respectively, and CEDIL − represents the maximum value in the study area at the end.  and   are calculated using the following equation [43]: The median values of   and   were used as thresholds in this study, because they can eliminate artificial interference in the analyzed results and provide a more balanced classification scheme.According to the classification results based on the spatiotemporal evolution model, the 110 cities in the YREB are classified into four types, with highcarbon and high-growth cities (HCHG cities) representing the   and   of the YREB CEDIL  , which are greater than the median, indicating a very bad state.Low-carbon and low-growth cities (LCLG cities) represent the   and   of the YREB CEDIL  and are less than the median, indicating a desirable state.High-carbon and low-growth cities (HCLG cities) show that the   of the YREB CEDIL  is less than the median, while the   is greater than the median, with a positive trend from bad to desirable conditions.Low-carbon and high-growth cities (LCHG cities) present a   for the YREB CEDIL  that is greater than the median, while the   is less than the median, with a tendency to change from desirable to bad conditions.

Global Moran's Index and Cold-Hot Spot Analysis
Spatial effect detection is a prerequisite for the scientific selection of regression analysis models, and if the dependent variable has strong spatial autocorrelation and heterogeneity, the spatial measurement model is selected to have a higher accuracy than the traditional statistical regression model.In this study, we use Moran's I to quantify the global autocorrelation of the total industrial land area, total industrial carbon emissions, industrial land carbon emission density, and their spatiotemporal evolution models.Zero is a key threshold for measuring spatial autocorrelation properties, where values greater than zero represent positive autocorrelation.The maximum and minimum values of the Moran index are 1 and −1, respectively, with a larger absolute value indicating stronger spatial correlation.To further visualize the local spatial pattern of spatial autocorrelation, this study uses the cold-hot spot analysis tool and introduces the  −    * index to identify the areas where high and low values are clustered.With   as the spatial weight matrix, 1 representing spatial adjacency, and 0 representing spatial disjacency, Moran′s I and   * () are calculated as follows [44,45]:

Geodetector
Geodetector is the most commonly used spatial econometric model.Its main advantage is that it can measure the driving effect of a single factor and analyze the correlations between different factors.We employed Geodetector to analyze the driving mechanism of the carbon emission intensity of industrial land, quantifying the influences of measurement factors, pathways, and interactions [46].Firstly, we discretized the source data of the independent variable (  ).During the research process, the quantile method using Python was employed to divide the 110 cities in the study area into multiple zones and form multiple zoning schemes.Secondly, we overlayed the geographical distribution pattern of the dependent variable with the independent variable zoning scheme and analyzed the similarity between the two as measured using Geodetector and characterized by the q-index.Thirdly, we took the scenario with the largest q-index as the final calculation for each factor and conducted the impact analysis of the different factors.It should be noted that the maximum value of the q-index is 1, with larger values representing a stronger single-factor influence or dual-factor interaction.For example, for the driving factors   and   , q(  ) and q(  ∩   ) represent their single-factor driving force and dualfactor interaction, while Min(q(  ), q(  )), Max(q(  ), q(  )), q(  ) + q(  )) represent the minimum, maximum, and sum of their single-factor driving forces.In the analysis of the interaction, the interactions were classified into five types by analyzing the dual-factor interaction in comparison with the maximum value of the one-factor interaction effect, the minimum value, and the sum of the two [47].When q(  ∩   ) < Min(q(  ) , q(  )) and Min(q(  ), q(  )) < q(  ∩   ) < Max(q(  )), q(  )), it is determined that the interaction between the driving factors   and   is nonlinear weakened and single nonlinear weakened, respectively.It indicates that there is antagonism between the two, and it is necessary to avoid the combination of the two in the future policy design.When q (  ∩   ) > Max(q(  ), q(  )) and q(  ∩   ) > q(  ) + q(  ), it is determined that the interaction is of bifactor enhancement and non-linear enhancement, respectively.It indicates that the two are synergistic, and it is necessary to encourage the combination of the two in the future policy design.When q(  ∩  ) = q(  ) + q(  ), it signifies that two factors are independent of each other and that factors can be freely combined with other factors during the policy design [48,49].The following is the calculation formula for the relevant parameters [50]: where ℎ represents the number of partitions in the discretization result of the independent variable data (the value range is an integer from 2 to 10),  ℎ represents the number of cities in partition h, N = 110,  ℎ 2 represents the variance of the dependent variable in partition h,  2 represents the variance of the dependent variable in the YREB study area, and SSW and SST represent partition variance and total variance, respectively.

Indicator Selection and Data Source
The geographical pattern of carbon emission density of urban industrial land and its change is a complex process, affected by the combined effect of economic scale, industrial structure and industrialization level, government intervention, degree of opening up, scientific and technological levels, and other factors.Based on the research experiences of related scholars, we chose the gross domestic product (GDP) to represent the impact of the scale effect of the urban economy in this paper [51,52], the proportion of the tertiary industry (PTI) and per capita GDP (PCGDP) to represent the industrial structure and the industrialization process [53,54], the fiscal expenditure (FE) to represent the impact of government intervention [55], the international trade goods export (ITGE) and foreign direct investment (FDI) to represent the influence of openness to the outside world and the degree of globalization [56], and the number of patent authorizations (NPA) to represent the impact of innovation [57].It is notable that the distribution and evolution of carbon density come close to the natural environment, especially the geomorphology [58], but they are not included in the analytical framework of this paper due to the lack of the necessary complete data.
Urban industrial carbon emission data came from the China City CO2 Emissions Dataset, released by the China City Greenhouse Gas Working Group (CCG) on its website.CCG is an urban greenhouse gas research and cooperation platform initiated and organized by the Center for Climate Change and Environmental Policy Research of the Academy of Environmental Planning of the Ministry of Ecology and Environment (MEP), established in June 2017.It regularly releases long time-series, full-caliber, and full-coverage data on greenhouse gas emissions from Chinese cities.The industrial carbon emission list covers carbon emissions from the ore, chemical, metal, and electronics industries, as well as carbon emissions related to the production and use of other products, including industrial energy consumption and industrial process emissions.The data on urban industrial land use came from the China Urban Construction Statistical Yearbook, released by the Ministry of Housing and Urban-Rural Development, which provides statistical data on urban construction land use over a long time series and across all land types.Urban industrial land is defined as land used for independently established factories, workshops, handicraft workshops, production sites for construction and installation, and slag (ash) discharge sites.There are three classes of industrial land according to the degree of disturbance to the residential and public environment, the degree of pollution, and the degree of safety hazards.And some cities, to meet the demand for the transformation and upgrading of traditional industries into a high-tech, synergistic production space, combined production space, and headquarters economy, have allocated land for new industries integrating research and development, creativity, design, pilot testing, non-polluting production, and other innovative industrial functions and their ancillary services to industrial land as well.The impact factor data came from the China Urban Statistical Yearbook and were released by the Department of Urban Socio-economic Survey of the National Bureau of Statistics (Appendix A).

Industrial Carbon Emissions
From 2010 to 2020, the city with the highest industrial carbon emissions in the YREB shifted from Shanghai to Suzhou-JS, with a decrease from 190.77 million to 143.67 million tons.Huangshan has been the city with the smallest industrial carbon emissions, with a decrease from 620,000 to 260,000 tons.Different from the decreasing trend of the maximum and minimum values, the average and median industrial carbon emissions in the YREB both increased to a certain extent, from 24.6636 and 17.2250 tons to 26.5923 and 20.3150 million tons, respectively.It should be noted that the coefficient of variation of industrial carbon emissions in the YREB decreased from 1.20 to 0.99, much higher than 0.36, showing high spatial heterogeneity despite the decline [59] (Table 1).The geographic distribution pattern of industrial carbon emissions in the YREB was visualized through quartile spatial clustering analysis using GIS software(ArcGIS 10.4 produced by Esri company).It was found that in 2010, high-carbon cities formed four agglomerations in the Yangtze River Delta (YRD) city cluster, Chengdu-Chongqing city cluster, and southern Hunan, including Shanghai, Suzhou-JS, Chongqing, Ningbo, Wuxi, Wuhan, Nanjing, Xuzhou, Qujing, Bijie, and Liupanshui.However, low-carbon cities were agglomerated in western Yunnan, Anhui, and northern Sichuan, including Lu'an, Bozhou, Ya'an, Ziyang, Lincang, Baoshan, Nanchong, Fuzhou, Suining, Suizhou, Bazhong, Zhangjiajie, and Huangshan.The geographical distribution patterns of high-carbon cities showed that in 2020, the YRD urban cluster rapidly expanded to northern Jiangsu, and the carbon emissions of cities in central Zhejiang decreased rapidly; the Chengdu-Chongqing urban cluster was connected to the northern Hubei, with cities such as Shiyan and Xiangyang having transformed from low-carbon cities to high-carbon cities.Most of the high-carbon cities in Hunan, Yunnan, and Guizhou were transformed into low-carbon cities.In 2020, low-carbon cities were mainly found in the peripheral areas of provinces, especially in their border areas, including southwestern Yunnan, the northeast of Sichuan, the west and south of Hunan, the Yunnan-Sichuan-Guizhou border area, and the Jiangsu-Zhejiang-Anhui border area (Figure 2).

Urban Industrial Land
From 2010 to 2020, Shanghai remained the largest user of industrial land in the YREB, despite a decrease from 739 to 538 km 2 .The city with the smallest industrial land use changed from Tongren to Pu'er, increasing from 0.14 to 0.49.The average and median industrial land in the YREB showed some growth, from 31.61 and 15.26 to 37.06 and 22.50 km 2 , respectively.The coefficient of variation of industrial land in cities in the YREB decreased from 2.41 to 1.74 from 2010 to 2020, showing a spatial heterogeneity trend consistent with industrial carbon emissions, but with a higher degree of spatial differentiation (Table 1).
Most of the cities with a high value of industrial land in 2010 were clustered in the YRD, while Chongqing and the Yunnan-Guizhou border area were also the centers of industrial land clustering in the region, including Shanghai, Chongqing, Nanjing, Suzhou-JS, Wuhan, Ningbo, Chengdu, Kunming, Hefei, Hangzhou, and Wuxi.Cities with low industrial land value were mostly in Yunnan, southwestern Hunan, and northern Sichuan, including Huaihua, Lishui, Yibin, Shaoyang, Ya'an, Yingtan, Pu'er, Yuxi, Zhaotong, Baoshan, Zhangjiajie, Lijiang, Bazhong, Lincang, and Tongren.During the 2010-2020 period, the YREB pushed hard for industrial structure adjustment and upgrading, and Shanghai even proposed a development strategy based on resource reduction, prompting major changes in the distribution pattern of urban industrial land.Most of the cities with a high value of industrial land in 2020 were concentrated in the border areas of Sichuan, Yunnan, and Guizhou, including Shanghai, Chongqing, Wuhan, Chengdu, Ningbo, Nanjing, Suzhou-JS, Hangzhou, Changzhou, Hefei, and Nantong.Most of the low-value cities with industrial land were clustered in the northern Jiangsu-Anhui border area and the Chengdu-Chongqing urban cluster, including Dazhou, Zhaotong, Shaoyang, Yuxi, Lishui, Liupanshui, Ezhou, Ziyang, Zhangjiajie, Huaihua, Lijiang, Baoshan, Bazhong, Lincang, and Pu'er (Figure 2).Of note, the comparison between Figures 2 and 3 shows that the distribution pattern of industrial carbon emissions is inconsistent with the industrial land in the YREB, and there is a large spatial mismatch between the two.The results of the independent analysis separating the two do not well reveal the real condition of carbon emissions from industrial land.

Carbon Emission Density of Industrial Land
From 2010 to 2020, the city with the largest industrial carbon emission density in the YREB changed from Tongren to Pu'er, with the value decreasing from 59.71 to 16.55 million tons/km 2 .Huangshan remained the city with the lowest industrial carbon emissions, dropping from 104,400 to 24,700 tons/km 2 .The average industrial carbon emission density in the YREB decreased from 2.0165 to 1.5706 million tons/km 2 , while the median increased from 918,700 to 970,800 tons/km 2 .It should be noted that the variation coefficient of industrial carbon emission density in the YREB decreased from 2.86 to 1.48, and the spatial heterogeneity level was much higher than that of industrial carbon emission and industrial land use (Table 1).

Spatial Effects Detection
Table 2 shows that Moran's I is greater than zero for all indicators, indicating that the geographic patterns and spatial and temporal evolution patterns of industrial carbon emissions, industrial land use, and industrial carbon intensity in the YREB are in significant positive spatial autocorrelation.The spatio-temporal evolution patterns of industrial carbon emissions from 2010 to 2020 and the geographic patterns of industrial carbon emissions in 2010 and 2020 showed a high autocorrelation index.From the industrial carbon perspective, the hot spots in 2010 and 2020 were concentrated in the Shanghai metropolitan area, with a shrinking spatial scope.In 2010, most of the cold spot cities were in eastern Sichuan and southwestern Yunnan, and in 2020, the cold spot cities expanded significantly, covering most of the YREB.They formed a large agglomeration in the middle reaches of the Yangtze River, centered on Hunan and Jiangxi, and extended into the lower reaches in western Anhui and southern Zhejiang.In addition, they formed a second agglomeration in the upper reaches of the Yangtze River, centered in Yunnan, Sichuan, and Guizhou and extending into southwestern Yunnan.The hot spots of the spatio-temporal evolution pattern from 2010 to 2020 were concentrated in the border area of Jiangsu province, Anhui province, and Zhejiang province (around Nanjing metropolitan area), while the cold spots formed three small clusters in eastern Sichuan, northern Hubei, and southern Hunan.The spatial pattern of "center-periphery" is formed by sub-cold and cold spots and sub-hot and hot spots in combination.It is worth noting that the sub-hot spots expanded in a beaded pattern by leaps to the middle and upper reaches of the Yangtze River region, forming a spatial structure of core-edge circles (Figure 6).According to the industrial land use, hot and sub-hot cities in 2010 were relatively concentrated in the lower reaches of the YREB, especially in the border areas of Jiangsu, Zhejiang, and Anhui.There were hot spot clusters formed in 2020, and most of the subhot spot cities were distributed in Yunnan, with hot spots and sub-hot spots concentrated in the upper reaches of the YREB.Cold spot cities in 2010 were scattered, forming small clusters only in eastern Sichuan.Cold spot cities were mostly concentrated in the YRD urban cluster in 2020, while sub-cold spots were mostly concentrated in the middle reaches of the YREB and extended upward to the Chengdu-Chongqing urban cluster.The spatial pattern of hot and cold spots was characterized by an upstream and downstream gradient rather than a circle-edge structure.Most of the hot spot cities in the evolution model of urban industrial land use from 2010 to 2020 were concentrated in the YRD urban cluster, especially in Jiangsu province, while the cold spots were relatively concentrated in Jiangxi and Yunnan provinces (Figure 6).
In terms of the carbon intensity of industrial land use, the hot spots in 2010 were fewer and located in and around the Shanghai metropolitan area, while the sub-hot spots were spread out in their neighboring areas.The hot spots and sub-hot spots shifted from the lower to the upper YREB in 2020, mostly concentrated in Yunnan and its border area with Guizhou, with a slight increase in the number of hot cities.Cold spots in 2010 covered most areas and were mostly concentrated in the middle reaches of the Yangtze River, centered in Jiangxi and Hunan, extending upstream to northern Guizhou and eastern Sichuan and downstream to western Anhui and the southwest corner of Zhejiang.Cold spot cities in 2020 shrank significantly, mostly concentrated in the YRD urban cluster, with a small number concentrated in eastern Sichuan.Most of the hot spots in the spatio-temporal evolution model from 2010 to 2020 were concentrated in Hunan in the middle reaches of the Yangtze River, with a small part in the upper reaches in Yunnan and the sub-hot cities distributed in their periphery.There were few cold spots, mostly in eastern Sichuan, with a few in the central part of Zhejiang and southeastern Jiangxi.

Single Factor Driving Force
In terms of the geographic pattern, the FDI and FE showed the strongest influence on the carbon emission intensity of industrial land as key factors.The NPA, ITGE, PCGDP, and GDP were not less influential, and they were important factors.The PTI had a weak and not statistically significant influence as an auxiliary factor.In terms of spatio-temporal evolution patterns, the FE and NPA had the strongest influence on the carbon emission density of industrial land as key factors.The ITGE, PCGDP, FDI, and GDP had comparable influences as important factors.The PTI had the lowest influence as an auxiliary factor.Of note is that the FDI was not statistically significant, and both the PCGDP and GDP were not significant (Table 3).To sum up, there were significant differences in the influences of different factors on the carbon emission density of industrial land.From a similarity perspective, the FE was a key factor common to both geographic and spatio-temporal evolution patterns, while the NPA, FE, ITGE, and GDP were important factors, and the PTI was an auxiliary factor common to them.From the perspective of differences, the PCGDP and FDI had more significant influence on the geographical pattern than spatio-temporal evolution, while the PTI was just the opposite.In addition, most factors, including the PCGDP, FE, ITGE, and FDI, showed greater influence on the geographical pattern than the spatio-temporal evolution pattern.The PTI was the opposite of them, while the influence of the GDP and NPA on both remained stable overall.Regarding the interaction type, most of the factor pairs were of nonlinear enhancement, with a few of bifactor enhancement, and there was a significant synergy between different factors.This suggests that the current pattern and change in the carbon emission density of industrial land is the result of a combination of factors.That is, different factors reinforce each other, and each factor has a significant direct impact on the carbon emission density.The additive effects produced by the interaction of different factors cannot be ignored.For geographical patterns, there was a large number of factor pairs of bifactor enhancement, including GDP and ITGE, GDP and FDI, GDP and NPA, PCGDP and ITGE, ITGE and FDI, ITGE and NPA, and NPA and FDI.However, for spatiotemporal evolution patterns, there was a very small number of factor pairs, only NPA and GDP and NPA and ITGE (Tables 4 and 5).In terms of interaction intensity, FE and FDI showed the strongest interaction on geographic pattern, followed by FE and ITGE and FE and GDP showing significant interactions, with a comprehensive force exceeding 0.8 or even 0.9, far higher than other factor pairs, considered as super-factor pairs.FE and NPA, PCGDP and FE, PCGDP and FDI, PCGDP and NPA, PTI and FE, and PTI and FDI also showed high levels of interaction and should therefore be taken into account together in future policy designs.Multifactor pairs showed weaker interactions on spatio-temporal evolutionary patterns than geographic patterns, with FE and FDI, FDI and PCGDP, PCGDP and FE, and PTI and FE leading the way in the interacting force, which can be considered as super-factor pairs.It is important to note that NPA and PCGDP, NPA and FE, NPA and FDI, and GDP and PCGDP also showed non-weak interaction forces, and they are factor combination options that should not be ignored in the policy design.Overall, multifactor interactions are dominated by nonlinear enhancements with significant synergistic effects, and more and more super-factor pairs are emerging.In particular, the interaction of FDI, FE, and NPA with other factors showed a leading advantage, and they are the factors that should be emphasized and prioritized during the portfolio policy design.

Discussion
The analysis shows that there is not a synchronous linear relationship between industrial land and its carbon emissions.Introducing the concept of carbon emission intensity can well reveal the coupling relationship between the two, while better meeting the needs of enabling carbon peaking and carbon neutrality through spatial planning.The comparison of Figures 2-4 shows a large difference between industrial carbon emissions, industrial land use, and industrial land carbon emission density.For example, in 2020, Shanghai was the city with the largest industrial land use in the Yangtze River Economic Belt, but it was not the city with the highest carbon emissions.Instead, Suzhou had the highest carbon emissions, and Pu'er had the largest carbon emission density.The difference is universal, not an exception.For example, in the YRD city cluster, except Shanghai, all other cities have a high level of industrial land use, generally a medium or low level of industrial carbon emissions, and low or medium level of carbon emission densities.These three are not only different in their current patterns, but also increasingly in their evolution.For example, the evolution of industrial carbon emissions, industrial land use, and industrial land carbon intensity in Lincang present LCHG, LCLG, and HCHG patterns, respectively, completely differently.This result suggests that for the goal of carbon peaking and carbon neutrality, it is not enough to limit spatial planning to industrial land and the management of carbon emission scale, and a new concept of carbon emission density needs is needed.
The analysis shows that there is a significant spatial difference in the carbon emission density of industrial land in the YREB, with strong correlations between neighboring cities, and the evolution of density emerges in different patterns, including HCHG, HCLG, LCHG, and LCLG.Additionally, the spatial pattern and evolutionary pattern of carbon emission density of industrial land is the result of the joint action of factors.The FDI and FE play a key role in shaping spatial patterns, while the FE and NPA are key factors influencing evolution patterns.Different factors interact with each other to amplify their effects, so the final determination of the force of the factors should take into account rather than the direct force of a single factor, as well as the influence of the association between other factors, especially the strong interaction of the FE and FDI with other factors.It provides a scientific basis for combining strategies in the policy design.
Urban density gained attention as a major issue for cities after the Industrial Revolution [60].To alleviate the problems of dirtiness and overcrowding in Western industrial cities, scholars from the early 20th century to the 1960s developed two schools of thought in study urban density: decentralization and centralization [61].The former, represented by Howard's Garden City, created the idea of low density, while the latter, represented by Corbusier's Radiant City, created the idea of high density by advocating higher building heights, more green space, and increased road widths.In the face of challenges such as suburban sprawl, inner-city decay, and urban disease, the 1970s and 1990s saw the rise of "new urbanism", "smart growth", and "compact city" concepts, with higher densities dominating urban development [62].City governments hope to address urban disease by promoting the rebirth and revitalization of urban centers through high-density development, mixed-use land use, and transit-oriented development to conserve land and reduce energy consumption [63].Since the beginning of the 21st century, scholars have begun to pay attention to and gradually attach importance to the negative effects of highdensity urban development, such as traffic congestion caused by excessive high-density urban development, the declining living environment quality, insufficient green space, and serious air pollution [64].Based on the value trade-off between "economy, environment and people" and the historical reflections on the development of cities with low and high densities, the inverted U-shaped relationship between urban density and urban development benefits has gradually become the consensus of academia and industry [65].
Carbon emission density is an important dimension of urban density, and based on the results of the above analysis, we believe that it is also consistent with parabolic theory.That is, the carbon emission density increases gradually in the early stage and begins to decrease after reaching a peak (critical value).Both too high and too low carbon emission densities pose problems for the sustainability of urban industrial land use, with a medium density being superior to high and low densities [66,67].Medium-suitable density can balance the weaknesses of both high density and low density to maximize the economic, environmental, and human benefits of urban industrial land management.There are two possible reasons for the low carbon emission density of industrial land.First, the small scale of industrial carbon emissions as a molecule indicates a high level of development of greening and low carbonization of the industrial system.Significantly, for Chinese cities at the current stage, the excessive pursuit of low carbonization or even zero carbonization in industrial development implies high development costs and prices, which goes against the healthy development of industries and cities.Second, the large size of industrial land use as the denominator indicates extensive industrial land development.The non-intensive use of industrial land will lead to a waste of land resources, as well as increased transportation costs, high construction costs of supporting facilities and low interest rates.On the contrary, a high carbon emission density in industrial land may also lead to two unhealthy situations.On the one hand, the large scale of industrial carbon emissions as a numerator indicates that high-emission, high-energy-consumption, and high-pollution industrial sectors account for a high proportion of the industrial system, hindering industrial upgrading and high-quality urban development.On the other hand, the small size of industrial land use as the denominator indicates that the supply of industrial land resources is insufficient.In this context, enterprises only seek to meet the demand by increasing the land plot ratio, and therefore, they put up high-rise buildings in large quantities, boosting the cost of development.
It is important to note that the medium density is not a specific and fixed number but a dynamic range, and the upper and lower limits may not be exactly the same across regions and cities (Figure 7a).The values at the two ends of the medium density can be determined in single-city sample studies using a time-series analysis; in large multi-city sample studies, they need to be obtained using regional comparative analysis.In addition to observing the pattern of change in the carbon intensity, the determination of values should also take into account the local development characteristics and stage needs of the study area.In addition, zoning can be traced back to the theory of locality and is a common tool for spatial planning of cities and territories [68].It is of great value to quantitatively draw the carbon emission density planning zoning of urban industrial land and design and implement differentiated management policies under the goal of "carbon neutral carbon peak" for controlling the carbon emissions of urban industrial land.For specific study areas, their internal members can be divided into several planning zonings using quantile spatial clustering analysis, and the government needs to focus on ultra-high or ultra-low density urban governance.It should be noted that the number of members in different zonings should be in an olive-shaped structure of "small at both ends and large in the middle", and if, on the contrary, there is a dumbbell-shaped structure of being "large at both ends and small in the middle", it indicates unreasonable planning zoning of the carbon emission density of industrial land (Figure 7a).In the process of introducing the concept and algorithm of carbon density into spatial planning, planners need to pay attention to the differences between planning systems of different countries and regions, as well as the regional characteristics of different study areas, especially the industrial supporting features (transportation, logistics, and other productive services) [69].The greatest contribution of this paper is to put forward the concept and measurement method of urban industrial land carbon emission density and discuss its connotation and application value from the theoretical level.Many studies in the past have discussed the relationship between urban density and carbon emissions [70,71].For example, Norman and Ishii argued that the medium density is most conducive to reducing carbon emissions [72,73].However, the variable used in these studies was total carbon emissions rather than density.From the previous analysis, it can be seen that total carbon emissions and density are two completely different concepts, and their geographical patterns and spatio-temporal evolution also have significant differences.In traditional studies, we pay more attention to the scale and intensity of carbon emissions (the relationship with economy) while neglecting the connection between carbon emissions and space (i.e., carbon emission density).Spatial planning requires us to find and solve problems from the perspective of space.The carbon emission scale and intensity indicators do not match the needs of spatial planning, while the carbon emission density solves the problem.To a certain extent, proposing the concept of carbon emission density and its measurement method and establishing the direct relationship between carbon emissions and space is of great academic innovation for spatial planning research.Therefore, the design of policies based on total carbon emissions alone is incomplete, as the insufficient consideration of carbon intensity is likely to cause failure in urban land-use planning under the "dualcarbon" goal.The introduction of spatial econometric modeling to improve the accuracy of the analytical results and conclusions is another marginal contribution of this paper.The industrial carbon emissions, urban industrial land, and carbon emission density of industrial land all have significant spatial effects, while most of the past studies chose linear regression and other traditional statistical models, ignoring the influence of spatial heterogeneity and autocorrelation on variables and their relationships.
As mentioned earlier, medium densities are more suitable for sustainable urban development than high and low densities, and therefore, carbon emission control on the YREB industrial sites should focus on both the low-and high-density ends.The identification of interval thresholds for medium density is a prerequisite for policy design.According to the numerical distribution structure of carbon emission density of industrial land in the YREB in 2020 and the balance of the number of members of different zoning schemes (small at both ends, large in the middle), the suggestion is to set the upper threshold to 1.5 million tons/km 2 and the lower threshold to 500,000 tons/km 2 .Therefore, the proportion of 28.18%, 42.73%, and 29.09% of high-, medium-, and low-density cities in the YREB suggests that the distribution structure conforms to the "olive shape", which is reasonable (Figure 8).For low-density cities including Nantong, Bozhou, Lishui, Bengbu, Jingzhou, Chongqing, Ya'an, Ziyang, Fuyang, Zhuzhou, Mianyang, Hefei, Shaoxing, and Xiangyang, they are faced up with a low-density spread of industrial land (Table 6).They need to be guided by the concept of compact cities and smart growth in the future to moderately control the spread and expansion of industrial land, coupled with industrial upgrading and transformation, to gradually increase the carbon emission density of industrial land.In particular, Chengdu, Zigong, Suining, Nanchong, Suizhou, and Huangshan have excessively low carbon densities, and they urgently need to control the spread of low densities on industrial land.In addition, low-density cities should be active in participating in the regional carbon emission trading for industrial land and enhance their capacity for sustainable development through inter-city ecological compensation transfers.
In addition, we suggest that all cities in the YREB include the carbon emission density of industrial land in their land transfer standards to strengthen their direct control over carbon emissions.The current criteria for the granting of industrial land in the YREB cities focuses mainly on economic efficiency and is not sufficiently responsive to "dual-carbon" objectives.For example, the policy provided in the "Standard Land" Control Indicators for Industrial Land in Mianyang City takes the intensity of fixed assets investment (10 4 yuan/mu), the average output value per mu (10 4 yuan/mu), the average tax per mu (10 4 yuan/mu), and the plot ratio of industrial land as mandatory indicators, and it has no environmental constraints related to carbon emissions.Therefore, it is recommended that all cities incorporate the density and total amount of carbon emissions from industrial land into land transfer standards and management policies as soon as possible, as a binding indicator for the selection of urban industries and structural optimization, so as to promote the low carbonization of the industrial system and economic development.
It is worth to note that close cooperation between different cities in the YREB has been established, and as a world-class urban agglomeration, the YREB is closely linked to other cities around the world [74].And the analysis above shows that the carbon emission density of industrial land in the YREB shows significant spatial heterogeneity and autocorrelation.Therefore, it is necessary to develop policies to promote carbon peaking and carbon neutrality in an orderly manner based on the differences in the carbon emission densities of industrial land in different cities and to emphasize collaboration between different cities in achieving the goals, including carbon emission trading between internal cities and exchange and cooperation between external cities. Regional cooperation to promote the low-carbon transformation of industrial structure is to fully leverage the comparative advantages of each region under the constraint of carbon emission density, promote inter-regional industrial transfer and integrated development, jointly cultivate innovative and green industries, and optimize the combination of production and innovation factors in a larger spatial range.
In the end, we suggest that all cities in the YREB should highlight the synergy among multiple factors in designing policy systems and organically combine the policies to strengthen the indirect management of carbon emissions.For example, the FDI, FE, and NPA, as super-interacting factors, have a leading advantage in interacting with other factors, and significantly amplify their influence on the carbon emission intensity of industrial land.Therefore, for globalization, government intervention, and technological innovation represented by the FDI, FE, and NPA [75,76], it is necessary to design dedicated supporting policies, and more importantly, to highlight their combined relationship with related policies such as investment density, urban planning, and plot ratio control to optimize management performance [77,78].

Conclusions
Industry is one of the largest sources of carbon emissions, and the management of carbon emissions from industrial land is a key part to achieve the goal of "double carbon".This study conducted an empirical research based on the 2010-2020 carbon emission dataset of urban industrial land use in the YREB and proposed the concept of carbon emission density of urban industrial land use and its measurement method from the perspective of urban and territorial spatial planning.The findings are as follows: First, the geographic patterns of industrial carbon emissions, urban industrial land, and carbon emission density of industrial land did not agree with their spatio-temporal evolution.The management of industrial land under the goal of "carbon neutrality and carbon peaking" is incomplete if it is only concerned with the total amount of carbon emissions, and it is necessary to incorporate the carbon emission intensity of industrial land into the control system as soon as possible.Second, the carbon emission density of urban industrial land showed significant spatial heterogeneity, with large differences between different cities.This suggests that the design of carbon emission intensity control systems for industrial land use needs to be city specific, with differentiated management policies for different zonings.Third, the carbon emission intensity of urban industrial land showed significant positive spatial autocorrelation.This demonstrates the need to strengthen intercity cooperation in carbon intensity management, with a focus on pushing spatial synergy in the upper, middle, and lower reaches of the YREB, especially in urban clusters, metropolitan areas, and interprovincial border areas.Fourth, the geographical pattern and spatio-temporal evolution model of urban industrial land carbon emission density had a very complex driving mechanism, with prominent features of diversification and ranking of influencing factors, and significant interaction effects among multiple factors.Therefore, in the scheme design of the control system, the influence of key factors should be emphasized, and attention should be paid to the coordination of multiple factors to realize the "organic combination of policies".
It should be noted that due to the difficulty of data quantification and the limitation of data availability, we neglected the empirical analysis of single-city cases at the micro-scale in this study when analyzing the geographic pattern, spatial, and temporal evolutionary patterns of carbon emission intensity of the YREB's multi-city industrial land at the macro-scale, as well as its spatial effects and driving mechanisms.Additionally, the dynamics of carbon reduction capacity are closely related to national, local, and urban strategies, but they have not been empirically analyzed due to the lack of suitable quantitative indicators.

Data Availability Statement:
The land data is sourced from the Statistical Yearbook of Urban Construction published by the Ministry of Housing and Urban Rural Development of the People's Republic of China, while the carbon emission data is sourced from the China City Greenhouse Gas Working Group's dataset on carbon dioxide emissions from Chinese cities (http://www.cityghg.com/toCauses?id=4, 2 November 2023), the remaining data is sourced from the China Urban Statistical Yearbook published by the National Bureau of Statistics.
Acknowledgments: Thank you to Kaixu Zhao, Weiwei Li and all the reviewers for their contributions and assistance in the process of improving the paper.

Figure 2 .
Figure 2. The geographical pattern of ICE in the YREB.

Figure 3 .
Figure 3.The geographical pattern of UIL in the YREB.

Figure 4 .
Figure 4.The geographical pattern of CEDIL in the YREB.

Figure 5 .
Figure 5.The spatio-temporal evolution model of ICE, UIL, and CEDIL in the YREB.

Figure 6 .
Figure 6.The spatial effect of ICE, UIL, and CEDIL in the YREB.

Figure 7 .
Figure 7.The functional relationship between carbon emission density and combined benefits of urban industrial land.(a) Carbon emission density benefit curve; (b) Carbon emission density zoning map.

Figure 8 .
Figure 8. Threshold selection of CEDIL in the YREB.

Table 1 .
Descriptive statistics on ICE, UIL, and CEDIL in the YREB.

Table 2 .
Spatial autocorrelation analysis on ICE, UIL, and CEDIL in the YREB.

Table 3 .
The impact of single factors on the geographical pattern and spatiotemporal evolution pattern of CEDIL in the YREB.

Table 4 .
The multifactor interaction effect on the geographical pattern of CEDIL in the YREB.

Table 5 .
The multifactor interaction effect on the spatio-temporal evolution pattern of CEDIL in the YREB.

Table 6 .
The planning zoning of CEDIL in the YREB.

Normalized Data on Industrial Land, Carbon Emissions, and Its Density Urban Industrial Land Industrial Carbon Emissions Carbon Emission Density 2010
Interest:The authors declare no conflicts of interest.