A double continuum hydrological model for glacier applications

Peer reviewedPublisher PD

Simulated retreat of Jakobshavn Isbræ since the Little Ice Age controlled by geometry

Rapid retreat of Greenland's marine-terminating glaciers coincides with regional warming trends, which have broadly been used to explain these rapid changes. However, outlet glaciers within similar climate regimes experience widely contrasting retreat patterns, suggesting that the local fjord geometry could be an important additional factor. To assess the relative role of climate and fjord geometry, we use the retreat history of Jakobshavn Isbræ, West Greenland, since the Little Ice Age (LIA) maximum in 1850 as a baseline for the parameterization of a depth- and width-integrated ice flow model. The impact of fjord geometry is isolated by using a linearly increasing climate forcing since the LIA and testing a range of simplified geometries.We find that the total length of retreat is determined by external factors – such as hydrofracturing, submarine melt and buttressing by sea ice – whereas the retreat pattern is governed by the fjord geometry. Narrow and shallow areas provide pinning points and cause delayed but rapid retreat without additional climate warming, after decades of grounding line stability. We suggest that these geometric pinning points may be used to locate potential sites for moraine formation and to predict the long-term response of the glacier. As a consequence, to assess the impact of climate on the retreat history of a glacier, each system has to be analyzed with knowledge of its historic retreat and the local fjord geometry.</p

Marine ice sheet dynamics: Hysteresis and neutral equilibrium

International audienceThe stability of marine ice sheets and outlet glaciers is mostly controlled by the dynamics of their grounding line, i.e., where the bottom contact of the ice changes from bedrock or till to ocean water. The last report of the Intergovernmental Panel on Climate Change has clearly underlined the poor ability of models to capture the dynamics of outlet glaciers. Here we present computations of grounding line dynamics on the basis of numerical solutions of the full Stokes equations for ice velocity, coupled with the evolution of the air ice– and sea ice–free interfaces. The grounding line position is determined by solving the contact problem between the ice and a rigid bedrock using the finite element code Elmer. Results of the simulations show that marine ice sheets are unstable on upsloping beds and undergo hysteresis under perturbation of ice viscosity, confirming conclusions from boundary layer theory. The present approach also indicates that a 2-D unconfined marine ice sheet sliding over a downsloping bedrock does not exhibit neutral equilibrium. It is shown that mesh resolution around the grounding line is a crucial issue. A very fine grid size (<100 m spacing) is needed in order to achieve consistent results

Capabilities and performance of Elmer/Ice, a new-generation ice sheet model

The Fourth IPCC Assessment Report concluded that ice sheet flow models, in their current state, were unable to provide accurate forecast for the increase of polar ice sheet discharge and the associated contribution to sea level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice flow models, and as a result a significant number of new ice sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the open-source multi-physics code Elmer. It was one of the first full-Stokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger-scale problems, earning the status of an ice sheet model. Here, we summarise almost 10 yr of development performed by different groups. Elmer/Ice solves the full-Stokes equations, for isotropic but also anisotropic ice rheology, resolves the grounding line dynamics as a contact problem, and contains various basal friction laws. Derived fields, like the age of the ice, the strain rate or stress, can also be computed. Elmer/Ice includes two recently proposed inverse methods to infer badly known parameters. Elmer is a highly parallelised code thanks to recent developments and the implementation of a block preconditioned solver for the Stokes system. In this paper, all these components are presented in detail, as well as the numerical performance of the Stokes solver and developments planned for the future

A modeling study of the effect of runoff variability on the effective pressure beneath Russell Glacier, West Greenland

Basal sliding is a main control on glacier flow primarily driven by water pressure at the glacier base. The ongoing increase in surface melting of the Greenland Ice Sheet warrants an examination of its impact on basal water pressure and in turn on basal sliding. Here we examine the case of Russell Glacier, in West Greenland, where an extensive set of observations has been collected. These observations suggest that the recent increase in melt has had an equivocal impact on the annual velocity, with stable flow on the lower part of the drainage basin but accelerated flow above the Equilibrium Line Altitude (ELA). These distinct behaviors have been attributed to different evolutions of the subglacial draining system during and after the melt season. Here we use a high-resolution subglacial hydrological model forced by reconstructed surface runoff for the period 2008 to 2012 to investigate the cause of these distinct behaviors. We find that the increase in meltwater production at low elevation yields a more efficient drainage system compatible with the observed stagnation of the mean annual flow below the ELA. At higher elevation, the model indicates that the drainage system is mostly inefficient and is therefore strongly sensitive to an increase in meltwater availability, which is consistent with the observed increase in ice velocity

Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP

Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving “shelfy stream” models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients

Ice Thickness and Bed Elevation of the Northern and Southern Patagonian Icefields

International audienc

Two independent methods for mapping the grounding line of an outlet glacier - An example from the Astrolabe Glacier, Terre Adélie, Antarctica

The grounding line is a key element of coastal outlet glaciers, acting on their dynamics. Accurately knowing its position is fundamental for both modelling the glacier dynamics and establishing a benchmark for later change detection. Here we map the grounding line of the Astrolabe Glacier in East Antarctica (66°41' S, 140°05' E), using both hydrostatic and tidal methods. The first method is based on new surface and ice thickness data from which the line of buoyant floatation is found. The second method uses kinematic GPS measurements of the tidal response of the ice surface. By detecting the transitions where the ice starts to move vertically in response to the tidal forcing we determine control points for the grounding line position along GPS profiles. Employing a two-dimensional elastic plate model, we compute the rigid short-term behaviour of the ice plate and estimate the correction required to compare the kinematic GPS control points with the previously determined line of floatation. These two approaches show consistency and lead us to propose a grounding line for the Astrolabe Glacier that significantly deviates from the lines obtained so far from satellite imagery. © Author(s) 2014. CC Attribution 3.0 License