Ice-Bed Coupling Beneath and Beyond Ice Streams: Byrd Glacier, Antarctica

Ice sheet thickness is determined mainly by the strength of ice-bed coupling that controls holistic transitions from slow sheet flow to fast streamflow to buttressing shelf flow. Byrd Glacier has the largest ice drainage system in Antarctica and is the fastest ice stream entering Ross Ice Shelf. In 2004 two large subglacial lakes at the head of Byrd Glacier suddenly drained and increased the terminal ice velocity of Byrd Glacier from 820 m yr(-1) to 900 m yr(-1). This resulted in partial ice-bed recoupling above the lakes and partial decoupling along Byrd Glacier. An attempt to quantify this behavior is made using flowband and flowline models in which the controlling variable for ice height above the bed is the floating fraction phi of ice along the flowband and flowline. Changes in phi before and after drainage are obtained from available data, but more reliable data in the map plane are required before Byrd Glacier can be modeled adequately. A holistic sliding velocity is derived that depends on phi, with contributions from ice shearing over coupled beds and ice stretching over uncoupled beds, as is done in state-of-the-art sliding theories.

Step-Wise Changes in Glacier Flow Speed Coincide With Calving and Glacial Earthquakes at Helheim Glacier, Greenland

Geodetic observations show several large, sudden increases in flow speed at Helheim Glacier, one of Greenland's largest outlet glaciers, during summer, 2007. These step-like accelerations, detected along the length of the glacier, coincide with teleseismically detected glacial earthquakes and major iceberg calving events. No coseismic offset in the position of the glacier surface is observed; instead, modest tsunamis associated with the glacial earthquakes implicate glacier calving in the seismogenic process. Our results link changes in glacier velocity directly to calving-front behavior at Greenland's largest outlet glaciers, on timescales as short as minutes to hours, and clarify the mechanism by which glacial earthquakes occur. Citation: Nettles, M., et al. (2008), Step-wise changes in glacier flow speed coincide with calving and glacial earthquakes at Helheim Glacier, Greenland, Geophys. Res. Lett., 35, L24503, doi: 10.1029/2008GL036127.

Oxygen Isotopic and Soluble Ionic Composition of a Shallow Firn Core, Inilchek Glacier, Central Tien Shan

Oxygen isotopic and soluble ionic measurements made on snow-pit (2 in depth) and firn-core (12.4 m depth samples recovered from the accumulation zone 5100 m) of Inilchek glacier 43degrees N, 79degrees E) provide information on recent (1992-98) climatic and environmental conditions in the central Tien Shan region of central Asia. The combined 14.4 m snow-pit/firn-core profile lies within the firn zone, arid contains only one observed melt feature (10 m temperature = - 12 degreesC), Although some post-depositional attenuation of the sub-seasonal delta(18)O record is possible, annual cycles are apparent throughout the isotope profile. We therefore use the preserved delta(18)O record to establish a depth/age scale for the core. Mean delta(18)O values for the entire core and for summer periods are consistent with delta(18)O/temperature observations, and suggest the delta(18)O record provides a means to reconstruct past changes in summer surface temperature at the site. Major-ion (Na(+), K(+), Mg(2+), Ca(2+), NH(4)(+), Cl(-), NO(3)(-), SO(4)(2-)) data from the core demonstrate the dominant influence of dust deposition on the soluble chemistry at the site, arid indicate significant interannual variability in atmospheric-dust loading during the 1900s. Anthropogenic impacts oil NH(4)(+) concentrations are observed at the site, and suggest a summer increase in atmospheric NH(4)(+) that may be related to regional agricultural (nitrogen-rich fertilizer use activities.

Surface "Waves" on Byrd Glacier, Antarctica

Byrd Glacier has one of the largest ice catchment areas in Antarctica, delivers more ice to the Ross Ice Shelf than any other ice stream, and is the fastest of these ice streams. A force balance, combined with a mass balance, demonstrates that stream flow in Byrd Glacier is transitional from sheet flow in East Antarctica to shelf flow in the Ross Ice Shelf. The longitudinal pulling stress, calculated along an ice flowband from the force balance, is linked to variations of ice thickness, to the ratio of the basal water pressure to the ice overburden pressure where Byrd Glacier is grounded, and is reduced by an ice-shelf buttressing stress where Byrd Glacier is floating. Longitudinal tension peaks at pressure-ratio maxima in grounded ice and close to minima in the ratio of the pulling stress to the buttressing stress in floating ice. The longitudinal spacing of these tension peaks is rather uniform and, for grounded ice, the peaks occur at maxima in surface slope that have no clear relation to the bed slope. This implies that the maxima in surface slope constitute a "wave train" that is related to regular variations in ice-bed coupling, not primarily to bed topography. It is unclear whether these surface "waves" are "standing waves" or are migrating either upslope or downslope, possibly causing the grounding line to either retreat or advance. Deciding which is the case will require obtaining bed topography in the map plane, a new map of surface topography, and more sophisticated modeling that includes ice flow linked to subglacial hydrology in the map plane.

Basal Melting Along the Floating Part of Byrd Glacier

A mass balance calculation was made for the floating part of Byrd Glacier, using 1978-79 ice elevation and velocity data, over the 45 km of Byrd Glacier from its grounding line to where it leaves its fjord and merges with the Ross Ice Shelf. Smoothed basal melting rates were relatively uniform over this distance and averaged 12 +/- 2 m yr(-1).

A New Velocity Map for Byrd Glacier, East Antarctica, from Sequential Aster Satellite Imagery

New ice-velocity measurements are obtained for the main trunk of Byrd Glacier, East Antarctica, using recently acquired Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. The velocities are derived from the application of a cross-correlation technique to sequential images acquired in 2000 and 2001. Images were co-registered and ortho-rectified with the aid of a digital elevation model (DEM) generated from ASTER stereo imagery. This paper outlines the process of DEM generation, image co-registration and correction, and the application of the cross-correlation technique to obtain ice velocities. Comparison of the new velocity map with earlier measurements of velocity from 1978 indicates that the glacier has undergone a substantial deceleration between observations. Portions of the glacier flowing at speeds of similar to 850 m a(-1) in 1978/79 were flowing at similar to 650 m a(-1) in 2000/01. The cause of this change in ice dynamics is not known, but the observation shows that East Antarctic outlet glaciers can undergo substantial changes on relatively short timescales.

Preliminary Results from the Chemical Records of an 80.4 M Ice Core Recovered from East Rongbuk Glacier, Qomolangma (Mount Everest), Himalaya

High-resolution chemical records from an 80.4 m ice core from the central Himalaya demonstrate climatic and environmental changes since 1844. The chronological net accumulation series shows a sharp decrease from the mid-1950s, which is coincident with the widely observed glacier retreat. A negative correlation is found between the ice-core delta(18)O record and the monsoon precipitation for Indian region 7. The temporal variation of the terrestrial ions (Ca2+ and Mg2+) is controlled by both the monsoon precipitation for Indian regions 3,7 and 8, located directly south and west of the Himalaya, and the dust-storm duration and frequency in the northern arid regions, such as the Taklimakan desert, China. The NH4+ profile is fairly flat until the 1940s, then substantially increases until the end of the 1980s, with a slight decrease during the 1990s which may reflect new agricultural practices. The SO42- and NO3- profiles show an apparent increasing trend, especially during the period 1940s-80s. Moreover, SO42- concentrations for the East Rongbuk Glacier core are roughly double that of the nearby Dasuopu core at Xixabangma, Himalaya, due to local human activity including that of climbing teams who use gasoline for cooking, energy and transport.

Evidence for a Frozen Bed, Byrd Glacier, Antarctica

Ice thickness, computed within the fjord region of Byrd Glacier on the assumptions that Byrd Glacier is in mass-balance equilibrium and that ice velocity is entirely due to basal sliding, are on average 400 m less than measured ice thicknesses along a radio-echo profile. We consider four explanations for these differences: (1) active glacier ice is separated from a zone of stagnant ice near the base of the glacier by a shear zone at depth; (2) basal melting rates are some 8 m/yr; (3) internal shear occurs with no basal sliding in much of the region above the grounding zone; or (4) internal creep and basal sliding contribute to the flow velocity in varying proportions above the grounding zone. Large gradients of surface strain rate seem to invalidate the first explanation. Computed values of basal shear stress (140 to 200 kPa) provide insufficient frictional heat to melt the ice demanded by the second explanation. Both the third and fourth explanations were examined by making simplifying assumptions that prevented a truly quantitative evaluation of their merit. Nevertheless, there is no escaping the qualitative conclusion that internal shear contributes strongly to surface velocities measured on Byrd Glacier, as is postulated in both these explanations.

Solar Forcing Recorded by Aerosol Concentrations in Coastal Antarctic Glacier Ice, McMurdo Dry Valleys

Ice-core chemistry data from Victoria Lower Glacier, Antarctica, suggest, at least for the last 50 years, a direct influence of solar activity variations on the McMurdo Dry Valleys (MDV) climate system via controls on air-mass input from two competing environments: the East Antarctic ice sheet and the Ross Sea. During periods of increased solar activity, when total solar irradiance is relatively high, the MDV climate system appears to be dominated by air masses originating from the Ross Sea, leading to higher aerosol deposition. During reduced solar activity, the Antarctic interior seems to be the dominant air-mass source, leading to lower aerosol concentration in the ice-core record. We propose that the sensitivity of the MDV to variations in solar irradiance is caused by strong albedo differences between the ice-free MDV and the ice sheet.

Tidal Flexure of Jakobshavns Glacier, West Greenland

Jakobshavns Glacier, a floating outlet glacier on the West Greenland coast, was surveyed during July 1976. The vertical displacements of targets along two profiles perpendicular to the fjord wall bounding the north margin of the glacier were analyzed to determine the effect of flexure caused by tidal oscillations within the fjord. An analysis based on the assumption that vertical displacements of the glacier reflected pure elastic bending yielded the conclusion that the effective thickness of the ice (i.e., the thickness which remained unaffected by surface and basal cracking and which behaved as a continuum) was ∼160 m 2.6 km upglacier from the calving front and ∼110 m 0.6 km from the calving front. An analysis based on the more realistic assumption that observed bending reflected elastic and viscoplastic deformation yielded the conclusion that the average effective thickness of the ice was 316 ± 74 m (∼40% of the estimated 800-m total thickness) 2.6 km from the calving front and 160 ± 48 m (∼21% of the estimated 750-m total) 0.6 km from the calving front. A constitutive relationship appropriate for hard glide during flexure was used.

Glacier Calving: A Numerical Model of Forces in the Calving-Speed/Water-Depth Relation

Empirical data suggest that the race of calving of grounded glaciers terminating in water is directly proportional to the water depth. Important controls on calving may be the extent to which a calving face tends to become oversteepened by differential flow within the ice and the extent to which bending moments promote extrusion and bottom crevassing at the base of a calving face. Numerical modelling suggests that the tendency to become oversteepened increases roughly linearly with water depth. In addition, extending longitudinal deviatoric stresses at the base of a calving face increase with water depth. These processes provide a possible physical explanation for the observed calving-rate/water-depth relation.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.