In this study, singular vector analysis was performed for the period from 1856 to 2003 using the latest Zebiak-Cane model version LDEO5. The singular vector, representing the optimal growth pattern of initial perturbations/errors, was obtained by perturbing the constructed tangent linear model of the Zebiak-Cane model. Variations in the singular vector and singular value, as a function of initial time, season, ENSO states, and optimal period, were investigated. Emphasis was placed on exploring relative roles of linear and nonlinear processes in the optimal perturbation growth of ENSO, and deriving statistically robust conclusions using long-term singular vector analysis. It was found that the first singular vector is dominated by a west-east dipole spanning most of the equatorial Pacific, with one center located in the east and the other in the central Pacific. Singular vectors are less sensitive to initial conditions, i.e., independence of seasons and decades; while singular values exhibit a strong sensitivity to initial conditions. The dynamical diagnosis shows that the total linear and nonlinear heating terms play opposite roles in controlling the optimal perturbation growth, and that the linear optimal perturbation is more than twice as large as the nonlinear one. The total linear heating causes a warming effect and controls two positive perturbation growth regions: one in the central Pacific and the other in the eastern Pacific; whereas the total linearized nonlinear advection brings a cooling effect controlling the negative perturbation growth in the central Pacific.
