This paper presents geomorphological and historical records of the surge-type behaviour of Hansbreen, one of the most studied tidewater glaciers in Svalbard. The surge-type behaviour of the glacier has not been considered before due to the lack of evidence of this phenomenon. We integrate geomorphological mapping of the terrestrial and submarine forefields with historical data from the 19th and 20th centuries to reconstruct the glacier dynamics and identify the possible timing of surging. Landform assemblages are representative of the surging glacier landsystem, including crevasse-squeeze ridges (CSRs) and submarine streamlined glacial lineations. Abundant CSRs in the outer part of the terrestrial forefield were also documented in the 1980s, but most have been obliterated since then. We suggest the identified surge landsystem was produced during a surge of Hansbreen detected from photographs taken during the Austro-Hungarian expedition in 1872. Historical photogrammetric photos from the Norwegian expedition in 1918 revealed surge-diagnostic features in the glacier surface, including a folded medial moraine and a dense, complex network of crevasses. A potential next surge remains questionable in the following decades due to the low-lying accumulation area of the main stream hindering the mass build-up, but potential surges of the tributary glaciers should not be excluded.

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.

Due to the limited landmasses in the Southern Hemisphere, we must rely on data from sub-Antarctic islands within the Southern Ocean to record historical climate patterns. Over the past few decades, glaciers throughout the Southern Ocean region have experienced a noticeable retreat, especially in the Kerguelen Archipelago, whose glacial landforms offer valuable insights into long-term climate fluctuations. Our comprehensive glacial geomorphological study conducted in its remote north-western region meticulously examines morainic complexes from smaller cirque glaciers and larger outlet glaciers stemming from the Cook Ice Cap. We mapped these landforms to reconstruct historical glacier extents during the Holocene. The surface area of the three main glaciers had decreased in 1962–1964 by only 35% compared to their maximum extents, whereas surface area changes across 12 time intervals spanning from 1962 to 2019 from aerial and satellite imagery reveal a cumulative reduction of 43.5%. Additionally, we modelled changes in glacier thickness and equilibrium-line altitude for the key glaciers at three distinct stages: 1) their maximum extent before 1962, 2) the early 1960s and 3) 2019. This multifaceted analysis contributes valuable insights into the dynamics of Kerguelen's glaciers and the broader implications for understanding past and ongoing climate dynamics in the Southern Hemisphere.

Spatial variability in bed topography, characterized as bed roughness, impacts ice-sheet flow and organization and can be used to infer subglacial conditions and processes, yet is difficult to quantify due to sparse observations. Paleo-subglacial beds of formerly expanded glaciers found across the Antarctic continental shelf are well preserved, have relatively limited post-glacial sediment cover and contain glacial landforms that can be resolved at sub-meter vertical scales. We analyze high-resolution bathymetry offshore of Pine Island and Thwaites glaciers in the Amundsen Sea to explore spatial variability of bed roughness where streamlined subglacial landforms allow for the determination of ice-flow direction. We quantify bed roughness using std dev. and Fast Fourier Transform methods, each employed at local (100 km) and regional (101–2 km) scales and in along- and across-flow orientations to determine roughness expressions across spatial scales. We find that the magnitude of roughness is impacted by the parameters selected – which are often not sufficiently reported in studies – to quantify roughness. Important spatial patterns can be discerned from high-resolution bathymetry, highlighting both its usefulness in identifying patterns of streaming ice flow and underscores the need for a standardized way of characterizing topographic variability.

Glaciers in the Russian High Arctic have undergone accelerated mass loss due to atmospheric and oceanic warming in the Barents–Kara Sea region. Most studies have concentrated on the western Barents–Kara sector, despite evidence of accelerating mass loss as far east as Severnaya Zemlya. However, long-term trends in glacier change on Severnaya Zemlya are largely unknown and this record may be complicated by surge-type glaciers. Here, we present a long-term assessment of glacier change (1965–2021) on Severnaya Zemlya and a new inventory of surge-type glaciers using declassified spy-satellite photography (KH-7/9 Hexagon) and optical satellite imagery (ASTER, Sentinel-2A, Landsat-4/5 TM and 8 OLI). Glacier area reduced from 17 053 km2 in 1965 to 16 275 in 2021 (−5%; mean: −18%, max: −100%), with areal shrinkage most pronounced at land-terminating glaciers on southern Severnaya Zemlya, where there is a recent (post-2010s) increase in summer atmospheric temperatures. We find that surging may be more widespread than previously thought, with three glaciers classified confirmed as surge-type, eight as likely to have surged and nine as possible, comprising 11% of Severnaya Zemlya's 190 glaciers (37% by area). Under continued warming, we anticipate accelerated retreat and increased likelihood of surging as basal thermal regimes shift.

This study presents the first light detection and ranging (LiDAR)-based morphometric description of a water blister from a past ice-sheet bed caused by rapid supraglacial drainage. The blister formed during the rapid early Holocene deglaciation of the Fennoscandian Ice Sheet (FIS). It is located in southwest Finland within a subglacial meltwater route interpreted to represent the transition from a distributed to a channelized drainage system. A LiDAR digital elevation model was supplemented with sedimentological and ground-penetrating radar data on blister outflow channels and sedimentology of downflow polymorphous mounds and ridges (PMRs). Unlike the water blisters recorded from the rapid drainage of supraglacial lakes on the Greenland Ice Sheet, the smaller blister size here was either due to crevasse or moulin drainage, or was a supraglacial lake drainage that tapped into a pre-existing, relatively efficient drainage system and related semi-sorted sediments, promoting rapid drainage and reworking of PMRs along the meltwater route. The preservation potential or exposure probability of blister marks is presumably low but they can provide important information about evolution of subglacial drainage systems that is of value to modern interpretations of glacial hydrology.

Extending the record of glacier area changes into the past improves our understanding of climate change impacts. Although analogue maps showing historic glacier extents are abundant, digital outlines from before the satellite era are sparse as the digitisation of moraines and trimlines on freely available satellite images is challenging. With the now available very high-resolution images provided by Web Map Services (WMS), new doors are open for the precise digitisation. Here, we used the ESRI WMS to digitise Little Ice Age (LIA) glacier extents and present area changes since the LIA in four selected regions along with a detailed uncertainty analysis. We used modern glacier outlines as a starting point and additionally consulted Sentinel-2 images, the ArcticDEM and historic maps for interpretation. Dating records from the literature allowed calculating area change rates. In total, 493 LIA glaciers (4640 km2, now 891 ice bodies with 3590 km2) were digitised, yielding relative area changes of −20% (−0.14% a−1), −15% (−0.10% a−1), −26% (−0.16% a−1) and −61% (−0.19% a−1) for Alaska, Baffin Island, Novaya Zemlya and the tropics, respectively. The ESRI WMS images are a great asset to precisely map moraines and trimlines, but information about the timing of the related extents requires further sources.

Ice streams are sites of ice-sheet drainage and together with other processes, such as calving, have an impact on deglaciation rates and ice-sheet mass balance. Proglacial lake deposits provide records of ice-sheet deglaciation and have the potential to supplement other paleoclimate records. Oneida Lake, northeastern USA, contains a thick proglacial lake sequence that buries evidence of ice streaming and a paleo-calving margin that developed during retreat of the Laurentide Ice Sheet. Previous high-resolution digital elevation models identified the Oneida Ice Stream from glacial landforms northwest of the lake. In this study, we utilize seismic refractions from a multichannel seismic (MCS) reflection dataset to estimate the thickness of glacial deposits using seismic tomography. With this method we constrain the depth to top of Paleozoic strata, especially in areas where the reflection data yielded poor outcomes and validate our reflection data in regions of good coverage. We demonstrate that where long offset seismic data are available, the first-arrival tomography method is useful in studies of formerly glaciated basins. Our study identifies a ~108 m thick sedimentary section and potentially long paleoclimate record in Oneida Lake, and identifies a paleotopographic low that likely encouraged formation of the Oneida Ice Stream.

We use a mathematical model to investigate the effect of basal topography and ice surface slope on transport and deposition of sediment within a water-filled subglacial channel. In our model, three zones of different behaviour occur. In the zone furthest upstream, variations in basal topography lead to sediment deposition under a wide range of conditions. In this first zone, even very small and gradually varying basal undulations (~5 m amplitude) can lead to the deposition of sediment within a modelled channel. Deposition is concentrated on the downstream gradient of subglacial ridges, and on the upstream gradient of subglacial troughs. The thickness and steepness of the ice sheet has a substantial impact on deposition rates, with shallow ice profiles strongly promoting both the magnitude and extent of sediment deposition. In a second zone, all sediment is transported downstream. Finally, a third zone close to the ice margin is characterised by high rates of sediment deposition. The existence of these zones has implications for esker formation and the dynamics of the subglacial environment.

Sub-Antarctic Marion Island's glacial history has acted as a control on abiotic terrestrial processes and the colonization and distribution of biotic species found on the island today. Recent chronological studies have shown an early deglaciation of the island and identified new geomorphological features associated with past ice dynamics. These permit a reassessment of ice extent during and after the island's last local glacial maximum. In this paper, we provide a revised reconstruction of the island's palaeo-ice extent by using a geomorphology-based approach to delineate palaeo-ice margins and demarcate possible glacial basins. The model presented here provides the needed spatial context for future studies on the variations in the distribution of species (e.g. microorganisms and plant species) and abiotic processes and forms (e.g. soil development and periglacial landforms). In addition, it highlights areas that require improved geophysical assessment in order to produce a more complete island-scale reconstruction of former ice extents (e.g. the west coast).

Palaeo-glacial landforms can give insights into bed roughness that currently cannot be captured underneath contemporary-ice streams. A few studies have measured bed roughness of palaeo-ice streams but the bed roughness of specific landform assemblages has not been assessed. If glacial landform assemblages have a characteristic bed-roughness signature, this could potentially be used to constrain where certain landform assemblages exist underneath contemporary-ice sheets. To test this, bed roughness was calculated along 5 m × 5 m resolution transects (NEXTMap DTM, 5 m resolution), which were placed over glacial landform assemblages (e.g. drumlins) in the UK. We find that a combination of total roughness and anisotropy of roughness can be used to define characteristic roughness signatures of glacial landform assemblages. The results show that different window sizes are required to determine the characteristic roughness for a wide range of landform types and to produce bed-roughness signatures of these. Mega scale glacial lineations on average have the lowest bed-roughness values and are the most anisotropic landform assemblage.

We describe the glacial geomorphology and initial geochronology of two ice-free valley systems within the Neptune Range of the Pensacola Mountains, Antarctica. These valleys are characterized by landforms associated with formerly more expanded ice sheet(s) that were at least 200 m thicker than at present. The most conspicuous features are areas of supraglacial debris, discrete debris accumulations separated from modern-day ice and curvilinear ridges and mounds. The landsystem bears similarities to debris-rich cold-based glacial landsystems described elsewhere in Antarctica and the Arctic where buried ice is prevalent. Geochronological data demonstrate multiple phases of ice expansion. The oldest, occurring > 3 Ma, overtopped much of the landscape. Subsequent, less expansive advances into the valleys occurred > 2 Ma and > ~1 Ma. An expansion of some local glaciers occurred < 250 ka. This sequence of glacial stages is similar to that described from the northernmost massif of the Pensacola Mountains (Dufek Massif), suggesting that it represents a regional signal of ice-sheet evolution over the Plio-Pleistocene. The geomorphological record and its evolution over millions of years makes the Neptune Range valleys an area worthy of future research and we highlight potential avenues for this.

Studying subglacial drainage networks is important for understanding the potential relationship between channel dynamics and rapid glacier recession as well as the role of subglacial channels in subglacial sediment evacuation. In order to delineate the planform geometry of snout marginal subglacial channels, densely spaced ground-penetrating radar (GPR) measurements at a frequency of ~70 MHz were carried out over the snout marginal zones of two temperate glaciers in the southwestern Swiss Alps, the Haut Glacier d'Arolla and the Glacier d'Otemma. Three-dimensional (3-D) data processing and amplitude analysis of the GPR reflection along the glacier bed was used to map the channels. At the Haut Glacier d'Arolla, two relatively straight channels of several meters in width were identified. The positions of these channels correspond well with the locations of channel outlets at the glacier terminus, as well as with fractures appearing on the glacier surface one month after the GPR data acquisition. The latter are believed to represent the beginning of ice collapse above the subglacial channels. At the Glacier d'Otemma, a major subglacial conduit was detected with similar dimensions to those identified at the Haut Glacier d'Arolla, but greater sinuosity. The position of this channel was confirmed by drone-based imagery acquired after glacier margin collapse. Our results confirm that high-density 3-D GPR surveys can be used to map subglacial channels near temperate alpine glacier margins.

Subglacial topography is an important feature in numerous ice-sheet analyses and can drive the routing of water at the bed. Bed topography is primarily measured with ice-penetrating radar. Significant gaps, however, remain in data coverage that require interpolation. Topographic interpolations are typically made with kriging, as well as with mass conservation, where ice flow dynamics are used to constrain bed geometry. However, these techniques generate bed topography that is unrealistically smooth at small scales, which biases subglacial water flowpath models and makes it difficult to rigorously quantify uncertainty in subglacial drainage patterns. To address this challenge, we adapt a geostatistical simulation method with probabilistic modeling to stochastically simulate bed topography such that the interpolated topography retains the spatial statistics of the ice-penetrating radar data. We use this method to simulate subglacial topography using mass conservation topography as a secondary constraint. We apply a water routing model to each of these realizations. Our results show that many of the flowpaths significantly change with each topographic realization, demonstrating that geostatistical simulation can be useful for assessing confidence in subglacial flowpaths.

The cross-sectional profile of a glacial valley can be obtained with a variational principle in which the friction against the valley walls and the glacier bed is extremized, subject to a Lagrangian constraint. We show that the actual valley profile maximizes the friction, thus settling an old debate.

The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component for the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The geometric boundary conditions and in particular the nature of the subglacial bed of the NEGIS are essential to understand its ice flow dynamics. We present a record of more than 8000 km of radar survey lines of multi-channel, ultra-wideband radio echo sounding data covering an area of 24 000 km2, centered on the drill site for the East Greenland Ice-core Project (EGRIP), in the upper part of the NEGIS catchment. Our data yield a new detailed model of ice-thickness distribution and basal topography in the region. The enhanced resolution of our bed topography model shows features which we interpret to be caused by erosional activity, potentially over several glacial–interglacial cycles. Off-nadir reflections from the ice–bed interface in the center of the ice stream indicate a streamlined bed with elongated subglacial landforms. Our new bed topography model will help to improve the basal boundary conditions of NEGIS prescribed for ice flow models and thus foster an improved understanding of the ice-dynamic setting.

Glacier thinning and retreat drives initial acceleration of glacier sliding and erosion, de-buttressing of steep valley walls, and destabilization of ice-marginal deposits and bedrock, which can lead to massive rock avalanching and accelerated incision of tributary watersheds. A compelling example of these changes occurred in Taan Fjord in SE Alaska due to the rapid thinning and retreat of Tyndall Glacier over the past half century. Increased glacier sliding speeds led to both increased rates of subglacial erosion and the evacuation of subglacially stored sediments into the proglacial basins. The shrinking glacier also exposed proglacial tributary watersheds to rapid incision and denudation driven by >350 m of baselevel fall in a few decades. Moreover, in October 2015 a large tsunamigenic landslide occurred at the terminus of Tyndall Glacier, largely due to thinning exposing oversteepened, unstable slopes. Sediment yields from the glacier, the landslide and the tributary watersheds, measured from surveys of the sediments in the fjord collected in 1999 and 2016, are compared to ongoing changes in glacier and fjord geometry to investigate the magnitude of glacial and paraglacial denudation in Taan Fjord during retreat. In the last 50 years, sediment yields from the glacier and non-glacial tributaries kept pace with the rapid rate of retreat, and were on par with each other. Notably, basin-averaged erosion rates from the paraglacial landscape were twice that from the glacier, averaging 58 ± 9 and 26 ± 5 mm a−1, respectively. The sharp increases in sediment yields during retreat observed from both the glacier and the adjacent watersheds, including the landslide, highlight the rapid evolution of landscapes undergoing glacier shrinkage.

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.