Beach groundwater and nearshore hydrodynamic data were collected during a field experiment along two dissipative beach transects on Galveston Island, Texas, in the fall of 2023. The monitored beaches serve as nesting habitat for the critically endangered Kemp’s ridley sea turtle. Conditions ranged from calm to stormy, with two storms occurring during the experiment, inundating the entire beach up to the dune toe. Collected hydrodynamic data include readings from pressure loggers submerged in the foreshore and mounted in groundwater wells in the backshore, data from two wave buoys about 1.5 km offshore, and GoPro timestacks of the instantaneous waterline (wave runup). Other collected data include bathymetry and topography surveys, subsurface temperature and moisture content readings, and sediment characteristics. This comprehensive dataset can be used to (1) study relevant beach inundation and groundwater processes, including their effect on the local ecosystem (e.g., repeated flooding of sea turtle nests), (2) study the propagation of nearshore hydrodynamic processes into the beach matrix and groundwater table, and (3) validate existing beach groundwater models.

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.

The swash zone is an important region for the coastal morphodynamics. Often, model studies of the swash zone use depth-averaged models. These models typically assume a vertically uniform velocity and sand concentration for calculating the sand transport flux. However, this assumption is not always accurate in the swash zone. In order to investigate the vertical distribution of velocity and sand, we use a depth-resolving model that is able to capture these vertical variations. We simulate the flow and suspended sediment transport induced by bichromatic waves using a 2DV depth-resolving RANS model. Our verification of the model shows that special care needs to be taken to deal with bubbles in 2DV simulations. Furthermore, we show that turning off the (Wilcox, 2006, 2008) limiter for turbulence, increases the modelled turbulent kinetic energy that is induced by wave-breaking, resulting in improved predictions of sediment concentrations. Using the depth-resolving model, we show that the vertical distribution of velocity and sand is far from uniform in the swash zone. The results show that if one assumes vertically uniform depth-averaged velocities and concentrations, one can overpredict the sediment flux by 50%.

Aeolian sand strips are bedforms that often form on the beach during high wind events. The sand strips can move across the beach resulting in sediment transport, but these sediment fluxes have not yet been quantified. Permanent laser scanning (PLS) provides an opportunity to gather long-term and detailed information about sand strips. Here, we present a method that detects sand strips in PLS reflectance images and extracts their wavelength, orientation, and height. The detection method is used on a month-long hourly dataset of Noordwijk beach, the Netherlands. The average sand strip wavelength (13 m), height (3 cm), and orientation (onshore-oblique) found in this dataset are consistent with previous research. Migration rates are obtained by examining successive images of sand strip patterns. Using the concurrent sand strip height, the sediment transport rate associated with sand strips is estimated. During a single 8-hour event, a total sediment flux of 0.16 m3/m is observed of which 0.07 m3/m in onshore direction. Furthermore, the detection method is shown to be applicable to other beaches by calibrating its thresholds. This opens a new possibility for future research to use PLS to investigate the behavior and sediment transport rates of aeolian sand strips.

In spectral wave models, the nonlinear triad source term accounts for the transfer of energy to the bound higher harmonics. This paper presents an extension to commonly used spectral models that resolves the evolution of the bound wave energy by keeping track of the energy that has been bound by the triad interactions. This extension is referred to as the bound wave evolution (BWE) model. From this, the spatial evolution of the bound wave height is obtained, which serves as a proxy for the nonlinear wave shape. The accuracy of these bound wave heights, and thus wave shape predictions, is highly dependent on the accuracy of the triad source term. Therefore, in this study, the capability of the LTA and SPB triad formulations to capture the growth of the bound wave height is evaluated. For both of these formulations, it is found that slope dependent calibration parameters are required. Overall, despite being computationally more expensive, the SPB method proves to be significantly more accurate in predicting the bound wave evolution. In the shoaling zone, where the bound wave energy is dominated by triads, the BWE model is well capable of predicting the nonlinear wave’s shape. In the surf zone, however, where a combination of triads and wave breaking control the spectral evolution, the BWE model over-predicts the bound wave height. Nevertheless, this paper shows the promising capabilities of spectral models to predict the nonlinear wave shape.

During extreme conditions, the transport of the wave-averaged suspended sediment concentrations in the inner surf zone affects dune erosion. Although large-scale laboratory experiments have provided insight in what drives these sediment concentrations, corresponding field data are lacking. To fill this gap, novel field observations of suspended sediment concentrations are compared to drivers that govern sediment suspension during storm conditions known from literature. A total of 128 time intervals of 20 min are analysed, spread over 10 different high water events with different hydrodynamic conditions. For each time interval, the wave-averaged (i.e. 20 min mean) suspended sediment concentration is computed and compared to three suspension drivers. The studied drivers are (1) bed shear due to near bed velocities that originate from mean currents in combination with wave-induced orbital flow, (2) the horizontal pressure gradients under steep wave fronts that increase the forces on the bed material, and (3) bore-induced turbulence that is generated at the free surface and reaches the bed. The derived bore-induced turbulence generates the greatest correlation with the mean suspended sediment concentrations (r = 0.74, p = 4.47E-23). Samples that deviate from this correlation correspond to time intervals with lower values of derived bore turbulence, less wave energy saturation in the inner surf zone, and stronger mean currents. The correlation with the mean suspended sediment concentrations increases when the shear stress originating from mean currents is used for these time intervals (r = 0.83, p = 1.63E-33). For time intervals during which more energetic conditions persist and the wave energy is saturated in the nearshore, bore turbulence was the dominant mechanism in stirring up sediment. The outcome of this study suggests that, based on the events analysed, dune erosion models may achieve more accurate results if computations of suspended sediment concentrations include a bore-induced turbulence term, or if already included, properly address the relative importance of bore-induced turbulence when compared to bed shearing.

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.

High-resolution wave measurements at intermediate water depth are required to improve coastal impact modeling. Specifically, such data sets are desired to calibrate and validate models, and broaden the insight on the boundary conditions that force models. Here, we present a wave data set collected in the North Sea at three stations in intermediate water depth (6–14 m) during the 2021/2022 storm season as part of the RealDune/REFLEX experiments. Continuous measurements of synchronized surface elevation, velocity and pressure were recorded at 2–4 Hz by Acoustic Doppler Profilers and an Acoustic Doppler Velocimeter for a 5-month duration. Time series were quality-controlled, directional-frequency energy spectra were calculated and common bulk parameters were derived. Measured wave conditions vary from calm to energetic with 0.1–5.0 m sea-swell wave height, 5–16 s mean wave period and W-NNW direction. Nine storms, i.e., wave height beyond 2.5 m for at least six hours, were recorded including the triple storms Dudley, Eunice and Franklin. This unique data set can be used to investigate wave transformation, wave nonlinearity and wave directionality for higher and lower frequencies (e.g., sea-swell and infragravity waves) to compare with theoretical and empirical descriptions. Furthermore, the data can serve to force, calibrate and validate models during storm conditions.

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).

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.

Wave nonlinearity plays an important role in cross-shore beach morphodynamics and is often parameterized in engineering-type morphodynamic models through a nonlinear relationship with the Ursell number. It is not evident that the relationship established in previous studies also holds for sheltered sites with fetch-limited seas as they are more prone to effects of local winds and currents, the waves are generally steeper, and the beaches are typically reflective. This study investigates near-bed orbital velocity nonlinearity from wave records collected at two sheltered beaches in The Netherlands and contrasts them to earlier observations made along the exposed, wave-dominated North Sea coast. Our observations at sheltered beaches show that the Ursell number has comparable skill in predicting wave nonlinearity as it has on previously studied exposed coasts. However, the orbital velocities at sheltered coasts are more asymmetric for the same Ursell number than on exposed coasts. When exposed coast data were examined for moments with comparable high-steepness waves, a similar effect on asymmetry was observed. In addition, following and opposing winds were found to have a clear relationship with total nonlinearity, while they did not affect the phase between skewness and asymmetry at the sheltered beaches. Refitting the free parameters of an Ursell-based predictor improved the bias for the asymmetry parameterization. Whether this has implications for modeling of the magnitude of wave-nonlinearity-driven sediment transport using engineering type models is strongly dependent on the sediment transport formulation used, as these formulations depend on additional calibration coefficients too.

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.

Traditionally, many numerical studies of sediment transport or morphodynamics in the swash zone use depth-averaged approach. This is relatively computationally cheap, however, the assumptions it makes on vertical structures of the flow velocity and suspended sediment concentration are not always fully met. Depth-resolving models, albeit more computationally expensive, offer the ability to investigate sediment transport without imposing assumptions on vertical structures. Therefore, these models are a great tool for investigating processes related to sediment transport in the swash zone. In this study, will look at the behavior of suspended sediment in the swash zone, and the role vertical structures therein.

Coastal Sediments 2023, pp. 1923-1932 (2023) No AccessSHAPING THE BEACH: CROSS-SHORE SAND TRANSPORT IN THE SWASH ZONEJEBBE VAN DER WERF, SARA DIONÍSIO ANTÓNIO, JOOST KRANENBORG, WEIQIU CHEN, GEERT CAMPMANS, ERIK HORSTMAN, AD RENIERS, and SUZANNE HULSCHERJEBBE VAN DER WERFDepartment of Marine and Fluvial Systems, University of Twente, Drienerlolaan 5, 7522 NB Enschede, The NetherlandsUnit of Marine & Coastal Systems, Deltares, P.O. Box 177, 2600, MH, Delft, The Netherlands, SARA DIONÍSIO ANTÓNIODepartment of Marine and Fluvial Systems, University of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands, JOOST KRANENBORGDepartment of Marine and Fluvial Systems, University of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands, WEIQIU CHENDepartment of Marine and Fluvial Systems, University of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands, GEERT CAMPMANSDepartment of Marine and Fluvial Systems, University of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands, ERIK HORSTMANDepartment of Marine and Fluvial Systems, University of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands, AD RENIERSFaculty of Civil Engineering and Geosciences, Delft University of Technology, P.O. Box 5048, 2600 GA Delft, The Netherlands, and SUZANNE HULSCHERDepartment of Marine and Fluvial Systems, University of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlandshttps://doi.org/10.1142/9789811275135_0176Cited by:0 (Source: Crossref) PreviousNext AboutSectionsPDF/EPUB ToolsAdd to favoritesDownload CitationsTrack CitationsRecommend to Library ShareShare onFacebookTwitterLinked InRedditEmail Abstract: This paper summarizes the main results from the three main activities of Shaping The Beach research project (2018-2023), and gives recommendations for further studies into swash zone sand transport. The project involved a unique and integral combination of detailed wave flume experiments with advanced measuring techniques, sophisticated intra-swash numerical modelling, and engineering, practical sand transport modeling. Shoreline retreat/advance is driven by high-/low-energy wave conditions, independent of the beach slope. Milder beach slope does decrease the shoreline retreat and support more accretion. Both effects are captured by relating beach volume change to the slope-corrected Dean number as proposed by Hattori and Kawamata (1980). Depth-resolving VoF RANS models like OpenFoam are key to unravel the vertical flow and suspended sand concentration structures within the swash zone. These models only work well if air-flow interactions and turbulence are carefully represented. Existing practical models have skill in predicting swash zone transport rates and bed level changes. However, they do not perform equally well for erosive and accretive conditions and require heavy data fitting, which makes them not fully predictive (yet). Knowledge on swash zone sand transport is essential to understand and cope with coastal erosion and sedimentation by natural processes and human interferences. This knowledge is used to assess climate change impacts on beach stability, to design cost-effective beach protection measures such as sand nourishments, and to design coastal structures. FiguresReferencesRelatedDetails Recommended Coastal Sediments 2023Metrics History PDF download

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.

The generation and propagation of infragravity waves at sandy coasts has received significant attention (e.g. Herbers et al., 1994). This is in contrast with tidal inlets where both observations and modeling studies are scarce. Here we present novel observations and modeling of infragravity waves over an ebb-tidal delta located at Ameland, the Netherlands. The objective is to understand the evolution of the free and bound infragravity waves as they propagate from the deeper shelf onto an ebb-tidal shoal and into the ebbchannel, establish their pathways and examine the potential in sediment transport. The latter is important in view of the anticipated large-scale ebb-tidal delta nourishment at this location to protect the coast from future storm impacts.

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.

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.

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.