Modelling fluxes of water and sediment between Venice Lagoon and the sea

To describe the exchange of water and sediment through the Venice Lagoon inlets a 3-D hydrodynamic and sediment transport model has been developed and applied to a domain comprising Venice Lagoon and a part of the Adriatic Sea. The model has been validated for both current velocities and suspended particle concentration against direct observations and from observations empirically derived fluxes from upward-looking acoustic Doppler current profiler probes installed inside each inlet. The model provides estimates of the suspended sediment transport in the lower 3 m of the water column that is not detected by acoustic Doppler current profiler sensors. The bedload model prediction has been validated against measured sand transport rates collected by sand traps deployed in the Lido and Chioggia inlets. Results indicate that, in the Lido inlet, 87% of the total load is in suspension, while the rest moves as bedload.     

"THREE DIMENSIONAL MATHEMATICAL MODEL OF SEDIMENT TRANSPORT IN ESTUARINE REGIONS-- A CASE STUDY OF THE HAIHE RIVER MOUTH" Yuchuan Bai,Huanting SHEN and Shixiong HU

The paper presents the 3D finite element simulation of tidal flow and Sediment transport in the estuarine region of the Haihe river. The proposed model adopts sigma-transformation of the hydrodynamic and sediment transport equations. The hydrodynamic and sediment transport models are verified in case of a simple test problem for which analytical solutions are available. Finally the models are applied to muddy Haihe river estuary of North China and it is claimed that hydrodynamic and sediment transport models give a reliable comparison with the observed field data. However, there are certain discrepancies, and some reasonable questions regarding the present state-of-art, in the modeling of three-dimensional multilevel hydrodynamics and sediment transport, which are provided below for answer.     

Calibration of an estuarine sediment transport model to sediment fluxes as an intermediate step for simulation of geomorphic evolution

Modeling geomorphic evolution in estuaries is necessary to model the fate of legacy contaminants in the bed sediment and the effect of climate change, watershed alterations, sea level rise, construction projects, and restoration efforts. Coupled hydrodynamic and sediment transport models used for this purpose typically are calibrated to water level, currents, and/or suspended-sediment concentrations. However, small errors in these tidal-timescale models can accumulate to cause major errors in geomorphic evolution, which may not be obvious. Here we present an intermediate step towards simulating decadal-timescale geomorphic change: calibration to estimated sediment fluxes (mass/time) at two cross-sections within an estuary. Accurate representation of sediment fluxes gives confidence in representation of sediment supply to and from the estuary during those periods. Several years of sediment flux data are available for the landward and seaward boundaries of Suisun Bay, California, the landward-most embayment of San Francisco Bay. Sediment flux observations suggest that episodic freshwater flows export sediment from Suisun Bay, while gravitational circulation during the dry season imports sediment from seaward sources. The Regional Oceanic Modeling System (ROMS), a three-dimensional coupled hydrodynamic/sediment transport model, was adapted for Suisun Bay, for the purposes of hindcasting 19th and 20th century bathymetric change, and simulating geomorphic response to sea level rise and climatic variability in the 21st century. The sediment transport parameters were calibrated using the sediment flux data from 1997 (a relatively wet year) and 2004 (a relatively dry year). The remaining years of data (1998, 2002, 2003) were used for validation. The model represents the inter-annual and annual sediment flux variability, while net sediment import/export is accurately modeled for three of the five years. The use of sediment flux data for calibrating an estuarine geomorphic model guarantees that modeled geomorphic evolution will not exceed the actual supply of sediment from the watershed and seaward sources during the calibration period. Decadal trends in sediment supply (and therefore fluxes) can accumulate to alter decadal geomorphic change. Therefore, simulations of future geomorphic evolution are bolstered by this intermediate calibration step.     

Issues in Eulerian–Lagrangian modeling of sediment transport under saltation regime

The saltation regime is very important for understanding the sediment transport mechanism. However,there is no consensus on a model for the saltation regime. This study answers several questions raised with respect to the Eulerian-Lagrangian modeling of sediment transport. The first question is why the previous saltation models that use different combinations of hydrodynamic forces yielded acceptable results? The second question is which shear lift model(i.e. a shear lift expression and its coefficient) is more appropriate? Another important question is which hydrodynamic forces have greater contributions to the saltation characteristics of a sediment particle? The last question is what are the contributions of the turbulence fluctuations as well as effects of using two-and three-dimensional(2 D and 3 D) models on the simulation results? In order to fairly answer these questions, a systematic study was done by considering different scenarios. The current study is the first attempt to clearly discuss these issues. A comprehensive 3 D saltation model for non-cohesive sediment was developed that includes all the hydrodynamic forces acting on the particle. The random nature of sediment transport was included using turbulent flow and bed-particle collision models. The eddy interaction model was applied to generate a3 D turbulent flow field. Bed-particle collisions were considered using the concept of a contact zone and a corresponding contact point. The validation of the model was done using the available experimental data for a wide range of sediment size(0.03 to 4.8 cm). For the first question, the results indicated that some of the hydrodynamic effects show opposing trends and some have negligible effects. With these opposing effects it is possible to adjust the coefficients of different models to achieve acceptable agreement with the same experimental data while omitting some aspects of the physics of the process. A suitable model for the shear lift force was developed by linking the lift coefficient to the drag coefficient and the contributions of the hydrodynamic forces and turbulence fluctuations as well as the consequences of using of 2 D and 3 D models were studied. The results indicate that the shear lift force and turbulent flow fluctuations are important factors for the saltation of both sand and gravel, and they cannot be ignored.     

Water balance prediction in stormwater infiltration basins using 2-D modeling: An application to evaluate the clogging process

Sustainable urban drainage systems are built along roads and in urban areas to collect urban runoff and avoid flooding, and to filter water pollutants. Sediment collected by runoff is deposited in the stormwater basin and progressively reduces water infiltration efficiency, leading to the clogging of the basin. To help stormwater basin managers and stakeholders better understand and predict clogging rates in order to elaborate maintenance plans and schedules, water transport prediction models are necessary. However,because of the heterogeneous sediment hydrodynamic properties inside the stormwater basin, a twodimensional(2-D) water flow model is required to predict water levels and possible overflow as accurately as possible. Saturated hydraulic conductivity(Ks) and sediment water retention curves were measured in the overall sediment layer of the stormwater basin, in addition to sediment layer thickness and organic matter content(11 sampling points). Sediment depth was used to predict organic matter(OM) content, and the OM was used to predict Ks. Water height in the basin was modeled with the HYDRUS-2 D model by taking into account the sediment hydrodynamic properties distribution. The HYDRUS-2 D model gave a satisfactory representation of the measured data. Scenarios of the hydraulic properties of stormwater basin sediment were tested over time, and hydraulic resistance, R, was calculated to assess the stormwater basin performance. Presently, after 20 years of functioning, the stormwater basin still ensures efficient water infiltration, but the first outflow(Hydraulic resistance,R 24 h)) is expected to appear in the next 5 years, and clogging(R 47 h) in the next 13 years. This 2-D water balance model makes it possible to integrate the hydrodynamic heterogeneity of a stormwater basin. It gives interesting perspectives to better predict 2-D/3-D contaminant transport.     

Numerical modeling of tidal currents, sediment transport and morphological evolution in Hangzhou Bay, China

A 2D depth-averaged numerical model is set up to simulate the macro-scale hydrodynamic characteristics, sediment transport patterns and morphological evolution in Hangzhou Bay, a large macro-tidal estuary on the eastern coast of China. By incorporating the shallow water equations, the suspended sediment transport equation and the mass-balance equation for sediment; short-term hydrodynamics, sediment transport and long-term morphological evolution for Hangzhou Bay are simulated and the underlying physical mechanisms are analyzed. The model reproduces the spatial distribution patterns of suspended sediment concentration (SSC) in Hangzhou Bay, characterized by three high SSC zones and two low SSC zones. It also correctly simulates the residual flow, the residual sediment transport and the sediment accumulation patterns in Hangzhou Bay. The model results are in agreement with previous studies based on field measurements. The residual flow and the residual sediment transport are landwards directed in the northern part of the bay and seawards directed in the southern part. Sediment accumulation takes place in most areas of the bay. Harmonic analysis revealed that the tide is flood-dominant in the northern part of the bay and ebb-dominant in the southern part of the bay. The strength of the flood-dominance increases landwards along the northern Hangzhou Bay. In turn sediment transport in Hangzhou Bay is controlled by this tidal asymmetry pattern. In addition, the direction of tidal propagation in the East China Sea, the presence of the archipelago in the southeast and the funnel-shaped geometry of the bay, play important roles for the patterns of sediment transport and sediment accumulation respectively.     

STUDY ON THE MATHEMATICAL MODEL OF HEAVY－METAL POLLUTANT TRANSPOPT－TRANSFORMATION IN RIVERS

ＳＴＵＤＹＯＮＴＨＥＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＯＦＨＥＡＶＹ－ＭＥＴＡＬＰＯＬＬＵＴＡＮＴＴＲＡＮＳＰＯＰＴ－ＴＲＡＮＳＦＯＲＭＡＴＩＯＮＩＮＲＩＶＥＲＳ￥Ｓｕｉ－ｌｉａｎｇＨｕａｎｇａｎｄＺｈａｏ－ｈｕｉＷａｎ（ＰｏｓｔｄｏｃｔｏｒａｌＦｅ...     

NUMERICAL SIMULATION OF SEDIMENT RELEASE FROM RESERVOIRS

1 INTRODUCTION Reservoirs sedimentation is a serious problem in many countries, including the I. R. of Iran. Accumulation of sediment deposits decreases worldwide reservoir storage capacity by one percent per year (Mahmood, 1987). The loss of reservoir st…     

Modeling Sediment Transport Dynamics in Thompson Island Pool,Upper Hudson River

Two-dimensional, vertically-averaged hydrodynamic and sediment transport models were developed and applied as part of a PCB fate and transport modeling study of Thompson Island Pool (TIP), Upper Hudson River. Mechanistic formulations were used to simulate cohesive and non-cohesive suspended load transport; site-specific data were extensively used to determine model inputs. This modeling approach is compared and contrasted to non-mechanistic solids transport sub-models used in other contaminant fate studies. A minimum number of model parameters were adjusted to calibrate the sediment transport model using data collected during the 1994 spring flood. The model was validated during the 1997 spring flood and for a 22-year (1977–1998) period. Successful calibration and validation of the model showed that: (1) deposition and resuspension processes were realistically and accurately formulated in the model; (2) the model is an effective diagnostic tool for quantitatively evaluating net deposition and erosion from various areas of TIP; and (3) sediment transport results can be coupled with a PCB fate model with a high degree of confidence.     

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.     

Storm-driven sediment transport in Massachusetts Bay

Massachusetts Bay is a semi-enclosed embayment in the western Gulf of Maine about 50 km wide and 100 km long. Bottom sediment resuspension is controlled predominately by storm-induced surface waves and transport by the tidal- and wind-driven circulation. Because the Bay is open to the northeast, winds from the northeast (‘Northeasters’) generate the largest surface waves and are thus the most effective in resuspending sediments. The three-dimensional oceanographic circulation model Regional Ocean Modeling System (ROMS) is used to explore the resuspension, transport, and deposition of sediment caused by Northeasters. The model transports multiple sediment classes and tracks the evolution of a multilevel sediment bed. The surficial sediment characteristics of the bed are coupled to one of several bottom-boundary layer modules that calculate enhanced bottom roughness due to wave–current interaction. The wave field is calculated from the model Simulating WAves Nearshore (SWAN). Two idealized simulations were carried out to explore the effects of Northeasters on the transport and fate of sediments. In one simulation, an initially spatially uniform bed of mixed sediments exposed to a series of Northeasters evolved to a pattern similar to the existing surficial sediment distribution. A second set of simulations explored sediment-transport pathways caused by storms with winds from the northeast quadrant by simulating release of sediment at selected locations. Storms with winds from the north cause transport southward along the western shore of Massachusetts Bay, while storms with winds from the east and southeast drive northerly nearshore flow. The simulations show that Northeasters can effectively transport sediments from Boston Harbor and the area offshore of the harbor to the southeast into Cape Cod Bay and offshore into Stellwagen Basin. This transport pattern is consistent with Boston Harbor as the source of silver found in the surficial sediments of Cape Cod Bay and Stellwagen Basin.     

Muskingum equation based downstream sediment flow simulation models for a river system

Applications of sediment transport and water flow characteristics based sediment transport simulation models for a river system are presented in this study. An existing water–sediment model and a new sediment–water model are used to formulate the simulation models representing water and sediment movement in a river system. The sediment–water model parameters account for water flow characteristics embodying sediment transport properties of a section. The models are revised formulations of the multiple water inflows model describing water movement through a river system as given by the Muskingum principle. The models are applied to a river system in Mississippi River basin to estimate downstream sediment concentration, sediment discharge, and water discharge. River system and the river section parameters are estimated using a revised and the original multiple water inflows models by applying the genetic algorithm. The models estimate downstream sediment transport rates on the basis of upstream sediment/water flow rates to a system. Model performance is evaluated by using standard statistical criteria;downstream water discharge resulting from the original multiple water inflows model using the estimated river system parameters indicate that the revised models satisfactorily describe water movement through a river system. Results obtained in the study demonstrate the applicability of the sediment transport and water flow characteristics-based simulation models in predicting downstream sediment transport and water flow rates in a river system.     

APPLICATION OF BED EVOLUTION MODELS OVER LOOSE AND RIGID AREAS

LINTRODUCTIONThesimulationofmorphodynamicsinestUariesandcoastalareasisbecauseofitSeconomicalmeritageononehandandofitSscientificcomplexityontheOtheranoutStalldingfieldofresearchanddeVelopment.TheFederalWaterwaysEngineeringandResearchinstitote(BAW)isthecentraladvisoryinstitutefortheGermanwaterwayauthority.ItisthereforereSPonsibletoprovidethescientificsuPPOrtforoptimalsedimentmanagemelltproceduresforWatCrwaymaill~ce.ThehydredynndcSystemanalysisisthemainaPProachtosupportengineeringd…     

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.     

Distribution of clay minerals in marine sediments off Chennai,Bay of Bengal,India： Indicators of sediment sources and transport processes

Clay mineralogy, texture size and statistical analyses were carried out on surface sediments from the continental shelf of Chennai, Bay of Bengal, India. The purpose of this study is to characterize the clay mineral distribution and its relation to the hydrodynamics off Chennai to identify the sources and transport pathways of the marine sediments. Characterization of clay minerals in coastal sediments by Fourier Transform Infrared （FTIR） spectroscopy has provided the association of quartz, feldspar, kaolinite, chlorite, illite and iron oxides （magnetite and hematite） derived from river catchments and coastal erosion. Kaolinite, chlorite, illite, iron oxides, and organic matter are the dominant minerals in Cooum, and Adayar region. High quartz and feldspar zones were identified in Marina, which are being confined the sand zone and paralleling the coast. The strong relationships among the wave energy density, sand, quartz and carbonate revealed that wave induced littoral drift system play a dominant role in transportation and deposition of sediments in the Chennai coast. The sediment texture and minerals data are in agreement well with the previous results of hydrodynamics and littoral drift models in this region. Multivariate statistical analyses （correlation, cluster and factor analyses） were carried out and obtained results suggested that clay minerals and organic matter are trapped in silt and clay particles, whereas quartz, feldspar and carbonate are associated with sand particles. Results of sediment sources and transport processes from this study will be useful to predict the fate of the pollutants released from land or the potential change in sediment delivery to coastal areas.     

Impact of offsite sediment transport and toxicity on remediation of a contaminated estuarine bay

Sediment of Ostrich Bay, an arm of Dyes Inlet on Puget Sound, was historically contaminated with ordnance compounds from an onshore US Navy facility. An initial recommendation for a sediment cover to mitigate benthic risks was followed by studies of sediment transport and deposition to determine whether contaminated sediment from Dyes Inlet or other offsite sources in Puget Sound may contribute to Ostrich Bay impacts. A Sediment Trend Analysis (STA) identified net sediment transport pathways throughout the bay and inlet by examining changes in grain size distributions in multiple adjacent samples. Results indicated that fine-grained sedimentary material transports into and deposits throughout the Dyes Inlet system, with no erosion or transport out of Ostrich Bay. Echinoderm larvae mortality bioassay results were elevated in fine-grained sediments of both Ostrich Bay and Dyes Inlet. Ordnance compounds were undetected, and although sediment mercury concentrations were elevated at 0.48-1.4 mg/kg in both waterbodies, the relationship with toxicity was weak. Results of the studies and sedimentation modeling indicate that impacted sedimentary material deposits throughout the Dyes Inlet/Ostrich Bay system from unknown sources and will prevent natural recovery of Ostrich Bay as well as negate long-term effectiveness of active remedial measures. Stakeholders have recognized that remediation of the bay can be achieved only after the toxicity of depositing sediment decreases.     

Effect of self-weight consolidation on a hydro-sedimentological model for the Río de la Plata estuary

The effect of the consolidation process on the morphodynamics and fine sediment dynamics of the Río de la Plata estuary is explored through a circulation-wave-sediment transport model. The consolidation model is calibrated based on settling column experimental data. Different simulations are done in order to initialize the mud layer distribution and to investigate the impact of different erosion parameter assumptions on the modeled sediment dynamics. Finally a two-year simulation is done with and without the consolidation process and realistic hydrodynamic forcings. Considering the consolidation process, the model correctly reproduces measured vertical density profiles in the Montevideo Bay access channel. The simulated suspended sediment dynamics behavior in Montevideo Bay with the consolidation process provides a more realistic deposition pattern in regard to the dredging activities.     

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.     

Dispersion dependence on retardation in a real fracture geometry using lattice-gas cellular automaton

Fractures have been recently identified in potential host rock for high level nuclear waste disposal, like indurated argilite formations. These fractures appear as potential rapid pathways for radionuclides transport and hydrodynamic properties of the transport inside these systems must thus be characterized. Miscible non-sorbing and sorbing tracers displacements were performed on a 2-D model derived from a real fracture geometry observed in the Tournemire argilite formation with a lattice-gas cellular automaton (LGA). LGA was shown to easily handle the complex geometry of such a fracture. The numerical breakthrough curves obtained were inverted with the 1-D CDE and MIM transport models. Two main conclusions were drawn: (i) at the length scale of the study, the non-sorbing tracer transport in our fracture geometry was more accurately interpreted in terms of the MIM model rather than in terms of the classical CDE model; (ii) in order to correctly model the sorbing tracers migration, the hydrodynamic dispersion coefficient value was found to increase with the increase of the retardation factor. A semi-empirical relation based on the Taylor–Aris theory was then used to describe this dependency.