Search

Abstract The accelerated melting of ice on the Antarctic Peninsula and islands in the sub-Antarctic suggests that the cryosphere is edging towards an irreversible tipping point. How unusual is this trend of ice loss within the frame of natural variability, and to what extent can it be explained by underlying climate dynamics? Here, we present new high-resolution reconstructions of long-term changes in the extents of three glaciers on the island of South Georgia (54°S, 36°W), combining detailed analyses of glacial-derived sediments deposited in distal glacier-fed lakes and cosmogenic exposure dating of moraines. We document that the glaciers of South Georgia have gradually retracted since the Antarctic cold reversal (ACR, 14.5–12.8 ka), culminating in the disappearance of at least one of the reconstructed glaciers. The glacier retreat pattern observed in South Georgia suggests a persistent link to summer insolation at 55°S, which intensified during the period from the ACR to approximately 2 ka. It also reveals multi-decadal to centennial climate shifts superimposed on this long-term trend that have resulted in at least nine glacier readvances during the last 10.5 ka. Accompanying meridional changes in the Southern Hemisphere westerlies and their interconnection with local topography may explain these glacier readvances.

Abstract Climate is currently warming due to anthropogenic impact on the Earth’s atmosphere. To better understand the processes and feedbacks within the climate system that underlie this accelerating warming trend, it is useful to examine past periods of abrupt climate change that were driven by natural forcings. Glaciers provide an excellent natural laboratory for reconstructing the climate of the past as they respond sensitively to climate oscillations. Therefore, we study glacier systems and their behavior during the transition from colder to warmer climate phases, focusing on the period between 15 and 10 ka. Using a combination of geomorphological mapping and beryllium-10 surface exposure dating, we reconstruct ice extents in two glaciated valleys of the Silvretta Massif in the Austrian Alps. The mountain glacier record shows that general deglaciation after the Last Glacial Maximum (LGM) was repeatedly interrupted by glacier stabilization or readvance, perhaps during the Oldest Dryas to Bølling transition (landform age: 14.4 ± 1.0 ka) and certainly during the Younger Dryas (YD; 12.9–11.7 ka) and the Early Holocene (EH; 12–10 ka). The oldest landform age indicates a lateral ice margin that postdates the ‘Gschnitz’ stadial (ca. 17–16 ka) and predates the YD. It shows that local inner-alpine glaciers were more extensive until the onset of the Bølling warm phase (ca. 14.6 ka), or possibly even into the Bølling than during the subsequent YD. The second age group, ca. 80 m below the (pre-)Bølling ice margin, indicates glacier extents during the YD cold phase and captures the spatial and temporal fine structure of glacier retreat during this period. The ice surface lowered approximately 50–60 m through the YD, which is indicative of milder climate conditions at the end of the YD compared to its beginning. Finally, the third age group falls into a period of more substantial warming, the YD–EH transition, and shows discontinuous glacier retreat during the glacial to interglacial transition. The new geochronologies synthesized with pre-existing moraine records from the Silvretta Massif evidence multiple cold phases that punctuated the general post-LGM warming trend and illustrate the sensitive response of Silvretta glaciers to abrupt climate oscillations in the past.

Abstract The Last Glacial Termination (T1) featured major changes in global circulation systems that led to a shift from glacial to interglacial climate. While polar ice cores attest to an antiphased thermal pattern at millennial timescales, recent well-dated moraine records from both hemispheres suggest in-phase fluctuations in glaciers through T1, which is inconsistent with the bipolar see-saw paradigm. Here, we present a glacier chronology based on 30 new 10Be surface exposure ages from well-preserved moraines in the Lago Palena/General Vintter basin in northern Patagonia (~ 44°S). We find that the main glacier lobe underwent profound retreat after 19.7 ± 0.7 ka. This recessional trend led to the individualization of the Cerro Riñón glacier by ~ 16.3 ka, which underwent minor readvances at 15.9 ± 0.5 ka during Heinrich Stadial 1, during the Antarctic Cold Reversal with successive maxima at 13.5 ± 0.4, 13.1 ± 0.4, and 13.1 ± 0.5 ka, and a minor culmination at 12.5 ± 0.4 ka during Younger Dryas time. We conclude that fluctuations of Patagonian glaciers during T1 were controlled primarily by climate anomalies brought by shifts in the Southern Westerly Winds (SWW) locus. We posit that the global covariation of mountain glaciers during T1 was linked to variations in atmospheric CO2 (atmCO2) promoted by the interplay of the SWW-Southern Ocean system at millennial timescales.

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 10Be 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.

Alpine glaciers, with their valuable combination of highly sensitive response to climate and near-global extent, are powerful tools for investigating previous and present climate changes. They also represent critical water resources for areas around the globe, with the potential for far-reaching effects in a warming world. Advancements to understand and model glacial changes and the variables influencing them are therefore paramount. Many glacier models fall into one of two endmembers: either highly complex transient models requiring careful tuning of multiple parameters to individual glaciers, or basic empirical correlations of glacier area and length with few considerations for local and regional variations in characteristics. Here we detail a physical steady-state model for alpine glaciers relating directly to glacier mass balance (via the equilibrium line altitude) while retaining the simplicity of other morphology methods, and simultaneously including error estimates. We provide custom MATLAB functions as a user-friendly and generally-applicable method to estimate glacier equilibrium line altitudes from only a limited number of glacier bed topography and glacier width measurements. As a test of the model’s efficacy, we compare the model results for present-day glaciers in the Swiss Alps with previously published estimates of equilibrium line altitudes and intermediate model outputs. • The method estimates glacier equilibrium line altitudes from a limited set of bed topography measurements and constraints on glacier width. • The method is based on continuity equations, reducing the need for empirical coefficients tuned with measured data. • The method uses Monte Carlo sampling and bootstrapping to generate uncertainty bounds on the equilibrium line altitude estimates.