Dune erosion
Definition of Dune erosion:
Sand loss from a dune under wave attack, mainly by notching, avalanching and slumping processes.
This is the common definition for Dune erosion, other definitions can be discussed in the article

This article deals with dune erosion by storms and provides some simple rules from which retreat of the dune front can be estimated. The focus is on the 'collision regime', where waves hit the dune without exceeding the dune crest (the socalled 'overwash regime').
Contents
Introduction
Coastal dunes have developed naturally along many shorelines worldwide during the last millennia when sealevel rise slowed down. Wave action and onshore winds are the main agents for coastal dune development and sufficient sand supply to the coast is a primary condition ^{[1]}, see also the article Dune development. The coastal dune belt in many cases protects lowlying hinterland from flooding by the sea. Dune erosion therefore can be a serious threat.
A coastal dune can suffer large losses when attacked by storm waves. The front dune can be taken away over several tens of meters, leaving a steep dune scarp, see Fig. 1. For the Dutch coast it has been estimated that the dune foot can recede as much as 80100 m under exceptional circumstances (extreme storms of very long duration, which may occur with a yearly probability of 1/100,000).
The actual dune loss depends on many factors:
 beach geometry: beach slope and width;
 grain size;
 dune geometry: dune face slope, dune crest elevation;
 nearshore morphology: nearshore slope, ripchannel embayments, nearshore sandbars, subtidal morphology;
 beach and dune vegetation; presence of foredunes;
 presence of hard engineered structures;
 alongshore variability in wave energy and dune height.
The dune toe elevation (beach width[math]\times[/math]slope) appears to be a major factor; observations show that vulnerability to erosion decreases with increasing dune toe elevation^{[2]}. However, there is also strong evidence that dunes are less impacted by storms when situated on dissipative flat beaches than on reflective steep beaches^{[3]}.
Dune loss also depends on storm characteristics: maximum sustained wind speed and duration; maximum water level and simultaneous wave height, period and incidence angle.
If at least part of the dune belt survives the storm, it can recover by natural processes – similar to the processes that created the original dune. See section Poststorm dune recovery.
As an example, Fig. 2 shows a model result of the relationship between dune retreat [math]RD[/math] (expressed as the distance between the initial dune foot location and the dune foot location after the storm surge) and exceedance frequency (the probability that in a particular year a storm occurs that produces greater erosion) for a location along the Dutch coast. According to this figure there is a 1/100,000 probability per year that dune erosion will exceed 85 m.
To establish such a plot, a proper insight is needed in the possibly occurring extreme storm conditions in the (very) small probability range. However, such insight is subject to great uncertainty, as reliable measurements of extreme water levels and wave conditions are available only over relatively short periods of the order of 100 years. With such a restricted data set, it is hard to make estimates of storm conditions with a yearly exceedance probability less than about 1/1000.
Impact of dune erosion
Buildings on the front dune situated close to the dune foot are at risk under severe storms, see Fig. 3. A sound estimate of potential dune erosion is required when issuing building permits. In cases where the hinterland is situated below sea level, the dune belt serves as sea defence. Breach of the dune belt may have catastrophic consequences. Dune and beach monitoring and dune management are of crucial importance. In some cases, where the dune belt consists of a single dune row, dune reinforcement or shore nourishment may be needed. To ensure safety, several methods have been developed for estimating dune loss during exceptional storm conditions.
Brief explanation of dune erosion processes
The initial crossshore beach profile, which might be considered to be in a more or less dynamic equilibrium condition with normally occurring hydrodynamic conditions, will be reshaped during a severe storm surge. The much higher water levels and much higher wave heights and peak periods call for a quite different shape of an equilibrium profile than the shape of the initial profile. Offshore directed sediment transports will occur, especially with sand eroded from the dunes.
Wave attack at the dune foot steepens the dune profile which may collapse by avalanching (Fig.4). Wave attack may also create a notch at the dune foot leading to mass failure: collapse of a dune slab, initiated by tensile cracking at the top surface of the dune followed by shear failure along an internal failure plane or overturning due to the weight of the overhang ^{[4]}, see Fig. 5.
The sand eroded from the dune is removed further down by the backwash and the seaward undercurrent ('undertow') under the breaking waves. Most of the sand is deposited on the lower beach or upper shoreface. In cases where strong rip currents are present or a strong ebb current from nearby inlets, some of the sand can even be transported far offshore^{[3]}. The slope of the crossshore profile gradually decreases, and consequently the rate of dune erosion will decrease with time during the storm surge. Usually, the limited storm surge duration does not allow for the development of a new equilibrium profile. The shape of the crossshore profile after the storm surge, which is a transient state between the initial profile and the storm equilibrium profile, is often called 'storm erosion profile'.
Quantification of dune erosion
Although 3D effects are generally important in the dune erosion process, often a 1D or 2DV approach is adopted, ignoring alongshore variability. In this case the dune erosion process is considered as a typically offshore directed crossshore sediment transport problem. Sand from the dunes is transported to deeper water and settles there.
A first approach consists of assuming a closed sediment balance in crossshore direction. The same volume of sand which is eroded from the dunes and the very upper part of a crossshore profile is accumulated lower in the crossshore profile. [Because of differences in porosity of the eroded dune material (often loosely packed) and the settled material (often a bit more densely packed), the volume balance is not always strictly closed.]
During severe storm surges, with a great water level increase (a water level increase of a few metres is possible along some coasts) huge volumes of sand from the dunes are transported in offshore direction. And because dune erosion is a rather shortlasting process, some computation methods take only offshoredirected transports into account.
Erosion profile method
A quick first order estimate of storm dune erosion can be obtained by assuming that the storm erosion profile is close to the storm equilibrium profile. This assumption is unrealistic for storms of short duration and strongly overestimates storm dune erosion in this case. However, for storms of very long duration, which produce great dune losses with low exceedance probability, the assumption of a poststorm equilibrium profile is not unreasonable. The method is still used in the Netherlands, besides other more advanced methods, to provide an upper bound of possible storm dune erosion. It is based on empirical poststorm equilibrium profiles established by laboratory experiments and validated by field data of the Dutch coast ^{[5]}. The poststorm equilibrium profile was established for different values of the significant wave height [math]H[/math], peak wave period [math]T[/math] and mean fall velocity of dune sand [math]w[/math]. Major assumption are:
 the storm duration is sufficient for establishment of a poststorm equilibrium profile [math]y(x)[/math];
 the poststorm equilibrium profile does not (strongly) depend on the initial coastal profile [math]y_0(x)[/math];
 eroded dune sediment is deposited within a zone delimited by a storm closure depth [math]y_{max}[/math];
 the eroded dune scarp has a slope 1:1 and the slope at the toe of the sand deposit is 1:12 (these are less crucial assumptions).
With these assumptions the eroded dune volume is given by the beach volume between [math]y(x)[/math] and [math]y_0(x)[/math], such that the total sand volume landward of the closure depth [math]y_{max}[/math] is preserved. The procedure is explained in Fig. 6, where the empirical parameterized functions for the poststorm equilibrium profile [math]y(x)[/math] and the storm closure depth [math]y_{max}[/math] are also indicated.
The slope of the poststorm profile decreases (i.e. increasing dune erosion) as the wave period [math]T[/math] increases and the fall velocity decreases. A serious drawback of this straightforward method is that hardly any physics is involved. The development with time of the storm profile is unknown. Effects of varying water levels and varying wave characteristics during the storm surge cannot be accounted for.
Analytic dune erosion model
A simple analytic model was developed by Larson et al. (2004^{[6]}, 2016^{[7]}) for estimating the time evolution of dune erosion. This model is based on the assumption that the eroded dune volume [math]\Delta V[/math] in a time interval [math]\Delta t[/math] is proportional to the force [math]F[/math] exerted by swash waves hitting the dune foot. By expressing the force [math]F[/math] as a function of the dune foot height [math]z_D[/math], the maximum water level [math]R[/math] including wave runup and the wave period [math]T[/math], the following equation was obtained for the loss rate [math]dV/dt[/math] of dune volume (see the appendix and Fig. B1):
[math]\Large\frac{dV}{dt}\normalsize=4C_s\Large\frac{(Rz_D)^2}{T}\normalsize , \qquad (1)[/math]
where the sitedependent constant [math]C_s[/math] has to be determined empirically. According to this equation, the rate of dune erosion is an increasing function of the wave height, the wave period and the beach slope (see appendix). This is consistent with the dune erosion estimate following from the storm profile of Fig. 6. Equation (1) also shows that dune erosion strongly depends on the height [math]z_D[/math] of the dune foot (i.e. on the beach steepness and width). This is consistent with observations at NarrabeenCollaroy Beach (East Australia, near Sidney), which showed that alongshore variability in dune erosion was strongly correlated with alongshore variability in [math]z_D[/math]^{[2]}. Assuming that the beach slope [math]\beta[/math] remains constant during dune erosion (i.e., dune retreat [math]\Delta x[/math] equal to [math]\Delta z_D / \tan \beta[/math]), Larson et al. (2004^{[6]}) found reasonable agreement with observed dune erosion rates when using the parameterization [math]C_s=1.34 \, 10^{3} \exp(3.19 \, 10^{4} H/d_{50})[/math], where [math]d_{50}[/math] is the median sediment grainsize and [math]H[/math] the rootmeansquare wave height during the storm.
Contrary to the erosion profile method, the analytic model only considers the maximum water level [math]R[/math] without distinction between the wave runup and storm level components. Observations suggest that the influence of the storm level component is more important for dune erosion than the wave runup component^{[8]}.
DUROSTA
Based on theoretical work and laboratory experiments, Steetzel (1993) ^{[9]} has developed the 2dimensional DUROSTA computation model, in which actual sediment transports are calculated based on waveintegrated velocity and sand concentration profiles. Although in the mathematical description of the crossshore sediment transport the intrawave component is neglected, the results of the model compared well with large scale model tests in the Delta Flume of Delft Hydraulics. During storm surge conditions, breaking waves in the surf zone cause a strong return flow in the lower part of the water column and high sediment concentrations throughout the water column. The crossshore profile is updated at each time step, based on gradients in the computed sediment fluxes. The DUROSTA model therefore provides estimates of the time evolution of the dune and beach profiles during a storm surge.
XBeach
XBeach ^{[10]} is an opensource processbased numerical model that was originally developed to simulate hydrodynamic and morphodynamic processes and impacts on sandy coasts with a domain size of kilometers and on the time scale of storms. Since then, the model has been applied to other types of coasts and purposes. The model includes the hydrodynamic processes of short wave transformation (refraction, shoaling and breaking), long wave transformation (generation, propagation and dissipation of infragravity waves), wave setup and unsteady currents, as well as overwash and inundation. The morphodynamic processes include bed load and suspended sediment transport, bed update and breaching. Effects of vegetation and hard structures have been included. The model has been validated with a series of analytical, laboratory and field test cases using a standard set of parameter settings. The avalanching of sandy material from the dune face to the foreshore during storm conditions is taken into account when updating the bed levels. This is modeled through introduction of a critical bed slope (default critical slope of 1 for dry zones and 0.3 for wet zones). When the critical slope is exceeded, material is exchanged between adjacent cells to the amount needed to bring the slope back to the critical slope.
Comparison with observed dune erosion events at the Dutch coast show that XBeach estimates the magnitude and pattern of alongshore variations in erosion volume reasonably well. The 2014version of XBeach overpredicted the erosion volume in the region where a dune scarp developed and underestimated the erosion volume where the whole dune face collapsed in a series of slumps ^{[11]}. XBeach simulations further illustrated that the observed alongshore variation in dune erosion was steered primarily by the prestorm dune topography i.e., the presence of embryonic dune fields and the steepness of the dune face. The importance of alongshore variability in intertidal beach topography was found to be secondary, but not negligible during the initial stage of the storm, when the surge level was still low. An application of the 2019 version of XBeach to a breach experiment at a largescale sand dune in Germany showed that the poststorm beach profile was well simulated, but dune erosion was overestimated with respect to temporal profile evolution. It was concluded that the XBeach model is capable of accurately predicting the evolution of storminduced erosion during a sufficiently long and overwash/inundationdominant surge event when current highresolution input data is available^{[12]}.
Poststorm dune recovery
After being eroded under storm conditions, the dune front can recover naturally without human intervention. This has been observed for many sandy coasts with either a gently sloping or a steep sloping shoreface ^{[13]}^{[14]}^{[15]}. Strong beach accretion often already occurs during the phase of waning storm intensity. Afterstorm dune recovery starts with the development of embryo dunes on the backshore stimulated by vegetation growth^{[16]}. An example of natural dune recovery is shown in Fig. 7 for a retreating coast (Walcheren, the Netherlands). The extreme storm surge of 1953 caused a 30 m retreat of the dune foot. Afterwards the dune foot advanced slightly over a period of several years, until it reached about 15 years later the position corresponding to the retreating trend prior to the storm surge. The ongoing coastal retreat was finally halted in 1990, when the Dutch government decided to stop coastal retreat along the entire coastline by coastal sand nourishments. Recovery of the shoreline precedes recovery of the dune foot and is much faster. Fig. 8 shows afterstorm recovery of the shoreline position at NarrabeenCollaroy Beach in East Australia, near Sidney. The recovery period after each storminduced shoreline retreat is only a few months at this location. Shoreline recovery at this location is significantly correlated with onshore migration of nearshore sandbars^{[15]}^{[17]}. It further appears that poststorm recovery of wide gently sloping beaches is promoted by infragravity swash^{[8]}.
Impact of sealevel rise on dune erosion
The XBeach model predicts a linear relationship between dune erosion volume and sea level rise ^{[19]}. This is primarily due to the higher water level in front of the dune and not to changes in the significant shortwave or infragravity wave height. Changes in the offshore angle of wave incidence also affect the dune erosion volume. A 30° shift from shorenormal influences the dune erosion volume to the same extent as a 0.4 m sealevel rise. This increase in dune erosion volume is related to strong alongshore currents, generated as a result of the obliquity of the waves, that enhance stirring and hence offshore transport.
The model simulations further show that the effectiveness of coastal sand nourishment to mitigate the impact of sealevel rise, strongly depends on the location in the profile where this sand is added. However, the reduction in dune erosion volume relative to the total added volume remains low (<0.3) for all mitigation options. This suggests that directly increasing the volume of sand in the dunes may be more efficient than beach or shoreface nourishment.
Appendix: Analytical dune erosion model
The analytic model of Larson et al. (2004^{[6]}) is based on the following assumptions and simplifications:
 The eroded dune volume [math]\Delta V[/math] in a time interval [math]\Delta t[/math] is proportional to the force [math]F[/math] exerted by swash waves hitting the dune foot, [math]\Delta V \propto F \Delta t .[/math]
 The force [math]F \Delta t [/math] is proportional to the number of swash bores [math]n[/math] hitting the dune foot ([math]n=\Delta t/T[/math]) and the force exerted by each individual swash bore, [math]f=m \Large\frac{du}{dt}\normalsize \propto \Large\frac{mu}{T}\normalsize [/math], where [math]T[/math] is the wave period, [math]m[/math] the mass of the swash bore per unit dune width and [math]u[/math] the velocity of the swash bore when hitting the dune foot.
 The bore height [math]h[/math] is related to the bore velocity by [math]u \propto \sqrt{gh}[/math], where [math]g[/math] is the gravitational acceleration. The bore mass is thus proportional to [math]m \propto huT \propto u^3 T[/math], yielding [math]F = n \,f / \Delta t \propto u^4 /T[/math].
 The bore velocity at the dune foot can be related to the swash runup [math]R[/math]. A ballistic swash excursion model (neglecting friction) gives an estimate for this relation, [math]R=z_D+gu^2/2[/math], where [math]z_D[/math] is the height of the dune foot compared to the beach level where waves collapse on the beach. Using this relation, the force [math] F [/math] can be expressed as a function of the wave runup [math] R [/math], [math]F \propto (R – z_D)^2 / T[/math].
 Substitution with [math]\Delta V / \Delta t \approx dV/dt [/math] gives Eq. (1), where [math]C_s[/math] is a proportionality constant, depending on characteristics of the site, such as sediment grainsize, wave incidence angle, shoreface profile, etc. The coefficient [math]C_s[/math] should be determined empirically; usual values are found in the range [math]10^{3}\;  \; 2.5\;10^{3}[/math].
 The wave runup [math]R[/math] depends on the beach slope [math]\beta[/math], the offshore wave height [math]H[/math] and wave period [math]T[/math]; an approximate empirical expression is [math]R \sim \tan \beta \sqrt{HL} [/math], where [math]L \approx g T^2/2 \pi[/math] is the offshore wavelength (see also Swash zone dynamics for more details).
Substitution in Eq. (1) gives [math]\Large\frac{dV}{dt}\normalsize = 4C_s \tan \beta (\sqrt{\Large\frac{gH}{2 \pi}\normalsize }  \Large\frac{l}{T}\normalsize) [/math], where [math]l[/math] is the beach width, which increases as the dune foot retreats.
This expression shows that the rate of dune erosion increases as the wave height [math]H[/math] increases, the wave period [math]T[/math] increases, the beach slope [math]\beta[/math] increases and the beach width [math]l[/math] decreases. Greater dune erosion for steeper beaches under the same wave conditions is consistent with observations^{[8]}.
The dune foot height [math]z_D[/math] (relative to the beach level of wave collapse [math]z_C[/math]) and the beach slope [math]\beta[/math] are time dependent. In order to solve Eq. (1) a prescription must be given how the eroded dune material is distributed over the beach – i.e. a model for the dependence of [math]z_D[/math] and [math]\beta[/math] on [math]dV/dt[/math]. Based on experiments in large wave flumes, Larson et al. (2004^{[6]})assumed that the beach slope [math]\beta[/math] remains constant during dune erosion. The dune retreat [math]\Delta x[/math] was therefore set equal to [math]\Delta z_D / \tan \beta[/math], assuming that the major part of the dune volume loss is removed to the lower beach (below the level of wave collapse).
Considering a dune height [math]D[/math] above the level of wave collapse [math]z_C[/math], the loss of dune volume is given by [math]\Delta V=(Dz_D) \Delta x = (D z_D) \Delta z_D / \tan \beta[/math]. Combining this result with Eq. (1) yields a differential equation for the dune foot height [math]z_D[/math]:
[math]\Large\frac{d}{dt}\normalsize z_D = \Large\frac{1}{\tau} \frac{(Rz_D)^2}{Dz_D}\normalsize , [/math] where [math]\tau=T/(4 C_s \tan \beta)[/math].
If it is further assumed that the level of wave collapse [math]z_C[/math] does not change during the storm, then the wave runup [math]R[/math] remains constant too (the slope [math]\beta[/math] was assumed constant). The differential equation can then be solved, with solution
[math]t = \tau (DR)\Large(\frac{1}{Rz_D(0)}  \frac{1}{Rz_D(t)})\normalsize +\tau \ln\Large\frac{rz_D(0)}{rz_D(t)}\normalsize [/math], from which [math]z_D(t)[/math] and [math]V(t)[/math] can be derived.
Related articles
 Natural causes of coastal erosion
 Dune stabilisation
 Light revetments builtin into artificial dunes
 Shoreface profile
 Bruun rule
 Dune development
 Swash zone dynamics
 Shore protection vegetation
 Dynamics, threats and management of dunes
 Risk and coastal zone policy: example from the Netherlands
References
 ↑ Anthony, E.J., MraniAlaoui. M. and Héquette, A. 2010. Shoreface sand supply and mid to late Holocene aeolian dune formation on the stormdominated macrotidal coast of the southern North Sea. Marine Geology 276: 100–104
 ↑ ^{2.0} ^{2.1} Splinter, K.D., Kearney, E.T. and Turner, I.L. 2018. Drivers of alongshore variable dune erosion during a storm event: Observations and modelling. Coastal Engineering 131: 31–41
 ↑ ^{3.0} ^{3.1} Masselink, G. and van Heteren, S. 2014. Response of wavedominated and mixedenergy barriers to storms. Marine Geology 352 (2014) 321–347
 ↑ Erikson, Li, H., Larson, M. and Hanson, H. 2007. Laboratory investigation of beach scarp and dune recession due to notching and subsequent failure. Marine Geology 245: 1–19.
 ↑ van Gent, M.R.A., van Thiel de Vries, J.S.M., Coeveld, E.M., de Vroeg, J.H. and van de Graaff, J. 2008. Largescale dune erosion tests to study the influence of wave periods Coastal Engineering 55: 1041–1051
 ↑ ^{6.0} ^{6.1} ^{6.2} ^{6.3} Larson, M., Erikson, L. and Hanson, H. 2004. An analytical model to predict dune erosion due to wave impact. Coastal Engineering 51: 675– 696
 ↑ Larson, M., Palalane, J., Fredriksson, F. and Hanson, H. 2016. Simulating crossshore material exchange at decadal scale. Theory and model component validation. Coastal Engineering 116: 57–66
 ↑ ^{8.0} ^{8.1} ^{8.2} Cohn, N., Ruggiero, P., GarcíaMedina, G., Anderson, D., Serafin, K.A. and Biel, R. 2019. Environmental and morphologic controls on waveinduced dune response. Geomorphology 329: 108–128 Cite error: Invalid
<ref>
tag; name "C19" defined multiple times with different content Cite error: Invalid<ref>
tag; name "C19" defined multiple times with different content  ↑ Steetzel, H.J. 1993. Crossshore transport during storm surges. Thesis Tech. Univ. Delft, Delft Hydraulics Communications 476.
 ↑ https://xbeach.readthedocs.io/en/latest/user_manual.html
 ↑ De Winter, R.C,. Gongriep, F. and Ruessink, B.G. 2014. Observations and modeling of alongshore variability in dune erosion at Egmond aan Zee, the Netherlands. Coastal Engineering 99: 167175.
 ↑ Schweiger, C., Kaehler, C., Koldrack, N. and Schuettrumpf, H. 2020. Spatial and temporal evaluation of storminduced erosion modelling based on a twodimensional field case including an artificial unvegetated research dune. Coastal Engineering 161, 103752
 ↑ Vousdoukas, M.I., Pedro, L., Almeida, M. and Ferreira, O. 2012. Beach erosion and recovery during consecutive storms at a steepsloping, mesotidal beach. Earth Surf. Process. Landf. 37: 583–593
 ↑ Scott. T., Masselink, G., O'Hare, T., Saulter, A., Poate, T,. Russell. P., Davidson, M. and Conley, D. 2016. The extreme 2013/2014 winter storms: Beach recovery along the southwest coast of England. Marine Geology 382: 224–241
 ↑ ^{15.0} ^{15.1} ^{15.2} Phillips, M.S., Harley, M.D., Turner, I.L., Splinter, K.D. and Cox , R.J. 2017. Shoreline recovery on wavedominated sandy coastlines: the role of sandbar morphodynamics and nearshore wave parameters. Marine Geology 385: 146–159
 ↑ Suanez, S., Cariolet, JM., Cancouët, R., Ardhuin, F. and Delacourt, C. 2012. Dune recovery after storm erosion on a highenergy beach: Vougot Beach, Brittany (France). Geomorphology 139140: 16–33
 ↑ Brooks, S.M., Spencer, T. and Christie, E.K. 2017. Storm impacts and shoreline recovery: Mechanisms and controls in the southern North Sea. Geomorphology 283: 48–60
 ↑ Jeuken, C., Ruessink, G. and Marchand, M. 2001. Ruimtelijke en temporele aspecten van de duinvoetdynamiek. Report WL/Delft Hydraulics Z2838 (in Dutch)
 ↑ De Winter, R,C. and Ruessink, B.G. 2017. Sensitivity analysis of climate change impacts on dune erosion: case study for the Dutch Holland coast. Climatic Change (2017) 141:685–701 DOI 10.1007/s1058401719223
Please note that others may also have edited the contents of this article.
