Ce travail de recherche porte sur les conditions d'alimentation et de fonctionnement des sources karstiques saumâtres de Port Miou et du Bestouan, principaux exutoires de l'aquifère karstique littoral de Port Miou. La méthodologie adoptée consiste en une approche multicritères sur l'hydrologie, l'hydrogéologie, la spéléomorphométrie, la structure et l'héritage paléogéographique du système karstique dans son ensemble. Le suivi hydrométrique des sources sous-marines se déroule sur un an et demi au cours du contexte climatique actuel de sécheresse touchant l'ensemble de la Basse Provence. L'analyse des mécanismes de mélange entre l'eau douce et l'eau de mer dans l'aquifère karstique montre que la contamination des sources s'effectue de manière pérenne par des intrusions d'eau de mer à grande profondeur dans le réseau. L'étude des crues révèle que la mise en charge du réseau karstique ne diminue pas la quantité d'eau de mer participant au mélange, chaque crue étant caractérisée par un apport supplémentaire d'eau douce qui contribue à faire chuter la salinité et la température de l'eau de la source. La révision des bilans hydriques des principaux massifs calcaires jurassico-crétacés en région marseillaise dévoile la possibilité d'importants écoulements karstiques profonds au profit de l'aquifère de Port Miou. Il convient alors de proposer que ce système se structure selon une karstification à grande profondeur et sur une superficie importante, aujourd'hui en partie noyé sous la mer. La genèse et l'organisation de ce karst, vaste et profond, s'expliquent par le modèle messinien de karstification adopté récemment pour la mise en place de nombreux aquifères karstiques périméditerranéens.

Cette thèse a pour but de contribuer à l'étude des réservoirs carbonatés karstiques par une approche intégrant modélisation géologique numérique 3D et hydrogéologie. Une méthodologie originale a été ainsi développée pour être appliquée à la caractérisation des ressources en eau souterraine de l'Unité du Beausset (Bouches-du-Rhône et Var, SE France). Cette méthodologie lie l'architecture géologique 3D du réservoir à ses propriétés hydrodynamiques extraites de son étude hydrogéologique. Cette zone est un synclinal de roches sédimentaires mésozoïques d'une épaisseur pouvant atteindre 3000m. Ces roches affleurent sur environ 650 km², et ont subi de nombreuses phases de fracturation et karstification. La synthèse des données existantes permet de définir sept intervalles lithostratigraphiques qui représentent la compartimentation verticale de l'unité. Ce découpage a servi de base à l'élaboration d'une série de coupes géologiques proposant de nouvelles interprétations structurales. Ces coupes, basées sur une récolte de données de terrain, ont permis de construire un géomodèle structural 3D du secteur. Le terrain a également fait l'objet de l'acquisition de nouvelles données hydrogéologiques : sept sources ont fait l'objet d'un suivi en continu à l'aide de sondes de conductivité électrique, température et pression (CTD) ; plusieurs campagnes hydrochimiques exhaustives ont été menées pour obtenir le contenu en ions majeurs et en isotopes de l'eau des exutoires ; deux campagnes de jaugeages ont également eu lieu. Les enregistrements CTD ont été analysés selon les techniques d'analyses des séries temporelles. Deux nouvelles méthodes de traitement du signal sont proposées pour qualifier la fonctionnalité karstique des aquifères à partir de ces enregistrements CTD. Les analyses hydrochimiques ont permis de définir des bassins hydrogéologiques d'alimentation de sources par traçage naturel. Et les jaugeages ont servis à discuter de l'établissement d'un bilan hydrique général. Ces données hydrogéologiques représentent les caractéristiques hydrodynamiques des aquifères " échantillonnés " à travers les sources. La connaissance de la géométrie 3D des couches, et de leurs propriétés géologiques et karstologiques, ouvre la voie à un changement d'échelle pour caractériser le réservoir dans son ensemble. Cette approche lie les observations de surface et de subsurface pour définir la compartimentation hydraulique de l'unité et cibler les intervalles les plus propices à l'exploitation de nouvelles ressources. La proposition de nouvelles ressources aquifères est discutée à l'échelle des sources ou à celle d'un forage profond. L'implantation de forages profonds se fait sur la base de carte de potentialité de présence de ressources aquifères. Ces cartes sont construites à l'aide du géomodèle 3D et d'une approche paramétrique de qualité des aquifères, intégrant données statiques et dynamiques. L'approche choisie pour cette étude permet donc d'intégrer dans un même cadre numérique des observations et des données, d'origines et de résolutions variées, pour caractériser un réservoir carbonaté karstique.