skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A probabilistic framework for quantifying the role of anthropogenic climate change in marine-terminating glacier retreats
Abstract. Many marine-terminating outlet glaciers have retreated rapidly in recent decades, but these changes have not been formally attributed to anthropogenic climate change. A key challenge for such an attribution assessment is that if glacier termini are sufficiently perturbed from bathymetric highs, ice-dynamic feedbacks can cause rapid retreat even without further climate forcing. In the presence of internal climate variability, attribution thus depends on understanding whether (or how frequently) these rapid retreats could be triggered by climatic noise alone. Our simulations with idealized glaciers show that in a noisy climate, rapid retreat is a stochastic phenomenon. We therefore propose a probabilistic approach to attribution and present a framework for analysis that uses ensembles of many simulations with independent realizations of random climate variability. Synthetic experiments show that century-scale climate trends substantially increase the likelihood of rapid glacier retreat. This effect depends on the timescales over which ice dynamics integrate forcing. For a population of synthetic glaciers with different topographies, we find that external trends increase the number of large retreats triggered within the population, offering a metric for regional attribution. Our analyses suggest that formal attribution studies are tractable and should be further pursued to clarify the human role in recent ice-sheet change. We emphasize that early-industrial-era constraints on glacier and climate state are likely to be crucial for such studies.  more » « less
Award ID(s):
1947882
PAR ID:
10407345
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
16
Issue:
7
ISSN:
1994-0424
Page Range / eLocation ID:
2725 to 2743
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Increasing ice flux from glaciers retreating over deepening (retrograde) 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 persistence of a glacier at a bed peak is caused by ice slowing as it flows up a reverse-sloping bed to the peak. Persistence at bed peaks may lead to two very different future behaviors for a glacier: one where it persists at a bed peak indefinitely, and another where it retreats from the bed peak after potentially long delays following climate forcing. However, it is nearly impossible to distinguish which of these two future behaviors 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. We recommend that further research is needed on seemingly stable glaciers to determine their likely future. 
    more » « less
  2. Recent observations indicate that many marine‐terminating glaciers in Greenland and Antarctica are currently retreating and thinning, potentially due to long‐term trends in climate forcing. In this study, we describe a simple two‐stage model that accurately emulates the response to external forcing of marine‐terminating glaciers simulated in a spatially extended model. The simplicity of the model permits derivation of analytical expressions describing the marine‐terminating glacier response to forcing. We find that there are two time scales that characterize the stable glacier response to external forcing, a fast time scale of decades to centuries, and a slow time scale of millennia. These two time scales become unstable at different thresholds of bed slope, indicating that there are distinct slow and fast forms of the marine ice sheet instability. We derive simple expressions for the approximate magnitude and transient evolution of the stable glacier response to external forcing, which depend on the equilibrium glacier state and the strength of nonlinearity in forcing processes. The slow response rate of marine‐terminating glaciers indicates that current changes at some glaciers are set to continue and accelerate in coming centuries in response to past climate forcing and that the current extent of change at these glaciers is likely a small fraction of the future committed change caused by past climate forcing. Finally, we find that changing the amplitude of natural fluctuations in some nonlinear forcing processes, such as ice shelf calving, changes the equilibrium glacier state. 
    more » « less
  3. Abstract. Observations of recent mass loss rates of the West Antarctic Ice Sheet (WAIS) raise concerns about its stability since a collapse would increase global sea levels by several meters. Future projections of these mass loss trends are often estimated using numerical ice sheet models, and recent studies have highlighted the need for models to be benchmarked against present-day observed mass change rates. Here, we present an improved initialization method that optimizes local agreement not only with observations of ice thickness and surface velocity but also with satellite-based estimates of mass change rates. This is achieved by a combination of tuned thermal forcing under the floating ice shelves and friction under the ice sheet. Starting from this improved present-day state, we generate an ensemble of future simulations of Antarctic mass change by varying model physical choices and parameter values while fixing the climate forcing at present-day values. The dynamical response shows slow grounding-line retreat over several centuries, followed by a phase of rapid mass loss over about 200 years with a consistent rate of ∼3 mm GMSL yr−1 (global mean sea level). We find that, for all ensemble members, the Thwaites Glacier and Pine Island Glacier collapse. Our results imply that present-day ocean thermal forcing, if held constant over multiple centuries, may be sufficient to deglaciate large parts of the WAIS, raising global mean sea level by at least a meter. 
    more » « less
  4. Projections of the growth and demise of ice sheets and glaciers require physical models of the processes governing flow and fracture of ice. The flow of glacier ice has been treated using increasingly sophisticated models. By contrast, fracture, the process ultimately responsible for half of the mass lost from ice sheets through iceberg calving, is often included using ad hoc parameterizations. In this study we seek to bridge this gap by introducing a model where ice obeys a power law rheology appropriate for intact ice below a yield strength. Above the yield strength, we introduce a separate rheology appropriate for the flow of heavily fractured ice, where ice deformation occurs more readily along faults and fractures. We show that, provided the motion of fractured ice is sufficiently rapid compared to that of intact ice, the behavior of glaciers depends solely on the rheology of intact ice and the yield strength of ice and is insensitive to the precise rheology of fractured ice. Moreover, assuming that glacier ice is unyielded allows us to bound the long‐term average rate of terminus advance, providing a first principles estimate of rates of retreat associated with the marine ice cliff instability. We illustrate model behavior using idealized geometries and climate forcing and show that the model not only exhibits realistic patterns of advance and retreat but also has the potential to exhibit hysteresis. This hysteresis could provide an explanation for the sudden onset of rapid retreat observed in marine‐terminating glaciers. 
    more » « less
  5. null (Ed.)
    Abstract. Changes in glacier length reflect the integrated response to local fluctuations in temperature and precipitation resulting from both external forcing (e.g., volcanic eruptions or anthropogenic CO2) and internal climate variability. In order to interpret the climate history reflected in the glacier moraine record, the influence of both sources of climate variability must therefore be considered. Here we study the last millennium of glacier-length variability across the globe using a simple dynamic glacier model, which we force with temperature and precipitation time series from a 13-member ensemble of simulations from a global climate model. The ensemble allows us to quantify the contributions to glacier-length variability from external forcing (given by the ensemble mean) and internal variability (given by the ensemble spread). Within this framework, we find that internal variability is the predominant source of length fluctuations for glaciers with a shorter response time (less than a few decades). However, for glaciers with longer response timescales (more than a few decades) external forcing has a greater influence than internal variability. We further find that external forcing also dominates when the response of glaciers from widely separated regions is averaged. Single-forcing simulations indicate that, for this climate model, most of the forced response over the last millennium, pre-anthropogenic warming, has been driven by global-scale temperature change associated with volcanic aerosols. 
    more » « less