Observed evidence of the IO-NA relationship
To quantify the strength of the IO-NA warming chain, we defined the IO SST index (IO index) and the NA Ocean SST index (NA index) by using the SST averaged over each box (Fig. 1a). From 1970–2019, the IO index showed a warming trend (1.40°C per century), with strong multi-decadal fluctuations (Fig. 1d). Local maximum and minimum IO indices were observed close to 1940 and 1970, respectively. The NA index also showed a significant warming trend (2.02°C per century) and a maximum inter-decadal fluctuation close to 1940. The warming trends and interdecadal-to-multidecadal fluctuations in the IO and NA indices tend to overlap. The lead-lag correlation coefficient between the IO and NA indices (Fig. 1f) was positive regardless of lag. The maximum correlation coefficient was recorded at zero lag (0.80), indicating concurrent NA and IO warming. Global warming caused by the anthropogenic increase in CO2 concentrations may explain warming trends in the two indices; however, global warming cannot fully explain the phase matching in interdecadal-to-multidecadal variations. The warming trend in both regions can be modified by the strong interactions between the two basins. Therefore, we conjecture that the IO and NA might serve as mutual pacemakers.
To further analyse the IO-NA relationship, we examined global SST anomalies and surface winds regressed onto the NA index by using observational data from 1970 to 2018 (Fig. 1b). The NA warming pattern resembles a developing positive phase of the Atlantic Multidecadal Oscillation (AMO), and strong warming is observed in the tropical IO. The Pacific SST anomalies are close to the negative phase of the inter-decadal Pacific Oscillation and resemble local SST trends (e.g. Figure 1a). Thus, the warming trend in the Pacific may be related to both global warming and NA warming. The warm SST anomalies over the NA generate ascending motions and strong upper-level divergent flows, which are connected to upper-level convergences in the equatorial central and eastern Pacific (Supplementary Fig. S1), where descending motions dominate. Double anomalous anticyclonic circulation occurs in the central North and South Pacific (Fig. 1b). Descending motions in the central Pacific induce easterly anomalies over the central to western Pacific. In addition, warmer SST in the western Pacific strengthens Walker circulation by enhancing ascending motion over the equatorial western Pacific. The anomalous easterlies over the western Pacific deepen the equatorial thermocline and increase the sea surface height (SSH) (Supplementary Figs. S2c and S2d), thereby increasing ocean heat advection22,27 to the eastern Indian Ocean through enhanced south China Sea and Indonesian throughflows (Supplementary Fig. S3a).
Similarly, the SST and surface wind anomalies associated with IO warming are linked to strong NA warming. The NA warming pattern resembled the AMO positive phase. IO warming induces relatively higher SSH and SST near South Africa (Supplementary Figs. S2a and S2b), thereby enhancing the transport of warmer ocean warmer to the tropical AO. The ascending motion also generates westerly anomalies in the tropical AO by strengthening the sinking motion over the eastern Pacific (Supplementary Fig. S1b), which may contribute to tropical Atlantic warming. The ascending motion in the IO generates a descending motion in the tropical Atlantic region, which affects precipitation in the tropical Atlantic region (Supplementary Fig. S1b). Because the simultaneous regression does not explain causality, we conducted numerical experiments to illuminate the cause and effects of IO and NA warming.
IO-NA warming chain simulated in historical experiments
To explore whether IO and NA warming affect each other, we conducted two groups of numerical experiments using an Earth system model (see the ‘Methods’ section). The first group represents a historical simulation with prescribed IO SST warming (observed pattern) trends (0.01–0.70°C per 50 yrs) (HIS_OBS_IO). The model was freely coupled in other areas, and external forcings were based on the CMIP6 historical simulation protocol. Similarly, we conducted model simulations for the other group by specifying different trends of observed SST anomalies over the NA region (0.01–1.02°C per 50 yrs) (HIS_OBS_NA).
The ensemble IO Index simulated by the observed SST anomalies over the NA region (HIS_OBS_NA) captures the observed trends reasonably well for most periods (Fig. 2a). The range of each member from the ensemble mean is ± 0.12°C (approximate biases of 10–15%), indicating that the impacts of internal variability on the simulated IO SST are moderate; thus, the response of the IO influenced by NA warming is realistic. We conducted the same historical simulation—except for the prescribed climatological SST of NA (HIS_CLIM_NA)—to eliminate the impact of NA warming on IO warming. The results show that the IO SST warming was reduced by 40–50% (green line in Fig. 2a), indicating that the observed IO warming arises from the remote effects of NA warming and local IO warming.
For another historical simulation with observed SST forcings over the IO region (Fig. 2b), the ensemble mean SST over the NA region shows realistic temporal evolution, with negative SST anomalies being observed from the 1970s to the 1990s along with rapidly increasing positive SST anomalies after the late 1990s. The ensemble mean SST is slightly less than that observed, which may be attributed to the systematic mean SST biases of the model south of Greenland34. The ensemble spread from each member is moderate (\(\pm\) 0.26\(℃\)), signifying that all ensemble simulations capture the NA warming trends induced by IO warming reasonably well. The same historical simulation—but with climatological SST nudged over the IO region—yields a 50% reduction in NA SST warming (green line, Fig. 2b). This result suggests that approximately half of NA warming may be attributed to IO warming.
Historical runs that nudged IO (or NA) SST warming (Fig. 2c and 2d) demonstrated that increasing IO SST could lead to a warmer NA; likewise, IO SST would increase with a warmer NA. The IO SST trend forced by NA SST is smaller than the NA SST trends (by approximately 40–50%). Thus, IO warming may be attributed to a remote NA impact and to local surface flux changes. We will further study this issue in the future. The linear relationship between IO and NA warming suggests that the IO-NA warming chain may occur without strong anthropogenic forcing.
IO-NA warming chain simulated without anthropogenic forcing effects
We hypothesised that the IO-NA warming chain could be generated without increasing the anthropogenic forcing. To verify this hypothesis, we conducted idealised numerical simulations with fixed CMIP6 preindustrial (PI) forcings to remove the warming effect caused due to increased CO2 forcing. Under fixed PI external forcings, we conducted a suite of 100-yr simulations by using observed IO SST anomalies averaged from 1980–2018 (PI_IO) and progressively increased uniform SST anomalies by 1°C (PI_IO + 1), 2°C (PI_IO + 2), and 3°C (PI_IO + 3). We repeated the same experiment with decreased uniform SST anomalies by 1°C (PI_IO-1), 2°C (PI_IO-2), and 3°C (PI_IO-2). Additionally, we conducted a similar suite of fixed-forcing PI experiments with the nudged SST over the NA region (PI_NA).
The SST anomalies induced by 1°C warming over the NA region captured the observed patterns (Fig. 1a) and showed significant IO warming (Fig. 4a). Warming also occurred in the northern and southern subtropical Pacific and equatorial western Pacific; however, relatively cool SST prevailed in the equatorial eastern Pacific. The anticyclonic flows in the north and south Pacific generate easterly anomalies in the western Pacific that increase ocean heat transport from the western Pacific to the eastern Indian Ocean. SST and surface wind changes induced by 1°C IO warming (Fig. S4a) largely resemble those in the observations—except for the Indian Ocean surface winds (Fig. 1c). The above-mentioned simulation results suggest that the IO-NA warming chain can occur without the global warming effect.
We further examined how the NA warming trends were linked to the prescribed IO warming under PI forcings (Fig. 3a). When the IO warming increases, NA SST is progressively warmer. The increasing rate of NA SST per 1°C IO warming is approximately 0.39°C. The correlation coefficient between relative IO warming and the NA response is 0.98, implying that the response of NA by IO forcings is extremely linear from the model simulations. The simulation with NA warming forcings yields a rate of increase of IO SST per 1°C NA warming of 0.33°C; the correlation between IO and NA SST changes is 0.97. The results show that the IO-NA warming chain mechanism could occur without greenhouse gas forcing. The historical simulation shows much higher trends of NA warming by IO warming forcing than those observed using PI simulation. Thus, anthropogenic forcing may intensify the increasing trends of NA warming induced by IO warming.
Change in IO-NA warming chain under anthropogenic forcing
We conducted a suite of 100-yr simulations under quadruple CO2 external forcings by using prescribed observed SST anomalies that were averaged from 1980–2018 in the Indian Ocean (4CO2_IO); subsequently, we progressively increased uniform SST by 1°C, 2°C, and 3°C to determine how the IO-NA warming chain would change with continued anthropogenic forcing. We repeated the same experiment but with decreased uniform SST anomalies at 1°C, 2°C, and 3°C. Additionally, we repeated 4CO2_IO experiments with observed SST anomalies over the NA region (4CO2_NA).
The NA response to IO warming forcing was extremely linear in the 4CO2-IO experiment. The rate of NA SST increase with 1°C of IO warming is approximately 0.56°C, which is 50% higher than that observed in the PI_IO (Fig. 3). Under strong anthropogenic forcing, the mean SST warmed considerably in South Africa and the South Indian Ocean, thereby increasing ocean heat transport via the Agulhas Current (Supplementary Fig. S3b). Similarly, the change in IO SST due to NA warming remains extremely linear, with a correlation of r = 0.98. The rate of increase in IO SST due to NA (0.39°C) warming under strong anthropogenic forcing was slightly higher than that of the PI (0.33°C). Under enhanced greenhouse gas forcing, the eastern Pacific is warmer than the western Pacific. NA warming generates strong anomalous easterlies in the central Pacific with the movement of warm seawater to the western Pacific, resulting in a warmer SST, a higher SSH, and more IO warming (Supplementary Figs. S2c and S2d). These results show that the IO-NA warming chain is intensified by anthropogenic forcing.
Reduction of North American precipitation by the IO-NA warming chain
The intensified IO-NA warming chain due to global warming may affect IO-induced (NA-induced) precipitation changes in the mid-latitudes. In the PI_IO run, precipitation in western North America due to IO warming of 1°C is reduced by 7–18% (Fig. 3c). IO warming generates anticyclonic circulation in the North Pacific, thereby transporting cold water around the Bering Sea to the west of North America. The NA warming induced by IO warming generates strong sinking motions in the western U.S. region, thereby inducing less precipitation. Under anthropogenic forcing, the decreasing rate of precipitation change due to IO warming is approximately 15–47% (average of approximately 30%), which is approximately double the corresponding PI run. Similarly, the warm NA SST caused by IO warming strengthened the descending motion and further reduced precipitation in the western U.S. by 10–35%.