Abstract: Based on the subsurface temperature parameterization theoretical framework and the Zebiak-Cane (ZC) model, using the monthly reanalysis of Phase 3.3.1 of the Simple Ocean Data Assimilation provided by the US University of Maryland and the monthly sea surface temperature dataset from the Hadley Center for the 1980—2015 period, the role of subsurface temperature (T_sub) and its key controlling factors in the ENSO amplitude asymmetry are investigated, with a focus on contrasting mechanisms across the eastern and central-western equatorial Pacific. It is find that the asymmetry extends from the surface to the subsurface and exhibits a dipole structure of "positive in the east and negative in the west". This dipole primarily arises from asymmetric subsurface responses to zonal wind stress (τ), but the dominant mechanisms differ across the equatorial Pacific. In the eastern equatorial Pacific, the positive asymmetry results from combined changes in thermocline depth (h) and thermocline sharpness (h^*). The stronger positive τ′ during the warm phase produces larger positive h′ and negative h^*′, leading to enhanced positive T'_\rmsub . In the central-western Pacific, the negative asymmetry is mainly controlled by h^*. The weaker negative τ′ during the cold phase actually induces stronger positive h^*′, which generates stronger negative T'_\rmsub . Whereas h′ contributes little in the western Pacific because the mean thermocline is much deeper than the mixed-layer depth. Overall, the asymmetric response of h^* to τ is the primary driver of the negative ENSO amplitude asymmetry in the central-western equatorial Pacific. By replacing the T_sub parameterization with a new scheme based on the h^* first-order approximation, the simulation of ENSO amplitude asymmetry in ZC model has been significantly improved.