Hostname: page-component-78c5997874-xbtfd Total loading time: 0 Render date: 2024-11-19T08:38:50.316Z Has data issue: false hasContentIssue false

Air and Water-Vapour Convection in Snow

Published online by Cambridge University Press:  20 January 2017

N. Klever*
Affiliation:
Federal Institute for Geosciences and Natural Resources, Alfred-Benk-Haus, Stilleweg 2, D-3000 Hannover 1, FDR
Rights & Permissions [Opens in a new window]

Abstract

Heat and mass transport by an air and water-vapour mixture (pore air) in snow due to thermal convection has been computed. It can be shown theoretically that thermal convection is always occuring in newly fallen snow, but never occurs in snow with rounded grains. Rayleigh numbers of depth hoar and partly settled snow are close to the critical Rayleigh number, therefore a small change in one or more of the snow parameters will decide whether there is convective or conductive heat transport.

The influence of thermal convection on the mass transport rates of water vapour in snow is computed based on the assumption of the ventilation approach. Although thermal convection occurs more often than previously assumed, its influence amounts only to about 2-3%.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1985

INTRODUCTION

Since the early papers of Paulcke (1934), Seligman (1936) and Bader (1939) several authors have tried to understand the role of the flux of water-vapour on the metamorphism in dry snow. The influence of thermal convection of the pore air on metamorphism has been a matter of question for a long time (de Quervain 1963). Yen (1962, 1963) and de Quervain (1972) have introduced a term for moisture exchange due to ventilation to the equation of mass balance. A detailed dimensional analysis of heat and mass transfer in a snowpack was published by Male and others (1972). They concluded that the development of general analytical, numerical or experimental models will prove impractical, and suggested that field studies of situations involving a restricted number of parameters is the only tractable approach. First measurements of air circulation in snow were made by Neher (1939). Trabant and Benson (1972, 1973) were among the first to attribute their high observed mass flux to thermal convection. Akitaya (1974) examined the possibility of thermal convection in snow by laboratory experiments. As far as I know only one significant quantitative study of theoretical models has been published until now (Palm and Tveitereid 1979). Both Akitaya and Palm and Tveitereid postulated that thermal convection only occurs in extreme situations like coarse grained old snow with strong vertical temperature gradients and large air permeabilities. A weakness of Palm and Tveitereid’s paper is that the models are based on permeability data of Kuroiwa (1968), which are only classified between fine and coarse grained snow, Further, these authors took the same thermal conductivity for all different kinds of snow.

In the last years the interest in problems connected with pore air in snow has increased, as shown by the papers of Powers and others (1985, this issue) and Klever (1983, 1984).

GOVERNING EQUATIONS AND NUMERICAL METHOD

Fluid flow in porous media is governed by the equations of conservation - energy, momentum and mass -in combination with the equation of state. Combining these equations and introducing the streamfunction Ψ by:

Lead to the following coupled equations

All variables in these equations have been made dimensionless by normalizing factors, indicated by primes (Table 1). This process defines the dimensionless Rayleigh number in the following way:

Table I. NORMALIZING FACTORS FOR THE SNOW MODELS. LENGTH AND TEMPERATURE ARE MEAN VALUES FOR SNOW LAYERS: PERMEABILITY AND THERMAL CONDUCTIVITY ARE MEAN VALUES DERIVED FROM THE EQUATIONS OF SHIMIZU AND DEVAUX. THE PROPERTIES OF THE FLUID ARE TAKEN FROM STANDARD AIR PROPERTIES

When simulating thermal convection of pore air in snow with the help of the above equations the influence of the following physical processes is neglected:

  • the volume change due to compressibility of a gas

  • the process of diffusion of water-vapour and

  • phase exchanges solid-fluid-gas.

Concerning the first simplification it can be argued that this is acceptable for velocities that are small compared with the speed of sound. The others imply that mass transport by metamorphism cannot be taken into consideration within the numerical calculations. Therefore no conclusions concerning net density changes are possible.

The numerical method used was a finite-difference approximation to the differential equations. To avoid numerical instabilities the Il’in scheme is used (Gorenflo and Kiesner 1982). This method is extended to simulate sloped regions with a constant or changing gradient. A detailed description of the theoretical and numerical procedure is given by Klever (1984).

RAYLEIGH NUMBER AND BOUNDARY CONDITIONS

An example for normalizing factors is given in Table I. From these factors a Rayleigh-number can be computed to Ra = 25.656. All models have the same Rayleigh number. The definite change between conductive and convective heat transport is described by the critical Rayleigh number. To compute models for different types of snow, permeability and thermal conductivity are introduced numerically by functions and not by Rayleigh number. This formulation leads to different critical Rayleigh numbers for the different types of snow. This procedure seems to be legitimized by its usefulness for comparing dimensionless results as well as for the introduction of lateral variation of permeability and conductivity (Klever 1984). The critical Rayleigh number Rayleigh in this paper is connected to the critical Rayleigh number Ranc normally used in literature by

Λt is the thermal conductivity of the whole system defined by

where λs, and λf refers to the thermal conductivity of snow and fluid (air).

Ranc = 27.1 is the critical Rayleigh number of the models with the following boundary conditions (Ribando and Torrance 1976): temperature is kept constant at the top and bottom and adiabatic conditions exist on both sides. The boundary conditions for the streamfunction are a permeable top and impereable boundaries elsewhere. This should simulate a normal snowcover under normal meteorological conditions.

For our models λ s and k were computed by

where d means grain size and λ f was assumed to be 0.0026 W/m°C (standard air).

Fig. 1a. Normalized stream function Ψ/ Ψ *[1] (where Ψ* = 0.139 g/m s) of model TSM-N-02 with newly fallen snow. The Rayleigh number of this model is Ra = 25.656 and the critical Rayleigh number is here Rac = 6.295 (see Table 2).

Fig. 1b. Normalized temperature T/T* (1] (where T* = 10°C) of model TSM-N-02.

TABLE 2.

MODELS

Mean values for density, porosity and crystal size (Table 2) have been taken from statistics of seasonal snowcover at Weißfluhjoch/Davos (de Quervain 1980/81). Permeability, thermal conductivity and the critical Rayleigh number (Table 2) have been computed from these values according to the equations given above. The presence of deviations from these mean values was considered by changing the normalizing factors and therefore the Rayleigh number Ra.

It is remarkable that in snow with rounded grains no sensible variation of the parameters can change Ra to exceed the critical Rayleigh number Rac. Therefore no convection will occur. Just the opposite is true in the case of newly fallen snow: here the Rayleigh number Ra will be greater than the critical Rayleigh number Rac for all variations of the parameters in the admissible range. Therefore convection always seems to be possible in newly fallen snow. However, in the case of partly settled snow and depth hoar slight changes in some parameters can have a great bearing on whether convection occurs or not.

In Figure I the normalized streamfunction and temperature are shown for a model with newly fallen snow. Using the factors from Table 2 the velocity of pore air in this model can be computed to 0.31 mm/s. The maximum velocity is computed from a depth hoar model (not shown here) to 4.2 mm/s. Assuming that the ventilation approach (Colbeck 1983) can give a first approximation to the mass transport by vapour, this velocity has only a small influence: from the ventilation factor an increase of the mass transport rate is computed only to 2.8%.

Fig. 2. Different types of convective cells in a sloped layer as functions of the experimental conditions (Ψ, Ra): A) unicellular flow, B) polyhedral cells, C) longitudinal stable coils, D) fluctuating regime, and E) oscillating longitudinal coils (from Combarnous and Bories 1975).

From theory and experiments it is known that convection always occurs in sloped regions (Combarnous and Bories 1975). These authors have shown that for convection in sloped layers, different types of convective cells are occurring (Figure 2). Because of the large extent of the cells in unicellular flow region, the experiments of Neher (1939) can therefore approximately be re-interpreted as a hint at convection in sloped regions.

DISCUSSION

It can be seen from Table 2 that thermal convection in snow occurs more often as previously assumed. The application of the ventilation approach as a first approximation to the mass transport shows that the computed velocities have only a mean influence of about 2-3% on the mass transport by vapour. On the other hand, based on a paper by Keller and Hallett (1982), Colbeck (1983) suggests an effect of thermal convection on crystal shape. Starting from this idea and comparing the schematic picture of Wakahama (1968) (Figure 3) with a convective cell, I believe that thermal convection may be the initial state of the occurrence of the drainage of melt water or water channels in the early melt season. There are some objections to this suggestion: (1) the invariance of the geometry of the percolation channels as compared to the possible variability of the geometry of the convection cells due to different snow layer depths and due to transient state of convection, (2) the occurrence of the channels in snow layers, where convection seems very unlikely, and (3) more plausible mechanisms such as radiation and melting effects. However the obvious similarity remains, especially in the case of a sloped region if the longitudinal stable coils of the cell-type classification of Combarnous and Bories (1975) are taken into consideration. Much more work has to be done to prove or disprove this suggestion.

Fig. 3. Schematic representation of the layered structure in a snow cover in the early melt season. Large “drains of melt water or water channel” are often found in the snow cover as shown in this figure (from Wakahama 1968)

References

Akitaya, E 1974 Studies on depth hoar. In International Association of Hydrologicat Sciences Publication 114, Symposium on Snow Mechanics: 42-47 Google Scholar
Bader, H 1939 Mineralogische und strukturelle Charakterisierung des Schnees und seiner Metamophose, In Bader, H and others (eds) Der Schnee und seine Metamorphose, Beiträge zur Geologie der Schweiz Geotechnische Serie Hydrologie, Lieferung 3. Bern, Kümmerly und Frey: 3-61 Google Scholar
Bader, H, Haefeli, R, Bucher, E, Neher, J, Eckel, O, Thomas, Chr 1939 Der Schnee und seine Metamorphose. Beitrage zur Geologie der Schweiz Hydrologie, Lieferung 3. Bern, Kümmerly und Frey Google Scholar
Benson, C S, Trabant, D C 1973 Field measurements on the flux of water vapour through dry snow. In Proceedings of the International Symposium on The Role of Snow and Ice in Hydrology, Banff. 6-13.9.1972. Banff, World Meteorological Organization (WMO), International Association of Hydrological Sciences HAHSl’ UNESCO: 291-298 Google Scholar
Colbeck, S C 1983 Theory of metamorphism of dry snow. Journal of Geophysical Research 88, C9: 5475-5482 Google Scholar
Combarnous, M, Bories, S A 1975 Hydrothermal convection in saturated porous media. Advances in Hydro-science 10: 231-307 Google Scholar
Devaux, J 1933 L’economie radio-thermique des Champs de neige et des glaciers. Annales de Pkysicae 10, 20: 5-67 Google Scholar
Gorenflo, R, Kiesner, St 1982 Über das ?Γ?η -Differenzenverfahren für die Diffusions- Konvektions-Gleichung, International Series for Numerical Mathematics (ISNM) 58: 73-89 Google Scholar
Keller, V W, Hatlett, J 1982 Influence of air velocity on the habit of ice crystal growth from the vapor. Journal of Crystal Growth 3 60: 91-106 CrossRefGoogle Scholar
Klever, N 1983 Konvektion von Wasserdampf im Schnee? Paper, held at 13. Int. Polartagung, Bamberg, 4-7 Oktober 1983, 15 pGoogle Scholar
Klever, N 1984 Stationäre Konvektion in porösen Medien numerische Untersuchungen an unterschiedlichen Fragestellungen aus der Hydrothermik und der Schneemetamorphose. Berliner Geowissen-schaflliche Abhandlungen, Reihe, B, Heft 11, 114 pGoogle Scholar
Kuroiwa, D 1968 Liquid permeability of snow. Low Temperature Science A, 26: 29-52 Google Scholar
Male, D H, Norum, D I, Besant, R W 1972 A dimension analysis of heat and mass transfer in a snowpack. Proceedings of the International Symposium on The Role of Snow and Ice in Hydrology, Banff 6-13,9.1972. Banff. World Meteorological Organization (WMO), International Asociation of Hydrological Sciences (IAHS), UNESCO.- 258-290 Google Scholar
Neher, J 1939 Schneeuntersuchungen im Gelände. In Bader and others (eds) Der Schnee und seine Metamorphose. Beiträge zur Geologie der Schweiz Geotechnische Serie Hydrologie, Lieferung 3. Bern, Kiimmerly und Frey: 255-272 Google Scholar
Palm, E, Tveitereid, M 1979 On heat and mass flux through dry snow. Journal of Geophysical Research 84: 745-749 Google Scholar
Paulcke, W 1934 Eisbildungen 1. Der Schnee und seine Diagenese. Zeitschrift fur Gletscherkunde XXI(4/5): 259-282 Google Scholar
Powers, D J, Colbeck, S C, O’Neill, K 1984 The role of thermal convection in snow metamorphism. Annals of Glaciology 6: 43-47 Google Scholar
de Quervain, M 1963 On the metamorphism of snow. In Kingery, W D (ed) Ice and Snow. Cambridge, Mass, MIT Press: 377-390 Google Scholar
de Quervain, M 1972 Snow structure, heat and mass flux through snow. Proceedings of the International Symposium on The Rote of Snow and Ice in Hydrology, Banff, 6-13.9.1972. Banff, World Meteorological Association (WMO), International Association of Hydrological Sciences (IAHS), UNESCO: 203-225 Google Scholar
de Quervain, M 1980/81 Schneekunde, Lawinenkunde, Lawinenschutz. Skript einer Einführungsvorlesung, gehalten an der ETH Zürich, im WS 80/81 (13. Ausgabe) Seligman G 1936 Snow structure and ski fields. London, MacMillan Google Scholar
Shimizu, H 1970 Air permeability of deposited snow. Contributions from the Institute of Low Temperature Science, A, 22: 1-32 Google Scholar
Trabant, D, Benson, C 1972 Field experiments on the development of depth hoar. In Doe, B R, Smith, D K, Gruner, J W (eds) Studies in Mineralogy and Precumbrian Geology, Geological Society of America Memoir 135: 309-322 Google Scholar
Wakahama, G 1968 The metamorphism of wet snow. General Assembly of Bern, 25.9-7.10.1967. Commission of Snow and Ice. Reports and Discussions. Word, W (ed), Gentbrugge. IAHS Pubi. 79: 370-379 Google Scholar
Yen, Y C 1962 Effective thermal conductivitv of ventilated snow. Journal of Geophysical Research 67(3): 1091-1098 Google Scholar
Yen, Y C 1963 Heat transfer by vapour transfer in ventilated snow. Journal of Geophysical Research 68(4): 1093-1101 Google Scholar
Figure 0

Table I. NORMALIZING FACTORS FOR THE SNOW MODELS. LENGTH AND TEMPERATURE ARE MEAN VALUES FOR SNOW LAYERS: PERMEABILITY AND THERMAL CONDUCTIVITY ARE MEAN VALUES DERIVED FROM THE EQUATIONS OF SHIMIZU AND DEVAUX. THE PROPERTIES OF THE FLUID ARE TAKEN FROM STANDARD AIR PROPERTIES

Figure 1

Fig. 1a. Normalized stream function Ψ/ Ψ *[1] (where Ψ* = 0.139 g/m s) of model TSM-N-02 with newly fallen snow. The Rayleigh number of this model is Ra = 25.656 and the critical Rayleigh number is here Rac = 6.295 (see Table 2).

Figure 2

Fig. 1b. Normalized temperature T/T* (1] (where T* = 10°C) of model TSM-N-02.

Figure 3

Fig. 2. Different types of convective cells in a sloped layer as functions of the experimental conditions (Ψ, Ra): A) unicellular flow, B) polyhedral cells, C) longitudinal stable coils, D) fluctuating regime, and E) oscillating longitudinal coils (from Combarnous and Bories 1975).

Figure 4

Fig. 3. Schematic representation of the layered structure in a snow cover in the early melt season. Large “drains of melt water or water channel” are often found in the snow cover as shown in this figure (from Wakahama 1968)