Articles | Volume 14, issue 9
Geosci. Model Dev., 14, 5843–5861, 2021
https://doi.org/10.5194/gmd-14-5843-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Geosci. Model Dev., 14, 5843–5861, 2021
https://doi.org/10.5194/gmd-14-5843-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Model description paper 24 Sep 2021
Model description paper | 24 Sep 2021
fenics_ice 1.0: a framework for quantifying initialization uncertainty for time-dependent ice sheet models
Conrad P. Koziol et al.
Related authors
Conrad P. Koziol and Neil Arnold
The Cryosphere, 12, 971–991, https://doi.org/10.5194/tc-12-971-2018, https://doi.org/10.5194/tc-12-971-2018, 2018
Short summary
Short summary
We model the summer acceleration of ice velocities at a land-terminating margin of the Greenland Ice Sheet. Model results compare favourably against GPS data, reflecting positively on the model components and the datasets used. When we run the model into the future, we find that summer velocities increase with increasing levels of surface melt but that changes in annual velocities may be limited.
Conrad P. Koziol and Neil Arnold
The Cryosphere, 11, 2783–2797, https://doi.org/10.5194/tc-11-2783-2017, https://doi.org/10.5194/tc-11-2783-2017, 2017
Short summary
Short summary
We develop a new ice sheet model and couple it to an existing subglacial hydrology model. A workflow for initializing the coupled model at the start of summer is proposed and demonstrated on the Russell Glacier area of Western Greenland. This is a first step towards modelling ice velocities during the summer.
Helen Ockenden, Robert G. Bingham, Andrew Curtis, and Daniel Goldberg
The Cryosphere Discuss., https://doi.org/10.5194/tc-2021-287, https://doi.org/10.5194/tc-2021-287, 2021
Preprint under review for TC
Short summary
Short summary
Hills and valleys hidden under the ice of Thwaites Glacier have an impact on ice flow and future ice loss, but there are not many three dimensional observations of their location or size. We apply a mathematical theory to new high resolution observations of the ice surface to predict the bed topography beneath the ice. There is a good correlation with ice-penetrating radar observations. The method may be useful in areas with few direct observations, or as a further constraint for other methods.
Alexander Robinson, Daniel Goldberg, and William H. Lipscomb
The Cryosphere Discuss., https://doi.org/10.5194/tc-2021-239, https://doi.org/10.5194/tc-2021-239, 2021
Revised manuscript accepted for TC
Short summary
Short summary
Here we investigate the numerical stability of several commonly used methods, in order to determine which of them are both capable of resolving the complex physics of the ice flow and are also computationally efficient. We find that the so-called DIVA solver outperforms the others. Its representation of the physics is consistent with more complex methods, while it remains computationally efficient at high resolution.
Jowan M. Barnes, Thiago Dias dos Santos, Daniel Goldberg, G. Hilmar Gudmundsson, Mathieu Morlighem, and Jan De Rydt
The Cryosphere, 15, 1975–2000, https://doi.org/10.5194/tc-15-1975-2021, https://doi.org/10.5194/tc-15-1975-2021, 2021
Short summary
Short summary
Some properties of ice flow models must be initialised using observed data before they can be used to produce reliable predictions of the future. Different models have different ways of doing this, and the process is generally seen as being specific to an individual model. We compare the methods used by three different models and show that they produce similar outputs. We also demonstrate that the outputs from one model can be used in other models without introducing large uncertainties.
Stephen L. Cornford, Helene Seroussi, Xylar S. Asay-Davis, G. Hilmar Gudmundsson, Rob Arthern, Chris Borstad, Julia Christmann, Thiago Dias dos Santos, Johannes Feldmann, Daniel Goldberg, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, Gunter Leguy, William H. Lipscomb, Nacho Merino, Gaël Durand, Mathieu Morlighem, David Pollard, Martin Rückamp, C. Rosie Williams, and Hongju Yu
The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, https://doi.org/10.5194/tc-14-2283-2020, 2020
Short summary
Short summary
We present the results of the third Marine Ice Sheet Intercomparison Project (MISMIP+). MISMIP+ is one in a series of exercises that test numerical models of ice sheet flow in simple situations. This particular exercise concentrates on the response of ice sheet models to the thinning of their floating ice shelves, which is of interest because numerical models are currently used to model the response to contemporary and near-future thinning in Antarctic ice shelves.
Conrad P. Koziol and Neil Arnold
The Cryosphere, 12, 971–991, https://doi.org/10.5194/tc-12-971-2018, https://doi.org/10.5194/tc-12-971-2018, 2018
Short summary
Short summary
We model the summer acceleration of ice velocities at a land-terminating margin of the Greenland Ice Sheet. Model results compare favourably against GPS data, reflecting positively on the model components and the datasets used. When we run the model into the future, we find that summer velocities increase with increasing levels of surface melt but that changes in annual velocities may be limited.
Conrad P. Koziol and Neil Arnold
The Cryosphere, 11, 2783–2797, https://doi.org/10.5194/tc-11-2783-2017, https://doi.org/10.5194/tc-11-2783-2017, 2017
Short summary
Short summary
We develop a new ice sheet model and couple it to an existing subglacial hydrology model. A workflow for initializing the coupled model at the start of summer is proposed and demonstrated on the Russell Glacier area of Western Greenland. This is a first step towards modelling ice velocities during the summer.
Daniel N. Goldberg, Sri Hari Krishna Narayanan, Laurent Hascoet, and Jean Utke
Geosci. Model Dev., 9, 1891–1904, https://doi.org/10.5194/gmd-9-1891-2016, https://doi.org/10.5194/gmd-9-1891-2016, 2016
Short summary
Short summary
Geophysical adjoint models are powerful tools, allowing sensitivity studies that are not possible otherwise, and enabling optimized fit of models to observing data sets. The complexity involved requires the use of algorithmic differentiation (AD) software, but AD adjoint calculation for ice models can be slow, with prohibitive memory requirements. In this paper, we present a method to improve the performance of ice model adjoint generation, in terms of timing, memory load, and accuracy.
D. N. Goldberg, P. Heimbach, I. Joughin, and B. Smith
The Cryosphere, 9, 2429–2446, https://doi.org/10.5194/tc-9-2429-2015, https://doi.org/10.5194/tc-9-2429-2015, 2015
Short summary
Short summary
We calibrate a time-dependent ice model through optimal fit to transient observations of surface elevation and velocity, a novel procedure in glaciology and in particular for an ice stream with a dynamic grounding line. We show this procedure gives a level of confidence in model projections that cannot be achieved through more commonly used glaciological data assimilation methods. We show that Smith Glacier is in a state of retreat regardless of climatic forcing for the next several decades.
D. N. Goldberg and P. Heimbach
The Cryosphere, 7, 1659–1678, https://doi.org/10.5194/tc-7-1659-2013, https://doi.org/10.5194/tc-7-1659-2013, 2013
Related subject area
Cryosphere
SNICAR-ADv3: a community tool for modeling spectral snow albedo
STEMMUS-UEB v1.0.0: integrated modeling of snowpack and soil water and energy transfer with three complexity levels of soil physical processes
A versatile method for computing optimized snow albedo from spectrally fixed radiative variables: VALHALLA v1.0
Ice Algae Model Intercomparison Project phase 2 (IAMIP2)
A Gaussian process emulator for simulating ice sheet–climate interactions on a multi-million-year timescale: CLISEMv1.0
SITool (v1.0) – a new evaluation tool for large-scale sea ice simulations: application to CMIP6 OMIP
Mapping high-resolution basal topography of West Antarctica from radar data using non-stationary multiple-point geostatistics (MPS-BedMappingV1)
WIFF1.0: A hybrid machine-learning-based parameterization of Wave-Induced sea-ice Floe Fracture
Development of adjoint-based ocean state estimation for the Amundsen and Bellingshausen seas and ice shelf cavities using MITgcm–ECCO (66j)
Sensitivity of Northern Hemisphere climate to ice–ocean interface heat flux parameterizations
icepack: a new glacier flow modeling package in Python, version 1.0
CAPS v1.0: An improved regional coupled modeling system for Arctic sea ice and climate simulation and prediction
Benefits of sea ice initialization for the interannual-to-decadal climate prediction skill in the Arctic in EC-Earth3
Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain
NEMO-Bohai 1.0: a high-resolution ocean and sea ice modelling system for the Bohai Sea, China
Performance of MAR (v3.11) in simulating the drifting-snow climate and surface mass balance of Adélie Land, East Antarctica
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18
The Utrecht Finite Volume Ice-Sheet Model: UFEMISM (version 1.0)
PERICLIMv1.0: a model deriving palaeo-air temperatures from thaw depth in past permafrost regions
Assessing the simulated soil hydrothermal regime of the active layer from the Noah-MP land surface model (v1.1) in the permafrost regions of the Qinghai–Tibet Plateau
CrocO_v1.0: a particle filter to assimilate snowpack observations in a spatialised framework
A fully coupled Arctic sea-ice–ocean–atmosphere model (ArcIOAM v1.0) based on C-Coupler2: model description and preliminary results
The Framework For Ice Sheet–Ocean Coupling (FISOC) V1.1
Comparison of sea ice kinematics at different resolutions modeled with a grid hierarchy in the Community Earth System Model (version 1.2.1)
Snow profile alignment and similarity assessment for aggregating, clustering, and evaluating snowpack model output for avalanche forecasting
Improvements in one-dimensional grounding-line parameterizations in an ice-sheet model with lateral variations (PSUICE3D v2.1)
Implementation of the RCIP scheme and its performance for 1-D age computations in ice-sheet models
COSIPY v1.3 – an open-source coupled snowpack and ice surface energy and mass balance model
Using Arctic ice mass balance buoys for evaluation of modelled ice energy fluxes
Impact of the ice thickness distribution discretization on the sea ice concentration variability in the NEMO3.6–LIM3 global ocean–sea ice model
Simulating the Early Holocene demise of the Laurentide Ice Sheet with BISICLES (public trunk revision 3298)
Extended enthalpy formulations in the Ice-sheet and Sea-level System Model (ISSM) version 4.17: discontinuous conductivity and anisotropic streamline upwind Petrov–Galerkin (SUPG) method
The Community Firn Model (CFM) v1.0
Description and validation of the ice-sheet model Yelmo (version 1.0)
The Whole Antarctic Ocean Model (WAOM v1.0): Development and Evaluation
Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD
On the calculation of normalized viscous–plastic sea ice stresses
Modelling thermomechanical ice deformation using an implicit pseudo-transient method (FastICE v1.0) based on graphical processing units (GPUs)
Version 1 of a sea ice module for the physics-based, detailed, multi-layer SNOWPACK model
A module to convert spectral to narrowband snow albedo for use in climate models: SNOWBAL v1.2
On the discretization of the ice thickness distribution in the NEMO3.6-LIM3 global ocean–sea ice model
Scientific workflows applied to the coupling of a continuum (Elmer v8.3) and a discrete element (HiDEM v1.0) ice dynamic model
A rapidly converging initialisation method to simulate the present-day Greenland ice sheet using the GRISLI ice sheet model (version 1.3)
Assessment of sub-shelf melting parameterisations using the ocean–ice-sheet coupled model NEMO(v3.6)–Elmer/Ice(v8.3)
CSIB v1 (Canadian Sea-ice Biogeochemistry): a sea-ice biogeochemical model for the NEMO community ocean modelling framework
LIVVkit 2.1: automated and extensible ice sheet model validation
ATAT 1.1, the Automated Timing Accordance Tool for comparing ice-sheet model output with geochronological data
The Open Global Glacier Model (OGGM) v1.1
Description and evaluation of the Community Ice Sheet Model (CISM) v2.1
Mark G. Flanner, Julian B. Arnheim, Joseph M. Cook, Cheng Dang, Cenlin He, Xianglei Huang, Deepak Singh, S. McKenzie Skiles, Chloe A. Whicker, and Charles S. Zender
Geosci. Model Dev., 14, 7673–7704, https://doi.org/10.5194/gmd-14-7673-2021, https://doi.org/10.5194/gmd-14-7673-2021, 2021
Short summary
Short summary
We present the technical formulation and evaluation of a publicly available code and web-based model to simulate the spectral albedo of snow. Our model accounts for numerous features of the snow state and ambient conditions, including the the presence of light-absorbing matter like black and brown carbon, mineral dust, volcanic ash, and snow algae. Carbon dioxide snow, found on Mars, is also represented. The model accurately reproduces spectral measurements of clean and contaminated snow.
Lianyu Yu, Yijian Zeng, and Zhongbo Su
Geosci. Model Dev., 14, 7345–7376, https://doi.org/10.5194/gmd-14-7345-2021, https://doi.org/10.5194/gmd-14-7345-2021, 2021
Short summary
Short summary
We developed an integrated soil–snow–atmosphere model (STEMMUS-UEB) dedicated to the physical description of snow and soil processes with various complexities. With STEMMUS-UEB, we demonstrated that the snowpack affects not only the soil surface moisture conditions (in the liquid and ice phase) and energy-related states (albedo, LE) but also the subsurface soil water and vapor transfer, which contributes to a better understanding of the hydrothermal implications of the snowpack in cold regions.
Florent Veillon, Marie Dumont, Charles Amory, and Mathieu Fructus
Geosci. Model Dev., 14, 7329–7343, https://doi.org/10.5194/gmd-14-7329-2021, https://doi.org/10.5194/gmd-14-7329-2021, 2021
Short summary
Short summary
In climate models, the snow albedo scheme generally calculates only a narrowband or broadband albedo. Therefore, we have developed the VALHALLA method to optimize snow spectral albedo calculations through the determination of spectrally fixed radiative variables. The development of VALHALLA v1.0 with the use of the snow albedo model TARTES and the spectral irradiance model SBDART indicates a considerable reduction in calculation time while maintaining an adequate accuracy of albedo values.
Hakase Hayashida, Meibing Jin, Nadja S. Steiner, Neil C. Swart, Eiji Watanabe, Russell Fiedler, Andrew McC. Hogg, Andrew E. Kiss, Richard J. Matear, and Peter G. Strutton
Geosci. Model Dev., 14, 6847–6861, https://doi.org/10.5194/gmd-14-6847-2021, https://doi.org/10.5194/gmd-14-6847-2021, 2021
Short summary
Short summary
Ice algae are tiny plants like phytoplankton but they grow within sea ice. In polar regions, both phytoplankton and ice algae are the foundation of marine ecosystems and play an important role in taking up carbon dioxide in the atmosphere. However, state-of-the-art climate models typically do not include ice algae, and therefore their role in the climate system remains unclear. This project aims to address this knowledge gap by coordinating a set of experiments using sea-ice–ocean models.
Jonas Van Breedam, Philippe Huybrechts, and Michel Crucifix
Geosci. Model Dev., 14, 6373–6401, https://doi.org/10.5194/gmd-14-6373-2021, https://doi.org/10.5194/gmd-14-6373-2021, 2021
Short summary
Short summary
Ice sheets are an important component of the climate system and interact with the atmosphere through albedo variations and changes in the surface height. On very long timescales, it is impossible to directly couple ice sheet models with climate models and other techniques have to be used. Here we present a novel coupling method between ice sheets and the atmosphere by making use of an emulator to simulate ice sheet–climate interactions for several million years.
Xia Lin, François Massonnet, Thierry Fichefet, and Martin Vancoppenolle
Geosci. Model Dev., 14, 6331–6354, https://doi.org/10.5194/gmd-14-6331-2021, https://doi.org/10.5194/gmd-14-6331-2021, 2021
Short summary
Short summary
This study introduces a new Sea Ice Evaluation Tool (SITool) to evaluate the model skills on the bipolar sea ice simulations by providing performance metrics and diagnostics. SITool is applied to evaluate the CMIP6 OMIP simulations. By changing the atmospheric forcing from CORE-II to JRA55-do data, many aspects of sea ice simulations are improved. SITool will be useful for helping teams managing various versions of a sea ice model or tracking the time evolution of model performance.
Zhen Yin, Chen Zuo, Emma J. MacKie, and Jef Caers
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2021-297, https://doi.org/10.5194/gmd-2021-297, 2021
Preprint under review for GMD
Short summary
Short summary
We provide a geospatial modeling approach to sample training images to fill large scale irregular geophysical data gaps in West Antarctica. Using a repository of global topographic training images, our approach models high-resolution basal topography and quantifies the geospatial uncertainty . It generated multiple high-resolution topographic realizations to investigate the impact of basal topographic uncertainty on critical subglacial hydrological flow patterns associated with ice velocity.
Christopher Horvat and Lettie A. Roach
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2021-281, https://doi.org/10.5194/gmd-2021-281, 2021
Revised manuscript accepted for GMD
Short summary
Short summary
Sea ice is a composite of individual pieces, called floes, ranging in horizontal size from meters to kilometers. Variations in sea ice geometry are often forced by ocean waves, a process that is an important target of global climate models as it affects the rate of sea ice melting. Yet directly simulating these interactions is computationally expensive. We present a neural network-based model of wave-ice fracture that allows models to incorporate their effect with no added computational cost.
Yoshihiro Nakayama, Dimitris Menemenlis, Ou Wang, Hong Zhang, Ian Fenty, and An T. Nguyen
Geosci. Model Dev., 14, 4909–4924, https://doi.org/10.5194/gmd-14-4909-2021, https://doi.org/10.5194/gmd-14-4909-2021, 2021
Short summary
Short summary
High ice shelf melting in the Amundsen Sea has attracted many observational campaigns in the past decade. One method to combine observations with numerical models is the adjoint method. After 20 iterations, the cost function, defined as a sum of the weighted model–data difference, is reduced by 65 % by adjusting initial conditions, atmospheric forcing, and vertical diffusivity. This study demonstrates adjoint-method optimization with explicit representation of ice shelf cavity circulation.
Xiaoxu Shi, Dirk Notz, Jiping Liu, Hu Yang, and Gerrit Lohmann
Geosci. Model Dev., 14, 4891–4908, https://doi.org/10.5194/gmd-14-4891-2021, https://doi.org/10.5194/gmd-14-4891-2021, 2021
Short summary
Short summary
The ice–ocean heat flux is one of the key elements controlling sea ice changes. It motivates our study, which aims to examine the responses of modeled climate to three ice–ocean heat flux parameterizations, including two old approaches that assume one-way heat transport and a new one describing a double-diffusive ice–ocean heat exchange. The results show pronounced differences in the modeled sea ice, ocean, and atmosphere states for the latter as compared to the former two parameterizations.
Daniel R. Shapero, Jessica A. Badgeley, Andrew O. Hoffman, and Ian R. Joughin
Geosci. Model Dev., 14, 4593–4616, https://doi.org/10.5194/gmd-14-4593-2021, https://doi.org/10.5194/gmd-14-4593-2021, 2021
Short summary
Short summary
This paper describes a new software package called "icepack" for modeling the flow of ice sheets and glaciers. Glaciologists use tools like icepack to better understand how ice sheets flow, what role they have played in shaping Earth's climate, and how much sea level rise we can expect in the coming decades to centuries. The icepack package includes several innovations to help researchers describe and solve interesting glaciological problems and to experiment with the underlying model physics.
Chao-Yuan Yang, Jiping Liu, and Dake Chen
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2021-220, https://doi.org/10.5194/gmd-2021-220, 2021
Revised manuscript accepted for GMD
Short summary
Short summary
We present an improved coupled modeling system for Arctic sea ice and climate prediction. We perform Arctic sea ice prediction experiments with improved/updated physical parameterizations, which show better skill in predicting sea ice state as well as atmospheric and oceanic state in the Arctic compared with its predecessor. The improved model also show extended predictive skill of Arctic sea ice after the summer season. This provides added values of this prediction system for decision-making.
Tian Tian, Shuting Yang, Mehdi Pasha Karami, François Massonnet, Tim Kruschke, and Torben Koenigk
Geosci. Model Dev., 14, 4283–4305, https://doi.org/10.5194/gmd-14-4283-2021, https://doi.org/10.5194/gmd-14-4283-2021, 2021
Short summary
Short summary
Three decadal prediction experiments with EC-Earth3 are performed to investigate the impact of ocean, sea ice concentration and thickness initialization, respectively. We find that the persistence of perennial thick ice in the central Arctic can affect the sea ice predictability in its adjacent waters via advection process or wind, despite those regions being seasonally ice free during two recent decades. This has implications for the coming decades as the thinning of Arctic sea ice continues.
Moritz Kreuzer, Ronja Reese, Willem Nicholas Huiskamp, Stefan Petri, Torsten Albrecht, Georg Feulner, and Ricarda Winkelmann
Geosci. Model Dev., 14, 3697–3714, https://doi.org/10.5194/gmd-14-3697-2021, https://doi.org/10.5194/gmd-14-3697-2021, 2021
Short summary
Short summary
We present the technical implementation of a coarse-resolution coupling between an ice sheet model and an ocean model that allows one to simulate ice–ocean interactions at timescales from centuries to millennia. As ice shelf cavities cannot be resolved in the ocean model at coarse resolution, we bridge the gap using an sub-shelf cavity module. It is shown that the framework is computationally efficient, conserves mass and energy, and can produce a stable coupled state under present-day forcing.
Yu Yan, Wei Gu, Andrea M. U. Gierisch, Yingjun Xu, and Petteri Uotila
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2021-100, https://doi.org/10.5194/gmd-2021-100, 2021
Revised manuscript accepted for GMD
Short summary
Short summary
In this study, we developed NEMO-Bohai, an ocean-ice model for the Bohai Sea, China. This study presented the scientific design and technical choices of the parameterizations for the NEMO-Bohai model. The model was calibrated and evaluated with in situ and satellite observations of ocean and sea ice. Nemo-Bohai is intended to be a useful tool for long-term ocean and ice simulations as well as the ocean and climate change studies.
Charles Amory, Christoph Kittel, Louis Le Toumelin, Cécile Agosta, Alison Delhasse, Vincent Favier, and Xavier Fettweis
Geosci. Model Dev., 14, 3487–3510, https://doi.org/10.5194/gmd-14-3487-2021, https://doi.org/10.5194/gmd-14-3487-2021, 2021
Short summary
Short summary
This paper presents recent developments in the drifting-snow scheme of the regional climate model MAR and its application to simulate drifting snow and the surface mass balance of Adélie Land in East Antarctica. The model is extensively described and evaluated against a multi-year drifting-snow dataset and surface mass balance estimates available in the area. The model sensitivity to input parameters and improvements over a previously published version are also assessed.
Thiago Dias dos Santos, Mathieu Morlighem, and Hélène Seroussi
Geosci. Model Dev., 14, 2545–2573, https://doi.org/10.5194/gmd-14-2545-2021, https://doi.org/10.5194/gmd-14-2545-2021, 2021
Short summary
Short summary
Numerical models are routinely used to understand the past and future behavior of ice sheets in response to climate evolution. As is always the case with numerical modeling, one needs to minimize biases and numerical artifacts due to the choice of numerical scheme employed in such models. Here, we assess different numerical schemes in time-dependent simulations of ice sheets. We also introduce a new parameterization for the driving stress, the force that drives the ice sheet flow.
Constantijn J. Berends, Heiko Goelzer, and Roderik S. W. van de Wal
Geosci. Model Dev., 14, 2443–2470, https://doi.org/10.5194/gmd-14-2443-2021, https://doi.org/10.5194/gmd-14-2443-2021, 2021
Short summary
Short summary
The largest uncertainty in projections of sea-level rise comes from ice-sheet retreat. To better understand how these ice sheets respond to the changing climate, ice-sheet models are used, which must be able to reproduce both their present and past evolution. We have created a model that is fast enough to simulate an ice sheet at a high resolution over the course of an entire 120 000-year glacial cycle. This allows us to study processes that cannot be captured by lower-resolution models.
Tomáš Uxa, Marek Křížek, and Filip Hrbáček
Geosci. Model Dev., 14, 1865–1884, https://doi.org/10.5194/gmd-14-1865-2021, https://doi.org/10.5194/gmd-14-1865-2021, 2021
Short summary
Short summary
We present a simple model that derives palaeo-air temperature characteristics related to the palaeo-active-layer thickness, which can be recognized using many relict periglacial features found in past permafrost regions. Its evaluation against modern temperature records and an experimental palaeo-air temperature reconstruction showed relatively high model accuracy, which suggests that it could become a useful tool for reconstructing Quaternary palaeo-environments.
Xiangfei Li, Tonghua Wu, Xiaodong Wu, Jie Chen, Xiaofan Zhu, Guojie Hu, Ren Li, Yongping Qiao, Cheng Yang, Junming Hao, Jie Ni, and Wensi Ma
Geosci. Model Dev., 14, 1753–1771, https://doi.org/10.5194/gmd-14-1753-2021, https://doi.org/10.5194/gmd-14-1753-2021, 2021
Short summary
Short summary
In this study, an ensemble simulation of 55296 scheme combinations for at a typical permafrost site on the Qinghai–Tibet Plateau (QTP) was conducted. The general performance of the Noah-MP model for snow cover events (SCEs), soil temperature (ST) and soil liquid water content (SLW) was assessed, and the sensitivities of parameterization schemes at different depths were investigated. We show that Noah-MP tends to overestimate SCEs and underestimate ST and topsoil SLW on the QTP.
Bertrand Cluzet, Matthieu Lafaysse, Emmanuel Cosme, Clément Albergel, Louis-François Meunier, and Marie Dumont
Geosci. Model Dev., 14, 1595–1614, https://doi.org/10.5194/gmd-14-1595-2021, https://doi.org/10.5194/gmd-14-1595-2021, 2021
Short summary
Short summary
In the mountains, the combination of large model error and observation sparseness is a challenge for data assimilation. Here, we develop two variants of the particle filter (PF) in order to propagate the information content of observations into unobserved areas. By adjusting observation errors or exploiting background correlation patterns, we demonstrate the potential for partial observations of snow depth and surface reflectance to improve model accuracy with the PF in an idealised setting.
Shihe Ren, Xi Liang, Qizhen Sun, Hao Yu, L. Bruno Tremblay, Bo Lin, Xiaoping Mai, Fu Zhao, Ming Li, Na Liu, Zhikun Chen, and Yunfei Zhang
Geosci. Model Dev., 14, 1101–1124, https://doi.org/10.5194/gmd-14-1101-2021, https://doi.org/10.5194/gmd-14-1101-2021, 2021
Short summary
Short summary
Sea ice plays a crucial role in global energy and water budgets. To get a better simulation of sea ice, we coupled a sea ice model with an atmospheric and ocean model to form a fully coupled system. The sea ice simulation results of this coupled system demonstrated that a two-way coupled model has better performance in terms of sea ice, especially in summer. This indicates that sea-ice–ocean–atmosphere interaction plays a crucial role in controlling Arctic summertime sea ice distribution.
Rupert Gladstone, Benjamin Galton-Fenzi, David Gwyther, Qin Zhou, Tore Hattermann, Chen Zhao, Lenneke Jong, Yuwei Xia, Xiaoran Guo, Konstantinos Petrakopoulos, Thomas Zwinger, Daniel Shapero, and John Moore
Geosci. Model Dev., 14, 889–905, https://doi.org/10.5194/gmd-14-889-2021, https://doi.org/10.5194/gmd-14-889-2021, 2021
Short summary
Short summary
Retreat of the Antarctic ice sheet, and hence its contribution to sea level rise, is highly sensitive to melting of its floating ice shelves. This melt is caused by warm ocean currents coming into contact with the ice. Computer models used for future ice sheet projections are not able to realistically evolve these melt rates. We describe a new coupling framework to enable ice sheet and ocean computer models to interact, allowing projection of the evolution of melt and its impact on sea level.
Shiming Xu, Jialiang Ma, Lu Zhou, Yan Zhang, Jiping Liu, and Bin Wang
Geosci. Model Dev., 14, 603–628, https://doi.org/10.5194/gmd-14-603-2021, https://doi.org/10.5194/gmd-14-603-2021, 2021
Short summary
Short summary
A multi-resolution tripolar grid hierarchy is constructed and integrated in CESM (version 1.2.1). The resolution range includes 0.45, 0.15, and 0.05°. Based on atmospherically forced sea ice experiments, the model simulates reasonable sea ice kinematics and scaling properties. Landfast ice thickness can also be systematically shifted due to non-convergent solutions to an
elastic–viscous–plastic (EVP) model. This work is a framework for multi-scale modeling of the ocean and sea ice with CESM.
Florian Herla, Simon Horton, Patrick Mair, and Pascal Haegeli
Geosci. Model Dev., 14, 239–258, https://doi.org/10.5194/gmd-14-239-2021, https://doi.org/10.5194/gmd-14-239-2021, 2021
Short summary
Short summary
The adoption of snowpack models in support of avalanche forecasting has been limited. To promote their operational application, we present a numerical method for processing multivariate snow stratigraphy profiles of mixed data types. Our algorithm enables applications like dynamical grouping and summarizing of model simulations, model evaluation, and data assimilation. By emulating the human analysis process, our approach will allow forecasters to familiarly interact with snowpack simulations.
David Pollard and Robert M. DeConto
Geosci. Model Dev., 13, 6481–6500, https://doi.org/10.5194/gmd-13-6481-2020, https://doi.org/10.5194/gmd-13-6481-2020, 2020
Short summary
Short summary
Buttressing by floating ice shelves at ice-sheet grounding lines is an
important process that affects ice retreat and whether structural failure
occurs in deep bathymetry. Here, we use a simple algorithm to better
represent 2-D grounding-line curvature in an ice-sheet model. Along with other
enhancements, this improves the performance in idealized-fjord intercomparisons
and enables better diagnosis of potential structural failure at future
retreating Antarctic grounding lines.
Fuyuki Saito, Takashi Obase, and Ayako Abe-Ouchi
Geosci. Model Dev., 13, 5875–5896, https://doi.org/10.5194/gmd-13-5875-2020, https://doi.org/10.5194/gmd-13-5875-2020, 2020
Short summary
Short summary
The present study introduces the rational function-based constrained interpolation profile (RCIP) method for use in 1 d dating computations in ice sheets and demonstrates the performance of the scheme. Comparisons are examined among the RCIP schemes and the first- and second-order upwind schemes. The results show that, in particular, the RCIP scheme preserves the pattern of input histories, in terms of the profile of internal annual layer thickness, better than the other schemes.
Tobias Sauter, Anselm Arndt, and Christoph Schneider
Geosci. Model Dev., 13, 5645–5662, https://doi.org/10.5194/gmd-13-5645-2020, https://doi.org/10.5194/gmd-13-5645-2020, 2020
Short summary
Short summary
Glacial changes play a key role from a socioeconomic, political, and scientific point of view. Here, we present the open-source coupled snowpack and ice surface energy and mass balance model, which provides a lean, flexible, and user-friendly framework for modeling distributed snow and glacier mass changes. The model provides a suitable platform for sensitivity, detection, and attribution analyses for glacier changes and a tool for quantifying inherent uncertainties.
Alex West, Mat Collins, and Ed Blockley
Geosci. Model Dev., 13, 4845–4868, https://doi.org/10.5194/gmd-13-4845-2020, https://doi.org/10.5194/gmd-13-4845-2020, 2020
Short summary
Short summary
This study calculates sea ice energy fluxes from data produced by ice mass balance buoys (devices measuring ice elevation and temperature). It is shown how the resulting dataset can be used to evaluate a coupled climate model (HadGEM2-ES), with biases in the energy fluxes seen to be consistent with biases in the sea ice state and surface radiation. This method has potential to improve sea ice model evaluation, so as to better understand spread in model simulations of sea ice state.
Eduardo Moreno-Chamarro, Pablo Ortega, and François Massonnet
Geosci. Model Dev., 13, 4773–4787, https://doi.org/10.5194/gmd-13-4773-2020, https://doi.org/10.5194/gmd-13-4773-2020, 2020
Short summary
Short summary
Climate models need to capture sea ice complexity to represent it realistically. Here we assess how distributing sea ice in discrete thickness categories impacts how sea ice variability is simulated in the NEMO3.6–LIM3 model. Simulations and satellite observations are compared by using k-means clustering of sea ice concentration in winter and summer between 1979 and 2014 at both poles. Little improvements in the modeled sea ice lead us to recommend using the standard number of five categories.
Ilkka S. O. Matero, Lauren J. Gregoire, and Ruza F. Ivanovic
Geosci. Model Dev., 13, 4555–4577, https://doi.org/10.5194/gmd-13-4555-2020, https://doi.org/10.5194/gmd-13-4555-2020, 2020
Short summary
Short summary
The Northern Hemisphere cooled by several degrees for a century 8000 years ago due to the collapse of an ice sheet in North America that released large amounts of meltwater into the North Atlantic and slowed down its circulation. We numerically model the ice sheet to understand its evolution during this event. Our results match data thanks to good ice dynamics but depend mostly on surface melt and snowfall. Further work will help us understand how past and future ice melt affects climate.
Martin Rückamp, Angelika Humbert, Thomas Kleiner, Mathieu Morlighem, and Helene Seroussi
Geosci. Model Dev., 13, 4491–4501, https://doi.org/10.5194/gmd-13-4491-2020, https://doi.org/10.5194/gmd-13-4491-2020, 2020
Short summary
Short summary
We present enthalpy formulations within the Ice-Sheet and Sea-Level System model that show better performance than earlier implementations. A first experiment indicates that the treatment of discontinuous conductivities of the solid–fluid system with a geometric mean produce accurate results when applied to coarse vertical resolutions. In a second experiment, we propose a novel stabilization formulation that avoids the problem of thin elements. This method provides accurate and stable results.
C. Max Stevens, Vincent Verjans, Jessica M. D. Lundin, Emma C. Kahle, Annika N. Horlings, Brita I. Horlings, and Edwin D. Waddington
Geosci. Model Dev., 13, 4355–4377, https://doi.org/10.5194/gmd-13-4355-2020, https://doi.org/10.5194/gmd-13-4355-2020, 2020
Short summary
Short summary
Understanding processes in snow (firn), including compaction and airflow, is important for calculating how much mass the ice sheets are losing and for interpreting climate records from ice cores. We have developed the open-source Community Firn Model to simulate these processes. We used it to compare 13 different firn compaction equations and found that they do not agree within 10 %. We also show that including firn compaction in a firn-air model improves the match with data from ice cores.
Alexander Robinson, Jorge Alvarez-Solas, Marisa Montoya, Heiko Goelzer, Ralf Greve, and Catherine Ritz
Geosci. Model Dev., 13, 2805–2823, https://doi.org/10.5194/gmd-13-2805-2020, https://doi.org/10.5194/gmd-13-2805-2020, 2020
Short summary
Short summary
Here we describe Yelmo v1.0, an intuitive and state-of-the-art hybrid ice sheet model. The model design and physics are described, and benchmark simulations are provided to validate its performance. Yelmo is a versatile ice sheet model that can be applied to a wide variety of problems.
Ole Richter, David E. Gwyther, Benjamin K. Galton-Fenzi, and Kaitlin A. Naughten
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2020-164, https://doi.org/10.5194/gmd-2020-164, 2020
Revised manuscript accepted for GMD
Short summary
Short summary
Here we present an improved model of the Antarctic continental shelf ocean and demonstrate that it is capable of reproducing present day conditions. The improvements are fundamental and regard the inclusion of tides and ocean eddies. We conclude that the model is well suited to gain new insights into processes that are important for Antarctic ice sheet retreat and global ocean changes. Hence, the model will ultimately help to improve projections of sea level rise and climate change.
Ahmad Jan, Ethan T. Coon, and Scott L. Painter
Geosci. Model Dev., 13, 2259–2276, https://doi.org/10.5194/gmd-13-2259-2020, https://doi.org/10.5194/gmd-13-2259-2020, 2020
Short summary
Short summary
Computer simulations are important tools for understanding the response of Arctic permafrost to a warming climate. To build confidence in an emerging class of permafrost simulators, we evaluated the Advanced Terrestrial Simulator against field observations from a frozen tundra site near Utqiaġvik (formerly Barrow), Alaska. The 3-year simulations agree well with observations of snow depth, summer water table, soil temperature at multiple locations, and spatially averaged evaporation.
Liz C. Logan, Sri Hari Krishna Narayanan, Ralf Greve, and Patrick Heimbach
Geosci. Model Dev., 13, 1845–1864, https://doi.org/10.5194/gmd-13-1845-2020, https://doi.org/10.5194/gmd-13-1845-2020, 2020
Short summary
Short summary
A new capability has been developed for the ice sheet model SICOPOLIS (SImulation COde for POLythermal Ice Sheets) that enables the generation of derivative code, such as tangent linear or adjoint models, by means of algorithmic differentiation. It relies on the source transformation algorithmic (AD) differentiation tool OpenAD. The reverse mode of AD provides the adjoint model, SICOPOLIS-AD, which may be applied for comprehensive sensitivity analyses as well as gradient-based optimization.
Jean-François Lemieux and Frédéric Dupont
Geosci. Model Dev., 13, 1763–1769, https://doi.org/10.5194/gmd-13-1763-2020, https://doi.org/10.5194/gmd-13-1763-2020, 2020
Short summary
Short summary
Sea ice dynamics plays an important role in shaping the sea cover in polar regions. Winds and ocean currents exert large stresses on the sea ice cover. This can lead to the formation of long cracks and ridges, which strongly impact the exchange of heat, momentum and moisture between the atmosphere and the ocean. It is therefore crucial for a sea ice model to be able to represent these features. This article describes how internal sea ice stresses should be diagnosed from model simulations.
Ludovic Räss, Aleksandar Licul, Frédéric Herman, Yury Y. Podladchikov, and Jenny Suckale
Geosci. Model Dev., 13, 955–976, https://doi.org/10.5194/gmd-13-955-2020, https://doi.org/10.5194/gmd-13-955-2020, 2020
Short summary
Short summary
Accurate predictions of future sea level rise require numerical models that predict rapidly deforming ice. Localised ice deformation can be captured numerically only with high temporal and spatial resolution. This paper’s goal is to propose a parallel FastICE solver for modelling ice deformation. Our model is particularly useful for improving our process-based understanding of localised ice deformation. Our solver reaches a parallel efficiency of 99 % on GPU-based supercomputers.
Nander Wever, Leonard Rossmann, Nina Maaß, Katherine C. Leonard, Lars Kaleschke, Marcel Nicolaus, and Michael Lehning
Geosci. Model Dev., 13, 99–119, https://doi.org/10.5194/gmd-13-99-2020, https://doi.org/10.5194/gmd-13-99-2020, 2020
Short summary
Short summary
Sea ice is an important component of the global climate system. The presence of a snow layer covering sea ice can impact ice mass and energy budgets. The detailed, physics-based, multi-layer snow model SNOWPACK was modified to simulate the snow–sea-ice system, providing simulations of the snow microstructure, water percolation and flooding, and superimposed ice formation. The model is applied to in situ measurements from snow and ice mass-balance buoys installed in the Antarctic Weddell Sea.
Christiaan T. van Dalum, Willem Jan van de Berg, Quentin Libois, Ghislain Picard, and Michiel R. van den Broeke
Geosci. Model Dev., 12, 5157–5175, https://doi.org/10.5194/gmd-12-5157-2019, https://doi.org/10.5194/gmd-12-5157-2019, 2019
Short summary
Short summary
Climate models are often limited to relatively simple snow albedo schemes. Therefore, we have developed the SNOWBAL module to couple a climate model with a physically based wavelength dependent snow albedo model. Using SNOWBAL v1.2 to couple the snow albedo model TARTES with the regional climate model RACMO2 indicates a potential performance gain for the Greenland ice sheet.
François Massonnet, Antoine Barthélemy, Koffi Worou, Thierry Fichefet, Martin Vancoppenolle, Clément Rousset, and Eduardo Moreno-Chamarro
Geosci. Model Dev., 12, 3745–3758, https://doi.org/10.5194/gmd-12-3745-2019, https://doi.org/10.5194/gmd-12-3745-2019, 2019
Short summary
Short summary
Sea ice thickness varies considerably on spatial scales of several meters. However, contemporary climate models cannot resolve such scales yet. This is why sea ice models used in climate models include an ice thickness distribution (ITD) to account for this unresolved variability. Here, we explore with the ocean–sea ice model NEMO3.6-LIM3 the sensitivity of simulated mean Arctic and Antarctic sea ice states to the way the ITD is discretized.
Shahbaz Memon, Dorothée Vallot, Thomas Zwinger, Jan Åström, Helmut Neukirchen, Morris Riedel, and Matthias Book
Geosci. Model Dev., 12, 3001–3015, https://doi.org/10.5194/gmd-12-3001-2019, https://doi.org/10.5194/gmd-12-3001-2019, 2019
Short summary
Short summary
Scientific workflows enable complex scientific computational scenarios, which include data intensive scenarios, parametric executions, and interactive simulations. In this article, we applied the UNICORE workflow management system to automate a formerly hard-coded coupling of a glacier flow model and a calving model, which contain many tasks and dependencies, ranging from pre-processing and data management to repetitive executions on heterogeneous high-performance computing (HPC) resources.
Sébastien Le clec'h, Aurélien Quiquet, Sylvie Charbit, Christophe Dumas, Masa Kageyama, and Catherine Ritz
Geosci. Model Dev., 12, 2481–2499, https://doi.org/10.5194/gmd-12-2481-2019, https://doi.org/10.5194/gmd-12-2481-2019, 2019
Short summary
Short summary
To provide reliable projections of the ice-sheet contribution to future sea-level rise, ice sheet models must be able to simulate the observed ice sheet present-day state. Using a low computational iterative minimisation procedure, based on the adjustment of the basal drag coefficient, we rapidly minimise the errors between the simulated and the observed Greenland ice thickness and ice velocity, and we succeed in stabilising the simulated Greenland ice sheet state under present-day conditions.
Lionel Favier, Nicolas C. Jourdain, Adrian Jenkins, Nacho Merino, Gaël Durand, Olivier Gagliardini, Fabien Gillet-Chaulet, and Pierre Mathiot
Geosci. Model Dev., 12, 2255–2283, https://doi.org/10.5194/gmd-12-2255-2019, https://doi.org/10.5194/gmd-12-2255-2019, 2019
Short summary
Short summary
The melting at the base of floating ice shelves is the main driver of the Antarctic ice sheet current retreat. Here, we use an ideal set-up to assess a wide range of melting parameterisations depending on oceanic properties with regard to a new ocean–ice-sheet coupled model, published here for the first time. A parameterisation that depends quadratically on thermal forcing in both a local and a non-local way yields the best results and needs to be further assessed with more realistic set-ups.
Hakase Hayashida, James R. Christian, Amber M. Holdsworth, Xianmin Hu, Adam H. Monahan, Eric Mortenson, Paul G. Myers, Olivier G. J. Riche, Tessa Sou, and Nadja S. Steiner
Geosci. Model Dev., 12, 1965–1990, https://doi.org/10.5194/gmd-12-1965-2019, https://doi.org/10.5194/gmd-12-1965-2019, 2019
Short summary
Short summary
Ice algae, the primary producer in sea ice, play a fundamental role in shaping marine ecosystems and biogeochemical cycling of key elements in polar regions. In this study, we developed a process-based numerical model component representing sea-ice biogeochemistry for a sea ice–ocean coupled general circulation model. The model developed can be used to simulate the projected changes in sea-ice ecosystems and biogeochemistry in response to on-going rapid decline of the Arctic.
Katherine J. Evans, Joseph H. Kennedy, Dan Lu, Mary M. Forrester, Stephen Price, Jeremy Fyke, Andrew R. Bennett, Matthew J. Hoffman, Irina Tezaur, Charles S. Zender, and Miren Vizcaíno
Geosci. Model Dev., 12, 1067–1086, https://doi.org/10.5194/gmd-12-1067-2019, https://doi.org/10.5194/gmd-12-1067-2019, 2019
Short summary
Short summary
A robust validation of ice sheet models is presented using LIVVkit, version 2.1. It targets ice sheet and coupled Earth system models, and handles datasets and operations that require high-performance computing and storage. We apply LIVVkit to a Greenland ice sheet simulation to show the degree to which it captures the surface mass balance. LIVVkit identifies a positive bias due to insufficient melting compared to observations that is focused largely around Greenland's southwest region.
Jeremy C. Ely, Chris D. Clark, David Small, and Richard C. A. Hindmarsh
Geosci. Model Dev., 12, 933–953, https://doi.org/10.5194/gmd-12-933-2019, https://doi.org/10.5194/gmd-12-933-2019, 2019
Short summary
Short summary
During the last 2.6 million years, the Earth's climate has cycled between cold glacials and warm interglacials, causing the growth and retreat of ice sheets. These ice sheets can be independently reconstructed using numerical models or from dated evidence that they leave behind (e.g. sediments, boulders). Here, we present a tool for comparing numerical model simulations with dated ice-sheet material. We demonstrate the utility of this tool by applying it to the last British–Irish ice sheet.
Fabien Maussion, Anton Butenko, Nicolas Champollion, Matthias Dusch, Julia Eis, Kévin Fourteau, Philipp Gregor, Alexander H. Jarosch, Johannes Landmann, Felix Oesterle, Beatriz Recinos, Timo Rothenpieler, Anouk Vlug, Christian T. Wild, and Ben Marzeion
Geosci. Model Dev., 12, 909–931, https://doi.org/10.5194/gmd-12-909-2019, https://doi.org/10.5194/gmd-12-909-2019, 2019
Short summary
Short summary
Mountain glaciers are one of the few remaining subsystems of the global climate system for which no globally applicable community-driven model exists. Here we present the Open Global Glacier Model (OGGM; www.oggm.org), developed to provide a modular and open-source numerical model framework for simulating past and future change of any glacier in the world.
William H. Lipscomb, Stephen F. Price, Matthew J. Hoffman, Gunter R. Leguy, Andrew R. Bennett, Sarah L. Bradley, Katherine J. Evans, Jeremy G. Fyke, Joseph H. Kennedy, Mauro Perego, Douglas M. Ranken, William J. Sacks, Andrew G. Salinger, Lauren J. Vargo, and Patrick H. Worley
Geosci. Model Dev., 12, 387–424, https://doi.org/10.5194/gmd-12-387-2019, https://doi.org/10.5194/gmd-12-387-2019, 2019
Short summary
Short summary
This paper describes the Community Ice Sheet Model (CISM) version 2.1. CISM solves equations for ice flow, heat conduction, surface melting, and other processes such as basal sliding and iceberg calving. It can be used for ice-sheet-only simulations or as the ice sheet component of the Community Earth System Model. Model solutions have been verified for standard test problems. CISM can efficiently simulate the whole Greenland ice sheet, with results that are broadly consistent with observations.
Cited articles
Alexanderian, A., Petra, N., Stadler, G., and Ghattas, O.: A-Optimal Design of
Experiments for Infinite-Dimensional Bayesian Linear Inverse Problems with
Regularized ℓ_0-Sparsification, SIAM J. Sci. Comp.,
36, A2122–A2148, https://doi.org/10.1137/130933381, 2014. a
Alnæs, M., Blechta, J., Hake, J., Johansson, A., Kehlet, B., Logg, A.,
Richardson, C., Ring, J., Rognes, M. E., and Wells, G. N.: The FEniCS
Project Version 1.5, Archive of Numerical Software, 3, 9–23, 2015. a
Alnæs, M. S., Logg, A., Ølgaard, K. B., Rognes, M. E., and Wells, G. N.:
Unified Form Language: A Domain-Specific Language for Weak Formulations of
Partial Differential Equations, ACM T. Math. Softw., 40, 1–37,
https://doi.org/10.1145/2566630, 2014. a
Arthern, R. J., Hindmarsh, R. C. A., and Williams, C. R.: Flow speed within the
Antarctic ice sheet and its controls inferred from satellite observations,
J. Geophys. Res.-Earth, 120, 1171–1188,
https://doi.org/10.1002/2014JF003239, 2015. a
Babaniyi, O., Nicholson, R., Villa, U., and Petra, N.: Inferring the basal sliding coefficient field for the Stokes ice sheet model under rheological uncertainty, The Cryosphere, 15, 1731–1750, https://doi.org/10.5194/tc-15-1731-2021, 2021. a
Bui-Thanh, T., Ghattas, O., Martin, J., and Stadler, G.: A
computational framework for infinite-dimensional Bayesian inverse problems.
Part I: The linearized case, with application to global seismic inversion,
arXiv e-prints, arXiv:1308.1313 2013. a, b, c, d
Chen, P.: Hessian Matrix vs. Gauss-Newton Hessian Matrix, SIAM J. Numer. Anal.,
49, 1417–1435, https://doi.org/10.1137/100799988, 2011. a
Cornford, S. L., Martin, D. F., Graves, D. T., Ranken, D. F., Le Brocq, A. M.,
Gladstone, R. M., Payne, A. J., Ng, E. G., and Lipscomb, W. H.: Adaptive
Mesh, Finite Volume Modeling of Marine Ice Sheets, J. Comput. Phys., 232,
529–549, https://doi.org/10.1016/j.jcp.2012.08.037, 2013. a
Cornford, S. L., Martin, D. F., Payne, A. J., Ng, E. G., Le Brocq, A. M., Gladstone, R. M., Edwards, T. L., Shannon, S. R., Agosta, C., van den Broeke, M. R., Hellmer, H. H., Krinner, G., Ligtenberg, S. R. M., Timmermann, R., and Vaughan, D. G.: Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate, The Cryosphere, 9, 1579–1600, https://doi.org/10.5194/tc-9-1579-2015, 2015. a, b, c, d
Cornford, S. L., Seroussi, H., Asay-Davis, X. S., Gudmundsson, G. H., Arthern, R., Borstad, C., Christmann, J., Dias dos Santos, T., Feldmann, J., Goldberg, D., Hoffman, M. J., Humbert, A., Kleiner, T., Leguy, G., Lipscomb, W. H., Merino, N., Durand, G., Morlighem, M., Pollard, D., Rückamp, M., Williams, C. R., and Yu, H.: Results of the third Marine Ice Sheet Model Intercomparison Project (MISMIP+), The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, 2020. a
Cuffey, K. and Paterson, W. S. B.: The Physics of Glaciers, Butterworth
Heinemann, Oxford, 4th Edn., 2010. a
Deconto, R. M. and Pollard, D.: Contribution of Antarctica to past and
future sea-level rise, Nature, 531, 591–597, https://doi.org/10.1038/nature17145,
2016. a
Dukowicz, J. K., Price, S. F., and Lipscomp, W. H.: Consistent approximations
and boundary conditions for ice-sheet dynamics from a principle of least
action, J. Glaciol., 56, 480–496, 2010. a
Fürst, J. J., Durand, G., Gillet-Chaulet, F., Merino, N., Tavard, L., Mouginot, J., Gourmelen, N., and Gagliardini, O.: Assimilation of Antarctic velocity observations provides evidence for uncharted pinning points, The Cryosphere, 9, 1427–1443, https://doi.org/10.5194/tc-9-1427-2015, 2015. a
Gagliardini, O., Durand, G., Zwinger, T., Hindmarsh, R. C. A., and Meur, E. L.:
Coupling of ice shelf melting and buttressing is a key process in ice sheet
dynamics, Geophys. Res. Lett., 37, L14501, https://doi.org/10.1029/2010GL043334, 2010. a, b
Gillet-Chaulet, F., Gagliardini, O., Seddik, H., Nodet, M., Durand, G., Ritz, C., Zwinger, T., Greve, R., and Vaughan, D. G.: Greenland ice sheet contribution to sea-level rise from a new-generation ice-sheet model, The Cryosphere, 6, 1561–1576, https://doi.org/10.5194/tc-6-1561-2012, 2012. a
Gladstone, R. M., Lee, V., Rougier, J., Payne, A. J., Hellmer, H.,
Le Brocq, A., Shepherd, A., Edwards, T. L., Gregory, J., and
Cornford, S. L.: Calibrated prediction of Pine Island Glacier retreat
during the 21st and 22nd centuries with a coupled flowline model, Earth
Planet. Sc. Lett., 333, 191–199, https://doi.org/10.1016/j.epsl.2012.04.022,
2012. a
Glen, J. W.: The creep of polycrystalline ice, P. Roy. Soc. Lond. A Mat., 228, 519–538, 1955. a
Goldberg, D. N., Heimbach, P., Joughin, I., and Smith, B.: Committed retreat of Smith, Pope, and Kohler Glaciers over the next 30 years inferred by transient model calibration, The Cryosphere, 9, 2429–2446, https://doi.org/10.5194/tc-9-2429-2015, 2015. a, b
Habermann, M., Truffer, M., and Maxwell, D.: Changing basal conditions during the speed-up of Jakobshavn Isbræ, Greenland, The Cryosphere, 7, 1679–1692, https://doi.org/10.5194/tc-7-1679-2013, 2013. a
Hernandez, V., Roman, J. E., and Vidal, V.: SLEPc: A scalable and flexible
toolkit for the solution of eigenvalue problems, ACM T. Math. Softw.,
31, 351–362, 2005. a
Higham, N. J.: Functions of Matrices: Theory and Computation, Society for
Industrial and Applied Mathematics, Philadelphia, PA, USA, 2008. a
Hindmarsh, R. C. A. and Payne, A. J.: Time-step limits for stable
solutions of the ice-sheet equation, Ann. Glaciol., 23, 74–85,
https://doi.org/10.1017/S0260305500013288, 1996. a
Isaac, T., Petra, N., Stadler, G., and Ghattas, O.: Scalable and efficient
algorithms for the propagation of uncertainty from data through inference to
prediction for large-scale problems, with application to flow of the
Antarctic ice sheet, J. Comput. Phys., 296, 348–368,
https://doi.org/10.1016/j.jcp.2015.04.047, 2015. a, b, c, d, e, f, g, h, i, j, k, l, m, n
Joughin, I., Smith, B., and Holland, D. M.: Sensitivity of 21st Century Sea
Level to Ocean-Induced Thinning of Pine Island Glacier, Antarctica,
Geophys. Res. Lett., 37, L20502, https://doi.org/10.1029/2010GL044819, 2010. a, b
Kalmikov, A. G. and Heimbach, P.: A Hessian-Based Method for Uncertainty
Quantification in Global Ocean State Estimation, SIAM J. Sci.
Comp., 36, S267–S295, https://doi.org/10.1137/130925311, 2014. a, b, c, d
Kaminski, T., Kauker, F., Eicken, H., and Karcher, M.: Exploring the utility of quantitative network design in evaluating Arctic sea ice thickness sampling strategies, The Cryosphere, 9, 1721–1733, https://doi.org/10.5194/tc-9-1721-2015, 2015. a, b
Keuthen, M. and Ulbrich, M.: Moreau–Yosida regularization in shape
optimization with geometric constraints, Comput. Optim.
Appl., 62, 181–216, 2015. a
Khodabakhshi, P., Willcox, K. E., and Gunzburger, M.: A multifidelity method
for a nonlocal diffusion model, Appl. Math. Lett., 121, 107361, https://doi.org/10.1016/j.aml.2021.107361,
2021. a
Larour, E., Rignot, E., Joughin, I., and Aubry, D.: Rheology of the Ronne Ice
Shelf, Antarctica, inferred from satellite radar interferometry data using
an inverse control method, Geophys. Res. Lett., 32, L05503,
https://doi.org/10.1029/2004GL021693, 2005. a
Larour, E., Utke, J., Csatho, B., Schenk, A., Seroussi, H., Morlighem, M., Rignot, E., Schlegel, N., and Khazendar, A.: Inferred basal friction and surface mass balance of the Northeast Greenland Ice Stream using data assimilation of ICESat (Ice Cloud and land Elevation Satellite) surface altimetry and ISSM (Ice Sheet System Model), The Cryosphere, 8, 2335–2351, https://doi.org/10.5194/tc-8-2335-2014, 2014. a
Logg, A., Mardal, K.-A., and Wells, G.: Automated Solution of Differential
Equations by the Finite Element Method: The FEniCS Book, Springer Publishing
Company, Incorporated, 2012. a
Loose, N., Heimbach, P., Pillar, H., and Nisancioglu, K.: Quantifying
Dynamical Proxy Potential through Oceanic Teleconnections in the
North Atlantic, Earth and Space Science Open Archive,
Earth and Space Science Open
Archive, https://doi.org/10.1002/essoar.10502065.1, 2020. a
MacAyeal, D. R.: Large-scale ice flow over a viscous basal sediment: Theory and
application to Ice Stream B, Antarctica, J. Geophys.
Res.-Sol. Ea., 94, 4071–4087, 1989. a
MacAyeal, D. R.: The basal stress distribution of Ice Stream E, Antarctica,
inferred by control methods, J. Geophys. Res., 97, 595–603,
1992. a
Maddison, J. R., Goldberg, D. N., and Goddard, B. D.: Automated Calculation of
Higher Order Partial Differential Equation Constrained Derivative
Information, SIAM J. Sci. Comp., 41, C417–C445,
https://doi.org/10.1137/18M1209465, 2019. a, b, c
Martin, J., Wilcox, L. C., Burstedde, C., and Ghattas, O.: A Stochastic Newton
MCMC Method for Large-Scale Statistical Inverse Problems with Application to
Seismic Inversion, SIAM J. Sci. Comp., 34, A1460–A1487,
https://doi.org/10.1137/110845598, 2012. a, b
Morales, J. L. and Nocedal, J.: Remark on “Algorithm 778: L-BFGS-B: Fortran
Subroutines for Large-Scale Bound Constrained Optimization”, ACM T.
Math. Softw., 38, 1–4, https://doi.org/10.1145/2049662.2049669, 2011. a
Morlighem, M., Rignot, E., Seroussi, G., Larour, E., Ben Dhia, H., and Aubry,
D.: Spatial patterns of basal drag inferred using control methods from a
full-Stokes and simpler models for Pine Island Glacier, West Antarctica,
Geophys. Res. Lett., 37, L14502, https://doi.org/10.1029/2010GL043853, 2010. a, b, c
Nias, I. J., Cornford, S. L., and Paybe, A. J.: Contrasting the modelled
sensitivity of the Amundsen Sea Embayment ice streams, J. Glaciol.,
62, 552–562, https://doi.org/10.1017/jog.2016.40, 2016. a
Pattyn, F., Perichon, L., Aschwanden, A., Breuer, B., de Smedt, B., Gagliardini, O., Gudmundsson, G. H., Hindmarsh, R. C. A., Hubbard, A., Johnson, J. V., Kleiner, T., Konovalov, Y., Martin, C., Payne, A. J., Pollard, D., Price, S., Rückamp, M., Saito, F., Souček, O., Sugiyama, S., and Zwinger, T.: Benchmark experiments for higher-order and full-Stokes ice sheet models (ISMIP–HOM), The Cryosphere, 2, 95–108, https://doi.org/10.5194/tc-2-95-2008, 2008. a, b
Petra, N., Martin, J., Stadler, G., and Ghattas, O.: A Computational Framework
for Infinite-Dimensional Bayesian Inverse Problems, Part II: Stochastic
Newton MCMC with Application to Ice Sheet Flow Inverse Problems, SIAM J. Sci. Comp., 36, A1525–A1555, https://doi.org/10.1137/130934805, 2014. a, b, c
Raymond, M. J. and Gudmundsson, G. H.: Estimating basal properties of ice streams from surface measurements: a non-linear Bayesian inverse approach applied to synthetic data, The Cryosphere, 3, 265–278, https://doi.org/10.5194/tc-3-265-2009, 2009. a
Rignot, E., Mouginot, J., and Scheuchl, B.: MEaSUREs InSAR-Based Antarctica
Ice Velocity Map, Version 2, Boulder, Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center [data set],
https://doi.org/10.5067/D7GK8F5J8M8R,
2017. a
Ritz, C., Edwards, T. L., Durand, G., Payne, A. J., Peyaud, V., and
Hindmarsh, R. C. A.: Potential sea-level rise from Antarctic ice-sheet
instability constrained by observations, Nature, 528, 115–118,
https://doi.org/10.1038/nature16147, 2015. a
Robel, A. A., Seroussi, H., and Roe, G. H.: Marine ice sheet instability
amplifies and skews uncertainty in projections of future sea-level rise,
P. Natl. Acad. Sci. USA, 116, 14887–14892,
https://doi.org/10.1073/pnas.1904822116, 2019. a
Rommelaere, V.: Large-scale rheology of the Ross Ice Shelf, Antarctica,
computed by a control method, J. Glaciol., 24, 694–712, 1997. a
Schoof, C.: A variational approach to ice stream flow, J. Fluid Mech., 556,
227–251, 2006. a
Sergienko, O. V., MacAyeal, D. R., and Thom, J. E.: Reconstruction of snow/firn
thermal diffusivities from observed temperature variation: Application to
iceberg C16, Ross Sea, Antarctica, 2004-07, Ann. Glaciol., 49, 91–95,
2008. a
Shapero, D. R., Badgeley, J. A., Hoffman, A. O., and Joughin, I. R.: icepack: a new glacier flow modeling package in Python, version 1.0, Geosci. Model Dev., 14, 4593–4616, https://doi.org/10.5194/gmd-14-4593-2021, 2021. a
Thacker, W. C.: The role of the Hessian matrix in fitting models to
measurements, J. Geophys. Res., 94, 6177–6196,
https://doi.org/10.1029/JC094iC05p06177, 1989. a
Tierney, L.: Markov Chains for Exploring Posterior Distributions, Ann.
Stat., 22, 1701–1728, https://doi.org/10.1214/aos/1176325750, 1994. a
Todd, J. A., Koziol, C. P., Goldberg, D. N., and Maddison, J. R.: EdiGlacUQ/fenics_ice: fenics_ice (v1.0.1), Zenodo [code], https://doi.org/10.5281/zenodo.5153231, 2021. a
Vieli, A. and Payne, A. J.: Application of control methods for modelling the
flow of Pine Island Glacier, West Antarctica, Ann. Glaciol., 36,
197–204, 2003. a
Villa, U., Petra, N., and Ghattas, O.: hIPPYlib: an Extensible Software
Framework for Large-scale Deterministic and Bayesian Inverse Problems,
J. Open Source Softw., 3, p. 940, https://doi.org/10.21105/joss.00940, 2018. a
Zhu, C., Byrd, R. H., Lu, P., and Nocedal, J.: Algorithm 778: L-BFGS-B: Fortran
Subroutines for Large-Scale Bound-Constrained Optimization, ACM T. Math.
Softw., 23, 550–560, https://doi.org/10.1145/279232.279236, 1997. a
Short summary
Sea level change due to the loss of ice sheets presents great risk for coastal communities. Models are used to forecast ice loss, but their evolution depends strongly on properties which are hidden from observation and must be inferred from satellite observations. Common methods for doing so do not allow for quantification of the uncertainty inherent or how it will affect forecasts. We provide a framework for quantifying how this
initialization uncertaintyaffects ice loss forecasts.
Sea level change due to the loss of ice sheets presents great risk for coastal communities....
