The notion that the Atlantic Meridional Overturning Circulation (AMOC) can have more than one stable equilibrium emerged in the 1980s as a powerful hypothesis to explain rapid climate variability during the Pleistocene. Ever since, the idea that a temporary perturbation of the AMOC—or a permanent change in its forcing—could trigger an irreversible collapse has remained a reason for concern. Here we review literature on the equilibrium stability of the AMOC and present a synthesis that puts our understanding of past and future AMOC behavior in a unifying framework. This framework is based on concepts from Dynamical Systems Theory, which has proven to be an important tool in interpreting a wide range of model behavior. We conclude that it cannot be ruled out that the AMOC in our current climate is in, or close to, a regime of multiple equilibria. But there is considerable uncertainty in the location of stability thresholds with respect to our current climate state, so we have no credible indications of where our present‐day AMOC is located with respect to thresholds. We conclude by identifying gaps in our knowledge and proposing possible ways forward to address these gaps.

Because of the pervasive role of water in the Earth system, the relative abundances of stable isotopologues of water are valuable for understanding atmospheric, oceanic, and biospheric processes, and for interpreting paleoclimate proxy reconstructions. Isotopologues are transported by both large‐scale and turbulent flows, and the ratio of heavy to light isotopologues changes due to fractionation that can accompany condensation and evaporation processes. Correctly predicting the isotopic distributions requires resolving the relationships between large‐scale ocean and atmospheric circulation and smaller‐scale hydrological processes, which can be accomplished within a coupled climate modeling framework. Here we present the water isotope‐enabled version of the Community Earth System Model version 1 (iCESM1), which simulates global variations in water isotopic ratios in the atmosphere, land, ocean, and sea ice. In a transient Last Millennium simulation covering the 850–2005 period, iCESM1 correctly captures the late‐twentieth‐century structure of δ¹⁸O and δD over the global oceans, with more limited accuracy over land. The relationship between salinity and seawater δ¹⁸O is also well represented over the observational period, including interbasin variations. We illustrate the utility of coupled, isotope‐enabled simulations using both Last Millennium simulations and freshwater hosing experiments with iCESM1. Closing the isotopic mass balance between all components of the coupled model provides new confidence in the underlying depiction of the water cycle in CESM, while also highlighting areas where the underlying hydrologic balance can be improved. The iCESM1 is poised to be a vital community resource for ongoing model development with both modern and paleoclimate applications.

Oxygen isotope speleothem records exhibit coherent variability over the pan-Asian summer monsoon (AM) region. The hydroclimatic representation of these oxygen isotope records for the AM, however, has remained poorly understood. Here, combining an isotope-enabled Earth system model in transient experiments with proxy records, we show that the widespread AM δ¹⁸O_c signal during the last deglaciation (20 to 11 thousand years ago) is accompanied by a continental-scale, coherent hydroclimate footprint, with spatially opposite signs in rainfall. This footprint is generated as a dynamically coherent response of the AM system primarily to meltwater forcing and secondarily to insolation forcing and is further reinforced by atmospheric teleconnection. Hence, widespread δ¹⁸O_p depletion in the AM region is accompanied by a northward migration of the westerly jet and enhanced southwesterly monsoon wind, as well as increased rainfall from South Asia (India) to northern China but decreased rainfall in southeast China.

Studying the El Niño–Southern Oscillation (ENSO) in the past can help us better understand its dynamics and improve its future projections. However, both paleoclimate reconstructions and model simulations of ENSO strength at the Last Glacial Maximum (LGM; 21 ka B.P.) have led to contradicting results. Here we perform model simulations using the recently developed water isotope‐enabled Community Earth System Model (iCESM). For the first time, model‐simulated oxygen isotopes are directly compared with those from ENSO reconstructions using the individual foraminifera analysis (IFA). We find that the LGM ENSO is most likely weaker comparing with the preindustrial. The iCESM suggests that total variance of the IFA records may only reflect changes in the annual cycle instead of ENSO variability as previously assumed. Furthermore, the interpretation of subsurface IFA records can be substantially complicated by the habitat depth of thermocline‐dwelling foraminifera and their vertical migration with a temporally varying thermocline.

The Beaufort Gyre (BG), the largest Arctic Ocean freshwater reservoir, has drastically increased its liquid freshwater content by 40% in the past two decades. If released within a short period, the excess freshwater could potentially impact the large-scale ocean circulation by freshening the upper subpolar North Atlantic. Here, we track BG-sourced freshwater using passive tracers in a global ocean sea-ice model and show that this freshwater exited the Arctic mostly through the Canadian Arctic Archipelago, rather than Fram Strait, during an historical release event in 1983–1995. The Labrador Sea is the most affected region in the subpolar North Atlantic, with a freshening of 0.2 psu on the western shelves and 0.4 psu in the Labrador Current. Given that the present BG freshwater content anomaly is twice the historical analog studied here, the impact of a future rapid release on Labrador Sea salinity could be significant, easily exceeding similar fluxes from Greenland meltwater.

Greenland ice-core δ ¹⁸ O-temperature reconstructions suggest a dramatic cooling during the Younger Dryas (YD; 12.9–11.7 ka), with temperatures being as cold as the earlier Oldest Dryas (OD; 18.0–14.6 ka) despite an approximately 50 ppm rise in atmospheric CO ₂ . Such YD cooling implies a muted Greenland climate response to atmospheric CO ₂ , contrary to physical predictions of an enhanced high-latitude response to future increases in CO ₂ . Here we show that North Atlantic sea surface temperature reconstructions as well as transient climate model simulations suggest that the YD over Greenland should be substantially warmer than the OD by approximately 5 °C in response to increased atmospheric CO ₂ . Additional experiments with an isotope-enabled model suggest that the apparent YD temperature reconstruction derived from the ice-core δ ¹⁸ O record is likely an artifact of an altered temperature-δ ¹⁸ O relationship due to changing deglacial atmospheric circulation. Our results thus suggest that Greenland climate was warmer during the YD relative to the OD in response to rising atmospheric CO ₂ , consistent with sea surface temperature reconstructions and physical predictions, and has a sensitivity approximately twice that found in climate models for current climate due to an enhanced albedo feedback during the last deglaciation.