The beds of many fast-flowing glaciers are underlain by sediment (till) whose physical properties regulate the mechanical and hydrological controls on glacier slip. Despite the importance of the till’s physical properties in affecting slip and hydrology, they are often ill-constrained, in part due to the inaccessible nature of glacier beds. Here, we assess a range of till properties for 11 different till samples, 10 of which were acquired from the Byrd Polar Rock Repository, collected from a range of locations along the Transantarctic Mountain Range, and one till sample collected from the paleo-grounding zone of Bindschadler Ice Stream (Ice Stream D). Using these sediment samples, we assess their coefficient of consolidation, friction angle, cohesion value, hydraulic conductivity, mineralogy, density and grain-size distribution. We provide an overview of the methods used and the resulting measurements but refrain from wide-scale interpretations of the results. In general, we find that some physical properties have a relatively small variance among the samples (e.g. friction angle, cohesion and grain-size maturity), whereas other properties vary more widely (e.g. hydraulic conductivity and coefficient of consolidation).

Water pressures in conduit networks (e.g. linked cavities, R-channels, subglacial water sheets) support a fraction of glacier weight; the remainder is carried by elevated normal stresses where basal ice is separated from the bed by extremely thin water films. We show how the requirements for melt flow toward conduits constrain film water pressure regimes. Where permeable bed materials facilitate efficient drainage, film pressures approach conduit water pressures, and deficits below the ice normal stress require the direct ice–bed stress transfer, termed thermomolecular pressure, that accompanies interfacial premelting. Beneath warm-based regions of the Greenland Ice Sheet, premelting is likely wherever meter-scale till depths are present; bulk melting conditions, with film water pressures balancing ice normal stresses, are anticipated where the ice and bed are both effectively impermeable. Englacial transport through vein networks can facilitate premelting over small film areas on impermeable beds; absent significant vein transport, such locations are prone to the formation of basal crevasses, which can buffer water pressures near the minimum compressive stress. Interfacial premelting is important for the preservation of unconsolidated bedforms, the strengthening of debris–bed frictional contacts and the enhancement of basal temperature fluctuations; anticorrelations between nearby borehole pressure records are expected when bulk melting prevails.

We located over one million icequakes that occurred during a 23 day period up-glacier of the grounding line of Rutford Ice Stream (RIS), West Antarctica. To investigate basal conditions and better understand the sliding process, we focused on sticky spots, spatially constrained clusters of regularly repeating basal icequakes. We developed a multistep workflow combining location and clustering to identify sticky spots with high confidence. The basal icequakes are consistent with low-angle slip in the ice-flow direction. Seismicity during our recording period showed weak tidal correlation and a minor sensitivity to atmospheric pressure, suggesting limited external forcing of basal processes. Icequake density was elevated on stoss sides of topographic highs likely characterized by exposed bedrock or compacted tills, and in low-lying moats eroded up-glacier of topographic highs, as well as in troughs extending down-glacier from the moats, likely reflecting hard beds that promote stick–slip. Local magnitudes followed a Gutenberg–Richter distribution with two distinct high (>2) b-values. Within sticky spots, repeating seismicity featured consistent waveforms and strong positive correlations between local magnitude and interevent time. During multihour active phases, seismic slip in sticky spots can fully accommodate ice flow, highlighting the role of stick–slip motion in basal sliding near the RIS grounding line. These results will help inform sliding relations used in numerical ice-flow models.

The basal thermal state of the Antarctic ice sheet strongly influences ice dynamics and mass balance. Although basal thermal states can be simulated using ice-sheet models, significant uncertainties persist, with poorly constrained geothermal heat flux being one important source of uncertainty. The substantial computational cost of ice-sheet models restricts large ensemble simulations essential for uncertainty quantification and sensitivity analyses. To address this limitation, we develop machine-learning emulators for rapidly simulating the basal temperature and melt rate of Totten Glacier. Three distinct machine-learning emulators were trained using outputs of full-Stokes model simulations with different geothermal heat flux fields, achieving accurate replication (R2 up to 0.92). At least five simulations are needed to train a reliable emulator for basal ice temperature. Model interpretability reveals complementary physics learning: Random forest and XGBoost prioritize local dynamic drivers (basal topography, ice velocity), while neural network captures lower-complexity statistical mapping that is more strongly influenced by globally predictive variables (surface temperature, ice thickness). For basal melt rate, surface velocity emerges as the primary driver. These findings highlight substantial potential for integrating machine learning with ice-sheet modeling to efficiently and better predict basal thermal states.

The basal thermal state of the Antarctic ice sheet (AIS)—whether the base is frozen or thawed—fundamentally underpins its flow and is an important factor in understanding its large-scale response to external forcings. Here, we present a first synthesis of the AIS basal thermal state combining two indirect and independent methods: (1) a compilation of nine three-dimensional thermomechanical simulations that calculate AIS basal temperature as part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) and (2) an estimate of the basal slip ratio, defined as the ratio of observed surface speed to deformational speed. This synthesis is evaluated against direct observations from deep boreholes and predicted flowpaths for water originating from subglacial lakes detected by altimetry and radar sounding. The synthesis predicts a thawed bed across most of West Antarctica and localized regions in East Antarctica. Most of the Antarctic Peninsula, the Transantarctic Mountains and several regions of East Antarctica are likely frozen at the bed. Overall, our synthesis suggests 46% of the AIS bed is likely thawed, 18% likely frozen and the remaining 36% is uncertain. Additional observations, particularly at the continental scale, are required to improve our understanding of Antarctica’s basal thermal state.

The aim of this paper is to characterise the internal structures and ice-flow history of representative valley glaciers in Svalbard and infer from them dynamic changes over centennial timescales. Three polythermal and one cold valley glacier are investigated using field- and laboratory-based techniques and remote sensing. Structures along flow-unit boundaries indicate that ice-flow configuration in three of the glaciers has remained stable spanning the residence time of the ice. Deformation of a flow-unit boundary in the fourth reveals an ice-flow instability, albeit one that has been maintained since its most recent advance. Macro-crystallographic, sedimentological and isotopic analyses indicate that basal ice is elevated to the glacier surface, as shown by entrained sediments and enrichment in heavy isotopes. In narrow zones of enhanced cumulative strain, new ice facies are generated through dynamic recrystallisation. The surface density of longitudinal foliation is shown to represent the relative magnitude of cumulative strain. Geometric similarities between flow-unit boundaries in Svalbard valley glaciers and larger scale longitudinal surface structures in ice sheets suggest that deformation mechanisms are common to both.

Small quantities of liquid water lining triple junctions in polycrystalline glacier ice form connected vein networks that enable material exchange with underlying basal environments. Diffuse debris concentrations commonly observed in ice marginal regions might be attributed to this mechanism. Following recent cryogenic ring-shear experiments, we observed emplacement along grain boundaries of loess particles several tens of microns in size. Here, we describe an idealized model of vein liquid flow to elucidate conditions favoring such particle transport. Gradients in liquid potential drive flow toward colder temperatures and lower solute concentrations, while deviations of the ice stress state from hydrostatic balance produce additional suction toward anomalously low ice pressures. Our model predicts particle entrainment following both modest warming along the basal interface resulting from anticipated natural changes in effective stress, and the interior relaxation of temperature and solute concentration imposed by our experimental protocols. Comparisons with experimental observations are encouraging, but suggest that liquid flow rates are somewhat higher and/or more effective at dragging larger particles than predicted by our idealized model with nominal parameter choices. Diffuse debris entrainment extending several meters above the glacier bed likely requires a more sophisticated treatment that incorporates effects of ice deformation or other processes.

We present new data from the debris-rich basal ice layers of the NEEM ice core (NW Greenland). Using mineralogical observations, SEM imagery, geochemical data from silicates (meteoric 10Be, εNd, 87Sr/86Sr) and organic material (C/N, δ13C), we characterize the source material, succession of previous glaciations and deglaciations and the paleoecological conditions during ice-free episodes. Meteoric 10Be data and grain features indicate that the ice sheet interacted with paleosols and eroded fresh bedrock, leading to mixing in these debris-rich ice layers. Our analysis also identifies four successive stages in NW Greenland: (1) initial preglacial conditions, (2) glacial advance 1, (3) glacial retreat and interglacial conditions and (4) glacial advance 2 (current ice-sheet development). C/N and δ13C data suggest that deglacial environments favored the development of tundra and taiga ecosystems. These two successive glacial fluctuations observed at NEEM are consistent with those identified from the Camp Century core basal sediments over the last 3 Ma. Further inland, GRIP and GISP2 summit sites have remained glaciated more continuously than the western margin, with less intense ice-substratum interactions than those observed at NEEM.

The geometry and thermal structure of western Greenland ice sheet are known to have undergone relatively substantial change over the Holocene. Evolution of the frozen and melted fractions of the bed associated with the ice-sheet retreat over this time frame remains unclear. We address this question using a thermo-mechanically coupled flowline model to simulate a 11 ka period of ice-sheet retreat in west central Greenland. Results indicate an episode of ~100 km of terminus retreat corresponded to ~16 km of upstream frozen/melted basal boundary migration. The majority of migration of the frozen area is associated with the enhancement of the frictional and strain heating fields, which are accentuated toward the retreating ice margin. The thermally active bedrock layer acts as a heat sink, tending to slow contraction of frozen-bed conditions. Since the bedrock heat flux in our region is relatively low compared to other regions of the ice sheet, the frozen region is relatively greater and therefore more susceptible to marginward changes in the frictional and strain heating fields. Migration of melted regions thus depends on both geometric changes and the antecedent thermal state of the bedrock and ice, both of which vary considerably around the ice sheet.

The importance of glacier sliding has motivated a rich literature describing the thermomechanical interactions between ice, liquid water and bed materials. Early recognition of the gradient in melting temperature across small bed obstacles led to focused studies of regelation. An appreciation for the limits on ice deformation rates downstream of larger obstacles highlighted a role for cavitation, which has subsequently gained prominence in descriptions of subglacial drainage. Here, we show that the changes in melting temperature that accompany changes in normal stress along a sliding ice interface near cavities and other macroscopic drainage elements cause appreciable supercooling and basal mass exchange. This provides the basis of a novel formation mechanism for widely observed laminated debris-rich basal ice layers.

Surface speeds in Greenland's ablation zone undergo substantial variability on an annual basis which are presumed to mainly be driven by changes in sliding. Yet, meltwater-forced changes in ice–bed coupling can also produce variable deformation motion, which impacts the magnitude of sliding changes inferred from surface measurements and provides important context to flow dynamics. We examine spatiotemporal changes in deformation, sliding and surface velocities over a 2-year period using GPS and a dense network of inclinometers installed in borehole grid drilled in western Greenland's ablation zone. We find time variations in deformation motion track sliding changes through the summer and entire measurement period. A distinct spatial deformation and sliding pattern is also observed within the borehole grid which remains similar during winter and summer flow. We suggest that positively covarying sliding and deformation across seasonal timescales is characteristic of passive areas that are coupled to regions undergoing transient forcing, and the spatial patterns are consistent with variations in the local bed topography. The covarying deformation and sliding result in a 1.5–17% overestimate of sliding changes during summer compared to that inferred from surface velocity changes alone. This suggests that summer sliding increases are likely overestimated in many locations across Greenland.

The Subglacial Antarctic Lakes Scientific Access (SALSA) Project accessed Mercer Subglacial Lake using environmentally clean hot-water drilling to examine interactions among ice, water, sediment, rock, microbes and carbon reservoirs within the lake water column and underlying sediments. A ~0.4 m diameter borehole was melted through 1087 m of ice and maintained over ~10 days, allowing observation of ice properties and collection of water and sediment with various tools. Over this period, SALSA collected: 60 L of lake water and 10 L of deep borehole water; microbes >0.2 μm in diameter from in situ filtration of ~100 L of lake water; 10 multicores 0.32–0.49 m long; 1.0 and 1.76 m long gravity cores; three conductivity–temperature–depth profiles of borehole and lake water; five discrete depth current meter measurements in the lake and images of ice, the lake water–ice interface and lake sediments. Temperature and conductivity data showed the hydrodynamic character of water mixing between the borehole and lake after entry. Models simulating melting of the ~6 m thick basal accreted ice layer imply that debris fall-out through the ~15 m water column to the lake sediments from borehole melting had little effect on the stratigraphy of surficial sediment cores.

Basal motion is the primary mechanism for ice flux in Greenland, yet a widely applicable model for predicting it remains elusive. This is due to the difficulty in both observing small-scale bed properties and predicting a time-varying water pressure on which basal motion putatively depends. We take a Bayesian approach to these problems by coupling models of ice dynamics and subglacial hydrology and conditioning on observations of surface velocity in southwestern Greenland to infer the posterior probability distributions for eight spatially and temporally constant parameters governing the behavior of both the sliding law and hydrologic model. Because the model is computationally expensive, characterization of these distributions using classical Markov Chain Monte Carlo sampling is intractable. We skirt this issue by training a neural network as a surrogate that approximates the model at a sliver of the computational cost. We find that surface velocity observations establish strong constraints on model parameters relative to a prior distribution and also elucidate correlations, while the model explains 60% of observed variance. However, we also find that several distinct configurations of the hydrologic system and stress regime are consistent with observations, underscoring the need for continued data collection and model development.

Basal motion of ice sheets depends in part on the roughness and material properties of the subglacial bed and the occurrence of water. To date, basal motion represents one of the largest uncertainties in ice-flow models. It is that component of the total flow velocity that can change most rapidly and can, therefore, facilitate rapid variations in dynamic behaviour. In this study, we investigate the subglacial properties of the East Antarctic Ice Sheet by statistically analysing the roughness of the bed topography, inferred from radio-echo sounding measurements. We analyse two sets of roughness parameters, one derived in the spatial and the other in the spectral domain, with two roughness parameters each. This enables us to compare the suitability of the four roughness parameters to classify the subglacial landscapes below the ice sheet. We further investigate the relationship of the roughness parameters with observed surface flow velocity and modelled basal temperatures of the ice sheet. We find that one of the roughness parameters, the Hurst exponent derived in the spatial domain, coincides with the thermal condition at the base of the ice sheet for slow flow velocities and varies with flow velocity.

Basal units – visibly distinct englacial structures near the ice-bed interface – warrant investigation for a number of reasons. Many are of unknown composition and origin, characteristics that could provide substantial insight into subglacial processes and ice-sheet history. Their significance, moreover, is not limited to near-bed depths; these units appear to dramatically influence the flow of surrounding ice. In order to enable improved characterization of these features, we develop and apply an algorithm that allows for the automatic detection of basal units. We use a tunable layer-optimized SAR processor to distinguish these structures from the bed, isochronous englacial layers and the ice-sheet surface, presenting a conceptual framework for the use of radio-echo character in the identification of ice-sheet features. We also outline a method by which our processor could be used to place observational constraints on basal units’ configuration, composition and provenance.

Subglacial rock friction is an important control on the sliding dynamics and erosive potential of hard-bedded glaciers, yet it remains largely unconstrained. To explore the relative influence of basal melt rate, effective stress and ice temperature on frictional resistance, we conducted abrasion experiments in which limestone beds were slid beneath a fixed slab of ice laden with granitic rock fragments. Shear stress scales linearly with melt rate and cryostatic stress, confirming that both viscous drag and effective stress are first-order controls on the contact force in drained conditions. Furthermore, temperature gradients in the ice increase the contribution of viscous drag on basal shear stress. In all experiments, the relationship between melt rate and shear stress is best explained by a model that accounts for the effects of regelation and viscous creep on the bed-normal drag force. We interpret this to mean fluid flow around entrained clasts contributed to basal drag even at subfreezing temperatures. Incorporating premelting dynamics into the Watts/Hallet model for subglacial rock friction, we find that the predicted debris-bed drag decreases by approximately an order of magnitude, with a corresponding ~3.5 × increase in the transition radius. This is lower than we observe for ice slightly below the pressure melting point.

Quarrying and abrasion are the two principal processes responsible for glacial erosion of bedrock. The morphologies of glacier hard beds depend on the relative effectiveness of these two processes, as abrasion tends to smooth bedrock surfaces and quarrying tends to roughen them. Here we analyze concentrations of bedrock discontinuities in the Tsanfleuron forefield, Switzerland, to help determine the geologic conditions that favor glacial quarrying over abrasion. Aerial discontinuity concentrations are measured from scaled drone-based photos where fractures and bedding planes in the bedrock are manually mapped. A Tukey honest significant difference test indicates that aerial concentration of bed-normal bedrock discontinuities is not significantly different between quarried and non-quarried areas of the forefield. Thus, an alternative explanation is needed to account for the spatial variability of quarried areas. To investigate the role that bed-parallel discontinuities might play in quarrying, we use a finite element model to simulate bed-normal fracture propagation within a stepped bed with different step heights. Results indicate that higher steps (larger spacing of bed-parallel discontinuities) propagate bed-normal fractures more readily than smaller steps. Thus, the spacing of bed-parallel discontinuities could exert strong control on quarrying by determining the rate that blocks can be loosened from the host rock.

The spatial distribution of basal water critically impacts the evolution of ice sheets. Current estimates of basal water distribution beneath the Greenland Ice Sheet (GrIS) contain large uncertainties due to poorly constrained boundary conditions, primarily from geothermal heat flux (GHF). The existing GHF models often contradict each other and implementing them in numerical ice-sheet models cannot reproduce the measured temperatures at ice core locations. Here we utilize two datasets of radar-detected basal water in Greenland to constrain the GHF at regions with a thawed bed. Using the three-dimensional ice-sheet model SICOPOLIS, we iteratively adjust the GHF to find the minimum GHF required to reach the bed to the pressure melting point, GHFpmp, at locations of radar-detected basal water. We identify parts of the central-east, south and northwest Greenland with significantly high GHFpmp. Conversely, we find that the majority of low-elevation regions of west Greenland and parts of northeast have very low GHFpmp. We compare the estimated constraints with the available GHF models for Greenland and show that GHF models often do not honor the estimated constraints. Our results highlight the need for community effort to reconcile the discrepancies between radar data, GHF models, and ice core information.

Johannes Weertman published his first glaciological paper in 1957 only 5 years after getting his DSc in metallurgy from the Carnegie Institute of Technology. The paper presented the very first sliding law developed quantitatively from first principles, and involved the unconventional idealization of bed roughness using cubic ‘tombstones’ of rock. Since 1957, there has been a great deal of progress in understanding glacier sliding, but few studies, if any, have preserved the original tombstone geometry that was a hallmark of this first theory. The current study presents a partial reanalysis of the sliding process over a bed with tombstone obstacles using modern numerical methods. The result confirms the enduring applicability of Weertman's model as a pedagogical tool and motivates new questions about (1) folding flow near bedrock obstacles that invert normal ice stratigraphy, (2) the presence and role of stress singularities on sharp edges of bedrock, and (3) the validity of a presumption that regelation flow can be plug-like.

Basal ice of glaciers and ice sheets frequently contains a well-developed stratification of distinct, semi-continuous, alternating layers of debris-poor and debris-rich ice. Here, the nature and distribution of shear within stratified basal ice are assessed through the anisotropy of magnetic susceptibility (AMS) of samples collected from Matanuska Glacier, Alaska. Generally, the AMS reveals consistent moderate-to-strong fabrics reflecting simple shear in the direction of ice flow; however, AMS is also dependent upon debris content and morphology. While sample anisotropy is statistically similar throughout the sampled section, debris-rich basal ice composed of semi-continuous mm-scale layers (the stratified facies) possesses well-defined triaxial to oblate fabrics reflecting shear in the direction of ice flow, whereas debris-poor ice containing mm-scale star-shaped silt aggregates (the suspended facies) possesses nearly isotropic fabrics. Thus, deformation within the stratified basal ice appears concentrated in debris-rich layers, likely the result of decreased crystal size and greater availability of unfrozen water associated with high debris content. These results suggest that variations in debris-content over small spatial scales influence ice rheology and deformation in the basal zone.