Characterization of bedload intermittency near the threshold of motion using a Lagrangian sediment transport model

At the smallest scales of sediment transport in rivers, the coherent structures of the turbulent boundary layer constitute the fundamental mechanisms of bedload transport, locally increasing the instantaneous hydrodynamic forces acting on sediment particles, and mobilizing them downstream. Near the critical threshold for initiating sediment motion, the interactions of the particles with these unsteady coherent structures and with other sediment grains, produce localized transport events with brief episodes of collective motion occurring due to the near-bed velocity fluctuations. Simulations of these flows pose a significant challenge for numerical models aimed at capturing the physical processes and complex non-linear interactions that generate highly intermittent and self-similar bedload transport fluxes. In this investigation we carry out direct numerical simulations of the flow in a rectangular flat-bed channel, at a Reynolds number equal to Re = 3632, coupled with the discrete element method to simulate the dynamics of spherical particles near the bed. We perform two-way coupled Lagrangian simulations of 48,510 sediment particles, with 4851 fixed particles to account for bed roughness. Our simulations consider a total of eight different values of the non-dimensional Shields parameter to study the evolution of transport statistics. From the trajectory and velocity of each sediment particle, we compute the changes in the probability distribution functions of velocities, bed activity, and jump lengths as the Shields number increases. For the lower shear stresses, the intermittency of the global bedload transport flux is described by computing the singularity or multifr actal spectrum of transport, which also characterizes the widespread range of transport event magnitudes. These findings can help to identify the mechanisms of sediment transport at the particle scale. The statistical analysis can also be used as an ingredient to develop larger, upscaled models for predicting mean transport rates, considering the variability of entrainment and deposition that characterizes the transport near the threshold of motion.     

Sediment processes and flow reversal in the undular tidal bore of the Garonne River (France)

A tidal bore is a series of waves propagating upstream as the tidal flow turns to rising, and the bore front corresponds to the leading edge of the tidal wave in a funnel shaped estuarine zone with macro-tidal conditions. Some field observations were conducted in the tidal bore of the Garonne River on 7 June 2012 in the Arcins channel, a few weeks after a major flood. The tidal bore was a flat undular bore with a Froude number close to unity: $\hbox {Fr}_{1} = 1.02$ and 1.19 (morning and afternoon respectively). A key feature of the study was the simultaneous recording of the water elevation, instantaneous velocity components and suspended sediment concentration (SSC) estimates, together with a detailed characterisation of the sediment bed materials. The sediment was some silty material ( $\hbox {d}_{50} \approx 13~\upmu \hbox {m}$ ) which exhibited some non-Newtonion thixotropic behaviour. The velocity and SSC estimate were recorded simultaneously at high frequency, enabling a quantitative estimate of the suspended sediment flux at the end of the ebb tide and during the early flood tide. The net sediment flux per unit area was directed upstream after the bore, and its magnitude was much larger than that at end of ebb tide. The field observations highlighted a number of unusual features on the morning of 7 June 2012. These included (a) a slight rise in water elevation starting about 70 s prior to the front, (b) a delayed flow reversal about 50 s after the bore front, (c) some large fluctuations in suspended sediment concentration (SSC) about 100 s after the bore front and (d) a transient water elevation lowering about 10 min after the bore front passage. The measurements of water temperature and salinity showed nearly identical results before and after the tidal bore, with no evidence of saline and thermal front during the study.     

Kohäsive Feinpartikel in fluvialen Systemen: einige Gedanken zu Erkenntnissen und Forschungsdefiziten

Fluvial cohesive sediment is of fundamental environmental and multidisciplinary concern owing to its significant impact on energy, nutrient and trace element fluxes, river sediment budgets as well as habitat quality. Consequently, numerous studies in geomorphology, hydraulics, hydrology and river ecology accentuate the importance of fluvial fine grained sediment. Increased deposition of fine grained particles for instance can negatively effect benthic habitats. Settling velocity and deposition rates of fines are key terms in sediment transport models. Thus, the knowledge of transport and quantitative deposition dynamics are necessary prerequisites for a sustainable sediment and habitat management as well as a reliable, physically based modelling of fine sediment conveyance and contaminant transfer in streams. However, it is generally a challenging task to characterise and trace fine-grained particle transport, deposition and resuspension dynamics in a fluvial environment. One inherent complexity of fluvial cohesive sediment dynamics is for instance the relevance of interparticle forces and flocculation/aggregation processes. Owing to this complexity and the lack of standardized methods cohesive particle transport mechanisms are hardly predictable in a specific field situation. Consequently, interdisciplinary research approaches are needed. This article discusses specific findings and research gaps in the field of cohesive sediment research.     

Bed load transport of sediments and morphodynamics in the Topolobampo coastal lagoon system,Mexico

The Topolobampo coastal lagoon system, located on the eastern side of the Gulf of California, is a marine zone of considerable economic importance with vessel traffic, dredging operations and aquaculture development. Despite its relevance as a conservation site, this ecosystem has been poorly studied. Since life in marine substrates is abundant, we investigated the capability of tidal hydrodynamics in the lagoon to erode and to accumulate sediment. We calculated the morphodynamics caused by bed load sediment transport applying a two dimensional non-linear hydrodynamic finite difference model. Bed erosion and accretion patterns of sediment (for specific grain size: 170 μm) were determined from the divergence of sediment transport. After a year of numerical simulation of sediment transport the areas of noticeable changes on the bottom of the lagoon have been revealed. Most of sediment accretion took place in the narrow steeped channel connecting the Topolobampo and Ohuira sections. This area appears characterized by the presence of high tidal velocity gradients. Another finding was that accretion areas were coupled with erosion zones in an alternating form. This outcome suggests that sediment do not travel long distances but is deposited nearby the erosion sites. The results are strong evidence of the influence of tidal hydrodynamics on the sediments distribution in the Topolobampo coastal lagoon system and on the generation of substrates where marine life may find protection.     

An engineering model for countercurrent flow under wind-induced waves and current

A general model for the phase-averaged velocity field in wind-induced countercurrent flow is proposed. The influence of waves on the time-averaged velocity is accounted for by introducing a skewness factor in a parabolic eddy viscosity model. The skewness factor represents the net effect of the wavy surface in the engineering model for velocity. The coherent velocity components are described separately by an orbital velocity obtained from linear wave theory and are added to the time-averaged components to give a complete model for the phase-averaged velocity field. The proposed model collapses to the standard model for deep-water conditions, but is also shown to yield the correct behavior for intermediate conditions. Moreover, the bed shear stress, derived from the proposed velocity model, is also shown to be in agreement with experiments.     

A model diagnostic study of age of river-borne sediment transport in the tidal York River Estuary

As nutrients and organic matters are transported preferentially in an adsorbed state and tend to bind to the sediments, sediment transport plays an important role on eutrophication processes in the estuaries. The timescale of sediment transport is of significance for studying the retention of pollutants and eutrophication processes in the estuaries. Unlike transport of dissolved substances that is mainly controlled by advection and diffusion processes, the sediment transport is significantly affected by the intermittent settling and resuspension processes. A three-dimensional model with suspended sediment transport was utilized to investigate the transport timescale of river-borne sediment in the tidal York River Estuary. The results indicate that river discharge dominantly determines the age of river-borne sediment in the estuary. High river discharge results in a low sediment age compared to that under mean flow. The intermittent effects of settling and resuspension events greatly affect the river-borne sediment age. Both settling velocity and critical shear stress are shown to be key parameters in determining the sediment transport timescale. The sediment age decreases as settling velocity and/or critical shear stress decrease, while it increases with the increase of settling velocity that prevents the sediment to be transported out of the estuary.     

One-dimensional numerical modelling of dam-break waves over movable beds: application to experimental and field cases

This paper reports a numerical study on dam-break waves over movable beds. A one-dimensional (1-D) model is built upon the Saint-Venant equations for shallow water waves, the Exner equation of sediment mass conservation and a spatial lag equation for non-equilibrium sediment transport. The set of governing equations is solved using an explicit finite difference scheme. The model is tested in various idealized experimental cases, with fairly good agreement between the numerical predictions and measurements. Discrepancies are observed at the earlier stage of the dam-break wave and around the dam location due to no vertical velocity component being taken into account. Sensitivity tests confirm that the friction coefficient is an important parameter for the evaluation of sediment transport processes operating during a dam-break wave. The influence of the non-equilibrium adaptation length (or the lag distance) is negligible on the wavefront celerity and weak on the free surface and bed profiles, which indicates that one may ignore the spatial lag effect in dam-break wave studies. Finally, the simulation of the Lake Ha!Ha! dyke-break flood event shows that the model can provide relevant results if a convenient formula for computing the sediment transport capacity and an appropriate median grain diameter of riverbed material are selected.     

Pier scour within long contraction in cohesive sediment bed

In the present study, an attempt has been made to study the bridge pier scour embedded long contraction for the bed sediment prepared by a mixture of clay and fine sand with varying proportions having a specific range of antecedent moisture content. Results particularly focused on the clay–sand mixed cohesive bed at varying clay fractions of the sediment bed, approach flow velocity, contraction ratio and different pier shapes, on maximum equilibrium scour depth for pier scour within long contraction. Further, regression based equations for the estimation of non-dimensional maximum scour depth for piers within long contraction in clay–sand mixed cohesive bed embedded were proposed.     

Environmental Impact of Undular Tidal Bores in Tropical Rivers

A tidal bore impacts significantly on the estuarine ecosystem, although little is known on the flow field, mixing and sediment motion beneath tidal bores. In the absence of detailed systematic field measurements, a quasi-steady flow analogy was applied to investigate undular tidal bores with inflow Froude numbers between 1.25 and 1.6. Experimental results indicated that rapid flow redistributions occur beneath the free-surface undulations, with significant variations in bed shear stress between wave crests and troughs. Dynamic similarity was used to predict detailed flow characteristics of undular tidal bores. The effects of periodic loading on river sediments, scour of river bed and flow mixing behind the bore are discussed. A better understanding of these processes will contribute to better management practices in tidal bore affected rivers, including the Styx and Daly rivers in tropical Australia.     

Large Eddy Simulation of turbulence generated by a weak breaking tidal bore

A tidal bore is a natural and fragile phenomenon, which is of great importance for the ecology of an estuary. The bore development is closely linked with the tidal range and the river mouth shape, and its existence is sensitive to any small change in boundary conditions. Despite their ecological and cultural value, little is known on the flow field, turbulent mixing and sediment motion beneath tidal bores. Indeed, some striking features can be highlighted in two-dimensional simulations, such as large velocity fluctuations and flow recirculation structures. Using Large Eddy Simulation method, the numerical results emphasised the complicated turbulent structures and their unsteadiness under a tidal bore.     

A neural network experiment on the simulation of daily nitrate-nitrogen and suspended sediment fluxes from a small agricultural catchment

In this paper, we report an application of neural networks to simulate daily nitrate-nitrogen and suspended sediment fluxes from a small 7.1 km² agricultural catchment (Melarchez), 70 km east of Paris, France. Nitrate-nitrogen and sediment losses are only a few possible consequences of soil erosion and biochemical applications associated to human activities such as intensive agriculture. Stacked multilayer perceptrons models (MLPs) like the ones explored here are based on commonly available inputs and yet are reasonably accurate considering their simplicity and ease of implementation. Note that the simulation does not resort on water quality flux observations at previous time steps as model inputs, which would be appropriate, for example, to predict the water chemistry of a drinking water plant a few time steps ahead. The water quality fluxes are strictly mapped to historical mean flux values and to hydro-climatic variables such as stream flow, rainfall, and soil moisture index (12 model input candidates in total), allowing its usage even when no flux observations are available. Self-organizing feature maps based on the network structure established by Kohonen were employed first to produce the training and the testing data sets, with the intent to produce statistically close subsets so that any difference in model performance between validation and testing has to be associated to the model and not to the data subsets. The stacked MLPs reached different levels of performance simulating the nitrate-nitrogen flux and the suspended sediment flux. In the first instance, 2-input stacked MLP nitrate-nitrogen simulations, based on the same-day stream flow and on the 80-cm soil moisture index, have a performance of almost 90% according to the efficiency index. On the other hand, the performance of 3-input stacked MLPs (same-day stream flow, same-day historical flux, and same-day stream flow increment) reached a little more than 75% according to the same criterion. The results presented here are deemed already promising enough, and should encourage water resources managers to implement simple models whenever appropriate.     

Mixing at a sediment concentration interface in turbulent open channel flow

In this work we address the role of turbulence on mixing of clear layer of fluid with sediment-laden layer of fluid at a sediment concentration interface. This process can be conceived as the entrainment of sediment-free fluid into the sediment-laden layer, or alternatively, as the transport of sediment into the top sediment-free flow. This process is governed by four parameters—Reynolds number of the flow \(Re_\tau\), non-dimensional settling velocity of the sediment (proxy for sediment size) \(\tilde{V}\), Richardson number \(Ri_\tau\) and Schmidt number Sc. For this work we have performed direct numerical simulations for fixed Reynolds and Schmidt numbers while varying the values of Richardson number and particle settling velocity. In the simple model considered here, the flow’s momentum and turbulence pre-exists over the entire layer of fluid, while the sediment is initially confined to a layer close to the bed. Mixing of sediment-free fluid with the sediment-laden layer is associated primarily with upward transport of sediment and buoyancy. There is no simultaneous upward transport of fluid momentum and turbulence into the sediment-free fluid layer, which is already in motion and turbulent. The analysis performed shows that the ability of the flow to transport a given sediment size decreases with the distance from the bottom, and thus only fine enough sediment particles are transported across the sediment concentration interface. For these cases, the concentration profiles evolve to a final steady state in good agreement with the well-known Rouse profile. The approach towards the Rouse profile happens through a transient self-similar state. This behavior of the flow is not seen for larger particles. Detailed analysis of the three dimensional structure of the sediment concentration interface shows the mechanisms by which sediment particles are lifted up by tongues of sediment-laden fluid with positive correlation between vertical velocity and sediment concentration. Finally, the mixing ability of the flow is addressed by monitoring the time evolution of the center of mass of the sediment-laden layer and the vertical location of the sediment-free/sediment-laden interface.     

A Constant Erosion Rate Model for Fine Sediment in the York River,Virginia

Results of field experiments using the Virginia Institute of Marine Science' Sea Carousel and tripod system reveal a highly dynamic sediment activity at the Clay Bank site in the York River. At the water-sediment interface, the critical bed shear stress for sediment erosion varied between 0.026 Pa and 0.1 Pa. For the well consolidated sediment below the interface, the critical bed shear stress increased significantly, to more than 0.6 Pa. The seasonal variation of erosion rate for the surficial sediment is significant. For the well consolidated sediment below the surficial sediment, however, the seasonal variation diminished and erosion rate approached a constant. Experimental results from the Carousel deployments imply that the erosion process for quasi-steady tidal flows is always near equilibrium. For this reason, a constant rate erosion model is proposed for the time during tidal acceleration phases. For other times, the erosion rate would be zero. The measured suspended sediment concentration at 10 cm above bed from the VIMS tripod system also Supports this model.     

A flume study of drift in marine infaunal amphipods (Haustoriidae)

Amphipods of the infaunal family Haustoriidae are characteristic of high-energy marine sands and occur both in the sediment and the overlying water column. Sediments in this habitat are subject to constant reworking by tidal currents, suggesting that resident amphipod populations are affected by this disturbance in a phenomenon similar to freshwater invertebrate drift. A controlled-velocity laboratory flume was used to examine the effect of haustoriid density, current velocity, illumination, and food availability on drift rates to determine a causal basis for drift. Drift is densityindependent and greatest at night and during high current flow; it is also greater at night from sterile sediment than from untreated sand. Flume transport was usually less than 10% of amphipods present in the sediment. Haustoriids captured downstream were mostly adults occurring in a 1:1 sex ratio, suggesting no obvious function of drift in reproduction. Current-induced displacement of haustoriids may produce the patchiness in distribution observed in nature. Disturbance of the bed could also function to keep amphipod densities below the carrying capacity of the local environment. In certain cases, food limitation may cause amphipods to actively leave the substrate. Under all conditions, greater drift rates in darkness are probably adaptive in avoidance of predators. Despite the nature of sediment movement in a high-energy environment, haustoriid drift may have an active as well as passive component.     

Fourier analysis of the roll-up and merging of coherent structures in shallow mixing layers

The current study investigates the role of nonlinearity in the development of two-dimensional coherent structures (2DCS) in shallow mixing layers. A nonlinear numerical model based on the depth-averaged shallow water equations is used to investigate temporal shallow mixing layers, where the mapping from temporal to spatial results is made using the velocity at the center of the mixing layer. The flow is periodic in the stream-wise direction and the transmissive boundary conditions are used in the cross-stream boundaries to prevent reflections. The numerical results are examined with the aid of Fourier decomposition. Results show that the previous success in applying local linear theory to shallow mixing layers does not imply that the flow is truly linear. Linear stability theory is confirmed to be only valid within a short distance from the inflow boundary. Downstream of this linear region, nonlinearity becomes important for the roll-up and merging of 2DCS. While the energy required for the merging of 2DCS is still largely provided by the velocity shear, the merging mechanism is one where nonlinear mode interaction changes the velocity field of the subharmonic mode and the gradient of the along-stream velocity profile which, in turn, changes the magnitude of the energy production of the subharmonic mode by the velocity shear implicitly. The nonlinear mode interaction is associated with energy up-scaling and is consistent with the inverse energy cascade which is expected to occur in shallow shear flows. Current results also show that such implicit nonlinear interaction is sensitive to the phase angle difference between the most unstable mode and its subharmonic. The bed friction effect on the 2DCS is relatively small initially and grows in tandem with the size of the 2DCS. The bed friction also causes a decrease in the velocity gradient as the flow develops downstream. The transition from unstable to stable flow occurs when the bed friction balances the energy production. Beyond this point, the bed friction is more dominant and the 2DCS are progressively damped and eventually get annihilated. The energy production by the velocity shear plays an important role from the upstream end all the way to the point of transition to stable flow. The fact that linear stability theory is valid only for a short distance from the inflow boundary suggests that some elements of nonlinearity is incorporated in the mean velocity profile in experiments by the averaging process. The implicit nature of nonlinear interaction in shallow mixing layers and the sensitivity of the nonlinear interaction to phase angle difference between the most unstable mode and its subharmonic allows local linear theory to be successful in reproducing features of the instability such as the dominant mode of the 2DCS and its amplitude.     

The wave-induced solute flux from submerged sediment

The issue of the transport of dissolved nutrients and contaminants between the sediment in the bottom of a lake or reservoir and the body of water above it is an important one for many reasons. In particular the biological and chemical condition of the body of water is intricately linked to these mass transport processes. As the review by Boudreau (Rev Geophys 38(3):389–416, 2000) clearly demonstrates those transport processes are very complex involving mechanisms as diverse as the wave-induced flux between the sediment and the overlying water and the effect of burrowing animals on the transport within the sediment as well as basic diffusion mechanisms. The present paper focuses on one facet of these transport processes; we re-examine the balance of diffusion and wave-induced advection and demonstrate that the wave-induced flux of a solute from submerged sediment is not necessarily purely diffusive as suggested by Harrison et al. (J Geophys Res 88:7617–7622, 1983) but can be dominated by a mean or time-averaged flux induced by the advective fluid motion into and out of the sediment caused by the fluctuating pressure waves associated with wave motion. Indeed along the subtidal shoreline where the fluctuating bottom pressures are greatest, wave-induced advection will dominate the mean, time-averaged transport of solute into or out of the sediment as suggested in the work of Riedl et al. (Mar Biol 13:210–221, 1972). However, the present calculations also indicate that this advective flux decreases rapidly with increasing depth so that further away from the shoreline the advective flux becomes negligible relative to the diffusive flux and therefore the latter dominates in deeper water.     

Model of turbulence penetration into a suspension layer on a sediment bed and effect on vertical solute transfer

The vertical diffusional mass (solute) transfer through a suspended sediment layer, e.g. at the bottom of a lake, reservoir or estuary, by the propagation of velocity fluctuations from above was investigated. The attenuation of the velocity fluctuations in the suspension layer and the associated effect on solute transfer through the suspension layer was simulated. To represent large eddies traveling downstream in water over a high-concentration suspended sediment layer, a streamwise velocity fluctuation moving in downstream direction was imposed along the upper boundary of the suspension layer. Velocity fluctuations and downstream velocity were normalized by the shearvelocity (U_*) at the top of the suspension layer. Streamwise and vertical velocity components inside the suspension layer, were obtained from the 2-D continuity and the Navier–Stokes equations. The persistence of turbulence with depth—as it penetrates from the overlying water into the suspension layer—was found to depend on its amplitude, its period, and on the apparent viscosity of the suspension. The turbulence was found to propagate efficiently into the suspension layer when its frequency is low, and the apparent viscosity of the suspension is high. Effects on vertical mass transfer were parameterized by penetration depth and effective diffusion coefficient, and related to apparent viscosity of the suspension, Schmidt number and shear velocity on top of the suspension layer. The enhancement of turbulence penetration by viscosity is similar to the flow near an oscillating flat plate (Stokes’ second problem), but is opposite to turbulence penetration into a stationary porous and permeable sediment bed. The information is applicable to water quality modeling mear the sediment/water interface of lakes, river impoundments and estuaries.     

Instantaneous transport of a passive scalar in a turbulent separated flow

The results of large-eddy simulations of flow and transient solute transport over a backward facing step and through a 180° bend are presented. The simulations are validated successfully in terms of hydrodynamics and tracer transport with experimental velocity data and measured residence time distribution curves confirming the accuracy of the method. The hydrodynamics are characterised by flow separation and subsequent recirculation in vertical and horizontal directions and the solute dispersion process is a direct response to the significant unsteadiness and turbulence in the flow. The turbulence in the system is analysed and quantified in terms of power density spectra and covariance of velocity fluctuations. The injection of an instantaneous passive tracer and its dispersion through the system is simulated. Large-eddy simulations enable the resolution of the instantaneous flow field and it is demonstrated that the instabilities of intermittent large-scale structures play a distinguished role in the solute transport. The advection and diffusion of the scalar is governed by the severe unsteadiness of the flow and this is visualised and quantified. The analysis of the scalar mass transport budget quantifies the mechanisms controlling the turbulent mixing and reveals that the mass flux is dominated by advection.     

Estimation of total sediment load concentration obtained by experimental study using artificial neural networks

Estimation of sediment concentration in rivers is very important for water resources projects planning and managements. The sediment concentration is generally determined from the direct measurement of sediment concentration of river or from sediment transport equations. Direct measurement is very expensive and cannot be conducted for all river gauge stations. However, sediment transport equations do not agree with each other and require many detailed data on the flow and sediment characteristics. The main purpose of the study is to establish an effective model which includes nonlinear relations between dependent (total sediment load concentration) and independent (bed slope, flow discharge, and sediment particle size) variables. In the present study, by performing 60 experiments for various independent data, dependent variables were obtained, because of the complexity of the phenomena, as a soft computing method artificial neural networks (ANNs) which is the powerful tool for input–output mapping is used. However, ANN model was compared with total sediment transport equations. The results show that ANN model is found to be significantly superior to total sediment transport equations.     

Estimation of methane emitted by natural wetlands of Orissa: Based on case studies on three representative wetlands

The rapid increase in the concentration of greenhouse gases in our atmosphere has led to significant changes in global climate. Methane is second to carbon dioxide in terms of its contribution to global warming. The wetlands serve as a major natural source of methane. Orissa State, located on the east coast of India, has a number of natural wetlands. Three representative wetlands: Chilka – the largest semi-saline lake of Asia, Anshupa – a land-locked fresh water lake and Gahirmatha – a tidal flood plain of Bhitara Kanika region that has sprawling mangrove vegetation (extension of the Sunderbans of the Gangetic delta) were chosen for intensive monitoring of water and sediment quality, along with methane flux over a period of three years during 1997–2000. The average seasonal methane fluxes were integrated to arrive at the annual flux from each wetland category. Finally, those results were extended to other natural wetlands of Orissa to arrive at the average annual methane flux of Orissa State. This is useful in calculating the total budget of greenhouse gases of India.