Observation of Phase-Locked States in the Atlantic Multi-Decadal Oscillation (AMO) ()
1. Introduction
The AMO index is the area weighted average sea surface temperature (SST) of the Atlantic Ocean from latitudes 0N to 70N. The multi-decadal oscillation period is reported as 65 - 70 years by Schlesinger et al. [1] .
Studies of the AMO on decadal time scales are few. Douglass [2] reported that the AMO and three other climate indices all had the same climate shift dates and that these dates were the same as have been reported in other climate indices. Muller, Curry et al. [3] (henceforth MC2013) in a study of land surface temperature data also compared various SST data sets to the AMO with a focus on the decadal range, which they characterize as 2 to 15 years.
In this paper we will point out that we already know that the AMO has abrupt climate shifts [2] at certain dates and that inclusion of other indices shows that phase-locked states are observed between these dates. Thus, one would expect to find phase-locked states of AMO between these dates. This paper will show that this expectation is borne out. For example, we find that between the AMO shift dates 1986-87 and 2001-02 there is a phase-locked state of period 3 years.
2. Background Material
We summarize results of prior investigations that will be discussed in this paper.
2.1. Climate Shifts
Abrupt shifts in various climate indices have been well documented. Trenberth [4] was among the first to report existence of an abrupt shift in a climate index in the Pacific Ocean during the mid-1970s.
Simultaneous climate shifts have been reported in a “global” set of indices [2] . The indices studied were:
• SST3.4, central tropical Pacific SST;
• NPC1 (PDO), Pacific Ocean between 20N and 65N;
• SPC1, southern Pacific SST; and
• AMO, the subject of this paper.
The climate shifts since 1950 are found to occur during the following years [2] :
• 1956-1959;
• 1964-66;
• 1969;
• 1976-77;
• 1986-87 and
• 2001-02.
Weaker shifts were found at other dates. This list is consistent with climate shifts reported in prior papers.
A question arises: From what to what does the climate shift refer? The answer is: from one phase-locked state to another phase- locked state. See next section.
2.2. Phase-Locked States
Phase-locked states in the Earth’s atmosphere/ocean climate system have been reported by Douglass [5] . In this paper the El Niño/La Niña climate index SST3.4 was studied from 1870 to 2011. It was noticed in the record that certain time segments had sinusoid-like oscillations whose period appeared to be 2 or 3 years. A scheme novel to the application of climate indices was used to examine these time segments.
This scheme is called the autocorrelation method (AM) in which the autocorrelation of the time segment containing the 2- or 3-year oscillations was calculated vs. the delay time τ. If this time segment contains a sinusoidal signal of period 2 or 3 years then the autocorrelation function will have a maximum at a delay of 24 or 36 months. See [5] for details. This scheme not only gave the periodicity, it also gave the beginning and the end dates. These dates corresponded closely to the dates of previous determined dates of climate shifts. The periodicity of 2 or 3 years was identified with the 2nd and 3rd subharmonic of an annual forcing. Ten phase-locked states were reported.
In a later paper Douglass [6] studied five Pacific climate indices that included SST3.4. The other four were:
• SOI, the Southern Oscillation Index;
• MSU-lt, equatorial lower troposphere temperature;
• PDO, the Pacific Decadal Oscillation (29N to 65N); and
• W3, an equatorial wind index.
All showed phase-locked states over nearly the same date ranges as SST3.4.
2.3. Historic El Niño Episodes
Douglass, Knox, Curtis, Giese and Ray [7] (henceforth DKCGR) studied the index SST3.4 in more detail. They defined and enumerated the historical El Niño episodes. This more complete study reported 18 phase-locked segments since 1872. Within each segment sinusoid-like oscillations were observed and characterized by the 3 discreet indices described in [2] . The positive maxima of these oscillations were identified with an El Niño episode. A total of 40 episodes were found and are listed in their Table 2.
3. Data
AMO
The AMO index data are monthly values of the area weighted sea surface temperature (SST) of the Atlantic Ocean from the equator to 70N. The data range from 1856 to the present. We use the unsmoothed data set that can be downloaded at NOAA/AMO [8] .
SST3.4
We also use the monthly values of El Ñino/La Nina index, SST3.4, for comparison. This data is the same that was used to defined and find 40 historical El Ñino episodes. See [7] .
The data in both sets are measured in degrees Centigrade.
4. Analyzing AMO
4.1. Annual Term
The first task is to remove the strong annual signal from the AMO. A very common method used to “remove” the annual signal is the “climatology” scheme. However, this scheme can lead to false signals in the derived anomaly index [9] . Also, in this paper a better scheme was proposed. In this better scheme the data was operated upon with a filter f which is a running 12-month average. This filter removes not only the 1.0-year signal exactly but also the harmonics. The data after application of f on AMO are labeled
. (1)
As in [5] we also define a “high” frequency index:
(2)
Figure 1(a) shows a representative time segment of hAMO from 1990 to 2000. One sees a 1-year periodic signal of amplitude about 2˚C. The fact that the
(a) (b) (c)
Figure 1. (a) Annual component hAMO of the Atlantic Multi-decadal Oscillation, 1990-1999. See text for definition of hAMO; (b) Average of the monthly values of hAMO. The maximum occurs in August; (c) Autocorrelation of AMO_12, which is the dataset AMO with the annual component removed. The maximum at a delay of 68 - 72 years gives a new estimate of the multi-decadal period.
magnitude of the positive cycle is larger than the magnitude of the negative cycle means that there is a harmonic component.
Figure 1(b) shows the average of the monthly values of hAMO over the same time interval. From this plot one sees that the maximum in hAMO occurs between August and September.
Figure 1(c) shows the autocorrelation of AMO-12. The maximum at a delay of 68/72 years gives a new estimate of the multi-decadal period.
4.2. Multi-Centennial Trend
After the annual signal has been removed one sees that there is a multi-centennial trend that must be removed before an anomaly index can be defined. At this point in the analysis of AMO_12 many investigators consider the global warming hypotheses as part of the trend and propose ways to model it. Some use a linear trend which is seen not to describe the trend. There is, in fact, no need to explain the trend. Just remove it imperially. Along with aAMO _12 (red) in Figure 2(a) a second order polynomial fit is shown (black).
We define the AMO anomaly as
(3)
4.3. The Multi-Decadal Oscillation Period
Figure 2(b) shows aAMO (in red). The multi-decadal oscillation is shown by a running average over 120 months and is labeled aAMO_120. This curve is sinusoid-like. The positive values of aAMO_120 are usually called “warm” time intervals (Figure 2(b), black) which are
• 1863 to 1901;
• 1927 to 1965;
• 1997 to ? (see discussion below on this value).
These values are listed in Table 1 along with the values from Klotzbach et al. [10] .
What is the value of the multi-decadal period?
Figure 2(c) shows the autocorrelation of aAMO_12 vs. delay time. The delay time at the maximum is the oscillation period which is estimated as 68 to 72
Table 1. Warm time intervals of AMO.
(a)(b)
Figure 2. (a) AMO_12 and a polynomial fit, 20.6057 − 0.00130x + 0.0000263x2, where x is the number of years past 1860. (b) aAMO (red) and 120-month-filtered version aAMO_12 (black). “Warm” time intervals are indicated by tinted rectangles.
years. This new estimate is close to the value of 65 - 70 years given by [11] . The structure seen on this plot is due to the decadal effects that are considered in the next section.
When does the third warm time interval end?
Using the 68 - 72-year value for the period then adding it to the end of warm segment 2, the estimate for the end of segment 3 is 2030 - 2034.
Some have suggested that the third warm time segment has already ended [10] . Our analysis indicates that it has not.
5. Phase-Locked States in AMO_12
Here we do not need to use anomalies to determine phase-locked states. We compare AMO_12 to the previously published results on SST3.4 [7] , which showed phase-locked states. We have studied AMO_12 from the mid-19th century and have observed the phase-locked phenomena. However, only results from data since 1990 are shown.
Figure 3 shows the SST3.4_12 data (black) and the date range of each SST3.4_12 phase-locked state is indicated by a tinted rectangle. Figures 3(a)-(c) shows phase-locked states 9, 10 and 11 of SST3.4_12. In addition, the various El Niño episodes listed in DKCGR are indicated by their assigned number.
5.1. Phase-Locked Segments of AMO_12
Figure 3(a) shows both AMO_12 (red) and SST3.4_12 (black) from 1991 to 2002. This time interval contains phase-locked segment 9 of SST3.4_12 with a period of 3 years, which was determined by the AM method. The range was from about June 1992 to about January 1999. The three maxima are the three El Niño episodes contained in this phase-locked state. These are labeled 32, 33 and 34 by DKCGR. Episode 34 is better known as “the El Niño of 1997/1998”. We now examine AMO_12. One sees that this plot also shows a signal with 3 maxima spaced about 3 years apart. These maxima occur after the maxima in SST3.4_12. The delays are 17, 11 and 11 months, respectively.
The AMO_12 data were also analyzed by the autocorrelation method. When the beginning date was about April 1992 and the ending date was about April 2002 the plot in Figure 4 was obtained. The maximum of the autocorrelation function occurs at a delay of 36 months. Thus, this time segment has an oscillation with a period of 3 years.
Figure 3(b) shows AMO_12 and SST3.4_12 from 1999 to 2008. This time segment includes phase-locked segment 9 of SST3.4_12 of period 2 years. There are 3 El Nino episodes labeled 35, 36 and 37 spaced 2 years apart. The date range of this SST3.4_12 phase-locked state is from about January 2001 to about March 2009. One sees that there are also 3 maxima in AMO_12 that occur after the maxima in SST3.4_12. The corresponding lags are 9, 12 and 14 months, respectively.
Figure 3(c) shows the time interval containing phase-locked segment 11 of SST3.4_12. This segment began at about March 2009 and is ongoing [11] . The signature of this aSST3.4_12 phase-locked state is the three maxima spaced 3 years apart. The 3 maxima are El Niño episodes 38, 39 and 40. Douglass has predicted that a fourth maximum may occur during the boreal winter of 2018 [11] .
(a)(b)(c)
Figure 3. (a) SST3.4_12 (black), left scale, and AMO_12 (red). The phase-locked segment 9 of period three years is indicated by a tinted rectangle. The three numbered maxima of SST3.4 are the three El Niño episodes as identified in reference 7; (b) Same as (a) except for phase-locked segment 10 of period two years; (c) Same as (a) except for phase-locked segment 11 of period three years.
Figure 4. Autocorrelation of AMO_12 from April 1992 to April 2002. The maximum at 36 months verifies that AMO_12 has a period of three years during this time.
One sees that there are 3 maxima in AMO_12 that occur after the 3 maxima of SST3.4_12. The lags are: 7, 0 and 12 months, respectively.
5.2. AMO_12 Lags SST3.4_12
From Figures 3(a)-(c) one sees that each maximum in SST_12 is followed by a maximum in AMO_12. Table 2 lists the dates of these maxima and their differences. The average and standard deviation of the delays are 11 ± 5 months.
We also calculated a delayed correlation plot of aAMO_12 vs. aSST3.4_12 from 1950 to 2016 (not shown). We found a maximum in the correlation coefficient of 0.77 when the aAMO_12 delay was 11 months. The associated relationship can be quantified as follows:
.
One sees that decadal changes in aAMO_12 are about 10 times smaller than changes in aSST3.4_12. Thus, these two methods of estimating delay are consistent with each other.
6. Discussion
Muller, Curry et al. [3] (MC2013) studied the correlation of land temperatures vs. various other surface temperature climate indices since 1950 also studied other indices. Two of them, AMO and SST3.4, are the same as studied in more detail here. Their findings that concern us are as follows:
1) A high correlation of SST3.4 and AMO with AMO lagging by a few months.
Table 2. Time Delay of AMO_12 after SST3.4_12.
2) A claim of a strong signal in AMO of a period of 9.1 ± 0.4 years.
MC2013 show a plot of spectral power vs. the frequency in cycles/year for AMO. This plot has a narrow maximum at about 0.11 cycles/year (9.1 years) and a second and third maxima but are slightly smaller maximum at about 0.28 and 0.34 cycles/year (3.5 and 2.9 years), which we will label maxima 1, 2 and 3. We will call this combination of three maxima the “signature”.
7. Comments on MC2013
1) Our result of AMO lagging SST3.4 by 11±5 months is consistent with the findings of MC2013.
2) Since our more detailed study of AMO found no signal of period 9.1 years we did a Fourier analysis of AMO. We calculated the Fourier spectrum of AMO using a fast Fourier transform (FFT) and obtained a plot with the same signature as reported by MC2013 and with the same 3 maxima. Thus, we agree. What is the explanation? While it is true that a periodic signal will result in a maximum in the Fourier spectrum the converse may not be true. In this case there is no 9.1 year periodic signal in AMO corresponding to the maximum in the Fourier spectrum. Also, if there were a signal of 9.1 years in ZMO there would be a maximum in the delayed autocorrecting function at a delay of 9.1 years. We calculated the autocorrelation function of AMO from 1950 to 2017 (the same range as in MC2013). We obtained a plot (not shown) almost the same as in Figure 1(c) with no maximum at a delay of 9.1 years.
We conclude that the claim of MC2013 of a 9.1 year periodic signal is an error in interpretation of a maximum in the Fourier spectrum of AMO.
8. Conclusions
We have measured the AMO multi-decadal period to be 68 - 72 years which is close to the value of 65 to 70 years found in the literature. We also estimate that the present “warm” time interval will not end until about 2031/2035.
More importantly, we find that the AMO has time segments showing the same phase-locked states of periods 2 and 3 years as reported in many other climate indices. However, the AMO lags the SST3.4 by about a year.
In addition, we conclude that interpreting a maximum in the Fourier spectrum of AMO as an oscillation in AMO of period of 9.1 years is incorrect.
Acknowledgements
We wish to acknowledge helpful comments from P Klotzbach and R.S. Knox.