skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on January 1, 2027

Title: A 3He-based Holocene glacial chronology from Villarrica volcano, Chile
Understanding alpine glacier extent during past climate variability is instructive for determining the glacier response to future climate change. Villarrica volcano is a late Pleistocene stratovolcano located in Chile's Southern Volcanic Zone that was covered by the Patagonian Ice Sheet during the last glacial period, and still retains small remnant glaciers today. Moraines preserved several kilometers from the summit on different flanks of the volcano record a history of expanded glacier lengths during the Holocene. However, the precise ages of these moraines are unknown, and the Holocene glacial history of Villarrica remains poorly constrained, limiting our understanding of how glaciers in this region responded to Holocene climate change. To constrain the timing of these moraines, we analyzed cosmogenic 3He in olivine from 25 basaltic andesite moraine boulders for cosmogenic surface exposure dating. Our new chronology reveals multiple late Holocene glacier advances from different flanks of the volcano, with the glaciers culminating and abandoning their moraines during the early Neoglacial period at ∼3355 ± 190 a and ∼1735 ± 215 a, and during the last millennium spanning the Little Ice Age period at ∼720 ± 225 a, ∼370 ± 75 a, and in the last ∼200 years. Our analysis of Holocene climate proxies from south-central Chile indicates that the early Neoglacial advances and subsequent retreat likely reflect increased effective moisture delivered by intensified Southern Westerly Winds and associated shifts in their latitudinal position. In contrast, we interpret the last millennium glacier advances as primarily driven by reduced summer ablation linked to regional cooling, followed by glacier retreat due to increased temperatures. Our chronology and closely spaced moraine positions suggest that glacier retreat on Villarrica, and possibly the broader Southern Volcanic Zone, has been gradual during the late Holocene and interrupted by short-lived advances driven by varying changes in temperature and moisture.  more » « less
Award ID(s):
2121372
PAR ID:
10650766
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Quaternary Science Reviews
Volume:
372
Issue:
C
ISSN:
0277-3791
Page Range / eLocation ID:
109707
Subject(s) / Keyword(s):
Villarrica volcano Southern Volcanic Zone (SVZ) Patagonian Ice Sheet (PIS) Late Holocene Cosmogenic isotopes Glacial chronology Southern Hemisphere Paleoclimate
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 ∘C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments. 
    more » « less
  2. The causes underlying Holocene glacier fluctuations remain elusive, despite decades of research efforts. Cosmogenic nuclide dating has allowed systematic study and thus improved knowledge of glacier-climate dynamics during this time frame, in part by filling in geographical gaps in both hemispheres. Here we present a new comprehensive Holocene moraine chronology from Mt. San Lorenzo (47°S) in central Patagonia, Southern Hemisphere. Twenty-four new 10 Be ages, together with three published ages, indicate that the Río Tranquilo glacier approached its Holocene maximum position sometime, or possibly on multiple occasions, between 9,860 ± 180 and 6,730 ± 130 years. This event(s) was followed by a sequence of slightly smaller advances at 5,750 ± 220, 4,290 ± 100 (?), 3,490 ± 140, 1,440 ± 60, between 670 ± 20 and 430 ± 20, and at 390 ± 10 years ago. The Tranquilo record documents centennial to millennial-scale glacier advances throughout the Holocene, and is consistent with recent glacier chronologies from central and southern Patagonia. This pattern correlates well with that of multiple moraine-building events with slightly decreasing net extent, as is observed at other sites in the Southern Hemisphere (i.e., Patagonia, New Zealand and Antarctic Peninsula) throughout the early, middle and late Holocene. This is in stark contrast to the typical Holocene mountain glacier pattern in the Northern Hemisphere, as documented in the European Alps, Scandinavia and Canada, where small glaciers in the early-to-mid Holocene gave way to more-extensive glacier advances during the late Holocene, culminating in the Little Ice Age expansion. We posit that this past asymmetry between the Southern and Northern hemisphere glacier patterns is due to natural forcing that has been recently overwhelmed by anthropogenic greenhouse gas driven warming, which is causing interhemispherically synchronized glacier retreat unprecedented during the Holocene. 
    more » « less
  3. Waitt RB, Thackray GD (Ed.)
    Mountain glacier moraine sequences and their chronologies allow us to evaluate the timing and climate conditions that underpin changes in the equilibrium line altitudes (ELAs), which can provide valuable information on the paleoclimatology of understudied regions such as tropical East Africa. However, moraine sequences are inherently discontinuous, and the precise climate conditions that they represent can be ambiguous due to the sensitivity of mountain glaciers to temperature, precipitation, and other environmental variables. Here, we used a two-dimensional (2-D) iceflow and mass-balance model to simulate glacier extents and ELAs in the Rwenzori Mountains in East Africa over the past 31,000 yr (31 k.y.), including the Last Glacial Maximum (LGM), late glacial period, and the Holocene Epoch. We drove the glacier model with two independent, continuous temperature reconstructions to simulate possible glacier length changes through time. Model input paleoclimate values came from branched glycerol dialkyl glycerol tetraether (brGDGT) temperature reconstructions from alpine lakes on Mount Kenya for the last ~31 k.y., and precipitation reconstructions for the LGM came from various East African locations. We then compared the simulated fluctuations with the positions and ages (where known) of the Rwenzori moraines. The simulated glacier extents reached within 1.1 km of the dated LGM moraines in one valley (93% of the full LGM extent) when forced by the brGDGT temperature reconstructions (maximum cooling of 6.1 °C) and a decrease in precipitation (-10% than modern amounts). These simulations suggest that the Rwenzori glaciers required a cooling of at least 6.1 °C to reach the dated LGM moraines. Based on the model output, we predict an age of 12–11 ka for moraines located halfway between the LGM and modern glacier extents. We also predict ice-free conditions in the Rwenzori Mountains for most of the early to middle Holocene, followed by a late Holocene glacier readvance within the last 2000 yr. 
    more » « less
  4. The reconstruction of former mountain glaciers from geomorphic mapping and cosmogenic-nuclide surface exposure dating provides a unique opportunity to infer patterns of past terrestrial climate variability. Tropical mountain glaciers are particularly valuable as there are comparatively few terrestrial climate proxies at equatorial latitudes relative to higher latitudes. As the single largest climate zone on Earth, the tropics play an outsized role in mediating global climate via the ocean-atmosphere transfer of latent heat and water vapor. Nonetheless, there remains a persistent gap in our understanding of how the tropics influenced – or were influenced by – the high-magnitude climate shifts of the Late Pleistocene, and whether this high-energy region simply responded to extratropical forcing or was itself a driver of global climatic change. To help address this knowledge gap, we analyzed geologic evidence for past glacial fluctuations in three adjacent valleys in the Sierra Nevada del Cocuy, the highest subrange of the Eastern Cordillera in the Colombian Andes, to provide a terrestrial record of atmospheric temperature during the latter part of Termination 1. Coupled with geomorphic mapping and paleo-snowline reconstructions, our beryllium-10 glacial chronology indicates that glaciers in the humid inner tropics underwent pronounced growth and gradual decay during the Antarctic Cold Reversal (14.5–12.8 ka) and Younger Dryas (12.8–11.7 ka) periods, respectively, following a trend that, according to directly dated moraine records from throughout both polar hemispheres, appears to have been global. While the specific mechanism(s) behind this large-scale behavior remains to be corroborated, we revisit the hypothesis that ocean atmosphere heat transfer and water vapor flux are key drivers of abrupt Lateglacial temperature fluctuations. Subsequent to the Lateglacial, deglaciation of the Sierra Nevada del Cocuy accelerated during the Early Holocene, a pattern also observed in other tropical glacier records. More recently, the magnitude of snowline rise and glacier retreat over the last two centuries supports the view that modern tropospheric warming is anomalously strong at least relative to the last ~16,000 years. 
    more » « less
  5. Abstract. Determining the timing and extent of Quaternary glaciations around the globe is critical to understanding the drivers behind climate change and glacier fluctuations. Evidence from the southern mid-latitudes indicates that local glacial maxima preceded the global Last Glacial Maximum (LGM), implying that feedbacks in the climate system or ice dynamics played a role beyond the underlying orbital forcings. To shed light on these processes, we investigated the glacial landforms shaped and deposited by the Lago Argentino glacier (50° S), an outlet lobe of the former Patagonian Ice Sheet, in southern Argentina. We mapped geomorphological features on the landscape and dated moraine boulders and outwash sediments using 10Be cosmogenic nuclides and feldspar infrared stimulated luminescence (IRSL) to constrain the chronology of glacial advance and retreat. We report that the Lago Argentino glacier lobe reached more extensive limits prior to the global LGM, advancing during the middle to late Pleistocene between 243–132 ka and during Marine Isotope Stage 3 (MIS 3), culminating at 44.5 ± 8.0 and at 36.6 ± 1.0 ka. Our results indicate that the most extensive advance of the last glacial cycle occurred during MIS 3, and we hypothesize that this was a result of longer and colder winters, as well as increased precipitation delivered by a latitudinal migration of the Southern Westerly Winds belt, highlighting the role of local and regional climate feedbacks in modulating ice mass changes in the southern mid-latitudes. 
    more » « less