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.

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.

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.

Mass Balance and Accumulation Rate Across Siple Dome, West Antarctica

Snow-accumulation rates and rates of ice-thickness change (mass balance) are studied at several sites on Siple Dome, West Antarctica. Accumulation rates are derived from analyses of gross beta radioactivity in shallow firn cores located along a 60 km transect spanning both flanks and the crest of the dome. There is a north-south gradient in snow-accumulation rate across the dome that is consistent with earlier radar mapping of internal stratigraphy. Orographic processes probably control this distribution. Mass balance is inferred from the difference between global positioning system (GPS)-derived vertical velocities and snow-accumulation rates for sites close to the firn-core locations. Results indicate that there is virtually no net thickness change at four of the five sites. The exception is at the northernmost site where a small amount of thinning is detected, that appears to be inconsistent with other studies. A possible cause of this anomalous thinning is recent retreat of the grounding line of Ice Stream D.

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.

GPR Reflection Profiles of Clark and Commonwealth Glaciers, Dry Valleys, Antarctica

Englacial horizons deeper than 100 m are absent within 100 MHz ground-penetrating radar (GPR) surface profiles we recorded on Clark and Commonwealth Glaciers in the Antarctic Dry Valleys region. Both glaciers show continuous bottom horizons to 280 m, with bottom signal-to-noise ratios near 30 dB. Density horizons should fade below 50 m depth because impermeable ice occurred by 36 m. Folding within Commonwealth Glacier could preclude radar strata beneath about 80 m depth, but there is no significant folding within Clark Glacier. Strong sulfate concentrations and contrasts exist in our shallow ice core. However, it appears that high background concentration levels, and possible decreased concentration contrasts with depth placed their corresponding reflection coefficients at the limit of, or below, our system sensitivity by about 77 m depth. Further verification of this conclusion awaits processing of our deep-core chemistry profiles.

A 200 Year Sub-Annual Record of Sulfate in West Antarctica, from Sixteen Ice Cores

Sixteen high-resolution ice-core records from West Antarctica and South Pole are used to examine the spatial and temporal distribution of sulfate for the last 200 years. The preservation of seasonal layers throughout the length of each record results in a dating accuracy of better than 1 year based on known global-scale volcanic events. A dual transport source for West Antarctic sea-salt (ss) SO42- and excess (xs) SO42- is observed: lower-tropospheric for areas below 1000m elevation and mid-/upper-tropospheric/stratospheric for areas located above 1000m. Our XsSO(4)(2-) records with volcanic peaks removed do not display any evidence of an anthropogenic impact on West Antarctic SO42- concentrations but do reveal that a major climate transition takes place over West Antarctica at similar to 1940. Global-scale volcanic eruptions appear as significant peaks in the robust-spline residual xsSO(4)(2-) records from sites located above 1000 m elevation but do not appear in the residual records from sites located below 1000 m.

Stratigraphic Continuity in 400 Mhz Short-Pulse Radar Profiles of Firn in West Antarctica

We track dated firn horizons within 400 MHz short-pulse radar profiles to find the continuous extent over which they can be used as historical benchmarks to study past accumulation rates in West Antarctica. The 30-40 cm pulse resolution compares with the accumulation rates of most areas. We tracked a particular set that varied from 30 to 90 m in depth over a distance of 600 km. The main limitations to continuity are fading at depth, pinching associated with accumulation rate differences within hills and valleys, and artificial fading caused by stacking along dips. The latter two may be overcome through multi-kilometer distances by matching the relative amplitude and spacing of several close horizons, along with their pulse forms and phases. Modeling of reflections from thin layers suggests that the - 37 to - 50 dB range of reflectivity and the pulse waveforms we observed are caused by the numerous thin ice layers observed in core stratigraphy. Constructive interference between reflections from these close, high-density layers can explain the maintenance of reflective strength throughout the depth of the firn despite the effects of compaction. The continuity suggests that these layers formed throughout West Antarctica and possibly into East Antarctica as well.

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.

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.

A High-Altitude Snow Chemistry Record from Amundsenisen, Dronning Maud Land, Antarctica

In this paper a detailed record of major ions from a 20 in deep firn core from Amundsenisen, western Dronning Maud Land, Antarctica, is presented. The core was drilled at 75degreesS, 2degrees E (2900 m.a.s.l.) during austral summer 1991/92. The following ions were measured at 3 cm resolution: Na+, Mg2+, Ca2+, Cl-, NO3-, SO42- and CH3SO3H (MSA). The core was dated back to 1865 using a combination of chemical records and volcanic reference horizons. The volcanic eruptions identified in this core are Mount Ngauruhoe, New Zealand (1974-75), Mount Agung, Indonesia (1963), Azul, Argentina (1932). and a broad peak that corresponds in time to Tarawera, New Zealand (1886), Falcon Island, South Shetlands, Southern Ocean (1885), and Krakatau, Indonesia (1883). There are no trends in any of the ion records, but the annual to decadal changes are large. The mean concentrations of the measured ions are in agreement with those from other high-altitude cores from the Antarctic plateau. At this core site there may be a correspondence between peaks in the MSA record and major El Nino-Southern Oscillation events.

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.

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.

A Record of Atmospheric Co2 During the Last 40,000 Years from the Siple Dome, Antarctica Ice Core

We have measured the CO2 concentration of air occluded during the last 40,000 years in the deep Siple Dome A ( hereafter Siple Dome) ice core, Antarctica. The general trend of CO2 concentration from Siple Dome ice follows the temperature inferred from the isotopic composition of the ice and is mostly in agreement with other Antarctic ice core CO2 records. CO2 rose initially at similar to 17.5 kyr B. P. ( thousand years before 1950), decreased slowly during the Antarctic Cold Reversal, rose during the Younger Dryas, fell to a local minimum at around 8 kyr B. P., and rose continuously since then. The CO2 concentration never reached steady state during the Holocene, as also found in the Taylor Dome and EPICA Dome C ( hereafter Dome C) records. During the last glacial termination, a lag of CO2 versus Siple Dome isotopic temperature is probable. The Siple Dome CO2 concentrations during the last glacial termination and in the Holocene are at certain times greater than in other Antarctic ice cores by up to 20 ppm (mumol CO2/mol air). While in situ production of CO2 is one possible cause of the sporadic elevated levels, the mechanism leading to the enrichment is not yet clear.

Monsoonal Circulation of the McMurdo Dry Valleys, Ross Sea Region, Antarctica: Signal from the Snow Chemistry

McMurdo Dry Valleys (MDV, Ross Sea region, Antarctica) precipitation exhibits extreme seasonality in ion concentration, 3-5 orders of magnitude between summer and winter precipitation. To identify aerosol sources and investigate causes for the observed amplitude in concentration variability, four snow pits were sampled along a coast-Polar Plateau transect across the MDV. The elevation of the sites ranges from 50 to 2400 m and the distance from the coast from 8 to 93 km. Average chemistry gradients along the transect indicate that most species have either a predominant marine or terrestrial source in the MDV. Empirical orthogonal function analysis on the snow-chemistry time series shows that at least 57% of aerosol deposition occurs concurrently. A conceptual climate model, based on meteorological observations, is used to explain the strong seasonality in the MDV. Our results suggest that radiative forcing of the ice-free valleys creates a surface low-pressure cell during summer which promotes air-mass flow from the Ross Sea. The associated precipitating air mass is relatively warm, humid and contains a high concentration of aerosols. During winter, the MDV are dominated by air masses draining off the East Antarctic ice sheet, that are characterized by cold, dry and low concentrations of aerosols. The strong differences between these two air-mass sources create in the MDV a polar version of the monsoonal flow, with humid, warm summers and dry, cold winters.