Modeling of turbidity dynamics caused by wind-induced waves and current in the Taihu Lake

A simple turbidity model was developed with a sound physical basis based on in situ high-frequency observations of short-term, strong wind-induced sediment suspension in Taihu Lake, China. The validation results show that the model could successfully simulate turbidity caused by strong wind events, despite the relatively poor simulation accuracy for high values of turbidity caused by the entrainment of cyanobacteria by turbulence. The in situ observations and model simulation results indicate that the wind waves were within a narrow spectral band, with spectral energy mainly distributed within the 0.28–0.75 Hz band on opposite sides of the peak frequency. These high-frequency and low-energy wind waves are sensitive to depth filtering. However, the average depth of the lake is only 1.9 m, and wind waves still represent the main force of sediment suspension at the sediment-water interface. Moreover, lake currents were of significance to the maintenance of background turbidity in calm waves or ripples and in the determination of critical shear stress. By considering the spatial distribution of hydrodynamics and sediment, the model can be used to simulate the turbidity of the entire lake as well as boundary conditions for three-dimensional numerical models.     

Fully coupled mathematical modeling of turbidity currents over erodible bed

Turbidity currents may feature active sediment transport and rapid bed deformation, such as those responsible for the erosion of many submarine canyons. Yet previous mathematical models are built upon simplified governing equations and involve steady flow and weak sediment transport assumptions, which are not in complete accordance with rigorous conservation laws. It so far remains unknown if these could have considerable impacts on the evolution of turbidity currents. Here a fully coupled modeling study is presented to gain new insights into the evolution of turbidity currents. The recent analysis of the multiple time scales of subaerial sediment-laden flows over erodible bed [Cao Z, Li Y, Yue Z. Multiple time scales of alluvial rivers carrying suspended sediment and their implications for mathematical modeling. Adv Water Resour 2007;30(4):715–29] is extended to subaqueous turbidity currents to complement the fully coupled modeling. Results from numerical simulations show the ability of the present coupled model to reproduce self-accelerating turbidity currents. Comparison among the fully and partially coupled and decoupled models along with the analysis of the relative time scale of bed deformation explicitly demonstrate that fully coupled modeling is essential for refined resolution of those turbidity currents featuring active sediment transport and rapid bed deformation, and existing models based on simplified conservation laws need to be reformulated.     

Forward-modeling of co-evolution of turbidity currents,sediment transport,and cyclic steps in the Rio Muni Basin

Previous quantitative studies of field-scale cyclic steps are mostly based on analysis of field data. Such studies have shed light on the erosion/deposition patterns over these morphological features as well as the magnitudes of the turbidity current parameters back estimated using the measured geometry data. However, it remains unclear to what extent such back estimated hydraulic features and erosion/deposition patterns can be numerically reproduced by process-based numerical models. Here, a tw...     

Spatio-temporal deposition profile of an experimentally produced turbidity current with a continuous suspension supply

A turbidity current is a turbulent, particle-laden gravity current that is driven by density differences resulting from the presence of suspended sediment particles. The current travels downslope, bearing a large amount of sediment over a great distance, and forms fluvial and submarine bedforms. Knowledge of the spatio-temporal deposition profile of turbidity-deposited sediment is important for a better understanding of sediment transport by turbidity currents. In the current study, the depositi...     

Numerical modeling of the propagation and morphological changes of turbidity currents using a cost-saving strategy of solution updating

Existing layer-averaged numerical models for turbidity currents have mostly adopted the global minimum time step (GMiTS) for solution updating, which confines their computational efficiency and limits their attractiveness for field applications. This paper presents a highly efficient layer-averaged numerical model for turbidity currents by implementing the combined approach of the local graded-time-step (LGTS) and the global maximum-time-step (GMaTS). The governing equations are solved for unstructured triangular meshes by the shock-capturing finite volume method along with a set of well-balanced evaluations of the numerical flux and geometrical slope source terms. The quantitative accuracy of the model, given reasonably estimated empirical and model parameters (e.g., bed friction, water entrainment, sediment deposition and erosion coefficients), is demonstrated by comparing the numerical solutions against laboratory data of the current front positions and deposition profiles, as well as field data of the current front positions. The improved computational efficiency is demonstrated by comparing the computational cost of the present model against that of a traditional model that uses a GMiTS. For the present simulated cases, the maximum reduction of the computational cost is approximately 80% (e.g., a simulation that cost 1 h before will only require 12 min with the new model).     

Effects of bottom slope, flocculation and hindered settling on the coupled dynamics of currents and suspended sediment in highly turbid estuaries, a simple model

This study aims at gaining basic understanding about two specific phenomena that are observed in the highly turbid estuaries tidal Ouse, Yangtze and Ems, i.e. (1) the accumulation of suspended matter in the deeper parts of the estuaries and (2) the relatively high values of turbidity near the surface in the area of the turbidity maximum. A semi-analytical model is analysed to verify the hypothesis that these phenomena result from bottom slope-induced turbidity currents and from hindered settling, respectively. The model governs the dynamics of residual flow, driven by fresh water discharge, salinity gradients and turbidity gradients. It further uses the condition of morphodynamic equilibrium (no divergence of net sediment transport) to compute the residual sediment concentration. New aspects are that depth variations on flow and mixing processes, as well as flocculation and hindered settling of sediment, are explicitly accounted for. Tides act as a source of mixing and erosion of sediment only, thus processes like tidal pumping are not considered. Model results show that the estuarine turbidity maximum (ETM) shifts in the down-slope direction, compared to the case of a constant depth. Slope-induced turbidity currents, which are directed down-slope near the bottom and up-slope near the surface, are responsible for this shift, thereby confirming the first part of the hypothesis above. The down-slope shift of the ETM is reduced by currents resulting from gradients in depth-dependent mixing, which counteract turbidity currents, but which are always weaker. Including flocculation and hindered settling yields increased surface sediment concentrations in the area of the turbidity maximum, compared to the situation of a constant settling velocity, thereby supporting the second part of the hypothesis. Sensitivity experiments reveal that the conclusions are not sensitive to the values of the model parameters.     

NUMERICAL SIMULATION OF TURBIDITY CURRENT IN RESERVOIR

I. INTRODUCTIONWhen a sediment--laden flow reaches the backwater zone of a reservoir, the suddenreduction of flow velocity causes sediment particles to settle toward the river bed. Undercertain circumstsnces, it will plunge and form a layer of sediment--water mixture flowingbeneath the water surface. This flowing layer is called the turbidity current. A turbiditycurrent is relatively stable and has important impacts on reservoir sedimentation.In the case of deep reservoirs, due to temper…     

Numerical investigation of the factors influencing the vertical profiles of current,salinity, and SSC within a turbidity maximum zone

The flow-sediment interaction plays a considerable role on the vertical (internal) profiles of current,salinity and suspended sediment concentration (SSC) within a turbidity maximum zone (TMZ).Numerical modeling provides valuable insights into the complex estuarine physical processes.By combining numerical modeling with field observations,the influencing factors of fine sediment dynamics within the TMZ of Yangtze Estuary have been explored in this study.Firstly,during the neap tide,the measured data present that the current is too weak to break the density stratification,and the vertical flow structure is effectively altered.Secondly,a three-dimensional numerical model based on the Delft3D has been developed and a range of numerical sensitivity analyses were carried out to distinguish the dominant mechanisms and physical processes responsible for the phenomena observed from the measurement data.The numerical investigation highlights the following findings.(1) The vertical profile of currents within the TMZ is largely affected by saltwater intrusion,especially during lower currents when the baroclinic pressure gradient can significantly reshape the local vertical profiles of velocity.(2) The baroclinic effects are primarily determined by the stratification of salinity.(3) In addition to salinity,SSC also influences the local density stratification when its contribution to fluid density is comparable to that of salinity.(4) The settling velocity determines the overall sediment distribution and vertical profiles of the SSC in the water column.The SSC-dependent settling velocity (including the flocculation-induced acceleration and hindered settling deceleration phases) affects the longitudinal movement of the sediments.(5) The vertical profiles of current,salinity and SSC within the TMZ are highly associated with the turbulence determined by the model.The approach to modulate the vertical eddy viscosity in the model,based on the empirical dependency between Rig and Prt,may lead to a numerical instability in the stratified flow.In order to improve the stratification of SSC,additional turbulence damping effect is suggested to be implemented in the model.     

Feedback between residual circulations and sediment distribution in highly turbid estuaries: An analytical model

Motivated by field studies of the Ems estuary which show longitudinal gradients in bottom sediment concentration as high as O(0.01 kg/m⁴), we develop an analytical model for estuarine residual circulation based on currents from salinity gradients, turbidity gradients, and freshwater discharge. Salinity is assumed to be vertically well mixed, while the vertical concentration profile is assumed to result from a balance between a constant settling velocity and turbulent diffusive flux. Width and depth of the model estuary are held constant. Model results show that turbidity gradients enhance tidally averaged circulation upstream of the estuarine turbidity maximum (ETM), but significantly reduce residual circulation downstream, where salinity and turbidity gradients oppose each other. We apply the condition of morphodynamic equilibrium (vanishing sediment transport) and develop an analytical solution for the position of the turbidity maximum and the distribution of suspended sediment concentration (SSC) along a longitudinal axis. A sensitivity study shows great variability in the longitudinal distribution of suspended sediment with the applied salinity gradient and six model parameters: settling velocity, vertical mixing, horizontal dispersion, total sediment supply, fresh water flow, and water depth. Increasing depth and settling velocity move the ETM upstream, while increasing freshwater discharge and vertical mixing move the ETM downstream. Moreover, the longitudinal distribution of SSC is inherently asymmetric around the ETM, and depends on spatial variations in the residual current structure and the vertical profile of SSC.     

Numerical model of local scour considering the unsteady sediment inflow and sediment sorting: Application to scour upstream of slit weir

The current study proposes a novel framework for the numerical model for estimating the temporal scour considering unsteady sediment inflow and the sediment sorting process. The framework was applied to local scour upstream of a slit weir. The scour model is based on an ordinary nonlinear differential equation derived from sediment continuity and scour rate equations. A one-dimensional(1-D)Boussinesq-type model coupled with nonequilibrium sediment transport was incorporated in the scour model to...     

Time‐varying suspended sediment‐discharge rating curves to estimate climate impacts on fluvial sediment transport

This study presents time‐varying suspended sediment‐discharge rating curves to model suspended‐sediment concentrations (SSCs) under alternative climate scenarios. The proposed models account for hysteresis at multiple time scales, with particular attention given to systematic shifts in sediment transport following large floods (long‐term hysteresis). A series of nested formulations are tested to evaluate the elements embedded in the proposed models in a case study watershed that supplies drinking water to New York City. To maximize available data for model development, a dynamic regression model is used to estimate SSC based on denser records of turbidity, where the parameters of this regression are allowed to vary over time to account for potential changes in the turbidity‐SSC relationship. After validating the proposed rating curves, we compare simulations of SSC among a subset of models in a climate change impact assessment using an ensemble of flow simulations generated using a stochastic weather generator and hydrologic model. We also examine SSC estimates under synthetic floods generated using a peaks‐over‐threshold model. Our results indicate that estimates of extreme SSC under new climate and hydrologic scenarios can vary widely depending on the selected model and may be significantly underestimated if long‐term hysteresis is ignored when simulating impacts under sequences of large storm event. Based on the climate change scenarios explored here, average annual maximum SSC could increase by as much as 2.45 times over historical values.     

3D Modeling of sediment movement by ships-generated wakes in confined shipping channel

Ship-generated waves and return currents are capable of re-suspending significant quantities of bottom and bank sediments.However,most of the previous studies done on the subject do not show how and where sediment is re-suspended by the wakes and the directions of net transport.In this paper,a 3D numerical model based on hydro-sedimentary coupling is presented to search the relationship between the sediment movement,and the pattern of ship-generated waves around and far away from the vessel and the return currents around the ships.The hydrodynamic model is based on 3D Navier-Stokes equations including the standard k-ε model for turbulence processes,and the sediment transport model is based on a 3D equation for the re-suspended sediment transport.The computation results show that the areas of sediment concentration and transport(whether by resuspension or by the bedload) depend mainly on the position,the speed of the ship in the waterways,the kinematics of ship-generated waves and on the return flows.Thus,a map of sediment distribution and the modes of sediment transport generated by the passage of the ship are presented.     

Quasi-stationary flow structure in turbidity currents

A turbidity current is a particle-laden current driven by density differences due to suspended sediment particles. Turbidity currents can transport large amounts of sediment down slopes over great distances, and play a significant role in fluvial, lake and submarine systems. To better understand the sediment transport process, the flow system of an experimentally produced turbidity current in an inclined flume was investigated using video processing. We observed that the current progresses with constant frontal velocity and maintains an unchanged global interface geometry. In addition, the spatio-temporal profiles of the inner mean and turbulence velocity obtained by ultrasound velocity profiler (UVP) showed that similar distributions were maintained, with low dissipation. The results indicate that the turbidity current progressed in a quasi-stationary state, which enabled long-distance sediment transport. To understand the mechanisms behind the quasi-stationary flow, we analyzed the forces acting on the turbidity current. We found that under particular densities of suspended particles, the gravitational force is balanced by the viscous forces along the slope direction. We conclude that this specific force balance induces the quasi-stationary flow structure, enabling the long-distance transport of a substantial amount of sediment downstream with low dissipation.     

Numerical modeling of lock-exchange gravity/turbidity currents by a high-order upwinding combined compact difference scheme

This study presents two-dimensional direct numerical simulations for sediment-laden current with higher density propagating forward through a lighter ambient water.The incompressible NavierStokes equations including the buoyancy force for the density difference between the light and heavy fluids are solved by a finite difference scheme based on a structured mesh.The concentration transport equations are used to explore such rich transport phenomena as gravity and turbidity currents.Within the framework of an Upwinding Combined Compact finite Difference(UCCD)scheme,rigorous determination of weighting coefficients underlies the modified equation analysis and the minimization of the numerical modified wavenumber.This sixth-order UCCD scheme is implemented in a four-point grid stencil to approximate advection and diffusion terms in the concentration transport equations and the first-order derivative terms in the Navier-Stokes equations,which can greatly enhance convective stability and increase dispersive accuracy at the same time.The initial discontinuous concentration field is smoothed by solving a newly proposed Heaviside function to prevent numerical instabilities and unreasonable concentration values.A two-step projection method is then applied to obtain the velocity field.The numerical algorithm shows a satisfying ability to capture the generation,development,and dissipation of the Kelvin-Helmholz instabilities and turbulent billows at the interface between the current and the ambient fluid.The simulation results also are compared with the data in published literatures and good agreements are found to prove that the present numerical model can well reproduce the propagation,particle deposition,and mixing processes of lock-exchange gravity and turbidity currents.     

Seasonal variation of water column structure and sediment transport in a mud depo-center off the Zhejiang-Fujian coast in China

Two surveys were conducted in December, 2008, and August, 2009, in the mud depo-center off the Zhejiang-Fujian coast (MDZFC) in the inner shelf of East China Sea to depict the seasonal variation of the water column structure and analyze the factors responsible for the variation. The results were also used to discuss the sediment transport process and formation mechanism of the MDZFC. The water column structures varied significantly between the two surveys, with respect to the temperature, salinity, and turbidity. The summer water body, with relatively high temperatures and salinities, was evidently stratified with respect to the temperature, whereas the salinity remained constant throughout the water column. The stratification restricts sediment resuspension and transport. From the north to the south, the temperature in the middle-bottom water layer slightly increased, whereas the salinity remained mostly constant. In winter, the water body, with relatively low temperatures and salinities, was well mixed vertically. The temperature and salinity both increased from the surface to the bottom toward the east (deep water) and the south. A wedge-shaped water mass, which appears as a coastal upwelling, with relatively low temperature and high salinity in summer and relatively high temperature and high salinity in winter, spread landward along the sea floor, from the sea deeper than 50 m, whereas the extension was relatively stronger in winter. The water turbidity in winter was clearly higher than in summer. In the surface layer, the turbidity was generally greater than 5 FTU in winter and less than 1 FTU in summer. In the bottom layer, the turbidity was much greater than 200 FTU in winter and slightly greater than 50 FTU in summer. Moreover, the turbid water layer close to the sea floor in winter can reach into an area deeper than 50 m with a thickness of over 10 m; however, it was only limited to only 30-m-deep water with a thickness of 5 m in summer. The differences of marine sedimentary environment in the MDZFC were attributed to the seasonal variations of hydrodynamics environment, weather conditions, sediment supplies, and seasonal circulations. The results suggest that winter is the key season for particle transportation and deposition. The bottom turbid layer is the primarily channel of sediment transport, and the upwelling currents and the oceanic front systems play an important role in the sediment deposit processes and the formation of the MDZFC.     

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…     

PHYSICAL MODELING OF HYDRAULIC DESILTATION IN TAPU RESERVOIR

1 INTRODUCTIONIn recent years, the concePt of long-term sustained use of reservoirs has been addressed because areservoir is very much considered to be a nonrenewable resource (Morris and Fan, l998). Technically,many options for reservoir sedAnentation control can be utlized to pursue the sustainable develoPment ofwater resources. In general, reduction of incoming sedimen yields from watersheds is often emPloyedin conjunction with hydraulic methods such as flushing or density currnt vot…     

Varved clay paleomagnetism: Sedimentogenesis and secular variation record

The genesis, varve, and paleomagnetic informativeness of glaciolacustrine sediments are considered based on paleomagnetic and lithological studies. The magnetic properties of varved clays from periglacial basins have been considered using the sedimentological model of formation of varved clays as products of glacial turbidity currents, which explains the differences in the magnetic characteristics and degree of paleomagnetic informativeness of these clays. The distal zone sediments (formed due to precipitation from glacial turbidity currents, suspension, and sometimes along the periphery of the proximal zone) are most successfully used to study secular variations. The possible errors of the secular variation records in varved clays and such errors in estimating annual bedding of these clays have been analyzed. It has been indicated that information about the composition, concentration, and magnetic fraction size should be completed with detailed lithological studies, which make it possible to reconstruct the character of clastic material income into a basin and to decipher a paleomagnetic record. An analysis of the magnetic susceptibility anisotropy parameters and the magnetization formation in a sediment makes it possible to take into account disturbances caused by the sedimentation process and, thereby, to detect a pure geomagnetic signal     

Numerical simulation of tidal current and erosion and sedimentation in the Yangshan deep-water harbor of Shanghai

Yangshan near-shore sea area is the multi-island and multi-channel area with strong flow velocity and high suspended sediment concentration. Based on the characteristics of tidal currents, waves, and sediment in the Yangshan area, a two-dimensional numerical model of tidal currents, sediment transport, and sea bed deformation is developed. In the model, the effects of tidal currents and wind waves on sediment transport are considered. According to characteristics of the study area, unstructured grids are applied to fit the boundaries of the near-shore sea area. The results show that the calculated values are in good agreement with the measured data. The field of tidal currents, suspended sediment concentrations, and the deformation of the seabed can be successfully simulated.     

Reducing the impacts of flood-induced reservoir turbidity on a regional water supply system

This paper proposed an integrated simulation model to incorporate the impact of flood-induced reservoir turbidity into water supply. The integrated model includes a regional water allocation model and a one-dimensional settling model of cohesive particles based on Kynch’s theory. It simulates the settling of sediment flocculation in a turbid reservoir. The restrictions of water supply during floods is mimicked by simulating turbidity profiles for control points and then quantifying the associated treatment capability of raw water in the regional water allocation model for each time step. This framework can simulate shortages caused by flood-induced high turbidity as well as extended droughts, thus provide a basis for comprehensive evaluations of emergent and regular water supply facilities. A case study of evaluating different measures to mitigate the impact of turbid reservoir on water supply in northern Taiwan is presented to demonstrate the efficacy of the proposed approach.