ABSTRACT
After the Last Glacial Maximum (LGM) global mean sea level (GMSL) rise was characterized by rapid increases over short (decadal to centennial) timescales superimposed on a longer term ...secular rise and these have been termed meltwater pulses (MWPs). In this paper we review the timing, impact and nature of these and the effects of rapid drainage of large post‐glacial MWPs into the world's oceans. We identify nine MWPs, four of which occurred in three periods of rapid sea level rise (19.5–18.8, 14.8–13.0 and 11.5–11.1 ka BP). The rest are dated to the period during the Early Holocene sea level rise after 11 ka BP. We show that drainage of the known post‐glacial lakes in total produced less than around 1.2 m of the 125 m of GMSLR since the LGM.
Many Himalayan glaciers are characterised in their lower reaches by a rock debris layer. This debris insulates the glacier surface from atmospheric warming and complicates the response to climate ...change compared to glaciers with clean-ice surfaces. Debris-covered glaciers can persist well below the altitude that would be sustainable for clean-ice glaciers, resulting in much longer timescales of mass loss and meltwater production. The properties and evolution of supraglacial debris present a considerable challenge to understanding future glacier change. Existing approaches to predicting variations in glacier volume and meltwater production rely on numerical models that represent the processes governing glaciers with clean-ice surfaces, and yield conflicting results. We developed a numerical model that couples the flow of ice and debris and includes important feedbacks between debris accumulation and glacier mass balance. To investigate the impact of debris transport on the response of a glacier to recent and future climate change, we applied this model to a large debris-covered Himalayan glacier—Khumbu Glacier in Nepal. Our results demonstrate that supraglacial debris prolongs the response of the glacier to warming and causes lowering of the glacier surface in situ, concealing the magnitude of mass loss when compared with estimates based on glacierised area. Since the Little Ice Age, Khumbu Glacier has lost 34% of its volume while its area has reduced by only 6%. We predict a decrease in glacier volume of 8–10% by AD2100, accompanied by dynamic and physical detachment of the debris-covered tongue from the active glacier within the next 150 yr. This detachment will accelerate rates of glacier decay, and similar changes are likely for other debris-covered glaciers in the Himalaya.
•Debris-covered glaciers show a differing response to climate change to clean-ice glaciers.•Many glaciers in the Himalaya are debris-covered, particularly in the Everest region.•We present the first dynamic model of debris transport and feedbacks with mass balance.•Debris-covered Khumbu Glacier in Nepal will lose 8–10% volume by AD2100.•The glacier tongue will detach from the active glacier and accelerate glacier decay.
We present PATICE, a GIS database of Patagonian glacial geomorphology and recalibrated chronological data. PATICE includes 58,823 landforms and 1,669 geochronological ages, and extends from 38°S to ...55°S in southern South America. We use these data to generate new empirical reconstructions of the Patagonian Ice Sheet (PIS) and subsequent ice masses and ice-dammed palaeolakes at 35 ka, 30 ka, 25 ka, 20 ka, 15 ka, 13 ka (synchronous with the Antarctic Cold Reversal), 10 ka, 5 ka, 0.2 ka and 2011 AD. At 35 ka, the PIS covered of 492.6 x103 km2, had a sea level equivalent of ~1,496 mm, was 350 km wide and 2090 km long, and was grounded on the Pacific continental shelf edge. Outlet glacier lobes remained topographically confined and the largest generated the suites of subglacial streamlined bedforms characteristic of ice streams. The PIS reached its maximum extent by 33 – 28 ka from 38°S to 48°S, and earlier, around 47 ka from 48°S southwards. Net retreat from maximum positions began by 25 ka, with ice-marginal stabilisation then at 21 – 18 ka, which was then followed by rapid, irreversible deglaciation. By 15 ka, the PIS had separated into disparate ice masses, draining into large ice-dammed lakes along the eastern margin, which strongly influenced rates of recession. Glacial readvances or stabilisations occurred at least at 14 – 13 ka, 11 ka, 6 – 5 ka, 2 – 1 ka, and 0.5 – 0.2 ka. We suggest that 20th century glacial recession (% a-1) is occurring faster than at any time documented during the Holocene.
Display omitted
Glaciers respond sensitively to climate change and paleoclimate conditions can be inferred from the geomorphic evidence by reconstructing equilibrium line altitude (ELA) and by modeling glacier mass ...balance and ice flow. However, such work has not been extensively carried out on the vast Tibetan Plateau, and thus the knowledge of glacier sensitivity to paleoclimate change and paleoglacier volume is still lacking on the plateau. Recent improvements in understanding glacial extent and chronology in the Payuwang Valley, southern Tibetan Plateau, present an opportunity to estimate the paleoclimate conditions during the Last Glacial cycle. Using a coupled mass-balance and ice-flow model, this study tests the glacier sensitivity to ELA change in the Payuwang Valley. By fitting modeled glaciers to the corresponding moraine positions, we reconstruct glacier extents probably occurred at oxygen isotope stage 3 (MIS 3) and Last Glacial Maximum (LGM) in the valley and infer the related climate conditions that could have supported these glacier advances. The glacier sensitivity tests indicate that in the Payuwang Valley, the glaciers responded to centennial-scale shifts in climate. With ELAs below 5600m asl, the tributary glaciers coalesce into the main valley and expand rapidly there. The model results also show that the Payuwang Valley contained ice volumes of ~1.58×109m3 and ~1.23×109m3 with ELA lowering values of 370m and 340m in the MIS 3 and LGM glacial advances, respectively. By examining other independent paleoclimatic reconstructions, we conclude that during the MIS 3 glacier advance, the temperature was only 1.5°C lower or even 0.2°C higher than the period of 2005–2008 with precipitation 40–100% higher than the 2005–2008 amount; and during the LGM glacier advance, the temperature was 3.3–4.4°C lower than the period of 2005–2008 with the precipitation 30–70% lower than the 2005–2008 amount. These estimates add to the knowledge of better understanding the relationship between glacier dynamic and climatic change on the Tibetan Plateau.
•A coupled mass-balance and ice-flow model was used.•Glacier sensitivity to ELA change was tested in a Tibet valley.•Glacier extents were reconstructed for the Last Glacial cycle.•Paleoclimates that support the glacier extents were inferred.
As the world’s glaciers recede in response to a warming atmosphere, a change in the magnitude and frequency of related hazards is expected. Among the most destructive hazards are glacier lake ...outburst floods (GLOFs), and their future evolution is concerning for local populations and sustainable development policy. Central to this is a better understanding of triggers. There is a long‐standing assumption that earthquakes are a major GLOF trigger, and seismic activity is consistently included as a key hazard assessment criterion. Here, we provide the first empirical evidence that this assumption is largely incorrect. Focusing on the Tropical Andes, we show that, of 59 earthquakes (1900–2021) the effects of which intersect with known glacier lakes, only one has triggered GLOFs. We argue that, to help develop climate resilient protocols, the focus for future assessments should be on understanding other key GLOF drivers, such as thawing permafrost and underlying structural geology.
Plain Language Summary
Climate change is increasing glacier melt in high mountains and increasing the size and number of glacial lakes. Over time, these lakes may drain catastrophically to generate glacier lake outburst floods (GLOFs), which can devastate downstream communities and destroy valuable infrastructure such as roads, bridges, and hydroelectric power facilities. As a result, many risk assessments have been carried out to understand what the triggers of GLOFs might be, in order to better predict their occurrence. One of the main triggers has been assumed to be earthquakes, but this association has not been properly tested. In this paper, we use earthquake and GLOF data from the Peruvian and Bolivian Andes to test this and find that there is little association between earthquakes and GLOFs. We conclude by arguing that focus needs to be on other GLOF triggers, and earthquake activity should be used as a secondary – not primary – indicator in GLOF susceptibility studies.
Key Points
Glacier lake outburst floods (GLOFs) are a major risk concern, with earthquakes being commonly considered effective triggers in hazard assessments
Of 59 earthquakes (1900–2021), the effects of which intersect lakes in the Tropical Andes, only one triggered GLOFs
We suggest that earthquake activity should be used as a secondary – not primary – susceptibility indicator for GLOFs
Subglacial erosion beneath glaciers occurs predominantly by abrasion and plucking, producing distinct erosional forms. The controls on the relative importance of abrasion vs. plucking are poorly ...understood. On the one hand, glacial conditions that favour or suppress cavity formation (ice velocity, ice thickness, and water pressure) are thought to favour plucking or abrasion, respectively. Conversely, bedrock properties are also known to control landforms, but this has rarely been analysed quantitatively. In this study we compare landforms and bedrock properties of sandstone and quartzite at the bed of a palaeo-ice stream near Ullapool in NW Scotland. The boundary between the rock types is at right angles to the westward palaeo-ice flow, and palaeoglacial conditions on both rock types were similar. We report quantitative parameters for bedrock properties (Schmidt hammer hardness and joint spacing) and use morphometric parameters to analyse the landforms. Torridon sandstone is soft but thick-bedded and with a wide joint spacing. Erosional bedforms include roche moutonnées with smoothed tops and concave stoss sides, whalebacks, and elongate p-forms, indicating a high proportion of abrasion over plucking. Cambrian quartzite is hard but thin-bedded with narrow joint spacing. Erosional landforms are angular to subangular with abundant plucked lee faces, suggesting a high proportion of plucking over abrasion. Hardness and joint spacing thus exert a strong control on subglacial erosional landforms and the mechanisms that formed them. Thus glacial conditions (ice velocity, ice thickness) can only be inferred from glacial erosional landforms if the effects of bedrock properties of the substrate are considered.
► We compare subglacial erosional landforms on quartzite and sandstone. ► Quartzite landforms are formed by plucking; sandstone landforms by abrasion. ► We also measured hardness and joint spacing of the two rock types. ► Abrasion is favoured by soft rocks and wide joint spacing. ► Plucking is favoured by hard rocks and close joint spacing.
The prevalence and increased frequency of high-magnitude Glacial Lake Outburst Floods (GLOFs) in the Chilean and Argentinean Andes suggests this region will be prone to similar events in the future ...as glaciers continue to retreat and thin under a warming climate. Despite this situation, monitoring of glacial lake development in this region has been limited, with past investigations only covering relatively small regions of Patagonia. This study presents new glacial lake inventories for 1986, 2000 and 2016, covering the Central Andes, Northern Patagonia and Southern Patagonia. Our aim was to characterise the physical attributes, spatial distribution and temporal development of glacial lakes in these three sub-regions using Landsat satellite imagery and image datasets available in Google Earth and Bing Maps. Glacial lake water volume was also estimated using an empirical area-volume scaling approach. Results reveal that glacial lakes across the study area have increased in number (43%) and areal extent (7%) between 1986 and 2016. Such changes equate to a glacial lake water volume increase of 65 km3 during the 30-year observation period. However, glacial lake growth and emergence was shown to vary sub-regionally according to localised topography, meteorology, climate change, rate of glacier change and the availability of low gradient ice areas. These and other factors are likely to influence the occurrence of GLOFs in the future. This analysis represents the first large-scale census of glacial lakes in Chile and Argentina and will allow for a better understanding of lake development in this region, as well as, providing a basis for future GLOF risk assessments.
•New glacial lake inventories covering the Central and Patagonian Andes presented for 1986, 2000 and 2016.•Glacial lake water volume was estimated using an empirical area-volume scaling approach.•Glacial lakes across the study area have increased in number (43%) and areal extent (7%) between 1986 and 2016.•Glacial lake growth and emergence was shown to vary sub-regionally according to localised factors.•21 previously unreported GLOF events were identified through the analysis of Landsat imagery.
South American glaciers, including those in Patagonia, presently contribute the largest amount of meltwater to sea level rise per unit glacier area in the world. Yet understanding of the mechanisms ...behind the associated glacier mass balance changes remains unquantified partly because models are hindered by a lack of knowledge of subglacial topography. This study applied a perfect-plasticity model along glacier centre-lines to derive a first-order estimate of ice thickness and then interpolated these thickness estimates across glacier areas. This produced the first complete coverage of distributed ice thickness, bed topography and volume for 617 glaciers between 41°S and 55°S and in 24 major glacier regions. Maximum modelled ice thicknesses reach 1631m±179m in the South Patagonian Icefield (SPI), 1315m±145m in the North Patagonian Icefield (NPI) and 936m±103m in Cordillera Darwin. The total modelled volume of ice is 1234.6km3±246.8km3 for the NPI, 4326.6km3±865.2km3 for the SPI and 151.9km3±30.38km3 for Cordillera Darwin. The total volume was modelled to be 5955km3±1191km3, which equates to 5458.3Gt±1091.6Gt ice and to 15.08mm±3.01mm sea level equivalent (SLE). However, a total area of 655km2 contains ice below sea level and there are 282 individual overdeepenings with a mean depth of 38m and a total volume if filled with water to the brim of 102km3. Adjusting the potential SLE for the ice volume below sea level and for the maximum potential storage of meltwater in these overdeepenings produces a maximum potential sea level rise (SLR) of 14.71mm±2.94mm. We provide a calculation of the present ice volume per major river catchment and we discuss likely changes to southern South America glaciers in the future. The ice thickness and subglacial topography modelled by this study will facilitate future studies of ice dynamics and glacier isostatic adjustment, and will be important for projecting water resources and glacier hazards.
Display omitted
•New outlines and new centrelines for 617 glaciers between 41°S and 55°S•Ice thickness statistics and ice volume per glacier reported•Ice below sea level and within overdeepenings quantified•Ice volume per major hydrological catchment determined
Mountain glaciers are key indicators of climate change. Their response is revealed by the environmental equilibrium-line altitude (ELA), i.e. the regional altitude of zero mass balance averaged over ...a long period of time. We introduce a simple approach for distributed modelling of the environmental ELA over the entire European Alps based on the parameterization of ELA in terms of summer temperature and annual precipitation at a glacier. We use 200 years of climate records and forecasts to model environmental ELA from 1901 to 2100 at 5 arcmin grid cell resolution. Historical environmental ELAs are reconstructed based on precipitation from the Long-term Alpine Precipitation reconstruction (LAPrec) dataset and temperature from the Historical Instrumental climatological Surface Time series of the greater Alpine region (HISTALP). The simulations of future environmental ELAs are forced with high-resolution EURO-CORDEX regional climate model projections for the European domain using three different greenhouse gas emissions scenarios (Representative Concentration Pathways, RCP). Our reconstructions yielded an environmental ELA across the European Alps of 2980 m above sea level for the period 1901−1930, with a rise of 114 m in the period 1971−2000. The environmental ELA is projected to exceed the maximum elevation of 69%, 81% and 92% of the glaciers in the European Alps by 2071−2100 under mitigation (RCP2.6), stabilization (RCP4.5) and high greenhouse gas emission (RCP8.5) scenarios, respectively.
This paper examines ice-sheet wide variations in subglacial thermal regime and ice dynamics using the landform record exposed on the beds of former mid-latitude ice sheets (the Laurentide, ...Cordilleran, Fennoscandian and British-Irish Ice Sheets). We compare the landform patterns beneath these former ice sheets to the flow organisation beneath parts of the contemporary Antarctic Ice Sheet inferred from RADARSAT-1 Antarctic Mapping Project (RAMP) data. The evidence preserved in the landform record and observed on contemporary ice masses can be grouped into four major ice-dynamical components that collectively define the subglacial thermal organisation (STO) of ice sheets. These ice-dynamical components are frozen-bed patches, ice streams, ice-stream tributaries and lateral shear zones. Frozen-bed patches appear at a wide range of spatial scales, spanning four orders of magnitude. In some areas, frozen-bed zones comprise large proportions of the bed (e.g. near the ice divide in continental areas), whilst in other areas they constitute isolated “islands” in areas dominated by thawed-bed conditions. Ice streams, narrow zones of fast flow in ice sheets that are otherwise dominated by slow sheet flow, are also common features of Quaternary ice sheets. Tributaries to ice streams flow at velocities intermediate between full ice-stream and sheet flow, and may divert ice drainage from one primary ice-stream corridor to an adjacent one. Sharp lateral boundaries between landforms indicate sliding and non-sliding conditions, respectively. These lateral boundaries represent important discontinuities in the glacial landscape and mark the location of shear zones between thawed-bed ice streams and intervening frozen-bed areas. We use the landform evidence in the area around Great Bear Lake, Canada to trace the evolution of an ice-stream web through time, demonstrating that frozen-bed patches are integral components of this complex system. We conclude that frozen-bed patches are important for the stability of ice sheets because they laterally constrain and isolate peripheral drainage basins and their ice streams.