SURFEX is a new externalized land and ocean surface platform that describes the surface fluxes and the evolution of four types of surfaces: nature, town, inland water and ocean. It is mostly based on pre-existing, well-validated scientific models that are continuously improved. The motivation for the building of SURFEX is to use strictly identical scientific models in a high range of applications in order to mutualise the research and development efforts. SURFEX can be run in offline mode (0-D or 2-D runs) or in coupled mode (from mesoscale models to numerical weather prediction and climate models). An assimilation mode is included for numerical weather prediction and monitoring. In addition to momentum, heat and water fluxes, SURFEX is able to simulate fluxes of carbon dioxide, chemical species, continental aerosols, sea salt and snow particles. The main principles of the organisation of the surface are described first. Then, a survey is made of the scientific module (including the coupling strategy). Finally, the main applications of the code are summarised. The validation work undertaken shows that replacing the pre-existing surface models by SURFEX in these applications is usually associated with improved skill, as the numerous scientific developments contained in this community code are used to good advantage.

Detailed studies of snow cover processes require models that offer a fine description of the snow cover properties. The detailed snowpack model Crocus is such a scheme, and has been run operationally for avalanche forecasting over the French mountains for 20 yr. It is also used for climate or hydrological studies. To extend its potential applications, Crocus has been recently integrated within the framework of the externalized surface module SURFEX. SURFEX computes the exchanges of energy and mass between different types of surface and the atmosphere. It includes in particular the land surface scheme ISBA (Interactions between Soil, Biosphere, and Atmosphere). It allows Crocus to be run either in stand-alone mode, using a time series of forcing meteorological data or in fully coupled mode (explicit or fully implicit numerics) with atmospheric models ranging from meso-scale models to general circulation models. This approach also ensures a full coupling between the snow cover and the soil beneath. Several applications of this new simulation platform are presented. They range from a 1-D stand-alone simulation (Col de Porte, France) to fully-distributed simulations in complex terrain over a whole mountain range (Massif des Grandes Rousses, France), or in coupled mode such as a surface energy balance and boundary layer simulation over the East Antarctic Ice Sheet (Dome C).

This paper presents the Meso-NH model version 5.4. Meso-NH is an atmospheric non hydrostatic research model that is applied to a broad range of resolutions, from synoptic to turbulent scales, and is designed for studies of physics and chemistry. It is a limited-area model employing advanced numerical techniques, including monotonic advection schemes for scalar transport and fourth-order centered or odd-order WENO advection schemes for momentum. The model includes state-of-the-art physics parameterization schemes that are important to represent convective-scale phenomena and turbulent eddies, as well as flows at larger scales. In addition, Meso-NH has been expanded to provide capabilities for a range of Earth system prediction applications such as chemistry and aerosols, electricity and lightning, hydrology, wildland fires, volcanic eruptions, and cyclones with ocean coupling. Here, we present the main innovations to the dynamics and physics of the code since the pioneer paper of Lafore et al. (1998) and provide an overview of recent applications and couplings.

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.

The Crocus snowpack model within the Interactions between Soil–Biosphere–Atmosphere (ISBA) land surface model was run over northern Eurasia from 1979 to 1993, using forcing data extracted from hydrometeorological datasets and meteorological reanalyses. Simulated snow depth, snow water equivalent, and density over open fields were compared with local observations from over 1000 monitoring sites, available either once a day or three times per month. The best performance is obtained with European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim). Provided blowing snow sublimation is taken into account, the simulations show a small bias and high correlations in terms of snow depth, snow water equivalent, and density. Local snow cover durations as well as the onset and vanishing dates of continuous snow cover are also well reproduced. A major result is that the overall performance of the simulations is very similar to the performance of existing gridded snow products, which, in contrast, assimilate local snow depth observations. Soil temperature at 20-cm depth is reasonably well simulated. The methodology developed in this study is an efficient way to evaluate different meteorological datasets, especially in terms of snow precipitation. It reveals that the temporal disaggregation of monthly precipitation in the hydrometeorological dataset from Princeton University significantly impacts the rain–snow partitioning, deteriorating the simulation of the onset of snow cover as well as snow depth throughout the cold season.

In alpine regions, wind-induced snow transport strongly influences the spatio-temporal evolution of the snow cover throughout the winter season. To gain understanding on the complex processes that drive the redistribution of snow, a new numerical model is developed. It directly couples the detailed snowpack model Crocus with the atmospheric model Meso-NH. Meso-NH/Crocus simulates snow transport in saltation and in turbulent suspension and includes the sublimation of suspended snow particles. The coupled model is evaluated against data collected around the experimental site of Col du Lac Blanc (2720 m a.s.l., French Alps). First, 1-D simulations show that a detailed representation of the first metres of the atmosphere is required to reproduce strong gradients of blowing snow concentration and compute mass exchange between the snowpack and the atmosphere. Secondly, 3-D simulations of a blowing snow event without concurrent snowfall have been carried out. Results show that the model captures the main structures of atmospheric flow in alpine terrain. However, at 50 m grid spacing, the model reproduces only the patterns of snow erosion and deposition at the ridge scale and misses smaller scale patterns observed by terrestrial laser scanning. When activated, the sublimation of suspended snow particles causes a reduction of deposited snow mass of 5.3% over the calculation domain. Total sublimation (surface + blowing snow) is three times higher than surface sublimation in a simulation neglecting blowing snow sublimation.