The development of tropical cyclones (TCs) strongly depends on the air-sea interactions that include heat fluxes and momentum transfer (Emanuel 1991(Emanuel , 1999)). The TC's strong wind and associated wave forcing drives upper-ocean currents that generate vigorous turbulence.
Turbulent eddies erode the thermocline by entraining deep cool water into the warmer upper layer, resulting in ocean surface boundary layer (OSBL) deepening and sea surface cooling (Price 1981). In turn, sea surface cooling reduces air-sea heat fluxes that drive the TC, resulting in a negative feedback between TC winds and sea surface temperature (Bender and Ginis 2000;Ginis 2002). The inertial resonance between the turning wind stress and surface currents is also a critical dynamical process under TCs because it increases the shear at the mixed layer base, leading to stronger mixing on the right-hand side of TCs (Price 1981;Skyllingstad et al. 2000;Sanford et al. 2011;Sullivan et al. 2012;Reichl et al. 2016b).
Recent studies indicate that wave-driven Langmuir turbulence (LT) plays an important role in upper-ocean Denotes content that is immediately available upon publication as open access.
a Current affiliation: Program in Atmospheric and Oceanic Sciences, Princeton University, Princeton, New Jersey. turbulence under TCs (Sullivan et al. 2012;Rabe et al. 2015;Reichl et al. 2016b,a). LT was originally observed in moderate wind conditions as parallel bands of floating material on the sea surface, which are due to strong surface current convergences of horizontal roll vortices in the OSBL (Langmuir 1938). Over the last decades, comprehensive field observations, mostly conducted in moderate wind conditions, have revealed characteristic features of LT, such as strong surface convergence regions, downwelling jets, and the spacing of roll vortices from several meters to kilometers (Thorpe 2004;Weller and Price 1988;Farmer and Li 1995;Plueddemann et al. 1996;Smith 1992;Gargett et al. 2004;Gargett and Grosch 2014). A systematic mathematical theory of LT is based on the wave-averaged Navier-Stokes equation, the so-called Craik-Leibovich (CL) equation, and suggests that LT is driven by the CL vortex force, which is the cross-product of Stokes drift and vorticity vectors (Craik and Leibovich 1976). Physically, the Stokes drift shear tilts vertical vorticity into the direction of wave propagation, generating LT. Today LT is recognized as a fundamental upper-ocean turbulent process (McWilliams et al. 1997;Thorpe 2004;Li et al. 2005;Sullivan and McWilliams 2010;Belcher et al. 2012;D'Asaro 2014) that contributes significantly to turbulent transport and windand wave-driven mixed layer deepening (Kukulka et al. 2009(Kukulka et al. , 2010;;Grant and Belcher 2011).
Previous investigations of LT involve turbulenceresolving large-eddy simulation (LES) that is based on the filtered CL equations with explicit wave effects, which resolves LT and associated relatively large vortical structures (Skyllingstad and Denbo 1995;McWilliams et al. 1997). LES studies show that LT enhances the vertical fluxes of momentum and heat, inducing stronger vertical velocity variance and mixed layer deepening (McWilliams et al. 1997;Li et al. 2005;Polton and Belcher 2007;Grant and Belcher 2009;Kukulka et al. 2009;Noh et al. 2009). Direct comparisons of LES results with ocean observations reveal that LES captures in detail many of the observed LT characteristics (Skyllingstad et al. 1999;Gargett et al. 2004;Li et al. 2009;Kukulka et al. 2009Kukulka et al. , 2013;;D'Asaro et al. 2014). However, most LES studies are conducted in moderate wind conditions with monochromatic waves, and only a few of them examine LT in extreme TC conditions. Sullivan et al. (2012) explored LT dynamics under a TC by forcing an LES model with realistic TC winds and waves, which were simulated by a spectral wave model. They contrasted time series of OSBL turbulence statistics at two stations: one on the right-hand side (rhs) with strong inertial resonance and the other one on the lefthand side (lhs) with weak inertial resonance. The intensity of LT strongly depends on location and time because of the TC's complex wind and wave forcing, yielding more energetic LT on the rhs of the TC. Furthermore, their results indicate that the direction of roll vortices due to LT is aligned with the wind direction and tracks the Lagrangian shear direction. As a result of the TC's transient forcing conditions, wind vector and wave propagation directions are misaligned, reducing LT intensity.
Motivated by OSBL observations of depth-averaged vertical velocity variance (VVV) obtained from Lagrangian floats under Hurricane Gustav (2008), Rabe et al. (2015) investigated OSBL turbulence with LES experiments forced by Gustav (2008)'s winds and waves. Simulated VVV is only consistent with the observed VVV with LT, that is, for simulations with CL vortex force, indicating LT's significant role in OSBL dynamics. LES results demonstrate that LT enhances VVV and varies with complex sea states found under TCs. Misaligned wind and wave fields near the TC eye are associated with an observed suppression of VVV, which is also predicted by the LES. Thus, wind-wave misalignment can reduce VVV and suppress LT to the levels close to shear turbulence (ST).
Building on this previous work (Sullivan et al. 2012;Rabe et al. 2015), we recently designed a series of LES experiments for the full spatial TC extent to develop a turbulence closure scheme with explicit sea-statedependent LT effects (Reichl et al. 2016b) and to investigate the role of sea-state-dependent LT in the OSBL response to TCs (Reichl et al. 2016a). In regional ocean models under TCs, which are commonly based on the Reynolds-averaged Navier-Stokes (RANS) equations, LT cannot be resolved, so that smaller-scale turbulent transport processes in the OSBL have to be parameterized. We modified the K-profile parameterization (KPP) model (Large et al. 1994) to match mean current and temperature profiles obtained from the LES model. In the KPP model, we replaced the Eulerian current with the Lagrangian current (Eulerian current plus Stokes drift) to compute the turbulent momentum flux following McWilliams et al. (2012) and also introduced turbulence enhancement factors following McWilliams and Sullivan (2000). Our new KPP model with explicit sea-state-dependent LT significantly improves estimations of LES temperature and currents compared to results of the standard (unmodified) KPP model (Reichl et al. 2016b). We next introduced this new KPP model in a regional three-dimensional coupled wave-ocean RANS model to demonstrate that sea-state-dependent LT substantially modifies the threedimensional OSBL response (Reichl et al. 2016a). Results indicate that LT reduces upwelling and horizontal advection as a result of enhanced near-surface mixing and that simulations without sea-state-dependent LT cannot accurately reproduce the sea surface cooling and horizontal transport. More recently, Blair et al. (2017) investigated the upper-ocean response under Hurricane Edouard (2014) with our new KPP model and satellite observations, suggesting the importance of sea-statedependent LT on the mixed layer depth evolution.
In this study, we use the same LES approach as in Reichl et al. (2016b) to comprehensively investigate LT dynamics, OSBL energetics, and the influence of inertial currents on the OSBL evolution for a wide range of realistic TC conditions. We first review our basic numerical approach and provide an overview of complex wind and wave conditions under TCs (section 2) and then illustrate that the LT-driven OSBL is not only seastate dependent but also influenced by inertial currents in TC conditions (section 3).
To analyze LT in TC conditions, we use the same wind, wave, and turbulence modeling approaches and datasets as in our previous study (Reichl et al. 2016b), which are briefly reviewed in the following subsections.
The TC wind field is constructed based on a Holland wind model (Holland 1980(Holland , 2008) ) with the radius of maximum wind (RMW) of 50 km, maximum wind speed at 10-m height of 65 m s 21 , and a translation speed of 5 m s 21 , which represent typical TC parameters (Reichl et al. 2016b). The output of the wind model is the wind velocity at 10-m height with speed U 10 over a domain with a length from 2648 to 648 km in the TC's propagation direction (along X) and a width from 2500 to 500 km across the TC's propagation direction (along Y) (Fig. 1, top left). In addition to the horizontally averaged turbulent statistics investigated by Reichl et al. (2016b), in this study we also output and investigate the fourdimensional spatiotemporal high-resolution LES temperature, pressure, and velocity fields.
The wind stress vector at the sea surface has the same direction as the wind velocity, and its magnitude is parameterized by t 5 r a C d U 2 10 , where r a is the air density and C d is the drag coefficient that depends on U 10 as follows (Sullivan et al. 2012 (1)
Recent studies show the complexity of drag coefficient under TCs, which varies with different storm quadrants and even wind directions (Holthuijsen et al. 2012;Hsu et al. 2017). Since there is substantial uncertainty of the drag coefficient, we adopt Eq. ( 1) with simpler assumptions than what has been done in previous studies (Sullivan et al. 2012;Rabe et al. 2015;Reichl et al. 2016b).
Following previous approaches (McWilliams et al. 1997;Skyllingstad and Denbo 1995;Kukulka et al. 2010), we use an LES model to simulate the upper-ocean response to TC's wind and wave fields. The LES model solves the grid-filtered CL equations (Craik and Leibovich 1976). LT is generated by the Craik-Leibovich vortex force, which is the cross-product of the Stokes drift u s and vorticity. When the wave effect is not considered (u s 5 0), the LES model simply solves gridfiltered Navier-Stokes equations that only generate ST. The Stokes drift is obtained through the wave model, introduced in the following section (section 2c). We only solve the density equation in the LES model and assume that temperature is linearly related to density by
where r is the water density, r 0 5 1024 kg m 23 is the reference density, T is water temperature, and g 5 2 3 10 24 K 21 is the thermal expansion coefficient. LES experiments are performed for 18 stations across the TC translation direction along Y from 2200 to 200 km, with a minimum and maximum spacing between stations of 20 and 50 km, respectively (Table 1). The coordinate Y is zero at the TC center, and Y . 0 and Y , 0 are on the rhs and lhs of the TC eye, respectively. The LES domain size spans horizontally 750 m 3 750 m and is 240 m deep with 256 grid points in each direction for experiments in Table 1. Inside the LES domain we use x i with index i 5 1, 2, 3 to denote the spatial coordinate (x, y, z) in east, north, and vertical directions, respectively. For all experiments the initial potential temperature profile features a 10-m homogeneous mixed layer with T 5 302:4 K and a constant stratification with dT/dz 5 0:04 K m 21 below the mixed layer, where z denotes the depth. Small constant surface cooling of 25 W m 22 is imposed (McWilliams et al. 1997) for all simulations, which facilitates the initial spinup of turbulence but is otherwise insignificant for the OSBL dynamics presented here. The LES model is forced at the surface with the modeled wind stress vector (section 2a), and the Stokes drift vector u s is imposed based on simulated 2D wave height spectra (section 2c).
The third-generation wave model WAVEWATCH III (Tolman 2009) is used to simulate the directional frequency spectra of surface gravity waves in TC conditions following Reichl et al. (2016b). Its computational domain is 3000 km long in the TC's translation direction and 1800 km wide across the TC's translation direction. The wave model has a horizontal grid spacing of 8.33 km, and the wave spectrum is discretized into 48 evenly spaced directions and 40 logarithmically spaced frequencies. The simulated wave field is stationary in a coordinate system translating with the TC. The skill of WAVEWATCH III physics parameterizations to simulate complex features of the wave field on the left of the translation direction remains a topic of research (see The dashed magenta lines denote two LES transects which will be examined in the main text. The two black circles on one transect, marked by capital letters A and B, are two locations that will be discussed in detail. Note that the x axis represents the translating distance of the TC, and the y axis indicates the distance of each station from the TC's path. Both distances are normalized by the RMW. The TC travels from right to left, and (X 2 X 0 ) 5 0 is where the TC's eye passes. Hsu et al. 2018). However, comparison of the version of WAVEWATCH III employed for this study with wave observations shows good skill on average to predict the sea state under extreme hurricane conditions (Fan et al. 2009).
The most energetic wave fields with large H s and long l m are found on the rhs of the TC near the RMW at (Y 2 Y 0 )/RMW 5 1 (Fig. 1, top center and top right). This is because waves travel with the TC and are thus exposed to greater wind forcing, creating favorable conditions for wave development. Note that the TC translates from right [(X 2 X 0 ) . 0] to left [(X 2 X 0 ) , 0], and the X axis is equivalent to time t 5 (X 2 X 0 )/(5 m s -1 ). In contrast, waves on the lhs are much weaker with smaller l m and H s because the wind direction opposes the TC direction. Longer waves far ahead of the TC near (X 2 X 0 )/RMW 5 -2 indicate the presence of fasttraveling swell waves, but their H s value is small so these waves do not significantly contribute to wave forcing (Fig. 1, top right). As waves are more developed on the rhs of the TC, the misalignment between the wind direction u w and the peak wave spectrum u P is smaller on the rhs (Fig. 1, bottom left). In other regions under the TC, however, the propagation direction of energetic waves differs significantly from the wind direction, in particular on the lhs.
The time-and space-dependent Stokes drift in the LES model is obtained by integrating the simulated wave spectra (Kenyon 1969):
where F(v, u) is the two-dimensional wave height spectrum, k is the horizontal wavenumber vector, k 5 jkj denotes its magnitude, v is the radian frequency, and u is the wave propagation direction. It is assumed that surface gravity waves satisfy the deep-water dispersion relation v 5 ffiffiffiffiffi ffi gk p , where g is the acceleration of gravity. Note that the largest parameterized wavenumber in the model is 400 m 21 , which specifies the upper bound of radian frequency in Eq. (3) (Reichl et al. 2016b).
Strong wind-wave misalignment in TC conditions reduces LT's intensity (Sullivan et al. 2012;Rabe et al. 2015;Reichl et al. 2016b). To scale LT for conditions with windwave misalignment, Van Roekel et al. (2012) modified the surface layer Langmuir number from Harcourt and D'Asaro (2008) that scales LT forced by equilibrium wind-wave spectra. The so-called projected surface layer Langmuir number projects the wind and wave forcing onto the Langmuir cell direction and is given by
where u w is the wind direction, u s is the direction of depth-averaged Stokes drift within 0:2H B , and a is the direction of the depth-averaged Lagrangian shear, which is an estimate for the direction of Langmuir cells and defined by tan(a) 5
The Eulerian current and the Stokes drift in the (x, y) direction are denoted as (u, y) and (u s , y s ), respectively, and hi denotes a horizontal spatial average inside the LES numerical domain. The spatial distribution of La SLu with relatively small La SLu values (Fig. 1, bottom center) indicates that LT plays an important role in OSBL dynamics and that LT is stronger on the rhs (Reichl et al. 2016b). Interestingly, small values of La SLu do not always coincide with large values of H s and l m , suggesting that LT is not necessarily more pronounced for more energetic or more developed wave fields.
Recent studies show that the Stokes drift decay length scale influences the dynamics and structure of LT (Sullivan et al. 2012;Gargett and Grosch 2014;Kukulka and Harcourt 2017). To account for the effects of Stokes drift shear for wave spectra, we compute the length scale d s for the depth-averaged Stokes drift shear as
In TC conditions, we find that d s is closely related to the penetration depth scale of the Stokes drift introduced by Sullivan et al. ( 2012), for which the misalignment between surface Stokes drift and Stokes drift at greater depth is considered. The differences between d s and l m in TC conditions further illustrate that l m is a poor estimator of Stokes drift decay length scale that drives LT (Fig. 1, top right and bottom right). In summary, traditional parameters that describe the sea state are not sufficient to characterize LT dynamics in TC conditions.
We first investigate the response of the OSBL to tropical cyclones and identify regions of relatively large entrainment of deep cold water into the OSBL for the LT and ST cases (section 3a). Then we investigate different mechanisms that induce greater entrainment in the LT case (section 3b) and the ST case (section 3c).
a. The response of OSBL to tropical cyclones
The cold wake under the TC is caused by upper-ocean mixing accompanied by mixed layer deepening (Price 1981;Sullivan et al. 2012). The mixed layer depth H M is defined as the depth of largest temperature gradient where the stratification is strongest. The spatial and temporal distribution of H M shows an unsymmetrical pattern with greater H M on the rhs of the TC due to resonant wind forcing that drives strong inertial currents (Price 1981;Skyllingstad et al. 2000;Sanford et al. 2011) (Fig. 2, top left). LT enhances mixed layer deepening, and the contribution of LT is spatially variant (Fig. 2, top right). Enhanced mixed layer deepening due to LT is most pronounced near the RMW, especially on the rhs of the TC where strong, wind-aligned waves create favorable conditions for LT (cf. the forcing shown in the Fig. 1). However, after the storm passes, H M in both LT and ST cases approaches a similar final value on the rhs of the TC (Fig. 2, top right). This suggests greater entrainment of cool water into the OSBL for the ST case when the wind subsides, which compensates the smaller mixed layer deepening near the RMW, finally resulting in a similar H M .
Because of the highly transient forcing conditions under TCs, it is useful to introduce an OSBL depth characterized by active turbulent mixing. Following previous approaches for atmospheric boundary layer turbulence (Zilitinkevich et al. 2007), we specify a turbulent boundary layer depth H B where the turbulent stress decays to 5% of its surface value. Figure 3 shows the normalized turbulent stress profiles at two example transects with two specific locations (transects are denoted by a dashed magenta line and locations are marked with A and B in the bottom left of Fig. 2). Location A is on the rhs of the TC near the RMW where LT rapidly changes H M . Location B is also on the rhs but about 2.5 times the RMW behind the TC's eye in a region where H M is similar for the LT and ST case.
As expected, H B is close to H M during rapid mixed layer deepening and strong wind and wave forcing (regions outside the solid magenta contour in the bottom right of Fig. 2). However, under the eye, where winds and turbulence levels are relatively weak, H M does not accurately characterize the boundary layer depth because H B is much smaller than H M (regions with positive values in the bottom right of Fig. 2). Across the TC eye, the winds decrease for 21 , (X 2 X 0 )/RMW , 0, and the normalized turbulent stress is relatively large owing to relatively strong residual turbulence (Fig. 3, bottom). After the eye passes for 0 , (X 2 X 0 )/RMW , 1, the normalized turbulent stress is relatively small as a result of relatively weak residual turbulence. These history effects are more pronounced for the ST case, which indicates that the turbulence in the LT case is more responsive to the transient TC forcing.
Behind the TC eye, we observe active turbulence (jt t j/u 2 * . 0:2) at great depth (z , 20:5H B ), although winds are weak (Fig. 3, top). The turbulent stress profile in the ST case shows a more homogenized structure in the upper OSBL (z . 20:5H B ) (Fig. 3, top right), while the turbulent stress in the LT case decays faster with depth after the TC passes (Fig. 3, top left). This difference in stress profiles is related to more energetic inertial currents in the ST case (discussed in section 3c).
To further examine the effect of TC's transient forcing on the OSBL depth over the extent of the whole TC domain the rhs with (X 2 X 0 )/RMW . 2 enclosed by the solid blue contour in the bottom left of Fig. 2]. For (X 2 X 0 )/RMW . 1:8, H B exceeds H E , suggesting the generation of turbulence that is due to inertial currents instead of local winds (Fig. 3). Such relatively large H B . H E and the persistent mixed layer deepening behind the TC eye suggest the importance of inertial currents in influencing OSBL dynamics.
To understand how turbulence drives OSBL deepening, we first examine the evolution of profiles of the resolved turbulent buoyancy fluxes 2hw 0 Bi at two transects (Fig. 4). Here B 5 r 0 g/r 0 is the buoyancy, w is the vertical velocity, and the prime denotes the deviation from the horizontal average denoted by symbol h i.
For both the LT and ST cases, the greatest buoyant fluxes occur below 0:6H M , which is in the region of the shear layer. The shear layer is defined as the layer with pronounced stratification and mean current shear and approximately locates below 0:6H M (temperature and current profiles will be shown in sections 3b and 3c) (Skyllingstad et al. 2000;Grant and Belcher 2011).
For the transect that is tangential to the RMW, the buoyancy fluxes in the LT case are larger than in the ST case near the maximum wind period (location A), suggesting that LT enhances buoyancy fluxes that contribute to enhanced mixed layer deepening (Fig. 4, top left).
However, buoyancy fluxes in the ST case increase behind the TC eye and are about one order larger than in the LT case at (X 2 X 0 )/RMW 5 2.5 (location B) (Fig. 4, top right). This is consistent with the relatively strong mixed layer deepening in the ST case behind the TC eye, which is due to the turbulence induced by strong inertial currents.
For the transect crossing the TC's eye region, we observe two periods with strong buoyancy fluxes for both the LT and ST case, which corresponds to the timing of two maximum winds as the TC eye passes (Fig. 4, bottom). The peaks of the buoyancy fluxes in the LT case are greater than in the ST case, indicating LT's significant role in enhancing the turbulent mixing even under transient winds. Greater buoyancy fluxes in the ST case are not found until (X 2 X 0 )/RMW . 2, where H B is deeper than H E , suggesting the presence of shear-driven turbulence due to energetic inertial currents (Fig. 2, bottom left).
To identify the influence of LT and inertial currents on mixed layer deepening over the full spatial extent of the TC, we examine the depth-integrated buoyancy fluxes, which drive changes in total (depth integrated) potential energy as
where SGS symbolizes subgridscale density fluxes (second term on the rhs), and H is the depth of LES domain. The greatest depth-integrated buoyancy flux occurs around the maximum wind radius where the maximum mixed layer deepening rates are also found (Fig. 5, left). The stronger turbulent buoyancy transport on the rhs agrees with the larger mixed layer deepening on the rhs (Fig. 2, top left). Pronounced differences of Ð 0 2H 2hw 0 Bi dz between the LT and the ST cases show that the enhanced entrainment by LT mainly occurs near the RMW (Fig. 5, right). In the ST case, entrainment is greater inside the TC eye as the result of enhanced leftover turbulence under transient winds (Fig. 3, bottom right). The greater entrainment on the rhs behind the TC eye is caused by the stronger turbulence generated by inertial currents (section 3c).
To explore the relation of wind forcing to entrainment, we scale the depth-integrated buoyancy fluxes Ð 0 2H 2hw 0 Bi dz by u 3 * , which is similar to the scaling proposed by Grant and Belcher (2009) (Fig. 6). In both the LT and the ST cases, significant deviations from this scaling are observed. Smaller values of Ð 0 2H 2hw 0 Bi dz are mainly found on the lhs of the storm, while greater values of Ð 0 2H 2hw 0 Bi dz are located mostly on its rhs (Fig. 6, left and center). In the ST case differences are caused by energetic inertial currents on the rhs, inducing shear-driven turbulent mixing. This is consistent with the stronger buoyancy fluxes and more homogenized turbulence stress profiles on the rhs, behind the TC eye.
In the LT case, differences are due to LT that is more vigorous on the rhs due to stronger wave forcing and weaker on the lhs due to wind-wave misalignment (refer to Fig. 1). To examine the influence of LT on entrainment under the TC's complicated wind and wave forcing, we consider a similar scaling approach as Grant and Belcher (2009) but replace La t with La SLu , which is more adequate in complex wind and wave conditions, so that Ð 0 2H 2hw 0 Bi dz is expected to scale with La 22 SLu u 3 * . This scaling significantly reduces the scatter, illustrating the importance of sea-state-dependent Langmuir turbulence in buoyancy entrainment (Fig. 6, right).
For aligned wind-wave conditions and relatively shallow mixed layers, LT enhances turbulent entrainment through a three-step process: first the large coherent structures in LT facilitate the transport of momentum through the boundary layer; then the shear-instability near the mixed layer base is locally enhanced, inducing greater erosion of thermocline; finally, the eroded colder water is transported upward and mixed by LT near the surface (Kukulka et al. 2010). In this section, we evaluate this process step by step in tropical cyclone conditions by examining the flow field at a location where enhanced mixed layer deepening and greater buoyancy fluxes due to LT are found (marked with A in the bottom of Fig. 2).
To address the enhanced shear instability by LT, we first examine the horizontally averaged Lagrangian current and temperature profiles at a location A (Fig. 7, top). The Lagrangian current is the sum of the Eulerian current and Stokes drift. To better illustrate the relation between wind and current directions, we project the horizontal Eulerian currents and Stokes drift into the along-wind and crosswind direction denoted by subscripts // and ?, respectively (Fig. 7, top left).
At location A, both the LT and the ST cases show predominant along-wind Eulerian currents compared to the crosswind currents. This is because the onset of the TC's maximum wind mainly accelerates currents and the inertial currents have not fully spun up yet. With LT, the vertical profile of hui // is more uniform in the upper OSBL with z . 20:5H B compared to the ST case (Fig. 7, top left). Large coherent structures in the presence of LT efficiently transport momentum downward, which homogenizes the upper OSBL and decreases the shear in the mean Eulerian currents. However, the near-surface Lagrangian current shear is more alike in both cases, illustrating the predominant Stokes drift shear in the LT case. In addition, the currents in the LT and the ST cases are also similar in the shear layer (z , 20:6H B ), although the buoyancy entrainment at the mixed layer base substantially differs in both cases.
To assess the influence of stratification on entrainment, we also investigate the temperature profiles at location A. The OSBL temperature with LT is about 0.5 K lower than without LT because of the enhanced entrainment (Fig. 7, top center). Within the shear layer, the temperature gradient in the LT case is slightly greater than the ST case, implying stronger stratification. To better understand the competition between destabilizing current shear and stabilizing stratification, we examine the gradient Richardson number:
In the ST case, Ri gradually increases with depth and approaches a value close to the critical value (Ri 5 0:25) near the mixed layer base (solid red line in the top right of Fig. 7). In the LT case, Ri significantly exceeds 0.25 in the upper half of the OSBL because currents are homogenized, illustrating the importance of LT in generating turbulence (solid black line in the top right of Fig. 7). In the lower half of the OSBL, Ri approaches a local minimum near H B , indicating the importance of locally generated shear instability that effectively erodes the thermocline and enhances entrainment of cool water (Kukulka et al. 2010).
To further address the locally enhanced shear instability by LT, we examine the velocity variance profile whose anisotropy and magnitude characterize the turbulence's type and intensity, respectively.
At location A with LT, the normalized hw 02 i has a peak value that is larger than both hu 02 i and hy 02 i at the depth of 0:2H B , which is induced by strong downwelling flows due to LT (Fig. 8, top left). Relatively large hy 02 i near the surface (z . 20:4H B ) is associated with strong LT crosswind convergence regions above downwelling flows (note that x is approximately aligned with the wind) (Thorpe 2004). Below 0.4H B , horizontal velocity variances are larger than hw 02 i, suggesting that shear-driven turbulence is dominant at greater depth (z , 20:4H B ). Enhanced hu 02 i near the bottom of the OSBL indicates enhanced shear-driven turbulence at the OSBL base due to a greater momentum transport, which is expected for rapid mixed layer deepening with LT (Kukulka et al. 2010). In the ST case, normalized hw 02 i in the upper half of the OSBL is much smaller than in the LT case and hu 02 i is always larger than hy 02 i throughout the OSBL as expected for shear-driven turbulence (Fig. 8, top right). In addition, the peak in horizontal velocity variances near the OSBL base is absent.
The profiles of velocity variances in Fig. 8 demonstrate that LT changes the anisotropy of velocity variances, inducing a larger vertical velocity variance in the upper OSBL (z . 20:4H B ). Therefore, the anisotropy of velocity variances is used as a metric to distinguish the regime dominated by LT or ST in the upper part of the OSBL. Here we introduce a turbulent anisotropy ratio R t as the ratio of vertical to horizontal velocity variances following (e.g., Polton and Belcher 2007; Rabe et al. 2015) R t 5 hw 02 i hu 02 i 1 hy 02 i . (9) In the ST case, R t is smaller (,0:3) without much depth variation (Fig. 9, bottom). With LT, we define the depth H LT with R t . 0:5 down to which LT is dominant (Fig. 9, top). Our results indicate that H LT ' 0:4H B for FIG. 7. The profiles of (left) along-wind (solid lines) and crosswind currents (dashed lines), (center) temperature, and (right) gradient Richardson number at (top) location A and (bottom) location B with LT (black lines) and without LT (ST) (red lines). The solid and dashed magenta lines in the left panels denote Lagrangian currents in the along-wind and crosswind direction, respectively. The thin black dashed lines in the right panels denote Ri 5 0:25, and the solid blue lines in the bottom panels indicate the Ekman layer depth H E .
sufficiently large d s /H B (black dashed line in the top panel of Fig. 9). For smaller d s /H B , the Stokes drift decay length controls H LT (Kukulka and Harcourt 2017), and we find H LT ' 6:5d s (magenta dashed line in the top panel of Fig. 9). We find that rapid mixed layer deepening due to LT only occurs for 0:4H B , 6:5d s , which is consistent with Sullivan et al. (2012). In our simulations, this criterion is met before the TC passes, suggesting that LT-enhanced entrainment mainly occurs before the TC passes.
To investigate enhanced turbulent kinetic energy (TKE) production and different mechanisms, we examine the budget of resolved TKE:
The Eulerian shear production P E [first term on the rhs of Eq. ( 10)] and the Stokes drift shear production P S [second term on the rhs of Eq. ( 10)] are two source terms that generate most of the TKE. The remaining terms in Eq. ( 10) from left to right are buoyancy flux (third term), divergence of the energy fluxes due to TKE advection (fourth term), rate of turbulent pressure work (fifth term), and TKE dissipation rate due to subgridscale motion (sixth term). The SGS is the sum of all remaining subgridscale terms (Kukulka et al. 2010;Skyllingstad et al. 2000).
At location A with LT, the Stokes drift shear production is dominant down to H LT (Fig. 10, top left). The Eulerian shear production first decreases with depth, reaching its minimum value at z 5 20:1H B , and then increases with depth and becomes a dominant term at the OSBL base. Augmented P E at the OSBL base is accompanied by greater buoyancy fluxes, consistent with the local peak of horizontal velocity variance (Fig. 8, top left). In contrast, P E for the ST case monotonically decays with depth (top right profile in Fig. 10), which is consistent with the relatively weak buoyancy fluxes and the decaying velocity variance profiles at greater depth (z . 0:6H B ) (top right profile, Fig. 8). In the absence of LT, neither the pronounced downward TKE fluxes in the surface layer nor the local peak of Eulerian shear production near the mixed layer base is present, even during the rapid mixed layer deepening period. Thus, the TKE budgets at location A reveal that LT induces a greater Eulerian shear production within the shear layer by locally enhanced shear instabilities, which overall drives greater buoyancy fluxes, enhancing the buoyancy entrainment.
The projected surface layer Langmuir number La SLu is effective in scaling normalized depth-integrated buoyancy fluxes (see above) and depth-averaged vertical velocity variance under TC conditions with misaligned winds and waves (Rabe et al. 2015;Reichl et al. 2016b). Let us next explore why the scaling with La SLu is consistent with a scaling of the dominant TKE budget terms between P S and h«i (Fig. 10). Following Kukulka and Harcourt (2017), an LT velocity scale w L for an idealized TKE balance between P S and h«i can be expressed as
where l(z) is a depth-dependent turbulent mixing length. Note that for aligned wind and waves, and for l 5 H B (nonlocal mixing), the rhs of Eq. ( 11) is approximately u 3 * La 22 t , consistent with the scaling proposed by Grant and Belcher (2009). Based on Eq. ( 11), Kukulka and Harcourt (2017) introduced a modified Langmuir number La f as
For scaling purposes, let us next assume that the turbulent stress is along the wind direction u w , so that Eq. ( 12) becomes after integration by parts
where g s is the depth-dependent angle of the Stokes drift vector direction. Since l is roughly proportional to z close to the surface and the Stokes drift shear production is predominantly important in the upper 20% of the OSBL (Fig. 10, left; Kukulka and Harcourt 2017), we can further approximate La f as (see appendix)
Thus, La f is proportional to La SLu in Reichl et al. (2016b) if the Lagrangian shear is aligned with the wind direction (a 5 u w ), which is approximately the case for strong LT if the Eulerian current is well mixed (Fig. 1, bottom right). Thus, the successful scaling of VVV and buoyancy fluxes with La SLu can be linked to a leadingorder TKE budget analysis.
To understand the influence of inertial currents on turbulence dynamics, we examine flow properties at location B that is 2.5 RMW behind the TC eye where the ST case shows more energetic buoyancy fluxes and greater mixed layer deepening (Fig. 5, right).
Unlike at location A, the magnitude of hui ? is comparable to hui == and current shear extends below H E , suggesting the presence of energetic inertial currents that drive turbulence at greater depth (Fig. 7, bottom left). As expected for LT, Eulerian currents are nearly uniform and only weakly sheared in the near-surface layer (2z , 6:5d s ' 0:2H B ) that is directly influenced by LT (Fig. 7, bottom left; cf. top panel of Fig. 9). However, the LT-affected surface layer at location B is much shallower than at location A, which does not directly influence the deeper shear-driven layer. At location B, the temperature differences between the ST and LT cases are smaller compared to location A (Fig. 7, middle) so that buoyancy fluxes for the ST case exceed those of the LT case in regions behind the TC eye (Fig. 5, right). Consistently, Ri profiles are similar for both cases in the shear layer with Ri , 0:25 (Fig. 7, bottom right), suggesting that, unlike location A, LT does not influence the deeper layer through locally enhanced shear instabilities.
At location B with LT, the profiles of hu 0 i 2 i near the LT-influenced surface layer (2z , H LT 5 6:5d s ' 0:2H B ) are similar to those at location A for 2z , H LT (Fig. 8, bottom left). At greater depth (2z . H LT ), the horizontal velocity variances exceed hw 02 i and the alongwind velocity variance (approximately hy 02 i) is largest, which is expected for shear-driven turbulence. However, all velocity variances gradually decrease with depth without peak at location B. This is different from location A with LT where LT locally enhances the shear instability near the mixed layer base, resulting in locally enhanced along-wind velocity variances. In the shear layer, velocity variances are significantly larger for the ST than for the LT case, suggesting relatively strong shear turbulence and Eulerian shear production that contributes to the greater entrainment in the ST case.
As for location A, the leading-order budget terms at location B are TKE production and dissipation (Fig. 10, bottom). Unlike location A, the normalized TKE production is now substantially enhanced for the ST case at greater depth, so that the relative magnitude of TKE production terms evolve over time (Fig. 11, bottom right). Near the surface, TKE production remains relatively large because of strong Lagrangian nearsurface shear for both ST and LT cases. In the LT case, the total TKE production (P S 1 P E ) becomes smaller than P E for the ST case for (X 2 X 0 )/RMW . 1 (Fig. 11, right). Since current shear is similar for both LT and ST cases (Fig. 7, bottom left), it is the relatively large turbulent stress in the shear layer (z , 20:6H M ) that causes this enhanced P E in the ST case (Fig. 3 Second, enhanced turbulent near-surface mixing reduces the surface currents and, thus, the surface flux of kinetic energy into the OSBL, which is discussed in the following section.
To understand the relation between turbulence generation and the evolution of inertial currents, we investigate the depth-integrated total energy budget at (Y 2 Y 0 )/RMW 5 1, based on the governing LES equations:
where « ikm is the Levi-Civita tensor and (f 1 , f 2 , f 3 ) 5 (0, 0, f ) is the Coriolis vector. The t i and hu i i 0 denote the wind stress vector and surface current vector, respectively. Note that we are considering the depthintegrated energy budget below the wave-breaking layer and assume that breaking energy input is all dissipated locally, consistent with previous LES models under TC (Sullivan et al. 2012;Rabe et al. 2015;Reichl et al. 2016b). The surface flux to mean energy F ME [first term on the rhs of Eq. ( 15)], the depth-integrated Stokes drift shear production I Stokes [second term on the rhs of Eq. ( 15)], and depth-integrated dissipation rate D [third term on the rhs of Eq. ( 15)] are source and sink terms that drive the temporal change of the depthintegrated energy of the system, which partitions into mean kinetic energy (MKE) [first term on the lhs of Eq. ( 15)], TKE [second term on the lhs of Eq. ( 15)], and potential energy (PE) [third term on the lhs of Eq. ( 15)]. The last lhs term of Eq. ( 15) T CS transfers energy between waves and Eulerian currents due to the Coriolis-Stokes force. This term is much smaller than other budget terms and, therefore, neglected in the following presentation. With these symbols, Eq. ( 15) is written as
Throughout the passage of the storm, F ME (black dashed line in Fig. 12a) and, with LT, I Stokes (solid cyan line in Fig. 12a) are predominantly balanced by D (solid and dashed green lines in Fig. 12a) and changes in mean kinetic energy (solid and dashed red lines in Fig. 12a), which are similar in magnitude. The changes in turbulent kinetic energy (solid and dashed blue lines in Fig. 12a) are small compared to other budget terms. Normalized F ME steadily increases as hu i i 0 increases owing to strong inertial currents. Because inertial currents dissipate relatively slowly, F ME remains relatively large even after the storm passes. With LT, I Stokes is comparable to F ME for (X 2 X 0 )/RMW , 0. 5, but I Stokes is much smaller than F ME behind the eye for (X 2 X 0 )/RMW . 1, when energetic waves propagate away from the storm but inertial currents are still strong.
With LT, the total energy budget dynamics changes at about (X 2 X 0 )/RMW 5 0:5, when the LT dynamics change from the rapid mixed layer deepening regime to the regime for which LT is relatively shallow and the LT-driven surface layer is decoupled from the deeper shear layer [refer to section 3c(1) and top panel of Fig. 9]. For (X 2 X 0 )/RMW , 0:5, that is, during rapid mixed layer deepening, energy is more efficiently used to entrain deeper cool water and to change the potential energy of the system more quickly compared to the ST case (solid magenta line in Fig. 12b). Furthermore, less energy is transferred to inertial currents because in the LT case currents are more mixed vertically and greater vertical shear in the ST case results in relatively large MKE and d(MKE)/dt (red line in Fig. 12b).
For (X 2 X 0 )/RMW . 0:5, on the other hand, the relatively shallow LT-influenced near-surface layer causes currents to be more sheared at greater depth inducing a relatively large transfer to MKE, so that MKE changes more in the LT case than the ST case during this period (solid red line in Fig. 12b). Nevertheless, LT mixes near-surface currents so that hu i i 0 is still relatively small, which is accompanied by relatively small F ME compared to the ST case (solid and dashed black lines in Fig. 12a). Thus, the presence of LT weakens the generation of turbulence at greater depth in the shear layer, so that the entrainment and the associated change in potential energy is greater for the ST case once the winds subside (solid magenta line in Fig. 12b).
Based on the analysis of high-resolution large-eddy simulation (LES) results, we have investigated the response of the ocean surface boundary layer (OSBL) with and without Langmuir turbulence (LT) to complex wind and wave forcing in tropical cyclone (TC) conditions. The Stokes drift vector that drives LT through the Craik-Leibovich vortex force is determined from spectral wave simulations obtained from the WAVEWATCH III model. Both the LES and the wave models are forced by winds from the TC Holland wind model. LES experiments are performed with and without LT (i.e., without Stokes drift) and for multiple, densely spaced stations across the TC translation direction, covering the full spatial extent of the TC.
The examination of OSBL depth and entrainment illustrates that LT substantially enhances entrainment with greater buoyancy fluxes, inducing rapid OSBL deepening during maximum TC winds. The mechanism that drives rapid OSBL deepening is that LT facilitates the downward momentum transport, locally enhancing shear instabilities near the mixed layer base, which occurs only when OSBL depth H B is relatively shallow compared to a Stokes drift decay length scale d s . LT enhances vertical velocity variances near a surface layer whose depth H LT is estimated from the enhanced anisotropy of velocity variances. During the rapid OSBL deepening, the LT-affected surface layer is tightly coupled with the deeper shear-driven layer so that H LT is proportional to H B . As a result, depth-integrated entrainment fluxes approximately scale with the projected surface layer Langmuir number La SLu , which is consistent with an idealized TKE budget analysis for misaligned wind and waves.
After the TC passes, H LT is proportional to d s and the LT-affected surface layer influences relatively weakly the deeper shear-driven layer. In the LT case, the OSBL is characterized by relatively small Lagrangian shear at middepth. At the same time, the normalized surface energy flux to mean currents increases substantially behind the TC eye, exceeding the energy input by the depth-integrated Stokes drift shear production. This is because of strong inertial currents that have been excited by TC winds. LT still homogenizes Eulerian currents in the LT-affected surface layer so that surface currents are reduced in the LT case. The reduced surface currents in the LT case, in turn, reduce the surface energy flux to mean currents, resulting in a smaller total energy input than for the ST case. Here we assume the wind stress and related surface energy fluxes are transferred through the wave breaking layer within which most of the energy input due to breaking waves is all dissipated locally. Together, the small Lagrangian shear at middepth and the reduced total energy input contribute to relatively small TKE levels and entrainment rates for the LT case behind the TC's passage. Therefore, our study illustrates that background currents need to be taken into account for a complete scaling of LT and its effects on the OSBL under TC conditions.
La 22 f ; ð 0 20:2H B ju s j cos(g s 2 u w ) 0:2H B u * dz 5 ð 0 20:2H B ju s j cos(g s 2 u w ) dz 0:2H B u * . (A4)
By comparing Eq. (A4) to Eq. (A3), we show that La SLu is proportional to La f .
J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 48
SEPTEMBER 2018W A N G E T A L .