The impact of submersed aquatic vegetation on the development of river mouth bars

This work is inspired by the sudden resurgence of the submersed aquatic vegetation (SAV) bed in the Chesapeake Bay (USA). Because the SAV bed occurs at the mouth of the Bay's main tributary (Susquehanna River), it plays a significant role in modulating sediment and nutrient inputs from the Susquehanna to the Bay. Previous model studies on the impact of submersed aquatic vegetation on the development of river mouth bars lacked a complete mechanistic understanding. This study takes advantage of new advances in 3D computational models that include explicit physical-sedimentological feedbacks to obtain this understanding. Specifically, we used Delft3D, a state-of-the-art hydrodynamic model that provides fine-scale computations of three-dimensional flow velocity and bed shear stress, which can be linked to sediment deposition and erosion. Vegetation is modeled using a parameterization of hydraulic roughness that depends on vegetation height, stem density, diameter, and drag coefficient. We evaluate the hydrodynamics, bed shear stresses, and sediment dynamics for different vegetation scenarios under conditions of low and high river discharge. Model runs vary the vegetation height, density, river discharge, and suspended-sediment concentration. Numerical results from the idealized model show that dense SAV on river mouth bars substantially diverts river discharge into adjacent channels and promotes sediment deposition at ridge margins, as well as upstream bar migration. Increasing vegetation height and density forms sandier bars closer to the river mouth and alteration of the bar shape. Thus, this study highlights the important role of SAV in shaping estuarine geomorphology, which is especially relevant for coastal management. © 2019 John Wiley & Sons, Ltd.     

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.     

A 1D model for tides waves and fine sediment in short tidal basins—Application to the Wadden Sea

In order to simulate the dynamics of fine sediments in short tidal basins, like the Wadden Sea basins, a 1D cross-sectional averaged model is constructed to simulate tidal flow, depth-limited waves, and fine sediment transport. The key for this 1D model lies in the definition of the geometry (width and depth as function of the streamwise coordinate). The geometry is computed by implementing the water level and flow data, from a 2D flow simulation, and the hypsometric curve in the continuity equation. By means of a finite volume method, the shallow-water equations and sediment transport equations are solved. The bed shear stress consists of the sum of shear stresses by waves and flow, in which the waves are computed with a depth-limited growth equation for wave height and wave frequency. A new formulation for erosion of fines from a sandy bed is proposed in the transport equation for fine sediment. It is shown by comparison with 2D simulations and field measurements that a 1D schematization gives a proper representation of the dynamics in short tidal basins.

     

Development of a river delta: a case study of Cimanuk river mouth,Indonesia

This paper describes delta development processes with particular reference to Cimanuk Delta in Indonesia. Cimanuk river delta, the most rapidly growing river delta in Indonesia, is located on the northern coast of Java Island. The delta is subject to ocean waves of less than 1 m height due to its position in the semi‐enclosed Java Sea in the Indonesian archipelago. The study has been carried out using a hydrodynamic model that accounts for sediment movement through the rivers and estuaries. As an advanced approach to management of river deltas, a numerical model, namely MIKE‐21, is used as a tool in the management of Cimanuk river delta. From calibration and verification of hydrodynamic model, it was found that the best value of bed roughness was 0·1 m. For the sediment‐transport model, the calibration parameters were adjusted to obtain the most satisfactory results of suspended sediment concentration and volume of deposition. By comparing the computed and observed data in the calibration, the best values of critical bed shear stress for deposition, critical bed shear stress for erosion and erosion coefficient were 0·05 N m^?2, 0·15 N m^?2, and 0·00001 kg m^?2 s^?1, respectively. The calibrated model was then used to analyse sensitivity of model parameters and to simulate delta development during the periods 1945–1963 and 1981–1997. It was found that the sensitive model parameters were bed shear stresses for deposition and erosion, while the important model inputs were river suspended sediment concentration, sediment characteristics and hydrodynamic. The model result showed reasonable agreement with the observed data. As evidenced by field data, the mathematical model proves that the Cimanuk river delta is a river‐dominated delta because of its protrusion pattern and very high sediment loads from the Cimanuk river. It was concluded that 86% of sediment load from the Cimanuk river was deposited in the Cimanuk delta. Copyright © 2008 John Wiley & Sons, Ltd.     

A 1D model for tides waves and fine sediment in short tidal basins—Application to the Wadden Sea

In order to simulate the dynamics of fine sediments in short tidal basins, like the Wadden Sea basins, a 1D cross-sectional averaged model is constructed to simulate tidal flow, depth-limited waves, and fine sediment transport. The key for this 1D model lies in the definition of the geometry (width and depth as function of the streamwise coordinate). The geometry is computed by implementing the water level and flow data, from a 2D flow simulation, and the hypsometric curve in the continuity equation. By means of a finite volume method, the shallow-water equations and sediment transport equations are solved. The bed shear stress consists of the sum of shear stresses by waves and flow, in which the waves are computed with a depth-limited growth equation for wave height and wave frequency. A new formulation for erosion of fines from a sandy bed is proposed in the transport equation for fine sediment. It is shown by comparison with 2D simulations and field measurements that a 1D schematization gives a proper representation of the dynamics in short tidal basins.     

River ice cover influence on sediment transportation at present and under projected hydroclimatic conditions

A large number of rivers are frozen annually, and the river ice cover has an influence on the geomorphological processes. These processes in cohesive sediment rivers are not fully understood. Therefore, this paper demonstrates the impact of river ice cover on sediment transport, i.e. turbidity, suspended sediment loads and erosion potential, compared with a river with ice‐free flow conditions. The present sediment transportation conditions during the annual cycle are analysed, and the implications of climate change on wintertime geomorphological processes are estimated. A one‐dimensional hydrodynamic model has been applied to the Kokemäenjoki River in Southwest Finland. The shear stress forces directed to the river bed are simulated with present and projected hydroclimatic conditions. The results of shear stress simulations indicate that a thermally formed smooth ice cover diminishes river bed erosion, compared with an ice‐free river with similar discharges. Based on long‐term field data, the river ice cover reduces turbidity statistically significantly. Furthermore, suspended sediment concentrations measured in ice‐free and ice‐covered river water reveal a diminishing effect of ice cover on riverine sediment load. The hydrodynamic simulations suggest that the influence of rippled ice cover on shear stress is varying. Climate change is projected to increase the winter discharges by 27–77% on average by 2070–2099. Thus, the increasing winter discharges and possible diminishing ice cover periods both increase the erosion potential of the river bed. Hence, the wintertime sediment load of the river is expected to become larger in the future. Copyright © 2015 John Wiley & Sons, Ltd.     

Size sorting of fine-grained sediments during erosion: Results from the western Gulf of Lions

Sediment cores from the western Gulf of Lions France were subject to known bottom shear stresses with the goal of understanding size-specific sediment erodibility. On cruises in October 2004, February and April 2005, cores with an undisturbed sediment–water interface were collected along a transect extending seaward from the Tet river mouth. The cores were exposed to increasing shear stresses (0.01–0.4 Pa) onboard the vessel shortly after collection by using a Gust erosion chamber. Samples of the suspensate were collected during the erosion experiments and analyzed for disaggregated inorganic grain size (DIGS) using a Coulter Multisizer IIe. Size-specific mobility plots were generated by dividing the proportion of each grain size in suspension at each shear stress by its proportion in the sediment before erosion. If all grain sizes that make up the bottom sediment are eroded equally from the bed, then mobility equals one for all grain sizes. Values >1 indicate that the suspended sediment is enriched in the size class and values <1 indicate that the size class is enriched in the bed. Results show that in non-cohesive, sandy silts, fine grains (clays and fine silts) are eroded preferentially from the bed at low shear stresses. With increasing bottom stress progressively larger grains are eroded from the bed. In cohesive silts, preferential erosion of the finer sizes no longer occurs, with all sizes up to medium silts eroding at approximately the same rate. Effectively, a sandy silt can be winnowed of its fine grain fraction during erosion while cohesive silts cannot. This difference in the sortability of cohesive and non-cohesive sediment during erosion may control the position and maintenance of the sand–mud transition and the sequestration of surface-adsorbed contaminants.     

A probabilistic model for sediment entrainment:The role of bed irregularity

A generalized probabilistic model is developed in this study to predict sediment entrainment under the incipient motion, rolling, and pickup modes. A novelty of the proposed model is that it incorporates in its formulation the probability density function of the bed shear stress, instead of the near-bed velocity fluctuations, to account for the effects of both flow turbulence and bed surface irregularity on sediment entrainment. The proposed model incorporates in its formulation the collective effects of three para-meters describing bed surface irregularity, namely the relative roughness, the volumetric fraction and relative position of sediment particles within the active layer. Another key feature of the model is that it provides a criterion for estimating the lift and drag coefficients jointly based on the recognition that lift and drag forces acting on sediment particles are interdependent and vary with particle protrusion and packing density. The model was validated using laboratory data of both fine and coarse sediment and was compared with previously published models. The study results show that all the examined models perform adequately for the fine sediment data, where the sediment particles have more uniform gra-dation and relative roughness is not a factor. The proposed model was particularly suited for the coarse sediment data, where the increased bed irregularity was captured by the new parameters introduced in the model formulation. As a result, the proposed model yielded smaller prediction errors and physically acceptable values for the lift coefficient compared to the other models in case of the coarse sediment data.     

THREE DIMENSIONAL NUMERICAL MODELLING OF FLOW AND SEDIMENT TRANSPORT IN RIVERS

The 3D numerical model, ECOMSED (open source code), was used to simulate flow and sediment transport in rivers. The model has a long history of successful applications to oceanic, coastal and estuarine waters. Improvements in the advection scheme, treatment of river roughness parameterization and shear stress partitioning were necessary to reproduce realistic and comparable results in a river application. To account for the dynamics of the mobile bed boundary, a model for the bed load transport was included in the code. The model reproduced observed secondary currents, bed shear stress distribution and erosion-deposition patterns on a curved channel. The model also successfully predicted the general flow patterns and sediment transport characteristics of a 1-km long reach of the River Klar?lven, located in the north of the county of V?rmland, Sweden.     

PROCESS-BASED MODEL OF GRAIN LIFTING FROM RIVER BED TO ESTIMATE SUSPENDED-SEDIMENT CONCENTRATION IN A SMALL HEADWATER BASIN

Suspended sediment is supplied from river bed sediment in Hiyamizusawa Brook, Hokkaido, Japan, during the early snowmelt season. The stirring up of fine grains from the river bed is an important control of the time variation of suspended-sediment flux. In this stream, about 10 per cent of the river bed is covered with sand sediment, 80 per cent with cobbles and/or pebbles and the remaining 10 per cent is exposed bedrock. A model previously used to explain the stirring up of fine grains within a cobble and pebble bed is applied to a sand bed, with the modification that fine grains in a sand bed are assumed to be stirred up from the tractive layer formed on the surface, whereas those in a cobble and pebble bed are assumed to be stirred up from the gaps formed by the selective movement of pebbles on the river bed. The lift force acting at the river bed is estimated from the bed shear stress, and the maximum grain size capable of being stirred up was calculated from the lift force. Consequently, the amount of fine material stirred up from the river bed is estimated from the grain size distribution of river bed sediment, and the suspended-sediment flux is thus calculated. All stirred-up fines are assumed to become suspended sediment. The simulated time variation of suspended-sediment concentration was similar to that obtained in the field study. The calculated grain size of suspended sediment was also equivalent to the field data.     

Numerical study of generation of the tidal shear front off the Yellow River mouth

Tidal shear front off the Yellow River mouth has been observed and modeled in the previous studies. However, a detailed investigation of the front generation has not been conducted. The aim of this paper is to use a three-dimensional tidal model coupled to a sediment transport module to examine the front formation. The model predicted a tidal shear front that propagated offshore and lasted 1–2 h at both flood and ebb phase off the Yellow River mouth. The sensitivity numerical experiments showed that the topography with a strong slope off the Yellow River mouth was a determining factor for the front generation, and a parallel orientation between the major axes of ellipses and co-tidal lines of maximum tidal current was a necessary condition. While the bottom friction and the river runoff had no effect on the front location but affected the front intensity, the front generation was not sensitive to the coastline variation. The study concluded that the bottom slope off the river mouth induces a strong variation in the bottom stress in a cross-shore direction, which produces both maximum phase gradient and sediment concentration variability across the tidal shear front. With the extending Yellow River delta, the tidal shear front under the new bathymetry of year 2003 has been strengthened and pushed further offshore due to an increased bottom slope.     

Sediment mobility and bed armoring in the St Clair River: insights from hydrodynamic modeling

The lake levels in Lake Michigan‐Huron have recently fallen to near historical lows, as has the elevation difference between Lake Michigan‐Huron compared to Lake Erie. This decline in lake levels has the potential to cause detrimental impacts on the lake ecosystems, together with social and economic impacts on communities in the entire Great Lakes region. Results from past work suggest that morphological changes in the St Clair River, which is the only natural outlet for Lake Michigan‐Huron, could be an appreciable factor in the recent trends of lake level decline. A key research question is whether bed erosion within the river has caused an increase in water conveyance, therefore, contributed to the falling lake level. In this paper, a numerical modeling approach with field data is used to investigate the possibility of sediment movement in the St Clair River and assess the likelihood of morphological change under the current flow regime. A two‐dimensional numerical model was used to study flow structure, bed shear stress, and sediment mobility/armoring over a range of flow discharges. Boundary conditions for the numerical model were provided by detailed field measurements that included high‐resolution bathymetry and three‐dimensional flow velocities. The results indicate that, without considering other effects, under the current range of flow conditions, the shear stresses produced by the river flow are too low to transport most of the coarse bed sediment within the reach and are too low to cause substantial bed erosion or bed scour. However, the detailed maps of the bed show mobile bedforms in the upper St Clair River that are indicative of sediment transport. Relatively high shear stresses near a constriction at the upstream end of the river and at channel bends could cause local scour and deposition. Ship‐induced propeller wake erosion also is a likely cause of sediment movement in the entire reach. Other factors that may promote sediment movement, such as ice cover and dredging in the lower river, require further investigation. Copyright © 2012 John Wiley & Sons, Ltd.     

Formation process of meandering channel by a 2D numerical simulation

A 2D depth-averaged model for hydrodynamic,sediment transport and river morphological adjustment is presented.The sediment transport submodel considers non-uniform sediment,bed surface armoring,impact of secondary flow on the direction of bed-load transport,and transverse slope of river bed.The bank erosion submodel incorporates a simple simulation method for updating bank geometry during either degradational or aggradational bed evolution.The model is applied to a 180°bend with a constant radius under unsteady flow conditions,and to Friedkin’s laboratory meander channels.The results are in acceptable agreement with measurements,confirming the two dimensional model’s potential in predicting the formation of river meandering and improving understanding of patterning processes.Future researches are needed to clarify some simplifications and limitations of the model.     

Relationship between bed shear stress and suspended sediment concentration:annular flume experiments

In this study,annular flume experiments were carried out,using the sediment samples collected from the lower part of the inter-tidal zone at Xiaoyangkou,Jiangsu coast,China.The Ariathurai-Partheniades equation was used to determine the bed shear stress,by evaluating variations in the suspended sediment concentration within the water column.The derived relation between the bed shear stress and suspended sediment concentration shows that,at various stages of seabed erosion, suspended sediment concentration increases rapidly when the flow velocity is increased,but the pattern of change in the bed shear stress does not follow suit.At low concentrations,bed shear stress initially increases markedly with increasing flow velocity.However,when the concentration reaches an apparently critical level around 0.55 kg m"3,the rate of change in the bed shear stress abruptly slows down,or becomes almost constant,in response to further increases in the flow velocity.Results of experiments indicate that,from a critical level onward,suspended sediment concentration has a strong influence on the bed shear stress.     

Effective macroscopic transport parameters between parallel plates with constant concentration boundaries

A macroscopic transport model is developed, following the Taylor shear dispersion analysis procedure, for a 2D laminar shear flow between parallel plates possessing a constant specified concentration. This idealized geometry models flow with contaminant dissolution at pore-scale in a contaminant source zone and flow in a rock fracture with dissolving walls. We upscale a macroscopic transient transport model with effective transport coefficients of mean velocity, macroscopic dispersion, and first-order mass transfer rate. To validate the macroscopic model the mean concentration, covariance, and wall concentration gradient are compared to the results of numerical simulations of the advection–diffusion equation and the Graetz solution. Results indicate that in the presence of local-scale variations and constant concentration boundaries, the upscaled mean velocity and macrodispersion coefficient differ from those of the Taylor–Aris dispersion, and the mass transfer flux described by the first-order mass transfer model is larger than the diffusive mass flux from the constant wall. In addition, the upscaled first-order mass transfer coefficient in the macroscopic model depends only on the plate gap and diffusion coefficient. Therefore, the upscaled first-order mass transfer coefficient is independent of the mean velocity and travel distance, leading to a constant pore-scale Sherwood number of 12. By contrast, the effective Sherwood number determined by the diffusive mass flux is a function of the Peclet number for small Peclet number, and approaches a constant of 10.3 for large Peclet number.     

Modelling of meander migration in an incised channel

An updated linear computer model for meandering rivers with incision has been developed. The model simulates the bed topography, flow field, and bank erosion rate in an incised meandering channel. In a scenario where the upstream sediment load decreases （e.g., after dam closure or soil conservation）, alluvial river experiences cross section deepening and slope flattening. The channel migration rate might be affected in two ways： decreased channel slope and steeped bank height. The proposed numerical model combines the traditional one-dimensional （1D） sediment transport model in simulating the channel erosion and the linear model for channel meandering. A non-equilibrium sediment transport model is used to update the channel bed elevation and gradations. A linear meandering model was used to calculate the channel alignment and bank erosion/accretion, which in turn was used by the 1D sediment transport model. In the 1D sediment transport model, the channel bed elevation and gradations are represented in each channel cross section. In the meandering model, the bed elevation and gradations are stored in two dimensional （2D） cells to represent the channel and terrain properties （elevation and gradation）. A new method is proposed to exchange information regarding bed elevations and bed material fractions between 1D river geometry and 2D channel and terrain. The ability of the model is demonstrated using the simulation of the laboratory channel migration of Friedkin in which channel incision occurs at the upstream end.     

On the use of horizontal acoustic Doppler profilers for continuous bed shear stress monitoring

Continuous monitoring of bed shear stress in large river systems may serve to better estimate alluvial sediment transport to the coastal ocean.Here we explore the possibility of using a horizontally deployed acoustic Doppler current profiler(ADCP) to monitor bed shear stress,applying a prescribed boundary layer model,previously used for discharge estimation.The model parameters include the local roughness length and a dip correction factor to account for sidewall effects.Both these parameters depend on river stage and on the position in the cross-section, and were estimated from shipborne ADCP data.We applied the calibrated boundary layer model to obtain bed shear stress estimates over the measuring range of the HADCP.To validate the results,co-located coupled ADCPs were used to infer bed shear stress,both from Reynolds stress profiles and from mean velocity profiles. From HADCP data collected over a period of 1.5 years,a time series of width profiles of bed shear stress was obtained for a tidal reach of the Mahakam River,East Kalimantan,Indonesia.A smaller dataset covering 25 hours was used for comparison with results from the coupled ADCPs.The bed shear stress estimates derived from Reynolds stress profiles appeared to be strongly affected by local effects causing upflow and downflow,which are not included in the boundary layer model used to derive bed shear stress with the horizontal ADCP.Bed shear stresses from the coupled ADCP are representative of a much more localized flow,while those derived with the horizontal ADCP resemble the net effect of the flow over larger scales.Bed shear stresses obtained from mean velocity profiles from the coupled ADCPs show a good agreement between the two methods,and highlight the robustness of the method to uncertainty in the estimates of the roughness length.