Actes - Climato.be

25 ème Colloque de l’Association Internationale de Climatologie, Grenoble 20121. Donnés utilisées et méthodologieLes données journalières de hauteur de neige enregistrées à 0700 UTC dans les 111stations automatiques et manuelles (figure 1) ont été utilisées pour calculer les valeursmoyennes mensuelles et saisonnières de précipitation neigeuse et du nombre de jours avecprécipitation solide dans le cas où au moins 80 % des données journalières était disponible,autrement la valeur moyenne a été rejetée (Klein Tank et al., 2002). La techniqued'interpolation par krigeage a été utilisée pour estimer les champs moyens saisonniers deprécipitation neigeuse cumulée et de jours neigeux en Piémont au cours de la dernièredécennie 2000-2010, pendant laquelle la plupart des stations était opérationnelle etfonctionnait en continu. Au pas de temps mensuel, on a calculé l'épaisseur moyenne de laneige au sol et on a utilisé les données satellitaires de couverture neigeuse MODIS Terra à500 m pour en déduire des informations sur la présence/absence de la neige au sol. Lesproduits de couverture neigeuse MODIS sur 8 jours ont été agrégés sur une période d'un mois(Terzago et al., 2010) afin de déterminer la couverture neigeuse maximale mensuelle et deminimiser le numéro des pixels non-classés à cause de la présence de nuages. Les pixels nonneigeuxdans le produit mensuel ont été considérés comme des « stations fictives » avechauteur de neige observée égale à 0. Pour des raisons de temps de calcul, un sous-ensemble deces points choisi aléatoirement a été retenu pour l'interpolation. Puisque la quantité de neigeaccumulée au sol dépend linéairement de l’altitude, l'interpolation spatiale a été effectuée enutilisant un krigeage avec dérive externe (KED) qui intègre dans le modèle une tendancelinéaire avec l'altitude. Les données saisonnières ont été lissées sur des blocs de 20 x 20kilomètres carrés et interpolées sur une grille de résolution 2,5 km.Figure 1 : Localisation des 111 stations nivo-méteorologiques manuelles et automatiques du réseau de l'AgenceRégionale de Protection Environnementale (ARPA) du Piémont (a) et couverture neigeuse maximale dedécembre 2006 dérivée du MODIS (b).L’interpolation a été évaluée par validation croisée. Les valeurs prédites ont été comparéesaux différentes valeurs observées disponibles pour la validation en utilisant trois statistiquesproposées par Hunter et Meentmeyer (2005) : i) le coefficient de détermination R² pourévaluer le degré de variabilité dans les valeurs observées expliquées par le modèle, ii) l'erreurabsolue moyenne (MAE) pour déterminer l'amplitude globale de l'erreur et iii) l'erreurmoyenne (ME) afin d'identifier les tendances de sous/surestimation ou la directionnalité del’erreur globale. Les résultats ont été comparés avec ceux obtenus en interpolant seulement lesdonnées des stations météorologiques pour déterminer la valeur ajoutée de l'utilisation des740