A new mixing equation for bed material composition in bed form dominant conditions

Vertical mixing beneath the bed surface and its effect on sediment composition has long been underestimated. This paper proposes a mixing equation to illustrate the temporal and spatial variations of bed material composition for bed form dominant rivers. A continuous mass conservation equation for elevation-specific nonuniform bed material composition is initially presented. A vertical mixing equation is finally derived as a diffusion-type partial differential equation with a source term charact...     

Lithology‐controlled evolution of stream bed sediment and basin‐scale sediment yields in adjacent mountain watersheds,Idaho, USA

The composition, grain‐size, and flux of stream sediment evolve downstream in response to variations in basin‐scale sediment delivery, channel network structure, and diminution during transport. Here, we document downstream changes in lithology and grain size within two adjacent ~300 km² catchments in the northern Rocky Mountains, USA, which drain differing mixtures of soft and resistant rock types, and where measured sediment yields differ two‐fold. We use a simple erosion–abrasion mass balance model to predict the downstream evolution of sediment flux and composition using a Monte Carlo approach constrained by measured sediment flux. Results show that the downstream evolution of the bed sediment composition is predictably related to changes in underlying geology, influencing the proportion of sediment carried as bedload or suspended load. In the Big Wood basin, particle abrasion reduces the proportion of fine‐grained sedimentary and volcanic rocks, depressing bedload in favor of suspended load. Reduced bedload transport leads to stronger bed armoring, and coarse granitic rocks are concentrated in the stream bed. By contrast, in the North Fork Big Lost basin, bedload yields are three times higher, the stream bed is less armored, and bed sediment becomes dominated by durable quartzitic sandstones. For both basins, the geology‐based mass balance model can reproduce within ~5% root‐mean‐square error the composition of the bed substrate using realistic erosion and abrasion parameters. As bed sediment evolves downstream, bedload fluxes increase and decrease as a function of the abrasion parameter and the frequency and size of tributary junctions, while suspended load increases steadily. Variable erosion and abrasion rates produce conditions of variable bed‐material transport rates that are sensitive to the distribution of lithologies and channel network structure, and, provided sufficient diversity in bedrock geology, measurements of bed sediment composition allow for an assessment of sediment source areas and yield using a simple modeling approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical modeling of bedload and suspended load contributions to morphological evolution of the Seine Estuary (France)

This numerical modeling study (i) assesses the influence of the sediment erosion process on the sediment dynamics and subsequent morphological changes of a mixed-sediment environment, the macrotidal Seine estuary, when non-cohesive particles are dominant within bed mixtures (non-cohesive regime), and (ii) investigates respective contributions of bedload and suspended load in these dynamics. A three dimensional (3D) process-based morphodynamic model was set up and run under realistic forcings (including tide, waves, wind, and river discharge) during a 1-year period. Applying erosion homogeneously to bed sediment in the non-cohesive regime, i.e., average erosion parameters in the erosion law (especially the erodibility parameter, E₀), leads to higher resuspension of fine sediment due to the presence of coarser fractions within mixtures, compared to the case of an independent treatment of erosion for each sediment class. This results in more pronounced horizontal sediment flux (two-fold increase for sand, +30% for mud) and erosion/deposition patterns (up to a two-fold increase in erosion over shoals, generally associated with some coarsening of bed sediment). Compared to observed bathymetric changes, more relevant erosion/deposition patterns are derived from the model when independent resuspension fluxes are considered in the non-cohesive regime. These results suggest that this kind of approach may be more relevant when local grain-size distributions become heterogeneous and multimodal for non-cohesive particles. Bedload transport appears to be a non-dominant but significant contributor to the sediment dynamics of the Seine Estuary mouth. The residual bedload flux represents, on average, between 17 and 38% of the suspended sand flux, its contribution generally increasing when bed sediment becomes coarser (can become dominant at specific locations). The average orientation of residual fluxes and erosion/deposition patterns caused by bedload generally follow those resulting from suspended sediment dynamics. Sediment mass budgets cumulated over the simulated year reveal a relative contribution of bedload to total mass budgets around 25% over large erosion areas of shoals, which can even become higher in sedimentation zones. However, bedload-induced dynamics can locally differ from the dynamics related to suspended load, resulting in specific residual transport, erosion/deposition patterns, and changes in seabed nature.     

A MATHEMATICAL MODEL FOR RESERVOIR SEDIMENTATION AND FLUVIAL PROCESSES

I. INTRODUCTIONAt present moot sediment transport models applied in engineering practice are based on equilibriumsediment transport approach, i. e. sediment--carrying Capacity is used to replace the actual sediment concentration (ref. 1 -- 9). However, the sediment--carrying capacity, in general, is not equal to sedimentconcentration, they may differ a lot especially for the case of reservoir sedimentation process and/orthe scouring process of river channel in the downstream of a reservoi…     

A three-dimensional numerical model of particulate transport for coastal waters

A Lagrangian marker particle in Eulerian finite difference cell solution to the three-dimensional incompressible mass transport equation was developed for predicting particulate transport in coastal and estuarine waters. Special features of the solution procedure include a finite difference grid network which translates horizontally and vertically with the mean particle motion and expands with the dispersive growth of the marker particle cloud. The cartesian vertical coordinate of the three-dimensional mass transport equation has been transformed, using instantaneous water column depth to allow adaptation to flow situations with a temporally and spatially varying bottom topography and free surface. Results from this model for turbulent diffusion and advection of a uniform plug flow of sediment in an unbounded uniform flow field with various sediment settling velocities were in excellent agreement with the corresponding analytic solutions. Using current information from a two-dimensional vertically averaged hydrodynamic's model, the model was utilized to predict the long term diffusion and advection of dilute neutrally and negatively buoyant suspended sediment clouds resulting from a hypothetical instantaneous release of dredge spoil waste at Brown's Ledge in Rhode Island Sound.     

Modelling the effects of Zostera noltei meadows on sediment dynamics: application to the Arcachon lagoon

A three-dimensional model has been modified to describe the complex interactions between hydrodynamics, sediment dynamics and biological parameters in the presence of Zostera noltei. The model treats seagrass leafs as flexible blades that bend under hydrodynamic forcing and alter the local momentum and turbulence fluxes and, therefore, the benthic shear conditions; these changes cause related changes to the mass balance at the boundary of the bed, in turn affecting the suspended matter in the column and ultimately primary productivity and the growth of the dwarf-grass. Modelling parameters related to the impact of Z. noltei to the local flow and to erosion and deposition rates were calibrated using flume experimental measurements; results from the calibration of the model are presented and discussed. The coupled model is applied in the Arcachon Bay, an area with high environmental significance and large abundance of dwarf-grass meadows. In the present paper, results from preliminary applications of the model are presented and discussed; the effectiveness of the coupled model is assessed comparing modelling results with available field measurements of suspended sediment concentrations and seagrass growth parameters. The model generally reproduces sediment dynamics and dwarf-grass seasonal growth in the domain efficiently. Investigations regarding the effects of the vegetation to the near-bed hydrodynamics and to the sediment suspension in the domain show that dwarf-grass meadows play an important part to velocity attenuation and to sediment stabilisation, with flow and suspended sediment concentrations damping, compared to an unvegetated state, to reach 35–50 and 65 %, respectively, at peak seagrass growth.     

AN ANALYTICAL SOLUTION FOR CALCULATING THE INITIATION OF SEDIMENT MOTION

This paper presents an analytical solution for calculating the initiation of sediment motion and the risk of river bed movement. It thus deals with a fundamental problem in sediment transport, for which no complete analytical solution has yet been found. The analytical solution presented here is based on forces acting on a single grain in state of initiation of sediment motion. The previous procedures for calculating the initiation of sediment motion are complemented by an innovative combination of optical surface measurement technology for determining geometrical parameters and their statistical derivation as well as a novel approach for determining the turbulence effects of velocity fluctuations. This two aspects and the comparison of the solution functions presented here with the well known data and functions of different authors mainly differ the presented solution model for calculating the initiation of sediment motion from previous approaches. The defined values of required geometrical parameters are based on hydraulically laboratory tests with spheres. With this limitations the derivated solution functions permit the calculation of the effective critical transport parameters of a single grain, the calculation of averaged critical parameters for describing the state of initiation of sediment motion on the river bed, the calculation of the probability density of the effective critical velocity as well as the calculation of the risk of river bed movement. The main advantage of the presented model is the closed analytical solution from the equilibrium of forces on a single grain to the solution functions describing the initiation of sediment motion.     

Bed-roughness investigation for a 2-D model calibration： the San Martin case study at Lower Parana

This paper illustrates how the acoustic Doppler current profiler （ADCP） and single-beam echo-sounder （SBES） recordings can be used for the calibration of existing software to assist in generalizing the morphodynamic processes in large rivers at key sites such as bifi.trcations and confluences. Calibration of the MIKE21C numerical model by the Danish Hydraulic Institute at the 25-km-long reach of Lower Paran~ near Rosario （Argentina） is presented. This reach includes two downstream confluences and two bifurcations. The model simulates a 2-D depth-averaged flow velocity and the related sediment fluxes to predict the bifurcation morphodynamics that affects the Paranh waterway. To investigate the river channel bathymetry, roughness, flow discharge allocation at bifurcations, suspended sediment concentration and grain size distributions, several instruments were used. These instruments included two ADCPs by Teledyne RDI working at frequencies of 600 and 1,200 kHz, a Sontek ADCP working at a frequency of 1,000 kHz and a SBES. The method to assess suspended sediment concentration and grain size distributions has been previously described. This paper focuses primarily on investigating dune morphology （by means of SBES depth measurements） and friction velocity （by means of ADCP profiling） to determine the river channel bed-roughness. The 2-D model results agree with observed values of bed-roughness, flow velocity and suspended sediment concentration distributions at the investigated sections, known data of water slope and total load of bed sediment are in good agreement with model results.     

Modeling equilibrium bed profiles of short tidal embayments

In many tidal embayments, bottom patterns, such as the channel-shoal systems of the Wadden Sea, are observed. To gain understanding of the mechanisms that result in these bottom patterns, an idealized model is developed and analyzed for short tidal embayments. In this model, the water motion is described by the depth- and width-averaged shallow water equations and forced by a prescribed sea surface elevation at the entrance of the embayment. The bed evolves due to the divergence and convergence of suspended sediment fluxes. To model this suspended-load sediment transport, the three-dimensional advection–diffusion equation is integrated over depth and averaged over the width. One of the sediment fluxes in the resulting one-dimensional advection–diffusion equation is proportional to the gradient of the local water depth. In most models, this topographically induced flux is not present. Using standard continuation techniques, morphodynamic equilibria are obtained for different parameter values and forcing conditions. The bathymetry of the resulting equilibrium bed profiles and their dependency on parameters, such as the phase difference between the externally prescribed M₂ and M₄ tide and the sediment fall velocity, are explained physically. With this model, it is then shown that for embayments that are dominated by a net import of sediment, morphodynamic equilibria only exist up to a maximum embayment length. Furthermore, the sensitivity of the model to different morphological boundary conditions at the entrance of the embayment is investigated and it is demonstrated how this strongly influences the shape and number of possible equilibrium bottom profiles. This paper ends with a comparison between the developed model and field data for the Wadden Sea’s Ameland and Frisian inlets. When the model is forced with the observed M₂ and M₄ tidal constituents, morphodynamic equilibria can be found with embayment lengths similar to those observed in these inlets. However, this is only possible when the topographically induced suspended sediment flux is included. Without this flux, the maximum embayment length for which morphodynamic equilibria can be found is approximately a third of the observed length. The sensitivity of the model to the topographically induced sediment flux is discussed in detail.     

NUMERICAL SIMULATIONS OF SCOUR AND DEPOSITION IN A CHANNEL NETWORK

1 INTRODUCTION Evolution of the river bed in alluvial channels has been studied by many researchers using analytical and numerical approaches. The use of analytical approach alone is insufficient for solving natural river engineering problems. With rapid growth in computer technology, numerical models have become a popular means for the study of mobile bed hydraulics. During the past decade, several numerical models have been developed. Most of the computer codes, such as HEC2SR (Si…     

Quantifying bed-related suspended load in gravel bed rivers through an analysis of the bedload-suspended load relationship

Suspended load transport can strongly impact ecosystems, dam filling and water resources. However, contrary to bedload, the use of physically based predicting equations is very challenging because of the complexity of interactions between suspended load and the river system. Through the analysis of extensive data sets, we investigated extent to which one or several river bed or flow parameters could be used as a proxy for quantifying suspended fluxes in gravel bed rivers. For this purpose, we gathered in the literature nearly 2400 instantaneous field measurements collected in 56 gravel bed rivers. Among all standard dimensionless parameters tested, the strongest correlation was observed between the suspended sediment concentration and the dimensionless bedload rate. An empirical relation between these two parameters was calibrated. Used with a reach average bedload transport formula, the approach allowed to successfully reproduce suspended fluxes measured during major flood events in seven gravel bed alpine rivers, morphodynamically active and distant from hillslope sources. These results are discussed in light of the complexity of the processes potentially influencing suspended load in a mountainous context. The approach proposed in this paper will never replace direct field measurements, which can be considered the only confident method to assess sediment fluxes in alpine streams; however, it can increment existing panel tools that help river managers to estimate even rough but not unrealistic suspended fluxes when measurements are totally absent. © 2019 John Wiley & Sons, Ltd.     

One-dimensional numerical simulation of non-uniform sediment transport under unsteady flows

One-dimensional numerical models are popularly used in sediment transport research because they can be easily programmed and cost less time compared with two- and three-dimensional numerical models. In particular, they possess greater capacity to be applied in large river basins with many tributaries. This paper presents a one-dimensional numerical model capable of calculating total-load sediment transport. The cross-section-averaged sediment transport capacity and recovery coefficient are addressed in the suspended load model. This one-dimensional model, therefore, can be applied to fine suspended loads and to hyperconcentrated flows in the Yellow River. Moreover, a new discretization scheme for the equation of unsteady non-uniform suspended sediment transport is proposed. The model is calibrated using data measured from the Yantan Reservoir on the Hongshui River and the Sanmenxia Reservoir on the Yellow River. A comparison of the calculated water level and river bed deformation with field measurements Shows that the improved numerical model is capable of predicting flow, sediment transport, bed changes, and bed-material sorting in various situations, with reasonable accuracy and reliability.     

Numerical Modeling of Non-Uniform Sediment Transport in River Channels

The mathematical model for simulating deformations of river channels composed of heterogeneous alluvium has been developed. The combination of shallow water equations and a three-layer model is used to describe the fluid flow and non-uniform sediment transport in bed (layer II) and suspended (layer III) loads. Changes in the fractional composition of unerodible bottom sediments (layer I) are also considered. The algorithm provides mass conservation for each fraction. The comparison of calculations results and experimental data (hydraulic washing of a desilting basin from sediments and armoring processes in heterogeneous soils) confirms the operability of the model. The model is applied to calculate the silting and hydraulic washes of the reservoir of a hydroelectric power station on a mountain river.     

Relationship between bed shear stress and suspended sediment concentration:annular flume experiments

In this study,annular flume experiments were carried out,using the sediment samples collected from the lower part of the inter-tidal zone at Xiaoyangkou,Jiangsu coast,China.The Ariathurai-Partheniades equation was used to determine the bed shear stress,by evaluating variations in the suspended sediment concentration within the water column.The derived relation between the bed shear stress and suspended sediment concentration shows that,at various stages of seabed erosion, suspended sediment concentration increases rapidly when the flow velocity is increased,but the pattern of change in the bed shear stress does not follow suit.At low concentrations,bed shear stress initially increases markedly with increasing flow velocity.However,when the concentration reaches an apparently critical level around 0.55 kg m"3,the rate of change in the bed shear stress abruptly slows down,or becomes almost constant,in response to further increases in the flow velocity.Results of experiments indicate that,from a critical level onward,suspended sediment concentration has a strong influence on the bed shear stress.     

Transport of suspended sediment under the dam‐break flow on an inclined plane bed of arbitrary slope

The problem of transport of suspended sediment after the break of a dam on an inclined bed is considered. To that end we use the shallow‐water approximation for arbitrary, constant slopes of the bottom, taking into consideration the effect of friction. The numerical technique and the frictional model are validated by comparison with available experimental data and asymptotic analytical solutions, with special attention to the numerical solution near the wetting front. The transport of suspended sediment down the inclined bed is obtained and discussed as a function of the slope of the bed for different values of the parameters characterizing the sediment and its transport properties. For sufficiently large times we always find the formation of roll waves near the water front, which affects the transport of sediments significantly. These strong oscillations are accurately computed with the numerical method used. The relative importance of the bed load (to the suspended load) sediment transport is also discussed as a function of the size of the sediment particles and the slope of the bed for different models on the initiation of sediment suspension from bed load. We also check the dilute sediment approach and characterize the conditions for its failure. Finally, the results of the present simplified model are intended to be used as tests of more complex numerical models. Copyright © 2007 John Wiley & Sons, Ltd.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

PREDICTION OF THE GRAIN SIZE OF SUSPENDED SEDIMENT： IMPLICATIONS FOR CALCULATING SUSPENDED SEDIMENT CONCENTRATIONS USING SINGLE FREQUENCY ACOUSTIC BACKSCATTER

Collection of samples of suspended sediment transported by streams and rivers is difficult and expensive. Emerging technologies, such as acoustic backscatter, have promise to decrease costs and allow more thorough sampling of transported sediment in streams and rivers. Acoustic backscatter information may be used to calculate the concentration of suspended sand-sized sediment given the vertical distribution of sediment size. Therefore, procedures to accurately compute suspended sediment size distributions from easily obtained river data are badly needed. In this study, techniques to predict the size of suspended sand are examined and their application to measuring concentrations using acoustic backscatter data are explored. Three methods to predict the size of sediment in suspension using bed sediment, flow criteria, and a modified form of the Rouse equation yielded mean suspended sediment sizes that differed from means of measured data by 7 to 50 percent. When one sample near the bed was used as a reference, mean error was reduced to about 5 percent. These errors in size determination translate into errors of 7 to 156 percent in the prediction of sediment concentration using backscatter data from 1 MHz single frequency acoustics.     

THREE-DIMENSIONAL NUMERICAL MODELING OF COHESIVE SEDIMENT TRANSPORT BY TIDAL CURRENT IN PEARL RIVER ESTUARV

IINTRODUCTIONEstUariesareprominentcoastalfeatUres.Estuariesareofgreateconomicssignificancetomankind.Attheseareas,manyharborsandwaterchannelshavetobebuiltforeconomicpurposes.ThedesignandconstrUctionofcoastalstrUctUresinestUariesrequireknowledgeofhydrodynamicsaswellassedimenttransportinsuchregions.ThenatUreofestuariesiscontrolledbyvariouscoastalhydrodynamicprocesses.Undertheactionofhydrodynamics,sedimentdepositionsorerosionswilloccurinestuariesornearcoastalstrUCtures.Tomaintainnavigati…     

ENTRAINMENT FUNCTION OF SUSPENDED SEDIMENT IN OPEN CHANNELS

1　ＩＮＴＲＯＤＵＣＴＩＯＮ　Ｒｉｖｅｒｆｌｏｗｓｗｉｔｈｓｕｓｐｅｎｄｅｄｓｅｄｉｍｅｎｔｓａｒｅａｍｏｎｇｔｈｅｍｏｓｔｉｍｐｏｒｔａｎｔｆａｃｔｏｒｓａｆｆｅｃｔｉｎｇｔｈｅｇｅｏｍｏｒｐｈｉｃｅｖｏｌｕｔｉｏｎ.Ａｋｅｙｆｅａｔｕｒｅｏｆｔｈｉ...     

Estimating sand concentrations using ADCP-based acoustic inversion in a large fluvial system characterized by bi-modal suspended-sediment distributions

Quantifying sediment flux within rivers is a challenge for many disciplines due, mainly, to difficulties inherent to traditional sediment sampling methods. These methods are operationally complex, high cost, and high risk. Additionally, the resulting data provide a low spatial and temporal resolution estimate of the total sediment flux, which has impeded advances in the understanding of the hydro-geomorphic characteristics of rivers. Acoustic technologies have been recognized as a leading tool for increasing the resolution of sediment data by relating their echo intensity level measurements to suspended sediment. Further effort is required to robustly test and develop these techniques across a wide range of conditions found in natural river systems. This article aims to evaluate the application of acoustic inversion techniques using commercially available, down-looking acoustic Doppler current profilers (ADCPs) in quantifying suspended sediment in a large sand bed river with varying bi-modal particle size distributions, wash load and suspended-sand ratios, and water stages. To achieve this objective, suspended sediment was physically sampled along the Paraná River, Argentina, under various hydro-sedimentological regimes. Two ADCPs emitting different sound frequencies were used to simultaneously profile echo intensity level within the water column. Using the sonar equation, calibrations were determined between suspended-sand concentrations and acoustic backscatter to solve the inverse problem. The study also analyzed the roles played by each term of the sonar equation, such as ADCP frequency, power supply, instrument constants, and particle size distributions typically found in sand bed rivers, on sediment attenuation and backscatter. Calibrations were successfully developed between corrected backscatter and suspended-sand concentrations for all sites and ADCP frequencies, resulting in mean suspended-sand concentration estimates within about 40% of the mean sampled concentrations. Noise values, calculated using the sonar equation and sediment sample characteristics, were fairly constant across evaluations, suggesting that they could be applied to other sand bed rivers. © 2018 John Wiley & Sons, Ltd.