1. Introduction

Mass transport, heat transfer and phase changes are the fundamental thermodynamic processes in the snow cover. Anim-portant phenomenon, metamorphism — the change of ice grain forms of the snow cover over time — is induced by these processes. Traditionally, metamorphism has been studied using the concepts of mass and energy conservation which require a detailed knowledge of the movement of water molecules on the ice surface or the diffusion of water vapor in the pore space. For an overview, see Reference ColbeckColbeck (1987). Rate equations describing grain growth have been proposed (Reference ColbeckColbeck, 1983; Reference GublerGubler, 1985; Reference Lehning, Bartelt, Brown, Fierz and SatyawaliLehning and others, 2002) with the goal of modeling the development of the seasonal snow cover (Reference Brun, Sudul, David and BrunotBrun and others, 1992; Reference Bartelt and LehningBartelt and Lehning, 2002).

The application of heat and mass transport principles is essential for understanding the processes involved in snow metamorphism; however, conservation principles do not include geometric aspects and cannot solve the problem of form. This is evident in the present generation of metamorphism models which predict growth rates and crystal size, but the form changes are still based on empirical par- ameterizations of field or laboratory observations (Reference Brun, Sudul, David and BrunotBrun and others, 1992; Reference Lehning, Bartelt, Brown, Fierz and SatyawaliLehning and others, 2002). Metamorphism can be described by such relationships, but not the underlying principle.

The question is whether there exists a simple thermodynamic function related directly to form. Such a function would provide a thermodynamic principle that would allow the determination of the natural snow structure under well- defined thermodynamic constraints such as the temperature state, the applied temperature gradient and mass supply. Heat-transfer engineers use such functions because they try to reduce the increase in entropy per unit time, the so-called entropy production (Reference BejanBejan, 1996). During the design of a thermodynamic system (e.g. a heat exchanger), their objective is to find an optimal geometric arrangement of system components such that entropy production is minimized. The minimum of entropy production corresponds to a thermodynamic state where the least work is dissipated (Reference PrigoginePrigogine, 1980, p. 88). In geophysics and biology the principle of minimum entropy production has been applied to see if the minimum entropy production yields the optimal or natural visible form of river basins or living trees (Reference BejanBejan, 2000). Natural scientists argue that the optimal form arises out of the competition between at least two mechanisms, or two driving forces, and their respective fluxes. A similar principle was used by Reference BejanBejan (2000) to predict the formation of snowflakes. Since the conjugate forces and fluxes are also the source of increasing entropy, the final form is the one that makes entropy production stationary. In a thermodynamic system, the driving forces can be temperature and pressure, and the corresponding fluxes are heat and mass flows.

In snow the mass fluxes are the reason for the changes in form of the ice grains. The driving force is the difference in water-vapor pressure over the ice surface which depends on temperature and the curvature radius of the ice grains. This force and flux produce entropy, together with other form- dependent and form-independent processes. In the following, we postulate that the form of the ice grains will be changed until the entropy production is a minimum within the constraints given by the thermodynamic variables (p, T, V). This will be a non-equilibrium steady state. In view of this, we must:

The irreversible processes are formulated using constitutive relations for heat conductivity and mass transfer. It is important to note that not all thermodynamic processes produce entropy. For instance, sublimation and subsequent deposition do not produce entropy because phase changes are reversible and thus isentropic processes. The model for the geometric arrangement of the ice grains will be a commonly used model for snow microstructure applications (Reference Lehning, Bartelt, Brown, Fierz and SatyawaliLehning and others, 2002; Reference MillerMiller, 2002): the grains are uniform spheres of polycrystalline ice that are bonded together by necks (Fig. 1). Such forms are found typically in dry, early-winter seasonal alpine snow packs. The single ice crystals are joined together at crystal boundaries; the grains are joined together at the necks. A grain can consist of many single crystals. The crystal boundaries can be seen as fine dark lines in Figure 1. Under sufficient magnification, they appear as grooves. The migration of molecules along these grooves is not the subject of our investigations. Hence, for a given density, our model contains two independent form parameters, the sphere’s curvature radius and the radius of the neck. We will choose the curvature radius of the grain surface as the appropriate form variable. We can easily find the minimum entropy production by differentiating the sum of the entropy-producing terms with respect to curvature radius.

Fig. 1. A detail of snow showing the grain bonding and crystal boundaries (SLF). Note the existence of both positive and negative curvatures.

In summary, we hypothesize that metamorphism develops toward a non-equilibrium steady state where the entropy production is a minimum and therefore stationary. This procedure provides unambiguous information about form changes under the given thermodynamic constraints. These are interpreted within the framework of isothermal- and temperature-gradient snow metamorphism.

2. Entropy Production in Snow

We consider snow to be a porous medium consisting of ice (subscript i) and air (superscript a). The ice phase consists of grains which are bonded together to form a complicated and irregular but load-bearing ice skeleton. The ice skeleton is considered to be rigid. The pore air contains water vapor (supscript v) which diffuses through the pore space of the granular skeleton. The diffusion rate is given by the mass flow rate of the air ṁ _a; the density of the dry air is ρ_a; the density of the water vapor ρ _v. The partial vapor pressure is p _v and is related to the density via the ideal gas law. The water vapor has the same temperature as the dry air, T _v = T _a. We assume no free (fluid) water is present on the ice, i.e. the snow is dry. The temperature of the ice, T _i, and that of the adjacent air, T _a, are not in thermal equilibrium, T _i ≠ T _a. In the presence of a temperature gradient, heat is being transferred in the ice (heat flux Q _i) by conduction, and air (heat flux Q _a) by conduction and advection. Temperature gradients are induced by meteorological forcing at the snow-cover surface; the heat flux and vapor fluxes are primarily in one direction, that is, in the direction of the temperature gradient (the z direction). A list of the thermodynamic variables is provided in Table 1.

Table 1. Nomenclature of thermodynamic variables. The subscripts are: i = ice, a = air, v = vapor, g = grain, n = neck, b = boundary (interfacial) layer, s = saturation. The double prime superscript denotes a quantity per unit area

At the interfacial ice-air boundary (subscript b), heat (free convection) and mass are exchanged between the ice lattice and pore air. The heat flux across the boundary per unit surface area we denote the mass flux across the boundary per unit surface area we denote (The double prime denotes a quantity per unit area; a triple prime a quantity per unit volume.) The interfacial boundary is saturated with water vapor (saturation pressure p _s, concentration density ρ _s) which is maintained by sublimation or deposition processes from or to the ice surface. The heat flux across the interfacial boundary is driven by the temperature difference between the ice and air, T _i – T _a; the mass flux by the density difference, ρ _s-ρ _v. The layer has an unknown width. A schematic representation of this system is shown in Figure 2.

Fig. 2. Definition of thermodynamic and microstructural variables. Snow is defined as a porous medium of volumetric ice content θ_i and volumetric air content θ_a. The grains have positive radius r _g; the bonds have negative radius r _c. The pore space is filled with dry air (subscript a) and water vapor (subscript v). The temperature of air and ice can differ. Surrounding the ice grains is a saturated boundary layer of pressure p _s. The unit microstructural cell defines the control volume boundaries of our open thermodynamic system.

In summary, we consider a volume of snow to be an open thermodynamic system undergoing both heat and mass transfer. Due to symmetry, we chose the control volume such that it contains half an ice grain and half a bond (see Fig. 2). The control volume contains both ice and air. The relationship between grain-size and pore diameter is found by imposing the volumetric constraint θ _i + θ _a = 1, where θ is the volumetric fraction of a component.

For this open system, the second law of thermodynamics can be stated as (Reference BejanBejan, 1997):

(1)

where Ṡ is the total entropy production, and s the entropy density (J K^-1 kg^-1). The first term in square brackets is the rate of entropy accumulation inside the control volume due to interfacial boundary heat and mass transfer, Ṡ _b; the second term in brackets is the entropy transfer due to heat conduction in the ice and air; the third term in brackets is the net entropy flow rate out of the control volume via mass flow in the pore space. The entropy variables are listed in Table 2.

Table 2. Entropy variables. The subscripts are: g = grain, n = neck, b = interfacial mass transfer, h = interfacial heat transfer, k = conductivity, f = airflow. A triple prime denotes a quantity per unit volume

The entropy production per unit volume due to the heat flux across the interfacial boundary is given by (Reference BejanBejan, 1996):

(2)

where A _s is the specific surface area of the boundary (m² m^-3), in this case the total surface area of the interface within the control volume. The entropy generation is positive so long as the temperature difference exists. (If T _a > T _i then the sign of the heat flux is negative.)

The entropy production per unit volume due to the mass flux across the interfacial layer over a grain and necks of the control volume is given by (Reference BejanBejan, 1996)

(3)

and

(4)

respectively. Both expressions are required since the curvature of the grains and necks differ, meaning that both the mass flux and pressure gradient across the interfacial layer differ and must be considered separately. In summary, the total entropy accumulation within the control volume (the first part of Equation (1) per unit volume) is

(5)

The entropy generation due to heat conduction across the control-volume boundaries in the ice and air is (Reference BejanBejan, 1996):

(6)

and

(7)

where k _i and k _a are the thermal conductivities of ice and air, respectively. The entropy generation in the pore space per unit volume, due to the net entropy flow rate out of the control volume via mass flow, is (Reference BejanBejan, 1996)

(8)

where μ is the viscosity of air and Φ is the viscous dissipation function. This defines the third term in the statement of the second law, Equation (1), per unit volume.

The specific total entropy production Ṡ′″ is the sum of the entropy production in the ice, interfacial boundary layer and pore air,

(9)

In the above derivation, several assumptions were made that were not explicitly mentioned. Is entropy production additive? In principle, entropy production is not additive because the processes are not independent, as given by Onsager’s reciprocity relations (Reference BejanBejan, 1997). By adding the entropy production terms, we are implicitly assuming, for instance, that Soret and Dufour diffusion effects can be neglected. This implies that the cross-diagonal terms of Onsager’s relations are zero. Another assumption is that surface tension effects are reversible and therefore do not contribute to entropy production. Two other entropy producing processes we did not include in our model are radiation and viscous (irreversible) deformation. We assume that the ice matrix is rigid.

3. The Microstructural Model

Consider Figure 1 which shows a thin cross-section of a snow sample. The sample consists of well-rounded grains which are bonded together at necks. Note that both positive (concave) and negative (convex) curvatures exist in the sample. Furthermore, the curvature radius of the grains is in general larger than the grain radius. In the following, we do not attempt to model the complex geometric arrangement of grains. Rather we concentrate on the basic microstructural elements: grains and necks. Thus, the microstructural model consists of two spherical grains of radius r _g (Fig. 2). The distance between grain centers is 2r _g. A concaveshaped neck is constructed by revolving the bond arc length of radius r _c around the axis of grain contact, the z axis. Due to symmetry, only the section consisting of a half-neck (of length l _n) and half-grain (of radius r _g) is considered in the entropy minimization analysis. The important feature of this model is that the geometry is defined by only two parameters, the grain radius r _g and the neck radius r _c. The model has been employed by other authors to describe snow metamorphism (Reference Lehning, Bartelt, Brown, Fierz and SatyawaliLehning and others, 2002). Other snow types (e.g. dendritic new snow) would require a different microstructural model. A complete list of the microstructural variables is contained in Table 3.

Table 3. Nomenclature of microstructural variables. The subscripts are: i = ice, a = air, g = grain, n = neck. The variables are depicted in Figures 2 and 3

The relationship between the grain radius r _g and neck radius r _c can be expressed in terms of the angle α (see Fig. 3)

(10)

Fig. 3. Definition of microstructural variables. Spherical ice grains of radius r_g are connected by bonds of radius r _b. The radius of curvature of the bonds is r _c. The model consists of only two free parameters, r _g and r _c. α is half the bond angle.

The neck radius r _b is

(11)

The cross-sectional neck radius r _b should not be confused with the neck curvature radius r _c which is used to define the vapor pressure near the neck surface. The cross-sectional radius of the neck at the point of grain contact is denoted r _n and can be found by solving the following quadratic equation:

(12)

The length of the neck is

(13)

The surface area of the ice, A _s, is given by the sum of the neck surface area

(14)

and the grain surface area

(15)

Likewise, the volume of the ice is given by the sum of the neck volume

(16)

and the grain volume

(17)

The total ice volume is

(18)

The total volume, V _t, is the sum of the ice and air volumes

(19)

where θ _i and θ _a are the volumetric fractions of ice and air, respectively.

4. Entropy Production Components

4.1. Interfacial heat and mass transfer

The vapor saturation pressure in the interfacial layer over the ice grain is denoted p _sg and given by the vapor equation of Laplace

(20)

where is saturation pressure over a flat surface and σ is the surface tension. The interfacial layer over the ice is assumed to be saturated at temperature T _i, and thus the water-vapor pressure over a flat ice surface, , is found from the well- known Clausius-Clapeyron relation. Similarly, the vapor saturation pressure over the neck is denoted p _sn and approximately given by

(21)

Note that the bond, unlike the grain, possesses both convex (radius r _n) and concave (radius r _c) surfaces. The air is assumed to be saturated at vapor temperature T _v, and thus the water-vapor pressure in the air, p _v, can likewise be determined from the Clausius-Clapeyron relation.

The mass exchange from the ice to the surrounding pore space is given by the difference in vapor concentration between the boundary layer and air (Reference Incropera and DeWittIncropera and Dewitt, 2002). Specifically, the mass-transfer rate across the interfacial grain layer is

(22)

where h _m is the interfacial mass-transfer coefficient, ρ _sg is the vapor density of the interfacial layer and ρ _v is the vapor density of the pore space. Assuming an ideal gas, the mass- transfer rate can be expressed in terms of the pressure difference (Reference Incropera and DeWittIncropera and DeWitt, 2002)

(23)

where the interfacial mass-transport coefficient is re lated to h _m according to

(24)

M _v is the molecular mass of water vapor (M _v = 1.80 kg mol^-1) and R is the molar gas constant. The pressure difference Δp is defined as

(25)

that is, the part of the pressure difference between the interfacial layer and air which is independent of curvature effects. In a similar manner, the mass transfer over the neck surface can be found

(26)

Note that are positive when mass sublimates from the ice. In order to find the total entropy production in the control volume, Equations (3) and (4) must be integrated over the volume. We find the total entropy production rates for the grains and necks:

(27)

(28)

The interfacial heat flux between the ice lattice and pore space can be given by

(29)

where h is the local convection coefficient (Reference Incropera and DeWittIncropera and DeWitt, 2002). The total entropy production per unit cell for this process can likewise be found by integrating Equation (2) over the total boundary area:

(30)

Reference KavianyKaviany (1995) provides a relationship between the dimensionless Nusselt number (Nu = (2hr _g/k _a)) and air speed based on packed-bed experiments with spherical particles. For low Reynolds number, Nu ≈ 2. Thus, the local heat convection coefficient can be approximated,

(31)

This formula is based on an evaluation of both steady-state and transient packed-bed experiments, for a wide range of Reynolds number. It should be noted that few if any data exist for the interfacial heat-transfer coefficient for snow and that at low Reynolds number the above equation has not been experimentally verified (Reference KavianyKaviany, 1995).

The relationship between the convection heat-transfer coefficient and the mass-transfer coefficient can be found in Reference Incropera and DeWittIncropera and DeWitt (2002). We make the following approximation:

(32)

where p ≈ 1 kg m^-3 is the density of air, c is the specific heat capacity of air c ≈ 1000 J kg^-1 K^-1, and Le^1-n ≈ 1 is the dimensionless Lewis number, the ratio of thermal to mass diffusivities (see Incropera and DeWitt, 2002).

4.2. Heat conduction

The entropy produced in the control volume by heat conduction in the ice and pore air is (Reference BejanBejan, 1996)

(33)

and

(34)

We are assuming that the conductivities, k _i and k _a, are independent of z.

4.3. Air current in pore space

Air movement in snow is in the laminar, creeping flow regime. Reynolds numbers are small, Re <10. For this flow regime, various porous-media airflow models can be applied (Reference KavianyKaviany, 1995). Reference Buser and GoodBuser and Good (1987) found that a capillary tube model fits the experimental data of Shimuzu (1970) as well as their own data. The entropy production in tubes in the laminar flow regime has been determined by Reference BejanBejan (1996). For Hagen-Poiseuille flow the velocity distribution in the pore space varies parabolically (Fig. 2), according to,

(35)

where u _avg is the average air velocity and d the tube diameter. In a Hagen-Poiseuille flow model, the flow velocity tangential to the ice surface (r = d/2) is zero and reaches a maximum at the tube middle, r = 0. The entropy production per unit volume at any position r for this case is (Reference BejanBejan, 1996)

(36)

where μ is the viscosity of air (μ = 2.2 × 10^-5 kgm^-1s^-1). For our purposes, the entropy production in the control volume must be found by integrating the above equation over the length (r _g) and tube area. We find:

(37)

Typical average velocities are u _avg ≈ 1.0 × 10 ^-6 ms^-1 arising from pressure differences induced by temperature gradients. Higher values due to wind pumping at the surface of the snow cover can be found, but it is estimated that these values are seldom over u _avg < 1.0 × 10^-4 m s^-1.

5. Entropy Production Minimization

The total entropy production in the control volume is given by the sum of the entropy production in the ice, due to heat conduction (Equation (33)),

(38)

the interfacial layer, due to heat and mass transfer (from Equations (27), (28) and (30)),

(39)

and, finally, the entropy production in the adjacent pore space, due to heat conduction and fluid friction (Equations (34) and (37)):

(40)

Given the air and ice temperature, the temperature gradient and the pore air velocity, the total entropy production

(41)

can be expressed in terms of the grain radius, r _g, assuming a bond angle α.

Figure 4 plots the total entropy production as a function of the grain curvature radius r _g for two bond angles (a) α = 5° and (b) α = 20°. The plot depicts the entropy production at four different temperatures varying between T = 253 K and T = 272 K. The system is placed in thermal equilibrium, i.e. T = T _a = T _i; no temperature difference exists between the ice and air. Thus, the figure shows the entropy production minima as a function of varying isothermal temperature. At lower temperatures, the mimima are located at larger grain radii. In addition, the minima at lower temperatures are associated with higher entropy production levels. Note that when the system reaches an isothermal state near T = 272 K, a small change in equilibrium temperature influences the entropy production significantly. The change in entropy production from T = 253 K to T = 263 K is small in comparison to the change from T = 263 K to T = 270 K. At T = 272 K the entropy production minimum is at a very large curvature radius, not seen on the graph.

Fig. 4. Entropy production as a function of grain curvature for different temperatures. (a) α = 5º and (b) α = 20º. Production minima are located at large curvature radii with increasing temperature. At T = 272 K the entropy production minimum is at a very large curvature radius, not seen on the graph.

Figure 5 displays the contributions to the total entropy production, classified according to the different thermo-dynamical processes: interfacial mass transport, interfacial heat transport, heat conduction in the ice and air and fluid friction. The graph shows that the primary producer of entropy is clearly interfacial mass transport (by a factor of 100 over convective heat transport). Entropy production due to fluid friction in the pore air is insignificant. Since snow metamorphism is defined as the process of mass change under thermodynamic actions, this result underscores the connection between entropy production and snow metamorphism.

Fig. 5. Components of entropy production. Ice temperature T_i = 265 K; air temperature T _a = 263 K; air velocity u _a= 1.0 ×10^-6; temperature gradient dT/dz = 10Km^-1. The plot depicts the importance of interfacial mass transfer in the total entropy production.

6. Discussion of Results

6.1. Minimum entropy production at equilibrium temperatures (isothermal metamorphism)

The isothermal temperature results depicted in Figure 4 can be visualized from a different perspective. As shown, the entropy production curve as a function of r _g contains a well- defined minimum. In order to find the grain radius at which the total entropy production is minimal, we solve

(42)

The curvature radius at which entropy production is minimum for a defined thermodynamic state is termed the optimal curvature radius, r _opt. Figure 6 plots the optimal curvature radius as a function of the difference in temperature between the ice and air, ΔT = T _i – T _a, for three different ice temperatures. The three curves intersect the ΔT = 0 K axis at three different locations. The colder the temperature, the smaller the curvature radius. This agrees with observations made by Reference MarboutyMarbouty (1980). Note also that the optimal curvature radius for the isothermal case is not the largest radius possible. Larger optimal curvature radii, for a given ice and air temperature, are predicted when the temperature is higher in the ice than in the air.

Fig. 6. Optimal grain curvature radius as a function of temperature difference. The maximum optimal radii do not occur on the ΔT = 0ºC isothermal axis. The lower the temperature the farther away the peaks are from the isothermal axis.

The influence of bond size on entropy production is depicted in Figure 7a, which shows the entropy production for two different bond angles α = 5° and α = 45°. Isothermal conditions were imposed on the ice and air, T _i = T _a = 268 K. The bond angle a influences the location of the minimum entropy production. Smaller bond angles induce smaller optimal curvature radii. Evidently, the geometry of the microstructural model plays an important role in the application of the principle of minimum entropy. Figure 7b shows the influence of super saturation (1.05 % of saturation) and partial saturation (0.95% of saturation) on the entropy production. The graph shows that ropt is not sensitive to the state of saturation; however, the magnitude of the entropy production changes significantly.

Fig. 7. (a) Entropy production for different bond configurations, α = 5° and α = 45°. (b) Entropy production for different degrees of saturation, α = 5°. Calculations assume isothermal conditions: T _i = T _a = 268 K.

6.2. Minimum entropy production at non-equilibrium temperatures (temperature-gradient metamorphism)

Several other important results can be ascertained from Figure 6. First, note that the optimal radius curves are not symmetric around the isothermal ΔT = 0 K axis. The maximum radii for a given temperature are located to the right of the ΔT = 0 K axis; that is, for ΔT > 0K, or when the temperature in the ice is higher than the air temperature, T _i > T _a. The peaks are extremely sharp, meaning that slight changes in the thermodynamic state will produce significant changes in optimal curvature radius. Note also that the colder the ice temperature, the farther away the peaks are located from the ΔT = 0 K axis. For example, for an ice temperature of T _i = 258 K, the maximum optimal radius is located at ΔT = 2 K; for a higher ice temperature T _i = 272 K, the peak is located very near the ΔT = 0K axis. This result suggests that in order to obtain flat grains at low temperatures, a large temperature difference between the ice and air (T _i – T _a) must be maintained. This will seldom occur since the energy transfer (convective heat exchanges) between the ice and air phases will ensure a state near ΔT = 0 K. Thus, the theory of minimum entropy production predicts that very flat grains at colder temperatures are possible, but probably rare. Alternatively, at temperatures near T _i = 273 K, smaller temperature differences of the order ≈0.1 K are required to obtain the optimal curvature radius.

One way to induce and maintain temperature differences between the ice and air is to apply a temperature gradient dT/dz to the system, a snow layer or the entire snowpack. This can be shown by considering the following very simple model in which the ice and air components of a control volume are divided into two parts as shown in Figure 8. The goal of the analysis is to find the steady-state temperature gradient is applied. We are interested in the magnitude of the ice-air temperature difference ΔT(z) = T _i(z) — T_a(z). For simplicity, the ice part consists of a cone. The upper radius is the neck radius, r _n; the lower radius is the grain radius r _g. The temperature distribution in the cone is calculated disregarding surface heat and mass exchanges. For steady-state conditions, the heat flux is

(43)

Fig. 8. A simple microstructural cell model consisting of ice and air. The ice is in the shape of a cone of length l with top radius r _n and bottom radius r _g. The air is in the shape of an inverse cone with radii d _n and d _p. No heat is transferred across the ice-air boundary. The model is used to demonstrate the temperature difference between the ice and air given a heat flux in the ice Q _i and air Q _a.

The temperature distribution at position z is found from the above equation:

(44)

The cross-sectional area is

(45)

with

(46)

The integral, Equation (44), can easily be evaluated (for this simple geometry):

(47)

The temperature difference between the upper (z = r _n) and lower (z = 0) cone surfaces can subsequently be found:

(48)

For air, the same calculation procedures apply, substituting d _p for r _g and d _n for r _n. The height of the microstructural cell is in both cases l.

An example calculation is presented in Figure 9. A temperature gradient of 100 K m^-1 is applied over a cell of 1 mm length. The base temperature at z = 0 is 268 K; the neck temperature is 268.1 K. The volumetric content of ice is θ _i = 0.3. The grain radius is r _g = 1mm; the neck radius is r _n = 0.4 mm. Figure 9a displays the temperature distribution in the ice and air over the cell length. Figure 9b displays the temperature difference between the ice and air. The important fact to observe is that a temperature difference of 0.1 K over the length of the grain produces a maximumtemperature difference of 0.03 K between the ice and air. More complicated and realistic geometries will strengthen this effect. If the bond and ice radius are equal, the temperature difference will vanish.

Fig. 9. A temperature gradient of 100 Km^-1 is applied over a microstructural cell of 1mm length. The base temperature at z = 0 is 268K; the neck temperature is 268.1K. The volumetric content of ice is θ_i = 0.3.The grain radius is r _g = 1mm; the neck radius is r _n = 0.4 mm. (a)Temperature distribution in the ice and air over the cell length. (b)Temperature difference. The important fact to observe is that a temperature difference of 0.1K over the length of the grain produces a maximum temperature difference of 0.03 K between the ice and air.

In summary, because the ice lattice consists of a non-uniform, irregular structure, global temperature gradients will induce local temperature differences between the ice and air. If the temperature of the ice lattice is 268 K < T _i < 273 K, these temperature differences must only be of the order ΔT < 0.1 K to produce massive changes in the optimal curvature radius. The fact that larger temperature differences are required at colder temperatures agrees well with the experimental observations made by Reference SokratovSokratov (2001). He noted that the relative importance of temperature gradient was large when the temperature was high, and was negligible for low temperatures.

It appears that temperature-gradient metamorphism can be defined within the framework of minimum entropy production: applied temperature gradients induce slight differences between the ice and air because of the highly irregular and non-uniform granular structure of snow. Depending on the temperature and temperature difference, flat surfaces are required to minimize entropy production.

An observed feature of temperature-gradient metamorphism is that recrystallization rates stop at a limiting density (Reference SokratovSokratov, 2001). Reference Bozhinskiy and LosevBozhinskiy and Losev (1998) state that the limiting density is 350 kg m^-3; Reference MarboutyMarbouty (1980) states a value of 400 kg m^-3. A possible explanation for this phenomenon is that at higher densities, large grains prevail with bond sizes about the size of the grain radius (large a angles). In such a situation, the temperature difference between the ice and air vanishes since the cross-sectional areas of the ice and air are constant. Without irregularities, the temperature gradient cannot induce temperature differences across the interface and, subsequently, temperature- gradient metamorphism stops.

6.3. Comparison of minimum entropy production with isothermal and temperature-gradient metamorphism

A direct comparison can be made between minimum entropy production and mass diffusion theory. The results are presented in Table 4. Diffusion theory divides snow metamorphism into two regimes: isothermal- and temperature- gradient metamorphism. We take Reference McClung and SchaererMcClung and Schaerer (1993) as the standard reference describing these two regimes. We do not consider new snow metamorphism or wet snow metamorphism, for these would require different microstructural or thermodynamic models.

Table 4. The principle of minimum entropy production predicts form changes as a function of ΔT. The widely accepted theory, from Reference McClung and SchaererMcClung and Schaerer (1993), is based on mass diffusion arising from temperature gradients. The two theories are not in contradiction

Table 4 shows the position of minimum entropy production within the framework of established metamorphism theory. It shows that the principle of minimum entropy production is not in contradiction with observations. Equilibrium conditions are likely to exist when the temperature gradients are small and the grains are tightly packed; non-equilibrium conditions are likely to exist under high temperature gradients and when the pore space is large. Isothermal metamorphism thus corresponds to minimum entropy production at equilibrium temperatures, and temperature-gradient metamorphism corresponds to minimum entropy production at non-equilibrium temperatures.

7. Summary and Conclusions

In summary, the theory of minimum entropy production predicts:

Because the controlling physical process of snow metamorphism appears to be interfacial mass transport which is driven by temperature and concentration differences between the ice and air phases, physical snowpack models should be based on a non-equilibrium treatment of heat and mass transfer (Reference Bartelt, Buser and SokratovBartelt and others, 2004).

These results are based on an elementary application of the principle of minimum entropy production. We were unable to derive a more complicated production function in which the control volume contained more than one grain and bond. Only the volumetric constraint condition between the ice and pore air was used in this work. In order to determine grain-size (not only grain curvature), additional (and non-trivial) mass constraints will probably have to be invoked with larger control volumes. The formulation was also based on a rigid, non-deforming ice matrix. Entropy principles including irreversible mechanical deformations (Reference Coleman and NollColeman and Noll, 1963) can and should be applied in future.

Finally, there is a fundamental assumption behind the invocation of the minimum entropy principle. We assume that non-equilibrium states can be deduced from principles of thermodynamic equilibrium (Reference JaynesJaynes, 1980). The practical consequence of this assumption is that it is impossible to say how a system behaves as a function of time. Therefore, it is impossible to state how fast the system attains the state of minimum entropy production. In principle, we cannot even say whether the system reaches this state or not. In our particular case, it is impossible for us to state how fast the grain curvature is changing; we can only suggest the general tendency or direction of the change. This is useful, but it is no substitute for rate equations. We must postulate: if a thermodynamic system can reduce entropy production by changing some variable(s), it will develop towards the minimum entropy production within the given constraints (Reference Glansdorf and PrigogineGlansdorff and Prigogine, 1974). To find the implicit time dependence of entropy production on form will be the next challenge.

According to our analysis, the snow cover is in the neighborhood of equilibrium (Reference PrigoginePrigogine, 1980) because the entropy production principle appears to hold. This will allow us to apply linear non-equilibrium thermodynamics to the snow cover in the future.