1. Introduction
Most of the midlatitudes are under the influence of storm tracks. The dynamics of storm tracks has been the central theme of climate dynamics, inspiring numerous studies (e.g., Hartmann 1974; Blackmon 1976; Blackmon et al. 1977; Hoskins et al. 1989; Hoskins and Valdes 1990; Chang 2001; Chang et al. 2002; Bengtsson et al. 2006). The storm tracks have a conspicuous seasonal cycle in their locations: they shift equatorward in step with the jet stream from fall to midwinter and then migrate poleward afterward (Hoskins et al. 1989; Nakamura 1992). In addition, at interannual time scales, ENSO can impact the latitudinal position of the storm track and the jet stream, with the warm phase leading to an equatorward contraction (Robinson 2002; Seager et al. 2003; Orlanski 2005; L’Heureux and Thompson 2006; Lu et al. 2008; Chen et al. 2008; Brayshaw et al. 2008) and vice versa. The issue of the latitudinal position of the storm-track and jet-stream shift is particularly topical, given that the poleward shift of the midlatitude storm tracks is deemed to be one of the most robust predicted features under global warming (e.g., Kushner et al. 2001; Yin 2005; Lorenz and DeWeaver 2007; Lu et al. 2008; Chen et al. 2008). Nevertheless, the mechanistic research on the storm track has hitherto been largely focused on the intensity (e.g., James 1987; Chang 2001; Nakamura 1992; Straus and Shukla 1988), and much is yet to be learned as to the mechanisms that control the location of the storm track.
Hoskins and Valdes (1990) perceived the storm track as being self-maintaining in the sense that storm-track eddies in general are most vigorous downstream of the region of peak baroclinicity, and the mixing of temperature by eddies is relatively benign where the baroclinicity is largest. Further, the enhanced baroclinicity over the storm-track entrance region is actively maintained by condensational heating, which in turn is caused by the cyclones themselves. The self-maintaining nature of the storm track implies that the transient component and the mean flow (and the associated thermal structure) are just integral parts of the same phenomenon: the storm track will always ensure that the baroclinicity is locally maximum nearby. The classic theory of linear baroclinic instability states that eddies grow by tapping the background available potential energy (e.g., Eady 1949; Charney 1947; Phillips 1954); thus, one may argue that the storm track is excited by the mean flow baroclinicity. Some have attempted to explain the shift of eddy field by a measure of mean baroclinicity such as the Eady growth rate (e.g., Inatsu et al. 2003; Brayshaw et al. 2008). The linear baroclinic instability theory no doubt has its value in understanding the existence of the storm track; however, in many cases one has to explain where the mean baroclinicity or the variation thereof comes from in the first place. As will be shown later for the cases considered in the current study, the baroclinicity evaluated based on the mean temperature in radiative–convective equilibrium gives little clue as to which way the storm track should move, whereas the latitudes of the maxima of eddy statistics vary in tandem with the mean wind when eddy adjustment is taken into account in the maintenance of the mean temperature.
In this study, we will use the eddy-driven surface westerly (EDSW) as the proxy for the storm track. The EDSW, as the term implies, results from the balance between the vertically integrated eddy momentum convergence and the surface drag; the surface drag acting on it is locally balanced by the Coriolis force associated with the poleward flow at the bottom of the Ferrel cell. Thus, intuitively, the maximum EDSW can be used as a proxy for the position of the storm track, and any predictive power obtained for the westerly wind can potentially be carried over to the other properties of the midlatitude storm track.
It is beyond the scope of the present study to construct a rigorous closure theory for the storm track. Our goal, instead, is to strive for the most primary dynamic factors that can be used to infer the movement of the storm-track proxy, the EDSW. The most important attribute of these factors should be that they are quantifiable from the conditions given for the problem, such as the sea surface temperatures (SSTs) so that through them one might be able to link the variation in the position of the storm track/EDSW to the given conditions. This exercise, which will be referred to as scaling for convenience, will be carried out for a suite of aquaplanet model simulations with two distinct GCMs under the same prescribed SST boundary conditions. If the proposed scaling can be validated, it should help shed light on the key processes that are responsible for the shift of the storm track/EDSW under different climate conditions. Two GCMs with distinct physics and parameterizations are used here to assess the robustness of our scaling.
The whole endeavor is also motivated by the observation that the best-measured quantity in both modern and paleoclimate records is the surface temperature, and one may potentially tap into this resource by forming a heuristic relationship that links surface temperature to the midlatitude mean flow and eddy statistics. Any success along this line may lead to some predictive power for the changes of storm track in the past and future.
The models and experiments will be detailed in section 2. Model simulations of the storm tracks and the EDSW will be reported in section 3. Section 4 lays out the rationale of the scaling, which is to be validated in section 5. We conclude with summary and discussion on the possible implications for climate change in the past and future.
2. Model description and simulations
a. Idealized GCM
The idealized GCM is an ice-free, land-free model (or so-called aquaplanet model) consisting of various simplified physical parameterizations coupled to a spectral dynamical core that solves the primitive equations (Frierson et al. 2006, 2007a). For longwave radiative transfer, a gray scheme is used with fluxes that are a function only of temperature and not of water vapor or other constituents. Water vapor itself is a prognostic variable, but there is no liquid water or clouds. Condensation and moist convection is represented by a simplified version of the Betts–Miller convection scheme (Betts 1986; Betts and Miller 1986; Frierson et al. 2007a) and a grid-scale condensation scheme. The simplified model physics also includes a simplified Monin–Obukhov surface flux scheme and a K-profile boundary layer scheme. All simulations in this study are run at a resolution of T42 in the horizontal and 25 levels in the vertical.
b. Full GCM
The full GCM simulations in this paper are the same simulations originally used to study poleward heat transports in the study of Caballero and Langen (2005). The model is a comprehensive GCM, with full representations of clouds, radiation, convection, and other physics. The atmospheric model used for these simulations is the Parallel Community Climate Model version 3 (PCCM3), which is the atmospheric component of the Fast Ocean-Atmosphere Model (FOAM; Jacob 1997). The model uses the physical parameterizations of the National Center for Atmospheric Research (NCAR) Community Climate System Model, version 3.6 (CCM3.6; Kiehl et al. 1996) and the dynamical core of the NCAR CCM2. The full GCM is run at T42 resolution, with 18 vertical levels. Sea ice is specified where SST is below 0°C.
c. Boundary conditions


3. Results of model simulations
We first examine the EDSW position in the suite of simulations with the two models. As shown in Fig. 1c, the axis of storm track [measured as vertically integrated eddy kinetic energy (EKE)] is in excellent alignment with the latitude of the maximum EDSW in the idealized model simulations. To a lesser extent, this relationship also holds for the full GCM (Fig. 1d). In general, it is justifiable that the shift of the EDSW is representative of that of the EKE of the storm track, with some cautionary discretion applied to the full model. This result, in a context of our aquaplanet simulations, corroborates the notion that mean wind moves in concert with the transients.
Both the idealized and full GCMs simulate a monotonic poleward shift of the EDSW/storm track with increasing Tm, probably the most robust feature in all the simulations.1 Another robust feature, not shown, is the intensification of storm activity and the EDSW with increasing ΔT (see Fig. 3a in Caballero and Langen 2005). However, the two GCMs show a distinct difference in the sensitivity of the storm-track/EDSW position to ΔT. For the idealized model, the position of the EDSW is also a monotonic function of ΔT, moving poleward with increasing ΔT. This is not the case for the full GCM. In the middle of the Tm–ΔT space, as ΔT increases from its lowest value, the EDSW first shifts equatorward and then poleward after passing its most equatorial position at the intermediate value similar to the conditions of present climate. The distinct behavior between the two models rules out a possible unifying scaling theory for the position of storm track/EDSW valid for the whole range of SST parameters and for both models. Therefore, in this study, we will only strive for an anomaly-wise scaling for the anomalous shift of EDSW resulting from perturbations with respect to a certain chosen reference state.
Further inspection on the zonal wind profiles of the full GCM suggests that this peculiar phenomenon of an equatorward shift with ΔT may be related to the large separation between the subtropical jet and eddy-driven jet when the global mean SST is warm but the SST gradient is weak (the upper-left domain in Fig. 1b). These cases with a zonal wind profile characterized by a distinct split-jet condition are marked with crosses in Fig. 1b. It is worth noting that the wind profiles in the split-jet regime bear a strong resemblance to those of the bifurcation discussed by Lee and Kim (2003) and Chen and Plumb (2009). In the case of Lee and Kim (2003), the bifurcation occurs when the subtropical jet is intensified by equatorial heating; as a result, the midlatitude disturbances are more subject to the waveguiding effect of a stronger subtropical jet (which itself is stable in position) and hence shift equatorward. The sharp regime transition toward a merged jet at ∼30° latitude never takes place in our case. Rather, the eddy-driven jet merely progresses gradually equatorward as the subtropical jet grows in strength with increasing ΔT. As ΔT approaches its middle range, the eddy-driven jet and the subtropical jet start to merge, and a different mechanism begins to come into play in the response of EDSW position to ΔT, a subject of focus of the present scaling study to be elaborated in the following sections. A typical transition from split- to single-jet regime in the full GCM with increasing ΔT is exemplified by Fig. 2, which shows the profiles of 400-hPa zonal wind for the cases of Tm = 15°C, with ΔT varying from 10° to 50°C at increments of 5°C.
4. Rationale for scaling


















Observations (Edmon et al. 1980) demonstrate that the meridional structure of the vertically integrated PV flux (















5. Validation
a. Relationship (9) based on model simulations
Figure 4 shows how well the relationship (9) holds for the idealized model between the fractional differential changes of ξ estimated from the model-simulated mean temperatures and the variation in the position of the EDSW (the total value is used). A qualitatively very similar result is found from the full model simulations and thus not shown. Note that the result of Fig. 4 is not very sensitive to which exact pixel point around (Tm = 20°C, ΔT = 40°C) in the Tm–ΔT space is chosen for the reference state. Although linearity does exist in general, it holds best for the cases neighboring the reference state, as highlighted by the circles in Fig. 4, but starts to deteriorate when (Tm, ΔT) values deviate markedly from the reference value. However, it is still encouraging to see that the scaling relationship (9) is valid for perturbations in both Tm and ΔT as long as they are small. Although this result should not be taken as evidence for the predictive power of (9), because ξ is estimated from the model simulations, this corroborates the anticipated connection between the differential change of baroclinicity (and hence the upward eddy EP fluxes) and the shift of the eddy-driven wind. For example, if the baroclinicity increases on the poleward side of the storm track relative to the equatorward side, the source of the baroclinic wave activity shifts poleward. Consequently, the divergent wave propagation in the upper troposphere, which corresponds to convergence of eddy momentum flux, should shift poleward accordingly, driving a poleward shift of the EDSW through momentum balance.
We next examine how ξ and its components [i.e., the static stability and the tropospheric temperature (meridional) gradient] vary under each SST condition. Figure 5 shows the profiles of ξ in the idealized simulations broken down into different ΔT groups. Overall, the magnitude of ξ increases with ΔT (with an intriguing exception for Tm = 0°C wherein the midlatitude static stability increases at a greater rate than the tropospheric temperature gradient with increasing ΔT) and decreases substantially with Tm because of the increase of static stability. A similar tendency is also found in the simulations with the full GCM. As a result, the midlatitude isentropic slope flattens considerably with increasing Tm, in contrast to the constant isentropic slopes as rationalized from the arguments of baroclinic adjustment (Stone 1978; Stone and Nemet 1996) or neutral supercriticality because of weak nonlinearity in the atmospheric eddy–eddy interactions (Schneider and Walker 2006). This result suggests that the existing dry theories for the midlatitude adjustment are inadequate to describe the midlatitude tropospheric thermal stratification under large SST boundary perturbations, at least for these idealized experiments examined here (see also Juckes 2000; Frierson 2008). In the meantime, the flattening of the isentropic slope with increasing mean temperature is consistent with previous studies on the effects of moisture on baroclinic eddies (Stone and Yao 1990; Held 1978b; Gutowski 1992). For example, Held (1978b) studied the effect of adding hydrological cycle to a dry two-layer climate model and calculated the ratio of the vertical eddy heat flux, weighted by the static stability, to the meridional eddy heat flux weighted by the meridional gradient of potential temperature, a ratio that approximates the mixing slope divided by the isentropic slope. Held found that this ratio was 0.55 in the dry model, close to the value given by classical theories of baroclinic instability, but increased to 0.90 in his moist case. The increase of the ratio of the mixing slope with respect to the mean isentropic slope was also found in Gutowski (1992) to be the key effect of condensation in the life cycle of midlatitude eddies. Given the fact that vertical transport of dry static energy by eddies
The most prominent aspect of the static stability is its increase with Tm as shown in Fig. 6 for the idealized model (see also Frierson 2008). A qualitatively similar result is also found for the full model. The increase is most prevalent in the tropics, where the vertical thermal profile approximately follows the moist adiabat. The increase also spreads poleward outside of the territory controlled by Hadley cell. From the coldest to the warmest Tm case, the average static stability near the edge of the Hadley cell or the equatorward flank of the reference storm track increases by a factor of 8 (see Fig. 9). For each ΔT group, the static stability increases preferentially over the subtropical latitudes relative to higher latitudes as Tm increases, pushing the peak of ξ poleward, a phenomenon that is typical of the response to the greenhouse gas (GHG)–induced global warming (Frierson 2006; Lu et al. 2008; Yin 2005).
The tropospheric temperature gradients in the simulations are depicted in Fig. 7, together with the SST gradient (black curves). For the warm Tm cases, the gradient profiles exhibit two peaks, a subtropical one associated with the thermally forced subtropical jet (a same suite of simulations running on an axisymmetric, eddy-free configuration elucidate unambiguously the subtropical peak, not shown) and a subpolar one associated with the eddy-driven jet. For the moderate and weak Tm cases, these two peaks are merged and indistinguishable. If one defines the edge of the Hadley cell as the subtropical peak of the temperature gradient, the width of the Hadley cell can hardly exceed 35° latitude, in accordance with the notion that the meridional scale of the Hadley cell is set by the location where the thermally driven subtropical jet first becomes baroclinically unstable (Held et al. 2000; Lu et al. 2007; Frierson et al. 2007b). However, the location of the EDSW associated with the eddy momentum convergence varies over a much wider range between 30° and 65° latitude (see also Figs. 1, 3). It is particularly notable that, for the cases with large Tm and ΔT, the EDSW maximum and the associated maximum temperature gradient are located significantly poleward of the maximum SST gradient (located at 45° for all cases). This poleward enhancement of the gradient reflects the notion of a self-maintaining eddy-driven jet: a self-maintaining jet preserves and sometimes reinforces the midlatitude gradient and places the gradient on the poleward side of the imposed baroclinicity through the transformed Eulerian mean overturning circulation (Robinson 2000, 2006; Chen and Plumb 2009). The creation of baroclinicity in the life cycle of eddies has also been explicitly demonstrated in the seminal work of Simons and Hoskins (1978) and Hoskins (1983).
The self-maintaining nature of the EDSW poses a challenge for us to predict the storm-track or jet position based on the atmospheric temperature gradient, because it is largely the result of the baroclinic eddy adjustment. It may be more desirable to predict the location of the storm track using the baroclinicity estimated from the temperature of the axisymmetric, eddy-free simulations, and in this approach the effects in both static stability and tropospheric temperature gradient of the eddy adjustment are excluded from the predicting factors. However, the evaluation of ξ using the axisymmetric simulations turns out to be fruitless. Figure 8 contrasts the locations of the maximum baroclinicity ξ estimated from the original eddy-permitted simulations and the eddy-free simulations by the same idealized model and their alignment/misalignment with locations of the storm tracks. Although ξ from the original simulations serves as a good indicator of the axis of the storm track, ξ from the eddy-free simulations just scatters between 50° and 60° latitude and shows little agreement with the storm track locations. It even fails to shift poleward with increasing Tm, the most robust behavior of the storm track in all our eddy-permitted model simulations. As a result, the baroclinicity from the eddy-free state provides almost no clue as to which way the storm track should shift. This could be partly due to the specific prescriptions of the SST conditions, of which the maximum gradients are always fixed at 45° latitude. On the other hand, the variation of static stability, to the extent that it can be predicted solely from the information of SST regardless of the position of the storm track, as will be demonstrated next, may be of some guiding value regarding the shift of the storm track and the tropospheric temperature gradient.
b. Implication from a static stability theory
Here, we estimate (δξ+ − δξ−)/ξpm making only use of the variation of static stability, which is predicted based on the theory of Juckes (2000) from the given SST boundary conditions. Note that, unlike the estimation in section 5a, no information of the simulated atmospheric temperature is used except that of the reference state.





Here, for our specific purpose of scaling, we predict the (dry) bulk static stability, which is defined as the potential temperature difference between the tropopause and the surface, and averaged over the 30°–50° latitudinal band (on the equatorward flank of the reference EDSW position): that is,
In view of the fact that the static stability does not vary as much over the poleward side as compared with the equatorward side of the jet (Fig. 6), we neglect the change in δξ+ and approximate δξ− with a theoretical estimate δξt− (with subscript t indicating theoretical estimation) in the evaluation of (δξ+ − δξ−)/ξpm. As such, (δξ+ − δξ−)/ξpm is simplified to be −δξt−/ξpm. In estimating −δξt−/ξpm, the temperature gradient component of δξt− is fixed at the reference value and theoretical estimates at the equatorial flank of the jet are used for the static stability part in δξt−; ξpm is approximated by an average between the reference value at the poleward flank of the jet (which is therefore fixed) and the theoretically estimated value at the equatorward flank of the jet (which varies with SSTs). The result is presented in Fig. 10 for the idealized model, showing a tight relationship between the quantity −δξt−/ξpm and the axis of the EDSW/storm track. The extended good linear fit to the large perturbation cases is not expected and should be interpreted with caution. In fact, when the baroclinicity changes (because of temperature gradients and/or static stability) on the poleward of the reference jet are taken into consideration, the linear fit deteriorates toward the situation as depicted in Fig. 4. One obvious reason, but not the only reason, is that, for the wide range of variation of the westerly jet/storm-track positions between 30° and 60° latitude, many may violate the constraint of small perturbations with respect to the 47° latitude reference jet position. Nevertheless, the overall results of this scaling effort corroborate the notion that eddy activity and the eddy-induced wind do feel the changes in subtropical static stability, with an increase in the latter causing a poleward shift of the former. Quantitatively, a unit (or 100%) increase in −δξt−/ξpm can lead to ∼7° shift of the EDSW/storm track.
A dynamical scaling that works would hopefully be only weakly dependent on the model configuration. Thus, to assess model dependency, we apply this same scaling to the simulations with the full GCM over the similar suite of SST boundary conditions. The result for all the cases turns out to be a significant deviation from a linear relationship (Fig. 11a), confirming our suspicion that the linearity between −δξt−/ξpm and the EDSW shift for large perturbation cases in the idealized model may be fortuitous. This may also be due to the regime behavior of the position of the EDSW over the wide range of ΔT, as discussed in section 3. Indeed, when applying the scaling to the Tm–ΔT domain wherein the wind profiles are characteristic of single jet and the EDSW shows a monotonic relationship to both Tm and ΔT, the linear relationship between −δξt−/ξpm and the actual shift resumes (Fig. 11b). It is especially encouraging to note that the slope of the linear relationship in Fig. 11b is very similar to that in the idealized model, a ∼6° shift for per unit increase of −δξt−/ξpm. In summary, the simple static stability-based scaling has shown to be valid and robust in the shift of the EDSW/storm track simulated by two rather different aquaplanet GCMs, insofar as the wind profile resides in a same dynamic regime and the SST perturbations are not so large as to cause wind regime transitions.
6. Summary and concluding remarks
The midlatitude storm track shifts poleward in simulations for the third Coupled Model Intercomparison Project (CMIP3) under the GHG-induced global warming. No widely accepted theory exists to rationalize this phenomenon, except that some studies suggest that it links to the rise of the tropopause (Williams 2006; Lorenz and DeWeaver 2007) and the attendant eddy feedbacks through the eddy phase speed changes (Chen and Held 2007; Chen et al. 2008; Lu et al. 2008). Here, by diagnosing a suite of aquaplanet simulations with two distinct GCMs under specified SST boundary conditions, we demonstrate that the variation of the zonal index (or the shift of the surface westerly winds) can be thought of as being constrained by the structure of the lower-tropospheric Eliassen–Palm activity flux (predominantly the heat flux) coming out of the lower troposphere. The latter, parameterized as a quantity measuring the tropospheric baroclinicity, shows a strong sensitivity to the perturbation in the subtropical static stability. In both an idealized GCM and a GCM with full physics, the preferential stabilization at the equatorial flank of the storm track/EDSW can, at least within the climate regime similar to the present climate in the Northern Hemisphere winter, shift poleward the baroclinicity and hence the production of the eddy activity flux. Because the source of the horizontal eddy activity flux (equivalent to the convergence of eddy momentum flux) in the upper troposphere determines the position of the EDSW through the momentum balance upon vertical integration, a poleward displacement of the source of the eddy activity flux heralds a poleward shift of EDSW. Through this chain of reasoning, the displacement in the position of EDSW may be predictable, provided the information of static stability can be derived exclusively from the given SST conditions. Finally, a good relationship results between the static stability-based estimation of baroclinicity and the position of EDSW in both the idealized and the full GCMs, suggesting the importance of static stability as a key dynamical factor for understanding the shift of storm track and the associated eddy-mean flow interaction.
This assertion should not be confused with the scaling theory for the width of the Hadley cell advocated by Frierson et al. (2007b), wherein the subtropical static stability has been found to scale with the Hadley cell width in the same aquaplanet simulations examined here (see their Figs. 1 and 2). We stress the distinction between the scaling for the Hadley cell width and the scaling constructed here for the storm-track position: the static stability in the former acts on the subtropical baroclinicity associated with the thermally forced subtropical jet within the Hadley cell, whereas the static stability in the latter acts through the midlatitude baroclinicity to influence the storm track. Moreover, the Hadley cell width scaling is much less susceptible to the jet regime transition, as one can discern from inspecting our Figs. 1a,b and Fig. 1 in Frierson et al. (2007b).
The importance of static stability in the shift of storm track and EDSW is further corroborated by the analysis of the simulations by the state-of-the-art CMIP3 models: as one can infer from Fig. 6b in Lu et al. (2008), the models with greater increase in static stability to the immediate equatorward side of the storm track tends to show a larger poleward displacement of EDSW. In the meantime, the scaling put forward in this study should not be considered contradictory but complementary to the mechanism proposed by Lorenz and DeWeaver (2007) for the case of global warming, which emphasizes the role of the tropopause rise. The rise of tropopause is associated with an upper-tropospheric warming and a stratospheric cooling under the direct and indirect effect of the greenhouse gas forcing; they can modulate the upper-tropospheric wave propagation and eddy momentum flux and consequently the position of the eddy-driven westerly wind (Chen et al. 2007; Chen and Held 2007; Chen et al. 2008). The relative contributions to the EDSW shift from the mechanism of static stability versus that of tropopause rise remain to be quantified, a topic of our ongoing investigation.
The anomaly-wise scaling proposed here is shown to be valid for perturbations in both the global mean and the equator-to-pole gradient of the SST up to an order of 10 K, it thus may have some implications on the change of the storm track during the past of the earth’s climate. For example, during the glacial periods when the global mean temperature was about 10 K cooler than the present, our scaling would suggest the storm track should be several degrees equatorward relative to today. This speculation is consistent with the higher-than-today loading of mineral dust from ice core records (e.g., Petit et al. 1999), attributable to the larger exposure of the subtropical desert to the uplift by midlatitude storm systems as they retreat equatorward in a colder climate (Chylek et al. 2001). For another application, during the mid-Holocene (∼6000 yr ago), the solar radiation in the Southern Hemisphere had a weaker equator-to-pole gradient during the austral spring followed by an overall dimming during summer, which could conceivably lead to an equatorward movement of the austral summer storm track. Whether this is true remains to be verified by paleoclimatic data.
Acknowledgments
The perceptive and constructive comments from Paul O’Gorman and other two anonymous reviewers have dramatically improved the clarity of the manuscript. We also thank David Straus for his internal review at COLA. JL is supported by the startup fund at George Mason University and the COLA omnibus fund from NSF Grant 830068, NOAA Grant NA09OAR4310058, and NASA Grant NNX09AN50G. GC is supported by the startup fund at Cornell University. DMF is supported by NSF Grants ATM-0846641 and AGS-0936069.
REFERENCES
Bengtsson, L., K. I. Hodges, and E. Roeckner, 2006: Storm tracks and climate change. J. Climate, 19 , 3518–3543.
Betts, A. K., 1986: A new convective adjustment scheme. Part I: Observational and theoretical basis. Quart. J. Roy. Meteor. Soc., 112 , 677–692.
Betts, A. K., and M. J. Miller, 1986: A new convective adjustment scheme. Part II: Single column tests using GATE wave, BOMEX, and arctic air-mass data sets. Quart. J. Roy. Meteor. Soc., 112 , 693–709.
Blackmon, M. L., 1976: A climatological spectral study of the 500 mb geopotential height of the Northern Hemisphere. J. Atmos. Sci., 33 , 1607–1623.
Blackmon, M. L., J. M. Wallace, N-C. Lau, and S. L. Mullen, 1977: An observational study of the Northern Hemisphere wintertime circulation. J. Atmos. Sci., 34 , 1040–1053.
Branscome, L. E., 1983: A parameterization of transient eddy heat flux on a beta-plane. J. Atmos. Sci., 40 , 2508–2521.
Brayshaw, D. J., B. Hoskins, and M. Blackburn, 2008: The storm-track response to idealized SST perturbations in an aquaplanet GCM. J. Atmos. Sci., 65 , 2842–2860.
Bretherton, F. P., 1966: The propagation of groups of internal gravity waves in a shear flow. Quart. J. Roy. Meteor. Soc., 92 , 466–480.
Caballero, R., and P. L. Langen, 2005: The dynamic range of poleward energy transport in an atmospheric general circulation model. Geophys. Res. Lett., 32 , L02705. doi:10.1029/2004GL021581.
Chang, E. K. M., 2001: GCM and observational diagnoses of the seasonal and interannual variations of the Pacific storm track during the cool season. J. Atmos. Sci., 58 , 1784–1800.
Chang, E. K. M., S. Lee, and K. L. Swanson, 2002: Storm track dynamics. J. Climate, 15 , 2163–2183.
Charney, J. G., 1947: The dynamics of long wave in a baroclinic westerly current. J. Meteor., 4 , 135–162.
Chen, G., and I. M. Held, 2007: Phase speed spectra and the recent poleward shift of Southern Hemisphere surface westerlies. Geophys. Res. Lett., 34 , L21805. doi:10.1029/2007GL031200.
Chen, G., and R. A. Plumb, 2009: Quantifying the eddy feedback and the persistence of the zonal index in an idealized atmospheric model. J. Atmos. Sci., 66 , 3707–3719.
Chen, G., I. M. Held, and W. A. Robinson, 2007: Sensitivity of the latitude of the surface westerlies to surface friction. J. Atmos. Sci., 64 , 2899–2915.
Chen, G., J. Lu, and D. M. W. Frierson, 2008: Phase speed spectra and the latitude of surface westerlies: Interannual variability and global warming trend. J. Climate, 21 , 5942–5959.
Chen, G., R. A. Plumb, and J. Lu, 2010: Sensitivities of zonal mean atmospheric circulation to SST warming in an aqua-planet model. Geophys. Res. Lett., 37 , L12701. doi:10.1029/2010GL043473.
Chylek, P., G. Lesins, and U. Lohmann, 2001: Enhancement of dust source area during past glacial periods due to change of the Hadley circulation. J. Geophys. Res., 106 , (D16). 18477–18485.
Eady, E. T., 1949: Long waves and cyclone waves. Tellus, 1 , 33–52.
Edmon, H. J. Jr., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen–Palm cross sections for the troposphere. J. Atmos. Sci., 37 , 2600–2616.
Frierson, D. M. W., 2006: Robust increases in midlatitude static stability in simulations of global warming. Geophys. Res. Lett., 33 , L24816. doi:10.1029/2006GL027504.
Frierson, D. M. W., 2008: Midlatitude static stability and simple and comprehensive general circulation models. J. Atmos. Sci., 65 , 1049–1062.
Frierson, D. M. W., I. M. Held, and P. Zurita-Gotor, 2006: A gray-radiation aquaplanet moist GCM. Part I: Static stability and eddy scale. J. Atmos. Sci., 63 , 2548–2566.
Frierson, D. M. W., I. M. Held, and P. Zurita-Gotor, 2007a: A gray-radiation aquaplanet moist GCM. Part II: Energy transports in altered climates. J. Atmos. Sci., 64 , 1680–1693.
Frierson, D. M. W., J. Lu, and G. Chen, 2007b: The width of Hadley circulation in simple and comprehensive general circulation models. Geophys. Res. Lett., 34 , L18804. doi:10.1029/2007GL031115.
Green, J. S., 1970: Transfer properties of the large scale eddies and the general circulation of the atmosphere. Quart. J. Roy. Meteor. Soc., 96 , 157–185.
Gutowski, W. J. Jr., 1992: Life cycles of moist baroclinic eddies. J. Atmos. Sci., 49 , 306–319.
Hartmann, D. L., 1974: Time spectral analysis of mid-latitude disturbances. Mon. Wea. Rev., 102 , 348–362.
Held, I. M., 1978a: The vertical scale of an unstable baroclinic wave and its importance for eddy heat flux parameterizations. J. Atmos. Sci., 35 , 572–576.
Held, I. M., 1978b: The tropospheric lapse rate and climate sensitivity: Experiments with a two-level atmospheric model. J. Atmos. Sci., 35 , 2083–2098.
Held, I. M. Coauthors 2000: The general circulation of the atmosphere. 2000 Program of Study in Geophysical Fluid Dynamics, Woods Hole, Massachusetts, Woods Hole Oceanographic Institute, 1–54.
Hoskins, B. J., 1983: Modelling of the transient eddies and their feedback on the mean flow. Large-Scale Dynamical Processes in the Atmosphere, B. J. Hoskins and R. P. Pearce, Eds., Academic Press, 169–199.
Hoskins, B. J., and P. J. Valdes, 1990: On the existence of storm tracks. J. Atmos. Sci., 47 , 1854–1864.
Hoskins, B. J., H-H. Hsu, I. N. James, M. Matsutani, P. D. Sardeshmukh, and G. H. White, 1989: Diagnostics of the global atmospheric circulation based on ECMWF analysis 1979-1989. World Climate Research Programme Rep. 27 and World Meteorological Organisation Tech. Doc. 326, 217 pp.
Inatsu, M., H. Mukougawa, and S-P. Xie, 2003: Atmospheric response to zonal variations in midlatitude SST: Transient and stationary eddies and their feedback. J. Climate, 16 , 3314–3329.
Jacob, R., 1997: Low frequency variability in a simulated atmosphere ocean system. Ph.D thesis, University of Wisconsin—Madison, 159 pp.
James, I. N., 1987: Suppression of baroclinic instability in horizontally sheared flows. J. Atmos. Sci., 44 , 3710–3720.
Juckes, M. N., 2000: The static stability of the midlatitude troposphere: The relevance of moisture. J. Atmos. Sci., 57 , 3050–3057.
Kiehl, J. T., J. J. Hack, G. B. Bonan, B. A. Boville, B. P. Brieglieb, D. L. Williamson, and P. J. Rasch, 1996: Description of the NCAR Community Climate Model (CCM3). National Center for Atmospheric Research Tech. Note TN-420, 159 pp.
Kushner, P. J., I. M. Held, and T. L. Delworth, 2001: Southern Hemisphere atmospheric circulation response to global warming. J. Climate, 14 , 2238–2249.
Lee, S., and H-K. Kim, 2003: The dynamical relationship between subtropical and eddy-driven jets. J. Atmos. Sci., 60 , 1490–1503.
L’Heureux, M. L., and D. W. J. Thompson, 2006: Observed relationships between the El Niño–Southern Oscillation and the extratropical zonal-mean circulation. J. Climate, 19 , 276–287.
Lorenz, D. J., and E. T. DeWeaver, 2007: Tropopause height and the zonal wind response to global warming in the IPCC scenario integrations. J. Geophys. Res., 112 , D10119. doi:10.1029/2006JD008087.
Lu, J., G. A. Vecchi, and T. Reichler, 2007: Expansion of the Hadley cell under global warming. Geophys. Res. Lett., 34 , L06805. doi:10.1029/2006GL028443.
Lu, J., G. Chen, and D. M. W. Frierson, 2008: Response of the zonal mean atmospheric circulation to El Niño versus global warming. J. Climate, 21 , 5835–5851.
Nakamura, H., 1992: Midwinter suppression of baroclinic wave activity in the Pacific. J. Atmos. Sci., 49 , 1629–1642.
Orlanski, I., 2005: A new look at the Pacific storm track variability: Sensitivity to tropical SSTs and to upstream seeding. J. Atmos. Sci., 62 , 1367–1390.
Pauluis, O., A. Czaja, and R. Korty, 2008: The global atmospheric circulation on moist isentropics. Science, 321 , 1075–1078.
Pauluis, O., A. Czaja, and R. Korty, 2010: The global atmospheric circulation in moist isentropic coordinates. J. Climate, 23 , 3077–3093.
Pavan, V., and I. M. Held, 1996: The diffusive approximation for eddy fluxes in baroclinically unstable jets. J. Atmos. Sci., 53 , 1262–1272.
Petit, J. R. Coauthors 1999: Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica. Nature, 399 , 429–436.
Phillips, N. A., 1954: Energy transformations and meridional circulations associated with simple baroclinic waves in a two-level, quasi-geostrophic model. Tellus, 6 , 273–286.
Ring, M. J., and R. A. Plumb, 2008: The response to a simplified GCM to axisymmetric forcings: Application of the fluctuation–dissipation theorem. J. Atmos. Sci., 65 , 3880–3898.
Robinson, W. A., 2000: A baroclinic mechanism for the eddy feedback on the zonal index. J. Atmos. Sci., 57 , 415–422.
Robinson, W. A., 2002: On the midlatitude thermal response to tropical warming. Geophys. Res. Lett., 29 , 1190. doi:10.1029/2001GL014158.
Robinson, W. A., 2006: On the self-maintenance of midlatitude jets. J. Atmos. Sci., 63 , 2109–2122.
Schneider, T., and C. C. Walker, 2006: Self-organization of atmospheric macroturbulence into critical states of weak nonlinear eddy–eddy interactions. J. Atmos. Sci., 63 , 1569–1586.
Schneider, T., P. A. O’Gorman, and X. Levine, 2010: Water vapor and the dynamics of climate changes. Rev. Geophys., 48 , RG3001. doi:10.1029/2009RG000302.
Seager, R., N. Harnik, Y. Kushnir, W. Robinson, and J. Miller, 2003: Mechanisms of hemispherically symmetric climate variability. J. Climate, 16 , 2960–2978.
Simons, A. J., and B. J. Hoskins, 1978: The life cycles of some nonlinear baroclinic waves. J. Atmos. Sci., 35 , 414–432.
Stone, P. H., 1972: A simplified radiative-dynamical model for the static stability of rotating atmospheres. J. Atmos. Sci., 29 , 405–418.
Stone, P. H., 1978: Baroclinic adjustment. J. Atmos. Sci., 35 , 561–571.
Stone, P. H., and M-S. Yao, 1990: Development of a two-dimensional zonally averaged statistical-dynamical model. Part III: The parameterization of the eddy fluxes of heat and moisture. J. Climate, 3 , 726–740.
Stone, P. H., and B. Nemet, 1996: Baroclinic adjustment: A comparison between theory, observations, and models. J. Atmos. Sci., 53 , 1663–1674.
Straus, D. M., and J. Shukla, 1988: A comparison of a GCM simulation of the seasonal cycle of the atmosphere with observations. Part II: Stationary waves and transient fluctuations. Atmos.—Ocean, 26 , 575–607.
Williams, G. P., 2006: Circulation sensitivity to tropopause height. J. Atmos. Sci., 63 , 1954–1961.
Yin, J. H., 2005: A consistent poleward shift of the storm tracks in simulations of 21st century climate. Geophys. Res. Lett., 32 , L18701. doi:10.1029/2005GL023684.
APPENDIX
Formulation of Shift
Empirically speaking, dipolar wind anomalies centered about the axis of a jet can lead to a shift of the jet. Here, we derive a functional relationship between the shift of the jet and the dipolar wind anomalies considering perturbations in both latitudinal position and the magnitude of the jet but with the shape kept the same.










To relate the shift of wind to the differential change of the baroclinicity ξ, we substitute (7′) into (A4) and make use of another assumption that fractional change of the peak wind speed is proportional to that of the baroclinicity: that is, Urm/Upm ∼ ξrm/ξpm, which is a reasonable first-order approximation to the actual model simulations. The result is simply the proportionality relation (9).

(top) Latitude of EDSW simulated in (a) the idealized and (b) the full GCMs. In (b), the crosses mark the cases of double-jet condition. (bottom) The latitude of the EDSW vs that of the axis of the storm track in (c) the idealized and (d) the full GCMs. The storm track is measured as the vertically integrated EKE. The dots in (c) and (d) are color coded based on the corresponding Tm
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

(top) Latitude of EDSW simulated in (a) the idealized and (b) the full GCMs. In (b), the crosses mark the cases of double-jet condition. (bottom) The latitude of the EDSW vs that of the axis of the storm track in (c) the idealized and (d) the full GCMs. The storm track is measured as the vertically integrated EKE. The dots in (c) and (d) are color coded based on the corresponding Tm
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
(top) Latitude of EDSW simulated in (a) the idealized and (b) the full GCMs. In (b), the crosses mark the cases of double-jet condition. (bottom) The latitude of the EDSW vs that of the axis of the storm track in (c) the idealized and (d) the full GCMs. The storm track is measured as the vertically integrated EKE. The dots in (c) and (d) are color coded based on the corresponding Tm
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

The profiles of 400-hPa zonal wind in the simulations of the full GCM for the cases of Tm = 15°C, ΔT varying from 10 to 50 K at increments of 5 K. The cases identified to be in split-jet regimes in Fig. 1b are highlighted by the thick lines.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

The profiles of 400-hPa zonal wind in the simulations of the full GCM for the cases of Tm = 15°C, ΔT varying from 10 to 50 K at increments of 5 K. The cases identified to be in split-jet regimes in Fig. 1b are highlighted by the thick lines.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
The profiles of 400-hPa zonal wind in the simulations of the full GCM for the cases of Tm = 15°C, ΔT varying from 10 to 50 K at increments of 5 K. The cases identified to be in split-jet regimes in Fig. 1b are highlighted by the thick lines.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Profiles of
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Profiles of
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
Profiles of
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

(δξ+ − δξ−)/ξpm estimated from model-simulated temperature vs the actual axis of the storm track in the idealized GCM. The circles highlight the cases neighboring the reference state.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

(δξ+ − δξ−)/ξpm estimated from model-simulated temperature vs the actual axis of the storm track in the idealized GCM. The circles highlight the cases neighboring the reference state.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
(δξ+ − δξ−)/ξpm estimated from model-simulated temperature vs the actual axis of the storm track in the idealized GCM. The circles highlight the cases neighboring the reference state.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Plot of ξ (m−1 s−1) in the set of simulations by the idealized model, broken down into six ΔT groups. Curves in each group/panel are color coded corresponding to the value of Tm as in Fig. 4.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Plot of ξ (m−1 s−1) in the set of simulations by the idealized model, broken down into six ΔT groups. Curves in each group/panel are color coded corresponding to the value of Tm as in Fig. 4.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
Plot of ξ (m−1 s−1) in the set of simulations by the idealized model, broken down into six ΔT groups. Curves in each group/panel are color coded corresponding to the value of Tm as in Fig. 4.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

As in Fig. 5, but for
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

As in Fig. 5, but for
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
As in Fig. 5, but for
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

As in Fig. 5, but for the tropospheric temperature gradients (K m−1).
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

As in Fig. 5, but for the tropospheric temperature gradients (K m−1).
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
As in Fig. 5, but for the tropospheric temperature gradients (K m−1).
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

(a) Latitude of the max ξ vs the location of the max EKE in the original eddy-permitted idealized GCM; (b) As in (a), but ξ values are estimated from the temperature simulated by the axisymmetric version of the idealized GCM. The cases in each panel are color coded based on the corresponding Tm
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

(a) Latitude of the max ξ vs the location of the max EKE in the original eddy-permitted idealized GCM; (b) As in (a), but ξ values are estimated from the temperature simulated by the axisymmetric version of the idealized GCM. The cases in each panel are color coded based on the corresponding Tm
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
(a) Latitude of the max ξ vs the location of the max EKE in the original eddy-permitted idealized GCM; (b) As in (a), but ξ values are estimated from the temperature simulated by the axisymmetric version of the idealized GCM. The cases in each panel are color coded based on the corresponding Tm
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Bulk static stability (K) estimated using the Juckes theory vs that simulated directly by the idealized GCM.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Bulk static stability (K) estimated using the Juckes theory vs that simulated directly by the idealized GCM.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
Bulk static stability (K) estimated using the Juckes theory vs that simulated directly by the idealized GCM.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

−δξt−/ξpm vs the displacement of the EDSW in the idealized GCM. The theoretical estimate of −δξt−/ξpm is not evaluated from GCM mean fields, but rather based on the Juckes theory for the bulk static stability and the meridional temperature gradient fixed at the reference value.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

−δξt−/ξpm vs the displacement of the EDSW in the idealized GCM. The theoretical estimate of −δξt−/ξpm is not evaluated from GCM mean fields, but rather based on the Juckes theory for the bulk static stability and the meridional temperature gradient fixed at the reference value.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
−δξt−/ξpm vs the displacement of the EDSW in the idealized GCM. The theoretical estimate of −δξt−/ξpm is not evaluated from GCM mean fields, but rather based on the Juckes theory for the bulk static stability and the meridional temperature gradient fixed at the reference value.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

(a) Theoretical estimate of −δξt−/ξpm vs the displacement of EDSW axes in the full GCM for all cases; (b) as in (a), but only for the cases of ΔT ≥ 30°C identified as the single-jet regime. The blue, green, yellow, and red colors represent Tm = 10°, 15°, 20°, and 25°C, respectively.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

(a) Theoretical estimate of −δξt−/ξpm vs the displacement of EDSW axes in the full GCM for all cases; (b) as in (a), but only for the cases of ΔT ≥ 30°C identified as the single-jet regime. The blue, green, yellow, and red colors represent Tm = 10°, 15°, 20°, and 25°C, respectively.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
(a) Theoretical estimate of −δξt−/ξpm vs the displacement of EDSW axes in the full GCM for all cases; (b) as in (a), but only for the cases of ΔT ≥ 30°C identified as the single-jet regime. The blue, green, yellow, and red colors represent Tm = 10°, 15°, 20°, and 25°C, respectively.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Schematic for the shift of the jet under perturbations in both the position and the amplitude of the jet.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1

Schematic for the shift of the jet under perturbations in both the position and the amplitude of the jet.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
Schematic for the shift of the jet under perturbations in both the position and the amplitude of the jet.
Citation: Journal of the Atmospheric Sciences 67, 12; 10.1175/2010JAS3477.1
In a set of simulations with a similar idealized GCM over an extended range of climate forcing (Schneider et al. 2010), the surface westerlies actually move equatorward with warming over part of the range of the simulations, despite the fact that the heat flux and near-surface EKE move poleward, an interesting counterexample to the cases studied here.