LETTER TO THE EDITOR
Unstable interaction of gravityinertial waves with Rossby waves with application to solar system atmospheres
J. F. McKenzie^{1,2}  T. B. Doyle^{3}
1  King's College, University of Cambridge, UK
2 
School of Mathematical Sciences, University of KwaZuluNatal, South Africa
3 
Materials Research Division, iThemba LABS, Somerset West, South Africa
Received 12 October 2009 / Accepted 12 January 2010
Abstract
This letter reports on the important features of an
analysis of the combined theory of gravity  inertial  Rossby
waves on a plane in the Boussinesq approximation. In
particular, it is shown that the coupling
between higher frequency gravity  inertial waves and lower
frequency Rossby waves, arising from the accumulated influences
of the
effect, stratification characterized by the Väisäla  Brunt frequency
N, the Coriolis frequency f,
and the component of vertical propagation wave
number k_{z}, may lead to an unstable coupling between buoyancy 
inertial modes with westward propagating Rossby waves.
``Supersonic'' fast rotators (such as Jupiter) are predicted to be unstable in a
fairly narrow band of latitudes around their equators. The Earth is moderately
supersonic and exhibits instability within about
of its
equator. Slow ``subsonic'' rotators (e.g. Mercury, Venus, and the
Sun's corona) are unstable at all latitudes except those very close to
the poles where the
effect vanishes.
Key words: hydrodynamics  instabilities  planets and satellites: atmospheres  waves
1 Introduction
A brief report is presented on an instability which may have important consequences in the dynamics of rotating planetary atmospheres. The structure of each planetary atmosphere is unique, with properties determined by momentum balance and heating and cooling processes in the presence of heat conduction/transport and viscous stresses. The complex basic state of each planetary atmosphere (exosphere) is beyond the scope of this letter. Instead we assume a highly idealized situation in which the atmosphere is characterized by a density scale height H (determined by g and the temperature T) rotating with a frequency which, together with the planetary radius R, define a dimensionless (rotational Mach or Froude) number , which plays a critical role in the stability problem. The present analysis is governed by the dispersion equation for coupled inertialgravityRossby waves on a plane in the Boussinesq approximation. Under certain circumstances the coupling between westward propagating ``higher'' frequency inertialgravity waves and ``lower'' frequency Rossby waves drives the system unstable, with growth rates of order days. The stability analysis indicates that rapidly rotating planets (e.g. Jupiter) are unstable within a fairly narrow belt of latitudes around their equators, whereas slowly rotating atmospheres (e.g. Venus and Mercury) are unstable at all latitudes except near the poles. Earth and Mars are ``transonic'' and are unstable within latitudes of about and , respectively.2 Dispersion equation and stability analysis
The instability condition and the growth rate, arising from the unstable
coupling between buoyancyinertial modes and Rossby waves,
follow from the existence of complex roots for
in the ``midlatitude'' dispersion equation
(LonguetHiggins 1968; Pedlosky 1987; Gill 1982; and more recently
McKenzie 2009) given by
(1) 
where
(2) 
and
(3) 
Here is the Coriolis frequency (with the latitude of the plane), is the normalized wave frequency, and is the normalized horizontal wave number, where ), and V is a characteristic gravity  inertial speed given by
V=N/k_{z}.  (4) 
Here k_{z} is the vertical component of the wave number vector and N is the Väisäla  Brunt frequency, defined (Lighthill 1980) by
(5) 
in which H is the density scale height, is the pressure, the density, the velocity of sound, the adiabatic constant, the gas constant and the temperature.
We emphasize that although the assumed basic state is highly idealized, the instability analysis may, nevertheless, capture conditions on spherical shells in which N and may be regarded as constants, and for horizontal length scales consistent with the plane approximation. The analysis may, moreover, also apply to latitudinally sheared zonal flows (as is evident in the banded structure of Jupiter) in which the wave frequency is replaced by its Doppler shifted counterpart (see e.g. Mekki and McKenzie1977).
In defining the coupling parameter m in Eqs. (1) and (3) we have
taken the modulus of
and assigned a minus sign in the
right hand side of the dispersion equation (Eq. (1)) to make it
applicable to waves propagating in the second or third quadrants,
corresponding to a westward component of propagation (since plane
waves are proportional to
and with ). The
coupling parameter m is, in fact, the quantity ,
as
defined by Gill (1982). This parameter is normally regarded as
small because the Rossby wave frequency is generally very much
smaller than the gravityinertial wave frequency. In the present
work this approximation is not used and situations where the wave
frequencies of these modes may overlap in certain bands of wave
number are included. It will be noted that in the absence of the
effect (
,
where R is the
radius of the body) Eq. (1) reduces to the inertial  gravity mode
dispersion relation, namely,
(6) 
for wave frequencies , but for . Whereas for very low frequencies, the dispersion equation goes over to that of the Rossby wave, namely,
(7) 
The coupled dispersion equation (Eq. (1)) gives two positive roots for ; one of which represents a Rossby wave, whose frequency is increased above its uncoupled value , and the other is an inertial  gravity wave, whose frequency is decreased below its uncoupled value . The third root yields a negative value for (with ) corresponding to an eastward propagating gravityinertial wave. As the coupling ( ) increases the two positive roots move towards each other and coalesce where the right hand side of Eq. (1) is equal to the minimum value of the left hand side, namely, , which occurs at . A further increase in leads to complex conjugate roots (corresponding to instability) given approximately by
(8) 
where
(9) 
This approximation to the complex conjugate roots follows directly from the Taylor expansion of the LHS of Eq. (1) in the neighbourhood of the double root, at [neglecting terms of O ]. Since has a minimum value of 1 at the condition for instability is simply
m>1.  (10) 
When this is satisfied there is a band of unstable wave numbers lying between the two roots of
(11) 
given approximately by
(12) 
for m just moderately greater than unity. The central frequency of the instability (occurring at ) is the local Coriolis frequency divided by , i.e. , and the corresponding growth rate is . The solutions (real roots) of Eq. (1) for westward propagation are shown in Fig. 1, for the stable regime (m<1), the unstable regime (m>1) cases, and for the critical case m=1. The dashed line in the figure is the locus of . The instability growth rate , as a function of , for different values of the coupling parameter m>1 is shown in Fig. 2, in which the solid curves (exact solutions)may be compared with the dashed curves (approximate solutions), showing that the latter is accurate only near marginal stability ( ).
Figure 1: The diagnostic diagram, ( , ) plot, for the coupled modes, including a stable case (m=1/2), an unstable case (m=2), and the critical case (m=1). The higher frequency curve segments (above the locus ) describe the gravity  inertial mode, and the lower frequency segments the Rossby mode. The coupling of these modes leads to a convective instability for m> 1. 

Open with DEXTER 
Figure 2: The normalized growth rate as a function of for m = 1.2, 2 and 3. The solid and dashed curves are derived, respectively, from the exact complex conjugate roots of Eq. (1)) and from the approximate form given by Eq. (8). 

Open with DEXTER 
The instability condition m>1 may be written as
(13) 
which, on squaring, translates into a biquadratic for . This may be cast in the interesting form:
(14) 
(15) 
and
(16) 
where again H is the density scale height. The expressions for follow from the definitions of Eq. (5) for N, H and . Note that the Lamb parameter introduced by LonguetHiggins (1968) is related to the rotational Mach number by . The Mach number plays a critical role in the instability condition. This is shown in Fig. 3 as the region below the curve of the critical latitudes as a function of an equivalent Mach number M, defined by relations (14)(16), in terms of and the configuration of the wave, through the parameter .
Figure 3: The instability condition, relations [1316], with the critical latitude as a function of the effective rotational Mach number, M, separating the parameter space into stable and unstable regions. 

Open with DEXTER 
3 Application to planetary atmospheres
The results of the present analysis, as shown in Fig. 3, predict that for supersonically rotating planets ( ), such as the outer giants (e.g. Jupiter), instability is confined to a fairly narrow belt of latitudes around their equators. For subsonically rotating bodies (e.g. the Sun, Mercury and Venus), for which , instability occurs at almost all latitudes, with the fastest growing mode near their equators. The Earth and Mars are ``transonic'' planets for which is around and , respectively. Note that in the ``hypersonic'' Jupiter instability is confined to approximately about its equator. The present analysis breaks down near the equator where the Coriolis frequency goes to zero. In this region a special treatment is required as a consequence of the development of a waveguide system (Moore & Philander 1977; Cane & Sarachik 1976; and Maas & Harlander 2007). These analyses indicate that ``globally'' the instability disappears, while ``locally'' it may still exist within a JWKB approximation. This point requires further study and clarification.
Finally we note that a similar analysis pertains to oceans except
that the effective Mach number M in the instability condition
(Eq. (14)) is now given by
(17) 
where h is the ocean depth and is the shallow water speed. In an Earth's ocean of depth h=6 km this speed is roughly 250 m/s, yielding . Purely westward propagating waves ( ) are, therefore, unstable in a band of latitudes of about of the equator.
Since this instability feeds off rotational kinetic energy and gravitational buoyancy its nonlinear evolution may play an important role in atmospheric and ocean dynamics over time scales of, or less than, a few planetary days.
AcknowledgementsJ.F. McK. wishes to thank the PeiLing Chan Chair, Department of Physics, University of Alabama, for support, and J.F. McK. and T.B.D. thank the Foundation for Research and Development of South Africa (under Grant number 2053776) for support.
References
 Cane, M. A., & Sarachik, E. S. 1976, J. Mar. Res., 35, 395 (In the text)
 Gill, A. E. 1982, AtmosphereOcean Dynamics, International Geophys. Ser. 30 (In the text)
 Lighthill, J. 1980, Waves in Fluids (Cambridge University Press) (In the text)
 LonguetHiggins, M. S. 1968, Phil. Trans. Roy. Soc., 262, 511 [NASA ADS] [CrossRef] (In the text)
 Maas, L. R. M., & Harlander, U. 2007, J. Fluid Mech., 570, 47 [NASA ADS] [CrossRef] (In the text)
 McKenzie, J. F. 2009, Ann. Geophys., 27, 4221 [NASA ADS] [CrossRef] (In the text)
 Mekki, O. M., & McKenzie, J. F. 1977, Phil. Trans. R. Soc. London Ser. A, 287, 115 [NASA ADS] [CrossRef] (In the text)
 Moore, D. W., & Philander, S. G. H. 1977, Modelling of the tropic ocean circulation. In The Sea (Chap. 8), ed. Goldberg et al. (New York: Inter science) (In the text)
 Pedlosky, J. 1987, Geophys. Fluid Dynamics (Springer Verlag) (In the text)
All Figures
Figure 1: The diagnostic diagram, ( , ) plot, for the coupled modes, including a stable case (m=1/2), an unstable case (m=2), and the critical case (m=1). The higher frequency curve segments (above the locus ) describe the gravity  inertial mode, and the lower frequency segments the Rossby mode. The coupling of these modes leads to a convective instability for m> 1. 

Open with DEXTER  
In the text 
Figure 2: The normalized growth rate as a function of for m = 1.2, 2 and 3. The solid and dashed curves are derived, respectively, from the exact complex conjugate roots of Eq. (1)) and from the approximate form given by Eq. (8). 

Open with DEXTER  
In the text 
Figure 3: The instability condition, relations [1316], with the critical latitude as a function of the effective rotational Mach number, M, separating the parameter space into stable and unstable regions. 

Open with DEXTER  
In the text 
Copyright ESO 2010