Addendum to: Nature https://doi.org/10.1038/s41586-020-2641-xPublished online 02 September 2020
In this Article, we demonstrated that the El Niño–Southern Oscillation (ENSO) response to greenhouse warming is self-modulating and depends on its historical variability, partly by a negative interexperiment correlation between initial ENSO variability and its initial-to-future change. Since publication, it has come to our attention that initial variability correlated with the difference between itself and its future value will tend to be, by definition, anti-correlated. Furthermore, an anti-correlation between initial and future variability can be amplified by using a normalized time series, such as the normalized E-index described in our Article. For example, when the standard deviation in the initial period is amplified by internal variability, this also amplifies the standard deviation for the whole time series (assuming an average level of variability after the initial period). This increases the normalization coefficient, and thus reduces the normalized standard deviation for the future period. As such, the strong anti-correlations in Fig. 1b, c and 4 of our original Article are, in part, contributed by the statistical formulation.
Here we reformulate the analysis to avoid these statistical issues. By using raw sea surface temperature (SST) anomalies from large ensemble experiments of the three different models—the Community Earth System Model (CESM), the Geophysical Fluid Dynamics Laboratory Climate Model version 3 (GFDL-CM3) and Geophysical Fluid Dynamics Laboratory Earth System Model second model (GFDL-ESM2M))—we show that the connection between the twentieth-century and twenty-first-century ENSO variability is still present (Fig. 1). The variability of quadratically detrended SST anomalies at the geographical centre of eastern Pacific ENSO shows that experiments that generate greater initial variability produce smaller future variability and therefore a smaller increase. The relationship is statistically significant above the 95% or 90% confidence level (Fig. 1). The correlation from the direct SST anomalies (r = −0.33 to −0.42) avoids any statistical artefacts and is fundamental to the result expressed in Fig. 1b, c and 4 of the original Article. The interexperiment relationship in our Article indicates that the initial variability explains about 60% of the future-minus-initial change in variability. With the revised approach, the initial variability explains 10–18% of future variability, with the remainder being attributed to other factors including the effect of stochastic forcing, such as atmospheric noise, which is compounded throughout the evolution of the twenty-first century. The butterfly effect is expected to operate among the members of the Coupled Model Intercomparison Project (CMIP) multi-model ensemble, despite the large inter-model differences in the ENSO variability due to various factors such as model physics and climate sensitivity. This is underscored by the comparable inter-model spread in the initial-to-future trend of ENSO variability between the CMIP multi-model ensemble and the single-model CESM large ensemble—the future-minus-initial standard deviation in Niño3.4 (E-index) variability for the CMIP ensemble is 0.27 °C (0.32 s.d.), and 0.18 °C (0.24 s.d.) for the CESM-LE—thus indicating a substantial contribution of internal variability to the CMIP inter-model differences in the ENSO variability trend.
The mechanisms presented in our original Article are also robust regarding the use of raw data. For example, the nonlinear relationship between the eastern Pacific SST anomalies and the net heat flux damping that underpins ENSO nonlinear thermodynamics shown in Fig. 2a of this Amendment is similar to that in Fig. 2a of the original Article. The relationship indicates that thermal damping during an El Niño is stronger than that during a La Niña, and represents an anomalous oceanic heat loss and atmospheric heat gain, respectively. Thus, experiments with stronger initial SST variability produce stronger cumulative heat loss in the central-to-eastern equatorial Pacific over the initial 50 years (see Fig. 2b, c; compare with Fig. 2d, e of the original Article), and the anomalous cooling eventually reaches the upper equatorial Pacific (see Fig. 2d; compare with Fig. 3c of the original Article), leading to a smaller warming. This anomalous cooling offsets the greenhouse-warming-induced intensification of stratification in the upper ocean that otherwise facilitates a stronger increase in ENSO variability.
In summary, although our main conclusion on the link between twentieth-century and twenty-first-century ENSO variability remains qualitatively consistent, the magnitude of the relationship is considerably reduced from what the original anti-correlations between initial variability and its future change might suggest. We thank J. Räisänen for bringing this to our attention, and are grateful for his open, detailed and constructive feedback. The original Article has not been corrected.
About this article
Cite this article
Cai, W., Ng, B., Geng, T. et al. Addendum: Butterfly effect and a self-modulating El Niño response to global warming.
Nature (2021). https://doi.org/10.1038/s41586-021-03261-4