Abstract. According to the living data process in ESSD, this publication presents extensions of a comprehensive hydrometeorological and glaciological data set for several research sites in the Rofental (1891–3772 m a.s.l., Ötztal Alps, Austria). Whereas the original dataset has been published in a first original version in 2018 (https://doi.org/10.5194/essd-10-151-2018), the new time series presented here originate from meteorological and snow-hydrological recordings that have been collected from 2017 to 2020. Some data sets represent continuations of time series at existing locations, others come from new installations complementing the scientific monitoring infrastructure in the research catchment. Main extensions are a fully equipped automatic weather and snow monitoring station, as well as extensive additional installations to enable continuous observation of snow cover properties. Installed at three high Alpine locations in the catchment, these include automatic measurements of snow depth, snow water equivalent, volumetric solid and liquid water content, snow density, layered snow temperature profiles, and snow surface temperature. One station is extended by a particular arrangement of two snow depth and water equivalent recording devices to observe and quantify wind-driven snow redistribution. They are installed at nearby wind-exposed and sheltered locations and are complemented by an acoustic-based snow drift sensor. The data sets represent a unique time series of high-altitude mountain snow and meteorology observations. We present three years of data for temperature, precipitation, humidity, wind speed, and radiation fluxes from three meteorological stations. The continuous snow measurements are explored by combined analyses of meteorological and snow data to show typical seasonal snow cover characteristics. The potential of the snow drift observations are demonstrated with examples of measured wind speeds, snow drift rates and redistributed snow amounts in December 2019 when a tragic avalanche accident occurred in the vicinity of the station. All new data sets are provided to the scientific community according to the Creative Commons Attribution License by means of the PANGAEA repository (https://www.pangaea.de/?q=%40ref104365).

Abstract Twenty-seven models participated in the Earth System Model–Snow Model Intercomparison Project (ESM-SnowMIP), the most data-rich MIP dedicated to snow modeling. Our findings do not support the hypothesis advanced by previous snow MIPs: evaluating models against more variables and providing evaluation datasets extended temporally and spatially does not facilitate identification of key new processes requiring improvement to model snow mass and energy budgets, even at point scales. In fact, the same modeling issues identified by previous snow MIPs arose: albedo is a major source of uncertainty, surface exchange parameterizations are problematic, and individual model performance is inconsistent. This lack of progress is attributed partly to the large number of human errors that led to anomalous model behavior and to numerous resubmissions. It is unclear how widespread such errors are in our field and others; dedicated time and resources will be needed to tackle this issue to prevent highly sophisticated models and their research outputs from being vulnerable because of avoidable human mistakes. The design of and the data available to successive snow MIPs were also questioned. Evaluation of models against bulk snow properties was found to be sufficient for some but inappropriate for more complex snow models whose skills at simulating internal snow properties remained untested. Discussions between the authors of this paper on the purpose of MIPs revealed varied, and sometimes contradictory, motivations behind their participation. These findings started a collaborative effort to adapt future snow MIPs to respond to the diverse needs of the community.

Abstract. Thirty-year simulations of seasonal snow cover in 22 physically based models driven with bias-corrected meteorological reanalyses are examined at four sites with long records of snow observations. Annual snow cover durations differ widely between models but interannual variations are strongly correlated because of the common driving data. No significant trends are observed in starting dates for seasonal snow cover, but there are significant trends towards snow cover ending earlier at two of the sites in observations and most of the models. A simplified model with just two parameters controlling solar radiation and sensible heat contributions to snowmelt spans the ranges of snow cover durations and trends. This model predicts that sites where snow persists beyond annual peaks in solar radiation and air temperature will experience rapid decreases in snow cover duration with warming as snow begins to melt earlier and at times of year with more energy available for melting.

Abstract. The 30-year simulations of seasonal snow cover in 22 physically based models driven with bias-corrected meteorological reanalyses are examined at four sites with long records of snow observations. Annual snow cover durations differ widely between models, but interannual variations are strongly correlated because of the common driving data. No significant trends are observed in starting dates for seasonal snow cover, but there are significant trends towards snow cover ending earlier at two of the sites in observations and most of the models. A simplified model with just two parameters controlling solar radiation and sensible heat contributions to snowmelt spans the ranges of snow cover durations and trends. This model predicts that sites where snow persists beyond annual peaks in solar radiation and air temperature will experience rapid decreases in snow cover duration with warming as snow begins to melt earlier and at times of year with more energy available for melting.

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.

In this paper, the hydrological impacts of future socio-economic and climatic development are assessed for a regional-scale Alpine catchment (Brixental, Tyrol, Austria). Therefore, coupled storylines of future land use and climate scenarios were developed in a transdisciplinary stakeholder process by means of questionnaire analyses and interviews with local experts from various relevant societal sectors. Resulting future land use maps for each decade were used as spatial input in the hydrological model WaSiM, to which a new module for the consideration of snow-canopy interaction processes has been added. Simulation results for three developed storylines, each combined with a moderate (A1B) and an extreme (RCP8.5) climate future, show that in a warmer and dryer climate the amount of annual simulated streamflow at the gauge of the catchment undergoes a significant reduction. The (mainly natural) reforestation of the catchment - caused by abandonment of previously cultivated areas - leads to additional losses of water by enhanced interception and evapotranspiration processes. Further cultivation of the current mountain pasture areas has a certain potential to attenuate undesirable long-term impacts of climate change on the catchment water balance.

This paper presents a novel data fusion technique for improving the snow cover monitoring for a mesoscale Alpine region, in particular in those areas where two information sources disagree. The presented methodological innovation consists in the integration of remote-sensing data products and the numerical simulation results by means of a machine learning classifier (support vector machine), capable to extract information from their quality measures. This differs from the existing approaches where remote sensing is only used for model tuning or data assimilation. The technique has been tested to generate a time series of about 1300 snow maps for the period between October 2012 and July 2016. The results show an average agreement between the fused product and the reference ground data of 96%, compared to 90% of the moderate-resolution imaging spectroradiometer (MODIS) data product and 92% of the numerical model simulation. Moreover, one of the most important results is observed from the analysis of snow cover area (SCA) time series, where the fused product seems to overcome the well-known underestimation of snow in forest of the MODIS product, by accurately reproducing the SCA peaks of winter season.

Mountain regions with complex orography are a particular challenge for regional climate simulations. High spatial resolution is required to account for the high spatial variability in meteorological conditions. This study presents a very high-resolution regional climate simulation (5 km) using the Weather Research and Forecasting Model (WRF) for the central part of Europe including the Alps. Global boundaries are dynamically downscaled for the historical period 1980–2009 (ERA-Interim and MPI-ESM), and for the near future period 2020–2049 (MPI-ESM, scenario RCP4.5). Model results are compared to gridded observation datasets and to data from a dense meteorological station network in the Berchtesgaden Alps (Germany). Averaged for the Alps, the mean bias in temperature is about −0.3 °C, whereas precipitation is overestimated by +14% to +19%. R 2 values for hourly, daily and monthly temperature range between 0.71 and 0.99. Temporal precipitation dynamics are well reproduced at daily and monthly scales (R 2 between 0.36 and 0.85), but are not well captured at hourly scale. The spatial patterns, seasonal distributions, and elevation-dependencies of the climate change signals are investigated. Mean warming in Central Europe exhibits a temperature increase between 0.44 °C and 1.59 °C and is strongest in winter and spring. An elevation-dependent warming is found for different specific regions and seasons, but is absent in others. Annual precipitation changes between −4% and +25% in Central Europe. The change signals for humidity, wind speed, and incoming short-wave radiation are small, but they show distinct spatial and elevation-dependent patterns. On large-scale spatial and temporal averages, the presented 5 km RCM setup has in general similar biases as EURO-CORDEX simulations, but it shows very good model performance at the regional and local scale for daily meteorology, and, apart from wind-speed and precipitation, even for hourly values.

Geochemical and isotopic tracers were often used in mixing models to estimate glacier melt contributions to streamflow, whereas the spatio‐temporal variability in the glacier melt tracer signature and its influence on tracer‐based hydrograph separation results received less attention. We present novel tracer data from a high‐elevation catchment (17 km2, glacierized area: 34%) in the Oetztal Alps (Austria) and investigated the spatial, as well as the subdaily to monthly tracer variability of supraglacial meltwater and the temporal tracer variability of winter baseflow to infer groundwater dynamics. The streamflow tracer variability during winter baseflow conditions was small, and the glacier melt tracer variation was higher, especially at the end of the ablation period. We applied a three‐component mixing model with electrical conductivity and oxygen‐18. Hydrograph separation (groundwater, glacier melt, and rain) was performed for 6 single glacier melt‐induced days (i.e., 6 events) during the ablation period 2016 (July to September). Median fractions (±uncertainty) of groundwater, glacier melt, and rain for the events were estimated at 49±2%, 35±11%, and 16±11%, respectively. Minimum and maximum glacier melt fractions at the subdaily scale ranged between 2±5% and 76±11%, respectively. A sensitivity analysis showed that the intraseasonal glacier melt tracer variability had a marked effect on the estimated glacier melt contribution during events with large glacier melt fractions of streamflow. Intra‐daily and spatial variation of the glacier melt tracer signature played a negligible role in applying the mixing model. The results of this study (a) show the necessity to apply a multiple sampling approach in order to characterize the glacier melt end‐member and (b) reveal the importance of groundwater and rainfall–runoff dynamics in catchments with a glacial flow regime.

Abstract A distributed snow model is applied to simulate the spatiotemporal evolution of the Austrian snow cover at 1 km × 1 km spatial and daily temporal resolution for the period 1948–2009. After a comprehensive model validation, changes in snow cover conditions are analyzed for all of Austria as well as for different Austrian subregions and elevation belts focusing on the change in snow cover days (SCDs). The comparison of SCDs for the period 1950–79 to those achieved for 1980–2009 for all of Austria shows a decrease in SCDs with a maximum of >35 SCDs near Villach (Carinthia). The analysis of SCD changes in different subregions of Austria reveals mean changes between −11 and −15 days with highest absolute change in SCDs for southern Austria. Two decrease maxima could be identified in elevations of 500–2000 m MSL (between −13 and −18 SCDs depending on the subregion considered) and above 2500 m MSL (over −20 SCDs in the case of central Austria). The temporal distribution of SCD change in the Austrian subregions is characterized by a reduction of SCDs in midwinter and at the end of winter rather than by fewer SCDs in early winter. With respect to the temporal distribution of SCD change in different elevation belts, changes in elevations below 1000 m MSL are characterized by a distinct reduction of SCDs in January. With increasing elevation the maximum change in SCDs shifts toward the summer season, reaching a maximum decrease in the months of June–August above 2500 m MSL.

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).