1. Introduction
Recent global warming attribution studies have shown that the global temperature trend is substantially due to changes in the value of anthropogenic forcings, while its interannual or decadal variability is influenced by the natural variability modes of the climate system, the most impactful of which appears likely to be the El Niño Southern Oscillation (ENSO) [
1,
2].
Although it is obviously being investigated what perturbative effects anthropogenic forcings may have on the course of these natural cycles (see, for instance, Chapter 3 of [
1]), their typical behaviour is always ascribed to the natural dynamics of the climate system. Recently, however, this has been called into question with regard to the Atlantic Multidecadal Oscillation (AMO) (see
Figure 1), a cycle in the Atlantic Ocean that influences climate and meteorological phenomena in several regions of the world: see, for instance, [
3], for a recent example of these impacts.
A good number of AMO “attribution” studies have been carried out in the past. Papers on this topic can be divided into two broad strands: those that claim a fundamental role for factors internal to climate dynamics (including ocean circulation) in driving the behaviour of the AMO time series and those that see external forcings, especially sulphate aerosols, as influential, particularly over the last century. Thus, is AMO a genuine representation of the internal natural variability of the climate system or is substantially driven by external factors? This issue is much debated and there is no consensus on the answer to the previous question in the scientific community.
In this framework, here, we employ data-driven methods for investigating this topic. After some preliminary applications to the study of the influences on AMO by other climate modes of internal variability, we focus on the possible driving roles of the external (natural or anthropogenic) forcings. In doing so, we present an analysis that can complement the past ones, performed almost exclusively through Global Climate Models (GCMs), and can contribute to the debate on AMO “attribution”.
Kushnir [
4] was probably the first to analyse quasi-periodic variations in the North Atlantic on a multi-decadal scale, using about a century of temperature and pressure observations. Meanwhile, Schlesinger and Ramankutty [
5], performing a spectral analysis of global surface air and sea temperature data, showed that the ~65-year oscillation is mainly the North Atlantic mode: it occurs only in this area. Thus, the term AMO implies not only the Atlantic origin of this oscillation but also its low-frequency character. Further details on the variability of low-frequency AMO have been described by Polonskii and Voskresenskaya [
6].
But by what mechanism are these quasi-periodic AMO oscillations produced? The first clue was certainly the meridional heat transport (MHT) in the North Atlantic. However, there are different views on the main causes of AMO generation and some authors have shown that AMO could be generated without any active ocean dynamics.
It is clear that, in the complex geophysical system, different mechanisms can act simultaneously. Here, however, our aim is to analyse which single influences can be the driving factors for AMO.
To our knowledge, the studies that have faced this problem of AMO “attribution” were published in [
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21]. The first papers focused on the role of natural changes in oceanic circulation [
7,
8]. Then, volcanic forcings were considered as the main drivers of AMO [
9] and a role for anthropogenic aerosol emissions came to be considered as very relevant for the twentieth century [
10]. Zhang et al. [
11] found discrepancies that cast doubts on previous claims in [
10] but without any alternative explanation for the AMO behaviour.
Later, Knudsen et al. [
12] suggested that the Atlantic Meridional Overturning Circulation (AMOC) could be crucial as a link between external forcing and North Atlantic sea-surface temperatures, a conjecture that is able to reconcile previous opposing theories concerning the origin of the AMO. Clement et al. [
13] showed that AMO is forced by mid-latitude atmospheric circulation, while Bellucci et al. [
14] studied the 1940–1975 North Atlantic cooling and found a key role for anthropogenic forcings in driving it. In the meantime, Cane et al. [
15] showed that the ocean is influent on AMO, even if it seems not so important for the simulation of the surface temperature climate variability. Murphy et al. [
16], via a modelling study, showed a clear role for external forcings in driving the observed AMO and argued that, very likely, internal variability is insufficient to drive the AMO multidecadal changes observed over the last century. A new paper, [
17], promptly confirmed this latter result and showed that greenhouse gases (GHGs) and tropospheric aerosols were probably the main drivers of the AMO in the latter part of the twentieth century.
On contrast, other results [
18,
19,
20] suggested internal sources for the AMO behaviour, e.g., oceanic variability and circulation. Finally, it should of course be noted that the AMO signal has been estimated for centuries and millennia well before the industrial era: see, for example, [
22]. Therefore, it is clear that the above-mentioned studies may seem rather limited in establishing the unforced or forced nature of AMO. In this regard, however, a paper by Mann et al. [
21] should be noted, where they addressed this issue and showed how, throughout the last millennium, this oscillation may be due to pulses of volcanic activity.
Here we stress again the different results found and that the vast majority of studies have been performed through dynamical analyses via GCMs. In this research framework, however, data-driven methods can contribute to analysing the problem via complementary approaches. This has been performed, for instance, in [
23], where one of us (A.P.) applied a neural network tool [
24] to this analysis, finding a clear role for anthropogenic forcings (especially sulphate aerosols) in driving AMO behaviour.
In any case, as already shown for global warming attribution studies, also other data-driven methods can be applied to analyse causality relationships in the climate system. See, for instance, some of our previous papers on this topic [
25,
26,
27]. Thus, here we perform an analysis of the influence of external (natural and anthropogenic) forcings on the AMO behaviour by means of a linear Granger causality analysis and by a nonlinear extension of this method, with the aim of contributing to the scientific debate on this topic.
3. Data Description
In this investigation, even for a consistent comparison with the previous paper by Pasini et al. [
23], we consider annual data for the AMO index and the following radiative forcings: RFWARM, which is the warming part of the anthropogenic forcing (GHGs + BC RFs), its cooling part represented by sulphates (RFSOX), the RF of the solar activity (RFSOLAR) and the RF of volcanic emissions (RFVOL). In the following analysis, we also use the sums RFNAT (RFSOLAR + RFVOL), RFANTH (RFWARM + RFSOX) and RFTOT (RFNAT + RFANTH).
In this paper, we consider the AMO index calculated with a standard method [
30]: its time series is freely available at
www.esrl.noaa.gov/psd/data/timeseries/AMO (accessed on 1 March 2024). Data about the anthropogenic radiative forcings are downloaded from the dataset collected at
http://www.sterndavidi.com/datasite.html (accessed on 1 March 2024, see also the paper by Stern and Kaufmann [
31]). In the present paper, the data on GHG concentrations come from the NASA/GISS website and their RF calculations are performed through well-known formulas [
32,
33]. The global estimates of sulphate emissions in the past [
34,
35] and the computation of direct and indirect RFs as in [
31]—which is based on slight modifications of previous studies [
36,
37]—supply us with data about the radiative forcing of sulphates. These last data are available until 2007: to consider a prolonged time series but without a too long extrapolation, we continue this series of RFSOX with constant data until 2011. Past data about the RF of black carbon come from the RCP8.5 scenario [
38]. For data concerning natural radiative forcings, solar irradiance is approximated by an index previously assembled [
39], and its data are downloadable from
https://data.giss.nasa.gov/modelforce/solar.irradiance/ (accessed on 1 March 2024). The conversion from solar irradiance to RFSOLAR is obtained in a standard way [
33]. Optical thickness data [
40], available at
https://data.giss.nasa.gov/modelforce/straaer/ (accessed on 1 March 2024), are considered as a proxy of the volcanic activity of dust emissions. RFVOL is set at 27 times the optical thickness [
31].
Preliminary Analysis
Preliminary to the Granger causality analysis, it is necessary to establish the order of integration of the involved time series. We remember that an integrated time series of order d, denoted as I(d), is a time series that have to be differenced d times to achieve stationarity.
To detect the order of integration of the variables, we use the Augmented Dickey–Fuller (ADF) test. To carry out the ADF tests, we estimate the following auxiliary regressions for each variable of interest y:
where Δ is the first difference operator, u
t∼
WN(0, σ
2u).
In this parametric framework,
yt is non-stationary (has a unit root) if
c = 0. We test
H0:
c = 0 against
H1:
c < 0 using the test statistic
where
is the ordinary least square estimate of
c and
is its standard error. Under
H0:
c = 0, this statistic follows asymptotically a Dickey–Fuller distribution. We reject
H0 if
ADFt is less than the critical value.
Applying the ADF test, we have obtained the following results: the series AMO, RFSOLAR, RFVOL, RFNAT, RFTOT are I(0), while the series RFSOX, RFWARM, RFANTH are I(1): see
Table 1.
4. Results
Although the main objective of this paper is to understand whether AMO is forced or unforced by external forcings, obviously a causality analysis between various modes of the natural variability of the coupled atmosphere–ocean system and AMO is of interest and is, as far as we know, a novelty when seen as an application of our methods.
Thus, referring to the data analysed in previous papers published in the literature, such as [
42], we applied our methods to the study of influences from AMOC to AMO and from NAO to AMO. The results (not explicitly reported here but available on request) show that these natural variability modes have no causal role on AMO, as resulting from both the linear Granger analysis and the nonlinear Diks–Panchenko test.
Our results about the influences of external forcings on AMO give quite clear indications too. Firstly, natural forcings do not show any causal role on AMO in both linear and nonlinear analyses. The
p-values are always very high, even in the case of the sum of RFSOLAR and RFVOL in RFNAT (see
Table 2,
Table 3,
Table 4,
Table 5,
Table 6,
Table 7,
Table 8,
Table 9 and
Table 10), and the results are insensitive to choices of
Lx,
Ly and
ϵ.
On the other hand, in the case of anthropogenic forcings the
p-values are generally smaller. However, RFSOX does not show a causal role on AMO (see
Table 11,
Table 12 and
Table 13), while RFWARM does, at 5% significance in the linear case (see
Table 14). In the nonlinear tests, RFWARM also does not show a causal effect, except in one case for
Lx =
Ly = 3 and epsilon equal to 1 (see
Table 15 and
Table 16). When RFWARM and RFSOX are summed up in RFANTH, the causality remains in the linear Granger causality (see
Table 17), and this is even the case when considering the total forcing RFTOT (see
Table 18). In the nonlinear tests, the null hypothesis of noncausality is always confirmed even in the latter cases.
How should we read these results in the framework of our previous knowledge and studies? Certainly, on the one hand, our work confirms the results of those previous studies, such as the only other data-driven investigation [
23], which showed that anthropogenic forcings have a certain influence on AMO trends. On the other hand, it is not the influence of RFSOX that is important here but that of RFWARM. In both [
23] and the present study, the negligible influence of natural forcings (and other modes of natural variability) and the importance of anthropogenic ones is confirmed, albeit in a different manner. Furthermore, whereas in [
23] it was the nonlinear model that showed the influence of anthropogenic variables on AMO, here causality is found almost exclusively in Granger’s linear analysis.