Le mont Covey Hill constitue une vaste zone de recharge pour l'aquifère rocheux régional et se distingue par l'abondance et la diversité des milieux naturels qu'il abrite. Ces milieux se caractérisent, entre autres, par la présence de plusieurs ruisseaux et résurgences qui servent d'habitats à deux espèces de salamandres susceptibles d'être menacées au Québec. L'alimentation en eau de ces cours d'eau et indirectement de ces sources est en partie assurée par la présence d'une tourbière naturelle située près du sommet de la colline. La capacité de rétention d'une tourbière permet de réduire les pertes d'eau par ruissellement dans les cours d'eau et de limiter les impacts d'une sécheresse estivale sur le réseau hydrique. Sur le mont Covey Hill, une modification de l'équilibre hydrique et écologique de la tourbière pourrait réduire la quantité d'eau qui y est emmagasinée et ainsi déstabiliser l'ensemble de l'écosystème de la colline. Ce projet de recherche a pour objectifs de comprendre la dynamique hydrologique de la tourbière et de définir son rôle dans l'alimentation en eau de la colline et ce, dans un contexte de conservation du milieu. Pour répondre à ces objectifs, le projet de recherche se divise en deux volets: 1) réalisation d'un bilan hydrique spécifique à la tourbière et 2) développement d'un modèle d'écoulement avec MODFLOW-2000 pour approfondir la compréhension du système hydrique tourbière/aquifère. Les deux volets ont pu être réalisés grâce à l'instrumentation permanente en place au Laboratoire naturel de Covey Hill où se trouve la tourbière à l'étude. Les observations hydrologiques montrent que l'apport d'eau provient principalement des précipitations (91-92 %) dont une grande partie retourne vers l'atmosphère sous forme d'évapotranspiration (47-62 %). Les résultats ne montrent pas une contribution significative provenant de l'aquifère à proximité de la tourbière (8-9 %) et aucun apport direct de la tourbière vers l'aquifère n'a été observé. La contribution de la tourbière à l'hydrologie de la colline provient donc principalement de l'écoulement par les ruisseaux exutoires (38-53 %) lequel est alimenté en grande partie par la couche de surface de la tourbière (acrotelme). L'écoulement souterrain actuel a d'abord été simulé sur l'ensemble de la colline afin de quantifier les échanges entre la tourbière, les cours d'eau et l'aquifère et de les comparer avec les résultats du bilan hydrique. Par la suite, deux scénarios de perturbation des dépôts tourbeux ont été simulés en retirant, d'une part, la couche représentant l'acrotelme et d'autre part l'ensemble des couches représentant la tourbière. Des simulations de réduction de la recharge basées sur les prédictions de changement climatique du Groupe d'experts intergouvernemental sur l'évolution du climat (GOEC) ont ensuite été appliquées sur les trois scénarios. Ces simulations ont permit de comprendre l'impact des différentes conditions climatiques sur la dynamique de la colline et le rôle de la tourbière dans cette dynamique. Les résultats de la modélisation montrent que 99% des écoulements au sein de la tourbière se produisent dans la couche acrotelme alimentée par les eaux de pluie ainsi que par l'aquifère superficiel environnant. La modélisation permet d'identifier plus précisément la portion de l'aquifère qui contribue à l'alimentation de la tourbière. Dans un contexte de conservation du milieu naturel, cette zone indique le périmètre minimal de protection de la tourbière où toute modification de la recharge et de la qualité de l'eau doit être évitée. Les simulations de perturbation des dépôts tourbeux montrent que l'absence de la couche acrotelme produit un impact majeur sur les échanges aquifère/tourbière et sur la contribution de la tourbière à ses exutoires. Les résultats des simulations de réduction de la recharge montrent cependant l'importance de la capacité de rétention de la couche catotelme de la tourbière sur le maintien de la nappe. Pour une réduction de la recharge de 50% (prévue d'ici 2050) le modèle avec tourbière et sans acrotelme montre un abaissement de la nappe de 38 à 40 m au sommet de la colline et de seulement 6 m en bordure de la tourbière. La réduction de la recharge sur le scénario sans tourbière montre un abaissement de la nappe plus important avec 43 m au sommet de la colline et 10 m en bordure de la tourbière. Ceci montre l'effet de la capacité de rétention de la couche catotelme de la tourbière. La préservation de l'intégralité de la tourbière apparaît donc être essentielle afin d'atténuer l'impact d'éventuelles modifications du climat sur le système hydrologique de la colline et donc sur l'habitat des espèces de salamandres de ruisseaux. Les conclusions de ce projet de recherche permettent de poser les bases scientifiques à partir desquelles pourront s'appuyer les mesures de protection et de gestion du site qui doivent être mises en oeuvre par Conservation de la nature Canada (CNC) afin de maintenir l'intégrité hydrologique et écologique du milieu. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Tourbière, Hydrologie, Hydrogéologie, Modélisation, Conservation.

Le présent mémoire s'insère dans un projet de plus grande envergure ayant eu lieu de 2003 à 2006 et qui visait à caractériser qualitativement et quantitativement les eaux souterraines de la vallée d'Annapolis en Nouvelle-Écosse. Le projet intitulé « Étude des aquifères de la vallée Annapolis-Cornwallis » (ÉAVAC) était dirigé par la Commission géologique du Canada en partenariat avec le ministère de l'environnement de la Nouvelle-Écosse (Nova Scotia Environment and Labour, NSEL). L'utilisation intensive de fertilisants pendant de nombreuses années a induit une augmentation des concentrations en nitrates dans les eaux de surface et dans les eaux souterraines. La présente étude avait comme objectif de mieux comprendre la contamination de l'eau souterraine par les nitrates sur le bassin du ruisseau Thomas. Ce petit bassin (8 km²) situé au nord-ouest de la municipalité de Berwick, est considéré représentatif de l'ensemble de la vallée de par sa géologie, sa topographie et ses activités agricoles. Plus de la moitié des terres du bassin sont cultivées et les principales cultures sont la fléole des prés, la luzerne, l'orge, les fraises et le maïs. Des travaux de terrain effectués en 2005 et 2006 ont permis de recueillir plusieurs échantillons de sol et d'eau sur lesquels des analyses granulométriques et géochimiques (pH, conductivité, ions majeurs), isotopiques (isotopes stables de l'eau et des nitrates) ont été effectuées. À partir de ces résultats et d'une carte d'occupation du territoire, il a été possible de simuler la migration des nitrates au travers de la zone non saturée à l'aide d'un logiciel de modélisation (AGRIFLUX). Finalement, une évaluation de la vulnérabilité intrinsèque et spécifique à la contamination des nitrates a été effectuée à l'aide de la méthode DRASTIC standard et modifiée. Les analyses en laboratoire ont révélé que la partie centrale de la zone à l'étude montre des concentrations en nitrates élevées (>10 mg N-NO₃/L) dans l'eau de l'aquifère rocheux. Les résultats de modélisation indiquent que la culture de la fraise et de l'orge ainsi que celle de la fléole des prés provoque un lessivage élevé en nitrates. Cependant, ces cultures ne sont pas situées dans les secteurs où les plus fortes concentrations en nitrates ont été décelées. L'utilisation de la méthode DRASTIC modifiée a permis d'obtenir une meilleure corrélation entre les concentrations mesurées et les zones de vulnérabilité par rapport à la méthode standard, mais la zone caractérisée par les plus fortes concentrations en nitrates n'est toujours pas classifiée vulnérable. Les caractéristiques géologiques de la région à l'étude et les directions d'écoulement souterrain pourraient expliquer ce résultat. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Hydrogéologie, Nitrates, Modélisation, AGRIFLUX, Vulnérabilité, DRASTIC.

L'étude des écoulements souterrains et du transport de contaminant en conditions de terrain est souvent très difficile, notamment en raison de l'hétérogénéité souvent importante des formations aquifères. Pour mieux comprendre les processus en jeu, les paramètres hydrodynamiques et hydrodispersifs sont souvent estimés en laboratoire dans des conditions entièrement contrôlées. L'objectif de ce mémoire est de vérifier s'il est possible de reproduire l'écoulement et le transport de contaminant dans une maquette hydrogéologique munie d'un simulateur de pluie en laboratoire dont les sables auraient été caractérisés en laboratoire. Cinq sables sont d'abord étudiés par des techniques classiques: des courbes granulométriques ainsi que des essais de Darcy et de traçage au NaCI dans un dispositif de Darcy de 43 cm de longueur. Ces techniques ont permis de déterminer les conductivités hydrauliques et les dispersivités des sables. Des essais d'écoulement et des traçages au NaCI sont ensuite réalisés en laboratoire dans une maquette hydrogéologique représentant un aquifère de dimensions 2 m x 1.32 m x 0.25 m. La maquette est d'abord remplie d'un seul des cinq sables pour représenter un aquifère homogène. Des expériences d'écoulement et de transport sont réalisés dans cet aquifère, avec et sans la recharge imposée au moyen d'un simulateur de pluie. La maquette hydrogéologique est par la suite vidée et remplie à nouveau d'une combinaison des cinq sables disposés en blocs de 4 cm x 4 cm x 25 cm suivant un variogramme exponentiel. Les expériences d'écoulement et de transport sont réalisées avec cet aquifère hétérogène, à nouveau avec et sans recharge. Les charges, les débits et les concentrations mesurés dans toutes ces expériences sont ensuite simulés avec un modèle Modflow et un modèle MT3D. Les paramètres calés pour reproduire les résultats de laboratoire sont comparés à ceux déterminés expérimentalement. Les résultats montrent que les conductivités hydrauliques estimées au moyen des courbes granulométriques sont peu représentatives des valeurs obtenues par l'expérience de Darcy. Les conductivités hydrauliques calées dans Modflow avec et sans recharge, pour les aquifères homogène et hétérogène sont quant à elles très similaires aux valeurs obtenues avec le dispositif de Darcy. Ceci indique que les paramètres hydrodynamiques mesurés en laboratoires sont fiables et représentatives des sables testés. Par ailleurs, les dispersivités mesurées dans le dispositif de Darcy sont jusqu'à 60% supérieures à celles calées dans MT3D pour reproduire les traçages réalisés en laboratoire. Ceci s'explique par la différence d'échelle entre le dispositif de Darcy (43 cm) et les blocs de l'aquifère hétérogène (4 cm). Ces résultats confirment qu'il est beaucoup plus difficile d'estimer en laboratoire les propriétés hydrodispersives d'un matériau aquifère que ses propriétés hydrodynamiques. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Traçage, Simulateur de pluie, Maquette hydrogéologique.

Au cours des 50 dernières années, l'Afrique de l'Ouest a subi une des plus fortes variations de la pluviométrie observées à l'échelle du Globe. Ces fluctuations climatiques ont non seulement affecté le régime des précipitations mais elles ont également engendré des conséquences sur les ressources hydrologiques et végétales. Cette thèse traite de la variabilité climatique et de ses conséquences sur les ressources en eau et la dynamique des zones humides notamment des Niayes du Sénégal au cours de la période comprise entre 1950 et 2004. La thèse s'articule autour de trois volets. Le premier est une caractérisation inter-annuelle et quotidienne de la variabilité climatique par l'analyse de tendance de quelques indices de variabilité et d'extrêmes de précipitations. Le second volet porte sur la réponse des systèmes hydrogéologiques des Niayes à la variabilité pluviométrique par une analyse inter-annuelle (1958-2002) et intra-annuelle du niveau de la nappe. Deux exemples extrêmes sont analysés pour illustrer le phénomène: l'année excédentaire de 1958 et l'année déficitaire de 1974. Enfin, l'impact de la variabilité hydro-climatologique sur les zones humides a été estimé à partir d'un suivi de l'évolution des zones humides via l'analyse diachronique de photographies aériennes des années 1954, 1974 et 2000. Les résultats de l'étude montrent que le déficit pluviométrique enregistré dans la région des Niayes depuis la fin des années 1960 est lié à une réduction du nombre de jours pluvieux (2j/10 ans en moyenne) et une hausse du nombre maximum de jours secs consécutifs (4j/10 ans en moyenne) au cours de la saison de mousson entre avril et octobre. Le nombre de jours dépassant les fortes valeurs de précipitations (par rapport au 90e centile) a également subi une baisse. Par contre, les précipitations cumulées sur trois jours atteignent souvent le 1/4 du cumul annuel ce qui confirme le rôle joué par les événements extrêmes sur l'intensité moyenne et les cumuls mensuels ou saisonniers des précipitations. En plus d'une réduction de la durée de la saison pluvieuse qui est passée de plus de cinq mois à près de trois mois entre 1950 et 2004, les précipitations ont davantage baissé au Sud qu'au Nord, induisant une réduction du gradient Sud-Nord. La Nappe des Sables Quaternaires est actuellement à un niveau extrêmement bas résultant d'une longue période pluviométrique déficitaire. Le maximum piézométrique est passé, sur l'ensemble de la région, de +55 m en moyenne en 1958 à +22 m en 1994 (au dessus du niveau 0 de la mer). Les plus fortes variations sont enregistrées au cours des années 1970 et 1980 -notamment au sud de la région des Niayes -avec -0.60 m tous les dix ans. Le bilan hydrologique de l'année 1974 suggère que les facteurs naturels sont les principales causes de destockage de la nappe avec 96% des sorties dont 89% par évapotranspiration. À l'échelle saisonnière, la recharge de la nappe dépend de la fréquence, de l'intensité et de la durée des précipitations et des séquences sèches à l'intérieur de la saison pluvieuse plus que du cumul annuel. En année de précipitations excédentaires (1958), le temps de recharge est de 5 à 7 mois contrairement en année «sèche» (1974) où les fluctuations piézométriques sont quasi nulles. La dégradation actuelle des conditions hydro-climatiques se manifeste par un asséchement progressif de zones humides. Sur l'ensemble de la région des Niayes, la superficie des zones inondées en permanence est passée de plus de 1000 ha en 1954 à moins de 170 ha en 1974 pour atteindre environ 50 ha en 1982. En 2000, aucune zone inondée en permanence n'est perceptible. Dans le secteur de Mboro (15°9 N -16°54 W) les zones inondées en permanence se sont transformées progressivement en sols à hydromorphie partielle. La «reprise pluviométrique» observée à la fin des années 1990 n'a pas inititié une réaparition en surface de la nappe qui subit un «effet mémoire» dû aux déficits pluviométriques répétés. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Sénégal, Niayes, Variabilité climatique, Précipitation, Hydrogéologie, Nappe des Sables Quaternaires, Zones humides, Changement climatique.

Identifier la géométrie des écoulements souterrains est un objectif majeur de l'hydrologie. L'objectif de ce travail de thèse est de montrer que la méthode du Potentiel Spontané (PS) peut fournir des informations spatialisées sur la géométrie des circulations hydriques et sur les propriétés hydrauliques du sous-sol. La méthode PS est basée sur la mesure du champ électrique naturellement présent dans le sous-sol. Des expériences de terrain ont été menées pour mesurer la signature électrique d'un écoulement hydrique lors d'expériences de pompage et d'infiltration. Pour une meilleure compréhension du problème physique sous-jacent, des simulations de pompage et d'infiltration ont été réalisées dans une cuve en Plexiglas. Les résultats obtenus sur le terrain et en laboratoire démontrent la capacité de la méthode PS pour suivre la progression d'un écoulement et confirment l'intérêt de celle-ci comme outil adapté à la détermination du niveau piézométrique et des propriétés hydrauliques du sous-sol.

La région de Katiola, située au centre de la Côte d’Ivoire est affectée par des variations climatiques. Cette région d’une superficie de 9452 kilomètres est à cheval sur deux bassins versants, à l’ouest le bassin du Bandama et à l’est le bassin du Comoé. Malgré le manque des ressources en eau, cette zone connaît une intense activité de mobilisation et d’utilisation des eaux : périmètres sucriers dans les localités de Tafiré, les parcelles de riz irriguées et activités pastorales à travers le département. Ce travail consiste à localiser les sites potentiels des puits et forages à l’aide de la géostatistique et de la télédétection. Des outils tel que les sondes géophysiques, permettent une compréhension des caractéristiques hydrogéologiques des formations. Une compréhension de la productivité des formations géologiques par l'étude des relevés de forages et de la télédétection peut permettre une meilleure compréhension des aquifères de ces formations. L’objectif de cette étude est de mettre en place des valeurs de référence devant servir comme outils de gestion, de surveillance et voire même d’aide à la prise de décision. Egalement rentabiliser le coût des forages en proposant le maximum d’ouvrages avec beaucoup de chance de succès et établir un lien entre les formations et la circulation de l'eau souterraine (espacement, longueur, connectivité). Cette étude nous a présenté les potentialités que regorgent les images SRTM. Leur traitement nous a permis de faire une cartographie assez détaillée des limites de bassin versant, du réseau hydrographique, et d’établir l’ordre des cours d’eau selon Strahler. Le calcul de l’Index Topographique d’Humidité montre bien la corrélation entre la topographie, les sols et la qualité des sols. Malgré l’étendue de notre zone d’étude, il a été possible de trouver une corrélation entre les indices topographiques, la répartition des champs et les caractéristiques du sol. Ces résultats nous ont permis d’avoir une nette amélioration de la classification des sols et l’occupation à partir d’image LANDSAT. Les images LANDSAT, nous ont également permis de suivre le phénomène d’eutrophisation et de disparition des plans d’eau de surface dont souffre le département. Nous avons pu établir l’ordre d’apparition des végétaux sur les plans d’eau, ceci jusqu'à la disparition totale du cours d’eau / The Katiola region, located in central Côte d'Ivoire, suffers climatic variations. This region has an area of 9,452 km². It is straddling two watersheds, the Bandama basin to the West, and the Comoé basin to the East. Despite the lack of water resources, this area is experiencing an intense mobilization and use of water: sugar stores in the localities of Tafiré, paddy fields and irrigated pastoral activities across the department. This work aims to locate potential sites for drilling wells using eostatistics and remote sensing. An understanding of the productivity of geologic formations by studying the records of boreholes and remote sensing can provide a better understanding of the aquifers in these formations. The objective of this study is also to establish benchmarks to serve as management tools, monitoring and even aid in decision making. The goal is also to make a maximum of successful drillings, and establish a link between these informations (spacing, length, connectivity) and groundwater movements. This study presents the potential of the DEM SRTM images. Their processing llowed us to make a fairly detailed mapping of the boundaries of catchments, drainage network, and establish the order of watercourses according to the Strahler classification. The calculation of the Topographic Index Humidity shows the correlation between topography, soil and soil quality. Despite the size of our study area, it was possible to find a correlation between the topographic indices, the distribution of fields and soil characteristics. These results enabled us improving soil classification and occupation from LANDSAT imagery. Using LANDSAT images, we also monitored the eutrophisation and isappearance of surficial water bodies in the department. We could follow the order of appearance of plants on the water, until the total disappearance of the watercourse

Comme bien d'autres grandes métropoles construites sur des formations géologiques aquifères, Paris présente une situation hydrogéologique complexe et perturbée.L'imperméabilisation du sol modifie notablement la recharge naturelle des nappes qui se fait alors, principalement, par les fuites des réseaux (voirie, eau non potable, eaux usée …). L'urbanisation intense du sous-sol, qui représente plus de 30 % de la surface parisienne (hors bois de Vincennes et Boulogne), influence les écoulements des aquifères. A cela, s'ajoute le rôle des 20 à 30 millions de mètres cubes d'eau qui son prélevés annuellement dans les nappes phréatiques.Les aquifères parisiens sont surveillés par l'Inspection Générale des Carrières de Paris (Ville de Paris) grâce à un réseau de 305 piézomètres, qui suit en particulier les ouvrages à proximité de la Seine pour anticiper les conséquences d'une crue. La réalisation d'un modèle hydrogéologique du sous-sol parisien permettant d'évaluer les effets de tels évènements et tenant compte des spécificités de l'anthropisation présentait donc un intérêt.Dans cet objectif, un SIG en trois dimensions du sous-sol parisien a, en premier lieu, été construit. Il intègre plusieurs éléments :- La géologie parisienne depuis la craie campanienne jusqu'aux formations tertiaires. 20 000 sondages ont été exploités pour cela.- Les infrastructures souterraines : 20 000 parkings résidentiels et constructions en sous-sol, plus de 140 parcs de stationnement concédés par la ville de Paris, 197 km de galeries du métro souterrain et 46 km du Réseau Express Régional.- Les prélèvements d'eau, en se basant sur les prélèvements annuels déclarés (12 millions de mètres cubes).- Les données piézométriques issues de 669 piézomètres répartis sur Paris et des communes avoisinantes.Avant ce travail, il n'existait aucune représentation synthétique de l'organisation générale du sous-sol de la capitale et de ses aménagements.Un modèle mathématique hydrodynamique s'étendant sur l'emprise de la capitale et de la proche banlieue a ensuite été construit. Ce modèle multicouche comporte 11 couches alternant 6 niveaux aquifères et 5 niveaux semi perméables, discrétisées selon 607 171 mailles carrées de 50 mètres de coté. Au stade actuel, ce modèle a permis la réalisation de simulations exploratoires reproduisant en régime permanent la structure de l'écoulement dans les conditions d'aujourd'hui. / Paris, as many large cities built on aquifers, shows a disturbed hydrogeological environment.Land waterproofing reduces the natural recharge of the water table, then mainly driven by sewer network leakage. Underground intensive urbanisation covering more than 30 % of the Parisian surfaces (without Vincennes and Boulogne woods) influences groundwater flow. In addition, had to be taken in consideration the 20-30 millions cubic meters of water yearly pumped from the aquifer system.Parisian aquifers are surveyed by IGC (Inspection General des Carrières) from the city of Paris, thanks to a 305 piezometers network, in particular close to the Seine river, in order to monitor and anticipate a flood event. In order to appraise such events a hydrogeological model of the Parisian underground had to be setup, including anthropic specification.To achieve this model a 3D GIS of the Parisian underground has been done, integrating several topics:- Parisian geology, including Campanian chalk and tertiary formation. 20 000 geological drillings have been used.- Underground infrastructures: 20 000 residential car parks and cellars more than 140 public car parks conceded by the city of Paris, 197 km of underground subway and 46 km of RER.- Water pumping. Only declared samples (12 millions cubic meters)- Piezometric data from 669 piezometers over Paris and neighbouring townsBefore this thesis, no synthetic representation of Paris geology and underground has been done.In a second time a mathematical model including Paris and suburbs has been built. This multilayer model consists of 11 layers interbedding, 6 aquifers level and 5 semi-permeable discretized through 607 171 square mesh (50 m side). Until now, this mathematical model allowed exploratory simulation in steady state showing flow structure of the present day hydrogeological situation.

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.

Une nouvelle approche de modélisation des milieux géologiques fracturés représentés par un modèle conceptuel de fractures discrètes et déterministes est présentée dans cette thèse. L'objectif principal de l'étude est de reproduire l'hétérogénéité et la complexité des milieux poreux fracturés dans un modèle géométrique tridimensionnel afin d'effectuer des simulations numériques dans le but d'améliorer les capacités de modélisation en hydrogéologie. Ceci est réalisé à travers le couplage entre une plateforme de modélisation géologique (GOCAD) et un code numérique (HydroGeoSphere). Les principaux défis à relever sont: la représentation géométrique du réseau de fractures, la sélection d'un maillage approprié pour la discrétisation spatiale du domaine de simulation et l'adaptation du code numérique à ce maillage. La nouvelle approche est basée sur une première phase de modélisation géologique 3D, suivie par la génération d'un maillage tétraédrique 3D et par la simulation numérique de l'écoulement souterrain en conditions saturées et du transport de solutés. En général, le maillage tétraédrique s'avère plus adéquat que les maillages de blocs ou de prismes pour discrétiser les geometries complexes telles que les milieux fracturés. De plus, une définition alternative du maillage "dual", qui est essentiel pour appliquer la méthode numérique élément finis - volume de contrôle utilisée par HydroGeoSphere, est analysée et intégrée dans le code numérique. Le code numérique proposé est d'abord vérifié par l'intermédiaire de simples scénarios de simulation dont les solutions, analytiques et numériques, sont déjà connues. La complexité des simulations est augmentée de façon graduelle. L'approche de modélisation est finalement appliquée au site Olkiluoto (Finlande) où un laboratoire de recherche souterrain est en construction afin d'évaluer la faisabilité du stockage géologique profond de déchets nucléaires à haute activité. Les techniques de modélisation géologique mises au point permettent de modéliser facilement la géométrie des fractures identifiées à travers la caractérisation géologique in situ. De plus, le modèle numérique s'avère adéquat pour la simulation de l'écoulement et du transport de solutés dans ce site complexe. Ce travail de recherche présente une contribution au développement des techniques de modélisation hydrogéologique des milieux fracturés.

Dans ce travail, une modélisation hydrogéologique du bassin-versant de la Somme dans le nord de la France a été menée avec une attention particulière sur l'interaction rivière-aquifère. Le modèle couplé a été utilisé. Un Système d'Information Géographique a été incorporé à toutes les étapes du processus de modélisation pour la préparation des données d'entrée et la visualisation des résultats des simulations. Le processus commence par l'analyse du Modèle Numérique de Terrain (MNT). Ensuite, les maillages de surface et de l'aquifère ont été crées. Une discrétisation plus fine a été effectuée sur le réseau des rivières et les limites des sous-bassins. Pour le modèle de surface, des données météorologiques, d'occupation du sol et de type de sol ont été acquises. Le modèle de surface permet de répartir les precipitations en évapotranspiration, infiltration et ruissellement de surface. Des niveaux piézométriques en régime permanents sont calculés par le modèle de nappe. L'écoulement dans la zone non saturée a été simulé par un modèle utilisant la cascade de Nash. Les simulations d'écoulement transitoire souterrain et de surface ont été réalisées en prenant en compte les interactions nappe-rivière sur un pas de temps journalier. La calibration et la validation des résultats ont été faites en utilisant les mesures des débits et des niveaux piézométriques dans le bassin. La forte influence de la nappe sur le regime hydrologique du bassin est bien représentée par le modèle. Des comparaisons de prédiction des zones inondées en 2001 ont été effectuées avec d'autres modèles. Enfin, des analyses de sensibilité ont été réalisées afin d'orienter de futures collectes de données terrain.