Abstract: Using the original setup described in Gallée et al. (2013), the MAR regional climate model including a coupled snowpack/aeolian snow transport parameterization, was run at a fine spatial (5 km horizontal and 2 m vertical) resolution over 1 summer month in coastal Adélie Land. Different types of feedback were taken into account in MAR including drag partitioning caused by surface roughness elements. Model outputs are compared with observations made at two coastal locations, D17 and D47, situated respectively 10 and 100 km inland. Wind speed was correctly simulated with positive values of the Nash test (0.60 for D17 and 0.37 for D47) but wind velocities above 10 m s−1 were underestimated at both D17 and D47; at D47, the model consistently underestimated wind velocity by 2 m s−1. Aeolian snow transport events were correctly reproduced with the right timing and a good temporal resolution at both locations except when the maximum particle height was less than 1 m. The threshold friction velocity, evaluated only at D17 for a 7-day period without snowfall, was overestimated. The simulated aeolian snow mass fluxes between 0 and 2 m at D47 displayed the same variations but were underestimated compared to the second-generation FlowCaptTM values, as was the simulated relative humidity at 2 m above the surface. As a result, MAR underestimated the total aeolian horizontal snow transport for the first 2 m above the ground by a factor of 10 compared to estimations by the second-generation FlowCaptTM. The simulation was significantly improved at D47 if a 1-order decrease in the magnitude of z0 was accounted for, but agreement with observations was reduced at D17. Our results suggest that z0 may vary regionally depending on snowpack properties, which are involved in different types of feedback between aeolian transport of snow and z0.

Abstract: In the SURFEX/ISBA-Crocus multi-layer snowpack model, the snow microstructure has up to now been characterised by the grain size and by semi-empirical shape variables which cannot be measured easily in the field or linked to other relevant snow properties. In this work we introduce a new formulation of snow metamorphism directly based on equations describing the rate of change of the optical diameter (dopt). This variable is considered here to be equal to the equivalent sphere optical diameter, which is inversely proportional to the specific surface area (SSA). dopt thus represents quantitatively some of the geometric characteristics of a porous medium. Different prognostic rate equations of dopt, including a re-formulation of the original Crocus scheme and the parameterisations from Taillandier et al. (2007) and Flanner and Zender (2006), were evaluated by comparing their predictions to field measurements carried out at Summit Camp (Greenland) in May and June 2011 and at Col de Porte (French Alps) during the 2009/10 and 2011/12 winter seasons. We focused especially on results in terms of SSA. In addition, we tested the impact of the different formulations on the simulated density profile, the total snow height, the snow water equivalent (SWE) and the surface albedo. Results indicate that all formulations perform well, with median values of the RMSD between measured and simulated SSA lower than 10 m2 kg-1. Incorporating the optical diameter as a fully fledged prognostic variable is an important step forward in the quantitative description of the snow microstructure within snowpack models, because it opens the way to data assimilation of various electromagnetic observations.