Les données utilisées pour ce projet ont été gracieusement mises à notre disposition par l’Agence Spatiale Européenne (ESA) au travers du projet C1P 8902 et traitées avec le logiciel GAMMA. Les données topographiques proviennent de l’Institut Géographique National (IGN) et les données météorologiques de MétéoFrance. Nous tenons également à remercier vivement les quatre relecteurs de cet article qui ont grandement contribué à son amélioration.
1Le pergélisol se définit comme un matériel de sub-surface dont la température est inférieure ou égale à 0 °C durant une année entière (Delaloye, 2004). Ce phénomène purement thermique se rencontre en France à des altitudes supérieures à environ 2 500 m, voire 3 000 m, et son étude se fait par l’intermédiaire de mesures de températures du sous-sol via des forages ponctuels et onéreux ou par des observations indirectes. Les formations superficielles des versants alpins en zone de pergélisol sont souvent affectées par des processus de fluage, donnant naissance à des formes telles que les glaciers rocheux qui sont des indices de la présence de pergélisol. Ces derniers sont définis comme l’expression visible, dans des matériaux non-consolidés, de la déformation lente et continue d’un corps de pergélisol sursaturé en glace (Delaloye, 2004). On distingue différents types de glaciers rocheux (Ikeda et Matsuoka, 2002 ; Delaloye, 2004 ; Delaloye et al., 2010) : 1) les glaciers rocheux actifs qui contiennent de la glace et présentent des mouvements principalement dans le sens de plus grande pente (quelques décimètres à 3 m/an) ; 2) les glaciers rocheux inactifs qui contiennent encore de la glace mais ne présentent plus de mouvements ; 3) les glaciers rocheux fossiles ou reliques qui ne contiennent plus de glace et ne présentent plus de mouvements ; 4) les glaciers rocheux déstabilisés qui sont des glaciers rocheux actifs ayant connu au cours de ces dernières années ou décennies une évolution rapide et une forte augmentation de leur vitesse moyenne (> 3 m/an).
2Dans le contexte de la remontée altitudinale de l’isotherme 0 °C, on s’attend à une modification des mouvements superficiels des glaciers rocheux. Quelques études ponctuelles en France tendent à montrer cette évolution (Bodin, 2007; Bodin et al., 2009 ; Krysiecki, 2009) mais l’absence de données actuelles de référence sur le degré d’activité des glaciers rocheux dans les Alpes françaises ne permet pas de quantifier ces modifications. Plusieurs auteurs se sont intéressés à la question de la répartition et de l’activité des glaciers rocheux mais leurs études ne couvrent que des régions limitées (Evin, 1987 ; Chardon, 1993, 1996) ou sont non spatialisées (Francou, 1988). Avec le développement des outils de la géomatique, des inventaires géolocalisés ont été réalisés sur divers massifs (Monnier, 2004, 2006 ; Bodin, 2007 ; Echelard, 2010). Un inventaire systématique et normalisé a été entrepris par le service RTM et devrait couvrir à terme l’ensemble des Alpes françaises (Bouvet et al., 2011). À l’exception de quelques données ponctuelles de mouvement ou de mesures géophysiques révélant la présence de glace (e.g., Evin et Assier, 1983), la distinction entre formes actives, inactives et fossiles ne se fonde, dans tous ces inventaires, que sur des critères géomorphologiques (fraîcheur des formes, absence de végétation, signes d’activité visibles…).
3À l’heure actuelle, le suivi des mouvements des glaciers rocheux en France se limite à quelques sites à dans les Alpes françaises via le réseau de mesure français du permafrost PermaFRANCE (Schoeneich et al., 2010). La méthode la plus utilisée pour suivre les déformations des glaciers rocheux est celle du GPS (Global Positioning System) différentiel (D-GPS), qui nécessite beaucoup de temps par site pour extraire des mesures de positions sur un ensemble de points. Il y a donc peu de données sur le nombre et sur la localisation des glaciers rocheux actifs dans les Alpes françaises. Il est alors intéressant d’utiliser l’imagerie satellitaire pour essayer de déterminer les glaciers rocheux actifs et déstabilisés des Alpes françaises. L’interférométrie radar différentielle (D-InSAR) utilise l’analyse de deux images radar de la même zone géographique pour détecter des variations de surfaces interprétables en termes de mouvements, voire d’établir des champs de vitesses. Cette méthode déjà éprouvée pour la détermination des vitesses de glaciers blancs (Mattar et al., 1998 ; Strozzi et al., 2002 ; Forster et al., 2003 ; Trouvé et al., 2007) est adaptée pour la détection de glaciers rocheux (Strozzi et al., 2004 ; Delaloye, 2005 ; Echelard, 2010 ; Echelard et al., 2011, 2012). En Suisse, elle a permis de compléter les données d’inventaire et de détecter plusieurs cas d’accélération voire de déstabilisation de glaciers rocheux (Delaloye et Lambiel, 2007 ; Delaloye et al., 2010). La taille des images radar ERS (European Remote Sensing) couvrant des zones de 100 km x 100 km, l’analyse des interférogrammes (images produites par le processus de l’interférométrie radar différentielle) permettrait d’obtenir des données de déplacements sur les glaciers rocheux à l’échelle de massifs entiers.
4L’objectif de cette étude est d’apporter une première contribution à l’évaluation du nombre de glaciers rocheux actifs dans les Alpes françaises, d’obtenir un maximum de données quantifiées sur les mouvements de ceux-ci et de détecter ceux dont les vitesses sont les plus rapides. Cela permettra par la suite d’évaluer leurs évolutions et les enjeux potentiels. Pour ce faire, la grande quantité d’images d’archive radar disponibles pour les satellites ERS 1 et 2 a été utilisée. Ces satellites ont acquis des données pendant 20 ans à partir de 1991 et nous permettent de disposer d’un jeu de données exceptionnel sur les Alpes françaises pour cette étude. Nous avons également comparé nos résultats avec des inventaires existants de glaciers rocheux pour vérifier, modifier et mettre à jour les données globales de mouvements des glaciers rocheux des Alpes françaises.
5Se fondant sur l’analyse de la différence de phase entre deux images radar acquises de manière aéroportée (avion ou satellite) ou par radar au sol, l’interférométrie radar permet de détecter des mouvements de l’ordre de la fraction de la longueur d’onde du signal émis (de quelques millimètres à quelques centimètres selon le type de capteur). Dans notre cas, nous avons utilisé des données satellitaires d’ERS 1 et 2 de 100 km par 100 km en bande C (λ = 5,66 cm) et dont la résolution spatiale en géométrie sol est de 25 m x 25 m. Si un changement (de relief par exemple) est survenu au sol entre le premier passage du satellite et un second passage, il est alors possible de mesurer une différence de phase entre les deux images, le temps de parcours de l’onde sera en effet très légèrement différent. Les détections de mouvements de surface entre deux acquisitions pourront être centimétriques dans la ligne de visée du satellite (Gabriel et al., 1989). Cependant, pour obtenir des interférogrammes (produit de la différence de phase de deux images de phases) fiables pour les analyses, de nombreux facteurs sont à prendre en compte. L’équation de la différence de phase (ΔΦ), composée de plusieurs termes, peut se synthétiser comme suit (Wegmuller et Strozzi, 1998) :
6ΔΦ = Φ topographie + Φ orbitale + Φ changement + Φ distorsion atmosphérique + Φ bruit (1)
7Dans cette équation, seul le terme « changement de surface » nous intéresse. Pour supprimer les autres, nous utilisons des données exogènes :
-
La composante topographique est liée à l’effet stéréographique induit par les orbites légèrement différentes du ou des satellite(s) au cours de l’acquisition des données. Cela crée des franges interférométriques représentant des lignes d’équialtitudes du relief. Cette composante est utilisée pour produire des modèles numériques de terrain à haute résolution (voir par exemple Zebker et al., 1994). Pour enlever cette composante, nous utilisons un modèle numérique de terrain de la zone d’étude (BD Alti de l’IGN, Institut Géographique National, à 25 m de résolution) pour simuler une image de phase factice de la topographie et ainsi la soustraire de l’interférogramme.
-
La composante orbitale, liée à des paramètres propres à l’orbite des satellites (position, inclinaison, etc.) est corrigée grâce aux données orbitales fournies par les agences spatiales avec les images radar.
-
Les distorsions atmosphériques (nébulosité, humidité, etc.) peuvent être la cause d’erreurs sur les interférogrammes. Cependant, ces distorsions se produisent à l’échelle de la troposphère, ce qui n’introduit pas de franges d’interférences nuisibles à l’échelle de glaciers rocheux mesurant quelques centaines de mètres (Kenyi et Kaufmann, 2003).
-
Le bruit, représentant une faible partie de l’équation, est fortement atténué par les filtres appliqués aux interférogrammes permettant d’affiner la continuité des franges interférométriques et de supprimer les pixels aberrants (filtre adaptatif proposé par R.M. Goldstein et C.L. Werner (1997).
8Théoriquement, après ces étapes, seule la composante liée aux changements de surface survenus entre les deux acquisitions radar dans la ligne de visée du satellite est présente sur l’interférogramme. Une fois cette étape réalisée l’interférogramme est reprojeté en « géométrie sol ». En effet, l’angle de visée des satellites n’est pas vertical (23° pour ERS) et les images radar sont dites en « géométrie radar ». Certaines zones de l’image sont compressées alors que d’autres sont dilatées ou non analysables (fig. 1). Pour cela, nous utilisons un modèle numérique de terrain (BD Alti de l’IGN, à 25 m de résolution), qui permet de transformer la « géométrie radar » en « géométrie sol ». Au final, la résolution des interférogrammes, pour des pixels homogènes en « géométrie sol », est de 25 m de côté par pixel pour les données provenant des satellites ERS. Ce changement de géométrie est effectué afin que les images soient intégrées puis analysées dans un SIG (système d’information géographique).
Fig. 1 – Distorsions géométriques observées sur les images de RSO en géométrie radar
Fig. 1 – Geometrics distortions observed on SAR data in radar geometry
La distance entre le satellite et la cible étant de plusieurs centaines de kilomètres, les arcs de cercle correspondant aux ondes émises par le radar peuvent être assimilés à leur tangente (traits pointillés).
As the distance between the satellite and the target is several hundred kilometres, arcs corresponding to waves emitted by the radar beam can be assimilated to their tangent (dashed line).
Modifié d’après Barboux et al. (2009).
Adapted from Barboux et al. (2009).
9Le choix des images utilisées joue également un rôle important pour la qualité des interférogrammes. Plusieurs critères doivent être respectés dans la sélection des images radar pour que les traitements d’images soient efficients (cf. supra, « Méthode D-InSAR »). Le premier critère est d’ordre technique car à partir de 2001, ERS 2 a rencontré des problèmes de gyroscopes rendant inexploitables pour l’interférométrie radar différentielle toutes les données à compter de cette date. D’autres critères sont spécifiques à notre étude sur les glaciers rocheux et sont listés ci-dessous :
-
Les glaciers rocheux actifs se situant pour la plupart au-dessus de 2 500 m, nous avons choisi de ne garder que les scènes acquises entre le 15 juillet et le 15 octobre, période autour de laquelle le couvert neigeux sur les glaciers rocheux est minimal.
-
L’humidité au sol fausse le retour de l’onde émise par le satellite et empêche toute exploitation de la donnée. Pour pallier cela, les données de précipitations de MétéoFrance ont été utilisées sur sept stations à travers les Alpes françaises afin de ne sélectionner que les scènes pour lesquelles il n’y avait pas eu de précipitation le jour de l’acquisition ainsi que la veille.
-
Les satellites ont une position dans l’espace légèrement différente entre la première et la deuxième acquisition. Cette différence s’appelle la « baseline » et doit théoriquement être inférieure à 600 m dans l’optique d’une utilisation pratique (Coulson et Guignard, 1994). Elle joue un rôle important sur le nombre de franges topographiques présentes sur les interférogrammes. En pratique, la valeur tolérable de la « baseline » est liée à la précision du MNT utilisé. Plus la résolution du MNT est bonne, plus les « baselines » pourront être longues. D. Walter et W. Busch (2012) ont calculé les valeurs de « baseline » à respecter en fonction des satellites et des précisions des MNT utilisés. Suivant cette méthodologie, nous avons choisi de ne garder que les scènes ayant une « baseline » inférieure à 200 m pour minimiser les franges topographiques et la dé-corrélation spatiale entre les images.
10Pour travailler sur l’ensemble des Alpes françaises, il a été nécessaire de travailler sur des données de quatre orbites différentes, deux descendantes et deux ascendantes, elles-mêmes subdivisées en trois sections d’orbites. La répétitivité des passages d’ERS au-dessus d’une même zone ont été de 1 jour, 3 jours ou 35 jours (selon les périodes) et 1 728 scènes prises au-dessus des Alpes françaises ont été archivées. Cependant, une fois pris en compte tous les critères cités précédemment, seules 29 scènes ont été utilisées (tous passages et numéros d’orbites confondus) et ont permis de produire 11 interférogrammes analysables après traitements. Le tableau 1 récapitule les caractéristiques des interférogrammes produits. Parmi les 29 scènes que nous avons sélectionnées, les interférogrammes produits n’ont permis d’obtenir que des données interprétables pour les versants ouest des Alpes françaises (180° à 360° environ). Ces interférogrammes proviennent d’images satellitaires descendantes. Les interférogrammes provenant d’images ascendantes furent de mauvaise qualité pour deux raisons qui tiennent à l’acquisition des données : des « baselines » spatiales trop longues et/ou de mauvais recalages des deux images lors du traitement.
Tab. 1 – Caractéristiques des interférogrammes produits.
Tab. 1 – Characteristics of generated interferograms.
N : Alpes du Nord ; S : Alpes du Sud.
N: Northern Alps; S: Southern Alps.
11Après les traitements d’images vus précédemment, une analyse visuelle des interférogrammes orthorectifiés et géoréférencés a été réalisée dans un SIG. Des masques ont été créés pour analyser uniquement les zones d’intérêt pour l’étude des glaciers rocheux :
-
Masque d’altitude créé à partir de la BD Alti à 25 m de résolution de l’IGN. Toutes les zones se situant en dessous de 2 000 m ont été écartées de l’analyse. Bien que les glaciers rocheux actifs se situent aux alentours de 2 500 m, nous avons préféré garder une marge d’ajustement car certains microclimats locaux rendent possibles la présence de permafrost sporadique, donc potentiellement des glaciers rocheux, en-dessous de cette valeur moyenne de 2 500 m.
-
Masque des zones analysables. L’acquisition des images radar d’ERS se fait à 23° par rapport à la verticale, inférant des zones de repliement, d’inversion des pixels et d’ombre du relief, ce qui limitent considérablement l’analyse de ces zones (Maître, 2001). Ces zones ont été supprimées des interférogrammes.
-
Masque ne faisant apparaître que les zones potentielles de pergélisol dans les Alpes françaises.
-
L’analyse visuelle des interférogrammes nous apporte trois types d’informations : 1) l’absence de variation de phase sur l’interférogramme signifie que l’on ne peut pas détecter de mouvements à ce pas de temps (fig. 2A) ; 2) la présence d’une variation de phase sur l’interférogramme avec un motif de frange d’interférence signifie une détection de mouvement (fig. 2B) ; 3) la présence de variations de phases aléatoires au sein d’une zone localisée de l’image signifie la détection d’un mouvement plus rapide que les capacités de l’interférogramme (demi-longueur d’onde/pas de temps) ou une erreur liée par exemple à un résidu de neige (fig. 2C). La figure 2 illustre les types d’informations évoqués ci-dessus.
12Pour l’étude, deux géomorphologues ont séparément analysé les données dans le but de limiter le plus possible les erreurs et/ou les omissions humaines. Cette analyse des interférogrammes s’est faite sur SIG et a combiné les données radar, les masques évoqués précédemment et des données optiques (orthophotographies). Les analyses des deux géomorphologues ont ensuite été fusionnées afin d’établir un inventaire des glaciers rocheux actifs détectés par D-InSAR. Les indices de mouvement détectés hors glacier rocheux n’ont pas été pris en compte.
13Enfin, une couche vecteur a été créée dans laquelle figure un point pour chaque glacier rocheux détecté à partir des interférogrammes. Pour chaque point, nous avons implémenté une base de données dans laquelle figurent les pas de temps des interférogrammes sur lesquels le glacier rocheux a été détecté. Cela a permis d’établir un classement des vitesses relatives des glaciers rocheux détectés. Par exemple, un glacier rocheux qui présente des franges sur un interférogramme à 3 jours d’intervalle et qui est dé-corrélé sur un interférogramme de 6 jours fluera plus rapidement qu’un glacier rocheux qui présente des franges sur un interférogramme à 6 jours d’intervalle et qui est dé-corrélé sur un interférogramme à 9 jours et ainsi de suite. Cela a abouti à la classification suivante : formes fluant à des vitesses dont l’ordre de grandeur est le cm/jour, le cm/semaine, le cm/2 semaines ou le cm/mois.
Fig. 2 – Illustration des différents types d’informations présents sur les interférogrammes.
Fig. 2 – Illustration of different types of information present in interferograms.
Les lettres se rapportent aux explications fournies dans la partie « Méthodes d’analyses des interférogrammes ».
The letters refer to the explanations in paragraph “Méthodes d’analyses des interférogrammes”.
14Pour évaluer l’efficacité et la pertinence des détections sur les interférogrammes, nous avons comparé nos résultats avec deux inventaires de glaciers rocheux réalisés dans le massif du Queyras (Echelard, 2010) et la haute vallée de l’Arc (Monnier, 2004, 2006) à l’aide de données topographiques et de visites sur le terrain. Ces inventaires font figurer dans leur base de données une variable sur l’aspect des glaciers rocheux basé sur l’observation d’indices géomorphologiques (altitude du front et de la racine ; raideur du front ; présence ou non de végétation ; analyse diachronique de deux orthophotographies) permettant de présumer de l’état d’un glacier rocheux (Barsch, 1996 ; Ikeda et Matsuoka, 2002). Cependant, les catégories proposées pour la variable d’aspect dans ces deux inventaires sont différentes. Pour celui du Queyras, l’auteur a distingué trois catégories : les glaciers rocheux actifs, les glaciers rocheux inactifs et les glaciers rocheux fossiles. En Vanoise, l’auteur a distingué deux catégories : les glaciers rocheux frais et les glaciers rocheux dégradés. Ces derniers sont considérés comme l’équivalent des glaciers rocheux fossiles et les glaciers rocheux frais comme l’équivalent des formes actives ou inactives (Monnier, 2006). Nous verrons par la suite que cette différence terminologique donne lieu à des différences de résultats lors de la comparaison entre inventaires et interférogrammes. Ainsi, en comparant les glaciers rocheux détectés avec les données D-InSAR et les glaciers rocheux des inventaires, nous avons pu chiffrer le pourcentage de détection concordantes et le nombre de formes non détectées qui auraient pu l’être. Cela a également permis de repérer d’éventuelles erreurs d’interprétation (oublis ou classe d’activité erronée) dans les inventaires.
15La double analyse des interférogrammes a abouti à une base de données de 152 glaciers rocheux actifs sur l’ensemble des Alpes françaises. La figure 3 présente quelques exemples de détection. La grande majorité de ces détections de mouvements (en %) se trouve sur des versants orientés de 180° à 360°. La figure 4, qui présente la répartition des détections dans les Alpes, montre la quasi absence de détection dans le département de Haute-Savoie. Cela est dû à l’absence de données exploitables dans ce département. Le classement des vitesses relatives décrit au paragraphe « Méthodes d’analyses des interférogrammes » a permis d’obtenir la répartition suivante des glaciers rocheux actifs détectés : les glaciers rocheux fluant à des vitesses dont l’ordre de grandeur est le cm/jour représentent 13,4 %, le cm/semaine 19,4 %, le cm/2 semaines 27,6 % et le cm/mois 39,6 %. Les résultats trouvés sont cohérents puisque les glaciers rocheux fluent généralement à des vitesses comprises entre quelques décimètres par an et plusieurs mètres par an dans des cas exceptionnels (Delaloye et al., 2010). Sur toutes ces formes détectées, nous avons identifié un glacier rocheux qui semble déstabilisé dans le massif de la Vanoise (fig. 3C). Il présente des fractures en son centre, attestant des dynamiques différentes entre l’amont et l’aval ; d’après l’interférogramme à un jour d’intervalle, sa vitesse de fluage est estimée entre 3 et 5 m/an dans sa partie aval.
Fig. 3 – Exemples de détections de mouvements de glaciers rocheux sur les interférogrammes générés.
Fig. 3 – Examples of rockglaciers movements detections in interferograms.
A : Glacier rocheux de Panestrel, massif du Queyras. B : Glacier rocheux du Brec de Chambeyron, massif de l’Ubaye. C : Glacier rocheux de Pierre Brune nord-ouest, massif de la Vanoise. Bt indique la baseline temporelle de l’interférogramme ; les traits-tillés délimitent les glaciers rocheux. Les orthophotographies datent respectivement de 2003, 2009 et 2006.
A: Panestrel rockglacier, Queyras massif. B: Brec de Chambeyron rockglacier, Ubaye massif. C: Pierre Brune nord-ouest rockglacier, Vanoise massif. Bt indicates the temporal baseline of interferogram; dotted lines mark out rockglaciers. Orthophotography from IGN (Institut Géographique National) from respectively 2003, 2009 and 2006.
Fig. 4 – Carte de la répartition des glaciers rocheux détectés par D-InSAR dans les Alpes françaises.
Fig. 4 – Map of detected rockglaciers by D-InSAR over the French Alps.
Chaque point représente un glacier rocheux détecté. 1 : glaciers rocheux détectés par les données D-InSAR ; 2 : localisation des exemples présentés dans la fig. 3.
Each point representing a rockglacier detected. 1: rockglaciers detected by D-InSAR data; 2: location of the examples in fig. 3.
16La comparaison avec les inventaires de glaciers rocheux n’a pris en compte que les formes potentiellement détectables par les données radar. Comme nous ne disposions que d’interférogrammes descendants, seuls les glaciers rocheux dont l’azimut moyen se situe entre 180° et 360° ont été pris en compte. Compte tenu de la résolution des images (pixel de 25 m), seuls les glaciers rocheux ayant au moins la longueur ou la largeur supérieure à 300 m (soit 12 pixels de long ou de large) ont été retenus. En deçà de cette taille, nous avons estimé que trop peu de pixels définiraient le glacier rocheux pour permettre une analyse pertinente.
17Sur 18 glaciers rocheux actifs inventoriés potentiellement détectables, 12 l’ont été sur les interférogrammes, soit 2/3 de détections concordantes. La figure 5 présente un exemple de la comparaison entre l’inventaire et les interférogrammes pour le secteur du lac Saint-Anne dans le Queyras. À l’inverse, sur les 24 glaciers rocheux localisés dans des zones analysables par D-InSAR mais qualifiés d’inactifs ou fossiles (i.e., sans mouvement), aucun n’a été détecté sur les interférogrammes. Ces résultats confirment donc pleinement les informations provenant de l’inventaire concernant ces formes. Au total, 85,7 % des informations sur l’activité des glaciers rocheux du massif du Queyras provenant de l’inventaire ont été confirmées par les données d’interférométrie radar différentielle. Cependant, le pas de temps le plus long dont nous disposons est de 70 jours pour des mouvements minimum détectable d’environ 1 cm. On ne peut donc pas exclure que certaines formes présentant des mouvements inférieurs à 5 cm/an n’aient pas été détectées. Par ailleurs, 7 formes ont été détectées bien que non répertoriées : elles seront à ajouter à l’inventaire réalisé en 2010. Il s’agit de formes localisées dans des zones d’ombre ou sous des nuages sur les orthophotographies utilisées pour l’inventaire. De plus, la géologie du Queyras fait que certains glaciers rocheux du massif sont formés de débris assez fins composés de schistes lustrés ou de flyschs à helminthoïdes, ce qui rend parfois difficile l’interprétation géomorphologique lors de la réalisation d’un inventaire à partir d’orthophotographies. En effet, les sillons et les bourrelets caractéristiques des glaciers rocheux y sont souvent moins francs.
18Sur 34 glaciers rocheux frais inventoriés potentiellement détectables, 17 l’ont été sur les interférogrammes, soit 50 % de détections concordantes. À l’inverse, sur les 45 glaciers rocheux localisés dans des zones analysables par D-InSAR mais qualifiés de dégradés (i.e., sans mouvement), un seul l’a été sur les interférogrammes. Ces résultats confirment presque la totalité des informations provenant de l’inventaire concernant l’absence de mouvement de ces formes. Comme pour l’inventaire dans le Queyras, on ne peut exclure que certaines formes présentant des mouvements inférieurs à 5 cm/an n’aient pas été détectées. Au total, 77 % des informations sur l’activité des glaciers rocheux de la haute vallée de l’Arc provenant de l’inventaire ont été confirmées par les données d’interférométrie radar différentielle. Par ailleurs, 36 glaciers rocheux ont été détectés bien que non répertoriées. Deux remarques peuvent être émises à partir de ces résultats : 1) la concordance entre l’inventaire et les données interférométriques n’est pas satisfaisante ; 2) un nombre de glaciers rocheux supérieur à l’échantillon de départ a été détecté (échantillon de départ : 34 ; détection non présente dans l’échantillon : 36). Cela s’explique en grande partie par des différences de définitions entre les auteurs, qui vont être discutées.
Fig. 5 – Comparaison entre glaciers rocheux actifs de l’inventaire et glaciers rocheux détectés par D-InSAR.
Fig. 5 – Comparison between active rockglaciers from the inventory and rockglaciers detected by D-InSAR.
1 : détectables par D-InSAR ; 2 : non détectables par D-InSAR ; 3 : glaciers rocheux détectés par d-InSAR.
1: detectable by D-InSAR; 2: undetectable by D-InSAR; 3: rockglaciers detected by D-InSAR.
19Les résultats obtenus à l’échelle des Alpes françaises ont permis d’obtenir une première carte faisant figurer 152 glaciers rocheux actifs répartis dans les Alpes françaises alors que jusqu’à présent, d’après le réseau de suivi national du permafrost en France (PermaFRANCE), seuls 7 glaciers rocheux disposaient de données de mouvements (par GPS ou théodolite ; Schoeneich et al., 2010). Les autres formes des Alpes françaises, qualifiées d’actives dans les inventaires et la littérature scientifique, le sont uniquement sur la base de données et d’observations indirectes. Cette étude a donc permis de multiplier par plus de 20 le nombre de glaciers rocheux que l’on peut qualifier d’actifs sur la base de données de mouvements. En revanche, plusieurs limites importantes sont à préciser :
-
Le manque de données en orbites ascendantes fait que seule une moitié des versants a été analysée (de 180° à 360°) ;
-
Les données de mouvements de surface obtenues sont des ordres de grandeurs relatifs et non des valeurs absolues ;
-
La différentiation entre mouvements latéraux et verticaux est délicate à estimer pour de très faible mouvements ;
-
Les images d’archives d’ERS utilisables pour de l’interférométrie radar différentielle datent des années 1990. Les résultats obtenus sont donc des données de mouvements des glaciers rocheux au cours de cette décennie. Cependant, on constate que la température du pergélisol, qui régit directement la vitesse des glaciers rocheux à l’échelle décennale (Delaloye et al., 2010 ; Kääb et al., 2007), a peu évolué en quelques décennies (Delaloye et al., 2010). On peut donc admettre que la très large majorité des résultats obtenus est toujours valable à l’heure de l’écriture de cet article.
-
La résolution au sol des images radar ERS (25 m) ne permet pas d’obtenir de robustes résultats pour certains glaciers rocheux actifs de petite taille (< 100 m voire 150 m).
20En ce qui concerne la comparaison avec les inventaires de glaciers rocheux, 77 % des données de l’inventaire de la haute vallée de l’Arc et 85,7 % de celles de l’inventaire du Queyras ont été confirmés par l’imagerie radar (pour les glaciers rocheux de plus de 300 m situés dans les zones visibles des interférogrammes). La question est alors de savoir quelles sont les causes des cas non confirmés. Rappelons que les 14,3 % pour le Queyras et les 33 % pour la Vanoise/Haute-Maurienne qui restent non confirmés par les données radar correspondent, respectivement, à 6 glaciers rocheux qualifiés d’actifs et à 17 qualifiés de frais dans les inventaires mais non détectés sur les interférogrammes. Nous proposons plusieurs hypothèses, communes pour les deux inventaires, pour expliquer cela :
-
Des artefacts atmosphériques ont rendu la détection impossible sur ses quelques formes ;
-
Les filtres utilisés pendant le traitement des données dans le but de réduire le bruit présent sur les interférogrammes (filtre adaptatif proposé par R.M. Goldstein et C.L. Werner, 1997) ont fait disparaître certaines informations de mouvements de petite ampleur ;
-
Ces glaciers rocheux ne fluaient pas suffisamment rapidement dans les années 1990 pour être détectés par des satellites ayant les caractéristiques techniques du type de ceux d’ERS 1 et 2.
21En ce qui concerne cette dernière hypothèse, elle semble peu probable au vu d’une récente publication (Delaloye et al., 2010) qui tend à montrer que la dynamique des glaciers rocheux des Alpes suisses n’a guère évolué entre les années 1990 et la situation moyenne actuelle. L’hypothèse selon laquelle ces erreurs seraient liées à des artefacts semble possible mais limitée car les artefacts atmosphériques laissent apparaître des taches caractéristiques sur les interférogrammes ; taches que nous avons pris le soin de vérifier avant d’analyser les interférogrammes. Quant aux filtres utilisés pendant le traitement des données, il est possible qu’ils aient fait disparaître une certaine quantité d’informations présentes originellement sur les données, donc potentiellement tout ou partie des 6 glaciers rocheux en question. Malgré tout, cet impact semble également limité car nous avons pris soin de ne pas considérer les glaciers rocheux de petite taille pour, justement, parer au mieux à cette éventualité.
22Pour l’inventaire dans le massif du Queyras, l’hypothèse qui nous semble la plus probable est celle qui suppose des erreurs d’interprétation lors de la création de l’inventaire des glaciers rocheux. Etant donné le pourcentage significatif d’informations concordantes entre l’inventaire et les données D-InSAR (85,7 %) ainsi que la difficulté de faire la différence entre un glacier rocheux actif ou inactif sans données directes de déplacement, il semble probable que la majorité des 6 cas particuliers qualifiés d’actifs dans l’inventaire soient en réalité inactifs.
23Pour l’inventaire dans la haute vallée de l’Arc, l’auteur a choisi de ne pas faire de distinction entre « glaciers rocheux actifs » et « glaciers rocheux inactifs », en les regroupant dans une catégorie « glaciers rocheux frais ». Il est vrai que sans donnée de mouvement, il est souvent difficile de faire la différence. Par conséquence, cela explique une grande partie des 17 glaciers rocheux frais de l’inventaire non détectés sur les interférogrammes. Il est fort probable que ces formes soient inactives et donc logiquement non détectées sur les interférogrammes.
24C’est également une différence de définition qui est à l’origine des 36 glaciers rocheux détectés sur les interférogrammes mais non répertoriés dans l’inventaire. Dans la littérature scientifique, certains auteurs estiment que la présence de bourrelets et de sillons est essentielle à l’identification d’un glacier rocheux (Capps, 1910 ; Vitek et Giardino, 1987, in Monnier, 2004) alors que pour d’autres, elle n’est pas nécessairement requise (Evin, 1987 ; Delaloye, 2004). Dans l’inventaire de la haute vallée de l’Arc, l’auteur prend comme élément essentiel de l’identification des glaciers rocheux la présence de bourrelets et de sillons, contrairement à notre étude durant l’analyse des interférogrammes. Cela explique les 36 formes détectées en plus de celles de l’inventaire qui correspondent bien, après vérification, à des glaciers rocheux dont la morphologie de surface ne présente pas, ou peu, de bourrelets et de sillons marqués.
25Les résultats obtenus au regard des objectifs de départ sont satisfaisants mais incomplets. Concernant le premier objectif, celui de l’évaluation du nombre de glaciers rocheux actifs dans les Alpes françaises, les interférogrammes ont permis de dresser une carte en faisant figurer 152 glaciers rocheux. Parmi eux, un cas anormalement rapide a été détecté dans le massif de la Vanoise (cf. supra « Résultats », fig. 3C) et sera suivi de près dans les années à venir par le Parc National de la Vanoise ainsi que dans le cadre du réseau de mesure français du permafrost (PermaFRANCE). Ces informations pourront être utiles aux services de gestion territoriale française pour l’évaluation des risques. Cependant, nous n’avons pu exploiter que des données provenant d’orbites descendantes, ce qui a limité les zones de détection. Le second objectif était d’obtenir des données quantifiées de mouvement, mais la volonté d’obtenir des résultats à l’échelle des Alpes françaises n’a pas été compatible avec le temps et la minutie nécessaire au déroulement de la phase des interférogrammes produits avec les données ERS (dans des zones de relief escarpée, les franges d’interférences sont trop souvent tronquées pour être déroulées de manière automatique).
26Cette étude a été réalisée avec des images ERS car l’accès aux données d’archives de ces satellites est aisé. Cela a permis de faire une première évaluation des zones les plus pertinentes à étudier et à suivre dans les Alpes françaises en termes de glaciers rocheux actifs. À partir de ces résultats, d’autres projets sont à envisager en utilisant des satellites de générations plus récentes pour des données mieux calibrées et des résultats de meilleure résolution. Bien que pertinente, cette méthode nécessite tout de même la présence d’utilisateurs compétents pour analyser les interférogrammes. L’automatisation des détections via des filtres de zonages et des recherches de « motifs » (ou patterns) dans les images permettraient un suivi simplifié et efficace des glaciers rocheux par des personnes spécialistes ou pas de l’imagerie radar et de la géomorphologie périglaciaire.