Analysis of the Spatiotemporal Evolution of Carbon Budget and Carbon Compensation Zoning in the Core Area of the Yangtze River Delta Urban Agglomeration

: As a world-class urban agglomeration, the Yangtze River Delta urban agglomeration is significant for China’s carbon neutrality strategy when calculating the carbon budget and dividing carbon compensation zones. This paper focused on 129 county-level cities in the core area of the Yangtze River Delta urban agglomeration, calculating the carbon budget for 2000–2020 using night-time light data and analyzing the evolution of spatiotemporal patterns. On this basis, a carbon compensation zoning model was constructed using the K-means algorithm; 129 cities were divided into different carbon compensation zones after combining this model with the main functional zones. The results showed that (1) the carbon emissions increased, with overall stabilization after 2012. The distribution of high carbon emission areas evolved from a “core-periphery” pattern into a “Z-shaped” pattern. Ecological carbon sinks showed a continuous decline, with the southern sinks performing better than the northern ones. (2) The carbon budget of the study area showed a deficit, gradually widening and exhibiting an unbalanced spatial distribution characterized by a “high in the south and low in the north” pattern. (3) Eleven types of carbon compensation zones were designated after overlaying the main functional zones. Low-carbon development suggestions were proposed for each zone type.


Introduction
Carbon emissions from human activities are contributing to global warming, leading to heightened climate extremes, rising sea levels, and food and water shortages, which present substantial threats to human life.Sea levels will rise by 0.6-3 m in 2030 under a high carbon emissions situation [1], and at least 50 major cities will be at risk of inundation globally, mainly in Asia [2].The probability of extreme events increases by 340% every 20 years (under 2 • C warming) [3].Countries accounting for more than 70% of the world's population and global GDP will be exposed to more than 500% of the risk of flooding (under conditions of a 4 • C rise in temperatures) [4], and the frequency of category 4 and 5 hurricanes in the Atlantic will double [5].These changes not only seriously threaten the safety of human life and property but also threaten the global ecological balance.Warming also poses a huge potential threat to China, where the frequency of heavy precipitation and droughts continues to increase [6], as do flood losses [7].China is the largest developing country and the largest emitter of carbon dioxide.In response to the climate crisis, China has committed to reaching the peak of carbon emissions around 2030 and striving to achieve "carbon neutrality" by 2060.Achieving carbon peaking and carbon neutrality is an important strategic decision for China.
Carbon budget is a measure of the balance between carbon emissions and carbon sinks, often expressed as the difference between the two.The traditional methods for Land 2024, 13, 747 3 of 23 erations of North America, the Pacific Coast urban agglomerations of Japan, the British urban agglomerations, and Northwest Europe urban agglomeration.The core area of the Yangtze River Delta urban agglomeration (CAYRDUA) is an area with a high concentration of the economy, population, industry, and carbon emissions within the YRDUA, along with strong economic radiation capability.It is an important engine of the YRDUA.The region accounted for 31.78% of the land area of the YRDUA, 52.23% of the total population, and 53.49% of the GDP of the YRDUA in 2020.It boasts a high urbanization rate of 79.17% and accounts for 63.96% of the total electricity consumption in the YRDUA.County-level cities in the CAYRDUA are actively developing, with significant differences in the economic base, functional positioning, and industrial structure.For example, the Chongming District in Shanghai is currently developing as a carbon-neutral demonstration area, while Kunshan City in Suzhou has a concentrated manufacturing industry, a highly developed economy, a solid industrial base, and a substantial energy demand, with its GDP surpassing that of three provincial-level administrative divisions in 2020.Therefore, the study of carbon budget in the CAYRDUA must be conducted at the county level to better reflect the actual situation and development pattern of the study area.Such studies are important in identifying low-carbon development pathways in both global and Chinese contexts.
Calculating carbon emissions often relies on energy consumption data, which are generally unavailable at the municipal and county levels in China and are not conducive to long-time-series carbon emission studies.Carbon emission calculation methods based on night-time light data can remedy this problem.Night-time light data have been shown to have a linear relationship with carbon emissions [42,43] and can explain 96% of the variation in carbon dioxide emissions among Chinese provinces [44].In this paper, based on two types of night-time light data (DMSP/OLS and NPP/VIIRS), land use data, energy consumption data, and other socio-economic data, we calculate the carbon budget of county-level cities in the CAYDRUA from 2000 to 2020 and analyze their spatiotemporal distribution and evolutionary patterns.The normalized revealed comparative advantage (NRCA) index is used to quantify the three attributes of socio-economy, ecology, and the scale of carbon emissions.Additionally, the K-means algorithm is used to classify the carbon compensation zones in the CAYRDUA, and reasonable low-carbon development suggestions are proposed for each zone in combination with the main functional zones.
The structure of this paper is as follows.Section 2 describes the study area, the data used, and the research framework.Section 3 presents the results of carbon emissions, carbon sinks, and carbon budget calculations for CAYDRUA and analyzes the temporal evolution and spatial evolution of the carbon budget.In addition, the result of carbon compensation zoning is shown, and low-carbon development recommendations are given for each zoning.Section 4 discusses the differences between the results of similar research and similar studies with similar urbanization patterns, as well as the trend of carbon budget development in CAYDRUA.Finally, Section 5 exhibits the conclusions of the research.

Study Area
The Yangtze River Delta region features a flat terrain, with the northern part comprising plains below ten meters above sea level and the southern part consisting of hilly and mountainous areas.The region has a subtropical monsoon climate with forest coverage of 32.93%, dense river networks, great development potential, and abundant natural resources [45].The CAYRDUA (118 • 22 ′ -123 • 25 ′ E, 28 • 01 ′ -33 • 25 ′ N) is an urban agglomeration consisting of 8 prefecture-level cities in Jiangsu Province, 7 prefecture-level cities in Zhejiang Province, and Shanghai Municipality, including 129 county-level cities (Figure 1).Forestland in the CAYRDUA is mainly distributed in the south, while arable land is mainly distributed in Jiangsu Province in the north, and the proportion of construction land in Shanghai is higher than in other areas [46].The region accounts for less than 2% of the national territory area, with a population of 122.95 million, contributing to a GDP of CNY 13.08 trillion, with a large population density of 1080 people per km² (Table S1).
Zhejiang Province, and Shanghai Municipality, including 129 county-level cities (Figure 1).Forestland in the CAYRDUA is mainly distributed in the south, while arable land is mainly distributed in Jiangsu Province in the north, and the proportion of construction land in Shanghai is higher than in other areas [46].The region accounts for less than 2% of the national territory area, with a population of 122.95 million, contributing to a GDP of CNY 13.08 trillion, with a large population density of 1080 people per km² (Table S1).With the rapid urban expansion and economic development, on the one hand, the region's demand for fossil energy has led to a rapid increase in carbon emissions, with per capita carbon emissions 1.6 times higher than the national average [40].On the other hand, the loss of carbon sinks due to the conversion of cropland and forestland into construction land has posed a huge challenge to the ecological carrying capacity of the region.In the context of the "dual carbon" target strategy, an analysis of this region's carbon budget can provide a reference for formulating carbon reduction plans and green development policies in the developed areas of China.

Data
In this study, the remote-sensing data used include DMSP/OLS night-time light data from 2000 to 2013 (https://ngdc.noaa.gov/,accessed on 17 July 2022) and monthly NPP/VIIRS night-time light data from 2012 to 2020 (https://eogdata.mines.edu/,accessed on 17 July 2022).Administrative division data were downloaded from the National Catalogue Service For Geographic Information (https://www.webmap.cn/,accessed on 9 September 2022).Land use data from 2000 to 2020 were obtained from the CLCD land use dataset (https://zenodo.org/records/5816591,accessed on 27 September 2023).Energy consumption data at the provincial level from 2000 to 2020 were taken from the China Energy Statistical Yearbook.Other socio-economic data were obtained from statistical yearbooks, government-issued data, and statistical bulletins of the cities in the study area.The main functional zoning came from the planning documents issued by the governments of Jiangsu Province, Zhejiang Province, and Shanghai Municipality.With the rapid urban expansion and economic development, on the one hand, the region's demand for fossil energy has led to a rapid increase in carbon emissions, with per capita carbon emissions 1.6 times higher than the national average [40].On the other hand, the loss of carbon sinks due to the conversion of cropland and forestland into construction land has posed a huge challenge to the ecological carrying capacity of the region.In the context of the "dual carbon" target strategy, an analysis of this region's carbon budget can provide a reference for formulating carbon reduction plans and green development policies in the developed areas of China.

Data
In this study, the remote-sensing data used include DMSP/OLS night-time light data from 2000 to 2013 (https://ngdc.noaa.gov/,accessed on 17 July 2022) and monthly NPP/VIIRS night-time light data from 2012 to 2020 (https://eogdata.mines.edu/,accessed on 17 July 2022).Administrative division data were downloaded from the National Catalogue Service For Geographic Information (https://www.webmap.cn/,accessed on 9 September 2022).Land use data from 2000 to 2020 were obtained from the CLCD land use dataset (https://zenodo.org/records/5816591,accessed on 27 September 2023).Energy consumption data at the provincial level from 2000 to 2020 were taken from the China Energy Statistical Yearbook.Other socio-economic data were obtained from statistical yearbooks, government-issued data, and statistical bulletins of the cities in the study area.The main functional zoning came from the planning documents issued by the governments of Jiangsu Province, Zhejiang Province, and Shanghai Municipality.

Method
The research process of this paper is shown in Figure 2.

Calculation of Carbon Emissions
In this paper, night-time light data were used in combination with the IPCC method to calculate county-level carbon emissions [47].Eight fossil fuels, including raw coal, coke, crude oil, gasoline, kerosene, diesel oil, fuel oil, and natural gas, were selected to calculate the provincial carbon emissions, and then, the carbon emissions at the county level were Land 2024, 13, 747 5 of 23 estimated based on the correlation between night-time light data and carbon emissions.The calculation formula is as follows: where CO 2 represents the carbon emissions from eight fossil fuels; E i represents the consumption of the type of ith fossil energy converted to standard coal consumption (SCE); K i represents the carbon emission factor of the type of ith fossil energy.The coefficients of conversion coefficients of the consumption of different energy sources to SCE, the percentage of usage for the eight fossil fuels, as well as the coefficients of the carbon emission factors, are shown in Table 1.

Method
The research process of this paper is shown in Figure 2.

Calculation of Carbon Emissions
In this paper, night-time light data were used in combination with the IPCC method to calculate county-level carbon emissions [47].Eight fossil fuels, including raw coal, coke, crude oil, gasoline, kerosene, diesel oil, fuel oil, and natural gas, were selected to calculate the provincial carbon emissions, and then, the carbon emissions at the county level were estimated based on the correlation between night-time light data and carbon emissions.The calculation formula is as follows: where CO2 represents the carbon emissions from eight fossil fuels; Ei represents the consumption of the type of ith fossil energy converted to standard coal consumption (SCE); Ki represents the carbon emission factor of the type of ith fossil energy.The coefficients of conversion coefficients of the consumption of different energy sources to SCE, the percentage of usage for the eight fossil fuels, as well as the coefficients of the carbon emission factors, are shown in Table 1.The processing of night-time light data consisted of three parts: the correction of DMSP/OLS data, the correction of NPP/VIIRS data, and the integration of the two nighttime light data.
(1) The correction of DMSP/OLS night-time light data referred to the results of Lv et al., including mutual correction, saturation correction, and continuous correction [49].(2) The correction of NPP/VIIRS data used the method of processing monthly data into annual data.The reference maximum for outlier removal was determined by selecting the brightest pixel among those of Beijing, Shanghai, and Guangzhou.Any outliers lighter than the brightest pixel were adjusted to the highest pixel value, while pixel values below 0.3 × 10 −9 W•cm −2 •SR −1 were identified as background noise [50].(3) A linear regression without intercept was chosen to integrate the two night-time light data [51], and a correla-tion between DMSP/OLS and NPP/VIIRS county-level total DN (TDN) in 2012 and 2013 was constructed.The integrated results are shown in Figure 3, with an R 2 value of 0.9471.Based on this, a night-time light dataset was constructed for 2000-2020.
(1) The correction of DMSP/OLS night-time light data referred to the results of Lv et al., including mutual correction, saturation correction, and continuous correction [49].(2) The correction of NPP/VIIRS data used the method of processing monthly data into annual data.The reference maximum for outlier removal was determined by selecting the brightest pixel among those of Beijing, Shanghai, and Guangzhou.Any outliers lighter than the brightest pixel were adjusted to the highest pixel value, while pixel values below 0.3 × 10 −9 W•cm −2 •SR −1 were identified as background noise [50].(3) A linear regression without intercept was chosen to integrate the two night-time light data [51], and a correlation between DMSP/OLS and NPP/VIIRS county-level total DN (TDN) in 2012 and 2013 was constructed.The integrated results are shown in Figure 3, with an R 2 value of 0.9471.Based on this, a night-time light dataset was constructed for 2000-2020.A linear regression without intercept was established between TDN and carbon emissions at the provincial level (e.g., Table 2) [52].A linear regression without intercept was established between TDN and carbon emissions at the provincial level (e.g., Table 2) [52].The carbon emissions at the county level were calculated based on the linear relationship between provincial TDN and carbon emissions [53], and an adjustment factor [52] was added to the formula as follows: where CO 2ce represents the county-level carbon emissions calculated using linear regression without the intercept model; A represents the regression coefficient; TDN c represents TDN at the county level; CO 2p represents the carbon emissions at the provincial level calculated using the IPCC method; CO 2pe represents the provincial-level carbon emissions calculated using linear regression without the intercept model; k represents the adjustment factor; CO 2c represents the results of carbon emissions at the county level.

Calculation of Carbon Sink
The calculation of carbon sequestration by vegetation in the aboveground portion was performed by multiplying the area of different vegetation types in the CAYRDUA by the amount of carbon sequestered per unit area.The amount of carbon sequestered per unit area by different vegetation types was derived from previous studies [54] (Table 3).The land use types and areas from 2000 to 2020 were obtained from the CLCD dataset [55].
where S i above represents carbon sequestration by the aboveground vegetation systems in the ith county-level city; j represents the jth vegetation type; A i j represents the area of the jth vegetation type in the ith county-level city; α represents the carbon sequestration coefficient per unit area of the vegetation system.Belowground carbon sequestration only considered plant roots and litter, and according to Li et al., the ratio of belowground-aboveground biomass (root-top ratio) was used to calculate belowground carbon sequestration per year [24].We referred to the results of root-top ratio of different vegetation types estimated by Huang et al. [56], as shown in Table 4.The vegetation in some areas south of the Huaihe River is dominated by broadleaved forests and coniferous forests, with root-top ratios ranging from 0.2 to 0.4 [56].In this paper, the root-top ratio of the forest was taken as 0.2.
S i = S i below + S i above (7) where S i below represents carbon sequestration by belowground vegetation systems in the ith county-level city; β represents the mean root-top ratio of the jth vegetation type; S i represents carbon sequestration of the ith county-level city.The relationship between the amount of carbon sequestered and the mass of carbon dioxide absorbed was established through the equation of the process of photosynthesis, from which the carbon dioxide absorbed by the vegetation in the CAYRDUA was calculated annually.

Standard Deviation Ellipse Analysis
The standard deviation ellipse (SDE) was employed to analyze the spatial distribution of attribute values.In this study, the method of SDE analysis was utilized to examine the spatiotemporal patterns of carbon budget and variations in carbon budget centers [57].
where M( x, ŷ) represent the coordinates of the center of SDE; x i y i represent the latitude and longitude coordinates of the center of the spatial observation unit i, respectively; w i represents the weight of the spatial observation unit i; θ represents the azimuth of the SDE; x i y i represent the latitude and longitude coordinate deviations between the center of spatial observation cell i and the mean center, respectively; δ x δ y represents the standard deviation of the extended x-axis and y-axis; n represents the total number of space observation units in the study area.

Spatial Autocorrelation Analysis
In this paper, the global and the local Moran's I indices were used to analyze the general trend of aggregation and the local spatial clustering pattern of the carbon budget from 2000 to 2020.The global Moran's I index reflects the spatial autocorrelation between different county-level cities, and the local Moran's I index reflects the local similarities and differences between county-level cities [48].The calculation formula is as follows: where I represents the global Moran's I index value; I i represents the local Moran's I index value; n represents the number of administrative units; x i and x j represent the carbon budget of the ith and jth county-level cities; x represents the average of the carbon budget; W ij represent the spatial weighting matrix of the ith and jth county-level cities.

Calculation of Carbon Budget
Carbon budget is a measure of the balance between carbon emissions and carbon sinks, often expressed as the difference between the two.The formula for calculating the carbon budget is as follows: where CB represents the carbon budget; CS represents the carbon sink; CE represents carbon emissions.

Construction of Carbon Compensation Zoning Model
In this paper, the three attributes of socio-economy, ecology, and the scale of carbon emission were considered to constitute the carbon compensation zoning index system.The economic contribution coefficient of carbon emission (ECC) and the ecological support coefficient of the carbon sink (ESC) were selected as indicators to characterize the socioeconomic and ecological attributes, respectively [58].The expression of each index is shown below.(1) ECC (16) where G i , G represent the GDP of the ith county-level city and the whole region; C i , C represent the carbon emissions of the ith county-level city and the whole region.The coefficient is used to measure the contribution of regional carbon emissions and the economy.When ECC < 1, it means that the county-level city contributes more to carbon emissions than to the economy.
(2) ESC (17) where CA i , CA represent the carbon sink of the ith county-level city and the whole region.This coefficient measures the region's contribution to the carbon sink and carbon emissions.When ESC > 1, it means that the county-level city contributes more to the carbon sink than to carbon emissions.
In this paper, the normalized revealed comparative advantage (NRCA) index was used to discern the comparative advantage of the different attributes of each county-level city in the CAYRDUA [58], with the following formula: (18) where X j i represents the indicator value of the jth attribute for the ith county-level city; X j represents the value of the total indicator for the jth attribute for all county-level cities; X i represents the total indicator values for all attributes for the ith county-level city.X represents the total indicator value for all county-level cities and attributes.If NRCA j i > 0, it means that the jth attribute of the ith county-level city has a comparative advantage in the study area.
This study constructed carbon compensation zoning by clustering the NRCA index for each attribute using the method of K-means.K-means is a division-based clustering algorithm that divides categories by calculating the Euclidean distance between the sample points and the cluster class center, recalculating the sample center, and iterating until the abort condition is met.

Analysis of Spatiotemporal Distribution of Carbon Emissions
From a temporal perspective (Figure 4), the change in total carbon emissions in the study area from 2000 to 2020 is observed three stages: (1) the rapid growth stage (2000)(2001)(2002)(2003)(2004), which was manifested in the continuous rapid growth of carbon emissions from 43,946.47 × 10 4 tons in 2000 to 65,611.25 × 10 4 tons in 2004, with an average annual growth rate of 10.54%; (2) the stable growth stage (2005)(2006)(2007)(2008)(2009)(2010)(2011), where the growth rate of carbon emissions in the CAYRDUA slowed down but still maintained a positive total growth rate, with an average annual growth rate of 6.8%; (3) the stability stage (2012-2020), where the total amount of carbon emissions in the CAYRDUA was stable in general, increasing from 114,234.87 × 10 4 tons in 2012 to 120,728.18 × 10 4 tons in 2020, with an average annual growth rate of 0.69%.
As shown in Figure 5, in 2000, the high carbon emission regions in the CAYRDUA showed a "core-periphery" distribution, with Shanghai as the center, which expanded in the central and northern parts of the region toward Nantong, Wuxi, and Nanjing, and in the southeastern part of the region toward Ningbo, with the area of high carbon emission gradually expanding.By 2020, the areas with higher carbon emissions exhibited a "Z-shaped" distribution in the central part of the CAYRDUA.As shown in Figure 5, in 2000, the high carbon emission regions in the CAYRDUA showed a "core-periphery" distribution, with Shanghai as the center, which expanded in the central and northern parts of the region toward Nantong, Wuxi, and Nanjing, and in the southeastern part of the region toward Ningbo, with the area of high carbon emission gradually expanding.By 2020, the areas with higher carbon emissions exhibited a "Zshaped" distribution in the central part of the CAYRDUA.As shown in Figure 5, in 2000, the high carbon emission regions in the CAYRDUA showed a "core-periphery" distribution, with Shanghai as the center, which expanded in the central and northern parts of the region toward Nantong, Wuxi, and Nanjing, and in the southeastern part of the region toward Ningbo, with the area of high carbon emission gradually expanding.By 2020, the areas with higher carbon emissions exhibited a "Zshaped" distribution in the central part of the CAYRDUA.

Analysis of Spatiotemporal Distribution of Carbon Sink
From a temporal perspective (Figure 6), the carbon sinks in the CAYRDUA showed a continuous downward trend, from 24,438.34 × 10 4 tons in 2000 to 21,783.84 × 10 4 tons in 2020, with an average annual growth rate of −0.57%.Spatially, the carbon sinks in the CAYRDUA exhibited a distribution pattern of "high in the south and low in the north", with the distribution pattern of carbon sinks in the study area remaining largely unchanged.Low-value carbon sink areas were primarily concentrated in the central and northern plain areas, while high-value areas were concentrated in the southern hilly regions.
a continuous downward trend, from 24,438.34 × 10 4 tons in 2000 to 21,783.84 × 10 4 tons in 2020, with an average annual growth rate of −0.57%.Spatially, the carbon sinks in the CAYRDUA exhibited a distribution pattern of "high in the south and low in the north", with the distribution pattern of carbon sinks in the study area remaining largely unchanged.Low-value carbon sink areas were primarily concentrated in the central and northern plain areas, while high-value areas were concentrated in the southern hilly regions.

Analysis of the Spatiotemporal Distribution and Evolution of the Carbon Budget
The carbon emissions in the CAYRDUA far exceed the carbon sinks, indicating a significant imbalance in the carbon budget.The carbon budget declined rapidly from −19,508.13 × 10 4 tons in 2000 to −98,944.34 × 10 4 tons in 2020, a reduction of 407.20%.By 2020, the carbon sink in the study area was only 18.04% of the carbon emission.The carbon budget change in the CAYRDUA is shown in two stages (Figure 7): (1) the rapid decline period (2000-2011), when the carbon budget deficit continued to expand, decreasing by 373.13% in 2011 compared to 2000.This period witnessed rapid industrial development in the study area, leading to increased energy demand.Simultaneously, rapid urbanization resulted in extensive conversion of cropland and forestland into construction land; (2) the relatively stable period (2012-2020), when the carbon budget in 2020 decreased by 7.62% compared with 2012, indicating a generally stable trend of change.This indicates that the continuous deterioration of the carbon budget in the study area has improved significantly.However, it is still negative, which is not conducive to achieving "carbon neutrality" in the study area.

Analysis of the Spatiotemporal Distribution and Evolution of the Carbon Budget
The carbon emissions in the CAYRDUA far exceed the carbon sinks, indicating a significant imbalance in the carbon budget.The carbon budget declined rapidly from −19,508.13 × 10 4 tons in 2000 to −98,944.34 × 10 4 tons in 2020, a reduction of 407.20%.By 2020, the carbon sink in the study area was only 18.04% of the carbon emission.The carbon budget change in the CAYRDUA is shown in two stages (Figure 7): (1) the rapid decline period (2000-2011), when the carbon budget deficit continued to expand, decreasing by 373.13% in 2011 compared to 2000.This period witnessed rapid industrial development in the study area, leading to increased energy demand.Simultaneously, rapid urbanization resulted in extensive conversion of cropland and forestland into construction land; (2) the relatively stable period (2012-2020), when the carbon budget in 2020 decreased by 7.62% compared with 2012, indicating a generally stable trend of change.This indicates that the continuous deterioration of the carbon budget in the study area has improved significantly.However, it is still negative, which is not conducive to achieving "carbon neutrality" in the study area.
As shown in Figure 8, there exists a significant spatial imbalance in carbon budget distribution, exhibiting a pattern of "higher in the south and lower in the north "in the CAYRDUA.Regions with severe carbon deficits are concentrated in the central and northern parts of the CAYRDUA, primarily in Jiangsu and Shanghai, gradually expanding northwestward and southeastward.Conversely, areas with carbon surpluses are predominantly distributed in southern cities, such as Hangzhou, Shaoxing, Taizhou, and Huzhou.The carbon budget in Zhejiang Province is significantly better than that in Jiangsu and Shanghai.In 2000, the carbon budget in Zhejiang Province exceeded that of Jiangsu Province by 13,887.73×10 4 tons, and this gap has widened over time.By 2020, this figure had increased to 32,077.07 × 10 4 tons.This spatially significant imbalance is attributed to a combination of various factors.On the one hand, Shanghai is a traditional financial center with intense land development and high population density.Meanwhile, Jiangsu is experiencing rapid development in the heavy industry [40], leading to a continuous increase in the proportion of carbon emissions from heavy industry, with Jiangsu province accounting for more than 50% of the total carbon emissions in the CAYRDUA from 2014 to 2019.On the other hand, Zhejiang Province in the southern part of the study area is characterized by mountainous terrain, favorable hydrothermal conditions, lush vegetation, and strong carbon sequestration capacity.In contrast, Jiangsu Province and Shanghai are characterized by flat terrain, with a predominant development of agricultural and construction land, resulting in relatively poorer carbon sequestration capacity.As shown in Figure 8, there exists a significant spatial imbalance in carbon budget distribution, exhibiting a pattern of "higher in the south and lower in the north "in the CAYRDUA.Regions with severe carbon deficits are concentrated in the central and northern parts of the CAYRDUA, primarily in Jiangsu and Shanghai, gradually expanding northwestward and southeastward.Conversely, areas with carbon surpluses are predominantly distributed in southern cities, such as Hangzhou, Shaoxing, Taizhou, and Huzhou.The carbon budget in Zhejiang Province is significantly better than that in Jiangsu and Shanghai.In 2000, the carbon budget in Zhejiang Province exceeded that of Jiangsu Province by 13,887.73×10 4 tons, and this gap has widened over time.By 2020, this figure had increased to 32,077.07 × 10 4 tons.This spatially significant imbalance is attributed to a combination of various factors.On the one hand, Shanghai is a traditional financial center with intense land development and high population density.Meanwhile, Jiangsu is experiencing rapid development in the heavy industry [40], leading to a continuous increase in the proportion of carbon emissions from heavy industry, with Jiangsu province accounting for more than 50% of the total carbon emissions in the CAYRDUA from 2014 to 2019.On the other hand, Zhejiang Province in the southern part of the study area is characterized by mountainous terrain, favorable hydrothermal conditions, lush vegetation, and strong carbon sequestration capacity.In contrast, Jiangsu Province and Shanghai are characterized by flat terrain, with a predominant development of agricultural and construction land, resulting in relatively poorer carbon sequestration capacity.
In terms of carbon budget changes at the prefectural level, Hangzhou exhibited the most favorable carbon budget situation from 2000 to 2003, with Taizhou following closely behind.However, after 2003, the carbon budget deficit in Hangzhou increased, and the carbon budget rapidly decreased from 2003 to 2008, with an annual average growth rate of −253.49%.Suzhou had the worst carbon budget situation among these cities.As a traditional industrial powerhouse, Suzhou's carbon emissions far exceed those of other prefectural cities, comprising 15.11% of the total carbon emissions in the CAYRDUA in 2020, while its carbon sink accounted for less than 3%.From 2000 to 2020, the shape of the SDE in the CAYRDUA exhibited a northwestsoutheast pattern (Figure 9), being predominantly concentrated in the southern part of the CAYRDUA.The center shifted from Zhuji in Shaoxing to Jiande in Hangzhou, indicating significantly less carbon budget loss in the southwestern part compared to the southeastern part of the study area.Shaoxing and the southeastern cities have a weakened influence on the carbon budget patterns.The lengths of the axes, the reduction in the lengths of the short and long axes, as well as the narrowing distribution range of the SDE further indicate the deterioration of the carbon budget in the southern region of the study area.The oblateness increased from 0.47 in 2000 to 0.61, indicating a strengthening directionality of the carbon budget, with a more concentrated distribution in the northwest-southeast direction.In terms of carbon budget changes at the prefectural level, Hangzhou exhibited the most favorable carbon budget situation from 2000 to 2003, with Taizhou following closely behind.However, after 2003, the carbon budget deficit in Hangzhou increased, and the carbon budget rapidly decreased from 2003 to 2008, with an annual average growth rate of −253.49%.Suzhou had the worst carbon budget situation among these cities.As a traditional industrial powerhouse, Suzhou's carbon emissions far exceed those of other prefectural cities, comprising 15.11% of the total carbon emissions in the CAYRDUA in 2020, while its carbon sink accounted for less than 3%.
From 2000 to 2020, the shape of the SDE in the CAYRDUA exhibited a northwest-southeast pattern (Figure 9), being predominantly concentrated in the southern part of the CAYRDUA.The center shifted from Zhuji in Shaoxing to Jiande in Hangzhou, indicating significantly less carbon budget loss in the southwestern part compared to the southeastern part of the study area.Shaoxing and the southeastern cities have a weakened influence on the carbon budget patterns.The lengths of the axes, the reduction in the lengths of the short and long axes, as well as the narrowing distribution range of the SDE further indicate the deterioration of the carbon budget in the southern region of the study area.The oblateness increased from 0.47 in 2000 to 0.61, indicating a strengthening directionality of the carbon budget, with a more concentrated distribution in the northwest-southeast direction.In the spatial autocorrelation analysis, the carbon budget data in CAYRDUA passed the 99% confidence test and the 95% confidence test, as shown in Table 5.The global Moran's I index of 0.32 in 2000 decreased to 0.09 in 2020, indicating that the spatial autocorrelation of CAYRDUA's carbon budget decreased.In 2000, there was a clear pattern of high-high clustering and low-low clustering in the carbon budget of the CAYRDUA (Figure 10).Low-low clustering was mainly observed in Shanghai, as well as in Kunshan and Taicang, two small parts of Suzhou.Highhigh clustering of the carbon budget was concentrated in the southern part of the CAYR-DUA, including some county-level cities in Hangzhou, Taizhou, Shaoxing, Huzhou, Ningbo, and other areas.Over time, the high-high clustering area in the southern region of the study area decreased, splitting into two parts concentrated in the southwest and southeast of the CAYRDUA, with a reduction in the high-high clustering areas of Hangzhou and Shaoxing.This may be attributed to Hangzhou's strategic goal of becoming a "Powerful Industrial City" in 2002, and the growth rate of industrial land in Hangzhou and Shaoxing exceeded the average level of Zhejiang Province from 2005 to 2020.Under the coordinated development of the CAYRDUA, the surrounding cities of Shanghai ac- In the spatial autocorrelation analysis, the carbon budget data in CAYRDUA passed the 99% confidence test and the 95% confidence test, as shown in Table 5.The global Moran's I index of 0.32 in 2000 decreased to 0.09 in 2020, indicating that the spatial autocorrelation of CAYRDUA's carbon budget decreased.In 2000, there was a clear pattern of high-high clustering and low-low clustering in the carbon budget of the CAYRDUA (Figure 10).Low-low clustering was mainly observed in Shanghai, as well as in Kunshan and Taicang, two small parts of Suzhou.High-high clustering of the carbon budget was concentrated in the southern part of the CAYRDUA, including some county-level cities in Hangzhou, Taizhou, Shaoxing, Huzhou, Ningbo, and other areas.Over time, the high-high clustering area in the southern region of the study area decreased, splitting into two parts concentrated in the southwest and southeast of the CAYRDUA, with a reduction in the high-high clustering areas of Hangzhou and

Carbon Compensation Zoning Results and Overlay Analysis
In this study, the NRCA index of three attributes was clustered by the K-means clustering algorithm, dividing the study area into 43 payment zones, 48 balanced zones, and 38 compensated zones.Overlaying with main functional zones, the zones are ultimately categorized into 11 types (Figure 11 and Table 6).

Carbon Compensation Zoning Results and Overlay Analysis
In this study, the NRCA index of three attributes was clustered by the K-means clustering algorithm, dividing the study area into 43 payment zones, 48 balanced zones, and 38 compensated zones.Overlaying with main functional zones, the zones are ultimately categorized into 11 types (Figure 11 and Table 6).

Payment Zone
The payment zone is distributed in the central part of the CAYRDUA, with a parallel distribution in the north and south.Concentrated in areas such as Shanghai, Suzhou, Nanjing, and Ningbo, the payment zone accounts for 34.02% of the total land area, with a highly developed economy (53.69% of GDP), high carbon emissions (63.70% of carbon emissions), and poor ecological carrying capacity (28.71% of carbon sinks).There is a serious mismatch between the economic contribution capacity of carbon emissions and the ecological carrying capacity (the ECC is 0.84, and the ESC is 0.45).This region is an important economic growth pole in the CAYRDUA.

Payment Zone
The payment zone is distributed in the central part of the CAYRDUA, with a parallel distribution in the north and south.Concentrated in areas such as Shanghai, Suzhou, Nanjing, and Ningbo, the payment zone accounts for 34.02% of the total land area, with a highly developed economy (53.69% of GDP), high carbon emissions (63.70% of carbon emissions), and poor ecological carrying capacity (28.71% of carbon sinks).There is a serious mismatch between the economic contribution capacity of carbon emissions and the ecological carrying capacity (the ECC is 0.84, and the ESC is 0.45).This region is an important economic growth pole in the CAYRDUA.
(1) Payment zone-optimized development zone.This region is mainly composed of 31 county-level cities, primarily in Suzhou, Wuxi, Hangzhou, and Ningbo, spatially surrounded by the northern shore of the Taihu Lake and the southern shore of the Yangtze River, accounting for 23.42% of the total area, with the carbon emission contribution slightly higher than the economic contribution (the ECC is 0.85).It is predominantly composed of economically developed and industrially concentrated county-level cities (37.07% of GDP share), facing prominent resource and environmental issues (carbon emissions account for 43.42%, carbon sinks account for 15.19%, the ESC is only 0.35).County-level cities with manufacturing industries as the pillar industry should focus on improving per-unit output and enhancing economic contribution capability in subsequent development.For regions with industries such as metal smelting and rolling processing with high carbon emissions, such as the Zhangjiagang City of Suzhou, the Wujin District of Changzhou, the Jingkou District of Zhenjiang, and the Baoshan District of Shanghai, efforts should be made to promote the consumption of clean energy and introduce efficient and green equipment in production.Furthermore, it is essential to adhere to the "red lines" of arable land and ecology, strengthen the protection of the surrounding environment of the Taihu Lake Basin, construct Taihu Lake ecological wetlands, and enhance carbon sink capacity.
(2) Payment zone-key development zone.This region is situated in the Shanghai Municipality and Jiangsu Province, spatially distributed along the north bank of the Yangtze River and the eastern coastal areas, comprising 11 county-level cities in areas such as Shanghai and Nanjing.Characterized by a solid economic foundation and substantial development potential, this area occupies 8.51% of the total land area, contributing 19.01% of carbon emissions and 10.07% of carbon sinks.Its ecological carrying capacity surpasses that of the Payment zone-optimized development zone (with an ECC of 0.86 and an ESC of 0.53).In subsequent development phases, this area, as the primary recipient of industrial spillover effects, should prioritize the development of low-emission, high-end manufacturing industries.Simultaneously, attention should be directed toward urban green space construction and suburban vegetation coverage.Coastal county-level cities ought to exercise rational control over the development and utilization of tidal flats, minimize human activity disturbance, prioritize the protection of tidal flat ecosystems, and ensure that coastal tidal wetlands fulfill their due carbon sink function.
(3) Payment zone-key ecological functional zone.This region comprises only one county-level city, the Chongming District of Shanghai.It occupies 2.09% of the land area of the CAYRDUA, exhibiting strong ecological functions (with carbon sinks accounting for 3.45%, an ECC of 0.19, and an ESC of 2.72).Following industrial restructuring in the Chongming District, the proportion of the secondary industry's GDP decreased from 55.7% to 26.1%.This area is endowed with several national forest parks, important water sources, nature reserves, and significant wetlands and tidal flats.In subsequent development phases, the area should further optimize its industrial structure and actively promote the development of ecotourism and ecological agriculture.

Balanced Zone
The balanced zone is primarily concentrated in the northern part of the research area and the northern counties of Zhejiang Province, spatially surrounding the outer edge of the payment zone.It accounts for 24.13% of the total land area, showing sound economic development (with a GDP share of 27.93%).There is a good alignment between carbon emissions and carbon sink contributions (carbon emissions account for 23.65%, while carbon sinks account for 20.98%).The difference between the carbon emission-economic contribution capability and the ecology carrying capacity is relatively small (with an ECC of 1.18 and an ESC of 0.89).
(1) Balanced zone-optimized development zone.This region accounts for 7.89% of the total land area, spatially mostly scattered along the Yangtze River, the Beijing-Hangzhou Canal, and the Qiantang River.It consists of 29 county-level cities, primarily in Nanjing, Zhenjiang, Shanghai, Jiaxing, and Ningbo, exhibiting a relatively high level of economic development (accounting for 17.60% of the GDP and 13.48% of the carbon emissions), while also possessing adequate carbon sink capabilities (8.36% of carbon sinks, an ECC of 1.31, and an ESC of 0.62) Most of the area belongs to the central urban area.In the subsequent development, it should optimize the spatial layout of the land, eliminate backward production capacity, and actively develop green industries.Ecologically, efforts should focus on constructing wetlands along the river basin and enhancing coastal soil and water conservation capacity.
(2) Balanced zone-key development zone.This region is mainly located in Jiangsu Province.It consists of 12 county-level cities, spatially distributed in the northern part of the Payment zone-key development zone and along the southeastern coast of the CAYRDUA.This area accounts for 9.40% of the total land, characterized by a relatively low level of economic development, contributing less to the total carbon emissions (7.09% of GDP and 6.73% of carbon emissions), yet exhibiting good ecological carbon sequestration capacity (accounting for 7.63% of carbon sinks).The ecological carrying capacity slightly exceeds the carbon emission-economic contribution capacity (with an ECC of 1.05 and an ESC of 1.13).In subsequent development phases, this area should actively collaborate with surrounding first-tier cities to develop distinctive industries and enhance economic contribution capabilities.Simultaneously, efforts should be made to promote energy conservation and emission reduction among the enterprises, reducing greenhouse gas emissions.Enhancing forest coverage and quality and increasing carbon sink capacity are essential.
(3) Balanced zone-agricultural products' main production zone.This region accounts for 6.80% of the CAYRDUA, comprising six county-level cities, exhibiting good ecological carrying capacity (with an ESC of 1.48) but insufficient economic development (with an ECC of 0.95 and a GDP share of 3.18%).In subsequent development phases, efforts should be made to optimize the spatial structure, establish high-standard farmland, minimize the fragmentation of farmland, and enhance agricultural carbon sinks.By leveraging local soil characteristics and regional features, the cultivation of distinctive agricultural products should be promoted to elevate economic development levels.
(4) Balanced zone-key ecological functional zone.This region comprises only one county-level city, the Shengsi County in Zhoushan City.This area occupies 0.04% of the total land area and is characterized by a low level of economic development (with a GDP share of 0.07%), exhibiting a minor disparity between carbon emission-economic contribution capability and ecological carrying capacity (with an ECC of 0.70 and an ESC of 0.50).In subsequent development phases, there should be active promotion of marine economic industries and the construction of low-carbon and efficient smart ports.Simultaneously, the development of clean energy, such as wind energy, solar energy, and tidal energy, should be promoted to reduce carbon emissions.Ecologically, rational restoration and regular monitoring of coastal shorelines and wetlands should be conducted.

Compensated Zone
The compensated zone is concentrated in the southern part of the CAYRDUA.It accounts for about 41.84% of the land area, with a high ecological carrying capacity (the ECC is 1.45, the ESC is 3.98, and carbon sinks account for 50.31%), a relatively slow economic development (the share of GDP is 18.37%), and low carbon emission (the carbon emission only accounts for 12.65%).
(1) Compensated zones-optimized development zone.This region includes 14 county-level cities, primarily from Nanjing, Shanghai, Hangzhou, and Huzhou, accounting for 6.72% of the CAYRDUA.Despite relatively good economic development, the carbon emissions are low (with a GDP share of 9.76% and carbon emissions accounting for 4.44%, with an ECC of 2.20).It exhibits strong carbon sink capabilities (with carbon sinks accounting for 9.24% and an ESC of 2.08).This area is mainly situated in the hinterland of the CAYRDUA and can be divided into two parts.The first part comprises the old urban areas of first-tier cities, such as Shanghai, Hangzhou, and Nanjing, characterized by severe industrial hollowing.In subsequent development phases, efforts should be made to promote green and low-carbon lifestyles and emphasize urban greening actively.The second part comprises county-level cities primarily engaged in the secondary industry, mainly concentrated in cities such as Huzhou, Jiaxing, and Shaoxing in Zhejiang Province.This area has relatively favorable ecological conditions and high carbon sink capacity, with low intensity of land development.In subsequent development phases, attention should be paid to adjusting the industrial structure and promoting the transformation of traditional manufacturing industries.Strengthening the protection of vital water sources and existing ecological resources and exploring the path of ecological economy should also be emphasized.
(2) Compensated zones-key development zone.This region is mainly located along the southeastern coast of the CAYRDUA, which is entirely within Zhejiang Province.It includes 10 county-level cities, primarily in Ningbo, Shaoxing, Taizhou, and Zhoushan, accounting for 11.38% of the total land area.Despite low levels of economic development and carbon emissions (with a GDP share of 3.73% and carbon emissions accounting for 3.58%), this area demonstrates strong carbon sink capabilities (accounting for 15.13% of the total).It possesses good ecology carrying capacity (with an ECC of 1.04 and an ESC of 4.23).In the subsequent development phases, it is essential to capitalize on coastal advantages and foster unique marine industries.Additionally, exploring carbon emission markets and trading is crucial to identifying new avenues for economic growth.To address issues such as soil erosion and coastal pollution, ecology restoration measures should be implemented to enhance the region's carbon sink capacity further.
(3) Compensated zones-agricultural products' main production zone.This region is mainly located in Jiangsu Province, which is characterized by a flat terrain.It includes six county-level cities, primarily in Nanjing, Nantong, Taizhou, Yangzhou, and Jiaxing, accounting for 7.67% of the total land area.It contributes relatively less to carbon emissions and GDP (with a GDP share of 2.93% and carbon emissions accounting for 2.53%, with an ECC of 1.16) but exhibits good carbon sink functionality (contributing to 5.38% of carbon sinks and possessing an ESC of 2.13).In the subsequent development phases, efforts should be intensified to promote significant agricultural product development, develop distinctive brands, and foster specialized breeding techniques.Optimization of the ratio and usage of agricultural pesticides and fertilizers to promote low-carbon agriculture should be prioritized.Exploring new pathways that integrate specialty agriculture with tourism to drive local economic development should also be considered.
(4) Compensated zones-key ecological functional zone.This region is primarily located in the southwestern mountainous and hilly areas of the CAYRDUA, mainly including eight county-level cities, primarily in Hangzhou, Huzhou, Shaoxing, and Taizhou, accounting for 16.08% of the total land area.It exhibits strong ecological carrying capacity (with carbon sinks accounting for 20.56% and an ESC of 9.79).However, it has a low level of economic development (with a GDP share of 1.95% and an ECC of 0.9286).In the later stages of development, particular attention should be given to protecting water sources, such as Qiandao Lake and its surrounding wetlands, controlling the intensity of shoreline development, and addressing soil erosion in hilly areas.Economically, fostering collaboration with the surrounding ecological function zones is important for nurturing a distinctive tourism industry chain, consequently bolstering the ecological contribution to the economy.

Discussion
In this study, we calculated the carbon budget of the CAYRDUA by utilizing nighttime light data, land use data, energy consumption data, and other socio-economic data; analyzed the distribution patterns and evolutionary trends; explored carbon compensation zones at the county level; and proposed low-carbon development recommendations.The study found that carbon emissions in the CAYRDUA have slowed down since 2011, which is consistent with previous studies.Wei et al. suggested that carbon emissions in the YRDUA were around 1.6 billion tons after 2011 [59], whereas this study indicates that the total carbon emissions in the study area stabilized at around 1.2 billion tons.This difference can be attributed to the different scopes of urban agglomeration.The carbon sinks in the CAYRDUA estimated in this study were higher than those in other similar studies [60], which was because the carbon sinks estimated in this study comprised the mass of carbon dioxide rather than the mass of carbon.In addition to considering the aboveground carbon sinks, the carbon dioxide fixed by the underground root system and litter was also estimated in this study, which better reflected the changes in carbon budget in the CAYRDUA in each year of the long time series.The Guangdong-Hong Kong-Macao Greater Bay Area urban agglomeration, the Beijing-Tianjin-Hebei urban agglomeration, and the YRDUA are China's three major urban agglomerations.In terms of time, all three urban agglomerations showed decreasing trends in carbon budget after 2000 [58,61].In 2000, the carbon budget of the Guangdong-Hong Kong-Macao Greater Bay Area was better than that of the YRDUA, and the worst carbon budget was found in the Beijing-Tianjin-Hebei urban agglomeration.The carbon budget of the YRDUA was the worst after 2017.The carbon budget of the Guangdong-Hong Kong-Macao Greater Bay Area was better [58,61].This situation may be related to the industrial structure and economic scale, with Beijing-Tianjin-Hebei being a highly polluting and energy-consuming heavy industrial base [62].At the same time, the YRDUA has a sizeable economic scale, with a concentration of manufacturing industries and a high demand for energy.The Guangdong-Hong Kong-Macao Greater Bay Area has diversified industries, with high-value-added and tertiary industries working together as the pillar industries.It has a higher carbon emission efficiency than the other two city clusters, playing a greater role in low-carbon development [63].The industrial structure and carbon emission reduction experience of the Guangdong-Hong Kong-Macao Greater Bay Area urban agglomeration can provide a reference for the YRDUA and encourage the YRDUA to achieve a balance in the carbon budget.
From a spatial perspective, areas with high carbon emission contributions and economic development are often located in important water source areas, major river basins, and coastal areas, consistent with the findings of Xiong et al. [64].The payment areas are faceted on the northern shore of Taihu Lake, the southern shore of the Yangtze River, and the eastern coast, a phenomenon related to economic, policy, and geographic factors.From 2001 to 2016, the financial center of the Yangtze River Delta shifted northwestward to Suzhou, and the intensity of financial ties between Shanghai and the surrounding cities continued to increase [65].Leveraging the natural advantage of the Yangtze River waterway, Suzhou and Shanghai radiated to cities in southern Jiangsu Province, facilitating industrial transfer.Research has shown that the degree of synergy in the manufacturing industry between Suzhou, Wuxi, and Changzhou is higher than in other cities [66].In 2016, China issued the "Outline of Yangtze River Economic Belt Development", proposing to build development axes along the Yangtze River.It can be inferred that the trend of payment zones developing continuously to the west will continue.
The balanced zone-key development zone is mainly concentrated on the northern bank of the Yangtze River and the coastal areas surrounding the peripheral areas of the payment zone.The balanced zone-optimized development zone is scattered along both sides of the Beijing-Hangzhou Grand Canal and the coastal areas.From the perspective of ports, the Yangtze River, as China's largest river, has a higher port density in its basin compared to the Beijing-Hangzhou Grand Canal basin [67], benefiting the county-level cities of Yangzhou, Taizhou, and Nantong adjacent to the Yangtze River, which have a better development potential and become the main industrial receiving sites for first-tier cities.In addition to the development axis along the Yangtze River, coastal cities have also proposed a development axis along the coast.Building upon major maritime ports, such as the Shanghai Port and Ningbo-Zhoushan Port, industries and regions with high carbon emissions will continue to migrate and expand to the north and south of the CAYRDUA.
This paper puts forward a low-carbon development proposal for each carbon compensation zone.It plans to realize the proposal through cooperation with the government and undertaking government projects, suggesting that the government incorporate it into the existing planning instruments.In the process, considering the coordinated development of carbon emission reduction plans and the economy, new jobs and economic growth in localities will increase by developing the carbon trading market and new energy industries, thereby increasing the benefits of the green economy.This will contribute to realizing the "dual carbon" target.

Conclusions
The study computed carbon emissions and carbon sinks at the county level within the CAYRDUA from 2000 to 2020.The carbon budget was then calculated, and its spatiotemporal patterns and evolution trends were analyzed.A carbon compensation zoning index system was developed based on three attributes: socio-economy, ecology, and the scale of carbon emissions.Different types of carbon compensation zones were delineated using NRCA and the K-means algorithm.Subsequently, low-carbon development recommendations were proposed for each type of carbon compensation zone, considering regional characteristics in conjunction with the main functional zones.The main conclusions are as follows: (1) The growth of carbon emissions in the CAYRDUA can be delineated into three stages: a period of continuous rapid increase from 2000 to 2004, followed by a slowdown in growth from 2005 to 2011 and eventual stabilization of total carbon emissions after 2012.Spatially, high carbon emission areas have transitioned from a "core-periphery" pattern to a "Z-shaped" distribution, being primarily concentrated in the central region of the CAYR-DUA, with an expanding spatial extent over time.Carbon sinks in the CAYRDUA have exhibited a sustained declining trend, with a spatial distribution pattern characterized by a "high in the south, low in the north" trend.High-carbon sink regions are predominantly located in the hilly southern areas, while low-carbon sink regions are situated in the central and northern plains.
(2) In the CAYRDUA, a significant imbalance exists in the carbon budget in terms of quantity and spatial distribution.Quantitatively, carbon emissions far exceed carbon sinks, leading to a severe imbalance in the carbon budget, with the deficit continuously increasing.From 2000 to 2011, there has been a rapid decline in the carbon budget, followed by a relative stabilization of the total carbon budget after 2012.Spatially, carbon budget distribution demonstrates a "high in the south, low in the north" pattern, with significant regional disparities.Zhejiang Province exhibits better carbon budget capabilities than Jiangsu and Shanghai.From 2000 to 2020, the number of cities with a carbon budget surplus has decreased, becoming highly concentrated in the southwest.The spatial distribution of carbon budgets in the CAYRDUA exhibits obvious spatial autocorrelation.As time progresses, the degree of spatial autocorrelation decreases, ecological disturbances increase, and the area of high-high clustering decreases.
(3) This study constructed a carbon compensation zoning index system based on three attributes: socio-economy, ecology, and the scale of carbon emission.Subsequently, the CAYRDUA, consisting of 129 county-level cities, was classified into 43 payment zones, 48 balanced zones, and 38 compensated zone.These zones were overlaid with the main functional zoning, and considering the regional characteristics and development levels, appropriate low-carbon development recommendations were provided for each zone.The research findings hold significance in providing insights for developed regions and urban agglomerations in China to explore low-carbon pathways and formulate emission reduction policies.
In this paper, we discussed the evolution of carbon budget and carbon compensation zoning at the county scale and gave suggestions for low-carbon development by combining the characteristics of different carbon compensation zones.However, the shortcoming is that we did not analyze in depth the driving factors affecting the carbon budget.Moreover, fewer indicators were used in carbon compensation zoning, and more indicators with zoning characteristics will be added in the later study for more scientific and accurate carbon compensation zoning.

Figure 1 .
Figure 1.Study area: (a) The location of the CAYRDUA; (b) The CAYRDUA consisting of 16 prefecture-level cities.

Figure 1 .
Figure 1.Study area: (a) The location of the CAYRDUA; (b) The CAYRDUA consisting of 16 prefecture-level cities.

Figure 3 .
Figure 3.The relationship between county-level TDN of two datasets in 2012 and 2013.

Figure 3 .
Figure 3.The relationship between county-level TDN of two datasets in 2012 and 2013.

Figure 7 .
Figure 7. Changes in carbon budget and growth rate in the CAYRDUA, 2000-2020.

Figure 11 .
Figure 11.Distribution of carbon compensation zones.

Figure 11 .
Figure 11.Distribution of carbon compensation zones.

Table 2 .
The linear regression model of carbon emissions.

Table 2 .
The linear regression model of carbon emissions.

Table 5 .
The parameter of the global Moran's I.

Table 5 .
The parameter of the global Moran's I.

Table 6 .
Main indicators of carbon compensation zones in CAYRDUA.

Table 6 .
Main indicators of carbon compensation zones in CAYRDUA.