Timothy C. Bartholomaus

Distributed subglacial discharge drives significant submarine melt at a Greenland tidewater glacier

This is the publisher’s final pdf. The article is copyrighted by American Geophysical Union and published by John Wiley & Sons, Inc. It can be found at: http://agupubs.onlinelibrary.wiley.com/agu/journal/10.1002/%28ISSN%291944-8007/

The impact of glacier geometry on meltwater plume structure and submarine melt in Greenland fjords

Meltwater from the Greenland Ice Sheet often drains subglacially into fjords, driving upwelling plumes at glacier termini. Ocean models and observations of submarine termini suggest that plumes enhance melt and undercutting, leading to calving and potential glacier destabilization. Here we systematically evaluate how simulated plume structure and submarine melt during summer months depends on realistic ranges of subglacial discharge, glacier depth, and ocean stratification from 12 Greenland fjords. Our results show that grounding line depth is a strong control on plume-induced submarine melt: deep glaciers produce warm, salty subsurface plumes that undercut termini and shallow glaciers produce cold, fresh surface-trapped plumes that can overcut termini. Due to sustained upwelling velocities, plumes in cold, shallow fjords can induce equivalent depth-averaged melt rates compared to warm, deep fjords. These results detail a direct ocean-ice feedback that can affect the Greenland Ice Sheet.

Contrasts in the response of adjacent fjords and glaciers to ice-sheet surface melt in West Greenland

ABSTRACT Neighboring tidewater glaciers often exhibit asynchronous dynamic behavior, despite relatively uniform regional atmospheric and oceanic forcings. This variability may be controlled by a combination of local factors, including glacier and fjord geometry, fjord heat content and circulation, and glacier surface melt. In order to characterize and understand contrasts in adjacent tidewater glacier and fjord dynamics, we made coincident ice-ocean-atmosphere observations at high temporal resolution (minutes to weeks) within a 10 000 km2 area near Uummannaq, Greenland. Water column velocity, temperature and salinity measurements reveal systematic differences in neighboring fjords that imply contrasting circulation patterns. The observed ocean velocity and hydrography, combined with numerical modeling, suggest that subglacial discharge plays a major role in setting fjord conditions. In addition, satellite remote sensing of seasonal ice flow speed and terminus position reveal both speedup and slow-down in response to melt, as well as differences in calving style among the neighboring glaciers. Glacier force budgets and modeling also point toward subglacial discharge as a key factor in glacier behavior. For the studied region, individual glacier and fjord geometry modulate subglacial discharge, which leads to contrasts in both fjord and glacier dynamics.

Subglacial discharge at tidewater glaciers revealed by seismic tremor

Abstract 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.

Subseasonal changes observed in subglacial channel pressure, size, and sediment transport

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.

An Empirical Approach for Estimating Stress-Coupling Lengths for Marine-Terminating Glaciers

Variability in the dynamic behavior of marine-terminating glaciers is poorly understood, despite an increase in the abundance and resolution of observations. When paired with ice thicknesses, surface velocities can be used to quantify the dynamic redistribution of stresses in response to environmental perturbations through computation of the glacier force balance. However, because the force balance is not purely local, force balance calculations must be performed at the spatial scale over which stresses are transferred within glacier ice, or the stress-coupling length (SCL). Here we present a new empirical method to estimate the SCL for marine-terminating glaciers using high-resolution observations. We use the empirically-determined periodicity in resistive stress oscillations as a proxy for the SCL. Application of our empirical method to two well-studied tidewater glaciers (Helheim Glacier, SE Greenland, and Columbia Glacier, Alaska, USA) demonstrates that SCL estimates obtained using this approach are consistent with theory (i.e., can be parameterized as a function of the ice thickness) and with prior, independent SCL estimates. In order to accurately resolve stress variations, we suggest that similar empirical stress-coupling parameterizations be employed in future analyses of glacier dynamics.

Tidal and seasonal variations in calving flux observed with passive seismology

The seismic signatures of calving events, i.e., calving icequakes, offer an opportunity to examine calving variability with greater precision than is available with other methods. Here using observations from Yahtse Glacier, Alaska, we describe methods to detect, locate, and characterize calving icequakes. We combine these icequake records with a coincident, manually generated record of observed calving events to develop and validate a statistical model through which we can infer iceberg sizes from the properties of calving icequakes. We find that the icequake duration is the single most significant predictor of an icebergs size. We then apply this model to 18 months of seismic recordings and find elevated iceberg calving flux during the summer and fall and a pronounced lull in calving during midwinter. Calving flux is sensitive to semidiurnal tidal stage. Large calving events are tens of percent more likely during falling and low tides than during rising and high tides, consistent with a view that deeper water has a stabilizing influence on glacier termini. Multiple factors affect the occurrence of mechanical fractures that ultimately lead to iceberg calving. At Yahtse Glacier, seismology allows us to demonstrate that variations in the rate of submarine melt are a dominant control on iceberg calving rates at seasonal timescales. On hourly to daily timescales, tidal modulation of the normal stress against the glacier terminus reveals the nonlinear glacier response to changes in the near-terminus stress field.

The periodic topography of ice stream beds: Insights from the Fourier spectra of mega-scale glacial lineations

Ice stream bed topography contains key evidence for the ways ice streams interact with, and are potentially controlled by, their beds. Here we present the first application of two–dimensional Fourier analysis to 22 marine and terrestrial topographies from 5 regions in Antarctica and Canada, with and without mega-scale glacial lineations (MSGLs). We find that the topography of MSGL-rich ice stream sedimentary beds is characterized by multiple, periodic wavelengths between 300 and 1200 m and amplitudes from decimeters to a few meters. This periodic topography is consistent with the idea that instability is a key element to the formation of MSGL bedforms. Dominant wavelengths vary among locations and, on one paleo ice stream bed, increase along the direction of ice flow by 1.7±0.52% km-1. We suggest that these changes are likely to reflect pattern evolution via downstream wavelength coarsening, even under potentially steady ice stream geometry and flow conditions. The amplitude of MSGLs is smaller than that of other fluvial and glacial topographies, but within the same order of magnitude. However, MSGLs are a striking component of ice stream beds because the topographic amplitude of features not aligned with ice flow is reduced by an order of magnitude relative to those oriented with the flow direction. This study represents the first attempt to automatically derive the spectral signatures of MSGLs. It highlights the plausibility of identifying these landform assemblages using automated techniques and provides a benchmark for numerical models of ice stream flow and subglacial landscape evolution.

Near-glacier surveying of a subglacial discharge plume: implications for plume parameterizations†

At tidewater glaciers, plume dynamics affect submarine melting, fjord circulation, and the mixing of meltwater. Models often rely on buoyant plume theory to parameterize plumes and submarine melting; however, these parameterizations are largely untested due to a dearth of near-glacier measurements. Here, we present a high-resolution ocean survey by ship and remotely-operated boat near the terminus of Kangerlussuup Sermia in west Greenland. These novel observations reveal the 3D structure and transport of a near-surface plume, originating at a large undercut conduit in the glacier terminus, that is inconsistent with axisymmetric plume theory, the most common representation of plumes in ocean-glacier models. Instead, the observations suggest a wider upwelling plume – a ‘truncated’ line plume of ∼200 m width – with higher entrainment and plume-driven melt compared to the typical axisymmetric representation. Our results highlight the importance of a subglacial outlets geometry in controlling plume dynamics, with implications for parameterizing the exchange flow and submarine melt in glacial fjord models.

Inversion of a glacier hydrology model

ABSTRACT The subglacial hydrologic system exerts strong controls on the dynamics of the overlying ice, yet the parameters that govern the evolution of this system are not widely known or observable. To gain a better understanding of these parameters, we invert a spatially averaged model of subglacial hydrology from observations of ice surface velocity and outlet stream discharge at Kennicott Glacier, Wrangell Mountains, AK, USA. To identify independent parameters, we formally non-dimensionalize the forward model. After specifying suitable prior distributions, we use a Markov-chain Monte Carlo algorithm to sample from the distribution of parameter values conditioned on the available data. This procedure gives us not only the most probable parameter values, but also a rigorous estimate of their covariance structure. We find that the opening of cavities due to sliding over basal topography and turbulent melting are of a similar magnitude during periods of large input flux, though turbulent melting also exhibits the greatest uncertainty. We also find that both the storage of water in the englacial system and the exchange of water between englacial and subglacial systems are necessary in order to explain both surface velocity observations and the relative attenuation in the amplitude of diurnal signals between input and output flux observations.