Simple stochastic dynamical models capturing the statistical diversity of El Nino Southern Oscillation

The El Nino Southern Oscillation (ENSO) has significant impact on global climate and seasonal prediction. A simple modeling framework is developed here that automatically captures the statistical diversity of ENSO. First, a stochastic parameterization of the wind bursts including both westerly and easterly winds is coupled to a simple ocean-atmosphere model that is otherwise deterministic, linear, and stable. Second, a simple nonlinear zonal advection with no ad hoc parameterization of the background sea-surface temperature (SST) gradient and a mean easterly trade wind anomaly representing the multidecadal acceleration of the trade wind are both incorporated into the coupled model that enables anomalous warm SST in the central Pacific. Then a three-state stochastic Markov jump process is used to drive the wind burst activity that depends on the strength of the western Pacific warm pool in a simple and effective fashion. It allows the coupled model to simulate the quasi-regular moderate traditional El Nino, the super El Nino, and the central Pacific (CP) El Nino as well as the La Nina with realistic features. In addition to the anomalous SST, the Walker circulation anomalies at different ENSO phases all resemble those in nature. In particular, the coupled model succeeds in reproducing the observed episode during the 1990s, where a series of 5-y CP El Ninos is followed by a super El Nino and then a La Nina. Importantly, both the variance and the non-Gaussian statistical features in different Nino regions spanning from the western to the eastern Pacific are captured by the coupled model.