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

2013 ◽  
Vol 6 (1) ◽  
pp. 1689-1741 ◽  
Author(s):  
O. Gagliardini ◽  
T. Zwinger ◽  
F. Gillet-Chaulet ◽  
G. Durand ◽  
L. Favier ◽  
...  
Keyword(s):  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Element Model ◽  
Current Increase ◽  
Ice Flow ◽  
Flow Models ◽  
Flow Problems ◽  
The Status ◽  
And Performance ◽  
New Generation

Abstract. The Fourth IPCC Assessment Report concluded that ice-sheet flow models are unable to forecast the current 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. We present the components already included in Elmer/Ice, its numerical performance, selected applications, as well as developments planned for the future.

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

2013 ◽  
Vol 6 (4) ◽  
pp. 1299-1318 ◽  
Author(s):  
O. Gagliardini ◽  
T. Zwinger ◽  
F. Gillet-Chaulet ◽  
G. Durand ◽  
L. Favier ◽  
...  
Keyword(s):  
Stokes Equations ◽  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Element Model ◽  
Ice Flow ◽  
Flow Models ◽  
The Status ◽  
Recent Developments ◽  
And Performance ◽  
New Generation

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

scholarly journals Sustained high basal motion of the Greenland ice sheet revealed by borehole deformation

2014 ◽  
Vol 60 (222) ◽  
pp. 647-660 ◽  
Author(s):  
Claudia Ryser ◽  
Martin P. Lüthi ◽  
Lauren C. Andrews ◽  
Matthew J. Hoffman ◽  
Ginny A. Catania ◽  
...  
Keyword(s):  
Dynamical Behavior ◽  
Ice Sheet ◽  
Stress Transfer ◽  
Greenland Ice Sheet ◽  
Good Choice ◽  
Surface Velocity ◽  
Ice Streams ◽  
Ice Flow ◽  
Flow Models ◽  
Dependent Flow

AbstractIce deformation and basal motion characterize the dynamical behavior of the Greenland ice sheet (GrIS). We evaluate the contribution of basal motion from ice deformation measurements in boreholes drilled to the bed at two sites in the western marginal zone of the GrIS. We find a sustained high amount of basal motion contribution to surface velocity of 44–73% in winter, and up to 90% in summer. Measured ice deformation rates show an unexpected variation with depth that can be explained with the help of an ice-flow model as a consequence of stress transfer from slippery to sticky areas. This effect necessitates the use of high-order ice-flow models, not only in regions of fast-flowing ice streams but in all temperate-based areas of the GrIS. The agreement between modeled and measured deformation rates confirms that the recommended values of the temperature-dependent flow rate factor A are a good choice for ice-sheet models.

scholarly journals Greenland ice sheet contribution to sea-level rise from a new-generation ice-sheet model

2012 ◽  
Vol 6 (6) ◽  
pp. 1561-1576 ◽  
Author(s):  
F. Gillet-Chaulet ◽  
O. Gagliardini ◽  
H. Seddik ◽  
M. Nodet ◽  
G. Durand ◽  
...  
Keyword(s):  
Mass Balance ◽  
Sea Level Rise ◽  
Sea Level ◽  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Surface Mass Balance ◽  
Ice Flow ◽  
Surface Mass ◽  
Continental Scale ◽  
New Generation

Abstract. Over the last two decades, the Greenland ice sheet (GrIS) has been losing mass at an increasing rate, enhancing its contribution to sea-level rise (SLR). The recent increases in ice loss appear to be due to changes in both the surface mass balance of the ice sheet and ice discharge (ice flux to the ocean). Rapid ice flow directly affects the discharge, but also alters ice-sheet geometry and so affects climate and surface mass balance. Present-day ice-sheet models only represent rapid ice flow in an approximate fashion and, as a consequence, have never explicitly addressed the role of ice discharge on the total GrIS mass balance, especially at the scale of individual outlet glaciers. Here, we present a new-generation prognostic ice-sheet model which reproduces the current patterns of rapid ice flow. This requires three essential developments: the complete solution of the full system of equations governing ice deformation; a variable resolution unstructured mesh to resolve outlet glaciers and the use of inverse methods to better constrain poorly known parameters using observations. The modelled ice discharge is in good agreement with observations on the continental scale and for individual outlets. From this initial state, we investigate possible bounds for the next century ice-sheet mass loss. We run sensitivity experiments of the GrIS dynamical response to perturbations in climate and basal lubrication, assuming a fixed position of the marine termini. We find that increasing ablation tends to reduce outflow and thus decreases the ice-sheet imbalance. In our experiments, the GrIS initial mass (im)balance is preserved throughout the whole century in the absence of reinforced forcing, allowing us to estimate a lower bound of 75 mm for the GrIS contribution to SLR by 2100. In one experiment, we show that the current increase in the rate of ice loss can be reproduced and maintained throughout the whole century. However, this requires a very unlikely perturbation of basal lubrication. From this result we are able to estimate an upper bound of 140 mm from dynamics only for the GrIS contribution to SLR by 2100.

scholarly journals Sensitivity of Greenland ice sheet projections to spatial resolution in higher-order simulations: the AWI contribution to ISMIP6-Greenland using ISSM

2020 ◽  
Author(s):  
Martin Rückamp ◽  
Heiko Goelzer ◽  
Angelika Humbert
Keyword(s):  
Sea Level Rise ◽  
Sea Level ◽  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Higher Order ◽  
Initial State ◽  
Ice Flow ◽  
Outlet Glacier ◽  
Flow Models ◽  
High Emission

Abstract. Projections of the contribution of the Greenland ice sheet to future sea-level rise include uncertainties primarily due to the imposed climate forcing and the initial state of the ice sheet model. Several state-of-the-art ice flow models are currently being employed on various grid resolutions to estimate future mass changes in the framework of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). Here we investigate the sensitivity to grid resolution on centennial sea-level contributions from the Greenland ice sheet and study the mechanism at play. To this end, we employ the finite-element higher-order ice flow model ISSM and conduct experiments with four different horizontal resolutions, namely 4, 2, 1 and 0.75 km. We run the simulation based on the ISMIP6 core GCM MIROC5 under the high emission scenario RCP8.5 and consider both atmospheric and oceanic forcing in full and separate scenarios. Under the full scenarios, finer simulations unveil up to ~5 % more sea-level rise compared to the coarser resolution. The sensitivity depends on the magnitude of outlet glacier retreat, which is implemented as a series of retreat masks following the ISMIP6 protocol. Without imposed retreat under atmosphere-only forcing, the resolution dependency exhibits an opposite behaviour with about ~ 5 % more sea-level contribution in the coarser resolution. The sea-level contribution indicates a converging behaviour

scholarly journals Greenland Ice Sheet contribution to sea-level rise from a new-generation ice-sheet model

2012 ◽  
Vol 6 (4) ◽  
pp. 2789-2826 ◽  
Author(s):  
F. Gillet-Chaulet ◽  
O. Gagliardini ◽  
H. Seddik ◽  
M. Nodet ◽  
G. Durand ◽  
...  
Keyword(s):  
Mass Balance ◽  
Sea Level Rise ◽  
Sea Level ◽  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Surface Mass Balance ◽  
Ice Flow ◽  
Surface Mass ◽  
Continental Scale ◽  
New Generation

Abstract. Over the last two decades, the Greenland Ice Sheet (GrIS) has been losing mass at an increasing rate, enhancing its contribution to sea-level rise. The recent increases in ice loss appear to be due to changes in both the surface mass balance of the ice sheet and ice discharge (ice flux to the ocean). Rapid ice flow directly affects the discharge, but also alters ice-sheet geometry and so affects climate and surface mass balance. The most usual ice-sheet models only represent rapid ice flow in an approximate fashion and, as a consequence, have never explicitly addressed the role of ice discharge on the total GrIS mass balance, especially at the scale of individual outlet glaciers. Here, we present a new-generation prognostic ice-sheet model which reproduces the current patterns of rapid ice flow. This requires three essential developments: the complete solution of the full system of equations governing ice deformation; an unstructured mesh to usefully resolve outlet glaciers and the use of inverse methods to better constrain poorly known parameters using observations. The modelled ice discharge is in good agreement with observations on the continental scale and for individual outlets. By conducting perturbation experiments, we investigate how current ice loss will endure over the next century. Although we find that increasing ablation tends to reduce outflow and on its own has a stabilising effect, if destabilisation processes maintain themselves over time, current increases in the rate of ice loss are likely to continue.

scholarly journals An approach to computing discrete adjoints for MPI-parallelized models applied to Ice Sheet System Model 4.11

2016 ◽  
Vol 9 (11) ◽  
pp. 3907-3918 ◽  
Author(s):  
Eric Larour ◽  
Jean Utke ◽  
Anton Bovin ◽  
Mathieu Morlighem ◽  
Gilberto Perez
Keyword(s):  
Ice Sheet ◽  
Parallel Execution ◽  
System Model ◽  
Ice Flow ◽  
Flow Models ◽  
Adjoint State ◽  
Expected Performance ◽  
Polar Ice Sheets ◽  
New Generation ◽  
Discrete Adjoints

Abstract. Within the framework of sea-level rise projections, there is a strong need for hindcast validation of the evolution of polar ice sheets in a way that tightly matches observational records (from radar, gravity, and altimetry observations mainly). However, the computational requirements for making hindcast reconstructions possible are severe and rely mainly on the evaluation of the adjoint state of transient ice-flow models. Here, we look at the computation of adjoints in the context of the NASA/JPL/UCI Ice Sheet System Model (ISSM), written in C++ and designed for parallel execution with MPI. We present the adaptations required in the way the software is designed and written, but also generic adaptations in the tools facilitating the adjoint computations. We concentrate on the use of operator overloading coupled with the AdjoinableMPI library to achieve the adjoint computation of the ISSM. We present a comprehensive approach to (1) carry out type changing through the ISSM, hence facilitating operator overloading, (2) bind to external solvers such as MUMPS and GSL-LU, and (3) handle MPI-based parallelism to scale the capability. We demonstrate the success of the approach by computing sensitivities of hindcast metrics such as the misfit to observed records of surface altimetry on the northeastern Greenland Ice Stream, or the misfit to observed records of surface velocities on Upernavik Glacier, central West Greenland. We also provide metrics for the scalability of the approach, and the expected performance. This approach has the potential to enable a new generation of hindcast-validated projections that make full use of the wealth of datasets currently being collected, or already collected, in Greenland and Antarctica.

scholarly journals Sensitivity of Greenland Ice Sheet Projections to Model Formulations

2013 ◽  
Vol 59 (216) ◽  
pp. 733-749 ◽  
Author(s):  
H. Goelzer ◽  
P. Huybrechts ◽  
J.J. Fürst ◽  
F.M. Nick ◽  
M.L. Andersen ◽  
...  
Keyword(s):  
Regional Climate Model ◽  
Sea Level ◽  
Flow Model ◽  
Climate Model ◽  
Regional Climate ◽  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Ice Flow ◽  
Outlet Glacier ◽  
Glacier Dynamics

AbstractPhysically based projections of the Greenland ice sheet contribution to future sea-level change are subject to uncertainties of the atmospheric and oceanic climatic forcing and to the formulations within the ice flow model itself. Here a higher-order, three-dimensional thermomechanical ice flow model is used, initialized to the present-day geometry. The forcing comes from a high-resolution regional climate model and from a flowline model applied to four individual marine-terminated glaciers, and results are subsequently extended to the entire ice sheet. The experiments span the next 200 years and consider climate scenario SRES A1B. The surface mass-balance (SMB) scheme is taken either from a regional climate model or from a positive-degree-day (PDD) model using temperature and precipitation anomalies from the underlying climate models. Our model results show that outlet glacier dynamics only account for 6–18% of the sea-level contribution after 200 years, confirming earlier findings that stress the dominant effect of SMB changes. Furthermore, interaction between SMB and ice discharge limits the importance of outlet glacier dynamics with increasing atmospheric forcing. Forcing from the regional climate model produces a 14–31 % higher sea-level contribution compared to a PDD model run with the same parameters as for IPCC AR4.

scholarly journals Winter drainage of surface lakes on the Greenland Ice Sheet from Sentinel-1 SAR imagery

2021 ◽  
Vol 15 (3) ◽  
pp. 1587-1606
Author(s):  
Corinne L. Benedek ◽  
Ian C. Willis
Keyword(s):  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Shape From Shading ◽  
Lake Area ◽  
Ice Flow ◽  
Surface Mass ◽  
Flow Acceleration ◽  
Optical Imagery ◽  
Sar Imagery ◽  
Lake Drainage

Abstract. Surface lakes on the Greenland Ice Sheet play a key role in its surface mass balance, hydrology and biogeochemistry. They often drain rapidly in the summer via hydrofracture, which delivers lake water to the ice sheet base over timescales of hours to days and then can allow meltwater to reach the base for the rest of the summer. Rapid lake drainage, therefore, influences subglacial drainage evolution; water pressures; ice flow; biogeochemical activity; and ultimately the delivery of water, sediments and nutrients to the ocean. It has generally been assumed that rapid lake drainage events are confined to the summer, as this is typically when observations are made using satellite optical imagery. Here we develop a method to quantify backscatter changes in satellite radar imagery, which we use to document the drainage of six different lakes during three winters (2014/15, 2015/16 and 2016/17) in fast-flowing parts of the Greenland Ice Sheet. Analysis of optical imagery from before and after the three winters supports the radar-based evidence for winter lake drainage events and also provides estimates of lake drainage volumes, which range between 0.000046 ± 0.000017 and 0.0200 ± 0.002817 km3. For three of the events, optical imagery allows repeat photoclinometry (shape from shading) calculations to be made showing mean vertical collapse of the lake surfaces ranging between 1.21 ± 1.61 and 7.25 ± 1.61 m and drainage volumes of 0.002 ± 0.002968 to 0.044 ± 0.009858 km3. For one of these three, time-stamped ArcticDEM strips allow for DEM differencing, which demonstrates a mean collapse depth of 2.17 ± 0.28 m across the lake area. The findings show that lake drainage can occur in the winter in the absence of active surface melt and notable ice flow acceleration, which may have important implications for subglacial hydrology and biogeochemical processes.

scholarly journals A comparison of the performance of depth-integrated ice-dynamics solvers

2021 ◽  
Author(s):  
Alexander Robinson ◽  
Daniel Goldberg ◽  
William H. Lipscomb
Keyword(s):  
Numerical Stability ◽  
Ice Sheet ◽  
Model Performance ◽  
Greenland Ice Sheet ◽  
Grid Resolution ◽  
Ice Dynamics ◽  
Fast Solvers ◽  
Ice Flow ◽  
Hybrid Solver ◽  
Constant Viscosity

Abstract. In the last decade, the number of ice-sheet models has increased substantially, in line with the growth of the glaciological community. These models use solvers based on different approximations of ice dynamics. In particular, several depth-integrated dynamics approximations have emerged as fast solvers capable of resolving the relevant physics of ice sheets at the continen- tal scale. However, the numerical stability of these schemes has not been studied systematically to evaluate their effectiveness in practice. Here we focus on three such solvers, the so-called Hybrid, L1L2-SIA and DIVA solvers, as well as the well-known SIA and SSA solvers as boundary cases. We investigate the numerical stability of these solvers as a function of grid resolution and the state of the ice sheet. Under simplified conditions with constant viscosity, the maximum stable timestep of the Hybrid solver, like the SIA solver, has a quadratic dependence on grid resolution. In contrast, the DIVA solver has a maximum timestep that is independent of resolution, like the SSA solver. Analysis indicates that the L1L2-SIA solver should behave similarly, but in practice, the complexity of its implementation can make it difficult to maintain stability. In realistic simulations of the Greenland ice sheet with a non-linear rheology, the DIVA and SSA solvers maintain superior numerical stability, while the SIA, Hybrid and L1L2-SIA solvers show markedly poorer performance. At a grid resolution of ∆x = 4 km, the DIVA solver runs approximately 15 times faster than the Hybrid and L1L2-SIA solvers. Our analysis shows that as resolution increases, the ice-dynamics solver can act as a bottleneck to model performance. The DIVA solver emerges as a clear outlier in terms of both model performance and its representation of the ice-flow physics itself.

Sign in / Sign up

Export Citation Format

Share Document