ENSO skewness hysteresis and associated changes in strong El Nino under a CO2 removal scenario

El Nino-Southern Oscillation (ENSO) sea surface temperature (SST) anomaly skewness encapsulates the nonlinear processes of strong ENSO events and affects future climate projections. Yet, its response to CO2 forcing remains not well understood. Here, we find ENSO skewness hysteresis in a large ensemble CO2 removal simulation. The positive SST skewness in the central-to-eastern tropical Pacific gradually weakens (most pronounced near the dateline) in response to increasing CO2, but weakens even further once CO2 is ramped down. Further analyses reveal that hysteresis of the Intertropical Convergence Zone migration leads to more active and farther eastward-located strong eastern Pacific El Nino events, thus decreasing central Pacific ENSO skewness by reducing the amplitude of the central Pacific positive SST anomalies and increasing the scaling effect of the eastern Pacific skewness denominator, i.e., ENSO intensity, respectively. The reduction of eastern Pacific El Nino maximum intensity, which is constrained by the SST zonal gradient of the projected background El Nino-like warming pattern, also contributes to a reduction of eastern Pacific SST skewness around the CO2 peak phase. This study highlights the divergent responses of different strong El Nino regimes in response to climate change.