Field studies of sediment transport by overland flow

Field studies on sandy soils of the Cottenham Series in mid-Bedfordshire show that the mean annual rate of sediment transport by overland flow on an 11° mid-slope is 98 g cm¹. The feasibility of using sediment transport equations to predict erosion by overland flow on a storm basis is examined by comparing the observed values of sediment yield with values predicted by four sediment transport equations and a regression equation which relates soil loss to runoff energy and rainfall energy. An expression combining Engelund's sediment transport capacity equation and the Manning equation for flow velocity, as modified by Savat for disturbed flow, best reflects field conditions. Although there is a significant correlation (r = 0.69; n 30) between the observed and predicted values using this expression, the coefficient of determination is too low for predictive purposes. Reasons for this are presented.     

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.     

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.     

Empirical relations for the sediment transport capacity of interrill flow

Experiments were carried out in order to measure the sediment transport capacity of interrill flow with and without rainfall and to relate the transport capacity to selected hydraulic parameters, such as effective stream power and shear velocity. Different sediments were used in order to study the effect of grain size. The proposed relationships show considerable variations with grain size and there is only a minor effect of rainfall on the transport capacity which also seems to be grain size dependent. The proposed relationships can be used to predict sediment transport capacity of interrill flow and can therefore contribute to the development of physically based erosion models.     

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.     

Effects of vegetation on flow conveyance and sediment transport capacity

In-stream and riparian vegetation may significantly affect flow and sediment transport in vegetated channels. A hydraulic model has been developed in this paper to compute the flow discharge in channels with rigid and flexible vegetation under emergent and submerged conditions. An empirical formula has also been presented to determine the bed-load discharge in vegetated channels. The model has been tested against experimental and field data available in the literature. The computed flow discharge and bed-load transport rate agree well with the measured data.     

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.     

Processes of accelerated pluvial erosion on desert hillslopes modified by vehicular traffic

Accelerated pluvial erosion on hillslopes modified by off-road vehicles (ORVs) is analysed using results from 50 rainfall simulation experiments conducted in the Mojave Desert, California. Sediment yield from 1 m² hillslope plots subjected to intense, 20-minute rainfalls is typically increased 10 to 20-fold following ORV use. Salient effects of vehicle traffic, which reduce infiltration, increase runoff sediment transport efficiency, and enhance gully formation, are further studied by combining simple theoretical relations with experimental data. This analysis helps identify factors controlling erosion on natural desert hillslopes, as well as those used by ORVs. Erosion of natural or vehicle-used desert surfaces is heavily influenced by runoff hydraulics. Calculated Darcy-Weisbach friction factors decrease by an average of 13-fold following vehicular slope modification, whereas runoff Reynolds numbers increase by an average of 5 1/2-fold. The capacity of overland flow to transport sediment is related to runoff power and its degree of localization, which usually increase considerably following ORV activity; however, the ability of overland flow to move large grains (competency) is related to a combination of factors not always systematically influenced by ORV use. Kinematic runoff routing, which is used to extrapolate experimental results to longer slope lengths, leads to the suggestion that the hydraulic roughness of desert hillslopes strongly influences their erosional behaviour.     

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.     

Relation of sediment transport capacity to stone cover and size in rain‐impacted interrill flow

This study examines how the sediment transport capacity of interrill overland flow varies with stone cover and stone size at two flow intensities. Six series of flume experiments were conducted on two slopes (2° and 10°) with stones of three sizes (28·0, 45·5 and 91·3 mm) serving as roughness elements. Bed sediment size, water discharge and simulated rainfall intensity were the same in all experiments. It was found (1) that transport capacity is positively related to stone size, with the relation becoming stronger as stone cover increases and flow intensity decreases; and (2) that transport capacity is negatively related to stone cover at the high flow intensity and curvilinearly related to stone cover at the low flow intensity. The curvilinear relations are concave‐upward with the lowest transport capacities occurring at stone covers between 0·40 and 0·60. The highest transport capacities are found at stone covers of 0 and 1, with the transport capacity being greater at the former stone cover than at the latter. Copyright © 2000 John Wiley & Sons, Ltd.     

LES-DEM simulations of sediment transport

In this work, a fully-coupled Computational Fluid Dynamics (CFD) model and Discrete Element Method (DEM) are used to simulate a unidirectional turbulent open-channel flow over the full range of sediment transport regimes. The fluid and particles are computed on separate grids using a dual-grid formulation to maintain consistency and avoid instability issues. The results of coupling the dispersed phase to a multiphase flow solver that uses volume-averaged Navier-Stokes equations are compared to those obtained from coupling through drag to a single flow solver. The current work also examines the applicability and limitations of lumping particles as a representative particle to reduce the cost of simulations. Insight to the impact of different turbulent events to the entrainment of particles is also given. The simulation results of sediment transport from both coupling techniques show good agreement with empirical formulas in the bedload regime, but under-predict sediment transport in the suspended load regime. In the suspended load regime, using partial coupling, the rate of sediment transport was found to be under-predicted as compared to full-coupling. The deviation in results in the suspended load regime was found to increase with increases in the applied shear stress. Both coupling methods revealed the same effect on the friction factor where friction increases in the bedload regime and decreases in the suspended load regime reaching a maximum at the transition between regimes. This result is contrary to past studies which have shown a discrete jump in the friction factor at the transition. Lumping particles as representative particles is shown to reduce the simulation cost by more than a factor of 5 when using a scaling factor of 2. By doing a quadrant analysis on information obtained from particle and flow field results, it was found that most of the particles are entrained by more frequent sweep events.     

Sediment-transport competence of rain-impacted interrill overland flow

Laboratory experiments to determine the maximum size of sediment transported in shallow, rain-impacted flow were conducted in a recirculating flume 4·80 m long and 0·50 m wide. Rainfall intensities were varied between 51 and 138 mm h⁻¹, flow was introduced from a header tank into the flume at rates ranging from 0 to 0·64 l s⁻¹, and experiments were conducted on gradients between 3·5 and 10°. The following equation was developed: ML = (REFE)^1·6363 in which M is particle mass, L is distance moved in unit time (cm min⁻¹), RE is rainfall energy (J m⁻² s⁻¹) and FE is flow energy (J m⁻² s⁻¹). This equation can be used to predict sediment-transport competence of interrill overland flow. The equation is limited in its utility insofar as it has been developed using quartz grains and takes no account of variations in absorption of rain energy by natural ground surfaces. © 1998 John Wiley & Sons, Ltd.     

Transporting capacity of overland flow on plane and on irregular beds

In this paper the transporting capacity of thin flows, in the laminar and transitional flow regime, is studied. Experiments were carried out on irregular as well as on plane beds, using two totally different set-ups. The results of these two types of experiment were convergent. In both cases, sediment concentration was clearly related to grain shear velocity and unit stream power, expressed as the product of mean velocity and slope (Yang, 1973). The data agreed with those of Kramer and Meyer (1969). For a sandy bed, the unit stream power relationship was able to predict reasonably well the sediment concentrations measured on a mulched surface. For laminar and transitional flows, both the unit stream power and the shear velocity are related in the same way to slope and unit discharge. The unit stream power is a parameter which in particular can be very easily measured and might therefore become useful in obtaining a quick estimate of the transporting capacity of a thin flow. However, before a sediment transport equation for thin flows can be developed, more information is needed about the influence of the flow regime and grain size and density.     

Field measurement of the velocity of overland flow using dye tracing

A new method of computing the mean velocity of overland flow using dye tracing is proposed in which a specified cross-section is divided into zones of relatively uniform flow characteristics, termed partial sections. The mean surface velocity for each partial section is determined by timing the passage of injected dye, and this figure is multiplied by 0.67, 0.70, or 0.80, depending on whether the flow is laminar, transitional, or turbulent, to give mean velocity. The mean velocity for the entire cross-section is calculated by multiplying the mean velocity of each partial section by its cross-sectional area, summing the products for all the partial sections, and dividing by the total area of the cross-section. A field test shows that mean velocity derived in this manner differs from mean velocity derived by the discharge method (i.e. by dividing discharge measured volumetrically by cross-sectional area) by an average of only 7.07 per cent. Thus the partial-section technique appears to provide a reliable method of estimating mean velocity of overland flow.     

Availability and performance of sediment detachment and transport functions for overland flow conditions

Abstract

Soil erosion is a global environmental problem. To quantify water erosion rates at the field, hillslope or catchment scale, several spatially-distributed soil erosion models have been developed. The accuracy of those models depends largely on the sediment detachment and transport functions used, many of which were developed from empirical research. In this paper, the physical basis of the available sediment detachment and transport functions is reviewed, and their application boundaries determined. Well-known and widely-used sediment detachment and transport functions are discussed on the basis of composite force predictors, i.e. shear stress, stream power, unit stream power and effective stream power, and their suitability is elucidated based on information in the literature. It was found that only a few sediment detachment functions are available, and those have been poorly tested. Most erosion models ignore direct calculation of sediment detachment, but use the sediment transport capacity deficit approach to estimate detachment rate. Many more sediment transport functions are available that also tested better for overland flow conditions. However, our tests did not result in a single function that appeared to perform best under a range of experimental conditions. The unit stream power-based functions developed by Govers seem to be the most promising ones for water erosion modelling. It is therefore recommended to evaluate the performance of existing sediment transport functions with more detailed field and laboratory datasets.

Editor Z.W. Kundzewicz     

Improved numerical solvers for implicit coupling of subsurface and overland flow

Due to complex dynamics inherent in the physical models, numerical formulation of subsurface and overland flow coupling can be challenging to solve. ParFlow is a subsurface flow code that utilizes a structured grid discretization in order to benefit from fast and efficient structured solvers. Implicit coupling between subsurface and overland flow modes in ParFlow is obtained by prescribing an overland boundary condition at the top surface of the computational domain. This form of implicit coupling leads to the activation and deactivation of the overland boundary condition during simulations where ponding or drying events occur. This results in a discontinuity in the discrete system that can be challenging to resolve. Furthermore, the coupling relies on unstructured connectivities between the subsurface and surface components of the discrete system, which makes it challenging to use structured solvers to effectively capture the dynamics of the coupled flow. We present a formulation of the discretized algebraic system that enables the use of an analytic form of the Jacobian for the Newton–Krylov solver, while preserving the structured properties of the discretization. An effective multigrid preconditioner is extracted from the analytic Jacobian and used to precondition the Jacobian linear system solver. We compare the performance of the new solver against one that uses a finite difference approximation to the Jacobian within the Newton–Krylov approach, previously used in the literature. Numerical results explores the effectiveness of using the analytic Jacobian for the Newton–Krylov solver, and highlights the performance of the new preconditioner and its cost. The results indicate that the new solver is robust and generally outperforms the solver that is based on the finite difference approximation to the Jacobian, for problems where the overland boundary condition is activated and deactivated during the simulation. A parallel weak scaling study highlights the efficiency of the new solver.     

Test of scale modelling of sediment transport in steady unidirectional flow

In two steady uniform flows at different physical scales in a small open channel, with variables characterizing flow, sediment, and fluid adjusted for dynamic similitude by means of four dimensionless modelling parameters (a Reynolds number, a Froude number, a density ratio, and a length ratio), measured frequency distributions of height, spacing, and migration rate of current ripples were almost identical when scaled, thus verifying that exact Reynolds-Froude modelling of loose-sediment transport is valid and workable. Modelling should be valid as well for a wide range of other transport conditions in the same kind of flow, because no additional kinds of forces or effects would be present in transport of loose grains in modes other than as ripples. In scaled-down modelling, a scale ratio of 2.5 is attainable without recourse to exotic fluids by use of water at 85°C to model natural flows at 10°C.     

Unsteady overland flow on flat surfaces induced by spatial permeability contrasts

Lateral redistribution of surface water in patchy arid ecosystems has been hypothesized to contribute to the maintenance of vegetation patches through the provision of a water subsidy from bare sites to vegetated sites. Such runon-runoff processes occur during Hortonian runoff events on topographically sloping ground. Surface flow redistribution may also occur on topographically flat ground if the presence of the vegetation patch creates a contrast in infiltration rate, leading to a free-surface gradient in ponded water. The precise dynamics and the eco-hydrologic role of this process has resisted complete theoretical treatment to date. Here the overland flow equations are modified to account for the presence of vegetation situated over a flat surface. The resulting model is solved numerically to determine whether this mechanism could influence the spatial partitioning of water in patchy arid ecosystems. Assumptions made about infiltration processes and overland flow in existing eco-hydrologic models of patchy and patterned arid ecosystems are evaluated in comparison to the solution of the ‘full’ coupled Saint-Venant equations with various infiltration models. The results indicate that the optimization of vegetation spatial patch scales with respect to water redistribution may be determined by the size of the infiltration redistribution length L over which the presence of an infiltration contrast perturbs baseline infiltration behavior.     

Linear and non-linear approaches to predict the Darcy-Weisbach friction factor of overland flow using the extreme learning machine approach

In modeling of overland flow and erosion, the overland flow friction factor (f), is a crucial factor. Due to the importance of a good understanding of f and its variability, the current study aimed to investigate the capability of non-linear approaches to estimate the Darcy-Weisbach friction factor of overland flow and its components (sediment transport, wave, form, and grain friction factors) through the Extreme Learning Machine (ELM) approach. Four datasets were used herein which were obtained from flume experiments done by different researchers. In order to investigate the effects of different parameters on the friction factor, numerous models consisting of various parameters were utilized to predict the friction factor using the ELM approach. The modeling procedure was established in two stages; the first stage aimed to model the overland flow friction factor and investigate the effect of the different parameters on the friction factor using non-linear separation via the ELM approach. In the second stage, the friction factor was linearly separated into different types of friction factors and then the separate components were estimated. Sensitivity analysis results confirmed the key role of Froude number (Fr) values for most of the models. On the other hand, the results obtained for estimated values of the friction factor were acceptable and outperformed available empirical approaches.     

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.