Quasi-geostrophic Flow Research Articles

Variational treatment of inertia–gravity waves interacting with a quasi-geostrophic mean flow

The equations for three-dimensional hydrostatic Boussinesq dynamics are equivalent to a variational principle that is closely analogous to the variational principle for classical electrodynamics. Inertia–gravity waves are analogous to electromagnetic waves, and available potential vorticity (i.e. the amount by which the potential vorticity exceeds the potential vorticity of the rest state) is analogous to electric charge. The Lagrangian can be expressed as the sum of three parts. The first part corresponds to quasi-geostrophic dynamics in the absence of inertia–gravity waves. The second part corresponds to inertia–gravity waves in the absence of quasi-geostrophic flow. The third part represents a coupling between the inertia–gravity waves and quasi-geostrophic motion. This formulation provides the basis for a general theory of inertia–gravity waves interacting with a quasi-geostrophic mean flow.

A Simple Model for Sheddies: Ocean Eddies Formed from Shed Vorticity

AbstractRecent studies show that vertical eddy diffusivity is sufficient on its own to introduce intense horizontal shear layers at sloping ocean margins (Molemaker et al.; Gula et al.; Dewar et al.). These layers influence mesoscale energy and potential vorticity budgets but cannot be fully represented in models without sloping boundaries, no-slip boundary conditions, and sufficiently high resolution. This paper investigates the detachment of these shear layers and their subsequent rolling up into concentrated eddies. These shed eddies, or “sheddies,” may have significant oceanographic impacts. Their growth is considered using a simple point vortex model that adapts the Brown–Michael model of vortex shedding to quasigeostrophic flow and allows detailed consideration of the vorticity fluxes. The model shows good qualitative agreement with observations and experimental and numerical results. It is applied to a number of examples of well-known cases of sheddy formation, including the Agulhas cyclones, California Undercurrent, and Canary Eddy Corridor, and also is used to investigate the effects of shed vorticity in the growth of the Cook Strait eddy and the interaction of the North Brazil Current rings with the islands of the Lesser Antilles.

A three-component model for the coupled evolution of near-inertial waves, quasi-geostrophic flow and the near-inertial second harmonic

We derive an asymptotic model that describes the nonlinear coupled evolution of (i) near-inertial waves (NIWs), (ii) balanced quasi-geostrophic flow and (iii) near-inertial second harmonic waves with frequency near $2f_{0}$, where $f_{0}$ is the local inertial frequency. This ‘three-component’ model extends the two-component model derived by Xie & Vanneste (J. Fluid Mech., vol. 774, 2015, pp. 143–169) to include interactions between near-inertial and $2f_{0}$ waves. Both models possess two conservation laws which together imply that oceanic NIWs forced by winds, tides or flow over bathymetry can extract energy from quasi-geostrophic flows. A second and separate implication of the three-component model is that quasi-geostrophic flow catalyses a loss of NIW energy to freely propagating waves with near-$2f_{0}$ frequency that propagate rapidly to depth and transfer energy back to the NIW field at very small vertical scales. The upshot of near-$2f_{0}$ generation is a two-step mechanism whereby quasi-geostrophic flow catalyses a nonlinear transfer of near-inertial energy to the small scales of wave breaking and diapycnal mixing. A comparison of numerical solutions with both Boussinesq and three-component models for a two-dimensional initial value problem reveals strengths and weaknesses of the model while demonstrating the extraction of quasi-geostrophic energy and production of small vertical scales.

Recent geomagnetic secular variation from Swarm and ground observatories as estimated in the CHAOS-6 geomagnetic field model

We use more than 2 years of magnetic data from the Swarm mission, and monthly means from 160 ground observatories as available in March 2016, to update the CHAOS time-dependent geomagnetic field model. The new model, CHAOS-6, provides information on time variations of the core-generated part of the Earth’s magnetic field between 1999.0 and 2016.5. We present details of the secular variation (SV) and secular acceleration (SA) from CHAOS-6 at Earth’s surface and downward continued to the core surface. At Earth’s surface, we find evidence for positive acceleration of the field intensity in 2015 over a broad area around longitude 90°E that is also seen at ground observatories such as Novosibirsk. At the core surface, we are able to map the SV up to at least degree 16. The radial field SA at the core surface in 2015 is found to be largest at low latitudes under the India–South-East Asia region, under the region of northern South America, and at high northern latitudes under Alaska and Siberia. Surprisingly, there is also evidence for significant SA in the central Pacific region, for example near Hawaii where radial field SA is observed on either side of a jerk in 2014. On the other hand, little SV or SA has occurred over the past 17 years in the southern polar region. Inverting for a quasi-geostrophic core flow that accounts for this SV, we obtain a prominent planetary-scale, anti-cyclonic, gyre centred on the Atlantic hemisphere. We also find oscillations of non-axisymmetric, azimuthal, jets at low latitudes, for example close to 40°W, that may be responsible for localized SA oscillations. In addition to scalar data from Orsted, CHAMP, SAC-C and Swarm, and vector data from Orsted, CHAMP and Swarm, CHAOS-6 benefits from the inclusion of along-track differences of scalar and vector field data from both CHAMP and the three Swarm satellites, as well as east–west differences between the lower pair of Swarm satellites, Alpha and Charlie. Moreover, ground observatory SV estimates are fit to a Huber-weighted rms level of 3.1 nT/year for the eastward components and 3.8 and 3.7 nT/year for the vertical and southward components. We also present an update of the CHAOS high-degree lithospheric field, making use of along-track differences of CHAMP scalar and vector field data to produce a new static field model that agrees well with the MF7 field model out to degree 110.

Balance dynamics in rotating stratified turbulence

If classical quasigeostrophic (QG) flow breaks down at smaller scales, it gives rise to questions of whether higher-order nonlinear balance can be maintained, to what scale and for how long. These are naturally followed by asking how this is affected by stratification and rotation. To address these questions, we perform non-hydrostatic Boussinesq simulations where the initial data is balanced using the Baer–Tribbia nonlinear normal mode initialization scheme (NNMI), which is accurate to second order in the Rossby number, as the next-order improvement to first-order QG theory. The NNMI procedure yields an ageostrophic contribution to the energy spectrum that has a very steep slope. However, as time passes, a shallow range emerges in the ageostrophic spectrum when the Rossby number is large enough for a given Reynolds number. It is argued that this shallow range is the unbalanced part of the motion that develops spontaneously in time and eventually dominates the energy at small scales. If the initial flow is not nonlinearly balanced, the shallow range emerges at even lower Rossby number and it appears at larger scales. Through numerous simulations at different rotation and stratification, this study gives a clear picture of how energy is cascaded in different initially balanced regimes of rotating stratified flow. We find that at low Rossby number the flow mainly consists of a geostrophic part and a balanced ageostrophic part with a steep spectrum. As the Rossby number increases, the unbalanced part of the ageostrophic energy increases at a rate faster than the balanced part. Hence, the total energy spectrum displays a shallow range above a transition wavenumber. This wavenumber evolves to smaller values as rotation weakens.

An ocean drum: quasi-geostrophic energetics from a Riemann geometry perspective

We revisit the discussion of the energetics of quasi-geostrophic flows from a geometric perspective based on the introduction of an effective metric, built in terms of the flow stratification and the Coriolis parameter. In particular, an appropriate notion of normal modes is defined through a spectral geometry problem in the ocean basin (a compact manifold with boundary) for the associated Laplace–Beltrami scalar operator. This spectral problem can be used to systematically encode non-local aspects of stratification and topography. As examples of applications we revisit the isotropy assumption in geostrophic turbulence, identify (a patch of) the hyperbolic space as the leading-order term in the effective geometry for the deep mesoscale ocean and, finally, discuss some diagnostic tools based on a simple statistical mechanics toy-model to be used in numerical simulations and/or observations of quasi-geostrophic flows.

Available potential vorticity and wave-averaged quasi-geostrophic flow

We derive a wave-averaged potential vorticity equation describing the evolution of strongly stratified, rapidly rotating quasi-geostrophic (QG) flow in a field of inertia-gravity internal waves. The derivation relies on a multiple-time-scale asymptotic expansion of the Eulerian Boussinesq equations. Our result confirms and extends the theory of Bühler & McIntyre (J. Fluid Mech., vol. 354, 1998, pp. 609–646) to non-uniform stratification with buoyancy frequency$N(z)$and therefore non-uniform background potential vorticity$f_{0}N^{2}(z)$, and does not require spatial-scale separation between waves and balanced flow. Our interest in non-uniform background potential vorticity motivates the introduction of a new quantity: ‘available potential vorticity’ (APV). Like Ertel potential vorticity, APV is exactly conserved on fluid particles. But unlike Ertel potential vorticity, linear internal waves have no signature in the Eulerian APV field, and the standard QG potential vorticity is a simple truncation of APV for low Rossby number. The definition of APV exactly eliminates the Ertel potential vorticity signal associated with advection of a non-uniform background state, thereby isolating the part of Ertel potential vorticity available for balanced-flow evolution. The effect of internal waves on QG flow is expressed concisely in a wave-averaged contribution to the materially conserved QG potential vorticity. We apply the theory by computing the wave-induced QG flow for a vertically propagating wave packet and a mode-one wave field, both in vertically bounded domains.

Existence, stability and formation of baroclinic tripoles in quasi-geostrophic flows

Hetons are baroclinic vortices able to transport tracers or species, which have been observed at sea. This paper studies the offset collision of two identical hetons, often resulting in the formation of a baroclinic tripole, in a continuously stratified quasi-geostrophic model. This process is of interest since it (temporarily or definitely) stops the transport of tracers contained in the hetons. First, the structure, stationarity and nonlinear stability of baroclinic tripoles composed of an upper core and two lower (symmetric) satellites are studied analytically for point vortices and numerically for finite-area vortices. The condition for stationarity of the point vortices is obtained and it is proven that the baroclinic point tripoles are neutral. Finite-volume stationary tripoles exist with marginal states having very elongated (figure-of-eight shaped) upper cores. In the case of vertically distant upper and lower cores, the latter can nearly join near the centre of the plane. These steady states are compared with their two-layer counterparts. Then, the nonlinear evolution of the steady states shows when they are often neutral (showing an oscillatory evolution); when they are unstable, they can either split into two hetons (by breaking of the upper core) or form a single heton (by merger of the lower satellites). These evolutions reflect the linearly unstable modes which can grow on the vorticity poles. Very tall tripoles can break up vertically due to the vertical shear mutually induced by the poles. Finally, the formation of such baroclinic tripoles from the offset collision of two identical hetons is investigated numerically. This formation occurs for hetons offset by less than the internal separation between their poles. The velocity shear during the interaction can lead to substantial filamentation by the upper core, thus forming small upper satellites, vertically aligned with the lower ones. Finally, in the case of close and flat poles, this shear (or the baroclinic instability of the tripole) can be strong enough that the formed baroclinic tripole is short-lived and that hetons eventually emerge from the collision and drift away.

Routes to Dissipation under Different Dynamical Conditions

AbstractIn this study, it is investigated how ageostrophic dynamics generate an energy flux toward smaller scales. Numerical simulations of baroclinic instability are used with varying dynamical conditions ranging from quasigeostrophic balance to ageostrophic flows. It turns out that dissipation at smaller scales by viscous friction is much more efficient if the flow is dominated by ageostrophic dynamics than in quasigeostrophic conditions. In the presence of ageostrophic dynamics, an energy flux toward smaller scales is observed while energy is transferred toward larger scales for quasigeostrophic dynamics. Decomposing the velocity field into its rotational and divergent components shows that only the divergent velocity component, which becomes stronger for ageostrophic flows, features a downscale flux. Variation of the dynamical conditions from ageostrophic dynamics to quasigeostrophic balanced flows shows that the forward energy flux and therefore the small-scale dissipation decreases as soon as the horizontal divergent velocity component decreases. A functional relationship between the small-scale dissipation and the local Richardson number is estimated. This functional relationship is used to obtain a global estimate of the small-scale dissipation of 0.31 ± 0.23 TW from a high-resolution realistic global ocean model. This emphasizes that an ageostrophic direct route to dissipation might be of importance in the ocean energy cycle.

Kraichnan–Leith–Batchelor similarity theory and two-dimensional inverse cascades

AbstractWe study the scaling properties and Kraichnan–Leith–Batchelor (KLB) theory of forced inverse cascades in generalized two-dimensional (2D) fluids (${\it\alpha}$-turbulence models) simulated at resolution $8192^{2}$. We consider ${\it\alpha}=1$ (surface quasigeostrophic flow), ${\it\alpha}=2$ (2D Euler flow) and ${\it\alpha}=3$. The forcing scale is well resolved, a direct cascade is present and there is no large-scale dissipation. Coherent vortices spanning a range of sizes, most larger than the forcing scale, are present for both ${\it\alpha}=1$ and ${\it\alpha}=2$. The active scalar field for ${\it\alpha}=3$ contains comparatively few and small vortices. The energy spectral slopes in the inverse cascade are steeper than the KLB prediction $-(7-{\it\alpha})/3$ in all three systems. Since we stop the simulations well before the cascades have reached the domain scale, vortex formation and spectral steepening are not due to condensation effects; nor are they caused by large-scale dissipation, which is absent. One- and two-point p.d.f.s, hyperflatness factors and structure functions indicate that the inverse cascades are intermittent and non-Gaussian over much of the inertial range for ${\it\alpha}=1$ and ${\it\alpha}=2$, while the ${\it\alpha}=3$ inverse cascade is much closer to Gaussian and non-intermittent. For ${\it\alpha}=3$ the steep spectrum is close to that associated with enstrophy equipartition. Continuous wavelet analysis shows approximate KLB scaling $\mathscr{E}(k)\propto k^{-2}~({\it\alpha}=1)$ and $\mathscr{E}(k)\propto k^{-5/3}~({\it\alpha}=2)$ in the interstitial regions between the coherent vortices. Our results demonstrate that coherent vortex formation (${\it\alpha}=1$ and ${\it\alpha}=2$) and non-realizability (${\it\alpha}=3$) cause 2D inverse cascades to deviate from the KLB predictions, but that the flow between the vortices exhibits KLB scaling and non-intermittent statistics for ${\it\alpha}=1$ and ${\it\alpha}=2$.

Dynamic Transitions of Quasi-geostrophic Channel Flow

The main aim of this paper is to study the dynamic transitions in flows described by the two-dimensional, barotropic vorticity equation in a periodic zonal channel. In [Z.-M. Chen et al., SIAM J. Appl. Math., 64 (2003), pp. 343--368], the existence of a Hopf bifurcation in this model as the Reynolds number crosses a critical value was proven. In this paper, we extend these results by addressing the stability problem of the bifurcated periodic solutions. Our main result is the explicit expression of a nondimensional parameter $\gamma$ which controls the transition behavior. We prove that depending on $\gamma$, the modeled flow exhibits either a continuous (Type I) or catastrophic (Type II) transition. Numerical evaluation of $\gamma$ for a physically realistic region of parameter space suggests that a catastrophic transition is preferred in this flow, which may lead to chaotic flow regimes.

Ribbon turbulence

We investigate the non-linear equilibration of a two-layer quasi-geostrophic flow in a channel with an initial eastward baroclinically unstable jet in the upper layer, paying particular attention to the role of bottom friction. In the limit of low bottom friction, classical theory of geostrophic turbulence predicts an inverse cascade of kinetic energy in the horizontal with condensation at the domain scale and barotropization in the vertical. By contrast, in the limit of large bottom friction, the flow is dominated by ribbons of high kinetic energy in the upper layer. These ribbons correspond to meandering jets separating regions of homogenized potential vorticity. We interpret these results by taking advantage of the peculiar conservation laws satisfied by this system: the dynamics can be recast in such a way that the initial eastward jet in the upper layer appears as an initial source of potential vorticity levels in the upper layer. The initial baroclinic instability leads to a turbulent flow that stirs this potential vorticity field while conserving the global distribution of potential vorticity levels. Statistical mechanical theory of the 112 layer quasi-geostrophic model predicts the formation of two regions of homogenized potential vorticity separated by a minimal interface. We explain that cascade phenomenology leads to the same result. We then show that the dynamics of the ribbons results from a competition between a tendency to reach the equilibrium state and baroclinic instability that induces meanders of the interface. These meanders intermittently break and induce potential vorticity mixing, but the interface remains sharp throughout the flow evolution. We show that for some parameter regimes, the ribbons act as a mixing barrier which prevents relaxation toward equilibrium, favouring the emergence of multiple zonal (eastward) jets.

Stable and symmetric quasi-geostrophic barotropic flows over topography

It is known that over a flat bottom, Taylor’s vortex is dynamically stable. The requirement of stability in the presence of a modulated topography with definite symmetry properties allows to modify the stream function of the vortex. The resulting solutions exhibit the same symmetry properties as the topography, and admit a variety of coupled vortices parametrized by the strength of the topography itself.

Comparison of the Effect of Parameterized Eddy Fluxes of Thickness and Potential Vorticity

AbstractParameterization of mesoscale eddies is an important problem of modern ocean dynamics and modeling. The most widely used scheme is the so-called Gent–McWilliams parameterization, which describes the eddy-induced transport of tracers, including temperature, density, and isopycnal thickness (TH). An alternative scheme, proposed by Green and Welander, deals with parameterizing eddy fluxes of potential vorticity (PV). Many recent studies propose using it, for it includes the effect of eddy Reynolds stresses that may influence mean flows. These two schemes are compared in the simplest configuration of two-layer quasigeostrophic channel flow, which enables analytical solutions for zonal-mean fields. It is shown how the parameterizations shape the zonally averaged zonal velocity profiles, with special attention paid to the role of the Reynolds stresses and momentum conservation. The zonally averaged zonal velocity profiles are sensitive to the amplitude and profiles of TH and PV diffusivities. For small enough diffusivities the TH parameterization may lead to solutions resembling those for the PV parameterization if it uses the diffusivity of the latter; that is, it may mimic the impact of the Reynolds stresses on the mean flow.

The stability of dipolar gyres on a beta-plane

When a source-sink dipole forces a fluid on a -plane limited by a western boundary, the linear steady solution can be obtained analytically and consists of zonally elongated gyres that extend west of the forcing and close as western boundary currents. The nondimensional parameter (with the zonal velocity of the flow and the distance between the source and sink) is used to characterize the nonlinearity of the flow. When reaches 0.1, the numerical shallow-water solution shows that the configuration with the source to the north of the sink becomes unstable, while the reverse configuration remains steady. Indeed, that reverse configuration remains steady for much larger values of the nonlinearity parameter , and begins to share some of the characteristics of a pure inertial circulation. The asymmetry of the stability properties of the two configurations, also found in the laboratory experiments of Colin de Verdière [Quasigeostrophic flows and turbulence in a rotating homogeneous fluid, 1977], is rationalized herein through the stability properties of the zonal central jet that flows between the source and sink. We consider, in turn, (i) the Kuo’s [J. Meteor. 1949, 6, 105–122] zero potential vorticity gradient necessary criteria (valid for an infinite zonal jet), (ii) enstrophy budgets and (iii) linear stability analysis of the mean flow. All three methods point out to the enhanced instability of the westward jet. We show that the transition regime has the characteristics of a super critical Hopf bifurcation.

Entropy Evolution and Uncertainty Estimation with Dynamical Systems

This paper presents a comprehensive introduction and systematic derivation of the evolutionary equations for absolute entropy H and relative entropy D, some of which exist sporadically in the literature in different forms under different subjects, within the framework of dynamical systems. In general, both H and D are dissipated, and the dissipation bears a form reminiscent of the Fisher information; in the absence of stochasticity, dH/dt is connected to the rate of phase space expansion, and D stays invariant, i.e., the separation of two probability density functions is always conserved. These formulas are validated with linear systems, and put to application with the Lorenz system and a large-dimensional stochastic quasi-geostrophic flow problem. In the Lorenz case, H falls at a constant rate with time, implying that H will eventually become negative, a situation beyond the capability of the commonly used computational technique like coarse-graining and bin counting. For the stochastic flow problem, it is first reduced to a computationally tractable low-dimensional system, using a reduced model approach, and then handled through ensemble prediction. Both the Lorenz system and the stochastic flow system are examples of self-organization in the light of uncertainty reduction. The latter particularly shows that, sometimes stochasticity may actually enhance the self-organization process.

Geostrophic adjustment with gyroscopic waves: stably neutrally stratified fluid without the traditional approximation

AbstractWe examine nonlinear geostrophic adjustment in a rapidly rotating (small Rossby number Ro) stably neutrally stratified (SNS) fluid consisting of a stratified upper layer with $N>f$ ($N$ is the buoyancy frequency, $f$ the Coriolis parameter) and a homogeneous lower layer, the density and other fields being continuous at the interface between the layers. The angular speed of rotation is non-parallel to gravity; the traditional and hydrostatic approximations are not used. The wave spectrum in the model consists of internal and gyroscopic waves. During the adjustment an arbitrary long-wave perturbation is split in a unique way into slow quasi-geostrophic (QG) and fast ageostrophic components with typical time scales $(Ro\, f)^{-1}$ and $f^{-1}$, respectively. The QG flow is governed by two coupled nonlinear equations of conservation of QG potential vorticity (PV) in the layers. The fast component is a sum of internal waves and inertial oscillations (long gyroscopic waves) confined to the homogeneous layer and modulated by an amplitude depending on coordinates and slow time. On times $t\sim (\, f\, Ro)^{-1}$ the slow component is not influenced by the fast one but the inertial oscillations amplitude is coupled to the QG flow and obeys an equation practically coinciding with that in the barotropic case (Reznik, J. Fluid Mech., vol. 743, 2014, pp. 585–605). A non-stationary boundary layer with large vertical gradients of horizontal velocity develops in the stratified layer near the interface to prevent penetration of the inertial oscillations into the stratified fluid; an analogous weaker boundary layer arises near the upper rigid lid. At large times the internal waves gradually decay because of dispersion and the resulting motion consists of the slow QG component and inertial oscillations confined to the barotropic lower layer.

Nonlinear energy transfers and phase diagrams for geostrophically balanced rotating-stratified flows.

Equilibrium statistical mechanics tools have been developed to obtain indications about the natural tendencies of nonlinear energy transfers in two-dimensional and quasi-two-dimensional flows like rotating and stratified flows in geostrophic balance. In this article we consider a simple model of such flows with a nontrivial vertical structure, namely, two-layer quasigeostrophic flows, which remain amenable to analytical study. We obtain the statistical equilibria of the system in the case of a linear vorticity-stream function relation, build the corresponding phase diagram, and discuss the most probable outcome of nonlinear energy transfers, both on the horizontal and on the vertical, in the presence of stratification and rotation.

Variability modes in core flows inverted from geomagnetic field models

The flow of liquid metal inside the Earth's core produces the geomagnetic field and its time variations. Understanding the variability of those deep currents is crucial to improve the forecast of geomagnetic field variations, which affect human spacial and aeronautic activities. Moreover, it may provide relevant information on the core dynamics. The main goal of this study is to extract and characterize the leading variability modes of core flows over centennial periods, and to assess their statistical robustness. To this end, we use flows that we invert from two geomagnetic field models (gufm1 and COV-OBS), and apply Principal Component Analysis and Singular Value Decomposition of coupled fields. The quasi geostrophic (QG) flows inverted from both geomagnetic field models show similar features. However, COV-OBS has a less energetic mean and larger time variability. The statistical significance of flow components is tested from analyses performed on subareas of the whole domain. Bootstrapping methods are also used to extract significant flow features required by both gufm1 and COV-OBS. Three main empirical circulation modes emerge, simultaneously constrained by both geomagnetic field models and expected to be robust against the particular a priori used to build them (large scale QG dynamics). Mode 1 exhibits three large vortices at medium/high latitudes, with opposite circulation under the Atlantic and the Pacific hemispheres. Mode 2 interestingly accounts for most of the variations of the Earth's core angular momentum. In this mode, the regions close to the tangent cylinder and to the equator are correlated, and oscillate with a period between 80 and 90 years. Each of these two modes is energetic enough to alter the mean flow, sometimes reinforcing the eccentric gyre, and other times breaking it up into smaller circulations. The three main circulation modes added together to the mean flow account for about 70% of the flows variability, 90% of the root mean square total velocities, and 95% of the secular variation induced by the total flows. Direct physical interpretation of the computed modes is not straightforward. Nonethe-less, similarities found between the two first modes and time/spatial features identified in different studies of core dynamics, suggest that our approach can help to pinpoint the relevant physical processes inside the core on centennial timescales.

Spectral non-locality, absolute equilibria and Kraichnan–Leith–Batchelor phenomenology in two-dimensional turbulent energy cascades

AbstractWe study the degree to which Kraichnan–Leith–Batchelor (KLB) phenomenology describes two-dimensional energy cascades in$\alpha $turbulence, governed by$\partial \theta / \partial t+ J(\psi , \theta )= \nu {\nabla }^{2} \theta + f$, where$\theta = {(- \Delta )}^{\alpha / 2} \psi $is generalized vorticity, and$\hat {\psi } (\boldsymbol{k})= {k}^{- \alpha } \hat {\theta } (\boldsymbol{k})$in Fourier space. These models differ in spectral non-locality, and include surface quasigeostrophic flow ($\alpha = 1$), regular two-dimensional flow ($\alpha = 2$) and rotating shallow flow ($\alpha = 3$), which is the isotropic limit of a mantle convection model. We re-examine arguments for dual inverse energy and direct enstrophy cascades, including Fjørtoft analysis, which we extend to general$\alpha $, and point out their limitations. Using an$\alpha $-dependent eddy-damped quasinormal Markovian (EDQNM) closure, we seek self-similar inertial range solutions and study their characteristics. Our present focus is not on coherent structures, which the EDQNM filters out, but on any self-similar and approximately Gaussian turbulent component that may exist in the flow and be described by KLB phenomenology. For this, the EDQNM is an appropriate tool. Non-local triads contribute increasingly to the energy flux as$\alpha $increases. More importantly, the energy cascade is downscale in the self-similar inertial range for$2. 5\lt \alpha \lt 10$. At$\alpha = 2. 5$and$\alpha = 10$, the KLB spectra correspond, respectively, to enstrophy and energy equipartition, and the triad energy transfers and flux vanish identically. Eddy turnover time and strain rate arguments suggest the inverse energy cascade should obey KLB phenomenology and be self-similar for$\alpha \lt 4$. However, downscale energy flux in the EDQNM self-similar inertial range for$\alpha \gt 2. 5$leads us to predict that any inverse cascade for$\alpha \geq 2. 5$will not exhibit KLB phenomenology, and specifically the KLB energy spectrum. Numerical simulations confirm this: the inverse cascade energy spectrum for$\alpha \geq 2. 5$is significantly steeper than the KLB prediction, while for$\alpha \lt 2. 5$we obtain the KLB spectrum.