Can multi-glacier moraine analysis constrain paleoclimate reconstructions? Evidence from 3D glacier modeling in the Cordillera Real, Bolivia
- 1Aix-Marseille Université (Marseille, France)
- 2Centre de Recherche et d’Enseignement de Géosciences de l’Environnement (Aix-en-Provence, France)
- 3Centre National de la Recherche Scientifique (Paris, France)
Abstract
Glacier moraines are widely used as proxies for past climate, as they are interpreted to reflect near-equilibrium glacier extents. However, a single glacier geometry can arise from multiple temperature-precipitation combinations and their quantification relies on glacier modeling that introduce additional uncertainties. Here, we investigate how this non-uniqueness can be constrained from a 3D full-Stokes ice-flow model (Elmer/Ice), applied to neighboring coeval glaciers in the Zongo-Charquini area of the Bolivian Cordillera Real. Using synthetic experiments, in which Temperature-Precipitation ( T - P ) anomaly curves are derived from several glacier extents generated under prescribed climate conditions, we show that differences in glacier geometry and hypsometry cannot be used to identify unique temperature and precipitation conditions. However, we also show that glacier thickness may provide a way to dissociate temperature and precipitation but with a limited sensitivity that would require highly accurate past volume estimates to be useful in practice, which is often not feasible. Using real moraine records, our results show that single-glacier reconstructions are strongly influenced by both, methodological choice (melt model), and site-specific uncertainties, resulting in a range of T - P anomaly curves describing the same climate. In this context, using a multi-glacier framework allows to better assess the temperature-precipitation condition and its associated uncertainties. Using this approach combined with independent tree-ring and hydrological lake-balance precipitation reconstruction in the region, we estimate temperature anomalies of −0.8 ± 0.3 °C for the Little Ice Age and −1.3 ± 0.4 °C for the Early Holocene (∼ 10 ka BP), relative to a present-day reference climate (1950–2023).
Versions
Related
- Inside the Extinction Envelope: A Systemic Review and Meta-Analysis of Paleoclimate Extinction Markers with Comments on the Necessary Array of Last-Threshold Interventions
- Coupled multiscale paleoclimate reconstruction with four-dimensional variational data assimilation
- Dating of Late Pleistocene sediments along the coast of Brittany (W France) by Optically Stimulated Luminescence, paleoclimate implications and Eemian highstand
- Volcanic controls on smectite authigenesis and paleoclimate signals across the Paleocene-Eocene Thermal Maximum (PETM) sedimentary record from the Norwegian-Greenland Sea
- Fault‐Related Carbonates as Quaternary Seismotectonic and Paleoclimate Archives (Monte Marine Fault, Central Apennines)
BibTeX
@article{reyes2026multiglacier,
title = {{Can multi-glacier moraine analysis constrain paleoclimate reconstructions? Evidence from 3D glacier modeling in the Cordillera Real, Bolivia}},
author = {Nicolás Acuña Reyes and Léo Martin and Adrien Gilbert and Vincent Jomelli and Antoine Rabatel and Simon Filhol and Deborah Verfaillie and Pierre-Henri Blard and Jérôme Lavé},
journal = {Climate of the Past},
year = {2026},
month = nov,
doi = {10.5194/cp-22-1941-2026},
url = {https://cp.copernicus.org/articles/22/1941/2026/},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Can multi-glacier moraine analysis constrain paleoclimate reconstructions? Evidence from 3D glacier modeling in the Cordillera Real, Bolivia
AU - Nicolás Acuña Reyes
AU - Léo Martin
AU - Adrien Gilbert
AU - Vincent Jomelli
AU - Antoine Rabatel
AU - Simon Filhol
AU - Deborah Verfaillie
AU - Pierre-Henri Blard
AU - Jérôme Lavé
PY - 2026
DA - 2026/11/05
JO - Climate of the Past
DO - 10.5194/cp-22-1941-2026
UR - https://cp.copernicus.org/articles/22/1941/2026/
AB - Glacier moraines are widely used as proxies for past climate, as they are interpreted to reflect near-equilibrium glacier extents. However, a single glacier geometry can arise from multiple temperature-precipitation combinations and their quantification relies on glacier modeling that introduce additional uncertainties. Here, we investigate how this non-uniqueness can be constrained from a 3D full-Stokes ice-flow model (Elmer/Ice), applied to neighboring coeval glaciers in the Zongo-Charquini area of the Bolivian Cordillera Real. Using synthetic experiments, in which Temperature-Precipitation ( T - P ) anomaly curves are derived from several glacier extents generated under prescribed climate conditions, we show that differences in glacier geometry and hypsometry cannot be used to identify unique temperature and precipitation conditions. However, we also show that glacier thickness may provide a way to dissociate temperature and precipitation but with a limited sensitivity that would require highly accurate past volume estimates to be useful in practice, which is often not feasible. Using real moraine records, our results show that single-glacier reconstructions are strongly influenced by both, methodological choice (melt model), and site-specific uncertainties, resulting in a range of T - P anomaly curves describing the same climate. In this context, using a multi-glacier framework allows to better assess the temperature-precipitation condition and its associated uncertainties. Using this approach combined with independent tree-ring and hydrological lake-balance precipitation reconstruction in the region, we estimate temperature anomalies of −0.8 ± 0.3 °C for the Little Ice Age and −1.3 ± 0.4 °C for the Early Holocene (∼ 10 ka BP), relative to a present-day reference climate (1950–2023).
SN - 1814-9332
ER -