Explainable El Niño predictability from climate mode interactions.

The El Niño-Southern Oscillation (ENSO) provides most of the global seasonal climate forecast skill1-3, yet, quantifying the sources of skilful predictions is a long-standing challenge4-7. Different sources of predictability affect ENSO evolution, leading to distinct global effects. Artificial intelligence forecasts offer promising advancements but linking their skill to specific physical processes is not yet possible8-10, limiting our understanding of the dynamics underpinning the advancements. Here we show that an extended nonlinear recharge oscillator (XRO) model shows skilful ENSO forecasts at lead times up to 16-18 months, better than global climate models and comparable to the most skilful artificial intelligence forecasts. The XRO parsimoniously incorporates the core ENSO dynamics and ENSO's seasonally modulated interactions with other modes of variability in the global oceans. The intrinsic enhancement of ENSO's long-range forecast skill is traceable to the initial conditions of other climate modes by means of their memory and interactions with ENSO and is quantifiable in terms of these modes' contributions to ENSO amplitude. Reforecasts using the XRO trained on climate model output show that reduced biases in both model ENSO dynamics and in climate mode interactions can lead to more skilful ENSO forecasts. The XRO framework's holistic treatment of ENSO's global multi-timescale interactions highlights promising targets for improving ENSO simulations and forecasts.

Increased occurrences of consecutive La Niña events under global warming.

Most El Niño events occur sporadically and peak in a single winter1-3, whereas La Niña tends to develop after an El Niño and last for two years or longer4-7. Relative to single-year La Niña, consecutive La Niña features meridionally broader easterly winds and hence a slower heat recharge of the equatorial Pacific6,7, enabling the cold anomalies to persist, exerting prolonged impacts on global climate, ecosystems and agriculture8-13. Future changes to multi-year-long La Niña events remain unknown. Here, using climate models under future greenhouse-gas forcings14, we find an increased frequency of consecutive La Niña ranging from 19 ± 11% in a low-emission scenario to 33 ± 13% in a high-emission scenario, supported by an inter-model consensus stronger in higher-emission scenarios. Under greenhouse warming, a mean-state warming maximum in the subtropical northeastern Pacific enhances the regional thermodynamic response to perturbations, generating anomalous easterlies that are further northward than in the twentieth century in response to El Niño warm anomalies. The sensitivity of the northward-broadened anomaly pattern is further increased by a warming maximum in the equatorial eastern Pacific. The slower heat recharge associated with the northward-broadened easterly anomalies facilitates the cold anomalies of the first-year La Niña to persist into a second-year La Niña. Thus, climate extremes as seen during historical consecutive La Niña episodes probably occur more frequently in the twenty-first century.

Publisher Correction: Butterfly effect and a self-modulating El Niño response to global warming.

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

Butterfly effect and a self-modulating El Niño response to global warming.

El Niño and La Niña, collectively referred to as the El Niño-Southern Oscillation (ENSO), are not only highly consequential1-6 but also strongly nonlinear7-14. For example, the maximum warm anomalies of El Niño, which occur in the equatorial eastern Pacific Ocean, are larger than the maximum cold anomalies of La Niña, which are centred in the equatorial central Pacific Ocean7-9. The associated atmospheric nonlinear thermal damping cools the equatorial Pacific during El Niño but warms it during La Niña15,16. Under greenhouse warming, climate models project an increase in the frequency of strong El Niño and La Niña events, but the change differs vastly across models17, which is partially attributed to internal variability18-23. Here we show that like a butterfly effect24, an infinitesimal random perturbation to identical initial conditions induces vastly different initial ENSO variability, which systematically affects its response to greenhouse warming a century later. In experiments with higher initial variability, a greater cumulative oceanic heat loss from ENSO thermal damping reduces stratification of the upper equatorial Pacific Ocean, leading to a smaller increase in ENSO variability under subsquent greenhouse warming. This self-modulating mechanism operates in two large ensembles generated using two different models, each commencing from identical initial conditions but with a butterfly perturbation24,25; it also operates in a large ensemble generated with another model commencing from different initial conditions25,26 and across climate models participating in the Coupled Model Intercomparison Project27,28. Thus, if the greenhouse-warming-induced increase in ENSO variability29 is initially suppressed by internal variability, future ENSO variability is likely to be enhanced, and vice versa. This self-modulation linking ENSO variability across time presents a different perspective for understanding the dynamics of ENSO variability on multiple timescales in a changing climate.

Twofold expansion of the Indo-Pacific warm pool warps the MJO life cycle.

The Madden-Julian Oscillation (MJO) is the most dominant mode of subseasonal variability in the tropics, characterized by an eastward-moving band of rain clouds. The MJO modulates the El Niño Southern Oscillation1, tropical cyclones2,3 and the monsoons4-10, and contributes to severe weather events over Asia, Australia, Africa, Europe and the Americas. MJO events travel a distance of 12,000-20,000 km across the tropical oceans, covering a region that has been warming during the twentieth and early twenty-first centuries in response to increased anthropogenic emissions of greenhouse gases11, and is projected to warm further. However, the impact of this warming on the MJO life cycle is largely unknown. Here we show that rapid warming over the tropical oceans during 1981-2018 has warped the MJO life cycle, with its residence time decreasing over the Indian Ocean by 3-4 days, and increasing over the Indo-Pacific Maritime Continent by 5-6 days. We find that these changes in the MJO life cycle are associated with a twofold expansion of the Indo-Pacific warm pool, the largest expanse of the warmest ocean temperatures on Earth. The warm pool has been expanding on average by 2.3 × 105 km2 (the size of Washington State) per year during 1900-2018 and at an accelerated average rate of 4 × 105 km2 (the size of California) per year during 1981-2018. The changes in the Indo-Pacific warm pool and the MJO are related to increased rainfall over southeast Asia, northern Australia, Southwest Africa and the Amazon, and drying over the west coast of the United States and Ecuador.

Author Correction: El Niño-Southern Oscillation complexity.

In this Review, the middle initial of author Kim M. Cobb was omitted. The original Review has been corrected online.

Increased variability of eastern Pacific El Niño under greenhouse warming.

The El Niño-Southern Oscillation (ENSO) is the dominant and most consequential climate variation on Earth, and is characterized by warming of equatorial Pacific sea surface temperatures (SSTs) during the El Niño phase and cooling during the La Niña phase. ENSO events tend to have a centre-corresponding to the location of the maximum SST anomaly-in either the central equatorial Pacific (5° S-5° N, 160° E-150° W) or the eastern equatorial Pacific (5° S-5° N, 150°-90° W); these two distinct types of ENSO event are referred to as the CP-ENSO and EP-ENSO regimes, respectively. How the ENSO may change under future greenhouse warming is unknown, owing to a lack of inter-model agreement over the response of SSTs in the eastern equatorial Pacific to such warming. Here we find a robust increase in future EP-ENSO SST variability among CMIP5 climate models that simulate the two distinct ENSO regimes. We show that the EP-ENSO SST anomaly pattern and its centre differ greatly from one model to another, and therefore cannot be well represented by a single SST 'index' at the observed centre. However, although the locations of the anomaly centres differ in each model, we find a robust increase in SST variability at each anomaly centre across the majority of models considered. This increase in variability is largely due to greenhouse-warming-induced intensification of upper-ocean stratification in the equatorial Pacific, which enhances ocean-atmosphere coupling. An increase in SST variance implies an increase in the number of 'strong' EP-El Niño events (corresponding to large SST anomalies) and associated extreme weather events.

El Niño-Southern Oscillation complexity.

El Niño events are characterized by surface warming of the tropical Pacific Ocean and weakening of equatorial trade winds that occur every few years. Such conditions are accompanied by changes in atmospheric and oceanic circulation, affecting global climate, marine and terrestrial ecosystems, fisheries and human activities. The alternation of warm El Niño and cold La Niña conditions, referred to as the El Niño-Southern Oscillation (ENSO), represents the strongest year-to-year fluctuation of the global climate system. Here we provide a synopsis of our current understanding of the spatio-temporal complexity of this important climate mode and its influence on the Earth system.

Seasonal sea surface cooling in the equatorial Pacific cold tongue controlled by ocean mixing.

Sea surface temperature (SST) is a critical control on the atmosphere, and numerical models of atmosphere-ocean circulation emphasize its accurate prediction. Yet many models demonstrate large, systematic biases in simulated SST in the equatorial 'cold tongues' (expansive regions of net heat uptake from the atmosphere) of the Atlantic and Pacific oceans, particularly with regard to a central but little-understood feature of tropical oceans: a strong seasonal cycle. The biases may be related to the inability of models to constrain turbulent mixing realistically, given that turbulent mixing, combined with seasonal variations in atmospheric heating, determines SST. In temperate oceans, the seasonal SST cycle is clearly related to varying solar heating; in the tropics, however, SSTs vary seasonally in the absence of similar variations in solar inputs. Turbulent mixing has long been a likely explanation, but firm, long-term observational evidence has been absent. Here we show the existence of a distinctive seasonal cycle of subsurface cooling via mixing in the equatorial Pacific cold tongue, using multi-year measurements of turbulence in the ocean. In boreal spring, SST rises by 2 kelvin when heating of the upper ocean by the atmosphere exceeds cooling by mixing from below. In boreal summer, SST decreases because cooling from below exceeds heating from above. When the effects of lateral advection are considered, the magnitude of summer cooling via mixing (4 kelvin per month) is equivalent to that required to counter the heating terms. These results provide quantitative assessment of how mixing varies on timescales longer than a few weeks, clearly showing its controlling influence on seasonal cooling of SST in a critical oceanic regime.

Late-twentieth-century emergence of the El Niño propagation asymmetry and future projections.

The El Niño/Southern Oscillation (ENSO) is the Earth's most prominent source of interannual climate variability, exerting profound worldwide effects. Despite decades of research, its behaviour continues to challenge scientists. In the eastern equatorial Pacific Ocean, the anomalously cool sea surface temperatures (SSTs) found during La Niña events and the warm waters of modest El Niño events both propagate westwards, as in the seasonal cycle. In contrast, SST anomalies propagate eastwards during extreme El Niño events, prominently in the post-1976 period, spurring unusual weather events worldwide with costly consequences. The cause of this propagation asymmetry is currently unknown. Here we trace the cause of the asymmetry to the variations in upper ocean currents in the equatorial Pacific, whereby the westward-flowing currents are enhanced during La Niña events but reversed during extreme El Niño events. Our results highlight that propagation asymmetry is favoured when the westward mean equatorial currents weaken, as is projected to be the case under global warming. By analysing past and future climate simulations of an ensemble of models with more realistic propagation, we find a doubling in the occurrences of El Niño events that feature prominent eastward propagation characteristics in a warmer world. Our analysis thus suggests that more frequent emergence of propagation asymmetry will be an indication of the Earth's warming climate.

More extreme swings of the South Pacific convergence zone due to greenhouse warming.

The South Pacific convergence zone (SPCZ) is the Southern Hemisphere's most expansive and persistent rain band, extending from the equatorial western Pacific Ocean southeastward towards French Polynesia. Owing to its strong rainfall gradient, a small displacement in the position of the SPCZ causes drastic changes to hydroclimatic conditions and the frequency of extreme weather events--such as droughts, floods and tropical cyclones--experienced by vulnerable island countries in the region. The SPCZ position varies from its climatological mean location with the El Niño/Southern Oscillation (ENSO), moving a few degrees northward during moderate El Niño events and southward during La Niña events. During strong El Niño events, however, the SPCZ undergoes an extreme swing--by up to ten degrees of latitude toward the Equator--and collapses to a more zonally oriented structure with commensurately severe weather impacts. Understanding changes in the characteristics of the SPCZ in a changing climate is therefore of broad scientific and socioeconomic interest. Here we present climate modelling evidence for a near doubling in the occurrences of zonal SPCZ events between the periods 1891-1990 and 1991-2090 in response to greenhouse warming, even in the absence of a consensus on how ENSO will change. We estimate the increase in zonal SPCZ events from an aggregation of the climate models in the Coupled Model Intercomparison Project phases 3 and 5 (CMIP3 and CMIP5) multi-model database that are able to simulate such events. The change is caused by a projected enhanced equatorial warming in the Pacific and may lead to more frequent occurrences of extreme events across the Pacific island nations most affected by zonal SPCZ events.

Central-Pacific El Niño-Southern Oscillation less predictable under greenhouse warming.

El Niño-Southern Oscillation (ENSO) is the dominant mode of interannual climate variability in the tropical Pacific, whose nature nevertheless may change significantly in a warming climate. Here, we show that the predictability of ENSO may decrease in the future. Across the models in the Coupled Model Intercomparison Project Phase 6 (CMIP6), we find a robust decrease of the persistence and predictability for the Central Pacific (CP) ENSO under global warming, notably in passing through the boreal spring. The strength of spring predictability barrier will be increased by 25% in the future. The reduced predictability of CP ENSO is caused by the faster warming over surface ocean in tropical Pacific and, in turn, the enhanced thermodynamical damping rate on CP ENSO in response to global warming. In contrast, the predictability of Eastern Pacific ENSO will not change. Our results suggest that future greenhouse warming will make the prediction of CP ENSO more challenging, with far-reaching implications on future climate predictions.

Enhanced risk of record-breaking regional temperatures during the 2023-24 El Niño.

In 2023, the development of El Niño is poised to drive a global upsurge in surface air temperatures (SAT), potentially resulting in unprecedented warming worldwide. Nevertheless, the regional patterns of SAT anomalies remain diverse, obscuring where historical warming records may be surpassed in the forthcoming year. Our study underscores the significant influence of El Niño and the persistence of climate signals on the inter-annual variability of regional SAT, both in amplitude and spatial distribution. The likelihood of global mean SAT exceeding historical records, calculated from July 2023 to June 2024, is estimated at 90%, contingent upon annual-mean sea surface temperature anomalies in the eastern equatorial Pacific exceeding 0.6 °C. Regions particularly susceptible to recording record-high SAT include coastal and adjacent areas in Asia such as the Bay of Bengal and the South China Sea, as well as Alaska, the Caribbean Sea, and the Amazon. This impending warmth heightens the risk of year-round marine heatwaves and escalates the threat of wildfires and other negative consequences in Alaska and the Amazon basin, necessitating strategic mitigation measures to minimize potential worst-case impacts.

Emergence of the Central Atlantic Niño.

The Atlantic Niño is characterized by sea surface warming in the equatorial Atlantic, which can trigger La Niña, the cold phase of El Niño-Southern Oscillation (ENSO). Although observations show that the Atlantic Niño has weakened by approximately 30% since the 1970s, its remote influence on ENSO remains strong. Here, we show that this apparent discrepancy is due to the existence of two types of Atlantic Niño with distinct patterns and climatic impacts, which we refer to as the central and eastern Atlantic Niño. Our results show that with equal strength, the central Atlantic Niño has a stronger influence on tropical climate than its eastern counterpart. Meanwhile, the eastern Atlantic Niño has weakened by approximately 50% in recent decades, allowing the central Atlantic Niño to emerge and dominate the remote impact on ENSO. Given the distinct climatic impacts of the two types, it is necessary to distinguish between them and investigate their behaviors and influences on climate in future studies.

Controls on tropical Pacific Ocean productivity revealed through nutrient stress diagnostics.

In situ enrichment experiments have shown that the growth of bloom-forming diatoms in the major high-nitrate low-chlorophyll (HNLC) regions of the world's oceans is limited by the availability of iron. Yet even the largest of these manipulative experiments represents only a small fraction of an ocean basin, and the responses observed are strongly influenced by the proliferation of rare species rather than the growth of naturally dominant populations. Here we link unique fluorescence attributes of phytoplankton to specific physiological responses to nutrient stress, and use these relationships to evaluate the factors that constrain phytoplankton growth in the tropical Pacific Ocean on an unprecedented spatial scale. On the basis of fluorescence measurements taken over 12 years, we delineate three major ecophysiological regimes in this region. We find that iron has a key function in regulating phytoplankton growth in both HNLC and oligotrophic waters near the Equator and further south, whereas nitrogen and zooplankton grazing are the primary factors that regulate biomass production in the north. Application of our findings to the interpretation of satellite chlorophyll fields shows that productivity in the tropical Pacific basin may be 1.2-2.5 Pg C yr(-1) lower than previous estimates have suggested, a difference that is comparable to the global change in ocean production that accompanied the largest El Niño to La Niña transition on record.

Emergence of changing Central-Pacific and Eastern-Pacific El Niño-Southern Oscillation in a warming climate.

El Niño-Southern Oscillation (ENSO) features strong warm events in the eastern equatorial Pacific (EP), or mild warm and strong cold events in the central Pacific (CP), with distinct impacts on global climates. Under transient greenhouse warming, models project increased sea surface temperature (SST) variability of both ENSO regimes, but the timing of emergence out of internal variability remains unknown for either regime. Here we find increased EP-ENSO SST variability emerging by around 2030 ± 6, more than a decade earlier than that of CP-ENSO, and approximately four decades earlier than that previously suggested without separating the two regimes. The earlier EP-ENSO emergence results from a stronger increase in EP-ENSO rainfall response, which boosts the signal of increased SST variability, and is enhanced by ENSO non-linear atmospheric feedback. Thus, increased ENSO SST variability under greenhouse warming is likely to emerge first in the eastern than central Pacific, and decades earlier than previously anticipated.

Sea level extremes and compounding marine heatwaves in coastal Indonesia.

Low-lying island nations like Indonesia are vulnerable to sea level Height EXtremes (HEXs). When compounded by marine heatwaves, HEXs have larger ecological and societal impact. Here we combine observations with model simulations, to investigate the HEXs and Compound Height-Heat Extremes (CHHEXs) along the Indian Ocean coast of Indonesia in recent decades. We find that anthropogenic sea level rise combined with decadal climate variability causes increased occurrence of HEXs during 2010-2017. Both HEXs and CHHEXs are driven by equatorial westerly and longshore northwesterly wind anomalies. For most HEXs, which occur during December-March, downwelling favorable northwest monsoon winds are enhanced but enhanced vertical mixing limits surface warming. For most CHHEXs, wind anomalies associated with a negative Indian Ocean Dipole (IOD) and co-occurring La Niña weaken the southeasterlies and cooling from coastal upwelling during May-June and November-December. Our findings emphasize the important interplay between anthropogenic warming and climate variability in affecting regional extremes.

Decadal climate variability in the tropical Pacific: Characteristics, causes, predictability, and prospects.

Climate variability in the tropical Pacific affects global climate on a wide range of time scales. On interannual time scales, the tropical Pacific is home to the El Niño­Southern Oscillation (ENSO). Decadal variations and changes in the tropical Pacific, referred to here collectively as tropical Pacific decadal variability (TPDV), also profoundly affect the climate system. Here, we use TPDV to refer to any form of decadal climate variability or change that occurs in the atmosphere, the ocean, and over land within the tropical Pacific. "Decadal," which we use in a broad sense to encompass multiyear through multidecadal time scales, includes variability about the mean state on decadal time scales, externally forced mean-state changes that unfold on decadal time scales, and decadal variations in the behavior of higher-frequency modes like ENSO.

On the interchangeability of sea-surface and near-surface air temperature anomalies in climatologies.

On global and hemispheric scales, sea-surface temperature (SST) anomalies are assumed to be good surrogates for near-surface marine air temperature (MAT) anomalies. In fact, global gridded temperature datasets commonly blend SST and near-surface air temperature anomalies to overcome the lack of geographically homogeneous and reliable MAT observations. Here, we show that SST and MAT anomalies differ regarding crucial statistical properties such as multiannual trends and probabilistic distributions of daily and monthly averages. We provide evidence of the lack of interchangeability from an array of moored buoys in the tropical Pacific Ocean. We identify statistically significant discrepancies between SST and MAT anomalies for single as well as groups of such buoys. Thus, caution is required when characterizing and interpreting MAT variability through SST observations, especially at shorter than decadal timescale.