Reply

George Mellor Princeton University, Princeton, New Jersey

Search for other papers by George Mellor in
Current site
Google Scholar
PubMed
Close
Full access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Corresponding author address: George Mellor, Program in Atmospheric and Oceanic Sciences, Princeton University, P.O. Box CN710, Sayre Hall Princeton, NJ 08544-0710. E-mail:glmellor@princeton.edu

The original article that was the subject of this comment/reply can be found at http://journals.ametsoc.org/doi/abs/10.1175/2008JPO3971.1.

Corresponding author address: George Mellor, Program in Atmospheric and Oceanic Sciences, Princeton University, P.O. Box CN710, Sayre Hall Princeton, NJ 08544-0710. E-mail:glmellor@princeton.edu

The original article that was the subject of this comment/reply can be found at http://journals.ametsoc.org/doi/abs/10.1175/2008JPO3971.1.

1. Introduction

The comments by Bennis and Ardhuin (2011, hereinafter BA11) are partially correct but, for the most part, are incorrect.

2. Discussion

The paper by Mellor (2008) contained some worthy aspects: in particular, the expression for the depth-dependent radiation stress,
e1
where the F functions are defined in Mellor (2003, 2008) and BA11. Equation (1) includes a Kronecker delta coefficient δαβ and a Dirac delta function , where is the mean elevation (, where h is the water depth). The problem with Mellor (2008) is that the radiation stress appears as ∂Sαβ/∂xβ in the Cartesian form of the phase-averaged momentum equation; additional terms were missing (see appendix A of Mellor 2005) that would cancel the Leibnitz terms in BA11’s Eq. (17). Thus, a vertical integration of the momentum equation did not equate to the corresponding equations in Phillips (1977) and Smith (2006). BA11 are correct about this issue. A recent corrigendum corrects the Mellor (2003) paper so that Eq. (1) does appear as ∂(DSαβ)/∂xβ but in the sigma form of the phase-averaged momentum equation. As seen below, it can be vertically integrated so that agreement with Longuet-Higgins and Stewart (1964), Phillips (1977), and Smith (2006) is obtained.
The delta function in Eq. (1) is defined in Mellor (2008) and in BA11’s Eq. (8) such that if but
eq1
for small ε and . The “approximated” form in Eq. (14) of BA11 is incorrect, however. A proper form might be
e2a
where the right side is a “generalized function” (Lighthill 1958; Greenberg 1978) or impulse function. Differentiation with respect to x, where , and then integration with respect to z′ yields zero, not infinity as BA11 would have it. A seemingly small change—replacing z′ by −|z′|—in Eq. (14) such that
e2b
can also be an acceptable impulse function. As in Eq. (2a), the integral of Eq. (2b) with respect to z′ after differentiation with respect to x is zero and is so by inspection.
In their section 3, I believe that BA11 solve the two-dimensional, sigma-coordinate equation
e3a
where τ(0) = τ(−1) = 0 and from Eq. (1)
e3c
Steady state is obtained. The “sigma” coordinate . Here, as derived in Mellor (2003), U is the current velocity plus Stokes drift. I use the symbol Ω instead of W to distinguish it from the Cartesian vertical velocity. Also, I write Eq. (3a) in flux form, which, with the continuity equation (Mellor 2003), is convertible to Eq. (18) of BA11. The form of Eq. (3a) is convenient and readily integrates from ς = −1 to ς = 0, yielding
e4
where
e5
with c and cg being the phase and group speeds, respectively. The two terms on the right of Eq. (5) correspond to those in Eq. (3c). The second term on the left of Eq. (4) was converted to ∂(M2/D)/∂x by Phillips with the help of some fairly restrictive assumptions as to the profile shape of U(ς). BA11 apparently accept Eq. (3a), derived in Mellor (2003), as a valid equation; it is not clear how BA11 regard the last term in Eq. (5) or its origin, the last term in Eq. (3c).
I can understand BA11’s concern over the singular term in Eq. (3c). It follows directly from Longuet-Higgins and Stewart (1964), however. It is a contribution to wave momentum due to the (hydrostatic) pressure imbedded in the crest and trough at the wave surface and given by the expression
e6
where, here, overbars represent phase averaging and is the instantaneous wave elevation relative to . In a depth-dependent context, the term is clearly concentrated at the surface; thus, the delta representation. When integrated, is subsumed in Eq. (5) and loses its singular identity.

Thus, if the surface of a flow is wavy then there must be a surface contribution to the phase-averaged momentum equation due to the wave’s intrinsic surface pressure field. This effect is missing in the solution labeled “exact” in BA11’s Fig. 3, which, therefore, is incorrect in my opinion. A more thorough discussion of the representation of wave-induced pressure is in Mellor (2011).

The existence of a concentrated surface forcing term in the momentum equation does dictate a realistic need to include a subsurface viscous or eddy viscosity stress term. Note that such forcing acts similar to a surface wind stress, which also requires a subsurface viscous or eddy viscosity term. Thus, instead of the exact solution, the long-dashed curve in BA11’s Fig. 3 may be correct, although I have no knowledge as to upstream and downstream boundary conditions and other details.

A question does arise as to whether the U in Eq. (3b) should include or exclude the Stokes drift in contrast to U terms in Eq. (3a), which do include Stokes drift. I hope to answer that question in a paper now in progress.

BA11 seem to suggest that models such as those that are based on Eq. (3) and some form of Sαβ that is based on Airy wave solution are not valid. Their position seems to be that a full Laplacian-type solution involving 10 vertical modes is needed. I disagree with that doomsday speculation.

Another question concerns the validity of Eqs. (3) and (1) in shallow water with a bottom slope; the question is, how shallow or how steep is the slope? Let k, a, and h be wavenumber, amplitude, and water depth, respectively. Whereas, in the development of Eqs. (3) and (1), terms of O(ka)2 were retained and terms of O(ka)4 were discarded, it can be shown (details available on request) that terms of O[ka(∂h/∂x)/sinh(kh)]2 should be less than or equal to O(ka)4; that is, (∂h/∂x)/sinh(kh) should be small: comparable to or less than the wave slope (ka).

3. Summary

The BA11 criticism of Mellor (2008) is correct and is partially the subject of the recently published corrigendum to Mellor (2003).

Their generalized definition of a Dirac delta function is incorrect, does not conform to standard usage (Lighthill 1958; Greenberg 1978) or my interpretation, and leads to an erroneous conclusion.

BA11 object to the existence of a surface-trapped, wave-related pressure contribution to the momentum balance. It is, however, implicit in the derivations of Longuet-Higgins and Stewart (1964) and Phillips (1977) and is explicitly part of their vertically integrated wave radiation stress.

The basis of BA11’s skepticism as to the utility of schemes that are based on a wave radiation stress term to couple wave dynamics to general circulation dynamics is unclear to me.

REFERENCES

  • Bennis, A.-C., and F. Ardhuin, 2011: Comments on “The depth-dependent current and wave interaction equations: A revision.” J. Phys. Oceanogr., 41, 20082012.

    • Search Google Scholar
    • Export Citation
  • Greenberg, M. D., 1978. Foundations of Applied Mathematics. Prentice Hall, 636 pp.

  • Lighthill, M. J., 1958: Introduction to Fourier Analysis and Generalized Functions. Cambridge University Press, 79 pp.

  • Longuet-Higgins, M. S., and R. W. Stewart, 1964: Radiation stresses in water waves; a physical discussion with applications. Deep-Sea Res., 11, 529562.

    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2003: The three-dimensional current and surface wave equations. J. Phys. Oceanogr., 33, 19781989; Corrigendum, 41, 14171418.

    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2005: Some consequences of the three-dimensional current and surface wave equations. J. Phys. Oceanogr., 35, 22912298.

    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2008: The depth-dependent current and wave interaction equations: A revision. J. Phys. Oceanogr., 38, 25872596.

  • Mellor, G. L., 2011: Wave radiation stress. Ocean Dyn., 61, 563568.

  • Phillips, O. M., 1977: The Dynamics of the Upper Ocean. Cambridge University Press, 336 pp.

  • Smith, J. A., 2006: Wave–current interactions in finite depth. J. Phys. Oceanogr., 36, 14031419.

Save
  • Bennis, A.-C., and F. Ardhuin, 2011: Comments on “The depth-dependent current and wave interaction equations: A revision.” J. Phys. Oceanogr., 41, 20082012.

    • Search Google Scholar
    • Export Citation
  • Greenberg, M. D., 1978. Foundations of Applied Mathematics. Prentice Hall, 636 pp.

  • Lighthill, M. J., 1958: Introduction to Fourier Analysis and Generalized Functions. Cambridge University Press, 79 pp.

  • Longuet-Higgins, M. S., and R. W. Stewart, 1964: Radiation stresses in water waves; a physical discussion with applications. Deep-Sea Res., 11, 529562.

    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2003: The three-dimensional current and surface wave equations. J. Phys. Oceanogr., 33, 19781989; Corrigendum, 41, 14171418.

    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2005: Some consequences of the three-dimensional current and surface wave equations. J. Phys. Oceanogr., 35, 22912298.

    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2008: The depth-dependent current and wave interaction equations: A revision. J. Phys. Oceanogr., 38, 25872596.

  • Mellor, G. L., 2011: Wave radiation stress. Ocean Dyn., 61, 563568.

  • Phillips, O. M., 1977: The Dynamics of the Upper Ocean. Cambridge University Press, 336 pp.

  • Smith, J. A., 2006: Wave–current interactions in finite depth. J. Phys. Oceanogr., 36, 14031419.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 303 221 122
PDF Downloads 50 20 2