A document by Enze Zhang. Click on the document to view its contents.

Increasing ice flux from glaciers retreating over deepening bed topography has been implicated in the recent acceleration of mass loss from the Greenland and Antarctic ice sheets. We show in observations that some glaciers have remained at peaks in bed topography without retreating despite enduring significant changes in climate. Observations also indicate that some glaciers which persist at bed peaks undergo sudden retreat years or decades after the onset of local ocean or atmospheric warming. Using model simulations, we show that glacier persistence may lead to two very different futures: one where glaciers persist at bed peaks indefinitely, and another where glaciers retreat from the bed peak suddenly without a concurrent climate forcing. However, it is difficult to distinguish which of these two futures will occur from current observations. We conclude that inferring glacier stability from observations of persistence obscures our true commitment to future sea-level rise under climate change.

The Greenland Ice Sheet discharges ice to the ocean through hundreds of outlet glaciers. Recent acceleration of Greenland outlet glaciers has been linked to both oceanic and atmospheric drivers. Here, we leverage temporally dense observations, regional climate model output, and newly developed time series analysis tools to assess the most important forcings causing ice flow variability at one of the largest Greenland outlet glaciers, Helheim Glacier, from 2009 to 2017. We find that ice speed correlates most strongly with catchment-integrated runoff at seasonal to interannual scales, while multi-annual flow variability correlates most strongly with multi-annual terminus variability. The multiple relevant time scales and the influence of subglacial topography on Helheim Glacier’s dynamics highlight different regimes that can inform modeling and forecasting of its future. Notably, our results suggest that the recent terminus history observed at Helheim is a response to, rather than the cause of, upstream changes.

Projections of the sea level contribution from the Greenland and Antarctic ice sheets rely on atmospheric and oceanic drivers obtained from climate models. The Earth System Models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6) generally project greater future warming compared with the previous CMIP5 effort. Here we use four CMIP6 models and a selection of CMIP5 models to force multiple ice sheet models as part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We find that the projected sea level contribution at 2100 from the ice sheet model ensemble under the CMIP6 scenarios falls within the CMIP5 range for the Antarctic ice sheet but is significantly increased for Greenland. Warmer atmosphere in CMIP6 models results in higher Greenland mass loss due to surface melt. For Antarctica, CMIP6 forcing is similar to CMIP5 and mass gain from increased snowfall counteracts increased loss due to ocean warming.

Numerical, process-based simulations of tidewater glacier evolution are necessary to project future sea-level change under various climate scenarios. Previous work has shown that nonlinearities in tidewater glacier and ice stream dynamics can lead to biases in simulated ice mass change in the presence of noisy forcings. Ice sheet modeling projections that will be used in the upcoming IPCC Assessment Report 6 (AR6) utilize atmospheric and oceanic forcings at annual temporal resolution, omitting any higher frequency forcings. Here, we quantify the effect of seasonal (<1 year) tidewater glacier terminus oscillations on decadal-scale (30 year) mass change. We use an idealized geometry to mimic realistic tidewater glacier geometries, and investigate the impact of the magnitude of seasonal oscillations, bed slope at the glacier terminus, and basal friction law. We find that omitting seasonal terminus motion results in biased mass change projections, with up to an 18% overestimate of mass loss when seasonality is neglected. The bias is most sensitive to the magnitude of the seasonal terminus oscillations and exhibits very little sensitivity to choice of friction law. Our results show that including seasonality is required to eliminate a potential bias in ice sheet mass change projections. In order to achieve this, seasonality in atmospheric and oceanic forcings must be adequately represented and observations of seasonal terminus positions and tidewater glacier thickness changes must be acquired to evaluate numerical models.