A JavaScript API for the Ice Sheet System Model: towards on online interactive model for the Cryosphere Community

Abstract. Earth System Models (ESMs) are becoming increasingly complex, requiring extensive knowledge and experience to deploy and use in an efficient manner. They run on high-performance architectures that are significantly different from the everyday environments that scientists use to pre and post-process results (i.e. MATLAB, Python). This results in models that are hard to use for non specialists, and that are increasingly specific in their application. It also makes them relatively inaccessible to the wider science community, not to mention to the general public. Here, we present a new software/model paradigm that attempts to bridge the gap between the science community and the complexity of ESMs, by developing a new JavaScript Application Program Interface (API) for the Ice Sheet System Model (ISSM). The aforementioned API allows Cryosphere Scientists to run ISSM on the client-side of a webpage, within the JavaScript environment. When combined with a Web server running ISSM (using a Python API), it enables the serving of ISSM computations in an easy and straightforward way. The deep integration and similarities between all the APIs in ISSM (MATLAB, Python, and now JavaScript) significantly shortens and simplifies the turnaround of state-of-the-art science runs and their use by the larger community. We demonstrate our approach via a new Virtual Earth System Laboratory (VESL) Web site

Linking glacially modified waters to catchment-scale subglacial discharge using autonomous underwater vehicle observations

© The Author(s), 2016. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Cryosphere 10 (2016): 417-432, doi:10.5194/tc-10-417-2016.Measurements of near-ice (<  200 m) hydrography and near-terminus subglacial hydrology are lacking, due in large part to the difficulty in working at the margin of calving glaciers. Here we pair detailed hydrographic and bathymetric measurements collected with an autonomous underwater vehicle as close as 150 m from the ice–ocean interface of the Saqqarliup sermia–Sarqardleq Fjord system, West Greenland, with modeled and observed subglacial discharge locations and magnitudes. We find evidence of two main types of subsurface glacially modified water (GMW) with distinct properties and locations. The two GMW locations also align with modeled runoff discharged at separate locations along the grounded margin corresponding with two prominent subcatchments beneath Saqqarliup sermia. Thus, near-ice observations and subglacial discharge routing indicate that runoff from this glacier occurs primarily at two discrete locations and gives rise to two distinct glacially modified waters. Furthermore, we show that the location with the largest subglacial discharge is associated with the lighter, fresher glacially modified water mass. This is qualitatively consistent with results from an idealized plume model.Support was provided by the National Science Foundation’s Office of Polar Programs (NSF-OPP) through PLR-1418256 to F. Straneo, S. B. Das and A. J. Plueddemann, PLR-1023364 to S. B. Das, and through the Woods Hole Oceanographic Institution Ocean and Climate Change Institute Arctic Research Initiative to F. Straneo, S. B. Das, and A. J. Plueddemann. L. A. Stevens was also supported by a National Science Foundation Graduate Research Fellowship. S. B. Das was also supported by the Woods Hole Oceanographic Institution James E. and Barbara V. Moltz Research Fellowship. M. Morlighem was supported by the National Aeronautics and Space Administration’s (NASA) Cryospheric Sciences Program through NNX15AD55G

Terminus-driven retreat of a major southwest Greenland tidewater glacier during the early 19th century : insights from glacier reconstructions and numerical modelling

Peer reviewedPublisher PD

Plastic bed beneath Hofsjökull Ice Cap, central Iceland, and the sensitivity of ice flow to surface meltwater flux

The mechanical properties of glacier beds play a fundamental role in regulating the sensitivity of glaciers to environmental forcing across a wide range of timescales. Glaciers are commonly underlain by deformable till whose mechanical properties and influence on ice flow are not well understood but are critical for reliable projections of future glacier states. Using synoptic-scale observations of glacier motion in different seasons to constrain numerical ice flow models, we study the mechanics of the bed beneath Hofsjökull, a land-terminating ice cap in central Iceland. Our results indicate that the bed deforms plastically and weakens following incipient summertime surface melt. Combining the inferred basal shear traction fields with a Coulomb-plastic bed model, we estimate the spatially distributed effective basal water pressure and show that changes in basal water pressure and glacier accelerations are non-local and non-linear. These results motivate an idealized physical model relating mean basal water pressure and basal slip rate wherein the sensitivity of glacier flow to changes in basal water pressure is inversely related to the ice surface slope.This research was conducted at the California Institute of Technology and the University of Iceland with funding provided by the NASA Crysopherice Sciences Program (Award NNX14AH80G). B. M. was partially funded by a NASA Earth and Space Sciences Fellowship and an Achievement Rewards for College Students (ARCS) fellowship. InSAR data are freely available from the Alaska Satellite Facility via the UAVSAR website (http://uavsar.jpl.nasa.gov).Peer Reviewe

Evaluation of four calving laws for Antarctic ice shelves

Many floating ice shelves in Antarctica buttress the ice streams feeding them, thereby reducing the discharge of icebergs into the ocean. The rate at which ice shelves calve icebergs and how fast they flow determine whether they advance, retreat, or remain stable, exerting a first-order control on ice discharge. To parameterize calving within ice sheet models, several empirical and physical calving “laws” have been proposed in the past few decades. Such laws emphasize dissimilar features, including along- and across-flow strain rates (the eigencalving law), a fracture yield criterion (the von Mises law), longitudinal stretching (the crevasse depth law), and a simple ice thickness threshold (the minimum thickness law), among others. Despite the multitude of established calving laws, these laws remain largely unvalidated for the Antarctic Ice Sheet, rendering it difficult to assess the broad applicability of any given law in Antarctica. We address this shortcoming through a set of numerical experiments that evaluate existing calving laws for 10 ice shelves around the Antarctic Ice Sheet. We utilize the Ice-sheet and Sea-level System Model (ISSM) and implement four calving laws under constant external forcing, calibrating the free parameter of each of these calving laws for each ice shelf by assuming that the current position of the ice front is in steady state and finding the set of parameters that best achieves this position over a simulation of 200 years. We find that, in general, the eigencalving and von Mises laws best reproduce observed calving front positions under the steady-state position assumption. These results will streamline future modeling efforts of Antarctic ice shelves by better informing the relevant physics of Antarctic-style calving on a shelf-by-shelf basis.</p

Calibrating calving parameterizations using graph neural network emulators: application to Helheim Glacier, East Greenland

Calving is responsible for the retreat, acceleration, and thinning of numerous tidewater glaciers in Greenland. An accurate representation of this process in ice sheet numerical models is critical to better predict the future response of the ice sheet to climate change. While traditional numerical models have been used to simulate ice dynamics and calving under specific parameterized conditions, the computational demand of these models makes it difficult to efficiently fine-tune these parameterizations, adding to the overall uncertainty in future sea level rise. In this study, we adopt three standard graph neural network (GNN) architectures, including graph convolutional network, graph attention network, and equivariant graph convolutional network (EGCN), to develop surrogate models for finite-element simulations from the Ice-sheet and Sea-level System Model. GNNs are particularly well-suited for this problem as they naturally capture the representation of unstructured meshes used by finite-element models. When these GNNs are trained with numerical simulations of Helheim Glacier, Greenland, for different calving stress thresholds, they successfully reproduce the observed evolution of ice velocity, ice thickness, and ice front migration between 2007 and 2020. Moreover, these emulators exhibit uncertainties of less than 10 %–20 % when extrapolating to out-of-sample calving parameterization cases. Among the three GNN architectures, EGCN outperforms the others by preserving the equivariance of graph structures. By leveraging the GPU-based GNN emulators, which are 30–34 times faster than traditional numerical simulations, we determine the temporal variations of the optimal calving threshold that minimizes the misfit between modeled and observed ice fronts. This fine-tuned calving parameterization, enabled by GNN emulators, can enhance the reliability of numerical models in capturing glacier mass loss driven by calving.</p

Source-transformation Differentiation of a C++-like Ice Sheet model

International audienceAlgorithmic Differentiation (AD) has become one of the most powerful tools to improve our understanding of theEarth System. If AD has been used by the ocean and atmospheric circulation modelingcommunity for almost 20 years, it is relatively new in the ice sheet modeling community. The Ice SheetSystem Model (ISSM) is a C++, object-oriented, massively parallelized, new generation ice sheet model that recentlyimplemented AD to improve its data assimilation capabilities. ISSM currently relies on Object Overloading throughADOL-C and AMPI. However, experience shows that Object Overloading AD on ISSM is significantly more memoryintensive compared to the primal code. We want to investigate other AD approaches to improve the performance ofthe AD adjoint. Yet, to our knowledge, there is no source-to-source AD tool that supports C++.To overcome this problem, we have developed a prototype of ISSM entirely in C, called Boreas, in order to testsource-to-source transformation and compare the performance of these two approaches to AD. Boreas is a clone ofISSM, the main difference with ISSM is that all the objects are converted to C-structures and some function nameshave been adapted in order to be unique, but the code architectures are identical. The programming style of Boreas isa first attempt at defining a programming style of (or a sub-language of) C++ that source-transformation AD couldhandle. We present here the first results of Source-Transformation AD of Boread with the AD tool Tapenade

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 thermal state of an ice sheet is an important control on its past and future evolution. Some parts of the ice sheet may be polythermal, leading to discontinuous properties at the cold–temperate transition surface (CTS). These discontinuities require a careful treatment in ice sheet models (ISMs). Additionally, the highly anisotropic geometry of the 3D elements in ice sheet modelling poses a problem for stabilization approaches in advection-dominated problems. Here, we present extended enthalpy formulations within the finite-element Ice-Sheet and Sea-Level System model (ISSM) that show a better performance than earlier implementations. In a first polythermal-slab experiment, we found that the treatment of the discontinuous conductivities at the CTS with a geometric mean produces more accurate results compared to the arithmetic or harmonic mean. This improvement is particularly efficient when applied to coarse vertical resolutions. In a second ice dome experiment, we find that the numerical solution is sensitive to the choice of stabilization parameters in the well-established streamline upwind Petrov–Galerkin (SUPG) method. As standard literature values for the SUPG stabilization parameter do not account for the highly anisotropic geometry of the 3D elements in ice sheet modelling, we propose a novel anisotropic SUPG (ASUPG) formulation. This formulation circumvents the problem of high aspect ratio by treating the horizontal and vertical directions separately in the stabilization coefficients. The ASUPG method provides accurate results for the thermodynamic equation on geometries with very small aspect ratios like ice sheets

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

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