scholarly journals Liquid water infiltration into a layered snowpack: evaluation of a 3D water transport model with laboratory experiments

2017 ◽  
Author(s):  
Hiroyuki Hirashima ◽  
Francesco Avanzi ◽  
Satoru Yamaguchi

Abstract. The heterogeneous movement of liquid water through snowpack during precipitation and snowmelt leads to complex liquid water distributions that are important for avalanche and runoff forecasting. We reproduced the formation of capillary barriers and the development of preferential flow through snow using a multi-dimensional water transport model, which was then validated using laboratory experiments of liquid water infiltration into layered, initially dry snow. Three-dimensional simulations assumed the same column shape and size, grain size, snow density, and water input rate as the laboratory experiments. Model evaluation focused on the timing of water movement, the thickness of the upper layer affected by ponding, and on water content profiles and the wet snow fraction. Simulation results showed that the model reconstructs some relevant features of capillary barriers including ponding in the upper layer, preferential infiltration far from the interface, and the timing of liquid water arrival at the snow base. In contrast, the area of preferential flow paths was usually underestimated and consequently the averaged water content in areas characterized by preferential flow paths was also underestimated. Improving the representation of water preferential infiltration into initially dry snow is necessary to reproduce the transition from a dry-snow-dominant condition to a wet-snow-dominant one, especially in long-period simulations.

2017 ◽  
Vol 21 (11) ◽  
pp. 5503-5515 ◽  
Author(s):  
Hiroyuki Hirashima ◽  
Francesco Avanzi ◽  
Satoru Yamaguchi

Abstract. The heterogeneous movement of liquid water through the snowpack during precipitation and snowmelt leads to complex liquid water distributions that are important for avalanche and runoff forecasting. We reproduced the formation of capillary barriers and the development of preferential flow through snow using a three-dimensional water transport model, which was then validated using laboratory experiments of liquid water infiltration into layered, initially dry snow. Three-dimensional simulations assumed the same column shape and size, grain size, snow density, and water input rate as the laboratory experiments. Model evaluation focused on the timing of water movement, thickness of the upper layer affected by ponding, water content profiles and wet snow fraction. Simulation results showed that the model reconstructs relevant features of capillary barriers, including ponding in the upper layer, preferential infiltration far from the interface, and the timing of liquid water arrival at the snow base. In contrast, the area of preferential flow paths was usually underestimated and consequently the averaged water content in areas characterized by preferential flow paths was also underestimated. Improving the representation of preferential infiltration into initially dry snow is necessary to reproduce the transition from a dry-snow-dominant condition to a wet-snow-dominant one, especially in long-period simulations.


2016 ◽  
Vol 10 (5) ◽  
pp. 2013-2026 ◽  
Author(s):  
Francesco Avanzi ◽  
Hiroyuki Hirashima ◽  
Satoru Yamaguchi ◽  
Takafumi Katsushima ◽  
Carlo De Michele

Abstract. Data of liquid water flow around a capillary barrier in snow are still limited. To gain insight into this process, we carried out observations of dyed water infiltration in layered snow at 0 °C during cold laboratory experiments. We considered three different finer-over-coarser textures and three different water input rates. By means of visual inspection, horizontal sectioning, and measurements of liquid water content (LWC), capillary barriers and associated preferential flow were characterized. The flow dynamics of each sample were also simulated solving the Richards equation within the 1-D multi-layer physically based snow cover model SNOWPACK. Results revealed that capillary barriers and preferential flow are relevant processes ruling the speed of water infiltration in stratified snow. Both are marked by a high degree of spatial variability at centimeter scale and complex 3-D patterns. During unsteady percolation of water, observed peaks in bulk volumetric LWC at the interface reached  ∼ 33–36 vol % when the upper layer was composed by fine snow (grain size smaller than 0.5 mm). However, LWC might locally be greater due to the observed heterogeneity in the process. Spatial variability in water transmission increases with grain size, whereas we did not observe a systematic dependency on water input rate for samples containing fine snow. The comparison between observed and simulated LWC profiles revealed that the implementation of the Richards equation reproduces the existence of a capillary barrier for all observed cases and yields a good agreement with observed peaks in LWC at the interface between layers.


2016 ◽  
Author(s):  
Nicolas R. Leroux ◽  
John W. Pomeroy

Abstract. Accurate estimation of the water flux through melting snowpacks is of primary importance for runoff prediction. Lateral flows and preferential flow pathways in porous media flow have proven critical for improving soil and groundwater flow models, but though many physically-based layered snowmelt models have been developed, only 1D matrix flow over level ground is currently accounted for in snow models. Snowmelt models that include these processes may improve snowmelt discharge timing and contributing area calculations in hydrological models. A two-dimensional snow model (SMPP – Snowmelt Model with Preferential flow Paths) is presented that simulates heat and water flows through both snowpack matrix and preferential flow paths, as well as snowmelt and refreezing of meltwater. The model assumes thermodynamic equilibrium between solid and liquid phases and uses the latest improvements made in snow science to estimate snow hydraulic and thermal properties. A finite volume method is applied to solve for the 2D heat and water equations. The use of a water entry pressure for dry snow combined with consideration of the impact of heterogeneities in surface fluxes and internal snow properties – density, grain size and layer thickness – allowed calculation of the formation of preferential flow paths in the snowpack. The simulation of water flow through preferential flow paths resulted in liquid water reaching the base of the snowpack earlier than for a homogeneous wetting front. Moreover, the preferential flow paths in the model increased the exchange of energy between the snow surface and the internal snowpack, resulting in faster warming of the snowpack. A sensitivity analysis, conducted on the snow internal properties showed that initial conditions such as density and temperature, should be carefully measured in the field to accurately estimate liquid water percolating through the snowpack. Furthermore, two empirical coefficients used in the water flow equation were showed to greatly impact model outputs. This heterogeneous flow model is an important tool to help understand snowmelt flow processes in complex and level terrains and to demonstrate how uncertainty in snowmelt-derived runoff calculations might be reduced.


2016 ◽  
Vol 10 (6) ◽  
pp. 2731-2744 ◽  
Author(s):  
Nander Wever ◽  
Sebastian Würzer ◽  
Charles Fierz ◽  
Michael Lehning

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.


2011 ◽  
Vol 52 (58) ◽  
pp. 201-208 ◽  
Author(s):  
Christoph Mitterer ◽  
Hiroyuki Hirashima ◽  
Jürg Schweizer

AbstarctWet-snow avalanches are difficult to forecast, as the change from stable to unstable snow conditions occurs rapidly in a wet snowpack, often in response to water production and movement. Snow stratigraphy plays a vital role in determining flux behaviour. Capillary barriers or melt–freeze crusts can impede and divert water horizontally over large areas and thus may act as a failure layer for wet-snow avalanches. We present a comparison of measured and modelled liquid water content, θw, and snow stratigraphy during periods of wet-snow instabilities. Special attention is given to the reproducibility of capillary barriers, ponding of water on melt–freeze crusts and the timing of first wetting and of water arrival at the bottom of the snowpack, because these factors are believed to play a major role in the formation of wet-snow avalanches. In situ measurements were performed in the vicinity of automatic weather stations or close to recent wet-snow avalanches in order to compare them with model results. The simulations are based on two different water flux models incorporated within the 1-D snow-cover model SNOWPACK. The comparison of the two model runs with observed θw and stratigraphy revealed that both water-transport models reproduced the ponding of water on melt–freeze crusts. However, in both models melt–freeze crusts were transformed to normal melt forms earlier than observed in nature, so still existing ponding was not captured by the models. Only one of the models was able to reproduce capillary barriers in agreement with observations. The time of the first wetting at the surface was well predicted, but the simulated arrival time of the wetting front at the bottom of the snowpack differed between the simulations; it was either too early or too late compared with the observation.


Sign in / Sign up

Export Citation Format

Share Document