Next Article in Journal
Long-Term Monitoring and Research in Forest Hydrology: Towards Integrated Watershed Management
Next Article in Special Issue
Policy Implementation for Water Resources Protection: Assessing Spatio-Temporal Trends of Results from Process-Based Outcomes of Resource-Directed Measures Projects in South Africa
Previous Article in Journal
Identifying Thresholds, Regime Shifts, and Early Warning Signals Using Long-Term Streamflow Data in the Transboundary Rio Grande–Rio Bravo Basin
Previous Article in Special Issue
Sewage Irrigation Fields—From Relict Landscape to Blue-Green Urban Infrastructure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ecological–Economic Assessment and Managerial Significance of Water Conservation in the Headwaters of the Yellow River

1
State Key Laboratory of Grassland Agro-Ecosystems, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730000, China
2
Key Laboratory of Grassland Livestock Industry Innovation, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730000, China
3
College of Life Science and Engineering, Hexi University, Zhangye 734000, China
4
College of Agriculture and Husbandry, Qinghai University, Xining 810016, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(16), 2553; https://doi.org/10.3390/w14162553
Submission received: 12 July 2022 / Revised: 10 August 2022 / Accepted: 13 August 2022 / Published: 19 August 2022
(This article belongs to the Special Issue Water ​Management and ​Environmental Protection)

Abstract

:
The water conservation function plays a vital role in the land–water cycle. As the “Chinese water tower”, the headwaters of the Yellow River are of great significance to the safety of the Yellow River basin and even the global ecosystem. Taking the grassland ecosystem in the Yellow River source area as the research object, the InVEST water yield model with modified parameters and the ecological value evaluation of the modified equivalent factor method were used to explore the simulated spatio-temporal changes and the value of grassland water conservation from 2001 to 2020. The results show that: (1) the average total amount of water conservation in the source area is 549 × 108 m3, which is 16% of the runoff in the Yellow River basin, with a growth rate of 7.5 mm/year 1 and a contribution rate of 30%; (2) the total ecological value of grassland water conservation in 2020 is USD 340.03 × 108. The proportion of improved grassland in ecological restoration and management is only 0.51%, while the proportion of original alpine meadow reaches 67% and its ecological function and value are irreplaceable; (3) based on the comprehensive indicators of water conservation capacity, value and importance, Qumalai, Chengduo and Maduo counties are ranked as priority areas for the ecological protection of water resources.

1. Introduction

Over the past decades, a substantial contradiction concerning the development of human society and the protection of the ecological environment has become increasingly exacerbated and 60% of global ecosystems have been degraded to varying degrees [1]. The headwater of the Yellow River (HYR) is a major ecological and functional area for water conservation in China. The safety of the ecosystem is related to the ecological, environmental and economic development of the Yellow River basin and even East Asia, focusing the interest of scientists from all walks of life [2]. However, the vegetation type of the ecosystem in the headwater is mainly grassland, which is complex, diverse and highly fragile. Due to the stress of global climate change and the influence of anthropological factors, the ecological environment has been damaged and the grasslands’ degradation is serious; this leads to a reduction in the water conservation capacity, deterioration of the soil and production loss, all of which impact the lifestyle of local herdsmen in the region [3,4,5]
For the regional rotational grazing and grass stand management, a comprehensive grazing ban was implemented on moderately degraded grassland in the HYR, under the Qinghai Sanjiangyuan Nature Reserve Construction Phase I and Phase II. The project achieved remarkable results and the ecological situation has improved significantly [6,7,8]. The existing ecological protection projects and policies are often based on the overall regional planning and the grasslands are not managed in a classified way. Therefore, a water conservation function should be evaluated based on the grassland type, the ecological function and the ecological value, while the priority conservation areas should be determined.
The assessment is based in the northeast of the Qinghai-Tibet Plateau location. Due to the high altitude and complex terrain, it is difficult to obtain scientific research data. With the continuous development of geographic information systems (GIS), hydrological model simulation becomes the preferred method for estimating water yield [9,10,11]. Common hydrologic models include: SWAT (Soil and Water Assessment Tool) [12,13], ARIES (the artistic intelligence for ecosystem services) [14] and InVEST (Integrated Value of Ecosystem Services and Tradeoffs). The InVEST model is a widely used and mature open-source ecosystem modeling tool. It is superior to other hydrological models in terms of spatial analysis function and valuation accuracy [15,16,17]. The climate season factor Z of the InVEST model, which enters the model, is an empirical constant. The value of Z is uncertain and ranges from 1 to 30 [18]. For example, when assessing the water yield of 22 catchments in Scotland and Wales, the value of Z is 30 [19]. For calculating the water retention of Ghana and Côte d’Ivoire in West Africa, the value of Z is 17 [20], whereas for the assessment of water production in Kentucky, the value of Z is 30 [21]. For the assessment of water yield in the Danjiang River basin, China, the value of Z is 12.6 [22]. The InVEST model simulates the water conservation quantity in the HYR, and the spatial and temporal dynamic analysis of the water storage quantity index in the source area still needs more rigorous and scientific methods. The Mann–Kendall (MK) trend analysis test method has been mostly reported to analyze hydrological–meteorological element changes [23,24,25]. It is a nonparametric statistical test that handles the environmental time series, judges the significance of the trends and estimates the mutation point time at which the simulated values appear. Lepeska et al., recognized that urbanization has a potential impact on water resources’ loss. From 1990 to 2018, water resources decreased year by year with the development of urbanization [26]. Irannezhad et al., used the Mann–Kendall test climate change effects on the hydrological process of snow cover in Finland over the past 100 years, and it is believed that the temperature change had the most significant effect on the snow water equivalent in Finland (p < 0.05) [27].
The ecological assessment of water yield is an important reference for investment in ecological conservation projects. Currently, the ecosystem service value (ESV) is mainly calculated by the functional value method and the equivalence factor method [28,29,30]. The equivalent factor method can be used to estimate the ecological value of different ecosystems and it is more reasonable to apply this method to the study of the ecosystem value of grassland [31,32,33,34,35]). The equivalence factor method is based on the IGBP classification method [36,37,38,39,40,41,42], making the evaluation results more accurate and conducive to the formulation and implementation of policy suggestions. At present, some scholars have adopted this method to evaluate the ecological value and to achieve reliable results. Xie et al., calculated that the ecological service value of the Qinghai-Tibet Plateau ecosystem is USD 596.48 × 1010/year−1, accounting for 17.68% of the annual service value of the national ecosystem, 0.61% of the global ecosystem and 16.5% of the water conservation value [29]. The water conservation value of grassland in the HYR is USD 15.15 × 108, grassland with medium cover (12.53%) and grassland with high cover (11.93%) [43]. Based on the improved equivalence factor method of the net primary productivity (NPP) of vegetation, the total value of different types of ecosystem services in China in 2010 was calculated as USD 242 × 1012, and the national equivalence factor value is USD 21,699.41/hm2 [44]. However, for the regional averaging of the national equivalence factor in the HYR, it is necessary to take into account the vegetation cover in the study area, to improve the equivalence factor value based on the NPP data and to develop the reliability of the equivalence factor method at the regional scale [45].
Using the alpine meadow in the HYR, three key problems will be solved using the modified parameter investment water yield model and the modified equivalent factor method: (1) the spatio-temporal pattern of water yield distribution; (2) the evaluation of the ecological function value of water conservation; (3) the comprehensive evaluation of the water conservation function based on grassland type, ecological function and ecological value, and the determination of an ecological priority conservation area to theoretically support the ecological conservation and management of the headwater area.

2. Materials and Methods

2.1. Study Area

The HYR (32°30′ N–35°0′ N—95°50′ E–103°30′ E) is the largest runoff area and accounts for 16.20% of the total Yellow River basin, with a drainage area of 12.37 × 104 km2 [46]. The study area is between 2680 and 6248 m above sea level (a.s.l.) and increases gradually with an average elevation of 4123 m (a.s.l.) from east to west. The area under study characterizes a typical continental highland climate and the annual average temperature ranges from −12.7 °C to 5.6 °C. Under the influence of the southwest monsoon, the annual average precipitation ranges from 281.8 mm to 1058.6 mm. However, the summer (June to September) generally accounts for 90% of the total annual precipitation. The average annual evaporation ranges from 730 mm to 1700 mm. In the HYR, there are various types of grassland, of which the alpine meadow accounts for 51.55% of the total area of the source area, while improved grassland accounts for only 0.01%. The data of vegetation cover types involved in this study generated the spatial distribution map of vegetation types in the study area based on the IGBP and the Chinese grassland type survey map (CERES land classification from the International Geosphere Biosphere Program IGBP) (Figure 1).

2.2. Data Source and Preprocessing

Environmental variables and socioeconomic data are needed to evaluate water yield and ecological value for the time period 2001–2020 (Table 1). The environmental variables’ data are divided into four categories: soil physical and chemical factors (S, 9), terrain factors (T, 5), meteorological factors (A, 4) and vegetation factors (B, 2). The raster images of all variable factors are projected and reclassified according to the geodetic world system 84 (WGS-1984), socioeconomic data (E, 2).

2.3. Methods

2.3.1. Spatial and Temporal Pattern and Trend Analysis of Water Yield

Water Yield Calculation

Using the water yield module of the InVEST model, the total annual water yield (mm) was calculated by formula (Equation (1)), where AET (x) is the annual actual evapotranspiration of grid unit X, and P (x) is the annual precipitation of grid unit X. The formula is as follows. AET (x)/P (x) is the approximate Budyko curve (Equation (2)) in terms of annual precipitation and actual evapotranspiration [47]. ω (x) is a nonphysical empirical adjustment parameter that describes the natural climate and soil properties of the watershed (Equation (3)) [48]. The Zhang coefficient is a climatic seasonal factor that represents local hydrogeological properties and is an empirical constant. The standard value ranges from 1 to 30 [18,19]. PAWC effective water retention capacity (Equation (4)), where sand, silt, clay and soil organic matter (om) are described in percent (%).
Y ( x ) = ( 1     AET ( x ) P ( X ) )   ×   P ( x )
AET ( x ) P ( X ) = 1 + ω x R xj 1 + ω x R xj + 1 R xj
ω x = z   ×   ( AWC x P x )
AWC = 54.509     0.132   ×   sand     0.003   ×   sand 2     0.055   ×   silt     0.006   × silt 2     0.738   ×   clay     0.007   ×   clay 2     2.688   ×   om + 0.501   ×   om 2
Calculation formula for water retention (mm) (Equation (5)), topographic index (TI), soil saturated hydraulic conductivity (Ksat) and velocity coefficient of total annual water storage [49]. The calculation formula is as follows:
Retention = min ( 1 ,   249 Velocity )   ×   min ( 1 ,   0.93   ×   TI 3 )   ×   min ( 1 ,   K sat 300 )   ×   Y ( x )

Water Retention Spatio-Temporal Dynamics and Trend Analysis

The InVEST model is based on ArcGIS and MATLAB 2017 software for analyzing the annual water retention and the spatio-temporal dynamic pattern of the HYR from 2001 to 2020 and for conducting geographic mapping. The trend analysis uses Theil Sen (SEN) median trend analysis and nonparametric (Mann–Kendall, MK) to test each pixel to determine water retention trend and mutation point detection. When the Sen slope > 0, it shows an upward trend, whereas in the opposite situation it shows a downward trend (SEN slope < 0). With a confidence level α = test the significance of the change trend of water retention at 0.05 and confidence level α = 0.05 to test the significance of the change trend of the water yield. If z > 1.96, the results increase significantly; if z < −1.96, the results decrease significantly; if −1.96 ≤ z ≤ 1.96, the results do not change significantly (Table 2).
The Mann–Kendall method (MK) estimates the point of annual change in the simulation value of the annual water production from 2000 to 2020. UFk is a standard normal distribution (Equation (6)), calculating the statistical sequence by x1, x2, …, xn; UBk= −UFk generates the inverse order xn, xn−1, …, x1 (Equation (7)).
UF k = S k E ( S k ) ] VAR ( S k )
UB k = UF k ,   k = n ,   n 1 ,   ,   1 ,   UB 1 = 0
Confidence level α = 0.05 to check the change trend of time series of water retention. When UF or UB > 0, it indicates that the series shows an upward trend, and when UF or UB < 0, it shows a downward trend. The time point corresponding to the intersection of UF and UB curves is the time point at which mutation starts.

2.3.2. Ecological Value Evaluation of Regional Water Yield

Based on the Chinese terrestrial ecological service value equivalent table [30], the ecological service value equivalent factor in the HYR is calculated using the national average grain and unit price data from 2001 to 2020, according to the principle that the economic value of a standard ecosystem service equivalence factor is 1/7 of the economic value of food production per unit area of farmland [35,42,50]. Based on the NPP data of grassland in the source area of the Yellow River, an ESV equivalent factor of USD 14,461.68/(hm2·year−1) was calculated for the study area and the total value of aquatic ecological services of grassland in the HYR was calculated in combination with the distribution area of grassland species in the headwaters.

2.3.3. Division of Ecological Priority Protection Areas

For the delineation method of the water conservation area, we followed the guidelines for the delineation of the ecological red line of protection [51,52,53,54]. Since the delineation of ecological red line certainty, only the water protection amount and importance level in 2020 were extracted and superimposed to delineate the ecological red line and the protection area.
According to the evaluation and classification steps of the guidelines, the grid values corresponding to 30%, 50%, 80% and 90% of the total water conservation capacity are classified as evaluation and classification nodes of water conservation capacity. The analytical hierarchy process (AHP) is used to build the hierarchical structure of importance rating. At the same time, the significance of the data is standardized using the correlation test between the variable environmental factors and the water production. The AHP is used to build the evaluation matrix. After the consistency test, the weighting of the influencing factors is determined. Combined with the comprehensive index algorithm, the water retention importance index (WRI) is obtained. The value range of the WRI ranges from 0–1 and the closer it is to 1, the higher the importance. Finally, the overlay analysis in ArcGIS is used and the natural segment points method is employed for classification.

3. Results

3.1. Spatial and Temporal Distribution of Water Yield

3.1.1. Spatio-Temporal Variation Characteristics

The annual average water retention from 2001 to 2020 in the study area was 443.86 mm, and the overall interannual variation trend of water conservation showed a slight increase, with a rate of 7.5 mm/year −1. The total annual average water retention was 5.49 × 1010 m3.
The estimated water retention was mainly between 300 and 700 mm, which accounted for 93.65% of HYR, and the areas with less than 300 mm and more than 700 mm accounted for 3.21% and 3.15%, respectively. The high value areas were located in Maduo, Toronto and Qumai counties, with heavy rainfall and narrow alpine grassland areas. After rainfall, it is easy to form surface runoff and the regional water retention capacity is high. The median water conservation was observed in Ruoergai, Aba and Hongyuan counties, with water retention of more than 600 mm/year −1, whereas Maqu and Jiuzhi counties also showed high water retention, at 540 mm and 587.79 mm/year −1, respectively. The median area has a high forest cover and the forest canopy, litter and surface soil effectively intercept the rainfall, thus reducing the water conservation capacity. The least water retention of Qumalai county was 229.69 mm/year- 1, with less rainfall, flat terrain and dense forest vegetation coverage intercepting rainfall and surface runoff, resulting in less water retention (Figure 2).

3.1.2. Trend Change Characteristics of Water Yield

The water conservation capacity of most regions in the HYR showed an upward trend from 2001 to 2020 (Figure 3). The area with significant growth accounted for 9.44%, mainly in the southeast of Zeku and Henan counties, the middle of Maqin and the north of Qumalai, with 90.56% of regional water conservation, which showed a slight upward trend.
The water retention in the source area in the last 20 years is statistically described in four time periods (Table 3) and the water production in each period shows an increasing trend. Among them, the annual water production from 2015 to 2020 is the largest. It should be noted that the water production of Chengduo, Gande, Maqu, Dari, Ruoergai, Jiuzhi, Aba and Hongyuan counties increased consistently from 2001 to 2020, while that of the other counties increased from 2001 to 2010, decreased from 2011 to 2015 and increased from 2016 to 2020.
Analysis of the water retention data by the MK test shows that 2005 is a mutation year. Therefore, we take 2005 as the cut-off point to analyze the change in water retention (Figure 3). The annual water retention increased from 347.14 mm in 2001 to 510.17 mm in 2005 and the lowest annual water retention was recorded in 2002 with 306.06 mm. From 2005 to 2015, the annual water production was 435.77 mm, with a growth rate of 2.71 mm/year −1. Another abrupt change occurred in 2015, which shows the water production increased significantly from 435.77 mm to 518.75 mm. In 2018, the water yield in the HYR was the highest since 2001. Overall, the simulation results show that the water yield in the HYR continues to increase (Figure 4).
In general, the spatial and temporal changes of water storage in grassland ecosystems are related to rainfall, vegetation cover, surface runoff and local ecological governance policies, and have seen high values in 2005 and 2015 [55,56,57,58]. The ecological protection project of the HYR is an important part of the Sanjiangyuan ecological protection program. The “Sanjiangyuan ecological protection project of China” is divided into two phases, every eight years; 2005 and 2014 are the two phases of the project, respectively. Affected by the policy, the ecological environment in the source area has been significantly improved, the vegetation coverage has increased, the overall production and water storage capacity has shown a steady growth trend, the ecological service function of water conservation in the source area has been improved and the ecological management has achieved remarkable results.

3.2. Ecological Service Value of Water Conservation in Grassland Ecosystem

The service value of grassland water conservation in the HYR is divided into five intervals: 2001, 2005, 2010, 2015 and 2020. The total ecological value was USD 173.06 × 108, USD 247.20 × 108, USD 296.17 × 108, USD 296.26 × 108 and USD 340.03 × 108 during 2001, 2005, 2010, 2015 and 2020, respectively, with the value of water conservation services increasing steadily over the past 20 years. The ecological value of water conservation apparently differs between different grassland types. The low-value area improves grassland, but the value of ecological service also shows a weak growth trend from USD 3.3 × 108 in 2001 to USD 5.74 × 108 in 2020, with an annual contribution rate of 5%, distributed in the north of Zeku County. Alpine grassland in the high-value area was from USD 146.15 × 108 in 2001 to USD 227.59 × 108 in 2020, with an annual contribution rate of 80%. From high to low, the average ecological values and contribution rates over the years are: alpine meadow (66.93%) > temperate grassland (11.30%) > mountain meadow (7.85%) > lowland meadow (7.19%) > temperate desert (5.17%) > alpine grassland (1.03%) > improved grassland (0.51%) (Figure 5).
The value of water conservation in grassland ecosystems is increasing year after year, which further illustrates the importance of ecological protection. Monetary value cannot only make people intuitively feel the ecological value and improve the awareness of ecological protection but can also provide a reference for the formulation of local ecological protection and management policies and planning and design according to ecological values and land uses. The first phase of the ecological protection project invested RMB 7.5 billion, reaching one third of the assessed value in 2005; the second phase investment reached RMB 16 billion, reaching one half of the assessed value in 2015. Ecological restoration planning was classified and managed according to different ecosystems and land use types with remarkable achievements.

3.3. Spatial Division of Water Conservation Capacity and Its Importance

Five levels of water conservation capacity were identified. The lowest value area (grade I) is concentrated in the northwest of Qumalai, Maduo and Chengduo counties; the low-value area (grade II) is in the southeast of Maduo and Chengduo counties and the northwest of Xinghai and Maqu counties. The medium-grade area (grade III) is in the northeast of Gande, Maqin, Dari and Xinghai counties; the high-grade area (grade IV) is in the northwest of Zeku, Henan, Maqu and Jiuzhi counties. The highest areas (grade V) are in Ruoergai, Aba and Hongyuan counties and account for 14.30%, 13.40%, 21.40%, 34.20% and 16.70% of the HYR, respectively. In general, the spatial distribution of water yield in the source area of the Yellow River is imbalanced and decreases from the southeast to the northwest, showing a step-like distribution (Figure 6).
Five levels of importance were identified: general important area (grade I; <127.5 mm), slightly important area (grade II; 127.5–202.5 mm), moderately important area (grade III; 202.5–340.0 mm), highly important area (grade IV; 340.0–382.5 mm) and extremely important area (grade V; >382.5 mm).
The patterns of the water conservation importance are: the northwest of Qumalai, Maduo and Chengduo counties (mainly grade I); the northwest of Qumalai, Maduo and Chengduo districts (mainly grade II); the southeast of Maduo and Chengduo counties and the northwest of Xinghai and Maqu counties (mainly grade III); the higher-grade area (grade IV) is located in the northwest of Zeku, Henan, Maqu and Jiuzhi counties. The highest areas (grade V) are concentrated in Ruoergai, Aba and Hongyuan counties (mainly grade V) and account for 14.30%, 13.40%, 21.40%, 34.20% and 16.70% of the HYR, respectively. In general, the spatial division of water conservation importance in the HYR is uneven and decreases from the southeast to the northwest, showing a step-like distribution (Figure 7).
Overall, the spatial division of water conservation capacity and its importance is an indicator of the ecosystem assessment of water conservation capacity. In this study, by dividing the importance level of the water conservation capacity, the spatial mapping of the water conservation function facilitates the division and management of the regional ecological reserves.

4. Discussion

4.1. Regionalization of Model Parameters

Regionalization of Z-Value Parameters in the Model

At present, only a few researchers have used the water production module of the InVEST model to analyze the ecological service function of water conservation in the HYR. Yang et al., used the InVEST model to obtain the total water production and to trend change dynamics of the grassland in the Yellow River basin from 1990 to 2018. Our results demonstrated that, when Z is taken as 3.6, the average annual water production is simulated as (534 × 108 m3) and the actual average water production in 1995, 2005 and 2018 was (596 × 108 m3), which is relatively close to their study results [59]. A recent study conducted by Zhao et al., calculated the water volume of grassland in the source area of the Yellow River, and the value of Z is 3.5. It is obtained that the average water production for many years from 2001 to 2020 is 509 × 108 m3, which is similar to the results of the present study (549 × 108 m3) [60]. Kang et al., analyzed the annual water production and its change trend in the source area of the Yellow River from 1990 to 2010. They calculated the value of Z as 5, and the simulated total water production in 2010 was 0.27 × 108 m3, lower than the actual total water production in 2010 (497.27 × 108 m3). The reason for this unreasonable value is that the simulated and measured values are not verified [61]; therefore, Z = 5 is not applicable to the source area of the Yellow River. In this study, based on the annual average water resources of the Tangnaihai hydrological station in the Yellow River Water Resources Bulletin [62], the water yield results of the InVEST model were repeatedly verified (Figure 8). The InVEST model is based on ArcGIS software for data processing and mapping of the environmental variable grid. We take Z = 3.5 as the starting value [60], design the difference (every 0.5 is an interval and gradually increase the difference to run the model. When Z = 5, the average annual water production in 2001, 2005, 2010 and 2020 is close to the value of the actual water production in the bulletin (Figure 8).

4.2. Improvement Measures of Water Conservation Function

The contribution rate of improved grassland to the value of water conservation function in the source area is only 0.51%, which plays a greater role in maintaining biodiversity and productivity, while the restoration ability of water conservation and ecological function is lower [63]. On the other hand, the original alpine grassland plays a significant role, ensuring 60% of the water conservation and ecological value of the source area. Ecological protection projects and policies play a positive role in the restoration of regional ecological functions [59]. The previous “one size fits all” management model is no longer applicable. Due to the change of management ideas, the goal of restoring grassland productivity is changed to the goal of restoring ecological functions. According to the classification standard of ecological service value, ecological management is being carried out on the grassland in the study area and the high-value area is fenced. Generally, the vegetation coverage has increased by 10% and the biomass has increased by 50g/m2 in the past three years, which can effectively enhance the water conservation function [64]. The median area is divided into rotational grazing and grassland livestock balance and the appropriate livestock carrying capacity of herdsman. The grasslands in the warm season (livestock carrying capacity = their own grassland area/10.62) and the cold season (livestock carrying capacity = their own grassland area/25.21) are checked and verified, respectively, in order to curb the phenomenon of overgrazing and to reduce the ecological pressure of the grasslands [65,66]. In low-value areas, grassland fertilization can increase the content of soil organic carbon and phosphorus, change the structure of soil particle size, increase soil water content, reduce soil erosion and improve the water conservation function of grassland [66].

4.3. Ecological Priority and Regional Management

From 2005 to 2010, the implementation of China’s “phase I project of ecological protection and construction in the HYR” effectively curbed the trend of ecological degradation in the HYR and the ecological situation in the source area showed a stable improvement. In 2018, the water production of the source area was 498.79 mm, reaching the highest value in the recent 20 years. In January 2018, Sanjiangyuan National Park, the source of the Yellow River Park, carried out key ecological restoration projects. According to the remote sensing monitoring results of 2020, the vegetation coverage of sandy land in the source area increased from 33.36% to nearly 40%. The wetland vegetation coverage increased by more than 4%, the grassland vegetation coverage reached 74%, the grass yield reached 3082 kg/hm2, the water yield reached 587.5 mm and the water conservation function was significantly improved [67]. However, the ecological environment in the source area of the Yellow River has a stage of overall improvement and local degradation. Based on the grassland type, ecological function and ecological value, evaluating the water conservation function determines the priority protection areas and determines that Qumalai, Chengduo and Maduo counties are in extremely sensitive areas. The National Nature Reserve of the Yellow River source park is established in Maqu County to expand the ecological reserve to Qumalai and Chengduo counties, with Maqu County as the center. These three counties and districts are set as the priority protection areas for the ecological functions of water conservation.

5. Conclusions

The stability of the water conservation function and the evaluation of the ecological value in the HYR are related to the protection and management of water resources in the Yellow River basin. The quantitative analysis of water retention and the management of ecological spatial planning are the focus of the research. Based on the calculation of water retention using improved parameters of the InVEST model, this study analyzes the spatio-temporal dynamic characteristics of water conservation in the headwaters of the Yellow River from 2001 to 2020; the modified NPP equivalent factor method is used to evaluate the ecological value of the water production services and to classify the ecological conservation areas. Finally, the following three conclusions are drawn:
From 2001 to 2020, the average annual water yield of the study area was 443.86 mm and the average annual total water yield was 4.29 × 1010 m3, with an increase in water retention in most areas. The area with significant growth accounted for 9.44% and is mainly located in the southeast of Zeku County, Henan, the middle of Maqin and the north of Qumalai; 90.56% of the regional water retention showed a slight upward trend, with an average annual growth rate of 7.5 mm/year−1. The water production capacity of the study area is divided into five levels. Low-value areas are in the northwest of Qumalai County, Maduo County and Chengduo County, while high value areas are in Ruoergai, Aba and Hongyuan counties.
The total value of water conservation ecological services in the grassland ecosystem increased from USD 146.15×108 in 2001 to USD 227.59 × 108 in 2020, showing a steady growth trend. The areas with low value are in the improved grasslands, while the areas with high value are distributed in the alpine grassland. The ecological value of alpine grassland increased from USD 146.15 × 108 in 2001 to USD 227.59 × 108 in 2020, with an average annual growth rate of 80%.
The Chinese government has carried out ecological protection projects in the HYR. The implementation of the projects has achieved remarkable results, the grassland degradation has been alleviated and the water conservation function has been significantly improved. Qumalai, Chengduo and Maduo counties are in extremely important areas. It is proposed to classify them into priority areas for the ecological protection of water resources. In areas with high ecological values such as alpine meadows, temperate grasslands and mountain meadows in the headwaters, regional rotational grazing and grass balance should be introduced, and a comprehensive grazing ban should be implemented for lowland meadows, temperate deserts and alpine grasslands. Our suggestions are based on the policies and measures implemented in previous ecological projects, which are reasonable for the restoration and management of the water conservation function in the HYR.
This study evaluates the ecological function and value of water conservation in the HYR, provides reference data and policy suggestions for the formulation of ecological conservation measures and technical project investments and is of great significance for the sustainable development of water resources in the headwaters and watershed of the Yellow River.

Author Contributions

Conceptualization, H.L. and D.W.; methodology, D.W., Y.Z., J.Z. and W.Y.; formal analysis, D.W., X.Z., R.T. and H.L.; resources, T.H. and Y.Z.; data curation, D.W., R.T., T.H. and J.Z.; writing—original draft preparation, D.W., T.H., K.M., Y.P. and K.G.M.; writing—review and editing, D.W., K.M., W.Y., Y.P., H.L., K.G.M. and J.Z.; visualization, H.L., D.W. and W.Y.; supervision, H.L.; project administration, H.L.; funding acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (32171680 and 31772666).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brett, A.B. China’s response to a national land-system sustainability emergency. Nature 2018, 559, 193–204. [Google Scholar] [CrossRef]
  2. Immerzeel, W.; Wanders, N.; Lutz, A.; Shea, J.; Bierkens, M. Reconciling high-altitude precipitation in the upper Indus basin with glacier mass balances and runoff. Hydrol. Earth Syst. Sci. 2015, 19, 4673–4687. [Google Scholar] [CrossRef]
  3. Zhang, C.F.; Bao, Z.X.; Yang, X.T.; Liu, Z.X.; Guan, X.X.; Wang, G.-Q. Evolution characteristics and response relationships of hydro- meteorological variables in the source region of the yellow river during 1961 to 2013. J. N. China Univ. Water Resour. Electr. Power Nat. Sci. Ed. 2019, 40, 15–19. (In Chinese) [Google Scholar] [CrossRef]
  4. Hu, Y.; Maskey, S.; Uhlenbrook, S.; Zhao, H.L. Streamflow trends and climate linkages in the source region of the Yellow River, China. Hydrol. Process. 2011, 25, 3399–3411. [Google Scholar] [CrossRef]
  5. Lan, Y.C.; Zhao, G.H.; Zhang, Y.N.; Wen, J.; Liu, J.Q.; Hu, X.L. Response of runoff in the source region of the Yellow River to climate warming. Quat. Int. 2010, 226, 60–65. [Google Scholar] [CrossRef]
  6. Choi, C.Y.; Shi, X.; Shi, J.B.; Gan, X.J.; Wen, C.J.; Zhang, J.W.; Jackson, M.V.; Fuller, R.A.; Gibson, L. China’s ecological conservation redline policy is a new opportunity to meet post-2020 protected area targets. Conserv. Lett. 2021, 15, e12853. [Google Scholar] [CrossRef]
  7. Gao, J.X.; Zou, C.X.; Zhang, K.; Xu, M.J.; Wang, Y. The establishment of Chinese ecological conservation redline and insights into improving international protected areas. J. Environ. Manag. 2020, 264, 110505. [Google Scholar] [CrossRef]
  8. Shao, Q.Q.; Zhao, Z.P.; Liu, J.Y.; Fan, J.W. The characteristics of land cover and macroscopical ecology changes in the source region of three rivers on Qinghai Tibet Plateau during last 30 years. Geogr. Res. 2010, 29, 1439–1451. (In Chinese) [Google Scholar]
  9. Ragettli, S.; Pellicciotti, F.; Immerzeel, W.W.; Miles, E.S.; Petersen, L.; Heynen, M.; Shea, J.M.; Stumm, D.; Joshi, S.; Shrestha, A. Unraveling the hydrology of a Himalayan catchment through integration of high resolution in situ data and remote sensing with an advanced simulation model. Adv. Water Resour. 2015, 78, 94–111. [Google Scholar] [CrossRef]
  10. 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]
  11. Yang, D.; Liu, W.; Tang, L.; Chen, L.; Li, X.Z.; Xu, X.L. Estimation of water provision service for monsoon catchments of south China: Applicability of the InVEST model. Landsc. Urban Plan. 2019, 182, 133–143. [Google Scholar] [CrossRef]
  12. Kusre, B.C.; Baruah, D.C.; Bordoloi, P.K.; Patra, S.C. Assessment of hydropower potential using GIS and hydrological modeling technique in Kopili River Basin in Assam (India). Appl. Energy 2010, 87, 298–309. [Google Scholar] [CrossRef]
  13. Baker, T.J.; Miller, S.N. Using the soil and water assessment tool (SWAT) to assess land use impact on water resources in an east African watershed. J. Hydrol. 2013, 486, 100–111. [Google Scholar] [CrossRef]
  14. Villa, F.; Bagstad, K.J.; Voigt, B.; Johnson, G.W.; Portela, R.; Honzák, M.; Batker, D. A Methodology for adaptable and robust ecosystem services assessment. PLoS ONE 2014, 9, e91001. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, X.; Zhang, G.S.; Long, X.; Zhang, Q.; Liu, D.S.; Wu, H.J.; Li, S. Identifying the drivers of water yield ecosystem service: A case study in the Yangtze River Basin, China. Ecol. Indic. 2021, 132, 108304. [Google Scholar] [CrossRef]
  16. Chen, J.H.; Wang, D.C.; Li, G.D.; Sun, Z.C.; Wang, X.; Zhang, X.; Zhang, W. Spatial and temporal heterogeneity analysis of water conservation in Beijing-Tianjin-Hebei urban agglomeration based on the geodetector and spatial elastic coefficient trajectory models. GeoHealth 2012, 4, e2020GH000248. [Google Scholar] [CrossRef]
  17. Wei, P.; Chen, S.; Wu, M.; Deng, Y.; Xu, H.; Jia, Y.; Liu, F. Using the InVEST model to assess the impacts of climate and land use changes on water yield in the upstream regions of the Shule River Basin. Water 2021, 13, 1250. [Google Scholar] [CrossRef]
  18. 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 Users Guide; The Natural Capital Project: Stanford, CA, USA, 2014. [Google Scholar]
  19. 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]
  20. Leh, M.D.K.; Matlock, M.D.; Cummings, E.C.; Nalley, L.L. Quantifying and mapping multiple ecosystem services change in West Africa. Agric. Ecosyst. Environ. 2013, 221, 285. [Google Scholar] [CrossRef]
  21. Bai, Y.; Ochuodho, T.O.; Yang, J. Impact of land use and climate change on water-related ecosystem services in Kentucky, USA. Ecol. Indic. 2019, 102, 51–64. [Google Scholar] [CrossRef]
  22. Li, M.Y.; Liang, D.; Xia, J.; Song, J.X.; Ceng, D.D.; Wu, J.T.; Cao, Y.L.; Sun, H.T.; Li, Q. Evaluation of water conservation function of Danjiang river basin in Qinling Mountains, China based on InVEST model. J. Environ. Manag. 2021, 286, 112212. [Google Scholar] [CrossRef] [PubMed]
  23. Ullah, S.; You, Q.; Ali, A.; Ullah, W.; Jan, M.A.; Zhang, Y. Observed changes in maximum and minimum temperatures over China-Pakistan economic corridor during 1980–2016. Atmos. Res. 2019, 216, 37–51. [Google Scholar] [CrossRef]
  24. Mallick, J.; Talukdar, S.; Alsubih, M.; Salam, R.; Ahmed, M.; Kahla, N.B.; Shamimuzzaman, M. Analysing the trend of rainfall in Asir region of Saudi Arabia using the family of Mann-Kendall tests, innovative trend analysis, and detrended fluctuation analysis. Theor. Appl. Climatol. 2021, 143, 823–841. [Google Scholar] [CrossRef]
  25. Gumus, V.; Avsaroglu, Y.; Simsek, O. Streamflow trends in the Tigris river basin using Mann−Kendall and innovative trend analysis methods. J. Earth Syst. Sci. 2022, 131, 34. [Google Scholar] [CrossRef]
  26. Lepeška, T.; Wojkowski, J.; Wałęga, A.; Młyński, D.; Radecki-Pawlik, A.; Olah, B. Urbanization—Its hidden impact on water losses: Prądnik river basin, Lesser Poland. Water 2020, 12, 1958. [Google Scholar] [CrossRef]
  27. Irannezhad, M.; Ronkanen, A.K.; Klove, B. Effects of climate variability and change on snowpack hydrological processes in Finland. Cold Reg. Sci. Technol. 2015, 118, 14–29. [Google Scholar] [CrossRef]
  28. Ouyang, Z.Y.; Song, C.S.; Zheng, H.; Polasky, S.; Xiao, Y.; Bateman, I.J.; Liu, J.G.; Ruckelshaus, M.; Shi, F.Q.; Xiao, Y. Using Gross Ecosystem Product (GEP) to value nature in decision making. Proc. Natl. Acad. Sci. USA 2020, 117, 14593–14601. [Google Scholar] [CrossRef]
  29. Xie, G.D.; LU, C.X.; Leng, Y.F. Ecological assets valuation of the Tibetan Plateau. J. Nat. Resour. 2004, 18, 189–196. (In Chinese) [Google Scholar]
  30. Xie, G.D.; Zhang, C.X.; Zhen, L.; Zhang, L. Dynamic changes in the value of China’s ecosystem services—Sciencedirect. Ecosyst. Serv. 2017, 26, 146–154. [Google Scholar] [CrossRef]
  31. Costanza, R.; Groot, R.D.; Sutton, P.; Ploeg, S.V.D.; Anderson, S.J.; Kubiszewski, I.; Farber, S.; Turner, R.K. Changes in the global value of ecosystem services. Glob. Environ. Chang.-Hum. Policy Dimens. 2014, 26, 152–158. [Google Scholar] [CrossRef]
  32. Wang, J.L.; Lu, Y.H.; Zeng, Y.; Zhao, Z.J.; Zhang, L.W.; Fu, B.J. Spatial heterogeneous response of land use and landscape functions to ecological restoration: The case of the Chinese loess hilly region. Environ. Earth Sci. 2014, 72, 2683–2696. [Google Scholar] [CrossRef]
  33. Li, P.H.; Zhang, R.Q.; Xu, L.P. Three-dimensional ecological footprint based on ecosystem service value and their drivers: A case study of Urumqi. Ecol. Indic. 2021, 131, 108117. [Google Scholar] [CrossRef]
  34. Guo, S.L.; Hu, Z.Y.; Ma, H.Z.; Xu, D.D.; He, R.W. Spatial and temporal variations in the ecological efficiency and ecosystem service value of agricultural land in China. Agriculture 2022, 12, 803. [Google Scholar] [CrossRef]
  35. Chen, W.X.; Zeng, Y.Y.; Zeng, J. Spatio-temporal variation of ecosystem service values adjusted by vegetation cover: A case study of Wuyishan National Park Pilot, China. J. For. Res. 2022, 33, 851–863. [Google Scholar] [CrossRef]
  36. Rosswall, T. The international geosphere-biosphere programme: A study of global change (IGBP). Environ. Geol. Water Sci. 1992, 20, 77–78. [Google Scholar] [CrossRef]
  37. Friedl, M.A.; Gopal, S.; Muchoney, D.; Strahler, A.H. Global land cover mapping from MODIS: Algorithm design and preliminary results. In Proceedings of the 2000 IEEE International Geoscience and Remote Sensing Symposium, Honolulu, HI, USA, 24–28 July 2000. [Google Scholar] [CrossRef]
  38. Chen, D.; Zhuang, Q.F.; Zhu, L.; Zhang, W.J. Comparison of methods for reconstructing MODIS land surface temperature under cloudy conditions. Appl. Sci. 2022, 12, 6068. [Google Scholar] [CrossRef]
  39. Krishnan, M.V.N.; Sulemana, A. Regional scale analysis of land cover dynamics in Kerala over last two decades through MODIS data and statistical techniques. J. Environ. Stud. Sci. 2022, 5, 13412. [Google Scholar] [CrossRef]
  40. Li, W.J.; Hu, J.; Feng, T.; Liu, Q. Suitability and sensitivity of the potential distribution of cyclobalanopsis glauca forests under climate change conditions in Guizhou Province, Southwestern China. Atmosphere 2022, 13, 456. [Google Scholar] [CrossRef]
  41. Fatras, C.; Parrens, M.; Peña Luque, S.; Al Bitar, A. Hydrological dynamics of the Congo basin from water surfaces based on L-band microwave. Water Resour. Res. 2020, 18, e2020WR027259. [Google Scholar] [CrossRef]
  42. Lin, H.; Ma, H.; Nyandwi, C.; Feng, Q.; Liang, T. A new net primary productivity model and new management strategy of grassland classification based on CSCS in China. Rangel. J. 2021, 43, 285–295. [Google Scholar] [CrossRef]
  43. Sun, L.; Zhang, Y.J.; Zhang, T. Valuation of ecosystem services based on land use and land cover in source region of the Yellow River. Geogr. Geo-Inf. Sci. 2014, 30, 99–104. (In Chinese) [Google Scholar]
  44. Ge, Q.Q.; Xu, W.J.; Fu, M.C.; Han, Y.X.; An, G.Q.; Xu, Y.T. Ecosystem service values of gardens in the Yellow River Basin, China. J. Arid Land. 2022, 14, 284–296. [Google Scholar] [CrossRef]
  45. Tang, R.; Zhao, Y.; Lin, H. Spatio-Temporal Variation Characteristics of Aboveground Biomass in the Headwater of the Yellow River Based on Machine Learning. Remote Sens. 2021, 13, 3404. [Google Scholar] [CrossRef]
  46. Wang, Y.; Sun, Y.; Niu, F.; Wu, Q. Technology using 137Cs measurements to investigate the impact of soil erosion on soil nutrients in alpine meadows within the Yangtze River Region, China. Cold Reg. Sci. Technol. 2017, 135, 28–33. [Google Scholar] [CrossRef]
  47. Zhang, L.; Hickel, K.; Dawes, W.R.; Chiew, F.; Western, A.W.; Briggs, P.R. A rational function approach for estimating mean annual evapotranspiration. Water Resour. Res. 2004, 40, 502. [Google Scholar] [CrossRef]
  48. 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]
  49. Shi, X.L.; Du, C.L.; Guo, X.D.; Shi, W.J. Heterogeneity of water-retention capacity of forest and its influencing factors based on meta-analysis in the Beijing-Tianjin-Hebei region. J. Geogr. Sci. 2021, 31, 69–90. [Google Scholar] [CrossRef]
  50. Zhang, T.T.; Lyu, X.G.; Zou, Y.C.; Liu, J.P.; Jiang, M.; Xu, C.G.; Zhou, C.C.; Xu, C.; Xue, Z.S. Value assessment of wetland ecosystem services in the Da Hinggan Mountains, China. Chin. Geogr. Sci. 2022, 32, 302–311. [Google Scholar] [CrossRef]
  51. Bai, Y.; Jiang, B.; Wang, M.; Li, H.; Alatalo, J.M.; Huang, S.F. New Ecological Redline Policy (ERP) to secure ecosystem services in China. Land Use Policy 2016, 55, 348–351. [Google Scholar] [CrossRef]
  52. Bai, Y.; Wong, C.P.; Jiang, B.; Hughes, A.C.; Wang, M.; Wang, Q. Developing china’s ecological redline policy using ecosystem services assessments for land use planning. Nat. Commun. 2018, 9, 3034. [Google Scholar] [CrossRef]
  53. Gao, J.X.; Wang, Y.; Zou, C.X.; Xu, D.L.; Lin, N.F.; Wang, L.X.; Zhang, K. China’s ecological conservation redline: A solution for future nature conservation. Ambio 2020, 49, 1519–1529. [Google Scholar] [CrossRef] [PubMed]
  54. Zhang, K.; Zou, C.; Lin, N.; Qiu, J.; Pei, W.; Yang, Y.; Bao, X.; Zhang, Z. The Ecological conservation redline program: A new model for improving China’s protected area network. Environ. Sci. Policy 2022, 131, 10–13. [Google Scholar] [CrossRef]
  55. Lv, L.T.; Ren, T.T.; Sun, C.Z.; Zheng, D.F.; Wang, H. Spatial and temporal changes of water supply and water conservation function in Sanjiangyuan National Park from 1980 to 2016. Acta Ecol. Sin. 2020, 3, 993–1003. (In Chinese) [Google Scholar] [CrossRef]
  56. Yin, Y.H.; Wu, S.H.; Zhao, D.H.; Dai, E.H. Ecosystem water conservation changes in response to climate change in the Source Region of the Yellow River from 1981 to 2010. Geogr. Res. 2016, 1, 49–57. (In Chinese) [Google Scholar]
  57. Xinhua News Agency. The central committee of the communist party of China and the state council issued the outline of the yellow river basin ecological protection and high-quality development plan. China Water Resour. 2021, 21, 3–16. (In Chinese) [Google Scholar]
  58. Xu, X.L.; Wang, L.; Li, J.; Cai, H.Y. Analysis of the grassland restoration trend and degradation situation in the “three-river headwaters” region since the implementation of the ecological project. Acta Geogr. Sin. 2017, 19, 50–58. (In Chinese) [Google Scholar]
  59. Yang, J. Study on Grassland Ecosystem Service and Its Trade-Off and Synergy in the Yellow River Basin. Master’s Thesis, Gansu Agricultural University, Lanzhou, China, 21 October 2021. [Google Scholar]
  60. Zhao, Y.T. Alpine Grassland Ecological Security Assessment for the Headwaters of the Yellow River, China. Master’s Thesis, Lanzhou University, Lanzhou, China, 3 June 2021. [Google Scholar]
  61. Kang, H.H. Evaluation and Prediction of Ecosystem Services Based on Land Cover in the Source Region of Yellow River. Master’s Thesis, Gansu Agricultural University, Lanzhou, China, 3 June 2017. [Google Scholar]
  62. Yellow River Water Resources Bulletin. Available online: http://www.yrcc.gov.cn/other/hhgb/ (accessed on 24 May 2022).
  63. Zhu, H.X.; Fu, B.J.; Wang, S.; Zhu, L.H.; Zhang, L.W.; Jiao, L.; Wang, C. Reducing soil erosion by improving community functional diversity in semi-arid grasslands. J. Appl. Ecol. 2015, 52, 1063–1072. [Google Scholar] [CrossRef]
  64. Zhang, M.N.; Liu, B.; Li, G.Y.; Kuang, Y.Y.; Yue, X.Q.; Jiang, S.C.; Liu, J.S.; Wang, L. The relative and combined effects of herbivore assemblage and soil nitrogen on plant diversity. Sci. China-Life Sci. 2021, 65, 830–837. [Google Scholar] [CrossRef]
  65. The Grassland Law of the People’s Republic of China. Available online: http://www.forestry.gov.cn/main/3949/20180918/114120127762082 (accessed on 1 June 2022).
  66. Huang, L.; Ning, J.; Zhu, P.; Zheng, Y.H.; Zhai, J. The conservation patterns of grassland ecosystem in response to the forage-livestock balance in North China. J. Geogr. Sci. 2021, 31, 518–534. [Google Scholar] [CrossRef]
  67. Wang, C.X.; Liang, W.; Yan, J.W.; Jin, Z.; Zhang, W.B.; Li, X.F. Effects of vegetation restoration on local microclimate on the Loess Plateau. J. Geogr. Sci. 2022, 32, 291–316. [Google Scholar] [CrossRef]
Figure 1. Location of the HRY, grassland types and county locations.
Figure 1. Location of the HRY, grassland types and county locations.
Water 14 02553 g001
Figure 2. Spatial distribution of water retention mean value in the HYR.
Figure 2. Spatial distribution of water retention mean value in the HYR.
Water 14 02553 g002
Figure 3. The trend of water retention Change in the HYR from 2001 to 2020.
Figure 3. The trend of water retention Change in the HYR from 2001 to 2020.
Water 14 02553 g003
Figure 4. The trend of water retention mean value change in the HYR from 2001 to 2020. The Y axis represents the statistical sequence of water yield calculated at time X, which includes time series UF and inverse time series UB, and x refers to the years 2001–2020.
Figure 4. The trend of water retention mean value change in the HYR from 2001 to 2020. The Y axis represents the statistical sequence of water yield calculated at time X, which includes time series UF and inverse time series UB, and x refers to the years 2001–2020.
Water 14 02553 g004
Figure 5. Water conservation function values of different grassland types (USD × 108).
Figure 5. Water conservation function values of different grassland types (USD × 108).
Water 14 02553 g005
Figure 6. Spatial division of water retention in the HYR (the pie chart shows the proportion of water retention distribution ranges).
Figure 6. Spatial division of water retention in the HYR (the pie chart shows the proportion of water retention distribution ranges).
Water 14 02553 g006
Figure 7. The importance of water conservation evaluation.
Figure 7. The importance of water conservation evaluation.
Water 14 02553 g007
Figure 8. InVEST modeled mean annual water yield against observed mean annual water yield. Blue line indicates a relationship with intercept = zero and slope = one.
Figure 8. InVEST modeled mean annual water yield against observed mean annual water yield. Blue line indicates a relationship with intercept = zero and slope = one.
Water 14 02553 g008
Table 1. The environment variables and data sources.
Table 1. The environment variables and data sources.
CategoryAbbreviationVariable InterpretationResolution/UnitData source
SCLAY 1Clay content (0–30 cm)250 mHWSD (https://www.fao.org/soils-portal)
The accessed date: 11 March 2021
CLAY 2Clay content (30–100 cm)250 mHWSD (https://www.fao.org/soils-portal)
The accessed date: 11 March 2021
SAND 1Sand content (0–30 cm)250 mHWSD (https://www.fao.org/soils-portal)
The accessed date: 11 March 2021
SAND 2Sand content (30–100 cm)250 mHWSD (https://www.fao.org/soils-portal)
The accessed date: 11 March 2021
SILTSilt content (0–30 cm)250 mHWSD (https://www.fao.org/soils-portal)
The accessed date: 11 March 2021
SOCSoil organic carbon (0–30 cm)250 mHWSD (https://www.fao.org/soils-portal)
The accessed date: 11 March 2021
LUCCLand use/land cover change30 m2005, 2010, 2015, 2020 four years’ data of Soil Science Database in China (http://sdb.casnw.net/portal/)
The accessed date: 20 March 2021
IGBPLand use classification data250 mMCD12Q1.V006(https://climatedataguide.ucar.edu/climate-data/ceres-igbp-land-classification)
The accessed date: 30 June 2020
SDSoil depth1000 mSoil Science Database in China (http://sdb.casnw.net/portal/)
The accessed date: 2 April 2021
VelocityVelocity coefficient--USDA-NRCS
The accessed date: 5 April 2021
APREAnnual precipitation1000 mMeteorological Station
The accessed date: 30 June 2020
TEMAnnual temperature1000 mMeteorological Station
The accessed date: 30 June 2020
VEvapotranspiration1000 mMeteorological Station
The accessed date: 30 June 2020
RARadiation from the top of solar atmosphere1000 mChina Meteorological Data Service Centre (http://data.cma.cn/)
The accessed date: 30 June 2020
TLONLongitude30 mArcGis Calculated
LATLatitude30 mArcGis Calculated
DEMElevation30 mGeospatial Data Cloud(http://www.gscloud.cn/)
The accessed date: 30 June 2020
SLPSlope30 mArcGis Calculated
TITopographic index30 mArcGis Calculated
BPAWCPlant available water content30 mArcGis Calculated
KsatVegetation evapotranspiration coefficient--FAO56
NPPNet Primary Productivity1000 mMOD16A3
ECropWheat growing areas; the per unit yield; unit-pricehm2; kg/hm2; USD/kgChina Statistical Yearbook 2002–2021, China Grain and Material Reserve Yearbook 2002–2021
The accessed date: 16 February 2022
GDPCounty GDPUSDQinghai Statistical Yearbook 2021; Gansu Statistical Yearbook 2021; Aba Prefecture Statistical Yearbook 2021
The accessed date: 17 February 2022
Table 2. Water retention change trend assessment table.
Table 2. Water retention change trend assessment table.
SENslopeZTrend
>0.001>1.96Significantly increasing
>0.001−1.96–1.96Increasing
−0.001–0.001−1.96–1.96Stable
<0.001−1.96–1.96Decreasing
<0.001<−1.96Significantly decreasing
Table 3. Statistical description of water yield of counties in the HYR (unit: mm).
Table 3. Statistical description of water yield of counties in the HYR (unit: mm).
Country2001–20052005–20102011–20152016–2020
MeanMeanMeanMean
Xinghai290.24339.6299.4424.32
Qumalai195.27240.17220.2263.12
Maduo262.43303.63303365.73
Tongde300.93349.41315.91447.75
Zeku362.96414.2381.63509.36
Maqin395.5432.09419.02516.88
Chengduo277.36327.05331.08356.71
Henan412.14446.51422.37540.13
Gande458.34483.5490.15569.31
Maqu505.19496.08529.65632.68
Dari451.44473.6509.8558.69
Ruoergai565.87530.2607.6710.77
Jiuzhi545.21548.92593.75663.29
Aba554.83538.09623.99710
Hongyuan603.22587.23686.87785.74
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, D.; Zhao, Y.; Yang, W.; Ma, K.; Hao, T.; Zhao, J.; Tang, R.; Pu, Y.; Zhang, X.; Mujtaba, K.G.; et al. Ecological–Economic Assessment and Managerial Significance of Water Conservation in the Headwaters of the Yellow River. Water 2022, 14, 2553. https://doi.org/10.3390/w14162553

AMA Style

Wang D, Zhao Y, Yang W, Ma K, Hao T, Zhao J, Tang R, Pu Y, Zhang X, Mujtaba KG, et al. Ecological–Economic Assessment and Managerial Significance of Water Conservation in the Headwaters of the Yellow River. Water. 2022; 14(16):2553. https://doi.org/10.3390/w14162553

Chicago/Turabian Style

Wang, Danni, Yuting Zhao, Wenxue Yang, Kexin Ma, Tianxing Hao, Jingwei Zhao, Rong Tang, Yanfei Pu, Xiujuan Zhang, Kalhoro Ghulam Mujtaba, and et al. 2022. "Ecological–Economic Assessment and Managerial Significance of Water Conservation in the Headwaters of the Yellow River" Water 14, no. 16: 2553. https://doi.org/10.3390/w14162553

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