“…Past analyses have used a Linear Inverse Model (LIM; Penland & Sardeshmukh, , hereafter PS95), to empirically extract a linear dynamical system forced with state‐ independent (i.e., additive) white noise from the covariance statistics of Tropical Pacific anomalies. Such a LIM can still explain ENSO irregularity (PS95), the (nonnormal) cycle of ENSO decay and growth (PS95, Vimont et al, ), the spectral characteristics of the main ENSO indices (Ault et al, ; Newman et al, ) and has also been shown to have forecast skill comparable to that of fully coupled, nonlinear General Circulation Models (Newman & Sardeshmukh, ), suggesting that it provides a good approximation to the Tropical Pacific deterministic dynamics at time scales relevant for ENSO (PS95; see also Penland, ). Moreover, inverse models augmented with deterministic nonlinearity do not seem to improve forecast skill over their linear counterparts (Chen et al, ; Kondrashov et al, ), which suggests that nonlinear processes contribute mostly to the unpredictable part of the system.…”