Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
The Cryosphere, 9, 2183-2200, 2015
© Author(s) 2015. This work is distributed under
the Creative Commons Attribution 3.0 License.
Research article
20 Nov 2015
A new methodology to simulate subglacial deformation of water-saturated granular material
A. Damsgaard1, D. L. Egholm1, J. A. Piotrowski1, S. Tulaczyk2, N. K. Larsen1,3, and C. F. Brædstrup1 1Department of Geoscience, Aarhus University, 8000 Aarhus C, Denmark
2Department of Earth and Planetary Sciences, University of California, Santa Cruz, California 95064, USA
3Centre for GeoGenetics, Natural History Museum of Denmark, University of Copenhagen, 1350 Copenhagen K, Denmark
Abstract. The dynamics of glaciers are to a large degree governed by processes operating at the ice–bed interface, and one of the primary mechanisms of glacier flow over soft unconsolidated sediments is subglacial deformation. However, it has proven difficult to constrain the mechanical response of subglacial sediment to the shear stress of an overriding glacier. In this study, we present a new methodology designed to simulate subglacial deformation using a coupled numerical model for computational experiments on grain-fluid mixtures. The granular phase is simulated on a per-grain basis by the discrete element method. The pore water is modeled as a compressible Newtonian fluid without inertia. The numerical approach allows close monitoring of the internal behavior under a range of conditions.

Our computational experiments support the findings of previous studies where the rheology of a slowly deforming water-saturated granular bed in the steady state generally conforms to the rate-independent plastic rheology. Before this so-called critical state, deformation is in many cases accompanied by volumetric changes as grain rearrangement in active shear zones changes the local porosity. For previously consolidated beds porosity increases can cause local pore-pressure decline, dependent on till permeability and shear rate. We observe that the pore-water pressure reduction strengthens inter-granular contacts, which results in increased shear strength of the granular material. In contrast, weakening takes place when shear deformation causes consolidation of dilated sediments or during rapid fabric development. Both processes of strengthening and weakening depend inversely on the sediment permeability and are transient phenomena tied to the porosity changes during the early stages of shear.

We find that the transient strengthening and weakening in turn influences the distribution of shear strain in the granular bed. Dilatant strengthening has the ability to distribute strain during early deformation to large depths, if sediment dilatancy causes the water pressure at the ice–bed interface to decline. Oppositely, if the ice–bed interface is hydrologically stable the strengthening process is minimal and instead causes shallow deformation. The depth of deformation in subglacial beds thus seems to be governed by not only local grain and pore-water feedbacks but also larger-scale hydrological properties at the ice base.

Citation: Damsgaard, A., Egholm, D. L., Piotrowski, J. A., Tulaczyk, S., Larsen, N. K., and Brædstrup, C. F.: A new methodology to simulate subglacial deformation of water-saturated granular material, The Cryosphere, 9, 2183-2200,, 2015.
Publications Copernicus
Short summary
This paper details a new algorithm for performing computational experiments of subglacial granular deformation. The numerical approach allows detailed studies of internal sediment and pore-water dynamics under shear. Feedbacks between sediment grains and pore water can cause rate-dependent strengthening, which additionally contributes to the plastic shear strength of the granular material. Hardening can stabilise patches of the subglacial beds with implications for landform development.
This paper details a new algorithm for performing computational experiments of subglacial...