State of the art numerical snowpack models essentially rely on observational data for initialization, forcing, parametrization and validation. Such data are available in increasing amount, but the propagation of related uncertainties on the simulation results has received rather limited attention so far. Depending on their complexity, even small errors can have a profound effect on simulations, which dilutes our confidence in the results. This paper aims at quantification of the overall and fractional contributions of some archetypical measurement uncertainties on snowpack simulations in Arctic environments.


Figure 1: The map shows the location of Kongsvegen glacier and the position of the automatic weather stations KNG8, KNG6 and KNG1 (Norwegian Polar Institute, 2014).


 

The sensitivity pattern is studied at two sites representing the accumulation and ablation area of Kongsvegen glacier (Svalbard) using the snowpack scheme Crocus. The contribution of measurement errors on model output variance, either alone or by interaction, is decomposed using Global Sensitivity Analysis. This allows investigation of the temporal evolution of the fractional contribution of different factors on key model output metrics, which provides a more detailed understanding of the model’s sensitivity pattern. The analysis demonstrates that the specified uncertainties in precipitation and longwave radiation forcings had an outstanding influence on the calculated surface height changes and surface energy balance components. The model output sensitivity patterns also reveal some characteristic seasonal imprints. For example, uncertainties in longwave radiation trace on calculated surface energy balance rather continuous throughout the year and at both study sites, while precipitation exerts most influence during winter and at the upper site. Such findings are valuable to identify critical parameters and improve their measurement and correspondingly updated simulations may shed new light on the confidence of results from snow or glacier mass and energy balance models. This matters for many applications in the fields of e.g. avalanche and hydrological forecasting.


Figure 2: Yearly averaged total-order effects (on model uncertainty) of factors on surface height change (SHC), surface energy balance (SEB), sensible heat and latent heat flux (HF) at KNG8 and KNG1. The whiskers show the 95% confidence interval derived from 1000 empirical bootstrap samples. The mean (taken over the whole period) of the 6 hourly first-order sums (linear effects) are given in the upper right corner.