1. Introduction
Recently it has been recognized that there is a possible link between the North Atlantic Oscillation (NAO) or Arctic Oscillation (AO), a dominant dipole mode in the Northern Hemisphere (NH), and the NH surface temperature, which has important implications in understanding the cause of global warming (Hurrell 1995, 1996; Thompson and Wallace 1998; Thompson et al. 2000). Thompson et al. (2000) noted that the observed positive trend in the NAO/AO index over the last three decades of the last century significantly contributed to the observed warming trend over Eurasia and North America. Although the NAO is a natural mode of atmospheric variability (Feldstein 2003; Vallis et al. 2004), its close relationship with warming trends makes it difficult to distinguish whether the positive trend of the NAO at the end of the twentieth century is due to natural variability or due to anthropogenically forced warming.
The NAO trend has been found to be related to ocean variability (Rodwell et al. 1999) or greenhouse gas concentrations (Shindell et al. 1999). However, recent observations show that the decadal trend of the NAO has reversed over the past several years and its link with global warming is not obvious (Cohen and Barlow 2005; Overland and Wang 2005). This means that the factor that is crucial for driving NAO trends is not clearly identified. Wittman et al. (2005) found in a simple stochastic model that a dipolar meridional structure of the NAO/AO is likely due to variations in the north–south position of the jet. However, why the phase of NAO depends on the meridional displacement of the jet is not sufficiently understood. The purpose of this paper is to provide a simple explanation for why the northward (southward) shift of a preexisting jet can excite the positive (negative) phase of the NAO by presenting an analytical solution of topographically forced stationary waves in a specified jet. It is also natural that anticyclonic (cyclonic) wave breaking is closely related to the northward (southward) displacement of a jet prior to the NAO.
This paper is organized as follows: In section 2, the meridional structure of an observed jet prior to the NAO is described. In section 3 we present an analytical solution of topographically forced stationary waves in a specified jet, and find that the northward (southward) shift of the zonal jet tends to produce a dipole component with a low-over-high (high-over-low) anomaly over the North Atlantic through the interaction with the topography of the land–sea contrast (LSC). This dipole component can be regarded as the initial state of subsequent positive (negative) phase NAO events. In addition, a comparison with the numerical solutions is also discussed in this section. A weakly nonlinear NAO model established by Luo et al. (2007c) is extended in section 4 to include the role of a specified jet. In section 5 we discuss how the south–north displacement of a preexisting jet influences the life cycle of positive (negative) phase NAO events. Section 6 presents the relationship between the preexisting jet and subsequent jet. Conclusions and a discussion are summarized in section 7.
2. Meridional profile of Atlantic jet prior to the NAO
The dataset used here is the daily mean, multilevel wind fields from the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR). To clearly see the meridional distribution of an Atlantic jet prior to the NAO, 20 negative and 10 positive phase events are chosen based upon the daily NAO index defined by Benedict et al. (2004). For each phase, the Atlantic jet prior to the NAO is defined as a time mean of zonal winds over the North Atlantic for 30 days before each NAO event begins, which is obtained through filtering out waves with periods from 2 to 21 days.
Figure 1 shows that there is a northward (southward) shift of a jet over the North Atlantic prior to the positive (negative) phase of NAO in the troposphere from the upper layer to the surface. However, what causes the north–south excursion of the preexisting jet is unclear and is beyond the scope of this paper. As we will demonstrate in the next section, the north–south displacement of the preexisting Atlantic jet is crucial for the phase of the subsequent NAO.
3. Asymptotic solution of linear topographically forced Rossby modes and a link to the meridional shift of a preexisting jet
a. Asymptotic solution


In the real atmosphere the zonal jet probably cannot be exactly represented by the mathematical expression given in (2), but can be decomposed into many meridional harmonic modes by using the Fourier expansion. In this case, this zonal jet can be reduced to (2) by neglecting modes higher than the first two modes.
For example, for u0 = 0.85 and α1 = 0.24,








It is found that in a jet with a horizontal shear the approximate solution for a planetary wave anomaly induced by the LSC topography exhibits a rather complicated spatial structure whose flow pattern depends strongly on the amplitude of each component, as given in (5a)–(5f). An interesting point is that ψA12 can have a dipole structure as the preexisting jet exhibits a meridional displacement. This shows that the meridional excursion of the jet plays an important role in the excitation of the dipole mode through an interaction with the LSC topography in the NH.
b. Spatial structures of topographic stationary modes and meridional shift of a preexisting jet
It is further shown that when the mean westerly wind u0 is weaker and is near the resonant westerly wind uc = β/k2 + m2, ψA12 can become dominant compared to ψA11 and ψA13, but the solution of equations is incorrect if u0 = uc is satisfied. This situation will be avoided if u0 is somewhat far from the value of uc. For example, for u0 = 0.85, both ψA11 and ψA13 may be negligible because they are rather small, but ψA12 is significant. In a uniform westerly wind (αi = 0, where i = 1, 2), the topographic planetary wave anomaly can only exhibit a monopole structure because there is only ψA = ψA0. A dominant dipole anomaly ψA12 tends to be excited through the interaction between the jet and LSC topography for a weak westerly wind u0 as the jet departs meridionally from its mean position (α2 ≠ 0). Nevertheless, both ψA11 and ψA13 can become more important as u0 increases further, because β/u0 − (k2 + m2/4) and β/u0 − (k2 + 9m2/4) become smaller. In any case, a dipole anomaly ψA12 can be excited by the LSC topography when the zonal jet undergoes a north–south excursion.
Because the LSC topography is dominated by zonal wavenumber 2 and because an NAO anomaly has a wavenumber 2 structure (Wallace 2000), we can choose k = 2/[6.371 cos(55°N)] in the whole paper (Luo et al. 2007a, b,c; Luo et al. 2008). In addition, xT = 0, h0 = 0.4, m = −2π/Ly, and Ly = 5 are fixed in this section, but u0 and αi (i = 1, 2) are allowed to vary to represent different cases of the zonal jet.
Here, the LSC topography and associated topographically forced wave (ψA = ψA0) in a uniform westerly wind (αi = 0, where i = 1, 2) are shown in Fig. 3 for u0 = 0.85 and u0 = 1.15. It is evident that when the background westerly wind is uniform, the LSC topography tends to induce a planetary wave with a meridional monopole structure as noted above, whose amplitude is dependent upon the setting of this uniform westerly wind. For u0 = 0.85 and α1 = 0.24, ψA and ψA12 are shown in Fig. 4 for α2 = 0.24 and α2 = −0.24. It is noted that when the jet exhibits a northward (southward) shift, the positive anomaly center of ψA can be in the south (north) side of the ocean (North Atlantic), reflecting that ψA12 can have a low-over-high (high-over-low) dipole. Of interest is that ψA12, which has a dipole meridional structure, is dominant, in that both ψA11 and ψA13 are negligible for u0 = 0.85. The spatial structure of this dipole is found to depend upon the direction of the jet departing from its mean position. When the jet has a poleward (equatorward) displacement this dipole anomaly ψA12 exhibits a low-over-high (high-over-low) anomaly in the North Atlantic. Once this dipole anomaly interacts with synoptic-scale eddies upstream, it will be amplified into a positive or negative phase NAO (NAO+ or NAO− hereafter) anomaly (Luo et al. 2007a). This can, to some extent, explain why the poleward (equatorward) displacement of the zonal mean jet is a precursor of NAO+ (NAO−) events, as noted in section 2. Of course, it is also important that the eddy fluxes from preexisting synoptic-scale eddies should have a matching spatial structure to provide a potential for the NAO occurrence (Luo et al. 2007a).
Figure 5 shows the anomaly fields ψA, ψA11, and ψA12 of the topographically induced stationary wave for u0 = 1.15. For this case, ψA13 is not shown here because its amplitude is rather weak. It is found that the amplitude of ψA for u0 = 1.15 is evidently larger than that for u0 = 0.85 as shown in Fig. 4, and at the same time ψA12 becomes weaker even though ψA11 can have a large amplitude. Thus, the setting of u0 does not affect the sign of the spatial structure of ψA12 unless it is particularly large. Certainly, the amplitude of ψA12 is also dominated by the extent of the north–south excursion of the jet.
The results found in this section indicate that a dipole anomaly can be excited by the interaction between the jet and the LSC topography when the jet departs from its mean position. This dipole component due to the meridional wobbling of the jet can provide an initial state of the NAO anomaly although Wittman et al. (2005) noted in a stochastic model that the occurrence of a meridional dipole is a simple consequence of the north–south excursion of the jet center. More recently, Luo et al. (2008) found that anticyclonic wave breaking (AWB) or cyclonic wave breaking (CWB) tend to accompany the occurrence of an NAO+ or NAO− event when an initial state with a low-over-high or high-over-low dipole is prespecified. Thus, it is not difficult to conclude that the AWB (CWB) is more likely to be naturally linked with the preexisting northward (southward) shift of a zonal jet.
c. A comparison with numerical solutions in a beta plane and on sphere
For the jet and topography given in (2) and (3), the numerical solution to (1) can be obtained by solving a Poisson equation using a finite-difference method if ∂ψA/(∂x) = ψAx is defined in (1). It is found that for a northward (southward) jet prior to the NAO+ (NAO−), the center of the topographically induced stationary wave anomaly tends to be in the south (north) side of the ocean (not shown), which is consistent with both the analytical solution ψA presented in (5a) and the observational characteristics of stationary wave anomalies for two phases of the NAO obtained by Luo et al. (2007c, their Fig. 2). This result also holds for spherical geometry, the treatment for which is given in the appendix.
4. Generalized weakly nonlinear NAO model
It needs to be pointed out that when the preexisting planetary wave is stationary, u0 = β/(k2 + m2) will be satisfied so that the solution (5e) is invalid. This difficulty cannot be avoided unless some assumptions are made. To analytically treat the role of the jet shift in the NAO life cycle, we assume αi = ε2αεi to avoid the above difficulty in treating the nonlinear evolution of NAO events. This is because the terms induced by the interaction between the jet shift and the topography can only appear in the higher-order equation if αi = ε2αεi is used. Such an assumption can allow us to obtain the weakly nonlinear asymptotic solutions of the NAO life cycle.













It is found that when the jet has no displacement from its mean position, that is to say, when α2 = 0, (7k) becomes identical to the forced nonlinear Schrödinger (NLS) equation that governs the evolution of an NAO anomaly forced jointly by synoptic-scale eddies and LSC topography derived by Luo et al. (2007c). Thus, (7k) can be considered as an extension of the previously obtained forced NLS equation. If the initial states of the planetary and synoptic waves prior to the NAO are specified, one can predict the life cycle of an NAO event by solving (7k) when the preexisting jet undergoes a north–south shift. It seems that the solutions (7a)–(7k) are not directly related to the solution (5e). In fact, the dipole component described by solution (5e) is assumed to provide an initial state of the NAO represented by (7), but the amplitude evolution of the subsequent NAO is still affected by the meridional shift of the zonal jet during its life cycle.
5. Life cycle of an NAO event and its link to the north–south shift of a jet
In this section, αi is relaxed to be large although the strict assumption of αi = ε2αεi may be violated. In addition, the parameters α1 = 0.18, a0 = 0.17, μ = 1.2, F = 1, and x0 = 2.87/2 are fixed, and without the loss of generality B(x, 0) = 0.4 is chosen as an initial amplitude of the NAO anomaly. The other parameters of synoptic-scale eddies are chosen the same as in Luo et al. (2007c). Such a treatment can allow us to conveniently discuss the impact of the strong north–south excursion of the jet on the phase of the NAO.
a. NAO+ event and meridional shift of a preexisting jet
The aim of this section is to look at the influence of the jet shift on the NAO occurrence if the initial state of the NAO anomaly is preexisting.
For three cases of α2 = 0, α2 = −0.18, and α2 = 0.18, the anomaly field ψNAOA of an eddy-driven NAO+ event (α = 1, m = 2π/Ly) is shown in Fig. 6 for xT = −2 and h0 = −0.4. It is found that for a symmetric jet an isolated NAO+ anomaly with a low-over-high dipole can be formed from a prespecified low-over-high anomaly under the forcing of synoptic-scale eddies (Fig. 6a). Once the zonal jet exhibits a meridional displacement, the northward shift of the jet will increase the amplitude of the NAO+ anomaly and strengthen its downstream energy dispersion during the decaying process (Fig. 6b), but the southward shift of the jet will play a reverse role (Fig. 6c). Thus, the northward displacement of a preexisting jet is a favorable environment for the occurrence of the NAO+ anomaly and strong European blocking events. In a synoptic weather map the CWB is characterized by the poleward shift of warm air and the southward intrusion of cold air (Benedict et al. 2004), but a reverse movement of air characterizes the AWB. As seen in the total field in Fig. 7, the AWB seems to accompany the occurrence of an NAO+ event, and the CWB occurs concurrently over Europe during the decaying phase (Fig. 7a). When the preexisting jet is located more northward, both the AWB in the Atlantic sector and subsequent CWB over the European continent are enhanced, and a strong European blocking is also observed (Fig. 7b). However, if the zonal jet is moved to the south, the NAO+ anomaly is strongly suppressed and no blocking event is seen over the European continent (Fig. 7c). Thus, the northward displacement of a preexisting jet is conducive to the occurrence of NAO+ events, AWB, and European blocking events.
b. Impact of the meridional shift of a preexisting jet on the NAO− event
Observations show that the center of action of the NAO− event is generally located more eastward than the region where the NAO+ event is (Benedict et al. 2004). As a result, we choose xT = −1 and h0 = 0.4 to discuss the evolution of an NAO− event (α = −1, m = −2π/Ly). Figure 8 shows the anomaly fields of ψNAOA for α2 = 0, α2 = 0.18, and α2 = −0.18. It is found that for a symmetric jet an NAO− event can be produced from an initial state with a high-over-low dipole through the forcing of synoptic-scale eddies (Fig. 8a). This NAO− anomaly is enhanced as the zonal jet moves to the south (Fig. 8b), but weakened as the zonal jet moves to the north (Fig. 8c). At the same time, the decay of an NAO− anomaly is also affected by the meridional shift of the preexisting jet. A prominent CWB is also seen in the total field (Fig. 9). As indicated by Luo et al. (2008), the type of synoptic wave breaking is controlled by the spatial structure of the initial planetary wave that determines the phase of the NAO. Thus, the CWB is essentially linked with the negative phase of NAO. When the preexisting jet shows a southward displacement, the CWB is enhanced, thus naturally concluding that the equatorward shift of an Atlantic jet is a favorable environment for the CWB (Fig. 9b). In contrast, the CWB is strongly suppressed once the jet moves to the north (Fig. 9c).
6. Positive feedback between a preexisting jet and its subsequent evolution
As noted in section 4, the poleward (equatorward) displacement of a zonal jet prior to the NAO implicates the more likely occurrence of subsequent NAO+ (NAO−) events. More recently, Luo et al. (2007c) found that an intensified jet can be established and can move to the north (south) during the life cycle of an NAO+ (NAO−) event. Thus, it is not difficult to conclude that there is a positive feedback between a preexisting jet and the subsequent intensified jet during the NAO life cycle. In other words, the northward (southward) shift of the zonal jet prior to the NAO implicates the likely emergence of a more northward (southward) intensified jet associated with an NAO+ (NAO−) event.
7. Conclusions and discussion
In this paper, we have presented an analytical solution of topographically forced stationary modes in a specified jet to clarify why the poleward (equatorward) shift of an Atlantic jet before an NAO event begins is more likely to excite positive (negative) phase NAO events. It is shown that a dipole component can arise from the interaction between a zonal jet and the LSC topography as the zonal jet has a north–south excursion. At the same time, the center of the topographically induced stationary wave anomaly shifts toward the north (south) as the preexisting jet is moved to the south (north). The spatial structure of the topographically induced dipole component depends strongly upon whether a preexisting jet moves to the north or the south. When the jet undergoes a northward (southward) shift, the induced dipole anomaly can have a low-over-high (high-over-low) structure over the North Atlantic, which can be regarded as an initial state or embryo of positive (negative) phase NAO events. As indicated in Benedict et al. (2004) and Luo et al. (2008), anticyclonic (cyclonic) wave breaking always accompanies the occurrence of positive (negative) phase NAO events. Thus, it is possible that during the life cycle of a positive (negative) phase NAO event the anticyclonic (cyclonic) wave breaking is inevitably enhanced when the preexisting jet exhibits a northward (southward) shift.
On the other hand, in the present work a weakly nonlinear NAO model proposed by Luo et al. (2007c) has been extended to examine the impact of the north–south shift of a specified jet on the NAO life cycle if the initial state of the NAO is preexisting. It is shown that the poleward (equatorward) shift of the jet prior to the NAO is indeed favorable for the occurrence of both the positive (negative) phase NAO event and anticyclonic (cyclonic) wave breaking, thus confirming previous observational results (Benedict et al. 2004; Franzke et al. 2004; Rivière and Orlanski 2007).
Simmons and Hoskins (1980) and Thorncroft et al. (1993) have shown that the type of wave breaking is very sensitive to the meridional shear of the basic state zonal winds. Lee and Feldstein (1996) found in an aquaplanet general circulation model that the meridional shear of the zonal wind in the upper troposphere, rather than the barotropic shear as in Thorncroft et al. (1993), is a crucial factor for determining the type of wave breaking. Orlanski (2005) recently noted that the type of wave breaking is in part controlled by SST anomalies created by ENSO because of their direct influence on the low-level baroclinicity. More recently, Woollings et al. (2008) conjectured that the wave breaking arises from the interaction between planetary and synoptic waves. Even so, we find here that the north–south excursion of the zonal jet prior to the NAO seems a crucial factor determining the type of wave breaking. However, what causes the north–south excursion of the zonal jet prior to the NAO is a very interesting problem and deserves further study.
Although the results here are gained in a beta-plane channel, they are also concluded to hold in a spherical geometry from numerical studies in the beta-plane channel model and in a spherical model performed by Charney and DeVore (1979) and Legras and Ghil (1985). Of course, it is impossible to get the analytical solutions presented in this paper if a spherical model is used. However, using the barotropic vorticity equation in a beta-plane channel can make it easier to obtain the analytical solutions.
The authors acknowledge support from the National Outstanding Youth Natural Science Foundation of China under Grant 40325016, the National Natural Science Foundation of China (4057016), the CityU Strategic Research Grant 7002231, and Taishan scholar funding. The authors thank two anonymous reviewers for their useful suggestions in improving this paper.
REFERENCES
Benedict, J. J., , S. Lee, , and S. B. Feldstein, 2004: Synoptic view of the North Atlantic Oscillation. J. Atmos. Sci., 61 , 121–144.
Charney, J. G., , and J. G. DeVore, 1979: Multiple flow equilibria in the atmosphere and blocking. J. Atmos. Sci., 36 , 1205–1216.
Cohen, J., , and M. Barlow, 2005: The NAO, the AO, and global warming: How closely related? J. Climate, 18 , 4498–4513.
Feldstein, S. B., 2003: The dynamics of NAO teleconnection pattern growth and decay. Quart. J. Roy. Meteor. Soc., 127 , 901–924.
Franzke, C., , S. Lee, , and S. B. Feldstein, 2004: Is the North Atlantic Oscillation a breaking wave? J. Atmos. Sci., 61 , 145–160.
Hart, J. E., 1979: Barotropic quasi-geostrophic flow over anisotropic mountains. J. Atmos. Sci., 36 , 1736–1746.
Hurrell, J. W., 1995: Decadal trends in the North Atlantic Oscillation: Regional temperatures and precipitation. Science, 269 , 676–679.
Hurrell, J. W., 1996: Influence of variations in extratropical wintertime teleconnections on Northern Hemisphere temperature. Geophys. Res. Lett., 23 , 665–668.
Lee, S., , and S. B. Feldstein, 1996: Two types of wave breaking in an aquaplanet GCM. J. Atmos. Sci., 53 , 842–857.
Legras, B., , and M. Ghil, 1985: Persistent anomalies, blocking and variations in atmospheric predictability. J. Atmos. Sci., 42 , 433–471.
Luo, D., , A. Lupo, , and H. Wan, 2007a: Dynamics of eddy-driven low-frequency dipole modes. Part I: A simple model of North Atlantic Oscillations. J. Atmos. Sci., 64 , 3–38.
Luo, D., , T. Gong, , and A. Lupo, 2007b: Dynamics of eddy-driven low-frequency dipole modes. Part II: Free mode characteristics of NAO and diagnostic study. J. Atmos. Sci., 64 , 29–51.
Luo, D., , T. Gong, , and Y. Diao, 2007c: Dynamics of eddy-driven low-frequency dipole modes. Part III: Meridional displacement of westerly jet anomalies during two phases of NAO. J. Atmos. Sci., 64 , 3232–3248.
Luo, D., , T. Gong, , and Y. Diao, 2008: Dynamics of eddy-driven low-frequency dipole modes. Part IV: Planetary and synoptic wave-breaking processes during the NAO life cycle. J. Atmos. Sci., 65 , 737–765.
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.
Overland, J. E., , and M. Wang, 2005: The Arctic climate paradox: The recent decrease of the Arctic Oscillation. Geophys. Res. Lett., 32 .L06701, doi:10.1029/2004GL021752.
Rivière, G., , and I. Orlanski, 2007: Characteristics of the Atlantic storm-track eddy activity and its relation with the North Atlantic Oscillation. J. Atmos. Sci., 63 , 241–266.
Rodwell, M. J., , D. P. Rowell, , and C. K. Folland, 1999: Oceanic forcing of the wintertime North Atlantic Oscillation and European climate. Nature, 398 , 320–323.
Shindell, D. T., , R. L. Miller, , G. Schmidt, , and L. Pandolfo, 1999: Simulation of recent northern winter climate trends by greenhouse-gas forcing. Nature, 399 , 452–455.
Simmons, A. J., , and B. J. Hoskins, 1980: Barotropic influences on the growth and decay of nonlinear baroclinic waves. J. Atmos. Sci., 37 , 1679–1684.
Thompson, D. W. J., , and J. M. Wallace, 1998: The Arctic Oscillation signature in the wintertime geopotential height and temperature fields. Geophys. Res. Lett., 25 , 1297–1300.
Thompson, D. W. J., , J. M. Wallace, , and G. C. Hegerl, 2000: Annular modes in the extratropical circulation. Part II: Trends. J. Climate, 13 , 1018–1036.
Thorncroft, C. D., , B. J. Hoskins, , and M. E. McIntyre, 1993: Two paradigms of baroclinic-wave life-cycle behaviour. Quart. J. Roy. Meteor. Soc., 119 , 17–55.
Vallis, G. K., , E. P. Gerber, , P. J. Kushner, , and B. A. Cash, 2004: A mechanism and simple dynamical model of the North Atlantic Oscillation and annular modes. J. Atmos. Sci., 61 , 264–280.
Wallace, J. M., 2000: North Atlantic Oscillation/annular mode: Two paradigms–one phenomenon. Quart. J. Roy. Meteor. Soc., 126 , 791–805.
Wittman, M. A. H., , A. J. Charlton, , and L. M. Polvani, 2005: On the meridional structure of annular modes. J. Climate, 18 , 2119–2122.
Woollings, T. J., , B. J. Hoskins, , M. Blackburn, , and P. Berrisford, 2008: A new Rossby wave-breaking interpretation of the North Atlantic Oscillation. J. Atmos. Sci., 65 , 609–626.
APPENDIX
Linear Topographically Forced Stationary Wave Equation on a Sphere and Its Numerical Solution


Here, we fix the parameters θS = 30°N, θN = 80°N, α1 = 0.24, and h0 = 0.4. The horizontal distribution of the topography with zonal wavenumber 2 on sphere can be seen in Fig. A1.
For the topography given in Fig. A1 and the prescribed parameters, a finite-difference method is used to obtain ψAλ through solving a Poisson equation that is satisfied by ψAλ = ∂ψA/(∂λ). On this basis, ψA(λ, θ) is obtained by integrating ψAλ with respect to λ using a Fourier expansion.
For u0 = 1.15 and α1 = 0.24, the topographically induced stationary wave anomaly is shown in Fig. A2 for α2 = 0, α2 = −0.24, and α2 = 0.24. It is found that for a symmetric jet (α2 = 0), the center of the topographically induced stationary wave anomaly is not in the middle part of the latitude band from 30° to 80°N because of the effect of spherical geometry. Instead, it is in the south of this latitude band. This feature is different from the result obtained in a beta-plane channel. However, as we will note below, on a sphere the characteristics of a topographically induced stationary wave anomaly depending upon the meridional shift of a zonal jet are basically consistent with the analytical and numerical solutions obtained in a beta-plane channel. This point is clear in Fig. A2. When the zonal jet is shifted to the north, the positive center of topographically induced stationary wave anomaly is in the more south of the latitude band in the ocean (trough) as well as a negative center seen in the north of the latitude band. However, when this jet is displaced toward the south, the positive center over the ocean is shifted to the north. The meridional shift of the center of the forced stationary wave anomaly reflects that a dipole component can be excited by the interaction between the jet shift and topography, and can be regarded as the initial states of NAO events. This conclusion also holds for other values of u0 (not shown). In conclusion, the stationary waves on the sphere are found to be basically consistent with those in a beta plane, and they do not propagate in the meridional direction. This reflects that the stationary waves in the real atmosphere can, to large extent, be captured by the analytical and numerical solutions in a beta plane.

Meridional profile of the zonal jet at (a) 300 and (b) 850 mb and (c) the vertical mean over the North Atlantic between 300 and 850 mb for the 1-month mean prior to the positive (negative) phase of NAO. The solid and dashed curves denote the positive and negative phases, respectively.
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

Meridional profile of a specified jet, in which the dotted–dashed, solid, and dotted curves correspond in turn to α2 = 0, α2 = 0.24, and α2 = −0.24.
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(a) Distribution of the land–sea contrast topography, in which the solid and dashed curves represent the topographic ridge and trough, respectively, and topographically induced stationary wave anomaly for different uniform westerly winds (b) u0 = 0.85 and (c) u0 = 1.15, in which the solid and dashed curves represent positive and negative anomalies, respectively. The contour interval (CI) is 0.2.
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

Anomaly fields ψA and ψA12 of topographically induced stationary wave components in a specified jet (CI = 0.1) in which the dashed and solid lines denote the negative and positive anomalies, respectively: (a) northward-shifting jet and (b) southward-shifting jet.
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

As in Fig. 4 but for anomaly fields ψA, ψA11, and ψA12.
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

Anomaly field ψNAO of an eddy-driven positive phase NAO event (CI = 0.2): (a) symmetric jet (α2 = 0), (b) northward-shifting jet (α2 = −0.18), and (c) southward-shifting jet (α2 = 0.18).
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

As in Fig. 6 but for the total field (CI = 0.3).
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

Anomaly field ψNAOA of an eddy-driven negative phase NAO event (CI = 0.2): (a) symmetric jet (α2 = 0), (b) southward-shifting jet (α2 = 0.18), and (c) northward-shifting jet (α2 = −0.18).
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

As in Fig. 8 but for the total field (CI = 0.3).
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

(Continued)
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

Fig. A1. Horizontal distribution of wavenumber 2 topography on a sphere (CI = 0.1), in which the solid and dashed lines represent the ridge and trough, respectively.
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1

Fig. A2. Anomalies of topographically induced stationary waves on sphere for u0 = 1.15 and different jets (CI = 0.5), in which the solid and dashed lines correspond to the positive and negative anomalies, respectively: (a) symmetric jet (α2 = 0), (b) northward-shifting jet (α2 = −0.24), and (c) southward-shifting jet (α2 = 0.24).
Citation: Journal of the Atmospheric Sciences 65, 6; 10.1175/2007JAS2560.1