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.     

Modeling the preservation potential of tidal flat sedimentary records,Jiangsu coast,eastern China

A forward modeling approach is proposed to simulate the preservation potential of tidal flat deposits. The preservation potential is expressed as a function of net deposition rate and a factor that represents the vertical flux of suspended load, or seabed lowering during erosion periods associated with bedload transport. The model takes into account a number of geometric parameters of a tidal flat sediment system and sediment dynamic processes. The former includes high water level, total sediment supply, the annual rate of the supply, the ratio of mud to bilk sediment in the supply, the bed slope of the tidal flat profile, as well as the slope of the stratigraphic boundary; the latter includes spring-neap cycles of tidal water level changes, boundary layer processes, resuspension of fine-grained sediments, bedload transport due to tidal currents, and bed elevation changes in response to sediment movement. Using this model, numerical experiments are carried out for a tidal flat system on the Jiangsu coast, eastern China, with the input data being derived from literature and from a series of sediment cores collected along an onshore–offshore transect. The results show that the preservation potential is highest over the upper part of the inter-tidal zone and in the lower part of the sub-tidal zone, and lowest near mean sea level and at low water on springs. The preservation potential tends to decrease with the advancement of the shoreline. The bed slope, tidal current direction and resuspension intensity influence the spatial distributions of the preservation potential. An implication of these results is that the temporal resolution of the tidal flat record depends upon the location and depth within the deposit; this should be taken into account in the interpretation of sedimentary records. Further studies are required to improve the model, on the hydrodynamic processes associated with extremely shallow water depths, sediment dynamic modeling of bed slope and profile shape, and the combined action of tides and waves for sediment transport on tidal flats.     

2D stream hydrodynamic,sediment transport and bed morphology model for engineering applications

A 2D depth‐averaged hydrodynamic, sediment transport and bed morphology model named STREMR HySeD is presented. The depth‐averaged sediment transport equations are derived from the 3D dilute, multiphase, flow equations and are incorporated into the hydrodynamic model STREMR. The hydrodynamic model includes a two‐equation turbulence model and a correction for the mean flow due to secondary flows. The suspended sediment load can be subdivided into different size classes using the continuum (two‐fluid) approach; however, only one bed sediment size is used herein. The validation of the model is presented by comparing the suspended sediment transport module against experimental measurements and analytical solutions for the case of equilibrium sediment‐laden in a transition from a rigid bed to a porous bed where re‐suspension of sediment is prevented. On the other hand, the bed‐load sediment transport and bed evolution numerical results are compared against bed equilibrium experimental results for the case of a meander bend. A sensitivity analysis based on the correction for secondary flow on the mean flow including the effect of secondary flow on bed shear stresses direction as well as the downward acceleration effect due to gravity on transverse bed slopes is performed and discussed. In general, acceptable agreement is found when comparing the numerical results obtained with STREMR HySeD against experimental measurements and analytical solutions. Copyright © 2007 John Wiley & Sons, Ltd.     

Modeling river dune development and dune transition to upper stage plane bed

Large asymmetric bedforms known as dunes commonly dominate the bed of sand rivers. Due to the turbulence generation over their stoss and lee sides, dunes are of central importance in predicting hydraulic roughness and water levels. During floods in steep alluvial rivers, dunes are observed to grow rapidly as flow strength increases, undergoing an unstable transition regime, after which they are washed out in what is called upper stage plane bed. This transition of dunes to upper stage plane bed is associated with high transport of bed sediment in suspension and large decrease in bedform roughness. In the present study, we aim to improve the prediction of dune development and dune transition to upper stage plane bed by introducing the transport of suspended sediment in an existing dune evolution model. In addition, flume experiments are carried out to investigate dune development under bed load and suspended load dominated transport regimes, and to get insight in the time scales related to the transition of dunes to upper stage plane bed. Simulations with the extended model including the transport of suspended sediment show significant improvement in the prediction of equilibrium dune parameters (e.g. dune height, dune length, dune steepness, dune migration rate, dune lee side slope) both under bed load dominant and suspended load dominant transport regimes. The chosen modeling approach also allows us to model the transition of dunes to upper stage plane bed which was not possible with the original dune evolution model. The extended model predicts change in the dune shapes as was observed in the flume experiments with decreasing dune heights and dune lee slopes. Furthermore, the time scale of dune transition to upper stage plane bed was quite well predicted by the extended model. Copyright © 2015 John Wiley & Sons, Ltd.     

QUANTITATIVE ANALYSIS OF UN-MEASURED BED LOAD AT DIVERSION INTAKES OF THE YELLOW RIVER IN SHANDONG PROVINCE

1 INTRODUCTION At present, most of the instruments used for sampling sediment loads at the canal intakes along the Yellow River in Shandong Province are horizontal samplers. These kinds of samplers cannot sample the suspended load 0.3 m from the riverbed. More recently use of bottle samplers has allowed sampling of the sediment load up to 0.1 m from the riverbed. However, bottle samplers are not able to measure suspended loads with high sediment concentrations and coarse sediment particl…     

Applying regression analysis to calculating suspended sediment runoff: Specific features of the method

Major theoretical prerequisites for the application of regression analysis are considered as applied to calculating suspended sediment given observational data on water discharge and sediment load in rivers. A concept of sediment transport curve as a calculation model based on regression analysis is introduced. The role of the probabilistic structure of raw data in the determination of the optimal form of sediment transport curve and the potential structure of errors of regression models is discussed. The most widely used theoretical models of curves, as well as the methods and algorithms for assessing regression coefficients are analyzed. The application limits and the accuracy of the method are derived from its stochastic nature, and the most common methodological mistakes in its use are discussed.     

Experimental study on the effects of artificial bed roughness on turbidity currents over abrupt bed slope change

A series of experimental observations are presented in the current study to discuss the effects of artificial bed roughness on the turbidity current flowing in a rectangular channel with an abrupt change in bed slope.For this purpose,two different types of elements,sinusoidal and trapezoidal,with various heights and arrangements are considered as artificial bed roughness.A Vectrino velocity meter was used to measure the velocity and sediment concentration profiles.The effects of inlet sediment concentration on front velocity,body velocity,unit discharge,sediment concentration,and suspended load transport rate also were investigated.Accurate equations were developed for estimation of the velocity of a turbidity current over smooth and rough beds.The unexpected experimental results showed that unlike the effect of roughness height,a change in the roughness arrangement has no significant influence on the velocity of a turbidity current.Also,the effect of bed roughness on the front velocity of a denser current is more significant.     

Optimization of the Sefid-Roud Dam desiltation process using a sophisticated one-dimensional numerical model

Although water and soil conservation activities reduce reservoir sedimentation, it is inevitable that reservoirs fed by rivers transporting high amounts of sediment will experience sedimentation. The Ghezel-Ozan and Shah-Roud rivers, which flow to the Sefld-Roud reservoir dam, are both highly sediment-laden and transport significant amounts of sediment in both bed load and suspended load forms to the reservoir. Hence, it seems that the only practical way to remove the sediment from the reservoir is to flush it out using the Chasse method. In the present paper, field measurements of Chasse operation characteristics taken in previous years are presented, and a numerical model that simulates this process is introduced. After calibrating the model using field measured data, the calculated results （for reservoir pressure flushing and released sediment volume） of the numerical model were compared with other measured data for the same Chasse operation and the results agree well. Finally, using the numerical simulation results, the best approaches to ensure highly effective flushing while conserving reservoir water are presented （at least for the Sefid-Roud dam）. The operation of the bottom outlet gates, the shape of the output hydrograph, and the reservoir water level variation during flushing were optimized. In addition, the numerical model and related parameters, which need to be calibrated, are discussed.     

Interannual variation of suspended sediment load and sediment yield in an alpine catchment

Abstract

The results are described of 16 years operation of a measuring station for the automatic recording of water discharge, bed load and suspended sediment transport in the Rio Cordon catchment, a small alpine basin (5 km²) located in northeastern Italy. Hillslope erosion processes were investigated by surveying individual sediment sources repeatedly. Annual and seasonal variations of suspended sediment load during the period 1986–2001 are analysed along with their contribution to the total sediment yield. The results show that suspended load accounted for 76% of total load and that most of the suspended sediment transport occurred during two flood events: an extreme summer flash flood in September 1994 (27% of the 16-years total suspended load) and a snowmelt-induced event in May 2001 accompanied by a mud flow which fed the stream with sediments. The role of active sediment source areas is discussed in relation to the changes in flood peak—suspended load trends which became apparent after both the 1994 and the 2001 events.     

2D NUMERICAL SIMULATION OF FLOOD AND FLUVIAL PROCESS IN THE MEANDERING AND ISLAND-BRAIDED MIDDLE YANGTZE RIVER

The characteristics of water flow and sediment transport in a typical meandering and island-braided reach of the middle Yangtze River is investigated using a two-dimensional （2D） mathematical model. The major problems studied in the paper include the carrying capacity for suspended load, the incipient velocity and transport formula of non-uniform sediment, the thickness of the mixed layer on the riverbed, and the partitioning of bed load and suspended load. The model parameters are calibrated using extensive field data. Water surface profiles, distribution of flow velocities, riverbed deformation are verified with site measurements. The model is applied to a meandering and island-braided section of the Wakouzi-Majiazui reach in the middle Yangtze River, which is about 200 km downstream from the Three Gorges Dam, to study the training scheme of the navigation channels. The model predicts the processes of sediment deposition and fiver bed erosion, changes of flow stage and navigation conditions for the first 20 years of impoundment of the Three Gorges Project.     

CASE STUDIES OF THREE-DIMENSIONAL NUMERICAL SIMULATION FOR TOTAL SEDIMENT TRANSPORT

1 INTRODUCTION With advances in computer technology and numerical methods, three-dimensional (3D) mathematical models for sediment transport are gradually applied more often and for more practical projects of hydraulic engineering. Three-dimensional mathematical river models can describe not only the secondary flow, but also the transport, deposition, and erosion of sediment in the river channel and the flood plain. However, at present the theories of sediment transport are not as well de…     

Physics‐based numerical modelling of large braided rivers dominated by suspended sediment

Physics‐based models have been increasingly developed in recent years and applied to simulate the braiding process and evolution of channel units in braided rivers. However, limited attention is given to lowland braided rivers where the transport of suspended sediment plays a dominant role. In the present study, a numerical model based on the basic physics laws of hydrodynamics and sediment transport is used to simulate the evolution process of a braided river dominated by suspended load transport. The model employs a fractional method to simulate the transport of graded sediments and uses a multiple‐bed‐layer approach to represent the sediment sorting process. An idealized braided river has been produced, with the hydrodynamic, sediment transport and morphological processes being analysed. In particular, the formation process of local pool–bar units in the predicted river has been investigated. A sensitivity analysis has also been undertaken to investigate the effects of grid resolution and an upstream perturbation on the model prediction. A variety of methods are applied to analyse the geometrical and topographical properties of the modelled river. Self‐organizing characteristics related to river geometry and topography are analysed by state‐space plots, which indicate a close relationship with the periodical erosion and deposition cycles of braiding. Cross‐sectional topography and slope frequency display similar geometries to natural rivers. Scaling characteristics are found by correlation analysis of bar parameters. Copyright © 2014 John Wiley & Sons, Ltd.     

Dispersal pattern of suspended sediment in the shear frontal zone off the Huanghe (Yellow River) mouth

The in situ records of a cruise in September 1995 off the Huanghe mouth and laboratory measurements indicate that the shear front off the river mouth results from the phase difference between the nearshore and offshore tides and plays significant role in the river-laden sediment dispersal. Two types of shear front, identified from the behaviors of currents inside and outside the shear front, alternate over tidal cycle, each of which lasts for ∼2–3 h. The dispersal patterns of suspended sediment at the stations inside and outside the shear front are distinctly different from each other. In addition, the gravity-driven hyperpycnal flow generated near the mouth is terminated within shallow water due to the barrier effect of shear front. A dispersal pattern of river-laden suspended sediment in the shear frontal zone is proposed to interpret the difference of sediment transport inside and outside the shear front. The fresh and highly turbid river effluents discharge to the sea during ebb tides and are transported northwestwards inside the shear front under the combined impacts of northward ebb currents, down-slope transport of hyperpycnal flow and confining action of shear front; after partially mixing with the ambient seawater the river effluents are then transported southeastwards outside the shear front along the flood currents, causing the intermittent increase in suspended sediment concentration and corresponding decrease in salinity outside the shear front. Over annual time scale the subaqueous slope has a geomorphological response to the ephemeral shear front. Most of the river-laden sediment deposit inside the shear front with a high accumulation rate, while erosion is dominant outside the shear front due to the lack of sediment supply.     

Multiple time scales of alluvial rivers carrying suspended sediment and their implications for mathematical modeling

Flow, sediment transport and bed deformation in alluvial rivers normally exhibit multiple time scales. Enhanced knowledge of the time scales can facilitate better approaches to the understanding of the fluvial processes. Yet prior studies of the time scales are based upon the concept of sediment transport capacity at low concentrations, which however is not generally applicable. This paper presents new formulations of the time scales of fluvial flow, suspended sediment transport and bed deformation, under the framework of shallow water hydrodynamics, non-capacity sediment transport and the theory of characteristics for the hyperbolic governing equations. The time scale of bed deformation in relation to that of flow depth is demonstrated to delimit the applicability region of mathematical river models, and the time scale of suspended sediment transport relative to that of the pertinent flow information is analyzed to address if the concept of sediment transport capacity is applicable. For shallow flows with high sediment concentrations, bed deformation may considerably affect the flow and a fully coupled model is normally required. In contrast, for deep flows at low sediment concentrations, a decoupled model is mostly justified. This pilot study of the time scales delivers a new theoretical basis, on which the interaction between flow, suspended sediment transport and bed deformation can be potentially better characterized.     

Temporal variations of sediment and morphological characteristics at a large confluence accounting for the effects of floodplain submergence

The large confluence between the Yangtze River and the outflow channel of Poyang Lake is receiving attention due to its importance in flood control and ecological protection in the Yangtze River basin. There is a large floodplain along the outflow channel of Poyang Lake, which is submerged during high flow and dry during low flow. The effects of the submergence of this floodplain on sediment and morphological characteristics at this large confluence have not been known. Hence, a field investigation was done in March 2019 (relatively high flow, Survey 3) to complement the previous field studies done in August (high flow, Survey 1) and December 2018 (low flow, Survey 2) to identify the temporal variations of sediment and morphological characteristics considering the submergence of this large floodplain. The predominant sediment transport modes were wash load for Poyang Lake and confluence particles and mixed bedload/suspended load for the Yangtze River particles. The sediment transport processes were largely affected by both the secondary flows and the water density contrast between the tributaries with a lock-exchange sediment rich, denser flow moving across the inclined mixing interface in Surveys 1 and 2. The sediment flux across the mixing interface was weakened in Survey 3 when the density contrast was very small. The stagnation zone near the confluence apex had a low sediment concentration and played a role in preventing the sediment flux exchange between the two flows, and its size, and, thus, its importance as a barrier to sediment mixing were related to the submergence of the floodplain. The bed morphology with the local scour holes at the confluence was largely affected by the large-size helical cells, and this kind of effect was weakened as the secondary flows got restricted in Survey 3. The current results expand the database and knowledge on the sediment transport and morphological features of large river confluences.     

Tidally-induced sediment transport patterns in the upper Bay of Fundy: A numerical study

The Minas Basin, the eastern end of the Bay of Fundy, is well known for its high tide ranges and strong tidal currents, which can be exploited to extract electricity power. The properties of the tidally-induced sediment transport in the Minas Basin, where significant changes in tidal processes may occur due to a recently proposed tidal power project, have been studied with a three-dimensional hydrodynamic model, an empirical bed load sediment transport model and surface sediment concentrations derived from the remotely-sensed images. The hydrodynamic model was evaluated against independent observational data, which include tidal elevation, tidal current (in the full water column and bottom layer), residual current profile and tidal asymmetry indicators. The evaluation shows that the model is in good agreement with the observations.The sediment transport includes two components, bed load and suspended particulate load. The bed load is calculated using the modelled bottom shear stress and the observed grain size data. The estimated features of bed load transport roughly agree with the observed patterns of the erosion and deposition in the Minas Basin and Cobequid Bay. The transport of the suspended load is estimated using the modelled velocity fields and the surface sediment concentration derived from remote-sensing images. The comparisons between the modelled results and the limited observations illustrate that the observed directions of suspended sediment transport are basically reproduced by the model. The modelled net suspended sediment input into the Minas Basin through Minas Passage is 2.4×10⁶ m³ yr^?1, which is comparable to the observed value of 1.6×10⁶ m³ yr^?1.The variations of the bed load and the suspended load in space and time are also presented. The total net transport, defined as the mean value of the sum of bed and suspended load transports during the tidal cycle, shows strong spatial variability. The magnitude of the transport flux ranges from 0.1 to 0.2 kg m^?1 s^?1 in Minas Channel and Minas Passage, 0.1 kg m^?1 s^?1 in Cobequid Bay, to 0.01 kg m^?1 s^?1 in the central Minas Basin and Southern Bight. In Minas Channel, the sediment transport follows the structure of the tidal residual circulation, which features a large anticlockwise gyre. The sediment in Minas Passage moves eastward and deposits into the central Minas Basin. However, the sediment from the eastern part of the Basin moves westward and deposits in the central Minas Basin as well. In the Cobequid Bay, sediment moves eastward and deposits in the upper bay.     

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.     

DEPTH-INTEGRATED MODELING FOR AGGRADING/DEGRADING MOBILE CHANNELS I： MODEL DEVELOPMENT

1 INTRODUCTION The study of sediment transport in alluvial river is one of the most important fields in hydraulic engineering. Sediment transport has direct influence on the evolution of riverbeds, estuaries and coastlines, and, in turn, affects decision-making of flood control, operating rule of reservoir, design of hydraulic structure and many other aspects. Models with different orders of dimensions were presented in the literature, and most of them had common basis that they were formu…     

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.