Ever since the first deep ice cores were drilled, it has been a challenge to determine their original, in-situ orientation. In general, the orientation of an ice core is lost as the drill is free to rotate during transport to the surface. For shallow ice cores, it is usually possible to match the adjacent core breaks, which preserves the orientation of the ice column. However, this method fails for deep ice cores, such as the EastGRIP ice core in Northeast Greenland. We provide a method to reconstruct ice core orientation using visual stratigraphy and borehole geometry. As the EastGRIP ice core is drilled through the Northeast Greenland Ice Stream, we use information about the directional structures to perform a full geographical re-orientation. We compared the core orientation with logging data from core break matching and the pattern of the stereographic projections of the crystals’ c-axis orientations. Both comparisons agree very well with the proposed orientation method. The method works well for 441 out of 451 samples from a depth of 1375–2120 m in the EastGRIP ice core. It can also be applied to other ice cores, providing a better foundation for interpreting physical properties and understanding the flow of ice.

Cryo-Raman tomography allows us, for the first time, to determine accurate morphologies and volumes of natural air hydrates in Antarctic ice cores. The measurements show complex growth structures that are not accounted for by the available models of hydrate formation.

A new densification model, which simulates the effect of impurities on the densification of polar firn, is presented. The classical densification models of Herron and Langway (1980) and Pimienta and Barnola (Barnola and others, 1991) are modified by assuming that the activation energy for deformation is reduced by the impurities. Motivated by recent observations, the impurity effect is formulated on an empirical basis using the seasonally varying Ca2+ ion concentration. Excellent agreement between simulated and measured high-resolution density profiles confirms the new approach. The same parameterization applies for Greenland and Antarctica despite the one order of magnitude difference in impurity concentration. The new models allow us, for the first time, to simulate the density layering in firn down to the firn–ice transition. Our results emphasize the importance of impurities and density layering for the air entrapment and for dating gas records of deep ice cores, in particular for glacial climate conditions where the impurity concentrations are 10–100-fold higher than in modern firn.

This study aims to demonstrate that deep ice cores can be synchronized using internal horizons in the ice between the drill sites revealed by airborne radio-echo sounding (RES) over a distance of >1000km, despite significant variations in glaciological parameters, such as accumulation rate between the sites. In 2002/03 a profile between the Kohnen station and Dome Fuji deep ice-core drill sites, Antarctica, was completed using airborne RES. The survey reveals several continuous internal horizons in the RES section over a length of 1217 km. The layers allow direct comparison of the deep ice cores drilled at the two stations. In particular, the counterpart of a visible layer observed in the Kohnen station (EDML) ice core at 1054 m depth has been identified in the Dome Fuji ice core at 575 m depth using internal RES horizons. Thus the two ice cores can be synchronized, i.e. the ice at 1560 m depth (at the bottom of the 2003 EDML drilling) is ∼49ka old according to the Dome Fuji age/depth scale, using the traced internal layers presented in this study.

We investigated the large-scale (10–1000 m) and small-scale (mm–cm) variations in size, number and arrangement of air bubbles in the EPICA Dronning Maud Land (EDML) (Antarctica) ice core, down to the end of the bubble/hydrate transition (BHT) zone. On the large scale, the bubble number density shows a general correlation with the palaeo-temperature proxy, δ18O, and the dust concentration, which means that in Holocene ice there are fewer bubbles than in glacial ice. Small-scale variations in bubble number and size were identified and compared. Above the BHT zone there exists a strong anticorrelation between bubble number density and mean bubble size. In glacial ice, layers of high number density and small bubble size are linked with layers with high impurity content, identified as cloudy bands. Therefore, we regard impurities as a controlling factor for the formation and distribution of bubbles in glacial ice. The anticorrelation inverts in the middle of the BHT zone. In the lower part of the BHT zone, bubble-free layers exist that are also associated with cloudy bands. The high contrast in bubble number density in glacial ice, induced by the impurities, indicates a much more pronounced layering in glacial firn than in modern firn.

A new radioscopic imaging technique has been developed to measure firn density in unprecedented resolution and accuracy even when the porosity is low or the geometry of a core or piece of core is not perfect. The technique is based on an X-ray microfocus computer tomograph (ICE-CT) designed especially for ice-core applications. Applied on an archive piece of the Antarctic firn core B32 drilled in Dronning Maud Land in 1998, the obtained density profile shows a strong correlation with the calcium ion concentration as found previously in Greenland. Given the impurity–density relationship found previously in Greenland, our result suggests both improved accuracy of the new density measurements and an impurity–density relationship with a similar magnitude in Greenland to that on the Antarctic plateau. Our measurements provide first evidence that the impurity–density relationship is a universal feature of polar firn and that the calcium ion concentration can serve as a proxy to describe quantitatively the effect of the impurities on densification.

Measurements of N2/O2 ratios inside individual air bubbles at various depths in the EDML (Antarctic) ice core are presented here. The small bubbles (diameter less than ~200 µm) in deeper ice are significantly enriched in O2 compared to the larger bubbles. The N2/O2 ratios show a systematic dependence on bubble size which is not the case for bubbles in shallower ice. This is interpreted as an effect of pressure relaxation during storage of the cores.

Ice in polar ice sheets undergoes deformation during its flow towards the coast. Deformation and recrystallization microstructures such as subgrain boundaries can be observed and recorded using high-resolution light microscopy of sublimation-edged sample surfaces (microstructure mapping). Subgrain boundaries observed by microstructure mapping reveal characteristic shapes and arrangements. As these arrangements are related to the basal plane orientation, full crystallographic orientation measurements are needed for further characterization and interpretation of the subgrain boundary types. X-ray Laue diffraction measurements validate the sensitivity of different boundary types with sublimation used by microstructure mapping for the classification. X-ray Laue diffraction provides misorientation values of all four crystal axes. Line scans across a subgrain boundary pre-located by microstructure mapping can determine the rotation axis and angle. Together with the orientation of the subgrain boundary this yields information on the dislocation types. Tilt and twist boundaries composed of dislocations lying in the basal plane, and tilt boundaries composed of nonbasal dislocations were found. A statistical analysis shows that nonbasal dislocations play a significant role in the formation of all subgrain boundaries.

Static (or ‘normal’) grain growth, i.e. grain boundary migration driven solely by grain boundary energy, is considered to be an important process in polar ice. Many ice-core studies report a continual increase in average grain size with depth in the upper hundreds of metres of ice sheets, while at deeper levels grain size appears to reach a steady state as a consequence of a balance between grain growth and grain-size reduction by dynamic recrystallization. The growth factor k in the normal grain growth law is important for any process where grain growth plays a role, and it is normally assumed to be a temperature-dependent material property. Here we show, using numerical simulations with the program Elle, that the factor k also incorporates the effect of the microstructure on grain growth. For example, a change in grain-size distribution from normal to log-normal in a thin section is found to correspond to an increase in k by a factor of 3.5.

Air bubbles in ice cores play an essential role in climate research, not only because they contain samples of the palaeoatmosphere, but also because their shape, size and distribution provide information about the past firn structure and the embedding of climate records into deep ice cores. In this context, we present profiles of average bubble size and bubble number for the entire EDML (Antarctica) core and the top 600 m of the EDC (Antarctica) core, and distributions of bubble sizes from selected depths. The data are generated with an image-processing framework which automatically extracts position, orientation, size and shape of an elliptical approximation of each bubble from thick-section micrographs, without user interaction. The presented software framework allows for registration of overlapping photomicrographs to yield accurate locations of bubble-like features. A comparison is made between the bubble parameterizations in the EDML and EDC cores and data published on the Vostok (Antarctica) ice core. The porosity at the firn/ice transition is inferred to lie between 8.62% and 10.48% for the EDC core and between 10.56% and 12.61 % for the EDML core.

Analyses of shallow cores obtained at the European Project for Ice Coring in Antarctica (EPICA) drilling site Kohnen station (75°00′ S, 00°04′ E; 2892 m a.s.l.) on the plateau of Dronning Maud Land reveal the presence of conserved snow dunes in the firn. In situ observations during three dune formation events in the 2005/06 austral summer at Kohnen station show that these periods were characterized by a phase of 2 or 3 days with snowdrift prior to dune formation which only occurred during high wind speeds of >10 m s-1 at 2 m height caused by the influence of a low-pressure system. The dune surface coverage after a formation event varied between 5% and 15%, with a typical dune size of (4 ± 2) m × (8 ± 3) m, a maximum height of 0.2 ± 0.1 m and a periodicity length of about 30 m. The mean density within a snow dune varied between 380 and 500 kg m-3, whereas the mean density at the surrounding surface was 330 ± 5 kgm-3. The firn cores covering a time-span of 22 ± 2 years reveal that approximately three to eight events per year occurred, during which snow dunes had been formed and were preserved in the firn.

Subgrain boundaries revealed as shallow sublimation grooves on ice sample surfaces are a direct and easily observable feature of intracrystalline deformation and recrystallization. Statistical data obtained from the EPICA Dronning Maud Land (EDML) deep ice core drilled in East Antarctica cannot detect a depth region of increased subgrain-boundary formation. Grain-boundary morphologies show a strong influence of internal strain energy on the microstructure at all depths. The data do not support the classical view of a change of dominating recrystallization regimes with depth. Three major types of subgrain boundaries, reflecting high mechanical anisotropy, are specified in combination with crystal-orientation analysis.

Clear evidence for the formation of mixed clathrate hydrates of air and hydrochlorofluorocarbon densifier (known as HCFC-141b, sometimes also called R-141b) is found by means of synchrotron X-ray diffraction and Raman spectroscopy on a sample recovered from the bottom of the EPICA Dronning Maud Land deep borehole in Antarctica. Subglacial water (SGW) appears to have reacted with the drilling liquid to build a large lump of clathrate hydrate. The hydrate growth may well have been accelerated by the stirring of the SGW–densifier mixture during drilling. Moreover, dissolved air in the SGW appears to have participated in the formation of mixed hydrates of air and HCFC-141b as evidenced by the concomitant appearance of Raman signals from both constituents. Our findings elucidate to some extent the meaning of earlier accounts of the formation of ‘heavy chips’ that may sink to the bottom of the borehole, possibly affecting or even impeding the drilling advance. These observations raise concerns with respect to the use of HCFC-141b densifiers in ice-core drilling liquids under warm ice conditions.

Results of laboratory uniaxial compression tests over the stress range 0.18–0.52 MPa and the strain range 0.5–8.6% at approximately –5 and –20°C are presented. Grain-size analysis and comparisons with annealing tests confirm that grain-growth reducing processes are active during deformation. Microstructural observations reveal that subgrain-rotation recrystallization and grain-shape changes due to strain-induced grain-boundary migration are the causes of the grain-growth deceleration. Further results from microstructural observations show that obstacle formation by dislocation walls and subgrain boundaries is the reason for isotropic hardening during creep. Subgrainboundary types that are likely to be relevant for studies on the activity of different dislocation types are described.

Automatic c-axes analyzers have been developed over the past few years, leading to a large improvement in the data available for analysis of ice crystal texture. Such an increase in the quality and quantity of data allows for stricter statistical estimates. The current textural parameters, i.e. fabric (crystallographic orientations) and microstructure (grain-boundary networks), are presented. These parameters define the state of the polycrystal and give information about the deformation undergone by the ice. To reflect the findings from automatic measurements, some parameter definitions are updated and new parameters are proposed. Moreover, a MATLAB® toolbox has been developed to extract all the textural parameters. This toolbox, which can be downloaded online, is briefly described.

This work presents a method of mapping deformation-related sublimation patterns, formed on the surface of ice specimens, at microscopic resolution (3–4 μm pixel−1). The method is based on the systematic sublimation of a microtomed piece of ice, prepared either as a thick or a thin section. The mapping system consists of an optical microscope, a CCD video camera and a computer-controlled xy-stage. About 1500 images are needed to build a high-resolution mosaic map of a 4.5 × 9 cm section. Mosaics and single images are used to derive a variety of statistical data about air inclusions (air bubbles and air clathrate hydrates), texture (grain size, shape and orientation) and deformation-related features (subgrain boundaries, slip bands, subgrain islands and loops, pinned and bulged grain boundaries). The most common sublimation patterns are described, and their relevance for the deformation of polar ice is briefly discussed.

The densification of dry polar snow and firn results in a continuous increase of density with depth accompanied by significant density fluctuations within seasonal layers. Density measurements of high spatial resolution reveal a persistent minimum of density fluctuations in the vicinity of the snow–firn transition (0.55–0.65 g cm-3) in firn-core records. In this study we give an explanation for the fluctuation minimum by applying a new method of X-ray microtomography to obtain three-dimensional (3-D) structural data of a Greenland firn core. At 13 different depths between 10 and 78 m a set of 16 samples of 40 cm total length for each depth interval was measured. A reconstructed firn segment of 40 cm covers 1–2 years of snow accumulation. Using digital image analysis techniques, different structural parameters are estimated including 3-D pore and particle sizes and specific surface areas. It is shown that the densification rates of snow and firn layers consisting of coarse particles are much higher than those consisting of fine particles within the same depth interval. This causes a density crossing of fine- and coarse-grained layers with a minimum of density variations at the crossover point. This crossing-over implies that formerly dense layers in the seasonal density signal are not of the same origin as dense layers in the deeper part of the firn column and that the seasonal density signal will totally change shape with depth. It is speculated that in coarse- and fine-grained firn the dominant mechanism of densification acts over different regimes of density.

Fabric analysis of the upper 1300 m of the Dome C (East Antarctica) ice core reveals a slight clustering tendency of c axes towards the vertical, which gradually enhances with depth from an initially isotropic orientational distribution of c axes at the free surface. Such a strain-induced anisotropy is compatible with the expected macroscale stress state in a dome, i.e. dominated by vertical compression. Yet, when one analyzes the orientational distribution of the visible gliding layers of individual crystallites (slip bands), the evidence is quite contrasting. Direct observation of slip bands in samples from the Dome C ice core taken from different depths (204–1291 m) indicates a higher slip activity in nearly horizontal planes, in such a manner that >60% of the detected slip bands have an inclination of <30˚ with respect to the horizontal. Furthermore, the observed slip activity is not symmetric, i.e. the number of slip bands discerned at 20˚, say, is usually not comparable with the number found at 160˚. Such features are not consistent with the predicted slip activity induced by compression and/or extension. In this work, we present evidence for this unexpected orientational distribution of slip bands and discuss some of the possible causes. Natural and artificial agents are investigated, together with their respective consequences for ice-sheet modeling and ice-core processing. Additionally, we show the occurrence of bent slip bands in certain crystallites. Such a bending represents an early stage of polygonization, and highlights the strong inhomogeneity of deformation at the crystal level. Moreover, it indicates that polygonization might be mathematically interpreted as a continuous process of rotation, characterized by the divergence of c axes from a common direction.

Detailed measurements of crystal outlines and fabrics have been performed on 35 000 crystals in fifteen 10 × 20 cm2 vertical thin sections from the North Greenland Icecore Project (NorthGRIP) ice core, evenly distributed in the depth interval 115–880m. The crystals exhibit important changes over this period. As the ice gets older the mean crystal area increases towards a constant value, the shape of the crystals becomes increasingly irregular, and the area distribution of crystals develops from a single log-normal distribution into a bimodal lognormal distribution. The c-axis fabric of the ice shows a smooth development of an increasingly stronger vertical fabric with depth, and the formation of a weak vertical girdle. Already in the younger samples the fabric is rather strongly oriented towards vertical. The fabric and the area of individual crystals are found not to correlate. A simple model, which takes into account the vertical strain of the ice, is applied in an attempt to determine the crystal growth rate at NorthGRIP.

A study of c-axis orientations in the upper 1500m of the Dome C (East Antarctica) deep ice core has been carried out using an automatic ice-fabric analyzer (AIFA). Twenty-nine vertical and a few horizontal thin sections from different depths in the core have been analyzed. Several statistical parameters describing fabric strength and fabric shapes have been calculated from the c-axis orientation data. The fabric diagrams display a near-random c-axis orientation distribution in the uppermost parts of the ice sheet. A tendency of c-axis rotation towards a broad single-maximum fabric is observed in the lowest part of the studied interval. The fabric development at Dome C thus appears typical for ice-sheet summit and dome sites. The fabric development at Dome C is compared with the fabric evolution in the Dome F and GRIP ice cores, and data on crystal size obtained with image-analysis techniques are presented. Studies of misorientation angles between the c axes of neighbouring crystals reveal little evidence for polygonization, but microscopic observations show that sub-grain boundaries are present in half of the grains at any depth.