Search

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.

Water level records at two stations in the Guadalquivir Estuary (Spain), one near the estuary mouth (Bonanza) and one about 77 km upstream (Sevilla), have been analysed to study the amplification of the tide in the estuary. The tidal amplification factor shows interesting temporal variation, including a spring-neap variation, some extreme low values, and especially the anomalous behaviour that the amplification factor is larger during a number of periods. These variations are explained by data analysis combined with numerical and analytical modelling. The spring-neap variation is due to the quadratic relation between the bottom friction and the tidal flow velocity. The river flood events are the direct causes of the extreme low values of the amplification factor, and they trigger the non-linear interaction between the tidal flow and suspended sediment transport. The fluvial sediment input during a river flood causes high sediment concentration in the estuary, up to more than 10 g/l. This causes a reduction of the effective hydraulic drag, resulting in stronger tidal amplification in the estuary for a period after a river flood. After such an event the tidal amplification in the estuary does not always fall back to the same level as before the event, indicating that river flood events have significant influence on the long-term development of this estuary.

This is Part II of two papers on man-induced regime shifts in small, narrow, and converging estuaries, with focus on the interaction between effective hydraulic drag, fine sediment import, and tidal amplification, induced by river engineering works, e.g., narrowing and deepening. Paper I describes a simple linear analytical model for the tidal movement in narrow, converging estuaries and a conceptual model on the response of tidal rivers to river engineering works. It is argued that such engineering works may set in motion a snowball effect bringing the river into an alternative steady state. Part II analyses the historic development in tidal range in four rivers, e.g., the Elbe, Ems, Loire, and Scheldt, all in northwest Europe; data are available for many decades, up to a century. We use the analytical model derived in Part I, showing that the effective hydraulic drag in the Ems and Loire has decreased considerably over time, as anticipated in Part I. We did not find evidence that the Upper Sea Scheldt is close to its tipping point towards hyperturbid conditions, but risks have been identified. In the Elbe, tidal reflections against the profound step in bed level around Hamburg seem to have affected the tidal evolution in the last decades. It is emphasized that the conceptual picture sketched in these papers is still hypothetical and needs to be validated, for instance through hind-cast modeling of the evolution of these rivers. This will not be an easy task, as historical data for a proper calibration of the models required are scarce.

This is Part I of two papers on man-induced regime shifts in small, narrow, and converging estuaries, with focus on the interaction between effective hydraulic drag, fine sediment import, and tidal amplification, induced by river engineering works, e.g., narrowing and deepening. In this part, a simple linear analytical model is derived, solving the linearized shallow water equations in exponentially converging tidal rivers. Distinguishing reflecting and non-reflecting conditions, a non-dimensional dispersion equation is derived which yields the real and imaginary wave numbers as a function of the estuarine convergence number and effective hydraulic drag. The estuarine convergence number describes the major geometrical features of a tidal river, e.g., intertidal area, convergence length, and water depth. This model is used in Part II analyzing the historical development of the tide in four rivers. Part I also presents a conceptual model on the response of tidal rivers to narrowing and deepening. It is argued that, upon the loss of intertidal area, flood-dominant conditions prevail, upon which fine sediments are pumped into the river, reducing its effective hydraulic drag. Then a snowball effect may be initiated, bringing the river into a hyper-turbid state. This state is self-maintaining because of entrainment processes, and favorable from an energetic point of view, and therefore highly stable. We may refer to an alternative steady state.

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.

Recent acoustic Doppler current profiler (ADCP)-measurements in the Scheldt estuary near Antwerp, Belgium, revealed anomalous, i.e. anti-clockwise circulations in a left bend during the major part of the flood period; these circulations were established shortly after the turn of the tide. During ebb, anti-clockwise circulations persisted, as predicted by classical theory. These data were analysed with a 3D and a 1DV-model. The 3D simulations reveal that the anomalous circulations are found when salinity is included in the computations—without salinity “normal” circulations were found. From analytical and 1DV simulations, it is concluded that a longitudinal salinity gradient % MathType!Translator!2!1!AMS LaTeX.tdl!TeX -- AMS-LaTeX! % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbbjxAHX % garmWu51MyVXgatuuDJXwAK1uy0HwmaeHbfv3ySLgzG0uy0Hgip5wz % aebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY-Hhbbf9v8qqaq % Fr0xc9pk0xbba9q8WqFfea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qq % Q8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeWaeaaakeaada % WcgaqaaiabgkGi2kaadofaaeaacqGHciITcaWG4baaaaaa!3C47! $${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-\nulldelimiterspace} {\partial x}$$ may induce a near-bed maximum in flow velocity reversing the direction of the secondary currents. The 1DV-model was then used to assess the contribution of various processes one by one. It was found that because of a reduction in vertical mixing, the vertical velocity profile is not at equilibrium during the first phase of accelerating tide, further enhancing the effects of % MathType!Translator!2!1!AMS LaTeX.tdl!TeX -- AMS-LaTeX!% MathType!MTEF!2!1!+-% feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbbjxAHX% garmWu51MyVXgatuuDJXwAK1uy0HwmaeHbfv3ySLgzG0uy0Hgip5wz% aebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY-Hhbbf9v8qqaq% Fr0xc9pk0xbba9q8WqFfea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qq% Q8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeWaeaaakeaada% WcgaqaaiabgkGi2kaadofaaeaacqGHciITcaWG4baaaaaa!3C47! $${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-\nulldelimiterspace} {\partial x}$$ . A small vertical salinity gradient % MathType!Translator!2!1!AMS LaTeX.tdl!TeX -- AMS-LaTeX! % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbbjxAHX % garmWu51MyVXgatuuDJXwAK1uy0HwmaeHbfv3ySLgzG0uy0Hgip5wz % aebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY-Hhbbf9v8qqaq % Fr0xc9pk0xbba9q8WqFfea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qq % Q8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeWaeaaakeaada % WcgaqaaiabgkGi2kaadofaaeaacqGHciITcaWG6baaaaaa!3C49! $${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial z}}} \right. \kern-\nulldelimiterspace} {\partial z}$$ appeared to have a very large effect as the crosscurrents of the secondary circulations induced by % MathType!Translator!2!1!AMS LaTeX.tdl!TeX -- AMS-LaTeX! % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbbjxAHX % garmWu51MyVXgatuuDJXwAK1uy0HwmaeHbfv3ySLgzG0uy0Hgip5wz % aebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY-Hhbbf9v8qqaq % Fr0xc9pk0xbba9q8WqFfea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qq % Q8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeWaeaaakeaada % WcgaqaaiabgkGi2kaadofaaeaacqGHciITcaWG4baaaaaa!3C47! $${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-\nulldelimiterspace} {\partial x}$$ became an order of larger magnitude. However, at the site under consideration, the effects of transverse salinity gradients, generated by differential advection in the river bend, were dominant: adverse directions of the secondary circulations were found even when the vertical velocity profile became more regular with a more or less logarithmic shape, i.e. when the effects of % MathType!Translator!2!1!AMS LaTeX.tdl!TeX -- AMS-LaTeX! % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbbjxAHX % garmWu51MyVXgatuuDJXwAK1uy0HwmaeHbfv3ySLgzG0uy0Hgip5wz % aebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY-Hhbbf9v8qqaq % Fr0xc9pk0xbba9q8WqFfea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qq % Q8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeWaeaaakeaada % WcgaqaaiabgkGi2kaadofaaeaacqGHciITcaWG4baaaaaa!3C47! $${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-\nulldelimiterspace} {\partial x}$$ and % MathType!Translator!2!1!AMS LaTeX.tdl!TeX -- AMS-LaTeX! % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbbjxAHX % garmWu51MyVXgatuuDJXwAK1uy0HwmaeHbfv3ySLgzG0uy0Hgip5wz % aebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY-Hhbbf9v8qqaq % Fr0xc9pk0xbba9q8WqFfea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qq % Q8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeWaeaaakeaada % WcgaqaaiabgkGi2kaadofaaeaacqGHciITcaWG6baaaaaa!3C49! $${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial z}}} \right. \kern-\nulldelimiterspace} {\partial z}$$ did not play a dominant role anymore. It is argued that data on the secondary velocity structure, which can be measured easily owing to today’s developments in ADCP equipment, may serve as an indicator for the accuracy at which the salinity field is computed with 3D numerical models. Moreover, the large effect of the salinity structure on the velocity field must have a large impact on the morphological development of estuaries, which should therefore be accounted for in morphological modelling studies.

The objective of the study presented in this paper is to investigate the predictive capabilities of a process-based sand–mud model in a quantitative way. This recently developed sand–mud model bridges the gap between noncohesive sand models and cohesive mud models. It explicitly takes into account the interaction between these two sediment fractions and temporal and spatial bed composition changes in the sediment bed [Van Ledden (2002) 5:577–594, Van Ledden et al. (2004a) 24:1–11, Van Ledden et al. (2004b) 54:385–391]. The application of this model to idealized situations has demonstrated a good qualitative agreement between observed and computed bed levels and bed composition developments. However, in real-life situations, a realistic quantitative prediction of the magnitude and timescale of this response is important to assess the short-term and long-term impacts of human interventions and/or natural changes. For this purpose, the Friesche Zeegat in the Wadden Sea (the Netherlands) is used as a reference to hindcast the morphological response in the period 1970–1994. Due to the closure of the Lauwerszee in 1969, the tidal prism of this tidal basin was reduced by about 30%. Significant changes in the bed level and bed composition have occurred in the decades following the closure to adjust to the new hydrodynamic conditions. We modeled the long-term bed level and bed composition development in the Friesche Zeegat in the period 1970–1994 starting with the geometry of 1970 by using a research version of Delft3D, which incorporates the sand–mud formulations proposed by [Van Ledden (2002) 5:577–594].The computed total net deposition in the tidal basin in the period 1970–1994 agrees well with the observations, but the observed decrease of the import rate with time is not predicted. The model predicts net deposition in the deeper parts and at the intertidal area in the basin and net erosion in between, which resembles the observations qualitatively. Furthermore, the computed distribution of sand and mud in the basin of the Friesche Zeegat appears to be realistic. Analysis of the results shows that the absence of the decreasing import rate in the basin is caused by a poor quantitative prediction of the changes in the hypsometry of the basin. Because of this, the computed velocity asymmetry in the main channel tends toward flood dominance, whereas the observations indicate that the system is ebb-dominant in 1992. Although the sand–mud model needs to be further improved and verified, the results presented in this paper indicate that the model can be applied as a first step to estimate the effects of human interventions on the large-scale bed level and bed composition changes in tidal systems with sand and mud.

The process-based morphodynamic model Delft3D-MOR and the idealized model of schuttelaars and De Swart (2000) are compared with each other. The differences between the two models in their mathematical-physical formulation as well as the boundary conditions are identified. Their effect on producing cross-sectionally averaged morphological equilibria of tidal inlets with arbitrary length and forced at the seaward boundary by a prescribed M 2 and M 4 sea- surface elevation is studied and an inventory is made of all relevant differences. The physical formulations in the source code of Delft3D-MOR are modified in various steps to resemble the formulations in the idealized model. The effect of each of the differences between the idealized and process-based model are studied by comparing the results of the idealized model to those of the adapted process-based model. The results of the idealized model can be qualitatively reproduced by the process-based model as long as the same morphological boundary condition is applied at the open sea end. This means that the simplifications concerning the mathematical formulation of the physical processes in the idealized model can be justified. Furthermore, it can be inferred that the morphological boundary condition at the open sea end is an essential element in controlling the behaviour of morphodynamic models for tidal inlets and estuaries.

The influence of sand and mud transport on the morphological behaviour of a short tidal basin is investigated in this paper. For this purpose, a morphological model is applied in which sand and mud transport are included and the temporal and spatial bed composition variations are taken into account. Initially, the morphological development shows a sand wave near the entrance of the basin and a mud deposition wave more landward. A quasi equilibrium bed level profile is found after a long period (order century) with a sandy bed surface over almost the entire basin and only a small muddy area near the landward end. The dimensionless ratio between the deposition and erosion flux turns out to be a crucial parameter for the understanding of the observed behaviour. Comparison with previous studies on short tidal basins for sand indicates only that the presence of mud in a combined sand mud model does not change the equilibrium bed level profile considerably for the applied parameter settings herein, but drastically decreases the morphological time scale. Comparison between model results and field data of the Wadden Sea suggests that the obtained bed level and bed composition profile are realistic, indicating that the process-based sand mud model is a first step towards a better understanding of sand mud distributions in tidal basins.