Testing slope effect on flow resistance equation for mobile bed rills

In this paper, a recently theoretically deduced rill flow resistance equation, based on a power‐velocity profile, is tested experimentally on plots of varying slopes in which mobile bed rills are incised. Initially, measurements of flow velocity, water depth, cross‐sectional area, wetted perimeter and bed slope conducted in 106 reaches of rills incised on an experimental plot having a slope of 14% were used to calibrate the flow resistance equation. Then, the relationship between the velocity profile parameter Γ, the channel slope, and the flow Froude number, which was calibrated using the 106 rill reach data, was tested using measurements carried out in plots having slopes of 22% and 9%. The measurements carried out in the latter slope conditions confirmed that (a) the Darcy–Weisbach friction factor can be accurately estimated using the proposed theoretical approach, and (b) the data were supportive of the slope independence hypothesis of rill velocity stated by Govers.     

Estimating flow resistance in steep slope rills

Recent research recognized that the slope of 18% can be used to distinguish between the ‘gentle slope’ case and that of ‘steep slope’ for the detected differences in hydraulic variables (flow depth, velocity, Reynolds number, Froude number) and those representatives of sediment transport (flow transport capacity, actual sediment load). In this paper, using previous measurements carried out in mobile bed rills and flume experiments characterized by steep slopes (i.e., slope greater than or equal to 18%), a theoretical rill flow resistance equation to estimate the Darcy-Weisbach friction factor is tested. The main aim is to deduce a relationship between the velocity profile parameter Γ, the channel slope, the Reynolds number, the Froude number and the textural classes using a data base characterized by a wide range of hydraulic conditions, plot or flume slope (18%–84%) and textural classes (clay ranging from 3% to 71%). The obtained relationship is also tested using 47 experimental runs carried out in the present investigation with mobile bed rills incised in a 18%—sloping plot with a clay loam soil and literature data. The analysis demonstrated that: (1) the soil texture affects the estimate of the Γ parameter and the theoretical flow resistance law (Equation 25), (2) the proposed Equation (25) fits well the independent measurements of the testing data base, (3) the estimate of the Darcy-Weisbach friction factor is affected by the soil particle detachability and transportability and (4) the Darcy-Weisbach friction factor is linearly related to the rill slope.     

Flow resistance equation for rills

In this paper, a new flow resistance equation for rill flow was deduced applying dimensional analysis and self‐similarity theory. At first, the incomplete self‐similarity hypothesis was used for establishing the flow velocity distribution whose integration gives the theoretical expression of the Darcy–Weisbach friction factor. Then the deduced theoretical resistance equation was tested by some measurements of flow velocity, water depth, cross section area, wetted perimeter, and bed slope carried out in 106 reaches of some rills shaped on an experimental plot. A relationship between the velocity profile, the channel slope, and the flow Froude number was also established. The analysis showed that the Darcy–Weisbach friction factor can be accurately estimated by the proposed theoretical approach based on a power–velocity profile.     

Comparing flow resistance law for fixed and mobile bed rills

Rills caused by run‐off concentration on erodible hillslopes have very irregular profiles and cross‐section shapes. Rill erosion directly depends on the hydraulics of flow in the rills, which may differ greatly from hydraulics of flow in larger and regular channels. In this paper, a recently theoretically deduced rill flow resistance equation, based on a power–velocity profile, was tested experimentally on plots of varying slopes (ranging from 9% to 26%) in which mobile and fixed bed rills were incised. Initially, measurements of flow velocity, water depth, cross‐section area, wetted perimeter, and bed slope, carried out in 320 reaches of mobile bed rills and in 165 reaches of fixed rills, were used for calibrating the theoretical flow resistance equation. Then the relationship between the velocity profile parameter Γ, the channel slope, and the flow Froude number was separately calibrated for the mobile bed rills and for the fixed ones. The measurements carried out in both conditions (fixed and mobile bed rills) confirmed that the Darcy–Weisbach friction factor can be accurately estimated using the proposed theoretical approach. For mobile bed rills, the data were supportive of the slope independence hypothesis of velocity, due to the feedback mechanism, stated by Govers. The feedback mechanism was able to produce quasicritical flow conditions. For fixed bed rills, obtained by fixing the rill channel, by a glue, at the end of the experimental run with a mobile bed rill, the slope independence of the flow velocity measurements was also detected. Therefore, an experimental run carried out by a rill bed fixed after modelling flow action is useful to detect the feedback mechanism. Finally, the analysis showed that, for the investigated conditions, the effect of sediment transport on the flow resistance law can be considered negligible respect to the grain roughness effect.     

Testing a new rill flow resistance approach using the Water Erosion Prediction Project experimental database

In this paper, a recently theoretically deduced rill flow resistance equation, based on a power‐velocity profile, was tested using the Water Erosion Prediction Project database. This database includes measurements of flow velocity, water depth, cross section area, wetted perimeter, and bed slope that were made in rills shaped on experimental sites distributed across the continental United States. In particular, three different experimental conditions (only rainfall, only flow, and rain with flow) were examined, and for each condition, the theoretically based relationship for estimating the Γ function of the power velocity profile was calibrated. The results established that (a) the Darcy‐Weisbach friction factor can be accurately estimated using the proposed theoretical approach, and (b) the flow resistance increases with the effect of rainfall impact.     

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.     

Testing a theoretical resistance law for overland flow on a stony hillslope

Overland flow, sediments, and nutrients transported in runoff are important processes involved in soil erosion and water pollution. Modelling transport of sediments and chemicals requires accurate estimates of hydraulic resistance, which is one of the key variables characterizing runoff water depth and velocity. In this paper, a new theoretical power–velocity profile, originally deduced neglecting the impact effect of rainfall, was initially modified for taking into account the effect of rainfall intensity. Then a theoretical flow resistance law was obtained by integration of the new flow velocity distribution. This flow resistance law was tested using field measurements by Nearing for the condition of overland flow under simulated rainfall. Measurements of the Darcy–Weisbach friction factor, corresponding to flow Reynolds number ranging from 48 to 194, were obtained for simulated rainfall with two different rainfall intensity values (59 and 178 mm hr⁻¹). The database, including measurements of flow velocity, water depth, cross-sectional area, wetted perimeter, and bed slope, allowed for calibration of the relationship between the velocity profile parameter Γ, the slope steepness s, and the flow Froude number F, taking also into account the influence of rainfall intensity i. Results yielded the following conclusions: (a) The proposed theoretical flow resistance equation accurately estimated the Darcy–Weisbach friction factor for overland flow under simulated rainfall, (b) the flow resistance increased with rainfall intensity for laminar overland flow, and (c) the mean flow velocity was quasi-independent of the slope gradient.     

A full-scale study of Darcy-Weisbach friction factor for channels vegetated by riparian species

In this article, an open channel flow resistance equation, deduced applying dimensional analysis and incomplete self-similarity condition for the flow velocity distribution, was tested using measurements carried out in a full-scale channel equipped with three types of riparian plants (Salix alba L., Salix caprea L. and Alnus glutinosa L.). In the experimental channel, having banks lined with boulders, the vegetation branches were anchored in a concrete bottom. For each species, the measurements were carried out with plants having different amounts of leaves, different plant density and plant area index. The relationship between the scale factor Γ of the velocity profile and the Froude number was separately calibrated by measurements carried out without and with vegetation. The component of Darcy-Weisbach friction factor corresponding to the riparian vegetation f_v was calculated as the difference between the measured friction factor value (channel grain roughness + vegetation) and that calculated for the channel without vegetation in the same hydraulic conditions. Using these f_v values, the relationship between the scale factor Γ and the Froude number was calibrated. In this last relationship, a scaling coefficient a varying with the investigated vegetation type was introduced. This coefficient, as expected, gives the highest friction factor values for vegetation having branches with leaves. The theoretical flow resistance law, coupled with the relationship for estimating the Γ function having a scaling coefficient different for each investigated vegetation type, allowed an accurate estimate of the Darcy-Weisbach friction factor (errors less than or equal to 20% for 82.6% of the investigated cases). Finally, for the investigated vegetation species that are characterized by a condition with few leaves or leafless, the scaling coefficient a resulted strongly related to the bending stiffness. This analysis demonstrated that the highest Darcy-Weisbach friction factors correspond to vegetation species characterized by the highest values of bending stiffness. The friction factor values calculated for this last condition are characterized by errors that were less than or equal to ±20% for 90.6% of cases.     

The influence of roughness geometry and Shields parameter on flow resistance in gravel-bed channels

The spatial variability of bed particles of a gravel-bed channel is analysed and treated experimentally in order to simulate the effects of the arrangement of coarse bed elements on the flow resistance law. For the studied bed patterns, characterized by the concentration Γ of coarser elements arranged on the bed layer, a particle arrangement parameter α is proposed. The α parameter is useful for estimating the intercept b₀ of the semi-logarithmic flow resistance law deduced by flume measurements carried out for the hydraulic condition of large-scale and transition roughness. The differences between the experimental friction factor parameter values and the ones calculated by the proposed semi-logarithmic relationship are explained by the ratio between the Shields parameter and its critical value. The analysis shows that estimates of the friction factor parameter are not improved by introducing the Froude number into the flow resistance law. © 1997 by John Wiley & Sons, Ltd.     

Evaluating overland flow sediment transport capacity

An equation for evaluating the sediment transport capacity of overland flow is a necessary part of a physically based soil erosion model describing sediment detachment and transport as distributed processes. At first, for the hydraulic conditions of small-scale and large-scale roughness, the sediment transport capacity relationship used in the WEPP model is calibrated by Yalin and Govers' equation. The analysis shows that the transport coefficient K_t depends on the Shields parameter, Y, according to a semi-logarithmic (Yalin) or a linear (Govers) equation. The reliability of the semi-logarithmic equation is verified by Smart's, and Aziz and Scott's experimental data. Then the Low's formula, whose applicability is also proved by Smart's, and Aziz and Scott's data, is transformed as a stream power equation in which a stream power coefficient, K_SP, depending on Shields parameter, slope, sediment and water-specific weight, appears. A relationship between transport capacity and effective stream power is also proposed. Finally, the influence of rainfall on sediment transport capacity and the prediction of critical shear stress corresponding to overland flow are examined. © 1998 John Wiley & Sons, Ltd.     

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.     

Roughness effect on the correction factor of surface velocity for rill flows

Flow velocity is one of the most important hydrodynamic variables for both channelized (rill and gullies) and interrill erosive phenomena. The dye tracer technique to measure surface flow velocity V_s is based on the measurement of the travel time of a tracer needed to cover a known distance. The measured V_s must be corrected to obtain the mean flow velocity V using a factor α_v = V/V_s which is generally empirically deduced. The V_s measurement can be influenced by the method applied to time the travel of the dye-tracer and α_v can vary in different flow conditions. Experiments were performed by a fixed bed small flume simulating a rill channel for two roughness conditions (sieved soil, gravel). The comparison between a chronometer-based (CB) and video-based (VB) technique to measure V_s was carried out. For each slope-discharge combination, 20 measurements of V_s, characterized by a sample mean V_m, were carried out. For both techniques, the frequency distributions of V_s/V_m resulted independent of slope and discharge. For a given technique, all measurements resulted normally distributed, with a mean equal to one, and featured by a low variability. Therefore, V_m was considered representative of surface flow velocity. Regardless of roughness, the V_m values obtained by the two techniques were very close and characterized by a good measurement precision. The developed analysis on α_v highlighted that it is not correlated with Reynolds number for turbulent flow regime. Moreover, α_v is correlated neither with the Froude number nor with channel slope. However, the analysis of the empirical frequency distributions of the correction factor demonstrated a slope effect. For each technique (CB, VB)-roughness (soil, gravel) combination, a constant correction factor was statistically representative even if resulted in less accurate V estimations compared to those yielded by the slope-specific correction factor.     

Longitudinal velocity patterns and bed morphology interaction in a rill

Present‐day understanding of rill dynamics is hampered by a lack of detailed data on velocity distributions in rills. The latter are dif?cult to collect with traditional techniques due to the very low water depths and the relatively high ?ow velocities in rills. The objectives of this paper were to investigate the feasibility of miniaturized acoustic Doppler velocimeter (mADV) measurements in rill ?ow and to explore longitudinal variations in ?ow velocities and their relationship with rill bed morphology. Detailed data on longitudinal ?ow velocity were required to achieve these objectives. A 1·8 m long rill was formed freely in a ?ume at 5° slope and 0·001 m³ s^?1 discharge. Rill topography was characterized by an alternation of steps and pools. The ?ume surface was then ?xed to preserve rill roughness. A topographical scanning of the entire ?ume surface was made. Velocity was measured with a mADV along the rill, and at different depths. Flow depth in a longitudinal direction was also measured using an elevation gauge. A strong relationship exists between rill topography and ?ow hydraulics. Over steps, ?ow was unidirectional and rapidly accelerating until a threshold Froude number (Fn) value between 1·3 and 1·7 was reached and a hydraulic jump occurred leading to the formation of a pool. In the pool, the ?ow pattern was multidirectional and complex. The ?ow was subcritical when leaving the pool and accelerated over the next step until the threshold Froude number value was again reached. Energy loss in the rill was concentrated in the pools, mainly due to the action of a hydraulic jump. This mechanism of energy dissipation appeared to be an essential factor in rill formation and bedform evolution. Copyright © 2004 John Wiley & Sons, Ltd.     

Variability of shallow overland flow velocity and soil aggregate transport observed with digital videography

Field experiments at Tiramoana station 30 km north of Christchurch, New Zealand using an erosion plot 16·5 m long, 0·6 m wide, and with a slope of 14–14·5° on rendzina soil aimed to measure the variability of flow velocity and of soil aggregates transport rate in shallow overland flow. Discharge/cross‐section area ratio was used to estimate mean velocity, and high‐speed digital video camera and image analysis provided information about flow and sediment transport variability. Six flow runs with 0·5–3·0 L s^?1 discharges were supercritical with Froude numbers close to or more than 1. Mean flow velocity followed Poiseuille law, float numbers were more than 1·5 and hydraulic resistance was an inverse proportional function of the Reynolds number, which is typical for laminar flows. Hence actual velocity varied through time significantly and the power spectrum was of ‘red‐noise’, which is typical for turbulent flow. Sediment transport rates had even higher variability, and soil aggregates transport was a compound Poisson process. Copyright © 2008 John Wiley & Sons, Ltd.     

Roughness coefficients and velocity estimation in well-inundated sheet and rilled overland flow without strongly eroding bed forms

Even with the flow of water over a soil surface in which roughness elements are well inundated, and in less erosive situations where erosional bed forms are not pronounced, the magnitude of resistance coefficients in equations such as those of Darcy–Weisbach, Chezy or Manning vary with flow velocity (at least). Using both original laboratory and field data, and data from the literature, the paper examines this question of the apparent variation of resistance coefficients in relation to flow velocity, even in the absence of interaction between hydraulics and resulting erosional bed forms. Resistance equations are first assessed as to their ability to describe overland flow velocity when tested against these data sources. The result is that Manning's equation received stronger support than the Darcy–Weisbach or Chezy equations, though all equations were useful. The second question addressed is how best to estimate velocity of overland flow from measurements of slope and unit discharge, recognizing that the apparent flow velocity variation in resistance coefficients is probably a result of shortcomings in all of the listed resistance equations. A new methodology is illustrated which gives good agreement between estimated and measured flow velocity for both well-inundated sheet and rill flow. Comments are given on the predictive use of this methodology. Copyright © 1999 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.     

Effect of particle density and inflow concentration of suspended sediment on bedload transport in rill flow

Laboratory flume experiments were carried out to evaluate the effect of particle density on bedload transport of sand‐sized particles and the effect of a suspended load of clay particles (kaolinite) on bedload transport of sand‐sized particles in rill flow conditions. Three materials in the range 400–600 µm were selected to simulate bedload transport of primary particles and aggregates: sand (2650 kg/m³), crushed brick (2450 kg/m³) and anthracite (1300–1700 kg/m³). In the two first experiments, two different methods were applied to determine bedload transport capacity of coarse particles for various conditions of flow discharge (from 2 to 15 L/min) and slope (2.2, 3 and 4%). In the third experiment, clear water was replaced with kaolinite–water mixture and bedload transport capacity of crushed brick particles was determined for a 4% slope and different concentrations of kaolinite (0, 7, 41 and 84 g/L). The results showed that bedload transport increased significantly with the decrease in particle density but the effect of particle density on transport rates was much less important than flow discharge. Velocity measurements of clear flow, flow mixed with coarse particles and coarse particles confirmed the existence of a differentiation between suspended load and bedload. In these experimental conditions, suspended load of kaolinite did not affect bedload rates of crushed brick particles. Three transport capacity formulae were tested against observed bedload rates. A calibration of the Foster formula revealed that the shear stress exponent should be greater than 1.5. The Low and the Govers unit stream power (USP) equations were then evaluated. The Low equation was preferred for the prediction of bedload rates of primary particles but it was not recommended in the case of aggregates of low density because of the limited experimental conditions applied to derive this equation. Copyright © 2008 John Wiley & Sons, Ltd.     

Comment on “Overland runoff erosion dynamics on steep slopes with forages under field simulated rainfall and inflow by C. Li and C. Pan”

Managing sloping landscapes to control soil erosion processes due to rainfall and runoff is a relevant problem, especially when the vegetation is absent or sparse. The aim of this paper was to investigate the applicability of a theoretically resistance law for overland flow under simulated rainfall, based on a power-velocity profile, using field measurements carried out by Li and Pan for three plots with planted forage species (Astragalus adsurgens, Medicago sativa and Cosmos bipinnatus).The relationship between the velocity profile parameter Γ, the flow Froude number and the rain Reynolds number was calibrated using the data by Li and Pan. The obtained overland flow resistance law was also verified by independent field measurements carried out by Li and Pan in the same plots with the same forage species subjected to three different treatments (intact grass control IG, no litter or leaves NLL (only the grass stems and roots were reserved) and only roots remaining OR). The theoretical approach and the measurements carried out in the investigated conditions allowed to state that a) the Darcy-Weisbach friction factor can be accurately estimated using the proposed theoretical approach, b) the Darcy-Weisbach friction factor varies with rainfall intensity and c) for the investigated forage species, the vegetation treatment (IG, NLL, OR) does not significantly affect the flow resistance in laminar regime.     

Hydraulic properties affected by litter and stem cover under overland flow

Vegetation stems and litter cover have different effects on sediment transport capacity under the same experimental conditions, which in essence, may be due to differences in their hydraulic properties, but the availability of comparative studies is limited. This study aimed to compare the hydraulic properties affected by litter and stem cover, compare differences in the drag forces exerted by litter and stems on overland flow, and develop new Manning's n and flow velocity equations for litter cover. Two series of flume experiments were conducted with the same slope gradients (8.8%, 17.6%, 26.8%) and flow discharge rates (0.5, 1.0 × 10⁻³ m³ s⁻¹). Artificial Gramineae stems with a 0%–30% cover level and Pinus tabulaeformis litter with a 0%–70% cover level were used in series 1 and series 2, respectively. The flow velocity and depth were measured. The results showed that the Froude number and flow velocity affected by stem cover were much lower than those affected by litter cover, while the opposite trend was observed in the relative magnitude of the Reynolds number, flow depth and shear stress. The form resistance caused by stems was 22–57 times greater than that caused by litter for the same cover level, which suggests that stem cover contributes more than litter cover to increasing the flow resistance and reducing the flow's ability for sediment detachment and transport. Two new equations for calculating Manning's n and flow velocity under the influence of litter cover were developed, with R² and NSE values of 0.96. The results of this study contribute to revealing the mechanisms of the differences of the effects of stem and litter cover on soil erosion.