1 

Residuele stromingen in een getijdenkanaal ten gevolge van de bodemwrijving
Met behulp van een wiskundig model het bepalen van de invloed van de bodemwrijving op de residuele stroming in een open kanaal.

[PDF]
[Abstract]

2 

Simulation of drying and ﬂooding in a tidal embayment using the level set method

[PDF]

3 

Waterbeweging in een estuarium: analytische modellering van getijde beweging voor diepteafhankelijke viscositeitsprofielen met behulp van de breedtegemiddelde ondiepwatervergelijkingen.
In een estuarium is het belangrijk te weten wat de kenmerken zijn van de stroming van het water om de ecologie te bestuderen. Hiervoor wordt een tweedimensionaal model met behulp van de breedtegemiddelde ondiepwatervergelijkingen (shallow water equations) afgeleid. De aandacht wordt daarna vervolgens gefocust op de leidende orde vergelijkingen. Bij het oplossen van deze vergelijkingen wordt de nadruk wordt gelegd op de modellering van de viscositeit in de diepterichting. Er wordt daarbij vooral gekeken naar de invloed van verschillende viscositeitsprofielen op de waterbeweging.

[PDF]
[Abstract]

4 

The Influence of Retention Basins on Tidal Dynamics in Estuaries
The Ems estuary suffers from high concentrations of suspended sediment, which is harmful for the life in the river. Sediment import occurs because of tidal asymmetry, a difference between the velocity of the current during high tide and low tide. The use of retention basins are a possible measure to reduce the tidal asymmetry.
An analytic model, using partial differential equations derived from the Shallow Water equations and NavierStokes equations, was used to investigate the influence of retention basins on the tidal dynamics and sediment transport in estuaries in general and the Ems estuary specifically. For these investigations, this analytic model is useful because system conditions can easily be changed.
Results show that the locations of retention basins are of great importance for their performance. In the Ems estuary the location minimizing the sediment import is at the end of the estuary, near the weir. Enlarging the retention basin is better than placing a second basin at a different location in the estuary.

[PDF]
[Abstract]

5 

Equilibrium and stability of a double inlet system
Barrier island coasts are a common feature in many parts of the world. An example is the Wadden coast of The Netherlands, Germany and Denmark. These coasts consist of barrier islands separated by tidal inlets with at the landward side tidal basins.Characteristic for the Wadden Sea is that the tidal basins are not completely separated,but are connected via topographic highs allowing exchange of water between the basins. As a result the tidal inlets that connect the basins to the North Sea will interact. The focus in this thesis is on the effect of this interaction on the crosssectional equilibrium and stability of tidal inlets that are part of a double inlet system. The knowledge gained in this study will help to develop rational management plans for this kind of system. In determining the equilibrium values and stability of crosssectional areas of the inlets use is made of flow diagrams. A flow diagram consists of the equilibrium flow curves of each inlet and a flow field showing the adaptation of the inlet crosssections after the system has been removed from equilibrium. Each intersection of the equilibrium flow curves represents a stable or unstable equilibrium. The equilibrium flow curve for each
inlet is the locus of the values of the crosssectional areas for which the velocity
amplitude in the inlet equals the equilibrium velocity i.e. approximately 1 m/s according to ESCOFFIER [1940]. As a start the double inlet system is schematized as a basin connected to the ocean by two channels. The water surface area of the basin is assumed constant and water levels are assumed to fluctuate uniformly. On the seaward side a simple harmonic, semidiurnal tide is used to force the system. Analyzing the double inlet system under these conditions by means of the flow diagrams leads to the conclusion that a stable equilibrium of the two inlets does not exist. Ultimately only one inlet remains open and the other will close. This confirms the earlier conclusions of VAN DE KREEKE [1990] and BORSJE [2003] concerning the crosssectional stability of multiple inlet systems.

[PDF]
[Abstract]

6 

Sediment transport patterns in the Dutch Western Wadden Sea: influence of forcing and bathymetry
The Dutch Western Wadden Sea, comprising the basins of Texel, Eierlandse Gat and Vlie, suffers a shortage of sediments. This implies that sediment is imported into the basins. The influence of the main forcing agents on the sediment transport is investigated. These forcing agents are the tide, waves and wind. In case of the wind forcing, both the largescale and local wind forcing are taken into account.
The tide induces an import of sediments through the Texel inlet, and export through the inlets of Eierlandse Gat and Vlie. Waves have most influence on the ebbtidal delta region of the inlets. Due to the sheltering effect of these deltas and the barrier islands they lose most of their energy when they enter the inlets. As a result the influence of the waves on the sediment transports inside the basins is low. Wind forcing generates a residual flow that is of the same order of magnitude as the residual flow generated by the tide. It induces residual sediment transports through the inlets that are of the same order of magnitude as the tide induced sediment transports.
A wave and wind climate is constructed to determine the annual residual sediment transports. With the bathymetry of the year 1998, in total a volume of 6306*103m3/yr is imported into the basins.
With the help of the semiempirical equilibrium relations, the required volume of sediments for the basins to be in equilibrium is determined. In total a volume of 13.78*108m3 is added, 88% of this is placed in the Texel basin. As a result the transport through the Texel inlet approaches the equilibrium (little to no residual transport), but through the other two inlets there is a strong residual export. The shift of the tidal divides is one reason for these results. These divides shift towards the Texel inlet, adjusting the tidal propagation in the basins. In the basins of the Eierlandse Gat and Vlie this results in an increase of the exports. This effect was not taken into account with the determination of the required volume of sediments. Another reason for the strong export through the inlets of the Eierlandse Gat and Vlie is the placement of the added sediments. In the basins of these inlets the additional sediments are predominantly placed in the ebbchannels. This increases the residual export of sediment through these channels.
This implies that the determined equilibrium bathymetry is not the exact equilibrium. Yet it gives a fair indication of the amount of sediments needed for the basins to be in equilibrium.

[PDF]
[Abstract]

7 

Tidal Embayments: Modelling and understanding their morphodynamics
Tidal basins are observed all over the world's coastline, for example along the north coast of The Netherlands. These basins are important both from an economic and ecological point of view. Complex channel and shoal patterns can be found in these inlets, as they develop due to the interaction of the currents generated by tides, wind and density differences with the erodible bed. External conditions, such as human interferences, can also strongly influence the morphologic behaviour. The aim of this thesis is to elucidate the physical mechanisms resulting in the formation of channel and shoal patterns in short tidal embayments.

[PDF]
[Abstract]

8 

Resonance characteristics of tides in branching channels

file embargo until: 20140709

9 

Subsurface Characterization Using a Cellular Automaton Approach
In this paper, a random Cellular Automaton model is developed to characterise heterogeneity of geological formations. The CAmodel is multilateral and can be easily applied in both two and three dimensions. We demonstrate that conditioning on well data is possible and the method is numerically efficient. To construct the model, the subsurface is subdivided into N cells, with an initial lithology assigned to each cell. Rules to update the current cell states are chosen from a set of rules, independently for each cell. The model converges typically in less than 50 iterations to a steady state or periodic solution. Within one period the realisations exhibit similar statistical properties. The final fraction of the various lithologies can be tuned by choosing the proper initial fractions. In this way, geological knowledge of those fractions can be satisfied.

[PDF]
[Abstract]

10 

Impact of estuarine convergence on residual circulation in tidally energetic estuaries and inlets
Estuarine convergence (landward reduction of width and/or depth) is known to have the potential to significantly enhance estuarine circulation, a result theoretically derived under the assumption of constant eddy viscosity. Recent studies of longitudinally uniform energetic tidal channels indicate that tidal straining, a process driven by tidally varying eddy viscosity, is a major driver of estuarine circulation. The combined effect of estuarine convergence and tidal straining is investigated, for the first time, in this paper. The present idealized numerical study shows that estuarine convergence is reducing or even reversing tidal straining circulation in such a way that estuarine circulation can be weakened. This is a counterintuitive hydrodynamic effect of estuarine convergence, which may reduce (rather than increase) upestuary particulate matter transport in estuaries and tidal inlets.

file embargo until: 20140804
[Abstract]

11 

The influence of basin geometry on equilibrium and stability of double inlet systems
This study investigates the influence of basin geometry on the crosssectional stability of double inlet systems. The inlet is in equilibrium when the amplitude of the inlet velocities equals the equilibrium velocity (~1 m s1). This equilibrium is stable when after a perturbation the crosssections of both inlets return to their original equilibrium value. The necessary amplitudes of the inlet velocities are obtained using an idealized 2DH hydrodynamic that calculates tidal elevation and flow in a geometry consisting of several adjacent rectangular compartments.
Model results suggest that regardless of the inclusion or exclusion of bottom friction in the basin, stable equilibrium states exist. Qualitatively, the influence of basin geometry does not change the presence of stable equilibrium. Quantitatively, however, taking a basin surface area of 1200 km2, equilibrium values can differ up to a factor 2 depending on the geometry of the basin.

[PDF]
[Abstract]

12 

Instability of timedependent winddriven ocean gyres
The winddriven ocean circulation at midlatitudes is susceptible to several types of instabilities. One of the simplest models of these flows is the quasigeostrophic barotropic potential vorticity equation in an idealized ocean basin. In this model, the route to complex spatio/temporal flows is through successive bifurcations. The aim of this study is to describe the physics of the destabilization process of a periodic winddriven flow associated with a secondary bifurcation. Although bifurcation theory has proven to be a valuable tool to determine the physical mechanisms of destabilization of fluid flows, the analysis of the stability of timedependent (for example, periodic) flows, using this methodology, is computationally unpractical, due to the large number of degreesoffreedom involved. The approach followed here is to construct a loworder model using numerical Galerkin projection of the full model equations onto the dynamically active eigenmodes. The resulting reduced model is shown to capture the local dynamics of the full model. The physical mechanism of the destabilization of the periodic winddriven flow is deduced from the reduced model. While there are several stabilizing processes, notably rectification, the destabilization occurs due to timedependent increase of the background horizontal shear in the flow.

[PDF]
[Abstract]

13 

Lateral entrapment of sediment in tidal estuaries: An idealized model study
Two physical mechanisms leading to lateral accumulation of sediment in tidally dominated estuaries are investigated, involving Coriolis forcing and lateral density gradients. An idealized model is used that consists of the three‐dimensional shallow water equations and sediment mass balance. Conditions are assumed to be uniform in the along‐estuary direction. A semidiurnal tidal discharge and tidally averaged density gradients are prescribed. The erosional sediment flux at the bed depends both on the bed shear stress and on the amount of sediment available in mud reaches for resuspension. The distribution of mud reaches over the bed is selected such that sediment transport is in morphodynamic equilibrium, that is, tidally averaged erosion and deposition of sediment at the bed balance. Analytical solutions are obtained by using perturbation analysis. Results suggest that in most estuaries lateral density gradients induce more sediment transport than Coriolis forcing. When frictional forces are small (Ekman number E < 0.02), the Coriolis mechanism dominates and accumulates sediment on the right bank (looking up‐estuary in the Northern Hemisphere). On the other hand, when frictional forces are moderate to high (E > 0.02), the lateral density gradient mechanism dominates and entraps sediment in areas with fresher water. Results also show that the lateral sediment transport induced by the semidiurnal tidal flow is significant when frictional forces are small (E ∼ 0.02). Model predictions are in good agreement with observations from the James River estuary.

[PDF]
[Abstract]

14 

Analysis of tidal straining as driver for estuarine circulation in wellmixed estuaries
Tidal straining, which can mathematically be described as the covariance between eddy viscosity and vertical shear of the alongchannel velocity component, has been acknowledged as one of the major drivers for estuarine circulation in channelized tidally energetic estuaries. In this paper, the authors investigate the role of lateral circulation for generating this covariance. Five numerical experiments are carried out, starting with a reference scenario including the full physics and four scenarios in which specific key physical processes are neglected. These processes are longitudinal internal pressure gradient forcing, lateral internal pressure gradient forcing, lateral advection, and the neglect of temporal variation of eddy viscosity. The results for the viscosity–shear covariance are correlated across different experiments to quantify the change due to neglect of these key processes. It is found that the lateral advection of vertical shear of the alongchannel velocity component and its interaction with the tidally asymmetric eddy viscosity (which is also modified by the lateral circulation) is the major driving force for estuarine circulation in wellmixed tidal estuaries.

[PDF]
[Abstract]

15 

Observations of barrier island length explained using an exploratory morphodynamic model
Barrier coasts display a chain of islands, separated by tidal inlets that connect a backbarrier basin to a sea or ocean. Observations show that barrier island length generally decreases for increasing tidal range and increasing basin area. However, this has neither been reproduced in model studies nor explained from the underlying physics. This is the aim of our study. Here we simulate barrier coast dynamics by combining a widely used empirical relationship for inlet dynamics with a processbased model of the tidal hydrodynamics. Our model results show stable inlet systems with more than one inlet open that support the observed qualitative relationships and fit in existing barrier coast classifications. To explain this, we identify a competition between a destabilizing mechanism (bottom friction in inlets, tending to reduce the number of open inlets) and a stabilizing one (spatially varying pressure gradients over the inlets, tending to keep the inlets open).

[PDF]
[Abstract]

16 

Channel and shoal development in a short tidal embayment: an idealized model study
In many tidal embayments, complex patterns of channels and shoals are observed. To gain a better understanding of these features, an idealized model, that describes the interaction of water motion, sediment transport and bed evolution in a semienclosed, rectangular basin, is developed and analysed. To explain the initial formation of channels and shoals, twodimensional perturbations superposed on a laterally uniform equilibrium bottom are studied. These perturbations evolve due to convergences of various residual suspended sediment fluxes: a diffusive flux, a flux related to the
bed topography, an advective flux resulting from internally generated overtides and an advective flux due to externally prescribed overtides. For most combinations of these fluxes, perturbations start to grow if the bottom friction is strong enough. Their growth is mainly a result of convergences of diffusive and topographically induced sediment fluxes. Advective contributions due to internally generated overtides enhance this growth. If only diffusive sediment fluxes are considered, the underlying equilibrium is always unstable. This can be traced back to the depth dependence of the deposition parameter. Contrary to the results of previous idealized models, the channels and shoals always initiate in the shallow, landward areas. This is explained by the enhanced generation (compared to that in previous models) of frictional torques in shallow regions. The resulting initial channel–shoal formation compares well with results found in complex numerical model studies. The instability mechanism and the location of the initial formation of bottom patterns do not change qualitatively when varying parameters. Changes are mainly related to differences in the underlying
equilibrium profile due to parameter variations.

[PDF]
[Abstract]

17 

Effect of bottom stress formulation on modelled flow and turbidity maxima in crosssections of tidedominated estuaries
A threedimensional numerical model with a prognostic salinity field is used to investigate the effect
of a partial slip bottom boundary condition on lateral flow and sediment distribution in a transect of a tidally dominated channel. The transect has a symmetrical Gaussian crosschannel bottom profile. For a deep, wellmixed, tidally dominated channel, partial slip decreases the relative importance of Coriolis deflection on the generation of crosschannel flow patterns. This has profound implications for the lateral distribution of residual salinity that drives the crosschannel residual circulation pattern. Transverse sediment transport, however, is always found to be governed by a balance between advection of residual sediment concentration by residual lateral flow on the one hand and crosschannel diffusion on the other hand. Hence, the changes in the crosschannel distribution of residual
salinity modify the lateral sediment distribution. For no slip, a single turbidity maximum occurs. In contrast, partial slip gives a gradual transition to a symmetrical density distribution with a turbidity maximum near each bank. For a more shallow, partially mixed tidal channel that represents the James River, a single turbidity maximum at the left bank is found irrespective of the nearbed slip condition. In this case, semidiurnal contributions to sediment distribution and lateral flow play an important role in crosschannel sediment transport. As vertical viscosity and diffusivity are increased, a second maximum at the right bank again exists for partial slip.

[PDF]
[Abstract]

18 

Influence of topography on tide propagation and amplification in semienclosed basins
An idealized model for tide propagation and amplification in semienclosed 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 depthaveraged 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 semidiurnal and diurnal
tide observations in the Gulf of California, the Adriatic Sea and the Persian Gulf.

[PDF]
[Abstract]

19 

Importance of crosschannel bathymetry and eddy viscosity parameterisation in modelling estuarine flow
For a proper understanding of flow patterns in curved tidal channels, quantification of contributions from
individual physical mechanisms is essential. We study quantitatively how such contributions are affected by crosschannel bathymetry and three alternative eddy viscosity parameterisations. Two models are presented for this purpose, both describing flow in curved but otherwise prismatic channels with an (almost) arbitrary transverse bathymetry. One is a numerical model based on the full threedimensional
shallow water equations. Special feature of this diagnostic model is that assumptions regarding the relative importance of particular physical mechanisms can be incorporated in the computations by switching corresponding terms in the model equations on or off. We also present an idealized model that provides semianalytical approximate solutions of the shallow water equations for all three considered
alternative eddy viscosity parameterisations. It forms an aid in explaining and theorising about results
obtained with the numerical model. Observations regarding Chesapeake Bay serve as a reference case for the present study. We find that the relative importance of both along channel advective forcing and transverse diffusive forcing depends on local characteristics of the crosssectional bottom profile rather than global ones. In our reference case, tideresidual alongchannel flow induced by these forcings
is not small compared to the total tidal residual. Building on this observation, we present an indicative test to judge whether advective processes should be included in leading order in modelling tidedominated estuarine flow. Furthermore, depending on the applied eddy viscosity parameterisation
(uniformly or parabolically distributed over the vertical), we find qualitatively different spatial patterns for
the alongchannel advective forcing.

[PDF]
[Abstract]

20 

Morphodynamic instabilities of planar beaches: Sensitivity to parameter values and process formulations
The initial growth of bed perturbations on planar sloping beaches under the forcing of obliquely incident, breaking waves is investigated using a state‐of‐the‐art numerical model. This allows for a systematic investigation of the sensitivity of the spatial structures of the bed perturbations and their growth and migration rates to different model formulations and parameterizations. If the sediment is only transported in the direction of the net current velocity and sediment stirring is taken proportional to the wave height squared, growing up‐current oriented crescentic bars are found with a preferred spacing of 800 m and a down‐current migration rate of 10 m d−1. Varying the angle of wave incidence, drag coefficient and bed slope results in qualitatively similar growing bed forms. Using an Engelund and Hansen transport formula, very oblique down‐current oriented bars are obtained that grow in time. No preferred wavelength, however, is found. Using the Bailard transport formula results in growing, up‐current oriented bars with a preferred spacing smaller than 300 m for wave angles smaller than 7°. When using either the Engelund and Hansen or Bailard sediment transport formulation, it is essential to take the transport in the direction of the wave orbital velocity into account in order to have growing bed perturbations.

[PDF]
[Abstract]
