Abstract. This paper describes ESM-SnowMIP, an international coordinated modelling effort to evaluate current snow schemes, including snow schemes that are included in Earth system models, in a wide variety of settings against local and global observations. The project aims to identify crucial processes and characteristics that need to be improved in snow models in the context of local- and global-scale modelling. A further objective of ESM-SnowMIP is to better quantify snow-related feedbacks in the Earth system. Although it is not part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6), ESM-SnowMIP is tightly linked to the CMIP6-endorsed Land Surface, Snow and Soil Moisture Model Intercomparison (LS3MIP).

In this study, we assess the impact of forcing data errors, model structure, and parameter choices on 1‐D snow simulations simultaneously within a global variance‐based sensitivity analysis framework. This approach allows inclusion of interaction effects, drawing a more representative picture of the resulting sensitivities. We utilize all combinations of a multiphysics snowpack model to mirror the influence of model structure. Uncertainty ranges of model parameters and input data are extracted from the literature. We evaluate a suite of 230,000 model realizations at the snow monitoring station Kühtai (Tyrol, Austria, 1,920 m above sea level) against snow water equivalent observations. The results show throughout the course of 25 winter seasons (1991–2015) and different model performance criteria a large influence of forcing data uncertainty and its interactions on the model performance. Mean interannual total sensitivity indices are in the general order of parameter choice < model structure < forcing error, with precipitation, air temperature, and the radiative forcings controlling the variance during the accumulation period and air temperature and longwave irradiance controlling the variance during the ablation period, respectively. Model skill is highly sensitive to the snowpack liquid water transport scheme throughout the whole winter period and to albedo representation during the ablation period. We found a sufficiently long evaluation period (>10 years) is required for robust averaging. A considerable interaction effect was revealed, indicating that an improvement in the knowledge (i.e., reduction of uncertainty) of one factor alone might not necessarily improve model results.

Abstract. In this study, the fully distributed, physically based hydroclimatological model AMUNDSEN is set up for catchments in the highly glacierized Ötztal Alps (Austria, 558 km2 in total). The model is applied for the period 1997–2013, using a spatial resolution of 50 m and a temporal resolution of 1 h. A novel parameterization for lateral snow redistribution based on topographic openness is presented to account for the highly heterogeneous snow accumulation patterns in the complex topography of the study region. Multilevel spatiotemporal validation is introduced as a systematic, independent, complete, and redundant validation procedure based on the observation scale of temporal and spatial support, spacing, and extent. This new approach is demonstrated using a comprehensive set of eight independent validation sources: (i) mean areal precipitation over the period 1997–2006 derived by conserving mass in the closure of the water balance, (ii) time series of snow depth recordings at the plot scale, (iii–iv) multitemporal snow extent maps derived from Landsat and MODIS satellite data products, (v) the snow accumulation distribution for the winter season 2010/2011 derived from airborne laser scanning data, (vi) specific surface mass balances for three glaciers in the study area, (vii) spatially distributed glacier surface elevation changes for the entire area over the period 1997–2006, and (viii) runoff recordings for several subcatchments. The results indicate a high overall model skill and especially demonstrate the benefit of the new validation approach. The method can serve as guideline for systematically validating the coupled components in integrated snow-hydrological and glacio-hydrological models.

Abstract. Seasonal snow cover is an important temporary water storage in high-elevation regions. Especially in remote areas, the available data are often insufficient to accurately quantify snowmelt contributions to streamflow. The limited knowledge about the spatiotemporal variability of the snowmelt isotopic composition, as well as pronounced spatial variation in snowmelt rates, leads to high uncertainties in applying the isotope-based hydrograph separation method. The stable isotopic signatures of snowmelt water samples collected during two spring 2014 snowmelt events at a north-and a south-facing slope were volume weighted with snowmelt rates derived from a distributed physicsbased snow model in order to transfer the measured plotscale isotopic composition of snowmelt to the catchment scale. The observed δ 18 O values and modeled snowmelt rates showed distinct inter-and intra-event variations, as well as marked differences between north-and south-facing slopes. Accounting for these differences, two-component isotopic hydrograph separation revealed snowmelt contributions to streamflow of 35 ± 3 and 75 ± 14 % for the early and peak melt season, respectively. These values differed from those determined by formerly used weighting methods (e.g., using observed plot-scale melt rates) or considering either the north-or south-facing slope by up to 5 and 15 %, respectively.

[1] Runoff generation in Alpine regions is typically affected by snow processes. Snow accumulation, storage, redistribution, and ablation control the availability of water. In this study, several robust parameterizations describing snow processes in Alpine environments were implemented in a fully distributed, physically based hydrological model. Snow cover development is simulated using different methods from a simple temperature index approach, followed by an energy balance scheme, to additionally accounting for gravitational and wind-driven lateral snow redistribution. Test site for the study is the Berchtesgaden National Park (Bavarian Alps, Germany) which is characterized by extreme topography and climate conditions. The performance of the model system in reproducing snow cover dynamics and resulting discharge generation is analyzed and validated via measurements of snow water equivalent and snow depth, satellite-based remote sensing data, and runoff gauge data. Model efficiency (the Nash-Sutcliffe coefficient) for simulated runoff increases from 0.57 to 0.68 in a high Alpine headwater catchment and from 0.62 to 0.64 in total with increasing snow model complexity. In particular, the results show that the introduction of the energy balance scheme reproduces daily fluctuations in the snowmelt rates that trace down to the channel stream. These daily cycles measured in snowmelt and resulting runoff rates could not be reproduced by using the temperature index approach. In addition, accounting for lateral snow transport changes the seasonal distribution of modeled snowmelt amounts, which leads to a higher accuracy in modeling runoff characteristics.