Hostname: page-component-8448b6f56d-c4f8m Total loading time: 0 Render date: 2024-04-23T07:29:21.103Z Has data issue: false hasContentIssue false

Bubbles and Bubble Pressures in Antarctic Glacier Ice

Published online by Cambridge University Press:  30 January 2017

Anthony J. Gow*
Affiliation:
Research Division, Snow and Ice Branch, U.S. Army Cold Regions Research and Engineering Laboratory, Hanover, New Hampshire 03755, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

Application of the gas law to fourth-place density measurements of ice samples from two deep drill holes at “Byrd” station and “Little America V”, Antarctica, shows that virtually all density increase beyond the pore close-off density (0.830 g cm−3) can be attributed to compression of the entrapped bubbles of air. Data from “Byrd” station also indicate that the lag between overburden pressure and bubble pressure, initially 4–5 kg cm−2 at pore close-off, diminishes to less than 1.0 kg cm−2 at about 200 m depth. By substituting the overburden pressure for the bubble pressure in the pressure-density relationship based on the gas law, one can determine ice densities below 200 m more accurately than they can be measured per se on cores, because of the relaxation that occurs in samples recovered from high confining pressures. This relaxation, resulting in a progressive increase in the bulk volume of the ice with time, is generally attributed to decompression of the entrapped air bubbles following removal of the ice from high confining pressures. However. calculations of the stress in ice due to bubble pressure, together with measurements of bubble sizes in cores from various depths at “Byrd” station, both tend to indicate that there has been negligible decompression of the inclosed bubbles. It is suggested that most of this relaxation may be due to the formation of micro-cracks in the ice. Anomalous bubble pressure–density relations at “Little America V” tend to confirm abundant petrographic evidence of the existence of considerable deformation in the upper part of the Ross Ice Shelf.

Studies of crystal–bubble relations at “Byrd” station revealed that the concentration of bubbles in ice remains remarkably constant at approximately 220 bubbles/cm3. Bubbles and crystals were found to be present in approximately equal numbers at pore close-off at 64 m depth, at which level the average bubble diameter was 0·95 mm, decreasing to 0.49 mm at 116 m and to 0·33 mm at 279 m. Despite a ten-fold increase in the size of crystals between 64 and 279 m, the bubbles showed no tendency to migrate to grain boundaries during recrystallization of the ice. The observation that most of the bubbles had assumed substantially spherical shapes by 120 m depth points to essentially hydrostatic conditions in the upper layers of the ice sheet at “Byrd” station.

Des bulles ed la pression des bulles dans la glace des glaciaires Antarctiques. L′application de la loi des gaz à la mesure jusqu′à la quatrième décimale de la densité de échantillons de glace prélevés dans des trous de forage profond effectués à la station “Byrd” et à “Little America V” dans l′Antarctique, montre que les augmentations de densité au delà de la densité de fermeture des pores (0,830 g cm−3) peuvent pratiquement toutes être attribuées à la compression des bulles d′air enfermées dans la glace. Des données obtenues à la station “Byrd” indiquent également que le retard entre la pression de surcharge et la pression des bulles, qui est à l′origine de 4 à 5 kg cm−2 à la fermeture des pores, diminue jusqu′à moins de 1 kg cm−2’ à une profondeur d′environ 200 m. La substitution de la pression des bulles par la pression de surcharge dans la relation pression-densité basée sur la loi des gaz permet de déterminer les densités de glace au delà de 200 m avec plus de précision que ne le permettrait une mesure à vue sur des échantillons prélevés à cause de la relaxation qui se produit dans des échantillons tirés des situations où ils sont soumis à des pressions élevées. Cette relaxation, qui provoque une augmentation progressive du volume total de la glace avec le temps, est généralement attribuée à la décompression des bulles d′air enfermées dans la glace après retrait de l′échantillons des conditions de pression élevées auxquelles il était soumis. Toutefois, des calculs des tensions présentes dans la glace et dues à la pression des bulles, ainsi que des mesures de dimensions des bulles dans des échantillons de glace prélevés à différentes profondeurs à la station “Byrd”, indiquent que la décompression à laquelle sont soumises les bulles est négligeable. Les résultats indiquent que la relaxation peut être causé par la formation de fentes microscopiques dans la glace. Des relations anormales entre la pression des bulles et la densité obtenues à “Little America V” confirment les indications petrographiques abondantes tendant à prouver l′existence de déformations considérables dans la partie supérieure du Ross Ice Shelf.

L′étude des relations bulles–cristaux à la station “Byrd” révèle que la concentration de bulles dans la glace reste remarquablement constante à environ 220 bulles par cm3. Des bulles et des cristaux ont été trouvés en nombre à peu près égal à des pressions de fermeture des pores à une profondeur de 64 m, profondeur à laquelle le diamètre moyen des bulles était de 0,95 mm, et diminuait jusqu′à 0,49 mm à 116 m et 0,33 mm à 279 m. Malgré le décuplement de la taille des cristaux entre 64 et 279 m, les bulles ne montrèrent aucune tendance à la migration vers les limites des grains durant la recristallisation de la glace. L′observation que la majeure partie des bulles avaient acquis des formes essentiellement sphériques à partir de 120 m de profondeur montre l′existence de phénomènes surtout hydrostatiques dans les niveaux supérieurs de la nappe de glace à la station “Byrd”.

Zusammenfassng

Zusammenfassng

Blasen und Blasendruck in antarktischen Gletschereis. Die Anwendung des Gasgesetzes auf viertstellige Dichtemessungen in Eisproben, aus zwei Tiefbohrlöchern im Gebiet der ‘Byrd” Station und bei “Little America V”, Antarktika zeigt, dass praktisch alle Steigerungen der Dichtezunahme über die Porenschliessdichte (0,80 g cm−3) hinaus auf die Kompression der eingeschlossenen Luftblasen zurückgeführt werden kann. Aus den von der “Byrd” Station gewonnenen Daten geht ebenfalls hervor, dass sich der Minderbetrag zwischen Überlagerungsdruck und Blasendruck, der anfangs 4–5 kg cm−2 bei Porenschliessung beträgt, auf weniger als 1,0 kg cm−2’ in etwa 200 m Tiefe verringert. Ersetzt man in der auf dem Gasgesetz basierenden Druck-Dichte-Beziehung den Blasendruck durch den Uberlagerungsdruck, so kann man Eisdichten unterhalb einer Tiefe von 200 in mit grösserer Genauigkeit bestimmen als durch Messung an den Bohrkernen selbst, und zwar wegen der EntspannunF, die in Proben eintritt, die hohen Einschliessungsdrücken entnommen werden. Diese Entspannung, die sich in einer progressiven Zunahme des Eisvolumens der Zeit äussert, wird im allgemeinen auf die Dekompression der eingeschlossenen Luftblassen nach der Entnahme des Eises aus hohen Einschliessungsdrücken zurückgeführt. Berechnungen der durch den Blasendruck im Eis erzeugten Spannungen sowie in Verbindung damit durchgeführte Messungen der Blasengrösse in Bohrkernen aus verschiedenen Tiefen im Gebiet der “Byrd” Station scheinen jedoch darauf hinzudeuten, dass die Dekompression in den eingeschlossenen Luftblasen nur unbedeutend ist. Es wird vermutet, dass diese Entspannung in erster Linie auf der Bildung von Mikrorissen im Eis beruht. Anomale Blasendruck-Dichte-Beziehungen im Gebiet von “Little America V” scheinen die in reichen Masse vorhandenen petrographischen Beweise für das Vorhandensein einer beträchtlichen Deformation im oberen Teil des Rosseisschelfes zu bestätigen.

Auf der “Byrd” Station durchgeführte Untersuchungen der Kristall-Blasen-Beziehungen haben ergeben, dass die Konzentration der Blasen im Eis mit 220 Blasen pro cm3 bemerkenswert konstant bleibt. Blasen und Kristalle wurden in etwa gleicher Anzahl bei Porenschliessdruck in 64 m Tiefe festgestellt; auf diesem Niveau betrug der mittlere Blasendurchmesser 0,95 mm, der sich in 116 m Tiefe auf 0,49 mm und in 279 m Tiefe auf 0,33 mm verringerte. Trotz einer zehnfachen Vergrösserung der Kristalle zwischen 64 und 279 m zeigten die Blasen keine Tendenz, während der Rekristallisation des Eises zu den Korngrenzen zu wandern. Die Beobachtung, dass die meisten Blasen bei 120 m Tiefe fast Kugelgestalt angenommen hatten, deutet auf im wesentlichen hydrostatische Verhältnisse in den oberen Schichten der Eisdecke im Gebiet der “Byrd” Station hin.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1968

Introduction

Polar glacier ice is characterized by the widespread occurrence of air bubbles that are trapped initially during the process of transformation of snow to ice. The structure of consolidated snow (firn) is essentially sponge-like in that much of the pore space is intercommunicating, i.e. permeable. With increasing depth of burial, this pore space is gradually eliminated until the stage is reached where the snow is no longer permeable and the firn is said to have transformed into glacier ice. The freshly formed ice still contains 10 per cent by volume of air trapped in the form of discrete bubbles which become compressed as they sink deeper into the ice sheet. It is this compression of bubbles in the ice that is primarily responsible for additional increase in the bulk density of the ice.

An exceptional opportunity to study changes both in the pressure and in the distribution and structure of bubbles as a function of depth in the Antarctic ice sheet arose during the I.G.Y. (1957–58) when ice cores were obtained to depths of 308and 254 m at “Byrd” station and “Little America V”, respectively. Approximately Ioo samples of ice of density >0.830 g cm−3 from both locations were used to obtain 4–place density measurements, and these data, in conjunction with observations of bubble structure in thin sections of ice, have been used to accomplish the following objectives:

  • i. To investigate the relationship between bubble pressure and the in situ overburden pressure, and to determine, if possible, the depth at which the difference in pressures is reduced to less than 1–2 kg cm−2.

  • ii. To investigate the pattern of changes in the size, shape, number and distribution of air bubbles as a function of depth in the Antarctic ice sheet.

Previous studies

Studies of bubbles in glacier ice, including measurements of bubble pressure, have been made by a number of workers including Reference Koch and WegenerKoch and Wegener (1930), Reference BaderBader (1950), Reference ShumskiyShumskiy (1955), Reference CoachmanCoachman and others (1956), Reference ScholanderScholander and others (1956), Reference LangwayLangway (1958), Reference SchyttSchytt (1958), Reference Scholander and NuttScholander and Nutt (1960), and Reference Matsuo and MiyakeMatsuo and Miyake (1966). Pressure measurements obtained by Schytt at Maudheim, and by Matsuo and Miyake at Syowa base, appear to be the only ones made so far on glacier ice from Antarctica. Schytt reasoned that, if all density increase beyond the pore close-off depth, i.e. for densities >0.830 g cm−3, can be attributed to bubble compression (compressibility of ice is negligibly small in the upper 200–300 m of the ice sheet), then the bubble pressure P of any particular ice sample can be expressed as P = P c V c/V, where P c and V c are the bubble pressure and void ratio respectively at pore close-off and V is the void ratio of the sample. Subsequent tests of this relationship by Schytt with samples from between 90.5 and 97 m at Maudheim revealed very satisfactory agreement between the measured and calculated pressures.

Direct measurements of bubble pressure, i.e. measurements of pressure within individual bubbles in glacier ice, have been made by Reference Scholander and NuttScholander and Nutt (1960) on samples from a number of Greenland icebergs. Most bubbles were found to have pressures of 4–5 atmospheres, though pressures as high as 20 atmospheres were recorded in some bubbles. However, most of the pressure measurements reported so far have been obtained indirectly from very precise determinations of the density of bubbly ice made in conjunction with measurements of the volume of air released upon melting the samples. Of these measurements, those made by Reference LangwayLangway (1958, Reference Langway1962) on core samples from “Site 2” in Greenland are probably the most accurate obtained so far on polar glacier ice. Langway reported pressures as high as 27 atmospheres in ice from 360·5 m depth but the major point of interest of his study was that in most of the samples tested the bubbles were seemingly compressed less than might have been expected from the weight of overlying snow and ice. Other conclusions reached by Langway as a result of this work were that the composition of the inclosed air in ice at “Site 2” has remained constant for the last goo years and that the mass of entrapped air also has remained unchanged and does not diffuse through the ice.

Analytical Procedures

In the present study no attempt was made to measure bubble pressures per se. However, it should be possible to calculate the average bubble pressure in a sample of ice if the density of that sample has been accurately determined, and if there are reasonably accurate data on the pore close-off density and pressure, and the variation of overburden pressure with depth.

If the density at pore close-off is denoted by γ c Footnote * and the pressure in the bubbles at this depth is P c (this is essentially equal to the average station pressure corrected for the depth at which the pores finally close off), then the average bubble pressure P is:

(1)

where γ it is density of bubble-free ice at the in situ temperature and γ is the densityof the sample at the same temperature. This is, of course, equivalent to the expression P = P c V c/V with the void ratios simply re-written in terms of density. Relationship (1) can be readily transformed to evaluate γ and γ c if the various other parameters are known, e.g.

(2)

(3)

The original densities were determined by hydrostatic weighing of ice cores using highly purified isooctane of density 0·7160 g cm−3 at –10°C as the immersion liquid. These data are tabulated fully by Reference GowGow (in press) and are accurate to at least ± 0·00015 g cm−3.

All densities were corrected for in situ temperatures using the simplified expression of Reference ButkovichButkovich (1957) for correcting for isopeizic density of ice between 0° and −30°C. At “Byrd” station the in situ temperature of the ice remains essentially constant at −28°C and the isopeizic density, γ it, of bubble-free ice at this temperature is 0.9204 g cm−3. At “Little America V”, however, temperatures in the ice shelf are not constant but increase from −21°C at 52 m to −2·0°C at the bottom (depth 258 m), causing a corresponding decrease of γ it from 0.9195 to 0.9168 g cm−3.

Results and Discussion

“Byrd” station

Rather than use relationship (1) to compute pressures from measured densities, relationship (2) was used instead to solve for γ over a range of pressures of from 1 to 25 atmospheres. The calculated densities were matched with depth–density data from cores and the ice loads (σ) appropriate to each density were then calculated. The absolute load was obtained by merely adding the average station pressure (0.82 kg cm−2) to the calculated ice loads. According to Langway’s bubble-pressure data at “Site 2”, pore close-off occurs at approximately 0.828 g cm−3. Observations on permeabilities of thin slices of core from “Byrd” station also indicated that the pores had closed off at much the same value, but in the absence of more precise data it was decided to evaluate (2) for three separate values of γ c, namely 0.820, 0.825 and 0.835 g cm−3. Pressure–density data for γ c = 0.825 g cm−3 are given in Table I to indicate the method of calculation and also to give the reader some idea of the average bubble pressure associated with any particular density. Curves of bubble pressure and pressure lag (ΔP) versus absolute load are presented in Figures 1 and 2.

Table I. Bubble-Pressure–Density Relationships in Ice Cores from “Byrd” Station, Antarctica. Bubble–pressure Calculations Based on Pore Close-off Density γ c = 0.825 g cm −3

Fig. 1. Bubble pressure versus absolute-load pressure at “Byrd” station, Antarctica. The broken line represents the probable trend in in situ ice below 200 m. Divergence at a σ+Pc>14−16 kg cm−2 due to expanded (relaxed) condition of samples

Fig. 2. Bubble-pressure lag (ΔP) versus absolute-load pressure at “Byrd” station. Data from Figure 1

Data obtained with P c = 0.820 g cm−3 are particularly instructive in that the computed bubble pressures actually begin to exceed in situ pressures at a depth of 165 m. This indicates immediately that porc close-off at ‘Byrd” station must generally occur at densities >0.820 g cm−3. A further feature of the data is the sudden and marked divergence of ΔP at absolute-load pressures >15 kg cm−2, i.e. deeper than 175 m. This can be due only to an actual expansion of the ice which, in the case of the very deepest samples, is caused by fracturing of the ice during drilling. However, some cores from between 175 and 240 m contained no detectable cracks whatsoever. These samples had “relaxed” to some degree according to density measurements, and the question is immediately raised as to what the mechanism of relaxation is in these particular samples. Langway has attributed identical relaxation of uncracked ice from “Site 2” in Greenland to an actual decompression of the entrapped bubbles (Fig. 3). An inspection of Langway’s data for samples tested immediately after drilling and later at intervals of 12 and 18 months indicates that this decompression is confined to samples with bubble pressures >10 kg cm−2. Less compressed bubbles apparently retain their pressures. In support of his decompression mechanism, Langway cited evidence of the increasing occurrence of small pressure cracks around bubbles in ice from depths >150 m. Some cracking of bubbles was also observed in the deeper samples at “Byrd” station and at the time it was believed that decompression of cracked bubbles was also the principal mode of relaxation of ice recovered from high confining pressures at “Byrd” station. However, subsequent examination of numerous thin sections of ice from below 170 m at “Byrd” station showed that most of the bubbles were not intersected by pressure cracks. Furthermore, actual measurements of the volumes of individual bubbles at various depths in the ice indicated that they were entirely commensurate with in situ pressures. In these samples at least it would appear that relaxation of the ice did not involve any significant decompression of bubbles. The evidence for this will be considered in a later section of this paper dealing with bubble structures in ice at “Byrd” station.

Fig. 3. Bubble pressure versus load pressure at “Site 2”, Greenland (from Langway, 1958). Broken line represents probable trend in situ ice below 200 m

Intuitively one would expect ΔP to diminish progressively with depth and for the bubble pressure to asymptotically approach the in situ pressure. This trend is certainly indicated in results of calculations at “Byrd” station where, for γ c = 0.825 g cm−3, ΔP diminishes to less than 1.0 kg cm−2 at 200 m depth. This result is encouraging in that it seems to validate the idea that densification beyond pore close-off occurs essentially as the result of compression of entrapped bubbles. A further test of the validity of the pressure–density relationship (1) is provided by the pressure-density data from “Site 2”, Greenland. Reference BaderBader (1965) has given the following parameters for “Site 2”.

Results of computations for some representative pressure-density values of samples tested by Langway immediately after cores were pulled are presented in Table II. From these data it can be seen that measured and computed pressures agree remarkably well and that ΔP diminishes progressively with depth except for the very deepest samples which, as is clearly demonstrated in Figure 3, had already begun to relax before the cores could be pulled to the surface and tested. This increase in ΔP below 261 m is almost certainly due to a volume expansion of samples, though an under-evaluation of as little as 0.003 g cm−3 in the original density measurement could account for most of the discrepancy between the calculated bubble pressures in the deepest samples and their in situ pressures.

Table II. Bubble-Pressure–Load-Pressure Relations in Ice Cores from “Site. 2”, Greenland. Data from Landway (1958, 1962) and Reference BaderBader (1965)

Because of the unavoidable relaxation of ice from depths greater than 200 m, it becomes increasingly difficult to obtain accurate measurements of in situ density. However, since the difference in pressure between bubbles and overburden is substantially eliminated at about 200 m, one can, by simply substituting the value of (σ+ Pc) for P in the pressure-density relationship

compute densities that are considerably more accurate than could ever be measured on deep cores per se. With due allowance for the compressibility of ice, one needs only to know the temperature distribution in an ice sheet in order to obtain very precise values of in situ density. For example, the density of all ice below 400 m at “Byrd” station can be shown to exceed 0·9180 g cm−3, the exact value depending on the temperature and in situ pressure. Ice at the bottom of the ice sheet at “Byrd” station (2500 m) would have a density of 0.9183 g cm−3, if it was at the pressure melting point (−2·0°C). This would be a minimum value as any decrease in the temperature of the bottom ice could only increase its density.

“Little America V”

Progressively increasing temperatures in the Ross Ice Shelf cause a corresponding decrease in γ it. In this case it is simpler to compute P for a range of densities using the englacial temperature corresponding to each selected density to determine γ it. For example, the density at 83 m is 0.9040 g cm−3. The temperature at this depth is −21°C and the corresponding value of γ it would be 0.9194 g cm−3. The average tation pressure at “Little America V” (elevation 43 m), corrected for the depth at which the pores finally close off (52 m), is 1.00 kg cm−2. Pressure–density relations were evaluated for γ c = 0.825 and 0.835 g cm−3, respectively (Fig. 4).

Fig. 4. Bubble pressure versus absolute-load pressure at “Little America V”, Ross Ice Shelf, Antarctica

These data tend to confirm abundant stratigraphie and petrographic evidence of the existence of considerable deformation in the Ross Ice Shelf (Reference GowGow, 1963). It is most unlikely that bubble pressures could ever exceed in situ pressures under normal hydrostatic conditions yet, even for γ c = 0.835 g cm−3, ΔP is all but eliminated at 118 m depth at “Little America V”. For γ c =0.825 g cm−3, the computed bubble pressures actually begin to exceed the in situ pressures at less than zoo m. It should be noted that an exceptional increase in firn densification was also observed below 36 m at “Little America V” and apparently this anomalous behavior has continued to some depth below the firn–ice transition.

Bubbles in Ice

Petrographic studies of Antarctic cores included the preparation of numerous thin sections, several of which were examined in some detail to determine changes in the size, shape and distribution of bubbles as a function of depth in the ice sheet. These studies were conducted primarily on samples from ‘Byrd” station but some comparative studies were also made of thin sections of cores from “Little America V”. The gross differences in crystal–bubble relations at these two locations are clearly demonstrated in Figures 5 and 6.

Fig. 5. Bubble—crystal relations in ice from “Byrd” station. Bubbles are the small gray and black inclusions located within crystals and at the intersections of crystals. The grid spacing is 1 cm.

a. 70 m depth; b. 116 m depth; c. 163 m depth; d. 305 m depth

Fig. 6. Bubble—crystal relations in ice from “Little America V”. Note the orientated elongate bubbles at 116 depth. The grid spacing is 1 cm.

a. 52 m depth; b. 116 m depth; c. 240 m depth

General properties

In newly formed ice at about 64. m depth at “Byrd” station the majority of bubbles are located at the intersections of crystals. Bubbles tend to be somewhat irregular at this stage (Fig. 5a) but as the depth of burial increases the bubbles become smaller and spherical. At the same time there is a slow but progressive increase in the size of crystals composing the ice. This growth is obviously to do with grain-boundary migration, but it is of considerable interest to note that during this process there is little if any tendency for the entrapped bubbles of air to be swept towards the grain boundaries. The observation that grain growth is not accompanied by any preferred orientation of the crystals, coupled with the fact that substantially spherical bubbles only are formed at depth, is a strong indication that stress conditions in the ice at “Byrd” station are essentially hydrostatic.

The number of bubbles in a unit volume of ice appears to remain constant which means that as the size of crystals increases so does the number of bubbles within a crystal. A count of the average number of bubbles in crystals in sections of given thickness from various depths in the ice thus affords a relatively rapid yet reliable means of estimating the rate of growth of crystals in the ice sheet.

Whereas crystal—bubble patterns in the ice at “Byrd” station appear to have evolved under essentially hydrostatic conditions, those at “Little America V” have not. Here, englacial deformation has caused an appreciable orientating of bubbles as well as crystals. Elongate and “waisted” bubbles are especially well developed in ice from between 100 and 130 m (Fig. 6b) and crystals generally are much larger than those at comparable depths at “Byrd” station. Numerous highly compressed spherical bubbles characterize ice in the deepest cores, but as can be observed in Figures 5d and 6c only a relatively small proportion of these bubbles is actually intersected by pressure cracks.

Bubble dimensions

Measurements of bubble dimensions were confined to sections of ice from “Byrd” station. These measurements were made on photomicrographs rather than on the thin sections themselves and, in order to further simplify the measurements, only those sections in which most of the bubbles had become spherical were used, i.e. sections from 116 m and deeper. All sections examined contained several hundred bubbles but only those bubbles located close to the upper surface of a section were used for measuring diameters. Fifty bubbles from each of five selected depths were measured and the results of these determinations are presented in Table III. From these data it can be seen that the average bubble diameter decreased from 0.49 mm at 116 m to 0.33 mm at 279 m. This corresponds to a decrease in the volume of individual bubbles of from 0.062 to 0.019 mm3.

Table III. Bubble Dimensions in Glacier Ice from “Byrd” Station, Antarctica

Since the porosity of a sample of ice can be readily ascertained from its density, these data on bubble volumes can also be used to determine the average number of bubbles present in a unit volume of ice. Results of these computations (Table III) show that the concentration of bubbles at various depths in the ice sheet at “Byrd” station remains remarkably constant at between 200 and 240 bubbles/cm3—further proof too that there is little if any tendency for entrapped air bubbles to migrate or coalesce during recrystallization of the ice.

The above data can also be used to obtain a reliable estimate of the average size of the equivalent spherical bubble at pore close-off. Because of the irregular shapes of bubbles at this depth, it is not practical to attempt to determine bubble dimensions directly. Since the porosity at close-off is ≈10 per cent, 1 cm3 of sample will contain 100 mm3 of entrapped air. Using a value of 220 bubbles/cm3 (the average count obtained from measurements on the five samples from 116 m and deeper), it is found that the mean volume of the equivalent spherical bubble at pore close-off at “Byrd” station is 0.45 mm3. Its diameter would be 0.95 mm.

Bubble–crystal relations

With data already obtained on crystal-size variations with depth at “Byrd” station (Reference Gow and BentleyGow, 1964), the relationship between the number of bubbles and crystals at pore close-off can now be determined. Average crystal cross-sections between 61 and 66 m at “Byrd” station varied from 3·1 to 3·7 mm2. If a value of 3·4 mm2 is taken and, it is assumed for simplicity, that the sections are hemispherical, then the volume of the corresponding spherical crystal will be 4·71 mm2. If crystals of these dimensions could be deformed-packed into the remaining 900 mm3 of space in 1 cm3 of ice, then the number of crystals at close-off would be 900/4·71 = 191. The number of bubbles in the same 1 cm3 of ice is ≈220, from which it can be concluded that at pore close-off there are approximately as many crystals present as bubbles. One interesting aspect of this result is that it does tend to demonstrate that the 1 : 1 relationship between voids and particles observed in systematically packed aggregates of equidimensional particles also holds for naturally dense-packed structures. As it turns out the crystal-pore structure at the firn–ice transition probably conforms as closely to a dense-packed aggregate as any in nature.

A further implication of the above result is that any change in the number of crystals per cubic centimeter at the firn–ice transition should be accompanied by a corresponding change in the number and size of bubbles. For example, fire crystal-growth studies to a depth of 50 m at the South Pole (Reference GowGow, unpublished) have shown that the mean crystal cross-section increases essentially linearly with depth at the rate of 0.005 mm2 m−1. At 50 m the mean cross-sectional area of crystals was 0.65 mm2 which, if extrapolated to the depth of the firn-ice transition 110 m (estimated) gives a mean crystal cross-section of 0.95 mm2. 1 cm3 of ice at 110 m would thus contain nearly 1400 crystals. This would represent a seven-fold increase in the number of crystals and bubbles over that observed at “Byrd” station but, since the total available pore space is still the same (≈10 per cent), the bubbles at pore close-off at the South Pole should be much smaller, of the order 0.5 mm in diameter compared with 0.95 mm for the diameter of the equivalent spherical bubble at “Byrd” station. On the other hand, estimates of bubble dimensions based on crystal-size data at pore close-off at a number of other locations in Antarctica indicate that at “Wilkes S-2” the number of bubbles should be much the same as at “Byrd” station (≈200 bubbles/cm3), increasing to about 300 bubbles/cm3 at Maudheim and to 350 bubbles/cm3 at “Southice”. Of course, the size and number of crystals and bubbles at pore close-off will depend upon the rate of densification, which in turn depends on the rate of accumulation and the mean annual temperature. In general, however, it would appear that the bubbliness of ice as determined by the number of bubbles/cm3 should increase towards the center of the Antarctic ice sheet.

Relaxation of ice

Measurements of bubble volume also provide an independent check on whether bubbles in the deepest ice have undergone any significant decompression. The original thin sections were prepared from samples that had lain in storage for more than a year—ample time for measurable relaxation to have occurred prior to sectioning. Since there is no evidence that bubbles with pressures <10 kg cm−2 undergo any significant decompression (Langway, 1958), pressure volume data for ice from either 116 or 135 m can be used to evaluate average bubble pressures in samples of deeper ice from “Byrd” station. Computations of bubble pressure P′ based on data from 116 m are presented in Table IV. The two pressure values at each depth agree remarkably well. All indications from the data in Table IV are that no significant decompression of entrapped air bubbles has occurred in the cores from “Byrd” station. Pressure cracking of bubbles could possibly occur in cores of very deep ice and, in view of the present attempt to drill to the bottom of the Antarctic ice sheet at “Byrd” station, this possibility seems worth investigating. The separation distance between neighboring bubbles in the ice as well as the stress associated with a particular bubble pressure both need to be considered. This stress in ice tending to split bubbles apart might be calculated as follows: Let the volume V of some bubble at pressure P be

, where r is the radius of the bubble. Now let this bubble be surrounded by a sphere of ice of radius R such that the average distance between neighboring bubbles is 2R (Fig. 7). The total force trying to split apart the two hemispheres of ice defined by the dashed line is F = πr 2 P.

Table IV. Bubble Volume-Pressure Relationships in Glacier Ice from “Byrd” Station, Antarctica

Fig. 7. Schematic diagram of two pressurized bubbles of radius r in glacier ice. Distance between bubble centers is 2R

The area of the interface across which this force acts is

Accordingly,

If the rupture stress of the ice is S = F /A, then the maximum allowable bubble pressure is

According to the data of Reference ButkovichButkovich (1959), the tensile strength of glacier ice of density >0.900 g cm−3 generally lies between 20 and 30 kg cm−2. If the lower value is taken, then

Before the above relationship can be evaluated it is necessary to know something of the manner of distribution of bubbles in the ice at “Byrd” station. A fairly uniform distribution of bubbles seems reasonable, at least on the basis of observations at pore close-off where the vast majority of bubbles were observed to be located at the intersections of essentially equidimensional crystals. The very simplest distribution would involve an arrangement of layers in which the distance between bubbles in the same layer is the same as the spacing between layers. In the particular case of 220 bubbles/cm3 (average concentration of bubbles observed at “Byrd” station), this spacing would be

, i.e. approximately six layers of 36 bubbles each in 1 cm3 of ice. At pore close-off, where the mean diameter (2r) of bubbles is 0.95 mm, R/r is 1.7, increasing to 5.1 at 279 m. These and intermediate values at “Byrd” station, together with the calculated values of P = 20[(R/r)2−1], are presented in Table V. From the data listed in Table V it would appear that the stress in ice due to bubble pressure is much too low to cause the ice to fracture.Footnote * In practice, of course, some small fraction of bubbles would be situated at less than the critical distance, in which case fracturing of the intervening ice could be expected to occur. The values R/r (critical) needed to attain this condition at various depths in the ice at “Byrd” station are also listed in Table V. It should be noted, however, that when bubbles are intersected by pressure cracks these cracks invariably occur parallel to the basal planes of crystals. Such a favored orientation demonstrates further the inherent weakness within the basal plane and it is likely that the rupture stress in this plane is appreciably smaller than in other directions in ice. Nevertheless the rupture stress would need to decrease by more than an order of magnitude for any widespread fracturing to occur under these conditions. In the case of those deep-core samples which are known to have undergone measurable expansion, but which have done so without any detectable cracking, it may be necessary to invoke some process of micro-fracturing to account for the expansion. Cracks formed in this way might be so small as to pass unnoticed, though it may only require a sharp knock or some other source of mechanical shock, e.g. undampened vibrations in the drill string for such micro-cracks to propagate into large visible fractures.

Table V. Estimate of Maximum Allowable Pressure (P ) in Bubbles in Glacier Ice at “Byrd” Station. Pressures Calculated According to the Relation P = 20[(R/r)2−1]. See Text for Derivation

Up to now the feeling has been that ice from deep within an ice sheet would literally explode in the core barrel because of rapid decompression of the entrapped air bubbles. However, from the results given above, it would now appear that with the proper precautions it should be possible to obtain usable cores to considerable depths in the Antarctic ice sheet. This would be best accomplished in a liquid-loaded hole such as the one that is now being drilled at “Byrd” station and which hopefully will penetrate to the bottom of the ice sheet. In this type of deep drilling it is necessary to load the hole with liquid of approximately the same density as ice so as to prevent the walls of the hole from closing. However, there are additional advantages to this system. A liquid-loaded hole dampens out virtually all vibrations in the drill string thus reducing the risk of mechanical breakage of ice cores during drilling. More important perhaps is the fact that the liquid load maintains pressure on the ice during the actual process of coring thus preventing any immediate relaxation of the ice. Since preparing the first draft of this manuscript a U.S. Army CRREL drill team has successfully cored to the bottom (1 390 m) of the Greenland ice sheet at “Camp Century” (Reference Hansen and LangwayHansen and Langway, 1966). A preliminary inspection of cores from near the bottom of the hole (liquid loaded) showed little if any fracturing of the ice and there were no immediate signs of any significant decompression of entrapped air bubbles even though in situ pressures near the bottom of the ice sheet at “Camp Century” probably exceed 100 kg cm−2.

Conclusions

Densification of bubbly ice after pore close-off is due almost entirely to compression of the entrapped air bubbles. This relationship between the bulk density γ of the ice and the pressure P of its inclosed bubbles can be expressed as follows:

where γ it is the density of bubble-free ice at the englacial temperature and P c and y c are the bubble pressure and density respectively of the ice at pore close-off. Computations of bubble pressure in ice cores from “Byrd” station, Antarctica, indicate that the pressure lag ΔP between bubbles and the absolute load, initially 4–5 kg cm−2 at pore close-off (64 m depth), diminishes to less than 1.0 kg cm−2 at 200 m. A direct consequence of this result is that by substituting the absolute-load pressure, σ+P c for P in the above relationship one can determine ice densities below 200 m more accurately than can be obtained with measurements per se on the cores because of the unavoidable relaxation of the latter following their removal from high confining pressures. This relaxation of the ice is generally attributed to decompression of the entrapped air bubbles. However, calculations of the stresses in the ice due to bubble pressure together with actual measurements of bubble volumes in “relaxed” cores from “Byrd” station tend to indicate that there is little if any significant decompression of bubbles even in ice recovered from very high confining pressures. Apparently the observed increase in bulk volume of ice samples with time does not depend upon prior expansion of bubbles. These data point up the need for further investigations of the mechanics of relaxation of ice recovered from high confining pressures.

Bubble pressure-density relations in ice cores at “Little America V” indicate anomalously high bubble pressures in the upper layers of the Ross Ice Shelf. These data tend to confirm the abundant stratigraphic and petrographic evidence for the existence of very considerable deformation near the seaward edge of the Ross Ice Shelf.

Antarctic glacier ice is characterized by an abundance of air bubbles. Measurements on ice from “Byrd” station show that the concentration of bubbles remains remarkably constant at 200–240 bubbles/cm3 of ice. At the pore close-off depth at 64 m the ice was observed to contain approximately the same number of crystals as bubbles. The mean diameter of air bubbles at this depth was 0.95 mm, decreasing to 0.49 mm at 116 m and to 0.33 mm at 279 m. A very uniform distribution of bubbles was observed at all depths. Bubbles showed no tendency to be swept towards grain boundaries during recrystallization of the ice. These observations and the fact that the majority of bubbles had become spherical at a depth of less than 120 m would also tend to indicate the existence of essentially hydrostatic conditions in the upper part of the ice sheet at “Byrd” station. By contrast, the formation of numerous elongate bubbles orientated symmetrically with respect to well-formed fabrics in ice between 90 and 130 m at “Little America V” is compatible only with the existence of considerable shearing in the ice shelf.

Acknowledgements

This research was accomplished in conjunction with a U.S. Army Cold Regions Research and Engineering Laboratory research project devoted to deep core-drilling studies in Antarctica. The project was sponsored by the U.S. Antarctic Research Program and supported by grants from the National Science Foundation. The support of these institutions, together with the assistance in the field of Task Force 43, U.S. Navy, is gratefully acknowledged. The author is indebted to Mr B. Lyle Hansen and Dr Chester C. Langway, Jr., for helpful discussions and for critically reviewing the manuscript.

Footnotes

* Bubble close-off and the firn—ice transition occur at essentially the same depth since impermeable firm is bubbly ice by definition. This transformation is observed to occur at densities of 0.82–0.84 g cm−3. It might be noted that this definition of ice is true only for samples of the dimensions of thin sections, Le. of one or two grain diameters at most. In larger samples of firn, zero permeability may effectively be reached at somewhat lower densities, in which case some compression of bubbles may have already occurred before all pore spaces had become sealed off.

* A similar conclusion was reached using the formula for so-called thick-walled spherical vessels given by Reference RoarkRoark (1954), though it is questionable whether this formula is valid for the very thick shells envisaged here.

References

Bader, H. 1950. The significance of air bubbles in glacier ice. Journal of Glaciology, Vol. 1, No. 8, p. 44351.Google Scholar
Bader, H. 1965. Theory of densification of dry, bubbly glacier ice. U.S. Cold Regions Research and Engineering Laboratory. Research Report 141.Google Scholar
Butkovich, T. R. 1957. Linear thermal expansion of ice. U.S. Snow, Ice and Permafrost Research Establishment. Research Report 40.Google Scholar
Butkovich, T. R. 1959. Some physical properties of ice from the Tuto Tunnel and Ramp, Thule. Greenland. U.S. Snow, Ice and Permafrost Research Establishment. Research Report 47.Google Scholar
Coachman, L. K., and others. 1956. Gas enclosures in a temperate glacier, by L. K. Coachman, E. Hemmingsen and P. F. Scholander. Tellus, Vol. 8, No. 4, p. 4153.Google Scholar
Gow, A. J. 1963. The inner structure of the Ross Ice Shelf at Little America V, Antarctica, as revealed by deep core drilling. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Assemblée générale de Berkeley, 19–.8–31–8 1963. Commission des Neiges et des Glaces, p. 27284.Google Scholar
Gow, A. J. 1964. Profile studies. (In Bentley, C. R., and others. Physical characteristics of the Antarctic ice sheet, by C. R. Bentley, R. L. Cameron, C. Bull, K. Kojima and A. J. Gow. Antarctic Map Folio Series (New York, American Geographical Society), Folio 2, p. 46.)Google Scholar
Gow, A. J. In press. Deep core studies of the accumulation and densification of snow at Byrd station and Little America V, Antarctica. U.S. Cold Regions Research and Engineering Laboratory. Research Report 197.Google Scholar
Gow, A. J. Unpublished. On the rates of growth of grains and crystals in South Polar firn.Google Scholar
Hansen, B. L. Langway, C. C. jr. 1966. Deep core drilling in ice and core analyses at Camp Century, Greenland, 1961–1966. Antarctic Journal of the United States, Vol. 1, No. 5, p. 20708.Google Scholar
Koch, J. P. Wegener, A. 1930. Wissenschaftliche Ergebnisse der dänischen Expedition nach Dronning Louises-Land und quer über das Inlandeis von Nordgrönland, 1912–13. Meddelelser om Grenland. Bd. 75.Google Scholar
Langway, C. C. jr. 1958. Bubble pressures in Greenland glacier ice. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Symposium de Chamonix, 16–24 sept. 1958, p. 33649.Google Scholar
Langway, C. C. jr. 1962. Some physical and chemical investigations of a 411 meter deep Greenland ice core and their relationship to accumulation. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Commission des Neiges et Glaces. Colloque d’Obergurgl, 10–9–18–9 1962, p. 10118.Google Scholar
Matsuo, S. Miyake, Y. 1966. Gas composition in ice samples from Antarctica. Journal of Geophysical Research, Vol. 71, No. 22, p. 523541.Google Scholar
Roark, R. J. 1954. Formulas for stress and strain. New York, McGraw-Hill.Google Scholar
Scholander, P. F. Nutt, D. C. 1960. Bubble pressure in Greenland icebergs. Journal of Glaciology. Vol. 3, No. 28, p. 67178.CrossRefGoogle Scholar
Scholander, P. F., and others. 1956. Gases in icebergs, by P. F. Scholander, J. W. Kanwisher and D. C. Nutt. Science, Vol. 123, No. 3185, p. 10405.CrossRefGoogle Scholar
Schytt, V. 1958. Glaciology. II. A. Snow studies at Maudheim.-B. Snow studies inland.-C. The inner structure of the ice shelf at Maudheim as shown by core drilling. Norwegian-British-Swedish Antarctic Expedition, 1949–52. Scientific Results, Vol. 4.Google Scholar
Shumskiy, P. A. 1955. Osnovy strukturnogo ledovedeniya. Petrografya presnogo l’da kak metod glyatsiologicheskogo issledovaniya. Moscow, lzdatel’stvo Akademii Nauk SSSR. [English translation; Principles of structural glaciology: the petrogrably of fresh-water ice as a method of glaciological investigation. Translated from the Russian by David Kraus. New York, Dover Publications, 1964.]Google Scholar
Figure 0

Table I. Bubble-Pressure–Density Relationships in Ice Cores from “Byrd” Station, Antarctica. Bubble–pressure Calculations Based on Pore Close-off Density γc = 0.825 g cm −3

Figure 1

Fig. 1. Bubble pressure versus absolute-load pressure at “Byrd” station, Antarctica. The broken line represents the probable trend in in situ ice below 200 m. Divergence at a σ+Pc>14−16 kg cm−2 due to expanded (relaxed) condition of samples

Figure 2

Fig. 2. Bubble-pressure lag (ΔP) versus absolute-load pressure at “Byrd” station. Data from Figure 1

Figure 3

Fig. 3. Bubble pressure versus load pressure at “Site 2”, Greenland (from Langway, 1958). Broken line represents probable trend in situ ice below 200 m

Figure 4

Table II. Bubble-Pressure–Load-Pressure Relations in Ice Cores from “Site. 2”, Greenland. Data from Landway (1958, 1962) and Bader (1965)

Figure 5

Fig. 4. Bubble pressure versus absolute-load pressure at “Little America V”, Ross Ice Shelf, Antarctica

Figure 6

Fig. 5. Bubble—crystal relations in ice from “Byrd” station. Bubbles are the small gray and black inclusions located within crystals and at the intersections of crystals. The grid spacing is 1 cm.a. 70 m depth; b. 116 m depth; c. 163 m depth; d. 305 m depth

Figure 7

Fig. 6. Bubble—crystal relations in ice from “Little America V”. Note the orientated elongate bubbles at 116 depth. The grid spacing is 1 cm.a. 52 m depth; b. 116 m depth; c. 240 m depth

Figure 8

Table III. Bubble Dimensions in Glacier Ice from “Byrd” Station, Antarctica

Figure 9

Table IV. Bubble Volume-Pressure Relationships in Glacier Ice from “Byrd” Station, Antarctica

Figure 10

Fig. 7. Schematic diagram of two pressurized bubbles of radius r in glacier ice. Distance between bubble centers is 2R

Figure 11

Table V. Estimate of Maximum Allowable Pressure (P) in Bubbles in Glacier Ice at “Byrd” Station. Pressures Calculated According to the Relation P = 20[(R/r)2−1]. See Text for Derivation