Reply to “Comments on ‘A Wave-Resolving Simulation of Langmuir Circulations with a Nonhydrostatic Free-Surface Model: Comparison with Craik–Leibovich Theory and an Alternative Eulerian View of the Driving Mechanism’”

Yasushi Fujiwara Graduate School of Science, Kyoto University, Kyoto, Japan

Search for other papers by Yasushi Fujiwara in
Current site
Google Scholar
PubMed
Close
,
Yutaka Yoshikawa Graduate School of Science, Kyoto University, Kyoto, Japan

Search for other papers by Yutaka Yoshikawa in
Current site
Google Scholar
PubMed
Close
, and
Yoshimasa Matsumura Atmosphere and Ocean Research Institute, The University of Tokyo, Tokyo, Japan

Search for other papers by Yoshimasa Matsumura in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

Fujiwara et al. explicitly simulated Langmuir circulations using a wave-resolving simulation (WRS) technique and found that the residual wave effect on vorticity was well represented by the vortex force of the Craik–Leibovich (CL) equation, at least in the simulated situation. In response to the simulation results, Mellor has proposed a view that ubiquitous applicability of the CL formulation is still questionable and that the three-dimensional radiation stress (3DRS) formulation that he has derived encompasses both of the vortex force effect and an effect that is lower order in terms of wave steepness. Here, these opinions are discussed in terms of the approximations used in the wave-averaged formulations. The asymptotic expansion of the Eulerian-averaged momentum equation allows the separate discussion of two different wave effects: pressure correction and torque. It is argued that the approximation adopted in Mellor’s 3DRS formulation is presumably not accurate enough to properly parameterize the wave torque effect, and possible approaches to examine its performance are proposed. We agree with the view that the applicability of the CL formulation needs further investigation. WRS will be a helpful tool for this purpose.

© 2019 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Research Fellow, Japan Society for the Promotion of Science, Tokyo, Japan.

Corresponding author: Yasushi Fujiwara, fujiwara@kugi.kyoto-u.ac.jp

The original article that was the subject of this comment/reply can be found at http://journals.ametsoc.org/doi/abs/10.1175/JPO-D-17-0199.1.

Abstract

Fujiwara et al. explicitly simulated Langmuir circulations using a wave-resolving simulation (WRS) technique and found that the residual wave effect on vorticity was well represented by the vortex force of the Craik–Leibovich (CL) equation, at least in the simulated situation. In response to the simulation results, Mellor has proposed a view that ubiquitous applicability of the CL formulation is still questionable and that the three-dimensional radiation stress (3DRS) formulation that he has derived encompasses both of the vortex force effect and an effect that is lower order in terms of wave steepness. Here, these opinions are discussed in terms of the approximations used in the wave-averaged formulations. The asymptotic expansion of the Eulerian-averaged momentum equation allows the separate discussion of two different wave effects: pressure correction and torque. It is argued that the approximation adopted in Mellor’s 3DRS formulation is presumably not accurate enough to properly parameterize the wave torque effect, and possible approaches to examine its performance are proposed. We agree with the view that the applicability of the CL formulation needs further investigation. WRS will be a helpful tool for this purpose.

© 2019 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Research Fellow, Japan Society for the Promotion of Science, Tokyo, Japan.

Corresponding author: Yasushi Fujiwara, fujiwara@kugi.kyoto-u.ac.jp

The original article that was the subject of this comment/reply can be found at http://journals.ametsoc.org/doi/abs/10.1175/JPO-D-17-0199.1.

1. Introduction

Using a latest nonhydrostatic free-surface numerical model, Fujiwara et al. (2018, hereinafter FYM) performed a wave-resolving simulation (WRS) to study the driving mechanism of simulated Langmuir circulations. In their analysis, FYM found that the simulated wave-induced torque is very well represented by the curl of the vortex force in the Craik–Leibovich (CL) equation (Craik and Leibovich 1976, hereinafter CL76) in that particular case.

Mellor (2019, hereinafter M19) recently submitted a comment on FYM. M19 reviews the numerical work of FYM and puts forward his view on the various wave-averaged formulations such as the three-dimensional radiation stress (3DRS) equation (Mellor 2003, 2015, hereinafter M03 and M15, respectively) and the CL equation. M19 shares the common question with FYM: when does the CL equation properly parameterize the residual wave effect, and when does it not? His answer is as follows: (i) its ubiquitous applicability is still questionable, and (ii) the 3DRS of M15 contains both the vortex force and an effect that is lower order in terms of wave steepness and is missing in CL76.

In this reply, we discuss the opinions of M19 through the consideration of the wave effect terms in the wave-averaged momentum equation. Since wave-averaged formulations essentially approximate the wave-associated fluctuation to obtain a useful parameterization of its residual effects, it is crucial to know the required accuracy of the approximation. Here we feature the Eulerian wave-averaged momentum equation to identify two different wave effects (pressure correction and torque) and the accuracy of the wave-associated fluctuation required to approximate these effects. This separation of effects clarifies the “lower-order effect” mentioned by M19.

2. Required accuracy of wave motion in the Eulerian framework

Many wave-averaged formulations are derived using the asymptotic expansion (e.g., Longuet-Higgins and Stewart 1962; Craik and Leibovich 1976; McWilliams et al. 2004; Aiki and Greatbatch 2014). While it might not be so physically intuitive as nonasymptotic derivations (e.g., Longuet-Higgins and Stewart 1964; Smith 2006), the asymptotic expansion allows us a solid discussion of accuracy and errors with the hierarchical series of equations. Here, we shall discuss the required accuracy of approximations based on the asymptotic approach, partly following the derivation by CL76.

In the Eulerian wave-averaged momentum equation, the residual wave effect is
e1
Here is velocity and is vorticity. The overline and prime denote the Eulerian temporal average over one (or several) wave period(s) and the deviation from the Eulerian temporal average, respectively. The first term is a gradient of a scalar and does not affect vorticity. It does modify the pressure field, whose horizontal variation changes mean surface elevation, the effect known as the wave setup/setdown. The second term is responsible for any forces accompanied by torque. The expression in Lamb’s form [Eq. (1)] is useful to distinguish different wave effects, namely, pressure correction and torque.

The use of different coordinate systems leads to different forms of wave-averaged equations, but they should essentially contain the same physics behind. Therefore, the wave-averaged equation based on the Eulerian average [Eq. (1)] and the one based on the ς coordinate (M03, M15) are equivalently accurate as long as we have the completely accurate field of the wave-associated fluctuation (quantities such as ). In practice, we usually do not have such information because there is no general analytic solution of waves in the presence of background current and topography. One possible solution is to explicitly simulate wave motions (i.e., WRS) and thereby numerically obtain all the information of both wave-averaged motion and fluctuation, but it is computationally demanding and thus often unreasonable for realistic simulations. Another approach is to use an asymptotically approximate form of , etc. to obtain useful parameterizations of the wave effect terms.

In the following paragraphs, we shall discuss the accuracy of the approximated wave orbital motions required to properly parameterize the wave residual effect terms in the Eulerian-averaged form [Eq. (1)]. We feature the Eulerian framework because it provides a clear view of the terms responsible for the pressure correction and torque effects, which act with different magnitudes. To facilitate the discussion, we follow the scaling assumptions adopted by CL76, namely, (i) the wave-averaged velocity is , and (ii) the spatial scale of is , where ε is wave steepness, c is wave phase speed, and k is wavenumber. Under these assumptions, the wave-averaged flow and the fluctuation (wave motion in the leading order) from the averaged flow are expanded in series:
e2
e3
The first- and the second-order fluctuation and follow the equations
e4
e5
where p denotes the kinematic pressure minus the hydrostatic pressure from a fixed level ( using the raw pressure P and the reference density ). Compressibility, viscosity, stratification, and system rotation are neglected, and is used assuming the temporal scale separation. Because only gradient forces are present in the rhs of Eqs. (4) and (5), the fluctuation is irrotational up to the second order (). This justifies the use of the Airy wave (the irrotational linear wave) solution to approximate . The third-order fluctuation follows
e6
Therefore, reflects the interaction of the background current and the leading-order wave motion , and it is accompanied by vorticity () due to the nongradient term .
Now let us consider which order of wave-associated fluctuation contributes to the wave torque, which is caused by the term . Substituting the series expansion (2) into the term using the fact , we obtain
e7
Therefore, the parameterization of the leading-order torque effect requires the information of the third-order fluctuation (fluctuating vorticity; see also Fig. 7 of FYM). In other words, the Airy wave solution (the first-order wave) alone is unable to represent the torque effect in the Eulerian derivation. CL76 dealt with the wave-averaged vorticity equation and showed that can be transformed to , where is the Stokes drift and is . After “uncurling” the wave-averaged vorticity equation, we obtain the wave-averaged momentum equation with the vortex force and an unknown gradient force (Leibovich 1977).1
Next, let us consider the pressure correction effect. Substituting the series expansion (2) into the term yields
e8
where the leading-order term is unaffected by the wave-averaged current . Therefore, the information of (i.e., the Airy wave) would be sufficiently accurate to capture the setup/setdown.2

To summarize, the asymptotic expansion of the Eulerian-averaged momentum equation clarifies the two different wave effects, namely, the pressure correction and torque, and the accuracy required to properly parameterize these effects. The Airy wave solution would be sufficiently accurate to parameterize the pressure correction, while the higher-order solution with the wave-current interaction effect must be used to parameterize the leading-order torque. The pressure correction is “the lower-order effect” mentioned by M19, and it cannot be parameterized through the “uncurling” of wave-averaged vorticity equation, as was done in CL76 and Leibovich (1977). Meanwhile, as noted before, the pressure correction does not affect vorticity.

3. Discussion and conclusions

M03 and M15 derived a wave-averaged momentum equation using the so-called ς coordinate, which follows the fluctuating upper surface. The wave motion in the residual wave effect terms was approximated with the Airy wave solution in the derivations [similar to truncating the second- and higher-order terms in Eq. (2)], obtaining the second-order accurate 3DRS ( in M19). This approximation is less accurate than the one required to derive the torque effect in the Eulerian formulation, that is, the interaction of the wave-averaged current and the lower-order wave motion that appears in the higher-order fluctuation is neglected.

Although the Eulerian discussion cannot be directly applied to the ς coordinate, we consider that the torque effect parameterization of M03/M15 needs to be carefully examined. As suggested by the opinion ii of M19, the horizontal component of the M03/M15 formulation does contain a vortex force term due to horizontal Stokes drift, which arises from expanding the advection term of the total momentum ( in M19). However, the 3DRS term is second-order accurate and may contain an error that contaminates the leading-order torque effect terms, which is typically higher-order than (e.g., the fourth-order in CL76). Also, while unmentioned in M19, the vertical component of the wave-averaged momentum equation needs to be discussed.

As the reliability of the torque effect parameterization in the Airy wave–based formulation is poorly supported, we suggest that solid experimental or theoretical evidence should be given. One simple way to give such supports would be to compare the vorticity evolution in a simulation based on the Airy wave–based formulation and a corresponding WRS. A reference result of a simple test (e.g., the development of laminar Langmuir circulations under monochromatic waves in a periodic domain) with a well-controlled WRS will be useful, and the authors are currently preparing such a result. Another more theoretical way might be to derive a wave-averaged vorticity equation using the asymptotic expansion and examine the leading-order wave torque term from the general 3DRS (i.e., including all higher-order motions). If it is properly represented with Airy wave approximation or less significant than the torque terms corresponding to the vortex force term in Eq. (12) of M19, the Airy wave approximation would be sufficiently accurate. The asymptotic discussion in the present manuscript could be used as a guideline for such discussion. These attempts will clarify whether the accuracy of the approximation is sufficient or not, and either conclusion will be largely beneficial for the users of the formulations.

Based on various assumptions and frameworks, various wave-averaged formulations have been proposed. Some formulations include both the wave setup/setdown effect and the torque effect (e.g., McWilliams et al. 2004; Smith 2006; Ardhuin et al. 2008). However, it is always worthwhile to pay attention to the validity of the formulations because the assumptions that they are based on may not be strictly satisfied all over the simulated domain and period,3 where overlooked terms in the derivation may become significant. Also, the models of wave momentum evolution become questionable when there is small-scale current due to turbulence or fronts because they are at best based on the wave action conservation and the ray equation (e.g., McWilliams et al. 2004), which are satisfied only when the spatial scale of the current is much larger than the wavelength. It is important to monitor the performance of formulations in such cases. In this sense, the opinion i of M19 (ubiquitous application of the CL equation is questionable) is important.

Again, well-controlled WRS will be a useful tool to examine the performance of wave-averaged formulations. Using WRS, FYM gave a support for the vortex force parameterization of the torque effect in the simulated particular conditions, which roughly satisfied the asymptotic assumptions of CL76.4 As mentioned in their conclusion, there are many other possible situations yet to be studied, and further investigation is required to expand our understanding of wave–current interaction in the upper ocean. When WRS is used in such attempts, careful examinations of numerical model performance and analysis method are necessary since the residual wave effect on current is quite subtle compared to the forces driving orbital motions.

Acknowledgments

This research is supported by Grant-in-Aid for Scientific Research on Innovative Areas (MEXT KAKENHI Grants JP15H05824 and JP15H05825) and Grant-in-Aid for JSPS Research Fellow (JSPS KAKENHI Grant JP17J07923).

REFERENCES

  • Aiki, H., and R. J. Greatbatch, 2014: A new expression for the form stress term in the vertically Lagrangian mean framework for the effect of surface waves on the upper ocean circulation. J. Phys. Oceanogr., 44, 323, https://doi.org/10.1175/JPO-D-12-0228.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ardhuin, F., N. Rascle, and K. A. Belibassakis, 2008: Explicit wave-averaged primitive equations using a generalized Lagrangian mean. Ocean Modell., 20, 3560, https://doi.org/10.1016/j.ocemod.2007.07.001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Craik, A. D. D., and S. Leibovich, 1976: A rational model for Langmuir circulations. J. Fluid Mech., 73, 401426, https://doi.org/10.1017/S0022112076001420.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fujiwara, Y., Y. Yoshikawa, and Y. Matsumura, 2018: A wave-resolving simulation of Langmuir circulations with a nonhydrostatic free-surface model: Comparison with Craik–Leibovich theory and an alternative Eulerian view of the driving mechanism. J. Phys. Oceanogr., 48, 16911708, https://doi.org/10.1175/JPO-D-17-0199.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Leibovich, S., 1977: Convective instability of stably stratified water in the ocean. J. Fluid Mech., 82, 561581, https://doi.org/10.1017/S0022112077000846.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Longuet-Higgins, M. S., and R. Stewart, 1962: Radiation stress and mass transport in gravity waves, with application to ‘surf beats.’ J. Fluid Mech., 13, 481504, https://doi.org/10.1017/S0022112062000877.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Longuet-Higgins, M. S., and R. Stewart, 1964: Radiation stresses in water waves; A physical discussion, with applications. Deep-Sea Res. Oceanogr. Abstr., 11, 529562, https://doi.org/10.1016/0011-7471(64)90001-4.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., J. M. Restrepo, and E. M. Lane, 2004: An asymptotic theory for the interaction of waves and currents in coastal waters. J. Fluid Mech., 511, 135178, https://doi.org/10.1017/S0022112004009358.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2003: The three-dimensional current and surface wave equations. J. Phys. Oceanogr., 33, 19781989, https://doi.org/10.1175/1520-0485(2003)033<1978:TTCASW>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2015: A combined derivation of the integrated and vertically resolved, coupled wave–current equations. J. Phys. Oceanogr., 45, 14531463, https://doi.org/10.1175/JPO-D-14-0112.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2016: On theories dealing with the interaction of surface waves and ocean circulation. J. Geophys. Res. Oceans, 121, 44744486, https://doi.org/10.1002/2016JC011768.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2019: Comments on “A wave-resolving simulation of Langmuir circulations with a nonhydrostatic free-surface model: Comparison with Craik–Leibovich Theory and an alternative Eulerian view of the driving mechanism.” J. Phys. Oceanogr., 49, 885888, https://doi.org/10.1175/JPO-D-18-0222.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Smith, J. A., 2006: Wave–current interactions in finite depth. J. Phys. Oceanogr., 36, 14031419, https://doi.org/10.1175/JPO2911.1.

1

Mellor (2016) and M19 argue that there is some missing term in the CL equation that was dropped in “curl–uncurl” processes. Recall that such missing information must take a gradient form.

2

This argument is somewhat crude because the term , whose divergence is generally nonzero, can also cause wave setup/setdown. When the horizontal variation of is small, the wave setup/setdown effect from can be smaller than from . In such cases, the higher-order fluctuation must be taken into account to correctly estimate the setup/setdown.

3

For example, in turbulent flows the energy spans over a wide range of spatiotemporal scales, and there is likely a scale at which the asymptotic assumption (e.g., varies over distance ) is not satisfied.

4

With the parameters listed at the end of the section 3 of FYM, , and . The conditions of CL76 (see section 2) are roughly satisfied with an error of one order of magnitude.

Save
  • Aiki, H., and R. J. Greatbatch, 2014: A new expression for the form stress term in the vertically Lagrangian mean framework for the effect of surface waves on the upper ocean circulation. J. Phys. Oceanogr., 44, 323, https://doi.org/10.1175/JPO-D-12-0228.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ardhuin, F., N. Rascle, and K. A. Belibassakis, 2008: Explicit wave-averaged primitive equations using a generalized Lagrangian mean. Ocean Modell., 20, 3560, https://doi.org/10.1016/j.ocemod.2007.07.001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Craik, A. D. D., and S. Leibovich, 1976: A rational model for Langmuir circulations. J. Fluid Mech., 73, 401426, https://doi.org/10.1017/S0022112076001420.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fujiwara, Y., Y. Yoshikawa, and Y. Matsumura, 2018: A wave-resolving simulation of Langmuir circulations with a nonhydrostatic free-surface model: Comparison with Craik–Leibovich theory and an alternative Eulerian view of the driving mechanism. J. Phys. Oceanogr., 48, 16911708, https://doi.org/10.1175/JPO-D-17-0199.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Leibovich, S., 1977: Convective instability of stably stratified water in the ocean. J. Fluid Mech., 82, 561581, https://doi.org/10.1017/S0022112077000846.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Longuet-Higgins, M. S., and R. Stewart, 1962: Radiation stress and mass transport in gravity waves, with application to ‘surf beats.’ J. Fluid Mech., 13, 481504, https://doi.org/10.1017/S0022112062000877.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Longuet-Higgins, M. S., and R. Stewart, 1964: Radiation stresses in water waves; A physical discussion, with applications. Deep-Sea Res. Oceanogr. Abstr., 11, 529562, https://doi.org/10.1016/0011-7471(64)90001-4.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., J. M. Restrepo, and E. M. Lane, 2004: An asymptotic theory for the interaction of waves and currents in coastal waters. J. Fluid Mech., 511, 135178, https://doi.org/10.1017/S0022112004009358.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2003: The three-dimensional current and surface wave equations. J. Phys. Oceanogr., 33, 19781989, https://doi.org/10.1175/1520-0485(2003)033<1978:TTCASW>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2015: A combined derivation of the integrated and vertically resolved, coupled wave–current equations. J. Phys. Oceanogr., 45, 14531463, https://doi.org/10.1175/JPO-D-14-0112.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2016: On theories dealing with the interaction of surface waves and ocean circulation. J. Geophys. Res. Oceans, 121, 44744486, https://doi.org/10.1002/2016JC011768.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., 2019: Comments on “A wave-resolving simulation of Langmuir circulations with a nonhydrostatic free-surface model: Comparison with Craik–Leibovich Theory and an alternative Eulerian view of the driving mechanism.” J. Phys. Oceanogr., 49, 885888, https://doi.org/10.1175/JPO-D-18-0222.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Smith, J. A., 2006: Wave–current interactions in finite depth. J. Phys. Oceanogr., 36, 14031419, https://doi.org/10.1175/JPO2911.1.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 818 536 144
PDF Downloads 254 77 43