Net sediment transport in tidal basins: quantifying the tidal barotropic mechanisms in a unified framework

Net sediment transport in tidal basins is a subtle imbalance between large fluxes produced by the flood/ebb alternation. The imbalance arises from several mechanisms of suspended transport. Lag effects and tidal asymmetries are regarded as dominant, but defined in different frames of reference (Lagrangian and Eulerian, respectively). A quantitative ranking of their effectiveness is therefore missing. Furthermore, although wind waves are recognized as crucial for tidal flats’ morphodynamics, a systematic analysis of the interaction with tidal mechanisms has not been carried out so far. We review the tide-induced barotropic mechanisms and discuss the shortcomings of their current classification for numerical process-based models. Hence, we conceive a unified Eulerian framework accounting for wave-induced resuspension. A new methodology is proposed to decompose the sediment fluxes accordingly, which is applicable without needing (semi-) analytical approximations. The approach is tested with a one-dimensional model of the Vlie basin, Wadden Sea (The Netherlands). Results show that lag-driven transport is dominant for the finer fractions (silt and mud). In absence of waves, net sediment fluxes are landward and spatial (advective) lag effects are dominant. In presence of waves, sediment can be exported from the tidal flats and temporal (local) lag effects are dominant. Conversely, sand transport is dominated by the asymmetry of peak ebb/flood velocities. We show that the direction of lag-driven transport can be estimated by the gradient of hydrodynamic energy. In agreement with previous studies, our results support the conceptualization of tidal flats’ equilibrium as a simplified balance between tidal mechanisms and wave resuspension.     

The links between entrance geometry,hypsometry and hydrodynamics in shallow tidally dominated basins

Geomorphological characteristics of tidal basins control hydrodynamics and sediment transport potential within such basins, for example, by adjusting the balance in tidal asymmetry. In this study we examine the effects of entrance geometry on tidal velocity asymmetry, slack water asymmetry, bed shear stress patterns and hypsometric profile shapes by comparison of six shallow meso-tidal basins of Tauranga Harbour, New Zealand. Numerical model results show how tidal distortion increases with distance from a basin entrance. A simple ratio between basin width and entrance width defines levels of basin dilation. Sub-basins with a constricted geometry and deep entrance channels are associated with small bed shear stress values and high rates of flood-directed tidal velocity asymmetry in the sheltered basin centres, indicating a large potential for sediment deposition of larger particles. Moreover, slack water asymmetry within these basins is weakly ebb-directed, indicating a small potential for transport of fine sediments out of the basins. The constricted depositional basins are characterized by convex hypsometric profiles with elevated intertidal regions. Unconstricted geometries are associated with larger bed shear stress values and more ebb-directed tidal velocity asymmetry within basin centres, suggesting limited potential for overall sediment deposition. The slack tide duration asymmetry is weakly flood-dominant indicating that limited input of fine sediment into the basins is possible. The comparatively high-energy conditions within these exposed basins are associated with a less convex hypsometric intertidal profile. The ability to estimate tidal asymmetries is advantageous when developing management strategies related to ecosystem functioning, navigability or coastal protection in specific geomorphic settings. © 2019 John Wiley & Sons, Ltd.     

Parameterization of wind turbine impacts on hydrodynamics and sediment transport

Monopile foundations of offshore wind turbines modify the hydrodynamics and sediment transport at local and regional scales. The aim of this work is to assess these modifications and to parameterize them in a regional model. In the present study, this is achieved through a regional circulation model, coupled with a sediment transport module, using two approaches. One approach is to explicitly model the monopiles in the mesh as dry cells, and the other is to parameterize them by adding a drag force term to the momentum and turbulence equations. Idealised cases are run using hydrodynamical conditions and sediment grain sizes typical from the area located off Courseulles-sur-Mer (Normandy, France), where an offshore windfarm is under planning, to assess the capacity of the model to reproduce the effect of the monopile on the environment. Then, the model is applied to a real configuration on an area including the future offshore windfarm of Courseulles-sur-Mer. Four monopiles are represented in the model using both approaches, and modifications of the hydrodynamics and sediment transport are assessed over a tidal cycle. In relation to local hydrodynamic effects, it is observed that currents increase at the side of the monopile and decrease in front of and downstream of the monopile. In relation to sediment transport effect, the results show that resuspension and erosion occur around the monopile in locations where the current speed increases due to the monopile presence, and sediments deposit downstream where the bed shear stress is lower. During the tidal cycle, wakes downstream of the monopile reach the following monopile and modify the velocity magnitude and suspended sediment concentration patterns around the second monopile.     

The effect of tidal asymmetry and temporal settling lag on sediment trapping in tidal estuaries

Over decades and centuries, the mean depth of estuaries changes due to sea-level rise, land subsidence, infilling, and dredging projects. These processes produce changes in relative roughness (friction) and mixing, resulting in fundamental changes in the characteristics of the horizontal (velocity) and vertical tides (sea surface elevation) and the dynamics of sediment trapping. To investigate such changes, a 2DV model is developed. The model equations consist of the width-averaged shallow water equations and a sediment balance equation. Together with the condition of morphodynamic equilibrium, these equations are solved analytically by making a regular expansion of the various physical variables in a small parameter. Using these analytic solutions, we are able to gain insight into the fundamental physical processes resulting in sediment trapping in an estuary by studying various forcings separately. As a case study, we consider the Ems estuary. Between 1980 and 2005, successive deepening of the Ems estuary has significantly altered the tidal and sediment dynamics. The tidal range and the surface sediment concentration has increased and the position of the turbidity zone has shifted into the freshwater zone. The model is used to determine the causes of these historical changes. It is found that the increase of the tidal amplitude toward the end of the embayment is the combined effect of the deepening of the estuary and a 37% and 50% reduction in the vertical eddy viscosity and stress parameter, respectively. The physical mechanism resulting in the trapping of sediment, the number of trapping regions, and their sensitivity to grain size are explained by careful analysis of the various contributions of the residual sediment transport. It is found that sediment is trapped in the estuary by a delicate balance between the M ₂ transport and the residual transport for fine sediment ( $\emph{w}_s=0.2$  mm s^???1) and the residual, M ₂ and M ₄ transports for coarser sediment ( $\emph{w}_s=2$  mm s^???1). The upstream movement of the estuarine turbidity maximum into the freshwater zone in 2005 is mainly the result of changes in tidal asymmetry. Moreover, the difference between the sediment distribution for different grain sizes in the same year can be attributed to changes in the temporal settling lag.     

Influence of topography on tide propagation and amplification in semi-enclosed basins

An idealized model for tide propagation and amplification in semi-enclosed rectangular basins is presented, accounting for depth differences by a combination of longitudinal and lateral topographic steps. The basin geometry is formed by several adjacent compartments of identical width, each having either a uniform depth or two depths separated by a transverse topographic step. The problem is forced by an incoming Kelvin wave at the open end, while allowing waves to radiate outward. The solution in each compartment is written as the superposition of (semi)-analytical wave solutions in an infinite channel, individually satisfying the depth-averaged linear shallow water equations on the f plane, including bottom friction. A collocation technique is employed to satisfy continuity of elevation and flux across the longitudinal topographic steps between the compartments. The model results show that the tidal wave in shallow parts displays slower propagation, enhanced dissipation and amplified amplitudes. This reveals a resonance mechanism, occurring when the length of the shallow end is roughly an odd multiple of the quarter Kelvin wavelength. Alternatively, for sufficiently wide basins, also Poincaré waves may become resonant. A transverse step implies different wavelengths of the incoming and reflected Kelvin wave, leading to increased amplitudes in shallow regions and a shift of amphidromic points in the direction of the deeper part. Including the shallow parts near the basin’s closed end (thus capturing the Kelvin resonance mechanism) is essential to reproduce semi-diurnal and diurnal tide observations in the Gulf of California, the Adriatic Sea and the Persian Gulf.     

Tidal and meteorological forcing of sediment transport in tributary mudflat channels

Field observations of flow and sediment transport in a tributary channel through intertidal mudflats indicate that suspended sediment was closely linked to advection and dispersion of a tidal salinity front. During calm weather when tidal forcing was dominant, high concentrations of suspended sediment advected up the mudflat channel in the narrow region between salty water from San Francisco Bay and much fresher runoff from the small local watershed. Salinity and suspended sediment dispersed at similar rates through each tidal inundation, such that during receding ebbs the sediment pulse had spread spatially and maximum concentrations had decreased. Net sediment transport was moderately onshore during the calm weather, as asymmetries in stratification due to tidal straining of the salinity front enhanced deposition, particularly during weaker neap tidal forcing. Sediment transport by tidal forcing was periodically altered by winter storms. During storms, strong winds from the south generated wind waves and temporarily increased suspended sediment concentrations. Increased discharge down the tributary channels due to precipitation had more lasting impact on sediment transport, supplying both buoyancy and fine sediment to the system. Net sediment transport depended on the balance between calm weather tidal forcing and perturbations by episodic storms. Net transport in the tributary channel was generally off-shore during storms and during calm weather spring tides, and on-shore during calm weather neap tides.     

The influence of changes in tidal asymmetry on residual sediment transport in the Western Scheldt

This study fits into a wider research program from RWS RIKZ concerning the exchange of sediment between the coast and the tidal basins, Western Scheldt and Wadden Sea, which are the largest basins along the Dutch Coast, over different time-scales. For both basins, questions about the evolution of the import/export at the mouth recently arose. In case of the Western Scheldt, which is the subject of this study, mainly the uncertainty about the future developments after the change from import to export at the mouth was noticed in the 1990s, which necessitated a more detailed study of this area.     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

Suspended sediment transport along an idealised tidal embayment: settling lag, residual transport and the interpretation of tidal signals

We present semi-analytical solutions for suspended sediment concentration (SSC) and residual sediment transport in a simple mathematical model of a short tidal embayment. These solutions allow us to investigate in some detail the characteristic tidal and semi-tidal variation of SSC and the processes leading to residual sediment transport, including settling and scour lags, the roles of ‘local’ and ‘advective’ contributions, and the presence of internally or externally generated overtides. By interpreting the transport mechanisms in terms of the classic conceptual models of settling lag we clarify how these models may be expressed in mathematical terms. Our results suggest that settling lag is usually a more important process than scour lag, and that a local model which neglects advection may predict the direction of net sediment transport incorrectly. Finally, we discuss our results in the context of other transport processes and morphodynamic feedback.     

Effect of tidal forcing on a subterranean estuary

Groundwater flow and chemical transport in subterranean estuaries are poorly understood despite their potentially important implications for chemical fluxes from aquifers to coastal waters. Here, a numerical study of the dynamics in a subterranean estuary subject to tidal forcing is presented. Simulations show that salt transport associated with tidally driven seawater recirculation leads to the formation of an upper saline plume in the intertidal region. Computed transit times and flow velocities indicate that this plume represents a more active zone for mixing and reaction than the dispersion zone of the lower, classical salt wedge. Proper conceptualisation of this surficial mixing zone extends our understanding of processes within the subterranean estuary. Numerical tracer simulations reveal that tidal forcing may reduce the threat of a land-derived contaminant discharging to the marine environment by modifying the subsurface transport pathway and local geochemical conditions. Mixing and stratification in the subterranean estuary are strongly affected by both inland and tidal forcing. Based on the estuarine analogy we present a systematic classification of subterranean estuaries.     

Morphological controls in sandy estuaries: the influence of tidal flats and bathymetry on sediment transport

The morphodynamics of shallow, vertically well-mixed estuaries, characterised by tidal flats and deeper channels, have been investigated. This paper examines what contributes to flood/ebb-dominant sediment transport in localised regions through a 2D model study (using the TELEMAC modelling system). The Dyfi Estuary in Wales, UK has been used as a case study and, together with idealised estuary shapes, shows that shallow water depths lead to flood dominance in the inner estuary whilst tidal flats and deep channels cause ebb dominance in the outer estuary. For medium sands and with an artificially ‘flattened’ bathymetry (i.e. no tidal flats), the net sediment transport switches from ebb-dominant to flood-dominant where the parameter a/h (local tidal amplitude ÷ local tidally averaged water depth) exceeds 1.2. Sea level rise will reduce this critical value of a/h and also reduce the ebb-directed sediment transport significantly, leading to a flood-dominated estuarine system. A similar pattern, albeit with greater transport, was simulated with tidal flats included and also with a reduced grain size. This suggests that analogous classifications for flood/ebb asymmetry of the tide in estuaries as a whole may not represent the local sediment transport in sufficient detail. Through the Dyfi simulations, the above criterion involving a/h is shown to be complicated further by augmented flow past a spit at the estuary mouth which gives rise to a self-maintaining scour hole. Simulations of one year of bed evolution in an idealised flat-bottomed estuary, including tidal flow past a spit, recreate the flood/ebb dominance on either side of the spit and the formation of a scour hole in between. The erosion rate at the centre of the hole is reduced as the hole deepens, suggesting the establishment of a self-maintaining equilibrium state.     

Suspended sediment fluxes at an intertidal flat: The shifting influence of wave,wind, tidal,and freshwater forcing

Using in situ, continuous, high frequency (8–16 Hz) measurements of velocity, suspended sediment concentration (SSC), and salinity, we investigate the factors affecting near-bed sediment flux during and after a meteorological event (cold front) on an intertidal flat in central San Francisco Bay. Hydrodynamic forcing occurs over many frequency bands including wind wave, ocean swell, seiching (500–1000 s), tidal, and infra-tidal frequencies, and varies greatly over the time scale of hours and days. Sediment fluxes occur primarily due to variations in flow and SSC at three different scales: residual (tidally averaged), tidal, and seiching. During the meteorological event, sediment fluxes are dominated by increases in tidally averaged SSC and flow. Runoff and wind-induced circulation contribute to an order of magnitude increase in tidally averaged offshore flow, while waves and seiching motions from wind forcing cause an order of magnitude increase in tidally averaged SSC. Sediment fluxes during calm periods are dominated by asymmetries in SSC over a tidal cycle. Freshwater forcing produces sharp salinity fronts which trap sediment and sweep by the sensors over short (∼30 min) time scales, and occur primarily during the flood. The resulting flood dominance in SSC is magnified or reversed by variations in wind forcing between the flood and ebb. Long-term records show that more than half of wind events (sustained speeds of greater than 5 m/s) occur for 3 h or less, suggesting that asymmetric wind forcing over a tidal cycle commonly occurs. Seiching associated with wind and its variation produces onshore sediment transport. Overall, the changing hydrodynamic and meteorological forcing influence sediment flux at both short (minutes) and long (days) time scales.     

Wave-induced steady streaming,mass transport and net sediment transport in rough turbulent ocean bottom boundary layers

The interaction between two important mechanisms which causes streaming has been investigated by numerical simulations of the seabed boundary layer beneath both sinusoidal waves and Stokes second order waves, as well as horizontally uniform bottom boundary layers with asymmetric forcing. These two mechanisms are streaming caused by turbulence asymmetry in successive wave half-cycles (beneath asymmetric forcing), and streaming caused by the presence of a vertical wave velocity within the seabed boundary layer as earlier explained by Longuet-Higgins. The effect of wave asymmetry, wave length to water depth ratio, and bottom roughness have been investigated for realistic physical situations. The streaming induced sediment dynamics near the ocean bottom has been investigated; both the resulting suspended load and bedload are presented. Finally, the mass transport (wave-averaged Lagrangian velocity) has been studied for a range of wave conditions. The streaming velocities beneath sinusoidal waves (Longuet-Higgins streaming) is always in the direction of wave propagation, while the streaming velocities in horizontally uniform boundary layers with asymmetric forcing are always negative. Thus the effect of asymmetry in second order Stokes waves is either to reduce the streaming velocity in the direction of wave propagation, or, for long waves relative to the water depth, to induce a streaming velocity against the direction of wave propagation. It appears that the Longuet-Higgins streaming decreases as the wave length increases for a given water depth, and the effect of wave asymmetry can dominate, leading to a steady streaming against the wave propagation. Furthermore, the asymmetry of second order Stokes waves reduces the mass transport (wave-averaged Lagrangian velocity) as compared with sinusoidal waves. The boundary layer streaming leads to a wave-averaged transport of suspended sediments and bedload in the direction of wave propagation.     

Effect of advective and diffusive sediment transport on the formation of local and global bottom patterns in tidal embayments

Bathymetric field data of tidal basins reveal two main classes of bottom patterns: (1) tidal bars, located near the entrance of the basin (length scale determined by the embayment width) and (2) global channel-shoal patterns which scale with the basin length. Previous models were able to describe only either one of these patterns. In this paper it is shown that both of them can be investigated within the framework of an idealised model of a rectangular tidal embayment, with fixed side walls and an erodible bed. The water motion is described by the depth-averaged shallow-water equations and is forced by a prescribed vertical tide at the seaward entrance. Sediment is transported as suspended load and only realistic values of the bottom friction parameter are considered. By assuming the ratio of embayment length over tidal wave length to be small, the model allows for a morphodynamic equilibrium, characterised by a spatially uniform tide moving over a bottom which slopes upwards toward the landward boundary. This equilibrium is unstable for a range of values of the model parameters, such that growth of bedforms occurs. Both global and local bottom patterns are found. In this study particular emphasis is laid on the mechanism governing the growth of a new type of localised bottom pattern. These patterns consist of small bars located near the entrance of the basin, resembling multiple row bars, and are found when advective sediment fluxes prevail over diffusive sediment fluxes. The formation process of these new bedforms is discussed in detail. The results agree well with field data. Comparison of the results with those obtained with a process-based, numerical model shows that, although the idealised model is strongly simplified, it is capable of producing the essential morphodynamics. Therefore, the idealised model is a useful tool to investigate mechanisms of bottom pattern growth.Responsible Editor: Iris Grabemann     

Numerical modeling of tidal currents, sediment transport and morphological evolution in Hangzhou Bay, China

A 2D depth-averaged numerical model is set up to simulate the macro-scale hydrodynamic characteristics, sediment transport patterns and morphological evolution in Hangzhou Bay, a large macro-tidal estuary on the eastern coast of China. By incorporating the shallow water equations, the suspended sediment transport equation and the mass-balance equation for sediment; short-term hydrodynamics, sediment transport and long-term morphological evolution for Hangzhou Bay are simulated and the underlying physical mechanisms are analyzed. The model reproduces the spatial distribution patterns of suspended sediment concentration (SSC) in Hangzhou Bay, characterized by three high SSC zones and two low SSC zones. It also correctly simulates the residual flow, the residual sediment transport and the sediment accumulation patterns in Hangzhou Bay. The model results are in agreement with previous studies based on field measurements. The residual flow and the residual sediment transport are landwards directed in the northern part of the bay and seawards directed in the southern part. Sediment accumulation takes place in most areas of the bay. Harmonic analysis revealed that the tide is flood-dominant in the northern part of the bay and ebb-dominant in the southern part of the bay. The strength of the flood-dominance increases landwards along the northern Hangzhou Bay. In turn sediment transport in Hangzhou Bay is controlled by this tidal asymmetry pattern. In addition, the direction of tidal propagation in the East China Sea, the presence of the archipelago in the southeast and the funnel-shaped geometry of the bay, play important roles for the patterns of sediment transport and sediment accumulation respectively.     

The effect of a small vegetation dieback event on salt marsh sediment transport

Vegetation is a critical component of the ecogeomorphic feedbacks that allow a salt marsh to build soil and accrete vertically. Vegetation dieback can therefore have detrimental effects on marsh stability, especially under conditions of rising sea levels. Here, we report a variety of sediment transport measurements associated with an unexpected, natural dieback in a rapidly prograding marsh in the Altamaha River Estuary, Georgia. We find that vegetation mortality led to a significant loss in elevation at the dieback site as evidenced by measurements of vertical accretion, erosion, and surface topography compared to vegetated reference areas. Below-ground vegetation mortality led to reduced soil shear strength. The dieback site displayed an erosional, concave-up topographic profile, in contrast to the reference sites. At the location directly impacted by the dieback, there was a reduction in flood dominance of suspended sediment concentration. Our work illustrates how a vegetation disturbance can at least temporarily reverse the local trajectory of a prograding marsh and produce complex patterns of sediment transport. © 2018 John Wiley & Sons, Ltd.     

On the stability relationships between tidal asymmetry and morphologies of tidal basins and estuaries

Simple stability relationships are practically useful to provide a rapid assessment of coastal and estuarine landforms in response to human interventions and long‐term climate change. In this contribution, we review a variety of simple stability relationships which are based on the analysis of tidal asymmetry (TA). Most of the existing TA‐based stability relationships are derived using the one‐dimensional tidal flow equations assuming a certain regular shape of the tidal channel cross‐sections. To facilitate analytical solutions, specific assumptions inevitably need to be made, for example by linearizing the friction term and dropping some negligible terms in the tidal flow equations. We find that three major types of TA‐based stability relationships have been proposed between three non‐dimensional channel geometric ratios (represented by the ratio of channel widths, ratio of wet surface areas and ratio of storage volumes) and the tide‐related parameter a/h (i.e. the ratio between tidal amplitude and mean water depth). Based on established geometric relations, we use these non‐dimensional ratios to restate the existing relationships so that they are directly comparable. Available datasets are further extended to examine the utility of these TA‐based relationships. Although a certain agreement is shown for these relationships, we also observe a large scatter of data points which are collected in different types of landscape, hydrodynamic and sedimentological settings over the world. We discuss in detail the potential reasons for this large scatter and subsequently elaborate on the limited applicability of the various TA‐based stability relationships for practical use. We highlight the need to delve further into what constitutes equilibrium and what is needed to develop more robust measures to determine the morphological state of these systems. Copyright © 2018 John Wiley & Sons, Ltd.