Nature-based solutions are increasingly recognized as effective and multifunctional components of climate-resilient flood protection. While tropical mangroves have received substantial attention, temperate riparian forests, particularly willow systems, offer comparable wave attenuation and biodiversity benefits, yet remain understudied. This study assesses the ecological and protective value of three types of willow floodplain forests: a so-called wild-grown willow forest, a pollard willow forest, and a willow plantation. Using field data from the Biesbosch National Park (the Netherlands), we quantified forest structure, ground-dwelling invertebrate diversity, and modelled wave attenuation under storm scenarios. Structural complexity and biodiversity were highest in the wild-grown forest, with significantly greater invertebrate order richness, larger body sizes, and more heterogeneous canopy architecture. The pollard forest showed the highest wave attenuation efficiency due to their dense, low-lying crown structures. The plantation forest showed lower values across both axes. We integrated these findings into a trade-off model evaluating ecological value, flood protection efficiency, and a 50-year simple cost analysis of each forest type as a hybrid solution alongside traditional dikes. While the pollard forest is the most spatially efficient for flood attenuation, the wild-grown system provides greater ecological value at lower lifecycle cost. Our results underscore the importance of tailoring hybrid flood defense strategies to local priorities - balancing biodiversity, spatial constraints, and economic feasibility. The framework developed here can inform ecosystem-based design in delta regions worldwide, supporting integrated climate adaptation that aligns safety with ecological resilience.

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.

Wave overtopping can cause severe erosion on the crest and landward slope of a dike. Accurate erosion prediction requires resolving the spatiotemporal evolution of overtopping flow, which remains insufficiently understood. This study investigates the behavior of overtopping flow along the crest and landward slope using a new, efficient numerical model based on the Steep-Slope Shallow Water Equations. The model was validated against measurements from field experiments using a wave overtopping simulator. The model was applied to a typical dike geometry (Bc = 5 m; tan(α) = 1V:3H), where the overtopping flow is imposed at the waterside crest line based on a schematization using empirical equations from literature and insights from small-scale overtopping experiments. The spatiotemporal evolution of the flow was analyzed along the crest and landward slope for different overtopping volumes. The flow was observed to stretch in time along both the crest and slope while the wavefront steepens, causing the peak flow thickness (hpeak) to decrease significantly. The peak flow velocity (upeak) decreases along the crest but initially accelerates on the slope due to gravitational forcing. As the flow becomes progressively thinner and faster downslope, frictional forcing increases, reducing the acceleration. Eventually, gravitational and frictional forces balance, causing upeak to decelerate and then decrease. Overall, the model captures key spatiotemporal dynamics such as flow stretching, wavefront steepening, and deceleration of upeak on long slopes, which are absent in time-independent analytical models. It offers a computationally efficient approach that provides a practical middle ground between simplified analytical methods and full CFD simulations.

Salt marshes provide critical coastal protection by attenuating waves, yet their performance under extreme storm conditions and subsequent recovery have lacked quantitative assessment. This study quantifies temporal variations in wave attenuation capacity across a complete storm cycle (before-during-after) at a Scirpus mariqueter marsh in the Yangtze Estuary, and evaluates cumulative impacts under storm sequences through numerical modeling. Our results show: (1) Storm landfall led to significant reduction in wave attenuation capacity of salt marshes, with wave damping coefficient (β) decreasing substantially in the post-storm period; (2) The weakened wave attenuation capacity of salt marshes was attributed to storm-induced vegetation damage in stem density, with recovery taking more than a few weeks; (3) Under the scenario of storm sequences, cumulative damage in stem density caused significant decline in wave attenuation capacity of salt marshes, with β gradually decreasing until the vegetation disappeared completely. This study reveals the vulnerability of salt marsh wave attenuation to storm disturbances, particularly under storm sequences, providing critical insights for coastal wetland management under increasing storm frequency.

The dynamics of coastal dunes are often affected by a combination of environmental conditions and human interventions. During storms and calm periods, marine and aeolian flows cause sediment transport that results in sedimentation and erosion. Dunes typically grow due to aeolian transport of sediment and erode due to marine forcing by waves and currents. Shoreface nourishments can increase the sediment budget in the coastal profile which may influence the marine and aeolian sediment transport. To what extent shoreface nourishments influence sediment transport and the connected dynamics of coastal dunes is yet unknown. In this paper we investigate coastal dune volume change from yearly profile measurements and relate this to a shoreface nourishment program. The measured dune growth along the entire Dutch sandy coastline is characterized at most places by a break in slope where the dune growth increases after a specific ‘breakpoint’ year. The derived breakpoint years generally correspond with the start of the nourishment program with a delay of several years (1–5 years). These results provide a starting point for further elaborating on the relationship between the dynamics of coastal dunes and nourishments. The next step would be to delineate the potential influence of nourishments on either the growth of dunes through aeolian sediment transport or the erosion due to marine events.

During severe storms, waves can overtop dikes, leading to erosion of the crest and landward slope, which may ultimately result in breaching. To accurately model this erosion, the overtopping flow needs to be described in a time-dependent manner for each individual wave overtopping event. The peak flow velocity (upeak) and peak flow thickness (hpeak) are critical boundary conditions in this context. Previous studies have shown that these flow characteristics are related to the overtopping volume, yet often propose deterministic models that overlook the variability and interdependency between these characteristics.

The goal of this study is to address these gaps by explicitly quantifying the variation and interdependence of upeak and hpeak, using data from small-scale FlowDike experiments. We propose generalized distributions to describe the variation in these flow characteristics, with upeak varying by 13% to 23%, depending on the waterside slope angle, and hpeak varying by approximately 20%. Furthermore, the interdependency between upeak and hpeak is modeled using a Student-t copula (ν=9.361, ρ=−0.497), revealing a moderate negative correlation. This suggests that overtopping events with a high upeak are less likely to have a large hpeak, and vice versa.

The findings of this study can be directly applied to improve models that describe the loading caused by overtopping waves and the resulting erosion. By incorporating the variation and interdependence of upeak and hpeak, these models can provide a more detailed representation of the peak flow characteristics of overtopping waves. Furthermore, these insights can be applied to the design of wave overtopping simulators, enabling the simulation of more realistic overtopping flows by incorporating more of their natural variation.

During extreme high-water events in river systems, the load on a levee section may exceed its resistance, initiating the breaching process which eventually leads to levee failure. The success of an emergency measure to intervene in the initial phases of levee failure is mainly dependent on its timely application. Quick action is required to prepare and deploy an emergency measure before damages to the levee section have become irreparable. In this study, we investigate the key parameters for successful application of an emergency measure, focusing on the BresDefender case study. The BresDefender is a floating pontoon used by the Dutch military, which is intended to avoid or postpone levee failure. A model has been developed taking the operational steps in the implementation of the emergency measure during a high water and the (uncertainty) in the duration of these processes into account. The model is used to quantify the probability of successful operation to prevent levee failure due to overflow or slope instability. The probability of successful application of the BresDefender has been simulated for river flood situations in the Netherlands. For the river Rhine, where the examined cases were prone to slope instabilities, the probability of arriving in time was found to be 70%. But for the Meuse case, where the examined cases were prone to overflow, the probability of arriving in time was found to be only 0%. The critical steps in the process after the occurrence of damage to the levee are damage detection, the decision to repair the damage, the transport of the emergency measure, and the placement of the measure. By incorporating emergency measures in emergency preparedness procedures, the time required for the critical steps will be decreased and the probability of successful application of the emergency measure, i.e., its contribution to flood risk reduction, will be enhanced.

This research investigates how salt marshes contribute to both wave energy dissipation and spectral period transformation, advancing their role as a nature-based solution for coastal protection. Using laboratory simulations with a scaled barren foreshore, salt marsh and dike model, we examine the interactions between vegetation, water depth, and wave properties under varied conditions, including storm scenarios with irregular waves. Results indicate a case specific threshold at which the salt marsh model attenuates energy optimally, as for very shallow water depths wave energy is predominantly dissipated by the barren foreshore. The spectral wave period T m − 1 , 0 increases when waves propagate from deep to shallow water depths, as a result of wave breaking and generation of infragravity waves. The presence of salt marsh vegetation further enhances this effect by preferentially damping high frequency components. This highlights that an increase in T m − 1 , 0 in vegetated environments may not always correspond to an increased hydrodynamic load on the dike.

The passage of ships in confined waterways creates a stern wave that can overflow bank protection structures such as groins. This overflow, due to the long-period primary ship-induced waves, can be high in velocity, especially at the lee-side slope of groins, potentially causing significant damage to the structure. This study derives an equation to express overflow velocities, intended as a design tool for groins exposed to these types of waves. A detailed experimental investigation was performed on four physical models of groins with different slopes and stone sizes in the armor layer under the influence of different hydraulic heads. Particle image velocimetry (PIV) was used to measure the flow velocities at the crest and lee sides of the structure. All PIV measurements were performed thrice under free-flow conditions with no initial water level at the lee side of the structure. The depth- and time-averaged flow velocities (Uavg) were extracted from four positions along the lee-side slope and accelerated from 0.7 to 2.2 m/s. A dimensionless equation of the overflow velocities was obtained as a function of the hydraulic head (h), slope (θ), freeboard (Rc), and nominal rock diameter (dn50).

Rock groins in the Elbe Estuary are constructed to maintain proper water levels for navigation and for embankment erosion protection. At certain localities, significant damages to rock groins have been observed due to the primary ship-generated waves. Primary waves are generated along the ship's hull and then propagate toward the river banks and groin fields, appearing in the interaction with the structures as a turbulent overflow phenomenon. Eventually, this overflowing may cause damages mainly to the crest and leeward side of the groins. Since this overflowing is the most pronounced with large primary waves at certain water levels, the estimation of the probabilities of extreme primary waves is a key element for a safe and reliable design of groins. For this goal, nonparametric Bayesian networks (NPBNs) are used here to infer the probability distribution function of the extreme primary wave heights at the tip of a groin in the Elbe Estuary. Results demonstrate the suitability of the NPBN in their prediction. The model framework allows the designer to predict the probabilities of primary ship-generated waves at groins when the information of ship dimensions, nautical parameters, and waterway geometry is available. These probabilities can later be used for design purposes for current and future conditions.

Mangrove forests are vital for flood reduction, yet their failure mechanisms during storms are poorly known, hampering their integration into engineered coastal protection. In this paper, we aimed to unravel the relationship between the resistance of mangrove trees to overturning and root distribution and the properties of the soil, while avoiding damage to natural mangrove forests. We therefore (i) tested the stability of 3D-printed tree mimics that imitate typical shallow mangrove root systems, mimicking both damaged and intact root systems, in sediments representing the soil properties of contrasting mangrove sites, and subsequently (ii) tested if the existing stability models for terrestrial trees are applicable for mangrove tree species, which have unique shallow root systems to survive waterlogged soils. Root systems of different complexities were modeled after Avicennia alba, Avicennia germinans, and Rhizophora stylosa, and printed at a 1:100 scale using material densities matching those of natural tree roots, to ensure the geometric scaling of overturning moments. The mimic stability increased with the soil shear strength and root plate surface area. The optimal root configuration for mimic stability depended on the sediment properties: spreading root systems performed better in softer sediments, while concentrating root biomass near the trunk improved stability in stronger sediments. An adapted terrestrial tree resistance model reproduced our measurements well, suggesting that such models could be adapted to predict the stability of shallow-rooted mangroves living in waterlogged soils. Field tree-pulling experiments are needed to further confirm our conclusions with real-world data, examine complicating factors like root intertwining, and consider mangrove tree properties like aerial roots. Overall, this work establishes a foundation for incorporating mangrove storm damage into hybrid coastal protection systems.

Mangroves have been used for coastal protection in Bangladesh since the 1960s, but their integration with embankment designs has not been fully explored. This paper investigates the effect of existing mangroves on required embankment performance, with a focus on the wave-damping effect of mangroves. Existing mangroves reduce the required thickness of embankment revetment by up to 16–30% in the west, 47–82% in the central region, and 53–77% in the east. Notable mangrove sites include the belt south of polder 45 (Amtali), with an average width of 1.77 km, and the Kukri-Mukri polder, with an average width of 1.82 km. These mangroves reduce the need for thick slope protection, allowing the replacement of concrete revetments with softer materials, such as clay or grass, combined with mangrove foreshore. Additional large mangrove belts are found in Sandwip and Mirersarai. By replacing or reducing revetment requirements, mangrove forests can minimize carbon emissions from construction while providing carbon sequestration and other ecosystem services. This study can inform future sustainable investments in coastal protection systems by identifying areas where mangroves offer the greatest wave-damping benefits, which could be focus of follow-up feasibility studies.

Saltmarshes are a promising nature-based alternative for conventional flood protection. However, saltmarshes can erode under storm conditions, whereby the seaward edge of the saltmarsh often forms a vertical cliff. Despite its importance, the effect of storm conditions on erosion at the saltmarsh cliff remains understudied, especially when waves traverse over a cliff. This research investigates the complex flow patterns around a saltmarsh cliff non-intrusively using Particle Image Velocimetry (PIV) conducted through a series of scaled monochromatic wave flume experiments. We adopted realistic foreshore configurations (e.g. cliff heights) and hydraulic loading conditions from the Dutch Wadden Sea. Results show two local near-bed velocity maxima on top of the saltmarsh, created during different wave phases by water depth contraction, wave transmission and interaction between flow and vortices that are shed from the cliff. Under the wave crest, high onshore-directed near-bed velocities were measured at approximately 2.5–4 times the cliff height onshore from the cliff. Under the wave trough, high offshore-directed velocities were found at the marsh edge. Both onshore- and offshore-directed velocities increase with increasing cliff height, larger wave height or lower water depth. Vegetation on top of the marsh reduces both the incoming and outgoing velocities in front of the cliff. Increasing the cliff height resulted in a greater reduction in velocities by the vegetation. These results demonstrate how local near-bed velocity maxima and location are influenced by the presence of a cliff and the interaction with vegetation on top of the saltmarsh. This research highlights the vulnerability of the cliff even during inundation of the cliff and will help to implement saltmarshes as nature-based solutions for flood defence.

Marine pollution is a major global environmental problem. The transport and dispersion of marine pollution is driven by a wide range of hydrodynamic processes, including wave-induced currents (e.g., Stokes drift) that are generated by free-surface and internal gravity waves in density-stratified fluids. While the (Lagrangian-mean) Stokes drift is known to fundamentally change transport patterns, wave-induced Eulerian-mean currents, such as those generated in the presence of the Coriolis force due to the Earth's rotation, are generally less well understood. To address this, the Delta Transport Processes Laboratory (DTP-Lab), a multi-purpose lab with novel facilities and state-of-the-art equipment, is being constructed in the Hydraulic Engineering Laboratory at TU Delft. The DTP-Lab combines multiple components: a 4.40-m diameter turntable, which can support a (removable) 5-m long flume; a 12.7-m long stainless steel flume; a piston-type, wet-back, force-controlled surface wave generator; a pumping system to create any type of density stratification; and a 3D Particle Tracking Velocimetry (PTV) system. The design and construction of these components along with technical validation and performance tests are presented in this technical note. A scaling analysis demonstrates the suitability of the laboratory to investigate wave-induced current under rotation. The DTP-Lab will pioneer the combined experimental study of surface waves, density stratification and Coriolis forces. The DTP-Lab is presented here with the objective of giving practical information to future users and to describe its novelty and range of applications.

This study treats a detached homogenous low-crested structure (HLCS) made of Cubipod concrete elements placed seaward of a vertical wall (forming a basin in between) to reduce overtopping. Assessing the complex hydrodynamics and effects of changing the geometry of such a system in relation to overtopping reduction is challenging. The numerical model OpenFOAM was applied to this end. Forchheimer coefficients for wave transmission and the flow through the HLCS were calibrated and validated using existing physical modeling data (α = 500 and β = 1.0, with varying porosity based on the Cubipod shape), while the effect of the basin and vertical seawall was determined fully numerically. The crest freeboard (R_c), crest width (B), and basin length (L_B) were the main geometrical parameters that influenced the performance of the HLCS in reducing overtopping. An exponential decay was observed in the overtopping discharge when the values of these geometrical parameters increased. As L_B increased, this decay was primarily due to the dissipation of the broken-wave bores. The largest gradient in the predicted overtopping discharge was noted at R_c/H_s,i ≈ 0, B/H_s,i ≈ 4.5, and L_B/L_p ≈ 1.2, where H_s,i is the incident significant wave height and L_p is the peak wavelength in the basin.

The severity of damages to riverine structures, such as groynes or revetments, across German estuaries has increased in the past years due to the increase in the ship-induced loads. However, few studies can be found in the literature focused on the damage of rock slopes under ship wave attack. In this study, the field data of a rock-armored groyne (lateral slope 1/4 and rocks with nominal diameter Dn50»12.6cm and high density ρs=3.7t/m3) tracked for a year by Melling et al. (2020) is analyzed; the field campaign began after the structure was rebuilt and finished when the structure already presented severe damage. During this field campaign, the incident ship-induced primary waves and water levels were recorded, and laser scans of the groyne armor were taken. Using those field laser scans, damage curves along the life of the structure were derived. Also, ship data from the AIS was retrieved. Then, each increment of the damage (increment of the dimensionless eroded area, ΔSe) was related to a ship-wave event and a passing ship. The most significant variable to describe ΔSe was found to be the primary wave height Hp, while the best explanatory variables for Hp were the partial blockage factor, the ship length and width and the relative velocity of the ship. The shape of the dependence between these variables is also analyzed by pairs using copula space. A clear tail dependence is observed between several pairs. For instance, the pair ΔSe and Hp presents upper tail dependence, meaning that the high values of ΔSe and Hp are more correlated than the smaller ones. This implies that models more complex than Gaussian copula, which is commonly used in Coastal Engineering applications, might be needed to model the probabilistic dependence between the variables.

Vegetation in front of levees, dikes and seawalls can decrease wave energy and therefore contribute to the safety against flooding. However, wave damping predictions by vegetation are still inaccurate due to measurement and modelling uncertainties. Many studies focused on finding reliable predictive tools using scaled flume tests with vegetation mimics. Scaling down vegetation can however lead to discrepancies with realistic scales, known as scale errors. In this work scaled tests were conducted on 3D-printed elastic replicas of willow trees and compared with real-scale experiments. We identified differences in measured wave dissipation between the scaled hydraulic model (1:10) and its real-scale prototype with 5m high live willow trees under storm conditions (1:1). The maximum measured wave damping (30%) was roughly 1.5 times higher than the real-scale experiments (20%). Following the same trend of the real-scale experiments, this amount of wave height damping declined for larger water levels. Largest effects are argued to be due to increased viscous damping (smaller branch Reynolds numbers) and non-exact flexibility scaling. These significant deviations illustrate that full-scale experiments, although expensive, are still needed to validate the results of scaled experiments for woody vegetation. Alternatively, accounting for these discrepancies can make scaled experiments more reliable and expensive real-scale experiments less needed for wave damping studies on woody vegetation.

Riparian forests in front of dikes can dampen incoming waves and thereby contribute to flood safety. In real-scale flume experiments with live pollard willow trees (forming a 40-m-long forest), it was observed that during storm conditions, a maximum reduction of 20 % in incoming wave height could be achieved (van Wesenbeeck, et al., 2022). Notably, this amount of wave damping occurred at a water depth of 3 meters, aligning with the section of the trees with the maximum frontal-surface area. For a larger water depth, measured wave damping however declined. This is potentially partly caused by the natural tapering form of the trees. Typically, trees are characterized by smaller branch diameters and more flexible branches higher up in the canopy (McMahon & Kronauer, 1976); and flexible vegetation mimics are known to dampen less than rigid mimics due to motion (Van Veelen, T, Reeve, & Karunarathna, 2020). Hence, both the frontal-surface area and branch rigidity decrease along the height of the willow trees, potentially leading to less wave damping by the forest when subject to large waves at higher water levels.

Physical model tests have been performed to study static stability of rock-armoured mild slopes. Current stability design formulae for steeper rock-armoured slopes focus on plunging and surging waves. Slopes of 1:6 and milder usually have more spilling breakers which decreases the load. Also, on mild slopes displaced rocks more often remain present in the wave attack zone, which increases the strength. These aspects lead to an overdesigned structure when existing formulae for steep rock-armoured slopes are used. The present wave flume tests were used to understand the processes and develop a design formula for rock-armoured mild slopes with an impermeable core. These tests were performed for statically stable rock-armoured slopes of 1:6 to 1:10. The tests confirmed that not all existing damage parameters are able to accurately describe the static stability on milder slopes. For mild slopes it is more accurate to describe the damage based on the eroded depth rather than on the eroded area or number of moved stones. In this study, a design formula and guidelines are provided for practicing engineers that design or evaluate the stability of mild rock-armoured slopes.

Hard stabilization methods have traditionally been employed to mitigate coastal erosion. Concrete armour is widely used due to its high level of dependence, robustness, ease of production and cost effectiveness (Cooke et al., 2020; Pikey and Cooper, 2012). It is inevitable that coastline ‘armouring’ will continue to rise because of the growing human population and urbanization, desire for and value of coastal property, opposed to predicted climate change (Chapman and Underwood, 2011). The environmental impact of such 'armouring' on coastal systems can be detrimental, resulting in a degradation or destruction of habitats and the loss of ecologically trivial species (Gittman et al., 2015).

The CoastalockTM, a single-layer armour unit, aims to blend coastal protection with marine habitat creation. This armour unit is designed to mimic inter- and sub-tidal habitats, with chemical composition of substrate and micro and macro features that provide niches for various species. The key feature of CoastalockTM is the cavity that is integrated into the design, that caters to diverse marine life needs depending on its orientation (ECOncrete Tech Ltd., 2019). CoastalockTM's hydraulic performance is under research. Preliminary tests conducted in the Hydraulic Engineering Laboratory (HEL) of the Technical University of Delft (TUD) on a 2V:3H impermeable slope in deep water conditions highlighted that with tight placement of the units significant pressure gradients across the top layer led to damage. The introduction of spacings between units for enhanced permeability improved stability significantly (Gutiérrez et al., 2023). A redesign of the unit was proposed incorporating protrusions to enforce the spacings between the blocks (Molenkamp, 2022).

This research focuses on evaluating the influence of a porous core on the hydraulic performance of a CoastalockTM armour layer, specifically assessing its stability, overtopping, and reflection on a 2V:3H breakwater slope in deep water conditions—from the toe to just below the crest. A pivotal aspect of this research is the investigation of the impact of protrusions on the hydraulic performance. Furthermore, the study explores the influence of different toe configurations, aiming to comprehend the vulnerability of the armour layer to sliding. Toe scour falls outside the scope of this study.