Spatial Patterns in Mass Balance of the Siple Coast and Amundsen Sea Sectors, West Antarctica

Local rates of change in ice-sheet thickness were calculated at IS sites in West Antarctica using the submergence velocity technique. This method entails a comparison of the vertical velocity of the ice sheet, measured using repeat global positioning system surveys of markers, and local long-term rates of snow accumulation obtained using firn-core stratigraphy. Any significant difference between these two quantities represents a thickness change with time. Measurements were conducted at sites located similar to 100-200 km apart along US ITASE traverse routes, and at several isolated locations. All but one of the sites are distributed in the Siple Coast and the Amundsen Sea basin along contours of constant elevation, along flowlines, across ice divides and close to regions of enhanced flow. Calculated rates of thickness change are different from site to site. Most of the large rates of change in ice thickness (similar to 10 cm a(-1) or larger) are observed in or close to regions of rapid flow, and are probably related to ice-dynamics effects. Near-steady-state conditions are calculated mostly at sites in the slow-moving ice-sheet interior and near the main West Antarctic ice divide. These results are consistent with regional estimates of ice-sheet change derived from remote-sensing measurements at similar locations in West Antarctica.

Stratigraphic Variation Within Polar Firn Caused by Differential Accumulation and Ice Flow: Interpretation of a 400 Mhz Short-Pulse Radar Profile from West Antarctica

We investigate causes of the stratigraphic variation revealed in a 177 km, 400 MHz short-pulse radar profile of firn from West Antarctica. The profile covers 56 m depth, and its direction was close to those of the ice flow and mean wind. The average, near-surface accumulation rates calculated from the time delays of one radar horizon consistently show minima on leeward slopes and maxima on windward slopes, confirming an earlier study based on stake observations. The stratigraphic variation includes up to 30 m depth variation in individual horizons over tens of km, fold limbs that become progressively steeper with depth, and fold-hinge loci that change direction or propagate down-ice with depth over distances far less than predicted by the ice speeds. We use an accumulation rate model to show how local rate anomalies and the effect of ice speed upon a periodic variation in accumulation rate cause these phenomena, and we reproduce two key features seen in the stratigraphic variations. We conclude that the model provides an explanation of changes in spatial stratigraphy and local measures of accumulation history given the constraints of surface topography, ice and wind velocities, and a general accumulation rate for an area.

Dating the Siple Dome (Antarctica) Ice Core By Manual and Computer Interpretation of Annual Layering

The Holocene portion of the Siple Dome (Antarctica) ice core was dated by interpreting the electrical, visual and chemical properties of the core. The data were interpreted manually and with a computer algorithm. The algorithm interpretation was adjusted to be consistent with atmospheric methane stratigraphic ties to the GISP2 (Greenland Ice Sheet Project 2) ice core, (BE)-B-10 stratigraphic ties to the dendrochronology C-14 record and the dated volcanic stratigraphy. The algorithm interpretation is more consistent and better quantified than the tedious and subjective manual interpretation.

CO(2) Diffusion in Polar Ice: Observations from Naturally Formed CO(2) Spikes in the Siple Dome (Antarctica) Ice Core

One common assumption in interpreting ice-core CO(2) records is that diffusion in the ice does not affect the concentration profile. However, this assumption remains untested because the extremely small CO(2) diffusion coefficient in ice has not been accurately determined in the laboratory. In this study we take advantage of high levels of CO(2) associated with refrozen layers in an ice core from Siple Dome, Antarctica, to study CO(2) diffusion rates. We use noble gases (Xe/Ar and Kr/Ar), electrical conductivity and Ca(2+) ion concentrations to show that substantial CO(2) diffusion may occur in ice on timescales of thousands of years. We estimate the permeation coefficient for CO(2) in ice is similar to 4 x 10(-21) mol m(-1) s(-1) Pa(-1) at -23 degrees C in the top 287 m (corresponding to 2.74 kyr). Smoothing of the CO(2) record by diffusion at this depth/age is one or two orders of magnitude smaller than the smoothing in the firn. However, simulations for depths of similar to 930-950m (similar to 60-70 kyr) indicate that smoothing of the CO(2) record by diffusion in deep ice is comparable to smoothing in the firn. Other types of diffusion (e.g. via liquid in ice grain boundaries or veins) may also be important but their influence has not been quantified.

A Coastal Transect of McMurdo Dry Valleys (Antarctica) Snow and Firn: Marine and Terrestrial Influences on Glaciochemistry

Samples of snow and firn from accumulation zones on Clark, Commonwealth, Blue and Victoria Upper Glaciers in the McMurdo Dry Valleys (similar to 77-78 degrees S, 161-164 degrees E), Antarctica, are evaluated chemically and isotopically to determine the relative importance of local (site-specific) factors vs regional-scale influences in defining glaciochemistry. Spatial variation in snow and firn chemistry confirms documented trends within individual valleys regarding major-ion deposition relative to elevation and to distance from the coast. Sodium and methylsulfonate (MS-), for example, follow a decreasing gradient with distance from the coast along the axis of Victoria Valley (350-119 mu gL(-1) for Na+; 33-14 mu gL(-1) for MS-); a similar pattern exists between Commonwealth and Newall Glaciers in the Asgaard Range. When comparing major-ion concentrations (e.g. Na-+,Na- MS-, Ca2+) or trace metals (e.g. Al, Fe) among different valleys, however, site-specific exposures to marine and local terrestrial chemical sources play a dominant role. Because chemical signals at all sites respond to particulates with varying mixtures of marine and terrestrial sources, each of these influences on site glaciochemistry must be considered when drawing temporal climate inferences on regional scales.

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.