A two-dimensional horizontal wave propagation and mud mass transport model

The present study offers a two-dimensional horizontal wave propagation and morphodynamic model for muddy coasts. The model can be applied on a general three-dimensional bathymetry of a soft muddy coast to calculate wave damping, fluid mud mass transport and resulting bathymetry change under wave actions. The wave propagation model is based on time-dependent mild slope equations including the wave energy dissipation due to the wave-mud interaction of bottom mud layers as well as the combined effects of the wave refraction, diffraction and breaking. The constitutive equations of the visco-elastic–plastic model are adopted for the rheological behavior of fluid mud. The mass transport velocity within the fluid mud layer is calculated combining the Stokes’ drift, the mean Eulerian velocity and the gravity-driven mud flow. The results of the numerical model are compared against a series of conducted wave basin experiments, wave flume experiments and field observations. Comparisons between the computed results with both the field and laboratory data reveal the capability of the proposed model to predict the wave transformation and mud mass transport.     

Adaptive criterion curves describing incipient motion of sediment under wave and current conditions

Sediment incipient motion is a fundamental issue in sediment transport theory and engineering practice. Whilst Shields curve often is used to determine the threshold of sediment movement under unidirectional current conditions, it is unclear whether it can be directly applied for the wave or combined wave-current conditions. The study developed adaptive criterion curves describing incipient motion of sediment under wave and current conditions based on the flow pattern around the sediment particles. Firstly, the flow pattern law for fixed particles was recognized based on the friction law under various dynamic conditions (wave, current, and their combinations), and the flow pattern demarcations for incipient sediment motion were obtained with the threshold conditions for sediment movement under various dynamic conditions combined. Secondly, the exact shape of the Shields curve in each flow regime was derived under the current condition. By combining the flow pattern demarcations for incipient sediment motion under the wave condition, the criterion curve under the wave condition was derived. By combining the flow pattern demarcations for incipient sediment motion under the combined current-wave condition, the criterion curve for sediment incipient motion under the combined current-wave condition was derived. The results indicated that the flow pattern around incipient particles includes laminar, laminar-rough turbulent transition, and rough turbulent regimes. The criterion curves for sediment incipient motion under various dynamic conditions stayed the same in the laminar and rough turbulent regimes, but different in the transition regime. Depending on the relative strengths of the currents and waves, the shape of the criterion curve under the combined current-wave condition transitions adaptively between the criterion curve under the current condition and the criterion curve under the wave conditions.     

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.     

A two-phase numerical model for suspended-sediment transport in estuaries

This paper presents a full 2-D X/Z numerical model for sediment transport in open channels and estuaries using a two-phase (fluid–solid particle) approach. The physical concept and the mathematical background of the model are given and test-cases have been carried out to validate the proposed model. In order to illustrate its feasibility for a real estuary, the model has been applied to simulate the suspended-sediment transport and the formation of turbidity maximum in the Seine estuary. The numerical results show that the main characteristics of estuarine hydro-sediment dynamics in the Seine estuary are in fact reproduced by the proposed model. A qualitative agreement between the numerical results and the actual observations has been obtained and is presented in this paper.     

Sensitivity and uncertainty analysis of a sediment transport model: a global approach

Computerized sediment transport models are frequently employed to quantitatively simulate the movement of sediment materials in rivers. In spite of the deterministic nature of the models, the outputs are subject to uncertainty due to the inherent variability of many input parameters in time and in space, along with the lack of complete understanding of the involved processes. The commonly used first-order method for sensitivity and uncertainty analyses is to approximate a model by linear expansion at a selected point. Conclusions from the first-order method could be of limited use if the model responses drastically vary at different points in parameter space. To obtain the global sensitivity and uncertainty features of a sediment transport model over a larger input parameter space, the Latin hypercubic sampling technique along with regression procedures were employed. For the purpose of illustrating the methodologies, the computer model HEC2-SR was selected in this study. Through an example application, the results about the parameters sensitivity and uncertainty of water surface, bed elevation and sediment discharge were discussed.     

Steep wave,turbulence, and sediment concentration statistics beneath a breaking wave field and their implications for sediment transport

A new methodology based on wavelet analysis is used to estimate steep wave statistics under depth-limited conditions and the corresponding high concentration sediment statistics. Steep waves here are defined as wave crests within the wavelet transform exceeding a root mean square derived acceleration threshold. The method is applied to laboratory data obtained in a large-scale wave-flume experiment conducted in 2005 at Oregon State University's O. H. Hinsdale Wave Research Laboratory from an acoustic Doppler velocimeter and a fiber optic backscatter sensor array above a mobile sand bed. The steep wave and high concentration statistical results for the erosive condition suggest that sand suspensions are intermittent when a wave-breaking timescale (the ratio of breaking wave height and rms wave velocity) is used to detect the concurrence among steep wave, high velocity turbulent fluctuations, and sand concentration events near the bed. More importantly, at 1 cm above the bed, though the accretive case has more steep wave events, the erosive case has more steep waves and concurrent high concentration events, suggesting a more intense breaking wave process near the sensors. The use of a longer time window, based on the dominant wave period in the detection process of steep wave and high concentration events at 1 cm above the bed, does not change the resulting statistics for the erosive condition. However, increased percentages of high concentration events correlated with steep wave and high velocity turbulence events for the accretive condition are obtained. These increased percentages are conjectured to be due to advection of non-local turbulent events and sediment concentration peaks from upstream.     

Parameterization of near-bed processes under collinear wave and current flows from a two-phase sheet flow model

Sediment transport models require appropriate representation of near-bed processes. We aim here to explore the parameterizations of bed shear stress, bed load transport rate and near-bed sediment erosion rate under the sheet flow regime. To that end, we employ a one-dimensional two-phase sheet flow model which is able to resolve the intrawave boundary layer and sediment dynamics at a length scale on the order of the sediment grain. We have conducted 79 numerical simulations to cover a range of collinear wave and current conditions and sediment diameters in the range 210–460 μm$\mathrm{\mu }\mathrm{m}$. The numerical results confirm that the intrawave bed shear stress leads the free stream velocity, and we assess an explicit expression relating the phase lead to the maximum velocity, wave period and bed roughness. The numerical sheet flow model is also used to provide estimates for the bed load transport rate and to inspect the near-bed sediment erosion. A common bed load transport rate formulation and two typical reference concentration approaches are assessed. A dependence of the bed load transport rate on the sediment grain diameter is observed and parameterized. Finally, the intrawave near-bed vertical sediment flux is further investigated and related to the time derivative of the bed shear stress.     

Fine sediment carrying capacity of combined wave and current flows

The so-called fine sediment in many coastal areas and estuaries in China is mostly referred to the mixture of cohesive sediment and non-cohesive sediment. To predict the mixed type time sediment transport, sediment carrying capacity formulae combined with the 2-D suspended sediment transport equation and morphologic equation have been widely used in China. In the present study, the sediment carrying capacity formula suggested by Dou et al. （1995） for wave conditions has been improved and implemented for the prediction of sediment transport in nearshore regions where wave activities are significant. The improvement is based on the wave energy dissipation principle inside and outside the surf zone. In the improved formula, sediment in suspension increases with the magnitude of the wave period and this feature complies with general observations. More than 300 laboratory and field measured data sets have been reviewed and 12 of them have been used to verify and determine the major coefficients in the improved formula. The application of the sediment carrying capacity model in combined wave and current situations shows that the model can faithfully reproduce the cross-shore sediment concentration distributions at the southwest coast of Bohai Bay.     

Modelling gravel-bed topography from sediment transport

A simple queueing model which generates bed topography consistent with the mechanics of gravel motion is presented. The equations on which it is based are derived from the application of simple theoretical ideas and an analysis of published flume data (Meland and Norrman, 1966; Francis, 1973; Abbott and Francis, 1977). The simulation produces an extremely variable bed topography in which at least two scales of bedform may be identified. Features of 5–10 grain wavelengths, similar to pebble clusters, are superimposed on other forms whose wavelength is up to 30–40 grains and which resemble step-pools systems and antidunes. Sediment transport shows many of the characteristics found in the field and the model is, therefore, thought to provide a useful starting point from which to examine the interaction between flow, sediment transport, and bedforms in gravel-bed and cobble streams.     

One-dimensional numerical simulation of non-uniform sediment transport under unsteady flows

One-dimensional numerical models are popularly used in sediment transport research because they can be easily programmed and cost less time compared with two- and three-dimensional numerical models. In particular, they possess greater capacity to be applied in large river basins with many tributaries. This paper presents a one-dimensional numerical model capable of calculating total-load sediment transport. The cross-section-averaged sediment transport capacity and recovery coefficient are addressed in the suspended load model. This one-dimensional model, therefore, can be applied to fine suspended loads and to hyperconcentrated flows in the Yellow River. Moreover, a new discretization scheme for the equation of unsteady non-uniform suspended sediment transport is proposed. The model is calibrated using data measured from the Yantan Reservoir on the Hongshui River and the Sanmenxia Reservoir on the Yellow River. A comparison of the calculated water level and river bed deformation with field measurements Shows that the improved numerical model is capable of predicting flow, sediment transport, bed changes, and bed-material sorting in various situations, with reasonable accuracy and reliability.     

Two-dimensional modeling of fine sediment transport with mixed sediment and consolidation: Application to the Gironde Estuary,France

In order to optimize ship navigation in the macrotidal Gironde Estuary, a recent project funded by the port of Bordeaux aims at better understand and forecast hydrodynamic and fine sediment transport within the estuary. In the framework of this project, a two-dimensional hydro-sedimentary model is built. The model includes hydrodynamic forcings, mixed-sediment transport, and consolidation processes. The harmonic analysis of the astronomical tides reveals a strong distortion of the tidal wave inducing the growth of overtide constituents and the non-significant effect of tide-surge interactions in annual-scale prediction. Depending on hydrological conditions, river discharge can considerably alter the model accuracy due to the migration of the turbidity maximum zone modifying the bottom roughness. Comparison with measurements shows the ability of the model to reproduce suspended-sediment concentrations in the central Estuary. Sensitivity of the model to sediment features has also been discussed in regard of suspended-sediment concentrations and fluid mud deposits. The model will be further coupled with ship squat predictions and a morphodynamic model.     

Developing a Lagrangian sediment transport model for open channel flows

A three-dimensional stochastic Lagrangian particle tracking sediment transport model is developed to solve the discrete advection-dispersion equation using a combination of empirical dispersion equations.The performance of three widely-used longitudinal dispersion coefficient equations was examined to select one of them as the primary dispersion equation term in the developed model. Also, a conditional empirical equation was used to consider the effect of vertical dispersion term in top layers n...     

Sediment transport and resuspension due to combined motion of wave and current in the northern Adriatic Sea during a Bora event in January 2001: A numerical modelling study

The Adriatic Sea general circulation model coupled to a third generation wave model SWAN and a sediment transport model was implemented in the Adriatic Sea to study the dynamics of the sediment transport and resuspension in the northern Adriatic Sea (NAS) during the Bora event in January 2001. The bottom boundary layer (BBL) was resolved by the coupled model with high vertical resolution, and the mechanism of the wave–current interaction in the BBL was also represented in the model. The study found that, during the Bora event of 13–17 January 2001, large waves with significant wave height 2 m and period of 5 s were generated by strong winds in the northwestern shelf of the Adriatic where the direction of wave propagation was orthogonal to the current. The combined motion of the wave and current in the BBL increased the bottom stress over the western Adriatic shelf, resulting in stronger sediment resuspension there. Combining stronger bottom resuspension and strong upward vertical flux of resuspended sediments due to turbulent mixing, the model predicted that sediment concentration near the Po River was much higher than that predicted by the model run without wave forcing. The study also shows that wave–current interaction in the BBL reduced the western Adriatic Coastal Currents (WACCs) in the shallower north. It is concluded that wave forcing significantly changed the sediment distributions and increased the total horizontal fluxes over the western shelf. These results signified wave effect on sediment flux and distribution in the NAS, and suggested that waves cannot be neglected in the study of dynamics of sediment transport and resuspension in the shallow coastal seas. By including the tidal forcing in the coupled model, we also examined the effect of tides on the sediment transport dynamics in the NAS.     

Fuzzy logic algorithm for runoff-induced sediment transport from bare soil surfaces

Utilizing the rainfall intensity, and slope data, a fuzzy logic algorithm was developed to estimate sediment loads from bare soil surfaces. Considering slope and rainfall as input variables, the variables were fuzzified into fuzzy subsets. The fuzzy subsets of the variables were considered to have triangular membership functions. The relations among rainfall intensity, slope, and sediment transport were represented by a set of fuzzy rules. The fuzzy rules relating input variables to the output variable of sediment discharge were laid out in the IF-THEN format. The commonly used weighted average method was employed for the defuzzification procedure.The sediment load predicted by the fuzzy model was in satisfactory agreement with the measured sediment load data. Predicting the mean sediment loads from experimental runs, the performance of the fuzzy model was compared with that of the artificial neural networks (ANNs) and the physics-based models. The results of showed revealed that the fuzzy model performed better under very high rainfall intensities over different slopes and over very steep slopes under different rainfall intensities. This is closely related to the selection of the shape and frequency of the fuzzy membership functions in the fuzzy model.     

A non‐equilibrium sediment transport model for rill erosion

Sediment transport in rill flows exhibits the characteristics of non‐equilibrium transport, and the sediment transport rate of rill flow gradually recovers along the flow direction by erosion. By employing the concept of partial equilibrium sediment transport from open channel hydraulics, a dynamic model of rill erosion on hillslopes was developed. In the model, a parameter, called the restoration coefficient of sediment transport capacity, was used to express the recovery process of sediment transport rate, which was analysed by dimensional analysis and determined from laboratory experimental data. The values of soil loss simulated by the model were in agreement with observed values. The model results showed that the length and gradient of the hillslope and rainfall intensity had different influences on rill erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Field studies of sediment transport by overland flow

Field studies on sandy soils of the Cottenham Series in mid-Bedfordshire show that the mean annual rate of sediment transport by overland flow on an 11° mid-slope is 98 g cm¹. The feasibility of using sediment transport equations to predict erosion by overland flow on a storm basis is examined by comparing the observed values of sediment yield with values predicted by four sediment transport equations and a regression equation which relates soil loss to runoff energy and rainfall energy. An expression combining Engelund's sediment transport capacity equation and the Manning equation for flow velocity, as modified by Savat for disturbed flow, best reflects field conditions. Although there is a significant correlation (r = 0.69; n 30) between the observed and predicted values using this expression, the coefficient of determination is too low for predictive purposes. Reasons for this are presented.     

A probabilistic/deterministic approach for the prediction of the sediment transport rate

In this paper we present a probabilistic/deterministic model for the evaluation of the sediment transport rate in a stream. Starting from Einstein’s theory, the approach was obtained by trying to overcome some of the intrinsic limitations. The approach is based on two distinct probability functions, one relevant to the detachment of grains and the second relevant to the length of particle jumps. The sediment transport rate is obtained by integrating the distribution of the ranges of the particle jumps multiplied by the average particle velocity. The relationship for the average ranges of particle jumps is an opportune combination of the Einstein and Yalin expressions. The final formulation was calibrated by means of a large number of experimental data and also by comparison with some of the most widely-used empirical formulas. The results show a better agreement between theory and experiments than do the other theories analyzed.     

Approximate velocity formula over mobile sediment bed induced by velocity-skewed waves and current

A velocity formula is proposed for flow over a mobile sediment bed induced by velocity-skewed waves and current. The formula is obtained by a separation of waves and current velocities and requires seven free variables related to free stream velocity and sediment characteristics. The formula includes two parts:(1) a wave part consisting of the free stream velocity and defect function, which considers phase lead, wave boundary layer thickness, and mobile bed level, and(2) a current part, which ch...     

LES-DEM simulations of sediment transport

In this work, a fully-coupled Computational Fluid Dynamics (CFD) model and Discrete Element Method (DEM) are used to simulate a unidirectional turbulent open-channel flow over the full range of sediment transport regimes. The fluid and particles are computed on separate grids using a dual-grid formulation to maintain consistency and avoid instability issues. The results of coupling the dispersed phase to a multiphase flow solver that uses volume-averaged Navier-Stokes equations are compared to those obtained from coupling through drag to a single flow solver. The current work also examines the applicability and limitations of lumping particles as a representative particle to reduce the cost of simulations. Insight to the impact of different turbulent events to the entrainment of particles is also given. The simulation results of sediment transport from both coupling techniques show good agreement with empirical formulas in the bedload regime, but under-predict sediment transport in the suspended load regime. In the suspended load regime, using partial coupling, the rate of sediment transport was found to be under-predicted as compared to full-coupling. The deviation in results in the suspended load regime was found to increase with increases in the applied shear stress. Both coupling methods revealed the same effect on the friction factor where friction increases in the bedload regime and decreases in the suspended load regime reaching a maximum at the transition between regimes. This result is contrary to past studies which have shown a discrete jump in the friction factor at the transition. Lumping particles as representative particles is shown to reduce the simulation cost by more than a factor of 5 when using a scaling factor of 2. By doing a quadrant analysis on information obtained from particle and flow field results, it was found that most of the particles are entrained by more frequent sweep events.     

Modeling the impact of climate change on sediment transport and morphology in coupled watershed-coast systems:A case study using an integrated approach

Climate change is an issue of major concern nowadays.Its impact on the natural and human environment is studied intensively,as the expected shift in climate will be significant in the next few decades.Recent experience shows that the effects will be critical in coastal areas,resulting in erosion and inundation phenomena worldwide.In addition to that,coastal areas are subject to ＂pressures＂ from upstream watersheds in terms of water quality and sediment transport.The present paper studies the impact of climate change on sediment transport and morphology in the aforementioned coupled system.The study regards a sandy coast and its upstream watershed in Chalkidiki,North Greece;it is based on：（a）an integrated approach for the quantitative correlation of the two through numerical modeling,developed by the authors,and（b）a calibrated application of the relevant models Soil and Water Assessment Tool（SWAT）and PELNCON-M,applied to the watershed and the coastal zone,respectively.The examined climate change scenarios focus on a shift of the rainfall distribution towards fewer and more extreme rainfall events,and an increased frequency of occurrence of extreme wave events.Results indicate the significance of climatic pressures in wide-scale sediment dynamics,and are deemed to provide a useful perspective for researchers and policy planners involved in the study of coastal morphology evolution in a changing climate.