The Younger Dryas stadial, a cold event spanning 12,800 to 11,500 years ago, during the last deglaciation, is thought to coincide with the last major glacial re-advance in the tropical Andes. This interpretation relies mainly on cosmic-ray exposure dating of glacial deposits. Recent studies, however, have established new production rates for cosmogenic (10)Be and (3)He, which make it necessary to update all chronologies in this region and revise our understanding of cryospheric responses to climate variability. Here we present a new (10)Be moraine chronology in Colombia showing that glaciers in the northern tropical Andes expanded to a larger extent during the Antarctic cold reversal (14,500 to 12,900 years ago) than during the Younger Dryas. On the basis of a homogenized chronology of all (10)Be and (3)He moraine ages across the tropical Andes, we show that this behaviour was common to the northern and southern tropical Andes. Transient simulations with a coupled global climate model suggest that the common glacier behaviour was the result of Atlantic meridional overturning circulation variability superimposed on a deglacial increase in the atmospheric carbon dioxide concentration. During the Antarctic cold reversal, glaciers advanced primarily in response to cold sea surface temperatures over much of the Southern Hemisphere. During the Younger Dryas, however, northern tropical Andes glaciers retreated owing to abrupt regional warming in response to reduced precipitation and land-surface feedbacks triggered by a weakened Atlantic meridional overturning circulation. Conversely, glacier retreat during the Younger Dryas in the southern tropical Andes occurred as a result of progressive warming, probably influenced by an increase in atmospheric carbon dioxide. Considered with evidence from mid-latitude Andean glaciers, our results argue for a common glacier response to cold conditions in the Antarctic cold reversal exceeding that of the Younger Dryas.

Abstract. Volcanic aerosols exert the most important natural radiative forcing of the last millennium. State-of-the-art paleoclimate simulations of this interval are typically forced with diverse spatial patterns of volcanic forcing, leading to different responses in tropical hydroclimate. Recently, theoretical considerations relating the intertropical convergence zone (ITCZ) position to the demands of global energy balance have emerged in the literature, allowing for a connection to be made between the paleoclimate simulations and recent developments in the understanding of ITCZ dynamics. These energetic considerations aid in explaining the well-known historical, paleoclimatic, and modeling evidence that the ITCZ migrates away from the hemisphere that is energetically deficient in response to asymmetric forcing.Here we use two separate general circulation model (GCM) suites of experiments for the last millennium to relate the ITCZ position to asymmetries in prescribed volcanic sulfate aerosols in the stratosphere and related asymmetric radiative forcing. We discuss the ITCZ shift in the context of atmospheric energetics and discuss the ramifications of transient ITCZ migrations for other sensitive indicators of changes in the tropical hydrologic cycle, including global streamflow. For the first time, we also offer insight into the large-scale fingerprint of water isotopologues in precipitation (δ18Op) in response to asymmetries in radiative forcing. The ITCZ shifts away from the hemisphere with greater volcanic forcing. Since the isotopic composition of precipitation in the ITCZ is relatively depleted compared to areas outside this zone, this meridional precipitation migration results in a large-scale enrichment (depletion) in the isotopic composition of tropical precipitation in regions the ITCZ moves away from (toward). Our results highlight the need for careful consideration of the spatial structure of volcanic forcing for interpreting volcanic signals in proxy records and therefore in evaluating the skill of Common Era climate model output.

Abstract. Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-Intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El Niño-Southern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850 CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium.An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records.Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the "amount effect". During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

How does climate sensitivity vary with the magnitude of climate forcing? This question was investigated with the use of a modified coupled atmosphere‐ocean model, whose stability was improved so that the model would accommodate large radiative forcings yet be fast enough to reach rapid equilibrium. Experiments were performed in which atmospheric CO2 was multiplied by powers of 2, from 1/64 to 256 times the 1950 value. From 8 to 32 times, the 1950 CO2, climate sensitivity for doubling CO2 reaches 8°C due to increases in water vapor absorption and cloud top height and to reductions in low level cloud cover. As CO2 amount increases further, sensitivity drops as cloud cover and planetary albedo stabilize. No water vapor‐induced runaway greenhouse caused by increased CO2 was found for the range of CO2 examined. With CO2 at or below 1/8 of the 1950 value, runaway sea ice does occur as the planet cascades to a snowball Earth climate with fully ice covered oceans and global mean surface temperatures near −30°C.