Cette étude s'inscrit dans un contexte de demande croissante en ressources en eau, conduisant à mieux contraindre et estimer la distribution des aquifères karstifiés qui représentent une ressource importante, notamment sur le pourtour méditerranéen. Le but de ce travail est d'élaborer une méthode prédictive de la distribution de la karstification et de l'épaisseur de massifs carbonatés karstifiés, notamment pour les karsts profonds. La méthodologie est développée à partir de l'étude des karsts de la plateforme carbonatée téthysienne du Languedoc montpelliérain; région qui a subi une histoire tectonique complexe et polyphasée, entrainant des phases de karstification successives. Elle est basée sur la compréhension des relations entre eustatisme, géodynamique et niveau de base, enregistrées par les paléokarsts de la région. Le remplissage sédimentaire des paléokarsts permet d'identifier les épisodes de karstification et de les dater ; c'est un marqueur des variations -que l'on cherche à quantifier- du niveau de base local responsable de la formation et du comblement des karsts. On analyse 3 paléokarsts d'extension régionale dont on peut étudier : le remplissage karstique, la distribution le long de profils amont-aval et l'extension verticale indicatrice de l'amplitude de l'abaissement du niveau de base dont il est issu. On en déduit les causes géodynamiques à l'origine de leur formation: - Le " paléokarst des bauxites " à remplissage caractéristique, résulte de l'évolution en 2 phases de " l'Isthme Durancien " : altération bauxitique de la couverture néocomienne sur le sud du Massif Central, liée à une inversion au Barrémien, suivi d'une surrection ≥ 1600m à l'Aptien-Albien, causant l'érosion des altérites en amont et le dépôt des bauxites détritiques en aval. - Le " paléokarst à remplissage marin paléocène inférieur " indique des variations du niveau de base local ≥350m dans la zone amont sud-cévenole. On propose que ces variations de niveau de base d'amplitude très supérieures aux variations eustatiques et dans une zone faiblement déformée, résultent d'assèchements et d'ennoiements d'un bassin endoréique à seuil. - Le " paléokarst à remplissage détritique d'origine cévenole ", enregistre la chute du niveau de base au Serravalien-Tortonien, qui se traduit par l'incision de canyons dans l'arrière-pays de la marge du Golfe du Lion. On détermine une surrection de plus de 400 mètres à l'amont, près de 250m dans les blocs faillés intermédiaires et quelques dizaines de mètres dans la zone avale. En parallèle, on établit la géométrie et la distribution actuelles des massifs carbonatés dans la zone d'étude (y compris sous couverture), en construisant un modèle géologique 3D, des Cévennes à la Méditerranée. On intègre les données de sismique réflexion, de forages et les données gravimétriques. La restauration tectonique séquentielle d'une coupe Nord-Sud extraite du modèle géologique permet d'appréhender la position des paléo-niveaux de base par rapport aux massifs carbonatés, pour chaque époque de karstification. On déduit l'organisation 3D des paléo-écoulements à travers la géométrie complexe des massifs carbonatés languedociens, notamment pendant l'évènement messinien. En terme d'évaluation de la ressource en eau, on déduit que: - Les surrections de l'arrière-pays (Crétacé " moyen ", Miocène) produisent des karsts de grande amplitude, mais ne constituent pas de réservoirs significatifs. - La karstification à l'origine de réservoirs exploitables importants et profonds est créée par l'assèchement de bassins endoréiques, qui ont affecté la région au Paléocène inférieur et au Messinien . - Ces deux évènements sont à l'origine de deux principaux réservoirs karstiques profonds qui ont une signature gravimétrique : (1) En amont du chevauchement de Montpellier, entre le Pic Saint Loup, la faille des Matelles et le bassin de l'Hérault, le karst profond se développe sur une épaisseur atteignant 900m, et est attribuée à l'événement du Paléocène. (2) Au Sud du chevauchement de Montpellier les karsts côtiers s'étendent jusqu'à plus de 10 km au large, scellés par une épaisse couverture Miocène et Plio-Quaternaire, et résultent de l'événement Messinien.

L'objectif premier de ce travail est l'étude puis la modélisation numérique 3D du fonctionnement hydrogéologique et salin actuel d'une partie de l'aquifère superficiel de la Moyenne Camargue afin de l'utiliser, dans un second temps, pour la prévision qualitative des modifications induites par le changement climatique. La caractérisation de l'aquifère a été mené selon deux approches : une étude structurale, en utilisant des sondages à la tarière et différentes techniques géophysiques (ERT, ES et cartographie EM), et une étude hydrogéologique, intégrée par plusieurs campagnes de pompages, suivi piézométrique et des prélèvements de l'eau souterraine. Cette méthodologie a servi pour caractériser la géométrie de l'aquifère superficiel ainsi que son fonctionnement actuel et les processus qui l'affectent. Divers phénomènes environnementaux ont été identifiés, tels que : un lessivage du sol grâce à l'action de l'irrigation, un endoréisme observé dans les zones basses produit par l'effet de l'évaporation, et la reproduction des fluctuations de l'étang de Vaccarès dans la piézométrie, jusqu'au moins 30 m de la rive. Le modèle numérique construit à partir des observations du terrain reproduit les principaux processus observés. Il met en évidence un possible flux sortant vers le Vaccarès et le rôle important du marais dans l'hydrodynamique et le comportement hydrosalin du système. Les résultats du modèle face aux différents scénarios de prédiction montrent une intrusion d'eau du Vaccarès et, donc, une salinisation du système induit par la remontée du niveau marin, ainsi que la présence de l'irrigation et le marais comme facteurs essentiels dans la neutralisation de cette salinisation.

L’objectif premier de ce travail est l’étude puis la modélisation numérique 3D du fonctionnement hydrogéologique et salin actuel d’une partie de l’aquifère superficiel de la Moyenne Camargue afin de l’utiliser, dans un second temps, pour la prévision qualitative des modifications induites par le changement climatique. La caractérisation de l’aquifère a été mené selon deux approches : une étude structurale, en utilisant des sondages à la tarière et différentes techniques géophysiques (ERT, ES et cartographie EM), et une étude hydrogéologique, intégrée par plusieurs campagnes de pompages, suivi piézométrique et des prélèvements de l’eau souterraine. Cette méthodologie a servi pour caractériser la géométrie de l’aquifère superficiel ainsi que son fonctionnement actuel et les processus qui l’affectent. Divers phénomènes environnementaux ont été identifiés, tels que : un lessivage du sol grâce à l’action de l’irrigation, un endoréisme observé dans les zones basses produit par l’effet de l’évaporation, et la reproduction des fluctuations de l’étang de Vaccarès dans la piézométrie, jusqu’au moins 30 m de la rive. Le modèle numérique construit à partir des observations du terrain reproduit les principaux processus observés. Il met en évidence un possible flux sortant vers le Vaccarès et le rôle important du marais dans l’hydrodynamique et le comportement hydrosalin du système. Les résultats du modèle face aux différents scénarios de prédiction montrent une intrusion d’eau du Vaccarès et, donc, une salinisation du système induit par la remontée du niveau marin, ainsi que la présence de l’irrigation et le marais comme facteurs essentiels dans la neutralisation de cette salinisation. / The main objective in this investigation is to construct a representative 3D numerical groundwater flow model of a local superficial aquifer in the Medium Camargue plain to simulate the actual hydrogeological and salt behavior. The 3D model will be used, in a second time, to simulate modifications of the system behavior provide by climate changes. The characterization of the aquifer was carried out by two approaches: a structural survey, using auger drilling technique and some geophysics techniques (EM mapping ERT, ES and), and a hydrogeological survey composed by several pumping tests, piezometric monitoring and groundwater samples. The methodological approach used allowed characterizing the geometry, behavior of the system and the process governing performance of the superficial aquifer. Several environmental phenomena were been identified, such as: a soil’s leaching; an endorheic phenomenon at the lower elevation areas; and the reproduction of the Vaccares pond’s fluctuations at the water table, until 30 m from shore. The numerical model, which has been constructed through field observations, reproduces principal process of the aquifer. This model highlights a possible outflow from the aquifer to the Vaccares pond and the important role of marsh in the hydrodynamic and hydrosaline behavior of the system. The results of the prediction scenarios tested with the numerical model highlights a possible intrusion from Vaccares water and, therefore, a salinization of the system produced by increase of sea level rise. The important presence of the irrigation and the marsh to neutralize a future salinization has been exposed.