A 12,000 Year Record of Explosive Volcanism in the Siple Dome Ice Core, West Antarctica

Air mass trajectories in the Southern Hemisphere provide a mechanism for transport to and deposition of volcanic products on the Antarctic ice sheet from local volcanoes and from tropical and subtropical volcanic centers. This study extends the detailed record of Antarctic, South American, and equatorial volcanism over the last 12,000 years using continuous glaciochemical series developed from the Siple Dome A (SDMA) ice core, West Antarctica. The largest volcanic sulfate spike ( 280 mu g/L) occurs at 5881 B. C. E. Other large signals with unknown sources are observed around 325 B. C. E. ( 270 mu g/L) and 2818 B. C. E. ( 191 mu g/L). Ages of several large equatorial or Southern Hemisphere volcanic eruptions are synchronous with many sulfate peaks detected in the SDMA volcanic ice chemistry record. The microprobe "fingerprinting'' of glass shards in the SDMA core points to the following Antarctic volcanic centers as sources of tephra found in the SDMA core: Balenny Island, Pleiades, Mount Berlin, Mount Takahe, and Mount Melbourne as well as Mount Hudson and possibly Mount Burney volcanoes of South America. Identified volcanic sources provide an insight into the poorly resolved transport history of volcanic products from source volcanoes to the West Antarctic ice sheet.

Variations of Ice Bed Coupling Beneath and Beyond Ice Streams: The Force Balance

A geometrical force balance that links stresses to ice bed coupling along a flow band of an ice sheet was developed in 1988 for longitudinal tension in ice streams and published 4 years later. It remains a work in progress. Now gravitational forces balanced by forces producing tensile, compressive, basal shear, and side shear stresses are all linked to ice bed coupling by the floating fraction phi of ice that produces the concave surface of ice streams. These lead inexorably to a simple formula showing how phi varies along these flow bands where surface and bed topography are known: phi = h(O)/h(I) with h(O) being ice thickness h(I) at x = 0 for x horizontal and positive upslope from grounded ice margins. This captures the basic fact in glaciology: the height of ice depends on how strongly ice couples to the bed. It shows how far a high convex ice sheet (phi = 0) has gone in collapsing into a low flat ice shelf (phi = 1). Here phi captures ice bed coupling under an ice stream and h(O) captures ice bed coupling beyond ice streams.

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.

Elastic Uplift in Southeast Greenland Due to Rapid Ice Mass Loss

The rapid unloading of ice from the southeastern sector of the Greenland ice sheet between 2001 and 2006 caused an elastic uplift of similar to 35 mm at a GPS site in Kulusuk. Most of the uplift results from ice dynamic-induced volume losses on two nearby outlet glaciers. Volume loss from Helheim Glacier, calculated from sequential digital elevation models, contributes about similar to 16 mm of the observed uplift, with an additional similar to 5 mm from volume loss of Kangerdlugssuaq Glacier. The remaining uplift signal is attributed to significant melt-induced ice volume loss from the ice sheet margin along the southeast coast between 62 degrees N and 66 degrees N.

High-Resolution Ice Cores from US ITASE (West Antarctica): Development and Validation of Chronologies and Determination of Precision and Accuracy

Shallow ice cores were obtained from widely distributed sites across the West Antarctic ice sheet, as part of the United States portion of the International Trans-Antarctic Scientific Expedition (US ITASE) program. The US ITASE cores have been dated by annual-layer counting, primarily through the identification of summer peaks in non-sea-salt sulfate (nssSO(4)(2-)) concentration. Absolute dating accuracy of better than 2 years and relative dating accuracy better than 1 year is demonstrated by the identification of multiple volcanic marker horizons in each of the cores, Tambora, Indonesia (1815), being the most prominent. Independent validation is provided by the tracing of isochronal layers from site to site using high-frequency ice-penetrating radar observations, and by the timing of mid-winter warming events in stable-isotope ratios, which demonstrate significantly better than 1 year accuracy in the last 20 years. Dating precision to 1 month is demonstrated by the occurrence of summer nitrate peaks and stable-isotope ratios in phase with nssSO(4)(2-), and winter-time sea-salt peaks out of phase, with phase variation of < 1 month. Dating precision and accuracy are uniform with depth, for at least the last 100 years.

Topographic Control of Regional Accumulation Rate Variability at South Pole and Implications for Ice-Core Interpretation

Snow-accumulation rates are known to be sensitive to local changes in ice-sheet surface slope because of the effect of katabatic winds. These topographic effects can be preserved in ice cores that are collected at non-ice-divide locations. The trajectory of an ice-core site at South Pole is reconstructed using measurements of ice-sheet motion to show that snow was probably deposited at places of different surface slope during the past 1000 years. Recent accumulation rates, derived from shallow firn cores, vary along this trajectory according to surface topography, so that on a relatively steep flank mean annual accumulation is similar to 18% smaller than on a nearby topographic depression. These modern accumulation rates are used to reinterpret the cause of accumulation rate variability with time in the long ice-core record as an ice-dynamics effect and not a climate-change signal. The results highlight the importance of conducting ancillary ice-dynamics measurements as part of ice-coring programs so that topographic effects can be deconvolved from potential climate signals.

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.

Thermal Convection in Ice Sheets: We Look But Do Not See

Thermal convection in the Antarctic and Greenland ice sheets has been dismissed on the grounds that radio-echo stratigraphy is undisturbed for long distances. However, the undisturbed stratigraphy lies, for the most part, above the density inversion in polar ice sheets and therefore does not disprove convection. An echo-free zone is widespread below the density inversion, yet nobody has cited this as a strong indication that convection is indeed present at d�pth. A generalized Rayleigh criterion for thermal convection in e1astic-viscoplastic polycrystalline solids heated from below is developed and applied to ice-sheet convection. An infinite Rayleigh number at the onset of primary creep decreases with time and becomes constant when secondary creep dominates, suggesting that any thermal buoyancy stress can initiate convection but convection cannot be sustained below a buoyancy stress of about 3 kPa. An analysis of the temperature profile down the Byrd Station core hole suggests that about 1000 m of ice below the density inversion will sustain convection. Creep along the Byrd Station strain network, radar sounding in East Antarctica, and seismic sounding in West Antarctica are examined for evidence of convective creep superimposed on advective creep. It is concluded that the evidence for convection is there, if we look for it with the intention offinding it.

Climate Variability in West Antarctica Derived from Annual Accumulation-Rate Records from ITASE Firn/Ice Cores

Thirteen annually resolved accumulation-rate records covering the last similar to 200 years from the Pine Island-Thwaites and Ross drainage systems and the South Pole are used to examine climate variability over West Antarctica. Accumulation is controlled spatially by the topography of the ice sheet, and temporally by changes in moisture transport and cyclonic activity. A comparison of mean accumulation since 1970 at each site to the long-term mean indicates an increase in accumulation for sites located in the western sector of the Pine Island-Thwaites drainage system. Accumulation is negatively associated with the Southern Oscillation Index (Sol) for sites near the ice divide, and periods of sustained negative Sol (1940-42, 1991-95) correspond to above-mean accumulation at most sites. Correlations of the accumulation-rate records with sea-level pressure (SLP) and the SOI suggest that accumulation near the ice divide and in the Ross drainage system may be associated with the midlatitudes. The post-1970 increase in accumulation coupled with strong SLP-accumulation-rate correlations near the coast suggests recent intensification of cyclonic activity in the Pine Island-Thwaites drainage system.

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.

Modeling Englacial Radar Attenuation at Siple Dome, West Antarctica, Using Ice Chemistry and Temperature Data

The radar reflectivity of an ice-sheet bed is a primary measurement for discriminating between thawed and frozen beds. Uncertainty in englacial radar attenuation and its spatial variation introduces corresponding uncertainty in estimates of basal reflectivity. Radar attenuation is proportional to ice conductivity, which depends on the concentrations of acid and sea-salt chloride and the temperature of the ice. We synthesize published conductivity measurements to specify an ice-conductivity model and find that some of the dielectric properties of ice at radar frequencies are not yet well constrained. Using depth profiles of ice-core chemistry and borehole temperature and an average of the experimental values for the dielectric properties, we calculate an attenuation rate profile for Siple Dome, West Antarctica. The depth-averaged modeled attenuation rate at Siple Dome (20.0 +/- 5.7 dB km(-1)) is somewhat lower than the value derived from radar profiles (25.3 +/- 1.1 dB km(-1)). Pending more experimental data on the dielectric properties of ice, we can match the modeled and radar-derived attenuation rates by an adjustment to the value for the pure ice conductivity that is within the range of reported values. Alternatively, using the pure ice dielectric properties derived from the most extensive single data set, the modeled depth-averaged attenuation rate is 24.0 +/- 2.2 dB km(-1). This work shows how to calculate englacial radar attenuation using ice chemistry and temperature data and establishes a basis for mapping spatial variations in radar attenuation across an ice sheet.

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.

Bending Shear: The Rate-Controlling Mechanism for Calving Ice Walls

Bending shear was observed to produce nearly vertical shear bands in a calving ice wall standing on dry land on Deception Island (Iat. 63.0 oS., long. 60.6 W.), and slabs calved straight downward when shear rupture occurred along these shear bands (Hughes, 1989). A formula for the calving rate was developed from the Deception Island data, and we have attempted to justify generalizing this formula to include ice walls standing along beaches or in water. These are environments in which a wave-washed groove develops along the base of the ice wall or along a water line above the base. The rate of wave erosion provides an alternative mechanism for controlling the calving rate in these environments. We have determined that the rate at which bending creep produces nearly vertical shear bands, along which shear r upture occurs, controls the calving rate in all environments. Shear rupture occurs at a calving shear stress of about I bar. Our results justify using the calving formula to compute the calving rate of ice walls in computer models of ice-sheet dynamics. This is especially important in simulating retreat of Northern Hemisphere ice sheets during the last deglaciation, when marine and lacustrine environments were common along retreating ice margins. These margins would have been ice walls standing along beaches or in water, because floating ice shelves are not expected in the ablation zone of retreating ice sheets.

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.

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.