Search

This paper provides the first description of version 1.0 of the open-source Glacier Energy and Mass Balance model. GEMB models the ice sheet and glacier surface-atmospheric energy and mass exchange, and firn state. It is a column model (no horizontal communication) of intermediate complexity that includes those processes deemed most relevant to glacier studies with the goal of retaining computational efficiency that can accommodate the very long (thousands of years) spin-ups necessary for initializing deep firn columns and for running sensitivity experiments to characterize model uncertainty on continental scales. The model is one-way coupled with the atmosphere which allows the model to be run off-line with a diversity of climate forcing but neglects feedback to the atmosphere. GEMB provides numerous parameterization choices for various key processes (e.g. albedo, subsurface shortwave absorption, and compaction), making it well suited for uncertainty quantification and model exploration. The model is evaluated against the current state-of-the-art and in situ observations and is shown to perform well. [ABSTRACT FROM AUTHOR]

Geologic archives constraining the variability of the Greenland ice sheet (GrIS) during the Holocene provide targets for ice sheet models to test sensitivities to variations in past climate and model formulation. Even as data–model comparisons are becoming more common, many models simulating the behavior of the GrIS during the past rely on meshes with coarse horizontal resolutions (≥10 km). In this study, we explore the impact of model resolution on the simulated nature of retreat across southwestern Greenland during the Holocene. Four simulations are performed using the Ice Sheet System Model (ISSM): three that use a uniform mesh and horizontal mesh resolutions of 20, 10, and 5 km, and one that uses a nonuniform mesh with a resolution ranging from 2 to 15 km. We find that the simulated retreat can vary significantly between models with different horizontal resolutions based on how well the bed topography is resolved. In areas of low topographic relief, the horizontal resolution plays a negligible role in simulated differences in retreat, with each model instead responding similarly to retreat driven by surface mass balance (SMB). Conversely, in areas where the bed topography is complex and high in relief, such as fjords, the lower-resolution models (10 and 20 km) simulate unrealistic retreat that occurs as ice surface lowering intersects bumps in the bed topography that would otherwise be resolved as troughs using the higher-resolution grids. Our results highlight the important role that high-resolution grids play in simulating retreat in areas of complex bed topography, but also suggest that models using nonuniform grids can save computational resources through coarsening the mesh in areas of noncomplex bed topography where the SMB predominantly drives retreat. Additionally, these results emphasize that care must be taken with ice sheet models when tuning model parameters to match reconstructed margins, particularly for lower-resolution models in regions where complex bed topography is poorly resolved. [ABSTRACT FROM AUTHOR]

Identifying fast and robust numerical solvers is a critical issue that needs to be addressed in order to improve projections of polar ice sheets evolving in a changing climate. This work evaluates the impact of using advanced numerical solvers for transient ice-flow simulations conducted with the JPL–UCI Ice Sheet System Model (ISSM). We identify optimal numerical solvers by testing a broad suite of readily available solvers, ranging from direct sparse solvers to preconditioned iterative methods, on the commonly used Ice Sheet Model Intercomparison Project for Higher-Order ice sheet Models benchmark tests. Three types of analyses are considered: mass transport, horizontal stress balance, and incompressibility. The results of the fastest solvers for each analysis type are ranked based on their scalability across mesh size and basal boundary conditions. We find that the fastest iterative solvers are  ∼ 1.5–100 times faster than the default direct solver used in ISSM, with speed-ups improving rapidly with increased mesh resolution. We provide a set of recommendations for users in search of efficient solvers to use for transient ice-flow simulations, enabling higher-resolution meshes and faster turnaround time. The end result will be improved transient simulations for short-term, highly resolved forward projections (10–100 year time scale) and also improved long-term paleo-reconstructions using higher-order representations of stresses in the ice. This analysis will also enable a new generation of comprehensive uncertainty quantification assessments of forward sea-level rise projections, which rely heavily on ensemble or sampling approaches that are inherently expensive. [ABSTRACT FROM AUTHOR]

Reducing the uncertainty in the past, present, and future contribution of ice sheets to sea-level change requires a coordinated effort between the climate and glaciology communities. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the Greenland and Antarctic ice sheets. In this paper, we describe the framework for ISMIP6 and its relationship with other activities within CMIP6. The ISMIP6 experimental design relies on CMIP6 climate models and includes, for the first time within CMIP, coupled ice-sheet-climate models as well as standalone ice-sheet models. To facilitate analysis of the multi-model ensemble and to generate a set of standard climate inputs for standalone ice-sheet models, ISMIP6 defines a protocol for all variables related to ice sheets. ISMIP6 will provide a basis for investigating the feedbacks, impacts, and sea-level changes associated with dynamic ice sheets and for quantifying the uncertainty in ice-sheet-sourced global sea-level change. [ABSTRACT FROM AUTHOR]

Identifying fast and robust numerical solvers is a critical issue that needs to be addressed in order to improve projections of polar ice sheets evolving in a changing climate. This work evaluates the impact of using sophisticated numerical solvers for transient ice flow simulations using the NASA-JPL/UCI Ice Sheet System Model (ISSM). We identify optimal numerical solvers by testing them on a commonly used ice flow benchmark test, the Ice Sheet Model Intercomparison Project for Higher-Order ice sheet Models (ISMIP-HOM) Experiment F. Three types of analyses are considered: mass transport, horizontal stress balance, and vertical stress balance. A broad suite of solvers is tested, ranging from direct sparse solvers to preconditioned iterative methods. The results of the fastest solvers for each analysis type are ranked based on their scalability across mesh size for each basal sliding conditions specified in Experiment F. We find that the fastest iterative solvers are ~1.5-100 times faster than the default direct solver used in ISSM with speed-ups improving rapidly with increased mesh resolution. We provide a set of recommendations for users in search of efficient solvers to use for transient ice flow simulations, enabling higher-resolution meshes and faster turnaround time. The end result will be improved transient simulations for short-term, highly resolved forward projections (10-100 year time scale) and also improved long-term paleo-reconstructions using higher-order representation of stresses in the ice. This analysis will also enable a new generation of comprehensive uncertainty quantification assessments of forward sea-level rise projections, which rely heavily on ensemble or sampling approaches that are inherently expensive. [ABSTRACT FROM AUTHOR]