Deciphering the climate changes that influenced the glacial fluctuations of the last millennium requires documenting the spatial and temporal patterns of these glacial events. Here, we estimate the change in equilibrium line altitudes (ELAs) between the most prominent glacial advance of the last millennium and the present for four alpine glaciers located in different climatic regimes along the Andes. For each glacier, we reconstruct scenarios of climatic conditions (temperature and precipitation anomalies) that accommodate the observed ELA changes. We focus on the following glaciers: an alpine glacier in the Cordillera Vilcanota (13°S), Tapado glacier (30°S), Cipreses glacier (34°S), and Tranquilo glacier (47°S). Our results show that the range of possible temperature and precipitation anomalies that accommodate the observed ELA changes overlap significantly at three of the four sites (i.e. Vilcanota, Cipreses, and Tranquilo). Only Tapado glacier exhibits a set of climate anomalies that differs from the other three sites. Assuming no change in precipitation, the estimated ELA changes require a cooling of at least 0.7°C in the Cordillera Vilcanota, 1.0°C at Tapado glacier, 0.6°C at Cipreses glacier, and 0.7°C at Tranquilo glacier. Conversely, assuming no change in temperature, the estimated ELA changes are explained by increases in precipitation exceeding 0.52 m yr−1 (64% of the annual precipitation) in the Cordillera Vilcanota, 0.31 m yr−1 (89%) at Tapado glacier, 0.22 m yr−1 (27%) at Cipreses glacier, and 0.3 m yr−1 (27%) at Tranquilo glacier. By mapping the ELA changes and modeling the potential climate forcing across diverse climate settings, we aim to contribute toward documenting the spatial variability of climate conditions during the last millennium, a key step to decipher the mechanisms underlying the glacial fluctuation that occurred during this period.

a b s t r a c tImproving the late Quaternary paleoclimate record through climate interpretations of low-latitude glacier length changes advances our understanding of past climate change events and the mechanisms for past, present, and future climate change. Paleotemperature reconstructions at low-latitude glaciers are uniquely fruitful because they can provide both site-specific information and enhanced understanding of regional-scale variations due to the structure of the tropical atmosphere. We produce Little Ice Age (LIA) and Younger Dryas (YD) paleoclimate reconstructions for the Huancan e outlet glacier of the Quelccaya Ice Cap (QIC) and low-latitude southern hemisphere regional sea surface temperatures (SSTs) using a coupled ice-flow and energy balance model. We also model the effects of long-term changes in the summit temperature and precipitiation rate and the effects of interannual climate variability on the Huancan e glacier length. We find temperature to be the dominant climate driver of glacier length change. Also, we find that interannual climate variability cannot adequately explain glacier advances inferred from the geomorphic record, necessitating that these features were formed during past colder climates. To constrain our LIA reconstruction, we incorporate the QIC ice core record, finding a LIA air temperature cooling at the ice cap of between~0.7 C and~1.1 C and~0.4 C and regional SSTs cooling of~0.6 C. For the YD paleoclimate reconstructions, we propose two limits on the precipitation rate, since the ice core record does not extend into the Pleistocene: 1) the precipitation rate scales with the Clausius-Clapeyron relationship (upper limit on cooling) and 2) the precipitation rate increases by 40% (lower limit on cooling), which is an increase about twice as great as the regional increases realized in GCM simulations for the period. The first limit requires~1.6 C cooling in ice cap air temperatures and 0.9 C cooling in SSTs, and the second limit requires~1.0 C cooling in ice cap air temperatures and 0.5 C cooling in SSTs. Our temperature reconstructions are in good agreement with the magnitude and trend of GCM simulations that incorporate the forcing mechanisms hypothesized to have caused these climate change events.

The Laurentide Ice Sheet (LIS) was the largest ice sheet during the last glacial period. An accurate representation of its behavior during the last deglaciation is critical to understanding its influence on and response to a changing climate. We use 10Be dating and Bayesian modeling to track the recession of the southwest sector of the Labrador Dome of the LIS along an ~500-km-long transect west of Lake Superior during the last deglaciation. This transect reflects terrestrial ice-margin retreat and crosses multiple moraine sets, with the southwestern part of the transect deglaciated by ca. 19 ka and the northeastern part deglaciated by ca. 10 ka. The predominant behavior of the ice margin during this interval is near-constant retreat with retreat rates varying between ~59 m/a and 38 m/a. The moraine sets mark standstills and/or readvances that in total constitute only ~17% of the retreat interval. The spatial and temporal pattern of ice-margin retreat tracked here differs from existing reconstructions that are based on using isochrons to define ice-margin positions. Acknowledging the uncertainties associated with the modeled ages of ice-margin retreat, we suggest that the overall retreat pattern is consistent with forcing by a gradual increase in Northern Hemisphere, high-latitude summer insolation. The pattern of ice-margin retreat is inconsistent with Greenland ice-core temperature records, and thus these records may not be suitable to drive models of the LIS.

Changes in glacier length and extent are indicators of contemporary and archives of past climate changes, but this common climate proxy presents a challenge for inferring a climate signal. Modeling studies suggest that length fluctuations can occur due to interannual climate variability within an unchanging mean climate and that changes in interannual climate variability can also drive changes in average length. This paper quantifies the impacts of interannual climate variability on average glacier length and mass balance, using a flowline model coupled to a simplified mass-balance model. Results illustrate that changes in the magnitude of interannual temperature variability can non-linearly affect the mean glacier length through a mass-balance asymmetry between warm and cold years. This asymmetry is present in models where melt only initiates after a temperature threshold is crossed. Glaciers susceptible to this asymmetry can be identified based on the shape of their mass-balance profiles. The presence of mass-balance asymmetries in glaciological databases is evaluated, but current records are too short for high statistical resolving power. While the asymmetry in this study can affect the average length and mass-balance, its impacts are small, and paleoclimate interpretations from glacier-length changes are likely not notably influenced by this process.

As tropical glaciers continue to retreat, we need accurate knowledge about where they are located, how large they are, and their retreat rates. Remote sensing data are invaluable for tracking these hard-to-reach glaciers. However, remotely identifying tropical glaciers is prone to misclassification errors due to ephemeral snow cover. We reevaluate the size and retreat rates of the two largest tropical ice masses, the Quelccaya Ice Cap (Peru) and Nevado Coropuna (Peru), using remote sensing data from the Landsat missions. To quantify their glacial extents more accurately, we expand the time window for our analysis beyond the dry season (austral winter), processing in total 529 Landsat scenes. We find that Landsat scenes from October, November, and December, which are after the dry season, better capture the glacial extent since ephemeral snow cover is minimized. We compare our findings to past studies of tropical glaciers, which have mainly analyzed scenes from the dry season. Our reevaluation finds that both tropical ice masses are smaller but retreating less rapidly than commonly reported. These findings have implications for these ice masses as sustained water resources for downstream communities.