Open Access
How to translate text using browser tools
1 July 2007 Predictive Formulas for Breaker Depth Index and Breaker Type
Benoît Camenen, Magnus Larson
Author Affiliations +

This study was undertaken to establish the most reliable breaker depth index formula that will yield satisfactory predictions for a wide range of hydraulic and beach conditions. The applicability of six existing formulas for computing the breaker depth index was examined with the use of a large amount of published laboratory data (524 cases collected from 22 sources). The behavior of these formulas was studied with respect to the main governing parameters, which are the beach slope and the deep-water wave steepness. It was found that most formulas show quite good predictions for cases including gentle slopes (0.01 < m ≤ 0.07). However, the predictions are typically not satisfactory for breaking waves on steep slopes (m > 0.1), and the formulas do not always present a physically correct behavior with respect to the two main parameters. A new formula is proposed to predict the breaker depth index with the best possible accuracy. Also, a discussion is included on the application of the formula to random waves and on the relationship between the breaker depth index and the Irribaren number to distinguish between different breaker types.


To analyze or simulate coastal processes (e.g., to estimate surf zone hydrodynamics in studies on nearshore morphology), the breaking wave height is one of the essential quantities to determine because it characterizes the most dynamic phenomenon in the nearshore zone. Thus, to compute wave height transformation in the surf zone, it is necessary to determine the initiation of breaking and the breaker type. A spilling breaker creates intense turbulence mainly at the surface, whereas a plunging breaker generates a large vortex that penetrates below the water surface. (Ting and Kirby, 1995; Wang, 1998; Wang and Kraus, 1999). Galvin (1968) observed a relationship between the breaker type and the Irribaren number ξ, defined as,

where m is the slope of the beach and λ = H/L is the deep-water (offshore) wave steepness (H and L are, respectively, the deep-water wave height and wavelength).

A simple way to estimate the location of the breaker line is to compute the wave height-to-water depth ratio where the wave will break (breaker depth index γb), which seems, according to different authors (Galvin, 1972; Goda, 1970; Miche, 1944; Weggel, 1972), to be a function mainly of the offshore wave steepness and the mean beach slope. Rattanapitikon and Shibayama (2000) compared many existing formulas and showed that they differed significantly and induced large uncertainties when slopes deviated from the typical range (i.e., m < 0.02 or m > 0.1). Because these formulas are semiempirical and often based on a limited data set, their accuracy over large data ranges might not be so good.

In this paper, we compare six existing formulas for calculating γb, with an extensive data set. On the basis of this comparison, a new semiempirical formula is proposed involving the beach slope and the offshore wave steepness that displays improved behavior and agreement with the data. Laboratory data on breaking wave height employed in this paper have been compiled from various sources, as summarized in Table 1. Also, because the Irribaren number is often used to distinguish between different breaker types, a discussion is included on the relationship between the breaker depth index and the Irribaren number.



The breaker depth index is commonly used to define the wave height at breaking (Horikawa, 1988, pp. 79–88; Van Rijn, 1990, pp. 298–305).

where h is the water depth and subscript b denotes the point of incipient breaking.

A second parameter, the breaker height index Ωb = Hb/H is also widely used (Smith and Kraus, 1990). However, according to the present authors, it induces a greater uncertainty in the prediction of Hb. Indeed, for deep-water wave steepnesses up to 0.02, this parameter is less sensitive than Equation (2) because predictions by existing formulas are approximately between 1 and 1.1 (cf. Rattanapitikon, Vivattanasirisak, and Shibayama, 2003; Smith and Kraus,1990).

Many different formulas for γb= have been proposed on the basis of monochromatic wave experiments. Rattanapitikon and Shibayama (2000) and Smith and Kraus (1990) compared some of these formulas with data. It appears that the main governing factors for γb are, as for the breaker type, the bottom slope m and the deep-water steepness λ. According to the results of Rattanapitikon and Shibayama (2000), six formulas are interesting to compare.

Weggel (1972)


Battjes (1974)


Ostendorf and Madsen (1979)


Singamsetti and Wind (1980)


Smith and Kraus (1990)


Goda (1970) modified by Rattanapitikon and Shibayama (2000)


It should be noted that Equations (3), (5), and (8) require an iterative technique, whereas Equations (4), (6), and (7) allow for direct computation. The range of validity for these formulas is implicitly assumed to be the range of the data used for their calibration. Both the formulas by Singamsetti and Wind (1980) and Smith and Kraus (1990) were calibrated with their own data set only (cf. Table 1). The formulas by Battjes (1974), Ostendorf and Madsen (1979), and Weggel (1972) are mainly based on the Galvin (1968, 1969) and Iversen (1952) data set and thus have a similar range of validity (0.02 < m ≤ 0.2). Finally, the Goda (1970) formula modified by Rattanapitikon and Shibayama (2000) was calibrated with a similar data set to that compiled for this paper (i.e., 0.00 < m ≤ 0.5).

Comparison with Compiled Data Set

Table 2 summarizes the errors in the predictions of γb for the six studied formulas and for the newly developed formula discussed later, including error estimates for all the data as well as for the cases with a plane-sloping beach, a step beach, and a barred beach. The quantities P10 and P20 denote the percentage of values obtained with an error of less than 10% and 20%, respectively, whereas Erms is the root mean square relative error (%), defined as

where γb,pred is the predicted breaker depth index and γb,exp is the breaker depth index observed experimentally. The results obtained for a step beach (Horikawa and Kuo, 1966) and for a barred beach (Smith and Kraus, 1990) were separated out because these cases present a different situation than the classical plane-sloping beach. Experiments on a step beach (m = 0) induce significant reflection, which can influence the wave breaking process. In the same way, Smith and Kraus (1990) observed that the bar, depending on its shape, can induce a strong seaward return flow, which could influence wave breaking. Thus, error estimates are presented for all data, as well as separately for the data with regard to the plane-sloping beach, the step beach, and the barred beach.

The mean value and the standard deviation (SD) of the difference Δb = γb,pred − γb,exp is also included (computed for all the data and the data for which m > 0). The figures in Appendix A show the comparison between predicted and observed breaker depth index obtained for the six formulas with the use of all the experimental data.

The best results for the existing formulas and the use of all the data are obtained with the Battjes (1974), Goda (1970), and the modified Ostendorf and Madsen (1979) formulas. The latter formula also yields the best results for the barred beach data. Finally, the Weggel (1972) formula shows the best results for the step beach experiments, although it gives a constant value γb = 0.78 independent of the deep-water wave conditions. Some formulas can reach an agreement with the data of close to 80% for a permitted error of 20%, whereas none reaches 50% for a permitted error of 10%. This deviation can partly be explained by the errors that come from the experimental measurements since the breaker height is typically visually estimated. Also, the break point can be defined in several ways. Singamsetti and Wind (1980) listed seven possible definitions, which could induce additional errors. Closer inspection of the graphs (see Appendix A) allows for the following comments.

  • Weggel's (1972) formula generally overestimates the breaker depth index and produces considerable dispersion of the results. This overestimation could occur because Weggel intended his formula to yield an envelope to the observed data to obtain a high estimate on the breaker height appropriate for design.

  • Ostendorf and Madsen's (1979) formula gives good agreement with the data but seems to exhibit a discrete behavior (same computed value for different experimental data), and no value can be computed for a beach with a slope equal to zero (γb = 0, if m = 0).

  • Singamsetti and Wind's (1980) formula generally overestimates the breaker depth index and produces considerable dispersion of the results.

  • Battjes' (1974) formula yields overall good results, but no value can be computed for a beach with a slope equal to zero (γb = 0, if m = 0).

  • Smith and Kraus' (1990) formula often underestimates the breaker depth index.

  • The modified Goda (1970) formula gives less dispersion in the results but tends to overestimate for small values on the breaker depth index and to underestimate for large values.

Dependence of Breaker Depth Index on Deep-Water Wave Steepness and Beach Slope

With regard to the dependence of γb on the deep-water wave steepness, the behavior of each of the studied formulas is not entirely satisfactory. The graphs in Appendix B.1 present the results of comparing data and formulas for varying λ under fixed bed slope. Although it is quite difficult to observe any trend in the experimental data of Horikawa and Kuo (1966) (m = 0, m = 0.0125; see graphs in Appendix B.1a and 1b), because significant dispersion exists in the data, a decrease in γb with an increase in λ remains the principal tendency in the other data sets. This decrease becomes clearer and quite pronounced in the data from Iwagaki et al. (1974), Singamsetti and Wind (1980), and Smith and Kraus (1990) for mild slopes (see Appendix B.1c–f). Only the Smith and Kraus (1990) formula displays a correct behavior in accordance with this observation. The five other formulas seem insensitive for high values on λ. For larger values of m (see Appendix B.1g), this sensitivity is not so apparent, mainly because of the lack of data. Finally, for smaller values of λ, most of the data indicate a value of γb that is nearly independent of λ, with a decrease in λ (see Appendix B.1c–e with the data of Galvin [1968, 1969], Iversen [1952], Van Dorn [1978] and Walker [1974]). Both the Battjes (1974) and the Singamsetti and Wind (1980) formulas show unsatisfactory predictions for those specific cases. To summarize the influence of the deep-water wave steepness on the breaker depth index, the best observed behavior is displayed by the Smith and Kraus (1990) formula for mild slopes and by the modified Goda (1970) formula for slopes up to 10%.

The beach slope, m, has often been the only parameter used for calibration of semiempirical formulas, mainly on the basis of Miche's (1944) results. Thus, good agreement is observed for most of the formulas for mild slopes (i.e., γb displays an exponential dependence on m; see graphs in Appendix B.2). According to the data (except for the step beach experiment), the Smith and Kraus (1990) formula yields the best results for slopes below 0.1. However, for larger slopes, the tendency changes, and γb becomes a decreasing function of m (see graphs B.2 with Smith and Kraus' data) after having reached a maximum for m ≈ 0.2. Above this value, the Battjes (1974) and the Singamsetti and Wind (1980) formulas display incorrect behavior. The Ostendorf and Madsen (1979), Smith and Kraus (1990), and Weggel (1972) formulas, which all assume a fixed value for m > 0.1, give better agreement with data, but it is not satisfactory. Only the Goda (1970) formula modified by Rattanapitikon and Shi-bayama (2000) shows a correct behavior for beach slopes up to 0.2. For lower values on the deep-water wave steepness (λ < 0.03), it seems like the maximum γb is reached for a milder slope (m ≈ 0.1), and the observed decrease with larger slopes is smoother.

A New Formula For Breaker Depth Index

One aim of this study was to propose a new formula for γb that exhibits the best possible behavior with respect to the dependence on the mean slope m and the offshore wave steepness γ. The new formula is based on Miche's (1944) results similar to the Ostendorf and Madsen (1979) formula, but it employs the offshore wave steepness instead of the ratio hb/Lb to allow for a direct computation of the breaker depth index.

The Miche (1944) formula can be written as

where Lb is the wavelength at the break point. The function f(m, λ) depends on the bottom slope and was empirically introduced by Ostendorf and Madsen (1979) with f(m, λ) = 0.8 + 5 min(m, 0.1) (f(m, λ) = 1 in case of the Miche [1944] formula).

Assuming shoaling follows the first-order Stokes theory, the wavelength can be expressed as a function of its deep-water value L = gT2w/(2π) and the water depth h.


(if kh « 1, k = i1551-5036-23-4-1028-ex01.gif). Because hb = αbH3b, where αb >1 is the shoaling coefficient at the break point.

where Ab = i1551-5036-23-4-1028-ex02.gif is assumed to be approximately constant. It appears from Figure 1 that Ab ≈ 0.5 is a reasonable value over a wide range of steepness and slope values.

With Equation (12), Equation (10) can be rewritten in an explicit form.


Following the method proposed by Rattanapitikon and Shibayama (2000), an empirical function can be fit to the data to obtain the function f* (m, λ)


Table 3 presents a comparison between the results from the iterative method (Equation [10]) and the explicit method (Equation [13]) with a special fitting of f*(m, λ), where P20 is the percentage of data points in agreement within a factor of 1.2. If f* = 1, Equation (13) yields a formula quite close to Miche's (1944) formula. Assuming that f(m, λ) = f(m), Ostendorf and Madsen (1979) fit a straight line using data in the range of 0.1 < m < 0.1. Rattanapitikon and Shibayama (2000) fit a polynomial function of m using a similar data set as the one compiled in this study. The fourth column presents a fit with the use of the present data set and Equation (13). Similar results are found for the formulas. The differences appear to be mainly due to the different data sets used by the authors for the fitting (P20 and Erms not improved). It appears that employing Equation (13) does not reduce the predictive skill of the formula compared with the iterative relationships. More than 95% of the predicted data points are in agreement within a factor of 1.2. Some differences occur for the small values of γb predicted by the iterative formula, which explains the apparent dispersion of the results (Erms > 2).

These expressions for f (and f*) are only dependent on the mean slope, although a clear influence from the deep-water wave steepness can be observed in the graphs (Appendix C). Moreover, the polynomial function proposed by Rattanapitikon and Shibayama (2000) diverges rapidly for large slopes. Instead, a sinusoidal function of the beach slope is introduced.

where mmax is the beach slope for which f* reaches its maximum value, and A1, A2, and α are fitting parameters that are functions of λ

The graphs in Appendix C present f* as a function of m together with the experimental data for three different intervals of λ The best polynomial fit (second order) and the final relationship obtained by the authors are included. A maximum of f* is clearly observed and seems to be a linear function of λ. The following relationship is proposed.


In the same way, the coefficient α, which represents the “curvature” of the sinusoidal function (for mmmax and m >mmax), can be related to the deep-water wave steepness according to the formulas in Equation (17).


Finally, the best fit is obtained (Equations [18] and [19]).

Figure 2 presents the result when plotting predictions with the proposed formula against all data, showing an improvement compared with the previous relationships.

Errors obtained with the new formula are summarized in Table 2. It is interesting to see that only 15% of the results for all beach types are found to have an error greater than 20%, and less than 50% of the results for an error greater than 10%. Predictions for plane, barred, and step beaches are the best among the studied formulas, but, as explained in the previous section, experimental results for step beaches (experiment of Horikawa and Kuo, 1966) can imply an overestimation of γb because of possible reflection. However, it can be noted that (as for the Battjes [1974], Goda [1970], Ostenford and Madsen [1979], and Smith and Kraus [1990] formulas) the maximum predicted value by Equation (15) is 1.15, although some observed values reach 1.6. It is a result of the formula being fit to a large amount of scattered data. Finally, because of the errors and uncertainties in the experimental data, it seems that the quadratic error of less than 1.9% is difficult to achieve. One of the great improvements with the new formula is that a correct behavior with regard to beach slope and deep-water wave steepness is obtained for most of the cases (unlike previous formulas).

The distribution of the difference Δγb = γb,pred − γb,exp (cf. Figure 3) follows a Gaussian probability density function (pdf), which implies that errors correspond mainly to the experimental uncertainties or randomness in the breaking process itself. Because Equations (13) and (15) yield a prediction of the mean value of the observed experimental data, they cannot predict extreme observed values. Looking at the figures in Appendix B, it can be observed that, for similar conditions (same slope and offshore wave steepness), a large scatter appears depending on the experimental data. Thus, in an application of Equations (13) and (15), a random variation might be employed to achieve more realistic predictions of γb.

Application to Random Wave Model

Modeling random waves in the surface zone follow basically two different approaches:

  1. assume a pdf, valid everywhere, and use this function in the wave transformation calculations to obtain the statistical wave properties;

  2. assume a pdf in the offshore, represented by a number of individual wave components, compute the evolution of each monochromatic wave, and determine the local pdf by adding together the different components.

The calculation of the wave transformation is determined by applying an energy flux balance. For shore-parallel depth contours, it reads,

where E is the wave energy, Cg the group velocity, and D the time-averaged energy dissipation due to broken waves and frictional losses as the bed (x-axis points onshore). Assuming that the friction at the bottom is negligible and that the energy dissipation at each wave front corresponds to that in a bore or a hydraulic jump, the energy dissipation can be written as
where H and T are the wave height and period, respectively,and Qb is the proportion of broken waves.

In case of approach 1 described above, H = Hrms (root mean square wave height) and T = Tp (peak wave period). Many of the random wave models are based on the approach of Battjes and Janssen (1978). The wave heights in the surf zone are characterized by a Rayleigh distribution, which is assumed to be valid for both broken and unbroken waves. The usual criterion employed to truncate the distribution is that a wave is breaking when its height Hi exceeds some fraction of the water depth. So, to a fixed water depth h, there corresponds a critical wave height Hib over which all the waves are broken (see Figure 4). The quantity Qb is found by integrating the Rayleigh distribution over all waves for which Hi > Hib (Equation [22]).


The height Hib is generally evaluated from calibration of the random wave model with data. A common expression developed by Battjes and Stive (1985) is based on the Miche (1944) formula.

where k is the wave number and γ the breaker parameter. As a first approximation, γ = 0.88. Nairn (1990) modified the value proposed by Battjes and Stive (1985).

More recently, Ruessink, Walstra, and Southgate (2003) found a linear relationship with the product of the local wave number and the water depth.


As Ruessink, Walstra, and Southgate (2003) discussed, these relationships include not only the physics of the breakers but also encapsulate structural errors in the model, which might explain the different behavior of these calibration parameters compared with “measured” parameters. However, because the definition of Hib corresponds exactly to the breaker depth index for a monochromatic wave, Equations (13) and (15) can be used (i.e., Hib = γb h), which gives acceptable results.

In the case of the second approach H = Hi (individual wave height from the pdf) and a breaker depth index formula is directly applicable because waves are modeled individually. So, Qb = 0 if the individual wave is not broken Hi/h ≤ γb and Qb = 1 if the individual wave is broken.

An example is given in Figure 5 of the calculated evolution of random waves according to both approaches, and two different offshore wave conditions and beach slopes are shown. For the first approach, Hib was estimated by Equations (13) and (15) or Equation (23) with Equations (24) or (25). For the second approach, only Equations (13) and (15) were used.

In Figure 5, ξ = 0.125 for both cases; however in Figure 5a, λ is one-fourth of that in Figure 5b, and the beach slope m is half. In Figure 5a, it appears that the wave-by-wave approach of Equation (15) gives similar results to the parametric approach with a calibrated coefficient. On the other hand, in Figure 5b, the parametric model gives similar results independently of the formula used to calculate Hib. The wave-by-wave approach yields a smaller wave height inside the surf zone. Comparing the wave-by-wave and parametric approaches with Equation (15), it appears that the wave height in the surf zone is always larger with the parametric method. If a coefficient 0.7 (∼1/ i1551-5036-23-4-1028-ex03.gif) is added for the calculation of Hib with Equation (15) similar results are obtained for all conditions. Therefore, assuming that the wave-by-wave approach produces results closer to reality, Equation (15) can be used in a parametric model after introducing a coefficient of 0.7.


Breaker Type

Assuming that an accurate estimation of the breaker depth index can be made for all cases, it is interesting to see how this index is related to the breaker type. Galvin (1968, 1972) defined the following breaker types:

  • spilling breakers occur if the wave crest becomes unstable and flows down the front face of the wave producing a foamy water surface;

  • plunging breakers occur if the crest curls over the front face and falls into the base of the wave, resulting in a high splash;

  • collapsing breakers occur if the crest remains unbroken while the lower part of front face steepens and then falls, producing an irregular turbulent water surface;

  • surging breakers occur if the crest remains unbroken and the front face of the wave advances up the beach with minor breaking.

Breaker type is controlled by the bottom slope m and the deep-water steepness λ. Battjes (1974) used the Irribaren number (or surf similarity parameter) to describe breaker type on the basis of previous results of Galvin (1968, 1972):


Relationship Between Breaker Depth Index and Irribaren Number

Figure 6 shows how measured γb varies with ξ. It is clearly seen that γb is an increasing function of ξ. The relationship proposed by Battjes (1974) is a good approximation on the basis of comparison with all the data. However, marked scatter exists, which is the weakness of this formula, as was observed in the previous section.

In Figure 7, a sensitivity study of the new formula (Equation [15]) is made for γb against ξ employing various beach slopes m (Figure 7a) or various deep-water wave steepnesses λ (Figure 7b). Curves are limited by imposed constraints, i.e., the deep-water wave steepness cannot be greater than 0.14 and the maximum beach slope is fixed as the maximum value from the experimental data (m ≤ 0.5).

The most interesting result to emerge from these graphs is that they show two different types of behavior of the breaker depth index depending on the breaker type. Indeed, for ξ < ξ∞, cr ≈ 0.5 (spilling waves, cr is critical value), γb is mainly sensitive to the offshore wave steepness λ, and for a fixed ξ, γb is nearly independent of the slope m. On the contrary, for ξ < ξ∞, cr (plunging and surging waves), γb is mainly sensitive to the beach slope m, and for a fixed ξ, γb is nearly independent of λ. Moreover, the maximum values on the breaker depth index (see Figure 7b) seem to correspond to plunging breakers (according to values of λ, this maximum appears for 0.8 < ξ < 4). Thus, surging and collapsing breakers only occur for conditions corresponding to the decreasing part of the computed curves. The physical mechanisms controlling these types of breakers do not allow high values of γb, even if they were observed for lower values of λ (≈0.01) and steep beaches (m > 0.2). In the same way, a separation point between spilling and plunging waves could more or less be the inflection point of the increasing part of the computed curves (see Figure 7b), which explains the critical ξ from Galvin (1968, 1972) varying between 0.2 and 0.5, which Smith and Kraus (1990) also noticed.

Figure 8 presents the final result concerning breaker type prediction on the basis of the Irribaren number and the breaker depth index with the use of data in which complete information is available. The data from Okazaki and Sunamura (1989) were also added because they performed an interesting study on step beaches. Unfortunately, the experimental results for the breaker depth index were not presented in their article. Thus, γb has been estimated by Equation (15). The symbols used for this data set have been made smaller.

The lower curve (separating of the nonbreaking wave region) corresponds to Equation (15) with λ = 0.14. The upper curve (separating of the unstable wave region) corresponds to the maximum values obtained with this equation with a factor of 1.25. The factor 1.25 was employed because the new equation does not compute the maximum values observed from the experimental data. Finally, the three other curves are empirically determined to separate the spilling breaker, plunging breaker, collapsing breaker, and surging breaker regions following the ideas presented in the previous paragraph. Large uncertainties still exist, but they seem to be mainly due to the uncertainties in the available experimental data. This graph provides an improvement for breaker type prediction compared with Galvin's (1968, 1972) relationships (Equations [26]). The separation lines in Figure 8 are defined by the formulas in Equations (27).



More than 500 experimental cases from 22 published sources covering a wide range of beach slopes and wave conditions were used to czompare predictions from six existing breaker depth index formulas. A sensitivity study against beach slope and offshore wave steepness allowed for exploring the limits of these different formulas. The modified Goda (1970) formula, adding a bottom slope effect proposed by Rattanapitikon and Shibayama (2000), improves the behavior of the original formula substantially from slopes up to 0.1 but gives an overestimation for smaller slopes. On the basis of the sensitivity study, a new breaker depth index formula (Equation [13]) was proposed starting with Miche's (1944) expression and a correction factor that is a function of both beach slope and offshore wave steepness (Equation [15]).

This relationship presents better results when compared with all the data, but above all, it shows an improved behavior with respect to the two studied parameters. Some questions still remain for step beaches because reflection has a great influence. A tentative method was proposed for the application of the new breaker depth index formula in random wave models. With the use of a wave-by-wave model, the formula can be used as it stands (because it corresponds to the summation of the effects of a number of monochromatic waves), whereas a coefficient of 0.7 should be employed with a parametric model. Finally, with this new formula, a more accurate prediction of breaker type (compared with the Galvin [1968, 1972] equations) was obtained employing the Irribaren number and the breaker depth index.


The following symbols and subscripts are used in this paper.


  • A1, A2 = empirical coefficients

  • C = wave speed

  • Cg = group speed

  • D = wave energy dissipation

  • E = wave energy

  • h = water depth

  • H = wave height

  • Hib = critical wave height for the truncation of the wave probability distribution

  • L = wavelength

  • m = beach slope

  • Qb = proportion of broken waves

  • T, Tp = wave period and peak wave period

  • x = horizontal coordinate perpendicular to the beach

  • α = empirical coefficient

  • γb = breaker depth index

  • λ = wave steepness

  • Ωb = breaker height index

  • ξ = Irribaren number


  • b = denotes the point of incipient breaking

  • ∞ = denotes a deep-water value

  • cr = denotes a critical value

  • exp = denotes an experimental result

  • max = denotes a maximum value

  • pred = denotes a predicted result

  • rms = denotes a root mean square value



J. Battjes 1974. Surf similarity. In: Proceedings of the 14th International Conference on Coastal Engineering (Copenhagen, Den-mark, ASCE), pp. 1050–1061. Google Scholar


J. Battjes and J. Janssen . 1978. Energy loss and set-up due to breaking of random waves. In: Proceedings of the 16th International Conference on Coastal Engineering (Hamburg, Germany, ASCE), pp. 569–587. Google Scholar


E. Christensen and R. Deigaard . 2001. Large eddy simulation of breaking waves. Coastal Engineering 42 1:53–86. Google Scholar


C. Galvin 1968. Breaker travel and choice of design wave height. Journal of Geophysical Research 73:3651–3659. Google Scholar


C. Galvin 1969. Breaker type classification on three laboratory beaches. Journal of Waterways, Port, Coastal and Ocean Engineering 95:(WW2). 175–200. [formerly Waterways and Harbors Division]. Google Scholar


C. Galvin 1972. Wave breaking in shallow water. In: Meyer, R.E. (ed.), Waves on Beaches and Resulting Sediment Transport. New York: Academic Press, pp. 413–456. Google Scholar


Y. Goda 1970. A synthesis of breaker indices. Transactions of JSCE 2:227–230. Google Scholar


K. Horikawa 1988. Nearshore Dynamics and Coastal Processes, Theory, Measurement and Predictive Models. Tokyo, Japan: University of Tokyo Press. Google Scholar


K. Horikawa and C. Kuo . 1966. A study of wave transformation inside the surf zone. In: Proceedings of the 10th International Conference on Coastal Engineering (Tokyo, Japan, ASCE), pp. 217–233. Google Scholar


H. Iversen 1952. Laboratory Study of Breakers. Washington, DC: National Bureau Standards Circular Technical Report 52. Google Scholar


Y. Iwagaki, T. Sakai, K. Tsukioda, and N. Sawai . 1974. Relationship between vertical distribution of water particle velocity and type of breaker on beaches. Coastal Engineering Journal 17:51–58. Google Scholar


K. Maruyama, T. Sakakiyama, R. Kajima, S. Saito, and T. Shimizu . 1983. Experimental Study on Wave Height and Particle Velocity Near the Surf Zone Using a Large Wave Flume. Chiba, Japan: Central Research Institute Electric Power Industry (CRIEPI) Technical Report [in Japanese]. Google Scholar


R. Miche 1944. Mouvements ondulatoires de l'oce´an pour une eau profonde constante et de´croissante. Annales des Ponts et Chaus-se´es, (mai–juin) 25–28, (juillet–aou⁁t) 131–164, 270– 292, 369–406 [in French]. Google Scholar


M. Mizuguchi 1980. A heuristic model of wave height distribution in the surf zone. In: Proceedings of the 17th International Conference on Coastal Engineering (Sydney, Australia, ASCE), pp. 278–289. Google Scholar


K. Nadaoka and T. Kondoh . 1982. Laboratory measurements of velocity field structure in the surf zone by ldv. Coastal Engineering Journal 25:125–145. Google Scholar


R. Nairn 1990. Prediction of Cross-Shore Sediment Transport and Beach Profile Evolution. London, UK: Imperial College, Department of Civil Engineering, Doctoral thesis. Google Scholar


A. Okayasu, T. Shibayama, and N. Nimura . 1986. Velocity field under plunging waves. In: Proceedings of the 20th International Conference on Coastal Engineering (Taipei, Taiwan, ASCE), pp. 660–674. Google Scholar


S. Okazaki and T. Sunamura . 1989. Wave-tank study of breaker types on steeply inclined beaches. Annual Report of the Institute of Geoscience, University of Tsukuba 15:33–35. Google Scholar


D. Ostendorf and O. Madsen . 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. 169. p. Google Scholar


W. Rattanapitikon and T. Shibayama . 2000. Verification and mod-ification of breaker height formulas. Coastal Engineering Journal 42 4:389–406. Google Scholar


W. Rattanapitikon, T. Vivattanasirisak, and T. Shibayama . 2003. A proposal of new breaker height formula. Coastal Engineering Journal 45 1:29–48. Google Scholar


B. Ruessink, D. Walstra, and H. Southgate . 2003. Calibration and verification of a parametric wave model on barred beaches. Coastal Engineering 48:139–149. Google Scholar


H. Saeki and M. Sasaki . 1973. A study of the deformation of waves after breaking. In: Proceedings of the 20th Japanese Conference on Coastal Engineering (Nagoya, Japan, JSCE), pp. 559–564. [in Jap-anese]. Google Scholar


S. Singamsetti and H. Wind . 1980. Characteristics of Breaking and Shoaling Periodic Waves Normally Incident on to Plane Beaches of Constant Slope. Delft, The Netherlands: Delft University of Technology Technical Report M1371. Google Scholar


E. Smith and N. Kraus . 1990. Laboratory Study on Macro-Features of Wave Breaking Over Bars and Artificial Reefs. Vicksburg, Virginia: U.S. Army Corps of Engineers, Waterways Experiment Station Technical Report CERC-90-12. Google Scholar


M. Stive and J. Battjes . 1984. Velocity and pressure field in spilling breakers. In: Proceedings of the 19th International Conference on Coastal Engineering (Houston, Texas, ASCE), pp. 1420–1436. Google Scholar


M. Stive and H. Wind . 1986. Cross-shore mean flow in the surf zone. Coastal Engineering 10:151–158. Google Scholar


F. Ting 2002. Laboratory of wave and turbulence characteristics in narrow-band irregular breaking waves. Coastal Engineering 46:291–313. Google Scholar


F. Ting and J. Kirby . 1995. Dynamics of surf-zone turbulence in a strong plunging breaker. Coastal Engineering 34:177–204. Google Scholar


W. Van Dorn 1978. Breaking invariants in shoaling waves. Journal of Geophysical Research 83:C6. 2981–2988. Google Scholar


L. Van Rijn 1990. Principles of Fluid Flow and Surface Waves in River, Estuaries, Seas and Ocean. Amsterdam: Aqua Publications. Google Scholar


P. Visser 1982. The Proper Longshore Current in a Wave Basin. Delft, The Netherlands: Delft University of Technology Technical Report 82-1. Google Scholar


J. Walker 1974. Wave Transformations Over a Sloping Bottom and Over a Three Dimensional Shoal. Honolulu: University of Hawaii Technical Report 11. Google Scholar


P. Wang 1998. Longshore sediment flux in the water column and across the surf zone. Journal of Waterways, Port, Coastal and Ocean Engineering 124:108–117. Google Scholar


P. Wang and N. Kraus . 1999. Longshore sediment transport rate measured by short-term impoundment. Journal of Waterways, Port, Coastal and Ocean Engineering 125:118–126. Google Scholar


J. Weggel 1972. Maximum breaker height for design. In: Proceedings of the 13th International Conference on Coastal Engineering (Vancouver, Canada, ASCE), pp. 419–432. Google Scholar



Comparison between Experimental Data and the Six Studied Formulas



Breaker Depth Index as a Function of Beach Slope or Deep-Water Wave Steepness

In this Appendix, the data on the breaker depth index are displayed as a function of the beach slope (deep-water wave steepness fixed) or as a function of the deep-water wave steepness (beach slope fixed). The curves correspond to the six studied formulas plus the new one (symbols for the experimental data correspond to those used in the plots of Appendix A).

B.1. Breaker Depth Index as a Function of Deep-Water Wave Steepness


B.2. Breaker Depth Index as a Function of the Beach Slope



Fit of the Function): Influence of Beach Slope m (a, c, and e) and Wave Steepness λ (b, d, and f)


Figure 1.

The ratio hb/Lb as a function of the deep-water wave steepness λ using all the data available (the straight line corresponds to Equation [12] with Ab = 0.5, and the dashed lines to a prediction within a factor of 1.25).


Figure 2.

Comparison between experimental data and the proposed formula (Equation [15]).


Figure 3.

Distribution of the difference Δγb = γb,pred − γb,exp with the new formula (the black line corresponds to the Gaussian function).


Figure 4.

Theoretical wave height distribution and fraction of breaking waves (breaking waves correspond to the shaded area).


Figure 5.

Comparison between the two random wave models with the use of different breaker depth index formulas (CL, Camenen and Larson, Equation [15]; Nai, Nairn [1990]; Rue, Ruessink, Walstra, and Southgate [2003]) for a test case involving a barred beach with two different mean slopes and offshore wave steepnesses. a) A spilling breaking case. b) A plunging breaking case.


Figure 6.

Comparison between experimental data on breaker depth index and the Irribaren number (relationship from Battjes [1974] is also shown, see Equation [4]).


Figure 7.

Comparison between the proposed formula (Equation [15]) and the Irribaren number for fixed beach slope (a) and fixed deep-water wave steepness (b).


Figure 8.

Observation of different breaker types with respect to the Irribaren number and the breaker depth index (1: spilling breaker area, 2: plunging breaker area, 3: collapsing breaker area, 4: surging area).


Table 1.

Summary of compiled laboratory data on breaker depth index and breaker types.


Table 2.

Summary of results (predictive errors) obtained with the different formulas.


Table 3.

Comparison between the iterative method (Equation [10]) and the explicit method (Equation [13]) with the use of the original f(m, λ) or the fit f*(m, λ).


[1] The major part of the research underlying this paper is being carried out in the HUMOR project supported by the European Community and the Japanese Society for the Promotion of Science.

Benoît Camenen and Magnus Larson "Predictive Formulas for Breaker Depth Index and Breaker Type," Journal of Coastal Research 2007(234), 1028-1041, (1 July 2007).
Received: 29 July 2005; Accepted: 15 February 2006; Published: 1 July 2007
beach slope
breaker depth index
breaker type
Breaking wave
deep-water wave steepness
random waves
Back to Top