Contrasting projections of the ENSO-driven CO2 flux variability in the equatorial Pacific under high-warming scenario

The El Niño-Southern Oscillation (ENSO) widely modulates the global carbon cycle. More specifically, it alters the net uptake of carbon in the tropical ocean. Indeed, over the tropical Pacific less carbon is released by oceans during El Niño, while the opposite is the case for La Niña. Here, the skill of Earth system models (ESMs) from the latest Coupled Model Intercomparison Project (CMIP6) to simulate the observed tropical Pacific CO2 flux variability in response to ENSO is assessed. The temporal amplitude and spatial extent of CO2 flux anomalies vary considerably among models, while the surface temperature signals of El Niño and La Niña phases are generally well represented. Under historical conditions followed by the high-warming Shared Socio-economic Pathway (SSP5-8.5) scenarios, about half the ESMs simulate a reversal in ENSO-CO2 flux relationship. This gradual shift, which occurs as early as the first half of the 21st century, is associated with a high CO2-induced increase in the Revelle factor that leads to stronger sensitivity of partial pressure of CO2 (pCO2) to changes in surface temperature between ENSO phases. At the same time, uptake of anthropogenic CO2 substantially increases upper-ocean dissolved inorganic carbon (DIC) concentrations (reducing its vertical gradient in the thermocline) and weakens the ENSO-modulated surface DIC variability. The response of the ENSO-CO2 flux relationship to future climate change is sensitive to the contemporary mean state of the carbonate ion concentration in the tropics. We present an emergent constraint between the simulated contemporary carbonate concentration with the projected cumulated CO2 fluxes. Models that simulate shifts in the ENSO-CO2 flux relationship simulate positive bias in surface carbonate concentrations.

Figure caption: . Average contemporary surface CO32- concentration (in µmol C L-1) plotted against the cumulated sea-air CO2 fluxes (in Pg C) in the EP region from 1985 to 2014 (a) and 1850 to 2100 (b). ρ is the correlation, and reversed ESMs are marked in yellow, while preserved ESMs are marked in blue. The observations are given as brown lines, with dashed lines being the carbonate observation error.