Introduction
Individual snow crystals form in clouds when supercooled, supersaturated water vapor nucleates as ice onto a particulate. In this environment, crystals may develop a wide variety of shapes, but all with a basic hexagonal structure, as demonstrated by Reference Bentley and HumphreysBentley and Humphreys (1931). In a classic study, Reference Maeno and EbinumaNakaya (1954) determined a relationship between the dominant crystallographic growth direction and temperature. From 0˚C to ~–4˚C plates develop (fastest growth in the plane of the a axes); then down to ~–10˚C columns grow (c axis fastest); a transition back to plates is observed down to ~–22˚C, below which columnar crystals again occur. Secondary features (e.g. branching and dendrites on plates, or needles and hollow sheaths on columns) develop when the supersaturation level exceeds that of saturation over liquid water at the given temperature (Reference LibbrechtLock, 1990). Branching growth on isolated crystals grown from the vapor begins when the crystal growth is dominated by water-vapor particle diffusion. This type of growth is most notably observed at –5˚C where needle (c-axis) growth is dominant and at –15˚C where branching along the a axis is most apparent. Complexity increases with the level of supersaturation (Reference Korolev, Bailey, Hallett and IsaacLibbrecht, 2005).
When the crystals collect as a snow cover, dynamic morphologic alterations ensue. In a dry snow cover, depending largely on temperature gradient, one of two general processes dominates (de Quervain, 1958; Reference AkitayaAkitaya, 1974). Larger temperature gradients (>20˚Cm–1) lead to what Reference ColbeckColbeck (1982) called the kinetic growth form. The process is particularly prevalent near a relatively warm base, in which case the resulting structure is termed depth hoar (Reference Rosenthal, Saleta and DozierSeligman, 1936). Lower temperature gradients yield a structure composed of rounded, sintered grains, referred to by Reference ColbeckColbeck (1982) as the equilibrium growth form. These processes are also often referred to as ‘temperature gradient’ and ‘equitemperature’ metamorphism, respectively, although actual equal temperature conditions, as implied by the term equitemperature, are unlikely. Macroscopically, the sintered snow is, in general, structurally stronger.
Sintering is a mechanism that is concerned with the consolidation of powdered materials, resulting in densification, grain growth and microstructural changes. It is commonly exploited in engineering applications in ceramics and powder metallurgy. Although sintering dynamics are quite complex, insight into the process is obtained by considering two spheres in contact. In this two-particle model, the amount of material supplied to the neck (or bond) region is largely based on geometry. Forces that govern sintering in this model result from surface curvature variation, which influences differences in bulk pressure, vacancy concentration and vapor pressure. Reference IvantsovKang (2005) lists six material transport mechanisms: (1) lattice diffusion from the grain boundary to the neck, (2) grain boundary diffusion from the grain boundary to the neck, (3) viscous flow from the bulk grain, (4) surface diffusion from the particle surface to the neck, (5) lattice diffusion from the particle surface to the neck and (6) gas phase transport by either evaporation/ condensation or gas diffusion. Material supplied from the grain boundary will result in the grains moving closer together, whereas the other sources yield geometric changes but not contraction.
Sintering in snow has most frequently been attributed to vapor transport (Reference De QuervainHobbs and Mason, 1964; Reference LockMaeno and Ebinuma, 1983). However, this mechanism does not account for shrinkage or densification. Reference ColbeckColbeck (1998) presents a grain boundary diffusion model for snow, and there is some physical evidence that grain boundary diffusion may play a role in snow sintering (Reference Adams, Miller and BrownAdams and others, 2001). Although Reference Adams, Miller and BrownAdams and others’ (2001) analysis was for a pure ice system, they indicated that impurities would tend to concentrate at grain boundaries and influence the diffusion. Reference BarnesBarnes (2003) suggested that the grain boundary ridges observed by Reference Adams, Miller and BrownAdams and others (2001) may indeed be impurities and offered that, while not required, the presence of impurities is in accord with grain boundary diffusion and that soluble impurities could increase the diffusion rate. Reference NakayaRosenthal and others (2007) demonstrated evidence that the ridges are impurities and discuss the fact that dopants are used in enhancing engineered sintering processes.
In the present study, we examine microstructural metamorphism of disaggregated snow particles in the absence of an imposed temperature gradient. Microstructure and intergranular bonding between snow grains is of importance to avalanche stability, over-snow mobility, conductivity and other thermomechanical properties. A strong macrostructure is often observed to rapidly develop in snow that has been mechanically worked by mechanisms such as wind, avalanches, snowplows and snowmillers, etc. (Reference AbeleAbele, 1990). Causal processes that have been attributed to this include mechanical heat dissipation, which produces slight melt followed by refreeze between grains, and particle sintering mechanisms.
Methods
The impetus of this study was to observe sintering in processed snow. Older, rounded grained natural snow that had been stored in an insulated container in a –10˚C freezer for several months comprised the initial constituent. Sample preparation was carried out in a cold laboratory near the temperature of the intended metamorphism study. Snow was disaggregated by sifting through a 0.99 mm sieve into a 16mm diameter, 2 mm tall crucible. The vessel used in most of the studies was constructed by affixing a stainless-steel ring to a thin glass plate. Steel was chosen to assist in achieving a uniform sample temperature. Immediately after snow was sifted into the crucible, it was transferred onto a temperature-controlled silver block, which is the core of the enclosed microscope cold chamber. A glass cover slip was placed over the crucible to facilitate sample temperature uniformity and minimize vapor exchange. In this paper, we refer to these cases as ‘equitemperature’ experiments.
For the short transfer distance to the microscope, which was outside the cold laboratory, samples were placed in a small insulated container filled with snow and cold packs. It should be noted that the cold-stage chamber must be briefly opened to place the sample, causing a perturbation in the stage environment.
The Linkam BCS196 microscope stage uses liquid nitrogen injected directly into the silver block by a LPN94/2 cooling pump which is controlled by a TMS94 temperature programmer. Temperature feedback is provided by a 1/10 DIN class B Pt100 platinum resistor to give <0.01 s response time. A second cryostage was also used in two of the experiments discussed here. The Linkam GS120 is a temperature gradient stage that can be used to program gradients across a sample by using two aligned Peltier-based thermal elements. These elements both lie on the focal plane of the microscope, separated by a 2.5 mm gap. Each individual element can be independently controlled to 0.1˚C from –25 to 120˚C to establish a temperature gradient across the stage, transverse to the microscope observation orientation.
Metamorphism was observed over the course of 2–4 hours, and time-lapse images recorded using a Nikon Eclipse 80i microscope with a Qimaging Retiga-2000R digital camera. The charge-coupled device has 1600 by 1200 square (7.4 μm×7.4 μm) pixels. The viewing area was selected to examine grain-to-grain and grain-to-pore interaction. This utilized a 100× microscope magnification to yield a 1.4 μm resolution. An effort was made to focus between the upper grains into a region below the top surface to achieve conditions more representative of the interior of a natural snow cover.
In general, the samples were back-lit through a small hole in the silver block or between the Peltier elements. The light intensity was kept to a minimum, to reduce its potential thermal influence, and compensated with exposure time. For the same reason, the exterior of the cryostage was preconditioned by placing cold packs on top, prior to beginning an observation sequence.
As is typical of experimental protocols, the design described evolved through nuanced refinements as efforts were made to eliminate conditions that might unintentionally influence the outcome. Consequently, although all of the experiments did not strictly adhere to all aspects of the protocol presented, the resultant unexpected metamorphism (described below) was reasonably consistent. The question remained, however, as to whether this development was an artifact of the scale of the experiment and irrelevant to larger systems.
To examine a somewhat more in situ type scenario, without boundary influence, a separate procedure was conducted. Instead of being sifted directly into the crucible, snow was sifted through the 0.99 mm sieve to 70 mm depth in an 85 mm diameter cylindrical paper cup that remained undisturbed in the cold laboratory for 4 hours. A sample was then taken from the center of the specimen, placed in the crucible, inserted into the cryostage and immediately observed, presumably prior to additional metamorphism. The stage was set to –30˚C in an effort to slow further metamorphism in this ‘cup’ experiment.
Observations and Hypotheses
As stated, the notion of this study was to observe sintering processes of mechanically processed snow in a controlled thermal environment. Instead of examining two-particle sintering, processes relevant to an assemblage of particles were considered. Unexpectedly, protrusions developed on some surfaces of individual ice grains that were not in contact with other grains. Furthermore, the growth of these was seemingly random in direction on different grains (Fig. 1). Reference KangKorolev and others (2004), using a thermal diffusion chamber in a similar temperature range to that employed here, have observed the growth of extended features onto isolated, generally smooth spherical ice drops.
The protrusions that developed in the present study are of particular interest because when they were sufficiently close to a neighboring grain they were observed to connect with the adjacent grain. Bonds were thus established, although the grains were not initially in contact. Why did these relatively slender protrusions grow in a preferential direction on a preferred face on individual crystals, yet develop a similar morphology but in different directions on other crystals?
Ice crystals grown from the vapor onto an ice substrate in the presence of a temperature gradient will assume the crystallographic orientation of the substrate crystal (Reference Adams and MillerAdams and others, 2003), and crystal growth within the snow cover must adhere to the same physics as atmospherically grown crystals. We hypothesize that the selective orientation and habit of the protrusions is governed by the preferred crystallographic morphology dependence on temperature. The needle-like projections observed in Figure 1 developed with the cold stage at –5˚C. The temperature-dependence hypothesis is supported in observations of a separate experiment run at –15˚C (Fig. 2) in which the new growth appears to be more plate-like in structure.
There is a potential that these observations might be influenced as the result of a temperature gradient that could develop in the microscope chamber. There is possible variation in the chamber since temperature is controlled by the stage and there is a glass observation port above the sample. It was not possible to determine any preferred growth orientation using the microscope. Therefore, the gradient stage was employed to examine the progressive metamorphism that might be expected in the presence of a prescribed large temperature gradient. The Peltier elements were set to –5˚C and –15˚C, and the ensuing microstructural changes observed (Fig. 3). A temperature gradient was established in the direction of the crucible diameter. In this environment, faceted striated crystals aligned with the imposed temperature gradient developed at the expense of the adjacent grains. This demonstrates classic temperature gradient metamorphism, described as the ‘hand-to-hand’ (Reference SeligmanYosida and others, 1955) transfer of water vapor across the pore space from the warmer to colder grain.
It is important to note that based on observation in the equitemperature experiments the grain surface to which the new growth joined did not measurably sublimate, nor grow, although on the same grain there was growth at other sections in different directions (Fig. 1).
In the equitemperature studies, the filling-in of necks joining grains which were initially in contact appears to follow the two-particle sintering concept simultaneously with the projection growth (Fig. 4). The column appears in the lower right section of the bottom image. This particular equitemperature experiment was performed using the gradient stage with both elements set to –5˚C. This observed hollow columnar prism may have been influenced by the 0.1˚C variations possible in the Peltier elements. This type of hollow faceted development was not noted on samples in the more precise stage.
Examinations of samples taken from the cup experiments indicate that the protrusion growth will also develop in regions further beneath the surface than that which was observed to form in the microscope stage (Fig. 5). However, these protrusions are smaller than those that were observed growing under the cryostage environment for a similar duration. Furthermore, since the protrusions may develop quickly, this demonstration does not preclude the potential that these formed during the first few minutes in the cryostage.
There were instances where the needles grew for a time, but when they did not connect with another particle they were observed to retreat. The shallow depth of field of the microscope stage coupled with the metamorphism process over time caused difficulty in maintaining image focus and illumination consistency, making longer-term microscope observations problematic.
Discussion
Development of the crystal growth from the vapor within a snow cover is examined in this approximation, assuming that the essential physics that determine the growth rate and habit are the same as for atmospherically grown crystals. Although this assumption is utilized here, it should be noted that there are substantive differences between the growth from the vapor of isolated crystals and that which develops in an aggregate such as snow. For example, the concept of the far field and uniform mean vapor distribution would not be the same. In the following analysis, we consider aspects of the summary presented by Reference Korolev, Bailey, Hallett and IsaacLibbrecht (2005). Crystal growth is considered in terms of the simple hexagonal prism and the development of the dendritic form.
The velocity of the ice growth normal to the surface, based purely on attachment kinetics-limited growth, is given as (Hertz–Knudsen formula),
where α ≤ 1 is the condensation coefficient (1 if an atom hitting the surface immediately sticks) and σ surf is the supersaturation just above the ice surface.
is described by Reference Korolev, Bailey, Hallett and IsaacLibbrecht (2005) as ‘the rate at which the surface is simultaneously evaporating and regrowing’. ρ sat and ρ solid are the densities for water-vapor equilibrium over a flat ice surface and solid ice respectively, k is Boltzmann’s constant, T is temperature and m is the mass of a water molecule. Faceted crystal growth is taken to be attachment kinetics-limited. However, water-vapor particle and heat diffusion also play a role in crystal morphology. When the ice growth velocity is ≪αv kin, it is considered to be a diffusion-limited process. At –5˚C, v kin = 496 μms–1, which for the approximately constant protrusion growth velocity observed (v tip ≈ 0:024 μms–1) implies that diffusion-limited growth is likely.
A precise development addressing dendritic growth within the snowpack is not yet available; however, some useful insight is presented by utilizing a number of assumptions in the development that is implemented below. Early work in the area of dendritic growth was established by Reference Hobbs and MasonIvantsov (1947), and a significant body of work has followed.
In air, water-vapor particle diffusion is taken to restrict growth more than does heat diffusion. Reference Korolev, Bailey, Hallett and IsaacLibbrecht (2005) considers a particle diffusion model where the velocity and radius of a cylindrically symmetric parabolic dendrite tip, v tip and R tip, are taken to remain constant. Dendritic tip velocity is approximated as
D = 2×10–5m2 s–1 is the temperature-dependent diffusion constant for water vapor in air, α is the condensation coefficient, s is a dimensionless solvability parameter, σ ∞ is the supersaturation at a distance ‘far’ from the growing surface, and B is taken by Libbrecht to be equal to 10 for experimental studies. R kin is presented as
For molecularly rough surfaces, α ≈ 1, although accurate values are difficult to determine. s, which is taken as a constant, dependent on the intrinsic attachment kinetics, has not yet been calculated. However, in air at 1 atm, R kin = 30 nm. Therefore if R tip is known, αs can be estimated from Equation (4). Then through the measured value for the tip velocity, v tip, an approximation of the supersaturation (σ ∞) may be ascertained using Equation (3).
The dendritic growth shown in Figure 1 was examined to obtain a coarse approximation of the tip velocity and radius. These were acquired using a scale (integral to the microscope software) to visually determine and manually mark positions on the images with a cursor. Unfortunately, there is considerable subjectivity associated with the method. The ice-boundary edges in the image are not precise, and determination of what constitutes the actual tip radius is inexact. In addition, any slope of the protrusion into or out of the plane will influence the apparent velocity, because the image is based on a projected length. Keeping in mind these difficulties and the limitations of the model, the level of supersaturation in the pores was calculated. The assumption that the tip velocity and radius remain constant appears to be generally consistent with these measurements. Taking as an average v tip = 0:024 μms–1 and R tip = 9 μm produces a value of σ ∞ = 0:015, which yields an excess vapor density (supersaturation density) of ≈0.05 gm–3. This supersaturation and temperature places us in a region consistent with solid prism growth and below the vapor saturation over supercooled water at –5˚C (Reference Korolev, Bailey, Hallett and IsaacLibbrecht, 2005).
The significance of the observations noted in this study is the potential relation to the strengthening of processed snow. However, extrapolation of the observations of this study to the large-scale field environment must be scrutinized. Among the potential unintended influences innate to the experimental design are: the constraint to optically focus on subjects slightly below the surface, with the potential of significant surface boundary effects; the perturbation of the cryostage environment when the sample is inserted; heating from the microscope illumination source; and slight temperature gradients that may result from the viewing window above the sample and the cold platform. Measures taken to address these issues, such as minimizing the light intensity, using steel for the vertical cylindrical walls of the crucible, and placing cold packs on the external stage environment, did not appear to alter the observations. Morphologic changes that developed in the presence of an imposed temperature gradient were markedly different from those that developed under the equitemperature conditions. The presence of growth protrusions in the cup experiment, designed to mitigate the influence of the surface boundary, adds credence to the prospect of the protrusion development in a natural environment.
Development of the observed metamorphism that occurs in the experimental conditions and which is potentially relevant to the full-scale environment suggests that, depending upon conditions, a number of processes may play a role in snow sintering. We might speculate that while radius-of-curvature driven evaporation/condensation processes dominate under some conditions, the dendritic growth discussed here may play a role in some specific situations, and grain boundary diffusion may perhaps play a more obvious role in the longer term.
Conclusion
A conglomeration of disaggregated snow subjected to constant temperature conditions can develop protrusions on individual grains that grow into the interstitial pore space. Growth occurs on specific regions of the ice surface, determined by the crystallographic orientation of the nucleation grain. The preferred growth plane and the crystallographic habit are a consequence of the recognized relationship of temperature with the morphology of ice crystal growth from vapor. It is likely that the observed dendritic development is particle-diffusion limited, with respect to water vapor in air. We hypothesize that the excess vapor density, which drives this process, is due to the availability of the high-surface-energy contours that result from the sieving process employed in disaggregation of the old sintered snow. Fine ice particles and sharp corners that result from fracturing will sublimate, thus providing a mass supply for crystal growth. Connections between two grains can develop when the temperature, the crystallographic orientation of the nucleate grain and the proximity of a neighboring grain are conducive. Protrusions that have connected with another grain may continue to sinter, filling in and strengthening the neck region. But those that do not may disappear when the supply of mass from the fine particles and corners is exhausted and the protrusion tips themselves become regions of relatively high vapor pressure. The dendritic growth assumes the orientation of the crystal from which it is derived; consequently, we conjecture, as sintering progresses, the two originally separate units may assume the morphology of two crystals joined at a single bond interface.
It is well known that snow that has been mechanically worked, such as avalanche debris, wind deposits and snowplows or snowmillers, rapidly strengthens macroscopically (Reference AbeleAbele, 1990). Results of the present study suggest that one mechanism in the process may be the growth of dendritic protrusions that grow between and join ice grains. These produce bonds in addition to those that form between grains that are initially in contact and which develop through mechanisms described by the two-particle sintering model. Although this mechanism has intriguing implications, it must be noted that conditions in this study were small-scale and may have been influenced by the boundary and initial conditions of the experimental protocol. Additional work is warranted to validate the relevance of the observed results to field-scale studies.
Acknowledgements
We thank the M.J. Murdock Charitable Trust and the US National Science Foundation (NSF) Award ANT-0521360 for instrumentation to support this work, as well as NSF grant EAR-0635977.