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.     

Response of soil detachment rate to the hydraulic parameters of concentrated flow on steep loessial slopes on the Loess Plateau of China

Using hydraulic parameters is essential for describing soil detachment and developing physically based erosion prediction models. Many hydraulic parameters have been used, but the one that performs the best for describing soil detachment on steep slopes when the lateral expansion (widening) of rills is not limited has not been identified. An indoor concentrated flow scouring experiment was performed on steep loessial slopes to investigate soil detachment rates for different flow rates and slope gradients. The experiments were conducted on a slope‐adjustable plot (5 m length, 1 m width, 0.5 m depth). Sixteen combinations of 4 flow rates (10, 15, 20, and 25 L/min) and 4 slope gradients (17.6%, 26.8%, 36.4%, and 46.6%) were investigated. The individual and combined effects of slope gradient and flow hydraulic parameters on soil detachment rate were analysed. The results indicated that soil detachment rate increased with flow rate and slope gradient. Soil detachment rate varied linearly and exponentially with flow rate and slope gradient, respectively. Multivariate, nonlinear regression analysis indicated that flow depth exerted the greatest influence on the soil detachment rate, followed by unit discharge per unit width, slope gradient, and flow rate in this study. Shear stress and stream power could efficiently describe the soil detachment rate using a power equation. However, the unit stream power and unit energy of the water‐carrying section changed linearly with soil detachment rate. Stream power was an optimal hydraulic parameter for describing soil detachment. These findings improve our understanding of concentrated flow erosion on steep loessial slopes.     

Water flow velocity over frozen and nonfrozen black soil slopes

Ice‐ and snow‐melted water flow over partially thawed frozen soil of cultivated slopes causes serious soil erosion, which results in soil degradation and affects productivity in Northeast China. Water flow velocity over frozen and nonfrozen soil shows importance in understanding meltwater erosion. In this work, a series of laboratory experiments were conducted to measure water flow velocity over frozen and nonfrozen soil slopes. Experiments were performed using the electrolyte trace method under the pulse boundary model, under conditions of 4 slope gradients (5°, 10°, 15°, and 20°), 3 flow rates (1, 2, and 4 L/min), and 7 sensors positioned at 0.1, 1.0, 2.0, 3.0, 4.0, 5.0, and 6.0 m away from the electrolyte injection point. Results showed that velocities over frozen soil slopes increased with flow rate and slope gradient. Flow velocities over nonfrozen soil slopes increased with flow rate and slope gradients from 5° to 15° and stabilized at 15°. Flow velocities over frozen soil slopes were 30%, 54%, 71%, and 91% higher than those over nonfrozen ones at slope gradients of 5°, 10°, 15°, and 20°. Flow velocities over frozen soil slopes under different flow rates of 1, 2, and 4 L/min were approximately 52%, 59%, and 79% higher than those over nonfrozen soil, respectively. This study can help in assessing the erosion of partially thawed frozen soil by meltwater flow.     

Rill erosion along the thalweg of a hillslope hollow: A case study from the cotswold hills,Central England

Rill erosion is an important erosional form on agricultural soils in England, causing large losses of soil, particularly on cultivated slopes. This paper describes a rill system that developed in a small agricultural catchment in north Oxfordshire during the winter of 1992–93. The rill system comprised two components: a system of ‘feeder rills’ along the valley-side slopes, which were the result of flow concentration and erosion along wheelings, and a thalweg rill, which formed along a dry valley bottom as a result of surface runoff concentration from the feeder rills. Total volumetric soil loss from the rill system was 32·28 m³, equivalent to 3·01 m³, ha^?1 for the rill catchment area, or 3·91 t ha^?1. Mean discharge for the thalweg rill and feeder rills, calculated during a storm event, was 31·101s^?1 and 1·171s^?1, respectively. All flows were fully turbulent and supercritical. We emphasize the need for a spatially distributed approach to the study of runoff and erosion at the catchment scale.     

Hydraulics of overland flow influenced by litter incorporation under extreme rainfall

Water erosion on hillslopes is a worldwide environmental problem, which is a rainfall‐induced process, especially extreme rainfall. The great intensity of extreme rainfall strongly enhances the power of overland flow to detach soil and transport sediment. Plant litter is one of the most important constituents of ecosystems that often covers the soil surface and can be incorporated into topsoil. However, little attention has been paid to its effect on flow hydraulics owing to the veiled nature. This study aimed to examine the effects of incorporated litter on the hydraulic properties under extreme rainfall condition. To reach this goal, six litter rates of 0, 0.05, 0.10, 0.20, 0.35, and 0.50 kg m^?2 and four litter types collected from deciduous trees, coniferous trees, shrubs, and herbs were incorporated into topsoil. Then, simulated rainfall experiments were performed on five slope gradients (5°, 10°, 15°, 20°, and 25°) with an extreme rainfall intensity of 80 mm h^?1. The results showed that Froude number and flow velocity of the overland flow decreased, whereas flow resistance increased exponentially with litter incorporation rate. Litter type had an influence on flow hydraulics, which can mainly be attributed to the variations in surface coverage of the exposed litter and the litter morphology. Flow velocity and Darcy–Weisbach coefficient increased markedly with slope gradient. However, the variation of slope gradient did not modify the relationships between flow hydraulics and incorporated litter rate. The random roughness, resulting from heterogeneous erosion due to the uneven protection of surface exposed litter, increased linearly with litter incorporated rate. As rainfall proceeded, flow hydraulics varied with incorporated litter rate and slope gradient complicatedly due to the increases in flow rate and coverage of the exposed litter and the modification of soil surface roughness.     

Effects of sediment size on transport capacity of overland flow on steep slopes

Abstract

Sediment transport capacity is a key concept in determining rates of detachment and deposition in process-based erosion models, yet limited studies have been conducted on steep slopes. We investigated the effects of sediment size on transport capacity of overland flow in a flume. Unit flow discharge ranged from 0.66 to 5.26?×?10^-3 m² s^-1, and slope gradient varied from 8.7 to 42.3%. Five sediment size classes (median diameter, d ₅₀, of 0.10, 0.22, 0.41, 0.69 and 1.16 mm) were used. Sediment size was inversely related to transport capacity. The ratios of average transport capacity of the finest class to those of the 0.22, 0.41, 0.69 and 1.16 mm classes were 1.09, 1.30, 1.55 and 1.92, respectively. Sediment transport capacity increased as a power function of flow discharge and slope gradient (R²?=?0.98), shear stress (R²?=?0.95), stream power (R²?=?0.94), or unit stream power (R²?=?0.76). Transport capacity generally decreased as a power function of sediment size (exponent?=??0.35). Shear stress and stream power predicted transport capacity better than unit stream power on steep slopes when transport capacity was <7 kg m^-1 s^-1. Sediment transport capacity increased linearly with mean flow velocity. Critical or threshold velocity increased as a power function of sediment size (R²?=?0.93). Further studies with fine soil particles are needed to quantify the effects of sediment size on transport capacity of overland flow on steep slopes.

Citation Zhang, G.-H., Wang, L.-L., Tang, K.-M., Luo, R.-T. & Zhang, X.C. (2011) Effects of sediment size on transport capacity of overland flow on steep slopes. Hydrol. Sci. J. 56(7), 1289–1299.     

Velocities of shallow overland flow on semiarid hillslopes

A series of 188 rainfall plot simulations was conducted on grass, shrub, oak savanna, and juniper sites in Arizona and Nevada. A total of 897 flow velocity measurements were obtained on 3.6% to 39.6% slopes with values ranging from 0.007 m s^‐1 to 0.115 m s^‐1. The experimental data showed that shallow flow velocity on rangelands was related to discharge and ground litter cover and was largely independent of slope gradient or soil characteristics. A power model was proposed to express this relationship. These findings support the slope–velocity equilibrium hypothesis. Namely, eroding soil surfaces evolve such that steeper areas develop greater hydraulic roughness. As a result overland flow velocity becomes independent of the slope gradient over time. Our findings have implications for soil erosion modeling suggesting that hydraulic friction is a dynamic, slope and discharge dependent property. Copyright © 2018 John Wiley & Sons, Ltd.     

Dynamics of fluvial and slope processes in the changing geomorphic environment of singapore

The nature and rates of fluvial and slope processes change over time and space as urbanized areas replace forested land in Singapore. Storm-based and time-based data, from undisturbed rainforests, heavily disturbed construction sites, urban grass-covered slopes and an experimental plot, are collected to observe the impact of rainwater on the soil moisture conditions, surface microtopography, runoff generation, sediment movement, and ground lowering in the three different categories of land use. The undisturbed forested environment is characterized by high throughfall (58% of total rainfall) and frequent negative soil moisture suctions. The slow and unconcentrated overland flow during heavy storms is restricted by the forest floor microtopography. No rills develop. Ground lowering is recorded as 3·2–3·4 mm a^?1. But sediment movement is episodic and suspended sediment concentrations in overland flow are 172–222 mg l^?1. During urban construction, gully development is rapid on the bare slopes, runoff generation, voluminous, and sediment-laden discharges (5200–75498 mg l^?1) lead to sediment plumes at channel mouths. Ground lowering rates are measured at 132·4 mm a^?1. Once grass-covered, runoff carries less suspended sediment (800 mg l^?1) and ground lowering rates are reduced, but depend on the condition of the cover, ranging from 0·2 to 8·2 mm a^?1. As urban development continues, environments are altered both in time as well as spatially.     

Establishment of a sediment transport capacity equation on loessal slope via experimental investigation

Areliable sediment transport capacity function provides response against challenges of soil erosion prediction on the Loess Plateau of China. The popular sediment transport capacity functions are questionable on loess slopes because the experimental conditions from which they were derived, like bed materials, gentle slopes,and surface roughness, are different from soil erosion processes on the loess slopes. Due to the foregoing uncertainty, a suitable sediment transport capacity function was dev...     

Effects of rainfall intensity and slope gradient on the application of vetiver grass mulch in soil and water conservation

The objective of this study is to investigate the effect of rainfall intensity and slope gradient on the performance ofvetiver grass mulch (VGM) in soil and water conservation.The study involved field ...     

Impacts of sediment load on Manning coefficient in supercritical shallow flow on steep slopes

The Manning equation is one of the most widely used formulae for calculating the velocity of shallow overland flow in hydrological and erosion models. Precise estimation of the Manning's friction coefficient (n) is critical to determining overland flow and soil erosion processes. Few studies have been conducted to quantify the effects of sediment load on Manning's n on steep slopes. This study was conducted to investigate the potential effects of sediment load on Manning's n in a flume with a fixed bed, under wide ranges of hydraulics and sediment loads. Slope gradient varied from 8·7 to 34·2%, unit flow rate from 0·66 to 5·26 × 10^?3 m² s^?1, and sediment load from 0 to 6·95 kg m^?1 s^?1. The Reynolds number ranged from 350 to 5899. Results showed that Manning's n varied in both sediment‐free and sediment‐laden flows ranging from 0·012 to 0·055. The apparent Manning's coefficients of sediment‐laden flow were much greater than those of sediment‐free flow. The mean Manning coefficient of sediment‐laden flow was 51·27% greater than the mean value of sediment‐free flow. For sediment‐laden flow, Manning's n could be estimated with a power function of unit flow discharge and sediment content. Further studies are needed to quantify the potential effects of sediment load on the Manning's n on erodible beds and in fields. Copyright © 2010 John Wiley & Sons, Ltd.     

Rill hydraulics on a semiarid hillslope,southern Arizona

Seventy field experiments were conducted in seven rills located on a semiarid rangeland hillslope underlain by gravelly soils at Walnut Gulch, Arizona. The rills, which are characterized by wide, shallow cross-sections and gravel-covered beds, have mean at-a-station hydraulic geometry exponents of b = 0·33, f = 0·34 and m = 0·33. Although the differences between these values and typical values of b = 0·30, f = 0·40 and m = 0·30 for cropland rills are not statistically significant, they are thought to be real, as cropland rills often have more rectangular cross-sections and steeper sides than the rangeland rills under study. For rills formed in silty loamy soils, Govers developed an empirical relation between mean flow velocity and discharge. Emphasizing the generality of this relation, he suggested that it may be used as a simple means of routing runoff through rills. He also noted that this relation appeared to be unaffected by either slope or soil materials. The present data represent rills underlain by coarser and somewhat more varied gravel-rich soils. These data do not conform to Govers' relation, and a multiple regression analysis reveals that slope and soil materials, either directly or indirectly through bed roughness, exert almost as much influence on flow velocity as does discharge. Three alternative methods are developed for predicting flow velocity in the rills under study. All three methods give good results with the largest root mean square deviation being 3·115 cm s⁻¹.     

Relationship between discharge,velocity and flow area for rills eroding loose,non-layered materials

A relationship between discharge, flow velocity and flow area in rills is established using data from four field and laboratory studies. The proposed relationship is shown to predict successfully flow velocities measured in six other studies. Although slopes range from 0.035 to 0.45 and soil materials range from stony sands over silt loams to vertisols, mean flow velocity can be well predicted from discharge alone. Thus, there is no important influence of slope and/or soil material characteristics on flow velocities in rills. The proposed relationship may be used to improve performance of deterministic flow routing models when applied to rilled catchments. Furthermore, it allows the calculation of unit stream power, which has been shown to be related to the transporting capacity of overland flow, in terms of slope and discharge.     

Hydrodynamic characteristics of rill flow on steep slopes

Rill erosion is a dominant sediment source on sloping lands. However, the amount of soil loss from rills on steep slopes is vastly more than that on gentle slopes because of differences in rill shape and hydraulic patterns. The aims of this paper are to determine the hydrodynamic characteristics of rills and the friction coefficients in steep slope conditions and to propose modifications of some hydraulic parameters used in soil loss prediction models. A series of inflow experiments was conducted on loess slopes. The results show that the geometric and hydraulic properties of rill on the steep loess slopes, which are characterized by the mean width of cross sections, mean velocity and mean depth of flow, are related to discharge and slope gradient in power functions. However, the related exponents to discharge are 0.26, 0.48 and 0.26, respectively, which are different from the exponents derived in previous studies, which were conducted on gentle slopes. The Manning roughness coefficient ranged from 0.035 to 0.071, with an average of 0.0536, and the Darcy–Weisbach friction coefficients varied from 0.4 to 1.9. The roughness coefficients are closely related to the Reynolds numbers and flow volumes; however, the correlations vary with slope gradient. The roughness coefficients are directly proportional to the Reynolds number and the flow volume on steep slopes, in contrast with the roughness coefficients found on gentle slopes, which decrease as the Reynolds number and flow volume increase. This difference is caused by the interactions among the hydraulics of the flow, the shape of the rills and the sediment concentrations on steep slopes. The results indicate that parameters used in models to predict rill erosion have to be modified according to slope gradient. Copyright © 2015 John Wiley & Sons, Ltd.     

Prediction of concentrated flow width in ephemeral gully channels

Empirical prediction equations of the form W = aQ^b have been reported for rills and rivers, but not for ephemeral gullies. In this study six experimental data sets are used to establish a relationship between channel width (W, m) and flow discharge (Q, m³ s^?1) for ephemeral gullies formed on cropland. The resulting regression equation (W = 2·51 Q^0·412; R² = 0·72; n = 67) predicts observed channel width reasonably well. Owing to logistic limitations related to the respective experimental set ups, only relatively small runoff discharges (i.e. Q < 0·02 m³s^?1) were covered. Using field data, where measured ephemeral gully channel width was attributed to a calculated peak runoff discharge on sealed cropland, the application field of the regression equation was extended towards larger discharges (i.e. 5 × 10^?4m³s^?1 < Q < 0·1 m³s^?1). Comparing W–Q relationships for concentrated flow channels revealed that the discharge exponent (b) varies from 0·3 for rills over 0·4 for gullies to 0·5 for rivers. This shift in b may be the result of: (i) differences in flow shear stress distribution over the wetted perimeter between rills, gullies and rivers, (ii) a decrease in probability of a channel formed in soil material with uniform erosion resistance from rills over gullies to rivers and (iii) a decrease in average surface slope from rills over gullies to rivers. The proposed W–Q equation for ephemeral gullies is valid for (sealed) cropland with no significant change in erosion resistance with depth. Two examples illustrate limitations of the W–Q approach. In a first example, vertical erosion is hindered by a frozen subsoil. The second example relates to a typical summer situation where the soil moisture profile of an agricultural field makes the top 0·02 m five times more erodible than the underlying soil material. For both cases observed W values are larger than those predicted by the established channel width equation for concentrated flow on cropland. For the frozen soils the equation W = 3·17 Q^0·368 (R² = 0·78; n = 617) was established, but for the summer soils no equation could be established. Copyright © 2002 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.     

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.     

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.     

Rill characteristics and sediment transport as a function of slope length during a storm event on loess soil

On the basis of detailed rill surveys carried out on bare plots of different lengths at slopes of 12 per cent, basic rill parameters were derived. Rill width and maximum depth increased with plot length, whereas rill amount and cross‐sectional area, expressed per unit length, remained similar. On smaller plots, all rills were connected in a continuous transport system reaching the plot outlet, whilst on larger plots (10 and 20 m long) part of the rills ended with a deposition areas inside the plots. Amounts of erosion, calculated from rill volume and soil bulk density, were compared with soil loss measured at the plot outlets. On plots 10 and 20 m long, erosion estimated from volume of all rills was larger than measured soil loss. The latter was larger than erosion estimated from volume of contributing rills. To identify contributing soil loss area on these plots, two methods were applied: (i) ratio of total soil loss to maximum soil loss per unit area, and (ii) partition of plot area according to the ratio of contributing to total rill volume. Both methods resulted in similar areas of 21·8–23·5 m² for the plot 10 m long and 31·2 m² for the plot 20 m long. Identification of contributing areas enabled rill (5·9 kg m^?2) and interrill (2·6 kg m^?2) erosion rate to be calculated, the latter being very close to the value predicted from the Universal Soil Loss Equation. Although rill and interrill rates seemed to be similar on all plots, their ratio increased slightly with plot length. Application of this ratio to compute slope length factor of the Revised Universal Soil Loss Equation resulted in similar values to those predicted with the model. The achieved balance of soil loss suggested that all the sediment measured at the plot outlet originated from contributing rills and associated contributing rill areas. The results confirmed the utility of different plot lengths as a research tool for analysing the dynamic response of soil to rainfall–runoff. Copyright © 2005 John Wiley & Sons, Ltd.     

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.