Transient Evolution of Langmuir Turbulence in Ocean Boundary Layers Driven by Hurricane Winds and Waves

Peter P. Sullivan National Center for Atmospheric Research, Boulder, Colorado

Search for other papers by Peter P. Sullivan in
Current site
Google Scholar
PubMed
Close
,
Leonel Romero Scripps Institution of Oceanography, La Jolla, California

Search for other papers by Leonel Romero in
Current site
Google Scholar
PubMed
Close
,
James C. McWilliams Department of Atmospheric and Oceanic Science, University of California, Los Angeles, Los Angeles, California

Search for other papers by James C. McWilliams in
Current site
Google Scholar
PubMed
Close
, and
W. Kendall Melville Scripps Institution of Oceanography, La Jolla, California

Search for other papers by W. Kendall Melville 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.

Abstract

A large-eddy simulation (LES) model, which adopts wave-averaged equations with vortex force, is used to investigate Langmuir turbulence and ocean boundary layer (OBL) dynamics in high-wind hurricane conditions. The temporally evolving spatially asymmetric wind and wave Stokes drift velocity imposed in the LES are generated by a spectral wave prediction model adapted to Hurricane Frances traveling at a speed of 5.5 m s−1. The potency of Langmuir turbulence depends on the turbulent Langmuir number, the wind–Stokes drift alignment, and the depth scale of the Stokes profile Ds relative to the OBL depth h. At the time of maximum winds, large-scale vigorous coherent cells develop on the right-hand side of the storm under the inertially rotating winds; the Stokes drift velocity is well tuned to the surface winds. Much weaker cells develop on the left-hand side of the storm, partly because of reduced Stokes production. With misaligned winds and waves the vertical momentum fluxes can be counter to the gradient of Stokes drift, and the cell orientation tracks the direction of the mean Lagrangian shear. The entrainment flux is increased by 20% and the sea surface temperature is 0.25 K cooler on the right-hand side of the storm in the presence of Langmuir turbulence. Wave effects impact entrainment when the ratio Ds/|h| > 0.75. Because of wind–wave asymmetry Langmuir cells add quantitatively to the left–right asymmetry already understood for hurricanes due to resonance. And the transient evolution of the OBL cannot be understood simply in terms of equilibrium snapshots.

Current affiliation: Marine Science Institute, University of California, Santa Barbara, Santa Barbara, California.

Corresponding author address: Peter P. Sullivan, MMM Division, NCAR, Boulder, CO 80307-3000. E-mail: pps@ucar.edu

Abstract

A large-eddy simulation (LES) model, which adopts wave-averaged equations with vortex force, is used to investigate Langmuir turbulence and ocean boundary layer (OBL) dynamics in high-wind hurricane conditions. The temporally evolving spatially asymmetric wind and wave Stokes drift velocity imposed in the LES are generated by a spectral wave prediction model adapted to Hurricane Frances traveling at a speed of 5.5 m s−1. The potency of Langmuir turbulence depends on the turbulent Langmuir number, the wind–Stokes drift alignment, and the depth scale of the Stokes profile Ds relative to the OBL depth h. At the time of maximum winds, large-scale vigorous coherent cells develop on the right-hand side of the storm under the inertially rotating winds; the Stokes drift velocity is well tuned to the surface winds. Much weaker cells develop on the left-hand side of the storm, partly because of reduced Stokes production. With misaligned winds and waves the vertical momentum fluxes can be counter to the gradient of Stokes drift, and the cell orientation tracks the direction of the mean Lagrangian shear. The entrainment flux is increased by 20% and the sea surface temperature is 0.25 K cooler on the right-hand side of the storm in the presence of Langmuir turbulence. Wave effects impact entrainment when the ratio Ds/|h| > 0.75. Because of wind–wave asymmetry Langmuir cells add quantitatively to the left–right asymmetry already understood for hurricanes due to resonance. And the transient evolution of the OBL cannot be understood simply in terms of equilibrium snapshots.

Current affiliation: Marine Science Institute, University of California, Santa Barbara, Santa Barbara, California.

Corresponding author address: Peter P. Sullivan, MMM Division, NCAR, Boulder, CO 80307-3000. E-mail: pps@ucar.edu

1. Introduction

The interaction between winds, waves, and currents ultimately determines the fate of turbulent mixing in the upper ocean boundary layer (OBL). There is a growing body of large-eddy simulation (LES) modeling studies supported by observational evidence that finds surface waves, and in particular phase-averaged wave–current interactions and stochastic wave breaking, play important roles in the dynamics of the upper ocean boundary layer (Sullivan and McWilliams 2010). The most compelling theoretical model of wave–current interactions is embodied in wave-averaged equations and their mathematical representation by a vortex force (Craik and Leibovich 1976; Leibovich 1983; McWilliams et al. 1997; McWilliams and Restrepo 1999; McWilliams et al. 2004). The latter is found to generate coherent Langmuir cells on a background of OBL turbulence that enhance turbulence kinetic energy and promote entrainment at the thermocline (Skyllingstad and Denbo 1995; McWilliams et al. 1997). Wave breaking can modulate the formation of these coherent structures, and breaking energetics increase dissipation rates in the upper OBL by several orders of magnitude compared to classic law-of-the-wall scaling (Terray et al. 1996; Sullivan et al. 2007; McWilliams et al. 2012). The regime where surface waves play an important role in the dynamics of the OBL is now referred to as Langmuir turbulence (McWilliams et al. 1997). Langmuir turbulence is created during the passage of storms, cold fronts, and tropical cyclones. These intermittent high-impact weather events, which help define the overall ocean structure and global sea surface temperature (SST), are subgrid-scale to the current generation of coupled-climate models, and thus improved parameterizations of OBLs with wave effects are generally needed (Sullivan and McWilliams 2010; Belcher et al. 2012).

Observations of wavy upper-ocean dynamics are conducted in the presence of broadband wind wave complexity, for example, wind wave disequilibrium, persistent swell, high winds with inertially rotating wind stress, and sharp fronts and submesoscale motions with breaking waves are common (Smith 1998; D’Asaro 2001; Melville et al. 2005; Capet et al. 2008; Sanford et al. 2011). LES studies of the OBL that include vortex force are, however, almost exclusively conducted for idealized steady flow regimes. Essentially, lengthy LES time integrations are carried out assuming a false independence of winds and waves; that is, there is no wind wave evolution over the long time periods used to gather turbulence statistics. This is an acceptable approximation for near wind–wave equilibrium conditions but is less certain in nonstationary wave states. Thus, a next step in complexity is exploring the dynamics of Langmuir turbulence in an evolving wind wave regime. Recent LES investigations consider wavy OBLs with a modest degree of nonequilibrium wind wave forcing, for example, idealized wind–wave misalignment (Van Roekle et al. 2012), nonstationary mixing in a shallow ocean (Kukulka et al. 2010), evolution of Langmuir turbulence over a range of wind speeds and wave age (Harcourt and D’Asaro 2008), and wind-driven mixing below the thermocline (Grant and Belcher 2011).

The present study expands our previous LES modeling investigations of Langmuir turbulence to the high wind (hurricane) driven ocean boundary layer. This regime features rapidly turning wind stress and temporally evolving wave fields that are a complex interacting mix of wind waves and swell with multiple spectral peaks propagating in directions misaligned with the surface winds (Wright et al. 2001; Fan et al. 2009). The hurricane-forced OBL then allows an assessment of Langmuir turbulence in a highly nonequilibrium wind wave regime. The classic study of Price (1981) and more recently Sanford et al. (2011) emphasize the importance of inertial resonance for ocean mixing, that is, wind stress rotating in the same direction as the underlying currents. These studies, however, do not consider OBL dynamics in the presence of the companion significant surface wave fields that high winds generate.

The goal of the present work is to examine the response of the wavy OBL to realistic time varying wind and wave forcing induced by an idealized moving hurricane. The previous LES investigation most closely related to the present work is that of Skyllingstad et al. (2000) who examined OBLs driven by simplified monochromatic wave fields and perfectly resonant wind stress, that is, a time-varying surface stress with period T = 2π/f in which f is the Coriolis parameter. A preliminary LES of a hurricane OBL with a simplified time-varying wave field is described by Sullivan and McWilliams (2010). There are numerous studies of high-wind OBLs carried out using 1D single-column parameterizations (e.g., Price et al. 1986; Large et al. 1994; Crawford and Large 1996; Zedler et al. 2002). These investigations emphasize the importance of inertially rotating wind stress for mixing at the thermocline but do not explicitly include surface wave effects; wave effects may be indirectly included in these models through empirically tuned constants and mixing rules. The computational elements of the present problem are a spectral wave prediction model WaveWatch III (Tolman 2002; Romero and Melville 2010b) and an LES model of the OBL with wave effects (Sullivan et al. 2007). The manuscript is arranged as follows: section 2 describes the large-scale hurricane wind and wave fields; sections 3 and 4 provide an overview of the problem posing and discuss the external wind and wave forcings applied to the LES model; sections 5 and 6 are devoted to a discussion of the results; and section 7 is a summary of the findings.

2. Hurricane wind and wave fields

Present computer power is insufficient to simultaneously resolve all dynamical scales in a hurricane-driven ocean basin where the largest scales of motion are hundreds of kilometers in horizontal directions while the smallest scales associated with Langmuir turbulence are about one meter. To make our small-scale turbulence computations tractable, we utilize a simplified one-way coupling strategy where the wave fields are first computed over a large ocean basin using surface winds that are an idealization of a translating hurricane vortex. In the present study, the wind fields come from hindcast products and are an idealization of Hurricane Frances. The winds are first used to drive the wave model, and then the Stokes drift, computed from the wave spectra, and the surface wind field are imposed as time and space varying external point forcings to a LES model of the upper OBL. This one-way coupling allows us to examine how the OBL dynamics, including Langmuir turbulence, respond to rapidly varying large-scale forcings. This coupling strategy with LES neglects some large-scale processes, namely, Ekman pumping and mean-field advection, which are most important for the evolution of the spatially 3D hurricane wake at late times (Price 1981; Sanford et al. 2011).

Hurricane Frances was a category 4 hurricane that developed in the Atlantic basin in 2004 and was one of the most heavily studied storms in the Coupled Boundary Layer Air-Sea Transfer (CBLAST) program (Black et al. 2007). Novel profiling floats developed by Sanford et al. (2007), D’Asaro et al. (2007), and Sanford et al. (2011) document the upper-ocean mixing induced by Frances over several days including the development of the hurricane cold wake. The impact of the storm on the large-scale currents and scalars is modeled by Zedler (2007) and Sanford et al. (2011).

a. Winds and surface fluxes

Surface winds, momentum, and heat flux parameterizations are needed to drive the LES and also the wave model. We make particular choices for our modeling but appreciate that there is uncertainty in specifying these external forcings owing to the variability in the available observations. An acceptable representation of the Hurricane Frances winds is derived from the NOAA hindcast data products described by Zedler (2007, see appendix). The specific formulas for the tangential and inflow radial winds (υφ, υr) expressed in terms of the radius r are
e1
which apply over the domain rmaxr ≤ 20 rmax, where rmax = 40 km is the radius of maximum winds. As suggested by Zedler (2007), these expressions are augmented with a decaying exponential function for rrmax. In (1), the winds are in units of meters per second, r is in units of kilometer, and the wind υr is positive inward. Based on the measurements of Sanford et al. (2007), we also add a constant translation speed VH = 5.5 m s−1 to (1), which results in an asymmetrical storm.
The parameterization of the surface stress τ at high winds is a key atmosphere–ocean coupling variable in hurricane prediction. There is vigorous debate concerning its parameterization but, as is customary, we adopt the bulk aerodynamic formula
e2
where ρa is the air density and U10 = (U10, V10) is the surface wind at a reference height za = 10 m. The drag coefficient is parameterized with the Large and Pond (1981) prescription
e3
This formula likely overpredicts the drag coefficient at very low winds |U10| ~ 5 m s−1 but closely follows the observations at wind speeds 10 ≤ |U10| ≤ 20 m s−1. Based on a growing body of experimental data collected in the field and laboratory (Powell et al. 2003; Donelan et al. 2004; French et al. 2007) we further limit (saturate) the maximum Cd = 1.8 × 10−3 at wind speeds beyond 20 m s−1. This saturated drag coefficient is near the estimates obtained from field observations but is lower than the value found by Donelan et al. (2004), Cd ~ 2.3 × 10−3, and Bell et al. (2012, see Fig. 19) who report a mean Cd ~ 2.4 × 10−3 above 50 m s−1 based on a budget analysis in the atmosphere. Sanford et al. (2011, see Fig. 7), using an integral momentum balance in the ocean, finds smaller values Cd ~ 1.5 × 10−3 in the wind speed range from 30 to 50 m s−1.
Surface latent Qlat and sensible Qsen fluxes are also imposed at the water surface in the LES. The parameterization of these fluxes is
e4
where is the specific heat, is the latent heat of vaporization, (ΔT, Δq) are the temperature and specific humidity differences between the sea surface and the reference height, and (Ch, Ce) are exchange coefficients. We choose Ch = Ce = 1.5 × 10−3, ΔT = 2.5 K, and compute Δq based on a relative humidity of 80% at the reference height and assuming saturation at the sea surface. These choices are representative values used in hurricane models (e.g., Price 1981; Jacob et al. 2000; Jacob and Shay 2003; D’Asaro 2003; Zedler et al. 2009) and result in a total heat flux from the ocean to the atmosphere of more than 1000 W m−2 at the time of maximum winds. Despite this seemingly high value the surface heat flux is ~0.1 of the dominant entrainment flux at the time of maximum winds. The sensitivity of the ocean boundary layer evolution to the surface exchange coefficients is described by Zedler et al. (2009).

b. Wave field

The surface wind and wave fields input to the LES are products derived from the wave prediction model WaveWatch III (Tolman 2002) adapted to hurricane conditions (Romero and Melville 2010b; L. Romero et al. 2012, manuscript in preparation). The hurricane wave simulations use nearly exact computations of the nonlinear energy transfers (van Vledder 2006), wind input by Snyder et al. (1981), and the energy dissipation by Romero and Melville (2010b) with the following modifications. For this study, the dissipation is improved to allow a smoother development of the spectrum with multiple peaks and to allow the spectral tail to follow the rapidly turning winds, while matching the empirical degree of saturation with weak dependence on the wind speed or friction velocity (Banner et al. 1989; Romero and Melville 2010a). A description of the code and its validation for high winds and hurricanes is given by Romero and Melville (2010b) and Romero et al. (2012).

For the present application, the wind–wave simulations of Hurricane Frances are carried out in an XY domain (±410, ±810) km with a horizontal resolution of 10 km, that is, 83 × 163 grid points. The spectral grid is discretized with a directional resolution of 10° using a constant relative bandwidth /ϖ = 0.078; ϖ is the radial frequency. The range of resolved frequencies is [0.157, 15.341] rad s−1 corresponding to the wavenumber range [0.0025, 24] rad m−1. The storm eye is initially placed at (X, Y) = (0, −510) km and is translated in time. Time series of winds and waves are collected at 61 discrete points distributed along the line (X, Y) = (±300, 0) km. The primary output of the wave model is the 2D wave height spectrum F(ϖ, φ), which varies with ϖ and angular position φ. The wave spectra are postprocessed to generate the Stokes drift velocity and wave age fields.

The Xt Hovmöller contour maps presented in Figs. 1 and 2 illustrate some of the variability in the wind and wave fields induced by our hypothetical storm over the horizontal domain X = ±300 km for a time period of more than 50 h. In these and subsequent figures, time is referenced to the time of maximum winds tm on the right-hand side of the storm. Here tm = 45.625 h is based on a constant storm translation speed Vh = 5.5 m s−1 and the initial storm center location. Thus, the cross section (X, Y = 0) is ahead of the storm for ttm < 0, cuts through the storm eye at ttm = 0, and is behind the storm (in the storm wake) for ttm > 0. Inspection of the contours in Fig. 1 shows that the wind speed |U10| is asymmetrical about X = 0 km, a consequence of the translation speed Vh, with the largest values ~50 m s−1 on the right-hand side of the storm track.

Fig. 1.
Fig. 1.

The hurricane storm propagates from left to right, and the fields are shown in the form of an Xt Hovmöller diagram; that is, the temporal variation of a particular field is shown for a chosen X cross section. Time is referenced to the time of maximum winds tm. Wind speed |U10| (color scale m s−1) imposed in the WaveWatch III simulations and (bottom) the wave age A = Cp/|U10| over the lifetime of the storm.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Fig. 2.
Fig. 2.

Stokes velocity components (top) us and (bottom) υs at the water surface for the imposed wind field (color scale units m s−1).

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Wave age, shown in the lower panel of Fig. 1, is defined as A = Cp/|U10| where Cp is the phase speed of the peak in the wave height spectrum: A varies widely in space and time across the storm path. Extensive periods of growing waves A ~ 0.5 are observed near the time of maximum winds with the lowest wave age on the right-hand side of the storm. For time periods well ahead of the storm center, A ~ 1.2 indicative of near wind–wave equilibrium conditions. At ttm ~ −22 h, the wave age jumps abruptly, A > 2, as fast-moving swell propagates rapidly ahead of the storm. In the wake of the storm ttm > 0, the wave age again slowly approaches wind wave equilibrium and at long times ttm > 20 h the wave state again becomes swell dominated. These estimates of wave age in a hurricane are consistent with values deduced from observations of 1D spectra documented by Young (1998, 2006) who reports values ranging from [0.5 − 2.0].

A key variable for our LES is the Stokes drift velocity us, which appears in the vortex force where is the resolved vorticity; see (8). For a broadband wave spectrum, us is obtained by the integration (Kenyon 1969):
e5
where k denotes the horizontal wavenumber vector [also see McWilliams and Restrepo (1999)]. The Stokes drift velocity us inherits the complexity of the wave field through (5); it varies with wind speed, spatial position, and time as depicted in Fig. 2. For our particular choice of storm propagation direction, the largest Stokes drift occurs in component υs on the right-hand side of the storm track. The distribution of the wave field is more asymmetrical than the wind field mainly due to the “extended” fetch on the right-hand side of a translating hurricane (Young 2003). The asymmetrical distribution of υs about X = 0 km reflects the stronger winds and better wind wave alignment on the storm right hand side. Clearly us is largest (most negative) in front of the storm center with contours roughly symmetrical about X = 0 km.
McWilliams et al. (1997) introduced a surface based turbulent Langmuir number
e6
as a measure of the relative importance of shear forcing to vortex force. In (6) the water friction velocity , where ρw is the density of water, us is evaluated at the water surface, and τ is given by (2). Grant and Belcher (2009) also adopt (6) in their scaling of the turbulent kinetic energy equation with wave effects. Harcourt and D’Asaro (2008) prefer a surface based Langmuir number , where 〈usSL is the average Stokes drift in the surface layer. These definitions of Langmuir number are, however, strictly only applicable in the situation of aligned winds and waves; Van Roekle et al. (2012) attempts to generalize (6) for wind–wave misalignment. Here we adopt (6) as our working definition of Langmuir number because of its simplicity and the ability to compare with previous studies; note that in (6) us is computed from the full wave spectrum according to (5).

Figure 3 shows contours of turbulent Langmuir number for our hypothetical storm. Outside the radius of maximum winds there are broad regions where 0.25 < Lat < 0.4. The very low values in the storm eye should be favorable for generating Langmuir cells with high aspect ratio, i.e., with their streamwise-to-crosswise length scale much greater than unity (Chini et al. 2009). Historically the criterion Lat < 0.4 (e.g., Li et al. 2005), used for aligned winds and waves, implies that often the hurricane driven OBL should be in a regime favorable to the generation of Langmuir turbulence. However, the applicability of the criterion Lat < 0.4 for nonequilibrium situations where weaker Langmuir cells can develop off-axis between the wind and wave directions is unknown (Van Roekle et al. 2012).

Fig. 3.
Fig. 3.

(top) The turbulent Langmuir number Lat (nondimensional color scale) and (bottom) Stokes depth scale Ds. In the bottom panel the color bar is in units of meters.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

The vertical profile of Stokes drift us(z) spirals with depth, a consequence of a multicomponent wave field composed of swell and wind waves that propagate in different directions. To characterize the depth penetration by the Stokes drift, an important variable in Langmuir turbulence, we define a vertical length scale
e7
where the integral is taken over the interval zs < z < 0 with zs chosen as the depth where the integrand remains positive [see Harcourt and D’Asaro (2008) for an alternate length scale for aligned winds and waves]. The above definition is robust and reduces to Ds = λ = 2π/k for a monochromatic wave field. The variation of Ds(X, t) is given in Fig. 3. Over the time interval (−20, 10) h, Ds ~ 60 m with values greater than 130 m on the right-hand side of the storm at the time of maximum winds. At ttm ~ 5 h, there is a small region to the left of the storm eye where Ds becomes small, a reflection of the confused sea state induced by opposing winds and waves. This is also a region where Lat > 0.6. Based on the results in Fig. 3 we conclude that the vertical profile of Stokes drift is track dependent and differs noticeably from a monochromatic prescription, but occupies a significant fraction of the OBL even at the time of maximum winds; that is, Ds/|h| > 1 where h is the OBL depth.

3. LES problem posing

The structure of the hurricane driven OBL depends on the large-scale forcings discussed in section 2. To thoroughly investigate all possible combinations of wind and wave forcing is beyond our computational reach with the LES and, thus, we restrict our attention to two locations X = ±60 km from the storm center shown in Fig. 1. These locations are of interest because they feature high wind conditions with different wind stress time histories and have a mixture of differently oriented winds and waves. As a further parameter variation, on each side of the storm, simulations are also conducted with and without surface wave effects, that is, with and without Stokes drift which generates the vortex force.

a. LES model

The details of our LES model for the OBL with wave effects are extensively discussed by McWilliams et al. (1997) and Sullivan et al. (2007). The equations for the resolved flow components are the Craik–Leibovich (CL) theory with crucial wave–current coupling through the vortex force. The LES model, with wave effects, is based on the incompressible Boussinesq equations with a single-point, second-moment turbulent kinetic energy (TKE) closure subgrid-scale parameterization and a flat ocean surface (Moeng 1984; McWilliams et al. 1997). The added wave effects are the vortex force, Stokes–Coriolis coupling, Lagrangian mean advection associated with Stokes drift, a wave-averaged increment to the pressure that arises through conservative wave–current interaction (McWilliams et al. 1997), and Stokes production in the subgrid-scale (SGS) TKE equation (Sullivan et al. 2007). Additional acceleration and energy generation due to nonconservative stochastic wave breaking, described by Sullivan et al. (2004, 2007), are not considered here.

We use standard LES notation where spatially filtered (or resolved) variables are denoted with an overbar: are the resolved-scale velocity and vorticity; e is the subgrid-scale TKE; is the temperature; and π is a generalized pressure variable. Other variables appearing in the LES equation set are the Coriolis frequency f = (0, 0, f); reference temperature ; gravitational acceleration g; and subgrid-scale momentum and density fluxes, (τij, τ), respectively. These subgrid-scale fluxes are modeled using the eddy viscosity prescription described in Moeng (1984) and Sullivan et al. (1994), which implies that the principal feedback of e on and occurs through subgrid-scale mixing with eddy viscosity and diffusivities (νt, νs) ∝ e1/2.

In the situation of misaligned winds and waves, Stokes drift velocity couples to the resolved vorticity to generate a three-dimensional vortex force
e8
that contributes to all three resolved momentum equations; in the case of aligned winds and waves, which is most often studied, the vortex force reduces to the two-component form us(0, −ζ, η).

b. LES setup

Two computational domains of varying size are used in the LES. First, the vertical dimension of the computational box is chosen to be about twice the expected maximum depth of the OBL; based on observations |h|max ~ 120 m (Sanford et al. 2011) and thus we choose a box of depth −240 m. The horizontal dimension is picked to be ~5|h|max so as to allow the interior motions to develop independent of the periodic sidewall boundary conditions (see Schmidt and Schumann 1989). On the left-hand side X = −60 km, which has slightly weaker forcing, a small domain (Lx, Ly, Lz) = (750, 750, −240) m was deemed acceptable. On the right-hand side of the storm, X = 60 km, the vigorous forcing from the inertially rotating wind stress was found to generate a system of large-scale random gravity waves in the hurricane wake. Test runs indicated that these structures were only adequately captured by increasing the horizontal area of the small domain LES by a factor of 4. Thus, a large domain is used on the right hand side with (Lx, Ly, Lz) = (1500, 1500, −240) m. The number of grid points in the small domain is (Nx, Ny, Nz) = (512, 512, 256) and is increased to (1024, 1024, 256) in the large domain; thus the horizontal spacing Δx = Δy = 1.46 m in all the computational boxes. In the vertical direction, we use a stretched vertical grid where the spacing between any neighboring cells is held fixed at the ratio K = 1.00451 with the first w grid point located at Δz1 = −0.5 m; thus the first u point is at z1/2 = −0.25 m in a staggered grid. Based on these choices, the cell aspect ratio Δxz = (2.93, 1.54) at z = (0, −100) m. The small and large computational boxes are able to capture small-scale turbulence at early times and larger scales of motion that develop during the rapid deepening of the OBL at the time of maximum winds.

The initial temperature sounding is simply constructed from two linear segments: over the depth hiz ≤ 0, θi = 302.4 K and for z < hi the temperature decreases at a rate of 0.04 K m−1. This is approximately the variation observed for Hurricane Frances (Sanford et al. 2011). The depth hi = −32 m is the initial OBL depth. The Coriolis parameter is f = 0.681 × 10−4 s−1, corresponding to an inertial period T = 25.6 h. The simulations are initiated from rest with small random temperature perturbations under the wind and wave forcing from hour 10 of the WaveWatch III simulation. At this time, the wind and wave fields are slowly varying with winds less than 10 m s−1. Next, these external forcings are held constant and the simulations are integrated for 10 physical hours to generate fully developed equilibrium turbulence. The resulting currents are ~ 0.05 m s−1, and the initial state of these weak currents has minimal impact on the long time ocean response under strong forcing (Crawford and Large 1996, p. 889). Finally, the winds and Stokes drift are allowed to follow the time history shown in Figs. 5 and 6. The time step Δt in the code is dynamically picked to obey a fixed Courant–Fredrichs–Lewy condition. Here Δt varies by an order of magnitude, from 0.3 to 5 s, and more than 300 000 time steps are needed to cover the 70 physical hours of the simulation. With a storm translation speed of 5.5 m s−1 the LES domains are initially located more than 900 km from the storm center. The parallelization of the code is described by Sullivan and Patton (2011).

4. Transient wind and wave forcing applied to LES

The temporal and spatial variability in the large-scale forcing conditions critically impacts the time evolution of the OBL and the development of Langmuir turbulence. To aid the interpretation of the OBL dynamics generated by the simulations we first discuss the time history of winds, waves, and cooling that are applied at the upper boundary of the LES; see Figs. 4, 5, 6, and 7. In these figures, the winds U10, surface cooling flux Q* = Qsen + Qlat, Stokes drift velocity us, wind–wave alignment, and surface wind stress τ are time slices at the horizontal locations X = ±60 km taken from Figs. 1 and 2.

Fig. 4.
Fig. 4.

External winds (U10, V10) and cooling Q* applied to the LES: (a) right-hand side of the storm X = 60 km and (b) left-hand side of the storm X = −60 km.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Fig. 5.
Fig. 5.

Surface wind stress and wave forcing on the right hand side of the storm, X = 60 km: (a) time variation of the wind stress τ where the (black, gray) lines are the components τ = (τx, τy), respectively. The (black, gray) dotted lines are resonant stress components τi(t) that rotate clockwise with inertial frequency f; that is, τi(t) = |τ(t)|eift. (b) The surface Stokes drift velocity components (us, υs) normalized by the local friction velocity u*(X, t) are depicted. (c) The wind–wave alignment given by (9) at the water surface (black line) and at depth z = −(5.8, 12.4, 40, 80.2) m gray (solid, dotted, dash-dot, long dash) lines, respectively.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Fig. 6.
Fig. 6.

Surface wind stress and wave forcing on the left-hand side of the storm, X = −60 km. The panels, labeling and variables are the same as in Fig. 5. Note the difference in the vertical scale in the lower panels of Figs. 6 and 5.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Fig. 7.
Fig. 7.

Vertical profiles of Stokes drift velocity normalized by the local friction velocity (us, υs)/u* at time instances ttm = (−25, −10, 0, 5) h: z is made dimensionless by the Stokes depth scale Ds given by (7). The solid (black, gray) lines are (us, υs)/u*, respectively, along the line X = 60 km. The dashed (black, gray) lines are their counterparts along the line X = −60 km. For the selected time instances the depth scale Ds = (14.5, 86.0, 123.5, 74.7) m and (12.6, 75.2, 45.7, 37.1) m for X = 60 (−60) km. Note at ttm = 5 h the solid black and gray lines overlap.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

The primary forcing to the LES (and also the wave model) are the surface winds, shown in Fig. 4, converted to a surface wind stress according to the prescription (2). The amplitude of the wind stress, depicted in Figs. 5 and 6, is a strong function of (X, t) within the storm track and, as expected, each component displays a strong pulse at the time of maximum winds mimicking the surface winds. Figures 5 and 6 also compare the time varying stress τ(t) with an “inertial” stress of equal magnitude that rotates clockwise with frequency f; that is, τi(t) = |τ(t)|eift. This exposes the direction and rate of rotation relative to an inertial motion. Both wind stress components on the right side of the storm track are partially resonant with the ideal inertial stress τi. The period of resonance is relatively short, about 10 h, but occurs during the time of maximum winds (−5, 5) h. Meanwhile, on the left-hand side, the wind stress magnitudes are modestly lower than on the right-hand side and component τy clearly rotates in a direction opposite to that of an inertial rotation. Transient wind forcing at frequency f with clockwise rotation is able to drive currents below the mixed layer and thereby efficiently mix the OBL (e.g., Price 1981; Crawford and Large 1996; Skyllingstad et al. 1999). The observations of Sanford et al. (2011) demonstrate how wind stress with limited resonance time is able to efficiently mix a hurricane OBL. At early and, more so, at late times the wind stress rotation applied at the upper boundary of the LES is noninertial on both sides of the storm track.

Surface cooling plays a minor role, despite the more than 1200 W m−2 applied at the time of maximum winds; it is offset by the rapid increase in surface shear stress. The ratio |hi|/L where is the Monin–Obukhov length is often used in atmospheric boundary layer analysis as a measure of the relative importance of shear and buoyancy forcing at the surface (e.g., Deardorff 1972; Moeng and Sullivan 1994). In this expression, the buoyancy parameter β = g/θo and κ = 0.4 is the von Kármán constant. Broadly, over the bulk of the time period the OBL surface shear forcing dominates surface cooling since −1 < |hi|/L < 0 and during the time period of maximum winds, (−10, 10) h, we find −0.1 < |hi|/L < 0. Cooling plays a modest role at early times and contributes to the well-mixed (vertically uniform) shape of the mean current profiles.

It is informative to discuss how the Stokes velocity varies in our simulations since it impacts the development of Langmuir turbulence. In the lower and middle panels of Figs. 5 and 6 we show the time variation of the wind–Stokes alignment
e9
at depths z = (0, −5.8, –12.4, −40, −80.7) m, and the components of the normalized Stokes drift velocity at the surface (us, υs)/u*. The surface based Lat, given by (6), can be inferred from this plot and its range is narrow; along X = 60 km, 0.28 < Lat < 0.32 and along X = −60 km, 0.28 < Lat < 0.41. First, we notice us(0, t)/u* has the same shape and amplitude on both sides of the storm track, X = ±60 km. This is expected based on the contours of Stokes drift velocity in Fig. 2. Thus, the important differences in us/u* between the right-hand and left-hand sides are due to component υs/u*. On the right-hand side of the track υs(0, t)/u* is nearly constant with time, less than a 25% variation over the period (−30, 5) h. This correlation between Stokes drift and wind stress is meaningful since u* varies by a factor of 7 over the time period of the simulations. Meanwhile on the left-hand side υs(0, t)/u* at first decreases linearly, then changes sign and finally reaches a minimum near ttm = 0 h. These differences are a result of the extended fetch on the right side of the storm versus the progressive misalignment between winds and waves on the left side as time advances.

The wind–Stokes alignment is higher at the surface than at depth, which results from a multimode surface wave spectrum composed of wind waves and swell: recall in wind wave equilibrium cosφ = 1 at all z. Short wavelength wind waves rapidly align themselves to varying surface shear forcing compared to the slow adjustment of long wavelength swell; this pattern of wave development in a hurricane is reported by Wright et al. (2001). Thus surface values of Stokes drift are only a partial reflection of the wave state that drives Langmuir turbulence. Swell first appears at ttm = −22 h, as indicated by wave age A > 1.2 in the lower panel of Fig. 1, and reduces the wind–Stokes alignment at depth throughout the passage of the storm. At the surface, on the right-hand side of the storm the wind–Stokes alignment is always within 10° at the surface and no more than 45° at z = −5.8 m. At the time of maximum winds the Stokes drift is near “resonant” with the winds; that is, the wind–Stokes alignment is near unity over a depth of more than 80 m and the Stokes depth scale exceeds 130 m. Thus we anticipate strong Langmuir turbulence to develop at this time. Over the time period (−5, 5) h, the wind–Stokes alignment on the left-hand side increases, 60° at the surface and almost 90° at z = −10.2 m, because the wind and dominant wave propagation directions are opposed.

Selected vertical profiles of the Stokes drift velocity are shown in Fig. 7 at times before, at, and after the time of maximum winds. These profiles illustrate that the normalized surface values of us/u* and, hence, Lat, have only a mild variation over the time period. The invariance of Lat however disguises the profile shape changes that occur with (X, t). For example, the profiles of υs(z, t) at ttm = (−25, 0) h have similar surface values of Lat but the profile at ttm = 0 extends much deeper into the OBL: at ttm = (−25, 0, 0) h Ds = (14.5, 123) m, respectively. Intuitively, we expect the formation of Langmuir turbulence to depend on Lat, the vertical shape of the Stokes profile and its depth penetration.

5. Identification of Langmuir turbulence

a. Temporal evolution of vertical velocity

The LES of steady aligned wind wave flows and a few field observations find vertical velocity is a robust signature of Langmuir turbulence. In the Langmuir turbulence regime, coherent Langmuir cells form and these structures can be detected in the instantaneous vertical velocity field and furthermore impact the statistical moments of w. Figure 8 compares the maximum total (resolved plus SGS) second-order moment 〈w2max and the minimum resolved third-order moment diagnosed from our simulations. Since the LES equations are closed at the second-moment level there is no estimate of third-order SGS quantities. Both moments are normalized by local u*(X, t). These temporally varying moments are obtained by first computing the vertical profiles of resolved and at all z for a given t: angle brackets denote an xy spatial average over the LES domain. The and are found by searching for the appropriate maximum or minimum along their respective vertical profiles. We then add an estimate of the SGS contribution, 2e/3, to the resolved w variance; that is, 〈w2max is the sum of resolved and SGS contributions. We note that the vertical locations of 〈w2max and are not the same and vary with time. These positions are always located in the upper surface layer of the OBL but are found deeper in the water as the winds and waves increase.

Fig. 8.
Fig. 8.

Temporal variation of the maximum total vertical velocity variance 〈w2max and minimum resolved third-order moment for cases with and without wave effects at locations left and right of the storm center. Simulations with wave effects for (a) X = 60 km and (b) X = −60 km and simulations with no wave effects for (c) X = 60 km and (d) X = −60 km. The moments are normalized appropriately by either or .

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

The time series depicted in Fig. 8 are compelling evidence for Langmuir turbulence in the hurricane driven OBL. The w moments are sensitive to the location within the storm track, depend on the time and space evolving wind and wave fields, and the inclusion of phase averaged wave–current interactions in the LES. Over the lifetime of the vortex-force simulation 〈w2max > 2.5 along the track X = 60 km. Notice also that during the time period with resonant wind forcing, (−5, 5) h, the maximum vertical velocity variance is near 3 and this level persists in the hurricane wake. Recall u* varies by about a factor of 7 and Ds varies from 20 to 130 m over this time period and, thus, the dynamic range of 〈w2〉 is broad. Surprisingly, simple scaling by u* captures much of the variation in w variance. On the left side of the storm track, the maximum w variance at first tracks its counterpart on the right-hand side, then gradually decays, reaches a minimum value at the time of maximum winds, and gradually recovers in the hurricane wake. In the simulations without vortex force, the maximum w variance is nearly independent of the track position and is similar to what is expected for a neutral flat-wall boundary layer flow.

The differences between the variations on the right and left hand side of the storm can be traced to the wave fields, the wind–Stokes alignment shown in Figs. 5 and 6, and the production of w variance by Stokes drift. The Stokes contribution to the ensemble average w variance budget consists of two terms
e10
where 〈uw〉 and 〈uw〉 are ensemble average (total) vertical momentum fluxes. Note there is no mean shear production in the w budget. The derivation of (10) follows the same sequence of steps used to derive the TKE budget with wave effects (e.g., McWilliams et al. 1997). The total TKE budget depends on the vertical gradients of the currents and Stokes drift, i.e., the mean Lagrangian shear
e11
Notice that in (10) the production term s depends on the correlation between turbulent fluxes and Stokes drift gradients. These terms can be positively correlated as in the case of wind wave equilibrium, uncorrelated or even negatively correlated: that is, the momentum fluxes can be countergradient to the Stokes drift gradient. A time series of Stokes production, given by (10), at the location of the maximum vertical velocity variance is shown in Fig. 9. There is a significant temporal correlation between the enhancement and reduction in w variance, shown in Fig. 8, with the time variation of Stokes production. In particular the variance on the right and left hand sides of the storm track first begin to depart at ttm = −12 h. At this time on the left-hand side of the storm, the production term reaches a positive maximum and switches sign and becomes a destruction term. An even steeper reduction in w variance occurs at the time of maximum winds where is a sink of w variance and is small and positive. On the right-hand side of the storm are always positive with reaching a maximum in the storm wake. Thus Stokes production is most often a significant source of w variance, but can in certain misaligned wind wave states be a sink. Wave–current interactions, which are key processes in the budgets of turbulence, play a role in the enhancement and reduction of w variance in the hurricane OBL. The results emphasize the importance of the vertical gradients of the Stokes drift profile and their alignment with the vertical momentum fluxes. It is worth noting that < 0 also implies that the currents and waves are moving in opposite directions near the surface (see section 6b), which can enhance wave breaking (Melville 1996).
Fig. 9.
Fig. 9.

Production of vertical velocity variance by Stokes drift. The terms and are computed at the vertical location of maximum vertical velocity variance along the right-hand side (solid lines) and left-hand side (dotted lines) of the storm track: are indicated by black (gray) line colors. The terms are made dimensionless by and are obtained by spatial averaging over an xy plane and by time averaging over 45 min. The momentum fluxes, which appear in , are the sum of the LES resolved and SGS contributions. Note the periods where are negative on the left-hand side of the storm track indicating that the turbulent fluxes and gradients of Stokes drift are countergradient.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

D’Asaro (2001), using drifting Lagrangian floats, reports measurements of vertical velocity variance that are “remarkably well parameterized by wind stress” with values ranging from 1.75 to 2 times that measured in solid-wall turbulent layers driven by shear. Our vortex-force simulations on the right hand side of the storm track are in good agreement with these measurements, and we speculate that the observed high correlation with wind stress is a consequence of the tight relationship between surface Stokes drift and u* for approximately aligned winds and waves, as shown in Fig. 5. Harcourt and D’Asaro (2008) and Grant and Belcher (2009) propose a velocity scale for Langmuir turbulence, . However our range of Lat is somewhat limited along X = 60 km because of the tight correlation between us and u*; hence along this track wL ~ u*. We mention that the applied wind forcing in our simulations varies smoothly with time (there is no assumption of wind wave equilibrium) and then the wave fields adjust naturally to the changing local wind conditions. Hence the Stokes drift velocity imposed in the LES also varies smoothly with time and changing winds.

The impact of Langmuir turbulence on the third-order vertical velocity moment is impressive. Figure 8 depicts the minimum based on the resolved vertical velocity. With vortex force the minimum third-order moment can be 10 times larger than its no vortex force counterpart at the time of maximum winds. Again we notice a significant reduction in on the left-hand side of the storm in the simulation that includes the vortex force mirroring the reduction in 〈w2max. The third-order moment of vertical velocity indicates an abundance of Langmuir cells, and is important as it appears in the turbulent kinetic energy equation. The changes in indicate that Langmuir turbulence has modified the vertical transport, the divergence , in the OBL in significant ways (Grant and Belcher 2009). We find that the moments 〈w2max and are stronger and clearer indicators of Langmuir circulations than the nondimensional skewness , which is often used to describe Langmuir turbulence (e.g., Li and Garrett 1993; McWilliams et al. 1997; Sullivan et al. 2007). In our simulation with vortex force on the right hand side of the storm minimum skewness varies over the range −1.1 < Sw < −0.6, which indicates that the downwelling velocities are narrower and stronger compared to the broader upward motions.

b. Flow structures

Flow visualization of the vertical velocity field provides insight into the formation and evolution of the coherent structures in a high-wind hurricane driven OBL. Extensive interrogation of the LES volumes shows that elongated streaky structures, that is, coherent Langmuir cells, are pervasive at all stages of the storm evolution. The elliptical shaped cells meander randomly in their long-axis direction with occasional merges at downstream junctions. Figure 10 is a typical sample that illustrates these processes. This image compares the instantaneous normalized vertical velocity at the time of maximum winds on the left and right hand sides of the storm. Also, included in this image are the directions of the surface and Stokes drift velocity, which provide a reference for the cell alignment. To allow a fair comparison of the spatial scales, the visualization is depicted over equal subareas of the computational domain. Notice the impressive size and strength of the cells on the right hand side of the storm. At times their long-axis is more than 100 m with pockets of downwelling velocity exceeding 5u* ~ 0.35 m s−1: similar large values were observed by vertical profiling floats on the right hand side of Hurricane Frances (T. B. Sanford 2011, personal communication). Meanwhile, on the left hand side coherent cells can be identified but they are clearly weaker in magnitude and smaller in scale. This flow visualization is consistent with the temporal evolution of vertical velocity variance and third-order moment discussed previously.

Fig. 10.
Fig. 10.

Contours of vertical velocity normalized by u* at the time of maximum winds at (top) z = −7.2 m along the track X = −60 km and (bottom) at z = −12.6 m along the track X = 60 km. The subarea depicted in the top (bottom) panels is 44 (11) % of the total xy area in their respective computational domains. The black (red) vectors are the directions of the surface wind (Stokes drift).

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Instantaneous images of w illustrate the chaotic interaction between cells and turbulence. However, in order to more clearly identify the average coherent structure (a Langmuir cell) in our simulations we employ the technique of linear stochastic estimation (LSE) (e.g., Adrian et al. 1989; Adrian 1996; Christensen and Adrian 2001). LSE is a robust technique for identifying structures in turbulent flows and is an excellent approximation of a conditional average that utilizes only unconditional two-point correlation data. For our flows, we estimate the deviation from the mean horizontal currents, that is, the horizontal perturbation velocities at x′ = x + Δx where Δx is the spatial displacement, by the conditional average
e12
We choose resolved vertical velocity at the vertical location of 〈w2max, that is, z = zmax, as the event indicator. LSE utilizes unconditional two-point spatial correlations
JPO-D-12-025.1-e13a
JPO-D-12-025.1-e13b
where σ denotes the root-mean square of the given quantity. In our implementation of (12) and (13) with LES data the angle brackets denote an average over an xy plane for each particular point (x′, y′).

Time sequences of the conditionally averaged field at z = zmax and contours of the normalized two-point spatial correlation Rww(x′, y′, zmax) are provided in Figs. 11 and 12. We define the spatial extent of the coherent Langmuir structure by contours of Rww; the positive contour level 0.1 encompasses an ellipsoidal region of strong downwelling with side lobes of broad weaker upwelling, contour level −0.05. Inside a downwelling region, the perturbation velocity field is observed to be closely aligned (parallel) to the major axis of the elliptical shaped structure. The coherent structure defined by the two-point correlation function Rww and the accompanying perturbation velocity fields are a consistent average representation of the instantaneous flow fields shown previously. The flow patterns are in good agreement with the findings of McWilliams et al. (1997) and Sullivan et al. (2007), and are consistent with the CL2 instability mechanism described by Leibovich (1983). Two of the most prominent features of the cells in the hurricane OBL are their rotation and scale expansion and contraction for varying time and track position. The flow patterns exhibit strong horizontal convergence perpendicular to the major axis of the cell center. There is vigorous inflow at the cell tail and weaker outflow at the cell head with the strongest perturbation velocities centered just above the largest downwelling motions. On the left-hand side, the cells rotate continuously in counterclockwise fashion. Close inspection of the flow pattern at ttm ~ 0 shows the development of a vertical rotational motion that displaces the perturbation velocity vectors slightly off center. Van Roekle et al. (2012), in idealized LES with misaligned winds and waves, finds a similar result and attributes this to advection of horizontal vorticity across the cell. On the right hand side of the storm, the cells first rotate counter clockwise then rapidly clockwise. The cells grow rapidly in scale near the time of maximum winds as they tend to align themselves with the surface winds. In contrast, the cell growth is much more modest on the left hand side of the storm, consistent with the temporal changes in the vertical velocity moments. On the left hand side, the cells are larger at ttm = (−12.7, −3.44) h than they are at the time of maximum winds.

Fig. 11.
Fig. 11.

Normalized two-point spatial correlation Rww(x′, y′, zmax) along the track X = 60 km. The solid (dotted) contour levels are 0.1 (−0.05). The positive contour level encloses a region of strong downwelling while the negative contour level encloses a broad region of weaker upwelling. The vectors (arrows) show the magnitude and direction of the horizontal perturbation velocities (u′, υ′) induced by a Langmuir cell. (a)–(f) ttm = (−30.2, −12.5, −3.04, −0.003, 3.37, 12.8) h, repectively. Note the change in the horizontal and vertical scales between the lower and upper panels.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Fig. 12.
Fig. 12.

The normalized two-point spatial correlation Rww(x′, y′, zmax) and perturbation horizontal velocity (u′, υ′) along the track X = −60 km: contour levels and labeling as in Fig. 11; (a)–(f) correspond to ttm = ( −30.9, −12.7, −3.44, −0.003, 2.55, 10.6) h, respectively. Note how the Langmuir cells rotate counterclockwise in contrast to those in Fig. 11.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

c. Langmuir cell alignment

One of the important questions is the orientation of a Langmuir cell for misaligned winds and waves. Cell orientation significantly complicates the development of OBL parameterizations with wave effects (McWilliams et al. 2012). In Fig. 13, we show the surface wind and Stokes drift velocities and the direction of the Lagrangian shear given by (11). The orientation of the cells in LES is deduced from the conditionally averaged velocity fields sampled at the cell center [u′(0, 0, zmax), υ′(0, 0, zmax)]. The Langmuir shear s is a key quantity in Langmuir turbulence as it appears as a production term in the turbulent kinetic energy equation [see Eq. (5.1) from McWilliams et al. (1997) and Grant and Belcher (2009)]. In addition, Holm (1996) shows how the Lagrangian shear contributes to vortex stretching, that is, the tilting of vertical vorticity into streamwise directions. The results in Fig. 13 indicate that the cells are oriented primarily with s, which is also broadly parallel to the surface wind direction. The deviation from s is largest along X = 60 km in the wake where the cell alignment tends more toward the wind direction. Figure 14 shows how s evolves with time at the water surface. Notice that there are periods when Stokes and current shear are of opposite sign, which leads to a reduction in turbulence variances (see discussion in section 5a).

Fig. 13.
Fig. 13.

Time variation of the surface wind (dotted black line), the surface Stokes drift velocity (gray line), and the Lagrangian shear (black line): (left) X = −60 km and (right) X = 60 km.The dominant Langmuir cell direction deduced from the LES are indicated by solid circles. Angles are measured counterclockwise from the x axis.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Fig. 14.
Fig. 14.

Mean Lagrangian shear components at the water surface on (top) the left-hand side and (bottom) right-hand side of the storm. The shears are made dimensionless by z1/u*, where z1 is the location of the first LES gridpoint below the water surface; the heavy black lines indicate time periods where the current and Stokes shears have opposite signs.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

6. Impacts of Langmuir turbulence on OBL statistics

In the hurricane driven OBL, coherent Langmuir cells interacting with background turbulence leave an imprint on the vertical velocity. Langmuir turbulence also effects the OBL low-order statistics, namely, mean currents, temperature profiles, momentum and heat fluxes, variances and dissipation. This is an overly abundant amount of information to present, and thus we focus on a few select statistics that demonstrate the role of wave–current interactions.

a. Horizontal variances

In simulations with steady aligned winds and waves the coherent Langmuir cells modulate the horizontal streamwise and spanwise variances compared to their flat-wall counterparts. The streamwise oriented cells induce strong spanwise convergence at downwelling sites and their vortical rotation about the x axis amplifies the spanwise variance 〈υ2〉 at the expense of streamwise variance 〈u2〉. For example, McWilliams et al. (1997) reports that, based on the resolved LES velocities, the ratio near the water surface; theoretical linearized rapid distortion results for a propagating surface wave are in qualitative agreement with this numerical result (Teixeira and Belcher 2002). Wave breaking, however, tends to weaken this variance anisotropy by disrupting the near-surface cell structure (Noh et al. 2004; Sullivan et al. 2007). This intercomponent exchange of horizontal variance is also present in our simulations and reflects the directional changes in the time varying winds and waves as depicted in Fig. 15. As shown previously, cells tend to align with the time dependent Lagrangian shear direction and the rotating cells induce a smooth exchange of variance between the horizontal velocity components. The cells also noticeably alter the phase relationship between the resolved horizontal variances so that and reach maximum values at times markedly different than simulations without vortex force. This is especially apparent at the time of maximum winds on the right hand side of the storm track. In the absence of vortex force is a minimum at ttm = 0, while it reaches a maximum when vortex force is included. Note that at ttm = 0 the cells are primarily oriented in the y direction, which results in an amplification of the u variance with a reduction in the υ variance. Our flow visualization shows that the cells are weaker and smaller in scale on the left-hand side of the storm track and this results in a lower value of u variance compared to its counterpart on the right-hand side near the water surface.

Fig. 15.
Fig. 15.

Temporal variation of the resolved horizontal variances at z = −0.5 m below the water surface for simulations with and without wave effects at locations left and right of the storm center: simulations with wave effects for (a) X = 60 km and (b) X = −60 km and with no wave effects for (c) X = 60 km and (d) X = −60 km. As the winds and waves veer with time the horizontal component with maximum variance also shifts tracking the rotation of the coherent Langmuir cells. At ttm = 0 the winds are primarily aligned with the y direction and then the υ component is a maximum as might be expected for a wall-bounded shear flow without wave effects. The inclusion of vortex force alters this relationship so that the u component is a maximum.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

b. Current and variance vertical profiles

Vertical profiles of the mean currents highlight the complex coupling that occurs between winds and waves and between waves and currents. Figure 16 compares the mean currents and Stokes drift velocity, normalized by local u*, on the left and right sides of the storm track for simulations at the time of maximum winds. Wave influences are pronounced. On the right-hand side the Stokes drift velocity is dominated by the υs component with the wind–Stokes alignment near unity over a sizeable fraction of the OBL (see lower panel of Fig. 5). The potent Langmuir cells induce a vertically uniform υ-current profile down to a depth z/|h| ~ 0.6. In contrast, in the simulation without wave effects the vertical gradient is surface intensified but varies continuously over the entire OBL. In cases with and without wave effects, there are steep current gradients near the base of the OBL.

Fig. 16.
Fig. 16.

Profiles of mean currents and Stokes drift velocity 〈us, υs〉 normalized by u* at the time of maximum winds: left (right) sides of the storm are shown in left (right) panels. The black (gray) colors are simulations with (without) vortex force. Averages are computed over an xy plane and 3600 s in time.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Inspection of the current profiles shows that Stokes drift induces a return Eulerian flow opposite to the direction of wave propagation: this is a consequence of Stokes–Coriolis coupling, , in the LES equations (Huang 1979; Gnanadesikan and Weller 1995; McWilliams et al. 1997; Polton et al. 2005). On the right-hand side of the storm this process leads to Eulerian current changes that are similar to simulations with aligned winds and waves. However, on the left-hand side of the storm with wind–wave misalignment, the wave-induced return flow leads to oppositely signed gradients and ∂us/∂z. As discussed in section 5a, the vertical momentum flux 〈uw〉 is then countergradient to the gradient of Stokes drift, which inhibits the production of w variance. Over the bulk of the OBL the currents with and without wave effects are broadly similar with the important changes occurring near the water surface.

The resolved velocity variance profiles displayed in Fig. 17 paint a consistent picture of the impact of Langmuir turbulence in the nonequilibrium hurricane regime. On the left-hand side of the storm at ttm = 0, the w variance in the simulation with vortex force displays a slightly elevated subsurface maximum compared to its counterpart without vortex force. The wave effects on the horizontal variances are confined to a shallow layer near the surface −0.15 < z/|h| < 0, and outside this region the profiles from the simulations with and without vortex force are nearly indistinguishable. This is consistent with the flow visualization and the interpretation of weak Langmuir turbulence on the left hand side of the storm discussed in section 5. The variance profiles along X = 60 km are in sharp contrast to their counterparts along X = −60 km. The wave effects are pronounced and the changes in the υ and w variances extend over a large fraction of the OBL, more than 50%. Also, it is apparent that vertical transport by the Langmuir cells contributes to the elevated υ variance in the region −1 < z/|h| < −0.5. This hints at a coupling with shear instabilities at the thermocline induced by inertial resonance.

Fig. 17.
Fig. 17.

(left to right) Vertical profiles of resolved variances normalized by for (top) left side and (bottom) right side of the storm track. The black (gray) solid lines in each figure are simulations with (without) vortex force. Results are at ttm = 0.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

c. Entrainment flux and SST

One of the primary bulk functions of the hurricane OBL is to entrain cool water from below the thermocline into the mixed layer. Over time, this lowers SST that can, depending on the translation speed of the storm, feedback onto the large-scale hurricane dynamics. Previous modeling and observational studies (e.g., Price 1981; Crawford and Large 1996; Skyllingstad et al. 2000; Sanford et al. 2011) find that the primary dynamical mechanism for entraining cool water is shear instability at the base of the OBL induced by inertially rotating wind stress (see the upper panel of Fig. 5). This critical dynamical process is also found to be dominant in our simulations but is augmented by additional mixing induced by the coherent Langmuir cells that develop in the OBL on both sides of the storm track. This is illustrated in Fig. 18 where we compare the minimum total heat flux at the base of the OBL, that is, the entrainment flux 〈′〉min, and the change in SST, Δθ = 〈θ〉 − θi at z = 0, for simulations with and without wave effects. The entrainment flux 〈′〉min is the sum of LES resolved and SGS contributions and is simply computed by sweeping through the vertical profile of the mean total temperature (or heat) flux and locating its minimum value. Vertical profiles of mean temperature and heat flux are obtained by averaging over an xy plane in the LES domain and 1000 s in time.

Fig. 18.
Fig. 18.

Temporal variation of (top) entrainment flux and (bottom) SST change. The groups of curves marked “L” (“R”) are conditions on the left (right) hand side of the storm. In each figure the black (gray) lines are simulations with (without) vortex force. The entrainment flux is normalized by Q*(X, t) the local value of the surface heat flux (see Fig. 4) and Δθ = 〈θ〉 − θi at z = 0.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

Far ahead of the storm center the normalized entrainment flux 〈′〉min/Q*(X, t) ~ −0.25, which is slightly more negative than the classical value of −0.2 for a daytime convective atmospheric boundary layer (e.g., Schmidt and Schumann 1989; Moeng and Sullivan 1994). This is expected since the Stokes drift is modest and the OBL is in a regime of mixed shear and convective forcing with the ratio |hi|/L ~ −1. With increasing time, the OBL becomes increasingly dominated by shear and wave forcing. At the time of maximum winds, the entrainment flux exceeds its surface value by a factor of (−10, −3) on the (right, left) sides of the storm, respectively. Langmuir turbulence noticeably enhances the entrainment flux compared to the situation without wave effects. For example, on the right-hand side at ttm = −4 h the normalized entrainment flux increases from −4 to −6, a 50% increase when wave effects are included. At ttm = 0, the increase is about 20% and the resulting maximum change in Δθ is about 0.25 K cooler. In absolute values, the changes in entrainment flux and SST on the left hand side of the storm are more modest than their counterparts on the right hand side. On both sides of the storm track wave effects cause the minimum entrainment flux to occur earlier in time, perhaps by as much as 3 h. This suggests enhanced coupling with the inertial currents at the thermocline. Our estimate that the normalized entrainment (cooling) flux on the right-hand side of the storm is greater than 10 times the surface heat flux is in good agreement with the observations by Sanford et al. (2011, see p. 1044). The LES predicted changes in SST on the left and right hand sides of the storm are in good agreement with the observations of D’Asaro et al. (2007), but those observations do not readily lend themselves to quantifying the importance of surface waves.

In the hurricane wake the entrainment flux, especially on the right-hand side of the storm, appears to retreat toward its prestorm value in a series of steep rises and declines that is noticeably different than in the time period ttm < 0. Flow visualization reveals packets of large-scale propagating gravity waves that overturn intermittently, creating large oscillations in the entrainment flux as well as other OBL statistics. These turbulent gravity waves appear to be a product of the vigorous resonant wind stress forcing. The large horizontal scale of these structures motivated our use of the wide LES horizontal domain described in section 3b.

d. Boundary layer depth and Stokes scale

It is interesting to notice that the augmentation of the bulk temperature entrainment and SST change induced by Langmuir turbulence occurs during the time period ttm < 0. This can be interpreted in terms of the surface Lat and the variations of the bulk OBL length scales h and Ds. Figure 19 depicts the time growth of the OBL and the ratio of the Stokes depth scale to the local OBL depth, that is, the ratio Ds/|h|, along the tracks X = ±60 km. The OBL depth h is estimated using the maximum temperature gradient method described in Sullivan et al. (1998); essentially h is the spatial average of the boundary layer depths h(x, y) deduced along every vertical column in the LES grid box. A comparison of Fig. 18 with Fig. 19 shows that the largest impact of Langmuir turbulence on the entrainment flux occurs in the deepening phase of the OBL, that is, when the ratio Ds/|h| > 0.75; a similar result is reported by Skyllingstad et al. (2000). Based on the definition of Ds given by (7), this empirical ratio corresponds to a nondimensional e-folding length scale for a monochromatic Stokes wave of zL/|h| = (3/16π) or about 6% of the OBL depth. We note that in the simulations of McWilliams et al. (1997) and Sullivan et al. (2007) Ds/|h| ~ 1.82 and they find significant wave effects on entrainment. In the present simulations when Ds/|h| > 0.75, the coherent cells are of sufficient strength to effectively transport momentum and kinetic energy to the base of the thermocline and thereby increase the entrainment of cool water. In the hurricane wake, ttm > 0, the ratio Ds/|h| < 0.5 because the OBL mixes to very deep levels by inertial resonance and Ds is reduced by the falling surface winds. The ratio Ds/|h| appears to be an empirical indicator of when Langmuir turbulence impacts entrainment in a stratified OBL.

Fig. 19.
Fig. 19.

Temporal variation of the OBL depth |h| (left ordinate) and ratio Ds/|h| (right ordinate). The solid black (gray) lines are simulations with vortex force on the right (left) side of the storm. The Stokes depth scale Ds is predicted by (7). The dashed black (gray) lines are estimates of |h| on the right (left) side for simulations with no vortex force. For reference, the LES by McWilliams et al. (1997) is carried out with the ratio Ds/|h| ~ 1.82.

Citation: Journal of Physical Oceanography 42, 11; 10.1175/JPO-D-12-025.1

7. Summary and conclusions

A turbulence resolving large-eddy simulation (LES) model is used to study ocean boundary layers (OBLs) driven by extreme (hurricane) winds and waves. The LES model adopts wave-averaged equations that include the vortex force , where us is the wave field Stokes drift velocity and is the vorticity. Large-scale time-varying wind and wave fields input to the LES are generated by WaveWatch III; the wind and wave fields are an idealization of Hurricane Frances traveling at a speed of 5.5 m s−1. The surface wind and wave fields are in a nonequilibrium state and, as a result, the Stokes drift velocity is a complex two-component vector that varies with time, storm track position, and spirals with depth. High-resolution LES using as many as 1024 × 1024 × 256 grid points integrated for more than 70 physical hours are conducted at fixed positions right and left (X = ±60 km) of the storm track.

The surface-based turbulent Langmuir number Lat, which is used to assess the relative importance of shear and wave forcing, varies from 0.2 to more than 0.6 but over wide regions Lat ~ 0.3. Thus the hurricane forced OBL is often in a Langmuir turbulence regime. The LES solutions exhibit important differences from previous simulations that adopt wind wave equilibrium. The potency of Langmuir turbulence and the size and strength of the coherent Langmuir cells depends on the locally evolving wind wave state. At the time of maximum winds on the right-hand side of the storm, the near inertial rotation of the winds creates a “resonant” Stokes drift that is well aligned with the surface winds over a large fraction of the OBL. Imagery shows that the resulting instantaneous Langmuir cells can be more than 100 m in length with downwelling velocities exceeding 5u* (u* is the friction velocity) and are depth filling. The cells are much weaker and smaller on the left-hand side because of wind wave misalignment. The latter results in oppositely signed vertical gradients of current and Stokes drift. Hence the vertical momentum fluxes can be countergradient to ∂us/∂z, which results in less production of vertical velocity variance near the water surface. Vertical momentum fluxes countergradient to vertical gradients of Stokes drift are one of many possible regimes in nonequilibrium winds and waves. The average coherent Langmuir cell, identified using linear stochastic estimation, rotates, grows, and shrinks depending on Lat, and the alignment between the surface wind and depth-varying Stokes profile. Over time on both sides of the storm track, the cells tend to closely align with the mean Lagrangian shear direction.

The vertical profiles of mean currents and temperature, turbulence variances, and momentum and heat fluxes, are impacted by Langmuir turbulence to varying degrees. At the time of maximum winds, Langmuir turbulence enhances the minimum heat flux (or entrainment flux) at the base of the OBL by 20% and lowers the sea surface temperature by 0.25 K on the right-hand side of the storm. Inertially rotating wind stress, however, remains the dominant dynamical process for entraining cold water into the OBL. Wave effects also produce changes to entrainment and SST on the left-hand side of the storm, but to a lesser degree. In the nonequilibrium wind wave regime studied here the surface based turbulent Langmuir number Lat appears to be insufficient in describing when wave effects impact the OBL. The shape of the Stokes drift vertical profile, wind wave alignment, and the depth scale of the Stokes profile Ds are important bulk parameters. When the ratio Ds/|h| > 0.75, where h is the OBL depth, Langmuir turbulence enhances entrainment at the thermocline. Wave breaking, which is a dominant feature of the high wind OBL, is omitted here but will be considered in future work.

Acknowledgments

We thank the anonymous reviewers for their comments and suggestions that improved this work. PPS and JCM were supported by the Physical Oceanography Program through the Office of Naval Research (ONR) Departmental Research Initiatives, Impact of Typhoons on the Ocean in the Pacific (ITOP) and High Resolution Air-Sea Interaction (HiRes). PPS is also supported by the National Science Foundation through the National Center for Atmospheric Research. LR and WKM were supported by grants to WKM from NSF (Physical Oceanography) and ONR HiRes. This research used computer resources provided by NCAR and the Department of Defense High Performance Computing Modernization Program.

REFERENCES

  • Adrian, R. J., 1996: Stochastic estimation of the structure of turbulent fields. Eddy Structure Identification, J. P. Bonnet, Ed., Springer Verlag, 145–196.

  • Adrian, R. J., B. G. Jones, M. K. Chung, Y. Hassan, C. K. Nithianandan, and A. T.-C. Tung, 1989: Approximation of turbulent conditional averages by stochastic estimation. Phys. Fluids, 1A, 992998.

    • Search Google Scholar
    • Export Citation
  • Banner, M. L., I. Jones, and J. C. Trinder, 1989: Wavenumber spectra of short gravity waves. J. Fluid Mech., 25, 321344.

  • Belcher, S. E., and Coauthors, 2012: A global perspective on Langmuir turbulence in the ocean surface boundary layer. Geophys. Res. Lett., doi:10.1029/2012GL052932, in press.

    • Search Google Scholar
    • Export Citation
  • Bell, M. M., M. T. Montgomery, and K. A. Emanuel, 2012: Air–sea enthalpy and momentum exchange at major hurricane wind speeds observed during CBLAST. J. Atmos. Sci., 69, 31973222.

    • Search Google Scholar
    • Export Citation
  • Black, P., and Coauthors, 2007: Air–sea exchange in hurricanes: Synthesis of observations from the Coupled Boundary Layers Air–Sea Transfer experiment. Bull. Amer. Meteor. Soc., 88, 357374.

    • Search Google Scholar
    • Export Citation
  • Capet, X., J. C. McWilliams, M. J. Molemaker, and A. Shchepetkin, 2008: Mesoscale to submesoscale transition in the California Current system. Part I: Flow structure, eddy flux, and observational tests. J. Phys. Oceanogr., 38, 2943.

    • Search Google Scholar
    • Export Citation
  • Chini, G. P., K. Julien, and E. Knobloch, 2009: An asympotically reduced model of turbulent Langmuir circulation. Geophys. Astrophys. Fluid Dyn., 103, 179197.

    • Search Google Scholar
    • Export Citation
  • Christensen, K. T., and R. J. Adrian, 2001: Statistical evidence of hairpin vortex packets in wall turbulence. J. Fluid Mech., 431, 433443.

    • Search Google Scholar
    • Export Citation
  • Craik, A., and S. Leibovich, 1976: A rational model for Langmuir circulations. J. Fluid Mech., 73, 401426.

  • Crawford, G. B., and W. G. Large, 1996: A numerical investigation of resonant inertial response of the ocean to wind forcing. J. Phys. Oceanogr., 26, 873891.

    • Search Google Scholar
    • Export Citation
  • D’Asaro, E., 2001: Turbulent vertical kinetic energy in the ocean mixed layer. J. Phys. Oceanogr., 31, 35303537.

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

  • D’Asaro, E., T. B. Sanford, P. P. Niiler, and E. J. Terrill, 2007: Cold wake of Hurricane Frances. Geophys. Res. Lett., 34, L13604, doi:10.1029/2007GL030160.

    • Search Google Scholar
    • Export Citation
  • Deardorff, J. W., 1972: Numerical investigation of neutral and unstable planetary boundary layers. J. Atmos. Sci., 29, 91115.

  • 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
  • Fan, Y., I. Ginis, T. Hara, C. W. Wright, and E. J. Walsh, 2009: Numerical simulations and observations of surface wave fields under an extreme tropical cyclone. J. Phys. Oceanogr., 39, 20972116.

    • Search Google Scholar
    • Export Citation
  • French, J. R., W. M. Drennan, J. A. Zhang, and P. G. Black, 2007: Turbulent fluxes in the hurricane boundary layer. Part I: Momentum flux. J. Atmos. Sci., 64, 10891102.

    • Search Google Scholar
    • Export Citation
  • Gnanadesikan, A., and R. A. Weller, 1995: Structure and instability of the Ekman spiral in the presence of surface gravity waves. J. Phys. Oceanogr., 25, 31483171.

    • Search Google Scholar
    • Export Citation
  • Grant, A. L. M., and S. E. Belcher, 2009: Characteristics of Langmuir turbulence in the ocean mixed layer. J. Phys. Oceanogr., 39, 18711887.

    • Search Google Scholar
    • Export Citation
  • Grant, A. L. M., and S. E. Belcher, 2011: Wind-driven mixing below the oceanic mixed layer. J. Phys. Oceanogr., 41, 15561575.

  • Harcourt, R. R., and E. A. D’Asaro, 2008: Large-eddy simulation of Langmuir turbulence in pure wind seas. J. Phys. Oceanogr., 38, 15421562.

    • Search Google Scholar
    • Export Citation
  • Holm, D. D., 1996: The ideal Craik-Leibovich equations. Physica D, 98, 415441.

  • Huang, N. E., 1979: On surface drift currents in the ocean. J. Fluid Mech., 91, 119208.

  • Jacob, S. D., and L. K. Shay, 2003: The role of oceanic mesoscale features on the tropical cyclone-induced mixed layer response: A case study. J. Phys. Oceanogr., 33, 649676.

    • Search Google Scholar
    • Export Citation
  • Jacob, S. D., L. K. Shay, A. J. Mariano, and P. G. Black, 2000: The 3D oceanic mixed layer response to Hurricane Gilbert. J. Phys. Oceanogr., 30, 14071429.

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

  • Kukulka, T., A. J. Plueddemann, J. H. Trowbridge, and P. P. Sullivan, 2010: Rapid mixed layer deepening by the combination of Langmuir and shear instabilities: A case study. J. Phys. Oceanogr., 40, 23812400.

    • Search Google Scholar
    • Export Citation
  • Large, W. G., and S. Pond, 1981: Open ocean 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., 1983: The form and dynamics of Langmuir circulations. Annu. Rev. Fluid Mech., 15, 391427.

  • Li, M., and C. Garrett, 1993: Cell merging and the jet/downwelling ratio in Langmuir circulations. 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 J. M. Restrepo, 1999: The wave-driven ocean circulation. J. Phys. Oceanogr., 29, 25232540.

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

  • McWilliams, J. C., J. R. Restrepo, and E. M. Lane, 2004: An asymptotic theory for the interaction of waves and currents in shallow coastal water. J. Fluid Mech., 511, 135178.

    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., E. Huckle, J.-H. Liang, and P. P. Sullivan, 2012: The wavy Ekman layer: Langmuir circulations, breaking waves, and Reynolds stress. J. Phys. Oceanogr., 42, 17931816.

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

  • Melville, W. K., L. Romero, and J. M. Kleiss, 2005: Extreme wave events in the Gulf of Tehuantepec. Rogue Waves: Proc. 14th ‘Aha Huliko’a Hawaiian Winter Workshop, Honolulu, HI, University of Hawaii at Manoa, 23–28.

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

    • Search Google Scholar
    • Export Citation
  • Moeng, C.-H., and P. P. Sullivan, 1994: A comparison of shear and buoyancy driven planetary-boundary-layer flows. J. Atmos. Sci., 51, 9991022.

    • 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
  • Polton, J. A., D. M. Lewis, and S. E. Belcher, 2005: The role of wave-induced Coriolis–Stokes forcing on the wind-driven mixed layer. J. Phys. Oceanogr., 35, 444457.

    • Search Google Scholar
    • Export Citation
  • Powell, M. D., P. J. Vickery, and T. A. Reinhold, 2003: Reduced drag coefficient for high wind speeds in tropical cyclones. Nature, 422, 279283.

    • Search Google Scholar
    • Export Citation
  • Price, J. F., 1981: Upper-ocean response to a hurricane. J. Phys. Oceanogr., 11, 153175.

  • Price, J. F., R. A. Weller, and R. Pinkel, 1986: Diurnal cycling: Observations and models of the upper ocean response to diurnal heating, cooling, and wind mixing. J. Geophys. Res., 91, 84118427.

    • Search Google Scholar
    • Export Citation
  • Romero, L., and W. K. Melville, 2010a: Airborne observations of fetch-limited waves in the Gulf of Tehuantepec. J. Phys. Oceanogr., 40, 441465.

    • Search Google Scholar
    • Export Citation
  • Romero, L., and W. K. Melville, 2010b: Numerical modeling of fetch-limited waves in the Gulf of Tehuantepec. J. Phys. Oceanogr., 40, 466486.

    • Search Google Scholar
    • Export Citation
  • Sanford, T. B., J. F. Price, J. B. Girton, and D. C. Webb, 2007: Highly resolved observations and simulations of the ocean response to a hurricane. Geophys. Res. Lett., 34, L13604, doi:10.1029/2007GL029679.

    • Search Google Scholar
    • Export Citation
  • Sanford, T. B., J. F. Price, and J. B. Girton, 2011: Upper-ocean response to Hurricane Frances (2004) observed by profiling EM-Apex floats. J. Phys. Oceanogr., 41, 10411056.

    • Search Google Scholar
    • Export Citation
  • Schmidt, H., and U. Schumann, 1989: Coherent structure of the convective boundary layer. J. Fluid Mech., 200, 511562.

  • 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. D., W. D. Smyth, J. N. Moum, and H. Wijesekera, 1999: Upper-ocean turbulence during a westerly wind burst: A comparison of large-eddy simulation results and microstructure measurements. J. Phys. Oceanogr., 29, 528.

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

    • Search Google Scholar
    • Export Citation
  • Smith, J. A., 1998: Evolution of Langmuir circulation during a storm. J. Geophys. Res., 103 (C6), 12 64912 668.

  • Snyder, R. L., F. W. Dobson, J. A. Elliott, and R. B. Long, 1981: Array measurements of atmospheric pressure fluctuations above surface gravity waves. J. Fluid Mech., 102, 159.

    • Search Google Scholar
    • Export Citation
  • Sullivan, P. P., and J. C. McWilliams, 2010: Dynamics of winds and currents coupled to surface waves. Annu. Rev. Fluid Mech., 42, 1942.

    • Search Google Scholar
    • Export Citation
  • Sullivan, P. P., and E. G. Patton, 2011: The effect of mesh resolution on convective boundary-layer statistics and structures generated by large-eddy simulation. J. Atmos. Sci., 68, 23952415.

    • Search Google Scholar
    • Export Citation
  • Sullivan, P. P., J. C. McWilliams, and C.-H. Moeng, 1994: A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows. Bound.-Layer Meteor., 71, 247276.

    • Search Google Scholar
    • Export Citation
  • Sullivan, P. P., C.-H. Moeng, B. Stevens, D. H. Lenschow, and S. D. Mayor, 1998: Structure of the entrainment zone capping the convective atmospheric boundary layer. J. Atmos. Sci., 55, 30423064.

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

    • Search Google Scholar
    • Export Citation
  • Sullivan, P. P., J. C. McWilliams, and W. K. Melville, 2007: Surface gravity wave effects in the oceanic boundary layer: Large-eddy simulation with vortex force and stochastic breakers. J. Fluid Mech., 593, 405452.

    • Search Google Scholar
    • Export Citation
  • Teixeira, M. A. C., and S. E. Belcher, 2002: On the distortion of turbulence by a progressive surface wave. J. Fluid Mech., 458, 229267.

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

    • Search Google Scholar
    • Export Citation
  • Tolman, H. L., 2002: User manual and system documentation of Wavewatch-III, version 2.22. NOAA/NWS/NCEP/MMAB Tech. Rep. Tech. Note 222, 133 pp.

  • Van Roekle, L. P., B. Fox-Kemper, P. P. Sullivan, P. E. Hamlington, and S. R. Haney, 2012: The form and orientation of Langmuir cells for misaligned winds and waves. J. Geophys. Res., 117, C05001, doi:10.1029/2011JC007516.

    • Search Google Scholar
    • Export Citation
  • van Vledder, G. P., 2006: The WRT method for the computation of non-linear four-wave interactions in discrete spectral wave models. Coastal Eng., 53, 223242.

    • Search Google Scholar
    • Export Citation
  • Wright, C., and Coauthors, 2001: Hurricane directional wave spectrum spatial variation in the open ocean. J. Phys. Oceanogr., 31, 24722488.

    • Search Google Scholar
    • Export Citation
  • Young, I. R., 1998: Observations of the spectra of hurricane generated waves. Ocean Eng., 25 (4–5), 261276.

  • Young, I. R., 2003: A review of the sea state generated by hurricanes. Mar. Struct., 16, 201218.

  • Young, I. R., 2006: Directional spectra of hurricane wind-waves. J. Geophys. Res., 111, C08020, doi:10.1029/2006JC003540.

  • Zedler, S. E., 2007: Strong wind forcing of the ocean. Ph.D. thesis, University of California, San Diego, 127 pp.

  • Zedler, S. E., T. D. Dickey, S. C. Doney, J. F. Price, X. Xu, and G. L. Mellor, 2002: Analyses and simulations of the upper ocean’s response to Hurricane Felix at the Bermuda Testbed Mooring site: 13–23 August 1995. J. Geophys. Res., 107, 129.

    • Search Google Scholar
    • Export Citation
  • Zedler, S. E., P. P. Niiler, D. Stammer, E. Terrill, and J. Morzel, 2009: Ocean’s response to Hurricane Frances and its implications for drag coefficient parameterizations at high wind speeds. J. Geophys. Res., 114, C04016, doi:10.1029/2008JC005205.

    • Search Google Scholar
    • Export Citation
Save
  • Adrian, R. J., 1996: Stochastic estimation of the structure of turbulent fields. Eddy Structure Identification, J. P. Bonnet, Ed., Springer Verlag, 145–196.

  • Adrian, R. J., B. G. Jones, M. K. Chung, Y. Hassan, C. K. Nithianandan, and A. T.-C. Tung, 1989: Approximation of turbulent conditional averages by stochastic estimation. Phys. Fluids, 1A, 992998.

    • Search Google Scholar
    • Export Citation
  • Banner, M. L., I. Jones, and J. C.