Radiocarbon dating of well-preserved, in-place vegetation exposed by the retreating Quelccaya Ice Cap of southeastern Peru constrains the last time the ice cap's extent was smaller than at present. Seventeen plant samples from two sites along the central western margin collectively date to 4700 and 5100 cal yr BP and strongly indicate that current ice cap retreat is unprecedented over the past ∼ 5 millennia. Seventeen vegetation samples interbedded in a nearby clastic sedimentary sequence suggest ice-free conditions at this site from ∼ 5200 to at least ∼ 7000 cal yr BP, and place minimum constraint on early- to mid-Holocene ice cap extent.

Moraine chronology is combined with digital topography to model deglacial rates of paleoglacier volumes in both the Huancané Valley on the west side of the Quelccaya Ice Cap and the Upismayo Valley on the northwest side of the Cordillera Vilcanota. The fastest rates of deglaciation (39×10−5 to 114×10−5 km3 yr−1 and 112×10−5 to 247×10−5 km3 yr−1 for each valley, respectively) were calculated for the most recent paleoglaciers, corresponding to the last few centuries. These results are consistent with observations in the Venezuelan Andes showing high rates of deglaciation since the Little Ice Age. These rates also fall within the range of 20th century rates of deglaciation measured on the Quelccaya Ice Cap (29×10−5 to 220×10−5 km3 yr−1, Brecher and Thompson, 1993; Thompson, 2000). These results imply that rates of deglaciation may fluctuate significantly over time and that high rates of deglaciation may not be exclusive to the late 20th century. Equilibrium line altitude (ELA) depressions for the ice volumes of the last glaciation modeled here were computed as 230 m for the Quelccaya Ice Cap and 170 m for the Cordillera Vilcanota. Maximum ELA depressions are lower than previously published: <500 m for the Cordillera Vilcanota and <400 m for the Quelccaya Ice Cap. These lower values could imply a topographic control over paleoglacier extent.

In this paper, we reanalyze the glacier mass balance on Glaciar Antisana 15α over the 1995–2012 period. Annual glacier mass balances were quantified on the basis of monthly glaciological measurements using an adaptation of Lliboutry's statistical approach. The geodetic mass balance was computed between 1997 and 2009 giving a cumulative balance of −1.39 ± 1.97 m w.e. and a slightly negative adjusted annual glaciological mass balance (−0.12 ± 0.16 m w.e. a−1). Despite a careful analysis of uncertainties, we found a large discrepancy between the cumulative glaciological and the geodetic mass balances over the common period, of 4.66 m w.e. This discrepancy can mainly be explained by underestimated net accumulation in the glacier upper reaches, which could be due to the peculiar climate conditions of the equatorial zone with year round accumulation, thereby preventing clear identification of annual layers. An increase of ~70% in measured rates of net accumulation would be needed to balance the glaciological and geodetic mass balances; a hypothesis confirmed by estimated ice flux in the vicinity of the ELA. Consequently, the vertical gradient of precipitation may be higher than previously estimated and the accumulation processes (including the role of frost deposition) need to be carefully analyzed.