Atlantic and Pacific multidecadal oscillations and Northern Hemisphere temperatures

See allHide authors and affiliations

Science  27 Feb 2015:
Vol. 347, Issue 6225, pp. 988-991
DOI: 10.1126/science.1257856

Is the end of the warming hiatus nigh?

Which recent climate changes have been forced by greenhouse gas emissions, and which have been natural fluctuations of the climate system? Steinman et al. combined observational data and a large collection of climate models to assess the Northern Hemisphere climate over the past 150 years (see the Perspective by Booth). At various points in time, the Pacific Decadal Oscillation and the Atlantic Multidecadal Oscillation have played particularly large roles in producing temperature trends. Their effects have combined to cause the apparent pause in warming at the beginning of the 21st century, known as the warming “hiatus.” This pause is projected to end in the near future as temperatures resume their upward climb.

Science, this issue p. 988; see also p. 952


The recent slowdown in global warming has brought into question the reliability of climate model projections of future temperature change and has led to a vigorous debate over whether this slowdown is the result of naturally occurring, internal variability or forcing external to Earth’s climate system. To address these issues, we applied a semi-empirical approach that combines climate observations and model simulations to estimate Atlantic- and Pacific-based internal multidecadal variability (termed “AMO” and “PMO,” respectively). Using this method, the AMO and PMO are found to explain a large proportion of internal variability in Northern Hemisphere mean temperatures. Competition between a modest positive peak in the AMO and a substantially negative-trending PMO are seen to produce a slowdown or “false pause” in warming of the past decade.

Distinguishing between forced and unforced variability in climate is critical for assessing the impact of anthropogenic forcing on temperature, drought, hurricane activity, weather extremes, and other climate phenomena. The North Atlantic and North Pacific oceans are the key drivers of internal variability in Northern Hemisphere temperatures on multidecadal time scales, but there is substantial uncertainty in their relative contributions to the observed variability. We applied a new semi-empirical method that uses a combination of observational data and a large ensemble of coupled climate model simulations to assess the relative roles of both forced and internal variability in the Northern Hemisphere over the historical period.

The Atlantic Multidecadal Oscillation (AMO) (1) is the leading mode of internal variability in North Atlantic sea surface temperature (SST) on multidecadal (~50 to 70 years) time scales (24). The Pacific Decadal Oscillation (PDO) (5, 6) is the leading mode of North Pacific internal SST variability but, as defined, consists of at least two distinct signals, one roughly bidecadal with a ~16- to 20-year period and the other multidecadal with a ~50- to 70-year period (4, 79). The PDO and AMO time series typically are defined in terms of the temporal pattern of temperature change in the north-central Pacific and North Atlantic, respectively. The multidecadal component of the PDO may in part be related to the AMO [although centered in the Atlantic, it appears (2, 3) to project at least weakly onto the Pacific] and in part reflective of low-frequency variability related to the El Niño–Southern Oscillation (ENSO) and its extratropical response (1016). We distinguish the multidecadal component from the conventionally defined PDO by terming it the “PMO,” and we term the multidecadal component of internal Northern Hemisphere mean temperature variability the “NMO.”

Prior methods used to define these internal variability components and their influence on Northern Hemisphere temperature include (i) a simple linear detrending of the mean North Atlantic SST time series (1721), (ii) estimating the forced trend based on regression of North Atlantic SST against global mean SST and removing the forced trend to yield an estimate of the internal variability (16, 22, 23), and (iii) defining the forced component as the mean of North Atlantic SST in an ensemble of climate model simulations and defining the internal variability component as the difference between the observed SST series and the multimodel mean (24, 25). These methods, as shown below, do not in general yield correct results. We estimated the Atlantic and Pacific-basin multidecadal internal variability components and their contribution to Northern Hemisphere temperature change on the basis of a new target region regression approach.

Our method is based on the principle that internal variability is uncorrelated among distinct realizations of a large ensemble. We therefore used the mean of the Coupled Model Intercomparison Project Phase 5 (CMIP5) ensemble (26) as an initial estimate of the forced component of surface temperature for the North Atlantic, North Pacific, and the entire (land + ocean) Northern Hemisphere region (Fig. 1). The estimated forced series is rescaled via linear regression against the actual temperature series so as to accommodate potential differences in the amplitude of the true forced response relative to the multimodel mean response (for example, because of disparities in climate sensitivity). We define the AMO, PMO, and NMO as the difference between the observations and estimated, regional forced temperature series for each of the three respective regions, low-pass filtered at a frequency of 40 years in order to isolate multidecadal variability (27).

Fig. 1 CMIP5-All ensemble means shown with individual model means.

(A) Northern Hemisphere SST+SAT. (B) North Atlantic SST. (C) North Pacific SST. Ensemble mean, black curves; individual model means, colored curves. Thin black line depicts the 95% confidence limits of the model mean determined via bootstrap resampling. Blue line depicts observed temperatures.

We analyzed both the subensemble of simulations (n = 24) of the GISS-E2-R model (28) (henceforth “CMIP5-GISS”); the subensemble of simulations (n = 45) of models (M = 15) with aerosol indirect effects (“CMIP5-AIE”); and the larger, full (n = 170 total realizations) ensemble of all (M = 44) models (“CMIP5-All”) (Fig. 1, fig. S1, and table S1). The three ensembles are complementary in their characteristics. The GISS-E2-R simulations (which comprise the largest CMIP5 ensemble for an individual model) are consistent in their forcings and include representation of the first aerosol indirect effect (cloud albedo). The CMIP5-AIE models all have full representations of both the first and second (cloud lifetime) indirect aerosol effects, which are potentially important contributions to the net radiative forcing (29). The CMIP5-All ensemble provides a much larger sample, but individual simulations vary in the forcings that were used and how they were implemented. Recent work (30) has explored the hypothesis that at least some of the difference between modeled and observed temperature changes arises from errors in the forcing estimates (for example, the accumulated effects of small volcanic eruptions over the past decade are not accounted for in the vast majority of CMIP5 simulations). Our assumption is that these three different ensembles mean estimates of the forced temperature signal span a representative range of uncertainty in the underlying forcing.

In defining the AMO, PMO, and NMO, we considered target regions spanning the equator to 60° north over the Atlantic (0° to 80°W) and Pacific (120°E to 100°W) oceans (the areal mean over all SST gridboxes in each basin), and the full Northern Hemisphere (ocean + land). The CMIP5-All multimodel ensemble mean series (latitude weighted) for each of the target regions, along with the ensemble of individual simulations, were compared with the actual historical observations over the interval 1854–2012 C.E. (Fig. 1 and fig. S1) (27). We used Goddard Institute for Space Studies (GISS) Surface Temperature (GISTEMP) (31) for the observational NH mean (ocean + land) series, owing to recent evidence (32) that other products may underestimate recent warming by under-sampling the Arctic. For the regional observational SST estimates, we used the mean of the Hadley Centre Global Sea Ice and Sea Surface Temperature (HadISST) (33), National Oceanic and Atmospheric Administration (NOAA) Extended Reconstructed Sea Surface Temperature (ERSST) (34, 35), and Kaplan (3638) products.

The results of the target-region regression analysis show for each of the three model ensembles that the estimated internal variability components derived from the various realizations are statistically independent, as they should be if the method is performing correctly, contrasting with what we find for the other previously used methods (Fig. 2, figs. S2 to S4, and table S2) (27). We next applied the methods in a semi-empirical setting in order to estimate the actual historical AMO, PMO, and NMO series. Under the assumption that the observational temperature series are the sum of a forced component and the real-word realization of internal variability, we estimate the true historical internal variability component as the residual series after the forced components are removed.

Fig. 2 CMIP5-All mean series and estimated 1-σ bounds for mean series under the assumption of statistical independence of internal variability among ensemble members.

(A) AMO. (B) PMO. (C) NMO. Solid lines indicate mean of N-1 realizations; dashed lines indicate estimated 1-σ bounds determined by using detrending (blue), global SST regression (red), and target region regression (black). Individual realizations of CMIP5-All internal variability as well as results for target region differencing are shown in the supplementary materials (fig. S2).

Our approach gives similar results whether CMIP5-All, CMIP5-GISS, CMIP5-AIE [or even individual models with a minimum of n ≥ 10 realizations (fig. S6)] ensemble means are used (39). The root mean square amplitude of the AMO and PMO are similar for all three ensembles (0.10/0.11/0.09°C for AMO and 0.09/0.09/0.11°C for PMO, for CMIP5-All/CMIP5-GISS/CMIP5-AIE, respectively). Unlike with the linear detrending approach, the PMO and AMO are not found to be significantly correlated. An analysis of the full multimodel ensemble reveals any putative correlation between the AMO and PMO [and arguments of a “stadium wave” climate signal (40)] to be an artifact of the linear detrending approach (fig. S7) (27). Shown also (Fig. 3) are the results of a simple bivariate regression demonstrating that the NMO can be very closely approximated [coefficient of determination (R2) = 0.86/0.88/0.91 for CMIP5-All/CMIP5-GISS/CMIP5-AIE, respectively] by a weighted combination of the AMO and PMO series (41). The amplitude of the NMO (0.07°C using either CMIP5-All or CMIP5-GISS, and 0.08°C using CMIP5-AIE) is consistent with results from long model control runs (3).

Fig. 3 Semi-empirical estimate of AMO, PMO, and NMO based on target region regression using historical climate model realizations.

(A) CMIP5-GISS. (B) CMIP5-AIE. (C) CMIP-All. In (A) to (C), blue, AMO; green, PMO; and black, NMO. Bivariate regression-based approximation of NMO (red) strongly correlates (R2 = 0.86/0.88/0.91 for CMIP5-All/CMIP5-GISS, CMIP5-AIE, respectively) with semi-empirical NMO estimate (black). 95% confidence limits of the AMO, PMO, and NMO CMIP5-All means were determined by using the ensemble of target region mean series resulting from bootstrap resampling (Fig. 1) and are shown as colored shading.

Our analysis shows the NMO to be decreasing at the end of the series (Fig. 3 and figs. S5 and S6). Mann et al. (42) assessed the recent decrease in the NMO in terms of a negative-trending AMO contribution. However, we reach a somewhat different conclusion in the present study, finding that the recent decrease in the NMO is instead a result of a sharply decreasing PMO (with a relatively flat AMO contribution). That observation is consistent with recent findings that the anomalous slowing of warming over the past decade is tied to subsurface heat burial in the tropical Pacific and a tendency for persistent “La Niña”–like conditions (4346). Our analysis attributes this trend to internal variability as a consequence of the failure of the CMIP5 models to identify a recent forced trend of this nature. However, there is paleoclimate evidence suggesting that a La Niña–like response might arise from positive radiative forcing (47), and the possibility remains that state-of-the-art climate models fail to capture such a dynamical response to anthropogenic radiative forcing.

Some recent work (18, 19, 21, 22, 25) has attributed a potentially large proportion of observed regional and hemispheric temperature changes to multidecadal internal variability related to the so-called “AMO” and/or “PDO.” Using the CMIP5 multimodel historical climate simulations, we have established that the methods used in these studies tend to inflate and distort the estimated internal variability owing to an incorrect partition of internal and forced variability. We have demonstrated that our target-region regression method correctly isolates the internal variability components.

Applying our method to observational surface temperature data, we find that internal variability is likely to have had a substantial influence on multidecadal Northern Hemisphere temperature changes over the historical period, contributing up to 0.15°C peak warming/cooling. The AMO appears to have been influential in the early and middle 20th century, but the PMO has played a more dominant role in recent decades. This result is consistent across the three ensembles (GISS, AIE, and All) (Fig. 3). Our findings (the AIE experiments, especially) suggest that natural internal variability has had a modest influence on Atlantic SST over the past half century and that multidecadal climate variability attributed to Atlantic SST changes (such as variations in tropical storm frequency and strength and Sahel and Midwestern North American drought) (4851) was largely driven by external forcing (as concluded in other recent work) (52). Our results also highlight the substantial uncertainties associated with the role of anthropogenic aerosol forcing in recent decades because the greatest discrepancies using the three different ensembles occur during that time period.

Our findings have strong implications for the attribution of recent climate changes. We find that internal multidecadal variability in Northern Hemisphere temperatures (the NMO), rather than having contributed to recent warming, likely offset anthropogenic warming over the past decade. This natural cooling trend appears to reflect a combination of a relatively flat, modestly positive AMO and a sharply negative-trending PMO. Given the pattern of past historical variation, this trend will likely reverse with internal variability instead, adding to anthropogenic warming in the coming decades.

Supplementary Materials

Materials and Methods

Supplementary Text

Figs. S1 to S7

Tables S1 to S3

References (5358)

References and Notes

  1. Materials and methods are available as supplementary materials on Science Online.
  2. The regression analyses applied to the full CMIP5 multimodel mean yields a scaling factor (“beta”) for Northern Hemisphere temperature changes that slightly exceeds unity (beta = 1.053 ± 0.0169), implying a real-world forced response that is slightly greater than that estimated by the CMIP5 multimodel mean. In contrast, North Atlantic mean temperatures yields a scaling factor slightly below unity (beta = 0.916 ± 0.0155), and North Pacific mean temperatures yield a scaling factor substantially below unity (beta = 0.629 ± 0.0182), suggesting that the CMIP5 multimodel mean substantially overestimates the amplitude of forced temperature changes over the North Pacific. Further details, including results for the two subensembles (CMIP5-A1E and CMIP5-GISS), are available in the supplementary materials (table S3).
  3. NMO = 0.35 AMO + 0.43 PMO for CMIP5-All; NMO = 0.42 AMO + 0.36 PMO for CMIP5-GISS; NMO = –0.06 AMO + 0.85 PMO for CMIP5-AIE; AMO and PMO regression coefficients are significant at the P << 0.05 level based on a one-sided test.
  4. Acknowledgments: All raw data, Matlab code, and results from our analysis are available at the supplementary website: We acknowledge the World Climate Research Programme’s Working Group on Coupled Modeling, which is responsible for CMIP, and we thank the climate modeling groups for producing and making available their model output. We thank K. Emanuel and G. Schmidt for helpful comments on earlier versions of the manuscript. B.A.S. acknowledges support by the U.S. National Science Foundation Atmospheric and Geospace Sciences–Postdoctoral Research Fellowships (AGS-PRF) (AGS-1137750). Kaplan SST V2 data were provided by the NOAA/Office of Oceanic and Atmospheric Research/Earth System Research Laboratory Physical Sciences Division, Boulder, Colorado, USA: HadISST data were provided by the Met Office Hadley Centre: ERSST data were provided by NOAA:
View Abstract

Navigate This Article