Abstract Variations in eastern Indian Ocean upper-ocean thermal properties are assessed for the period 1970–2004, with a particular focus on asymmetric features related to opposite phases of Indian Ocean dipole events, using high-resolution ocean model hindcasts. Sensitivity experiments, where interannual atmospheric forcing variability is restricted to the Indian or Pacific Ocean only, support the interpretation of forcing mechanisms for large-scale asymmetric behavior in eastern Indian Ocean variability. Years are classified according to eastern Indian Ocean subsurface heat content (HC) as proxy of thermocline variations. Years characterized by an anomalous low HC feature a zonal gradient in upper-ocean properties near the equator, while high events have a meridional gradient from the tropics into the subtropics. The spatial and temporal characteristics of the seasonal evolution of HC anomalies for the two cases is distinct, as is the relative contribution from Indian Ocean atmospheric forcing versus remote influences from Pacific wind forcing: low events develop rapidly during austral winter/spring in response to Indian Ocean wind forcing associated with an enhanced southeasterly monsoon driving coastal upwelling and a shoaling thermocline in the east; in contrast, formation of an anomalous high eastern Indian Ocean HC is more gradual, with anomalies earlier in the year expanding from the Indonesian Throughflow (ITF) region, initiated by remote Pacific wind forcing, and transmitted through the ITF via coastal wave dynamics. Implications for seasonal predictions arise with high HC events offering extended lead times for predicting thermocline variations and upper-ocean properties across the eastern Indian Ocean.