‘Detachment’ of icefield outlet glaciers – catastrophic thinning and retreat of the Columbia Glacier (Canada)

We present an investigation of changes taking place on the Columbia Glacier – a lake-terminating outlet of the Columbia Icefield in the Canadian Rockies. The Columbia Icefield is the largest, and one of the most important, ice bodies in the Canadian Rockies. Like other ice masses, it stores water as snow and ice during the winter and releases it during warmer summer months, sustaining river flows and the ecosystems that rely on them. However, the Columbia Glacier and Icefield is shrinking. We use Landsat and Sentinel-2 imagery to show that the Columbia Glacier has retreated increasingly rapidly in recent years, and suggest that this looks set to continue. Importantly, we identify a previously undocumented process that appears to be playing an important role in the retreat of this glacier. This process involves the ‘detachment’ of the glacier tongue from its accumulation area in the Columbia Icefield. This process is important because the tongue is cut off from the accumulation area and there is no replenishment of ice that melts in the glacier's ablation area by flow from upglacier. As a consequence, for a given rate of ablation, the ice in the tongue will disappear much faster than it would if the local mass loss by melting/calving was partly offset by mass input by glacier flow. Such a change would alter the relationship between rates of surface melting and rates of glacier frontal retreat. We provide evidence that detachment has already occurred elsewhere on the Columbia Icefield and that it is likely to affect other outlet glaciers in the future. Modelling studies forecast this detachment activity, which ultimately results in a smaller ‘perched’ icefield without active outlets

The internal layering of Pine Island Glacier, West Antarctica, from airborne radar-sounding data

This paper presents an overview of internal layering across Pine Island Glacier, West Antarctica, as measured from airborne-radar data acquired during a survey conducted by the British Antarctic Survey and the University of Texas in the 2004/05 season. Internal layering is classified according to type (continuous/discontinuous/missing) and the results compared with InSAR velocities. Several areas exhibit disruption of internal layers that is most likely caused by large basal shear stresses. Signs of changes in flow were identified in a few inter-tributary areas, but overall the layering classification and distribution of layers indicate that only minor changes in ice-flow regime have taken place. This is supported by bed-topography data that show the main trunk of the glacier, as well as some of the tributaries, are topographically controlled and located in deep basins

Inland extent of the Weddell Sea Rift imaged by new aerogeophysical data

Peer reviewedPostprin

Subglacial controls on the flow of Institute Ice Stream, West Antarctica

The Institute Ice Stream (IIS) rests on a reverse-sloping bed, extending >150 km upstream into the ~1.8 km deep Robin Subglacial Basin, placing it at the threshold of marine ice-sheet instability. Understanding IIS vulnerability has focused on the effect of grounding-line melting, which is forecast to increase significantly this century. Changes to ice-flow dynamics are also important to IIS stability, yet little is known about them. Here we reveal the trunk of the IIS occurs downstream of the intersection of three discrete subglacial features; a large ‘active’ subglacial lake, a newly-discovered sharp transition to a zone of weak basal sediments, and a major tectonic rift. The border of IIS trunk flow is confined by the sediment on one side, and by a transition between basal melting and freezing at the border with the Bungenstock Ice Rise. By showing how basal sediment and water dictate present-day flow of IIS, we reveal that ice-sheet stability here is dependent on this unusual arrangement

Using bed-roughness signatures to characterise glacial landform assemblages beneath contemporary and palaeo ice-sheets

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

Seasonal evolution of the supraglacial drainage network at Humboldt Glacier, North Greenland, between 2016 and 2020

Supraglacial rivers and lakes are important for the routing and storage of surface meltwater during the summer melt season across the Greenland Ice Sheet (GrIS), yet remain poorly mapped and quantified across the northern part of the ice sheet, which is rapidly losing mass. Here we produce, for the first time, a high-resolution record of the supraglacial drainage network (including both rivers and lakes) and its seasonal behaviour at Humboldt Glacier, a wide-outlet glacier draining a large hydrologic catchment (13,488 km2), spanning the period 2016 to 2020 using 10 m spatial resolution Sentinel-2 imagery. Our results reveal a perennially extensive yet interannually-variable supraglacial network extending from an elevation of 200 m a.s.l to a maximum of ~1440 m a.s.l recorded in 2020, with limited development of the network observed in the low melt years of 2017 and 2018. The supraglacial drainage network is shown to cover an area ranging between 965.7 km2 (2018) and 1566.3 km2 (2019) at its maximum seasonal extent, with spatial coverage of up to 2685 km2 recorded during the early phases of the melt season when a slush zone is most prominent. Up-glacier expansion and the development of an efficient supraglacial drainage network as surface runoff increases and the snowline retreats is clearly visible. Preconditioning of the ice surface following a high melt year is also observed, with the earlier widespread exposure of the supraglacial drainage network in 2020 compared to other years; a finding that may become representative with persistent warmer years into the future. Overall, this study provides evidence of a persistent, yet dynamic, supraglacial drainage network at this prominent northern GrIS outlet glacier and advances our understanding of such hydrologic processes, particularly under ongoing climatic warming and enhanced runoff

Aerodynamic roughness of glacial ice surfaces derived from high-resolution topographic data

This paper presents new methods of estimating the aerodynamic roughness (z0) of glacier ice directly from three-dimensional point clouds and digital elevation models (DEMs), examines temporal variability of z0, and presents the first fully distributed map of z0 estimates across the ablation zone of an Arctic glacier. The aerodynamic roughness of glacier ice surfaces is an important component of energy balance models and meltwater runoff estimates through its influence on turbulent fluxes of latent and sensible heat. In a warming climate these fluxes are predicted to become more significant in contributing to overall melt volumes. Ice z0 is commonly estimated from measurements of ice surface microtopography, typically from topographic profiles taken perpendicular to the prevailing wind direction. Recent advances in surveying permit rapid acquisition of high-resolution topographic data allowing revision of assumptions underlying conventional z0 measurement. Using Structure from Motion (SfM) photogrammetry with Multi-View Stereo (MVS) to survey ice surfaces with millimeter-scale accuracy, z0 variation over 3 orders of magnitude was observed. Different surface types demonstrated different temporal trajectories in z0 through 3 days of intense melt. A glacier-scale 2 m resolution DEM was obtained through terrestrial laser scanning (TLS), and subgrid roughness was significantly related to plot-scale z0. Thus, we show for the first time that glacier-scale TLS or SfM-MVS surveys can characterize z0 variability over a glacier surface potentially leading to distributed representations of z0 in surface energy balance models

Quantifying bed roughness beneath contemporary and palaeo-ice streams

Bed roughness is an important control on ice-stream location and dynamics. The majority of previous bed roughness studies have been based on data derived from radio echo sounding (RES) transects across Antarctica and Greenland. However, the wide spacing of RES transects means that the links between roughness and flow are poorly constrained. Here, we use Digital Terrain Model (DTM)/bathymetry data from a well-preserved palaeo-ice stream to investigate basal controls on the behaviour of contemporary ice streams. Artificial transects were set up across the Minch Palaeo-Ice Stream (NW Scotland) to mimic RES flight lines over Institute and Möller Ice Streams (Antarctica). We then explored how different data-resolution, transect orientation and spacing, and different methods, impact upon roughness measurements. Our results show that fast palaeo-ice flow can occur over a rough, hard bed, not just a smooth, soft bed, which much previous work has suggested. Smooth areas of the bed occur over both bedrock and sediment covered regions. Similar trends in bed roughness values were found using Fast Fourier Transform analysis and standard deviation methods. Smoothing of bed roughness results can hide important details. We propose that the typical spacing of RES transects is too wide to capture different landform assemblages, and that transect orientation influences bed roughness measurements in both contemporary and palaeo-ice-stream setting