Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Continuation of tropical Pacific Ocean temperature trend may weaken extreme El Niño and its linkage to the Southern Annular Mode

### Subjects

A Publisher Correction to this article was published on 04 February 2020

## Abstract

Observational records show that occurrences of the negative polarity of the Southern Annular Mode (low SAM) is significantly linked to El Niño during austral spring and summer, potentially providing long-lead predictability of the SAM and its associated surface climate conditions. In this study, we explore how this linkage may change under a scenario of a continuation of the ocean temperature trends that have been observed over the past 60 years, which are plausibly forced by increasing greenhouse gas concentrations. We generated coupled model seasonal forecasts for three recent extreme El Niño events by initialising the forecasts with observed ocean anomalies of 1 September 1982, 1997 and 2015 added into (1) the current ocean mean state and into (2) the ocean mean state updated to include double the recent ocean temperature trends. We show that the strength of extreme El Niño is reduced with the warmer ocean mean state as a result of reduced thermocline feedback and weakened rainfall-wind-sea surface temperature coupling over the tropical eastern Pacific. The El Niño-low SAM relationship also weakens, implying the possibility of reduced long-lead predictability of the SAM and associated surface climate impacts in the future.

## Introduction

The Southern Annular Mode (SAM) is the leading mode of variability of the Southern Hemisphere (SH) extratropical circulation on weekly and longer timescales that describes a meridional vacillation of the eddy-driven jet and associated storm track1,2,3,4,5. The spatial pattern of the positive polarity of SAM (high SAM) is characterised by a nearly zonally symmetric annular pattern of positive anomalies of pressure/geopotential height in the SH midlatitudes and negative anomalies in the Antarctic region, associated with a poleward shift of the eddy-driven jet and storm track3. The negative polarity of SAM (low SAM) is characterised by an equatorward shift of the eddy-driven jet and associated storm track. SAM, which is intrinsic to the troposphere, exhibits a decorrelation time of about two weeks6 and can be realistically simulated in an atmospheric general circulation model forced with climatological sea surface temperatures (SSTs) at the lower boundary7,8. In addition to the well established trend toward the higher polarity phase of SAM in response to the Antarctic ozone depletion, year-to-year variations of seasonal mean SAM during spring and summer can be promoted by the El Niño-Southern Oscillation (ENSO)4,9,10,11,12,13,14,15 and by downward coupling of anomalous conditions in the Antarctic stratospheric polar vortex that develops as early as austral winter16,17,18.

The relationship between ENSO and SAM is particularly important for a long-lead prediction of  SAM and its surface climate anomalies because ENSO can be skillfully predictable at lead times of 2–3 seasons and beyond4,19. During austral spring and summer, between 10–36% of the variance of SAM is explained by its relationship with eastern Pacific-type (EP) ENSO4,9,10,13. The warm phase of ENSO (El Niño) is associated with low SAM and the cold phase of ENSO (La Niña) is associated with high SAM. Seasonal mean SAM during austral spring and early summer is predictable with an atmosphere-ocean coupled forecast system at up to 6 month lead time, and the long-lead predictive skill is likely to be attributed to the relationship of SAM with eastern Pacific ENSO4.

The predictability of SAM stemming from ENSO leads to predictability of SH extratropical atmosphere, ocean and sea-ice variations in regions which are strongly influenced by SAM13,20,21,22,23,24,25. This is in addition to the predictability arising from direct impacts of ENSO26. Because of the potential benefits of predicting extratropical climate as a result of the co-variation of SAM with ENSO, it seems natural to question how the ENSO-SAM relationship will change in a future climate given that there has been much focus on how ENSO and its impacts might change in warmer climte27,28. Unfortunately, there is great uncertainty in how this relationship may change in the future because, for instance, the models used in the Climate Coupled Model Intercomparison Project (CMIP5)29 fall significantly short in simulating the observed linkage between ENSO and SAM in the current climate, demonstrating large inter-model spread all year round30. Therefore, an idealised model experiment would be useful to explore how SAM might respond to ENSO in a future climate.

In this study, we specifically focus on the response of SAM to extreme El Niño in a future waremer climate because characteristics of extreme El Niño and extreme La Niña31 are not entirely symmetric and their impacts on the atmosphere are not exactly opposite to each other32,33. This focus on extreme El Niño is also to increase the signal in the model experiments. We explore our research question by superimposing 1982, 1997 and 2015 El Niño oceanic conditions on a hypothetically warmer mean state, which we derive as a continuation of the observed ocean temperature trends since 1960. This is done by using a series of seasonal forecast sensitivity experiments with the global coupled model seasonal prediction system POAMA (Predictive Ocean and Atmosphere Model for Australia)34.

The pattern of the observed ocean surface temperature trends of 1960–2014 is characterised by significant surface warming trends in the tropical western Pacific and Indian Oceans together with a slight surface cooling trend along the equatorial eastern Pacific35,36,37,38,39 (Fig. 1a). In the subsurface along the equator, there has been a significant cooling trend along and below the thermocline in the equatorial Pacific (Fig. 1b). This overall trend pattern is often described as a “La Niña-like” mean state change for convenience40. This trend pattern and magnitude in SSTs is well captured by the 1st empirical orthogonal function (EOF) eigenvector of decadally smoothed SST variability, which is well separated from the 2nd EOF eigenvector that represents the Inter-decadal Pacific Oscillation (IPO)/Pacific Decadal Oscillation (PDO)41 (Supplementary Fig. S1). This observed trend pattern is distinctively different from the trend patterns simulated by the CMIP multi-model mean historical and future projections in responses to increasing greenhouse gases42. The projected trend shows greater warming over the tropical eastern Pacific compared to its surrounding oceans and flattening of the thermocline in the equatorial Pacific32,42,43 (i.e. a “El Niño-like” mean state change).

An El Niño-like mean state change in response to increasing greenhouse gases is explained by weakening of the Pacific Walker circulation and the zonal SST gradient in the tropical Pacific, consistent with the energetic and hydrological balances under increasing greenhouse gas forcing44,45,46. The amplitude of El Niño is closely linked to the degree of mean warming in the eastern equatorial Pacific in the CMIP5 models, the majority of which show El Niño-like mean state change patterns42. Consequently, a significant increase in the frequency of extreme El Niño is projected in the 21st century compared to the 20th century28. However, the fidelity of this El Niño-like simulated mean state change in response to greenhouse gas forcing and its potential impact on ENSO characteristics have been vigorously debated36,39,47,48,49,50,51. The El Niño-like mean state change in some models has been attributed to systematic model biases such as a too regular and symmetric El Niño-La Niña cycle together with ocean mixed layers that are too deep in the tropical Pacific Ocean51,52; too weak inter-basin teleconnections in the tropical oceans53; and a too-cold cold tongue in the equatorial Pacific with too high relative humidity and too low wind speed39. On the other hand, the ocean dynamical thermostat mechanism40,54, the non-linear ENSO warming suppression mechanism55, and/or the inter-basin warming contrast mechanisms56,57 suggest that equatorial eastern Pacific SST would warm more slowly than the equatorial western Pacific SST as the Earth warms. Therefore, the possibility cannot be ruled out that the observed La Niña-like ocean temperature trend has been in part forced by global warming39,56,58, and it may continue into the future. Thus, we suggest that it is a valid and valuable question to address how the characteristics and teleconnection of extreme El Niño will change in the future if the observed La Niña-like long-term mean state warming continues, a scenario which has not yet been given much attention.

We address this premise by conducting four forecast sensitivity experiments whereby the ocean initial conditions are altered in a manner following the approaches of earlier studies using POAMA59,60. These experiments are summarised in Table 1. The present climate El Niño experiment (pElNiño) was conducted by initialising the coupled forecast model POAMA with observed ocean initial conditions for 00 UTC 1 September of 1982, 1997 and 2015 that represent the developing stages of the three extreme El Niños observed in the modern instrumental record43. We also initialised POAMA with the mean ocean conditions for 00 UTC 1 September computed over 1981–2013 in order to produce a set of climatological forecasts for the present climate (the present climatology experiment; pClim). To explore changes to extreme El Niño events and their linkages to SAM in the climate warmed up by enhanced observed ocean temperature trends, we computed the trends in ocean temperatures and salinity at all available vertical levels, latitudes, and longitudes on 00 UTC 1 September for the period of 1960–2014. We then doubled their magnitudes in order to produce a stronger impact in our modeling framework (Supplementary Fig. S2). We added these doubled trends to the observed ocean initial conditions used in pElNiño to generate forecasts of El Niño that occur with a warmer ocean mean state (the warmer climate El Niño experiment; wElNiño). Finally, we generated the climatological forecasts of the warmer climate by adding the doubled ocean temperature and salinity trends to the climatological ocean initial conditions used in pClim. We refer to this as the warmer climatology experiment, wClim. We denote present and warmer climate predicted El Niño with respect to their respective present and warmer climatologies as pElNiño’ and wElNiño’, respectively.

The atmosphere and land component models of POAMA were initialised with 33 different conditions drawn from observed states during 1981–2013 in order to scramble atmosphere and land initial conditions within the observed range so to generate ensemble forecasts. The CO2 concentration was fixed to 345 ppm, which is the default value of POAMA retrospective forecast set, and the ozone concentration was prescribed by the observed monthly climatology34. Therefore, differences in the predicted atmospheric circulation in our experiments can be interpreted as responses to the differences in ocean forecasts in the experiments. We initialised all experiments on 1 September and limited our interest to the three month mean forecasts for October to December (OND) because ENSO and SAM are better correlated on a seasonal time scale than on shorter time scales and the three month mean observed correlation is maximum in OND4.

## Revisiting the Observed Relationship Between ENSO and SAM

Austral spring-summer SAM is influenced by both ENSO and Antarctic stratospheric vortex variation, the latter of which is considered to be a stronger driver16. SAM has also shown a strong positive trend from the 1980s to the late 1990s due to the anthropogenically-driven Antarctic ozone depletion61. To highlight the relationship of SAM with ENSO, which is largely an interannual variation, we de-trended the SAM index for the period of 1979–2016. Then, we regressed out the component of SAM related to the stratospheric polar vortex variation from the de-trended SAM index, using the SH stratosphere-troposphere coupled mode index18 as a predictor representing anomalous stratospheric conditions (i.e. de-trended residual SAM; see Methods for the climate indices used in the study).

Earlier studies18,62 reported no significant relationship between the Niño3.4 SST index, which represents eastern Pacific ENSO variability, and the Antarctic stratospheric vortex variation during the last 40 years, but SAM is significantly better correlated with Niño3.4 SST after the removal of the influence of the SH stratospheric polar vortex variation from the SAM index (r = −0.53, statistically significant at the 0.1% level (i.e. p < 0.001); Fig. 2c). The three extreme El Niño events of 1982, 1997 and 2015, which were accompanied by a substantial spread in the magnitude of raw SAM (e.g. SAM was neutral during the late spring of the strong 2015 El Niño), all exhibited strong low SAM after the removal of the stratospheric influence. This result supports the connection between ENSO and SAM in OND.

Figure 3 shows the SST anomalies of the three extreme El Niño events and de-trended residual MSLP anomalies of OND. Although there are differences in the spatial details and magnitudes of anomalies of SSTs, all three El Niño events are associated with low pressure anomalies equatorward of 60°S with a wavenumber 3 structure and high pressure anomalies poleward of 60°S, which typify low SAM. Strong stationary Rossby wave propagation from the Maritime Continent poleward and eastward to the Amundsen-Bellingshausen Seas is a distinctive teleconnection driven by El Niño, which is likely to further contribute to the amplitude of low SAM by placing a strong low pressure anomaly centre north of 60°S and a strong high pressure anomaly centre slightly south of 60°S. As noted in Fig. 2a, SAM was observed to be neutral despite the extreme strength of El Niño in OND 2015 because anomalous Antarctic stratospheric vortex strengthening18 and associated ozone depletion were acting to produce high SAM (https://ozonewatch.gsfc.nasa.gov/meteorology/figures/merra/ozone/toms_areas_2015_omi+merra.pdf). After removing this influence of the stratospheric polar vortex, the residual SAM is strongly negative with zonally symmetric pressure anomalies like those of 1982 and 1997 (Fig. 3 right panels). This implies that in 2015 the forcing of high SAM from a strengthening of the polar stratospheric vortex and associated depletion of Antarctic ozone may have completely countered the forcing of low SAM from El Niño in the OND season, thereby lowering the predictability of SAM and the associated surface climate.

The promotion of low SAM by El Niño has previously been shown to result from intensification and equatorward contraction of the the Hadley circulation because the El Niño SST anomalies act to increase near equatorial diabatic heating (latent heat release as a result of moist atmospheric convection), which is flanked by enhanced diabatic cooling (longwave radiation to space as a result of enhanced subsidence and reduced latent heat release) in the subtropics (Fig. 4a). Thus, westerlies on the equatorward side of the climatological subtropical jet strengthen (Fig. 4b). This increase in westerly winds closer to the equator allows deeper penetration of extratropical baroclinic Rossby waves into the tropics11,12. As a result, eddy momentum flux divergence in the tropics shifts closer to the equator, and momentum flux convergence anomalously increases in the midlatitudes (30–40°S), while decreasing in the higher latitudes (50–70°S; Fig. 4c). Extratropical baroclinicity and the associated storm track thus shift equatorward, which is manifest as low SAM. These general features resulting from the eastern Pacific SST variations during El Niño, which were proposed in the literature, are confirmed by the composites of the three extreme El Niños in Fig. 4 and by the anomalies of the individual extreme El Niño in Supplementary Fig. S3.

## Experiment Results

### El Niño and SAM in the present climate

Figure 5 displays 33-member ensemble mean anomalies of SSTs and MSLP averaged during OND for the present climate 1982, 1997 and 2015 El Niños. Anomalies are formed relative to the present climate as pElNiño’ = pElNiño-pClim. The strong El Niño events are skillfully predicted at this short lead time of 1 month, although the amplitude of 1982 El Niño is substantially underpredicted (by ~0.7 °C; Supplementary Fig. S4). In all three El Niño cases, overall patterns of forecast MSLP anomalies feature low SAM with low pressure anomalies being dominant equatorward of 60°S and high pressure anomalies over the polar cap, and a stationary Rossby wave train extending from the Maritime Continent poleward and eastward to the south eastern Pacific, which are consistent with the observed anomalies shown in Fig. 3. However, a strong low-pressure center observed over the southern Atlantic Ocean is missed in all three of the simulated El Niño cases (Fig. 5 right panels). Reproducibility of low SAM in pElNiño’, which is initialised with realistic ocean conditions but random atmosphere and land conditions, confirms that strong El Niño can be an important forcing of low SAM, and the polarity of SAM would have been more negative during El Niño of 2015 if the Antarctic polar vortex were not significantly stronger than usual.

### Warmer climate mean state

For the results of the forecasts of extreme El Niño in the idealised future warmer climate (wElNiño) to be scientifically reliable and attributable to the change in the mean state of the initial conditions, key features of the observed ocean temperature trends added in the initial conditions and resultant ocean circulation changes should be faithfully maintained through the October to December verification period. This is demonstrated by comparing mean differences in SSTs, equatorial Pacific upper ocean temperatures and tropical Pacific mixed layer circulations between wClim and pClim to the observed trends in OND over 1960–2014 in Fig. 6. We note that although we added doubled the observed ocean temperature trends to the climatological ocean initial conditions for wClim, the magnitudes of the ocean surface and subsurface warming over the following four months are more or less comparable to the magnitudes of the observed trends of the 55 years of 1960–2014. This loss of amplitude in the mean state trend may be partly due to the scrambled atmosphere initial conditions and their damping effect until the atmosphere and the ocean are brought to a balance in the first month of forecasts24, or the lack of changing the greenhouse gas forcing in the warmer climate runs. Nevertheless, it is encouraging to see the overall similarities between the warmer minus present climate differences and the observed trends for the spatial structures of SSTs, equatorial Pacific subsurface temperatures and tropical Pacific mixed layer ocean circulations. For instance, the pattern correlation between Figs. 6a,b is 0.55 (p < 0.001 with 5760 grid values), and the enhanced zonal SST gradient between the tropical western Pacific (120–160°E, 5°S-5°N)40 and eastern Pacific (Niño3.4 region) is maintained in OND forecasts in wClim minus pClim although the gradient is simulated with only about one third of the observed magnitude (0.3 °C compared to 0.9 °C). The observed reduction of upwelling into the mixed layer (averaged 0–45 m depth) in the equatorial eastern Pacific is also faithfully captured by wClim forecasts (Figs. 6e,f). On the other hand, forecasts of SSTs and subsurface temperatures are significantly different from the observed trends in the far eastern Pacific east of 120°W, which should be borne in mind when we diagnose below how the change in the mean state acts to change El Niño growth.

The atmospheric response to the warmer climate (Fig. 7a) shows more intense warming in the upper troposphere in the tropics than in the higher latitudes of the SH, steepening the meridional temperature gradient between the tropics and the south pole. This change results in a poleward shift of the eddy-driven jet with anomalously increased eddy momentum flux convergence around 60°S but decreased eddy momentum flux convergence around 45°S (Figs. 7b,c). The associated MSLP change projects onto high SAM (standardized SAM index = 0.9), although the high-pressure anomalies in the midlatitudes are not as annular as that of conventional SAM (Fig. 7d). These atmospheric changes are similar to those projected by climate models in response to increasing greenhouse gases63 despite the contrasting tropical Pacific zonal SST gradients between wClim and climate change simulations. The similarity in the atmospheric responses despite differences in the degree of warming in the tropical eastern SSTs is because equatorial waves quickly spread anomalous temperature change across the entire tropics59.

### El Niño on the warmer ocean mean state

We now turn to the warmer climate El Niño compared to the present climate El Niño. Because forecasts of the SST and MSLP patterns of the three El Niño years in the present climate appear to be all very similar to one another, we will present results of a grand ensemble of 99 members by combining the anomalies for the 1982, 1997 and 2015 simulations.

When El Niño of the same extreme strength occurs on the “La Niña-like” warmer ocean mean state (wElNiño’ = wElNiño-wClim), its strength is simulated to be weaker than that in the current climate (pElNiño’ = pElNiño-pClim) over the tropical eastern Pacific (p < 0.05, see Methods for statistical significance calculation; Fig. 8a, Supplementary Fig. S4). Therefore, the maximum warming of El Niño in the warmer mean state appears to be more confined to the central Pacific, giving it some central Pacific flavour. Analysis of advective feedback terms in the ocean mixed layer heat budget (averaged over 0–45 m, 5°S-5°N; see Methods) suggests that the reduced mean upwelling in the eastern Pacific in the warmer ocean mean state (Fig. 6f) is partly responsible for the weaker growth of eastern Pacific SST anomalies in the warmer climate (Fig. 8b). Across the central to western Pacific, the reduced strength of mean westward surface currents in the warmer ocean mean state (Fig. 6f) also contributes to weaker growth of warm anomalies over and west of the dateline (Fig. 8c), while weaker anomalous eastward currents during El Niño in the warmer climate contributes to weaker warming over the central Pacific (Fig. 8d). The cause of the weaker anomalous eastward currents during El Niño in the warmer climate is addressed below. In contrast, the enhanced mean zonal temperature gradient of the warmer ocean mean state over the equatorial Pacific (e.g. Figs. 6b,d) positively contributes to the growth of positive SSTs during El Niño over the central to western Pacific (Fig. 8e), thereby compensating for some of the negative contributions caused by the weakening of the mean westward zonal currents and the weaker El Niño-generated anomalous eastward currents there.

The cause of the weaker El Niño-generated anomalies in the waremer mean state can be traced to westward shifts of the maximum rainfall and wind responses to the El Niño SST anomalies while reducing their responses in the tropical eastern Pacific. These changes are quantified by the regressions of rainfall and 10-m zonal wind anomalies onto the Niño3.4 SST anomalies during OND (all 99 forecast members are used for the regression calculation). Figures 9a–c show that El Niño-driven rainfall in the warmer climate is shifted westward relative to that in the present climate, and over the Niño3.4 region the rainfall anomaly is weaker (p < 0.05). This likely occurs due to the cooler equatorial mean SSTs (Fig. 6b) and weaker El Niño SST anomalies caused by the reduced thermocline feedback in the eastern Pacific (Figs. 8a,b), which makes a SST anomaly less efficient at driving a rainfall/wind response60,64,65. In contrast, there is an increased rainfall response over and west of the dateline during wElNiño’ compared to pElNiño’ likely boosted by the warmer SST mean state of the tropical central to western Pacific, shifting the location of the maximum rainfall associated with the Niño3.4 SSTs about ~10° westward (p < 0.05). Likewise, SST anomalies over the Niño3.4 region induce weaker westerly response over the eastern Pacific but stronger westerly response west of the dateline with the warmer mean state compared to the present mean state (p < 0.05; Fig. 9d–f). Together these westward shifts of the zonal wind and rainfall responses during El Niño on the warmer mean state feed back into producing a weaker El Niño in the eastern Pacific.

Interestingly, these 1) weakened thermocline feedback and 2)air-sea coupling strength and resultant weakened El Niño with its maximum SST warming and rainfall response concentrated in the central Pacific are consistent with what has been diagnosed to have occurred during the early 2000s as compared to the 1980s and 1990s and is thought to reflect the shift to the cold phase of the IPO60,64,65,66,67. The cold phase of the IPO is also characterised by an enhanced zonal SST gradient across the tropical Pacific and shallower thermocline in the equatorial eastern Pacific, but there are many differences between the cold phase of the IPO and our warmer ocean mean state (Supplementary Fig. S1). Hence, the consistency in the changes in the ocean and atmosphere feedback and El Niño properties between our warmer ocean mean state and the cold phase of IPO highlight the role of the equatorial Pacific mean state in modulating the amplitude of El Niño and the longitude of maximum warming.

### El Niño and SAM on the warmer ocean mean state

Although El Niño weakens on the warmer ocean mean state, it still promotes low SAM. However, the El Niño-low SAM connection is substantially weaker in the warmer climate compared to the current climate as judged by standardized SAM being −0.7 as derived from the pressure anomalies of Fig. 10d compared to standardized SAM being −2.1 as derived from those of Fig. 10c (this weakening of the SAM strength is statistically significant at the 1% level). The Rossby wave train emanating from the Maritime Continent is also weaker in the warmer climate. This weakening of teleconnection is likely due to both weakening of El Niño and weakening of the tropical convective heating response to El Niño on the warmer mean state. As shown in Figs. 10a,b, wElNiño’ induces less intense warming over the tropical upper troposphere, and therefore, the meridional temperature gradient between the tropics and the SH midlatitudes is not as steep as that caused by pElNiño’, resulting in weakened westerly anomalies associated with wElNiño’ compared to those with pElNiño’s (Figs. 10e,f). This reduced change in the mean winds in the upper troposphere leads to reduced changes in eddy momentum flux convergence anomalies in the subtropics to the high latitudes of the SH (Fig. 10g,h); which feed less of anomalous momentum back to the mean flow, resulting in weaker dipole wind anomalies in the SH extratropics, leading to weaker low SAM. The strengths of the extratropical zonal-mean zonal wind dipole and the eddy momentum flux convergence dipole in the upper troposphere are statistically significantly different between pElNiño’ and wElNiño’ at the 5% level (Supplementary Fig. S8). Consequently, predictability of SAM in OND would be substantially reduced during El Niño under our scenario of continuation of the observed ocean temperature trends to the future.

## Concluding Remarks

This study was motivated by earlier research findings that ENSO is an important source of predictability of SH extratropical climate through its connection to SAM during austral spring and summer; yet this ENSO-SAM relationship is not skillfully simulated by climate models, therefore making it hard to foresee any possible change in this relationship in a future warmer climate. Furthermore, climate models have not been able to reach a strong consensus on how ENSO amplitude will change43, which would be a key determinant in predictability of extratropical climate via its teleconnections. The strength of ENSO and its spatial characteristics influence and are influenced by the tropical ocean mean state, especially the tropical Pacific Ocean mean state60,66. There is still vigorous debate, though, about whether the response of the tropical Pacific to global warming will be more “La Niña-like” with greater warming over the tropical western Pacific than over the eastern Pacific, which is what has been observed over the past 60 years39, or more “El Niño-like” with the opposite pattern of warming, which is projected by the majority of the CMIP5 models47.

In this study, we limited our focus to El Niño and attempted to address this question: How will extreme El Niño and its relationship with SAM change in a warmer climate if the “La Niña-like” ocean temperature trends that have occurred over the past 60 years continue into the future? To do so, we conducted a series of forecast sensitivity experiments for three extreme El Niños (1982, 1997 and 2015). The forecast simulations were initialised on 1 September using the present ocean mean state and using a hypothetical future ocean mean state that was created by addition of the doubled observed ocean temperature trends over the past 60 years (we have referred to this as the warmer ocean mean state). Present and warmer climatological forecasts were also produced with the climatological conditions of 1 September of 1981–2013 without and with the addition of the doubled observed ocean temperature trends, respectively, as reference forecast sets.

Statistical analysis of observations and the forecast experiments confirmed that extreme El Niño is a key driver of strong low SAM in the October to December season. The observed El Niño-low SAM relationship is robust once we remove the influence of the Antarctic stratospheric vortex on the SAM by regressing it out in the observational analysis. This conclusion is also supported by the forecast experiments that use scrambled atmospheric initial conditions so that the promotion of low SAM during the El Niño forecasts can be attributed to the El Niño forcing.

The El Niño experiments with the present and the hypothetical warmer ocean mean states revealed that extreme El Niño is likely to lose some strength particularly over the tropical eastern Pacific as a result of reduced mean upper ocean upwelling on the warmer ocean mean state which causes a reduction in the thermocline feedback. In our experiments, the increased zonal temperature gradient of the observed La Niña-like upper ocean mean state change contributes to anomalous SST warming associated with extreme El Niño, as highlighted by the recent study of Wang et al.68. However, this enhanced SST growth appears to be limited to the central and western Pacific65 and is offset by weaker anomaly growth caused by the reduced strength of mean westward surface currents in the same region. The air-sea coupling strength appears to also significantly weaken in the eastern Pacific due to the cooler equatorial mean SST and weaker SST anomalies associated El Niño and shift westward in the warmer climate, further contributing to the weakening of extreme El Niño in the eastern Pacific. This delicate balance of processes in our model framework tips the scale in favor of weaker extreme El Niño events in a warmer world with the continuation of the observed long-term ocean temperature trends.

The weakened El Niño in the warmer climate appears to result in a significant weakening of the low SAM response. Our present climate El Niño experiment confirmed the findings of earlier studies10,11,12 that tropical upper tropospheric warming caused by El Niño increases westerlies on the equatorward side of the SH subtropical jet, shifting the critical latitude equatorward. This shift induces an equatorward shift of the momentum flux convergence-divergence dipole, resulting in low SAM. This chain of processes appears to continue to operate in the warmer climate but with significantly weaker strength of zonal wind, eddy momentum transport and pressure anomalies because El Niño itself becomes weaker and the tropical convective heating response to El Niño becomes weaker as well. Therefore, if the ocean temperature trends that have been observed over the past 60 years continue into the future, predictability of low SAM and its associated impacts on SH surface climate during extreme El Niño is likely to be substantially reduced. A previous study60 indicated that La Niña may also be expected to weaken in response to a “La Niña-like” temperature trend, but a further study is required to assess whether there will be a similar reduction in predictability of high SAM. An idealised experiment like ours but to examine the ENSO-SAM relationship with El Nino-like warming trends in the tropical Pacific, as suggested by most CMIP models, would also be valuable in understanding the mechanisms of the resultant changes in ENSO and SAM.

## Methods

### Data for observational analyses

We used reanalysis data of the European Centre for Medium-Range Weather Forecasts Interim project (ERA-Interim)69 and SST analyses of Hurrell et al.70 and Reynolds OI v271 for the period of 1979–2016. Anomalies were computed against the base period of 1981–2013 for comparisons with POAMA experimental forecasts. The 55 year observed temperature and upper ocean circulation trends were calculated with POAMA ocean data assimilation system reanalysis data (PEODAS)72 for the period 1960–2014.

### Climate indices

The SAM index was obtained by following Gong and Wang’s definition73, which is the difference between normalised zonally averaged MSLP anomalies at 40°S and 65°S. The strength of extreme El Niños was determined by the Niño3.4 index, which was obtained by averaging SSTs over the domain of 5°S-5°N, 190–240°E. The stratosphere-troposphere (S-T) coupled mode index was obtained following the method of Lim et al.18 by applying height-time domain EOF analysis to anomalies (seasonal cycle removed) of monthly mean zonal-mean zonal winds averaged over 55–65°S. The input data to the EOF was ordered from April to March each year for pressure level data extending from 1000 to 1 hPa. The resultant 1st principal component time series consists of one value each year and depicts the year-to-year variations of the SH spring polar vortex strength and its downward coupling. To obtain the de-trended residual SAM independent of the influence of the Antarctic stratospheric polar vortex, we removed from the SAM index the components linearly related to time and the S-T coupled mode index variability.

### Model and initial conditions for forecast sensitivity experiments

For the model experiments, we used the Bureau of Meteorology’s atmosphere-ocean fully coupled dynamical seasonal climate forecast system, POAMA version a34. Its atmosphere and ocean component models are the Bureau’s Atmosphere Model version 3 (T47/L17)74, and the Australian Community Ocean Model version 2 (2° longitude by 0.5–1° latitude from the tropics to the pole)75, which are coupled by OASIS76 coupler.

High quality observed atmosphere, land and ocean conditions are generated from the Bureau of Meteorology’s atmosphere and land initialisation scheme (ALI)77 and PEODAS, respectively.

### Statistical significance tests

Statistical significance of the observed trends in Figs 1 and 6 was assessed by a two-sided Student t-test with 55 samples (i.e. data of 1960–2014), using the incomplete beta function available in NCAR Command Language. In Fig. 4, 1 standard deviation threshold was used to indicate if air temperature, zonal wind and eddy momentum flux convergence anomalies associated with extreme El Niño are significantly different from the climatological conditions in the ERA-Interim reanalysis set. A two-sided Student t-test was used to estimate the statistical significance on the difference of the two means (Figs 8, 9a,b,d,e, 10) and of the two regression coefficients (Fig. 9c,f) with 99 samples for p’ElNiño vs. w’ElNiño and and the difference of the two means with 33 samples for pClim vs wClim (Figs 6b,d,f and 7).

### Ocean mixed layer heat advection analysis

Based on the ocean mixed layer heat budget and heat advection analyses of Zhao et al.60 and Abellan et al.78, we computed contributions of advective feedback terms to the change of El Niño growth as follows:

$$\frac{\partial T}{\partial t}\approx -\,{\boldsymbol{u}}\nabla T$$
(1)
$${\boldsymbol{u}}\nabla T=\bar{u}T{^{\prime} }_{x}+u^{\prime} {\bar{T}}_{x}+u^{\prime} T{^{\prime} }_{x}+\bar{v}T{^{\prime} }_{y}+v^{\prime} {\bar{T}}_{y}+v^{\prime} T{^{\prime} }_{y}+\bar{w}T{^{\prime} }_{z}+w^{\prime} {\bar{T}}_{z}+w^{\prime} T{^{\prime} }_{z}$$
(2)

where u denotes zonal, meridional and vertical velocities, and Tx, Ty and Tz indicate zonal (x), meridional (y) and vertical (z) temperature gradients, respectively. Overbar and prime signs denote a temporal mean and a departure from the mean, respectively.

From Eq. (2), we further limited our interest to linear terms to tease out the contributions from the mean state changes vs the anomalous changes to changes in the growth of El Niño. To estimate the contributions of the mean state changes, we used the mean state terms $$\bar{u}$$, $$\bar{v}$$, $$\bar{w}\,$$, $${\bar{T}}_{x}$$, $${\bar{T}}_{y}$$, and $${\bar{T}}_{z}$$ from the two different climatologies, pClim and wClim, together with the anomalies from the present climate u′, v′, w′, Tx, Ty and Tz. Then, to estimate the contributions of the anomalous changes in El Niño, we used the anomalies u′, v′, w′, Tx, Ty and Tz of El Niño for the two different mean states, wElNiño’ and pElNiño’, together with the mean state $$u^{\prime}$$, $$v^{\prime}$$, $$w^{\prime}$$, $${T^{\prime} }_{x}$$, $${T^{\prime} }_{y}$$, and $${T^{\prime} }_{z}$$ from the present climatology (pClim).

Figure 8 and Supplementary Figs. S5–7 show that the most important term that explains the reduction of El Niño strength in the tropical eastern Pacific is the thermocline feedback change by the reduced mean upwelling on the warmer ocean, $$\bar{w}{{T}^{^{\prime} }}_{z}$$, in Eq. (2). The reduced mean westward currents of the warmer mean state in $$\bar{u}T{\text{'}}_{x}\,\,$$and the reduced anomalous eastward current of warmer climate El Niño in $$u^{\prime} {\bar{T}}_{x}\,\,$$also suggest some contributions to the reduction of SST warming associated with warmer climate El Niño over the tropical central Pacific. However, as the enhanced zonal temperature gradient of the warmer mean state in $$u^{\prime} {\bar{T}}_{x}\,\,$$somewhat compensates the temperature reduction, an overall contribution from the zonal advective feedback change to the change in the SSTs of warmer climate El Niño seems small.

## Change history

• ### 04 February 2020

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

## References

1. 1.

Kidson, J. W. Indices of the Southern Hemisphere Zonal Wind. Journal of Climate 1, 183–194 (1988).

2. 2.

Hartmann, D. L. & Lo, F. Wave-Driven Zonal Flow Vacillation in the Southern Hemisphere. J. Atmos. Sci. 55, 1303–1315 (1998).

3. 3.

Thompson, D. W. J. & Wallace, J. M. Annular Mode in the extratropical circulation. Part I: Month-to-month variability. J. Clim. 13, 1000–1016 (2000).

4. 4.

Lim, E.-P., Hendon, H. H. & Rashid, H. Seasonal predictability of the Southern Annular Mode due to its association with ENSO. J. Clim. 26, 8037–8054 (2013).

5. 5.

Karoly, D. J., Hope, P. & Jones, P. D. Decadal Variations of the Southern Hemisphere Circulation. Int. J. Climatol. 16, 723–738 (1996).

6. 6.

Lorenz, D. J. & Hartmann, D. L. Eddy – Zonal Flow Feedback in the Southern Hemisphere. J. Atmos. Sci. 3312–332, 10.11757 (2001).

7. 7.

Limpasuvan, V. & Hartmann, D. L. Eddies and the annular modes of climate variability. Geophys. Res. Lett. 26, 3133–3136 (1999).

8. 8.

Chen, G., Lu, J. & Frierson, D. M. W. Phase Speed Spectra and the Latitude of Surface Westerlies: Interannual Variability. J. Clim. 5942–5959, https://doi.org/10.1175/2008JCLI2306.1 (2008).

9. 9.

Zhou, T. J. & Yu, R. C. Sea-surface temperature induced variability of the Southern Annular Mode in an atmospheric general circulation model. Geophys Res Lett 31, L24206 (2004).

10. 10.

L’Heureux, M. L. & Thompson, D. W. J. Observed relationships between the El Niño–Southern Oscillation and the extratropical zonal-mean circulation. J. Clim. 19, 276–287 (2006).

11. 11.

Seager, R., Harnik, N., Kushnir, Y., Robinson, W. A. & Miller, J. A. Mechanisms of hemispherically symetric climate variability. J. Clim. 16, 2960–2978 (2003).

12. 12.

Lu, J., Chen, G. & Frierson, D. M. W. Response of the zonal mean atmospheric circulation to El Niño versus global warming. J. Clim. 21, 5835–5851 (2008).

13. 13.

Silvestri, G. E. & Vera, C. S. Antarctic Oscillation signal on precipitation anomalies over southeastern South America. Geophys. Res. Lett. 30, 2115 (2003).

14. 14.

Adames, Á. F. & Wallace, J. M. On the Tropical Atmospheric Signature of El Niño. J. Atmos. Sci. 74, 1923–1939 (2017).

15. 15.

Ding, Q., Steig, E. J., Battisti, D. S. & Wallace, J. M. Influence of the tropics on the southern annular mode. J. Clim. 25, 6330–6348 (2012).

16. 16.

Seviour, W. J. M. et al. Skillful seasonal prediction of the Southern Annular Mode and Antarctic ozone. J. Clim. 27, 7462–7474 (2014).

17. 17.

Byrne, N. J. & Shepherd, T. G. Seasonal persistence of circulation anomalies in the Southern Hemisphere stratosphere and its implications for the troposphere. J. Clim. 31, 3467–3483 (2018).

18. 18.

Lim, E.-P., Hendon, H. H. & Thompson, D. W. J. Seasonal evolution of stratosphere-troposphere coupling in the Southern Hemisphere and implications for the predictability of surface climate. J. Geophys. Res. Atmos. 123(12), 002–12,016 (2018).

19. 19.

Barnston, A. G., Tippett, M. K., L’Heureux, M. L., Li, S. & Dewitt, D. G. Skill of real-time seasonal ENSO model predictions during 2002–11: Is our capability increasing? Bull. Am. Meteorol. Soc. 93, 631–651 (2012).

20. 20.

Reason, C. J. C. & Rouault, M. Links between the Antarctic Oscillation and winter rainfall over western South Africa. Geophys. Res. Lett. 32, 1–4 (2005).

21. 21.

Gillett, N. P., Kell, T. D. & Jones, P. D. Regional climate impacts of the Southern Annular Mode. Geophys. Res. Lett. 33, L23704 (2006).

22. 22.

Sen Gupta, A. & England, M. H. Coupled ocean-atmosphere-ice response to variations in the southern annular mode. J. Clim. 19, 4457–4486 (2006).

23. 23.

Ciasto, L. M., Alexander, M. A., Deser, C. & England, M. H. On the persistence of cold-season SST anomalies associated with the annular modes. J. Clim. 24, 2500–2515 (2011).

24. 24.

Lim, E. P. & Hendon, H. H. Understanding and predicting the strong Southern Annular Mode and its impact on the record wet east Australian spring 2010. Clim. Dyn. 44, 2807–2824 (2015).

25. 25.

Wang, G. et al. Compounding tropical and stratospheric forcing of the record low Antarctic sea-ice in 2016. Nat. Commun. 10, 13 (2019).

26. 26.

Troccoli, A., Harrison, M., Anderson, D. L. T. & Mason, S. J. Seasonal Climate: Forecasting and Managing Risk [electronic resource]/edited by Alberto Troccoli, Mike Harrison, David L. T. Anderson, Simon J. Mason. NATO Science Series, Vol. 82 (2008).

27. 27.

Power, S., Delage, F., Chung, C., Kociuba, G. & Keay, K. Robust twenty-first-century projections of El Niño and related precipitation variability. Nature 502, 541–545 (2013).

28. 28.

Cai, W. et al. Increased variability of eastern Pacific El Niño under greenhouse warming. Nature 564, 201–206 (2018).

29. 29.

Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An Overview of CMIP5 and the Experiment Design. Bull. Am. Meteorol. Soc. 93, 485–498 (2012).

30. 30.

Lim, E.-P. et al. The impact of the Southern Annular Mode on future changes in Southern Hemisphere rainfall. Geophys. Res. Lett. 43, 7160–7167 (2016).

31. 31.

Dommenget, D., Bayr, T. & Frauen, C. Analysis of the non-linearity in the pattern and time evolution of El Niño southern oscillation. Clim. Dyn. 40, 2825–2847 (2013).

32. 32.

Chung, C. T. Y., Power, S. B., Arblaster, J. M., Rashid, H. A. & Roff, G. L. Nonlinear precipitation response to El Niño and global warming in the Indo-Pacific. Clim. Dyn. 42, 1837–1856 (2014).

33. 33.

Chung, C. T. Y. & Power, S. B. Precipitation response to La Niña and global warming in the Indo-Pacific. Clim. Dyn., https://doi.org/10.1007/s00382-014-2105-9 (2014).

34. 34.

Cottrill, A. et al. Seasonal Forecasting in the Pacific Using the Coupled Model POAMA-2. Weather Forecast. 28, 668–680 (2013).

35. 35.

Deser, C., Phillips, A. S. & Alexander, M. A. Twentieth century tropical sea surface temperature trends revisited. Geophys. Res. Lett. 37, 1–6 (2010).

36. 36.

Chen, C., Cane, M. A., Wittenberg, A. T. & Chen, D. ENSO in the CMIP5 Simulations: Life Cycles, Diversity, and Responses to Climate Change. J. Clim. 30, 775–801 (2017).

37. 37.

Ashok, K., Sabin, T. P., Swapna, P. & Murtugudde, R. G. Is a global warming signature emerging in the tropical Pacific? Geophys. Res. Lett. 39, 1–5 (2012).

38. 38.

Seager, R. et al. Strengthening tropical Pacific zonal sea surface temperature gradient consistent with rising greenhouse gases. Nat. Clim. Chang. 9, 517–522 (2019).

39. 39.

Jiang, N. & Zhu, C. Asymmetric Changes of ENSO Diversity Modulated by the Cold Tongue Mode Under Recent Global Warming. Geophys. Res. Lett. 45, 12,506–12,513 (2018).

40. 40.

Cane, M. A. et al. Twentieth-century sea surface temperature trends. Science (80-.). 275, 957–960 (1997).

41. 41.

Meehl, G. A. et al. Decadal prediction: Can it be skillful? Bull. Am. Meteorol. Soc. 90, 1467–1485 (2009).

42. 42.

Zheng, X. T., Xie, S. P., Lv, L. H. & Zhou, Z. Q. Intermodel uncertainty in ENSO amplitude change tied to Pacific Ocean warming pattern. J. Clim. 29, 7265–7279 (2016).

43. 43.

Santoso, A. et al. Dynamics and Predictability of El Niño–Southern Oscillation: An Australian Perspective on Progress and Challenges. Bull. Am. Meteorol. Soc. 100, 403–420 (2019).

44. 44.

Boer, G. J. Climate change and the regulation of the surface moisture and energy budgets. Clim. Dyn. 8, 225–239 (1993).

45. 45.

Knutson, T. R. & Manabe, S. Time-Mean Response over the Tropical Pacific to Increased CO2 in a Coupled Ocean-Atmosphere Model. J. Clim. 8, 2181–2199 (1995).

46. 46.

Held, I. M. & Soden, B. J. Robust responses of the hydrological cycle to global warming. J. Clim. 19, 5686–5699 (2006).

47. 47.

Collins, M. et al. The impact of global warming on the tropical Pacific Ocean and El Niño. Nat. Geosci. 3, 391 (2010).

48. 48.

Luo, J. J., Wang, G. & Dommenget, D. May common model biases reduce CMIP5’s ability to simulate the recent Pacific La Niña-like cooling? Clim. Dyn. 50, 1335–1351 (2018).

49. 49.

Cai, W. et al. Increasing frequency of extreme El Niño events due to greenhouse warming. Nat. Clim. Chang. 4, 111–116 (2014).

50. 50.

Dinezio, P. N. et al. Mean climate controls on the simulated response of ENSO to increasing greenhouse gases. J. Clim. 25, 7399–7420 (2012).

51. 51.

Kohyama, T., Hartmann, D. L. & Battisti, D. S. La Niña–like Mean-State Response to Global Warming and Potential Oceanic Roles. J. Clim. 30, 4207–4225 (2017).

52. 52.

Kohyama, T., Hartmann, D. L. & Battisti, D. S. Weakening of Nonlinear ENSO Under Global Warming. Geophys. Res. Lett. 45, 8557–8567 (2018).

53. 53.

Cai, W. et al. Pantropical climate interactions. Science (80-.). 363 (2019).

54. 54.

Clement, A. C., Seager, R., Cane, M. A. & Zebiak, S. E. An Ocean Dynamical Thermostat. J. Clim. 9, 2190–2196 (1996).

55. 55.

Kohyama, T. & Hartmann, D. L. Nonlinear ENSO warming suppression (NEWS). J. Clim. 30, 4227–4251 (2017).

56. 56.

Luo, J.-J., Sasaki, W. & Masumoto, Y. Indian Ocean warming modulates Pacific climate change. Proc. Natl. Acad. Sci. USA 109, 18701–6 (2012).

57. 57.

McGregor, S. et al. Recent walker circulation strengthening and pacific cooling amplified by atlantic warming. Nat. Clim. Chang. 4, 888–892 (2014).

58. 58.

Solomon, A. & Newman, M. Reconciling disparate twentieth-century Indo-Pacific ocean temperature trends in the instrumental record. Nat. Clim. Chang. 2, 691–699 (2012).

59. 59.

Lim, E. P. et al. Interaction of the recent 50 year SST trend and La Niña 2010: amplification of the Southern Annular Mode and Australian springtime rainfall. Clim. Dyn. 47, 2273–2291 (2016).

60. 60.

Zhao, M., Hendon, H., Oscar, A., Liu, G. & Guomin, W. Weakened Eastern Pacific El Niño Predictability in the Early Twenty-First Century. J. Clim. 29, 6805–6822 (2016).

61. 61.

Thompson, D. W. J. & Solomon, S. Interpretation of recent Southern Hemisphere climate change. Science (80-.). 296, 895–899 (2002).

62. 62.

Hurwitz, M. M., Newman, P. A., Oman, L. D. & Molod, A. M. Response of the Antarctic Stratosphere to Two Types of El Niño Events. J. Atmos. Sci. 68, 812–822 (2011).

63. 63.

Collins, M. et al. Long-term Climate Change: Projections, Commitments and Irreversibility. Clim. Chang. 2013 Phys. Sci. Basis. Contrib. Work. Gr. I to Fifth Assess. Rep. Intergov. Panel Clim. Chang. 1029–1136, https://doi.org/10.1017/CBO9781107415324.024 (2013).

64. 64.

Chung, P. H. & Li, T. Interdecadal relationship between the mean state and El Niño types. J. Clim. 26, 361–379 (2013).

65. 65.

N. C. Johnson, M. L. L’Heureux, C.‐H. & Chang, Z.‐Z. Hu, On the Delayed Coupling Between Ocean and Atmosphere in Recent Weak El Niño Episodes. Geophysical Research Letters https://doi.org/10.1029/2019GL084021 (2019).

66. 66.

McPhaden, M. J., Lee, T. & McClurg, D. El Niño and its relationship to changing background conditions in the tropical Pacific Ocean. Geophys. Res. Lett. 38, 2–5 (2011).

67. 67.

Lübbecke, J. F. & Mcphaden, M. J. Assessing the twenty-first-century shift in enso variability in terms of the bjerknes stability index. J. Clim. 27, 2577–2587 (2014).

68. 68.

Bin W. et al. Historical change of El Niño properties sheds light on future changes of extreme El Niño. Proceedings of the National Academy of Sciences 116(45), 22512–22517 (2019).

69. 69.

Dee, D. et al. The ERA - Interim reanalysis: Configuration and performance of the data assimilation system. Quaterly J. R. Meteorol. Soc. 137, 553–597 (2011).

70. 70.

Hurrell, J. W., Hack, J. J., Shea, D., Caron, J. M. & Rosinski, J. A new sea surface temperature and sea ice boundary dataset for the community atmosphere model. J. Clim. 21, 5145–5153 (2008).

71. 71.

Reynolds, R. W., Rayner, N. A., Smith, T. M., Stokes, D. C. & Wang, W. An improved in situ and satellite SST analysis for climate. J. Clim. 15, 1609–1625 (2002).

72. 72.

Yin, Y., Alves, O. & Oke, P. R. An Ensemble Ocean Data Assimilation System for Seasonal Prediction. Mon. Weather Rev. 139, 786–808 (2011).

73. 73.

Gong, D. & Wang, S. Definition of Antarctic Oscillation index. Geophys. Res. Lett. 26, 459–462 (1999).

74. 74.

Colman, R. et al. BMRC Atmospheric Model (BAM) version 3.0: comparison with mean climatology (2005).

75. 75.

Oke, P. R., Schiller, A., Griffin, D. A. & Brassington, G. B. Ensemble data assimilation for an eddy-resolving ocean model of the Australian region. Q. J. R. Meteorol. Soc. 131, 3301–3311 (2005).

76. 76.

Valke, S., Terray, L. & Piacentini, A. The OASIS coupled user guide version 2.4. Tech. Rep. TR/CMGC/00-10, CERFACS (2000).

77. 77.

Hudson, D., Alves, O., Hendon, H. & Wang, G. The impact of atmospheric initialisation on seasonal prediction of tropical Pacific SST. Clim. Dyn. (2011).

78. 78.

Abellán, E., McGregor, S., England, M. H. & Santoso, A. Distinctive role of ocean advection anomalies in the development of the extreme 2015–16 El Niño. Clim. Dyn. 51, 2191–2208 (2018).

## Acknowledgements

This study was supported in part by the Australian Government’s National Environmental Science Programme. Michael J. McPhaden was supported by NOAA. The authors are grateful to Drs Irina Rudeva and Blair Trewin at the Bureau of Meteorology and two anonymous reviewers for their constructive feedback on the manuscript. This research was undertaken on the NCI National Facility in Canberra, Australia, which is supported by the Australian Commonwealth Government. The NCAR Command Language (NCL; http://www.ncl.ucar.edu) version 6.4.0 was used for data analysis and visualization of the results. We also acknowledge NCAR/UCAR, NOAA and ECMWF for producing and providing Hurrell et al. (2008) SST analysis, Reynolds OI v2 SST analysis, and ERA-Interim reanalysis, respectively. This is PMEL contribution no. 4994.

## Author information

Authors

### Contributions

E.L. and P.H. conceived the idea, and E.L. conducted the experiments, analysed the results and wrote the first draft with the help of C.C. and F.D.; H.H.H. and M.J.M. contributed to the interpretation of the results, and all authors contributed to the writing of the manuscript.

### Corresponding author

Correspondence to Eun-Pa Lim.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Lim, EP., Hendon, H.H., Hope, P. et al. Continuation of tropical Pacific Ocean temperature trend may weaken extreme El Niño and its linkage to the Southern Annular Mode. Sci Rep 9, 17044 (2019). https://doi.org/10.1038/s41598-019-53371-3

• Accepted:

• Published:

• ### Long‐Term Trend of Equatorial Atlantic Zonal Sea Surface Temperature Gradient Linked to the Tropical Pacific Cold Tongue Mode Under Global Warming

• Yang Li
• , Quanliang Chen
• , Nan Xing
• , Zhigang Cheng
• , Yulei Qi
• , Fan Feng
•  & Minggang Li

Journal of Geophysical Research: Oceans (2021)

• ### Exploring atmospheric circulation leading to three anomalous Australian spring heat events

• Roseanna C. McKay
• , Julie M. Arblaster
• , Pandora Hope
•  & Eun-Pa Lim

Climate Dynamics (2021)

• ### Pantropical Response to Global Warming and the Emergence of a La Niña‐Like Mean State Trend

• Sang‐Ki Lee
• , Dongmin Kim
• , Gregory R. Foltz
•  & Hosmay Lopez

Geophysical Research Letters (2020)

• ### Role of Tropical Variability in Driving Decadal Shifts in the Southern Hemisphere Summertime Eddy-Driven Jet

• Dongxia Yang
• , Julie M. Arblaster
• , Gerald A. Meehl
• , Matthew H. England
• , Eun-Pa Lim
• , Susan Bates
•  & Nan Rosenbloom

Journal of Climate (2020)

• ### Emerging arboviruses in the urbanized Amazon rainforest

• Rachel Lowe
• , Sophie Lee
• , Raquel Martins Lana
• , Cláudia Torres Codeço
• , Marcia C Castro
•  & Mercedes Pascual

BMJ (2020)
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 2, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0, "math_score": 0.8254827857017517, "perplexity": 5723.461981236543}, "config": {"markdown_headings": true, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "plain"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "warc_path": "s3://commoncrawl/crawl-data/CC-MAIN-2021-25/segments/1623487636559.57/warc/CC-MAIN-20210618104405-20210618134405-00258.warc.gz"}