ALBERT

Description: Ocean ventilation is the process by which climatically important tracers such as heat and carbon are exchanged between the atmosphere and ocean interior. In this paper a series of numerical simulations are used to study the interaction of submesoscales and a topographically steered jet in driving rapid ventilation. The ventilation is found to increase both as a function of wind stress and model resolution, with a submesoscale-resolving 1/120° model exhibiting the largest ventilation rate. The jet in this simulation is found to be persistently unstable to submesoscale instabilities, which are known to feature intense vertical circulations. The vertical tracer transport is found to scale as a function of the eddy kinetic energy and mean isopycnal slope, whose behaviors change as a function of the wind stress and due to the emergence of a strong potential vorticity gradient due to the lateral shear of the jet. These results highlight the importance of jet–submesoscale interaction as a bridge between the atmosphere and the ocean interior.

Description: Despite slow rates of ocean mixing, observational and modeling studies suggest that buoyancy is redistributed to all depths of the ocean on surprisingly short interannual to decadal time scales. The mechanisms responsible for this redistribution remain poorly understood. This work uses an Earth system model to evaluate the global steady-state ocean buoyancy (and related steric sea level) budget, its interannual variability, and its transient response to a doubling of CO2 over 70 years, with a focus on the deep ocean. At steady state, the simple view of vertical advective–diffusive balance for the deep ocean holds at low to midlatitudes. At higher latitudes, the balance depends on a myriad of additional terms, namely mesoscale and submesoscale advection, convection and overflows from marginal seas, and terms related to the nonlinear equation of state. These high-latitude processes rapidly communicate anomalies in surface buoyancy forcing to the deep ocean locally; the deep, high-latitude changes then influence the large-scale advection of buoyancy to create transient deep buoyancy anomalies at lower latitudes. Following a doubling of atmospheric carbon dioxide concentrations, the high-latitude buoyancy sinks are suppressed by a slowdown in convection and reduced dense water formation. This change is accompanied by a slowing of both upper and lower cells of the global meridional overturning circulation, reducing the supply of dense water to low latitudes beneath the pycnocline and the commensurate flow of light waters to high latitudes above the pycnocline. By this mechanism, changes in high-latitude buoyancy are communicated to the global deep ocean on relatively fast advective time scales.

Description: Recent work on the global overturning circulation and its energetics assumes that processes caused by nonlinearities of the equation of state of seawater are negligible. Nonlinear processes such as cabbeling and thermobaricity cause diapycnal motion as a consequence of isopycnal mixing. The nonlinear equation of state also causes the helical nature of neutral trajectories; as a consequence of this helical nature, it is not possible to define a continuous “density” surface that aligns with neutral tangent planes. The result is an additional diapycnal advection, which needs to be accounted for in water-mass analysis. In this paper, the authors take advantage of new techniques in constructing very accurate continuous density surfaces to more precisely estimate isopycnal and diapycnal processes caused by the nonlinear equation of state. They then quantify the diapycnal advection due to each of these nonlinear processes and show that they lead in total to a significant downward diapycnal advection, particularly in the Southern Ocean. The nonlinear processes are therefore another source of dense water formation in addition to high-latitude convection. To maintain the abyssal stratification in the global ocean, these dense water masses have to be brought back toward surface layers, and this can occur by either diabatic or adiabatic processes. Including these nonlinear processes into the advection–diffusion balance, the authors show that observed diapycnal diffusivities are unlikely to be able to support the amount of dense water produced in the global ocean, thus placing more importance on the adiabatic way of bringing the deep waters back to the surface.

Description: In the absence of diapycnal mixing processes, fluid parcels move in directions along which they do not encounter buoyant forces. These directions define the local neutral tangent plane. Because of the nonlinear nature of the equation of state of seawater, these neutral tangent planes cannot be connected globally to form a well-defined surface in three-dimensional space; that is, continuous “neutral surfaces” do not exist. This inability to form well-defined neutral surfaces implies that neutral trajectories are helical. Consequently, even in the absence of diapycnal mixing processes, fluid trajectories penetrate through any “density” surface. This process amounts to an extra mechanism that achieves mean vertical advection through any continuous surface such as surfaces of constant potential density or neutral density. That is, the helical nature of neutral trajectories causes this additional diasurface velocity. A water-mass analysis performed with respect to continuous density surfaces will have part of its diapycnal advection due to this diasurface advection process. Hence, this additional diasurface advection should be accounted for when attributing observed water-mass changes to mixing processes. Here, the authors quantify this component of the total diasurface velocity and show that locally it can be the same order of magnitude as diasurface velocities produced by other mixing processes, particularly in the Southern Ocean. The magnitude of this diasurface advection is proportional to the ocean’s neutral helicity, which is observed to be quite small in today’s ocean. The authors also use a perturbation experiment to show that the ocean rapidly readjusts to its present state of small neutral helicity, even if perturbed significantly. Additionally, the authors show how seasonal (rather than spatial) changes in the ocean’s hydrography can generate a similar vertical advection process. This process is described here for the first time; although the vertical advection due to this process is small, it helps to understand water-mass transformation on density surfaces.

Description: We introduce a simple algorithm to improve existing density surfaces to ensure that the resulting surfaces are as close to neutral as possible. This means the slopes at any point on the surfaces are close to neutral tangent planes – the directions along which layered stirring and mixing occurs – minimizing the fictitious diapycnal diffusivity. Inverse techniques and layered models have been used for decades to understand ocean circulation. The most-used density surfaces are potential density or neutral density surfaces. Both these density surfaces and all others produce a fictitious diapycnal diffusivity to some degree due to the helical nature of neutral trajectories – with the magnitude of this artificial diffusivity in some cases being larger than the values measured in the ocean. Here we show how this error can be reduced by up to four orders of magnitude and therefore becomes insignificant compared to measured values, thus providing surfaces which would produce more accurate results when used for inverse techniques.

Description: We introduce a simple algorithm to improve existing density surfaces to ensure that the resulting surfaces are as close to neutral as possible. This means the slopes at any point on the surfaces are close to neutral tangent planes – the directions along which layered stirring and mixing occurs – minimizing the fictitious diapycnal diffusivity. Inverse techniques and layered models have been used for decades to understand ocean circulation. The most-used density surfaces are potential density or neutral density surfaces. Both these density surfaces and all others produce a fictitious diapycnal diffusivity to some degree due to the helical nature of neutral trajectories – with the magnitude of this artificial diffusivity in some cases being larger than the values measured in the ocean. Here we show how this error can be reduced by up to four orders of magnitude and therefore becomes insignificant compared to measured values, thus providing surfaces which would produce more accurate results when used for inverse techniques.

Description: The first direct estimate of the rate at which geostrophic turbulence mixes tracers across the Antarctic Circumpolar Current is presented. The estimate is computed from the spreading of a tracer released upstream of Drake Passage as part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). The meridional eddy diffusivity, a measure of the rate at which the area of the tracer spreads along an isopycnal across the Antarctic Circumpolar Current, is 710 ± 260 m2 s−1 at 1500-m depth. The estimate is based on an extrapolation of the tracer-based diffusivity using output from numerical tracers released in a one-twentieth of a degree model simulation of the circulation and turbulence in the Drake Passage region. The model is shown to reproduce the observed spreading rate of the DIMES tracer and suggests that the meridional eddy diffusivity is weak in the upper kilometer of the water column with values below 500 m2 s−1 and peaks at the steering level, near 2 km, where the eddy phase speed is equal to the mean flow speed. These vertical variations are not captured by ocean models presently used for climate studies, but they significantly affect the ventilation of different water masses.

Description: On planetary scales, surface wind stress and differential buoyancy forcing act together to produce isopycnal surfaces that are relatively flat in the tropics/subtropics and steep near the poles, where they tend to outcrop. Tilted isopycnals in a rapidly rotating fluid are subject to baroclinic instability. The turbulent, mesoscale eddies generated by this instability have a tendency to homogenize potential vorticity (PV) along density surfaces. In the Southern Ocean (SO), the tilt of isopycnals is largely maintained by competition between the steepening effect of surface forcing and the flattening effect of turbulent, spatially inhomogeneous eddy fluxes of PV. Here quasigeostrophic theory is used to investigate the influence of a planetary–geometric constraint on the equilibrium slope of tilted density/buoyancy surfaces in the SO. If the meridional gradients of relative vorticity and PV are small relative to β, then quasigeostrophic theory predicts ds/dz = β/f0 = cot(ϕ0)/a, or equivalently r ≡ |∂zs/(β/f0)| = 1, where f is the Coriolis parameter, β is the meridional gradient of f, s is the isopycnal slope, ϕ0 is a reference latitude, a is the planetary radius, and r is the depth-averaged criticality parameter. It is found that the strict r = 1 condition holds over specific averaging volumes in a large-scale climatology. A weaker r = O(1) condition for depth-averaged quantities is generally satisfied away from large bathymetric features. The r = O(1) constraint is employed to derive a depth scale to characterize large-scale interior stratification, and an idealized sector model is used to test the sensitivity of this relationship to surface wind forcing. Finally, the possible implications for eddy flux parameterization and for the sensitivity of SO circulation/stratification to changes in forcing are discussed.

Description: Mesoscale eddies play a major role in the transport of tracers in the ocean. Focusing on a sector in the east Pacific, the authors present estimates of eddy diffusivities derived from kinematic tracer simulations using satellite-observed velocity fields. Meridional diffusivities are diagnosed, and how they are related to eddy properties through the mixing length formulation of Ferrari and Nikurashin, which accounts for the suppression of diffusivity due to eddy propagation relative to the mean flow, is shown. The uniqueness of this study is that, through systematically varying the zonal-mean flow, a hypothetical “unsuppressed” diffusivity is diagnosed. At a given latitude, the unsuppressed diffusivity occurs when the zonal-mean flow equals the eddy phase speed. This provides an independent estimate of eddy phase propagation, which agrees well with theoretical arguments. It is also shown that the unsuppressed diffusivity is predicted very well by classical mixing length theory, that is, that it is proportional to the rms eddy velocity times the observed eddy size, with a spatially constant mixing efficiency of 0.35. Then, the suppression factor is estimated and it is shown that it too can be understood quantitatively in terms of easily observed mean flow properties. The authors then extrapolate from these sector experiments to the global scale, making predictions for the global surface eddy diffusivity. Together with a prognostic equation for eddy kinetic energy and a theory explaining observed eddy sizes, these concepts could potentially be used in a closure for eddy diffusivities in coarse-resolution ocean climate models.