Abstract. The runoff from a snow cover during spring snowmelt or rain-on-snow events is an important factor in the hydrological cycle. In this study, three water balance schemes for the 1 dimensional physically-based snowpack model SNOWPACK are compared to lysimeter measurements at two alpine sites with a seasonal snow cover, but with different climatological conditions: Weissfluhjoch (WFJ) and Col de Porte (CDP). The studied period consists of 14 and 17 yr, respectively. The schemes include a simple buckettype approach, an approximation of Richards Equation (RE), and the full RE. The results show that daily sums of snowpack runoff are strongly related to a positive energy balance of the snow cover and therefore, all water balance schemes show very similar performance in terms of Nash-Sutcliffe efficiency (NSE) coefficients (around 0.63 and 0.72 for WFJ and CDP, respectively) and r 2 values (around 0.83 and 0.72 for WFJ and CDP, respectively). An analysis of the runoff dynamics over the season showed that the bucket-type and approximated RE scheme release meltwater slower than in the measurements, whereas RE provides a better agreement. Overall, solving RE for the snow cover yields the best agreement between modelled and measured snowpack runoff, but differences between the schemes are small. On sub-daily time scales, the water balance schemes behave very differently. In that case, solving RE provides the highest agreement between modelled and measured snowpack runoff in terms of NSE coefficient (around 0.48 at both sites). At WFJ, the other water balance schemes loose most predictive power, whereas at CDP, the bucket-type scheme has an NSE coefficient of 0.39. The shallower and less stratified snowpack at CDP likely reduces the differences between the water balance schemes. Accordingly, it can be concluded that solving RE for the snow cover improves several aspects of modelling snow cover runoff, especially for deep, sub-freezing snow covers and in particular on the sub-daily time scales. The additional computational cost was found to be in the order of a factor of 1.5-2.

No abstract

No abstract

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 SNOW-PACK. 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.

No abstract