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: The demise of the world's largest piedmont glacier: a probabilistic forecast
Abstract. Sít' Tlein, located in the St. Elias Range, which straddles Alaska's Wrangell–St. Elias National Park and Kluane National Park in the Yukon, is the world's largest piedmont glacier. Sít' Tlein has thinned considerably over 30 years of altimetry, yet its low-elevation piedmont lobe has remained intact in contrast to the glaciers that once filled neighboring Icy and Disenchantment bays. In an effort to forecast changes to Sít' Tlein over decadal to centennial timescales, we take a data-constrained dynamical modeling approach in which we infer the parameters of a higher-order model of ice flow – the bed elevation, basal traction, and surface mass balance – with a diverse but spatiotemporally sparse set of observations including satellite-derived, time-varying velocity fields; radar-derived bed and surface elevation measurements; and in situ and remotely sensed observations of accumulation and ablation. Nonetheless, such data do not uniquely constrain model behavior, so we adopt an approximate Bayesian approach based on the Laplace approximation and facilitated by low-rank parametric representations to quantify uncertainty in the bed, traction, and mass balance fields alongside the induced uncertainty in model-based predictions of glacier change. We find that Sít' Tlein is considerably out of balance with contemporary (and presumably future) climate, and we expect its piedmont lobe to largely disappear over the coming centuries. If warming ceases, and surface mass balance remains at 2023 levels, then by 2073 (2173) we forecast a mass loss (expressed in terms of 95 % credible interval) of 323–444 km3 (546–728 km3). If instead surface mass balance continues to change at the same rate as inferred over the historical period, then we forecast a 2073 (2173) mass loss of 383–505 km3 (740–900 km3). In either case, the resulting retreat and subsequent replacement of glacier ice with a marine embayment or lake will yield a significant modification to the regional landscape and ecosystem.  more » « less
Award ID(s):
1929718 1929566
PAR ID:
10656488
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
The Cryosphere
Date Published:
Journal Name:
The Cryosphere
Volume:
19
Issue:
6
ISSN:
1994-0424
Page Range / eLocation ID:
2321 to 2353
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ### Access GPS, photo, and position data can be accessed and downloaded from the directory via: [http://arcticdata.io/data/10.18739/A2MS3K37K](http://arcticdata.io/data/10.18739/A2MS3K37K). ### Overview The following files gather Surface Mass Balance stakes, timelapse camera photos, Hobo TidBits temperature sensors data, Global Positioning System (GPS) measurements, and debris thickness measurements on the surface of Malaspina Glacier, specifically in its lobe. The data was collected in Summer 2021. The glacier is located on the southwestern coast of Alaska, in the Wrangell St-Elias National Park. The study was mostly focused on the lobe of this piedmont glacier. The goal of the study was to measure several glaciological variables to be later compared to satellite-derived datasets, and incorporated into an ice flow model. The research topic was the retreat of Malaspina glacier over the next decades. We identified regions of interests within the glacier's lobe, that could potentially be markers of a starting retreat. The data gathered was used for publications currently in preparation as of October 2025. This study was led under the NSF grant number #1929566. 
    more » « less
  2. Abstract Sít’ Tlein (Malaspina Glacier), located in Southeast Alaska, has a complex flow history. This piedmont glacier, the largest in the world, is fed by three main tributaries that all exhibit similar flow patterns, yet with varying surge cycles. The piedmont lobe is dramatically reshaped by surges that occur at approximately decadal timescales. By combining historical accounts with modern remote sensing data, we derive a surge history over the past century. We leverage the Stochastic Matrix Factorization, a novel data analysis and interpolation technique, to process and interpret large datasets of glacier surface velocities. A variant of the Principal Component Analysis allows us to uncover spatial and temporal patterns in ice dynamics. We show that Sít’ Tlein displays a wide range of behaviors, spanning quiescence to surge with seasonal to decadal variations of ice flow direction and magnitude. We find that in the lobe, surges dominate the velocity dataset’s variance (spanning 1984–2021), while seasonal variations represent a much smaller part of the variance. However, despite the regular surge pulses, the glacier lobe is far from equilibrium, and widespread retreat of the glacier is inevitable, even without further climate warming. 
    more » « less
  3. null (Ed.)
    Abstract Glacier surges are periodic episodes of mass redistribution characterized by dramatic increases in ice flow velocity and, sometimes, terminus advance. We use optical satellite imagery to document five previously unexamined surge events of Sít’ Kusá (Turner Glacier) in the St. Elias Mountains of Alaska from 1983 to 2013. Surge events had an average recurrence interval of ~5 years, making it the shortest known regular recurrence interval in the world. Surge events appear to initiate in the winter, with speeds reaching up to ~25 m d −1 . The surges propagate down-glacier over ~2 years, resulting in maximum thinning of ~100 m in the reservoir zone and comparable thickening at the terminus. Collectively, the rapid recurrence interval, winter initiation and down-glacier propagation suggest Sít’ Kusá's surges are driven by periodic changes in subglacial hydrology and glacier sliding. Elevation change observations from the northern tributary show a kinematic disconnect above and below an icefall located 23 km from the terminus. We suggest the kinematic disconnect inhibits drawdown from the accumulation zone above the icefall, which leads to a steady flux of ice into the reservoir zone, and contributes to the glacier's exceptionally short recurrence interval. 
    more » « less
  4. Abstract We used measurements of radar-detected stratigraphy, surface ice-flow velocities and accumulation rates to investigate relationships between local valley-glacier and regional ice-sheet dynamics in and around the Schmidt Hills, Pensacola Mountains, Antarctica. Ground-penetrating radar profiles were collected perpendicular to the long axis of the Schmidt Hills and the margin of Foundation Ice Stream (FIS). Within the valley confines, the glacier consists of blue ice, and profiles show internal stratigraphy dipping steeply toward the nunataks and truncated at the present-day ablation surface. Below the valley confines, the blue ice is overlain by firn. Data show that upward-progressing overlap of actively accumulating firn onto valley-glacier ice is slightly less than ice flow out of the valleys over the past ∼1200 years. The apparent slightly negative mass balance (-0.25 cm a -1 ) suggests that ice-margin elevations in the Schmidt Hills may have lowered over this time period, even without a change in the surface elevation of FIS. Results suggest that (1) mass-balance gradients between local valley glaciers and regional ice sheets should be considered when using local information to estimate regional ice surface elevation changes; and (2) interpretation of shallow ice structures imaged with radar can provide information about local ice elevation changes and stability. 
    more » « less
  5. Abstract Malaspina Glacier, located on the coast of southern Alaska, is the world's largest piedmont glacier. A narrow ice‐cored foreland zone undergoing rapid thermokarst erosion separates the glacier from the relatively warm waters of the Gulf of Alaska. Glacier‐wide thinning rates for Malaspina are greater than 1 m/yr, and previous geophysical investigations indicated that bed elevation exceeds 300 m below sea level in some places. These observations together give rise to the question of glacial stability. To address this question, glacier evolution models are dependent upon detailed observations of Malaspina's subglacial topography. Here, we map 2,000 line‐km of the glacier's bed using airborne radar sounding data collected by NASA's Operation IceBridge. When compared to gridded radar measurements, we find that glaciological models overestimate Malaspina's volume by more than 30%. While we report a mean bed elevation 100 m greater than previous models, we find that Malaspina inhabits a broad basin largely grounded below sea level. Several subglacial channels dissect the glacier's bed: the most prominent of these channels extends at least 35 km up‐glacier from the terminus toward the throat of Seward Glacier. Provided continued foreland erosion, an ice‐ocean connection may promote rapid retreat along these overdeepened subglacial channels, with a global sea‐level rise potential of 1.4 mm. 
    more » « less