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.     

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.     

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.     

DYNAMIC AND KINEMATIC WAVE MODELS FOR 1-D RILL FLOW: A COMPARISON

1INTRODUCTIONInprocess-basedoverlandflowandsoilerosionmodels,surfacerunoffonahillslopeisoftenrepresentedaseitherbroadsheetfloworflowinrillswithassumedrectangularchannelcrosssections(e.g.,Bairdelal1992,NSERL1995).Inmostcasesthehydraulicsofoverlandflowiscalculatedbyusingthekinematicwavemodel,whichisasirnplificationofthedynamicwavemodel(theequationsystemofSaint-Venantequationandequationofcontinuity).Forabroadsheetoverlandflowonhillslopeduetorainfallexcess,theequationofcontinuityiswherehis…     

Rill development and soil erosion: a laboratory study of slope and rainfall intensity

A total of 15 rainfall simulation experiments were conducted in a 1 m by 2 m box varying slope (10, 20, 30%) and rainfall intensity (60, 90, 120 mm h^?1). The experiments were performed to study how rill networks initiate and evolve over time under controlled conditions with regard to the treatment variables considered, and to allow for input in a computer simulation model. Runoff and sediment yield samples were collected. Digital elevation models were calculated by means of photogrammetry for several time steps of most experiments. The soil used in the experiments was a basal till derived Cambisol typical for the Swiss Plateau. While significant differences were found for sediment yield, runoff did not vary significantly with treatment combinations. Increasing rainfall intensity had a larger effect on sediment yield than increasing slope. Rill density and energy expenditure decreased with time, suggesting that energy expenditure was a useful parameter to describe the emergence of rill network at the laboratory scale. Copyright © 2010 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.     

Physically based modelling of sheet erosion (detachment and deposition processes) in complex hillslopes

A one‐dimensional uncoupled model governed by this research is a physics‐based modelling of the rainfall‐runoff induced erosion process. The presented model is composed of three parts of a three‐dimensional (3D) hillslope geometry, a nonlinear storage (kinematic wave) model for hillslope hydrological response, and an unsteady physically based surface erosion model. The 3D hillslope geometry model allows describing of the hillslope morphology by defining their plan shape and profile curvature. By changing these two topographic parameters, nine basic hillslope types are derived. The modelling of hillslope hydrological response is based on a flow continuity equation as the relation of discharge and flow depth is passed on kinematic wave approximation. The erosion model is based on a mass conservation equation for unsteady flow. The model assumes that suspended sediment does not affect flow dynamics. The model also accounts for the effect of flow depth plus loose soil depth on soil detachment. The presented model was run for two different precipitations, slope content, and length, and results were plotted for sediment detachment/deposition rate. Based on the obtained results, in hillslopes with convex and straight profile curvatures, sediment detachment only occurred in the whole length of the hillslope. However, in concave ones, sediment detachment and deposition only occurred together in hillslope. The hillslopes with straight profiles and convergent plans have the highest rate of detachment. Also, results show that most detachment rates occur in convex profile curvatures, which are about 15 times more than in straight profiles. Copyright © 2015 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.     

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.     

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.     

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.     

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.     

A non‐equilibrium sediment transport model for rill erosion

Sediment transport in rill flows exhibits the characteristics of non‐equilibrium transport, and the sediment transport rate of rill flow gradually recovers along the flow direction by erosion. By employing the concept of partial equilibrium sediment transport from open channel hydraulics, a dynamic model of rill erosion on hillslopes was developed. In the model, a parameter, called the restoration coefficient of sediment transport capacity, was used to express the recovery process of sediment transport rate, which was analysed by dimensional analysis and determined from laboratory experimental data. The values of soil loss simulated by the model were in agreement with observed values. The model results showed that the length and gradient of the hillslope and rainfall intensity had different influences on rill erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Experimental investigation of the morphodynamic response of riffles and pools to unsteady flow and increased sediment supply

Gravel-bed rivers characteristically exhibit shallow riffles in wide sections and deeper pools where the channel becomes constricted and narrow. While rivers can adjust to changing flow and sediment supply through some combination of adjustments of channel slope, bed-surface sorting, and channel shape, the degree to which riffle-pools may adopt these changes in response to changing flows and sediment supplies remains unclear. This article presents results from a flume experiment investigating how constant- and variable-width channels adjust their morphology in response to changing flow and increased sediment supply. Two flume geometries were used: (1) constant-width and (2) variable-width, characterized by a sinusoidal pattern with a mean width equal to that of the first channel. The variable-width channel developed bed undulations in phase with the width, representing riffle-pools. The experiment consisted of three phases for each flume geometry: (1) steady flow, constant sediment supply; (2) unsteady flow, constant sediment supply; and (3) unsteady flow, doubled sediment supply. Unsteady flow was implemented in the form of repeated symmetrical stepped hydrographs, with a mean discharge equal to that in the steady flow phase. In all phases the bed and sediment supply were composed of a sand/gravel mixture ranging from 1 to 8 mm. In both the straight and variable-width channels, transitioning from steady flow to repeated hydrographs did not result in significant changes in bed morphology. The two channel geometries had different responses to increased sediment supply: the slope of the straight channel increased nearly 40%, while the variable-width channel reduced the relief between bars and pools and decreased the variability in cross-sectional elevation with a slight slope increase. Bar-pool relief varied with repeat discharge hydrographs. Pool elevation changed twice the distance of bar elevations, emphasizing the relevance of pool scour for riffle-pool self-maintenance in channels with width variations.     

Modelling the process of soil detachment by rill flow on steep loessial hillslopes

Soil detachment by rill flow is a key process of rill erosion, modelling this process can help in understanding rill erosion mechanisms. However, many soil detachment models are established on conceptual assumptions rather than experimental data. The objectives of this study were to establish a model of soil detachment by rill flow based on flume experimental data and to quantitatively verify the model. We simulated the process of soil detachment by rill flow in flume experiments with a soil-feeding hopper using loessial soil on steep slopes. Seven flow discharges, six slopes and five sediment loads were combined. Soil detachment capacity, sediment transport capacity, and soil detachment rate by rill flow under different sediment loads were measured. The process of soil detachment by rill flow can be modelled by a dual power function based on soil detachment capacity and transport capacity deficit as variables. The established model exhibited high credibility (NSE=0.97; R²=0.97). The contributions of soil detachment capacity and transport capacity deficit to soil detachment rate by rill flow reached 60% and 36%, respectively. Soil detachment capacity exerted more influence on soil detachment rate than did transport capacity deficit. The performance of the WEPP rill erosion equation is also favourable (NSE=0.95; R²=0.97). The two power exponents in the model we established strengthen the role of soil detachment capacity in soil detachment rate and weaken that for transport capacity deficit. Soil detachment capacity and transport capacity deficit played important roles in the determination of soil detachment rate by rill flow. The results can be applied to implement the numerical modeling and prediction of rill erosion processes on steep loessial hillslopes. © 2019 John Wiley & Sons, Ltd.     

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.     

Using a thermal tracer to estimate overland and rill flow velocities

Flow velocity is a basic hydraulic property of surface flows and its precise calculation is necessary for process based hydrological models, such as soil erosion and rill development models, as well as for modelling sediment and solute transport by runoff. This study presents a technique based on infrared thermography to visualize very shallow flows and allow a quantitative measurement of overland flow and rill flow velocities. Laboratory experiments were conducted to compare the traditional dye tracer technique with this new thermal tracer technique by injecting a combined tracer (heated dye) into shallow flowing surface water. The leading edge tracer velocities estimated by means of infrared video and by the usual real imaging video were compared. The results show that thermal tracers can be used to estimate both overland and rill flow velocities, since measurements are similar to those resulting from using dye tracers. The main advantage of using thermography was the higher visibility of the leading edge of the injected tracer compared with the real image videos. Copyright © 2013 John Wiley & Sons, Ltd.     

A MATHEMATICAL QUASI-2D OVERLAND FLOW MODEL

1 INTRODUCTIONThe prediction of future impacts on terrestrial ecosystems by atmospheric, climatic and land-usechanges is the aim of watershed management. Meeting these requirements scientists, managers and policymakers try to achieve the sustainable management of the vitally important resources of watersheds due toan integrated ecosystem approach at the catchment scale. As composite landscapes often have a highdegree of contingency between its elements, the transport over these landscape s…     

MODELING SEDIMENT PROCESSES ON SMALL WATERSHEDS:A CONCEPTUAL FRAMEWORK Ⅰ. BROAD SHALLOW FLOW PROCESSE

A conceptual modeling framework for developing process-based mathematical models of sediment generation, transport, and deposition on broad shallow flow areas is presented. The governing equations relevant to process-based modeling of broad shallow flow sediment processes on small watersheds are presented and the effects of space and time averaging on the predictive equations are described. Starting from the most general one-dimensional, unsteady model of sediment processes, simpler model structures are obtained and the successive simplifications made on the governing equations in order to obtain simpler and less detailed formulations are described. Specific model formulations are given for illustrative purposes and applications of these models to erosion and sediment yield prediction from broad shallow flow areas are shown using sediment data from rainfall simulator plots. In spite of some progress made in the development of process-based erosion and sediment yield models from broad shallow flow areas, further developmental modeling efforts must be based on a clear separation between hydrologic and hydraulic processes, and the soil properties which are significant for each.