RAMMS ou le miroir aux alouettes

Dans un article qui a été récemment accepté (septembre 2017) par la revue scientifique Journal of Glaciology [1], trois scientifiques ont décidé de porter l’estocade contre le modèle numérique d’avalanche RAMMS (acronyme pour Rapid Mass Movement) commercialisé par le WSL/SLF depuis le début des années 2010. RAMMS comporte trois modules (avalanche, chutes de pierre, et lave torrentielle). Le module avalanche a connu un certain succès commercial au cours des sept dernières années au point que l’on voit fleurir des affirmations dithyrambiques le citant comme LA référence en la matière. Le présent article revient sur les causes de la polémique et les conséquences pratiques de ce qui s’annonce comme un scandale dans un milieu d’ordinaire bien policé.

Les protagonistes

Plantons le décor. Les trois scientifiques sont

  • Dieter Issler, un consultant en avalanche pour le NGI en Norvège, physicien de formation (ETHZ) et ancien du WSL/SLF (il a notamment participé à la mise en place du site de la Sionne),
  • Jim Jenkins (prof., université de Cornell, un des pères de la théorie cinétique granulaire), et
  • Jim McElwaine (prof., université de Durham, un mathématicien avec une passion coupable pour la neige et les avalanches, et un partenaire de longue date du SLF).

La revue scientifique est Journal of Glaciology, dont un des éditeurs est Perry Bartelt du SLF et papa de RAMMS. Il y a publié l’essentiel des articles consacrés à son modèle numérique.

Hum… cela sent le linge sale qu’on lave en famille d’habitude. Alors qu’en est-il ?

Des tabous en science ?

Il y a en science des frontières que l’on se garde de franchir, des quasi-tabous. Le plus souvent, ce sont des frontières tacites acceptées par la communauté, mais parfois, le législateur est passé par là et met des garde-fous légaux (comme l’interdiction du clonage humain). Comme l’avait anticipé le philosophe français Jacques Ellul, ce type de frontière finit toujours par être franchi :

« parce que tout ce qui est technique, sans distinction de bien et de mal, s’utilise forcément quand on l’a en main. »

Dans le domaine des modèles numériques en général, et en dynamique des fluides géophysiques en particulier, les développeurs (le plus souvent des universitaires) peuvent mettre à disposition leur modèle pour des applications ponctuelles ou des comparaisons avec des données de terrain ou de laboratoire. La diffusion large d’un produit de la recherche nécessite un gros travail d’interface (souvent de pré- et post-traitement de l’information) et de documentation. Et le travail est tellement conséquent qu’il faut des années pour que de petites équipes académiques puissent y parvenir. Dans bien des cas, les outils restent donc sous forme de version beta sans interface et c’est l’utilisateur qui doit se débrouiller. Rendre un code librement accessible ou le commercialiser sous-tend que le produit correspond à un état de l’art accepté par la communauté scientifique et qu’il a été validé à différents niveaux.

Valider un code numérique

Le travail de validation d’un code numérique comporte plusieurs étapes :

  • vérifier que les équations utilisées sont pertinentes à décrire le phénomène étudié,
  • contrôler que la méthode de résolution numérique est stable et converge vers la bonne solution physique,
  • appliquer le code à des cas idéalisés (pour montrer que le code capture la physique du phénomène avec la précision souhaitée) puis des cas réels.

Généralement, le processus de développement et de validation fait l’objet d’articles dans des revues scientifiques. Le travail doit être examiné par des scientifiques (rapporteurs) indépendants. C’est ce qu’on appelle le processus de revue par les pairs (peer-reviewing). Et naturellement, pour éviter les conflits d’intérêt, les auteurs ne soumettent pas leurs articles dans des revues dans lesquelles ils ont des responsabilités éditoriales ou un lien étroit avec les responsables de la publication.

Et alors quel est le problème avec les avalanches ?

Tous les modèles actuels sont fondés sur des équations de Saint-Venant, un jeu d’équations tirées de l’hydraulique pour décrire les rivières en crue. Le premier code numérique moderne d’avalanche a été écrit par Jean-Paul Vila dans les années 1980 [2-3], et tous les codes actuels en sont inspirés à des degrés divers. Cependant, même après 30 ans de tests et développements, les progrès ont été limités. En 2000, dans le cadre d’un programme européen de recherche, différents groupes ont calibré les paramètres de frottement de leur modèle numérique (tous fondés sur le même jeu d’équations de Saint-Venant) et ont abouti à des valeurs très différentes [4]. Cela a souligné la grande sensibilité des résultats aux détails de la résolution numérique.

D’autres problèmes se sont posés et touchent le cœur de l’approche actuelle. Par exemple, depuis Voellmy [5], le frottement de l’avalanche est supposé inclure du frottement solide (de type Coulomb) et du frottement turbulent (fonction de la vitesse quadratique). Avec la multiplication de sites expérimentaux permettant de suivre le mouvement des avalanches dans le détail, il est devenu clair que le modèle de Voellmy ne permet guère de décrire fidèlement le frottement réel [6-7]. Cependant, pour des raisons tant pratiques qu’historiques, le modèle de Voellmy continue à être utilisé par les praticiens. Ceux-ci appliquent donc leur propre technique pour ajuster les paramètres du modèle de Voellmy.

Malgré les progrès, les méthodes numériques peinent parfois à résoudre correctement les équations. Par exemple, la localisation du front pose souvent problème (en effet c’est le point où la hauteur écoulement devient nulle, mais avec des modèles de frottement de type Voellmy, le frottement devient localement infini (car il s’écrit en 1/h avec h la hauteur d’écoulement) et même avec un frottement faible, le front est une singularité qui nécessite un traitement particulier, voir [8-9]). Pour plusieurs codes, l’avalanche ne s’arrête jamais. La vitesse tend vers 0, sans toutefois l’atteindre. Il faut donc des « patches » pour corriger ces erreurs numériques.

Et RAMMS alors ?

Plusieurs groupes dans le monde se sont intéressés à décrire des écoulements de fluides complexes à l’aide d’équations de type Saint-Venant. Il faut prendre conscience que l’étalement d’un fluide ordinaire comme de la mayonnaise sur une table de cuisine est déjà un problème d’une grande complexité, sans solution actuellement ! Alors, des fluides complexes comme la neige offrent un défi considérable.

Le groupe de Perry Bartelt au SLF a tenté de relever le défi en puisant des idées dans la thermodynamique des systèmes hors équilibre de Prigogine et la théorie cinétique des gaz dissipatifs. Le point d’achoppement est que le modèle (ou devrait-on dire « les modèles » tant le modèle initial a subi des évolutions au cours des 10 dernières années) n’est pas cohérent avec ce que l’on sait des fluides complexes tels que les milieux granulaires. Les démonstrations sont par ailleurs émaillées d’erreurs. C’est en substance ce que dit l’article très détaillé d’Issler, Jenkins, et McElwaine. Iverson et George notaient, de leur côté, des problèmes similaires [10].

À l’absence de validation théorique des concepts développés s’est ajouté le manque de validation expérimentale. Indépendamment de toutes autres considérations, un problème de taille surgit pour toute validation avec la neige : il est très difficile de mesurer les propriétés mécaniques de la neige en écoulement par des expériences indépendantes et contrôlées. Les paramètres de tous les modèles sont donc ajustés ou calibrés en faisant correspondre des observations de terrain et des prédictions du modèle. C’est ce qu’on appelle caler un modèle. Mais caler les paramètres d’un modèle sur un événement ne dit rien sur la pertinence du modèle. Cela ne permet pas non plus de dire grand-chose sur un autre événement (même sur le même site) car a priori les valeurs des paramètres changent sans qu’on sache comment. Dans sa dernière version, RAMMS possède 16 paramètres. Avec un modèle sur-paramétré (c.-à-d. qui possède plus de paramètres qu’il n’existe de variables indépendantes que l’on puisse observer et utiliser pour tester le modèle), il est plus facile de caler le modèle, mais il devient quasiment impossible de l’utiliser en pratique de façon rigoureuse. Une erreur de détermination de quelques % sur un paramètre peut se propager et s’amplifier, aboutissant finalement à des erreurs considérables dans la prédiction. Un problème bien connu des hydrologues… [11].

Le miroir aux alouettes

Malgré l’absence de fondement théorique reconnu au sein de la communauté scientifique et l’insuffisance du travail de validation à partir d’observations de terrain, RAMMS a réussi à séduire autorités locales et bureaux d’études. Pourtant, des avertissements sur l’illusion de précision de l’approche numérique ont été émis. En 2004, à la fin de sa carrière au SLF, Bruno Salm (qui a beaucoup contribué aux méthodes quantitatives pour les avalanches et peut être considéré comme un des grands acteurs de la nivologie moderne) a mis en garde contre les voix de sirène : avec l’augmentation croissante de la puissance des ordinateurs, les modèles gagnent en réalisme, mais non en précision [12] :

« tous les modèles présentés sont – jusqu’à présent – d’une façon ou d’une autre entachés d’incertitudes. Un accroissement de la complexité des modèles n’implique pas nécessairement une précision plus grande ou une meilleure stratégie de lutte contre les dangers naturels. »

Pour les autorités locales, le numérique semble offrir non seulement un extraordinaire gain en précision (c.-à-d. prédire les distances d’arrêt pour le zonage d’avalanche avec une précision métrique), mais également un outil objectif permettant un traitement homogène du zonage de risque (implicitement il est supposé que le dire d’expert est subjectif). C’est ainsi qu’en 2011, les deux ingénieurs généraux de l’administration française Le Gallou et Guignard ont préconisé le recours systématique au numérique dans le zonage d’avalanche pour éviter les écarts qu’ils ont notés d’une commune à l’autre [13]. RAMMS est jusqu’à présent le seul candidat en lice, puisque personne d’autre que le SLF n’a commercialisé de code de calcul avec une une large diffusion. Il existe bien un candidat autrichien (SAMOS), mais il a une diffusion confidentielle [14-15], et il y a aussi quelques codes académiques, mais ils sont utilisés par leurs propres développeurs (en France, il y a l’IRSTEA de Grenoble, qui est le pionnier mondial du domaine, mais également l’ENS de Lyon et l’Université de Savoie). Étonnamment, personne dans l’administration ne s’est demandé pourquoi aucun des autres groupes de recherche dans le monde ne s’est lancé dans le créneau du code commercial. Les successeurs de Vila seraient-ils moins bons ou audacieux que leur père spirituel ? Ou bien suffisamment prudents pour ne pas se lancer sur une pente glissante ? Un responsable du RTM évoquait l’idée que, certes, un modèle comme RAMMS n’est pas rigoureux sur le plan scientifique, mais il est socialement acceptable car il répond à un besoin de la société (il veut dire à un besoin d’un groupe de personnes de l’administration, la synecdoque étant une figure de style prisée dans l’administration). La situation n’est pas sans rappeler les financiers qui avant la crise financière de 2007 s’étaient engagés dans des placements hasardeux mais tellement prometteurs en termes de rendement, et qui une fois le crack boursier advenu, se sont écriés que jamais au grand jamais ils n’auraient pu se douter du caractère extravagant des montages financiers. Personne ne les en avait avertis. Pourtant les voix s’étaient bien élevés avant 2007 pour dénoncer l’emballement financier, mais comme le dit le proverbe populaire, il n’y a pas plus sourd que celui qui ne veut pas entendre. Dans la finance, entre des promesses de rendement à 3 % et d’autres à 15 %, le choix a été vite fait, avec les conséquences que l’on mesure encore aujourd’hui. Nous vivons une situation similaire avec les modèles hydrauliques et d’avalanches. L’État a dépensé plus de 450 k€ pour des simulations numériques d’une trentaine de couloirs à Chamonix, simulations réalisées à l’aide de RAMMS, en dépit des remarques sur la faible plus-value de l’exercice. Des rapports ont été rendus récemment par les administrations pour « encadrer » l’usage des modèles numériques. Étonnamment, là encore, aucun scientifique (du développeur de code au mécanicien en passant par le géographe, le sociologue, l’historien ou le spécialiste de terrain) n’a été consulté. On mesure en France et en Suisse le fossé croissant qui sépare milieu académique et administrations en charge de la gestion territoriale.

Pour les bureaux d’études, le numérique est une formidable porte d’entrée à un marché qui pendant longtemps nécessitait des années de formation et de reconnaissance par les pairs. Il suffit désormais d’acquérir la licence d’un code numérique pour du même coup bénéficier d’un savoir numérique. Or, si dans quelques domaines de l’ingénieur (les fluides industriels ou le calcul des structures), une grande partie du travail de calcul est réalisée de nos jours par des ordinateurs, avec une partie humaine réduite à la portion congrue, tant les codes sont fiables, cela n’est pas le cas dans d’autres domaines (dont l’hydraulique torrentielle et les avalanches) où les modèles sont beaucoup plus rustiques et leurs résultats nécessitent un examen critique.

Et demain ?

Comme l’a documenté l’excellent ouvrage de Nicolas Chevassus-au-Louis [16], on note qu’avec la pression croissante qu’ont les scientifiques à publier, il y a en retour un accroissement du nombre d’articles litigieux (selon Chevassus-au-Louis, cela concerne quand même une fraction non négligeable de l’ordre de quelques % de l’ensemble des articles publiés). Cela va d’erreurs d’inadvertance (erreurs de mesure, de calcul, etc.) à des problèmes plus graves de fraude ou de plagiat. Quand les errata ne suffisent pas à corriger des erreurs manifestes, on voit des scientifiques retirer leurs articles. Dans le cas présent, il est encore trop tôt pour savoir quelle position sera adoptée par Bartelt et ses coauteurs, mais une réponse à l’article d’Issler, Jenkins et McElwaine est nécessaire (comme cela a été le cas avec la critique d’Iverson et George [10] sur un autre article de Bartelt). La rétraction est une possibilité ouverte pour clore la polémique. Les conditions d’acceptation des articles de Bartelt dans Journal of Glaciology font planer un doute et ne sont pas sans rappeler le scandale qui a secoué le groupe de Claude Allègre à l’IPGP et leurs articles dans une revue publiée par Elsevier [17].

Si on dépasse cette querelle interne au milieu académique, on peut se demander si le scandale va ébranler la diffusion du code numérique, s’il va être retiré du marché ou gelé dans une version plus cohérente avec l’état de l’art scientifique réel (et non idéalisé). Tout du moins, un sérieux rappel sur les limites de la modélisation s’impose. Dans la méthode dite VSG de 1990 (VSG pour Voellmy Salm Gubler) qui est en quelque sorte l’ancêtre de RAMMS, les auteurs prévenaient l’utilisateur [18]

« Il faut mettre en garde sur l’utilisation de cette méthode comme un « livre de cuisine » (on portera une attention particulière à l’usage des mots « environ » et « valeur indicative »). Bien au contraire, chaque avalanche requiert un traitement propre et une connaissance de la matière dans son ensemble. À cela s’ajoute aussi l’expérience qui joue un rôle important et que chacun doit acquérir. Il faut en dernier lieu avoir à l’esprit que les avalanches échappent aux calculs trop précis. »

C’est un avertissement qui ne mange pas de pain. Toutefois, si on examine les scandales dans le domaine sanitaire, on constate une grande inertie des intervenants. En pratique, il faut beaucoup de temps avant que les choses ne bougent…

Le scandale ne sera pas sans incidence plus large. La probabilité qu’un zonage tout numérique aboutisse prochainement à un accident dans le monde est grande. Comment cela sera-t-il jugé ? Qui sera responsable ? Ne risque-t-on pas une défiance envers les approches scientifiques ? Les sondages d’opinion montrent qu’une part croissante de citoyens éprouve de la défiance envers la science et les « experts » (combien parmi nous font confiance aux experts des OGM ou du nucléaire ?) alors que la fascination technologique est toujours bien présente. Comment redonner de la confiance aux utilisateurs et autorités car si l’approche numérique n’est pas la panacée, elle a ses mérites ?

Pour les maîtres d’ouvrage et maîtres d’œuvre, un regard critique doit être porté sur la plus-value d’une simulation numérique (par rapport à une expertise classique) et sur la capacité du bureau d’études à manipuler un code de type RAMMS. Il faut indiquer qu’un bon praticien saura déjouer la plupart des pièges numériques alors qu’un praticien inexpérimenté prendra les solutions numériques pour argent comptant. Le rendu sera souvent proche (de jolies couleurs sur une carte), mais le diable se cache dans les détails.

Mise à jour (27 janvier 2018)

Perry Bartelt a usé de son droit de réponse en publiant un article pour nier toute erreur dans ses articles [18]. Dans la conclusion, il dit créer une nouvelle mécanique :

« Ici, nous devons être francs : nous n’avons pas négligé les solutions offertes par la mécanique granulaire et l’étude des courants gravitaires de particules. La situation est beaucoup plus grave : nous les avons abandonnées pour la simple raison que nous avons constaté qu’elles ne fonctionnaient pas. Confrontés à la tâche de développer un modèle d’avalanche pour la pratique en ingénierie, nous nous sommes posé une question simple : qu’est-ce que l’application des lois de Newton et Stokes implique lorsqu’on s’attelle à la description des conditions aux limites au fond, du frottement de l’écoulement, et de l’aérosol. Nous voulions éviter les formulations de modèle qui abondent dans le jargon inutile et trompeur, l’inexactitude, la contradiction et l’illogisme général. Nous voulions avant tout un modèle qui puisse être testé en utilisant des mesures expérimentales et qui soit compréhensible pour les ingénieurs et les technologues. [19] »

Références

  • [1] Issler, D, J. T. Jenkins, and J. McElwaine, Comments on avalanche flow models based on extensions of the concept of random kinetic energy, à paraître dans Journal of Glaciology. Article accessible depuis le dépôt arxiv : arxiv.org/abs/1710.00524
  • [2] Vila, J.-P., Modélisation mathématique et simulation d’écoulements à surface libre . La Houille Blanche, 6/7, 485-489, 1984.
  • [3] Vila, J.P., Sur la théorie et l’approximation numérique des problèmes hyperboliques non-linéaires, application aux équations de Saint-Venant et à la modélisation des avalanches denses, Ph.D. thesis thesis, Paris VI, 1986.
  • [4] Barbolini, M., U. Gruber, C.J. Keylock, M. Naaim, and F. Savi, Application of statistical and hydraulic-continuum dense-snow avalanche models to five European sites, Cold Regions Science and Technology, 31, 133-149, 2000.
  • [5] Voellmy, A., Über die Zerstörungskraft von Lawinen, Schweizerische Bauzeitung, 73, 159–165, 1955.
  • [6] Ancey, C., and M. Meunier, Estimating bulk rheological properties of flowing snow avalanches from field data, Journal of Geophysical Research, 109, F01004, 2004.
  • [7] Gauer, P., and K. Kristensen, Four decades of observations from NGI’s full-scale avalanche test site Ryggfonn—Summary of experimental results, Cold Regions Science and Technology, 125, 162-176, 2016.
  • [8] Ancey, C., S. Cochard, S. Wiederseiner, and M. Rentschler, Existence and features of similarity solutions for supercritical non-Boussinesq gravity currents, Physica D, 226, 32-54, 2007.
  • [9] Ancey, C., S. Cochard, S. Wiederseiner, and M. Rentschler, Front dynamics of supercritical non-Boussinesq gravity currents, Water Resources Research, 42, W08424, 2006
  • [10] Iverson, R.M., and D.L. George, Discussion of “The relation between dilatancy, effective stress and dispersive pressure in granular avalanches” by P. Bartelt and O. Buser, Acta Geotechnica, 11, 1465–1468, 2016.
  • [11] Beven, K., Environmental Modelling: An Uncertain Future, Routledge (Taylor & Francis), Oxon, 2009.
  • [12] Salm, B., A short and personal history of snow avalanche dynamics, Cold Regions Science and Technology, 39, 83-92, 2004.
  • [12] Le Gallou, J.-Y., and P. Guignard, Modalités de prise en compte des avalanches exceptionnelles pour améliorer la prévention des risques et renforcer la sécurité des personnes, Inspection générale de l’administration, Paris, 2011.
  • [14] Fischer, J.-T., A. Kofler, W. Fellin, M. Granig, and K. Kleemayr, Optimization of Computational Snow Avalanche Simulation Tools, in International Snow Science Workshop 2014 Proceedings, pp. 665-669, Banff, Canada, 2014.
  • [15] Sampl, P., and M. Granig, Avalanche simulation with SAMOS-AT, in Proceedings of the International Snow Science Workshop, 27 September–2 October 2009, pp. 519-523, Davos, Switzerland, 2009.
  • [16] Chevassus-au-Louis, N., Malscience, Éditions du Seuil, Paris, 2016.
  • [17] Foucart, S et H. Morin, Petits arrangements entre géologues, Le Monde du 26 décembre 2008.
  • [18] Salm, B., A. Burkard, and H. Gubler, Berechnung von Fliesslawinen, eine Anleitung für Praktiker mit Beispielen, Eidgenössisches Institut für Schnee- und Lawinenforschung (Davos), 1990.
  • [19] Bartelt, P., and O. Buser, Avalanche dynamics by Newton. Reply to comments on avalanche flow models based on the concept of random kinetic energy, Journal of Glaciology, doi: 10.1017/jog.2018.1, 2018.
  • [20] “Here we must be blunt: we have not disregarded the solutions offered by granular mechanics/particulate gravity currents. The situation is far more serious: we have abandoned them for the simple reason we found that they do not work. Confronted with the task of developing an avalanche model for engineering practice we asked ourselves a simple question: what does the consequential application of the laws of Newton and Stokes imply for the all-important descriptions of the basal boundary conditions, avalanche flow friction and the formation of the powder cloud. We wanted to avoid model formulations that abound in needless and misleading jargon, misstatement, contradiction and general illogic. Above all we wanted a model that can be tested using experimental measurements and that is comprehensible to engineers and technologists.” Tiré de la réf. [19], Sec. 5.



Avalanche

La première génération

En 1923, un ingénieur du RTM Paul Mougin proposait le premier modèle dynamique d’avalanche (une analogie simple entre un écoulement de neige et un bloc solide). Dans une série d’articles de 1955, l’ingénieur suisse Adolf Voellmy généralisait ce modèle et offrait un cadre de calcul plus complet des caractéristiques dynamiques (vitesse moyenne, pression cinétique) d’une avalanche. C’est le premier d’une longue série de modèles, qui ont été utilisés jusque dans les années 2000.

La seconde génération

En 1966, le chercheur suisse Bruno ébauche la première esquisse de ce qui allait être l’approche hydraulique des avalanches. L’analogie était fondée non plus sur un bloc glissant, mais une rivière en crue. En URSS, à la fin des années 1960 et dans les années 1970, des mécaniciens comme Sergei Grigorian et Margarita Eglit élaborent un cadre théorique relativement complet exploitant cette analogie. Pour lui donner une dimension opérationnelle encore fallait-il pouvoir calculer des solutions sur des cas d’intérêt pratique. À la fin des années 1970, ce sont deux ingénieurs du CTGREF (ce qui est aujourd’hui IRSTEA) de Grenoble, Rémi Pochat et Gérard Brugnot, développent les premières solutions numériques sur les toutes premières stations de calcul disponibles. Il faut toutefois attendre le travail de thèse de Jean-Paul Vila en 1986 pour voir apparaître le premier schéma numérique moderne (fondé sur la technique des volumes finis) capable de résoudre des modèles d’avalanche. Tous les modèles d’avalanche sont peu ou prou basés sur la même approche.

À quand la troisième génération ?

Et depuis 1986 quels sont les progrès ? Si la puissance de calcul a été considérablement améliorée, les rendus graphiques sont plus réalistes, force est de reconnaître que les progrès en matière de modélisation n’ont pas permis de déboucher sur une troisième génération de modèle. Et bien futé qui pourra dire à quoi ressemblera cette troisième génération…

Alors pourquoi tant de difficultés ? Contrairement à l’eau, la neige est un matériau compressible, au comportement mécanique complexe et variable dans le temps. Les conditions initiales, les caractéristiques du manteau neigeux, sa variabilité spatiale ou temporelle, sont autant de paramètres peu connus. Le travail réalisé sur le terrain (sur des sites comme la Sionne en Valais) ou sur des matériaux modèles (p. ex., les écoulements granulaires en laboratoire) ont permis de mieux appréhender la dynamique d’écoulements de fluides complexes, mais les progrès restent lents.




Hiver 2017

Un hiver précoce

L’hiver 2016-2017 commence précocement avec de belles chutes de neige jusqu’à basse altitude. Heureuse surprise pour moi qui monte depuis les Paccots au col du Lys le samedi 12 novembre 2016. Bonne neige et épaisseur suffisante. Mais une chute précoce n’est jamais bon signe pour la suite.

Préalpes fribourgeoises depuis le col du Lys
Figure 1 : le Pays d’en Haut vu depuis le col du Lys (FR) le 12 novembre 2016.

Dans les Alpes du Nord, l’isotherme 0 °C remonte rapidement de 1000 à 3000 m entre le 10 et le 16 novembre [1]. Il n’y a plus qu’à voir la neige fondre et attendre les prochaines chutes. Du 20 au 26 novembre, le flux passe au sud et amène de fortes précipitations sur la chaîne frontalière entre la France, la Suisse, et l’Italie. Les chutes de neige sont abondantes, principalement en haute Maurienne en France et dans la région du Simplon en Suisse. Pour le Simplon, le cumul de neige sur 6 jours se situe entre 1 et 2 m à 2400 m [1]. Selon Météo-Suisse [2], l’épisode de flux de sud du 20 au 25 novembre 2016 est le troisième épisode le plus long depuis 1981. Il est notamment plus long que celui de novembre 1996, qui avait également occasionné de fortes chutes de neige le long de la frontière italienne et dans les Alpes du Sud (avalanche causant des dommages à Orcières Merlette). En revanche, si les quantités de neige ont été très importantes, elles n’ont pas été exceptionnelles pour les Alpes Suisses puisque aucun cumul supérieur à 250 cm n’a été observé.


Figure 2 : le cirque de l’Arcelle vu depuis hélicoptère après l’avalanche du 25 novembre 2016. Source : Service des pistes de Val Cenis.

En haute Maurienne, les fortes précipitations provoquent une activité avalancheuse d’ampleur en haute altitude. Le télésiège de l’Arcelle (Val Cenis) est touché le 25 novembre 2016 par la phase aérosol d’une grosse avalanche. Il est difficile de donner des valeurs précises d’enneigement car la station était fermée à cette époque.
Ce que l’on sait c’est que la limite pluie/neige oscille probablement autour de 2500 m le lundi 21, puis s’abaisse à 2000 m les 22 et 23. En l’absence de données nivométriques sur l’Arcelle, on peut se faire une idée de la vigueur de l’épisode en reprenant les valeurs fournies par Météo-France sur les postes voisins :

  • à Bessans (1715 m), il est tombé 287 mm de précipitations du 18 au 25 novembre 2016, donc 202 mm le 24 novembre (veille de l’accident) ;
  • au col du Mont-Cenis (1715 m), il est tombé 284 mm de précipitations du 21 au 26 novembre 2016, donc 200 mm le 24 novembre ;
  • au col du Mont-Cenis (1280 m), il est tombé 245 mm de précipitations du 21 au 26 novembre 2016, donc 188 mm le 24 novembre.

Ce sont des valeurs tout à fait exceptionnelles et assez homogènes pour la haute Maurienne puisque la veille de l’accident, tous les postes ont mesuré une lame d’eau d’environ 200 mm en 24 h (du 24 au 25 novembre 7 h TU). Pour donner un ordre de grandeur du caractère exceptionnel, le cumul de pluie sur Termignon est de 1100 mm ; il est donc tombé en 24 h ce qu’il tombe en moyenne en deux mois. La valeur record de précipitation au col du Mont-Cenis (pluviographe EDF) était de 146 mm. L’analyse statistique des précipitations donne une période de retour d’environ 300 ans pour la pluie journalière du 24 novembre 2016 !

Pour le cirque de l’Arcelle, on peut estimer que le cumul sur 5 jours devait être voisin de 2 m, mais l’oscillation de l’isotherme du 0 °C et le régime de lombarde associé aux précipitations ont entraîné des variations très importantes dans la répartition des précipitations de neige tandis que l’effet de foehn a conduit à une hausse importante des températures dès qu’on s’éloigne de la zone frontalière.


Figure 3 : gare d’arrivée du TSD l’Arcelle vue depuis hélicoptère après l’avalanche du 25 novembre 2016. Source : Service des pistes de Val Cenis.


Figure 4 : dommage à la gare (endommagement des surfaces en plexiglas du tympan). Source : Service des pistes de Val Cenis.

Quand on examine les données disponibles pour la Suisse [2], on note qu’au cours des 35 dernières années, il y a eu 12 événements significatifs, soit un en moyenne tous les 3 ans, mais leur distribution est très variable : 2 épisodes pour les années 1980, 5 pour les années 1990, 2 pour les années 2000, et 3 déjà pour l’actuelle décennie. Pour la haute Maurienne, nous n’avons pas de statistiques, mais la chronique des gros accidents d’avalanche montre également qu’il y a une distribution très variable (les hivers 1972 et 1978 n’ont pas connu d’équivalents au cours des deux dernières décennies même si localement comme en 2013 des fortes chutes de neige ont pu se produire). Dans notre compréhension actuelle du changement climatique, nous redoutons toutefois la multiplication de tels épisodes : la mer Méditerranée est plus chaude, elle peut pulser de l’air chaud et humide vers les Alpes, et si cet air rencontre un front froid ou des masses d’air froid, de fortes chutes de neige sont attendues.

Les Pyrénées dans la tourmente

Le mois de décembre 2016 est un mois de disette. La neige se fait attendre un peu partout. Il faut attendre le 10 janvier pour commencer à voir la neige tomber en abondance [4]. À Cauterets dans les Pyrénées (65), les chutes de neige sont importantes : du mardi 10 janvier 2017 au lundi 16 janvier 2017, on observe 280 cm de cumul en 7 jours, l’épisode est le troisième plus important (assez loin du record de 440 cm de février 2013, mais sur 14 jours). La période de retour de tout l’épisode est de T = 15 ans. Plus remarquable est la valeur des chutes de neige du 14 au 16 janvier, avec 175 cm sur 3 jours, qui se rapproche du record de 2013 (220 cm). La période de retour est T = 15 à 20 ans d’après l’analyse nivo-météorologique statistique. Sur cette base, on peut estimer que l’épisode de neige de janvier 2017, est parmi les plus forts que la station a subis au cours des 35 dernières années.


Figure 5 : tapis du Grand Yéti endommagé après l’avalanche du 16 janvier 2017. Source : Service des pistes de Cauterets.

L’analyse des chutes de neige montre que (i) les trois plus gros cumuls se sont produits à une date récente (2013 et 2016) ; (ii) on compte 17 événements ayant amené plus de 150 cm de neige depuis 183 (en moyenne un tous les deux ans), mais parmi ceux-ci 8 se sont produit au cours des 5 dernières saisons. Il y a donc eu un accroissement très sensible (d’un facteur 7) de la fréquence d’occurrence des gros épisodes de neige puis que sur la période 1983-2012, on compte 8 événements en 30 ans, soit une probabilité annuelle de 0,26 alors que pour les 5 dernières années, on en compte 9, soit une probabilité annuelle de 1,8.

Tableau 1 : liste des épisodes de neige ayant amené plus de 150 cm sur Cauterets à 1800 m depuis déc. 1983.
DateDurée de l'épisodeCumul de neigePériode de retour (ans)
5 avril 198691532
3 février 1994131532
1 mars 199581592
24 mars 201441622
18 janvier 201471652
18 mars 200781763
15 mars 1985101843
26 février 201471953
23 mars 199282175
21 février 201572325
21 décembre 199392517
13 janvier 201362598
28 mars 198792608
26 février 20161530814
9 janvier 20171031315
20 janvier 2013932116
2 février 20131444068

L’origine de cette soudaine augmentation de la fréquence des gros cumuls n’est pas connue avec certitude. Il peut s’agir d’une « série noire ». Le cas se produit lorsque des scénarios météorologiques se maintiennent durant quelques années. Il y a ici peu d’éléments en faveur de cette interprétation, mais statistiquement cela reste dans le domaine du possible. Il peut s’agir d’une conséquence du réchauffement climatique en cours. Le conflit entre masses d’air chaud et humide et air froid entraîne, dans certains secteurs, un accroissement sensible de l’intensité des précipitations. Ce scénario est plausible, notamment car il se voit sur plusieurs zones sous influence maritime. Dans ce cas de figure, la station doit se préparer à une multiplication sensible des « gros coups de neige », donc des situations avalancheuses sur le site.

Un flux de sud dévastateur

Après une fin janvier et un mois de février où petites chutes de neige et coups de foehn alternent, le mois de mars voit une offensive de l’hiver avec de fortes chutes de neige dans les Alpes suisses ; on estime à un peu près 1 m de neige le cumul sur la Suisse occidentale durant cet épisode [5]. Avec le fort vent de sud, il y a eu un fort transport de neige et donc localement les quantités de neige cumulées ont pu être significativement plus importantes. Comme en témoigne l’avalanche de la Croix de Fer au-dessus du Châtelard (VS), les épaisseurs des cassures sont parfois impressionnantes. Le SLF déclenche le 8 mars une des plus grosses avalanches observées au cours de la dernière décennie [6].

avalanche du lac de Catogne
Figure 6 : avalanche du 7 ou 8 mars 2017 dans la combe des Jeurs, lac de Catogne sous la Croix de Fer, au-dessus du Châtelard (VS). Source : Jérôme Descomes-Sevoie.

Le hameau (non habité durant l’hiver) du Van d’en Haut sur le territoire de la commune de Salvan (VS) a été touché sévèrement par une avalanche issue due la pointe de la Djoua le mercredi 8 mars 2017 [7]. Dans le val Ferret et la vallée de la Dranse, il y a eu une activité avalancheuse de grande ampleur comme le montrent les dommages à la forêt à Champex ou au chalet d’alpage dans le val Ferret.

Van d'en Haut, commune de Salvan
Figure 8 : avalanche du mercredi 8 mars 2017 au Van d’en Haut, commune de Salvan (VS). Source : PiX4D.

alpage du mont Percé
Figure 9 : alpage du mont Percé (2023 m) dans le Val Ferret (commune d’Orsières, VS) endommagé par une avalanche venue du mont Fourchon à une date inconnue, sans doute le 8 mars 2017. Cliché pris en avril 2017.

combe du Barmay
Figure 10 : combe du Barmay sous les clochers d’Arpette (commune de Champex, VS) en avril 2017.

Avalanche du Van-d’en-Haut

Je montre ici quelques photographies prises le vendredi 7 avril 2017 en compagnie du bureau Geoformer de Brigue et Silva Plus en charge de la requalification du risque sur le Van d’en Haut.

.portelement_2 .play-icon.youtube-icon { background: url(http://www.toraval.ch/wp-content/plugins/portfolio-gallery/assets/images/admin_images/play.youtube.png) center center no-repeat; background-size: 30% 30%; } .portelement_2 .play-icon.vimeo-icon { background: url(http://www.toraval.ch/wp-content/plugins/portfolio-gallery/assets/images/admin_images/play.vimeo.png) center center no-repeat; background-size: 30% 30%; } .portelement_2 .play-icon { position: absolute; top : 0; left : 0; width: 100%; height: 100%; } .portelement_2 .add-H-relative { position: relative; } .portelement_2 { width: 277px !important; max-width: calc(100% - 10px); height:auto; margin: 5px 0 5px 0; float: left; overflow: hidden; position: relative; outline:none; background:#f7f7f7; border:1px solid #D0D0D0; box-sizing: border-box; } .default-block_2 { position:relative;; width:100%; } .default-block_2 .image-block_2 { margin :0; padding :0; line-height :0; border-bottom:1px solid #D0D0D0; } .default-block_2 img { margin :0 !important; padding :0 !important; width: 100%; border-radius :0; display: inline-block; } .default-block_2 .title-block_2 { display:block; height:auto; padding:10px 0 10px 0; width:100%; text-overflow: ellipsis; } .default-block_2 .title-block_2 h3 { position:relative; margin :0 !important; padding :0 5px 0 5px !important; max-width:275px !important; width: 100%; text-overflow: ellipsis; overflow: hidden; /*white-space:nowrap;*/ font-weight:normal; color:#555555; font-size:15px !important; line-height:19px !important; } .wd-portfolio-panel_2 { position: relative; display:block; width: calc(100% - 10px); margin:0 5px 0 5px; padding :0; text-align:left; } .wd-portfolio-panel_2 .description-block_2 p,.wd-portfolio-panel_2 .description-block_2 * { text-align:justify; font-size:14px !important; color:#5b5b5b; margin :0 !important; padding :0 !important; } .wd-portfolio-panel_2 .description-block_2 h1, .wd-portfolio-panel_2 .description-block_2 h2, .wd-portfolio-panel_2 .description-block_2 h3, .wd-portfolio-panel_2 .description-block_2 h4, .wd-portfolio-panel_2 .description-block_2 h5, .wd-portfolio-panel_2 .description-block_2 h6, .wd-portfolio-panel_2 .description-block_2 p, .wd-portfolio-panel_2 .description-block_2 strong, .wd-portfolio-panel_2 .description-block_2 span { padding:2px !important; margin :0 !important; } .wd-portfolio-panel_2 .description-block_2 ul, .wd-portfolio-panel_2 .description-block_2 li { padding:2px 0 2px 5px; margin :0 0 0 8px; } .wd-portfolio-panel_2 .thumbs-list_2 { list-style:none; clear:both; display:table; width:100%; padding :0; margin:3px 0 0 0; text-align:center; } .wd-portfolio-panel_2 .thumbs-list_2 li { display:inline-block; margin :0 3px 0 2px; } .wd-portfolio-panel_2 .thumbs-list_2 li a { display:block; width:75px; height:75px; opacity:0.7; display:table; border: none; box-shadow: none; } .wd-portfolio-panel_2 .thumbs-list_2 li a:hover { opacity:1; } .wd-portfolio-panel_2 img { margin :0 !important; padding :0 !important; display:table-cell; vertical-align:middle; width:75px !important; max-height:75px !important; width:100%; height:100%; } .wd-portfolio-panel_2 > div { padding-top:10px; margin-bottom:10px; background:url('http://www.toraval.ch/wp-content/plugins/portfolio-gallery/assets/images/admin_images/divider.line.png') center top repeat-x; } .wd-portfolio-panel_2 .button-block { padding-top:10px; margin-bottom:10px; } .wd-portfolio-panel_2 .button-block a, .wd-portfolio-panel_2 .button-block a:link, .wd-portfolio-panel_2 .button-block a:visited { padding:10px; display:inline-block; font-size:14px; background:#e74c3c; color:#ffffff; padding:6px 12px; text-decoration:none; } .wd-portfolio-panel_2 .button-block a:hover, .wd-portfolio-panel_2 .button-block a:focus, .wd-portfolio-panel_2 .button-block a:active { background:#df2e1b; color:#ffffff; text-decoration:none; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 { position: relative; overflow: hidden; float:;margin-top:5px;max-width:180px;width:20%;display:inline-block;margin-left:1%; opacity: 0; margin-bottom: 10px; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul { margin : 0 !important; padding : 0 !important; list-style: none; } #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul { margin : 0 !important; padding : 0 !important; overflow: hidden; width: 100%; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul { float: left; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul li { border-radius: 0px; list-style-type: none; margin : 0 !important; padding: 0; border: 1px solid #ccc;} #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul li a { background-color: #F7F7F7 !important; border-radius: 0px; font-size:14px !important; color:#555555 !important; text-decoration: none; cursor: pointer; margin : 0 !important; display: block; padding:3px; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul li a:hover { background-color: #FF3845 !important; color:#ffffff !important; cursor: pointer; } #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 { position: relative; overflow: hidden; float:;margin-top:5px;max-width:180px;width:20%;display:inline-block;margin-left:1%; opacity: 0;} #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li { list-style-type: none; border-radius: 0px; border: 1px solid #ccc;} #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a { font-size:14px !important; color:#555555 !important; background-color: #F7F7F7 !important; border-radius: 0px; padding:3px; display: block; text-decoration: none; } #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a:hover { color:#ffffff !important; background-color: #FF3845 !important; cursor: pointer } #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li.active a, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li.active a:link, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li.active a:visited, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li.active a:hover, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li.active a:focus, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li.active a:active { color:#ffffff !important; background-color: #FF3845 !important; cursor: pointer; } #huge_it_portfolio_content_2 section { position:relative; display:block; } #huge_it_portfolio_content_2 #huge_it_portfolio_container_2 { width: 79%; max-width: 100% !important; opacity: 0;} #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 #port-sort-direction{ position: static; } @media screen and (max-width: 768px) { #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a:link, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a:visited { font-size: 2vw !important; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul li a { font-size:2vw !important; } } @media screen and (max-width: 480px) { #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 { float: left; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 #sort-by{ float: left; width: 100% !important; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 #port-sort-direction{ float: left; width: 100% !important; position: relative; padding-left: 31% !important; right: 31%; } #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a:link, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a:visited { font-size: 3vw !important; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul li a { line-height: 3vw; font-size:3vw !important; } } @media screen and (max-width: 420px) { #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a:link, #huge_it_portfolio_content_2 #huge_it_portfolio_filters_2 ul li a:visited { font-size: 4vw !important; } #huge_it_portfolio_content_2 #huge_it_portfolio_options_2 ul li a { font-size:4vw !important; } } @media screen and (max-width: 317px) { .portelement_2 { width: calc(98% - 2px); margin: 1% !important; float: left; overflow: hidden; outline:none; } } #huge-it-container-loading-overlay_2 { width: 100%; height: 100%; position: absolute; z-index: 1; background: url(http://www.toraval.ch/wp-content/plugins/portfolio-gallery/assets/images/loading/loading-1.svg) center top ; background-repeat: no-repeat; background-size: 60px auto; } #huge_it_portfolio_options_and_filters_2{ position: relative; float: left; width: 20%; max-width: 180px; float:; margin-left:1%;}

Remerciements

  • Philippe Dupla, directeur Espaces Cauterets
  • Yves Flammier, service des pistes Val Cenis
  • André Burkard, Geoformer

Sauf mention contraire, les photographies sont de Christophe Ancey.