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.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Vertical Structure and Energetic Constraints for a Backscatter Parameterization of Ocean Mesoscale Eddies
Abstract Mesoscale eddies modulate the stratification, mixing, tracer transport, and dissipation pathways of oceanic flows over a wide range of spatiotemporal scales. The parameterization of buoyancy and momentum fluxes associated with mesoscale eddies thus presents an evolving challenge for ocean modelers, particularly as modern climate models approach eddy‐permitting resolutions. Here we present a parameterization targeting such resolutions through the use of a subgrid mesoscale eddy kinetic energy budget (MEKE) framework. Our study presents two novel insights: (a) both the potential and kinetic energy effects of eddies may be parameterized via a kinetic energy backscatter, with no Gent‐McWilliams along‐isopycnal transport; (b) a dominant factor in ensuring a physically‐accurate backscatter is the vertical structure of the parameterized momentum fluxes. We present simulations of 1/2° and 1/4° resolution idealized models with backscatter applied to the equivalent barotropic mode. Remarkably, the global kinetic and potential energies, isopycnal structure, and vertical energy partitioning show significantly improved agreement with a 1/32° reference solution. Our work provides guidance on how to parameterize mesoscale eddy effects in the challenging eddy‐permitting regime.  more » « less
Award ID(s):
1912357
PAR ID:
10648397
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
16
Issue:
7
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Atmospheric research)
    The surface wind structure and vertical turbulent transport processes in the eyewall of hurricane Isabel (2003) are investigated using six large-eddy simulations (LESs) with different horizontal grid spacing and three-dimensional (3D) sub-grid scale (SGS) turbulent mixing models and a convection permitting simulation that uses a coarser grid spacing and one-dimensional vertical turbulent mixing scheme. The mean radius-height distribution of storm tangential wind and radial flow, vertical velocity structure, and turbulent kinetic energy and momentum fluxes in the boundary layer generated by LESs are consistent with those derived from historical dropsonde composites, Doppler radar, and aircraft measurements. Unlike the convection permitting simulation that produces storm wind fields lacking small-scale disturbances, all LESs are able to produce sub-kilometer and kilometer scale eddy circulations in the eyewall. The inter-LES differences generally reduce with the decrease of model grid spacing. At 100-m horizontal grid spacing, the vertical momentum fluxes induced by the model-resolved eddies and the associated eddy exchange coefficients in the eyewall simulated by the LESs with different 3D SGS mixing schemes are fairly consistent. Although with uncertainties, the decomposition in terms of eddy scales suggests that sub-kilometer eddies are mainly responsible for the vertical turbulent transport within the boundary layer (~1 km depth following the conventional definition) whereas eddies greater than 1 km become the dominant contributors to the vertical momentum transport above the boundary layer in the eyewall. The strong dependence of vertical turbulent transport on eddy scales suggests that the vertical turbulent mixing parameterization in mesoscale simulations of tropical cyclones is ultimately a scale-sensitive problem. 
    more » « less
  2. ; ; ; ; (, Journal of Advances in Modeling Earth Systems)
    Ocean mesoscale eddies are often poorly represented in climate models, and therefore, their effects on the large scale circulation must be parameterized. Traditional parameterizations, which represent the bulk effect of the unresolved eddies, can be improved with new subgrid models learned directly from data. Zanna and Bolton (ZB20) applied an equation‐discovery algorithm to reveal an interpretable expression parameterizing the subgrid momentum fluxes by mesoscale eddies through the components of the velocity‐gradient tensor. In this work, we implement the ZB20 parameterization into the primitive‐equation GFDL MOM6 ocean model and test it in two idealized configurations with significantly different dynamical regimes and topography. The original parameterization was found to generate excessive numerical noise near the grid scale. We propose two filtering approaches to avoid the numerical issues and additionally enhance the strength of large‐scale energy backscatter. The filtered ZB20 parameterizations led to improved climatological mean state and energy distributions, compared to the current state‐of‐the‐art energy backscatter parameterizations. The filtered ZB20 parameterizations are scale‐aware and, consequently, can be used with a single value of the non‐dimensional scaling coefficient for a range of resolutions. The successful application of the filtered ZB20 parameterizations to parameterize mesoscale eddies in two idealized configurations offers a promising opportunity to reduce long‐standing biases in global ocean simulations in future studies. 
    more » « less
  3. ; ; ; ; ; ; ; (, Journal of Advances in Modeling Earth Systems)
    Abstract There are two distinct parameterizations for the restratification effect of mesoscale eddies: the Greatbatch and Lamb (1990, GL90,https://journals.ametsoc.org/view/journals/phoc/20/10/1520-0485_1990_020_1634_opvmom_2_0_co_2.xml?tab_body=abstract-display) parameterization, which mixes horizontal momentum in the vertical, and the Gent and McWilliams (1990, GM90,https://journals.ametsoc.org/view/journals/phoc/20/1/1520-0485_1990_020_0150_imiocm_2_0_co_2.xml) parameterization, which flattens isopycnals adiabatically. Even though these two parameterizations are effectively equivalent under the assumption of quasi‐geostrophy, GL90 has been used much less than GM90, and exclusively inz‐coordinate models. In this paper, we compare the GL90 and GM90 parameterizations in an idealized isopycnal coordinate model, both from a theoretical and practical perspective. From a theoretical perspective, GL90 is more attractive than GM90 for isopycnal coordinate models because GL90 provides an interpretation that is fully consistent with thickness‐weighted isopycnal averaging, while GM90 cannot be entirely reconciled with any fully isopycnal averaging framework. From a practical perspective, the GL90 and GM90 parameterizations lead to extremely similar energy levels, flow and vertical structure, even though their energetic pathways are very different. The striking resemblance between the GL90 and GM90 simulations persists from non‐eddying through eddy‐permitting resolution. We conclude that GL90 is a promising alternative to GM90 for isopycnal coordinate models, where it is more consistent with theory, computationally more efficient, easier to implement, and numerically more stable. Assessing the applicability of GL90 in realistic global ocean simulations with hybrid coordinate schemes should be a priority for future work. 
    more » « less
  4. ; ; ; (, Journal of Physical Oceanography)
    Abstract Surface and upper-ocean measurements of mesoscale eddies have revealed the central role they play in ocean transport, but their interior and deep ocean characteristics remain undersampled and underexplored. In this study, mooring arrays, sampling with high vertical resolution, and a high-resolution global atmosphere–ocean coupled simulation are used to characterize full-depth mesoscale eddy vertical structure. The vertical structure of eddy kinetic energy, e.g., partitioning of barotropic to baroclinic eddy kinetic energy or vertical modal structure, is shown to depend partly on bathymetric slope and roughness. This influence is contextualized alongside additional factors, such as latitude and vertical density stratification, to present a global landscape of vertical structure. The results generally reveal eddy vertical structure to decay with increasing depth, consistent with theoretical expectations relating to the roles of surface-intensified stratification and buoyancy anomalies. However, at high latitudes and where the seafloor is markedly flat and smooth (approximately 20% of the ocean’s area), mesoscale eddy vertical structures are significantly more barotropic by an approximate factor of 2–5. From a climate modeling perspective, these results can inform the construction, implementation, and improvement of energetic parameterizations that account for the underrepresentation of mesoscale eddies and their effects. They also offer expectation as to a landscape of eddy vertical structure to be used in inferring vertical structure from surface measurements. Significance StatementThis work addresses the question of how do ocean seafloor features (bathymetry) affect the vertical structure of ocean currents and eddies? Seafloor features modify eddies in complex ways not often accounted for in global ocean simulations. We analyze high-resolution velocity observations, find diverse structures at four mooring sites, and consider how sloping and rough bathymetry change distributions of eddy kinetic energy throughout the water column. Comparison to theory and model output reveals a relationship between vertical structure and bathymetry. These results show that vertical structures vary significantly with bathymetry, density stratification, and latitude and contribute to model development efforts to reproduce the effects of eddy turbulence without explicit representation. These results also enhance interpretations of more numerous surface observations. 
    more » « less
  5. ; ; (, Journal of Physical Oceanography)
    null (Ed.)
    Abstract Theories of the Beaufort Gyre (BG) dynamics commonly represent the halocline as a single layer with a thickness depending on the Eulerian-mean and eddy-induced overturning. However, observations suggest that the isopycnal slope increases with depth, and a theory to explain this profile remains outstanding. Here we develop a multilayer model of the BG, including the Eulerian-mean velocity, mesoscale eddy activity, diapycnal mixing, and lateral boundary fluxes, and use it to investigate the dynamics within the Pacific Winter Water (PWW) layer. Using theoretical considerations, observational data, and idealized simulations, we demonstrate that the eddy overturning is critical in explaining the observed vertical structure. In the absence of the eddy overturning, the Ekman pumping and the relatively weak vertical mixing would displace isopycnals in a nearly parallel fashion, contrary to observations. This study finds that the observed increase of the isopycnal slope with depth in the climatological state of the gyre is consistent with a Gent–McWilliams eddy diffusivity coefficient that decreases by at least 10%–40% over the PWW layer. We further show that the depth-dependent eddy diffusivity profile can explain the relative magnitude of the correlated isopycnal depth and layer thickness fluctuations on interannual time scales. Our inference that the eddy overturning generates the isopycnal layer thickness gradients is consistent with the parameterization of eddies via a Gent–McWilliams scheme but not potential vorticity diffusion. This study implies that using a depth-independent eddy diffusivity, as is commonly done in low-resolution ocean models, may contribute to misrepresentation of the interior BG dynamics. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email