Subsurface Heat Content Research Articles

Interhemispheric synchrony of Mid-Late Holocene SST in the Indo-Pacific warm pool linked to upper water column dynamics

The sea surface temperature (SST) of the Indo-Pacific warm pool (IPWP) plays a crucial role in global climate system. Despite its importance, knowledge of its SST changes over century to millennial timescales remains limited and controversial, hindering our comprehensive understanding of climate change mechanisms. Our study bridges this gap with a novel SST reconstruction spanning approximately 7000 years, derived by aggregating 75 coral Sr/Ca-SST records from the South China Sea. These records reveal a notable synchronization in SST variations with foraminiferal reconstructions across the IPWP, reflecting a prevalent pattern of interhemispheric cooling and multi-century variability in the IPWP during the Mid-Late Holocene. Further analysis reveals that these (sub)millennial SST variations are closely linked to subsurface water temperature and ocean heat content (OHC) changes, supporting the theory that the release of deeper OHC acted as a heat source for low-latitude processes during the Holocene.

Investigating Multidecadal Trends in Ice Cover and Subsurface Temperatures in the Laurentian Great Lakes Using a Coupled Hydrodynamic–Ice Model

Abstract While changing lake surface conditions have received significant research scrutiny, changes in subsurface conditions, including stratification and heat content, remain largely unexplored. In this work, we highlight changes in thermal structure, stratification dynamics, and ice characteristics in the Laurentian Great Lakes (Lakes Superior, Huron, Michigan, Erie, and Ontario) as simulated between 1979 and 2021. Three-dimensional lake hydrodynamics and ice cover were modeled using the Finite Volume Community Ocean Model (FVCOM) coupled with the Los Alamos sea ice model (CICE). Analysis revealed significant increases in surface (0.4°–0.6°C decade−1) and subsurface (0.1°–0.4°C decade−1) temperatures as well as dramatic losses in ice cover (1%–8% decade−1) and ice volume (0–3 km3 decade−1) over the last 40 years. Estimated surface heating rates were strongest during the summer and fall, while subsurface warming was most rapid during the nearly isothermal winter and spring. Intensified (decreased) summer (winter) stratification led to shifts in lake turnover dynamics, with delayed fall turnover dates (2–6 days decade−1) and earlier spring overturn dates (2–9 days decade−1). Modeled surface temperatures (LST), bottom temperatures (LBT), and annual averaged ice cover (AAIC) were used to estimate low-frequency climate signals, which were highly correlated with the Atlantic multidecadal oscillation. Warming trends fit to residual climate signals (LST: 0.1°C decade−1; LBT: 0.03°C decade−1; AAIC: −1% decade−1), calculated by removing low-frequency variability from the raw climate signal, were lower than those fit to associated low-frequency components, suggesting that recent climate change in the Great Lakes may be strongly influenced by natural multidecadal climate variability.

A meridional dipole mode in the Indian Ocean subsurface ocean heat content and its multidecadal variability

The dominant mode of the Indian Ocean variability can be impacted by changes in the background state (multidecadal timescale) of the subsurface heat content. The multidecadal variability of subsurface ocean heat content (sub-OHC) in the Indian Ocean is examined using four reanalysis products from 1958 to 2017. The analysis reveals a meridional basin-wide dipole mode in the subsurface OHC until the late 1980s, followed by the mode embedded in the uniform basin-wide patterns. These patterns are also observed in the trends of thermocline and sea surface height. The observed patterns in the Indian Ocean are explained by two distinct mechanisms. Firstly, the multidecadal variability of dipole patterns over the Indian Ocean is influenced by local wind forcing. Wind stress trends and Ekman pumping velocity trends favor downwelling (upwelling) in the off-equatorial southern Indian Ocean region, leading to thermocline depth deepening (shallowing) during 1958–1975 and 1976–1987, respectively. Secondly, the combined effect of heat transport from the western Pacific through Indonesian Through Flow and local wind forcing accounts for the basin-wide cooling and warming trends observed during 1988–2000, and 2001–2014, respectively.

ENSO phase space dynamics with an improved estimate of the thermocline depth

The recharge oscillator model of the El Niño Southern Oscillation (ENSO) describes the ENSO dynamics as an interaction between the eastern tropical Pacific sea surface temperatures (T) and subsurface heat content (thermocline depth; h), defining a dynamical cycle with different phases. h is often approximated on the basis of the depth of the 20 °C isotherm (Z20). In this study we will address how the estimation of h affects the representation of ENSO dynamics. We will compare the ENSO phase space with h estimated based on Z20 and based on the maximum gradient in the temperature profile (Zmxg). The results illustrate that the ENSO phase space is much closer to the idealised recharge oscillator model if based on Zmxg than if based on Z20. Using linear and non-linear recharge oscillator models fitted to the observed data illustrates that the Z20 estimate leads to artificial phase dependent structures in the ENSO phase space, which result from an in-phase correlation between h and T. Based on the Zmxg estimate the ENSO phase space diagram show very clear non-linear aspects in growth rates and phase speeds. Based on this estimate we can describe the ENSO cycle dynamics as a non-linear cycle that grows during the recharge and El Nino state, and decays during the La Nina states. The most extreme ENSO states are during the El Nino and discharge states, while the La Nina and recharge states do not have extreme states. It further has faster phase speeds after the El Nino state and slower phase speeds during and after the La Nina states. The analysis suggests that the ENSO phase speed is significantly positive in all phases, suggesting that ENSO is indeed a cycle. However, the phase speeds are closest to zero during and after the La Nina state, indicating that the ENSO cycle is most likely to stall in these states.

Unexpected limitation of tropical cyclone genesis by subsurface tropical central-north Pacific during El Niño

The vast tropical Pacific is home to the majority of tropical cyclones (TCs) which threaten the rim countries every year. The TC genesis is nourished by warm sea surface temperatures (SSTs). During El Niño, the western Pacific warm pool extends eastward. However, the number of TCs does not increase significantly with the expanding warm pool and it remains comparable between El Niño and La Niña. Here, we show that the subsurface heat content change counteracts the favorable SSTs in the tropical central-north Pacific. Due to the anomalous positive wind stress curl, the 26 °C isotherm shoals during El Niño over this region and the heat content diminishes in the tropical central-north Pacific, even though warm SST anomalies prevail. This negative correlation between SST and 26 °C isotherm depth anomalies is opposite to the positive correlation in the tropical eastern and western Pacific. This is critical because quantifying the dynamics of the subsurface ocean provides insight into TC genesis. The trend in TC genesis continues to be debated. Future projections must account for the net effect of the surface-subsurface dynamics on TCs, especially given the expected El Niño-like pattern over the tropical Pacific under global warming.

Processes Controlling the Southern Ocean Temperature Change

Abstract Observations show that since the 1950s, the subsurface Southern Ocean has been rapidly warming in the south and cooling in the north. However, what drives these opposing latitudinal patterns of change is not well understood. By analyzing the outputs of atmosphere–ocean general circulation models, we examine and quantify the processes controlling the temperature change in the Southern Ocean. We find that subsurface cooling in the north results from the combined effect of weak warming driven by greenhouse gases and cooling forced by aerosols. Contrary to previous work, we find that the influence of stratospheric ozone forcing is small and unnecessary to explain the observed patterns. Analysis of heat budget diagnostics shows that ocean heat content change between 400 and 1200 m in the Southern Ocean is primarily driven by slight imbalances between the large and opposing changes in residual mean advection and isopycnal diffusion. We show that while surface heat flux anomalies are the main contributor to the warming in the Southern Ocean, changes in wind stress are also required to produce the observed pattern of temperature change. These results highlight the importance of ocean dynamics in subsurface ocean heat content change. The counterbalancing effects of aerosols and greenhouse gases on ocean heat content change in the Southern Ocean will change with aerosols declining in the future, while our results suggest that changes in stratospheric ozone forcing will have only a minor effect on future patterns of Southern Ocean thermohaline structure.

Thermal coupling of the Indo-Pacific warm pool and Southern Ocean over the past 30,000 years

The role of the tropical Pacific Ocean and its linkages to the southern hemisphere during the last deglacial warming remain highly controversial. Here we explore the evolution of Pacific horizontal and vertical thermal gradients over the past 30 kyr by compiling 340 sea surface and 7 subsurface temperature records, as well as one new ocean heat content record. Our records reveal that La Niña-like conditions dominated during the deglaciation as a result of the more intense warming in the western Pacific warm pool. Both the subsurface temperature and ocean heat content in the warm pool rose earlier than the sea surface temperature, and in phase with South Pacific subsurface temperature and orbital precession, implying that heat exchange between the tropical upper water column and the extratropical Southern Ocean facilitated faster warming in the western Pacific. Our study underscores the key role of the thermal coupling between the warm pool and the Southern Ocean and its relevance for future global warming.

The Different Relationships between the ENSO Spring Persistence Barrier and Predictability Barrier

Abstract In this paper, we investigate the relationship between the El Niño–Southern Oscillation (ENSO) spring persistence barrier (PB) and predictability barrier (PD) and apply it to explain the interdecadal modulation of ENSO prediction skill using the anomaly correlation coefficient (ACC). Previous studies showed that a longer persistence (i.e., autocorrelation) tends to produce a higher prediction skill. Using the recharge oscillator model of ENSO, both analytical and numerical solutions suggest that the predictability (i.e., ACC) is related to the persistence of sea surface temperature (SST) and cross correlation between SST and subsurface ocean heat content in the tropical Pacific. In particular, a larger damping rate in SST anomalies will lead to a lower persistence and ACC and a stronger PD. However, a shortened ENSO period, which controls the cross correlation, will lead to a lower persistence but a higher ACC associated with a weaker PD. Finally, we apply our solutions to observations and suggest that a higher ACC associated with a weaker PD after 1960 is caused by the shortened ENSO period.

Asymmetries in the ENSO phase space

El Niño Southern Oscillation (ENSO) dynamics are best described by the recharge oscillator model, in which the eastern tropical Pacific sea surface temperatures (T) and subsurface heat content (thermocline depth; h) have an out-of-phase relationship. This defines a 2-dimensional phase space diagram between T and h. In an idealized, stochastically forced damped oscillator, the mean phase space diagram should be a perfectly symmetrical circle with a clockwise propagation over time. However, the observed phase space shows strong asymmetries. In this study we illustrate how the ENSO phase space can be used to discuss the phase-dependency of ENSO dynamics. A normalized spherical coordinate system allows the definition of phase-depending ENSO growth rates and phase transition speeds. Based on these we discuss the implications of the observed asymmetries with regards to the dynamics and predictability of ENSO; with a particular focus on the variations in the growth rate and coupling of ENSO along the oscillation cycle. Using linear and non-linear recharge oscillator models we will show how dynamics and noise are driving ENSO at different phases of the ENSO cycles. The results illustrate that the ENSO cycle with positive phase transitions is present in all phases but has strong variations in its strength. Much of these variations result from presenting the ENSO phase space with estimates of h based on the iso-thermal depth, that is not ideal as it is not out-of-phase with T. Future work should address how h can be estimated better, including aspects such as the vertical temperature gradients and the meridional or zonal range. We further illustrated that a non-linear growth rate of T can explain most of the observed non-linear phase space characteristics.

Pacific subsurface temperature as a long‐range indicator of El Niño, regional precipitation, and fire

AbstractThe SubNiño4 index based on the subsurface potential temperature around the thermocline beneath the west Pacific warm pool, the Niño 4 region, is examined as a long‐range indicator of the surface El Niño–Southern Oscillation (ENSO) and ENSO‐driven atmospheric response. The SubNiño4 index captures the evolution of subsurface ocean heat content between the El Niño and La Niña phases of the ENSO cycle, allowing it to serve as a long‐range indicator of surface ENSO and hence also many ENSO‐driven atmospheric anomalies. The SubNiño4 index has more temporally stable correlations with Niño 3.4 than the widely used western equatorial Pacific warm‐water volume indicator. For a lead time of the order of 12 months, Niño 3.4 correlations afforded by the lead observed SubNiño4 index become similar to and can exceed those produced by typical dynamical ENSO predictions. The value and viability of the SubNiño4 index as a simple statistical long‐range indicator of ENSO‐driven atmospheric response is shown for regional precipitation anomalies throughout the Tropics and fires in Continental and Maritime Southeast Asia.

Reassessing seasonal sea ice predictability of the Pacific-Arctic sector using a Markov model

Abstract. In this study, a regional linear Markov model is developed to assess seasonal sea ice predictability in the Pacific-Arctic sector. Unlike an earlier pan-Arctic Markov model that was developed with one set of variables for all seasons, the regional model consists of four seasonal modules with different sets of predictor variables, accommodating seasonally varying driving processes. A series of sensitivity tests are performed to evaluate the predictive skill in cross-validated experiments and to determine the best model configuration for each season. The prediction skill, as measured by the sea ice concentration (SIC) anomaly correlation coefficient (ACC) between predictions and observations, increased by 32 % in the Bering Sea and 18 % in the Sea of Okhotsk relative to the pan-Arctic model. The regional Markov model's skill is also superior to the skill of an anomaly persistence forecast. SIC trends significantly contribute to the model skill. However, the model retains skill for detrended sea ice extent predictions for up to 7-month lead times in the Bering Sea and the Sea of Okhotsk. We find that subsurface ocean heat content (OHC) provides a crucial source of prediction skill in all seasons, especially in the cold season, and adding sea ice thickness (SIT) to the regional Markov model has a substantial contribution to the prediction skill in the warm season but a negative contribution in the cold season. The regional model can also capture the seasonal reemergence of predictability, which is missing in the pan-Arctic model.

Ocean Dynamics are Key to Extratropical Forcing of El Niño

AbstractEl Niño–Southern Oscillation (ENSO) has been recently linked with extratropical Pacific Ocean atmospheric variability. The two key mechanisms connecting the atmospheric variability of the extratropical Pacific with ENSO are the heat flux–driven “seasonal footprinting mechanism” (SFM) and the ocean dynamics–driven “trade wind charging” (TWC) mechanism. However, their relative contributions to ENSO are still unknown. Here we present modeling evidence that the positive phase of the SFM generates a weaker, short-lived central Pacific El Niño–like warming pattern in the autumn, whereas the TWC positive phase leads to a wintertime eastern Pacific El Niño–like warming. When both mechanisms are active, a strong, persistent El Niño develops. While both mechanisms can trigger equatorial wind anomalies that generate an El Niño, the strength and persistence of the warming depends on the subsurface heat content buildup by the TWC mechanism. These results suggest that while dynamical coupling associated with extratropical forcing is crucial to maintain an El Niño, thermodynamical coupling is an extratropical source of El Niño diversity.

Subsurface Evolution and Persistence of Marine Heatwaves in the Northeast Pacific

AbstractThe reappearance of a northeast Pacific marine heatwave (MHW) sounded alarms in late summer 2019 for a warming event on par with the 2013–2016 MHW known as The Blob. Despite these two events having similar magnitudes in surface warming, differences in seasonality and salinity distinguish their evolutions. We compare and contrast the ocean's role in the evolution and persistence of the 2013–2016 and 2019–2020 MHWs using mapped temperature and salinity data from Argo floats. An unusual near‐surface freshwater anomaly in the Gulf of Alaska during 2019 increased the stability of the water column, preventing the MHW from penetrating deep as strongly as the 2013–2016 event. This freshwater anomaly likely contributed to the intensification of the MHW by increasing the near‐surface buoyancy. The gradual buildup of subsurface heat content throughout 2020 in the region suggests the potential for persistent ecological impacts.

Testing the Trade Wind Charging Mechanism and Its Influence on ENSO Variability

AbstractDuring the positive phase of the North Pacific Oscillation, westerly wind anomalies over the subtropical North Pacific substantially increase subsurface heat content along the equator by “trade wind charging” (TWC). TWC provides a direct pathway between extratropical atmospheric circulation and El Niño–Southern Oscillation (ENSO) initiation. Previous model studies of this mechanism lacked the ocean–atmospheric coupling needed for ENSO growth, so it is crucial to examine whether TWC-induced heat content anomalies develop into ENSO events in a coupled model. Here, coupled model experiments, forced with TWC favorable (+TWC) or unfavorable (−TWC) wind stress, are used to examine the ENSO response to TWC. The forcing is imposed on the ocean component of the model through the first winter and then the model evolves in a fully coupled configuration through the following winter. The +TWC (−TWC) forcing consistently charges (discharges) the equatorial Pacific in spring and generates positive (negative) subsurface temperature anomalies. These subsurface temperature anomalies advect eastward and upward along the equatorial thermocline and emerge as like-signed sea surface temperature (SST) anomalies in the eastern Pacific, creating favorable conditions upon which coupled air–sea feedback can act. During the fully coupled stage, warm SST anomalies in +TWC forced simulations are amplified by coupled feedbacks and lead to El Niño events. However, while −TWC forcing results in cool SST anomalies, pre-existing warm SST anomalies in the far eastern equatorial Pacific persist and induce local westerly wind anomalies that prevent consistent development of La Niña conditions. While the TWC mechanism provides adequate equatorial heat content to fuel ENSO development, other factors also play a role in determining whether an ENSO event develops.

A method for constructing synthetic parallel tracks for augmenting in-situ data with satellite observations

In-situ observations have significant importance for calibration and validation of satellite data, processes studies etc, however in-situ observations are often available with less spatio-temporal coverage due to practical limitations in making observations, as well as cost effectiveness. Albeit satellite observations have better spatial and temporal coverage, they have certain limitations like obstacles of cloud coverage, retrieval algorithms etc. Hence worthy inferences on the processes (that was intended to study) can be obtained only through skillful utilization of both in-situ observations, satellite and/or reanalysis products. We have illustrated a method for deriving synthetic parallel tracks to in-situ track/survey. Using this method, the satellite and/or reanalysis products can be extracted over these synthetic parallel tracks to derive worthy understanding/conclusions on the processes. The points on the original track are reproduced in terms of slope and distance, which can be utilized to construct a new synthetic parallel track initiating at a predetermined distance to the starting point of original track. As an example, we have demonstrated variations in oceanic physical properties along synthetic parallel tracks to the track of the Khai-muk cyclone. It is easy to understood using this technique that the variations in sub-surface heat content (as evidenced by variations in sea level along Khai-muk cyclone track) have a relatively major role in controlling the intensity of the system, compared to the changes in surface temperature. The method illustrated here is applicable to any geophysical track/survey with varying spatial and/or temporal extension.

The Pacific Meridional Mode and ENSO: a Review

This paper reviews recent progress in understanding of the North Pacific Meridional Mode (NPMM) and its influence on the timing, magnitude, flavor, and intensity of the El Nino-Southern Oscillation (ENSO). The NPMM is a seasonally evolving mode of coupled climate variability and features several distinct opportunities to influence ENSO. They include: (1) A Wind-Evaporation-SST (WES) feedback-driven propagation of surface anomalies onto the equator during boreal spring, (2) Trade Wind Charging (TWC) of equatorial subsurface heat content by NPMM-related surface wind stress curl anomalies in boreal winter and early spring, (3) The reflection of NPMM-forced ocean Rossby waves off the western boundary in boreal summer, and (4) A Gill-like atmospheric response associated with anomalous deep convection in boreal summer and fall. The South Pacific Meridional Mode (SPMM) also significantly modulates ENSO, and its interactions with the NPMM may contribute to ENSO diversity. Together, the NPMM and SPMM are also important components of Tropical Pacific Decadal Variability; however, future research is needed to improve understanding on these timescales. Since 1950, the boreal spring NPMM skillfully predicts about 15–30% of observed winter ENSO variability. Improving simulated NPMM-ENSO relationships in forecast models may reduce ENSO forecasting error. Recent studies have begun to explore the influence of anthropogenic climate change on the NPMM-ENSO relationship; however, the results are inconclusive.

Assessment of the interannual variability of local atmospheric and ITF contribution to the subsurface heat content of southern tropical Indian Ocean in GECCO2 and ORAS4 using ROMS

Abstract Atmospheric and Ocean reanalysis products are handy tools for studying the past climate. However, the choice of these products could influence the outcome of the study. In this study, we compared two popular Ocean reanalysis products (ORAS4 and GECCO2) in the context of the Indian Ocean (IO) heat content variability. The region of interest is the southern tropical IO (STIO) as the region features a seasonal upwelling zone and the Indonesian Throughflow (ITF) is nearly aligned with it. The comparison mainly includes the trends in sea surface height anomaly, IO response to Interdecadal Pacific Oscillation (IPO), and the recent changing behavior of Indo-Pacific Oceans. Among the drivers of southern tropical Indian Ocean subsurface heat content (STIOSHC), the local atmospheric wind forcing seems to be contributing more to its variability in ORAS4 where as ITF seems to be contributing more in GECCO2. To get a deeper understanding of the local and remote influences, we distinguish the interannual variability in the STIOSHC due to local atmospheric contribution and ITF contribution by setting up an Indian Ocean model using ROMS and carrying out some experiments. We show that the interannual variability in STIOSHC contributed by local winds is higher in ORAS4 when compared to GECCO2 and it is mainly reflected in the western IO. In the eastern IO, in ORAS4, the local wind and ITF contributions to the STIOSHC variability are found to be comparable whereas in GECCO2, the ITF contribution is greater than that of ORAS4. In GECCO2, the local wind and ITF contributions are found to be more or less same whereas in ORSA4, the local wind contribution is noticeably higher than the ITF contribution. Also, the local wind and ITF contribution to STIOSHC are found to be out of phase with dominant time periods of 8 to 15 years. The high STIOSHC during first (second) half of the warming hiatus period is found to be contributed by ITF (local surface forcing). The sea surface height anomalies and sea surface temperature anomalies (SSTA) in ORAS4 at few key locations in the Indian Ocean are found to be more in line with observations than those of GECCO2. We posit that ORAS4 is more reliable than GECCO2 over the IO.

Distinct Mechanisms of Decadal Subsurface Heat Content Variations in the Eastern and Western Indian Ocean Modulated by Tropical Pacific SST

Decadal variability of the subsurface ocean heat content (OHC) in the Indian Ocean is investigated using a coupled climate model experiment, in which observed eastern tropical Pacific sea surface temperature (EPSST) anomalies are specified. This study intends to understand the contributions of external forcing relative to those of internal variability associated with EPSST, as well as the mechanisms by which the Pacific impacts Indian Ocean OHC. Internally generated variations associated with EPSST dominate decadal variations in the subsurface Indian Ocean. Consistent with ocean reanalyses, the coupled model reproduces a pronounced east–west dipole structure in the southern tropical Indian Ocean and discontinuities in westward-propagating signals in the central Indian Ocean around 100°E. This implies distinct mechanisms by which the Pacific impacts the eastern and western Indian Ocean on decadal time scales. Decadal variations of OHC in the eastern Indian Ocean are attributed to 1) western Pacific surface wind anomalies, which trigger oceanic Rossby waves propagating westward through the Indonesian Seas and influence Indonesian Throughflow transport, and 2) zonal wind anomalies over the central tropical Indian Ocean, which trigger eastward-propagating Kelvin waves. Decadal variations of OHC in the western Indian Ocean are linked to conditions in the Pacific via changes in the atmospheric Walker cell, which trigger anomalous wind stress curl and Ekman pumping in the central tropical Indian Ocean. Westward-propagating oceanic Rossby waves extend the influence of this anomalous Ekman pumping to the western Indian Ocean.

Eastern Boundary Circulation and Hydrography Off Angola: Building Angolan Oceanographic Capacities

AbstractThe eastern boundary region off Angola encompasses a highly productive ecosystem important for the food security of the coastal population. The fish-stock distribution, however, undergoes large variability on intraseasonal, interannual, and longer time scales. These fluctuations are partly associated with large-scale warm anomalies that are often forced remotely from the equatorial Atlantic and propagate southward, reaching the Benguela upwelling off Namibia. Such warm events, named Benguela Niños, occurred in 1995 and in 2011. Here we present results from an underexplored extensive in situ dataset that was analyzed in the framework of a capacity-strengthening effort. The dataset was acquired within the Nansen Programme executed by the Food and Agriculture Organization of the United Nations and funded by the Norwegian government. It consists of hydrographic and velocity data from the Angolan continental margin acquired biannually during the main downwelling and upwelling seasons over more than 20 years. The mean seasonal changes of the Angola Current from 6° to 17°S are presented. During austral summer the southward Angola Current is concentrated in the upper 150 m. It strengthens from north to south, reaching a velocity maximum just north of the Angola Benguela Front. During austral winter the Angola Current is weaker, but deeper reaching. While the southward strengthening of the Angola Current can be related to the wind forcing, its seasonal variability is most likely explained by coastally trapped waves. On interannual time scales, the hydrographic data reveal remarkable variability in subsurface upper-ocean heat content. In particular, the 2011 Benguela Niño was preceded by a strong subsurface warming of about 2 years’ duration.

Influences of Pacific Climate Variability on Decadal Subsurface Ocean Heat Content Variations in the Indian Ocean

Abstract Decadal variabilities in Indian Ocean subsurface ocean heat content (OHC; 50–300 m) since the 1950s are examined using ocean reanalyses. This study elaborates on how Pacific variability modulates the Indian Ocean on decadal time scales through both oceanic and atmospheric pathways. High correlations between OHC and thermocline depth variations across the entire Indian Ocean Basin suggest that OHC variability is primarily driven by thermocline fluctuations. The spatial pattern of the leading mode of decadal Indian Ocean OHC variability closely matches the regression pattern of OHC on the interdecadal Pacific oscillation (IPO), emphasizing the role of the Pacific Ocean in determining Indian Ocean OHC decadal variability. Further analyses identify different mechanisms by which the Pacific influences the eastern and western Indian Ocean. IPO-related anomalies from the Pacific propagate mainly through oceanic pathways in the Maritime Continent to impact the eastern Indian Ocean. By contrast, in the western Indian Ocean, the IPO induces wind-driven Ekman pumping in the central Indian Ocean via the atmospheric bridge, which in turn modifies conditions in the southwestern Indian Ocean via westward-propagating Rossby waves. To confirm this, a linear Rossby wave model is forced with wind stresses and eastern boundary conditions based on reanalyses. This linear model skillfully reproduces observed sea surface height anomalies and highlights both the oceanic connection in the eastern Indian Ocean and the role of wind-driven Ekman pumping in the west. These findings are also reproduced by OGCM hindcast experiments forced by interannual atmospheric boundary conditions applied only over the Pacific and Indian Oceans, respectively.