1. Introduction
The Antarctic Circumpolar Wave (ACW) is identified by Reference White and PetersonWhite and Peterson (1996) as anomalies in sea-level pressure (SLP), meridional wind stress (MWS), sea-surface temperature (SST) and sea-ice extent (SIE) over the Southern Ocean, propagating eastward with a period of 4−5 years and taking 8−10 years to encircle the globe. Thus, it is a wavenumber-2 phenomenon. At about the same time, Reference Jacobs and MitchellJacobs and Mitchell (1996) found a similar ACW mode in anomalous sea-surface height from TOPEX/POSEIDON altimeter observations. The datasets used by Reference White and PetersonWhite and Peterson (1996) consist of in situ and satellite measurements for SST and SIE, and 1985−94 European Centre for Medium-range Weather Forecasts (ECMWF) operational analyses for SLP and MWS. Utilizing 1979−94 ECMWF re-analysis (ERA) (a longer dataset as compared to the period covered by Reference White and PetersonWhite and Peterson (1996)), Reference Bonekamp, Sterl and KomenBonekamp and others (1999) identified both propagating and standing oscillations in SLP and MWS anomalies. For the period 1985−94, coinciding with the period studied by Reference White and PetersonWhite and Peterson (1996), eastward propagation is present in more or less the same way as in the operational analysis used by Reference White and PetersonWhite and Peterson (1996). However, during the earlier period 1979−84, the eastward propagation is absent, and the standing oscillations exist instead.
Unlike the wavenumber-2 ACW mode described by Reference White and PetersonWhite and Peterson (1996), coupled climate-model studies by Reference Christoph, Barnett and RoecknerChristoph and others (1998) and Reference Cai, Baines and GordonCai and others (1999) found a dominant wavenumber-3 pattern in atmospheric, oceanic and sea-ice variables. Although Reference White and PetersonWhite and Peterson (1996) reported the phase-locking of observed ACW components (i.e. all anomalies propagate eastward at the same rate), the atmospheric anomalies in the studies of Reference Christoph, Barnett and RoecknerChristoph and others (1998) and Reference Cai, Baines and GordonCai and others (1999) showed little or no eastward propagation, while the oceanic anomalies propagated eastward. The period of their oceanic ACW mode is 4−5 years, but it takes about 12−16 years for an anomaly to encircle the globe.
Several studies have tried to understand the dynamics of the ACW. Reference Qiu and JinQiu and Jin (1997) and Reference Peterson and WhitePeterson and White (1998) concluded that the ACW is a mode of the coupled atmosphere-ocean sea-ice system and that atmosphere-ocean interaction plays a dominant role in maintaining the ACW from dissipation as it travels around the globe. In the simulations of Reference Qiu and JinQiu and Jin (1997) with a simple coupled atmosphere-ocean model, the atmosphere reacts to a warm SST anomaly with a high (low) SLP anomaly to the east (west) of it (see Reference Qiu and JinQiu and Jin, 1997, plate 2). The corresponding wind-stress anomalies deepen (shallow) the ocean mixed layer by Ekman pumping, and the warm (cold) SST anomaly is strengthened by a warm (cold) geostrophic flow from the north (south). The whole anomaly system is advected with the Antarctic Circumpolar Current (ACC), and the SLP anomaly is phase-locked to the SSTanomaly. Reference Peterson and WhitePeterson and White (1998) suggested that the ACW has a major source in the western subtropical South Pacific, where interannual anomalies in SST develop in response to the El Niño-Southern Oscillation (ENSO) near the Equator (see Reference Peterson and WhitePeterson and White, 1998, fig. 1). Subsequently, these ENSO-scale SSTanomalies propagate southward into the Southern Ocean in tandem with SLP anomalies. The system then migrates east around the globe through a combination of oceanic advection with the ACC and ocean-atmosphere coupling. In contrast to Reference Peterson and WhitePeterson and White (1998), Reference Christoph, Barnett and RoecknerChristoph and others (1998) and Reference Cai, Baines and GordonCai and others (1999) found that ENSO is not a significant source for the ACW. They argued that the ACW has its origins in the southern mid-to high latitudes.
Modeling studies cast doubt on the coupled atmosphere-ocean explanation of the ACW. The breakdown of phase-locking in atmospheric and oceanic components of the ACW in the coupled model studies of Reference Christoph, Barnett and RoecknerChristoph and others (1998) and Reference Cai, Baines and GordonCai and others (1999) suggests that the feedback of SST anomalies on the atmospheric circulation is weak. Reference Christoph, Barnett and RoecknerChristoph and others (1998) concluded that the ACW represents the net result of propagating oceanic anomalies interacting with a spatially fixed resonant atmospheric background pattern such that the oceanic anomalies are selectively amplified or dissipated. Reference Weisse, Mikolajewicz, Sterl and DrijfhoutWeisse and others (1999) used an ocean-only model driven by stochastic atmospheric forcing with a spatially fixed zonal wavenumber-3 pattern, and showed that the irregular atmospheric fluctuations in time are sufficient to excite propagating signals in the ocean surface layer, generating similar ACW patterns described by Reference Christoph, Barnett and RoecknerChristoph and others (1998). The large-scale spatial coherence of the atmospheric forcing and the average zonal velocity of the ACC determine the time-scale of the oceanic variability.
In this study, we examine the ACW using a global coupled ocean-ice model forced with the ERA data. Section 2 contains a brief description of the sea-ice and ocean models and the atmospheric forcing data. Section 3 presents the propagation characteristics of the ACW. Section 4 discusses the temporal phase relation among variables in atmosphere, ocean and sea ice. Conclusions are given in section 5.
2. Model Description
The global coupled ocean-ice model consists of two pieces: (1) an almost-global, primitive-equation, z-coordinate ocean model (Reference Semtner and ChervinSemtner and Chervin, 1992) with added free-surface (Reference Killworth, Stainforth, Webb and PatersonKillworth and others, 1991) and Antarctic sea ice; and (2) an Arctic ice-ocean model. The almost-global ocean model utilizes a Mercator grid with an average horizontal grid size of 1/4° (∼20 km at 60° S or 60° N) and 20 vertical levels. The model bathymetry is obtained without applying any smoothing, which provides higher accuracy in mapping the true geometry onto the model domain. Due to numerical singularity at the Poles, the domain of the almost-global ocean model extends only from 75° S to 65 °N, thus omitting the Arctic Ocean. The Antarctic sea-ice model employs a polar stereographic projection with a horizontal grid size of 18 km. It uses an elastic-viscous-plastic ice rheology (Reference Hunke and DukowiczHunke and Dukowicz, 1997; Reference Hunke and ZhangHunke and Zhang, 1999) and the zero-layer approximation for heat conduction through ice (Reference SemtnerSemtner, 1976). The surface heat budget follows that of Reference Parkinson and WashingtonParkinson and Washington (1979).
To complete the global model domain, the almost-global ocean model is adapted to the Arctic Ocean (Reference ParsonsParsons, 1995; Reference MaslowskiMaslowski, 1997) by transforming the Mercator grid into a rotated coordinate system, in which the model Equator passes through the geographic North Pole. It has a 1/6° (-18 km) horizontal grid size and 30 vertical levels. The model domain includes the Arctic Ocean, Barents and Kara Seas, Nordic Seas, Baffin and Hudson Bays, Canadian Archipelago, and subpolar North Atlantic Ocean (for more information on the model, see Reference Zhang, Maslowski and SemtnerZhang and others (1999)). The Arctic sea-ice model uses the same grid as the Arctic ocean model and has the same dynamics and thermodynamics as the Antarctic sea-ice model described earlier. The almost-global model is joined to the Arctic model by means of interpolating between the different model grids both horizontally and vertically, and exchanging lateral boundary information in the subpolar North Atlantic.
Following initialization with Reference LevitusLevitus (1982) climatological temperature and salinity and uniform 2 m sea ice, the almost-global and the Arctic models were spun up separately for several decades. The 3−daily averaged ECMWF operational analysis was then used to force the two models for a decade. Finally, the two pieces of the ocean-ice model were run interactively for 15 years with daily-averaged ERA (1979−93) data (surface net heat flux, 10 m winds, precipitation and evaporation) as atmospheric forcing. The last 15 year model output is analyzed and presented here.
3. Eastward Propagation
In this section, time-longitude diagrams of the monthly anomalies of net heat flux (NHF) at the surface of sea ice or open ocean, MWS at 10 m, SIE, and ocean temperature and salinity are shown (Fig. 1) to illustrate the eastward propagation characteristic of the ACW. As in Reference White and PetersonWhite and Peterson (1996), the monthly anomalies in each variable at each grid-point are determined relative to the 15 year monthly mean values, thus removing the mean seasonal cycles. To further suppress the high- and low-frequency signals, the monthly anomalies are then filtered with a bandpass filter having a 3−7 year admittance window. All monthly anomalies are along 56° S except the SIE, which is for 5° longitude increments. No space filtering is employed. The sign convention used here is that positive values denote northward SIE anomalies, warm-ocean temperature anomalies, saltier than normal salinity, northward MWS anomalies and more than normal downward NHF.
We first examine the time-longitude diagrams of monthly anomalies of NHF and MWS from ERA (Fig. 1), which are used to drive the ocean-ice model. The NHF is the sum of net sensible- and latent-heat fluxes and net short- and longwave radiations at the surfaces of sea ice or open ocean. The MWS is computed using the bulk aerodynamic formula with ERA 10 m winds and stability-dependent drag coefficient (Reference Large and PondLarge and Pond, 1982). Consistent with the results from Reference Bonekamp, Sterl and KomenBonekamp and others (1999), eastward propagations of NHF and MWS anomalies are present from about 1983. Prior to 1983, standing oscillations exist in both NHF and MWS anomalies. For the propagating mode, it takes about 8−9 years for an individual phase to travel around the globe, and variability at a particular location occurs on a time-scale of 4−5 years. Thus, the eastward-propagating NHF and MWS anomalies have a spatial length scale characterized by wavenumber 2, consistent with the findings by Reference White and PetersonWhite and Peterson (1996). For the standing oscillation, a higher wavenumber (3−4) is present, resembling the atmospheric variability described by Reference Christoph, Barnett and RoecknerChristoph and others (1998).
Both the simulated and observed time-longitude diagrams of SIE anomalies are displayed in Figure 1. The observed SIE is obtained from the sea-ice concentration derived from the Scanning Multichannel Microwave Radiometer (SMMR) and Special Sensor Microwave/Imager (SSM/I)brightness temperatures using the NASA Team algorithm (Reference Gloersen, Cavalieri, Campbell, Comiso, Parkinson and ZwallyGloersen and others, 1992; Reference Comiso, Cavalieri, Parkinson and GloersenComiso and others, 1997). Similar to NHF and MWS anomalies, both standing and eastward-propagating oscillations are present in the SIE anomaly. The fact that the standing oscillation is also present in the satellite passive-microwave observations indicates that the sudden switching-on of the ERA forcing (see section 2) has little impact on the modeled SIE oscillation in the beginning of the time series. The propagating mode of the SIE anomaly has the same oscillation period and phase speed as the NHF and MWS anomalies, suggesting the presence of eastward propagation in the SIE anomaly is a response to that in the atmosphere. Notice that the simulated propagating and standing anomalies of the SIE are in very good agreement with the observed.
Like the SIE anomaly, the anomalies of the first-level (mid-depth: 12.5 m) ocean temperature and salinity (hereafter referred to as SST and SSS, respectively) also have a pronounced eastward propagation during the period 1983−93 (Fig. 1). There are indications of standing oscillation from 1979 to 1982, but the signal is not as clear as for the SIE. Their phase speeds and oscillation periods are the same as those of SIE, NHF and MWS. The predominant wavenumber-2 pattern and the dominant time-scale of 4−5 years in the SSTanomaly is illustrated further by its first complex empirical orthogonal function (CEOF) pattern (Fig. 2) and the corresponding temporal normalized amplitude (Fig. 3). The CEOF analysis (Reference PreisendorferPreisendorfer, 1988) was conducted at our request by J. L. Annis of the Scripps Institution of Oceanography on the SST anomaly filtered with a bandpass filter having a 3−6 year admittance window. A wave partitioning is then performed to separate the dominant modes of CEOF into standing and propagating waves. Figure 2 shows the phase cycle of the propagating-wave component of the first CEOF mode, which explains 59% of the total variance. One sees not only the predominant wavenumber-2 pattern but also the eastward propagation of the SSTanomaly in Figure 2. A warm anomaly is initially seen over the central South Pacific, with cold anomalies to the west and east of it. In the following phases, the warm anomaly drifts eastward along with the cold anomalies. At the end of the sequence, part of the warm anomaly passes through Drake Passage.
The Hovmoeller diagrams of first-level temperature and salinity anomalies (Fig. 1) are illustrative for the other ocean levels down to 1000 m except that more coherent and intense anomalies are present over the eastern Pacific and Atlantic sectors than over the Indian and the western Pacific sectors in the interior ocean. The anomalies of the eighth-level (mid-depth: 310 m) ocean temperature and salinity (hereafter referred to as T310 and S310, respectively) are shown in Figure 1 to demonstrate this. The same dominant oscillation period and phase speed are found at these levels. Interestingly, the speed of eastward propagation of temperature and salinity anomalies appears to be independent of depth. As we know, the zonal velocity decreases with depth. In our model, the mean zonal velocity at 56° S is 0.097 m s"1 in the upper 100 m, and 0.085 m s"1 at 310 m depth. This suggests that, below the surface, advection of the ACC does not play as significant a role in the eastward movement of the anomalies. This seems to support the notion proposed by Reference Bonekamp, Sterl and KomenBonekamp and others (1999) that the circumpolar character of the marginal stratification in the Southern Ocean allows convection along the trajectory of the surface anomalies, which are advected with the surface ACC. In this way the anomalies can be maintained in time against the turbulent dissipation in the top and intermediate layers. The stronger temperature and salinity anomalies over the eastern Pacific and Atlantic sectors in the interior ocean imply stronger anomalous oceanic convections there.
The eastward propagations of ocean temperature and salinity anomalies in the interior ocean are also found in the model studies of Reference Bonekamp, Sterl and KomenBonekamp and others (1999) and Reference Cai, Baines and GordonCai and others (1999). However, the predominant wavenumber in their models is 3 whereas in our model it is 2. One significant difference between our ocean model and theirs is the higher horizontal and vertical resolutions in our model, allowing us to reasonably simulate the local strength of the ACC which consists of many eddies and narrow jet streams. The coarse-resolution models tend to underestimate the strength of the ACC, although they obtained reasonable Drake Passage transport due to the too wide Drake Passage. Assuming the oceanic part of the ACW is strongly governed by advection as claimed by Reference Christoph, Barnett and RoecknerChristoph and others (1998) and Reference Weisse, Mikolajewicz, Sterl and DrijfhoutWeisse and others (1999), the strength of the ACC would largely determine the time-scale of the oceanic component of the ACW. This may account for the slower phase speed of the oceanic ACW in the model of Reference Bonekamp, Sterl and KomenBonekamp and others (1999) than in ours, since their simulated average ACC strength is about one-half of ours even though the two models are forced by the same ERA atmospheric forcing.
4. Phase Relationship
In this section, the phase relationship among the various variables discussed in section 3 is examined to understand the dynamics of the ACW. Note that the cold (warm) ocean temperature anomalies, northward (southward) SIE anomalies and saltier (fresher) than normal ocean salinity anomalies are in phase (Fig. 1). This phase relationship is confirmed by the almost zero-lag cross-correlations between SST and SIE anomalies, and between SSTand SSS anomalies (Fig. 4). The zero-lag cross-correlation between SST and T310 anomalies indicates that the temperature anomalies at the two depths are in phase. The lagged cross-correlations of NHF anomalies vs SSTand SIE anomalies show warm (cold) SSTand southward (northward) SIE anomalies trailing more (less) than normal downward NHF by about 10 months (or ∼ 90° out of phase). The lagged cross-correlation of MWS vs SIE indicates that northward (southward) MWS anomalies and northward (southward) SIE anomalies are in phase.
This phase relation supports the schematic view of the ACW proposed by Reference White and PetersonWhite and Peterson (1996, fig. 3) if one overlays less (more) than normal NHF on their high (low) SLP anomalies. (We did not analyze the ERA SLP anomaly because we did not use ERA SLP to force the model.) It can be summarized as follows: less (more) than normal NHF, coinciding with high (low) SLP anomalies, induces northward (southward) MWS anomalies 90° downstream as required by the geostrophic balance; the anomalous advection of cold (warm) air associated with the anomalous northward (southward) MWS cools (heats) the ocean surface, resulting in a cold (warm) SST anomaly; at the same time, the northward (southward) MWS anomaly advects more (less) sea ice northward, and the cooler (warmer) air and ocean allows more (less) sea-ice growth, inducing northward (southward) SIE anomalies; the surface anomalies are advected with the ACC, and the anomalous oceanic convection along the path transports the ACW surface signal downwards and sustains the anomalies against the turbulent dissipation in the top and intermediate layers, especially over the eastern South Pacific and South Atlantic Oceans.
5. Conclusions
This study demonstrates that a one-way atmosphere-to-ocean-ice coupling is sufficient to generate the ACW signature in ocean and sea ice, providing the ACW signature is present in the atmosphere. We cannot determine whether the presence of the ACW in the atmosphere depends on positive feedback from the ocean and sea ice since our model does not interact with the atmosphere. We also do not know whether the existence of the ACW in ERA depends on feedback from the ocean since the SST, used to drive the ECMWF model, is mainly derived from the observational dataset of Reference Reynolds and SmithReynolds and Smith (1994), which contains the ACW signature (Reference White and PetersonWhite and Peterson, 1996). However, a recent analysis by S. G. Warren and J.L. Annis (personal communication, 2000) of 100 years from a coupled climate simulation having the same atmospheric resolution as in Reference Christoph, Barnett and RoecknerChristoph and others (1998) but higher oceanic resolution does show the ACW with wavenumber-2 anomalies propagating in both the atmosphere and ocean. This suggests that oceanic feedback to the atmosphere is essential for the correct phase propagation in the atmospheric component of the ACW. Ongoing analysis of this simulation from the Parallel Climate Model of Reference WashingtonWashington and others (2000) is being pursued by Warren and Annis.
Acknowledgements
We offer sincere thanks to W. Maslowski and R.Tokmakian at the Naval Postgraduate School, Monterey, for their effort in joining the almost-global model to the Arctic model, and kudos especially to J.L. Annis at the Scripps Institution of Oceanography, University of California, San Diego, for performing the CEOF analysis. We gratefully acknowledge support by the DOE CHAMMP and Climate Change Prediction programs and the U.S. National Science Foundation ARCSS and OPP (grant OPP-9980607) programs. Model simulations were performed mainly on Climate Simulation Laboratory computers at the U.S. National Center for Atmospheric Research. The sea-ice concentration derived from the satellite passive-microwave observations was provided by the Earth Observing System Distributed Active Archive Centre at the U. S. National Snow and Ice Data Center.