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

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

Abstract.A long term data acquisition effort of profile soil moisture is under way in southwestern France at 13 automated weather stations. This ground network was developed in order to validate remote sensing and model soil moisture estimates. In this paper, both those in situ observations and a synthetic data set covering continental France are used to test a simple method to retrieve root zone soil moisture from a time series of surface soil moisture information. A recursive exponential filter equation using a time constant, T , is used to compute a soil water index. The Nash and Sutcliff coefficient is used as a criterion to optimise the T parameter for each ground station and for each model pixel of the synthetic data set. In general, the soil water indices derived from the surface soil moisture observations and simulations agree well with the reference root-zone soil moisture. Overall, the results show the potential of the exponential filter equation and of its recursive formulation to derive a soil water index from surface soil moisture estimates. This paper further investigates the correlation of the time scale parameter T with soil properties and climate conditions. While no significant relationship could be determined between T and the main soil properties (clay and sand fractions, bulk density and organic matter content), the modelled spatial variability and the observed inter-annual variability of T suggest that a weak climate effect may exist.

[1] Thirty-three snowpack models of varying complexity and purpose were evaluated across a wide range of hydrometeorological and forest canopy conditions at five Northern Hemisphere locations, for up to two winter snow seasons. Modeled estimates of snow water equivalent (SWE) or depth were compared to observations at forest and open sites at each location. Precipitation phase and duration of above-freezing air temperatures are shown to be major influences on divergence and convergence of modeled estimates of the subcanopy snowpack. When models are considered collectively at all locations, comparisons with observations show that it is harder to model SWE at forested sites than open sites. There is no universal ''best'' model for all sites or locations, but comparison of the consistency of individual model performances relative to one another at different sites (and vice versa). Calibration of models at forest sites provides lower errors than uncalibrated models at three out of four locations. However, benefits of calibration do not translate to subsequent years, and benefits gained by models calibrated for forest snow processes are not translated to open conditions.

Abstract:Two downscaling methods designed for the study of the hydrological impact of climate change on the Seine basin in France are tested for present climate. First, a multivariate statistical downscaling (SD) methodology based on weather typing and conditional resampling is described. Then, a bias correction technique for dynamical downscaling based on quantile-quantile mapping is introduced. To evaluate the end-to-end SD methodology, the atmospheric forcing derived from the large-scale circulation (LSC) of the ERA40 reanalysis by SD is used to force a hydrological model. Simulated discharges reproduce historical values reasonably well. Next, the dynamical and statistical approaches are compared using the Météo-France ARPEGE general circulation model in a variable resolution configuration (resolution around 60 km over France). The ARPEGE simulation is downscaled using the two methodologies, and hydrological simulations are performed. Regarding downscaled temperature and precipitation, the statistical approach is more efficient in reproducing the temporal and spatial autocorrelation properties. The simulated river discharges from the two approaches are nevertheless very similar: the two methods reproduce well the seasonal cycle and the daily distribution of streamflows. Finally, the results of the study are discussed from a practical impact study perspective.

[1] The hydrometeorological model SIM consists of a meteorological analysis system (SAFRAN), a land surface model (ISBA), and a hydrogeological model (MODCOU). It generates atmospheric forcing at an hourly time step, and it computes water and surface energy budgets, the river flow at more than 900 river-gauging stations, and the level of several aquifers. SIM was extended over all of France in order to have a homogeneous nationwide monitoring of the water resources: it can therefore be used to forecast flood risk and to monitor drought risk over the entire nation. The hydrometeorological model was applied over a 10-year period from 1995 to 2005. In this paper the databases used by the SIM model are presented; then the 10-year simulation is assessed by using the observations of daily streamflow, piezometric head, and snow depth. This assessment shows that SIM is able to reproduce the spatial and temporal variabilities of the water fluxes. The efficiency is above 0.55 (reasonable results) for 66% of the simulated river gauges, and above 0.65 (rather good results) for 36% of them. However, the SIM system produces worse results during the driest years, which is more likely due to the fact that only few aquifers are simulated explicitly. The annual evolution of the snow depth is well reproduced, with a square correlation coefficient around 0.9 over the large altitude range in the domain. The streamflow observations were used to estimate the overall error of the simulated latent heat flux, which was estimated to be less than 4%.

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

Abstract. The overall objective of the present study is to introduce the new ECOCLIMAP-II database for Europe, which is an upgrade for this region of the former initiative, ECOCLIMAP-I, already implemented at global scale. The ECOCLIMAP programme is a dual database at 1 km resolution that includes an ecosystem classification and a coherent set of land surface parameters that are primarily mandatory in meteorological modelling (notably leaf area index and albedo). Hence, the aim of this innovative physiography is to enhance the quality of initialisation and impose some surface attributes within the scope of weather forecasting and climate related studies. The strategy for implementing ECOCLIMAP-II is to depart from prevalent land cover products such as CLC2000 (Corine Land Cover) and GLC2000 (Global Land Cover) by splitting existing classes into new classes that possess a better regional character by virtue of the climatic environment (latitude, proximity to the sea, topography). The leaf area index (LAI) from MODIS and normalized difference vegetation index (NDVI) from SPOT/Vegetation (a global monitoring system of vegetation) yield the two proxy variables that were considered here in order to perform a multi-year trimmed analysis between 1999 and 2005 using the K-means method. Further, meteorological applications require each land cover type to appear as a partition of fractions of 4 main surface types or tiles (nature, water bodies, sea, urban areas) and, inside the nature tile, fractions of 12 plant functional types (PFTs) representing generic vegetation types -principally broadleaf forest, needleleaf forest, C3 and C4 crops, grassland and bare land -as incorporated by the SVAT model ISBA (Interactions Surface Biosphere Atmosphere) developed at Météo France. This landscape division also forms the cornerstone of a validation exercise. The new ECOCLIMAP-II can be verified with auxiliary land cover products at very fine and coarse resolutions by means of versatile land occupation nomenclatures.