Abstract. Physically based multilayer snowpack models suffer from various modelling errors. To represent these errors, we built the new multiphysical ensemble system ES-CROC (Ensemble System Crocus) by implementing new representations of different physical processes in the deterministic coupled multilayer ground/snowpack model SUR-FEX/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.

Abstract. Physically based multilayer snowpack models suffer from various modelling errors. To represent these errors, we built the new multi-physical ensemble system ESCROC (Ensemble System Crocus) by implementing new representations of different physical processes in the deterministic coupled multi-layer 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 dataset. 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 observations uncertainties is necessary to obtain an optimal dispersion, but this issue can also be a consequence of the fact that meterorological forcing uncertainties were not accounted for. ESCROC is a promising system to integrate numerical snow modelling errors in ensemble forecasting and ensemble assimilation systems in support of avalanche hazard forecasting and other snowpack modelling applications.

Abstract. Monitoring the evolution of the snowpack properties in mountainous areas is crucial for avalanche hazard forecasting and water resources management. In-situ and remotely sensed observations provide precious information on the snowpack but usually offer a limited spatio-temporal coverage of bulk or surface variables only. In particular, visible-near infrared (VIS-NIR) reflectance observations can inform on the snowpack surface properties but are limited by terrain shading and clouds. Snowpack modelling enables to estimate any physical variable, virtually anywhere, but is affected by large errors and uncertainties. Data assimilation offers a way to combine both sources of information, and to propagate information from observed areas to non observed areas. Here, we present CrocO, (Crocus-Observations) an ensemble data assimilation system able to ingest any snowpack observation (applied as a first step to the height of snow (HS) and VIS-NIR reflectances) in a spatialised geometry. CrocO uses an ensemble of snowpack simulations to represent modelling uncertainties, and a Particle Filter (PF) to reduce them. The PF is known to collapse when assimilating a too large number of observations. Two variants of the PF were specifically implemented to ensure that observations information is propagated in space while tackling this issue. The global algorithm ingests all available observations with an iterative inflation of observation errors, while the klocal algorithm is a localised approach performing a selection of the observations to assimilate based on background correlation patterns. Experiments are carried out in a twin experiment setup, with synthetic observations of HS and VIS-NIR reflectances available in only a 1/6th of the simulation domain. Results show that compared against runs without assimilation, analyses exhibit an average improvement of snow water equivalent Continuous Rank Probability Score (CRPS) of 60 % when assimilating HS with a 40-member ensemble, and an average 20 % CRPS improvement when assimilating reflectance with a 160-member ensemble. Significant improvements are also obtained outside the observation domain. These promising results open a way for the assimilation of real observations of reflectance, or of any snowpack observations in a spatialised context.

Uncertainties of snowpack models and of their meteorological forcings limit their use by avalanche hazard forecasters, or for glaciological and hydrological studies. The spatialized simulations currently available for avalanche hazard forecasting are only assimilating sparse meteorological observations. As suggested by recent studies, their forecasting skills could be significantly improved by assimilating satellite data such as snow reflectances from satellites in the visible and the near-infrared spectra. Indeed, these data can help constrain the microstructural properties of surface snow and light absorbing impurities content, which in turn affect the surface energy and mass budgets. This paper investigates the prerequisites of satellite data assimilation into a detailed snowpack model. An ensemble version ofMétéo-France operational snowpack forecasting system (named S2M) was built for this study. This operational system runs on topographic classes instead of grid points, so-called 'semi-distributed' approach. Each class corresponds to one of the 23 mountain massifs of the French Alps (about 1000km 2 each), an altitudinal range (by step of 300m) and aspect (by step of 45 o ). We assess the feasability of satellite data assimilation in such a semi-distributed geometry. Ensemble simulations are compared with satellite observations from MODIS and Sentinel-2, and with in-situ reflectance observations. The study focuses on the 2013-2014 and 2016-2017 winters in the Grandes-Rousses massif. Substantial Pearson R 2 correlations (0.75-0.90) of MODIS observations with simulations are found over the domain. This suggests that assimilating it could have an impact on the spatialized snowpack forecasting system. However, observations contain significant biases (0.1-0.2 in reflectance) which prevent their direct assimilation. MODIS spectral band ratios seem to be much less biased. This may open the way to an operational assimilation of MODIS reflectances into the Météo-France snowpack modelling system. Highlights -Ensemble simulations of the snowpack are compared with satellite reflectances -Spatial aggregation into the semi-distributed geometry filters the observation noises -Satellite reflectances carry useful information worth to assimilate -MODIS reflectances can not be directly assimilated because they are biased -Ratios of MODIS reflectances show no evidence of bias and could be assimilated