1. Introduction
Soil conservation is a vital regulating service in regional ecosystems through ecosystem structure and processes [
1] that is defined as the ability of an ecosystem to prevent soil erosion [
2,
3]. Soil erosion is an environmental problem that exists worldwide that can destroy soil resources and lead to land deterioration and serious ecological degradation [
4,
5,
6]. The Status of the World’s Soil Resources report shows that most of the world’s soil resources are in fair or poor condition [
4], with more than 2 billion hectares affected by soil erosion, costing the global GDP about USD 8 billion per year [
7]. In China, the threat of soil erosion is equally serious, especially in the Yellow River Basin (YRB). In 2021, the total area of soil and water loss in the YRB reached about 2.6 × 10
5 km
2 [
8]. Soil conservation is the prevention of soil loss through erosion or the loss of fertility induced by overuse, acidification, salinization, or other chemical soil pollution [
9,
10]. Improving soil conservation functions is important for the ecological protection and high-quality development of the YRB [
11].
According to the Soil and Water Conservation Monitoring Center of the Ministry of Water Resources of China, which issued the
Technical Guidelines for Monitoring the Dynamics of Soil and Water Loss in 2020, of the 76 tributaries in the YRB with watersheds larger than 1000 km
2, there are 11 tributaries with serious soil erosion. Among them, the most serious and difficult to manage are the Ten Kongduis (kongdui is the translation of “short-term flood gullies” in Mongolian), which are ten tributaries in the upper Inner Mongolia section of the YRB. Pisha sandstone, known as “Earth’s ecological cancer”, is located in the upper stream of the Ten Kongduis, which are highly susceptible to soil and water loss. The average annual sediment delivered from the Ten Kongduis to the YRB is about 2.7 × 10
7 t, which is more than one-tenth of the total amount of sediment entering the YRB in China. The Ten Kongduis are one of the direct producers of the “hanging river on the ground” of the Inner Mongolia section and lower reaches of the YRB [
12]. In order to implement the strategy of ecological protection and high-quality development of the YRB, the Chinese Inner Mongolia Autonomous Region has listed the comprehensive management of the Ten Kongduis as a major initiative of the development strategy [
13]. Knowledge of sediment initiation and transport processes allows management to intervene in the transport process by, for example, changing how land is used and managed to reduce sediment loads [
14]. Therefore, the study of the spatial and temporal variability in soil conservation in the Ten Kongduis is of extraordinary scientific significance both in terms of the discipline and for the ecological and environmental management of the region.
In general, field observations of soil erosion and sediment yield are most accurate. However, this is impractical for large study areas because of financial and human resource constraints. The literature shows that some models can be used to understand the processes of soil erosion and sediment yield in watersheds [
15], including the Universal Soil Loss Equation (USLE), Revised Universal Soil Loss Equation (RUSLE), Water Erosion Predict Project (WEPP), European Soil Erosion Model (EUROSEM) [
16], Water and Tillage Erosion Model and Sediment Delivery Model (WaTEM/SEDEM) [
17], Soil and Water Assessment Tool (SWAT) [
18], Agricultural Non-Point Source Pollution (AGNPS) model [
19], and Limburg Soil Erosion Model (LISEM) [
20]. Meanwhile, because physically based soil loss and sediment transport models require large amounts of data, it is not easy to apply them in areas where data are scarce [
15]. So far, for the Ten Kongduis, there are only three hydrological stations. Therefore, the main study on the Ten Kongduis is based on the measured data at the three hydrological stations to analyze the changes in the runoff and sediment discharge and less on the whole region to carry out research on the spatial and temporal dynamics of the sediment yield. Therefore, the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) model, which works with few climate and soil data [
21], was chosen to evaluate the soil conservation in the Ten Kongduis at different spatial and temporal scales. The soil conservation module of the InVEST model considers the sum of soil erosion reduction and retention as soil conservation. The former is the difference between the potential and actual erosion of the soil, and the latter is the ability of the plot to intercept sediment uphill, which can be expressed as the product of the incoming sand volume and the sediment retention rate. Compared with the Universal Soil Loss Equation (USLE), this model integrates the erosion mitigation and sediment retention effects of vegetation, and the evaluation results are more scientific [
16]. Soil erosion and sediment yields have been simulated with the InVEST model in many countries around the world [
22,
23,
24,
25,
26,
27].
Thus, based on the InVEST model, the objective of this paper was to quantify and map the characteristics of the spatial and temporal variations in soil conservation services in the Ten Kongduis since 2000, obtain the effects of rainfall and land use on soil conservation, assess the spatial and temporal distribution of sediment retention services, and identify where sediment retention occurs in the landscape. The results can allow water and land managers to investigate the spatial and temporal changes in ecosystem soil and water conservation functions and their driving factors in typical management projects carrying out ditching and land reclamation and design better strategies to reduce sediment loads by changing the land use and management practices in the Ten Kongduis and YRB through understanding where sediments are generated and transported and provide a scientific basis for soil and water conservation and the ecosystem safety management of watersheds.
4. Discussions
It is worth mentioning that a common feature of these changes is that when the value of rainfall erosivity in the base year was greater than in the simulation year and the land use data were changed, the impact of land use on avoided export was negative, and when the value of rainfall erosivity in the base year was less than in the simulation year and the land use data were changed, the impact of land use on avoided export was positive. The impact of land use on avoided export was positive. Meanwhile, the effect of land use on avoided erosion was exactly the opposite of that of avoided export. The preliminary inference is that because barren land showed significant changes between years, while other land uses did not change significantly between years, when the rainfall erosivity was large, the high reduction in vegetation at this time eroded soil from the slopes, and the amount of soil erosion was much larger than the amount of sediment export. When the rainfall erosivity was small, the vegetation at this time mainly reduced the eroded soil entering the stream, and the difference between the soil erosion and sediment export decreased.
Based on the above analysis, it can be seen that the rain erosivity factor had a positive effect on soil conservation values, both for the years with the largest and smallest avoided erosion and avoided export values. The rain erosivity factor is a reflection of the ability of rainfall to strip, move, and scour the surface of the soil, showing the potential ability of rainfall to cause soil erosion and loss [
6]; thus, it seems that the greater the rain erosivity factor, the greater the amount of soil erosion. Meanwhile, according to the definition of soil conservation, it can be seen that avoided erosion and avoided export are the avoidance of erosion and reduction in erosion, so they can be used to identify where sediment is trapped/retained in the landscape. In this way, it seems that the effect of rain erosivity on both should be negative. However, from the formulas for avoided erosion and avoided export, the two are sized by the difference between the RKLS and USLE. It appears that both the RKLS and USLE are directly proportional to the rainfall erosivity, i.e., the greater the rainfall erosivity, the greater the values of both the RKLS and USLE, when other factors are held constant. Therefore, it can be inferred that for the Ten Kongduis, the extent of the increase in the RKLS is greater than the extent of the increase in the USLE when the rainfall erosivity increases, and, at this time, the benefit of vegetation and good management practices is more pronounced. Where more erosion is created, more erosion is retained. Soil that has been eroded away from slopes under vegetation or other water conservation measures may not necessarily be transported to the stream channel in the current year; thus, focusing solely on protecting high-retention areas does not necessarily address where erosion occurs in the first place.