Knowledge of how glaciations formed landscapes is particularly important as receding glaciers currently uncover subglacial landscapes that are prone to a series of natural hazards, but that also bring opportunities for hydropower or water resources. We present high-resolution (1 × 1 m) swath bathymetric data of a proglacial lake in front of the Rhonegletscher (Swiss Alps) that started to form in the early 2000s allowing a look into a freshly uncovered glacier bed and its characterized morphology in an overdeepened setting. The comparison of two surveys from 2015 and 2021 allows an unprecedented quantification of the accumulation and erosion processes in the central lake basin. This highly dynamic environment is characterized by iceberg calving, fluctuating outflow conduits, rapid sedimentation due to particle-laden meltwaters and dumped glacial debris. Assuming constant sediment yield, the Rhone Lake would persist for ~300 years. However, as intense glacier retreat continues in the coming decades, a chain of overdeepened lakes will be revealed that will act as long persisting sediment traps.

Reconstructing Pleistocene glaciation timing and extent is vital for understanding paleoclimate. Whereas late Pleistocene glaciation has been studied extensively in western North American mountain ranges, the glacial history of the western Elk Range in Colorado remains understudied, particularly in the East River watershed, a site of intense scientific focus. Here we use cosmogenic nuclide exposure and depth–profile dating methods to determine the timing of glaciation in the East River watershed. We use glacier modeling to reconstruct paleoglacier extents and quantify past climate conditions. Our findings indicate that the East River glacier retreated from its maximum position approximately 17–18 ka, moving to recessional positions between 13 and 15 ka, before experiencing more substantial retreat to high-elevation cirques around 13 ka. Glacier modeling suggests that the maximum ice extents at 17–18 ka could have been sustained by temperature depressions of approximately −6.5°C compared to modern conditions, assuming consistent precipitation. Additionally, the ice position at 13–15 ka could have been supported by temperature depressions of around −4.0°C. These results offer insights into the deglaciation timeline in the East River watershed and broader western Elk Range as well as paleoclimate conditions during the late Pleistocene, which may aid future research on critical zone evolution in the East River watershed.

We review successes and challenges from five recent subglacial bedrock drilling campaigns intended to find evidence for Antarctic Ice Sheet retreat during warm periods in the geologic past. Insights into times when the polar ice sheets were smaller than present serve as guiding information for modeling efforts that aim to predict the rate and magnitude of future sea level rise that would accompany major retreat of the Antarctic Ice Sheet. One method to provide direct evidence for the timing of deglaciations and minimum extent of prior ice sheets is to extract subglacial bedrock cores for cosmogenic nuclide analysis from beneath the modern ice sheet surface. Here we summarize the lessons learned from five field seasons tasked with obtaining bedrock cores from shallow depths (<120 m beneath ice surface) across West Antarctica since 2016. We focus our findings on drilling efforts and technology and geophysical surveys with ground-penetrating radar. Shallow subglacial drilling provides a high risk, high reward means to test for past instabilities of the Antarctic Ice Sheet, and we highlight key challenges and solutions to increase the likelihood of success for future subglacial drilling efforts in polar regions.

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change.

Cosmogenic nuclide measurements in glacial deposits extend our knowledge of glacier chronologies beyond the observational record. The short half-life of in situ cosmogenic 14C makes it particularly useful for studying glacier chronologies, as resulting exposure ages are less sensitive to nuclide inheritance when compared with more commonly measured, long-lived nuclides. An increasing number of laboratories using an automated process to extract carbon from quartz has led to in situ 14C measurements in Antarctic samples at an accelerating rate over the past decade, shedding light on deglaciation in Antarctica. In situ 14C has had the greatest impact in the Weddell Sea Embayment, where inferences on the thickness of ice and timing of deglaciation were limited by inheritance in other cosmogenic nuclide systems. Future subglacial measurements of the nuclide hold much potential as they can provide direct evidence of proposed Holocene thinning and subsequent re-thickening of parts of the Antarctic ice sheets.

The Autonomous Pinger Unit (APU) is an electro-thermal drill with acoustic instrumentation developed for the project EnEx-RANGE in view of a future space mission for the sub-surface exploration of Saturn's moon Enceladus. A main goal is the development of navigation technology for an acoustic guidance system allowing maneuvering a probe through glacial ice. In total 13 APUs were built and tested in terrestrial analog scenarios on alpine glaciers. The APUs form a spatially distributed network that defines a system of reference for the navigation of the maneuverable probe to a point of interest. The APUs have a novel melting head, slow control systems, and a modern system-on-chip (SoC) module that controls the probe and processes the recorded data. The APUs use acoustic emitters and receivers to measure the transit time of acoustic signals between them, allowing for the position reconstruction of all APUs by trilateration. Several auxiliary sensors monitor the internal state of the probe and assist the position estimation. With this instrumentation, the APUs have the ability of dynamically optimizing themselves within the network by changing their position. This paper gives an overview of the developed APU hardware and presents performance results from the field tests.

Ramparted depressions (doughnut-shaped debris-cored ridges with peat- and/or sediment-filled central basins) are commonly perceived to represent the relict collapsed forms of permafrost ground-ice mounds (i.e. pingos or lithalsas). In Wales, UK, ramparted depressions of Late Pleistocene age have been widely attributed to permafrost-related processes. However, a variety of alternative glacial origins for these enigmatic landforms are also consistent with the available geological and geomorphological evidence, although previous studies have barely considered such alternative processes of formation. From detailed geophysical, sedimentological and remote-sensing studies at two field sites, we demonstrate that: (i) the wastage of stagnating glacier ice is a viable alternative explanation for the formation of ramparted depressions in Wales; (ii) the glacial geomorphology and geology of these landforms is analogous to supraglacial and subglacial landforms from the last Laurentide and Fennoscandian ice sheets; (iii) these landforms have significant potential for characterising the nature of deglaciation around the margins of the Irish Sea during the last glacial cycle, and may record evidence for the overextension and stagnation of the south-eastern margin of the Irish Sea Ice Stream; and (iv) investigations of ramparted depressions within formerly glaciated terrains must consider both glacial and periglacial mechanisms of formation.

Glacial erosion produces vast quantities of fine-grained sediment that has a far-reaching impact on Earth surface processes. To gain a better understanding of the production of glacial silt and clay, we use automated mineralogy to quantify the microstructure and mineralogy of rock and sediment samples from 20 basins in the St. Elias Mountains, Yukon, Canada. Sediments were collected from proglacial streams, while rock samples were collected from ice marginal outcrops and fragmented using electrical pulse disaggregation. For both rock fragments and sediments, we observe a log-normal distribution of grain sizes and a sub-micrometer terminal grain size. We find that the abrasion of silt and clay results in both rounding and the exploitation of through-going fractures. The abundance of inter- versus intragranular fractures depends on mineralogy and size. Unlike the relatively larger grains, where crushing and abrasion are thought to exploit and produce discrete populations of grain sizes, the comminution of fines leads to a grain size, composition and rounding that is continuously distributed across size, and highly dependent on source-rock properties.

Glacier-erosion rates range across orders of magnitude, and much of this variation cannot be attributed to basal sliding rates. Subglacial till acts as lubricating ‘fault gouge’ or ‘sawdust’, and must be removed for rapid subglacial bedrock erosion. Such erosion occurs especially where and when moulin-fed streams access the bed and are unconstrained by supercooling or other processes. Streams also may directly erode bedrock, likely with strong time-evolution. Erosion is primarily by quarrying, aided by strong fluctuations in the water system driven by variable surface melt and by subglacial earthquakes. Debris-bed friction significantly affects abrasion, quarrying and general glacier flow. Frost heave drives cirque headwall erosion as winter cold air enters bergschrunds, creating temperature gradients to drive water flow along premelted films to growing ice lenses that fracture rock, and the glacier removes the resulting blocks. Recent subglacial bedrock erosion and sediment flux are in many cases much higher than long-term averages. Over glacial cycles, evolution of glacial-valley form feeds back strongly on erosion and deposition. Most of this is poorly quantified, with parts open to argument. Glacial erosion and interactions are important to tectonic and volcanic processes as well as climate and biogeochemical fluxes, motivating vigorous research.

At its late Pleistocene maximum, the Laurentide Ice Sheet was the largest ice mass on Earth and a key player in the modulation of global climate and sea level. At the same time, this temperate ice sheet was itself sensitive to climate, and high-magnitude fluctuations in ice extent, reconstructed from relict glacial deposits, reflect past changes in atmospheric temperature. Here, we present a cosmogenic 10Be surface-exposure chronology for the Berlin moraines in the White Mountains of northern New Hampshire, USA, which supports the model that deglaciation of New England was interrupted by a pronounced advance of ice during the Bølling-Allerød. Together with recalculated 10Be ages from the southern New England coast, the expanded White Mountains moraine chronology also brackets the timing of ice sheet retreat in this sector of the Laurentide. In conjunction with existing chronological data, the moraine ages presented here suggest that deglaciation was widespread during Heinrich Stadial 1 event (~18–14.7 ka) despite apparently cold marine conditions in the adjacent North Atlantic. As part of the White Mountains moraine system, the Berlin chronology also places a new terrestrial constraint on the former glacial configuration during the marine incursion of the St. Lawrence River valley north of the White Mountains.

High-resolution digital elevation models of Finland and Sweden based on LiDAR (Light Detection and Ranging) reveal subglacial landforms in great detail. We describe the ice-sheet scale distribution and morphometric characteristics of a glacial landform that is distinctive in morphology and occurs commonly in the central parts of the former Scandinavian Ice Sheet, especially up-ice of the Younger Dryas end moraine zone. We refer to these triangular or V-shaped landforms as murtoos (singular, ‘murtoo’). Murtoos are typically 30–200 m in length and 30–200 m in width with a relief of commonly <5 m. Murtoos have straight and steep edges, a triangular tip oriented parallel to ice-flow direction, and an asymmetric longitudinal profile with a shorter, but steeper down-ice slope. The spatial distribution of murtoos and their geomorphic relation to other landforms indicate that they formed subglacially during times of climate warming and rapid retreat of the Scandinavian Ice Sheet when large amounts of meltwater were delivered to the bed. Murtoos are formed under warm-based ice and may be associated with a non-channelized subglacial hydraulic system that evacuated large discharges of subglacial water.

Stratigraphic evidence and extensive optically stimulated luminescence (OSL) geochronology from an 18-km-long reach of the middle Rangitata Valley, South Island, New Zealand, provide evidence for at least six distinct glacial advances during the last glacial cycle. These include four well-constrained Marine Oxygen Isotope Stage (MIS) 3 and 2 advances at ca. 38 ka, ca. 27 ka, ca. 21 ka and at 18 ka, as well as less well-constrained advances in MIS 4 and/or early MIS 3. Ice occupied a farther downvalley reach of the Rangitata from 38 ka to after 18 ka, indicating that near-full glacial conditions persisted for most of the last 20 ka of the last glaciation, though the glacier still fluctuated significantly, as reflected by the numerous distinguishable advances. Global or regional cooling alone cannot explain the persistence of near-maximum glacial conditions for this extended period, nor can it explain the occurrence of the largest advances ca. 32 ka. Instead, we invoke the northward expansion of the westerlies during MIS 3 as the cause for the early widespread glaciation, wherein enhanced westerly flow under moderate cooling maximised glacial extents. Local insolation favoured extended MIS 3 glaciation until ca. 32 ka. Increasing summer insolation gradually reduced glacial extents after ca. 28 ka.

Our study combines new geological and paleoecological information to reconstruct the glacial history and terrestrial paleoenvironments on Haida Gwaii during the advance phase of the Fraser glaciation (Marine Isotope Stage 2). At Cape Ball on eastern Graham Island, five accelerator mass spectrometry radiocarbon ages ranging from 23,200±280 to 26,650±390 14C yr BP (ca. 27,000–31,400 cal yr BP) record the earliest approach of mainland glaciers to Haida Gwaii. Abundant marine dinoflagellate cysts indicate isostatic depression by glacial ice in Hecate Strait to the east. At Mary Point on the north coast of Graham Island, similar outwash of a piedmont lobe advancing westward along Dixon Entrance preserves plant remains dated from 19,270±360 to 23,740±300 14C yr BP (22,500–28,600 cal yr BP). These sediments also contain marine indicators. Plant macrofossils, pollen, and invertebrates support the geological evidence of a proglacial environment under a colder-than-present macroclimate. Although some trees were likely present on Graham Island at this time, tundra-like plant communities dominated low-lying areas. A large area that appears to have been ice-free during this time is a portion of the continental shelf off the east coast of Moresby Island, referred to provisionally as the “Hecate Refugium.”

The Glasgow area has a combination of highly variable superficial deposits and a legacy of heavy industry, quarrying and mining. These factors create complex foundation and hydrological conditions, influencing the movement of contaminants through the subsurface and giving rise locally to unstable ground conditions. Digital geological three-dimensional models developed by the British Geological Survey are helping to resolve the complex geology underlying Glasgow, providing a key tool for planning and environmental management. The models, covering an area of 3200km2 to a depth of 1.2km, include glacial and post-glacial deposits and the underlying, faulted Carboniferous igneous and sedimentary rocks. Control data, including 95,000 boreholes, digital mine plans and published geological maps, were used in model development. Digital outputs from the models include maps of depth to key horizons, such as rockhead or depth to mine workings. The models have formed the basis for the development of site-scale high-resolution geological models and provide input data for a wide range of other applications from groundwater modelling to stochastic lithological modelling.

Isla de los Estados is a mountainous island southeast of Tierra del Fuego, in southernmost South America. Its central and eastern parts have an alpine topography, transected by U-shaped valleys, small, partly over-deepened fjords, and a multitude of abandoned cirques, all associated with extensive former local glaciations. Traces of glacial erosion generally reach 400–450 m a.s.l., and above that trimline a distinct sharp-edged nunatak derived landscape is present. The westernmost part of the island has a lower, more subdued topography, reflecting its "softer" geology but possibly also over-running and erosion by mainland-derived ice streams. The present study concentrated on glacigenic sediment sequences exposed along coastal erosional cliffs. A combination of OSL and 14C datings show that these sediments mostly date from the latest (Wisconsinan/Weichselian) glacial cycle, i.e. from the last ca. 100 ka with the oldest (glaciolacustrine) deposits possibly as old as 90–80 ka. The upper parts of overlying tills, with associated lateral and terminal moraines from glaciers that expanded onto an eustatically exposed dry shelf north of the island, date from the last global glacial maximum (LGM). Radiocarbon ages of peat and lake sediments indicate that deglaciation began 17–16 cal ka BP.

A large ice sheet still covered almost all of Maine and eastern New England until ca. 15 cal ka BP, reaching south of 45 °S, despite rising summer insolation intensity and major ice recession elsewhere outside the North Atlantic region. Furthermore, the well-studied moraine belt along eastern coastal Maine, including the prominent Pineo Ridge delta/moraine complex and Pond Ridge moraine, indicates repeated readvances and stillstands between ca. 16 and 15 cal ka BP. This moraine belt reflects a considerable ice sheet response over eastern North America during this time period, coeval with the latter half of the European Oldest Dryas period. Moraine deposition was concurrent with reduction or elimination of North Atlantic meridional overturning, starting with the earlier onset of peak IRD and Heinrich Event 1 (HE-1). The existing 14C chronology suggests that the coastal moraine belt and the persistence of the ice sheet until ∼ 15 cal ka BP was a response to the severe cooling of the North Atlantic region after ∼ 17 cal ka BP.

At least three sets of moraines mark distinct glacial stands since the last glacial maximum (LGM) in the Three Sisters region of the Oregon Cascade Range. The oldest stand predates 8.1 ka (defined here as post-LGM), followed by a second between ∼ 2 and 8 ka (Neoglacial) and a third from the Little Ice Age (LIA) advance of the last 300 years. The post-LGM equilibrium line altitudes were 260 ± 100 m lower than that of modern glaciers, requiring 23 ± 9% increased winter snowfall and 1.4 ± 0.5°C cooler summer temperatures than at present. The LIA advance had equilibrium line altitudes 110 ± 40 m lower than at present, implying 10 ± 4% greater winter snowfall and 0.6 ± 0.2°C cooler summer temperatures.

Late Pleistocene glaciation of the Australian mainland was restricted to a small area of the southeastern highlands. Geomorphic mapping of the area and exposure dating using the in situ produced cosmogenic isotope 10Be provides evidence for at least two distinct glaciations. The Early Kosciuszko glaciation consisted of a single glacier advance before 59,300 ± 5400 years ago (Snowy River Advance). The Late Kosciuszko glaciation comprised three glacier advances 32,000 ± 2500 (Headley Tarn Advance), 19,100 ± 1600 (Blue Lake Advance), and 16,800 ± 1400 years ago (Mt. Twynam Advance). The Early Kosciuszko glaciation was the most extensive and the Late Kosciuszko advances were progressively less extensive. These periods of glaciation in the highlands correspond to episodes of periglacial activity and peaks in lake levels and river discharge at lower elevations in southeastern Australia. Glacier advances on the Kosciuszko Massif correlate with advances in Tasmania, South America, and New Zealand and are broadly representative of hemispheric climate changes during the last glacial cycle.

Stagnation of the terminal region of a glacier occurs in response to sufficiently large and rapid climatic warming, so the presence of stagnation deposits provides quantitative information on the climate change that forced retreat. Here we use a simple flow-line glacier model to investigate the relationship between stagnation, climate forcing and aspects of the glacier bed. For climatic warming greater than the threshold to cause stagnation, larger or faster warming events cause longer regions of a glacier to stagnate. Smaller or slower warming episodes, below the threshold for stagnation, cause retreat while active flow persists along the entire glacier length. The threshold for stagnation depends not only on the climatic forcing but also on many other aspects of the glacier, with stagnation favored by factors including a lower mean bed slope with greater roughness. Quantitative determination of the climatic forcing consistent with the occurrence or absence of stagnation deposits requires that these site-specific characteristics be incorporated in modeling.

Of numerous theories of drumlin genesis, none has been widely accepted. It seems evident, however, that some form of positive feedback process is involved. Under certain circumstances perturbations are amplified. Herein we suggest that patchy areas of frozen bed provide the initial perturbation. Such frozen patches may occur in local areas underlain by material of lower thermal conductivity or on slight topographic highs. Drag exerted by the frozen patch deflects ice flow into its lee, dragging with it mobile till eroded from the thawed area. The energy balance is such that this till likely refreezes, either producing a topographic perturbation or amplifying an existing one. The resulting topography then deflects more of the geothermal heat away from the developing hill and into the adjacent trough, resulting in a positive feedback. Once the thermal perturbation exceeds a critical (though as yet undefined) level, melting may decouple the ice from the bed, preventing further entrainment of till from thawed areas, and thus limiting the height and length of the drumlin.