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.

Variability in Accumulation Rates from GPR Profiling on the West Antarctic Plateau

lsochronal layers in firn detected with ground-penetrating radar (GPR) and dated using results from ice-core analyses are used to calculate accumulation rates along a 100 km across-flow profile in West Antarctica. Accumulation rates are shown to be highly variable over short distances. Elevation measurements from global positioning system surveys show that accumulation rates derived from shallow horizons correlate well with surface undulations, which implies that wind redistribution of snow is the leading cause of this variability. Temporal changes in accumulation rate over 25-185 year intervals are smoothed to along-track length scales comparable to surface undulations in order to identify trends in accumulation that are likely related to changes in climate. Results show that accumulation rates along this profile have decreased in recent decades, which is consistent with core-derived time series of annual accumulation rates measured at the two ends of the radar profile. These results suggest that temporal variability observed in accumulation-rate records from ice cores and GPR profiles can be obscured by spatial influences, although it is possible to resolve temporal signals if the effects of local topography and ice flow are quantified and removed.

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.

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.

Sulfur Isotopic Measurements from a West Antarctic Ice Core: Implications for Sulfate Source and Transport

Measurements of delta(34)S covering the years 1935-76 and including the 1963 Agung (Indonesia) eruption were made on a West Antarctic firn core, RIDSA (78.73 degrees S, 116.33 degrees W; 1740m a.s.l.), and results are used to unravel potential source functions in the sulfur cycle over West Antarctica. The delta(34)S values Of SO42- range from 3.1 parts per thousand to 9.9 parts per thousand. These values are lower than those reported for central Antarctica, from near South Pole station, of 9.3-18.1 parts per thousand (Patris and others, 2000). While the Agung period is isotopically distinct at South Pole, it is not in the RIDSA dataset, suggesting differences in the source associations for the sulfur cycle between these two regions. Given the relatively large input of marine aerosols at RIDSA (determined from Na+ data and the seasonal SO42- cycle), there is likely a large marine biogenic SO42- influence. The delta(34)S values indicate, however, that this marine biogenic SO42-, with a well-established delta(34)S of 18 parts per thousand, is mixing with SO42- that has extremely negative delta(34)S values to produce the measured isotope values in the RIDSA core. We suggest that the transport and deposition of stratospheric SO42- in West Antarctica, combined with local volcanic input, accounts for the observed variance in delta(34)S values.

Ice Elevation Near the West Antarctic Ice Sheet Divide During the Last Glaciation

Interior ice elevations of the West Antarctic Ice Sheet (WAIS) during the last glaciation, which can serve as benchmarks for ice-sheet models, are largely unconstrained. Here we report past ice elevation data from the Ohio Range, located near the WAIS divide and the onset region of the Mercer Ice Stream. Cosmogenic exposure ages of glacial erratics that record a WAIS highstand similar to 125 m above the present surface date to similar to 11.5 ka. The deglacial chronology prohibits an interior WAIS contribution to meltwater pulse 1A. Our observational data of ice elevation changes compare well with predictions of a thermomechanical ice-sheet model that incorporates very low basal shear stress downstream of the present day grounding line. We conclude that ice streams in the Ross Sea Embayment had thin, low-slope profiles during the last glaciation and interior WAIS ice elevations during this period were several hundred meters lower than previous reconstructions.

Boron-Bearing Kornerupine from Fiskenaesset, West Greenland: A Reexamination of Specimens from the Type Locality

In 1884, Lorenzen proposed the formula MgAI2SiO6 for his new mineral kornerupine from Fiskenæsset and did not suspect it to contain boron. Lacroix and de Gramont (1919) reported boron in Fiskenæsset kornerupine, while Herd (1973) found none. New analyses (ion microprobe mass analyser and spectrophotometric) of kornerupine in three specimens from the type locality, including the specimens analysed by Lorenzen and Herd, indicate the presence of boron in all three, in amounts ranging from 0.50 to 1.44 wt.% B203, e.g. (Li0.04 Na0.01 Ca0.01) (Mg3.49 Mn0.01 Fe0.17 Ti0.01 Al5.64)Σ9.30 (Si3.67 Al1.02 B0.31)Σ5 O21 (OH0.99 F0.01) for Lorenzen's specimen. Textures and chemical compositions suggest that kornerupine crystallized in equilibrium in the following assemblages, all with anorthite (An 92-95) and phlogopite (XFe = atomic Fe/(Fe + Mg) = 0.028-0.035): (1) kornerupine (0.045)-gedrite (0.067); (2) kornerupine (0.038-0.050)-sapphirine (0.032-0.035); and (3) kornerupine (0.050)-hornblende. Fluorine contents of kornerupine range from 0.01 to 0.06%, of phlogopite, from 0.09 to 0.10%. In the first assemblage, sapphirine (0.040) and corundum are enclosed in radiating bundles of kornerupine; additionally sapphirine, corundum, and/or gedrite occur with chlorite and pinite (cordierite?) as breakdown products of kornerupine. Kornerupine may have formed by reactions such as: gedrite + sapphirine + corundum + B203 (in solution) + H20 = kornerupine + anorthite + Na-phlogopite under conditions of the granulite facies. Boron for kornerupine formation was most likely remobilized by hydrous fluids from metasedimentary rocks occurring along the upper contact of the Fiskenæsset gabbro-anorthosite complex with amphibolite.

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.

Relating Crevassing to Non-Linear Strain in the Floating Part of Jakobshavn Isbrae, West Greenland

Jakobshavn Isbrae is a major ice stream that drains the west-central Greenland ice sheet and becomes afloat in Jakobshavn Isfiord (69degreesN, 49degreesW), where it has maintained the world's fastest-known sustained velocity and calving rate (7 km a(-1)) for at least four decades. The floating portion is approximately 12 km long and 6 km wide. Surface elevations and motion vectors were determined photogrammetrically for about 500 crevasses on the floating ice, and adjacent grounded ice, using aerial photographs obtained 2 weeks apart in July 1985. Surface strain rates were computed from a mesh of 399 quadrilateral elements having velocity measurements at each corner. It is shown that heavy crevassing of floating ice invalidates the assumptions of linear strain theory that (i) surface strain in the floating ice is homogeneous in both space and time, (ii) the squares and products of strain components are nil, and (iii) first- and second-order rotation components are small compared to strain components. Therefore, strain rates and rotation rates were also computed using non-linear strain theory. The percentage difference between computed linear and non-linear second invariants of strain rate per element were greatest (mostly in the range 40-70%) where crevassing is greatest. Isopleths of strain rate parallel and transverse to flow and elevation isopleths relate crevassing to known and inferred pinning points.