L'acquisition récente de 3000 km de données de sismique sparker haute résolution et morphobathymétriques EM 1000 dans le domaine côtier sud armoricain apporte des éléments qui éclairent l'histoire cénozoïque du Sud Bretagne et notre compréhension de la genèse des vallées et estuaires quaternaires. L'étude est basée sur la cartographie des éléments morphologiques, stratigraphiques et structuraux de la prolongation en mer des principaux fleuves de l'Odet, du Blavet, d'Etel et de la Vilaine. Elle a permis la reconnaissance de cinq paléo-réseaux hydrographiques. L'étude de l'évolution des profils longitudinaux et transversaux de ces réseaux, des surfaces d'abrasion marines qui les recoupent, de l'architecture interne de leur remplissage, et leur comparaison avec les réseaux contemporains préservés le long de la marge atlantique ont permis de discuter l'impact de la lithologie, de la tectonique et du climat sur la dynamique des paléo-réseaux hydrographiques.

L'objectif que nous avons essayé d'atteindre dans le cadre de cette thèse est l 'évaluation de l'aléa sismique au Maroc par l'application de différentes étapes. La première partie de ce travail consiste en l'élaboration d'un fichier informatisé (chapitre 1) de la sismicité de la région comprise entre 21°N-38°N et O°W-20°W à partir du début de ce siècle. Un peu plus de quatre mille six cent événements ont été recensés faisant l objet d'une révision complète, 34% des événements ont été relocalisés. Nous avons effectué un nouveau dépouillement des archivés au D.P.G., ce qui nous a permis de traçer une cartes macrosismiques des séismes largement ressentis. Nous avons essayé d'utiliser ce fichier pour calculer l'aléa sismique à l'aide du programme de calcul EQRISK développé par McGuire (première partie: chapitre 2). Les résultats obtenus sont donnés sous forme de cartes des valeurs d'intensités susceptibles d'être ressenties pour les deux valeurs de probabilité annuelle: 0.01 et 0.005. Dans le premier chapitre de la deuxième partie de cette tIhèse, nous avons procédé à une révision complète des données macrosismiques et instrumentales du séisme d'Agadir, ce qui nous a permis de dresser une nouvelle carte macrosismique et de mieux préciser le mécanisme au foyer. Une étude de l'effet de site a été effectuée pour mettre en évidence l'influence des couches sédimentaires, sous la ville d'Agadir, sur les constructions pendant le tremblement de terre. Au dernier chapitre, de la deuxième partie, une étude sismotectonique détaillée a été réalisée dans la province d'Al Hoceima dans le cadre de l'Action Inteqrée 359/80. Cette région est l'une des plus sismiques au Maroc, de nombreuses secousses y sont ressenties chaque année . Une compagne de microsismicité a été effectuée à proximité de l'accident de Nekor. Un réseau, composé de 16 stations sismiques portables et de 9 stations numériques à 3 composantes télécommandées, a été installé pour une durée de cinq semaines . L'objectif de cette expérience était de localiser précisément la sismicité, de calculer des mécanismes à la source, d'entreprendre une étude du mouvement du sol en réponse à une sollicitation sismique ..

L'ophiolite d'Oman permet d'étudier un morceau de lithosphère océanique formé à une dorsale rapide et échoué sur le continent. La structure générale permet de reconnaître la présence de plusieurs diapirs dont certains alimentaient la dorsale d'origine, partiellement obductée avec l'ophiolite. Le diapir de Maqsad est l'exemple type d'un diapir à l'axe. Les linéations verticales couvrent 100 km2 et tournent à l'horizontale au sein d'une zone de transition riche en magma, épaisse de 500 m, juste sous le Moho. Un flux forcé et radial est issu du diapir. Un calcul prédit que ce flux divergent devient passif à 10 km du diapir, comme le suggèrent nos données de terrain. Le diapir de Nakhl a des caractères similaires, tandis que le diapir de Mansah, éloigné de l'axe de la paleodorsale, est bordé par des zones de cisaillement et des injections de diabase qui attestent qu'il poinçonnait une lithosphère refroidie. Dans une seconde partie, les traces de magma dans les péridotites sont étudiées par analyse d'image sur des lames minces. On montre que les poches de liquide ont un allongement préférentiellement parallèle à la linéation et proportionnel à la force de la fabrique cristallographique de l'olivine. Ces données sont utilisées pour calculer les propriétés sismiques des échantillons. Nos résultats suggèrent que la quantité de liquide au Moho de la dorsale Est Pacifique peut être sous estimée à cause de l'anisotropie des roches. Enfin, une modélisation sismique d'un diapir, à partir de nos données structurales et des vitesses sismiques calculées permet de comparer les données d'Oman et de l'EPR.

La tectonique des plaques repose sur le postulat que les plaques sont rigides et que l'essentiel des déformations se concentre sur leurs frontières, généralement étroites en domaine océanique. Plusieurs plaques océaniques semblent y déroger. La plaque Indo-Australienne montre par exemple une forte sismicité en dehors de ses frontières, qui traduit une déformation intraplaque intense. Les mécaniques au foyer expriment une déformation compressive dans la partie Est de la plaque, et une déformation extensive sur le flanc Est de la dorsale centrale Indienne et sur le flanc Nord de la dorsale Sud Est Indienne. L'inversion de données cinématiques, qui supposent des plaques rigides, abouti au concept de plaques composites composées d'éléments rigides séparés par des limites larges et diffuses. En utilisant l'exemple de la plaque Indo-Australienne, ce travail cherche à savoir comment se comporte la lithosphère océanique quand elle est soumise à des contraintes horizontales et le rôle joué par les structures préexistantes comme les plateaux sous-marins ou les zones de fracture fossiles sur la déformation. La première partie de ce travail s'est intéressée à la sismicité de faible magnitude détectée par un réseau temporaire de 3 hydrophones déployés dans l'océan Indien Austral. L'objectif était de mieux identifier et localiser les structures actives des zones déformées. En raison de la géométrie du réseau, ces données apportent peu d'informations nouvelles mais confirment le rôle des zones de fracture fossiles dans la déformation des domaines en extension. Par ailleurs, la découverte, sur des zones de fracture actives des dorsales Est Pacifique et Juan de Fuca, de précurseurs sismiques nous a incités à en rechercher sur nos données. Un précurseur sismique est défini comme un événement se produisant dans un rayon de 100 km autour d'un événement majeur, quelques heures à quelques jours avant celui-ci. L'analyse d'une sélection d'évènements détectés par le réseau permanent du traité international de non-prolifération nucléaire (OTICE), dont la localisation est très précise, n'a pas permis de reconnaître de tels précurseurs. En revanche, notre réseau temporaire détecte de nombreux évènements qui pourraient y être assimilés, si leur localisation, basée sur seulement 3 hydrophones, pouvait être améliorée. La combinaison des données de ces 2 réseaux serait nécessaire pour progresser sur ces deux questions. La seconde partie de ce mémoire est consacrée à la modélisation de la déformation de la plaque Indo-Australienne, afin notamment de préciser le rôle de structures préexistantes, reliefs ou zones de fracture fossiles remobilisées, comme en atteste la sismicité. Pour aborder cette question, nous avons testé deux outils " ADELI " (Hassani et al., 1997) et " SHELLS " (Bird, 1989; Bird et Kong, 1994; Kong et Bird, 1995; Bird, 1999). Le second code numérique s'est avéré le plus adapté. Ce modèle décompose les plaques en éléments finis triangulaires sur la sphère, tient compte de l'âge de la lithosphère considérée, et peut prendre en compte des zones de faiblesse (zones de fracture) et des reliefs. A partir de conditions aux limites en vitesse et d'hypothèses sur la rhéologie de la lithosphère, le modèle calcule l'état des contraintes et des déformations de la plaque considérée. La plaque Indo-Australienne présentant des zones déformées disjointes, nous avons d'abord modélisé chacune d'elles séparément : la zone en extension, entre les plaques Inde et Capricorne est découpée par un réseau de zones de fracture oblique par rapport à l'orientation générale Nord-Sud de la dorsale Centrale Indienne et comprend le relief des Chagos ; la zone en extension entre les plaques Capricorne et Australie est découpée par un réseau de zones de fracture perpendiculaires à l'orientation générale Ouest-Nord-Ouest Est-Sud-Est de la dorsale Sud-Est Indienne et comprend le plateau océanique de Saint-Paul et Amsterdam ; enfin, la zone en compression entre les plaques Inde, Capricorne et Australie, est caractérisée par une lithosphère ancienne découpées par des zones de fracture Nord-Sud dans les bassins de Wharton et Indien Central et la ride du 90ème Est. Les conditions aux limites de ces modèles sont issues d'une approche cinématique de plaque rigide. Pour intégrer les limites de plaques diffuses sur les bords des maillages, nous avons utilisé la notion de bord libre, cela nous a permis de reproduire les déformations dans ces trois régions sans imposer de limites particulières entre les plaques. Enfin, nous avons réuni les trois maillages afin de modéliser la plaque Indo-Australienne dans son ensemble. Les modélisations effectuées ont mis en évidence l'importance des structures préexistantes comme le banc des Chagos qui concentre la déformation au niveau de la limite diffuse entre les plaques Inde et Capricorne où des séismes de magnitude 8 se sont produits. Les zones de fracture jouent par leur réactivation un rôle essentiel, notamment à l'Est de la ride du 90ème Est, où elles provoquent des comportements différents entre le bassin Indien Central où la déformation est marquée essentiellement par de la compression dans une direction Nord-Sud et le bassin de Wharton, où du décrochement sénestre se manifeste en plus de la compression. Nos modèles prédisent des taux d'extension et de raccourcissement proches d'études précédentes avec un éloignement de l'Inde par rapport à la plaque Capricorne d'environ 2 mm/an dans une direction Nord-Sud, et un raccourcissement de l'ordre de 11,5 mm/an entre l'Australie et l'Inde dans une direction Nord-Ouest Sud-Est. Au premier ordre, le concept de plaque composite semble valide pour la plaque Indo-Australienne, cependant la répartition des zones diffuses s'avère plus complexe que prévue.

Nous étudions ici différents aspects de la sismicité du Sud-Est de la France, en nous aidant des données nouvelles disponibles aujourd'hui. L'apparition récente de nouveaux réseaux sismologiques (L.O.G., S.E.O., I.G.G.) dans les Alpes, et la multiplication d'études géophysiques (grands profils, aéromagnétisme,paléomagnetisme,photos satellites, nivellement) et géologiques devraient aboutir dans les prochaines années à une synthèse géodynamique des Alpes. Notre travail apporte certains éléments Sismologiques pour cette future synthèse. Nous présentons d'autre part une étude sur une nouvelle méthode de zonage sismique que nous appliquons en particulier au ca's du S.E. de la France. Les trois premiers chapitres concernent surtout l'acquisition ou le traitement de nouvelles données de sismologie. Dans le Chapitre l, on constitue un fichier de séismes historiques et récents, avec lequel on trace des cartes de sismicité et on étudie la loi magnitude-fréquence. Le Chapitre II livre des solutions focales nouvelles. Le Chapitre III donne les résultats d'une campagne d'enregistrement de petits séismes en Ubaye-Queyras en septembre-octobre 1977 . Le Chapitre IV est consacré à la synthèse des données acquises dans les trois premiers chapitres et à leur comparaison avec d'autres données (géophysiques et géologiques). Enfin, le Chapitre V expose notre méthode de zonation sismique. Nous présentons en annexe le fichier de séismes, les cartes de sismicité, ainsi qu'un exposé sur la loi magnitude-fréquence.

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.

Ce rapport résume les recherches que j'ai entreprises depuis la fin de mon doctorat de thèse à l'automne 2001. Le travail présenté a été fait dans le cadre de divers projets et en collaboration avec de nombreux chercheurs : l'Imperial College London, Royaume-Uni (2001–2004) ; Bureau de Recherches Géologiques et Minières (BRGM), France (2004–présent) ; le Centre de recherches pour le génie parasismique (Earthquake Engineering Research Centre) à l'Université d'Islande (2009– 2010). En outre, ce rapport énumère les activités d'enseignement, de supervision et de conseils dans lesquels j'ai été impliqué depuis 2001. Ma recherche s'est principalement focalisée sur la prédiction de mouvements fort du sol à des fins d'ingénierie, soit pour des projets de conception ou de rénovation, soit pour l'évaluation de l'aléa et du risque sismique. La plupart de mes études portent sur l'estimation des mouvements sismiques empiriques par le biais des équations de prédiction de mouvements du sol (GMPEs, aussi appelées modèles de mouvements du sol ou relations d'atténuation). Cette recherche met l'accent sur : l'amélioration des prévisions des mouvements du sol médians et de la variabilité associée ; la quantification, la compréhension et potentiellement la réduction de la variabilité ; l'évaluation et la modélisation de la dépendance régionale des mouvements du sol ; et la confrontation de simulations et d'estimations empiriques. Mes recherches montrent que, bien que des progrès significatifs ont été accomplis au cours des deux dernières décennies dans l'amélioration de la précision des estimations de mouvements du sol médians pour un scénario donné, l'incertitude épistémique reste élevée et cela doit être pris en compte aussi bien dans l'évaluation de l'aléa que du risque sismique. En outre, toutes les méthodes permettant de réduire les écarts-types de GMPEs proposées jusqu'à présent, même si elles semblaient prometteuses, se sont révélées largement inefficaces. Mes études ont démontré que les mouvements sismique du sol peuvent varier considérablement selon les séismes ou les sites d'étude. Cela aussi doit être pris en compte lors de l'évaluation de l'aléa et du risque sismique. Toute tentative de réduire les incertitudes épistémiques et de dériver les GMPEs avec des écartstypes plus faibles sont tributaires d'une part de l'augmentation de la densité des réseaux de mouvements forts ; et d'autre part et peut-être de manière plus importante encore, de l'amélioration de la precision des métadonnées associées aux enregistrements de mouvements forts. Une telle base de données ainsi améliorée devrait conduire à une meilleure compréhension des phénomènes physiques et, de ce fait, à des résultats empiriques concernant l'effet de source, de trajet et de site sur les mouvements sismiques.

Mes travaux de recherche s'inscrivent dans l'imagerie sismique des premiers kilomètres du sous-sol, en particulier dans le domaine pétrolier. Il s'agit de reconstruire les propriétés de la sub-surface, par exemple la vitesse de propagation des ondes de compression, à partir d'enregistrements en surface de déplacements ou de variations de la pression liés au passage d'un train d'ondes acoustiques ou élastiques. Je me suis attaché, dans ce domaine, à regarder différemment les données sismiques. Les enregistrements sismiques dépendent de la position de source et d'un ensemble de récepteurs, ainsi que du temps d'enregistrement. Ils sont couramment analysés trace par trace, par exemple pour filtrer les données ou bien pour les migrer, c'est-à-dire retrouver en profondeur les perturbations du milieu, typiquement avec des migrations de type Kirchhoff. Dans le cas de la migration dite "reverse-time migration", les points de tir sont vus comme un ensemble : l'information de tous les récepteurs est rétro-propagée depuis la surface jusqu'en profondeur, également pour retrouver les propriétés du milieu. Je propose de considérer des groupes de traces adjacentes, autour d'une fenêtre en temps, et d'examiner différents algorithmes d'imagerie sous cet aspect. L'aspect cohérence locale est justifié dans le domaine pétrolier par les acquisitions actuelles qui sont de plus en plus denses, et aussi par la notion de zone de Fresnel. La cohérence latérale des signaux sismiques est traitée sous deux angles. Le premier donne une vision très épurée : les événements localement cohérents sont décrits en 2D par la position de la trace centrale, une position dans la trace, et une pente qui indique la cohérence. Je montre qu'à partir de ce type d'approche, il est possible de retrouver les grandes longueurs d'onde du modèle de vitesse par analyse de vitesse dans le domaine migré. J'ai pu aussi établir le lien entre cette méthode d'analyse de vitesse et la tomographie de pente. Le second angle d'attaque prend en compte la signature des données sismiques, et en particulier la bande passante limitée. Les applications sont alors beaucoup plus nombreuses : migration, analyse de vitesse, sensibilité de l'image migrée par rapport au modèle de vitesse, et d'autres tâches liées au pré-traitement comme le débruitage, la prédiction des multiples, ... Dans chacun des cas, j'analyse les avantages et limitations de l'utilisation de la cohérence locale. Il s'avère que le choix du code de décomposition en événements locaux (``curvelets'', ...) est étroitement lié à l'application qui en est faite : si l'objectif est de comprimer les données sismiques, alors les curvelets ne sont pas adaptées. S'il s'agit au contraire de savoir comment l'image migrée dépend du choix du modèle de vitesse pris pour la migration, alors les curvelets ont beaucoup d'avantages. Plus particulièrement, les curvelets offrent la flexibilité de la décomposition des données, très utile pour la suppression des bruits cohérents. De plus, elles diagonalisent presque l'opérateur de démigration/migration qui donne la sensibilité de l'image migrée par rapport au choix du modèle de vitesse. Enfin, elles permettent d'exploiter l'aspect multi-échelle des données, avec des applications sur la suppression de l'aliasing. Je montre aussi qu'il est possible de développer d'autres schémas de décomposition et de reconstruction, comme par exemple les ``circlets''. Dans les perspectives de recherche, je tiens à dépasser le cadre asymptotique hautes fréquences dans lequel les curvelets trouvent un cadre naturel. Je propose en particulier une nouvelle formulation du problème d'inversion des formes d'onde (DFWI, Differential Full WaveForm Inversion). Cette méthode se veut générale pour retrouver les propriétés du sous-sol et permettre une optimisation locale, tout en s'affranchissant de la détermination d'un modèle de départ très proche de la solution. Cette idée doit être approfondie. Elle exploite la cohérence locale des données sismiques et suppose que le signal est bien échantillonné en temps et en espace. Les applications sont très nombreuses et touchent le domaine pétrolier, la géotechnique, la sismologie globale et régionale.

Cette thèse vise à mettre en évidence expérimentalement le mécanisme d'interaction entre une ossature en béton armé et son remplissage en maçonnerie de briques creuses sous sollicitations quasi statiques alternées. L'étude expérimentale montre que la maçonnerie est capable d'absorber une bonne quantité d'énergie et que la rigidité et la résistance de la structure sont augmentées par rapport au cas du cadre vide. Le calcul montre que le cadre rempli peut être convenablement modélisé en remplaçant la maçonnerie par un système de bielles non-linéaires. Le modèle à 3 bielles diagonales parallèles rend mieux compte des efforts internes dans l'ossature.

La compréhension du fonctionnement hydrodynamique des grands complexes aquifères nécessite une bonne connaissance des réservoirs sollicités. Au sein de ces ensembles, l'influence des structures géologiques de grande extension reste mal connue. La structure de Villagrains-Landiras, limite méridionale du système aquifère nord aquitain, sert ici de cadre pour étudier l'impact de ces structures sur l'alimentation de systèmes aquifères profonds.<br />Différentes méthodes, faisant appel à plusieurs domaines des géosciences (géologie, géophysique, hydrogéologie, hydrologie, géomorphologie, géochimie isotopique, ...) ont été employées pour caractériser le fonctionnement hydrodynamique de cette entité et établir son rôle au sein du complexe aquifère nord aquitain.<br /><br />L'interprétation des données issues des prospections effectuées (sondages de reconnaissance, sismique réflexion) a aboutit à la révision de la géologie du secteur et des conditions de gisement des réservoirs rencontrés. Les travaux réalisés montrent l'importance des phénomènes érosifs qui ont affecté l'ensemble étudié induisant l'absence de nombreux horizons aquifères du Tertiaire (Éocène, Oligocène et Miocène). Les formations du Crétacé supérieur, largement fissurées et karstifiées, occupent quant à elles une place majeure au sein de la structure.<br /><br />L'étude du fonctionnement hydrodynamique de la structure s'est appuyée sur une étude transversale intégrant diverses analyses.<br />Tout d'abord, l'examen des modelés et des différentes chroniques piézométriques a montré l'importance de la structure dans la mise en charge des aquifères du complexe nord aquitain.<br />Ensuite, les différentes mesures effectuées sur le terrain ont mis en évidence l'existence, dans certains secteurs, d'une relation étroite entre les eaux de surface et de sub-surface et les eaux souterraines. Les cours d'eau recoupant l'axe de la structure sont ainsi caractérisés par des échanges avec les réservoirs souterrains (pertes et/ou émergences).<br />Par ailleurs, les nombreuses dépressions circulaires immergées présentes sur l'axe de la structure de Villagrains-Landiras, sont des marques de l'existence d'une communication privilégiée entre les eaux de la nappe de surface et les eaux du réservoir de la base du Crétacé supérieur. Ces formes caractéristiques ont une origine complexe combinant plusieurs processus (thermokarstique et cryptokarstique), la karstification des formations du secondaire ayant abouti à la suffosion des terrains meubles du Plio-Quaternaire.<br />En définitive, les différentes analyses isotopiques effectuées viennent appuyer les précédentes observations. Les secteurs dans lesquels des traces de tritium ont été retrouvées coïncident avec les zones d'échanges mises en évidence par les mesures en rivières et l'étude des dépressions ennoyées. D'une manière générale, la structure de Villagrains-Landiras apparaît comme une zone de recharge, les temps de résidence étant plus courts sur l'axe de la structure et ce pour différentes nappes.<br /><br />Enfin, les écoulements au sein du complexe aquifère nord aquitain ont été appréhendés à l'aide d'un modèle hydrodynamique régional représentatif des échanges entre les 6 aquifères composant ce système. La géométrie des réservoirs a complètement été reprise sur la zone de Villagrains-Landiras de manière à restituer plus fidèlement le fonctionnement hydrodynamique de cet ensemble.<br />Ainsi, ce modèle a pu être utilisé pour mesurer l'impact à grande échelle de l'implantation de champs captants exploitant la nappe de la base du Crétacé supérieur sur la commune de Saint-Magne.