Using surface velocities to calculate ice thickness and bed topography: A case study at Columbia Glacier, Alaska, USA

R.W. McNabb; R. Hock; Shad O’Neel; Lowell A. Rasmussen; Y. Ahn; M. Braun; H. Conway; S. Herreid; I. Joughin; W.T. Pfeffer; B.E. Smith; M. Truffer

doi:10.3189/2012jog11j249

USGS · 70041058

Using surface velocities to calculate ice thickness and bed topography: A case study at Columbia Glacier, Alaska, USA

R.W. McNabb·R. Hock·Shad O’Neel·Lowell A. Rasmussen·Y. Ahn·M. Braun·H. Conway·S. Herreid·I. Joughin·W.T. Pfeffer·B.E. Smith·M. Truffer

View original source ↗Read source document ↗Export metadata ↓

Abstract

Information about glacier volume and ice thickness distribution is essential for many glaciological applications, but direct measurements of ice thickness can be difficult and costly. We present a new method that calculates ice thickness via an estimate of ice flux. We solve the familiar continuity equation between adjacent flowlines, which decreases the computational time required compared to a solution on the whole grid. We test the method on Columbia Glacier, a large tidewater glacier in Alaska, USA, and compare calculated and measured ice thicknesses, with favorable results. This shows the potential of this method for estimating ice thickness distribution of glaciers for which only surface data are available. We find that both the mean thickness and volume of Columbia Glacier were approximately halved over the period 1957–2007, from 281m to 143 m, and from 294 km 3 to 134 km 3 , respectively. Using bedrock slope and considering how waves of thickness change propagate through the glacier, we conduct a brief analysis of the instability of Columbia Glacier, which leads us to conclude that the rapid portion of the retreat may be nearing an end.

Explore related subjects

Journal of Glaciology Article

Source-reported bounding extent: 52.68° to 62.54° latitude; -155.11° to -141.02° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

R.W. McNabb, R. Hock, Shad O’Neel, Lowell A. Rasmussen, Y. Ahn, M. Braun, H. Conway, S. Herreid, I. Joughin, W.T. Pfeffer, B.E. Smith, M. Truffer. 2017-09-08. Using surface velocities to calculate ice thickness and bed topography: A case study at Columbia Glacier, Alaska, USA. https://doi.org/10.3189/2012jog11j249

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

The Python Energy Balance model for Snow and Ice (PEBSI): Application and tradeoff analysis on Gulkana Glacier, Alaska

Glacier energy-balance models offer mechanistic insights into glacier mass balance under a changing climate, yet their considerable data requirements hinder large-scale applications. Here we present the open-source Python Energy Balance model for Snow and Ice (PEBSI), which includes physically based albedo evolution using the Snow, Ice and Aerosol Radiative (SNICAR) model. PEBSI is calibrated and validated using robust in situ data from Gulkana Glacier, Alaska from 2000 to 2024. Simulations forced with original and bias-corrected climate reanalysis data show that statistically downscaling reanalysis data with in situ observations is necessary to reproduce summer mass balance (mean absolute error [MAE] = 0.75 m w.e. vs 0.22 m w.e., respectively). A grid search across two parameters, a precipitation factor and a densification parameter, reveals tradeoffs in performance compared to seasonal mass balance and end-of-winter snow density and depth. No single combination of parameters minimizes all errors, underscoring the inherent overparameterization of energy-balance models and challenges with translating coarse climate data to the glacier scale. The calibrated model successfully simulates the 2024 melt season, agreeing with surface-height change (MAE = 0.48 m) and albedo (MAE = 0.066) observations. Moving forward, PEBSI provides unique opportunities to quantify albedo feedbacks and their impact on present and future glacier mass loss.

Alaska↗

USGS2025

Equilibrium line altitudes, accumulation areas, and the vulnerability of glaciers in Alaska

The accumulation area ratio (AAR) of a glacier reflects its current state of equilibrium, or disequilibrium, with climate and its vulnerability to future climate change. Here, we present an inventory of glacier-specific annual accumulation areas and equilibrium line altitudes (ELAs) for over 3000 glaciers in Alaska and northwest Canada (88% of the regional glacier area) from 2018 to 2022 derived from Sentinel-2 imagery. We find that the 5 year average AAR of the entire study area is 0.41, with an inter-annual range of 0.25–0.49. More than 1000 glaciers, representing 8% of the investigated glacier area, were found to have effectively no accumulation area. Summer temperature and winter precipitation from ERA5-Land explained nearly 50% of the inter-annual ELA variability across the entire study region ( R 2 = 0.47 ). An analysis of future climate scenarios (SSP2-4.5) projects that ELAs will rise by ∼170 m on average by the end of the 21st century. Such changes would result in a loss of 25% of the modern accumulation area, leaving a total of 1900 glaciers (22% of the investigated area) with no accumulation area. These results highlight the current state of glacier disequilibrium with modern climate, as well as glacier vulnerability to projected future warming.

Alaska↗

USGS2024

GNSS reflectometry from low-cost sensors for continuous in situ contemporaneous glacier mass balance and flux divergence

Recent advances in remote sensing have produced global glacier surface elevation change data. Parsing these elevation change signals into contributions from the climate (i.e. climatic mass balance) and glacier dynamics (i.e. flux divergence) is critical to enhance our process-based understanding of glacier change. In this study, we evaluate three approaches for direct, continuous measurements of the climatic mass balance, flux divergence, and elevation change at a site on Gulkana Glacier in Alaska using low-cost GNSS sensors, GNSS interferometric reflectometry (GNSS-IR), banded ablation stakes with time-lapse cameras, and combinations thereof. Cumulative climatic mass balance over the season was 4.85 m and the three approaches were within 0.08 m through early July before the snowpack melted, and within 0.28 m through mid-August. The flux divergence increased from 0.52 ± 0.03 cm d -1 before June 3 to roughly 0.73 cm d -1 after June 27. We demonstrate a single GNSS system fixed atop an ablation stake can measure contemporaneous climatic mass balance, flux divergence, and elevation change based on the antenna’s position and GNSS-IR techniques. The ability of these systems to measure glacier mass balance and flux divergence offers unique opportunities for year-round observations on mountain glaciers in the future.

Alaska↗