Effect of saltating sediment on flow resistance and bed roughness in overland flow

The acceleration of saltating grains by overland flow causes momentum to be transferred from the flow to the grains, thereby increasing flow resistance and bed roughness. To assess the impact of saltating sediment on overland flow hydraulics, velocity profiles in transitional and turbulent flows on a fixed sand-covered bed were measured using hot-film anemometry. Five discharges were studied. At each discharge, three flows were measured: one free of sediment, one with a relatively low sediment load, and one with a relatively high sediment load. In these flows from 83 to 90 per cent of the sediment was travelling by saltation. As a result, in the sediment-laden flows the near-bed velocities were smaller and the velocity profiles steeper than those in the equivalent sediment-free flows. Sediment loads ranged up to 87·0 per cent of transport capacity and accounted for as much as 20·8 per cent of flow resistance (measured by the friction factor) and 89·7 per cent of bed roughness (measured by the ratio of the roughness length to median grain diameter). It is concluded that saltating sediment has a considerable impact on overland flow hydraulics, at least on fixed granular beds. Saltation is likely to have a relatively smaller effect on overland flow on natural hillslopes and agricultural fields where form and wave resistance dominate. Still, saltation is generally of greater significance in overland flow than in river flow, and for this reason its effect on overland flow hydraulics is deserving of further study. © 1998 John Wiley & Sons, Ltd.     

Partitioning resistance to overland flow on rough mobile beds

For overland flows transporting predominantly bed load over rough mobile beds without rainfall, resistance to flow f may be divided into four components: surface resistance f_s, form resistance f_f, wave resistance f_w, and bed‐mobility resistance f_m. In this study it is assumed that f = f_s + f_f + f_w + f_m, and an equation is developed for each component. The equations for f_s and f_f are borrowed from the literature, while those for f_w and f_m are developed from two series of flume experiments in which the beds are covered with various concentrations of large‐scale roughness elements. The first series consists of 65 experiments on fixed beds, while the second series contains 194 experiments on mobile beds. All experiments were performed on the same slope (S = 0·114) and with the same size of sediment (D = 0·00074 m). The equations for f_w and f_m are derived by a combination of dimensional analysis and regression analysis. The analyses reveal that the major controls of f_w and f_m are the Froude number F and the concentration of the roughness elements C_r. When the equations for f_w and f_m are summed, the C_r terms cancel out, leaving f_w+m = 0·63F^?2. An equation is developed that predicts total f, and the contributions of f_s, f_f, f_w and f_m to f are computed from the series 1 and 2 experiments. An analysis of the first series reveals that in clear‐water flows over fixed beds, f_w accounts for 52 per cent of f. A similar analysis of the second series indicates that in sediment‐laden flows over mobile beds f_w comprises 37 per cent and f_m 32 per cent of f, so that together f_w and f_m account for almost 70 per cent of f. Finally, regression analyses indicate that where F > 0·5, f_w and f_m each vary with F^?2 and f_w/f_m = 1·18. The equation developed here for predicting total f applies only to the range of hydraulic, sediment, and bed roughness conditions represented by the experimental data. With additional data from a broader range of conditions the same methodology as employed here could be used to develop a more general equation. Copyright © 2006 John Wiley & Sons, Ltd.     

Terrestrial laser scanning soil surfaces: a field methodology to examine soil surface roughness and overland flow hydraulics

Conventional roughness–resistance relationships developed for pipe and open‐channel flows cannot accurately describe shallow overland flows over natural rough surfaces. This paper develops a new field methodology combining terrestrial laser scanning (TLS) and overland flow simulation to provide a high‐resolution dataset of surface roughness and overland flow hydraulics as simulated on natural bare soil surfaces. This method permits a close examination of the factors controlling flow velocity and a re‐evaluation of the relationship between surface roughness and flow resistance. The aggregate effect of flow dynamics, infiltration and depression storage on retarding the passage of water over a surface is important where runoff‐generating areas are distant from well‐defined channels. Experiments to separate these effects show that this ‘effective resistance’ is dominated by surface roughness. Eight measurements of surface roughness are found to be related to flow resistance: standard deviation of elevations, inundation ratio, pit density (measured both perpendicular and parallel to the flow direction), slope, median depth, skewness of the depth distribution and frontal area. Hillslope position is found to affect the significant roughness measures. In contrast, infiltration rate has little effect on the velocity of water fronts advancing over the soil surfaces examined here and the effect of depression storage is limited. Overland flow resistance is depth dependent where complex microtopographic structures are progressively inundated. Copyright © 2010 John Wiley & Sons, Ltd.     

Predicting sediment transport by interrill overland flow on rough surfaces

Modelling soil erosion requires an equation for predicting the sediment transport capacity by interrill overland flow on rough surfaces. The conventional practice of partitioning total shear stress into grain and form shear stress and predicting transport capacity using grain shear stress lacks rigour and is prone to underestimation. This study therefore explores the possibility that inasmuch as surface roughness affects flow hydraulic variables which, in turn, determine transport capacity, there may be one or more hydraulic variables which capture the effect of surface roughness on transport capacity suffciently well for good predictions of transport capacity to be achieved from data on these variables alone. To investigate this possibility, regression analyses were performed on data from 1506 flume experiments in which discharge, slope, water temperature, rainfall intensity, and roughness size, shape and concentration were varied. The analyses reveal that 89·8 per cent of the variance in transport capacity can be accounted for by excess flow power and flow depth. Including roughness size and concentration in the regression improves that explained variance by only 3·5 per cent. Evidently, flow depth, when used in combination with excess flow power, largely captures the effect of surface roughness on transport capacity. This finding promises to simplify greatly the task of developing a general sediment equation for interrill overland flow on rough surfaces. Copyright © 1998 John Wiley & Sons, Ltd.     

A simple model of hillslope response for overland flow generation

This paper deals with the derivation of the hydrological response of a hillslope on the assumption of quick runoff by surface runoff generation. By using the simple non‐linear storage based model, first proposed by Horton, an analytical solution of the overland flow equations over a plane hillslope was derived. This solution establishes a generalization for different flow regimes of Horton's original solution, which is valid for the transitional flow regime only. The solution proposed was compared successfully with that of Horton and, for the turbulent flow regime, to the one derived from kinematic wave theory. This solution can be applied easily to both stationary and non‐stationary rainfall excess events. An analytical solution for the instantaneous response function (IRF) was also derived. Finally, simple expressions to compute peak and time to peak of IRF are proposed. Copyright © 2001 John Wiley & Sons, Ltd.     

Rill flow resistance law under equilibrium bed‐load transport conditions

In this paper, a recently deduced flow resistance equation for open channel flow was tested under equilibrium bed‐load transport conditions in a rill. First, the flow resistance equation was deduced applying dimensional analysis and the incomplete self‐similarity condition for the flow velocity distribution. Then, the following steps were carried out for developing the analysis: (a) a relationship (Equation  13 ) between the Γ function of the velocity profile, the rill slope, and the Froude number was calibrated by the available measurements by Jiang et al.; (b) a relationship (Equation  17 ) between the Γ function, the rill slope, the Shields number, and the Froude number was calibrated by the same measurements; and (c) the Darcy–Weisbach friction factor values measured by Jiang et al. were compared with those calculated by the rill flow resistance equation with Γ estimated by Equations  13 and 17 . This last comparison demonstrated that the rill flow resistance equation, in which slope and Shields number, representative of sediment transport effects, are introduced, is characterized by the lowest values of the estimate errors.     

The effect of bed mobility on resistance to overland flow

When sediment grains are transported as bed load in overland flow, there is a net transfer of momentum from the flow to the grains. When these grains collide with other grains, whether on the bed or in the flow, streamwise flow velocity decreases and resistance to flow increases. Resistance to flow generated in this manner is termed bed‐load transport resistance. Resistance to flow f over a plane bed may be partitioned into grain resistance f_g and bed‐load transport resistance f_bt. We use the symbols f_btf and f_btm to denote f_bt for flows over fixed beds and over mobile beds, respectively, and we compute the effect of bed mobility on flow resistance f_mob by subtracting f_btf from f_btm. The data for this study come from 54 flume experiments with fixed beds and 38 with mobile beds. On average f_mob is approximately equal to half of f_btm, which is about one‐quarter of f. Hence, f_mob is about one‐tenth of f. Predictive equations are developed for f_btf, f_btm and f_mob using dimensional analysis to identify the relevant independent variables and regression analysis to evaluate the coefficients associated with these variables. Values of f_mob are always positive which implies that mobile beds offer greater resistance to flow than do fixed beds. Evidently bed‐load grains colliding with mobile beds lose more momentum to the bed than do grains colliding with fixed beds. In other words, grain collisions with mobile beds are less elastic than those with fixed beds. Copyright © 2005 John Wiley & Sons, Ltd.     

Size characteristics of sediment in interrill overland flow on a semiarid hillslope,Southern Arizona

This study examines the size characteristics of sediment removed from a semiarid hillslope by interrill overland flow. Rainfall simulation experiments were conducted on a runoff plot 18 m wide and 35 m long established on a piedmont hillslope in southern Arizona. The top of the plot coincided with the hillslope divide, and its outlet was located within a shallow rill. Samples of runoff were obtained from two cross-sections located in the interrill portion of the plot upslope of the rill and from a calibrated flume through which was directed interrill overland flow reaching the bottom of the plot. Analyses of sediment contained in these samples showed that sediment in interrill flow is finer than the matrix soil. The fineness of the interrill sediment compared to the matrix soil appears to be due to the inability of interrill overland flow to transport the coarser fraction of the sediment supplied to it by raindrop detachment. This finding implies that the rate of soil erosion in interrill areas is not. as is commonly supposed, limited by the rate at which raindrops can detach sediment but by the rate at which they detach sediment of a size that the overland flow is competent to transport. The relative fineness of sediment eroded from this hillslope is consistent with other evidence for the recent evolution of shrub-covered hillslopes in southern Arizona.     

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.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

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.     

Tracing sediment movement in interrill overland flow on a semi-arid grassland hillslope using magnetic susceptibility

Experiments were undertaken to determine the feasibility of tracing sediment movement in interrill overland flow. Crushed magnetite was introduced as a source-line 10 cm wide by 8 m long on a runoff plot 18 m wide by 29 m long located in southern Arizona. Initial magnetic susceptibilities along this source line, and along three transects located 0·25, 2·95 and 5 m downslope of the source-line, were measured. Movement of the magnetite in response to three rainfall simulation experiments was monitored. During the first two experiments, overland flow discharge was sampled at miniature flumes located along two cross sections on the plot downslope of the source-line, and at a supercritical flume at the plot outlet. Magnetic susceptibilities along the source-line and transects were measured after all three experiments. Results show that the magnetite moves very early in the experiments and that it reaches one of the flumes 2 m downslope of the source-line in 3 min. Most of the tracer moves a very short distance: 29·7 per cent is deposited within 25 cm of the source-line and only 2·2 per cent is deposited 2·95 m away. The deposition rate appears to decrease exponentially away from the source-line. Very little magnetite is recorded in the flow through the miniature flumes: in general it makes up less than 1 per cent of the total sediment load. No temporal pattern in these percentages is observed. Magnetite appears to be an effective tracer of sediment movement in interrill overland flow, though its higher density than natural soil may affect its detachment and transport.     

Effects of patterned Artemisia capillaris on overland flow velocity under simulated rainfall

Vegetation cover is an important factor for erosion control. Laboratory‐simulated rainfall experiments were conducted to quantify the effectiveness of patchy distributed Artemisia capillaris in retarding overland flow velocity. Simulated storms (60, 90, 120, and 150 mm h^?1) were applied on a bare plot (CK) and four different plant patterns, a banded pattern perpendicular to the slope direction (BP), a single long strip parallel to slope direction (LP), small patches distributed like a checkerboard (SP1), and small patches distributed like a letter “X” (SP2). All treatments had three replicates. Each plot underwent two sets of experiments, intact plant plots and root plots (the above‐ground parts were removed, only roots were reserved), respectively. Results showed that flow velocity increased with rainfall intensity, and the lower slope velocity (V_l) was higher than the upper slope velocity (V_u). The removal of grass shoots increased flow velocity. Compared with bare soil plot, intact plants reduced mean flow velocity by 14%–60%, whereas the reduction declined to <40% for the root plots. BP and both SP treatments performed more effectively than LP in retarding flow velocity, whereas no significant differences were identified between BP and SP. The contributions of A. capillaris shoots and roots to the reductions in flow velocity under different rainfall intensities were different. The shoots made greater contribution of 53%–97% at 60 and 90 mm h^–1, and the roots contributed more (51%–81%) at 120 and 150 mm h^–1. Runoff and sediment rate had significant (p < 0.05) linear correlations with mean flow velocity. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

The performance of grass filter strips in controlling high‐concentration suspended sediment from overland flow under rainfall/non‐rainfall conditions

There is little information on the performance of vegetative filter strips (VFS) in filtering high‐concentration sediment from subcritical overland flow. Flume experiments on simulated grass strips were conducted using combinations of three slope gradients (3°, 9° and 15°), five 1‐m‐wide slope positions (from upslope to downslope), two flow rates (60 and 20 L min^‐1 m^‐1) and sediment concentrations of 100–300 kg m^‐3 under simulated rainfall and non‐rainfall conditions. The results showed that sediment deposition efficiency increased with VFS width as a power function. Rainfall significantly reduced sediment deposited within VFS. Higher sediment concentration corresponded to a larger sediment deposition load but reduced deposition efficiency. Flow rate had a negative effect on deposition efficiency but no effect on deposition load. Sediments were more easily deposited at the upper slope position than downslope, and the upper slope position had a higher percentage of coarse sediments. The deposited sediment had significantly greater median diameters (D₅₀) than the inflow sediment. A greater proportion of coarse sediments larger than 25 µm in diameter were deposited, and particles smaller than 1 µm and of 10–25 µm had a better deposition performance than particles of 1–10 µm. Rainfall reduced the deposited sediment D₅₀ at a slope gradient of 3° and had no significant influence on it at 9° or 15°. A higher sediment concentration led to a smaller D₅₀ of the deposited sediment. Rainfall had no significant effect on overland flow velocity. Both the deposited sediment load and D₅₀ decreased with increasing flow velocity, and flow velocity was the most sensitive factor impacting sediment deposition. The results from this study should be useful to control sediment flowing into rivers in areas with serious soil erosion. Copyright © 2013 John Wiley & Sons, Ltd.     

How nature uses overland flow regime to maximize equilibrium detention storage and flood attenuation

Equilibrium detention storage is an important parameter as it has a proportional effect on flood attenuation. In this paper, based on the kinematic wave theory, a working formula for the equilibrium detention storage of an overland plane with upstream inflow has been derived. Since the flow regime over a concrete plane can vary throughout the entire range laminar to turbulent, this case has been selected to examine the effect of flow regime on the equilibrium detention storage. In the examination, the derived formula has been applied to four flow regimes: (a) laminar, (b) transitional, (c) near turbulent, and (d) turbulent. The examination shows that for planes with a small discharge, laminar flow gives the maximum detention storage. For planes with a medium discharge, transitional flow gives the maximum detention storage, and for planes with a large discharge, near turbulent flow gives the maximum detention storage. The flow regime can cause more than two‐fold increase in detention storage. All these results can be attributed to the respective flow resistance, and have been endorsed with analyses of the water surface profile and the rising limb of the hydrograph. Finally, relating the results to real‐life situations, it shows that the flow regime that gives the maximum detention storage is also the dominating flow regime in nature. Hence, extraordinarily, the flow regimes that exist in nature in fact provide maximum flood attenuation. Copyright © 2008 John Wiley & Sons, Ltd.     

An experimental study on the evolution law of roll wave parameters in overland flow

Roll waves commonly occur in overland flow and have an important influence on the progress of soil erosion on slopes. This study aimed to explore the evolution and mechanism of roll waves on steep slopes. The potential effects of flow rate, rainfall intensity and bed roughness on the laws controlling roll wave parameters were investigated. The flow rates, rainfall intensities and bed roughness varied from 5 to 30 L/min, 0 to 150 mm/h, and 0.061 to 1.700 mm, respectively. The results indicate that roll waves polymerize significantly along the propagation path, and bed roughness and rainfall affect the generation and evolution of roll waves. The wave velocity, length and height decreased with bed roughness, whereas the wave frequency increased with increasing bed roughness under fixed flow rate and rainfall intensity conditions. Rainfall increased the wave velocity and wavelength and decreased the wave frequency. The wave velocity, height and wavelength tended to increase with an increasing flow rate. Rainfall promoted the generation of roll waves, whereas bed roughness had the opposite effect. The generation of roll waves is closely related to the Froude number (Fr) and flow resistance. In this experiment, the range of the Reynolds number for the roll waves generated in the laminar region was 142–416, and the range of the flow resistance coefficient was 0.64–4.85. The critical value of the Fr for flow instability in the laminar region was approximately 0.57. Exploring the generation and evolution law of roll waves is necessary for understanding the processes and dynamic mechanisms of slope soil erosion.     

High spatial resolution mapping of surface velocities and depths for shallow overland flow

Point measurements of flow rate, depth or velocity are not sufficient to validate overland flow models, particularly when the interaction of the water with the soil surface creates a complex flow geometry. In this study, we present the coupling of two techniques obtaining spatial data of flow depths and surface velocity measurements for water depths as low as 1 mm. Overland flow experiments were performed in the laboratory at various flow rates and slopes on two surfaces. The first surface was 120 cm by 120 cm showing three undulations of sinusoidal shape with an amplitude of 1 cm and a wavelength of 20 cm, while the second was a 60 cm by 60 cm moulded reproduction of a seedbed with aggregates up to 2 cm in size. Large scale particle image velocimetry (LSPIV) was used for velocity measurements with a sub‐centimetre spatial resolution. An instantaneous‐profile laser scanner was used to map flow depths with a sub‐millimetre spatial resolution. A sensitivity analysis of the image processing of the LSPIV showed good robustness of the method. Comparison with measurements performed with hot film anemometer and salt velocity gauge showed that LSPIV surface velocities were representative of the flow. Water depths measured with the laser scanner were also in good agreement with single‐point measurements performed with a dial indicator. Spatially‐distributed flow rates could be computed by combining both presented techniques with a mean relative error less than 20%. Copyright © 2012 John Wiley & Sons, Ltd.