skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Deep-Cycle Turbulence in the Upper Pacific Equatorial Ocean: Characterization by LES and Heat Flux Parameterization
Abstract 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
Award ID(s):
1851390 2048631
PAR ID:
10534513
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
https://journals.ametsoc.org/view/journals/phoc/54/2/JPO-D-23-0015.1.xml
Date Published:
Journal Name:
Journal of Physical Oceanography
Volume:
54
Issue:
2
ISSN:
0022-3670
Page Range / eLocation ID:
577 to 599
Format(s):
Medium: X
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. 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
  3. null (Ed.)
    Abstract This study investigates the dynamics of velocity shear and Reynolds stress in the ocean surface boundary layer for idealized misaligned wind and wave fields using a large-eddy simulation (LES) model based on the Craik–Leibovich equations, which captures Langmuir turbulence (LT). To focus on the role of LT, the LES experiments omit the Coriolis force, which obscures a stress–current-relation analysis. Furthermore, a vertically uniform body force is imposed so that the volume-averaged Eulerian flow does not accelerate but is steady. All simulations are first spun-up without wind-wave misalignment to reach a fully developed stationary turbulent state. Then, a crosswind Stokes drift profile is abruptly imposed, which drives crosswind stresses and associated crosswind currents without generating volume-averaged crosswind currents. The flow evolves to a new stationary state, in which the crosswind Reynolds stress vanishes while the crosswind Eulerian shear and Stokes drift shear are still present, yielding a misalignment between Reynolds stress and Lagrangian shear (sum of Eulerian current and Stokes drift). A Reynolds stress budgets analysis reveals a balance between stress production and velocity–pressure gradient terms (VPG) that encloses crosswind Eulerian shear, demonstrating a complex relation between shear and stress. In addition, the misalignment between Reynolds stress and Eulerian shear generates a horizontal turbulent momentum flux (due to correlations of along-wind and crosswind turbulent velocities) that can be important in producing Reynolds stress (due to correlations of horizontal and vertical turbulent velocities). Thus, details of the Reynolds stress production by Eulerian and Stokes drift shear may be critical for driving upper-ocean currents and for accurate turbulence parameterizations in misaligned wind-wave conditions. 
    more » « less
  4. The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance StatementThe traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence. 
    more » « less
  5. AMS (Ed.)
    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 StatementOnly 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