Long-Term Dynamics of Ecosystem Services and Their Inﬂuencing Factors in Ecologically Fragile Southwest China

: Southwest China has one of the largest karst landscapes on the Earth and an ecologically fragile environment. A better understanding of how ecosystem services function in karst areas helps ecological preservation and policy implementation. However, little effort has been made to evaluate the long-term dynamics of ecosystem services across Southwest China. This study systematically analyzed the spatio-temporal patterns and the values of three typical ecosystem service functions, including water conservation, carbon sequestration, and soil conservation, as well as the effects of precipitation and land use changes between 2000 and 2020. The results showed that water conservation exhibited an overall decrease from southeast to northwest, while soil conservation showed the opposite trend. The regions with an increasing trend in carbon sequestration were mainly distributed in Guizhou, Guangxi, and Sichuan. Compared to the year 2000, the percentage of water conservation and soil conservation decreased by 4.50% and 0.76%, respectively. However, carbon sequestration increased by 94.35%. The total value of ecosystem services in Southwest China showed a 90.00% increase in 2020 relative to 2000. Water conservation and carbon sequestration had a much closer correlation with precipitation. The impact of soil conservation was more signiﬁcantly inﬂuenced by land use changes.


Introduction
Ecosystem functions provided for humans are called ecosystem services, and they include water conservation, carbon sequestration, and soil conservation, which are essential for human survival and development [1,2]. The entire ecosystem value is estimated to be about tens of trillions of yuan per year, far exceeding traditional economic values [3]. However, increasing conflicts between economic development and environmental protection, as well as population growth and urbanization, are placing ecosystems under immense pressure [2,4]. Therefore, it is imperative to evaluate the functions and values of ecosystem services in order to quantify the impact of ecological conditions on human habitats and societal benefits, as well as to implement appropriate strategies for ecosystem management and ecological policies [5].
The karst landform, which covers approximately 15% of the Earth's surface, is predominantly utilized for agricultural and rural activities [6,7]. Southwest China is known as China's ecological barrier, owing to its rich biodiversity and vital ecosystem services [8,9]. Southwest China also owns the largest continuous distribution of karst areas in the country.
higher distribution in the west and lower in the east. The main soil types include red soil, purplish soil, and yellow soil [35]. Yunnan-Kweichow Plateau, Sichuan Basin, etc. [34]. The terrain is undulating and complex, with a higher distribution in the west and lower in the east. The main soil types include red soil, purplish soil, and yellow soil [35].

Precipitation Data
The precipitation data are a 1-km monthly precipitation dataset for China  and were obtained from the National Tibetan Plateau Data Center (http://data.tpdc.ac.cn, accessed on 6 August 2022). The coordination system of the data is WGS 84. This dataset was generally based on the 0.5° global climate dataset released by Climate Research Unit and the global high-resolution climate dataset released by World-Clim [36,37]. Additionally, 496 independent meteorological observation data points were utilized for verification, yielding reliable results [38,39]. The precipitation data were used to calculate water conservation and rainfall erosivity factor (R) in the Revised Universal Soil Loss Equation (RUSLE).

Soil Properties Data
The soil properties data, including sand, silt, clay, and organic carbon contents, were downloaded from National Tibetan Plateau Data Center (http://data.tpdc.ac.cn, accessed on 14 November 2022). The coordinate system of the data was WGS 84, with a spatial resolution of 30 arc seconds. In order to maintain the consistency of spatial resolution, the soil properties data were resampled to 1 km using ArcGIS 10.4.1. The source data for this dataset were obtained from 8595 soil profiles collected during the second soil census and 1:1 million Chinese soil maps [40]. The previous study demonstrated that this dataset was applied in Southwest China [41]. Therefore, we used soil properties data to calculate the soil erodibility factor (K) in the RUSLE.

Precipitation Data
The precipitation data are a 1-km monthly precipitation dataset for China  and were obtained from the National Tibetan Plateau Data Center (http://data. tpdc.ac.cn, accessed on 6 August 2022). The coordination system of the data is WGS 84. This dataset was generally based on the 0.5 • global climate dataset released by Climate Research Unit and the global high-resolution climate dataset released by WorldClim [36,37]. Additionally, 496 independent meteorological observation data points were utilized for verification, yielding reliable results [38,39]. The precipitation data were used to calculate water conservation and rainfall erosivity factor (R) in the Revised Universal Soil Loss Equation (RUSLE).

Soil Properties Data
The soil properties data, including sand, silt, clay, and organic carbon contents, were downloaded from National Tibetan Plateau Data Center (http://data.tpdc.ac.cn, accessed on 14 November 2022). The coordinate system of the data was WGS 84, with a spatial resolution of 30 arc seconds. In order to maintain the consistency of spatial resolution, the soil properties data were resampled to 1 km using ArcGIS 10.4.1. The source data for this dataset were obtained from 8595 soil profiles collected during the second soil census and 1:1 million Chinese soil maps [40]. The previous study demonstrated that this dataset was applied in Southwest China [41]. Therefore, we used soil properties data to calculate the soil erodibility factor (K) in the RUSLE.

Digital Elevation Model (DEM) Data
The DEM is produced by ALOS Global Digital Surface Model with a point spacing of 3pprox.. 30-m mesh and was downloaded from the Japan Aerospace Exploration Agency (https://earth.jaxa.jp/en/, accessed on 2 March 2023). The dataset has been collated with images acquired by the Advanced Land Observing Satellite (ALOS) called "DAICHI". The DEM data were generated based on the digital surface model dataset with a 5-m mesh version of the "World 3D Topographic Data" dataset of the "World 3D Topographic Data", which represented the most accurate global scale elevation data [42,43]. Its results in Southwest China are reliable [44]. In this study, the spatial resolution and coordinate system were adjusted to 1km and WGS 84, respectively. The DEM was applied to obtain the slope length and slope steepness factor (LS) in the RUSLE.

Land Use and Cover Data
The land use data were derived from China's National Land Use and Cover Change (CNLUCC) dataset during 2000-2020. The CNLUCC dataset was developed by the Data Center for Resources and Environment Sciences (https://www.resdc.cn, accessed on 28 October 2022), Chinese Academy of Sciences. It was based on visual interpretation of Landsat TM/ETM images with a spatial resolution of 1 km × 1 km. The coordinate system is WGS 84. Many researchers have proven that classification techniques above 92% were achieved through random sampling verification of survey sites in the field and assessments of the accuracy of the Kappa coefficient, which were applied in Southwest China [45]. The original CNLUCC dataset comprised 23 land use classes, based on a well-established hierarchical land cover classification system [46]. To facilitate the calculation of ecosystem service functions in this study, the 23 land use classes were further consolidated into six aggregated land use categories, namely construction land, farmland, forestland, grassland, unused land, and water. The land use data were used to calculate the crop management factor (C) and conservation support practice factor (P) in the RUSLE.

Net Ecosystem Productivity (NEP) Data
The NEP data are provided by the National Ecosystem Science Data Center, National Science & Technology Infrastructure of China (http://www.nesdc.org.cn, accessed on 29 August 2022). The dataset comprises simulated daily net ecosystem productivity over the globe during 1981-2019, which has been verified in some previous studies to be applicable in Southwest China [47]. In particular, the data for 2020 were not publicly available and were obtained directly from the research team. The coordinate system is WGS 84, with a spatial resolution of 0.072727 • , and was resampled to 1 km by ArcGIS 10.4.1. NEP was simulated using the process-based Boreal Ecosystem Productivity Simulator (BEPS) model, driven by meteorological data, remotely sensed vegetation parameters, nitrogen deposition atmospheric, and carbon dioxide (CO 2 ) concentration [47,48]. The amount of annual carbon sequestration was represented by NEP in the study.

Evapotranspiration Data
The evapotranspiration data were derived from A Big Earth Data Platform for Three Poles (https://data.tpdc.ac.cn, accessed on 29 February 2023). The data have a spatialtemporal resolution of daily and 500 m from 26 February 2000 to 31 December 2020. They have been widely applied in ecology studies, especially in Southwest China [49]. In the study, the resolution was resampled to 1 km using ArcGIS 10.4.1 and the lacking parts of 2000 were replaced by contemporaneous data from 2001. The actual evapotranspiration was the sum of the five types of evapotranspiration listed earlier and was used to calculate the water yield.

Water Conservation
The water conservation was calculated using the InVEST (Integrated Valuation of Ecosystem Services and Tradeoffs) Water Yield model. The model is founded on the principle of water balance, utilizing the difference between the precipitation within the grid cells and the actual losses from evapotranspiration to determine the water yield [50,51]. The formula is as follows: where Y(x) is the annual water yield of each image element x (mm); AET(x) is the annual actual evapotranspiration of each image element x (mm); P(x) is the annual precipitation of each image element x (mm). The shadow engineering method and the market price in the direct market method are commonly employed to evaluate the value of water conservation services [52,53]. The market price method was used in this study as follows [53]: where V x is the value of water yield (yuan); WR x is the amount of water yield in each raster (m 3 ); P is the price of water resources (yuan·m −3 ). This study used the price of water supply for residential domestic water from the waterworks in that year (with adjacent years substituted for years in which the price was not acquired).

Carbon Sequestration
NEP was used as a proxy for the carbon sequestration in the study area, with positive values indicating carbon sinks and negative values indicating carbon sources. The value of carbon sequestration was calculated by the market value method, which multiplied the physical amount of carbon by the carbon trading price [54]. The formula is as follows: where V c f is the value of carbon sequestration (yuan); Q CO 2 is the total amount of carbon sequestration (t); C c is the average price of the national carbon trading market transaction each year and was 64.10 yuan·t −1 before 2013, 26.38 yuan·t −1 in 2015, and 29.19 yuan·t −1 in 2020.

Soil Conservation
The reduction in soil erosion in the presence of vegetation cover and land management is defined as soil conservation [55]. Soil erosion was calculated based on the RUSLE, which mainly considered the effects of rainfall, topography, and vegetation on soil erosion [56,57]. The formula is as follows: where Q sr is soil erosion (t·hm −2 ·a −1 ); R is rainfall erosivity factor (MJ·mm·hm −2 ·a −1 ); K is soil erodibility factor (t·hm 2 ·h·hm −2 ·MJ −1 ·mm −1 ); LS is the slope length and slope steepness factor (dimensionless); C is the crop management factor (dimensionless); P is the conservation support practice factor (dimensionless).
(1) The rainfall erosivity factor (R) reflects the magnitude of the potential erosive capacity of the soil caused by rainfall [58]. The formula is as follows: where EI 30i is the average monthly rainfall erosivity factor (MJ·mm·hm −2 ·h −1 ); P i is the rainfall in month i (mm); P a is the average annual rainfall (mm); and the average annual rainfall erosivity factor is the sum of the average monthly rainfall erosivity force in a year.
(2) The soil erodibility factor (K) is an index to evaluate the sensitivity of soil erosion and is the amount of soil erosion per unit of rainfall erosion force [59]. The calculation was conducted in accordance with the "Technical Specification for Investigation and Assessment of National Ecological Status-Ecosystem Services Assessment" (HJ 1173-2021) issued by the Ministry of Ecology and Environment of the People's Republic of China. The formula is as follows: where K is the soil erodibility factor (t·hm 2 ·h·hm −2 ·MJ −1 ·mm −1 ); K EPIC is the soil erodibility factor calculated using the erosion-productivity evaluation model (t·hm 2 ·h·hm −2 ·MJ −1 · mm −1 ); m s , m silt , m c and orgC is sand, silt, clay, and organic carbon contents, respectively (%).
(3) The slope length and slope steepness factor (LS) accounts for the effect of the topography of the site on erosion [60]. It was calculated by Tool LS factor in SAGA GI [61]. The formula is as follows: where L is the slope length factor (dimensionless); S is the slope steepness factor (dimensionless); m is the slope length index (dimensionless); θ is the slope ( • ); λ is the slope length (m).
(4) Crop management factor (C) reflects the impact of cropping and management practices, as well as disturbances to the natural environment, including agriculture, productivity, crop rotation, growing season length, and sub-surface biomass, on soil erosion [60]. According to the relevant literature [62], it was classified into various land use types, as presented in Table 1. Table 1. Crop management factor (C) and conservation support practice factor (P) of different types of land use. (5) Conservation support practice factor (P) is the ratio of soil loss after specific practice measures to soil loss when planted downhill. It was assigned to different land use types based on the relevant literature [62], as shown in Table 1.

Land Use C P
Soil erosion has two primary impacts on humans: soil fertility loss and sedimentation [63]. Therefore, this study estimated the amount of soil conservation by considering both soil fertility conservation and sedimentation reduction values. The formula is as follows: where E is the value of soil conservation (yuan·hm −2 ); E a is the value of soil fertility loss (yuan·hm −2 ); E b is the value of reduced sedimentation (yuan·hm −2 ). Soil fertility maintenance encompasses both nutrient and organic matter values [64]. The former is determined by the market prices of chemical fertilizers, soil conservation practices, and nutrient content. The latter is estimated by converting the organic matter value into an equivalent amount of fuelwood, based on market prices. The formula is as follows: where E 1 is the nutrient value (yuan·hm −2 ); E 2 is the organic matter value (yuan·hm −2 ); Q is the annual soil conservation (t·hm −2 ); i is nitrogen, phosphorus, and potassium, respectively; m i is the average percentage content of nitrogen, phosphorus, and potassium in the soil, taken as 0.06%, 0.06%, and 1.40%, respectively; n i are the factors for converting alkaline nitrogen to urea, fast-acting phosphorus to calcium superphosphate, and fastacting potassium to potassium chloride, taken as 2.16, 4.07, and 1.92, respectively; p i is the national average market price of urea, calcium superphosphate, and potassium chloride in 2000, 2005, 2010, 2015, and 2020 (with adjacent years substituted for years when the prices are not available); D i is the soil organic carbon content (%); P S is the opportunity cost of fuelwood, which was taken as 195 yuan·t −1 in 2010 and before, 210 yuan·t −1 in 2015, and 140 yuan·t −1 in 2020 according to the average market price of fuelwood in that year; S is the coefficient of conversion of fuelwood into soil organic matter, which was taken as 0.5. The shadow engineering method is used to estimate the value of sedimentation by the cost required to dig and transport a unit volume of earth, which is replaced by the reservoir construction cost [65]. The formula is as follows: where Q is the annual soil conservation (t·hm −2 ); C is the construction cost per unit reservoir capacity (yuan·m −3 ), obtained by dividing the total reservoir investment (10 8 yuan) by the total reservoir construction capacity (

Statistical Analysis
The trend analysis method refers to analyzing the long-term sequence of relevant indicators to identify trends [66]. In this study, trend analysis was used to investigate the trends of three ecosystem service functions at both regional and pixel scales. The formula is as follows: where θ slope is the slope of the trend of a certain indicator; n is the total number of years of trend change, which was taken as 5; i is the year, and i took 1, 2, . . ., 5, representing 2000, 2005, . . ., 2020, respectively; y i is the value of the indicator in the year i. The larger the absolute value of θ slope , the more obvious the change of y i . The higher the absolute value of θ slope , the more obvious the change of y i . When θ slope < 0, it indicates a decreasing trend of an ecosystem service indicator; when θ slope > 0, it indicates an increasing trend of an ecosystem service indicator. If the regression coefficient passes the significance test (F test, p < 0.05), it indicates a significant ascending or descending trend [67].

Partial Correlation Analysis
Precipitation and land use changes have significant impacts on ecosystem service functions. However, linear correlation analysis may introduce bias when examining the relationship between influencing factors (precipitation and land use changes) and ecosystem service functions (water conservation, carbon sequestration, and soil conservation). Partial correlation analysis offers an advantage by excluding the impact of a third variable when two variables are simultaneously associated with it, allowing for the estimation of the correlation between the other two variables [66,68]. This method effectively overcomes the limitations of conventional correlation analysis. The formula for partial correlation analysis is as follows: where r xy·z is the partial correlation coefficient between variables x and y after fixing variable z; and r xy , r xz , and r yz are the correlation coefficients between variables x and y, x and z, and y and z, respectively. The correlation coefficient is calculated as follows: where r xy is the correlation coefficient for x and y, whose value ranges from −1 to 1. If one of the influencing factors passed the significance test (p < 0.05) but the other did not, the former was considered the dominant factor. If both factors passed the test, then both were considered dominant. However, if neither factor passed, it was assumed that these pixels were not influenced by either precipitation or land use changes.

Spatio-Temporal Variability of Ecosystem Service Functions
As illustrated in Figure 2, water yield exhibited a decreasing trend along the direction of southeast to the northwest from 2000 to 2020. Meanwhile, this study found that the marginal region of Yunnan province was apparently larger than the central zone for water yield, and West Sichuan has the least water yield. Strong spatial-temporal heterogeneity can be found for carbon sequestration in this region ( Figure 3). Overall, Guizhou, Guangxi, and Sichuan provinces showed a broad trend of rising carbon sequestration, whereas Chongqing and Yunnan provinces fluctuated throughout the study period. It was interesting to note that the spatial patterns of soil conservation contrasted with the water yield in Southwest China (Figure 4). Most areas had lower amounts of soil conservation, except the western part of Sichuan province. Table 2 illustrates that the dynamic trends of water conservation and carbon sequestration were similar, exhibiting a decrease, followed by an increase, and then a subsequent decrease after 2015. The minimum and maximum levels were observed in 2005 and 2015, respectively. Notably, the water conservation level in 2020 decreased by 4.50% compared to 2000, while carbon sequestration increased by 94.35%. Although the amount of soil conservation continued to increase after reaching its minimum in 2005, it decreased by 0.76% in 2020 compared to 2000. Particularly, the values for water conservation and soil conservation no longer decreased as indicated by the rates of change presented in Table 2. Instead, they increased by 137.33% and 77.71%, respectively. Carbon sequestration changed from the only improved function (Table 2) to the lowest change rate of 6.95% (Table 3). These results suggest that the impact of economic factors on the values of ecosystem service functions should not be neglected. As a consequence, the total value of ecosystem service functions in Southwest China exhibited a 90.00% increase in 2020 in comparison with 2000.         Table 2 illustrates that the dynamic trends of water conservation and carbon sequestration were similar, exhibiting a decrease, followed by an increase, and then a subsequent decrease after 2015. The minimum and maximum levels were observed in 2005 and 2015, respectively. Notably, the water conservation level in 2020 decreased by 4.50% compared to 2000, while carbon sequestration increased by 94.35%. Although the amount of soil conservation continued to increase after reaching its minimum in 2005, it decreased by 0.76% in 2020 compared to 2000. Particularly, the values for water conservation and soil conservation no longer decreased as indicated by the rates of change presented in Table 2. Instead, they increased by 137.33% and 77.71%, respectively. Carbon sequestration changed from the only improved function (Table 2) to the lowest change rate of 6.95% (Table 3). These results suggest that the impact of economic factors on the values of ecosystem service functions should not be neglected. As a consequence, the total value of ecosystem service functions in Southwest China exhibited a 90.00% increase in 2020 in comparison with 2000.   Figure 5 reveals that the trend analysis indicates a declining trend in the annual water yield of approximately 43.47% in Southwest China, with these regions being widely distributed in the western and northeastern parts of Southwest China. The minimum of the trend was observed in the south of Yunnan, whereas the southern parts of Chongqing and Guizhou were the areas where the maximum was distributed. Approximately 5.83% of the entire region, mainly located in West Sichuan and Southwest Yunnan, exhibited a significant decreasing trend. Conversely, there was only 4.62% (p < 0.05) in Southwest China, primarily in southern Chongqing, eastern Guangxi, and eastern Sichuan. Upon comparison with Figure 2, it is apparent that the regions displaying growth and decline trends in Figure 5 are roughly associated with areas characterized by high and low amounts of water conservation, respectively.  During 2000-2020, around 40.50% of the study area exhibited a negative trend in carbon sequestration, with notable increases observed in the north and east of Southwest China ( Figure 6). Meanwhile, approximately 59.50% of the area showed a positive trend, primarily in the southern region. Regions with significant trends were uniformly distributed, with 5.87% and 5.15% of the study area displaying a significant decreasing and increasing trend, respectively.

Long-Term Trend of Ecosystem Service Functions
The soil conservation trend was concentrated within the range of −20 t·hm −2 ·yr −1 to 20 t·hm −2 ·yr −1 , as shown in Figure 7. The areas exhibiting a declining trend were predominantly located in the west and south of the study area, covering roughly 44.21% of the overall area. Meanwhile, 55.79% of the study area demonstrated an increasing trend, with a highly significant increase (p < 0.05) observed in 7.52% of the region, mainly distributed in the central areas of Southwest China.

Effect of Environmental Factors
Spatial partial correlation analysis was utilized to investigate the response of three key ecosystem service functions to precipitation and land use changes (Figures 8 and 9). The results showed that water conservation was mainly controlled by precipitation, accounting for approximately 45.71% (p < 0.05) of the total area, with the highest impact observed in Guizhou, southern Chongqing, and northern Sichuan. Land use changes had relatively little impact, accounting for only 1.50% (p < 0.05) of the total area, while the combined effect of precipitation and land use changes accounted for 4.06% (p < 0.05), and no dominant factors accounted for 48.73% (p < 0.05).
In addition, carbon sequestration and soil conservation were less affected by precipitation, with only 4.19% (p < 0.05) and 2.37% (p < 0.05) of the total area, respectively. Regarding the influencing factors of carbon sequestration, land use changes had less impact than precipitation. Spatial distribution controlled by precipitation was mainly concentrated in the northeast of Yunnan. Land use changes had the largest impact on soil conservation, particularly in southern Chongqing and eastern Yunnan, accounting for approximately 4.04% (p < 0.05) of the total area. The area controlled by precipitation in soil conservation was distributed in western Yunnan, eastern Guangxi, and northern Chongqing. However, for both carbon sequestration and soil conservation, approximately 91.67% (p < 0.05) and 92.79% (p < 0.05) of the study area, respectively, were not affected by either precipitation or land use changes over the past two decades (Figures 8 and 9). Gray dots indicate areas that passed the significance test (p < 0.05). The bar graphs represent the frequency of the p value in the study area.

Effect of Environmental Factors
Spatial partial correlation analysis was utilized to investigate the response of three key ecosystem service functions to precipitation and land use changes (Figures 8 and 9). The results showed that water conservation was mainly controlled by precipitation, accounting for approximately 45.71% (p < 0.05) of the total area, with the highest impact observed in Guizhou, southern Chongqing, and northern Sichuan. Land use changes had relatively little impact, accounting for only 1.50% (p < 0.05) of the total area, while the combined effect of precipitation and land use changes accounted for 4.06% (p < 0.05), and no dominant factors accounted for 48.73% (p < 0.05).
In addition, carbon sequestration and soil conservation were less affected by precipitation, with only 4.19% (p < 0.05) and 2.37% (p < 0.05) of the total area, respectively. Regarding the influencing factors of carbon sequestration, land use changes had less impact than precipitation. Spatial distribution controlled by precipitation was mainly concentrated in the northeast of Yunnan. Land use changes had the largest impact on soil conservation, particularly in southern Chongqing and eastern Yunnan, accounting for approximately 4.04% (p < 0.05) of the total area. The area controlled by precipitation in soil conservation was distributed in western Yunnan, eastern Guangxi, and northern Chongqing. However, for both carbon sequestration and soil conservation, approximately 91.67% (p < 0.05) and 92.79% (p < 0.05) of the study area, respectively, were not affected by either precipitation or land use changes over the past two decades (Figures 8 and 9).   Figure 8 shows that precipitation was the dominant factor for water conserv which indicates that precipitation effectively accounted for the spatial and tem trends in water yield [27,69]. The increasing trends of water conservation in north chuan, southern Chongqing, and Guizhou were mainly controlled by precipitatio namics (Figures 5 and 8). Furthermore, it was noteworthy that extreme precipitatio to climate change under human activities was closely related to water yield [26]. Th tial distribution of water yield during anomalous extreme precipitation years exhib similar pattern of variation to that of total extreme precipitation [26].

Impact of Precipitation on Ecosystem Services
As shown in Table 2, NEP increased rapidly in the study area post-2010, con with the findings of Chen et al. [70], with more than 60% attributed to ecological eng ing with significant results. However, the enhancement of ecosystem carbon uptak ciency through ecological engineering was contingent upon the increases of precip and radiation [70]. It can be inferred that the significant increase in precipitation northeast of Yunnan and Sichuan during the last two decades was responsible for t nificant increase in NEP (Figures 6 and 8). In addition, as with most of the previous s [71,72], carbon sequestration exhibited a significant increase but decreased during 2010 (Table 2), probably due to drought and the reduction in solar radiation in Sout China [73]. Drought can lead to decreased plant production, increased mortality, a duced growth [74,75].
Overall, the annual precipitation and land use changes were not the predom factors in soil conservation (Figures 8 and 9). Areas of high amounts of soil conser were primarily concentrated in regions with lower susceptibility to extreme precip and higher fractional vegetation cover [26]. Previous studies [76][77][78][79][80] showed that th tial pattern of fractional vegetation coverage did not align well with that of soil con  Figure 8 shows that precipitation was the dominant factor for water conservation, which indicates that precipitation effectively accounted for the spatial and temporal trends in water yield [27,69]. The increasing trends of water conservation in northern Sichuan, southern Chongqing, and Guizhou were mainly controlled by precipitation dynamics (Figures 5 and 8). Furthermore, it was noteworthy that extreme precipitation due to climate change under human activities was closely related to water yield [26]. The spatial distribution of water yield during anomalous extreme precipitation years exhibited a similar pattern of variation to that of total extreme precipitation [26].

Impact of Precipitation on Ecosystem Services
As shown in Table 2, NEP increased rapidly in the study area post-2010, consistent with the findings of Chen et al. [70], with more than 60% attributed to ecological engineering with significant results. However, the enhancement of ecosystem carbon uptake efficiency through ecological engineering was contingent upon the increases of precipitation and radiation [70]. It can be inferred that the significant increase in precipitation in the northeast of Yunnan and Sichuan during the last two decades was responsible for the significant increase in NEP (Figures 6 and 8). In addition, as with most of the previous studies [71,72], carbon sequestration exhibited a significant increase but decreased during 2000-2010 (Table 2), probably due to drought and the reduction in solar radiation in Southwest China [73]. Drought can lead to decreased plant production, increased mortality, and reduced growth [74,75].
Overall, the annual precipitation and land use changes were not the predominant factors in soil conservation (Figures 8 and 9). Areas of high amounts of soil conservation were primarily concentrated in regions with lower susceptibility to extreme precipitation and higher fractional vegetation cover [26]. Previous studies [76][77][78][79][80] showed that the spatial pattern of fractional vegetation coverage did not align well with that of soil conservation in the study area, except for Sichuan. Based on the research by Eekhout et al. [81], it can be inferred that soil conservation was more affected by extreme precipitation than vegetation in Southwest China. Extreme precipitation enhanced soil erosion through the impact of raindrops on soil particles [82] and the infiltration of excess surface runoff that can promote gully erosion [83]. It should be noted that karst areas in Southwest China are characterized by the special above-ground and underground dual structure and porous media. Therefore, a large amount of soil can be lost during rainfall through bedrock fissures, waterfall holes, and solution channels [84], which was also an important factor affecting the spatial and temporal distribution of soil conservation in Southwest China.

Impact of Land Use Changes on Ecosystem Services
Land use changes affected the structure and function of ecosystems, which influenced the supply of ecosystem services [85,86]. The expansion of farmland and grassland [87] can simplify ecosystems [86,88] and alter local environmental characteristics [89], thereby affecting the provision of ecosystem services. Xiong et al. [90] suggested that the conversion of large amounts of unused land into farmland and grassland with higher ecological value coefficients was the main reason for the improvement of the ecosystem service functions in karst regions. However, urbanization promoted the conversion of a large amount of farmland into construction land, resulting in the deterioration of ecosystem service functions [91].
The IPCC (Intergovernmental Panel on Climate Change) identified land use changes as a major factor contributing to alterations in carbon sequestration in terrestrial ecosystems [45], owing to the varying carbon sequestration capacities of different land use types [29]. Urbanization-induced shifts from high to low carbon density land use types were identified as the second largest cause of carbon sequestration increase after CO 2 emissions from fossil fuel combustion [45]. Thus, it was plausible to hypothesize that the decline in carbon sequestration in the study area was linked to urbanization ( Figures 6, 8 and 10). Moreover, land use can impact water conservation function, such as infiltration, runoff, and evapotranspiration, through the water cycle [92,93]. In the south of Chongqing, a significant decrease in soil conservation was observed (Figure 7), which was attributed to land use changes as the dominant factor ( Figure 8). Li [94] reported that this region was one of the most economically developed areas in Southwest China, and the accelerated urbanization and relocation of villages and towns ( Figure 10) caused a decrease in soil conservation functions [95]. Additionally, ecologically well-conditioned land with high biomass and production had strong soil respiration, which resulted in high CO 2 concentrations [96]. In karst regions, the interaction among precipitation, soil CO 2 , and carbonate rocks led to dissolved inorganic carbon, which was one of the most crucial mechanisms of carbon sequestration [97].
Ecological policies can lead to land use changes. Since 2000, Southwest China has undertaken various ecological restoration projects such as the "Food for Green", "Natural Forest Protection Project", "Karst Rock Desertification Restoration Project", and other ecological restoration projects [66]. In 2008, the Rock Desertification Control Project (RDCP) was started in Southwest China, which significantly improved carbon sequestration and soil conservation in the region [14]. Numerous studies have demonstrated that the afforestation and conservation projects implemented in Southwest China have a significant positive impact on fractional vegetation coverage [14,66]. The increase in fractional vegetation coverage led to significant increases in carbon sequestration and soil conservation, but water conservation did not change significantly or even decreased [25,98]. In karst areas, irrational land use and reduced surface vegetation cover had an amplifying effect on soil conservation [99]. Ecological policies can lead to land use changes. Since 2000, Southwest China has undertaken various ecological restoration projects such as the "Food for Green", "Natural Forest Protection Project", "Karst Rock Desertification Restoration Project", and other ecological restoration projects [66]. In 2008, the Rock Desertification Control Project (RDCP) was started in Southwest China, which significantly improved carbon sequestration and soil conservation in the region [14]. Numerous studies have demonstrated that the afforestation and conservation projects implemented in Southwest China have a significant positive impact on fractional vegetation coverage [14,66]. The increase in fractional vegetation coverage led to significant increases in carbon sequestration and soil conservation, but water conservation did not change significantly or even decreased [25,98]. In karst areas, irrational land use and reduced surface vegetation cover had an amplifying effect on soil conservation [99].

Management Implications
Ecological restoration projects have made significant contributions to ecosystem services. Precipitation, as a representative of climatic conditions, is one of the most fundamental factors in enhancing ecosystem services [70]. Implementing ecological restoration projects in areas where the climate was predicted to be suitable for vegetation growth (such as due to precipitation) can reap half the benefits. Deforestation bans and afforestation were recognized as activities to mitigate climate change [70,100]. These measures led to an increase in forest cover, accelerated the accumulation of soil organic matter, improved soil biology, and reduced soil erosion [16]. However, they also brought a series of biophysical effects. The carbon cycle warming effects of deforestation were offset by net

Management Implications
Ecological restoration projects have made significant contributions to ecosystem services. Precipitation, as a representative of climatic conditions, is one of the most fundamental factors in enhancing ecosystem services [70]. Implementing ecological restoration projects in areas where the climate was predicted to be suitable for vegetation growth (such as due to precipitation) can reap half the benefits. Deforestation bans and afforestation were recognized as activities to mitigate climate change [70,100]. These measures led to an increase in forest cover, accelerated the accumulation of soil organic matter, improved soil biology, and reduced soil erosion [16]. However, they also brought a series of biophysical effects. The carbon cycle warming effects of deforestation were offset by net cooling associated with changes in albedo and evapotranspiration. Afforestation mitigates warming in the tropics, but probably not in subtropical, temperate, and boreal regions [101,102]. Planting large numbers of trees in non-forested areas might also not contribute to climate change mitigation [103]. Single-purpose land use and unrealistic land use policies often upset the balance between ecosystem services [5], resulting in reduced overall services [15]. Spatially oriented land use policies can serve to enhance the value of ecosystem services [5]. For example, after the 1998 floods in China, the overall ecosystem service functions were enhanced by 8.5% in the Poyang Lake area by implementing the project of returning fields to the lake for ecological restoration [104]. Local governments should take these into full consideration when formulating policies to mitigate climate change and tailor them to local conditions, rather than blindly pursuing the expansion of vegetation quantity and area.
However, compared to other studies, this study differs in the amounts and values of ecosystem service functions, probably due to the uncertainty of applied assessment datasets, models, and methods.

Conclusions
In this study, we quantified ecosystem services and their trends in Southwest China between 2000 and 2020. Further analysis was also conducted to explore the driving factors of ecosystem service changes. The main conclusions were as follows: (1) Over the past two decades, the spatial distribution of water conservation showed a decrease from southeast to northwest, the opposite of soil conservation. The regions with an increasing trend in carbon sequestration were concentrated in Guizhou, Guangxi, and Sichuan.
(2) Compared to 2000, the carbon sequestration level in 2020 increased by 94.35%, while water conservation and soil conservation decreased by 4.50% and 0.76%, respectively. The total value of three ecosystem services in Southwest China in 2020 increased by 90.00% compared to 2000.
(3) Next, 4.62% of Southwest China showed a significant increase in water conservation, which was mainly distributed in southern Chongqing, eastern Guangxi, and eastern Sichuan. For carbon storage, 5.83% and 5.15% of the study area showed significant downward and upward trends, respectively. Meanwhile, the trend of soil conservation exhibited a significant increase of 7.52%, which was primarily concentrated in the central areas of Southwest China.
(4) Water conservation was mainly controlled by precipitation, accounting for 45.71% of the total area. Extreme precipitation was also an important factor affecting water yield. Among the factors studied in this study, precipitation had a greater contribution to carbon sequestration (accounting for 4.19%). Annual precipitation only had a low impact (2.37%) on soil conservation. We suggested that extreme precipitation might have a greater impact on it through the impact of raindrops on soil particles.
(5) Land use changes contributed only 1.50% to water conservation functions through the water cycle. It contributed 2.78% to carbon storage, with the increase due to the expansion of vegetation area and its decrease related to urbanization. Moreover, land use changes had a tight negative impact (4.04%) on soil conservation due to accelerated urbanization and village relocation. (6) We suggested that appropriate ecological restoration projects should be projected based on climate predictions and that ecological policies should be tailored to local conditions, rather than blindly pursuing the expansion of vegetation quantity and area.