ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods

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

Rapid submarine ice melting in the grounding zones of ice shelves in West Antarctica.

Enhanced submarine ice-shelf melting strongly controls ice loss in the Amundsen Sea embayment (ASE) of West Antarctica, but its magnitude is not well known in the critical grounding zones of the ASE's major glaciers. Here we directly quantify bottom ice losses along tens of kilometres with airborne radar sounding of the Dotson and Crosson ice shelves, which buttress the rapidly changing Smith, Pope and Kohler glaciers. Melting in the grounding zones is found to be much higher than steady-state levels, removing 300-490 m of solid ice between 2002 and 2009 beneath the retreating Smith Glacier. The vigorous, unbalanced melting supports the hypothesis that a significant increase in ocean heat influx into ASE sub-ice-shelf cavities took place in the mid-2000s. The synchronous but diverse evolutions of these glaciers illustrate how combinations of oceanography and topography modulate rapid submarine melting to hasten mass loss and glacier retreat from West Antarctica

Complex Basal Thermal Transition Near the Onset of Petermann Glacier, Greenland

The basal thermal regime of ice sheets exerts a strong control on ice‐sheet stability and the onset of rapidly streaming flow. However, the nature of this thermal transition where sliding initiates is largely unconstrained by geophysical observations. In the Greenland Ice Sheet, topographic troughs or elevated geothermal heat fluxes typically define the onset of outlet glaciers. In contrast, Petermann Glacier in Northern Greenland does not have any distinct bed troughs or localized geothermal heating associated with its onset, making it an ideal site to investigate the basal thermal state and examine its role in the onset of Petermann Glacier. Here we use radar bed reflectivity and an ice‐sheet thermomechanical model to examine the basal thermal regime beneath Petermann Glacier. Our results reveal a complex thermal transition near the onset of Petermann Glacier. As the bed shifts from largely frozen to largely thawed with increasing distances from the ice divide, our results show that this thermal transition happens through alternating bands of frozen and thawed bands. The complex thermal state across the onset region suggests that lateral meltwater injection and local meltwater production determine the location of Petermann Glacier. Given the lack of topographic pinning at the onset location, the upstream margin of Petermann is vulnerable to migrate depending on a combination of advective cooling and meltwater supply from the interior of the ice sheet

Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14)

Accurate projections of the evolution of ice sheets in a changing climate require a fine mesh/grid resolution in ice sheet models to correctly capture fundamental physical processes, such as the evolution of the grounding line, the region where grounded ice starts to float. The evolution of the grounding line indeed plays a major role in ice sheet dynamics, as it is a fundamental control on marine ice sheet stability. Numerical modeling of a grounding line requires significant computational resources since the accuracy of its position depends on grid or mesh resolution. A technique that improves accuracy with reduced computational cost is the adaptive mesh refinement (AMR) approach. We present here the implementation of the AMR technique in the finite element Ice Sheet System Model (ISSM) to simulate grounding line dynamics under two different benchmarks: MISMIP3d and MISMIP+. We test different refinement criteria: (a) distance around the grounding line, (b) a posteriori error estimator, the Zienkiewicz-Zhu (ZZ) error estimator, and (c) different combinations of (a) and (b). In both benchmarks, the ZZ error estimator presents high values around the grounding line. In the MISMIP+ setup, this estimator also presents high values in the grounded part of the ice sheet, following the complex shape of the bedrock geometry. The ZZ estimator helps guide the refinement procedure such that AMR performance is improved. Our results show that computational time with AMR depends on the required accuracy, but in all cases, it is significantly shorter than for uniformly refined meshes. We conclude that AMR without an associated error estimator should be avoided, especially for real glaciers that have a complex bed geometry.121215232CONSELHO NACIONAL DE DESENVOLVIMENTO CIENTÍFICO E TECNOLÓGICO - CNPQ140186/2015-

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

Pathways of ocean heat towards Pine Island and Thwaites grounding lines

In the Amundsen Sea, modified Circumpolar Deep Water (mCDW) intrudes into ice shelf cavities, causing high ice shelf melting near the ice sheet grounding lines, accelerating ice flow, and controlling the pace of future Antarctic contributions to global sea level. The pathways of mCDW towards grounding lines are crucial as they directly control the heat reaching the ice. A realistic representation of mCDW circulation, however, remains challenging due to the sparsity of in-situ observations and the difficulty of ocean models to reproduce the available observations. In this study, we use an unprecedentedly high-resolution (200 m horizontal and 10 m vertical grid spacing) ocean model that resolves shelf-sea and sub-ice-shelf environments in qualitative agreement with existing observations during austral summer conditions. We demonstrate that the waters reaching the Pine Island and Thwaites grounding lines follow specific, topographically-constrained routes, all passing through a relatively small area located around 104°W and 74.3°S. The temporal and spatial variabilities of ice shelf melt rates are dominantly controlled by the sub-ice shelf ocean current. Our findings highlight the importance of accurate and high-resolution ocean bathymetry and subglacial topography for determining mCDW pathways and ice shelf melt rates

An Interdisciplinary Perspective on Greenland’s Changing Coastal Margins

Greenland’s coastal margins are influenced by the confluence of Arctic and Atlantic waters, sea ice, icebergs, and meltwater from the ice sheet. Hundreds of spectacular glacial fjords cut through the coastline and support thriving marine ecosystems and, in some places, adjacent Greenlandic communities. Rising air and ocean temperatures, as well as glacier and sea-ice retreat, are impacting the conditions that support these systems. Projecting how these regions and their communities will evolve requires understanding both the large-scale climate variability and the regional-scale web of physical, biological, and social interactions. Here, we describe pan-Greenland physical, biological, and social settings and show how they are shaped by the ocean, the atmosphere, and the ice sheet. Next, we focus on two communities, Qaanaaq in Northwest Greenland, exposed to Arctic variability, and Ammassalik in Southeast Greenland, exposed to Atlantic variability. We show that while their climates today are similar to those of the warm 1930s­–1940s, temperatures are projected to soon exceed those of the last 100 years at both locations. Existing biological records, including fisheries, provide some insight on ecosystem variability, but they are too short to discern robust patterns. To determine how these systems will evolve in the future requires an improved understanding of the linkages and external factors shaping the ecosystem and community response. This interdisciplinary study exemplifies a first step in a systems approach to investigating the evolution of Greenland’s coastal margins

Multisystem synthesis of radar sounding observations of the Amundsen Sea sector from the 2004‐2005 field season

The Amundsen Sea Embayment of the West Antarctic Ice Sheet contains Thwaites and Pine Island Glaciers, two of the most rapidly changing glaciers in Antarctica. To date, Pine Island and Thwaites Glaciers have only been observed by independent airborne radar sounding surveys, but a combined cross-basin analysis that investigates the basal conditions across the Pine Island-Thwaites Glaciers boundary has not been performed. Here, we combine two radar surveys and correct for their differences in system parameters to produce unified englacial attenuation and basal relative reflectivity maps spanning both Pine Island and Thwaites Glaciers. Relative reflectivities range from -24.8 to +37.4 dB with the highest values beneath fast-flowing ice at the ice sheet margin. By comparing our reflectivity results with previously derived radar specularity and trailing bed echoes at Thwaites Glacier, we find a highly diverse subglacial landscape and hydrologic condition that evolve along flow. Together, these findings highlight the potential for joint airborne radar analysis with ground-based seismic and geomorphological observations to understand variations in the bed properties and cross-catchment interactions of ice streams and outlet glaciers