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

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, multi-layer 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. We found that preferential flow paths arriving at a layer transition in the snowpack may lead to ponding conditions. Subsequent refreezing then can form dense layers in the snowpack, that regularly exceed 700 kg m−3. A comparison of simulated density profiles with bi-weekly 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 multi-dimensional model.

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.


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.


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.


2019 ◽  
Vol 13 (7) ◽  
pp. 1819-1842 ◽  
Author(s):  
Vincent Verjans ◽  
Amber A. Leeson ◽  
C. Max Stevens ◽  
Michael MacFerrin ◽  
Brice Noël ◽  
...  

Abstract. As surface melt is increasing on the Greenland Ice Sheet (GrIS), quantifying the retention capacity of the firn layer is critical to linking meltwater production to meltwater runoff. Firn-densification models have so far relied on empirical approaches to account for the percolation–refreezing process, and more physically based representations of liquid water flow might bring improvements to model performance. Here we implement three types of water percolation schemes into the Community Firn Model: the bucket approach, the Richards equation in a single domain and the Richards equation in a dual domain, which accounts for partitioning between matrix and fast preferential flow. We investigate their impact on firn densification at four locations on the GrIS and compare model results with observations. We find that for all of the flow schemes, significant discrepancies remain with respect to observed firn density, particularly the density variability in depth, and that inter-model differences are large (porosity of the upper 15 m firn varies by up to 47 %). The simple bucket scheme is as efficient in replicating observed density profiles as the single-domain Richards equation, and the most physically detailed dual-domain scheme does not necessarily reach best agreement with observed data. However, we find that the implementation of preferential flow simulates ice-layer formation more reliably and allows for deeper percolation. We also find that the firn model is more sensitive to the choice of densification scheme than to the choice of water percolation scheme. The disagreements with observations and the spread in model results demonstrate that progress towards an accurate description of water flow in firn is necessary. The numerous uncertainties about firn structure (e.g. grain size and shape, presence of ice layers) and about its hydraulic properties, as well as the one-dimensionality of firn models, render the implementation of physically based percolation schemes difficult. Additionally, the performance of firn models is still affected by the various effects affecting the densification process such as microstructural effects, wet snow metamorphism and temperature sensitivity when meltwater is present.


2020 ◽  
Vol 24 (6) ◽  
pp. 3271-3288
Author(s):  
Anne Hartmann ◽  
Ekaterina Semenova ◽  
Markus Weiler ◽  
Theresa Blume

Abstract. Preferential flow strongly controls water flow and transport in soils. It is ubiquitous but difficult to characterize and predict. This study addresses the occurrence and the evolution of preferential flow during the evolution of landscapes and here specifically during the evolution of hillslopes. We targeted a chronosequence of glacial moraines in the Swiss Alps to investigate how water flow paths evolve along with the soil-forming processes. Dye tracer irrigation experiments with a Brilliant Blue FCF solution (4 g L−1) were conducted on four moraines of different ages (30, 160, 3000, and 10 000 years). At each moraine, three dye tracer experiments were conducted on plots of 1.5 m ×1.0 m. The three plots at each moraine were characterized by different vegetation complexities (low, medium, and high). Each plot was further divided into three equal subplots for the application of three different irrigation amounts (20, 40, and 60 mm) with an average irrigation intensity of 20 mm h−1. The day after the experiment five vertical soil sections were excavated, and the stained flow paths were photographed. Digital image analysis was used to derive average infiltration depths and flow path characteristics such as the volume and surface density of the dye patterns. Based on the volume density, the observed dye patterns were assigned to specific flow type categories. The results show a significant change in the type of preferential flow paths along the chronosequence. The flow types change from a rather homogeneous matrix flow in coarse material with high conductivities and a sparse vegetation cover at the youngest moraine to a heterogeneous infiltration pattern at the medium-age moraines. Heterogeneous matrix and finger flow are dominant at these intermediate age classes. At the oldest moraine only macropore flow via root channels was observed in deeper parts of the soil, in combination with a very high water storage capacity of the organic top layer and low hydraulic conductivity of the deeper soil. In general, we found an increase in water storage with increasing age of the moraines, based on our observations of the reduction in infiltration depth as well as laboratory measurements of porosity. Preferential flow is, however, not only caused by macropores, but especially for the medium-age moraine, it seems to be mainly initiated by soil surface characteristics (vegetation patches and microtopography).


2016 ◽  
Vol 28 (2) ◽  
pp. 648-663 ◽  
Author(s):  
Yinghu Zhang ◽  
Jianzhi Niu ◽  
Mingxiang Zhang ◽  
Zixing Xiao ◽  
Weili Zhu

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.


2000 ◽  
Vol 151 (6) ◽  
pp. 198-204 ◽  
Author(s):  
Peter Waldner ◽  
Martin Schneebeli ◽  
Hans Wunderli

The evolution of the release of nutrients from a snow cover may essentially determine their availability for plants. In this study, the nitrate release from a snow pack on an open field in Alptal (canton of Schwyz, Switzerland) was measured in the winter of 1998/99. With the diminution of the nitrate concentration in single snow layers, a hint for ionic fractionation has been found. The measurements allowed to estimate the variability of the water and nitrate release from the snow pack. Its high spatial and temporal heterogeneity is explained with preferential flow paths of the melt water.


2020 ◽  
Author(s):  
Anne Hartmann ◽  
Ekaterina Semenova ◽  
Markus Weiler ◽  
Theresa Blume

Abstract. The presence or absence of preferential flow strongly controls water flow and transport in soils. It is ubiquitous, but difficult to characterize and predict. This study addresses the occurrence and the evolution of preferential flow during the evolution of landscapes, and here specifically during the evolution of hillslopes. We targeted a chronosequence of glacial moraines in the Swiss Alps to investigate how water flow paths evolve along with the soil forming processes. Dye tracer irrigation experiments with Brilliant Blue solution (4 g/l) were conducted on four moraines of different ages (30, 160, 3000, and 10 000 yrs). At each moraine, three dye tracer experiments were conducted on plots of 1.5 × 1.0 m. The three plots at each moraine were characterized by different vegetation complexities (low, medium, high). Each plot was further divided into three equal subplots for the application of three different irrigation amounts (20, 40, 60 mm) with an average irrigation intensity of 20 mm/h. The day after the experiment five vertical soil sections were excavated and the stained flow paths were photographed. Digital image analysis was used to derive average infiltration depths and flow path characteristics such as the volume and surface density of the dye patterns. Based on the volume density, the observed dye patterns were assigned to specific flow type categories. The results show a significant change in type of preferential flow paths along the chronosequence. The flow types change from a rather homogeneous gravity driven matrix flow in coarse material with high conductivities and a sparse vegetation cover at the youngest moraine to a heterogeneous infiltration pattern at the medium-age moraines. Heterogeneous matrix and finger flow are dominant at these intermediate age classes. At the oldest moraine only macro pore flow via root channels was observed in deeper parts of the soil, in combination with a very high water storage capacity of the organic top layer and low hydraulic conductivity of the deeper soil. In general, we found an increase in water storage with increasing age of the moraines, based on our observations of the reduction in infiltration depth as well as laboratory measurements of porosity. Preferential flow is, however, not only caused by macropores, but especially for the medium age moraine seems to be initiated mainly by soil surface characteristics (vegetation patches and micro-topography).


2014 ◽  
Vol 18 (1) ◽  
pp. 121-138 ◽  
Author(s):  
J. Wienhöfer ◽  
E. Zehe

Abstract. Rapid flow processes in connected preferential flow paths are widely accepted to play a key role in the rainfall–runoff response at the hillslope scale, but a quantitative description of these processes is still a major challenge in hydrological research. This paper investigates the approach of incorporating preferential flow paths explicitly in a process-based model for modelling water flow and solute transport at a steep forested hillslope. We conceptualise preferential flow paths as spatially explicit structures with high conductivity and low retention capacity, and evaluate simulations with different combinations of vertical and lateral flow paths in conjunction with variable or constant soil depths against measured discharge and tracer breakthrough. Out of 122 tested realisations, six set-ups fulfilled our selection criteria for the water flow simulation. These set-ups successfully simulated infiltration, vertical and lateral subsurface flow in structures, and allowed predicting the magnitude, dynamics and water balance of the hydrological response of the hillslope during successive periods of steady-state sprinkling on selected plots and intermittent rainfall on the entire hillslope area. The number of equifinal model set-ups was further reduced by the results of solute transport simulations. Two of the six acceptable model set-ups matched the shape of the observed breakthrough curve well, indicating that macrodispersion induced by preferential flow was captured well by the topology of the preferential flow network. The configurations of successful model set-ups suggest that preferential flow related to connected vertical and lateral flow paths is a first-order control on the hydrology of the study hillslope, whereas spatial variability of soil depth is secondary especially when lateral flow paths are present. Virtual experiments for investigating hillslope controls on subsurface processes should thus consider the effect of distinctive flow paths within the soil mantle. The explicit representation of flow paths in a hydrological process model was found to be a suitable approach for this purpose.


Sign in / Sign up

Export Citation Format

Share Document