Are you an EPFL student looking for a semester project?
Work with us on data science and visualisation projects, and deploy your project as an app on top of Graph Search.
Many methods exist to model snow densification in order to calculate the depth of a single snow layer or the depth of the total snow cover from its mass. Most of these densification models need to be tightly integrated with an accumulation and melt model and need many forcing variables at high temporal resolution. However, when trying to model snow depth (HS) on climatological timescales, which is often needed for winter tourism-related applications, these preconditions can cause barriers. Often, for these types of applications, empirical snow models are used. These can estimate snow accumulation and snowmelt based on daily precipitation and temperature data only. To convert the resultant snow water equivalent (SWE) time series into snow depth, we developed the empirical model SWE2HS. SWE2HS is a multilayer densification model which uses daily snow water equivalent as sole input. A constant new snow density is assumed and densification is calculated via exponential settling functions. The maximum snow density of a single layer changes over time due to overburden and SWE losses. SWE2HS has been calibrated on a data set derived from a network of manual snow stations in Switzerland. It has been validated against independent data derived from automatic weather stations (AWSs) in the European Alps (Austria, France, Germany, Switzerland) and against withheld data from the Swiss manual observer station data set which was not used for calibration. The model fits the calibration data with root mean squared error (RMSE) of 8.4 cm, coefficient of determination (R2) of 0.97, and bias of -0.3 cm; it is able to achieve RMSE of 20.5 cm, R2 of 0.92, and bias of 2.5 cm on the validation data set from automatic weather stations and RMSE of 7.9 cm, R2 of 0.97, and bias of -0.3 cm on the validation data set from manual stations. The temporal evolution of the bulk density can be reproduced reasonably well on all three data sets. Due to its simplicity, the model can be used as post-processing tool for output of any other snow model that provides daily snow water equivalent output. Owing to its empirical nature, SWE2HS should only be used in regions with a similar snow climatology as the European Alps or has to be recalibrated for other snow climatological conditions. The SWE2HS model is available as a Python package which can be easily installed and used.
Michael Lehning, Dylan Stewart Reynolds, Michael Haugeneder