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.

Evolution of Melt Pond Volume on the Surface of the Greenland Ice Sheet

The presence of surface meltwater on ice caps and ice sheets is an important glaciological and climatological characteristic. We describe an algorithm for estimating the depth and hence volume of surface melt ponds using multispectral ASTER satellite imagery. The method relies on reasonable assumptions about the albedo of the bottom surface of the ponds and the optical attenuation characteristics of the ponded meltwater. We apply the technique to sequences of satellite imagery acquired over the western margin of the Greenland Ice Sheet to derive changes in melt pond extent and volume during the period 2001 - 2004. Results show large intra- and interannual changes in ponded water volumes, and large volumes of liquid water stored in extensive slush zones.

Calculating Basal Thermal Zones Beneath the Antarctic Ice Sheet

A procedure is presented for using a simple flowline model to calculate the fraction of the bed that is thawed beneath present-day ice sheets, and therefore for mapping thawed, frozen, melting and freezing basal thermal zones. The procedure is based on the proposition, easily demonstrated, that variations in surface slope along ice flowlines are due primarily to variations in bed topography and ice-bed coupling, where ice-bed coupling for sheet flow is represented by the basal thawed fraction. This procedure is then applied to the central flowlines of flow bands on the Antarctic ice sheet where accumulation rates, surface elevations and bed topography are mapped with sufficient accuracy, and where sheet flow rather than stream flow prevails. In East Antarctica, the usual condition is a low thawed fraction in subglacial highlands, but a high thawed fraction in subglacial basins and where ice converges on ice streams. This is consistent with a greater depression of the basal melting temperature and a slower rate of conducting basal heat to the surface where ice is thick, and greater basal frictional heat production where ice flow is fast, as expected for steady-state flow. This correlation is reduced or even reversed where steady-state flow has been disrupted recently, notably where ice-stream surges produced the Dibble and Dalton Iceberg Tongues, both of which are now stagnating. In West Antarctica, for ice draining into the Pine Island Bay polynya of the Amundsen Sea, the basal thawed fraction is consistent with a prolonged and ongoing surge of Pine Island Glacier and with a recently initiated surge of Thwaites Glacier. For ice draining into the Ross Ice Shelf, long ice streams extend nearly to the West Antarctic ice divide. Over the rugged bed topography near the ice divide, no correlation consistent with steady-state sheet flow exists between ice thickness and the basal thawed fraction. The bed is wholly thawed beneath ice streams, even where stream flow is slow. This is consistent with ongoing gravitational collapse of ice entering the Ross Sea embayment and with unstable flow in the ice streams.

Solar Forcing of the Polar Atmosphere

We present highly resolved, annually dated, calibrated proxies for atmospheric circulation from several Antarctic ice cores (ITASE (International Trans-Antarctic Scientific Expedition), Siple Dome, Law Dome) that reveal decadal-scale associations with a South Pole ice-core Be-10 proxy for solar variability over the last 600 years and annual-scale associations with solar variability since AD 1720. We show that increased (decreased) solar irradiance is associated with increased (decreased) zonal wind strength near the edge of the Antarctic polar vortex. The association is particularly strong in the Indian and Pacific Oceans and as such may contribute to understanding climate forcing that controls drought in Australia and other Southern Hemisphere climate events. We also include evidence suggestive of solar forcing of atmospheric circulation near the edge of the Arctic polar vortex based on ice-core records from Mount Logan, Yukon Territory, Canada, and both central and south Greenland as enticement for future investigations. Our identification of solar forcing of the polar atmosphere and its impact on lower latitudes offers a mechanism for better understanding modern climate variability and potentially the initiation of abrupt climate-change events that operate on decadal and faster scales.

On the Pulling Power of Ice Streams

Gravity wants to pull an ice sheet to the center of the Earth, but cannot because the Earth's crust is in the way, so ice is pushed out sideways instead. Or is it? The ice sheet "sees" nothing preventing it from spreading out except air, which is much less massive than ice. Therefore, does not ice rush forward to fill this relative vacuum; does not the relative vacuum suck ice into it, because Nature abhors a vacuum? If so, the ice sheet is not only pulled downward by gravity, it is also pulled outward by the relative vacuum. This pulling outward will be most rapid where the ice sheet encounters least resistance. The least resistance exists along the bed of ice streams, where ice-bed coupling is reduced by a basal water layer, especially if the ice stream becomes afloat and the floating part is relatively unconfined around its perimeter and unpinned to the sea floor. Ice streams are therefore fast currents of ice that develop near the margins of an ice sheet where these conditions exist. Because of these conditions, ice streams pull ice out of ice sheets and have pulling power equal to the longitudinal gravitational pulling force multiplied by the ice-stream velocity. These boundary conditions beneath and beyond ice streams can be quantified by a basal buoyancy factor that provides a life-cycle classification of ice streams into inception, growth, mature, declining and terminal stages, during which ice streams disintegrate the ice sheet. Surface profiles of ice streams are diagnostic of the stage in a life cycle and, hence, of the vitality of the ice sheet.

ICESat Measurements Reveal Complex Pattern of Elevation Changes on Siple Coast Ice Streams, Antarctica

We compare ICESat data (2003-2004) to airborne laser altimetry data (1997-98 and 1999-2000) to monitor surface changes over portions of Van der Veen (VdVIS), Whillans (WIS) and Kamb ice streams (KIS) in the Ross Embayment of the West Antarctic Ice Sheet. The spatial pattern of detected surface changes is generally consistent with earlier observations. However, important changes have occurred during the past decade. For example, areas on the VdVIS and WIS, where large thinning was detected by the airborne surveys, are now closer to being in balance. The upper trunk of KIS continues to build up with thickening rates reaching 0.4 m/year. Our results provide new evidence that the overall mass balance of the region is becoming more positive, but a significant spatial variability exists. They also demonstrate the potential of ICESat data for detecting spatial patterns of surface elevation change in Antarctica.

Fracture and Back Stress Along the Byrd Glacier Flowband on the Ross Ice Shelf

East Antarctic ice discharged by Byrd Glacier continues as a flowband to the calving front of the Ross Ice Shelf. Flow across the grounding line changes from compressive to extensive as it leaves the fjord through the Transantarctic Mountains occupied by Byrd Glacier. Magnitudes of the longitudinal compressive stress that suppress opening of transverse tensile cracks are calculated for the flowband. As compressive back stresses diminish, initial depths and subsequent growth of these cracks, and their spacing, are calculated using theories of elastic and ductile fracture mechanics. Cracks are initially about one millimeter wide, with approximately 30 in depths and 20 in spacings for a back stress of 83 kPa at a distance of 50 kin beyond the fjord, where floating ice is 600 in thick. When these crevasses penetrate the whole ice thickness, they release tabular icebergs 20 kin to 100 kin wide, spaced parallel to the calving front of the Ross Ice Shelf

Multi-Decadal Record of Ice Dynamics on Daugaard Jensen Gletscher, East Greenland, from Satellite Imagery and Terrestrial Measurements

The history of ice velocity and calving front position of Daugaard Jensen Gletscher, a large outlet glacier in East Greenland, is reconstructed from field measurements, aerial photography and satellite imagery for the period 1950-2001. The calving terminus of the glacier has remained in approximately the same position over the past similar to 50 years. There is no evidence of a change in ice motion between 1968 and 2001, based on a comparison of velocities derived from terrestrial surveying and feature tracking using sequential satellite images. Estimates of flux near the entrance to the fjord vs snow accumulation in the interior catchment show that Daugaard Jensen Gletscher has a small negative mass balance. This result is consistent with other mass-balance estimates for the inland region of the glacier.

Variations in Ice Rafted Detritus on Beaches in the South Shetland Islands: A Possible Climate Proxy

Raised beach ridges on Livingston Island of the South Shetland Islands display variations in both quantity and source of ice rafted detritus (IRD) received over time. Whereas the modem beach exhibits little IRD, all of which is of local origin, the next highest beach (similar to250 C-14 yr BP) has large amounts, some of which comes from as far away as the Antarctic Peninsula. Significant quantities of IRD also were deposited similar to 1750 C-14 yr BP. Both time periods coincide with generally cooler regional conditions and, at least in the case of the similar to250 yr old beach, local glacial advance. We suggest that the increases in ice rafting may reflect periods of greater glacial activity, altered ocean circulation, and/or greater iceberg preservation during the late Holocene. Limited IRD and lack of far-travelled erratics on the modem beach are both consistent with the ongoing warming trend in the Antarctic Peninsula region.

The Effects of Flowline Length Evolution on Chemistry-Delta O-18 Profiles from Penny Ice Cap, Baffin Island, Canada

The isotopic and chemical signatures for ice-age and Holocene ice from Summit, Greenland and Penny Ice Cap, Baffin Island, Canada, arc compared. The usual pattern of low delta(18)O, high Ca2+ and high Cl- is presented in the Summit records, but Penny Ice Cap has lower than present Cl- in its ice-age ice. A simple extension of the Hansson model (Hansson, 1994) is developed and used to simulate these signatures. The low ice-age Cl- from Penny Ice Cap is explained by having the ice-age ice originating many thousands of km inland near the centre of the Laurentide ice sheet and much further from the marine sources. Summit's flowlines all start close to the present site. The Penny Ice Cap early-Holocene delta(18)O's had to be corrected to offset the Laurentide meltwater distortion. The analysis suggests that presently the Summit and Penny Ice Cap marine impurity originates about,500 km away, and that presently Penny Ice Cap receives a significant amount of local continental impurity.

Thickness Changes on Whillans Ice Stream and Ice Stream C, West Antarctica, Derived from Laser Altimeter Measurements

Repeat airborne laser altimeter measurements are used to derive surface elevation changes on parts of Whillans Ice Stream and Ice Stream C, West Antarctica. Elevation changes are converted to estimates of ice equivalent thickness change using local accumulation rates, surface snow densities and vertical bedrock motions. The surveyed portions of two major tributaries of Whillans Ice Stream are found to be thinning almost uniformly at an average rate of similar to 1 m a(-1). Ice Stream C has a complicated elevation-change pattern, but is generally thickening. These results are used to estimate the contribution of each surveyed region to the current rate of global sea-level rise.

Possible Effects of Anthropogenically-Increased CO2 on the Dynamics of Climate: Implications for Ice-Age Cycles

A dynamical model, developed to account for the observed major variations of global ice mass and atmospheric CO2 during the late Cenozoic, is used to provide a quantitative demonstration of the possibility that the anthropogenically-forced increase of atmospheric CO2, if maintained over a long period of time (perhaps by tectonic forcing), could displace the climatic system from an unstable regime of oscillating ice ages into a more stable regime representative of the pre-Pleistocene. This stable regime is characterized by orbitally-forced oscillations that are of much weaker amplitude than prevailed during the Pleistocene.

LGM Summer Climate on the Southern Margin of the Laurentide Ice Sheet: Wet Or Dry?

Regional climate simulations are conducted using the Polar fifth-generation Pennsylvania State University (PSU)-NCAR Mesoscale Model (MM5) with a 60-km horizontal resolution domain over North America to explore the summer climate of the Last Glacial Maximum (LGM: 21 000 calendar years ago), when much of the continent was covered by the Laurentide Ice Sheet (LIS). Output from a tailored NCAR Community Climate Model version 3 (CCM3) simulation of the LGM climate is used to provide the initial and lateral boundary conditions for Polar MM5. LGM boundary conditions include continental ice sheets, appropriate orbital forcing, reduced CO2 concentration, paleovegetation, modified sea surface temperatures, and lowered sea level. The simulated LGM summer climate is characterized by a pronounced low-level thermal gradient along the southern margin of the LIS resulting from the juxtaposition of the cold ice sheet and adjacent warm ice-free land surface. This sharp thermal gradient anchors the midtropospheric jet stream and facilitates the development of synoptic cyclones that track over the ice sheet, some of which produce copious liquid precipitation along and south of the LIS terminus. Precipitation on the southern margin is orographically enhanced as moist southerly low-level flow (resembling a contemporary, Great Plains low-level jet configuration) in advance of the cyclone is drawn up the ice sheet slope. Composites of wet and dry periods on the LIS southern margin illustrate two distinctly different atmospheric flow regimes. Given the episodic nature of the summer rain events, it may be possible to reconcile the model depiction of wet conditions on the LIS southern margin during the LGM summer with the widely accepted interpretation of aridity across the Great Plains based on geological proxy evidence.

On the Disintegration of Ice Shelves: The Role of Fracture

Crevasses can be ignored in studying the dynamics of most glaciers because they are only about 20 m deep, a small fraction of ice thickness. In ice shelves, however, s urface crevasses 20 m deep often reach sealevel and bottom crevasses can move upward to sea-level (Clough, 1974; Weertman, 1980). The ice shelf is fractured completely through if surface and basal crevasses meet (Barrett, 1975; Hughes, 1979). This is especially likely if surface melt water fills surface crevasses (Weertman, 1973; Pfeffer, 1982; Fastook and Schmidt, 1982). Fracture may therefore play an important role i n the disintegration of ice shelves. Two fracture criteria which can be evaluated experimentally and applied to ice shelves, are presented. Fracture is then examined for the general strain field of an ice shelf and for local strain fields caused by shear rupture alongside ice streams entering the ice shelf, fatigue rupture along ice shelf grounding lines, and buckling up-stream from ice rises. The effect of these fracture patterns on the stability of Antarctic ice shelves and the West Antarctic ice sheet is then discussed.

Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.