Understanding the dynamic response of glaciers to climate change is vital for assessing water resources and hazards, and subglacial hydrology is a key player in glacier systems. Traditional observations of subglacial hydrology are spatially and temporally limited, but recent seismic deployments on and around glaciers show the potential for comprehensive observation of glacial hydrologic systems. We present results from a high-density seismic deployment spanning the surface of Lemon Creek Glacier, Alaska. Our study coincided with a marginal lake drainage event, which served as a natural experiment for seismic detection of changes in subglacial hydrology. We observed glaciohydraulic tremor across the surface of the glacier that was generated by the subglacial hydrologic system. During the lake drainage, the relative changes in seismic tremor power and water flux are consistent with pressurization of the subglacial system of only the upper part of the glacier. This event was not accompanied by a significant increase in glacier velocity; either some threshold necessary for rapid basal motion was not attained, or, plausibly, the geometry of Lemon Creek Glacier inhibited speedup. This pressurization event would have likely gone undetected without seismic observations, demonstrating the power of cryoseismology in testing assumptions about and mapping the spatial extent of subglacial pressurization.

The stability of the world’s largest glaciers and ice sheets depends on mechanical and thermodynamic processes occurring at the glacier–ocean boundary. A buoyant agglomeration of icebergs and sea ice, referred to as ice mélange, often forms along this boundary and has been postulated to affect ice-sheet mass losses by inhibiting iceberg calving. Here, we use terrestrial radar data sampled every 3 min to show that calving events at Jakobshavn Isbræ, Greenland, are preceded by a loss of flow coherence in the proglacial ice mélange by up to an hour, wherein individual icebergs flowing in unison undergo random displacements. A particle dynamics model indicates that these fluctuations are likely due to buckling and rearrangements of the quasi-two-dimensional material. Our results directly implicate ice mélange as a mechanical inhibitor of iceberg calving and further demonstrate the potential for real-time detection of failure in other geophysical granular materials.

Lemon Creek Glacier, a temperate valley glacier in the Juneau Icefield of Southeast Alaska, is the site of long running (>60 years) glaciological studies. However, the most recent published estimates of its thickness and subglacial topography come from two ∼50 years old sources that are not in agreement and do not account for the effects of years of negative mass balance. We collected a 1-km long active-source seismic line on the upper section of the glacier parallel and near to the centerline of the glacier, roughly straddling the equilibrium-line altitude. We used these data to perform joint reflection-refraction velocity modeling and reflection imaging of the glacier bed. We find that this upper section of Lemon Creek Glacier is as much as 150 m (∼65%) thicker than previously suggested with a large overdeepening in an area previously believed to have a uniform thickness. Our results lead us to reinterpret the impact of basal motion on ice flow and have a significant impact on expectations of subglacial hydrology. We suggest that further efforts to develop a whole-glacier model of subglacial topography are necessary to support studies that require accurate models of ice thickness and subglacial topography.

Sedimentary processes are known to help facilitate tidewater glacier advance, but their role in modulating retreat is uncertain and poorly quantified. In this study we use repeated seafloor bathymetric surveys and satellite‐derived terminus positions from LeConte Glacier, Alaska, to evaluate the evolution of a morainal bank and related changes in terminus dynamics over a 17‐year period. The glacier experienced a rapid retreat between 1994 and 1999, before stabilizing at a constriction in the fjord. Since then, the glacier terminus has remained stabilized while constructing a morainal bank up to 140 m high in water depths of 240–260 m, with rates of sediment delivery of 3.3 Å~ 105 to 3.8 Å~ 105 m3 a−1. Based on repeated interannual surveys between 2016 and 2018, the moraine is a dynamic feature characterized by push ridges, evidence of active gravity flows, and bulldozing by the glacier at rates of up to meters per day. Beginning in 2016, the summertime terminus has become increasingly retracted, revealing a newly emerging basin potentially signaling the onset of renewed retreat. Between 2000 and 2016, the growing moraine reduced the exposed submarine area of the terminus by up to 22%, altered the geometry of the terminus during seasonal advances, and altered the terminus stress balance. These feedbacks for calving, melting, and ice flow likely represent mechanisms whereby moraine growth may delay glacier retreat, in a system where readvance is unlikely.

Suicide Basin is a partly glacierized marginal basin of Mendenhall Glacier, Alaska, that has released glacier lake outburst floods (GLOFs) annually since 2011. The floods cause inundation and erosion in the Mendenhall Valley, impacting homes and other infrastructure. Here, we utilize in-situ and remote sensing data to assess the recent evolution and current state of Suicide Basin. We focus on the 2018 and 2019 melt seasons, during which we collected most of our data, partly using unmanned aerial vehicles (UAVs). To provide longer-term context, we analyze DEMs collected since 2006 and model glacier surface mass balance over the 2006–2019 period. During the 2018 and 2019 outburst flood events, Suicide Basin released ∼ 30 Å~ 106 m3 of water within approximately 4–5 days. Since lake drainage was partial in both years, these ∼ 30 Å~ 106 m3 represent only a fraction (∼ 60%) of the basin’s total storage capacity. In contrast to previous years, subglacial drainage was preceded by supraglacial outflow over the ice dam, which lasted ∼ 1 day in 2018 and 6 days in 2019. Two large calving events occurred in 2018 and 2019, with submerged ice breaking off the main glacier during lake filling, thereby increasing the basin’s storage capacity. In 2018, the floating ice in the basin was 36 m thick on average. In 2019, ice thickness was 29 m, suggesting rapid decay of the ice tongue despite increasing ice inflow from Mendenhall Glacier. The ice dam at the basin entrance thinned by more than 5 m a–1 from 2018 to 2019, which is approximately double the rate of the reference period 2006–2018. While ice-dam thinning reduces water storage capacity in the basin, that capacity is increased by declining ice volume in the basin and longitudinal lake expansion, with the latter process challenging to predict. The potential for premature drainage onset (i.e., drainage before the lake’s storage capacity is reached), intermittent drainage decelerations, and early drainage termination further complicates prediction of future GLOF events.

In most models of glacial erosion, glacier sliding velocity is hypothesized to control rates of bedrock erosion. If this hypothesis is correct, then the elevation difference between hanging and trunk valley floors, the hang height, should be dictated by the relative sliding velocities of the glaciers that occupied these valleys. By considering sliding velocity to be proportional to balance velocity and using mass continuity, hang height is expressed in terms of glacier catchment areas, slopes, and widths, which can be estimated for past glaciers from the morphology of glacial valleys. These parameters were estimated for 46 hanging valleys and their trunk valleys in three adjacent regions of Jasper National Park. The variability in valley morphology can account for 55–85% of the hang height variability if erosion rate scales with balance velocity raised to a power of 1/3. This correspondence is in spite of spatial variations in glaciation duration, snow accumulation rates, and other variables that likely affected hang heights but cannot be readily estimated and so are not included in our formulation. Thus it appears that balance velocity, and by extension, sliding velocity if the two are proportional, may be a reasonable control variable for assessing erosion rate.

Jakobshavn Isbræ, a tidewater glacier that produces some of Greenland’s largest icebergs and highest speeds, reached record-high flow rates in 2012 (Joughin and others, 2014). We use terrestrial radar interferometric observations from August 2012 to characterize the events that led to record-high flow. We find that the highest speeds occurred in response to a small calving retreat, while several larger calving events produced negligible changes in glacier speed. This non-linear response to calving events suggests the terminus was close to flotation and therefore highly sensitive to terminus position. Our observations indicate that a glacier’s response to calving is a consequence of two competing feedbacks: (1) an increase in strain rates that leads to dynamic thinning and faster flow, thereby promoting desta- bilization, and (2) an increase in flow rates that advects thick ice toward the terminus and promotes restabilization. The competition between these feedbacks depends on temporal and spatial variations in the glacier’s proximity to flotation. This study highlights the importance of dynamic thinning and advective processes on tidewater glacier stability, and further suggests the latter may be limiting the current retreat due to the thick ice that occupies Jakobshavn Isbræ’s retrograde bed.

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.