Effects of rainfall and slope on runoff,soil erosion and rill development: an experimental study using two loess soils

Runoff generation and soil loss from slopes have been studied for decades, but the relationships among runoff, soil loss and rill development are still not well understood. In this paper, rainfall simulation experiments were conducted in two neighbouring plots (scale: 1 m by 5 m) with four varying slopes (17.6%, 26.8%, 36.4% and 46.6%) and two rainfall intensities (90 and 120 mm h^?1) using two loess soils. Data on rill development were extracted from the digital elevation models by means of photogrammetry. The effects of rainfall intensity and slope gradient on runoff, soil loss and rill development were different for the two soils. The runoff and soil loss from the Anthrosol surface were generally higher than those from the Calcaric Cambisol surface. Higher rainfall intensity produced less runoff and more sediment for almost each treatment. With increasing slope gradient, the values of cumulative runoff and soil loss peaked, except for the treatments with 90 mm h^?1 rainfall on the slopes with Anthrosol. With rainfall duration, runoff discharge decreased for Anthrosol and increased for Calcaric Cambisol for almost all the treatments. For both soils, sediment concentration was very high at the onset of rainfall and decreased quickly. Almost all the sediment concentrations increased on the 17.6% and 26.8% slopes and peaked on the 36.4% and 46.6% slopes. Sediment concentrations were higher on the Anthrosol slopes than on the Calcaric Cambisol slopes. At 90 mm h^?1 rainfall intensity, increasingly denser rills appeared on the Anthrosol slope as the slope gradient increased, while only steep slopes (36.4% and 46.6%) developed rills for the Calcaric Cambisol soil. The contributions of rill erosion ranged from 36% to 62% of the cumulative soil losses for Anthrosol, while the maximum contribution of rill erosion to the cumulative soil loss was only 37.9% for Calcaric Cambisol. Copyright © 2014 John Wiley & Sons, Ltd.     

The effect of slope length on sediment yield by rainfall impact under different land use types

To evaluate the impact of slope length on sediment yield under different rainfall intensities and land use types on low hill gentle slope, the characteristics of sediment yield process were analyzed based on the field artificial rainfall simulation. For the study, grassland and capsicum slope were taken from Anji county of Zhejiang province, China. Results indicated that rainfall intensity had stronger influence than slope length on sediment yield in south region. For capsicum slope, sediment yield increased quickly with increasing slope length when rainfall intensity greater than 90 mm h^–1. The slope length had no significant effect on sediment yield when rainfall intensity less than or equaled to 90 mm h^–1. For grassland, data from experiments indicated that sediment yield increased slowly with increasing slope length under rainfall intensity less than 120 mm h^–1. There was a decreasing tendency of sediment yield at 6 m slope length under all rainfall events. It was concluded from particle size analysis of erosional sediment that silt and clay particles <0.02 mm were always preferentially transported on both capsicum slope (silt 47.1%, clay 40.9%) and grassland (silt 38.3%, clay 35.9%). We hope these results are useful for soil and water conservation and land management.     

The effect of soil type,meteorological forcing and slope gradient on the simulation of internal erosion processes at the local scale

Numerical simulation experiments of water erosion at the local scale (20 × 5 m) using a process‐based model [Plot Soil Erosion Model_2D (PSEM_2D)] were carried out to test the effects of various environmental factors (soil type, meteorological forcing and slope gradient) on the runoff and erosion response and to determine the dominant processes that control the sediment yield at various slope lengths. The selected environmental factors corresponded to conditions for which the model had been fully tested beforehand. The use of a Green and Ampt model for infiltration explained the dominant role played by rainfall intensity in the runoff response. Sediment yield at the outlet of the simulated area was correlated positively with rainfall intensity and slope gradient, but was less sensitive to soil type. The relationship between sediment yield (soil loss per unit area) and slope length was greatly influenced by all environmental factors, but there was a general tendency towards higher sediment yield when the slope was longer. Contribution of rainfall erosion to gross erosion was dominant for all surfaces with slope lengths ranging from 4 to 20 m. The highest sediment yields corresponded to cases where flow erosion was activated. An increase in slope gradient resulted in flow detachment starting upstream. Sediment exported at the outlet of the simulated area came predominantly from the zone located near the outlet. The microrelief helped in the development of a rill network that controlled both the ratio between rainfall and flow erosion and the relationship between sediment yield and slope length. Copyright © 2010 John Wiley & Sons, Ltd.     

Changes of soil surface roughness under water erosion process

The objective of this study was to determine the changing characteristics of soil surface roughness under different rainfall intensities and examine the interaction between soil surface roughness and different water erosion processes. Four artificial management practices (raking cropland, artificial hoeing, artificial digging, and contour tillage) were used according to the local agriculture customs of the Loess Plateau of China to simulate different types of soil surface roughness, using an additional smooth slope for comparison purposes. A total of 20 rainfall simulation experiments were conducted in five 1 m by 2 m boxes under two rainfall intensities (0.68 and 1.50 mm min^?1) on a 15° slope. During splash erosion, soil surface roughness decreased in all treatments except raking cropland and smooth baseline under rainfall intensity of 0.68 mm min^?1, while increasing for all treatments except smooth baseline under rainfall intensity of 1.50 mm min^?1. During sheet erosion, soil surface roughness decreased for all treatments except hoeing cropland under rainfall intensity of 0.68 mm min^?1. However, soil surface roughness increased for the artificial hoeing and raking cropland under rainfall intensity of 1.50 mm min^?1. Soil surface roughness has a control effect on sheet erosion for different treatments under two rainfall intensities. For rill erosion, soil surface roughness increased for raking cropland and artificial hoeing treatments, and soil surface roughness decreased for artificial digging and the contour tillage treatments under two rainfall intensities. Under rainfall intensity of 0.68 mm min^?1, the critical soil surface roughness was 0.706 cm for the resistance control of runoff and sediment yield. Under rainfall intensity of 1.50 mm min^?1, the critical soil surface roughness was 1.633 cm for the resistance control of runoff, while the critical soil surface roughness was 0.706 cm for the resistance control of sediment yield. These findings have important implications for clarifying the erosive nature of soil surface roughness and harnessing sloped farmland. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydraulic conditions for rill incision under simulated rainfall: A laboratory experiment

A series of controlled laboratory experiments were conducted in order to obtain precise data on the hydraulic and sediment transport conditions during rill formation. Tests were carried out using a crusting-prone binary mixed soil in a 15 m long flume at an average slope of 0·087 under simulated rainfall. Rainfall intensities varied from 30–35 mm h^?1 and developed about 70 per cent of the kinetic energy of natural rainfall of similar intensity. Runoff and sediment discharge measured at the downstream weir were strongly influenced by rill forming processes. Essentially, rill incision reduced runoff discharge as a result of increased percolation in rill channels but increased sediment discharge. Secondary entrainment processes, such as bank collapse, also increased sediment discharge at the weir. Knickpoint bifurcation and colluvial deposition, however, decreased sediment discharge. Rills always developed through the formation of a knickpoint. The critical condition for knickpoint initiation was the development of supercritical flow and waves which mould and incise the bed. Prior smoothing of the soil surface by entrainment and redistribution of sediment facilitated supercritical flow. Statistical analysis showed that hydraulic and sediment transport conditions differed significantly in rilled and unrilled flows. The relationship between sediment discharge, rill erosion, and flow hydraulics was found to be nonlinear, conforming to a standard power function in the form y = ax^b. Rills were also associated with significantly increased sediment transport capacities. However, rill initiation was not clearly defined by any specific hydraulic threshold. Instead, rilled and unrilled flows were separated by zones of transition within which both types of flow occur.     

The changing dynamics of rill erosion on sloping farmland during the different growth stages of a maize crop

Rill erosion is a serious concern in the hilly region of China with purple soil, and maize is extensively cultivated in this region. Evaluations of the dynamic mechanisms of rill erosion in sloping farmland areas are particularly important during the maize growing season to determine whether rill erosion can occur. A new ridge tillage (RT) system was designed using local agricultural methods in China. Twelve artificial rainfall experiments were conducted in three 1 × 2 m experimental plots with a slope of 15°, which is a typical slope in the study area. The rainfall intensities were designated as 1.0, 1.5, and 2.0 mm min^?1. The rainfall experiments were performed in the field to determine the characteristics of run‐off and sediment transport related to rill erosion processes during different stages of maize growth and to analyse how hydraulic parameters and the sediment yield of the rill erosion process are related. The results showed that rill flow patterns were mainly classified as subcritical transition flow during all the growth stages of maize. The effects of hydrodynamic parameters on the sediment yield were ordered as follows: Reynolds number > stream power > Froude number > shear stress. The total sediment yield varied by stage as follows: seedling stage > jointing stage > mature stage > tasseling stage. The sediment yield and run‐off rate exhibited a linear relationship that was well described at the hillslope scale. To initiate soil loss in sloping farmland areas with purple soil during the maize growing season, the critical hydrodynamic shear stress and stream power must be at least 46.505 Pa and 1.541 N m^?1 s^?1, respectively.     

Rill network development on loessial hillslopes in China

Rill network development not only potentially affects hillslope and drainage network evolution, but also causes severe soil degradation. However, the studies on rill network development remain inconclusive. This study aimed to investigate the temporal and spatial development of hillslope rill networks and their characteristics based on rainfall simulations and field observations. A soil pan (10.0 m long × 3.0 m wide × 0.5 m deep) on a 20° slope was applied three successive simulated rains at two intensities of 50 and 100 mm h^–1. The field observations were performed on two bare hillslope runoff plots (10.0 m long × 3.0 m wide) at 20°. Three typical erosive natural rainfall events were observed in the field, and rills were measured in detail, similar to the laboratory rainfall simulation. The results indicated that with increases in rainfall events, the rill network morphology varied from incipient formation to the maximum drainage network density. Four rill network development indicators (rill distribution density, distance between rills, rill bifurcation number, and confluence point number) exhibited different changes over time and space. Among the four indicators, the rill bifurcation number was the best indicator for describing rill network development. Rill flow energy increased and decreased cyclically on a slope ranging between ~3 and 4 m. Moreover, rill networks on loessial hillslopes generally evolved into dendritic rather than parallel forms. The development characteristics of the rill network were relatively similar between the laboratory simulation and natural field conditions. Over time, rill erosion control measures become increasingly difficult to implement as the rill network develops. The morphology of eroding rills evolved over time and space, which led to corresponding rill network development. Further study should quantify the impacts of rill network development on soil degradation and land development. © 2020 John Wiley & Sons, Ltd.     

Laboratory experiments of sediment transport from bare soil with a rill

Abstract

Mathematical models developed for quantification of sediment transport in hydrological watersheds require data collected through field or laboratory experiments, but these are still very rare in the literature. This study aims to collect such data at the laboratory scale. To this end, a rainfall simulator equipped with nozzles to spray rainfall was constructed, together with an erosion flume that can be given longitudinal and lateral slopes. Eighty experiments were performed, considering microtopographical features by pre-forming a rill on the soil surface before the start of each experiment. Medium and fine sands were used as soil, and four rainfall intensities (45, 65, 85 and 105 mm h^-1) were applied in the experiments. Rainfall characteristics such as uniformity, granulometry, drop velocity and kinetic energy were evaluated; flow and sediment discharge data were collected and analysed. The analysis shows that the sediment transport rate is directly proportional to rainfall intensity and slope. In contrast, the volumetric sediment concentration stays constant and does not change with rainfall intensity unless the slope changes. These conclusions are restricted to the conditions of experiments performed under rainfall intensities between and 105 mm h^-1 for medium and fine sands in a 136-cm-wide, 650-cm-long and 17-cm-deep erosion flume with longitudinal and lateral slopes varying between 5 and 20%.

Editor Z.W. Kundzewicz; Associate editor G. Mahé

Citation Aksoy, H., Unal, N.E., Cokgor, S., Gedikli, A., Yoon, J., Koca, K., Inci, S.B., Eris, E., and Pak, G., 2013. Laboratory experiments of sediment transport from bare soil with a rill. Hydrological Sciences Journal, 58 (7), 1505–1518.     

An experimental investigation of the effects of thawed soil depth on rill flow velocity

Water flow velocity is an important hydraulic variable in hydrological and soil erosion models, and is greatly affected by freezing and thawing of the surface soil layer in cold high-altitude regions. The accurate measurement of rill flow velocity when impacted by the thawing process is critical to simulate runoff and sediment transport processes. In this study, an electrolyte tracer modelling method was used to measure rill flow velocity along a meadow soil slope at different thaw depths under simulated rainfall. Rill flow velocity was measured using four thawed soil depths (0, 1, 2 and 10 cm), four slope gradients (5°, 10°, 15° and 20°) and four rainfall intensities (30, 60, 90 and 120 mm·h⁻¹). The results showed that the increase in thawed soil depth caused a decrease in rill flow velocity, whereby the rate of this decrease was also diminishing. Whilst the rill flow velocity was positively correlated with slope gradient and rainfall intensity, the response of rill flow velocity to these influencing factors varied with thawed soil depth. The mechanism by which thawed soil depth influenced rill flow velocity was attributed to the consumption of runoff energy, slope surface roughness, and the headcut effect. Rill flow velocity was modelled by thawed soil depth, slope gradient and rainfall intensity using an empirical function. This function predicted values that were in good agreement with the measured data. These results provide the foundation for a better understanding of the effect of thawed soil depth on slope hydrology, erosion and the parameterization scheme for hydrological and soil erosion models.     

Rill length and plot‐scale effects on the hydrogeomorphologic response of gravelly roadbeds

Although unpaved roads are well‐recognized as important sources of Hortonian overland flow and sediment in forested areas, their role in agriculturally‐active rural settings still lacks adequate documentation. In this study, we assessed the effect of micro‐catchment size, slope, and ground cover on runoff and sediment generation from graveled roadbeds servicing a rural area in southern Brazil. Fifteen replications based on 30‐min‐long simulated rainfall experiments were performed at constant rainfall intensities of 22–58 mm h^?1 on roadbeds with varying characteristics including ~3–7 m² micro‐catchment areas, 2–11° slopes, 2–9.7‐m‐long shallow rill features, and 30–100% gravel cover. The contributions of micro‐catchment size and rill length were the most important physical characteristics affecting runoff response and sediment production; both the size of the micro‐catchment and the length of the rills were inversely related to sediment loss and this contradicts most of the rill erosion literature. The effect of micro‐catchment size on runoff and sediment response suggests a potentially problematic spatial‐scale subjectivity of experimental plot results. The inverse relationship between rill length and sediment generation is interpreted here as related to the predominance of coarse fragments within rills, the inability of the shallow flows generated during the simulations to erode this sediment, and their role as zones of net sediment storage. Copyright © 2015 John Wiley & Sons, Ltd.     

Using beryllium‐7 to monitor the relative proportions of interrill and rill erosion from loessal soil slopes in a single rainfall event

Quantifying the relative proportions of soil losses due to interrill and rill erosion processes during erosion events is an important factor in predicting total soil losses and sediment transport and deposition. Beryllium‐7 (⁷Be) can provide a convenient way to trace sediment movement over short timescales providing information that can potentially be applied to longer‐term, larger‐scale erosion processes. We used simulated rainstorms to generate soil erosion from two experimental plots (5 m × 4 m; 25° slope) containing a bare, hand‐cultivated loessal soil, and measured ⁷Be activities to identify the erosion processes contributing to eroded material movement and/or deposition in a flat area at the foot of the slope. Based on the mass balance of ⁷Be detected in the eroded soil source and in the sediments, the proportions of material from interrill and rill erosion processes were estimated in the total soil losses, the deposited sediments in the flat area, and in the suspended sediments discharged from the plots. The proportion of interrill eroded material in the discharged sediment decreased over time as that of rill eroded material increased. The amount of deposited material was greatly affected by overland flow rates. The estimated amounts of rill eroded material calculated using ⁷Be activities were in good agreement with those based on physical measurements of total plot rill volumes. Although time lags of 45 and 11 minutes existed between detection of sediment being removed by rill erosion, based on ⁷Be activities, and observed rill initiation times, our results suggest that the use of ⁷Be tracer has the potential to accurately quantify the processes of erosion from bare, loessal cultivated slopes and of deposition in flatter, downslope areas that occur in single rainfall events. Such measurements could be applied to estimate longer‐term erosion occurring over larger areas possessing similar landforms. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical modeling of gravitational erosion in rill systems

A self-organizing model was developed for simulating rill erosion process on slopes with particular attention to the role of gravitational erosion.For a complete simulation circle,processes such as precipitation,infiltration,runoff,scouring,gravitational erosion and elevation variation were fully considered.Precipitation time(or runoff time) was regarded as iteration benchmark in the model.To specify the contribution of gravitational erosion to the process of rill formation and development,a gravitational erosion module was inserted into the model.Gravitational erosion in rill development was regarded as a Gaussian random process.A model was calibrated by our experimental data,and further validated satisfactorily with 22 runs of experimental results from different investigators. Systematic comparison was made between sediment yields with and without consideration of gravitational erosion module.It was demonstrated that the model could reasonably simulate the rill erosion process under a variety of slope gradients,rainfall intensities and soil conditions upon the gravitational erosion being considered.However,the role of gravitational erosion on sediment yields in rill systems varies significantly under different conditions,although it is of the utmost importance in steeper slopes.The process of gravitational erosion in rill development was studied by a newly-defined parameter a>,which is defined as the volume ratio of gravitational erosion over hydraulic-related erosion.The gravitational contribution to the total erosion could be over 50%for the rill systems with higher rainfall intensity and steeper slopes.     

Flow energy and channel adjustments in rills developed in loamy sand and sandy loam soils

The storms usually associated with rill development in nature are seldom prolonged, so development is often interrupted by interstorm disturbances, e.g. weathering or tillage. In laboratory simulated rainfall experiments, active rill development can be prolonged, and under these conditions typically passes through a period of intense incision, channel extension and bifurcation before reaching quasi‐stable conditions in which little form change occurs. This paper presents laboratory experiments with coarse textured soils under simulated rainfall which show how channel adjustment processes contribute to the evolution of quasi‐stability. Newly incised rills were stabilized for detailed study of links between rill configuration and flow energy. On a loamy sand, adjustment towards equilibrium occurred due to channel widening and meandering, whereas on a sandy loam, mobile knickpoints and chutes, pulsations in flow width and flow depth and changes in stream power and sediment discharge occurred as the channel adjusted towards equilibrium. The tendency of rill systems towards quasi‐stability is shown by changes in stream power values which show short‐lived minima. Differences in energy dissipation in stabilized rills indicate that minimization of energy dissipation was reached locally between knickpoints and at the downstream ends of rills. In the absence of energy gradients in knickpoints and chutes, stabilized rill sections tended toward equilibrium by establishing uniform energy expenditure. The study confirmed that energy dissipation increased with flow aspect ratio. In stabilized rills, flow acceleration reduced energy dissipation on the loamy sand but not on the sandy loam. On both soils flow deceleration tended to increase energy dissipation. Understanding how rill systems evolve towards stability is essential in order to predict how interruptions between storms may affect long‐term rill dynamics. This is essential if event‐based physical models are to become effective in predicting sediment transport on rilled hillslopes under changing weather and climatic conditions. Copyright © 2008 John Wiley and Sons, Ltd.     

Effects of rock fragments incorporated in the soil matrix on concentrated flow hydraulics and erosion

Rock fragments can act as a controlling factor for erosional rates and patterns in the landscape. Thus, the objective of this study is to better understand the role that rock fragments incorporated into the soil matrix play in concentrated flow hydraulics and erosion. Laboratory flume experiments were conducted with soil material that was mixed with rock fragments. Rock fragment content ranged from 0 to 40 per cent by volume. Other treatments were slope (7 and 14%) and flow discharge (5·7 and 11·4 l min^?1). An increase in rock fragment content resulted in lower sediment yield, and broader width of flow. Rock fragment cover at the soil surface, i.e. surface armour, increased with time in experiments with rock fragments. Flow energy was largely dissipated by rock fragment cover. For more turbulent flow conditions, when roughness elements were submerged in the flow, hydraulic roughness was similar for different rock fragment contents. In experiments with few or no rock fragments a narrow rill incised. Flow energy was dissipated by headcuts. Total sediment yield was much larger than for experiments with rock fragments in the soil. Adding just a small number of rock fragments in the soil matrix resulted in a significant reduction of sediment yield. Copyright © 2007 John Wiley & Sons, Ltd.     

RILL EROSION PROCESS AND RILL FLOW HYDRAULIC PARAMETERS

In the rill erosion process, run-on water and sediment from upslope areas, and rill flow hydraulic parameters have significant effects on sediment detachment and transport. However, there is a lack of data to quantify the effects of run-on water and sediment and rill flow hydraulic parameters on rill erosion process at steep hillslopes, especially in the Loess Plateau of China. A dual-box system, consisting of a 2-m-long feeder box and a 5-m-long test box with 26.8% slope gradient was used to quantify the effects of upslope runoff and sediment, and of rill flow hydraulic parameters on the rill erosion process. The results showed that detachment-transport was dominated in rill erosion processes; upslope runoff always caused the net rill detachment at the downslope rill flow channel, and the net rill detachment caused by upslope runoff increased with a decrease of runoff sediment concentration from the feeder box or an increase of rainfall intensity. Upslope runoff discharging into the rill flow channel or an increase of rainfall intensity caused the rill flow to shift from a stratum flow into a turbulent flow. Upslope runoff had an important effect on rill flow hydraulic parameters, such as rill flow velocity, hydraulic radius, Reynolds number, Froude number and the Darcy-Weisbach resistance coefficient. The net rill detachment caused by upslope runoff increased as the relative increments of rill flow velocity, Reynolds number and Froude number caused by upslope runoff increased. In contrast, the net rill detachment decreased with an increase of the relative decrement of the Darcy-Weisbach resistance coefficient caused by upslope runoff. These findings will help to improve the understanding of the effects of run-on water and sediment on the erosion process and to find control strategies to minimize the impact of run-on water.     

Rill erosion processes on a constantly saturated slope

Rill erosion processes on saturated soil slopes are important for understanding erosion hydrodynamics and determining the parameters of rill erosion models. Saturated soil slopes were innovatively created to investigate the rill erosion processes. Rill erosion processes on saturated soil slopes were modelled by using the sediment concentrations determined by sediment transport capacities (STCs) measurement and the sediment concentrations at different rill lengths. Laboratory experiments were performed under varying slope gradients (5°, 10°, 15°, and 20°) and unit-width flow rates (0.33, 0.67, and 1.33 × 10⁻³ m³ s⁻¹ m⁻¹) to measure sediment concentrations at different rill lengths (1, 2, 4, and 8 m) on saturated soil slopes. The measured sediment concentrations along saturated rills ranged from 134.54 to 1,064.47 kg/m³, and also increased exponentially with rill length similar to non-saturated rills. The model of the rill erosion process in non-saturated soil rills was applicable to that in saturated soil rills. However, the sediment concentration of the rill flow increased much faster, with the increase in rill length, to considerably higher levels at STCs. The saturated soil rills produced 120–560% more sediments than the non-saturated ones. Moreover, the former eroded remarkably faster in the beginning section of the rills, as compared with that on the non-saturated soil slopes. This dataset serves as the basis for determining the erosion parameters in the process-based erosion models on saturated soil slopes.     

The transport of aggregates associated with soil organic carbon under the rain-induced overland flow on the Chinese Loess Plateau

Organic carbon (OC) is easily enriched in sediment particles of different sizes due to aggregate breakdown and selective transport for sheet erosion. However, the transport of aggregate-associated OC has not been thoroughly investigated. To address this issue, 27 simulated rainfall experiments were conducted in a 1 m × 0.35 m box on slope gradients of 15°, 10°, and 15°and under three rainfall intensities of 45 mm h⁻¹, 90 mm h⁻¹ and 120 mm h⁻¹. The results showed that OC was obviously enriched in sediment particles of different sizes under sheet erosion. The soil organic carbon (SOC) concentrations of each aggregate size class in sediments were different from those in the original soil, especially when the rainfall intensity or slope was sufficiently low, such as 45 mm h^–1 or 5°, respectively. Under a slope of 5°, the SOC enrichment ratios (ERocs) of small macroaggregates and microaggregates were high but decreased over time. As rainfall intensity increased, OC became enriched in increasingly fine sediment particles. Under a rainfall intensity of 45 mm h^–1, the ERocs of the different aggregate size classes were always high throughout the entire erosion process. Under a rainfall intensity of > 45 mm h^–1 and slope of > 5°, the ERocs of the different aggregate size classes were close to 1.0, especially those of clay and silt. Therefore, the high ERocs in sediments resulted from the first transport of effective clay. Among total SOC loss, the proportion of OC loss caused by the transport of microaggregates and silt plus clay-sized particles was greater than 50%. We also found that low stream power and low water depth were two requirements for the high ERocs in aggregates. Stream power was closely related to sediment particle distribution. Flow velocity was significantly and positively related to the percentage of OC-enriched macroaggregates in the sediments (P > 0.01). Our study will provide important information for understanding the fate of SOC and building physical-based SOC transport models. © 2019 John Wiley & Sons, Ltd.     

A methodology to study rain splash and wash processes under natural rainfall

The complex interactions between rainfall‐driven erosion processes and rainfall characteristics, slope gradient, soil treatment and soil surface processes are not very well understood. A combination of experiments under natural rainfall and a consistent physical theory for their interpretation is needed to shed more light on the underlying processes. The present study demonstrates such a methodology. An experimental device employed earlier in laboratory studies was used to measure downslope rain splash and ‘splash‐creep’, lateral splash, upslope splash and rainfall‐driven runoff transport (wash) from a highly aggregated clay‐rich oxisol exposed to natural rainfall in West Java, Indonesia. Two series of measurements were made: the first with the soil surface at angles of 0°, 5°, 15° and 40°; and the second all at an angle of 5° but with different tillage and mulching treatments. A number of rainfall erosivity indices were calculated from rainfall intensity measurements and compared with measured transport components. Overall storm kinetic energy correlated reasonably well with sediment transport, but much better agreement was obtained when a threshold rainfall intensity (20 mm h^?1) was introduced. Rain splash transport measurements were interpreted using a recently developed theory relating detachment to sediment transport. Furthermore, a conceptually sound yet simple wash transport model is advanced that satisfactorily predicted observed washed sediment concentrations. The lack of replication precluded rigorous assessment of the effect of slope and soil treatment on erosion processes, but some general conclusions could still be drawn. The results stress the importance of experiments under conditions of natural rainfall. Copyright © 2002 John Wiley & Sons, Ltd.     

Sediment delivery on rill and interrill areas

Equations which relate sediment delivery to a power function of flow rate and slope gradient were evaluated in this study. The data used to parameterize the equations were obtained from sites where crop residues had been removed, and moldboard plowing and disking had occurred. Measurements of sediment delivery resulting from simulated rainfall were obtained from preformed rills and interrill areas. The equations provided reliable sediment delivery estimates for selected soils located throughout the United States. To use the sediment delivery equations, soil-related parameter values must be identified. Multiple regression analyses were performed to relate parameter values used in the equations to selected soil properties. Equations were also developed for estimating rill sediment delivery under rainfall conditions from rill soil loss and discharge data collected without the addition of rainfall. The equations identified in this study, and appropriate soils information, can be used to predict sediment delivery on both rill and interrill areas.     

The European Soil Erosion Model (EUROSEM): a dynamic approach for predicting sediment transport from fields and small catchments

The European Soil Erosion Model (EUROSEM) is a dynamic distributed model, able to simulate sediment transport, erosion and deposition over the land surface by rill and interill processes in single storms for both individual fields and small catchments. Model output includes total runoff, total soil loss, the storm hydrograph and storm sediment graph. Compared with other erosion models, EUROSEM has explicit simulation of interill and rill flow; plant cover effects on interception and rainfall energy; rock fragment (stoniness) effects on infiltration, flow velocity and splash erosion; and changes in the shape and size of rill channels as a result of erosion and deposition. The transport capacity of runoff is modelled using relationships based on over 500 experimental observations of shallow surface flows. EUROSEM can be applied to smooth slope planes without rills, rilled surfaces and surfaces with furrows. Examples are given of model output and of the unique capabilities of dynamic erosion modelling in general. © 1998 John Wiley & Sons, Ltd.