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

Abstract

<div class="page" title="Page 1"> <div class="section"> <div class="layoutArea"> <div class="column"> <p>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. Over the tropical Pacific, less carbon is released 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 CO<sub>2</sub> 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–CO<sub>2</sub> flux relationship. This gradual shift, which occurs as early as the first half of the 21st century, is associated with a high CO<sub>2</sub>-induced increase in the Revelle factor that leads to stronger sensitivity of partial pressure of CO<sub>2</sub> (pCO<sub>2</sub>) to changes in surface temperature between ENSO phases. At the same time, uptake of anthropogenic CO<sub>2</sub> 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–CO<sub>2</sub> 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 CO<sub>2</sub> fluxes. Models that simulate shifts in the ENSO–CO<sub>2</sub> flux relationship simulate positive bias in surface carbonate concentrations.</p> </div> </div> </div> </div>

Equatorial Pacific pCO2 Interannual Variability in CMIP6 Models

“…For example, Vaittinada Ayar et al. (2022) shows that IPSL‐CM6A‐LR doesn't have a strong alkalinity bias, however, we find that its pCO 2, nonT : pCO 2, T variance ratio is weaker than the ratio in MRI‐ESM2‐0 (Figure 7b), which is a model they show with a strong alkalinity bias.…”

“…Previous studies of historical simulations from the earlier CMIP5 found that FCO 2 IAV were also underestimated in models (Dong et al., 2016, 2017). Results from another CMIP6 study also find that most models simulate weak FCO 2 anomalies while overestimating SST IAV (Vaittinada Ayar et al., 2022).…”

“…The latest generation of ESMs from CMIP6 have demonstrated progress in representing the mean state of ocean biogeochemistry (Séférian et al., 2020). However, as in CMIP5, weak FCO 2 IAV were also found in CMIP6 (Vaittinada Ayar et al., 2022). Identifying sources of model biases in FCO 2 IAV for the contemporary period, where some data constraints exist, is a first step toward model improvements.…”

“…Vaittinada Ayar et al. (2022) find that models with strong alkalinity biases have weak surface DIC biases (i.e., weak surface DIC variability), which leads to a reduction in pCO 2, nonT variability. They find that for some models (CanESM5, GFDL‐CM4 and MRI‐ESM2‐0), pCO 2, nonT variability is weak enough that pCO T variability can dominate total pCO 2 anomalies.…”

“…Vaittinada Ayar et al. (2022) proposed that models without a strong alkalinity bias may be better predictors of future ENSO‐CO 2 flux dynamics. However, we find that these models underestimate equatorial Pacific pCO 2 IAV and ENSO‐related covariability.…”