L'objectif de ce mémoire est de développer et d'étudier un nouveau type de filtre particulaire appelé le filtre de Kalman-particulaire à noyaux (KPKF). Le KPKF modélise la densité conditionnelle de l'état comme un mélange de gaussiennes centrées sur les particules et ayant des matrices de covariance de norme petite. L'algorithme du KPKF utilise une correction de l'état de type Kalman et une correction de type particulaire modifiant les poids des particules. Un nouveau type de ré-échantillonnage permet de préserver la structure de cette mixture. Le KPKF combine les avantages du filtre particulaire régularisé en terme de robustesse et du filtre de Kalman étendu en terme de précision. Cette nouvelle méthode de filtrage a été appliquée au recalage de navigation inertielle d'un aéronef disposant d'un radio altimètre. Les résultats obtenus montrent que le KPKF fonctionne pour de grandes zones d'incertitude initiale de position.

Ce mémoire de thèse est constitué d'une thématique régionale et d'une thématique générale. La première thématique correspond à une étude du Rift Ouest-Européen (ROE) par l'analyse du Rift du Massif Central (RMC). Dans le Massif Central (MC), la répartition du volcanisme associé à l'extension cénozoïque permet de différencier trois épisodes: le Volcanisme Pré-rift, le Volcanisme Syn-rift et la Phase Volcanique Majeure. Le Volvanisme Pré-rift (Paléo-Eocène) est associé à une surrection du MC. Le Volcanisme Syn-rift est lié à la période d'extension Eo-Miocène à l'origine des grabens du MC. Durant cette période, la chronologie des événements plaide en faveur d'une évolution de type rift passif. La Phase Volcanique Majeure, responsable des principales provinces volcaniques, est caractérisée par deux pics de magmatisme contem.porains de périodes de surrection sans création de graben et sans sédimentation. Cette évolution est compatible avec le stade initial d'un rift actif. Cette évolution du RMC est commune à l'ensemble du ROE et la répartition du ROE concentrique' autour des Alpes suggère un lien étroit entre ces deux systèmes géologiques. La phase éo-alpine (== 65 Ma) serait responsable du flambage litho sphérique en avant des Alpes et du Volcanisme Pré-rift. A l'Eocène, la création de la racine lithosphérique alpine permet d'expliquer le rift passif contemporain au niveau du ROE, puis le "rift actif' à partir du Miocène. La thématique générale a pour but d'étudier la dynamique de l'extension par une approche analogique. Les expériences à l'échelle crustale montrent le rôle majeur de la vitesse d'extension sur les structures formées. Dimensionnées à partir du rapport de résistance entre les parties fragile et ductile, les modèles montrent que la géométrie des structures varie en fonction du nombre de discontinuités de vitesse et de la vitesse d'extension. Appliqués à la nature, ces modèles permettent d'expliquer la différence de géométrie entre le graben du Rhin et le RMC.

Ce travail vise à étendre les méthodes d'analyse d'images satellitaires de crues au-delà de la détection des limites d'inondation afin d'estimer des niveaux d'eau distribués dans l'espace et d'aider la modélisation hydraulique.Inspirée des travaux de Raclot (2003) sur photographies aériennes qui fournissent des incertitudes moyennes de ±20cm, la méthode d'estimation des niveaux d'eau utilise des images satellitaires RADAR de crue et un MNT fin. Elle repose sur 1) une phase de télédétection pour cartographier l'inondation et analyser la pertinence de ses limites pour l'estimation des niveaux d'eau, 2) une phase d'analyse spatiale dans laquelle les niveaux d'eau sont estimés par croisement entre les limites pertinentes et un MNT fin, puis contraints par le schéma de circulation des eaux. Les estimations de niveaux d'eau obtenues ont une incertitude moyenne de ±38cm pour une image RADARSAT-1 d'une crue de la Moselle (France, 1997). Des travaux de validation ont permis de calculer une RMSE de l'ordre de 13cm pour une image ENVISAT de l'Alzette (Luxembourg, 2003).Pour aider la modélisation hydraulique, la démarche proposée vise à réduire le phénomène d'équifinalité grâce aux images satellitaires. Pour cela, un calage "traditionnel" à partir d'hydrogrammes observés est complété par la comparaison entre résultats du modèle et surfaces inondées ou niveaux d'eau extraits des images. Pour cerner les incertitudes du calage, des simulations Monte-Carlo ont été mises en place. En perspective, la prévision de l'évolution d'une crue après acquisition d'une image devrait bénéficier de modèles mieux contraints grâce à l'utilisation des images comme condition initiale ou donnée de calage.

Ce travail porte sur la morphodynamique des bancs de galets alternés dans une rivière de montagne aménagée. Le site expérimental est un tronçon long de 8 km de la rivière Arc située dans les Alpes françaises. Une analyse des cartes historiques, des photos aériennes et des relevés topographiques montre que le lit de la rivière a évolué d’un style en tresse vers un lit rectiligne avec des bancs alternés contraints entre deux digues. Au milieu des années 1990, le lit de la rivière a été aplani pour former des profils trapézoïdaux. Cependant, les bancs alternés ont réapparu en moins de dixans. Une analyse des facteurs contrôlant leur formation est fournie. Les observations de terrain et l’analyse numérique 1D suggèrent que le motif de bancs alternés représente un état de quasi équilibre du tronçon. D’autre part, la topographie d’un sous-tronçon incluant un banc de galets aété mesurée 18 fois entre 2006 et 2012. Une méthodologie de construction des MNT basée sur leslignes directrices et incluant une estimation de l’erreur a été proposée pour une application aux litsde rivières. La technique de différenciation des MNT a été utilisée pour détecter les changements morphologiques du banc. Ainsi, la morphodynamique du sous-tronçon a été étudiée en la reliant aux principaux processus qui la contrôlent tels que l’hydrologie, les formes des hydrogrammes, la répartition des contraintes de cisaillement, les apports amont en sédiments et l’énergie du coursd’eau. La forme des hydrogrammes de crue ainsi que les apports en sédiments semblent être lesprincipaux facteurs influençant le bilan sédimentaire sur le site étudié. En outre, des modélisationsnumériques 2D hydro-sédimentaires ont été utilisées pour appréhender la dynamique du banc. Lessimulations ont révélé l’importance de l’interaction entre les sédiments fins et les galets sur ladynamique globale du banc. / The current research investigates the morphodynamics of gravel alternate bars in an engineeredmountainous river. The experimental site corresponds to an 8 km long reach of the Arc River, in the French Alps. An analysis of historical maps and recent aerial photographs and topographic campaigns shows that the original braided river bed has evolved to a straight bed with alternate bar system constrained within two embankments following two centuries of man-made engineeringprojects. In the mid 1990s the river bed was flattened to fit a trapezoidal cross-sectional design. However, alternate bars re-appeared in less than ten years. An analysis of the controlling factors of bar formation and evolution is provided. Field observation and 1D numerical analysis suggest that the alternate-bar pattern represents a quasi-equilibrium state of the river bed. A singlegravel bar within the reach was monitored 18 times from 2006 to 2012. A methodology of DTM construction based on breaklines is proposed for river beds and an error estimation is presented. DTM differencing technique was used to detect morphologic changes of the gravel bar. Therefore, the gravel bar and main channel dynamics were analyzed relative to main controlling processes suchas extreme hydrologic events, sediment supply, hydrograph shape, bed shear stress distribution andriver energy. Sediment supply and hydrograph shape appeared to be the main factors influencing the study site sediment budget. Moreover, 2D numerical modelling of water flow and sediment transport was performed. Sensitivity analysis was conducted to yield results in agreement with observations. Simulations revealed the relevance of the interaction between fine sediments and gravels on the bar dynamics.

Le domaine maritime littoral constitue un milieu fragile, placé à l'interface Terre/Mer où s'exercent de nombreuses pressions issues des milieux terrestres, maritimes et atmosphériques. Les approches globales de Gestion Intégrée des Zones Côtières (GIZC) ont pour but de pallier les lacunes des politiques précédentes au niveau des interactions simultanées des différents acteurs du littoral. Cette approche requiert un Modèle Numérique de Terrain (MNT) littoral continu ainsi que des données maritimes et terrestres cohérentes entre elles. En raison d'acquisitions ou de productions différentes pour les MNT terrestres et maritimes, l'ensemble du relief n'est pas traité de manière ni homogène ni continue. De plus les techniques d'intégration des MNT actuelles entraînent des lissages des reliefs, voire la perte des éléments caractéristiques. Des outils et une nouvelle méthodologie d'intégration de MNT ont donc été définis. Les outils utilisent la distance de Fréchet pour calculer une distance maximale d'écartement entre les points homologues de deux lignes. La valeur de cette distance donne une mesure de ressemblance de forme qui sert de base à leur appariement ultérieur. La méthodologie de fusion proposée dans cette thèse se veut générique tout en s'appuyant sur les lignes caractéristiques spécifiques du paysage. Elle effectue un enchaînement de traitements répartis en trois étapes majeures : 1) transformations cartographiques et analyse des zones de recouvrement ; 2) segmentation et extraction de lignes caractéristiques ; 3) appariement et fusion à l'aide de déformations élastiques. Ces outils d'intégration ont été testés sur des données réelles en 2D pour la comparaison du trait de côte, puis validés sur des MNT simulés 3D.

Les modèles numériques de terrain, cas particulier de modèles numériques de surfaces, n'ont pas la même erreur quadratique moyenne en planimétrie qu'en altimétrie. Différentes solutions ont été envisagées pour déterminer séparément l'erreur en altimétrie et l'erreur planimétrique, disposant, bien entendu, d'un modèle numérique plus précis comme référence. La démarche envisagée consiste à déterminer les paramètres des ellipsoïdes d'erreur, centrées dans la surface de référence. Dans un premier temps, l'étude a été limitée aux profils de référence avec l'ellipse d'erreur correspondante. Les paramètres de cette ellipse sont déterminés à partir des distances qui séparent les tangentes à l'ellipse du centre de cette même ellipse. Remarquons que cette distance est la moyenne quadratique des distances qui séparent le profil de référence des points du modèle numérique à évaluer, c'est à dire la racine de la variance marginale dans la direction normale à la tangente. Nous généralisons à l'ellipsoïde de révolution. C'est le cas ou l'erreur planimétrique est la même dans toutes les directions du plan horizontal (ce n'est pas le cas des MNT obtenus, par exemple, par interférométrie radar). Dans ce cas nous montrons que le problème de simulation se réduit à l'ellipse génératrice et la pente du profil correspondant à la droite de pente maximale du plan appartenant à la surface de référence. Finalement, pour évaluer les trois paramètres d'un ellipsoïde, cas où les erreurs dans les directions des trois axes sont différentes (MNT obtenus par Interférométrie SAR), la quantité des points nécessaires pour la simulation doit être importante et la surface tr ès accidentée. Le cas échéant, il est difficile d'estimer les erreurs en x et en y. Néanmoins, nous avons remarqué, qu'il s'agisse de l'ellipsoïde de révolution ou non, que dans tous les cas, l'estimation de l'erreur en z (altimétrie) donne des résultats tout à fait satisfaisants. / A Digital Surface Model (DSM) is a numerical surface model which is formed by a set of points, arranged as a grid, to study some physical surface, Digital Elevation Models (DEM), or other possible applications, such as a face, or some anatomical organ, etc. The study of the precision of these models, which is of particular interest for DEMs, has been the object of several studies in the last decades. The measurement of the precision of a DSM model, in relation to another model of the same physical surface, consists in estimating the expectancy of the squares of differences between pairs of points, called homologous points, one in each model which corresponds to the same feature of the physical surface. But these pairs are not easily discernable, the grids may not be coincident, and the differences between the homologous points, corresponding to benchmarks in the physical surface, might be subject to special conditions such as more careful measurements than on ordinary points, which imply a different precision. The generally used procedure to avoid these inconveniences has been to use the squares of vertical distances between the models, which only address the vertical component of the error, thus giving a biased estimate when the surface is not horizontal. The Perpendicular Distance Evaluation Method (PDEM) which avoids this bias, provides estimates for vertical and horizontal components of errors, and is thus a useful tool for detection of discrepancies in Digital Surface Models (DSM) like DEMs. The solution includes a special reference to the simplification which arises when the error does not vary in all horizontal directions. The PDEM is also assessed with DEM's obtained by means of the Interferometry SAR Technique

Les Systèmes d'Information Géographique (SIG) gèrent des informations complexes localisées dans un espace géographique. Durant la dernière décennie, ces systèmes ont vu l'émergence d'une part des SIG en trois dimensions et plus particulièrement des modèles numériques de terrain (MNT) et d'autre part de l'intégration de données temporelles. Les travaux de recherche en informatique présentés dans ce manuscrit s'attaquent à ces deux thématiques et sont appliqués au domaine maritime. En termes de MNT, ces recherches ont essentiellement porté sur l'intégration de MNT terrestre et maritime. Une méthode de fusion générant un MNT continu a été définie. Elle est basée sur l'appariement d'éléments caractéristiques linéaires (crêtes, talwegs, lignes de rupture de pente). Elle emploie des distances linéaires discrètes introduites par Fréchet. Ces distances fournissent des vecteurs d'appariement autorisant la fusion de MNT à l'aide d'une technique de déformation élastique. En ce qui concerne les SIG spatio-temporels, ces travaux concernent les déplacements d'objets mobiles et plus particulièrement les navires. Dans le cadre de la navigation maritime, les SIG ont une place prépondérante via les aides à la navigation, les systèmes de surveillance du trafic et les simulateurs de formation. Dans ce cadre, un nouveau modèle de représentation des déplacements, la vue relative, a été formalisé. Ce formalisme s'appuie sur la perception humaine de l'espace spatio-temporel et des objets mobiles. Cette représentation est complémentaire de la représentation absolue. Les déplacements relatifs vis-à-vis d'un objet référent sont ainsi plus facilement analysables. Les évolutions des distances et des vitesses relatives sont mieux perçues. Parallèlement, une simulation des activités de navires à l'aide de systèmes multi-agents a été conçue. Les déplacements simultanés d'un grand nombre de navires dans un environnement maritime sont ainsi simulés. Elle prend en compte les différentes activités concurrentes des objets mobiles et les contraintes géographiques. Qui plus est, elle intègre des raisonnements à base de patrons pour reproduire les raisonnements humains. Finalement, à partir des bases de données historiques des positions d'objets mobiles dans un environnement ouvert, des techniques d'extractions des trajectoires et de fouille de données spatio-temporelles ont été réalisées. Elles permettent pour chaque itinéraire et chaque type d'objets de définir des routes types spatio-temporelles et de détecter à la volée des comportements inhabituels. Ces travaux, bien qu'appliqués au domaine maritime, sont génériques. Les méthodes de fusion de MNT peuvent être appliquées à tous types de MNT ayant des zones de recouvrement. De même, les travaux de modélisation, d'analyses et de simulation sur les objets mobiles sont adaptables à tous les déplacements d'objets mobiles dans des espaces ouverts.

L'objectif de cette thèse est de contraindre la nature pétrologique et géométrique de la partie profonde d'un domaine orogénique en convergence. Dans ce but, nous avons créé dans un système d'information géographique projeté (Madagascar Laborde 1923), une base de données géologique du tiers Sud de l'île de Madagascar. Les données utilisées sont variées : observations de terrain, images satellites, données topographiques, mesures gravimétriques, radiométriques et aéromagnétiques. L'analyse, le traitement, la confrontation et l'intégration de ces différentes informations (rendues compatibles) nous a permis de réaliser avec précision de nouvelles cartes géologiques de la région, ainsi que des représentations cartographiques nouvelles. Nous mettons notamment en évidence des zones de cisaillements verticales (ZC) majeures et mineures, des fractures, des dykes et des limites lithologiques.<br />Nous précisons la nature des alignements sur l'axe Bongolava – Ranotsara : (i) La structure ductile de direction N 140°, cartographiée dans plusieurs littératures, n'est qu'une impression régionale (à petite échelle) due à la déflexion N 157° des ZC de direction Nord – Sud. La faille de Ranotsara (de direction N 140°) n'existe pas en termes de faille régionale, elle est une impression (à petite échelle) de la restitution cartographique causée par le relais de failles de taille ≤60 km de long et qui sont orientées N 115°, N 140° et N 160°. <br />Nous identifions que les ZC sont les lieux de localisation d'anomalie en éléments incompatibles (U-Th). Elles sont également le lieu de concentration d'oxydes de magnétite et titanomagnétite, associé à la disparition de la biotite. Ces observations soulignent la concentration des magmas et de la fusion au sein des ZC. Ainsi, les ZC sont le lieu de concentration de magmas qui permettent l'advection de chaleur. En effet, dans le Sud de Madagascar les conditions de haute température (800°C) sont homogénéisées sur de grande surface. <br />Nous établissons une nouvelle carte des ensembles lithologiques qui met en évidence des contrastes de compétence: les ZC se développent préférentiellement dans l'ensemble d'orthogneiss peu alumineux; les ensembles de métasédiments alumineux et les gneiss migmatitiques à amphibole sont plus résistants et forment des boudins à différentes échelles ; l'ensemble de granites et syénites stratoïdes et gneiss migmatitiques à amphibole constituent les blocs les plus rigides. Nous identifions des dykes kilométriques de basaltes, de dolérites et de microgranites qui sont orientés préférentiellement sur les directions N 25° et N 120°.<br />Par croisement des données, nous montrons que la position des minéralisations économiques est contrôlée par les anisotropies kilométriques (ZC, charnières de plis, limites lithologiques) qui localisent les circulations fluides et magmatiques. Les minéralisations primaires (métamorphiques ou métasomatiques : corindons, fer et graphite) se localisent au sein des ZC. Les minéralisations secondaires (uranothorianite, or) s'expliquent par des remobilisations au sein des ZC. La position finale des minéralisations s'explique cependant par la géométrie des réseaux hydrographiques et des failles. A l'échelle du Gondwana, la géométrie des ZC majeures néoprotérozoïques contrôle le développement des systèmes de failles et de dykes (170Ma 90Ma) orientés N 160°. <br />Dans 4 régions : Satrokala, Fenoarivo, Analavoka, Lavaraty ; existent des anomalies magnétiques et gravimétriques de type «bulls eyes» qui se localisent à l'intersection des dykes N 160° et des dykes N 30°, N 140° et N 170° ; et parfois avec les ZC. Ces points correspondent à des creux topographiques circulaires (10 m - 2,5 km de diamètre) isolés ou associés en champ de point. Ces zones sont potentiellement des intrusions kimberlitiques. La présence de ces pipes à l'intersection d'anisotropies kilométriques est favorable à l'ascension rapide des magmas silicatés, soit la présence de kimberlites diamantifères.

La production de Modèles Numériques de Terrain (MNT) a subi d’importantes évolutions durant les deux dernières décennies en réponse à une demande croissante pour des besoins scientifiques et industriels. De nombreux satellites d’observation de la Terre, utilisant des capteurs tant optiques que radar, ont permis de produire des MNT couvrant la plupart de la surface terrestre. De plus, les algorithmes de traitement d’images et de nuages de points ont subi d’importants développements. Ces évolutions ont fourni des MNT à différentes échelles pour tout utilisateur. Les applications basées sur la géomorphologie ont profité de ces progrès. En effet, ces applications exploitent les formes du terrain dont le MNT constitue une donnée de base. Cette étude a pour objectif d’évaluer l’impact des paramètres de production de MNT par photogrammétrie et par InSAR sur la qualité de position et de forme de ces modèles. La qualité de position, évaluée par les producteurs de MNT, n’est pas suffisante pour évaluer la qualité des formes. Ainsi, nous avons décrit les méthodes d’évaluation de la qualité de position et de forme et la différence entre elles. Une méthode originale de validation interne, qui n’exige pas de données de référence, a été proposée. Ensuite, l’impact des paramètres de l’appariement stéréoscopique, du traitement interférométrique ainsi que du rééchantillonnage, sur l’altitude et les formes, a été évalué. Finalement, nous avons conclu sur des recommandations pour choisir correctement les paramètres de production, en particulier en photogrammétrie.Nous avons observé un impact négligeable de la plupart des paramètres sur l’altitude, à l’exception de ceux de l’InSAR. Par contre, un impact significatif existe sur les dérivées de l’altitude. L’impact des paramètres d’appariement présente une forte dépendance avec la morphologie du terrain et l’occupation du sol. Ainsi, le choix de ces paramètres doit être effectué en prenant en considération ces deux facteurs. L’effet des paramètres du traitement interférométrique se manifeste par des erreurs de déroulement de phase qui affectent principalement l’altitude et peu les dérivées. Les méthodes d’interpolation et la taille de maille présentent un impact faible sur l’altitude et important sur ses dérivées. En effet, leur valeur et leur qualité dépendent directement de la taille de maille. Le choix de cette taille doit s’effectuer selon les besoins de l’application visée. Enfin, nous avons conclu que ces paramètres sont interdépendants et peuvent avoir des effets similaires. Leur choix doit être effectué en prenant en considération à la fois l’application concernée, la morphologie du terrain et son occupation du sol afin de minimiser l’erreur des résultats finaux et des conclusions. / The production of Digital Elevation Models (DEMs) has undergone significant evolution duringthe last two decades resulting from a growing demand for scientific as well as industrial purposes.Many Earth observation satellites, using optical and radar sensors, have enabled the production ofDEMs covering most of the Earth’s surface. The algorithms of image and point cloud processing havealso undergone significant evolution. This progress has provided DEMs on different scales, which canfulfill the requirements of many users. The applications based on geomorphology have benefitted fromthis evolution. Indeed, these applications concentrate specifically on landforms for which the DEMconstitutes a basic data.The aim of this study is to assess the impact of the parameters of DEM production byphotogrammetry and InSAR on position and shape quality. The position quality, assessed by DEMproducers, is not sufficient for the evaluation of shape quality. Thus, the evaluation methods ofposition and shape quality and the difference between them are described. A novel method of internalvalidation, which does not require reference data, is proposed. Then, the impact of image matchingand interferometric processing parameters as well as resampling, on elevation and shapes, is assessed.Finally, we conclude on recommendations on how to choose the production parameters correctly,particularly for photogrammetry.We observe little impact from most of the parameters on the elevation, except InSAR parameters.On the other hand, there is a significant impact on the elevation derivatives. The impact of matchingparameters presents a strong dependence on the terrain morphology and the landcover. Therefore,these parameters have to be selected by taking into account these two factors. The effect ofinterferometric processing manifests by phase unwrapping errors that mainly affect the elevation andless the derivatives. The interpolation methods and the mesh size present a small impact on theelevation and a significant impact on the derivatives. Indeed, the value of the derivatives and theirquality depend directly on the mesh size. The selection of this size has to be made according to theforeseen application. Finally, we conclude that these parameters are interdependent and can havesimilar effects. They must be selected according to the foreseen application, the terrain morphologyand the landcover in order to minimize the error in the final results and the conclusions.

L'étude conjointe des parties immergée et émergée de l'Ile de La Réunion (Océan Indien), basée sur la valorisation des données marines multifaisceaux acquises ces 20 dernières années et confrontée à l'ensemble des connaissances géologiques et géophysiques, démontre l'importance considérable des processus de déstabilisations de flanc dans l'évolution des édifices de ce système volcanique océanique. Au total, une cinquantaine d'événements multi-échelles, primaires ou secondaires, d'origine aérienne, côtière ou sous-marine, ont démantelé de façon récurrente les volcans du Piton des Neiges, des Alizés et du Piton de La Fournaise au cours des 2 derniers millions d'années au moins. La plupart correspondent à des phénomènes instantanés et catastrophiques de type avalanche de débris. Leurs dépôts sont principalement accumulés en mer, au sein de quatre gigantesques excroissances topographiques, les Reliefs Sous-Marins, localisés à l'est, au nord, à l'ouest et au sud de l'île. Ces éperons triangulaires, de 20-30 km de large au niveau de la côte et 100-150 km à leur extrémité, à 70-80 km du littoral, se sont construits progressivement par la superposition et la juxtaposition des produits déplacés dans des épisodes de déstabilisation pouvant impliquer plusieurs centaines de km3 de matériaux. Une très faible proportion est retrouvée à terre, au cœur et en bordure littorale du massif du Piton des Neiges et constitue le faciès proximal des dépôts immergés. Les principales structures de la morphologie sub-aérienne de l'île sont interprétées comme (1) les zones sources de déstabilisations relativement récentes, souvent remodelées par l'érosion et/ou comblées par l'activité volcanique ultérieure, ou (2) les traces rémanentes d'amphithéâtres d'avalanches de débris anciens dont les failles bordières ont pu être réactivées dans des processus volcano-tectoniques (effondrements caldériques, étalement gravitaire, subsidence du complexe hypovolcanique du Piton des Neiges,...) et/ou lors d'épisodes ultérieurs de déstabilisation. Ces failles ont probablement guidé le démantèlement par érosion des massifs et ainsi facilité le creusement des dépressions centrales du Piton des Neiges (les trois cirques de Mafate, Salazie, Cilaos et le paléo-cirque des Marsouins) et des principales vallées de l'île. L'étude du volcan Mount Haddington (Ile de James Ross, Antarctique) et les expériences analogiques montrent que la déformation ductile d'horizons de faible résistance inclus dans la structure des volcans boucliers, et en particulier de deltas volcanoclastiques essentiellement constitués de dépôts d'avalanches de débris, est un élément clé dans la génération de glissements rotationnels et d'avalanches de débris ultérieures. Ce modèle est privilégié à La Réunion, en considérant un rôle prépondérant des processus gravitaires dans le déclenchement des déstabilisations des volcans de l'île.