USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Glaciology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
USGS2017

Röthlisberger channel theory: its origins and consequences

The theory of channelized water flow through glaciers, most commonly associated with the names of Hans Röthlisberger and Ron Shreve and their 1972 papers in the Journal of Glaciology, was developed at a time when interest in glacier-bed processes was expanding, and the possible relationship between glacier sliding and water at the bed was becoming of keen interest. The R-channel theory provided for the first time a physically based conceptual model of water flow through glaciers. The theory also marks the emergence of glacier hydrology as a glaciological discipline with goals and methods distinct from those of surface-water hydrology.

Journal of Glaciology↗
USGS2007

Distributed energy-balance modeling of snow-cover evolution and melt in rugged terrain: Tobacco Root Mountains, Montana, USA

A distributed energy-balance model was developed for simulating snowpack evolution and melt in rugged terrain. The model, which was applied to a 43-km2 watershed in the Tobacco Root Mountains, Montana, USA, used measured ambient data from nearby weather stations to drive energy-balance calculations and to constrain the model of Liston and Sturm [Liston, G.E., Sturm, M., 1998. A snow-transport model for complex terrain. Journal of Glaciology 44 (148), 498-516] for calculating the initial snowpack thickness. Simulated initial snow-water equivalent ranged between 1 cm and 385 cm w.e. (water equivalent) with high values concentrated on east-facing slopes below tall summits. An interpreted satellite image of the snowcover distribution on May 6, 1998, closely matched the simulated distribution with the greatest discrepancy occurring in the floor of the main trunk valley. Model simulations indicated that snowmelt commenced early in the melt season, but rapid meltout of snow cover did not occur until after the average energy balance of the entire watershed became positive about 45 days into the melt season. Meltout was fastest in the lower part of the watershed where warmer temperatures and tree cover enhanced the energy income of the underlying snow. An interpreted satellite image of the snowcover distribution on July 9, 1998 compared favorably with the simulated distribution, and melt curves for modeled canopy-covered cells mimicked the trends measured at nearby snow pillow stations. By the end of the simulation period (August 3), 28% of the watershed remained snow covered, most of which was concentrated in the highest parts of the watershed where initially thick accumulations had been shaded by surrounding summits. The results of this study provide further demonstration of the critical role that topography plays in the timing and magnitude of snowmelt from high mountain watersheds. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Hydrology↗
USGS2022

Beyond glacier-wide mass balances: Parsing seasonal elevation change into spatially resolved patterns of accumulation and ablation at Wolverine Glacier, Alaska

We present spatially distributed seasonal and annual surface mass balances of Wolverine Glacier, Alaska, from 2016 to 2020. Our approach accounts for the effects of ice emergence and firn compaction on surface elevation changes to resolve the spatial patterns in mass balance at 10 m scale. We present and compare three methods for estimating emergence velocities. Firn compaction was constrained by optimizing a firn model to fit three firn cores. Distributed mass balances showed good agreement with mass-balance stakes (RMSE = 0.67 m w.e., r = 0.99, n = 41) and ground-penetrating radar surveys (RMSE = 0.36 m w.e., r = 0.85, n = 9024). Fundamental differences in the distributions of seasonal balances highlight the importance of disparate physical processes, with anomalously high ablation rates observed in icefalls. Winter balances were found to be positively skewed when controlling for elevation, while summer and annual balances were negatively skewed. We show that only a small percent of the glacier surface represents ideal locations for mass-balance stake placement. Importantly, no suitable areas are found near the terminus or in elevation bands dominated by icefalls. These findings offer explanations for the often-needed geodetic calibrations of glaciological time series.

Alaska↗
USGS2017

Geometry, mass balance and thinning at Eklutna Glacier, Alaska: an altitude-mass-balance feedback with implications for water resources

We analyzed glacier surface elevations (1957, 2010 and 2015) and surface mass-balance measurements (2008–2015) on the 30 km 2 Eklutna Glacier, in the Chugach Mountains of southcentral Alaska. The geodetic mass balances from 1957 to 2010 and 2010 to 2015 are −0.52 ± 0.46 and −0.74 ± 0.10 m w.e. a −1 , respectively. The glaciological mass balance of −0.73 m w.e. a −1 from 2010 to 2015 is indistinguishable from the geodetic value. Even after accounting for loss of firn in the accumulation zone, we found most of the mass loss over both time periods was from a broad, low-slope basin that includes much of the accumulation zone of the main branch. Ice-equivalent surface elevation changes in the basin were −1.0 ± 0.8 m a −1 from 1957 to 2010, and −0.6 ± 0.1 m a −1 from 2010 to 2015, shifting the glacier hypsometry downward and resulting in more negative mass balances: an altitude-mass-balance feedback. Net mass loss from Eklutna Glacier accounts for 7 ± 1% of the average inflow to Eklutna Reservoir, which is entirely used for water and power by Anchorage, Alaska's largest city. If the altitude-mass-balance feedback continues, this ‘deglaciation discharge dividend’ is likely to increase over the short-term before it eventually decreases due to diminishing glacier area.

Alaska↗
USGS2019

Reanalysis of the U.S. Geological Survey Benchmark Glaciers: Long-term insight into climate forcing of glacier mass balance

Mountain glaciers integrate climate processes to provide an unmatched signal of regional climate forcing. However, extracting the climate signal via intercomparison of regional glacier mass balance records can be problematic when methods for extrapolating and calibrating direct glaciological measurements are mixed or inconsistent. To address this problem, we reanalyzed and compared long-term mass balance records from the U.S. Geological Survey Benchmark Glaciers. These five glaciers span maritime and continental climate regimes of the western United States and Alaska. Each glacier exhibits cumulative mass loss since the mid-20th century, with average rates ranging from –0.58 to –0.30 m water equivalent (w.e.) a-1. We produced a set of solutions using different extrapolation and calibration methods to inform uncertainty estimates, which range from 0.22–0.44 m w.e. a-1. Mass losses are primarily driven by increasing summer warming. Continentality exerts a stronger control on mass loss than latitude. Similar to elevation, topographic shading, snow redistribution, and glacier surface features often exert first-order control on mass balance. The reanalysis underscores the value of geodetic calibration to resolve mass balance magnitude, as well as the irreplaceable value of direct measurements in contributing to process-based understanding of glacier mass balance.

Journal of Glaciology↗
USGS2016

Observations and modeling of fjord sedimentation during the 30 year retreat of Columbia Glacier, AK

To explore links between glacier dynamics, sediment yields and the accumulation of glacial sediments in a temperate setting, we use extensive glaciological observations for Columbia Glacier, Alaska, and new oceanographic data from the fjord exposed during its retreat. High-resolution seismic data indicate that 3.2 × 10 8 m 3 of sediment has accumulated in Columbia Fjord over the past three decades, which corresponds to ~5 mm a −1 of erosion averaged over the glaciated area. We develop a general model to infer the sediment-flux history from the glacier that is compatible with the observed retreat history, and the thickness and architecture of the fjord sediment deposits. Results reveal a fivefold increase in sediment flux from 1997 to 2000, which is not correlated with concurrent changes in ice flux or retreat rate. We suggest the flux increase resulted from an increase in the sediment transport capacity of the subglacial hydraulic system due to the retreat-related steepening of the glacier surface over a known subglacial deep basin. Because variations in subglacial sediment storage can impact glacial sediment flux, in addition to changes in climate, erosion rate and glacier dynamics, the interpretation of climatic changes based on the sediment record is more complex than generally assumed.

Alaska↗
USGS2017

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

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.

Alaska↗
USGS2017

Analysis of a GRACE global mascon solution for Gulf of Alaska glaciers

We present a high-resolution Gravity Recovery and Climate Experiment (GRACE) mascon solution for Gulf of Alaska (GOA) glaciers and compare this with in situ glaciological, climate and other remote-sensing observations. Our GRACE solution yields a GOA glacier mass balance of –65 ± 11 Gt a –1 for the period December 2003 to December 2010, with summer balances driving the interannual variability. Between October/November 2003 and October 2009 we obtain a mass balance of –61 ± 11 Gt a –1 from GRACE, which compares well with –65 ± 12 Gt a –1 from ICESat based on hypsometric extrapolation of glacier elevation changes. We find that mean summer (June–August) air temperatures derived from both ground and lower-troposphere temperature records were good predictors of GRACE-derived summer mass balances, capturing 59% and 72% of the summer balance variability respectively. Large mass losses during 2009 were likely due to low early melt season surface albedos, measured by the Moderate Resolution Imaging Spectroradiometer (MODIS) and likely associated with the 31 March 2009 eruption of Mount Redoubt, southwestern Alaska. GRACE data compared well with in situ measurements at Wolverine Glacier (maritime Alaska), but poorly with those at Gulkana Glacier (interior Alaska). We conclude that, although GOA mass estimates from GRACE are robust over the entire domain, further constraints on subregional and seasonal estimates are necessary to improve fidelity to ground observations.

Alaska↗
USGS2019

Formation of pedestalled, relict lakes on the McMurdo Ice Shelf, Antarctica

Surface debris covers much of the western portion of the McMurdo Ice Shelf and has a strong influence on the local surface albedo and energy balance. Differential ablation between debris-covered and debris-free areas creates an unusual heterogeneous surface of topographically low, high-ablation, and topographically raised (‘pedestalled’), low-ablation areas. Analysis of Landsat and MODIS satellite imagery from 1999 to 2018, alongside field observations from the 2016/2017 austral summer, shows that pedestalled relict lakes (‘pedestals’) form when an active surface meltwater lake that develops in the summer, freezes-over in winter, resulting in the lake-bottom debris being masked by a high-albedo, superimposed, ice surface. If this ice surface fails to melt during a subsequent melt season, it experiences reduced surface ablation relative to the surrounding debris-covered areas of the ice shelf. We propose that this differential ablation, and resultant hydrostatic and flexural readjustments of the ice shelf, causes the former supraglacial lake surface to become increasingly pedestalled above the lower topography of the surrounding ice shelf. Consequently, meltwater streams cannot flow onto these pedestalled features, and instead divert around them. We suggest that the development of pedestals has a significant influence on the surface-energy balance, hydrology and flexure of the ice shelf.

Journal of Glaciology↗
USGS2020

Explaining mass balance and retreat dichotomies at Taku and Lemon Creek Glaciers, Alaska

We reanalyzed mass balance records at Taku and Lemon Creek Glaciers to better understand the relative roles of hypsometry, local climate and dynamics as mass balance drivers. Over the 1946–2018 period, the cumulative mass balances diverged. Tidewater Taku Glacier advanced and gained mass at an average rate of +0.25±0.28 m w.e. a–1, contrasting with retreat and mass loss of –0.60±0.15 m w.e. a-1 at land-terminating Lemon Creek Glacier. The uniform influence of regional climate is demonstrated by strong correlations among annual mass balance and climate data. Regional warming trends forced similar statistically significant decreases in surface mass balance after 1989: –0.83 m w.e. a–1 at Taku Glacier and –0.81 m w.e. a–1 at Lemon Creek Glacier. Divergence in cumulative mass balance arises from differences in glacier hypsometry and local climate. Since 2013 negative mass balance and glacier-wide thinning prevailed at Taku Glacier. These changes initiated terminus retreat, which could increase dramatically if calving begins. The future mass balance trajectory of Taku Glacier hinges on dynamics, likely ending the historic dichotomy between these glaciers.

Alaska↗
USGS2022

Hydrologic modeling of a perennial firn aquifer in southeast Greenland

A conceptual model, based on field observations and assumed physics of a perennial firn aquifer near Helheim Glacier (southeast Greenland), is evaluated via steady-state 2-D simulation of liquid water flow and energy transport with phase change. The simulation approach allows natural representation of flow and energy advection and conduction that occur in vertical meltwater recharge through the unsaturated zone and in lateral flow within the saturated aquifer. Agreement between measured and simulated aquifer geometry, temperature, and recharge and discharge rates confirms that the conceptual field-data-based description of the aquifer is consistent with the primary physical processes of groundwater flow, energy transport and phase change. Factors that are found to control simulated aquifer configuration include surface temperature, meltwater recharge rate, residual total-water saturation and capillary fringe thickness. Simulation analyses indicate that the size of perennial firn aquifers depends primarily on recharge rates from surface snowmelt. Results also imply that the recent aquifer expansion, likely due to a warming climate, may eventually produce lakes on the ice-sheet surface that would affect the surface energy balance.

Journal of Glaciology↗
USGS2023

How to handle glacier area change in geodetic mass balance

Innovations in geodesy enable widespread analysis of glacier surface elevation change and geodetic mass balance. However, coincident glacier area data are less widely available, causing inconsistent handling of glacier area change. Here we quantify the bias introduced into meters water equivalent (m w.e.) specific geodetic mass balance results when using a fixed, maximum glacier area, and illustrate the bias for five North American glaciers. Sites span latitudes from the northern U.S. Rocky Mountains (48°N) to the Central Alaska Range (63°N) between 1948 and 2021. Results show that fixed (maximum) area treatment subdues the m w.e. mass change signal, underestimating mass balance by up to 19% in our test cases. This bias scales with relative glacier area change and the mass balance magnitude. Thus, the bias for specific geodetic mass balances will be most pronounced across rapidly deglaciating regions. Our analysis underscores the need for temporally resolved glacier area in geodetic mass balance studies.

Journal of Glaciology↗
USGS2024

Direct measurements of firn-density evolution from 2016 to 2022 at Wolverine Glacier, Alaska

Knowledge of snow and firn-density change is needed to use elevation-change measurements to estimate glacier mass change. Additionally, firn-density evolution on glaciers is closely connected to meltwater percolation, refreezing and runoff, which are key processes for glacier mass balance and hydrology. Since 2016, the U.S. Geological Survey Benchmark Glacier Project has recovered firn cores from a site on Wolverine Glacier in Alaska's Kenai Mountains. We use annual horizons in repeat cores to track firn densification and meltwater retention over seasonal and interannual timescales, and we use density measurements to quantify how the firn air content (FAC) changes through time. The results suggest the firn is densifying due primarily to compaction rather than refreezing. Liquid-water retention in the firn is transient, likely due to gravity-fed drainage and irreducible-water-content decreases that accompany decreasing porosity. We show that the uncertainty (±60 kg m −3 ) in the commonly used volume-to-mass conversion factor of 850 kg m −3 is an underestimation when glacier-wide FAC variability exceeds 12% of the glacier-averaged height change. Our results demonstrate how direct measurements of firn properties on mountain glaciers can be used to better quantify the uncertainty in geodetic volume-to-mass conversions.

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↗
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↗
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↗
USGS2017

Mapping tide-water glacier dynamics in east Greenland using landsat data

Landsat multispectral scanner and thematic mapper images were co-registered For the Kangerdlugssuaq Fjord region in East Greenland and were used to map glacier drainage-basin areas, changes in the positions of tide-water glacier termini and to estimate surface velocities of the larger tide-water glaciers. Statistics were compiled to document distance and area changes to glacier termini. The methodologies developed in this study are broadly applicable to the investigation of tide-water glaciers in other areas. The number of images available for consecutive years and the accuracy with which images are co-registered are key factors that influence the degree to which regional glacier dynamics can be characterized using remotely sensed data. Three domains of glacier state were interpreted: net increase in terminus area in the southern part of the study area, net loss of terminus area for glaciers in upper Kangerdlugssuaq Fjord and a slight loss of glacier terminus area northward from Ryberg Fjord. Local increases in the concentrations of drifting icebergs in the fjords coincide with the observed extension of glacier termini positions Ice-surface velocity estimates were derived for several glaciers using automated image cross-correlation techniques The velocity determined for Kangerdlugssuaq Gletscher is approximately 5.0 km a −1 and that for Kong Christian IV Gletscher is 0.9 km a −1 . The continuous presence of icebergs and brash ice in front of these glaciers indicates sustained rates of ice-front calving.

Journal of Glaciology↗
USGS2017

A laboratory study of particle ploughing and pore-pressure feedback: A velocity-weakening mechanism for soft glacier beds

If basal-water discharge and pressure are sufficiently high, a soft-bedded glacier will slip over its bed by ploughing, the process in which particles that span the ice-bed interface are dragged across the bed surface. Results of laboratory experiments indicate that resistance to ploughing can decrease with increasing ploughing velocity (velocity weakening). During ploughing at various velocities (15-400 ma-1), till was compacted in front of idealized particles, causing pore pressures there that were orders of magnitude higher than the ambient value. This excess pore pressure locally weakened the till in shear, thereby decreasing ploughing resistance by a factor of 3.0-6.6 with a six-fold increase in ploughing velocity. Characteristic timescales of pore-pressure diffusion and compaction down-glacier from ploughing particles depend on till diffusivity, ploughing velocity and sizes of ploughing particles. These timescales accurately predict the ranges of these variables over which excess pore pressure and velocity weakening occurred. Existing ploughing models do not account for velocity weakening. A new ploughing model with no adjustable parameters predicts ploughing resistance to no worse than 38% but requires that excess pore pressures be measured. Velocity weakening by this mechanism may affect fast glacier flow, sediment transport by bed deformation and basal seismicity.

Journal of Glaciology↗