Coupled Model of Two-phase Debris Flow,Sediment Transport and Morphological Evolution

The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.

Evaluation of an erosion-sediment transport model for a hillslope using laboratory flume data

Climate change can escalate rainfall intensity and cause further increase in sediment transport in arid lands which in turn can adversely affect water quality. Hence, there is a strong need to predict the fate of sediments in order to provide measures for sound erosion control and water quality management. The presence of micro- topography on hillslopes influences processes of runoff generation and erosion, which should be taken into account to achieve more accurate modelling results. This study presents a physically based mathematical model for erosion and sediment transport coupled to one-dimensional overland flow equations that simulate rainfall-runoff generation on the rill and interrill areas of a bare hillslope. Modelling effort at such a fine resolution considering the flow con- nection between Jnterrill areas and rills is rarely verified. The developed model was applied on a set of data gath- ered from an experimental setup where a 650 cm×136 cm erosion flume was pre-formed with a longitudinal rill and interrJll having a plane geometry and was equipped with a rainfall simulator that reproduces natural rainfall characteristics. The flume can be given both longitudinal and lateral slope directions. For calibration and validation, the model was applied on the experimental results obtained from the setup of the flume having 5% lateral and 10% longitudinal slope directions under rainfall intensities of 105 and 45 mm/h, respectively. Calibration showed that the model was able to produce good results based on the R2 （0.84） and NSE （0.80） values. The model performance was further tested through validation which also produced good statistics （R2=0.83, NSE=0.72）. Results in terms of the sedigraphs, cumulative mass curves and performance statistics suggest that the model can be a useful and an important step towards verifying and improving mathematical models of erosion and sediment transport.

Coupled 2D Hydrodynamic and Sediment Transport Modeling of Megaflood due to Glacier Dam-break in Altai Mountains,Southern Siberia

One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approaches to enhancing the understanding of the flood events.The computational domain of this flood is over 9460 km2 and about 3.784 × 106 cells are involved as a 50 m × 50 m mesh is used,which necessitates a computationally efficient model.Here the Open MP(Open Multiprocessing) technique is adopted to parallelize the code of a coupled 2D hydrodynamic and sediment transport model.It is shown that the computational efficiency is enhanced by over 80% due to the parallelization.The floods over both fixed and mobile beds are well reproduced with specified discharge hydrographs at the dam site.Qualitatively,backwater effects during the flood are resolved at the bifurcation between the Chuja and Katun rivers.Quantitatively,the computed maximum stage and thalweg are physically consistent with the field data of the bars and deposits.The effects of sediment transport and morphological evolution on the flood are considerable.Sensitivity analyses indicate that the impact of the peak discharge is significant,whilst those of the Manningroughness,medium sediment size and shape of the inlet discharge hydrograph are marginal.

Modeling the circulation and sediment transport in the Beibu Gulf

Water circulation and sediment transport in the Beibu Gulf are important for its environmental protection and resource exploitation.By employing the Regional Ocean Modeling System（ROMS）,we studied the seasonal variation of circulation,sediment transport and long-term morphological evolution in the Beibu Gulf.The simulation results show that the circulation induced by tide and wind is cyclonic both in winter and summer in the gulf and that the wind-driven circulation is stronger in winter than that in summer.The sediment concentration is higher in the Qiongzhou Strait,west of the Hainan Island and the coast of Vietnam and the Leizhou Peninsula.The sediment is transported westwards in winter and eastwards in summer in the Qiongzhou Strait.The west entrance of the Qiongzhou Strait is dominated by westward transport all the year round.The sediment discharged by rivers is deposited near the river mouths.The simulated result demonstrates that the sediment transport is mainly controlled by tidal induced bottom resuspension in the Beibu Gulf.Four characteristics are summarized for the distribution patterns of erosion and deposition.（1） The erosion and deposition are insignificant in most area of the gulf.（2） Sediment deposition is more significant in the mouths of Qiongzhou Strait.（3） The erosion is observed in the seabed of Qiongzhou Strait.（4） Erosion and deposition occur alternatively in the west of Hainan Island.

A Higher-Efficient Non-Hydrostatic Finite Volume Model for Strong Three-Dimensional Free Surface Flows and Sediment Transport

In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional （2-D） grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.

Comparison of depth-averaged concentration and bed load flux sediment transport models of dam-break flow

This paper presents numerical simulations of dam-break flow over a movable bed. Two different mathematical models were compared: a fully coupled formulation of shallow water equations with erosion and deposition terms(a depth-averaged concentration flux model), and shallow water equations with a fully coupled Exner equation(a bed load flux model). Both models were discretized using the cell-centered finite volume method, and a second-order Godunov-type scheme was used to solve the equations. The numerical flux was calculated using a Harten, Lax, and van Leer approximate Riemann solver with the contact wave restored(HLLC). A novel slope source term treatment that considers the density change was introduced to the depth-averaged concentration flux model to obtain higher-order accuracy. A source term that accounts for the sediment flux was added to the bed load flux model to reflect the influence of sediment movement on the momentum of the water. In a onedimensional test case, a sensitivity study on different model parameters was carried out. For the depth-averaged concentration flux model,Manning's coefficient and sediment porosity values showed an almost linear relationship with the bottom change, and for the bed load flux model, the sediment porosity was identified as the most sensitive parameter. The capabilities and limitations of both model concepts are demonstrated in a benchmark experimental test case dealing with dam-break flow over variable bed topography.

Dynamic simulation and modeling of PCP transport between sediment and water

A model to calculate the concentration of organic pollutant in the water column is presented. The inflow, outflow, adsorption of pollutant to the suspended particles and settling of particles, diffusion across the sediment water interface, and volatilization to the atmosphere were considered in this model. Dynamic simulation experiment of pentachlorophenol(PCP) has been made in the laboratory. The model parameters were also estimated with experimental data(diffusion mass transfer coefficient, 0.00112m/d; net settling velocity of particles, 0.323 m.d -1 ; volatilization rate constant, 1.94×10 -4 d -1 ). There are little differences between calculated values and measured values. Adsorption of pollutant to the particles and settling of particles are very important mechanisms for removal of organic pollutant from the water. The principles and methods of this model are also applicable to the prediction of the concentration of other trace elements and organic constituents in aquatic systems.

Linking Fish Habitat Modelling and Sediment Transport in Running Waters

The assessment of ecological status for running waters is one of the major issues within an integrated river basin management and plays a key role with respect to the implementation of the European Water Frame- work Directive (WFD).One of the tools supporting the development of sustainable river management is physi- cal habitat modeling,e.g.,for fish,because fish population are one of the most important indicators for the e- colngical integrity of rivers.Within physical habitat models hydromorphological ...

Applicative limitations of sediment transport on predictive modeling in geomorphology

Sources of uncertainty or error that arise in attempting to scale up the results of laboratory-scale sediment transport studies for predictive modeling of geomorphic systems include: (i) model imperfection, (ii) omission of important processes, (iii) lack of knowledge of initial conditions, (iv) sensitivity to initial conditions, (v) unresolved heterogeneity, (vi) occurrence of external forcing, and (vii) inapplicability of the factor of safety concept. Sources of uncertainty that are unimportant or that can be controlled at small scales and over short time scales become important in large-scale applications and over long time scales. Control and repeatability, hallmarks of laboratory-scale experiments, are usually lacking at the large scales characteristic of geomorphology. Heterogeneity is an important concomitant of size, and tends to make large systems unique. Uniqueness implies that prediction cannot be based upon first-principles quantitative modeling alone, but must be a function of system history as well. Periodic data collection, feedback, and model updating are essential where site-specific prediction is required.

A three-dimensional,wave-current coupled,sediment transport model for POM

In the high-energy environment of coastal seas and estuaries,strong sediment resuspension/ deposition events are driven by surface waves,tides,winds and buoyancy driven currents.In recent years,A POM based three-dimensional,wave-current coupled,sediment transport model has been developed by the University of New South Wales.This paper presents several examples of the model applications to study sediment dynamics in the environments where forcings such as waves,tides,and winds are equally important to affect sediment fluxes and distributions.Firstly,the sediment transport model coupled to the Yellow Sea general circulation model and a third generation wave model SWAN was implemented in the Yellow Sea to study the dynamics of the sediment transport and resuspension in the northern Jiangsu shoal-wate(rNJSW).The sediment distributions and fluxes and their inter-annual variability were studied by realistic numerical simulations.The study found that the surface waves played a dominant role over the tides to form the turbidity maxima along the muddy coast of NJSW. Secondly,the sediment transport model was used to explore the effect of suspended sediment-induced stratification in the bottom boundary laye(rBBL).The model uses a re-parameterized bottom drag coefficient Cd that incorporates a linear stability function of flux Richardson number Rf.The study has shown that the sediment induced stratification in the BBL reduces the vertical eddy viscosity and bottom shear stress in comparison with the model prediction in a neutrally stratified BBL.In response to these apparent reductions,the tidal current shear is increased and sediments are abnormally concentrated within a thin wall layer that is overlain by a thicker layer with much smaller concentration.The formation of this fluid-mud layer near the seabed has led to a significant reduction in the total sediment transport.This study contributes to the understanding of formations of tidal flats along the coasts of turbid seas and estuaries.

Developing a 2D vertical flow and sediment transport model for open channels using the Youngs-VOF method

A 2D vertical （2DV） numerical model, without o-coordinate transformation in the vertical direction, is developed for the simulation of flow and sediment transport in open channels. In the model, time-averaged Reynolds equations are closed by the k-e nonlinear turbulence model. The modified Youngs- VOF method is introduced to capture free surface dynamics, and the free surface slope is simulated using the ELVIRA method. Based on the power-law scheme, the k-e model and the suspended-load transport model are solved numerically with an implicit scheme applied in the vertical plane and an explicit scheme applied in the horizontal plane. Bedload transport is modeled using the Euler-WENO scheme, and the grid-closing skill is adopted to deal with the moving channel bed boundary. Verification of the model using laboratory data shows that the model is able to adequately simulate flow and sediment transport in open channels, and is a good starting point for the study of sediment transport dynamics in strong nonlinear flow scenarios.

On the Existence of Global Weak Solutions to 1D Sediment Transport Model

This paper is devoted to the study of the existence of weak solution in time with a periodic domain of sediment transport model. We consider a one-dimensional viscous sediment transport model which combines a viscous Shallow-Water system with a transport equation that describes the bottom evolution. The model studied does not take into account all the regularizing terms used by Roamba Brahima, Zongo Yacouba and Jean de Dieu Zabsonré (2017) and we use a better transport equation than that used by Zabsonré (2012).

Numerical Modeling of Sediment Transport and Its Effect on Algal Biomass Distribution in Lake Pontchartrain Due to Flood Release from Bonnet CarréSpillway

In order to protect the city of New Orleans from the Mississippi River flooding, the Bonnet Carré Spillway (BCS) was constructed from 1929 to 1936 to divert flood water from the river into Lake Pontchartrain and then into the Gulf of Mexico. During the BCS opening for flood release, large amounts of freshwater, nutrients, sediment, etc. were discharged into Lake Pontchartrain, and caused a lot of environmental problems. To evaluate the environmental impacts of the flood water on lake ecosystems, a two-dimensional numerical model was developed based on CCHE2D and applied to simulate the flow circulation, sediment transport and algal biomass distribution in Lake Pontchartrain. The effect of sediment concentration on the growth of algae was considered in the model. The numerical model was calibrated using field measured data provided by USGS, and then it was validated by the BCS Opening Event in 1997. The simulated results were generally in good agreement with filed data and satellite imagery. The field observation and numerical model show that during the spillway opening for flood release, the sediment concentration is very high, which greatly restricts the growth of algae, so there is no algal bloom observed in the lake. After the closure of BCS, the sediment concentration in the lake reduces gradually, and the nutrient concentration of the lake is still high. Under these conditions, numerical results and satellite imagery showed that the chlorophyll concentration was high and algal bloom might occur.

Modeling of suspended sediment by coupled wave-current model in the Zhujiang(Pearl) River Estuary

A three-dimensional wave-current-sediment coupled numerical model is developed to understand the sediment transport dynamics in the Zhujiang(Pearl)River Estuary(ZRE),China.The model results are in good agreement with observed data,and statistics show good model skill scores.Numerical studies are conducted to assess the scenarios of suspended sediment in the ZRE under the effects of different forcing(river discharges,waves,and winds).The model results indicate that the estuarine gravitational circulation plays an important role in the development of estuarine turbidity maximum in the ZRE,particularly during neap tides.The increased river discharge can result in a seaward sediment transport.The suspended sediment concentration(SSC)in the bottom increases with both wave bottom orbital velocity and wave height.Because of the shallow water depth,the effect of waves on sediment in the west shoal is greater than that in the east channel.The southwesterly wind-induced wave affects the SSC more than those resulting from the northeasterly wind,while the northeasterly wind-driven circulation has a slightly greater influence on the SSC than that of the southwesterly wind.However,a steady southwesterly wind condition favors the increase of the SSC in the Lingding Bay more so than a steady northeasterly wind condition.If the other forcings are same,the averaged SSC under a steady southwesterly wind condition is about 1.1 times that resulting from a steady northeasterly wind.

The reverse sediment transport trend between abandoned Huanghe River(Yellow River) Delta and radial sand ridges along Jiangsu coastline of China——an evidence from grain size analysis

To reveal the sediment transporting mechanism between the abandoned Huanghe River （Yellow River） Delta and radial sand ridges, “End Member” Model and grain size trend analysis have been employed to separate the “dynamic populations” in the surficial sediment particle spectra and to determine the possible sediment transporting pathway. The results reveal four “dynamic subpopulations”（EM1 to EM4） and two reverse sediment transporting directions： a northward transport tend from the radial sand ridges to mud patch, and a southward transport trend in deep water area outside the mud patch. Combined with the published hydrodynamic information, the transporting mechanism of dynamic populations has been discussed, and the main conclusion is that the transporting of finer subpopulations EM1 and EM2 is controlled by the “anticlockwise residual current circulation” forming during tidal cycle, which favor a northward transporting trend and the forming of mud patch on the north of radial sand ridges, while the transporting of coarser EM3 is mainly controlled by wind driven drift in winter, which favors a southward transporting direction.

Numerical Modelling Sediment-Bacteria Interaction Processes in the Severn Estuary

Faecal bacteria exist in both free-living and attached forms in estuarine waters. The deposition of sediments can take faecal bacteria out of the water column and to the bed. The sediments can subsequently be re-suspended to the water column, which can then lead to re-suspension of the faecal bacteria of the attached forms back to the water column. Therefore, the fate and transport of faecal bacteria is highly related to the governing sediment transport processes, particularly where these processes are significant, such as the Severn Estuary, UK. However, little attempt has been made to model such processes in terms of predicting the impact of the sediment fluxes on faecal bacteria levels. Details are given of the refinement of a numerical model of faecal bacteria transport, where the sediment transport processes are significant. After testing the sediment-bacteria interaction model favourably against known results in previous study, the model was applied to the Severn Estuary and Bristol Channel, UK, to investigate the impact of suspended sediment fluxes on the corresponding faecal bacteria transport processes. The model predictions have proved to be encouraging, with the results being compared to a traditional faecal bacteria modelling approach, where sediment bacteria interactions were not included. The new model provides improved predictions of faecal bacteria concentrations when sediment transport is included and for the Bristol Channel Severn Estuary it can be seen that the effects of the sediments on the bacterial levels in the water column can be significant.

Numerical simulation of sediment transport in coastal waves and wave-induced currents

Prediction of coastal sediment transport is of particularly importance for analyzing coast erosion accurately and solving the corresponding coast protection engineering problems.The present study provided a numerical scheme for sediment transport in coastal waves and wave-induced currents.In the scheme,the sand transport model was implemented with wave refraction-diffraction model and near-shore current model.Coastal water wave was simulated by using the parabolic mild-slope equation in which wave refraction,diffraction and breaking effects are considered.Wave-induced current was simulated by using the nonlinear shallow water equations in which wave provides radiation stresses for driving current.Then,sediment transport in waves and wave-induced currents was simulated by using the two-dimensional suspended sediment transport equations for suspended sediment and the bed-load transport equation for bed load.The numerical scheme was validated by experiment results from the Large-scale Sediment Transport Facility at the US Army Corps of Engineer Research and Development Center in Vicksburg.The numerical results showed that the present scheme is an effective tool for modeling coastal sediment transport in waves and near-shore currents.

Applications of GSTARS Computer Models for Solving River and Reservoir Sedimentation Problems

GSTARS (Generalized Sediment Transport model for Alluvial River Simulation) is a series of computer models developed by the U.S. Bureau of Reclamation while the author was employed by that agency. The stream tube concept is used in all GSTARS models which allow us to solve one-dimensional equations for each stream tube independently and obtain semi-two-dimensional variation of the hydraulic conditions along and across stream tubes for rivers and reservoirs. Sedi-ment transport, scour, and deposition processes are simulated along each stream tube independently to give us a semi-three-dimensional variation of the bed geometry. Most sediment transport computer models assume that channel width is given and cannot change during the simulation process. GSTARS models apply the theory of minimum stream power to the determination of optimum channel width and channel geometry. The concepts of channel side stability, and active, inactive, and armoring layers are used in all GSTARS models for realistic long-term simulation and prediction of the scour and deposition processes in rivers and reservoirs.GSTARS models have been applied in many countries for solving a wide range of river and reservoir sedimentation problems. Case studies will be used to illustrate the applications of GSTARS computer models.

Numerical Analysis of Adaptation-to-capacity Length for Fluvial Sediment Transport

Over the last several decades,various sediment transport capacity formulations have been used by geomorphologists and engineers to calculate fluvial morphological changes.However,it remains poorly understood if the adaptation to capacity could be fulfilled instantly in response to differing inflow discharges and sediment supplies,and thus if the calculation of morphological changes in rivers based on the assumed capacity status is fully justified.Here we present a numerical investigation on this issue.The distance required for sediment transport to adapt to capacity(i.e.,adaptation-to-capacity length) of both bed load and suspended sediment transport is computationally studied using a coupled shallow water hydrodynamic model,in line with varied inlet sediment concentrations.It is found that the adaptation-to-capacity length generally decreases as the Rouse number increases,irrespective of whether the inlet sediment concentration increases or reduces.For cases with vanishing inlet sediment concentration a unified relationship is found between the adaptation-to-capacity length and the Rouse number.Quantitatively,the adaptation-to-capacity length of bed load sediment is limited to tens of times of the flow depth,whilst that of suspended sediment increases substantially with decreasing Rouse number and can be up to hundreds of times of the flow depth.The present finding concurs that bed load sediment transport can adapt to capacity much more rapidly than suspended sediment transport,and it facilitates a quantitative criterion on which the applicability of bed load or suspended sediment transport capacity for natural rivers can be readily assessed.

Numerical simulations of flow and sediment transport within the Ning-Meng reach of the Yellow River,northern China

Effective management of a river reach requires a sound understanding of flow and sediment transport generated by varying natural and artificial runoff conditions. Flow and sediment transport within the Ning-Meng reach of the Yellow River（NMRYR）, northern China are controlled by a complex set of factors/processes, mainly including four sets of factors：（1） aeolian sediments from deserts bordering the main stream;（2） inflow of water and sediment from numerous tributaries;（3） impoundment of water by reservoir/hydro-junction; and（4） complex diversion and return of irrigation water. In this study, the 1-D flow ＆ sediment transport model developed by the Yellow River Institute of Hydraulic Research was used to simulate the flow and sediment transport within the NMRYR from 2001 to 2012. All four sets of factors that primarily control the flow and sediment transport mentioned above were considered in this model. Compared to the measured data collected from the hydrological stations along the NMRYR, the simulated flow and sediment transport values were generally acceptable, with relative mean deviation between measured and simulated values of 〈15%. However, simulated sediment concentration and siltation values within two sub-reaches（i.e., Qingtongxia Reservoir to Bayan Gol Hydrological Station and Bayan Gol Hydrological Station to Toudaoguai Hydrological Station） for some periods exhibited relatively large errors（the relative mean deviations between measured and simulated values of 18% and 25%, respectively）. These errors are presumably related to the inability to accurately determine the quantity of aeolian sediment influx to the river reach and the inflow of water from the ten ephemeral tributaries. This study may provide some valuable insights into the numerical simulations of flow and sediment transport in large watersheds and also provide a useful model for the effective management of the NMRYR.