Decoupling the Effect of Climate and Land-Use Changes on Carbon Sequestration of Vegetation in Mideast Hunan Province, China

Urbanization and global climate change are two important global environmental phenomena in the 21st century. Human activities and climate changes usually increase the uncertainties of the ecosystem functions and structures and can greatly affect regional landscape patterns and the carbon cycle. Consequently, it is critical to understand how various climate and land-use changes influence the carbon dynamics at a regional scale. In this study, we quantitatively analyzed the spatial and temporal changes of net primary productivity (NPP) and the effects of climate factors and human disturbance factors (i.e., land-use changes) on the “Chang–Zhu–Tan” (CZT) urban agglomeration region from 2000 to 2015. The Carnegie–Ames–Stanford Approach (CASA) model was combined with spatially explicit land-use maps, monthly climate data, and MODIS NDVI images to simulate the carbon dynamics in the CZT area. Based on our six different scenarios, we also analyzed the relative roles of climate change and land-use change in total production. Our results indicated that the annual NPP of the study area maintained an upward trend by 7.31 gC•m−2•yr−1 between 2000 and 2015. At the same time, the average annual NPP was 628.99 gC•m−2 in the CZT area. We also found that the NPP was lower in the middle of the north region than in others. In addition, land-use changes could contribute to a positive effect on the total production in the study area by 3.42 T gC. Meanwhile, the effect of climate changes on the total production amounted to −1.44 T gC in the same region and period. Temperature and precipitation had negative effects on carbon sequestration from 2000 to 2015. As forest land constituted over 62.60% of the total land use from 2000 to 2015, the negative effect of carbon sequestration caused by urbanization could be ignored in the CZT area. Although climate and land-use changes had simultaneously positive and negative effects during the period 2000–2015, prioritizing the protection of existing forest land could contribute to increasing carbon sequestration and storage at the regional scale. Our study assists in understanding the impact of climate changes and land-use changes on carbon sequestration while providing a scientific basis for the rational and effective protection of the ecological environment in mid-east Hunan Province, China.


Introduction
At present, the increase in carbon emissions caused by human urbanization and climate changes has become a global economic and environmental problem [1]. Urbanization and global climate changes can greatly affect regional landscape patterns and the carbon cycle [2]. Therefore, decoupling and understanding the complex impacts of climate change and human-induced land-use change on carbon sequestration in the ecosystem is becoming more and more important [3]. In China, urban agglomerations have become very important boosters of urbanization [4]. Due to the rapid urban expansion and population growth, land-use change in these urban agglomerations is very common [5,6]. Quantitative analysis of the regional carbon sequestration capacity caused by vegetation dynamic change can improve the management of local land resources in the changing environment.
Net primary productivity (NPP) plays an important role in the global carbon cycle [7]. It is the key indicator revealing ecosystem carbon dynamics [8]. Previous studies mostly only focused on the spatial and temporal variations in NPP for a whole region [9,10]. However, the change of land type in the process of urbanization depends on the special conditions and specific environmental characteristics of different cities [11,12]. Consequently, the effects of various land-use changes on carbon sequestration for different cities in the same region are ignored. Moreover, cities close to each other often have similar climatic characteristics, so few studies pay attention to whether the effect of similar climate conditions on the NPP of cities with different land-use policies will be different. In this context, it is necessary to analyze the impact of climate change on NPP for different cities in the same region.
As it is impossible to make a comprehensive measurement of NPP on a large scale, relevant models to indirectly estimate NPP have been widely applied in regional or global studies [13]. At present, many NPP estimation models, which are based on remote sensing methods, are widely used in global and regional NPP evaluation. For example, the TRIPLEX-Flux [12] and TRIPLEX-GHG [14] models, which are based on ecosystem process design, are used to simulate and predict the carbon dynamics of forest ecosystems and wetland ecosystems, respectively. CITYgreen [15] and i-tree [16] are two useful urban forest-specific models for simulating the carbon dynamics of urban green space ecosystems. The Carnegie-Ames-Stanford Approach (CASA) model, which is a process-based and light-use efficiency model, has been applied to estimate the carbon dynamics of terrestrial ecosystems [17]. Comparing these models, the CASA model is easier to master, and its parameters are easy to obtain by remote sensing technology. Therefore, the CASA model has been widely used in NPP estimation in different regions and scales in China [18][19][20].
The "Chang-Zhu-Tan" (CZT) [21] urban agglomeration is an important part of the urban agglomeration in mid-east Hunan Province, China [4]. Although the structures of these three cities are relatively close, they maintain the largest urban ecological green center (i.e., 528.32 km 2 ) between them [22]. The rapid expansion of urbanization and the implementation of ecological protection policies in the last 20 years provide unique research value when assessing the relationship between urbanization and carbon balance. In this study, we combined spatially explicit land-use maps, monthly climate data, and MODIS NDVI images with the CASA model to estimate the annual NPP of the CZT urban agglomeration between 2000 and 2015. Specifically, the objectives of this study were to (1) simulate the vegetation NPP and characterize the spatiotemporal variation in CZT from 2000 to 2015; (2) reveal the climate (i.e., temperature, precipitation, and solar radiation) and land-use (i.e., afforestation and urbanization) changes and their link with NPP variation in the CZT area; and (3) evaluate the relative effects of climate and land-use changes on the total production of each city based on our six different scenarios. This study will provide some fundamental support for the management of local landscapes under the background of climate and land-use changes and contribute to improving the understanding of the driving mechanism of urban agglomeration development and ecosystem carbon dynamics in subtropical humid areas.

Study Area
The CZT urban agglomeration has an area of about 28,121 km 2 and lies between 26 • 04 -28 • 65 N and 111 • 90 -114 • 21 in China. There are three cities in this study area: Changsha, Zhuzhou, and Xiangtan, with a total population of 12,000,000-13,000,000. The distance between the core urban areas of each city is no more than 40 km. The study area contains many natural geographical elements, such as rivers, hills, forest land, grassland, wetland, and cropland, so its ecosystem structure is complex [23,24]. The region has a midsubtropical monsoon climate, with a historical annual mean temperature of 16.8-17.3 • C and an annual accumulated temperature of about 5457 • C. The maximum monthly mean temperature is 32.0 • C in July while the minimum monthly mean temperature is 5.0 • C in January. Meanwhile, in this region, rain and heat occur in the same period. The annual precipitation ranges from 936.4 to 1954.2 mm, with a mean annual precipitation of 1416.4 mm. The soil type is mainly well-drained clay loam red soil and paddy soil and is classified as an Alliti-Udic Ferrosol in the Chinese Soil Taxonomy (CST) [25], which corresponds to Acrisol in the World Reference Base for Soil Resources [26]. Subtropical evergreen broadleaved forest is the climax vegetation in CZT [27].

Data
The climate data for 2000-2015 were obtained from the National Tibetan Plateau Third Pole Environment Data Center (TPDC) [28]. They include monthly mean temperature, monthly mean precipitation, and monthly mean solar radiation. Land-use maps of the study area, derived from Landsat TM/ETM+/OLI imageries in 2000, 2005, 2010, and 2015, were obtained from the Resource and Environment Science and Data Center (https://www.resdc.cn/, accessed on 07 November 2021). The land use is divided into nine classes: paddy field, dry land, arboreal forest, shrub, open forest land, grassland, water, built-up land, and bare land (see Table 1) ( Figure 1). Moreover, monthly NDVI values were obtained using the maximum-value compositing algorithm based on time-series MODIS NDVI images at the pixel level.

Scenarios
In this study, we calculated the total production (TP) in an area based on Equation (1):

Scenarios
In this study, we calculated the total production (TP) in an area based on Equation (1): To decouple and analyze the effects of climate changes (i.e., temperature, precipitation, and solar radiation) and land-use changes (i.e., afforestation, urbanization, and others) on the total production in CZT, we designed six scenarios to estimate the potential TP in 2000 and 2015 (see Table S1) [9].
In this study, we kept the climate conditions (i.e., precipitation, temperature, and solar radiation) at the 2000 level in scenario A, and then we calculated the potential TP in 2000, indicated as TP 2000 . We also estimated TP A based on the land-use map and NDVI images from 2015. The effects of land-use change (∆LUCC) could be calculated using Equation (2): Forests 2021, 12, 1573

of 14
The overall effects on TP (∆ALL) were calculated based on the difference between the actual TP of 2000 and 2015 (see Equation (3)): Based on the above information, we calculated the effects of climate changes (∆Climate) using Equation (4): To decouple the effects of different climate conditions, we designed three different scenarios: keeping the precipitation at the 2000 level with no change (scenario B), keeping the temperature at the 2000 level with no change (scenario C), and keeping the solar radiation at the 2000 level with no change (scenario D). Based on these scenarios, we estimated the effects of precipitation (∆Precipitation), temperature (∆Temperature), and solar radiation (∆Radiation) on TP, respectively (Table S1) [9].
To decouple the effects of different land-use changes, we used Equation (5) as follows: where ∆Afforestation represents the effects of afforestation on TP, estimated based on scenario E (Table S1); ∆Urbanization represents the effects of urbanization on TP, estimated based on scenario F (Table S1); and ∆Others represents the effects of other land-use changes.

Modeling 2.3.1. CASA Model
The CASA model is based on process development and is driven by light-use efficiency [29]. It was used to estimate the NPP from absorbed photosynthetic active solar radiation (APAR) and light-use efficiency (ε) [17]. APAR was estimated using Equation (6): where SOL is the total solar radiation; FPAR is the fraction of the photosynthetically active solar radiation absorbed by vegetation canopy; t is time; x is the spatial location. The constant 0.5 indicates the proportion of effective solar radiation available to vegetation from the total solar radiation [30]. In this study, FPAR was calculated using Equation (7): where FPAR(x, t) NDVI is the FPAR calculated by NDVI in pixel x and month t (see Equation (8)); SR is the simple ratio of NDVI and is calculated by NDVI(x, t) (see Equation (9)); FPAR(x, t) SR is the FPAR calculated by SR in pixel x and month t (see Equation (10)) [31].
where FPAR max and FPAR min are 0.950 and 0.001, respectively, in this study; NDVI i,min is the minimum value of the NDVI for the land-use type i in month t; NDVI i,max is the maximum value of the NDVI for the land-use type i in month t; SR i,min is the minimum SR value for the land-use type i in month t; SR i,max is the maximum SR value for the land-use type i in month t; and NDVI(x, t) is the NDVI in pixel x and month t [30].
Light-use efficiency (ε), which is the efficiency of vegetation to convert absorbed photosynthetic effective solar radiation into organic carbon, is a key parameter for calculating NPP [20]. It can be affected by various environmental conditions (e.g., temperature and moisture). In the CASA model, ε is estimated using Equation (11): where T ε1 (x, t) and T ε2 (x, t) are temperature stress coefficients, which indicate the reduction in light-use efficiency caused by temperature; W ε (x, t) is the moisture stress coefficient, which reflects the reduction in light-use efficiency caused by moisture; and ε max is the maximum light-use efficiency and is determined by the empirical method [20]. For details of the above parameters, see Table S2.

Coefficient of Variation
In this study, the stability of the interannual variation of the NDVI in the CZT urban agglomeration is represented as the coefficient of variation (C v ) from 2000 to 2015. C v is a useful indicator for analyzing the interannual variability of an ecosystem. It is a statistical index, which means the ratio of a variable relative to its average, and is estimated using Equation (12): To illuminate and analyze the spatial dynamics of the NPP in CZT, we evaluated the mean annual NPP and its linear trend between 2000 and 2015. According to the linear trends of annual NPP and their significance, we defined three NPP change types (i.e., no change, increase, and decrease). Our simulations demonstrated that the mean annual NPP presented a distinct heterogeneity in the study areas (Figure 3a). The areas with a very low annual NPP were mainly distributed in the middle of the north region, which also comprises the core built-up area of the three cities (Figure 3a). The maximum increasing rate was 36.52 gC•m −2 •yr −1 in the middle of the east of the study areas, while the maximum decreasing rate was −24.33 gC•m −2 •yr −1 in the north-central part of the study areas (Figure 3b). We also found that the mean annual NPP indicated a decreasing trend around the built-up areas during 2000-2015. In addition, the mean annual NPP changed significantly in the central region and the middle of the east region in the study area, as we found that the p-values of the t-test for the modeled slope of these pixels were less than 0.05 in these regions (Figure 3c).  Our simulations demonstrated that the mean annual NPP presented a distinct heterogeneity in the study areas (Figure 3a). The areas with a very low annual NPP were mainly distributed in the middle of the north region, which also comprises the core built-up area of the three cities (Figure 3a). The maximum increasing rate was 36.52 gC•m −2 •yr −1 in the middle of the east of the study areas, while the maximum decreasing rate was −24.33 gC•m −2 •yr −1 in the north-central part of the study areas (Figure 3b). We also found that the mean annual NPP indicated a decreasing trend around the built-up areas during 2000-2015. In addition, the mean annual NPP changed significantly in the central region and the middle of the east region in the study area, as we found that the p-values of the t-test for the modeled slope of these pixels were less than 0.05 in these regions (Figure 3c). The simulated results show that the annual NPP had an increasing trend in 37.51% of the study area and a decreasing trend in 11.43% of the study area (Figure 3d). The decrease areas (i.e., the red color in Figure 3d) were mainly distributed in the middle of the north region. Meanwhile, the increase areas (i.e., the green color in Figure 3d) were mainly distributed in the middle of the east region.

The Effects of Climate Changes on Total Production
In the study area, temperature, precipitation, and solar radiation indicated different linear change rates of 0.035 • C•yr −1 , −0.241 mm•yr −1 , and 3.886 MJ•m −2 •yr −1 , respectively, from 2000 to 2015. Figure 4 shows that the spatial distribution characteristics of them are quite different. The temperature had an increasing trend in most regions, except the southern region of CZT (Figure 4a). The increased precipitation was mainly concentrated in the northeast, while the southern area had decreased precipitation (Figure 4b). For solar radiation, the northeast of the study area showed a decreasing trend, but the western and southern regions showed an increasing trend from 2000 to 2015 (Figure 4c). Based on the scatters and fitted linear correlations analysis, we found that the annual NPP positively correlated with temperature ( Figure 5a) and solar radiation (Figure 5c) but negatively correlated with precipitation ( Figure 5b). The correlation coefficient between NPP and temperature was lower in the north central part of the study area, but a higher correlation coefficient was sparsely distributed in the central parts of the CZT area (Figure 6a). For precipitation, higher correlation coefficients were distributed in the northeast region ( Figure 6b). However, the opposite distribution could be observed for the correlation coefficient between NPP and solar radiation in the study area (Figure 6c). The simulated results show that the annual NPP had an increasing trend in 37.51% of the study area and a decreasing trend in 11.43% of the study area (Figure 3d). The decrease areas (i.e., the red color in Figure 3d) were mainly distributed in the middle of the north region. Meanwhile, the increase areas (i.e., the green color in Figure 3d) were mainly distributed in the middle of the east region.

The Effects of Climate Changes on Total Production
In the study area, temperature, precipitation, and solar radiation indicated different linear change rates of 0.035 °C•yr −1 , −0.241 mm•yr −1 , and 3.886 MJ•m −2 •yr −1 , respectively, from 2000 to 2015. Figure 4 shows that the spatial distribution characteristics of them are quite different. The temperature had an increasing trend in most regions, except the southern region of CZT (Figure 4a). The increased precipitation was mainly concentrated in the scatters and fitted linear correlations analysis, we found that the annual NPP positively correlated with temperature ( Figure 5a) and solar radiation (Figure 5c) but negatively correlated with precipitation (Figure 5b). The correlation coefficient between NPP and temperature was lower in the north central part of the study area, but a higher correlation coefficient was sparsely distributed in the central parts of the CZT area (Figure 6a). For precipitation, higher correlation coefficients were distributed in the northeast region (Figure 6b). However, the opposite distribution could be observed for the correlation coefficient between NPP and solar radiation in the study area (Figure 6c).   related with precipitation (Figure 5b). The correlation coefficient between NPP and temperature was lower in the north central part of the study area, but a higher correlation coefficient was sparsely distributed in the central parts of the CZT area (Figure 6a). For precipitation, higher correlation coefficients were distributed in the northeast region (Figure 6b). However, the opposite distribution could be observed for the correlation coefficient between NPP and solar radiation in the study area (Figure 6c).    Based on our results, we found that the effect of climate changes on the total production was −1.44 T gC in the whole area of the CZT urban agglomeration between 2000 and 2015 ( Table 2). Our scenario simulations showed that Changsha, Zhuzhou, and Xiangtan faced negative effects of climate changes on the total production by −0.50, −0.72, and −0.22 T gC from 2000 to 2015, respectively. We also found that precipitation was the main meteorological driving factor influencing the NPP of the CZT areas, followed by solar radia- Based on our results, we found that the effect of climate changes on the total production was −1.44 T gC in the whole area of the CZT urban agglomeration between 2000 and 2015 ( Table 2). Our scenario simulations showed that Changsha, Zhuzhou, and Xiangtan faced negative effects of climate changes on the total production by −0.50, −0.72, and −0.22 T gC from 2000 to 2015, respectively. We also found that precipitation was the main meteorological driving factor influencing the NPP of the CZT areas, followed by solar radiation and temperature (Table 3). Solar radiation had a positive effect on the total production by 1.59 T gC, but precipitation and temperature had negative effects on the total production by −1.62 and −1.41 T gC, respectively, during the period 2000-2015 (Table 2). In the same period, we found that the total production of Zhuzhou was the most affected by temperature (i.e., −0.73 T gC), precipitation (i.e., −0.76 T gC), and solar radiation (i.e., 0.77 T gC) ( Table 2). Table 2. The effects of climate changes, precipitation, temperature, and radiation on the total production (T gC) in the CZT urban agglomeration from 2000 to 2015. ∆Climate, the effect of climate changes; ∆Precipitation, the effect of precipitation changes; ∆Temperature, the effect of temperature changes; ∆Radiation, the effect of radiation changes.  Table 3. The effects of land-use changes on the total production (T gC) in the CZT urban agglomeration from 2000 to 2015. ∆LUCC, the effect of all land-use changes; ∆Afforestation, the effect of nonforest land changes to forest land; ∆Urbanization, the effect of nonconstruction land changes to construction land; ∆Others, the effect of other land-use changes.

The Effects of Land-Use Changes on Total Production
Our simulations indicated that the effect of land-use changes on the total production was 3.42 T gC in the whole area of the CZT urban agglomeration between 2000 and 2015 (Table 3). We found that the afforestation area was 163.72 km 2 in the study areas between 2000 and 2015, and it had a positive effect on the total production by 0.05 T gC (Table 3). At the same time, the area of built-up land increased from 650.34 to 1387.83 km 2 in CZT ( Table 1). The simulated results show that the negative effect of urbanization (i.e., −0.05 T gC) was almost equal to the positive effect of urbanization on the total production in CZT during the period 2000-2015. In fact, urbanization also had slight effects on the total production for each city in the CZT region from 2000 to 2015. The effects of urbanization on the total production were −0.03 T gC in Changsha, −0.01 T gC in Xiangtan, and −0.01 T gC in Zhuzhou (Table 3).

Validation
As it is very difficult to monitor NPP in a large area, it is often difficult to verify the accuracy of simulation results on a large scale [10,32]. Here, we compare and verify our simulation results with relevant studies. Zhang and Zeng [33] reported that the annual NPP was 260.88-323.01 gC•m −2 •yr −1 in the core urban areas of the CZT urban agglomeration in 2015, and our results on the annual NPP in the same region were basically consistent.
The annual NPP of the core urban areas was 250.33-301.22 gC•m −2 •yr −1 in this study in 2015. Although uncertainties caused by the different spatiotemporal scales and land-use type changes in different regions always exist, other previous studies demonstrated that the NPP estimated by the CASA model combined with MODIS datasets was reliable in China [10,18,34].

Response of Carbon Dynamics to Climate Changes
Climatic conditions have different effects on the NPP of terrestrial ecosystems of different regions [35,36]. In this study, climate changes had a negative effect on the total production of the CZT area from 2000 to 2015. Our investigation demonstrated that this negative effect mainly came from the increasing trend of temperature (0.035 • C•yr −1 ) and the decreasing trend of precipitation (−0.241 mm•yr −1 ) in the study areas. Lauerwald et al. [37] pointed out that the NPP of a forest ecosystem could decrease with an increase in temperature if there is enough precipitation in the region, and, as mentioned earlier, the CZT urban agglomeration is located in a humid area with a subtropical monsoon climate. Although precipitation showed a decreasing trend from 2000 to 2015, the average annual precipitation was still more than 1500 mm•yr −1 . Another reason might be the increase in soil respiration (i.e., autotrophic and heterotrophic soil respiration) caused by the increase in temperature, which leads to carbon losses [38].
The mechanism behind the effect of climate factors on regional carbon dynamics is complicated [10]. In this study, although the effect of precipitation on total production was greater than that of both temperature and solar radiation, there was little difference between them. Our results are similar to those of Nemani et al. [39], who reported that about 40% of vegetation growth is mainly limited by water resources, while 33% and 27% is limited by temperature and solar radiation, respectively.

Response of Carbon Dynamics to Land-Use Changes
Land-use change is one of the most important human disturbances to regional carbon dynamics [40]. Our observation indicates that urban expansions and ecological land losses (i.e., cropland, forest land, and grassland) had a highly similar spatial distribution in the CZT areas (Figure 1). We found that land-use change in CZT was driven by urbanization expansion and the influence of ecological protection from 2000 to 2015. In this process, the increased built-up land (i.e., 737.49 km 2 ) mainly came from cropland (i.e., 399.02 km 2 ) and forest land (i.e., 351.79 km 2 ) from 2000 to 2015 (Table 4). This result is consistent with that of Zhu and He [4]. All these caused a reduction in productivity (i.e., −0.05 T gC) in the CZT areas between 2000 and 2015. Moreover, because ∆Others was 3.42 T gC in our study areas from 2000 to 2015, the effects of other land-use changes played a key positive role in carbon sequestration in the study areas. The cities most and least affected by ∆Others were Zhuzhou (i.e., 1.87 T gC) and Xiangtan (i.e., 0.51 T gC), respectively (Table 3). We think that the main source of ∆Others is the natural growth brought on by the protection of original vegetation. The natural growth of original forest land has been the main driver of the increasing trend of NPP in the CZT areas in the past decade. During this period, the afforested area was 163.72 km 2 , while 489.89 km 2 of forest land was converted to other land types. However, the majority of forest land (i.e., 17,452.78 km 2 ) was protected, with no land-use change. From the spatial distribution map of NPP change (Figure 3d) and the spatially explicit land-use maps (Figure 1), we also found that the areas with increased NPP were highly coincident with well-protected forest areas (such as urban green core regions, forest parks, etc.). These areas may benefit from the impact of the recent ecological restoration project, the "overall planning of ecological core green area of CZT urban agglomeration" [22]. In addition, our observation showed that about 37.51% of the study area had an increasing trend of annual NPP from 2000 to 2015 (Figure 3d), while 51.06% of the study area had no change trend for annual NPP (Figure 3d). These results demonstrate that the overall status of land vegetation in CZT became better during the period 2000-2015.

Limitation and Future Work
Some uncertainties in this study must also be considered. The first uncertainty comes from the downscaling of MODIS NDVI products, because the NPP simulated by the CASA model needs more precise NDVI data to improve the simulation accuracy [41]. Secondly, the input of climate data, which involved a fusion of remote sensing products, reanalysis datasets, and in situ observation data from weather stations, might bring uncertainty to the data accuracy [28]. Thirdly, we did not consider the uncertainty from natural disturbances (e.g., fire, insects, and wind). Ignoring these natural disturbances may lead to overestimation of the NPP at the regional scale [42,43]. In addition, some previous studies indicated that the increase in atmospheric carbon dioxide concentration and the acceleration of nitrogen deposition may affect the simulation of regional NPP [34,[44][45][46]. Overall, the impacts of the above shortcomings on the estimation of NPP certainly need to be further considered and addressed in future work.

Conclusions
In this study, the CASA model was used to estimate the annual NPP in the CZT urban agglomeration from 2000 to 2015. Although urbanization has obviously occurred in the past sixteen years, the annual NPP in this region has maintained an upward trend. Zhuzhou had the highest increasing trend of NPP at 9.70 gC•m −2 •yr −1 from 2000 to 2015. The maximum annual NPP (i.e., 780.45 gC•m −2 ) occurred in Zhuzhou in 2014. Meanwhile, we found that the mean annual NPP presented distinct heterogeneity in the study areas. The mean annual NPP was lower in the middle of the north region than elsewhere, which indicated a decreasing trend around the built-up areas.
Moreover, this study revealed the effects of climate changes and land-use changes on the carbon sequestration capacity in CZT between 2000 and 2015. According to our scenario simulations, land-use changes had a positive effect on the total production of the CZT area of 3.42 T gC-greater than the negative effects from climate changes of −1.44 T gC during the period 2000-2015. We found that both temperature and precipitation had significantly negative effects on the total production in the study areas. Contrarily, solar radiation had a positive effect on the total production in the same period. Our results also showed that the positive effect of afforestation offset the negative effect of urbanization in CZT from 2000 to 2015. As the proportion of forest land in the region remained relatively high (i.e., >62.60%), the negative effect on carbon sequestration of vegetation caused by urbanization and climate changes could be ignored in the CZT urban agglomeration from 2000 to 2015.
Our study suggests that protecting existing forests and vegetation is the most effective way to increase regional carbon sink.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10.3 390/f12111573/s1, Table S1: The information of six different scenarios in the "Chang-Zhu-Tan" urban agglomeration from 2000 to 2015; Table S2: Parameters of the CASA model for different land-use types in the CZT urban agglomeration.