Future Antarctic bed topography and its implications for ice sheet dynamics

The Antarctic bedrock is evolving as the solid Earth responds to the past and ongoing evolution of the ice sheet. A recently improved ice loading history suggests that the Antarctic Ice Sheet (AIS) has generally been losing its mass since the Last Glacial Maximum. In a sustained warming climate, the AIS is predicted to retreat at a greater pace, primarily via melting beneath the ice shelves. We employ the glacial isostatic adjustment (GIA) capability of the Ice Sheet System Model (ISSM) to combine these past and future ice loadings and provide the new solid Earth computations for the AIS. We find that past loading is relatively less important than future loading for the evolution of the future bed topography. Our computations predict that the West Antarctic Ice Sheet (WAIS) may uplift by a few meters and a few tens of meters at years AD 2100 and 2500, respectively, and that the East Antarctic Ice Sheet is likely to remain unchanged or subside minimally except around the Amery Ice Shelf. The Amundsen Sea Sector in particular is predicted to rise at the greatest rate; one hundred years of ice evolution in this region, for example, predicts that the coastline of Pine Island Bay will approach roughly 45 mm yr−1 in viscoelastic vertical motion. Of particular importance, we systematically demonstrate that the effect of a pervasive and large GIA uplift in the WAIS is generally associated with the flattening of reverse bed slope, reduction of local sea depth, and thus the extension of grounding line (GL) towards the continental shelf. Using the 3-D higher-order ice flow capability of ISSM, such a migration of GL is shown to inhibit the ice flow. This negative feedback between the ice sheet and the solid Earth may promote stability in marine portions of the ice sheet in the future

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

Polarization analysis of CuXX-lines emitted from X-pinch

International audienceSoft x-ray emission from CuXX L-shell lines emitted by a dense X-pinch plasma have been investigated with high-resolution curved Bragg crystals at different angles of orientation. Single shot time integrated spectra show clear evidences of polarization for the Ne-like spectral lines 2s22p6 1S0 → 2s22p53s 1P1 (λ = 12.570 Å), 2s22p6 1S0 → 2s22p53s 3P1 (λ = 12.8277 Å). The variation of the intensity ratio of these two well-separated L-shell lines is discussed in view of its application for suprathermal electron characterization under real experimental conditions of pinch plasmas. We demonstrated that the simultaneous use of two different polarization spectrometers (means 4 Bragg crystals) permitted a high level of confidence for the analysis of the variation of the line ratios due to polarization

Representation of sharp rifts and faults mechanics in modeling ice shelf flow dynamics: Application to Brunt/Stancomb-Wills Ice Shelf, Antarctica

Ice shelves play a major role in buttressing ice sheet flow into the ocean, hence the importance of accurate numerical modeling of their stress regime. Commonly used ice flow models assume a continuous medium and are therefore complicated by the presence of rupture features (crevasses, rifts, and faults) that significantly affect the overall flow patterns. Here we apply contact mechanics and penalty methods to develop a new ice shelf flow model that captures the impact of rifts and faults on the rheology and stress distribution of ice shelves. The model achieves a best fit solution to satellite observations of ice shelf velocities to infer the following: (1) a spatial distribution of contact and friction points along detected faults and rifts, (2) a more realistic spatial pattern of ice shelf rheology, and (3) a better representation of the stress balance in the immediate vicinity of faults and rifts. Thus, applying the model to the Brunt/Stancomb-Wills Ice Shelf, Antarctica, we quantify the state of friction inside faults and the opening rates of rifts and obtain an ice shelf rheology that remains relatively constant everywhere else on the ice shelf. We further demonstrate that better stress representation has widespread application in examining aspects affecting ice shelf structure and dynamics including the extent of ice mélange in rifts and the change in fracture configurations. All are major applications for better insight into the important question of ice shelf stability

High-resolution ice-thickness mapping in South Greenland

Airborne radar sounding is difficult in South Greenland because of the presence of englacial water, which prevents the signal from reaching the bed. Data coverage remains suboptimal for traditional methods of ice-thickness and bed mapping that rely on geostatistical techniques, such as kriging, because important features are missing. Here we apply two alternative approaches of highresolution (̃300 m) ice-thickness mapping, that are based on the conservation of mass, to two regions of South Greenland: (1) Qooqqup Sermia and Kiattuut Sermiat, and (2) Ikertivaq. These two algorithms solve optimization problems, for which the conservation of mass is either enforced as a hard constraint, or as a soft constraint. For the first region, very few measurements are available but there is no gap in ice motion data, whereas for Ikertivaq, more ice-thickness measurements are available, but there are gaps in ice motion data. We show that mass-conservation algorithms can be used as validation tools for radar sounding. We also show that it is preferable to apply mass conservation as a hard constraint, rather than a soft constraint, as it better preserves elongated features, such as glacial valleys and ridges

Ice Sheet Model Intercomparison Project (ISMIP6) contribution to CMIP6

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

Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11)

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

Sensitivity of the dynamics of Pine Island Glacier, West Antarctica, to climate forcing for the next 50 years

Pine Island Glacier, a major contributor to sea level rise in West Antarctica, has been undergoing significant changes over the last few decades. Here, we employ a three-dimensional, higher-order model to simulate its evolution over the next 50 yr in response to changes in its surface mass balance, the position of its calving front and ocean-induced ice shelf melting. Simulations show that the largest climatic impact on ice dynamics is the rate of ice shelf melting, which rapidly affects the glacier speed over several hundreds of kilometers upstream of the grounding line. Our simulations show that the speedup observed in the 1990s and 2000s is consistent with an increase in sub-ice-shelf melting. According to our modeling results, even if the grounding line stabilizes for a few decades, we find that the glacier reaction can continue for several decades longer. Furthermore, Pine Island Glacier will continue to change rapidly over the coming decades and remain a major contributor to sea level rise, even if ocean-induced melting is reduced

Simulating ice thickness and velocity evolution of Upernavik Isstrom 1849-2012 by forcing prescribed terminus positions in ISSM

Abstract. Tidewater glacier velocity and mass balance are known to be highly responsive to terminus position change. Yet it remains challenging for ice flow models to reproduce observed ice margin changes. Here, using the Ice Sheet System Model (Larour et al., 2012), we simulate the ice velocity and thickness changes of Upernavik Isstrøm (north-western Greenland) by prescribing a collection of 27 observed terminus positions spanning 164 years (1849–2012). The simulation shows increased ice velocity during the 1930s, the late 1970s and between 1995 and 2012 when terminus retreat was observed along with negative surface mass balance anomalies. Three distinct mass balance states are evident in the reconstruction: (1849–1932) with near zero mass balance, (1932–1992) with ice mass loss dominated by ice dynamical flow, and (1998–2012), when increased retreat and negative surface mass balance anomalies led to mass loss that was twice that of any earlier period. Over the multi-decadal simulation, mass loss was dominated by thinning and acceleration responsible for 70 % of the total mass loss induced by prescribed change in terminus position. The remaining 30 % of the total ice mass loss resulted directly from prescribed terminus retreat and decreasing surface mass balance. Although the method can not explain the cause of glacier retreat, it enables the reconstruction of ice flow and geometry during 1849–2012. Given annual or seasonal observed terminus front positions, this method could be a useful tool for evaluating simulations investigating the effect of calving laws. </jats:p