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

Physically based multilayer snowpack models suffer from various modelling errors. To represent these errors, we built the new multiphysical ensemble system ESCROC (Ensemble System Crocus) by implementing new representations of different physical processes in the deterministic coupled multilayer ground/snowpack model SURFEX/ISBA/Crocus. This ensemble was driven and evaluated at Col de Porte (1325 m a.s.l., French alps) over 18 years with a high-quality meteorological and snow data set. A total number of 7776 simulations were evaluated separately, accounting for the uncertainties of evaluation data. The ability of the ensemble to capture the uncertainty associated to modelling errors is assessed for snow depth, snow water equivalent, bulk density, albedo and surface temperature. Different sub-ensembles of the ESCROC system were studied with probabilistic tools to compare their performance. Results show that optimal members of the ESCROC system are able to explain more than half of the total simulation errors. Integrating members with biases exceeding the range corresponding to observational uncertainty is necessary to obtain an optimal dispersion, but this issue can also be a consequence of the fact that meteorological forcing uncertainties were not accounted for. The ESCROC system promises the integration of numerical snow-modelling errors in ensemble forecasting and ensemble assimilation systems in support of avalanche hazard forecasting and other snowpack-modelling applications.

This work introduces the S2M (SAFRAN–SURFEX/ISBA–Crocus–MEPRA) meteorological and snow cover reanalysis in the French Alps, Pyrenees and Corsica, spanning the time period from 1958 to 2021. The simulations are made over elementary areas, referred to as massifs, designed to represent the main drivers of the spatial variability observed in mountain ranges (elevation, slope and aspect). The meteorological reanalysis is performed by the SAFRAN system, which combines information from numerical weather prediction models (ERA-40 reanalysis from 1958 to 2002, ARPEGE from 2002 to 2021) and the best possible set of available in situ meteorological observations. SAFRAN outputs are used to drive the Crocus detailed snow cover model, which is part of the land surface scheme SURFEX/ISBA. This model chain provides simulations of the evolution of the snow cover, underlying ground and the associated avalanche hazard using the MEPRA model. This contribution describes and discusses the main climatological characteristics (climatology, variability and trends) and the main limitations of this dataset. We provide a short overview of the scientific applications using this reanalysis in various scientific fields related to meteorological conditions and the snow cover in mountain areas. An evaluation of the skill of S2M is also displayed, in particular through comparison to 665 independent in situ snow depth observations. Further, we describe the technical handling of this open-access dataset, available at https://doi.org/10.25326/37#v2020.2. The S2M data are provided by Météo-France – CNRS, CNRM, Centre d'Études de la Neige, through AERIS (Vernay et al., 2022).

This paper describes the impact of the various changes made to the Safran–Isba–Modcou (SIM) hydrometeorological system and demonstrates that the new version of the model performs better than the previous one by making comparisons with observations of daily river flows and snow depths. SIM was developed and put into operational service at Météo-France in the early 2000s. The SIM application is dedicated to the monitoring of water resources and can therefore help in drought monitoring or flood risk forecasting on French territory. This complex system combines three models: SAFRAN, which analyses meteorological variables close to the surface, the ISBA land surface model, which aims to calculate surface fluxes at the interface with the atmosphere and ground variables, and finally MODCOU, a hydrogeological model which calculates river flows and changes in groundwater levels. The SIM model has been improved first by reducing the infrared radiation bias of SAFRAN and then by using the more advanced ISBA multi-layer surface diffusion scheme to have a more physical representation of surface and ground processes. In addition, more accurate and recent databases of vegetation, soil texture, and orography were used. Finally, in mountainous areas, a sub-grid orography representation using elevation bands was adopted, as was the possibility of adding a reservoir to represent the effect of aquifers in mountainous areas. The numerical simulations carried out with the SIM model covered the period from 1958 to 2018, thereby providing an extensive historical analysis of the water resources over France.

This article investigates the climatic response of a series of indicators for characterizing annual snow conditions and corresponding meteorological drivers at 1500 m altitude in the Chartreuse mountain range in the Northern French Alps. Past and future changes were computed based on reanalysis and observations from 1958 to 2016, and using CMIP5–EURO-CORDEX GCM–RCM pairs spanning historical (1950–2005) and RCP2.6 (4), RCP4.5 and RCP8.5 (13 each) future scenarios (2006–2100). The adjusted climate model runs were used to drive the multiphysics ensemble configuration of the detailed snowpack model Crocus. Uncertainty arising from physical modeling of snow accounts for 20 % typically, although the multiphysics is likely to have a much smaller impact on trends. Ensembles of climate projections are rather similar until the middle of the 21st century, and all show a continuation of the ongoing reduction in average snow conditions, and sustained interannual variability. The impact of the RCPs becomes significant for the second half of the 21st century, with overall stable conditions with RCP2.6, and continued degradation of snow conditions for RCP4.5 and 8.5, the latter leading to more frequent ephemeral snow conditions. Changes in local meteorological and snow conditions show significant correlation with global temperature changes. Global temperature levels 1.5 and 2 ∘C above preindustrial levels correspond to a 25 and 32 % reduction, respectively, of winter mean snow depth with respect to the reference period 1986–2005. Larger reduction rates are expected for global temperature levels exceeding 2 ∘C. The method can address other geographical areas and sectorial indicators, in the field of water resources, mountain tourism or natural hazards.

A multireplicate multimodel ensemble of hydrological simulations covering the 1860–2099 period has been produced for the Upper Durance River basin (French Alps). An original quasi‐ergodic analysis of variance was applied to quantify uncertainties related to General Circulation Models (GCMs), Statistical Downscaling Models (SDMs) and the internal variability of each GCM/SDM simulation chain. For temperature, GCM uncertainty prevails and SDM uncertainty is nonnegligible. Significant warming and in turn significant changes are predicted for evaporation, snow cover and seasonality of discharges. For precipitation, GCM and SDM uncertainty components are of the same order. A high contribution of the large and small‐scale components of internal variability is also obtained, inherited, respectively, from the GCMs and the different replicates of a given SDM. The same applies for annual discharge. The uncertainty in values that could be experienced for any given future period is therefore very high. For both discharge and precipitation, even the sign of future realizations is uncertain at a 90% confidence level. These findings have important implications. Similarly to GCM uncertainty, SDM uncertainty cannot be neglected. The same applies for both components of internal variability. Climate change impact studies based on a single SDM realization are likely to be no more relevant than those based on a single GCM run. They may lead to poor decisions for climate change adaptation.