La Zone Non Saturée (ZNS) peut jouer un rôle majeur pour la circulation et le stockage de l'eau dans les hydrosystèmes karstiques mais ces mécanismes restent peu connus. L'hydrogéophysique peut-elle permettre d'évaluer les variations spatiales et temporelles du stock de l'eau de la ZNS pour améliorer la compréhension de son fonctionnement ? En 2011 et 2012 différents outils ont été utilisés sur un site expérimental à l'aplomb du Laboratoire Souterrain à Bas Bruit de Rustrel (LSBB) dans le bassin d'alimentation de la Fontaine de Vaucluse afin d'en étudier la structure et le fonctionnement hydrogéologique. Ces investigations se sont déroulées selon une méthodologie comprenant quatre phases : (1) phase préliminaire d'étude géologique et hydrogéologique (2) caractérisation géophysique de la stucture géologique (3) caractérisation hydrogéophysique du fonctionnement hydrodynamique (4) synthèse par géomodélisation 3 D. Les résultats obtenus sur le site expérimental soulignent le potentiel de l'hydrogéophysique en milieu karstique et valident la méthodologie proposée. Les investigations ont permis d'identifier rapidement des structures géologiques pouvant jouer un rôle majeur dans la régulation des flux d'eau. Des variations saisonnières du stock d'eau ont été mesurées jusqu'à plusieurs dizaines de mètres dans la ZNS. Le suivi géophysique évènementiel a permis d'imager l'hétérogénéité spatiale et temporelle de la recharge. Enfin, la géomodélisation s'avère essentielle dans la méthologie proposée pour synthétiser les connaissances du site, finaliser le schéma de fonctionnement et poser les bases de futurs travaux de modélisation hydrodynamique dans la ZNS du karst. / Karst Unstatured Zone can play a major hydrodynamic role but water flow and storage mechanisms remains poorly known. Can hydrogeophysics help assessing the spatial and temporal waterstorage variation in UZ and improve understanding karst UZ functioning ? In 2011 and 2012, various tools have been employed on an experimental site above the low noise underground research laboratory of Rustrel (LSBB) within the Fontaine de Vaucluse catchment, in order to study the hydrogeological structure and functioning. These investigations were conducted using a methodology including four phases : (1) Preliminary study of geological and hydrogeological contexte, (2) Geophysical characterization of the geological features, (3) Hydrogeophysical characterization of the hydrodynamic functioning, (4) 3D geomodeling of the karst UZ. The acquired results highligth the potential of hydrogeophysical investigation in karst media and validated the proposed methology. These investigations allowed a fast identification of geological structures which may play a major role in water flow regulation. Seasonal variations in water storage were measured until several tens of meters deep within the UZ. Hydrogeophysical monitoring allowed imaging spatial and temporal heterogeneityof recharge; Finally, the geomodeling is essential within the proposed methology to synthesize site knowledge, finalize the conceptual functioning schema and lay the foundation for future work of hydrodynamicmodeling in karst UZ / La Zona No Saturada (ZNS) puede jugar un papel importante en la circulación y el almacenamiento del agua en los hidrosistemas kársticos, sin embargo sus mecanismos son poco conocidos. ¿La hidrogeofísica podría evaluar las variaciones espaciales y temporales del agua subterránea de la ZNS para mejorar la comprensión de su comportamiento?En el año 2011 y 2012, diferentes herramientas fueron utilizadas en un sitio experimental del Laboratorio Subterráneo de Bajo nivel de Ruido (LSBB, siglas en francés) ubicada en la cuenca de alimentación de la Fontaine de Vaucluse (Francia) con el propósito de estudiar su estructura y comportamiento hidrogeológico. Estas investigaciones fueron desarrolladas según una metodología compuesta de cuatro fases: (1) una fase preliminar de estudio geológico e hidrogeológico, (2) una caracterización geofísica de la estructura geológica, (3), una caracterización hidrogeofísica del comportamiento hidrodinámico, y (4) una síntesis por geomodelización 3D.Los resultados obtenidos en el campo validan la metodología utilizada y muestran la importancia del uso de la hidrogeofísica en el estudio del medio kárstico. Las investigaciones permitieron identificar rápidamente las estructuras geológicas que juegan un rol importante en la regulación de los flujos de agua. Las variaciones temporales del agua subterránea fueron medidas hasta varias decenas de metros de profundidad en la ZNS. El monitoreo geofísico de eventos puntuales permitió realizar una representación grafica de la heterogeneidad espacial y temporal de la recarga. Finalmente, la geomodelización es considerada como parte esencial en la metodología propuesta para sintetizar los conocimientos del área de estudio, finalizar el esquema del funcionamiento y proponer las bases de los futuros trabajos en modelización hidrodinámica de la ZNS del karst. / يمكن للطبقة غير المشبعة٬ في نظم المياه الجوفية ذات الطبيعة الجيرية (الكارست) أن تلعب دورا مهما في حركة و تخزين المياه الجوفية٬ لكن آليات عمل هذه الطبقة تبقى إلى حد ما مبهمة. إن السؤال الذي يبقى مطروحا : هل بإمكان علم الهيدروجيوفيزياء أن يمكننا من تقييم التغيرات في البعدين المكاني و الزمني لمخزون المياه في الطبقة غير المشبعة بغية تطوير فهمنا لآلية عملها؟خلال الفترة الممتدة بين سنتي 2011 – 2012 م٬ تم إستخدام العديد من الأدوات في حقل تجريبي على مستوى المختبر التحت أرضي عديم التشويش بروستال بفرنسا (LSBB de Rustel)٬ والذي يتواجد في حوض التغذية لمنبع الفوكلوز(Fontaine de Vaucluse) بغية دراسة بنية الطبقة غير المشبعة للكارست و كيفية أدائها فيما يتعلق بالمياه الجوفية. إن هذه الأبحاث نفذت بمنهجية مكونة من أربعة مراحل٬ و هن كالآتي: (1) المرحلة الابتدائية و تخص الدراسة الجيولوجية و الهيدروجيولوجية٬ (2) التوصيف الجيوفيزيائي للبنية الجيولوجية٬ (3) التوصيف الهيدروجيوفيزيائي لحركة الطبقة المائية٬ (4) الخلاصة من خلال إنشاء نموذج جيولوجي ثلاثي الأبعاد.إن النتائج المتحصل عليها في الموقع التجريبي أظهرت أهمية الهيدروجيوفيزياء في الوسط الكارستي و أكدت نجاعة المنهجية المقترحة. إن الأبحاث المنجزة٬ مكنت في وقت وجيز٬ من تحديد البنيات الجيولوجية التي يمكن أن تلعب دورا مهما في معادلة التدفقات المائية. من منظور آخر٬ مكنتنا الطريقة المقترحة من قياس التغيرات الموسمية للمخزون المائي في الطبقة غير المشبعة (ZNS) و ذلك على عمق عشرات الأمتار. إن الرصد الجيوفيزيائي الحدثي مكن من تصوير التغيرات المكانية و الزمنية لإعادة التغذية. أخيرا يظهر جليا أن النمذجة الجيولوجية ثلاثية الأبعاد ذات الأسس الجيولوجية٬ الهيدروجيولوجية و الجيوفيزيائية٬ أمر أساسي بالنسبة للمنهجية المقترحة من أجل حوصلة معلوماتنا بخصوص الموقع المدروس٬ تكملة مخطط التشغيل و وضع أسس الأعمال المستقبلية لنماذج حركة المياه خلال الطبقة غير المشبعة لأحواض المياه الجوفية ذات الطبيعة الجيرية (الكارست Karst).المفردات الرئيسية:الهيدروجيوفيزياء (علم جيوفيزياء المياه الجوفية)، الكارست (حوض مائي جوفي ذو طبيعة جيرية)، إعادة تغذية الطّبقة المائيّة، الطّبقة غير المشبّعة (ZNS)، النّمذجة الجيولوجيةGéomodélisation) )، منبع الفوكلوز (ontaine deF Vaucluse)، المختبر التحت أرضي عديم التشويش بروستال بفرنسا (LSBB de Rustel). / 非饱和带对喀斯特地区水分的运动和存储具有重要作用，但是这一机制仍然不太清楚。水文地球物理学方法是否有助于评估非饱和带水分存储的时空变异，并深化我们对喀斯特非饱和带功能的认识？为了研究沃克吕兹泉流域的水文地质结构和功能，分别于2011和2012年，在吕斯特勒地区地下低噪声试验室(lsbb)上方的试验点进行了各种试验。研究内容包括四个部分：(1)：初步的地质和水文地质调查；(2)：地质结构的地球物理方法勘探；(3)：水动力过程的水文地球物理监测；(4) 喀斯特非饱和带结构的三维模拟。研究结果表明，水文地球物理学方法在喀斯特地区的研究中具有很大应用潜力，并检验了我们提出的方法的可靠性。利用这种方法，对对水流调节起主要作用的地质结构进行了快速识别；对深达几十米的非饱和带储水量的季节性变化进行了监测；对地下补给过程的时空异质性进行了成像。最后，综合研究区的各种信息，建立地学模型，确定了水文地质结构和功能的概念模式，为后续喀斯特非饱和带水动力学过程模拟研究奠定了基础。

Les résultats présentés dans ce document portent sur les anisotropies sismique et électrique d'un massif calcaire fracturé dans des galeries parallèles souterraines. Les fractures étant orientées verticalement, leurs propriétés sont approximées par isotropie transverse horizontale (HTI). Plusieurs méthodes d'inversion ont été étudié pour traiter les données sismiques (temps d’arrivée des ondes P et S); et pour traiter les données de résistivité électrique mesuré dans une configuration pôle-pôle (valeur de potentiel mesuré).Pour traiter les données sismiques (quatre campagnes), les méthodes étudiées dans ce document sont: une tomographie isotrope, une approximation de cosinus, une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de la matrice de rigidité (de l'isotropie transverse horizontale) et une tomographie anisotrope pour un milieu avec une isotropie transverse inclinée. Toutes les méthodes conduisent à la conclusion qu'il y a effectivement une anisotropie présente dans le massif rocheux et confirment la direction de la vélocité maximale parallèle à la direction de la fracturation. Une forte anisotropie de 15 % est présente dans la zone étudiée.Pour traiter les données électriques (deux campagnes), les méthodes étudiées sont une approximation de cosinus et une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de résistivité transversale ρT, de résistivité longitudinale ρL et de l'angle d'orientation du modèle par rapport au système de référence. La présence à 180° de deux maxima de résistivité apparente par rapport à l'azimut de mesure ne peuvent être modélisés par des approches conventionnelles. Une modélisation conceptuelle d'un maillage des fils conducteurs dans un milieu non conducteur a été développé et montre des résultats prometteurs, avec le succès de la modélisation des deux maxima.Des mesures répétées montrent de légères variations de résistivité apparentes et des variations des ondes P. Les variations des ondes S ne sont pas prononcés, ce qui reflète un changement de la saturation en eau. La porosité de 10-15% est estimé à partir des mesures sismiques. / Les résultats présentés dans ce document portent sur les anisotropies sismique et électrique d'un massif calcaire fracturé dans des galeries parallèles souterraines. Les fractures étant orientées verticalement, leurs propriétés sont approximées par isotropie transverse horizontale (HTI). Plusieurs méthodes d'inversion ont été étudié pour traiter les données sismiques (temps d’arrivée des ondes P et S); et pour traiter les données de résistivité électrique mesuré dans une configuration pôle-pôle (valeur de potentiel mesuré).Pour traiter les données sismiques (quatre campagnes), les méthodes étudiées dans ce document sont: une tomographie isotrope, une approximation de cosinus, une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de la matrice de rigidité (de l'isotropie transverse horizontale) et une tomographie anisotrope pour un milieu avec une isotropie transverse inclinée. Toutes les méthodes conduisent à la conclusion qu'il y a effectivement une anisotropie présente dans le massif rocheux et confirment la direction de la vélocité maximale parallèle à la direction de la fracturation. Une forte anisotropie de 15 % est présente dans la zone étudiée.Pour traiter les données électriques (deux campagnes), les méthodes étudiées sont une approximation de cosinus et une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de résistivité transversale ρT, de résistivité longitudinale ρL et de l'angle d'orientation du modèle par rapport au système de référence. La présence à 180° de deux maxima de résistivité apparente par rapport à l'azimut de mesure ne peuvent être modélisés par des approches conventionnelles. Une modélisation conceptuelle d'un maillage des fils conducteurs dans un milieu non conducteur a été développé et montre des résultats prometteurs, avec le succès de la modélisation des deux maxima.Des mesures répétées montrent de légères variations de résistivité apparentes et des variations des ondes P. Les variations des ondes S ne sont pas prononcés, ce qui reflète un changement de la saturation en eau. La porosité de 10-15% est estimé à partir des mesures sismiques.

Les résultats présentés dans ce document portent sur les anisotropies sismique et électrique d'un massif calcaire fracturé dans des galeries parallèles souterraines. Les fractures étant orientées verticalement, leurs propriétés sont approximées par isotropie transverse horizontale (HTI). Plusieurs méthodes d'inversion ont été étudié pour traiter les données sismiques (temps d'arrivée des ondes P et S); et pour traiter les données de résistivité électrique mesuré dans une configuration pôle-pôle (valeur de potentiel mesuré).Pour traiter les données sismiques (quatre campagnes), les méthodes étudiées dans ce document sont: une tomographie isotrope, une approximation de cosinus, une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de la matrice de rigidité (de l'isotropie transverse horizontale) et une tomographie anisotrope pour un milieu avec une isotropie transverse inclinée. Toutes les méthodes conduisent à la conclusion qu'il y a effectivement une anisotropie présente dans le massif rocheux et confirment la direction de la vélocité maximale parallèle à la direction de la fracturation. Une forte anisotropie de 15 % est présente dans la zone étudiée.Pour traiter les données électriques (deux campagnes), les méthodes étudiées sont une approximation de cosinus et une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de résistivité transversale ρT, de résistivité longitudinale ρL et de l'angle d'orientation du modèle par rapport au système de référence. La présence à 180° de deux maxima de résistivité apparente par rapport à l'azimut de mesure ne peuvent être modélisés par des approches conventionnelles. Une modélisation conceptuelle d'un maillage des fils conducteurs dans un milieu non conducteur a été développé et montre des résultats prometteurs, avec le succès de la modélisation des deux maxima.Des mesures répétées montrent de légères variations de résistivité apparentes et des variations des ondes P. Les variations des ondes S ne sont pas prononcés, ce qui reflète un changement de la saturation en eau. La porosité de 10-15% est estimé à partir des mesures sismiques.

Le principal objectif de ce travail est le développement, l'application et la validation de nouveaux traceurs hydrogéologiques sur la base de signatures spectrales spécifiques des matières organiques dissoutes (MOD) et du suivi de leur évolution dans les hydrosystèmes karstiques méditerranéens. Pour cela, de nombreuses analyses hydrodynamiques et hydrochimiques ont été réalisées, et plusieurs méthodes de caractérisation de la MOD ont été utilisées. Le travail s'est alors fondé sur l'étude de différents niveaux du bassin expérimental de la Fontaine de Vaucluse (sol, épikarst et zone non saturée), grâce à un accès privilégié rendu possible par le Laboratoire Souterrain à Bas Bruit de Rustrel-Pays d'Apt (Vaucluse, France). Un bon traceur hydrogéologique requiert des caractéristiques particulières (solubilité dans l'eau, stabilité, tendance faible à l'adsorption, limite de détection faible et des variations bien identifiables). Les composés organiques correspondant le mieux à ces critères sont des substances dérivées de la décomposition de la lignine. Leurs différentes caractéristiques et leurs faibles concentrations dans l'hydrosystème poussent à utiliser des modes de détection suffisamment sensibles tels que la spectrométrie de fluorescence et la Résonance Paramagnétique Électronique (RPE). L’application conjointe des analyses en RPE et en spectrofluorimétrie sur les lixiviats de sols et les eaux souterraines du LSBB permet de proposer deux protocoles de traçages hydrogéologiques. Le premier concerne la détermination des zones superficielles de sol appartenant au bassin d’alimentation des écoulements suivis. Le second, s'appuyant notamment sur l'évolution et la dynamique des MOD dans l'hydrosystème, correspond au développement de traceurs quantitatifs du temps de transit de l'eau pouvant s’appliquer de manière évènementielle par le marquage d’une période particulière de leur variation, ou en continu permettant ainsi le suivi des variations du temps de transit. Enfin, le rapprochement des résultats obtenus par ces nouveaux traceurs organiques avec ceux d'un traceur isotopique (l'oxygène 18) a permis de tester la validité des calculs du temps de transit, mais aussi de démontrer leur complémentarité / The aim of this study is to develop, apply and validate new hydrogeological tracers, based on specific spectral fingerprints of Dissolved Organic Matter (DOM), and a monitoring of their evolution in Mediterranean karstic hydrosystems. Many hydrodynamic and hydrochemical analyses were made, and several DOM characterization methods were used. The basis of this work is therefore the study of different levels of Fontaine de Vaucluse experimental basin (soil, epikarst and unsaturated zone), thanks to a privileged access : the Low-Noise Underground Laboratory of Rustrel - Pays d’Apt (Vaucluse, France). Good hydrogeological tracers require specific characteristics (water solubility, stability, a low tendency to adsorption, a low detection limit and identifiable variations). Organic compounds, corresponding to these criterions at best, are lignin-derived compounds. Their different characteristics and their low concentrations in the hydrosystem drive us to use sensitive enough detection methods, such as fluorescence spectrometry and Electronic Paramagnetic Resonance (EPR). Linked application of EPR and spectrofluorimetric analyses of leachate soil and LSBB groundwater allow to develop two hydrogeological tracing protocols. The first one concerns the determination of upper soil zone, which belongs to monitoring flow catchment area. The second one, using the evolution and the dynamics of DOM in the hydrosystem, corresponds to the development of transit time qualitative tracers, which apply either in an event-driven way by the marking of specific period of their variations, or continuously allowing the monitoring of transit time variations. Finally, the comparison between the results obtained by these new organic tracers and by an isotopic tracer (18 Oxygen) allows to validate transit time calculations and prove their complementarity

Le principal objectif de ce travail est le développement, l'application et la validation de nouveaux traceurs hydrogéologiques sur la base de signatures spectrales spécifiques des matières organiques dissoutes (MOD) et du suivi de leur évolution dans les hydrosystèmes karstiques méditerranéens. Pour cela, de nombreuses analyses hydrodynamiques et hydrochimiques ont été réalisées, et plusieurs méthodes de caractérisation de la MOD ont été utilisées. Le travail s'est alors fondé sur l'étude de différents niveaux du bassin expérimental de la Fontaine de Vaucluse (sol, épikarst et zone non saturée), grâce à un accès privilégié rendu possible par le Laboratoire Souterrain à Bas Bruit de Rustrel-Pays d'Apt (Vaucluse, France). Un bon traceur hydrogéologique requiert des caractéristiques particulières (solubilité dans l'eau, stabilité, tendance faible à l'adsorption, limite de détection faible et des variations bien identifiables). Les composés organiques correspondant le mieux à ces critères sont des substances dérivées de la décomposition de la lignine. Leurs différentes caractéristiques et leurs faibles concentrations dans l'hydrosystème poussent à utiliser des modes de détection suffisamment sensibles tels que la spectrométrie de fluorescence et la Résonance Paramagnétique Électronique (RPE). L'application conjointe des analyses en RPE et en spectrofluorimétrie sur les lixiviats de sols et les eaux souterraines du LSBB permet de proposer deux protocoles de traçages hydrogéologiques. Le premier concerne la détermination des zones superficielles de sol appartenant au bassin d'alimentation des écoulements suivis. Le second, s'appuyant notamment sur l'évolution et la dynamique des MOD dans l'hydrosystème, correspond au développement de traceurs quantitatifs du temps de transit de l'eau pouvant s'appliquer de manière évènementielle par le marquage d'une période particulière de leur variation, ou en continu permettant ainsi le suivi des variations du temps de transit. Enfin, le rapprochement des résultats obtenus par ces nouveaux traceurs organiques avec ceux d'un traceur isotopique (l'oxygène 18) a permis de tester la validité des calculs du temps de transit, mais aussi de démontrer leur complémentarité