Search

Slowly evolving stratified flow over rough topography is subject to substantial drag due to internal motions, but often numerical simulations are carried out at resolutions where this “wave” drag must be parameterized. Here we highlight the importance of internal drag from topography with scales that cannot radiate internal waves, but may be highly nonlinear, and we propose a simple parameterization of this drag that has a minimum of fit parameters compared to existing schemes. The parameterization smoothly transitions from a quadratic drag law () for low Nh/u0 (linear wave dynamics) to a linear drag law () for high Nh/u0 flows (nonlinear blocking and hydraulic dynamics), where N is the stratification, h is the height of the topography, and u0 is the near-bottom velocity; the parameterization does not have a dependence on Coriolis frequency. Simulations carried out in a channel with synthetic bathymetry and steady body forcing indicate that this parameterization accurately predicts drag across a broad range of forcing parameters when the effect of reduced near-bottom mixing is taken into account by reducing the effective height of the topography. The parameterization is also tested in simulations of wind-driven channel flows that generate mesoscale eddy fields, a setup where the downstream transport is sensitive to the bottom drag parameterization and its effect on the eddies. In these simulations, the parameterization replicates the effect of rough bathymetry on the eddies. If extrapolated globally, the subinertial topographic scales can account for 2.7 TW of work done on the low-frequency circulation, an important sink that is redistributed to mixing in the open ocean.

The thickness-weighted average (TWA) framework, which treats the residual-mean flow as the prognostic variable, provides a clear theoretical formulation of the eddy feedback onto the residual-mean ...

We investigate the role of small‐scale, high‐frequency motions on lateral transport in the ocean, by using velocity fields and particle trajectories from an ocean general circulation model (MITgcm‐llc4320) that permits submesoscale flows, inertia‐gravity waves, and tides. Temporal averaging/filtering removes most of the submesoscale turbulence, inertia‐gravity waves, and tides, resulting in a largely geostrophic flow, with a rapid drop‐off in energy at scales smaller than the mesoscales. We advect two types of Lagrangian particles: (a) 2‐D particles (surface restricted) and (b) 3‐D particles (advected in full three dimensions) with the filtered and unfiltered velocities and calculate Lagrangian diagnostics. At large length/time scales, Lagrangian diffusivity is comparable for filtered and unfiltered velocities, while at short scales, unfiltered velocities disperse particles much faster. We also calculate diagnostics of Lagrangian coherent structures:rotationally coherent Lagrangian vortices detected from closed contours of the Lagrangian‐averaged vorticity deviation and material transport barriers formed by ridges of maximum finite‐time Lyapunov exponent. For temporally filtered velocities, we observe strong material coherence, which breaks down when the level of temporal filtering is reduced/removed, due to vigorous small‐scale mixing. In addition, for the lowest temporal resolution, the 3‐D particles experience very little vertical motion, suggesting that degrading temporal resolution greatly reduces vertical advection by high‐frequency motions. Our study suggests that Lagrangian diagnostics based on satellite‐derived surface geostrophic velocity fields, even with higher spatial resolutions as in the upcoming Surface Water and Ocean Topography mission, may overestimate the presence of mesoscale coherent structures and underestimate dispersion.

With the advent of high-performance computing, we are now capable of simulating the ocean and climate system on decadal to centennial timescales. However, global and basin-scale simulations still lack the spatial resolution necessary to resolve the mesoscales (hereon referred to as mesoscale-permitting simulations), a scale roughly on the order of O(100 km). Here, we provide a first step towards a potential vorticity (PV) based mesoscale closure scheme in order to improve the representation of mesoscale eddies in such simulations by taking advantage of the thickness-weighted averaged (TWA) framework. In the TWA framework the total eddy feedback can be encapsulated in the Eliassen-Palm (E-P) flux divergence. This implies that mesoscale closure schemes aimed at representing the total eddy feedback should therefore be representing the E-P flux divergence. The TWA framework further elucidates that its divergence is equivalent to the eddy Ertel PV flux. In other words, if one is to parametrize the eddy Ertel PV flux, one has parametrized the total eddy feedback onto the mean flow. Using a 1/12° North Atlantic ensemble simulation with 24 members, which allows us to decompose the mesoscale variability from the forced dynamics, we show that the eddy Ertel PV flux can be related to the local-gradient of mean Ertel PV as an active tracer via an anisotropic eddy diffusivity tensor. What follows is that not only does the tensor bring together the isopycnal thickness skew diffusivity and isopycnic tracer diffusivity, the former known as the Gent-McWilliams (GM) parametrization and latter the Redi parametrization, but also incorporates the eddy momentum fluxes. Although the Redi parametrization has existed longer than GM, there has been much more development in the latter, leaving the Redi diffusivity poorly constrained. Being able to treat GM and Redi simultaneously is another strength of our framework.

The thickness-weighted average (TWA) framework, which treats the residual-mean flow as the prognostic variable, has provided us with a clear theoretical understanding of the eddy feedback onto the residual-mean flow. The averaging operator involved in the TWA framework, although in theory being an ensemble mean, in practice has often been approximated by a temporal mean, which conflates the temporal variability with the eddies. Here, we analyze an ensemble of North Atlantic simulations at mesoscale resolving resolution (1/12{degree sign}). We therefore recognize means and eddies in terms of ensemble means and fluctuations about those means, in keeping with the TWA formalism proposed by Young (2012). Eddy-mean flow feedbacks are encapsulated in the Eliassen-Palm (E-P) flux tensor and its divergence indicates that the eddies contribute to the zonal meandering of the Gulf Stream and its deceleration in the meridional direction. We also show that the eddy Ertel potential vorticity (PV) flux can be parametrized as an isopycnic local-gradient flux of the residual-mean Ertel PV via an anisotropic eddy diffusivity tensor. As the E-P flux divergence and eddy Ertel PV flux are directly related to one another, this provides a new pathway forward for a unified mesoscale eddy closure scheme.

It has been hypothesized that submesoscale flows play an important role in the vertical transport of climatically important tracers, due to their strong associated vertical velocities. However, the multi-scale, non-linear, and Lagrangian nature of transport makes it challenging to attribute proportions of the tracer fluxes to certain processes, scales, regions, or features. Here we show that criteria based on the surface vorticity and strain joint probability distribution function (JPDF) effectively decomposes the surface velocity field into distinguishable flow regions, and different flow features, like fronts or eddies, are contained in different flow regions. The JPDF has a distinct shape and approximately parses the flow into different scales, as stronger velocity gradients are usually associated with smaller scales. Conditioning the vertical tracer transport on the vorticity-strain JPDF can therefore help to attribute the transport to different types of flows and scales. Applied to a set of idealized Antarctic Circumpolar Current simulations that vary only in horizontal resolution, this diagnostic approach demonstrates that small-scale strain dominated regions that are generally associated with submesoscale fronts, despite their minuscule spatial footprint, play an outsized role in exchanging tracers across the mixed layer base and are an important contributor to the large-scale tracer budgets. Resolving these flows not only adds extra flux at the small scales, but also enhances the flux due to the larger-scale flows.

Mixing by mesoscale eddies in the ocean plays a major role in setting the distribution of oceanic tracers, with important implications for physical and biochemical systems at local to global scales. Roach et al. (2016; https://doi.org/10.1002/2015JC011440) demonstrated that a two‐particle analysis of Argo trajectories produces robust estimates of horizontal mixing in the Southern Ocean. Here we extend this analysis to produce global 1° × 1° maps of eddy diffusivity at the nominal Argo parking depth of 1,000 m. We also applied this methodology to estimate surface eddy diffusivities from Global Drifter Program (GDP) surface drifters. The global mean eddy diffusivity was 543 ± 155 m2/s at 1,000 m and 2637 ± 311 m2/s at the surface, with elevated diffusivities in regions of enhanced eddy kinetic energy, such as western boundary currents and along the Antarctic Circumpolar Current. The eddy kinetic energy at the equator is high at both the surface and depth, but the eddy diffusivity is only enhanced near the surface. At depth the eddy diffusivity is strongly suppressed due to the presence of mean flow. We used our observational estimates to test the validity of an eddy diffusivity parameterization that accounts for mixing suppression in the presence of zonal mean flows. Our results indicated that this parameterization generally agrees with the directly observed eddy diffusivities in the midlatitude and high‐latitude oceans.

We provide the first observational estimate of the circumpolar distribution of cross‐stream eddy diffusivity at 1000 m in the Southern Ocean using Argo float trajectories. We show that Argo float trajectories, from the float surfacing positions, can be used to estimate lateral eddy diffusivities in the ocean and that these estimates are comparable to those obtained from RAFOS floats, where they overlap. Using the Southern Ocean State Estimate (SOSE) velocity fields to advect synthetic particles with imposed behavior that is "Argo‐like" and "RAFOS‐like" diffusivity estimates from both sets of synthetic particles agreed closely at the three dynamically very different test sites, the Kerguelen Island region, the Southeast Pacific Ocean, and the Scotia Sea, and support our approach. Observed cross‐stream diffusivities at 1000 m, calculated from Argo float trajectories, ranged between 300 and 2500 m2 s−1, with peaks corresponding to topographic features associated with the Scotia Sea, the Kerguelen Plateau, the Campbell Plateau, and the Southeast Pacific Ridge. These observational estimates agree with previous regional estimates from the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES) near the Drake Passage, and other estimates from natural tracers (helium), inverse modeling studies, and current meter measurements. These estimates are also compared to the suppressed eddy diffusivity in the presence of mean flows. The comparison suggests that away from regions of strong topographic steering suppression explains both the structure and magnitude of eddy diffusivity but that eddy diffusivities in the regions of topographic steering are greater than what would be theoretically expected and the ACC experiences localized enhanced cross‐stream mixing in these regions.

The large-scale middepth circulation and eddy diffusivities in the southeast Pacific Ocean and Scotia Sea sectors between 110° and 45°W of the Antarctic Circumpolar Current (ACC) are described based on a subsurface quasi-isobaric RAFOS-float-based Lagrangian dataset. These RAFOS float data were collected during the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). The mean flow, adjusted to a common 1400-m depth, shows the presence of jets in the time-averaged sense with speeds of 6 cm s−1 in the southeast Pacific Ocean and upward of 13 cm s−1 in the Scotia Sea. These jets appear to be locked to topography in the Scotia Sea but, aside from negotiating a seamount chain, are mostly free of local topographic constraints in the southeast Pacific Ocean. The eddy kinetic energy (EKE) is higher than the mean kinetic energy everywhere in the sampled domain by about 50%. The magnitude of the EKE increases drastically (by a factor of 2 or more) as the current crosses over the Hero and Shackleton fracture zones into the Scotia Sea. The meridional isopycnal stirring shows lateral and vertical variations with local eddy diffusivities as high as 2800 ± 600 m2 s−1 at 700 m decreasing to 990 ± 200 m2 s−1 at 1800 m in the southeast Pacific Ocean. However, the cross-ACC diffusivity in the southeast Pacific Ocean is significantly lower, with values of 690 ± 150 and 1000 ± 200 m2 s−1 at shallow and deep levels, respectively, due to the action of jets. The cross-ACC diffusivity in the Scotia Sea is about 1200 ± 500 m2 s−1.

As part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES), 210 subsurface floats were deployed west of the Drake Passage on two targeted density surfaces. Absolute (single particle) diffusivities are calculated for the floats. The focus is on the meridional component, which is less affected by the mean shear. The diffusivities are estimated in several ways, including a novel method based on the probability density function of the meridional displacements. This allows the determination of the range of possible lateral diffusivities, as well as the period over which the spreading can be said to be diffusive. The method is applied to the float data and to synthetic trajectories generated with the Massachusetts Institute of Technology General Circulation Model (MITgcm). Because of ballasting problems, many of the floats did not remain on their targeted density surface. However, the float temperature records suggest that most occupied a small range of densities, so the floats were grouped together for the analysis. The latter focuses on a subset of 109 of the floats, launched near 105°W. The different methods yield a consistent estimate for the diffusivity of 800 ± 200 m2 s−1. The same calculations were made with model particles deployed on 20 different density surfaces and the result for the particles deployed on the neutral density surface γ = 27.7 surface was the same within the errors. The model was then used to map the variation of the diffusivity in the vertical, near the core of the Antarctic Circumpolar Current (ACC). The results suggest mixing is intensified at middepths, between 1500 and 2000 m, consistent with several previous studies.