Three-dimensional numerical modeling of cohesive sediment transport and wind wave impact in a shallow oxbow lake

It was observed that in some closed inland lakes sediment transport was dominated by wind-induced currents, and the sediment resuspension was primarily driven by wind-induced waves. This paper presents the development and application of a three-dimensional numerical model for simulating cohesive sediment transport in water bodies where wind-induced currents and waves are important. In the model, the bottom shear stresses induced by currents and waves were calculated, and the processes of resuspension (erosion), deposition, settling, etc. were considered. This model was first verified by a simple test case consisting of the movement of a non-conservative tracer in a prismatic channel with uniform flow, and the model output agreed well with the analytical solution. Then it was applied to Deep Hollow Lake, a small oxbow lake in Mississippi. Simulated sediment concentrations were compared with available field observations, with generally good agreement. The transport and resuspension processes of cohesive sediment due to wind-induced current and wave in Deep Hollow Lake were also discussed.     

Numerical study of three-dimensional shelf circulation on the Scotian Shelf using a shelf circulation model

A numerical shelf circulation model was developed for the Scotian Shelf, using a nested-grid setup consisting of a three-dimensional baroclinic inner model embedded inside a two-dimensional barotropic outer model. The shelf circulation model is based on the Princeton Ocean Model and driven by three-hourly atmospheric forcing provided by a numerical weather forecast model and by tidal forcing specified at the inner model's open boundaries based on pre-calculated tidal harmonic constants. The outer model simulates the depth-mean circulation forced by wind and atmospheric pressure fields over the northwest Atlantic Ocean with a horizontal resolution of 1/12°. The inner model simulates the three-dimensional circulation over the Gulf of St. Lawrence, the Scotian Shelf, and the adjacent slope with a horizontal resolution of 1/16°. The performance of the shelf circulation model is assessed by comparing model results with oceanographic observations made along the Atlantic coast of Nova Scotia and in the vicinity of Sable Island (on the Scotian Shelf) during two periods: October 2000–March 2001 and April–June 2002. Analysis of model results on Sable Island Bank indicates that tidal currents account for as much as ∼80% of the total variance of near-bottom currents, and currents driven by local winds account for ∼30% of the variance of the non-tidal near-bottom currents. Shelf waves generated remotely by winds and propagating into the region also play an important role in the near-bottom circulation on the bank.     

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.     

Numerical prediction of turbidity currents on unsteady sediment-supply

A three-dimensional(3D) numerical model of unstable turbidity currents is developed based on the mechanism of sediment transport and turbulence theory.In this model,numerical simulation of turbidity currents without subsequent supply of muddy water was conducted using the same parameters as were used in the flume experiments.The evolution process of turbidity currents of completely losing supplies observed in the experiment was simulated by the model;validation of the numerical model and the algorithm was conducted.If momentarily interrupted process is regarded as a special case of the gradually interrupted,based on the preceding numerical simulation validity,it is feasible to simulate the motion law of turbidity currents under losing gradually supplies.By this method in this article,the characteristic of sediment-laden flow of losing gradually supplies was obtained,as well as its relationship between front velocity and sediment concentration.     

Contemporary hydrodynamics and morphological change of a microtidal estuary: a numerical modelling study

Contemporary hydrodynamics and morphological change are examined in a shallow microtidal estuary, located on a wave-dominated coast (Port Stephens, NSW, Australia). Process-based numerical modelling is undertaken by combining modules for hydrodynamics, waves, sediment transport and bathymetry updates. Model results suggest that the complex estuarine bathymetry and geometry give rise to spatial variations in the tidal currents and a marked asymmetry between ebb and flood flows. Sediment transport paths correspond with tidal asymmetry patterns. The SE storms significantly enhance the quantities of sediment transport, while locally generated waves by the westerly strong winds also are capable of causing sediment entrainment and contribute to the delta morphological change. The wave/wind-induced currents are not uniform with flow over shoals driven in the same direction as waves/winds while a reverse flow occurring in the adjacent channel. The conceptual sediment transport model developed in this study shows flood-directed transport occurs on the flood ramp while ebb-directed net transport occurs in the tidal channels and at the estuary entrance. Accretion of the intertidal sand shoals and deepening of tidal channels, as revealed by the model, suggest that sediment-infilling becomes advanced, which may lead to an ebb-dominated estuary. It is likely that a switch from flood- to ebb-dominance occurs during the estuary evolution, and the present-day estuary acts as a sediment source rather than sediment sink to the coastal system. This is conflictive to the expectation drawn from the estuarine morphology; however, it is consistent with previous research suggesting that, in an infilling estuary, an increase in build-up of intertidal flats/shoals can eventually shift an estuary towards ebb dominance. Thus, field data are needed to validate the result presented here, and further study is required to investigate a variety of estuaries in the Australian area.     

Newtonian nudging for a Richards equation-based distributed hydrological model

The objective of data assimilation is to provide physically consistent estimates of spatially distributed environmental variables. In this study a relatively simple data assimilation method has been implemented in a relatively complex hydrological model. The data assimilation technique is Newtonian relaxation or nudging, in which model variables are driven towards observations by a forcing term added to the model equations. The forcing term is proportional to the difference between simulation and observation (relaxation component) and contains four-dimensional weighting functions that can incorporate prior knowledge about the spatial and temporal variability and characteristic scales of the state variable(s) being assimilated. The numerical model couples a three-dimensional finite element Richards equation solver for variably saturated porous media and a finite difference diffusion wave approximation based on digital elevation data for surface water dynamics. We describe the implementation of the data assimilation algorithm for the coupled model and report on the numerical and hydrological performance of the resulting assimilation scheme. Nudging is shown to be successful in improving the hydrological simulation results, and it introduces little computational cost, in terms of CPU and other numerical aspects of the model’s behavior, in some cases even improving numerical performance compared to model runs without nudging. We also examine the sensitivity of the model to nudging term parameters including the spatio-temporal influence coefficients in the weighting functions. Overall the nudging algorithm is quite flexible, for instance in dealing with concurrent observation datasets, gridded or scattered data, and different state variables, and the implementation presented here can be readily extended to any of these features not already incorporated. Moreover the nudging code and tests can serve as a basis for implementation of more sophisticated data assimilation techniques in a Richards equation-based hydrological model.     

MODELING OF THE HIGH CONCENTRATION LAYER OF COHESIVE SEDIMENT UNDER THE ACTION OF WAVES AND CURRENTS

MODELING OF THE HIGH CONCENTRATION LAYEROF COIIESIVE SEDMNT UNDER Tus ACTIONOF WAVES AND CURRENTSQinghe ZHANG', Yongsheng WU', Jiian LIAN1 and Pingxing DING3Abstract:High concenhation layer of cohesive sediment frequenhy occurs in muddy estUaries and coastalzones, and causes raPid siltation of the waterways. A one dimensional vertical coupled modeldescribing the interactions betWeen waves, currentS and suspended cohesive sediment is develoPed inthe pre…     

Sediment transport event analysis on the western Adriatic continental shelf

The sediment-transport mechanisms that contribute to and redistribute the modern sediment deposits on the western Adriatic continental shelf were evaluated utilizing data collected from two instrumented benthic tripods deployed at 12-m water depth, one in the northern Adriatic basin on the Po River subaqueous delta, and the other in the central Adriatic basin on the Pescara River shelf. Sediment-resuspension events driven by cold, northeasterly Bora winds dominate the along-shelf transport climatology at both tripod locations, but at the Po delta site, the southwesterly Scirocco wind events also play a significant role. At the Pescara shelf site, interaction between Bora wind-driven currents and the Western Adriatic Coastal Current strongly contributes to the resuspension and advection of suspended sediment. Interannual variability of the forcing mechanisms (including strength, frequency, and relative mix of Bora and Scirocco wind events) is evident in the three winters of data collected on the Po River subaqueous delta. In both types of wind events, and throughout all years of data collection, the net along-shelf sediment transport is significantly larger than the net across-shelf transport at the 12-m sites. This may be characteristic of low-energy environments, where sediment resuspension and transport occurs in such shallow water that it is not subjected to strong downwelling features characteristic of higher-energy environments.     

Post-glacial sediment dynamics in the Irish Sea and sediment wave morphology: Data–model comparisons

The irregular seafloor of the narrow Irish Sea on the NW European Shelf has been documented over several decades. From recently collected swath bathymetry data, very large trochoidal, nearly symmetrical sediment waves are observed in many parts of the Irish Sea and appear similar to those described from other continental shelf seas in North America that were covered by glacigenic sediments during the Last Glacial Maximum. Swath multibeam and single beam bathymetry data, backscatter intensity, shallow seismic imagery, video footage and sediment cores from the Irish Sea high sediment waves have been integrated to identify their genesis with reference to present and past hydrodynamic variability. From cross-sectional profiles over asymmetrical sediment waves in the Irish Sea the direction of asymmetry is used to map residual bed stress directions and associated bedload transport paths. Irish Sea peak bed stress vectors were generated using a two-dimensional palaeo-tidal model for the NW European shelf seas and compare well with the observations. Tidally induced bed stresses are modelled to have increased between 7–10 ka BP, to be nearly symmetrical in magnitude and to have reversed in dominant direction on a millennial scale. These environmental conditions during the post-glacial marine transgression are suggested here to help comprehend the construction of the very large sediment waves, with local variations due to differences in sediment grain size, sediment supply, water depth and intensified currents due to seafloor slopes. Model parameterisation using an open ocean boundary with time-dependent tidal changes and the implementation of high-resolution bathymetric information will improve future models of small-scale bed shear stress patterns and improve the predictive value of such modelling efforts.     

Dispersal of Mississippi and Atchafalaya sediment on the Texas-Louisiana shelf: Model estimates for the year 1993

A three-dimensional coupled hydrodynamic-sediment transport model for the Texas-Louisiana continental shelf was developed using the Regional Ocean Modeling System (ROMS) and used to represent fluvial sediment transport and deposition for the year 1993. The model included water and sediment discharge from the Mississippi River and Atchafalaya Bay, seabed resuspension, and suspended transport by currents. Input wave properties were provided by the Simulating WAves Nearshore (SWAN) model so that ROMS could estimate wave-driven bed stresses, critical to shallow-water sediment suspension. The model used temporally variable but spatially uniform winds, spatially variable seabed grain size distributions, and six sediment tracers from rivers and seabed.At the end of the year 1993, much of the modeled fluvial sediment accumulation was localized with deposition focused near sediment sources. Mississippi sediment remained within 20-40 km of the Mississippi Delta. Most Atchafalaya sediment remained landward of the 10-m isobath in the inner-most shelf south of Atchafalaya Bay. Atchafalaya sediment displayed an elongated westward dispersal pattern toward the Chenier Plain, reflecting the importance of wave resuspension and perennially westward depth-averaged currents in the shallow waters (<10 m). Due to relatively high settling velocities assumed for sediment from the Mississippi River as well as the shallowness of the shelf south of Atchafalaya Bay, most sediment traveled only a short distance before initial deposition. Little fluvial sediment could be transported into the vicinity of the “Dead Zone” (low-oxygen area) within a seasonal-annual timeframe. Near the Mississippi Delta and Atchafalaya Bay, alongshore sediment-transport fluxes always exceeded cross-shore fluxes. Estimated cumulative sediment fluxes next to Atchafalaya Bay were episodic and “stepwise-like” compared to the relatively gradual transport around the Mississippi Delta. During a large storm in March 1993, strong winds helped vertically mix the water column over the entire shelf (up to 100-m isobath), and wave shear stress dominated total bed stress. During fair-weather conditions in May 1993, however, the freshwater plumes spread onto a stratified water column, and combined wave-current shear stress only exceeded the threshold for suspending sediment in the inner-most part of the shelf.     

Storm-driven shelf-to-canyon suspended sediment transport at the southwestern Gulf of Lions

Shelf-to-canyon suspended sediment transport during major storms was studied at the southwestern end of the Gulf of Lions. Waves, near-bottom currents, temperature and water turbidity were measured on the inner shelf at 28-m water depth and in the Cap de Creus submarine canyon head at 300 m depth from November 2003 to March 2004. Two major storm events producing waves H_s>6 m coming from the E–SE sector took place, the first on 3–4 December 2003 (max H_s: 8.4 m) and the second on 20–22 February 2004 (max H_s: 7 m). During these events, shelf water flowed downcanyon producing strong near-bottom currents on the canyon head due to storm-induced downwelling, which was enhanced by dense shelf water cascading in February 2004. These processes generated different pulses of downcanyon suspended sediment transport. During the peak of both storms, the highest waves and the increasing near-bottom currents resuspended sediment on the canyon head and the adjacent outer shelf causing the first downcanyon sediment transport pulses. The December event ended just after these first pulses, when the induced downwelling finished suddenly due to restoration of shelf water stratification. This event was too short to allow the sediment resuspended on the shallow shelf to reach the canyon head. In contrast, the February event, reinforced by dense shelf water cascading, was long enough to transfer resuspended sediment from shallow shelf areas to the canyon head in two different pulses at the end of the event. The downcanyon transport during these last two pulses was one order of magnitude higher than those during the December event and during the first pulses of the February event and accounted for more than half of the total downcanyon sediment transport during the fall 2003 and winter 2004 period. Major storm events, especially during winter vertical mixing periods, produce major episodes of shelf-to-canyon sediment transport at the southwestern end of the Gulf of Lions. Hydrographic structure and storm duration are important factors controlling off-shelf sediment transport during these events.     

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.     

A 3-D finite volume model for sediment transport in coastal waters

Ocean Dynamics - A 3-D model has been developed to simulate sediment transport and bed change induced by currents and waves in coastal waters. The currents are calculated with the 3-D...     

Non-linear response of shoreface-connected sand ridges to interventions

A non-linear morphodynamic model of a microtidal coastal shelf is used to study the response of shoreface-connected sand ridges and the net sand balance of the shelf to large-scale interventions. The model describes the interaction between storm-driven currents and the erodible bottom. The transport of sediment comprises both bedload and suspended load contributions and is due to the joint action of waves (stirring of sediment from the bed) and net currents (causing transport). Three basic types of interventions are studied: extracting sand from ridges, nourishing sand at the shelf and constructing navigation channels. The model results indicate that for all interventions studied a relatively fast local recovery (time scale of decades to centuries) of the disturbed bathymetry to its original pattern takes place. Readjustment of the global system to its original equilibrium state (the saturation process) occurs on a longer time scale (several centuries). During the adjustment stage, significant net sand exchanges between inner shelf and adjacent outer shelf and near-shore zone occur. The results further suggest that extraction of sand from the shelf and dredging of navigation channels have negative implications for the stability of the beach (its sand volume decreases).Responsible Editor: Iris Grabemann     

Modeling undertow due to random waves

A numerical model of undertow due to random waves is developed. The model includes three sub-models: (i) a model for multi-directional and multi-frequency random wave transformation, (ii) a surface roller evolution model, and (iii) a model for calculating the vertical distribution and the mean value of the undertow velocity. The calculation of wave trough level is performed based on a theory for the wave asymmetry. The model was successfully validated against small- and large-scale laboratory experiments. Thus, the model is expected to provide reliable input for the modeling of sediment transport and morphological change due to waves and currents.     

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.     

Monitoring the physical forcing of resuspension events in the Thermaikos Gulf—NW Aegean during 2001–2003

This work constitutes an attempt to assess the relative importance of natural versus anthropogenic forcing for sediment resuspension on the shelf of the Thermaikos Gulf (NW Aegean) through a 2-year experiment in the framework of the E.U. project INTERPOL. Four periods of different hydrographic characteristics were identified, including two periods of stronger vertical homogenisation and two of stratified conditions. The former periods were characterized by stronger near-bed currents, while relatively weak internal wave motions dominated the periods of stratification. The near-bed currents showed strong coherence with the wind during the period of full homogenisation, whereas during stratification periods the wind provided indirect forcing evident mostly in the internal-wave bands. The site was too deep for the surface waves to cause any sediment resuspension. In conclusion, the observed near-bed currents at the site of interest did not appear to produce significant local resuspension of sediment; notably, the current-turbidity correlation suggested a shoreward transport of suspended material. The large-scale turbidity variability was rather related to the onset of the trawling period in mid-October 2001 and the increased riverine discharge in late summer/autumn 2002.     

Currents induced by vertical varied radiation stress in standing waves and evolution of the bed composed of fine sediments

This paper extends the conventional concept of radiation stress （Longuet-Higgins and Stewart, 1964） in progressive water waves to standing waves, so that its vertical profile could be defined and calculated in a new technical way. The hydrodynamic numerical model being coupled with the vertically varying radiation stress in standing waves is used to simulate the currents being induced by standing waves in the vertical section. Numerical modeling of suspended sediment transport is then carried out to simulate the evolution of the bed composed of fine sediments by the currents. The scour and deposition patterns simulated are in qualitative agreement with prior laboratory and field observations.     

Suspended sediment transport in the Gulf of Lions (NW Mediterranean): Impact of extreme storms and floods

In situ observations were combined with 3D modeling to gain understanding of and to quantify the suspended sediment transport in the Gulf of Lions (NW Mediterranean Sea). The outputs of a hydrodynamic–sediment transport coupled model were compared to near-bottom current and suspended sediment concentration measurements collected at the head of seven submarine canyons and at a shallow shelf site, over a 6-month period (November 2003–May 2004). The comparisons provide a reasonable validation of the model that reproduces the observed spatial and time variations. The study period was marked by an unusual occurrence of marine storms and high river inputs. The major water and sediment discharges were supplied by the Rhone, the largest Mediterranean river, during an exceptional flood accompanying a severe marine storm in early December 2003. A second major storm, with moderate flooding, occurred in February 2004. The estimate of river input during the studied period was 5.9 Mt. Our study reveals (i) that most of the particulate matter delivered by the Rhone was entrapped on the prodelta, and (ii) that marine storms played a crucial role on the sediment dispersal on the shelf and the off-shelf export. The marine storms occurring in early December 2003 and late February 2004 resuspended a very large amount of shelf sediment (>8 Mt). Erosion was controlled by waves on the inner shelf and by energetic currents on the outer shelf. Sediment deposition took place in the middle part of the shelf, between 50 and 100 m depth. Resuspended sediments and river-borne particles were transported to the southwestern end of the shelf by a cyclonic circulation induced by these onshore winds and exported towards the Catalan shelf and into the Cap de Creus Canyon which incises the slope close to the shore. Export taking place mostly during marine storms was estimated to reach 9.1 Mt during the study period.     

Physical controls on the dynamics of inlet sandbar systems

Knowledge of the physical processes acting at inlet systems and their interaction with sediments and sediment bodies is important to the understanding of such environments. The objectives of this study are to identify and assess the relative importance of the controlling processes across the complex sandbar system at the Teign inlet (Teignmouth, UK) through the combined application of a numerical model, field data and Argus video images. This allows the determination of the regions dominated by wave processes or by tidal processes and definition of the variability of these regions under different wave, tide and river-discharge conditions. Modelling experiments carried out for one stage of the evolution of the system show that the interaction between tidal motion and waves generates complex circulation patterns that drive the local sediment transport and sandbar dynamics, producing a cyclic morphological behaviour of the sandbars that form the ebb-tidal delta. The relative importance of each physical process on the sediment transport and consequent morphodynamics varies across the region. The main inlet channel is dominated by tidal action that directs the sediment transport as a consequence of the varying tidal flow asymmetry, resulting in net offshore transport. Sediment transport over the shoals and secondary channels at both sides of the main channel is dominated by wave-related processes, displacing sediment in the onshore direction. The interaction between waves and tide-generated currents controls the transport over the submerged sandbar that defines the channels seaward extend. High river discharge events are also proven to be important in this region, as they can change sediment-transport patterns across the area.Responsible Editor: Iris Grabemann