Advances in Meteorology

Volume 2013 (2013), Article ID 981271, 11 pages

http://dx.doi.org/10.1155/2013/981271

Research Article

The Sensitivity of Characteristics of Large Scale Baroclinic Unstable Waves in Southern Hemisphere to the Underlying Climate

Sergei Soldatenko and Chris Tingwell

Centre for Australian Weather and Climate Research, 700 Collins Street, Docklands, VIC 3008, Australia

Received 26 September 2013; Accepted 18 November 2013

Academic Editor: Luis Gimeno

Copyright © 2013 Sergei Soldatenko and Chris Tingwell. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The sensitivity of the main characteristics of baroclinically unstable waves with respect to fundamental parameters of the atmosphere (the static stability parameter and vertical shear of a zonal wind ) is theoretically explored. Two types of waves are considered: synoptic scale waves and planetary scale (ultralong) waves based on an Eady-type model and model with vertically averaged primitive equations. Sensitivity functions are obtained that estimate the impact of variations in and on the growth rate and other characteristics of unstable waves and demonstrate that waves belonging to the short-wave part of the spectrum of unstable waves are more sensitive to changes in the static stability parameter than waves belonging to the long-wave part of the spectrum. The obtained theoretical results show that the increase of the static stability and decrease of the meridional temperature gradient in midlatitude baroclinic zones in some areas of the southern hemisphere lead to a slowing of the growth rate of baroclinic unstable waves and an increasing wavelength of baroclinic unstable wave maximum growth rate, that is, a spectrum shift of unstable waves towards longer wavelengths. These might affect the favorable conditions for the development of baroclinic instability and, therefore, the intensity of cyclone generation activity.

1. Introduction

Static stability and the meridional temperature gradient (MTG) are among the most important fundamental parameters characterizing the state of the atmosphere and, in particular, midlatitude large-scale eddy dynamics [1, 2]. Static stability and MTG play a significant role in the development of baroclinic instability which is the dominant mechanism for generating large-scale atmospheric eddies (cyclones) that form the storm tracks in midlatitudes. The physical nature of baroclinic instability is well understood and explained in the scientific literature, including text books on dynamic meteorology (e.g., [1–4]). Baroclinic instability can be viewed as sloping convection where growing perturbations draw upon the available potential energy which is proportional to a meridional temperature gradient. Since the publication of the pioneering theoretical works of Charney [5] and Eady [6], in which the fundamental baroclinic mechanism of the atmospheric large-scale instability was first described, many scientific papers have been published that examine the growth of initially infinitesimal perturbations in the atmosphere and ocean caused by baroclinic effects. Both linear theory for the onset of baroclinic instability and its nonlinear saturation have been explored in many research articles. A thorough review of baroclinic instability is presented in [7] and to some extent in [8, 9]. Theoretical models of baroclinic instability typically represent linearized dynamics equations and the instability problem is examined as an eigenvalue problem. The nonmodal instability analysis in both linear and nonlinear formulations is more general than eigenanalysis. This technique suggests the solution of a Cauchy (initial value) problem (e.g., [10–12]). Analytical weakly nonlinear theories of baroclinic instability represent a further extension of research (e.g., [13–16]) which focuses on the finite-amplitude behaviour of unstable baroclinic waves. Other areas of studies have explored the life-cycle behaviour of baroclinic unstable waves to understand how initially infinitesimal perturbations grow to large but finite amplitude modifying the mean flow. This kind of research requires numerical integration of the atmospheric nonlinear model equations (e.g., [17–20]).

It is important to underline the significant role of the two-layer model of Philips [21, 22] and its modifications (e.g., [13, 23–28]) in developing baroclinic instability theories. However, for this purpose a two-layer model was usually used in the quasigeostrophic approximation. A generalized baroclinic instability analysis based on two-level primitive-equation model has been described in [29] and applied to detailed consideration of the baroclinic instability mechanism.

There is evidence of an increase over the past decades of the static stability in the extratropics [30] and a pole-ward movement of the midlatitude precipitation zones and storm tracks which are uniquely linked to zones with strong MTG-baroclinic zones (e.g., [31–35]). Recent studies have also indicated an increase in intensity of extratropical cyclones and a decrease in frequency [36–44]. Changes in geographical locations and intensity of storm tracks are more distinct in the southern hemisphere (SH) (e.g., [45]) affecting the essential features of weather patterns over large territories such as Australia [46–54] and indicating a change of favourable conditions for baroclinic instability because the key role of baroclinic instability in the development of midlatitude cyclones is well established.

To study how baroclinic instability has changed in the SH in recent decades, a comparative analysis of climates for periods 1949–1968 and 1975–1994 has been carried out in [50, 54] using National Centres for Environmental Prediction (NCEP) reanalyses and the European Centre for Medium-Range Weather Forecasting (ECMWF) Reanalysis (ERA40) data. As a measure of the baroclinic instability the Phillips criterion [55] adapted to spherical geometry was chosen:

where and are the zonal wind velocities at 300 hPa and 700 hPa isobaric levels, respectively, is the static stability measure for a given reference state, is the specific heat of air at constant pressure, is the angular rotation speed of the earth, is the earth’s radius, is a dimensionless constant, and is a latitude. J. S. Frederiksen and C. S. Frederiksen [50] obtained a significant decreasing trend in baroclinic instability over the middle latitudes of the SH. The most significant negative trend is localized between 30S and 40S. For instance, the difference between the two periods, 1949–1968 and 1975–1994, reached a maximum of 17% in the SH subtropical jet stream (for the July climate). In addition it was found that further poleward the baroclinic instability intensified creating favourable conditions for cyclogenesis. Using a two-level linearized primitive equation model, Frederiksen et al. [51] studied the growth rate of unstable modes for the July reference state averaged over the above mentioned periods. They identified a 30% reduction in the growth rate of cyclone-scale modes between the two periods.

It should be pointed out that, since the beginning of the 1980s, exploration of the essential features of the atmospheric general circulation in the SH has been facilitated by the availability of plentiful satellite upper-air data. In particular, the importance of baroclinic instability and large-scale eddies in the formation of zonal wind characteristics, such as a tropospheric double-jet phenomenon, was studied in [56, 57].

Taking into account the significant role of baroclinic instability in the development of midlatitude cyclones and the intensity changes of baroclinic instability in some areas of the SH in recent decades, this paper examines the sensitivity of the main characteristics of baroclinically unstable waves (e.g., the growth rates of unstable waves as function of wavelength) to fundamental atmospheric parameters: the static stability parameter which is characterized by temperature lapse rate , and zonal wind vertical shear which, by thermal wind balance, characterizes the meridional temperature gradient. Two classes of waves [55] are considered. The first one is characterized by the Rossby number and the representative horizontal length scale which is smaller than the Earth’s radius , that is, . Here is a horizontal velocity scale and is the Coriolis parameter. For this class of waves , where is the vertical component of relative vorticity and is the horizontal divergence. Synoptic scale waves belong to this type of motion [58]. The second class, which includes planetary, or ultralong, waves [59], is characterized by and . For these waves the magnitude of horizontal divergence is comparable with the vertical component of vorticity; that is, . To study synoptic scale waves, the Eady-type model is used with uniform zonal wind shear between upper and lower boundaries on an -plane. In this context parameters and are considered to be variables that control the development of baroclinic instability in the atmosphere. Ultralong waves are investigated based on a model with vertically averaged equations [60, 61] with a -plane approximation. Simplified models such as the Eady model of baroclinic instability and models with vertically averaged equations, despite their simplicity, allow solutions to be obtained that clearly illustrate real physical processes in the atmosphere.

2. Synoptic Scale Baroclinically Unstable Waves

2.1. The Model Equations

We consider the inviscid primitive equation atmospheric model in normalized isobaric coordinates on an -plane in the following form:

The state variables of the model are the horizontal velocity vector , the vertical pressure velocity , where is pressure, the geopotential , and the temperature . The operator is applied to the horizontal coordinates and , directed eastward and northward, respectively. The normalized pressure , where  hPa is a “standard” pressure approximately equal to the surface pressure, is taken as the vertical coordinate, while the time is denoted by . Other notations are the diabatic heating rate per unit time per unit mass , the gas constant , the unit vector in the vertical direction , and the reference state static stability measure in the normalized isobaric coordinate system:

where is a reference temperature, is the gravity acceleration, is the dry adiabatic lapse rate, and is the reference state lapse rate. We employ the -plane approximation so that the Coriolis parameter is a constant with being the latitude of interest. Hereafter we consider only adiabatic process and thus assume zero heating rate . The following boundary conditions are used for the pressure velocity:

The atmospheric reference state, defined by , , , , and , is steady and satisfies the following relations:

where . By substituting (5) into the set of (2), one can see that (5) is a solution of (2) that describes the zonal flow:

which matches the specified distribution of the zonally averaged temperature and represents thermal wind balance. To consider only the baroclinic mechanism of the atmospheric instability, meridional variability of the basic zonal flow is excluded. In other words the barotropic impact on the instability of the basic zonal flow is not taken into account. Therefore we assume that the velocity of the basic zonal flow does not depend on the horizontal -coordinate; that is, . Thus, the problem now is the study of the instability of the basic zonal flow (6) with respect to infinitesimal perturbations. For this purpose the system (2) is linearized around the basic state (5). Representing the state variables as , where is a basic state and is an infinitesimal perturbation, and taking into account the hydrostatic equations and the thermal wind relationship (6), the linearized system can be written as

The static stability parameter is expressed as

where is a specific volume and is a reference potential temperature. Suppose that . Applying the method of separation of variables, we assume the solutions of the form

where is a function of only, is a wave number, and is a phase velocity of perturbations which in general is a complex value: . After substituting (9) into (7), we can finally obtain the following single equation for :

Similar equations were considered in a number of publications (e.g., [4, 62, 63]). The analytical solution of this equation can be obtained by eliminating gravity waves from the consideration and assuming that . For this particular case, (10) can be transformed into

The boundary conditions for are specified as

Equation (11), together with boundary conditions (12), represents the eigenvalue problem for the complex phase speed . The fundamental solutions of (11) are expressed through the Bessel functions of the first and second kinds. With the homogeneous boundary conditions (12), the following two discrete eigenvalues can be obtained [4]:

where

The phase velocity will amplify exponentially if has an imaginary part . From (13) we can see that this will occur if the discriminant in (13) is less than zero:

which gives by the Newton’s iteration algorithm the necessary condition for instability: . Besides two discrete eigenvalues (13), the eigenvalue problem (11)-(12) has a continuous spectrum of eigenvalues that are real and therefore can be neglected in the problem of baroclinic instability [4]. The growth rate of unstable waves is calculated by the following expression:

As we can see from (16), at a given latitude the growth rate is a function of (static stability), (wind shear), and quantity that depends on the horizontal wavenumber and .

2.2. Impact of Static Stability and Vertical Wind Shear on Baroclinic Instability

Within the Eady problem framework, the static stability parameter and the vertical wind shear represent the main control variables. By varying and , one can obtain estimates of the impact of these parameters on the development of baroclinic instability in the atmosphere. In this research, parameters corresponding to the basic state are given the following values:  ms⁻¹ [50] and  m² Pa⁻² s⁻² [1]. These parameter values can be used as an approximation to describe the zonal-averaged atmospheric conditions for JJA (June, July, and August) in the SH [50]. The latitude of interest is assumed to be  S which gives  s⁻¹.

Figure 1 shows plot of Eady growth rate versus zonal wavenumber obtained with (16). The growth rate has a short-wave instability cutoff beyond which waves are stable. Let be the wavelength that corresponds to a short-wave cutoff. Value of can be obtained from (14) when which gives  km. To calculate the wavelength of maximum growth rate , one can take and set the result equal to zero which gives . Then by using (14) we can obtain  km.

Figure 1: Growth rate versus zonal wavenumber for  m² Pa⁻² s⁻² and  m s⁻¹.

The influence of the static stability parameter on the wavelength of maximum growth rate and the short wave cut-off are shown in Figure 2. In general, an increase in the parameter leads to an increase in both and . The functional dependences between and , and between and are almost linear: departure of static stability parameter from its nominal value  m² Pa⁻² s⁻² results in about % change for both and with respect to the nominal value . For instance, if , then and  km, and if , then and  km.

Figure 2: Length of waves with maximum growth rates and short-wave cutoff , as functions of the static stability parameter .

Figure 3 illustrates the growth rate of unstable waves versus the static stability parameter at different values of . Parameters and influence the growth rate in the opposite direction: growth rate decreases if increases and if decreases. Note that the decrease of the parameter indicates the weakening of the intensity of the baroclinic zone, that is, reduction of the MTG. In nature both of these processes take place, which leads to a synergistic effect. For instance, if decreases by 10% and the static stability parameter increases by 10%, the growth rate decreases by 14%.

Figure 3: Growth rate versus static stability parameter for different values of parameter (units m s⁻¹).

Since is a nonlinear function of (16), to estimate the influence of infinitesimal perturbations in on variations in , the sensitivity function

and the relative sensitivity function

can be used. The function shows changes in due to variations in . The relative sensitivity function is used to compare model parameters to find out what parameter is the most important for a certain percent change in the parameters. Sensitivity functions (17) and (18) are evaluated in the vicinity of some nominal value of the parameter . We can select several nominal values to cover some range of changes in . Differentiating (16) with respect to control parameter , we can obtain the expression for

Sensitivity versus zonal wavenumbers for different values of with  m s⁻¹ are shown in Figure 4. The absolute value of the sensitivity of with respect to exponentially increases with decreasing wavelength. For planetary scale waves (zonal wave numbers 1–4), the absolute value of the sensitivity of with respect to is palpably less than sensitivity for synoptic scale waves (zonal wave numbers ≥5). Absolute and relative sensitivity functions and calculated for different values of for various zonal wave numbers are shown in Tables 1 and 2, respectively.

Table 1: Absolute sensitivity as a function of zonal wavenumber for different values of static stability parameter .

 m² Pa⁻² s⁻²
Zonal wave number

1
2
3
4
5
6
7

−0.0157
−0.1251
−0.4175
−0.9772
−1.8879
−3.2470
−5.2005

−0.0157
−0.1245
−0.4139
−0.9675
−1.8794
−3.3014
−5.5927

−0.0157
−0.1240
−0.4110
−0.9626
−1.8979
−3.4964
−6.9371

−0.0157
−0.1234
−0.4088
−0.9629
−1.9509
−3.9665
−23.0031

−0.0157
−0.1230
−0.4071
−0.9691
−2.0540
−5.3057

Table 2: Relative sensitivity as a function of zonal wavenumber for different values of static stability parameter .

 m² Pa⁻² s⁻²
Zonal wave number

1
2
3
4
5
6
7

−0.0062
−0.0250
−0.0575
−0.1057
−0.1735
−0.2692
−0.4095

−0.0093
−0.0379
−0.0881
−0.1656
−0.2837
−0.4749
−0.8371

−0.0124
−0.0509
−0.1201
−0.2375
−0.4220
−0.8004
−1.9939

−0.0156
−0.0642
−0.1539
−0.3087
−0.6071
−1.4399
−31.4243

−0.0187
−0.0777
−0.1898
−0.3974
−0.8759
−3.4322

Figure 4: Sensitivity function versus zonal wavenumber for different nominal values of the static stability parameter (units 10⁶ m² Pa⁻² s⁻²).

The expression for sensitivity function can be easily obtained by differentiating (16) with respect to control parameter :

The function versus zonal wave number for different is plotted in Figure 5. It is clear to see that for a given value of the parameter the graph of function is very much like the classic picture of the growth rates versus zonal wavenumber [5]. It is interesting that the relative sensitivity function does not depend on the wavelength (wavenumber) and for all of the unstable waves is equal to unity:

Since relative sensitivity functions allow direct comparison of the importance of model parameters on the growth rate , we can see that because the parameter (i.e., the meridional temperature gradient) is more important than the static stability parameter except for  m² Pa⁻² s⁻² for waves with (see Table 2). Wherever the midlatitude values in Table 2 are less than one, the growth rate is more sensitive to the meridional temperature gradient (i.e., ) than the static stability ().

Figure 5: Sensitivity function versus zonal wavenumber for different nominal values of the static stability parameter (units 10⁶ m² Pa⁻² s⁻²).

The obtained results are consistent with observations [50, 51, 54]: an increase in static stability and a decrease of the MTG have occurred over the past few decades in some areas of the SH, which has led to a decrease in the growth rate of baroclinic unstable waves, a shift of the spectrum of unstable waves in the long wavelength part of spectrum, and a weakened intensity of cyclogenesis. Naturally, these changes affect favourable conditions for the development of baroclinic instability and the essential features of weather patterns over large territories, particularly over Australia.

3. Planetary Scale Waves

To study the influence of the static stability parameter on the dynamics of planetary scale (ultralong) waves, a thin film approximation is applied. This approximation employs a specific averaging technique over the vertical coordinate to the system of primitive equations [60]. As a result, a two-dimensional set of equations can be obtained that describes the dynamics of a two-dimensional baroclinic film. These equations reproduce all the wavelike solutions that correspond to the main weather-forming modes of three-dimensional models and therefore can be used in theoretical studies of large-scale dynamic processes in the atmosphere. The system of vertically averaged equations can be written as [60]

where . For instance, if the original primitive equations are written in the Phillips’ vertical coordinate system [64], the operator for vertical averaging is introduced as and state variables are represented as . Equations (22) are obtained by neglecting the orography and terms , , and [65]. A detailed linear analysis of the vertically averaged equation (22) is represented in [60]. In particular, two types of wave solutions were found: fast waves that propagate westward and slow waves that move eastward. Within the framework of this model, ultralong waves are always neutral for any vertically averaged zonal wind velocity. Indeed, linearizing (22) around the following basic state:

and assuming the beta plane approximation , where , and representing the solution in the form (9), one can finally obtain under different asymptotics the following expressions for four wave solutions [60]:

(a)acoustic waves

(b)Rossby wave

(c)baroclinic wave

Here and . These results however are valid only for the specific case of a neutral atmosphere with [60, 61]. To take into account the atmospheric static stability on the behaviour of ultralong waves the polytropic model of the atmosphere can be used for which

where is the temperature at the surface and is a vertical temperature gradient. Integrating (27) with respect to vertical coordinate, we can obtain [60]. Assuming the geostrophic approximation on a -plane, the set of vertically averaged equations can be written as [61]

Here , , , , and the operator . The basic state is defined as a stationary solution of system (28) for which

or in other words

Linearizing (28) around the basic state (30), the following cubic characteristic equation can be obtained in which the second order terms are neglected [61]:

where and . If the discriminant of this equation is positive, then the wave solution is unstable. The domain of zonal flow instability can be found numerically (see diagram in [61]). In Figure 6, we reproduce only for the 1st quadrant of a Cartesian plane the domain of instability calculated as a function of vertically averaged zonal wind velocity and dimensionless temperature lapse rate .

Figure 6: Domain of instability (filled) as a function of dimensionless temperature lapse rate and vertically averaged zonal wind velocity .

The imaginary part of phase velocity which characterises the growth rate of unstable waves is displayed in Figure 7 as a function of dimensionless temperature lapse rate for different values of vertically averaged zonal wind velocity . A maximum phase velocity exists for given values of ; that is, dependent on the ratio of . For instance, if  m s⁻¹ then the maximum value  m s⁻¹ is reached at . Figure 6 shows that increasing vertically averaged zonal wind is associated with increasing . This is further evident in Figure 8, which shows as a function of for a range of values. The lower and the larger , that is, , increases with decreasing static stability.

Figure 7: Imaginary part of phase speed versus dimensionless temperature lapse rate for different values of vertically averaged zonal wind velocity .

Figure 8: Imaginary part of phase speed versus vertically averaged zonal wind velocity for different values of dimensionless temperature lapse rate .

4. Concluding Remarks

We have studied theoretically the impact of variations in the static stability parameter and zonal wind shear on the characteristics of baroclinically unstable waves of synoptic scales using Eady-type model with the uniform between upper and lower boundaries on an -plane. Quantitative estimates of variations in and on the growth rate , wavelength of maximum growth rate , and short-wave cutoff were obtained.

Analytical expressions are derived for sensitivity functions for the growth rate with respect to variations in static stability parameter and wind shear velocity. These expressions allow estimating to a first-order approximation the influence of changes in and on . Analytical expressions for relative sensitivity functions allow estimating the significance of variations in and on the growth rate of baroclinically unstable waves with a given zonal wave number.

To study the impact of variations in atmospheric static stability and zonal wind velocity on the instability of planetary scale waves, the model with vertically averaged primitive equations with -plane approximation was applied. As control parameters, we have used dimensionless temperature lapse rate and vertically averaged zonal wind velocity . We have estimated the influence of and on the imaginary part of phase speed , which was used as a measure of instability.

The obtained results are qualitatively consistent with changes in the essential weather patterns that occurred over the last several decades in some areas of the SH and, in particular, over Australia (e.g., [49, 50, 52–54]). Climate change results suggest SH midlatitude static stability may increase and the MTG (the vertical wind shear ) may decrease, which according to our linear theoretical models leads to a slowing of the growth rate of baroclinic unstable waves and an increasing wavelength of baroclinic unstable wave with maximum growth rate , that is, a spectrum shift of unstable waves towards longer wavelengths. These might affect the favourable conditions for the development of baroclinic instability and, therefore, the rate of cyclogenesis and a reduction in cyclone intensity. The obtained sensitivity functions demonstrate that waves belonging to the short-wave part of the spectrum of unstable waves are more sensitive to changes in the static stability parameter than waves belonging to the long-wave part of the spectrum.

To obtain more realistic estimates of the sensitivity of the growth rate of unstable waves with respect to static stability parameter and MTG, numerical modeling based on a full GCM is required. It is hoped to carry out such work in the future.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors thank Dr. I. Pisnichenko for clarification of certain questions relevant to this paper.

References

J. R. Holton, An Introduction to Dynamic Meteorology, Academic Press, 3rd edition, 1992.

G. K. Vallis, Atmospheric and Oceanic Fluid Dynamics, Cambridge University Press, 2006.

M. Mak, Atmospheric Dynamics, Cambridge University Press, 2011.

V. P. Dymnikov, Stability and Predictability of the Large Scale Atmospheric Processes, Institute of Numerical Mathematics of the Russian Academy of Sciences, Moscow, Russia, 2007.

J. G. Charney, “The dynamics of long waves in baroclinic westerly current,” Journal of Meteorology, vol. 4, pp. 135–162, 1947.

E. T. Eady, “Long waves and cyclone waves,” Tellus, vol. 1, pp. 33–52, 1949.

R. T. Pierrehumbert and K. L. Swanson, “Baroclinic instability,” Annual Review of Fluid Mechanics, vol. 27, no. 1, pp. 419–467, 1995.

R. Gall, “A comparison of linear instability theory with the eddy statistics of a general circulation model,” Journal of the Atmospheric Sciences, vol. 33, no. 3, pp. 349–373, 1976.

I. M. Held, “Progress and problems in large-scale atmospheric dynamics,” in The Global Circulation of the Atmosphere, T. Schneider and A. Sobel, Eds., pp. 1–21, Princeton University Press, Princeton, NJ, USA, 2007.

B. Farrell, “Modal and non-modal baroclinic waves,” Journal of the Atmospheric Sciences, vol. 41, no. 4, pp. 668–673, 1984.

M. K. Tippett, “Transient moist baroclinic instability,” Tellus A, vol. 51, no. 2, pp. 273–288, 1999.

D. Hodyss and R. Grotjahn, “Nonmodal and unstable normal mode baroclinic growth as a function of horizontal scale,” Dynamics of Atmospheres and Oceans, vol. 37, no. 1, pp. 1–24, 2003.

J. Pedlosky, “Finite amplitude baroclinic waves,” Journal of the Atmospheric Sciences, vol. 27, no. 1, pp. 15–30, 1970.

R. D. Romea, “The effects of friction and beta on finite amplitude baroclinic waves,” Journal of the Atmospheric Sciences, vol. 34, pp. 1689–1695, 1977.

J. Pedlosky, “Finite-amplitude baroclinic waves at minimum critical shear,” Journal of the Atmospheric Sciences, vol. 39, no. 3, pp. 555–562, 1982.

T. Warn and P. Gauthier, “Potential vorticity mixing by marginally unstable baroclinic disturbances,” Tellus A, vol. 41, no. 2, pp. 115–131, 1989.

A. J. Simmons and B. J. Hoskins, “The lifecycles of some nonlinear waves,” Journal of the Atmospheric Sciences, vol. 35, pp. 414–432, 1978.

A. J. Simmons and B. J. Hoskins, “Barotropic influences of the growth and decay of nonlinear baroclinic waves,” Journal of the Atmospheric Sciences, vol. 37, pp. 1679–1684, 1980.

S. B. Feldstein and I. M. Held, “Barotropic decay of baroclinic waves,” Journal of the Atmospheric Sciences, vol. 46, pp. 1679–1684, 1989.

C. D. Thorncroft, B. J. Hoskins, and M. E. McIntyre, “Two paradigms of baroclinic-wave life-cycle behaviour,” Quaterly Journal, vol. 119, no. 509, pp. 17–55, 1993.

N. A. Phillips, “A simple three-dimensional model for the study of largescale extratropical flow patterns,” Journal of Meteorology, vol. 8, pp. 381–394, 1951.

N. A. Phillips, “Energy transformations and meridional circulations associated with simple baroclinic waves in a two-level, quasi-geostrophic model,” Tellus, vol. 6, pp. 273–286, 1954.

L. M. Polvani and J. Pedlosky, “The effect of dissipation on spatially growing nonlinear baroclinic waves,” Journal of the Atmospheric Sciences, vol. 45, no. 14, pp. 1977–1989, 1988.

T. G. Shepherd, “Nonlinear saturation of baroclinic instability. Part I: the two-layer model,” Journal of the Atmospheric Sciences, vol. 45, no. 14, pp. 2014–2025, 1988.

G. Balasubramanian and M. K. Yau, “Baroclinic instability in a two-layer model with parameterized slantwise convection,” Journal of the Atmospheric Sciences, vol. 51, no. 7, pp. 971–990, 1994.

A. Wiin-Nielsen, “On the structure of atmospheric waves in middle latitudes,” Atmosfera, vol. 16, no. 2, pp. 83–102, 2003.

D. D. Holm and B. A. Wingate, “Baroclinic instabilities of the two-layer quasigeostrophic alpha model,” Journal of Physical Oceanography, vol. 35, no. 7, pp. 1287–1296, 2005.

J. Egger, “Baroclinic instability in the two-layer model: interpretations,” Meteorologische Zeitschrift, vol. 18, no. 5, pp. 559–565, 2009.

J. Lamberts, G. Lapeyre, and V. Zeitlin, “Moist versus dry baroclinic instability in a simplified two-layer atmospheric model with condensation and latent heat release,” Journal of the Atmospheric Sciences, vol. 69, pp. 1405–1426, 2012.

D. M. W. Frierson, “Robust increases in midlatitude static stability in simulations of global warming,” Geophysical Research Letters, vol. 33, no. 24, Article ID L24816, 2006.

N. M. J. Hall, B. J. Hoskins, P. J. Valdes, and C. A. Senior, “Storm tracks in a high-resolution GCM with doubled carbon dioxide,” Quarterly Journal, vol. 120, no. 519, pp. 1209–1230, 1994.

S. B. Feldstein, “The recent trend and variance increase of the annular mode,” Journal of Climate, vol. 15, no. 3, pp. 88–94, 2002.

G. J. Marshall, “Trends in the southern annular mode from observations and reanalyses,” Journal of Climate, vol. 16, pp. 4134–4143, 2003.

G. M. Ostermeier and J. M. Wallace, “Trends in the North Atlantic Oscillation-Northern Hemisphere annular mode during the twentieth century,” Journal of Climate, vol. 16, no. 2, pp. 336–341, 2003.

D. W. J. Thompson and S. Solomon, “Interpretation of recent Southern Hemisphere climate change,” Science, vol. 296, no. 5569, pp. 895–899, 2002.

S. J. Lambert, “The effect of enhanced greenhouse warming on winter cyclone frequencies and strengths,” Journal of Climate, vol. 8, no. 5, pp. 1447–1462, 1995.

F. Lunkeit, K. Fraedrich, and S. E. Bauer, “Storm tracks in a warmer climate: sensitivity studies with a simplified global circulation model,” Climate Dynamics, vol. 14, no. 11, pp. 813–826, 1998.

D. J. Lorenz and D. L. Hartmann, “Eddy-zonal flow feedback in the Southern Hemisphere,” Journal of the Atmospheric Sciences, vol. 58, no. 21, pp. 3312–3327, 2001.

Q. Geng and M. Sugi, “Possible change of extratropical cyclone activity due to enhanced greenhouse gases and sulphate aerosols: study with a high-resolution AGCM,” Journal of Climate, vol. 16, pp. 2262–2274, 2003.

J. C. Fyfe, “Extratropical Southern Hemisphere cyclones: harbingers of climate change,” Journal of Climate, vol. 16, pp. 2802–2805, 2003.

J. H. Yin, “A consistent poleward shift of the storm tracks in simulations of 21st century climate,” Geophysical Research Letters, vol. 32, no. 18, Article ID L18701, pp. 1–4, 2005.

S. J. Lambert and J. C. Fyfe, “Changes in winter cyclone frequencies and strengths simulated in enhanced greenhouse warming experiments: results from the models participating in the IPCC diagnostic exercise,” Climate Dynamics, vol. 26, no. 7-8, pp. 713–728, 2006.

G. Gastineau and B. J. Soden, “Model projected changes of extreme wind events in response to global warming,” Geophysical Research Letters, vol. 36, no. 10, Article ID L10810, 2009.

Y. Wu, M. Ting, R. Seager, H.-P. Huang, and M. A. Cane, “Changes in storm tracks and energy transports in a warmer climate simulated by the GFDL CM2.1 model,” Climate Dynamics, vol. 37, no. 1, pp. 53–72, 2011.

L. Bengtsson, K. I. Hodges, and E. Roeckner, “Storm tracks and climate change,” Journal of Climate, vol. 19, no. 15, pp. 3518–3543, 2006.

I. Smith, “An assessment of recent trends in Australian rainfall,” Australian Meteorological Magazine, vol. 53, no. 3, pp. 163–173, 2004.

N. Nicholls, Detecting, Understanding and Attributing Climate Change, Australian Greenhouse Office Publication, 2007.

B. C. Bates, P. Hope, B. Ryan, I. Smith, and S. Charles, “Key findings from the Indian Ocean Climate Initiative and their impact on policy development in Australia,” Climatic Change, vol. 89, no. 3-4, pp. 339–354, 2008.

J. S. Frederiksen and C. S. Frederiksen, “Decadal changes in Southern Hemisphere winter cyclogenesis,” CSIRO Marine and Atmospheric Research Paper 002, 2005.

J. S. Frederiksen and C. S. Frederiksen, “Interdecadal changes in southern hemisphere winter storm track modes,” Tellus A, vol. 59, no. 5, pp. 599–617, 2007.

J. S. Frederiksen, C. S. Frederiksen, and S. L. Osbrough, “Modelling of changes in Southern Hemisphere weather systems during the 20th century,” in Proceedings of the 18th World IMACS Congress and International Congress on Modelling and Simulation, pp. 2562–2568, Cairns, Australia, July 2009.

J. S. Frederiksen, C. S. Frederiksen, S. L. Osbrough, and J. M. Sisson, “Causes of changing Southern Hemispheric weather systems,” in Managing Climate Change, I. Jupp, P. Holper, and W. Cai, Eds., pp. 85–98, CSIRO Publishing, 2010.

J. S. Frederiksen, C. S. Frederiksen, S. L. Osbrough, and J. M. Sisson, “Changes in Southern Hemisphere rainfall, circulation and weather systems,” in Proceedings of the 19th International Congress on Modelling and Simulation (MODSIM '11), pp. 2712–2718, Perth, Australia, December 2011.

C. S. Frederiksen, J. S. Frederiksen, J. M. Sisson, and S. L. Osbrough, “Observed and projected changes in the annual cycle of Southern Hemisphere baroclinicity for storm formation,” in Proceedings of the 19th International Congress on Modelling and Simulation (MODSIM' 11), pp. 2719–2725, Perth, Australia, December 2011.

N. A. Phillips, “Geostrophic motion,” Reviews of Geophysics, vol. 1, no. 2, pp. 123–176, 1963.

I. Bordi, K. Fraedrich, F. Lunkeit, and A. Sutera, “Tropospheric double jets, meridional cells, and eddies: a case study and idealized simulations,” Monthly Weather Review, vol. 135, no. 9, pp. 3118–3133, 2007.

I. Bordi, K. Fraedrich, M. Ghil, and A. Sutera, “Zonal flow regime changes in a GCM and in a simple quasigeostrophic model: the role of stratospheric dynamics,” Journal of the Atmospheric Sciences, vol. 66, no. 5, pp. 1366–1383, 2009.

J. G. Charney, “On the scale of atmospheric motions,” Geofysiske Publikasjoner, vol. 17, no. 2, pp. 3–17, 1948.

A. Burger, “Scale consideration of planetary motions of the atmosphere,” Tellus, vol. 10, pp. 195–205, 1958.

I. A. Pisnichenko, “Ultralong-wave dynamics in a two-dimensional baroclinic atmosphere model,” Atmospheric and Oceanic Physics, vol. 16, no. 9, pp. 883–892, 1980.

I. A. Pisnichenko, “Influence variable static stability on the dynamics of ultralong waves in two-dimensional baroclinic model of the atmosphere,” Atmospheric and Oceanic Physics, vol. 19, no. 11, pp. 1223–1226, 1983.

A. Wiin-Nielsen, “On barotropic and baroclinic models, with special emphasis on ultra-long waves,” Monthly Weather Review, vol. 87, pp. 171–183, 1959.

B. W. Bin Wang and A. Barcilon, “Moist stability of a baroclinic zonal flow with conditionally unstable stratification,” Journal of the Atmospheric Sciences, vol. 43, no. 7, pp. 705–719, 1986.

N. A. Phillips, “A coordinate system having some special advantages for numerical forecasting,” Journal of Meteorology, vol. 14, pp. 184–185, 1957.

D. M. Alishaev, “Dynamics of a two-dimensional baroclinic atmosphere,” Atmospheric and Oceanic Physics, vol. 16, no. 2, pp. 99–107, 1980.