The optimally growing initial errors (OGEs) of El Niño events are found in the Community Earth System Model (CESM) by the conditional nonlinear optimal perturbation (CNOP) method. Based on the characteristics of low-dimensional attractors for ENSO (El Niño Southern Oscillation) systems, we apply singular vector decomposition (SVD) to reduce the dimensions of optimization problems and calculate the CNOP in a truncated phase space by the differential evolution (DE) algorithm. In the CESM, we obtain three types of OGEs of El Niño events with different intensities and diversities and call them type-1, type-2 and type-3 initial errors. Among them, the type-1 initial error is characterized by negative SSTA errors in the equatorial Pacific accompanied by a negative west–east slope of subsurface temperature from the subsurface to the surface in the equatorial central-eastern Pacific. The type-2 initial error is similar to the type-1 initial error but with the opposite sign. The type-3 initial error behaves as a basin-wide dipolar pattern of tropical sea temperature errors from the sea surface to the subsurface, with positive errors in the upper layers of the equatorial eastern Pacific and negative errors in the lower layers of the equatorial western Pacific. For the type-1 (type-2) initial error, the negative (positive) temperature errors in the eastern equatorial Pacific develop locally into a mature La Niña (El Niño)-like mode. For the type-3 initial error, the negative errors in the lower layers of the western equatorial Pacific propagate eastward with Kelvin waves and are intensified in the eastern equatorial Pacific. Although the type-1 and type-3 initial errors have different spatial patterns and dynamic growing mechanisms, both cause El Niño events to be underpredicted as neutral states or La Niña events. However, the type-2 initial error makes a moderate El Niño event to be predicted as an extremely strong event.