## 1. Introduction

The space–time evolution of interannual sea surface temperature (SST) anomalies associated with the Antarctic circumpolar wave (ACW) displays general eastward propagation around the Southern Ocean (Fig. 1). White and Peterson (1996) found the phase relationship between SST and sea level pressure (SLP) anomalies such that associated poleward (equatorward) meridional surface wind (MSW) anomalies approximately coincide with warm (cool) SST anomalies over most of the Southern Ocean. In the presence of an apparent advection of SST anomalies by the eastward flow of the Antarctic Circumpolar Current (ACC), the persistent phase relationships between SST, SLP, and MSW anomalies suggested to White and Peterson (1996) that maintenance of the ACW involves coupling between ocean and atmosphere. Here, we propose that this occurs because SST anomalies drive MSW anomalies in the overlying atmosphere and corresponding zonal surface wind (ZSW) anomalies drive an anomalous SST tendency through anomalous meridional Ekman heat advection.

Recently, White and Chen (1998; hereafter WC98) began testing this coupling hypothesis by simulating the atmospheric response to interannual SST anomalies in the Southern Ocean with an equilibrium climate model for the lower troposphere. From analysis of European Centre for Medium-Range Weather Forecasts (ECMWF) and National Centers for Environmental Prediction (NCEP) datasets, they found anomalous vortex stretching of the lower troposphere over the Southern Ocean occurring in response to SST-induced midlevel diabatic heating anomalies on interannual timescales (see the appendix). They found this leading to low-level anomalous convergence balanced by anomalous meridional advection of planetary vorticity, yielding poleward (equatorward) MSW anomalies over warm (cool) SST anomalies, as observed. Since the relationship between extratropical SST anomalies and midlevel diabatic heating anomalies holds for winter and not summer, this relationship extends to interannual timescales because of the former’s dominance over the latter in the annual average. In the present study, we demonstrate how this atmospheric response to anomalous SST forcing can produce a positive feedback onto the ocean that yields the characteristic phasing of the observed ACW.

We begin by examining the NCEP Reanalysis for 15 years from 1982 to 1996 (Kalnay et al. 1996) to determine more precise phase relationships between SST, SLP, and ZSW anomalies associated with the ACW. This gives us clues as to the nature of coupling in the ACW. Second, we conduct a scale analysis that isolates specific processes believed responsible for the anomalous upper-ocean heat balance associated with the ACW. This simplifies the coupling, allowing us to construct an analytical model that yields a dispersion relation for free coupled waves of the system that match many of the wave characteristics observed in the ACW. Finally, we construct a numerical coupled model by linking a global equilibrium climate model for the lower troposphere (WC98) to a global heat budget model for the upper ocean. The latter ocean model balances the anomalous SST tendency with the advection of SST anomalies by the mean broadscale flow of the ACC and with the advection of the background SST field by anomalous meridional Ekman flow, the latter driven by the anomalous ZSW response to SST anomalies. The numerical coupled model allows us to simulate much of the space–time evolution exhibited by the ACW in Fig. 1.

## 2. Spiral structure in global SST and SLP anomaly patterns associated with the ACW

Recently, Robertson (1996) demonstrated how coupling between the extratropical ocean and atmosphere on interdecadal timescales requires the spatial scales of covarying SST and SLP anomalies to become anisotropic and their spatial orientation to become tilted at some angle from the zonal direction. He finds this tilt arising from the maintenance of covarying ocean and atmosphere signals against dissipation, as we shall demonstrate here in this study. As such, this tilt in the spatial orientation of covarying oceanic and atmospheric anomaly patterns becomes a signature not only for coupling but for the type of coupling between ocean and atmosphere. A similar tilt can be seen in the spatial orientation of SST anomalies associated with the ACW in Fig. 1, where the spatial orientation of elongated positive and negative SST anomalies can be seen aligned in the southwest–northeast direction as the ACW propagates eastward across the Indian and Pacific sectors of the Southern Ocean. A more realistic view of this is given on a polar stereographic projection (Fig. 2), where this tilt in spatial orientation in Fig. 1 can be seen to be a manifestation of a spiral structure in global spatial patterns of covarying SST and SLP anomalies around the Southern Ocean. Within this spiral structure, SLP anomalies can be seen straddling the boundary between warm and cool SST anomalies, consistent with poleward (equatorward) MSW anomalies aligned with warm (cool) SST anomalies, as explained by WC98. Most importantly, the spiral structure allows westward (eastward) ZSW anomalies to come into approximate alignment with the warm (cool) SST anomalies as well.

If the spiral structure in global spatial patterns of covarying SST and SLP anomalies in the Southern Ocean is a signature for ocean–atmosphere coupling in the ACW, we should be able to use it to gain some insight into how covarying SST and SLP anomalies interact to produce ACW characteristics. Closer inspection in Fig. 2 of the SLP anomalies along the spiraling boundaries between warm and cool SST anomalies finds their extrema displaced poleward from the latitude connecting maximum SST anomalies. In particular, this can be seen with regard to the high SLP anomaly in the Pacific sector and with the two low SLP anomalies south of Africa and Australia. In these cases, associated ZSW anomalies between 45° and 65°S are in position to generate anomalous meridional Ekman advection of the background SST field (in addition to anomalous sensible-plus-latent heat flux and/or anomalous vertical mixing at the base of the near-surface mixed layer), yielding an anomalous SST tendency displaced poleward and eastward from SST anomalies of the same sign. Not only do these ZSW anomalies suggest a source for eastward (and poleward) propagation of the ACW, but that portion of the resultant anomalous SST tendency that is in phase with the SST anomalies suggests a source for maintaining the ACW against dissipation (i.e., anomalous heat loss to the atmosphere).

Therefore, we now ask whether an ocean–atmosphere coupled model can be constructed that will simulate this observed displacement between covarying SST and SLP anomalies. At the same time, we ask whether these displacements can explain the eastward propagation of the ACW observed in Fig. 1, while accounting for maintenance of the ACW against dissipation and its attendant global spiral structure in Fig. 2.

## 3. Scale analysis for the ACW

A necessary approach toward understanding coupling between ocean and atmosphere in the ACW lies in identifying the salient terms in anomalous heat and vorticity budgets in both the upper ocean and lower troposphere. White and Chen have already conducted a scale analysis for the equilibrium response of the lower troposphere to SST-induced midlevel diabatic heating on large zonal spatial scales like those of the ACW (see the appendix). Here we do the same for oceanic heat and vorticity budgets, where scales and model parameters for both are summarized in Table 1.

*h*

*t*

*c*

_{R}

*h*

*x*

*τ*

*f,*

*h*is the depth anomaly of the interface representing the mean depth of the main pycnocline, positive downward;

*c*

_{R}is the nondispersive Rossby wave speed, negative westward (i.e., −

*β*g′

*H*

_{O}/

*f*

^{2});

*β*is the meridional derivative of the Coriolis parameter (

*f*); g′ is the scale of reduced gravity for the upper layer of the ocean (i.e., 0.01 m s

^{−2});

*H*

_{O}is the mean depth scale of the main pycnocline (i.e., 500 m); and

*τ*is anomalous kinematic horizontal wind stress at the sea surface (i.e., in units of m

^{2}s

^{−2}). The much thicker lower layer is assumed to be at rest. We obtain characteristic wavelength scales of the ACW from inspection of Figs. 1 and 2, yielding a zonal wavelength scale of 180° (i.e.,

*L*

_{X}= 1.1 × 10

^{7}m near 55°S) and a meridional wavelength scale of 40° (i.e.,

*L*

_{Y}= 4.4 × 10

^{6}m). As such, the average slope of the spiral structure in covarying global SST and SLP anomaly patterns from 45° to 65° latitude in Fig. 2 is approximately Δ

*y*/Δ

*x*= −

*k*/

*λ*= (

*L*

_{Y}/

*L*

_{X}) = 0.4, where

*k*is the zonal wavenumber (i.e., positive eastward) and

*λ*is the meridional wavenumber (negative southward). We obtain the characteristic period scale

*L*

_{T}= 4.5 yr from White and Peterson (1996), representing observed periodicity ranging from 4 to 5 yr over the 16 years of record from 1979 to 1994. These scales, together with −1.2 × 10

^{−4}s

^{−1}for

*f*and −0.004 m s

^{−1}for

*c*

_{R}near 55°S, finds the relative scale of the stretching term (i.e., ∂

*h*/∂

*t*) on the lhs of Eq. (3.1) dominating that of the meridional advection of planetary vorticity (i.e.,

*C*

_{PX}∂

*h*/∂

*x*) by more than an order of magnitude. Therefore, Eq. (3.1) reduces to the expression for Ekman pumping; that is,

*h*

*t*

*τ*

^{X}

*y*

*f,*

*τ*

^{X}/∂

*y*also dominates ∂

*τ*

^{Y}/∂

*x*in curl(

*τ*) = (∂

*τ*

^{Y}/∂

*x*− ∂

*τ*

^{X}/∂

*y*), with

*τ*

^{X}and

*τ*

^{Y}zonal and meridional components, respectively, of the anomalous surface wind stress. So, whatever influence SST anomalies have upon the atmosphere in the ACW, the feedback upon the main pycnocline in the ocean is through Ekman pumping in Eq. (3.2), principally through ∂

*τ*

^{X}/∂

*y,*the latter proportional to the meridional gradient of the anomalous ZSW.

*T*

*t*

*KT*

*u*

_{O}

*T*

*x*

*T*/∂

*y*

*υ*

_{G}

*υ*

_{E}

*T*/∂

*z*

*w*

_{E}

*T*is the uniform anomalous temperature in this layer, representing anomalous SST;

*u*

_{O}represents the mean zonal geostrophic flow of the broadscale ACC, uniform over this layer;

*KT*is dissipation representing the anomalous sensible-plus-latent heat loss to the atmosphere;

*υ*

_{G}and

*υ*

_{E}represent anomalous meridional geostrophic and Ekman flow, respectively, both uniform over the near-surface mixed layer, the latter proportional to anomalous ZSW; and

*w*

_{E}represents the vertical velocity anomaly at the base of the near-surface mixed layer. By assuming that climatological atmospheric forcing during winter deepens the near-surface mixed layer to the top of the main pycnocline at a scale depth,

*H*

_{p}, of 100 m, then fluctuations from year to year in the depth of the main pycnocline from anomalous Ekman pumping will yield cool (warm) SST anomalies when the pycnocline is displaced upward (downward) from its long-term mean depth. For scaling purposes, we propose that

*w*

_{E}in Eq. (3.3) is less than or equal to ½∂

*h*/∂

*t*from Eq. (3.2) (where

*w*

_{O}= 0 at the sea surface) and that (

^{−2}°C m

^{−1}).

In subsequent development of the coupled model, we will continue to assume winter conditions because WC98 found anomalous midtropospheric level diabatic heating in the Southern Ocean responding to SST anomalies exclusively during winter, not during summer. This derives from the influence that SST anomalies have upon the location and/or intensity of aggregate synoptic storm activity during winter. White and Chen inferred from this that the phase relationship between SST and MSW anomalies observed on interannual timescales by White and Peterson (1996) represents an annual average of an inherent winter (and possibly autumn and spring) mechanism.

Additional assumptions have been made in the heat budget model in Eq. (3.3). In the feedback from atmosphere to ocean, ZSW anomalies influence more than just anomalous meridional Ekman heat advection represented in this model; they also influence anomalous wind mixing at the base of the near-surface mixed layer and the anomalous sensible-plus-latent heat loss at the sea surface, both of which are in phase with the meridional Ekman heat advection in their collective influence on the anomalous SST tendency (Miller et al. 1994). But, for the sake of simplicity, we have ignored the latter two mechanisms in favor of the former, recognizing that they may contribute to effects attributed in this model solely to the anomalous meridional Ekman heat advection. Also, we have ignored the influence that anomalous cloud cover has upon the anomalous net radiation flux in its contribution to the anomalous SST tendency. Furthermore, we have ignored the near-surface Ekman deflection of anomalous geostrophic winds in the near-surface planetary boundary layer of the lower troposphere, assuming that anomalous oceanic Ekman flow (i.e., vertically averaged over the near-surface mixed layer of the ocean) is directed to the left of anomalous surface wind stresses in the Southern Hemisphere, the latter aligned with anomalous surface geostrophic winds.

From the description of horizontal and vertical temperature distributions and zonal geostrophic surface currents (0/4000 db) in the Southern Ocean (Olbers et al. 1992), we estimate the scales and model parameters used in Eq. (3.3) from spatial averages taken over the Southern Ocean from 45° to 65°S, summarized in Table 1. The dissipation parameter *K* corresponds to an *e*-folding timescale for anomalous sensible-plus-latent heat loss to the atmosphere of 3 months. This estimate for *K*^{−1} is obtained from 4(*ρ*_{A}/*ρ*_{O})(*C*_{PA}/*C*_{PO})(*C*_{D}*u*_{A}/*H*_{P}), balancing the dissipation of anomalous heat storage in the near-surface mixed layer of depth *H*_{P} with the anomalous sensible-plus-latent heat loss to the atmosphere, where the factor of 4 assumes that anomalous latent heat loss is three times the anomalous sensible heat loss, the latter directly related to anomalous SST. This factor could be larger or smaller depending upon location, season, and year of interest.

*π*radians. As such, each term is important in the anomalous upper ocean heat budget of the ACW and cannot be ignored. On the rhs of Eq. (3.3), relative scaling between

*υ*

_{G}

*T*/∂

*y*

*w*

_{E}

*T*/∂

*z*

*h*) for the two-layer model ocean; that is,

*w*

_{E}

*T*/∂

*z*

*υ*

_{G}

*T*/∂

*y*

*υ*

_{E}

*T*/∂

*y*

*w*

_{E}

*T*/∂

*z*

*u*

_{A}); that is,

*γ*is the coefficient (i.e., 1.3 × 10

^{−5}m s

^{−1}) linearly relating anomalous ZSW to anomalous zonal surface wind stress, estimated from (

*ρ*

_{A}/

*ρ*

_{O})(

*C*

_{D}/

*u*

_{A}). In Eq. (3.5),

*υ*

_{E}is meridional Ekman flow (e.g., Gill 1982) while the vertical velocity

*w*

_{E}is half the Ekman pumping velocity given in Eq. (3.2). Substituting Eq. (3.5) into Eq. (3.3), and applying the scales and parameters evaluated previously and summarized in Table 1, finds relative scales for

*υ*

_{E}

*T*/∂

*y*

*w*

_{E}

*T*/∂

*z*

*π*radians. This indicates that meridional Ekman heat advection dominates both vertical Ekman heat advection and meridional geostrophic heat advection on the rhs of Eq. (3.3) in the feedback from atmosphere to ocean in the ACW. This is consistent with the meridional Ekman heat advection scenario discussed previously in the examination of the phase differences between ACW SST and SLP anomaly patterns in Fig. 2.

## 4. An analytic coupled model for the ACW

The foregoing scale analysis given for the upper ocean and that given by WC98 for the lower troposphere greatly simplifies modeling of ocean–atmosphere coupling in the ACW. In the lower troposphere, it allows us to relate SST anomalies to vertical-averaged zonal geostrophic wind anomalies (assumed equivalent to ZSW anomalies), and in the upper ocean it allows us to relate ZSW anomalies to the anomalous SST tendency through anomalous meridional Ekman heat advection. The question now becomes whether such a simple coupled model can simulate the wave characteristics of the ACW in a manner consistent with observed phase relationships in Figs. 1 and 2.

*T*

*t*

*KT*

*u*

_{O}

*T*

*x*

*T*/∂

*y*

*γ*

*fH*

_{P}

*u*

_{A}

*T*/∂

*t*) following a fluid parcel traveling zonally with the mean speed

*u*

_{O}of the broadscale ACC and damped by the anomalous turbulent loss of heat to the atmosphere (i.e., parameterized by

*KT*), is balanced by anomalous heat advection stemming from anomalous meridional Ekman flow [see Eq. (3.5a)] operating on the mean meridional SST gradient (i.e.,

*T*/∂

*y*

*γ*and

*K*to be constant, neglecting their dependency upon wind speed.

*KT*) in Eq. (4.1) influences the overlying atmosphere and instigates the feedback to the ocean. How this works has been examined by White and Chen and is summarized here in the appendix. There, we learn that some portion or all of

*KT*instigates midlevel anomalous diabatic heating (

*Q*) in Eq. (A.2), which leads to the anomalous MSW (

*υ*

_{A}) response in Eq. (A.5). Yet coupling between ocean and atmosphere in association with the ACW requires understanding of the anomalous ZSW (

*u*

_{A}) response to SST anomalies. This is determined by substituting anomalous MSW from Eq. (A.5) into the expression for the conservation of mass in the lower troposphere in Eq. (A.4), yielding

*u*

_{A}

*x*

*fα*

*βC*

^{2}

_{A}

*T*

*y*

*β*

*f*

*T*

*f*/

*β*and the low-level convergence is proportional to SST in Eq. (A.4). The resulting expression indicates that ∂

*u*

_{A}/∂

*x*derives both from SST(

*T*) and its meridional derivative. Therefore, the influence of anomalous ZSW upon the anomalous heat budget of the upper ocean in Eq. (4.1) can only be achieved by taking the zonal derivative of the latter, yielding

*T*/∂

*y*

*T*) constitutes the coupled partial differential equation for the model ACW.

*T*=

*T*

_{O}exp(

*k*

_{M}

*x*+

*λ*

_{M}

*y*−

*σ*

_{M}

*t*) into it, yielding

*σ*

_{M}is the complex frequency,

*k*

_{M}is the zonal wavenumber, and

*λ*

_{M}is the meridional wavenumber of the model ACW. The zonal phase speed of the model ACW (i.e.,

*c*

_{ACWX}=

*σ*

_{M}/

*k*

_{M}) depends upon the broadscale flow of the ACC, augmented by a zonal coupling speed, (

*c*

_{CX}); that is,

*T*/∂

*y*

*f*] < 0; so, the zonal coupling speed (

*c*

_{CX}) is eastward, in the same direction as the broadscale mean flow (

*u*

_{O}) of the ACC.

*K*+ [

*T*/∂

*y*

*γα*)/(

*βC*

^{2}

_{A}

*H*

_{P})](

*λ*

_{M}/

*k*

_{M})}. For the coupled model of the ACW to maintain its amplitude against dissipation, this expression is set to zero, thereby defining the magnitude and direction of the meridional wavenumber for stable coupled waves; that is,

*λ*

_{M}

*k*

_{M}

*T*/∂

*y*

*γα*

*βC*

^{2}

_{A}

*H*

_{P}

*K*

^{−1}

*T*/∂

*y*

*k*

_{M}> 0, then

*λ*

_{M}< 0. So, as the coupled model ACW propagates eastward around the Southern Ocean, its amplitude remains stable if it has a poleward wavenumber component

*λ*

_{M}given by Eq. (4.6). Therefore, the meridional phase speed of the ACW is obtained by substituting Eq. (4.6) for

*λ*

_{M}into the expression for

*c*

_{ACWY}=

*σ*

_{M}/

*λ*

_{M}in the real part of Eq. (4.4), yielding

*c*

_{CY}) is highly dispersive depending upon

*λ*

^{−2}

_{M}

*β*/

*f*) and the dissipation timescale K

^{−1}. Note that this coupling phase speed is poleward and decreases with latitude through the inverse Coriolis parameter (i.e., negative in the Southern Hemisphere).

The magnitude of the zonal coupling speed *c*_{CX} can be estimated from Eq. (4.5) for scales and model parameters estimated previously and summarized in Table 1. We can see that *c*_{CX} is highly dispersive, depending upon ^{−2}_{M}*k*_{M} = 0.5(*σ*_{M}/*u*_{O})[1 + (1 + 4*T*/∂*y**γα*)(*u*_{O}/*σ*^{2}_{M}*fC*^{2}_{A}*H*_{P}))^{1/2}]}, yielding *k*_{M} = 0.7(*σ*_{M}/*u*_{O}) = 6 × 10^{−7} m^{−1} (i.e., *L*_{XM} = 2*π*/*k*_{M} = 1.1 × 10^{7} m) at 55°S for a model frequency (*σ*_{M} = 4.5 × 10^{−8} s^{−1}) corresponding to a period of 4.5 yr. This model zonal wavelength *L*_{XM} is the same as that observed earlier (i.e., *L*_{X} = 1.1 × 10^{7} m at 55°S) from consideration of Figs. 1 and 2. From knowledge of this model zonal wavenumber, the zonal coupling speed (*c*_{CX}) can be computed from Eq. (4.5b), yielding *c*_{CX} = 0.04 m s^{−1}. Take note that both the zonal wavenumber and the zonal coupling speed are independent of the dissipation timescale *K*^{−1}. The zonal coupling speed is approximately the same as the scale magnitude for the broadscale flow of the ACC (i.e., *u*_{O} = 0.05 m s^{−1}). Summing the two speeds in Eq. (4.5a) yields an estimate for *c*_{ACWX} of approximately 0.09 m s^{−1}, similar to the observed speed of 0.08 m s^{−1} (i.e., *L*_{X}/*L*_{T} in Table 1).

The meridional coupling speed (*c*_{CY}) in Eq. (4.7) and its associated meridional wavenumber (*λ*_{M}) in Eq. (4.6) depend upon the zonal wavenumber (*k*_{M}). The meridional wavelength (*L*_{YM}*K*^{−1}. For scales and model parameters estimated previously and summarized in Table 1, we find stable coupled waves requiring *λ*_{M} = −2.2*k*_{M} = −1.2 × 10^{−6} m^{−1} (i.e., *L*_{YM} = 2*π*/*λ*_{M} = −5.2 × 10^{6} m), directed poleward. This model meridional wavelength is similar to that observed earlier (i.e., *L*_{Y} = −4.4 × 10^{6} m) from inspection of Figs. 1 and 2. Moreover, it yields a southwest–northeast alignment in crests and troughs of the model ACW, with an average slope of Δ*y*/Δ*x* = −*k*_{M}/*λ*_{M} = 0.5, similar to that observed earlier (i.e., −*k*/*λ* = 0.4) in the spiral structure displayed in Figs. 1 and 2. This value for *λ*_{M} allows meridional coupling speed (*c*_{CY}) to be computed from Eq. (4.7b), yielding *c*_{CY} = −0.02 m s^{−1} at 55°S. This is approximately one-half the zonal coupling phase speed of *c*_{CX} in Eq. (4.5b). When added to the advection of the coupled wave by the broadscale flow of the ACW [i.e., *u*_{O}*k*_{M}/*λ*_{M} in Eq. (4.7a)], their sum *c*_{ACWY} = −0.04 m s^{−1} again approximately one-half that of *c*_{ACWX} in Eq. (4.5a). Therefore, for the model ACW to remain stable it must propagate poleward at approximately one-half its eastward phase speed.

## 5. Model phase relationship between pycnocline depth and SST anomalies

*T*) anomalies can now be determined by taking the zonal derivative of Eq. (3.2), allowing ∂

*u*

_{A}/∂

*x*to be represented by Eq. (4.2), yielding

*T*=

*T*

_{O}cos(

*k*

_{M}

*x*+

*λ*

_{M}

*y*−

*σ*

_{M}

*t*), then the solution for

*h*can be written

*h*

_{O}=

*T*

_{O}[(

*γαλ*

^{2})/(

*βC*

^{2}

_{A}

*σ*

_{M}

*k*

_{M})]. With

*λ*

_{M}defined by Eq. (4.6) to be −2.2

*k*

_{M}for stable coupled waves, then the ratio of the cosine to the sine term in Eq. (5.2) is approximately 5 to 1. Therefore, since (2

*β*/

*f*)(1/

*λ*

_{M}) > 0 in the Southern Hemisphere, positive

*h*anomalies (i.e., positive downward) are displaced approximately 10° of phase to the east of positive SST anomalies. For magnitudes of

*T*

_{O}= 0.5°C, the magnitude of

*h*

_{O}is approximately 10 m, corresponding to a sea-level height amplitude [i.e.,

*h*

_{O}(

*ρ*′/

*ρ*

_{O})] of approximately 1 cm. This amplitude is consistent with that observed by Jacobs and Mitchell (1996), who also observed no significant lag between SST and sea-level height anomalies associated with the ACW. Moreover, this anomalous vertical displacement of the pycnocline yields an anomalous SST change in Eq. (3.3) of ΔT = −(

*T*/∂

*z*

*h*/2) = 0.1°C. As already discussed in section 3 and as we shall demonstrate in the next section, this influence is too small to significantly influence SST anomalies associated with the ACW.

## 6. Model phase relationships leading to eastward and poleward phase propagation

To understand how phase relationships between model SST, MSW, SLP, and ZSW anomalies are related to eastward and poleward phase propagation of the model ACW in Eqs. (4.4)–(4.7), we employ a graphical approach. This allows us to examine phase relationships as SST anomalies generate ZSW anomalies and as the latter, in turn, generate an anomalous SST tendency through anomalous meridional Ekman heat advection. Moreover, this graphical approach allows us to demonstrate the importance of the beta effect [i.e., (2*β*/*f*) in Eqs. (4.5) and (4.7)] in this coupled interaction.

*T*=

*T*

_{O}sin(

*k*

_{M}

*x*+

*λ*

_{M}

*y*), with

*k*

_{M}= 2

*π*/

*L*

_{XM}and

*λ*

_{M}= −2

*π*/

*L*

_{YM}assumed to be uniform over a schematic ocean extending from 30° to 60°S, 30° to 120°E with no continental boundaries. We evaluate

*L*

_{XM}= 180° longitude (i.e., 1.3 × 10

^{7}m at 45°S),

*L*

_{YM}= 60° latitude (i.e., 6.0 × 10

^{6}m), and

*T*

_{O}= 0.5°C, similar to model scales computed in the previous section, with constant

*f*and

*β*evaluated at 45°S, and model parameters and scales otherwise given in Table 1. These parameters and scales are taken as uniform over the domain, consistent with conditions allowing for the idealized wave solution to Eq. (4.3); yet,

*f*and

*β*are allowed to vary with latitude when we choose. Next, we use the relationship between MSW and SST anomalies in Eq. (A.5) in the appendix to compute the corresponding anomalous SLP(

*P*) wave;that is,

*P*

*x*

*ρ*

_{A}

*f*

^{2}

*α*

*βC*

^{2}

_{A}

*P*= −

*T*

_{O}

*ρ*

_{A}[

*f*

^{2}

*α*/(

*βC*

^{2}

_{A}

*k*

_{M}) cos(

*k*

_{M}

*x*+

*λ*

_{M}

*y*), which is displayed (second panel, Fig. 3) for variable Coriolis parameter (left) and constant Coriolis parameter (right). Since SLP anomalies in Eq. (6.1) vary as the square of the Coriolis parameter (

*f*), their response to SST anomalies is poleward intensified, as discussed by WC98. By comparing left and right panels in Fig. 3, we find the SLP amplitude {i.e.,

*P*

_{O}=

*T*

_{O}

*ρ*

_{A}[

*f*

^{2}

*α*/(

*βC*

^{2}

_{A}

*k*

_{M}) = 6 hPa} at approximately 55°S computed with variable Coriolis parameter. So, maximum SLP anomalies occur along the sloping boundary between warm and cool SST anomalies but are displaced southward by approximately 10° latitude from that centered at 45°S computed on the basis of constant Coriolis parameter. This southward displacement is similar to that observed in Fig. 2.

Negative ZSW anomalies near 45°S, 90°W (i.e., geostrophically inferred from peak SLP anomalies near 55°S) hold the key to understanding how the ACW propagates to the east. To see how this works, we compute an idealized ZSW wave analytically from Eq. (4.2), yielding *u*_{A} = *T*_{O}[*fα*/(*βC*^{2}_{A}*k*_{M})[(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)], which again is displayed (third panel, Fig. 3) for variable Coriolis parameter (left) and constant Coriolis parameter (right). This finds the ZSW wave computed with variable Coriolis parameter occurring nominally along the sloping boundary between warm and cool SST anomalies, but displaced southward and eastward by approximately 5° latitude and longitude, respectively, from that computed on the basis of constant Coriolis parameter. This displacement stems from the cosine term in the expression for the ZSW wave, which is located 90° of phase eastward and poleward of peak SST anomalies. As such, the amplitude of the ZSW wave occurring in phase with the SST wave (i.e., the sine term in the solution) is *u*_{AO} = *T*_{O}[*fα*/(*βC*^{2}_{A}*λ*_{M}/*k*_{M}) = −6 m s^{−1}, while that shifted 90° of phase to the east of the SST wave (i.e., the cosine term in the solution) is *u*_{AO} = −*T*_{O}[*fα*/(*βC*^{2}_{A}*β*/*f*)(1/*k*_{M}) = −1.8 m s^{−1}. The former is responsible for growth of the SST wave in Eq. (4.3), while the latter is responsible for its poleward and eastward phase propagation. This propagation depends upon the beta effect stemming from the coupling through the coefficient 2*β*/*f.*

To see this, we next compute the anomalous SST tendency due strictly to the influence of the ZSW wave upon anomalous meridional Ekman heat advection in Eq. (4.1), yielding Δ*T*/Δ*t* = *T*_{O}[*γα*/(*βC*^{2}_{A}*H*_{P})](*T*/∂*y**k*_{M})[(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)], which is displayed (bottom panel, Fig. 3) for variable Coriolis parameter (left) and constant Coriolis parameter (right). For constant Coriolis parameter, the Δ*T*/Δ*t* wave is in phase with the SST wave in the top panel, in which case the feedback from atmosphere to ocean solely produces a growth tendency for the SST wave. On the other hand, for variable Coriolis parameter, the Δ*T*/Δ*t* wave is displaced eastward and poleward by approximately 18° of phase and, therefore, not only gives the SST wave a growth tendency but causes it to propagate poleward and eastward as well. This propagation stems from the cosine term in the expression for the Δ*T*/Δ*t* wave (i.e., in concert with that in the expression for the ZSW wave) and is due to the poleward intensification of the anomalous SLP response to anomalous SST in Eq. (6.1), the latter translated into a beta effect (i.e., 2*β*/*f*) when the ocean and atmosphere are coupled.

The corresponding anomalous meridional Ekman flows are relatively weak on these interannual timescales, computed to be *υ*_{E} = −(*γ*/*fH*_{P})*u*_{A} = −2.3 × 10^{−3} m s^{−1} for *u*_{A} = −1.8 m s^{−1}. However, over a one-year period even this relatively weak magnitude of Ekman flow will create an anomalous SST tendency of Δ*T* = −*υ*_{E} (*T*/∂*y**t* = 0.5°C displaced 90° of phase eastward and poleward of the original SST wave. Since this portion of the anomalous SST change is the same as the amplitude of the SST wave, this magnitude of anomalous Ekman flow can account for eastward and poleward phase propagation of the ACW. The total anomalous SST change (i.e., due to anomalous meridional Ekman heat transport both in phase and shifted 90° of phase eastward and poleward of original SST wave) is approximately 2.1°C, with that portion in phase with the SST wave (i.e., 2.0°C) accounting for a growth tendency in the amplitude of the wave. This growth tendency is balanced by the anomalous sensible-plus-latent heat loss to the atmosphere, thereby keeping the wave amplitude stable.

*υ*

_{E}

*T*=

*T*

_{O}sin(

*k*

_{M}

*x*+

*λ*

_{M}

*y*) and following the sequence of interactions given above to arrive at the meridional Ekman flow wave

*υ*

_{E}= −(

*γ*/

*fH*

_{P})

*u*

_{A}= −

*T*

_{O}[

*αγ*/(

*βH*

_{P}

*C*

^{2}

_{A}

*k*

_{M})[(2

*β*/

*f*) cos(

*k*

_{M}

*x*+

*λ*

_{M}

*y*) −

*λ*

_{M}sin(

*k*

_{M}

*x*+

*λ*

_{M}

*y*)]. This has a constituent common with

*T*(i.e., the sine term), which together yield a net poleward eddy heat flux EHF =

*ρ*

_{O}

*C*

_{PO}

*H*

_{P}2

*L*

_{XM}

*υ*

_{E}

*T*

*ρ*

_{O}

*C*

_{PO}

*T*

^{2}

_{O}

*αγ*

*βC*

^{2}

_{A}

*λ*

_{M}

*k*

_{M}

*L*

_{XM}

^{13}

^{15}watts), established indirectly from budget considerations (Niiler 1992). The portion of the

*υ*

_{E}wave responsible for this net poleward eddy heat flux is that leading to the anomalous meridional Ekman heat advection that maintains the ACW against dissipation, not that leading to eastward and poleward phase propagation. This is also the portion of the

*υ*

_{E}wave that leads to the tilt in the orientation of the covarying SST and SLP anomaly patterns in Fig. 1. As such, the spiral structure in the ACW anomaly patterns in Fig. 2 can be considered a signature not only of the coupling between ocean and atmosphere, but also of the net poleward eddy heat flux in the upper ocean by the ACW, weak as it is.

## 7. Numerical simulation of the coupled ACW

To demonstrate that this simple coupled model for the ACW is capable of simulating most of its observed wave characteristics in Figs. 1 and 2, we construct a numerical version of this coupled model by linking the numerical version of the global equilibrium climate model of WC98 with a numerical version of the anomalous heat budget model for the near-surface mixed layer in the Southern Ocean given by Eq. (4.1). We extend the latter by including both meridional and zonal components of the mean background flow of the ACC in the Southern Ocean. The finite difference analogue for this coupled model is constructed on a 2° by 5° latitude–longitude grid with realistic coastline topography of the Southern Ocean and spatial-varying estimates of background surface currents and SST gradients. We choose numerical model parameters to be the same as analytic model parameters given above with one exception; that is, *K*^{−1} is evaluated with an *e*-folding timescale of 1 year. This is larger than 3 months utilized in the analytic model required to keep the numerical coupled model stable. In the analytic solution to Eq. (4.3), increasing *K*^{−1} in this way would have had no influence upon the eastward phase speed of the model ACW, but it should have increased the poleward phase speed by a factor of 4. Yet, as we shall see, the poleward phase speed in the numerical model is that predicted by the analytic coupled model. This suggests that increasing *K*^{−1} to balance some numerical noise and/or dispersion preserved the necessary thermodynamic balance. We initialize this numerical coupled model with an observed pattern of SST anomalies from April 1987 (Fig. 2) and determine whether it can simulate the space–time evolution of the ACW and maintain its amplitude against dissipation. We compare the observed ACW to the model ACW after 13 years of model integration, after which the model generates its own repeating evolution of the ACW independent of the April 1987 initial conditions.

In this numerical model, the mean surface flow is composed of both geostrophic and Ekman components. The former is made relative to 2000 db, computed using the dynamic method from the available hydrographic data from the period 1950–89, gridded onto the model 2° lat. by 5° long. grid. The annual-mean Ekman flow is computed from 15 years (1982–96) of surface wind stress estimates from the NCEP Reanalysis Project (Kalnay et al. 1996) utilizing a mean near-surface mixed-layer depth scale in autumn–winter of H_{P} = 100 m, with vertical-averaged Ekman flow directed 90° to the left of the surface wind stress in the Southern Hemisphere. When this Ekman flow is added onto the geostrophic flow, the total surface flow (Fig. 4) displays the relatively narrow core of the eastward flowing ACC (i.e., with speeds ranging from 0.05 to 0.12 m s^{−1}) extending around the Southern Ocean between 45° and 60°S, embedded in the general eastward broadscale flow of the ACC (i.e., with speeds of approximately 0.05 m s^{−1}) occurring throughout the Southern Hemisphere from approximately 40° to 70°S. The Ekman flow produces a broadscale equatorward flow directed away from Antarctica over most of the Southern Ocean.

Given this broadscale surface flow pattern, we conduct two numerical experiments: one in the absence of coupling and another in its presence. In the absence of coupling (not shown), anomalous SST initial conditions from April 1987 are advected eastward by the broadscale flow of the ACC and equatorward by mean Ekman flow in all three oceans, consistent with Fig. 4, with amplitudes decreasing over time to insignificance after 6–8 years. Even so, a remnant of the initial SST signal is advected by the ACC completely around the Southern Ocean after 12–14 years, nearly twice the time it takes (i.e., approximately 8 years) for the observed ACW to circle the globe (White and Peterson 1996). This indicates that the broadscale flow of the ACC is too slow (i.e., by nearly a factor of 2) to account for the eastward propagation of the ACW.

In the presence of coupling and beginning 13 years after initiation of model integration, an animation of SST anomalies from the coupled model (Fig. 5) reveals a model ACW with a 4-yr period and the required two wavelengths circling the globe in the Southern Ocean, similar to that observed in Fig. 1. This arises because the zonal coupling speed in Eq. (4.5b) approximately doubles the broadscale flow of the ACC. Moreover, magnitudes of the interannual SST anomalies are maintained against dissipation and dispersion of the initial SST signal. This occurs in response to the balance between anomalous meridional Ekman heat advection and anomalous sensible-plus-latent heat loss to the atmosphere expressed in Eq. (4.3). In addition, model SST anomalies in the Indian and Pacific sectors of the Southern Ocean display the characteristic elongation and tilt of their orientation into the southwest–northeast direction observed in Fig. 1 and predicted by Eq. (4.6). It also displays significant poleward propagation, where SST anomalies in the vicinity of New Zealand can be seen propagating eastward and poleward toward Drake Passage as observed in Fig. 1. The speed of this southward propagation is approximately 0.04 m s^{−1}, the same as that evaluated analytically in Eq. (4.7a). However, the model simulation differs significantly from that observed in the Atlantic sector of the Southern Ocean where the mean geostrophic flow of the ACC in Fig. 4 is underresolved, unable to simulate the strong advection of observed SST anomalies through and east of Drake Passage observed in Fig. 1.

We also integrated a version of the numerical coupled model where in the model atmosphere the Coriolis parameter southward of 30°S was set to a constant value (i.e., *f* = *f*_{O} at 45°S) except where its meridional gradient is taken explicitly. From the analytical model in connection with discussion of Fig. 3, this step would be expected to eliminate eastward and southward coupling speeds but retain maintenance of the initial signal against dissipation. Indeed, this is what happened, with the model ACW maintained against dissipation and characterized by the southwest–northeast tilt in the alignment of the anomalies but propagating eastward only with the broadscale flow of the ACC in Fig. 4. This model experiment demonstrates that the fundamental physics of coupling in the numerical model is similar to that in the analytical model, with eastward and poleward propagation stemming from the *f*-squared dependency in the SLP response to SST anomalies in Eq. (6.1), which translates into a beta effect when the ocean and atmosphere are coupled.

We can also demonstrate the effect of this *f*-squared dependency in the numerical simulation by examining the phase relationship between covarying SST and SLP anomalies in the model ACW (Fig. 6) as was done with observations in Fig. 2 and with the analytical coupled model in Fig. 3. With the latter we showed that the impact of the *f*-squared dependency upon SLP anomalies along the sloping boundary between warm and cool SST anomalies is to displace them southward from the latitude line connecting maximum SST anomalies, as observed in Fig. 2. This is also observed in the snapshot from the numerical simulation in Fig. 6.

This numerical coupled model is also capable of simulating the global spiral structure in covarying SST and SLP anomal patterns associated with the ACW observed in the polar-stereographic projection in Fig. 6. This indicates that the global spiral structure is a signature for coupling between ocean and atmosphere. It is a manifestation of the southwest–northeast tilt in the orientation of the anomalies in the ACW in Eq. (4.6) and derives from the balance between anomalous meridional Ekman heat advection and anomalous sensible-plus-latent heat loss to the atmosphere. It does not depend upon the beta effect due to coupling. The latter is only responsible for eastward and poleward coupling speeds.

## 8. Discussion and conclusions

In the present study, covarying SST and SLP anomaly patterns in the ACW were observed to take on a global spiral structure in the Southern Ocean around Antarctica (Fig. 2). Earlier, White and Peterson (1996) had found poleward (equatorward) MSW anomalies to be aligned with warm (cool) SST anomalies in the ACW. Recently, WC98) had explained this association in terms of a balance between anomalous meridional planetary vorticity advection and SST-induced anomalous vortex stretching of the lower troposphere (see appendix). This indicated that SLP anomalies should straddle the boundary between warm and cool SST anomalies and be poleward intensified. Here, we observed this to be so, with the resulting pattern of SLP anomalies yielding ZSW anomalies that are displaced eastward of peak SST anomalies. Therefore, we hypothesized that associated anomalous meridional Ekman heat advection, and/or anomalous vertical mixing and anomalous sensible-plus-latent heat loss with which it is in phase, induces an anomalous SST tendency that accounts for the eastward propagation of the ACW.

To test this hypothesis, we constructed analytic and numerical models of this coupled system, with MSW anomalies in balance with underlying SST anomalies (WC98) and with associated ZSW anomalies providing a positive feedback to the ocean, forcing the anomalous SST tendency with anomalous meridional Ekman advection of the background SST field. The analytic model yielded the dispersion relation for free waves of the coupled system [i.e., Eqs. (4.4)–(4.7)], while the numerical model yielded a simulation for the evolution of the model ACW in the presence of spatially varying background SST, geostrophic plus Ekman flow, and coastline topography (Figs. 5 and 6). Both modeling efforts established the presence of eastward and poleward phase speeds of approximately 0.09 cm s^{−1} and 0.04 cm s^{−1}, respectively, with zonal wavelengths of 180° longitude (i.e., 1.1 × 10^{7} m at 55°S) and meridional wavelengths of 50° latitude (i.e., 5.2 × 10^{6} m), all similar to those observed (summarized in Table 1). In the analytic model, these speeds and wavelengths stem from requiring the coupled system to fluctuate with a period of 4.5 yr observed in the ACW by Peterson and White (1996). In the numerical model, a nominal 4–5 yr periodicity was generated intrinsically, with both frequency and wavenumber arising together from model scale and parameter estimation. We conclude that this particular frequency and wavenumber set arises as a global normal mode of the particular thermodynamics contained in this coupled ocean–atmosphere model of the Southern Ocean.

The thermodynamics of coupling leading to the eastward and poleward phase propagation of the ACW, and to its maintenance against dissipation, were explained through examination of the phase relationship between SST and SLP anomalies in both observations (Fig. 2) and numerical simulations (Fig. 6). These relationships are summarized in the schematic diagram given in Fig. 7, which is a significant improvement over that lacking the spiral structure introduced initially by White and Peterson (1996). In this schematic diagram, maximum SLP anomalies occur along the sloping boundary between warm and cool SST anomalies, displaced poleward from the latitude connecting maximum SST anomalies. In the analytic coupled model, we found this displacement arising from the *f*-squared dependency in the anomalous SLP response to SST anomalies, yielding poleward intensification of the response. The spiral structure in this SLP anomaly pattern acts to bring anomalous westward (eastward) ZSW anomalies into phase with warm (cool) SST anomalies, allowing anomalous meridional Ekman heat advection to balance the anomalous sensible-plus-latent heat loss to the overlying atmosphere and thereby maintain the ACW against this form of dissipation. The poleward intensification of this SLP anomaly pattern yields a portion of the ZSW anomalies that drives anomalous meridional Ekman heat advection, and its corresponding anomalous SST tendency, eastward and poleward of SST anomalies. This accounts for the eastward and poleward propagation of the covarying SST and SLP anomalies.

We also demonstrated that the model ACW is associated with a net poleward eddy heat flux from subtropical to polar latitudes, deriving from the coincidence of the anomalous poleward (equatorward) Ekman flow and warm (cool) SST anomalies in the near-surface mixed layer of the Southern Ocean. Yet, it is two orders of magnitude smaller than the net poleward transport of heat by the ocean (Niiler 1992). We demonstrated that it derives from the same coupled interactions that give rise to the maintenance of the ACW against dissipation, both of which are linked to the spiral structure in the global SST and SLP anomaly patterns in Fig. 7. Therefore, the latter becomes a signature not only for the coupling between ocean and atmosphere but also for the net poleward eddy flux of heat associated with the ACW, albeit weak.

Even though this coupling is restricted to the near-surface mixed layer of the ocean, the analytical model found downward displacements in the main pycnocline nearly in phase with warm SST anomalies, displaced only 10° of phase to the east. This is consistent with that observed by Jacobs and Mitchell (1996), where pycnocline displacements were inferred from altimetric sea-level height anomalies. But the interpretation given here is different from that given by Jacobs and Mitchell (1996). They proposed that anomalous vertical Ekman pumping and associated meridional geostrophic heat advection is the dominant advective influence upon the anomalous SST tendency. In our analytic model, we found this mechanism too weak; instead, we found anomalous meridional Ekman heat advection responsible for the feedback from atmosphere to ocean.

While the numerical version of the coupled model was able to simulate much of the observed evolution of the ACW in the Indian and Pacific sectors of the Southern Ocean, it was unable to simulate the observed jetlike advection of the ACW through Drake Passage and east of there in the Atlantic sector (Fig. 1). We ascribed this to the inability of the mean geostrophic flow (0/2000 db) used in the model to resolve the jetlike core of the ACC throughout much of this domain. To simulate the behavior of the ACW in this domain correctly will require high-resolution modeling of the ACC there. Moreover, the numerical model did not simulate maximum interannual variability in SST and SLP in the vicinity of the sea ice edge in the Pacific sector of the Southern Ocean (Peterson and White 1998). This may require the inclusion of model interactions between ocean, atmosphere, and sea ice anomalies. Furthermore, this model did not contain remote wind stress forcing arising from meridional atmospheric teleconnections (i.e., quasi-stationary Rossby waves in the atmosphere), which may have projected tropical ENSO signals into the Southern Ocean domain (Karoly 1989).

A number of assumptions need to be confirmed before the thermodynamics expressed in this coupled model can be accepted. Principal among these is the assumption that the response of the atmosphere to anomalous SST forcing produces MSW and ZSW anomalies that follow the geostrophic winds in the lower troposphere. In this study, we have ignored the Ekman deflection of geostrophic wind anomalies in the near-surface planetary boundary layer of the lower troposphere. This might have yielded a component of anomalous Ekman heat advection in the near-surface mixed layer that significantly alters the brand of coupling between ocean and atmosphere presented here. Moreover, for the sake of simplicity this model assumed permanent winter conditions, when the near-surface mixed layer coincides with the top of the main pycnocline and when SST anomalies instigate midlevel diabatic heating in the overlying troposphere. We might expect that, had the numerical model included the annual cycle of atmospheric feedback and mixed layer dynamics, the resulting ACW would have exhibited at least a winter dominance and at most a chaotic element that would have transformed the periodic behavior of the model ACW into the quasiperiodic behavior of the limit cycle.

Finally, while anomalous meridional Ekman heat advection (i.e., linear relation to ZSW anomalies) was chosen here to represent the principal feedback mechanism from the atmosphere to the ocean, this feedback could also have derived from anomalous wind mixing at the base of the near-surface mixed layer and/or from anomalous sensible-plus-latent heat flux, both of which are also linearly related to ZSW anomalies. Also, anomalous cloud cover might have affected the anomalous SST tendency through its influence on anomalous air–sea radiation fluxes. These additional mechanisms, and others, will need to be sorted out quantitatively before we understand fully the coupling mechanisms responsible for the ACW.

## Acknowledgments

This research was partially supported by the Office of Global Programs of NOAA (NOAA NA 47GP0-188) in concert with the Scripps–Lamont Consortium on the Ocean’s Role in Climate, and by the National Science Foundation, Division of Ocean Sciences (OCE.96-33474). Warren White is also supported by the National Aeronautics and Space Administration (NASA) under Contract NA27GPO-539. Our thanks extend to Ted Walker, who provided most of the computational and visualization support.

## REFERENCES

Gill, A. E., 1982:

*Atmosphere–Ocean Dynamics.*Int. Geophys. Ser., Vol. 30 Academic Press, 662 pp.Hoskins, B. J., and D. Karoly, 1981: The steady linear response of a spherical atmosphere to thermal and orographic forcing.

*J. Atmos. Sci.,***38,**1179–1196.Jacobs, G. A., and J. L. Mitchell, 1996: Ocean circulation variations associated with the Antarctic Circumpolar Wave.

*Geophys. Res. Lett.,***23,**2947–2950.Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project.

*Bull. Amer. Meteor. Soc.,***77,**437–471.Karoly, D. J., 1989: Southern Hemisphere circulation features associated with El Niño–Southern Oscillation events.

*J. Climate,***2,**1239–1252.Kaylor, R. E., 1977: Filtering and decimation of digital time series. Inst. Phys. Sci. and Technol. University of Maryland Tech. Note BN 850, 14 pp.

Kraus, E. B., and J. A. Businger, 1994:

*Atmosphere–Ocean Interaction.*Oxford University Press, 362 pp.Miller, A. J., D. R. Cayan, T. P. Barnett, N. E. Graham, and J. M. Oberhuber, 1994: Interdecadal variability of the Pacific Ocean: Model response to observed heat flux and wind stress anomalies.

*Climate Dyn.,***10,**287–302.Namias, J., 1985: Remarks of the potential for long-range forecasting.

*Bull. Amer. Meteor. Soc.,***66,**165–173.Niiler, P. P., 1992: The ocean circulation.

*Climate System Modeling,*K. Trenberth, Ed., Cambridge University Press, 117–144.Olbers, D., V. Gouretski, G. Seiss, and J. Schröter, 1992:

*Hydrographic Atlas of the Southern Ocean.*Alfred Wegener Inst., 82 plates, 17 pp.Palmer, T. N., and Z. Sun, 1985: A modeling and observational study of the relationship between sea surface temperature in the north-west Atlantic and the atmospheric general circulation.

*Quart. J. Roy. Meteor. Soc.,***111,**947–975.Pedlosdky, J., 1996:

*Ocean Circulation Theory.*Springer, 453 pp.Peterson, R., and W. B. White, 1997: Slow oceanic teleconnections linking tropical ENSO and the Antarctic Circumpolar Wave.

*J. Geophys. Res.,*in press.Robertson, A. W., 1996: Interdecadal variability over the North Pacific in a multi-century climate simulation.

*Climate Dyn.,***12,**227–241.White, W. B., and R. Peterson, 1996: An Antarctic circumpolar wave in surface pressure, wind, temperature, and sea ice extent.

*Nature,***380,**699–702.

## APPENDIX

### Lower-Tropospheric Response to Extratropical SST Anomalies

*βυ*

_{A}

*fD*

_{A}

*β*is the meridional derivative of the Coriolis parameter

*f,*and

*υ*

_{A}and −

*D*

_{A}are meridional wind and convergence anomalies, respectively, both vertically averaged over the lower troposphere. Second, they found the low-level convergence anomalies associated with warm SST anomalies over most of the global ocean.

*Q*

*α*

_{A}

*T,*

*Q*is the midlevel diabatic heating anomaly,

*T*is the SST anomaly, and

*α*

_{A}is an empirically derived coefficient equating the two variables. Moreover, they found these year-to-year changes in winter midlevel diabatic heating arising from changes in the midlevel latent heat release stemming from year-to-year changes in aggregate synoptic storm activity. This was consistent with an earlier hypothesis proposed by Namias (1985) that extratropical SST anomalies destabilize the overlying troposphere during the autumn–winter, thereby modifying the location and/or intensity of aggregate synoptic storm activity.

*w*

_{A}

*Q*

*g*

*ρ*

_{A}

*N*

^{2}

*C*

_{PA}

*H*

_{A}

*w*

_{A}is the anomalous vertical velocity at midlevel;

*N*

^{2}is the mean Väisälä frequency at midlevel [i.e., (

*g*/

*z*

*H*

_{A}is the mean height of the lower troposphere;

*ρ*

_{A}is the mean density of air; and

*C*

_{PA}is the specific heat of air. By applying conservation of mass to the lower troposphere (i.e., where

*w*

_{A}/

*H*

_{A}= −

*D*

_{A}= −∂

*u*

_{A}/∂

*x*− ∂

*υ*

_{A}/∂

*y*), the potential temperature balance in Eq. (A.3) can be linked with the empirical relationship in Eq. (A.2), yielding

*D*

_{A}

*u*

_{A}

*x*

*υ*

_{A}

*y*

*α*

*C*

^{2}

_{A}

*T,*

*C*

_{A}is the internal gravity wave speed (i.e.,

*NH*

_{A}= 50 m s

^{−1}) and

*α*is a coefficient [i.e., (

*α*

_{A}

*g*)/(

*ρ*

_{A}

*C*

_{PA}

^{−3}m

^{2}s

^{−3}°C

^{−1}]. This relationship is consistent with the association observed between warm SST anomalies and low-level anomalous convergence.

*u*

_{A}) in the lower troposphere to be computed from knowledge of extratropical SST anomalies (

*T*) and vertical-averaged meridional wind anomalies (

*υ*

_{A}). The latter can be related directly to SST anomalies by linking anomalous mass and potential temperature conservation in Eq. (A.4) to anomalous vorticity conservation in Eq. (A.1), yielding

*υ*

_{A}

*fα*

*βC*

^{2}

_{A}

*T.*

*Quart. J. Roy. Meteor. Soc.*) demonstrated that this latter diagnostic relationship allows both magnitude and phase of global SLP anomaly patterns [i.e.,

*υ*

_{A}=

*f*

^{−1}∂

*P*/∂

*x,*where the vertical-averaged low-level pressure anomaly (

*P*) is assumed equivalent to the SLP anomaly] to be diagnosed from observed global SST anomaly patterns for the 8 years of analysis from 1985 to 1992. This included the diagnoses of MSW anomalies (i.e., assumed equivalent to

*υ*

_{A}) from SST anomalies in the ACW observed by White and Peterson (1996).

Interannual SST and SLP anomalies for April 1987 plotted contemporaneously on a polar-stereographic projection centered at the South Pole and extending to 30°S. Anomalies were formed from SST and SLP variables provided by the NCEP reanalysis project for the 15-yr period from 1982 to 1996 (Kalnay et al. 1996). The zonal mean was removed from both SST and SLP anomalies and spatial smoothing applied to enhance the structure of the ACW. Contours of SST are given in the color bar, ranging from ±1.0°C. Contours of SLP are 0.5 hPa, with positive anomalies given by thin solid lines and negative anomalies given by thin dashed lines, with the zero contour given by a thick solid line.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Interannual SST and SLP anomalies for April 1987 plotted contemporaneously on a polar-stereographic projection centered at the South Pole and extending to 30°S. Anomalies were formed from SST and SLP variables provided by the NCEP reanalysis project for the 15-yr period from 1982 to 1996 (Kalnay et al. 1996). The zonal mean was removed from both SST and SLP anomalies and spatial smoothing applied to enhance the structure of the ACW. Contours of SST are given in the color bar, ranging from ±1.0°C. Contours of SLP are 0.5 hPa, with positive anomalies given by thin solid lines and negative anomalies given by thin dashed lines, with the zero contour given by a thick solid line.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Interannual SST and SLP anomalies for April 1987 plotted contemporaneously on a polar-stereographic projection centered at the South Pole and extending to 30°S. Anomalies were formed from SST and SLP variables provided by the NCEP reanalysis project for the 15-yr period from 1982 to 1996 (Kalnay et al. 1996). The zonal mean was removed from both SST and SLP anomalies and spatial smoothing applied to enhance the structure of the ACW. Contours of SST are given in the color bar, ranging from ±1.0°C. Contours of SLP are 0.5 hPa, with positive anomalies given by thin solid lines and negative anomalies given by thin dashed lines, with the zero contour given by a thick solid line.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

A schematic representation is given of the model phase relationships between anomalous SST, SLP, ZSW, and SST tendency, extending over the domain from 30° to 60°S, 30° to 120°E. Top panel: The model SST wave is assumed to be proportional to sin(*k*_{M}*x* + *λ*_{M}*y*), with wavenumbers and amplitudes given in the text. Second panel: The model SLP wave in response to the SST wave then is proportional to −*f*^{2} cos(*k*_{M}*x* + *λ*_{M}*y*), with amplitudes given in the text. Third panel: The model ZSW wave is related to the SLP wave through the geostrophic relationship, proportional to (*f*/*β*)[(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)], with amplitudes given in the text. Bottom panel: Finally, the model SST tendency wave due only to anomalous meridional Ekman heat advection, driven by the model ZSW wave, becomes proportional to [(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)]*β*^{−1}, with amplitudes given in the text. In the left column are results based upon a variable Coriolis parameter in the atmosphere; in the right column are results based upon a constant Coriolis parameter where terms within brackets containing *β* have been set to zero. Negative anomalies are shaded; positive anomalies are unshaded. In each panel, amplitudes have been weighted by latitude with a Gaussian filter that allows zonal and meridional phase shifts between right and left panels to be discerned more easily. The sloping heavy black line is the same in each panel, used as a reference to discern shifts in phase between right and left panels. See text describing how model waves associated with each variable can be followed from one panel to the next and how wave amplitudes are computed.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

A schematic representation is given of the model phase relationships between anomalous SST, SLP, ZSW, and SST tendency, extending over the domain from 30° to 60°S, 30° to 120°E. Top panel: The model SST wave is assumed to be proportional to sin(*k*_{M}*x* + *λ*_{M}*y*), with wavenumbers and amplitudes given in the text. Second panel: The model SLP wave in response to the SST wave then is proportional to −*f*^{2} cos(*k*_{M}*x* + *λ*_{M}*y*), with amplitudes given in the text. Third panel: The model ZSW wave is related to the SLP wave through the geostrophic relationship, proportional to (*f*/*β*)[(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)], with amplitudes given in the text. Bottom panel: Finally, the model SST tendency wave due only to anomalous meridional Ekman heat advection, driven by the model ZSW wave, becomes proportional to [(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)]*β*^{−1}, with amplitudes given in the text. In the left column are results based upon a variable Coriolis parameter in the atmosphere; in the right column are results based upon a constant Coriolis parameter where terms within brackets containing *β* have been set to zero. Negative anomalies are shaded; positive anomalies are unshaded. In each panel, amplitudes have been weighted by latitude with a Gaussian filter that allows zonal and meridional phase shifts between right and left panels to be discerned more easily. The sloping heavy black line is the same in each panel, used as a reference to discern shifts in phase between right and left panels. See text describing how model waves associated with each variable can be followed from one panel to the next and how wave amplitudes are computed.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

A schematic representation is given of the model phase relationships between anomalous SST, SLP, ZSW, and SST tendency, extending over the domain from 30° to 60°S, 30° to 120°E. Top panel: The model SST wave is assumed to be proportional to sin(*k*_{M}*x* + *λ*_{M}*y*), with wavenumbers and amplitudes given in the text. Second panel: The model SLP wave in response to the SST wave then is proportional to −*f*^{2} cos(*k*_{M}*x* + *λ*_{M}*y*), with amplitudes given in the text. Third panel: The model ZSW wave is related to the SLP wave through the geostrophic relationship, proportional to (*f*/*β*)[(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)], with amplitudes given in the text. Bottom panel: Finally, the model SST tendency wave due only to anomalous meridional Ekman heat advection, driven by the model ZSW wave, becomes proportional to [(2*β*/*f*) cos(*k*_{M}*x* + *λ*_{M}*y*) − *λ*_{M} sin(*k*_{M}*x* + *λ*_{M}*y*)]*β*^{−1}, with amplitudes given in the text. In the left column are results based upon a variable Coriolis parameter in the atmosphere; in the right column are results based upon a constant Coriolis parameter where terms within brackets containing *β* have been set to zero. Negative anomalies are shaded; positive anomalies are unshaded. In each panel, amplitudes have been weighted by latitude with a Gaussian filter that allows zonal and meridional phase shifts between right and left panels to be discerned more easily. The sloping heavy black line is the same in each panel, used as a reference to discern shifts in phase between right and left panels. See text describing how model waves associated with each variable can be followed from one panel to the next and how wave amplitudes are computed.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Global distribution of the surface geostrophic plus Ekman flow over the Southern Hemisphere from 30° to 70°S. Vectors indicate flow direction and intensity, while contours indicate speed of the zonal component. The geostrophic component is made relative to 2000 db, computed from the long-term annual-mean geopotential anomaly computed from all available hydrograhic data collected from 1950 to the present. The vertical-averaged Ekman component (i.e., directed 90° to the left of the surface wind stress in the Southern Hemisphere) was computed from the long-term annual mean wind stress (i.e., 1982–96) obtained from the NCEP Reanalysis Project (Kalnay et al. 1996), acting over the winter near-surface mixed layer depth scale of 100 m. Contour intervals are 0.02 m s^{−1}. Shading is for effect.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Global distribution of the surface geostrophic plus Ekman flow over the Southern Hemisphere from 30° to 70°S. Vectors indicate flow direction and intensity, while contours indicate speed of the zonal component. The geostrophic component is made relative to 2000 db, computed from the long-term annual-mean geopotential anomaly computed from all available hydrograhic data collected from 1950 to the present. The vertical-averaged Ekman component (i.e., directed 90° to the left of the surface wind stress in the Southern Hemisphere) was computed from the long-term annual mean wind stress (i.e., 1982–96) obtained from the NCEP Reanalysis Project (Kalnay et al. 1996), acting over the winter near-surface mixed layer depth scale of 100 m. Contour intervals are 0.02 m s^{−1}. Shading is for effect.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Global distribution of the surface geostrophic plus Ekman flow over the Southern Hemisphere from 30° to 70°S. Vectors indicate flow direction and intensity, while contours indicate speed of the zonal component. The geostrophic component is made relative to 2000 db, computed from the long-term annual-mean geopotential anomaly computed from all available hydrograhic data collected from 1950 to the present. The vertical-averaged Ekman component (i.e., directed 90° to the left of the surface wind stress in the Southern Hemisphere) was computed from the long-term annual mean wind stress (i.e., 1982–96) obtained from the NCEP Reanalysis Project (Kalnay et al. 1996), acting over the winter near-surface mixed layer depth scale of 100 m. Contour intervals are 0.02 m s^{−1}. Shading is for effect.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Time sequence of interannual SST anomalies from the numerical coupled model of the ACW in the Southern Ocean, extending from 30°S to 70°S and zonally around the globe. The time sequence displays each January and June for 4 years beginning 13 years after the start of model integration from SST initial conditions given for April 1987 in Fig. 2. In each figure, the zonal mean was removed to enhance the structure of the ACW. Contours are given by the color bar ranging from ±1.0°C.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Time sequence of interannual SST anomalies from the numerical coupled model of the ACW in the Southern Ocean, extending from 30°S to 70°S and zonally around the globe. The time sequence displays each January and June for 4 years beginning 13 years after the start of model integration from SST initial conditions given for April 1987 in Fig. 2. In each figure, the zonal mean was removed to enhance the structure of the ACW. Contours are given by the color bar ranging from ±1.0°C.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Time sequence of interannual SST anomalies from the numerical coupled model of the ACW in the Southern Ocean, extending from 30°S to 70°S and zonally around the globe. The time sequence displays each January and June for 4 years beginning 13 years after the start of model integration from SST initial conditions given for April 1987 in Fig. 2. In each figure, the zonal mean was removed to enhance the structure of the ACW. Contours are given by the color bar ranging from ±1.0°C.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Model SST and SLP anomalies for June in year 14 of the numerical model integration plotted contemporaneously on a polar stereographic projection centered at the South Pole and extending to 30°S. Anomalies were formed intrinsically in the model, but the zonal mean was removed to enhance the structure of the ACW. Contours of SST anomalies are given in the color bar, ranging from ±1.0°C. Contours of SLP anomalies are at 0.5 hPa intervals, with positive anomalies given by thin solid lines and negative anomalies given by thin dashed lines, with the zero contour given by a thick black line.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Model SST and SLP anomalies for June in year 14 of the numerical model integration plotted contemporaneously on a polar stereographic projection centered at the South Pole and extending to 30°S. Anomalies were formed intrinsically in the model, but the zonal mean was removed to enhance the structure of the ACW. Contours of SST anomalies are given in the color bar, ranging from ±1.0°C. Contours of SLP anomalies are at 0.5 hPa intervals, with positive anomalies given by thin solid lines and negative anomalies given by thin dashed lines, with the zero contour given by a thick black line.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Model SST and SLP anomalies for June in year 14 of the numerical model integration plotted contemporaneously on a polar stereographic projection centered at the South Pole and extending to 30°S. Anomalies were formed intrinsically in the model, but the zonal mean was removed to enhance the structure of the ACW. Contours of SST anomalies are given in the color bar, ranging from ±1.0°C. Contours of SLP anomalies are at 0.5 hPa intervals, with positive anomalies given by thin solid lines and negative anomalies given by thin dashed lines, with the zero contour given by a thick black line.

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Schematic diagram of interannual SST anomalies (red: warm, blue: cool) and SLP anomalies (black lines: high; white lines: low) observed in Figs. 1 and 2, representing their spatial patterns and phase relationships in the ACW over the Southern Ocean. Schematic meridional geostrophic wind anomalies are given by the broad arrows; anomalous sea ice extent is given by the gray lines. Information about sea ice extent is based upon a 13-yr dataset analyzed by White and Peterson (1996).

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Schematic diagram of interannual SST anomalies (red: warm, blue: cool) and SLP anomalies (black lines: high; white lines: low) observed in Figs. 1 and 2, representing their spatial patterns and phase relationships in the ACW over the Southern Ocean. Schematic meridional geostrophic wind anomalies are given by the broad arrows; anomalous sea ice extent is given by the gray lines. Information about sea ice extent is based upon a 13-yr dataset analyzed by White and Peterson (1996).

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Schematic diagram of interannual SST anomalies (red: warm, blue: cool) and SLP anomalies (black lines: high; white lines: low) observed in Figs. 1 and 2, representing their spatial patterns and phase relationships in the ACW over the Southern Ocean. Schematic meridional geostrophic wind anomalies are given by the broad arrows; anomalous sea ice extent is given by the gray lines. Information about sea ice extent is based upon a 13-yr dataset analyzed by White and Peterson (1996).

Citation: Journal of Physical Oceanography 28, 12; 10.1175/1520-0485(1998)028<2345:TACWAB>2.0.CO;2

Summary of scales and model parameters.