Shoreface profile

From MarineBiotech Infopages
Jump to: navigation, search
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


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[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). Other criteria for the outer closure depth can be derived from the condition that the wave-induced bed shear stress should not exceed the threshold of sediment transport under extreme storms[3]. Typical values of the outer closure depth are in the range 15 – 40 m.

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

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.

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, 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[9][10][11]. 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[12][13]. 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 [14]. 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 armouring'.

Figure 2. The barred coast of North Holland, The Netherlands., 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) 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 [15][16] - 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[17], vortex shedding at ripple crests[18] and non-linear ripple-ripple interactions[19]. 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[20]. 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 [21][22]:

[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. 3. Beach morphologies according to the parameters [math]\Omega[/math] and [math]RTR[/math], based on beach surveys in England, Wales, Northern Ireland and Australia[23][24][25]. Reflective beaches are found in the red zone and dissipative beaches in the blue zone. Almost no natural beaches are found in the white zone.

Wright and Short (1984) [26] 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[23][24][25] 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. 3. 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[23].

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.

Cross-shore sediment transport

Direct observation of net cross-shore sediment transport and its separate components in the field is very difficult, in particular in the surf zone. Experimental data seldom yield conclusive results, for several reasons: the hydrodynamic conditions for the installation and operation of measuring devices are harsh, especially in periods where sediment transport is strongest, the temporal and spatial variability of the transport processes is very high, and sediment concentration sensors cannot be easily calibrated because of their sensitivity to air bubbles (especially acoustic backscatter sensors) and grainsize (especially optical backscatter sensors)[27]. 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. We focus on processes relevant for 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 sheetflow layer, of about 10 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 flow velocity; increase (decrease) of flow 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. Bedload transport can dominate in situations where unbroken waves approach the shore, e.g. on the lower shoreface and on steep, reflective beaches[27].

Onshore wave-induced sediment transport

Wave transformation in shallow water is a major factor for onshore sediment transport [28]. 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 [29][30][31], see also Shallow-water wave theory. Wave skewness and asymmetry are most pronounced on gently sloping dissipative coasts; they hardly develop on steep reflective coastal profiles.

Phase lags between bed shear stress and suspended sediment concentration determine to which degree wave skewness induces onshore transport. In some situations (rippled bed, sheet flow) 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 [32][33][29][34]). Onshore transport is stimulated by acceleration skewness[35] and by wave breaking that generates turbulent fluid motions and sediment resuspension under the wave crest, especially for plunging waves[36]. But at the same time, wave breaking removes most of the wave skewness[37], 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 also contributes to onshore transport [38], 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 [39][40] (see also Sediment transport formulas for the coastal environment). Another, but generally lesser contribution to onshore sediment transport is further due to net mass transport in the upper part of the water column between wave trough and wave crest (in-phase variation of wave water level and wave orbital velocity).

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 upper water column between wave trough and wave crest [41][42][43]. The offshore undertow sand transport is particularly strong just landward of the location where waves are breaking and where much sediment is brought into suspension by wave energy dissipation and by the longshore current[44][45]. Although the undertow current is smaller than the wave orbital velocities, it produces substantial offshore transport across the shoreface, and in particular in the surf zone, where it can dominate the onshore transport related to wave-asymmetry effects. Bound infragravity waves may also contribute to net offshore sediment transport[46], 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[47]. However, small transport contributions cannot be ignored because the large transport components related to wave skewness / asymmetry and undertow largely cancel each other out[48]. For this reason, the downslope transport component must also be considered, see next section.

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 [1] (we do not consider here downslope turbidity currents or mud avalanching occurring at the outer edge of river deltas that receive large amounts of sediment) . Gravity-induced sediment transport is generally 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[1]. When a sediment particle is lifted from the bed it travels in onshore or in offshore direction depending on the phase of the wave orbital motion. When it settles higher on the shoreface it is more easily resuspended than when it settles lower on the shoreface, because it is subject to greater bed shear stresses in the former case. This implies that settling high on the shoreface will occur more frequently for coarser sediment particles with high fall velocity, while finer sediment particles with small fall velocity will preferentially settle lower on the shoreface. The slope effect thus requires a correction to the net sediment transport over a flat bottom, which is an increasing function of the bed slope [math]\tan \beta[/math] and a decreasing function of the fall velocity [math]w[/math]. 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 [49]

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

where [math]u[/math] is the instantaneous wave orbital velocity, [math]w[/math] the sediment fall velocity and [math]\alpha[/math] an efficiency parameter. This is consistent with laboratory experiments, in which an approximate correspondence was found between coastal profiles and an adjusted Dean parameter[50]

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

Experiments by Hattori and Kawamata (1980[51]) 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 wave-induced onshore and offshore transport components) in situations close to equilibrium[52]. 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 (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] [53], 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 XBEACH [54], 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)

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. 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 [55] and the US Atlantic and Gulf coasts [56] Bruun (1954) and Dean (1977) derived an equilibrium profile (often called the Dean/Bruun profile) of the form (Fig. 4)

Figure 4. The Bruun/Dean profile (Eq. 7) and the equivalent Bodge profile (Eq. 8).

[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[57]) empirically determined [math]A \approx 0.5 \; w^{0.44} [/math] [m1/3]. Another estimate was given by Kriebel et al. (1991[58]) [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 [56] 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) [59],

[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].

From an analysis of terrace-shaped shoreface profiles in northern Spain, Bernabeu et al. (2003) [60] found that the Bruun/Dean profile was adequate only for the upper shoreface; for the lower shoreface they proposed the form

[math]h(x)=A' \; (x-x_0)^{2/3} , \qquad(9)[/math]

where [math]x_0[/math] is determined such that the two profiles (7) and (9) match at the inflexion point between the upper and lower shoreface.

Inman et al. (1993) [61] 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, which is more consistent with uniform wave energy dissipation in the surf zone per unit surface rather than per unit volume.

The cross-shore profile of gravel beaches is not the same as for sandy beaches, see: Gravel Beaches.

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. Waves also play an important role on the lower shoreface in shaping the cross-shore profile while wave interaction with the seabed generates wave skewness and asymmetry and resulting net sediment transport. But transport rates on the lower shoreface are typically small, 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][62]. 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 in the nearshore zone, as explained in the article Coastal mud belt. However, fine sediment is more generally transported from the upper to the lower shoreface, as the medium sediment grainsize is typically larger on the upper shoreface than on 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[62]. 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[63].

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. 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 [52] [1][64]. 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[65] 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 [64].

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.

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

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), \quad \mathcal{H}cos(2x)=sin(2x), \quad \mathcal{H}sin(2x)=-cos(2x)[/math]).

By analysing a large number of field observations of the wave orbital velocity on the shoreface, Ruessink et al. (2012[66]) 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.

Abreu et al. (2010[30]) 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.


  1. 1.0 1.1 1.2 1.3 1.4 1.5 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. 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
  10. 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
  11. 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
  12. Cornaglia, P. 1889. Delle Spiaggie. Accademia Nazionale dei Lincei, Atti. Cl. Sci. Fis., Mat. e Nat. Mem. 5: 284-304
  13. Johnson, D.W. 1919. Shore processes and shoreline development. Prentice Hall, N-Y, 584 pp.
  14. Van Rijn, L.C. 1998. Principles of Coastal Morphology. Aqua Publications, The Netherlands (aqua
  15. 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
  16. 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
  17. Sleath, J.F.A. 1984. Sea bed mechanics. Wiley, New York.
  18. Fredsøe, J., Andersen, K.H. and Sumer, M.B. 1999. Wave plus current over a ripple-covered bed. Coastal Eng. 38: 177-221
  19. 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
  20. 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
  21. Gourlay, M. R. 1968. Beach and dune erosion tests. Rep. m935/m936, Delft Hydraul. Lab., Delft
  22. Dean, R.G. 1973. Heuristic models of sand transport in the surf zone. Proc. Conf. Eng. Dynamics in the Surf Zone, Sydney: 208-214.
  23. 23.0 23.1 23.2 Scott, T., Masselink, G. and Russell, P. 2011. Morphodynamic characteristics and classification of beaches in England and Wales. Mar. Geol. 286: 1–20
  24. 24.0 24.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
  25. 25.0 25.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
  26. Wright, L.D. and Short, A.D. 1984. Morphodynamic variability of surf zones and beaches: a synthesis. Mar.Geol. 56: 93-118
  27. 27.0 27.1 Aagaard, T. and Hughes, M.G. 2022. Sediment Transport. Chapter 8.6, Treatise on Geomorphology 2nd edition. Elsevier
  28. Elgar, S., Gallagher, E.L. and Guza, R.T. 2001.Nearshore sandbar migration. J.Geophys.Res. 106: 11,623-11,727
  29. 29.0 29.1 Nielsen, P. 1992. Coastal bottom boundary layers and sediment transport. In: Advanced Series on Ocean Engineering, IV. World Scientific.
  30. 30.0 30.1 Abreu, T., Silva, P.A., Sancho, F. and Temperville, A. 2010. Analytical approximate wave form for asymmetric waves. Coastal Eng. 57: 656-667
  31. 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
  32. 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
  33. Ribberink, J.S. and Al-Salem, A.A. 1995. Sheet flow and suspension of sand in oscillatory boundary layer. Coastal Eng. 25: 205-225
  34. 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
  35. Dubarbier, B., Castelle, B., Marieu, V. and Ruessink, G. 2015. Process-based modelling of cross-shore sandbar behavior. Coastal Engineering 95: 35–50
  36. 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
  37. Guza, R.T. and Thornton, E.B. 1985. Velocity moments in the nearshore. J. Waterways Port Coastal Ocean Div. ASCE III: 235-256
  38. Longuet-Higgins, M.S. 1953. Mass transport in water waves. Royal Soc. London, Phil.Trans. 245A: 535-581
  39. 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
  40. 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
  41. Stive, M.J.F. and Wind, H.G. 1986. Cross-shore mean flow in the surf zone. Coastal Eng. 10: 325-340
  42. 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
  43. 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
  44. Thornton, E.B., Humiston, R.T. and Birkemeier, W. 1996. Bar/trough generation on a natural beach. J. Geophys. Res. 101: 21-97-12110
  45. 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
  46. Roelvink, J.A. and Stive, M.J.F. 1989. Bar generating cross-shore flow mechanisms on a beach. J. Geophys. Res. 94: 4785-4800
  47. Aagaard, T., Greenwood, B. and Hughes, M. 2013. Sand transport on dissipative, intermediate and reflective beaches. Earth Sci. Rev. 124: 32–50
  48. 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
  49. Bagnold, R.A. 1963. Mechanics of marine sedimentation. In: The sea, vol.3. Ed. M.N.Hill, Wiley-Interscience: 507-528
  50. 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
  51. Hattori, M. and Kawamata, R. 1980. Onshore–offshore transport and beach profile change. In: Proc. 17th ASCE Coastal Engineering Conference, 2, pp. 1175–1193
  52. 52.0 52.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
  53. 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
  55. Bruun, P. 1954. Coast erosion and the development of beach profiles. Beach Erosion Board, US Army Corps of Eng., Tech.Mem. 44: 1-79
  56. 56.0 56.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.
  57. 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
  58. 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
  59. Bodge, K.R. 1992. Representing equilibrium beach profiles with an exponential expression. J. Coastal Research 8: 47-55.
  60. 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
  61. 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
  62. 62.0 62.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
  63. 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
  64. 64.0 64.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
  65. Aagaard, T. and Huges, M.G. 2017. Equilibrium shoreface profiles: A sediment transport approach. Marine Geology 390: 321–330
  66. 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

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 (2023): Shoreface profile. Available from [accessed on 21-03-2023]