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.

The seasonal cycle of fast ice thickness in Prydz Bay, East Antarctica, was observed between March and December 2012. In March, we observed a 0.16 m thickness gain of 0.22 m-thick first-year ice (FYI), while 1.16 m-thick second-year ice (SYI) nearby simultaneously ablated by 0.59 m. A 1-D thermodynamic sea-ice model was applied to identify the factors that led to the simultaneous growth of FYI and melt of SYI. The different evolutions were explained by the difference in the conductive heat flux between the FYI and SYI. As the FYI was thin, there was a large temperature gradient between the ice base and the colder ice surface. This generated an upward conductive heat flux, which was larger than the heat flux from the ocean to the ice base, yielding basal growth of ice. In the case of the thicker SYI the temperature gradient and, hence, the conductive heat flux were smaller, and not sufficient to balance the oceanic heat flux at the ice base, yielding basal ablation. Penetration of solar radiation affected the conductive heat flux in both cases, and the model results were sensitive to the initial ice temperature profile and the uncertainty of the oceanic heat flux.

After the NEEM (Greenland) deep ice-core drilling was declared terminated with respect to developing stratigraphic climate reconstructions, efforts were turned toward collecting basal ice-sheet debris and, if possible, drilling into the bedrock itself. In 2010, several meters of banded debris-rich ice were obtained under normal ice-drilling operations with the NEEM version of the Hans Tausen (HT) drill, but further penetration was obstructed by a rock in the path of the drill head at 2537.36 m. During short campaigns in 2011 and 2012, attempts were made to penetrate further using various reinforced ice cutters mounted on the HT drill head, tailored to cut through rock. These had some success in penetrating coarse material, but produced severely damaged cutters. Additionally a 51 mm diameter diamond cutting tipped rock drill was adapted to fit the NEEM drill. With this device, several additional meters of core containing subglacial sediments, rocks and rock fragments were collected. With these tools 1.39m of additional material were obtained during the 2011 field season, and 7.1 m during 2012. Subglacial water refreezing into the newly formed borehole hindered further penetration, and the bedrock interface was not reached before final closure of the NEEM Camp.

Despite their importance as freshwater reservoirs for downstream river systems, few glaciers in central Chile have been comprehensively surveyed. This study presents ground-penetrating radar (GPR) and field-based observations for characterizing the englacial and basal conditions of Glaciar Olivares Alfa (33°110 S, 70°130 W), central Chilean Andes. Using a 50 MHz radar mounted onto a helicopter platform, data were collected covering large portions of the glacier accumulation and ablation zones. The radar data revealed boundaries of a temperate-ice layer at the base of the eastern body of Glaciar Olivares Alfa which appears to be covered by colder ice that extends throughout large parts of the glacier. The thickness of the temperate ice layer is highly variable across the glacier, being on average 40% of the total ice thickness. Radar data analyses reveal regions of cold ice at the bottom/base of the glacier and also patterns of highly saturated sediments beneath the glacier. Using GPR data, this study represents the most exhaustive analysis of glacier ice structure performed in the central Chilean Andes. The results will enable improved estimations of the glacier’s mass balance and ice dynamics, helping us to understand its further development and its impact on water availability.

Glacier dynamics are inextricably linked to the basal conditions of glaciers. Seismic reflection methods can image the glacier bed under certain conditions. However, where a seismically thin layer of material is present at the bed, traditional analyses may fail to fully characterize bed properties. We use a targeted full-waveform inversion algorithm to quantify the basal-layer parameters of a mountain glacier: thickness (d), P-wave velocity (α) and density (ρ). We simultaneously invert for the seismic quality factor (Q) of the bulk glacier ice. The inversion seeks to minimize the difference between the data and a one-dimensional reflectivity algorithm using a gradient-based search with starting values initialized from a Monte Carlo scheme. We test the inversion algorithm on four basal layer synthetic data traces with 5% added Gaussian noise. The inversion retrieved thin-layer parameters within 10% of synthetic test parameters with the exception of seismic Q. For the seismic dataset from Bench Glacier, Alaska, USA, inversion results indicate a thin basal layer of debris-rich ice within the study area having mean velocity 4000 ± 700 m s–1, density 1900 ± 200 kg m–3 and thickness 6 ± 1.5 m.

Large numbers of small valley glaciers on Svalbard were thicker and more extensive during the Little Ice Age (LIA), demonstrated by prominent ice-cored moraines up to several kilometres beyond present-day margins. The majority of these glaciers have since experienced a long period of strongly negative mass balance during the 20th century and are now largely frozen to their beds, indicating they are likely to have undergone a thermal transition from a polythermal to a cold-based regime. We present evidence for such a switch by reconstructing the former flow dynamics and thermal regime of Tellbreen, a small cold-based valley glacier in central Spitsbergen, based on its basal sequence and glaciological structures. Within the basal sequence, the underlying matrix-supported diamict is interpreted as saturated subglacial traction till which has frozen at the bed, indicating that the thermal switch has resulted in a cessation of subglacial sediment deformation due to freezing of the former deforming layer. This is overlain by debris-poor dispersed facies ice, interpreted to have formed through strain-induced metamorphism of englacial ice. The sequential development of structures includes arcuate fracture traces, interpreted as shear planes formed in a compressive/transpressive stress regime; and fracture traces, interpreted as healed extensional crevasses. The formation of these sediment/ice facies and structures is indicative of dynamic, warm-based flow, most likely during the LIA when the glacier was significantly thicker.

Semi-permanent snow is part of the continuum between seasonal snow and glacier ice. Although ubiquitous in the High Arctic, most late-lying snow banks and snow beds have lost their perennial status over the past decade as the summers have become progressively warmer. The loss over the past decade is the most unprecedented since aerial photography of the Canadian Arctic islands was first undertaken over half a century ago, and it has produced observable thermal, hydrological and ecological impacts. Where the ground has become exposed beneath the perennial snow cover, seasonal ground thaw has deepened. Tundra ponds and patchy wetlands fed principally by meltwater in the summer have suffered water-level decline or desiccation. The water balance of headwater basins is also affected, losing a buffering vehicle that accumulates storage surplus from the wet cool years to support streamflow and evaporation in the dry warm years. The tundra vegetation, already sparse, undergoes changes in the long term. As an essential source of water in the polar desert environment, the widespread distribution of semi-permanent snow magnifies its Arctic-wide importance.