skip to main content

Title: Deep-Cycle Turbulence in the Upper Pacific Equatorial Ocean: Characterization by LES and Heat Flux Parameterization

Observations in the Pacific Equatorial Undercurrents (EUC) show that the nighttime deep-cycle turbulence (DCT) in the marginal-instability (MI) layer of the EUC exhibits seasonal variability that can modulate heat transport and sea surface temperature. Large-eddy simulations (LES), spanning a wide range of control parameters, are performed to identify the key processes that influence the turbulent heat flux at multiple time scales ranging from turbulent (minutes to hours) to daily to seasonal. The control parameters include wind stress, convective surface heat flux, shear magnitude, and thickness of the MI layer. In the LES, DCT occurs in discrete bursts during the night, exhibits high temporal variability within a burst, and modulates the mixed layer depth. At the daily time scale, turbulent heat flux generally increases with increasing wind stress, MI-layer shear, or nighttime convection. Convection is found to be important to mixing under weak wind, weak shear conditions. A parameterization for the daily averaged turbulent heat flux is developed from the LES suite to infer the variability of heat flux at the seasonal time scale. The LES-based parameterized heat flux, which takes into account the effects of all control parameters, exhibits a seasonal variability that is similar to the observed heat flux from theχ-pods.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Physical Oceanography
Medium: X Size: p. 577-599
["p. 577-599"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Microstructure observations in the Pacific cold tongue reveal that turbulence often penetrates into the thermocline, producing hundreds of watts per square meter of downward heat transport during nighttime and early morning. However, virtually all observations of this deep-cycle turbulence (DCT) are from 0°, 140°W. Here, a hierarchy of ocean process simulations, including submesoscale-permitting regional models and turbulence-permitting large-eddy simulations (LES) embedded in a regional model, provide insight into mixing and DCT at and beyond 0°, 140°W. A regional hindcast quantifies the spatiotemporal variability of subsurface turbulent heat fluxes throughout the cold tongue from 1999 to 2016. Mean subsurface turbulent fluxes are strongest (∼100 W m −2 ) within 2° of the equator, slightly (∼10 W m −2 ) stronger in the northern than Southern Hemisphere throughout the cold tongue, and correlated with surface heat fluxes ( r 2 = 0.7). The seasonal cycle of the subsurface heat flux, which does not covary with the surface heat flux, ranges from 150 W m −2 near the equator to 30 and 10 W m −2 at 4°N and 4°S, respectively. Aseasonal variability of the subsurface heat flux is logarithmically distributed, covaries spatially with the time-mean flux, and is highlighted in 34-day LES of boreal autumn at 0° and 3°N, 140°W. Intense DCT occurs frequently above the undercurrent at 0° and intermittently at 3°N. Daily mean heat fluxes scale with the bulk vertical shear and the wind stress, which together explain ∼90% of the daily variance across both LES. Observational validation of the scaling at 0°, 140°W is encouraging, but observations beyond 0°, 140°W are needed to facilitate refinement of mixing parameterization in ocean models. Significance Statement This work is a fundamental contribution to a broad community effort to improve global long-range weather and climate forecast models used for seasonal to longer-term prediction. Much of the predictability on seasonal time scales is derived from the slow evolution of the upper eastern equatorial Pacific Ocean as it varies between El Niño and La Niña conditions. This study presents state-of-the-art high-resolution regional numerical simulations of ocean turbulence and mixing in the eastern equatorial Pacific. The results inform future planning for field work as well as future efforts to refine the representation of ocean mixing in global forecast models. 
    more » « less
  2. Abstract

    Multiyear turbulence measurements from oceanographic moorings in equatorial Atlantic and Pacific cold tongues reveal similarities in deep cycle turbulence (DCT) beneath the mixed layer (ML) and above the Equatorial Undercurrent (EUC) core. Diurnal composites of turbulence kinetic energy dissipation rate,ϵ, clearly show the diurnal cycles of turbulence beneath the ML in both cold tongues. Despite differences in surface forcing, EUC strength and core depth DCT occurs, and is consistent in amplitude and timing, at all three sites. Time‐mean values ofϵat 30 m depth are nearly identical at all three sites. Variations of averaged values ofϵin the deep cycle layer below 30 m range to a factor of 10 between sites. A proposed scaling in depth that isolates the deep cycle layers and ofϵby the product of wind stress and current shear collapses vertical profiles at all sites to within a factor of 2.

    more » « less
  3. This study utilizes a large-eddy simulation (LES) approach to systematically assess the directional variability of wave-driven Langmuir turbulence (LT) in the ocean surface boundary layer (OSBL) under tropical cyclones (TCs). The Stokes drift vector, which drives LT through the Craik–Leibovich vortex force, is obtained through spectral wave simulations. LT’s direction is identified by horizontally elongated turbulent structures and objectively determined from horizontal autocorrelations of vertical velocities. In spite of a TC’s complex forcing with great wind and wave misalignments, this study finds that LT is approximately aligned with the wind. This is because the Reynolds stress and the depth-averaged Lagrangian shear (Eulerian plus Stokes drift shear) that are key in determining the LT intensity (determined by normalized depth-averaged vertical velocity variances) and direction are also approximately aligned with the wind relatively close to the surface. A scaling analysis of the momentum budget suggests that the Reynolds stress is approximately constant over a near-surface layer with predominant production of turbulent kinetic energy by Stokes drift shear, which is confirmed from the LES results. In this layer, Stokes drift shear, which dominates the Lagrangian shear, is aligned with the wind because of relatively short, wind-driven waves. On the contrary, Stokes drift exhibits considerable amount of misalignments with the wind. This wind–wave misalignment reduces LT intensity, consistent with a simple turbulent kinetic energy model. Our analysis shows that both the Reynolds stress and LT are aligned with the wind for different reasons: the former is dictated by the momentum budget, while the latter is controlled by wind-forced waves.

    more » « less
  4. Abstract

    Several years of moored turbulence measurements fromχpods at three sites in the equatorial cold tongues of Atlantic and Pacific Oceans yield new insights into proxy estimates of turbulence that specifically target the cold tongues. They also reveal previously unknown wind dependencies of diurnally varying turbulence in the near-critical stratified shear layers beneath the mixed layer and above the core of the Equatorial Undercurrent that we have come to understand as deep cycle (DC) turbulence. Isolated by the mixed layer above, the DC layer is only indirectly linked to surface forcing. Yet, it varies diurnally in concert with daily changes in heating/cooling. Diurnal composites computed from 10-min averaged data at fixedχpod depths show that transitions from daytime to nighttime mixing regimes are increasingly delayed with weakening wind stressτ. These transitions are also delayed with respect to depth such that they follow a descent rate of roughly 6 m h−1, independent ofτ. We hypothesize that this wind-dependent delay is a direct result of wind-dependent diurnal warm layer deepening, which acts as the trigger to DC layer instability by bringing shear from the surface downward but at rates much slower than 6 m h−1. This delay in initiation of DC layer instability contributes to a reduction in daily averaged values of turbulence dissipation. Both the absence of descending turbulence in the sheared DC layer prior to arrival of the diurnal warm layer shear and the magnitude of the subsequent descent rate after arrival are roughly predicted by laboratory experiments on entrainment in stratified shear flows.

    Significance Statement

    Only recently have long time series measurements of ocean turbulence been available anywhere. Important sites for these measurements are the equatorial cold tongues where the nature of upper-ocean turbulence differs from that in most of the world’s oceans and where heat uptake from the atmosphere is concentrated. Critical to heat transported downward from the mixed layer is the diurnally varying deep cycle of turbulence below the mixed layer and above the core of the Equatorial Undercurrent. Even though this layer does not directly contact the surface, here we show the influence of the surface winds on both the magnitude of the deep cycle turbulence and the timing of its descent into the depths below.

    more » « less
  5. Large-eddy simulation (LES) is used to model turbulent winds in a nominally neutral atmospheric boundary layer at varying mesh resolutions. The boundary layer is driven by wind shear with zero surface heat flux and is capped by a stable inversion. Because of entrainment the boundary layer is in a weakly stably stratified regime. The simulations use meshes varying from 1282× 64 to 10242× 512 grid points in a fixed computational domain of size (2560, 2560, 896) m. The subgrid-scale (SGS) parameterizations used in the LES vary with the mesh spacing. Low-order statistics, spectra, and structure functions are compared on the different meshes and are used to assess grid convergence in the simulations. As expected, grid convergence is primarily achieved in the middle of the boundary layer where there is scale separation between the energy-containing and dissipative eddies. Near the surface second-order statistics do not converge on the meshes studied. The analysis also highlights differences between one-dimensional and two-dimensional velocity spectra; differences are attributed to sampling errors associated with aligning the horizontal coordinates with the vertically veering mean wind direction. Higher-order structure functions reveal non-Gaussian statistics on all scales, but are highly dependent on the mesh resolution. A generalized logarithmic law and a k−1spectral scaling regime are identified with mesh-dependent parameters in agreement with previously published results.

    more » « less