Hostname: page-component-76fb5796d-zzh7m Total loading time: 0 Render date: 2024-04-25T19:41:16.480Z Has data issue: false hasContentIssue false

Raman tomography of natural air hydrates

Published online by Cambridge University Press:  10 July 2017

Christian Weikusat*
Affiliation:
Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany
Sepp Kipfstuhl
Affiliation:
Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany
Ilka Weikusat
Affiliation:
Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany
*
Correspondence: Christian Weikusat <christian.weikusat@awi.de>
Rights & Permissions [Opens in a new window]

Abstract

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.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2015

Introduction

Polar ice is a composite of ice, air and various impurities of atmospheric origin. This makes ice cores unique climate archives, because the air inclusions in the polar ice sheets are the only direct source for measuring and reconstructing the composition of the atmosphere of the past (e.g. Reference Raynaud, Barnola, Chappellaz, Blunier, Indermühle and StaufferRaynaud and others, 2000). To correctly interpret the findings, a firm understanding of the processes that may be able to change the gas content of the air inclusions is mandatory.

The precursors of polar ice, snow and firn, form a porous permeable material up to 120 m thick on top of the polar ice sheets. Atmospheric gases circulate through the connected pores and are entrapped at the transition from firn to ice, where the network of open pores closes off to form individual gas inclusions (air bubbles (ABs)) (Reference Schwander and StaufferSchwander and Stauffer, 1984). With time, the ABs move deeper into the ice sheet and start shrinking as the hydrostatic pressure of the overlying ice increases. When the gas pressure inside the AB exceeds the dissociation pressure of hydrate formation the ABs begin to convert to air hydrates (AHs). AHs are crystalline structures composed of water molecules and containing cavities suitable for housing one or two gas molecules at a time. In contrast to the known phase diagrams (Reference Miller, Whalley, Jones and GoldMiller, 1973), AB–AH conversion in polar ice takes place over a large depth range (transition zone) of several hundred meters (Reference Uchida, Yasuda, Oto, Shen and OhmuraUchida and others, 2014). Microstructure mapping images of the EPICA Dronning Maud Land (EDML) (Antarctica) ice core show first AHs at ∼720 m depth; ABs were observed to coexist with AHs down to ∼1225 m depth (Reference BendelBendel and others, 2013). The widely accepted reason for the thickness of the transition zone is the slow nucleation of hydrates, after which the remaining transition proceeds comparatively quickly (e.g. Reference Shoji and LangwayShoji and Langway, 1982). This is supported by the fact that no intermediate stages of the AB–AH transition have been found in ice cores (Reference NaritaNarita and others, 1999; Reference Lipenkov and HondohLipenkov, 2000; Reference Kipfstuhl, Pauer, Kuhs and ShojiKipfstuhl and others, 2001).

The influence of the AB–AH conversion processes on the measured gas contents has been indicated previously: Reference IkedaIkeda and others (1999) and Reference Ikeda-Fukazawa, Hondoh, Fukumura, Fukazawa and MaeIkeda-Fukazawa and others (2001) presented results of Raman measurements on individual air inclusions that clearly show gas fractionation in the transition zone, where initial AHs have lower N2/O2 ratios (i.e. they are oxygen-enriched) and the late ABs much higher N2/O2 ratios (i.e. nitrogen-enriched) than anticipated. Furthermore, anomalies in mixing ratios and isotopic composition of O2, N2 and argon were reported by Reference Bender, Sowers and LipenkovBender and others (1995). Reference LüthiLüthi and others (2010) showed significant deviations of CO2 and O2/N2 ratios in and just below the transition zone, that are caused by diffusion and fractionation processes during the AB–AH conversion. Little is known about these processes, and while a reduction of resolution in the transition zone helps with measuring CO2, the processes involved have to be understood to ascertain the correct storage of CO2 in the climate archive.

The AHs produced by a direct conversion from ABs are called primary air hydrates (pAHs) (Reference Kipfstuhl, Pauer, Kuhs and ShojiKipfstuhl and others, 2001). Two types of this transformation have been proposed by Reference Salamatin, Lipenkov, Hondoh and IkedaSalamatin and others (1999):

  • Type A: The hydrate nucleates and rapidly forms a closed hydrate shell, covering the entire bubble wall. In a second, much slower, step the remaining bubble volume is spherically filled with hydrate, where water molecules diffusing through the closed shell towards the growing AH are the rate-limiting factor.

  • Type B: The AH does not cover the bubble perimeter but proceeds as a flat crystallization front through the AB. Water diffusion through AH is not needed in this case, which indicates this process is much faster.

Additionally, Reference Kuhs, Klapproth, Chazallon and HondohKuhs and others (2000) discussed the production of monocrystalline or polycrystalline AHs through a single nucleus or polynucleation respectively. At the beginning of the transition zone, the packing density inside an AH is much larger than inside a bubble. Therefore, if the AH growth uses the gas content of the transforming bubble, the resulting AH will have a significantly smaller volume than the originating bubble. Differences in bubble pressure and diffusion of gases from neighboring bubbles or the ice matrix to the newly forming AH can influence the actual volume mismatch between AH and the original bubble (Reference Salamatin, Lipenkov, Hondoh and IkedaSalamatin and others, 1999), complicating the deduction of bubble radii from AH sizes.

While most AHs appear smooth-surfaced and very clear (smooth air hydrates, sAHs), AHs with visible structure have been found in all ice cores and have been described as ‘cloudy’, ‘irregular’ or ‘graupel-like’ (Reference Shoji and LangwayShoji and Langway, 1982; Reference Nakahara, Shigesato, Higashi, Hondoh and LangwayNakahara and others, 1988; Reference Uchida, Hondoh, Mae, Lipenkov and DuvalUchida and others, 1994a; Reference NaritaNarita and others, 1999; Reference Pauer, Kipfstuhl, Kuhs and ShojiPauer and others, 1999, Reference Pauer, Kipfstuhl, Kuhs and Shoji2000; Reference Lipenkov and HondohLipenkov, 2000; Reference Kipfstuhl, Pauer, Kuhs and ShojiKipfstuhl and others, 2001). An example from the NorthGRIP ice core is shown in Figure 1a. Carefully chosen illumination conditions reveal three-dimensional (3-D) structures inside the pAH (Fig. 1b). The genesis of these structures, which have up to now only been accessible to visual inspection, is still unclear. Reference Salamatin, Lipenkov, Hondoh and IkedaSalamatin and others (1999) suggested that clusters of structured AHs could represent an AH shell fragmented by hydrostatic pressure during the type A AH growth. The amount of structural data retrievable by microscopic methods is very limited, as no information about the actual 3-D structure can be obtained. While timescales of several thousand years are not adequately reproducible in the laboratory, numerous researchers have performed laboratory experiments on bubble–hydrate conversions (e.g. Reference Uchida, Hondoh, Mae, Duval, Lipenkov, Maeno and HondohUchida and others, 1992, Reference Uchida, Hondoh, Mae, Lipenkov and Duval1994a; Reference Ikeda, Uchida and MaeIkeda and others, 1993; Reference Shimada and HondohShimada and Hondoh, 2004) and also Raman spectroscopy on synthetic air hydrates (Reference Chazallon, Champagnon, Panczer, Pauer, Klapproth and KuhsChazallon and others, 1998). We present for the first time measured data of the intricate 3-D structure of natural air hydrates.

Fig. 1. AH from the NorthGRIP core: (a) transmitted light; (b) asymmetric illumination reveals 3-D structure.

Methods

Samples

We used samples from the EDML ice core from 1005 m, 1048 m, 1083 m and 1084 m (all depths are inside the transition zone). The EDML core was drilled between 2001 and 2006 at Kohnen Station, Dronning Maud Land, East Antarctica (75°00′ S, 0°04′ E; 2882 m a.s.l.) within the framework of the European Project for Ice Coring in Antarctica (Reference Oerter, Graf, Meyer and WilhelmsOerter and others, 2004). The samples studied were drilled in the 2002/03 season and have been stored at −30°C ever since.

The samples were cut, microtomed and sublimated to produce clean surfaces (see Reference Weikusat, Freitag and KipfstuhlWeikusat and others, 2012, for details). Overview images taken with a 10× lens and the automatic image stitching function of the WITec control software were used to locate suitable AHs for measurements.

Raman tomography

The novel Alfred Wegener Institute cryo-Raman system is equipped with fibre optics that allow for a split set-up where the microscope unit is in the cold laboratory and the controller and spectrometer are at room temperature. This set-up eliminates the need for sample transfer between the cold laboratory and a cooling cell mounted to the microscope. As such transfers inevitably expose the samples to large temperature gradients and moisture from the air they have been the most critical point in cryo-analysis. The cryo-Raman set-up allows for sample preparation in the same laboratory where the microscope unit is situated, and transfer to the system is virtually riskless. The set-up also provides much more stable surfaces and reduced sublimation conditions, allowing for long-time measurements of up to several days.

The system consists of a WITec alpha 300 M+, UHTS 300 spectrometer using a grating with 600 grooves mm−1, Peltier-cooled CCD detector, frequency-doubled Nd:YAG laser (532 nm), and Zeiss Epiplan 50× LWD objective (NA 0.35). The laser spot size for this set-up is ∼1μm, which defines the spatial resolution limit of the measurements. The confocal hole is provided by the diameter of the light-transmitting fibre from microscope to spectrometer. We used a 100 μm fibre which gives a good compromise between confocality and signal intensity. The optical fibre transmitting the laser to the microscope is polarization-maintaining, while the fibre to the spectrometer depolarizes the light for a better detection efficiency. The temperature in the cold laboratory, and thus the temperature of the sample, were held at −15°C for all measurements.

Hyperspectral large-area stack scans map a rectangular area of the sample at a given step size, where a full spectrum is acquired for each pixel. The microscope turret is then moved a predefined step downwards and another scan with the same parameters is acquired for the new depth. One image is generated from each scan in the stack. These images can then be used to create a 3-D reconstruction of the measured AH.

Data analysis

The Raman spectra of ice and air hydrate are very similar, both being crystalline structures made from H2O molecules (Reference Schicks, Erzinger and ZiemannSchicks and others, 2005). Additionally, due to the choice of the 100 μm optical fibre and the implied reduced confocality, there is always a small contribution from the ice matrix in every recorded spectrum. The differences used here to distinguish between ice- and air-hydrate phase are shown in Figure 2. The low-energy mode of ice at ∼216 cm−1 is shifted ∼3 cm−1 to lower values (∼213 cm−1) in the hydrate while the band broadens considerably (FWHM (full width at half-maximum) in ice: ∼27 cm−1; FWHM in AH: ∼37 cm−1). The shoulder at ∼290 cm−1 in the ice spectrum is absent from the AH spectrum. A broadband weighted center-of-mass filter from 180 to 340 cm−1 was used to extract the phase boundaries between ice and air hydrate (WITec Project software). To obtain binary images (Fig. 3), pixels with a value >250 cm−1 were defined as hydrate, and pixels with lower values as ice. It is important to note that other parameters (e.g. the height of the O2 and N2 bands (∼1552 cm−1 and ∼2326 cm−1) or shifts in the OH region (∼3000–3600 cm−1)) yield the same results qualitatively, albeit with worse signal/noise ratio. The resulting maps were exported as binary images, with black representing AH and white, ice. The obtained slices were combined to produce a 3-D reconstruction of the measured AH. Automatic segmentation of the slices and initial mesh generation were performed with ImageJ (Reference Abràmoff, Magalhães and RamAbràmoff and others, 2004). Smoothing of the initial mesh and calculation of geometric parameters were done with Meshlab (Visual Computing Lab ISTI – CNR, 2013).

Fig. 2. Differences in the Raman spectra of ice and air hydrate that were used for imaging (averages of 100 spectra with accumulation time of 0.3 s per spectrum). a.u.: arbitrary units.

Fig. 3. Binary slices generated from the Raman stack scan of pAH-1 (black: air hydrate; white: ice). The relative depth is given for each image; the scale is given on the photograph (lower right). See Figure 4 for the corresponding 3-D reconstruction.

Peak fitting and calibration

Background subtraction and fitting of Voigt functions to the Raman bands were done with the open-source software package FITYK 1.2.1 (Reference WojdyrWojdyr, 2010). The intensities of the main vibrational bands of N2 (∼2326 cm−1) and O2 (∼1552 cm−1) were used to calculate N2/O2 ratios. To calibrate the instrument response, several measurements of room air were averaged, yielding a raw N2/O2 ratio of 2.97 ± 0.04 (mean value and standard deviation of all air measurements). Assuming the laboratory air composition to be N2/O2 = 3.73, every calculated ratio was multiplied by a correction factor of 1.26.

N2/O2 ratios and crystal orientation

It has been shown that the intensity of the Raman bands is a function of the orientation of the AH crystal to the polarization direction of the incident laser (Reference Ikeda, Fukazawa, Mae, Hondoh and LangwayIkeda and others, 1997). Ikeda and co-workers found values between 8.06 and 10.33 for the N2/O2 ratios measured at different laser polarization angles for measurements on a well-defined crystal plane of a monocrystalline AH from the Vostok (East Antarctica) ice core. However, the peaks for N2 and O2 scale homogeneously for most polarization angles (which is reasonable, as they inhabit the same cavities). Additionally, the laser used for this study was only ∼90% polarized and we measured on an arbitrary plane of hydrates, so we assume that the effect is much less pronounced. To be on the safe side we should assume an additional margin of uncertainty of ∼±1.1.

Results

Structure and volumes

Figure 4a presents a microscopic image of a structured pAH (pAH-1) from 1005 m depth. It is almost perfectly round, with a diameter of ∼130 μm. The typical ‘cloudy’ appearance and rough surface are directly visible. With a step size of 3 μm in the x-y plane and a z-step size of 10 μm, the whole stack scan is comprised of 48 355 full spectra. An acquisition time of 0.2 s per spectrum was used, amounting to ∼3 hours total measurement time. The 3-D reconstruction is shown in Figure 4bd. pAH-1 is not solid but incorporates a large ice intrusion that is connected to the surrounding matrix ice by two holes in the AH shell (Fig. 4c). The ice volume is asymmetrically distributed inside the AH (Fig. 4d). The volume of pAH-1 is 9.89 × 105 μm3, and the volume of the ice intrusions is 2.12 × 105 μm3 (Fig. 4d), giving an ice: AH ratio of 0.2. The combined area of the holes in the AH shell is ∼1687 μm2.

Fig. 4. pAH-1 (depth 1005 m). (a) Microscopic image (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation (steps are an artifact of large z-step size). (c) Side view; large holes in the shell are evident. (d) Semi-transparent side view (the holes are to the left) with visible ice intrusion.

Figure 5a shows a microscopic image of a structured pAH (pAH-2) from 1048 m depth, exhibiting similar features to pAH-1 (Fig. 4a). The mapped section of the AH is indicated. The diameter of pAH-2 is ∼134 μm. A volume of 55 μm 165 μm × 150 μm was scanned with an x-y step size of 2 μm and a z-step size of 10 μm, totaling 34 031 full spectra. With an acquisition time of 0.2 s per spectrum the whole stack took ∼2 hours to complete. Compared to pAH-1, the ice intrusion inside pAH-2 has a more complex structure. Holes in the AH shell that connect the ice intrusion to the matrix can be seen in Figure 5b and c. The AH volume in this section is 3.45 × 105 μm3, and the ice intrusion volume is 1.07 × 105 μm3, giving an ice:AH ratio of 0.31.

Fig. 5. pAH-2 (depth 1048 m). (a) Microscopic image with highlighted scan area (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation (steps are an artifact of large z-step size). (c) Side view of the section with evident ice inclusion. (d) Side view from the opposite side.

Figure 6a shows two connected AHs from 1084 m depth: a very clear oval sAH with smooth boundaries (sAH-1) above an irregular and structured pAH (pAH-3). Several solid inclusions (black dots) can be seen inside sAH-1. The longest diameter of sAH-1 is ∼95 μm, and the diameter of pAH-3 is ∼120 μm. The scan was done with a step size of 3 μm for all three directions, amounting to 212 333 full spectra. With an acquisition time of 0.3 s per spectrum, the total measurement time for the stack was ∼18 hours. The 3-D reconstruction in Figure 6b exhibits a continouus ice intrusion through pAH-3 (the top part of sAH-1 was outside the scan range). Figure 6c and d clearly show that both AHs are connected. Figure 6d visualizes the dimensions of the ice intrusion through pAH-3. The volume of the ice intrusion is 9.17 × 104 μm3. The volume of pAH-3 could not be calculated reliably due to the connection of both AHs. The volume of the measured connected structure is 1.59 × 106 μm3. As the top part of sAH-1 was outside the scan range, this calculation only yields a lower limit for the whole structure.

Fig. 6. Connected AHs pAH-3 and sAH-1 (depth 1084 m). (a) Microscopic images with focus on sAH-1 (left) and pAH-3 (right). (b) 3-D reconstruction viewed in the same orientation (top part of sAH-1 out of scan range; scale bar for (b–d)). Note the hole through pAH-3. (c) Side view reveals a connection of both AHs. (d) Semi-transparent side view with visible ice intrusion.

Figure 7a shows a smooth AH from 1084 m depth (sAH-2). The scan has a step width of 4 μm for all directions and consists of 24 000 full spectra. With an acquisition time of 0.3 s per spectrum the measurement time was ∼2 hours. The 3-D reconstructions (Fig. 7bd) show a slightly flattened AH with smooth outlines and no internal structuration. The volume of sAH-2 is 4.11 × 106 μm3.

Fig. 7. sAH-2 (depth 1084 m). (a) Microscopic image. (b) 3-D reconstruction viewed in the same orientation (scale bar for (b–d)). (c) View perpendicular to (b). (d) View perpendicular to (c).

Figure 8a shows a faceted sAH from 1084 m depth (sAH-3). The scan has a step size of 3 μm for all directions and the stack consists of 23 466 full spectra. With an acquisition time of 0.3 s per spectrum the measurement time was ∼2 hours. Multiple twinned octahedra are visible in the 3-D reconstruction (Fig. 8bd). The volume of sAH-3 is 1.12 × 105 μm3.

Fig. 8. Faceted AH (sAH-3, depth 1084 m). (a) Microscopic image (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation. (c) View perpendicular to (b). Several twinned octahedra are discernible. (d) View perpendicular to (c).

Figure 9a shows a plate-like sAH from 1083 m depth (sAH-4). The scan step size was 3 μm for all directions. With an acqusition time of 0.3 s per spectrum and a total of 341 185 spectra the measurement time was ∼28 hours. The dimensions of sAH-4 are 245 μm × 125 μm × 20 μm; the volume is 6.5 × 105 μm3. The lower left part of sAH-4 that appears to be slightly tilted in Figure 9a is only partially connected to the rest of the AH (Fig. 9c).

Fig. 9. Plate-like sAH-4 (depth 1083 m). (a) Microscopic image (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation. (c) View perpendicular to (b).

N2/O2 ratios

To determine N2/O2 ratios, several spectra were acquired on different points for each AH and the resulting ratios were averaged. Table 1 shows the results for all presented AHs. The ratio for pAH-1 is 3.60 ± 0.07, which is close to the atmospheric value of 3.7. pAH-2, pAH-3 and sAH-1 are distinctly enriched in nitrogen with ratios of 5.51 ± 0.15, 4.84 ± 0.10 and 4.91 ± 0.15, respectively. The connected AHs, pAH-3 and sAH-1, exhibit very similar ratios. The solid and smooth sAH-2 has a distinctly higher ratio of 8.30 ± 0.20; the plate-like sAH-4 exhibits a similarly high nitrogen content, with an N2/O2 ratio of 8.09 ± 0.20. The highest relative nitrogen content was found in the faceted sAH-3 with an N2/O2 ratio of 17.49 ± 0.70.

Table 1. Averaged N2/O2 ratios and standard deviations

The homogeneity of the mixing ratios inside an AH was checked for pAH-1 by averaging all spectra for each slice and determining the N2/O2 ratios from the averaged spectra. The values in Table 2 are rather homogeneous and show no systematic alterations with relative depth in the AH.

Table 2. N2/O2 ratios from averaged spectra for each slice of pAH-1 (cf. Fig. 3)

Discussion

Structure and formation processes

The studied pAHs with complex internal structures and incorporated ice bodies are very round and seem to have clearly retained the shape of the originating bubble. A destructive collapse of the AH shell during growth is certainly more likely to form either irregular agglomerations of splinter-like pieces (as observed by, e.g., Reference Salamatin, Lipenkov, Hondoh and IkedaSalamatin and others, 1999), or solid, volume-optimized AH bodies (and eventually sAHs), if the process is combined with subsequent recrystallization. We therefore suggest that the intrusion of ice into the former AB volume is a syngenetic and non-destructive process. Reference Salamatin, Hondoh, Uchida and LipenkovSalamatin and others (1998) found that the rate-limiting step of post-nucleation AH conversion is the plastic deformation of the ice-AH volume around the bubble. They also deduced that AHs are harder than the surrounding ice. This was later confirmed by Durham and co-workers in a series of ice-hydrate rheology experiments, where both type I and type II hydrates with various guest molecules proved to be up to 20 times harder than ice under identical strain conditions (Reference Durham, Kirby and SternDurham and others, 2000, Reference Durham, Kirby, Stern and Zhang2003, Reference Durham, Stern, Kirby and Circone2005, Reference Durham, Prieto-Ballesteros, Goldsby and Kargel2010). We propose the following path for the AB–AH transition yielding structured pAHs: The AH nuclei start growing into the bubble, without forming a closed hydrate shell around the whole bubble. With proceeding AH growth the pressure inside the remnant bubble decreases, causing the bubble volume to shrink and the surrounding ice to flow simultaneously into the former bubble volume. In this scenario the diffusion of water molecules is not inhibited by a closed hydrate shell and the whole growth process can happen very fast, similar to the B-type formation described by Reference Salamatin, Lipenkov, Hondoh and IkedaSalamatin and others (1999), yet the B-type transformation maintains a flattened interface between growing AH and gas phase throughout the process and should therefore result in solid AH bodies without internal structures. As AH growth can only continue as long as the bubble pressure is above the dissociation pressure of the AH, the viscosity of the surrounding ice is an important factor for the growth velocity. (If we assume a spherical geometry for the ice intrusion in pAH-1, we can gain a rough estimate for the timescales needed to fill this volume with ice: as a first approximation, we set the driving stress to the hydrostatic pressure, P = 9 MPa. The strain rate for single crystal ice can then be extracted from figure 2 of Reference Duval, Ashby and AndermanDuval and others (1983). As the orientation of the crystallite is unknown, we chose a medium rate of . To obtain the velocity, we integrated the strain rate, considering only the x direction for further simplification. With x=110 μm, ice intrusion volume V =2.12 × 105 μm3 and opening area of the hydrate shell A = 1687 μm2, a time of t = 3.2 hours can be readily calculated. While this may give a lower boundary for the formation of pAH-1, it is by no means intended to represent a mathematical treatment of the microstructural flow of ice, which is well outside the scope of this publication.)

After all gas has been used up for the AH, an ice intrusion remains inside the AH. This process has been witnessed in laboratory experiments, even though the intrusion of ice into the forming hydrate was not recognized at that time (Reference Uchida, Hondoh, Mae, Duval and LipenkovUchida and others, 1994b; Reference Shimada and HondohShimada and Hondoh, 2004). Driven by a thermodynamic urge to reduce the surface energy, the AH can then recrystallize to form a smooth AH without internal structuration. The connected AHs seen in Figure 6 are an example of this process, where the pAH volume is diffusively redistributed in favour of a smooth and compact sAH body. This is further supported by the fact that the N2/O2 ratios of pAH-3 and sAH-1 are identical in the given error range. The impurities visible in sAH-1 could have acted as nucleation points for the sAH (Reference Ohno, Lipenkov and HondohOhno and others, 2010). Microscopic images of more AHs undergoing such a metamorphism are provided by Reference Pauer, Kipfstuhl, Kuhs and ShojiPauer and others (2000). It is important to keep in mind that this is not the only process capable of forming smooth AHs; a direct conversion from AB to sAH has been demonstrated in laboratory experiments (Reference Ikeda, Uchida and MaeIkeda and others, 1993).

However, the shapes of sAH-3 and sAH-4 are evidently not a product of the above-mentioned processes. sAH-3 exhibits several twinned octahedra that clearly follow the cubic crystal structure of the AH, while sAH-4 closely resembles the so-called plate-like inclusions (PLIs), relaxation features that are known to form on the basal planes of ice crystals, following the hexagonal structure of ice Ih (Reference Nedelcu, Faria and KuhsNedelcu and others, 2009).

The structured pAHs with ice intrusions are more susceptible to mechanical destruction than solid AHs. The frequently seen clusters of small AHs, described by Reference Salamatin, Lipenkov, Hondoh and IkedaSalamatin and others (1999) as a product of a collapsing AH shell during transition, could also have been formed by post-genetic separation due to ice deformation.

N2/O2 ratios

With the exception of pAH-1, which has a ratio close to the atmospheric value, all AHs measured for this study are enriched in nitrogen. This is contrary to other works that found the AHs in the transition zone of Dome Fuji, Antarctica, to be oxygen-enriched (Reference IkedaIkeda and others, 1999; Reference Ikeda-Fukazawa, Hondoh, Fukumura, Fukazawa and MaeIkeda-Fukazawa and others, 2001). However, it has also been shown that ABs in the transition zone are enriched in nitrogen: Reference Weikusat, Freitag and KipfstuhlWeikusat and others (2012) reported N2/O2 values up to 6 for EDML while Reference Ikeda-Fukazawa, Hondoh, Fukumura, Fukazawa and MaeIkeda-Fukazawa and others (2001) found extreme values up to 14 in the Vostok core. As all ABs will have transformed by the end of the transition zone, it seems likely that the N2/O2 ratios measured in the AHs are due to the conversion from already nitrogen-enriched ABs. Exceptionally high N2/O2 ratios as measured in sAH-2 (8.30), sAH-4 (8.09) and sAH-3 (17.49) could also indicate a different path of hydrate formation, not necessarily implying equally nitrogen-enriched ABs. For sAH-3 and sAH-4, this is supported by the unusual morphologies, which are clearly not derived from a common air bubble. It is also worth noting that the solid AHs without ice intrusions investigated in this study show higher N2/O2 ratios than the cloudy pAHs, yet more data are needed for a statistical evaluation. Selective gas losses during the long storage time could also have altered the measured mixing ratios. A detailed interpretation of these unusually high values is outside the scope of this work and will be addressed in a dedicated publication.

Conclusion

The structured pAHs are a result of a fast conversion process, where the AH shell is not closed at any time and thus a fast conversion is possible during the whole process. Ice from the surrounding matrix flows to intrude the former AB volume, keeping the pressure in the remnant bubble above the AH dissociation pressure. The flow velocity is then an important factor for the conversion. On the microscale, the flow can be very fast (hours) and would therefore be on the same timescales as the formation processes seen in laboratory experiments (Reference Uchida, Hondoh, Mae, Duval, Lipenkov, Maeno and HondohUchida and others, 1992, Reference Uchida, Hondoh, Mae, Lipenkov and Duval1994a; Reference Ikeda, Uchida and MaeIkeda and others, 1993; Reference Shimada and HondohShimada and Hondoh, 2004). The constant contact to the ice matrix also eliminates the need for (slow) diffusion of water molecules through the AH layer.

Detailed investigations and modelling of the processes with a suitable tool (e.g. ELLE; Reference Bons, Köhn and JesselBons and others, 2008) could yield information about deformation and strain distribution on the microscale.

Acknowledgements

We thank Daniela Jansen for fruitful discussions, and three anonymous reviewers for helpful comments and suggestions. This work was supported by the Deutsche Forschungsgemeinschaft (DFG) in the framework of the priority programme ‘Antarctic Research with comparative investigations in Arctic ice areas’ by grant DFG WE4771 (CW) and by the Helmholtz Association by YIG grant VH-NG-802 (IW).

References

Abràmoff, MD, Magalhães, PJ and Ram, SJ (2004) Image processing with ImageJ. Biophoton. Int, 11, 3642 Google Scholar
Bendel, V and 6 others (2013) High-resolution variations in size, number and arrangement of air bubbles in the EPICA DML (Antarctica) ice core. J. Glaciol., 59(217), 972980 (doi: 10.3189/2013JoG12J245)CrossRefGoogle Scholar
Bender, ML, Sowers, T and Lipenkov, V (1995) On the concentrations of O2, N2, and Ar in trapped gases from ice cores. J. Geophys. Res, 100, 18 65118 660 Google Scholar
Bons, P, Köhn, D and Jessel, M (2008) Microdynamic simulation – from microprocess to patterns in rocks. (Volume 106 of Lecture Notes in Earth Sciences) Springer, Berlin Google Scholar
Chazallon, B, Champagnon, B, Panczer, G, Pauer, F, Klapproth, A and Kuhs, W (1998) Micro-Raman analysis of synthetic air clathrates. Eur. J. Mineral, 10, 11251134 Google Scholar
Durham, WB, Kirby, SH and Stern, LA (2000) Comparative rheologies of solid H2O, CO2, and CO2 clathrate hydrate. In Proceedings of the 2nd International Conference on Mars Polar Science and Exploration, 21–25 August 2000, Reykjavik, Iceland. Lunar and Planetary Institute, Houston, TX Google Scholar
Durham, W, Kirby, S, Stern, L and Zhang, W (2003) The strength and rheology of methane clathrate hydrate. J. Geophys. Res, 108 (doi: 10.1029/2002JB001872)Google Scholar
Durham, W, Stern, L, Kirby, S and Circone, S (2005) Rheological comparisons and structural imaging of sl and sll endmember gas hydrates and hydrate/sediment aggregates. In Proceedings of the 5th International Conference on Gas Hydrates, Trondheim, Norway, June 2005. Tapir Academic Press, Trondheim Google Scholar
Durham, W, Prieto-Ballesteros, O, Goldsby, DL and Kargel, JS (2010) Rheological and thermal properties of icy materials. Space Sci. Rev., 153(1–4), 273298 (doi: 10.1007/s11214-009-9619-1)Google Scholar
Duval, P, Ashby, MF and Anderman, I (1983) Rate-controlllng processes in the creep of polycrystalline ice. J. Phys. Chem, 87, 40664074 Google Scholar
Ikeda, T, Uchida, T and Mae, S (1993) The effect of hydrostatic pressure on the formation of air-hydrate crystals. In Proceedings of the NIPR Symposium on Polar Meteorology and Glaciology 7. National Institute of Polar Research, Tokyo Google Scholar
Ikeda, T, Fukazawa, H, Mae, S, Hondoh, T and Langway, CC Jr (1997) Crystal-orientation dependence of Raman spectra of natural air hydrate single crystal. J. Phys. Chem. B, 101, 61806183 Google Scholar
Ikeda, T and 7 others (1999) Extreme fractionation of gases caused by formation of clathrate hydrates in Vostok Antarctic ice. Geophys. Res. Lett, 26, 9194 Google Scholar
Ikeda-Fukazawa, T, Hondoh, T, Fukumura, T, Fukazawa, H and Mae, S (2001) Variation in N2/O2 ratio of occluded air in Dome Fuji Antarctic ice. J. Geophys. Res, 106, 17 799–17 810Google Scholar
Kipfstuhl, S, Pauer, F, Kuhs, WF and Shoji, H (2001) Air bubbles and clathrate hydrates in the transition zone of the NGRIP deep ice core. Geophys. Res. Lett, 28, 591594 Google Scholar
Kuhs, WF, Klapproth, A and Chazallon, B (2000) Chemical physics of air clathrate hydrates. In Hondoh, T ed. Physics of ice core records. Hokkaido University Press, Sapporo, 373392 Google Scholar
Lipenkov, VY (2000) Air bubbles and air-hydrate crystals in the Vostok ice core. In Hondoh, T ed. Physics of ice core records. Hokkaido University Press, Sapporo, 327358 Google Scholar
Lüthi, D and 11 others (2010) CO2 and O2/N2 variations in and just below the bubble–clathrate transformation zone of Antarctic ice cores. Earth Planet. Sci. Lett, 297, 226233 (doi: 10.1016/j.epsl. 2010.06.023)Google Scholar
Miller, SL (1973) The clathrate hydrates – their nature and occurrence. In Whalley, E, Jones, SJ and Gold, LW eds Physics and chemistry of ice: papers presented at the Symposium on the Physics and Chemistry of Ice, held in Ottawa, Canada, 14–18 August 1972. Royal Society of Canada, Ottawa Google Scholar
Nakahara, J, Shigesato, Y, Higashi, A, Hondoh, T and Langway, CC Jr (1988) Raman spectra of natural clathrates in deep ice cores. Philos. Mag. B, 57, 421430 (doi: 10.1080/13642818808208514)Google Scholar
Narita, H and 8 others (1999) Characteristics of air bubbles and hydrates in the Dome Fuji ice core, Antarctica. Ann. Glaciol, 29, 207210 CrossRefGoogle Scholar
Nedelcu, AF, Faria, SH and Kuhs, W (2009) Correspondence. Raman spectra of plate-like inclusions in the EPICA-DML (Antarctica) ice core. J. Glaciol, 55, 183184 (doi: 10.3189/002214309788609010)Google Scholar
Oerter, H, Graf, W, Meyer, H and Wilhelms, F (2004) The EPICA ice core from Dronning Maud Land: first results from stable-isotope measurements. Ann. Glaciol, 39, 307312 (doi: 10.3189/172756404781814032)Google Scholar
Ohno, H, Lipenkov, VY and Hondoh, T (2010) Formation of air clathrate hydrates in polar ice sheets: heterogeneous nucleation induced by micro-inclusions. J. Glaciol, 56, 917921 (doi: 10.3189/002214310794457317)CrossRefGoogle Scholar
Pauer, F, Kipfstuhl, S, Kuhs, WF and Shoji, H (1999) Air clathrate crystals from the GRIP deep ice core, Greenland: a number-, size- and shape-distribution study. J. Glaciol, 45, 2230 Google Scholar
Pauer, F, Kipfstuhl, S, Kuhs, WF and Shoji, H (2000) Classification of air clathrates found in polar ice sheets. Polarforschung, 66, 3138 Google Scholar
Raynaud, D, Barnola, JJ-M, Chappellaz, J, Blunier, T, Indermühle, A and Stauffer, B (2000) The ice record of greenhouse gases: a view in the context of future changes. Quat. Sci. Rev, 19, 917 Google Scholar
Salamatin, AN, Hondoh, T, Uchida, T and Lipenkov, VY (1998) Post-nucleation conversion of an air bubble to clathrate air-hydrate crystal in ice. J. Cryst. Growth, 193, 197218 Google Scholar
Salamatin, AN, Lipenkov, VY, Hondoh, T and Ikeda, T (1999) Simulated features of the air-hydrate formation process in the Antarctic ice sheet at Vostok. Ann. Glaciol, 29, 191201 Google Scholar
Schicks, JM, Erzinger, J and Ziemann, MA (2005) Raman spectra of gas hydrates – differences and analogies to ice Ih and (gas saturated) water. Spectrochim. Acta, 61, 23992403 (doi: 10.1016/j.saa.2005.02.019)Google Scholar
Schwander, J and Stauffer, B (1984) Age difference between polar ice and the air trapped in its bubbles. Nature, 311, 4547 Google Scholar
Shimada, W and Hondoh, T (2004) In situ observation of the transformation from air bubbles to air clathrate hydrate crystals using a Mizuho ice core. J. Cryst. Growth, 265, 309317 (doi: 10.1016/j.jcrysgro.2004.01.040)Google Scholar
Shoji, H and Langway, CC Jr (1982) Air hydrate inclusions in fresh ice core. Nature, 298, 548550 Google Scholar
Uchida, T, Hondoh, T, Mae, S, Duval, P and Lipenkov, VY (1992) In-situ observations of growth process of clathrate air-hydrates under hydrostatic pressure. In Maeno, N and Hondoh, T eds Physics and chemistry of ice. Hokkaido University Press, Sapporo Google Scholar
Uchida, T, Hondoh, T, Mae, S, Lipenkov, VY and Duval, P (1994a) Air-hydrate crystals in deep ice-core samples from Vostok Station, Antarctica. J. Glaciol, 40, 7986 CrossRefGoogle Scholar
Uchida, T, Hondoh, T, Mae, S, Duval, P and Lipenkov, VY (1994b) Effects of temperature and pressure on the transformation rate from air bubbles to air-hydrate crystals in ice sheets. Ann. Glaciol, 20, 143147 Google Scholar
Uchida, T, Yasuda, K, Oto, Y, Shen, R and Ohmura, R (2014) Natural supersaturation conditions needed for nucleation of air-clathrate hydrates in deep ice sheets. J. Glaciol, 60(224), 11111116 (doi: 10.3189/2014JoG13J232)Google Scholar
Visual Computing Lab ISTI – CNR (2013) Meshlab. http://meshlab.sourceforge.net/ Google Scholar
Weikusat, C, Freitag, J and Kipfstuhl, S (2012) Raman spectroscopy of gaseous inclusions in EDML ice core: first results – microbubbles. J. Glaciol, 58, 761766 (doi: 10.3189/2012JoG11J222)Google Scholar
Wojdyr, M (2010) Fityk: a general-purpose peak fitting program. J. Appl. Cryst, 43, 11261128 (doi: 10.1107/ S0021889810030499)Google Scholar
Figure 0

Fig. 1. AH from the NorthGRIP core: (a) transmitted light; (b) asymmetric illumination reveals 3-D structure.

Figure 1

Fig. 2. Differences in the Raman spectra of ice and air hydrate that were used for imaging (averages of 100 spectra with accumulation time of 0.3 s per spectrum). a.u.: arbitrary units.

Figure 2

Fig. 3. Binary slices generated from the Raman stack scan of pAH-1 (black: air hydrate; white: ice). The relative depth is given for each image; the scale is given on the photograph (lower right). See Figure 4 for the corresponding 3-D reconstruction.

Figure 3

Fig. 4. pAH-1 (depth 1005 m). (a) Microscopic image (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation (steps are an artifact of large z-step size). (c) Side view; large holes in the shell are evident. (d) Semi-transparent side view (the holes are to the left) with visible ice intrusion.

Figure 4

Fig. 5. pAH-2 (depth 1048 m). (a) Microscopic image with highlighted scan area (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation (steps are an artifact of large z-step size). (c) Side view of the section with evident ice inclusion. (d) Side view from the opposite side.

Figure 5

Fig. 6. Connected AHs pAH-3 and sAH-1 (depth 1084 m). (a) Microscopic images with focus on sAH-1 (left) and pAH-3 (right). (b) 3-D reconstruction viewed in the same orientation (top part of sAH-1 out of scan range; scale bar for (b–d)). Note the hole through pAH-3. (c) Side view reveals a connection of both AHs. (d) Semi-transparent side view with visible ice intrusion.

Figure 6

Fig. 7. sAH-2 (depth 1084 m). (a) Microscopic image. (b) 3-D reconstruction viewed in the same orientation (scale bar for (b–d)). (c) View perpendicular to (b). (d) View perpendicular to (c).

Figure 7

Fig. 8. Faceted AH (sAH-3, depth 1084 m). (a) Microscopic image (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation. (c) View perpendicular to (b). Several twinned octahedra are discernible. (d) View perpendicular to (c).

Figure 8

Fig. 9. Plate-like sAH-4 (depth 1083 m). (a) Microscopic image (scale bar for all images). (b) 3-D reconstruction viewed in the same orientation. (c) View perpendicular to (b).

Figure 9

Table 1. Averaged N2/O2 ratios and standard deviations

Figure 10

Table 2. N2/O2 ratios from averaged spectra for each slice of pAH-1 (cf. Fig. 3)