Modeling Nonuniform Suspended Sediment Transport in Alluvial Rivers

Problems and difficulties in modeling sediment transport in alluvial rivers arise when one uses the theory of equilibrium transport of uniform sediment to simulate riverbed variation. A two-dimensional mathematical model for nonuniform suspended sediment transport is presented to simulate riverbed deformation. Through dividing sediment mixture into several size groups in which the sediment particles are thought to be uniform, the nonuniformity and the exchange between suspended sediment and bed material are considered. The change of concentration along the flow direction, size redistribution, and cross-sectional bed variation can then be described reasonably well by the model. In simulating the flow field with big dry-wet flats, moving boundary problems are solved very well by introducing a so-called finite-slot technique. Verification with laboratory data shows that the model has a good ability to simulate channel bed variations. Last, the model was applied to a real alluvial river system. Variables such as water level, sediment concentration, suspended sediment size distribution, and riverbed variation were obtained with encouraging results.     

Monitoring Suspended Sediment Dynamics Using MBES

This technical note describes use of a multibeam echosounder (MBES) to quantify the dynamics of suspended sediment in a large open channel. A methodology is detailed that uses the backscatter magnitude from the MBES water-column data to adjust the magnitude of sonar returns for the various sonar settings, spatially and temporally average the data to account for the random nature of acoustic backscatter from the suspended sediment, and calibrate the processed data with direct samples. A case study of flow at the confluence of the Rio Paraná and Rio Paraguay, Argentina, where there is a distinct turbidity difference along the mixing interface of the two flows, is used to demonstrate the unique capabilities of MBES to quantify sediment concentrations and dynamics within the water column.     

Aspect Ratio to Maximize Sediment Transport in Rigid Bank Channels

The continuity equation, Manning’s equation, Einstein’s wall correction procedure and sediment transport equations are combined to indicate channel aspect ratios which maximize sediment transport for a given water discharge in rigid-bank trapezoidal and rectangular channels with fixed slope. Higher aspect ratios are required to maximize sediment transport for channels conveying bed load than for those with a dominant suspended load. A total load equation predicts optimum aspect ratios lying in between those for bed load and suspended load channels. The equations imply that the optimum aspect ratio increases markedly as the channel bank to channel bed roughness ratio increases. The resulting optimum ratios are smaller than the aspect ratios of many natural rivers.     

Prediction of Sediment Load Concentration in Rivers using Artificial Neural Network Model

An artificial neural model is used to estimate the natural sediment discharge in rivers in terms of sediment concentration. This is achieved by training the network to extrapolate several natural streams data collected from reliable sources. The selection of water and sediment variables used in the model is based on the prior knowledge of the conventional analyses, based on the dynamic laws of flow and sediment. Choosing an appropriate neural network structure and providing field data to that network for training purpose are addressed by using a constructive back-propagation algorithm. The model parameters, as well as fluvial variables, are extensively investigated in order to get the most accurate results. In verification, the estimated sediment concentration values agree well with the measured ones. The model is evaluated by applying it to other groups of data from different rivers. In general, the new approach gives better results compared to several commonly used formulas of sediment discharge.     

Lateral Variations in Suspended Sediment Concentration over Dunes

A series of laboratory flume experiments were conducted to characterize lateral variations of suspended sediment over dunes. Forty experimental runs were made using flow depths of 0.328 and 0.127?m in a 1.2-m-wide channel. The Froude number was 0.47 and the median diameter of the bed material was 0.52?mm. Point samples of suspended sediment and depth-integrated values calculated using acoustic backscatter data were collected simultaneously at two lateral positions at spacings of 0.41, 0.20, 0.10, and 0.05?m. Root mean square differences between the paired samples increased from 30 to 94% and from 20 to 45% of the mean concentration for point samples and depth-integrated samples, respectively, as the lateral spacing between samples was increased. Lateral variations in the concentration of suspended sediment should be considered when designing sediment sampling strategies over dune beds.     

Fluctuations of Suspended Sediment Concentration and Turbulent Sediment Fluxes in an Open-Channel Flow

A complex problem of turbulent-sediment interactions in an open-channel flow is approached experimentally, using specially designed field experiments in an irrigation canal. The experimental design included synchronous measurements of instantaneous three-dimensional (3D) velocities and suspended sediment concentration using acoustic Doppler velocimeters (ADV) and a water sampling system. Various statistical measures of sediment concentration fluctuations, turbulent sediment fluxes, and diffusion coefficients for fluid momentum and sediment are considered. Statistics, fractal behavior, and contributions of bursting events to vertical fluxes of fluid momentum and sediment are evaluated using quadrant analysis. It has been found that both turbulence and sediment events are organized in fractal clusters which introduce additional characteristic time and spatial scales into the problem and should be further explored. It is also shown that Barenblatt’s theory of sediment-laden flows appears to be a good approximation of experimental data.     

Total Load Transport Formula for Flow in Alluvial Channels

A user-friendly total bed-material load transport formula for flow in alluvial channels under equilibrium transport conditions has been developed based on dimensional analysis. The main advantages of this formula are its ease of computation, accuracy in prediction, and the wide range of application. The total sediment discharge gt is computed directly and is linearly related to the new total load transport parameter, TT. The latter involves variables that can be easily measured in field conditions, i.e., flow depth, mean flow velocity, energy slope, median sediment size and density, and water temperature. The factor of proportionality k in the formula has been checked for a wide range of hydraulic conditions and it remains a constant equal to 12.5. Comparisons between the computed and measured total sediment discharge indicate that the predictions are good.     

Diagonal Cartesian Method for the Numerical Simulation of Flow and Suspended Sediment Transport over Complex Boundaries

There is increasing demand for simulation tools of flow and suspended sediment transport over complex boundaries in hydraulic engineering. The diagonal Cartesian method, which approximates complex boundaries using both Cartesian grid lines and diagonal lines segments, is presented in the paper to simulate the complex boundaries of two-dimensional shallow-water turbulence equations and nonequilibrium suspended sediment transport equation. The method, which utilizes cell-centered nodes on a nonstaggered grid, uses boundary velocity information at the wall boundary to avoid the specification of water level. An enlarged finite-difference method is introduced for momentum and suspended sediment equations on the complex boundary. This paper describes an application of the diagonal Cartesian method to calculate the tidal current and suspended sediment concentration of Quanzhou Bay in the Fujian province of China. The results show that the method predicts the flow and suspended sediment concentration well, and the calculations agree well with the measurement.     

Critical Shear Stress of Bimodal Sediment in Sand-Gravel Rivers

A new model for the critical shear stress and the transport of graded sediment is presented. The model is based on the size distribution of the bed surface and can be used to compute sediment transport rates in numerical simulations with an active layer model. This model makes a distinction between unimodal and bimodal sediments. It is assumed that all size fractions of unimodal sediments have the same critical shear stress while there is selective transport for the gravel fractions of bimodal sediments. A recently published laboratory transport data set is used to calibrate our model.     

Density Stratification Effects in Sand-Bed Rivers

In this paper the effects of density stratification in sand-bed rivers are studied by the application of a model of vertical velocity and concentration profiles, coupled through the use of a turbulence closure that retains the buoyancy terms. By making the governing equations dimensionless, it is revealed that the slope is the additional dimensionless parameter introduced by inclusion of the buoyancy terms. The primary new finding is that in general density stratification effects tend to be greater in large, low-slope rivers than in their smaller, steeper brethren. Under high flow conditions the total suspended load and size distribution of suspended sediment can be significantly affected by density stratification, and should be accounted for in any general theory of suspended transport.     

Influence of Coherent Flow Structures on the Dynamics of Suspended Sediment Transport in Open-Channel Flow

The influence of suspended sediments on coherent flow structures has been studied by simultaneously measuring the longitudinal and vertical components of the instantaneous velocity vector and the instantaneous suspended particle concentration with an acoustic particle flux profiler. The measurements were carried out in clear water and in particle-laden open-channel flows. In both cases, they clearly show the predominance of ejection and sweep phases that are part of a burst cycle. The analysis further demonstrates the importance of the ejection and sweep phases in sediment resuspension and transport. Ejections pick up the sediment at the bed and carry it up through the water column close to the surface. It is shown that ejections and sweeps are in near equality in the near-bottom layer, whereas ejections clearly dominate in the remaining water column. The implications of these results for sediment transport dynamics are discussed.     

Variability in Suspended-Sediment Concentration over Mobile Sand Beds

Lateral and longitudinal variability in suspended-sediment concentration was investigated over dune and upper-regime plane beds under steady, uniform flow in a laboratory flume. A much greater sampling time was required to reliably predict mean concentration in a dune bed than for an upper-regime plane bed. Laterally separated, simultaneous samples over both dune and upper-regime plane beds were poorly correlated.     

Modeling Suspended Sediment during Construction in Great Barrier Reef World Heritage Area

A marina was constructed in the Great Barrier Reef World Heritage Area in close proximity to coral reefs that could be damaged by excess turbidity generated during construction. Since there was uncertainty about both the fate of suspended sediments and their effect on corals, initial water quality constraints were set very conservatively. In order to better understand the movement of suspended sediment during construction, a numerical model study was commissioned using three-dimensional, numerical, hydrodynamic, and Lagrangian particle tracking models. The study was successful in: (1) increasing the understanding of and reducing the uncertainty of sediment dispersal patterns under a range of common forcing conditions; (2) testing the variation in suspended sediment concentrations over sensitive areas for two different outfall locations; (3) offering evidence that a good choice in outfall locations will reduce the threat to corals; and importantly (4) presenting the results in a way that enhanced understanding by nontechnical reef managers. This final result was achieved by creating movies of sediment movement that clearly demonstrated the complex hydrodynamic processes involved with near-coastal water currents. Specific model results showed: (1) that a more seaward outfall increases effluent dispersal away from sensitive areas; (2) the highest concentrations of effluent over sensitive sites occur during no wind and neap tide conditions; and (3) prevailing southeast winds advect effluent offshore, away from sensitive sites.     

Formation and Evolution of the Jiuduansha Shoal over the Past 50 Years

The formation and geomorphic evolution of the Jiuduansha Shoal are investigated over the past 50 years in this paper. The sedimentary processes, formation, and geomorphic evolution of the Jiuduansha Shoal in the Yangtze River Estuary, eastern China, are analyzed based on digitized sea maps from 1945 to 2001, satellite images from 1975 to 2001, and field survey data in the spring-neap tidal cycle in the dry and flood seasons in 2003. The suspended sediment concentration and hydrodynamics of the North Passage and South Passage during the dry and flood seasons in the spring-neap tidal cycle in 2003 were investigated, and relations between the North Passage and South Passage of the Jiuduansha Shoal are analyzed. Results show that seasonal and spring-neap tidal cycle variations in sediment and water discharge from the drainage basin correspond to erosion and deposition on the Jiuduansha Shoal. The results and data can provide useful information for the management of the Yangtze River Estuary and restoration of the estuarine tidal flat ecology.     

Numerical Model of Sediment Pulses and Sediment-Supply Disturbances in Mountain Rivers

Sediment pulses in rivers can result from many mechanisms including landslides entering from side slopes and debris flows entering from tributaries. Artificial sediment pulses can be caused by the removal of a dam. This paper presents a numerical model for the simulation of gravel bedload transport and sediment pulse evolution in mountain rivers. A combination of the backwater and quasi-normal flow formulations is used to calculate flow parameters. Gravel bedload transport is calculated with the surface-based bedload equation of Parker in 1990. The Exner equation of sediment continuity is used to express the mass balance at different grain size groups and lithologies, as well as the abrasion of gravel. The river is assumed to have no geological controls such as bedrock outcrops and immobile boulder pavements. The results of nine numerical experiments designed to study various key parameters relevant to the evolution of sediment pulses are reported here. Results of the numerical runs indicate that the evolution of sediment pulses in mountain rivers is dominated by dispersion rather than translation. Here dispersion is an expression for the observation that a sediment pulse aggrades both upstream and downstream of its apex whereas its amplitude decreases in time. The results also indicate that grain abrasion is an important and yet often neglected mechanism in removing the excess sediment associated with pulse inputs from some mountain rivers.     

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.     

Effect of Shallowness on Sediment Deposition in Recirculation Zones

Various hydrodynamic and sediment deposition experiments were conducted on a purpose-built, wide, shallow water table using a range of incoming flow and bed roughness conditions, including five flow rates, different inflow sediment concentrations, and different artificial roughness coefficients generated by four transverse string intervals. A recirculation zone was created behind a 0.5-m-wide backward-facing step on the left side of the table. Surface velocities, sediment concentrations, and processes and features of sediment deposition in recirculation zones were measured. Experimental results indicate that (1)?as the value of the shallow stability parameter increased from 0.033 to 0.47, the length of the recirculation zone decreased to less than half, but the corresponding area and volume of sediment deposition increased about four times, and (2)?it is the strength of the mixing layer rather than the size of the recirculation zone that dominates the process of sediment deposition as well as the final deposition feature.     

Modeling Noncohesive Suspended Sediment Transport in Stream Channels Using an Ensemble-Averaged Conservation Equation

The governing conservation equation for the transport of noncohesive suspended sediment in erodible channels is recognized as a stochastic partial differential equation due to the uncertainties in the parameters, and a deterministic ensemble-averaged equation is developed. Variables in this one-dimensional equation are represented as averaged quantities, and their covariances are also taken into account. Lateral inflows and deposition and entrainment of sediment are incorporated in the formulation. A hypothetical test problem is constructed to examine the model behavior. Manning’s coefficient, bed slope and bottom width are taken as the primary random parameters. Results from the solution of the ensemble-averaged equation are compared to results from Monte Carlo simulations. For comparison purposes, predicted values are also obtained by solving the deterministic transport equation without the covariance terms. It is found that predictions obtained from this latter approach deviate significantly from Monte Carlo simulation results. On the other hand, the ensemble-averaged predictions compare favorably to the Monte Carlo simulation results indicating that this promising technique needs further exploration.     

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.     

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.