Difference between revisions of "Breaker index"

From Coastal Wiki
Jump to: navigation, search
(One intermediate revision by the same user not shown)
Line 1: Line 1:
{{Definition|title= Breaker index
 
|definition= The ratio of wave height and still water depth at the shoreface location where waves start breaking.  }}
 
  
  
==Notes==
 
The hydrodynamics in the surf zone and the associated sediment transport largely depend on the process of wave breaking. Waves break when the steepness of the wave front exceeds a certain limit. The wave breaking process is very complex. Hence, instead of solving the underlying dynamical equations, it is represented in most models by empirical formulas. These models therefore make use of the breaker index concept. This concept is based on observations showing that the wave height progressively decreases with depth in the surf zone. The strongest assumption is that the ratio of wave height and still water depth is constant throughout the surf zone,
 
  
<math>\large\frac{H(x)}{h(x)}\normalsize = \gamma = constant , \qquad (1)</math>
+
The process of wave breaking plays a crucial role in the morphodynamics of the shoreface and the beach, such as accretion or erosion trends and the development and dynamics of nearshore bars and rip currents. The breaker index concept has been introduced to estimate the location on the shoreface where wave breaking occurs based on empirical formulas, without the need for a detailed simulation of the breaker process.
  
where <math>H</math> is the wave height, <math>h</math> the still water depth and <math>\gamma</math> the breaker index. This assumption is called '''saturated wave breaking'''.
+
==A practical wave breaking criterion==
  
This assumption is too rough for most applications. In a more refined approximation the two parameters wave steepness <math>H/L</math> and the shoreface bedslope <math>m=\tan \beta</math> are taken into consideration and the breaker index is evaluated at the breakpoint indicated by the subscript <math>b</math>,
+
According to numerical simulations<ref>Barthelemy, X., Banner, M., Peirson, W., Fedele, F., Allis, M. and Dias, F. 2018. On a unified breaking onset threshold for gravity waves in deep and intermediate depth water. Journal of Fluid Mechanics 841: 463-488</ref> and flume experiments<ref>Saket, A., Peirson, W., Banner, M., Barthelemy, X. and Allis, M. 2017. Wave breaking onset of two-dimensional deep-water wave groups in the presence and absence of wind. J. Fluid Mech. 811: 642–658 </ref>, wave breaking starts on the shoreface when the horizontal surface water particle velocity at the wave crest comes close (within about 15%) to the wave crest celerity of the tallest wave in the incident wave train. However, it is not easy to derive from this criterion at which location on the shoreface wave breaking will occur. Even if the characteristics of the incident waves and the shoreface profile are well known, one has to face the problem that wave breaking is highly sensitive to fine-scale processes. It cannot be easily modelled, although significant progress has been achieved with numerical techniques such as [https://en.wikipedia.org/wiki/Computational_fluid_dynamics Computational Fluid Dynamics (CFD)]<ref name=A>Aggarwal, A., Chella, M.A., Bihs, H. and Myrhaug, D. 2020. Properties of breaking irregular waves over slopes. Ocean Engineering 216, 108098</ref> or [https://en.wikipedia.org/wiki/Smoothed-particle_hydrodynamics Smoothed Particle Hydrodynamics (SPH)]<ref> Lowe, R.J., Buckley, M.L., Altomare, C., Rijnsdorp, D.P., Yao, Y., Suzuki, T. and Bricker, J.D. 2019  Numerical simulations of surf zone wave dynamics using Smoothed Particle Hydrodynamics. Ocean Modelling 144, 101481</ref>. For this reason, many investigations have been conducted to establish empirical rules for wave breaking from which the location of wave breaking can be easily derived. These empirical rules refer to the so-called breaker index. The breaker index relates the wave height at breaking, <math>H_b</math>, to the water depth <math>h_b</math> at the breakpoint; it is defined as
  
<math>\large\frac{H_b}{h_b}\normalsize = \gamma_b .   \qquad (2)</math>
+
<math>\gamma_b=H_b/h_b . \qquad (1)</math>
  
It may be expected that shorter waves will break earlier on the shoreface than longer waves (i.e. at greater depth for the same wave height), and that wave running over a steep slope will break later (i.e. at smaller depth) than waves running over a gentle slope. The breaker index should thus be an increasing function of the wavelength and the bed slope.  
+
Important parameters for the wave breaking process are the steepness <math>H_0/L_0</math> of the incident waves (<math>H_0, L_0</math> are the offshore wave height and wavelength) and the average shoreface slope <math>m</math>. Empirical rules establish a relationship between the breaker index and these parameters.  
  
 +
==Wave breaking process==
  
[[File:BreakerIndex.jpg|thumb|left|350px|Fig. 1. The area covered by experimental data for the breaker index as a function of the surf similarity parameter. The density of the experimental data is represented by varying degrees of blueness. The black line is the Battjes formula Eq. 4.  Figure adapted from Ostendorf and Madsen (1979<ref name=O>Ostendorf, D. and Madsen, O.  1979. An Analysis of Longshore Current and Associated Sediment Transport in the Surf Zone. Boston: Massachusetts Institute of Technology, Department of Civil Engineering Technical Report 241, 169p.</ref>).]]
+
[[File:WaveBreakingTypes.jpg|thumb|center|800px|Fig. 1. Schematic of wave evolution towards breaking, for swell waves (left panel) and sea waves (right panel).]]
  
Experiments by Battjes (1974<ref name=B>Battjes, J.A. 1974. Surf similarity. Proceedings 14th International Conference on Coastal Engineering, pp. 466–480</ref>) suggested that this dependence on wavelength and bed slope could be represented by the [[surf similarity parameter]]
+
The process of wave breaking starts on the shoaling zone where incident waves become progressively skewed by interaction with the seabed: the wave crest becomes sharper and higher, while the wave trough widens and flattens (Fig. 1). Incident waves with short wavelengths (locally generated 'sea') are steeper than incident waves with longer wavelengths (remotely generated swell); the former category of waves will break earlier, thus at greater depth <math>h_b</math> than the latter category<ref name=A/><ref name=G10>Goda, Y. 2010. Reanalysis of regular and random breaking wave statistics. Coastal Engineering Journal 52: 71-106</ref>, implying a larger value of  <math>\gamma_b</math>. These short sea waves are also less skewed before breaking; they break when the crest becomes unstable and flows down the front face of the wave, producing spilling breaker bores that surf towards the shoreline (Fig. 2). This breaking mode is also favoured on gently sloping shorefaces when the [[surf similarity parameter]] <math>\xi=m/\sqrt{H_0/L_0}</math> has values typically smaller than about 0.4. Incident waves with longer wavelengths will break later, thus at smaller depth <math>h_b</math>. While shoaling, they are strongly skewed and break when the crest curls over the front face and falls into the base of the wave, producing a so-called plunging breaker<ref name=CL>Camenen, B. and Larson, M. 2007. Predictive Formulas for Breaker Depth Index and Breaker Type. Journal of Coastal Research 23: 1028–1041</ref> (Fig. 3). This breaking mode is also favoured on steep sloping shorefaces when the [[surf similarity parameter]] has values typically larger than 0.4.
  
<math>\xi=m / \sqrt{H/L} ,  \qquad (3)</math>
 
  
according to the formula
+
{| border="0" align="center"
 +
|-
 +
| valign="top"|
 +
[[File:SpillingWave2.jpg|thumb|left|300px|Fig. 2. Spilling waves. Photo credit Andrew Dawley Flickr Creative Commons.]]
 +
| valign="top"|
 +
[[File:PlungingWave4.jpg|thumb|left|300px|Fig. 3. Plunging wave. Photo credit Kernowfile Flickr Creative Commons.]]
 +
|}
 +
 
 +
 
 +
 
 +
==Wave breaking formulas==
 +
 
 +
Analytical formulas for the breaker index have been derived from laboratory experiments, field observations and numerical simulation models. Reviews of the various methods and formulas have been given by Camenen and Larson (2007<ref name=CL>Camenen, B. and Larson, M. 2007. Predictive Formulas for Breaker Depth Index and Breaker Type. Journal of Coastal Research 23: 1028–1041</ref>) and Robertson et al. (2013<ref>Robertson, B., Hall, K., Zytner, R. and Nistor, I. 2013. Breaking waves: review of characteristic relationships. Coastal Engineering Journal 55(1), 1350002</ref>). Most formulas involve the steepness <math>H_0/L_0</math> of the incident wave and the average seabed slope <math>m</math> and reproduce the qualitative tendencies described above, i.e. <math>\gamma_b</math> is an increasing function of the wave steepness and the seabed slope.
 +
 
 +
 
 +
[[File:BreakerIndex.jpg|thumb|right|300px|Fig. 4. The area covered by experimental data for the breaker index as a function of the surf similarity parameter. The density of the experimental data is represented by varying degrees of redness. The black line is the Battjes formula Eq. 4.  Figure adapted from Ostendorf and Madsen (1979<ref name=OM>Ostendorf, D. and Madsen, O.  1979. An Analysis of Longshore Current and Associated Sediment Transport in the Surf Zone. Boston: Massachusetts Institute of Technology, Department of Civil Engineering Technical Report 241, 169p.</ref>).]]
 +
 
 +
A theoretical limit for <math>\gamma_b</math> was established by Miche (1944<ref>Miche, R. 1944. Mouvements ondulatoires de l’océan pour une eau profonde constante et décroissante. Annales des Ponts et Chaussées 114: 369-406</ref>), who analysed the wave motion over a horizontal seabed to determine the maximum crest steepness of the incident wave before collapse (occurring when the surface fluid velocity equals the wave crest celerity). The result is
 +
 
 +
<math>\gamma_b < \Large\frac{0.88}{k_b \, h_b}\normalsize \tanh(k_b\, h_b) , \quad k_b=2 \pi / L_b , \qquad (2)</math>
 +
 
 +
where <math>L_b</math> is the wavelength at the breakpoint. This relation was empirically refined by Ostendorf and Madsen (1979<ref name=OM>Ostendorf, D. and Madsen, O.  1979. An Analysis of Longshore Current and Associated Sediment Transport in the Surf Zone. Boston: Massachusetts Institute of Technology, Department of Civil Engineering Technical Report 241, 169p.</ref>) who proposed the formula
 +
 
 +
<math>\gamma_b =  \Large\frac{0.88}{k_b h_b}\normalsize \tanh(p \, k_b \, h_b) , \quad p = 0.8 + 5 \; min(m, 0.1)  . \qquad (3)</math>
 +
 
 +
Experiments by Battjes (1974<ref name=B>Battjes, J.A. 1974. Surf similarity. Proceedings 14th International Conference on Coastal Engineering, pp. 466–480</ref>) suggested that the dependence of <math>\gamma_b</math> on wave steepness and bed slope could be represented by the [[surf similarity parameter]] <math>\xi</math> at the breakpoint,
  
 
<math>\gamma_b = 1.06 + 0.14 \ln \xi_b . \qquad (4)</math>  
 
<math>\gamma_b = 1.06 + 0.14 \ln \xi_b . \qquad (4)</math>  
  
Many other more elaborate formulas have been proposed. For example, the formulas by Ostendorf and Madsen (1979<ref name=O>Ostendorf, D. and Madsen, O. 1979. An Analysis of Longshore Current and Associated Sediment Transport in the Surf Zone. Boston: Massachusetts Institute of Technology, Department of Civil Engineering Technical Report 241, 169p.</ref>)
+
Another frequently used formula was established by Goda (2010<ref name=G10>Goda, Y. 2010. Reanalysis of regular and random breaking wave statistics. Coastal Engineering Journal 52: 71-106</ref>),
 +
 
 +
<math>\gamma_b = 0.17 \large\frac{L_0}{h_b}\normalsize \Big[1 - \exp[-\large\frac{3 \pi h_b}{2L_0}\normalsize (1+11m^{4/3})] \Big] .  \qquad (5)</math>
 +
 
 +
The formulas (3) and (5) both require an iterative solution because of the dependence on <math>h_b</math>.
 +
 
 +
 
 +
The results of experiments designed to determine the breaker index yield widely scattered results, even for experiments with similar bed slope and incident waves. Fig. 4 gives an impression of the range covered by experimental data when the breaker index is plotted as a function of the surf similarity parameter. Although the scatter can be partly explained by the dependence of the breaker index on other parameters besides <math>\xi</math>, it also suggests that the breaking process should be considered a stochastic process that is very sensitive to small factors that are not accounted for in the gross characteristics of the experimental setup<ref name=G10/>.
 +
<br clear=all>
 +
 
 +
==Irregular waves==
 +
 
 +
The formulas (2-5) were derived for regular monochromatic waves. However, incident wave trains in nature usually have an irregular, random character. An irregular wave field can be characterized by a wave spectrum with significant wave height <math>H_s</math> (see [[Statistical description of wave parameters]]). Since individual waves in a wave train have different breakpoints in the surf zone, the breaker index concept cannot be applied as such for <math>H_s</math>. This problem can be addressed in two ways. The first way is to consider the breaking of individual waves in an incident wave train, represented by a wave-by-wave breaker index. The second way is to consider the location in the surf zone where <math>H_s</math> starts to decay (the intersection of the shoaling zone where <math>H_s</math> increases and the surf zone where <math>H_s</math> decreases<ref>Kamphuis, J. W. 1991. Incipient wave breaking. Coastal Engineering 15: 185–203</ref>); this location is called the incipient breakpoint, indicated by the subscript <math>ib</math>. The wave-by-wave breaker index is similar the breaker index for regular waves<ref name=X>Xu, J., Liu, S., Li, J. and Jia, W. 2020. Experimental study of breaker index of normal and oblique incident unidirectional and multidirectional irregular waves on slope. Ocean Engineering 213, 107792</ref>. The incipient breaker index <math>\gamma_{ib}</math> is significantly smaller than the breaker index for regular waves, but the dependence on wave steepness and bed slope is similar. The empirical formulas for regular waves can be applied for irregular waves if the value of <math>\gamma_b</math> is reduced by about 30% <ref name=OM/><ref name=G10/><ref name=X/>.
 +
 
 +
==Wave decay in the surf zone==
 +
 
 +
After initial breaking the wave propagates shoreward, often as a bore, while decreasing in height. The decay of the wave height <math>H(x)</math> can be derived from the equation:
 +
 
 +
<math> gradient \; of \; wave \; energy \; flux = wave \; energy \; dissipation, \qquad (6)</math>
 +
 
 +
in which for both terms explicit expressions as function <math>H(x)</math> have to be substituted. The wave energy dissipation depends on several factors, including the proportion of broken waves in the surf zone. If all waves are broken, the surf zone is said to be saturated, but for an irregular (random) wave field this is generally not the case. Explicit expressions for the terms in Eq. (6) have been proposed by Battjes and Janssen (1978<ref name=BJ>Battjes, J.A. and Janssen, J.P.F.M. 1978. Energy loss and set-up due to breaking of random waves. Proc. 16th Int. Conf. on Coastal Engineering. ASCE, New York, pp. 570–587</ref>) and Baldock et al. (1998<ref name=B98>Baldock, T.E., Holmes, P., Bunker, S. and Van Weert, P. 1998. Crossshore hydrodynamics within an unsaturated surf zone. Coast. Eng. 34: 173–196</ref>), but they rely on many assumptions (see the Appendix).
 +
 
 +
Many laboratory and field investigations have therefore been conducted to establish empirical relationships for the dependence of wave height on water depth in the surf zone,
 +
 
 +
<math>H(x)/h(x) = \gamma . \qquad (7)</math>.  
 +
   
 +
The simplest relationship is a constant ratio of wave height and water depth in the surf zone,
 +
 
 +
<math>\gamma (x) = \gamma_b  . \qquad (8)</math>
 +
 
 +
Flume experiments suggest that this relationship is satisfied to some extent for rather steep shoreface profiles (<math>m \approx 1/30</math>), but not for more gently sloping shorefaces<ref>Dally, W. R., Dean, R. G. and Dalrymple, R. A. 1985. A model for breaker decay on beaches. Proc. 19th Int. Conf. on Coastal Engineering, Vol. 1, pp. 82–98</ref>. For more gently sloping profiles the ratio <math>H/h</math> displays a notable increase at small water depths<ref>Raubenheimer, B., Guza, R. T. and Elgar, S. 1996. Wave transformation across the inner surf zone. J. Geophys. Res. 101(C11): 25,589–25,597</ref><ref name=G10/><ref name=X/><ref>Power, H. E., Hughes, M. G., Aagaard, T. and Baldock, T. E. 2010. Nearshore wave height variation in unsaturated surf. J. Geophys. Res. 115, C08030</ref>. In shallow water, regeneration of waves can take place, which propagate shoreward with less energy loss than wave bores. For irregular waves, the increase of <math>H_s/h</math> is also due to a large percentage of shoaling waves (still increasing in height) after incipient breaking, that gradually decreases as the wave train travels shoreward<ref name=G10/><ref name=X/>.
 +
 
 +
The use of Eq. (7) can be justified if the surf zone depth decreases monotonically towards the shoreline, with an approximately constant slope <math>m</math>. In other situations there is no reason to assume that the wave height is related to the local depth<ref name=BJ/>. For example, for a wave propagating across a surf zone with [[nearshore sandbars]], the wave height at the trough between sandbars depends on the crest height of the previous sandbar.
 +
 
 +
 
 +
 
 +
==Appendix==
 +
 
 +
Wave energy dissipation in the surf zone can be described by Eq. (6),
 +
 
 +
<math>\Large\frac{dP}{dx}\normalsize + D =0 , \qquad (A1)</math>
 +
 
 +
where the <math>x</math>-axis points towards the local shoreline. In linear wave theory the wave energy flux <math>P</math>  can be approximated by
 +
 
 +
<math>P= ¼ \rho g H_s^2 c_g \cos \overline \theta  , \qquad (A2)</math>
 +
 
 +
where <math>g</math> is the gravitational acceleration, <math>\rho</math> the seawater density, <math>\overline \theta</math> the mean incident wave direction and <math>c_g</math> the wave group velocity, which in shallow water is about equal to the wave propagation speed <math>c \approx \sqrt{gh}</math>.  
 +
 
 +
An expression for the wave energy dissipation <math>D</math> was derived by Battjes and Janssen (1978<ref name=BJ/>) based on analogy with the power dissipated by a bore
  
<math>\gamma_b =  0.14 \large\frac{L_b}{h_b} \normalsize \tanh(p \, k_b \, h_b) , \quad k=2 \pi / L ,  \; p = 0.8 + 5 \; min(m, 0.1) ,  \qquad (5)</math>  
+
<math>D \approx ¼ Q_b \rho g \Large\frac{H^3}{Th}\normalsize ,  \qquad (A3)</math>
  
and by Goda (2010<ref name=G>Goda, Y. 2010. Reanalysis of regular and random breaking wave statistics. Coastal Engineering Journal 52: 71-106</ref>)
+
where <math>T</math> is the peak wave period, <math>H</math> the wave height and <math>Q_b</math> the fraction of broken waves. Baldock et al. (1998<ref name=B98/>), by assuming a Rayleigh distribution of the incident waves (see [[Statistical description of wave parameters]]), derived for <math>D</math> the approximate expression
  
<math>\gamma_b = 0.17 \large\frac{L}{h_b}\normalsize [1-exp(-\large\frac{3 \pi h_b}{2L}\normalsize (1+11m^{4/3}))] ,  \qquad (6)</math>
+
<math>D = \Large\frac{\rho g}{4 T}\normalsize \; (H_b^2 + ½ H_s^2 ) \; \exp[-2(H_b/H_s)^2] \, . \qquad (A4)</math>
  
both requiring an iterative solution because of the dependence on <math>h_b</math>.
+
An expression of the breaker wave height <math>H_b</math> was derived by Battjes and Janssen (1978<ref name=BJ/>) based on Miche's criterion,
  
The results of experiments designed to determine the breaker index yield widely scattered results, even for experiments with similar bed slope and incident waves. Fig. 1 gives an impression of the range covered by experimental data when the breaker index is plotted as a function of the surf similarity parameter.   
+
<math>H_b=\Large\frac{0.88}{k}\normalsize \tanh(\large\frac{\gamma}{0.88}\normalsize kh) , \qquad (A5)</math>
<br clear=all>
 
 
 
  
The formulas (4, 5, 6) were derived for regular monochromatic waves. Analysis of experiments <ref name=O/><ref name=G/> indicates that these formulas can also be applied for irregular waves with <math>H</math> replaced by <math>H_s</math> and with <math>\gamma</math> multiplied by 0.7.
+
where <math>k=2 \pi / (cT)</math> is the wave number and <math>\gamma</math> is a breaker index (<math>H_b \approx \gamma h</math> in shallow water). By solving Eq. (A1) and comparing the result with data from field surveys on beaches with nearshore sandbars, Ruessink et al. (2003<ref>Ruessink, B.G., Walstra, D.J.R. and Southgate, H.N. 2003. Calibration and verification of a parametric wave model on barred beaches. Coastal Eng. 48: 139–149</ref>) found for the breaker index <math>\gamma</math> defined by equation (A5) the empirical relationship
  
For gently sloping shorefaces, typically <math>m < 0.02</math>, it may be assumed that wave breaking is saturated throughout the surf zone, i.e. Eq. (1) approximately holds with <math>\gamma (x) = \gamma_b \; </math><ref>LeMehaute, B . 1962. On non-saturated breakers and the wave run-up. Proceedings of the 8th International Conference Coastal Engineering, pp. 77-92, ASCE</ref>. For steeper slopes the ratio <math>H(x)/h(x)</math> increases towards the shoreline<ref name=G/>.
+
<math>\gamma \approx  0.29+0.76 \, kh  . \qquad (A6)</math>
 +
 
 +
Zhang et al. (2021<ref>Zhang, C., Li, Y., Cai, Y., Shi, J., Zheng, J., Cai, F. and Qi, H. 2021. Parameterization of nearshore wave breaker index. Coast Eng. 168, 103914</ref>) and Myrhaug et al. (2022<ref>Myrhaug, D., Wang, H. and Holmedal, L.E. 2022. Discussion/comments of “Parameterization of nearshore wave breaker index” by Chi Zhang, Yuan Li, Yu Cai, Jian Shi, Jinhai Zheng, Feng Cai and Hong Shuai. Coastal Engineering 173, 104042</ref>) provided evidence that the dependence of the breaker index <math>\gamma</math> on <math>kh</math> is in fact a function of the offshore wave steepness.
 +
 
  
  
Line 49: Line 119:
 
:[[Shoreface profile]]
 
:[[Shoreface profile]]
 
:[[Shallow-water wave theory]]
 
:[[Shallow-water wave theory]]
 +
:[[Wave transformation]]
 
:[[Surf similarity parameter]]
 
:[[Surf similarity parameter]]
 
:[[Waves]]
 
:[[Waves]]
Line 57: Line 128:
 
==References==
 
==References==
 
<references/>
 
<references/>
 +
 +
 +
 +
{{author
 +
|AuthorID=120
 +
|AuthorFullName=Job Dronkers
 +
|AuthorName=Dronkers J}}
 +
 +
 +
[[Category:Physical coastal and marine processes]]

Revision as of 08:33, 22 May 2022


The process of wave breaking plays a crucial role in the morphodynamics of the shoreface and the beach, such as accretion or erosion trends and the development and dynamics of nearshore bars and rip currents. The breaker index concept has been introduced to estimate the location on the shoreface where wave breaking occurs based on empirical formulas, without the need for a detailed simulation of the breaker process.

A practical wave breaking criterion

According to numerical simulations[1] and flume experiments[2], wave breaking starts on the shoreface when the horizontal surface water particle velocity at the wave crest comes close (within about 15%) to the wave crest celerity of the tallest wave in the incident wave train. However, it is not easy to derive from this criterion at which location on the shoreface wave breaking will occur. Even if the characteristics of the incident waves and the shoreface profile are well known, one has to face the problem that wave breaking is highly sensitive to fine-scale processes. It cannot be easily modelled, although significant progress has been achieved with numerical techniques such as Computational Fluid Dynamics (CFD)[3] or Smoothed Particle Hydrodynamics (SPH)[4]. For this reason, many investigations have been conducted to establish empirical rules for wave breaking from which the location of wave breaking can be easily derived. These empirical rules refer to the so-called breaker index. The breaker index relates the wave height at breaking, [math]H_b[/math], to the water depth [math]h_b[/math] at the breakpoint; it is defined as

[math]\gamma_b=H_b/h_b . \qquad (1)[/math]

Important parameters for the wave breaking process are the steepness [math]H_0/L_0[/math] of the incident waves ([math]H_0, L_0[/math] are the offshore wave height and wavelength) and the average shoreface slope [math]m[/math]. Empirical rules establish a relationship between the breaker index and these parameters.

Wave breaking process

Fig. 1. Schematic of wave evolution towards breaking, for swell waves (left panel) and sea waves (right panel).

The process of wave breaking starts on the shoaling zone where incident waves become progressively skewed by interaction with the seabed: the wave crest becomes sharper and higher, while the wave trough widens and flattens (Fig. 1). Incident waves with short wavelengths (locally generated 'sea') are steeper than incident waves with longer wavelengths (remotely generated swell); the former category of waves will break earlier, thus at greater depth [math]h_b[/math] than the latter category[3][5], implying a larger value of [math]\gamma_b[/math]. These short sea waves are also less skewed before breaking; they break when the crest becomes unstable and flows down the front face of the wave, producing spilling breaker bores that surf towards the shoreline (Fig. 2). This breaking mode is also favoured on gently sloping shorefaces when the surf similarity parameter [math]\xi=m/\sqrt{H_0/L_0}[/math] has values typically smaller than about 0.4. Incident waves with longer wavelengths will break later, thus at smaller depth [math]h_b[/math]. While shoaling, they are strongly skewed and break when the crest curls over the front face and falls into the base of the wave, producing a so-called plunging breaker[6] (Fig. 3). This breaking mode is also favoured on steep sloping shorefaces when the surf similarity parameter has values typically larger than 0.4.


Fig. 2. Spilling waves. Photo credit Andrew Dawley Flickr Creative Commons.
Fig. 3. Plunging wave. Photo credit Kernowfile Flickr Creative Commons.


Wave breaking formulas

Analytical formulas for the breaker index have been derived from laboratory experiments, field observations and numerical simulation models. Reviews of the various methods and formulas have been given by Camenen and Larson (2007[6]) and Robertson et al. (2013[7]). Most formulas involve the steepness [math]H_0/L_0[/math] of the incident wave and the average seabed slope [math]m[/math] and reproduce the qualitative tendencies described above, i.e. [math]\gamma_b[/math] is an increasing function of the wave steepness and the seabed slope.


Fig. 4. The area covered by experimental data for the breaker index as a function of the surf similarity parameter. The density of the experimental data is represented by varying degrees of redness. The black line is the Battjes formula Eq. 4. Figure adapted from Ostendorf and Madsen (1979[8]).

A theoretical limit for [math]\gamma_b[/math] was established by Miche (1944[9]), who analysed the wave motion over a horizontal seabed to determine the maximum crest steepness of the incident wave before collapse (occurring when the surface fluid velocity equals the wave crest celerity). The result is

[math]\gamma_b \lt \Large\frac{0.88}{k_b \, h_b}\normalsize \tanh(k_b\, h_b) , \quad k_b=2 \pi / L_b , \qquad (2)[/math]

where [math]L_b[/math] is the wavelength at the breakpoint. This relation was empirically refined by Ostendorf and Madsen (1979[8]) who proposed the formula

[math]\gamma_b = \Large\frac{0.88}{k_b h_b}\normalsize \tanh(p \, k_b \, h_b) , \quad p = 0.8 + 5 \; min(m, 0.1) . \qquad (3)[/math]

Experiments by Battjes (1974[10]) suggested that the dependence of [math]\gamma_b[/math] on wave steepness and bed slope could be represented by the surf similarity parameter [math]\xi[/math] at the breakpoint,

[math]\gamma_b = 1.06 + 0.14 \ln \xi_b . \qquad (4)[/math]

Another frequently used formula was established by Goda (2010[5]),

[math]\gamma_b = 0.17 \large\frac{L_0}{h_b}\normalsize \Big[1 - \exp[-\large\frac{3 \pi h_b}{2L_0}\normalsize (1+11m^{4/3})] \Big] . \qquad (5)[/math]

The formulas (3) and (5) both require an iterative solution because of the dependence on [math]h_b[/math].


The results of experiments designed to determine the breaker index yield widely scattered results, even for experiments with similar bed slope and incident waves. Fig. 4 gives an impression of the range covered by experimental data when the breaker index is plotted as a function of the surf similarity parameter. Although the scatter can be partly explained by the dependence of the breaker index on other parameters besides [math]\xi[/math], it also suggests that the breaking process should be considered a stochastic process that is very sensitive to small factors that are not accounted for in the gross characteristics of the experimental setup[5].

Irregular waves

The formulas (2-5) were derived for regular monochromatic waves. However, incident wave trains in nature usually have an irregular, random character. An irregular wave field can be characterized by a wave spectrum with significant wave height [math]H_s[/math] (see Statistical description of wave parameters). Since individual waves in a wave train have different breakpoints in the surf zone, the breaker index concept cannot be applied as such for [math]H_s[/math]. This problem can be addressed in two ways. The first way is to consider the breaking of individual waves in an incident wave train, represented by a wave-by-wave breaker index. The second way is to consider the location in the surf zone where [math]H_s[/math] starts to decay (the intersection of the shoaling zone where [math]H_s[/math] increases and the surf zone where [math]H_s[/math] decreases[11]); this location is called the incipient breakpoint, indicated by the subscript [math]ib[/math]. The wave-by-wave breaker index is similar the breaker index for regular waves[12]. The incipient breaker index [math]\gamma_{ib}[/math] is significantly smaller than the breaker index for regular waves, but the dependence on wave steepness and bed slope is similar. The empirical formulas for regular waves can be applied for irregular waves if the value of [math]\gamma_b[/math] is reduced by about 30% [8][5][12].

Wave decay in the surf zone

After initial breaking the wave propagates shoreward, often as a bore, while decreasing in height. The decay of the wave height [math]H(x)[/math] can be derived from the equation:

[math] gradient \; of \; wave \; energy \; flux = wave \; energy \; dissipation, \qquad (6)[/math]

in which for both terms explicit expressions as function [math]H(x)[/math] have to be substituted. The wave energy dissipation depends on several factors, including the proportion of broken waves in the surf zone. If all waves are broken, the surf zone is said to be saturated, but for an irregular (random) wave field this is generally not the case. Explicit expressions for the terms in Eq. (6) have been proposed by Battjes and Janssen (1978[13]) and Baldock et al. (1998[14]), but they rely on many assumptions (see the Appendix).

Many laboratory and field investigations have therefore been conducted to establish empirical relationships for the dependence of wave height on water depth in the surf zone,

[math]H(x)/h(x) = \gamma . \qquad (7)[/math].

The simplest relationship is a constant ratio of wave height and water depth in the surf zone,

[math]\gamma (x) = \gamma_b . \qquad (8)[/math]

Flume experiments suggest that this relationship is satisfied to some extent for rather steep shoreface profiles ([math]m \approx 1/30[/math]), but not for more gently sloping shorefaces[15]. For more gently sloping profiles the ratio [math]H/h[/math] displays a notable increase at small water depths[16][5][12][17]. In shallow water, regeneration of waves can take place, which propagate shoreward with less energy loss than wave bores. For irregular waves, the increase of [math]H_s/h[/math] is also due to a large percentage of shoaling waves (still increasing in height) after incipient breaking, that gradually decreases as the wave train travels shoreward[5][12].

The use of Eq. (7) can be justified if the surf zone depth decreases monotonically towards the shoreline, with an approximately constant slope [math]m[/math]. In other situations there is no reason to assume that the wave height is related to the local depth[13]. For example, for a wave propagating across a surf zone with nearshore sandbars, the wave height at the trough between sandbars depends on the crest height of the previous sandbar.


Appendix

Wave energy dissipation in the surf zone can be described by Eq. (6),

[math]\Large\frac{dP}{dx}\normalsize + D =0 , \qquad (A1)[/math]

where the [math]x[/math]-axis points towards the local shoreline. In linear wave theory the wave energy flux [math]P[/math] can be approximated by

[math]P= ¼ \rho g H_s^2 c_g \cos \overline \theta , \qquad (A2)[/math]

where [math]g[/math] is the gravitational acceleration, [math]\rho[/math] the seawater density, [math]\overline \theta[/math] the mean incident wave direction and [math]c_g[/math] the wave group velocity, which in shallow water is about equal to the wave propagation speed [math]c \approx \sqrt{gh}[/math].

An expression for the wave energy dissipation [math]D[/math] was derived by Battjes and Janssen (1978[13]) based on analogy with the power dissipated by a bore

[math]D \approx ¼ Q_b \rho g \Large\frac{H^3}{Th}\normalsize , \qquad (A3)[/math]

where [math]T[/math] is the peak wave period, [math]H[/math] the wave height and [math]Q_b[/math] the fraction of broken waves. Baldock et al. (1998[14]), by assuming a Rayleigh distribution of the incident waves (see Statistical description of wave parameters), derived for [math]D[/math] the approximate expression

[math]D = \Large\frac{\rho g}{4 T}\normalsize \; (H_b^2 + ½ H_s^2 ) \; \exp[-2(H_b/H_s)^2] \, . \qquad (A4)[/math]

An expression of the breaker wave height [math]H_b[/math] was derived by Battjes and Janssen (1978[13]) based on Miche's criterion,

[math]H_b=\Large\frac{0.88}{k}\normalsize \tanh(\large\frac{\gamma}{0.88}\normalsize kh) , \qquad (A5)[/math]

where [math]k=2 \pi / (cT)[/math] is the wave number and [math]\gamma[/math] is a breaker index ([math]H_b \approx \gamma h[/math] in shallow water). By solving Eq. (A1) and comparing the result with data from field surveys on beaches with nearshore sandbars, Ruessink et al. (2003[18]) found for the breaker index [math]\gamma[/math] defined by equation (A5) the empirical relationship

[math]\gamma \approx 0.29+0.76 \, kh . \qquad (A6)[/math]

Zhang et al. (2021[19]) and Myrhaug et al. (2022[20]) provided evidence that the dependence of the breaker index [math]\gamma[/math] on [math]kh[/math] is in fact a function of the offshore wave steepness.


Related articles

Shoreface profile
Shallow-water wave theory
Wave transformation
Surf similarity parameter
Waves
Wave set-up
Littoral drift and shoreline modelling


References

  1. Barthelemy, X., Banner, M., Peirson, W., Fedele, F., Allis, M. and Dias, F. 2018. On a unified breaking onset threshold for gravity waves in deep and intermediate depth water. Journal of Fluid Mechanics 841: 463-488
  2. Saket, A., Peirson, W., Banner, M., Barthelemy, X. and Allis, M. 2017. Wave breaking onset of two-dimensional deep-water wave groups in the presence and absence of wind. J. Fluid Mech. 811: 642–658
  3. 3.0 3.1 Aggarwal, A., Chella, M.A., Bihs, H. and Myrhaug, D. 2020. Properties of breaking irregular waves over slopes. Ocean Engineering 216, 108098
  4. Lowe, R.J., Buckley, M.L., Altomare, C., Rijnsdorp, D.P., Yao, Y., Suzuki, T. and Bricker, J.D. 2019 Numerical simulations of surf zone wave dynamics using Smoothed Particle Hydrodynamics. Ocean Modelling 144, 101481
  5. 5.0 5.1 5.2 5.3 5.4 5.5 Goda, Y. 2010. Reanalysis of regular and random breaking wave statistics. Coastal Engineering Journal 52: 71-106
  6. 6.0 6.1 Camenen, B. and Larson, M. 2007. Predictive Formulas for Breaker Depth Index and Breaker Type. Journal of Coastal Research 23: 1028–1041
  7. Robertson, B., Hall, K., Zytner, R. and Nistor, I. 2013. Breaking waves: review of characteristic relationships. Coastal Engineering Journal 55(1), 1350002
  8. 8.0 8.1 8.2 Ostendorf, D. and Madsen, O. 1979. An Analysis of Longshore Current and Associated Sediment Transport in the Surf Zone. Boston: Massachusetts Institute of Technology, Department of Civil Engineering Technical Report 241, 169p.
  9. Miche, R. 1944. Mouvements ondulatoires de l’océan pour une eau profonde constante et décroissante. Annales des Ponts et Chaussées 114: 369-406
  10. Battjes, J.A. 1974. Surf similarity. Proceedings 14th International Conference on Coastal Engineering, pp. 466–480
  11. Kamphuis, J. W. 1991. Incipient wave breaking. Coastal Engineering 15: 185–203
  12. 12.0 12.1 12.2 12.3 Xu, J., Liu, S., Li, J. and Jia, W. 2020. Experimental study of breaker index of normal and oblique incident unidirectional and multidirectional irregular waves on slope. Ocean Engineering 213, 107792
  13. 13.0 13.1 13.2 13.3 Battjes, J.A. and Janssen, J.P.F.M. 1978. Energy loss and set-up due to breaking of random waves. Proc. 16th Int. Conf. on Coastal Engineering. ASCE, New York, pp. 570–587
  14. 14.0 14.1 Baldock, T.E., Holmes, P., Bunker, S. and Van Weert, P. 1998. Crossshore hydrodynamics within an unsaturated surf zone. Coast. Eng. 34: 173–196
  15. Dally, W. R., Dean, R. G. and Dalrymple, R. A. 1985. A model for breaker decay on beaches. Proc. 19th Int. Conf. on Coastal Engineering, Vol. 1, pp. 82–98
  16. Raubenheimer, B., Guza, R. T. and Elgar, S. 1996. Wave transformation across the inner surf zone. J. Geophys. Res. 101(C11): 25,589–25,597
  17. Power, H. E., Hughes, M. G., Aagaard, T. and Baldock, T. E. 2010. Nearshore wave height variation in unsaturated surf. J. Geophys. Res. 115, C08030
  18. Ruessink, B.G., Walstra, D.J.R. and Southgate, H.N. 2003. Calibration and verification of a parametric wave model on barred beaches. Coastal Eng. 48: 139–149
  19. Zhang, C., Li, Y., Cai, Y., Shi, J., Zheng, J., Cai, F. and Qi, H. 2021. Parameterization of nearshore wave breaker index. Coast Eng. 168, 103914
  20. Myrhaug, D., Wang, H. and Holmedal, L.E. 2022. Discussion/comments of “Parameterization of nearshore wave breaker index” by Chi Zhang, Yuan Li, Yu Cai, Jian Shi, Jinhai Zheng, Feng Cai and Hong Shuai. Coastal Engineering 173, 104042


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 (2022): Breaker index. Available from http://www.coastalwiki.org/wiki/Breaker_index [accessed on 29-03-2024]