# Twentieth Century Winter Changes in Southern Hemisphere Synoptic Weather Modes

Authors: Jorgen S. Frederiksen; Carsten S. Frederiksen
Journal: Advances in Meteorology
 (2011)

---

## Abstract

During the last sixty years, there have been large changes in the southern hemisphere winter circulation and reductions in rainfall particularly in the southern Australian region. Here we examine the corresponding changes in dynamical modes of variability ranging from storm tracks, onset-of-blocking modes, northwest cloud-band disturbances, Antarctic low-frequency modes, intraseasonal oscillations, and African easterly waves. Our study is performed using a global two-level primitive equation instability-model with reanalyzed observed July three-dimensional basic states for the periods 1949–1968, 1975–1994, and 1997–2006. We relate the reduction in the winter rainfall in the southwest of Western Australia since the mid-1970s and in south-eastern Australia since the mid-1990s to changes in growth rate and structures of leading storm track and blocking modes. We find that cyclogenesis and onset-of-blocking modes growing on the subtropical jet have significantly reduced growth rates in the latter periods. On the other hand there is a significant increase in the growth rate of northwest cloud-band modes and intraseasonal oscillation disturbances that cross Australia and are shown to be related to recent positive trends in winter rainfall over northwest Western Australia and central Australia, in general. The implications of our findings are discussed.


---

## Body

## 1. Introduction

Major shifts in the structure of the large-scale circulation of both the northern and southern hemispheres have occurred during the last sixty years or so (see [1, 2], for overview). A number of studies [3–16], have also shown that, since the early to mid-1970s, there has been a dramatic and continuing reduction, as well as negative trend, in the winter rainfall over southern Australia with very large reductions (approximately 20%) occurring first in the southwest of Western Australia (SWWA). More recently, during the last decade or so, there have also been large rainfall reductions in south-eastern Australia [8–11, 17–21].As reviewed by Frederiksen [22], normal mode instability theory has been widely used to understand the dynamics of weather systems [23] and their predictability [24]. Of course nonlinear effects will cause modifications of the disturbances as they increase [25]. Nevertheless, normal modes of climatological basic states and finite-time normal modes of time evolving basic states capture many of the important properties of atmospheric disturbances and their predictability [22].Recently, J. S. Frederiksen and C. S. Frederiksen [1, 2] investigated whether the impact of the dramatic changes in the mean July climate of the atmospheric circulation and thermal structure, for the periods 1949–1968 and 1975–1994, on changes in the southern hemisphere (SH) storm track modes could explain the dramatic decrease in rainfall between these periods. A significant difference between the two periods was a reduction of 17% in the peak strength of the SH subtropical jetstream over the southern Indian Ocean upstream of Australia. They found also that there was a hemisphere-wide reduction in the upper-level zonal winds and vertical wind shear near 30S, and increases in a hemispheric band near 50S. There was also a significant warming of the SH atmosphere south of 30S and a reduced equator-to-pole temperature gradient particularly in the eastern hemisphere. These changes were shown to be associated with about a 33% reduction in growth rates of leading July SH storm track modes, or weather systems, crossing Australia, resulting in a reduction in the intensity of storm formation and the southward deflection of some storms. This has been the primary cause of the rainfall reduction over southwest Western Australia and southern Australia, in general, since 1975.Frederiksen et al. [9] used the Phillips criterion [26] as a diagnostic for evaluating the ability of climate models to capture these changes in the growth rates of storm track modes and the preferred location of storm development. Over the last sixty years, there have been large reductions in the Phillips criterion, or baroclinic instability, near 30S in a hemispheric band consistent with the reductions in the growth rate of storm development in the Australian region. Frederiksen et al. [10, 11] showed that there was a close relationship between the reduction in the Phillips criterion and reductions in winter rainfall during the twentieth century. They also showed, using climate models that were able to reproduce observed changes in the Phillips criterion, that the impact of further increases in anthropogenic CO2 concentrations can lead to further large reductions and significant negative trends in baroclinic instability, by the end of the 21st century, resulting in further reductions in the growth rates of the storm track modes.In this paper, we extend the studies of J. S. Frederiksen and C. S. Frederiksen [1, 2] and Frederiksen et al. [9] and examine other important dynamical modes of variability including onset-of-blocking modes, northwest cloud-band disturbances, Antarctic low-frequency modes, intraseasonal oscillations, and African easterly waves. As in these studies, we look at the leading eigenmodes of the linearized atmospheric equations of motion using the primitive equation model of J. S. Frederiksen and C. S. Frederiksen [2, 27]. We focus on the month of July for three periods 1949–1968, 1975–1994, and 1997–2006, with the aim to see whether the properties of the leading dynamical modes can explain some of the observed changes and trends in Australian rainfall (Figures 1(a) and 1(b)) derived from the Australian Bureau of Meteorological National Climate Centre mean gridded rainfall dataset [28]. Figure 1(a) shows the difference in rainfall between the periods 1975–1994 and 1949–1968 with large reductions over SWWA and eastern and south-eastern Australia of around 20%, 30%, and 10%, respectively. The linear trend in rainfall, over the period 1950–1999, is shown in Figure 1(b). Consistent with Figure 1(a), this shows substantial negative trends over SWWA and eastern and south-eastern Australia and substantial positive trends over northwest Western Australia and over parts of central Australia.Observed July Australian rainfall differences (1975–94) minus (1949–68) (mm/month) (a) and 1950 to 1999 trend in July Australian rainfall (mm/month/year) (b).
(a)(b)The plan of this work is as follows. In Section2, we describe briefly the details of the instability model and some of the climatological fields that constitute the basic states used in the analysis. In Section 3, we discuss the properties of the instability modes in each of the three periods and discuss some of the more significant changes. In Sections 4 to 8, we discuss the five categories of SH dynamical modes that have changed most dramatically over the last sixty or so years, and in Section 9 we examine the sensitivity of our results to the choice of dissipation. Our conclusions are in Section 10.


## 2. Model Details and Basic States

The transient instabilities analyzed in this study have been obtained using the two-level linearized primitive equation model that was developed and applied in a series of studies by J. S. Frederiksen and C. S. Frederiksen [27, 29, 30] and by Frederiksen [31]. The current version of the model includes a generalized Kuo-type heating parameterization that incorporates closures for both convection and evaporation-wind feedback, as outlined in our appendix. We use the values of the model parameters specified in Appendix A that are the same as the standard values used by J. S. Frederiksen and C. S. Frederiksen [1, 2] except in Section 9 where the effect of the choice of dissipation is examined.For the results presented in this paper, the source of the observed fields is the National Centres for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis project [32]. Detailed comparisons of studies of transient modes with NCEP/NCAR and ERA40 basic states were made in J. S. Frederiksen and C. S. Frederiksen [2], and it was found that for the same periods of time the structures and growth rates of the modes considered were generally very similar.We focus on the global Southern Hemisphere (SH) mean July climate fields for the periods 1949–68 and 1975–94, as in J. S. Frederiksen and C. S. Frederiksen [2], and also consider the more recent period 1997–2006. We recall that between the periods 1948–68 and 1975–94 there is a very significant  reduction of 17% in the peak strength of the SH subtropical jet stream (Figure 1 of [2]). A similar reduction occuring between 1948–68 and 1997–2006 is seen in Figure 2(a) which shows the height (in pressure units) and latitudinal cross-section of the zonal wind, averaged over 110°E–160°E longitude. The latitudinal and longitudinal changes between the latter periods at 300 hPa are also shown in Figure 2(b). There we note the reduction in the subtropical jet strength near 30°S is largely zonally symmetric and is accompanied by a corresponding increase in the zonal wind near 50°S. These changes also extend to lower levels with reduced magnitudes as shown for 700 hPa in Figure 2(c).Difference between (1997–2006) and (1949–68) basic states of vertical cross-section of July zonal wind (ms−1) averaged between 110° and 160°E as a function of latitude and pressure (in hPa) (a), 300 hPa July zonal wind (ms−1) (b) and 700 hPa July zonal wind (ms−1) (c).
(a)(b)(c)The changes in the zonal wind shown in Figure2 are associated with changes in the thermal structure of the SH atmosphere. By 1997–2006 there has been a significant warming south of 30°S, as shown by the 500 hPa temperature difference between 1997–2006 and 1949–69 in Figure 3(a). This has                                                                                                                        reduced the equator-to-high-latitude temperature gradient across the whole Southern Hemisphere. More over, the difference in the 500 hPa geopotential height between the two periods (Figure 3(b)) is associated with a change in the strength of the Antarctic Oscillation or southern hemisphere annular mode (SAM), [33–36]. The changes in SAM extend to the surface pressure as shown in Figure 3(c). These large-scale global changes would be expected to have a significant effect on the stability of the Southern Hemisphere circulation and hence on the nature of the Southern Hemisphere storms and other modes of weather variability.Difference between (1997–2006) and (1949–68) basic states of vertical cross-section of July 500 hPa temperature (K) (a), 500 hPa geopotential height (m) (b), and sea-level pressure (hPa) (c).
(a)(b)(c)The flow fields for the basic states that we focus on in this paper are the monthly averaged July global fields for the periods 1949–68, 1975–94, and 1997–2006 from the NCEP/NCAR reanalyses. We use global 300 hPa and 700-hPa fields as representative of the upper and lower levels, respectively and the fields of streamfunction, velocity potential, and potential temperature are represented as rhomboidal wavenumberR=15 spectral fields as described in J. S. Frederiksen and C. S. Frederiksen [29].


## 3. Calculation of Weather Modes

The synoptic weather modes and lower-frequency disturbances for the 1949–68, 1975–94, and 1997–2006 basic states described in Section2 have been calculated using the methods described in Appendix A. As in our previous studies, a wide variety of modal disturbances that can be related to many aspects of atmospheric dynamics have been obtained. Here we consider some of the faster growing modes that show significant differences in growth rates between the three time periods. As noted in J. S. Frederiksen and C. S. Frederiksen [2], Southern Hemisphere July storm track modes have significantly smaller growth rates during 1975–94 compared with 1949–68. We find that growth rates of leading Southern Hemisphere July storm track modes are also reduced in 1997–2006 compared with 1949–68. A major focus of this paper is to examine growth rate changes in onset-of-blocking modes, Antarctic modes, northwest cloud band modes, and intraseasonal oscillation modes between these three periods. Here the modes are ordered by decreasing growth rate, and the mode number specifies that order.Table1 gives the mode number, period, and growth rates of modes of particular interest for the 1949–68 and 1975–94 basic states. Also shown are the mode types and the largest pattern correlations between corresponding modes for the two basic states. Corresponding results for the 1949–68 and 1997–2006 basic states are shown in Table 2. The modes have been ordered in terms of classes with generally increasing periods followed by the purely tropical African easterly wave modes. There are additional modes whose properties are not summarized in the Tables. These modes include northern hemisphere modes, which will not be considered here, and subdominant southern hemisphere modes, many of which are storm track modes, with properties similar to the leading modes or that do not affect the Australian region. In these two tables, the pattern correlations shown are the maximum over the phases of the modes and the averages taken over the five model fields ψ, τ,χ,θ, and σ. The pattern correlations of the perturbation fields are calculated as described in Appendix C of Frederiksen and Bell [37]. For a good fit, the pattern correlations as calculated require that both the real and imaginary parts of the two modes of interest are similar; this ensures that two travelling modes are similar throughout their evolution. It is, however, a more stringent criterion than just calculating pattern correlations between the same type of fields at just one corresponding phase, and therefore our correlations may tend to be smaller.Table 1
PeriodsTrd and growth rates ωid of selected modes for the 1949–68 NCEP/NCAR basic state with standard parameterizations (K′d=2.338×1016 m4 s−1, QFd=1500K, CFd=5×10-4K m−1), their mode number, and class of mode.  Also shown are the largest pattern correlations Ac between a given mode and the modes for the 1975–94 NCEP/NCAR basic state with standard parameterizations and the 1975–94 mode for which the pattern correlation is largest as well as the change in the growth rate Δωid.
ClassModeTrdωidTrdωidAcModeΔωid1949–68(days)(day−1)(days)(day−1)1975–94SH cyclogenesis11.200.42501.300.28240.91489−33.5%ML blocking306.110.25666.250.22790.459721−11.2%HL blocking87.400.34958.020.33170.83424−5.1%HL blocking115.220.31845.080.31010.84686−2.6%NWCB438.010.21008.190.26240.803612+24.9%Antarctic613.080.362417.480.36060.94163−0.5%Antarctic2652.10.2746121.30.27000.854211−1.7%ISO5759.00.144441.850.18580.850635+28.7%African EW224.140.28813.830.30710.79597+6.6%Table 2
PeriodsTrd and growth rates ωid of selected modes for the 1949–68 NCEP/NCAR basic state with standard parameterizations (K′d=2.338×1016m4 s−1, QFd=1500K, CFd=5×10-4K m−1), their mode number, and class of mode.  Also shown are the largest pattern correlations Ac between a given mode and the modes for the 1997–2006 NCEP/NCAR basic state with standard parameterizations and the 1997–2006 mode for which the pattern correlation is largest as well as the change in the growth rate Δωid.
ClassModeTrdωidTrdωidAcModeΔωid1949–68(days)(day−1)(days)(day−1)1997–2006SH cyclogenesis11.200.42501.330.26630.896012−37.3%ML blocking306.110.25666.270.20740.639926−19.2%HL blocking87.400.34958.010.35000.76445+0.1%HL blocking115.220.31844.990.30600.80286−3.9%NWCB438.010.21007.970.30540.63747+45.4%Antarctic613.080.362420.00.38100.90713+4.9%Antarctic2652.10.2746∞0.28500.647510+3.8%ISO5759.00.144425.40.18860.574433+30.6%African EW224.140.28813.860.30350.76388+5.3%The leading mode for the 1949–68 basic state is a SH cyclogenesis mode with large pattern correlation (approximately 0.9) with mode 9 for the 1975–94 basic state (Table1) and with mode 12 for the 1997–2006 basic state (Table 2).   For the later periods the growth rates of the SH cyclogenesis modes are reduced by 33% (Table 1) and 37% (Table 2), respectively. Thus the growth rate reduction has continued to the most recent period.  As noted in J. S. Frederiksen and C. S. Frederiksen [2], on average the 10 leading Southern Hemisphere cyclogenesis modes for 1975–94, most of which cross Australia, have growth rates that are 32% smaller than for the corresponding modes for 1949–68. This property of leading SH cyclogenesis modes crossing Australia to have reduced growth rates of approximately 30% continues to the most recent period of 1997–2006.From Tables1 and 2 we note that there is little systematic change in the growth rates of leading Antarctic modes in the latter periods 1975–94 and 1997–2006 compared with the early period. For the midlatitude SH onset-of-blocking modes, there is considerable reduction in growth rates in the latter periods, and pattern correlations between mode 30 for 1949–68 and those for the latter periods are less than the 0.7578 (not shown) between mode 21 for 1975–94 and mode 26 for 1997–2006. The growth rates of northwest cloud-band and intraseasonal oscillation modes increase considerably in the latter time intervals. For the modes with both high latitude and midlatitude multipole structures the changes in growth rates are small while the African easterly waves show moderate increases in growth rate during the latter time periods. All the modes considered are eastward propagating apart from the African easterly waves.In the following sections, we examine the structures of the main classes of Southern Hemisphere modes in Tables1 and 2 during the three periods of interest and relate them, in particular, to the observed changes in Australian rainfall (Figure 1).


## 4. Southern Hemisphere Storm Track Modes

### 4.1. Growth Rate and Periods

The leading Southern Hemisphere storm track modes [38, 39] for the time intervals 1949–68, 1975–94, and 1997–2006 are modes 1, 9, and 12. As shown in Tables 1 and 2, they all have periods of 1.3 days but growth rates of 0.423 day−1, 0.282 day−1, and 0.266 day−1 and pattern correlations with mode 1 for 1949–68 of 1.0, 0.915, and 0.896, respectively. This corresponds to reductions in growth rates of 33.5 per cent and 37.3 per cent during 1975–94 and 1997–2006, respectively.


### 4.2. Structures of Storm Track Modes

Figure4 shows the upper level streamfunction at a particular instance, its random-phase ensemble average (RPEA), the upper level divergence, and its RPEA for mode 12 for the 1997–2006 basic state, which is the fastest growing storm track mode crossing Australia. The RPEA of these fields is defined as in Equation (2.7) of Frederiksen [40]. The RPEA of a disturbance field is essentially the amplitude envelope within which the disturbance propagates with changing time or phase and determines in this case the structure of the storm track. Mode 12 for the 1997–2006 basic state is very similar to mode 1 for 1949–68 and mode 9 for 1975–94, as expected from their high pattern correlations. Like mode 1 for 1949–68 and mode 9 for 1975–94, the disturbance has a wavenumber of around 12 in the regions of large amplitude over southern Australia and consists of a series of eastward-propagating troughs (cool colours) and ridges (warm colours). As the troughs and ridges move eastward they amplify to reach a maximum in preferred regions of large RPEA amplitude. In Figures 4(b) and 4(d), the amplitude variations of the storm mode are more evident with both the streamfunction and divergence fields having largest amplitude focused over south-western Australia and secondary maxima over the central Southern Pacific Ocean. We note that precipitation is essentially determined by the lower level vertical velocity in pressure coordinates, and in our model this is proportional to the magnitude of the upper level divergence. Thus, the pattern of RPEA divergence in Figure 4(c) is also the expected pattern of rainfall associated with this mode, and as we have noted above, the maximum is focused over southern Australia. Thus the significant reduction in growth rates of storm track modes in the latter two time intervals (1975–94 and 1997–2006) has been consistent with the reduction in rainfall (Figure 1) over southern Australia and, in particular southwest Western Australia, since the 1970s.Shown in arbitrary units are the upper-level streamfunction (a), the RPEA of the upper-level streamfunction (b), the upper-level divergence (c), and the RPEA of the upper-level divergence (d) for mode 12, the fastest growing SH storm track mode for the July 1997–2006 basic state.
(a)(b)(c)(d)Mode 12, shown in Figure4, is the fastest growing of a class of cyclogenesis modes with zonal wavenumbers between 6 and 12 and large amplitudes in the Australian region.


## 4.1. Growth Rate and Periods

The leading Southern Hemisphere storm track modes [38, 39] for the time intervals 1949–68, 1975–94, and 1997–2006 are modes 1, 9, and 12. As shown in Tables 1 and 2, they all have periods of 1.3 days but growth rates of 0.423 day−1, 0.282 day−1, and 0.266 day−1 and pattern correlations with mode 1 for 1949–68 of 1.0, 0.915, and 0.896, respectively. This corresponds to reductions in growth rates of 33.5 per cent and 37.3 per cent during 1975–94 and 1997–2006, respectively.


## 4.2. Structures of Storm Track Modes

Figure4 shows the upper level streamfunction at a particular instance, its random-phase ensemble average (RPEA), the upper level divergence, and its RPEA for mode 12 for the 1997–2006 basic state, which is the fastest growing storm track mode crossing Australia. The RPEA of these fields is defined as in Equation (2.7) of Frederiksen [40]. The RPEA of a disturbance field is essentially the amplitude envelope within which the disturbance propagates with changing time or phase and determines in this case the structure of the storm track. Mode 12 for the 1997–2006 basic state is very similar to mode 1 for 1949–68 and mode 9 for 1975–94, as expected from their high pattern correlations. Like mode 1 for 1949–68 and mode 9 for 1975–94, the disturbance has a wavenumber of around 12 in the regions of large amplitude over southern Australia and consists of a series of eastward-propagating troughs (cool colours) and ridges (warm colours). As the troughs and ridges move eastward they amplify to reach a maximum in preferred regions of large RPEA amplitude. In Figures 4(b) and 4(d), the amplitude variations of the storm mode are more evident with both the streamfunction and divergence fields having largest amplitude focused over south-western Australia and secondary maxima over the central Southern Pacific Ocean. We note that precipitation is essentially determined by the lower level vertical velocity in pressure coordinates, and in our model this is proportional to the magnitude of the upper level divergence. Thus, the pattern of RPEA divergence in Figure 4(c) is also the expected pattern of rainfall associated with this mode, and as we have noted above, the maximum is focused over southern Australia. Thus the significant reduction in growth rates of storm track modes in the latter two time intervals (1975–94 and 1997–2006) has been consistent with the reduction in rainfall (Figure 1) over southern Australia and, in particular southwest Western Australia, since the 1970s.Shown in arbitrary units are the upper-level streamfunction (a), the RPEA of the upper-level streamfunction (b), the upper-level divergence (c), and the RPEA of the upper-level divergence (d) for mode 12, the fastest growing SH storm track mode for the July 1997–2006 basic state.
(a)(b)(c)(d)Mode 12, shown in Figure4, is the fastest growing of a class of cyclogenesis modes with zonal wavenumbers between 6 and 12 and large amplitudes in the Australian region.


## 5. Midlatitude Onset-of-Blocking Modes

### 5.1. Growth Rates and Periods

For the time intervals 1949–68, 1975–94, and 1997–2006, the leading southern hemisphere midlatitude onset-of-blocking modes [23, 30] are modes 30, 21, and 26 as shown in Tables 1 and 2. All have periods of 6 days but growth rates of 0.257 day−1, 0.228 day−1, and 0.207 day−1 and pattern correlations with mode 30 for 1949–68 of 1.0, 0.460, and 0.640, respectively. Thus there are reductions in growth rates of 11.2 per cent and 19.2 per cent in 1975–94 and 1997–2006 compared with 1949–68.


### 5.2. Structures of Midlatitude Onset-of-Blocking Modes

Midlatitude dipole onset-of-blocking modes are eastward propagating, have longer periods than cyclogenesis modes (approximately 6 days in Tables1 and 2), are of somewhat larger scale, and have dipole structures in the latitudinal direction in the major regions of observed Southern Hemisphere blocking. Figure 5 shows an example of a mid-latitude onset-of-blocking mode (mode 26) for the 1997–2006 basic state; the depicted upper-level streamfunction has dipole structures between Australia and New Zealand, extending downstream over the central Pacific Ocean and there are dipole structures over the southern Atlantic Ocean as well. We note from Tables 1 and 2 that there are somewhat similar onset-of-blocking modes for the periods 1949–68 and 1975–94 but that the different basic states cause changes in the structures since the pattern correlations are at best moderate. The primary differences in structures in the midlatitude onset-of-blocking modes for the three periods relate to differences in the longitude of maximum amplitude of upper-level RPEA streamfunction (not shown) near 30S.  For mode 26 for the 1997–2006 basic state, this maximum occurs near 120W; for mode 21 for the 1975–94 basic state, it occurs near 135W; and for mode 30 for the 1949–68 basic state it occurs near 170W with relatively larger amplitude blocks in the region of New Zealand. We also note that the growth rates in 1997–2006 and 1975–94 are approximately 19% and approximately 11% less than for 1949–68, respectively. We expect that the growth rate reductions in 1975–94 and 1997–2006 compared with 1949–68 at least partly reflect the reductions in baroclinicity of the subtropical jet that cause even greater reductions in the growth rate of storm track modes. The blocking highs in the region of New Zealand are also frequently associated with rain-bearing east coast lows on the east coast of Australia. Again, the changes and trends in rainfall over Australia in Figure 1 show reductions and negative trends in rainfall along the east coast of Australia consistent with the smaller average growth rates of these larger-scale modes since the 1970s.Figure 5
As in Figure4(a) but for the SH midlatitude onset-of-blocking mode, mode 26, for the July 1997–2006 basic state.


## 5.1. Growth Rates and Periods

For the time intervals 1949–68, 1975–94, and 1997–2006, the leading southern hemisphere midlatitude onset-of-blocking modes [23, 30] are modes 30, 21, and 26 as shown in Tables 1 and 2. All have periods of 6 days but growth rates of 0.257 day−1, 0.228 day−1, and 0.207 day−1 and pattern correlations with mode 30 for 1949–68 of 1.0, 0.460, and 0.640, respectively. Thus there are reductions in growth rates of 11.2 per cent and 19.2 per cent in 1975–94 and 1997–2006 compared with 1949–68.


## 5.2. Structures of Midlatitude Onset-of-Blocking Modes

Midlatitude dipole onset-of-blocking modes are eastward propagating, have longer periods than cyclogenesis modes (approximately 6 days in Tables1 and 2), are of somewhat larger scale, and have dipole structures in the latitudinal direction in the major regions of observed Southern Hemisphere blocking. Figure 5 shows an example of a mid-latitude onset-of-blocking mode (mode 26) for the 1997–2006 basic state; the depicted upper-level streamfunction has dipole structures between Australia and New Zealand, extending downstream over the central Pacific Ocean and there are dipole structures over the southern Atlantic Ocean as well. We note from Tables 1 and 2 that there are somewhat similar onset-of-blocking modes for the periods 1949–68 and 1975–94 but that the different basic states cause changes in the structures since the pattern correlations are at best moderate. The primary differences in structures in the midlatitude onset-of-blocking modes for the three periods relate to differences in the longitude of maximum amplitude of upper-level RPEA streamfunction (not shown) near 30S.  For mode 26 for the 1997–2006 basic state, this maximum occurs near 120W; for mode 21 for the 1975–94 basic state, it occurs near 135W; and for mode 30 for the 1949–68 basic state it occurs near 170W with relatively larger amplitude blocks in the region of New Zealand. We also note that the growth rates in 1997–2006 and 1975–94 are approximately 19% and approximately 11% less than for 1949–68, respectively. We expect that the growth rate reductions in 1975–94 and 1997–2006 compared with 1949–68 at least partly reflect the reductions in baroclinicity of the subtropical jet that cause even greater reductions in the growth rate of storm track modes. The blocking highs in the region of New Zealand are also frequently associated with rain-bearing east coast lows on the east coast of Australia. Again, the changes and trends in rainfall over Australia in Figure 1 show reductions and negative trends in rainfall along the east coast of Australia consistent with the smaller average growth rates of these larger-scale modes since the 1970s.Figure 5
As in Figure4(a) but for the SH midlatitude onset-of-blocking mode, mode 26, for the July 1997–2006 basic state.


## 6. Northwest Cloud-Band Modes

### 6.1. Growth Rates and Periods

In Figure1 we also see that there are regions of increased rainfall stretching southeast from northwest and central Western Australia. These are regions often affected by propagating cloud-bands originating off the northwest coast of Australia. Here we examine the extent to which these results may be related to changes in northwest cloud-band (NWCB) disturbances [41]. The leading northwest cloud band modes for the time intervals 1949–68, 1975–94, and 1997–2006 are modes 43, 12, and 7, and all have periods of 8 days but growth rates of 0.210 day−1, 0.262 day−1, and 0.305 day−1 and pattern correlations with mode 43 for 1949–68 of 1.0, 0.840, and 0.637, respectively. This corresponds to increases in growth rates of 24.9 per cent and 45.4 per cent during 1975–94 and 1997–2006, respectively.


### 6.2. Structures of Northwest Cloud-Band Modes

The northwest cloud-band modes in Tables1 and 2 are eastward propagating with periods of approximately 8 days although there are also slower-growing NWCB modes with longer periods (as discussed in a different context by C. S. Fredriksen and J. S. Fredriksen [41]). Mode 12 for the 1975–94 basic state is a NWCB mode that emanates from the region of the Indian Ocean, crosses Australia, and then propagates into the central Pacific and south across South America. It again is of larger scale than the storm track modes such as that shown in Figure 4. It also differs significantly from the storm track and onset-of-blocking modes in that its growth rate is approximately 25% greater than the corresponding mode for 1949–68. This increases to approximately 45% for NWCB mode 7 for the 1997–2006 basic state compared to the corresponding NWCB mode for the 1949–68 basic state.Figure6 shows some of the features of mode 12 for the 1975–94 basic state. The upper-level disturbance streamfunction is shown in Figure 6(a) with a characteristic wave train emanating from the Indian Ocean off the northwest coast. The RPEA of this field for this mode is shown in Figure 6(b) and illustrates the large impact (i.e., maxima in amplitude) of the mode over the Indian Ocean and central Australia. The upper-level divergence field depicted in Figure 6(c) shows large amplitude over central Western Australia and northern Queensland, where it would be expected to have a large impact on the rainfall. The fact that the growth rate continues to increase compared with the earliest period is consistent with the increased tendency for these modes to develop and consequently for the increase in rainfall shown over northwest Western Australia and northern Queensland (Figure 1).As in Figure4 but for the SH upper-level streamfunction (a), the RPEA of the upper-level streamfunction (b), and the upper-level divergence in the Australian region (c) for the leading NWCB mode, mode 12, for the July 1975–94 basic state.
(a)(b)(c)


## 6.1. Growth Rates and Periods

In Figure1 we also see that there are regions of increased rainfall stretching southeast from northwest and central Western Australia. These are regions often affected by propagating cloud-bands originating off the northwest coast of Australia. Here we examine the extent to which these results may be related to changes in northwest cloud-band (NWCB) disturbances [41]. The leading northwest cloud band modes for the time intervals 1949–68, 1975–94, and 1997–2006 are modes 43, 12, and 7, and all have periods of 8 days but growth rates of 0.210 day−1, 0.262 day−1, and 0.305 day−1 and pattern correlations with mode 43 for 1949–68 of 1.0, 0.840, and 0.637, respectively. This corresponds to increases in growth rates of 24.9 per cent and 45.4 per cent during 1975–94 and 1997–2006, respectively.


## 6.2. Structures of Northwest Cloud-Band Modes

The northwest cloud-band modes in Tables1 and 2 are eastward propagating with periods of approximately 8 days although there are also slower-growing NWCB modes with longer periods (as discussed in a different context by C. S. Fredriksen and J. S. Fredriksen [41]). Mode 12 for the 1975–94 basic state is a NWCB mode that emanates from the region of the Indian Ocean, crosses Australia, and then propagates into the central Pacific and south across South America. It again is of larger scale than the storm track modes such as that shown in Figure 4. It also differs significantly from the storm track and onset-of-blocking modes in that its growth rate is approximately 25% greater than the corresponding mode for 1949–68. This increases to approximately 45% for NWCB mode 7 for the 1997–2006 basic state compared to the corresponding NWCB mode for the 1949–68 basic state.Figure6 shows some of the features of mode 12 for the 1975–94 basic state. The upper-level disturbance streamfunction is shown in Figure 6(a) with a characteristic wave train emanating from the Indian Ocean off the northwest coast. The RPEA of this field for this mode is shown in Figure 6(b) and illustrates the large impact (i.e., maxima in amplitude) of the mode over the Indian Ocean and central Australia. The upper-level divergence field depicted in Figure 6(c) shows large amplitude over central Western Australia and northern Queensland, where it would be expected to have a large impact on the rainfall. The fact that the growth rate continues to increase compared with the earliest period is consistent with the increased tendency for these modes to develop and consequently for the increase in rainfall shown over northwest Western Australia and northern Queensland (Figure 1).As in Figure4 but for the SH upper-level streamfunction (a), the RPEA of the upper-level streamfunction (b), and the upper-level divergence in the Australian region (c) for the leading NWCB mode, mode 12, for the July 1975–94 basic state.
(a)(b)(c)


## 7. Intraseasonal Oscillation Modes

### 7.1. Growth Rates and Periods

Next, we consider the leading intraseasonal oscillation modes [29, 31] for the time intervals 1949–68, 1975–94, and 1997–2006. As seen from Tables 1 and 2, they are modes 57, 35, and 33 and have periods of 59 days, 35 days, and 25 days, growth rates of (0.144 day−1, 0.186 day−1, 0.189 day−1) and pattern correlations with mode 57 for 1949–68 of (1.0, 0.851, 0.574) respectively. This corresponds to increases in growth rates of 28.7% and 30.6% during 1975–94 and 1997–2006, respectively.


### 7.2. Structures of Intraseasonal Modes

The July intraseasonal oscillation modes for basic states during the second half of the twentieth century have streamfunctions that extend from the equator into both hemispheres. In the Southern Hemisphere, eastward propagating mode 35 for the 1975–94 basic state has a wave train extending from the Indian Ocean across Australia into the central and south Pacific as shown in Figure7(a) for the 300 hPa disturbance streamfunction. The 300 hPa divergence field at the phase depicted in Figure 7(b) has large amplitude over central Western Australia. It is interesting to note how the period and growth rate of this mode change between the time intervals 1949–68, 1975–94, and 1997–2006. We see from Tables 1 and 2 that the periods of this intraseasonal oscillation mode are reduced from 57 days to 42 days to 25 days. At the same time the growth rate increases by approximately 29% between 1949–68 and 1975–94 and by approximately 31% between 1949–68 and 1997–2006. We also note that while the correlation between mode 33 for 1997–2006 and mode 57 for 1949–68 is 0.5744, it is considerably larger with mode 35 for 1975–94 at 0.8030 (not shown). Thus there is a systematic change in the properties of the intraseasonal oscillation modes with time.As in Figure4 but for the global upper-level streamfunction (a) and the upper-level divergence in the Australian region (b) for the leading intraseasonal oscillation mode, mode 35, for the July 1975–94 basic state.
(a)(b)The 300 hPa divergence field depicted in Figure7(b) has large amplitude over central Western Australia, much like the NWCB mode (Figure 6(c)). Again, the increase in the growth rate of this mode is consistent with the increased development of these systems and hence the increase in rainfall over this part of Western Australia.


## 7.1. Growth Rates and Periods

Next, we consider the leading intraseasonal oscillation modes [29, 31] for the time intervals 1949–68, 1975–94, and 1997–2006. As seen from Tables 1 and 2, they are modes 57, 35, and 33 and have periods of 59 days, 35 days, and 25 days, growth rates of (0.144 day−1, 0.186 day−1, 0.189 day−1) and pattern correlations with mode 57 for 1949–68 of (1.0, 0.851, 0.574) respectively. This corresponds to increases in growth rates of 28.7% and 30.6% during 1975–94 and 1997–2006, respectively.


## 7.2. Structures of Intraseasonal Modes

The July intraseasonal oscillation modes for basic states during the second half of the twentieth century have streamfunctions that extend from the equator into both hemispheres. In the Southern Hemisphere, eastward propagating mode 35 for the 1975–94 basic state has a wave train extending from the Indian Ocean across Australia into the central and south Pacific as shown in Figure7(a) for the 300 hPa disturbance streamfunction. The 300 hPa divergence field at the phase depicted in Figure 7(b) has large amplitude over central Western Australia. It is interesting to note how the period and growth rate of this mode change between the time intervals 1949–68, 1975–94, and 1997–2006. We see from Tables 1 and 2 that the periods of this intraseasonal oscillation mode are reduced from 57 days to 42 days to 25 days. At the same time the growth rate increases by approximately 29% between 1949–68 and 1975–94 and by approximately 31% between 1949–68 and 1997–2006. We also note that while the correlation between mode 33 for 1997–2006 and mode 57 for 1949–68 is 0.5744, it is considerably larger with mode 35 for 1975–94 at 0.8030 (not shown). Thus there is a systematic change in the properties of the intraseasonal oscillation modes with time.As in Figure4 but for the global upper-level streamfunction (a) and the upper-level divergence in the Australian region (b) for the leading intraseasonal oscillation mode, mode 35, for the July 1975–94 basic state.
(a)(b)The 300 hPa divergence field depicted in Figure7(b) has large amplitude over central Western Australia, much like the NWCB mode (Figure 6(c)). Again, the increase in the growth rate of this mode is consistent with the increased development of these systems and hence the increase in rainfall over this part of Western Australia.


## 8. Other Dynamically Important Weather Modes

There are groups of other leading modes that also could potentially play important roles in the changing climate during the second half of the 20th century. These include high-latitude blocking modes, Antarctic modes, and African easterly waves.


### 8.1. High-Latitude Blocking Modes

High-latitude blocking modes, such as modes 4 and 6 for 1975–94, have significant amplitudes at high latitudes between New Zealand and Antarctica, as well as in midlatitudes; they have tripole structures in latitude at some longitudes (not shown). In mid-latitudes there are dipole structures downstream of Australia with additional centres around Antarctica similar to those for the mode shown in Figure7 of J. S. Frederiksen and C. S. Frederiksen [30]. For the three time intervals considered throughout this paper, these modes have periods between 5 and 8 days (Tables 1 and 2) and are eastward propagating but show little systematic change in growth rates. Perhaps this is because they also grow on both the polar jet that has strengthened and on the subtropical jet that has weakened with time.


### 8.2. Antarctic Modes

For the 1975–94 basic state the two leading Antarctic modes are mode 3 with a period of approximately 17 days and mode 11 with a period of approximately 52 days. These two lower frequency modes are located principally over the Antarctic regions and consist of superpositions of zonal wavenumbers 2 and 3. They have close similarities with corresponding Antarctic teleconnection pattern modes of Frederiksen and Frederiksen [30, Figure 11]. Again, there are closely similar teleconnection pattern modes with the basic states for 1949–68 and for 1997–2006. These teleconnection patterns are little changed in growth rate or structure during the three periods of interest. This indicates that while there has been a change in the strength of the mean Antarctic oscillation (Section 2), its variability about the mean has changed less.


### 8.3. African Easterly Waves

African easterly waves are westward-propagating tropical waves that are generated by the African easterly jet and originate on the east coast of North Africa. They are frequent precursors of Atlantic hurricanes with about 85% of intense hurricanes having their origins as African easterly waves ([42] and references therein). The upper-level streamfunction for mode 7 for the 1975–94 basic state is shown in Figure 8. It is a baroclinic mode with a period of 4 days and with growth rates that are between 5% and 7% larger after 1975 than for the earlier period of 1949–68 (Tables 1 and 2).Figure 8
As in Figure4(a) but for the African easterly wave, mode 7, for the July 1975–94 basic state.


## 8.1. High-Latitude Blocking Modes

High-latitude blocking modes, such as modes 4 and 6 for 1975–94, have significant amplitudes at high latitudes between New Zealand and Antarctica, as well as in midlatitudes; they have tripole structures in latitude at some longitudes (not shown). In mid-latitudes there are dipole structures downstream of Australia with additional centres around Antarctica similar to those for the mode shown in Figure7 of J. S. Frederiksen and C. S. Frederiksen [30]. For the three time intervals considered throughout this paper, these modes have periods between 5 and 8 days (Tables 1 and 2) and are eastward propagating but show little systematic change in growth rates. Perhaps this is because they also grow on both the polar jet that has strengthened and on the subtropical jet that has weakened with time.


## 8.2. Antarctic Modes

For the 1975–94 basic state the two leading Antarctic modes are mode 3 with a period of approximately 17 days and mode 11 with a period of approximately 52 days. These two lower frequency modes are located principally over the Antarctic regions and consist of superpositions of zonal wavenumbers 2 and 3. They have close similarities with corresponding Antarctic teleconnection pattern modes of Frederiksen and Frederiksen [30, Figure 11]. Again, there are closely similar teleconnection pattern modes with the basic states for 1949–68 and for 1997–2006. These teleconnection patterns are little changed in growth rate or structure during the three periods of interest. This indicates that while there has been a change in the strength of the mean Antarctic oscillation (Section 2), its variability about the mean has changed less.


## 8.3. African Easterly Waves

African easterly waves are westward-propagating tropical waves that are generated by the African easterly jet and originate on the east coast of North Africa. They are frequent precursors of Atlantic hurricanes with about 85% of intense hurricanes having their origins as African easterly waves ([42] and references therein). The upper-level streamfunction for mode 7 for the 1975–94 basic state is shown in Figure 8. It is a baroclinic mode with a period of 4 days and with growth rates that are between 5% and 7% larger after 1975 than for the earlier period of 1949–68 (Tables 1 and 2).Figure 8
As in Figure4(a) but for the African easterly wave, mode 7, for the July 1975–94 basic state.


## 9. Sensitivity Studies

In this study a rhomboidalR=15 truncation is used for both the perturbation and basic-state fields. This resolution, corresponding to (3R+1) by (5R+1)/2 grid points, or circa 750 km by 500 km at 30S in longitude by latitude, has been used in many simulation and instability studies and is adequate for resolving the leading storm track modes and lower-frequency modes of variability. For example, the storm track in Figure 4(b), associated with the smallest scale wavenumber 12 mode considered here, has a longitudinal extent of around 4000 km and a latitudinal extent of circa 2000 km in the Australian region. Studies of the relationship between resolution and the dissipation required to obtain the same large-scale mean and transient flow properties and predictability have been carried out by Frederiksen and Kepert [43] and Zidikheri and Frederiksen [44]. Here we examine the sensitivity of the instability properties for a range of dissipations to draw general conclusions.Our standard dissipationK′d=2.338×1016 m4 s−1 has also been the standard dissipation used in many general circulation models at rhomboidal R=15 resolution going back to the 1970s ([25] and references therein). We expect that this is the most realistic choice, but in this section we also present results for mid-latitude blocking, northwest cloud band and intraseasonal oscillation modes for twice this value and with zero dissipation to gauge the sensitivity of our results. Tables 3 and 4 show the mode number, period Trd, growth rate ωid, and e-folding time τid of these modes with twice and zero dissipation and the largest pattern correlations Ac with the corresponding mode with standard dissipation for the 1949–68 and 1975–94 basic states, respectively. Table 5 shows the change in the growth rate Δωid for these modes between 1949–68 and 1975–94 and how it depends on the strength of the dissipation.Table 3
PeriodsTrd, growth rates ωid, and e-folding times τid, of modes for the 1949–68 basic state with QFd=1500K and CFd=5×10-4K m−1 and for different K′d,  which have the largest pattern correlations Ac with leading ML blocking, NWCB, and ISO modes for the 1949–68 basic state with standard parameterizations. Also shown are the mode numbers.
ClassModeK′dTrdωidτidAcModem4 s−1(days)(day−1)(days)standardML blocking5806.0160.30693.2580.795630ML blocking154.676×10166.1880.22144.5170.937530NWCB4009.8330.38332.6090.741243NWCB454.676×10168.8000.10969.1260.704943ISO95053.880.17775.6280.837557ISO564.676×101661.410.086811.520.915557Table 4
As in Table3 for modes for 1975–94 basic states with different dissipation parameterizations, K′d, that have the largest pattern correlations, Ac, with leading ML blocking, NWCB, and ISO modes for the 1975–94 basic state with standard parameterizations.
ClassModeK′dTrdωidτidAcModem4 s−1(days)(day−1)(days)StandardML blocking5506.2050.28313.5320.770721ML blocking114.676×10166.2180.18625.3710.823721NWCB608.8060.48942.0430.633812NWCB194.676×10167.5140.15796.3340.785012ISO75034.370.22944.3580.865435ISO224.676×101643.280.14926.7010.958435Table 5
Changes in growth rates,Δωid, for leading blocking, NWCB, and ISO modes when basic state is changed from 1949–68 basic state to 1975–94 basic, for different dissipation parameters, K′d.
K′dBasic stateModeBasic stateModeΔωidm4 s−12.338×10161949–6830  ML blocking1975–9421 ML blocking−11.2%01949–6815 ML blocking1975–9411 ML blocking−15.9%4.676×10161949–6858 ML blocking1975–9455 ML blocking−7.7%2.338×10161949–6843 NWCB1975–9412 NWCB+24.9%01949–6840 NWCB1975–946 NWCB+27.7%4.676×10161949–6845 NWCB1975–9419 NWCB+44.1%2.338×10161949–6857 ISO1975–9435 ISO+28.7%01949–6895 ISO1975–9475 ISO+29.1%4.676×10161949–6856 ISO1975–9422 ISO+71.9%We see from Table5 that the changes in Δωid are generally quite similar for standard and zero dissipation and thus that the results should depend only weakly on dissipation in this range. For twice the dissipation, the differences in Δωid increase for NWCB and ISO modes that are close to the instability threshold while Δωid is more similar to that for the standard dissipation for the mid-latitude blocking modes that are further from threshold. This is to be expected. For all our choices of dissipation, the changes in growth rate Δωid for mid-latitude blocking, northwest cloud band, and intraseasonal oscillation modes between 1949–68 and 1975–94 are significant.


## 10. Discussion and Conclusions

In this paper, we have sought to understand some of the changes and trends in Australian winter rainfall over the last fifty or so years in the context of changes in the SH large-scale circulation and the associated properties of relevant SH weather systems and modes of variability. We have examined, in particular, the changes in the structures and growth rates of the leading SH instability modes derived from an atmospheric primitive equation model. Decreased, or increased growth rates are generally an indication of whether the particular mode is less or more likely to develop into mature systems and impact on the rainfall.We have conducted this study with July climatological basic states averaged over the periods 1949–1968, 1975–1994, and 1997–2006. As in J. S. Frederiksen and C. S. Frederiksen [1, 2] we find that the largest decreases in growth rates occur for the cyclogenesis, or storm track, modes with reductions of 33% and 37% during 1975–1994 and 1997–2006, respectively, compared to 1949–1968. These modes were shown to have largest impact on rainfall over southwestern Australia and southern Australia in general. The reduction in growth rate of these storm track modes implies a reduced tendency for their development and would explain the reduced rainfall and negative rainfall trends over much of southern Australia.Our analysis has also shown a significant decrease in the growth rate of modes associated with the onset of mid-latitude blocking. For these modes we have found reductions in growth rate of 11% and 19% in 1975–1994 and 1997–2006, respectively, compared with the earliest period. The structures of these modes, at appropriate phases, show a blocking high over New Zealand and a low pressure cell off the eastern seaboard of Australia. This low pressure cell is conducive to the development of rain bearing east coast lows on the east coast of Australia. This region has experienced quite large reductions in rainfall over many decades consistent with the reduced growth rate of these modes and hence the reduced tendency for the development of east coast lows.Not all of Australia has experienced reductions in winter rainfall. Over the central northwest of Western Australia, central eastern Australia, and northern Queensland, there have been increases and positive trends in rainfall. Our study suggests that these positive changes are attributable to the changing nature of two types of weather systems commonly associated with these regions. The first of these is associated with eastward-propagating northwest cloud-bands originating in the Indian Ocean off northwest Australia. Over the last fifty or so years, these NWCB modes have increased their growth rates by around 25% and 45%, respectively, for the 1975–1994 and 1997–2006 basic states compared with 1949–1968. The increasing likelihood of the development of these systems would have a positive impact over the three regions mentioned above and can partly explain the observational changes.Another category of mode that has an impact on rainfall over the central northwest of Western Australia is associated with variability at the intraseasonal timescale. For these modes we see an increase in growth rate of around 29% and 31% for 1975–1994 and 1997–2006 compared to 1949–1968. However, we have also noticed a continual decrease in the period of these modes from 59 days to 35 days and finally 25 days in the last period. Again, the increased growth rate of these modes is consistent with the increased rainfall trends over central northwestern Western Australia.Other modes that show some significant change in their growth rate are those associated with African easterly waves, which are westward-propagating tropical waves generated by the African easterly jet and originating off the east coast of North Africa. These have increased growth rates of 5% and 7% in the latter two periods compared with 1949–1968. However, these have no relevance for Australian rainfall. The other modes, such as high-latitude blocking modes and Antarctic modes, do not show any significant changes in their growth rates or horizontal structure.In subsequent papers, we plan to analyze the changing properties of SH atmospheric modes in other seasons to try to better understand the observed changes in rainfall and climate, in general.



---

*Source: 353829-2011-10-13.xml*