Search

The mean, trend and variability of net snow accumulation in firn cores are often used to validate model output, develop remote-sensing algorithms and quantify ice-sheet surface mass balance. Thus, accurately defining uncertainties associated with these in situ measurements is critical. In this study, we apply statistical simulation methods to quantify the uncertainty in firn-core accumulation data due to the uncertainty in depth–age scales. The methods are applied to a suite of firn cores from central West Antarctica. The results show that uncertainty in depth–age scales can give rise to spurious trends in accumulation that are the same order of magnitude as accumulation trends reported in West Antarctica. The depth–age scale uncertainties also significantly increase the apparent interannual accumulation variability, so these uncertainties must first be accounted for before using firn-core data to assess such processes as small-spatial-scale variability. Better quantification of error in accumulation will improve our ability to meaningfully compare firn-core data across different regions of the ice sheet, and provide appropriate targets for calibration and/or validation of model output and remote-sensing data.

The Wind River Range in Wyoming, USA, contains the largest concentration of glacial mass in the Rocky Mountains of the contiguous USA. Despite this distinction, only a few field or remotely sensed studies providing glacier volume changes have been published. The current study focuses on Continental Glacier located on the northern end of the range and uses two field datasets (high-accuracy GPS surface elevation points and ice-penetrating radar transects of the glacier bed) to create a three-dimensional model of glacier volume. Current surface elevations are compared with historical elevation data to calculate surface elevation change over time. An average thinning rate of 13.8 ± 7.8 m (0.30 ± 0.17 m a–1) between 1966 and 2012 was found. Surface elevation change rates varied across the glacier, ranging from +0.30 to –0.98 m a–1. Taking into account variable melt rates across the glacier, along with a glacial volume of 72.1 × 106 ± 10.8 × 106 m3, we estimate that Continental Glacier will be reduced in volume by 43% over the next 100 years and will disappear completely over the next 300–400 years, if current climatic conditions persist.

A pulse is a type of unstable glacier flow intermediate between normal flow and surging. Using Landsat MSS, TM and ETM+ imagery and feature-tracking software, a time series of mostly annual velocity maps from 1973 to 2012 was produced that reveals five pulses of Ruth Glacier, Alaska. Peaks in ice velocity were found in 1981, 1989, 1997, 2003 and 2010, approximately every 7 years. During these peak years the ice velocity increased 300%, from approximately 40 m a–1 to 160 m a–1. Based on the spatio-temporal behavior of Ruth Glacier during the pulse cycles, we suggest the pulses are due to enhanced basal motion via deformation of a subglacial till. The cyclical nature of the pulses is interpreted to be due to a thin till, with low permeability, that causes incomplete drainage of the till between the pulses, followed by eventual recharge and dilation of the till. These findings suggest care is needed when attempting to correlate changes in regional climate with decadal-scale changes in velocity, because in some instances basal conditions may have a greater influence on ice dynamics than climate.

Increased mass losses from the Greenland ice sheet and inferred contributions to sea-level rise have heightened the need for hydrologic observations of meltwater exiting the ice sheet. We explore whether temporal variations in ice-sheet surface hydrology can be linked to the development of a downstream sediment plume in Kangerlussuaq Fjord by comparing: (1) plume area and suspended sediment concentration from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery and field data; (2) ice-sheet melt extent from Special Sensor Microwave/Imager (SSM/I) passive microwave data; and (3) supraglacial lake drainage events from MODIS. Results confirm that the origin of the sediment plume is meltwater release from the ice sheet. Interannual variations in plume area reflect interannual variations in surface melting. Plumes appear almost immediately with seasonal surface-melt onset, provided the estuary is free of landfast sea ice. A seasonal hysteresis between melt extent and plume area suggests late-season exhaustion in sediment supply. Analysis of plume sensitivity to supraglacial events is less conclusive, with 69% of melt pulses and 38% of lake drainage events triggering an increase in plume area. We conclude that remote sensing of sediment plume behavior offers a novel tool for detecting the presence, timing and interannual variability of meltwater release from the ice sheet.

This study adjusts and compares digital elevation models (DEMs) created from photogrammetric and interferometric synthetic aperture radar techniques to determine volume and surface elevation changes of five icefields in a remote region of southwest British Columbia, Canada, between the mid-1980s and 1999. Preliminary differences between the DEMs in ice-free and vegetation-free areas indicated variable elevation offsets with increasing altitude (11 m km−1) and with increasing slope (2.7 m (10°)−1). Results indicate a surface elevation change of −6.0 ± 2.7 m (−0.5 ± 0.2 m a−1) and a total volume loss of −19.4 ± 8.8 km3 (−1.5 ± 0.7 km3 a−1), which represents a potential sea-level rise contribution of 0.004 ± 0.002 mm a−1. Temperature and snowfall data from four nearby meteorological stations indicate that increased temperatures and decreased snowfall throughout the late 1980s and 1990s are a likely cause of the thinning. Glacier terminus positions were compared between a historical map (1927) and satellite images (1974, 1990/91 and 2000/01). All observed glaciers were in retreat between 1927 and 1974, as well as between 1990/91 and 2000/01, but many glaciers advanced or significantly slowed in their retreat between 1974 and 1990/91.

Repeat-pass L-band interferometric synthetic aperture radar (InSAR) data for part of Hielo Patagónico Sur, Chile, were collected by the space-shuttle-based Spaceborne Imaging Radar C (SIR-C) over a 4 day span in October 1994. Three co-registered complex SAR images are used to generate phase-coherence maps, a digital elevation model (DEM) and an ice-velocity map. The phase-coherence maps indicate low coherence in the 5–15 km approaching the termini due to large displacements, ice deformation and melting. However, the coherence is high over nearly all of the remaining imaged icefield. Ice-velocity precision is greater than 2 cm d−1, while the DEM is good to about 25 m. A flow divide between two of the glaciers is mapped by locating a narrow band of near-zero ice velocity. Horizontal ice-surface velocity profiles calculated along flowlines show there is a high degree of spatial variability reaching a peak value of 5.5 m d−1 located 3.5 km from the terminus of Glaciar Europa. Longitudinal strain rates along the center lines calculated from these velocities at the locations of the initiation of crevassing are used to compute the tensile strength of ice (169–224 kPa).

The evolution of four dynamic radar glacier zones at the surface of an alpine icefield in British Columbia is monitored using a time series of 35 First European Remote Sensing Satellite (ERS-1) synthetic aperture radar (SAR) images acquired from 1992 to 1994. These zones result from changing wetness and textural properties, and appear to represent: (1) cold snow with no liquid water present; (2) an initial melt front with an upper boundary near the elevation of the 0° isotherm; (3) metamorphosed, rapidly melting first-year snow with a rough or pitted surface; and (4) bare ice. This interpretation is aided by temperature and runoff data, air photographs and field measurements of snowpack properties acquired with two ERS-1 SAR scenes, ice-surface elevations derived from 1:50 000 topographic maps and simulations of radar backscatter from a geometric optics model of surface scattering. Meltwater production is affected by the development of zones (2), (3) and (4), which form, migrate up-elevation and disappear each year between April and September.

The evolution of four dynamic radar glacier zones at the surface of an alpine icefield in British Columbia is monitored using a time series of 35 First European Remote Sensing Satellite (ERS-1) synthetic aperture radar (SAR) images acquired from 1992 to 1994. These zones result from changing wetness and textural properties, and appear to represent: (1) cold snow with no liquid water present; (2) an initial melt front with an upper boundary near the elevation of the 0° isotherm; (3) metamorphosed, rapidly melting first-year snow with a rough or pitted surface; and (4) bare ice. This interpretation is aided by temperature and runoff data, air photographs and field measurements of snowpack properties acquired with two ERS-1 SAR scenes, ice-surface elevations derived from 1:50 000 topographic maps and simulations of radar backscatter from a geometric optics model of surface scattering. Meltwater production is affected by the development of zones (2), (3) and (4), which form, migrate up-elevation and disappear each year between April and September.

Interferometric radar observations of Glaciar San Rafael, Chile, were collected in October 1994 by NASA’s Spaceborne Imaging Radar C (SIR-C) at both L- (24 cm) and C-band frequency (5.6 cm), with vertical transmit and receive polarization. The C-band data did not yield good geophysical products, because the temporal coherence of the signal was significantly reduced alter 24 h. The L-band data were, however, successfully employed to map the surface topography of the icefield with a 10 m uncertainty in height, and measure ice velocity with a precision of 4 mm d−1 or 1.4 m a−1. The corresponding error in strain rates is 0.05 a −1 at a 30 m horizontal spacing. The one-dimensional interferometric velocities were subsequently converted to horizontal displacements by assuming a flow direction and complemented by feature-tracking results near the calving front. The results provide a comprehensive view of the ice-flow dynamics of Glaciar San Rafael. The glacier has a core of rapid flow, 4.5 km in width and 3.5° in average slope, surrounded by slower-moving ice, not by rock. Ice velocity is 2.6 m d−1 or 0.95 km a−1 near the equilibrium-line altitude (1200 m), increasing rapidly before the glacier enters the narrower terminal valley, to reach 17.5 md−1 or 6.4 k ma−1 at the calving front. Strain rates are dominated by lateral shearing at the glacier margins (0.4–0.7 a−1), except for the terminal-valley section, where longitudinal strain rates average close to 1 a−1. This spectacular longitudinal increase in ice velocity in the last few kilometers may be a fundamental feature of tidewater glaciers.

An X-band FM-CW radar was used to determine the feasibility of observing annual snow-accumulation layers in Antarctica with a high-resolution inexpensive radar system. The formation of layering boundaries, their resultant electromagnetic discontinuity and their detection by reflected energy are presented. Large returns from depths corresponding to reasonable positions for annual layers were found. The average accumulation rates calculated from the radar returns agree with those measured in a previous pit study done in the same area. The detection of the annual accumulation layers with this system implies a simple, inexpensive mobile radar could be used to profile large areas allowing the distorting effects of local topography to be removed.This type of system with a concurrent pit study could provide insight into the effect of sub-surface strata on spaceborne or airborne microwave remote sensing.