Next Article in Journal
Train-Induced Building Vibration and Radiated Noise by Considering Soil Properties
Next Article in Special Issue
Rural Sustainable Environmental Management
Previous Article in Journal
Sustainable Business Model Based on Digital Twin Platform Network: The Inspiration from Haier’s Case Study in China
Previous Article in Special Issue
Safety Guarantee System of Drinking Water Source in Three Gorges Reservoir Area and Its Application in Huangjuedu Drinking Water Source Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Soil and Water Conservation Practices in the Loess Hilly Region Using a Coupled Rainfall-Runoff-Erosion Model

1
Department of Building, Civil and Environmental Engineering, Concordia University, Montreal, QC H3G 1M8, Canada
2
Department of Chemical and Materials Engineering, Concordia University, Montreal, QC H3G 1M8, Canada
3
Institute for Energy, Environment and Sustainable Communities, University of Regina, Regina, SK S4S 0A2, Canada
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(3), 934; https://doi.org/10.3390/su12030934
Submission received: 6 January 2020 / Revised: 21 January 2020 / Accepted: 22 January 2020 / Published: 27 January 2020
(This article belongs to the Special Issue Rural Sustainable Environmental Management)

Abstract

:
Soil and water conservation practices (SWCPs) are widely used to control soil and water loss. Quantifying the effect of SWCPs and climate change on soil and water erosion is important for regional environmental management. In this study, the Soil Conservation Service Curve Number (SCS-CN) and the Modified Universal Soil Loss Equation (MUSLE) were employed to investigate the patterns of surface runoff and soil erosion with different SWCPs in the hilly region on the Loess Plateau of China. The impact of climate change under RCP4.5 and RCP8.5 emission scenarios was considered from 2020 to 2050. Surface runoff grew with the increased rainfall and rainfall erosivity, while soil erosion presented large variations between years due to uneven distribution of rainfall and rainfall erosivity under two scenarios. Different SWCPs significantly reduced surface soil and water loss. Compared with bare slopes, the reduction rates were 15–40% for surface runoff and 35–67% for soil erosion under RCP4.5 and RCP8.5 emission scenarios, respectively. The combination of shrub and horizontal terracing was recommended due to its low water cost for sediment control among seven SWCPs.

1. Introduction

Soil erosion is a global problem of land degradation which may adversely influence natural and human-affected ecosystems [1]. Around 25–40 billion tons of surface soil and 10 million hm2 of farmland worldwide are eroded annually due to soil erosion, causing a direct economic loss of nearly $400 billion USD [2,3,4]. In particular, approximately five billion tons of topsoil and 67 thousand hm2 of cultivated land are lost each year in China as a result of soil erosion [5]. The Loess Plateau is one of the most important agricultural regions in China, with serious soil erosion [6]. The Grain-to-Green project implemented since 1999 aimed to improve soil conservation in this region. The soil erosion in the Loess Plateau, however, was still higher than 3600 t km−2 yr−1, exceeding the negligible erosion rate of 1000 t km−2 yr−1 [7]. The use of fertilizers in the farmlands with soils prone to erosion risk may affect groundwater and surface water quality. For example, the dissolved reactive phosphorus from fertilizers will transport with surface runoff from farmlands, accelerating eutrophication in receiving water bodies [8]. The release of heavy metals such as Zn, Cr, Cd, Hg, and Pb from soil erosion may also pose a risk for human health [9].
Soil erosion can be affected by some key factors such as climate change, land use, and anthropogenic activities [10,11,12]. The high-intensity rainfall events have important effects on soil erosion [13]. The mean surface temperature as well as the intensity and frequency of extreme rainfall events will increase during the twenty-first century in the Loess Plateau [14]. Due to climate change, global soil erosion in the next century is expected to increase by approximately 25–50% [15]. Moreover, soil and water erosion may be influenced by other factors (e.g., soil water content and vegetation cover), which will also be affected by climate change. Field monitoring is easily conducted at small scales (e.g., patches, hillslopes and plots), but it is a challenge to implement over large areas (e.g., watershed and region) considering the required labor and cost [16]. Therefore, evaluations are required to estimate soil erosion at a large scale (e.g., the Loess Plateau) for decision making and regional management considering the future climate change [17]. The Universal Soil Loss Equation (USLE) and its modifications are widely used empirical models for water erosion prediction because they can be applied at large scales with various inputs to consider the future climatic uncertainties [18,19].
Soil and water conservation practices (SWCPs) are one of the most effective measures for erosion control. It can also help adapt to climate change and improve rural development [20]. Over the past 60 years, various SWCPs, including terraces, natural restoration, afforestation, and engineering construction (e.g., check-dams), have been applied in more than 30% of the Loess Plateau [21]. These practices can significantly reduce soil and water loss and improve environmental quality [22,23,24]. For example, the use of a permeable reactive barrier filled with materials such as activated carbon and leonardite can intercept the contaminant transport and even transform the contaminant into low or non-toxic substances [25,26]. Liu et al. [27] summarized 53 types of SWCPs used in China. These SWCPs can significantly reduce soil erosion and runoff, but the efficiency of these SWCPs was also affected by the topography, soil characteristics, climatic regime, etc. [28]. Meanwhile, the regional water scarcity in the Yellow River Basin and the conflict between the upper and downstream of the Yellow River in water usage are serious problems. The average runoff may decrease with the increase of vegetation evapotranspiration due to the application of SWCPs [29]. With the growing water scarcity and increasing competition among water users, more efficient water use has increased in importance for sustainable water management and environmental decision-making [30,31,32]. Therefore, this study aims to develop an approach for the evaluation of SWCPs under future climate scenarios. Improved assessment of SWCPs will be achieved through the use of water cost which can help analyze the balance of soil conservation and runoff reduction, and the incorporation of climate change into soil erosion prediction, which can help understand the performance of SWCPs in a changing environment. Specific goals include evaluating the impact of climate change on soil erosion and surface runoff in a watershed and assessing the water cost for sediment control of seven SWCPs under the future climate scenarios. Climate change, surface runoff, and soil erosion simulated by a coupled rainfall-runoff-erosion model will be investigated to select appropriate SWCPs for the hilly regions on the Loess Plateau. This study regarding the effect of climate change and SWCPs on soil and water loss can be used to support watershed management.

2. Methodology

2.1. Erosion Measurement

The experimental plots locate at the Dongshan watershed (38°02′20″–38°03′48″ N, 111°46′48″–111°48′15″ E) of Shanxi Province in China, which is in the southeast of Loess Plateau. This area has a typical semi-arid continental monsoon climate and a fragile ecological environment, with a mean annual rainfall of 500 mm (420–540 mm), and a mean air temperature of 9.5 °C (8.80 to 10.2 °C). The rainfall is seasonally unsteady, with intensive and short-duration events during the rainy season (June–September), accounting for 75–90% of total annual rainfall [33]. The average effective evaporation is about 2160 mm per year, nearly 5 times more than rainfall [34]. Soils in this region are mainly loose and porous loess, containing 64–73% silt and 17–20% clay, which are vulnerable to erosive activities (e.g., raindrop, runoff, gravity, freeze-thaw, and wind) [33,35]. The watershed, with area around 433 hm2 and elevation from 1135 to 1775 m above sea level, has been affected by human disturbances (e.g., urbanization, deforestation, and mining) and climate change. The reduction of vegetation cover in this region resulted in a gradual increase in soil erosion, flooding and sediment yield. The water-induced erosion models can be modified based on erosion measurements since the soil erosion on the Loess Plateau mainly comes from steep, sparsely vegetated, and highly erosive loess slopes with intensive rainfall and water flow [18,21].
Based on the major SWCPs in the Loess Plateau, eight standard experimental plots with the dimension of 5 m in width and 20 m in length each were built to study the impact of SWCPs on soil and water erosion. They were A (Sea-buckthorn + horizontal terracing), B (grass + horizontal terracing), C (forest/Chinese Pine + fish-scale pitting), D (bare slope), E (forest/Arborvitae + fish-scale pitting), F (forest/Arborvitae + horizontal terracing), G (forest/Apricot + horizontal terracing), and H (forest/Chinese Pine + horizontal terracing). The basic characteristics of these plots are shown in Table 1. These experimental plots were located on the upper and middle positions of the natural semi-shady slope, with a slope of 18°. The soil type at this site was cinnamon soil and the soil thickness was 130 cm. Each plot was bordered with a 35-cm-wide cement wall. A drainage ditch on the upper part and 1.2 m-wide isolation zones on both sides of each plot were designed to prevent the raindrop splashing from the surrounding areas. The lowest part of each plot was built with an outlet (30 × 30 cm) linked to the trough, water pipe, and a cement reservoir (2 × 1.5 × 2 m).
The reduction rates of surface runoff ( f Q , %) and soil erosion ( f A , %) of different SWCPs were calculated as follows [36]:
f Q = Q 0 Q i Q 0 ,
f A = A 0 A i A 0 ,
where Q 0 and A 0 are surface runoff (mm) and soil erosion (t hm−2 yr−1) of the reference plot (D, natural and bare slope), respectively; Q i and A i are surface runoff (mm) and soil erosion (t hm−2 yr−1) of experimental plots with different SWCPs (e.g., grass, fish-scale pitting, horizontal terracing), respectively.
Surface runoff and soil erosion on the slopes can be reduced by SWCPs at the same time, but the reduction effects vary in water and soil erosion among different practices. Water cost for sediment control was proposed by Wang et al. [37] to compare the differences in the decrease of runoff when reducing the same amount of sediment among different practices. It can be calculated by the ratio of reduced runoff and soil erosion (Rrs, mm (t hm−2 yr−1)−1), which is defined as the ratio of the water-reduced runoff to the reduced soil erosion by SWCPs. The formula was as follows:
R r s = Q d A d ,
where Qd and Ad are the reduced surface runoff (mm) and soil erosion (t hm−2 yr−1) of each SWCP, respectively. The SWCPs with a small Rrs have smaller impact on runoff when reducing the same amount of soil erosion.

2.2. Rainfall-Runoff-Erosion Model

2.2.1. Regional Climate Model

In this study, climate simulations conducted by regional climate models (RCMs) were used to dynamically downscale climate predictions of global climate models (GCMs) [38]. The RCM system (RegCM4.5) with the hydrostatic core was used in this study, coupled with the Community Land Model (CLM4.5). The microphysics scheme and Emanuel scheme were applied as the moisture scheme and cumulus convective scheme, respectively. This parameterization scheme combination showed satisfying performance in previous studies. For example, Gao et al. [39] found the combination of the Emanuel scheme and CLM3.5 could better simulate the precipitation over China than other scheme combinations. Aerosol processes and atmosphere-ocean coupling, however, were not considered in this study. Besides RCM and scheme options, GCM, scenario, and boundary conditions may cause uncertainties [14]. Ideally, more reliable and consistent climate projections can be produced through a multi-model ensemble [40]. RCM precipitation post-processed by nonparametric transformations could reduce biases and produce reliable estimates of local-scale climate [41]. However, due to the high requirement of computational resources, it is difficult to address all sources of uncertainty and bias. Thus, RegCM and the abovementioned schemes were selected in this study. The area described in the RegCM4.5 is a region of 110 by 135 grid points with a resolution of 50 km centered at E111°46′ and N38°3′. To predict the future climate, two Representative Concentration Pathways (RCPs) published by the Inter-governmental Panel on Climate Change (IPCC) were selected to drive the simulation with the highest (RCP8.5: 8.5W/m2 by 2100) and moderate (RCP4.5: 4.5W/m2 by 2100) scenarios. The simulations were conducted between 1981 and 2005, with the assumption of accelerated development in the initial five years and exceeding the baseline after 1986. The future climate simulations were run from 2006 to 2050, and the analysis was conducted with the results from 2020 to 2050. The values of monthly rainfall on a 0.5° grid were available for the hilly Loess Plateau between 2001 and 2050. In this study, climate change scenarios were generated using downscaled average monthly rainfall inputs during 1986–2005. All scenarios were run for 30 years, and rainfall changed during the entire period of 2020–2050. The simulations were presented as annual and monthly averages. The simulated rainfall will provide input to the following models.

2.2.2. SCS-CN Model

The Soil Conservation Service curve number (SCS-CN) is an empirical rainfall-runoff model with assumptions and a few inputs to estimate runoff caused by a rainfall event. It is based on the principle of the water balance and two basic assumptions—the proportion of actual runoff to potential maximum runoff is equivalent to the ratio of infiltration to potential maximum interception, and the initial abstraction corresponds to the potential maximum interception. This method was presented as follows [39]:
Q =   { ( P λ S ) 2 P λ S + S P > λ S 0 P < λ S ,
S = 254 × ( 100 C N 1 ) ,
where P is rainfall (mm); Q is surface runoff (mm); λ is the initial abstraction rate, which is assumed to be 0.2 in the SCS-CN method. Various values of initial abstraction rate ranging from 0.0 to 0.3 were reported in previous studies, and it was more reasonable for λ varying between 0.00–0.05 [42,43]. Combined with the sensitivity analysis conducted in the Loess Plateau by Wang [44], λ = 0.00 was selected for this study. S is the maximum retention (mm), which is directly linked to the curve number (CN). CN is a dimensionless variable, which can be determined by two methods: (1) using the CN value-based lookup table based on soil type and land use provided by SCS and (2) calculating the CN value according to the measured data of the experimental plots based on the backward formula from SCS-CN equation (λ = 0.00) as follows:
C N = 100 P 254 × ( P Q 1 ) + 1 .

2.2.3. MULSE Model

The USLE model was widely applied to predict soil erosion by water forces on the Loess Plateau and it has six important aspects including rainfall erosivity (R), soil erodibility (K), slope length (L), slope gradient (S), crop management (C), and erosion control practice (P) [45]. USLE has some modifications and revisions, mainly including the revised USLE (RUSLE), Chinese Soil Loss Equation (CSLE), and the modified USLE (MUSLE) [19,46,47]. Li et al. [18] compared 11 commonly used soil erosion models on the Loess Plateau and found that empirical models produced sound and quick predictions of soil erosion and sediment yield. On the Loess Plateau, soil erosion was dominated by heavy and intensive rains while low-rainfall events were less important for it. Thus, considering the direct impact of surface runoff on soil erosion, the rainfall-runoff erosivity index (QRR), reflecting the erosive force under the combined rainfall and runoff, was introduced into the MUSLE model for Loess Plateau [39]. The runoff-erosion model was formulated as follows:
A = a ( Q R R ) b × K × L S × C × P ,
where A is soil erosion (t hm−2 yr−1) on a storm basis; R represents rainfall erosivity factor (MJ mm hm−2 h−1 yr−1); Q R is the runoff coefficient (%), which is the ratio of future runoff and rainfall predicted by the SCS-CN model and the RegCM4 model, respectively; K is the soil erodibility factor (t hm2 h hm−2 MJ−1 mm−1); L is the slope length factor; S stands for the slope factor; C is the dimensionless vegetation cover factor; and P refers to a dimensionless erosion control practice factor. Recommended values for empirical coefficients a and b in this model on the Loess Plateau were 1.723 and 1.548, respectively [39]. The measured data from experimental plots were used to determine the application conditions and model variables for MUSLE.
Rainfall erosivity factor (R) is one of the crucial factors in the MULSE model and it has a direct link with the potential of soil erosion caused by rainfall [48]. Some algorithms of rainfall erosivity have been studied, primarily including classical algorithms and simple algorithms. Wu et al. [49] reviewed 23 R-factor algorithms and their application scope from 1958 to 2006. The study found that the two simple algorithms—the Wischmeier model and the half-month model—performed better for the hilly Loess Plateau. Although the half-month model was of higher accuracy and more appropriate to be used in the hilly areas, input data for the Wischmeier model were more easily available [49,50]. Thus, monthly rainfall (Pi) and annual rainfall (P) were used to estimate the R factor with the Wischmeier model in this study [45]:
R = 17.02 × i = 1 12 1.735 × 10 ( 1.5 l o g ( P i 2 P ) 0.08188 ) ,
where R is the annual rainfall erosivity factor (MJ mm hm−2 h−1 yr−1); Pi is the average monthly rainfall of i-th month from January to December (mm); P is the annual rainfall (mm).
Soil erodibility factor (K) is the key parameter to assess the soil susceptibility to erosion, reflecting the soil vulnerability to raindrop splash and runoff flow. The widely used K estimators so far include USLE, the Revised Universal Soil Loss Equation (RUSLE), the Erosion Productivity Impact Calculator (EPIC), and the Geometric Mean Diameter based (Dg) model [51]. These models were developed according to American soil-erosion databases, so they might be ineffective for other regions, e.g., the Loess Plateau. Therefore, Zhang et al. [51] combined several methods for K factor estimation based on the data from the Second Nation Soil Survey (China) and provided the optimally estimated K factor (Kopt) combining the revised Knomo, Kepic, and KDg methods. The new estimator performed well in their study and was adopted in the present study as follows:
K o p t = 0.01044 + 1.7530 3 [ ( 0.03336 + 0.7449   K n o m o ) + ( 0.01383 + 0.5158   K e p i c ) + ( 0.00911 + 0.5507   K D g ) ] ,
K n o m o = 0.1317 × 2.1 × 10 4 ( 12 S O M ) M 1.14 + 3.25 ( S 2 ) + 2.5 ( P 3 ) 100 ,   M = SIL × ( SAN + SIL ) ,
K e p i c = 0.1317 × [ 0.2 + 0.3 exp ( 0.0256 S A N ( 1 S I L 100 ) ) ] ( S I L C L A + S I L ) 0.3 ( 1.0 0.25 S O M S O M + exp ( 3.72 2.95 S O M ) ) ( 1.0 0.7 S N 1 S N 1 + exp ( 5.51 + 22.9 S N 1 ) ) ,   SN 1 = 1 S A N 100 ,
K D g = 7.594 { 0.0017 + 0.0494 exp [ 0.5 ( log ( D g ) + 1.675 0.6986 ) 2 ] } ,
where SAN, SIL, and CLA are the fraction of sand (Dg = 0.05-0.1 mm, %), silt (Dg = 0.002-0.05 mm, %) and clay (Dg < 0.002 mm, %) in the US Soil Texture System, respectively; SOM is the soil organic matter (%); S is the soil-texture code used in soil classification; P is the profile-permeability class; Dg is the geometric mean diameter (mm).
Topographical factors, including slope length factor (L) and slope gradient factor (S), are the factors that directly affect soil erosion. In the original erosion model, slope length is regarded as the horizontal distance from the origin of runoff to the end where soil begins to deposit or runoff starts to concentrate in the defined channel with decreased slope gradient [52]. Rainfall redistribution due to interception and evaporation of vegetation canopies can preserve and reduce surface runoff, so the surface runoff is affected by the specific area and land-use changes along the flow paths [50]. Some calculating methods for the L-factor have been applied in different models at different scales and scopes [18,50,53]. Qin et al. [50] reviewed four L-factor calculating methods which were currently widely used for USLE and RUSLE in many countries and proposed a new method to calculate the L-factor, considering the comprehensive effect of upslope vegetation coverage and topography on the downslope surface runoff and soil erosion. However, soil erosion increased linearly with the rise of slope length in a uniform slope with limited length [45]. The formulas for L-factor and S-factor were primarily based on experimental data from natural runoff plots with different slope lengths and slope gradients, respectively, which were as follows [50,54]:
L = ( λ 22.13 ) m     ,   m =   {     0.5             θ > 5 ° 0.4 3 ° < θ 5 ° 0.3 1 ° < θ 3 ° 0.2 θ 1 ° ,
S =   { 21.9 sin ( θ / 180 ) 0.96                 θ 10 ° 16.8 sin ( θ / 180 ) 0.5                 5 ° θ < 10 ° 10.8 sin ( θ / 180 ) + 0.03                 θ < 5 ° ,
where λ and θ are the horizontal slope length (m) and the slope gradient (°), respectively; m is a dimensionless empirical coefficient depending on the slope gradient.
Cover management factor (C) refers to the ratio of soil erosion from a slope with and without cover management [45]. This factor can reflect the importance of land cover, crops and crop management on soil erosion [45,55]. C-factor is closely relevant to vegetation cover which is one sensitive factor influencing soil erosion. Therefore, an equation describing the relationship between C-factor and vegetation coverage based on the observations on runoff plots under natural and simulated rainfall was introduced, which was expressed as follows [55]:
C = { 1 f = 0 0.6508 0.3436 lg f 0 < f 78.3 % 0 f > 78.3 % ,
where f is the vegetation coverage (%).
Erosion control practice factor (P) is the ratio between soil erosion in a slope with and without one or more specific control practices [45]. These practices will probably change the flow path, direction, and amount of surface runoff, resulting in the reduction of soil erosion [56]. SWCPs in China could be classified into the following three major categories: biological, engineering, and tillage practices [27,47]. By far, contour farming and no-tillage were the most widely applied practices among different tillage options; horizontal terracing was the most frequently used engineering practice; and hedgerow or shrub planting was the most usually tested biological practice [27,28]. The efficacy of 15 most frequently tested SWCPs on reducing soil and water erosion was assessed and used to estimate the P-factor [28]. If there is no significant control practice, the P factor is regarded as 1. Based on the above results, the values of the P-factor are listed in Table 1.

2.3. Calibration and Validation

The calibration and validation method was based on the differences between measured and modelled values. The Nash-Sutcliffe efficiency (NSE) and Root Mean Square Error (RMSE) were adopted to assess the performance of SCN-CN and MUSLE as follows [57]:
R M S E = 1 n i = 1 n ( O i P i ) 2 ,
N S E = 1 i = 1 n ( O i P i ) 2 i = 1 n ( O i O ¯ ) 2 ,
where Oi is the i-th observed value; Pi is the i-th predicted value; O ¯ is the average value of the entire observations. RMSE is always positive, and a lower value close to zero means a good fit for the data [58]. NSE ranges from negative infinity to one and the value closer to one denotes a higher model accuracy. NSE > 0.4 is required for this study as the model performance evaluation standard.

3. Results and Discussion

3.1. CN and K in the Rainfall-Runoff-Erosion Model

Runoff occurs when rainfall is greater than soil absorption or infiltration. Rainfall resulting in surface runoff and soil erosion is called erosive rainfall [59]. The observed rainfall and surface runoff of 44 erosive rainfall events between 2014 and 2017 were used as reference data to calculate the CN of the experimental plots with different SWCPs from January to December (Table 2). The CN values averaged 26.6 among eight plots, varying from 19.2 in plots with Chinese Pine and fish-scale pitting to 38.8 in bare slope. They were generally higher in summers or rainy seasons (May–October) with intensive and frequent rainfalls. These results were consistent with those in some previous studies. Soulis et al. [60] applied the SCS-CN method to analyze the direct runoff generation mechanism in a Mediterranean experimental watershed in Greece and found a strong correlation between the CN values obtained from measured runoff and the rainfall depth. The proposed two-CN system well described the CN-rainfall variation observed in natural watersheds, which could be attributed to the spatial variability of soils and land cover in the watershed [61]. Xiao et al. [62] found CN values had large differences, ranging from 15 in rangeland to 100 in bare rock land, for various land uses in the Liudaogou watershed, which was also within the hilly region of the Loess Plateau. The bare slope with large CN value indicated bare slope had more surface runoff and less infiltration [63]. Wang [44] also reported CN values in the Loess Plateau, varying between 8.81 (λ = 0.00) and 75.27 (λ = 0.20), which was closely related to the selected initial abstraction rate. The model with λ = 0.00 had the smallest CN value and the largest potential maximum interception (S) in the watershed because this revision ignored the initial abstraction before runoff and the antecedent soil moisture content.
Table 3 shows the values of soil erodibility factor (Kopt, t hm2 h hm−2 MJ−1 mm−1) of experimental plots with different SWCPs, varying between 0.001 in plot B and 0.012 in Plot E. In previous studies, Zhang et al. [64] calculated a set of erodibility values for main soils on the Loess Plateau based on the data from several field stations and found that the K-factor for loess soils varied from 0.3 to 0.6 in the American system and from 0.008 to 0.04 in the metric system. Zhang et al. [51] also proposed a linear combination of the revised Kepic and KDg methods for soil erodibility estimation based on the data from the Second Nation Soil Survey and the K value of 0.0163 t hm2 h hm−2 MJ−1 mm−1 was obtained for loessial soils with slope gradient of 9%. Wang et al. [65] reported that previous equations overestimated most of the K-factor and the KDg, and Knomo could be the only reliable method to be applied directly to the Loess Plateau. The mean K values were 0.0321 t hm2 h hm−2 MJ−1 mm−1 in Chinese water erosion regions and 0.0384 in Shanxi, which were close to the values of KDg (0.0411) and Knomo (0.0451) but lower than Kepic values of 0.053 t hm2 h hm−2 MJ−1 mm−1 in this study. Soil characteristics, such as soil texture and SOM, had significant impacts on soil erodibility [51,66]. The soil in this study area had low SOM with an average of 1.01%, thus, the values of soil erodibility factor could be overestimated for these plots when Kepic estimator was adopted to estimate K-factor.

3.2. Observed and Predicted Soil and Water Erosion

The calibration and validation results for the surface runoff and soil erosion are presented in Table 4. The minimum RMSE between the observed (Oi) and predicted (Pi) soil erosion was achieved with the empirical coefficients (a and b) of SWCPs in the range of 1.609–3.409 and 1.150–1.631, respectively. The modeling performance of MULSE for soil erosion and SCS-CN for the surface runoff was good, and the NSE values showed a positive efficiency for all experimental plots (NSE > 0.4). The surface runoff was overestimated for the majority of severe runoff events on all experimental plots, and soil erosion was underestimated for smaller runoff values and overestimated for heavy erosive events. Model parameters and input data may lead to uncertainty in the simulation process. Slope gradient, antecedent soil moisture, and rainfall intensity are crucial factors to predict surface runoff and soil erosion [67]. In this study, only monthly rainfall was used as climatic input to the model. Rainfall intensity and rainfall duration greatly influenced the amount of runoff, but they were not included in this SCS-CN model [68]. The rainfall erosivity calculated using the Wischmeier model based on monthly rainfall was higher than that calculated using a half-month model based on daily rainfall. That could be attributed to a certain degree of deviation between high intensity and short duration rainfall [49], which overestimated the soil erosion especially for heavy erosive events. Consequently, the MUSLE model may overestimate the total soil erosion from eight experimental plots in this study. Notably, the gravitational soil erosion was not considered in the MUSLE model of this study because water-induced erosion was dominating in this area.

3.3. Impact of Climate Change on Soil and Water Erosion

The impact of climate change on soil and water erosion is shown in Figure 1 in terms of the relative difference in rainfall, rainfall erosivity factor, surface runoff, and soil erosion between the historical climate (1986–2005) and the two RCP scenarios (RCP4.5 and RCP8.5) during 2020–2050. Projections by RegCM4 showed that the annual average rainfall for the two scenarios during that period were 909.6 and 922.6 mm, respectively, with an increase of 16% and 18% compared to historical climate (785.0 mm), respectively. The result was consistent with the previous studies. Yellow River Basin would be likely to experience rainfall increase in the projected scenarios [14]. Compared with historical climate (1950–1999), GCMs projected a 4–18% and 23–37% increase in annual rainfall for the region during 2010–2039 and 2070–2099, respectively [69]. Rainfall erosivity is influenced by the amount, intensity, and duration of rainfall. Its temporal distributions have a strong relationship with the erosive rainfall. Generally, more rainfall can cause larger rainfall erosivity [49]. The annual averages of the rainfall erosivity for RCP4.5 and RCP8.5 in the years of 2020–2050 were 20,518.72 and 20,606.60 MJ mm hm−2 h−1 yr−1, respectively, with a slight increase of 1% and 7% compared to historical climate (20,226.18 MJ mm hm−2 h−1 yr−1). These values were higher than those in previous studies probably because the impact of high intensity and short duration rainfall on the Wischmeier model was greater than that on the half-month model [49]. Annual surface runoff for RCP4.5 and RCP8.5 emission scenarios showed a similar average, approximately 141 mm, during the period 2020–2050, which was 28% higher than that of the historical climate. However, annual soil erosion was projected with the changes in the range of −96% to 235% and −92% to 448% for the two scenarios compared to the historical climate, respectively. The large variations with years might be due to the large intra-annual variability of rainfall and rainfall erosivity. For example, the rainfall depth (around 725 mm) was approximately equal for 2040, 2048, and 2050 under RCP8.5 scenario, but the rainfall erosivity and soil erosion in 2040 were over twice as much as the others. This was because heavy rainfalls were more concentrated in rainy seasons and 57%, 40%, and 31% of total rainfall depth occurred in July and August in 2040, 2048, and 2050, respectively. The surface runoff and soil erosion in the rainy season (May–October) under future scenarios accounted for nearly 90% and 99% of the total amount in a year, respectively. This was consistent with the results of Azari et al. [70] that nearly 75% of the soil erosion occurred in the rainy season because of heavy rainfall. The frequency and intensity of large storms may increase rainfall erosivity force and cause more surface runoff and soil erosion [69,70]. Water-induced soil erosion mainly resulted from soil particle detachment by the shearing force of raindrops and overland flow, and the movement and transport by rain splash and runoff [71]. Surface runoff and soil erosion generally increased with rainfall, but many other factors such as antecedent rainfall and vegetation coverage also affected them [72]. Antecedent rainfall had a strong impact on the generation of surface runoff and soil erosion on slopes during rainfall events. Surface runoff and soil erosion on slopes with large antecedent rainfall may be higher than that without antecedent rainfall because the antecedent rainfall could recharge and saturate surface fractures which usually reduce surface runoff. Saturation excess runoff generally occurred when the rainfall intensity was higher than soil infiltration rates [73]. Therefore, more heavy rainfall and runoff events in the rainy season will result in more soil and water erosion. However, climate change may contribute to 14.3–31% of streamflow change in north China [74,75], while human activities, such as land-use change and SWCPs, would also be an important factor influencing runoff and soil erosion [68,76].

3.4. The Effect of Different SWCPs on Soil and Water Conservaion

Seven SWCPs and a bare slope without any SWCPs were studied using SCS-CN and MULSE with RCP4.5 and RCP8.5 during 2020–2050 to quantify the effect of different SWCPs on soil and water erosion under future climate change. The predicted surface runoff and soil erosion of experimental plots are presented in Figure 2. The results showed the reduction rates of 15–40% for surface runoff and 35–67% for soil erosion under RCP4.5 and RCP8.5 scenarios during 2020–2050, respectively. This means that SWCPs performed well for the slopes in terms of soil and water conservation. Previous studies also reported that more than 90% of the erosion resulted from the bare slopes without vegetation cover, which was significantly higher than that on the slopes with vegetation cover [7,77]. Compared with historical climate, under traditional tillage, the WEPP model (Watershed Erosion Prediction Project) predicted a change of 49–112% and 29–79% for runoff, 31–167% and 2–81% for soil erosion during 2010–2039 and 2070–2099, respectively [78,79]. However, conservation farming produced low runoff and soil erosion. Therefore, it could reduce the adverse effect of future climate changes on agro-ecosystems. Compared to traditional farming under the historical climate, conservation farming could reduce runoff by −34–71% and 18–38% and decrease soil erosion by 26–77% and 56–68% during 2010–2039 and 2070–2099, respectively [69,78]. Moreover, the reduction rates of soil erosion by different SWCPs were mostly higher than those of surface runoff because surface runoff was simulated by the SCS-CN method which was closely related to rainfall. The response of soil erosion was nonlinear with rainfall and was estimated by MULSE, which was affected by rainfall amount and intensity, the soil moisture content in the early stage, soil physical properties, SWCPs, land use, and vegetation coverage [77].
The reduction rates of surface runoff and soil erosion varied among experimental plots with different SWCPs (Figure 2). These SWCPs could reduce surface runoff and soil erosion because the topography of the watershed could be changed to greatly delay surface runoff, improve rainfall infiltration, and reduce the runoff discharge. A large proportion of rainfall can be retained by vegetation and evaporate into the atmosphere, causing a significant decrease in the effective rainfall for runoff generation [68]. Under the same rainfall conditions, plot A (shrub/Seabuckthorn + horizontal terracing), B (grass + horizontal terracing), and H (forest/Chinese Pine + horizontal terracing) showed high reduction rates of soil erosion (about 54, 68, and 67%, respectively), which were due to the high vegetation cover (50.5–71.0%). It can be seen that plot A showed a lower reduction effect (about 15%) on surface runoff than other plots, which was probably because a part of retained rainfall by shrubs flew down the trunks or dropped from leaves to the surface quickly and produced secondary surface runoff and soil erosion. Previous studies have shown that ecological factors (e.g., tree canopy, shrub and grass coverage, litter depth, and woody debris) may affect the surface runoff and soil erosion. Specifically, the size and erosive power of raindrops were influenced by tree canopy. Ground cover, such as short shrub, grass, woody debris, and litter cover, offered a rough surface that reduced soil particle detachment, movement, and transport down the slope [72,80]. Therefore, stemflow and throughfall from tall forest may cause more secondary soil erosion than short shrubs because of strong erosive power of raindrops and less ground cover in this study.
Water costs for sediment control (Rrs, mm (t hm−2 yr−1)-1) of seven SWCPs under two scenarios (RCP4.5 and RCP8.5) during the period 2020–2050 are presented in Figure 3. The results showed that Rrs was in the range of 0.35–3.13 and 0.34–1.90 under RCP4.5 and RCP8.5 scenarios during 2020–2050, respectively. When reducing the soil erosion by similar quantities, plot A with the lowest Rrs, nearly 0.41 mm (t hm−2 yr−1)−1, had smaller impacts on runoff, while plot C (forest/Chinese Pine + fish-scale pitting) with the largest Rrs of 1.29 mm (t hm−2 yr−1)−1 reduced larger runoff compared to other SWCPs. Considering the Rrs of individual SWCPs in the hilly loess region, Yan [81] showed that the Rrs slightly varied among different SWCPs, with the order: shrub > deciduous broadleaved forest > evergreen coniferous forest > fish-scale pitting > horizontal terracing > grass. It suggested that the combination of forest and fish-scale pitting could result in a higher Rrs than a single SWCP, while horizontal terracing would reduce Rrs when combined with shrub because Rrs simultaneously decreased with the increase of rainfall, runoff, and soil erosion [82].

4. Conclusions and Recommendations

SWCPs have obvious effects on reducing soil and water loss. The results of this study showed that surface runoff and soil erosion were closely related to rainfall depth, rainfall intensity, and rainfall erosivity. Surface runoff showed a 28% increase compared to historical climate, while soil erosion presented large variations with years due to uneven distribution of rainfall and rainfall erosivity in a year. Compared with bare slopes, the reduction rates were 15–40% for surface runoff and 35–67% for soil erosion under RCP4.5 and RCP8.5 scenarios during 2020–2050, respectively. When reducing the soil erosion by similar quantities, the combination of shrub and horizontal terracing had smaller impacts on runoff compared to other SWCPs. Growing vegetation captured eroded sediment from the upper part of the slopes. Moreover, vegetation growth stabilized the slope, which decreased gravitational erosion [71]. In the long run, vegetation can help decrease surface runoff and soil erosion to a safe level through improving soil chemical-physical properties and anti-erodibility [30]. However, compared to engineering practices, the effect of biological practices on soil erosion conservation might be postponed because vegetation like trees, shrubs, and grass require a long period to adapt to the changes in soil features, topography, and local climate [17]. Therefore, the combination of shrub and horizontal terracing is the desired recommendation for controlling surface runoff and soil erosion because this combination has smaller impacts on runoff when reducing the soil erosion by similar quantities. The results of this study can be used to assess soil and water erosion, support small watershed management, and preserve the ecosystem of erosive areas.

Author Contributions

These authors contributed equally to this work. Their individual contributions were as follows: conceptualization, M.C. and C.A.; methodology, M.C. and C.A.; software, M.C. and C.L.; validation, M.C.; formal analysis, M.C.; investigation, M.C.; resources, M.C. and C.L.; writing—original draft preparation, M.C.; writing—review and editing, M.C., C.A. and C.G.; supervision, C.A. and C.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Natural Sciences and Engineering Research Council of Canada.

Acknowledgments

The authors are grateful to the editor and the anonymous reviewers for their insightful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lal, R. Soil erosion and the global carbon budget. Environ. Int. 2003, 29, 437–450. [Google Scholar] [CrossRef]
  2. Pennock, D.; McKenzie, N.; Montanarella, L. Status of the world’s soil resources; FAO Rome: Rome, Italy, 2015. [Google Scholar]
  3. Montanarella, L. Agricultural policy: govern our soils. Nat. News 2015, 528, 32. [Google Scholar] [CrossRef] [PubMed]
  4. Pimentel, D.; Burgess, M. Soil erosion threatens food production. Agriculture 2013, 3, 443–463. [Google Scholar] [CrossRef] [Green Version]
  5. MWR of PRC. The proclamation of soil and water loss in China; MWR: Beijing, China, 2002. [Google Scholar]
  6. Wang, Y.; Gao, X.; Li, T.; Yue, Y.; Fang, H.; Ni, J. Geocode-based aquatic habitats in hierarchical system of the Yellow River Basin. J. Environ. Inform. 2018, 32, 69–81. [Google Scholar] [CrossRef]
  7. Fu, B.; Liu, Y.; Lü, Y.; He, C.; Zeng, Y.; Wu, B. Assessing the soil erosion control service of ecosystems change in the Loess Plateau of China. Ecol. Complex. 2011, 8, 284–293. [Google Scholar] [CrossRef]
  8. Wang, Z.; Zhang, T.; Tan, C.; Wang, X.; Taylor, R.; Qi, Z.; Yang, J. Modeling the impacts of manure on phosphorus loss in surface runoff and subsurface drainage. J. Environ. Qual. 2019, 48, 39–46. [Google Scholar] [CrossRef] [Green Version]
  9. Sawut, R.; Kasim, N.; Maihemuti, B.; Hu, L.; Abliz, A.; Abdujappar, A.; Kurban, M. Pollution characteristics and health risk assessment of heavy metals in the vegetable bases of northwest China. Sci. Total Environ. 2018, 642, 864–878. [Google Scholar] [CrossRef]
  10. Zhang, P.; Huang, G.; An, C.; Fu, H.; Gao, P.; Yao, Y.; Chen, X. An integrated gravity-driven ecological bed for wastewater treatment in subtropical regions: Process design, performance analysis, and greenhouse gas emissions assessment. J. Clean. Prod. 2019, 212, 1143–1153. [Google Scholar] [CrossRef]
  11. Colman, C.B.; Oliveira, P.T.S.; Almagro, A.; Soares-Filho, B.S.; Rodrigues, D.B. Effects of climate and land-cover changes on soil erosion in Brazilian pantanal. Sustainability 2019, 11, 7053. [Google Scholar] [CrossRef] [Green Version]
  12. Xu, Q.; Chen, Q.; Zhao, S.; Liu, K.; Ma, J. Saving water and associated energy from distribution networks by considering landscape factors in pressure management and use of district metered areas. J. Environ. Inform. 2018, 31, 65–73. [Google Scholar] [CrossRef]
  13. Routschek, A.; Schmidt, J.; Kreienkamp, F. Impact of climate change on soil erosion—a high-resolution projection on catchment scale until 2100 in Saxony/Germany. Catena 2014, 121, 99–109. [Google Scholar] [CrossRef]
  14. Lu, C.; Huang, G.; Wang, X. Projected changes in temperature, precipitation, and their extremes over China through the RegCM. Clim. Dyn. 2019, 53, 5859–5880. [Google Scholar] [CrossRef]
  15. Delgado, J.A.; Nearing, M.A.; Rice, C.W. Conservation practices for climate change adaptation. In Advances in Agronomy; Elsevier: Amsterdam, The Netherlands, 2013; Volume 121, pp. 47–115. [Google Scholar]
  16. Fu, B.; Li, S.; Yu, X.; Yang, P.; Yu, G.; Feng, R.; Zhuang, X. Chinese ecosystem research network: progress and perspectives. Ecol. Complex. 2010, 7, 225–233. [Google Scholar] [CrossRef]
  17. Millington, J.D.; Wainwright, J.; Mulligan, M. Representing human decision-making in environmental modelling. Environ. Model. 2013, 291. [Google Scholar]
  18. Li, P.; Mu, X.; Holden, J.; Wu, Y.; Irvine, B.; Wang, F.; Gao, P.; Zhao, G.; Sun, W. Comparison of soil erosion models used to study the Chinese Loess Plateau. Earth-Sci. Rev. 2017, 170, 17–30. [Google Scholar] [CrossRef] [Green Version]
  19. Sadeghi, S.; Gholami, L.; Khaledi Darvishan, A.; Saeidi, P. A review of the application of the MUSLE model worldwide. Hydrol. Sci. J. 2014, 59, 365–375. [Google Scholar] [CrossRef] [Green Version]
  20. Abdul-Rahim, A.; Sun, C.; Noraida, A. The impact of soil and water conservation on agricultural economic growth and rural poverty reduction in China. Sustainability 2018, 10, 4444. [Google Scholar] [CrossRef] [Green Version]
  21. Zhao, G.; Mu, X.; Wen, Z.; Wang, F.; Gao, P. Soil erosion, conservation, and eco-environment changes in the Loess Plateau of China. Land Degrad. Dev. 2013, 24, 499–510. [Google Scholar] [CrossRef]
  22. Cai, M.; Xin, Z.; Yu, X. Spatio-temporal variations in PM leaf deposition: A meta-analysis. Environ. Pollut. 2017, 231, 207–218. [Google Scholar] [CrossRef]
  23. Morar, F.; Iantovics, L.; Gligor, A. Analysis of phytoremediation potential of crop plants in industrial heavy metal contaminated soil in the upper Mures River Basin. J. Environ. Inform. 2018, 31, 1–14. [Google Scholar]
  24. Cai, M.; Xin, Z.; Yu, X. Particulate matter transported from urban greening plants during precipitation events in Beijing, China. Environ. Pollut. 2019, 252, 1648–1658. [Google Scholar] [CrossRef]
  25. Santonastaso, G.; Erto, A.; Bortone, I.; Chianese, S.; Di Nardo, A.; Musmarra, D. Experimental and simulation study of the restoration of a thallium (I)-contaminated aquifer by Permeable Adsorptive Barriers (PABs). Sci. Total Environ. 2018, 630, 62–71. [Google Scholar] [CrossRef] [Green Version]
  26. Raptis, S.; Gasparatos, D.; Economou-Eliopoulos, M.; Petridis, A. Chromium uptake by lettuce as affected by the application of organic matter and Cr (VI)-irrigation water: Implications to the land use and water management. Chemosphere 2018, 210, 597–606. [Google Scholar] [CrossRef] [PubMed]
  27. Liu, B.; Liu, Y.; Zhang, K.; Xie, Y. Classification for soil conservation practices in China. J. Soil Water Conserv. 2013, 27, 80–84. [Google Scholar]
  28. Zhao, J.; Yang, Z.; Govers, G. Soil and water conservation measures reduce soil and water losses in China but not down to background levels: evidence from erosion plot data. Geoderma 2019, 337, 729–741. [Google Scholar] [CrossRef]
  29. Peng, H.; Jia, Y.; Niu, C.; Gong, J.; Hao, C.; Gou, S. Eco-hydrological simulation of soil and water conservation in the Jinghe River Basin in the Loess Plateau, China. J. Hydro-Environ. Res. 2015, 9, 452–464. [Google Scholar] [CrossRef]
  30. Zhou, Y.; Yang, B.; Han, J.; Huang, Y. Robust linear programming and its application to water and environmental decision-making under uncertainty. Sustainability 2019, 11, 33. [Google Scholar] [CrossRef] [Green Version]
  31. Huang, J.; Huang, G.; An, C.; He, Y.; Yao, Y.; Zhang, P.; Shen, J. Performance of ceramic disk filter coated with nano ZnO for removing Escherichia coli from water in small rural and remote communities of developing regions. Environ. Pollut. 2018, 238, 52–62. [Google Scholar] [CrossRef]
  32. Zhou, Y.; Huang, G.; Baetz, B.W. Multilevel factorial fractional programming for sustainable water resources management. J. Water Resour. Plan. Manag. 2016, 142, 04016063. [Google Scholar] [CrossRef]
  33. Wang, T.; Wu, J.; Kou, X.; Oliver, C.; Mou, P.; Ge, J. Ecologically asynchronous agricultural practice erodes sustainability of the Loess Plateau of China. Ecol. Appl. 2010, 20, 1126–1135. [Google Scholar] [CrossRef]
  34. Li, S.; Di, X.; Wu, D.; Zhang, J. Effects of sewage sludge and nitrogen fertilizer on herbage growth and soil fertility improvement in restoration of the abandoned opencast mining areas in Shanxi, China. Environ. Earth Sci. 2013, 70, 3323–3333. [Google Scholar] [CrossRef]
  35. Shi, H.; Shao, M. Soil and water loss from the Loess Plateau in China. J. Arid. Environ. 2000, 45, 9–20. [Google Scholar] [CrossRef]
  36. Yin, Z.; Gou, J.; Li, Y. Design strategies and soil erosion reduction effect of soil and water conservation measures in hilly area of central Sichuan Basin. Syst. Sci. Compr. Stud. Agric. 2009, 25, 369–374. (In Chinese) [Google Scholar]
  37. Wang, F.; Li, R.; Mu, X.; Yang, X. Characteristics and regional differences of runoff cost of sediment control of conservation measures and their simulation in Weihe River Basin. Sci. Soil Water Conserv. 2004, 2, 12–17. (In Chinese) [Google Scholar]
  38. Change, I.C. The Physical Science Basis. In Contribution of working group i to the fifth assessment report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2013; Volume 2013, pp. 33–118. [Google Scholar]
  39. Gao, G.Y.; Fu, B.J.; Lü, Y.H.; Liu, Y.; Wang, S.; Zhou, J. Coupling the modified SCS-CN and RUSLE models to simulate hydrological effects of restoring vegetation in the Loess Plateau of China. Hydrol. Earth Syst. Sci. 2012, 16, 2347–2364. [Google Scholar] [CrossRef] [Green Version]
  40. Kendon, E.J.; Rowell, D.P.; Jones, R.G.; Buonomo, E. Robustness of future changes in local precipitation extremes. J. Clim. 2008, 21, 4280–4297. [Google Scholar] [CrossRef]
  41. Gudmundsson, L.; Bremnes, J.; Haugen, J.; Engen-Skaugen, T. Downscaling RCM precipitation to the station scale using statistical transformations–a comparison of methods. Hydrol. Earth. Syst. Sci. 2012, 16, 3383–3390. [Google Scholar] [CrossRef] [Green Version]
  42. Baltas, E.; Dervos, N.; Mimikou, M. Determination of the SCS initial abstraction ratio in an experimental watershed in Greece. Hydrol. Earth Syst. Sci. Discuss. 2007, 11, 1825–1829. [Google Scholar] [CrossRef] [Green Version]
  43. Mishra, S.K.; Singh, V.P. Validity and extension of the SCS-CN method for computing infiltration and rainfall-excess rates. Hydrol. Process. 2004, 18, 3323–3345. [Google Scholar] [CrossRef]
  44. Wang, H. Estimating watersheds runoff by using improved Runoff Curve Number Model (SCS-CN) on the Loess Plateau of China. Master’s Thesis, Beijing Forestry University, Bejing, China, 2016. (In Chinese). [Google Scholar]
  45. Wischmeier, W.H.; Smith, D.D. Predicting rainfall erosion losses: A guide to conservation planning; Science and Education Administration, US Department of Agriculture: Washington, DC, USA, 1978.
  46. Renard, K.G.; Foster, G.R.; Weesies, G.A.; Porter, J.P. RUSLE: Revised universal soil erosion equation. J. Soil Water Conserv. 1991, 46, 30–33. [Google Scholar]
  47. Liu, B.; Zhang, K.; Xie, Y. An empirical soil erosion equation. In Proceedings of the 12th International Soil Conservation Organization Conference, Beijing, China, 26–31 May 2002; Tsinghua University Press: Beijing, China, 2002; III, p. 15. [Google Scholar]
  48. Angulo-Martínez, M.; Beguería, S. Estimating rainfall erosivity from daily precipitation records: A comparison among methods using data from the Ebro Basin (NE Spain). J. Hydrol. 2009, 379, 111–121. [Google Scholar] [CrossRef] [Green Version]
  49. Wu, L.; Liu, X.; Ma, X. Spatiotemporal distribution of rainfall erosivity in the Yanhe River watershed of hilly and gully region, Chinese Loess Plateau. Environ. Earth Sci. 2016, 75, 315. [Google Scholar] [CrossRef]
  50. Qin, W.; Guo, Q.; Cao, W.; Yin, Z.; Yan, Q.; Shan, Z.; Zheng, F. A new RUSLE slope length factor and its application to soil erosion assessment in a Loess Plateau watershed. Soil Till. Res. 2018, 182, 10–24. [Google Scholar] [CrossRef]
  51. Zhang, K.; Yu, Y.; Dong, J.; Yang, Q.; Xu, X. Adapting & testing use of USLE K factor for agricultural soils in China. Agricul. Ecosyst. Environ. 2019, 269, 148–155. [Google Scholar]
  52. Smith, D.D.; Wischmeier, W.H. Factors affecting sheet and rill erosion. Eos Trans. Am. Geophys. Union 1957, 38, 889–896. [Google Scholar] [CrossRef]
  53. Malagó, A.; Vigiak, O.; Bouraoui, F.; Pagliero, L.; Franchini, M. The hillslope length impact on SWAT streamflow prediction in large basins. J. Environ. Inform. 2018, 32, 82–97. [Google Scholar] [CrossRef] [Green Version]
  54. Sun, W.; Shao, Q.; Liu, J.; Zhai, J. Assessing the effects of land use and topography on soil erosion on the Loess Plateau in China. Catena 2014, 121, 151–163. [Google Scholar] [CrossRef]
  55. Kinnell, P. Event soil loss, runoff and the Universal Soil Loss Equation family of models: A review. J. Hydrol. 2010, 385, 384–397. [Google Scholar] [CrossRef]
  56. Renard, K.G.; Foster, G.R. Soil conservation: principles of erosion by water. Dryland. Agric. 1983, 23, 155–176. [Google Scholar]
  57. Nash, J.; Sutcliffe, J. River forcasting using conceptual models, 1. A discussion of principles. J. Hydrol. 1970, 10, 280–290. [Google Scholar] [CrossRef]
  58. Remesan, R.; Bray, M.; Mathew, J. Application of PCA and clustering methods in input selection of hybrid runoff models. J. Environ. Inform. 2018, 31, 137–152. [Google Scholar] [CrossRef] [Green Version]
  59. Oğuz, I. Rainfall erosivity in north-central anatolia in turkey. Appl. Ecol. Environ. Res. 2019, 17, 2719–2731. [Google Scholar]
  60. Soulis, K.; Valiantzas, J.; Dercas, N.; Londra, P. Analysis of the runoff generation mechanism for the investigation of the SCS-CN method applicability to a partial area experimental watershed. Hydrol. Earth Syst. Sci. Discuss. 2009, 6, 373–400. [Google Scholar] [CrossRef]
  61. Soulis, K.; Valiantzas, J. SCS-CN parameter determination using rainfall-runoff data in heterogeneous watersheds–the two-CN system approach. Hydrol. Earth Syst. Sci. 2012, 16, 1001–1015. [Google Scholar] [CrossRef] [Green Version]
  62. Xiao, B.; Wang, Q.; Fan, J.; Han, F.; Dai, Q. Application of the SCS-CN model to runoff estimation in a small watershed with high spatial heterogeneity. Pedosphere 2011, 21, 738–749. [Google Scholar] [CrossRef]
  63. Li, C.; Qin, J.; Li, J. Application of computational curve number to precipitation-runoff simulation in a typical watershed in Chinese Loess Plateau. J. Arid. Land. Resour. Environ. 2008, 22, 67–70. (In Chinese) [Google Scholar]
  64. Zhang, K.; Li, S.; Peng, W.; Yu, B. Erodibility of agricultural soils on the Loess Plateau of China. Soil Till. Res. 2004, 76, 157–165. [Google Scholar] [CrossRef] [Green Version]
  65. Wang, B.; Zheng, F.; Guan, Y. Improved USLE- K factor prediction: A case study on water erosion areas in China. Int. Soil Water Conserv. Res. 2016, 4, 168–176. [Google Scholar] [CrossRef] [Green Version]
  66. Auerswald, K.; Fiener, P.; Martin, W.; Elhaus, D. Use and misuse of the K factor equation in soil erosion modeling: An alternative equation for determining USLE nomograph soil erodibility values. Catena 2014, 118, 220–225. [Google Scholar] [CrossRef]
  67. Cao, L.; Zhang, K.; Dai, H.; Liang, Y. Modeling interrill erosion on unpaved roads in the Loess Plateau of China. Land Degrad. Dev. 2015, 26, 825–832. [Google Scholar] [CrossRef]
  68. Guo, Q.; Han, Y.; Yang, Y.; Fu, G.; Li, J. Quantifying the impacts of climate change, coal mining and soil and water conservation on streamflow in a coal mining concentrated watershed on the Loess Plateau, China. Water 2019, 11, 1054. [Google Scholar] [CrossRef] [Green Version]
  69. Li, Z.; Liu, W.; Zhang, X.; Zheng, F. Assessing the site-specific impacts of climate change on hydrology, soil erosion and crop yields in the Loess Plateau of China. Clim. Chang. 2011, 105, 223–242. [Google Scholar] [CrossRef]
  70. Azari, M.; Saghafian, B.; Moradi, H.R.; Faramarzi, M. Effectiveness of soil and water conservation practices under climate change in the Gorganroud Basin, Iran. CLEAN - Soil Air Water 2017, 45, 1700288. [Google Scholar] [CrossRef]
  71. Barthes, B.; Roose, E. Aggregate stability as an indicator of soil susceptibility to runoff and erosion; validation at several levels. Catena 2002, 47, 133–149. [Google Scholar] [CrossRef] [Green Version]
  72. Hartanto, H.; Prabhu, R.; Widayat, A.S.; Asdak, C. Factors affecting runoff and soil erosion: plot-level soil loss monitoring for assessing sustainability of forest management. For. Ecol. Manag. 2003, 180, 361–374. [Google Scholar] [CrossRef]
  73. Peng, T.; Wang, S. Effects of land use, land cover and rainfall regimes on the surface runoff and soil loss on karst slopes in southwest China. Catena 2012, 90, 53–62. [Google Scholar] [CrossRef]
  74. Wang, X.; He, K.; Dong, Z. Effects of climate change and human activities on runoff in the Beichuan River Basin in the northeastern Tibetan Plateau, China. Catena 2019, 176, 81–93. [Google Scholar] [CrossRef]
  75. Li, Y.; Liu, C.; Yu, W.; Tian, D.; Bai, P. Response of streamflow to environmental changes: A Budyko-type analysis based on 144 river basins over China. Sci. Total Environ. 2019, 664, 824–833. [Google Scholar] [CrossRef]
  76. Zare, M.; Mohammady, M.; Pradhan, B. Modeling the effect of land use and climate change scenarios on future soil loss rate in Kasilian watershed of northern Iran. Environ. Earth Sci. 2017, 76, 305. [Google Scholar] [CrossRef]
  77. Zheng, F. Effect of vegetation changes on soil erosion on the Loess Plateau. Pedosphere 2006, 16, 420–427. [Google Scholar] [CrossRef]
  78. Zhang, X.; Liu, W.; Li, Z.; Zheng, F. Simulating site-specific impacts of climate change on soil erosion and surface hydrology in southern Loess Plateau of China. Catena 2009, 79, 237–242. [Google Scholar] [CrossRef]
  79. Zhang, X.; Liu, W. Simulating potential response of hydrology, soil erosion, and crop productivity to climate change in Changwu tableland region on the Loess Plateau of China. Agric. For. Meteorol. 2005, 131, 127–142. [Google Scholar] [CrossRef]
  80. Montenegro, A.D.A.; Abrantes, J.; De Lima, J.; Singh, V.; Santos, T. Impact of mulching on soil and water dynamics under intermittent simulated rainfall. Catena 2013, 109, 139–149. [Google Scholar] [CrossRef]
  81. Yan, L. Runoff cost for sediment control of soil and water conservation practices at the plot’s scale on the Loess Plateau. Master’s Thesis, Northwest A&F University, Yangling, China, 2008. (In Chinese). [Google Scholar]
  82. Mo, L. Runoff cost of sediment control of soil and water conservation practices with hydrologic methods in Beiluo River Basin. Master’s Thesis, Northwest A&F University, Yangling, China, 2008; pp. 15–18. (In Chinese). [Google Scholar]
Figure 1. Relative differences in rainfall (P), rainfall erosivity factor (R), surface runoff (Q), and soil erosion (A) between the historical climate (1986–2005) and the two RCP scenarios (RCP4.5 and RCP8.5) during 2020–2050.
Figure 1. Relative differences in rainfall (P), rainfall erosivity factor (R), surface runoff (Q), and soil erosion (A) between the historical climate (1986–2005) and the two RCP scenarios (RCP4.5 and RCP8.5) during 2020–2050.
Sustainability 12 00934 g001
Figure 2. The reduction rates of different SWCPs for two scenarios (RCP4.5 and RCP8.5) during the period of 2020–2050: (a) surface runoff (fQ, %) under RCP4.5 scenario; (b) surface runoff (fQ, %) under RCP8.5 scenario; (c) soil erosion (fA, %) under RCP4.5 scenario; (d) soil erosion (fA, %) under RCP8.5 scenario.
Figure 2. The reduction rates of different SWCPs for two scenarios (RCP4.5 and RCP8.5) during the period of 2020–2050: (a) surface runoff (fQ, %) under RCP4.5 scenario; (b) surface runoff (fQ, %) under RCP8.5 scenario; (c) soil erosion (fA, %) under RCP4.5 scenario; (d) soil erosion (fA, %) under RCP8.5 scenario.
Sustainability 12 00934 g002
Figure 3. Water cost for sediment control (Rrs, mm (t hm−2 yr−1)−1) of different SWCPs during the period of 2020–2050: (a) Rrs under RCP4.5 scenario; (b) Rrs under RCP8.5 scenario.
Figure 3. Water cost for sediment control (Rrs, mm (t hm−2 yr−1)−1) of different SWCPs during the period of 2020–2050: (a) Rrs under RCP4.5 scenario; (b) Rrs under RCP8.5 scenario.
Sustainability 12 00934 g003
Table 1. Basic characteristics of eight experimental plots.
Table 1. Basic characteristics of eight experimental plots.
SWCPs 1ABCDEFGH
Land useshrubgrassforestbareforestforestforestforest
Vegetation cover (%)7161.738.51028.524.616.650.5
P-factor0.0270.0500.2031.0000.2030.1810.1880.181
Erosive events4234254438403739
Duration (min)9771078119510171054105910501053
Intensity (mm h−1)2.32.42.62.32.32.42.32.3
Rainfall (P, mm)26.329.834.827.829.129.628.329.0
Surface runoff (mm)1.21.01.41.71.41.21.11.1
Soil erosion (t hm−2 yr−1)7.668.1020.3329.2218.4216.1619.1811.41
1 SWCPs: A (Sea-buckthorn + horizontal terracing); B (grass + horizontal terracing); C (forest/Chinese Pine + fish-scale pitting); D (bare slope); E (forest/Arborvitae + fish-scale pitting); F (forest/Arborvitae + horizontal terracing); G (forest/Apricot + horizontal terracing); H (forest/Chinese Pine + horizontal terracing).
Table 2. Curve number (CN) value of experimental plots with different soil and water conservation practices (SWCPs) from January to December.
Table 2. Curve number (CN) value of experimental plots with different soil and water conservation practices (SWCPs) from January to December.
SWCPsABCDEFGH
126.821.915.738.316.17.522.929.3
226.821.915.738.316.17.522.929.3
326.821.915.738.316.17.522.929.3
426.821.915.738.316.17.522.929.3
526.821.915.738.316.17.522.929.3
639.029.110.628.923.524.222.213.3
722.219.822.330.826.725.217.422.3
826.820.323.228.622.818.323.623.9
935.417.934.431.620.420.523.525.7
1026.829.420.651.238.738.845.141.8
1126.829.420.651.238.738.845.141.8
1226.829.420.651.238.738.845.141.8
Average28.923.719.238.824.220.228.029.8
Table 3. Soil erodibility factor (K, t hm2 h hm−2 MJ−1 mm−1) of experimental plots with different SWCPs.
Table 3. Soil erodibility factor (K, t hm2 h hm−2 MJ−1 mm−1) of experimental plots with different SWCPs.
SWCPsABCDEFGH
SAN (%)20.419.224.315.313.517.122.117.9
SIL (%)53.752.75258.765.759.156.158.5
CLA (%)27.320.29.116.817.816.910.918.3
SOM (%)0.90.80.651.21.11.160.991.24
Dg (mm)0.91.10.70.80.70.80.90.8
p23422232
s33333333
Knomo0.0400.0410.0470.0430.0530.0450.0470.044
Kepic0.0520.0540.0540.0540.0560.0530.0530.052
KDg0.0380.0310.0480.0420.0480.0420.0380.042
Kopt0.0020.0010.0090.0050.0120.0060.0060.005
Table 4. Root mean square error (RMSE) and Nash-Sutcliffe efficiency (NSE) of observed and predicted surface runoff (Q, mm) and soil erosion (A, t hm−2 yr−1) from experimental plots with different SWCPs.
Table 4. Root mean square error (RMSE) and Nash-Sutcliffe efficiency (NSE) of observed and predicted surface runoff (Q, mm) and soil erosion (A, t hm−2 yr−1) from experimental plots with different SWCPs.
SWCPsMULSECalibrationValidation
QAQA
abObservedPredictedRMSENSEObservedPredictedRMSENSEObservedPredictedRMSENSEObservedPredictedRMSENSE
A3.1091.6211.381.350.650.724.283.525.040.541.011.210.660.659.9710.758.310.65
B3.4091.6310.820.840.280.812.041.551.590.931.141.440.820.5811.8612.248.390.46
C2.8491.2311.331.210.660.5811.58.666.60.751.381.940.870.5828.4827.428.050.88
D1.6811.151.61.760.470.9215.915.8619.510.531.722.151.110.7438.4446.4314.530.89
E1.6091.2331.191.260.650.826.676.836.710.731.541.740.880.7625.2824.315.440.57
F3.0231.28811.180.490.837.447.547.220.651.341.540.970.6721.426.9814.560.7
G2.5091.3190.961.150.40.7411.77.5410.730.541.161.421.110.4322.3428.519.510.43
H2.5591.2581.011.150.480.834.684.274.60.721.181.641.120.4115.6115.1512.80.48
Average 1.181.260.530.818.067.189.710.591.311.621.000.6321.5524.1914.30.72

Share and Cite

MDPI and ACS Style

Cai, M.; An, C.; Guy, C.; Lu, C. Assessment of Soil and Water Conservation Practices in the Loess Hilly Region Using a Coupled Rainfall-Runoff-Erosion Model. Sustainability 2020, 12, 934. https://doi.org/10.3390/su12030934

AMA Style

Cai M, An C, Guy C, Lu C. Assessment of Soil and Water Conservation Practices in the Loess Hilly Region Using a Coupled Rainfall-Runoff-Erosion Model. Sustainability. 2020; 12(3):934. https://doi.org/10.3390/su12030934

Chicago/Turabian Style

Cai, Mengfan, Chunjiang An, Christophe Guy, and Chen Lu. 2020. "Assessment of Soil and Water Conservation Practices in the Loess Hilly Region Using a Coupled Rainfall-Runoff-Erosion Model" Sustainability 12, no. 3: 934. https://doi.org/10.3390/su12030934

APA Style

Cai, M., An, C., Guy, C., & Lu, C. (2020). Assessment of Soil and Water Conservation Practices in the Loess Hilly Region Using a Coupled Rainfall-Runoff-Erosion Model. Sustainability, 12(3), 934. https://doi.org/10.3390/su12030934

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