Isopycnal mixing of tracers is important for ocean dynamics and biogeochemistry. Previous studies have primarily focused on the horizontal structure of mixing, but what controls its vertical structure is still unclear. This study investigates the vertical structure of the isopycnal tracer diffusivity diagnosed by a multiple‐tracer inversion method in an idealized basin circulation model. The first two eigenvalues of the symmetric part of the 3D diffusivity tensor are approximately tangent to isopycnal surfaces. The isopycnal mixing is anisotropic, with principal directions of the large and small diffusivities generally oriented along and across the mean flow direction. The cross‐stream diffusivity can be reconstructed from the along‐stream diffusivity after accounting for suppression of mixing by the mean flow. In the circumpolar channel and the upper ocean in the gyres, the vertical structure of the along‐stream diffusivity follows that of the rms eddy velocity times a depth‐independent local energy‐containing scale estimated from the sea surface height. The diffusivity in the deep ocean in the gyres instead follows the profile of the eddy kinetic energy times a depth‐independent mixing time scale. The transition between the two mixing regimes is attributed to the dominance of nonlinear interactions and linear waves in the upper and deep ocean, respectively, distinguished by a nonlinearity parameter. A formula is proposed that accounts for both regimes and captures the vertical variation of diffusivities better than extant theories. These results inform efforts to parameterize the vertical structure of isopycnal mixing in coarse‐resolution ocean models.
- Award ID(s):
- 1829969
- PAR ID:
- 10298804
- Date Published:
- Journal Name:
- Journal of Physical Oceanography
- Volume:
- 51
- Issue:
- 3
- ISSN:
- 0022-3670
- Page Range / eLocation ID:
- 845 to 860
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Mixing along isopycnals plays an important role in the transport and uptake of oceanic tracers. Isopycnal mixing is commonly quantified by a tracer diffusivity. Previous studies have estimated the tracer diffusivity using the rate of dispersion of surface drifters, subsurface floats, or numerical particles advected by satellite‐derived velocity fields. This study shows that the diffusivity can be more efficiently estimated from the dispersion of coherent mesoscale eddies. Coherent eddies are identified and tracked as the persistent sea surface height extrema in both a two‐layer quasigeostrophic (QG) model and an idealized primitive equation (PE) model. The Lagrangian diffusivity is estimated using the tracks of these coherent eddies and compared to the diagnosed Eulerian diffusivity. It is found that the meridional coherent eddy diffusivity approaches a stable value within about 20–40 days in both models. In the QG model, the coherent eddy diffusivity is a good approximation to the upper‐layer tracer diffusivity in a broad range of flow regimes, except for small values of bottom friction or planetary vorticity gradient, where the motions of same‐sign eddies are correlated over long distances. In the PE model, the tracer diffusivity has a complicated vertical structure and the coherent eddy diffusivity is correlated with the tracer diffusivity at the e‐folding depth of the energy‐containing eddies where the intrinsic speed of the coherent eddies matches the rms eddy velocity. These results suggest that the oceanic tracer diffusivity at depth can be estimated from the movements of coherent mesoscale eddies, which are routinely tracked from satellite observations.
-
Abstract The middepth ocean temperature profile was found by Munk in 1966 to agree with an exponential profile and shown to be consistent with a vertical advective–diffusive balance. However, tracer release experiments show that vertical diffusivity in the middepth ocean is an order of magnitude too small to explain the observed 1-km exponential scale. Alternative mechanisms suggested that nearly all middepth water upwells adiabatically in the Southern Ocean (SO). In this picture, SO eddies and wind set SO isopycnal slopes and therefore determine a nonvanishing middepth interior stratification even in the adiabatic limit. The effect of SO eddies on SO isopycnal slopes can be understood via either a marginal criticality condition or a near-vanishing SO residual deep overturning condition in the adiabatic limit. We examine the interplay between SO dynamics and interior mixing in setting the exponential profiles of σ 2 and ∂ z σ 2 . We use eddy-permitting numerical simulations, in which we artificially change the diapycnal mixing only away from the SO. We find that SO isopycnal slopes change in response to changes in the interior diapycnal mixing even when the wind forcing is constant, consistent with previous studies (that did not address these near-exponential profiles). However, in the limit of small interior mixing, the interior ∂ z σ 2 profile is not exponential, suggesting that SO processes alone, in an adiabatic limit, do not lead to the observed near-exponential structures of such profiles. The results suggest that while SO wind and eddies contribute to the nonvanishing middepth interior stratification, the exponential shape of the ∂ z σ 2 profiles must also involve interior diapycnal mixing.more » « less
-
Due to their limited resolution, numerical ocean models need to be interpreted as representing filtered or averaged equations. How to interpret models in terms of formally averaged equations, however, is not always clear, particularly in the case of hybrid or generalized vertical coordinate models, which limits our ability to interpret the model results and to develop parameterizations for the unresolved eddy contributions. We here derive the averaged hydrostatic Boussinesq equations in generalized vertical coordinates for an arbitrary thickness‐weighted average. We then consider various special cases and discuss the extent to which the averaged equations are consistent with existing ocean model formulations. As previously discussed, the momentum equations in existing depth‐coordinate models are best interpreted as representing Eulerian averages (i.e., averages taken at fixed depth), while the tracer equations can be interpreted as either Eulerian or thickness‐weighted isopycnal averages. Instead we find that no averaging is fully consistent with existing formulations of the parameterizations in semi‐Lagrangian discretizations of generalized vertical coordinate ocean models such as MOM6. A coordinate‐following average would require “coordinate‐aware” parameterizations that can account for the changing nature of the eddy terms as the coordinate changes. Alternatively, the model variables can be interpreted as representing either Eulerian or (thickness‐weighted) isopycnal averages, independent of the model coordinate that is being used for the numerical discretization. Existing parameterizations in generalized vertical coordinate models, however, are not always consistent with either of these interpretations, which, respectively, would require a three‐dimensional divergence‐free eddy tracer advection or a form‐stress parameterization in the momentum equations.more » « less
-
more » « less
Abstract Paleoproxy observations suggest that deep-ocean water-mass distributions were different at the Last Glacial Maximum than they are today. However, even modern deep-ocean water-mass distributions are not completely explained by observations of the modern ocean circulation. This paper investigates two processes that influence deep-ocean water-mass distributions: 1) interior downwelling caused by vertical mixing that increases in the downward direction and 2) isopycnal mixing. Passive tracers are used to assess how changes in the circulation and in the isopycnal-mixing coefficient impact deep-ocean water-mass distributions in an idealized two-basin model. We compare two circulations, one in which the upper cell of the overturning reaches to 4000-m depth and one in which it shoals to 2500-m depth. Previous work suggests that in the latter case the upper cell and the abyssal cell of the overturning are separate structures. Nonetheless, high concentrations of North Atlantic Water (NAW) are found in our model’s abyssal cell: these tracers are advected into the abyssal cell by interior downwelling caused by our vertical mixing profile, which increases in the downward direction. Further experiments suggest that the NAW concentration in the deep South Atlantic Ocean and in the deep Pacific Ocean is influenced by the isopycnal-mixing coefficient in the top 2000 m of the Southern Ocean. Both the strength and the vertical profile of isopycnal mixing are important for setting deep-ocean tracer concentrations. A 1D advection–diffusion model elucidates how NAW concentration depends on advective and diffusive processes.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JPO-D-20-0077.1
