Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Effects of Crop Rotation and Topography on Soil Erosion and Nutrient Loss under Natural Rainfall Conditions on the Chinese Loess Plateau

Effects of Crop Rotation and Topography on Soil Erosion and Nutrient Loss under Natural Rainfall... land Article Effects of Crop Rotation and Topography on Soil Erosion and Nutrient Loss under Natural Rainfall Conditions on the Chinese Loess Plateau 1 , 2 3 2 , Chenhui Li , Wenhai Shi and Mingbin Huang * College of Resources and Environment, Northwest A & F University, Xianyang 712100, China State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Northwest A & F University, Xianyang 712100, China School of Water and Environment, Chang’an University, Xi’an 710054, China * Correspondence: hmbd@nwsuaf.edu.cn Abstract: Erosive rainfall results in the loss of both soil and nutrients, which indirectly triggers soil deterioration and a reduction in land productivity. However, how rainfall affects runoff, soil erosion, and nutrient loss under different crop rotation patterns and topographic factors remains unclear. This experiment observed nine runoff-erosion plots on the Chinese Loess Plateau (CLP) from 2019 to 2020 to determine the effects of crop type, rotation pattern, and slope gradient and length on runoff, soil erosion, and nutrient loss. Runoff, soil erosion, and nutrient loss were highest for the fallow plots; values for these variables for spring corn and winter wheat plots were not significantly different. Crop rotation generated greater runoff, soil erosion, and nutrient loss compared to non-rotation. Soil erosion and associated nutrient loss increased, but not significantly, with slope for gradients of 0.5 , 1 , and 3 , while runoff and associated nutrient loss did not increase. In addition, soil erosion and associated nutrient loss were significantly greater for slope lengths of 20 m vs. 50 m. A structural equation model showed rainfall characteristics significantly impacted runoff and soil erosion and subsequently affected nutrient loss. This study increases the understanding of runoff, soil erosion, and nutrient loss from cropland with gentle slopes on the CLP. Keywords: crop rotation; runoff and sediment; nitrogen and phosphorus losses; Chinese Loess Plateau Citation: Li, C.; Shi, W.; Huang, M. Effects of Crop Rotation and Topography on Soil Erosion and Nutrient Loss under Natural Rainfall 1. Introduction Conditions on the Chinese Loess Soil erosion driven by runoff and sediment generally takes nutrients away from soil, Plateau. Land 2023, 12, 265. https:// reduces soil fertility, and then causes soil degradation [1,2]. Therefore, soil erosion has doi.org/10.3390/land12020265 attracted widespread attention as a worldwide problem [3–5], especially in severe soil Academic Editor: Vincent Chaplot erosion areas, such as the Chinese Loess Plateau (CLP). The CLP is the most erosive area in the world, generating approximately 90% of the sediment and nutrients that enter the Received: 17 December 2022 Yellow River [6], as well as a large amount of soil nitrogen and phosphorus that cause Revised: 6 January 2023 Accepted: 14 January 2023 water pollution and eutrophication [7]. Soil nutrients are critical to the productivity of Published: 17 January 2023 terrestrial ecosystems and have received a large amount of attention [8,9]. Some research shows inappropriate land use and tillage practices can increase soil nutrient loss [10], which usually occurs via two mechanisms—sediment-bound nutrients lost as particles and soluble nutrients lost with runoff [11]. Runoff and soil erosion from farmland are considered the Copyright: © 2023 by the authors. main sources of soil nutrients entering water bodies [12,13]. Therefore, investigating the soil Licensee MDPI, Basel, Switzerland. erosion and nutrient loss from farmland is important for better management of agricultural This article is an open access article production and the ecological environment. distributed under the terms and Soil nutrient loss processes are complex and controlled by various factors, such as conditions of the Creative Commons rainfall, soil tillage system, and topography [14–16]. The main driver of soil erosion in Attribution (CC BY) license (https:// semi-arid systems is rainfall, with characteristics such as duration, amount, and intensity creativecommons.org/licenses/by/ controlling hydrological processes [17,18]. Rainfall characteristics are also the main factors 4.0/). Land 2023, 12, 265. https://doi.org/10.3390/land12020265 https://www.mdpi.com/journal/land Land 2023, 12, 265 2 of 16 that affect soil nutrient loss and have been extensively studied [19,20]. The sediment yield rate decreases with increasing rainfall duration [21], and even the form of nutrient loss varies with rainfall duration [22]. Sediment concentration and soil loss increase as rainfall amount increases [23]. In general, heavy rainfalls of short duration and high intensity cause the most erosion [24] and also drive the redistribution of nutrients [25]. Rainfall reaching an intensity of 15 mm/h can cause severe soil erosion [26]. However, previous studies have mainly focused on simulated rainfall [27–29] and long-term natural rainfall experiments are relatively rare [30,31]. Therefore, soil nutrient migration in farmland ecosystems needs to be characterized under natural precipitation conditions. Crop type and rotation are other major factors that affect soil erosion and nutrient loss from farmland [32,33]. The direct impact of cultivating different crops on soil loss and water erosion has been widely studied. Carroll et al. (1997) show that, in central Queensland, the average annual runoff from wheat is lower and the soil loss higher than that from sorghum and sunflower [34]. Prasuhn (2012) demonstrates that, among 203 crop fields in the Swiss Midlands, the greatest total soil erosion occurs in winter wheat fields, followed by potato and fallow fields [35]. Based on 198 rainfall-runoff events, Fiene and Auerswald (2007) show the obvious difference in soil loss based on the previous season crop [36]. However, inadequate information is available about the effect of winter wheat and spring corn (common crops on the CLP) on soil and nutrient loss in surface runoff. Moreover, research on crop rotation has mainly focused on crop yield and soil moisture [37,38], while relatively little has been carried out to study soil erosion and nutrient loss in the CLP region [33,39]. Therefore, conditions of soil erosion and nutrient loss under common crop types and rotation need to be clarified to inform the development of more effective cropping systems to replace those in current use. Topographic factors such as slope gradient and length also affect soil erosion and nutrient loss. In some soil erosion models (such as USLE, CSLE, and EUROSEM), slope gradient and length are important parameters [40–42]. Many studies show the amount of soil erosion and associated nutrient loss are positively related to the gradient of culti- vated land [43,44]. The occurrence of soil erosion will especially increase when crops are cultivated on farmland with steep slopes [45]. However, some reports indicate total nitro- gen and total phosphorus in surface runoff decrease with increasing slope gradient [46]. Slope length also plays a key role in soil erosion and nutrient loss. Lal (1997) shows total runoff nutrient losses decrease with increasing slope length under conventional tillage practices, whereas the opposite result is obtained under no-till practices [47]. Xing et al. (2016) indicate that sediment-associated total nitrogen loss increases with slope length under simulated rainfall conditions [48]. Overall, studies to date show slope gradient and length can have either positive or negative impacts on runoff, erosion, and nutrient loss and have mainly focused on soil erosion and nutrient loss for steep slopes. The influence of slope gradient and length on soil erosion from gently sloping agricultural land has not received sufficient attention, and studies are needed to facilitate the development of new technologies for erosion control on gently sloping agricultural fields. Overall, the influence of rainfall, crop type and rotation, and topography on soil erosion and nutrient loss under gently sloping agricultural land in the CLP is not well understood. Therefore, the objectives of this study were to: (1) identify the runoff, soil erosion, and nutrient loss for different crops and crop rotation systems; (2) explain the effects of slope gradient and length on runoff, soil erosion, and nutrient loss for gently sloping agricultural land; (3) investigate how rainfall factors affect runoff, soil erosion, and nutrient loss. The results of this study can clarify the nutrient loss pathways of runoff and sediment from cropland and may also provide a theoretical basis for solving problems caused by nutrient loss. Land 2023, 12, x FOR PEER REVIEW 3 of 16 Land 2023, 12, 265 3 of 16 off and sediment from cropland and may also provide a theoretical basis for solving problems caused by nutrient loss. 2. Materials and Methods 2. Materials and Methods 2.1. Study Site 2.1. Study Site 0 00 The study was conducted in the Wangdonggou watershed (latitude: 107 40 30 – The study was conducted in the Wangdonggou watershed (latitude: 107°40′30″– 0 00  0 00  0 00 107 42 30 E; longitude: 35 12 16 –35 16 00 N; altitude: 946–1226 m) at the Changwu 107°42′30″ E; longitude: 35°12′16″–35°16′00″ N; altitude: 946–1226 m) at the Changwu Agro-ecological Experiment Station, Chinese Academy of Sciences (CAS) and Ministry Agro-ecological Experiment Station, Chinese Academy of Sciences (CAS) and Ministry of Water Resources (MWR), in Shaanxi Province, China. The watershed is a typical of Water Resources (MWR), in Shaanxi Province, China. The watershed is a typical loess 2 2 loess tableland-gully region with an area of 8.3 km and a gully density of 2.78 km/km 2 2 tableland-gully region with an area of 8.3 km and a gully density of 2.78 km/km (Fig- (Figure 1). The gradient of tableland varies from 0–3 . The study area has a warm-temperate ure 1). The gradient of tableland varies from 0–3°. The study area has a warm-temperate semi-humid continental monsoon climate with an average annual temperature of 9.2 C semi-humid continental monsoon climate with an average annual temperature of 9.2 °C (1957–2014) [49]. The average annual precipitation is 579 mm (1960–2016), with more than (1957–2014) [49]. The average annual precipitation is 579 mm (1960–2016), with more 70% falling from April to September [50]. The soil in this watershed is dominated by Heilu than 70% falling from April to September [50]. The soil in this watershed is dominated soil [51], and the parent material is medium loamy Malan loess. The thickness of the loess by Heilu soil [51], and the parent material is medium loamy Malan loess. The thickness layer is more than 100 m, and the soil in the whole profile is uniform and loose with good of the loess layer is more than 100 m, and the soil in the whole profile is uniform and air permeability. loose with good air permeability. Figure 1. Location of the study site and experimental design. (a) the Loess Plateau of China, (b) Wangdonggou watershed on the Loess Plateau of China, (c) runoff and sediment collecting equipment, and (d) nine runoff plots. Numbers of 1–9 show the locations of all nine runoff plots. Land 2023, 12, 265 4 of 16 2.2. Experimental Design Nine loess experimental plots were arranged at Changwu station with three crop/land use types: winter wheat (Triticum aestivum L.), spring corn (Zea mays L.) and fallow (Table 1). Winter wheat was sown in September and harvested in early June; spring corn was sown in April and harvested in September. The boundaries of each plot were constructed by reinforced concrete baffles with 35 cm buried and 15 cm above ground to separate runoff from inside and outside the plots. A self-recording rainfall gauge was placed around the experimental plots to record rainfall. Each plot had a set of equipment for collecting runoff and sediment, including a water collection tank, water conveyance tank, water-retaining dike, and runoff barrel. Table 1. Basic information about the study plots. Plot 1 2 3 4 5 6 7 8 9 Slope length (m) 20 20 20 50 50 50 20 20 20 Slope gradient ( ) 0.5 0.5 0.5 0.5 0.5 0.5 0.5 1 3 Crop rotation Yes Yes Yes Yes Yes Yes No No No 2019 SC WW FA SC WW FA WW WW WW Crop/land type 2020 FA SC WW FA SC WW WW WW WW 2019 78.0 80.0 47.3 74.7 78.0 48.0 79.3 76.7 78.0 Coverage (%) 2020 67.3 75.0 45.0 69.7 73.3 45.7 76.7 75.0 76.7 Note: FA: Fallow; SC: Spring corn; WW: Winter wheat. Coverage was measured in October of 2019 and 2020, respectively. The coverage of winter wheat fields took into account the presence of weeds and stubble at that time. Plots 1 to 3 had a slope length of 20 m and gradient of 0.5 with a wheat-corn-fallow rotation treatment; plots 4 to 6 had a slope length of 50 m and gradient of 0.5 with a wheat-corn-fallow rotation treatment; plots 7 to 9 had a slope length of 20 m and gradients of 0.5 , 1 , and 3 , respectively, with a winter wheat monoculture (Figure 1). Table 1 shows crop types planted during the study period (2019 and 2020). Data from plots 2 and 7 in 2019 and plots 3 and 7 in 2020 for winter wheat were used to compare the effects of rotation and no-rotation on runoff, soil erosion, and nutrient loss. Data from plots 1 and 4 in 2019 and plots 2 and 5 in 2020 for spring corn, data from plots 2 and 5 in 2019 and plots 3 and 6 in 2020 for winter wheat, and data from plots 3 and 6 in 2019 and plots 1 and 4 in 2020 for fallow were used to analyze the effects of different crop/land uses on runoff, soil erosion, and nutrient loss. Data from plots 7, 8, and 9 were used to compare the effect of slope gradient on runoff, soil erosion, and nutrient loss with precipitation as a random variable. Data from plots 1 to 6 were used to contrast the effects of slope length and crop type on runoff, soil erosion, and nutrient loss. 2.3. Data Collection Data were collected from May to October of 2019 and 2020 because most of the rainfall and runoff on the CLP occurs during these months. The total rainfall, rainfall duration, average rainfall intensity (ARI), and the maximum rainfall intensity (MRI) for each rainfall event were obtained from the self-recording rainfall gauge and are given in Table 2. The runoff volume of each plot was obtained by measuring the height of the inside of the runoff barrel according to the known barrel bottom area after the occurrence of erosive rainfall. The suspension (runoff + sediment) in the runoff barrel was thoroughly stirred and samples were collected using 5000-mL plastic bottles. Some samples were filtered through medium-speed phosphorus-free filter paper, and approximately 300-mL samples were immediately refrigerated (4 C) and brought back to the laboratory for determination of total nitrogen and total phosphorus in 24 h. The total nitrogen of water samples was measured by alkaline potassium persulfate ablation UV spectrophotometry and the total phosphorus of water samples was measured by ammonium molybdate spectrophotometry. The remaining samples were separated by filtering through filter paper after settling for Land 2023, 12, 265 5 of 16 24 h. Total sediment was calculated by air drying sediment samples and weighing them to obtain the sediment sample mass. Subsequently, the sediment samples were brought back to the laboratory for determination of total nitrogen and total phosphorus. Soil total nitrogen was measured by the semi-micro Kjeltec method, and soil total phosphorus was measured by the HClO –H SO elimination method. The vegetation cover, plant spacing, 4 2 4 and height of the plots were measured, and soil samples were collected in September and October 2019 for determination of their physical and chemical properties. The basic physical and chemical indicators are shown in Table 3. Table 2. Erosive rainfall characteristics during the study period. Date Rainfall Rainfall Duration ARI MRI (d/m/y) (mm) (h) (mm/h) (mm/h) 05/06/2019 56.5 18.5 3.1 6.8 28/06/2019 49.3 28.2 1.8 3.0 23/07/2019 28.6 4.0 6.3 7.2 29/07/2019 34.2 4.7 7.3 9.8 04/08/2019 20.5 13.3 1.5 3.7 26/08/2019 71.6 12.5 5.7 25.2 15/09/2019 70.0 8.0 3.6 8.8 19/06/2020 58.5 53.0 1.1 3.6 26/06/2020 66.0 17.3 3.8 14.8 11/07/2020 41.4 24.0 1.7 5.0 26/07/2020 47.0 32.5 1.4 4.6 07/08/2020 31.0 6.0 5.2 7.8 13/08/2020 20.8 9.5 2.2 4.6 19/08/2020 88.3 40.0 2.2 6.0 Note: MRI: maximum rainfall intensity; ARI: average rainfall intensity. Table 3. Soil properties of experimental plots. SOC TN TP AP AK pH BD Clay Silt Sand Soil Texture Plot (g/kg) (g/kg) (g/kg) (mg/kg) (mg/kg) – (g/cm ) (%) (%) (%) – 1 9.33 0.98 0.99 14.51 139.82 7.93 1.21 27.92 37.1 34.98 Loamy clay 2 8.67 1.02 0.87 10.01 125.66 7.95 1.28 27.31 37.28 35.41 Loamy clay 3 8.07 0.97 0.81 7.87 117.56 7.96 1.27 27.78 37.42 34.8 Loamy clay 4 7.78 0.91 0.99 10.85 112.84 8.08 1.25 28.08 39.04 32.88 Loamy clay 5 8.28 0.96 0.97 13.19 136.62 7.96 1.24 27.75 37.35 34.9 Loamy clay 6 8.15 0.98 0.94 7.91 134.93 7.79 1.48 26.63 36.77 36.6 Loamy clay 7 9.03 1.11 0.86 7.57 117.73 7.77 1.29 26.3 37.76 35.94 Loamy clay 8 8.85 1.07 0.84 11.46 107.95 7.72 1.38 26.9 39.42 33.68 Loamy clay 9 8.1 0.97 0.84 9.38 113.35 7.74 1.3 25.29 36.29 38.42 Loamy clay Note: SOC: soil organic carbon; TN: total nitrogen; TP: total phosphorus; AP: available phosphorus; AK: available potassium; BD: bulk density. 2.4. Data Analysis All statistical analyses were performed using R version 4.0.2. The mixed linear model and multiple comparisons based on Tukey’s test were used to identify the differences in the effects of tillage practices (rotation and no-rotation), crop type (spring corn, winter wheat, and fallow), and topography (slope gradient and length) on runoff, soil erosion, and nutrient loss; significance differences were judged based on the R package lmerTest [52]. The main factors influencing nutrient loss were decided using a random forest model [53], and the significance of each predictor was assessed using the R package rfPermute [54]. A structural equation model was constructed in combination with professional knowledge using R package piecewiseSEM [55], which was used to quantify the direct and indirect effects of explanatory variables on response variables. Land 2023, 12, x FOR PEER REVIEW 6 of 16 winter wheat, and fallow), and topography (slope gradient and length) on runoff, soil erosion, and nutrient loss; significance differences were judged based on the R package lmerTest [52]. The main factors influencing nutrient loss were decided using a random forest model [53], and the significance of each predictor was assessed using the R pack- age rfPermute [54]. A structural equation model was constructed in combination with professional knowledge using R package piecewiseSEM [55], which was used to quantify the direct and indirect effects of explanatory variables on response variables. Land 2023, 12, 265 6 of 16 3. Results 3.1. Runoff, Soil Erosion, and Nutrient Loss for Different Crops/Land Uses 3. Results Runoff, soil erosion, and nutrient loss for the spring corn, winter wheat, and fallow 3.1. Runoff, Soil Erosion, and Nutrient Loss for Different Crops/Land Uses plots are compared in Figure 2. The fallow plots had more runoff, soil erosion, and nu- Runoff, soil erosion, and nutrient loss for the spring corn, winter wheat, and fallow trient loss than the spring corn or winter wheat plots. Specifically, the fallow plots were plots are compared in Figure 2. The fallow plots had more runoff, soil erosion, and significantly different from the spring corn and winter wheat plots in terms of runoff nutrient loss than the spring corn or winter wheat plots. Specifically, the fallow plots were dept significantly h, soil erosion, different tot fraom l phosphorus the spring cor in n runo and winter ff (RTP wheat ), toplots tal nit inrterms ogen in of r so unof il flos depth, s (STN), soil erosion, total phosphorus in runoff (RTP), total nitrogen in soil loss (STN), and total and total phosphorus in soil loss (STP). For total nitrogen in runoff (RTN), the fallow phosphorus in soil loss (STP). For total nitrogen in runoff (RTN), the fallow plots were plots were only significantly different from the spring corn plots and not from winter only significantly different from the spring corn plots and not from winter wheat plots wheat plots (Figure 2c). Spring corn plots had a lower runoff depth and associated RTP (Figure 2c). Spring corn plots had a lower runoff depth and associated RTP and RTN but and RTN but higher soil loss and associated STP and STN compared to winter wheat higher soil loss and associated STP and STN compared to winter wheat plots; however, the plots; however, the differences were not significant. Overall, runoff, soil erosion, and differences were not significant. Overall, runoff, soil erosion, and nutrient loss from the nutrient loss from the fallow plots were significantly higher than for spring corn and fallow plots were significantly higher than for spring corn and winter wheat plots, but the winter wheat plots, but the differences between the latter two crops were not significant. differences between the latter two crops were not significant. Figure 2. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 2. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil (STP) for spring corn, winter wheat, and fallow plots. Values followed by different letters above loss (STP) for spring corn, winter wheat, and fallow plots. Values followed by different letters a column are significantly different at p  0.05. Different letters indicate significant differences above a column are significantly different at p ≤ 0.05. Different letters indicate significant differ- between treatments. ences between treatments. 3.2. Runoff, Soil Erosion, and Nutrient Loss for Different Rotation Systems The runoff, soil erosion, and nutrient loss for winter wheat cropped under rotation and no-rotation systems are compared in Figure 3. Runoff depth, soil erosion, RTN, and STP were significantly higher under the rotation system than under the no-rotation system in 2019; only RTP was significantly higher in 2020. Combining data from the two years, the rotation system resulted in significantly greater runoff depth, RTN, and Land 2023, 12, x FOR PEER REVIEW 7 of 16 3.2. Runoff, Soil Erosion, and Nutrient Loss for Different Rotation Systems The runoff, soil erosion, and nutrient loss for winter wheat cropped under rotation and no-rotation systems are compared in Figure 3. Runoff depth, soil erosion, RTN, and STP were significantly higher under the rotation system than under the no-rotation sys- Land 2023, 12, 265 7 of 16 tem in 2019; only RTP was significantly higher in 2020. Combining data from the two years, the rotation system resulted in significantly greater runoff depth, RTN, and STN than the no-rotation system by 72.74, 99.41, and 182.06%, respectively. For soil erosion, STN than the no-rotation system by 72.74, 99.41, and 182.06%, respectively. For soil RTP, and STP, the rotation system yielded higher values than the no-rotation system, erosion, RTP, and STP, the rotation system yielded higher values than the no-rotation but the differences were not significant. Overall, the use of a rotation system for winter system, but the differences were not significant. Overall, the use of a rotation system wheat resulted in more runoff, soil erosion, and nutrient loss compared to a no-rotation for winter wheat resulted in more runoff, soil erosion, and nutrient loss compared to a system. no-rotation system. Figure 3. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 3. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss (STP) in rotation and no-rotation systems for 2019 and 2020. ALL: both years (2019 and 2020); *, loss (STP) in rotation and no-rotation systems for 2019 and 2020. ALL: both years (2019 and 2020); p < 0.05; **, p < 0.01; *** p < 0.001. *, p < 0.05; **, p < 0.01; *** p < 0.001. 3.3. Impact of Slope Gradient and Length on Runoff, Soil Erosion, and Nutrient Loss 3.3. Impact of Slope Gradient and Length on Runoff, Soil Erosion, and Nutrient Loss A comparison of runoff, soil erosion, and nutrient loss for gradients of 0.5 , 1 , and 3 A comparison of runoff, soil erosion, and nutrient loss for gradients of 0.5°, 1°, and is illustrated in Figure 4. Soil loss and associated STP and STN losses increased with the 3° is illustrated in Figure 4. Soil loss and associated STP and STN losses increased with gradient but differences between the various slopes were not significant. Runoff depth and the gradient but differences between the various slopes were not significant. Runoff associated RTP and RTN losses did not increase with slope gradient. In general, slope had depth and associated RTP and RTN losses did not increase with slope gradient. In gen- no significant effect on soil erosion and nutrient loss for the gentle slope range considered era (0.5–3 l, slope h ). ad no significant effect on soil erosion and nutrient loss for the gentle slope range considered (0.5–3°). Runoff, soil erosion, and nutrient loss for different slope lengths for the three land use types considered are shown in Figure 5. Plots with a slope length of 50 m had lower runoff depth, soil erosion, RTN, STN, RTP, and STP than plots with a slope length of 20 m for all crop/land use types, which implies runoff, soil erosion, and nutrient loss de- Land 2023, 12, x FOR PEER REVIEW 8 of 16 crease with slope length. For spring corn and winter wheat, runoff depth, soil erosion, RTN, STN, and STP were significantly higher for the 20- vs. 50-m slope length; RTP was also higher for the 20-m slope, but the difference was not significant. For fallow plots, runoff depth and RTN were significantly lower for the 50-m slope; soil loss, RTP, STN, and STP were also lower, but the differences were not significant. When all crop/land use types were considered, increasing slope length resulted in significant decreases in Land 2023, 12, 265 8 of 16 runoff depth, RTN, and STN and insignificant decreases in soil loss, RTP, and STP. These results indicate that slope length significantly affected runoff and nitrogen loss. Figure 4. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 4. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss loss (STP) for slope gradients of 0.5°  , 1°, and 3  °. Values followed by different letters above a column (STP) for slope gradients of 0.5 , 1 , and 3 . Values followed by different letters above a column are are significantly different at p ≤ 0.05. significantly different at p  0.05. Runoff, soil erosion, and nutrient loss for different slope lengths for the three land use types considered are shown in Figure 5. Plots with a slope length of 50 m had lower runoff depth, soil erosion, RTN, STN, RTP, and STP than plots with a slope length of 20 m for all crop/land use types, which implies runoff, soil erosion, and nutrient loss decrease with slope length. For spring corn and winter wheat, runoff depth, soil erosion, RTN, STN, and STP were significantly higher for the 20- vs. 50-m slope length; RTP was also higher for the 20-m slope, but the difference was not significant. For fallow plots, runoff depth and RTN were significantly lower for the 50-m slope; soil loss, RTP, STN, and STP were also lower, but the differences were not significant. When all crop/land use types were considered, increasing slope length resulted in significant decreases in runoff depth, RTN, and STN and insignificant decreases in soil loss, RTP, and STP. These results indicate that slope length significantly affected runoff and nitrogen loss. Land 2023, 12, x FOR PEER REVIEW 9 of 16 Land 2023, 12, 265 9 of 16 Figure 5. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 5. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss loss (STP) for slope lengths of 20 and 50 m. The significance value in the upper left corner repre- (STP) for slope lengths of 20 and 50 m. The significance value in the upper left corner represents the sents the results of the analysis of all crops. *, p < 0.05; **, p < 0.01; *** p < 0.001. results of the analysis of all crops. *, p < 0.05; **, p < 0.01; *** p < 0.001. 3.4. Impact of Natural Rainfall Characteristics on Runoff, Soil Erosion, and Nutrient Loss 3.4. Impact of Natural Rainfall Characteristics on Runoff, Soil Erosion, and Nutrient Loss The influence of total rainfall, rainfall duration, MRI, ARI, soil condition, runoff depth, The influence of total rainfall, rainfall duration, MRI, ARI, soil condition, runoff and soil loss on RTN, RTP, STN, and STP losses were studied using a random forest model depth, and soil loss on RTN, RTP, STN, and STP losses were studied using a random (Figure 6). The pathways and coefficients affecting runoff, soil erosion, and nutrient loss forest model (Figure 6). The pathways and coefficients affecting runoff, soil erosion, and are presented in a structural equation model for each crop/land use type (Figure 7). The nutrient loss are presented in a structural equation model for each crop/land use type results show rainfall was the main driving force related to the occurrence of runoff and soil erosion, which indirectly leads to nutrient loss. For runoff, the structural equation model (Figure 7). The results show rainfall was the main driving force related to the occurrence shows the total rainfall and rainfall duration significantly altered runoff depth in the spring of runoff and soil erosion, which indirectly leads to nutrient loss. For runoff, the struc- corn plots (path coefficients of 0.714 and 0.335, respectively; Figure 7a), while the total tural equation model shows the total rainfall and rainfall duration significantly altered rainfall only significantly influenced runoff depth in the winter wheat and fallow plots runoff depth in the spring corn plots (path coefficients of 0.714 and −0.335, respectively; (path coefficients of 0.593 and 0.668, respectively; Figure 7b,c). Importantly, MRI did not Figure 7a), while the total rainfall only significantly influenced runoff depth in the winter significantly affect runoff depth for any of the three crop/land use types. Soil loss was wheat and fallow plots (path coefficients of 0.593 and 0.668, respectively; Figure 7b,c). significantly affected by runoff depth in spring corn and fallow plots (path coefficients of Importantly, MRI did not significantly affect runoff depth for any of the three crop/land 0.403 and 0.845, respectively; Figure 7a,c), and significantly influenced by MRI in the winter use types. Soil loss was significantly affected by runoff depth in spring corn and fallow wheat plots (path coefficient of 0.404; Figure 7b). plots (path coefficients of 0.403 and 0.845, respectively; Figure 7a,c), and significantly in- fluenced by MRI in the winter wheat plots (path coefficient of 0.404; Figure 7b). Land 2023, 12, x FOR PEER REVIEW 10 of 16 Land 2023, 12, x FOR PEER REVIEW 10 of 16 Land 2023, 12, 265 10 of 16 Figure 6. Random forest mean predictor importance (% increase in mean square error) of the fac- Figure Figure 6. 6. Random Random forest mean predictor im forest mean predictor importance portance (% (% incr increase in m ease in mean ean squar squ eaerr re error) or) of the of the factors fac- tors influencing (a) total nitrogen in runoff (RTN), (b) total phosphorus in runoff (RTP), (c) total tors influencing (a) total nitrogen in runoff (RTN), (b) total phosphorus in runoff (RTP), (c) total influencing (a) total nitrogen in runoff (RTN), (b) total phosphorus in runoff (RTP), (c) total nitrogen nitrogen in soil loss (STN), and (d) total phosphorus in soil loss (STP) for spring corn, (e–h) winter nitrogen in soil loss (STN), and (d) total phosphorus in soil loss (STP) for spring corn, (e–h) winter in soil loss (STN), and (d) total phosphorus in soil loss (STP) for spring corn, (e–h) winter wheat, and wheat, and (i–l) fallow plots, respectively. TR: total rainfall; SC: soil condition (including soil nu- wheat, and (i–l) fallow plots, respectively. TR: total rainfall; SC: soil condition (including soil nu- (i–l) fallow plots, respectively. TR: total rainfall; SC: soil condition (including soil nutrients, bulk trients, bulk density, and saturated hydraulic conductivity); RDe: runoff depth; RDu: rainfall du- trients, bulk density, and saturated hydraulic conductivity); RDe: runoff depth; RDu: rainfall du- density, and saturated hydraulic conductivity); RDe: runoff depth; RDu: rainfall duration; MRI: ration; MRI: maximum rainfall intensity; ARI: average rainfall intensity; SL: soil erosion. *, p < 0.05; ration; MRI: maximum rainfall intensity; ARI: average rainfall intensity; SL: soil erosion. *, p < 0.05; maximum rainfall intensity; ARI: average rainfall intensity; SL: soil erosion. *, p < 0.05; **, p < 0.01. **, p < 0.01. **, p < 0.01. Figure Figure 7. 7. Structural Structural equati equation on mod models els shshow ow the the pathways pathways and coeff and coef icificients ents affecting affecting runoff runof , sofi,l e soil ro- Figure 7. Structural equation models show the pathways and coefficients affecting runoff, soil ero- sion, and nutrient loss. Insignificant paths (p ≥ 0.05) are semi-transparent dashed lines, while sig- erosion, and nutrient loss. Insignificant paths (p  0.05) are semi-transparent dashed lines, while sion, and nutrient loss. Insignificant paths (p ≥ 0.05) are semi-transparent dashed lines, while sig- nificant paths (p < 0.05) are shown with solid lines. Black and red arrows indicate positive and significant paths (p < 0.05) are shown with solid lines. Black and red arrows indicate positive and nificant paths (p < 0.05) are shown with solid lines. Black and red arrows indicate positive and negative pathways, respectively. Conditional R denotes the proportion of variance explained. The negative negative pathways, respe pathways, respectively ctively. Conditio . Conditional nal R R denote denotes s the the proportion of proportion of variance variance explained. The explained. The standardized path coefficient represents the extent of the effect of the independent variable on the standardized path coefficient represents the extent of the effect of the independent variable on the standardized path coefficient represents the extent of the effect of the independent variable on the dependent variable. Significant pathways are marked by: *, p < 0.05; **, p < 0.01; ***, p < 0.001. dependent variable. Significant pathways are marked by: *, p < 0.05; **, p < 0.01; ***, p < 0.001. dependent variable. Significant pathways are marked by: *, p < 0.05; **, p < 0.01; ***, p < 0.001. Land 2023, 12, 265 11 of 16 For nutrient loss, the random forest model results indicate total rainfall, MRI, ARI, runoff depth, and soil erosion all had important effects (Figure 6). The results of structural equation modeling further demonstrate rainfall indirectly causes nutrient loss by affecting runoff depth and soil erosion (Figure 7). Runoff depth explains 83, 71, and 91% of the total variation in RTN loss and 49, 44, and 91% of the total variation in RTP loss for spring corn, winter wheat, and fallow plots, respectively. Soil erosion explains 97, 94, and 100% of the total variation in STN loss and 96, 97, and 96% of the total variation in STP loss for spring corn, winter wheat, and fallow plots, respectively (Figure 7). However, some differences were evident among three plots. For the fallow plots (Figure 7c), runoff depth and soil loss significantly impacted RTN, RTP, STN, and STP losses. For the winter wheat plots, runoff depth only had a significantly positive effect on RTN and RTP losses (path coefficients of 0.905 and 0.554, respectively; Figure 7b), while soil erosion had a significantly positive effect on STN and STP losses (path coefficients of 0.938 and 0.999, respectively; Figure 7b). For spring corn plots, runoff depth only had a significantly positive effect on RTN loss (path coefficient of 0.877; Figure 7a), while soil loss had a significantly negative effect on RTP loss and a significantly positive effect on STN and STP losses (path coefficients of 0.421, 0.938, and 0.999, respectively; Figure 7a). 4. Discussion 4.1. Impact of Crop Type on Runoff, Soil Erosion, and Nutrient Loss Our results indicate natural rainfall causes greater runoff, soil erosion, and nutrient loss in fallow plots compared to spring corn and winter wheat plots (Figure 2). Many studies also report that the choice of crop can effectively reduce surface runoff and soil and nutrient losses [34,35]. In the study region, the growing season of spring corn is from April to September, which aligns with the rainy season [56,57]. The coverage of spring corn can reduce soil erosion by intercepting rainfall and decreasing raindrop kinetic energy [57]. In addition, the active root system can increase erosion resistance. However, the growing season of winter wheat is from September to May of the next year, which does not align with the rainy season. Therefore, runoff, soil erosion, and nutrient losses are higher but not significantly greater in spring corn vs. winter wheat fields [58]. In our study, the spring corn and winter wheat plots both had greater biomass and coverage than the fallow plot, which likely resulted in the reduced soil erosion and nutrient loss. 4.2. Impact of Crop Rotation on Runoff, Soil Erosion, and Nutrient Loss Crop rotation is an important management practice [59]. Our results show rotation can result in greater runoff, soil erosion, and nutrient loss than no-rotation management in winter wheat plots (Figure 3). Crop rotation changes the coverage time and duration on the soil surface as well as the root distribution in the topsoil. In addition, farmers use different methods to harvest winter wheat and spring corn. Plots in the no-rotation system are left as wheat stubble after harvest, which provides better surface cover and roughness and reduces runoff and soil erosion [60,61]. Moreover, wheat residuals can increase soil organic matter content, improve soil physical characteristics, and subsequently reduce runoff, soil erosion, and nutrient loss [62,63]. For RTP, the major fraction of soil phosphorus is tightly adsorbed to mineral particles and bound within organic matter, so only a very small fraction leaches as dissolved soluble phosphate [64,65]. This might be the reason for the insignificant difference in phosphorus loss noted in our results for the different rotation systems. In this study, the results show that crop rotation generates greater runoff, which is different from the conclusion of Jiao et al. (2011) [66]. Jiao et al. (2011) compares surface runoff for three double cropping systems of wheat/maize, wheat/cotton, and wheat/soybean with wheat/fallow and indicates that all double cropping systems increases the ground cover and significantly reduces runoff [66,67]. However, our result is drawn for the same crop type and coverage of winter wheat with rotation and no-rotation. Land 2023, 12, 265 12 of 16 4.3. Impact of Slope Gradient and Length on Runoff, Soil Erosion, and Nutrient Loss Sloping farmland is one of the main types of agricultural land on the CLP. The gradient affects runoff, soil erosion, and nutrient loss by altering runoff flow rates and infiltration rates [68,69]. Our results show soil erosion and associated STP and STN losses increased with slope gradient, but the differences among the three gentle slope gradients (0.5 , 1 , and 3 ) were not significant (Figure 4b–f). Runoff depth and associated RTP and RTN losses did not increase with slope gradient, perhaps because the gentle gradients considered in this study were not sufficiently large to create stronger effects [68,70]. Slope length also impacts runoff, soil erosion, and nutrient loss [48,71]. Our study results show all indicators for the 20-m slope length were greater than those for the 50-m slope length. Among all indicators, runoff depth, RTN and STN were significantly larger for the 20-m slope length than for the 50-m slope length when all crop/land use types were considered (Figure 5a,c,e); these differences were not significant for soil erosion, RTP, and STP (Figure 5b,d,f). These findings confirm previous results reported by Xing et al. (2016) [48], where both runoff and associated total nitrogen loss rate decreased with increasing slope length. This is attributed to the longer slope increasing surface runoff infiltration and deposition of eroded soil along its length [72,73], which subsequently reduces runoff-associated RTN and RTP losses and sediment-associated STN and STP losses [74,75]. Furthermore, some studies suggest soil loss decreases with slope length according to a power-law trend [76,77], which might explain why the soil loss was less for the 50-m vs. 20-m slope. However, the two years of data from this study show significant nitrogen loss from runoff and soil but insignificant phosphorus loss for the different slope lengths. This might be due to RTN values being much higher than RTP values (Figure 5c,d). The latter were usually insignificant, likely due to phosphorus being predominantly and firmly associated with soil particle surfaces [78,79]; the amount of STP supports this suggestion. 4.4. Impact of Rainfall on Runoff, Soil Erosion, and Nutrient Loss The results of our study show rainfall characteristics are the dominant driver of runoff, soil erosion, and associated nutrient loss (Figures 6 and 7). This result is consistent with those of many other studies [14,20,80]. Runoff depth was mainly positively associated with total rainfall and negatively associated with rainfall duration (Figure 7), which were in agreement with earlier studies [81,82]. These results suggest total rainfall and rainfall duration determine the runoff depth, which indirectly causes runoff-associated nutrient loss. Our results are also consistent with assertions that runoff and MRI are the main factors influencing soil loss [83]. This is because MRI can reflect rainfall impact on soil particle detachment, whereas surface runoff depth can be used to represent the transport of sediments [11]. Our study also found rainfall indirectly causes nutrient loss by affecting runoff depth and soil erosion, with nutrient loss due to both N and P dissolved in runoff as well as partially bonded to soil particles [11,84]. The forms of nutrient loss also differed for the three crops/land uses studied (Figure 7). This result aligns with those of previous studies that show the intensity of rainfall can impact the form of nutrient loss [85]. At higher intensity the nutrient loss is attributed more to soil nutrients being transported with sediments, while at lower intensity the nutrient loss is attributed more soluble elements carried with runoff. Overall, the responses of runoff, soil erosion, and nutrient loss to rainfall intensity in plots of different crops differed and, thus, resulted in different forms of nutrient loss. 5. Conclusions This paper reports on an experiment featuring nine runoff-erosion plots on the CLP studied in 2019 to 2020 to determine the effect of crop type, rotation pattern, slope gradient, and slope length on runoff, soil erosion, and nutrient loss. Runoff, soil erosion, and nutrient loss were higher in fallow plots than in spring corn and winter wheat plots, with no significant differences noted between the latter two. Plots subject to crop rotation Land 2023, 12, 265 13 of 16 generated greater runoff, soil erosion, and nutrient loss compared to those managed with a no-rotation system. Soil loss and associated nutrient loss tended to increase with slope gradient, while runoff and associated nutrient loss did not show any trends. Runoff, soil erosion, and associated nutrient loss were significantly greater in plots with a slope length of 20- vs. 50-m under gentle gradient conditions. Total rainfall and rainfall duration were the primary factors determining runoff, soil erosion, and associated nutrient loss, and maximum rainfall intensity partially explained soil loss. Our results suggest that the effect of rainfall characteristics on soil nutrient loss varies with rotation system and crop type. These results are useful for developing targeted farmland nutrient conservation strategies. Author Contributions: Data collection, data analysis, original draft preparation and writing, C.L.; manuscript of validation and revision, W.S.; manuscript review, editing, and project supervision, M.H. All authors have read and agreed to the published version of the manuscript. Funding: This study was supported by the National Natural Science Foundation of China (42107351); the China Postdoctoral Science Foundation (2019M663917XB); the Natural Science Fund of Shaanxi Province (2021JQ-227). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: The vector boundary of the Wangdonggou watershed was downloaded from the Loess Plateau Science Data Centre, National Earth System Science Data Sharing Infrastruc- ture, National Science and Technology Infrastructure of China (http://loess.geodata.cn (accessed on 18 March 2021). The rest of the data are available upon request from the corresponding author. Acknowledgments: The authors thank the members of the Changwu Agri-Ecological Experimental Station, Chinese Academy of Sciences, for their assistance. Conflicts of Interest: The authors declare no conflict of interest. References 1. Lal, R. Soil erosion impact on agronomic productivity and environment quality. Crit. Rev. Plant Sci. 1998, 17, 319–464. [CrossRef] 2. Gomes, L.; Simoes, S.J.C.; Dalla Nora, E.L.; de Sousa-Neto, E.R.; Forti, M.C.; Ometto, J.P.H.B. Agricultural expansion in the Brazilian Cerrado: Increased soil and nutrient losses and decreased agricultural productivity. Land 2019, 8, 12. [CrossRef] 3. Liu, C.; Sui, J.; Wang, Z. Sediment load reduction in Chinese rivers. Int. J. Sediment Res. 2008, 23, 44–55. [CrossRef] 4. Ochoa, P.A.; Fries, A.; Mejia, D.; Burneo, J.I.; Ruiz-Sinoga, J.D.; Cerda, A. Effects of climate, land cover and topography on soil erosion risk in a semiarid basin of the Andes. Catena 2016, 140, 31–42. [CrossRef] 5. Martinez-Mena, M.; Carrillo-Lopez, E.; Boix-Fayos, C.; Almagro, M.; Franco, N.G.; Diaz-Pereira, E.; Montoya, I.; de Vente, J. Long-term effectiveness of sustainable land management practices to control runoff, soil erosion, and nutrient loss and the role of rainfall intensity in Mediterranean rainfed agroecosystems. Catena 2020, 187, 104352. [CrossRef] 6. Shao, F.; Tao, W.; Wang, Q.; Wu, J.; Su, L.; Sun, Y.; Zeng, S. A modified model for predicting nutrient loss in runoff using a time-varying mixing layer. J. Hydrol. 2021, 603, 127091. [CrossRef] 7. Chen, M.; Chen, J.; Sun, F. Agricultural phosphorus flow and its environmental impacts in China. Sci. Total Environ. 2008, 405, 140–152. [CrossRef] 8. Tan, Z.; Lal, R.; Wiebe, K.D. Global soil nutrient depletion and yield reduction. J. Sustain. Agric. 2005, 26, 123–146. [CrossRef] 9. Sharpley, A.; Kleinman, P. Effect of rainfall simulator and plot scale on overland flow and phosphorus transport. J. Environ. Qual. 2003, 32, 2172–2179. [CrossRef] 10. Fu, B.; Chen, L.; Ma, K.; Zhou, H.; Wang, J. The relationships between land use and soil conditions in the hilly area of the loess plateau in northern Shaanxi, China. Catena 2000, 39, 69–78. [CrossRef] 11. Shi, W.; Huang, M.; Wu, L. Prediction of storm-based nutrient loss incorporating the estimated runoff and soil loss at a slope scale on the Loess Plateau. Land Degrad. Dev. 2018, 29, 2899–2910. [CrossRef] 12. Bechmann, M.E.; Boe, F. Soil tillage and crop growth effects on surface and subsurface runoff, loss of soil, phosphorus and nitrogen in a cold climate. Land 2021, 10, 77. [CrossRef] 13. Guo, S.; Zhai, L.; Liu, J.; Liu, H.; Chen, A.; Wang, H.; Wu, S.; Lei, Q. Cross-ridge tillage decreases nitrogen and phosphorus losses from sloping farmlands in southern hilly regions of China. Soil Tillage Res. 2019, 191, 48–56. [CrossRef] 14. Zheng, F. Effects of accelerated soil erosion on soil nutrient loss after deforestation on the Loess Plateau. Pedosphere 2005, 15, 707–715. 15. Stevens, C.J.; Quinton, J.N.; Bailey, A.P.; Deasy, C.; Silgram, M.; Jackson, D.R. The effects of minimal tillage, contour cultivation and in-field vegetative barriers on soil erosion and phosphorus loss. Soil Tillage Res. 2009, 106, 145–151. [CrossRef] Land 2023, 12, 265 14 of 16 16. Zhang, X.; Hu, M.; Guo, X.; Yang, H.; Zhang, Z.; Zhang, K. Effects of topographic factors on runoff and soil loss in Southwest China. Catena 2018, 160, 394–402. [CrossRef] 17. Weltzin, J.F.; Tissue, D.T. Resource pulses in arid environments—Patterns of rain, patterns of life. New Phytol. 2003, 157, 171–173. [CrossRef] 18. Yu, Y.; Wei, W.; Chen, L.; Feng, T.; Daryanto, S. Quantifying the effects of precipitation, vegetation, and land preparation techniques on runoff and soil erosion in a Loess watershed of China. Sci. Total Environ. 2019, 652, 755–764. [CrossRef] 19. Jung, J.W.; Lim, S.S.; Kwak, J.H.; Park, H.J.; Yoon, K.S.; Kim, H.Y.; Baek, W.J.; Choi, W.J. Further Understanding of the impacts of rainfall and agricultural management practices on nutrient loss from rice paddies in a monsoon area. Water Air Soil Pollut. 2015, 226, 283. [CrossRef] 20. Lemma, B.; Kebede, F.; Mesfin, S.; Fitiwy, I.; Abraha, Z.; Norgrove, L. Quantifying annual soil and nutrient lost by rill erosion in continuously used semiarid farmlands, North Ethiopia. Environ. Earth Sci. 2017, 76, 190. [CrossRef] 21. Ma, L.; Li, J.; Liu, J. Effects of antecedent soil water content on infiltration and erosion processes on loessial slopes under simulated rainfall. Hydrol. Res. 2020, 51, 882–893. [CrossRef] 22. He, Z.; Weng, H.; Wu, T.; Ho, H.; Lin, Y.; Ran, Q. Impact of moving rainfall events on hillslope pollutant transport. Environ. Earth Sci. 2015, 74, 5989–5999. [CrossRef] 23. Nishigaki, T.; Sugihara, S.; Kilasara, M.; Funakawa, S. Surface runoff generation and soil loss under different soil and rainfall properties in the Uluguru Mountains, Tanzania. Land Degrad. Dev. 2017, 28, 283–293. [CrossRef] 24. Zheng, M.; Cai, Q.; Cheng, Q. Modelling the runoff-sediment yield relationship using a proportional function in hilly areas of the Loess Plateau, North China. Geomorphology 2008, 93, 288–301. 25. Fiener, P.; Dlugoss, V.; Van Oost, K. Erosion-induced carbon redistribution, burial and mineralisation—Is the episodic nature of erosion processes important? Catena 2015, 133, 282–292. 26. Foster, G.R.; Young, R.A.; Römkens, M.J.M.; Onstad, C.A. Processes of Soil Erosion by Water. In Soil Erosion and Crop Productivity; American Society of Agronomy, Inc.: Madison, WI, USA; Crop Science Society of America, Inc.: Madison, WI, USA; Soil Science Society of America, Inc.: Madison, WI, USA, 1985; Volume 9, pp. 137–162. 27. Shi, Z.; Fang, N.; Wu, F.; Wang, L.; Yue, B.; Wu, G. Soil erosion processes and sediment sorting associated with transport mechanisms on steep slopes. J. Hydrol. 2012, 454, 123–130. [CrossRef] 28. Cerda, A.; Gonzalez-Pelayo, O.; Gimenez-Morera, A.; Jordan, A.; Pereira, P.; Novara, A.; Brevik, E.C.; Prosdocimi, M.; Mahmood- abadi, M.; Keesstra, S.; et al. Use of barley straw residues to avoid high erosion and runoff rates on persimmon plantations in Eastern Spain under low frequency-high magnitude simulated rainfall events. Soil Res. 2016, 54, 154–165. [CrossRef] 29. Keesstra, S.D.; Rodrigo-Comino, J.; Novara, A.; Gimenez-Morera, A.; Pulido, M.; Di Prima, S.; Cerda, A. Straw mulch as a sustainable solution to decrease runoff and erosion in glyphosate-treated clementine plantations in Eastern Spain. An assessment using rainfall simulation experiments. Catena 2019, 174, 95–103. [CrossRef] 30. Hussein, M.H.; Kariem, T.H.; Othman, A.K. Predicting soil erodibility in northern Iraq using natural runoff plot data. Soil Tillage Res. 2007, 94, 220–228. [CrossRef] 31. Huo, J.; Liu, C.; Yu, X.; Chen, L.; Zheng, W.; Yang, Y.; Yin, C. Direct and indirect effects of rainfall and vegetation coverage on runoff, soil loss, and nutrient loss in a semi-humid climate. Hydrol. Process. 2021, 35, e13985. [CrossRef] 32. Martinez, J.R.F.; Zuazo, V.H.D.; Raya, A.M. Environmental impact from mountainous olive orchards under different soil- management systems (SE Spain). Sci. Total Environ. 2006, 358, 46–60. [CrossRef] [PubMed] 33. Yan, L.; Xue, L.; Petropoulos, E.; Qian, C.; Hou, P.; Xu, D.; Yang, L. Nutrient loss by runoff from rice-wheat rotation during the wheat season is dictated by rainfall duration. Environ. Pollut. 2021, 205, 117382. [CrossRef] [PubMed] 34. Carroll, C.; Halpin, M.; Burger, P.; Bell, K.; Sallaway, M.M.; Yule, D.F. The effect of crop type, crop rotation, and tillage practice on runoff and soil loss on a Vertisol in central Queensland. Aust. J. Soil Res. 1997, 35, 925–939. [CrossRef] 35. Prasuhn, V. On-farm effects of tillage and crops on soil erosion measured over 10 years in Switzerland. Soil Tillage Res. 2012, 120, 137–146. [CrossRef] 36. Fiene, P.; Auerswald, K. Rotation effects of potato, maize, and winter wheat on soil erosion by water. Soil Sci. Soc. Am. J. 2007, 71, 1919–1925. [CrossRef] 37. Zhang, Y.; Wang, R.; Wang, H.; Wang, S.; Wang, X.; Li, J. Soil water use and crop yield increase under different long-term fertilization practices incorporated with two-year tillage rotations. Agric. Water Manag. 2019, 221, 362–370. [CrossRef] 38. Cui, Z.; Yan, B.; Gao, Y.; Wu, B.; Wang, Y.; Xie, Y.; Xu, P.; Wang, H.; Wen, M.; Wang, Y.; et al. Crop yield and water use efficiency in response to long-term diversified crop rotations. Front. Plant. Sci 2022, 13, 1024898. [CrossRef] 39. Shen, Y.; Li, L.; Chen, W.; Robertson, M.; Unkovich, M.; Bellotti, W.; Probert, M. Soil water, soil nitrogen and productivity of lucerne-wheat sequences on deep silt loams in a summer dominant rainfall environment. Field Crop Res. 2009, 111, 97–108. [CrossRef] 40. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Index Losses: A Guide to Conservation Planning. In USDA Agricultural Handbook No. 537; US Department of Agriculture: Washington, DC, USA, 1981. 41. Liu, B.; Zhang, K.; Xie, Y. An empirical soil loss equation. In Process of Soil Erosion and Its Environment Effect, Proceedings of the 12th International Soil Conservation Organization Conference, Beijing, China, 26–31 May 2002; Tsinghua University Press: Beijing, China, 2002; Volume 2, pp. 21–25. Land 2023, 12, 265 15 of 16 42. Morgan, R.P.C.; Quinton, J.N.; Smith, R.E.; Govers, G.; Poesen, J.W.A.; Auerswald, K.; Chisci, G.; Torri, D.; Styczen, M.E. The European soil erosion model (EUROSEM): A process-based approach for predicting sediment transport from fields and small catchments. Earth Surf. Process. Landf. 1998, 23, 527–544. [CrossRef] 43. El Kateb, H.; Zhang, H.; Zhang, P.; Mosandl, R. Soil erosion and surface runoff on different vegetation covers and slope gradients: A field experiment in Southern Shaanxi Province, China. Catena 2013, 105, 1–10. [CrossRef] 44. Wang, S.; Feng, X.; Wang, Y.; Zheng, Z.; Li, T.; He, S.; Zhang, X.; Yu, H.; Huang, H.; Liu, T.; et al. Characteristics of nitrogen loss in sloping farmland with purple soil in southwestern China during maize (Zea mays L.) growth stages. Catena 2019, 182, 104169. [CrossRef] 45. Liu, B.; Nearing, M.A.; Risse, L.M. Slope gradient effects on soil loss for steep slopes. Trans. ASABE 1994, 37, 1835–1840. [CrossRef] 46. Peng, X.; Dai, Q.; Ding, G.; Li, C. Role of underground leakage in soil, water and nutrient loss from a rock-mantled slope in the karst rocky desertification area. J. Hydrol. 2019, 578, 124086. [CrossRef] 47. Lal, R. Soil degradative effects of slope length and tillage method on alfisols in western Nigeria. II. Soil chemical properties, plant nutrient loss and water quality. Land Degrad. Dev. 1997, 8, 221–244. [CrossRef] 48. Xing, W.; Yang, P.; Ren, S.; Ao, C.; Li, X.; Gao, W. Slope length effects on processes of total nitrogen loss under simulated rainfall. Catena 2016, 139, 73–81. [CrossRef] 49. Duan, L.; Huang, M.; Zhang, L. Differences in hydrological responses for different vegetation types on a steep slope on the Loess Plateau, China. J. Hydrol. 2016, 537, 356–366. [CrossRef] 50. Suo, L.; Huang, M. Stochastic modelling of soil water dynamics and sustainability for three vegetation types on the Chinese Loess Plateau. Soil Res. 2019, 57, 500–512. [CrossRef] 51. Liu, Q.; Zhu, B.; Tang, J.; Huang, W.; Zhang, X. Hydrological processes and sediment yields from hillslope croplands of regosol under different slope gradients. Soil Sci. Soc. Am. J. 2017, 81, 1517–1525. [CrossRef] 52. Kuznetsova, A.; Brockhoff, P.B.; Christensen, R.H.B. lmerTest package: Tests in linear mixed effects models. J. Stat. Softw. 2017, 82, 1–26. [CrossRef] 53. Breiman, L. Random forests. Mach. Learn 2001, 45, 5–32. [CrossRef] 54. Archer, E. Estimate Permutation p-Values for Random Forest Importance Metrics; R Package Version 2.5.1. Available online: https://github.com/EricArcher/rfPermute (accessed on 11 October 2022). 55. Lefcheck, J.S. PIECEWISESEM: Piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol. Evol. 2016, 7, 573–579. [CrossRef] 56. Zhuang, H.; Wang, Y.; Liu, H.; Wang, S.; Zhang, W.; Zhang, S.; Dai, Q. Large-scale soil erosion estimation considering vegetation growth cycle. Land 2021, 10, 473. [CrossRef] 57. Zhou, J.; Fu, B.; Gao, G.; Lu, Y.; Liu, Y.; Lu, N.; Wang, S. Effects of precipitation and restoration vegetation on soil erosion in a semi-arid environment in the Loess Plateau, China. Catena 2016, 137, 1–11. [CrossRef] 58. Tao, W.; Wang, Q.; Guo, L.; Lin, H. A new analytical model for predicting soil erosion and nutrient loss during crop growth on the Chinese loess plateau. Soil Tillage Res. 2020, 199, 104585. [CrossRef] 59. Zhao, J.; Yang, Y.; Zhang, K.; Jeong, J.; Zeng, Z.; Zang, H. Does crop rotation yield more in China? A meta-analysis. Field Crops Res. 2020, 245, 107659. [CrossRef] 60. Lehtinen, T.; Schlatter, N.; Baumgarten, A.; Bechini, L.; Krüger, J.; Grignani, C.; Zavattaro, L.; Costamagna, C.; Spiegel, H. Effect of crop residue incorporation on soil organic carbon and greenhouse gas emissions in European agricultural soils. Soil Use Manag. 2014, 30, 524–538. [CrossRef] 61. Lu, F. How can straw incorporation management impact on soil carbon storage? A meta-analysis. Mitig. Adapt. Strat. Gl. 2015, 20, 1545–1568. [CrossRef] 62. Powlson, D.S.; Riche, A.B.; Coleman, K.; Glendining, M.J.; Whitmore, A.P. Carbon sequestration in European soils through straw incorporation: Limitations and alternatives. Waste Manag. 2008, 28, 741–746. [CrossRef] 63. Osborne, S.L.; Johnson, J.M.F.; Jin, V.L.; Hammerbeck, A.L.; Varvel, G.E.; Schumacher, T.E. The impact of corn residue removal on soil aggregates and particulate organic matter. Bioenergy Res. 2014, 7, 559–567. [CrossRef] 64. Smil, V. Phosphorus in the environment: Natural flows and human interferences. Annu. Rev. Environ. Resour. 2000, 25, 53–88. [CrossRef] 65. Alewell, C.; Ringeval, B.; Ballabio, C.; Robinson, D.A.; Panagos, P.; Borrelli, P. Global phosphorus shortage will be aggravated by soil erosion. Nat. Commun. 2020, 11, 4546. [CrossRef] 66. Jiao, P.J.; Xu, D.; Wang, S.L.; Zhang, T.Q. Phosphorus loss by surface runoff from agricultural field plots with different cropping systems. Nutr. Cycl. Agroecosyst. 2011, 90, 23–32. [CrossRef] 67. Khokhar, A.; Yousuf, A.; Singh, M.; Sharma, V.; Sandhu, P.S.; Chary, G.R. Impact of land configuration and strip-intercropping on runoff, soil loss and crop yields under rainfed conditions in the Shivalik Foothills of North-West, India. Sustainability 2021, 13, 6282. [CrossRef] 68. Liu, D.; She, D.; Yu, S.E.; Shao, G.; Chen, D. Rainfall intensity and slope gradient effects on sediment losses and splash from a saline–sodic soil under coastal reclamation. Catena 2015, 128, 54–62. [CrossRef] 69. Guo, Z.; Ma, M.; Cai, C.; Wu, Y. Combined effects of simulated rainfall and overland flow on sediment and solute transport in hillslope erosion. J. Soil Sediments 2018, 18, 1120–1132. [CrossRef] Land 2023, 12, 265 16 of 16 70. Kumar, S.; Mishra, A. Critical Erosion Area Identification Based on Hydrological Response Unit Level for Effective Sedimentation Control in a River Basin. Water Resour. Manag. 2015, 29, 1749–1765. [CrossRef] 71. Liu, B.Y.; Nearing, M.A.; Shi, P.J.; Jia, Z.W. Slope length effects on soil loss for steep slopes. Soil Sci. Soc. Am. J. 2000, 64, 1759–1763. [CrossRef] 72. Sadeghi, S.H.R.; Seghaleh, M.B.; Rangavar, A.S. Plot sizes dependency of runoff and sediment yield estimates from a small watershed. Catena 2013, 102, 55–61. [CrossRef] 73. van de Giesen, N.; Stomph, T.J.; de Ridder, N. Surface runoff scale effects in West African watersheds: Modeling and management options. Agric. Water Manag. 2005, 72, 109–130. [CrossRef] 74. Stomph, T.J.; De Ridder, N.; Van de Giesen, N.C. A flume design for the study of slope length effects on runoff. Earth Surf. Proc. Land 2001, 26, 647–655. [CrossRef] 75. Young, R.A.; Olness, A.E.; Mutchler, C.K.; Moldenhauer, W.C. Chemical and physical enrichments of sediment from cropland. Trans. ASABE 1986, 29, 165–169. [CrossRef] 76. Prats, S.A.; Wagenbrenner, J.W.; Martins, M.A.S.; Malvar, M.C.; Keizer, J.J. Mid-term and scaling effects of forest residue mulching on post-fire runoff and soil erosion. Sci. Total Environ. 2016, 573, 1242–1254. [CrossRef] [PubMed] 77. Wu, S.B.; Chen, L.; Wang, N.L.; Zhang, J.; Wang, S.Q.; Bagarello, V.; Ferro, V. Variable scale effects on hillslope soil erosion during rainfall-runoff processes. Catena 2021, 207, 105606. [CrossRef] 78. Regan, J.T.; Rodgers, M.; Healy, M.G.; Kirwan, L.; Fenton, O. Determining phosphorus and sediment release rates from five Irish tillage soils. J. Environ. 2010, 39, 185–192. [CrossRef] 79. Brennan, R.B.; Fenton, O.; Grant, J.; Healy, M.G. Impact of chemical amendment of dairy cattle slurry on phosphorus, suspended sediment and metal loss to runoff from a grassland soil. Sci. Total Environ. 2011, 409, 5111–5118. [CrossRef] 80. Ziadat, F.M.; Taimeh, A.Y. Effect of rainfall intensity, slope, land use and antecedent soil moisture on soil erosion in an arid environment. Land Degrad. Dev. 2013, 24, 582–590. [CrossRef] 81. Rodriguez-Blanco, M.L.; Taboada-Castro, M.M.; Taboada-Castro, M.T. Rainfall-runoff response and event-based runoff coefficients in a humid area (northwest Spain). Hydrolog. Sci. J. 2012, 57, 445–459. [CrossRef] 82. Huo, J.; Yu, X.; Liu, C.; Chen, L.; Zheng, W.; Yang, Y.; Tang, Z. Effects of soil and water conservation management and rainfall types on runoff and soil loss for a sloping area in North China. Land Degrad. Dev. 2020, 31, 2117–2130. [CrossRef] 83. Zheng, H.J.; Nie, X.F.; Liu, Z.; Mo, M.H.; Song, Y.J. Identifying optimal ridge practices under different rainfall types on runoff and soil loss from sloping farmland in a humid subtropical region of Southern China. Agric. Water Manag. 2021, 255, 107043. [CrossRef] 84. Kato, T.; Kuroda, H.; Nakasone, H. Runoff characteristics of nutrients from an agricultural watershed with intensive livestock production. J. Hydrol. 2009, 368, 79–87. [CrossRef] 85. Bai, Y.H.; Zha, X.; Zhang, J.; Chen, S.F. The threshold of nitrogen and phosphorus loss in runoff on degraded Ferralsols of Fujian province, southern China. Environ. Earth Sci. 2020, 79, 395. [CrossRef] Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Land Multidisciplinary Digital Publishing Institute

Effects of Crop Rotation and Topography on Soil Erosion and Nutrient Loss under Natural Rainfall Conditions on the Chinese Loess Plateau

Land , Volume 12 (2) – Jan 17, 2023

Loading next page...
 
/lp/multidisciplinary-digital-publishing-institute/effects-of-crop-rotation-and-topography-on-soil-erosion-and-nutrient-FTvrZsQNmt
Publisher
Multidisciplinary Digital Publishing Institute
Copyright
© 1996-2023 MDPI (Basel, Switzerland) unless otherwise stated Disclaimer Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. Terms and Conditions Privacy Policy
ISSN
2073-445X
DOI
10.3390/land12020265
Publisher site
See Article on Publisher Site

Abstract

land Article Effects of Crop Rotation and Topography on Soil Erosion and Nutrient Loss under Natural Rainfall Conditions on the Chinese Loess Plateau 1 , 2 3 2 , Chenhui Li , Wenhai Shi and Mingbin Huang * College of Resources and Environment, Northwest A & F University, Xianyang 712100, China State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Northwest A & F University, Xianyang 712100, China School of Water and Environment, Chang’an University, Xi’an 710054, China * Correspondence: hmbd@nwsuaf.edu.cn Abstract: Erosive rainfall results in the loss of both soil and nutrients, which indirectly triggers soil deterioration and a reduction in land productivity. However, how rainfall affects runoff, soil erosion, and nutrient loss under different crop rotation patterns and topographic factors remains unclear. This experiment observed nine runoff-erosion plots on the Chinese Loess Plateau (CLP) from 2019 to 2020 to determine the effects of crop type, rotation pattern, and slope gradient and length on runoff, soil erosion, and nutrient loss. Runoff, soil erosion, and nutrient loss were highest for the fallow plots; values for these variables for spring corn and winter wheat plots were not significantly different. Crop rotation generated greater runoff, soil erosion, and nutrient loss compared to non-rotation. Soil erosion and associated nutrient loss increased, but not significantly, with slope for gradients of 0.5 , 1 , and 3 , while runoff and associated nutrient loss did not increase. In addition, soil erosion and associated nutrient loss were significantly greater for slope lengths of 20 m vs. 50 m. A structural equation model showed rainfall characteristics significantly impacted runoff and soil erosion and subsequently affected nutrient loss. This study increases the understanding of runoff, soil erosion, and nutrient loss from cropland with gentle slopes on the CLP. Keywords: crop rotation; runoff and sediment; nitrogen and phosphorus losses; Chinese Loess Plateau Citation: Li, C.; Shi, W.; Huang, M. Effects of Crop Rotation and Topography on Soil Erosion and Nutrient Loss under Natural Rainfall 1. Introduction Conditions on the Chinese Loess Soil erosion driven by runoff and sediment generally takes nutrients away from soil, Plateau. Land 2023, 12, 265. https:// reduces soil fertility, and then causes soil degradation [1,2]. Therefore, soil erosion has doi.org/10.3390/land12020265 attracted widespread attention as a worldwide problem [3–5], especially in severe soil Academic Editor: Vincent Chaplot erosion areas, such as the Chinese Loess Plateau (CLP). The CLP is the most erosive area in the world, generating approximately 90% of the sediment and nutrients that enter the Received: 17 December 2022 Yellow River [6], as well as a large amount of soil nitrogen and phosphorus that cause Revised: 6 January 2023 Accepted: 14 January 2023 water pollution and eutrophication [7]. Soil nutrients are critical to the productivity of Published: 17 January 2023 terrestrial ecosystems and have received a large amount of attention [8,9]. Some research shows inappropriate land use and tillage practices can increase soil nutrient loss [10], which usually occurs via two mechanisms—sediment-bound nutrients lost as particles and soluble nutrients lost with runoff [11]. Runoff and soil erosion from farmland are considered the Copyright: © 2023 by the authors. main sources of soil nutrients entering water bodies [12,13]. Therefore, investigating the soil Licensee MDPI, Basel, Switzerland. erosion and nutrient loss from farmland is important for better management of agricultural This article is an open access article production and the ecological environment. distributed under the terms and Soil nutrient loss processes are complex and controlled by various factors, such as conditions of the Creative Commons rainfall, soil tillage system, and topography [14–16]. The main driver of soil erosion in Attribution (CC BY) license (https:// semi-arid systems is rainfall, with characteristics such as duration, amount, and intensity creativecommons.org/licenses/by/ controlling hydrological processes [17,18]. Rainfall characteristics are also the main factors 4.0/). Land 2023, 12, 265. https://doi.org/10.3390/land12020265 https://www.mdpi.com/journal/land Land 2023, 12, 265 2 of 16 that affect soil nutrient loss and have been extensively studied [19,20]. The sediment yield rate decreases with increasing rainfall duration [21], and even the form of nutrient loss varies with rainfall duration [22]. Sediment concentration and soil loss increase as rainfall amount increases [23]. In general, heavy rainfalls of short duration and high intensity cause the most erosion [24] and also drive the redistribution of nutrients [25]. Rainfall reaching an intensity of 15 mm/h can cause severe soil erosion [26]. However, previous studies have mainly focused on simulated rainfall [27–29] and long-term natural rainfall experiments are relatively rare [30,31]. Therefore, soil nutrient migration in farmland ecosystems needs to be characterized under natural precipitation conditions. Crop type and rotation are other major factors that affect soil erosion and nutrient loss from farmland [32,33]. The direct impact of cultivating different crops on soil loss and water erosion has been widely studied. Carroll et al. (1997) show that, in central Queensland, the average annual runoff from wheat is lower and the soil loss higher than that from sorghum and sunflower [34]. Prasuhn (2012) demonstrates that, among 203 crop fields in the Swiss Midlands, the greatest total soil erosion occurs in winter wheat fields, followed by potato and fallow fields [35]. Based on 198 rainfall-runoff events, Fiene and Auerswald (2007) show the obvious difference in soil loss based on the previous season crop [36]. However, inadequate information is available about the effect of winter wheat and spring corn (common crops on the CLP) on soil and nutrient loss in surface runoff. Moreover, research on crop rotation has mainly focused on crop yield and soil moisture [37,38], while relatively little has been carried out to study soil erosion and nutrient loss in the CLP region [33,39]. Therefore, conditions of soil erosion and nutrient loss under common crop types and rotation need to be clarified to inform the development of more effective cropping systems to replace those in current use. Topographic factors such as slope gradient and length also affect soil erosion and nutrient loss. In some soil erosion models (such as USLE, CSLE, and EUROSEM), slope gradient and length are important parameters [40–42]. Many studies show the amount of soil erosion and associated nutrient loss are positively related to the gradient of culti- vated land [43,44]. The occurrence of soil erosion will especially increase when crops are cultivated on farmland with steep slopes [45]. However, some reports indicate total nitro- gen and total phosphorus in surface runoff decrease with increasing slope gradient [46]. Slope length also plays a key role in soil erosion and nutrient loss. Lal (1997) shows total runoff nutrient losses decrease with increasing slope length under conventional tillage practices, whereas the opposite result is obtained under no-till practices [47]. Xing et al. (2016) indicate that sediment-associated total nitrogen loss increases with slope length under simulated rainfall conditions [48]. Overall, studies to date show slope gradient and length can have either positive or negative impacts on runoff, erosion, and nutrient loss and have mainly focused on soil erosion and nutrient loss for steep slopes. The influence of slope gradient and length on soil erosion from gently sloping agricultural land has not received sufficient attention, and studies are needed to facilitate the development of new technologies for erosion control on gently sloping agricultural fields. Overall, the influence of rainfall, crop type and rotation, and topography on soil erosion and nutrient loss under gently sloping agricultural land in the CLP is not well understood. Therefore, the objectives of this study were to: (1) identify the runoff, soil erosion, and nutrient loss for different crops and crop rotation systems; (2) explain the effects of slope gradient and length on runoff, soil erosion, and nutrient loss for gently sloping agricultural land; (3) investigate how rainfall factors affect runoff, soil erosion, and nutrient loss. The results of this study can clarify the nutrient loss pathways of runoff and sediment from cropland and may also provide a theoretical basis for solving problems caused by nutrient loss. Land 2023, 12, x FOR PEER REVIEW 3 of 16 Land 2023, 12, 265 3 of 16 off and sediment from cropland and may also provide a theoretical basis for solving problems caused by nutrient loss. 2. Materials and Methods 2. Materials and Methods 2.1. Study Site 2.1. Study Site 0 00 The study was conducted in the Wangdonggou watershed (latitude: 107 40 30 – The study was conducted in the Wangdonggou watershed (latitude: 107°40′30″– 0 00  0 00  0 00 107 42 30 E; longitude: 35 12 16 –35 16 00 N; altitude: 946–1226 m) at the Changwu 107°42′30″ E; longitude: 35°12′16″–35°16′00″ N; altitude: 946–1226 m) at the Changwu Agro-ecological Experiment Station, Chinese Academy of Sciences (CAS) and Ministry Agro-ecological Experiment Station, Chinese Academy of Sciences (CAS) and Ministry of Water Resources (MWR), in Shaanxi Province, China. The watershed is a typical of Water Resources (MWR), in Shaanxi Province, China. The watershed is a typical loess 2 2 loess tableland-gully region with an area of 8.3 km and a gully density of 2.78 km/km 2 2 tableland-gully region with an area of 8.3 km and a gully density of 2.78 km/km (Fig- (Figure 1). The gradient of tableland varies from 0–3 . The study area has a warm-temperate ure 1). The gradient of tableland varies from 0–3°. The study area has a warm-temperate semi-humid continental monsoon climate with an average annual temperature of 9.2 C semi-humid continental monsoon climate with an average annual temperature of 9.2 °C (1957–2014) [49]. The average annual precipitation is 579 mm (1960–2016), with more than (1957–2014) [49]. The average annual precipitation is 579 mm (1960–2016), with more 70% falling from April to September [50]. The soil in this watershed is dominated by Heilu than 70% falling from April to September [50]. The soil in this watershed is dominated soil [51], and the parent material is medium loamy Malan loess. The thickness of the loess by Heilu soil [51], and the parent material is medium loamy Malan loess. The thickness layer is more than 100 m, and the soil in the whole profile is uniform and loose with good of the loess layer is more than 100 m, and the soil in the whole profile is uniform and air permeability. loose with good air permeability. Figure 1. Location of the study site and experimental design. (a) the Loess Plateau of China, (b) Wangdonggou watershed on the Loess Plateau of China, (c) runoff and sediment collecting equipment, and (d) nine runoff plots. Numbers of 1–9 show the locations of all nine runoff plots. Land 2023, 12, 265 4 of 16 2.2. Experimental Design Nine loess experimental plots were arranged at Changwu station with three crop/land use types: winter wheat (Triticum aestivum L.), spring corn (Zea mays L.) and fallow (Table 1). Winter wheat was sown in September and harvested in early June; spring corn was sown in April and harvested in September. The boundaries of each plot were constructed by reinforced concrete baffles with 35 cm buried and 15 cm above ground to separate runoff from inside and outside the plots. A self-recording rainfall gauge was placed around the experimental plots to record rainfall. Each plot had a set of equipment for collecting runoff and sediment, including a water collection tank, water conveyance tank, water-retaining dike, and runoff barrel. Table 1. Basic information about the study plots. Plot 1 2 3 4 5 6 7 8 9 Slope length (m) 20 20 20 50 50 50 20 20 20 Slope gradient ( ) 0.5 0.5 0.5 0.5 0.5 0.5 0.5 1 3 Crop rotation Yes Yes Yes Yes Yes Yes No No No 2019 SC WW FA SC WW FA WW WW WW Crop/land type 2020 FA SC WW FA SC WW WW WW WW 2019 78.0 80.0 47.3 74.7 78.0 48.0 79.3 76.7 78.0 Coverage (%) 2020 67.3 75.0 45.0 69.7 73.3 45.7 76.7 75.0 76.7 Note: FA: Fallow; SC: Spring corn; WW: Winter wheat. Coverage was measured in October of 2019 and 2020, respectively. The coverage of winter wheat fields took into account the presence of weeds and stubble at that time. Plots 1 to 3 had a slope length of 20 m and gradient of 0.5 with a wheat-corn-fallow rotation treatment; plots 4 to 6 had a slope length of 50 m and gradient of 0.5 with a wheat-corn-fallow rotation treatment; plots 7 to 9 had a slope length of 20 m and gradients of 0.5 , 1 , and 3 , respectively, with a winter wheat monoculture (Figure 1). Table 1 shows crop types planted during the study period (2019 and 2020). Data from plots 2 and 7 in 2019 and plots 3 and 7 in 2020 for winter wheat were used to compare the effects of rotation and no-rotation on runoff, soil erosion, and nutrient loss. Data from plots 1 and 4 in 2019 and plots 2 and 5 in 2020 for spring corn, data from plots 2 and 5 in 2019 and plots 3 and 6 in 2020 for winter wheat, and data from plots 3 and 6 in 2019 and plots 1 and 4 in 2020 for fallow were used to analyze the effects of different crop/land uses on runoff, soil erosion, and nutrient loss. Data from plots 7, 8, and 9 were used to compare the effect of slope gradient on runoff, soil erosion, and nutrient loss with precipitation as a random variable. Data from plots 1 to 6 were used to contrast the effects of slope length and crop type on runoff, soil erosion, and nutrient loss. 2.3. Data Collection Data were collected from May to October of 2019 and 2020 because most of the rainfall and runoff on the CLP occurs during these months. The total rainfall, rainfall duration, average rainfall intensity (ARI), and the maximum rainfall intensity (MRI) for each rainfall event were obtained from the self-recording rainfall gauge and are given in Table 2. The runoff volume of each plot was obtained by measuring the height of the inside of the runoff barrel according to the known barrel bottom area after the occurrence of erosive rainfall. The suspension (runoff + sediment) in the runoff barrel was thoroughly stirred and samples were collected using 5000-mL plastic bottles. Some samples were filtered through medium-speed phosphorus-free filter paper, and approximately 300-mL samples were immediately refrigerated (4 C) and brought back to the laboratory for determination of total nitrogen and total phosphorus in 24 h. The total nitrogen of water samples was measured by alkaline potassium persulfate ablation UV spectrophotometry and the total phosphorus of water samples was measured by ammonium molybdate spectrophotometry. The remaining samples were separated by filtering through filter paper after settling for Land 2023, 12, 265 5 of 16 24 h. Total sediment was calculated by air drying sediment samples and weighing them to obtain the sediment sample mass. Subsequently, the sediment samples were brought back to the laboratory for determination of total nitrogen and total phosphorus. Soil total nitrogen was measured by the semi-micro Kjeltec method, and soil total phosphorus was measured by the HClO –H SO elimination method. The vegetation cover, plant spacing, 4 2 4 and height of the plots were measured, and soil samples were collected in September and October 2019 for determination of their physical and chemical properties. The basic physical and chemical indicators are shown in Table 3. Table 2. Erosive rainfall characteristics during the study period. Date Rainfall Rainfall Duration ARI MRI (d/m/y) (mm) (h) (mm/h) (mm/h) 05/06/2019 56.5 18.5 3.1 6.8 28/06/2019 49.3 28.2 1.8 3.0 23/07/2019 28.6 4.0 6.3 7.2 29/07/2019 34.2 4.7 7.3 9.8 04/08/2019 20.5 13.3 1.5 3.7 26/08/2019 71.6 12.5 5.7 25.2 15/09/2019 70.0 8.0 3.6 8.8 19/06/2020 58.5 53.0 1.1 3.6 26/06/2020 66.0 17.3 3.8 14.8 11/07/2020 41.4 24.0 1.7 5.0 26/07/2020 47.0 32.5 1.4 4.6 07/08/2020 31.0 6.0 5.2 7.8 13/08/2020 20.8 9.5 2.2 4.6 19/08/2020 88.3 40.0 2.2 6.0 Note: MRI: maximum rainfall intensity; ARI: average rainfall intensity. Table 3. Soil properties of experimental plots. SOC TN TP AP AK pH BD Clay Silt Sand Soil Texture Plot (g/kg) (g/kg) (g/kg) (mg/kg) (mg/kg) – (g/cm ) (%) (%) (%) – 1 9.33 0.98 0.99 14.51 139.82 7.93 1.21 27.92 37.1 34.98 Loamy clay 2 8.67 1.02 0.87 10.01 125.66 7.95 1.28 27.31 37.28 35.41 Loamy clay 3 8.07 0.97 0.81 7.87 117.56 7.96 1.27 27.78 37.42 34.8 Loamy clay 4 7.78 0.91 0.99 10.85 112.84 8.08 1.25 28.08 39.04 32.88 Loamy clay 5 8.28 0.96 0.97 13.19 136.62 7.96 1.24 27.75 37.35 34.9 Loamy clay 6 8.15 0.98 0.94 7.91 134.93 7.79 1.48 26.63 36.77 36.6 Loamy clay 7 9.03 1.11 0.86 7.57 117.73 7.77 1.29 26.3 37.76 35.94 Loamy clay 8 8.85 1.07 0.84 11.46 107.95 7.72 1.38 26.9 39.42 33.68 Loamy clay 9 8.1 0.97 0.84 9.38 113.35 7.74 1.3 25.29 36.29 38.42 Loamy clay Note: SOC: soil organic carbon; TN: total nitrogen; TP: total phosphorus; AP: available phosphorus; AK: available potassium; BD: bulk density. 2.4. Data Analysis All statistical analyses were performed using R version 4.0.2. The mixed linear model and multiple comparisons based on Tukey’s test were used to identify the differences in the effects of tillage practices (rotation and no-rotation), crop type (spring corn, winter wheat, and fallow), and topography (slope gradient and length) on runoff, soil erosion, and nutrient loss; significance differences were judged based on the R package lmerTest [52]. The main factors influencing nutrient loss were decided using a random forest model [53], and the significance of each predictor was assessed using the R package rfPermute [54]. A structural equation model was constructed in combination with professional knowledge using R package piecewiseSEM [55], which was used to quantify the direct and indirect effects of explanatory variables on response variables. Land 2023, 12, x FOR PEER REVIEW 6 of 16 winter wheat, and fallow), and topography (slope gradient and length) on runoff, soil erosion, and nutrient loss; significance differences were judged based on the R package lmerTest [52]. The main factors influencing nutrient loss were decided using a random forest model [53], and the significance of each predictor was assessed using the R pack- age rfPermute [54]. A structural equation model was constructed in combination with professional knowledge using R package piecewiseSEM [55], which was used to quantify the direct and indirect effects of explanatory variables on response variables. Land 2023, 12, 265 6 of 16 3. Results 3.1. Runoff, Soil Erosion, and Nutrient Loss for Different Crops/Land Uses 3. Results Runoff, soil erosion, and nutrient loss for the spring corn, winter wheat, and fallow 3.1. Runoff, Soil Erosion, and Nutrient Loss for Different Crops/Land Uses plots are compared in Figure 2. The fallow plots had more runoff, soil erosion, and nu- Runoff, soil erosion, and nutrient loss for the spring corn, winter wheat, and fallow trient loss than the spring corn or winter wheat plots. Specifically, the fallow plots were plots are compared in Figure 2. The fallow plots had more runoff, soil erosion, and significantly different from the spring corn and winter wheat plots in terms of runoff nutrient loss than the spring corn or winter wheat plots. Specifically, the fallow plots were dept significantly h, soil erosion, different tot fraom l phosphorus the spring cor in n runo and winter ff (RTP wheat ), toplots tal nit inrterms ogen in of r so unof il flos depth, s (STN), soil erosion, total phosphorus in runoff (RTP), total nitrogen in soil loss (STN), and total and total phosphorus in soil loss (STP). For total nitrogen in runoff (RTN), the fallow phosphorus in soil loss (STP). For total nitrogen in runoff (RTN), the fallow plots were plots were only significantly different from the spring corn plots and not from winter only significantly different from the spring corn plots and not from winter wheat plots wheat plots (Figure 2c). Spring corn plots had a lower runoff depth and associated RTP (Figure 2c). Spring corn plots had a lower runoff depth and associated RTP and RTN but and RTN but higher soil loss and associated STP and STN compared to winter wheat higher soil loss and associated STP and STN compared to winter wheat plots; however, the plots; however, the differences were not significant. Overall, runoff, soil erosion, and differences were not significant. Overall, runoff, soil erosion, and nutrient loss from the nutrient loss from the fallow plots were significantly higher than for spring corn and fallow plots were significantly higher than for spring corn and winter wheat plots, but the winter wheat plots, but the differences between the latter two crops were not significant. differences between the latter two crops were not significant. Figure 2. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 2. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil (STP) for spring corn, winter wheat, and fallow plots. Values followed by different letters above loss (STP) for spring corn, winter wheat, and fallow plots. Values followed by different letters a column are significantly different at p  0.05. Different letters indicate significant differences above a column are significantly different at p ≤ 0.05. Different letters indicate significant differ- between treatments. ences between treatments. 3.2. Runoff, Soil Erosion, and Nutrient Loss for Different Rotation Systems The runoff, soil erosion, and nutrient loss for winter wheat cropped under rotation and no-rotation systems are compared in Figure 3. Runoff depth, soil erosion, RTN, and STP were significantly higher under the rotation system than under the no-rotation system in 2019; only RTP was significantly higher in 2020. Combining data from the two years, the rotation system resulted in significantly greater runoff depth, RTN, and Land 2023, 12, x FOR PEER REVIEW 7 of 16 3.2. Runoff, Soil Erosion, and Nutrient Loss for Different Rotation Systems The runoff, soil erosion, and nutrient loss for winter wheat cropped under rotation and no-rotation systems are compared in Figure 3. Runoff depth, soil erosion, RTN, and STP were significantly higher under the rotation system than under the no-rotation sys- Land 2023, 12, 265 7 of 16 tem in 2019; only RTP was significantly higher in 2020. Combining data from the two years, the rotation system resulted in significantly greater runoff depth, RTN, and STN than the no-rotation system by 72.74, 99.41, and 182.06%, respectively. For soil erosion, STN than the no-rotation system by 72.74, 99.41, and 182.06%, respectively. For soil RTP, and STP, the rotation system yielded higher values than the no-rotation system, erosion, RTP, and STP, the rotation system yielded higher values than the no-rotation but the differences were not significant. Overall, the use of a rotation system for winter system, but the differences were not significant. Overall, the use of a rotation system wheat resulted in more runoff, soil erosion, and nutrient loss compared to a no-rotation for winter wheat resulted in more runoff, soil erosion, and nutrient loss compared to a system. no-rotation system. Figure 3. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 3. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss (STP) in rotation and no-rotation systems for 2019 and 2020. ALL: both years (2019 and 2020); *, loss (STP) in rotation and no-rotation systems for 2019 and 2020. ALL: both years (2019 and 2020); p < 0.05; **, p < 0.01; *** p < 0.001. *, p < 0.05; **, p < 0.01; *** p < 0.001. 3.3. Impact of Slope Gradient and Length on Runoff, Soil Erosion, and Nutrient Loss 3.3. Impact of Slope Gradient and Length on Runoff, Soil Erosion, and Nutrient Loss A comparison of runoff, soil erosion, and nutrient loss for gradients of 0.5 , 1 , and 3 A comparison of runoff, soil erosion, and nutrient loss for gradients of 0.5°, 1°, and is illustrated in Figure 4. Soil loss and associated STP and STN losses increased with the 3° is illustrated in Figure 4. Soil loss and associated STP and STN losses increased with gradient but differences between the various slopes were not significant. Runoff depth and the gradient but differences between the various slopes were not significant. Runoff associated RTP and RTN losses did not increase with slope gradient. In general, slope had depth and associated RTP and RTN losses did not increase with slope gradient. In gen- no significant effect on soil erosion and nutrient loss for the gentle slope range considered era (0.5–3 l, slope h ). ad no significant effect on soil erosion and nutrient loss for the gentle slope range considered (0.5–3°). Runoff, soil erosion, and nutrient loss for different slope lengths for the three land use types considered are shown in Figure 5. Plots with a slope length of 50 m had lower runoff depth, soil erosion, RTN, STN, RTP, and STP than plots with a slope length of 20 m for all crop/land use types, which implies runoff, soil erosion, and nutrient loss de- Land 2023, 12, x FOR PEER REVIEW 8 of 16 crease with slope length. For spring corn and winter wheat, runoff depth, soil erosion, RTN, STN, and STP were significantly higher for the 20- vs. 50-m slope length; RTP was also higher for the 20-m slope, but the difference was not significant. For fallow plots, runoff depth and RTN were significantly lower for the 50-m slope; soil loss, RTP, STN, and STP were also lower, but the differences were not significant. When all crop/land use types were considered, increasing slope length resulted in significant decreases in Land 2023, 12, 265 8 of 16 runoff depth, RTN, and STN and insignificant decreases in soil loss, RTP, and STP. These results indicate that slope length significantly affected runoff and nitrogen loss. Figure 4. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 4. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss loss (STP) for slope gradients of 0.5°  , 1°, and 3  °. Values followed by different letters above a column (STP) for slope gradients of 0.5 , 1 , and 3 . Values followed by different letters above a column are are significantly different at p ≤ 0.05. significantly different at p  0.05. Runoff, soil erosion, and nutrient loss for different slope lengths for the three land use types considered are shown in Figure 5. Plots with a slope length of 50 m had lower runoff depth, soil erosion, RTN, STN, RTP, and STP than plots with a slope length of 20 m for all crop/land use types, which implies runoff, soil erosion, and nutrient loss decrease with slope length. For spring corn and winter wheat, runoff depth, soil erosion, RTN, STN, and STP were significantly higher for the 20- vs. 50-m slope length; RTP was also higher for the 20-m slope, but the difference was not significant. For fallow plots, runoff depth and RTN were significantly lower for the 50-m slope; soil loss, RTP, STN, and STP were also lower, but the differences were not significant. When all crop/land use types were considered, increasing slope length resulted in significant decreases in runoff depth, RTN, and STN and insignificant decreases in soil loss, RTP, and STP. These results indicate that slope length significantly affected runoff and nitrogen loss. Land 2023, 12, x FOR PEER REVIEW 9 of 16 Land 2023, 12, 265 9 of 16 Figure 5. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total Figure 5. The (a) runoff depth (RDe), (b) soil loss (SL), (c) total nitrogen in runoff (RTN), (d) total phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil phosphorus in runoff (RTP), (e) total nitrogen in soil loss (STN), and (f) total phosphorus in soil loss loss (STP) for slope lengths of 20 and 50 m. The significance value in the upper left corner repre- (STP) for slope lengths of 20 and 50 m. The significance value in the upper left corner represents the sents the results of the analysis of all crops. *, p < 0.05; **, p < 0.01; *** p < 0.001. results of the analysis of all crops. *, p < 0.05; **, p < 0.01; *** p < 0.001. 3.4. Impact of Natural Rainfall Characteristics on Runoff, Soil Erosion, and Nutrient Loss 3.4. Impact of Natural Rainfall Characteristics on Runoff, Soil Erosion, and Nutrient Loss The influence of total rainfall, rainfall duration, MRI, ARI, soil condition, runoff depth, The influence of total rainfall, rainfall duration, MRI, ARI, soil condition, runoff and soil loss on RTN, RTP, STN, and STP losses were studied using a random forest model depth, and soil loss on RTN, RTP, STN, and STP losses were studied using a random (Figure 6). The pathways and coefficients affecting runoff, soil erosion, and nutrient loss forest model (Figure 6). The pathways and coefficients affecting runoff, soil erosion, and are presented in a structural equation model for each crop/land use type (Figure 7). The nutrient loss are presented in a structural equation model for each crop/land use type results show rainfall was the main driving force related to the occurrence of runoff and soil erosion, which indirectly leads to nutrient loss. For runoff, the structural equation model (Figure 7). The results show rainfall was the main driving force related to the occurrence shows the total rainfall and rainfall duration significantly altered runoff depth in the spring of runoff and soil erosion, which indirectly leads to nutrient loss. For runoff, the struc- corn plots (path coefficients of 0.714 and 0.335, respectively; Figure 7a), while the total tural equation model shows the total rainfall and rainfall duration significantly altered rainfall only significantly influenced runoff depth in the winter wheat and fallow plots runoff depth in the spring corn plots (path coefficients of 0.714 and −0.335, respectively; (path coefficients of 0.593 and 0.668, respectively; Figure 7b,c). Importantly, MRI did not Figure 7a), while the total rainfall only significantly influenced runoff depth in the winter significantly affect runoff depth for any of the three crop/land use types. Soil loss was wheat and fallow plots (path coefficients of 0.593 and 0.668, respectively; Figure 7b,c). significantly affected by runoff depth in spring corn and fallow plots (path coefficients of Importantly, MRI did not significantly affect runoff depth for any of the three crop/land 0.403 and 0.845, respectively; Figure 7a,c), and significantly influenced by MRI in the winter use types. Soil loss was significantly affected by runoff depth in spring corn and fallow wheat plots (path coefficient of 0.404; Figure 7b). plots (path coefficients of 0.403 and 0.845, respectively; Figure 7a,c), and significantly in- fluenced by MRI in the winter wheat plots (path coefficient of 0.404; Figure 7b). Land 2023, 12, x FOR PEER REVIEW 10 of 16 Land 2023, 12, x FOR PEER REVIEW 10 of 16 Land 2023, 12, 265 10 of 16 Figure 6. Random forest mean predictor importance (% increase in mean square error) of the fac- Figure Figure 6. 6. Random Random forest mean predictor im forest mean predictor importance portance (% (% incr increase in m ease in mean ean squar squ eaerr re error) or) of the of the factors fac- tors influencing (a) total nitrogen in runoff (RTN), (b) total phosphorus in runoff (RTP), (c) total tors influencing (a) total nitrogen in runoff (RTN), (b) total phosphorus in runoff (RTP), (c) total influencing (a) total nitrogen in runoff (RTN), (b) total phosphorus in runoff (RTP), (c) total nitrogen nitrogen in soil loss (STN), and (d) total phosphorus in soil loss (STP) for spring corn, (e–h) winter nitrogen in soil loss (STN), and (d) total phosphorus in soil loss (STP) for spring corn, (e–h) winter in soil loss (STN), and (d) total phosphorus in soil loss (STP) for spring corn, (e–h) winter wheat, and wheat, and (i–l) fallow plots, respectively. TR: total rainfall; SC: soil condition (including soil nu- wheat, and (i–l) fallow plots, respectively. TR: total rainfall; SC: soil condition (including soil nu- (i–l) fallow plots, respectively. TR: total rainfall; SC: soil condition (including soil nutrients, bulk trients, bulk density, and saturated hydraulic conductivity); RDe: runoff depth; RDu: rainfall du- trients, bulk density, and saturated hydraulic conductivity); RDe: runoff depth; RDu: rainfall du- density, and saturated hydraulic conductivity); RDe: runoff depth; RDu: rainfall duration; MRI: ration; MRI: maximum rainfall intensity; ARI: average rainfall intensity; SL: soil erosion. *, p < 0.05; ration; MRI: maximum rainfall intensity; ARI: average rainfall intensity; SL: soil erosion. *, p < 0.05; maximum rainfall intensity; ARI: average rainfall intensity; SL: soil erosion. *, p < 0.05; **, p < 0.01. **, p < 0.01. **, p < 0.01. Figure Figure 7. 7. Structural Structural equati equation on mod models els shshow ow the the pathways pathways and coeff and coef icificients ents affecting affecting runoff runof , sofi,l e soil ro- Figure 7. Structural equation models show the pathways and coefficients affecting runoff, soil ero- sion, and nutrient loss. Insignificant paths (p ≥ 0.05) are semi-transparent dashed lines, while sig- erosion, and nutrient loss. Insignificant paths (p  0.05) are semi-transparent dashed lines, while sion, and nutrient loss. Insignificant paths (p ≥ 0.05) are semi-transparent dashed lines, while sig- nificant paths (p < 0.05) are shown with solid lines. Black and red arrows indicate positive and significant paths (p < 0.05) are shown with solid lines. Black and red arrows indicate positive and nificant paths (p < 0.05) are shown with solid lines. Black and red arrows indicate positive and negative pathways, respectively. Conditional R denotes the proportion of variance explained. The negative negative pathways, respe pathways, respectively ctively. Conditio . Conditional nal R R denote denotes s the the proportion of proportion of variance variance explained. The explained. The standardized path coefficient represents the extent of the effect of the independent variable on the standardized path coefficient represents the extent of the effect of the independent variable on the standardized path coefficient represents the extent of the effect of the independent variable on the dependent variable. Significant pathways are marked by: *, p < 0.05; **, p < 0.01; ***, p < 0.001. dependent variable. Significant pathways are marked by: *, p < 0.05; **, p < 0.01; ***, p < 0.001. dependent variable. Significant pathways are marked by: *, p < 0.05; **, p < 0.01; ***, p < 0.001. Land 2023, 12, 265 11 of 16 For nutrient loss, the random forest model results indicate total rainfall, MRI, ARI, runoff depth, and soil erosion all had important effects (Figure 6). The results of structural equation modeling further demonstrate rainfall indirectly causes nutrient loss by affecting runoff depth and soil erosion (Figure 7). Runoff depth explains 83, 71, and 91% of the total variation in RTN loss and 49, 44, and 91% of the total variation in RTP loss for spring corn, winter wheat, and fallow plots, respectively. Soil erosion explains 97, 94, and 100% of the total variation in STN loss and 96, 97, and 96% of the total variation in STP loss for spring corn, winter wheat, and fallow plots, respectively (Figure 7). However, some differences were evident among three plots. For the fallow plots (Figure 7c), runoff depth and soil loss significantly impacted RTN, RTP, STN, and STP losses. For the winter wheat plots, runoff depth only had a significantly positive effect on RTN and RTP losses (path coefficients of 0.905 and 0.554, respectively; Figure 7b), while soil erosion had a significantly positive effect on STN and STP losses (path coefficients of 0.938 and 0.999, respectively; Figure 7b). For spring corn plots, runoff depth only had a significantly positive effect on RTN loss (path coefficient of 0.877; Figure 7a), while soil loss had a significantly negative effect on RTP loss and a significantly positive effect on STN and STP losses (path coefficients of 0.421, 0.938, and 0.999, respectively; Figure 7a). 4. Discussion 4.1. Impact of Crop Type on Runoff, Soil Erosion, and Nutrient Loss Our results indicate natural rainfall causes greater runoff, soil erosion, and nutrient loss in fallow plots compared to spring corn and winter wheat plots (Figure 2). Many studies also report that the choice of crop can effectively reduce surface runoff and soil and nutrient losses [34,35]. In the study region, the growing season of spring corn is from April to September, which aligns with the rainy season [56,57]. The coverage of spring corn can reduce soil erosion by intercepting rainfall and decreasing raindrop kinetic energy [57]. In addition, the active root system can increase erosion resistance. However, the growing season of winter wheat is from September to May of the next year, which does not align with the rainy season. Therefore, runoff, soil erosion, and nutrient losses are higher but not significantly greater in spring corn vs. winter wheat fields [58]. In our study, the spring corn and winter wheat plots both had greater biomass and coverage than the fallow plot, which likely resulted in the reduced soil erosion and nutrient loss. 4.2. Impact of Crop Rotation on Runoff, Soil Erosion, and Nutrient Loss Crop rotation is an important management practice [59]. Our results show rotation can result in greater runoff, soil erosion, and nutrient loss than no-rotation management in winter wheat plots (Figure 3). Crop rotation changes the coverage time and duration on the soil surface as well as the root distribution in the topsoil. In addition, farmers use different methods to harvest winter wheat and spring corn. Plots in the no-rotation system are left as wheat stubble after harvest, which provides better surface cover and roughness and reduces runoff and soil erosion [60,61]. Moreover, wheat residuals can increase soil organic matter content, improve soil physical characteristics, and subsequently reduce runoff, soil erosion, and nutrient loss [62,63]. For RTP, the major fraction of soil phosphorus is tightly adsorbed to mineral particles and bound within organic matter, so only a very small fraction leaches as dissolved soluble phosphate [64,65]. This might be the reason for the insignificant difference in phosphorus loss noted in our results for the different rotation systems. In this study, the results show that crop rotation generates greater runoff, which is different from the conclusion of Jiao et al. (2011) [66]. Jiao et al. (2011) compares surface runoff for three double cropping systems of wheat/maize, wheat/cotton, and wheat/soybean with wheat/fallow and indicates that all double cropping systems increases the ground cover and significantly reduces runoff [66,67]. However, our result is drawn for the same crop type and coverage of winter wheat with rotation and no-rotation. Land 2023, 12, 265 12 of 16 4.3. Impact of Slope Gradient and Length on Runoff, Soil Erosion, and Nutrient Loss Sloping farmland is one of the main types of agricultural land on the CLP. The gradient affects runoff, soil erosion, and nutrient loss by altering runoff flow rates and infiltration rates [68,69]. Our results show soil erosion and associated STP and STN losses increased with slope gradient, but the differences among the three gentle slope gradients (0.5 , 1 , and 3 ) were not significant (Figure 4b–f). Runoff depth and associated RTP and RTN losses did not increase with slope gradient, perhaps because the gentle gradients considered in this study were not sufficiently large to create stronger effects [68,70]. Slope length also impacts runoff, soil erosion, and nutrient loss [48,71]. Our study results show all indicators for the 20-m slope length were greater than those for the 50-m slope length. Among all indicators, runoff depth, RTN and STN were significantly larger for the 20-m slope length than for the 50-m slope length when all crop/land use types were considered (Figure 5a,c,e); these differences were not significant for soil erosion, RTP, and STP (Figure 5b,d,f). These findings confirm previous results reported by Xing et al. (2016) [48], where both runoff and associated total nitrogen loss rate decreased with increasing slope length. This is attributed to the longer slope increasing surface runoff infiltration and deposition of eroded soil along its length [72,73], which subsequently reduces runoff-associated RTN and RTP losses and sediment-associated STN and STP losses [74,75]. Furthermore, some studies suggest soil loss decreases with slope length according to a power-law trend [76,77], which might explain why the soil loss was less for the 50-m vs. 20-m slope. However, the two years of data from this study show significant nitrogen loss from runoff and soil but insignificant phosphorus loss for the different slope lengths. This might be due to RTN values being much higher than RTP values (Figure 5c,d). The latter were usually insignificant, likely due to phosphorus being predominantly and firmly associated with soil particle surfaces [78,79]; the amount of STP supports this suggestion. 4.4. Impact of Rainfall on Runoff, Soil Erosion, and Nutrient Loss The results of our study show rainfall characteristics are the dominant driver of runoff, soil erosion, and associated nutrient loss (Figures 6 and 7). This result is consistent with those of many other studies [14,20,80]. Runoff depth was mainly positively associated with total rainfall and negatively associated with rainfall duration (Figure 7), which were in agreement with earlier studies [81,82]. These results suggest total rainfall and rainfall duration determine the runoff depth, which indirectly causes runoff-associated nutrient loss. Our results are also consistent with assertions that runoff and MRI are the main factors influencing soil loss [83]. This is because MRI can reflect rainfall impact on soil particle detachment, whereas surface runoff depth can be used to represent the transport of sediments [11]. Our study also found rainfall indirectly causes nutrient loss by affecting runoff depth and soil erosion, with nutrient loss due to both N and P dissolved in runoff as well as partially bonded to soil particles [11,84]. The forms of nutrient loss also differed for the three crops/land uses studied (Figure 7). This result aligns with those of previous studies that show the intensity of rainfall can impact the form of nutrient loss [85]. At higher intensity the nutrient loss is attributed more to soil nutrients being transported with sediments, while at lower intensity the nutrient loss is attributed more soluble elements carried with runoff. Overall, the responses of runoff, soil erosion, and nutrient loss to rainfall intensity in plots of different crops differed and, thus, resulted in different forms of nutrient loss. 5. Conclusions This paper reports on an experiment featuring nine runoff-erosion plots on the CLP studied in 2019 to 2020 to determine the effect of crop type, rotation pattern, slope gradient, and slope length on runoff, soil erosion, and nutrient loss. Runoff, soil erosion, and nutrient loss were higher in fallow plots than in spring corn and winter wheat plots, with no significant differences noted between the latter two. Plots subject to crop rotation Land 2023, 12, 265 13 of 16 generated greater runoff, soil erosion, and nutrient loss compared to those managed with a no-rotation system. Soil loss and associated nutrient loss tended to increase with slope gradient, while runoff and associated nutrient loss did not show any trends. Runoff, soil erosion, and associated nutrient loss were significantly greater in plots with a slope length of 20- vs. 50-m under gentle gradient conditions. Total rainfall and rainfall duration were the primary factors determining runoff, soil erosion, and associated nutrient loss, and maximum rainfall intensity partially explained soil loss. Our results suggest that the effect of rainfall characteristics on soil nutrient loss varies with rotation system and crop type. These results are useful for developing targeted farmland nutrient conservation strategies. Author Contributions: Data collection, data analysis, original draft preparation and writing, C.L.; manuscript of validation and revision, W.S.; manuscript review, editing, and project supervision, M.H. All authors have read and agreed to the published version of the manuscript. Funding: This study was supported by the National Natural Science Foundation of China (42107351); the China Postdoctoral Science Foundation (2019M663917XB); the Natural Science Fund of Shaanxi Province (2021JQ-227). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: The vector boundary of the Wangdonggou watershed was downloaded from the Loess Plateau Science Data Centre, National Earth System Science Data Sharing Infrastruc- ture, National Science and Technology Infrastructure of China (http://loess.geodata.cn (accessed on 18 March 2021). The rest of the data are available upon request from the corresponding author. Acknowledgments: The authors thank the members of the Changwu Agri-Ecological Experimental Station, Chinese Academy of Sciences, for their assistance. Conflicts of Interest: The authors declare no conflict of interest. References 1. Lal, R. Soil erosion impact on agronomic productivity and environment quality. Crit. Rev. Plant Sci. 1998, 17, 319–464. [CrossRef] 2. Gomes, L.; Simoes, S.J.C.; Dalla Nora, E.L.; de Sousa-Neto, E.R.; Forti, M.C.; Ometto, J.P.H.B. Agricultural expansion in the Brazilian Cerrado: Increased soil and nutrient losses and decreased agricultural productivity. Land 2019, 8, 12. [CrossRef] 3. Liu, C.; Sui, J.; Wang, Z. Sediment load reduction in Chinese rivers. Int. J. Sediment Res. 2008, 23, 44–55. [CrossRef] 4. Ochoa, P.A.; Fries, A.; Mejia, D.; Burneo, J.I.; Ruiz-Sinoga, J.D.; Cerda, A. Effects of climate, land cover and topography on soil erosion risk in a semiarid basin of the Andes. Catena 2016, 140, 31–42. [CrossRef] 5. Martinez-Mena, M.; Carrillo-Lopez, E.; Boix-Fayos, C.; Almagro, M.; Franco, N.G.; Diaz-Pereira, E.; Montoya, I.; de Vente, J. Long-term effectiveness of sustainable land management practices to control runoff, soil erosion, and nutrient loss and the role of rainfall intensity in Mediterranean rainfed agroecosystems. Catena 2020, 187, 104352. [CrossRef] 6. Shao, F.; Tao, W.; Wang, Q.; Wu, J.; Su, L.; Sun, Y.; Zeng, S. A modified model for predicting nutrient loss in runoff using a time-varying mixing layer. J. Hydrol. 2021, 603, 127091. [CrossRef] 7. Chen, M.; Chen, J.; Sun, F. Agricultural phosphorus flow and its environmental impacts in China. Sci. Total Environ. 2008, 405, 140–152. [CrossRef] 8. Tan, Z.; Lal, R.; Wiebe, K.D. Global soil nutrient depletion and yield reduction. J. Sustain. Agric. 2005, 26, 123–146. [CrossRef] 9. Sharpley, A.; Kleinman, P. Effect of rainfall simulator and plot scale on overland flow and phosphorus transport. J. Environ. Qual. 2003, 32, 2172–2179. [CrossRef] 10. Fu, B.; Chen, L.; Ma, K.; Zhou, H.; Wang, J. The relationships between land use and soil conditions in the hilly area of the loess plateau in northern Shaanxi, China. Catena 2000, 39, 69–78. [CrossRef] 11. Shi, W.; Huang, M.; Wu, L. Prediction of storm-based nutrient loss incorporating the estimated runoff and soil loss at a slope scale on the Loess Plateau. Land Degrad. Dev. 2018, 29, 2899–2910. [CrossRef] 12. Bechmann, M.E.; Boe, F. Soil tillage and crop growth effects on surface and subsurface runoff, loss of soil, phosphorus and nitrogen in a cold climate. Land 2021, 10, 77. [CrossRef] 13. Guo, S.; Zhai, L.; Liu, J.; Liu, H.; Chen, A.; Wang, H.; Wu, S.; Lei, Q. Cross-ridge tillage decreases nitrogen and phosphorus losses from sloping farmlands in southern hilly regions of China. Soil Tillage Res. 2019, 191, 48–56. [CrossRef] 14. Zheng, F. Effects of accelerated soil erosion on soil nutrient loss after deforestation on the Loess Plateau. Pedosphere 2005, 15, 707–715. 15. Stevens, C.J.; Quinton, J.N.; Bailey, A.P.; Deasy, C.; Silgram, M.; Jackson, D.R. The effects of minimal tillage, contour cultivation and in-field vegetative barriers on soil erosion and phosphorus loss. Soil Tillage Res. 2009, 106, 145–151. [CrossRef] Land 2023, 12, 265 14 of 16 16. Zhang, X.; Hu, M.; Guo, X.; Yang, H.; Zhang, Z.; Zhang, K. Effects of topographic factors on runoff and soil loss in Southwest China. Catena 2018, 160, 394–402. [CrossRef] 17. Weltzin, J.F.; Tissue, D.T. Resource pulses in arid environments—Patterns of rain, patterns of life. New Phytol. 2003, 157, 171–173. [CrossRef] 18. Yu, Y.; Wei, W.; Chen, L.; Feng, T.; Daryanto, S. Quantifying the effects of precipitation, vegetation, and land preparation techniques on runoff and soil erosion in a Loess watershed of China. Sci. Total Environ. 2019, 652, 755–764. [CrossRef] 19. Jung, J.W.; Lim, S.S.; Kwak, J.H.; Park, H.J.; Yoon, K.S.; Kim, H.Y.; Baek, W.J.; Choi, W.J. Further Understanding of the impacts of rainfall and agricultural management practices on nutrient loss from rice paddies in a monsoon area. Water Air Soil Pollut. 2015, 226, 283. [CrossRef] 20. Lemma, B.; Kebede, F.; Mesfin, S.; Fitiwy, I.; Abraha, Z.; Norgrove, L. Quantifying annual soil and nutrient lost by rill erosion in continuously used semiarid farmlands, North Ethiopia. Environ. Earth Sci. 2017, 76, 190. [CrossRef] 21. Ma, L.; Li, J.; Liu, J. Effects of antecedent soil water content on infiltration and erosion processes on loessial slopes under simulated rainfall. Hydrol. Res. 2020, 51, 882–893. [CrossRef] 22. He, Z.; Weng, H.; Wu, T.; Ho, H.; Lin, Y.; Ran, Q. Impact of moving rainfall events on hillslope pollutant transport. Environ. Earth Sci. 2015, 74, 5989–5999. [CrossRef] 23. Nishigaki, T.; Sugihara, S.; Kilasara, M.; Funakawa, S. Surface runoff generation and soil loss under different soil and rainfall properties in the Uluguru Mountains, Tanzania. Land Degrad. Dev. 2017, 28, 283–293. [CrossRef] 24. Zheng, M.; Cai, Q.; Cheng, Q. Modelling the runoff-sediment yield relationship using a proportional function in hilly areas of the Loess Plateau, North China. Geomorphology 2008, 93, 288–301. 25. Fiener, P.; Dlugoss, V.; Van Oost, K. Erosion-induced carbon redistribution, burial and mineralisation—Is the episodic nature of erosion processes important? Catena 2015, 133, 282–292. 26. Foster, G.R.; Young, R.A.; Römkens, M.J.M.; Onstad, C.A. Processes of Soil Erosion by Water. In Soil Erosion and Crop Productivity; American Society of Agronomy, Inc.: Madison, WI, USA; Crop Science Society of America, Inc.: Madison, WI, USA; Soil Science Society of America, Inc.: Madison, WI, USA, 1985; Volume 9, pp. 137–162. 27. Shi, Z.; Fang, N.; Wu, F.; Wang, L.; Yue, B.; Wu, G. Soil erosion processes and sediment sorting associated with transport mechanisms on steep slopes. J. Hydrol. 2012, 454, 123–130. [CrossRef] 28. Cerda, A.; Gonzalez-Pelayo, O.; Gimenez-Morera, A.; Jordan, A.; Pereira, P.; Novara, A.; Brevik, E.C.; Prosdocimi, M.; Mahmood- abadi, M.; Keesstra, S.; et al. Use of barley straw residues to avoid high erosion and runoff rates on persimmon plantations in Eastern Spain under low frequency-high magnitude simulated rainfall events. Soil Res. 2016, 54, 154–165. [CrossRef] 29. Keesstra, S.D.; Rodrigo-Comino, J.; Novara, A.; Gimenez-Morera, A.; Pulido, M.; Di Prima, S.; Cerda, A. Straw mulch as a sustainable solution to decrease runoff and erosion in glyphosate-treated clementine plantations in Eastern Spain. An assessment using rainfall simulation experiments. Catena 2019, 174, 95–103. [CrossRef] 30. Hussein, M.H.; Kariem, T.H.; Othman, A.K. Predicting soil erodibility in northern Iraq using natural runoff plot data. Soil Tillage Res. 2007, 94, 220–228. [CrossRef] 31. Huo, J.; Liu, C.; Yu, X.; Chen, L.; Zheng, W.; Yang, Y.; Yin, C. Direct and indirect effects of rainfall and vegetation coverage on runoff, soil loss, and nutrient loss in a semi-humid climate. Hydrol. Process. 2021, 35, e13985. [CrossRef] 32. Martinez, J.R.F.; Zuazo, V.H.D.; Raya, A.M. Environmental impact from mountainous olive orchards under different soil- management systems (SE Spain). Sci. Total Environ. 2006, 358, 46–60. [CrossRef] [PubMed] 33. Yan, L.; Xue, L.; Petropoulos, E.; Qian, C.; Hou, P.; Xu, D.; Yang, L. Nutrient loss by runoff from rice-wheat rotation during the wheat season is dictated by rainfall duration. Environ. Pollut. 2021, 205, 117382. [CrossRef] [PubMed] 34. Carroll, C.; Halpin, M.; Burger, P.; Bell, K.; Sallaway, M.M.; Yule, D.F. The effect of crop type, crop rotation, and tillage practice on runoff and soil loss on a Vertisol in central Queensland. Aust. J. Soil Res. 1997, 35, 925–939. [CrossRef] 35. Prasuhn, V. On-farm effects of tillage and crops on soil erosion measured over 10 years in Switzerland. Soil Tillage Res. 2012, 120, 137–146. [CrossRef] 36. Fiene, P.; Auerswald, K. Rotation effects of potato, maize, and winter wheat on soil erosion by water. Soil Sci. Soc. Am. J. 2007, 71, 1919–1925. [CrossRef] 37. Zhang, Y.; Wang, R.; Wang, H.; Wang, S.; Wang, X.; Li, J. Soil water use and crop yield increase under different long-term fertilization practices incorporated with two-year tillage rotations. Agric. Water Manag. 2019, 221, 362–370. [CrossRef] 38. Cui, Z.; Yan, B.; Gao, Y.; Wu, B.; Wang, Y.; Xie, Y.; Xu, P.; Wang, H.; Wen, M.; Wang, Y.; et al. Crop yield and water use efficiency in response to long-term diversified crop rotations. Front. Plant. Sci 2022, 13, 1024898. [CrossRef] 39. Shen, Y.; Li, L.; Chen, W.; Robertson, M.; Unkovich, M.; Bellotti, W.; Probert, M. Soil water, soil nitrogen and productivity of lucerne-wheat sequences on deep silt loams in a summer dominant rainfall environment. Field Crop Res. 2009, 111, 97–108. [CrossRef] 40. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Index Losses: A Guide to Conservation Planning. In USDA Agricultural Handbook No. 537; US Department of Agriculture: Washington, DC, USA, 1981. 41. Liu, B.; Zhang, K.; Xie, Y. An empirical soil loss equation. In Process of Soil Erosion and Its Environment Effect, Proceedings of the 12th International Soil Conservation Organization Conference, Beijing, China, 26–31 May 2002; Tsinghua University Press: Beijing, China, 2002; Volume 2, pp. 21–25. Land 2023, 12, 265 15 of 16 42. Morgan, R.P.C.; Quinton, J.N.; Smith, R.E.; Govers, G.; Poesen, J.W.A.; Auerswald, K.; Chisci, G.; Torri, D.; Styczen, M.E. The European soil erosion model (EUROSEM): A process-based approach for predicting sediment transport from fields and small catchments. Earth Surf. Process. Landf. 1998, 23, 527–544. [CrossRef] 43. El Kateb, H.; Zhang, H.; Zhang, P.; Mosandl, R. Soil erosion and surface runoff on different vegetation covers and slope gradients: A field experiment in Southern Shaanxi Province, China. Catena 2013, 105, 1–10. [CrossRef] 44. Wang, S.; Feng, X.; Wang, Y.; Zheng, Z.; Li, T.; He, S.; Zhang, X.; Yu, H.; Huang, H.; Liu, T.; et al. Characteristics of nitrogen loss in sloping farmland with purple soil in southwestern China during maize (Zea mays L.) growth stages. Catena 2019, 182, 104169. [CrossRef] 45. Liu, B.; Nearing, M.A.; Risse, L.M. Slope gradient effects on soil loss for steep slopes. Trans. ASABE 1994, 37, 1835–1840. [CrossRef] 46. Peng, X.; Dai, Q.; Ding, G.; Li, C. Role of underground leakage in soil, water and nutrient loss from a rock-mantled slope in the karst rocky desertification area. J. Hydrol. 2019, 578, 124086. [CrossRef] 47. Lal, R. Soil degradative effects of slope length and tillage method on alfisols in western Nigeria. II. Soil chemical properties, plant nutrient loss and water quality. Land Degrad. Dev. 1997, 8, 221–244. [CrossRef] 48. Xing, W.; Yang, P.; Ren, S.; Ao, C.; Li, X.; Gao, W. Slope length effects on processes of total nitrogen loss under simulated rainfall. Catena 2016, 139, 73–81. [CrossRef] 49. Duan, L.; Huang, M.; Zhang, L. Differences in hydrological responses for different vegetation types on a steep slope on the Loess Plateau, China. J. Hydrol. 2016, 537, 356–366. [CrossRef] 50. Suo, L.; Huang, M. Stochastic modelling of soil water dynamics and sustainability for three vegetation types on the Chinese Loess Plateau. Soil Res. 2019, 57, 500–512. [CrossRef] 51. Liu, Q.; Zhu, B.; Tang, J.; Huang, W.; Zhang, X. Hydrological processes and sediment yields from hillslope croplands of regosol under different slope gradients. Soil Sci. Soc. Am. J. 2017, 81, 1517–1525. [CrossRef] 52. Kuznetsova, A.; Brockhoff, P.B.; Christensen, R.H.B. lmerTest package: Tests in linear mixed effects models. J. Stat. Softw. 2017, 82, 1–26. [CrossRef] 53. Breiman, L. Random forests. Mach. Learn 2001, 45, 5–32. [CrossRef] 54. Archer, E. Estimate Permutation p-Values for Random Forest Importance Metrics; R Package Version 2.5.1. Available online: https://github.com/EricArcher/rfPermute (accessed on 11 October 2022). 55. Lefcheck, J.S. PIECEWISESEM: Piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol. Evol. 2016, 7, 573–579. [CrossRef] 56. Zhuang, H.; Wang, Y.; Liu, H.; Wang, S.; Zhang, W.; Zhang, S.; Dai, Q. Large-scale soil erosion estimation considering vegetation growth cycle. Land 2021, 10, 473. [CrossRef] 57. Zhou, J.; Fu, B.; Gao, G.; Lu, Y.; Liu, Y.; Lu, N.; Wang, S. Effects of precipitation and restoration vegetation on soil erosion in a semi-arid environment in the Loess Plateau, China. Catena 2016, 137, 1–11. [CrossRef] 58. Tao, W.; Wang, Q.; Guo, L.; Lin, H. A new analytical model for predicting soil erosion and nutrient loss during crop growth on the Chinese loess plateau. Soil Tillage Res. 2020, 199, 104585. [CrossRef] 59. Zhao, J.; Yang, Y.; Zhang, K.; Jeong, J.; Zeng, Z.; Zang, H. Does crop rotation yield more in China? A meta-analysis. Field Crops Res. 2020, 245, 107659. [CrossRef] 60. Lehtinen, T.; Schlatter, N.; Baumgarten, A.; Bechini, L.; Krüger, J.; Grignani, C.; Zavattaro, L.; Costamagna, C.; Spiegel, H. Effect of crop residue incorporation on soil organic carbon and greenhouse gas emissions in European agricultural soils. Soil Use Manag. 2014, 30, 524–538. [CrossRef] 61. Lu, F. How can straw incorporation management impact on soil carbon storage? A meta-analysis. Mitig. Adapt. Strat. Gl. 2015, 20, 1545–1568. [CrossRef] 62. Powlson, D.S.; Riche, A.B.; Coleman, K.; Glendining, M.J.; Whitmore, A.P. Carbon sequestration in European soils through straw incorporation: Limitations and alternatives. Waste Manag. 2008, 28, 741–746. [CrossRef] 63. Osborne, S.L.; Johnson, J.M.F.; Jin, V.L.; Hammerbeck, A.L.; Varvel, G.E.; Schumacher, T.E. The impact of corn residue removal on soil aggregates and particulate organic matter. Bioenergy Res. 2014, 7, 559–567. [CrossRef] 64. Smil, V. Phosphorus in the environment: Natural flows and human interferences. Annu. Rev. Environ. Resour. 2000, 25, 53–88. [CrossRef] 65. Alewell, C.; Ringeval, B.; Ballabio, C.; Robinson, D.A.; Panagos, P.; Borrelli, P. Global phosphorus shortage will be aggravated by soil erosion. Nat. Commun. 2020, 11, 4546. [CrossRef] 66. Jiao, P.J.; Xu, D.; Wang, S.L.; Zhang, T.Q. Phosphorus loss by surface runoff from agricultural field plots with different cropping systems. Nutr. Cycl. Agroecosyst. 2011, 90, 23–32. [CrossRef] 67. Khokhar, A.; Yousuf, A.; Singh, M.; Sharma, V.; Sandhu, P.S.; Chary, G.R. Impact of land configuration and strip-intercropping on runoff, soil loss and crop yields under rainfed conditions in the Shivalik Foothills of North-West, India. Sustainability 2021, 13, 6282. [CrossRef] 68. Liu, D.; She, D.; Yu, S.E.; Shao, G.; Chen, D. Rainfall intensity and slope gradient effects on sediment losses and splash from a saline–sodic soil under coastal reclamation. Catena 2015, 128, 54–62. [CrossRef] 69. Guo, Z.; Ma, M.; Cai, C.; Wu, Y. Combined effects of simulated rainfall and overland flow on sediment and solute transport in hillslope erosion. J. Soil Sediments 2018, 18, 1120–1132. [CrossRef] Land 2023, 12, 265 16 of 16 70. Kumar, S.; Mishra, A. Critical Erosion Area Identification Based on Hydrological Response Unit Level for Effective Sedimentation Control in a River Basin. Water Resour. Manag. 2015, 29, 1749–1765. [CrossRef] 71. Liu, B.Y.; Nearing, M.A.; Shi, P.J.; Jia, Z.W. Slope length effects on soil loss for steep slopes. Soil Sci. Soc. Am. J. 2000, 64, 1759–1763. [CrossRef] 72. Sadeghi, S.H.R.; Seghaleh, M.B.; Rangavar, A.S. Plot sizes dependency of runoff and sediment yield estimates from a small watershed. Catena 2013, 102, 55–61. [CrossRef] 73. van de Giesen, N.; Stomph, T.J.; de Ridder, N. Surface runoff scale effects in West African watersheds: Modeling and management options. Agric. Water Manag. 2005, 72, 109–130. [CrossRef] 74. Stomph, T.J.; De Ridder, N.; Van de Giesen, N.C. A flume design for the study of slope length effects on runoff. Earth Surf. Proc. Land 2001, 26, 647–655. [CrossRef] 75. Young, R.A.; Olness, A.E.; Mutchler, C.K.; Moldenhauer, W.C. Chemical and physical enrichments of sediment from cropland. Trans. ASABE 1986, 29, 165–169. [CrossRef] 76. Prats, S.A.; Wagenbrenner, J.W.; Martins, M.A.S.; Malvar, M.C.; Keizer, J.J. Mid-term and scaling effects of forest residue mulching on post-fire runoff and soil erosion. Sci. Total Environ. 2016, 573, 1242–1254. [CrossRef] [PubMed] 77. Wu, S.B.; Chen, L.; Wang, N.L.; Zhang, J.; Wang, S.Q.; Bagarello, V.; Ferro, V. Variable scale effects on hillslope soil erosion during rainfall-runoff processes. Catena 2021, 207, 105606. [CrossRef] 78. Regan, J.T.; Rodgers, M.; Healy, M.G.; Kirwan, L.; Fenton, O. Determining phosphorus and sediment release rates from five Irish tillage soils. J. Environ. 2010, 39, 185–192. [CrossRef] 79. Brennan, R.B.; Fenton, O.; Grant, J.; Healy, M.G. Impact of chemical amendment of dairy cattle slurry on phosphorus, suspended sediment and metal loss to runoff from a grassland soil. Sci. Total Environ. 2011, 409, 5111–5118. [CrossRef] 80. Ziadat, F.M.; Taimeh, A.Y. Effect of rainfall intensity, slope, land use and antecedent soil moisture on soil erosion in an arid environment. Land Degrad. Dev. 2013, 24, 582–590. [CrossRef] 81. Rodriguez-Blanco, M.L.; Taboada-Castro, M.M.; Taboada-Castro, M.T. Rainfall-runoff response and event-based runoff coefficients in a humid area (northwest Spain). Hydrolog. Sci. J. 2012, 57, 445–459. [CrossRef] 82. Huo, J.; Yu, X.; Liu, C.; Chen, L.; Zheng, W.; Yang, Y.; Tang, Z. Effects of soil and water conservation management and rainfall types on runoff and soil loss for a sloping area in North China. Land Degrad. Dev. 2020, 31, 2117–2130. [CrossRef] 83. Zheng, H.J.; Nie, X.F.; Liu, Z.; Mo, M.H.; Song, Y.J. Identifying optimal ridge practices under different rainfall types on runoff and soil loss from sloping farmland in a humid subtropical region of Southern China. Agric. Water Manag. 2021, 255, 107043. [CrossRef] 84. Kato, T.; Kuroda, H.; Nakasone, H. Runoff characteristics of nutrients from an agricultural watershed with intensive livestock production. J. Hydrol. 2009, 368, 79–87. [CrossRef] 85. Bai, Y.H.; Zha, X.; Zhang, J.; Chen, S.F. The threshold of nitrogen and phosphorus loss in runoff on degraded Ferralsols of Fujian province, southern China. Environ. Earth Sci. 2020, 79, 395. [CrossRef] Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Journal

LandMultidisciplinary Digital Publishing Institute

Published: Jan 17, 2023

Keywords: crop rotation; runoff and sediment; nitrogen and phosphorus losses; Chinese Loess Plateau

There are no references for this article.