Modelling snowpack bulk density using snow depth,cumulative degree-days,and climatological predictor variables

Snowpack water equivalent (SWE) is a key variable for water resource management in snow-dominated catchments. While it is not feasible to quantify SWE at the catchment scale using either field surveys or remotely sensed data, technologies such as airborne LiDAR (light detection and ranging) support the mapping of snow depth at scales relevant to operational water management. To convert snow depth to water equivalent, models have been developed to predict SWE or snowpack density based on snow depth and additional predictor variables. This study builds upon previous models that relate snowpack density to snow depth by including additional predictor variables to account for (1) long-term climatologies that describe the prevailing conditions influencing regional snowpack properties, and (2) the effect of intra- and inter-year variability in meteorological conditions on densification through a cumulative degree-day index derived from North American Regional Reanalysis products. A non-linear model was fit to 114 506 snow survey measurements spanning 41 years from 1166 snow courses across western North America. Under spatial cross-validation, the predicted densities had a root-mean-square error of 47.1 kg m⁻³, a mean bias of −0.039 kg m⁻³, and a Nash-Sutcliffe Efficiency of 0.70. The model developed in this study had similar overall performance compared to a similar regression-based model reported in the literature, but had reduced seasonal biases. When applied to predict SWE from simulated depths with random errors consistent with those obtained from LiDAR or Structure-from-Motion, 50% of the SWE estimates for April and May fell within −45 to 49 mm of the observed SWE, representing prediction errors of −15% to 20%.