Difference between revisions of "Shoreface profile"

From Coastal Wiki
Jump to: navigation, search
(40 intermediate revisions by the same user not shown)
Line 1: Line 1:
 +
 +
{{Definition|title=Shoreface
 +
|definition= The shoreface is the nearshore zone of the inner continental shelf that is bounded landward by the mean low-water line and that extends seaward to where the influence of wave action on cross-shore sediment transport is on average minor compared to other influences. A similar definition is: The shoreface is the zone seaward of the shoreline where offshore generated waves interact with the upward sloping seabed. }}
  
 
{{Definition|title=Shoreface profile
 
{{Definition|title=Shoreface profile
|definition= The shoreface profile, often called beach profile, is the cross-shore coastal depth profile extending from the low-water line to the [[closure depth]]. }}
+
|definition= The shoreface profile, often called beach profile, is the cross-shore coastal depth profile of the shoreface.  }}
 +
 
 +
 
 +
 
 +
This article focuses on the shoreface of sandy coasts. The shoreface of rocky shores is described in the article [[Rocky shore morphology]] and the shoreface of mud coasts in the article [[Coastal mud belt]].
 +
 
 +
 
 +
 
  
  
==Introduction==
+
==Closure depth==
  
 +
[[File:ShorefaceUpperLower.jpg|thumb|400px|right|Fig. 1. Schematic representation of a shoreface profile. The vertical/horizontal scale ratio is greatly exaggerated. ]]
  
[[Image: ShorefaceProfile.jpg|thumb|400px|right|Figure 1. Schematic representation of a shoreface profile. The vertical/horizontal scale ratio is greatly exaggerated. ]]
+
In the literature one can find various definitions of the shoreface. The definition given here, which is most commonly used, follows the definition given in Stive and De Vriend (1995<ref name=S>Stive, M.J.F. and De Vriend, H.J. 1995. Modelling shoreface profile evolution. Mar. Geol. 126: 235–248</ref>). The shoreface can be divided in two zones, the upper shoreface and the lower shoreface. The upper shoreface is the zone where under average conditions waves break and the lower shoreface the zone where waves shoal. The lower part of the shoreface extends to the so-called outer [[closure depth]] <math>h_{out}</math>; beyond this depth the seabed is hardly influenced by waves and wave-induced sediment transport is (on average) insignificant. The lower part of the shoreface is also called shoaling zone; wave shoaling is the process of wave amplification when waves travel from deep to shallow water before breaking (see [[Shallow-water wave theory]]). The outer [[closure depth]] <math>h_{out}</math> has been empirically determined by Hallermeier (1981<ref>Hallermeier, R.J. 1981. A profile zonation for seasonal sand beaches from wave climate. Coast. Eng. 4: 253–277</ref>), yielding the formula (see also [[Closure depth]])
  
 +
<math>h_{out} \approx (H_s – 0.3 \sigma) \, T_s \sqrt{\large\frac{g}{5000 d_{50}}\normalsize} \approx 0.013 H_s T_s \sqrt{\large\frac{g}{d_{50}}\normalsize}  , \qquad (1)</math>
  
[[Image: DutchCoastNorthHolland.jpg|thumb|400px|right|Figure 2. The barred coast of North Holland, The Netherlands. https://beeldbank.rws.nl, Rijkswaterstaat / Rens Jacobs]]
+
where <math>H_s, T_s</math> are the average significant wave height and period, <math>\sigma</math> is the wave height standard deviation, <math>d_{50}</math> the medium sediment grainsize and <math>g</math> the gravitational acceleration. The second approximation in Eq. (1) follows from the assumption of a Rayleigh wave distribution (see [[Statistical description of wave parameters]]). Typical values of the outer closure depth are in the range 15 – 40 m.
  
The shoreface is the zone where offshore generated waves dissipate their energy, thereby suspending and transporting large quantities of sediment. Most energy is dissipated in the upper part of the shoreface, the [[Definitions of coastal terms#Breaker zone or Surf zone|surf zone]] (or breaker zone), where waves overturn and break. The lower part of the shoreface extends to the so-called [[closure depth]], the depth where the seabed is hardly influenced by waves and where wave-induced sediment transport is insignificant. The lower part of the shoreface is also called shoaling zone; wave shoaling is the process of wave amplification when waves travel from deep to shallow water before breaking (see [[Shallow-water wave theory]]). The division between upper and lower shoreface profiles is generally marked by an inflexion point, a bar or a terrace, indicating transition of  morphodynamics <ref>Wright, L.D., Chappel, J., Thom, B.G., Bradshow, M.P. and Cowell, P. 1979.  Morphodynamics of reflective and dissipative beach and inshore systems: south-eastern Australia. Mar.Geol. 32: 105-140</ref> (Fig. 1). Another transition exists at the shoreline where waves collapse and run up the beach, the so-called [[Swash zone dynamics|swash zone]]. The surf zone profile is often undulated due to the presence of one or several bars. These bars have mainly a shore parallel orientation, but more complex three-dimensional patterns are ubiquitous (see [[Rhythmic shoreline features]] and Fig. 2).   
+
The definition of the outer closure depth implies that wave-driven sediment transport across this boundary is very small on average. It does not rule out significant sediment transport across this boundary due to e.g. tidal currents, wind-driven currents, ocean currents or internal waves. For example, the occurrence of shoreface-connected ridges (see [[Sand ridges in shelf seas]]) on the lower shoreface can be theoretically related to sediment transport by strong wind-driven currents. Other definitions of the outer boundary of the lower shoreface can also be found in the literature, for example the so-called Depth of Transport. This is the depth at which the total bed shear stress does not exceed the threshold of sediment transport under extreme storms<ref name=V19>Valiente, N.G., Masselink, G., Scott, T., Conley, D. and McCarroll, R.J. 2019. Role of waves and tides on depth of closure and potential for headland bypassing. Mar. Geol. 407: 60–75</ref>. The Depth of Transport is typically further seaward than the outer closure depth Eq. (1).
 +
 
 +
Most energy is dissipated in the upper part of the shoreface, where wave overturning and breaking occurs. The depth where waves start breaking depends on the wave climate and is therefore highly variable. An often used definition of the transition depth <math>h_{in}</math> between upper and lower shoreface was also given by Hallermeier (1978<ref>Hallermeier, R.J. 1978. Uses for a calculated limit depth to beach erosion. In: 16th Coastal Engineering Conference, American Society of Civil Engineers, pp. 1493- 1512</ref>),
 +
 
 +
<math>h_{in} \approx 2.28 H_{12h/y}- 68.5 \Large\frac{ H_{12h/y}^2}{g T_{12h/y}^2}\normalsize , \qquad (2)</math>
 +
 
 +
where <math>H_{12h/y}</math> is the effective wave height just seaward of the breaker zone that is exceeded on average for 12 hours per year (i.e. the significant wave height with a probability of yearly exceedance of 0.137%) and <math>T_{12h/y}</math> the associated wave period. Later refinements of this formula are discussed in the article [[Closure depth]]. Typical values of <math>h_{in}</math> are in the range 3 – 10 m. The width of the upper shoreface typically ranges between about 100 - 1000 meters and the average slope between 1/10 and 1/100. The upper shoreface is sometimes also called breaker zone or surf zone, although onshore surfing wave bores usually only occur in the shoreward part of the upper shoreface. It is therefore preferable to use the designation upper shoreface as a fixed geographical indication, and the designations surf zone or breaker zone as a physical indication of the zone where broken waves surf onshore. This latter zone is narrow and close to the shoreline in a gentle wave climate and can be very wide under storm conditions, extending from the seaward boundary of the upper shoreface to the dune foot. Other subdivisions of the coastal profile include the [[active coastal zone]], the beach face or foreshore and the backshore, see [[Definitions of coastal terms]].
 +
 
 +
The division between upper and lower shoreface profiles is often marked by an inflexion point, a bar or a terrace, indicating transition of  morphodynamics <ref>Wright, L.D., Chappel, J., Thom, B.G., Bradshow, M.P. and Cowell, P. 1979.  Morphodynamics of reflective and dissipative beach and inshore systems: south-eastern Australia. Mar.Geol. 32: 105-140</ref><ref>Niedoroda, A.W., Swift, D.J., Hopkins, T.S. and Ma, C.-M. 1984. Shoreface morphodynamics on wave-dominated coasts. Mar. Geol. 60: 331–354</ref><ref>Larson, M., Kraus, N.C. and Wise, R.A. 1999. Equilibrium beach profiles under breaking and non-breaking waves. Coast. Eng. 36: 59–85</ref> (Fig. 1). Another indicator of the transition between upper and lower shoreface is a break in the trend of the sediment grainsize with depth<ref>Aragones, L., Pagan, J.I., Lopez, I. and Serra, J.C. 2018. Depth of closure: New calculation method based on sediment data. International Journal of Sediment Research 33:  198-207</ref>.  The shoreface is bordered on the landward side by the [[Definitions of coastal terms#Beach face|beach face]]. This is the zone where under average conditions waves collapse and run up the beach (see [[Swash zone dynamics]]). The upper shoreface profile is often undulated due to the presence of one or several bars, see [[Nearshore sandbars]]. These bars have mainly a shore parallel orientation, but more complex three-dimensional patterns are ubiquitous (see [[Rhythmic shoreline features]] and Fig. 3).   
  
  
 
==Morphodynamic feedback==
 
==Morphodynamic feedback==
  
Sediment transport on the shoreface depends on the local wave climate. By "wave climate" is meant the description of incoming waves in terms of a statistical distribution of incident wave directions, wave heights and wave periods (see [[Statistical description of wave parameters]]). Sediment transport also depends on the characteristics of bed sediments and on the seabed bathymetry. For small wave incidence angles (close to the shore-normal direction) the relevant seabed bathymetry is the cross-shore shoreface profile.  
+
Sediment transport on the shoreface depends on the local wave climate. By 'wave climate' is meant the description of incoming waves in terms of a statistical distribution of incident wave directions, wave heights and wave periods (see [[Statistical description of wave parameters]]). Sediment transport also depends on the [[Coastal and marine sediments|characteristics of bed sediments]] and on the seabed bathymetry. For small wave incidence angles (close to the shore-normal direction) the relevant seabed bathymetry is the cross-shore shoreface profile.  
  
The shoreface profile influences cross-shore sediment transport, but the inverse also holds: changes of the shoreface profile result from gradients in the cross-shore sediment transport. This is called morphodynamic feedback: for a given incident wave field the shoreface profile determines the transformation of the wave field on the shoreface and therefore the resulting sediment transport, whereas the resulting sediment transport determines the shoreface profile. As a result of this morphodynamic feedback, the shoreface profile tends to an equilibrium shape for a given stable wave climate. A dynamic equilibrium is reached when onshore directed sediment transport and offshore directed transport are equal on average over a period exceeding the morphodynamic adaptation time. This holds for alongshore uniform coasts; the equilibrium requirement for non-uniform coasts is a vanishing divergence of the long-term average total sediment transport vector (cross-shore and alongshore). It was realised already long ago that the shoreface of a sedimentary coast should not be regarded as a geologically inherited feature, but as the result of natural adaptation to modern hydrodynamic conditions. <ref>Cornaglia, P. 1889. Delle Spiaggie. Accademia Nazionale dei Lincei, Atti.Cl.Sci.Fis., Mat.e Nat.Mem. 5: 284-304</ref> <ref>Johnson, D.W. 1919. Shore processes and shoreline development. Prentice Hall, N-Y, 584 pp.</ref>.
+
[[File:GoldCoastCyclone.jpg|thumb|450px|right|Fig. 2. Gold Coast, Australia, before and after the June 1967 cyclone. Observed shoreface profiles showing recovery of the initial close-to-equilibrium profile after strong perturbation by an extreme event. Redrawn after Patterson and Nielsen (2016<ref>Patterson, D.C. and Nielsen, P. 2016. Depth, bed slope and wave climate dependence of long term average sand transport across the lower shoreface. Coastal Engineering 117: 113–125</ref>).]]
 +
 
 +
The shoreface profile influences cross-shore sediment transport, but the inverse is also true: changes of the shoreface profile result from gradients in the cross-shore sediment transport. This is called morphodynamic feedback: for a given incident wave field the shoreface profile determines the transformation of the wave field on the shoreface and therefore the resulting sediment transport, whilst the resulting sediment transport determines the shoreface profile. As a result of this morphodynamic feedback, the shoreface profile tends to an equilibrium shape for a given stable wave climate<ref>Wright, L.D., Short, A.D. and Green, M.O. 1985. Short-term changes in the morphodynamic states of beaches and surf zones: an empirical predictive model. Marine Geology 62: 339—364</ref><ref>Ludka, B. C., Guza, R. T. O’Reilly, W. C. and Yates M. L. 2015. Field evidence
 +
of beach profile evolution toward equilibrium, J. Geophys. Res. Oceans 120: 7574–7597</ref><ref>Eichentopf, S., Van der Zanden, J., Caceres, I. and Alsina, J.M. 2019. Beach profile evolution towards equilibrium from varying initial morphologies. J. Mar. Sci. Eng. 7, 406</ref>. An illustration is presented in Fig. 2, which shows the recovery after strong perturbation of a shoreface profile that is initially close to equilibrium.
 +
 
 +
A dynamic equilibrium is reached when onshore directed sediment transport and offshore directed transport are equal on average over a period exceeding the morphodynamic adaptation time. This holds for alongshore uniform coasts; the equilibrium requirement for non-uniform coasts is a vanishing divergence of the long-term average total sediment transport vector (cross-shore and alongshore). It was realised already long ago that the shoreface of a sedimentary coast should not be regarded as a geologically inherited feature, but as the result of natural adaptation to modern hydrodynamic conditions and sediment supply<ref>Cornaglia, P. 1889. Delle Spiaggie. Accademia Nazionale dei Lincei, Atti. Cl. Sci. Fis., Mat. e Nat. Mem. 5: 284-304</ref><ref>Johnson, D.W. 1919. Shore processes and shoreline development. Prentice Hall, N-Y, 584 pp.</ref>. In real field situations, a coast never assumes its equilibrium profile; it is always in a transient state due to fluctuations in the wave climate and the mean water level. The coastal profile deviates from the equilibrium profile, even in the hypothetical case of a stationary wave climate, due to the cyclical or quasi-cyclical migration of seabed structures along the profile (see [[Nearshore sandbars]]).
  
 
Knowing how to determine the equilibrium shoreface profile has great practical interest. It provides quantitative insight in the response of the shoreface to changes in the local wave climate by human and natural causes; such insight is important for testing and improving the effectiveness of measures to combat coastal erosion and shoreline retreat.  
 
Knowing how to determine the equilibrium shoreface profile has great practical interest. It provides quantitative insight in the response of the shoreface to changes in the local wave climate by human and natural causes; such insight is important for testing and improving the effectiveness of measures to combat coastal erosion and shoreline retreat.  
Line 28: Line 53:
 
Morphodynamic feedback involves not only shoreface profile development, but also sediment grainsize and bed roughness. Grainsize and bed roughness strongly influence sediment transport <ref> Van Rijn, L.C. 1998. Principles of Coastal Morphology. Aqua Publications, The Netherlands  (aqua publications.nl.</ref>. Inversely, sediment transport strongly influences the grainsize distribution of bed sediments. Bed roughness depends on the presence of bedforms (ripples, dunes), which in turn depend on sediment transport and sediment grainsize. Consequently, shoreface profile development, wave transformation on the shoreface, sediment transport, grainsize sorting and bedform emergence are all dynamically coupled processes.   
 
Morphodynamic feedback involves not only shoreface profile development, but also sediment grainsize and bed roughness. Grainsize and bed roughness strongly influence sediment transport <ref> Van Rijn, L.C. 1998. Principles of Coastal Morphology. Aqua Publications, The Netherlands  (aqua publications.nl.</ref>. Inversely, sediment transport strongly influences the grainsize distribution of bed sediments. Bed roughness depends on the presence of bedforms (ripples, dunes), which in turn depend on sediment transport and sediment grainsize. Consequently, shoreface profile development, wave transformation on the shoreface, sediment transport, grainsize sorting and bedform emergence are all dynamically coupled processes.   
  
Grainsize sorting is primarily due to differences in transport modes of different grainsize fractions. The finer fractions are more easily brought and maintained in suspension than the coarser fractions. They are transported further away from their initial location and tend to settle on the lower shoreface zone where wave-induced bed shear stresses are lower than on the upper shoreface. The coarse fraction is transported partly as suspended load and partly as bedload in the wave boundary layer. Skewed shallow-water waves favour net onshore transport in the wave boundary layer. Sediments on the upper shoreface are therefore coarser than on the lower shoreface. Fine sediment particles can still be present on the upper shoreface, but they are hidden by overlying coarser sediment particles – a phenomenon called 'bed armouring'.  
+
Grainsize sorting is primarily due to differences in transport modes of different grainsize fractions. The finer fractions are more easily brought and maintained in suspension than the coarser fractions. They are transported further away from their initial location and tend to settle on the lower shoreface zone where wave-induced bed shear stresses are lower than on the upper shoreface. The coarse fraction is transported partly as suspended load and partly as bedload in the wave boundary layer. Wave transformation in shallow water leading to wave skewness (stronger onshore than offshore wave orbital velocity) and wave asymmetry (stronger onshore than offshore wave orbital acceleration) favour net onshore transport in the wave boundary layer. Fine sediments are resuspended from the upper shoreface due to [[undertow]] currents, longshore currents and turbulent shear stresses generated by breaking waves. Sediments on the upper shoreface are therefore coarser than on the lower shoreface. Fine sediment particles can still be present on the upper shoreface if they are 'hidden' by overlying coarser sediment particles – a phenomenon called 'bed armoring' (see [[Seabed armoring]]).
 +
 
 +
[[Image: DutchCoastNorthHolland.jpg|thumb|400px|right|Fig. 3. The barred coast of North Holland, The Netherlands. https://beeldbank.rws.nl, Rijkswaterstaat / Rens Jacobs]]
 +
 
 +
The largest bedforms on the shoreface, such as shoreface-connected ridges (scale of several kilometers) and shore-oblique bars (scale of the order of hundred meters) are mainly related to morfodynamic feedback processes that interact with longshore transport, see the articles [[Sand ridges in shelf seas]] and [[Rhythmic shoreline features]]. Other large bedforms, the so-called nearshore sandbars or breaker bars (scale also of the order of 100 meters, Fig. 3) arise from morphodynamic coupling with the wave breaking process, through the resulting flow structure at the breaker bar and interaction with cross-shore and longshore currents <ref>Reniers, A.J.H.M., Thornton, E.B., Stanton, T.P. and Roelvink, J.A. 2004. Vertical flow structure during Sandy Duck: observations and modelling. Coast.Eng. 51: 237-260</ref><ref>Jacobsen, N.G. and Fredsoe, J. 2014. Formation and development of a breaker bar under regular waves. Part 2: Sediment transport and morphology. Coastal Eng. 88: 55-68</ref> - see [[Nearshore sandbars]] for a more detailed discussion. Bedforms at much smaller scales – ripples with wavelengths of the order of 0.3 - 1 m – are formed in the less exposed parts of the shoreface under low to moderate wave conditions. They are often found on the lower shoreface and in the troughs between nearshore breaker bars. Morphodynamic feedback processes, leading to regular ripple patterns, involve ripple-induced near-bed flow circulation<ref> Sleath, J.F.A. 1984. Sea bed mechanics. Wiley, New York. </ref>, vortex shedding at ripple crests<ref> Fredsøe, J., Andersen, K.H. and Sumer, M.B. 1999. Wave plus current over a ripple-covered bed. Coastal Eng. 38: 177-221 </ref> and non-linear ripple-ripple interactions<ref> Marieu, V., Bonneton, P., Foster, D.L. and Ardhuin, F.  2008. Modeling of vortex ripple morphodynamics. J. Geophys. Res. 113, C09007, doi:10.1029/2007JC004659 </ref>. For a more detailed discussion, see the articles [[Wave ripples]] and [[Wave ripple formation]]. The time lag of sediment suspension and settling related to vortex shedding at the ripple crests strongly influences the direction of residual cross-shore sediment transport on the shoreface, see section [[#Onshore wave-induced sediment transport]].  
 +
<br clear=all>
 +
 
  
The largest bedforms on the shoreface, such as shoreface-connected ridges (scale of several kilometers) and shore-oblique bars (scale of the order of hundred meters) are mainly related to longshore transport through morfodynamic feedback processes, see the articles [[Sand ridges in shelf seas]] and [[Rhythmic shoreline features]]. Other large bedforms, the so-called breaker bars (scale also of the order of 100 meters) arise from morphodynamic coupling with the wave breaking process, through the resulting flow structure at the breaker bar <ref>Reniers, A.J.H.M., Thornton, E.B., Stanton, T.P. and Roelvink, J.A. 2004. Vertical flow structure during Sandy Duck: observations and modelling. Coast.Eng. 51: 237-260</ref><ref>Jacobsen, N.G. and Fredsoe, J. 2014. Formation and development of a breaker bar under regular waves. Part 2: Sediment transport and morphology. Coastal Eng. 88: 55-68</ref>. Bedforms at much smaller scales – ripples with wavelengths of the order of 1 m – are formed in the less exposed parts of the shoreface under low to moderate wave conditions. Morphodynamic feedback processes, leading to regular ripple patterns, involve ripple-induced near-bed flow circulation <ref> Sleath, J.F.A. 1984. Sea bed mechanics. Wiley, New York. </ref>, vortex shedding at ripple crests <ref> Fredsøe, J., Andersen, K.H. and Sumer, M.B. 1999. Wave plus current over a ripple-covered bed. Coastal Eng. 38: 177-221 </ref> and non-linear ripple-ripple interactions <ref> Marieu, V., Bonneton, P., Foster, D.L. and Ardhuin, F.  2008. Modeling of vortex ripple morphodynamics. J. Geophys. Res. 113, C09007, doi:10.1029/2007JC004659 </ref>. For a more detailed explanation, see the articles [[Wave ripples]] and [[Wave ripple formation]]. The time lag of sediment suspension and settling related to vortex shedding at the ripple crests strongly influences the direction of residual sediment transport on the shoreface.
+
==Beach classification; dissipative and reflective==
  
 +
The beach profile in the upper shoreface zone can be highly variable. The profile tends to steepen in periods of low-energy waves (especially long-period swell waves) and to flatten under high-energy storm waves (see [[Active coastal zone]]). Sediment characteristics, which largely depend on the supply source (see [[Coastal and marine sediments]]), also play an important role. Coarse-grained beaches typically have steeper slopes than fine-grained beaches<ref>Bujan, N., Cox, R. and Masselink, R. 2019. From fine sand to boulders: Examining the relationship between beach-face slope and sediment size. Marine Geology 417, 106012 </ref>. Gently sloping beaches dissipate almost completely the energy of incident waves, whereas steep-sloping beaches tend to reflect (partly) incident waves.  For distinguishing between reflective and dissipative beaches the Dean-parameter <math>\Omega</math> is often used <ref> Gourlay, M. R. 1968. Beach and dune erosion tests. Rep. m935/m936, Delft Hydraul. Lab., Delft </ref><ref> Dean, R.G. 1973. Heuristic models of sand transport in the surf zone. Proc. Conf. Eng. Dynamics in the Surf Zone, Sydney: 208-214. </ref>:
  
==Beach classification==
+
<math>\Omega = \Large \frac{H_b}{wT_p} \normalsize , \qquad(3)</math>
  
The beach profile in the upper shoreface zone can be highly variable. The profile tends to steepen in periods of low-energy waves (especially long-period swell waves) and to flatten under high-energy storm waves. Sediment characteristics also play an important role. Coarse-grained beaches, typical for high-energy coasts, have steeper slopes than fine-grained beaches, typical for moderate-to-low-energy coasts. Gently sloping beaches dissipate almost completely the energy of incident waves, whereas steep sloping beaches tend to reflect incident waves - at least partly.  For distinguishing between reflective and dissipative beaches the Dean-parameter <math>\Omega</math> is often used <ref> Gourlay, M. R. 1968. Beach and dune erosion tests. Rep. m935/m936, Delft Hydraul. Lab., Delft </ref><ref> Dean, R.G. 1973. Heuristic models of sand transport in the surf zone. Proc. Conf. Eng. Dynamics in the Surf Zone, Sydney: 208-214. </ref>:
+
where <math>H_b</math> is the significant wave height at breaking, <math>T_p</math> is the peak spectral wave period and <math>w</math> is the mean sediment fall velocity.
  
<math>\Omega = \Large \frac{H}{wT} \normalsize , \qquad(1)</math>
+
[[File:BeachClassification.jpg|thumb|right|400px|Fig. 4. Beach morphologies according to the parameters <math>\Omega</math> and <math>RTR</math>, based on beach surveys in England, Wales, Northern Ireland, Australia and China<ref name=S11>Scott, T., Masselink, G. and Russell, P. 2011. Morphodynamic characteristics and classification of beaches in England and Wales. Mar. Geol. 286: 1–20</ref><ref name=S13>Short, A.D. and Jackson, D.W.T. 2013. Beach morphodynamics. In: Shroder, J.F. (Ed.), Treatise in Geomorphology, vol. 10. Academic Press, San Diego</ref><ref name=J22>Jackson, D.W.T., Short, A.D., Loureiro, C. and Cooper, J.A.G. 2022. Beach morphodynamic classification using high-resolution nearshore bathymetry and process-based wave modelling. Estuarine, Coastal and Shelf Science 268, 107812</ref><ref>Li, Y., Zhang, C., Zhao, S., Qi, H., Cai, F. and Zheng, J. 2024. Equilibrium configurations of sandy-muddy transitional beaches on South China coasts: Role of waves in formation of sand-mud transition boundary. Coastal Engineering 187 (2024) 104401</ref>. Reflective beaches are found in the red zone and dissipative beaches in the blue zone.]] 
  
where <math>H</math> is the offshore wave height before breaking, <math>T</math> is the peak spectral wave period and <math>w</math> is the sediment fall velocity. By analysing beach profiles in Australia, Wright and Short (1984) <ref>Wright, L.D. and Short, A.D. 1984. Morphodynamic variability of surf zones and beaches: a synthesis. Mar.Geol. 56: 93-118</ref> found that beach types could be characterised by the parameter <math>\Omega</math>:
+
Wright and Short (1984) <ref>Wright, L.D. and Short, A.D. 1984. Morphodynamic variability of surf zones and beaches: a synthesis. Mar.Geol. 56: 93-118</ref> analysed beach profiles in Australia by evaluating the parameter <math>\Omega</math>, considering the local long-term average wave climate. They found that beach types could be classified as follows :
  
 
* reflective beaches correspond to <math>\Omega < 1 </math>;
 
* reflective beaches correspond to <math>\Omega < 1 </math>;
Line 45: Line 77:
 
* dissipative beaches correspond to <math>\Omega > 6 </math>.
 
* dissipative beaches correspond to <math>\Omega > 6 </math>.
  
The tidal range also influences the beach profile, because the surf zone is shifted along the shoreface with tidal level. For large tides (tidal range much larger than wave height) the surf zone covers during the tidal cycle a much larger area than for small tides. Beaches with large tides therefore have smaller upper shoreface slopes compared to beaches with the same <math>\Omega</math> and small tides.     
+
The tidal range <math>H_{tide}</math> also influences the beach profile, because the surf zone is shifted across the shoreface with tidal level. For large tides (tidal range much larger than wave height) the surf zone covers during the tidal cycle a much larger area than for small tides. Beaches with large tides therefore have smaller upper shoreface slopes compared to beaches with the same <math>\Omega</math> and small tides.     
 +
 
 +
Several field studies<ref name=S11/><ref name=S13/><ref name=J22/> have been performed to relate beach profiles to different combinations of the Dean parameter <math>\Omega</math> and the relative tidal range <math>RTR=H_{tide}/H_b</math>. Distinction is made between different beach morphologies: reflective beaches with steep graveled slopes, reflective/intermediate beaches with [[Nearshore sandbars|subtidal bars]], intermediate beaches with [[Rip current|rip cells]] and low-tide terrace or low-tide bar, dissipative beaches without significant bar formations, dissipative beaches with multiple inter-tidal bars and dissipative beaches with tidal flats, see Fig. 4. The figure only gives a rough indication of beach morphologies, because the natural variability is not fully captured by the two dimensionless parameters <math>\Omega</math> and <math>RTR</math>. Other factors play a role, in particular the temporal variability in  morphology and environmental conditions and the role of geology<ref name=S11/>. 
 +
 
 +
Another parameter that is often used for analysing the beach response to variations in the wave climate is the so-called [[surf similarity parameter]] that relates the wave steepness (wave height to wave length ratio) to the surf zone slope. 
 +
 
 +
 
 +
==Difference between upper and lower shoreface==
 +
The upper shoreface is the zone where the energy of incident waves is converted into water level gradients, currents, seabed agitation and sediment transport. This goes with strong morphodynamic activity manifested in profile change, generation of bedforms and shoreline advance or retreat.
 +
In contrast, wave energy dissipation on the lower shoreface is small and hardly affects the wave energy flux toward the coast, <math>F^{(x)} = c_g H_s^2 \cos \alpha</math>, which can therefore be considered constant (Green's law; see Table Symbols for the symbols). Assuming approximate validity of linear wave theory, <math>c_g \approx c \approx \sqrt{gh}</math>, the height of waves propagating to shallower waters increases as <math>H_s \propto h^{-1/4}</math>. While shoaling, these waves develop wave skewness and asymmetry on the lower shoreface (see Appendix). Wave skewness and asymmetry generate net sediment transport (see next paragraph) and therefore play an important role in shaping the cross-shore profile.
 +
 
 +
Sediment transport rates on the lower shoreface are much smaller than on the upper shoreface, implying a long adaptation time scale to changes in sea level and wave climate. On the upper shoreface these time scales are on the order of day-year, compared to decade-millennium time scales on the lower shoreface<ref name=S/><ref name=A20> Anthony, E.J. and Aagaard, T. 2020. The lower shoreface: morphodynamics and sediment connectivity with the upper shoreface and beach. Earth Sci. Rev. 210: 103334</ref>. Besides, waves are not the only sediment transport agent on the lower shoreface; tidal currents and wind-driven currents are important as well. The strongest currents are associated with extreme events (storm surges, cyclones), but persistent wind stress over longer periods can produce up or downwelling and associated up or downslope sediment transport. At places where river plumes flowing along the coast induce a cross-shore salinity gradient, fine sediment can be trapped on the lower shoreface and even on the upper shoreface, as explained in the article [[Coastal mud belt]]. However, on exposed coasts, energetic waves will transport fine sediment from the upper to the lower shoreface. The lower shoreface can also receive sediment directly from nearby river deltas. When sediment sources are absent or sparse, the seabed sediment of the lower shoreface has mainly a relict origin with possibly rocky outcrops, indurated clay or coarse terrigenous deposits<ref name=A20/>. Biotic activity on the lower shoreface can be intense with a high production of calcareous bioclasts. The lower shoreface can therefore act as a source of calcareous sand to the upper shoreface and the beach<ref>Short, A.D. 2013. Australia’s temperate carbonate coast: Sources, depositional environments and implications. In: Martini, I.P., Wanless, H.R. (Eds.), Sedimentary Coastal Zones from High to Low Latitudes: Similarities and Differences. London, Special Publications, Geological Society</ref>.
  
  
 
==Cross-shore sediment transport==
 
==Cross-shore sediment transport==
  
Sediment transport on the shoreface depends on many processes. We focus on processes relevant for cross-shore transport. It should be realised, however, that longshore transport processes can also play an important role. Shifts in the cross-shore position of the shoreface profile and the generation of bars on the shoreface are strongly influenced by longshore transport processes. See for example the articles: [[Natural causes of coastal erosion]], [[Littoral drift and shoreline modelling]] and [[Rhythmic shoreline features]].
+
Direct observation of net cross-shore sediment transport and its separate components in the field is difficult. This is particularly true for the surf zone, where the temporal and spatial variability of transport processes is very high. The hydrodynamic conditions for the installation and operation of measuring devices are harsh, especially in periods where sediment transport is most significant. Bottom-mounted pressure sensors are most often employed for the determination of wave characteristics. This can easily yield erroneous (underestimated) values if the nonlinearity of wave dynamics is not accounted for<ref>Martins, K., Bonneton, P., Mouragues, A. and Castelle, B. 2020. Non-hydrostatic, non-linear processes in the surf zone. Journal of Geophysical Research: Oceans 125,
In the following we focus on sandy beaches, with grain sizes typically between 0.2 and 0.5 mm (fine to medium sand). Coarse sedimentary beaches are dealt with in the article [[Gravel Beaches]].
+
e2019JC015521</ref>. Acoustic and optical sensors are sensitive to air bubbles. Their use for determining suspended sediment concentrations and grainsizes is hampered as breaking waves entrain large volumes of air<ref name=AH>Aagaard, T. and Hughes, M.G. 2022. Sediment Transport. Chapter 8.6, Treatise on Geomorphology 2nd edition. Elsevier</ref>. For these reasons, most knowledge is derived from laboratory experiments. Knowledge from the field is often obtained indirectly, through interpretation of the observed erosion and accretion patterns produced by cross-shore sediment transport. Mathematical models  are often used as interpretation tools.
 +
 
 +
Sediment transport on the shoreface depends on many processes. The focus is on processes relevant to cross-shore transport. It should be realised, however, that longshore transport processes can also play an important role, for example by influencing the sediment balance of the shoreface or by generating various kinds of bedforms, see for example the articles: [[Natural causes of coastal erosion]], [[Littoral drift and shoreline modelling]] and [[Rhythmic shoreline features]].  
 +
In the following we consider sandy beaches, with grain sizes typically between 0.2 and 0.5 mm (fine to medium sand). Coarse sedimentary beaches are dealt with in the article [[Gravel Beaches]].
 +
 
 +
A distinction is often made between bedload transport and suspended load transport. Bedload transport represents the motion of sediment particles that are dragged by the flow along the seabed, either by rolling or successive jumps (saltation) or by entrainment of a thin seabed top layer (so-called [[sheet flow]] layer, of typically 10-30 grains thick). Suspended load represents the transport of sediment particles that have been brought into suspension by turbulent fluid motions. An important difference between the two transport modes is the phase between the load of the transported sediment and the instantaneous flow velocity, especially in the case of wave-induced flow. Bedload responds almost instantaneously to changes in the wave orbital velocity; increase (decrease) of velocity implies immediate increase (decrease) of bedload. This implies that an onshore (offshore ) flow asymmetry (so-called velocity skewness) with higher peak velocity in onshore (offshore) flow induces a net bedload transport in onshore (offshore) direction. Suspended load generally displays a fast response to an increase in flow velocity; sediment particles can be lifted quickly from the seabed and brought into suspension when the flow velocity increases. However, unlike bedload, the response of the suspended sediment concentration to a decrease in flow velocity depends on the settling velocity of sediment particles, which induces a response delay that can be substantial compared to the wave period. Asymmetries in onshore (offshore) wave orbital flow are therefore less decisive for the net direction of suspended load than for bedload. In the surf zone where the waves break, suspended load transport dominates. Flume experiments suggest that bedload transport dominates in situations where unbroken waves approach the shore, e.g. on the lower shoreface and on steep, reflective beaches<ref>van der Zanden, J., van der A, D. A., Hurther, D., Cáceres, I., O'Donoghue, T., Hulscher, S. J. M. H. and Ribberink, J. S. 2017. Bedload and suspended load contributions to breaker bar morphodynamics. Coastal Engineering 129: 74–92</ref>.
 +
 
 +
As explained below, onshore transport on the shoreface is mainly related to bedload transport and offshore transport to suspended load transport. As bedload transport involves mainly the coarse sediment fraction whilst suspended load transport involves mainly the fine sediment fraction, it follows that under average conditions coarser sediment is preferentially transported onshore and finer sediment offshore<ref>Hassan, W.N. and Ribberink, J.S. 2005. Transport processes of uniform and mixed sands in oscillatory sheet flow. Coastal Engineering 52: 745– 770</ref>.  
  
 
===Onshore wave-induced sediment transport===
 
===Onshore wave-induced sediment transport===
Wave transformation in shallow water (see [[Shallow-water wave theory]]) is a major factor for onshore sediment transport <ref>Elgar, S., Gallagher, E.L. and Guza, R.T. 2001.Nearshore sandbar migration. J.Geophys.Res. 106: 11,623-11,727</ref>. Waves become asymmetric and skewed, producing near-bed sheer stresses which are stronger in forward (onshore) direction than in backward (offshore) direction <ref>Nielsen, P. 1992. Coastal bottom boundary layers and sediment transport. In: Advanced Series on Ocean Engineering, IV. World Scientific.</ref> <ref>Abreu, T., Silva, P.A., Sancho, F. and Temperville, A. 2010. Analytical approximate wave form for asymmetric waves. Coastal Eng. 57: 656-667</ref><ref>Van der A, D.A., O’Donoghue, T., Davies, A.G. and Ribberink, J.S. 2011. Experimental study of the turbulent boundary layer in acceleration-skewed oscillatory flow. J. of Fluid Mech. 684: 251-283</ref>. Phase lags between bed sheer stress and sediment suspension determine to which degree wave skewness induces onshore transport; in some situations (rippled bed, sheet flow) a large phase lag can even reverse the net sediment transport direction <ref>Vincent, C.E. and Green. M.O. 1990. Field measurements of the suspended sand concentration profiles and fluxes and of the resuspensionm coefficient over a rippled bed. J.Geophys.Res. 95: 11591-11601</ref><ref>Ribberink, J.S. and Al-Salem, A.A. 1995. Sheet flow and suspension of sand in oscillatory boundary layer. Coastal Eng. 25: 205-225</ref><ref>Ruessink, B.G., Michallet, H., Abreu, T., Sancho, F., van der A, D.A., van der Werf, J.J. and Silva, P.A., 2011. Observations of velocities, sand concentrations, and fluxes under velocity-asymmetric oscillatory flows. J. Geophys. Res. 116, C03004, doi:10.1029/2010JC006443</ref>. The development of forward streaming at the top of the wave boundary layer also contributes to onshore transport <ref>Longuet-Higgins, M.S. 1953. Mass transport in water waves. Royal Soc. London, Phil.Trans. 245A: 535-581</ref>, but this only holds for smooth beds and not in cases where the seabed is rippled, in cases of strongly skewed waves  or for sheet flow conditions <ref>Trowbridge, J.H. and Madsen, O.S. 1984. Turbulent wave boundary layers: 2. Second-order theory and mass transport. J.Geophys.Res. 89: 7999-8007</ref><ref>Kranenburg, W. M., Ribberink, J. S., Uittenbogaard, R. E. and Hulscher, S. J. M. H. 2012. Net currents in the wave bottom boundary layer: on wave shape streaming and progressive wave streaming. J. Geophys. Res. 117(F03005), doi:10.1029/2011JF002070</ref>. A lesser contribution to onshore transport is further due to net mass transport in the zone between wave trough and wave crest.  
+
Wave transformation in shallow water is a major factor for onshore sediment transport <ref>Elgar, S., Gallagher, E.L. and Guza, R.T. 2001.Nearshore sandbar migration. J.Geophys.Res. 106: 11,623-11,727</ref>. Due to shoaling effects, starting on the lower shoreface, incident waves become increasingly skewed (asymmetry in onshore / offshore peak flow velocity) and asymmetric (asymmetry in onshore / offshore flow acceleration, also called acceleration skewness), see the Appendix for details. These asymmetries in the wave orbital velocity produce shear stresses on the seabed which are stronger in onshore direction than in offshore direction <ref name=N>Nielsen, P. 1992. Coastal bottom boundary layers and sediment transport. In: Advanced Series on Ocean Engineering, IV. World Scientific.</ref><ref name=A10>Abreu, T., Silva, P.A., Sancho, F. and Temperville, A. 2010. Analytical approximate wave form for asymmetric waves. Coastal Eng. 57: 656-667</ref><ref>Van der A, D.A., O’Donoghue, T., Davies, A.G. and Ribberink, J.S. 2011. Experimental study of the turbulent boundary layer in acceleration-skewed oscillatory flow. J. of Fluid Mech. 684: 251-283</ref>, see also [[Shallow-water wave theory]]. Wave skewness and asymmetry are primarily effective for inducing onshore bedload transport rather than suspended load transport, since bedload transport responds almost instantaneously to changes in near-bed velocities. Wave skewness and asymmetry are most pronounced on gently sloping dissipative coasts; they hardly develop on steep reflective coastal profiles. 
 +
 
 +
[[File:RippleRegimeTransportCompressed.gif|thumb|right|400px|Fig. 5. Sand transport over a steep vortex bed ripple in the Aberdeen oscillatory flow tunnel. The horizontal orbital motion is skewed: the velocity to the right (wave propagation direction) is stronger than to the left, but the duration is shorter. The high orbital velocity to the right creates a sand-laden vortex on the lee-side of the ripple crest, which is swept back over the ripple right after flow reversal. Courtesy: Jebbe van der Werf.]]
 +
 
 +
Phase lags between bed shear stress and suspended sediment concentration determine to which degree wave skewness induces onshore transport. In some situations (rippled bed, fine sands) a large phase lag can reverse the net direction of suspended sediment transport. A large time lag implies that sediment particles that are stirred up by turbulent eddies during onshore orbital motion are transported offshore after reversal of the wave orbital velocity<ref>Vincent, C.E. and Green. M.O. 1990. Field measurements of the suspended sand concentration profiles and fluxes and of the resuspension coefficient over a rippled bed. J.Geophys.Res. 95: 11591-11601</ref><ref>Ribberink, J.S. and Al-Salem, A.A. 1995. Sheet flow and suspension of sand in oscillatory boundary layer. Coastal Eng. 25: 205-225</ref><ref name=N></ref><ref>Ruessink, B.G., Michallet, H., Abreu, T., Sancho, F., van der A, D.A., van der Werf, J.J. and Silva, P.A., 2011. Observations of velocities, sand concentrations, and fluxes under velocity-asymmetric oscillatory flows. J. Geophys. Res. 116, C03004, doi:10.1029/2010JC006443</ref>, see Fig. 5. Although strongly skewed waves promote onshore ripple migration, the net sediment transport can be directed offshore. Onshore transport is stimulated by acceleration skewness<ref>Dubarbier, B., Castelle, B., Marieu, V. and Ruessink, G. 2015. Process-based modelling of cross-shore sandbar behavior. Coastal Engineering 95: 35–50</ref>. It creates a sediment-laden vortex on the stoss side of bed ripples, which travels toward the onshore direction as it is carried by the onshore ambient flow<ref>Salimi-Tarazouj, A., Tian-Jian Hsu, T-J., Traykovski, P. and Chauchat, J. 2024. Investigating wave shape effect on sediment transport over migrating ripples using an Eulerian two-phase model. Coastal Engineering 189, 104470</ref>. Wave breaking that generates turbulent fluid motions and sediment resuspension under the wave crest also contributes to onshore transport, especially in the case of plunging waves<ref>Aagaard, T., Hughes, M.G. and Ruessink, B.G. 2018. Field observations of turbulence, sand suspension and cross-shore transport under spilling and plunging breakers. Journal of Geophysical Research - Earth Surface 123: 2844–2862</ref>. But at the same time, wave breaking removes most of the wave skewness<ref>Guza, R.T. and Thornton, E.B. 1985. Velocity moments in the nearshore. J. Waterways Port Coastal Ocean Div. ASCE III: 235-256</ref>, thus reducing the associated onshore sediment transport component in the surf zone.  
 +
 
 +
The development of forward streaming at the top of the wave boundary layer can provide a further contribution to onshore transport <ref>Longuet-Higgins, M.S. 1953. Mass transport in water waves. Royal Soc. London, Phil.Trans. 245A: 535-581</ref>, but this only holds for smooth beds and not in cases where the seabed is rippled, in cases of strongly skewed waves  or for sheet flow conditions<ref>Trowbridge, J.H. and Madsen, O.S. 1984. Turbulent wave boundary layers: 2. Second-order theory and mass transport. J.Geophys.Res. 89: 7999-8007</ref><ref>Kranenburg, W. M., Ribberink, J. S., Uittenbogaard, R. E. and Hulscher, S. J. M. H. 2012. Net currents in the wave bottom boundary layer: on wave shape streaming and progressive wave streaming. J. Geophys. Res. 117(F03005), doi:10.1029/2011JF002070</ref> (see also [[Sediment transport formulas for the coastal environment]]). Due to the high sediment concentration in the streaming layer (and in spite of the small vertical extent of the layer and the low flow velocity), the contribution of wave driven near-bed streaming to onshore transport is significant, especially in the shoaling zone and surf zone prior the plunging point<ref>Zheng, P., Gumbira, G., Li, M., van der Zanden, J., van der A, D., van der Werf, J., Chen, X. and Tang, X. 2023. Development, calibration and validation of a phase-averaged model for cross-shore sediment transport and morphodynamics on a barred beach. Continental Shelf Research 258, 104989</ref>.  
  
 
===Offshore wave-induced sediment transport===
 
===Offshore wave-induced sediment transport===
Offshore transport is due to so-called undertow: the compensating return current for the net onshore mass transport in the zone between wave trough and wave crest <ref> Cox, D.T. and Kobayashi, N. 1998. Application of an undertow model to irregular waves on plane and barred beaches. J.Coast.Res. 14: 1314-1324 </ref>. The undertow current is particularly strong just in front of the location where waves are breaking and where much sediment is brought into suspension<ref> Reniers, A.J.H.M., Thornton, E.B., Stanton, T.P. and Roelvink, J.A. 2004. Vertical flow structure during Sandy Duck: observations and modelling. Coast.Eng. 51: 237-260 </ref>. Bound infragravity waves may also contribute to net offshore sediment transport, as discussed in the article [[Infragravity waves]].  
+
Offshore transport is mainly due to so-called [[undertow]]: the compensating return current for the net onshore mass transport in the upper water column between wave trough and wave crest <ref>Stive, M.J.F. and Wind, H.G. 1986. Cross-shore mean flow in the surf zone. Coastal Eng. 10: 325-340</ref><ref> Cox, D.T. and Kobayashi, N. 1998. Application of an undertow model to irregular waves on plane and barred beaches. J.Coast.Res. 14: 1314-1324 </ref><ref> Gallagher, E.L., Elgar, S. and Guza, R.T. 1998. Observations of sand bar evolution on a natural beach. J. Geophys. Res. 103: 3203–3215</ref>. The offshore suspended sediment transport by the undertow current is particularly strong just landward of the location where waves are breaking (plunging waves) and where much sediment is brought into suspension by wave energy dissipation and by the longshore current<ref>Thornton, E.B., Humiston, R.T. and Birkemeier, W. 1996. Bar/trough generation on a natural beach. J. Geophys. Res. 101: 21-97-12110</ref><ref> Reniers, A.J.H.M., Thornton, E.B., Stanton, T.P. and Roelvink, J.A. 2004. Vertical flow structure during Sandy Duck: observations and modelling. Coast.Eng. 51: 237-260 </ref>. Although the undertow velocity is smaller than the amplitude of the wave orbital velocity, it produces substantial offshore transport across the shoreface, in particular in the surf zone, where it often dominates the onshore transport related to wave-asymmetry effects. The undertow current is strongest near the bottom in the surf zone, whereas it peaks near the surface on the lower shoreface. On the lower shoreface it is largely offset by onshore-directed Stokes drift and therefore less effective for offshore transport<ref>Lentz, S.J. and Fewings, M. 2012. The wind- and wave-driven inner-shelf circulation. Annual Review of Marine Science 4: 317–341</ref><ref>Larsen, B.E. and Fuhrman, D.R. 2023. Simulation of cross-shore breaker bar development utilizing a stabilized two-equation turbulence model. Coastal Engineering 180, 104269</ref>.
 +
 
 +
Bound infragravity waves may also contribute to net offshore sediment transport<ref>Roelvink, J.A. and Stive, M.J.F. 1989. Bar generating cross-shore flow mechanisms on a beach. J. Geophys. Res. 94: 4785-4800</ref>, as discussed in the article [[Infragravity waves]]. The transport component related to infragravity waves can be quite substantial on dissipative beaches, but the net contribution is subdued after wave breaking because the phase relationship with incident waves is lost and the net sand transport can be either onshore or offshore directed<ref name=A13>Aagaard, T., Greenwood, B. and Hughes, M. 2013. Sand transport on dissipative, intermediate and reflective beaches. Earth Sci. Rev. 124: 32–50</ref>.  However, small transport contributions cannot be ignored because the large transport components related to wave skewness / asymmetry and undertow largely cancel each other out<ref>Ruessink, B.G., Houwman, K.T. and Hoekstra, P. 1998. The systematic contribution of transporting mechanisms to the cross-shore sediment transport in water depths of 3 to 9 m. Mar. Geol. 152: 295–324</ref>. For this reason, the downslope transport component must also be considered, see next section. 
  
 
===Downslope sediment transport===
 
===Downslope sediment transport===
Downslope transport occurs when sediment particles roll down the shoreface slope. However, this is not the primary cause of downslope sediment transport for fine to medium sandy beaches <ref name=S>Stive, M. J. F. and de Vriend, H. J. 1995. Modelling shoreface profile evolution. Marine Geology 126: 235–248</ref>. Sediment transport by downrolling is a minor effect for two reasons: (1) the shoreface slope is small and (2) sediment particles are mainly transported in suspension, more than as bedload, especially on the upper shoreface. When a sediment particle is lifted from the bed to a certain level <math>z</math> into the fluid, it moves with the wave orbital flow over a net distance <math>l</math> until settling again onto the bed. The net distance <math>l</math> travelled by the sediment particle depends on the settling time; the longer the settling time, the greater the distance. The settling time is inversely proportional to the sediment fall velocity <math>w</math>. Because the shoreface has a seaward dipping slope, the net travelled distance by an offshore moving sediment particle is larger than for a flat bottom; the inverse holds for an onshore moving sediment particle. The slope effect thus requires a correction to the net sediment transport over a flat bottom, which depends on the slope <math>\beta</math> and the fall velocity <math>w</math>. Sediment transport models assume in general that the slope effect can be represented by a correction factor to the wave-induced transport components of the form <ref> Bagnold, R.A. 1963. Mechanics of marine sedimentation. In: The sea, vol.3. Ed. M.N.Hill, Wiley-Interscience: 507-528 </ref>
+
Several mechanisms can induce downslope sand transport. We do not consider here downslope turbidity currents or mud avalanching that may occur at the outer edge of deltas loaded with large amounts of fluvial sediment.
 +
 
 +
One mode of downslope transport consists of sand particles rolling down the shoreface slope. However, modeling suggests that this does not contribute significantly to downslope transport for fine to medium sandy beaches<ref name=S>Stive, M. J. F. and de Vriend, H. J. 1995. Modelling shoreface profile evolution. Marine Geology 126: 235–248</ref>.  
 +
 
 +
A second transport mode is so-called wave-supported gravity flow (WSGF). Observations on the Dutch coast show that wave stirring of the seabed under storm conditions can generate and sustain a high-density layer of suspended fine sand (thickness 15-20 cm, concentration exceeding 40 g/l) flowing down the shoreface slope (1/100 - 1/200) with a velocity of a few cm/s <ref> Flores, R. P., Rijnsburger, S., Meirelles, S., Horner-Devine, A. R., Souza, A. J., Pietrzak, J. D., Henriquez, M. and Reniers, A. 2018. Wave Generation of Gravity-Driven Sediment Flows on a Predominantly Sandy Seabed. Geophysical Research Letters 45: 7634-7645</ref>. Modeling suggests that the high-density sediment layer is kept in suspension by the combined action of wave-induced turbulent stresses and hindered settling. The wave-supported gravity flow was observed during a storm with near-bed wave orbital velocities of about 0.5 m/s, lasted for more than a day and far exceeded in volume other cross-shore sand transport processes. However, the frequency of occurrence of WSGF events is not well known. Wave-supported gravity flows have also been observed on the shoreface of muddy coasts<ref>Peng, Y., Yu, Q., Du, Z., Wang, L, Wang, Y and Gao, S. 2022, Gravity-driven sediment flows on the shallow sea floor of a muddy open coast. Marine Geology 445, 106759</ref>.
 +
 
 +
A third mechanism of downslope transport is related to the wave orbital excursion that moves sediment particles between zones of nearshore high bed shear stress and offshore low bed shear stress. Since resuspension after settling takes more time lower on the shoreface than higher on the shoreface, the downslope distance traveled by settling sediment particles is on average longer than the onshore distance. The average difference between offshore and onshore travel distance depends not only on the orbital excursion and shoreface slope, but also on the particle settling velocity. During the offshore excursion, sediment particles with a lower settling velocity settle later, thus lower, on the shoreface where the bed shear stresses are lower and the resuspension time longer. Fine sediment particles thus experience greater downslope transport than coarse particles. The slope effect requires a correction to the net sediment transport over a flat bottom, which is an increasing function of the product of bed slope <math>\tan \beta</math> and wave orbital excursion, and a decreasing function of the fall velocity <math>w</math>. Empirical formulas for suspended sediment transport assume in general that the slope effect can be represented by a correction factor to the instantaneous wave-induced transport components of the form<ref> Bagnold, R.A. 1963. Mechanics of marine sedimentation. In: The sea, vol.3. Ed. M.N.Hill, Wiley-Interscience: 507-528 </ref>
 +
 
 +
<math>(1- \gamma_s \Large \frac{u}{w}  \normalsize ) \, \tan \beta, \qquad (4)</math>
 +
 
 +
where <math>u</math>  is the instantaneous wave orbital velocity, <math>w</math> the sediment fall velocity and <math>\gamma_s</math> an adjustable parameter depending, among other things, on the resuspension processes. This is consistent with laboratory experiments, in which an approximate correspondence  was found between coastal profiles and an adjusted Dean parameter<ref>Bayle, P.M., Beuzen, T., Blenkinsopp, C.E., Baldock, T.E. and Turner, I.L. 2021. A new approach for scaling beach profile evolution and sediment transport rates in distorted laboratory models. Coastal Engineering 163, 103794</ref>
 +
 
 +
<math>\Omega_{\beta} = \Large \frac{H_s \tan \beta}{wT_p} \normalsize . \qquad(5)</math>
 +
 
 +
Experiments by Hattori and Kawamata (1980<ref>Hattori, M. and Kawamata, R. 1980. Onshore–offshore transport and beach profile change. In: Proc. 17th ASCE Coastal Engineering Conference, 2, pp. 1175–1193</ref>) with different profile slopes show that onshore transport dominates for <math>\Omega_{\beta} < 0.08</math> and offshore transport dominates for <math>\Omega_{\beta} > 0.08</math>. An equilibrium slope corresponds to (order of magnitude)
 +
 
 +
<math>\Omega_{\beta} \sim 0.08. \qquad(6)</math>
  
<math>(1-\gamma \Large \frac{\beta u }{w}  \normalsize ) , \qquad(2)</math>
+
Comparison with Eq. (3) shows that reflective beaches have a steep slope and dissipative beaches a gentle slope. The presence of the fall velocity <math>w</math> in the denominator of the downslope correction term in Eq. (4) also implies that net wave-induced transport can be downslope (offshore) for fine sediment even if upslope transport (onshore) prevails for coarse sediment. 
  
where <math>\beta</math> is the bed slope, <math>u</math>  the wave orbital velocity, <math>w</math> the sediment fall velocity and <math>\gamma</math> an efficiency parameter.
+
Although not proven, it is often assumed that the slope effect is a major factor compensating for net wave-induced onshore transport (sum of all wave-induced onshore and offshore transport components including undertow) in situations close to equilibrium<ref name=B>Bowen, A.J. 1980. Simple models of nearshore sedimentation; beach profiles and longshore bars. In: The coastline of Canada, Ed. by S.B. McCann, Geological Survey of Canada, Ottawa, pp. 1-11</ref>. This implies that the equilibrium shoreface slope <math>\tan \beta </math>  will be steeper (larger <math>\tan \beta </math>) for a coarser seabed (large fall velocity <math>w</math>) than for a finer sediment bed (small fall velocity <math>w</math>, smaller <math>\tan \beta </math>). This is expressed by the correction factor (Eq. 4) and is in general agreement with field observations. An analysis of long-term shoreface profile observations at the Gold Coast (Australia) suggests that the bed-slope effect on the total net sediment transport averaged over the wave period (or, more generally, the wave climate) can be related to the equilibrium bed slope <math>\tan \beta_{eq}</math> by the correction factor <math>(1 – \tan \beta / \tan \beta_{eq})</math> <ref>Patterson, D.C. and Nielsen, P. 2016. Depth, bed slope and wave climate dependence of long term average sand transport across the lower shoreface. Coastal Engineering 117: 113–125</ref>, which is equivalent to Eq. (4).
Although not proven, it is often assumed that the slope effect is the primary factor compensating for net wave-induced onshore transport in situations close to equilibrium. This implies that the equilibrium shoreface slope <math>\beta</math>  will be steeper (larger <math>\beta</math>) for a coarser seabed (large fall velocity <math>w</math>) than for a finer sediment bed (small fall velocity <math>w</math>, smaller <math>\beta</math>), according to the expression (2) for the correction factor. This is in general agreement with field observations.
 
  
===Sediment transport modelling===
+
===Shoreface modelling===
The processes involved in sediment transport on the shoreface, as described above, are complicated and not yet fully understood. Even the most sophisticated process-based mathematical models, like XBEACH <ref>https://oss.deltares.nl/web/xbeach/home</ref>, make use of approximate empirical formulations for several processes contributing to sediment transport and involve parameters that need to be tuned. Empirical formulas for wave-induced sediment transport on the shoreface that are often used in practice are given in the articles [[Sediment transport formulas for the coastal environment]] and [[Sand transport]]. These formulas are based on laboratory and field measurements for a large range of conditions. Rough estimates can be obtained with the empirical default parameters given in these articles. For better estimates it is necessary to tune these parameters to field data.  
+
The processes involved in sediment transport on the shoreface, as described above, are complicated and not yet fully understood. Even the most sophisticated numerical morphodynamic models, like Delft3D and XBEACH <ref>https://oss.deltares.nl/web/xbeach/home</ref>, make use of approximate empirical formulations for several processes contributing to sediment transport and involve parameters that need to be tuned. Empirical formulas for wave-induced sediment transport on the shoreface that are often used in practice are given in the articles [[Sediment transport formulas for the coastal environment]] and [[Sand transport]]. These formulas are based on laboratory and field measurements for a large range of conditions. Rough estimates can be obtained with the empirical default parameters given in these articles. For more accurate site-specific estimates it is necessary to tune these parameters to field data. A more detailed discussion of shoreface modelling is presented in the article [[Modelling coastal hydrodynamics]].
  
  
==Equilibrium shoreface profile==
+
==Equilibrium shoreface profile (sandy beach)==
  
For a given stable wave climate the shoreface profile tends after a long period to an equilibrium profile. This equilibrium profile can be computed by applying sediment transport formulas to an initial profile for a period comprising all wave conditions according to local wave climate statistics. As this requires tedious computations and produces results with large uncertainty margins, equilibrium profiles derived from field observations are more often used. From analysis of a large number of beach profiles of the Californian and Danish coasts <ref> Bruun, P. 1954. Coast erosion and the development of beach profiles. Beach Erosion Board, US Army Corps of Eng., Tech.Mem. 44: 1-79 </ref><ref> Dean, R.G. 1977. Equilibrium beach profiles: US Atlantic coast and Gulf coasts. Ocean Eng.Tech.Rep. 12, Univ. of Delaware, Newark, 45 pp .</ref> and the US Atlantic and Gulf coasts <ref  name=D>Dean, R.G. 1977. Equilibrium beach profiles: US Atlantic coast and Gulf coasts. Ocean Eng.Tech.Rep. 12, Univ. of Delaware, Newark, 45 pp.</ref> Bruun (1954) and Dean (1977) derived an equilibrium profile of the form (Fig. 3)
+
[[Image: DeanAndBodgeProfiles.jpg|thumb|300px|right|Fig. 6. The Bruun/Dean profile (Eq. 7) and the equivalent Bodge profile (Eq. 8).]]
  
[[Image: DeanAndBodgeProfiles.jpg|thumb|300px|right|Figure 3. The Bruun/Dean profile (Eq. 3) and the equivalent Bodge profile (Eq. 4).]]
+
For a given stable wave climate the shoreface profile tends after a certain (long) period to an equilibrium profile, as pointed out before. This equilibrium profile can be computed by applying sediment transport formulas to an initial profile for a period comprising all wave conditions according to local wave climate statistics. This requires heavy calculations that are moreover subject to large uncertainty margins. Therefore, equilibrium profiles derived from field observations are often used. By analyzing a large number of beach profiles of the Californian and Danish coasts <ref> Bruun, P. 1954. Coast erosion and the development of beach profiles. Beach Erosion Board, US Army Corps of Eng., Tech.Mem. 44: 1-79 </ref> and the US Atlantic and Gulf coasts <ref  name=D>Dean, R.G. 1977. Equilibrium beach profiles: US Atlantic coast and Gulf coasts. Ocean Eng.Tech.Rep. 12, Univ. of Delaware, Newark, 45 pp.</ref> Bruun (1954) and Dean (1977) derived an equilibrium profile (often called the Dean/Bruun profile) of the form (Fig. 6)
  
<math>h(x)=A \; x^{2/3}, \qquad(3)</math>
+
<math>h(x)=A \; x^{2/3}, \qquad(7)</math>
  
where <math>h</math> is the water depth [m], <math>x</math> the cross-shore distance [m] (<math>x=0</math> at the shoreline) and <math>A </math> a coefficient depending on the sediment fall velocity <math>w</math> [m/s], <math> \quad A \approx 0.5 \; w^{0.44} [m^{1/3}]</math>.  
+
where <math>h</math> is the water depth [m], <math>x</math> the cross-shore distance [m] (<math>x=0</math> at the shoreline) and <math>A</math> a coefficient depending on the sediment fall velocity <math>w</math> [m/s]. Dean (1987<ref>Dean, R. 1987. Coastal sediment processes: Toward coastal engineering solutions. Spec. Conf. Understanding Coastal Sediment Processes. N.C. Kraus, ed. ASCE, New York, pp. 1-24</ref>) empirically determined <math>A \approx 0.5 \; w^{0.44} </math> [m<sup>1/3</sup>]. Another estimate was given by Kriebel et al. (1991<ref>Kriebel, D.L., Kraus, N.C. and Larson, M. 1991. Engineering Methods for Predicting Beach Profile Response. Procs. Conf. Quantitative Approaches to Coastal Sediment Processes. N.C. Kraus, K.J. Gingerich, D.L. Kriebel eds. ASCE, New York, pp. 557-571</ref>) <math>A \approx 2.25 \; (w^2/g)^{1/3}</math> [m<sup>1/3</sup>]. In Eq. (7) it is implicitly assumed that the fall velocity is uniform across the shoreface profile.  
  
 
Dean <ref name=D></ref> showed that the exponent 2/3 is consistent with the assumption of constant wave energy dissipation per unit volume throughout the surf zone.
 
Dean <ref name=D></ref> showed that the exponent 2/3 is consistent with the assumption of constant wave energy dissipation per unit volume throughout the surf zone.
Line 85: Line 158:
 
Because the Bruun/Dean profile has an infinite slope at the shoreline, an alternative form was proposed by Bodge (1992) <ref>Bodge, K.R. 1992. Representing equilibrium beach profiles with an exponential expression. J. Coastal Research 8: 47-55.</ref>,
 
Because the Bruun/Dean profile has an infinite slope at the shoreline, an alternative form was proposed by Bodge (1992) <ref>Bodge, K.R. 1992. Representing equilibrium beach profiles with an exponential expression. J. Coastal Research 8: 47-55.</ref>,
  
<math>h(x)=B \; (1 – e^{-kx}) , \qquad(4)</math>
+
<math>h(x)=B \; (1 – e^{-kx}) , \qquad(8)</math>
 +
 
 +
where <math>k</math> is an empirical coefficient in the range <math>500 < k^{-1} < 2000 </math> [m] and where the coefficient <math>B</math> is related to the beach slope <math>\tan \beta=Bk</math> at the shoreline <math>x=0</math>.
  
where <math>k</math> is an empirical coefficient in the range <math>3 \; 10^{-5} - 1.16 \; 10^{-3} \; [m^{-1}]</math>.
+
[[File:UpperLowerShoreface.jpg|thumb|right|300px|Fig. 7. Theoretical equilibrium profile for the upper shoreface (Eq. 7) and lower shoreface (Eq. 9).]]
  
From an analysis of terrace-shaped shoreface profiles in northern Spain, Bernabeu et al. (2003) <ref>Bernabeu, A.M., Medina, R. and Vidal, C. 2003. A morphological model of the beach profile integrating wave and tidal influences. Mar.Geol. 197: 95-116</ref> found that the Bruun/Dean profile was adequate only for the upper shoreface; for the lower shoreface they proposed the form
+
The Bruun/Dean profile provides in general a reasonable description of observed long-term averaged profiles of micro-mesotidal sandy coasts<ref>Bernabeu, A.M., Medina, R. and Vidal, C. 2003. A morphological model of the beach profile integrating wave and tidal influences. Mar.Geol. 197: 95-116</ref>. This holds in particular for the upper shoreface where breaking waves are the dominant driving agent for sediment transport. The 2/3 power law is generally less adequate for describing the profile of the lower shoreface. This is not surprising, as sediment transport here is mainly driven by non-breaking waves. In a simplified ideal situation, based on realistic assumptions, a theoretical profile was derived through application of the principle of minimum energy dissipation<ref>Larson, M., Kraus, N. C. and Wise, R. A. 1999. Equilibrium beach profiles under breaking and non-breaking waves. Coastal Engineering 36: 59–85</ref><ref name=M20>Maldonado, S. 2020. Do beach profiles under nonbreaking waves minimize energy dissipation? Journal of Geophysical Research: Oceans 125, e2019JC015876</ref>. The theoretical profile also follows a power law (see Fig. 7):
  
<math>h(x)=A' \; (x-x_0)^{2/3} , \qquad(5)</math>
+
<math>h(x) = \Big[ h_b^{1/n} + (h_c^{1/n} - h_b^{1/n}) \Large\frac{x - x_b}{x_c - x_b}\normalsize \Big]^n , \quad x_b < x < x_c , \qquad (9)</math>
  
where <math>x_0</math> is determined such that the two profiles (3) and (5) match at the inflexion point between the upper and lower shoreface.
+
where <math>h_b=h(x_b)</math> is the depth at the transition from upper to lower shoreface and <math>x_b</math> the corresponding location, and <math>h_c=h(x_c)</math> is the closure depth and <math>x_c</math> the corresponding location. The exponent <math>n</math> has a value in the range 0.3 - 0.4, depending on the formulation of the bottom friction<ref name=M20/>.
 +
Inman et al. (1993) <ref>Inman, D.L., Elwany, M.H. and Jenkins, S.A. 1993. Shorerise and bar-berm profiles on ocean beaches. J.Geophy.Res. 98: 18181-18199</ref> found best fits with exponents ranging between 0.3 and 0.45 from an analysis of beach profiles for the coasts of south California, North Carolina and the Nile delta.  
  
Inman et al. (1993) <ref>Inman, D.L., Elwany, M.H. and Jenkins, S.A. 1993. Shorerise and bar-berm profiles on ocean beaches. J.Geophy.Res. 98: 18181-18199</ref> used a similar procedure for the analysis of beach profiles for the coasts of south California, North Carolina and the Nile delta. They found best fits with exponents ranging between 0.3 and 0.45 instead of 2/3.
+
For gravel beaches other cross-shore profiles are found, see [[Gravel Beaches]].
  
  
 
==Analytical equilibrium models==
 
==Analytical equilibrium models==
  
Because of the earlier mentioned difficulties to derive model estimates of the equilibrium shoreface profile taking into account the full wave climate, analytical models have been developed in which the wave climate is replaced by a single representative monochromatic incident wave. An alongshore uniform shoreface is assumed with perpendicular wave incidence. A simple wave transformation model is used based on depth saturation of broken waves and [[shallow-water wave theory]]. An analytical expression for the equilibrium shoreface slope can be derived by requiring that the total wave-integrated cross-shore sediment flux equals zero throughout the shoreface zone <ref>Bowen, A.J. 1980. Simple models of nearshore sedimentation; beach profiles and longshore bars. In: The coastline of Canada, Ed. by S.B. McCann, Geological Survey of Canada, Ottawa, pp. 1-11</ref> <ref name=S></ref> <ref name =O> Ortiz, A.C. and Ashton, A.D. 2016. Exploring shoreface dynamics and a mechanistic explanation for a morphodynamic depth of closure. J. Geophys. Res. Earth Surf 121: 442–464 doi:10.1002/2015JF003699</ref>. Although these models reproduce certain properties also observed in the field, such as increasing profile steepness with increasing grain size and with increasing wave period, a more detailed comparison of computed profiles with observations exhibits large discrepancies. The computed concavity of the upper shoreface profile is much stronger than generally observed in the field and much stronger than for the Bruun/Dean profile. This is perhaps not surprising in view of the simplicity of the model for wave-induced sediment transport, which ignores many of the processes described in the previous paragraphs. The model also ignores the influence of major storms, which flatten the profile out to great depths. The profile adaptation time scale for the entire shoreface being quite long (order of centuries <ref name =S></ref>), it is plausible that observed shoreface profiles often correspond to transient profiles <ref name = O></ref>.   
+
Because of the earlier mentioned difficulties to derive model estimates of the equilibrium shoreface profile, taking wave climate statistics fully into account, analytical models have been developed in which the wave climate is replaced by a single representative monochromatic incident wave. In these models an alongshore uniform shoreface is assumed with perpendicular wave incidence, and a simple wave transformation model is used based on [[Breaker index|depth saturation of broken waves]] (wave height limited by depth) and [[shallow-water wave theory]]. An analytical expression for the equilibrium shoreface slope can be derived by requiring that the total wave-integrated cross-shore sediment flux equals zero throughout the shoreface zone <ref name=B>Bowen, A.J. 1980. Simple models of nearshore sedimentation; beach profiles and longshore bars. In: The coastline of Canada, Ed. by S.B. McCann, Geological Survey of Canada, Ottawa, pp. 1-11</ref> <ref name=S></ref><ref name =O>Ortiz, A.C. and Ashton, A.D. 2016. Exploring shoreface dynamics and a mechanistic explanation for a morphodynamic depth of closure. J. Geophys. Res. Earth Surf 121: 442–464 doi:10.1002/2015JF003699</ref>. Although these models reproduce certain properties also observed in the field, such as increasing profile steepness with increasing grain size and with increasing wave period, a more detailed comparison of computed profiles with observations exhibits large discrepancies. This is perhaps not surprising in view of the simplicity of the model used for wave-induced sediment transport, which assumes zero phase difference between wave-orbital velocity and sediment load, and which ignores many other processes described in the previous paragraphs. It has been shown<ref>Aagaard, T. and Huges, M.G. 2017. Equilibrium shoreface profiles: A sediment transport approach. Marine Geology 390: 321–330</ref> that some discrepancies can be repaired by using expressions for the sediment transport components derived from field observations.
 +
 
 +
Analytical models also ignore the influence of extreme storms, which can flatten the profile out to great depths within a short period. The profile adaptation time scale for the entire shoreface being quite long<ref name =S/>, it is plausible that observed shoreface profiles generally correspond to transient profiles <ref name = O></ref>.   
  
 
   
 
   
  
 
==Related articles==
 
==Related articles==
 +
:[[Nearshore sandbars]]
 +
:[[Bruun rule]]
 +
:[[Active coastal zone]]
 +
:[[Shoreline retreat and recovery]]
 +
:[[Coastal Hydrodynamics And Transport Processes]]
 +
:[[Modelling coastal hydrodynamics]]
 +
:[[Coastal and marine sediments]]
 +
:[[Gravel Beaches]]
 +
 +
 +
==Further reading==
 +
:Komar, P.D. 1998. Beach processes and sedimentation. Prentice Hall, London, pp. 544.
 +
:Woodroffe, C.D. 2002. Coasts, form, processes and evolution. Cambridge Univ.Press, 623 pp.
 +
:Dronkers, J. 2017. Dynamics of Coastal Systems. World Scientific Publ. Co. 753 pp.
 +
 +
 +
 +
==Table Symbols==
 +
{| class="wikitable"
 +
|-
 +
! Variable !! Description !! Variable !! Description
 +
|-
 +
| <math>c</math>|| wave propagation speed (<math>[m/s]</math>) ||<math>c_g</math> || wave group propagation speed (<math>[m/s]</math>)
 +
|-
 +
| <math>d_{50}</math> || mean sediment grainsize (<math>[m]</math>) ||<math>g</math>|| gravitational acceleration (<math>\approx 9.8 \; m/s^2</math> )
 +
|-
 +
| <math>h</math> ||  water depth (<math>[m]</math>) ||<math>h_b</math> ||  water depth at breaking location (<math>[m]</math>)
 +
|-
 +
| <math>H_b</math>|| significant wave height at breaking location (<math>[m]</math>)||<math>H_s</math>|| significant wave height (<math>[m]</math>)
 +
|-
 +
| <math>k</math>|| wave number (<math>[m^{-1}]</math>) ||<math>t</math>|| time (<math>[s]</math>)
 +
|-
 +
| <math>T_p</math>|| spectral peak wave period (<math>[s]</math>) ||<math>u</math> || velocity (<math>[m/s]</math>)
 +
|-
 +
| <math>w</math>|| sediment fall velocity  (<math>[m/s]</math>) ||<math>x</math>|| distance from shoreline (seaward >0, <math>[m]</math>)
 +
|-
 +
| <math>\alpha</math>|| wave incidence angle ||<math>\tan \beta</math>|| shoreface slope
 +
|-
 +
| <math>\theta</math>|| wave phase <math>\theta = \omega t - kx</math> ||<math>\rho</math>|| seawater density (<math>[kg/m^3]</math>)
 +
|-
 +
| <math>\omega</math>|| radial wave frequency  (<math>[s^{-1}]</math>) ||<math>\Omega</math>|| Dean parameter (<math>= \Large \frac{H_b}{wT_p} \normalsize </math>)
 +
|}
 +
 +
 +
==Appendix: Wave skewness==
 +
 +
When a sinusoidal wave travels into shallow water the wave shape changes. The equations describing wave motion in deep water are almost linear; the solution is a simple sinus-shaped wave. The nonlinear terms in the wave equations become important in shallow water. These nonlinear terms introduce higher harmonic wave components in the solution. Initially, the first harmonic components grow the fastest. The wave orbital velocity can then be written as
 +
 +
<math>u(\theta) = U_1 \cos(\theta) + U_2 \cos(2\theta) + V_2 \sin(2\theta) , \qquad (A1)</math>
 +
 +
where <math>\theta = \omega t - kx</math> is the wave phase, with <math>t=</math>time, <math>x=</math>the distance travelled in the wave propagation direction, <math>\omega =</math>the wave angular frequency, <math>k=\omega / c=</math>the wave number and <math>c=</math>the wave propagation speed. The <math>U_2</math> component represents the velocity skewness <math>S_u</math> and the component <math>V_2</math> represents the wave asymmetry or acceleration skewness <math>A_u</math>  (Fig. A1).
 +
 +
 +
[[File:WaveAsymmetries.jpg|thumb|center|600px|Fig. A1. Wave orbital velocity as function of the phase <math>\theta</math>. Left panel: Velocity skewness. Right panel: Acceleration skewness. Blue curves: Eq. (A1). Red curves: Eq. (A5). Positive acceleration skewness stimulates onshore bedload transport, because the fast increase of the onshore orbital motion leaves the boundary layer barely time to grow, leading to strong near-bed velocity gradients and shear stresses<ref>Van der A, D.A., O’Donoghue, T. and Ribberink, J.S. 2010. Measurements of sheet flow transport in acceleration-skewed oscillatory flow and comparison with practical formulations. Coastal Engineering 57: 331–342</ref>. It has been suggested that high pressure forces under the steep wave front also play a role by mobilizing bed grains around flow reversal<ref>Calantoni, J. and Puleo, J.A. 2006. Role of pressure gradients in sheet flow of coarse sediments under sawtooth waves. Journal of Geophysical Research 111, 2005JC002875</ref>.]]
 +
  
[[Bruun rule]]
+
We define <math> <…> \equiv \Large\frac{1}{2 \pi}\normalsize \int_0^{2 \pi} … d \theta </math> and assume <math>|U_2|,|V_2| << |U_1| .</math>
  
[[Active coastal zone]]
+
The velocity and acceleration skewness are defined as
  
[[Definitions of coastal terms]]
+
<math>S_u \equiv \Large\frac{<u^3>}{<u^2>^{3/2}}\normalsize  \approx \Large\frac{3 U_2}{\sqrt{2} U_1}\normalsize, \quad  A_u \equiv \Large\frac{ < (\mathcal{H}u)^3 >}{<(\mathcal{H}u)^2>^{3/2} }\normalsize \approx \Large\frac{3 V_2}{\sqrt(2) U_1}\normalsize . \qquad (A2)</math>
  
[[Coastal and marine sediments]]
+
(the Hilbert transform  <math>\mathcal{H}</math> applied to sinusoidal functions gives: <math>\mathcal{H}cos(x)=-sin(x), \; \mathcal{H}sin(x)=cos(x), \; \mathcal{H}cos(2x)=-sin(2x), \; \mathcal{H}sin(2x)=cos(2x)</math>,  etc).
  
 +
By analysing a large number of field observations of the wave orbital velocity on the shoreface, Ruessink et al. (2012<ref>Ruessink, B.G., Ramaekers, G. and van Rijn, L.C. 2012. On the parameterization of the free-stream non-linear wave orbital motion in nearshore morphodynamic models. Coastal Engineering 65: 56–63</ref>) found empirical evidence that the velocity and acceleration skewness mainly depend on the Ursell number
  
+
<math>U_r = \Large\frac{3 g H}{8 h^2 \omega^2}\normalsize , \qquad (A3)</math>
==Further reading==
+
 
 +
where <math>g \approx </math>9.8 m s<sup>-2</sup> is the gravitational acceleration, <math>H</math> the significant wave height and <math>h</math> the local depth. They established the empirical expressions
 +
 
 +
<math> (S_u^2+A_u^2)^{-1/2} = 1.17 \, (1 + 0.205 \, U_r^{-3.37}), \quad \Large\frac{A_u}{S_u}\normalsize = -\cot \Big(\Large\frac{\pi}{2}\normalsize \tanh \big(0.815 U_r^{-0.672}\big)\Big) \normalsize, \qquad (A4)</math>
 +
 
 +
from which <math>S_u</math> and <math>A_u</math> can be solved. Besides the Ursell number, also other factors influencing the velocity and acceleration skewness have been identified, such as the offshore wave steepness, spectral band width and beach slope<ref>Rocha, M.V.L., Michallet, H. and Silva, P.A. 2017. Improving the parameterization of wave nonlinearities – The importance of wave steepness, spectral bandwidth and beach slope. Coastal Engineering 121: 77–89</ref>.
 +
 
 +
Abreu et al. (2010<ref name=A10/>) derived a more general expression than Eq. (A1) for the wave orbital velocity, that includes velocity skewness and acceleration skewness,
 +
 
 +
<math>u(\theta)=U_1 \Large\frac{\sqrt{1-r^2}}{1+\sqrt{1-r^2}} \, \frac{r \sin{\phi}+(1+\sqrt{1-r^2}) \sin(\theta)}{1 - r \cos(\theta +\phi)}\normalsize . \qquad (A5)</math>
  
Komar, P.D. 1998. Beach processes and sedimentation. Prentice Hall, London, pp. 544.
+
The parameters <math>r</math> and <math>\phi</math> are related to the velocity skewness and acceleration skewness by
  
Woodroffe, C.D. 2002. Coasts, form, processes and evolution. Cambridge Univ.Press, 623 pp.
+
<math>r = \Large\frac{2 \sqrt{p+1}}{p+2}\normalsize, \quad p = \Large\frac{9}{2 (S_u^2+A_u^2)}\normalsize , \quad \phi = -\Large\frac{\pi}{2}\normalsize -\arctan(A_u/S_u) . \qquad (A6)</math>
  
Dronkers, J. 2017. Dynamics of Coastal Systems. World Scientific Publ. Co. 753 pp.
+
A predictive model for wave-induced cross-shore sediment transport can be obtained by using the expression (A5) for the wave orbital velocity in formulas for cross-shore sediment transport, for example the formulas presented in the article [[Sediment transport formulas for the coastal environment]]. This yields only the transport contribution of the oscillatory wave component. Other components of cross-shore sediment transport discussed above have to be considered as well.   
   
 
  
  
Line 135: Line 277:
 
|AuthorName=Dronkers J}}
 
|AuthorName=Dronkers J}}
  
[[Category:Land and ocean interactions‏‎]]
+
[[Category:Physical coastal and marine processes]]
[[Category:Geomorphological processes and natural coastal features‏‎]]
+
[[Category:Beaches]]
[[Category:Coastal and marine natural environment]]
+
[[Category:Morphodynamics]]

Revision as of 10:16, 19 April 2024

Definition of Shoreface:
The shoreface is the nearshore zone of the inner continental shelf that is bounded landward by the mean low-water line and that extends seaward to where the influence of wave action on cross-shore sediment transport is on average minor compared to other influences. A similar definition is: The shoreface is the zone seaward of the shoreline where offshore generated waves interact with the upward sloping seabed.
This is the common definition for Shoreface, other definitions can be discussed in the article
Definition of Shoreface profile:
The shoreface profile, often called beach profile, is the cross-shore coastal depth profile of the shoreface.
This is the common definition for Shoreface profile, other definitions can be discussed in the article


This article focuses on the shoreface of sandy coasts. The shoreface of rocky shores is described in the article Rocky shore morphology and the shoreface of mud coasts in the article Coastal mud belt.



Closure depth

Fig. 1. Schematic representation of a shoreface profile. The vertical/horizontal scale ratio is greatly exaggerated.

In the literature one can find various definitions of the shoreface. The definition given here, which is most commonly used, follows the definition given in Stive and De Vriend (1995[1]). The shoreface can be divided in two zones, the upper shoreface and the lower shoreface. The upper shoreface is the zone where under average conditions waves break and the lower shoreface the zone where waves shoal. The lower part of the shoreface extends to the so-called outer closure depth [math]h_{out}[/math]; beyond this depth the seabed is hardly influenced by waves and wave-induced sediment transport is (on average) insignificant. The lower part of the shoreface is also called shoaling zone; wave shoaling is the process of wave amplification when waves travel from deep to shallow water before breaking (see Shallow-water wave theory). The outer closure depth [math]h_{out}[/math] has been empirically determined by Hallermeier (1981[2]), yielding the formula (see also Closure depth)

[math]h_{out} \approx (H_s – 0.3 \sigma) \, T_s \sqrt{\large\frac{g}{5000 d_{50}}\normalsize} \approx 0.013 H_s T_s \sqrt{\large\frac{g}{d_{50}}\normalsize} , \qquad (1)[/math]

where [math]H_s, T_s[/math] are the average significant wave height and period, [math]\sigma[/math] is the wave height standard deviation, [math]d_{50}[/math] the medium sediment grainsize and [math]g[/math] the gravitational acceleration. The second approximation in Eq. (1) follows from the assumption of a Rayleigh wave distribution (see Statistical description of wave parameters). Typical values of the outer closure depth are in the range 15 – 40 m.

The definition of the outer closure depth implies that wave-driven sediment transport across this boundary is very small on average. It does not rule out significant sediment transport across this boundary due to e.g. tidal currents, wind-driven currents, ocean currents or internal waves. For example, the occurrence of shoreface-connected ridges (see Sand ridges in shelf seas) on the lower shoreface can be theoretically related to sediment transport by strong wind-driven currents. Other definitions of the outer boundary of the lower shoreface can also be found in the literature, for example the so-called Depth of Transport. This is the depth at which the total bed shear stress does not exceed the threshold of sediment transport under extreme storms[3]. The Depth of Transport is typically further seaward than the outer closure depth Eq. (1).

Most energy is dissipated in the upper part of the shoreface, where wave overturning and breaking occurs. The depth where waves start breaking depends on the wave climate and is therefore highly variable. An often used definition of the transition depth [math]h_{in}[/math] between upper and lower shoreface was also given by Hallermeier (1978[4]),

[math]h_{in} \approx 2.28 H_{12h/y}- 68.5 \Large\frac{ H_{12h/y}^2}{g T_{12h/y}^2}\normalsize , \qquad (2)[/math]

where [math]H_{12h/y}[/math] is the effective wave height just seaward of the breaker zone that is exceeded on average for 12 hours per year (i.e. the significant wave height with a probability of yearly exceedance of 0.137%) and [math]T_{12h/y}[/math] the associated wave period. Later refinements of this formula are discussed in the article Closure depth. Typical values of [math]h_{in}[/math] are in the range 3 – 10 m. The width of the upper shoreface typically ranges between about 100 - 1000 meters and the average slope between 1/10 and 1/100. The upper shoreface is sometimes also called breaker zone or surf zone, although onshore surfing wave bores usually only occur in the shoreward part of the upper shoreface. It is therefore preferable to use the designation upper shoreface as a fixed geographical indication, and the designations surf zone or breaker zone as a physical indication of the zone where broken waves surf onshore. This latter zone is narrow and close to the shoreline in a gentle wave climate and can be very wide under storm conditions, extending from the seaward boundary of the upper shoreface to the dune foot. Other subdivisions of the coastal profile include the active coastal zone, the beach face or foreshore and the backshore, see Definitions of coastal terms.

The division between upper and lower shoreface profiles is often marked by an inflexion point, a bar or a terrace, indicating transition of morphodynamics [5][6][7] (Fig. 1). Another indicator of the transition between upper and lower shoreface is a break in the trend of the sediment grainsize with depth[8]. The shoreface is bordered on the landward side by the beach face. This is the zone where under average conditions waves collapse and run up the beach (see Swash zone dynamics). The upper shoreface profile is often undulated due to the presence of one or several bars, see Nearshore sandbars. These bars have mainly a shore parallel orientation, but more complex three-dimensional patterns are ubiquitous (see Rhythmic shoreline features and Fig. 3).


Morphodynamic feedback

Sediment transport on the shoreface depends on the local wave climate. By 'wave climate' is meant the description of incoming waves in terms of a statistical distribution of incident wave directions, wave heights and wave periods (see Statistical description of wave parameters). Sediment transport also depends on the characteristics of bed sediments and on the seabed bathymetry. For small wave incidence angles (close to the shore-normal direction) the relevant seabed bathymetry is the cross-shore shoreface profile.

Fig. 2. Gold Coast, Australia, before and after the June 1967 cyclone. Observed shoreface profiles showing recovery of the initial close-to-equilibrium profile after strong perturbation by an extreme event. Redrawn after Patterson and Nielsen (2016[9]).

The shoreface profile influences cross-shore sediment transport, but the inverse is also true: changes of the shoreface profile result from gradients in the cross-shore sediment transport. This is called morphodynamic feedback: for a given incident wave field the shoreface profile determines the transformation of the wave field on the shoreface and therefore the resulting sediment transport, whilst the resulting sediment transport determines the shoreface profile. As a result of this morphodynamic feedback, the shoreface profile tends to an equilibrium shape for a given stable wave climate[10][11][12]. An illustration is presented in Fig. 2, which shows the recovery after strong perturbation of a shoreface profile that is initially close to equilibrium.

A dynamic equilibrium is reached when onshore directed sediment transport and offshore directed transport are equal on average over a period exceeding the morphodynamic adaptation time. This holds for alongshore uniform coasts; the equilibrium requirement for non-uniform coasts is a vanishing divergence of the long-term average total sediment transport vector (cross-shore and alongshore). It was realised already long ago that the shoreface of a sedimentary coast should not be regarded as a geologically inherited feature, but as the result of natural adaptation to modern hydrodynamic conditions and sediment supply[13][14]. In real field situations, a coast never assumes its equilibrium profile; it is always in a transient state due to fluctuations in the wave climate and the mean water level. The coastal profile deviates from the equilibrium profile, even in the hypothetical case of a stationary wave climate, due to the cyclical or quasi-cyclical migration of seabed structures along the profile (see Nearshore sandbars).

Knowing how to determine the equilibrium shoreface profile has great practical interest. It provides quantitative insight in the response of the shoreface to changes in the local wave climate by human and natural causes; such insight is important for testing and improving the effectiveness of measures to combat coastal erosion and shoreline retreat.


Bedforms and sediment sorting

Morphodynamic feedback involves not only shoreface profile development, but also sediment grainsize and bed roughness. Grainsize and bed roughness strongly influence sediment transport [15]. Inversely, sediment transport strongly influences the grainsize distribution of bed sediments. Bed roughness depends on the presence of bedforms (ripples, dunes), which in turn depend on sediment transport and sediment grainsize. Consequently, shoreface profile development, wave transformation on the shoreface, sediment transport, grainsize sorting and bedform emergence are all dynamically coupled processes.

Grainsize sorting is primarily due to differences in transport modes of different grainsize fractions. The finer fractions are more easily brought and maintained in suspension than the coarser fractions. They are transported further away from their initial location and tend to settle on the lower shoreface zone where wave-induced bed shear stresses are lower than on the upper shoreface. The coarse fraction is transported partly as suspended load and partly as bedload in the wave boundary layer. Wave transformation in shallow water leading to wave skewness (stronger onshore than offshore wave orbital velocity) and wave asymmetry (stronger onshore than offshore wave orbital acceleration) favour net onshore transport in the wave boundary layer. Fine sediments are resuspended from the upper shoreface due to undertow currents, longshore currents and turbulent shear stresses generated by breaking waves. Sediments on the upper shoreface are therefore coarser than on the lower shoreface. Fine sediment particles can still be present on the upper shoreface if they are 'hidden' by overlying coarser sediment particles – a phenomenon called 'bed armoring' (see Seabed armoring).

Fig. 3. The barred coast of North Holland, The Netherlands. https://beeldbank.rws.nl, Rijkswaterstaat / Rens Jacobs

The largest bedforms on the shoreface, such as shoreface-connected ridges (scale of several kilometers) and shore-oblique bars (scale of the order of hundred meters) are mainly related to morfodynamic feedback processes that interact with longshore transport, see the articles Sand ridges in shelf seas and Rhythmic shoreline features. Other large bedforms, the so-called nearshore sandbars or breaker bars (scale also of the order of 100 meters, Fig. 3) arise from morphodynamic coupling with the wave breaking process, through the resulting flow structure at the breaker bar and interaction with cross-shore and longshore currents [16][17] - see Nearshore sandbars for a more detailed discussion. Bedforms at much smaller scales – ripples with wavelengths of the order of 0.3 - 1 m – are formed in the less exposed parts of the shoreface under low to moderate wave conditions. They are often found on the lower shoreface and in the troughs between nearshore breaker bars. Morphodynamic feedback processes, leading to regular ripple patterns, involve ripple-induced near-bed flow circulation[18], vortex shedding at ripple crests[19] and non-linear ripple-ripple interactions[20]. For a more detailed discussion, see the articles Wave ripples and Wave ripple formation. The time lag of sediment suspension and settling related to vortex shedding at the ripple crests strongly influences the direction of residual cross-shore sediment transport on the shoreface, see section #Onshore wave-induced sediment transport.


Beach classification; dissipative and reflective

The beach profile in the upper shoreface zone can be highly variable. The profile tends to steepen in periods of low-energy waves (especially long-period swell waves) and to flatten under high-energy storm waves (see Active coastal zone). Sediment characteristics, which largely depend on the supply source (see Coastal and marine sediments), also play an important role. Coarse-grained beaches typically have steeper slopes than fine-grained beaches[21]. Gently sloping beaches dissipate almost completely the energy of incident waves, whereas steep-sloping beaches tend to reflect (partly) incident waves. For distinguishing between reflective and dissipative beaches the Dean-parameter [math]\Omega[/math] is often used [22][23]:

[math]\Omega = \Large \frac{H_b}{wT_p} \normalsize , \qquad(3)[/math]

where [math]H_b[/math] is the significant wave height at breaking, [math]T_p[/math] is the peak spectral wave period and [math]w[/math] is the mean sediment fall velocity.

Fig. 4. Beach morphologies according to the parameters [math]\Omega[/math] and [math]RTR[/math], based on beach surveys in England, Wales, Northern Ireland, Australia and China[24][25][26][27]. Reflective beaches are found in the red zone and dissipative beaches in the blue zone.

Wright and Short (1984) [28] analysed beach profiles in Australia by evaluating the parameter [math]\Omega[/math], considering the local long-term average wave climate. They found that beach types could be classified as follows :

  • reflective beaches correspond to [math]\Omega \lt 1 [/math];
  • intermediate beaches correspond to [math]1\lt \Omega \lt 6 [/math];
  • dissipative beaches correspond to [math]\Omega \gt 6 [/math].

The tidal range [math]H_{tide}[/math] also influences the beach profile, because the surf zone is shifted across the shoreface with tidal level. For large tides (tidal range much larger than wave height) the surf zone covers during the tidal cycle a much larger area than for small tides. Beaches with large tides therefore have smaller upper shoreface slopes compared to beaches with the same [math]\Omega[/math] and small tides.

Several field studies[24][25][26] have been performed to relate beach profiles to different combinations of the Dean parameter [math]\Omega[/math] and the relative tidal range [math]RTR=H_{tide}/H_b[/math]. Distinction is made between different beach morphologies: reflective beaches with steep graveled slopes, reflective/intermediate beaches with subtidal bars, intermediate beaches with rip cells and low-tide terrace or low-tide bar, dissipative beaches without significant bar formations, dissipative beaches with multiple inter-tidal bars and dissipative beaches with tidal flats, see Fig. 4. The figure only gives a rough indication of beach morphologies, because the natural variability is not fully captured by the two dimensionless parameters [math]\Omega[/math] and [math]RTR[/math]. Other factors play a role, in particular the temporal variability in morphology and environmental conditions and the role of geology[24].

Another parameter that is often used for analysing the beach response to variations in the wave climate is the so-called surf similarity parameter that relates the wave steepness (wave height to wave length ratio) to the surf zone slope.


Difference between upper and lower shoreface

The upper shoreface is the zone where the energy of incident waves is converted into water level gradients, currents, seabed agitation and sediment transport. This goes with strong morphodynamic activity manifested in profile change, generation of bedforms and shoreline advance or retreat. In contrast, wave energy dissipation on the lower shoreface is small and hardly affects the wave energy flux toward the coast, [math]F^{(x)} = c_g H_s^2 \cos \alpha[/math], which can therefore be considered constant (Green's law; see Table Symbols for the symbols). Assuming approximate validity of linear wave theory, [math]c_g \approx c \approx \sqrt{gh}[/math], the height of waves propagating to shallower waters increases as [math]H_s \propto h^{-1/4}[/math]. While shoaling, these waves develop wave skewness and asymmetry on the lower shoreface (see Appendix). Wave skewness and asymmetry generate net sediment transport (see next paragraph) and therefore play an important role in shaping the cross-shore profile.

Sediment transport rates on the lower shoreface are much smaller than on the upper shoreface, implying a long adaptation time scale to changes in sea level and wave climate. On the upper shoreface these time scales are on the order of day-year, compared to decade-millennium time scales on the lower shoreface[1][29]. Besides, waves are not the only sediment transport agent on the lower shoreface; tidal currents and wind-driven currents are important as well. The strongest currents are associated with extreme events (storm surges, cyclones), but persistent wind stress over longer periods can produce up or downwelling and associated up or downslope sediment transport. At places where river plumes flowing along the coast induce a cross-shore salinity gradient, fine sediment can be trapped on the lower shoreface and even on the upper shoreface, as explained in the article Coastal mud belt. However, on exposed coasts, energetic waves will transport fine sediment from the upper to the lower shoreface. The lower shoreface can also receive sediment directly from nearby river deltas. When sediment sources are absent or sparse, the seabed sediment of the lower shoreface has mainly a relict origin with possibly rocky outcrops, indurated clay or coarse terrigenous deposits[29]. Biotic activity on the lower shoreface can be intense with a high production of calcareous bioclasts. The lower shoreface can therefore act as a source of calcareous sand to the upper shoreface and the beach[30].


Cross-shore sediment transport

Direct observation of net cross-shore sediment transport and its separate components in the field is difficult. This is particularly true for the surf zone, where the temporal and spatial variability of transport processes is very high. The hydrodynamic conditions for the installation and operation of measuring devices are harsh, especially in periods where sediment transport is most significant. Bottom-mounted pressure sensors are most often employed for the determination of wave characteristics. This can easily yield erroneous (underestimated) values if the nonlinearity of wave dynamics is not accounted for[31]. Acoustic and optical sensors are sensitive to air bubbles. Their use for determining suspended sediment concentrations and grainsizes is hampered as breaking waves entrain large volumes of air[32]. For these reasons, most knowledge is derived from laboratory experiments. Knowledge from the field is often obtained indirectly, through interpretation of the observed erosion and accretion patterns produced by cross-shore sediment transport. Mathematical models are often used as interpretation tools.

Sediment transport on the shoreface depends on many processes. The focus is on processes relevant to cross-shore transport. It should be realised, however, that longshore transport processes can also play an important role, for example by influencing the sediment balance of the shoreface or by generating various kinds of bedforms, see for example the articles: Natural causes of coastal erosion, Littoral drift and shoreline modelling and Rhythmic shoreline features. In the following we consider sandy beaches, with grain sizes typically between 0.2 and 0.5 mm (fine to medium sand). Coarse sedimentary beaches are dealt with in the article Gravel Beaches.

A distinction is often made between bedload transport and suspended load transport. Bedload transport represents the motion of sediment particles that are dragged by the flow along the seabed, either by rolling or successive jumps (saltation) or by entrainment of a thin seabed top layer (so-called sheet flow layer, of typically 10-30 grains thick). Suspended load represents the transport of sediment particles that have been brought into suspension by turbulent fluid motions. An important difference between the two transport modes is the phase between the load of the transported sediment and the instantaneous flow velocity, especially in the case of wave-induced flow. Bedload responds almost instantaneously to changes in the wave orbital velocity; increase (decrease) of velocity implies immediate increase (decrease) of bedload. This implies that an onshore (offshore ) flow asymmetry (so-called velocity skewness) with higher peak velocity in onshore (offshore) flow induces a net bedload transport in onshore (offshore) direction. Suspended load generally displays a fast response to an increase in flow velocity; sediment particles can be lifted quickly from the seabed and brought into suspension when the flow velocity increases. However, unlike bedload, the response of the suspended sediment concentration to a decrease in flow velocity depends on the settling velocity of sediment particles, which induces a response delay that can be substantial compared to the wave period. Asymmetries in onshore (offshore) wave orbital flow are therefore less decisive for the net direction of suspended load than for bedload. In the surf zone where the waves break, suspended load transport dominates. Flume experiments suggest that bedload transport dominates in situations where unbroken waves approach the shore, e.g. on the lower shoreface and on steep, reflective beaches[33].

As explained below, onshore transport on the shoreface is mainly related to bedload transport and offshore transport to suspended load transport. As bedload transport involves mainly the coarse sediment fraction whilst suspended load transport involves mainly the fine sediment fraction, it follows that under average conditions coarser sediment is preferentially transported onshore and finer sediment offshore[34].

Onshore wave-induced sediment transport

Wave transformation in shallow water is a major factor for onshore sediment transport [35]. Due to shoaling effects, starting on the lower shoreface, incident waves become increasingly skewed (asymmetry in onshore / offshore peak flow velocity) and asymmetric (asymmetry in onshore / offshore flow acceleration, also called acceleration skewness), see the Appendix for details. These asymmetries in the wave orbital velocity produce shear stresses on the seabed which are stronger in onshore direction than in offshore direction [36][37][38], see also Shallow-water wave theory. Wave skewness and asymmetry are primarily effective for inducing onshore bedload transport rather than suspended load transport, since bedload transport responds almost instantaneously to changes in near-bed velocities. Wave skewness and asymmetry are most pronounced on gently sloping dissipative coasts; they hardly develop on steep reflective coastal profiles.

Fig. 5. Sand transport over a steep vortex bed ripple in the Aberdeen oscillatory flow tunnel. The horizontal orbital motion is skewed: the velocity to the right (wave propagation direction) is stronger than to the left, but the duration is shorter. The high orbital velocity to the right creates a sand-laden vortex on the lee-side of the ripple crest, which is swept back over the ripple right after flow reversal. Courtesy: Jebbe van der Werf.

Phase lags between bed shear stress and suspended sediment concentration determine to which degree wave skewness induces onshore transport. In some situations (rippled bed, fine sands) a large phase lag can reverse the net direction of suspended sediment transport. A large time lag implies that sediment particles that are stirred up by turbulent eddies during onshore orbital motion are transported offshore after reversal of the wave orbital velocity[39][40][36][41], see Fig. 5. Although strongly skewed waves promote onshore ripple migration, the net sediment transport can be directed offshore. Onshore transport is stimulated by acceleration skewness[42]. It creates a sediment-laden vortex on the stoss side of bed ripples, which travels toward the onshore direction as it is carried by the onshore ambient flow[43]. Wave breaking that generates turbulent fluid motions and sediment resuspension under the wave crest also contributes to onshore transport, especially in the case of plunging waves[44]. But at the same time, wave breaking removes most of the wave skewness[45], thus reducing the associated onshore sediment transport component in the surf zone.

The development of forward streaming at the top of the wave boundary layer can provide a further contribution to onshore transport [46], but this only holds for smooth beds and not in cases where the seabed is rippled, in cases of strongly skewed waves or for sheet flow conditions[47][48] (see also Sediment transport formulas for the coastal environment). Due to the high sediment concentration in the streaming layer (and in spite of the small vertical extent of the layer and the low flow velocity), the contribution of wave driven near-bed streaming to onshore transport is significant, especially in the shoaling zone and surf zone prior the plunging point[49].

Offshore wave-induced sediment transport

Offshore transport is mainly due to so-called undertow: the compensating return current for the net onshore mass transport in the upper water column between wave trough and wave crest [50][51][52]. The offshore suspended sediment transport by the undertow current is particularly strong just landward of the location where waves are breaking (plunging waves) and where much sediment is brought into suspension by wave energy dissipation and by the longshore current[53][54]. Although the undertow velocity is smaller than the amplitude of the wave orbital velocity, it produces substantial offshore transport across the shoreface, in particular in the surf zone, where it often dominates the onshore transport related to wave-asymmetry effects. The undertow current is strongest near the bottom in the surf zone, whereas it peaks near the surface on the lower shoreface. On the lower shoreface it is largely offset by onshore-directed Stokes drift and therefore less effective for offshore transport[55][56].

Bound infragravity waves may also contribute to net offshore sediment transport[57], as discussed in the article Infragravity waves. The transport component related to infragravity waves can be quite substantial on dissipative beaches, but the net contribution is subdued after wave breaking because the phase relationship with incident waves is lost and the net sand transport can be either onshore or offshore directed[58]. However, small transport contributions cannot be ignored because the large transport components related to wave skewness / asymmetry and undertow largely cancel each other out[59]. For this reason, the downslope transport component must also be considered, see next section.

Downslope sediment transport

Several mechanisms can induce downslope sand transport. We do not consider here downslope turbidity currents or mud avalanching that may occur at the outer edge of deltas loaded with large amounts of fluvial sediment.

One mode of downslope transport consists of sand particles rolling down the shoreface slope. However, modeling suggests that this does not contribute significantly to downslope transport for fine to medium sandy beaches[1].

A second transport mode is so-called wave-supported gravity flow (WSGF). Observations on the Dutch coast show that wave stirring of the seabed under storm conditions can generate and sustain a high-density layer of suspended fine sand (thickness 15-20 cm, concentration exceeding 40 g/l) flowing down the shoreface slope (1/100 - 1/200) with a velocity of a few cm/s [60]. Modeling suggests that the high-density sediment layer is kept in suspension by the combined action of wave-induced turbulent stresses and hindered settling. The wave-supported gravity flow was observed during a storm with near-bed wave orbital velocities of about 0.5 m/s, lasted for more than a day and far exceeded in volume other cross-shore sand transport processes. However, the frequency of occurrence of WSGF events is not well known. Wave-supported gravity flows have also been observed on the shoreface of muddy coasts[61].

A third mechanism of downslope transport is related to the wave orbital excursion that moves sediment particles between zones of nearshore high bed shear stress and offshore low bed shear stress. Since resuspension after settling takes more time lower on the shoreface than higher on the shoreface, the downslope distance traveled by settling sediment particles is on average longer than the onshore distance. The average difference between offshore and onshore travel distance depends not only on the orbital excursion and shoreface slope, but also on the particle settling velocity. During the offshore excursion, sediment particles with a lower settling velocity settle later, thus lower, on the shoreface where the bed shear stresses are lower and the resuspension time longer. Fine sediment particles thus experience greater downslope transport than coarse particles. The slope effect requires a correction to the net sediment transport over a flat bottom, which is an increasing function of the product of bed slope [math]\tan \beta[/math] and wave orbital excursion, and a decreasing function of the fall velocity [math]w[/math]. Empirical formulas for suspended sediment transport assume in general that the slope effect can be represented by a correction factor to the instantaneous wave-induced transport components of the form[62]

[math](1- \gamma_s \Large \frac{u}{w} \normalsize ) \, \tan \beta, \qquad (4)[/math]

where [math]u[/math] is the instantaneous wave orbital velocity, [math]w[/math] the sediment fall velocity and [math]\gamma_s[/math] an adjustable parameter depending, among other things, on the resuspension processes. This is consistent with laboratory experiments, in which an approximate correspondence was found between coastal profiles and an adjusted Dean parameter[63]

[math]\Omega_{\beta} = \Large \frac{H_s \tan \beta}{wT_p} \normalsize . \qquad(5)[/math]

Experiments by Hattori and Kawamata (1980[64]) with different profile slopes show that onshore transport dominates for [math]\Omega_{\beta} \lt 0.08[/math] and offshore transport dominates for [math]\Omega_{\beta} \gt 0.08[/math]. An equilibrium slope corresponds to (order of magnitude)

[math]\Omega_{\beta} \sim 0.08. \qquad(6)[/math]

Comparison with Eq. (3) shows that reflective beaches have a steep slope and dissipative beaches a gentle slope. The presence of the fall velocity [math]w[/math] in the denominator of the downslope correction term in Eq. (4) also implies that net wave-induced transport can be downslope (offshore) for fine sediment even if upslope transport (onshore) prevails for coarse sediment.

Although not proven, it is often assumed that the slope effect is a major factor compensating for net wave-induced onshore transport (sum of all wave-induced onshore and offshore transport components including undertow) in situations close to equilibrium[65]. This implies that the equilibrium shoreface slope [math]\tan \beta [/math] will be steeper (larger [math]\tan \beta [/math]) for a coarser seabed (large fall velocity [math]w[/math]) than for a finer sediment bed (small fall velocity [math]w[/math], smaller [math]\tan \beta [/math]). This is expressed by the correction factor (Eq. 4) and is in general agreement with field observations. An analysis of long-term shoreface profile observations at the Gold Coast (Australia) suggests that the bed-slope effect on the total net sediment transport averaged over the wave period (or, more generally, the wave climate) can be related to the equilibrium bed slope [math]\tan \beta_{eq}[/math] by the correction factor [math](1 – \tan \beta / \tan \beta_{eq})[/math] [66], which is equivalent to Eq. (4).

Shoreface modelling

The processes involved in sediment transport on the shoreface, as described above, are complicated and not yet fully understood. Even the most sophisticated numerical morphodynamic models, like Delft3D and XBEACH [67], make use of approximate empirical formulations for several processes contributing to sediment transport and involve parameters that need to be tuned. Empirical formulas for wave-induced sediment transport on the shoreface that are often used in practice are given in the articles Sediment transport formulas for the coastal environment and Sand transport. These formulas are based on laboratory and field measurements for a large range of conditions. Rough estimates can be obtained with the empirical default parameters given in these articles. For more accurate site-specific estimates it is necessary to tune these parameters to field data. A more detailed discussion of shoreface modelling is presented in the article Modelling coastal hydrodynamics.


Equilibrium shoreface profile (sandy beach)

Fig. 6. The Bruun/Dean profile (Eq. 7) and the equivalent Bodge profile (Eq. 8).

For a given stable wave climate the shoreface profile tends after a certain (long) period to an equilibrium profile, as pointed out before. This equilibrium profile can be computed by applying sediment transport formulas to an initial profile for a period comprising all wave conditions according to local wave climate statistics. This requires heavy calculations that are moreover subject to large uncertainty margins. Therefore, equilibrium profiles derived from field observations are often used. By analyzing a large number of beach profiles of the Californian and Danish coasts [68] and the US Atlantic and Gulf coasts [69] Bruun (1954) and Dean (1977) derived an equilibrium profile (often called the Dean/Bruun profile) of the form (Fig. 6)

[math]h(x)=A \; x^{2/3}, \qquad(7)[/math]

where [math]h[/math] is the water depth [m], [math]x[/math] the cross-shore distance [m] ([math]x=0[/math] at the shoreline) and [math]A[/math] a coefficient depending on the sediment fall velocity [math]w[/math] [m/s]. Dean (1987[70]) empirically determined [math]A \approx 0.5 \; w^{0.44} [/math] [m1/3]. Another estimate was given by Kriebel et al. (1991[71]) [math]A \approx 2.25 \; (w^2/g)^{1/3}[/math] [m1/3]. In Eq. (7) it is implicitly assumed that the fall velocity is uniform across the shoreface profile.

Dean [69] showed that the exponent 2/3 is consistent with the assumption of constant wave energy dissipation per unit volume throughout the surf zone.

Because the Bruun/Dean profile has an infinite slope at the shoreline, an alternative form was proposed by Bodge (1992) [72],

[math]h(x)=B \; (1 – e^{-kx}) , \qquad(8)[/math]

where [math]k[/math] is an empirical coefficient in the range [math]500 \lt k^{-1} \lt 2000 [/math] [m] and where the coefficient [math]B[/math] is related to the beach slope [math]\tan \beta=Bk[/math] at the shoreline [math]x=0[/math].

Fig. 7. Theoretical equilibrium profile for the upper shoreface (Eq. 7) and lower shoreface (Eq. 9).

The Bruun/Dean profile provides in general a reasonable description of observed long-term averaged profiles of micro-mesotidal sandy coasts[73]. This holds in particular for the upper shoreface where breaking waves are the dominant driving agent for sediment transport. The 2/3 power law is generally less adequate for describing the profile of the lower shoreface. This is not surprising, as sediment transport here is mainly driven by non-breaking waves. In a simplified ideal situation, based on realistic assumptions, a theoretical profile was derived through application of the principle of minimum energy dissipation[74][75]. The theoretical profile also follows a power law (see Fig. 7):

[math]h(x) = \Big[ h_b^{1/n} + (h_c^{1/n} - h_b^{1/n}) \Large\frac{x - x_b}{x_c - x_b}\normalsize \Big]^n , \quad x_b \lt x \lt x_c , \qquad (9)[/math]

where [math]h_b=h(x_b)[/math] is the depth at the transition from upper to lower shoreface and [math]x_b[/math] the corresponding location, and [math]h_c=h(x_c)[/math] is the closure depth and [math]x_c[/math] the corresponding location. The exponent [math]n[/math] has a value in the range 0.3 - 0.4, depending on the formulation of the bottom friction[75]. Inman et al. (1993) [76] found best fits with exponents ranging between 0.3 and 0.45 from an analysis of beach profiles for the coasts of south California, North Carolina and the Nile delta.

For gravel beaches other cross-shore profiles are found, see Gravel Beaches.


Analytical equilibrium models

Because of the earlier mentioned difficulties to derive model estimates of the equilibrium shoreface profile, taking wave climate statistics fully into account, analytical models have been developed in which the wave climate is replaced by a single representative monochromatic incident wave. In these models an alongshore uniform shoreface is assumed with perpendicular wave incidence, and a simple wave transformation model is used based on depth saturation of broken waves (wave height limited by depth) and shallow-water wave theory. An analytical expression for the equilibrium shoreface slope can be derived by requiring that the total wave-integrated cross-shore sediment flux equals zero throughout the shoreface zone [65] [1][77]. Although these models reproduce certain properties also observed in the field, such as increasing profile steepness with increasing grain size and with increasing wave period, a more detailed comparison of computed profiles with observations exhibits large discrepancies. This is perhaps not surprising in view of the simplicity of the model used for wave-induced sediment transport, which assumes zero phase difference between wave-orbital velocity and sediment load, and which ignores many other processes described in the previous paragraphs. It has been shown[78] that some discrepancies can be repaired by using expressions for the sediment transport components derived from field observations.

Analytical models also ignore the influence of extreme storms, which can flatten the profile out to great depths within a short period. The profile adaptation time scale for the entire shoreface being quite long[1], it is plausible that observed shoreface profiles generally correspond to transient profiles [77].


Related articles

Nearshore sandbars
Bruun rule
Active coastal zone
Shoreline retreat and recovery
Coastal Hydrodynamics And Transport Processes
Modelling coastal hydrodynamics
Coastal and marine sediments
Gravel Beaches


Further reading

Komar, P.D. 1998. Beach processes and sedimentation. Prentice Hall, London, pp. 544.
Woodroffe, C.D. 2002. Coasts, form, processes and evolution. Cambridge Univ.Press, 623 pp.
Dronkers, J. 2017. Dynamics of Coastal Systems. World Scientific Publ. Co. 753 pp.


Table Symbols

Variable Description Variable Description
[math]c[/math] wave propagation speed ([math][m/s][/math]) [math]c_g[/math] wave group propagation speed ([math][m/s][/math])
[math]d_{50}[/math] mean sediment grainsize ([math][m][/math]) [math]g[/math] gravitational acceleration ([math]\approx 9.8 \; m/s^2[/math] )
[math]h[/math] water depth ([math][m][/math]) [math]h_b[/math] water depth at breaking location ([math][m][/math])
[math]H_b[/math] significant wave height at breaking location ([math][m][/math]) [math]H_s[/math] significant wave height ([math][m][/math])
[math]k[/math] wave number ([math][m^{-1}][/math]) [math]t[/math] time ([math][s][/math])
[math]T_p[/math] spectral peak wave period ([math][s][/math]) [math]u[/math] velocity ([math][m/s][/math])
[math]w[/math] sediment fall velocity ([math][m/s][/math]) [math]x[/math] distance from shoreline (seaward >0, [math][m][/math])
[math]\alpha[/math] wave incidence angle [math]\tan \beta[/math] shoreface slope
[math]\theta[/math] wave phase [math]\theta = \omega t - kx[/math] [math]\rho[/math] seawater density ([math][kg/m^3][/math])
[math]\omega[/math] radial wave frequency ([math][s^{-1}][/math]) [math]\Omega[/math] Dean parameter ([math]= \Large \frac{H_b}{wT_p} \normalsize [/math])


Appendix: Wave skewness

When a sinusoidal wave travels into shallow water the wave shape changes. The equations describing wave motion in deep water are almost linear; the solution is a simple sinus-shaped wave. The nonlinear terms in the wave equations become important in shallow water. These nonlinear terms introduce higher harmonic wave components in the solution. Initially, the first harmonic components grow the fastest. The wave orbital velocity can then be written as

[math]u(\theta) = U_1 \cos(\theta) + U_2 \cos(2\theta) + V_2 \sin(2\theta) , \qquad (A1)[/math]

where [math]\theta = \omega t - kx[/math] is the wave phase, with [math]t=[/math]time, [math]x=[/math]the distance travelled in the wave propagation direction, [math]\omega =[/math]the wave angular frequency, [math]k=\omega / c=[/math]the wave number and [math]c=[/math]the wave propagation speed. The [math]U_2[/math] component represents the velocity skewness [math]S_u[/math] and the component [math]V_2[/math] represents the wave asymmetry or acceleration skewness [math]A_u[/math] (Fig. A1).


Fig. A1. Wave orbital velocity as function of the phase [math]\theta[/math]. Left panel: Velocity skewness. Right panel: Acceleration skewness. Blue curves: Eq. (A1). Red curves: Eq. (A5). Positive acceleration skewness stimulates onshore bedload transport, because the fast increase of the onshore orbital motion leaves the boundary layer barely time to grow, leading to strong near-bed velocity gradients and shear stresses[79]. It has been suggested that high pressure forces under the steep wave front also play a role by mobilizing bed grains around flow reversal[80].


We define [math] \lt …\gt \equiv \Large\frac{1}{2 \pi}\normalsize \int_0^{2 \pi} … d \theta [/math] and assume [math]|U_2|,|V_2| \lt \lt |U_1| .[/math]

The velocity and acceleration skewness are defined as

[math]S_u \equiv \Large\frac{\lt u^3\gt }{\lt u^2\gt ^{3/2}}\normalsize \approx \Large\frac{3 U_2}{\sqrt{2} U_1}\normalsize, \quad A_u \equiv \Large\frac{ \lt (\mathcal{H}u)^3 \gt }{\lt (\mathcal{H}u)^2\gt ^{3/2} }\normalsize \approx \Large\frac{3 V_2}{\sqrt(2) U_1}\normalsize . \qquad (A2)[/math]

(the Hilbert transform [math]\mathcal{H}[/math] applied to sinusoidal functions gives: [math]\mathcal{H}cos(x)=-sin(x), \; \mathcal{H}sin(x)=cos(x), \; \mathcal{H}cos(2x)=-sin(2x), \; \mathcal{H}sin(2x)=cos(2x)[/math], etc).

By analysing a large number of field observations of the wave orbital velocity on the shoreface, Ruessink et al. (2012[81]) found empirical evidence that the velocity and acceleration skewness mainly depend on the Ursell number

[math]U_r = \Large\frac{3 g H}{8 h^2 \omega^2}\normalsize , \qquad (A3)[/math]

where [math]g \approx [/math]9.8 m s-2 is the gravitational acceleration, [math]H[/math] the significant wave height and [math]h[/math] the local depth. They established the empirical expressions

[math] (S_u^2+A_u^2)^{-1/2} = 1.17 \, (1 + 0.205 \, U_r^{-3.37}), \quad \Large\frac{A_u}{S_u}\normalsize = -\cot \Big(\Large\frac{\pi}{2}\normalsize \tanh \big(0.815 U_r^{-0.672}\big)\Big) \normalsize, \qquad (A4)[/math]

from which [math]S_u[/math] and [math]A_u[/math] can be solved. Besides the Ursell number, also other factors influencing the velocity and acceleration skewness have been identified, such as the offshore wave steepness, spectral band width and beach slope[82].

Abreu et al. (2010[37]) derived a more general expression than Eq. (A1) for the wave orbital velocity, that includes velocity skewness and acceleration skewness,

[math]u(\theta)=U_1 \Large\frac{\sqrt{1-r^2}}{1+\sqrt{1-r^2}} \, \frac{r \sin{\phi}+(1+\sqrt{1-r^2}) \sin(\theta)}{1 - r \cos(\theta +\phi)}\normalsize . \qquad (A5)[/math]

The parameters [math]r[/math] and [math]\phi[/math] are related to the velocity skewness and acceleration skewness by

[math]r = \Large\frac{2 \sqrt{p+1}}{p+2}\normalsize, \quad p = \Large\frac{9}{2 (S_u^2+A_u^2)}\normalsize , \quad \phi = -\Large\frac{\pi}{2}\normalsize -\arctan(A_u/S_u) . \qquad (A6)[/math]

A predictive model for wave-induced cross-shore sediment transport can be obtained by using the expression (A5) for the wave orbital velocity in formulas for cross-shore sediment transport, for example the formulas presented in the article Sediment transport formulas for the coastal environment. This yields only the transport contribution of the oscillatory wave component. Other components of cross-shore sediment transport discussed above have to be considered as well.


References

  1. 1.0 1.1 1.2 1.3 1.4 Stive, M.J.F. and De Vriend, H.J. 1995. Modelling shoreface profile evolution. Mar. Geol. 126: 235–248 Cite error: Invalid <ref> tag; name "S" defined multiple times with different content
  2. Hallermeier, R.J. 1981. A profile zonation for seasonal sand beaches from wave climate. Coast. Eng. 4: 253–277
  3. Valiente, N.G., Masselink, G., Scott, T., Conley, D. and McCarroll, R.J. 2019. Role of waves and tides on depth of closure and potential for headland bypassing. Mar. Geol. 407: 60–75
  4. Hallermeier, R.J. 1978. Uses for a calculated limit depth to beach erosion. In: 16th Coastal Engineering Conference, American Society of Civil Engineers, pp. 1493- 1512
  5. Wright, L.D., Chappel, J., Thom, B.G., Bradshow, M.P. and Cowell, P. 1979. Morphodynamics of reflective and dissipative beach and inshore systems: south-eastern Australia. Mar.Geol. 32: 105-140
  6. Niedoroda, A.W., Swift, D.J., Hopkins, T.S. and Ma, C.-M. 1984. Shoreface morphodynamics on wave-dominated coasts. Mar. Geol. 60: 331–354
  7. Larson, M., Kraus, N.C. and Wise, R.A. 1999. Equilibrium beach profiles under breaking and non-breaking waves. Coast. Eng. 36: 59–85
  8. Aragones, L., Pagan, J.I., Lopez, I. and Serra, J.C. 2018. Depth of closure: New calculation method based on sediment data. International Journal of Sediment Research 33: 198-207
  9. Patterson, D.C. and Nielsen, P. 2016. Depth, bed slope and wave climate dependence of long term average sand transport across the lower shoreface. Coastal Engineering 117: 113–125
  10. Wright, L.D., Short, A.D. and Green, M.O. 1985. Short-term changes in the morphodynamic states of beaches and surf zones: an empirical predictive model. Marine Geology 62: 339—364
  11. Ludka, B. C., Guza, R. T. O’Reilly, W. C. and Yates M. L. 2015. Field evidence of beach profile evolution toward equilibrium, J. Geophys. Res. Oceans 120: 7574–7597
  12. Eichentopf, S., Van der Zanden, J., Caceres, I. and Alsina, J.M. 2019. Beach profile evolution towards equilibrium from varying initial morphologies. J. Mar. Sci. Eng. 7, 406
  13. Cornaglia, P. 1889. Delle Spiaggie. Accademia Nazionale dei Lincei, Atti. Cl. Sci. Fis., Mat. e Nat. Mem. 5: 284-304
  14. Johnson, D.W. 1919. Shore processes and shoreline development. Prentice Hall, N-Y, 584 pp.
  15. Van Rijn, L.C. 1998. Principles of Coastal Morphology. Aqua Publications, The Netherlands (aqua publications.nl.
  16. Reniers, A.J.H.M., Thornton, E.B., Stanton, T.P. and Roelvink, J.A. 2004. Vertical flow structure during Sandy Duck: observations and modelling. Coast.Eng. 51: 237-260
  17. Jacobsen, N.G. and Fredsoe, J. 2014. Formation and development of a breaker bar under regular waves. Part 2: Sediment transport and morphology. Coastal Eng. 88: 55-68
  18. Sleath, J.F.A. 1984. Sea bed mechanics. Wiley, New York.
  19. Fredsøe, J., Andersen, K.H. and Sumer, M.B. 1999. Wave plus current over a ripple-covered bed. Coastal Eng. 38: 177-221
  20. Marieu, V., Bonneton, P., Foster, D.L. and Ardhuin, F. 2008. Modeling of vortex ripple morphodynamics. J. Geophys. Res. 113, C09007, doi:10.1029/2007JC004659
  21. Bujan, N., Cox, R. and Masselink, R. 2019. From fine sand to boulders: Examining the relationship between beach-face slope and sediment size. Marine Geology 417, 106012
  22. Gourlay, M. R. 1968. Beach and dune erosion tests. Rep. m935/m936, Delft Hydraul. Lab., Delft
  23. Dean, R.G. 1973. Heuristic models of sand transport in the surf zone. Proc. Conf. Eng. Dynamics in the Surf Zone, Sydney: 208-214.
  24. 24.0 24.1 24.2 Scott, T., Masselink, G. and Russell, P. 2011. Morphodynamic characteristics and classification of beaches in England and Wales. Mar. Geol. 286: 1–20
  25. 25.0 25.1 Short, A.D. and Jackson, D.W.T. 2013. Beach morphodynamics. In: Shroder, J.F. (Ed.), Treatise in Geomorphology, vol. 10. Academic Press, San Diego
  26. 26.0 26.1 Jackson, D.W.T., Short, A.D., Loureiro, C. and Cooper, J.A.G. 2022. Beach morphodynamic classification using high-resolution nearshore bathymetry and process-based wave modelling. Estuarine, Coastal and Shelf Science 268, 107812
  27. Li, Y., Zhang, C., Zhao, S., Qi, H., Cai, F. and Zheng, J. 2024. Equilibrium configurations of sandy-muddy transitional beaches on South China coasts: Role of waves in formation of sand-mud transition boundary. Coastal Engineering 187 (2024) 104401
  28. Wright, L.D. and Short, A.D. 1984. Morphodynamic variability of surf zones and beaches: a synthesis. Mar.Geol. 56: 93-118
  29. 29.0 29.1 Anthony, E.J. and Aagaard, T. 2020. The lower shoreface: morphodynamics and sediment connectivity with the upper shoreface and beach. Earth Sci. Rev. 210: 103334
  30. Short, A.D. 2013. Australia’s temperate carbonate coast: Sources, depositional environments and implications. In: Martini, I.P., Wanless, H.R. (Eds.), Sedimentary Coastal Zones from High to Low Latitudes: Similarities and Differences. London, Special Publications, Geological Society
  31. Martins, K., Bonneton, P., Mouragues, A. and Castelle, B. 2020. Non-hydrostatic, non-linear processes in the surf zone. Journal of Geophysical Research: Oceans 125, e2019JC015521
  32. Aagaard, T. and Hughes, M.G. 2022. Sediment Transport. Chapter 8.6, Treatise on Geomorphology 2nd edition. Elsevier
  33. van der Zanden, J., van der A, D. A., Hurther, D., Cáceres, I., O'Donoghue, T., Hulscher, S. J. M. H. and Ribberink, J. S. 2017. Bedload and suspended load contributions to breaker bar morphodynamics. Coastal Engineering 129: 74–92
  34. Hassan, W.N. and Ribberink, J.S. 2005. Transport processes of uniform and mixed sands in oscillatory sheet flow. Coastal Engineering 52: 745– 770
  35. Elgar, S., Gallagher, E.L. and Guza, R.T. 2001.Nearshore sandbar migration. J.Geophys.Res. 106: 11,623-11,727
  36. 36.0 36.1 Nielsen, P. 1992. Coastal bottom boundary layers and sediment transport. In: Advanced Series on Ocean Engineering, IV. World Scientific.
  37. 37.0 37.1 Abreu, T., Silva, P.A., Sancho, F. and Temperville, A. 2010. Analytical approximate wave form for asymmetric waves. Coastal Eng. 57: 656-667
  38. Van der A, D.A., O’Donoghue, T., Davies, A.G. and Ribberink, J.S. 2011. Experimental study of the turbulent boundary layer in acceleration-skewed oscillatory flow. J. of Fluid Mech. 684: 251-283
  39. Vincent, C.E. and Green. M.O. 1990. Field measurements of the suspended sand concentration profiles and fluxes and of the resuspension coefficient over a rippled bed. J.Geophys.Res. 95: 11591-11601
  40. Ribberink, J.S. and Al-Salem, A.A. 1995. Sheet flow and suspension of sand in oscillatory boundary layer. Coastal Eng. 25: 205-225
  41. Ruessink, B.G., Michallet, H., Abreu, T., Sancho, F., van der A, D.A., van der Werf, J.J. and Silva, P.A., 2011. Observations of velocities, sand concentrations, and fluxes under velocity-asymmetric oscillatory flows. J. Geophys. Res. 116, C03004, doi:10.1029/2010JC006443
  42. Dubarbier, B., Castelle, B., Marieu, V. and Ruessink, G. 2015. Process-based modelling of cross-shore sandbar behavior. Coastal Engineering 95: 35–50
  43. Salimi-Tarazouj, A., Tian-Jian Hsu, T-J., Traykovski, P. and Chauchat, J. 2024. Investigating wave shape effect on sediment transport over migrating ripples using an Eulerian two-phase model. Coastal Engineering 189, 104470
  44. Aagaard, T., Hughes, M.G. and Ruessink, B.G. 2018. Field observations of turbulence, sand suspension and cross-shore transport under spilling and plunging breakers. Journal of Geophysical Research - Earth Surface 123: 2844–2862
  45. Guza, R.T. and Thornton, E.B. 1985. Velocity moments in the nearshore. J. Waterways Port Coastal Ocean Div. ASCE III: 235-256
  46. Longuet-Higgins, M.S. 1953. Mass transport in water waves. Royal Soc. London, Phil.Trans. 245A: 535-581
  47. Trowbridge, J.H. and Madsen, O.S. 1984. Turbulent wave boundary layers: 2. Second-order theory and mass transport. J.Geophys.Res. 89: 7999-8007
  48. Kranenburg, W. M., Ribberink, J. S., Uittenbogaard, R. E. and Hulscher, S. J. M. H. 2012. Net currents in the wave bottom boundary layer: on wave shape streaming and progressive wave streaming. J. Geophys. Res. 117(F03005), doi:10.1029/2011JF002070
  49. Zheng, P., Gumbira, G., Li, M., van der Zanden, J., van der A, D., van der Werf, J., Chen, X. and Tang, X. 2023. Development, calibration and validation of a phase-averaged model for cross-shore sediment transport and morphodynamics on a barred beach. Continental Shelf Research 258, 104989
  50. Stive, M.J.F. and Wind, H.G. 1986. Cross-shore mean flow in the surf zone. Coastal Eng. 10: 325-340
  51. Cox, D.T. and Kobayashi, N. 1998. Application of an undertow model to irregular waves on plane and barred beaches. J.Coast.Res. 14: 1314-1324
  52. Gallagher, E.L., Elgar, S. and Guza, R.T. 1998. Observations of sand bar evolution on a natural beach. J. Geophys. Res. 103: 3203–3215
  53. Thornton, E.B., Humiston, R.T. and Birkemeier, W. 1996. Bar/trough generation on a natural beach. J. Geophys. Res. 101: 21-97-12110
  54. Reniers, A.J.H.M., Thornton, E.B., Stanton, T.P. and Roelvink, J.A. 2004. Vertical flow structure during Sandy Duck: observations and modelling. Coast.Eng. 51: 237-260
  55. Lentz, S.J. and Fewings, M. 2012. The wind- and wave-driven inner-shelf circulation. Annual Review of Marine Science 4: 317–341
  56. Larsen, B.E. and Fuhrman, D.R. 2023. Simulation of cross-shore breaker bar development utilizing a stabilized two-equation turbulence model. Coastal Engineering 180, 104269
  57. Roelvink, J.A. and Stive, M.J.F. 1989. Bar generating cross-shore flow mechanisms on a beach. J. Geophys. Res. 94: 4785-4800
  58. Aagaard, T., Greenwood, B. and Hughes, M. 2013. Sand transport on dissipative, intermediate and reflective beaches. Earth Sci. Rev. 124: 32–50
  59. Ruessink, B.G., Houwman, K.T. and Hoekstra, P. 1998. The systematic contribution of transporting mechanisms to the cross-shore sediment transport in water depths of 3 to 9 m. Mar. Geol. 152: 295–324
  60. Flores, R. P., Rijnsburger, S., Meirelles, S., Horner-Devine, A. R., Souza, A. J., Pietrzak, J. D., Henriquez, M. and Reniers, A. 2018. Wave Generation of Gravity-Driven Sediment Flows on a Predominantly Sandy Seabed. Geophysical Research Letters 45: 7634-7645
  61. Peng, Y., Yu, Q., Du, Z., Wang, L, Wang, Y and Gao, S. 2022, Gravity-driven sediment flows on the shallow sea floor of a muddy open coast. Marine Geology 445, 106759
  62. Bagnold, R.A. 1963. Mechanics of marine sedimentation. In: The sea, vol.3. Ed. M.N.Hill, Wiley-Interscience: 507-528
  63. Bayle, P.M., Beuzen, T., Blenkinsopp, C.E., Baldock, T.E. and Turner, I.L. 2021. A new approach for scaling beach profile evolution and sediment transport rates in distorted laboratory models. Coastal Engineering 163, 103794
  64. Hattori, M. and Kawamata, R. 1980. Onshore–offshore transport and beach profile change. In: Proc. 17th ASCE Coastal Engineering Conference, 2, pp. 1175–1193
  65. 65.0 65.1 Bowen, A.J. 1980. Simple models of nearshore sedimentation; beach profiles and longshore bars. In: The coastline of Canada, Ed. by S.B. McCann, Geological Survey of Canada, Ottawa, pp. 1-11
  66. Patterson, D.C. and Nielsen, P. 2016. Depth, bed slope and wave climate dependence of long term average sand transport across the lower shoreface. Coastal Engineering 117: 113–125
  67. https://oss.deltares.nl/web/xbeach/home
  68. Bruun, P. 1954. Coast erosion and the development of beach profiles. Beach Erosion Board, US Army Corps of Eng., Tech.Mem. 44: 1-79
  69. 69.0 69.1 Dean, R.G. 1977. Equilibrium beach profiles: US Atlantic coast and Gulf coasts. Ocean Eng.Tech.Rep. 12, Univ. of Delaware, Newark, 45 pp.
  70. Dean, R. 1987. Coastal sediment processes: Toward coastal engineering solutions. Spec. Conf. Understanding Coastal Sediment Processes. N.C. Kraus, ed. ASCE, New York, pp. 1-24
  71. Kriebel, D.L., Kraus, N.C. and Larson, M. 1991. Engineering Methods for Predicting Beach Profile Response. Procs. Conf. Quantitative Approaches to Coastal Sediment Processes. N.C. Kraus, K.J. Gingerich, D.L. Kriebel eds. ASCE, New York, pp. 557-571
  72. Bodge, K.R. 1992. Representing equilibrium beach profiles with an exponential expression. J. Coastal Research 8: 47-55.
  73. Bernabeu, A.M., Medina, R. and Vidal, C. 2003. A morphological model of the beach profile integrating wave and tidal influences. Mar.Geol. 197: 95-116
  74. Larson, M., Kraus, N. C. and Wise, R. A. 1999. Equilibrium beach profiles under breaking and non-breaking waves. Coastal Engineering 36: 59–85
  75. 75.0 75.1 Maldonado, S. 2020. Do beach profiles under nonbreaking waves minimize energy dissipation? Journal of Geophysical Research: Oceans 125, e2019JC015876
  76. Inman, D.L., Elwany, M.H. and Jenkins, S.A. 1993. Shorerise and bar-berm profiles on ocean beaches. J.Geophy.Res. 98: 18181-18199
  77. 77.0 77.1 Ortiz, A.C. and Ashton, A.D. 2016. Exploring shoreface dynamics and a mechanistic explanation for a morphodynamic depth of closure. J. Geophys. Res. Earth Surf 121: 442–464 doi:10.1002/2015JF003699
  78. Aagaard, T. and Huges, M.G. 2017. Equilibrium shoreface profiles: A sediment transport approach. Marine Geology 390: 321–330
  79. Van der A, D.A., O’Donoghue, T. and Ribberink, J.S. 2010. Measurements of sheet flow transport in acceleration-skewed oscillatory flow and comparison with practical formulations. Coastal Engineering 57: 331–342
  80. Calantoni, J. and Puleo, J.A. 2006. Role of pressure gradients in sheet flow of coarse sediments under sawtooth waves. Journal of Geophysical Research 111, 2005JC002875
  81. Ruessink, B.G., Ramaekers, G. and van Rijn, L.C. 2012. On the parameterization of the free-stream non-linear wave orbital motion in nearshore morphodynamic models. Coastal Engineering 65: 56–63
  82. Rocha, M.V.L., Michallet, H. and Silva, P.A. 2017. Improving the parameterization of wave nonlinearities – The importance of wave steepness, spectral bandwidth and beach slope. Coastal Engineering 121: 77–89


The main author of this article is Job Dronkers
Please note that others may also have edited the contents of this article.

Citation: Job Dronkers (2024): Shoreface profile. Available from http://www.coastalwiki.org/wiki/Shoreface_profile [accessed on 31-10-2024]