Glacier equilibrium line altitude variations during the “Little Ice Age” in the Mediterranean Andes (30◦n37◦tS)

Abstract. The "Little Ice Age" (LIA; 1500–1850 Common Era (CE)), has long been recognized as the last period when mountain glaciers in many regions of the Northern Hemisphere (NH) recorded extensive growth intervals in terms of their ice mass and frontal position. The knowledge about this relevant paleoclimatic interval is vast in mountainous regions such as the Alps and Rocky Mountains in North America. However, in extra-tropical Andean sub-regions such as the Mediterranean Andes of Chile and Argentina (MA; 30o–37o S), the LIA has been poorly documented. Paradoxically, the few climate reconstructions performed in the MA based on lake sediments and tree rings do not show clear evidence of a LIA climate anomaly as observed in the NH. In addition, recent studies have demonstrated temporal differences between mean air temperature variations across the last millennium between both hemispheres. This motivates our hypothesis that the LIA period was not associated with a significant climate perturbation in the MA region. Considering this background, we performed an experiment using daily climatic variables from three Global Climate Models (GCMs) to force a novel glaciological model. In this way, we simulated temporal variations of the glacier equilibrium-line altitude (ELA) to evaluate the glacier response during the period 1500–1848 CE. Overall, each GCM shows temporal changes in annual ELA, with anomalously low elevations during 1640–1670 and 1800–1848 CE. An interval with high ELA values was identified during 1550–1575 CE. The spectral properties of the mean annual ELA in each GCM present significant periodicities between 2–7 years, and also significant decadal to multi-decadal signals. In addition, significant and coherent cycles at interannual to multi-decadal scales were detected between modeled mean annual ELAs and the first EOF1 extracted from Sea Surface Temperature (SST) within the El Nino 3.4 of each GCM. Finally, significant Pearson correlation coefficients were obtained between the mean annual ELA and Pacific SST on interannual to multi-decadal timescales. According to our findings, we propose that Pacific SST variability was the main modulator of temporal changes of the ELA in the MA region of South America during 1500–1848 CE.