Baroclinic Turbulence in the Ocean: Analysis with Primitive Equation and Quasigeostrophic Simulations

2011 ◽  
Vol 41 (9) ◽  
pp. 1605-1623 ◽  
Author(s):  
Antoine Venaille ◽  
Geoffrey K. Vallis ◽  
K. Shafer Smith

Abstract This paper examines the factors determining the distribution, length scale, magnitude, and structure of mesoscale oceanic eddies in an eddy-resolving primitive equation simulation of the Southern Ocean [Modeling Eddies in the Southern Ocean (MESO)]. In particular, the authors investigate the hypothesis that the primary source of mesoscale eddies is baroclinic instability acting locally on the mean state. Using local mean vertical profiles of shear and stratification from an eddying primitive equation simulation, the forced–dissipated quasigeostrophic equations are integrated in a doubly periodic domain at various locations. The scales, energy levels, and structure of the eddies found in the MESO simulation are compared to those predicted by linear stability analysis, as well as to the eddying structure of the quasigeostrophic simulations. This allows the authors to quantitatively estimate the role of local nonlinear effects and cascade phenomena in the generation of the eddy field. There is a modest transfer of energy (an “inverse cascade”) to larger scales in the horizontal, with the length scale of the resulting eddies typically comparable to or somewhat larger than the wavelength of the most unstable mode. The eddies are, however, manifestly nonlinear, and in many locations the turbulence is fairly well developed. Coherent structures also ubiquitously emerge during the nonlinear evolution of the eddy field. There is a near-universal tendency toward the production of grave vertical scales, with the barotropic and first baroclinic modes dominating almost everywhere, but there is a degree of surface intensification that is not captured by these modes. Although the results from the local quasigeostrophic model compare well with those of the primitive equation model in many locations, some profiles do not equilibrate in the quasigeostrophic model. In many cases, bottom friction plays an important quantitative role in determining the final scale and magnitude of eddies in the quasigeostrophic simulations.

2009 ◽  
Vol 66 (4) ◽  
pp. 837-863 ◽  
Author(s):  
Pablo Zurita-Gotor ◽  
Geoffrey K. Vallis

Abstract This paper investigates the equilibration of baroclinic turbulence in an idealized, primitive equation, two-level model, focusing on the relation with the phenomenology of quasigeostrophic turbulence theory. Simulations with a comparable two-layer quasigeostrophic model are presented for comparison, with the deformation radius in the quasigeostrophic model being set using the stratification from the primitive equation model. Over a fairly broad parameter range, the primitive equation and quasigeostrophic results are in qualitative and, to some degree, quantitative agreement and are consistent with the phenomenology of geostrophic turbulence. The scale, amplitude, and baroclinicity of the eddies and the degree of baroclinic instability of the mean flow all vary fairly smoothly with the imposed parameters; both models are able, in some parameter ranges, to produce supercritical flows. The criticality in the primitive equation model, which does not have any convective parameterization scheme, is fairly sensitive to the external parameters, most notably the planet size (i.e., the f /β ratio), the forcing time scale, and the factors influencing the stratification. In some parameter settings of the models, although not those that are most realistic for the earth’s atmosphere, it is possible to produce eddies that are considerably larger than the deformation scales and an inverse cascade in the barotropic flow with a −5/3 spectrum. The vertical flux of heat is found to be related to the isentropic slope.


2010 ◽  
Vol 40 (2) ◽  
pp. 257-278 ◽  
Author(s):  
Andrew F. Thompson

Abstract Satellite altimetry and high-resolution ocean models indicate that the Southern Ocean comprises an intricate web of narrow, meandering jets that undergo spontaneous formation, merger, and splitting events, as well as rapid latitude shifts over periods of weeks to months. The role of topography in controlling jet variability is explored using over 100 simulations from a doubly periodic, forced-dissipative, two-layer quasigeostrophic model. The system is forced by a baroclinically unstable, vertically sheared mean flow in a domain that is large enough to accommodate multiple jets. The dependence of (i) meridional jet spacing, (ii) jet variability, and (iii) domain-averaged meridional transport on changes in the length scale and steepness of simple sinusoidal topographical features is analyzed. The Rhines scale, ℓβ = 2πVe/β, where Ve is an eddy velocity scale and β is the barotropic potential vorticity gradient, measures the meridional extent of eddy mixing by a single jet. The ratio ℓβ /ℓT, where ℓT is the topographic length scale, governs jet behavior. Multiple, steady jets with fixed meridional spacing are observed when ℓβ ≫ ℓT or when ℓβ ≈ ℓT. When ℓβ < ℓT, a pattern of perpetual jet formation and jet merger dominates the time evolution of the system. Zonal ridges systematically reduce the domain-averaged meridional transport, while two-dimensional, sinusoidal bumps can increase transport by an order of magnitude or more. For certain parameters, bumpy topography gives rise to periodic oscillations in the jet structure between purely zonal and topographically steered states. In these cases, transport is dominated by bursts of mixing associated with the transition between the two regimes. Topography modifies local potential vorticity (PV) gradients and mean flows; this can generate asymmetric Reynolds stresses about the jet core and can feed back on the conversion of potential energy to kinetic energy through baroclinic instability. Both processes contribute to unsteady jet behavior. It is likely that these processes play a role in the dynamic nature of Southern Ocean jets.


2003 ◽  
Vol 33 (8) ◽  
pp. 1579-1591 ◽  
Author(s):  
S. S. Drijfhout

Abstract The question of whether anticyclones can split and break up is readdressed using a numerical, multilayer, primitive equation model. Applying the conservation of integrated angular momentum (IAM) to barotropic and baroclinic vortices, it has been argued that anticyclones can never split, no matter what their structure is. When an anticyclone splits, the IAM has to increase as the newly formed eddies are pushed away from their original center. Conservation of IAM prohibits such an increase. Several numerical simulations, however, have shown anticyclonic splitting. In a multilayer model, a vertical transport of IAM is possible. For counterrotating eddies (an anticyclone on top of a cyclone) it is easy to see that a vertical exchange of IAM allows the eddy to break up. For a compensated or weakly corotating eddy, breakup is only possible when, in addition to a vertical transport of IAM, in the deep layer(s) IAM is exchanged between the core of the vortex and the surrounding fluid. In the presence of a tilting interface, the pressure gradient associated with the sea surface height (SSH) anomaly, in particular its non-equivalent-barotropic part, drives the required exchanges. The non-equivalent-barotropic SSH anomaly is associated with the vertical phase lag of the most unstable eigenmode (m = 2), which develops when this mode gains energy by baroclinic energy conversion. The previous conclusion that anticyclones cannot split on their own should be revised to the following: anticyclones cannot split by barotropic processes alone—baroclinic instability is a necessary ingredient for splitting to occur.


2018 ◽  
Vol 48 (4) ◽  
pp. 937-957 ◽  
Author(s):  
Benjamin A. Storer ◽  
Francis J. Poulin ◽  
Claire Ménesguen

AbstractThe stability of lens-shaped vortices is revisited in the context of an idealized quasigeostrophic model. We compute the stability characteristics with higher accuracy and for a wider range of Burger numbers (Bu) than what was previously done. It is found that there are four distinct Bu regions of linear instability. Over the primary region of interest (0.1 < Bu < 10), we confirm that the first and second azimuthal modes are the only linearly unstable modes, and they are associated with vortex tilting and tearing, respectively. Moreover, the most unstable first azimuthal mode is not precisely captured by the linear stability analysis because of the extra condition that is imposed at the vortex center, and accurate calculations of the second azimuthal mode require higher resolution than was previously considered. We also study the nonlinear evolution of lens-shaped vortices in the context of this model and present the following results. First, vortices with a horizontal length scale a little less than the radius of deformation (Bu > 1) are barotropically unstable and develop a wobble, whereas those with a larger horizontal length scale (Bu < 1) are baroclinically unstable and often split. Second, the transfer of energy between different horizontal scales is quantified in two typical cases of barotropic and baroclinic instability. Third, after the instability the effective Bu is closer to unity.


2011 ◽  
Vol 68 (4) ◽  
pp. 823-838 ◽  
Author(s):  
Pablo Zurita-Gotor ◽  
Geoffrey K. Vallis

Abstract This paper investigates the factors that determine the equilibrium state, and in particular the height and structure of the tropopause, in an idealized primitive equation model forced by Newtonian cooling in which the eddies can determine their own depth. Previous work has suggested that the midlatitude tropopause height may be understood as the intersection between a radiative and a dynamical constraint. The dynamical constraint relates to the lateral transfer of energy, which in midlatitudes is largely effected by baroclinic eddies, and its representation in terms of mean-flow properties. Various theories have been proposed and investigated for the representation of the eddy transport in terms of the mean flow, including a number of diffusive closures and the notion that the flow evolves to a state marginally supercritical to baroclinic instability. The radiative constraint expresses conservation of energy and so must be satisfied, although it need not necessarily be useful in providing a tight constraint on tropopause height. This paper explores whether and how the marginal criticality and radiative constraints work together to produce an equilibrated flow and a tropopause that is internal to the fluid. The paper investigates whether these two constraints are consistent with simulated variations in the tropopause height and in the mean state when the external parameters of an idealized primitive equation model are changed. It is found that when the vertical redistribution of heat is important, the radiative constraint tightly constrains the tropopause height and prevents an adjustment to marginal criticality. In contrast, when the stratification adjustment is small, the radiative constraint is only loosely satisfied and there is a tendency for the flow to adjust to marginal criticality. In those cases an alternative dynamical constraint would be needed in order to close the problem and determine the eddy transport and tropopause height in terms of forcing and mean flow.


2006 ◽  
Vol 63 (5) ◽  
pp. 1548-1557 ◽  
Author(s):  
G. P. Williams

Abstract Altering the tropospheric static stability changes the nature of the equatorial superrotation associated with unstable, low-latitude, westerly jets, according to calculations with a dry, global, multilevel, spectral, primitive equation model subject to a simple Newtonian heating function. For a low static stability, the superrotation fluxes with the simplest structure occur when the stratospheric extent and horizontal diffusion are minimal. Barotropic instability occurs on the jet's equatorward flank and baroclinic instability occurs on the jet's poleward flank. Systems with a high static stability inhibit the baroclinic instability and thereby reveal more clearly that the barotropic instability is the primary process driving the equatorial superrotation. Such systems produce a flatter equatorial jet and also take much longer to equilibrate than the standard atmospheric circulation.


2012 ◽  
Vol 42 (1) ◽  
pp. 18-38 ◽  
Author(s):  
G. Roullet ◽  
J. C. McWilliams ◽  
X. Capet ◽  
M. J. Molemaker

Abstract High-resolution simulations of β-channel, zonal-jet, baroclinic turbulence with a three-dimensional quasigeostrophic (QG) model including surface potential vorticity (PV) are analyzed with emphasis on the competing role of interior and surface PV (associated with isopycnal outcropping). Two distinct regimes are considered: a Phillips case, where the PV gradient changes sign twice in the interior, and a Charney case, where the PV gradient changes sign in the interior and at the surface. The Phillips case is typical of the simplified turbulence test beds that have been widely used to investigate the effect of ocean eddies on ocean tracer distribution and fluxes. The Charney case shares many similarities with recent high-resolution primitive equation simulations. The main difference between the two regimes is indeed an energization of submesoscale turbulence near the surface. The energy cycle is analyzed in the (k, z) plane, where k is the horizontal wavenumber. In the two regimes, the large-scale buoyancy forcing is the primary source of mechanical energy. It sustains an energy cycle in which baroclinic instability converts more available potential energy (APE) to kinetic energy (KE) than the APE directly injected by the forcing. This is due to a conversion of KE to APE at the scale of arrest. All the KE is dissipated at the bottom at large scales, in the limit of infinite resolution and despite the submesoscales energizing in the Charney case. The eddy PV flux is largest at the scale of arrest in both cases. The eddy diffusivity is very smooth but highly nonuniform. The eddy-induced circulation acts to flatten the mean isopycnals in both cases.


Sign in / Sign up

Export Citation Format

Share Document