Search

A field campaign was carried out at a sheltered sandy beach with the aim of gaining new insights into the driving processes behind sheltered beach morphodynamics. Detailed measurements of the local hydrodynamics, bed-level changes and sediment composition were collected at a man-made beach on the leeside of the barrier island Texel, bordering the Marsdiep basin that is part of the Dutch Wadden Sea. The dataset consists of (1) current, wave and turbidity measurements from a dense cross-shore array and a 3 km alongshore array; (2) sediment composition data from beach surface samples; (3) high-temporal-resolution RTK-GNSS beach profile measurements; (4) a pre-campaign spatially covering topobathy map; and (5) meteorological data. This paper outlines how these measurements were set up and how the data have been processed, stored and can be accessed. The novelty of this dataset lies in the detailed approach to resolve forcing conditions on a sheltered beach, where morphological evolution is governed by a subtle interplay between tidal and wind-driven currents, waves and bed composition, primarily due to the low-energy (near-threshold) forcing. The data are publicly available at 4TU Centre for Research Data at: https://doi.org/10.4121/19c5676c-9cea-49d0-b7a3-7c627e436541 (Van der Lugt et al., 2023)., Coastal Engineering, Lab Hydraulic Engineering, Environmental Fluid Mechanics, Civil Engineering & Geosciences

Coastal wave forecasting over large spatial scales is essential for many applications (e.g., coastal safety assessments, coastal management and developments, etc.). This demand explains the necessity for accurate yet effective models. A well-known efficient modelling approach is the quadratic approach (often referred to as frequency-domain models, nonlinear mild-slope models, amplitude models, etc.). The efficiency of this approach stems from a significant modelling reduction of the original governing equations (e.g., Euler equations). Most significantly, the description of wave nonlinearity essentially collapses into a single mode coupling term determined by the quadratic interaction coefficients. As a result, it is expected that the efficiency achieved by the quadratic approach is accompanied by a decrease in prediction accuracy. In order to gain further insight into the predictive capabilities of this modelling approach, this study examines six different quadratic formulations, three of which are of the Boussinesq type and the other three are referred to as fully dispersive. It is found that while the Boussinesq formulations reliably predict the evolution of coastal waves, the predictions by the fully dispersive formulations tend to be affected by false developments of modulational instability. Consequently, the predicted wave fields by the fully dispersive formulations are characterized by unexpectedly strong modulations of the sea-swell part and associated unexpected infragravity response. The impact of the modulational instability on wave prediction based on the quadratic approach is further demonstrated using existing laboratory results of bichromatic and irregular wave conditions., Environmental Fluid Mechanics

Spectral information of coastal waves and the associated statistical parameters (e.g., the significant wave height and mean wave period) over large spatial scales is essential for many applications (e.g., coastal safety assessments, coastal management and developments, etc.). This demand explains the necessity for accurate yet effective models. A well-known efficient modelling approach is the quadratic approach (often referred to as frequency-domain models, weakly nonlinear mild-slope models, amplitude models, etc.). The efficiency of this approach is achieved through modelling reduction of the original governing equations (e.g., Euler equations). Most significantly, wave nonlinearity is described solely by a single quadratic mode-coupling term. Therefore, doubts arise with regard to the predictive capabilities of the quadratic approach to reliably describe the nonlinear development of waves in the coastal environment where nonlinearity is typically significant. This study attempts to push the limit of the prediction capabilities of nonlinear coastal waves based on the quadratic approach. To this end, an optimization process is proposed, striving to extract the quadratic formulation which describes most adequately nonlinear wave developments over water depths and bathymetrical structures which characterize the coastal environment. The outcome is the model QuadWave1D: a fully dispersive quadratic model for coastal wave prediction in one-dimension. Based on a wide set of examples (including monochromatic, bichromatic and irregular wave conditions) and comparing to other representative quadratic formulations, it is found that QuadWave1D presents superior predictive capabilities of both the sea-swell components and the infragravity field., Environmental Fluid Mechanics

Nearshore hydro- and morphodynamic data were collected during a field experiment under calm conditions, moderate conditions, and storm conditions with dune erosion in the collision regime. The experiment was conducted on the Sand Engine near Kijkduin, the Netherlands, from October 18, 2021, to January 7, 2022. Two artificial unvegetated dunes were constructed just above the high water line to measure storm erosion and dune impacts from higher water levels and waves. During the experiment, three storms occurred that resulted in significant erosion of both dunes. The collected hydrodynamic data include pressure sensor and velocimeter data along two cross-shore transects. The collected morphodynamic data include bathymetry and topography surveys, optical backscatter sensor data in the inner surf zone, and a continuous cross-shore line-scanning lidar data set of the dune face. This comprehensive data set can be used to (1) study relevant nearshore hydrodynamic and morphodynamic processes that occur during calm conditions, moderate conditions, and storm conditions with dune erosion in the collision regime, and (2) validate existing dune erosion models., Coastal Engineering, Environmental Fluid Mechanics

Infragravity (IG) waves are key drivers for coastal erosion and thus need to be properly included in process-based modelling of coastal hazards. Uncertainties remain regarding the offshore boundary conditions for these long waves. Typically, only bound IG waves are included at the boundary, which means that the possible contribution of free IG waves, such as those radiated from distant coastlines, is neglected. Recent studies however suggest that incoming free IG waves could be significant, particularly in semi-enclosed basins such as the North Sea where they could contribute to coastal hazards (e.g., Reniers et al., 2021, Rijnsdorp et al. 2021). The objective of this work is to improve the understanding of the incoming IG wave field along the Dutch coast. We will quantify how bound and free IG waves develop in intermediate water depths and assess in which conditions (onshore directed) free IG waves become significant., Environmental Fluid Mechanics, Coastal Engineering

The present work presents physical laboratory measurements of surface elevation and pore water pressures in a fine sand bed under bichromatic waves in a large-scale laboratory experiment. This was done at three cross-shore locations in the swash zone, with pressures being measured at different depths in the bed. The measurements show that the pore pressure signal decays and shifts with increased depth. These measurements are used to validate a practical model, based on the theory of Yamamoto et al. (1978, https://doi.org/10.1017/S0022112078003006) and Guest and Hay (2017, https://doi.org/10.1002/2016JC012257). The model corresponds well with the measurements (nRMSE < 0.2 and R2 > 0.95 for most probes) and shows that a frequency-based model can reproduce the pressures in the bed, despite the bed being exposed during dry periods. Furthermore, the model provides the opportunity to calculate pressure gradients, both throughout the bed and at the bed surface. These modeled pressure gradients at the bed surface show that the vertical pressure gradient can have an important impact on the Shields parameter, thereby influencing sediment transport., Environmental Fluid Mechanics

Currents can affect the evolution of waves in nearshore regions through altering their wavenumber and amplitude. Including the effect of ambient currents (e.g., tidal and wind-driven) on waves in phase-resolving wave models is not straightforward as it requires appropriate boundary conditions in combination with a large domain size and long simulation duration. In this paper, we extended the non-hydrostatic wave-flow model SWASH with additional terms that account for the influence of a depth-uniform ambient current on the wave dynamics, in which the current field can be taken from an external source (e.g., from observations or a circulation model). We verified the model ability by comparing predictions to results from linear theory, laboratory experiments and a spectral wave model that accounts for wave interference effects. With this extension, the model was able to account for current-induced changes to the wave field (i.e., changes to the wave amplitude, length and direction) due to following and opposing currents, and two classical examples of sheared currents (a jet-like current and vortex ring). Furthermore, the model captured the wave dynamics in the presence of strong opposing currents. This includes reflections of relatively small amplitude waves at the theoretical blocking point, and transmission of breaking waves beyond the theoretical blocking point for larger wave amplitudes. The proposed model extension allows phase-resolving models to more accurately and efficiently simulate the wave dynamics in coastal regions with tidal and/or wind-driven flows., Environmental Fluid Mechanics

The coastline of Demak, Indonesia, has been eroding during the last 15 years. Coastal retreat in Demak is caused by a combination of mangrove deforestation and local subsidence due to groundwater extraction in the nearby city of Semarang. To restore the lost mangrove forest, permeable dams, consisting of bamboo poles with a brushwood filling, have been built to attenuate the waves, facilitate sedimentation at their land side, and thus create a suitable habitat for mangroves. However, existing designs required frequent brushwood maintenance. Therefore, a new type of design is proposed, consisting of only vertical bamboo poles without a filling of brushwood. Nevertheless, the hydrodynamic performance of this type of structure is not known. This study assesses the wave transformation through structures formed by bamboo poles for the physical conditions of Demak, Indonesia, with the numerical wave model SWASH. Field measurements and WaveWatch III data are analyzed to obtain the design conditions for the structures in Demak. SWASH is validated against laboratory experiments, and applied to investigate different structure designs. The model shows that for a structure consisting of two rows of bamboo poles, the transmission rate Et/Ei decreases from 75% to 55% when the row spacing in the wave direction is increased from sx = 0.42 m to sx =5.8 m. Even larger spacings do not result in less transmission, and at least three rows are needed to have a transmission rate lower than 50 % - a common wave reduction target used in restoration efforts with structures. This study thus identifies potential strategies to maximize wave attenuation by bamboo structures, which can be used to reduce wave attack along muddy coasts without the need of a brushwood filling. Hereby it provides an economically and user friendly alternative with respect to the previous brushwood struct, Environmental Fluid Mechanics, Hydraulic Structures and Flood Risk

Dune erosion during storm surges can lead to excessive damage to the dune system with devastating floods as a potential consequence. A risk assessment of areas protected by dunes can be facilitated by an understanding and description of the physical processes that take place. Field measurements, knowledge of underlying processes and numerical modelling have developed with time, which enabled a more comprehensive description and new predictive techniques. This review concerns dune erosion in the collision regime, and summarises relevant observations, describes underlying processes and explains existing models predicting dune erosion. Observations of dune erosion consist of field observations, laboratory experiments and manipulative field campaigns. The underlying physical processes that contribute to dune erosion are divided into processes that contribute to sediment transport due to hydrodynamic forcing, which occurs in the surf and swash zone, and sediment transport due to avalanching, which occurs in the swash zone, on the dune face and on the dune crest. The existing dune erosion models that are discussed here contain (empirical) equilibrium profile models and process-based models, which can both be a valuable tool for the risk assessment of storm surges. However, model uncertainties still remain, as specific processes are not yet fully understood and described. Examples are the influences of wave obliquity, sediment grain size, and vegetation on the dune face. By improving our knowledge through research and reducing these uncertainties, we can further improve our predictive models. This could eventually lead to more accurate predictions, more complete risk assessments, and sandy coastlines which are more resilient to excessive dune erosion and possible floods., Coastal Engineering, Environmental Fluid Mechanics, Hydraulic Engineering

Many sand spits are morphodynamically complex landforms, that are either analysed with complex and expensive computational models or at a conceptual level. Therefore, most case studies on spits in different environments are descriptive. A novel method based on the use of polar coordinates was devised to quantitatively analyse spit morphodynamics in a non-tidal, wind-dominated lake environment, using the Marker Wadden islands in Lake Markermeer, the Netherlands, as a case study. A high-resolution morphological data set allowed for the quantification of sedimentation processes around two spits, in two distinctive depth zones. Spit-platform growth is governed by alongshore currents that transport sediment over the spit-platform into deeper waters; the size of the spit-platform in turn affects the growth of the spit around the mean water level. Insight in this complex interplay of processes is crucial to understand spit behaviour in low-energy lake environments. At the Marker Wadden the submerged spit-platform grows during high energy wind events while the emerged spit part grows under mild to moderate energy conditions. With this new method we can quantitatively explore the role of different wave and flow conditions and predict spit growth direction in non-tidal, wind-dominated environments, beyond the level of conceptual descriptions., Coastal Engineering, Environmental Fluid Mechanics, Hydraulic Engineering

Weyl rule of association, proposed by Hermann Weyl for quantum mechanics applications (Weyl, 1931), can be used to associate between the dispersion relation of water waves and a non-local pseudo-differential operator. The central result of this study is that this operator correctly approximates the Dirichlet-to-Neumann operator derived for linear waves over a slowly varying bathymetry. This opens the door to a formal use of Weyl's operational calculus, and consequently, allowing straightforward derivations and generalizations of water waves’ models over mild slopes. Specifically, within the framework of linear wave theory, the formulation based on Weyl rule of association provides a generalized mild-slope model which does not impose a limit on the spectral width. Most significantly, the mild-slope formulation based on Weyl rule of association allows to derive a general linear kinetic equation for which the widely used energy balance equation (the central equation of forecasting models such as SWAN and WAVEWATCH) serves as a special case. This result not only provides a formal link between the deterministic description (i.e., Euler equations) and the stochastic description (i.e., the energy balance equation), but also establishes the theoretical foundations for the statistical description of bathymetry-induced wave interferences. Such a statistical description is especially important over coastal waters, where through the interaction with the bathymetry, waves are rapidly scattered and tend to form focal zones and associated interference patterns., Environmental Fluid Mechanics

Grain size affects the rates of aeolian sediment transport on beaches. Sediment in coastal environments typically consists of multiple grain-size fractions and exhibits spatiotemporal variations. Still, conceptual and numerical aeolian transport models are simplified and often only include a single fraction that is constant over the model domain. It is unclear to what extent this simplification is valid and if the inclusion of multi-fraction transport and spatial grain-size variations affects aeolian sediment transport simulations and predictions of coastal dune development. This study applies the numerical aeolian sediment transport model AeoLiS to compare single-fraction to multi-fraction approaches for a range of grain-size distributions and spatial grain-size scenarios. The results show that on timescales of days to years, single-fraction simulations with the median grain size, D50, often give similar results to multi-fraction simulations, provided the wind is able to mobilize all fractions within that time frame. On these timescales, vertical variability in grain size has a limited effect on total transport rates, but it does influence the simulation results on minute timescales. Horizontal grain-size variability influences both the total transport rates and the downwind bed grain-size composition. The results provide new insights into the influence of beach sediment composition and spatial variability on total transport rates toward the dunes. The findings of this study can guide the implementation of grain-size variability in numerical aeolian sediment transport models., Coastal Engineering, Environmental Fluid Mechanics

The Simulating WAves Nearshore (SWAN) model has been extended with an infragravity module to predict the Wave-Group-Forced (WGF) infragravity response to a frequency-directional sea-swell spectrum at a mildly sloping alongshore uniform beach. To that end the SWAN model has been extended with an WGF-infragravity source term denoted Ssb where the subscript denotes surfbeat. The corresponding WGF infragravity energy model has been verified with a set of benchmark tests using the infragravity amplitude model of Reniers et al. (2002). Next the implementation of the energy balance in SWAN has been validated with both prototype-scale laboratory experiments and field observations, showing a good comparison with observations not affected by the nodal structure of the (partially) standing infragravity waves. This suggests that the model is capable of providing improved infragravity boundary conditions in relatively shallow water compared to the typical assumption of equilibrium forcing conditions using for instance Hasselmann's equilibrium theory (Hasselmann, 1962). These infragravity boundary conditions can subsequently can be used by other more sophisticated models to compute runup, overtopping and dune erosion., Environmental Fluid Mechanics

A chenier is a beach ridge, consisting of sand and/or shells, overlying a muddy substrate. In this paper, we explore the cross-shore dynamics of cheniers in their ‘active’ phase, i.e. the phase between their formation and their landing on the shore and can no longer be reached by daily wave and tidal influences. While cheniers described in literature are known to only migrate onshore until they reach a stable position with their crest level above tidal influences, observations in Demak suggest the existence of an alternative stable state, highly dynamic on the short term, but stable on the longer term. To explore this alternative stable state, we developed an idealised chenier model to investigate cross-shore chenier dynamics under daily wave and tidal influences. The model is able to predict both onshore and offshore migration; onshore migration is mainly driven by wave action, while offshore migration is induced by a tidal phase lag, or the effect of the storm season. For certain combinations of waves, tide (incl. phase lag) and a storm season effect, the model predicts a dynamically stable chenier. In absence of a phase lag and storm season effect, the model yields a ‘classic’ stable chenier that welds onto the shoreline by onshore migration., Environmental Fluid Mechanics

In sandy beach systems, the aeolian sediment transport can be governed by the vertical structure of the sediment layers at the bed surface. Here, data collected with a newly developed sand scraper is presented to determine high-resolution vertical grain size variability and how it is affected by marine and aeolian processes. Sediment samples at up to 2 mm vertical resolution down to 50 mm depth were collected at three beaches: Waldport (Oregon, USA), Noordwijk (the Netherlands) and Duck (North Carolina, USA). The results revealed that the grain size in individual layers can differ considerably from the median grain size of the total sample. The most distinct temporal variability occurred due to marine processes that resulted in significant morphological changes in the intertidal zone. The marine processes during high water resulted both in fining and coarsening of the surface sediment. Especially near the upper limit of wave runup, the formation of a veneer of coarse sediment was observed. Although the expected coarsening of the near-surface grain size during aeolian transport events was observed at times, the opposite trend also occurred. The latter could be explained by the formation and propagation of aeolian bedforms within the intertidal zone locally resulting in sediment fining at the bed surface. The presented data lays the basis for future sediment sampling strategies and sediment transport models that investigate the feedbacks between marine and aeolian transport, and the vertical variability of the grain size distribution., Coastal Engineering, Environmental Fluid Mechanics

Understanding directional spectra of infragravity (IG) waves composed of free and bound components is required due to their impacts on various coastal processes (e.g., coastal inundation and morphological change). However, conventional reconstruction methods of directional spectra relying on linear wave theory are not applicable to IG waves in intermediate water depths (20–30 m) due to the presence of bound waves. Herein, a novel method is proposed to reconstruct directional spectra of IG waves in intermediate depth based on weakly nonlinear wave theory. This method corrects cross-spectra among observed wave signals by taking account of the nonlinearity of bound waves in order to reconstruct directional spectra of free IG waves. Numerical experiments using synthetic data representing various directional distributions show that the proposed method reconstructs free IG wave directional spectra more accurately than the conventional method. The method is subsequently applied to observations of severe sea-states at two field sites. At these sites, free IG waves are not isotropic and have clear peak directions. Numerical modeling of the wave fields shows that these peak directions correspond to the reflection of IG waves from the shore and/or coastal structures. Additionally, the validity of the underlying weakly nonlinear wave theory of the present method is assessed by a newly proposed method employing bispectral analysis. The bound wave response generally agrees with the theory at the field sites but deviates slightly for energetic sea states. The applicability of the present method on a sloping bottom is further discussed by an analytical solution., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Environmental Fluid Mechanics

Cheniers are ridges consisting of coarse-grained sediments, resting on top of the fine sediment that forms the otherwise muddy coast. In this paper, we use Delft3D to explore how cheniers are formed through wave winnowing. We identify three phases of chenier development: (a) a winnowing phase, during which mud is washed out of the seabed initially consisting of a mixture of sand and mud, (b) a sand transport phase, when the sand in the upper layer is transported onshore, and (c) a crest formation phase, during which a chenier crest rapidly develops at the landward limit of onshore sediment transport. The main mechanism driving onshore sand transport is wave asymmetry. During calm conditions, sand transport takes place within a narrow band limiting the volume of sand delivered nearshore, and therefore no chenier develops. In contrast, average storm conditions mobilize sufficient sand for a crest to develop. Our results thus reveal that chenier formation through wave winnowing does not require extreme storm conditions. Furthermore, our study showed that chenier formation through wave winnowing is a relatively slow process, with the largest time scales associated with the winnowing and sand transport. Once sufficient sand is available in the intertidal zone, the crest develops rapidly., Environmental Fluid Mechanics

We present a fully coupled 2DV morphodynamic model, implemented in OpenFOAM® that is capable of simulating swash-zone morphodynamics of sandy beaches. The hydrodynamics are described by the Reynolds-averaged Navier–Stokes (RANS) equations with a (Formula presented.) turbulence model and the Volume of Fluid (VoF) approach for discriminating between air and water. Sediment transport is described in terms of bedload and suspended load transport. We show that the default divergence scheme in OpenFOAM can become numerically unstable and lead to negative sediment concentrations, and propose a solution to avoid this problem. The model performance is assessed in terms of surface elevation, flow velocities, runup, suspended sediment concentrations, bed profile evolution and sediment transport volumes by comparing with measurements of field-scale (wave height of 0.6 m) solitary waves. The model shows reasonable agreement in terms of hydrodynamics and predicts the correct sediment transport volumes, although the deposition is predicted more onshore compared to the measurements. This is partially attributed to an overprediction of the runup. The model shows that the suspended sediment concentration displays a strong vertical dependence. These results show the potential of depth-resolving models in providing more insight into morphodynamic processes in the swash zone, particularly with respect to vertical structures in the flow and suspended sediment transport., Environmental Fluid Mechanics

To address the important research question of whether implicit (bottom friction) or explicit (stem drag) dissipation models are most appropriate for the prediction of wave attenuation due to aquatic vegetation, the Simulating Waves Nearshore (SWAN) spectral wave model has been extended with an explicit frequency-dependent dissipation model for submerged and emergent vegetation. The new explicit model is compared to existing explicit and implicit dissipation models in SWAN, and the distinguishing features of each of the dissipation models are quantified. The present work verifies the implementation of the new and existing dissipation models, outlines their distinguishing features, and compares model predictions against experimental data. The emphasis is on the transformation of the spectral wave periods Tm0;1 and Tm 1;0 over a canopy. Model evaluation based on academic and laboratory cases allows for recommendations regarding applicability of the three dissipation models, where the new method has the broadest applicability, since it bridges the gap in applicability between the other two dissipation models. The implementation of Jacobsen, McFall, and van der A (2019; A frequency distributed dissipation model for canopies; Coastal Engineering, 150, 135-146) is publicly available in SWAN version 41.31B., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Coastal Engineering, Environmental Fluid Mechanics

A chenier is a beach ridge, consisting of sand and/or shells, overlying a muddy substrate. In this paper, we explore the cross-shore dynamics of cheniers in their ‘active’ phase, i.e. the phase between their formation and their landing on the shore and can no longer be reached by daily wave and tidal influences. While cheniers described in literature are known to only migrate onshore until they reach a stable position with their crest level above tidal influences, observations in Demak suggest the existence of an alternative stable state, highly dynamic on the short term, but stable on the longer term. To explore this alternative stable state, we developed an idealised chenier model to investigate cross-shore chenier dynamics under daily wave and tidal influences. The model is able to predict both onshore and offshore migration; onshore migration is mainly driven by wave action, while offshore migration is induced by a tidal phase lag, or the effect of the storm season. For certain combinations of waves, tide (incl. phase lag) and a storm season effect, the model predicts a dynamically stable chenier. In absence of a phase lag and storm season effect, the model yields a ‘classic’ stable chenier that welds onto the shoreline by onshore migration., Environmental Fluid Mechanics

The Simulating WAves Nearshore (SWAN) model has been extended with an infragravity module to predict the Wave-Group-Forced (WGF) infragravity response to a frequency-directional sea-swell spectrum at a mildly sloping alongshore uniform beach. To that end the SWAN model has been extended with an WGF-infragravity source term denoted Ssb where the subscript denotes surfbeat. The corresponding WGF infragravity energy model has been verified with a set of benchmark tests using the infragravity amplitude model of Reniers et al. (2002). Next the implementation of the energy balance in SWAN has been validated with both prototype-scale laboratory experiments and field observations, showing a good comparison with observations not affected by the nodal structure of the (partially) standing infragravity waves. This suggests that the model is capable of providing improved infragravity boundary conditions in relatively shallow water compared to the typical assumption of equilibrium forcing conditions using for instance Hasselmann's equilibrium theory (Hasselmann, 1962). These infragravity boundary conditions can subsequently can be used by other more sophisticated models to compute runup, overtopping and dune erosion., Environmental Fluid Mechanics

In this paper, a new approach to model wave‐driven, cross‐shore shoreline change incorporating multiple timescales is introduced. As a base, we use the equilibrium shoreline prediction model ShoreFor that accounts for a single timescale only. High‐resolution shoreline data collected at three distinctly different study sites is used to train the new data‐driven model. In addition to the direct forcing approach used in most models, here two additional terms are introduced: a time-upscaling and a time‐downscaling term. The upscaling term accounts for the persistent effect of short‐term events, such as storms, on the shoreline position. The downscaling term accounts for the effect of long‐term shoreline modulations, caused by, for example, climate variability, on shorter event impacts. The multi‐timescale model shows improvement compared to the original ShoreFor model (a normalized mean square error improvement during validation of 18 to 59%) at the three contrasted sandy beaches. Moreover, it gains insight in the various timescales (storms to inter‐an-nual) and reveals their interactions that cause shoreline change. We find that extreme forcing events have a persistent shoreline impact and cause 57–73% of the shoreline variability at the three sites. Moreover, long‐term shoreline trends affect short‐term forcing event impacts and determine 20– 27% of the shoreline variability., Coastal Engineering, Environmental Fluid Mechanics

Seaport operability is key to the economic viability of ports. Metocean conditions (e.g., wind, short waves, and infragravity waves) affect this operability when certain thresholds are exceeded. This paper describes a method for the global mapping of seaport operability risk indicators using open-source metocean data. This global-scale assessment provides a geographic overview of operability risks and first-order insights into the most relevant metocean risk indicators at each location. The results show that locations around the equator and inland seas have lower operability risk than locations farther away from the equator. “Hotspots” are mainly located along the southern capes (Cape of Good Hope, Leeuwin, Horn), around the ‘Roaring Forties’, and at exposed locations along the oceans. Of the metocean parameters considered, short waves are found to be the most critical risk indicator for port operability at most locations. Using (the insights of) this study, port authorities, operators, and designers can prepare for metocean risks at an early stage and effectively respond with mitigation measures and layout adjustments to improve port operability., Materials and Environment, Rivers, Ports, Waterways and Dredging Engineering, Environmental Fluid Mechanics

Infragravity waves are low-frequency surface waves that can impact a variety of nearshore and oceanic processes. Recent measurements in the North Sea showed that significant bursts of infragravity energy occurred during storm events. Using a spectral wave model, we show that a substantial part of this energy was radiated from distant shorelines where it was generated by the incident sea-swell waves. These radiated infragravity waves can cross the North Sea basin and reach distant shorelines. The origin of the infragravity wave energy varied between the different storms, and particularly depends on where largest sea-swell waves made landfall. Along the coastlines of the North Sea, shoreward directed infragravity waves that originate from a remote source were non-negligible during storm events. This suggests that radiated infragravity waves can potentially contribute to coastal dynamics and hazards away from their region of generation., Environmental Fluid Mechanics

In numerical ocean models, the effect of waves on currents is usually expressed by either vortex-force or radiation stress representations. In this paper, the differences and similarities between those two representations are investigated in detail in conditions of both conservative and nonconservative waves. In addition, comparisons between different sets of equations of mean motion that apply different representations of wave-induced forcing terms are included. The comparisons are useful for selecting a suitable numerical ocean model to simulate the mean current in conditions of waves combined with currents., Coastal Engineering, Environmental Fluid Mechanics

A 2-layer non-hydrostatic model with improved dispersive behaviour is presented. Due to the assumption of a constant non-hydrostatic pressure distribution in the lower layer, the dispersive behaviour is improved without much additional computational time. A comparison with linear wave theory showed that this 2-layer model gives a better result for the dispersion relation and shoaling of waves in intermediate water. This means that the 2-layer model is applicable in shallow and intermediate water depths (up to relative depths kh equals 4), whereas the 1-layer model is only applicable in shallow water depths (kh smaller than 1). Three laboratory experiments, including a fringing reef and a barred beach, were used to validate the presented mode for different hydrodynamic conditions. Based on these results, it can be concluded that the 2-layer model can be applied to accurately simulate the bulk wave height and spectral properties. The low frequency wave height, the setup and in particular the second order statistics contain more scatter, but the model accurately captured the general trend. Furthermore, the model showed good results for complex bathymetries in shallow to intermediate water., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Environmental Fluid Mechanics

Mangrove vegetation constitutes a natural coastal defence against waves and erosion. Despite their protective role, mangrove ecosystems have experienced continuous degradation over the last decades due to human causes. At retreating mangrove coastlines, bamboo structures are built to create new habitat for mangrove colonization. Existing structures have experienced mixed rates of success due to the lack of a scientific basis in their design. Optimizing future structure designs requires investigating the effect of the bamboo poles on waves. We consequently conducted laboratory experiments to measure wave transformation, hydrodynamic forces, and flow velocities inside cylinder arrays, mimicking bamboo poles, with varying cylinder configurations and orientations. The experiments provided relationships for wave transmission, wave reflection, and the drag coefficients for configurations with volumetric porosities between n = 0.64 − 0.9. Configurations with a small lateral spacing (causing higher blockage) and a relatively longer streamwise spacing (causing less sheltering) exhibit larger forces and dissipation per element. Such arrangements enable optimizing wave dissipation at locations where the wave direction has low variability over the year. Placing the poles horizontally instead of vertically increases the forces and wave dissipation per element in relatively deeper water. Based on the experiments, we developed a conceptual analytical model that predicts wave reflection and dissipation through cylinder arrays, including blockage and sheltering. The model can reproduce the influence of cylinder arrangement on wave transformation, and it suggests that accurate predictions of sheltering and wave reflection are important to find optimal designs. Overall, these results provide useful insights on how to model and optimize the design of structures for mangrove restoration., Hydraulic Structures and Flood Risk, Environmental Fluid Mechanics

The prediction of wave runup, as well as its components, time-averaged setup and the time-varying swash, is a key element of coastal storm hazard assessments, as wave runup controls the transitions between morphodynamic response types such as dune erosion and overwash, and the potential for flooding by wave overtopping. While theoretically able to simulate the dominant low-frequency swash, previous studies using the infragravity-wave–resolving model XBeach (XBSB) have shown an underestimation of the observed swash variance and wave runup, which was in part related to the absence of incident-band swash motions in the model. Here, we use an incident-band wave-resolving, non-hydrostatic version of the XBeach model (XBNH) to simulate wave runup observed during the SandyDuck '97 experiment on an intermediate–reflective sandy beach. The results show that the XBNH model describes wave runup and the individual setup and swash components well. We subsequently examine differences in wave runup prediction between the XBSB and XBNH models and find that the XBNH model is a better predictor of wave runup than XBSB for this beach, which is due to better predictions of both the incident-band and infragravity-band swash. For a range of beach states from reflective to dissipative it is shown that incident-band swash is underestimated by XBSB relative to XBNH, in particular for reflective conditions. Infragravity-band swash is shown to be lower in XBSB than XBNH for most conditions, including dissipative conditions for which the mean difference is 16% of the deep water wave height. The difference in infragravity-band swash in XBNH relative to XBSB is shown to mainly be the result of processes occurring outside the swash zone, but approximately 15% of the difference is caused by explicitly resolving incident-band wave motions within the swash zone, such as swash-swash interactions, which inherently cannot be simulated by wave-averaged models., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Environmental Fluid Mechanics

Coastal safety assessments with wave-resolving storm impact models require a proper offshore description for the incoming infragravity (IG) waves. This boundary condition is generally obtained by assuming a local equilibrium between the directionally-spread incident sea-swell wave forcing and the bound IG waves. The contribution of the free incident IG waves is thus ignored. Here, in-situ observations of IG waves with wave periods between 100 s and 200 s at three measurement stations in the North Sea in water depths of O(30) m are analyzed to explore the potential contribution of the free and bound IG waves to the total IG wave height for the period from 2010 to 2018. The bound IG wave height is computed with the equilibrium theory of Hasselmann using the measured frequency-directional sea-swell spectra as input. The largest IG waves are observed in the open sea with a maximum significant IG wave height of O(0.3) m at 32 m water depth during storm Xaver (December 2013) with a concurrent significant sea-swell wave height in excess of 9 m. Along the northern part of the Dutch coast, this maximum has reduced to O(0.2) m at a water depth of 28 m with a significant sea-swell wave height of 7 m and to O(0.1) m at the most southern location at a water depth of 34 m with a significant sea-swell wave height of 5 m. These appreciable IG wave heights in O(30) m water depth represent a lower bound for the expected maximum IG wave heights given the fact that in the present analysis only a fraction of the full IG frequency range is considered. Comparisons with the predicted bound IG waves show that these can contribute substantially to the observed total IG wave height during storm conditions. The ratio between the predicted bound-and observed total IG variance ranges from 10% to 100% depending on the location of the observations and the timing during the storm. The ratio is typically high at the peak of the storm and is lower at both the onset and waning of the storm. There i, Environmental Fluid Mechanics, Coastal Engineering

High quality laboratory measurements of nearshore waves and morphology change at, or near prototype-scale are essential to support new understanding of coastal processes and enable the development and validation of predictive models. The DynaRev experiment was completed at the GWK large wave flume over 8 weeks during 2017 to investigate the response of a sandy beach to water level rise and varying wave conditions with and without a dynamic cobble berm revetment, as well as the resilience of the revetment itself. A large array of instrumentation was used throughout the experiment to capture: (1) wave transformation from intermediate water depths to the runup limit at high spatio-temporal resolution, (2) beach profile change including wave-by-wave changes in the swash zone, (3) detailed hydro and morphodynamic measurements around a developing and a translating sandbar., Corrigendum DOI: 10.1038/s41597-021-00874-2, Environmental Fluid Mechanics

In numerical ocean models, the effect of waves on currents is usually expressed by either vortex-force or radiation stress representations. In this paper, the differences and similarities between those two representations are investigated in detail in conditions of both conservative and nonconservative waves. In addition, comparisons between different sets of equations of mean motion that apply different representations of wave-induced forcing terms are included. The comparisons are useful for selecting a suitable numerical ocean model to simulate the mean current in conditions of waves combined with currents., Coastal Engineering, Environmental Fluid Mechanics

Infragravity waves are low-frequency surface waves that can impact a variety of nearshore and oceanic processes. Recent measurements in the North Sea showed that significant bursts of infragravity energy occurred during storm events. Using a spectral wave model, we show that a substantial part of this energy was radiated from distant shorelines where it was generated by the incident sea-swell waves. These radiated infragravity waves can cross the North Sea basin and reach distant shorelines. The origin of the infragravity wave energy varied between the different storms, and particularly depends on where largest sea-swell waves made landfall. Along the coastlines of the North Sea, shoreward directed infragravity waves that originate from a remote source were non-negligible during storm events. This suggests that radiated infragravity waves can potentially contribute to coastal dynamics and hazards away from their region of generation., Environmental Fluid Mechanics

Seaport operability is key to the economic viability of ports. Metocean conditions (e.g., wind, short waves, and infragravity waves) affect this operability when certain thresholds are exceeded. This paper describes a method for the global mapping of seaport operability risk indicators using open-source metocean data. This global-scale assessment provides a geographic overview of operability risks and first-order insights into the most relevant metocean risk indicators at each location. The results show that locations around the equator and inland seas have lower operability risk than locations farther away from the equator. “Hotspots” are mainly located along the southern capes (Cape of Good Hope, Leeuwin, Horn), around the ‘Roaring Forties’, and at exposed locations along the oceans. Of the metocean parameters considered, short waves are found to be the most critical risk indicator for port operability at most locations. Using (the insights of) this study, port authorities, operators, and designers can prepare for metocean risks at an early stage and effectively respond with mitigation measures and layout adjustments to improve port operability., Materials and Environment, Rivers, Ports, Waterways and Dredging Engineering, Environmental Fluid Mechanics

In this paper, a new approach to model wave‐driven, cross‐shore shoreline change incorporating multiple timescales is introduced. As a base, we use the equilibrium shoreline prediction model ShoreFor that accounts for a single timescale only. High‐resolution shoreline data collected at three distinctly different study sites is used to train the new data‐driven model. In addition to the direct forcing approach used in most models, here two additional terms are introduced: a time-upscaling and a time‐downscaling term. The upscaling term accounts for the persistent effect of short‐term events, such as storms, on the shoreline position. The downscaling term accounts for the effect of long‐term shoreline modulations, caused by, for example, climate variability, on shorter event impacts. The multi‐timescale model shows improvement compared to the original ShoreFor model (a normalized mean square error improvement during validation of 18 to 59%) at the three contrasted sandy beaches. Moreover, it gains insight in the various timescales (storms to inter‐an-nual) and reveals their interactions that cause shoreline change. We find that extreme forcing events have a persistent shoreline impact and cause 57–73% of the shoreline variability at the three sites. Moreover, long‐term shoreline trends affect short‐term forcing event impacts and determine 20– 27% of the shoreline variability., Coastal Engineering, Environmental Fluid Mechanics

A 2-layer non-hydrostatic model with improved dispersive behaviour is presented. Due to the assumption of a constant non-hydrostatic pressure distribution in the lower layer, the dispersive behaviour is improved without much additional computational time. A comparison with linear wave theory showed that this 2-layer model gives a better result for the dispersion relation and shoaling of waves in intermediate water. This means that the 2-layer model is applicable in shallow and intermediate water depths (up to relative depths kh equals 4), whereas the 1-layer model is only applicable in shallow water depths (kh smaller than 1). Three laboratory experiments, including a fringing reef and a barred beach, were used to validate the presented mode for different hydrodynamic conditions. Based on these results, it can be concluded that the 2-layer model can be applied to accurately simulate the bulk wave height and spectral properties. The low frequency wave height, the setup and in particular the second order statistics contain more scatter, but the model accurately captured the general trend. Furthermore, the model showed good results for complex bathymetries in shallow to intermediate water., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Environmental Fluid Mechanics

Mangrove vegetation constitutes a natural coastal defence against waves and erosion. Despite their protective role, mangrove ecosystems have experienced continuous degradation over the last decades due to human causes. At retreating mangrove coastlines, bamboo structures are built to create new habitat for mangrove colonization. Existing structures have experienced mixed rates of success due to the lack of a scientific basis in their design. Optimizing future structure designs requires investigating the effect of the bamboo poles on waves. We consequently conducted laboratory experiments to measure wave transformation, hydrodynamic forces, and flow velocities inside cylinder arrays, mimicking bamboo poles, with varying cylinder configurations and orientations. The experiments provided relationships for wave transmission, wave reflection, and the drag coefficients for configurations with volumetric porosities between n = 0.64 − 0.9. Configurations with a small lateral spacing (causing higher blockage) and a relatively longer streamwise spacing (causing less sheltering) exhibit larger forces and dissipation per element. Such arrangements enable optimizing wave dissipation at locations where the wave direction has low variability over the year. Placing the poles horizontally instead of vertically increases the forces and wave dissipation per element in relatively deeper water. Based on the experiments, we developed a conceptual analytical model that predicts wave reflection and dissipation through cylinder arrays, including blockage and sheltering. The model can reproduce the influence of cylinder arrangement on wave transformation, and it suggests that accurate predictions of sheltering and wave reflection are important to find optimal designs. Overall, these results provide useful insights on how to model and optimize the design of structures for mangrove restoration., Hydraulic Structures and Flood Risk, Environmental Fluid Mechanics

The prediction of wave runup, as well as its components, time-averaged setup and the time-varying swash, is a key element of coastal storm hazard assessments, as wave runup controls the transitions between morphodynamic response types such as dune erosion and overwash, and the potential for flooding by wave overtopping. While theoretically able to simulate the dominant low-frequency swash, previous studies using the infragravity-wave–resolving model XBeach (XBSB) have shown an underestimation of the observed swash variance and wave runup, which was in part related to the absence of incident-band swash motions in the model. Here, we use an incident-band wave-resolving, non-hydrostatic version of the XBeach model (XBNH) to simulate wave runup observed during the SandyDuck '97 experiment on an intermediate–reflective sandy beach. The results show that the XBNH model describes wave runup and the individual setup and swash components well. We subsequently examine differences in wave runup prediction between the XBSB and XBNH models and find that the XBNH model is a better predictor of wave runup than XBSB for this beach, which is due to better predictions of both the incident-band and infragravity-band swash. For a range of beach states from reflective to dissipative it is shown that incident-band swash is underestimated by XBSB relative to XBNH, in particular for reflective conditions. Infragravity-band swash is shown to be lower in XBSB than XBNH for most conditions, including dissipative conditions for which the mean difference is 16% of the deep water wave height. The difference in infragravity-band swash in XBNH relative to XBSB is shown to mainly be the result of processes occurring outside the swash zone, but approximately 15% of the difference is caused by explicitly resolving incident-band wave motions within the swash zone, such as swash-swash interactions, which inherently cannot be simulated by wave-averaged models., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Environmental Fluid Mechanics

Coastal safety assessments with wave-resolving storm impact models require a proper offshore description for the incoming infragravity (IG) waves. This boundary condition is generally obtained by assuming a local equilibrium between the directionally-spread incident sea-swell wave forcing and the bound IG waves. The contribution of the free incident IG waves is thus ignored. Here, in-situ observations of IG waves with wave periods between 100 s and 200 s at three measurement stations in the North Sea in water depths of O(30) m are analyzed to explore the potential contribution of the free and bound IG waves to the total IG wave height for the period from 2010 to 2018. The bound IG wave height is computed with the equilibrium theory of Hasselmann using the measured frequency-directional sea-swell spectra as input. The largest IG waves are observed in the open sea with a maximum significant IG wave height of O(0.3) m at 32 m water depth during storm Xaver (December 2013) with a concurrent significant sea-swell wave height in excess of 9 m. Along the northern part of the Dutch coast, this maximum has reduced to O(0.2) m at a water depth of 28 m with a significant sea-swell wave height of 7 m and to O(0.1) m at the most southern location at a water depth of 34 m with a significant sea-swell wave height of 5 m. These appreciable IG wave heights in O(30) m water depth represent a lower bound for the expected maximum IG wave heights given the fact that in the present analysis only a fraction of the full IG frequency range is considered. Comparisons with the predicted bound IG waves show that these can contribute substantially to the observed total IG wave height during storm conditions. The ratio between the predicted bound-and observed total IG variance ranges from 10% to 100% depending on the location of the observations and the timing during the storm. The ratio is typically high at the peak of the storm and is lower at both the onset and waning of the storm. There i, Environmental Fluid Mechanics, Coastal Engineering

High quality laboratory measurements of nearshore waves and morphology change at, or near prototype-scale are essential to support new understanding of coastal processes and enable the development and validation of predictive models. The DynaRev experiment was completed at the GWK large wave flume over 8 weeks during 2017 to investigate the response of a sandy beach to water level rise and varying wave conditions with and without a dynamic cobble berm revetment, as well as the resilience of the revetment itself. A large array of instrumentation was used throughout the experiment to capture: (1) wave transformation from intermediate water depths to the runup limit at high spatio-temporal resolution, (2) beach profile change including wave-by-wave changes in the swash zone, (3) detailed hydro and morphodynamic measurements around a developing and a translating sandbar., Corrigendum DOI: 10.1038/s41597-021-00874-2, Environmental Fluid Mechanics

Tide is influenced due to not only mainly tide generating force but also local wind and weather patterns. The East Asian monsoons cause strong seasonal climatic variations in the Mekong Delta. A two-dimensional, barotropic numerical model was employed to investigate the dynamics of tidal wave propagation in the South China Sea with a particular interest for its characteristics along the Mekong deltaic coast under wind monsoon climate. The results reveal that wind monsoon climate could causes damped or amplified tidal amplitudes around Mekong deltaic coast approximately 2–3 cm due to the changing atmospheric pressure, the tangential stress of wind over the water surface, and wind enhanced bottom friction. The monsoon climate influences rather strongly on the M2 semidiurnal tide system in the eastern Mekong deltaic coast, meanwhile the monsoon climate controls K1 diurnal tide in the western region of Mekong delta., Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Coastal Engineering, Environmental Fluid Mechanics

Cheniers are important for stabilising mud-dominated coastlines. A chenier is a body of wave-reworked, coarse-grained sediment consisting of sand and shells overlying a muddy substrate. In this paper we present and analyse a week of field observations of the dynamics of a single chenier along the coast of Demak, Indonesia. Despite relatively calm hydrodynamics during the one-week observational period, the chenier migrated surprisingly fast in the landward direction. The role of the tide and waves on the cross-shore chenier dynamics is explored using velocity moments as a proxy for the sediment transport. This approach shows that both tide and waves are capable of transporting the sediment of the chenier system. During calm conditions (representative for the south-east monsoon season), the tides generate a landward-directed sediment transport when the chenier crest is high relative to mean sea level. Waves only generate substantial sediment transport (direct, via skewness, and indirect, via stirring) when the chenier is submerged during periods with higher waves. The cross-shore chenier dynamics are very sensitive to the timing of tide and waves: most transport takes place when high water levels coincide with (relatively) high waves., Environmental Fluid Mechanics

Intertidal bars are naturally occurring morphological features along the waterline of sandy beaches. Present quantitative knowledge on intertidal bar behavior is limited, due to the scarcity of data resources and the limitations of traditional survey techniques. To investigate and quantify the cross-shore morphologic behavior of intertidal bars, hourly terrestrial laser scans of Kijkduin beach (The Netherlands) are used and a conceptual evolution intertidal bar model is constructed. In a six-week period in January and February 2017, a pronounced intertidal bar formed at Kijkduin beach and migrated onshore during mild wave conditions and eroded again during storm conditions. The observed maximum shoreward migration was 30 m horizontally with a maximum growth of about 1 m in the vertical direction. Onshore sediment transport fluxes peaked around 2 m3 per m width per day. In the conceptual model proposed here, run-up and overwash processes are dominant for shoreward growth and migration of the bar and submersion processes are responsible for bar destruction., Coastal Engineering, Environmental Fluid Mechanics

A large-scale field campaign was carried out on the ebb-tidal delta (ETD) of Ameland Inlet, a basin of the Wadden Sea in the Netherlands, as well as on three transects along the Dutch lower shoreface. The data have been obtained over the years 2017-2018. The most intensive campaign at the ETD of Ameland Inlet was in September 2017. With this campaign, as part of KustGenese2.0 (Coastal Genesis 2.0) and SEAWAD, we aim to gain new knowledge on the processes driving sediment transport and benthic species distribution in such a dynamic environment. These new insights will ultimately help the development of optimal strategies to nourish the Dutch coastal zone in order to prevent coastal erosion and keep up with sea level rise. The dataset obtained from the field campaign consists of (i) single-and multi-beam bathymetry; (ii) pressure, water velocity, wave statistics, turbidity, conductivity, temperature, and bedform morphology on the shoal; (iii) pressure and velocity at six back-barrier locations; (iv) bed composition and macrobenthic species from box cores and vibrocores; (v) discharge measurements through the inlet; (vi) depth and velocity from X-band radar; and (vii) meteorological data. The combination of all these measurements at the same time makes this dataset unique and enables us to investigate the interactions between sediment transport, hydrodynamics, morphology and the benthic ecosystem in more detail. The data provide opportunities to calibrate numerical models to a high level of detail. Furthermore, the open-source datasets can be used for system comparison studies. The data are publicly available at 4TU Centre for Research Data at https://doi.org/10.4121/collection:seawad (Delft University of Technology et al., 2019) and https://doi.org/10.4121/collection:kustgenese2 (Rijkswaterstaat and Deltares, 2019). The datasets are published in netCDF format and follow conventions for CF (Climate and Forecast) metadata. The http://data.4tu.nl (last access: 11 Nove, Environmental Fluid Mechanics, Coastal Engineering, Management Support

Sediment transport over intertidal flats is driven by a combination of waves, tides, and wind-driven flow. In this study we aimed at identifying and quantifying the interactions between these processes. A five week long dataset consisting of flow velocities, waves, water depths, suspended sediment concentrations, and bed level changes was collected at two locations across a tidal flat in the Wadden Sea (The Netherlands). A momentum balance was evaluated, based on field data, for windy and non-windy conditions. The results show that wind speed and direction have large impacts on the net flow, and that even moderate wind can reverse the tidal flow. A simple analytical tide–wind interaction model shows that the wind-induced reversal can be predicted as a function of tidal flow amplitude and wind forcing. Asymmetries in sediment transport are not only related to the tide–wind interaction, but also to the intratidal asymmetries in sediment concentration. These asymmetries are influenced by wind-induced circulation interacting with the large scale topography. An analysis of the shear stresses induced by waves and currents revealed the relative contributions of local processes (resuspension) and large-scale processes (advection) at different tidal flat elevations., Environmental Fluid Mechanics, Coastal Engineering

Considering the likely increase in coastal flooding in small island developing states (SIDSs) due to climate change, coastal managers at the local and global levels have been developing initiatives aimed at implementing disaster risk reduction (DRR) and adaptation measures. Developing science-based adaptation policies requires accurate coastal flood risk (CFR) assessments, which in the case of insular states are often subject to input uncertainty. We analysed the impact of a number of uncertain inputs on coastal flood damage estimates: (i) significant wave height, (ii) storm surge level and (iii) sea level rise (SLR) contributions to extreme sea levels, as well as the error-driven uncertainty in (iv) bathymetric and (v) topographic datasets, (vi) damage models, and (vii) socioeconomic changes. The methodology was tested through a sensitivity analysis using an ensemble of hydrodynamic models (XBeach and SFINCS) coupled with a direct impact model (Delft-FIAT) for a case study of a number of villages on the islands of São Tomé and Príncipe. Model results indicate that for the current time horizon, depth damage functions (DDFs) and digital elevation models (DEMs) dominate the overall damage estimation uncertainty. When introducing climate and socioeconomic uncertainties to the analysis, SLR projections become the most relevant input for the year 2100 (followed by DEM and DDF). In general, the scarcity of reliable input data leads to considerable predictive uncertainty in CFR assessments in SIDSs. The findings of this research can help to prioritize the allocation of limited resources towards the acquisitions of the most relevant input data for reliable impact estimation., Coastal Engineering, Hydraulic Structures and Flood Risk, Environmental Fluid Mechanics

The nonlinear wave shape, expressed by skewness and asymmetry, can be calculated from surface elevation or pressure time series using bispectral analysis. Here, it is shown that the same analysis technique can be used to calculate the bound superharmonic wave height. Using measured near-bed pressures from three different field experiments, it is demonstrated that there is a clear relationship between this bound wave height and the nonlinear wave shape, independent of the measurement time and location. This implies that knowledge on the spatially varying bound wave height can be used to improve wave shape-induced sediment transport predictions. Given the frequency-directional sea-swell wave spectrum, the bound wave height can be predicted using second order wave theory. This paper shows that in relatively deep water, where conditions are not too nonlinear, this theory can accurately predict the bispectrally estimated bound superharmonic wave height. However, in relatively shallow water, the mismatch between observed and predicted bound wave height increases significantly due to wave breaking, strong currents, and increased wave nonlinearity. These processes are often included in phase-averaged wind-wave models that predict the evolution of the frequency-directional spectrum over variable bathymetry through source terms in a wave action balance, including the transfer of energy to bound super harmonics. The possibility to calculate and compare with the observed bound super harmonic wave height opens the door to improved model predictions of the bound wave height, nonlinear wave shape and associated sediment transport in large-scale morphodynamic models at low additional computational cost., Environmental Fluid Mechanics

Although commonly used for the validation of morphological predictions, point-wise accuracy metrics, such as the root-mean-squared error (RMSE), are not well suited to demonstrate the quality of a high-variability prediction; in the presence of (often inevitable) location errors, the comparison of depth values per grid point tends to favour predictions that underestimate variability. In order to overcome this limitation, this paper presents a novel diagnostic tool that defines the distance between predicted and observed morphological fields in terms of an optimal sediment transport field, which moves the misplaced sediment from the predicted to the observed morphology. This optimal corrective transport field has the “cheapest” quadratic transportation cost and is relatively easily found through a parameter-free and symmetric solution procedure solving an elliptic partial differential equation. Our method, which we named effective transport difference (ETD), is a variation to a partial differential equation approach to the Monge–Kantorovich L2 optimal transport problem. As a new error metric, we propose the root-mean-squared transport error (RMSTE) as the root-mean-squared value of the optimal transport field. We illustrate the advantages of the RMSTE for simple 1D and 2D cases as well as for more realistic morphological fields, generated with Delft3D, for an idealized case of a tidal inlet developing from an initially highly schematized geometry. The results show that by accounting for the spatial structure of morphological fields, the RMSTE, as opposed to the RMSE, is able to discriminate between predictions that differ in the misplacement distance of predicted morphological features, and avoids the consistent favouring of the underprediction of morphological variability that the RMSE is prone to., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Coastal Engineering, Environmental Fluid Mechanics

Many coral reef-lined coasts are low-lying with elevations <4 m above mean sea level. Climate-change-driven sea-level rise, coral reef degradation, and changes in storm wave climate will lead to greater occurrence and impacts of wave-driven flooding. This poses a significant threat to their coastal communities. While greatly at risk, the complex hydrodynamics and bathymetry of reef-lined coasts make flood risk assessment and prediction costly and difficult. Here we use a large (>30,000) dataset of measured coral reef topobathymetric cross-shore profiles, statistics, machine learning, and numerical modeling to develop a set of representative cluster profiles (RCPs) that can be used to accurately represent the shoreline hydrodynamics of a large variety of coral reef-lined coasts around the globe. In two stages, the large dataset is reduced by clustering cross-shore profiles based on morphology and hydrodynamic response to typical wind and swell wave conditions. By representing a large variety of coral reef morphologies with a reduced number of RCPs, a computationally feasible number of numerical model simulations can be done to obtain wave runup estimates, including setup at the shoreline and swash separated into infragravity and sea-swell components, of the entire dataset. The predictive capability of the RCPs is tested against 5,000 profiles from the dataset. The wave runup is predicted with a mean error of 9.7–13.1%, depending on the number of cluster profiles used, ranging from 312 to 50. The RCPs identified here can be combined with probabilistic tools that can provide an enhanced prediction given a multivariate wave and water level climate and reef ecology state. Such a tool can be used for climate change impact assessments and studying the effectiveness of reef restoration projects, as well as for the provision of coastal flood predictions in a simplified (global) early warning system., Environmental Fluid Mechanics, Coastal Engineering

Wave forecasting in ocean and coastal waters commonly relies on spectral models based on the spectral action balance equation. These models assume that different wave components are statistically independent and as a consequence cannot resolve wave interference due to statistical correlation between crossing waves, as may be found in, for instance, a focal zone. This study proposes a statistical model for the evolution of wave fields over non-uniform currents and bathymetry that retains the information on the correlation between different wave components. To this end, the quasi-coherent model (Smit & Janssen, J. Phys. Oceanogr., vol. 43, 2013, pp. 1741-1758) is extended to allow for wave-current interactions. The outcome is a generalized action balance model that predicts the evolution of the wave statistics over variable media, while preserving the effect of wave interferences. Two classical examples of wave-current interaction are considered to demonstrate the statistical contribution of wave interferences: (1) swell field propagation over a jet-like current and (2) the interaction of swell waves with a vortex ring. In both examples cross-correlation terms lead to development of prominent interference structures, which significantly change the wave statistics. Comparison with results of the SWAN model demonstrates that retention of cross-correlation terms is essential for accurate prediction of wave statistics in shear-current-induced focal zones., Environmental Fluid Mechanics

Waves are important drivers for reef hydrodynamics, and therefore strongly contribute to flooding over reef-lined coasts. While high-frequency waves are largely dissipated when they propagate over the reef flat due to breaking and friction, low-frequency (LF) waves are generally able to reach the back-reef beach. There, they can reflect and form (quasi-) standing wave patterns, which under resonant conditions can lead to disproportionally high run-up on the beach (e.g., Pequignet et al., 2009; Gawehn et al., 2016). The probability of this phenomenon is expected to increase due to sea-level rise (e.g., Pequignet et al., 2009). In this study, we numerically investigate long wave resonance and the processes enhancing or limiting the resonant amplification of long waves over the reef flat. Besides the role of frictional dissipation (e.g., Pomeroy et al. 2012), we investigate how the nonlinear transformation of long waves influences the amplification rate., Environmental Fluid Mechanics, Coastal Engineering, Hydraulic Structures and Flood Risk