Zonal Sea Surface Temperature Gradient Research Articles

Zonal control on Holocene precipitation in northwestern Madagascar based on a stalagmite from Anjohibe

The Malagasy Summer Monsoon is an important part of the larger Indian Ocean and tropical monsoon region. As the effects of global warming play out, changes to precipitation in Madagascar will have important ramifications for the Malagasy people. To help understand how precipitation responds to climate changes we present a long-term Holocene speleothem record from Anjohibe, part of the Andranoboka cave system in northwestern Madagascar. To date, it is the most complete Holocene record from this region and sheds light on the nature of millennial and centennial precipitation changes in this region. We find that over the Holocene, precipitation in northwestern Madagascar is actually in phase with the Northern Hemisphere Asian monsoon on multi-millennial scales, but that during some shorter centennial-scale events such as the 8.2 ka event, Anjohibe exhibits an antiphase precipitation signal to the Northern Hemisphere. The ultimate driver of precipitation changes across the Holocene does not appear to be the meridional migration of the monsoon. Instead, zonal sea surface temperature gradients in the Indian Ocean seem to play a primary role in precipitation changes in northwestern Madagascar.

Emergent Constraint on Projection of the North Pacific Monsoon Trough and Its Implications for Typhoon Activity Using CMIP6 Models

AbstractThe North Pacific monsoon trough (NPMT) is an imperative large‐scale circulation pattern influencing tropical cyclone (TC) activity in the western North Pacific (WNP), thus its future change has great implications for the WNP TC activity. Here future change in the NPMT and its uncertainty are examined by using 35 climate models from Phase six of Coupled Model Intercomparison Project (CMIP6). Due to the El Niño‐like sea surface temperature (SST) pattern, the multi‐model ensemble (MME) mean projects an eastward extension of the NPMT in the latter half of the 21st century. However, considerable inter‐model uncertainty exists in the projection, which is associated with the diversity in the zonal SST gradient across the tropical Pacific. This diversity in the projected SST is largely rooted in the simulated precipitation biases in the tropical Pacific in the historical experiments of CMIP6 models, which can modulate future Pacific SST distribution through the shortwave‐SST feedback. This linkage between biases in historical climate and future climate allows us to constrain the projection by using observational data sets. Emergent constraint using observational precipitation data set reduces the projection uncertainty by 48% and projects an eastward extension of the north Pacific NPMT by 6.2°. This eastward extension of the NPMT indicates an eastward migration of TC genesis location, and elongated TC track and thus strengthened TC intensity, suggesting an increased TC‐related disaster for residents near the WNP.

Uncertainty in the Past and Future Changes of Tropical Pacific SST Zonal Gradient: Internal Variability versus Model Spread

Abstract The zonal sea surface temperature (SST) gradient across the tropical Pacific is a pacemaker of the variable rates of global warming. In both historical simulations and future projections, the current state-of-the-art climate models show evident spread in the changes of zonal SST gradient, but the reasons remain unknown. Here, we quantify the contributions of internal variability and model spread to the uncertainty of zonal SST gradient changes by analyzing 342 realizations from CMIP5 and CMIP6 models and several sets of large-ensemble simulations. We found that the internal variability dominates the total uncertainty at multidecadal time scales (∼31-yr trends). Although the ratio of internal uncertainty to the total uncertainty declines along with higher emission of greenhouse gases under global warming, it is still over 80% at the multidecadal time scales in the future. The Pacific decadal oscillation is identified as the key internal mode responsible for the multidecadal uncertainty. For the future projections at centurial time scales, the uncertainty of zonal SST gradient changes is mainly from the intermodel spread in response to external forcing, accounting for about 70% of the uncertainty based on the difference between 2070–99 and 1970–99. The model spread in the cloud–shortwave radiation–SST feedback over the tropical Pacific is important in the uncertainty of zonal SST gradient changes. In particular, the intensity of negative convective cloud feedback in the western Pacific dominates the spread in CMIP5 models, while the intensity of stratocumulus cloud feedback over the southeastern Pacific is the primary process influencing the uncertainty in CMIP6 models.

Responses of the Madden–Julian Oscillation to Global Warming: Impacts from Tropical Sea Surface Temperature Changes

Abstract The impacts of tropical sea surface temperature (SST) changes on the Madden–Julian oscillation (MJO) are investigated using the large-ensemble simulation from the Community Earth System Model version 2 (CESM2-LE) under the shared socioeconomic pathway (SSP370 scenario). Three SST change patterns are featured, distinguished by the zonal gradient of the change in the equatorial Pacific warming. MJO characteristics and its teleconnections responses are composited for the clusters, and their relationships to the zonal SST gradient changes are examined. Results show that the anomalously strong El Niño–like SST change pattern significantly intensifies the MJO amplitude and enhances its eastward extension compared to the anomalously weak El Niño–like SST change pattern. These changes in MJO amplitude are further interpreted through the α framework. We also found no statistically different extratropical geopotential height responses to MJO between the three SST warming patterns, possibly due to strong internal climate variability. Changes in Rossby wave source between clusters also show a weak relationship with the MJO teleconnections. Our results highlight the importance of Indo-Pacific zonal SST gradient changes on the changes of MJO but limited impacts on MJO teleconnections to the midlatitudes.

High‐Latitude, Indian Ocean, and Orbital Influences on Eastern African Hydroclimate Across the Plio‐Pleistocene Boundary

AbstractTerrestrial‐marine dust fluxes, pedogenic carbonate δ13C values, and various paleovegetation proxies suggest that Africa experienced gradual cooling and drying across the Pliocene‐Pleistocene (Plio‐Pleistocene) boundary (2.58 million years ago [Ma]). However, the timing, magnitude, resolution, and relative influences of orbitally‐driven changes in high latitude glaciations and low latitude insolation differ by region and proxy. To disentangle these forcings and investigate equatorial eastern African climate across the Plio‐Pleistocene boundary, we generated a high‐resolution (∼3,000‐year) data set of compound‐specific n‐alkane leaf wax δ2H values—a robust proxy for atmospheric circulation and precipitation amount—from the HSPDP‐BTB13‐1A core, which spans a ∼3.3–2.6 Ma sequence in the Baringo‐Tugen Hills‐Barsemoi Basin of central Kenya. In combination with the physical sedimentology, our data indicate that precipitation varied strongly with orbital obliquity, not precession, during the late Pliocene, perhaps imparted by variations in the cross‐equatorial insolation gradient. We also observe a marked shift toward wetter conditions beginning ∼3 Ma that corresponds with global cooling, drying in western Australia, and a steepening of the west‐east zonal Indian Ocean (IO) sea surface temperature (SST) gradient. We propose that northward migration of the Subtropical Front reduced Agulhas current leakage, warming the western IO and causing changes in the IO zonal SST gradient at 3 Ma, a process that has been observed in the latest Pleistocene‐Holocene but not over longer timescales. Thus, the late Cenozoic moisture history of eastern Africa is driven by a complex mixture of low‐latitude insolation, the IO SST gradient, and teleconnections to distal high‐latitude cooling.

Zonal Indian Ocean Variability Drives Millennial‐Scale Precipitation Changes in Northern Madagascar

AbstractThe low latitude Indian Ocean is warming faster than other tropical basins, and its interannual climate variability is projected to become more extreme under future emissions scenarios with substantial impacts on developing Indian Ocean rim countries. Therefore, it has become increasingly important to understand the drivers of regional precipitation in a changing climate. Here we present a new speleothem record from Anjohibe, a cave in northwest (NW) Madagascar well situated to record past changes in the Intertropical Convergence Zone (ITCZ). U‐Th ages date speleothem growth from 27 to 14 ka. δ18O, δ13C, and trace metal proxies reconstruct drier conditions during Heinrich Stadials 1 and 2, and wetter conditions during the Last Glacial Maximum and Bølling–Allerød. This is surprising considering hypotheses arguing for southward (northward) ITCZ shifts during North Atlantic cooling (warming) events, which would be expected to result in wetter (drier) conditions at Anjohibe in the Southern Hemisphere tropics. The reconstructed Indian Ocean zonal (west‐east) sea surface temperature (SST) gradient is in close agreement with hydroclimate proxies in NW Madagascar, with periods of increased precipitation correlating with relatively warmer conditions in the western Indian Ocean and cooler conditions in the eastern Indian Ocean. Such gradients could drive long‐term shifts in the strength of the Walker circulation with widespread effects on hydroclimate across East Africa. These results suggest that during abrupt millennial‐scale climate changes, it is not meridional ITCZ shifts, but the tropical Indian Ocean SST gradient and Walker circulation driving East African hydroclimate variability.

Future Projections of the El Niño—Southern Oscillation and Tropical Pacific Mean State in CMIP6

AbstractThe El Niño—Southern Oscillation (ENSO) is an important mode of tropical Pacific atmosphere‐ocean variability that drives teleconnections with weather and climate globally. However, prior studies using state‐of‐the‐art climate models lack consensus regarding future ENSO projections and are often impacted by tropical Pacific sea‐surface temperature (SST) biases. We used 173 simulations from 29 climate models participating in the Coupled Model Intercomparison Project, version 6 (CMIP6) to analyze model biases and future ENSO projections. We analyzed two ENSO indices, namely the ENSO Longitude Index (ELI), which measures zonal shifts in tropical Pacific deep convection and accounts for changes in background SST, and the Niño 3.4 index, which measures SST anomalies in the central‐eastern equatorial Pacific. We found that the warm eastern tropical‐subtropical Pacific SST bias typical of previous generations of climate models persists into many of the CMIP6 models. Future projections of ENSO shift toward more El Niño‐like conditions based on ELI in 48% of simulations and 55% of models, in association with a future weakening of the zonal equatorial Pacific SST gradient. On the other hand, none of the models project a significant shift toward La Niña‐like conditions. The standard deviation of the Niño 3.4 index indicates a lack of consensus on whether an increase or decrease in ENSO variability is expected in the future. Finally, we found a possible relationship between historical SST and low‐level cloud cover biases in the ENSO region and future changes in ELI; however, this result may be impacted by limitations in data availability.

The role of in situ ocean data assimilation in ECMWF subseasonal forecasts of sea‐surface temperature and mixed‐layer depth over the tropical Pacific ocean

AbstractThe tropical Pacific plays an important role in modulating the global climate through its prevailing sea‐surface temperature spatial structure and dominant climate modes like El Niño–Southern Oscillation (ENSO), the Madden–Julian Oscillation (MJO), and their teleconnections. These modes of variability, including their oceanic anomalies, are considered to provide sources of prediction skill on subseasonal timescales in the Tropics. Therefore, this study aims to examine how assimilating in situ ocean observations influences the initial ocean sea‐surface temperature (SST) and mixed‐layer depth (MLD) and their subseasonal forecasts. We analyze two subseasonal forecast systems generated with the European Centre for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS), where the ocean states were initialized using two Observing‐System Experiment (OSE) reanalyses. We find that the SST differences between forecasts with and without ocean data assimilation grow with time, resulting in a reduced cold‐tongue bias when assimilating ocean observations. Two mechanisms related to air–sea coupling are considered to contribute to this growth of SST differences. One is a positive feedback between the zonal SST gradient, pressure gradient, and surface wind. The other is the difference in Ekman suction and mixing at the Equator due to surface wind‐speed differences. While the initial mixed‐layer depth (MLD) can be improved through ocean data assimilation, this improvement is not maintained in the forecasts. Instead, the MLD in both experiments shoals rapidly at the beginning of the forecast. These results emphasize how initialization and model biases influence air–sea interaction and the accuracy of subseasonal forecasts in the tropical Pacific.

On the Influence of the Bay of Bengal’s Sea Surface Temperature Gradients on Rainfall of the South Asian Monsoon

Abstract The southwest monsoon delivers over 70% of India’s annual rainfall and is crucial to the success of agriculture across much of South Asia. Monsoon precipitation is known to be sensitive to sea surface temperature (SST) in the Bay of Bengal (BoB). Here, we use a configuration of the Unified Model of the Met Office coupled to an ocean mixed layer model to investigate the role of upper-ocean features in the BoB on southwest monsoon precipitation. We focus on the pronounced zonal and meridional SST gradients characteristic of the BoB; the zonal gradient in particular has an as-yet unknown effect on monsoon rainfall. We find that the zonal SST gradient is responsible for a 50% decrease in rainfall over the southern BoB (approximately 5 mm day−1), and a 50% increase in rainfall over Bangladesh and northern India (approximately 1 mm day−1). This increase is remotely forced by a strengthening of the monsoon Hadley circulation. The meridional SST gradient acts to decrease precipitation over the BoB itself, similarly to the zonal SST gradient, but does not have comparable effects over land. The impacts of barrier layers and high-salinity subsurface water are also investigated, but neither has significant effects on monsoon precipitation in this model; the influence of barrier layers on precipitation is felt in the months after the southwest monsoon. Models should accurately represent oceanic processes that directly influence BoB SST, such as the BoB cold pool, in order to faithfully represent monsoon rainfall.

Near-term tropical cyclone risk and coupled Earth system model biases

Most current climate models predict that the equatorial Pacific will evolve under greenhouse gas–induced warming to a more El Niño-like state over the next several decades, with a reduced zonal sea surface temperature gradient and weakened atmospheric Walker circulation. Yet, observations over the last 50 y show the opposite trend, toward a more La Niña-like state. Recent research provides evidence that the discrepancy cannot be dismissed as due to internal variability but rather that the models are incorrectly simulating the equatorial Pacific response to greenhouse gas warming. This implies that projections of regional tropical cyclone activity may be incorrect as well, perhaps even in the direction of change, in ways that can be understood by analogy to historical El Niño and La Niña events: North Pacific tropical cyclone projections will be too active, North Atlantic ones not active enough, for example. Other perils, including severe convective storms and droughts, will also be projected erroneously. While it can be argued that these errors are transient, such that the models’ responses to greenhouse gases may be correct in equilibrium, the transient response is relevant for climate adaptation in the next several decades. Given the urgency of understanding regional patterns of climate risk in the near term, it would be desirable to develop projections that represent a broader range of possible future tropical Pacific warming scenarios—including some in which recent historical trends continue—even if such projections cannot currently be produced using existing coupled earth system models.

Abrupt change in tropical Pacific climate mean state during the Little Ice Age

The mean state of the tropical Pacific ocean-atmosphere climate, in particular its east-west asymmetry, has profound consequences for regional climates and for the El Niño/Southern Oscillation variability. Here we present a new high-resolution paleohydrological record using the stable-hydrogen-isotopic composition of terrestrial-lipid biomarkers (δDwax) from a 1,400-year-old lake sedimentary sequence from northern Philippines. Results show a dramatic and abrupt increase in δDwax values around 1630 AD with sustained high values until around 1900 AD. We interpret this change as a shift to significantly drier conditions in the western tropical Pacific during the second half of the Little Ice Age as a result of a change in tropical Pacific mean state tied to zonal sea surface temperature (SST) gradients. Our findings highlight the prominent role of abrupt shifts in zonal SST gradients on multidecadal to multicentennial timescales in shaping the tropical Pacific hydrology of the last millennium, and demonstrate that a marked transition in the tropical Pacific mean state can occur within a period of a few decades.

Interdecadal Change of Tropical Cyclone Translation Speed during Peak Season in South China Sea: Observed Evidence, Model Results, and Possible Mechanism

Abstract Long-term variations in the translation speed of tropical cyclones (TCs) in the South China Sea (SCS) are examined based on five TC datasets from different institutions. TC translation speed during the TC peak season in the SCS shows an evident rhythm of interdecadal change throughout 1977–2020. This interdecadal change in TC translation speed in the SCS can be well reproduced by a newly developed trajectory model. The model results indicate that the interdecadal change in TC translation speed is primarily due to an interdecadal change in the steering flow in the SCS. Such an interdecadal change in the steering flow is closely related to an east–west shift of the subtropical high in the western North Pacific (WNP) ocean basin, which may be driven by the zonal sea surface temperature (SST) gradient between the north Indian (NI) and WNP ocean basins. A new index of the zonal SST gradient is proposed, which is shown to be effective for indicating the interdecadal change in east–west shift of subtropical high, and thus, the TC translation speed in the SCS. Significance Statement Identification of the prevailing trend of tropical cyclone activities due to climate change has been a challenging subject of scientific studies in recent years. This research focuses on the long-term variability of the tropical cyclone translation speed in the South China Sea. It is demonstrated that there is evidently an interdecadal rhythm for variation of the tropical cyclone translation speed in the South China Sea, and this variation can be well reproduced using a newly developed trajectory model. Based on the simulation results, the possible mechanism for this variation is also elucidated.

Influence of Indian Ocean warming on rainfall of Sri Lanka

AbstractRecent studies have pointed out an increased warming over the tropical Indian Ocean (IO) and how such warming could alter the convection and rainfall in the region. In this study, using observational/reanalysis datasets from 1981 to 2020, we show that the early Maha rainfall of Sri Lanka (SL) has been increasing, and it is considerably influenced by the North‐western IO sea surface temperature (SST) increase in the second inter‐monsoon (SIM). A significant warming in the North‐western IO alters the zonal SST gradient and strengthens the low‐level circulation, which induces large‐scale convergence over the western IO during SIM, contributing to the enhancement of rainfall over SL. However, despite the substantial impact of IO warming, the moderate correlation suggests the importance of examining other factors, which may influence the rainfall variability of SL.

Biases and improvements of the boreal winter–spring equatorial undercurrent in the Indian Ocean in the CMIP5 and CMIP6 models

We assessed the performance of state-of-the-art coupled models in reproducing the equatorial undercurrent (EUC) in the Indian Ocean based on the outputs of the Coupled Model Intercomparison Project Phase 6 (CMIP6) models and compared with the Phase 5 (CMIP5) models. Our results showed that the CMIP6 models reproduced the boreal winter–spring Indian Ocean EUC more realistically than the CMIP5 models, although both generations of models underestimated the strength of the Indian Ocean EUC compared with the observations. This underestimation of the Indian Ocean EUC can be attributed to the excessively strong and westward-extended cold tongue in the equatorial Pacific. In the CMIP models, a stronger winter-mean cold tongue favors a stronger zonal sea surface temperature gradient, which forces a strong easterly wind bias over the equatorial western Pacific. This, in turn, contributes to an acceleration of the Walker circulation. This enhanced Walker circulation over the Indo-Pacific Ocean directly causes a lower level westerly wind bias over the equatorial Indian Ocean and drives a shallow west–deep east thermocline tilt bias, ultimately leading to an excessively weak EUC in the Indian Ocean via wind-induced thermocline processes. Compared with the CMIP5 models, the overall improvement in the strength of the winter–spring Indian Ocean EUC in the CMIP6 models can be traced back to the improvement in the degree of the strong and westward-extended cold tongue bias. Our results suggest that efforts should be made to reduce the bias in the mean-state equatorial Pacific sea surface temperature to further improve the simulation and projection of the atmospheric and oceanic circulations in the Indian Ocean.

The future of the El Niño–Southern Oscillation: using large ensembles to illuminate time-varying responses and inter-model differences

Abstract. Future changes in the El Niño–Southern Oscillation (ENSO) are uncertain, both because future projections differ between climate models and because the large internal variability of ENSO clouds the diagnosis of forced changes in observations and individual climate model simulations. By leveraging 14 single model initial-condition large ensembles (SMILEs), we robustly isolate the time-evolving response of ENSO sea surface temperature (SST) variability to anthropogenic forcing from internal variability in each SMILE. We find nonlinear changes in time in many models and considerable inter-model differences in projected changes in ENSO and the mean-state tropical Pacific zonal SST gradient. We demonstrate a linear relationship between the change in ENSO SST variability and the tropical Pacific zonal SST gradient, although forced changes in the tropical Pacific SST gradient often occur later in the 21st century than changes in ENSO SST variability, which can lead to departures from the linear relationship. Single-forcing SMILEs show a potential contribution of anthropogenic forcing (aerosols and greenhouse gases) to historical changes in ENSO SST variability, while the observed historical strengthening of the tropical Pacific SST gradient sits on the edge of the model spread for those models for which single-forcing SMILEs are available. Our results highlight the value of SMILEs for investigating time-dependent forced responses and inter-model differences in ENSO projections. The nonlinear changes in ENSO SST variability found in many models demonstrate the importance of characterizing this time-dependent behavior, as it implies that ENSO impacts may vary dramatically throughout the 21st century.

Accelerated mafic weathering in Southeast Asia linked to late Neogene cooling.

Arc-continent collision in Southeast Asia during the Neogene may have driven global cooling through chemical weathering of freshly exposed ophiolites resulting in atmospheric CO2 removal. Yet, little is known about the cause-and-effect relationships between erosion and the long-term evolution of tectonics and climate in this region. Here, we present an 8-million-year record of seawater chemistry and sediment provenance from the eastern Indian Ocean, near the outflow of Indonesian Throughflow waters. Using geochemical analyses of foraminiferal shells and grain size-specific detrital fractions, we show that erosion and chemical weathering of ophiolitic rocks markedly increased after 4 million years (Ma), coincident with widespread island emergence and gradual strengthening of Pacific zonal sea-surface temperature gradients. Together with supportive evidence for enhanced mafic weathering at that time from re-analysis of the seawater 87Sr/86Sr curve, this finding suggests that island uplift and hydroclimate change in the western Pacific contributed to maintaining high atmospheric CO2 consumption throughout the late Neogene.

Future Indian Ocean warming patterns

Most future projections conducted with coupled general circulation models simulate a non-uniform Indian Ocean warming, with warming hotspots occurring in the Arabian Sea (AS) and the southeastern Indian Ocean (SEIO). But little is known about the underlying physical drivers. Here, we are using a suite of large ensemble simulations of the Community Earth System Model 2 to elucidate the causes of non-uniform Indian Ocean warming. Strong negative air-sea interactions in the Eastern Indian Ocean are responsible for a future weakening of the zonal sea surface temperature gradient, resulting in a slowdown of the Indian Ocean Walker circulation and the generation of southeasterly wind anomalies over the AS. These contribute to anomalous northward ocean heat transport, reduced evaporative cooling, a weakening in upper ocean vertical mixing and an enhanced AS future warming. In contrast, the projected warming in the SEIO is related to a reduction of low-cloud cover and an associated increase in shortwave radiation. Therefore, the regional character of air-sea interactions plays a key role in promoting future large-scale tropical atmospheric circulation anomalies with implications for society and ecosystems far outside the Indian Ocean realm.

Record-breaking 2020 summer marine heatwaves in the western North Pacific

Record-breaking marine heatwaves (MHWs) occurred in the western North Pacific during the summer of 2020. These unprecedented MHWs were consistent with favorable large-scale conditions that are linked to an anomalous western North Pacific subtropical high (WNPSH), resulting mainly from sea surface temperature (SST) anomalies across the tropical oceans. In addition, a moderate La Niña-like pattern was also conducive to transporting warm seawater to the western North Pacific. Mixed-layer heat budgets suggest that surface heat flux contributed to the SST anomaly in the western subtropical Pacific. In contrast, oceanic heat advection dominanted in the South China Sea and the western equatorial Pacific. Numerical model experiments indicated that the tropical Indian Ocean SST anomalies were responsible for the enhanced WNPSH. The increased zonal SST gradient across the tropical Pacific also played an important role. Inter-ocean interactions can modulate climate variability through ocean-atmospheric coupling and deserve more attention when predicting MHWs within the context of global warming. In addition, it is critical to consider MHWs as a powerful tool in detecting acute, intense thermal stress events in the coral bleaching pre-warning system.

Storylines of Maritime Continent dry period precipitation changes under global warming

The dry half of the year from May to October over the Maritime Continent (MC) has experienced unprecedented damages from forest fires in recent decades. The observed interannual rainfall variability during this period is closely tied to sea surface temperature (SST) variability over the equatorial Pacific (EP). Therefore, the future evolution of EP SST can be expected to influence the climatological precipitation over the MC. Whilst multi-model means (MMMs) suggest a future drying trend over the south-western part of the MC, there is considerable model uncertainty. Here, using a storyline approach with the 38 climate models from Coupled Model Intercomparison Project Phase 6, we distinguish the model uncertainty associated with changes in the zonal EP SST gradient from that associated with the basin-wide EP (BEP) warming. We find that an increase in east-to-west EP SST gradient would bring more rainfall over the north-eastern regions including northern Borneo, Sulawesi and New Guinea. In contrast, the intensity of the basin-wide warming of EP SST is directly linked with the drying response seen over the south-western MC in the MMM. This drying affects the highly vulnerable regions of Sumatra and Kalimantan for forest fires. Our results suggest that a storyline under higher BEP warming accompanied by an El-Niño like change in zonal SST gradient would lead to even drier climatic conditions over these key regions. However, the observed record of more than one hundred years favours a storyline of lower BEP warming accompanied by a La-Niña like change in zonal SST gradient, which would lead to minimal drying over the south-western MC and wetter conditions over the north-eastern parts of the MC.

Greenhouse warming and internal variability increase extreme and central Pacific El Niño frequency since 1980

El Niño has been recorded to change its properties since the 1980s, characterized by more common extreme El Niño and Central Pacific (CP) El Niño events. However, it is still unclear whether such change is externally forced or part of the natural variability. Here, we find that the frequency of the extreme and CP El Niño events also increased during the period 1875–1905, when the anthropogenic CO2 concentration was relatively lower, but with a positive phase of the Atlantic Multidecadal Oscillation (AMO). Models and palaeoclimate proxies reveal that a positive AMO enhances the zonal sea surface temperature gradient in the CP, which strengthens zonal advective feedback, favoring extreme and CP El Niño development. Moreover, we estimate that internal variability contributed to ~65% of the increasingly extreme and CP El Niño events, while anthropogenic forcing has made our globe experience ~1 more extreme and ~2 more CP events over the past four decades.