A subgrid‐scale eddy parameterization is developed, which makes use of an explicit eddy kinetic energy budget and can be applied at both “non‐eddying” and “eddy‐permitting” resolutions. The subgrid‐scale eddies exchange energy with the resolved flow in both directions via a parameterization of baroclinic instability (based on the established formulation of Gent and McWilliams) and biharmonic and negative Laplacian viscosity terms. This formulation represents the turbulent cascades of energy and enstrophy consistent with our current understanding of the turbulent eddy energy cycle. At the same time, the approach is simple and general enough to be readily implemented in ocean climate models, without adding significant computational cost. The closure has been implemented in the Modular Ocean Model Version 6 and tested in the “Neverworld” configuration, which employs an idealized analytically defined topography designed as a testbed for mesoscale eddy parameterizations. The parameterization performs well over a range of resolutions, seamlessly connecting non‐eddying and eddy‐resolving regimes.
Diagnosing Scale-Dependent Energy Cycles in a High-Resolution Isopycnal Ocean Model
Abstract Energy exchanges between large-scale ocean currents and mesoscale eddies play an important role in setting the large-scale ocean circulation but are not fully captured in models. To better understand and quantify the ocean energy cycle, we apply along-isopycnal spatial filtering to output from an isopycnal 1/32° primitive equation model with idealized Atlantic and Southern Ocean geometry and topography. We diagnose the energy cycle in two frameworks: 1) a non-thickness-weighted framework, resulting in a Lorenz-like energy cycle, and 2) a thickness-weighted framework, resulting in the Bleck energy cycle. This paper shows that framework 2 is more useful for studying energy pathways when an isopycnal average is used. Next, we investigate the Bleck cycle as a function of filter scale. Baroclinic conversion generates mesoscale eddy kinetic energy over a wide range of scales and peaks near the deformation scale at high latitudes but below the deformation scale at low latitudes. Away from topography, an inverse cascade transfers kinetic energy from the mesoscales to larger scales. The upscale energy transfer peaks near the energy-containing scale at high latitudes but below the deformation scale at low latitudes. Regions downstream of topography are characterized by a downscale kinetic energy transfer, in which mesoscale eddies are generated through barotropic instability. The scale- and flow-dependent energy pathways diagnosed in this paper provide a basis for evaluating and developing scale- and flow-aware mesoscale eddy parameterizations. Significance Statement Blowing winds provide a major energy source for the large-scale ocean circulation. A substantial fraction of this energy is converted to smaller-scale eddies, which swirl through the ocean as sea cyclones. Ocean turbulence causes these eddies to transfer part of their energy back to the large-scale ocean currents. This ocean energy cycle is not fully simulated in numerical models, but it plays an important role in transporting heat, carbon, and nutrients throughout the world’s oceans. The purpose of this study is to quantify the ocean energy cycle by using fine-scale idealized numerical simulations of the Atlantic and Southern Oceans. Our results provide a basis for how to include unrepresented energy exchanges in coarse global climate models.
more »
« less
- PAR ID:
- 10426875
- Date Published:
- Journal Name:
- Journal of Physical Oceanography
- Volume:
- 53
- Issue:
- 1
- ISSN:
- 0022-3670
- Page Range / eLocation ID:
- 157 to 176
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Abstract Oceanic mesoscale motions including eddies, meanders, fronts, and filaments comprise a dominant fraction of oceanic kinetic energy and contribute to the redistribution of tracers in the ocean such as heat, salt, and nutrients. This reservoir of mesoscale energy is regulated by the conversion of potential energy and transfers of kinetic energy across spatial scales. Whether and under what circumstances mesoscale turbulence precipitates forward or inverse cascades, and the rates of these cascades, remain difficult to directly observe and quantify despite their impacts on physical and biological processes. Here we use global observations to investigate the seasonality of surface kinetic energy and upper-ocean potential energy. We apply spatial filters to along-track satellite measurements of sea surface height to diagnose surface eddy kinetic energy across 60–300-km scales. A geographic and scale-dependent seasonal cycle appears throughout much of the midlatitudes, with eddy kinetic energy at scales less than 60 km peaking 1–4 months before that at 60–300-km scales. Spatial patterns in this lag align with geographic regions where an Argo-derived estimate of the conversion of potential to kinetic energy is seasonally varying. In midlatitudes, the conversion rate peaks 0–2 months prior to kinetic energy at scales less than 60 km. The consistent geographic patterns between the seasonality of potential energy conversion and kinetic energy across spatial scale provide observational evidence for the inverse cascade and demonstrate that some component of it is seasonally modulated. Implications for mesoscale parameterizations and numerical modeling are discussed. Significance Statement This study investigates the seasonality of upper-ocean potential and kinetic energy in the context of an inverse cascade, consisting of energy transfers to and through the mesoscale. Observations show a scale-dependent cycle in kinetic energy that coincides with temporal variability in mixed layer potential energy and progresses seasonally from smaller to larger scales. This pattern appears dominant over large regions of the ocean. Results are relevant to ocean and climate models, where a large fraction of ocean energy is often parameterized. A customizable code repository and dataset are provided to enable comparisons of model-based resolved and unresolved kinetic energy to observational equivalents. Implications result for a range of processes including mixed layer stratification and vertical structure of ocean currents.more » « less
-
Abstract The Beaufort Gyre (BG) is hypothesized to be partially equilibrated by those mesoscale eddies that form via baroclinic instabilities of its currents. However, our understanding of the eddy field’s dependence on the mean BG currents and the role of sea ice remains incomplete. This theoretical study explores the scales and vertical structures of eddies forming specifically due to baroclinic instabilities of interior BG flows. An idealized quasi-geostrophic model is used to show that flows driven only by the Ekman pumping contain no interior potential vorticity (PV) gradients and generate weak and large eddies, ℴ(200km) in size, with predominantly barotropic and first baroclinic mode energy. However, flows containing realistic interior PV gradients in the Pacific halocline layer generate significantly smaller eddies of about 50 km in size, with a distinct second baroclinic mode structure and a subsurface kinetic energy maximum. The dramatic change in eddy characteristics is shown to be caused by the stirring of interior PV gradients by large-scale barotropic eddies. The sea ice-ocean drag is identified as the dominant eddy dissipation mechanism, leading to realistic sub-surface maxima of eddy kinetic energy for drag coefficients higher than about 2×10 −3 . A scaling law is developed for the eddy potential enstrophy, demonstrating that it is directly proportional to the interior PV gradient and the square root of the barotropic eddy kinetic energy. This study proposes a possible formation mechanism of large BG eddies and points to the importance of accurate representation of the interior PV gradients and eddy dissipation by ice-ocean drag in BG simulations and theory.more » « less
-
Abstract. We describe an idealized primitive-equation model for studying mesoscale turbulence and leverage a hierarchy of grid resolutions to make eddy-resolving calculations on the finest grids more affordable.The model has intermediate complexity, incorporating basin-scale geometry with idealized Atlantic and Southern oceans and with non-uniform ocean depth to allow for mesoscale eddy interactions with topography.The model is perfectly adiabatic and spans the Equator and thus fills a gap between quasi-geostrophic models, which cannot span two hemispheres, and idealized general circulation models, which generally include diabatic processes and buoyancy forcing.We show that the model solution is approaching convergence in mean kinetic energy for the ocean mesoscale processes of interest and has a rich range of dynamics with circulation features that emerge only due to resolving mesoscale turbulence.more » « less
-
Abstract The Atlantic meridional overturning circulation (AMOC) plays a key role in climate due to uptake and redistribution of heat and carbon anomalies. This redistribution takes place along several main pathways that link the high-latitude North Atlantic with midlatitudes and the Southern Ocean and involves currents on a wide range of spatial scales. This numerical study examines the importance of mesoscale currents (“eddies”) in these AMOC pathways and associated time scales, using a highly efficient offline tracer model. The study uses two boundary impulse response (BIR) tracers, which can quantify the importance of the Atlantic tracer exchanges with the high-latitude atmosphere in the north and with the Southern Ocean in the south. The results demonstrate that mesoscale advection leads to an increase in the overall BIR inventory during the first 100 years and results in a more efficient and spatially uniform ventilation of the deep Atlantic. Mesoscale currents also facilitate meridional spreading of the BIR tracer and thus assist the large-scale advection. The results point toward the importance of spatial inhomogeneity and anisotropy of the eddy-induced mixing in several mixing “hotspots,” as revealed by an eddy diffusivity tensor. Conclusions can be expected to assist evaluations of eddy-permitting simulations that stop short of full resolution of mesoscale, as well as development of eddy parameterization schemes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JPO-D-22-0083.1
