Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

Evaluating competing hypotheses for the origin and dynamics of river anastomosis

Anastomosing rivers have multiple interconnected channels that enclose flood basins. Various theories potentially explain this pattern, including an increased discharge conveyance and sediment transport capacity of multiple channels, deltaic branching, avulsion forced by base‐level rise, or a tendency to avulse due to upstream sediment overloading. The former two imply a stable anabranching channel pattern, whereas the latter two imply disequilibrium and evolution towards a single‐channel pattern in the absence of avulsion. Our objective is to test these hypotheses on morphodynamic scenario modelling and data of a well‐documented case study: the upper Columbia River. Proportions of channel and floodplain sediments along the river valley were derived from surface mapping. Initial and boundary conditions for the modelling were derived from field data. A 1D network model was built based on gradually varied flow equations, sediment transport prediction, mass conservation, transverse slope and spiral meander flow effects at the bifurcations. The number of channels and crevasse splays decreases in a downstream direction. Also, measured sediment transport is higher at the upstream boundary than downstream. These observations concur with bed sediment overloading from upstream, which can have caused channel aggradation above the surrounding floodplain and subsequent avulsion. The modelling also indicates that avulsion was likely caused by upstream overloading. In the model, multi‐channel systems inevitably evolve towards single‐channel systems within centuries. The reasons are that symmetric channel bifurcations are inherently unstable, while confluenced channels have relatively less friction than two parallel channels, so that more discharge is conveyed through the path with more confluences and less friction. Furthermore, the present longitudinal profile curvature of the valley could only be reproduced in the model by temporary overfeeding. We conclude that this anastomosing pattern is the result of time‐varying sediment overloading and is not an equilibrium pattern feature, and suggest this is valid for many anastomosing rivers. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling dam-break flows over mobile beds using a 2D coupled approach

Dam-break flows usually propagate along rivers and floodplains, where the processes of fluid flow, sediment transport and bed evolution are closely linked. However, the majority of existing two-dimensional (2D) models used to simulate dam-break flows are only applicable to fixed beds. Details are given in this paper of the development of a 2D morphodynamic model for predicting dam-break flows over mobile beds. In this model, the common 2D shallow water equations are modified, so that the effects of sediment concentrations and bed evolution on the flood wave propagation can be considered. These equations are used together with the non-equilibrium transport equations for graded sediments and the equation of bed evolution. The governing equations are solved using a matrix method, thus the hydrodynamic, sediment transport and morphological processes can be jointly solved. The model employs an unstructured finite volume algorithm, with an approximate Riemann solver, based on the Roe-MUSCL scheme. A predictor–corrector scheme is used in time stepping, leading to a second-order accurate solution in both time and space. In addition, the model considers the adjustment process of bed material composition during the morphological evolution process. The model was first verified against results from existing numerical models and laboratory experiments. It was then used to simulate dam-break flows over a fixed bed and a mobile bed to examine the differences in the predicted flood wave speed and depth. The effects of bed material size distributions on the flood flow and bed evolution were also investigated. The results indicate that there is a great difference between the dam-break flow predictions made over a fixed bed and a mobile bed. At the initial stage of a dam-break flow, the rate of bed evolution could be comparable to that of water depth change. Therefore, it is often necessary to employ the turbid water governing equations using a coupled approach for simulating dam-break flows.     

Transport of moderately sorted gravel at low bed shear stresses: The role of fine sediment infiltration

A reliable estimation of sediment transport in gravel‐bed streams is important for various practical engineering and biological studies (e.g., channel stability design, bed degradation/aggradation, restoration of spawning habitat). In the present work, we report original laboratory experiments investigating the transport of gravel particles at low bed shear stresses. The laboratory tests were conducted under unsteady flow conditions inducing low bed shear stresses, with detailed monitoring of the bed topography using a laser scanner. Effects of bed surface arrangements were documented by testing loose and packed bed configurations. Effects of fine sediments were examined by testing beds with sand, artificial fine sand or cohesive silt infiltrated in the gravel matrix. Analysis of the experimental data revealed that the transport of gravel particles depends upon the bed arrangement, the bed material properties (e.g., size and shape, consolidation index, permeability) and the concentration of fine sediments within the surface layer of moving grains. This concentration is directly related to the distribution of fine particles within the gravel matrix (i.e., bottom‐up infiltration or bridging) and their transport mode (i.e., bedload or suspended load). Compared to loose beds, the mobility of gravel is reduced for packed beds and for beds clogged from the bottom up with cohesive fine sediments; in both cases, the bed shear stress for gravel entrainment increases by about 12%. On the other hand, the mobility of gravel increases significantly (bed shear stress for particle motion decreasing up to 40%) for beds clogged at the surface by non‐cohesive sand particles. Copyright © 2017 John Wiley & Sons, Ltd.     

Tidal straining effect on the suspended sediment transport in the Huanghe (Yellow River) Estuary, China

Tidal straining effect on sediment transport dynamics in the Huanghe (Yellow River) estuary was studied by field observations and numerical simulations. The measurement of salinity, suspended sediment concentration, and current velocity was conducted during a flood season in 1995 at the Huanghe river mouth with six fishing boats moored at six stations for 25-h hourly time series observations. Based on the measurements, the intra-tidal variations of sediment transport in the highly turbid river mouth was observed and the tidal straining effect occurred. Our study showed that tidal straining of longitudinal sediment concentration gradients can contribute to intra-tidal variability in sediment stratification and to asymmetries in sediment distribution within a tidal cycle. In particular, the tidal straining effect in the Huanghe River estuary strengthened the sediment-induced stratification at the flood tide, thus producing a higher bottom sediment concentration than that during the ebb. A sediment transport model that is capable of simulating sediment-induced stratification effect on the hydrodynamics in the bottom boundary layers and associated density currents was applied to an idealized estuary to demonstrate the processes and to discuss the mechanism. The model-predicted sediment processes resembled the observed characteristics in the Huanghe River estuary. We concluded that tidal straining effect is an important but poorly understood mechanism in the transport dynamics of cohesive sediments in turbid estuaries and coastal seas.     

A GPU-based numerical model coupling hydrodynamical and morphological processes

Sediment transport simulations are important in practical engineering. In this study, a graphics processing unit (GPU)-based numerical model coupling hydrodynamical and morphological processes was developed to simulate water flow, sediment transport, and morphological changes. Aiming at accurately predicting the sediment transport and sediment scouring processes, the model resolved the realistic features of sediment transport and used a GPU-based parallel computing technique to the accelerate calculation. This model was created in the framework of a Godunov-type finite volume scheme to solve the shallow water equations (SWEs). The SWEs were discretized into algebraic equations by the finite volume method. The fluxes of mass and momentum were computed by the Harten, Lax, and van Leer Contact (HLLC) approximate Riemann solver, and the friction source terms were calculated by the proposed a splitting point-implicit method. These values were evaluated using a novel 2D edge-based MUSCL scheme. The code was programmed using C++ and CUDA, which could run on GPUs to substantially accelerate the computation. The aim of the work was to develop a GPU-based numerical model to simulate hydrodynamical and morphological processes. The novelty is the application of the GPU techniques in the numerical model, making it possible to simulate the sediment transport and bed evolution in a high-resolution but efficient manner. The model was applied to two cases to evaluate bed evolution and the effects of the morphological changes on the flood patterns with high resolution. This indicated that the GPU-based high-resolution hydro-geomorphological model was capable of reproducing morphological processes. The computational times for this test case on the GPU and CPU were 298.1 and 4531.2 s, respectively, indicating that the GPU could accelerate the computation 15.2 times. Compared with the traditional CPU high-grid resolution, the proposed GPU-based high-resolution numerical model improved the reconstruction speed more than 2.0–12.83 times for different grid resolutions while remaining computationally efficient.     

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.     

Effects of vegetation on flow and sediment transport: comparative analyses and validation of predicting models

The presence of vegetation modifies flow and sediment transport in alluvial channels and hence the morphological evolution of river systems. Plants increase the local roughness, modify flow patterns and provide additional drag, decreasing the bed‐shear stress and enhancing local sediment deposition. For this, it is important to take into account the presence of vegetation in morphodynamic modelling. Models describing the effects of vegetation on water flow and sediment transport already exist, but comparative analyses and validations on extensive datasets are still lacking. In order to provide practical information for modelling purposes, we analysed the performance of a large number of models on flow resistance, vegetation drag, vertical velocity profiles and bed‐shear stresses in vegetated channels. Their assessments and applicability ranges are derived by comparing their predictions with measured values from a large dataset for different types of submerged and emergent vegetation gathered from the literature. The work includes assessing the performance of the sediment transport capacity formulae of Engelund and Hansen and van Rijn in the case of vegetated beds, as well as the value of the drag coefficient to be used for different types of vegetation and hydraulic conditions. The results provide a unique comparative overview of existing models for the assessment of the effects of vegetation on morphodynamics, highlighting their performances and applicability ranges. Copyright © 2014 John Wiley & Sons, Ltd.     

The settling dynamics of flocculating mud and sand mixtures: part 2—numerical modelling

Estuarine and coastal sediment transport is characterised by the transport of both sand-sized particles (of diameter greater than 63?μm) and muddy fine-grained sediments (silt, diameter less than 63?μm; clay, diameter less than 2?μm). These fractions are traditionally considered as non-cohesive and cohesive, respectively, because of the negligible physico-chemical attraction that occurs between sand grains. However, the flocculation of sediment particles is not only caused by physico-chemical attraction. Cohesivity of sediment is also caused by biology, in particular the sticky extra-cellular polymeric substances secreted by diatoms, and the effect of biology in binding sediment particles can be much larger than that of physico-chemical attraction. As demonstrated by Manning (2008) and further expanded in part 1 of this paper (Manning et al., submitted), the greater binding effect of biology allows sand particles to flocculate with mud. In many estuaries, both the sand and fine sediment fractions are transported in significant quantities. Many of the more common sediment transport modelling suites now have the capability to combine mud and sand transport. However, in all of these modelling approaches, the modelling of mixed sediment transport has still essentially separated the modelling of sand and mud fractions assuming that these different fractions do not interact except at the bed. However, the use of in situ video techniques has greatly enhanced the accuracy and reliability of settling velocity measurements and has led to a re-appraisal of this widely held assumption. Measurements of settling velocity in mixed sands presented by Manning et al. (2009) have shown strong evidence for the flocculation of mixed sediments, whilst the greater understanding of the role of biology in flocculation has identified mechanisms by which this mud-sand flocculation can occur. In the first part of this paper (Manning et al., submitted), the development of an empirical flocculation model is described which represents the interaction between sand and mud particles in the flocculation process. Measurements of the settling velocity of varying mud-sand mixtures are described, and empirical algorithms governing the variation of settling velocity with turbulence, suspended sediment concentration and mud-sand content are derived. The second part of this paper continues the theme of examination of the effects of mud-sand interaction on flocculation. A 1DV mixed transport model is developed and used to reproduce the vertical transport of mixed sediment fractions. The 1DV model is used to reproduce the measured settling velocities in the laboratory experiments described in the part 1 paper and also to reproduce measurements of concentration of mixed sediments in the Outer Thames. In both modelling exercises, the model is run using the algorithms developed in part 1 and repeated using an assumption of no interaction between mud and sand in the flocculation process. The results of the modelling show a significant improvement in the ability of the 1DV to reproduce the observed sediment behaviour when the empirical equations are used. This represents further strong evidence of the interaction between sand and mud in the flocculation process.     

Tools and cover effects in bedload transport observations in the Pitzbach,Austria

We report bedload data and acoustic impulse measurements due to particle impact from the Pitzbach in Austria. Impulse counts can be viewed as a measure of the energy delivered to the bed by moving particles. Impulse counts show a large scatter even for the same discharge and bedload supply. This scatter is due to varying grain size distribution, grain shape, mode of transport of the sediment particles and spatial and temporal distribution of the sediment load. The mean impulse count at given hydraulic conditions may increase or decrease with increasing sediment supply, suggesting that both tools and cover effects are active on the channel bed. Dependent on the local balance between sediment supply and transport capacity, either effect may be dominant at different locations along the cross‐section at the same time. Furthermore, the same bed location may respond to increasing sediment supply as tools‐dominated at some discharges and cover‐dominated at other discharges. Our observations may have implications for modelling of bedrock erosion in landscape evolution models and of bedrock channel morphology. Erosion models that do not incorporate both tools and cover effects are not sufficient to describe observations. Furthermore, a local erosion law cannot in general be used to describe erosion averaged over the channel cross‐section. The changing balance between sediment supply and transport capacity with increasing discharge highlights that a single representative discharge is not sufficient to capture the full erosion dynamics. Copyright © 2008 John Wiley & Sons, Ltd.     

Influence of bed clusters and size gradation on operational time distribution for non‐uniform bed‐load transport

The operational time distribution (OTD) defines the time for bed‐load sediment spent in motion, which is needed to characterize the random nature of sediment transport. This study explores the influence of bed clusters and size gradation on OTD for non‐uniform bed‐loads. First, both static and mobile bed armouring experiments were conducted in laboratorial flumes to monitor the transport of mixed sand/gravel sediments. Only in the mobile armouring experiment did apparent bed clusters develop, because of stable feeding and a longer transport period. Second, a generalized subordinated advection (GSA) model was applied to quantify the observed dynamics of tracer particles. Results show that for the static armour layer (without sediment feed), the best‐fit OTD assigns more weight to the large displacement of small particles, likely because of the size‐selective entrainment process. The capacity coefficient in the GSA model, which affects the width of the OTD, is space dependent only for small particles whose dynamics can be significantly affected by larger particles and whose distribution is more likely to be space dependent in a mixed sand and gravel system. However, the OTD for the mobile armour layer (with sediment recirculation) exhibited longer tails for larger particles. This is because the trailing edge of larger particles is more resistant to erosion, and their leading front may not be easily trapped by self‐organized bed clusters. The strong interaction between particle–bed may cause the capacity coefficient to be space‐dependent for bed‐load transport along mobile armour layers. Therefore, the combined laboratory experiments and stochastic model analysis show that the OTD may be affected more by particle–bed interactions (such as clusters) than by particle–particle interactions (e.g. hiding and exposing), and that the GSA model can quantify mixed‐size sand/gravel transport along river beds within either static or mobile armour layers. Copyright © 2016 John Wiley & Sons, Ltd.     

Spatial heterogeneity in estuarine mud dynamics

The fate of mud in an estuary over an entire year was unravelled using complementary, independent, spatially explicit techniques. Sequential ERS-2 SAR and Envisat MERIS-FR data were used to derive synoptic changes in intertidal bottom mud and suspended particulate matter (SPM) in the top of the water column, respectively. These satellite data were combined with in situ measurements and with a high resolution three-dimensional cohesive sediment model, simulating mud transport, resuspension, settling and deposition under the influence of tides, wind, waves and freshwater discharge. The spatial distribution of both bottom mud and SPM as observed by in situ and satellite techniques was largely explained by modelled estuarine circulation, tidal and wind-induced variations in vertical mixing and horizontal advection. The three data sources also showed similar spring-neap and seasonal variations in SPM (all factor 1.5 to 2), but semi-diurnal tidal variations were underestimated by the model. Satellite data revealed that changes in intertidal bottom mud were spatially heterogeneous, but on average mud content doubled during summer, which was confirmed by in situ data. The model did not show such seasonal variation in bed sediment, suggesting that seasonal dynamics are not well explained by the physical factors presently implemented in the model, but may be largely attributed to other (internal) factors, including increased floc size in summer, temporal stabilisation of the sediment by microphytobenthos and a substantially lower roughness of the intertidal bed in summer as observed by the satellite. The effects of such factors on estuarine mud dynamics were evaluated.     

Numerical modeling of bedload and suspended load contributions to morphological evolution of the Seine Estuary (France)

This numerical modeling study (i) assesses the influence of the sediment erosion process on the sediment dynamics and subsequent morphological changes of a mixed-sediment environment, the macrotidal Seine estuary, when non-cohesive particles are dominant within bed mixtures (non-cohesive regime), and (ii) investigates respective contributions of bedload and suspended load in these dynamics. A three dimensional (3D) process-based morphodynamic model was set up and run under realistic forcings (including tide, waves, wind, and river discharge) during a 1-year period. Applying erosion homogeneously to bed sediment in the non-cohesive regime, i.e., average erosion parameters in the erosion law (especially the erodibility parameter, E₀), leads to higher resuspension of fine sediment due to the presence of coarser fractions within mixtures, compared to the case of an independent treatment of erosion for each sediment class. This results in more pronounced horizontal sediment flux (two-fold increase for sand, +30% for mud) and erosion/deposition patterns (up to a two-fold increase in erosion over shoals, generally associated with some coarsening of bed sediment). Compared to observed bathymetric changes, more relevant erosion/deposition patterns are derived from the model when independent resuspension fluxes are considered in the non-cohesive regime. These results suggest that this kind of approach may be more relevant when local grain-size distributions become heterogeneous and multimodal for non-cohesive particles. Bedload transport appears to be a non-dominant but significant contributor to the sediment dynamics of the Seine Estuary mouth. The residual bedload flux represents, on average, between 17 and 38% of the suspended sand flux, its contribution generally increasing when bed sediment becomes coarser (can become dominant at specific locations). The average orientation of residual fluxes and erosion/deposition patterns caused by bedload generally follow those resulting from suspended sediment dynamics. Sediment mass budgets cumulated over the simulated year reveal a relative contribution of bedload to total mass budgets around 25% over large erosion areas of shoals, which can even become higher in sedimentation zones. However, bedload-induced dynamics can locally differ from the dynamics related to suspended load, resulting in specific residual transport, erosion/deposition patterns, and changes in seabed nature.     

A continuous model for sand dunes: Review,new developments and application to barchan dunes and barchan dune fields

Basically, sand dunes are patterns resulting from the coupling of hydrodynamic and sediment transport. Once grains move, they modify the surface topography which in turns modifies the flow. This important feedback mechanism lies at the core of continuous dune modelling. Here we present an updated review of such a model for aeolian dunes, including important modifications to improve its predicting power. For instance, we add a more realistic wind model and provide a self‐consistent set of parameters independently validated. As an example, we are able to simulate realistic barchan dunes, which are the basic solution of the model in the condition of unidirectional flow and scarce sediments. From the simulation, we extract new relations describing the morphology and dynamics of barchans that compare very well with existing field data. Next, we revisit the problem of the stability of barchan dunes and argue that they are intrinsically unstable bed‐forms. Finally, we perform more complex simulations: first, a barchan dune under variable wind strength and, second, barchan dune fields under different boundary conditions. The latter has important implications for the problem of the genesis of barchan dunes. Copyright © 2010 John Wiley & Sons, Ltd.     

A coupled hierarchical modeling approach to simulating the geomorphic response of river systems to anthropogenic climate change

Anthropogenic climate change is expected to change the discharge and sediment transport regime of river systems. Because rivers adjust their channels to accommodate their typical inputs of water and sediment, changes in these variables can potentially alter river morphology. In this study, a hierarchical modeling approach was developed and applied to examine potential changes in reach‐averaged bedload transport and spatial patterns of erosion and deposition for three snowmelt‐dominated gravel‐bed rivers in the interior Pacific Northwest. The modeling hierarchy was based on discharge and suspended‐sediment load from a basin‐scale hydrologic model driven by a range of downscaled climate‐change scenarios. In the field, channel morphology and sediment grain‐size data for all three rivers were collected. Changes in reach‐averaged bedload transport were estimated using the Bedload Assessment of Gravel‐bedded Streams (BAGS) software, and the Cellular Automaton Evolutionary Slope and River (CAESAR) model was used to simulate the spatial pattern of erosion and deposition within each reach to infer potential changes in channel geometry and planform. The duration of critical discharge was found to control bedload transport. Changes in channel geometry were simulated for the two higher‐energy river reaches, but no significant morphological changes were found for a lower‐energy reach with steep, cohesive banks. Changes in sediment transport and river morphology resulting from climate change could affect the management of river systems for human and ecological uses. Copyright © 2015 John Wiley & Sons, Ltd.     

The response of a gravel‐bed river planform configuration to flow variations and bed reworking: a modelling study

A 2D depth‐averaged model has been developed for simulating water flow, sediment transport and morphological changes in gravel‐bed rivers. The model was validated with a series of laboratory experiments and then applied to the Nove reach of the Brenta River (Northern Italy) to assess its bed material transport, interpret channel response to a series of intensive flood events (R.I. ≈ 10 years) and provide a possible evolutionary scenario for the medium term. The study reach is 1400 m long with a mean slope of 0.0039 m m^?1. High‐resolution digital terrain models were produced combining LiDAR data with colour bathymetry techniques. Extensive field sedimentological surveys were also conducted for surface and subsurface material. Data were uploaded in the model and the passage of two consecutive high intensity floods was simulated. The model was run under several hypotheses of sediment supply: one considering substantial equilibrium between sediment input and transport capacity, and the others reducing the sediment supply. The sediment supply was then calibrated comparing channel morphological changes as observed in the field and calculated by the model. Annual bed material transport was assessed and compared with other techniques. Low‐frequency floods (R.I. ≈ 1.5 years) are expected to produce negligible changes in the channel while high floods may erode banks rather than further incising the channel bed. Location and distribution of erosion and deposition areas within the Nove reach were predicted with acceptable biases stemming from imperfections of the model and the specified initial, boundary and forcing conditions. A medium‐term evolutionary scenario simulation underlined the different response to and impact of a consecutive sequence of floods. Copyright © 2015 John Wiley & Sons, Ltd.     

On the sensitivity of the sedimentary system in the East Frisian Wadden Sea to sea-level rise and wave-induced bed shear stress

The paper addresses the individual and collective contribution of different forcing factors (tides, wind waves, and sea-level rise) to the dynamics of sediment in coastal areas. The results are obtained from simulations with the General Estuarine Transport Model coupled with a sediment transport model. The wave-induced bed shear stress is formulated using a simple model based on the concept that the turbulent kinetic energy (TKE) associated with wind waves is a function of orbital velocity, the latter depending on the wave height and water depth. A theory is presented explaining the controls of sediment dynamics by the TKE produced by tides and wind waves. Several scenarios were developed aiming at revealing possible trends resulting from realistic (observed or expected) changes in sea level and wave magnitude. The simulations demonstrate that these changes not only influence the concentration of sediment, which is very sensitive to the magnitude of the external forcing, but also the temporal variability patterns. The joint effect of tides and wave-induced bed shear stress revealed by the comparison between theoretical results and simulations is well pronounced. The intercomparison between different scenarios demonstrates that the spatial patterns of erosion and deposition are very sensitive to the magnitude of wind waves and sea-level rise. Under a changing climate, forcing the horizontal distribution of sediments adjusts mainly through a change in the balance of export and import of sediment from the intertidal basins. The strongest signal associated with this adjustment is simulated North of the barrier islands where the evolution of sedimentation gives an integrated picture of the processes in tidal basins.     

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…     

A one‐dimensional model for simulating armouring and erosion on hillslopes: 1. model development and event‐scale dynamics

This paper presents an erosion model, ARMOUR, which simulates time‐varying runoff, erosion, deposition and surface armour evolution down a hillslope either as a result of a single erosion event or as the cumulative impact of many events over periods up to decades. ARMOUR simulates sediment transport for both cohesive and non‐cohesive soil and dynamically differentiates between ‘transport‐limited’ and ‘source‐limited’ processes. A variety of feasible processes for entrainment of different size classes can be modelled and evaluated against data. The generalized likelihood of uncertainty estimation (GLUE) technique was used to calibrate and validate ARMOUR using data collected during rainfall simulator experiments at two contrasting sites: (1) non‐cohesive stony sediments at Ranger Uranium Mine, Northern Territory, Australia; and (2) cohesive silty sediments at Northparkes Gold Mine, NSW, Australia. The spatial and temporal variations of model predictions within the individual runoff events showed that some entrainment processes could not model the spikes in concentration and subsequent depletion, while the hiding model of Andrews and Parker best simulated the concentration trends for both calibrated and independent runoff events. ARMOUR also successfully captured the coarsening of the surface material, though small, over the duration of the rainfall simulator trials. This was driven by the depletion of the finest size class of the soil. For a constant discharge, ARMOUR simulated higher sediment flux at the start of the storm with the sediment flux and concentration diminishing with time. For natural rainfall a power law relationship between sediment flux and discharge was observed. The calibration exercise showed that sediment concentration and discharge alone are insufficient to calibrate all aspects of the physics, in particular the armour depth. This appears to be because the armouring during the short duration events is driven by depletion of the finest classes of the sediments (diameters less then 62·5 mm), which are not normally measured. Copyright © 2006 John Wiley & Sons, Ltd.