Influence of Turbulence on Bed Load Sediment Transport

This paper summarizes the results of an experimental study on the influence of an external turbulence field on the bed load sediment transport in an open channel. The external turbulence was generated by (1) a horizontal pipe placed halfway through the depth h; (2) a series of grids with a clearance of about one-third of the depth from the bed, and extending over a finite length of the flume; and (3) a series of grids with a clearance in the range (0.1–1.0)h from the bed, but extending over the entire length of the flume. Two kinds of experiments were conducted: plane-bed experiments and ripple-covered-bed experiments. In the former case, the flow in the presence of the turbulence generator was adjusted so that the mean bed shear stress was the same as in the case without the turbulence generator in order to single out the effect of the external turbulence on the sediment transport. In the ripple-covered-bed case, the mean and turbulence quantities of the streamwise component of the velocity were measured, and the Shields parameter, due to skin friction, was determined. The Shields parameter, together with the RMS value of the streamwise velocity fluctuations, was correlated with the sediment transport rate. The sediment transport increases markedly with increasing turbulence level.     

Evaluation and Improvement of Bed Load Discharge Formulas based on Helley–Smith Sampling in an Alpine Gravel Bed River

Bed load discharge formulas have been evaluated by analyzing them in relation to measured Helley–Smith data for the gravel-bedded armored Drau River, Austria. Comparison of calculations with measurements leads to ranking of the formulas that depends on the evaluation parameters. The choice of formula is made with respect to our specific aims: the investigation of individual floods requires a different approach from that of long-term budgets. Formula performance is consistently improved when conditions for the threshold of motion are modified according to data measured up on the initiation of motion. Formulas such as those reported by Parker in 1990, Zanke in 1999, and Sun and Donahue in 2000 are capable of coping with partial transport, which is commonly found in Alpine rivers. These formulas therefore provide encouraging results, particularly after the introduction of modifications. The augmentation of field measurements, even if limited in scope, considerably improves the performance of bed load discharge formulas.     

Measurement of Bed Load Velocity using an Acoustic Doppler Current Profiler

A new technique has been developed to measure the apparent velocity of bed load (va) using an acoustic Doppler current profiler. The technique involves estimating the bias in bottom tracking due to a moving bottom. Mean va measured at sampling stations in the gravel-bed Fraser River correlated well (r2 = 0.93,?n = 9) with mean bed load transport rates measured using conventional samplers. Mean va was also correlated (r2 = 0.44,?n = 19) with boundary shear stress estimated by a log-law fit to the mean velocity profile. Estimates of va from individual 5 s ensemble averages were extremely variable: the coefficient of variation for a sampling station ranged from 1.0 to 6.4, and 25 min of sampling were required to achieve stable estimates of the mean and coefficient of variation (within 5% error). Variance was due to both real temporal variability of transport and measurement error. The mechanisms that produce this variability are discussed and preliminarily quantified.     

Multiple Linear Regression Model for Total Bed Material Load Prediction

A new total bed material load equation that is applicable for rivers in Malaysia was developed using multiple linear regression analyses. A total of 346 hydraulic and sediment data were collected from nine natural and channelized rivers having diverse catchment characteristics in Malaysia. The governing parameters were carefully selected based on literature survey and field experiments, examined and grouped into five categories namely mobility, transport, sediment, shape, and flow resistance parameters. The most influential parameters from each group were selected by using all possible regression model method. The suitable model selection criteria namely the R-square, adjusted R-square, mean square error, and Mallow’s Cp statistics were employed. The accuracy of the derived model is determined using the discrepancy ratio, which is a ratio of the calculated values to the measured values. The best performing models that give the highest percentage of prediction from the validation data were chosen. In general, the newly derived model is best suited for rivers with uniform sediment size distribution with a d50 value within the range of 0.37–4.0 mm and performs better than the commonly used Graf, Yang, and Ackers–White total bed material load equations.     

Unified View of Sediment Transport by Currents and Waves. I: Initiation of Motion, Bed Roughness, and Bed-Load Transport

Attention is given to the properties of sediment beds over the full range of conditions (silts to gravel), in particular the effect of fine silt on the bed composition and on initiation of motion (critical conditions) is discussed. High-quality bed-load transport data sets are identified and analyzed, showing that the bed-load transport in the sand range is related to velocity to power 2.5. The bed-load transport is not much affected by particle size. The prediction of bed roughness is addressed and the prediction of bed-load transport in steady river flow is extended to coastal flow applying an intrawave approach. Simplified bed-load transport formulas are presented, which can be used to obtain a quick estimate of bed-load transport in river and coastal flows. It is shown that the sediment transport of fine silts to coarse sand can be described in a unified model framework using fairly simple expressions. The proposed model is fully predictive in the sense that only the basic hydrodynamic parameters (depth, current velocity, wave height, wave period, etc.) and the basic sediment characteristics (d10, d50, d90, water temperature, and salinity) need to be known. The prediction of the effective bed roughness is an integral part of the model.     

Validation of Existing Bed Load Transport Formulas Using In-Sewer Sediment

Granular sediment in pipe inverts has been reported in a number of sewer systems in Europe. Given the range of flow conditions and particle characteristics of inorganic sewer sediments the mode of transport may normally be considered as bed load. Current commercial software for modeling the erosion and transport of sediments in sewer pipes still utilizes well-known, or modified versions of transport equations that were derived for transport of noncohesive sediment in alluvial streams. In this paper the performances of the equations of Ackers and White (originally developed for the transport of river sediments) and of May (derived from laboratory pipe experiments) are examined against two separate data sets. One set is from laboratory erosion experiments on sewer sediment obtained in Paris. A second data set has bed load transport rate measurements recorded in a sewer inlet pipe. The formulas were selected because of their widespread use in the prediction of in-sewer sediment transport both in commercial software and in the latest United Kingdom design guidance for new sewers. The results indicated that both the relationships performed poorly, even in such well-controlled conditions. These formulas have significant difficulties in predicting the erosion thresholds and fractional transport rates for non-uniformly sized in-sewer sediments. An empirical formula to adjust the threshold of motion for individual grain size fractions was developed which significantly improved predictions. Although such techniques have been used in gravel bed rivers, the threshold adjustment function for in-sewer deposits was significantly different from these previously published for fluvial gravels, indicating that a direct transfer of fluvial relationships to sewers may be inappropriate without further research.     

Wash Load and Bed-Material Load Transport in the Yellow River

It has been the conventional assumption that wash load is supply limited and is only indirectly related to the hydraulics of a river. Hydraulic engineers also assumed that bed-material load concentration is independent of wash load concentration. This paper provides a detailed analysis of the Yellow River sediment transport data to determine whether the above assumptions are true and whether wash load concentration can be computed from the original unit stream power formula and the modified unit stream power formula for sediment-laden flows. A systematic and thorough analysis of 1,160 sets of data collected from 9 gauging stations along the Middle and Lower Yellow River confirmed that the method suggested by the conjunctive use of the two formulas can be used to compute wash load, bed-material load, and total load in the Yellow River with accuracy.     

Effect of Seepage-Induced Nonhydrostatic Pressure Distribution on Bed-Load Transport and Bed Morphodynamics

Bed-load transport is commonly evaluated in the condition of a hydrostatic pressure distribution of the flow field; while this condition is reasonable for quasi-steady, quasi-uniform rectilinear flows, it cannot be satisfied in a large variety of flow conditions, i.e., near an obstacle as in the case of a bridge pier. The dimensionless Shields number, which contains the assumption of a hydrostatic pressure distribution in its denominator, therefore cannot be strictly applied to evaluate bed-load transport in all the configurations where nonhydrostatic pressure distributions are observed. In the present work, a generalization of the Shields number is proposed for the case of nonhydrostatic pressure distribution produced by groundwater flow. Experiments showing the effects of vertical groundwater flow on the bed morphodynamics are presented. The comparison between the experimental observations and numerical results, obtained by means of a morphodynamic model which employs the new formulation of the Shields number, suggests that the proposed generalization of the Shields number is able to account the effect of the nonhydrostatic pressure distribution on the bed-load transport.     

Bed-Load Transport Equation for Sheet Flow

When open-channel flows become sufficiently powerful, the mode of bed-load transport changes from saltation to sheet flow. Where there is no suspended sediment, sheet flow consists of a layer of colliding grains whose basal concentration approaches that of the stationary bed. These collisions give rise to a dispersive stress that acts normal to the bed and supports the bed load. An equation for predicting the rate of bed-load transport in sheet flow is developed from an analysis of 55 flume and closed conduit experiments. The equation is ib = ω where ib = immersed bed-load transport rate; and ω = flow power. That ib = ω implies that eb = tan?α = ub/u, where eb = Bagnold’s bed-load transport efficiency; ub = mean grain velocity in the sheet-flow layer; and tan?α = dynamic internal friction coefficient. Given that tan?α ≈ 0.6 for natural sand, ub ≈ 0.6u, and eb ≈ 0.6. This finding is confirmed by an independent analysis of the experimental data. The value of 0.60 for eb is much larger than the value of 0.12 calculated by Bagnold, indicating that sheet flow is a much more efficient mode of bed-load transport than previously thought.     

Bedload Transport in Alluvial Channels

Hydraulic, sediment, land-use, and rock-erosivity data of 22 alluvial streams were used to evaluate conditions of bedload transport and the performance of selected bedload-transport equations. Transport categories of transport-limited (TL), partially transport-limited (PTL), and supply-limited (SL) were identified by a semiquantitative approach that considers hydraulic constraints on sediment movement and the processes that control sediment availability at the basin scale. Equations by Parker et al. in 1982, Schoklitsch in 1962, and Meyer-Peter and Muller in 1948 adequately predicted sediment transport in channels with TL condition, whereas the equations of Bagnold in 1980, and Schoklitsch, in 1962, performed well for PTL and SL conditions. Overall, the equation of Schoklitsch predicted well the measured bedload data for eight of 22 streams, and the Bagnold equation predicted the measured data in seven streams.     

Stochastic Prediction of Sediment Transport in Sand-Gravel Bed Rivers

Classical deterministic bedload transport predictors are applied to sand-gravel bed rivers. The turbulent bed shear stress is modeled according to a probability distribution to obtain realistic bedload transport rates at incipient motion. In extending the predictors to stochastic predictors for nonuniform sediment, many parameters that represent near-bed turbulence and the particle size distribution must be chosen. The parameters that give realistic results are chosen by analyzing the results of a new experimental flume dataset with relatively large water depths. Choosing other combinations of parameters may give equal total bedload transport rates, but at the cost of large errors in fractional transport rates. Attention is given to the hiding-exposure phenomenon and a hindrance effect related to nonuniform sediment. Validation based on two independent field datasets shows that successful predictions of particle sizes near the threshold for motion are feasible using the stochastic approach, while the deterministic approach gives successful predictions well above incipient motion.     

Microforms in Gravel Bed Rivers: Formation, Disintegration, and Effects on Bedload Transport

This research aims to advance current knowledge on cluster formation and evolution by tackling some of the aspects associated with cluster microtopography and the effects of clusters on bedload transport. The specific objectives of the study are (1) to identify the bed shear stress range in which clusters form and disintegrate, (2) to quantitatively describe the spacing characteristics and orientation of clusters with respect to flow characteristics, (3) to quantify the effects clusters have on the mean bedload rate, and (4) to assess the effects of clusters on the pulsating nature of bedload. In order to meet the objectives of this study, two main experimental scenarios, namely, Test Series A and B (20 experiments overall) are considered in a laboratory flume under well-controlled conditions. Series A tests are performed to address objectives (1) and (2) while Series B is designed to meet objectives (3) and (4). Results show that cluster microforms develop in uniform sediment at 1.25 to 2 times the Shields parameter of an individual particle and start disintegrating at about 2.25 times the Shields parameter. It is found that during an unsteady flow event, effects of clusters on bedload transport rate can be classified in three different phases: a sink phase where clusters absorb incoming sediment, a neutral phase where clusters do not affect bedload, and a source phase where clusters release particles. Clusters also increase the magnitude of the fluctuations in bedload transport rate, showing that clusters amplify the unsteady nature of bedload transport. A fourth-order autoregressive, autoregressive integrated moving average model is employed to describe the time series of bedload and provide a predictive formula for predicting bedload at different periods. Finally, a change-point analysis enhanced with a binary segmentation procedure is performed to identify the abrupt changes in the bedload statistic characteristics due to the effects of clusters and detect the different phases in bedload time series using probability theory. The analysis verifies the experimental findings that three phases are detected in the bedload rate time series structure, namely, sink, neutral, and source.     

Two-Dimensional Total Sediment Load Model Equations

An unsteady total load equation is derived for use in depth-averaged sediment transport models. The equation does not require the load to be segregated a priori into bed and suspended but rather automatically switches to suspended load, bed load, or mixed load depending on a transport mode parameter consisting of local flow hydraulics. Further, the sediment transport velocity, developed from available data, is explicitly tracked, and makes the equation suitable for unsteady events of sediment movement. The equation can be applied to multiple size fractions and ensures smooth transition of sediment variables between bed load and suspended load for each size fraction. The new contributions of the current work are the consistent treatment of sediment concentration in the model equation and the empirical definition of parameters that ensure smooth transitions of sediment variables between suspended load and bed load.     

Bedload Layer Thickness and Disturbance Depth in Gravel Bed Streams

A field investigation in ten gravel bed stream reaches determined that substrate disturbance depth associated with a moving bedload layer was a small multiple of the bed surface D90. Disturbance depth during plane bed transport of coarse, heterogeneous mixtures appeared similar in magnitude to particle exchange depth and moving layer thickness. Maximum disturbance depth was distributed approximately uniformly over the most active areas of the streambed when local scour and fill were negligible. The distribution upper bound was the smaller of approximately 1.5 times the competent grain size or twice the surface D90, and was invariant with flow strength once the largest grains present were mobilized. Disturbance depth did not scale with grain sizes smaller than D50 when larger grains were mobilized. Thicker traction carpets were not predicted to occur because much larger shear stresses then observed naturally were needed to mobilize two or more layers of the bed simultaneously. Bedload transport rate in coarse streambeds is suggested to increase primarily with mobile fraction of bed surface area and grain velocity, than with layer thickness.     

Exponential Formula for Bedload Transport

An exponential formula that does not involve the concept of the critical shear stress is derived in this study for computing bedload transport rates. The formula represents well various experimental data sets ranging from the weak transport to high shear conditions. Comparisons of the present study are also made with many previous bedload formulas commonly cited in the literature.     

Measurements of Sediment Erosion and Transport with the Adjustable Shear Stress Erosion and Transport Flume

Soil and sediments play an important role in water management and water quality. Issues such as water turbidity, associated contaminants, reservoir sedimentation, undesirable erosion and scour, and aquatic habitat are all linked to sediment properties and behaviors. In situ analysis is necessary to develop an understanding of the erosion and transport of sediments. Sandia National Laboratories has recently patented the Adjustable Shear Stress Erosion and Transport (ASSET) Flume that quantifies in situ erosion of a sediment core with depth while affording simultaneous examination of transport modes (bedload versus suspended load) of the eroded material. Core erosion rates and ratios of bedload to suspended load transport of quartz sediments were studied with the ASSET Flume. The erosion and transport of a fine-grained natural cohesive sediment were also observed. Experiments using quartz sands revealed that the ratio of suspended load to bedload sediment transport is a function of grain diameter and shear stress at the sediment surface. Data collected from the ASSET Flume were used to formulate a novel empirical relation for predicting the ratio of bedload to suspended load as a function of shear stress and grain diameter for noncohesive sediments.     

Suspended Wash Load Transport of Nonuniform Sediments

The results of an experimental study on transport of suspended wash load through a coarse-bed stream are presented. The experiments were conducted under different concentrations of fine suspended sediment (wash load of uniform size, 0.064 mm diameter) and with three different coarse-bed sediments: two having uniform sizes and one with nonuniform size distribution. For any equilibrium concentration of wash load in suspension, a definite proportion of the wash material was observed to be present within the bed material. No difference is found in this regard between wash load and suspended load transport. Therefore, the relationship, as stated by Samaga et al., for the parameter representing sheltering—exposure and interference effects in the suspended load transport of nonuniform sediments was applied in a modified form by using the present data and the data collected from the literature.     

Transition between Two Bed-Load Transport Regimes: Saltation and Sheet Flow

In the saltation regime where bed-shear stress is low, bed load moves by sliding, rolling, and saltating along the bed, while in the sheet-flow regime where bed-shear stress is high, it travels by a combination of saltation and sheet flow. In this paper a theoretical model is developed for predicting the onset of the sheet-flow regime as shear stress increases. This model is based on a new variable Pb representing the proportion of grains on the bed that are entrained as bed load. The model yields the equation Pb = 2.56θG3 in which G = 1?θc/θ, θ = dimensionless bed-shear stress; and θc = critical value of θ at which grains begin to move. The equation shows that θt, which is the value of θ at the onset of the sheet-flow regime and is assumed to occur when Pb = 1, is around 0.5 with the exact value controlled by θc. For example, when θc = 0.045, θt = 0.52. The theoretical model is verified by performing a nonlinear regression analysis on data from 285 flume experiments. Additional flume experiments with a high-speed video (HSV) system result in consistent values of θ for the onset of the sheet-flow regime, which support the theoretical model. The HSV images further reveal that: (1) the sheet-flow regime is characterized by granular sheets or laminations; and (2) a zone of mixed saltation and rolling grains exists not only in the saltation regime but also in the sheet-flow regime.     

Field Assessment of Alternative Bed-Load Transport Estimators

Measurement of near-bed sediment velocities with acoustic Doppler current profilers (ADCPs) is an emerging approach for quantifying bed-load sediment fluxes in rivers. Previous investigations of the technique have relied on conventional physical bed-load sampling to provide reference transport information with which to validate the ADCP measurements. However, physical samples are subject to substantial errors, especially under field conditions in which surrogate methods are most needed. Comparisons between ADCP bed velocity measurements with bed-load transport rates estimated from bed-form migration rates in the lower Missouri River show a strong correlation between the two surrogate measures over a wide range of mild to moderately intense sediment transporting conditions. The correlation between the ADCP measurements and physical bed-load samples is comparatively poor, suggesting that physical bed-load sampling is ineffective for ground-truthing alternative techniques in large sand-bed rivers. Bed velocities measured in this study became more variable with increasing bed-form wavelength at higher shear stresses. Under these conditions, bed-form dimensions greatly exceed the region of the bed ensonified by the ADCP, and the magnitude of the acoustic measurements depends on instrument location with respect to bed-form crests and troughs. Alternative algorithms for estimating bed-load transport from paired longitudinal profiles of bed topography were evaluated. An algorithm based on the routing of local erosion and deposition volumes that eliminates the need to identify individual bed forms was found to give results similar to those of more conventional dune-tracking methods. This method is particularly useful in cases where complex bed-form morphology makes delineation of individual bed forms difficult.     

Effect of Coarse Surface Layer on Bed-Load Transport

Existing bed-load transport formulas may overestimate the transport rate in mountain rivers by two orders of magnitude or more. Recently published field data sets provide an opportunity to take a fresh look at the bed-load transport relationship and it is hypothesized that the overestimate is due to a failure to account for the effect of a coarse surface layer of bed material inhibiting the release of fine subsurface material. Bed-load transport is determined as gs = aρ(q?qc) where q=water discharge per unit width; qc=critical value for initiation of bed material movement; ρ=water density; and a=coefficient. The gs/q relationship is typically piecewise linear, characterized by two transport phases with, respectively, low and high rates of change. Twenty-one flume and 25 field data sets were used to quantify the relationship for Phase 2. The flume data confirm the dependence of a on S1.5, where S=channel slope, in agreement with earlier studies. The field data additionally show that a varies inversely with the degree of bed armoring, given by the ratio of surface to subsurface bed material size. The finding is consistent with the hypothesis and suggests the need to account for the bed material supply limitation in the bed-load transport formula. However, the available data are not entirely sufficient to rule out an alternative dependency, or codependency, on flow resistance. The critical conditions for initiation of Phase 2 transport are also quantified as a function of bed material size and channel slope. The resulting set of equations allows a more accurate estimation of Phase 2 bed-load transport rates. However, the equations are empirical and should be restricted for use within the range of conditions used in their development, to determine mean rather than instantaneous transport rates and to determine bulk transport rates, not transport by size fraction.