Seasonal lake Gornersee forms at the confluence of Gornergletscher and Grenzgletscher, Switzerland, and experiences outburst floods annually in midsummer. To study the interplay between lake drainage, glacier movement and crevasse activity, high-frequency seismometers and GPS receivers were deployed in networks near Gornersee during the summer ablation seasons of 2004, 2006 and 2007. We use a Rayleigh wave coherence method to locate 3289, 7939 and 4087 icequakes, respectively, primarily along well-defined surface crevasses. We calculate two-dimensional strains from triads of GPS stations and find mean differential strain rates of ~300 × 10−6 d−1 with diurnal variations up to 800 × 10−6 d−1. Crevasse icequake activity and glacial velocity are highest during early season, then decrease as meltwater channels erode and subglacial water pressure decreases. Glacial response to Gornersee drainage varied year-to-year, with icequake activity promoted at some crevasses and inhibited at others, suggesting syn-drainage icequakes may be indicative of local drainage patterns and small-scale features of the stress field. Diurnal pulses in icequake activity exhibit peak activity at different times of day in different locations, coincident with a southeast-to-northwest trending concentrated shear zone near the Gornergletscher–Grenzgletscher confluence, likely due to differences in the timing of peak strain rate in these regions.

Glacier-erosion rates range across orders of magnitude, and much of this variation cannot be attributed to basal sliding rates. Subglacial till acts as lubricating ‘fault gouge’ or ‘sawdust’, and must be removed for rapid subglacial bedrock erosion. Such erosion occurs especially where and when moulin-fed streams access the bed and are unconstrained by supercooling or other processes. Streams also may directly erode bedrock, likely with strong time-evolution. Erosion is primarily by quarrying, aided by strong fluctuations in the water system driven by variable surface melt and by subglacial earthquakes. Debris-bed friction significantly affects abrasion, quarrying and general glacier flow. Frost heave drives cirque headwall erosion as winter cold air enters bergschrunds, creating temperature gradients to drive water flow along premelted films to growing ice lenses that fracture rock, and the glacier removes the resulting blocks. Recent subglacial bedrock erosion and sediment flux are in many cases much higher than long-term averages. Over glacial cycles, evolution of glacial-valley form feeds back strongly on erosion and deposition. Most of this is poorly quantified, with parts open to argument. Glacial erosion and interactions are important to tectonic and volcanic processes as well as climate and biogeochemical fluxes, motivating vigorous research.

In January 2016, static GPS measurements were carried out in a 30 × 30 km2 area centered around Kunlun station at Dome Argus (Dome A), East Antarctica, to acquire high-precision 3-D geodetic coordinates at 49 sites. By comparing the coordinates with previous GPS measurements in 2008 and 2013 at the same sites, we constructed a detailed and long-term record of the ice-surface velocity field, 2008–2016, around Dome A. During this time span, the estimated ice-surface velocity ranges from 0.8 ± 0.3 to 28.7 ± 1.6 cm a−1, with a mean of 10.4 ± 0.3 cm a−1. From 2013 to 2016, the surface elevation of most Dome A areas exhibits a rising trend, and the maximum increase of snow surface elevation is 84.8 cm. The mean snow surface elevation change rate at Dome A is estimated to be 6.6 ± 0.7 cm a−1. The difference of 1.0 cm a−1 between the snow surface change rate derived from GPS and pole-height change rate from surface mass balance is suspected to be a result of a combination of firn densification and basal melt under Dome A.

Changing rates of water input can affect both the flow of glaciers and ice sheets and their propensity to crevasse. Here we examine geodetic and seismic observations during two substantial (10–18-times background velocity) rain-induced glacier accelerations at Haupapa/Tasman Glacier, New Zealand. Changes in rain rate result in glacier acceleration and associated uplift, which propagate down-glacier. This pattern of acceleration results in a change to the strain rate field, which correlates with an order of magnitude increase in the apparent seismicity rate and an overall down-glacier migration in located seismicity. After each acceleration event the apparent seismicity rate decreases to below the pre-acceleration rate for 3 days. This suggests that seismic events associated with surface crevasse growth occur early during phases of glacier acceleration due to elevated extensional stresses, and then do not occur again until stresses recover.

Measurements of a sample from ~580 m depth in the WAIS Divide (WDC06A) ice core reveal that bubbles are preferentially elongated in the basal plane of their parent grain, as expected if bubble shape preserves the record of dominant basal glide. This suggests that a method using bubbles as strain gauges could provide insights to grain-scale ice deformation. We introduce a technique using fabric and image analyses of paired thin and thick sections. Comparison of the crystallographic orientations of 148 grains and the shape orientations of 2377 intragrain bubbles reveals a strongly preferred elongation of bubbles in the grain basal planes (R2 = 0.96). Elongation magnitudes are consistent with a balance between ice flow deformation and diffusive restoration, with larger bubbles more elongated. Assuming bubbles record ice strain, grains with greater resolved stress on their basal planes from the far-field ice flow stresses show greater deformation, but with large variability suggesting that heterogeneity of the local stress field causes deformation even in unfavorably oriented grains. A correlation is also observed among bubble elongation, grain size, and bubble size, explaining a small but significant fraction of the variance ( P< 0.05), with implications for controls on ice deformation, as discussed here.

Late Holocene sediment deposits in Pine Island Bay, West Antarctica, are hypothesized to be linked to intensive meltwater drainage during the retreat of the paleo-Pine Island Ice Stream after the Last Glacial Maximum. The uppermost sediment units show an abrupt transition from ice-proximal debris to a draped silt during the late Holocene, which is interpreted to coincide with rapid deglaciation. The small scale and fine sorting of the upper unit could be attributed to origins in subglacial meltwater; however the thickness and deposition rate for this unit imply punctuated- rather than continuous-deposition. This, combined with the deposit's location seaward of large, bedrock basins, has led to the interpretation of this unit as the result of subglacial lake outbursts in these basins. However, the fine-scale sorting of the silt unit is problematic for this energetic interpretation, which should mobilize and deposit a wider range of sediment sizes. To resolve this discrepancy, we present an alternative mechanism in which the silt was sorted by a distributed subglacial water system, stored in bedrock basins far inland of the grounding line, and subsequently eroded at higher flow speeds during retreat. We demonstrate that this mechanism is physically plausible given the subglacial conditions during the late Holocene. We hypothesize that similar silt units observed elsewhere in Antarctica downstream of bedrock basins could be the result of the same mechanism.

Glacial seismicity provides important insights into glacier dynamic processes. We study the temporal distribution of cryogenic seismic signals (icequakes) at Holtedahlfonna, Svalbard, between April and August 2016 using a single three-component sensor. We investigate sources of observed icequakes using polarization analysis and waveform modeling. Processes responsible for five icequake categories are suggested, incorporating observations of previous studies into our interpretation. We infer that the most dominant icequake type is generated by surface crevasse opening through hydrofracturing. Secondly, bursts of high-frequency signals are presumably caused by repeated near-surface crevassing due to high strain rates during glacier fast-flow episodes. Furthermore, signals related to resonance in water-filled cracks, fracturing or settling events in dry firn or snow before the melt season, and processes at the glacier bed are observed. Amplitude of seismic background noise is clearly related to glacier runoff. We process ambient seismic noise to invert horizontal-to-vertical spectral ratios for a sub-surface seismic velocity model used to model icequake signals. Our study shows that a single seismic sensor provides useful information about seasonal ice dynamics in case deployment of a network is not feasible.

We present a highly detailed study of calving dynamics at Tunabreen, a tidewater glacier in Svalbard. A time-lapse camera was trained on the terminus and programmed to capture images every 3 seconds over a 28-hour period in August 2015, producing a highly detailed record of 34 117 images from which 358 individual calving events were distinguished. Calving activity is characterised by frequent events (12.8 events h−1) that are small relative to the spectrum of calving events observed, demonstrating the prevalence of small-scale calving mechanisms. Five calving styles were observed, with a high proportion of calving events (82%) originating at, or above, the waterline. The tidal cycle plays a key role in the timing of calving events, with 68% occurring on the falling limb of the tide. Calving activity is concentrated where meltwater plumes surface at the glacier front, and a ~ 5 m undercut at the base of the glacier suggests that meltwater plumes encourage melt-under-cutting. We conclude that frontal ablation at Tunabreen may be paced by submarine melt rates, as suggested from similar observations at glaciers in Svalbard and Alaska. Using submarine melt rate to calculate frontal ablation would greatly simplify estimations of tidewater glacier losses in prognostic models.

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.

Ice shelves regulate the rate of ice-sheet discharge along much of the Antarctic coastline. Pinning points, sites of localised grounding within floating ice, can in turn, regulate the flow and thickness of an ice shelf. While the net resistive effect of ice shelves has been quantified in a systematic way, few extant pinning points have been examined in detail. Here, complete force budgets are calculated and examined for ice rises and rumples in the Ross Ice Shelf, West Antarctica. The diverse features have different effects on ice shelf mechanics that do not depend simply on their size but may, we conclude, depend on the properties of seafloor materials.

Theory and experiments indicate that ice–bed separation during glacier slip over 2-D hard beds causes basal shear stress to reach a maximum at a particular slip velocity and decrease at higher velocities. We use the sliding theory of Lliboutry (1968) to explore how friction between debris particles in sliding ice and a rock bed affects this relationship between shear stress and slip velocity. Particle–bed contact forces and associated debris friction increase with increasing slip velocity, owing to increased rates of ice convergence with up-glacier facing surfaces. However, debris friction on diminished areas of the bed counteracts this effect as cavities grow. Thus, friction from debris alone increases only slightly with slip velocity, and for sediment particles larger than ~60 mm in diameter, debris friction peaks and decreases with increasing slip velocity. The effect on the sliding relationship is to steepen its rising limb and shift its shear stress peak to a slightly higher velocity. These results, which exclude the effect of debris friction on cavity size and debris concentrations above ~15%, indicate that the effect of debris in ice is to increase basal shear stress but not significantly change the form of the sliding relationship.

GPS measurements of tidal modulation of ice flow and seismicity within the grounding zone of Beardmore Glacier show that tidally induced fluctuations of horizontal flow are largest near the grounding line and decrease downstream. Seismic activity is continuous, but peaks occur on falling and rising tides. Beamforming methods reveal that most seismic events originate from two distinct locations, one on the grid-north side of the grounding zone, and one on the grid-south side. The broad pattern of deformation generated as Beardmore Glacier merges with the Ross Ice Shelf results in net extension along the grid-north side of the grounding zone and net compression along the grid-south side. During falling tides, seismic activity peaks on both sides because of increased vertical flexure across the grounding line. During rising tides, seismic activity in the region of extension on the grid-north side is relatively low because the tidal influence on both horizontal strain rate and vertical flexure is small. On the grid-south side during rising tides, however, tidally induced horizontal strain rates promote increased seismicity in regions of long-term compressional flow paths. Our study highlights how concurrent geodetic and seismic measurements provide insight into grounding-zone mechanics and their influence on ice-shelf buttressing.

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.

Over a year of seismic observations, ~5000 short duration icequakes were detected by a permanent broadband station installed at the Princess Elisabeth base, located ~180 km inland in eastern Dronning Maud Land, East-Antarctica. Icequake detection via seismic waveform pattern recognition indicates the presence of two dominating clusters of events, totalizing ~1500 icequakes. The corresponding icequake locations point towards two distinct zones of outcropping blue ice areas (BIAs) located respectively at 4 and 1 km from the seismic station, both on the leeward side of a nunatak protruding through the ice sheet. The temporal occurrence of these icequakes suggests a close genetic link with thermal contraction of ice caused by significant surface cooling controlled, in summer by variations in diurnal solar radiation and in winter by strong cooling during cold katabatic regimes. Further analysis demonstrates the dependence of these icequakes on the absolute surface temperature and on its temporal change. Besides providing information on the ice fracture mechanics and rheology, investigations of thermal icequakes may be regarded as a ground-based proxy for the monitoring of the thermal state of BIAs, and characterization of ice-sheet ablation zones.

We recorded the ice motion and icequakes on the floating part of Langhovde Glacier in East Antarctica to better understand the dynamic behavior of ice shelves and floating tongues. Diurnal and semi-diurnal variations in ice motion and seismicity were simultaneously observed at all four global navigation satellite system and three seismic stations over 2 weeks. The short-term along-flow ice motion is explained by the elastic response of the glacier to ocean tide-induced hydrostatic stress variations, which decayed at a rate of 0.8 km−1 toward the grounding line. We observed a large number of icequakes during mid-rising and high tides that covered a broad frequency range and formed two major groups of events centered at 10 and 120 Hz, respectively. The hourly occurrence rates were ~500 events h−1, with the observed seismicity consistent with fracture due to floating tongue bending. We also observed minor secondary peaks at high ice speeds, which could reflect surface cracking due to stretching or basal friction. Our observation demonstrates that tidal-modulation was the main factor to fracture the floating tongue of Langhovde Glacier.

Glacial erosion produces vast quantities of fine-grained sediment that has a far-reaching impact on Earth surface processes. To gain a better understanding of the production of glacial silt and clay, we use automated mineralogy to quantify the microstructure and mineralogy of rock and sediment samples from 20 basins in the St. Elias Mountains, Yukon, Canada. Sediments were collected from proglacial streams, while rock samples were collected from ice marginal outcrops and fragmented using electrical pulse disaggregation. For both rock fragments and sediments, we observe a log-normal distribution of grain sizes and a sub-micrometer terminal grain size. We find that the abrasion of silt and clay results in both rounding and the exploitation of through-going fractures. The abundance of inter- versus intragranular fractures depends on mineralogy and size. Unlike the relatively larger grains, where crushing and abrasion are thought to exploit and produce discrete populations of grain sizes, the comminution of fines leads to a grain size, composition and rounding that is continuously distributed across size, and highly dependent on source-rock properties.

Forced ice-shelf vibration modeling is performed using a full 3-D finite-difference elastic model, which also takes into account sub-ice seawater flow. The sea water flow is described by the wave equation. Ice-shelf flexure therefore results from hydrostatic pressure perturbations in the sub-ice seawater layer. Numerical experiments were undertaken for idealized rectangular ice-shelf geometry. The ice-plate vibrations were modeled for harmonic incoming pressure perturbations and for a wide range of incoming wave frequencies. The spectra showed distinct resonant peaks, which demonstrate the ability of the model to simulate a resonant-like motion. The spectra obtained by the full 3-D model are compared with exact solutions for the elastic thin plate with two fixed edges and two free edges. The spectra are also compared with the spectra modeled by the thin-plate Holdsworth and Glynn model (1978).

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.