Ice-bed Interface Research Articles

Multi-scale variations of subglacial hydro-mechanical conditions at Kongsvegen glacier, Svalbard

Abstract. The flow of glaciers is largely controlled by changes at the ice–bed interface, where basal slip and sediment deformation drive basal glacier motion. Determining subglacial conditions and their responses to hydraulic forcing remains challenging due to the difficulty of accessing the glacier bed. Here, we monitor the interplay between surface runoff and hydro-mechanical conditions at the base of the Kongsvegen glacier in Svalbard. From July 2021 to August 2022, we measured both subglacial water pressure and till strength. Additionally, we derived median values of subglacial hydraulic gradient and radius of channelized subglacial drainage system from seismic power, recorded at the glacier surface. To characterize the variations in the subglacial conditions caused by changes in surface runoff, we investigate the variations of the following hydro-mechanical properties: measured water pressure, measured sediment ploughing forces, and derived hydraulic gradient and radius, over seasonal, multi-day, and diurnal timescales. We discuss our results in light of existing theories of subglacial hydrology and till mechanics to describe subglacial conditions. We find that during the short, low-melt-rate season in 2021, the subglacial drainage system evolved at equilibrium with runoff, increasing its capacity as the melt season progressed. In contrast, during the long and high-melt-rate season in 2022, the subglacial drainage system evolved transiently to respond to the abrupt and large water supply. We suggest that in the latter configuration, the drainage capacity of the preferential drainage axis was exceeded, promoting the expansion of hydraulically connected regions and local weakening of ice–bed coupling and, hence, enhanced sliding.

Evaluating different geothermal heat-flow maps as basal boundary conditions during spin-up of the Greenland ice sheet

Abstract. There is currently poor scientific agreement on whether the ice–bed interface is frozen or thawed beneath approximately one third of the Greenland ice sheet. This disagreement in basal thermal state results, at least partly, from differences in the subglacial geothermal heat-flow basal boundary condition used in different ice-flow models. Here, we employ seven widely used Greenland geothermal heat-flow maps in 10 000-year spin-ups of the Community Ice Sheet Model (CISM). We perform two spin-ups: one nudged toward thickness observations and the other unconstrained. Across the seven heat-flow maps, and regardless of unconstrained or nudged spin-up, the spread in basal ice temperatures exceeds 10 ∘C over large areas of the ice–bed interface. For a given heat-flow map, the thawed-bed ice-sheet area is consistently larger under unconstrained spin-ups than nudged spin-ups. Under the unconstrained spin-up, thawed-bed area ranges from 33.5 % to 60.0 % across the seven heat-flow maps. Perhaps counterintuitively, the highest iceberg calving fluxes are associated with the lowest heat flows (and vice versa) for both unconstrained and nudged spin-ups. These results highlight the direct, and non-trivial, influence of the heat-flow boundary condition on the simulated equilibrium thermal state of the ice sheet. We suggest that future ice-flow model intercomparisons should employ a range of basal heat-flow maps, and limit direct intercomparisons with simulations using a common heat-flow map.

Modes of Antarctic tidal grounding line migration revealed by Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) laser altimetry

Abstract. Tide-forced short-term migration of the grounding line (GL) of Antarctic ice shelves can impact ice dynamics at the ice sheet margins and obscures assessments of long-term GL advance or retreat. However, the magnitude of tidally induced GL migration is poorly known, and the spatial patterns and modes of variability are not well characterised. Here we develop and apply a technique that uses Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) repeat-track laser altimetry to locate the inland limit of tidal ice shelf flexure for each sampled tide, enabling the magnitude and temporal variability of tidal GL migration to be resolved. We demonstrate its application at an ice plain north of Bungenstockrücken, in a region of the southern Ronne Ice Shelf subject to large ocean tides. We observe a 1300 km2 area of ephemeral grounding over which the GL migrates by up to 15 km between low and high tide and identify four distinct modes of migration: linear, asymmetric, threshold and hysteresis. The short-term movement of the GL dominates any long-term migration signal in this location, and the distribution of GL positions and modes contains information about spatial variability in the ice–bed interface. We discuss the impact of extreme tidal GL migration on ice shelf–ocean–subglacial systems in Antarctica and make recommendations for how GLs should be more precisely defined and documented in future by the community.

Design and performance of the Hotrod melt-tip ice-drilling system

Abstract. We introduce the design and performance of an electrothermal ice-drilling system designed to insert a temperature sensor cable into ice. The melt tip is relatively simple and low-cost, designed for a one-way trip to the ice–bed interface. The drilling system consists of a melt tip, umbilical cable, winch, interface, power supply, and support items. The melt tip and the winch are the most novel elements of the drilling system, and we make the hardware and electrical designs of these components available open-access. Tests conducted in a laboratory indicate that the melt tip has an electrical energy to forward melting heat transfer efficiency of ∼35 % with a theoretical maximum penetration rate of ∼12 m h−1at maximum 6.0 kW power. In contrast, ice-sheet testing suggests the melt tip has an analogous heat transfer efficiency of ∼15 % with a theoretical maximum penetration rate of ∼6 m h−1. We expect the efficiency gap between laboratory and field performance to decrease with increasing operator experience. Umbilical freeze-in due to borehole refreezing is the primary depth-limiting factor of the drilling system. Enthalpy-based borehole refreezing assessments predict refreezing below critical umbilical diameter in ∼4 h at −20 ∘C ice temperatures and ∼20 h at −2 ∘C. This corresponds to a theoretical depth limit of up to ∼200 m, depending on firn thickness, ice temperature, and operator experience.

Highly variable friction and slip observed at Antarctic ice stream bed

The slip of glaciers over the underlying bed is the dominant mechanism governing the migration of ice from land into the oceans, with accelerating slip contributing to sea-level rise. Yet glacier slip remains poorly understood, and observational constraints are sparse. Here we use passive seismic observations to measure both frictional shear stress and slip at the bed of the Rutford Ice Stream in Antarctica using 100,000 repetitive stick-slip icequakes. We find that basal shear stresses and slip rates vary from 104 to 107 Pa and 0.2 to 1.5 m per day, respectively. Friction and slip vary temporally over the order of hours, and spatially over 10s of metres, due to corresponding variations in effective normal stress and ice–bed interface material. Our findings suggest that the bed is substantially more complex than currently assumed in ice stream models and that basal effective normal stresses may be significantly higher than previously thought. Our observations can provide constraints on the basal boundary conditions for ice-dynamics models. This is critical for constraining the primary contribution of ice mass loss in Antarctica and hence for reducing uncertainty in sea-level rise projections.

Reconstructing terrestrial ice sheet retreat dynamics from hummocky topography using multiscale evidence: An example from central Ireland

The research reported here combines high-resolution digital elevation models (DEMs) derived from airborne LiDAR with ground-penetrating radar (GPR) surveys and macro- and micro-scale sedimentological analyses to study the geomorphology and internal structure and composition of ridges and mounds within an area of hummocky topography in the Brosna basin, central Ireland. Our evidence indicates that much of the hummocky topography consists of fragmented mega-scale glacial lineations (MSGLs) indicating a phase of accelerated ice flow, overlain by groups of small ridges composed of subglacially derived till and sediment gravity flow deposits. Subtle differences in ridge morphology indicate they may be multi-genetic and formed both ice-marginally as moraines and subglacially as small ribbed moraine. We interpret MSGL-ridge associations as a subglacial bedform continuum reflecting the evolution from a deforming bed to brittle deformation, due to changing thermal and hydrological conditions at the ice-bed interface during ice sheet retreat. Adjacent glaciofluvial landforms indicate initial en- and supra-glacial meltwater drainage, possibly directly related to changing basal thermal/hydrological characteristics and subglacial ridge formation. Subsequently a subglacial conduit system evolved, and later retreat involved formation of ice-marginal ridges. The partially preserved landsystem reflects formation at a polythermal to active temperate ice marginal zone, rather than the stagnating temperate margin previously assumed for this area. Our work demonstrates the usefulness of hummocky topography in identifying changes in ice sheet bed thermal/hydrological characteristics during deglaciation and the importance of combining multiple evidence strands in reconstructing the processes involved in glacial landform construction.

The impact of spatially varying ice sheet basal conditions on sliding at glacial time scales

AbstractSpatially variable basal conditions are thought to govern how ice sheets behave at glacial time scales (>1000 years) and responsible for changes in dynamics between the core and peripheral regions of the Laurentide and Fennoscandian ice sheets. Basal motion is accomplished via the deformation of unconsolidated sediments, or via sliding of the ice over an undeformable bed. We present an ice sheet sliding module for the Parallel Ice Sheet Model (PISM) that takes into account changes in sediment cover and incorporates surface meltwater. This model routes meltwater, produced at the surface and base of the ice sheet, toward the margin of the ice sheet. Basal sliding is accomplished through the deformation of water saturated sediments, or sliding at the ice-bed interface. In areas with continuous, water saturated sediments, sliding is almost always accomplished through sediment deformation. In areas with incomplete cover, sliding has a stronger dependence on the supply of water. We find that the addition of surface meltwater to the base is a more important factor for ice sheet evolution than the style of sliding. In a glacial cycle simulation, our model causes a more rapid buildup of the Laurentide Ice Sheet.

Numerical modelling of subglacial ribs, drumlins, herringbones, and mega‐scale glacial lineations reveals their developmental trajectories and transitions

AbstractInitially a matter of intellectual curiosity, but now important for understanding ice‐sheet dynamics, the formation of subglacial bedforms has been a subject of scientific enquiry for over a century. Here, we use a numerical model of the coupled flow of ice, water, and subglacial sediment to explore the formation of subglacial ribs (i.e., ribbed moraine), drumlins and mega‐scale glacial lineations (MSGLs). The model produces instabilities at the ice–bed interface, which result in landforms resembling subglacial ribs and drumlins. We find that a behavioural trajectory is present. Initially subglacial ribs form, which can either develop into fields of organized drumlins, or herringbone‐type structures misaligned with ice flow. We present potential examples of these misaligned bedforms in deglaciated landscapes, the presence of which means caution should be taken when interpreting cross‐cutting bedforms to reconstruct ice flow directions. Under unvarying ice flow parameters, MSGLs failed to appear in our experiments. However, drumlin fields can elongate into MSGLs in our model if low ice–bed coupling conditions are imposed. The conditions under which drumlins elongate into MSGLs are analogous to those found beneath contemporary ice streams, providing the first mechanism, rather than just an association, for linking MSGLs with ice stream flow. We conclude that the instability theory, as realized in this numerical model, is sufficient to explain the fundamental mechanics and process‐interactions that lead to the initiation of subglacial bedforms, the development of the distinctive types of bedform patterns, and their evolutionary trajectories. We therefore suggest that the first part of the longstanding ‘drumlin problem’ – how and why they come into existence – is now solved. However, much remains to be discovered regarding the exact sedimentary and hydrological processes involved.

What can deep ice, water, sediments, and bedrock at the ice–bed interface tell us?

What can deep ice, water, sediments, and bedrock at the ice–bed interface tell us?

Inverting ice surface elevation and velocity for bed topography and slipperiness beneath Thwaites Glacier

Abstract. There is significant uncertainty over how ice sheets and glaciers will respond to rising global temperatures. Limited knowledge of the topography and rheology of the ice–bed interface is a key cause of this uncertainty as models show that small changes in the bed can have a large influence on predicted rates of ice loss. Most of our detailed knowledge of bed topography comes from airborne and ground-penetrating radar observations. However, these direct observations are not spaced closely enough to meet the requirements of ice-sheet models, so interpolation and inversion methods are used to fill in the gaps. Here we present the results of a new inversion of surface elevation and velocity data over Thwaites Glacier, West Antarctica, for bed topography and slipperiness (i.e. the degree of basal slip for a given level of drag). The inversion is based on a steady-state linear perturbation analysis of the shallow-ice-stream equations. The method works by identifying disturbances to surface flow which are caused by obstacles or sticky patches in the bed and can therefore be applied wherever the shallow-ice-stream equations hold and where surface data are available, even where the ice thickness is not well known. We assess the performance of the inversion for topography with the available radar data. Although the topographic output from the inversion is less successful where the bed slopes steeply, it compares well with radar data from the central trunk of the glacier for medium-wavelength features (5–50 km). This method could therefore be useful as an independent test of other interpolation methods such as mass conservation and kriging. We do not have data to allow us to assess the success of the slipperiness results from our inversions, but we provide maps that may guide future seismic data collection across Thwaites Glacier. The methods presented here show significant promise for using high-resolution satellite datasets, calibrated by sparser field datasets, to generate high-resolution bed topography products across the ice sheets and therefore contribute to reduced uncertainty in predictions of future sea-level rise.

Hysteretic evolution of ice rises and ice rumples in response to variations in sea level

Abstract. Ice rises and ice rumples are locally grounded features found in coastal Antarctica and are surrounded by otherwise freely floating ice shelves. An ice rise has an independent flow regime, whereas the flow regime of an ice rumple conforms to that of the ice shelf and merely slows the flow of ice. In both cases, local highs in the bathymetry are in contact with the ice shelf from below, thereby regulating the large-scale ice flow, with implications for the upstream continental grounding line position. This buttressing effect, paired with the suitability of ice rises as a climate archive, necessitates a better understanding of the transition between ice rise and ice rumple, their evolution in response to a change in sea level, and their dynamic interaction with the surrounding ice shelf. We investigate this behaviour using a three-dimensional full Stokes ice flow model with idealised ice rises and ice rumples. The simulations span end-member basal friction scenarios of almost stagnant and fully sliding ice at the ice–bed interface. We analyse the coupling with the surrounding ice shelf by comparing the deviations between the non-local full Stokes surface velocities and the local shallow ice approximation (SIA). Deviations are generally high at the ice divides and small on the lee sides. On the stoss side, where ice rise and ice shelf have opposing flow directions, deviations can be significant. Differences are negligible in the absence of basal sliding where the corresponding steady-state ice rise is larger and develops a fully independent flow regime that is well described by SIA. When sea level is increased, and a transition from ice rise to ice rumple is approached, the divide migration is more abrupt the higher the basal friction. In each scenario, the transition occurs after the stoss-side grounding line has moved over the bed high and is positioned on a retrograde slope. We identify a hysteretic response of ice rises and ice rumples to changes in sea level, with grounded area being larger in a sea-level-increase scenario than in a sea-level-decrease scenario. This hysteresis shows not only irreversibility following an equal increase and subsequent decrease in sea level but also that the perturbation history is important when the ice rise or ice rumple geometry is not known. The initial grounded area needs to be carefully considered, as this will determine the formation of either an ice rise or an ice rumple, thereby causing different buttressing effects.

High-resolution subglacial topography around Dome Fuji, Antarctica, based on ground-based radar surveys over 30 years

Abstract. The retrieval of continuous ice core records of more than 1 Myr is an important challenge in palaeo-climatology. For identifying suitable sites for drilling such ice, knowledge of the subglacial topography and englacial layering is crucial. For this purpose, extensive ground-based ice radar surveys were carried out over Dome Fuji in the East Antarctic plateau during the 2017/18 and 2018/19 austral summers by the Japanese Antarctic Research Expedition, on the basis of ground-based radar surveys conducted over the previous ∼ 30 years. High-gain Yagi antennae were used to improve the antenna beam directivity, thereby significantly decreasing hyperbolic features of unfocused along-track diffraction hyperbolae in the echoes from mountainous ice–bedrock interfaces. We combined the new ice thickness data with the previous ground-based data, recorded since the 1980s, to generate an accurate high-spatial-resolution (up to 0.5 km between survey lines) ice thickness map. This map revealed a complex landscape composed of networks of subglacial valleys and highlands. Based on the new map, we examined the roughness of the ice–bed interface, the bed surface slope, the driving stress of ice and the subglacial hydrological condition. These new products and analyses set substantial constraints on identifying possible locations for new drilling. In addition, our map was compared with a few bed maps compiled by earlier independent efforts based on airborne radar data to examine the difference in features between datasets. Our analysis suggests that widely available bed topography products should be validated with in situ observations where possible.

Estimation of surface melt induced melt pond depths over Amery Ice Shelf, East Antarctica using Multispectral and ICESat-2 data

Surface melting is increasing over the Antarctica continent, with unknown impacts on glacier and ice sheet dynamics at the ice-bed interface. Surface meltwater drains to the bed of Antarctica ice sheets and outlet glaciers, causing accelerations in ice flow (up to 100% greater than the annual mean). Because melt occurrences are currently erratic (also Antarctica surface melt refreezes), effective subglacial drainage cannot be maintained, resulting in repeated short-lived ice flow disturbances and calving of ice shelves. Longterm exposure to strong surface melting promotes the formation of melt ponds and subsequently Supraglacial lakes and therefore there is possibility of hydrofracturing. The present study aims to estimate the volume of surface melt-induced melt ponds over the Amery Ice Shelf (AIS) using a physically-based radiative transfer model. The AIS, located on Antarctica's east coast, is one of the world's biggest glacier drainage basins. The ICESat-2 geolocated photon data were used to validate the depth estimates of melt ponds derived from Landsat-8 and Sentinel-2 multispectral datasets based on the model. The average depth estimated is around 1.9±0.6 m for the melt ponds adopted for the present study. The present study aims to quantify surface melt over the Antarctica continent (AIS region) in the context of India’s scientific credentials in polar science. We recommend monitoring extensive surface melt and subsurface refreezing which might result in more ponding above impermeable ice layers in the coming years.

The Impact of Basal Roughness on Inland Thwaites Glacier Sliding.

Swath radar technology enables three‐dimensional mapping of modern glacier beds over large areas at resolutions that are higher than those typically used in ice‐flow models. These data may enable new understanding of processes at the ice‐bed interface. Here, we use two densely surveyed swath‐mapped topographies (<50 m2 resolution) of Thwaites Glacier to investigate the sensitivity of inferred basal friction proxies to bed roughness magnitude and orientation. Our work suggests that along‐flow roughness influences inferred friction more than transverse‐flow roughness, which agrees with analytic form‐drag sliding theory. Using our model results, we calculate the slip length (the ratio of internal shear to basal slip). We find excellent agreement between the numerically derived slip lengths and slip lengths predicted by analytic form‐drag sliding theory, which suggests that unresolved short wavelength bed roughness may control sliding in the Thwaites interior.

A dynamic saline groundwater system mapped beneath an Antarctic ice stream

Antarctica's fast-flowing ice streams drain the ice sheet, with their velocity modulated by subglacial water systems. Current knowledge of these water systems is limited to the shallow portions near the ice-bed interface, but hypothesized deeper groundwater could also influence ice streaming. Here, we use magnetotelluric and passive seismic data from Whillans Ice Stream, West Antarctica, to provide the first observations of deep sub-ice stream groundwater. Our data reveal a volume of groundwater within a >1-kilometer-thick sedimentary basin that is more than an order of magnitude larger than the known subglacial system. A vertical salinity gradient indicates exchange between paleo seawater at depth and contemporary basal meltwater above. Our results provide new constraints for subglacial water systems that affect ice streaming and subglacial biogeochemical processes.

Streamlined subglacial bedform sensitivity to bed characteristics across the deglaciated Northern Hemisphere

AbstractStreamlined subglacial bedforms observed in deglaciated landscapes provide the opportunity to assess the sensitivity of glacier dynamics to bed characteristics across broader spatiotemporal scales than is possible for contemporary glacial systems. While many studies of streamlined subglacial bedforms rely on manual mapping and qualitative (i.e., visual) assessment, we semi‐automatically identify 11,628 sedimentary and bedrock bedforms, created during and following the Last Glacial Maximum across nine geologically and topographically diverse deglaciated sites in the Northern Hemisphere. Using this large dataset of landforms and associated morphometrics, we empirically test the importance of subglacial terrain on bedform morphology and ice‐flow behavior. A minimum bedform length–width ratio threshold provides a constraint on minimum morphometrics needed for streamlined bedforms to develop. Similarities in bedform metric distribution regardless of bed properties indicate that all bed types may support similar distributions of warm‐based ice flow conditions. Ice flow within valleys with easily erodible beds host the most elongate bedforms yet the widest range in bedform elongation and bedform surface relief. The presence of these highly elongate bedforms suggest high ice‐flow velocities occur within valley settings despite spatially heterogeneous landform‐generating processes. In contrast, lithified sedimentary beds within regions not constrained by topography on the scale of 1–102 km contain bedforms with high density and packing, low change in surface relief and low elongation, indicating spatially uniform and organized interactions at the ice–bed interface and consistency in ice‐flow velocity. Regardless of genesis, we find a sensitivity of bedform elongation (i.e., used to interpret ice‐flow speed or persistence) to topographic conditions on the scale of 1–102 km, while bedform density is sensitive to bed lithology. The findings presented in this study provide analogues for processes of subglacial erosion and deposition, ice–bed interactions and warm‐based ice flow within contemporary glacial systems.

An augmented lagrangian algorithm for recovery of ice thickness in unidirectional flow using the shallow ice approximation

A key parameter in ice flow modelling is basal slipping at the ice-bed interface as it can have a large effect on the resultant ice thickness. Unfortunately, its contribution to surface observations can be hard to distinguish from that of bed undulations. Therefore, inferring the ice thickness from surface measurements is an interesting and non-trivial inverse problem. This paper presents a method for recovering dually the ice thickness, in turn the bed topography, and the basal slip using only surface elevation and speed measurements. The unidirectional shallow ice approximation is first implemented to model steady state ice flow for given bedrock and basal slip profiles. This surface is then taken as synthetic observed data. To invert the data, an optimisation approach is taken, employing a basic Uzawa type algorithm to minimise the augmented Lagrangian objective functional. This optimisation recovers the diffusion coefficient which gives the best fit to observations. Combining this recovered diffusion with observed surface velocity, a simple Newton’s method is used to recover both the ice thickness and basal slip. The method was successful in each test case and this implies that it should be possible to recover both of these parameters in two-dimensional cases also.

Glaciotectonic disintegration of roches moutonnées during glacial ripping in east Sweden

ABSTRACT Roches moutonnées are typical landforms of glacial erosion developed in hard rocks, with an asymmetric profile caused by abrasion and lee-side plucking. In eastern Sweden, some roches moutonnées show extensive damage, including open fractures, disintegration into blocks, fracture caves and short boulder trains. Disintegration increases along ice-flow directions during deglaciation of the last Weichselian Fennoscandian Ice Sheet, indicating a subglacial origin: limited edge rounding can be explained by a combination of hard rock, slow abrasion rates and disintegration just prior to deglaciation. The roches moutonnées initially developed in kernels of gneissic rocks with a wide fracture spacing (large block size) and interlocking fracture pattern, and hence high overall rock mass strength. Dilated fractures and ‘fracture caves’ occur up to 15 m below the ice-bed interface. It is proposed that hydraulic jacking by overpressured water opened up the rock mass along pre-existing fractures. Jacking reduced rock mass strength, allowing glaciotectonic deformation of the roches moutonnées. Uneven hydraulic jacking led to uplift of individual fracture-bound blocks above the pre-existing smooth, abraded surface of the roches moutonnées, creating blunt, step-like edges. These edges allowed high ice pushing forces to act on large blocks: where blocks extend into the deeper rock mass, they further aided the disintegration of the rock mass. The disintegrated roches moutonnées can be regarded as transient features between intact bedrock and complete disintegration into boulders. The jacking-disintegration-transport sequence is characteristic of glacial ripping and very different from classic lee-side plucking.

Glacial processes and landforms

AbstractFrom 1965 to 2000 glacial geomorphology became increasingly specialized and developed significantly due to technological improvements, particularly in remote sensing, surveying and field-based glaciological process studies. The better understanding of basal thermal regimes in ice sheets and glaciers led to the development of concepts such as spatial and temporal migration of ice divides in dynamic ice sheets that could overprint subglacial landform assemblages, debris entrainment processes related to polythermal glacier systems, and glacier and ice sheet beds composed of cold and warm based mosaics. Process observations at the ice–bed interface led to the discovery of the third glacier flow mechanism, substrate deformation, which provided the impetus to reconstruct the genesis of subglacial bedforms such as drumlins and to evaluate the origins and potential flow law for till. Numerical evaluations of glacial erosion led to a better understanding of abrasion and quarrying as well as the erection of genetic models and erosion rates for larger-scale features such as U-shaped valleys and cirques. Linkages were made between debris transport pathways and moraine construction in supraglacial environments, with the role of glacier structure being linked to specific landforms, such as medial, lateral, hummocky and ice-cored moraines as well as rock glaciers. Our appreciation of the erosional and depositional impacts of glacifluvial systems was enhanced significantly with the advent of process observations on the hydrology of modern glaciers as well as the final vindication of J.H. Bretz and his proposed jökulhlaup origins of the Channelled Scablands and the Missoula Floods. In addition to the increasing numbers of studies at modern glacier snouts, the embracing of sedimentology by glacial geomorphologists was to result in significant developments in understanding the process-form regimes of subglacial, marginal and proglacial landforms, particularly the recognition of landform continua and hybrids. Advances resulting from this included the recognition of different modes of moraine and glacitectonic thrust mass development, lithofacies models of the varied glacifluvial depositional environments, and the initial expansion of work on the sediments and depocentres of glacimarine settings, the latter being the result of glacial research taking to submersibles and ice-strengthened ships for the first time. A similarly new frontier was the expansion of research on the increasingly higher resolution images returning from Mars, where extraterrestrial glaciations were recognized based on comparisons with Earth analogues. Holistic appreciations of glaciation signatures using landform assemblages were developed, initially as process-form models and later as glacial landsystems, providing an ever-expanding set of templates for reconstructing palaeoglaciology in the wide variety of topographic and environmental settings, which also acknowledge spatial and temporal change in glacier and ice sheet systems.

Seasonal Cold-Wave Propagation Into the Near-Surface Ice of Debris-Covered Khumbu Glacier, Nepal

Meltwater from high-elevation debris-covered glaciers—particularly those located in the greater Himalaya and Andes—shapes the water supply of major rivers and nourishes substantial terrestrial, estuarine, and marine habitats (Kraaijenbrink et al., 2017; Immerzeel et al., 2020). However, the relative inaccessibility and high elevation of such glaciers results in a paucity of data relating to their fundamental physical properties and processes, limiting the information available to constrain and evaluate numerical models of their behaviour and project future change. Knowledge of the subsurface properties of such glaciers is particularly deficient because it is largely obscured to satellite and airborne remote sensing; englacial investigations therefore commonly require direct access (Miles et al., 2020). Of the physical properties of glaciers, ice temperature exerts an important control over glaciological processes, such as glacier motion, and their modelled behaviour. For example, ice viscosity is sensitive to temperature such that, under the same stress, ice approaching the melting point deforms 5–10 times more rapidly than it would at −10°C (Deeley and Woodward, 1908; Cuffey and Paterson, 2010). Basal motion depends on lubrication facilitated by the presence of meltwater at the ice-bed interface and/or within the pore space of a subglacial sediment layer. Measurements of near-surface ice temperatures are important for modelling the surface energy balance and projecting the future mass-balance response of glaciers to anticipated climate change. This is especially the case for glaciers with a thick supraglacial debris layer that insulates the underlying ice (according to debris layer thickness and lithology), reducing ablation and potentially extending glacier longevity (Nicholson and Benn, 2006; Nicholson and Benn, 2013; Anderson and Anderson, 2016). Here, we present a one-year time series of near-surface ice temperatures, measured between 1.5 and 7.0 m below the ice surface, in a borehole drilled by hot water into the debris-covered tongue of Khumbu Glacier, Nepal.