# Impact of water vapor diffusion and latent heat on the effective thermal conductivity of snow

Kévin Fourteau, Florent Domine, Pascal Hagenmuller

► To cite this version:

Kévin Fourteau, Florent Domine, Pascal Hagenmuller. Impact of water vapor diffusion and latent heat on the effective thermal conductivity of snow. The Cryosphere, 2021, 15 (6), pp.2739-2755. 10.5194/tc-15-2739-2021. insu-03326421

HAL Id: insu-03326421

https://insu.hal.science/insu-03326421

Submitted on 26 Aug 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Impact of water vapor diffusion and latent heat on the effective thermal conductivity of snow

Kévin Fourteau¹, Florent Domine²,³, and Pascal Hagenmuller¹

¹Univ. Grenoble Alpes, Université de Toulouse, Météo-France, CNRS, CNRM, Centre d'Études de la Neige, Grenoble, France

²Takuvik Joint International Laboratory, Université Laval (Canada) and CNRS-INSU (France), Québec, QC, G1V 0A6, Canada

³Centre d'Études Nordiques (CEN) and Department of Chemistry, Université Laval, Québec, QC, G1V 0A6, Canada

**Correspondence:** Kévin Fourteau (kfourteau@protonmail.com)

Received: 27 October 2020 – Discussion started: 10 November 2020

Revised: 8 April 2021 – Accepted: 26 May 2021 – Published: 18 June 2021

**Abstract.** Heat transport in snowpacks is understood to occur through the two processes of heat conduction and latent heat transport carried by water vapor, which are generally treated as decoupled from one another. This paper investigates the coupling between both these processes in snow, with an emphasis on the impacts of the kinetics of the sublimation and deposition of water vapor onto ice. In the case when kinetics is fast, latent heat exchanges at ice surfaces modify their temperature and therefore the thermal gradient within ice crystals and the heat conduction through the entire microstructure. Furthermore, in this case, the effective thermal conductivity of snow can be expressed by a purely conductive term complemented by a term directly proportional to the effective diffusion coefficient of water vapor in snow, which illustrates the inextricable coupling between heat conduction and water vapor transport. Numerical simulations on measured three-dimensional snow microstructures reveal that the effective thermal conductivity of snow can be significantly larger, by up to about 50% for low-density snow, than if water vapor transport is neglected. A comparison of our numerical simulations with literature data suggests that the fast kinetics hypothesis could be a reasonable assumption for modeling heat and mass transport in snow. Lastly, we demonstrate that under the fast kinetics hypothesis the effective diffusion coefficient of water vapor is related to the effective thermal conductivity by a simple linear relationship. Under such a condition, the effective diffusion coefficient of water vapor is expected to lie in the narrow 100% to about 80% range of

the value of the diffusion coefficient of water vapor in air for most seasonal snows. This may greatly facilitate the parameterization of water vapor diffusion of snow in models.

## 1 Introduction

Thermal conductivity is one of the major physical properties of snow. It governs the magnitude of the thermal energy flux through the snowpack when subjected to a thermal gradient, and it thus plays an integral role in the energy budgets of the ground (Zhang et al., 1996), ice caps and glaciers (Gilbert et al., 2012), and sea ice (Lecomte et al., 2013), as well as in the temperature of the snow surface and therefore in meteorology (Domine et al., 2019). Moreover, variations of thermal conductivity between snow layers impact the temperature gradients at the layer scale and thus in part govern snow metamorphism (Vionnet et al., 2012). In light of its importance for the understanding of snow and environmental physics, snow thermal conductivity has been actively studied and measured for several decades (Yosida et al., 1955; Jaafar and Picot, 1970; Sturm and Johnson, 1992; Morin et al., 2010; Calonne et al., 2011; Riche and Schneebeli, 2013; Domine et al., 2015).

One of the peculiarities of snow is that energy transport does not solely occur through heat conduction. Indeed, when a snowpack is subjected to a thermal gradient, a macroscopic water vapor flux is also present (Sturm and Benson, 1997;
Pinzer et al., 2012). This vapor flux carries latent heat in par-
allel to heat conduction. Several studies have investigated the
influence of vapor transport on the total energy flux through
snow under a thermal gradient. Among others, Sturm and
Johnson (1992) report that the heat transport in snow is char-
acterized by an effective thermal conductivity $\mathbf{K}^{\text{eff}}$ encom-
passing both the effects of heat conduction and vapor trans-
port. In their framework, one can decompose the effective
thermal conductivity as $\mathbf{K}^{\text{eff}} = \mathbf{K}^{\text{cond}} + \mathbf{K}^{\text{vap}}$, where $\mathbf{K}^{\text{cond}}$
is “the hypothetical conductivity in the absence of any va-
por transport” (Sturm and Johnson, 1992) and $\mathbf{K}^{\text{vap}}$ corre-
sponds to the latent heat transported with water vapor. In op-
position to the idea of merging conduction and vapor trans-
port in a single effective thermal conductivity, Calonne et al.
(2011) “recommend purely conductive effects (i.e. conduc-
tion through ice and interstitial air) to be considered sep-
arately from non-conductive processes”, therefore treating
heat conduction as decoupled from vapor transport.

The aim of this article is to provide a simplified analy-
sis of the contribution of latent heat to the thermal energy
flux in snow and notably to quantify the coupling between
these processes at the macroscopic scale. For this we focus on
two limiting cases, considering the kinetics of deposition and
sublimation of water vapor to be either very fast or very slow.
We start by providing theoretical considerations on the rela-
tionship between water vapor transport and the effective ther-
mal conductivity. We then perform numerical simulations to
quantify the contribution of latent heat to the effective ther-
mal conductivity.

## 2 Theory

Let us consider a snow sample of volume $V$, subjected to
a macroscopic thermal gradient denoted $\nabla T^M$ (potentially
accompanied by a macroscopic vapor concentration gradi-
ent $\nabla C^M$). We also make the simplifying assumption that
convection in the pore space does not occur or can be ne-
glected (similarly to Riche and Schneebeli, 2013; Calonne
et al., 2014). Furthermore, let us assume that the sample is
taken large enough to be larger than its representative ele-
mentary volume (REV; Auriault et al., 2010; Calonne et al.,
2014). Moreover, we assume that the sample is small enough
that it does not span several snow layers and that the macro-
scopic thermal and water vapor gradients can be considered
constant over the sample. Under these conditions, the vol-
ume $V$ of snow is representative of the entire snow layer,
and both share the same physical properties. The existence of
this intermediate size between the microscopic and macro-
scopic scales is guaranteed by the separation of scales be-
tween the microscopic and macroscopic descriptions, which
is a necessary condition to treat snow as an equivalent homo-
geneous medium with well-defined physical properties (Au-
riault, 1991; Auriault et al., 2010). An illustration of the mi-

<figure><img src="image_4.png" /></figure>

**Figure 1.** Illustration of the microscopic and macroscopic points of view of snow. At the microscopic scale, snow is composed of an ice space ($\Omega_i$) and a pore space ($\Omega_a$), separated by a boundary ($\Gamma$). Heat conduction occurs through the ice and pore spaces, while vapor diffusion is limited to the pore space. At the macroscopic scale, snow is treated as an equivalent homogeneous medium, with an effective thermal conductivity $\mathbf{K}^{\text{eff}}$ and an effective water vapor diffusion coefficient $\mathbf{D}^{\text{eff}}$, and is subjected to a macroscopic thermal gradient $\nabla T^M$.

macroscopic and macroscopic points of view and of the separa-
tion of scale between them is given in Fig. 1.

The effective thermal conductivity $\mathbf{K}^{\text{eff}}$ of snow relates the macroscopic heat flux $\mathbf{Q}$, transporting thermal energy at the macroscopic scale, to the thermal gradient $\nabla T^M$ with $\mathbf{Q} = -\mathbf{K}^{\text{eff}}\nabla T^M$ (e.g., Yosida et al., 1955; Sturm and Johnson, 1992; Riche and Schneebeli, 2013). Similarly, an effective vapor diffusion coefficient $\mathbf{D}^{\text{eff}}$ can be defined, which relates the macroscopic vapor flux $\mathbf{F}$ to the macroscopic concentration gradient $\nabla C^M$ with $\mathbf{F} = -\mathbf{D}^{\text{eff}}\nabla C^M$ (e.g., Shertzer and Adams, 2018; Fourteau et al., 2021). In snow sciences, it is usually expected that the effective thermal conductivity and vapor diffusion coefficient depends only on the snow microstructure and on the physical properties of the underlying materials (the ice and the air) but not on the macroscopic thermal and water vapor concentration gradients (Yosida et al., 1955; Jaafar and Picot, 1970; Sturm and Johnson, 1992; Colbeck, 1993; Morin et al., 2010; Calonne et al., 2011; Riche and Schneebeli, 2013; Domine et al., 2015). One should however keep in mind that it might not necessarily be true, depending on the nature of the mechanisms at play at the microscopic scale (for instance in the case of a dependence of the sticking coefficient of water molecules onto ice on the local saturation of water vapor, as discussed in Fourteau et al., 2021).

Finally, the effective thermal conductivity is represented
by a 3×3 tensor. However, snow can be considered as a trans-
verse isotropic material (Löwe et al., 2013), and this tensor
is thus fully characterized by two scalar values, namely the vertical and horizontal thermal conductivities. These scalar values are respectively denoted $K_z^{eff}$ and $K_{xy}^{eff}$ in this article to differentiate them from the tensor. Similarly, $D^{eff}$ is represented by a 3 × 3 tensor, but for snow it reduces to a horizontal component and a vertical component. Again, these scalar components are respectively denoted $D_z^{eff}$ and $D_{xy}^{eff}$. We define the normalized effective diffusion coefficient $D^{norm}$ as $D^{norm} = D^{eff}/D_0$, where $D_0$ is the diffusion coefficient of water vapor in air. The definitions and symbols of the variables used for this work are summarized in Appendix A.

In this article, the effective thermal conductivity of snow will be obtained starting from the physics at the microscopic scale. The relevant microscopic physical mechanisms for heat transport are (i) heat conduction in the ice, (ii) heat conduction in the air, (iii) vapor diffusion in the air, and (iv) vapor deposition/sublimation at ice surfaces (Calonne et al., 2014). Moreover, we assume that the physics at the microscopic scale can be treated in a steady state. From our understanding this is justified as the timescale governing the microscopic scale is much shorter than the macroscopic timescale at which snow observations are made. Indeed, Hansen and Foslien (2015) report that the characteristic times at the macroscopic and microscopic scales differ by a factor of $10^6$. Consistent with this, Calonne et al. (2014) report that when expressed in a non-dimensional form, the time derivatives in the heat and mass equations are negligible compared to the flux terms. The microscopic equations governing energy and vapor transport are thus

$$ 
\left\{ 
\begin{array}{l}
\text{div}(-k_i \nabla T_i) = 0 \quad (\Omega_i) \\
\text{div}(-k_a \nabla T_a) = 0 \quad (\Omega_a) \\
\text{div}(-D_0 \nabla c) = 0 \quad (\Omega_a) \\
T_i = T_a \quad (\Gamma) \\
-k_i \nabla T_i \cdot \mathbf{n} = -k_a \nabla T_a \cdot \mathbf{n} - LD_0 \nabla c \cdot \mathbf{n} \quad (\Gamma) \\
-D_0 \nabla c \cdot \mathbf{n} = \alpha v_{\text{kin}}(c - c_{\text{sat}}) \quad (\Gamma)
\end{array}
\right. 
, \quad (1)
$$

where $\Omega_i$, $\Omega_a$, $\Gamma$, and $\mathbf{n}$ represent the ice space, the pore space, the ice-pore interface, and the normal vector to $\Gamma$ pointing toward the ice, respectively. The geometry of the microscopic problem is exemplified in Fig. 1. In Eq. (1), $k_i$ and $k_a$ are the thermal conductivities of ice and air, $T_i$ and $T_a$ are the ice and air temperatures, $D_0$ is the diffusion coefficient of vapor in the air, $c$ is the concentration of vapor in the pores, $c_{sat}$ is the saturation concentration of vapor at the ice interface, $L$ is the latent heat of sublimation of ice, $v_{\text{kin}} = \sqrt{(kT)/(2\pi m)}$ is referred to as the kinetic velocity and is related to the velocity of water molecules in the gas phase (with $k$ the Boltzmann's constant and $m$ the mass of a water molecule), and $\alpha$ is a coefficient less than or equal to unity referred to as the sticking coefficient (or sometimes the accommodation coefficient) of water vapor molecules on ice surfaces. The last equation of the system is referred to as the Hertz-Knudsen equation, and it governs the sublimation

and deposition of water molecules on the ice surfaces (Saito, 1996). The penultimate equation represents the impact of vapor sublimation/deposition on the continuity of the heat flux at the ice-pore interface. Finally, when a sufficiently large thermal gradient is imposed on a snow sample, the variations of $c_{sat}$ due to differences in curvature of the ice surface become negligible compared to variations due to temperature differences (Colbeck, 1983). This large temperature gradient condition corresponds to the regime of temperature gradient metamorphism, usually observed for thermal gradients of $10 \text{ K m}^{-1}$ and above (Sommerfeld and LaChapelle, 1970; Colbeck, 1982). In this case, the saturation concentration of vapor can be treated as a function of the ice surface temperature only.

The system of Eq. (1) shows that there exists a two-way coupling between heat and vapor transport in snow. Indeed, the ice and air temperatures are impacted by the phase change in water vapor and the release/absorption of latent heat, while the water vapor concentration in the pores is impacted by temperature through the value of $c_{sat}$ at the ice surfaces. This implies that the heat flux through a snow sample depends on the sublimation and deposition processes happening in the snow and that the magnitude of the coupling between the heat flux and the vapor transport depends on the kinetics of the adsorption and desorption of water molecules. This kinetics is encapsulated in the parameter $\alpha$ of the Hertz-Knudsen equation. A general treatment of the system of Eq. (1) in the case of an arbitrary $\alpha$ is however out of the scope of this article. In this work, we limit ourselves to two limiting cases, namely very slow (small $\alpha$) and very fast (large $\alpha$) surface kinetics. Moreover, we only focus on quantifying the energy and water vapor fluxes and their associated effective thermal conductivity and effective diffusion coefficient of vapor without deriving the complete macroscopic temperature and vapor equations (contrary to the work of Calonne et al., 2014, for instance). Finally, note that the notion of slow and fast kinetics is related to the notion of kinetics-limited (small $\alpha$) and diffusion-limited (large $\alpha$) metamorphism in snow (e.g., Krol and Löwe, 2016).

## 2.1 The slow kinetics case

In the slow kinetics case, we consider that $\alpha$ is sufficiently small that the sublimation/deposition of water vapor does not strongly impact the temperature field in the snow microstructure. In this case, the coupling between heat and water vapor can be neglected, and snow can be viewed as an inert medium for heat conduction. This slow kinetics case has been treated in detail by Calonne et al. (2011) and case 3 of Calonne et al. (2014). Here, the snow sample is characterized by an effective thermal conductivity $K_{\text{slow}}^{\text{eff}}$ that only accounts for the heat conduction through the ice and the air as if the snow medium were inert for water vapor. The subscript slow, used in $K_{\text{slow}}^{\text{eff}}$ and elsewhere in the paper, is used to emphasize the slow kinetics assumption. Calonne et al. (2011) showed
that the effective thermal conductivity depends on the snow microstructure and on the ice and air thermal conductivities but not on the macroscopic thermal gradient. It can be obtained with microscale numerical simulations of heat conduction, which do not include vapor transport (e.g., Calonne et al., 2011; Riche and Schneebeli, 2013). Following a similar decomposition of the effective conductivity as that reported by Sturm and Johnson (1992), one has $K_{\text{slow}}^{\text{eff}} = K_{\text{slow}}^{\text{cond}}$ and $K_{\text{slow}}^{\text{vap}} = 0$. Note that the fact that $K_{\text{slow}}^{\text{vap}} = 0$ does not imply that the vapor flux in snow is null. A macroscopic vapor flux might be present, but there is simply not enough mass that changes phase and release/absorption of latent heat to meaningfully impact the heat flux in the snow.

## 2.2 The fast kinetics case

In the fast kinetics case, $\alpha$ is sufficiently large that the adsorption/desorption of water molecules is fast enough to impose vapor saturation at the ice-air interface. Mathematically, this case can be treated by letting $\alpha \to \infty$. While this mathematical treatment is purely theoretical (as $\alpha \le 1$), it helps in apprehending the effect of fast kinetics. This case corresponds to the diffusion-limited case, and the Hertz-Knudsen equation is replaced by the saturation of vapor at the ice-air interface. Moreover, it can be shown that in this case the vapor concentration equals its saturation value not only at the ice surface but also throughout the entire pore space (see Yosida et al., 1955, or Fourteau et al., 2021, for demonstrations). The system of Eq. (1) can therefore be rewritten as

$$ 
\left\{ 
\begin{aligned}
\operatorname{div}(-k_i \nabla T_i) &= 0 && (\Omega_i) \\
\operatorname{div}(-k_a \nabla T_a) &= 0 && (\Omega_a) \\
\operatorname{div}(-D_0 \nabla c_{\text{sat}}) &= 0 && (\Omega_a) \\
T_i &= T_a && (\Gamma) \\
-k_i \nabla T_i \cdot \mathbf{n} &= -k_a \nabla T_a \cdot \mathbf{n} - L D_0 \nabla c_{\text{sat}} \cdot \mathbf{n} && (\Gamma)
\end{aligned}
\right. 
. \quad (2) $$

Using the chain rule, one has $\nabla c_{\text{sat}} = \beta \nabla T_a$, where $\beta = \frac{dc_{\text{sat}}}{dT}$. Re-injecting this equality into Eq. (2) yields

$$ 
\left\{ 
\begin{aligned}
\operatorname{div}(-k_i \nabla T_i) &= 0 && (\Omega_i) \\
\operatorname{div}(-k_a \nabla T_a) &= 0 && (\Omega_a) \\
\operatorname{div}(-\beta D_0 \nabla T_a) &= 0 && (\Omega_a) \\
T_i &= T_a && (\Gamma) \\
-k_i \nabla T_i \cdot \mathbf{n} &= -k_a \nabla T_a \cdot \mathbf{n} - \beta L D_0 \nabla T_a \cdot \mathbf{n} && (\Gamma)
\end{aligned}
\right. 
. \quad (3) $$

Multiplying the third line by $L$ and summing with the second line, one finally gets

$$ 
\left\{ 
\begin{aligned}
\operatorname{div}(-k_i \nabla T_i) &= 0 && (\Omega_i) \\
\operatorname{div}(-k_a + \beta L D_0) \nabla T_a) &= 0 && (\Omega_a) \\
T_i &= T_a && (\Gamma) \\
-k_i \nabla T_i \cdot \mathbf{n} &= -(k_a + \beta L D_0) \nabla T_a \cdot \mathbf{n} && (\Gamma)
\end{aligned}
\right. 
. \quad (4) $$

which is the system of equations governing the temperature and heat conduction in a microstructure without any explicit

vapor transport and where the conductivity of the air has been replaced by an apparent conductivity $k_v$, defined as

$$ k_v = k_a + \beta L D_0. \quad (5) $$

An equivalent demonstration of this result was proposed by Yosida et al. (1955). A similar result was also derived by Moyne et al. (1988) for a triphasic medium composed of water vapor, liquid water, and an inert solid phase, and it is consistent with the effective thermal conductivity model of soil proposed by De Vries (1958) and De Vries (1987).

As this system of equations is equivalent to the one of an inert medium with an increased air thermal conductivity, one can show using methods of homogenization (e.g., Auriault et al., 2010; Calonne et al., 2011) that at the macroscopic scale snow can be treated as an equivalent medium with a well-defined tensorial thermal conductivity $K_{\text{fast}}^{\text{eff}}$. This effective thermal conductivity depends on the snow microstructure and on the physical properties of ice, air, and vapor (through $k_i$, $k_a$, and $\beta L D_0$) but not on the macroscopic thermal gradient. Again, the subscript fast is used to stress that we are working under the fast kinetics assumption.

We now investigate the individual contributions of conduction and vapor transport to $K_{\text{fast}}^{\text{eff}}$. The macroscopic heat flux $Q$, which equals the volume average of the microscopic heat flux (Batchelor and Brien, 1977), can be decomposed as

$$ 
\begin{aligned}
Q &= -\frac{1}{V} \left( \int_{V_i} k_i \nabla T_i dV + \int_{V_a} k_v \nabla T_a dV \right) \\
&= -(1-\phi) \frac{1}{V_i} \int_{V_i} k_i \nabla T_i dV - \phi \frac{1}{V_a} \int_{V_a} k_a \nabla T_a dV \\
&\quad - \phi \frac{1}{V_a} \int_{V_a} \beta L D_0 \nabla T_a dV \\
&= -(1-\phi) k_i \langle \nabla T_i \rangle - \phi k_a \langle \nabla T_a \rangle - \phi \beta L D_0 \langle \nabla T_a \rangle,
\end{aligned} 
\quad (6) $$

where $V_i$ and $V_a$ are the ice and air volumes in the snow, $\phi$ is the porosity, and $\langle \nabla T_i \rangle$ and $\langle \nabla T_a \rangle$ stand for the spatial averages of the thermal gradients in the individual ice and air spaces. Note that the average thermal gradient in the ice (respectively the air) is defined by performing the volume average in the ice space only (respectively the air space only) and not in the entire snow volume. The first two terms of the last line of Eq. (6) respectively correspond to the contribution of ice and air heat conduction to the energy flux, while the last term corresponds to an additional contribution of latent heat transported with water vapor. Moreover, recalling that with fast kinetics $\nabla c = \nabla c_{\text{sat}} = \beta \nabla T_a$, the contribution of water vapor is given by
$$ \begin{aligned} \phi \beta L D_0 \langle \nabla T_a \rangle &= L \frac{1}{V} \int_{V_a} \beta D_0 \nabla T_a dV \\ &= L \frac{1}{V} \int_{V_a} D_0 \nabla c dV \\ &= L F, \end{aligned} \quad (7) $$

where $F = \frac{1}{V} \int_{V_a} D_0 \nabla c dV$ is the macroscopic vapor flux (Shertzer and Adams, 2018; Fourteau et al., 2021) which is linked to the macroscopic vapor gradient $\nabla C^M$ through an effective diffusion coefficient $D_{\text{fast}}^{\text{eff}}$ such that $F = -D_{\text{fast}}^{\text{eff}} \nabla C^M$ (Calonne et al., 2014; Fourteau et al., 2021). Since in the fast kinetics case vapor is at saturation throughout the pore space, the macroscopic vapor gradient is related to the macroscopic temperature gradient $\nabla T^M$ by $\nabla C^M = \beta \nabla T^M$. Therefore, the macroscopic vapor flux is given by

$$ F = -\beta D_{\text{fast}}^{\text{eff}} \nabla T^M. \quad (8) $$

The effective thermal conductivity of snow can thus be decomposed in

$$ K_{\text{fast}}^{\text{eff}} = K_{\text{fast}}^{\text{cond}} + \beta L D_{\text{fast}}^{\text{eff}}, \quad (9) $$

where $K_{\text{fast}}^{\text{cond}}$ is the contribution due to heat conduction in the ice and air, and $\beta L D_{\text{fast}}^{\text{eff}} = K_{\text{fast}}^{\text{vap}}$ is the contribution of latent heat transported with water vapor. Contrary to the slow kinetics case, we now find that latent heat has an impact on the thermal properties of snow and that the vapor flux directly contributes to the effective thermal conductivity. A similar version of Eq. (9), with the contribution of vapor transport being $\beta L D_{\text{fast}}^{\text{eff}}$, has notably been reported by Jordan (1991) and Sturm and Johnson (1992), which directly consider energy balance and transport at the macroscopic scale and suppose the existence of well-defined (i.e., independent of the macroscopic thermal gradient) $K_{\text{fast}}^{\text{cond}}$ and $D_{\text{fast}}^{\text{eff}}$.

It is important to note that $K_{\text{fast}}^{\text{cond}}$ in Eq. (9) is different from the thermal conductivity when latent heat effects are neglected, i.e., $K_{\text{fast}}^{\text{cond}} \neq K_{\text{slow}}^{\text{cond}}$ (Moyne et al., 1988). Indeed, the heat conduction in the microstructure is determined by the distribution of the local thermal gradients in the two phases, and latent heat modifies the microscopic thermal gradients compared to the case without latent heat and thus modifies the heat conduction. The presence of latent heat increases the apparent thermal conductivity of the pore space and thus reduces the thermal conductivity contrast between the two phases. In turn, this reduced thermal contrast increases the average temperature gradient of the ice phase (the highly conducting phase) and decreases the average temperature gradient of the gas phase (the poorly conducting phase). The increase in heat conduction in the ice is larger than the decrease in the air, and the contribution of heat conduction alone is therefore greater with the presence of latent heat than

without. Such an effect is illustrated and quantified with numerical simulations in Sect. 3.1.

Finally, we want to point out that in the fast kinetics case, the effective thermal conductivity and the effective water vapor diffusion coefficient are linearly related. Indeed, starting from the fact that the effective diffusion coefficient is given by the ratio of the magnitude of the vapor flux over the magnitude of the vapor concentration gradient, one has

$$ \begin{aligned} D_{\text{fast}}^{\text{eff}} &= \frac{\| \frac{1}{V} \int_{V_a} D_0 \nabla c dV \|}{\| \nabla C \|} = \frac{\| D_0 \beta \frac{1}{V} \int_{V_a} \nabla T_a dV \|}{\| \beta \nabla T^M \|} \\ &= D_0 \frac{V_a}{V} \frac{\| \frac{1}{V_a} \int_{V_a} \nabla T_a dV \|}{\| \nabla T^M \|} = D_0 \phi \frac{\| \langle \nabla T_a \rangle \|}{\| \nabla T^M \|}, \end{aligned} \quad (10) $$

where we used the facts that $\nabla c = \beta \nabla T_a$ and $\nabla C^M = \beta \nabla T^M$ and where $D_{\text{fast}}^{\text{eff}}$ either stands for the vertical or for the horizontal effective diffusion coefficient depending on the orientation of the macroscopic thermal gradient. The ratio of the average thermal gradient in the pore space over the macroscopic thermal gradient is governed by the effective thermal conductivity and the thermal conductivity of the ice and the air. Indeed, we have

$$ K_{\text{fast}}^{\text{eff}} \nabla T^M = (1 - \phi) k_i \langle \nabla T_i \rangle + \phi k_v \langle \nabla T_a \rangle \quad (11) $$

and

$$ \nabla T^M = (1 - \phi) \langle \nabla T_i \rangle + \phi \langle \nabla T_a \rangle, \quad (12) $$

where similarly $K_{\text{fast}}^{\text{eff}}$ either stands for the vertical or for the horizontal effective thermal conductivity depending on the orientation of the macroscopic thermal gradient. Equation (11) follows from the definition of the macroscopic energy flux (Batchelor and Brien, 1977), and Eq. (12) follows from the application of Stokes theorem and has notably been previously reported by Hansen and Foslien (2015). Combining Eqs. (11) and (12) we have

$$ \phi \langle \nabla T_a \rangle = \frac{k_i - K_{\text{fast}}^{\text{eff}}}{k_i - k_v} \nabla T^M. \quad (13) $$

Finally, injecting Eq. (13) into Eq. (10) we have

$$ D_{\text{fast}}^{\text{eff}} = D_0 \frac{k_i - K_{\text{fast}}^{\text{eff}}}{k_i - k_v}. \quad (14) $$

Since the effective thermal conductivity is larger than the conductivity of the least conducting phase, i.e., $K_{\text{fast}}^{\text{eff}} > k_v$, one finds that $D_{\text{fast}}^{\text{eff}} \le D_0$, as reported by Giddings and LaChapelle (1962) and Fourteau et al. (2021). The linear relationship between the effective thermal conductivity and the normalized effective water vapor diffusion coefficient at 263 K is displayed in Fig. 2 for effective thermal conductivities $K_{\text{fast}}^{\text{eff}}$ ranging from $k_v = 0.0336$ to $0.5 \text{ W K}^{-1} \text{ m}^{-1}$, as typically encountered with seasonal snow (e.g., Sturm et al., 1997; Calonne et al., 2011; Riche and Schneebeli, 2013, and
the numerical values computed in Sect. 3.2 of this paper). The application of Eq. (14) in Fig. 2 reveals that the normalized effective diffusion coefficient ranges from 1 to about 0.8 for most seasonal snow. Moreover, as the macroscopic vapor flux decreases with slower kinetics (Pinzer et al., 2012; Fourteau et al., 2021), this curve represents an upper limit for the effective diffusion coefficient.

## 2.3 Intermediate cases

Numerous works indicate that $\alpha$ depends on temperature, the local vapor saturation, and the crystallographic properties of the underlying ice surface (e.g., Saito, 1996; Libbrecht and Rickerby, 2013), but for snow it remains unclear what value or expression should be used for $\alpha$ in the Hertz–Knudsen equation (Legagneux and Domine, 2005). However, a recent study suggests that the very slow kinetics and fast surface kinetics cases correspond to the minimum and maximum macroscopic vapor diffusion in snow, respectively (Fourteau et al., 2021). We can therefore expect the energy flux $Q$ to be maximal in the fast kinetics case since this corresponds to the situation with maximal vapor flux and the fastest adsorption and desorption of water molecules onto the ice surface. Similarly, the energy flux is minimal in the slow kinetics case as latent heat effects are absent in this case. The energy flux in snow $Q$ can thus be bounded by the slow kinetics and the fast kinetics cases:

$$K_{\text{slow}}^{\text{eff}} \| \nabla T^{\text{M}} \| \le Q \le K_{\text{fast}}^{\text{eff}} \| \nabla T^{\text{M}} \|, \quad (15)$$

where this inequality applies both for the vertical and horizontal components of the effective thermal conductivities, depending on the orientation of the macroscopic thermal gradient.

## 3 Numerical simulations

To exemplify and quantify the points raised in Sect. 2, we performed finite element simulations of steady-state thermal conduction through several snow microstructures obtained experimentally with computed microtomography. The simulations were performed using the open-source Elmer FEM software (Malinen and Råback, 2013) and the readily available solver for the heat equation.

In each simulation, the temperatures of two opposite sides of the microstructure were imposed in order to obtain a thermal gradient of $50 \text{ K m}^{-1}$, with adiabatic conditions on the remaining sides. Similarly to Riche and Schneebeli (2013), the effective thermal conductivities are estimated by computing the ratio of the macroscopic heat flux $Q$ to the macroscopic thermal gradient $\nabla T^{\text{M}}$. The macroscopic heat flux $Q$ is computed as the volume average of the microscopic heat fluxes using the ParaView software. Note that the chosen value of $50 \text{ K m}^{-1}$ for the imposed gradient is purely arbitrary and does not have an impact on our results as the resulting effective thermal conductivity does not depend on the magnitude of the gradient.

In order to test the influence of temperature on the effective thermal conductivity of snow, the simulations were run for different mean temperatures, ranging from 223 to 273 K. The temperature dependence of the thermal conductivities of ice and air ($k_i$ and $k_a$) were respectively taken from Lide (2006) (based on Slack, 1980) and Kadya et al. (1985). The parameter $\beta = \frac{dC_{\text{sat}}}{dT}$ was obtained by assuming that water vapor follows the Clausius–Clapeyron and ideal gas laws (Eq. 11 of Fourteau et al., 2021). We set the diffusion coefficient of water vapor in air to $D_0 = 2 \times 10^{-5} \text{ m}^2 \text{ s}^{-1}$ (Calonne et al., 2014) and the latent heat of sublimation of ice $L = 28 \times 10^5 \text{ J kg}^{-1}$ (Lide, 2006), independent of temperature. Finally, we assume that the density of ice $\rho_{\text{ice}}$ is constant with temperature and equal to $917 \text{ kg m}^{-3}$ (Calonne et al., 2014). The densities of the samples reported in this article are computed using $\rho_{\text{ice}}$, and the ice volume fractions are deduced from the tomography images.

For the different microstructures and mean temperatures, two types of simulations were performed, one in which we assumed no impact of latent heat on the heat conduction (thus obtaining $K_{\text{slow}}^{\text{eff}}$), and the other in which we increased the apparent thermal conductivity of air by a $\beta L D_0$ term (thus obtaining $K_{\text{fast}}^{\text{eff}}$). Moreover, as seen in Sect. 2.2 with Eq. (14), under the fast kinetics assumption the effective diffusion coefficient of water vapor $D_{\text{fast}}^{\text{eff}}$ can be directly obtained from the effective thermal conductivity. Finally, we recall that the normalized effective diffusion coefficient is defined as the ratio of the effective diffusion coefficient with the diffusion coefficient of water vapor in free air, i.e., $D^{\text{norm}} = D^{\text{eff}}/D_0$.

In total we used 34 measured snow microstructures covering several types of seasonal snow. The particular snow types used are, according to the terminology of Fierz et al. (2009), precipitation particles (PP), decomposing and fragmented precipitation particles (DF), rounded grains (RG), faceted crystals (FC), depth hoar (DH), and melt forms (MF). We used sample sizes larger than the REV sizes reported by Calonne et al. (2011). The tetrahedral meshes used in the finite element simulations were produced using the Computational Geometry Algorithms Library (CGAL) and contain between 20 and 90 million elements. The samples are described in the Supplement and include both previously published snow samples (from Hagenmuller et al., 2016, 2019; Peinke et al., 2020) and new snow samples.

### 3.1 Effect of temperature on the effective thermal conductivity

In this section we analyze the influence of the mean temperature on the effective thermal conductivity. For simplicity, we limit ourselves to vertical temperature gradients and thus only deal with vertical effective thermal conductivities and vertical diffusion coefficients of water vapor. As all scalar
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Normalized effective water vapor diffusion coefficient as a function of the effective thermal conductivity under the fast kinetics hypothesis at 263 K. The shaded area covers the typical range of thermal conductivity values (from $k_v$ up to 0.5 W K$^{-1}$ m$^{-1}$) and the corresponding range of the normalized effective diffusion coefficient of water vapor (from 1 to about 0.8). At 263 K, $k_v = k_a + \beta LD_0 = 0.0336$ W K$^{-1}$ m$^{-1}$.</figcaption></figure>

components are vertical, we do not use the subscript z in order to lighten the notation.

The temperature dependence of $K^{\text{eff}}$ is due to the temperature dependence of the underlying materials. Indeed, an increasing temperature results in the decrease in the ice thermal conductivity $k_i$ and the increase in the apparent thermal conductivity of the air $k_v$ due to both the increase in the intrinsic thermal conductivity of air $k_a$ and the increase in the contribution of water vapor latent heat $\beta LD_0$, all displayed in Fig. 3.

Furthermore, we define for our analysis

$$ 
\begin{aligned}
K^{\text{ice}} &= (1-\phi)k_i \frac{\|\nabla T_i\|}{\|\nabla T\|}, \\
K^{\text{air}} &= \phi k_a \frac{\|\nabla T_a\|}{\|\nabla T\|},
\end{aligned} 
\qquad (16) 
$$

where $K^{\text{ice}}$ (not to be mistaken with $k_i$; see Appendix A) corresponds to the contribution of the ice heat conduction to the total effective thermal conductivity and $K^{\text{air}}$ (not to be mistaken with $k_a$) to the contribution of the air heat conduction. We have by construction

$$ K^{\text{cond}} = K^{\text{ice}} + K^{\text{air}}, \qquad (17) $$

where $K^{\text{cond}}$ is the vertical component of $K^{\text{cond}}$. With our simulations, the values of $K^{\text{cond}}$, $K^{\text{ice}}$, and $K^{\text{air}}$ are computed directly from $K^{\text{eff}}$, $k_i$, $k_a$, and $k_v$ and using Eqs. (11) and (12).

We first focus on only two snow samples: a low-density sample and a high-density sample. The low-density sample is composed of decomposing and fragmented precipitation particles (DF) with a density of 125 kg m$^{-3}$ and a specific surface area of 40 m$^2$ kg$^{-1}$. The high-density sample is composed of melt forms (MF) with a density of 380 kg m$^{-3}$ and

a specific surface area of 5 m$^2$ kg$^{-1}$. The three-dimensional microstructures of both samples are displayed in Fig. 4. The results of the finite element simulations for these two samples are reported in Fig. 5.

We start by analyzing the low-density sample (left column of Fig. 5). Under the fast kinetics hypothesis, the effective thermal conductivity of the low-density sample shows an exponential-like increase with increasing temperature. This increase in $K^{\text{eff}}$ is due to the combined effects of (i) an increase in $K^{\text{vap}}$, the vertical component of the contribution of latent heat transport, and (ii) an increase in $K^{\text{cond}}$, the heat conduction through the ice and air spaces. The increase in $K^{\text{cond}}$ is principally due to the increase in $K^{\text{ice}}$, the heat conduction in the ice. With increasing temperature, the increase in the apparent thermal conductivity of air reduces the contrast between the two phases, and the average ice thermal gradient increases. This increase more than offsets the decrease in the ice thermal conductivity $k_i$, and the net effect is an increase in $K^{\text{ice}}$. Under the slow kinetics hypothesis, however, the effective thermal conductivity only barely decreases over the range of temperature studied, consistent with the results of Calonne et al. (2011). In this case, the increase in the thermal conductivity of the air is not as pronounced, and the increase in the thermal gradient in the ice does not compensate for the decrease in the ice thermal conductivity. Overall $K^{\text{ice}}$ decreases with temperature in the slow kinetics case.

Contrary to the low-density sample, the high-density sample (right column of Fig. 5) shows a decrease in the effective thermal conductivity in the fast kinetics case. The increase in $K^{\text{vap}}$ with temperature does not counter the decrease in $K^{\text{cond}}$. This decrease in $K^{\text{cond}}$ can be attributed to the decrease in $K^{\text{ice}}$ with temperature. Here, the increase in the ice thermal gradient is not large enough to offset the decrease in $K^{\text{ice}}$.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Temperature dependence of the thermal conductivity of ice ($k_i$ in blue), of the thermal conductivity of air ($k_a$ in orange), of the contribution of latent heat to the apparent thermal conductivity of the air ($\beta LD_0$ in green), and of the apparent thermal conductivity of air including latent heat effect ($k_v = k_a + \beta LD_0$ in black). Note the break in the y axis.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Tetrahedral meshes (ice phase only) of the low-density DF sample (a) and high-density MF sample (b).</figcaption></figure>

crease in the ice thermal conductivity $k_i$, and overall $K^{\text{ice}}$ decreases. Under the slow kinetics hypothesis, the effective thermal conductivity of the high-density sample decreases with temperature slightly more rapidly than in the fast kinetics case.

For both samples, the difference between $K_{\text{slow}}^{\text{eff}}$ and $K_{\text{fast}}^{\text{eff}}$ is maximal near the melting point, where it reaches more than 50% for low-density snow. Moreover, neglecting the effect of water vapor transport on heat conduction under the fast kinetics case can lead to an underestimation of about 20% of the conduction contribution. Thus in the fast kinetics case, the effect of latent heat can only reasonably be neglected for low temperatures or high-density samples.

In order to better quantify the difference between the fast and slow kinetics cases, we computed the vertical effective thermal conductivity for the totality of our 34 snow samples under both hypotheses, at 248 and 273 K. The ratios of the

effective thermal conductivity in the fast kinetics case over the slow kinetics case are displayed in Fig. 6. They confirm that the relative difference is more important for low-density snow and for higher temperatures. Near the melting point (Fig. 6b), the fast kinetics effective thermal conductivity is between 10% and 50% higher than in the slow kinetics case. For colder snow at 248 K (Fig. 6a), the relative difference is less marked and ranges from 1% to 10%. When expressed in absolute terms, however, the difference between the fast and slow kinetics thermal conductivity is more marked for high-density snow.

Finally, Fig. 7 shows the variation of the vertical normalized effective diffusion coefficient of water vapor with temperature under the fast kinetics hypothesis and for the low and high-density samples shown in Fig. 4. The numerical values are consistent with the recent study of Fourteau et al. (2021), who obtained effective diffusion coefficients using
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Vertical effective thermal conductivity ($K^{\text{eff}}$), with the contributions of ice heat conduction ($K^{\text{ice}}$), air heat conduction ($K^{\text{air}}$), and vapor transport ($K^{\text{vap}}$) for a low-density snow sample (a, c) and a high-density snow sample (b, d), as well as in the fast (upper line) and slow (lower line) kinetics cases. $K^{\text{cond}}$ stands for the purely conductive part of $K^{\text{eff}}$ and is given by $K^{\text{cond}} = K^{\text{ice}} + K^{\text{air}}$.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Ratio of the fast kinetics over the slow kinetics vertical effective thermal conductivity for various snow samples as a function of density. Computations performed at 248 K in (a) and 273 K in (b). Note the different y-axis scales in both panels.</figcaption></figure>

finite element simulations explicitly representing vapor diffusion in the pores. Figure 7 reveals a slight decrease in the effective diffusion coefficient with temperature for both low- and high-density snow. This can be explained by the decrease in the air thermal gradient as the apparent conductivity of air increases with temperature. A lower air temperature gradient

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Vertical normalized effective diffusion coefficient $D_{\text{fast}}^{\text{norm}}$ of a low-density snow sample (blue) and a high-density snow sample (orange) as a function of temperature and under the hypothesis of fast kinetics. Note the break in the y axis.</figcaption></figure>

leads to a lower vapor concentration gradient in the pores and thus to a lower vapor flux and a lower $D_{\text{fast}}^{\text{norm}}$.

## 3.2 Effective thermal conductivity and diffusion coefficient as a function of snow density

The slow kinetics effective thermal conductivities of snow samples covering a broad range of densities and microstructures have been reported by Calonne et al. (2011) and Riche and Schneebeli (2013). Similarly, numerical values of the effective diffusion coefficient of water vapor in snow under limited kinetics have been provided by Calonne et al. (2014). Here, we provide numerical estimates of the effective thermal conductivities and effective diffusion coefficients of a broad range of snow samples, this time under the fast kinetics hypothesis. For each sample we computed the vertical and horizontal effective thermal conductivities and water vapor
effective diffusion coefficients, at five different temperatures (223, 248, 263, 268, and 273 K). The thermal conductivities and diffusion coefficients of each simulated sample are available in the Supplement.

The thermal conductivities computed at 263 K are displayed in Fig. 8 as a function of density. Similar to the work of Calonne et al. (2011) and Riche and Schneebeli (2013), we observe that density and thermal conductivity are well correlated, with denser snow samples presenting higher thermal conductivity values. For the low-density samples, for which the conduction of air plays a determinant role in the effective thermal conductivity, we report thermal conductivity values higher than the polynomial fits of Calonne et al. (2011) and Riche and Schneebeli (2013), both based on the slow kinetics hypothesis. This difference can be explained by the increased apparent thermal conductivity of the air due to latent heat effects. At higher density, our data lie above the reported data and polynomial fit of Calonne et al. (2011). As the relative difference between the fast and slow kinetics cases is small for high-density samples, one can expect the slow and fast kinetics simulations to yield similar values for high-density samples. The scatter between our values and the study of Calonne et al. (2011) is likely due to the inherent variability between snow samples, even for equal densities. Note that the fit proposed by Riche and Schneebeli (2013) was based on faceted crystals and depth hoar snow only and at 253 K. On the contrary the fit proposed by Calonne et al. (2011) was based on their entire sample set at 271 K.

We adjusted second-order polynomial functions to derive parameterizations of thermal conductivity as a function of density and this for each of the five temperatures studied. Our parametrization for the vertical effective thermal conductivity at 263 K is displayed as a solid line in Fig. 8. The parameterizations of the vertical effective thermal conductivity for the five different temperatures are given by

$$K_z^{\text{eff}} = \begin{cases} 2.564 \left( \frac{\rho}{\rho_{\text{ice}}} \right)^2 - 0.059 \frac{\rho}{\rho_{\text{ice}}} + 0.0205 & \text{for } T = 223 \, \text{K} \\ 2.172 \left( \frac{\rho}{\rho_{\text{ice}}} \right)^2 + 0.015 \frac{\rho}{\rho_{\text{ice}}} + 0.0252 & \text{for } T = 248 \, \text{K} \\ 1.985 \left( \frac{\rho}{\rho_{\text{ice}}} \right)^2 + 0.073 \frac{\rho}{\rho_{\text{ice}}} + 0.0336 & \text{for } T = 263 \, \text{K} \\ 1.883 \left( \frac{\rho}{\rho_{\text{ice}}} \right)^2 + 0.107 \frac{\rho}{\rho_{\text{ice}}} + 0.0386 & \text{for } T = 268 \, \text{K} \\ 1.776 \left( \frac{\rho}{\rho_{\text{ice}}} \right)^2 + 0.147 \frac{\rho}{\rho_{\text{ice}}} + 0.0455 & \text{for } T = 273 \, \text{K} \end{cases} \quad (18)$$

where $\rho/\rho_{\text{ice}}$ is the volume fraction of ice, and the constant terms in the polynomial equations correspond to $k_v$. Similar parameterizations for the horizontal thermal conductivity and for the geometric mean of the vertical and horizontal thermal conductivities are available in the Supplement of this article. This parametrization can be extended to other temperatures by first computing the thermal conductivity at the desired density for the five proposed temperatures and then performing an interpolation to the desired temperature.

These vertical effective thermal conductivity parameterizations, displayed in Fig. 9, show a decrease in the slopes of

the thermal conductivity versus density curves with increasing temperatures. This is consistent with the observations made in Sect. 3.1 that the thermal conductivity of the low-density sample increases with temperature, while the thermal conductivity of the high-density sample decreases with temperature. The transition between these two behaviors lies around 350 to 400 kg m⁻³. Note that Calonne et al. (2019) report that a similar transition between the low- and high-density samples also exists under limited kinetics but occurs at a much lower density of about 100 kg m⁻³.

Finally, the estimated normalized effective diffusion coefficients of water vapor are displayed in Fig. 10 as a function of density at 263 K. The normalized effective diffusion coefficients obtained by application of Eq. (14), together with the polynomial fit of the vertical effective thermal conductivity, are shown as a solid black line in Fig. 10. The normalized effective diffusion coefficient decreases with density and mostly remains in the 0.98 to 0.8 range. Notably, detailed seasonal snow models working under the fast kinetics assumption could thus make the reasonable simplifying assumption that $D^{\text{norm}} = 0.90$, independent of snow type. This would result in a less than 10% error on the effective diffusion coefficient of water vapor.

# 4 Discussion

## 4.1 Does the fast kinetics hypothesis apply for heat and mass transport in snow?

This paper studied two limiting cases, considering either that the kinetics of water vapor deposition/sublimation is sufficiently fast to impose saturated water vapor at the ice interface (very large $\alpha$) or that the kinetics is sufficiently slow so that latent heat does not have an impact on either the temperature gradients or the heat conduction in the snow microstructure (very small $\alpha$). It remains however unclear if one of these two limiting cases applies for snow modeling. For example, based on the observation of snow crystal growth with computed tomography, Krol and Löwe (2016) suggest that isothermal metamorphism is slightly better represented by a slow kinetics, while temperature gradient metamorphism data appear consistent with fast kinetics.

As seen in Sect. 3.1, the effective thermal conductivity of low-density snow displays a fundamentally different dependence on temperature depending on whether the slow or the fast kinetics hypothesis applies. In the slow kinetics case, the effective thermal conductivity slightly decreases with increasing temperature, while it increases in the fast kinetics case. Using the needle probe method, Sturm and Johnson (1992) measured the variation of the effective thermal conductivity of a low-density sample of depth hoar with temperature. Even though recent studies have highlighted a potential bias of the needle probe method when used with snow (Calonne et al., 2011; Riche and Schneebeli, 2013), this re-
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Effective thermal conductivity of snow as a function of density under the fast kinetics assumption at 263 K. The horizontal bar of a symbol marks the horizontal effective thermal conductivity value of a snow sample, while the tip of the vertical bar marks its vertical value. Snow classification according to Fierz et al. (2009). Black dot: apparent thermal conductivity of air at 263 K. Solid black line: second-order polynomial fit of the vertical effective thermal conductivity. Dashed grey line: polynomial fit proposed by Calonne et al. (2011) under the slow kinetics assumption at 271 K. Dotted grey line: polynomial fit proposed by Riche and Schneebeli (2013) under the slow kinetics assumption at 253 K.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Temperature dependence of the vertical effective thermal conductivity parameterizations under the fast kinetics hypothesis.</figcaption></figure>

ported bias does not have an impact on the trend of thermal conductivity measured at different temperatures in similar snow samples, as performed by Sturm and Johnson (1992). These data can thus be expected to reflect the variation of the effective thermal conductivity with temperature. These measurements, displayed in Fig. 7 of Sturm and Johnson (1992), clearly indicate an exponential-like increase in thermal conductivity with temperature, consistent with the fast kinetics hypothesis but not with the slow kinetics hypothesis.

The differences between the slow and fast kinetics cases on the effective diffusion coefficient of water vapor were also

studied by Fourteau et al. (2021). While direct measurements of the effective diffusion coefficient are difficult and should therefore be analyzed with caution, the reported experimental values of Sokratov and Maeno (2000) report an average normalized diffusion coefficient of 0.64 for snow densities of about 475 kg m⁻³, while Calonne et al. (2014) report a value of about 0.35 for limited kinetics, and extrapolation of our results suggests a value of about 0.70 under the fast kinetics. Finally, the numerical simulations of Fourteau et al. (2021) indicate that for water vapor diffusion the transition between the slow and fast kinetics regimes occurs for stick-
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Normalized effective diffusion coefficient as a function of density under the fast kinetics assumption at 263 K. Snow classification according to Fierz et al. (2009). Solid black line: normalized effective diffusion coefficient deduced from the application of Eq. (14) with the effective thermal conductivity polynomial fit of Fig. 8.</figcaption></figure>

ing coefficients $\alpha$ around $10^{-3}$. Based on data by Libbrecht (2006), Kaempfer and Plapp (2009) report that $\alpha$ is likely to be within the $10^{-3}$ to $10^{-1}$ range, thus within the fast kinetics regime.

All the above reasons suggest that the effective thermal conductivity and diffusion coefficient of water vapor in snow could be well represented under the fast kinetics hypothesis, at least during temperature gradient metamorphism. Further experimental work should be performed to confirm that the fast kinetics assumption generally applies for modeling mass and heat transport in snow and to highlight its potential limitations. Also, the derivation of a theoretical model able to describe heat and mass transfer with arbitrary surface kinetics would allow one to investigate intermediate kinetics in an effort to ultimately select the best modeling assumptions for snow. At the same time, this model could be formulated to explicitly take into account macroscopic convection as this phenomenon has been observed in sub-arctic shallow snow-packs (Trabant and Benson, 1972; Sturm and Johnson, 1991). Its derivation could be achieved using standard homogenization methods, such as the two-scale asymptotic expansion (e.g., Municchi and Icardi, 2020) or volume averaging methods (e.g., Whitaker, 1977).

## 4.2 The coupling of heat conduction with vapor transport

We showed that in the fast kinetics case, the pure conduction part $K^{\text{cond}}$ of the effective thermal conductivity is influenced by the presence of water vapor and its latent heat. Therefore, the definition of $K^{\text{cond}}$ given by Sturm and Johnson (1992), i.e., that it is “the hypothetical conductivity in the absence of

any vapor transport”, should be clarified to emphasize that $K^{\text{cond}}$ corresponds to the pure conduction occurring through the ice and pore spaces but in response to the actual microscopic thermal gradients that are influenced by the latent heat effects. Furthermore, the dependence of the pure conduction part on temperature is different from what would be expected from variations of the ice and air thermal conductivity only. This means that under the fast kinetics hypothesis a strong two-way coupling exists between heat conduction and water vapor transport, and the heat conduction process cannot be fully considered without latent heat processes. One should therefore be careful when treating heat conduction as decoupled from vapor transport (e.g., Calonne et al., 2011; Riche and Schneebeli, 2013). While this approximation is justified if the effects of latent heat are small, one should be aware of the potential limit of this approximation. Finally, in such a case it is not possible to experimentally decouple the measurement of $K^{\text{cond}}$ from $K^{\text{vap}}$ by performing measurements at low temperature (where $K^{\text{vap}} \simeq 0$). The inferred value of $K^{\text{cond}}$ at low-temperature does not hold at higher temperatures, at which the effect of latent heat is no longer negligible and thus has an impact on $K^{\text{cond}}$. A similar conclusion was reached by Moyne et al. (1988) for the thermal conductivity of a humid triphasic medium.

## 5 Conclusions

This paper investigates the effective thermal conductivity of snow and its relationship to the diffusion of water vapor and its associated latent heat. Using theory, we show that the kinetics of the sublimation and deposition processes at the ice
surfaces plays a significant role on the transport of heat in snow. In particular, if the kinetics is slow, we recall that snow can be treated as an inert medium and that heat transport only occurs through conduction in the ice and in the air. In contrast, if the kinetics is fast, vapor transport and latent heat effects become an integral part of heat transport, and the effective thermal conductivity of snow is composed of a purely conductive term and a term proportional to the water vapor diffusivity. Moreover, we show that under the latter hypothesis there is a simple linear relationship between the effective diffusion coefficient of water vapor in snow and the effective thermal conductivity. Since the effective thermal conductivity of snow rarely exceeds $0.5 \text{ W K}^{-1} \text{ m}^{-1}$, we conclude that under fast kinetics the normalized effective diffusion coefficient of water vapor ranges between 1 and about 0.80 for most seasonal snow.

We complemented this theoretical work by finite element simulations of heat conduction through snow microstructures obtained with computed tomography. The simulations were performed on a total of 34 samples, covering the typical seasonal snow types, under both the slow and fast kinetics hypotheses and for temperatures ranging from 223 to 273 K. The simulations were performed on large samples in order to ensure the representativeness of the results.

Using this new set of numerical simulations, we show that the influence of vapor transport in the fast kinetics case can lead to a significant increase in the effective thermal conductivity compared to the slow kinetics case, up to 50% for low-density snow near the melting point. Moreover, we show that under the fast kinetics hypothesis the purely conductive term of the effective thermal conductivity is influenced by the presence of water vapor and differs from the effective thermal conductivity in the absence of any vapor transport. Indeed, sublimation and deposition processes modify the ice surface temperature through latent heat effect, therefore affecting thermal gradients throughout the snow microstructure. This observation illustrates the coupled nature of heat and water vapor transport in snow, where one cannot be fully understood and quantified without the other. We also compared our numerical simulations to published experimental data of the dependence of the effective thermal conductivity of snow on temperature. This suggests that the fast kinetics option might be well suited to model heat and mass transport in snow during temperature gradient metamorphism. Finally, we provide our new numerical values of the effective thermal conductivity and of the effective diffusion coefficient of water vapor under the fast kinetics hypothesis, derived from snow microstructures measured with computed tomography, as well as parameterizations with snow density. These new data and parameterizations are primarily meant to be used in detailed snow physics models.
Appendix A

Symbols and definitions of the major variables used in this article. The convention followed is that a tensorial variable is denoted with a bold capitalized letter and its scalar components with the non-bold capitalized letter.

<table><thead><tr><th>Symbol</th><th>Definition</th></tr></thead><tbody><tr><td>K<sup>eff</sup></td><td>Effective thermal conductivity of snow</td></tr><tr><td>K<sup>eff</sup><sub>xy</sub></td><td>Horizontal component of the effective thermal conductivity of snow</td></tr><tr><td>K<sup>eff</sup><sub>z</sub></td><td>Vertical component of the effective thermal conductivity of snow</td></tr><tr><td>K<sup>eff</sup></td><td>Scalar component of the effective thermal conductivity of snow (either horizontal or vertical)</td></tr><tr><td>K<sup>cond</sup></td><td>Purely conductive part of the thermal conductivity of snow</td></tr><tr><td>K<sup>vap</sup></td><td>Vapor transport part of the thermal conductivity of snow</td></tr><tr><td>K<sup>air</sup></td><td>Contribution of air heat conduction to the vertical effective thermal conductivity of snow (Eq. 16)</td></tr><tr><td>K<sup>ice</sup></td><td>Contribution of ice heat conduction to the vertical effective thermal conductivity of snow (Eq. 16)</td></tr><tr><td>D<sup>eff</sup></td><td>Effective diffusion coefficient of water vapor in snow</td></tr><tr><td>D<sup>eff</sup><sub>xy</sub></td><td>Horizontal component of the effective diffusion coefficient of water vapor in snow</td></tr><tr><td>D<sup>eff</sup><sub>z</sub></td><td>Vertical component of the effective diffusion coefficient of water vapor in snow</td></tr><tr><td>D<sup>eff</sup></td><td>Scalar component of the effective diffusion coefficient of water vapor in snow (either horizontal or vertical)</td></tr><tr><td>D<sup>norm</sup></td><td>Normalized effective diffusion coefficient of water vapor in snow</td></tr><tr><td>D<sub>0</sub></td><td>Diffusion coefficient of water vapor in air</td></tr><tr><td>·slow</td><td>Subscript pertaining to the slow kinetics hypothesis</td></tr><tr><td>·fast</td><td>Subscript pertaining to the fast kinetics hypothesis</td></tr><tr><td>β</td><td>Derivative of the saturated water vapor concentration with respect to temperature, β = dc<sub>sat</sub>/dT</td></tr><tr><td>L</td><td>Latent heat of sublimation of ice</td></tr><tr><td>k<sub>i</sub></td><td>Thermal conductivity of ice</td></tr><tr><td>k<sub>a</sub></td><td>Thermal conductivity of air</td></tr><tr><td>k<sub>v</sub></td><td>Apparent thermal conductivity of air, k<sub>v</sub> = k<sub>a</sub> + βLD<sub>0</sub></td></tr></tbody></table>
*Code availability.* The codes for the numerical simulations and their analysis will be provided upon direct request to the corresponding author.

*Data availability.* The computed values of effective thermal conductivity of snow and of effective diffusion coefficient of water vapor in snow are available in the Supplement of the article.

*Supplement.* The supplement related to this article is available online at: https://doi.org/10.5194/tc-15-2739-2021-supplement.

*Author contributions.* The research was designed by FD, KF, and PH. FD obtained funding. KF performed research and wrote the paper with inputs from FD and PH.

*Competing interests.* The authors declare that they have no conflict of interest.

*Acknowledgements.* We acknowledge Marion Reveillet, Marie Dumont, François Tuzet, Neige Calonne, Anne Dufour, and the ANR JCJC EBONI (grant no. ANR-16-CE01-006) for providing the tomography scan of a melt forms sample. The rest of the unpublished samples were obtained thanks to the help of Jacques Roulle and the tomography apparatus with funding from the INSU-LEFE, the LabEx OSUG, and the CNRM. We thank Neige Calonne and Marie Dumont for their valuable inputs for the article. We are thankful to Charles Fierz and the two anonymous reviewers for reviewing the manuscript and to Carrie Vuyovich for editing it.

*Financial support.* This research has been supported by the Climate Initiative program of the BNP-Paribas Fondation (Acceleration of Permafrost Thaw program grant).

*Review statement.* This paper was edited by Carrie Vuyovich and reviewed by Charles Fierz and two anonymous referees.

## References

Auriault, J.: Heterogeneous medium. Is an equivalent macroscopic description possible?, International J. Engin. Sci., 29, 785–795, https://doi.org/10.1016/0020-7225(91)90001-J, 1991.

Auriault, J.-L., Boutin, C., and Geindreau, C.: Homogenization of coupled phenomena in heterogenous media, vol. 149, John Wiley & Sons, 2010.

Batchelor, G. K. and Brien, R. W.: Thermal or electrical conduction through a granular material, Proc. Royal Soc. Lond. A. Math. Phys. Sci., 355, 313–333, https://doi.org/10.1098/rspa.1977.0100, 1977.

Calonne, N., Flin, F., Morin, S., Lesaffre, B., du Roscoat, S. R., and Geindreau, C.: Numerical and experimental investigations of the

effective thermal conductivity of snow, Geophys. Res. Lett., 38, L23501, https://doi.org/10.1029/2011GL049234, 2011.

Calonne, N., Geindreau, C., and Flin, F.: Macroscopic modeling for heat and water vapor transfer in dry snow by homogenization, J. Phys. Chem. B, 118, 13393–13403, https://doi.org/10.1021/jp5052535, 2014.

Calonne, N., Milliancourt, L., Burr, A., Philip, A., Martin, C. L., Flin, F., and Geindreau, C.: Thermal Conductivity of Snow, Firn, and Porous Ice From 3-D Image-Based Computations, Geophys. Res. Let., 46, 13079–13089, https://doi.org/10.1029/2019GL085228, 2019.

Colbeck, S. C.: An overview of seasonal snow metamorphism, Revi. Geophys. 20, 45–61, https://doi.org/10.1029/RG020i001p00045, 1982.

Colbeck, S. C.: Theory of metamorphism of dry snow, J. Geophys. Res.-Oceans, 88, 5475–5482, https://doi.org/10.1029/JC088iC09p05475, 1983.

Colbeck, S. C.: The vapor diffusion coefficient for snow, Water Resour. Res., 29, 109–115, https://doi.org/10.1029/92WR02301, 1993.

De Vries, D. A.: Simultaneous transfer of heat and moisture in porous media, Eos, Trans. Am. Geophys. Union, 39, 909–916, https://doi.org/10.1029/TR039i005p00909, 1958.

De Vries, D. A.: The theory of heat and moisture transfer in porous media revisited, Int. J. Heat Mass Transf., 30, 1343–1350, https://doi.org/10.1016/0017-9310(87)90166-9, 1987.

Domine, F., Barrere, M., Sarrazin, D., Morin, S., and Arnaud, L.: Automatic monitoring of the effective thermal conductivity of snow in a low-Arctic shrub tundra, The Cryosphere, 9, 1265–1276, https://doi.org/10.5194/tc-9-1265-2015, 2015.

Domine, F., Picard, G., Morin, S., Barrere, M., Madore, J.-B., and Langlois, A.: Major Issues in Simulating Some Arctic Snow-pack Properties Using Current Detailed Snow Physics Models: Consequences for the Thermal Regime and Water Budget of Permafrost, J. Adv. Model. Earth Syst., 11, 34–44, 2019.

Fierz, C., Armstrong, R. L., Durand, Y., Etchevers, P., Greene, E., McClung, D. M., Nishimura, K., Satyawali, P. K., and Sokratov, S. A.: The International Classification for Seasonal Snow on the Ground, The International Classification for Seasonal Snow on the Ground, IHP-VII Technical Documents in Hydrology No. 83, IACS Contribution No. 1, UNESCO-IHP, Paris, 2009.

Fourteau, K., Domine, F., and Hagenmuller, P.: Macroscopic water vapor diffusion is not enhanced in snow, The Cryosphere, 15, 389–406, https://doi.org/10.5194/tc-15-389-2021, 2021.

Giddings, J. C. and LaChapelle, E.: The formation rate of depth hoar, J. Geophys. Res., 67, 2377–2383, https://doi.org/10.1029/JZ067i006p02377, 1962.

Gilbert, A., Vincent, C., Wagnon, P., Thibert, E., and Rabatel, A.: The influence of snow cover thickness on the thermal regime of Tête Rousse Glacier (Mont Blanc range, 3200 m asl): Consequences for outburst flood hazards and glacier response to climate change, J. Geophys. Res.-Ea. Surf., 117, F04018, https://doi.org/10.1029/2011JF002258, 2012.

Hagenmuller, P., Matzl, M., Chambon, G., and Schneebeli, M.: Sensitivity of snow density and specific surface area measured by microtomography to different image processing algorithms, The Cryosphere, 10, 1039–1054, https://doi.org/10.5194/tc-10-1039-2016, 2016.
Hagenmuller, P., Flin, F., Dumont, M., Tuzet, F., Peinke, I., Lal-
palus, P., Dufour, A., Roulle, J., Pézard, L., Voisin, D., Ando, E.,
Rolland du Roscoat, S., and Charrier, P.: Motion of dust parti-
cles in dry snow under temperature gradient metamorphism, The
Cryosphere, 13, 2345–2359, https://doi.org/10.5194/tc-13-2345-
2019, 2019.

Hansen, A. C. and Foslien, W. E.: A macroscale mixture the-
ory analysis of deposition and sublimation rates during heat
and mass transfer in dry snow, The Cryosphere, 9, 1857–1878,
https://doi.org/10.5194/tc-9-1857-2015, 2015.

Jaafar, H. and Picot, J. J. C.: Thermal conductivity of snow by a transient state probe method, Water Resour. Res., 6, 333–335, https://doi.org/10.1029/WR006i001p00333, 1970.

Jordan, R.: A one-dimensional temperature model for a snow cover: Technical documentation for SNTHERM. 89., Tech. rep., Cold Regions Research and Engineering Lab Hanover NH, 1991.

Kadoya, K., Matsunaga, N., and Nagashima, A.: Viscosity and Thermal Conductivity of Dry Air in the Gaseous Phase, J. Phys. Chem. Ref. Data, 14, 947–970, https://doi.org/10.1063/1.555744, 1985.

Kaempfer, T. U. and Plapp, M.: Phase-field modeling of dry snow metamorphism, Phys. Rev. E, 79, 031502, https://doi.org/10.1103/PhysRevE.79.031502, 2009.

Krol, Q. and Löwe, H.: Analysis of local ice crystal growth in snow, J. Glaciol., 62, 378–390, https://doi.org/10.1017/jog.2016.32, 2016.

Lecomte, O., Fichefet, T., Vancoppenolle, M., Domine, F., Mas-
sonnet, F., Mathiot, P., Morin, S., and Barriat, P.-Y.: On
the formulation of snow thermal conductivity in large-scale
sea ice models, J. Adv. Model. Earth Syst., 5, 542–557,
https://doi.org/10.1002/jame.20039, 2013.

Legagneux, L. and Domine, F.: A mean field model of the de-
crease of the specific surface area of dry snow during isother-
mal metamorphism, J. Geophys. Res. Earth Surf., 110, F04011,
https://doi.org/10.1029/2004JF000181, 2005.

Libbrecht, K. G.: Precision Measurements of Ice Crystal Growth Rates, Tech. rep., Department of Physics, California Institute of Technology, Pasadena, California 91125, US, 2006.

Libbrecht, K. G. and Rickerby, M. E.: Measure-
ments of surface attachment kinetics for faceted
ice crystal growth, J. Crystal Growth, 377, 1–8,
https://doi.org/10.1016/j.jcrysgro.2013.04.037, 2013.

Lide, D. R.: CRC handbook of chemistry and physics, chap. Prop-
erties of ice and supercooled water, pp. 6–5, CRC press, Taylor
and Francis, Boca Raton, FL, 85 edn., 2006.

Löwe, H., Riche, F., and Schneebeli, M.: A general treatment of snow microstructure exemplified by an improved relation for thermal conductivity, The Cryosphere, 7, 1473–1480, https://doi.org/10.5194/tc-7-1473-2013, 2013.

Malinen, M. and Råback, P.: Elmer Finite Element Solver for Multiphysics and Multiscale Problems, in: Multiscale Modelling Methods for Applications in Materials Science, edited by Kovdov, I. and Sutmann, G., Forschungszentrum Jülich GmbH, pp. 101–113, 2013.

Morin, S., Domine, F., Arnaud, L., and Picard, G.: In-situ monitoring of the time evolution of the effective thermal conductivity of snow, Cold Reg. Sci. Tech., 64, 73–80, https://doi.org/10.1016/j.coldregions.2010.02.008, 2010.

Moyne, C., Batsale, J.-C., and Degiovanni, A.: Approche ex-
périmentale et théorique de la conductivité thermique des mi-
lieux poreux humides – II. Théorie, Int. J. Heat Mass Transf.,
31, 2319–2330, https://doi.org/10.1016/0017-9310(88)90163-9,
1988.

Municchi, F. and Icardi, M.: Macroscopic models for filtration and heterogeneous reactions in porous media, Adv. Water Resour., 141, 103605, https://doi.org/10.1016/j.adwwatres.2020.103605, 2020.

Peinke, I., Hagenmuller, P., Andò, E., Chambon, G., Flin, F., and Roulle, J.: Experimental Study of Cone Penetration in Snow Using X-Ray Tomography, Front. Earth Sci., 8, 63, https://doi.org/10.3389/feart.2020.00063, 2020.

Pinzer, B. R., Schneebeli, M., and Kaempfer, T. U.: Va-
por flux and recrystallization during dry snow metamor-
phism under a steady temperature gradient as observed by
time-lapse micro-tomography, The Cryosphere, 6, 1141–1155,
https://doi.org/10.5194/tc-6-1141-2012, 2012.

Riche, F. and Schneebeli, M.: Thermal conductivity of snow mea-
sured by three independent methods and anisotropy considera-
tions, The Cryosphere, 7, 217–227, https://doi.org/10.5194/tc-7-
217-2013, 2013.

Saito, Y.: Statistical physics of crystal growth, World Scientific,
1996.

Shertzer, R. H. and Adams, E. E.: A Mass Diffusion Model for Dry Snow Utilizing a Fabric Tensor to Characterize Anisotropy, J Adv. Model. Earth Sys., 10, 881–890, https://doi.org/10.1002/2017MS001046, 2018.

Slack, G. A.: Thermal conductivity of ice, Phys. Rev. B, 22, 3065–3071, https://doi.org/10.1103/PhysRevB.22.3065, 1980.

Sokratov, S. A. and Maeno, N.: Effective water vapor diffusion co-efficient of snow under a temperature gradient, Water Resour. Res., 36, 1269–1276, https://doi.org/10.1029/2000WR900014, 2000.

Sommerfeld, R. A. and LaChapelle, E.: The Classification of Snow Metamorphism, J. Glaciol., 9, 3–18, https://doi.org/10.3189/S0022143000026757, 1970.

Sturm, M. and Benson, C. S.: Vapor transport, grain growth and depth-hoar development in the subarctic snow, J. Glaciol., 43, 42–59, https://doi.org/10.3189/S0022143000002793, 1997.

Sturm, M. and Johnson, J. B.: Natural convection in the subar-
tic snow cover, J. Geophys. Res.-Sol. Ea., 96, 11657–11671,
https://doi.org/10.1029/91JB00895, 1991.

Sturm, M. and Johnson, J. B.: Thermal conductivity measure-
ments of depth hoar, J. Geophys. Res.-Sol. Ea., 97, 2129–2139,
https://doi.org/10.1029/91JB02685, 1992.

Sturm, M., Holmgren, J., König, M., and Morris, K.: The ther-
mal conductivity of seasonal snow, J. Glaciol., 43, 26–41,
https://doi.org/10.3189/S0022143000002781, 1997.

Trabant, D. and Benson, C.: Field experiments on the development of depth hoar, Geol. Soc. Am. Mem., 135, 309–322, 1972.

Vionnet, V., Brun, E., Morin, S., Boone, A., Faroux, S., Le Moigne, P., Martin, E., and Willemet, J.-M.: The detailed snow-pack scheme Crocus and its implementation in SURFEX v7.2, Geosci. Model Dev., 5, 773–791, https://doi.org/10.5194/gmd-5-773-2012, 2012.

Whitaker, S.: Simultaneous Heat, Mass, and Momentum Transfer in Porous Media: A Theory of Drying, Adv. Heat Transf., 13, 119–203, https://doi.org/10.1016/S0065-2717(08)70223-5, 1977.
Yosida, Z., Oura, H., Kuroiwa, D., Huzioka, T., Kojima, K., Aoki, S. I., and Kinosita, S.: Physical Studies on Deposited Snow. I. Thermal Properties, Contributions from the Institute of Low Temperature Science, Hokkaido, Japan, 7, 19–74, available at: http://hdl.handle.net/2115/20216 (last access: 15 June 2021), 1955.

Zhang, T., Osterkamp, T. E., and Stamnes, K.: Influence of the depth hoar layer of the seasonal snow cover on the ground thermal regime, Water Resour. Res., 32, 2075–2086, https://doi.org/10.1029/96WR00996, 1996.