Mathieu Morlighem

Widespread, rapid grounding line retreat of Pine Island, Thwaites, Smith, and Kohler glaciers, West Antarctica, from 1992 to 2011

We measure the grounding line retreat of glaciers draining the Amundsen Sea sector of West Antarctica using Earth Remote Sensing (ERS-1/2) satellite radar interferometry from 1992 to 2011. Pine Island Glacier retreated 31 km at its center, with most retreat in 2005-2009 when the glacier ungrounded from its ice plain. Thwaites Glacier retreated 14 km along its fast flow core and 1 to 9 km along the sides. Haynes Glacier retreated 10 km along its flanks. Smith/Kohler glaciers retreated the most, 35 km along its ice plain, and its ice shelf pinning points are vanishing. These rapid retreats proceed along regions of retrograde bed elevation mapped at a high spatial resolution using a mass conservation technique that removes residual ambiguities from prior mappings. Upstream of the 2011 grounding line positions, we find no major bed obstacle that would prevent the glaciers from further retreat and draw down the entire basin.

Spatial patterns of basal drag inferred using control methods from a full-Stokes and simpler models for Pine Island Glacier, West Antarctica

[1] Basal drag is a fundamental control on ice stream dynamics that remains poorly understood or constrained by observations. Here, we apply control methods on ice surfacevelocities ofPine IslandGlacier, WestAntarctica toinfer the spatial pattern of basal drag using a full‐Stokes (FS) model of ice flow and compare the results obtained with two commonly‐used simplified solutions: the MacAyeal shelfy stream model and the Blatter‐Pattyn model. Over most of the model domain, the three models yield similar patterns of basal drag, yet near the glacier grounding‐line, the simplified models yield high basal drag while FS yields almost no basal drag. The simplified models overestimate basal drag because they neglect bridging effects in an ice stream region of rapidly varying ice thickness. This result reinforces theoretical studies that a FS treatment of ice flow is essential near glacier grounding lines. Citation: Morlighem, M., E. Rignot, H. Seroussi, E. Larour, H. Ben Dhia, and D. Aubry (2010), 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, doi:10.1029/2010GL043853.

Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM)

Ice flow models used to project the mass balance of ice sheets in Greenland and Antarctica usually rely on the Shallow Ice Approximation (SIA) and the Shallow-Shelf Approximation (SSA), sometimes combined into so-called hybrid models. Such models, while computationally efficient, are based on a simplified set of physical assumptions about the mechanical regime of the ice flow, which does not uniformly apply everywhere on the ice sheet/ice shelf system, especially near grounding lines, where rapid changes are taking place at present. Here, we present a new thermomechanical finite element model of ice flow named ISSM (Ice Sheet System Model) that includes higher-order stresses, high spatial resolution capability and data assimilation techniques to better capture ice dynamics and produce realistic simulations of ice sheet flow at the continental scale. ISSM includes several approximations of the momentum balance equations, ranging from the two-dimensional SSA to the three-dimensional full-Stokes formulation. It also relies on a massively parallelized architecture and state-of-the-art scalable tools. ISSM employs data assimilation techniques, at all levels of approximation of the momentum balance equations, to infer basal drag at the ice-bed interface from satellite radar interferometry-derived observations of ice motion. Following a validation of ISSM with standard benchmarks, we present a demonstration of its capability in the case of the Greenland Ice Sheet. We show ISSM is able to simulate the ice flow of an entire ice sheet realistically at a high spatial resolution, with higher-order physics, thereby providing a pathway for improving projections of ice sheet evolution in a warming climate. Copyright 2012 by the American Geophysical Union.

A mass conservation approach for mapping glacier ice thickness

The traditional method for interpolating ice thickness data from airborne radar sounding surveys onto regular grids is to employ geostatistical techniques such as kriging. While this approach provides continuous maps of ice thickness, it generates products that are not consistent with ice flow dynamics and are impractical for high resolution ice flow simulations. Here, we present a novel approach that combines sparse ice thickness data collected by airborne radar sounding profilers with high resolution swath mapping of ice velocity derived from satellite synthetic-aperture interferometry to obtain a high resolution map of ice thickness that conserves mass and minimizes the departure from observations. We apply this approach to the case of Nioghalvfjerdsfjorden (79North) Glacier, a major outlet glacier in northeast Greenland that has been relatively stable in recent decades. The results show that our mass conserving method removes the anomalies in mass flux divergence, yields interpolated data that are within about 5% of the original data, and produces thickness maps that are directly usable in high spatial-resolution, high-order ice flow models. We discuss the application of this method to the broad and detailed radar surveys of ice sheet and glacier thickness. Copyright 2011 by the American Geophysical Union.

Insights into spatial sensitivities of ice mass response to environmental change from the SeaRISE ice sheet modeling project I: Antarctica

Sophie Nowicki, Robert A. Bindschadler, Ayako Abe-Ouchi, Andy Aschwanden, Ed Bueler, Hyeungu Choi, Jim Fastook, Glen Granzow, Ralf Greve, Gail Gutowski, Ute Herzfeld, Charles Jackson, Jesse Johnson, Constantine Khroulev, Eric Larour, Anders Levermann, William H. Lipscomb, Maria A. Martin, Mathieu Morlighem, Byron R. Parizek, David Pollard, Stephen F. Price, Diandong Ren, Eric Rignot, Fuyuki Saito, Tatsuru Sato, Hakime Seddik, Helene Seroussi, Kunio Takahashi, Ryan Walker, and Wei Li Wang

Radiostratigraphy and age structure of the Greenland Ice Sheet

Several decades of ice-penetrating radar surveys of the Greenland and Antarctic ice sheets have observed numerous widespread internal reflections. Analysis of this radiostratigraphy has produced valuable insights into ice sheet dynamics and motivates additional mapping of these reflections. Here we present a comprehensive deep radiostratigraphy of the Greenland Ice Sheet from airborne deep ice-penetrating radar data collected over Greenland by The University of Kansas between 1993 and 2013. To map this radiostratigraphy efficiently, we developed new techniques for predicting reflection slope from the phase recorded by coherent radars. When integrated along track, these slope fields predict the radiostratigraphy and simplify semiautomatic reflection tracing. Core-intersecting reflections were dated using synchronized depth-age relationships for six deep ice cores. Additional reflections were dated by matching reflections between transects and by extending reflection-inferred depth-age relationships using the local effective vertical strain rate. The oldest reflections, dating to the Eemian period, are found mostly in the northern part of the ice sheet. Within the onset regions of several fast-flowing outlet glaciers and ice streams, reflections typically do not conform to the bed topography. Disrupted radiostratigraphy is also observed in a region north of the Northeast Greenland Ice Stream that is not presently flowing rapidly. Dated reflections are used to generate a gridded age volume for most of the ice sheet and also to determine the depths of key climate transitions that were not observed directly. This radiostratigraphy provides a new constraint on the dynamics and history of the Greenland Ice Sheet. Key Points Phase information predicts reflection slope and simplifies reflection tracing Reflections can be dated away from ice cores using a simple ice flow model Radiostratigraphy is often disrupted near the onset of fast ice flow

A damage mechanics assessment of the Larsen B ice shelf prior to collapse: Toward a physically-based calving law

Calving is a primary process of mass ablation for glaciers and ice sheets, though it still eludes a general physical law. Here, we propose a calving framework based on continuum damage mechanics coupled with the equations of viscous deformation of glacier ice. We introduce a scalar damage variable that quantifies the loss of load-bearing surface area due to fractures and that feeds back with ice viscosity to represent fracture-induced softening. The calving law is a standard failure criterion for viscous damaging materials and represents a macroscopic brittle instability quantified by a critical or threshold damage. We constrain this threshold using the Ice Sheet System Model (ISSM) by inverting for damage on the Larsen B ice shelf prior to its 2002 collapse. By analyzing the damage distribution in areas that subsequently calved, we conclude that calving occurs after fractures have reduced the load-bearing capacity of the ice by 60 ± 10%.

Fast retreat of Zachariæ Isstrøm, northeast Greenland

Shrinking shelf and faster flow Zachariæ Isstrøm, a large glacier in northeast Greenland, began a rapid retreat after detaching from a stabilizing sill in the late 1990s. Mouginot et al. report that between 2002 and 2014, the area covered by the glaciers ice shelf shrank by 95%; since 1999, the glaciers flow rate has nearly doubled; and its acceleration increased threefold in the fall of 2012. These dramatic changes appear to be the result of a combination of warmer air and ocean temperatures and the topography of the ocean floor at the head of the glacier. Rising sea levels should continue to destabilize the marine portion of Zachariæ Isstrøm for decades. Science, this issue p. 1357 A large glacier in northeast Greenland is retreating rapidly as air and ocean warm. After 8 years of decay of its ice shelf, Zachariæ Isstrøm, a major glacier of northeast Greenland that holds a 0.5-meter sea-level rise equivalent, entered a phase of accelerated retreat in fall 2012. The acceleration rate of its ice velocity tripled, melting of its residual ice shelf and thinning of its grounded portion doubled, and calving is now occurring at its grounding line. Warmer air and ocean temperatures have caused the glacier to detach from a stabilizing sill and retreat rapidly along a downward-sloping, marine-based bed. Its equal-ice-volume neighbor, Nioghalvfjerdsfjorden, is also melting rapidly but retreating slowly along an upward-sloping bed. The destabilization of this marine-based sector will increase sea-level rise from the Greenland Ice Sheet for decades to come.

BedMachine v3: Complete Bed Topography and Ocean Bathymetry Mapping of Greenland From Multibeam Echo Sounding Combined With Mass Conservation

Abstract Greenlands bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenlands marine‐terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine‐terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine‐based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

Ice flow sensitivity to geothermal heat flux of Pine Island Glacier, Antarctica

Model projections of ice flow in a changing climate are dependent on model inputs such as surface elevation, bedrock position or surface temperatures, among others. Of all these inputs, geothermal heat flux is the one for which uncertainty is greatest. In the area of Pine Island Glacier, Antarctica, available data sets differ by up to a factor of 2.5. Here, we evaluate the impact of such uncertainty on ice flow, using sampling analyses based on the Latin-Hypercube method. First, we quantify the impact of geothermal heat flux errors on ice hardness, a thermal parameter that critically controls the magnitude of ice flow. Second, we quantify the impact of the same errors on mass balance, specifically on the mass flux advecting through thirteen fluxgates distributed across Pine Island Glacier. We contrast our Results with similar uncertainties generated by errors in the specification of ice thickness. Model outputs indicate that geothermal heat flux errors yield uncertainties on ice hardness on the order of 5-7%, with maximum uncertainty reaching 15%. Resulting uncertainties in mass balance remain however below 1%. We discuss the uncertainty distribution and its relationship to the amount of heat available at the base of the ice sheet from friction, viscous and geothermal heating. We also show that comparatively, errors in ice thickness contribute more to model uncertainty than errors in geothermal heat flux, especially for fast-flowing ice streams.