Diagnosing the Quasi-Equilibrium Response of ENSO Variability under a Range of CO2 Levels

Several recent studies have highlighted differences in simulated properties of El Nino-Southern Oscillation (ENSO) under transient and equilibrium responses to increasing CO2. However, the reasons behind these disparate responses and the extent to which they are robust to different scales of CO2 forcing remain unclear. In this study, we adopt a climate system model with reduced SST bias in the eastern tropical Pacific and incrementally apply abrupt increases in CO2, analyzing outputs after each simulation reaches quasi-equilibrium with the imposed forcing. The results suggest that ENSO activity under quasi-equilibrium first increases and then decreases with increasing CO2, peaking in simulations with CO2 concentrations similar to the present day. Bjerknes-Jin stability analysis indicates that changes in the ENSO growth rate result primarily from changes in the thermocline feedback and thermodynamic damping terms. While thermodynamic damping increases monotonically with increasing CO2, the positive thermocline feedback varies within the range of internal variability up to twice the preindustrial value of CO2 and then weakens sharply with further increases. The mechanisms behind these changes include weaker mean ocean upwelling and weaker dynamical coupling between the atmosphere and subsurface ocean associated with substantial near-surface freshening at higher levels of CO2. These changes steepen the thermodynamic barrier to mixing between the surface and subsurface, weakening the east-west temperature gradient in the mean state and suppressing variability in the cold tongue. Analysis of similar model simulations from the Coupled Model Intercomparison Project (CMIP6) archive indicates that changes in the Bjerknes-Jin stability index are robust but do not establish a consensus as to the mechanisms behind them.