L'objectif de ce travail de thèse est de comprendre le rôle de la déglaciation et l'influence de l'évolution géomorphologique de la basse vallée glaciaire de la Romanche sur l'initiation et la localisation des instabilités gravitaires de ses versants rocheux.<br />Une étude multidisciplinaire combinant des données géologiques, morphologiques, géophysiques, et chronologiques a ainsi été menée à deux échelles spatiales différentes : celle de la basse vallée de la Romanche et celle du mouvement de terrain de Séchilienne.<br />La détermination de la fracturation du massif, obtenues par morphologie de surface, indiquent que les caractéristiques lithologiques et structurales sont des facteurs importants dans la localisation et l'initiation des instabilités gravitaires. De plus, le surcreusement glaciaire très important de la vallée, mis en évidence par prospection géophysique, a probablement engendré une modification importante de l'état des contraintes dans le massif rocheux. Cependant, ce processus érosif, associé à la décompression élastique générée par la fonte des glaciers et aux facteurs litho-structuraux, ne semble pas avoir été suffisant pour déclencher des instabilités gravitaires de manière immédiate après le retrait glaciaire. Les données chronologiques obtenues sur l'escarpement sommital de Séchilienne suggèrent en effet que ce mouvement a été déclenché avec un délai d'au moins 5400 ans par rapport au retrait glaciaire, et probablement suite à une modification brusque du climat correspondant à l'optimum climatique de l'Holocène. En particulier, le régime pluviométrique apparaît être un facteur majeur dans l'initiation et dans la dynamique du mouvement de terrain de Séchilienne.

Les risques naturels constituent une contrainte forte pour l'aménagement des zones de montagne. Cette thèse s'intéresse aux méthodes mises en œuvre par les pouvoirs publics pour gérer les risques d'inondations, d'avalanches et de mouvements de terrain dans les montagnes du Valais en Suisse, de la Vallée d'Aoste en Italie et des Hautes-Pyrénées en France. Une attention particulière est accordée au système politique en tant que facteur de différenciation.<br />Il existe, depuis deux siècles, une certaine convergence des discours et des politiques publiques ayant trait aux risques : le discours sur le changement climatique a succédé au discours sur le déboisement des montagnes, et le zonage réglementaire, a remplacé le reboisement et la protection des forêts. <br />Une analyse plus fine révèle cependant des différences importantes dans la répartition des compétences entre les différents acteurs, dans l'importance accordée à l'information préventive et surtout dans la pratique du zonage des risques.<br />L'organisation politique, qui confère aux collectivités locales un rôle plus ou moins important dans l'élaboration du zonage, constitue un facteur d'explication de ces différences. Entrent également en ligne de compte la culture politique au sens large ainsi que la philosophie du risque, qui associe une certaine tradition héritée de plusieurs siècles d'histoire, le contexte politique, social, économique et culturel actuel ainsi que les représentations dominantes dans les administrations.

La vitesse des ondes de cisaillement (Vs) est un parametre elastique permettant de caracteriser la compacite des formations superficielles. Il joue un role essentiel aussi bien en sismologie appliquee, dans l'estimation des effets de site ou du potentiel de liquefaction des sols, qu'en geotechnique pour la reconnaissance des glissements de terrain. Ce travail de these s'interesse aux differentes techniques sismiques permettant d'obtenir les valeurs de Vs a l'echelle de la dizaine a la centaine de metres de profondeur. Il vise a evaluer les capacites de ces methodes et a definir une strategie de reconnaissance des sites en fonction de leur complexite. La premiere partie du travail est consacree a l'evaluation de la fiabilite de l'inversion des ondes de surface pour des milieux stratifies horizontalement. En collaboration avec d'autres laboratoires europeens, les methodes de sismique active et passive ont ete appliquees sur 20 sites accelerometriques europeens pour lesquels des profils verticaux de Vs etaient disponibles a partir d'essais en forage. La comparaison des resultats montre que les profils de Vs obtenus par les techniques sismiques de surface sont comparables a ceux deduits des methodes en forage pour les sites mous (Vs30 < 600 m/s) mais fournissent des vitesses plus faibles pour les sites compacts. De plus, l'information fournie par les courbes de dispersion obtenues est apparue suffisante pour la caracterisation des sites selon la classification EC8. Sur la base des profils Vs mesures, une etude theorique sur l'inversion de la courbe de dispersion du mode fondamental de Rayleigh a ete realisee pour des sites presentant un contraste marque entre socle et sediments. Une strategie d'inversion en deux etapes, permettant de contraindre la profondeur du substratum durant l'inversion, a ete definie sur des cas synthetiques et a ete appliquee avec succes sur des cas reels stratifies horizontalement. Le cas de sites de geometrie plus complexe est aborde dans la deuxieme partie. L'inversion des courbes de dispersion obtenues par sismique active (donnees synthetiques et reelles) et passive (donnees synthetiques) au-dessus d'une interface en pente a fourni des profils de vitesse presentant un gradient de vitesse et une mauvaise estimation des vitesse et/ou profondeur du demi-espace. Une procedure a ete definie afin de detecter les sites 2D/3D, a partir des donnees des 20 sites accelerometriques europeens et des donnees du glissement de terrain d'Avignonet (Isere, France). Cette procedure, basee sur les observables de la sismique active (sismogrammes, spectrogramme spatial, comparaison de courbes de dispersion) et de la sismique passive (forme du rapport H/V, variations laterales de la frequence du pic H/V), vise a adapter le choix des techniques d'imagerie lors de la caracterisation de site. Dans la troisieme partie, la technique d'inter-correlation du bruit ambiant est appliquee au glissement de terrain d'Avignonet. Treize stations sismologiques ont ete deployees durant deux semaines sur le site et une inversion tomographique des temps d'arrivees mesures par inter-correlation a fourni des cartes frequentielles de vitesse de groupe des ondes de Rayleigh. L'inversion de ces donnees a permis d'obtenir une image 3D en Vs du glissement de terrain jusqu'a 100 m de profondeur. Enfin, l'evolution avec le temps des correlations calculees entre 2 stations permanentes sur le glissement suggere une legere diminution de la vitesse sismique au cours des trois ans d'enregistrement. Ces resultats montrent l'interet des techniques d'inter-correlation du bruit ambiant pour la reconnaissance tridimensionnelle de milieux heterogenes et le suivi temporel de l'activite des glissements de terrain.

La compréhension et la prévision de la dynamique d'évolution des mouvements de terrain peuvent être améliorées par un apport de données quantitatives sur leur géométrie tridimensionnelle, leur état de déstructuration et les circulations préférentielles de fluide. L'objet de ce travail de thèse, à la fois méthodologique et expérimental, est de mener une étude critique des apports et limites des différentes techniques géophysiques (tomographie électrique, méthodes sismiques, traitement du bruit de fond sismique et mesure de la polarisation spontanée) pour la caractérisation des mouvements de terrain. <br /><br />Nous avons mené des campagnes géophysiques sur cinq sites tests. Ces sites ont été choisis en milieu rocheux et de grande ampleur (Séchilienne, Isère, et La Clapière, Alpes Maritimes) ou de taille modérée (Ravin de l'Aiguille, Isère) et en contexte argileux (Super Sauze, Alpes de Hautes Provence, et Léaz, Ain). Nous montrons que, selon le degré de remaniement, l'étude des paramètres géophysiques tels que $\rho$ (tomographie électrique), $V_P$ (sismique réfraction et tomographie) et $V_S$ (sismique réfraction et analyse des ondes de surface) peut permettre de caractériser la géométrie du mouvement. Nous préconisons pour l'analyse des ondes de surface de contraindre la paramétrisation par $V_P$ si elle est connue en surface et par la fréquence du pic H/V. En outre, pour une géométrie régulière et des faibles gradients latéraux de $V_S$, nous montrons que la méthode H/V peut aider à obtenir rapidement des informations quantitatives sur la variation spatiale de la géométrie du substratum. <br /><br />%Nous préconisons pour l'étude de la géométrie des mouvements rocheux, le déploiement de méthodes tomographiques (électrique et sismique) et l'analyse des ondes de surface enregistrées lors de prospection sismique et par des réseaux de bruit de fond, en tenant compte des informations sur $V_P$ en surface et de la fréquence de résonance du pic H/V. La caractérisation des écoulements au sein des mouvements est améliorée lorsque l'on associe l'interprétation des mesures de résistivité et de polarisation spontanée. <br /><br />%Selon le degré de remaniement des argiles en mouvement, l'étude des paramètres géophy\-siques tels que $\rho$ (tomographie électrique), $V_P$ (sismique réfraction et tomographie) et $V_S$ (sismique réfraction et analyse des ondes de surface) peut permettre de caractériser la géométrie du mouvement. Pour une géométrie régulière et des faibles gradients latéraux de $V_S$, l'utilisation de la méthode H/V peut aider à obtenir rapidement des informations quantitatives sur la variation spatiale de la géométrie du substratum. <br /><br />Enfin, nous avons déployé un réseau d'électrodes impolarisables pour le suivi temporel de la polarisation spontanée sur le mouvement de grande ampleur de Séchilienne. Ce réseau a permis de mettre en évidence les trajets préférentiels de circulation de fluide au sein du massif et de rendre compte de leur complexité.

Les escarpements des massifs calcaires alpins sont particulièrement sensibles aux instabilités à cause de leur structure stratifiée alternant des couches de compétences différentes et des nombreux facteurs favorisant leur érosion, tels que la présence d'eau ou de glace dans les fissures de la roche. L'instabilité des escarpements calcaires se manifeste par des éboulements ou des mouvements progressifs de masses rocheuses. Ces derniers mouvements ont été quantifiés sur deux escarpements à lithologie comparable situés en Haute-Savoie, grâce à deux types de méthodes. Le G. P. S. (Global Positionning System) il permis de mesurer des mouvements significatifs sur le versant d'Allèves entre 1996 et 1998, grâce à un réseau de 18 repères. Des prismes implantés sur l'escarpement Valanginien des Tours Saint-Jacques ont permis d'effectuer des visées distancemétriques à partir de certains points du réseau G.P. S., et de quantifier les mouvements de l'escarpement grâce à la combinaison du G. P. S. et de la distancemétrie. Un traitement fin des données G. P. S. a été réalisé avec le logiciel BERNESE, afin de corriger les biais liés essentiellement aux effets de masque et troposphériques. Les coordonnées horizontales des repères ont ainsi été obtenues à 5 millimètres près. Sur l'escarpement de Passy, des inclinomètres en silice fournissent des mesures de basculement de l'abrupt au pas horaire, L'analyse de six années d'enregistrement (1991 -1996) a nécessité l'emploi de techniques de traitement du signal, afin d'isoler l'influence de la température, et les périodicités présentes dans les données. Le versant d'Allèves présente un raccourcissement horizontal révélé par G. P. S. de 1 cm/an orienté au SSE, par rapport au village d'Allèves, situé au pied du versant. Ce déplacement suggère un glissement sur une surface structurale présentant une inflexion le long d'un axe plongeant vers l'W. Une accélération significative des mouvements a été observée sur le versant durant la période 1997-98 par rapport à 1996-97. Elle pourrait traduire l'influence du séisme d'Epagny (1996) avec retard. Alors que les escarpements sommitaux montrent des déplacements horizontaux suggérant la déstabilisation du rebord des abrupts, les Tours Saint-Jacques révèlent des mouvements traduisant leur dislocation actuelle. La tour amont semble divisée en deux compartiments au moins séparés par des fractures verticales d'azimut N140°. Un décrochement dextre peut rendre compte des déplacements relatifs des prismes. Les mouvements des prismes implantés sur la tour du milieu, de part et d'autre d'une fracture verticale d'azimut N120° parcourant toute la face NW, traduisent également le comportement actuel de deux compartiments évoluant indépendemment. Enfin, les déplacements des prismes de la tour aval peuvent traduire l'ouverture actuelle de fractures verticales d'azimut N135°, parallèles aux faces NE et SW de la tour. Sur le site de Passy, aucune corrélation n'a été observée entre les précipitations et les signaux inclinométrigues, mais de fortes corrélations lient la température et les données inclinométriques dans le domaine des basses fréquences. Malgré ces corrélations, la dérive du signal montre que l'escarpement subit un basculement vers l'amont. Un axe de rotation horizontal orienté N60° peut rendre compte de ces basculements. Ceux-ci suggèrent un glissement de type rotationne!.

L’eau, par l’intermédiaire de la pression de fluides, est un phénomène déclencheur majeur de la déstabilisation des mouvements de terrain profonds. En conséquence, la caractérisation des mécanismes de déformation nécessite une bonne compréhension des processus hydrogéologiques contrôlant la déstabilisation. Les milieux fissurés et de surcroit les milieux instables présentent de fortes hétérogénéités, ce qui rend les études hydrogéologiques classiques peu adaptées. De plus, les mouvements de terrain profonds présentent des relations hydromécaniques complexes avec des évolutions significatives dépendantes du temps (déformation de type fluage). Cette thèse s’attache à caractériser les relations précipitations-déplacement du mouvement de terrain profond de Séchilienne. Un suivi saisonnier de traceurs naturels et artificiels a permis de définir un schéma conceptuel d’écoulement de l’eau souterraine sur l’ensemble du massif malgré un nombre limité de points d’intérêt hydrogéologiques. Les données de pression de fluides étant rarement mesurées, les paramètres indirects, tels que la recharge, sont souvent les seules données hydrogéologiques qui permettent de caractériser la relation précipitations-déstabilisation. Une méthode d’estimation de la recharge basée sur un calcul de bilan du sol a été développée afin d’estimer la recharge avec précision. En se basant sur le schéma conceptuel d’écoulement et le calcul de la recharge, une analyse en ondelettes couplée à un modèle numérique a permis de caractériser la relation précipitations-vitesse de déplacement. Cette modélisation tient compte de paramètres dépendant du temps et permet de simuler une déformation de type fluage (tendance pluriannuelle des vitesses de déplacement), conséquence des couplages hydro- mécaniques indirects. La caractérisation des processus hydrogéologiques contrôlant la déstabilisation a permis de définir un seuil statistique d’activation de la déstabilisation, basé sur une approche multi-dimensionnelle innovante. / Pore water pressure build-up by recharge of underground hydrosystems is one of the main triggering factors ofdeep-seated landslides. Consequently, the characterization of landslide deformation mechanisms requires athorough knowledge of the hydrogeological processes triggering the destabilization. Anisotropic andheterogeneous media combined with landslide deformation render classical hydrogeological investigationsunsuitable. Hydro-mechanical processes which lead to slope failure of deep-seated landslides are complex andare influenced by the evolution of the landslide deformation through time. This thesis aims at improving theunderstanding of the relationships between precipitation and displacement velocity based on the study of theSéchilienne deep-seated landslide. A time-related monitoring of natural and artificial tracers allows to define aconceptual groundwater flow model despite a limited number of hydrogeological points of interest. Fluid porepressures are rarely measured on landslide sites and, instead the groundwater recharge is generally used as themost relevant parameter. A parsimonious, yet robust, guideline workflow to calculate time series of groundwaterrecharge is developed. Based on the conceptual groundwater flow model and the recharge calculation, a waveletanalysis coupled to a numerical model integrating time-dependent parameters allows to characterize therelationship between precipitation and displacement velocity and to simulate the creep deformation resulting ofindirect hydro-mechanical coupling (multi-year trend of displacement velocities). The characterization of thehydrogeological processes controlling the destabilization allowed to define a statistical rainfall threshold basedon an innovative multi-dimensional approach.

Les changements du niveau de la mer sont un des indicateurs des changements climatiques qui affectent le monde depuis les dernière décennies. Dans les zones côtières, ce phénomène est à l'origine de l'érosion côtière, l'intrusion d'eau salée dans la nappe phréatique et les inondations. Les côtes du Canada ne sont pas épargnées. En effet, le niveau de la mer dans ces régions augmente actuellement et suit les tendances mondiales. Parmi ces zones côtières, celles du golfe du Saint-Laurent nous ont intéressées. Elles sont caractérisées par des propriétés physiques spécifiques comme l'isolement relatif par rapport l'océan Atlantique, la couverture de glace en hiver et la coulée de l'eau douce au printemps qui leurs confèrent un comportement particulier en termes de changements du niveau de la mer. Les changements du niveau de la mer peuvent être déterminés en intégrant les données des marégraphes et de l'altimétrie satellitaire que le résultat de cette procédure conduit à une détermination relative du niveau moyen de la mer. Cependant, parce que la structure de mesure de ces données est différente, certaines corrections doivent être appliquées. L'une des corrections parmi les plus essentielles est la longueur du Mouvement de Terrain Vertical (MTV) qui doit être prise en compte dans les données des marégraphes, car ces derniers sont connectés directement à la croûte terrestre. Différents phénomènes peuvent causer le MTV. À une échelle locale, la tectonique, l'affaissement et la sédimentation sont les facteurs qui créent le mouvement vertical et à une méso ou macro-échelle, la fonte des inlandsis qui entraîne une modification de la charge de masse de la surface de la terre, provoquent le mouvement des terres. Il existe divers modèles et méthodes de correction du MTV globaux et locaux, tels que l'Ajustement Isostatique Glaciaire (AIG), le Système Mondial de Navigation par Satellite (SMNS), l'Interférométrie par Radar Synthèse d'Ouverture (InRSO), Gravity Recovery and Climate Experiment (GRACE) et le calcul de la différence entre données marégraphiques et altimétriques par satellite. Parmi ces méthodes, l'utilisation de la vitesse verticale des données GPS est la plus précise. Cependant, en raison du faible nombre de stations GPS dans ou proximité des stations marégraphiques, l'utilisation d'une approche plus globale est inévitable. En utilisant cette correction, le niveau moyen absolu de la mer sera obtenu. Dans cette étude, l'évaluation du MTV se fait par un algorithme d'InRSO, appelé Diffuseurs Persistants InRSO (DPInRSO en anglais PSInSAR). Les images du RSO du satellite Sentinel-1 sont soumises à l'analyse DPInRSO. Afin d'évaluer la précision de cette approche, les résultats de DPInRSO (obtenus du logiciel Stanford Method for Persistent Scatterers (StaMPS)) pour deux régions de l'est du Canada (Halifax et St. John's) sont comparés à la GPS vitesse verticale en analysant les séries temporelles de données. Également, en utilisant la même méthode, le mouvement vertical des terres dans 24 stations marégraphiques du golfe du Saint-Laurent sera corrigé. L'étude du niveau moyen de la mer est mise en oeuvre en utilisant l'analyse spectrale (analyse harmonique et analyse spectrale des moindres carrés (ASMC)) afin de supprimer les effets à long terme des constituants de la marée (constituants annuels et semi-annuels). Aussi, pour évaluer le niveau moyen absolu de la mer (NMAM) dans la région considérée, le niveau moyen de la mer calculé est corrigé sur la base des résultats du MTV. En conclusion, la comparaison des résultats de la technique DPInRSO avec les données GPS indique le comportement similaire des tendances au cours de la période d'étude. En d'autres termes, la méthode DPInRSO mesure le mouvement vertical des terres en précision millimétrique dans la région d'étude (à l'exception d'une station) et elle peut être appliquée pour les régions qu'elles ne se connectent pas sur le positionnement GPS continu des marégraphes. Dans le cas d'un changement du niveau de la mer, l'incompatibilité entre les données marégraphiques et altimétriques en termes de temps rend la conclusion plus diffcile. Cependant, sur la base d'observations continues dans des séries temporelles de données altimétriques, nous avons obtenu une tendance plus homogène du niveau de la mer (niveau moyen de la mer) dans toutes les stations, mais pour évaluer le changement à long terme du niveau de la mer et en raison de l'emplacement exact de la station marégraphique, ces dernières données doivent être prises en compte (à la fois dans le calcul du niveau moyen de la mer et des constituants de la marée). Enfin, en raison du faible changement de MTV, les valeurs moyennes absolues du niveau de la mer sont très proches de celles relatives, ce qui indique que nous pouvons évaluer les changements du niveau de la mer dans le golfe du Saint-Laurent, en utilisant uniquement les données marégraphiques et altimétriques. En plus de faire des corrections, des cartes de la tendance du niveau de la mer et des différences de phase entre les stations marégraphiques seront représentées. La qualité des données des marégraphes et de l'altimétrie satellitaire sera discutée et rapportée sur les résultats tirés de l'étude. / Sea level changes are one of the indicators of climate change that has affected the world in recent decades. In coastal areas, this phenomenon has caused coastal erosion, saltwater intrusion into the groundwater and floods. These environmental changes can also be seen in Canada's coasts. Indeed, the sea level in these regions is currently increasing and follows the global trends. Among these coastal areas, those of the Gulf of St. Lawrence interested us. This area is characterized by specific physical properties such as relative isolation from the Atlantic Ocean, winter ice cover and spring freshwater flow that confer particular behavior in terms of sea level changes. Sea level changes can be determined by integrating the tide gauge and satellite altimetry data that the result of this procedure leads to a relative determination of mean sea level. However, because the measurement structures of these data are different, some corrections should be applied. One of the most essential corrections is the length of the Vertical Land Motion (VLM) which must be taken into account in tide gauge data, as these are connected directly to the Earth crust. Different phenomena can cause VLM. In the local scale, tectonics, subsidence and sedimentation are the factors that create vertical movement and in global scale, melting of ice sheets, which leads to change the mass loading of the earth's surface, cause the land motion. There are various global and local VLM correction models and methods such as, Glacial Isostatic Adjustment (GIA), Global Navigation Satellite Systems (GNSS), Interferometry Synthetic Aperture Radar (InSAR), Gravity Recovery and Climate Experiment (GRACE) and calculation of difference between satellite altimetry and tide gauge data. Between these methods, using the vertical velocity of GPS data is the most accurate one. However, due to the low number of GPS stations in or near the tide gauge stations, using a more comprehensive approach is inevitable. By applying this correction, the absolute mean sea level will be obtained. In this study, VLM is examined by one algorithm of InSAR which is called Persistent Scatterer InSAR (PSInSAR). The Sentinel-1 single-look SAR images are used for PSInSAR analysis. In order to accuracy assessment of this approach, the results of PSInSAR (obtained from Stanford Method for Persistent Scatterers software (StaMPS)) for two regions in east of Canada (Halifax and St. John's) are compared with GPS vertical velocity by analyzing their time series data. Equally, by using the same method, vertical land movement in 24 tide gauge station in Gulf of St. Lawrence will be corrected. Investigation of mean sea level is implemented by using the spectral analysis (Harmonic analysis and Least Square Spectral Analysis) in order to remove the long-term effects of tidal constituent (annual and semiannual constituents). Also, to assess the Absolute Mean Sea Level (AMSL) in considered region, calculated mean sea level is corrected based on the results of VLM. In conclusion, comparing of the results of PSInSAR technique with GPS data indicate the similar behavior of trends during the period of study. In other words, PSInSAR method measures vertical land motion in millimeter precision in study region (except of one station) and it can be applied for regions that they do not connect on continuous GPS positioning of tide gauges. In the case of sea level change, incompatibility between tide gauge and altimetry data in terms of time makes the conclusion harder. However, based on continues observations in time series of altimetry data, we obtained more homogeneous sea level trend (mean sea level) in all stations, but to evaluate the long-term sea level change and because of exact location of tide gauge station, this latter data must be considered (both in calculation of mean sea level and tidal constituents). Finally, because of small values of VLM, absolute mean sea level values are very close to relative one which indicate that, we can evaluate the sea level changes in Gulf of Saint Lawrence by using only the tide gauge and altimetry data. In addition to making corrections, the maps of sea level trend and phase differences between tide gauge stations will be represented. The quality of the tide gauge and satellite altimetry records will be discussed and reported over the drawn results of the study.