Hostname: page-component-78c5997874-j824f Total loading time: 0 Render date: 2024-11-05T05:11:44.233Z Has data issue: false hasContentIssue false

Characteristic Growth Processes of Ice Crystals on the Antarctic Ice Sheet

Published online by Cambridge University Press:  20 January 2017

Hiroshi Nishimura
Affiliation:
MTS Institute Inc., Maritime Building, 14–7, Sanbancho Chiyoda-ku, Tokyo 102, Japan
Norikazu Maeno
Affiliation:
Institute of Low Temperature Science, Hokkaido University, Nishi 8, Kita 19, Sapporo 060, Japan
Rights & Permissions [Opens in a new window]

Abstract

Characteristic growth processes were investigated by measuring cross-sectional areas of ice crystals for four 30 m snow cores drilled in Mizuho Plateau, Antarctica. Considerable difference was found in the growth rate of crystals between a temperature-gradient layer above 6 m depth and an isothermal layer below 10 m depth: the growth rate in the temperature-gradient layer was much larger than that in the isothermal layer. In the isothermal layer, temperature dependence of the growth rate K was expressed by an equation K = K0 exp(−E/RT), where R and T are the gas constant and absolute temperature respectively. The apparent activation energy E is 44.7 kJ mol−1.

On the other hand, in the temperature-gradient layer, the apparent activation energy was as small as 12 kJ mol−1: the difference was explained as due to the temperature gradient. Using the temperature profiles in snow that have been estimated from the meteorological data from several stations, the growth rates in the temperature-gradient layer were calculated. The calculated temperature dependence of the growth rates, taking into consideration vertical flux of water vapor between ice particles caused by the temperature gradient, showed good agreement with measured results. It is concluded that the growth process in the layer above 6 m depth is mainly due to vapor transport under the vertical temperature gradient.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1988

Introduction

Snow deposited on polar regions densifies and undergoes a change in structure and physical properties. During the densification, which is usually accompanied by a sintering process, ice crystals in snow grow in size. Measurements of such crystal growth have been made by many authors for different types of snow in Antarctica and Greenland: Reference Stephenson and OuraStephenson (1967) and Reference GowGow (1969) demonstrated a close relationship between the growth rate and temperature which was similar to that in a sintering process of ceramic materials under isothermal conditions.

In a surface-snow layer in polar regions, however, there are often large vertical temperature gradients, which result in a considerable amount of vertical transport of water vapor. Consequently, the crystal growth there is expected to be quite different from that which occurs under isothermal conditions. Such a discrepancy, between the growth rate of the isothermal layer and the temperature-gradient layer, has been found by Reference Nakaya, Kuroiwa and OuraNakaya and Kuroiwa (1967) at Site 2 in Greenland, and by Gow (1969) at South Pole Station in Antarctica.

Various authors (e.g. Reference Quervainde Quervain 1958, Reference AkitayaAkitaya 1974) have studied the crystal growth in the temperature-gradient layer in seasonal snow fields and laboratories, but the details of the growth processes in the temperature-gradient layer in polar regions were not investigated.

In the present study, crystal size in snow was measured in 30 m snow cores that had been drilled at four sites on Mizuho Plateau, East Antarctica, and regional difference in the crystal-growth process was examined with special reference to the temperature-gradient layer.

Coring Sites and Measuring Methods

The samples used for the measurements were the four 30 m snow cores, which were also used for the measurement of densities (Reference Nishimura, Maeno and SatowNishimura and others 1983), air permeabilities and specific areas of internal free surfaces (Reference Nishimura and MaenoNishimura and Maeno 1984, 1985). The cores were drilled at the four sites shown in Figure 1: S18 (69°02′S, 40°07′ E), W′200 (69°35′S, 48°50′E), U234 (71°01′ S, 47°29′E) and V142 (72°32′S, 51°57′E). Snow temperatures at 10 m depth are −15.9°C at S18, −33.1°C at W′200, −38.5°C at U234 and −48.1°C at V142; these are considered to be close to the mean annual air temperatures at each site. Mean annual accumulation rates of snow were estimated from stratigraphie analyses of the cores as 210kg m2 a−1 at S18, 290 kg m−2a−1 at W′200, 190 kg m−2 a−1 at U234 and 90 kg m−2 a−1 at V142.

Vertical thin sections of snow were prepared by cutting aniline-reinforced samples. Thin sections planed to thicknesses about 0.2−0.4 mm were photographed under polarized light. Mean cross-sectional areas of individual ice crystals were measured in photographs roughly 20 mm × 14 mm in area.

Fig. 1. Map of Mizuho Plateau, showing locations of core-sampling sites (S18, W′200, U234, V142), represented by open circles.

Results

Mean cross-sectional areas of crystals at the four sites are plotted against depth in Figure 2; they increase with depth, showing large rates of increase at more than 6 m depth at the four sites. The overall tendency is similar to that found at South Pole Station (Reference GowGow 1969).

Fig. 2. Mean value of cross-sectional areas of ice crystals, plotted against depth at S18, W′200, U234 and V142.

Figure 3 shows relations between the crystal areas and pressure of overburden snow at the four sites; pressures were calculated from the density-depth relations of the cores (Reference Nishimura, Maeno and SatowNishimura and others 1983). It should be noted that each relation between the crystal area and pressure is expressed roughly by two straight lines. The rates of increase are greater in shallower layers, less than about 3 × 104 Pa, compared with those in the deeper layers. The pressure of 3 × 104 Pa corresponds to depths from 6 to 10 m.

Assuming a constant accumulation rate of snow, the time that had elapsed since deposition at the surface can be estimated from the pressure of overburden snow divided by the mean annual accumulation rate. Then it was found that the crystal area increases linearly with time; the linear relationship is similar to the results obtained by Reference Stephenson and OuraStephenson (1967) and Gow (1969).

Fig. 3. Mean value of the cross-sectional areas of ice crystals, plotted against overburden pressure at the four sites.

Fig. 4. Relations between the growth rate of the cross-sectional areas of crystals and the reciprocal of temperature in the layer above 6 m depth (open circles) and below 10 m depth (solid circles).

Figure 4 shows the relation between the growth rate K and the reciprocal of temperature 1/T, where K is the proportionality constant appearing in the linear relation and T is the absolute temperature. In the deeper layer, below 10 m depth, the growth rate decreases significantly with decreasing temperature. This relation, which is almost like that of Gow (1969), is expressed as

(1)

where K 0 is constant (12.9 m2 a−1), E and R (8.314 J mol−1 K−1) are the apparent activation energy of crystal growth and the gas constant respectively. The activation energy obtained from the slope of the straight line for snow below 10 m in depth is 44.7 kj mol−1, which is slightly smaller than the value of 48.6 kJ mol−1 given by Gow (1969).

In the surface layer above about 6 m depth, the apparent activation energy is extremely small — almost four times smaller: E = 12.0 kj mol−1 and K 0 = 4.7 × 10−6 m2 a−1. This result suggests that the growth processes of crystals are different in the layers above 6 m and those below 10 m. In the deep layer, below about 10 m, it is considered to be an isothermal crystal-growth process, but in the surface layer, above 6 m, it seems to be necessary to take account of the effect of the temperature gradient, since polar snow in a surface layer above 6 m depth is usually subject to a large temperature gradient (Reference Nishimura and MaenoNishimura and Maeno 1984). In the following we try to explain the small apparent activation energy in the surface layer by the effect of the vertical temperature gradient in snow.

Thermal Condition in Snow

In Antarctica amplitudes of annual variations in air temperatures are usually greater in the interior regions than in the coastal regions. According to Nishimura and Maeno (1984), the amplitude of annual variations in air temperatures A a in Antarctica, which is defined as half the difference between the maximum and minimum of monthly mean air temperatures, is described by the mean annual air temperature θ ao as

(2)

Using this equation, A a at the four sites are estimated to be 9.1°C at S18, 12.7°C at W′200, 14.1°C at U234, and 16.9°C at V142.

Diurnal temperature variations also differ among the four sites. Like the amplitude of annual variations in air temperature, the average amplitudes of diurnal temperature variations A d are also related to the mean annual air temperature θ ao, which is expressed as the following equation:

(3)

Estimated values of A d are 3.3°C at S18, 4.3°C at W′200, 4.6°C at U234, and 5.3°C at V142.

We assume that the surface temperature θ0 of a snow layer changes with time t as

(4)

where λ is the period, and A and θ 0 are the amplitude and mean temperature of the snow surface respectively. Then the temperature at a depth z is given from a common theory of thermal conduction (Reference Carslaw and JaegerCarslaw and Jaeger 1967):

(5)

where a is the thermal diffusivity of snow.

The temperature gradient in snow at a depth z is given by differentiation of Equation (5):

(6)

Assuming that the mean annual temperature of the snow surface is identical with the mean annual air temperature, we can obtain the temperature and temperature gradient in snow by Equations (5) and (6) respectively. Summing up the profiles for the annual and diurnal variations, the temperature profiles in the snow layer are obtained as shown in Figure 5. This figure shows examples from three seasons at V142.

Fig. 5. Examples of temperature profiles for three different seasons, calculated for annual and diurnal variation.

Crystal Growth in the Temperature-Gradient Layer

To discuss the growth of ice particles in the temperature-gradient layer, we assume a simple arrangement of particles as shown in Figure 6. Cubic ice particles with a side length of a are arranged vertically with an interval of la, l being the distance between the centers of the particles; temperature gradient T′ is applied vertically, the temperature being higher at lower levels. Since ice is a

Fig. 6. Idealized particle arrangement and schematic profiles of average temperature in snow (broken line) and local temperature in the particles and spaces (solid line).

much better conductor of heat than air, the temperature gradient in ice particles is smaller than the average temperature gradient of snow T′. The temperature gradient in ice particles is described as

(7)

where k s and k i are the thermal conductivities of snow and ice respectively. The temperature gradient in air space between ice particles dT/dz, which is larger than that of snow, is expressed as

dT

(8)

where

(9)

When we assume that air in the pore space of snow is saturated by water vapor, the flux of water vapor between ice particles f is given by

(10)

where ρv is the saturated density of water vapor. D is the diffusion coefficient of water vapor given by the following equation (Reference Hall and PruppacherHall and Pruppacher 1977):

(11)

where x = 2.11 × 10−5 (P ao/P) and y = 1.94. P is air pressure, P ao = 1.013 × 105 Pa, and T 0 = 273 K.

The value of ρv on the ice surface is obtained from the Clausius-Clapeyron equation as

(12)

where P 0 is the saturated vapor pressure at temperature T 0, M is the molecular weight of water, and R and L are the gas constant and latent heat of the sublimation of ice respectively. Equation (12) can be rewritten with good approximation in the temperature region of the problem as

(13)

where A = MP 0 and

.

The gradient of vapor density dρv/dz is given by

(14)

Then the flux between ice particles is given by

(15)

The vapor density over the bottom surface of a particle is smaller than that over the upper surface of the lower neighbor. Then water vapor is transported upward according to Equation (15). The ice particle grows at the bottom and sublimes at the top. When we assume that the vapor flux between ice particles given by Equation (15) is supplied by water vapor from the upper surface of the particle and that the vapor is condensed on the bottom surface of the higher neighbor, the growth rate K of the particle for the cross-sectional area is given by

(16)

where f 1 and f 2 are the fluxes between particles in the lower and upper spaces of the particle. As the difference in temperature gradients T′ between the upper and the next lower spaces can be negligible, the growth rate is expressed as

(17)

Because B >> (y − 1)/T, Equation (17) can be simplified to

(18)

This equation shows that the growth rate K is related to temperature T, temperature gradient T′, and ala 2. ∝ is also a function of the structural parameters of snow, l and a as expressed by Equation (9). The values of a, l, k s/k i, and ala 2 at the four coring sites are given in Table I. In the present calculations, the value of a was assumed to be the length of the side of a square with the same area as a mean cross-sectional area of ice crystals in the layer above 6 m depth. The value of l, which is also the mean value in the layer above 6 m depth, is calculated by using a relation

, where ρ is snow density. The ratio k s/k i was obtained by using the value k i = 2.48 J m−1 s−1 °C−1 and the values of k s were obtained by the calculation (Reference Nishimura and MaenoNishimura and Maeno 1984). The values of alα2 at the four sites lie in the range from 1.59 mm2 to 3.10 mm2.

TABLE I. Structural parameters in the layer above 6 m depth at the four sites

Equation (18) was elevated at time steps of 1 h for the diurnal variation and 10 d for the annual one. The mean growth rates obtained for cross-sectional areas in the layer above 6 m depth at the four sites are represented in Figure

7 by open circles; measured values, represented by solid circles, are also shown. The calculated relation between the growth rates and temperature gives a smaller value for the apparent activation energy: 5.6 kJ mol−1, quite comparable to the observed value.

Fig. 7. Relation between the growth rate of the cross-sectional areas and the reciprocal of temperature in the layer above 6 m depth. Calculated mean growth rates K are represented by open circles; solid circles indicate measured results (Fig. 4).

From the theoretical calculation, it is clear that the small value of the apparent activation energy is determined mainly by the existence of temperature gradients at the four sites. Though the value of the activation energy is even smaller than the measured value, 12.0 kJ mol−1, it is suggested that the growth process in the layer above 6 m depth is mainly due to vapor transport under the vertical temperature gradient. But the calculated growth rates are smaller by about two orders of magnitude than the measured growth rates. This difference might be explained by considering the three-dimensional structure, the distribution of the sizes of ice particles and intervals between particles, and so on.

In the temperature-gradient layer, ice particles grow mainly by vertical vapor flux, which results in an anisotropic structure (i.e. a vertically elongated structure); the growth driven by surface energy in an isothermal process is rather weak. The measured results of the growth rates show that the discrepancy in the growth rates in the temperature-gradient layer (above 6 m depth) and the isothermal layer (below 10 m depth) is smaller in the regions where there are higher temperatures, and greater at lower temperatures. The growth rates in the temperature-gradient layers at U234 and V142 are considered to be affected significantly by the vertical vapor flux; this view is also supported by the observation that vertically elongated structures were found frequently at U234 and V142 (Reference Nishimura and MaenoNishimura and Maeno 1984). The above result agrees with that obtained from the study of regional differences in the development of vertical structure (Nishimura and Maeno 1984). Vertical vapor flux is considered to be a main reason for the regional differences in the development of vertical structure.

Acknowledgements

The authors would like to express their thanks to Professor S. Kobayashi of Niigata University and Dr R. Naruse of the Institute of Low Temperature Science, Hokkaido University, for their valuable discussions. The authors are also grateful to Professors G. Wakahama and Y. Suzuki and to Dr T. Kuroda of the Institute for their helpful discussions.

References

Akitaya, E. 1974 Studies on depth hoar. Inst. Low Temp. Sci. Contrib., Ser. A., 26.Google Scholar
Carslaw, H. S. Jaeger, J. C. 1967 Conduction of heat in solid Second edition. Oxford, Clarendon Press. Google Scholar
Gow, A. J. 1969 On the rates of growth of grains and crystals in south polar firn. J. Glaciol., 8(53), 241252.Google Scholar
Hall, W. D. Pruppacher, H. R. 1977 The survival of ice particles falling from cirrus clouds in subsaturated air. J. Almos. Sci., 33, 19952006.Google Scholar
Nakaya, U. Kuroiwa, D. 1967 Physical properties and internal structure of Greenland snow. In Oura, H., ed. Physics of Snow and Ice. International Conference on Low Temperature Science … 1966 … Proceedings. Vol. 1. Pt 2. Sapporo, Hokkaido University. Institute of Low Temperature Science, 953971.Google Scholar
Nishimura, H. Maeno, N. 1984 Snow structure and depth hoar formation in Mizuho Plateau,Antarctica. Mem. Natl Inst. Polar Res. Special Issue, 34, 137146.Google Scholar
Nishimura, H. Maeno, N. 1985 Studies on structures and physical properties of snow on Mizuho Plateau, Antarctica. Ann. Glaciol., 6, 105107.Google Scholar
Nishimura, H. Inoue, J. Satow, K. 1982 Meteorological data at Mizuho Station, Antarcticain 1981 JARE Data Rep., 77. (Meteorology 12.)Google Scholar
Nishimura, H. Maeno, N. Satow, K. 1983 Initial stage of densification of snow in Mizuho Plateau, Antarctica. Mem. Natl Inst. Polar Res. Special Issue, 29, 149158.Google Scholar
Quervain, M.R. de 1958 On metamorphism and hardening of snow under constant pressure and temperature gradient. IASH Publ. 46 (General Assembly of Toronto 1957 – Snow and Ice), 225239.Google Scholar
Stephenson, P. J. 1967 Some considerations of snow metamorphism in the Antarctic ice sheet in the light of ice crystal studies. In Oura, H., ed. Physics of Snow and Ice. International Conference on Low Temperature Science … 1966 … Proceedings. Vol. I, Pt 2. Sapporo, Hokkaido University. Institute of Low Temperature Science, 725740.Google Scholar
Figure 0

Fig. 1. Map of Mizuho Plateau, showing locations of core-sampling sites (S18, W′200, U234, V142), represented by open circles.

Figure 1

Fig. 2. Mean value of cross-sectional areas of ice crystals, plotted against depth at S18, W′200, U234 and V142.

Figure 2

Fig. 3. Mean value of the cross-sectional areas of ice crystals, plotted against overburden pressure at the four sites.

Figure 3

Fig. 4. Relations between the growth rate of the cross-sectional areas of crystals and the reciprocal of temperature in the layer above 6 m depth (open circles) and below 10 m depth (solid circles).

Figure 4

Fig. 5. Examples of temperature profiles for three different seasons, calculated for annual and diurnal variation.

Figure 5

Fig. 6. Idealized particle arrangement and schematic profiles of average temperature in snow (broken line) and local temperature in the particles and spaces (solid line).

Figure 6

TABLE I. Structural parameters in the layer above 6 m depth at the four sites

Figure 7

Fig. 7. Relation between the growth rate of the cross-sectional areas and the reciprocal of temperature in the layer above 6 m depth. Calculated mean growth rates K are represented by open circles; solid circles indicate measured results (Fig. 4).