Large-Eddy Simulation of Langmuir Turbulence in Pure Wind Seas

Ramsey R. Harcourt Applied Physics Laboratory, University of Washington, Seattle, Washington

Search for other papers by Ramsey R. Harcourt in
Current site
Google Scholar
PubMed
Close
and
Eric A. D’Asaro Applied Physics Laboratory, University of Washington, Seattle, Washington

Search for other papers by Eric A. D’Asaro in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

The scaling of turbulent kinetic energy (TKE) and its vertical component (VKE) in the upper ocean boundary layer, forced by realistic wind stress and surface waves including the effects of Langmuir circulations, is investigated using large-eddy simulations (LESs). The interaction of waves and turbulence is modeled by the Craik–Leibovich vortex force. Horizontally uniform surface stress τ0 and Stokes drift profiles uS(z) are specified from the 10-m wind speed U10 and the surface wave age CP/U10, where CP is the spectral peak phase speed, using an empirical surface wave spectra and an associated wave age–dependent neutral drag coefficient CD. Wave-breaking effects are not otherwise included. Mixed layer depths HML vary from 30 to 120 m, with 0.6 ≤ CP/U10 ≤ 1.2 and 8 m s−1 < U10 < 70 m s−1, thereby addressing most possible oceanic conditions where TKE production is dominated by wind and wave forcing.

The mixed layer–averaged “bulk” VKE 〈w2〉/u*2 is equally sensitive to the nondimensional Stokes e-folding depth D*S/HML and to the turbulent Langmuir number Lat = u*/US, where u* = |τ0|/ρw in water density ρw and US = |uS|z=0. Use of a D*S scale-equivalent monochromatic wave does not accurately reproduce the results using a full-surface wave spectrum with the same e-folding depth. The bulk VKE for both monochromatic and broadband spectra is accurately predicted using a surface layer (SL) Langmuir number LaSL = u*/〈uSSL, where 〈uSSL is the average Stokes drift in a surface layer 0 > z > − 0.2HML relative to that near the bottom of the mixed layer. In the wave-dominated limit LaSL → 0, turbulent vertical velocity scales as wrmsu*La−2/3SL. The mean profile (z) of VKE is characterized by a subsurface peak, the depth of which increases with D*S/HML to a maximum near 0.22HML as its relative magnitude /〈w2〉 decreases. Modestly accurate scalings for these variations are presented. The magnitude of the crosswind velocity convergence scales differently from VKE. These results predict that for pure wind seas and HML ≅ 50 m, 〈w2〉/u*2 varies from less than 1 for young waves at U10 = 10 m s−1 to about 2 for mature seas at winds greater than U10 = 30 m s−1. Preliminary comparisons with Lagrangian float data account for invariance in 〈w2〉/u*2 measurements as resulting from an inverse relationship between U10 and CP/U10 in observed regimes.

Corresponding author address: R. R. Harcourt, Applied Physics Laboratory, University of Washington, 1013 NE 40th St., Seattle, WA 98105. Email: harcourt@apl.washington.edu

Abstract

The scaling of turbulent kinetic energy (TKE) and its vertical component (VKE) in the upper ocean boundary layer, forced by realistic wind stress and surface waves including the effects of Langmuir circulations, is investigated using large-eddy simulations (LESs). The interaction of waves and turbulence is modeled by the Craik–Leibovich vortex force. Horizontally uniform surface stress τ0 and Stokes drift profiles uS(z) are specified from the 10-m wind speed U10 and the surface wave age CP/U10, where CP is the spectral peak phase speed, using an empirical surface wave spectra and an associated wave age–dependent neutral drag coefficient CD. Wave-breaking effects are not otherwise included. Mixed layer depths HML vary from 30 to 120 m, with 0.6 ≤ CP/U10 ≤ 1.2 and 8 m s−1 < U10 < 70 m s−1, thereby addressing most possible oceanic conditions where TKE production is dominated by wind and wave forcing.

The mixed layer–averaged “bulk” VKE 〈w2〉/u*2 is equally sensitive to the nondimensional Stokes e-folding depth D*S/HML and to the turbulent Langmuir number Lat = u*/US, where u* = |τ0|/ρw in water density ρw and US = |uS|z=0. Use of a D*S scale-equivalent monochromatic wave does not accurately reproduce the results using a full-surface wave spectrum with the same e-folding depth. The bulk VKE for both monochromatic and broadband spectra is accurately predicted using a surface layer (SL) Langmuir number LaSL = u*/〈uSSL, where 〈uSSL is the average Stokes drift in a surface layer 0 > z > − 0.2HML relative to that near the bottom of the mixed layer. In the wave-dominated limit LaSL → 0, turbulent vertical velocity scales as wrmsu*La−2/3SL. The mean profile (z) of VKE is characterized by a subsurface peak, the depth of which increases with D*S/HML to a maximum near 0.22HML as its relative magnitude /〈w2〉 decreases. Modestly accurate scalings for these variations are presented. The magnitude of the crosswind velocity convergence scales differently from VKE. These results predict that for pure wind seas and HML ≅ 50 m, 〈w2〉/u*2 varies from less than 1 for young waves at U10 = 10 m s−1 to about 2 for mature seas at winds greater than U10 = 30 m s−1. Preliminary comparisons with Lagrangian float data account for invariance in 〈w2〉/u*2 measurements as resulting from an inverse relationship between U10 and CP/U10 in observed regimes.

Corresponding author address: R. R. Harcourt, Applied Physics Laboratory, University of Washington, 1013 NE 40th St., Seattle, WA 98105. Email: harcourt@apl.washington.edu

Save
  • Agrawal, Y. C., E. A. Terray, M. A. Donelan, P. A. Hwang, A. J. Williams III, W. M. Drennan, K. K. Kahma, and S. A. Kitaigorodskii, 1992: Enhanced dissipation of kinetic energy beneath breaking waves. Nature, 359 , 219220.

    • Search Google Scholar
    • Export Citation
  • Andrews, D. G., and M. E. McIntyre, 1978: An exact theory of nonlinear waves on a Lagrangian-mean flow. J. Fluid Mech., 89 , 609646.

  • Anis, A., and J. N. Moum, 1995: Surface wave–turbulence interaction: Scaling ε(z) near the sea surface. J. Phys. Oceanogr., 25 , 20252045.

    • Search Google Scholar
    • Export Citation
  • Banner, M. L., 1990: Equilibrium spectra of wind waves. J. Phys. Oceanogr., 20 , 966984.

  • Craik, A. D. D., and S. Leibovich, 1976: A rational model for Langmuir circulation. J. Fluid Mech., 73 , 401426.

  • D’Asaro, E. A., 2001: Turbulent vertical kinetic energy in the ocean mixed layer. J. Phys. Oceanogr., 31 , 35303537.

  • D’Asaro, E. A., 2003: The ocean boundary layer below Hurricane Dennis. J. Phys. Oceanogr., 33 , 561579.

  • D’Asaro, E. A., and C. McNeil, 2006: Air-sea gas exchange at extreme wind speeds. J. Mar. Syst., 66 , 92109. doi:10.1016/j.jmarsys.2006. 06.007.

    • Search Google Scholar
    • Export Citation
  • Deardorff, J. W., 1980: Stratocumulus-capped mixed layers derived from a three-dimensional model. Bound.-Layer Meteor., 18 , 494527.

  • Donelan, M. A., J. Hamilton, and W. H. Hui, 1985: Directional spectra of wind-generated waves. Philos. Trans. Roy. Soc. London, A315 , 509562.

    • Search Google Scholar
    • Export Citation
  • Donelan, M. A., F. W. Dobson, S. D. Smith, and R. J. Anderson, 1993: On the dependence of sea surface roughness on wave development. J. Phys. Oceanogr., 23 , 21432149.

    • Search Google Scholar
    • Export Citation
  • Donelan, M. A., B. K. Haus, N. Reul, W. J. Plant, M. Stiassnie, H. C. Graber, O. B. Brown, and E. S. Saltzman, 2004: On the limiting aerodynamic roughness of the ocean in very strong winds. Geophys. Res. Lett., 31 .L18306, doi:10.1029/2004GL019460.

    • Search Google Scholar
    • Export Citation
  • Drennan, W. M., M. A. Donelan, E. A. Terray, and K. B. Katdaros, 1996: Oceanic turbulence dissipation measurements in SWADE. J. Phys. Oceanogr., 26 , 808815.

    • Search Google Scholar
    • Export Citation
  • Drennan, W. M., H. C. Graber, D. Hauser, and C. Quentin, 2003: On the wave age dependence of wind stress over pure wind seas. J. Geophys. Res., 108 .8062, doi:10.1029/2000JC000715.

    • Search Google Scholar
    • Export Citation
  • Duncan, J. H., 1981: An experimental investigation of breaking waves produced by a towed hydrofoil. Proc. Roy. Soc. London, A377 , 331348.

    • Search Google Scholar
    • Export Citation
  • Garwood JR, R. W., S. M. Isakari, and P. C. Gallacher, 1994: Thermobaric convection. The Polar Oceans and Their Role in Shaping the Global Environment,Geophys. Monogr., Vol. 85, Amer. Geophys. Union, 199–209.

    • Search Google Scholar
    • Export Citation
  • Harcourt, R. R., 1999: Numerical simulation of deep convection and the response of drifters in the Laborador Sea. Ph.D. thesis, University of California, Santa Cruz, 367 pp.

  • Harcourt, R. R., 2005: Thermobaric cabbeling near Maud Rise: Theory and large eddy simulation. Prog. Oceanogr., 67 , 186244.

  • Harcourt, R. R., E. L. Steffen, R. W. Garwood, and E. A. D’Asaro, 2002: Fully Lagrangian floats in Labrador Sea deep convection: Comparison of numerical and experimental results. J. Phys. Oceanogr., 32 , 493509.

    • Search Google Scholar
    • Export Citation
  • Hasselmann, K., and Coauthors, 1973: Measurements of wind wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Herausgegeben von Deutsch. Hydrograph. Institut., Reihe A, 12, 95 pp.

  • Jessup, A. T., C. J. Zappa, and H. Yeh, 1997: Defining and quantifying microscale wave breaking with infrared imagery. J. Geophys. Res., 102 , 2314523153.

    • Search Google Scholar
    • Export Citation
  • Kenyon, K. E., 1969: Stokes drift for random gravity waves. J. Geophys. Res., 74 , 69916994.

  • Klemp, J. B., and D. R. Duran, 1983: An upper boundary condition permitting internal gravity wave radiation in numerical mesoscale models. Mon. Wea. Rev., 111 , 430444.

    • Search Google Scholar
    • Export Citation
  • Large, W. G., and S. Pond, 1981: Open ocean momentum flux measurements in moderate to strong winds. J. Phys. Oceanogr., 11 , 324336.

  • Large, W. G., J. C. McWilliams, and S. C. Doney, 1994: Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parameterization. Rev. Geophys., 32 , 363403.

    • Search Google Scholar
    • Export Citation
  • Leibovich, S., 1980: On wave-current interaction theories of Langmuir circulations. J. Fluid Mech., 99 , 715724.

  • Leibovich, S., 1983: The form and dynamics of Langmuir circulations. Annu. Rev. Fluid Mech., 15 , 391427.

  • Lewis, A. W., and R. N. Allos, 1990: JONSWAP’s parameters: Sorting out inconsistencies. Ocean Eng., 17 , 409415.

  • Lherminier, P., R. R. Harcourt, R. W. Garwood Jr., and J-C. Gascard, 2001: Interpretation of mean vertical velocities measured by isobaric floats during deep convection. J. Mar. Syst., 29 , 221237.

    • Search Google Scholar
    • Export Citation
  • Li, M., and C. Garrett, 1993: Cell merging and the jet/downwelling ratio in Langmuir circulation. J. Mar. Res., 51 , 737769.

  • Li, M., C. Garrett, and E. Skyllingstad, 2005: A regime diagram for classifying turbulent large eddies in the upper ocean. Deep-Sea Res. I, 52 , 259278.

    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., and P. P. Sullivan, 2000: Vertical mixing by Langmuir circulations. Spill Sci. Technol. Bull., 6 , 225237.

  • McWilliams, J. C., P. P. Sullivan, and C-H. Moeng, 1997: Langmuir turbulence in the ocean. J. Fluid Mech., 334 , 130.

  • Mellor, G. L., and T. Yamada, 1982: Development of a turbulence closure model for geophysical fluid problems. Rev. Geophys. Space Phys., 20 , 851875.

    • Search Google Scholar
    • Export Citation
  • Melville, W. K., 1996: The role of surface-wave breaking in air–sea interaction. Annu. Rev. Fluid Mech., 28 , 279321.

  • Min, H. S., and Y. Noh, 2004: Influence of the surface heating on Langmuir circulation. J. Phys. Oceanogr., 34 , 26302641.

  • Moeng, C-H., 1984: A large-eddy-simulation model for the study of planetary boundary-layer turbulence. J. Atmos. Sci., 41 , 20522062.

  • Moeng, C-H., and J. C. Wyngaard, 1988: Spectral analysis of large-eddy simulations of the convective boundary layer. J. Atmos. Sci., 45 , 35733587.

    • Search Google Scholar
    • Export Citation
  • Noh, Y., H. S. Min, and S. Raasch, 2004: Large eddy simulation of the ocean mixed layer: The effects of wave breaking and Langmuir circulation. J. Phys. Oceanogr., 34 , 720735.

    • Search Google Scholar
    • Export Citation
  • Phillips, O. M., 1958: The equilibrium range in the spectrum of wind-generated waves. J. Fluid Mech., 4 , 426434.

  • Plueddemann, A. J., J. A. Smith, D. M. Farmer, R. A. Weller, W. R. Crawford, R. Pinkel, S. Vagle, and A. Gnanadesikan, 1996: Structure and variability of Langmuir circulation during the Surface Waves Processes Program. J. Geophys. Res., 101 , 35253543.

    • Search Google Scholar
    • Export Citation
  • Rapp, R. J., and W. K. Melville, 1990: Laboratory measurements of deep-water breaking waves. Philos. Trans. Roy. Soc. London, 331A , 735800.

    • Search Google Scholar
    • Export Citation
  • Skyllingstad, E. D., and D. W. Denbo, 1995: An ocean large-eddy simulation of Langmuir circulations and convection in the surface mixed layer. J. Geophys. Res., 100 , 85018522.

    • Search Google Scholar
    • Export Citation
  • Skyllingstad, E., W. Smyth, and G. Crawford, 2000: Resonant wind-driven mixing in the ocean boundary layer. J. Phys. Oceanogr., 30 , 18661890.

    • Search Google Scholar
    • Export Citation
  • Smith, J. A., 1992: Observed growth of Langmuir circulation. J. Geophys. Res., 97 , 56515664.

  • Smith, J. A., 2001: Observations and theories of Langmuir circulation: A story of mixing. Fluid Mechanics and the Environment: Dynamical Approaches, J. L. Lumley, Ed., Springer, 295–314.

    • Search Google Scholar
    • Export Citation
  • Smith, J. A., 2006: Observed variability of ocean wave Stokes drift, and the Eulerian response to passing groups. J. Phys. Oceanogr., 36 , 13811402.

    • Search Google Scholar
    • Export Citation
  • Smith, S. D., 1988: Coefficients for sea surface wind stress, heat flux, and wind profiles as a function of wind speed and temperature. J. Geophys. Res., 93 , 1546715472.

    • Search Google Scholar
    • Export Citation
  • Smyth, W. D., E. D. Skyllingstad, C. B. Crawford, and H. Wijesekera, 2002: Non-local fluxes in strokes drift effects in the K-profile parametrization. Ocean Dyn., 52 , 104115.

    • Search Google Scholar
    • Export Citation
  • Sullivan, P. P., J. C. McWilliams, and W. K. Melville, 2004a: The oceanic boundary layer driven by wave breaking with stochastic variability. Part 1. Direct numerical simulations. J. Fluid Mech., 507 , 143174.

    • Search Google Scholar
    • Export Citation
  • Sullivan, P. P., J. C. McWilliams, and W. K. Melville, 2004b: Impacts of breaking waves and Langmuir circulations on the ocean mixed layer in high winds. Preprint, 26th Conf. on Hurricanes and Tropical Meteorology, Miami, FL, Amer. Meteor. Soc., 3A.2. [Available online at http://ams.confex.com/ams/pdfpaers/75679.pdf.].

  • Sullivan, P. P., J. C. McWilliams, and W. K. Melville, 2005: Surface waves and ocean mixing: Insights from numerical simulations with stochastic surface forcing. Rogue Waves: Proc. 14th ‘Aha Huliko‘a Hawaiian Winter Workshop, Honolulu, HI, University of Hawaii at Manoa, 147–155.

  • Terray, E. A., A. J. Williams, M. A. Donelan, W. M. Drennan, Y. C. Agrawal, K. K. Kahma, S. A. Kitaigorodskii, and P. A. Hwang, 1996: Estimates of kinetic energy dissipation under breaking waves. J. Phys. Oceanogr., 26 , 792807.

    • Search Google Scholar
    • Export Citation
  • Thorpe, S. A., 2004: Langmuir circulation. Annu. Rev. Fluid Mech., 36 , 5579.

  • Toba, Y., 1973: Local balance in the air-sea boundary process. J. Oceanogr. Soc. Japan, 29 , 209220.

  • Tseng, R-S., and E. A. D’Asaro, 2004: Measurements of turbulent vertical kinetic energy in the ocean mixed layer from Lagrangian floats. J. Phys. Oceanogr., 34 , 19841990.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1873 1160 309
PDF Downloads 610 152 10