We present a workflow to track icebergs in proglacial fjords using oblique time-lapse photos and the Lucas-Kanade optical flow algorithm. We employ the workflow at LeConte Bay, Alaska, where we ran five time-lapse cameras between April 2016 and September 2017, capturing more than 400 000 photos at frame rates of 0.5–4.0 min−1 . Hourly to daily average velocity fields in map coordinates illustrate dynamic currents in the bay, with dominant downfjord velocities (exceeding 0.5 m s−1 intermittently) and several eddies. Comparisons with simultaneous Acoustic Doppler Current Profiler (ADCP) measurements yield best agreement for the uppermost ADCP levels (∼ 12 m and above), in line with prevalent small icebergs that trace near-surface currents. Tracking results from multiple cameras compare favorably, although cameras with lower frame rates (0.5 min−1 ) tend to underestimate high flow speeds. Tests to determine requisite temporal and spatial image resolution confirm the importance of high image frame rates, while spatial resolution is of secondary importance. Application of our procedure to other fjords will be successful if iceberg concentrations are high enough and if the camera frame rates are sufficiently rapid (at least 1 min−1 for conditions similar to LeConte Bay).

Water that pressurizes the base of glaciers and ice sheets enhances glacier velocities and modulates glacial erosion. Predicting ice flow and erosion therefore requires knowledge of subglacial channel evolution, which remains observationally limited. Here we demonstrate that detailed analysis of seismic ground motion caused by subglacial water flow at Mendenhall Glacier (Alaska) allows for continuous measurement of daily to subseasonal changes in basal water pressure gradient, channel size, and sediment transport. We observe intermittent subglacial water pressure gradient changes during the melt season, at odds with common assumptions of slowly varying, low-pressure channels. These observations indicate that changes in channel size do not keep pace with changes in discharge. This behavior strongly affects glacier dynamics and subglacial channel erosion at Mendenhall Glacier, where episodic periods of high water pressure gradients enhance glacier surface velocity and channel sediment transport by up to 30% and 50%, respectively. We expect the application of this framework to future seismic observations acquired at glaciers worldwide to improve our understanding of subglacial processes.

Subglacial discharge influences glacier basal motion and erodes and redeposits sediment. At tidewater glacier termini, discharge drives submarine terminus melting, affects fjord circulation, and is a central component of proglacial marine ecosystems. However, our present inability to track subglacial discharge and its variability significantly hinders our understanding of these processes. Here we report observations of hourly to seasonal variations in 1.5–10 Hz seismic tremor that strongly correlate with subglacial discharge but not with basal motion, weather, or discrete icequakes. Our data demonstrate that vigorous discharge occurs from tidewater glaciers during summer, in spite of fast basal motion that could limit the formation of subglacial conduits, and then abates during winter. Furthermore, tremor observations and a melt model demonstrate that drainage efficiency of tidewater glaciers evolves seasonally. Glaciohydraulic tremor provides a means by which to quantify subglacial discharge variations and offers a promising window into otherwise obscured glacierized environments.

We investigate the dynamics of ice mélange by analyzing rapid motion recorded by a time-lapse camera and terrestrial radar during several calving events that occurred at Jakobshavn Isbræ, Greenland. During calving events (1) the kinetic energy of the ice mélange is 2 orders of magnitude smaller than the total energy released during the events, (2) a jamming front propagates through the ice mélange at a rate that is an order of magnitude faster than the motion of individual icebergs, (3) the ice mélange undergoes initial compaction followed by slow relaxation and extension, and (4) motion of the ice mélange gradually decays before coming to an abrupt halt. These observations indicate that the ice mélange experiences widespread jamming during calving events and is always close to being in a jammed state during periods of terminus quiescence. We therefore suspect that local jamming influences longer timescale ice mélange dynamics and stress transmission.

The loss of the floating ice tongue on Jakobshavn Isbræ, Greenland, in the early 2000s has been concurrent with a pattern of thinning, retreat and acceleration leading to enhanced contribution to global sea level. These changes on decadal timescales have been well documented. Here we identify how the glacier responds to forcings on shorter timescales, such as from variations in surface melt, the drainage of supraglacial lakes and seasonal fluctuations in terminus position. Ice motion and surface melt were monitored intermittently from 2006 to 2008. Dual-frequency GPS were deployed 20–50 km upstream of the terminus along the glacier center line. Gaps in surface melt measurements were filled using a temperature-index model of ablation driven by surface air temperatures recorded during the same time period. Our results corroborate the premise that the primary factors controlling speeds on Jakobshavn Isbræ are terminus position and geometry. We also observe that surface speeds demonstrate a complex relationship with meltwater input: on diurnal timescales, velocities closely match changes in water input; however, on seasonal timescales a longer, more intense melt season was observed to effectively reduce the overall ice flow of the glacier for the whole year.

Full-glacier-thickness icebergs are frequently observed to capsize as they calve into the ocean. As they capsize they may collide with the glaciers’ termini; previous studies have hypothesized that such collisions are the source of teleseismic ‘glacial earthquakes’. We use laboratory-scale experiments, force-balance modeling and theoretical arguments to show that (1) the contact forces during these collisions are strongly influenced by hydrodynamic forces and (2) the associated glacial earthquake magnitudes (expressed as twice-integrated force histories) are related to the energy released by the capsizing icebergs plus a hydrodynamic term that is composed of drag forces and hydrodynamic pressure. Our experiments and first-order modeling efforts suggest that, due to hydrodynamic forces, both contact force and glacial earthquake magnitudes may not be directly proportional to the energy released by the capsizing icebergs (as might be expected). Most importantly, however, our results highlight the need to better understand the hydrodynamics of iceberg capsize prior to being able to accurately interpret seismic signals generated by iceberg collisions.

We propose a general framework for iceberg-calving models that can be applied to any calving margin. The framework is based on mass continuity, the assumption that calving rate and terminus velocity are not independent and the simple idea that terminus thickness following a calving event is larger than terminus thickness at the event onset. The theoretical, near steady-state analysis used to support and analyze the framework indicates that calving rate is governed, to first order, by ice thickness, thickness gradient, strain rate, mass-balance rate and backwards melting of the terminus; the analysis furthermore provides a physical explanation for a previously derived empirical relationship for ice-shelf calving (Alley and others, 2008). In the calving framework the pre- and post-calving terminus thicknesses are given by two unknown but related functions. The functions can vary independently of changes in glacier flow and geometry, and can therefore account for variations in calving behavior due to external forcings and/or self-sustaining calving processes (positive feedbacks). Although the calving framework does not constitute a complete calving model, any thickness-based calving criterion can easily be incorporated into the framework. The framework should be viewed as a guide for future attempts to parameterize calving.

Observations of surface motion and ice deformation from 2002–03 were used to infer mean stress fields in a cross-section of Black Rapids Glacier, Alaska, USA, over seasonal timescales. Basal shear stresses in a well-defined zone north of the center line (orographic left) were approximately 7% and 16% lower in spring and summer, respectively, than in winter. Correspondingly higher stresses were found near the margins. These changes in the basal shear stress distribution were sufficiently large to cause mean surface velocities to be 1.2 and 1.5 times larger in spring and summer than in winter. These results were inferred with a simple inverse finite-element flow model that can successfully reproduce bulk surface velocities and tiltmeter data. Stress redistribution between the well-defined zone and the margins may also occur over much shorter time periods as a result of rapidly changing basal conditions (ice–bed decoupling or enhanced till deformation), thereby causing large variations in surface velocity and strongly influencing the glacier’s net motion during summer.

We assess how the propagation of high-frequency elastic-flexural waves through an ice shelf is modified by the presence of spatially periodic crevasses. Analysis of the normal modes supported by the ice shelf with and without crevasses reveals that a periodic crevasse distribution qualitatively changes the mechanical response. The normal modes of an ice shelf free of crevasses are evenly distributed as a function of frequency. In contrast, the normal modes of a crevasse-ridden ice shelf are distributed unevenly. There are ‘band gaps’, frequency ranges over which no eigenmodes exist. A model ice shelf that is 50 km in lateral extent and 300 m thick with crevasses spaced 500 m apart has a band gap from 0.2 to 0.38 Hz. This is a frequency range relevant for ocean-wave/ice-shelf interactions. When the outermost edge of the crevassed ice shelf is oscillated at a frequency within the band gap, the ice shelf responds very differently from a crevasse-free ice shelf. The flexural motion of the crevassed ice shelf is confined to a small region near the outermost edge of the ice shelf and effectively ‘blocked’ from reaching the interior.

We investigated seismic signals generated during a large-scale, multiple iceberg calving event that occurred at Jakobshavn Isbræ, Greenland, on 21 August 2009. The event was recorded by a high-rate time-lapse camera and five broadband seismic stations located within a few hundred kilometers of the terminus. During the event two full-glacier-thickness icebergs calved from the grounded (or nearly grounded) terminus and immediately capsized; the second iceberg to calve was two to three times smaller than the first. The individual calving and capsize events were well-correlated with the radiation of low-frequency seismic signals (<0.1 Hz) dominated by Love and Rayleigh waves. In agreement with regional records from previously published ‘glacial earthquakes’, these low-frequency seismic signals had maximum power and/or signal-to-noise ratios in the 0.05–0.1 Hz band. Similarly, full waveform inversions indicate that these signals were also generated by horizontal single forces acting at the glacier terminus. The signals therefore appear to be local manifestations of glacial earthquakes, although the magnitudes of the signals (twice-time integrated force histories) were considerably smaller than previously reported glacial earthquakes. We thus speculate that such earthquakes may be a common, if not pervasive, feature of all full-glacier-thickness calving events from grounded termini. Finally, a key result from our study is that waveform inversions performed on low-frequency, calving-generated seismic signals may have only limited ability to quantitatively estimate mass losses from calving. In particular, the choice of source time function has little impact on the inversion but dramatically changes the earthquake magnitude. Accordingly, in our analysis, it is unclear whether the smaller or larger of the two calving icebergs generated a larger seismic signal.