Next Article in Journal
New Observations of the Meteorological Conditions Associated with Particulate Matter Air Pollution Episodes in Santiago, Chile
Next Article in Special Issue
Sensitivity Analysis of the Land Surface Characteristic Parameters in Different Climatic Regions of the Loess Plateau
Previous Article in Journal
Vectorial EM Propagation Governed by the 3D Stochastic Maxwell Vector Wave Equation in Stratified Layers
Previous Article in Special Issue
Composition, Distribution, and Attribution of Hydrochemistry in Drainage Systems in the North of Tianshan Mountains, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial and Temporal Evolution Characteristics of Water Conservation in the Three-Rivers Headwater Region and the Driving Factors over the Past 30 Years

1
Key Laboratory of Land Surface Pattern and Simulation, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Atmosphere 2023, 14(9), 1453; https://doi.org/10.3390/atmos14091453
Submission received: 14 August 2023 / Revised: 11 September 2023 / Accepted: 12 September 2023 / Published: 18 September 2023
(This article belongs to the Special Issue Land-Atmosphere Interactions over the Tibetan Plateau)

Abstract

:
The Three-Rivers Headwater Region (TRHR), located in the hinterland of the Tibetan Plateau, serves as the “Water Tower of China”, providing vital water conservation (WC) services. Understanding the variations in WC is crucial for locally tailored efforts to adapt to climate change. This study improves the Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) water yield model by integrating long-term time series of vegetation data, emphasizing the role of interannual vegetation variation. This study also analyzes the influences of various factors on WC variations. The results show a significant increase in WC from 1991 to 2020 (1.4 mm/yr, p < 0.05), with 78.17% of the TRHR showing improvement. Precipitation is the primary factor driving the interannual variations in WC. Moreover, distinct interactions play dominant roles in WC across different eco-geographical regions. In the north-central and western areas, the interaction between annual precipitation and potential evapotranspiration has the highest influence. Conversely, the interaction between annual precipitation and vegetation has the greatest impact in the eastern and central-southern areas. This study provides valuable insights into the complex interactions between the land and atmosphere of the TRHR, which are crucial for enhancing the stability of the ecosystem.

1. Introduction

The Three-Rivers Headwater Region (TRHR), located in the hinterland of the Tibetan Plateau, plays a key role in ensuring ecosystem security [1]. Renowned as the “Water Tower of China”, the TRHR is the birthplace of the Yellow River, Yangtze River, and Lancang River and is essential for water storage and conservation [2]. Water conservation (WC) is a crucial ecosystem regulation service, encompassing the interception and maintenance of precipitation by various components of the ecosystem, including vegetation canopy, litter, soil, lakes, water reservoirs, and more [3]. The WC services of the TRHR satisfy the water needs of the local ecosystem and provide water resources for external and downstream areas [4,5]. However, the TRHR experiences an uneven and relatively low precipitation distribution. Rising temperatures also contribute to the degradation of frozen soil and increased evapotranspiration, which have a negative impact on WC [6,7]. Consequently, accurately assessing the WC in the TRHR and determining its driving factors is crucial to understanding of the ecosystem’s ability to maintain the stability and sustainability of water resources.
Previous studies have revealed an increasing trend in WC within the TRHR since the 1990s [8,9,10]. The increase in WC results from the complex interplay between ecosystem and hydrological processes, influenced by various factors such as ecosystem type, surface characteristics, climate conditions, and human activities [11]. However, differences exist in understanding the attribution of WC variations in this region. Studies suggest the rise in WC can be attributed to factors such as precipitation [3], potential evapotranspiration [12], the combined effect of increased precipitation and vegetation growth [13,14], or the combined effect of increased precipitation and decreased potential evapotranspiration [10]. Climate, vegetation, and other factors, such as slope and altitude, also play a significant role in shaping the spatial distribution of WC [15]. However, comprehensive investigations into the spatial responses of WC services in the TRHR under the synergistic influence of environmental elements are lacking. Assessing the changes in WC services and understanding the driving factors behind temporal and spatial variations are crucial. Moreover, further exploration is necessary for quantitatively examining of the interplay among distinct influencing factors [16].
Currently, the main methods used to evaluate WC include the water storage capacity [5] and the water balance methods [17]. The water storage capacity method is primarily used to assess the scale of water storage capacity in sample plots, considering the contributions of vegetation, dead branches, and soil in intercepting precipitation [18]. By contrast, the assessment of large-scale WC is typically accomplished through remote sensing inversion and ecological model simulation based on water balance [19]. Remote sensing products provide information on WC at larger time and spatial scales but may not accurately capture local parameters such as soil water storage capacity [20]. The Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) model, based on the water balance, is widely used for long-term WC simulation due to its relatively straightforward data preparation and its ability to capture the physical mechanisms of soil–vegetation–atmosphere interactions in terrestrial ecosystems [21]. Previous research has highlighted that neglecting the influence of interannual changes in vegetation may result in inaccurate simulations of hydrological models [22,23]. However, most studies that employ the InVEST model to assess WC have overlooked the role of interannual variation in vegetation, potentially leading to inaccuracies in estimating actual evapotranspiration [24]. Therefore, adequately considering vegetation dynamics is crucial to accurately representing the WC response.
This study aims to simulate and analyze the temporal and spatial characteristics of WC in the TRHR over the past 30 years. Correlation analysis and geographical detectors are used to investigate the combined impact of climate, vegetation, and other driving factors on WC variations. A key aspect of this study is improving the InVEST water yield model by integrating long-term vegetation dynamics and ground-based meteorological observations, thus highlighting the impact of interannual vegetation variation. Specifically, considering its importance as a parameter controlling land–atmosphere interactions and water cycling, a regional correction of actual evapotranspiration on alpine grasslands was conducted [25]. The differences in the interaction between influencing factors on WC in different eco-geographical regions were also emphasized. Overall, this research deepens the understanding of land–atmosphere interactions in the TRHR, provides valuable insights for maintaining ecosystem stability in alpine areas and offers scientific support for ecosystem management and decision making processes.

2. Materials and Methods

2.1. Research Area

The TRHR (31°39′–36°12′ N, 89°45′–102°23′ E) is located in the northwestern part of Qinghai Province and is renowned as the birthplace of the Yellow River, Yangtze River and Lancang River (Figure 1a). This region has an average elevation of over 4400 m (Figure 1a) and is characterized by predominantly low temperatures. A west-to-east elevation trend is exhibited in this region, with altitudes ranging from 2300 to 6600 m. Slopes in the area range from 0 to 26 degrees, with steeper slopes (>10 degrees) concentrated in the southeastern and central-southern parts of the region (Figure 1e). The annual average temperature fluctuates between –8 °C and 17 °C, leading to widespread permafrost in the TRHR. Precipitation in the area is typically concentrated between June and September, with a gradient of decreasing rainfall from southeast to northwest, ranging from over 840 mm to less than 160 mm (Figure 1b). The region also experiences high potential evapotranspiration, with a multi-year average of 600 mm and a range between 400 mm and 1000 mm (Figure 1c). The ecosystem in the TRHR is relatively simple in structure and fragile in terms of its environmental conditions, making it highly vulnerable to the direct impacts of climate change. These unique climatic characteristics in the region have resulted in the formation of distinct eco-geographical regions (Figure 1f, where I and II represent sub-cold and temperate, and B, C, and D represent semi-humid, semi-arid, and arid zones, respectively) [26]. The dominant ecosystem type in the TRHR is alpine grassland, which covers 71.18% of the total area (Figure 1a). Thus, the region has extensive coverage of alpine shrub meadows, alpine grasslands, alpine meadow steppes, and alpine meadows (Figure 1f), collectively forming the alpine grassland. Additionally, the land use distribution in the TRHR includes farmland, woodland, water bodies, built-up land, and unused land, accounting for 0.66%, 4.26%, 5.63%, 0.1%, and 18.17% of the total land area in 2020, respectively. Unused land is located primarily in the western part of the study area, while forest land is in the central-southern and eastern regions. The leaf area index (LAI) in the TRHR exhibits a spatial gradient decreasing from southeast to northwest. The annual average LAI value ranges from 0 to 2 m2 m−2. The eastern region exhibits an average LAI value exceeding 1 m2 m−2, while the western region has an average LAI value below 0.2 m2 m−2 (Figure 1d).

2.2. Data Source and Preprocessing

Monthly observational data were gathered from 20 weather stations located in the TRHR. These data cover the period from 1991 to 2020 and include parameters such as temperature, precipitation, wind speed, and relative humidity. A thin plate spline function was employed to spatially interpolate these site-based observational data, utilizing elevation as a covariate, to obtain regional-scale data across the entire region. Potential evapotranspiration data were derived from the 1 km monthly potential evapotranspiration dataset in China provided by the National Tibetan Plateau Data Center [27]. Previous studies have demonstrated the effectiveness of this dataset in simulating the WC in the Tibetan Plateau [16].
The vegetation dataset used in the current study included the LAI and the homogenized vegetation index (NDVI). The LAI datasets were utilized to calculate the annual vegetation evapotranspiration coefficient (Kc) of the region. In contrast, the NDVI datasets were employed to assess the explanatory power of vegetation factors regarding the spatial distribution of WC. The Global Mapping (GLOBMAP) LAI version 3 dataset was acquired to obtain the required LAI data [28]. This dataset provided the necessary LAI data with a temporal resolution of 16 days before 2001 and 8 days after. The 1 km spatial resolution NDVI data were downloaded from Zenodo (https://zenodo.org/record/6295928, accessed on 1 September 2022). Monthly vegetation index data from 1991 to 2020 were generated using the maximum synthesis method.
Spatial data for the velocity coefficient were obtained based on empirical values associated with different land cover types to calculate water yield using the InVEST model [29]. Root depths for different vegetation types were determined by referring to relevant data from the Food and Agriculture Organization (FAO). Woodland areas were assigned a root depth of 3 m, grassland areas were assigned a root depth of 1 m, and non-vegetated areas were not assigned a root depth.
The annual total water yields (1994–2020) in the Qinghai Water Resources Bulletin were obtained to assess the accuracy of the model simulation. All the data used in this study are listed in Table 1.

2.3. Methods

2.3.1. Water Conservation Quantification

The WC for each grid was calculated using the water yield obtained from the InVEST model. This calculation considers various factors, including soil permeability, topography, and flow velocity coefficients, which were specific to different land-use types [29]:
W C i j = M i n 1 , 249 V e l o c i t y j · M i n 1 , 0.9 T I i 3 · M i n 1 , K s a t i 300 · Y i j
where W C i j (mm) is the WC for each land use type ( j ) in each grid ( i ). The influence of different land use types on surface runoff was considered by using the flow velocity coefficient ( V e l o c i t y j ). The topographic index ( T I i ), expressed as a dimensionless parameter, was used to assess the topographic characteristics of each grid ( i ). The Cosby model [30] was utilized to calculate the soil-saturated hydraulic conductivity ( K s a t i ) (cm/d) in centimeters per day. Additionally, the water yield ( Y i j ) was determined using the InVEST model and measured in millimeters (mm).
T I = l o g ( D a r e a S o i l d e p t h P s l o p e )
where D a r e a is the number of grids in the catchment area (dimensionless), S o i l d e p t h is the soil depth (mm), and P s l o p e is the percentage slope (%) calculated using the suite of Hydrology Tools (Arc Hydro) in ArcGIS 10.5.
K s a t = 60.96 10 ( 0.6 + 0.0126 s a n d 0.0064 c l a y )
where s a n d and c l a y represent the content of sand and clay in the soil (%), respectively.
Water yield in InVEST is defined as the amount of water that runs off the landscape. It calculates water yield by applying the water balance principle at the sub-watershed level. The conceptual diagram of the water balance model is illustrated in Figure 2.
The water yield model is based on the Budyko curve and annual average precipitation. The annual water yield for each pixel was calculated as follows:
Y x = 1 A E T x P x P x
where A E T x is the annual actual evapotranspiration of grid x and P x is the annual precipitation of grid x .
In the water balance formula, the vegetation evapotranspiration is calculated based on the Budyko hypothesis proposed by Zhang et al. [32]:
A E T x P x = 1 + K c l x · E T 0 x P x 1 + K c l x · E T 0 x P x w 1 / w
where E T 0 x is the potential evapotranspiration of grid x. K c l x represents the evapotranspiration coefficient of vegetation for land use type j and grid x. An estimation of Kc values for farmland and woodland was obtained using the relationship Kc = LAI/3 (Kc equal to 1 for LAI > 3). The Kc values for water bodies, built-up land, and unused land were set to 1.1, 0.3, and 0.5, respectively [31]. Kc for grassland was estimated using the method in Section 2.3.2. w x is an empirical parameter proposed by Donohue et al. [33]:
w x = Z A W C x P x + 1.25
where Z is the Zhang coefficient, which corresponds to the number of precipitation occurrences every year. When Z = 3.5, the error between the simulated water yield and actual runoff is the smallest (7.9%) [34]. A W C x represents the available water content of soil (mm), which is used to determine the total water stored and provided by soil for plant growth. AWC(x) is calculated as follows:
A W C x = M i n r o o t   r e s t r i c t i n g   l a y e r   d e p t h ,   r o o t   d e p t h × P A W C
where AWC(x) is estimated as the product of the plant available water capacity (PAWC) and the minimum of root-restricting layer depth and vegetation rooting depth. The calculation of AWC is performed using the following formula for PAWC:
P A W C = 54.509   0.132 s a n d % 0.003 ( s a n d % ) 2 0.055 s i l t %   0.006 ( s i l t % ) 2 0.738 c l a y % + 0.007 ( c l a y % ) 2   2.688 O M % + 0.501 ( O M % ) 2
Actual evapotranspiration (AET) for other land use types (water bodies, built-up land, unused land) is calculated directly with E T ( x ) :
A E T x = M i n ( K c l x · E T x , P x
where ET represents potential evapotranspiration and K c l x represents the evapotranspiration coefficient for the land use type without vegetation cover. The value of K c l x is taken in accordance with a fixed value (x) [35]: the value of the water body is 1.1, the value of built-up land is 0.3, and the value of unused land is 0.5.

2.3.2. Calibration Method of Evapotranspiration Coefficient in Grassland Based on Interannual Variation in Vegetation

Alpine grassland covers more than 70% of the total area of the TRHR, making it the dominant terrestrial ecosystem type. Therefore, the main objective of this study is to calibrate the Kc coefficient specifically for alpine grasslands. Capturing the dynamic evapotranspiration patterns of alpine grasslands within the InVEST model accurately is crucial. The growth process of alpine grassland can be categorized into different stages, which exhibit interannual variability influenced by climatic changes.
Monthly grid data of LAI were employed to estimate the Kc coefficient for each grid. Subsequently, the monthly average Kc coefficient was calculated for the entire study area by averaging the grid-based Kc coefficients. Finally, the 12 monthly regional average Kc coefficients for each year were combined to generate the annual value of Kc coefficients for grassland, following the approach recommended by the InVEST model. The specific procedure is as follows.
First, the monthly average Kc coefficients of grassland were computed based on the different growth stages of vegetation. Previous research and actual meteorological conditions identify three distinct stages for the growth season of alpine grasslands: the initial growth season (May), the middle growth season (June, July, and August), and the late growth season (September and October) [36]. Following the FAO-56 guidelines, the recommended values for Kc during the initial, middle, late, and non-growing stages of grassland are 0.4, 1.05, 0.85, and 0.4, respectively [35]. However, these coefficients can be adjusted to reflect the actual climatic conditions of the study area. The Kc coefficients for the middle and late growth stages were modified. Notably, the Kc coefficients for the early growth and non-growth stages remained unadjusted. The adjustment equation employed in this study is as follows:
K c = K c r e c o + [ 0.04 u 2 2 0.004 ( R H m i n 45 ) ] ( h / 3 ) 0.3
where R H m i n is the mean value for the daily minimum relative humidity (%), u 2 is the wind speed at 2 m height (m/s), and K r e c o is the recommended Kc value in FAO-56. h is the mean plant height (m), fitted by wind speed, relative humidity, and LAI [37]:
h = e x p [ 2 3 L A I 5.5 ]
Second, the annual vegetation evapotranspiration coefficients ( K c _ a n n u a l ) was determined using the calculation formula recommended by the InVEST model. Each month’s regional average Kc coefficients was calculated by averaging the grid-based monthly Kc coefficients across the entire study area. Then, each month’s regional average potential evapotranspiration was computed. Finally, the annual Kc coefficient was derived by using the following formula:
K c _ a n n u a l = m = 1 12 K o m × E T o m m = 1 12 E T o m
where K o m represents an average Kc coefficient of month m (1–12) and E T o m is the corresponding ETo.
The differences between the observed yearly total water yield (TWY) and the simulated TWY were evaluated to investigate the accuracy of the simulations by the improved InVEST model using the regionally calibrated Kc coefficient, including the coefficient of determination (R2) and the root mean square error (RMSE). The differences between the observed and simulated TWY obtained through three different methods of calculating Kc are further compared to illustrate the effectiveness of the improved model (Figure 3).
The approach proposed in Section 2.3.2 was utilized in the first method to estimate the Kc coefficients (Figure 3a). In the second method, the Kc coefficient was determined using the relationship Kc = LAI/3 (where Kc equals 1 when LAI > 3) (Figure 3b). This method is mentioned in the InVEST model manual. The calculation process involves synthesizing the monthly values of LAI using the maximum value synthesis method, followed by averaging the values for each month across each year to derive the annual average LAI. The Kc value for each grid is computed based on this relationship Kc = LAI/3, and the overall Kc value required by the model operation is obtained by averaging across the entire region. The third method uses a fixed Kc coefficient of 0.65, corresponding to the recommended Kc value for grassland in the InVEST model (Figure 3c). The observed TWY data are obtained from the Qinghai Water Resources Bulletin, covering 1994 to 2020. Comparatively, the simulation results of TWY based on the improved model display better performance, with an increase in the R2 ranging from 0.51% to 14.6% and a reduction in RMSE ranging from 0.478 × 1010 m3 to 2.193 × 1010 m3. The improved model enhances the performance of the InVEST water yield model by incorporating information on interannual vegetation changes.

2.3.3. Geographical Detector

The geographical detector method is a powerful tool for analyzing the driving forces by detecting the heterogeneity of the spatial stratification of elements [38]. Factor detectors and interaction detectors were utilized to examine the driving factors of WC. This study focused on two important aspects: identifying the factors that significantly influence WC and investigating the interactions between these factors.
The factor detector is employed to detect the spatial differentiation of dependent variables and assess the explanatory power of independent variables for WC. The q-value quantifies the strength of this explanatory power. The specific calculation formula is as follows:
q = 1 h = 1 L N h σ h 2 N σ 2 = 1 S S W S S T
S S W = h = 1 L N h σ h 2   S S T = N σ 2
where h is the stratification of the dependent or independent variable, with values ranging from 1 to L. N h   a n d   N are the number of cells within stratum h and the whole region, respectively. σ h 2 and σ 2 are the variances of the values of the dependent variable in the stratum and the whole region, respectively. S S W is the sum of the variances within the stratum, and S S T is the total variance of the whole region. q indicates the explanatory power of independent variables for dependent variables, with a value range between 0 and 1. A higher value of q indicates a stronger explanatory power of independent variables for dependent variables. The main influencing factors of WC can be determined by examining the value of q . The q value can be tested for significance using the geographic detector method because it follows a non-central distribution after a simple transformation.
This study assesses whether the two factors have an interaction. It determines the strength and characteristics of the interaction by comparing the values obtained from individual factors and the interaction values obtained from the geographic detector. For example, if q ( X 1 X 2 ) < M i n ( q ( X 1 ) , q ( X 2 ) ) , it indicates that the interaction is characterized by nonlinear weakening. If M i n ( q ( X 1 ) , q ( X 2 ) ) < q ( X 1 X 2 ) < M a x ( q ( X 1 ) , q ( X 2 ) ) , the interaction shows the characteristics of one-factor nonlinearity weakening. If q ( X 1 X 2 ) > M a x ( q ( X 1 ) , q ( X 2 ) ) , the interaction displays the characteristics of two-factor enhancement. If q ( X 1 X 2 ) = q ( X 1 ) + q ( X 2 ) , it signifies that the interaction is characterized by two-factor independence. Finally, if q ( X 1 X 2 ) > q ( X 1 ) + q ( X 2 ) , it indicates the interaction is characterized by nonlinear enhancement.
The objective of this study is to examine the factors that influence the spatial heterogeneity of WC. These factors include climatic factors such as annual average precipitation and annual average potential evapotranspiration, topographic factors such as slope and digital elevation model, and vegetation factors such as annual average growing season LAI. Continuous variables need to be converted into classified data to meet the requirements of the geographic detectors. Therefore, the natural breakpoint method in ArcGIS 10.5 was utilized to divide altitude, slope, annual average precipitation, annual average potential evapotranspiration, and LAI into nine categories. The analysis process of the geographic detector using NDVI as the vegetation factor was repeated to ensure the reliability of the results.

2.3.4. Trend Changes and Correlation Analysis

Simple linear regression analysis was initially employed to examine the interannual variation trends within the TRHR at regional and grid scales. Subsequently, Pearson correlation analysis was used to investigate the impact of individual environmental factors on the interannual variability in WC.

2.3.5. Total Research Approach

The specific research schemes are outlined as follows: (1) first, the monthly Kc coefficients are calculated using the monthly LAI grid data to correct for the grassland’s middle and late growth stages. The method recommended in the InVEST model manual is used to synthesize the monthly Kc coefficients into an annual value. Other parameter values necessary for model execution can be found in Section 2.3.1 and Section 2.3.2. (2) Next, the InVEST model is utilized to simulate the water yield of the TRHR from 1991 to 2020. The TWYs obtained from the model output are compared with the recorded TWYs in the water resources bulletin to verify the accuracy of the simulation results. Once the verification is complete, we estimate the WC based on the grid-scale water yield results. (3) The influence of vegetation and climate on the temporal changes in WC is analyzed using correlation analysis. Additionally, the geographic detector method is employed to explore the effects of climate, vegetation, and topography on the spatial differentiation of WC and the interactions among these factors (Figure 4).

3. Results

3.1. Characteristics of Interannual Variation in Water Conservation

After regional correction, the multi-year average Kc coefficient for the grassland is 0.72, which is higher than the recommended value of 0.65 suggested by the InVEST model (Figure 5). Overall, the Kc coefficient showed a slight downward trend over the years, which was not very pronounced. The coefficient of variation of Kc from 1990 to 2000 was 0.018, indicating significant fluctuation of the Kc coefficient of grassland during that period. By contrast, the overall coefficient of variation from 1990 to 2020 decreased to 0.013, indicating decreased fluctuation in the 2000s and 2010s. Among the years analyzed, the evapotranspiration coefficient of grassland reached its maximum value of 0.75 in 1997 and its minimum value of 0.71 in 1993.
WC in the TRHR exhibited a significant increase (p < 0.05) from 1991 to 2020, with an annual growth rate of 1.4 mm (Figure 6). A series of ecological projects have been implemented, and the ecological environment in the TRHR has improved since 2000. Accordingly, WC has changed from a weak increase to a strong increase since 2000. Specifically, WC demonstrated a nonsignificant upward trend (1.16 mm/yr) from 1991 to 1999. However, a significant increase was observed in WC, with a growth rate of 1.74 mm/yr (p < 0.05) from 2000 to 2020. Over the past 30 years, the average annual WC in the TRHR was 77.82 mm. During the 1990s, the average annual WC in this region was 65 mm. The annual average WC has increased to 83.31 mm since 2000, representing an overall growth of 28.17%. The change trend of WC is similar to that of precipitation, with a change rate of 3.71 mm/yr before 2000 and 3.8 mm/yr from 2000 to 2020. The lowest value for WC occurred in 2015 (46.31 mm), while the highest was observed in 2019 (114.76 mm). The average annual total WC over the years was 2.79 × 1010 m3. The overall trend of total WC aligns with the trend of WC over time.

3.2. Spatial and Temporal Distributions of Water Conservation

The spatial distribution of WC in the TRHR shows significant heterogeneity (Figure 7a). The areas with high WC were found mainly in the southeast of the Yellow River and the source area of the Lancang River, specifically in Banma County, Nangqian County, and Jiuzhi County. The average value of WC over the years was 77.82 mm. Among the sub-regions, the mean value of WC was highest in the HIB1 region (125.32 mm), followed by the HIIC2 (112.97 mm) and the HIC1 regions (71.97 mm). The HIB1 region is characterized by lower altitude and a more humid climate favorable conditions for WC. Moreover, this area boasts relatively high vegetation coverage, pivotal in regulating runoff and increasing water yield, making it a promising region for WC.
WC exhibited an increasing trend in 96.12% of the region during the study period, and among these trends, 78.17% were statistically significant (p < 0.05) (Figure 7b). These findings highlight substantial WC improvements across most TRHR over the past three decades. Notably, the southeastern region of the TRHR exhibited an even more pronounced trend, demonstrating a rate exceeding 2 mm/yr. However, approximately 3.87% of the region still faces challenges maintaining WC efforts. Moreover, approximately 1.72% of the total area of the TRHR experienced a significant decrease in WC throughout the study period, with a rate of –2.77 mm/yr (p < 0.05).

3.3. Spatial and Temporal Driving Factors of Water Conservation

3.3.1. Drivers of Interannual Change

The TRHR experienced an increase in average annual precipitation during the period from 1991 to 2020, with a rate of 3.6 mm/yr (p < 0.5) (Figure 8a). Before 2000, the increase in precipitation was not statistically significant, with an annual increase of only 3.71 mm/yr. However, the rate of precipitation increase has noticeably accelerated since 2000, with an annual increase of 3.80 mm/yr. By contrast, potential evapotranspiration demonstrated no significant change throughout the study period, with a rate of 0.48 mm/yr (Figure 8b). The rate of increase in potential evapotranspiration slowed after 2000. Additionally, the average annual LAI from 1991 to 2020 exhibited an upward trend, indicating improved vegetation growth (Figure 8c). The TRHR generally had a low baseline LAI, but a fluctuating upward trend was observed from 1991 to 2020, with a growth rate of 0.0008 m2 m−2/yr (p < 0.05). In particular, the annual fluctuation in LAI was more pronounced during 2000–2020 than in the 1991–2000 period.
Table 2 shows a significant positive correlation between precipitation and WC. The increase in precipitation in the study area over the past 30 years has played an important role in improving WC. However, no significant correlation between potential evapotranspiration and WC was observed throughout the study. Before 2000, potential evapotranspiration and WC had a positive correlation, but this relationship was not statistically significant. However, this correlation changed to a significant negative correlation after 2000. The increase in potential evapotranspiration indicates that meteorological factors such as temperature, radiation, and wind speed have a stronger ability to demand atmospheric water, ultimately leading to an increase in actual evapotranspiration and a decrease in WC. The increase in LAI also promotes greater WC. After 2000, the influence of LAI on WC increased.

3.3.2. Drivers of Spatial Heterogeneity

The geographical detector method was employed to identify the primary factors influencing the spatial distribution of the annual average WC in the region, given the significant spatial variability of WC within the TRHR. Analysis of the q values of each factor indicated that annual average precipitation had the highest explanatory power for WC (0.73), followed by LAI (0.582), annual average potential evapotranspiration (0.207), digital elevation model (DEM) (0.221), and slope (0.109), as determined. These factors were statistically significant (p < 0.01) (Table 3). The annual average precipitation accounted 73% of the spatial distribution of WC, which was the primary driving factor in the region. LAI demonstrated moderate explanatory power for WC, while potential evapotranspiration and altitude had weak influences. The slope had the weakest influence on the spatial heterogeneity of WC, contributing only 10.9%. A multi-year average NDVI from 1991 to 2020 was used as a vegetation factor, yielding consistent findings with those obtained using LAI to ensure the reliability of the results.
An interaction detector was used to analyze the interaction between two factors in the TRHR and various ecological regions to explore further the impact of different factors on the spatial variation in WC (Figure 9). The NDVI data were applied to replicate the entire interaction detection process and validate the reliability of the interaction detection findings. The entire interaction detection process was also replicated using NDVI data to validate the reliability of the findings. The interaction between two factors provided a better explanation for the spatial differentiation of WC compared to a single factor. In particular, the interaction between precipitation and other factors played a key role in the spatial differentiation pattern of WC. The interaction between precipitation and potential evapotranspiration had the strongest influence on WC in the north-central and western regions (HIC1) (explaining at least 63.1% of the spatial distribution of WC), where the background value of LAI was relatively low. Precipitation and vegetation together accounted for 72.8% and 60% of the distribution of WC in regions HIIC2 and HIB1, respectively.

4. Discussion

4.1. Interannual Variation in Water Conservation

Vegetation dynamics affect the regional WC through their influence on rainfall-runoff capacity [39] and evapotranspiration [40]. Thus, considering the impact of vegetation dynamics is crucial in quantifying WC services. Some hydrological models, including the Soil and Water Assessment Tool (SWAT), Water and Energy transfer Processes in Large river basins (WEP-L), and Water and Energy Budget Distributed Hydrological Model (WEB-DHM), also consider the effects of vegetation growth on precipitation distribution, soil moisture, and groundwater recharge [41,42,43]. This study focused specifically on the influence of vegetation on WC and improved the InVEST model to reflect interannual variation in vegetation and enhance the accuracy of the simulation.
The grassland growth and actual evapotranspiration in the Tibetan Plateau exhibit evident stage and interannual variations [44]. An empirical coefficient called the Kc coefficient is commonly used to determine vegetation water requirements. This coefficient is typically calculated by dividing actual evapotranspiration (AET) by reference evapotranspiration (ETo). Previous studies have employed this ratio to estimate the annual Kc coefficient. However, obtaining AET data in advance is necessary for the corresponding research, which limits the generalizability of this approach [45]. The Kc coefficient for a specific region can be determined by establishing an empirical relationship between vegetation indexes, such as LAI and NDVI, and the corresponding Kc coefficients. However, collecting and obtaining the observed Kc coefficients beforehand to establish this relationship is crucial [46]. This study first categorized grassland Kc coefficients into three growth stages to better understand the variations in grassland evapotranspiration in this research. Then, the Kc coefficients recommended by FAO-56 for each specific stage of grassland growth were utilized. The Kc coefficients for the middle and late growth stages were also refined by incorporating local climate data. These adjusted coefficients were combined to derive annual-scale grassland Kc coefficients. The methodology developed in this study enables rapid and accurate estimation of Kc coefficients across a wide range of temporal and spatial scales. However, when estimating the Kc coefficient, the absence of extensive plant height data requires estimating plant height using LAI and an empirical formula. However, this approach is specific to alpine grasslands because of the regional nature of the empirical formula employed and may not be apply to other regions. If simulated or observed plant height data become available for a broader range of times and locations, the method proposed in the current research could potentially be applied to grassland areas on a global scale.

4.2. Factors Affecting Temporal and Spatial Differentiation of Water Conservation

In the past 30 years, the WC in the TRHR has increased significantly, which is consistent with previous studies [8,13,47]. WC is a complex process influenced by various factors. Climatic factors, including precipitation, temperature, and humidity directly impact on the hydrological cycle and the ability to conserve water [48,49]. The relationship between interannual variations in WC and climatic factors was examined in the TRHR. The results indicated a significant positive correlation between precipitation and WC over the past 30 years, in line with previous research [3,13,16]. Precipitation emerged as the primary factor influencing the spatial variation in WC in the TRHR, followed by vegetation (measured by NDVI or LAI) consistent with previous studies in the region [14,16]. The evapotranspiration capacity of the terrestrial ecosystem, indicated by potential evapotranspiration, has a great effect on WC in the TRHR, along with vegetation. The TRHR displays significant spatial variations in WC, demonstrating high WC values found in areas with abundant precipitation and well-developed vegetation. Topography regulates WC by impacting vegetation structure, soil properties, and water and heat conditions [50]. Nevertheless, slope has a weaker effect on WC than other influencing factors.
The implementation of ecological projects in the TRHR also contributes to enhancing WC. Grassland degradation has been observed in this region since the 1980s, reducing effectiveness in intercepting runoff and storing water [51]. Since 2000, the government has initiated several ecological projects to address this issue, including the grassland ecological reward mechanism, the grazing withdrawal project, and the first and second stages of ecological conservation and restoration projects [52,53]. These projects transformed approximately 1.98 × 104 km2 of unused land into grassland in the TRHR by 2020. Furthermore, a remarkable increase in vegetation coverage has been observed in the entire region over the past 30 years, stabilizing surface temperature and enhancing soil moisture conditions [54]. Implementing ecological projects has effectively improved WC in the TRHR, strengthening its ecological security barrier function [55].
A single factor alone is often insufficient to comprehend the entire ecosystem’s complexity fully. The complexity of geographical processes typically involves the interplay of multiple factors. While the interaction between precipitation and vegetation enhances their explanatory power for WC, the degree of interaction may vary across eco-geographical regions. In the TRHR, HIC1 experiences relatively low precipitation and has a low background value of vegetation. Therefore, the interaction between precipitation and potential evapotranspiration effectively explains the distribution of WC in this area. By contrast, HIB1 and HIIC2 are relatively moist and have abundant vegetation. Vegetation plays a crucial role in regulating the distribution of precipitation and ensuring WC through transpiration and interception [56]. Thus, the interaction between vegetation and precipitation in these areas has a high explanatory capacity for WC.

4.3. Adaptation to Future Climate Change

The Tibetan Plateau is predicted to experience increased precipitation in the future, which will positively impact WC in most areas of the TRHR, such as the Yangtze River and Lancang River sources [3,57]. However, certain areas, particularly some counties at the source of the Yellow River, face a decline in WC because of reduced precipitation and increased potential evapotranspiration resulting from a weakening summer monsoon. This increase in rainfall is accompanied by higher rainfall erosivity, leading to soil erosion and uncertainties in WC [58]. The Tibetan Plateau is recognized as one of the most sensitive regions to global climate change, with its warming rate twice exceeding the global average [59]. The TRHR, located in the hinterland of the Tibetan Plateau, has observed even more significant warming changes. The increasing frequency of extreme heat events in the future will have a profound impact on the terrestrial ecosystem balance in this region [60]. Rising air temperatures lead to a decrease in soil moisture [61] and an increase in the saturated vapor pressure deficit [62], thereby increasing the vulnerability of WC services in this region. Climate change is anticipated to intensify vegetation expansion (greening) in the TRHR [63]. The development of grasslands is closely linked to WC services [64]. This large-scale vegetation greening will cause changes in evapotranspiration and precipitation patterns [65]. Permafrost in the TRHR has also been declining since 1980, while the thickness of the active layer has been increasing [51,66]. Hence, future research should monitor the impact intensity of vegetation greening and permafrost degradation on WC. Changes in WC also affect the terrestrial ecosystem in alpine regions [67]. Therefore, understanding WC services and their interactions with various environmental factors is crucial to developing effective measures to adapt to climate change.

4.4. Limitations and Future Work

WC requires a comprehensive understanding of the ecohydrological processes within a basin, including factors such as precipitation interception and evaporation. It also involves exploring the interaction between multiple hydrological processes and the ecosystem, such as runoff timing, water quality, flood protection, and temperature regulation. Runoff timing can indicate the capacity for WC. If the surface and groundwater systems have a strong water storage capacity, runoff from rainfall or snowmelt will slowly enter rivers or lakes, delaying runoff formation. In such cases, the runoff time will be relatively late [68]. WC processes reduce soil erosion and prevent pollutants from entering water bodies. For example, vegetation slows down the flow of rainfall, allowing more time for water to permeate the soil, which reduces surface runoff and helps decrease the risk of soil erosion and the introduction of pollutants into water bodies [69]. Additionally, a basin with a good WC capacity mitigates the likelihood of natural disasters such as floods and droughts by intercepting rainfall and storing soil water [70,71]. These aspects of WC function will be considered to further enhance the understanding of WC in future research.
Research methods employed in this study have several limitations. First, the WC simulation does not differentiate between the Kc coefficients at grid scales. Instead, a regional average Kc coefficient is utilized, which may introduce certain uncertainties to the results. Second, the InVEST model does not account for the influence of runoff in the calculation of regional WC, which is another limitation [31,72]. the model also fails to consider the impact of freezing and thawing on WC [14]. Therefore, future research should focus on developing more precise methods for simulating vegetation evapotranspiration coefficients at the grid scale and exploring ways to optimize WC models in plateau regions.
The environment in the TRHR is highly complex, and this study has not delved deeply into the various factors contributing to changes in WC. Human activities such as overgrazing, expansion of agricultural land, road construction projects, and ecological project initiatives alter the surface, directly or indirectly affecting the WC services of the region [73,74,75]. Future studies should focus on quantitatively analyzing the positive and negative impacts of human activities on WC to gain a more comprehensive understanding. Differentiating between the effects of human activities and natural environmental factors on WC is crucial.

5. Conclusions

In this study, we improved the InVEST model by integrating the long-term series of LAI data to capture the interannual variation in vegetation. The temporal and spatial distribution characteristics of WC in the TRHR were assessed quantitatively from 1991 to 2020. Trend and correlation analyses were used to examine the temporal variation trend and driving factors of WC. We also employed the geographical detector method to explore the influencing factors and their interactions in the spatial distribution of WC. The main findings are as follows:
(1)
WC in the TRHR exhibits a spatial pattern with high values in the south and low values in the north. WC significantly increased (1.4 mm/yr, p < 0.05) from 1991 to 2020. Compared with 1991–1999, the annual average value of WC in 2000–2020 increased by 28.17%. Over the past 30 years, WC improved in 78.17% of the regions (p < 0.5).
(2)
The increase in precipitation over the past three decades has had a significant positive impact on WC (R = 0.97, p < 0.01). However, the growth in potential evapotranspiration has had a significant inhibitory effect on WC since 2000 (R = –0.5, p < 0.05).
(3)
The spatial variation in WC is influenced primarily by precipitation, followed by vegetation and potential evapotranspiration. The interaction among these factors has a stronger explanatory power for WC than individual factors alone. The interaction between precipitation and other influencing factors demonstrates the greatest explanatory power. The combined influence of precipitation and vegetation accounts for approximately 79.1% of the WC distribution across the study area. However, the dominant interaction factors for WC vary in different eco-geographical regions. In the north-central and western regions (HIC1) with low vegetation, the interaction between annual precipitation and potential evapotranspiration explains 65% of the variation in WC, making it the dominant interaction factor in the region. In the eastern and central-southern areas (HIB1 and HIIC2), the interaction between annual precipitation and vegetation exhibits the strongest influence.
The method proposed provides valuable insights into simulating WC. Incorporating processes such as permafrost degradation is crucial for the accurate assessment of WC in plateau regions in future work.

Author Contributions

Conceptualization, data curation, formal analysis, methodology, writing—original draft, Y.P.; methodology, conceptualization, funding acquisition, supervision, writing—review and editing, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Second Tibetan Plateau Scientific Expedition Program (2019QZKK0403).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yao, Y.; Fu, B.J.; Liu, Y.X.; Li, Y.; Wang, S.; Zhan, T.Y.; Wang, Y.J.; Gao, D.X. Evaluation of ecosystem resilience to drought based on drought intensity and recovery time. Agric. For. Meteorol. 2022, 314, 108809. [Google Scholar] [CrossRef]
  2. Yao, T.D.; Wu, G.J.; Xu, B.Q.; Wang, W.C.; Gao, J.; An, B.S. Asian Water Tower Change and Its Impacts. Bull. Chin. Acad. Sci. 2019, 34, 1203–1209. (In Chinese) [Google Scholar]
  3. Xue, J.; Li, Z.X.; Feng, Q.; Gui, J.; Zhang, B.J. Spatiotemporal variations of water conservation and its influencing factors in ecological barrier region, Qinghai-Tibet Plateau. J. Hydrol. Reg. Stud. 2022, 42, 101164. [Google Scholar] [CrossRef]
  4. Xu, W.H.; Xiao, Y.; Zhang, J.J.; Yang, W.; Zhang, L.; Hull, V.; Wang, Z.; Zheng, H.; Liu, J.G.; Polasky, S. Strengthening protected areas for biodiversity and ecosystem services in China. Proc. Natl. Acad. Sci. USA 2017, 114, 1601–1606. [Google Scholar] [CrossRef] [PubMed]
  5. Hu, W.; Li, G.; Gao, Z.; Jia, G.; Wang, Z.; Li, Y. Assessment of the impact of the Poplar Ecological Retreat Project on water conservation in the Dongting Lake wetland region using the InVEST model. Sci. Total Environ. 2020, 733, 139423. [Google Scholar] [CrossRef] [PubMed]
  6. Biskaborn, B.K.; Smith, S.L.; Noetzli, J.; Matthes, H.; Vieira, G.; Streletskiy, D.A.; Schoeneich, P.; Romanovsky, V.E.; Lewkowicz, A.G.; Abramov, A.; et al. Permafrost is warming at a global scale. Nat. Commun. 2019, 10, 264. [Google Scholar] [CrossRef] [PubMed]
  7. Ma, N.; Zhang, Y.Q. Increasing Tibetan Plateau terrestrial evapotranspiration primarily driven by precipitation. Agric. For. Meteorol. 2022, 317, 108887. [Google Scholar] [CrossRef]
  8. Xiao, Y.; Xiong, Q.L.; Liang, P.H.; Xiao, Q. Potential risk to water resources under eco-restoration policy and global change in the Tibetan Plateau. Environ. Res. Lett. 2021, 16, 094004. [Google Scholar] [CrossRef]
  9. Shao, Q.Q.; Cao, W.; Fan, J.W.; Huang, L.; Xu, X.L. Effects of an ecological conservation and restoration project in the Three-River Source Region, China. J. Geogr. Sci. 2017, 27, 183–204. [Google Scholar] [CrossRef]
  10. Pan, T.; Wu, S.H.; Liu, Y. Relative Contributions of Land Use and Climate Change to Water Supply Variations over Yellow River Source Area in Tibetan Plateau during the Past Three Decades. PLoS ONE 2015, 10, e0123793. [Google Scholar] [CrossRef] [PubMed]
  11. Šatalová, B.; Kenderessy, P. Assessment of water retention function as tool to improve integrated watershed management (case study of Poprad river basin, Slovakia). Sci. Total Environ. 2017, 599, 1082–1089. [Google Scholar] [CrossRef]
  12. Liu, H.; Liu, S.L.; Wang, F.F.; Liu, Y.X.; Han, Z.L.; Wang, Q.B.; Yu, L.; Dong, Y.H.; Sun, J. Multilevel driving factors affecting ecosystem services and biodiversity dynamics on the Qinghai-Tibet Plateau. J. Clean. Prod. 2023, 396, 136448. [Google Scholar] [CrossRef]
  13. Liu, G.B.; Shao, Q.Q.; Fan, J.W.; Ning, J.; Huang, H.B.; Liu, S.C.; Zhang, X.Y.; Niu, L.; Liu, J.Y. Spatio-Temporal Changes, Trade-Offs and Synergies of Major Ecosystem Services in the Three-River Headwaters Region from 2000 to 2019. Remote Sens. 2022, 14, 5349. [Google Scholar] [CrossRef]
  14. Wang, Y.F.; Ye, A.Z.; Peng, D.Z.; Miao, C.Y.; Di, Z.H.; Gong, W. Spatiotemporal variations in water conservation function of the Tibetan Plateau under climate change based on InVEST model. J. Hydrol. Reg. Stud. 2022, 41, 101064. [Google Scholar] [CrossRef]
  15. Kurzweil, J.R.; Metlen, K.; Abdi, R.; Strahan, R.; Hogue, T.S. Surface water runoff response to forest management: Low-intensity forest restoration does not increase surface water yields. For. Ecol. Manag. 2021, 496, 119387. [Google Scholar] [CrossRef]
  16. Wen, X.; Shao, H.Y.; Wang, Y.; Lv, L.F.; Xian, W.; Shao, Q.F.; Shu, Y.; Yin, Z.Q.; Liu, S.H.; Qi, J.G. Assessment of the Spatiotemporal Impact of Water Conservation on the Qinghai–Tibet Plateau. Remote Sens. 2023, 15, 3175. [Google Scholar] [CrossRef]
  17. Yang, M.Y.; Dong, S.K.; Dong, Q.M.; Xu, Y.D.; Zhi, Y.L.; Liu, W.T.; Zhao, X.Y. Trade-offs in ecological, productivity and livelihood dimensions inform sustainable grassland management: Case study from the Qinghai-Tibetan Plateau. Agric. Ecosyst. Environ. 2021, 313, 107377. [Google Scholar] [CrossRef]
  18. Wu, H.F.; Hu, B.A.; Yan, J.X.; Cheng, X.Q.; Yi, P.H.; Kang, F.F.; Han, H.R. Mixed plantation regulates forest floor water retention and temperature sensitivity in restored ecosystems on the Loess Plateau, China. Catena 2023, 222, 106838. [Google Scholar] [CrossRef]
  19. Wu, X.; Shi, W.J.; Guo, B.; Tao, F.L. Large spatial variations in the distributions of and factors affecting forest water retention capacity in China. Ecol. Indic. 2020, 113, 106152. [Google Scholar] [CrossRef]
  20. de Figueiredo, T.; Royer, A.C.; Fonseca, F.; de Araujo Schutz, F.C.; Hernandez, Z. Regression Models for Soil Water Storage Estimation Using the ESA CCI Satellite Soil Moisture Product: A Case Study in Northeast Portugal. Water 2021, 13, 37. [Google Scholar] [CrossRef]
  21. Redhead, J.W.; Stratford, C.; Sharps, K.; Jones, L.; Ziv, G.; Clarke, D.; Oliver, T.H.; Bullock, J.M. Empirical validation of the InVEST water yield ecosystem service model at a national scale. Sci. Total Environ. 2016, 569, 1418–1426. [Google Scholar] [CrossRef] [PubMed]
  22. Yu, S.Q.; Qin, H.P.; Deng, W. Modeling the effects of vegetation dynamics on the hydrological performance of a bioretention system. J. Hydrol. 2023, 620, 129473. [Google Scholar] [CrossRef]
  23. Chen, S.Z.; Fu, Y.H.; Wu, Z.F.; Hao, F.H.; Hao, Z.C.; Guo, Y.H.; Geng, X.J.; Li, X.Y.; Zhang, X.; Tang, J. Informing the SWAT model with remote sensing detected vegetation phenology for improved modeling of ecohydrological processes. J. Hydrol. 2023, 616, 128817. [Google Scholar] [CrossRef]
  24. Trautmann, T.; Koirala, S.; Carvalhais, N.; Güntner, A.; Jung, M. The importance of vegetation to understand terrestrial water storage variations. Hydrol. Earth Syst. Sci. 2022, 26, 1089–1109. [Google Scholar] [CrossRef]
  25. Zhang, K.; Kimball, J.S.; Running, S.W. A review of remote sensing based actual evapotranspiration estimation. Wiley Interdiscip. Rev. Water 2016, 3, 834–853. [Google Scholar] [CrossRef]
  26. Zheng, D.; Fao, R. A Study on the Eco-Geographic Regional System of China. In Proceedings of the FAO FRA2000 Global Ecological Zoning Workshop, Cambridge, UK, 28–30 January 1999. [Google Scholar]
  27. Peng, S.Z.; Ding, Y.X.; Liu, W.Z.; Li, Z. 1 km monthly temperature and precipitation dataset for China from 1901 to 2017. Earth Syst. Sci. Data 2019, 11, 1931–1946. [Google Scholar] [CrossRef]
  28. Liu, Y.; Liu, R.G.; Chen, J.M. Retrospective retrieval of long-term consistent global leaf area index (1981–2011) from combined AVHRR and MODIS data. J. Geophys. Res. Biogeosci. 2012, 117, 1–14. [Google Scholar] [CrossRef]
  29. Wang, Y.C.; Zhao, J.; Fu, J.W.; Wei, W. Effects of the Grain for Green Program on the water ecosystem services in an arid area of China-Using the Shiyang River Basin as an example. Ecol. Indic. 2019, 104, 659–668. [Google Scholar] [CrossRef]
  30. Cosby, B.J.; Hornberger, G.M.; Clapp, R.B.; Ginn, T.R. A Statistical Exploration of the Relationships of Soil Moisture Characteristics to the Physical Properties of Soils. Water Resour. Res. 1984, 20, 682–690. [Google Scholar] [CrossRef]
  31. Sharp, R.; Tallis, H.T.; Ricketts, T.; Guerry, A.D.; Wood, S.A.; Chaplin-Kramer, R.; Nelson, E.; Ennaanay, D.; Wolny, S.; Olwero, N.; et al. Invest Version 3.2. 0 User’s Guide; The Natural Capital Project: Stanford, CA, USA, 2015. [Google Scholar]
  32. Zhang, L.; Hickel, K.; Dawes, W.R.; Chiew, F.H.S.; Western, A.W.; Briggs, P.R. A rational function approach for estimating mean annual evapotranspiration. Water Resour. Res. 2004, 40, W02502. [Google Scholar] [CrossRef]
  33. Donohue, R.J.; Roderick, M.L.; McVicar, T.R. Roots, storms and soil pores: Incorporating key ecohydrological processes into Budyko’s hydrological model. J. Hydrol. 2012, 436, 35–50. [Google Scholar] [CrossRef]
  34. Jiang, C.; Li, D.Q.; Wang, D.W.; Zhang, L.B. Quantification and assessment of changes in ecosystem service in the Three-River Headwaters Region, China as a result of climate variability and land cover change. Ecol. Indic. 2016, 66, 199–211. [Google Scholar] [CrossRef]
  35. Allen, R.; Pereira, L.; Raes, D.; Smith, M. Crop Evapotranspiration Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper, 56; FAO: Rome, Italy, 1998; pp. 1–289. [Google Scholar]
  36. Li, P.; Peng, C.H.; Wang, M.; Luo, Y.P.; Li, M.P.; Zhang, K.R.; Zhang, D.L.; Zhu, Q.A. Dynamics of vegetation autumn phenology and its response to multiple environmental factors from 1982 to 2012 on Qinghai-Tibetan Plateau in China. Sci. Total Environ. 2018, 637, 855–864. [Google Scholar] [CrossRef] [PubMed]
  37. Zhao, J.; Liu, C.Y.; Pan, J.H.; Liu, Y.Y.; Yang, D.H. Spatial and Temporal Distributions of Evapotranspiration in Gannan Grassland Based on MODIS Data. Resour Sci. 2011, 33, 341–346. (In Chinese) [Google Scholar]
  38. Wang, J.F.; Li, X.H.; Christakos, G.; Liao, Y.L.; Zhang, T.Z.; Gu, X.; Zheng, X.Y. Geographical Detectors-Based Health Risk Assessment and its Application in the Neural Tube Defects Study of the Heshun Region, China. Int. J. Geogr. Inf. Sci. 2010, 24, 107–127. [Google Scholar] [CrossRef]
  39. Sun, M.L.; Hu, J.; Chen, X.L.; Lü, Y.H.; Yang, L.X. Comparison of Five Models for Estimating the Water Retention Service of a Typical Alpine Wetland Region in the Qinghai–Tibetan Plateau. Remote Sens. 2022, 14, 6306. [Google Scholar] [CrossRef]
  40. Lian, J.J.; Li, D.F.; Huang, M.B.; Chen, H.S. Evaluation of remote sensing-based evapotranspiration estimates using a water transfer numerical simulation under different vegetation conditions in an arid area. Hydrol. Process. 2018, 32, 1801–1813. [Google Scholar] [CrossRef]
  41. Luan, J.K.; Miao, P.; Tian, X.Q.; Li, X.J.; Ma, N.; Faiz, M.A.; Xu, Z.W.; Zhang, Y.Q. Estimating hydrological consequences of vegetation greening. J. Hydrol. 2022, 611, 128018. [Google Scholar] [CrossRef]
  42. Peng, H.; Jia, Y.W.; Niu, C.W.; Gong, J.G.; Hao, C.F.; Gou, S. Eco-hydrological simulation of soil and water conservation in the Jinghe River Basin in the Loess Plateau, China. J. Hydro Environ. Res. 2015, 9, 452–464. [Google Scholar] [CrossRef]
  43. Wang, Y.W.; Wang, L.; Zhou, J.; Chai, C.H.; Hu, Z.D.; Zhao, L.; Wang, S.F.; Fan, M.T. Impacts of frozen ground degradation and vegetation greening on upper Brahmaputra runoff during 1981–2019. Int. J. Clim. 2023, 43, 3768–3781. [Google Scholar] [CrossRef]
  44. Shen, M.G.; Wang, S.P.; Jiang, N.; Sun, J.P.; Cao, R.Y.; Ling, X.F.; Fang, B.; Zhang, L.; Zhang, L.H.; Xu, X.Y. Plant phenology changes and drivers on the Qinghai–Tibetan Plateau. Nat. Rev. Earth Environ. 2022, 3, 633–651. [Google Scholar] [CrossRef]
  45. Dennedy-Frank, P.J.; Muenich, R.L.; Chaubey, I.; Ziv, G. Comparing two tools for ecosystem service assessments regarding water resources decisions. J. Environ. Manag. 2016, 177, 331–340. [Google Scholar] [CrossRef] [PubMed]
  46. Rawat, K.S.; Singh, S.K.; Bala, A.; Szabó, S. Estimation of crop evapotranspiration through spatial distributed crop coefficient in a semi-arid environment. Agric. Water Manag. 2019, 213, 922–933. [Google Scholar] [CrossRef]
  47. Jiang, C.; Zhang, L.B. Effect of ecological restoration and climate change on ecosystems: A case study in the Three-Rivers Headwater Region, China. Environ. Monit. Assess. 2016, 188, 382. [Google Scholar] [CrossRef]
  48. Huntington, T.G. Evidence for intensification of the global water cycle: Review and synthesis. J. Hydrol. 2006, 319, 83–95. [Google Scholar] [CrossRef]
  49. Khorchani, M.; Nadal-Romero, E.; Lasanta, T.; Tague, C. Carbon sequestration and water yield tradeoffs following restoration of abandoned agricultural lands in Mediterranean mountains. Environ. Res. 2022, 207, 112203. [Google Scholar] [CrossRef]
  50. Maurya, S.; Srivastava, P.K.; Gupta, M.; Islam, T.; Han, D. Integrating soil hydraulic parameter and microwave precipitation with morphometric analysis for watershed prioritization. Water Resour. Manag. 2016, 30, 5385–5405. [Google Scholar] [CrossRef]
  51. Wang, Y.F.; Lv, W.W.; Xue, K.; Wang, S.P.; Zhang, L.i.; Hu, R.H.; Zeng, H.; Xu, X.L.; Li, Y.M.; Jiang, L.L.; et al. Grassland changes and adaptive management on the Qinghai-Tibetan Plateau. Nat. Rev. Earth Environ. 2022, 3, 668–683. [Google Scholar] [CrossRef]
  52. Hou, L.L.; Xia, F.; Chen, Q.H.; Huang, J.K.; He, Y.; Rose, N.; Rozelle, S. Grassland ecological compensation policy in China improves grassland quality and increases herders’ income. Nat. Commun. 2021, 12, 4683. [Google Scholar] [CrossRef] [PubMed]
  53. Lu, F.; Hu, H.F.; Sun, W.J.; Zhu, J.J.; Liu, G.B.; Zhou, W.M.; Zhang, Q.F.; Shi, P.L.; Liu, X.P.; Wu, X.; et al. Effects of national ecological restoration projects on carbon sequestration in China from 2001 to 2010. Proc. Natl. Acad. Sci. USA 2018, 115, 4039–4044. [Google Scholar] [CrossRef] [PubMed]
  54. Guo, B.Y.; Wang, J.S.; Mantravadi, V.S.; Zhang, L.; Liu, G.Z. Effect of climate and ecological restoration on vegetation changes in the “Three-River Headwaters” region based on remote sensing technology. Environ. Sci. Pollut. Res. 2022, 29, 16436–16448. [Google Scholar] [CrossRef] [PubMed]
  55. Shao, Q.Q.; Liu, S.C.; Ning, J.; Liu, G.B.; Yang, F.; Zhang, X.; Niu, L.N.; Huang, H.B.; Fan, J.W.; Liu, J.Y. Assessment of ecological benefits of key national ecological projects in China in 2000–2019 using remote sensing. Acta Geogr. Sin. 2022, 77, 2133–2153. [Google Scholar]
  56. Huang, K.; Zhang, Y.J.; Zhu, J.; Liu, Y.J.; Zu, J.X.; Zhang, J. The Influences of Climate Change and Human Activities on Vegetation Dynamics in the Qinghai-Tibet Plateau. Remote Sens. 2016, 8, 876. [Google Scholar] [CrossRef]
  57. Chen, R.; Duan, K.; Shang, W.; Shi, P.; Meng, Y.; Zhang, Z. Increase in seasonal precipitation over the Tibetan Plateau in the 21st century projected using CMIP6 models. Atmos. Res. 2022, 277, 106306. [Google Scholar] [CrossRef]
  58. Teng, H.; Liang, Z.; Chen, S.; Liu, Y.; Rossel, R.A.V.; Chappell, A.; Yu, W.; Shi, Z. Current and future assessments of soil erosion by water on the Tibetan Plateau based on RUSLE and CMIP5 climate models. Sci. Total Environ. 2018, 635, 673–686. [Google Scholar] [CrossRef]
  59. Han, D.L.; Huang, J.P.; Ding, L.; Zhang, G.L.; Liu, X.; Li, C.; Yang, F. Breaking the Ecosystem Balance Over the Tibetan Plateau. Earth’s Future 2022, 10, EF002890. [Google Scholar] [CrossRef]
  60. Li, X.Y.; Long, D.; Scanlon, B.R.; Mann, M.E.; Li, X.D.; Tian, F.Q.; Sun, Z.L.; Wang, G.Q. Climate change threatens terrestrial water storage over the Tibetan Plateau. Nat. Clim. Change 2022, 12, 801–807. [Google Scholar] [CrossRef]
  61. Zhang, Q.; Fan, K.K.; Singh, V.P.; Song, C.Q.; Xu, C.Y.; Sun, P. Is Himalayan-Tibetan Plateau “drying”? Historical estimations and future trends of surface soil moisture. Sci. Total Environ. 2019, 658, 374–384. [Google Scholar] [CrossRef]
  62. Liu, Y.L.; Kumar, M.; Katul, G.G.; Feng, X.; Konings, A.G. Plant hydraulics accentuates the effect of atmospheric moisture stress on transpiration. Nat. Clim. Change 2020, 10, 691–695. [Google Scholar] [CrossRef]
  63. Bai, Y.F.; Guo, C.C.; Degen, A.A.; Ahmad, A.A.; Wang, W.Y.; Zhang, T.; Li, W.Y.; Ma, L.; Huang, M.; Zeng, H.J.; et al. Climate warming benefits alpine vegetation growth in Three-River Headwater Region, China. Sci. Total Environ. 2020, 742, 140574. [Google Scholar] [CrossRef]
  64. Gao, Q.; Guo, Y.; Xu, H.; Ganjurjav, H.; Li, Y.; Wan, Y.; Qin, X.; Ma, X.; Liu, S. Climate change and its impacts on vegetation distribution and net primary productivity of the alpine ecosystem in the Qinghai-Tibetan Plateau. Sci. Total Environ. 2016, 554, 34–41. [Google Scholar] [CrossRef]
  65. Cui, J.P.; Lian, X.; Huntingford, C.; Gimeno, L.; Wang, T.; Ding, J.Z.; He, M.Z.; Xu, H.; Chen, A.P.; Gentine, P. Global water availability boosted by vegetation-driven changes in atmospheric moisture transport. Nat. Geosci. 2022, 15, 982–988. [Google Scholar] [CrossRef]
  66. Wang, J.; Liu, D.S. Vegetation green-up date is more sensitive to permafrost degradation than climate change in spring across the northern permafrost region. Glob. Change Biol. 2022, 28, 1569–1582. [Google Scholar] [CrossRef]
  67. Di, Y.P.; Zhang, Y.J.; Zeng, H.; Tang, Z. Effects of changed Asian water tower on Tibetan Plateau ecosystem: A review. Bull. Chin. Acad. Sci. 2019, 34, 1322–1331. [Google Scholar]
  68. Zhang, J.L.; Zhou, L.L.; Ma, R.M.; Jia, Y.F.; Yang, F.; Zhou, H.Y.; Cao, X.Y. Influence of soil moisture content and soil and water conservation measures on time to runoff initiation under different rainfall intensities. Catena 2019, 182, 104172. [Google Scholar] [CrossRef]
  69. Gao, J.; Li, F.; Gao, H.; Zhou, C.B.; Zhang, X.L. The impact of land-use change on water-related ecosystem services: A study of the Guishui River Basin, Beijing, China. J. Clean. Prod. 2017, 163, S148–S155. [Google Scholar] [CrossRef]
  70. Zhang, B.; Li, W.H.; Xie, G.D.; Xiao, Y. Water conservation of forest ecosystem in Beijing and its value. Ecol. Econ. 2010, 69, 1416–1426. [Google Scholar]
  71. Li, S.H.; Yu, D.Y.; Huang, T.; Hao, R.F. Identifying priority conservation areas based on comprehensive consideration of biodiversity and ecosystem services in the Three-River Headwaters Region, China. J. Clean. Prod. 2022, 359, 132082. [Google Scholar] [CrossRef]
  72. Jia, G.Y.; Hu, W.M.; Zhang, B.; Li, G.; Shen, S.Y.; Gao, Z.H.; Li, Y. Assessing impacts of the Ecological Retreat project on water conservation in the Yellow River Basin. Sci. Total Environ. 2022, 828, 154483. [Google Scholar] [CrossRef] [PubMed]
  73. Liu, G.; Shao, Q.; Fan, J.; Huang, H.; Liu, J.; He, J. Assessment of Restoration Degree and Restoration Potential of Key Ecosystem-Regulating Services in the Three-River Headwaters Region Based on Vegetation Coverage. Remote Sens. 2023, 15, 523. [Google Scholar] [CrossRef]
  74. Wang, S.Z.; Fan, J.W.; Li, Y.Z.; Wu, D.; Zhang, Y.X.; Huang, L. Dynamic response of water retention to grazing activity on grassland over the Three River Headwaters region. Agric. Ecosyst. Environ. 2019, 286, 106662. [Google Scholar] [CrossRef]
  75. Ning, X.G.; Zhu, N.; Liu, Y.F.; Wang, H. Quantifying impacts of climate and human activities on the grassland in the Three-River Headwater Region after two phases of Ecological Project. Geogr. Sustain. 2022, 3, 164–176. [Google Scholar] [CrossRef]
Figure 1. Study area location and natural environment in the Three−Rivers Headwater Region ((a) The study area location and land use types. (b) Annual average precipitation. (c) Annual average potential evapotranspiration. (d) Annual average leaf area index. (e) Slope, (f) Ecosystem zones: sub−cold sub−humid Guoluo−Naqu Plateau mountain alpine shrub–meadow region (HIB1), sub−cold semi−arid Qiangtang Plateau lake basin alpine steppe region (HIC2), sub−cold semi−arid southern Qinghai alpine meadow–steppe (HIC1), temperate semi−arid eastern Qinghai−Qilian montane steppe zone (HIIC2), sub−cold arid Kunlun high mountain and plateau alpine desert region (HID1), temperate arid North Kunlun mountain desert region (HIID2), temperate arid Ngali mountain desert region (HIID3), temperate humid/sub−humid Western Sihuan and Eastern Xizang high mountain and deep valley coniferous forest region (HIIAB1)).
Figure 1. Study area location and natural environment in the Three−Rivers Headwater Region ((a) The study area location and land use types. (b) Annual average precipitation. (c) Annual average potential evapotranspiration. (d) Annual average leaf area index. (e) Slope, (f) Ecosystem zones: sub−cold sub−humid Guoluo−Naqu Plateau mountain alpine shrub–meadow region (HIB1), sub−cold semi−arid Qiangtang Plateau lake basin alpine steppe region (HIC2), sub−cold semi−arid southern Qinghai alpine meadow–steppe (HIC1), temperate semi−arid eastern Qinghai−Qilian montane steppe zone (HIIC2), sub−cold arid Kunlun high mountain and plateau alpine desert region (HID1), temperate arid North Kunlun mountain desert region (HIID2), temperate arid Ngali mountain desert region (HIID3), temperate humid/sub−humid Western Sihuan and Eastern Xizang high mountain and deep valley coniferous forest region (HIIAB1)).
Atmosphere 14 01453 g001
Figure 2. Conceptual diagram of the water balance model used in the InVEST water yield model [31]. Only parameters shown in color are included, and parameters shown in grey are ignored.
Figure 2. Conceptual diagram of the water balance model used in the InVEST water yield model [31]. Only parameters shown in color are included, and parameters shown in grey are ignored.
Atmosphere 14 01453 g002
Figure 3. Evaluation of simulated total water yield (TWY) versus observed TWY (TWY was simulated by the InVEST model using (a) calibration Kc values with interannual variation in vegetation, (b) Kc values obtained using the relationship Kc = LAI/3, or (c) Kc coefficients estimated at Kc = 0.65. Shadow represents a 95% confidence interval).
Figure 3. Evaluation of simulated total water yield (TWY) versus observed TWY (TWY was simulated by the InVEST model using (a) calibration Kc values with interannual variation in vegetation, (b) Kc values obtained using the relationship Kc = LAI/3, or (c) Kc coefficients estimated at Kc = 0.65. Shadow represents a 95% confidence interval).
Atmosphere 14 01453 g003
Figure 4. Implementation flow chart of this study (Formula a is derived from Zhao, et al. [37]. Formula b is derived from Allen, et al. [35]. Formula c is derived from Sharp, et al. [31]).
Figure 4. Implementation flow chart of this study (Formula a is derived from Zhao, et al. [37]. Formula b is derived from Allen, et al. [35]. Formula c is derived from Sharp, et al. [31]).
Atmosphere 14 01453 g004
Figure 5. Regional corrected Kc coefficient of grassland in the Three-Rivers Headwater Region from 1990 to 2020.
Figure 5. Regional corrected Kc coefficient of grassland in the Three-Rivers Headwater Region from 1990 to 2020.
Atmosphere 14 01453 g005
Figure 6. Temporal variation in water conservation in the Three-Rivers Headwater Region from 1991 to 2020 (relative to the 1991–2020 anomaly).
Figure 6. Temporal variation in water conservation in the Three-Rivers Headwater Region from 1991 to 2020 (relative to the 1991–2020 anomaly).
Atmosphere 14 01453 g006
Figure 7. The spatial distribution (a) and trend (b) of water conservation in the Three-Rivers Headwater Region from 1991 to 2020 (the inset panels in the bottom right of (a) display the water conservation values in different ecosystem zones using a violin diagram. The inset panels in the bottom right of (b) indicate the significance level (p < 0.05). The percentages of increasing (I) and decreasing (D) trends (percentage of significant trends in parentheses) are shown at the bottom of (b)).
Figure 7. The spatial distribution (a) and trend (b) of water conservation in the Three-Rivers Headwater Region from 1991 to 2020 (the inset panels in the bottom right of (a) display the water conservation values in different ecosystem zones using a violin diagram. The inset panels in the bottom right of (b) indicate the significance level (p < 0.05). The percentages of increasing (I) and decreasing (D) trends (percentage of significant trends in parentheses) are shown at the bottom of (b)).
Atmosphere 14 01453 g007
Figure 8. Interannual variability of climatic and vegetation elements from 1991 to 2000 ((a) Average annual precipitation. (b) Average annual potential evapotranspiration. (c) Annual average leaf area index).
Figure 8. Interannual variability of climatic and vegetation elements from 1991 to 2000 ((a) Average annual precipitation. (b) Average annual potential evapotranspiration. (c) Annual average leaf area index).
Atmosphere 14 01453 g008
Figure 9. Explanatory power of interactive influencing factors on water conservation in the Three-Rivers Headwater Region and its sub-regions ((ad): LAI represents vegetation factor, (eh): NDVI represents vegetation factor).
Figure 9. Explanatory power of interactive influencing factors on water conservation in the Three-Rivers Headwater Region and its sub-regions ((ad): LAI represents vegetation factor, (eh): NDVI represents vegetation factor).
Atmosphere 14 01453 g009
Table 1. Basic information regarding the datasets used in this research.
Table 1. Basic information regarding the datasets used in this research.
Data TypeData FormatData Source
Meteorological data (precipitation, temperature, relative humidity, wind speed)Point-based, daily, from 1991 to 2020National Climate Center of the China Meteorological Administration (http://data.cma.cn/, accessed on 1 September 2022). A 500 m resolution grid is interpolated by using the professional meteorological interpolation software ANUSPLINA version 4.4 (http://fennerschool.anu.edu.au/files/anusplin44.pdf, accessed on 1 September 2022)
LAIGrid, 5 km resolution, from 1991 to 2020Global Mapping (GLOBMAP) LAI version 3 dataset (https://zenodo.org/record/4700264, accessed on 1 September 2022)
NDVIGrid, 1 km resolution, from 1991 to 2020https://zenodo.org/record/6295928, accessed on 1 September 2022
Soil propertiesGrid, 1 km resolution, 2000Harmonized World Soil Database (HWSD) v1.2 (http://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/, accessed on 1 September 2022)
Soil depthGrid, 1 km resolution, 2000Soil Data Center, National Earth System Science Data Sharing Infrastructure, National Science and Technology Infrastructure of China (http://soil.geodata.cn, accessed on 1 September 2022)
Land use (1990, 1995, 2000, 2005, 2010, 2015, 2020)Grid, 30 m resolution, 2000, 2005, 2010, 2015, and 2020Resources and Environmental Sciences Data Center (RESDC), Chinese Academy of Sciences (http://www.resdc.cn, accessed on 1 September 2022)
Digital elevation model (DEM)Grid, 1 km resolution, 2020Resources and Environmental Sciences Data Center (RESDC), Chinese Academy of Sciences (http://www.resdc.cn, accessed on 1 September 2022)
Table 2. Correlation between climate, vegetation, and water conservation.
Table 2. Correlation between climate, vegetation, and water conservation.
Driving Element1991–20002000–20201991–2020
Precipitation0.92 **0.97 **0.97 **
Potential evapotranspiration0.13–0.5 *–0.19
LAI0.070.360.33
Note: ** indicates a significant correlation at the 0.01 level (bilateral); * indicates a significant correlation at the 0.05 level (bilateral).
Table 3. The q values of influencing water conservation factors in the Three-Rivers Headwater Region.
Table 3. The q values of influencing water conservation factors in the Three-Rivers Headwater Region.
(1) Influencing FactorsAnnual Average PrecipitationAnnual Average Potential EvapotranspirationLAISlopeAltitude
q value0.73 **0.207 **0.582 **0.109 **0.221 **
(2) Influencing FactorsAnnual Average PrecipitationAnnual Average Potential EvapotranspirationNDVISlopeAltitude
q value0.746 **0.185 **0.568 **0.116 **0.242 **
Note: ** represents significance level p < 0.001.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pan, Y.; Yin, Y. Spatial and Temporal Evolution Characteristics of Water Conservation in the Three-Rivers Headwater Region and the Driving Factors over the Past 30 Years. Atmosphere 2023, 14, 1453. https://doi.org/10.3390/atmos14091453

AMA Style

Pan Y, Yin Y. Spatial and Temporal Evolution Characteristics of Water Conservation in the Three-Rivers Headwater Region and the Driving Factors over the Past 30 Years. Atmosphere. 2023; 14(9):1453. https://doi.org/10.3390/atmos14091453

Chicago/Turabian Style

Pan, Yao, and Yunhe Yin. 2023. "Spatial and Temporal Evolution Characteristics of Water Conservation in the Three-Rivers Headwater Region and the Driving Factors over the Past 30 Years" Atmosphere 14, no. 9: 1453. https://doi.org/10.3390/atmos14091453

APA Style

Pan, Y., & Yin, Y. (2023). Spatial and Temporal Evolution Characteristics of Water Conservation in the Three-Rivers Headwater Region and the Driving Factors over the Past 30 Years. Atmosphere, 14(9), 1453. https://doi.org/10.3390/atmos14091453

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop