Observation of thin sections of the WAIS (West Antarctic Ice Sheet) Divide ice core in cross-polarized light reveals a wealth of microstructures and textural characteristics indicative of strain and recovery in an anisotropic crystalline substance undergoing high-temperature plastic deformation. The appearance of abundant subgrain domains-relatively strain-free regions inside crystals (grains) surrounded by walls of dislocations across which small structural orientation changes occur-is particularly noticeable in the depth range associated with the brittle ice (similar to 650-1300 m). Here we describe a subgrain texture, not previously reported in ice, that resembles chessboard-pattern subgrains in beta-quartz. This chessboard texture at WAIS Divide is strongly associated with the presence of bubbles. We hypothesize that chessboard-subgrain development may affect grain-size evolution, the fracture of ice cores recovered from the brittle ice zone and perhaps grain-boundary sliding as well.

We present analyses of bubble number-density (BND) data from the South Pole Ice Core (SPC14) showing warming of similar to 7.5 degrees C from the Late Glacial (similar to 19.5 ka), then relatively stable temperatures during the Holocene (<0.5 degrees C warming), in close agreement with results of independent paleothermometers. The BND data span from similar to 160 m just below pore close-off, to similar to 1200 m, where bubble loss by clathrate formation is significant. Measurements were made with standard bubble 'thick'-section techniques and a new application of three-dimensional micro-computed tomography (CT) imagery; the nearly identical results recommend the faster, nondestructive micro-CT. The very high BND at South Pole, typically 800 and 900 bubbles cm(-3), reflects the joint effects of the relatively low mean-annual temperature (-49 degrees C) and high accumulation rate (similar to 7.5 cm w.e. a(-1)). High BND is physically linked to small grain sizes at pore close-off, which in turn helps explain the near-absence of brittle-ice behavior at the site, contributing to the high quality of the recovered core with implications for siting of future ice cores. The accumulation history, derived from delta N-15-N-2 firn-column thickness estimates, correlates with the temperature history but varies somewhat more than saturation vapor pressure, suggesting dynamic controls including upstream slope variability.

Using samples from the South Pole Ice Core (SPC14), we present new bubble number-density (BND) measurements and a modeled temperature history reconstruction for the South Pole site back through ˜18.5 ka.Additionally, we show that 3D micro-CT sample imagery can accurately quantify BND, enabling more rapid and efficient future analyses.Using sampling and imaging techniques previously established for analyses of the WAIS Divide ice core (Spencer et al., 2006; Fegyveresi et al.,

Automated c-axis analyzers are a critical tool for harvesting large-scale ice crystal orientation data from thin section analysis, but existing examples are not designed for deployment into the field. We demonstrate a possible solution to this need with the Automated Lightweight Portable Analyzer for C-Axes (ALPACA), which implements the established four-sequence measurement algorithm of Wilen (2000) using three motorized axes of motion in a unit with volume about 0.034 cubic meters. In tests of eight easily-visible grains in thin section WDC-06A 420 VTS, ALPACA’s polarizer rotation angles of extinction agreed with previously published data to within 5º on all but one sequence of a single grain, with average sequence errors of 1.6º, 1.3º, 1.4º, and 2.9º. The information produced by these automated measurements enable production of complete grain-by-grain orientation data.

Thwaites Glacier (TG), West Antarctica, is losing mass in response to oceanic forcing. Future evolution could lead to deglaciation of the marine basins of the West Antarctic Ice Sheet, depending on ongoing and future climate forcings, but also on basal topography/bathymetry, basal properties, and physical processes operating within the grounding zone. Hence, it is important to know the distribution of bed types of TG's interior and grounding zone, and to incorporate them accurately in models in order to improve estimates of retreat rates and stability. Here we estimate properties of the bed by determining its acoustic impedance from amplitude analysis of reflection seismic data. We report on the results from two lines – a longitudinal (L-Line) and a transverse (N-Line) – on a central flowline of TG ∼100 km inland from the grounding zone. The data show considerable spatial variability in bed forms and properties, similar to results from a comparable survey farther inland. Notably, we find the same pattern here of hard (presumed bedrock) material on the stoss side of bumps and soft (presumed till) on the lee side. Physical understanding indicates the basal flow law describing motion over different regions of TG's bed likely varies from nearly-viscous over the bedrock regions to nearly-plastic over till regions, providing guidance for modeling.

Measurements of a sample from similar to 580 m depth in the WAIS Divide (WDC06A) ice core reveal that bubbles are preferentially elongated in the basal plane of their parent grain, as expected if bubble shape preserves the record of dominant basal glide. This suggests that a method using bubbles as strain gauges could provide insights to grain-scale ice deformation. We introduce a technique using fabric and image analyses of paired thin and thick sections. Comparison of the crystallo-graphic orientations of 148 grains and the shape orientations of 2377 intragrain bubbles reveals a strongly preferred elongation of bubbles in the grain basal planes (R-2 = 0.96). Elongation magnitudes are consistent with a balance between ice flow deformation and diffusive restoration, with larger bubbles more elongated. Assuming bubbles record ice strain, grains with greater resolved stress on their basal planes from the far-field ice flow stresses show greater deformation, but with large variability suggesting that heterogeneity of the local stress field causes deformation even in unfavorably oriented grains. A correlation is also observed among bubble elongation, grain size, and bubble size, explaining a small but significant fraction of the variance (P < 0.05), with implications for controls on ice deformation, as discussed here.

ABSTRACT We describe methods for measuring crystal orientation fabric with sonic waves in an ice core borehole, with special attention paid to vertical-girdle fabrics that are prevalent at the WAIS Divide. The speed of vertically propagating compressional waves in ice is influenced by vertical clustering of the ice crystal c-axes. Shear-wave speeds – particularly the speed separation between fast and slow shear polarizations – are sensitive to azimuthal anisotropy. Sonic data from the WAIS Divide complement thin-section measurements of fabric. Thin sections show a steady transition to strong girdle fabrics in the upper 2000 m of ice, followed by a transition to vertical-pole fabrics below 2500 m depth. Compressional-wave sonic data are inconclusive in the upper ice, due to noise, as well as the method's inherent insensitivity to girdle fabrics. Compared with available thin sections, sonic data provide better resolution of the transition to pole fabrics below 2500 m, notably including an abrupt increase in vertical clustering near 3000 m. Our compressional-wave measurements resolve fabric changes occurring over depth ranges of a few meters that cannot be inferred from available thin sections, but are sensitive only to zenithal anisotropy. Future logging tools should be designed to measure shear waves in addition to compressional waves, especially for logging in regions where ice flow patterns favor the development of girdle fabrics.

Bubble number densities from the West Antarctic Ice Sheet (WAIS) Divide deep core in West Antarctica record relatively stable temperatures during the middle Holocene followed by late Holocene cooling. We measured bubble number density, shape, size, and arrangement on new samples of the main WAIS Divide deep core WDC06A from similar to 580m to similar to 1600 depth. The bubble size, shape, and arrangement data confirm that the samples satisfy the requirements for temperature reconstructions. A small correction for cracks formed after core recovery allows extension of earlier work through the brittle ice zone, and a site-specific calibration reduces uncertainties. Using an independently constructed accumulation rate history and a steady state bubble number density model, we determined a temperature reconstruction that agrees closely with other independent estimates, showing a stable middle Holocene, followed by a cooling of similar to 1.25 degrees C in the late Holocene. Over the last similar to 5 millennia, accumulation has been higher during warmer times by similar to 12%degrees C-1, somewhat stronger than for thermodynamic control alone, suggesting dynamic processes.

Abstract A method has been devised and tested for measuring the c-axis orientation of crystal grains in thin sections of glacier ice. The crystal orientation and grain size of ice are of great interest to glaciologists since these parameters contain information on the prior thermal and flow history of the ice. The traditional method of determining c-axis orientation involves a transmission measurement through an ice sample, a process that is time-consuming and therefore impractical for obtaining a continuous record. A reflection- or backscatter-based method could potentially be used inside boreholes, with bubbles as reflectors to avoid such drawbacks. The concept demonstration of this paper is performed on ice slices, enabling a direct comparison of accuracy with traditional methods. Measurements of the crystal orientations (θ, ϕ) in 11 grains showed an average error of ±0.8° in ϕ, with no grain error >1.4°. Measurements of θ showed an average error of ±8.2° on ten grains, with unexplained disagreement on the remaining grain. Although the technique is applied specifically to glacier ice, it should be generally applicable to any transparent birefringent polycrystalline material.

Abstract The WAIS (West Antarctic Ice Sheet) Divide deep ice core was recently completed to a total depth of 3405 m, ending 50 m above the bed. Investigation of the visual stratigraphy and grain characteristics indicates that the ice column at the drilling location is undisturbed by any large-scale overturning or discontinuity. The climate record developed from this core is therefore likely to be continuous and robust. Measured grain-growth rates, recrystallization characteristics, and grain-size response at climate transitions fit within current understanding. Significant impurity control on grain size is indicated from correlation analysis between impurity loading and grain size. Bubble-number densities and bubble sizes and shapes are presented through the full extent of the bubbly ice. Where bubble elongation is observed, the direction of elongation is preferentially parallel to the trace of the basal (0001) plane. Preferred crystallographic orientation of grains is present in the shallowest samples measured, and increases with depth, progressing to a vertical-girdle pattern that tightens to a vertical single-maximum fabric. This single-maximum fabric switches into multiple maxima as the grain size increases rapidly in the deepest, warmest ice. A strong dependence of the fabric on the impurity-mediated grain size is apparent in the deepest samples.

AbstractOn 1 December 2011 the West Antarctic Ice Sheet (WAIS) Divide ice-core project reached its final depth of 3405 m. The WAIS Divide ice core is not only the longest US ice core to date, but is also the highest-quality deep ice core, including ice from the brittle ice zone, that the US has ever recovered. The methods used at WAIS Divide to handle and log the drilled ice, the procedures used to safely retrograde the ice back to the US National Ice Core Laboratory (NICL) and the methods used to process and sample the ice at the NICL are described and discussed.

Abstract Active seismic imaging of glaciers and ice sheets is important for constraining inputs to climate models, such as englacial ice fabric and the nature of the basal interface. However, acquiring high-quality seismic data is time-consuming and resource-intensive. Using traditional single-element geophones requires ideal weather conditions (e.g. light winds) and excellent source coupling. In addition, deploying and retrieving these geophones is slow and cumbersome. We have developed a four-element ‘georod’ that enhances signal levels by 20–30dB in a variety of conditions, including blowing snow and poorly coupled source detonations. The long, slender design of these georods makes them easy to deploy and retrieve, allowing researchers to acquire greater line-kilometers of seismic data during field campaigns that are commonly time-constrained.

Seismic attenuation α, or internal friction Q−1, in glacial ice is highly sensitive to temperature, particularly near the melting point. Here we detail a technique to estimate Q and apply it to active source seismic data from Jakobshavn Isbrae, Greenland. We compare our results to measured and modeled temperature profiles of the ice in the region. We find an excellent match, with differences between seismically estimated and modeled temperatures of less than 2°C. Mapping variations in seismic Q through glacial ice thus is shown to allow detailed estimation of englacial temperature profiles, which may be of special value in regions where in situ measurements are logistically difficult.

The thickness of the mantle transition zone beneath Ross Island, and parts of the Transantarctic Mountains and East Antarctic Craton has been mapped using data from the 2000–2003 Transantarctic Mountain Seismic Experiment to determine if, as indicated by some tomographic images, an upper mantle thermal anomaly centered beneath Ross Island is a deep‐seated feature extending into the mantle transition zone. Some 2700 receiver functions have been stacked using a 3D velocity model, revealing Ps conversions from the mantle transition zone discontinuities at depths of 410 and 660 km. Results yield an average nearly uniform transition zone thickness (266 ± 10 km) that is slightly larger than the global average, implying that the upper mantle thermal anomaly does not likely extend into the transition zone. This finding favors explanations for the upper mantle thermal anomaly invoking a plume head or small‐scale convection.

High‐resolution reflection seismic data from Jakobshavn Isbræ, Greenland, reveal complex fabric development. Abundant englacial reflectivity occurs for approximately half the thickness of the ice (the lower half), and disruption of the englacial reflectors occurs in the lower 10–15% of the ice‐thickness. These depths correspond to the higher impurity‐content, and more easily deformed, ice from the Younger Dryas and Last Glacial Maximum to Stage‐3. We conclude that the reflectivity results from contrasting seismic velocities due to changes in the crystal orientation fabric of the ice, and suggest that these fabric changes are caused by variations in impurity loading and subsequent deformation history. These findings emphasize the difference between ice‐divide and ice‐stream crystal orientation fabrics and have implications for predictive ice sheet modeling.