Speedup of Pine Island Glacier over the past several decades has made it Antarctica’s largest contributor to sea-level rise. The past speedup is largely due to grounding-line retreat in response to ocean-induced thinning that reduced ice-shelf buttressing. While speeds remained fairly steady from 2009 to late 2017, our Copernicus Sentinel 1A/B–derived velocity data show a >12% speedup over the past 3 years, coincident with a 19-km retreat of the ice shelf. We use an ice-flow model to simulate this loss, finding that accelerated calving can explain the recent speedup, independent of the grounding-line, melt-driven processes responsible for past speedups. If the ice shelf’s rapid retreat continues, it could further destabilize the glacier far sooner than would be expected due to surface- or ocean-melting processes.

Resistance at the ice-bed interface provides a strong control on the response of ice streams and outlet glaciers to external forcing, yet it is not observable by remote sensing. We used inverse methods constrained by satellite observations to infer the basal resistance to flow underneath three of the Greenland Ice Sheet's largest outlet glaciers. In regions of fast ice flow and high (>250 kPa) driving stresses, ice is often assumed to flow over a strong bed. We found, however, that the beds of these three glaciers provide almost no resistance under the fast-flowing trunk. Instead, resistance to flow is provided by the lateral margins and stronger beds underlying slower-moving ice upstream. Additionally, we found isolated patches of high basal resistivity within the predominantly weak beds. Because these small-scale (<1 ice thickness) features may be artifacts of overfitting our solution to measurement errors, we tested their robustness to different degrees of regularization.

Abstract. We introduce a new software package called “icepack” for modeling the flow of glaciers and ice sheets. The icepack package is built on the finite element modeling library Firedrake, which uses the Unified Form Language (UFL), a domain-specific language embedded into Python for describing weak forms of partial differential equations. The diagnostic models in icepack are formulated through action principles that are specified in UFL. The components of each action functional can be substituted for different forms of the user's choosing, which makes it easy to experiment with the model physics. The action functional itself can be used to define a solver convergence criterion that is independent of the mesh and requires little tuning on the part of the user. The icepack package includes the 2D shallow ice and shallow stream models. We have also defined a 3D hybrid model based on spectral semi-discretization of the Blatter–Pattyn equations. Finally, icepack includes a Gauss–Newton solver for inverse problems that runs substantially faster than the Broyden–Fletcher–Goldfarb–Shanno (BFGS) method often used in the glaciological literature. The overall design philosophy of icepack is to be as usable as possible for a wide a swath of the glaciological community, including both experts and novices in computational science.