Channel Sediment Transport Research Articles

Alongshore rip channel migration and sediment transport

The relationship between alongshore rip channel migration and sediment transport is investigated using time-averaged video images to identify the positions of rip channels at three sites along southern Monterey Bay, California, over 1–3 year periods. Daily migration rates are calculated for each site. The basic Coastal Engineering Research Center (CERC) formula for bulk alongshore sediment transport is used to estimate alongshore transport rates at each site, based on hourly input of directional wave spectra measured at an offshore buoy and shoaled and refracted to the nearshore. The wave transformation model is validated by comparison with measured nearshore directional wave spectra in 13 m water depth. An expanded CERC formulation that accounts for the additional forcing of alongshore wave height gradients is also tested but rejected. Correlation coefficients between daily rates of calculated alongshore transport and rip channel migration are low: 0.01, 0.44, and 0.38 at the Sand City, Stilwell, and Marina sites, respectively. The problem of digitization noise, associated with video-based rip channel identification, is discussed and several filtering methods are used to increase the signal to noise ratio and obtain improved correlations. Results indicate that higher frequency migration events (on time scales shorter than 8 days) will likely be obscured below the digitization “noise floor”, but that longer period oscillations (such as tidal and seasonal cycles) show up clearly in both rip migration and sediment transport datasets. Daily rates are cumulatively summed to generate time series of representative overall rip migration distance and net alongshore transport, and the cumulative summing procedure is shown to be similar to applying a low-pass filter. Mean migration is demonstrated to lag behind net transport in the summertime, suggesting that a minimum level of alongshore forcing is required to generate detectable rip migration. Rip channel migration and net alongshore transport have correlation coefficients that range from 0.76 to 0.94 at the three sites.

Flow movement and sediment transport in compound channels

There are many studies on the flow movement in compound channels, yet few are concerned with sediment transport. An experimental study on the flow movement and sediment transport in compound channels is presented. The experimental results indicate that the distribution of longitudinal velocity with depth in the main channel and the floodplains has a logarithmic component. The longitudinal velocity with flow depth in the interactive region does not obey a logarithmic distribution, but involves a wake function. In the boundary region, the longitudinal velocity obeys a parabolic distribution. In addition, based on the suspended sediment diffusion equation and the flow interaction between main channel and floodplains, expressions are derived to predict the lateral eddy viscosity and the sediment diffusion coefficients. Finally, an analytical solution for the lateral distribution of the depth-averaged velocity and sediment concentration in a compound channel is obtained. The results from the analytical solution agree well with experimentation.

PRELIMINARY RESULT OF DEVELOPMENT OF 3D NUMERICAL MODEL FOR MORPHOLOGICAL PROCESS

The paper presented a 3D numerical model for calculation of flow and sediment transport in channels. The flow is calculated by solving the Reynolds equation with hydrostatic assumption. Suspended-load is simulated by transport equation while bed-load is calculated using empirical formula. Bed transformation is obtained by solving an overall mass-balance equation. All equations are solved using finite volume method. The model is applied for study of bed deformation of VamNao river in the Mekong delta.

Lateral distributions of streamwise velocity in compound channels with partially vegetated floodplains

Natural rivers are commonly characterized by a main channel for primary flow conveyance and a floodplain, often partially covered with vegetation such as shrubs or trees, to carry extra flow during floods. The hydraulic resistance due to vegetation on the floodplain typically causes a further reduction of flow velocity and increases the velocity difference between the main channel and the floodplain. As a consequence a strong lateral shear layer leads to the exchange of mass and momentum between the main channel and floodplain, which in turn affects the overall channel conveyance and certain fluvial processes. The prediction of the lateral velocity distribution is important for many flood alleviation schemes, as well as for studies on sediment transport and dispersion in such channels. The present paper proposes a method for predicting the depth-averaged velocity in compound channels with partially vegetated floodplains, based on an analytical solution to the depth-integrated Reynolds-Averaged Navier-Stokes equation with a term included to account for the effects of vegetation. The vegetation is modelled via an additional term in the momentum equation to account for the additional drag force. The method includes the effects of bed friction, drag force, lateral turbulence and secondary flows, via four coefficients f, CD, λ & Γ respectively. The predicted lateral distributions of depth-averaged velocity agree well with the experimental data. The analytical solutions can also be used to predict the distribution of boundary shear stresses, which adds additional weight to the method proposed.

A balanced treatment of secondary currents, turbulence and dispersion in a depth-integrated hydrodynamic and bed deformation model for channel bends

This work deals with the formulation and numerical implementation of a two-dimensional mathematical and numerical model describing open channel hydrodynamics, sediment and/or scalar transport and riverbed evolution in curved channels. It is shown that a well balanced 2D model can predict flow features, sediment and scalar concentration, and bed elevation with an accuracy that is suitable for practical river engineering. The term “balanced” implies that important physical processes are modeled with a similar degree of complexity and exhaustiveness. The starting point of the model formulation is the assumption of self-similarity of vertical velocity profiles (horizontal velocities in the longitudinal and transverse directions), that are scaled by shear velocity and streamline curvature, both resolved by the model. The former is scaled by a bed-resistance coefficient that must be estimated or calibrated – as usual – on a application-specific basis, and the latter is computed by a new, grid-based but grid orientation independent, scheme that acts on the discrete solution. All processes, including bottom shear, momentum dispersion, scalar dispersion, turbulent diffusion, bed load, and suspended load, are modeled using physically based, averaged values of empirical or semi-empirical constants, and consistently with the assumed velocity profiles (and bed-generated turbulence). Bed deformation modeling can be implemented with either an equilibrium or non-equilibrium formulation of the Exner equation, depending on the adaptation length scale, which must be taken under consideration if significantly larger than the length scale of the spatial discretization. The governing equations are solved by high-resolution, unstructured-grid Godunov method, which is elsewhere tested and shown to be reliable and second-order accurate. Application of the model to laboratory test cases, using standard parameter values and previously reported bed-resistance coefficients, gives results comparable to many 2D and 3D models previously applied to the same cases, most part of which benefit from case-specific parameter tuning. There are obviously intrinsic limits to the descriptive ability of 2D models in river modeling, but the results of this study point to the utility and cost-effectiveness of a well-designed 2D model.

From flumes to rivers: Can sediment transport in natural alluvial channels be predicted from observations at the laboratory scale?

Doubt regarding the applicability of laboratory results to alluvial streams has led some to develop sediment transport predictors based solely on field data, and most current sediment transport formulae have typically been calibrated at least partially on field data. This paper examines the transferability of flume results to the field by exploring the extent to which a unified approach to the prediction of (1) flow regime, (2) depth, and (3) total sediment transport can be developed entirely with laboratory data. Relevance vector machine (RVM)‐based probabilistic models were constructed with only laboratory data, and their performances were tested against field data and found to be comparable with or better than currently available methods. Comparison of a laboratory‐trained RVM with a field‐trained RVM suggests that the prediction performances of the two models for unseen field data are not statistically different given the prediction uncertainty. For transferability, the choice of predictor variables is important with successful predictors being characterized by similar probability distribution in the laboratory and field data, e.g., as quantified by the Kullback‐Leibler divergence.

One‐dimensional numerical modeling of sediment transport and bed deformation in open channels

A one‐dimensional numerical model for simulating unsteady flow and sediment transport in open channels is presented and tested. The flow hydrodynamics is represented by the shallow water equations, and the bed morphodynamics is represented by the Exner equation and an additional equation describing the nonequilibrium sediment transport. Sediment size distribution is represented by the median grain diameter and the standard deviation, instead of the usual modeling with multiple particle size classes. Various methods for computing bed elevation changes at a cross section due to erosion or deposition of sediment are proposed and tested, including an innovative approach that relates the spatial pattern of erosion and deposition rates to boundary shear stress distribution, which is calculated by the Merged Perpendicular Method. An explicit finite difference scheme is employed for solving the water and sediment governing equations. The pertinence of the model is examined for two hypothetical cases. The model is then tested on one set of laboratory experiments on bed degradation under steady flow, showing excellent model data fit, and indicating that incorporating a nonequilibrium sediment transport equation into the model structure is an important element in reproducing the bed degradation process. Finally, the model is applied to simulate the morphological changes taking place in the Ha!Ha! River (Quebec) after the failure of the Ha!Ha! Dyke on July 1996. Relevant results can be obtained in terms of changes in longitudinal bed profile and cross‐sectional geometry as well as water levels, although some discrepancies are obtained between the simulated and surveyed cross‐sectional geometries, mainly because bank failure and channel widening are not modeled.

Two-Dimensional Nonequilibrium Noncohesive and Cohesive Sediment Transport Model

The purpose of this paper is to develop an unsteady 2D depth-averaged model for nonuniform sediment transport in alluvial channels. In this model, the orthogonal curvilinear coordinate system is adopted; the transport mechanisms of cohesive and noncohesive sediment are both embedded; the suspended load and bed load are treated separately. In addition, the processes of hydraulic sorting, armoring, and bed consolidation are also included in the model. The implicit two-step split-operator approach is used to solve the flow governing equations and the coupling approach with iterative method are used to solve the mass-conservation equation of suspended sediment, mass-conservation equation of active-layer sediment, and global mass-conservation equation for bed sediment simultaneously. Three sets of data, including suspension transport, degradation and aggradation cases for noncohesive sediment, and aggradation, degra- dation, and consolidation cases for cohesive sediment, have been demonstrated to show the rationality and accuracy of the model. Finally, the model is applied to evaluate the desilting efficiency for Ah Gong Diann Reservoir located in Taiwan to show its applicability.

A two-phase numerical model for suspended-sediment transport in estuaries

This paper presents a full 2-D X/Z numerical model for sediment transport in open channels and estuaries using a two-phase (fluid–solid particle) approach. The physical concept and the mathematical background of the model are given and test-cases have been carried out to validate the proposed model. In order to illustrate its feasibility for a real estuary, the model has been applied to simulate the suspended-sediment transport and the formation of turbidity maximum in the Seine estuary. The numerical results show that the main characteristics of estuarine hydro-sediment dynamics in the Seine estuary are in fact reproduced by the proposed model. A qualitative agreement between the numerical results and the actual observations has been obtained and is presented in this paper.

Geomorphic applications of stream‐gage information

AbstractIn the United States, several thousand stream gages provide what typically is the only source of continuous, long‐term streamflow and channel‐geometry information for the locations being monitored. In this paper, the geomorphic content of stream‐gage information, previous and potential applications of stream‐gage information in fluvial geomorphic research and various possible limitations are described. Documented applications include studies of hydraulic geometry, channel bankfull characteristics, sediment transport and channel geomorphic response to various types of disturbance. Potential applications include studies to determine the geomorphic effectiveness of large floods and in‐stream habitat change in response to disturbance. For certain applications, various spatial, temporal and data limitations may render the stream‐gage information of limited use; however, such information often is of considerable value to enable or enhance geomorphic investigations. Published in 2008 by John Wiley & Sons, Ltd.

A 3D non-hydrostatic model to predict flow and sediment transport in loose-bed channel bends

The validation and application of a model to solve for the flow, the sediment transport, and bed morphology changes in open channels is discussed. The non-hydrostatic model solves the three-dimensional (3D) incompressible, Reynolds-Averaged Navier–Stokes (RANS) equations in generalized curvilinear coordinates. The model uses adaptive grids in the vertical direction needed to account for changes in the free surface elevation and bed levels due to erosion/deposition at the bed as the code converges toward steady state (equilibrium conditions). The model is used to predict flow and sediment transport through straight and curved open channels, including the 140? open channel bend studied by Olesen (1985) and the 193? very sharp bend studied by Blanckaert (2002). An estimate of the peak non-dimensional shear stress at the banks is deduced function of the ratio between the radius of curvature of the bend and the channel width

Analysis of Suspended Sediment Transport in Open-Channel Flows: Kinetic-Model-Based Simulation

This paper presents a comprehensive analysis of suspended sediment transport in open channels under various flow conditions through a kinetic-model-based simulation. The kinetic model, accounting for both sediment-turbulence and sediment-sediment interactions, successfully represents experimentally observed diffusion and transport characteristics of suspended sediments with different densities and sizes. Without tuning any model coefficients, the nonmonotonic concentration distribution and the noticeable lag velocity with a negative value close to the wall are reasonably reproduced. Examination of flow conditions typical of suspension dominated rivers shows that the conventional method may overestimate or underestimate unit suspended-sediment discharge, depending on the Rouse number, sediment size, as well as shear velocity. The error may be less than 20% for dp 1.0 mm under typical flow conditions where shear velocity ranges from 1.0 to 12.5 cm∕s and flow depth ranges...

Influences of bedrock projections and large rocks on sediment transport in a mountainous gorge channel

Influences of bedrock projections and large rocks on sediment transport in a mountainous gorge channel

Channel flow competence and sediment transport in upland streams in southeast Australia

AbstractBedload transport data from planebed and step‐pool reach types are used to determine grain size transport thresholds for selected upland streams in southeast Australia. Morphological differences between the reach types allow the effects of frictional losses from bedforms, microtopography and bed packing to be incorporated into the dimensionless critical shear stress value. Local sediment transport data are also included in a regime model and applied to mountain streams, to investigate whether empirical data improve the delineation of reach types on the basis of dimensionless discharge per unit width (q*) and dimensionless bedload transport (qb*).Instrumented planebed and step‐pool sites are not competent to transport surface median grains (D50s) at bankfull discharge (Qbf). Application of a locally parametrized entrainment equation to the full range of reach types in the study area indicates that the majority of cascades, cascade‐pools, step‐pools and planebeds are also not competent at Qbf and require a 10 year recurrence interval flood to mobilize their D50s. Consequently, the hydraulic parameters of the regime diagram, which assume equilibrium conditions at bankfull, are ill suited to these streams and provide a poor basis of channel delineation. Modifying the diagram to better reflect the dominant transported bedload size (equivalent to the D16 of surface sediment) made only slight improvements to reach delineation and had greatest effect on the morphologies with smaller surface grain sizes such as forced pool‐riffles and planebeds. Likewise, the Corey shape factor was incorporated into the regime diagram as an objective method for adjusting a base dimensionless critical shear stress (τ*c50b) to account for lithologically controlled grain shape on bed packing and entrainment. However, it too provided only minor adjustments to reach type delineation. Copyright © 2007 John Wiley & Sons, Ltd.

A method for improving predictions of bed‐load discharges to reservoirs

AbstractEffective management options for mitigating the loss of reservoir water storage capacity to sedimentation depend on improved predictions of bed‐load discharges into the reservoirs. Most predictions of bed‐load discharges, however, are based on the assumption that the rates of bed‐load sediment availability equal the transport capacity of the flow, ignoring the spatio‐temporal variability of the sediment supply. This paper develops a semiquantitative method to characterize bed‐load sediment transport in alluvial channels, assuming a channel reach is non‐supply limited when the bed‐load discharge of a given sediment particle‐size class is functionally related to the energy that is available to transport that fraction of the total bed‐load. The method was applied to 22 alluvial stream channels in the USA to determine whether a channel reach had a supply‐limited or non‐supply‐limited bed‐load transport regime. The non‐supply‐limited transport regime was further subdivided into two groups on the basis of statistical tests. The results indicated the pattern of bed‐load sediment transport in alluvial channels depends on the complete spectrum of sediment particle sizes available for transport rather than individual particle‐size fractions represented by one characteristic particle size. The application of the method developed in this paper should assist reservoir managers in selecting bed‐load sediment transport equations to improve predictions of bed‐load discharge in alluvial streams, thereby significantly increasing the efficiency of management options for maintaining the storage capacity of waterbodies.

Spatially distributed modelling of soil erosion and sediment yield at regional scales in Spain

Initiated by the need to quantify erosion rates and the impacts of global changes on erosion, several attempts have been made to apply erosion models at regional scales. However, these models have often been directed towards on-site soil erosion estimates, emphasising sheet and rill erosion processes, and disregarding gully erosion, channel erosion and sediment transport. These models are therefore of limited use for the assessment of sediment yield, off-site impacts of erosion, and for the development of environmental management to control these impacts at regional scale. This study analyses and compares three spatially distributed models for the prediction of soil erosion and/or sediment yield at regional scales: the WATEM-SEDEM model that is based on the empirical Revised Universal Soil Loss Equation (RUSLE) in combination with a sediment transport equation, the physics-based Pan European Soil Erosion Risk Assessment model (PESERA), and a newly developed Spatially Distributed Scoring model (SPADS). The three models were applied to 61 Spanish drainage basins and model predictions were evaluated against data on measured reservoir sedimentation rates. Global data sets on land use, climate, elevation and soil characteristics were used as model input for WATEM-SEDEM and SPADS, whereas published soil erosion estimates of PESERA at 1 km 2 resolution were used directly. SPADS and WATEM-SEDEM provided best results after separate calibration for basins with a Sediment Delivery Ratio (SDR) higher than 5% and those with an SDR lower than 5%. In this way, SPADS explained 67% of variation in sediment yield, while WATEM-SEDEM explained 48% of the variation. PESERA represents a promising alternative to the use of empirical models at the regional scale as it can be applied to very diverse environments with little calibration. However, PESERA provides soil erosion rates and not sediment yield estimates. For most basins PESERA soil erosion rates vary between fifty and close to zero percent of total sediment yield. Two major factors may explain this discrepancy between modelled soil erosion rates and measured sediment yield. First, it may be that PESERA underestimates soil erosion under Mediterranean conditions, although PESERA soil erosion rates are of the same order of magnitude as erosion rates measured in erosion plot studies. Second, gully-, river channel erosion and sediment transport processes may be much more important than sheet- and rill erosion for regional scale sediment yield in these environments. These issues therefore require further attention in future model development. Although spatially lumped models provide better predictions of sediment yield at the basin scale, and while validation of the predicted spatial patterns of sources and sinks of sediment requires further research, spatially distributed models are expected to be of value to support management decisions regarding the assessment of on-site and off-site impacts of erosion at the regional scale.

Understanding confluence dynamics in the alluvial Ganga–Ramganga valley, India: An integrated approach using geomorphology and hydrology

Confluence dynamics in the Ganga–Ramganga valley in the western Ganga plains of India has been studied through systematic mapping of channel configuration using multi-date remote sensing images and topographic sheets for a period spanning nearly 100 years (1911–2000). The study has been supplemented with a detailed analysis of the channel morphology, hydrology and sediment transport characteristics of the different rivers. Our study indicates that new confluences have been created during this period and that the confluence points have moved both upstream and downstream on a historical time scale. Apart from major avulsions, other processes that have controlled the confluence movements include river capture, cut-offs and aggradation in the confluence area. River capture occurs through lateral bank erosion and migration, encroachment by the master stream and beheading of smaller rivers resulting in upstream movement of the confluence point. Another process which influences the upstream migration of the confluence is an increase in sinuosity of one of the channels near the confluence and then a cut-off. Aggradation in the confluence area and local avulsions of the primary channel in a multi-channel system seem to be the major process controlling the downstream movement of the confluence point. Analysis of channel morphology, hydrology and sediment budget for the study period supports our interpretations.

Dynamics of suspended sediment transport at field‐scale drain channels of irrigation‐dominated watersheds in the Sonoran Desert, southeastern California

AbstractSuspended sediment is a major source of pollution in irrigation‐dominated watersheds. However, little is known about the process and mechanisms of suspended sediment transport in drain channels directly connected to agricultural fields. This paper explains sediment dynamics using averaged 5 min flow discharge Q (m3 s−1) and suspended sediment concentration C (mg l−1) collected during one crop season in a small catchment containing a first‐order drain channel and its connected six agricultural fields within the Salton Sea watershed. The statistical properties and average trends of Q and C were investigated for both early (i.e. November) and late (i.e. January) stages of a crop season. Further in‐depth analysis on sediment dynamics was performed by selecting two typical single‐field irrigation events and two multiple‐field irrigation events. For each set of irrigation events, the process of suspended sediment transport was revealed by examining hydrograph and sediment graph responses. The mechanisms underlying suspended sediment transport were investigated by analysing the types of corresponding hysteresis loop. Finally, sediment rating curves for both hourly and daily data at early and late stages and for the entire crop season were established to seek possible sediment‐transport predictive model(s). The study suggests that the complicated processes of suspended sediment transport in irrigation‐dominated watersheds require stochastic rather than deterministic forecasting. Copyright © 2007 John Wiley & Sons, Ltd.

Tidal and meteorological forcing of sediment transport in tributary mudflat channels

Field observations of flow and sediment transport in a tributary channel through intertidal mudflats indicate that suspended sediment was closely linked to advection and dispersion of a tidal salinity front. During calm weather when tidal forcing was dominant, high concentrations of suspended sediment advected up the mudflat channel in the narrow region between salty water from San Francisco Bay and much fresher runoff from the small local watershed. Salinity and suspended sediment dispersed at similar rates through each tidal inundation, such that during receding ebbs the sediment pulse had spread spatially and maximum concentrations had decreased. Net sediment transport was moderately onshore during the calm weather, as asymmetries in stratification due to tidal straining of the salinity front enhanced deposition, particularly during weaker neap tidal forcing. Sediment transport by tidal forcing was periodically altered by winter storms. During storms, strong winds from the south generated wind waves and temporarily increased suspended sediment concentrations. Increased discharge down the tributary channels due to precipitation had more lasting impact on sediment transport, supplying both buoyancy and fine sediment to the system. Net sediment transport depended on the balance between calm weather tidal forcing and perturbations by episodic storms. Net transport in the tributary channel was generally offshore during storms and during calm weather spring tides, and onshore during calm weather neap tides.

Feedbacks between erosion and sediment transport in experimental bedrock channels

AbstractNatural bedrock rivers flow in self‐formed channels and form diverse erosional morphologies. The parameters that collectively define channel morphology (e.g. width, slope, bed roughness, bedrock exposure, sediment size distribution) all influence river incision rates and dynamically adjust in poorly understood ways to imposed fluid and sediment fluxes. To explore the mechanics of river incision, we conducted laboratory experiments in which the complexities of natural bedrock channels were reduced to a homogenous brittle substrate (sand and cement), a single sediment size primarily transported as bedload, a single erosion mechanism (abrasion) and sediment‐starved transport conditions. We find that patterns of erosion both create and are sensitive functions of the evolving bed topography because of feedbacks between the turbulent flow field, sediment transport and bottom roughness. Abrasion only occurs where sediment impacts the bed, and so positive feedback occurs between the sediment preferentially drawn to topographic lows by gravity and the further erosion of these lows. However, the spatial focusing of erosion results in tortuous flow paths and erosional forms (inner channels, scoops, potholes), which dissipate flow energy. This energy dissipation is a negative feedback that reduces sediment transport capacity, inhibiting further incision and ultimately leading to channel morphologies adjusted to just transport the imposed sediment load. Copyright © 2007 John Wiley & Sons, Ltd.