Next Article in Journal
The Influence of the South-to-North Water-Diversion Project on Terrestrial Water-Storage Changes in Hebei Province
Previous Article in Journal
Particle Size Distribution and Composition of Soil Sample Analysis in a Single Pumping Well Using a Scanning Electron Microscope Coupled with an Energy Dispersive X-ray (SEM-EDX) and the Laser Diffraction Method (LDM)
Previous Article in Special Issue
Research on the Features of Rainfall Regime and Its Influence on Surface Runoff and Soil Erosion in the Small Watershed, the Lower Yellow River
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of the Plateau Pika Mound Numbers on Soil Water Erosion Properties in Alpine Meadows of the Yellow River Source Zone, Western China

1
Geological Engineering Department, Qinghai University, Xining 810016, China
2
Key Lab of Cenozoic Resource and Environment in North Margin of the Tibetan Plateau, Xining 810016, China
3
State Key Laboratory of Plateau Ecology and Agriculture, Qinghai University, Xining 810016, China
4
College of Agriculture and Animal Husbandry, Qinghai University, Xining 810016, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(17), 3111; https://doi.org/10.3390/w15173111
Submission received: 4 July 2023 / Revised: 14 August 2023 / Accepted: 28 August 2023 / Published: 30 August 2023
(This article belongs to the Special Issue Soil Erosion Monitoring and Modeling)

Abstract

:
The plateau pika (Ochotona curzoniae) actively contributes to soil erosion and meadow degradation in western China’s Yellow River source zone. This study aimed to elucidate the effects of the pika mound numbers on the hydrodynamic characteristics and soil water erosion through simulated rainfall experiments. The inhibition effects of restored vegetation growth on meadow degradation and soil erosion were explored using a revegetated pika mound as a control. The results showed that at a rainfall intensity of 30 mm/h, the soil loss per unit time increased and then decreased with rainfall time and that 15–20 min of rainfall duration was the sensitive period for soil loss in different pika mound patch lands. The degradation of meadows due to the activities of plateau pika is an essential factor influencing soil erosion, and the soil erosion rate is positively correlated with both the pika mound numbers and slope. The mean flow velocity can better describe the process of the soil erosion, and its value increased with the number of pika mounds and slope. The Reynolds number ranged from 57.85 to 153.63 (Re <500), and it was preliminarily determined that all slope runoff was laminar flow. The Froude number increased linear function with pika mound numbers (p < 0.01), and both the pika mound numbers and slope were significant factors affecting it (p < 0.05). The Darcy–Weisbach resistance coefficient instead decreased with the pika mound numbers and slope, and the inhibitory effect of vegetation on soil erosion was probably limited when the mound numbers reached a certain level. According to the results of the grey correlation and Pearson correlation analysis, changes in the number of pika mounds led to variability in the soil water erosion properties by altering the landscape scale effect. The number of pika mound patches (NP), edge length index (TE), area index (AREA), and volume index (V) were the key influencing factors on the soil erosion rate and hydrodynamic parameters. We conclude that plateau pika activities intensify meadow degradation, and the continuous increase in pika mounds decreases the vegetation cover and flow resistance and increases soil erodibility. Controlling the intensity of plateau pika activities will effectively prevent soil erosion in degraded areas.

1. Introduction

The Yellow River source zone is in the northeast of the Qinghai–Tibet Plateau, the major flow-producing area of the upper Yellow River and an essential water-conserving region in China [1,2]. Alpine meadows, as the most important cover type in the source zone of the Yellow River [3], account for about 80% of the total area, and their environmental changes directly affect the soil and water conservation and ecological security of the entire Yellow River Basin [4]. Recently, influenced by rodent activity, livestock production, and climate change, alpine meadows have been degraded to different degrees [5] and show a trend of increasing degradation yearly [6]. Plateau rodents are the local dominant biological species; the small rodents in the Yellow River source zone mainly include plateau zokor (Eospalax baileyi) and plateau pika (Ochotona curzoniae) [7]. At present, along with the phenomenon of meadow vegetation degradation and soil degradation, which also provides a good vision and site for the activities of plateau rodents, the local government’s early rodent extermination behavior has seriously affected the number of their predators, resulting in a continuous surge in their populations and an expansion of their activity range in recent years. In addition, although the activities of these rodents may have some positive effects on soil air circulation, nitrogen cycling, and vegetation community restoration, their burrowing and gnawing behaviors disrupt the natural vegetation community features and soil structures of alpine meadows [8,9,10], and the accumulated loose mounds indirectly form secondary bald patches of land and supply a material source for soil erosion [11,12]. So far, many scholars have concluded that rodent activity-induced meadow degradation areas stimulate soil erosion under various wind, water, or freeze–thaw dynamics and generate or exacerbate other natural hazards [13,14], negatively affecting the sustainability and security of the whole alpine meadow ecosystem [4,15].
Landscape patterns are the consequence of coupling diverse ecological processes [16], and studying the spatial pattern changes in landscape patches is significant to understand the soil degradation and vegetation degradation mechanisms [17,18]. The source zone of the Yellow River is influenced by rodent activities [19], and the landscape fragmentation is severe, among which the rodent mound patches formed by plateau pika activities lead to the reverse succession of the meadow ecosystem, which has a specific impact on community habitat texture heterogeneity, species diversity, and soil erodibility [9,20,21]. Landscape indices are quantitative indicators that reflect the structural composition and spatial configuration of the landscape, which serve as a bridge between landscape patterns and ecological processes [22], and currently, with the advancement of 3S technology and landscape ecology, studies on the effects of soil erosion based on landscape indices have been widely carried out [23,24,25]. For example, Zhao [26] et al. showed in riparian zone areas that vegetation fragmentation distribution increased soil erodibility; the landscape pattern influence on the soil erosion process depended on the landscape scale. Sun et al. [27] used the Markov model to associate the relationship between landscape indices and soil sand production and suggested that the landscape index could be an essential indicator to influence the soil erosion sand production process. However, due to the soil erosion process’s complexity and the landscape index’s limitations [28,29], selecting an appropriate landscape index will effectively advance the intrinsic association between changes in landscape patterns and soil erosion processes [30]. Furthermore, hydrodynamic characteristics are essential to reflect the soil erosion process on slopes and its kinetic mechanism [31]. However, there are still significant gaps regarding the response mechanisms of the slope surface flow hydrodynamic characteristics to landscape patterns, and few studies have reported the role of landscape patterns on hydrodynamic characteristics. Therefore, the study of the number, shape, volume, and spatial distribution pattern of plateau pika mound patches on soil erosion and hydrodynamic characteristics based on microhabitat scale is vital to investigate in depth the response mechanism of meadow patches to hydraulic soil erosion.
Soil erosion is the process by which surface soil and its parent material are separated, transported by a combination of internal and external forces, closely related to rainfall, topography, soil, and vegetation [32,33,34,35]. Rainfall is a major driver of soil erosion, and many scholars have shown that changes in the rainfall-runoff can directly affect the soil erosion conditions and hydrodynamic characteristics [26,36]; hydrodynamic characteristics serve as visual reflections of soil erosion processes, which are related to landscape patterns, slope, and vegetation cover [37,38]. For example, Chen et al. [39] suggested that landscape fragmentation can increase soil erosion by altering the hydrological connectivity and sediment transport pathways. Li et al. [40] showed that as the vegetation cover decreased from 75% to 25%, the runoff velocity increased, the runoff resistance decreased, and the soil anti-erosion capacity decreased. Sun et al. [41] found that the flow resistance decreased when the slope and precipitation intensity increased, and that the Reynolds number can accurately describe the process of soil erosion. Furthermore, in the context of the continuous deterioration of the ecological environment of alpine meadows, so far, although many scholars have made significant progress in research on soil erosion and hydrodynamic characteristics, few studies have been conducted to quantify the characteristics of soil water erosion of the pika mound patch lands in the Yellow River source zone.
This research attempts to fill the above knowledge gap by conducting simulated rainfall experiments in areas with frequent activities of plateau pika with the aim to clarify the influence of pika mound numbers on the soil erosion and hydrodynamic characteristics. The targets were to (1) clarify the soil erosion regimes under different pika mound numbers; (2) reveal the hydrodynamic characteristics of surface flow under different pika mound numbers; (3) elucidate the influence of pika mound numbers on the soil erosion and hydrodynamic characteristics; (4) explore the inhibitory effect of restored vegetation growth on meadow degradation and soil erosion. The study results help reveal the soil erosion mechanism in the degraded meadow regions by rodent activities and provide a theoretical basis for restoring and managing meadow ecological vegetation.

2. Materials and Methods

2.1. Study Area

This study area is situated in Henan County, Huangnan Prefecture, Qinghai Province (34°47′ N, 101°32′ E), which is part of the alpine meadow region in the source zone of the Yellow River (Figure 1a). The altitude is about 3600 m, the average annual wind speed is 2.6 m s−1, and the annual precipitation is 597.1~615.5 mm, which is a typical continental plateau type climate. Rains are mainly centered in June–September of each year, and heavy rainfall events occur primarily in July–August, with rainfall durations ranging from 30 min to 180 min, characterized by short duration, suddenness, and small occurrence range [42,43]. The type of the meadow is mainly alpine dwarf tarragon meadow; the major plant species are Kobresia capillifolia, Kobresia humilis, Kobresia pygmaea, and Elymus nutans. The type of soil is humic cambisols (WRB), in which the undegraded areas are dominated by chalk and clay grains, while the degraded areas are dominated by sand grains and show significant coarse-graining and sanding. In addition, rodent activity is rampant in the study area, and many mounds have accumulated on the meadow surface (Figure 1b), among which plateau pika mounds are hemispherical and relatively gentle, with an average height of about 8 cm and a diameter of about 45 cm. The density of pika mounds is 0.36 mounds m−2 (Figure 1c). The soil properties of the pika mounds are highly susceptible to more severe soil erosion, and eventually intensify meadow degradation due to the apparent differences that occur before and after (Table 1), having negative impacts on the regional ecological environment and development of livestock [44].

2.2. Experimental Design

This study mainly focused on the erosive effects of rainfall scouring on the different pika mound patch lands. Through the field survey, we found that the plateau pika is widely distributed in the study area, and there is no apparent requirement for slope in its habitat choice, resulting in some of the distribution areas of the pika mounds having slopes reaching 33–35°, and the percentage of the area with slopes exceeding 30° was about 9.24%. Therefore, we chose to set up three removable water erosion test chambers in the field that were 1.5 m × 3 m in area and 0.5 m in depth, while the slope was designed as 10°, 20°, and 30°, respectively. The bottom of the box had natural drainage holes, and the front end was equipped with a “V”-shaped converging sediment deflector hole to collect and test the sediment runoff volume (Figure 2a). The soil in the test box was filled with native meadow soil from the study area that had been closed for more than 4 years and contained a vegetative root–soil structure. In situ meadow squares of 50 cm × 50 cm × 50 cm were cut from healthy meadow areas in June 2020. These squares were then sequentially pieced together and laid into the corresponding water erosion test chamber. The edge joints between the meadow squares were artificially filled and repaired to produce a complete test site. All test boxes were left in the field for 1–2 years after completion of the above preparatory work, and soil water erosion tests were carried out after the seams between the meadow squares were fully healed and the vegetation cover was greater than 90% (Figure 2b).
The mounds in our study were fresh and uncrusted mounds formed within a short period after the soil was scattered and buried by vegetation death due to the mounding behavior of the plateau pika. The number of mounds in the test box was controlled to be 1–5, and the number of mounds was set to 1, 2, 3, 4, and 5 by artificially simulating mounds of the same shape, size, weight, and water content in the test boxes. The soil of the simulated mounds was selected from fresh mound soil accumulated by highland pika activities in the study area, and the more disturbed mounds were artificially restored twice. Since the initial vegetation cover in each test chamber was unequal, the meadow vegetation was artificially treated with equal cover (i.e., the vegetation cover was all reduced to about 70%). Immediately after completing the three sets of replicate tests, the number of mounds was increased to two by artificially simulating additional mound patches and so on, until a whole set of hydraulic erosion tests were completed one after another. In addition, we found during the field survey that the recovery of the dominant plant species always accompanied the succession process of the pika mounds. This study conducted a control test on the old pika mound with about 60% of the surface vegetation restored, which was comparable in terms of vegetation growth, external morphology, and basic soil properties.
The hydraulic erosion test used a field portable artificial rainfall simulator (Figure 3a). The device’s working principle is as follows: when connecting the rainfall device, extend the inlet pipe of the pump into the water tank, set up a layer of filtering net at the water intake to avoid sediment impurities into the pump, and finally, start the pump so that the water extracted by the pump is evenly sprinkled on the ground through the nozzle. The vertical distance of the rainfall nozzle from the center of the meadow surface in the test chamber is 3 m, which distributes rainfall evenly from the center of the 1 mm diameter hole. The effective radius of low rainfall is 1.5 m, and the rainfall intensity controller can regulate the precipitation intensity between 0 and 60 mm h−1 with a simulated degree of ±0.5 mm h−1. Based on the real-time monitoring data of intense rainfall weather in 2018–2021 from small meteorological stations in the test area and the test demand, the average rainfall intensity for all trials was set to 30 mm h−1. In addition, to prevent the wind from affecting the trajectory and direction of the raindrops, the testing area was enclosed by a tarpaulin during the tests (Figure 3b).
This experiment recorded the time when runoff appeared from the inflow hole, the total precipitation duration was set to 1 h, and the runoff sediment samples were taken at 5-min intervals using catchment buckets. Meanwhile, soil erosion due to rain splash was included as a component of the measured sediment flux, considering that soil erosion due to rain splash is not only related to soil redistribution but is also challenging to separate from the lost soil. All of the samples collected in the field were taken back to the laboratory, and the sediment separation process was carried out by sedimentation and filtration methods. Furthermore, the separated samples were dried in an 105 °C oven to measure the amount and the rate of soil loss. In addition, each slope’s water temperature, runoff width, and runoff velocity were measured every 5 min during the sustained rainfall events. The runoff flow rate was determined by the dye (KMnO4) traced method, and three measurement zones of 1 m, 2 m, and 3 m were set up from the top of the slope downward, with the length of the zones designed to be 0.8 m, and the mean flow rate of each zone was used as the real-time runoff flow rate of the slope.

2.3. Indicator Testing and Analysis

The different hydrodynamic parameters representing the runoff flow process on the slope were taken to accurately describe the effect of changing the number of pika mounds on the hydrodynamic characteristics. The hydrodynamic parameters selected for this study included the runoff depth (H), mean flow velocity (v), Froude number (Fr), Reynolds number (Re), and Darcy–Weisbach resistance coefficient (f); the specific formulae and their applicability are given in the citations [45,46].
The landscape pattern index is an essential means to quantitatively describe the landscape structure, function, and ecological process, which can highly concentrate the landscape pattern information including fundamental landscape elements such as patch number, shape, size, and distribution combination [47]. In this study, a total of six landscape pattern indices including the number index (NP), edge length index (TE), density index (PD), area index (AREA), volume index (V), and separation index (Ni) were selected at the patch, type, and landscape level to accurately describe the influence of the pika mound numbers on the soil water erosion and hydrodynamic characteristics. The calculation and ecological significance of these indices are listed in Table 2 [26,48].
This study used Microsoft Excel 2021 software for data reduction and drawing charts, etc. In addition, based on the average value of the hydrodynamic parameters and soil erosion rate, we used IBM SPSS Statistics 26.0 software to analyze the statistical and correlation between the landscape pattern of the plateau pika mound patches and soil erosion and hydrodynamic characteristics.

3. Results

3.1. Soil Erosion Regulation under Different Pika Mound Numbers

3.1.1. Comparison of Total Soil Loss under Different Pika Mound Numbers

During this rainfall event (60 min for example), the total soil loss per unit area increased with the number of pika mounds and slopes (Figure 4). At the 10° slope, the soil loss from degraded patch lands with one, two, three, four, and five mounds was 2.12, 3.58, 5.11, 6.98, and 8.71 times higher than that of the revegetated pika mound, respectively. The increase in the number of mounds had a more significant impact on soil loss, and the total soil loss per unit area adding in a multiplicative power function with the number of pika mounds (p < 0.01). Compared to the fresh mounds under these different quantities, the amount of soil erosion was significantly reduced because the vegetation protected and ameliorated the revegetated mounds, resulting in a fundamental improvement in their soil characteristics. In addition, when the slope was increased from 10° to 30°, the total soil loss per unit area of the revegetated mound and patch lands under one, two, three, four, and five mound numbers increased by 63.97%, 66.41%, 59.70%, 56.44%, 43.93%, and 40.01%, respectively. The statistical analysis results showed significant differences in the total soil loss between different numbers of mounds on the same slope (p < 0.05). Between the same number of mounds on different slopes, except for the revegetated mound, which did not show any significant difference, the total soil loss under the other mounds showed significant differences (p < 0.05). Furthermore, we also found that as the number of pika mounds increased, the average slope of the different test plots increased by 0.3 to 0.6°, indicating that both the slope and number of pika mounds were significant factors influencing soil erosion.

3.1.2. Soil Loss Characteristics under Different Pika Mound Numbers

In Figure 5, the soil loss under the three gradients followed an upward and then downward trend with the time of rainfall erosion and reached a peak at about 20 min. The change in soil loss showed a significant increase (p < 0.05) during the period of rainfall until 15–20 min, and the soil loss generated during this period contributed the most to the total soil loss, indicating that 15–20 min of rainfall erosion was a sensitive period for soil loss in these experiments. It can also be seen from the figure that the magnitude of soil loss variation in the pika mound patch lands at different slopes were significantly different. Compared to the revegetated pika mound, soil loss under the remaining number of mounds increased significantly with the number of mounds (p < 0.05). In addition, combined with the observation results, we found that the height of the mound settled with 1.36 cm h−1, the radius of the mound expanded with 2.45 cm h−1, and the shape of the mound changed from a hemispherical mound before rainfall to an oblate spherical mound after rainfall. The surface layer of the mound was gradually exposed to more sand and gravel during the rainfall process, but the overall content was relatively small. The enrichment of the sand and gravel phenomenon did not appear until the end of the rainfall event, where there were only a few rain traces and imbricated texture by rain splash. The sediment size changed significantly before and after rainfall, for example, on a 20° slope, the average particle size of the sediment particles were 0.31 mm in 20 min before rainfall and 0.18 mm in 20 min after rainfall. Therefore, we suggest that the soil loss status of degraded meadow areas due to plateau pika activities is related to the slope, rainfall duration, and the number of pika mounds. Under the continuous high-intensity rainfall conditions, the soil loss was serious when the slope was steeper and contained more mounds.

3.1.3. Soil Erosion Regimes under Different Pika Mound Numbers

Under the precipitation strength of 30 mm h−1, the soil erosion occurred in all test plots after 6.5 min of sustained rainfall, regardless of the slope and number of mounds. As shown in Figure 6, the soil erosion rates varied significantly between different numbers of pika mounds at the same slope. On a 10° slope, for example, the soil erosion rate varied from 0.61 to 4.68 g m−2 min−1 in the revegetated mound and the degraded patch lands with one mound; from 2.07 to 9.93 g m−2 min−1 in the degraded patch lands with two and three mounds; and from 4.87 to15.42 g m−2 min−1 in the degraded patch lands with four and five mounds. Furthermore, as the gradient reached 20° or 30°, the soil erosion rate increased significantly (p < 0.05), except for the revegetated mounds, and was positively correlated with the slope as a linear function. Overall, when the mounds were about the same size and weight, the soil erosion patterns were composed of two phases: a sharp rise and gentle decline, and the soil erosion rate peaked at about 25 min. It is also worth mentioning that an increase in mound numbers significantly enhanced the growth trend in the first phase and the fluctuation range in the second phase compared to the vegetation restoration of mouse mounds. This effect became more evident with the increasing number of pika mounds and slopes.

3.1.4. Relationship between the Pika Mound Numbers and Soil Erosion in a Degraded Meadow Area

By establishing a functional relationship between the pika mound numbers and vegetation cover with the soil erosion rate (Figure 7), we concluded that the soil erosion rate showed a linear positive correlation with the number of pika mounds (p < 0.01) while there was a linear negative correlation with vegetation coverage (p < 0.01). The results indicate that soil erosion in the pika mound patch lands was influenced by the number of pika mounds and related to the vegetation cover of the meadows destroyed by the activities of plateau pika.

3.2. Hydrodynamic Characteristics of Alpine Meadows under Different Pika Mound Numbers

3.2.1. Average Flow Rate

Flow velocity is a primary parameter to study the hydrodynamic characteristics of slope runoff, which is closely related to slope, vegetation cover, and runoff volume. As shown in Table 3, the mean flow velocity ranged between 2.29–4.26 cm s−1, 2.65–5.19 cm s−1, and 2.98–5.82 cm s−1 for different pika mound patch lands at 10°, 20°, and 30° slopes, respectively, and the mean flow velocity increased linearly with the number of mounds (Table 4). In addition, as the gradient rose from 10° to 30°, the mean flow velocity increased by 30.13%, 36.40%, 36.25%, 28.37%, 32.15%, and 36.62% for the revegetated mound and different mound patch lands, respectively. The statistical analysis showed that the mean flow velocities differed significantly (p < 0.05) for all slope and pika mound number combination patterns. We suggest that an increase in the number of mounds mainly increases the mean flow velocity by weakening the vegetation’s inhibitory and interception effect on surface runoff, and the mean flow velocity tends to be higher on steeper slopes due to inertia forces.

3.2.2. Flow Regime

The criterion value of the Reynolds number (Re = 500) is taken to determine whether the slope flow is laminar or turbulent. [27]. The experiment values of Re at 10°, 20°, and 30° slopes ranged between 57.85–93.15, 86.55–131.92, and 105.12–153.63, respectively (Table 3), and were positively correlated with the pika mound numbers as an exponential function (Table 4). Depending on the basis for the discrimination of open channel flows, we tentatively determined that all slope runoff was laminar. Furthermore, from the results of the statistical analysis, the variability in Re was not significant between different numbers of mounds at the same slope, but the difference in Re between different slopes at the same number of mounds was significant (p < 0.05), concluding that the effect of slope on Re was much more significant than the change in the number of pika mounds.
The Froude number (Fr) represents the flow state of the runoff. As shown in Table 3, the values of Fr were 0.143–0.322, 0.157–0.361, and 0.161–0.385 (none of them exceeded 1) for different slopes, so it was determined that the relationship of Fr with slope and the number of mounds was consistent with that of Re (i.e., it increased with the slope and the number of mounds), and Fr increased as a linear function of the pika mound numbers (Table 4). The analytical results showed that Fr had significant differences (p < 0.05) for different combinations of slope and number of mounds; both were significant factors affecting Fr (two-way ANOVA, p < 0.05).

3.2.3. Flow Resistance

The Darcy–Weisbach resistance coefficient (f) characterizes the magnitude of the resistance to runoff generated by the bedding surface as water flows along the slope. In Table 3, f reduced with an increasing number of mounds and slope, where the range of f was 34.89–159.07 for different pika mound patch lands at a 10° slope, and f decreased to 34.51–125.80 and 16.99–71.47 when the slope was increased to 20° and 30°, respectively. Statistical analysis showed that f showed significant differences (p < 0.05) between different numbers of mounds in the same slope excluding revegetated pika mounds, while f showed significant differences (p < 0.05) between different slopes for the same number of mounds. In addition, f was negatively correlated with the number of mounds by an exponential function (Table 4), indirectly reflecting that the steeper topography and the outbreak of pika mounds in the area are very likely to cause a significant reduction in flow resistance, thus increasing the degree of soil erosion.

3.3. Relationship between Landscape Pattern with Soil Water Erosion Properties under Different Pika Mound Numbers

3.3.1. Analysis of Landscape Pattern Indicators of Pika Mound Patches

We analyzed the spatial distribution pattern of the pika patch landscape by selecting the appropriate landscape index (Table 5). As the pika mounds increased, the TE, AREA, V, and Ni significantly differed (p < 0.05) at different slopes, while the PD, r, and h did not show significant differences. The analysis also showed that TE and Ni increased exponentially (p < 0.01) and AREA and V increased linearly (p < 0.01) as a function of the pika mound numbers, while PD, r, and h did not show significant functions with the pika mound numbers. In addition, we also found that NP, TE, AREA, and V composed the first principal component (PCA), so their interpretation was as high as 81.32%. There was a highly significant positive correlation between NP and TE, AREA, V, and a highly significant negative correlation with Ni, but no significant correlation with r and h.

3.3.2. Relationship between Landscape Pattern and Soil Erosion under Different Pika Mound Patches

The correlation analysis showed that NP, TE, AREA, and V were extremely significantly correlated with RSE at different slopes, with correlation values ranging from 0.883 to 0.958, while PD, Ni, r, and h were less correlated with the soil erosion rate, with correlation values ranging from 0.547 to 0.731 (Table 6). This table also shows a highly significant positive correlation between RSE and NP, TE, AREA, and V, a significant negative correlation with Ni, and no significant correlation with PD, r, and h at different slopes. The correlation analysis results concluded that the interaction effect between the landscape pattern of pika patches and soil erosion was evident, in which NP, TE, AREA, and V were important factors affecting the RSE, and NP could be used as the main control factor to determine the intensity of soil erosion. This result indirectly shows that the spatial configuration of landscape patterns can be effectively optimized by controlling the number of pika mounds and restoring a certain degree of vegetation growth in the plateau pika mounds outbreak area, suppressing soil erosion.
By establishing the regression equation between the landscape index of the pika mound patches and the soil erosion rate at different slopes (Figure 8), we concluded that the soil erosion rate was positively correlated with the edge length index (TE) and volume index (V) as a linear function (p < 0.01), positively correlated with the area index (AREA) as a power function (p < 0.01), and negatively correlated with the separation index (Ni) as a power function (p < 0.01), while there was no significant regular variation with the density (PD), radius (r), and height (h) of the pika mound patches.

3.3.3. Relationship between Landscape Pattern and Hydrodynamic Characteristics under Different Pika Mound Patches

Based on the results of the correlation analysis at different slopes (Table 7), there were highly significant positive correlations between NP, TE, AREA, and V with v and Fr, significant negative correlations with H and f, and significant positive correlations with Re. Ni showed a highly significant positive correlation with H and f, and a significant negative correlation with v, Re, and Fr. No significant correlations were found between PD, r, and h and the hydrodynamic parameters. It shows that there is indeed an interaction between the landscape pattern and hydrodynamic characteristics of pika mound patches, in which NP, TE, AREA, V, and Ni affect the connectivity and mobility of the slope surface flow by changing the landscape scale effect, which leads to significant variability in the hydrodynamic characteristics, and they can all be significant factors affecting the hydrodynamic characteristics of slope runoff.

3.3.4. Relationship between Soil Erosion and Hydrodynamic Characteristics under the Different Pika Mound Patches

Based on the results of the correlation analysis (Table 8), v, Re, and Fr were most significantly associated with the RSE at different slopes, with correlation values around 0.8; whereas H and f were relatively less associated with the RSE, and the correlation values ranged from 0.569 to 0.690. In addition, a highly significant positive correlation was found between v, Re, and Fr with the RSE, while there was a significant negative correlation between H and f with the RSE. We initially judged that the average flow rate could be better to describe the soil erosion process.

4. Discussion

4.1. Comparison with Other Similar Studies

The power source of soil water erosion is mainly rainfall and slope runoff, and the process is influenced by many factors such as the vegetation cover, slope, and soil properties [11,32,33,49]. In order to further characterize the soil water erosion of the pika mounds, the results of the present study were compared with other similar reports. We found that the amount of soil loss per unit time produced showed an increase and then a decrease with rainfall duration (Figure 5), consistent with the findings of Li et al. [12] on plateau zokor mounds. However, from the information on localized changes in soil loss, the peak points of soil loss in the two studies showed apparent differences, which were considered likely to be related to the rainfall intensity and number of rodent mounds. In addition, our research concluded that the soil erosion rate of the vegetation restoration pika mound was reduced by 1.03 to 6.74 times compared with the different numbers of mounds. However, the soil erosion rate increased with the number of mounds and slope (Figure 6). The analyzed reason mainly lies in the fact that the increase in the number of mounds also reduces the vegetation cover, weakening the inhibition and interception function of meadow vegetation [50] and ultimately exacerbating soil erosion under the superimposed effect with slope. Our results agreed with the results of Zhao et al. [26] and Chen et al. [39] in the degraded vegetation area, and we believe that the vegetation restoration growth in rodent-infested degraded areas can effectively prevent soil erosion. At the same time, it is essential to reasonably control the activity range of plateau pika and the density of pika mounds. Marzen et al. [51] suggested that the influence of rainfall on slope erosion is largely the result of the potential erosive force of rain splash and runoff on loosening topsoil. We also observed that the rain splash would lead to the rain traces and stacked tile texture on the surface of the pika mounds. Rain splash is also an essential element in the study of soil erosion mechanisms, so studying the effect of rain splash on the soil erosion of pika mounds with different rainfall intensities is also a scientific problem that needs to be investigated. Hydrodynamic characteristics are an essential area to reflect the soil erosion process; Zhang et al. [37] and Yuan et al. [52] found that the increase in flow velocity contributed to the runoff connectivity and soil erodibility, which led to an increase in the surface runoff capacity and soil erosion, and this indirectly confirmed our viewpoint. We also found that runoff velocity was the hydrodynamic parameter most closely related to soil erosion through the correlation analysis (Table 8). We found that the runoff resistance decreased with the increase in mound numbers (Table 3), mainly because the increase in mounds significantly reduced the vegetation cover, which in turn led to a decline in runoff resistance, which agrees with the results of Zhao et al. [26] and Sun et al. [49] Furthermore, although increasing mound numbers may disrupt the runoff connectivity and increase landscape roughness, it also provides runoff channels for overland flow, and with an increase in the average slope, the runoff resistance tends to decrease. In addition, we observed that the height of the plateau pika mounds settled by 1.36 cm h−1 during the experiment, but Li et al. [12] found that the height of zokor mounds decreased at a rate of 1.8 cm h−1 in a previous hydraulic erosion test on plateau zokor mounds. We speculate that this might be related to the soil properties, slope, and rainfall devices on the one hand, and on the other hand, it was mainly affected by the shape and size of the two types of mounds. Together with the experimental results and the actual landscape scale, the erosion of the pika mound soil in the degraded meadow area was not only affected by the number of mounds, slope, and vegetation cover, but was also closely related to the landscape pattern of the patches, topographic, and micro-geomorphic features as well as the adaptability of vegetation replacement, etc. The erosion mechanism and ecological succession of the pika mound patches will be scientific issues that need to be deeply investigated.

4.2. Response Mechanisms of Plateau Pika Activities to Soil Erosion and Meadow Degradation

Plateau rodent activity significantly drives soil degradation and meadow degradation in the source zone of the Yellow River [53,54]. Degraded meadow areas induced by plateau pika activities can provide a large number of material sources for soil erosion, while repeated and continuous soil erosion can, in turn, exacerbate the degree of meadow degradation [55]. During the field investigation, we found that in areas with high activity of plateau pika, their digging behavior not only destroyed the vegetation community and soil structure of native meadows [56], but also caused the meadow vegetation to rot and wither due to the accumulation and burial of loose pika mounds after turning out the shallow soil, and eventually formed pika mound patch lands that were significantly different from the native meadows [11,57]. The degraded bald patches are also subject to solar radiation, toxic grass invasion, and climate change, which constantly cause an increase in the number and area of patches, which leads to an expansion in alpine meadow degradation. It has been shown that meadow degradation and soil erosion are intrinsically related [58,59], in particular, the mound patches accumulated by the activities of plateau pika are prone to eroding under the action of wind, water, and freeze–thaw due to their loose structure, fine particles, and particular shape. In our study, degraded patches of land with different numbers of pika mounds responded positively to the soil water erosion conditions. The total soil loss of the degraded patch lands with one, two, three, four, and five pika mounds were significantly different (p < 0.05). Their soil erosion rates were significantly higher than the revegetated pika mounds (p < 0.01). Furthermore, soil erosion is the process by which fine soil particles are separated and transported and includes the loss of soil nutrients [60]. The succession of pika mound patches, due to the lack of nutrients in the vegetation, leads to the expansion of pika mound patches in all directions and even causes patch-to-patch connectivity, which eventually forms a bidirectional vicious cycle between soil erosion and meadow degradation, continuously increasing the level of meadow degradation [12]. Soil erosion and meadow degradation are extraordinarily complex and closely interrelated processes, although they are closely related to the meteorological and hydrological conditions, geographic environment, and the adaptation of patch micro-geomorphic development and vegetation replacement. However, due to the random nature of plateau pika activities and the particular soil structure of the pika mound, it is also more sensitive in influencing meadow ecological processes and soil erosion, so to a certain extent, the frequent activities of plateau pika are the primary initiation of meadow degradation and soil degradation in the source zone of the Yellow River.

4.3. Properties of Alpine Pika Mound Patches and Their Effects on Ecological Processes in Alpine Meadows

Alpine meadows in the source zone of the Yellow River are an important protective barrier for the entire Yellow River Basin and the environmental eco-systems of the Tibetan Plateau [61]. As a basic configuration of the vegetation landscape pattern in this region, the different quantitative distribution patterns of pika mound patches will change and reshape the landscape pattern of the meadow area, impacting the stability and sustainability of the meadow ecological processes. In this study, we simulated fresh uncrusted pika mounds of the same size, weight, and shape but with different number classes, aiming to clarify the effects of different pika mound numbers on soil water erosion and its hydrodynamic characteristics. Our results show that plateau pika activities intensify meadow degradation and increase soil erodibility by reducing vegetation cover and runoff resistance. However, the timing of mound formation, number of mounds, and the shape and size of mounds vary among different periods of plateau pika behaviors, so this study also quantified and compared the mounds with the revegetated mounds in the following years. We found significant differences (p < 0.05) in the landscape index, soil erosion and sand production, and hydrodynamic characteristics of the revegetated pika mounds compared to other pika mound patches. This reason could be related to the spatial scale instead of the temporal scale (i.e., at the early stage of plateau pika mound patch succession, the soil particles of the pika mound patches were dispersed and loose, and the landscape was severely fragmented, so that the edge length index, area index, and volume index were more extensive). The destructive mechanical effects of plateau pika on meadows involve not only meadow degradation but also soil degradation (i.e., a massive reversal of the soil properties of pika mounds before and after [62]), leading to more intense soil and nutrient loss under the action of external erosion, which has a direct negative impact on meadow ecological processes. During the field investigation, it was also found that the succession of plateau pika mounds in different years was a spatially dynamic process, with the radius, height, area, and volume of the following year’s mound soil decreasing to different degrees with time succession, but the soil erosion resistance increased significantly. The main reason for this is that although the mounding behavior of plateau pika destroys meadow vegetation and soil, the microtopographic environment formed by the pika mounds can renew the subsurface organic matter composition and create more ecological niches for the subsequent invasion and colonization of other dominant species. To a certain extent, it has contributed to the improvement in the properties of pika mound patches and their benign cycle with the ecological processes of the meadow [63]. In summary, the presence of pika mound patches has both advantages and disadvantages for meadow ecological processes, but in the context of local climatic and environmental characteristics and the current status of meadow degradation, it is still necessary to artificially control the population density and mound numbers of plateau pika in the outbreak area. In addition, we concluded that aside from the number of pika mounds, the slope and vegetation cover also are essential factors influencing the soil water erosion processes of the pika mounds. We suggest that a combination of different vegetation cover and slope patterns can be integrated into the restoration of meadow vegetation to enhance the soil and water conservation function of degraded bald patches, which is fundamentally important in promoting sustainable management and the utilization of meadow resources.

5. Conclusions

We aimed to clarify the effects of different pika mound numbers on the soil erosion and hydrodynamic characteristics by conducting simulated rainfall experiments. The results showed that the amount of soil loss increased and then decreased with rainfall time and that 15–20 min of rainfall was a sensitive period for soil loss in different pika mound patch lands. The degradation of meadows due to the activities of plateau pika is an essential factor influencing soil erosion, and the soil erosion rate is positively correlated with both the pika mound number and slope. The mean flow velocity can better describe the soil erosion process, and the magnitude of its value increases with the number of pika mounds and slope. The Reynolds number ranged from 57.85 to 153.63 (Re < 500), and it was preliminarily determined that all slope runoff was laminar. The Froude number increased with the pika mound numbers and slope, and they were both significant factors (p < 0.05). The resistance coefficient instead decreased with the pika mound numbers and slope, and the inhibitory effect of vegetation on soil erosion was likely to be limited when the number of mounds reached a certain level. Based on the results of the grey correlation method and Pearson correlation analysis, changes in the number of pika mounds caused variability in the slope surface flow hydrodynamic characteristics and soil erosion conditions by changing the landscape scale effect, where the number of mound patches (NP), edge length index (TE), area index (AREA), and volume index (V) were the key influencing factors on the soil erosion rate and each hydrodynamic parameter.
We can conclude that outbreaks of plateau pika activity exacerbate meadow degradation, and that the development of pika mound patches leads to severe landscape fragmentation through altered landscape-scale effects, and ultimately causes decreases in the vegetation cover and runoff resistance and increases in the soil erodibility. To effectively restore the meadow ecosystem, we suggest that reasonable control of plateau pika activities should be combined with a reduction in the density of pika mounds to realize the dual purpose of soil erosion control.

Author Contributions

Conceptualization, G.L., Y.L. and H.Z.; Supervision, X.H. and X.L.; Field investigations, S.T., G.L., J.Z., C.J. and J.L.; Data analysis, W.C. and J.Z.; Writing the manuscript, S.T. and G.L.; Review and editing the manuscript, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Basic Research Project of Qinghai Provincial Science and Technology Department (2021-ZJ-701); the National Natural Science Foundation of China (U21A20191, 41662023, 42161068).

Data Availability Statement

Not applicable.

Acknowledgments

All authors thank Y.L. (Yurong Li), Y.L. (Yong Li), and Y.L. (Yuxiang Long) for the field data collection; the anonymous reviewers and editorial staff for their guidance and assistance in improving this manuscript.

Conflicts of Interest

All authors declare no conflict of interest.

References

  1. Ma, Q.; Jin, H.-J.; Bense, V.F.; Luo, D.-L.; Marchenko, S.S.; Harris, S.A.; Lan, Y.-C. Impacts of degrading permafrost on streamflow in the source area of Yellow River on the Qinghai-Tibet Plateau, China. Adv. Clim. Change Res. 2019, 10, 225–239. [Google Scholar] [CrossRef]
  2. Wang, X.; Yang, T.; Yong, B.; Krysanova, V.; Shi, P.; Li, Z.; Zhou, X. Impacts of climate change on flow regime and sequential threats to riverine ecosystem in the source region of the Yellow River. Environ. Earth Sci. 2018, 77, 465. [Google Scholar] [CrossRef]
  3. Qian, D.; Du, Y.; Li, Q.; Guo, X.; Fan, B.; Cao, G. Impacts of alpine shrub-meadow degradation on its ecosystem services and spatial patterns in Qinghai-Tibetan Plateau. Ecol. Indic. 2022, 135, 108541. [Google Scholar] [CrossRef]
  4. Xu, M.; Liu, M.; Xue, X.; Zhai, D. Warming effects on plant biomass allocation and correlations with the soil environment in an alpine meadow, China. J. Arid. Land 2016, 8, 773–786. [Google Scholar] [CrossRef]
  5. Zhang, Y.; Ganjurjav, H.; Dong, S.; Gao, Q. Excessive plant compensatory growth: A potential endogenous driver of meadow degradation on the Qinghai-Tibetan Plateau. Ecosyst. Health Sustain. 2020, 6, 1816500. [Google Scholar] [CrossRef]
  6. Yi, X.; Li, G.; Yin, Y. The impacts of grassland vegetation degradation on soil hydrological and ecological effects in the source region of the Yellow River—A case study in Junmuchang region of Maqin country. Procedia Environ. Sci. 2012, 13, 967–981. [Google Scholar] [CrossRef]
  7. Pang, X.P.; Guo, Z.G. Plateau pika disturbances alter plant productivity and soil nutrients in alpine meadows of the Qinghai-Tibetan Plateau, China. Rangel. J. 2017, 39, 133–144. [Google Scholar] [CrossRef]
  8. Yang, P.N.; Li, X.L.; Li, C.Y.; Duan, C.W. Response of Soil Microbial Diversity to Long-term Enclosure in Degraded Patches of Alpine Meadow in the Source Zone of the Yellow River. Huanjing Kexue 2023, 44, 2293–2303. [Google Scholar]
  9. Chen, Y.Y.; Yang, H.; Bao, G.S.; Pang, X.P.; Guo, Z.G. Effect of the presence of plateau pikas on the ecosystem services of alpine meadows. Biogeosciences 2022, 19, 4521–4532. [Google Scholar] [CrossRef]
  10. Wei, X.; Guo, Z. Effect of Plateau Pika Disturbance on Plant Aboveground Biomass of Alpine Meadows at Two Different Scales. Plants 2022, 11, 2266. [Google Scholar] [CrossRef]
  11. Li, J.-X.; Li, X.-L.; Gao, J.; Shi, Y.; Ma, G.-L.; Rang, K.Z.C. Influences of pika and simulated grazing disturbances on bare patches of alpine meadow in the Yellow River Source Zone. J. Mt. Sci. 2021, 18, 1307–1320. [Google Scholar] [CrossRef]
  12. Li, G.; Li, X.; Li, J.; Chen, W.; Zhu, H.; Zhao, J.; Hu, X. Influences of Plateau Zokor Burrowing on Soil Erosion and Nutrient Loss in Alpine Meadows in the Yellow River Source Zone of West China. Water 2019, 11, 2258. [Google Scholar] [CrossRef]
  13. Liu, J.; Chen, J.; Qin, Q.; You, H.; Han, X.; Zhou, G. Patch Pattern and Ecological Risk Assessment of Alpine Grassland in the Source Region of the Yellow River. Remote Sens. 2020, 12, 3460. [Google Scholar] [CrossRef]
  14. Li, X.; Zhang, Y.; Yan, H.; Salahou, M.K. Watershed-level spatial pattern of degraded alpine meadow and its key influencing factors in the Yellow River Source Zone of West China. Ecol. Indic. 2023, 146, 109865. [Google Scholar] [CrossRef]
  15. Jiang, X.-J.; Zhu, X.; Yuan, Z.-Q.; Li, X.G.; Liu, W.; Zakari, S. Lateral flow between bald and vegetation patches induces the degradation of alpine meadow in Qinghai-Tibetan Plateau. Sci. Total Environ. 2020, 751, 142338. [Google Scholar] [CrossRef]
  16. Hou, L.; Wu, F.; Xie, X. The spatial characteristics and relationships between landscape pattern and ecosystem service value along an urban-rural gradient in Xi’an city, China. Ecol. Indic. 2020, 108, 105720. [Google Scholar] [CrossRef]
  17. Xie, S.; Yan, D.; Li, J.; Liu, Y.; Sheng, Y.; Luan, Z. GEE-Based Spatial-Temporal Dynamics in a Ramsar Wetland, Honghe National Nature Reserve, Northeast China from 1985 to 2021. Land 2022, 11, 2137. [Google Scholar] [CrossRef]
  18. Romano, N.; Lignola, G.; Brigante, M.; Bosso, L.; Chirico, G. Residual life and degradation assessment of wood elements used in soil bioengineering structures for slope protection. Ecol. Eng. 2016, 90, 498–509. [Google Scholar] [CrossRef]
  19. Li, Z.; Lyu, S.; Chen, H.; Ao, Y.; Zhao, L.; Wang, S.; Zhang, S.; Meng, X. Changes in climate and snow cover and their synergistic influence on spring runoff in the source region of the Yellow River. Sci. Total Environ. 2021, 799, 149503. [Google Scholar] [CrossRef]
  20. Zhang, J.; Liu, D.; Meng, B.; Chen, J.; Wang, X.; Jiang, H.; Yu, Y.; Yi, S. Using UAVs to assess the relationship between alpine meadow bare patches and disturbance by pikas in the source region of Yellow River on the Qinghai-Tibetan Plateau. Glob. Ecol. Conserv. 2021, 26, e01517. [Google Scholar] [CrossRef]
  21. Ji, J.; Zhao, Y.; Zhang, W.; Gu, K.; Li, W. Quantitative assessment of biocrust distribution patterns using landscape indices benefits the study of their soil conservation functions. Geoderma 2023, 429, 116257. [Google Scholar] [CrossRef]
  22. LiLi, J.; Qi, H.H.; Duan, Y.Y.; Guo, Z.G. Effects of Plateau Pika Disturbance on the Spatial Heterogeneity of Vegetation in Alpine Meadows. Front. Plant Sci. 2021, 12, 771058. [Google Scholar] [CrossRef] [PubMed]
  23. Hancock, G.; Wells, T. Predicting soil organic carbon movement and concentration using a soil erosion and Landscape Evolution Model. Geoderma 2021, 382, 114759. [Google Scholar] [CrossRef]
  24. Zhang, S.; Fan, W.; Li, Y.; Yi, Y. The influence of changes in land use and landscape patterns on soil erosion in a watershed. Sci. Total Environ. 2017, 574, 34–45. [Google Scholar] [CrossRef] [PubMed]
  25. Sun, L.; Guo, H.; Liu, B.; Wu, S.; Weckler, P.R.; Yang, J. Characterizing erosion processes on a convex slope based on 3D reconstruction method. Geoderma 2021, 402, 115364. [Google Scholar] [CrossRef]
  26. Zhao, Q.H.; Zhang, Y.F.; Xu, S.S.; Ji, X.Y.; Wang, S.Q.; Ding, S.Y. Relationships between Riparian Vegetation Pattern and the Hydraulic Characteristics of Upslope flow. Sustainability 2019, 11, 2966. [Google Scholar] [CrossRef]
  27. Sun, T.C.; Liu, T.T.; Chu, L.; Li, Z.X.; Wang, T.W.; Cai, H.F. Effects of temporal and spatial variations in source-sink landscape patterns on soil erosion and sediment yield from typical watershed in the Three Gorges Reservoir area, China. Acta Ecol. Sin. 2019, 39, 7476–7492. [Google Scholar]
  28. Liu, X.; Huang, J.; Liu, W.; Ding, Y.; Ge, L. Research on the Coastal Landscape Pattern Index in the District of Nansha. Sustainability 2023, 15, 3378. [Google Scholar] [CrossRef]
  29. Flanagan, D.C.; Ascough, J.C., II; Nieber, J.L.; Misra, D.; Douglas-Mankin, K.R. Advances in Soil Erosion Research: Processes, Measurement, and Modeling. Trans. ASABE 2013, 56, 455–463. [Google Scholar] [CrossRef]
  30. Liu, Y.; Lv, Y.H.; Fu, B.J. Implication and limitation of landscape metrics in delineating relationship between landscape pattern and soil erosion. Acta Ecol. Sin. 2017, 37, 4923–4935. [Google Scholar]
  31. Li, C.; Pan, C. Overland runoff erosion dynamics on steep slopes with forages under field simulated rainfall and inflow. Hydrol. Process. 2020, 34, 1794–1809. [Google Scholar] [CrossRef]
  32. Wu, D.; Peng, R.; Huang, L.; Cao, W.; Huhe, T. Spatio-Temporal Analysis and Driving Factors of Soil Water Erosion in the Three-River Headwaters Region, China. Water 2022, 14, 4127. [Google Scholar] [CrossRef]
  33. Wang, F.; Xu, G.; Li, L.; Li, Z.; Li, P.; Zhang, J.; Cheng, Y. Response Relationship between Microtopographic Variation and Slope Erosion under Sand-Cover. Water 2019, 11, 2488. [Google Scholar] [CrossRef]
  34. Liu, Y.; Gao, G.; Li, H.; Liu, L.; Fan, Z.; Wen, T. Spatiotemporal Variations and Causes of Wind/Rainfall Erosion Climatic Erosivity in Qinghai Province, China. Atmosphere 2022, 13, 1649. [Google Scholar] [CrossRef]
  35. Alaboz, P.; Dengiz, O.; Demir, S.; Şenol, H. Digital mapping of soil erodibility factors based on decision tree using geostatistical approaches in terrestrial ecosystem. CATENA 2021, 207, 105634. [Google Scholar] [CrossRef]
  36. SuSu, Y.; Zhang, Y.; Wang, H.; Lei, N.; Li, P.; Wang, J. Interactive Effects of Rainfall Intensity and Initial Thaw Depth on Slope Erosion. Sustainability 2022, 14, 3172. [Google Scholar] [CrossRef]
  37. Zhang, S.Y.; Liang, Z.Q.; Xie, Z.Y.; Zhuo, M.N.; Liao, Y.S.; Guo, T.L.; Li, D.Q. Mechanisms of grass in slope erosion control in red soil region of southern China. Soil Water Conserv. 2016, 30, 1–5. [Google Scholar]
  38. Ding, W.; Li, M. Effects of grass coverage and distribution patterns on erosion and overland flow hydraulic characteristics. Environ. Earth Sci. 2016, 75, 477. [Google Scholar] [CrossRef]
  39. Chen, C.; Zhao, G.; Zhang, Y.; Bai, Y.; Tian, P.; Mu, X.; Tian, X. Linkages between soil erosion and long-term changes of landscape pattern in a small watershed on the Chinese Loess Plateau. CATENA 2023, 220, 106659. [Google Scholar] [CrossRef]
  40. Li, X.; Zhang, Y.; Ji, X.; Strauss, P.; Zhang, Z. Effects of shrub-grass cover on the hillslope overland flow and soil erosion under simulated rainfall. Environ. Res. 2022, 214, 113774. [Google Scholar] [CrossRef]
  41. Sun, J.; Fan, D.; Yu, X.; Li, H. Hydraulic characteristics of varying slope gradients, rainfall intensities and litter cover on vegetated slopes. Hydrol. Res. 2018, 49, 506–516. [Google Scholar] [CrossRef]
  42. Jin, J.; Wang, G.; Zhang, J.; Yang, Q.; Liu, C.; Liu, Y.; Bao, Z.; He, R. Impacts of climate change on hydrology in the Yellow River source region, China. J. Water Clim. Change 2020, 11, 916–930. [Google Scholar] [CrossRef]
  43. Jiang, P.; Yu, Z.; Yuan, F.; Acharya, K. The Multi-Scale Temporal Variability of Extreme Precipitation in the Source Region of the Yellow River. Water 2019, 11, 92. [Google Scholar] [CrossRef]
  44. Paruchuri, S.; Smith, A.T.; Fan, Z.; Dobson, F.S. Microhabitat use by plateau pikas: Living on the edge. J. Mammal. 2019, 100, 1221–1228. [Google Scholar] [CrossRef]
  45. Tong, S.; Li, G.; Li, X.; Li, J.; Zhai, H.; Zhao, J.; Zhu, H.; Liu, Y.; Chen, W.; Hu, X. Soil Water Erosion and Its Hydrodynamic Characteristics in Degraded Bald Patches of Alpine Meadows in the Yellow River Source Area, Western China. Sustainability 2023, 15, 8165. [Google Scholar] [CrossRef]
  46. An, J.; Zheng, F.; Lu, J.; Li, G. Investigating the Role of Raindrop Impact on Hydrodynamic Mechanism of Soil Erosion Under Simulated Rainfall Conditions. Soil Sci. 2012, 177, 517–526. [Google Scholar] [CrossRef]
  47. Reza, D.B.; Hamidreza, K.; Mohammadreza, Y. Investigation on ecosystem degradation induced by LU/LC changes using landscape pattern indices analysis. Arab. J. Geosci. 2018, 11, 443. [Google Scholar]
  48. Kwon, O.S.; Ra, J.H.; Ku, J.N.; Kim, J.H. Temporal-Spatial Analysis of Landscape Diversity using FRAGSTATS. J. Korea Soc. Environ. Restor. Technol. 2015, 18, 39–50. [Google Scholar] [CrossRef]
  49. Sun, R.; Ma, L.; Zhang, S.; Yu, Y.; Shen, M.; Zhang, H.; Wang, D.; Yang, Y.; Zhang, J.; Zhang, Y.; et al. Study on Landscape Patches Influencing Hillslope Erosion Processes and Flow Hydrodynamics in the Loess Plateau of Western Shanxi Province, China. Water 2020, 12, 3201. [Google Scholar] [CrossRef]
  50. Yang, Y.; Yang, S.; Wang, X.; Wang, C.; Wang, J. Study on the influence of ecological vegetation revetment on river flow and sediment environment. Water Supply 2020, 20, 3141–3155. [Google Scholar] [CrossRef]
  51. Marzen, M.; Iserloh, T.; Casper, C.M.; Ries, B.J. Quantification of particle detachment by rain splash and wind-driven rain splash. CATENA 2015, 127, 135–141. [Google Scholar] [CrossRef]
  52. Yuan, H.D.; Xin, Z.B.; Jiang, Q.L.; Yu, X.X.; Fan, D.X. Slope erosion and its hydrodynamic characteristic of Cinnamon soil under continuous rainfall. Soil Water Conserv. 2020, 34, 14–20+30. [Google Scholar]
  53. Hu, L.; Zi, H.; Ade, L.; Lerdau, M.; Wang, C. Effects of zokors (Myospalax baileyi) on plant, on abiotic and biotic soil characteristic of an alpine meadow. Ecol. Eng. 2017, 103, 95–105. [Google Scholar] [CrossRef]
  54. Qin, Y.; Yi, S.H.; Ding, Y.J.; Qin, Y.; Zhang, W.; Sun, Y.; Hou, X.M.; Yu, H.Y.; Meng, B.P.; Zhang, H.F.; et al. Effects of plateau pikas’ foraging and burrowing activities on vegetation biomass and soil organic carbon of alpine grasslands. Plant Soil 2020, 458, 201–216. [Google Scholar] [CrossRef]
  55. Li, G.R.; Li, X.L.; Li, J.F.; Zhu, H.L.; Chen, W.T.; Zhao, J.Y.; Hu, X.S.; Cui, Y.P.; Li, C.L. Soil wind erosion law in Ochotona Curzoniae mound of alpine meadow in the Yellow River. Soil Water Conserv. 2019, 33, 110–114+168. [Google Scholar]
  56. Wei, W.R.; He, J.D.; Zheng, Q.Y. Plateau pikas (Ochotona curzoniae) at low densities have no destructive effect on winter pasture in alpine meadows. Rangel. J. 2020, 42, 55–61. [Google Scholar] [CrossRef]
  57. Brierley, G.; Li, X.-L.; Fryirs, K.; Gao, J.; Shi, Y.; Perry, G.L.W.; Cullum, C. Degradation and recovery of alpine meadow catenas in the source zone of the Yellow River, Western China. J. Mt. Sci. 2022, 19, 2487–2505. [Google Scholar] [CrossRef]
  58. Hou, J.; Fu, B.; Wang, S.; Zhu, H. Comprehensive analysis of relationship between vegetation attributes and soil erosion on hillslopes in the Loess Plateau of China. Environ. Earth Sci. 2014, 72, 1721–1731. [Google Scholar] [CrossRef]
  59. Zhao, X.; Zhang, B.; Wu, P. Changes in key driving forces of soil erosion in the Middle Yellow River Basin: Vegetation and climate. Nat. Hazards 2013, 70, 957–968. [Google Scholar] [CrossRef]
  60. Sun, H.; Li, X.; Jin, L.; Su, X.; Li, C.; Kou, J. Effects of biological soil crusts on soil labile organic carbon of patchy alpine meadows in the Source Zone of the Yellow River, West China. CATENA 2023, 220, 106715. [Google Scholar] [CrossRef]
  61. Chen, Q.; Lu, Q.S.; Gao, Z.Q. Changes in Herdsmen’s Pastoral Behaviour Triggered by Rangeland Degradation in the Source Region of the Yellow River, Qinghai-Tibetan Plateau. Outlook Agric. 2015, 44, 81–87. [Google Scholar] [CrossRef]
  62. Li, G.; Li, X.; Chen, W.; Li, J.; Zhu, H.; Hu, X.; Zhou, H.; Sun, H. Effects of degradation severity on the physical, chemical and mechanical properties of topsoil in alpine meadow on the Qinghai-Tibet Plateau, west China. CATENA 2019, 187, 104370. [Google Scholar] [CrossRef]
  63. Xiang, Z.; Bhatt, A.; Tang, Z.; Peng, Y.; Wu, W.; Zhang, J.; Wang, J.; Gallacher, D.; Zhou, S. Disturbance of plateau zokor-made mound stimulates plant community regeneration in the Qinghai-Tibetan Plateau, China. J. Arid. Land 2021, 13, 1054–1070. [Google Scholar] [CrossRef]
Figure 1. (a) Geographical location of the study area. (b) Plateau pika mound and zokor mound outbreak area. (c) Plateau pika mound creation process and basic morphology of mound patches.
Figure 1. (a) Geographical location of the study area. (b) Plateau pika mound and zokor mound outbreak area. (c) Plateau pika mound creation process and basic morphology of mound patches.
Water 15 03111 g001
Figure 2. (a) Diagram of the water erosion test box. (b) Vegetation growth on the surface of the field water erosion test boxes.
Figure 2. (a) Diagram of the water erosion test box. (b) Vegetation growth on the surface of the field water erosion test boxes.
Water 15 03111 g002
Figure 3. (a) Principle of artificial rainfall simulation device. (b) Diagram of the field hydraulic erosion test setup.
Figure 3. (a) Principle of artificial rainfall simulation device. (b) Diagram of the field hydraulic erosion test setup.
Water 15 03111 g003
Figure 4. A comparison of the total soil loss per unit area under different gradients and the pika mound numbers. All tests were repeated thrice, and the average values were used instead. Error bars are the graphical displays of the data variability and used here to represent the standard error (SE). Different capital and lowercase letters represent significant differences in the total soil loss at the same levels (p < 0.05).
Figure 4. A comparison of the total soil loss per unit area under different gradients and the pika mound numbers. All tests were repeated thrice, and the average values were used instead. Error bars are the graphical displays of the data variability and used here to represent the standard error (SE). Different capital and lowercase letters represent significant differences in the total soil loss at the same levels (p < 0.05).
Water 15 03111 g004
Figure 5. Regulation of soil loss per unit time with rainfall time under different pika mound numbers.
Figure 5. Regulation of soil loss per unit time with rainfall time under different pika mound numbers.
Water 15 03111 g005
Figure 6. Variation pattern of the soil erosion rate with rainfall duration under different gradients of pika mound number.
Figure 6. Variation pattern of the soil erosion rate with rainfall duration under different gradients of pika mound number.
Water 15 03111 g006
Figure 7. (a) Relationship between the pika mound numbers and soil erosion rate. (b) Relationship between the vegetation cover and soil erosion rate.
Figure 7. (a) Relationship between the pika mound numbers and soil erosion rate. (b) Relationship between the vegetation cover and soil erosion rate.
Water 15 03111 g007
Figure 8. (a) Function between the edge length index and soil erosion rate. (b) Function between the area index and soil erosion rate. (c) Function between the volume index and soil erosion rate. (d) Function between the separation index and soil erosion rate.
Figure 8. (a) Function between the edge length index and soil erosion rate. (b) Function between the area index and soil erosion rate. (c) Function between the volume index and soil erosion rate. (d) Function between the separation index and soil erosion rate.
Water 15 03111 g008aWater 15 03111 g008b
Table 1. Basic physio-mechanical characteristics and the soil grain size class composition of pika mounds.
Table 1. Basic physio-mechanical characteristics and the soil grain size class composition of pika mounds.
Pika Mounds TypeMoisture ContentDensityPorosityCohesionFirmnessBelow a Certain Soil Grain Size (%)
(%)(g m−3)(%)(kPa)(kPa)d < 2 mmd < 0.5 mmd < 0.075 mm
Fresh pika mounds9.35 ± 1.750.87 ± 0.1117.34 ± 0.1611.56 ± 1.0210.01 ± 1.9588.7959.728.88
Revegetated pika mounds15.16 ± 2.081.51 ± 0.1511.38 ± 0.1426.61 ± 3.2459.64 ± 2.3875.2648.978.15
Healthy meadows38.25 ± 2.161.68 ± 0.084.75 ± 0.0735.86 ± 3.53153.50 ± 5.1868.5545.397.70
Table 2. Landscape pattern indicators of pika mound patches.
Table 2. Landscape pattern indicators of pika mound patches.
Landscape IndexEquationsDescription and Meaning of Indicators
NP N P = N The total number of all pika mound patches (N) to reflect the landscape’s spatial pattern and the whole landscape’s qualitative heterogeneity.
TE T E = E The sum of the edge lengths (E) of all pika mound patches can describe the degree of agglomeration and extension trends of a given patch type.
PD P D = N A i The overall heterogeneity and fragmentation of the landscape and the degree of fragmentation of a given patch type are examined.
AREA A R E A = 2 π r h The area of all pika mound patches is the basis for calculating other landscape indicators with implications for species abundance, energy nutrients, and productivity levels.
V V = π h 6 ( 3 r 2 + h 2 ) The volume of all pika mound patches can describe the overall spatial characteristics of the landscape and accurately reflect the degree of landscape disturbance by external or human factors.
Ni N i = A 2 A i N A The level of dispersion or aggregation of the individual spatial distribution of a specific type of patch in the landscape.
Note: A (m2) is the total landscape area, Ai (m2) is the area of landscape type “i”, r (m) is the pika mound patches radius, and h (m) is the height of the pika mound patches.
Table 3. Hydrodynamic characteristics at different pika mound numbers.
Table 3. Hydrodynamic characteristics at different pika mound numbers.
Slope
(°)
The Number of Pika Mounds (pcs)Runoff Depth
(h, mm)
Flow Velocity
(v, cm s−1)
Reynolds Number
(Re)
Froude Number (Fr)Darcy–Weisbach Resistance (f)
10°Revegetated pika mound3.57 ± 0.74 Aa2.29 ± 0.27 Cf57.85 ± 18.45 Cd0.143 ± 0.009 Be159.07 ± 31.77 Aa
1 pika mound3.25 ± 0.73 Ab2.61 ± 0.26 Ce66.08 ± 20.63 Ccd0.168 ± 0.016 Cd99.10 ± 13.41 Ab
2 pika mounds3.06 ± 0.75 Ac3.20 ± 0.32 Cd75.01 ± 25.91 Cbc0.236 ± 0.024 Bc71.56 ± 13.54 Ac
3 pika mounds2.83 ± 0.76 Ad3.56 ± 0.33 Cc80.79 ± 31.53 Cb0.281 ± 0.037 Bb51.89 ± 9.71 Ad
4 pika mounds2.75 ± 0.74 Ade3.95 ± 0.34 Cb76.36 ± 35.42 Cbc0.289 ± 0.049 Cb38.86 ± 8.56 Ae
5 pika mounds2.85 ± 0.77 Ad4.26 ± 0.29 Ca93.15 ± 39.01 Ca0.322 ± 0.073 Ca34.89 ± 6.56 Ae
20°Revegetated pika mound3.33 ± 0.72 Ba2.65 ± 0.32 Bf86.55 ± 25.82 Bd0.157 ± 0.010 Af125.80 ± 17.73 Ba
1 pika mound3.00 ± 0.72 Bb3.03 ± 0.35 Be89.75 ± 27.57 Bd0.180 ± 0.018 Be86.71 ± 16.54 Bb
2 pika mounds2.82 ± 0.76 Bc3.56 ± 0.33 Bd99.19 ± 33.12 Bc0.271 ± 0.021 Ac57.98 ± 9.70 Bc
3 pika mounds2.66 ± 0.79 Bd4.05 ± 0.42 Bc114.18 ± 38.56 Bb0.253 ± 0.040 Cd42.56 ± 9.68 Bd
4 pika mounds2.45 ± 0.70 Be4.63 ± 0.37 Bb105.39 ± 37.84 Bc0.310 ± 0.045 Bb29.81 ± 6.33 Be
5 pika mounds2.61 ± 0.71 Bd5.19 ± 0.42 Ba131.92 ± 42.18 Ba0.361 ± 0.065 Ba34.51 ± 5.41 Ade
30°Revegetated pika mound2.79 ± 0.61 Ca2.98 ± 0.40 Af105.12 ± 28.83 Ad0.161 ± 0.017 Af71.47 ± 8.71 Ca
1 pika mound2.56 ± 0.63 Cb3.56 ± 0.38 Ae116.64 ± 34.92 Acd0.203 ± 0.016 Ae50.29 ± 7.68 Cb
2 pika mounds2.37 ± 0.66 Cc4.36 ± 0.39 Ad134.71 ± 37.62 Abc0.270 ± 0.026 Ad30.55 ± 5.33 Cc
3 pika mounds2.29 ± 0.73 Ccd4.57 ± 0.35 Ac125.36 ± 41.09 Ac0.293 ± 0.034 Ac23.66 ± 5.11 Ccd
4 pika mounds2.38 ± 0.79 Bc5.22 ± 0.38 Ab141.86 ± 45.26 Ab0.329 ± 0.047 Ab17.87 ± 4.76 Cd
5 pika mounds2.18 ± 0.82 Cd5.82 ± 0.44 Aa153.63 ± 51.49 Aa0.385 ± 0.055 Aa16.99 ± 4.76 Bd
Note: All tests were repeated thrice, and 12 datasets were taken each time. The values in the table indicate the mean and standard deviation (SD). Different capital and lowercase letters represent significant differences in the individual hydrodynamic parameters at the same level (p < 0.05).
Table 4. Regression equation of the number of pika mounds with the hydrodynamic parameters at different slopes.
Table 4. Regression equation of the number of pika mounds with the hydrodynamic parameters at different slopes.
Hydrodynamic ParametersThe Number of Pika Mounds (pcs)
10° Slope20° Slope30° Slope
vy = 0.4066x + 2.2952(R2 = 0.9917)y = 0.514x + 2.5667 (R2 = 0.9963)y = 0.554x + 3.0333 (R2 = 0.9875)
Rey = 60.234e0.0826x (R2 = 0.8717)y = 85.092e0.078x (R2 = 0.8532)y = 108.2e0.0689x (R2 = 0.8825)
Fry = 0.0433x + 0.1663 (R2 = 0.9782)y = 0.0398x + 0.1559 (R2 = 0.9305)y = 0.0374x + 0.1474 (R2 = 0.9418)
fy = 140.93e–0.306x (R2 = 0.9738)y = 112.19e–0.285x (R2 = 0.9653)y = 64.547e–0.301x (R2 = 0.9733)
Note: y is the hydrodynamic parameter, x is the number of pika mounds.
Table 5. Landscape pattern indicators of the plateau pika mound patches.
Table 5. Landscape pattern indicators of the plateau pika mound patches.
Slope (°)The Number of Pika Mounds (pcs)NPTEAREAVPDNirh
10°Revegetated pika mound1 e1.130 f0.0712 f0.0033 f14.042 a14.894 a0.180 c0.063 b
1 pika mound1 e1.413 e0.1130 e0.0066 e8.846 c9.383 b0.225 ab0.080 a
2 pika mounds2 d2.700 d0.2106 d0.0118 d9.495 b7.121 c0.215 b0.078 a
3 pika mounds3 c4.145 c0.3357 c0.0193 c8.936 c5.472 d0.220 b0.081 a
4 pika mounds4 b5.903 b0.4841 b0.0296 b8.263 d4.382 e0.235 a0.082 a
5 pika mounds5 a7.065 a0.5864 a0.0345 a8.527 cd4.045 e0.225 ab0.083 a
20°Revegetated pika mound1 e1.036 f0.0632 f0.0027 f15.821 a16.780 a0.165 c0.061 b
1 pika mound1 e1.382 e0.1147 e0.0066 e8.721 c9.249 b0.220 b0.083 a
2 pika mounds2 d2.889 d0.2311 d0.0138 d8.654 cd6.491 c0.230 a0.080 a
3 pika mounds3 c4.239 c0.3476 c0.0204 c8.631 c5.285 d0.225 ab0.082 a
4 pika mounds4 b5.903 b0.4664 b0.0284 b8.577 d4.549 e0.235 a0.079 a
5 pika mounds5 a7.222 a0.5561 a0.0332 a8.991 b4.265 e0.230 a0.077 a
30°Revegetated pika mound1 e1.005 f0.0553 f0.0022 f18.095 a19.193 a0.160 c0.055 b
1 pika mound1 e1.444 e0.1156 e0.0069 e8.654 c9.179 b0.230 ab0.080 a
2 pika mounds2 d2.700 d0.2133 d0.0120 d9.375 b7.031 c0.215 b0.079 a
3 pika mounds3 c4.427 c0.3675 c0.0225 c8.164 d4.999 d0.235 a0.083 a
4 pika mounds4 b6.029 b0.4582 b0.0284 b8.730 c4.630 de0.240 a0.076 a
5 pika mounds5 a6.908 a0.5388 a0.0309 a9.279 b4.402 e0.220 b0.078 a
Note: NP = the number of pika mound patches; TE = edge length index; AREA = area index; V = volume index; PD = patch density index; Ni = separability index; r = radius of pika mound patches; h = height of pika mound patches. Different lowercase letters in the same column indicate significant differences between the same indices (p < 0.05).
Table 6. Correlation analysis of the soil erosion rate with the landscape pattern indicators of pika mound patches.
Table 6. Correlation analysis of the soil erosion rate with the landscape pattern indicators of pika mound patches.
Slope (°)CorrelationNPAREATEVPDNirh
10°Grey correlation0.9470.9320.9380.9080.6260.5470.6810.685
Pearson correlation0.992 **0.997 **0.995 **0.995 **−0.693−0.896 *0.6910.741
Sig. (two-tailed)00000.1270.0160.1280.092
20°Grey correlation0.9480.9580.9480.9490.6610.5840.7100.702
Pearson correlation0.990 **0.998 **0.995 **0.997 **−0.595−0.843 *0.7060.416
Sig. (two-tailed)00000.2130.0350.1170.413
30°Grey correlation0.9190.9080.9210.8830.6670.5940.7420.731
Pearson correlation0.984 **0.994 **0.990 **0.987 **−0.624−0.840 *0.6290.555
Sig. (two-tailed)00000.1860.0360.1810.253
Note: RSE = soil erosion rate. The definition of landscape pattern index in the table is the same as Table 5. * is a significant correlation (0.01 < p < 0.05), and ** is a highly significant correlation (p < 0.01).
Table 7. Correlation analysis of the hydrodynamic parameters with the landscape pattern indicators of pika mound patches.
Table 7. Correlation analysis of the hydrodynamic parameters with the landscape pattern indicators of pika mound patches.
Slope (°) NPAREATEVPDNirh
10°v0.977 **0.991 **0.986 **0.987 **−0.634−0.853 *0.6380.566
H−0.841 *−0.908 *−0.884 *−0.928 **0.817 *0.952 **−0.842 *−0.758
Re0.902 *0.916 *0.893 *0.895 *−0.658−0.841 *0.5910.650
Fr0.943 **0.972 **0.955 **0.970 **−0.663−0.873 *0.6560.622
f−0.851 *−0.910 *−0.885 *−0.920 **0.854 *0.979 **−0.846 *−0.802
20°v0.991 **0.998 **0.996 **0.998 **−0.592−0.842 *0.7060.408
H−0.848 *−0.899 *−0.873 *−0.910 *0.8050.954 **−0.885 *−0.672
Re0.907 *0.903 *0.899 *0.892 *−0.509−0.7590.5830.368
Fr0.954 **0.959 **0.961 **0.960 **−0.564−0.819 *0.6990.356
f−0.854 *−0.904 *−0.878 *−0.914 *0.829 *0.979 **−0.904 *−0.695
30°v0.970 **0.977 **0.973 **0.968 **−0.643−0.855 *0.6310.576
H−0.833 *−0.877 *−0.844 *−0.869 *0.7960.934 **−0.730−0.792
Re0.907 *0.902 *0.903 *0.884 *−0.626−0.823 *0.5820.557
Fr0.958 **0.957 **0.958 **0.944 **−0.617−0.831 *0.5990.543
f−0.862 *−0.915 *−0.892 *−0.923 **0.8100.960 **−0.799−0.76
Note: The hydrodynamic parameters in this table are as defined in Table 3, and the landscape pattern indices in this table are as defined in Table 5. * is a significant correlation (0.01 < p < 0.05), and ** is a highly significant correlation (p < 0.01).
Table 8. The correlation analysis between soil erosion and the hydrodynamic parameters at different slopes.
Table 8. The correlation analysis between soil erosion and the hydrodynamic parameters at different slopes.
Slope (°)CorrelationHvReFrf
10°Grey correlation0.6720.7980.7600.8170.571
Pearson correlation−0.887 *0.991 **0.919 **0.957 **−0.905 *
Sig. (two-tailed)0.01800.0110.0030.013
20°Grey correlation0.6900.8360.7930.8340.587
Pearson correlation−0.907 *0.998 **0.919 **0.953 **−0.914 *
Sig. (two-tailed)0.01300.0120.0040.011
30°Grey correlation0.6740.7860.7290.8170.569
Pearson correlation−0.910 *0.993 **0.938 **0.994 **−0.933 **
Sig. (two-tailed)0.01200.0060.0010.007
Note: The abbreviations representing the hydrodynamic parameters are the same as in Table 7. * is a significant correlation (0.01 < p < 0.05), and ** is a highly significant correlation (p < 0.01).
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

Tong, S.; Li, G.; Li, J.; Li, X.; Jiang, C.; Zhao, J.; Zhu, H.; Liu, Y.; Chen, W.; Hu, X. Influence of the Plateau Pika Mound Numbers on Soil Water Erosion Properties in Alpine Meadows of the Yellow River Source Zone, Western China. Water 2023, 15, 3111. https://doi.org/10.3390/w15173111

AMA Style

Tong S, Li G, Li J, Li X, Jiang C, Zhao J, Zhu H, Liu Y, Chen W, Hu X. Influence of the Plateau Pika Mound Numbers on Soil Water Erosion Properties in Alpine Meadows of the Yellow River Source Zone, Western China. Water. 2023; 15(17):3111. https://doi.org/10.3390/w15173111

Chicago/Turabian Style

Tong, Shengchun, Guorong Li, Jinfang Li, Xilai Li, Chengdong Jiang, Jianyun Zhao, Haili Zhu, Yabin Liu, Wenting Chen, and Xiasong Hu. 2023. "Influence of the Plateau Pika Mound Numbers on Soil Water Erosion Properties in Alpine Meadows of the Yellow River Source Zone, Western China" Water 15, no. 17: 3111. https://doi.org/10.3390/w15173111

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