Difference between revisions of "Dune erosion"

From Coastal Wiki
Jump to: navigation, search
(Large scale safety problem)
(25 intermediate revisions by 4 users not shown)
Line 1: Line 1:
 +
 
{{Definition|title=Dune erosion
 
{{Definition|title=Dune erosion
|definition= [[Dune]] erosion involves that, during a severe [[storm surge]], sediments from the mainland and upper parts of the [[beach]] are eroded and settled at deeper water within a short time period; this is a typical cross-shore sediment transport process.}}
+
|definition= Sand loss from a dune under wave attack, mainly by notching, avalanching and slumping processes.}}
 +
 
 +
 
 +
This article deals with dune erosion by storms and provides some simple rules from which retreat of the dune front can be estimated.
 +
 +
 
 +
==Introduction==
 +
 
 +
[[Image:DuneScarpEgmond.jpg|400px|thumbnail|right|Fig. 1. Post-storm eroded dune, Egmond aan Zee.]]
 +
 
 +
Coastal dunes have developed naturally along many shorelines worldwide during the last millennia when sea-level 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 <ref>Anthony, E.J., Mrani-Alaoui. M. and Héquette, A. 2010. Shoreface sand supply and mid- to late Holocene aeolian dune formation on the storm-dominated macrotidal coast of the southern North Sea. Marine Geology 276: 100–104</ref>. The coastal dune belt in many cases protects low-lying 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 80-100 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 local conditions such as beach width, beach height and dune height. 
 +
 
 +
If at least part of the dune belt survives the storm, it can recover by natural processes – the same processes that created the original dune. This will happen only if the conditions under which the dune belt was formed – sand supply, wave and wind climate – have not significantly changed.
 +
 
 +
[[Image:imageDu002.JPG|350px|thumbnail|left|Fig. 2. Storm-induced dune foot retreat as a function of frequency of exceedance for the Dutch coast.]]
 +
 
 +
 
 +
 
 +
 
 +
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==
 +
 
 +
[[Image:BadhotelSchiermonnikoog_1924.jpg|350px|thumbnail|right|Fig. 3. Badhotel at Schiermonnikoog after the 1924 storm.]]
 +
 
 +
 
 +
 
 +
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 cross-shore 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 <ref> 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 (2007) 1–19.</ref>, see Fig. 5.
 +
 
 +
{| border="0"
 +
|-
 +
| valign="top"|
 +
[[File: DuneFailureMechanisms.jpg|450px|thumbnail|left|Fig. 4. Dune failure mechanisms.]]
 +
| valign="top"|
 +
[[File: DuneErosionEgmond.jpg|450px|thumbnail|left|Fig. 5. Dune scarp after collapse. The people standing on the dune scarp have not read this article.]]
 +
|}
 +
 
 +
 
 +
The sand deposited at the dune foot is removed by the seaward undercurrent ('undertow') under the breaking waves and deposited on the shoreface. The slope of the cross-shore profile gradually decreases, and consequently the rate of dune erosion will decrease with time during the storm surge. It is, however, not expected that a new equilibrium profile will develop because of the limited storm surge duration. The shape of the cross-shore 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 undoubtedly important in the dune erosion process, often a 2D approach is adopted. In this case the dune erosion process is considered as a typically offshore directed cross-shore 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 cross-shore direction. The same volume of sand which is eroded from the dunes and the very upper part of a cross-shore profile is accumulated lower in the cross-shore 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 sediment 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 short lasting process, some computation methods take only offshore directed 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 post-storm 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 post-storm equilibrium profiles established by laboratory experiments and validated by field data of the Dutch coast <ref>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. Large-scale dune erosion tests to study the influence of wave periods Coastal Engineering 55: 1041–1051</ref>. The post-storm 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 post-storm equilibrium profile <math>y(x)</math>;
 +
# the post-storm 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 parameterized functions for the post-storm equilibrium profile <math>y(x)</math> and the storm closure depth <math>y_{max}</math> are also indicated.
 +
 
 +
 
 +
 
 +
[[Image:DUROSmodel.jpg|700px|thumbnail|center|Fig. 6. Erosion profile method for estimating dune erosion as given by the beach volume between  <math>y(x)</math> and <math>y_0(x)</math>. Underlying assumptions are explained in the text. ]]
 +
 
 +
 
 +
 
 +
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, 2016)<ref>Larson, M., Erikson, L. and Hanson, H. 2004. An analytical model to predict dune erosion due to wave impact. Coastal Engineering 51: 675– 696</ref><ref>Larson, M., Palalane, J., Fredriksson, F. and Hanson, H. 2016. Simulating cross-shore material exchange at decadal scale. Theory and model component validation. Coastal Engineering 116: 57–66 </ref> 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 zone dynamics|swash waves]] hitting the dune foot. This force 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 du/dt \propto mu/T </math>, where <math>T</math> is the wave period, <math>m</math> the mass of the [[Swash zone dynamics|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> (<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=nf \propto u^4 \Delta t/T</math>. The bore velocity at the dune foot can be related to the swash [[Wave run-up|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. The estimate of the dune volume loss <math>dV/dt</math> per unit time obtained by the model is finally given by
 +
 
 +
<math>\Large\frac{dV}{dt}\normalsize=-4C_s\Large\frac{(R-z_D)^2}{T}\normalsize. \qquad (1)</math>   
  
The basis of this article is especially written for the Coastal Wiki by the main author referred to at the bottom of this page.
+
The [[Wave run-up|runup]] <math>R</math> and the coefficient <math>C_s</math> depend on local conditions. The 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 \propto \tan \beta \sqrt{HL} </math>, where <math>L \approx g T^2/2 \pi</math> is the offshore wavelength. 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 approximate linear dependence of the [[Wave run-up|runup]] <math>R</math> on <math>\beta</math> and <math>T</math> implies that dune erosion is an increasing function of the beach slope and the wave period. 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). This is corroborated by a field study at Narrabeen-Collaroy Beach (East Australia, near Sidney), which showed that alongshore variability in dune erosion was strongly correlated with alongshore variability in <math>z_D</math><ref>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</ref>.  
  
This article is under construction; not yet finished.
+
===DUROSTA===
  
 +
Based on theoretical work and laboratory experiments, Steetzel (1993) <ref>Steetzel, H.J. 1993. Cross-shore transport during storm surges. Thesis Tech. Univ. Delft, Delft Hydraulics Communications 476.</ref> has developed the 2-dimensional DUROSTA computation model, in which actual sediment transports are calculated based on wave-integrated velocity and sand concentration profiles. Although in the mathematical description of the cross-shore sediment transport the intra-wave 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 cross-shore 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===
  
==Dunes as sea defense==
+
XBeach <ref>https://xbeach.readthedocs.io/en/latest/user_manual.html</ref> is an open-source process-based 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 (infragravity wave) transformation (generation, propagation and dissipation), wave-induced 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.  
River and sea dikes are good examples of structures to protect low-lying areas from flooding. Also a dune area (dune row) serves that aim in some cases. E.g. the safety of large parts of The Netherlands, often with ground levels even below Mean Sea Level (MSL), relies on dikes, but also on dunes for their protection against flooding.
 
  
Visiting the beach and the coastal zone in e.g. The Netherlands under normal weather conditions would easily give the impression that the dunes are certainly strong enough to properly protect the hinterland. However, during a severe storm surge, with under design conditions water levels at sea which are approximately 5 - 6 m above MSL and together with the much more severe wave conditions than normal (cf. wave heights Hs ≈ 7 - 9 m and peak periods Tp ≈ 12 -18 s), the dunes will be eroded in a very short period of time.
 
  
Existing design rules in The Netherlands yield erosion rates of 80 - 100 m of the dunes during design storm conditions. (The rates of 80 - 100 m are given as an order of magnitude value only to facilitate the further discussion; the actual erosion rates under design conditions depend on the specific local conditions; e.g. shape of initial cross-shore profile and particle size of the dune material.) It must be realized that because of the specific Dutch conditions, the design conditions in The Netherlands are very strict.
+
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 2014-version 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 <ref>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: 167-175.</ref>.  XBeach simulations further illustrated that the observed alongshore variation in dune erosion was steered primarily by the pre-storm 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.
  
===Rather wide dune area is sea defense===
 
Often the dune areas in The Netherlands are wide enough to accommodate 80 - 100 m of dune erosion during a single severe storm surge. In some cases, however, the row of dunes is rather slender; a careful judgement has to be passed whether the dunes provide the required rate of protection. Is a break-through expected under design conditions?, and if yes: what reinforcement is necessary to fulfil the requirements?
 
  
So far in the discussion the safety problem of people living behind the dunes was raised as the main issue. In the Sections 'Large scale safety problem' and 'Small scale safety problem'of this article it will be shown that more issues related to dune erosion are relevant to a coastal zone manager. Aspects like the safety of single houses and hotels in the erosion zone are dealt with.
+
==Post-storm dune recovery==
  
Also topics like how to deal with structural erosion and global sea level rise are relevant topics for a coastal zone manager. They are briefly discussed in this article with the present Dutch policy and insights as starting points.
+
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 <ref>Vousdoukas, M.I., Pedro, L., Almeida, M. and Ferreira, O. 2012. Beach erosion and recovery during consecutive storms at a steep-sloping, meso-tidal beach. Earth Surf. Process. Landf. 37: 583–593</ref><ref>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</ref><ref name=Ph>Phillips, M.S., Harley, M.D., Turner, I.L., Splinter, K.D. and Cox , R.J. 2017. Shoreline recovery on wave-dominated sandy coastlines: the role of sandbar morphodynamics and nearshore wave parameters. Marine Geology 385: 146–159</ref>. After-storm recovery starts with the development of embryo dunes on the backshore and is stimulated by vegetation growth on these embryo dunes<ref name=Su>Suanez, S., Cariolet, J-M., Cancouët, R., Ardhuin, F. and Delacourt, C. 2012. Dune recovery after storm erosion on a high-energy beach: Vougot Beach, Brittany (France). Geomorphology 139-140: 16–33</ref>. 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 30m 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 after-storm recovery of the shoreline position at Narrabeen-Collaroy Beach in East Australia, near Sidney. The recovery period after each storm-induced shoreline retreat is only a few months at this location. Shoreline recovery at this location is significantly correlated with onshore migration of nearshore sand bars<ref name=Ph></ref><ref> 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</ref>.
  
From the discussions it will become clear that for many reasons a proper insight in the rates of dune erosion as a function of the boundary conditions is necessary. Section 'Quantification of rate of dune erosion' will deal with the dune erosion process and the methods to quantify the rates of erosion during a severe storm surge.
+
{| border="0"
 +
|+Overall caption (optional)
 +
|-
 +
| valign="top"|
 +
[[File:DuneFootRecovery1953.jpg|thumb|left|450px|Fig. 7. Variation of the dune foot position in a beach transect of the island Walcheren (Netherlands) over the period 1900-2010 (solid white line, based on annual surveys of Rijkswaterstaat <ref>Jeuken, C., Ruessink, G. and Marchand, M. 2001. Ruimtelijke en temporele aspecten van de duinvoetdynamiek. Report WL/Delft Hydraulics Z2838 (in Dutch)</ref>). The trend is indicated by the dashed line, showing ongoing retreat of the dune foot from 1900-1990. The large dune erosion by the exceptional 1953 storm surge has no influence on the long-term trend of dune foot retreat. Stabilization after 1990 is due to  implementation of the current 'hold-the-line' coastal policy by sand nourishments.]]
 +
| valign="top"|
 +
[[File:NarrabeenCollaroyBeachShorelineRecovery.jpg|thumb|left|450px| Figure 8. Variation of the shoreline position in a transect of Narrabeen-Collaroy Beach (East Australia, near Sidney) over the period October 2004-May 2005. Strong shoreline retreat occurs during the storm periods indicated in light blue. After-storm recovery occurs during the periods in between. Figure redrafted after <ref name=Ph></ref>.]]
 +
|}
  
===Brief explanation of dune erosion process===
 
Figure 1 shows schematically what happens with a (in this case: wide enough) dune during a severe storm surge.
 
[[Image:image004.JPG|350px|thumbnail|right|Figure 1 Dune erosion due to storm surge]]
 
The initial cross-shore profile, which might be considered to be in a more or less dynamic equilibrium condition with the normal occurring boundary conditions, will be reshaped during the severe storm surge. The  much higher water levels and the 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 sediments from the dunes. Sand is eroded from the dunes and is settled at the foreshore again. During these reshaping processes the slopes of the cross-shore profile gradually decrease, and consequently it can be understood that the rate of dune erosion will decrease with time during the storm surge. It is, however, not expected that a real equilibrium profile will develop during the storm surge. The time available during the storm is too short to achieve such a real equilibrium profile, but the developments are in the direction of achieving equilibrium. The shape of a cross-shore profile as encountered after the storm surge is often called 'erosion profile'.
 
  
Right after a severe storm surge, when the boundary conditions are normal again, the shape of the (erosion) profile does not fit with these normal conditions. Onshore directed sediment transports will occur; wind will blow sand from the beach to the dunes; the 'old' situation will gradually be restored. Dune erosion because of severe storm surges is thus a temporary and a reversible process.
+
==Impact of sea-level rise on dune erosion==
  
===Large scale safety problem===
+
The XBeach model predicts a linear relationship between dune erosion volume and sea level rise <ref>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/s10584-017-1922-3</ref>. This is primarily due to the higher water level in front of the dune and not to changes in the significant short-wave or [[Infragravity waves|infragravity]] wave height. Changes in the offshore angle of wave incidence also affect the dune erosion volume. A 30° shift from shore-normal influences the dune erosion volume to the same extent as a 0.4 m sea-level 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.  
If, like in some cases, the dunes are really slender and landward of the row of dunes the ground levels are even below MSL (like in some places in The Netherlands) a large scale safety problem might occur. A break-through of the dunes causes flooding and will cause loss of lives and results in a lot of damage in the densely populated low-lying areas of The Netherlands.
 
Given the dimensions of a row of dunes and the specific design conditions, the Coastal Zone Manager has to judge whether the dunes are 'safe' or not. 'Safe' is in this respect in fact a relative notion. 'Safe' means that the strength (width, height) of the dunes fulfils the requirements. Absolute safety does not exist in many cases; often a set of even more extreme boundary conditions is conceivable which will result in a break-through of a row of dunes which was judged just 'safe' enough. The chance that such a set of boundary conditions will occur, is then, however, apparently smaller than has been agreed for the design conditions.
 
For the judgement of the safety of a row of dunes, a proper computation model is necessary and (a set of) design conditions. Section 'Quantification of rate of dune erosion' deals with these computation models.
 
  
For the further discussion in the present section it is good to realize that because of the many (stochastic) parameters which ultimately determine the rate of dune erosion, a probabilistic approach seems to be appropriate (see e.g. Van de Graaff, 1986). In a probabilistic approach an acceptable chance of failure must be the starting point instead of a single set of design conditions. Because of the high importance of the low-lying hinterland, a probability of failure for dunes of 10-5 per year (return period 100,000 years) has been agreed in The Netherlands for the most important parts of the country. This chance seems rather small, but compared to other threats (e.g. accidents with nuclear power plants, air planes or large industrial plants) the probability of failure is not that small taking into account the number of human lives and the high investments that are at stake.
+
The model simulations further show that the effectiveness of coastal sand nourishment to mitigate the impact of sea-level rise, strongly depends on the location in the profile where this sand is added. However, the ratio of the reduction in dune erosion volume 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 from a morphological perspective.
  
The probability of failure of 10-5 per year holds for the most important parts of The Netherlands; in some more rural areas larger probabilities of failure are the (legal) norm.
 
  
Figure 2 shows a schematic plot of the relationship between the rate of erosion RD (the distance RD is in Figure 1 the distance between the initial edge of the dune and the edge of the dune after the storm surge) and the frequency of exceedance. According to Figure 2 there is a chance of 10-5 per year that a rate of erosion of 85 m is reached or will be surpassed.
+
==Related articles==
 +
* [[Dune stabilisation]]
 +
* [[Light revetments built-in into artificial dunes]]
 +
* [[Natural causes of coastal erosion]]
 +
* [[Dealing with coastal erosion]]
 +
* [[Shoreface profile]]
 +
* [[Bruun rule]]
 +
* [[Risk and coastal zone policy: example from the Netherlands]]
 +
* [[Swash zone dynamics]]
  
The safety problem of large parts of The Netherlands (as far as the problem depends on the protection by dunes (large scale problem)), seems thus solvable if a proper insight is available in the possible occurring conditions in the (very) small chances range. However, this also introduces many uncertainties. E.g. reliable measurements of (maximum) water levels of the sea in measuring stations near the coast (cf. in ports) are available since approximately 1850. So for only 150 years. With such an in fact restricted data set, it is hard to estimate water levels which are associated with frequencies of exceedance in the range of 10-3 to 10-5 per year.
 
  
Data sets of reliable measurements of wave characteristics cover even shorter periods of time. Nevertheless, these kinds of uncertainties have been taken into account to arrive at figures like Figure 2.
+
==References==
Based on a legal norm, with the help of a proper computation method, and for a given situation (cross-section of cross-shore profile and dune area) one is next able to judge whether the cross-section meets the requirements; i.e. whether the cross-section is safe enough or not. If not: iteratively an improvement scheme for the dune area can be developed; e.g. a strengthening of the dunes at the landward side.
+
<references/>
  
===Small scale safety problem===
 
  
===Quantification of rate of dune erosion===
 
  
==See also==
+
{{2Authors
* [[Types and background of coastal erosion]]: article on the background of erosion, dune erosion and [[structural erosion]] (''Jan van de Graaff is also planning to write a separate article on dune erosion'').
+
|AuthorID1=11540
* [[Natural Causes of Coastal Erosion]]: Effects of e.g. transport gradient, loss of sand, protruding areas, marine deposit shorelines, down stream erosion, sea level rise, subsidence and natural variation on coastal erosion.
+
|AuthorFullName1= Jan van de Graaff
 +
|AuthorName1= Jan van de Graaff
 +
|AuthorID2=120
 +
|AuthorFullName2=Job Dronkers
 +
|AuthorName2=Dronkers J
 +
}}
  
{{author
+
[[Category:Coastal protection]]
|AuthorID=11540
+
[[Category:Physical coastal and marine processes]]
|AuthorName= Jan van de Graaff
+
[[Category:Beaches]]
|AuthorFullName= Jan van de Graaff}}
+
[[Category:Climate change, impacts and adaptation]]
[[Category:Theme 8]]
+
[[Category:Sea level rise]]

Revision as of 11:58, 5 July 2020

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.


Introduction

Fig. 1. Post-storm eroded dune, Egmond aan Zee.

Coastal dunes have developed naturally along many shorelines worldwide during the last millennia when sea-level 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]. The coastal dune belt in many cases protects low-lying 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 80-100 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 local conditions such as beach width, beach height and dune height.

If at least part of the dune belt survives the storm, it can recover by natural processes – the same processes that created the original dune. This will happen only if the conditions under which the dune belt was formed – sand supply, wave and wind climate – have not significantly changed.

Fig. 2. Storm-induced dune foot retreat as a function of frequency of exceedance for the Dutch coast.



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

Fig. 3. Badhotel at Schiermonnikoog after the 1924 storm.


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 cross-shore 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 [2], see Fig. 5.

Fig. 4. Dune failure mechanisms.
Fig. 5. Dune scarp after collapse. The people standing on the dune scarp have not read this article.


The sand deposited at the dune foot is removed by the seaward undercurrent ('undertow') under the breaking waves and deposited on the shoreface. The slope of the cross-shore profile gradually decreases, and consequently the rate of dune erosion will decrease with time during the storm surge. It is, however, not expected that a new equilibrium profile will develop because of the limited storm surge duration. The shape of the cross-shore 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 undoubtedly important in the dune erosion process, often a 2D approach is adopted. In this case the dune erosion process is considered as a typically offshore directed cross-shore 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 cross-shore direction. The same volume of sand which is eroded from the dunes and the very upper part of a cross-shore profile is accumulated lower in the cross-shore 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 sediment 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 short lasting process, some computation methods take only offshore directed 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 post-storm 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 post-storm equilibrium profiles established by laboratory experiments and validated by field data of the Dutch coast [3]. The post-storm 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:

  1. the storm duration is sufficient for establishment of a post-storm equilibrium profile [math]y(x)[/math];
  2. the post-storm equilibrium profile does not (strongly) depend on the initial coastal profile [math]y_0(x)[/math];
  3. eroded dune sediment is deposited within a zone delimited by a storm closure depth [math]y_{max}[/math];
  4. 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 parameterized functions for the post-storm equilibrium profile [math]y(x)[/math] and the storm closure depth [math]y_{max}[/math] are also indicated.


File:DUROSmodel.jpg
Fig. 6. Erosion profile method for estimating dune erosion as given by the beach volume between [math]y(x)[/math] and [math]y_0(x)[/math]. Underlying assumptions are explained in the text.


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, 2016)[4][5] 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. This force 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 du/dt \propto mu/T [/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] ([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=nf \propto u^4 \Delta t/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. The estimate of the dune volume loss [math]dV/dt[/math] per unit time obtained by the model is finally given by

[math]\Large\frac{dV}{dt}\normalsize=-4C_s\Large\frac{(R-z_D)^2}{T}\normalsize. \qquad (1)[/math]

The runup [math]R[/math] and the coefficient [math]C_s[/math] depend on local conditions. The 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 \propto \tan \beta \sqrt{HL} [/math], where [math]L \approx g T^2/2 \pi[/math] is the offshore wavelength. 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 approximate linear dependence of the runup [math]R[/math] on [math]\beta[/math] and [math]T[/math] implies that dune erosion is an increasing function of the beach slope and the wave period. 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). This is corroborated by a field study at Narrabeen-Collaroy Beach (East Australia, near Sidney), which showed that alongshore variability in dune erosion was strongly correlated with alongshore variability in [math]z_D[/math][6].

DUROSTA

Based on theoretical work and laboratory experiments, Steetzel (1993) [7] has developed the 2-dimensional DUROSTA computation model, in which actual sediment transports are calculated based on wave-integrated velocity and sand concentration profiles. Although in the mathematical description of the cross-shore sediment transport the intra-wave 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 cross-shore 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 [8] is an open-source process-based 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 (infragravity wave) transformation (generation, propagation and dissipation), wave-induced 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 2014-version 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 [9]. XBeach simulations further illustrated that the observed alongshore variation in dune erosion was steered primarily by the pre-storm 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.


Post-storm 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 [10][11][12]. After-storm recovery starts with the development of embryo dunes on the backshore and is stimulated by vegetation growth on these embryo dunes[13]. 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 30m 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 after-storm recovery of the shoreline position at Narrabeen-Collaroy Beach in East Australia, near Sidney. The recovery period after each storm-induced shoreline retreat is only a few months at this location. Shoreline recovery at this location is significantly correlated with onshore migration of nearshore sand bars[12][14].

Overall caption (optional)
Fig. 7. Variation of the dune foot position in a beach transect of the island Walcheren (Netherlands) over the period 1900-2010 (solid white line, based on annual surveys of Rijkswaterstaat [15]). The trend is indicated by the dashed line, showing ongoing retreat of the dune foot from 1900-1990. The large dune erosion by the exceptional 1953 storm surge has no influence on the long-term trend of dune foot retreat. Stabilization after 1990 is due to implementation of the current 'hold-the-line' coastal policy by sand nourishments.
Figure 8. Variation of the shoreline position in a transect of Narrabeen-Collaroy Beach (East Australia, near Sidney) over the period October 2004-May 2005. Strong shoreline retreat occurs during the storm periods indicated in light blue. After-storm recovery occurs during the periods in between. Figure redrafted after [12].


Impact of sea-level rise on dune erosion

The XBeach model predicts a linear relationship between dune erosion volume and sea level rise [16]. This is primarily due to the higher water level in front of the dune and not to changes in the significant short-wave or infragravity wave height. Changes in the offshore angle of wave incidence also affect the dune erosion volume. A 30° shift from shore-normal influences the dune erosion volume to the same extent as a 0.4 m sea-level 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 sea-level rise, strongly depends on the location in the profile where this sand is added. However, the ratio of the reduction in dune erosion volume 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 from a morphological perspective.


Related articles


References

  1. Anthony, E.J., Mrani-Alaoui. M. and Héquette, A. 2010. Shoreface sand supply and mid- to late Holocene aeolian dune formation on the storm-dominated macrotidal coast of the southern North Sea. Marine Geology 276: 100–104
  2. 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 (2007) 1–19.
  3. 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. Large-scale dune erosion tests to study the influence of wave periods Coastal Engineering 55: 1041–1051
  4. Larson, M., Erikson, L. and Hanson, H. 2004. An analytical model to predict dune erosion due to wave impact. Coastal Engineering 51: 675– 696
  5. Larson, M., Palalane, J., Fredriksson, F. and Hanson, H. 2016. Simulating cross-shore material exchange at decadal scale. Theory and model component validation. Coastal Engineering 116: 57–66
  6. 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
  7. Steetzel, H.J. 1993. Cross-shore transport during storm surges. Thesis Tech. Univ. Delft, Delft Hydraulics Communications 476.
  8. https://xbeach.readthedocs.io/en/latest/user_manual.html
  9. 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: 167-175.
  10. Vousdoukas, M.I., Pedro, L., Almeida, M. and Ferreira, O. 2012. Beach erosion and recovery during consecutive storms at a steep-sloping, meso-tidal beach. Earth Surf. Process. Landf. 37: 583–593
  11. 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
  12. 12.0 12.1 12.2 Phillips, M.S., Harley, M.D., Turner, I.L., Splinter, K.D. and Cox , R.J. 2017. Shoreline recovery on wave-dominated sandy coastlines: the role of sandbar morphodynamics and nearshore wave parameters. Marine Geology 385: 146–159
  13. Suanez, S., Cariolet, J-M., Cancouët, R., Ardhuin, F. and Delacourt, C. 2012. Dune recovery after storm erosion on a high-energy beach: Vougot Beach, Brittany (France). Geomorphology 139-140: 16–33
  14. 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
  15. Jeuken, C., Ruessink, G. and Marchand, M. 2001. Ruimtelijke en temporele aspecten van de duinvoetdynamiek. Report WL/Delft Hydraulics Z2838 (in Dutch)
  16. 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/s10584-017-1922-3


The main authors of this article are Jan van de Graaff and Job Dronkers
Please note that others may also have edited the contents of this article.

Citation: Jan van de Graaff; Job Dronkers; (2020): Dune erosion. Available from http://www.coastalwiki.org/wiki/Dune_erosion [accessed on 28-03-2024]