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

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.

[1]  M. Lehning,et al.  Use of a multilayer snow model to assess grazing conditions for reindeer , 2013, Annals of Glaciology.

[2]  James P. McNamara,et al.  An evaluation of the hydrologic relevance of lateral flow in snow at hillslope and catchment scales , 2013 .

[3]  O. Eisen,et al.  Verification of the multi-layer SNOWPACK model with different water transport schemes , 2015 .

[4]  E. Brun,et al.  A numerical model to simulate snow-cover stratigraphy for operational avalanche forecasting , 1992, Journal of Glaciology.

[5]  A. Sato,et al.  Experimental analysis of preferential flow in dry snowpack , 2013 .

[6]  Philip Marsh,et al.  Meltwater Movement in Natural Heterogeneous Snow Covers , 1985 .

[7]  M. Lehning,et al.  Modelling liquid water transport in snow under rain-on-snow conditions – considering preferential flow , 2016 .

[8]  S. Colbeck Water flow through heterogeneous snow , 1979 .

[9]  S. Yamaguchi,et al.  A Multi-Dimensional Water Transport Model to Reproduce the Preferential Flow in a Snowpack , 2014 .

[10]  Mark F. Meier,et al.  Retention of Greenland runoff by refreezing: Implications for projected future sea level change , 1991 .

[11]  Hiroyuki Hirashima,et al.  Observations of capillary barriers and preferential flow in layered snow during cold laboratory experiments , 2016 .

[12]  R. Kattelmann Snowmelt lysimeters in the evaluation of snowmelt models , 2000, Annals of Glaciology.

[13]  Chris Derksen,et al.  Observed and modelled effects of ice lens formation on passive microwave brightness temperatures over snow covered tundra , 2010 .

[14]  M. Lehning,et al.  Influence of Initial Snowpack Properties on Runoff Formation during Rain-on-Snow Events , 2016 .

[15]  M. Lehning,et al.  Solving Richards Equation for snow improves snowpack meltwater runoff estimations in detailed multi-layer snowpack model , 2014 .

[16]  A. Sato,et al.  Dependence of the water retention curve of snow on snow characteristics , 2012, Annals of Glaciology.

[17]  I. Baker,et al.  Investigating the thermophysical properties of the ice–snow interface under a controlled temperature gradient: Part I: Experiments & Observations , 2015 .

[18]  E. Mosley‐Thompson,et al.  Greenland meltwater storage in firn limited by near-surface ice formation , 2016 .

[19]  Michael Lehning,et al.  snowpack model calculations for avalanche warning based upon a new network of weather and snow stations , 1999 .

[20]  W. T. Pfeffer,et al.  Ice Columns and Frozen Rills in a Warm Snowpack, Green Lakes Valley, Colorado, U.S.A. , 2000 .

[21]  P. Bartelt,et al.  A physical SNOWPACK model for the Swiss avalanche warning Part III: meteorological forcing, thin layer formation and evaluation , 2002 .

[22]  Rolf Weingartner,et al.  Retrospective analysis of a nonforecasted rain-on-snow flood in the Alps – a matter of model limitations or unpredictable nature? , 2013 .

[23]  A. Roy,et al.  Microwave snow emission modeling uncertainties in boreal and subarctic environments , 2015 .

[24]  D. DiCarlo Stability of gravity‐driven multiphase flow in porous media: 40 Years of advancements , 2013 .

[25]  Y. Takeuchi,et al.  A multiple snow layer model including a parameterization of vertical water channel process in snowpack , 2008 .

[26]  Martin Schneebeli,et al.  A constant-speed penetrometer for high-resolution snow stratigraphy , 1998, Annals of Glaciology.

[27]  Tammo S. Steenhuis,et al.  Wetting front instability: 1. Theoretical discussion and dimensional analysis , 1989 .

[28]  Tammo S. Steenhuis,et al.  Wetting front instability: 2. Experimental determination of relationships between system parameters and two‐dimensional unstable flow field behavior in initially dry porous media , 1989 .

[29]  Daqing Yang,et al.  Winter rain on snow and its association with air temperature in northern Eurasia , 2008 .

[30]  W. T. Pfeffer,et al.  Thermal tracking of meltwater retention in Greenland's accumulation area , 2012 .

[31]  Philip Marsh,et al.  Wetting front advance and freezing of meltwater within a snow cover: 2. A simulation model , 1984 .

[32]  F. Perron,et al.  Radar investigations of melt pathways in a natural snowpack , 1999 .

[33]  E. Mosley‐Thompson,et al.  Changes in the firn structure of the western Greenland Ice Sheet caused by recent warming , 2015 .

[34]  M. Phillips,et al.  Seasonally intermittent water flow through deep fractures in an Alpine Rock Ridge: Gemsstock, Central Swiss Alps , 2016 .

[35]  M. Zambrano-Bigiarini Goodness-of-fit functions for comparison of simulated andobserved hydrological time series , 2014 .

[36]  P. Bartelt,et al.  A physical SNOWPACK model for the Swiss avalanche warning: Part II. Snow microstructure , 2002 .

[37]  Desiree Tullos,et al.  Variability in effect of climate change on rain-on-snow peak flow events in a temperate climate , 2013 .

[38]  M. Schneebeli Development and stability of preferential flow paths in a layered snowpack , 1995 .

[39]  N. Wever,et al.  Assessing wet snow avalanche activity using detailed physics based snowpack simulations , 2016 .

[40]  P. Bartelt,et al.  A physical SNOWPACK model for the Swiss avalanche warning: Part I: numerical model , 2002 .

[41]  P. Marsh Snowcover formation and melt: recent advances and future prospects , 1999 .

[42]  C. Derksen,et al.  Brief communication: Improved measurement of ice layer density in seasonal snowpacks , 2016 .

[43]  R. Kattelmann,et al.  Macropores in Snowpacks of Sierra Nevada , 1985, Annals of Glaciology.

[44]  C. Willmott ON THE VALIDATION OF MODELS , 1981 .

[45]  P. Marsh,et al.  Flow-finger continuity in serial thick-sections in a melting Sierran snowpack , 1995 .

[46]  M. Lehning,et al.  Model simulations of the modulating effect of the snow cover in a rain-on-snow event , 2014 .

[47]  Mark F. Meier,et al.  Modeling of Meltwater Infiltration in Subfreezing Snow , 1990 .

[48]  K. Tonnessen,et al.  Biogeochemistry of seasonally snow-covered catchments , 1995 .

[49]  Christian Geindreau,et al.  3-D image-based numerical computations of snow permeability: links to specific surface area, density, and microstructural anisotropy , 2012 .

[50]  S. Colbeck,et al.  The layered character of snow covers , 1991 .

[51]  Pratap Singh,et al.  Hydrological response of snowpack under rain-on-snow events: a field study , 1997 .