How to Realize Synergistic Emission Reduction in Future Urban Agglomerations: Spatial Planning Approaches to Reducing Carbon Emissions from Land Use: A Case Study of the Beijing–Tianjin–Hebei Region

: Land use changes in rapidly urbanizing regions around the world constitute a principal anthropogenic element fueling the surge in carbon emissions. Here, land use patterns within the Beijing–Tianjin–Hebei (BTH) urban agglomeration under low-carbon development (LCD) scenarios were simulated. Additionally, social network analysis was employed to formulate carbon balance planning guidelines for various administrative regions. (1) In the ecological protection scenario, carbon emissions from land use were 643.42 × 10 4 tons lower compared to the natural development scenario. Counties with high ecological support coefficients accounted for 22%, making them better suited for predicting outcomes related to low-carbon-oriented land use. (2) The spatial connections of carbon emissions in BTH were closely related, forming the three main carbon emission spatial linkage areas. (3) A carbon balance zoning plan for the BTH in 2035 under the LCD scenario was formulated. Furthermore, key areas for the implementation of carbon peak and carbon neutrality projects were delineated, and targeted measures for carbon reduction and sink increase were proposed. This study provides a new perspective for implementing territorial spatial planning in Chinese urban agglomerations and can aid the government in formulating a reasonable low-carbon-oriented regional planning policy.


Introduction
Global warming is a growing concern and has become one of the major environmental challenges confronting the world presently.Carbon dioxide (CO 2 ) is a greenhouse gas that is sourced mainly from man-made emissions, and its excessive emissions put significant pressure on the ecological environment [1,2].Cities are densely populated and have become hot spots for CO 2 emissions, with China's urban areas accounting for 75% of the country's overall carbon emissions [3].Urban agglomeration is a key factor for promoting urbanrural integration and coordinated regional development [4].In the process of achieving a low-carbon transformation of urban and rural areas, it is crucial to formulate sink increase and emission reduction strategies at the urban agglomeration level.
Carbon emissions are limited by land use change, which can affect global climate change and the carbon cycle by altering the surface carbon process and the intensity of human energy consumption [5].Carbon emissions resulting from land use/cover changes (LUCCs) surpass 30% of global man-made carbon emissions [6,7].Therefore, scientifically simulating LUCCs and exploring land use patterns under future low-carbon development (LCD) scenarios both aid in finding ways to reduce carbon emissions.
Land use carbon emissions can be calculated in several ways.Many scholars have characterized CO 2 emissions by calculating the carbon emissions from energy consumption [8].This method has some limitations, however.The calculation results of this method can only be expressed numerically, making it difficult to locate specific geographic spaces and to demonstrate the similarities and differences in the spatial distribution of carbon emissions.Some studies have used a combination of land use layout and carbon emission factors to compute carbon emissions [9].This method can effectively express the spatial differences in carbon emissions and provide a planning basis for carbon reduction measures in different regions, filling the gap in calculating carbon emissions based on energy.The future land use layout in different scenarios can be used to guide the use of various methods, including the cellular automata (CA) model [10,11], CLUE-S model [12], FLUS model [13], and PLUS model [14].Among these methods, the PLUS model provides benefits of high simulation accuracy and rapid data processing speed.By introducing land expansion analysis strategies and multitype random patch-seed generation mechanisms, it effectively simulates the intricate evolution process of multiple land types.At present, multiple studies have made notable contributions to carbon reduction in Chinese urban agglomerations from the perspective of emission reduction technologies.Although the development space for technology emission reduction is constrained, the importance of a top-down systematic planning methodology is gradually becoming prominent [15].As a hot spot for greenhouse gas emissions, urban agglomerations exhibit the uneven distribution of CO 2 because of their imbalanced socioeconomic factors [16].Exploring collaborative approaches to emission reduction in urban agglomerations can better demonstrate the importance of the urban network structure in reducing CO 2 emissions in the process of urbanization, thereby better guiding emission reduction planning.The economic development levels and natural environments of different cities and counties vary, and as a result, their carbon emissions also vary.These differences directly affect different regions' ability to achieve carbon peak targets.Therefore, policy makers must develop different carbon reduction strategies based on specific circumstances and the inherent emission differences of each city and county [17].The economic and population mobility relationships among counties are complex.The carbon emissions of each county within the urban cluster are greatly affected by surrounding counties, making it difficult to separate a county from the urban cluster and address its carbon issues separately.Therefore, when studying carbon emission reduction issues in urban agglomerations, considering the connections of carbon emissions between cities is essential.Social network analysis (SNA) serves as an effective methodology for examining the spatial structure of urban agglomerations and the spatial connections among cities.SNA can be employed to explore the spatial correlation characteristics of carbon emissions in urban agglomerations and the specificity of different cities and counties in carbon emission networks.Lately, numerous scholars have utilized SNA to examine spatial coupling, spatial correlation, and dynamic evolution of the spatiotemporal patterns of carbon emissions among regions in terms of the spatial scale [18,19].
Still, many gaps exist in current research on carbon emissions within urban agglomerations.Because of a lack of data, most studies can focus only on the carbon emissions of some large cities with recorded data, and the research locations are individual, not universal.In addition, due to limited interannual records of historical data, the data analysis results are applicable only to the policies and economic conditions at that time, making it difficult to align with recent carbon emission reduction governance policies.Moreover, most carbon emission simulations focus only on simulating future land use patterns or providing macropolicy recommendations to reduce emissions and enhance foreign exchange in the research area as a whole [20].Few studies have simulated the land use patterns of LCD scenarios or formulated detailed and feasible planning guidelines for different administrative regions informed by the location and socioeconomic attributes of cities [21].The classification types and methods of carbon balance zoning, as well as the planning guidelines for different Land 2024, 13, 554 3 of 26 regions to undertake the functions of increasing emissions and reducing emissions, remain unclear, and the theoretical basis needs to be improved.
Combining the PLUS model with the SNA method can compensate for this lack of data and can simulate the spatial carbon emission linkage in the context of formulating LCD policies in the future.Therefore, these methods can be used to more effectively formulate planning strategies and achieve the goal of increasing foreign exchange and reducing emissions.
The Beijing-Tianjin-Hebei (BTH) urban amalgamation stands as one of China's four national-level urban agglomerations, and it has a high proportion of heavy industry and significant deterioration of ecological functions [22], which puts significant pressure on economic development, industrial upgrading, and integrated development.The significant energy consumption and exhaust emissions in the BTH urban agglomeration have caused serious air pollution [23].Thus, coordinating the development of the BTH integration is of overriding significance to plan for carbon balance zoning and formulate a scientific linkage emission reduction scheme.Given this imperative, our study delved into the land use pattern under the future LCD scenario and the spatial association of carbon emissions.We also developed a low-carbon-oriented carbon balance zoning planning and management scheme.

Study Area
The BTH includes 13 cities (113 1).Its total area spans 21.5 million km 2 and is situated in the northeastern coastal plain of China.As the BTH is one of four national urban agglomerations, its coordinated and integrated development can effectively improve the layout and form of urban agglomerations and contribute to the construction of an environmentally friendly society.In 2015, the state released the Outline for Coordinated Development of the Beijing-Tianjin-Hebei Region, which clarifies the significance of ecological environmental protection and the synergy between emission reduction and sink increase.
Land 2024, 13, 554 3 of 28 different administrative regions informed by the location and socioeconomic attributes of cities [21].The classification types and methods of carbon balance zoning, as well as the planning guidelines for different regions to undertake the functions of increasing emissions and reducing emissions, remain unclear, and the theoretical basis needs to be improved.
Combining the PLUS model with the SNA method can compensate for this lack of data and can simulate the spatial carbon emission linkage in the context of formulating LCD policies in the future.Therefore, these methods can be used to more effectively formulate planning strategies and achieve the goal of increasing foreign exchange and reducing emissions.
The Beijing-Tianjin-Hebei (BTH) urban amalgamation stands as one of China's four national-level urban agglomerations, and it has a high proportion of heavy industry and significant deterioration of ecological functions [22], which puts significant pressure on economic development, industrial upgrading, and integrated development.The significant energy consumption and exhaust emissions in the BTH urban agglomeration have caused serious air pollution [23].Thus, coordinating the development of the BTH integration is of overriding significance to plan for carbon balance zoning and formulate a scientific linkage emission reduction scheme.Given this imperative, our study delved into the land use pattern under the future LCD scenario and the spatial association of carbon emissions.We also developed a low-carbon-oriented carbon balance zoning planning and management scheme.

Study Area
The BTH includes 13 cities (113°27-119°50 E, 36°05-42°40 N) (Figure 1).Its total area spans 21.5 million km 2 and is situated in the northeastern coastal plain of China.As the BTH is one of four national urban agglomerations, its coordinated and integrated development can effectively improve the layout and form of urban agglomerations and contribute to the construction of an environmentally friendly society.In 2015, the state released the Outline for Coordinated Development of the Beijing-Tianjin-Hebei Region, which clarifies the significance of ecological environmental protection and the synergy between emission reduction and sink increase.

Data Sources
The data utilized in this study comprise land use data, soil type data, annual mean temperature data, annual mean precipitation data, DEM data, nighttime light data, GDP data, population data, railway data, highway data, national highway data, river water data, national nature reserve data, and carbon emission calculation data (Table 1).

Methodology
The proposed planning simulation framework in this study consists of three main components (Figure 2).

1.
Land use pattern under the future LCD scenario: We used the PLUS model to import a dataset that includes land use data spanning three periods from 2010 to 2020, as well as climate, environment, and future driving factors to investigate the land use pattern under both the natural development (ND) scenario and ecological protection (EP) scenario in the BTH in 2035.We analyzed the carbon emission characteristics under the different scenarios and used the carbon emission and ecological support coefficient (ESC) as the criteria to screen the development scenarios that are more conducive to LCD scenarios.

2.
Spatial connection analysis of carbon emissions under the LCD scenario: We introduced the carbon emission and socioeconomic data into the gravity model to formulate the carbon emission gravity matrix for the BTH region.We used SNA to assess the spatial network characteristics of carbon emissions in the BTH in 2035 under the LCD scenario.

3.
Carbon balance zoning planning under the 2035 LCD scenario in the BTH: Drawing from the land use pattern under the LCD scenario and the spatial network characteristics of future carbon emissions, we analyzed the status and role of each region in increasing sinks and reducing emissions.On this basis, we formulated a future carbon balance zoning plan for the BTH to pinpoint critical areas for implementing the emission peak and carbon neutrality projects.
late the carbon emission gravity matrix for the BTH region.We used SNA to assess the spatial network characteristics of carbon emissions in the BTH in 2035 under the LCD scenario.3. Carbon balance zoning planning under the 2035 LCD scenario in the BTH: Drawing from the land use pattern under the LCD scenario and the spatial network characteristics of future carbon emissions, we analyzed the status and role of each region in increasing sinks and reducing emissions.On this basis, we formulated a future carbon balance zoning plan for the BTH to pinpoint critical areas for implementing the emission peak and carbon neutrality projects.

Multi-Scenario Land Use Pattern Projections
The PLUS model was derived from the FLUS model [26].Drawing from the prediction of land use demand, the Markov module developed two modules: (1) land expansion strategy analysis (LEAS) and ( 2) the CA model of multitype random patch seeds (CARS).CARS can be used for refined land use simulation and prediction and can be coupled with a multi-objective optimization algorithm.The simulation outcomes can provide enhanced support for formulating planning policies to attain sustainable development goals.
To explore the land use patterns under various development scenarios, in the PLUS model, we established two development scenarios (ND and EP) to simulate the land use patterns in the BTH region for the year 2035.The ND scenario continued the development trend from 2010 to 2020 and maintained the transfer probability and neighborhood weight of the original land use types.Under the EP scenario, the primary goal orientation was the protection of the ecological environment.It aimed to curb the unordered expansion of urban areas, decrease the likelihood of converting forest land, shrub, grassland, and water bodies into built-up land by 20%, elevate the likelihood of converting cultivated land into forest land by 60%, and augment the likelihood of converting cultivated land into shrub and grassland by 50% when setting the neighborhood weight [27].These two scenarios both took the national nature reserves as the limited expansion area.

Carbon Emission Calculation Model
Carbon emissions from the BTH region were computed based on land use patterns under various scenarios derived from the PLUS simulation.Built-up land and cultivated land served as the primary carbon sources in the urban agglomeration, whereas forest land, grassland, water bodies, and unused land acted as the primary carbon sinks in the

Multi-Scenario Land Use Pattern Projections
The PLUS model was derived from the FLUS model [26].Drawing from the prediction of land use demand, the Markov module developed two modules: (1) land expansion strategy analysis (LEAS) and ( 2) the CA model of multitype random patch seeds (CARS).CARS can be used for refined land use simulation and prediction and can be coupled with a multi-objective optimization algorithm.The simulation outcomes can provide enhanced support for formulating planning policies to attain sustainable development goals.
To explore the land use patterns under various development scenarios, in the PLUS model, we established two development scenarios (ND and EP) to simulate the land use patterns in the BTH region for the year 2035.The ND scenario continued the development trend from 2010 to 2020 and maintained the transfer probability and neighborhood weight of the original land use types.Under the EP scenario, the primary goal orientation was the protection of the ecological environment.It aimed to curb the unordered expansion of urban areas, decrease the likelihood of converting forest land, shrub, grassland, and water bodies into built-up land by 20%, elevate the likelihood of converting cultivated land into forest land by 60%, and augment the likelihood of converting cultivated land into shrub and grassland by 50% when setting the neighborhood weight [27].These two scenarios both took the national nature reserves as the limited expansion area.

Carbon Emission Calculation Model
Carbon emissions from the BTH region were computed based on land use patterns under various scenarios derived from the PLUS simulation.Built-up land and cultivated land served as the primary carbon sources in the urban agglomeration, whereas forest land, grassland, water bodies, and unused land acted as the primary carbon sinks in the urban agglomeration [28][29][30][31][32][33].Based on the China carbon emission inventory coefficient, and a coefficient from the relevant literature [29], we formulated the carbon emission coefficient of land use in the BTH region, as presented in Table 2.The formula for the land use carbon emission estimation model is as follows [9]: where C e is carbon emissions, e i is carbon emissions from a specific land use type, A i indicates the area of the ith land use type, and a i is the carbon emission coefficient of the ith land use type.

Ecological Support Coefficient
We used the ecological support coefficient (ESC) to assess the carbon sink capacity under different development scenarios [34,35].ESC can gauge the ecological capacity's contribution to carbon balance in each region, and its calculation formula is as follows: where C i and C represent the carbon emissions from each county-level unit and the entire region, respectively, and CA i and CA are the carbon sequestration values of each county-level unit and the entire region, respectively.If ESC > 1, the contribution rate of carbon absorption exceeds the contribution rate of carbon emissions, and this indicates a robust ecological carrying capacity for the county.If ESC < 1, the region's ecological carrying capacity is considered weak.Taking ESC as one of the evaluation bases, we selected the low-carbonoriented development scenario in the multi-scenario land use pattern prediction study described in Sections 4.1.3and 4.1.4,which encompassed a higher proportion of counties with lower carbon emissions and ESC values > 1, as the LCD scenario.Then, we analyzed the spatial connection of carbon emissions under the subsequent low-carbon scenario.

Spatial Correlation Analysis of Carbon Emissions Based on SNA
In this study, we employed the gravity model to measure the economic and social links between regions.GDP, POP size, and carbon emissions were also incorporated into the gravity model to adjust and measure the intensity of the carbon emissions linkage between the counties in the BTH [36].The calculation formula is as follows: where R is the gravitational value; d is the spherical geographical distance between counties; G, P, and C are GDP, POP, and total carbon emissions, respectively; k is the adjustment factor of the carbon emission ratio of a county in the two-way relationship; and i and j represent two different counties.
The GDP and POP data for 2035 were predicted using the gray (1, 1) model (GM).As the growth and change trends of the GDP and POP data are complex and uncertain, GM is suitable for predicting the GDP and POP data for 2035 [37,38].We selected the GDP data and population data from 2014 to 2020 for prediction.We obtained the GDP data from a previous study [39] and the population data from the LandScan population dataset The LandScan population data produced by the U.S. Department of Energy at Oak Ridge National Laboratory (ORNL) represents ambient population count (averaged over 24 h) at approximately 1 km 2 resolution.
The main characteristics of the spatial network were the network spatial pattern and the centrality of each county.The network spatial pattern was used primarily to visualize the overall structure of the network, indicating the ability to achieve coordinated development among regions.Centrality was employed to examine the individual roles of each county in collaborative emission reduction within the network.It reflects the significance of carbon emission exchange between counties for the coordinated development of the entire network.As a result, we clarified the key counties with different carbon emission reduction attributes.The centrality selection indicators are given in Table 3 [17].
X indicates whether there is a direct connection between counties, where 1 indicates a direct connection and 0 indicates no connection; n is the number of cities in the network Directly reflects the status of nodes in the network Indegree Centrality

Closeness of Centrality
In-closeness Centrality ) denotes the shortest spherical distance from county i to county j Reflects the degree to which each city in the entire network is not dominated by the other cities Out-closeness Centrality g jk represents the total number of shortest paths from county j to county k that generate links; g jk (i) is the number of shortest paths from county j to county k that generate links through county i Reflects the extent to which a node can control the relationships between other nodes 3.4.Carbon Balance Planning for the BTH Urban Agglomeration in 2035 under the LCD Scenario

Carbon Balance Zoning Planning
Based on the simulation results of land use under the future LCD scenario, a carbon balance zoning plan for the BTH region in 2035 was developed.ESC was used to characterize the carbon sink potential of a region, while carbon emission linkage regions were used to characterize the potential for regional collaborative emission reduction.We believe that when ESC > 1, the region exhibits carbon-sink-oriented characteristics; when ESC < 1, the region exhibits carbon-emission-oriented characteristics.When the county is in a carbon emission linkage zone, the correlation between networks is high, making it easy to carry out regional collaborative emission reduction work; when the county is not in the carbon emission linkage zone, the correlation between networks is low, making it difficult to carry out regional collaborative emission reduction work.Utilizing the calculated spatial distribution of ESC and the three major carbon emission linkage zones, we divided the carbon balance zoning according to this standard.

Key Carbon Balance and County Planning
Based on carbon balance zoning planning, the centrality of each county was calculated.The study identified the characteristics of each county in the carbon emission network.Building upon the centrality analysis of each county within the carbon emission network outlined in Section 3.3, four key counties were determined based on the scores of urban agglomerations in the above four indicators: core county, promotion county, bridge county, and improvement county.Among them, counties with higher indegree centrality were designated as core counties, counties with higher in-closeness centrality were designated as promotion counties, counties with higher betweenness centrality were designated as bridge counties, and counties with higher outdegree centrality were designated as improvement counties.4. Regarding the accuracy validation of the PLUS model, the kappa coefficient stood at 0.9149, with the overall simulation accuracy reaching 0.9402.This outcome indicates that the land use predictions generated by the model in this study had high accuracy [27,40].Consequently, the model could be extended for simulating land use types and their distributions in 2035 in future studies.The land use distribution results for 2035 under the two scenarios obtained from the simulation are shown in Figure 3 and Table 5.According to the simulation results, in comparison to land use in 2020, forest land, water bodies, and built-up land would continue to increase, while cultivated land, shrub, grassland, and unused land would decrease in 2035 under the ND scenario.Under the EP scenario, in contrast to the ND scenario, built-up land and cultivated land would decrease by 52,284 ha and 175,628 ha, respectively.The area of ecological land (forest land, shrub, and grassland) would increase by 223,247 ha.Therefore, under the constraints of the EP policy, the disorderly expansion of built-up land in the BTH urban agglomeration in 2035 could be effectively restricted, while the area of ecological land could be augmented.

Ecological Support Coefficient
ESC was employed to assess the carbon sink capacity under the different development scenarios (Figure 5).The spatial distribution trends of ESC under the two development scenarios were similar, and the overall characteristics were as follows: High ESC values were observed in the northwest, while lower ESC values were noted in the central, southern, and eastern regions.Under the ND scenario, 39 counties had an ESC > 1, accounting for 19.5%.Under the EP scenario, 44 counties had an ESC > 1, accounting for 22%.
We selected the development scenarios with lower carbon emissions and more counties with higher ESC values for the land use prediction results for the LCD scenarios.The spatial distributions of the carbon emissions under the two scenarios exhibited higher concentrations in the central and eastern regions, while they were lower in the western and northern regions.The main areas of carbon emission accumulation included Beijing, Tianjin, Tangshan, Cangzhou, and the Binhai New Area.These accumulation areas are characterized by developed industries, high-density population distributions, highly concentrated distributions of built-up land, and a relatively small proportion of ecological land (forest land and grassland) retaining the same characteristics of carbon emission agglomeration observed in 2020.

Ecological Support Coefficient
ESC was employed to assess the carbon sink capacity under the different development scenarios (Figure 5).The spatial distribution trends of ESC under the two development scenarios were similar, and the overall characteristics were as follows: High ESC values were observed in the northwest, while lower ESC values were noted in the central, southern, and eastern regions.Under the ND scenario, 39 counties had an ESC > 1, accounting for 19.5%.Under the EP scenario, 44 counties had an ESC > 1, accounting for 22%.
We selected the development scenarios with lower carbon emissions and more counties with higher ESC values for the land use prediction results for the LCD scenarios.

Carbon Emission Network in the BTH Urban Agglomeration
We used the gravity model to draw a carbon emission correlation intensity map of the BTH under the LCD scenario in 2035 (Figure 6).To illustrate the network structure more clearly, we opted to map the carbon emission correlations corresponding to the top 1% of the gravity values [41].Under the LCD scenario, a distinct spatial correlation pattern of carbon emissions emerged in the BTH region.Beijing, Tianjin, and the neighboring counties of Langfang, Baoding, Tangshan, and Cangzhou in the central area of the study region exhibited close relationships.Conversely, the carbon emission network in the southwestern part of the study area was mainly concentrated in Shijiazhuang and displayed less connectivity with other regions.In the southern region of the study area, the spatial network of carbon emissions was composed mainly of Xingtai and Handan.
According to the gravity values (threshold gravity value = 153,000), the counties with close carbon emission spatial correlations were clustered.We obtained the linkage areas  We used the gravity model to draw a carbon emission correlation intensity map of the BTH under the LCD scenario in 2035 (Figure 6).To illustrate the network structure more clearly, we opted to map the carbon emission correlations corresponding to the top 1% of the gravity values [41].Under the LCD scenario, a distinct spatial correlation pattern of carbon emissions emerged in the BTH region.Beijing, Tianjin, and the neighboring counties of Langfang, Baoding, Tangshan, and Cangzhou in the central area of the study region exhibited close relationships.Conversely, the carbon emission network in the southwestern part of the study area was mainly concentrated in Shijiazhuang and displayed less connectivity with other regions.In the southern region of the study area, the spatial network of carbon emissions was composed mainly of Xingtai and Handan.
According to the gravity values (threshold gravity value = 153,000), the counties with close carbon emission spatial correlations were clustered.We obtained the linkage areas of carbon emissions.Under the LCD scenario, in 2035, the BTH region formed three main groups of carbon emission spatial correlation networks: (1) the Beijing-Tianjin-Tangshan-Baoding-Cangzhou carbon emission aggregation network, with Beijing and Tianjin serving as the dual core hubs; (2) the Shijiazhuang carbon emission aggregation network, which was concentrated in Shijiazhuang Province; and (3) the Xingtai-Handan carbon emission aggregation network, which was composed primarily of Xingtai and Handan.These three carbon emission networks were grouped into three carbon emission linkage areas.

Analysis of the Centrality of Counties in the Carbon Emission Networks
Analyzing and visualizing the centrality of the counties in the networks further clarifies the key counties in the carbon emission networks in the BTH urban agglomeration and offers practical guidance for subsequent planning (Figure 7).The number of counties with the top 10% of the centrality in each city is shown in Table 7.In the centrality evaluation, the counties in Beijing, Tianjin, Baoding, Hengshui, and Cangzhou ranked relatively high, including 46 counties with characteristic nodes.The spatial correlation network of the carbon emissions in the BTH region was centered on the counties included in these

Analysis of the Centrality of Counties in the Carbon Emission Networks
Analyzing and visualizing the centrality of the counties in the networks further clarifies the key counties in the carbon emission networks in the BTH urban agglomeration and offers practical guidance for subsequent planning (Figure 7).The number of counties with the top 10% of the centrality in each city is shown in Table 7.In the centrality evaluation, the counties in Beijing, Tianjin, Baoding, Hengshui, and Cangzhou ranked relatively high, including 46 counties with characteristic nodes.The spatial correlation network of the carbon emissions in the BTH region was centered on the counties included in these cities, which had a significant impact on the spatial spillover of the carbon emissions in urban agglomerations.
Land 2024, 13, 554 13 of 28    Based on the multi-scenario carbon emission simulation, we selected the EP scenario with lower carbon emissions and more counties with high ESC values as the LCD scenario for the BTH urban agglomeration in 2035.Utilizing the ESC analysis, carbon emission spatial network characteristics, and network centrality of each county under this scenario, we formulated low-carbon-oriented carbon balance zoning planning and carbon balance key county planning.

Carbon Balance Zoning Planning
The carbon balance zoning planning was based mainly on the ESC results discussed in Section 4.1.4and the carbon emission linkage area summarized in the carbon emission network within the urban agglomeration discussed in Section 4.2.1.According to the planning method of carbon balance zoning plan in Section 3.4.1,we divided the zones into four types of zones, and the division process is shown in Figure 8.We divided the BTH urban agglomeration into 33 carbon sink balance areas, six ecological conservation areas, 88 carbon emission linkage control areas, and 74 low-carbon control areas (Table 8).
Land 2024, 13, 554 14 of 28 spatial network characteristics, and network centrality of each county under this scenario, we formulated low-carbon-oriented carbon balance zoning planning and carbon balance key county planning.

Carbon Balance Zoning Planning
The carbon balance zoning planning was based mainly on the ESC results discussed in Section 4.1.4and the carbon emission linkage area summarized in the carbon emission network within the urban agglomeration discussed in Section 4.2.1.According to the planning method of carbon balance zoning plan in Section 3.4.1,we divided the zones into four types of zones, and the division process is shown in Figure 8.We divided the BTH urban agglomeration into 33 carbon sink balance areas, six ecological conservation areas, 88 carbon emission linkage control areas, and 74 low-carbon control areas (Table 8).

Ecological conservation area ESC > 1, belongs to the carbon emission linkage area
The ecological conservation area primarily resided in the northwestern region of the Beijing-Tianjin-Tangshan-Baoding-Cangzhou carbon emission linkage area.It accounted for 4.699% of the urban agglomeration area.The region has strong carbon sink capacity and is situated within the carbon emission linkage area, so it can give full play to the carbon sink regulation capacity of the ecological region and drive the overall coordinated emission reduction within the carbon emission network.

Carbon sink balance area ESC > 1, does not belong to the carbon emission linkage area
The carbon sink balance area was located mainly in the northwestern part of Hebei Province.It accounted for 43.322% of the urban agglomeration area and was composed mainly of 33 counties.The region has a strong ecological function, high carbon sink benefit, and low carbon emissions, so it is suitable to introduce economically developed county industries appropriately and continue to carry out ecological environmental construction.

Carbon emission linkage control area
ESC ≤ 1, belongs to the carbon emission linkage area The carbon emission linkage control area was composed mainly of counties in three carbon emission linkage areas.It accounted for 22.207% of the urban agglomeration area and was composed mainly of

Ecological conservation area ESC > 1, belongs to the carbon emission linkage area
The ecological conservation area primarily resided in the northwestern region of the Beijing-Tianjin-Tangshan-Baoding-Cangzhou carbon emission linkage area.It accounted for 4.699% of the urban agglomeration area.The region has strong carbon sink capacity and is situated within the carbon emission linkage area, so it can give full play to the carbon sink regulation capacity of the ecological region and drive the overall coordinated emission reduction within the carbon emission network.

Carbon sink balance area ESC > 1, does not belong to the carbon emission linkage area
The carbon sink balance area was located mainly in the northwestern part of Hebei Province.It accounted for 43.322% of the urban agglomeration area and was composed mainly of 33 counties.The region has a strong ecological function, high carbon sink benefit, and low carbon emissions, so it is suitable to introduce economically developed county industries appropriately and continue to carry out ecological environmental construction.

Carbon-Balance Zoning Division Criteria Functional Features
Carbon emission linkage control area ESC ≤ 1, belongs to the carbon emission linkage area The carbon emission linkage control area was composed mainly of counties in three carbon emission linkage areas.It accounted for 22.207% of the urban agglomeration area and was composed mainly of 87 counties.The region exhibits significant carbon emissions and is integrated into the carbon emission network, enabling it to effectively utilize the influence mechanism of the carbon emission network and establish a regional carbon emission reduction linkage system [20].

Low-carbon control area ESC ≤ 1, does not belong to the carbon emission linkage area
The low-carbon control area was located mainly in the central and southern parts of the urban agglomeration.It accounted for 25.369% of the urban agglomeration area and was composed of 74 counties.The region has high carbon emissions, but it has less contact with other counties in the BTH region and is difficult for the other to affect.It is necessary to strengthen communication with other regions and gradually carry out project cooperation and technical exchanges with provinces and cities.

Key Carbon Balance and County Planning
Based on the centrality analysis of each county in the carbon emission network described in Section 4.2.2, we divided four types of counties according to the planning method of Key Carbon Balance and County Planning in Section 3.4.2.The division process is shown in Figure 9. Based on indegree centrality, outdegree centrality, in-closeness centrality, and betweenness centrality, we identified four types of key counties: core counties, promotion counties, bridge counties, and improvement counties (Table 9 and Figure 10).
Land 2024, 13, 554 15 of 28 87 counties.The region exhibits significant carbon emissions and is integrated into the carbon emission network, enabling it to effectively utilize the influence mechanism of the carbon emission network and establish a regional carbon emission reduction linkage system [20].
Low-carbon control area ESC ≤ 1, does not belong to the carbon emission linkage area The low-carbon control area was located mainly in the central and southern parts of the urban agglomeration.It accounted for 25.369% of the urban agglomeration area and was composed of 74 counties.The region has high carbon emissions, but it has less contact with other counties in the BTH region and is difficult for the other counties to affect.It is necessary to strengthen communication with other regions and gradually carry out project cooperation and technical exchanges with provinces and cities.

Key Carbon Balance and County Planning
Based on the centrality analysis of each county in the carbon emission network described in Section 4.2.2, we divided four types of counties according to the planning method of Key Carbon Balance and County Planning in Section 3.4.2.The division process is shown in Figure 9. Based on indegree centrality, outdegree centrality, in-closeness centrality, and betweenness centrality, we identified four types of key counties: core counties, promotion counties, bridge counties, and improvement counties (Table 9 and Figure 10).

Core county Counties with higher indegree centrality
These counties dominate the central position within the carbon emission network, deserving top priority in shaping regional policies for carbon emission reduction.The optimization of socioeconomic development and energy development top-level design in these core counties is essential.The construction of a carbon emission trading market should be strengthened with market-oriented low-carbon mechanisms.

Counties with higher incloseness centrality
Counties that are more connected to other counties can change their car-bon emission behaviors and the carbon emission behaviors of their neighboring counties more quickly, so they should be used as promotion counties for carbon emission reduction and sink increase, further promoting emission peak and carbon neutrality policies according to local conditions.The government can learn from the experience

Core county Counties with higher indegree centrality
These counties dominate the central position within the carbon emission network, deserving top priority in shaping regional policies for carbon emission reduction.The optimization of socioeconomic development and energy development top-level design in these core counties is essential.The construction of a carbon emission trading market should be strengthened with market-oriented low-carbon mechanisms.

Promotion county Counties with higher in-closeness centrality
Counties that are more connected to other counties can change their car-bon emission behaviors and the carbon emission behaviors of their neighboring counties more quickly, so they should be used as promotion counties for carbon emission reduction and sink increase, further promoting emission peak and carbon neutrality policies according to local conditions.The government can learn from the experience of advanced cities to promote the effective implementation of coordinated emission reduction planning.

Bridge county Counties with higher betweenness centrality
Counties that play a communication role in the network can control the coordinated emission reduction process in other cities in a planned manner.These counties provide network communication and can carry out key planning control.By limiting some functions, the government can minimize the exchange of ineffective carbon emissions, impede unnecessary carbon emission connections.

Improvement county
Counties with higher outdegree centrality These counties exhibit a greater number of spatial spillover relationships within the carbon emission network and bear a higher burden of carbon emissions.Their ecological conservation function is strong, and their spillover effect can be used to carry out ecological environmental construction and avoid ecological degradation caused by overurbanization and industrial development.
Land 2024, 13, 554 16 of 28 of advanced cities to promote the effective implementation of coordinated emission reduction planning.

Bridge county Counties with higher betweenness centrality
Counties that play a communication role in the network can control the coordinated emission reduction process in other cities in a planned manner.These counties provide network communication and can carry out key planning control.By limiting some functions, the government can minimize the exchange of ineffective carbon emissions, impede unnecessary carbon emission connections.

Improvement county
Counties with higher outdegree centrality These counties exhibit a greater number of spatial spillover relationships within the carbon emission network and bear a higher burden of carbon emissions.Their ecological conservation function is strong, and their spillover effect can be used to carry out ecological environmental construction and avoid ecological degradation caused by overurbanization and industrial development.

Multi-Scenario Simulation of Land Use
By setting different constraints, we conducted simulations of the land use patterns within the BTH urban agglomeration in 2035 under different scenarios.In comparison to the ND scenario, the EP scenario saw a transformation of built-up land in Zhangjiakou and Baoding into forest land and grassland, leading to an increase in their carbon sinks and a reduction in their carbon emissions.In addition, under the EP scenario, carbon emissions were low, and the counties with ESC > 1 accounted for a higher percentage of these counties.Carbon emissions in Zhangjiakou and Hengshui decreased significantly.The ecological carrying capacity of the western region of Zhangjiakou improved, transforming from a low-carbon-emission contribution and low-carbon-sink contribution area to a lowcarbon-sink contribution and high-carbon-sink contribution area.After comparing the two simulated scenarios, it is evident that the EP scenario imposed stricter limitations on the expansion of built-up land; retained more ecological land, including forest land, shrub, and grassland; produced fewer carbon emissions; eased the problem of the gathering situation of high-carbon-emission areas in spatial patterns, and had stronger carbon sink benefits and ecological benefits.Therefore, we selected the EP scenario as the prediction result for the LCD scenario in the BTH in 2035 for further research and analysis.These research findings demonstrate that through the preservation of ecological functional spaces, controlling unplanned construction land expansion, and execution of policies promoting low-carbon environmental protection development, it is possible to effectively decrease carbon emission intensity and enhance regional carbon sink efficiency.In future planning and construction, these cities will be more conducive to carrying out EP planning actions and maintaining the carbon sequestration benefits of existing forests and grasslands.
We found that carbon emissions were concentrated primarily in the central and eastern counties of the BTH urban agglomeration under the LCD scenario.These findings align with earlier research [42,43].The ecological support coefficient exhibited higher values in the northwestern region and lower values in the central, southern, and eastern regions.The distribution trend of the overall evaluation results was consistent with the results of a previous study [41], but there were some differences in the distribution of the counties with an ESC < 1. Zhang et al. used the Comprehensive Environmental Data Archive's (CEDA) database as the source of carbon emission data to calculate the ESCs of Beijing, Tianjin, and Hebei from 2005 to 2017.In this study, we calculated carbon emissions by setting multiple development scenarios, which was more conducive to evaluating carbon balance under different development scenarios and providing suggestions for the zoning planning of low-carbon transformation and development in BTH urban agglomerations.

Characteristics of Carbon Emission Spatial Networks
In this study, we utilized SNA to analyze the carbon emission network of the BTH urban agglomeration.The use of SNA to analyze the spatial association of carbon emissions can break through traditional geographic restrictions and identify the intricate spatial network connections in each region [19].The spatial connections of carbon emissions revealed a network structure with Beijing, Tianjin, Tangshan, Baoding, and Cangzhou at its core, alongside clusters located in Shijiazhuang, as well as in Xingtai and Handan.The spatial arrangement of this clustering aligned with findings from prior research [41].This network of carbon emission spatial clusters exhibited denser connections compared to other urban areas.The collaborative ability of carbon emission networks could be expressed by the number of network connections.Clustered areas with a large number of connections had stronger collaborative capabilities.Further analysis of typical cities that share the characteristics of this network, as well as exploration of the degree of centralization and its characteristics, could help identify the factors that affect the coordinated emission reduction of urban agglomerations.As a result, it would be possible to better formulate low-carbon plans with different goal orientations based on these characteristics.
The centrality degree indicated the crucial role of the nodes within the network.The top 10% of counties in terms of indegree centrality were located mainly in Beijing and Tianjin, where the economy is more developed.These counties had the following characteristics: they undertook a large number of external carbon emission relationships, their industrial structure was more advanced, and their dependence on basic raw materials and energy imports was strong.Furthermore, these counties reside within the core of the spatial network of carbon emissions, situated centrally within the Beijing-Tianjin-Tangshan-Baoding-Cangzhou carbon emission aggregation network.These findings indicate that augmenting the node count in counties in terms of indegree centrality helped promote carbon emission network collaboration.The focus of coordinating carbon emission networks should be on counties with high indegree centrality.The top 10% of the counties in terms of outdegree centrality were all located in Hebei, and these counties had small carbon emissions and strong ecological conservation functions.This result may have been due to the ecological resource advantages of these counties, but the technology level was low and the energy utilization rate was low.Although the core position of the network in these counties was higher than that of other resource-based cities with better ecological benefits, it was still in the middle to low end of the value range in the carbon emission network.These results prove that remote geographical locations, underdeveloped economies, and population levels acted as impediments to the coordinated development of carbon emission networks.Therefore, special consideration should be given to these counties when devising policies, given their relatively passive position in the carbon emission network.As it would be difficult for these counties to achieve positive cooperative effects, more policy support should be given to them.
The closeness of centrality primarily characterized the spatial transfer distance of carbon emissions from counties and determined whether or not carbon transfer links with other counties could be produced more quickly.The top 10% of counties with incloseness centrality are located in Beijing, Tianjin, and Baoding.These counties shared the following characteristics: they were highly involved in the spatial network of carbon emission integration and established extensive carbon emission links with the surrounding counties, thus becoming central actors in the network.It remains crucial to enhance the connection between hub cities and surrounding cities, enhance their communication, and enhance the synergy of carbon emission networks, thus enabling them to quickly take advantage of regional linkage mechanisms.The top 10% of counties with out-closeness centrality were concentrated in Hengshui, Baoding, and Cangzhou in Hebei Province.These counties occupied a subordinate position in the network and were easily influenced by other regions to make corresponding changes.
Betweenness centrality represented the communication function of one county with other counties.The top 10% of counties are located primarily in Baoding, Beijing, Cangzhou, Hengshui, and Tianjin.The distribution of these counties was relatively scattered and these counties were located in the outskirts areas of key cities.Counties with higher betweenness centrality were more likely to establish carbon emission networks, which meant that carbon emissions between provinces could be exchanged mainly through these cities.This result proved that enhancing the centrality of betweenness could strengthen the transfer of carbon emissions between provinces and promote regional collaboration in carbon emission networks.The counties with high betweenness centrality formed the central relationships connecting surrounding counties, served as good communication bridges, and held great significance in improving the carbon emission network structure.
In summary, fostering collaboration within the carbon emission network in the BTH urban agglomeration requires boosting the centrality of county nodes and internal connectivity.In counties with higher centrality of entry, key policies should be implemented to increase foreign exchange and reduce emissions.Cities with higher inward proximity should be used as policy pilots to strengthen the connections between counties with higher centrality of intermediaries to drive carbon exchange between provinces.Counties with higher centrality of entry should be given more policy assistance.To optimize coordinated regional emission reduction planning, it is essential to comprehensively manage the spatial structure characteristics of urban agglomerations, regulate the association mode, internal groups, and regional status of carbon emissions, and align with the principles of social and economic development [18].

Recommendations for Key Areas for the Implementation of Emission Peak and Carbon Neutrality Projects
In crafting and executing carbon emission reduction policies, our attention was directed towards the spatial interconnection impact of carbon emissions within the urban agglomeration.We devised a low-carbon-centric carbon balance zoning strategy, leveraging insights from the ESC, spatial network properties of carbon emissions, and network centrality assessment.On this basis, we could further adopt differentiated development strategies for different counties according to the specifics of carbon balance zoning across various regions of the BTH and the functions of key counties in carbon balance.Taking four characteristic counties and two groups of cross-regional cooperation counties as examples, we propose the suggestions given in Table 10 and shown in Figure 11.
The Fangshan District, located in Beijing, is in the ecological conservation area and undertakes the functions of promoting counties and bridge counties.The eastern part of the region is projected to experience a rise in built-up land, while the western part is anticipated to witness an expansion of forest land and grassland.In terms of carbon reduction, we should concentrate on the construction of three high-tech groups in the east and advocate for the development direction of low carbon.To foster compact urban development, efforts will be intensified to promote the construction of bus lanes and dedicated roads [44][45][46].In terms of increasing foreign exchange, we should construct a recreational and cultural green belt in the plain area [47][48][49].

A Low-Carbon-Oriented Urban Agglomeration Planning Framework
Numerous studies have consistently demonstrated the pivotal role of regional spatial planning in mitigating regional carbon emissions and attaining climate objectives [83][84][85].First, spatial planning is comprehensive and can regulate the spatial division of urban agglomerations on a macroscopic scale [86].Second, spatial planning helps determine the characteristics of cities and formulate differentiated development policies based on the social and economic disparities within the urban agglomeration [87].Additionally, spatial planning assists in reorganizing spatial relationships within the region, ensuring the balanced and coordinated development of spaces, including residential, commercial, transportation, service, and ecological areas [84].Given the significance of spatial planning in reducing regional carbon emissions, this study proposes a low-carbon-oriented carbon The Haidian District, located in Beijing, is in the carbon emission linkage control area and serves the functions of core counties and promotion counties.In the forthcoming years, the northern region is expected to witness an augmentation in both built-up and cultivated land, while forest land is projected to expand in the northwest.For carbon reduction endeavors, our focal point should center on Zhongguancun Science City situated in the southern region.It will spearhead the synchronized advancement of low-carbon technology across the carbon emission linkage area [50][51][52] and promote active and public transportation by improving street connectivity [53][54][55].In terms of increasing sinks, in the future, we should establish the ecological barrier of Xishan Mountain in the northwest as the core position [56][57][58].Pingquan City, located in Chengde, is in the carbon sink balance area and serves as an improvement county.In the forthcoming years, we anticipate a rise in forest land, grassland, and cultivated land across the entire region.It is crucial to take advantage of these additional regional ecological benefits.In terms of carbon reduction, we should strengthen Pingquan's support of Beijing's noncapital functions and actively accept carbon emission spillover within the Beijing-Tianjin region [59][60][61].In terms of increasing sinks, we should cooperate with the BTH and Liaoning-Inner Mongolia regions to carry out ecological protection [62][63][64][65].
Land 2024, 13, 554 20 of 26 Hejian City, located in Cangzhou, is in the low-carbon control area and serves as an improvement country and bridge country.In the future, there will be a trend of growth in both built-up and cultivated land across the entire region.Therefore, it is essential to prioritize efforts towards promoting carbon reduction initiatives.In terms of carbon reduction, we should decrease reliance on fossil fuels while concurrently elevating the utilization of nuclear and renewable energy sources [66,67].To augment sinks, it is imperative to establish an ecological security pattern grounded in agriculture, forestry, and water networks.Improve agricultural technology [68], advocate for improving the dietary structure of residents [69][70][71], and reduce the indiscriminate use of fertilizers, organic fertilizers, and chemical fertilizers [72].
The Taocheng District and Jizhou District, located in Hengshui, are in the carbon emission linkage control area and the low-carbon control area, respectively, and serve as bridge counties.In the future, there will be a surge in built-up land on the periphery of the region and an expansion of water bodies in the central area and cultivated land.To advance carbon reduction, fostering interregional collaboration and establishing a national demonstration zone for the conversion of scientific and technological advancements in the southern part of Beijing is crucial [73,74].We should conform to the expansion trend of the central waterbody in the aspect of increasing the sinks and should carry out ecological restoration in the Fuyang River water ecological comprehensive renovation area around the regional central Hengshui Lake National Nature Reserve [75,76].We should also promote scientific planning and the construction of tourist areas [77,78].
Yu County and Zhuolu County are located in Zhangjiakou, Laishui County is located in Baoding, and the Mentougou District and Fangshan District are located in Beijing.These counties are in the carbon sink balance area and ecological conservation area.In the future, they will follow a trend of increasing forest land, grassland, and cultivated land.In terms of increasing sinks, we should establish a way of putting regional coordinated ecological restoration planning into action; take advantage of the characteristic county role of bridge, promotion, and improvement counties [79,80]; carry out industry structure adjustment; and coordinate the implementation of emission peak and carbon neutrality policy projects.For the important nature reserves, we should carry out regional coordinated habitat protection optimization and ecological corridor construction.We should also increase the number of nature reserves [81,82].

A Low-Carbon-Oriented Urban Agglomeration Planning Framework
Numerous studies have consistently demonstrated the pivotal role of regional spatial planning in mitigating regional carbon emissions and attaining climate objectives [83][84][85].First, spatial planning is comprehensive and can regulate the spatial division of urban agglomerations on a macroscopic scale [86].Second, spatial planning helps determine the characteristics of cities and formulate differentiated development policies based on the social and economic disparities within the urban agglomeration [87].Additionally, spatial planning assists in reorganizing spatial relationships within the region, ensuring the balanced and coordinated development of spaces, including residential, commercial, transportation, service, and ecological areas [84].Given the significance of spatial planning in reducing regional carbon emissions, this study proposes a low-carbon-oriented carbon balance spatial planning method for future urban agglomerations (Figure 12).The overall planning method can be split into four steps: (1) By simulating multiple scenarios of land use, we calculated the carbon emissions of urban agglomerations under different scenarios for the future, selected development scenarios align more favorably with low-carbon principles, and optimized the overall spatial scale of regional carbon emissions from the standpoint of urban agglomerations.This approach has been validated in most urban studies [88][89][90], demonstrating that this top-down carbon reduction strategy can be used to adjust future urban land demand, formulate urban expansion plans with minimal ecological losses, and promote compact and green urban development.
(2) Taking into account the carbon emission characteristics and spatial correlation of urban agglomerations, we utilized carbon balance zoning planning to delineate four distinct carbon balance-oriented zones at the internal spatial correlation level of urban agglomerations.This approach clarifies the spatial development pattern of carbon emissions within urban agglomerations and addresses the question of how to perceive spatial carbon effects.Regional planning with different orientations helps divide the "zoning entities" under the concept of functional control of the main body.The "zoning" control emphasizes respecting and balancing the different roles played by different zoning entities in regulating regional carbon balance, to better protect ecological space and coordinate urban spatial development [43,91].(3) Through the analysis of individual network characteristics in social network analysis, we developed the key carbon balance and county planning based on the centrality indicators of each county and analyzed their position within the carbon emission network.Clear guidance and development goals for carbon reduction actions in different counties were established.In the field of urban planning, the literature has shown that governments and policymakers worldwide intend to utilize various data and measurement methods to ascertain urban characteristics; summarize regional political, economic, and social characteristics; and plan cities based on these factors to better address the needs of people regarding housing, economic opportunities, and social system changes [87].(4) Finally, we developed key areas for the implementation of emission peak and carbon neutrality projects by analyzing the land features within the county.By scaling down from the level of urban agglomerations, we analyzed the composition of land parcel elements within the county.By optimizing the structure of industrial land, energy land, ecological land, agricultural land, and transportation land, we further optimized the structure of the carbon-source and carbon-sink land and implemented urban spatial planning accordingly.We guided the practice of county-level urban spatial planning in terms of low-carbon layout, low-carbon industries, low-carbon transportation, low-carbon human settlements, and ecological green spaces to better regulate carbon spatial layout and formulate a carbon intensity design.Urban and transportation planners have conducted various studies on achieving carbon neutrality goals at this level, and combined with the low-carbon planning tasks of county-level urban spaces, they have screened and extracted planning and control elements from both qualitative and quantitative perspectives, and implemented them in practice [45,[92][93][94].

Limitations
While this study examined the interplay between future land use changes and social networks within spatial planning, it is important to acknowledge certain limitations.Specifically, during the process of simulating land use carbon emissions, variations in carbon emission intensities across different land use types and years might have existed, particularly concerning built-up land.In the process of comparing future development scenarios, In summary, we started from four levels: the overall space of the urban agglomeration, the spatial correlation areas within the urban agglomeration, each county, and the internal land features within the county.We optimized the carbon emission scale through scenario simulation, coordinated the carbon emission spatial development pattern through carbon balance zoning planning, clarified the carbon emission action goals through key carbon balance and county planning, and finally, established key areas for the implementation of emission peak and carbon neutrality projects to regulate the carbon spatial layout and formulate carbon intensity designs in response to the carbon effects of land types, space, urban areas, and spatial elements.By starting from a multi-scale spatial pattern, developing multilevel action plans, constructing a planning policy framework, and integrating carbon reduction and sink enhancement goals into multilevel and multi-class spatial planning, we aim to achieve the leading, comprehensive, and fundamental nature of regional carbon reduction policies, thereby more widely guiding their application in major urban agglomerations around the world.

Limitations
While this study examined the interplay between future land use changes and social networks within spatial planning, it is important to acknowledge certain limitations.Specifically, during the process of simulating land use carbon emissions, variations in carbon emission intensities across different land use types and years might have existed, particularly concerning built-up land.In the process of comparing future development scenarios, however, this research method can still effectively screen more low-carbon simulation results, and these results can be used as an optimization guide to build a carbon-neutral territorial spatial pattern system.Additionally, although this study developed a multi-scale and multilevel research framework through calculations and references to articles from various planning fields worldwide, it has not yet been validated through virtual simulation models to evaluate the feasibility of planning policies.In future research, it would be valuable to conduct additional evaluation studies to assess the effectiveness of carbon emission reduction planning initiatives.

Conclusions and Contribution
In this study, we focused on the urban agglomeration comprising BTH county-level units as our research subject.We established a theoretical framework for carbon balance planning from the standpoint of low-carbon development.Future land use patterns were simulated under different development scenarios, and we selected low-carbon-oriented development models.On this basis, we introduced the socioeconomic factors of the urban agglomeration and constructed a spatial association network of the carbon emissions in the BTH urban agglomeration under LCD scenarios.Based on the spatial network characteristics and centrality characteristics, we carried out carbon balance zoning planning, clarified the functional characteristics of the different zoning areas and key counties, delineated the key areas for carbon neutrality projects, and proposed suggestions.The primary findings of the research are outlined below: 1.
We simulated land use patterns under the ND and EP scenarios in the BTH in 2035.Carbon emissions from land use under the EP scenario decreased by 643.42 × 10 4 t.The spatial distribution was characterized by lower carbon emission and higher ESC values in the northwest and higher carbon emissions and lower ESC values in the central, southern, and eastern regions.We selected the EP scenario as the prediction result of the LCD scenario for the BTH in 2035 for further analysis; 2.
Carbon emissions within the BTH were closely related, forming three major carbon emission linkage areas: Beijing-Tianjin-Tangshan-Baoding-Cangzhou, Shijiazhuang, and Xingtai-Handan.The counties of Beijing, Tianjin, Baoding, Hengshui, and Cangzhou and other cities in the central region had different levels of individual centrality.

3.
Based on the characteristics of the ESC and carbon emission network, we formulated a carbon balance zoning plan for the BTH urban agglomeration in 2035 under the LCD scenario.The urban agglomeration underwent division into four areas: eco-

Figure 1 .
Figure 1.Geographical location of the study area.Figure 1. Geographical location of the study area.

Figure 1 .
Figure 1.Geographical location of the study area.Figure 1. Geographical location of the study area.

Figure 2 .
Figure 2. Workflow of this study.

Figure 2 .
Figure 2. Workflow of this study.
Analysis of the BTH Urban Agglomeration in 2035 4.1.1.Simulation Accuracy Verification Ensuring the precision of land use model simulation, according to BTH land use in 2010 and 2015, we conducted simulations and compared the expected results of land use in 2020 with the real land use distribution in the same year.The results are displayed in Table

Figure 3 .
Figure 3. Land use patterns under the (a) ND scenario and (b) EP scenario in 2035.4.1.3.Spatial Distribution of Carbon Emissions We calculated future carbon emissions in 2035 based on land use patterns under different scenarios (Figure 4 and Table 6).Carbon emissions in 2035 would be 41,408.18× 10 4 tons under the ND scenario and 40,764.76× 10 4 tons under the EP scenario.In comparison to the ND scenario, the EP scenario would reduce the emission of greenhouse gases by 643.42 × 10 4 tons.The net carbon sink was larger, leading to lower overall carbon emissions

Figure 3 .
Figure 3. Land use patterns under the (a) ND scenario and (b) EP scenario in 2035.

Figure 4 .
Figure 4. Spatial distributions of carbon emissions under the (a) ND scenario and (b) EP scenario in 2035.

Figure 4 .
Figure 4. Spatial distributions of carbon emissions under the (a) ND scenario and (b) EP scenario in 2035.

28 Figure 5 .
Figure 5. ESC under the (a) ND scenario and (b) EP scenario in 2035.

Figure 5 .
Figure 5. ESC under the (a) ND scenario and (b) EP scenario in 2035.

4. 2 .
Spatial Correlation Analysis of Carbon Emissions in the BTH Urban Agglomeration under the LCD Scenario 4.2.1.Carbon Emission Network in the BTH Urban Agglomeration

Land 2024, 13 , 554 12 of 28 Figure 6 .
Figure 6.Carbon emission networks under the LCD scenario in the BTH urban agglomeration.

Figure 6 .
Figure 6.Carbon emission networks under the LCD scenario in the BTH urban agglomeration.

Figure 8 .
Figure 8. Flowchart of carbon balance zoning planning.

Figure 8 .
Figure 8. Flowchart of carbon balance zoning planning.

Figure 9 .
Figure 9. Flowchart of key carbon balance and county planning.

Table 9 .
Planning of key counties for carbon balance in 2035 under the LCD scenario in the BTH urban agglomeration.

Figure 9 .
Figure 9. Flowchart of key carbon balance and county planning.

Figure 10 .
Figure 10.Carbon balance zoning and key county planning for the BTH urban agglomeration in 2035 under the LCD scenario.Figure 10.Carbon balance zoning and key county planning for the BTH urban agglomeration in 2035 under the LCD scenario.

Figure 10 .
Figure 10.Carbon balance zoning and key county planning for the BTH urban agglomeration in 2035 under the LCD scenario.Figure 10.Carbon balance zoning and key county planning for the BTH urban agglomeration in 2035 under the LCD scenario.

Figure 11 .
Figure 11.Key regional proposals for emission peak and carbon neutrality engineering for the BTH urban agglomeration in 2035 under the LCD scenario.

Figure 11 .
Figure 11.Key regional proposals for emission peak and carbon neutrality engineering for the BTH urban agglomeration in 2035 under the LCD scenario.

Table 1 .
Details of the datasets.

Table 2 .
Carbon emission coefficients of different land use types.

Table 3 .
Spatial network indicators of carbon emissions.

Table 4 .
Comparison between the simulated land use types and the actual land use types in 2020 at pixel level.

Table 5 .
Changes in the regional area under different scenarios in 2020-2025 (ha).

Table 5 .
Changes in the regional area under different scenarios in 2020-2025 (ha).Spatial Distribution of Carbon EmissionsWe calculated future carbon emissions in 2035 based on land use patterns under different scenarios (Figure4and Table6).Carbon emissions in 2035 would be 41,408.18× 10 4 tons under the ND scenario and 40,764.76× 10 4 tons under the EP scenario.In comparison to the ND scenario, the EP scenario would reduce the emission of greenhouse gases by 643.42 × 10 4 tons.The net carbon sink was larger, leading to lower overall carbon emissions and reduced net carbon emissions.

Table 6 .
Carbon emissions under the ND scenario and EP scenario in 2035 (10 4 t).

Table 6 .
Carbon emissions under the ND scenario and EP scenario in 2035 (10 4 t).

Table 7 .
Number of counties in the top 10% of the centrality indicators for each city.
4.3.Carbon Balance Planning for the BTH Urban Agglomeration in 2035 under the LCD Scenario

Table 7 .
Number of counties in the top 10% of the centrality indicators for each city.Carbon Balance Planning for the BTH Urban Agglomeration in 2035 under the LCD Scenario

Table 8 .
Carbon balance zoning planning for the BTH urban agglomeration in 2035 under the LCD scenario.

Table 8 .
Carbon balance zoning planning for the BTH urban agglomeration in 2035 under the LCD scenario.

Table 9 .
Planning of key counties for carbon balance in 2035 under the LCD scenario in the BTH urban agglomeration.

Table 10 .
Key regional proposals for emission peak and carbon neutrality projects in the BTH urban agglomeration in 2035 under LCD scenarios.