Detection, attribution and projection of mass loss from the Greenland Ice Sheet has been a central focus of the glaciological community, with surface meltwater thought to play a key role in feedbacks ...that could accelerate sea-level rise. While the prospect of runaway sliding has faded, much remains uncertain when it comes to the role of surface runoff and subglacial discharge in Greenland's future.
Modelling water flow under glaciers and ice sheets Flowers, Gwenn E.
Proceedings of the Royal Society. A, Mathematical, physical, and engineering sciences,
04/2015, Letnik:
471, Številka:
2176
Journal Article
Recenzirano
Odprti dostop
Recent observations of dynamic water systems beneath the Greenland and Antarctic ice sheets have sparked renewed interest in modelling subglacial drainage. The foundations of today's models were laid ...decades ago, inspired by measurements from mountain glaciers, discovery of the modern ice streams and the study of landscapes evacuated by former ice sheets. Models have progressed from strict adherence to the principles of groundwater flow, to the incorporation of flow 'elements' specific to the subglacial environment, to sophisticated two-dimensional representations of interacting distributed and channelized drainage. Although presently in a state of rapid development, subglacial drainage models, when coupled to models of ice flow, are now able to reproduce many of the canonical phenomena that characterize this coupled system. Model calibration remains generally out of reach, whereas widespread application of these models to large problems and real geometries awaits the next level of development.
Long‐term records of the flow patterns and dynamics of surge‐type glaciers improve our understanding of their underlying dynamic processes, and are critical to better resolve their contribution to a ...changing cryosphere. We adapt a modeling approach designed to emulate glacier surging and fold kinematics using the full Stokes ice‐flow model Elmer/Ice to simulate surging of the Dusty Glacier, located in the St. Elias Mountains, Canada. We combine distributed mass‐balance and numerical ice‐flow models to reconstruct the fold kinematics of the 2001–2003 surge of the Dusty Glacier by comparing model results to Landsat‐7 and Sentinel‐2 imagery, and assess the sensitivity of centennial‐scale modeled glacier structure to different mass balance and sliding parameterizations. This study demonstrates the feasibility of using the approach to reconstruct the surface structure kinematics of a surge‐type glacier in nature, highlighting its potential application to other surge‐type glaciers and regions.
Plain Language Summary
Glaciers can exhibit irregular flow patterns that complicate predictions of their evolution over the coming decades and centuries. We present a method for reproducing iconic surface structures known as folded medial moraines using glaciological modeling. These moraines are wavy flow patterns found on surge‐type glaciers, highlighted by sediment deposited onto the ice that traces their path. By reconstructing these patterns, the underlying climate and sliding conditions that contributed to the glacier's past flow can be identified. Improving our knowledge of these conditions can help improve glacier flow models. We demonstrate that our methodology successfully reconstructs the flow patterns present on a large surge‐type glacier in Yukon, Canada, and explore its past flow history, and possible future, based on these results.
Key Points
We use a distributed mass‐balance model and Elmer/Ice to reconstruct the 2001–2003 surge kinematics of the Dusty Glacier, Yukon, Canada
We explore the centennial‐scale sensitivity of glacier surface fold geometry to mass balance and sliding parameterizations
This study is a proof‐of‐concept for further model reconstructions of the past dynamics of surge‐type glaciers
Five decades of radioglaciology Schroeder, Dustin M.; Bingham, Robert G.; Blankenship, Donald D. ...
Annals of glaciology,
04/2020, Letnik:
61, Številka:
81
Journal Article
Recenzirano
Odprti dostop
Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, ...airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects.
Abstract One fifth of Earth's volcanoes are covered by snow or ice and many have active geothermal systems that interact with the overlying ice. These glaciovolcanic interactions can melt voids into ...glaciers, and are subject to controls exerted by ice dynamics and geothermal heat output. Glaciovolcanic voids have been observed to form prior to volcanic eruptions, which raised concerns when such features were discovered within Job Glacier on Qw̓elqw̓elústen (Mount Meager Volcanic Complex), British Columbia, Canada. In this study we model the formation, evolution, and steady-state morphology of glaciovolcanic voids using analytical and numerical models. Analytical steady-state void geometries show cave height limited to one quarter of the ice thickness, while numerical model results suggest the void height h scales with ice thickness H and geothermal heat flux $\dot {Q}$ as $h/H = a H^b \dot {Q}^c$ , with exponents b = − n /2 and c = 1/2 where n is the creep exponent. Applying this scaling to the glaciovolcanic voids within Job Glacier suggests the potential for total geothermal heat flux in excess of 10 MW. Our results show that relative changes in ice thickness are more influential in glaciovolcanic void formation and evolution than relative changes in geothermal heat flux.
Abstract
We investigate unusual discontinuous glacier motion on Thompson Glacier, Umingmat Nunaat, Arctic Canada, using synthetic aperture radar (SAR) images and ice-flow modeling. A novel ...intensity-rescaling scheme is developed to reduce errors in high-resolution speckle tracking, resulting in a ~25% improvement in accuracy. Interferometric SAR (InSAR) and speckle tracking using high resolution RADARSAT-2 data indicate velocity discontinuities of up to 1 cm d
−1
across deep and longitudinally extensive supraglacial channels on Thompson Glacier. We use a cross-sectional finite-element ice-flow model to determine the conditions under which velocity discontinuities of the observed magnitude and signature are possible. The modeling suggests that discontinuous motion across (long and straight) supraglacial channels can occur without ice fracture and under a wide variety of glacier thermal structures, including in fully temperate glaciers. Despite the wide range of conditions conducive to discontinuous motion, the form we observe requires that the associated channels be deep, longitudinally extensive and located in regions of lateral shearing. We speculate that these combined conditions are rare except on polythermal glaciers, where drainage features such as moulins are comparatively scarce and lower deformation rates allow channels to incise consistently and persist over many years.
The Kaskawulsh Glacier is an iconic outlet draining the icefields of the St. Elias Mountains in Yukon, Canada. We determine and attempt to interpret its catchment-wide mass budget since 2007. Using ...SPOT5/6/7 data we estimate a 2007–18 geodetic balance of −0.46 ± 0.17 m w.e. a−1. We then compute balance fluxes and observed ice fluxes at nine flux gates to examine the discrepancy between the climatic mass balance and internal mass redistribution by glacier flow. Balance fluxes are computed using a fully distributed mass-balance model driven by downscaled and bias-corrected climate-reanalysis data. Observed fluxes are calculated using NASA ITS_LIVE surface velocities and glacier cross-sectional areas derived from ice-penetrating radar data. We find the glacier is still in the early stages of dynamic adjustment to its mass imbalance. We estimate a committed terminus retreat of ~23 km under the 2007–18 climate and a lower bound of 46 km3 of committed ice loss, equivalent to ~15% of the total glacier volume.
Field data and numerical modeling show that glaciations have the potential either to enhance relief or to dampen topography. We aim to model the effect of the subglacial hydraulic system on ...spatiotemporal patterns of glacial erosion by abrasion and quarrying on time scales commensurate with drainage system fluctuations (e.g., seasonal to annual). We use a numerical model that incorporates a dual-morphology subglacial drainage system coupled to a higher-order ice-flow model and process-specific erosion laws. The subglacial drainage system allows for a dynamic transition between two morphologies: the distributed system, characterized by an increase in basal water pressure with discharge, and the channelized system, which exhibits a decrease in equilibrium water pressure with increasing discharge. We apply the model to a simple synthetic glacier geometry, drive it with prescribed meltwater input variations, and compute sliding and erosion rates over a seasonal cycle. When both distributed and channelized systems are included, abrasion and sliding maxima migrate ~20% up-glacier compared to simulations with distributed drainage only. Power-law sliding generally yields to a broader response of abrasion to water pressure changes along the flowline compared to Coulomb-friction sliding. Multi-day variations in meltwater input elicit a stronger abrasion response than either diurnal- or seasonal variations alone for the same total input volume. An increase in water input volume leads to increased abrasion. We find that ice thickness commensurate with ice sheet outlet glaciers can hinder the up-glacier migration of abrasion. Quarrying patterns computed with a recently published law differ markedly from calculated abrasion patterns, with effective pressure being a stronger determinant than sliding speeds of quarrying rates. These variations in calculated patterns of instantaneous erosion as a function of hydrology-, sliding-, and erosion-model formulation, as well as model forcing, may lead to significant differences in predicted topographic profiles on long time scales.
•Hydrologically-coupled ice-flow model was used to calculate seasonal erosion patterns.•Process-specific erosion laws were implemented for abrasion and quarrying.•Channelized drainage leads to up-glacier migration of abrasion rates.•Multi-day variations in meltwater input are more important than diurnal or seasonal.•Quarrying patterns with Iverson's (2012) law differ from abrasion patterns.