Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 4.790 IF 4.790
  • IF 5-year value: 5.921 IF 5-year
    5.921
  • CiteScore value: 5.27 CiteScore
    5.27
  • SNIP value: 1.551 SNIP 1.551
  • IPP value: 5.08 IPP 5.08
  • SJR value: 3.016 SJR 3.016
  • Scimago H <br class='hide-on-tablet hide-on-mobile'>index value: 63 Scimago H
    index 63
  • h5-index value: 51 h5-index 51
Volume 10, issue 6
The Cryosphere, 10, 2731–2744, 2016
https://doi.org/10.5194/tc-10-2731-2016
© Author(s) 2016. This work is distributed under
the Creative Commons Attribution 3.0 License.
The Cryosphere, 10, 2731–2744, 2016
https://doi.org/10.5194/tc-10-2731-2016
© Author(s) 2016. This work is distributed under
the Creative Commons Attribution 3.0 License.

Research article 16 Nov 2016

Research article | 16 Nov 2016

Simulating ice layer formation under the presence of preferential flow in layered snowpacks

Nander Wever1,2, Sebastian Würzer1,2, Charles Fierz2, and Michael Lehning2,1 Nander Wever et al.
  • 1École Polytechnique Fédérale de Lausanne (EPFL), School of Architecture, Civil and Environmental Engineering, Lausanne, Switzerland
  • 2WSL Institute for Snow and Avalanche Research SLF, Davos, Switzerland

Abstract. For physics-based snow cover models, simulating the formation of dense ice layers inside the snowpack has been a long-time challenge. Their formation is considered to be tightly coupled to the presence of preferential flow, which is assumed to happen through flow fingering. Recent laboratory experiments and modelling techniques of liquid water flow in snow have advanced the understanding of conditions under which preferential flow paths or flow fingers form. We propose a modelling approach in the one-dimensional, multilayer snow cover model SNOWPACK for preferential flow that is based on a dual domain approach. The pore space is divided into a part that represents matrix flow and a part that represents preferential flow. Richards' equation is then solved for both domains and only water in matrix flow is subjected to phase changes. We found that preferential flow paths arriving at a layer transition in the snowpack may lead to ponding conditions, which we used to trigger a water flow from the preferential flow domain to the matrix domain. Subsequent refreezing then can form dense layers in the snowpack that regularly exceed 700 kg m−3. A comparison of simulated density profiles with biweekly snow profiles made at the Weissfluhjoch measurement site at 2536 m altitude in the Eastern Swiss Alps for 16 snow seasons showed that several ice layers that were observed in the field could be reproduced. However, many profiles remain challenging to simulate. The prediction of the early snowpack runoff also improved under the consideration of preferential flow. Our study suggests that a dual domain approach is able to describe the net effect of preferential flow on ice layer formation and liquid water flow in snow in one-dimensional, detailed, physics-based snowpack models, without the need for a full multidimensional model.

Publications Copernicus
Download
Short summary
The study presents a dual domain approach to simulate liquid water flow in snow using the 1-D physics based snow cover model SNOWPACK. In this approach, the pore space is separated into a part for matrix flow and a part that represents preferential flow. Using this approach, water can percolate sub-freezing snow and form dense (ice) layers. A comparison with snow pits shows that some of the observed ice layers were reproduced by the model while others remain challenging to simulate.
The study presents a dual domain approach to simulate liquid water flow in snow using the 1-D...
Citation