Increasing numerical stability of mountain valley glacier simulations: implementation and testing of free-surface stabilization in Elmer/Ice

This paper concerns a numerical stabilization method for free-surface ice flow called the free-surface stabilization algorithm (FSSA). In the current study, the FSSA is implemented into the numerical ice-flow software Elmer/Ice and tested on synthetic two-dimensional (2D) glaciers, as well as on the real-world glacier of Midtre Lovénbreen, Svalbard. For the synthetic 2D cases it is found that the FSSA method increases the largest stable time-step size at least by a factor of 5 for the case of a gently sloping ice surface (∼ 3°) and by at least a factor of 2 for cases of moderately to steeply inclined surfaces (∼ 6° to 12°) on a fine mesh. Compared with other means of stabilization, the FSSA is the only one in this study that increases largest stable time-step sizes when used alone. Furthermore, the FSSA method increases the overall accuracy for all surface slopes. The largest stable time-step size is found to be smallest for the case of a low sloping surface, despite having overall smaller velocities. For an Arctic-type glacier, Midtre Lovénbreen, the FSSA method doubles the largest stable time-step size; however, the accuracy is in this case slightly lowered in the deeper parts of the glacier, while it increases near edges. The implication is that the non-FSSA method might be more accurate at predicting glacier thinning, while the FSSA method is more suitable for predicting future glacier extent. A possible application of the larger time-step sizes allowed for by the FSSA is for spin-up simulations, where relatively fast-changing climate data can be incorporated on short timescales, while the slow-changing velocity field is updated over larger timescales.</p

Linking bedrock discontinuities to glacial quarrying

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

High-resolution ice thickness and bed topography of a land-terminating section of the Greenland Ice Sheet

We present ice thickness and bed topography maps with a high spatial resolution (250–500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology. The combined data set is freely available at doi:10.1594/pangaea.830314

High resolution ice thickness, bed topography, and roughness of a land terminating section of the western Greenland Ice Sheet

We present ice thickness and bed topography maps with a high spatial resolution (250–500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology. The combined data set is freely available at doi:10.1594/pangaea.830314

Theoretical, contemporary observational and palaeo-perspectives on ice sheet hydrology: processes and products

publisher: Elsevier articletitle: Theoretical, contemporary observational and palaeo-perspectives on ice sheet hydrology: Processes and products journaltitle: Earth-Science Reviews articlelink: http://dx.doi.org/10.1016/j.earscirev.2016.01.010 content_type: article copyright: Copyright © 2016 Elsevier B.V. All rights reserved

Linking bedrock discontinuities to glacial quarrying

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.This article is published as Woodard, J. B., L. K. Zoet, N. R. Iverson, and C. Helanow. "Linking bedrock discontinuities to glacial quarrying." Annals of Glaciology 60, no. 80 (2019): 66-72. doi: 10.1017/aog.2019.36.</p