Next Article in Journal
Characterization of Surface Evidence of Groundwater Flow Systems in Continental Mexico
Next Article in Special Issue
The Use of Various Rainfall Simulators in the Determination of the Driving Forces of Changes in Sediment Concentration and Clay Enrichment
Previous Article in Journal
Challenges in Harmonized Environmental Impact Assessment (EIA), Monitoring and Decommissioning Procedures of Offshore Platforms in Adriatic-Ionian (ADRION) Region
Previous Article in Special Issue
Can Lumped Characteristics of a Contributing Area Provide Risk Definition of Sediment Flux?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Slope Erosion and Hydraulics during Thawing of the Sand-Covered Loess Plateau

1
State Key Laboratory of Eco-hydraulics in Northwest Arid Region of China, Xi’an University of Technology, Xi’an 710048, China
2
Key Laboratory of National Forestry Administration on Ecological Hydrology and Disaster Prevention in Arid Regions, Xi’an 710048, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(9), 2461; https://doi.org/10.3390/w12092461
Submission received: 20 July 2020 / Revised: 25 August 2020 / Accepted: 31 August 2020 / Published: 1 September 2020
(This article belongs to the Special Issue Soil Water Erosion)

Abstract

:
Seasonal freeze-thaw processes have led to severe soil erosion globally. Slopes are particularly susceptible to changes in runoff, it can be useful to study soil erosion mechanisms. We conducted meltwater flow laboratory experiments to quantify the temporal and spatial distribution of hydraulic parameters on sandy slopes in relation to runoff and sediment yield under constant flow, different soil conditions (unfrozen slope: US; frozen slope: FS), and variable sand thickness. The results showed that sand can prolong initial runoff time, and US and FS have significantly different initial runoff times. There was a significant linear relationship between the cumulative runoff and the cumulative sediment yield. Additionally, hydrodynamic parameters of US and FS varied with time and spatially, as the distance between US and FS is linearly related to the top of the slope. We found that the main runoff flow pattern was composed of laminar flow and supercritical flow. There was a significant linear relationship between flow velocity and hydraulic parameters. The flow velocity is the best hydraulic parameter to simulate the trend of slope erosion process. This study can provide a scientific basis for a model of slope erosion during thawing for the Loess Plateau.

1. Introduction

The hydrodynamic properties of slopes have a decisive effect on runoff and sediment yield. Their study can help in understanding the process and mechanism of slope soil erosion and understanding the parameters of slope water dynamics, which are helpful in the construction of predictive models of slope soil erosion [1,2,3].
The Loess Plateau is a sand-covered landform that experiences substantial wind and water erosion [4,5,6]. Due to the differences in physical characteristics, infiltration, hydraulic conductivity, and water holding capacity between the surface sand layer and loess layer, a distinct sand-soil interface is formed between the sand layer and loess layer, and then the typical sand-soil dual structure is formed [7,8,9,10]. Zhang et al. [11] found in field rainfall experiments that the runoff and sediment yield processes on the sand-covered slopes are significantly different from those on loess slopes. Under light rain, the sand-covered slopes store rainfall, minimizing runoff, and the sediment content in any runoff that does occur is very large. Wu et al. [12] conducted a qualitative description of the sand-soil interface flow on the sand-covered slope through field surveys. Many scholars have studied the relationship between runoff and sediment yield, erosion processes, and the influence of sand layer size composition on runoff and sediment yield process through laboratory simulated rainfall experiments [6,7,8]. Tang et al. [13,14] quantitatively studied the spatiotemporal distribution of hydraulic parameters under different rainfall intensities and different sand thickness and their relationship with runoff and sediment yield. This study showed that the Reynolds Number can characterize the process of runoff and sediment yield on sand-covered slopes well [15]. However, studies on hydraulic characteristics of this particular landform at the slope scale during soil erosion is relatively limited.
The area of wind-water erosion on the Loess Plateau is in the middle latitude of the temperate zone, with an annual average precipitation of 300–600 mm. About 1/3 of days are below 0 °C each year, and it is windy and sandy in winter and spring with heavy rain in summer [16,17,18,19]. The hilly-gullied loess region is significantly affected by freezing and thawing, and snowmelt runoff erosion is a major manifestation of freeze-thaw erosion during thawing [20,21]. This erosion process is also found elsewhere throughout the world. In inland northeastern Oregon, 86% of soil erosion events are caused by freeze-thaw processes and snowmelt runoff [22]. In the northwest coastal region of the United States along the Pacific, over 90% of the total annual snow erosion is caused by melting snow [20]. When the frozen soil thaws, its shear strength decreases, and its erodibility increases, thus making the soil in the thawing period more susceptible to erosion [23,24,25,26]. The results of rainfall experiments by Sharratt et al. [27] showed that the impermeable frozen layer of soil prevents water infiltration during thawing, resulting in increased soil surface moisture content, surface runoff, with a high sediment content and soil erosion. The freeze-thaw process changes the structure of the topsoil and thereby influences the water erosion process as well [28,29]. At the same time, part of the hilly-gullied loess region subject to freeze-thaw erosion overlaps with the flakes of sand. Although this area of overlap is not large, it is widely distributed. Due to the overlapping of several types of erosion, “wind erosion, water erosion and freeze-thaw erosion”, this superposition effect causes very serious soil erosion [30]. At present, scholars have done a lot of research on single-force erosion and wind-water composite erosion on the Loess Plateau. However, little research has been done on the problems of soil and water loss caused by multiple erosion types in the above areas.
Therefore, we studied sand-covered loess slopes using laboratory scouring experiments. We analyzed soil erosion characteristics and the spatiotemporal variation in hydraulic parameters of the sand-covered slopes when frozen or not, and described the relationship between hydraulic parameters and runoff and sediment yield. This study can provide a scientific basis for the construction of erosion forecast models for sand-covered loess slopes during thawing periods.

2. Materials and Methods

2.1. Material and Device

Two soil types were gathered from the Wangmaogou watershed (37°34′13′′–37°36′03′′ N, 110°20′46′′–110°22′46′′ E) of the Loess Plateau in Suide county, Shaanxi Province, China (Figure 1). Soil particle size was measured using a Mastersizer 2000 sediment particle size analyzer. The loess was comprised of 0.20% clay, 72.01% silt, and 27.79% sand. The sandy soil was comprised of 0.72% clay, 14.38% silt, and 84.9% sand. The soil was identified as a silt loam according to the soil classification standard of the United States Department of Agriculture. The dry bulk density of the soil was 1.3 g/cm3, its organic matter content was 2.0 g/kg, and its saturated water content was 46.41%.
To measure soil erosion characteristics, we used a two-part experimental device that consisted of the frozen soil system and the scour experiment (Figure 2). The frozen soil system adopts a freeze-thaw test system implemented by the Xi’an University of Technology. The internal dimensions of the freeze-thaw were 4.5 m (length) × 2.5 m (width) × 2.5 m (height). Its internal temperature varied from −30–40 ± 1 °C. The temperature error was less than 2.0 °C, and there was a refrigeration and heating system, to maintain the experimental conditions. The system contained a runoff collection unit, soil box, sink, steady flow flume, and a water tank. The size of the soil box is 2 m long, 0.2 m wide, and 0.2 m deep with marks at 0.5 m increments on the side of the box to measure section velocities. From the top of the slope to the slope are S1, S2, S3, and S4. A sink of 2 m length, 0.2 m width, and 0.05 m depth was joined on the top of the soil box to make a stable concentrated flow.

2.2. Experimental Design

The scour experiments were conducted at the State Key Laboratory of Eco-hydraulics in the northwest arid region of China (Xi’an University of Technology) in Xi’an in April 2017. We used a local standard runoff plot and calculated the experimental flow rate from the mean precipitation during thawing as 1 L/min after correction. The scour experiment included unfrozen (US) and frozen slopes (FS) treatments and four sand thicknesses (0, 1, 2, and 3 cm). To ensure that the initial conditions of the scour experiments were consistent, the slope grade was set to 12°. The rainfall temperature was maintained at 10 °C and remained constant. Each experiment was replicated three times and we present the average of the three replicates. Each scour experiment lasted 15 min after the flow started (Table 1). Before the experiments, the soil samples were air-dried and passed through a 10 × 10 mm sieve to remove impurities such as plant roots. The soil samples were then moistened to a water content of 15% and covered with plastic film for 24 h to evenly distribute the soil moisture. A layer of gauze was laid on the bottom of the soil tank before the tank was filled with a 5 cm layer of sand. To mimic the soil’s dry bulk density, the box was filled with soil in 5 cm layer intervals, and layers were mixed. To avoid confounds of being introduced by the box itself, the slope was designed to ensure erosion would flow through the middle of the box by lowering the slope and raising the sides. Once the box was filled, sand of different thicknesses (0, 1, 2, 3 cm) was layered above the soil in the box, along with water to increase the sand’s water content. The soil tank was then frozen at −20 °C for 24 h and then placed on a scouring device bracket for testing. Since the room’s temperature was higher than that of the freeze-thaw experiment system, the scour experiment also thawed.
The flow rate was determined before the experiment began to ensure the difference between the actual flow and the intended treatment flow was less than 5% for three consecutive trials. After error testing, a formal test was performed by recording the initial flow time and collecting runoff and sediment samples every minute. Dye tracing (KMnO₄) was used to measure the flow velocity of different sections of the box by dividing the travel distance by the mean traveling time multiplied by an adjustment coefficient of 0.65. Samples were collected and the sediment was separated and then dried at 105 °C for 24 h and subsequently weighed.

2.3. Hydraulics Parameter Calculation and Methods

2.3.1. Calculation of Hydraulics Parameters

In this study, flow velocity (V), the Reynolds number (Re), Froude number (Fr), and the Darcy–Weisbach roughness coefficient (f) were selected as the research objects [31], then were calculated by the expression:
V = V m × 0.65
where V is the mean flow velocity (m/s), Vm is the observed velocity (m/s), and the flow travel distance is divided by the mean travelling time:
R e = V R ν
where ν is the kinematic viscosity (m2/s), and R is the hydraulic radius (m), which can be replaced by the value of average flow depth h:
h = Q V b T
where Q is the total runoff in time T (m3/s) and b is the width of water surface (m):
ν (m2/s) is the kinematic viscosity, it was calculated as follows:
v = 0.01775 1 + 0.0337 t + 0.000221 t 2
where t is the water temperature (°C):
F r = V g h
where g is the acceleration due to gravity (m/s2), being 9.8 m/s2:
f = 8 g h sin α V 2
where α is the slope (°).
The coefficient of variation (CV) indicates the degree of data dispersion, the formula is as follows:
C V = ( S D / M e a n ) × 100 %
where CV is the coefficient of variation (%), SD is the standard deviation, and Mean is the average value.

2.3.2. Methods

The Photoshop (Adobe Photoshop CS4 Extended 11.0.1) was applied to design the experimental system. All statistical analyses were conducted using SPSS (IBM SPSS Statistics Version 21). Figures were generated in Origin 8.5.

3. Results

3.1. Erosion, Runoff, and Sediment Yield

3.1.1. Characteristics Values of Runoff and Sediment Yield

Table 2 showed the characteristic values of runoff and sediment yield under different treatments. Based on the initial runoff time of U0, the change in the initial runoff time under different treatments was calculated. It is calculated that the initial runoff time of U1, U2, and U3 is significantly longer than that of U0, and the initial runoff time of the slope surface under different sand thicknesses has been extended by 3.5 (U1), 4.73 (U2), and 6.36 (U3) times. The initial runoff time of F0 is 37.9% earlier than U0. The initial runoff time of F1, F2, and F3 did not change much compared with U0, but compared with U1, U2, and U3, the initial runoff time was much longer. Sand-covered slopes prolong initial runoff times and the effects become longer as the thickness of sand-cover increases. The initial runoff time of FS was significantly shorter than US. The total runoff under different treatments increased in the following order: U0 < U3 < U1 < U2 < F3 < F2 < F0 < F1. The CV of runoff under different treatments was between 2.48% and 22.14%, and the fluctuation range of the runoff was small, indicating that the impact of sand cover and soil freezing on the slope runoff process is small. The total sediment yield across different treatments declined in the order U0 < U1 < U2 < U3 < F0 < F1 < F2 < F3. Based on the total sediment yield of U0, the total sediment yield under different treatments are 3.37 (U1), 4.35 (U2), 4.96 (U3), 8.38 (F0), 8.88 (F1), 10.85 (F2), and 10.98 (F3) times. The CV of the sediment yield of US was between 27.8% and 63.3%, which indicates that the sediment yield process of US had a large degree of fluctuation. The increasing sand thickness increased the CV which indicates that the sediment yield of the slope varies drastically. The CV of FS was between 3.33% and 35.22%. Under the same conditions of sand cover thickness, the CV of FS was much smaller than that of US, which indicates that the sediment yield of FS was relatively stable.
Figure 2 showed the eroded topography under different treatments, with significant differences in surface morphology. On the unfrozen slopes and frozen slopes with different sand thicknesses, a rill appeared during the runoff process. However, under the same hydraulic conditions, the rill appeared in different shapes. For the US, the rill initially developed on the slope top and bottom, which extended to the slope middle at the same time. In U0, the rill showed a discontinuous distribution, and the development of the rill was primarily on the top (S1) and bottom (S4) of the slope. In U1, U2, and U3, the rills were continuous but shallow in depth. For the FS, the rill had the same characteristics. During the experimental processes, the rill initially only developed on the slope top and gradually extended to the slope bottom. The connected rill gradually appeared on the frozen slope (Figure 3).

3.1.2. Correlation between Accumulative Runoff and Accumulative Sediment Yield

A function is fitted to the relationship between cumulative runoff and cumulative sediment yield for each experiment. The fitted equation is M = CQ + D, where M (kg) is the cumulative sediment yield, Q (L) is the cumulative runoff, and C and D are regression coefficients. All fitting equations were significant at p < 0.001 (Table 3). The coefficient C was defined as the sediment yield coefficient. MC was the mean of regression coefficient C. This coefficient obeys a certain change law. The MC of the US and FS were 0.31 and 0.67, respectively. The MC of the FS was 2.16 times than that of the US, indicating that the dependency of sediment yield on runoff was stronger for the FS than for the US.

3.2. Hydraulics of Slope Runoff

3.2.1. Spatiotemporal Variations of Flow Velocity

For the US, the flow velocity varies from 0.23 to 0.35 m/s during the tests, and its fluctuation range is small (Figure 4a). For the FS, the flow velocity varies from 0.18 to 0.35 m/s during the test, and its fluctuation range is relatively large (Figure 4). Under the condition of the same sand thickness, the mean values of flow velocity on the FS were 85.92% (F0/U0), 96.13% (F1/U1), 84.84% (F2/U2), and 88.47% (F3/U3) of the US, respectively. During the entire experiment, the flow velocity of the US and FS generally showed a downward trend (Figure 4). However, due to the conversion from erosion between inter-rill erosion to rill erosion during the erosion phase, the flow velocity changed due to the occurrence of rills. In the early stage of runoff, erosion was mainly between inter-rill erosion, the slope was relatively smooth, the runoff resistance was small, and the flow velocity was large. When the rill was generated on the slope, the flow velocity decreased significantly, and it occasionally rose with the backwater and reaches unconnected rills. In later stages of the experiment, the flow velocity tended to stabilize because the rill no longer developed. The variability of the flow velocity was mainly due to the increased resistance caused by the ground surface, the runoff energy consumption caused by the water flow down-cut, and the collapse of the soil on the side of the rill. Due to the abundance and looseness of sand, the fluctuation of U3 flow velocity is more severe.
In space, the V of the US and FS and the distance from the top of the slope can be represented by a linear function (US: R2 = 0.893, p > 0.05; FS: R2 = 0.952, p < 0.05). For the US, the V increases continuously as the distance from the top of the slope increases; for the FS, the V increases first to the maximum and then decreases as the distance from the top of the slope increases. At the same section, the V of the US decreases with the increase of the sand thickness, and then increases; the V of the FS increases with the increase in the sand thickness. Regression analysis showed that the V can be described by a linear function of the distance from the top of the slope (Table 4). The fitted equation showed that there is a certain range in which the V increases continuously as the slope length increases. Therefore, the slope length is very important. In some cases, precautions should be taken to mitigate the harm of downhill scour. Setting intercepting trenches and terraces on the slope can greatly reduce the slope length and thus reduce the erosion of runoff.

3.2.2. Spatiotemporal Variations of the Reynolds Number

From open channel hydraulics, the runoff flow is laminar when the Re < 500, the runoff flow is turbulent when the Re > 2000, and the runoff is transitional when the Re is between 500 and 2000. During the entire experiment, the Re ranged from 149.2 to 533.69 under US and FS, indicating that most of the runoff is laminar (Figure 5). For the US, the average value of the Re increased across conditions in this order: U0 < U1 < U2 < U3. There was little change in Re in the U0 treatment over time (range 155.63 to 204.82). The Re of U1, U2, and U3 increased over time, and the Re of U3 at the end of the experiment exceeded 500. For the FS, the Re in F1, F2, and F3 were approximately identical, all increasing slowly over time. The F0 increased rapidly from 0 to 5 min, and then showed a slow downward trend.
In space, the Re of the US and FS and the distance from the top of the slope can be represented by a linear function (US: R2 = 0.712, p > 0.05; FS: R2 = 0.998, p < 0.01). For the US, the Re was largest (Re > 500 for U1, U2, and U3) at 0.5 m from the top of the slope. As the distance from the top of the slope increased, the Re decreased. In the same section, the Re of the sand-covered slope was significantly larger than U0 (Table 5). For the FS, the Re was largest at 0.5 m from the top of the slope, and runoff was a transitional flow; only the Re of F0 is greater than 500. As the distance from the top of the slope increased, the Re decreased. As the distance from the top of the slope increased, the Re decreased. In the same section, the Re of F0 was greater than that of the sand-covered slope. At the same sand thickness, the Re of F0 was greater than U0. The Re of sand-covered slopes of US was larger than that of FS.

3.2.3. Spatiotemporal Variations of the Froude Number

The critical value of subcritical flow and supercritical flow is 1. If the Fr is greater than 1, it is a supercritical transition of flow, otherwise, it is a subcritical flow. The Fr gradually decreased over time (Figure 6). The Fr in the US ranged from 3.88 to 5.24, and the Fr of the FS ranged from 2.93 to 4.85. This indicated that the runoff on the slope is supercritical during the experiment. For the US, the Fr of U0 was greater than the Fr of the sand-covered slope. The changes in the Fr of U1, U2, and U3 over time were roughly the same. For the FS, the Fr decreased faster with time as compared to the US, and showed a clear layering phenomenon.
In space, the Fr of the US and FS and the distance from the top of the slope can be represented by a linear function (US: R2 = 0.899, p > 0.05; FS: R2 = 0.956, p < 0.05). The Fr under different treatments increased with the distance from the top of the slope (except U2 and F3). The maximum Fr of U2 and F3 occurred at 1.5 m from the top of the slope. For the US, the Fr of U0 in different sections was greater than the Fr of sand-covered slopes. Under the same section, the Fr had no obvious change law with the increase of sand thickness (Table 6). Under the same sand thickness, the Fr of the US was larger than the Fr of the FS.

3.2.4. Spatiotemporal Variations of the Darcy-Weisbach Roughness Coefficient

As shown in Figure 7, the f has volatility, but generally increases gradually with time. The range of the f for the US and the FS is 0.06~0.56 and 0.08~1.81, respectively. For the US, there was little change in f of U0, and the change of the f of sand-covered slope with time showed strong fluctuation. This phenomenon may be caused by sediment pick-up. This occurred due to the back water and increased resistance, and forced the f value to increase during the experiment. For the FS, the f of F0 increased sharply with time, and the change of the f of F3 with time was relatively gentle. In the case of frozen soil, the average value of the f decreased with the increasing sand thickness.
In space, the f of the US and FS and the distance from the top of the slope can be represented by a linear function (US: R2 = 0.719, p > 0.05; FS: R2 = 0.972, p < 0.05). Under different treatments, the f decreased as the slope length increased. Given the same thickness of sand cover and the same section, f was greater in FS than that in US (Table 7).

3.3. Quantification of Hydrodynamic Parameters of Slope Erosion under US and FS

3.3.1. Relationship between Flow Velocity and Hydraulic Parameters under US and FS

The flow velocity is one of the basic factors that affect the hydraulic parameters such as Reynolds number (Re), Froude number (Fr) and Darcy-Weisbach roughness coefficient (f). Figure 8 showed the relationship between flow velocity and hydraulic parameters. The Re decreased with increasing V in both the US (R2 = 0.797) and the FS (R2 = 0.871), while the Fr increased with increasing V in the US (R2 = 0.913) and the FS (R2 = 0.977). The f decreased with increasing V on the US (R2 = 0.857) and the FS (R2 = 0.946). Comparing the fitting coefficients of the various relations, it can be found that due to the influence of soil freezed, the flow velocity has less influence on the hydrodynamic parameters. All the determination coefficients (R2) were high (from 0.797 to 0.977), and there was a significant linear relationship between flow velocity and hydraulic parameters (p < 0.01). This may be due to the fact that hydrodynamic parameters are also affected by flow depth and sediment concentration.

3.3.2. Interrelations of Flow Velocity and Hydraulic Parameters with Runoff Rate Response under US and FS

Table 8 showed the correlation between runoff rate with flow velocity (V) and hydraulics parameters (Re, Fr, and f) under different treatments. It can be seen from the equation that the above flow velocity and hydraulic parameters can be used to describe the runoff process under experimental conditions to a certain extent (Table 8). In terms of fitting effect for the runoff process of the US, the test parameters can be arranged in the order of V > Fr > Re > f. For the runoff process on the FS, the test parameters can be arranged in the order of Fr > Re > V > f. The runoff rate, flow velocity, and hydraulics parameters have a significant linear relationship (p < 0.01), and R2 is above 65%. By fitting the data of flow velocity, hydraulic parameters, and sediment yield rate, it is found that although there is a certain relationship between them, this relationship is not significant (p > 0.05).

4. Discussion

4.1. Effects of Slope, Sand Cover, and Soil Freezing on Soil Erosion

For the US, the initial runoff time of the sand-covered slope increased, and the effect was clearer with the increasing sand thickness (Table 2). This is consistent with the research results of Zhang et al. and Tang et al. [8,13]. This relationship occurred due to the high porosity of aeolian sand soil [6]. The greater the sand thickness, the greater the water storage effect, and ultimately the initial runoff time increases greatly. For the FS, the initial runoff time was significantly reduced under different sand thicknesses because the water present in the soil surface layer and the water in the soil pores condense to form an “ice cap”. In early stages, the “ice cap” hindered the inflow and infiltration, resulting in a significant reduction in the initial runoff time [32,33].
The total runoff and total sediment yield on the slope are related to the degree of erosion on the slope during soil erosion [7]. We found the total runoff under different treatments was 1.02 to 1.28 times than that of U0, and the total sediment yield under different treatments was 1.97 to 10.94 times than that of U0 (Table 2). Frozen soil and sand cover on the slope both lead to changes in the total runoff and total sediment yield. The reasons are as follows: (1) When the temperature drops below 0 °C, the water stored in the sand layer and the water in the soil pores freeze into ice, and the volume expands, which reduces soil stability. The bottom layer of the water-tight layer forms an impervious layer, the runoff flows along the contact surface, the friction between the water flow and the slope surface is reduced, and a small ditch is quickly formed, which increases the amount of erosion and eventually leads to increased soil erosion [34,35,36,37,38]; (2) frozen soil greatly shortens the initial runoff time and the appearance time of rill, which makes it easier to generate runoff on the slope surface, and also more likely to generate a fine ditch, leading to increased erosion [39,40,41,42]. The upper-most sand layer can prolong the initial runoff time and store more water. When the slope begins to produce water, the stored water will drain along with the runoff and carry more sediment, increasing sediment yield [4,5,6]. For the US and FS, the sediment yield of sand covered slope is 3.37~4.96 and 1.06~1.13 times of U0 and F0, respectively. In this study, the sediment yield under different treatments increased with the increasing runoff. There was a linear relationship between the cumulative runoff and cumulative sediment yield under different treatments (Table 3). This study further supported the previously reported relationships [43].

4.2. Effects of Slope Sand Covered and Soil Freezed on Hydraulic Parameters

Runoff on the slope is the driving force of soil erosion [44]. The movement of sediment particles will be affected by the runoff [45]. Therefore, there is a close relationship between slope runoff, sediment movement, and hydraulic parameters. In this study, the underlying condition is the main factor affecting hydraulic properties. The hydraulic characteristics of the slope are mainly affected by various factors such as sand cover and soil frozen [15,46]. Sand cover on the slope changed the infiltration capacity of the soil, which in turn changed the runoff of the slope. Therefore, the topography of the sand cover slope changed greatly during the experiment (Figure 2). Compared with the rill formed by U0, the sand-covered slope in the US formed wide and deep due to runoff erosion, and the runoff depth is increased, resulting in changes in hydraulic characteristics. The research results of Tang et al. [15] showed that the flow pattern of water greatly influences the erosion of sand-covered slopes. This study showed that flow velocity and hydraulics parameters can describe the runoff process under different treatments. Among them, R2 of V, Re, and Fr all reached more than 70% (Table 8). Since the data of flow velocity can be obtained directly by the experiment, the flow velocity can be used to better describe the runoff process of the US and FS. However, some seemingly effective hydraulic parameters cannot explain the process of sediment yield on the FS. Frozen soil condenses the water in the soil surface layer and pores into an “ice cap” that hinders runoff infiltration. The infiltration capacity of the slope is reduced, the runoff on the slope is larger, and the erosion is greater [47,48,49]. Therefore, during the experiment, narrow and deep ditches were quickly formed on the FS, resulting in changes in hydraulic characteristics. In addition, the limited observation technology during the erosion process leads to poor fitting of hydraulic parameters and sediment yield on the FS. Despite these shortcomings, the results of this study can still provide a reference for the establishment of a model of erosion on sand-covered loess slopes during thawing.

4.3. Implications for the Relationship between Hydraulic Parameters and Slope Erosion

Global climate change will cause the local permafrost area to melt in advance in the seasonal freeze-thaw area, thus changing the erosion situation in this area. Therefore, in the past few decades, soil erosion resulting from climate warming has captured strong attention in cold regions [49,50,51]. During the thawing period of winter and spring, the soil erosion in the wind water erosion crisscross zone of the Loess Plateau is usually the result of the combined action of water erosion, wind erosion, and freeze-thaw erosion. However, the problem of soil erosion caused by the combined action is far more than the harm of single action erosion itself. The superposition of different types of erosion has led to huge changes in soil erosion [30,52,53,54,55]. In future research, the analysis and quantitative description of composite erosion and single erosion should be encouraged. While it is necessary to measure the impact of each type of erosion on total erosion, it is also necessary to analyze the relationship between hydraulic parameters, runoff, and sediment yield. According to the different effectiveness of hydraulic parameters, selecting the appropriate hydraulic parameters to establish an evaluation model is key [56,57]. This information will help us clarify the feedback relationship between the soil erosion process and the hydrodynamic process from a scientific perspective. At the same time, erosion changes in local areas can also be predicted more comprehensively.

5. Conclusions

Flume tests were performed to study the mechanism of hydrodynamics erosion on the steep sand-covered Loess slopes during the thawing period to improve our understanding of the mechanisms of runoff and sediment yield and to establish a model of soil erosion during thawing. The results showed that the initial runoff time increases with the increase of sand thickness. Under the same sand thickness, the initial runoff time of FS is significantly shorter than the US. For the US, the total sediment yield of different sand thicknesses was significantly higher than that of U0. The cumulative runoff and sediment yield of different treatments can be expressed as a function of M = CQ + D. During the entire experiment, the flow velocity in the US and FS treatments generally showed a downward trend. The distance between the hydraulic parameters of US and FS and the top of the slope can be expressed as a linear function. The main flow pattern of runoff was composed of laminar flow and supercritical flow. Linear equations can be used to describe the relationship between flow velocity and the main hydraulic parameters including Reynolds number, Froude number, and Darcy-Weisbach roughness coefficient. Different hydrodynamic parameters show varying degrees of effectiveness in describing slope erosion processes. Flow velocity is the best hydraulic parameter to simulate the trend of slope erosion process.

Author Contributions

Y.S. and P.L. conceived the main idea of the paper. Z.R., L.X., T.W., and Y.Z. designed and performed the experiment. Y.S. wrote the manuscript and all authors contributed in improving the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded and supported by the National Natural Science Foundation of China (grant no. 51779204), the National Key Research and Development Program of China (grant no. 2017YFC0504501), and the Shaanxi Province Innovation Talent Promotion Plan Project Technology Innovation Team (grant no. 2018TD-037).

Acknowledgments

We thank the reviewers for their useful comments and suggestions. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
USunfrozen slope
U0 (1,2,3)unfrozen slope, the thickness of sand covering is 0 (1,2,3) cm
FSfrozen slope
F0 (1,2,3)frozen slope, the thickness of sand covering is 0 (1,2,3) cm

References

  1. Cao, L.X.; Zhang, K.L.; Dai, H.L.; Liang, Y. Modeling interrill erosion on unpaved roads in the loess plateau of China. Land Degrad. Dev. 2015, 26, 825–832. [Google Scholar] [CrossRef]
  2. Wang, G.Y.; Joho, I.; Yang, Y.S.; Chen, S.M.; Judi, K.; Xie, J.S.; Lin, W.L. Extent of soil erosion and surface runoff associated with large-scale infrastructure development in Fujian Province, China. Catena 2012, 89, 22–30. [Google Scholar] [CrossRef]
  3. Cheng, B.; Lv, Y.; Zhan, Y.; Su, D. Constructing China’s roads as works of art: A case study of “Esthetic Greenway” construction in the Shennongjia Region of China. Land Degrad. Dev. 2013, 26, 324–330. [Google Scholar] [CrossRef]
  4. Xu, J.X. The wind-water two-phase erosion and sediment-producing processes in the middle Yellow River basin, China. Sci. China 2000, 43, 176–186. [Google Scholar] [CrossRef]
  5. Xu, J.X.; Yang, J.S.; Yan, Y.X. Erosion and sediment yields as influenced by coupled eolian and fluvial processes: The Yellow River, China. Geomorphology 2006, 73, 1–15. [Google Scholar] [CrossRef]
  6. Xie, L.Y.; Bai, Y.J.; Zhang, F.B.; Yang, M.Y.; Li, Z.B. Effects of thickness and particle size composition of overlying sand laye on runoff and sediment yield on sand-covered loess slopes. Acta Pedol. Sin. 2017, 54, 60–72. [Google Scholar]
  7. Xu, G.C.; Tang, S.S.; Lu, K.X.; Li, P.; Li, Z.B.; Gao, H.D.; Zhao, B.H. Runoff and sediment yield under simulated rainfall on sand-covered slopes in a region subject to wind-water erosion. Environ. Earth Sci. 2015, 74, 2523–2530. [Google Scholar] [CrossRef]
  8. Zhang, F.B.; Bai, Y.J.; Xie, L.Y.; Yang, M.Y.; Li, Z.B.; Wu, X.R. Runoff and soil loss characteristics on loess slopes covered with aeolian sand layers of different thicknesses under simulated rainfall. J. Hydrol. 2017, 549, 244–251. [Google Scholar] [CrossRef]
  9. Zhang, F.B.; Yang, M.Y.; Li, B.B.; Li, Z.B.; Shi, W.Y. Effects of slope gradient on hydro-erosional processes on an aeolian sand-covered loess slope under simulated rainfall. J. Hydrol. 2017, 553, 447–456. [Google Scholar] [CrossRef]
  10. Zhang, X.; Li, Z.B.; Li, P.; Tang, S.S.; Wang, T.; Zhang, H. Influences of sand cover on erosion processes of loess slopes based on rainfall simulation experiments. J. Arid. Land 2018, 10, 39–52. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, L.P.; Tang, K.L.; Zhang, C.P. Research on soil wind erosion laws in loess hilly-gully region covered by Sheet Sand. J. Soil Water Conserv. 1999, 5, 40–45. [Google Scholar]
  12. Wu, X.R.; Zhang, F.B.; Wang, Z.L. Variation of sand and loess properties of binary structure profile in Hilly Region covered by sand of the Loess Plateau. J. Soil Water Conserv. 2014, 28, 190–193. [Google Scholar]
  13. Tang, S.S.; Li, Z.B.; Ren, Z.P.; Yao, J.W.; Tang, H. Experimental study on the process of runoff and sediment yield on sand-covered slope. J. Soil Water Conserv. 2015, 29, 25–28. [Google Scholar]
  14. Tang, S.S.; Li, Z.B.; Li, C.; Zhao, B.H. Runoff and Sediment yield process on sand covered slope under simulated rainfall. J. Northwest A F Univ. Nat. Sci. Ed. 2016, 44, 139–146. [Google Scholar]
  15. Tang, S.S.; Li, Z.B.; Lu, K.X.; Liu, Y.; Su, Y.Y.; Ma, Y.Y. Relationship between hydrodynamic parameters and runoff and sediment yield on sand-covered slope in rainfall simulation study. Trans. Chin. Soc. Agric. Eng. 2017, 33, 136–143. [Google Scholar]
  16. Zhao, B.H.; Li, Z.B.; Li, P.; Xu, G.C.; Gao, H.D.; Cheng, Y.T.; Chang, E.H.; Yuan, L.; Zhang, Y.; Feng, Z.H. Spatial distribution of soil organic carbon and its influencing factors under the condition of ecological construction in a hilly-gully watershed of the Loess Plateau, China. Geoderma 2017, 296, 10–17. [Google Scholar] [CrossRef]
  17. Shi, P.; Qin, Y.; Liu, Q.; Zhu, T.T.; Li, Z.B.; Li, P.; Ren, Z.P.; Liu, Y.; Wang, F.C. Soil respiration and response of carbon source changes to vegetation restoration in the Loess Plateau, China. Sci. Total Environ. 2019, 707, 135507. [Google Scholar] [CrossRef]
  18. Shi, P.; Feng, Z.H.; Gao, H.D.; Li, P.; Xiao, L. Has “Grain for Green” threaten food security on the Loess Plateau of China? Ecosyst. Health Sustain. 2020, 6, 1709560. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, Y.; Li, P.; Liu, X.J.; Xiao, L.; Shi, P.; Zhao, B.H. Effects of farmland conversion on the stoichiometry of carbon, nitrogen, and phosphorus in soil aggregates on the Loess Plateau of China. Geoderma 2019, 351, 188–196. [Google Scholar] [CrossRef]
  20. Wang, T.; Li, P.; Li, Z.B.; Hou, J.M.; Xiao, L.; Ren, Z.P.; Xu, G.C.; Yu, K.X.; Su, Y.Y. The effects of freeze–thaw process on soil water migration in dam and slope farmland on the Loess Plateau, China. Sci. Total Environ. 2019, 666, 721–730. [Google Scholar] [CrossRef]
  21. Wang, T.; Li, P.; Liu, Y.; Hou, J.M.; Li, Z.B.; Ren, Z.P.; Cheng, S.D.; Zhao, J.H.; Hinkelmann, R. Experimental investigation of freeze-thaw meltwater compound erosion and runoff energy consumption on loessal slopes. Catena 2020, 185, 104310. [Google Scholar] [CrossRef]
  22. Kirkby, M.J. Modeling Water Erosion Processes. In Soil Erosion; Kirkby, M.J., Morgan, R.P.C., Eds.; Wiley: Chichester, UK, 1980; pp. 183–196. [Google Scholar]
  23. Zuzel, J.; Pikul, J. Effects of straw mulch on runoff and erosion from small agricultural plots in northeastern oregon. Soil Sci. 1993, 156, 111–117. [Google Scholar] [CrossRef]
  24. Edwards, L.M.; Burney, J.R. The effect of antecedent freeze-thaw frequency on runoff and soil loss from frozen soil with and without subsoil compaction and ground cover. Can. J. Soil Sci. 1989, 69, 799–811. [Google Scholar] [CrossRef]
  25. Fan, H.; Liu, Y.; Xu, X.; Wu, M.; Zhou, L. Simulation of rill erosion in black soil and albic soil during the snowmelt period. Acta Agric. Scand. Sect. B Soil Plant. Sci. 2017, 67, 510–517. [Google Scholar] [CrossRef]
  26. Ting, J.M.; Torrence, M.R.; Ladd, C.C. Mechanisms of strength for frozen sand. J. Geotech. Eng. 1983, 109, 1286–1302. [Google Scholar] [CrossRef]
  27. Sharratt, B.S.; Lindstrom, M.J.; Benoit, G.R.; Young, R.A.; Wilts, A. Runoff and soil erosion during spring thaw in the northern US Corn Belt. J. Soil Water Conserv. 2000, 55, 487–494. [Google Scholar]
  28. Henry, H.A.L. Soil freeze—Thaw cycle experiments: Trends, methodological weaknesses and suggested improvements. Soil Biol. Biochem. 2007, 39, 977–986. [Google Scholar] [CrossRef]
  29. Pawluk, S. Freeze-thaw effects on granular structure reorganization for soil materials of varying texture and moisture content. Can. J. Soil Sci. 1988, 68, 485–494. [Google Scholar] [CrossRef]
  30. Shi, Z.H.; Fang, N.F.; Wu, F.Z. Soil erosion processes and sediment sorting associated with transport mechanisms on steep slopes. J. Hydrol. Amst. 2012, 454, 123–130. [Google Scholar] [CrossRef]
  31. Zhang, L.T.; Gao, Z.L.; Yang, S.W.; Yang, S.W.; Li, Y.H.; Tian, H.W. Dynamic processes of soil erosion by runoff on engineered landforms derived from expressway construction: A case study of typical steep spoil heap. Catena 2015, 128, 108–121. [Google Scholar] [CrossRef]
  32. Sharratt, B.S.; Lindstrom, M.J. Laboratory simulation of erosion from a partially frozen soil. In Soil Erosion; American Society of Agricultural and Biological Engineers: Honolulu, HI, USA, 2001; pp. 159–162. [Google Scholar]
  33. Pikul, J.L.; Aase, J.K. Fall contour ripping increases water infiltration into frozen soil. Soil Sci. Soc. Am. J. 1998, 62, 1017. [Google Scholar] [CrossRef] [Green Version]
  34. Ma, Q.; Zhang, K.; Jabro, J.D.; Ren, L.; Liu, H. Freeze–thaw cycles effects on soil physical properties under different degraded conditions in Northeast China. Environ. Earth Sci. 2019, 78, 321. [Google Scholar] [CrossRef]
  35. Kimaro, D.N.; Poesen, J.; Msanya, B.M.; Deckers, J.A. Magnitude of soil erosion on the northern slope of the Uluguru Mountains, Tanzania: Interrill and rill erosion. Catena 2008, 75, 38–44. [Google Scholar] [CrossRef]
  36. Wischmeier, W.H.; Smith, D.D. Predicting rainfall erosion losses: A guide to conservation planning. In Agriculture Handbook (USA); Department of Agriculture, Science and Education Administration: Washington, DC, USA, 1978; Volume 537, pp. 5–8. [Google Scholar]
  37. Xiao, L.; Zhang, Y.; Li, P.; Xu, G.C.; Shi, P.; Zhang, Y. Effects of freeze-thaw cycles on aggregate-associated organic carbon and glomalin-related soil protein in natural-succession grassland and Chinese pine forest on the Loess Plateau. Geoderma 2019, 334, 1–8. [Google Scholar] [CrossRef]
  38. Xiao, L.; Yao, K.H.; Li, P.; Liu, Y.; Zhang, Y. Effects of freeze-thaw cycles and initial soil moisture content on soil aggregate stability in natural grassland and Chinese pine forest on the Loess Plateau of China. J. Soils Sediments 2019, 20, 1222–1230. [Google Scholar] [CrossRef]
  39. Øygarden, L. Rill and gully development during extreme winter runoff event in Norway. Catena 2003, 50, 217–242. [Google Scholar] [CrossRef]
  40. Saxton, K.E.; Mccool, D.K.; Papendick, R.I. Slot mulch for runoff and erosion control. J. Soil Water Conserv. 1981, 36, 44–47. [Google Scholar]
  41. Sun, B.Y.; Xiao, J.B.; Li, Z.B.; Ma, B.; Zhang, L.T.; Huang, Y.L.; Bai, L.F. An analysis of soil detachment capacity under freeze-thaw conditions using the Taguchi method. Catena 2018, 162, 100–107. [Google Scholar] [CrossRef]
  42. Ferrick, M.G.; Gatto, L.W. Quantifying the effect of a freeze–thaw cycle on soil erosion: Laboratory experiments. Earth Surf. Process. Landf. 2005, 30, 1305–1326. [Google Scholar] [CrossRef]
  43. Wang, T.; Li, P.; Ren, Z.P.; Xu, G.C.; Li, Z.B.; Yang, Y.Y.; Yao, J.W. Effects of freeze-thaw on soil erosion processes and sediment selectivity under simulated rainfall. J. Arid Land 2017, 9, 234–243. [Google Scholar] [CrossRef] [Green Version]
  44. Asadi, H.; Ghadiri, H.; Rose, C.W. An investigation of flow-driven soil erosion processes at low streampowers. J. Hydrol. 2007, 342, 134–142. [Google Scholar] [CrossRef] [Green Version]
  45. Asadi, H.; Moussavi, A.; Ghadiri, H. Flow-driven soil erosion processes and the size selectivity of sediment. J. Hydrol. Amst. 2011, 406, 73–81. [Google Scholar] [CrossRef] [Green Version]
  46. Bao, Y.X.; Wang, X.; Zhou, L.L.; Chen, Z.Q.; Zhang, K. A study on the temporal and spatial evolutionary processes of the dynamic responses of hydrodynamic erosion parameters on freeze-thaw slopes. J. Soil Water Conserv. 2017, 31, 103–110. [Google Scholar]
  47. Layton, J.B.; Skidmore, E.L.; Thompson, C.A. Winter-associated changes in dry-soil aggregation as influenced by management. Soil Sci. Soc. Am. J. 1993, 57, 1568. [Google Scholar] [CrossRef]
  48. Bullock, M.S.; Nelson, S.D.; Kemper, W.D. Soil Cohesion as affected by freezing, water content, time and tillage. Soil Sci. Soc. Am. J. 1988, 52, 770. [Google Scholar] [CrossRef] [Green Version]
  49. Bochove, E.V.; Danielle, P.; Pelletier, F. Effects of freeze-thaw and soil structure on nitrous oxide produced in a clay Soil. Soil Sci. Soc. Am. J. 2000, 64, 1638–1643. [Google Scholar] [CrossRef]
  50. Dagesse, D.F. Freezing-induced bulk soil volume changes. Can. J. Soil Sci. 2010, 90, 389–401. [Google Scholar] [CrossRef] [Green Version]
  51. Sahin, U.; Anapali, O. Short communication: The effect of freeze-thaw cycles on soil aggregate stability in different salinity and sodicity conditions. Span. J. Agric. Res. 2007, 5, 431–434. [Google Scholar] [CrossRef]
  52. Cheng, Y.T.; Li, P.; Xu, G.C. The effect of soil water content and erodibility on losses of available nitrogen and phosphorus in simulated freeze-thaw conditions. Catena 2018, 166, 21–33. [Google Scholar] [CrossRef]
  53. Wang, S.J. Characteristics of freeze and thaw weathering and its contribution to sediment yield in middle Yellow River basin. Bull. Soil Water Conserv. 2004, 24, 1–5. [Google Scholar]
  54. Zhang, R.F.; Wang, X.; Fan, H.M.; Zhou, L.L.; Wu, M.; Liu, Y.H. Study on the regionalization of freeze-thaw zones in China and the erosion characteristics. Sci. Soil Water Conserv. 2009, 7, 24–28. [Google Scholar]
  55. Wang, F.; Fan, H.M.; Guo, C.J.; Zhou, L.L.; Wu, M.; Liu, Y.H.; Chen, Y.B. Comparison and analysis of climate and environment variation in two main freeze-thaw erosion regions in China. Ecol. Environ. 2008, 1, 173–177. [Google Scholar]
  56. Bryan, R.B. Soil erodibility and processes of water erosion on hillslope. Geomorphology 2000, 32, 385–415. [Google Scholar] [CrossRef]
  57. Knapen, A.; Poesen, J.; Govers, G.; Gyssels, G.; Nachtergaele, J. Resistance of soils to concentrated flow erosion: A review. Earth Sci. Rev. 2007, 80, 75–109. [Google Scholar] [CrossRef]
Figure 1. The geographical location of tested soil.
Figure 1. The geographical location of tested soil.
Water 12 02461 g001
Figure 2. The scour experimental system.
Figure 2. The scour experimental system.
Water 12 02461 g002
Figure 3. Eroded topography under different treatments. Note: S1 (2, 3, and 4), Section 1 (2, 3, and 4).
Figure 3. Eroded topography under different treatments. Note: S1 (2, 3, and 4), Section 1 (2, 3, and 4).
Water 12 02461 g003
Figure 4. Variations in flow velocity under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Figure 4. Variations in flow velocity under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Water 12 02461 g004
Figure 5. Variations in Reynolds number under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Figure 5. Variations in Reynolds number under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Water 12 02461 g005
Figure 6. Variations in the Froude number under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Figure 6. Variations in the Froude number under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Water 12 02461 g006
Figure 7. Variations in the Darcy-Weisbach roughness coefficient under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Figure 7. Variations in the Darcy-Weisbach roughness coefficient under different treatments over time. (a) Unfrozen slope, (b) frozen slope.
Water 12 02461 g007
Figure 8. Relationship between flow velocity and hydraulic parameters. (a) Re with V, (b) Fr with V, (c) f with V.
Figure 8. Relationship between flow velocity and hydraulic parameters. (a) Re with V, (b) Fr with V, (c) f with V.
Water 12 02461 g008
Table 1. Design table of scour experimental.
Table 1. Design table of scour experimental.
State of SlopeTreatmentDepth of Sand
(cm)
Flow Rate
(L/min)
Slope
(°)
Initial Soil Moisture Content (%)Time
(min)
Unfrozen slopes
(US)
U001121515
U111121515
U221121515
U331121515
Frozen slopes
(FS)
F001121515
F111121515
F221121515
F331121515
Table 2. Initial runoff time, runoff and sediment yield under different treatments.
Table 2. Initial runoff time, runoff and sediment yield under different treatments.
TreatmentInitial Runoff Time/(s)RunoffSediment Yield
Total Runoff/LCV/(%)Total Sediment Yield/kgCV/(%)
U039.389.6422.140.9327.80
U1138.0111.0114.053.1331.67
U2186.3611.315.514.0351.16
U3250.5810.6920.674.5963.3
F024.4513.822.487.763.33
F131.9114.186.178.2326.34
F234.0213.4613.8610.0426.22
F332.4112.2411.9310.1735.22
Table 3. Cumulative runoff and cumulative sediment yield fitted equation.
Table 3. Cumulative runoff and cumulative sediment yield fitted equation.
TreatmentFitted EquationMC
U0M = 0.093 Q + 0.029R2 = 0.997, p < 0.0010.31
U1M = 0.283 Q − 0.016R2 = 0.990, p < 0.001
U2M = 0.397 Q − 0.035R2 = 0.972, p < 0.001
U3M = 0.465 Q + 0.263R2 = 0.930, p < 0.001
F0M = 0.565 Q + 0.001R2 = 0.999, p < 0.0010.67
F1M = 0.593 Q − 0.361R2 = 0.987, p < 0.001
F2M = 0.740 Q + 0.566R2 = 0.986, p < 0.001
F3M = 0.778 Q + 1.003R2 = 0.993, p < 0.001
Table 4. Variations of flow velocity with the distance from the top of the slope (m/s).
Table 4. Variations of flow velocity with the distance from the top of the slope (m/s).
TreatmentDistance from the Top of the Slope E/mFitted Equation
0.511.52
U00.260.290.290.32VU = 0.071E + 0.212
R2 = 0.893, p > 0.05
U10.210.270.330.29
U20.220.310.310.35
U30.230.330.370.39
F00.190.200.260.33VF = 0.091E + 0.146
R2 = 0.952, p < 0.05
F10.180.210.310.37
F20.190.220.290.28
F30.240.250.330.32
Note: The parameter “E” means the distance from the top of the slope, the same to below.
Table 5. Variations of Reynolds number with the distance from the top of the slope.
Table 5. Variations of Reynolds number with the distance from the top of the slope.
TreatmentDistance from the Top of the Slope E/mFitted Equation
0.511.52
U0283.55194.56172.46181.68ReU = −135.01E + 484.59
R2 = 0.712, p > 0.05
U1570.51380.66433.42309.84
U2636.61375.33371.24392.56
U3587.71463.96388.73410.48
F0518.21499.82444.29430.38ReF = −57.09E + 489.25
R2 = 0.998, p < 0.01
F1433.68410.35340.95350.17
F2396.99499.04399.28347.70
F3488.91328.52326.88371.09
Table 6. Variations of Froude number with the distance from the top of the slope.
Table 6. Variations of Froude number with the distance from the top of the slope.
TreatmentDistance from the Top of the Slope E/mFitted Equation
0.511.52
U03.374.474.675.32FrU = 1.488E + 1.656
R2 = 0.899, p > 0.05
U11.862.893.673.66
U21.403.743.764.34
U31.762.984.074.29
F01.351.451.923.18FrF = 1.362E + 0.656
R2 = 0.956, p < 0.05
F11.251.683.394.03
F21.371.622.722.86
F31.742.383.563.24
Table 7. Variations of Darcy-Weisbach roughness coefficient with the distance from the top of the slope.
Table 7. Variations of Darcy-Weisbach roughness coefficient with the distance from the top of the slope.
TreatmentDistance from the Top of the Slope E/mFitted Equation
0.511.52
U00.210.090.080.06fU = −0.348E + 0.694
R2 = 0.719, p > 0.05
U10.720.250.130.13
U21.040.140.130.09
U30.650.240.110.09
F01.711.570.850.26fF = −0.754E + 1.695
R2 = 0.972, p < 0.05
F11.570.950.330.16
F21.181.090.450.36
F30.780.430.170.19
Table 8. Correlation between runoff rate with flow velocity and hydraulic parameters.
Table 8. Correlation between runoff rate with flow velocity and hydraulic parameters.
Runoff/LHydraulic Parameters
V/(m/s)ReFrf
USR = −6.394V + 2.559R = 0.0033Re − 0.443R = −0.278Fr + 1.788R = 1.207f + 0.445
R2 = 0.813, p < 0.01R2 = 0.725, p < 0.01R2= 0.767, p < 0.01R2 = 0.655, p < 0.01
FSR = −1.474V + 1.273R = 0.0013Re + 0.437R = −0.084Fr + 1.141R = 0.167f + 0.791
R2 = 0.787, p < 0.01R2 = 0.85, p < 0.01R2 = 0.866, p < 0.01R2 = 0.668, p < 0.01
Note: R was the runoff (L).

Share and Cite

MDPI and ACS Style

Su, Y.; Li, P.; Ren, Z.; Xiao, L.; Wang, T.; Zhang, Y. Slope Erosion and Hydraulics during Thawing of the Sand-Covered Loess Plateau. Water 2020, 12, 2461. https://doi.org/10.3390/w12092461

AMA Style

Su Y, Li P, Ren Z, Xiao L, Wang T, Zhang Y. Slope Erosion and Hydraulics during Thawing of the Sand-Covered Loess Plateau. Water. 2020; 12(9):2461. https://doi.org/10.3390/w12092461

Chicago/Turabian Style

Su, Yuanyi, Peng Li, Zongping Ren, Lie Xiao, Tian Wang, and Yi Zhang. 2020. "Slope Erosion and Hydraulics during Thawing of the Sand-Covered Loess Plateau" Water 12, no. 9: 2461. https://doi.org/10.3390/w12092461

APA Style

Su, Y., Li, P., Ren, Z., Xiao, L., Wang, T., & Zhang, Y. (2020). Slope Erosion and Hydraulics during Thawing of the Sand-Covered Loess Plateau. Water, 12(9), 2461. https://doi.org/10.3390/w12092461

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