Sea turtles are key species in many coastal ecosystems worldwide, particularly coral reef and seagrass habitats. Yet, six of seven species are endangered. Their nests, which incubate in beach sand and rely on specific climatic conditions for egg viability, face significant threats from inundation, for example through wave runup. This paper examines a method to rapidly predict wave runup in low-data coral reef environments, and the implications thereof on the inundation of sea turtle nests. The study uses two metamodels, BEWARE-2 and HyCReWW, to predict wave runup at Ras Baridi, Saudi Arabia, a key nesting site of the Red Sea green turtle population. The models were used to analyze runup events and inundation durations and provide a first estimate of a safe nesting elevation. Despite data limitations, the study provides valuable insights for coastal managers to protect sea turtle nests, suggesting that a 5-year return period runup elevation could serve as a threshold for nest relocation. However, the findings also highlight the importance of more accurate hydrodynamic predictions and the need for in-situ data to validate models and improve conservation strategies.

Threatened sea turtles rely on sandy beaches for nesting, linking their long-term survival to global beach availability. However, beaches worldwide are increasingly threatened by anthropogenic stressors and sea level rise (SLR). Reliable vulnerability assessments require understanding beach dynamics across multiple time scales, informed by long-term coastal change records. While many nesting beaches lie in remote, data-poor environments, recent advances in coastal remote sensing now allow us to monitor coastal change worldwide. Here, we combine satellite-derived shorelines (CoastSat), shoreline modeling (CoSMoS-COAST), and global data sets to investigate shoreline evolution and future vulnerability at nine globally important sea turtle nesting sites. We investigate seasonal and long-term shoreline change, hindcast (1980–2024) and forecast (2025–2100) shoreline positions under various SLR scenarios, and quantify available accommodation space based on backbeach elevation and infrastructure footprints. We find that shoreline evolution and vulnerability vary considerably, with three sites showing historical accretion trends and four sites showing erosion. This demonstrates that the previously widely applied bathtub approach—adding SLR to a static beach profile—is not suitable to assess the vulnerability of sea turtle nesting beaches to erosion. Three eroding beaches emerge as particularly vulnerable due to projected shoreline retreat coupled with limited accommodation space. Despite significant uncertainties arising from long-term shoreline projections, our results provide important insights into seasonal and long-term morphodynamics, identify vulnerable nesting sites, and offer a comprehensive, transferable framework for assessing shoreline evolution and relative erosion vulnerability at other sites. Understanding these dynamics is crucial to inform conservation and management strategies to future-proof these critical nesting habitats.

Plain Language Summary
Sea turtles depend on sandy beaches for nesting, which means their survival is closely linked to how these beaches change over time. Today, many beaches are increasingly pressured by human activity and rising sea levels, putting turtle nesting habitats at risk. To better understand which beaches are most vulnerable, we used satellite images, computer models, and global data to study nine of the world's most important nesting sites. We looked at how the shoreline has moved since 1980, how it might change through 2100 under different sea level rise (SLR) scenarios, and how much space may remain for turtles to nest given local terrain and development. Our results show that some beaches are naturally building up while others are eroding, and that vulnerability is not the same everywhere. In particular, three beaches appear especially at risk because they are eroding and have little room for turtles to nest further inland. These findings highlight the importance of moving beyond simple “bathtub” estimates of SLR, and instead considering the complex, long-term behavior of beaches. This approach can help identify priority sites for conservation and guide strategies to protect sea turtle nesting habitats in a changing world.

Accurate prediction of Wave-Group-Forced (WGF) InfraGravity (IG) waves depends on resolving the corresponding phase shift, typically achieved through a coupled phase – amplitude equation. However, this approach requires a grid resolution that resolves the structure of the wave groups making it computationally expensive at regional scales. To address this limitation, an existing local expression for the phase shift of normally incident WGF-IG waves has been extended to account for directional seas. The extended formulation is verified against predictions from the coupled phase – amplitude model using bichromatic wave forcing over a uniformly sloping beach for a wide range of sea-swell conditions. Results show that the local approach performs well in the off-resonant region for obliquely incident waves. When applied outside this regime, however, its accuracy decreases, with performance varying depending on sea-swell and bathymetric conditions. The coupled and local phase shift approaches are also validated with observations obtained during the Coast3D field experiment. The total, incoming and outgoing IG waves are predicted with comparable skill and root mean square error for both methods. The good match using the local expression is attributed to the fact that the conditions during Coast3D correspond to directionally broad sea-swell spectra with relative short peak periods propagating over moderately sloping bathymetry for which the verification showed significant skill. Additional validation with field observations at other locations are necessary to firmly determine the limitations of the use of a local phase shift.

Hybrid dune-dike structures are innovative developments creating coastal defense systems which are more conveniently integrated with the natural environment. In this study, a numerical study was conducted to investigate the temporal evolution of wave overtopping, with the changing profile of the dune under extreme storm conditions with a constant water level, of two types of hybrid dune-dike structures in Katwijk (dike-in-dune type) and Raversijde (dune-in-front-of-dike type). XBeach 1DH was used to firstly calculate bed profiles for different time steps during a 10-h storm duration using the Surfbeat mode and then, in a second step, mean wave overtopping rates were modelled for each calculated bed profile using the Non-hydrostatic mode. According to the simulation results, most of the dune erosion occurs during the first two hours of the storm, and then continues at a slower rate as the sand deposits in front of the dune. Once the hybrid structure is eroding (so for t > 0), the significant wave height at the dike toe and the mean overtopping discharge increase in time for both Katwijk and Raversijde, although it quickly reaches a plateau for Raversijde. The first simulations with the original non-eroded profiles deviate from this trend. The reason for this deviation needs to be further investigated.

Tracking coastal sediments can provide useful information about coastal dynamics, thereby helping coastal management. However, the highly dynamic conditions of the coasts makes analyzing the trajectories of a huge number of particles challenging. To solve this limitation, the framework of coastal sediment connectivity is designed. In this framework, recent advances in graph theory are used to quantify coastal systems as complex networks. In this context, sediment sinks/sources and pathways represent the graphical nodes and links, respectively. In this work, we take the first step to evaluate the ability of this newly-developed framework in quantifying the basic processes on a sandy beach. Firstly, we used Delft3D to obtain the velocity field and bed-level changes. Then, the Eulerian results were fed into SedTRAILS to simulate the sediment pathways. We show that the current version of the model can correctly calculate the basic metrics of the sediment-connectivity network (e.g., network link strength which is a proxy for sediment fluxes). More specifically, we show that this framework is capable of exploring the initiation of the rip channel formation.

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.

This study examines the importance of free infragravity (FIG) waves in the North Sea using a recent collection of wave measurements and a newly developed unstructured SWAN model. The measurements include new observations of infragravity waves at offshore (30–40 m water depth) and nearshore (10–20 m water depth) locations in the southern North Sea. These observations serve as the basis for model optimization and verification. Good agreement is obtained between model predictions and measurements during two recent storm periods, including severe storms with unusual wind directions and high wind speeds (e.g., “Storm Babet”). Model investigation along the coasts of Belgium and the Netherlands demonstrated a strong dependence between nearshore FIG conditions (i.e., energy intensity and sources) and storm characteristics (i.e., alongshore wind pattern and storm track). Specifically, several storms have demonstrated significant contributions of FIG energy originating from remote sources (e.g., the coasts of UK and Denmark). This suggests that nearshore FIG conditions in the North Sea cannot be determined based on the local sea-swell conditions alone and may be significantly underestimated if non-local contributions are ignored. Finally, modelled and measured results at nearshore locations along the Dutch coast revealed that under storm conditions FIG energy can be an order of magnitude higher than energy due to bound infragravity (BIG) waves. This result, augmented with estimated ratios of free and forced infragravity energy at the shoreline, emphasizes the necessity of considering the FIG waves as an integral part of coastal safety assessments along the coasts of the North Sea.

Temporary permeable structures of bamboo and brushwood have been implemented for mangrove restoration along retreating coastlines worldwide. However, deriving lessons from previous studies is inhibited by their lack of morphodynamic context, with missing bathymetric data or control comparisons. In this paper, we present a low-tech, low-cost, data collection methodology to support morphodynamic system understanding and modeling of mangrove coastlines. This method was applied to monitor a mangrove restoration project featuring temporary permeable structures of bamboo and PVC, installed in late 2021 on the subsiding muddy coast of Demak, Indonesia. Seabed level changes were regularly tracked with bathymetric surveys and monitoring poles across structures and at a nearby control site. Structures were positioned landward of a chenier, at −0.7 m to −0.9 m relative to mean sea level (MSL), and 30–70 m seaward of the mangrove fringe. Measurements from August 2021 to December 2022 revealed seabed erosion (−0.33 m to −0.4 m) seaward of the structures, with mixed responses landward: two sections eroded (−0.04 m to −0.05 m), one remained stable, and a creek-adjacent section eroded by −0.43 m. At the nearby control site, chenier migration and vertical growth promoted landward accretion, though elevations remained below MSL and thus unsuitable for mangrove colonization. The bathymetric and monitoring pole measurements presented in this study constitute valuable datasets for modeling studies aiming to unravel the dominant processes driving morphodynamic changes. Such models could also inform integrated approaches to mangrove restoration in subsiding coastlines, considering sediment supply, subsidence management, and structure integrity.

Beach groundwater dynamics play a critical role in coastal ecosystem functions, particularly in low-lying beach habitats used for nesting by endangered species like sea turtles. Incubating nests are susceptible to prolonged inundation below the groundwater table (GWT), as flooding duration critically affects egg viability. Understanding how oceanic processes and rain drive GWT fluctuations in the nesting area is essential for evaluating nest relocation strategies and designing nature-based solutions that mitigate nest flooding. Here, we analyze how infragravity waves, tides, storm surge, and rainfall drive short-term fluctuations (hourly to weekly) in the beach GWT on Galveston Island, Texas—a dissipative, mild-sloping barrier island system along the northwestern Gulf of Mexico coast. Applying tailored spectral analyses to field observations collected in 2023, we show that surge and rainfall dominate short-term GWT response in the nesting area, while higher-frequency wave and tidal signals are increasingly damped landward. To facilitate this analysis, we classify observed water levels into groundwater, mixed, and submerged regimes based on estimated wave runup. A flooding threshold analysis reveals multiple, prolonged nest inundation events (exceeding 12 h) across the backshore, even for the shallowest nests. This strongly suggests that Galveston Island’s beaches are currently unsuitable for turtle nesting, underscoring the need to continue the ongoing nest relocation program and further research nature-based solutions that enable sea turtle nesting (e.g., turtle-friendly nourishments).

Wave runup observations are important for coastal management providing data to validate predictive models of inundation frequencies and erosion rates, which are vital for assessing the vulnerability of coastal ecosystems and infrastructure. Automated algorithms to extract the instantaneous water line from video imagery struggle under dissipative conditions, where the presence of a seepage face and the lack of contrast between the sand and the swash impede proper extraction, requiring time-intensive data quality control or manual digitization. This study introduces two novel methods, based on color contrast (CC) and machine learning (ML). The CC method combines texture roughness — local entropy — with saturation. Images are first binarized using entropy values and then refined through noise reduction by binarization of the saturation channel. The ML method uses a convolutional neural network (CNN) informed by five channels: the grayscale intensity and its time gradient, the saturation channel, and the entropy and its time gradient. Both methods were validated against nine manually labeled, 80 min video time series. The CC method demonstrated strong agreement with manually digitized water lines (RMSE = 0.12 m, r=0.94 for the vertical runup time series; RMSE = 0.08 m, r=0.97 for the 2% runup exceedance (R_2%); and RMSE = 3.88 s, r=0.70 for the mean period (T_m−1,0)). The ML model compared well with the manually labeled time series (RMSE = 0.10 m, r=0.96 for the vertical runup time series; RMSE = 0.09 m, r=0.97 for R_2%; and RMSE = 3.51 s, r=0.79 for T_m−1,0). Furthermore, the computed R_2% values of both methods show a good agreement with the formula proposed by Stockdon et al. (2006) for extremely dissipative conditions, with RMSE-values lower than 0.13 m and correlations exceeding 0.70 for manual, CC, and ML estimates. While the CC method is deemed applicable for wave-by-wave analysis under similar dissipative conditions with a smooth seepage face and sufficient turbulent swash, the ML method still struggles with new, unseen data. However, it shows promise for a broader application and serves as a viable proof of concept. Together, these methods reduce the need for manual processing and enhance real-time coastal monitoring, contributing to more accurate predictive modeling of runup events and a better understanding of nearshore processes.

Climate change and human activity pose increasing challenges to endangered sea turtles, which are key species in many marine ecosystems worldwide. Among these challenges are the flooding and erosion of nesting beaches. In this perspective, we argue that existing methods and tools from coastal science and management hold significant, yet underused, potential for sea turtle conservation. We introduce a stepwise framework for integrating sea turtle ecology and coastal management to address these coastal threats. The framework follows an Observe–Understand–Predict–Intervene cycle and links ecological thresholds, coastal processes, and management interventions across scales, from Regional Management Units (RMUs) to individual beaches. We illustrate how state-of-the-art monitoring, modeling, and nature-based solutions (NBS) can be embedded within this framework to inform when and how to intervene. Increased in-situ data collection and interdisciplinary collaboration will be critical to apply and refine this approach, thereby enhancing the long-term resilience of nesting habitats.

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

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

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.

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.

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. Dataset: https://doi.org/10.4121/233f11ff-7804-4777-8b32-92c4606e56d8 Dataset License: CC-BY 4.0.

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.

Infragravity (IG) waves are relatively long waves with typical periods of several tens of seconds to several minutes. The energy at the IG band plays an important role in nearshore areas. For example, IG waves can significantly contribute to dune erosion and sediment transport (e.g., Roelvink et al., 2009), and may excite harbor oscillations (e.g., Bowers, 1977). Furthermore, IG waves may result in destructive inundation events (e.g., Roeber and Bricker, 2015). These documentations of IG waves' impacts emphasise the necessity to account for IG contributions as part of coastal hazard assessments, especially under storm conditions.

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.