The Pacific decadal oscillation (PDO) is the leading mode of sea surface temperature (SST) variability over the North Pacific (north of 20°N). Its South Pacific counterpart (south of 20°S) is the South Pacific decadal oscillation (SPDO). The effects of tropical eastern Pacific (TEP) SST forcing and internal atmospheric variability are investigated for both the PDO and SPDO using a 10-member ensemble tropical Pacific pacemaker experiment. Each member is forced by the historical radiative forcing and observed SST anomalies in the TEP region. Outside the TEP region, the ocean and atmosphere are fully coupled and freely evolve. The TEP-forced PDO (54% variance) and SPDO (46% variance) are correlated in time and exhibit a symmetric structure about the equator, driven by the Pacific–North American (PNA) and Pacific–South American teleconnections, respectively. The internal PDO resembles the TEP-forced component but is related to internal Aleutian low (AL) variability associated with the Northern Hemisphere annular mode and PNA pattern. The internal variability is locally enhanced by barotropic energy conversion in the westerly jet exit region around the Aleutians. By contrast, barotropic energy conversion is weak associated with the internal SPDO, resulting in weak geographical preference of sea level pressure variability. Therefore, the internal SPDO differs from the TEP-forced component, featuring SST anomalies along ~60°S in association with the Southern Hemisphere annular mode. The limitations on isolating the internal component from observations are discussed. Specifically, internal PDO variability appears to contribute significantly to the North Pacific regime shift in the 1940s.
The Pacific decadal oscillation (PDO) is the leading mode of monthly sea surface temperature (SST) anomalies in the North Pacific (poleward of 20°N) (Mantua et al. 1997; Zhang et al. 1997). It exhibits a horseshoe-like SST pattern, with anomalies of one sign in the Kuroshio–Oyashio extension (KOE) region (Nakamura et al. 1997; Nonaka et al. 2006) and of the opposite sign along the west coast of North America. The PDO predominantly varies on decadal-to-multidecadal time scales (Minobe 1997, 1999). It impacts U.S. summertime hydrologic variability (Barlow et al. 2001), North Pacific marine ecosystems, fisheries, and global climate (Miller and Schneider 2000; Mantua and Hare 2002).
Numerous studies have investigated the PDO mechanisms. A recent consensus is that the PDO is not a single mode but a combination of multiple physical processes (Schneider and Cornuelle 2005; Newman et al. 2016). Newman et al. (2016) decomposed the PDO into two components: tropical Pacific remote forcing and local North Pacific variability, both of which are associated with atmospheric Aleutian low (AL) variability. The warm phase of El Niño–Southern Oscillation (ENSO) excites the Pacific–North American (PNA) pattern, extending across the North Pacific to North America with a center of action at the AL (Horel and Wallace 1981). The strengthened AL over the North Pacific drives the positive PDO pattern through surface heat fluxes, wind-driven mixing, and Ekman advection (Alexander et al. 2002). AL variability also originates from atmospheric internal dynamics, forcing a reddened ocean response (Frankignoul and Hasselmann 1977). SST in the KOE region can vary through dynamical adjustment of upper-ocean gyre circulations excited by anomalous wind stress curl. The SST changes in turn feed back to the atmosphere, forming a coupled ocean–atmosphere process in the North Pacific (Zhang and Delworth 2015; Newman et al. 2016). Here we use a 10-member ensemble of partially coupled simulations to separate tropical forced and internal components of the PDO. Each member run shares the common radiative forcing and prescribed tropical Pacific SST variability and differs only in the initial condition. Here internal variability refers to the deviations from the ensemble mean and is unrelated to tropical Pacific variability. This internal component includes not only the atmospheric internal variability but also the coupled air–sea interaction process discussed above.
The South Pacific counterpart of the PDO (poleward of 20°S) is termed the South Pacific decadal oscillation (SPDO; e.g., Chen and Wallace 2015). The observed SPDO shows an equatorially symmetric SST structure with the PDO. Tropical Pacific forcing is suggested to be important for the SPDO (Shakun and Shaman 2009). The tropically excited Pacific–South American (PSA) pattern (Karoly 1989) is important in forcing the SPDO, mainly through the surface heat fluxes (Okumura 2013). It remains unclear if there is an SPDO component independent of tropical Pacific variability.
Previous studies have attempted to separate the PDO into tropical Pacific-forced and local internal variability components by means of statistical methods (e.g., Zhang et al. 1996; An and Wang 2005) and a hierarchy of coupled models. Ocean mixed layer models were used to study the ENSO-forced North Pacific variability, but they consider only thermodynamic coupling processes in the North Pacific (Alexander et al. 2002; Di Lorenzo et al. 2015). Partial coupling experiments were conducted to study the intrinsic North Pacific variability by prescribing the climatological SST in the tropics but with the ocean and atmosphere fully coupled in the extratropics (Zhang and Delworth 2015). Some studies investigated the impact of atmospheric noise on the North Pacific SST variability using an interactive ensemble coupling strategy (Kirtman and Shukla 2002; Yeh et al. 2007). Wang et al. (2012) compared the partial coupling experiment with a fully coupled run to investigate the influence of ENSO on the PDO. All the above studies have suggested that the PDO is made up of the tropical Pacific forced and intrinsic North Pacific components, both related to AL variability. However, not much attention has been given to possible differences in the atmospheric teleconnection pattern between the tropical forced and internal PDO. Even fewer studies examine the tropical-forced and internal components for the SPDO. It is unclear whether the SPDO components are similar in spatial structure to each other.
The present study uses a 10-member ensemble of the tropical Pacific pacemaker experiment (Kosaka and Xie 2016) to investigate the tropical eastern Pacific (TEP)-forced and internal components of the PDO and SPDO, with a focus on the associated atmospheric teleconnection patterns. We identify distinct features between the TEP-forced and internal components of the SPDO. We suggest that the small degrees of freedom limit our ability to isolate the internal component from observations. The comparison with the pacemaker experiment indicates that the PDO regime shift in the 1940s (Minobe 1997, 1999) is related to the internal variability.
The rest of the paper is organized as follows. Section 2 describes the pacemaker experiment, the method of isolating the TEP-forced and internal components, and the observational dataset used in this study. Section 3 compares the simulated PDO and SPDO variabilities by the pacemaker experiment with observations. Section 4 investigates the distinct characteristics between the TEP-forced and internal components for both the PDO and SPDO. Furthermore, barotropic energy conversion is diagnosed for atmospheric anomalies to explain the similarity (distinction) of the two components for the PDO (SPDO). Section 5 discusses the limitations on isolating the internal component from observations and identifies the relative contribution of TEP-forced and internal components to the PDO regime shifts. Section 6 is a summary with some discussion.
2. Data and methods
A long-term tropical Pacific pacemaker experiment, named Pacific Ocean–Global Atmosphere (POGA) (Kosaka and Xie 2016), is used to isolate the TEP-forced and internal components based on the Geophysical Fluid Dynamics Laboratory Coupled Model, version 2.1 (CM2.1; Delworth et al. 2006). The POGA ensemble has 10-member runs, each of which is forced by the identical historical radiative forcing and observed SST anomalies in the TEP domain (from 180° to the American coast, 15°S–15°N, with a 5° buffer zone north, south, and west of the domain), but with different atmospheric and oceanic initial conditions. Outside the TEP-restoring region, the ocean and atmosphere are fully coupled and freely evolve. The horizontal resolution of the ocean is approximately 1° longitude × 1° latitude, with finer meridional resolution equatorward of 30°S and 30°N (⅓° at the equator). The atmospheric resolution is about 2.5° longitude × 2° latitude. The time period of 1890–2014 is analyzed. We refer details on the model framework to Kosaka and Xie (2013, 2016).
To obtain the TEP-forced component, we calculate the ensemble mean of 10 POGA members and then subtract the ensemble mean of 20-member historical runs, each of which is forced by the historical and representative concentration pathway 4.5 (RCP4.5) radiative forcing and based on CM2.1, to remove the radiatively forced response. The internal component is obtained by subtracting the POGA ensemble mean from each POGA member. Note that the internal component here refers to variability independent of the prescribed TEP SST. Although 10 members may be not sufficient to extract the TEP-forced and internal components in a single event, the PDO and SPDO are obtained from empirical orthogonal function (EOF) analysis, which also exploits the degrees of freedom in the time dimension. Our analyses suggest that 10 members are sufficient for statistically isolating the TEP-forced and internal variability of the PDO and SPDO [see the strong (weak) connection of SST anomalies with the tropical Pacific in Figs. 3a and 5a (Figs. 3d and 5d)]. In addition to CM2.1 POGA, we also examine the POGA results based on the Community Earth System Model, version 1 (CESM POGA; Deser et al. 2017; see the appendix).
We use the Extended Reconstructed Sea Surface Temperature product, version 5 (ERSSTv5; Huang et al. 2017). Other SST datasets show similar results. The analysis period begins from 1900 (until 2014) because of large uncertainties in observations prior to the early twentieth century. We also use the Southern Oscillation index (SOI) from 1890 to 2014 to represent long-term tropical Pacific variability. This index is based on annual standardization, which has shown to maximize the signal of large-scale variability (Trenberth 1984).
We focus on the PDO and SPDO variabilities on interannual-to-multidecadal time scales. In the internal component, we subtract climatology in each ocean variable of each member to remove initial perturbations that persist throughout the integration in the Southern Ocean. The annual mean centered on the local wintertime is calculated, from July to the following June in the Northern Hemisphere (NH), and January–December in the Southern Hemisphere (SH). The reemergence process occurs in local winter and is important for the low-frequency PDO and SPDO (Newman et al. 2003; Shakun and Shaman 2009).
To isolate the TEP-forced covariability mode, we conduct joint EOF analysis of annual-mean (July–June) SST anomalies between the North and South Pacific. The analyses based on the annual mean of January–December are similar (not shown). SST anomalies are normalized by the spatially averaged SST standard deviation in each basin. This normalization method assigns equal weight on SST variability between the North and South Pacific.
3. Evaluating POGA simulation
We first evaluate the POGA’s skill in simulating the spatial–temporal structures of the observed PDO and SPDO. Figure 1 shows the spatial patterns of the PDO and SPDO along with the corresponding SLP anomalies. In general, POGA well simulates the PDO, with a horseshoe-like SST pattern associated with strengthened AL variability, although some biases exist, such as a slight westward shift in the main center of the AL anomaly, stronger SST anomalies in the KOE region, and weaker SST anomalies in the subtropical northeastern Pacific (Fig. 1a vs Fig. 1b). POGA also well simulates the SPDO pattern, characterized by SST cooling below the South Pacific convergence zone (SPCZ) region and around New Zealand and SST warming in the southeastern South Pacific and Southern Ocean. It is associated with low pressure anomalies in the midlatitudes, with maxima south of New Zealand, and high pressure anomalies centered around 65°S, 120°W (Fig. 1c vs Fig. 1d).
We further compare the temporal characteristics of the PDO and SPDO between POGA and observations (Fig. 2). Overall, POGA well simulates the PDO (Fig. 2a) and SPDO spectra (Fig. 2b), although the latter of which shows strong interannual and weak decadal-to-multidecadal variability. The model also captures the correlation between ENSO and the PDO or SPDO, albeit biased somewhat in lead/lag (Figs. 2c,d). The correlation between the PDO and SPDO in POGA is lower than in observations (Fig. 2e). This may arise from too strong internal variability in POGA (see further discussion in section 5a).
Nevertheless, POGA well simulates both the spatial and temporal characteristics of the observed PDO and SPDO. Thus, we can further investigate the TEP-forced and internal components based on POGA.
4. TEP-forced versus internal variability
This section explores the characteristics of the TEP-forced and internal components separately for the PDO and SPDO. We identify distinct spatial patterns between these two components, especially for the SPDO.
Figure 3 shows the TEP-forced and internal components of the PDO. Generally, both the TEP-forced and internal PDOs exhibit a horseshoe-like SST pattern in the North Pacific, although the internal component features strong SST cooling over the KOE region (Fig. 3a vs Fig. 3d). The TEP-forced PDO is closely related to the ENSO-like SST pattern in the tropics, while the internal PDO is confined to the extratropical North Pacific. Both the TEP-forced and internal PDOs are associated with the anomalous AL over the North Pacific (Figs. 3b,e), consistent with previous findings (e.g., Newman et al. 2016). However, these two components are related to different NH atmospheric teleconnections. The TEP-forced PDO is driven by the PNA pattern (Wallace and Gutzler 1981) excited from the tropical Pacific (Fig. 3c), which differs from the tropical/Northern Hemisphere pattern in Barnston and Livezey (1987). From the SLP field, the internal PDO is not only associated with the anomalous AL, but also with an anomalously high SLP system over the Arctic polar cap region (Fig. 3e). This SLP pattern is reminiscent of the negative phase of the Northern Hemisphere annular mode (NAM), with the center of action over the entire Arctic Ocean surrounded by a zonal ring of the opposite sign centered near 45°N (Thompson and Wallace 1998, 2000). From the 200-hPa geopotential height (Z200) anomaly pattern, the internal PDO appears to be a hybrid of the NAM and PNA pattern (Fig. 3f).
To examine whether the atmospheric teleconnection patterns associated with the internal PDO are the dominant modes of internal atmospheric circulation in the NH, we obtain EOFs of the internal component of the annual-mean (July–June) SLP anomalies poleward of 20°N. Figures 4a and 4b show the regressions of SLP anomalies against the normalized principal components (PCs). The first EOF (SLP1) is the NAM (Fig. 4a), and the second EOF (SLP2) shows an anomalously low pressure over the Aleutian Islands (Fig. 4b). The SLP2 corresponds to the PNA pattern, which is the second EOF mode of the NH atmospheric circulation anomaly in observations, as reported by Wallace and Thompson (2002). Spatial correlations of SLP1 and SLP2 with the SLP pattern related to the internal PDO (Fig. 3e) demonstrate that the internal PDO is primarily associated with the internal NAM teleconnection, and secondarily with the internal PNA teleconnection (Table 1).
Figure 5 shows the TEP-forced and internal components of the SPDO. The TEP-forced SPDO connects to the ENSO-like variability and exhibits a rather symmetric structure about the equator with the TEP-forced PDO (Figs. 3a, 5a). It is characterized by an SST cooling below the SPCZ region surrounded by SST warming on both the equatorward and poleward sides (Fig. 5a). Analogous to the TEP-forced PDO, the TEP-forced SPDO is driven by the PSA pattern (Mo and Ghil 1987) excited from the tropical Pacific (Figs. 5b,c).
The internal SPDO is very different from the TEP-forced counterpart (Fig. 5a vs Fig. 5d). Notably, the internal SPDO exhibits strong SST warming along ~60°S with weak cooling around New Zealand, but without the cooling in the SPCZ region (Fig. 5d). The internal SPDO is not associated with the PSA pattern but with the Southern Hemisphere annular mode (SAM), which features a primary center of action over Antarctica that is regionally enhanced over the Amundsen Sea and western Antarctica, accompanied by a zonal ring of opposite sign along ~45°S (Figs. 5e,f) (Thompson and Wallace 2000). The SAM is the leading mode of internal SLP variability in the SH winter (Fig. 4c). High spatial correlation with the internal SAM (Fig. 5e) indicates a strong association with the internal SPDO (Table 1).
We evaluate the TEP-forced and internal components by projecting POGA anomalies onto the PDO (Fig. 1a) and SPDO (Fig. 1c) patterns. The variance of the PDO and SPDO explained by TEP forcing is 54% and 46%, respectively. This result suggests that both the TEP-forced and internal components are important for the PDO and SPDO variability.
c. Atmospheric forcing of the PDO and SPDO
This subsection examines the mechanism for SST variability in the PDO and SPDO. As both surface turbulent heat flux and ocean Ekman advection are related to surface wind variability, we perform 1-yr-lead regressions of the sum of these two terms against PC1s of the PDO and SPDO (contours in Fig. 6). Ekman advection is given by , where is the specific heat of seawater, is the Coriolis parameter, is the unit normal vector, and is the surface wind stress that is converted from the surface wind using quadratic stress law. Our calculation suggests that the Ekman advection is weak compared to heat flux. Thus, Fig. 6 can be mainly regarded as the surface heat flux effect.
For both the TEP-forced and internal PDOs (Figs. 6a,b), the overall SST pattern is consistent with surface heat flux forcing and Ekman advection. In the KOE, surface heat flux and Ekman advection damp SST anomalies, indicating that ocean dynamics plays a primary role. In the central-to-eastern North Pacific on the boundary between the SST cooling and warming regions, vertical entrainment due to the intensified wind speed (not shown) might also contribute to the SST cooling. For the TEP-forced and internal SPDOs (Figs. 6c,d), SST anomalies are basically driven by the PSA pattern and SAM, respectively. Thus, SST anomalies in the open ocean are driven by atmospheric wind variability through surface processes (heat flux and Ekman advection). In section 4d, we address the question of why the atmospheric circulation patterns associated with the TEP-forced and internal SPDOs are distinct but similar for the PDOs.
d. Barotropic energy conversion
For the PDO, both components are related to similar AL anomaly patterns in SLP (Fig. 7a). The SLP anomalies in the midlatitudes associated with the TEP-forced SPDO are geographically anchored south of Australia (negative value of shading in Fig. 7b). However, there is no geographical preference for SLP pattern in the midlatitudes associated with the internal SPDO (negative value of contour in Fig. 7b).
Over the North Pacific, previous studies suggested that the stationary AL variability is due partly to barotropic energy conversion (Simmons et al. 1983; Trenberth et al. 1998). It refers to a process through which disturbances can gain kinetic energy from the climatological zonally varying flow. This preferentially occurs in the westerly jet exit region as west of the Aleutians. Here we calculate the local barotropic energy conversion (CK) associated with the internal PDO and SPDO in Cartesian coordinates, following the formulation (cf. Simmons et al. 1983):
Here is the zonal velocity of the basic state; and and are the zonal and meridional velocity anomalies, respectively, which are obtained by regression against normalized PC1 of the internal PDO or SPDO. Positive CK indicates the conversion of kinetic energy from the mean flow to atmospheric circulation perturbations. The calculation is based on the wintertime-centered annual mean at the 200-hPa level.
Over the North Pacific, positive CK is large in the westerly jet exit region (Fig. 8a), indicating that strong barotropic energy conversion anchors AL variability, which in turn induces the internal PDO. Because the AL anomaly pattern is similar, the TEP-forced and internal components of the PDO both feature a horseshoe-like SST pattern (Figs. 3a,d). By contrast, CK over the South Pacific is weak despite a comparable magnitude of the circulation anomaly (Fig. 8b), resulting in a weak zonal preference for the atmospheric circulation anomaly associated with the internal SPDO. Consequently, the internal SPDO features pronounced SST warming along ~60°S, associated with the negative phase of the atmospheric internal SAM (Figs. 5d–f). The TEP-forced component shows strong zonal variations, with low SLP anomalies of the PSA pattern (shading in Fig. 7b) driving the SST cooling in the SPCZ region and SST warming in the northeast–southwest direction on the poleward side (Fig. 5a).
Ultimately, the westerly jets vary more strongly in the Northern than in the Southern Hemisphere. The resultant interhemispheric difference in barotropic energy conversion explains why the internal PDO resembles the TEP-forced component while the two components of the SPDO are markedly distinct. The results based on the CESM POGA are similar (see the appendix).
5. Discussion: Observed internal variability
a. Limitations on isolating the internal component
In the POGA ensemble mean, the temporal correlation is moderate between the PDO and SPDO (Table 2). This could be because they are sensitive to details in the structure and seasonality of TEP SST or because of remaining internal variability in the 10-member average. Alternatively, we can perform joint SST EOF between the North and South Pacific (see the method in section 2). The leading PC captures interhemispheric coherent variability and is highly correlated with both the TEP-forced PDO and SPDO (Table 2). This leading EOF resembles the so-called interdecadal Pacific oscillation (IPO; Power et al. 1999; Henley et al. 2015).
We apply the joint EOF method over the extratropical North and South Pacific basins to observations. Figure 9 shows the regressions of annual-mean (July–June) SST anomalies against the PCs of joint EOF, along with the POGA ensemble mean results. The results based on the January–December mean are similar (not shown). The joint EOF1 is highly correlated spatially and temporally between observations and the POGA ensemble mean (Figs. 9a and 9d; Table 3). This mode explains 38% and 24% of the variance in the POGA ensemble mean and observations, respectively. The second and higher modes disagree in spatial pattern between the POGA ensemble mean and observations (Figs. 9b,c vs Figs. 9e,f), indicating a limitation on reliably isolating higher modes by removing the TEP-forced variability over the North and South Pacific due to the small degrees of freedom in the instrumental record.
By design, POGA does not address the origin of TEP variability. The extratropical anomalies propagating toward the tropical Pacific are eventually restored by the observed TEP SST anomalies. The internal component might induce the TEP variability via the tropical–extratropical interactions. For example, the extratropical Pacific meridional mode can propagate into the tropics via the seasonal footprinting mechanism (Vimont et al. 2001, 2003; Di Lorenzo et al. 2015; Ma et al. 2017). The South Pacific is suggested to be coupled with the tropical Pacific on decadal time scales via the oceanic tunnel (Luo and Yamagata 2001). Okumura (2013) pointed out that atmospheric internal variability over the subtropical South Pacific can thermodynamically force the tropical Pacific.
b. PDO regime shifts
The PDO features marked regime shifts in the 1940s, the late 1970s, and the late 1990s (Minobe 1997, 1999). The joint EOF1 of the POGA ensemble mean suggests that the covariability between the PDO and SPDO is dominated by TEP forcing (Fig. 9a). Indeed, temporal correlation between the internal PDO and SPDO is nearly zero (r = −0.01). In observations, the PDO and SPDO exhibit strong covariability in the late 1970s and 1990s, and the POGA ensemble mean captures these regime shifts (Figs. 10a, A3a). This suggests that the 1970s and 1990s regime shifts are dominated by TEP forcing, consistent with previous analyses (e.g., Nitta and Yamada 1989; Kosaka and Xie 2016). By contrast, the PDO–SPDO covariability is weak during the 1940s regime shift (black line in Fig. 10a vs Fig. 10b). The PDO exhibits a regime shift but the POGA ensemble mean fails to capture it (Figs. 10a, A3a). As a caveat, the TEP SST itself is somewhat uncertain during World War II because of sparse observations (Thompson et al. 2008) as discussed in Kosaka and Xie (2016). The SOI, which reflects large-scale tropical Pacific variability, exhibits a consistent phase shift with the PDO, but over a much shorter period (Fig. 10a). This suggests that the internal variability plays an important role in the 1940s regime shift of the PDO.
We have used a 10-member ensemble of a long-term tropical Pacific pacemaker experiment and showed that the TEP-forced and internal components are distinct for both the PDO and SPDO. The TEP-forced PDO and SPDO exhibit an interhemispherically symmetric structure about the equator, forced by the PNA and PSA patterns excited from the tropical Pacific, respectively. The internal PDO resembles the TEP-forced component but is related to internal AL variability associated with the NAM and PNA pattern. The internal variability is locally enhanced by barotropic energy conversion in the westerly jet exit region around the Aleutians, driving a horseshoe-like SST pattern similar to the TEP-forced PDO. By contrast, barotropic energy conversion is rather weak and ineffective to form the preferred spatial pattern associated with the internal SPDO, resulting in weak geographical preference of SLP variability. The resultant SST anomaly pattern differs from the TEP-forced SPDO, featuring pronounced SST variability along ~60°S in association with the SAM teleconnection. These atmospheric patterns associated with the internal PDO and SPDO are the dominant modes of the internal atmospheric variability in each hemisphere. The more efficient barotropic energy conversion over the North than South Pacific explains why the internal PDO resembles the TEP-forced component while the internal and TEP-forced SPDOs are markedly distinct. The internal variability in the PDO appears to contribute significantly to the North Pacific climate regime shift in the 1940s, although short observations limit the confidence in the attribution.
As discussed in the introduction, coupled ocean–atmosphere interaction could be important for the PDO. The atmospheric resolution of CM2.1 is insufficient to adequately represent the crucial mesoscale ocean–atmosphere interaction in the KOE region (Ma et al. 2016). Future research needs to examine detailed physical processes of the PDO and SPDO—both the TEP-forced and internal components—with higher-resolution POGA experiments.
The ERSST, version 5, product is provided by the NOAA/Earth System Research Laboratory (www.esrl.noaa.gov/psd/data/gridded/data.noaa.ersst.v5.html), and the SOI time series is provided by Climate and Global Dynamics Division, NCAR (www.cgd.ucar.edu/cas/catalog/climind/soiAnnual.html). CESM1 historical, RCP8.5, and POGA runs have been obtained from the Earth System Grid (www.earthsystemgrid.org). We sincerely thank the three anonymous reviewers for their constructive comments and suggestions. This work is supported by the Natural Science Foundation of China (41490640 and 41490641), China Scholarship Council (201706330016), U.S. National Science Foundation (AGS 1637450), Japanese Ministry of Education, Culture, Sports, Science and Technology (the Integrated Research Program for Advancing Climate Models and Arctic Challenge for Sustainability), the Japan Society for the Promotion of Science (15H05466 and 18H01278), and the Japan Science and Technology Agency (Belmont Forum CRA “InterDec”).
The TEP-Forced and Internal Components of the PDO and SPDO in CESM POGA
We also use 10-member runs of CESM POGA to investigate the characteristics of the TEP-forced and internal components of the PDO and SPDO. The prescribed SST dataset is identical to that in CM2.1 POGA. The horizontal resolutions of both atmosphere and ocean are approximately 1° longitude × 1° latitude. The time period of the experiment is from 1920 to 2013. We use CESM 40-member historical and RCP8.5 runs to remove the response to the radiative forcing.
Both the TEP-forced and internal PDOs exhibit similar horseshoe-like SST patterns and are associated with AL variability (Fig. A1), consistent with the results based on the CM2.1 POGA. The TEP-forced SPDO is also characterized by an interhemispherically symmetric SST pattern with the TEP-forced PDO and is driven by the PSA pattern (Figs. A2a,b). As in CM2.1, the internal SPDO in CESM1 markedly differs from the TEP-forced SPDO, without SST cooling in the SPCZ region (Fig. A2c) and is also associated with the zonally symmetric SAM teleconnection (Fig. A2d). Likewise, the internal SPDO is notably distinct from the TEP-forced SPDO, associated with little geographical preference of SLP anomalies because of the weak barotropic energy conversion. Regarding the observed PDO regime shifts, the CESM POGA ensemble mean captures the late 1970s and 1990s regime shifts but not the one in the 1940s (Fig. A3a). This result is consistent with CM2.1 POGA, further suggesting that the internal variability in the PDO partly contributes to the 1940s regime shift.