L’Afrique de l'Ouest est sujette à une variabilité climatique importante et elle connaît une longue période de sécheresse depuis les années 1970 dont l’impact sur les ressources en eau souterraine reste peu étudié. La présente étude porte sur le bassin versant du Kou au Burkina Faso qui abrite une grande réserve d’eau souterraine à l’origine de sources exceptionnelles. Cette ressource, essentielle pour les besoins de la population, présente pourtant des signes d'épuisement qui soulèvent une question essentielle : quelle est la part du climat et celle des prélèvements dans l’épuisement de cette ressource ? Pour répondre à cette question, nous avons utilisé des méthodes d’analyse statistique et la modélisation hydrogéologique. Les méthodes statistiques ont été appliquées aux chroniques de données climatiques, hydrométriques et piézométriques. Elles ont permis de décrire la réponse des eaux souterraines aux forçages du climat et d’identifier les facteurs explicatifs du comportement de la nappe. La modélisation hydrogéologique a pallié certaines limites de ces méthodes statistiques par une meilleure prise en compte de l’hétérogénéité spatiale de l’aquifère. Elle a en particulier montré que les pompages ont un impact significatif sur la ressource en eau. En somme, les déficits de précipitation, combinés aux pompages en croissance continue, mettent en péril la ressource. La gestion de cette dernière dans une perspective durable devrait donc revêtir une attention particulière. / West Africa is prone to critical climate variability and since the 1970s has been experiencing a long period of drought whose impact on groundwater resources is little studied. This study focuses on the Kou catchment in Burkina Faso which holds a large groundwater resource that lies at the root of exceptional springs. This resource, essential to fulfill the populations’ needs, shows signs of depletion that raises a critical question: what is the share of climate and that of withdrawals in the depletion of the resource? To answer this question, we used statistical analysis methods and hydrogeological modeling. The statistical methods were applied to climatic, hydrometric and groundwater level time series to describe the groundwater’s response to climate forcing and to identify explanatory factors of the groundwater’s behaviour. The hydrogeological modeling overcame some of the statistical methods limitations by better heeding the spatial heterogeneity of the aquifer’s characteristics. Particularly, it revealed that pumping have significant impact on the resource. To sum up, rainfall deficits combined with growing pumping endanger the resource. Hence, sustainable management of the resource should be a critical issue to pay attention to.

La compréhension du cycle de l'eau d'une région où les ressources en eau sont limitées est fondamentale pour assurer une gestion durable de celles-ci, c'est le cas de Galápagos. Cette thèse présente la première étude intégrée du fonctionnement hydrogéologique des aquifères perché de l'île San Cristóbal. Pour ce faire, une approche pluridisciplinaire, fondée sur la mise en place d'un site expérimental sur le bassin versant de Cerro Gato (CG), a été conduite. Pour quantifier les entrées d'eau dans le système hydrologique une approche canopée-transfert hydrique du sol a été utilise. La recharge est principalement affectée par l'altitude en raison de l'effet orographique sur les gradients de pluie et d'évapotranspiration. Néanmoins, à haute altitude, la différence est principalement associée à la couverture végétale en raison de l'apport additionnel fourni par l'interception du brouillard sous la forêt. Une analyse hydrologique montre que les pertes des bassins versants situés à haute-altitude deviennent les entrées d'eau souterraine des bassins versants situés à moyenne altitude (comme CG). Les données du méthode électromagnétique héliportée SkyTEM permettent définir l'extension du bassin hydrogéologique de CG qui est plus grande que le bassin hydrologique. Des approches précédemment décrites, découlent un modèle conceptuel des sources de CG. Ainsi, les eaux souterraines de ces sources sont issues d'un aquifère perché qui s'est formé à la faveur d'une fine couche peu perméable. Ce modèle est testé par simulations numériques qui montrent cette plausibilité. Cette thèse fournit les fondements scientifiques d'une gestion durable des ressources en eau. / The understanding of the hydrogeological functioning in regions such as San Cristobal Island (Galapagos), where water is limited, is fundamental for a suitable management of its resources. This work is the first study of this type on high-level aquifers in San Cristobal using a multi-disciplinary approach, based in the implementation of an experimental site and modeling strategies. For this purpose, a hydrological network was installed in Cerro Gato (CG) and surrounding watersheds. Inputs to the watershed are estimated using the joint modeling of a canopy and soil water transfer. Recharge rates are mainly affected by altitude in mid-elevation watersheds, whereas land cover is the main controlling factor at high-elevation watersheds because of the additional input of fog interception in forests. A hydrological analysis shows that losses from the high-elevation basins become groundwater inputs in the mid-elevation basins, while others have inputs from watersheds at the same altitude, such as CG. The detailed geometry of its hydrogeological watershed is obtained from the dataset of a high resolution AEM SkyTEM survey, which confirms the assumption that its hydrogeological watershed is bigger than its hydrological one. Results from these approaches allow proposing a hydrogeological conceptual model for the springs of CG, where the groundwater flow of springs is fed by a perched aquifer suspended by a low permeability thin layer. This model is tested with numerical simulations, which confirm the plausibility of the existence of a perched aquifer. This thesis provides scientific basis for an effective water management strategy.

Située à l'interface entre les eaux de surface et les eaux souterraines, la zone hyporhéique (ZH) est depuis maintenant plusieurs décennies considérée comme une zone cruciale pour la préservation des milieux aquatiques. Elle constitue souvent un indicateur fiable de la bonne qualité des eaux et une niche écologique primordiale pour de nombreuses espèces. Mais elle est aussi le lieu d'interaction entre deux masses d'eau de signature différente, ce qui conduit à la formation d'un milieu extrêmement fragile et siège d'un grand nombre de réactions biogéochimiques. L'objectif principal de cette thèse est de parvenir à une meilleure compréhension de la dynamique des échanges au sein de la zone hyporhéique. L'approche de cette problématique s'est faite sous un aspect innovant en couplant une démarche hydrogéologique "classique" à l'aide de mesures hydrodynamiques et géochimiques, et l'utilisation de la tomographie de résistivité électrique (ERT). Plusieurs campagnes de terrain ont été menées sur la rivière Essonne, choisie comme lieu d'expérimentation. Différents outils de prélèvement et/ou de mesure ont été mis en place et un grand nombre de mesures à différentes périodes de l'année ont été réalisées. Des expériences assez techniques et innovantes de suivi d'un abaissement et relèvement de barrage, ainsi qu'un traçage artificiel au sel ont pu être effectués grâce à la collaboration avec le syndicat chargé de la gestion et l'aménagement d'une partie du réseau hydrographique de l'Essonne (SIARCE). En parallèle avec cette étude expérimentale, une maquette numérique 3D de la zone d'étude a été réalisée à l'aide du logiciel HydroGeoSphere. Des tests de sensibilité ont permis d'identifier les paramètres hydrodynamiques les plus importants et de quantifier leur impact sur la formation et l'évolution de la zone hyporhéique. Finalement, les premières simulations des expériences menées sur le terrain ont permis de confronter l'approche expérimentale et l'approche théorique. / The Hyporheic Zone (HZ) is located at the interface between surface water and groundwater. For several decades it is considered as a hotspot for the development of a rich aquatic environment in rivers. Its system is often considered as a reliable indicator for water quality and a primary ecological niche for many species. From a hydrological point of view, it is also the place of interaction between two distinct water bodies with different geochemical signatures. This place of mixing forms a very fragile equilibrium where many biogeochemical reactions can occur. The objective of this thesis is to reach a better understanding of mixing and water fluxes in a dynamic context within the hyporheic zone. An innovative method was used by coupling a "classic" hydrogeological approach with hydrodynamic and geochemical measurements with Electrical Resistivity Tomography (ERT). Several field campaigns were done on the Essonne river as experimental site. A large number of measurements were done at various periods of the year and field equipment for water sampling and measurements were installed during these three years. Technical and innovative experiments were conducted such as a dam lowering and rising and an artificial salt tracer test in collaboration with the federation in charge of organization and management of the Essonne network. Finally, a 3D-model of the studied area was built with the HydroGeosphere software. The main hydrodynamic parameters have been tested in order to understand their impact and their variation in a static or dynamic environment on the hyporheic system and its development. In addition, field experiments were reproduced to compare the experimental and theoretical approach.

Le territoire burkinabè est constitué à plus de 80% de formations de socle. Ces formations renferment des ressources en eau souterraine qui constitue la seule source en eau potable de certaines populations. Cependant, le quart des forages réalisés en vue d’exploiter ces eaux sont négatifs. Une bonne compréhension du fonctionnement hydrogéologique et hydrodynamique des aquifères de socle qui sont très hétérogènes peut aider à résoudre ce problème. L’objectif de la thèse est de proposer un modèle conceptuel hydrogéologique en milieu de socle fracturé soudano-sahélien en combinant différentes méthodes hydrogéophysiques et hydrogéologiques classiques. Le site expérimental de Sanon (environ 14 km2) a été choisi comme site d’étude en raison de sa géologie qui est représentative de certaines formations de socle de l’Afrique de l’Ouest. Les résultats indiquent que le profil d’altération de Sanon est cohérent avec d'autres modèles proposés ces dernières années, suggérant que la conductivité hydraulique des aquifères de socle est une conséquence des processus d'altération plutôt que de la fracturation tectonique. Concernant les propriétés hydrodynamiques, une corrélation linéaire positive d’environ 80% a été obtenue entre les transmissivités et les paramètres RMP. La présence d’un dôme piézométrique dans la vallée centrale a été observée en toute saison. La recharge se fait directement au niveau de ce dôme piézométrique (11 et 13% de la pluie annuelle) et l’alimentation des zones périphériques se ferait plutôt par des transferts latéraux. L’ensemble des informations géologiques et hydrogéologiques a permis de faire une simulation numérique qui reproduit l’écoulement observé. / More than 80% of Burkina Faso subsoil is composed of hard rocks. These hard rocks contain ground water resources which are the only drinking water for some people. Yet, the quarter of the boreholes drilled to exploit these waters are negative. A good understanding of the hydrogeological and hydrodynamic functioning of hard rock aquifers which are very heterogeneous can help to solve that problem. The objective of the thesis is to propose a hydrogeological conceptual model of hard rock aquifer under Sudano-Sahelian climate by combining different hydrogeophysical and hydrogeological methods. The Sanon experimental site (14 km² catchment area) was chosen as the study site because of its geology, which is representative of some hard rocks in West Africa. The results show that Sanon weathering profile model is consistent with other models proposed in recent years, suggesting that the hydraulic conductivity of hard rock aquifers is a consequence of weathering processes, rather than tectonic fracturing. Concerning hydrodynamic properties, a positive linear correlation of about 80% was obtained between the transmissivities and the MRS parameters. The presence of a piezometric dome in the central valley on all the piezometric maps produced in both low and high waters was observed. As for recharging, it is done directly in the zone of the piezometric dome (between 11 and 13% of annual rainfall) and the feeding of the peripheral zones would be rather by lateral transfers. The set of geological and hydrogeological data allowed the performance of a numerical simulation in steady state which reproduces well the observed flow.

Le bassin crétacé de Tsagaan Els est situé dans le désert de Gobi en Mongolie. L’objectif de cette thèse est de mieux comprendre le fonctionnement hydrogéologique et hydrochimique de ce bassin endoréique afin de faciliter l’implantation de futurs projets miniers. Dans cette optique, un modèle hydrogéologique régional 3D a tout d’abord été développé sous MODFLOW. La calibration à partir des niveaux piézométriques mesurés sur les ouvrages présents au sein des licences et quelques puits nomades, aboutit à plusieurs solutions. Elles présentent chacune des conductivités hydrauliques en accord avec les tests de pompages, et des taux de recharge entre 0.6 et 3.1 mm/an, cohérents avec le climat actuel. Ce modèle a ainsi permis d’étendre nos connaissances piézométriques au-delà des licences et d’effectuer un bilan hydrique régional. Ensuite, l’âge des eaux a été estimé numériquement via ce modèle puis par des datations 14C. Une comparaison des résultats de ces méthodes a indiqué que les solutions présentant des taux de recharges élevés étaient plus probables que celles en présentant de plus faibles à conditions que les écoulements soient restés les mêmes durant les trente derniers milliers d’années. En parallèle, une étude hydrochimique portant sur les ions majeurs, le brome et les isotopes stable de l’eau (18O et 2H) a montré que l’origine de la salinité des eaux était principalement liée à la dissolution d’évaporites superficielles. En dernier lieu, un modèle 3D local, centré sur la playa terminale du bassin a permis de mieux comprendre le fonctionnement actuel et passé du bassin en s’appuyant sur l’évolution spatiotemporelle de la salinité. / The Tsagaan Els basin is located in the Gobi desert, in Mongolia. This thesis aims to better understand the hydrogeology and hydro-chemistry of the basin in order to optimize the future mining projects. To this end, a 3D groundwater flow model of the basin has first been developed under MODFLOW. The latter was calibrated using piezometric network data and levels in a few nomadic wells. Several solutions were found with hydraulic conductivity values in agreement with pumping tests and recharge rates between 0.6-3.1 mm/yr in agreement with the current climate. The model helped us to extend our piezometric knowledge beyond AREVA’s mining licenses and to establish a regional water balance model. Groundwater age was then estimated with this model and with 14C dating. A comparison of the results of these two methods showed that calibrated solutions with higher recharge are more likely than those with lower recharge under the hypothesis that groundwater flow has not changed significantly during the last 30 kyrs. Meanwhile, a hydrochemical study using major ions, bromines and the stable isotopes of water (18O and 2H) showed that groundwater salinity mainly comes from dissolution of surficial evaporites: likely gypsum and halite and possibly thenardite. Even though evaporation is the main driving force of the groundwater flow, this process, when taking place in subsurface, seems to have no real influence on δ18O and δ2H values of the groundwater. Last a 3D local groundwater flow model centered on the terminal playa and based on the spatiotemporal evolution of the salinity, led us to a better understanding of the past and current hydrogeological functioning of the basin.

Le bassin molassique de Carpentras est constitué de formations miocènes sablo-argileuses présentant de nombreuses variations de faciès. Il renferme un réservoir aquifère complexe aux caractéristiques hydrodynamiques intrinsèques médiocres (perméabilité comprise entre 10-4 et 10-7 m.s-1) mais dont la puissance est importante (600 m d'épaisseur). La qualité chimique de l'eau est excellente (naturellement potable) ; l'exploitation est facile et bon marché. Depuis les années 1960, le nombre de forages a augmenté rapidement et sans contrôle. En 2006, on compte des milliers d'ouvrages sur 600 km2. Dans un contexte de rationalisation des usages de l'eau, il est apparu nécessaire de faire un état des connaissances et d'actualiser celles-ci pour construire un schéma de fonctionnement du système aquifère.<br />Ce travail aborde la géologie, l'hydrodynamique et l'hydrochimie du bassin. Chaque élément fait l'objet d'une revue bibliographique critique et est approfondi par l'acquisition et l'analyse de données nouvelles. L'étude géophysique des bordures du bassin par tomographie de résistivité électrique permet de proposer un schéma de l'organisation des corps sédimentaires. Le suivi mensuel du niveau de l'eau dans 51 ouvrages pendant 20 mois révèle l'influence saisonnière très nette des pompages. Enfin deux campagnes de prélèvements et d'analyses des ions caractéristiques des eaux du bassin (100 éch.) ont été menées en 2004 et 2005. La comparaison des résultats avec ceux acquis en 1985 et 1996 permet une approche spatio-temporelle des concentrations (évolution de la contamination nitratée, temps de séjour).<br />On aboutit à un modèle hydrogéologique simple mais complet. La modélisation est alors utilisée comme un outil de synthèse des informations et des hypothèses de fonctionnement. Elle permet ainsi une visualisation générale du fonctionnement de l'aquifère et met en évidence les éléments qui nécessiteraient des investigations supplémentaires pour améliorer la compréhension.

Les tourbières holocènes des marais du Cotentin reposent sur des bassins sédimentaires sableux aquifères (K = 10-3 – 10-4 m.s-1) dont l'exploitation AEP crée des conditions hydrogéologiques variables dans le temps et dans l'espace sur un site donné. Deux années de suivi hydrogéologique et hydrochimique, ainsi qu'une modélisation hydrogéologique du fonctionnement de la zone humide ont permis de montrer que :<br />- Le fonctionnement hydrologique de la tourbière est contrôlé par celui de l'aquifère des sables et l'existence même de la tourbière est directement liée à sa présence ; <br />- Les conditions hydrologiques influencent directement la variabilité spatio-temporelle des conditions redox du milieu et donc les réactions biogéochimiques mises en jeu dans la tourbière ;<br />- La modélisation de la zone humide a mis en évidence la sensibilité et la vulnérabilité de ces écosystèmes à la pression anthropique (pompage) et/ou l'évolution globale du climat.

Les interactions entre la nappe phréatique peu profonde et les réseaux d’assainissement sont en général très complexes. Dans le but d’étudier l’impact des variations du niveau de la nappe phréatique sur les dégradations des conduites, un modèle hydrogéologique numérique tridimensionnel (3D) de la nappe phréatique ayant un aquifère multicouche au sein d’une zone de 83 km2 de l’agglomération bordelaise (Mérignac, Pessac, Talence) a été développé, calibré et validé. Pour ce faire, un modèle conceptuel de cet aquifère multicouche a dans un premier temps été défini sur la base des descriptions lithostratigraphiques des forages. Puis un travail géostatistique a été entrepris pour construire les limites des formations géologiques et les cartes piézométriques issues de campagnes de mesures de hauteurs d’eau de la nappe phréatique. Nous avons par la suite développé une méthodologie originale de transformation de l’information lithologique au droit des sondages en perméabilités, puis en champs de perméabilités en utilisant des méthodes géostatistiques. Dans un second temps, les conditions hydrogéologiques ont été simulées en utilisant le logiciel Visual Modflow ©. Les paramètres hydrogéologiques ont été calibrés manuellement et automatiquement sur la base des niveaux d’eau mesurés lors des campagnes de juillet 2010, octobre 2010 et mars 2011. La validation du modèle qui s’est effectuée en simulant les cotes piézométriques au pas mi-décadaire allant de juin 1993 à mai 2012 montre que les chroniques piézométriques observées sur les 35 piézomètres du domaine étudié sont correctement restituées et la dynamique de la nappe phréatique est reproduite. Dans un dernier temps, les distributions spatiales et temporelles des cotes piézométriques simulées ont permis d’identifier les zones d’ennoiement des réseaux d’assainissements en comparant les cotes piézométriques aux cotes des réseaux. Ces zones d’ennoiement ont été comparées aux principaux secteurs de casses de canalisations. Les secteurs de casses de canalisations sont localisés dans les zones à faible probabilité d’ennoiement des conduites. L’environnement naturel hydrogéologique pris comme paramètre isolé ne constitue donc pas un facteur explicatif au phénomène de casses des canalisations. / The interactions between shallow groundwater and sewer network are generally complex. The aim of this research is to study the impact of the water table level variation on the damages of the buried pipes. For this purpose, a three-dimensional (3D) numerical groundwater flow model of the shallow multi-aquifer of the Bordeaux urban area, within a 83 km2 area (Mérignac, Pessac, Talence), is developed, calibrated and validated. In order to do this, in a first step, a conceptual model was built based on stratigraphic descriptions of the boreholes. A geostastistical study was carried out to build the limits of the geological layers and the piezometrical maps from the measurement campaigns of water table levels. Then, an original methodology to transform the description of lithological units to local hydraulic conductivity values is proposed and the reconstruction of hydraulic conductivity fields at the urban scale was performed based on geostatistical methods. In a second step, hydrogeological conditions were simulated using the finite-difference groundwater flow model MODFLOW-2005 (Visual Modflow ©). The hydrogeological parameters were then calibrated manually and automatically based on water table elevation data measured in July 2010, October 2010 and March 2010. The model validation covered 6890 days, from June 1993 to may 2012. The results showed that the simulated heads are quite accurate and reproduce the main dynamics of the system at 35 piezometers. The spatial and temporal distribution of the heads of the piezometers simulated permitted to identify the flooding zones of the sewer network by comparing the piezometric head with the altitude of the networks. The sectors where the pipes are damaged were localized in the zone wherein there is low probability of flooding. Additionally, the natural hydrologic environment, taken as a sole factor, cannot explain the pipe breakage phenomena.

Les eaux souterraines constituent une réserve d’eau potable non négligeable, leur alimentation se fait en majeure partie par les précipitations, appelée recharge des nappes. Du fait de leur grande importance, la compréhension du fonctionnement de ces ressources en eau est plus que jamais indispensable. Celle-ci passe par l’élaboration de modèles mathématiques. Ces outils nous offrent une meilleure appréhension et une bonne prévision des phénomènes physiques. Les systèmes hydrogéologiques sont généralement très complexes et caractérisés par des dynamiques hydriques très variables dans le temps et dans l’espace. Cette complexité a attiré l’attention de nombreux hydrogéologues et un grand nombre de modèles très sophistiqués a été développé afin de décrire ces systèmes avec précision. Cependant, la prise en compte de la recharge de ces réservoirs reste toujours un défi dans la modélisation hydrogéologique. En effet, le plus souvent, les modèles hydrogéologiques simulent l’écoulement dans la nappe tout en considérant la recharge comme une constante sur le domaine et indépendante du système. De plus, elle est souvent calculée de façon simplifiée. Or, la recharge traduisant la quantité des précipitations atteignant les nappes est une composante hydrologique complexe et variable car elle interagit avec les nappes et dépend des conditions climatiques, du couvert végétal et du transfert de l’eau dans le sol. Ce présent travail vise à intégrer cette recharge variable et complexe aux modèles hydrogéologiques. À cet effet, un modèle couplé a été développé. Une première partie de ce modèle permet de calculer la recharge des nappes en modélisant les interactions précipitations-sol et l’hydrodynamique dans le sol. Cette modélisation a été effectuée en utilisant des modèles conceptuels simples basés sur des lois empiriques (Gardénia, Nash) et des modèles physiques résolvant l’équation de Richards. La recharge ainsi calculée est intégrée à la deuxième partie du modèle simulant l’hydrodynamique dans les nappes décrite par l’équation de diffusivité. Des méthodes numériques précises et robustes ont été utilisées pour résoudre les équations du modèle mathématique : les éléments finis non conformes ont été utilisés pour résoudre l’équation de diffusivité et l’équation de Richards est résolue sous sa forme mixte par une méthode itérative en temps. En somme, ce modèle couplé permet de décrire les variations de niveaux de nappe à partir des données météorologiques connaissant les paramètres caractéristiques de cet aquifère. [...] / Groundwater is the main available water resource for many countries; they are mainly replenished by water from precipitation, called groundwater recharge. Due to its great importance, management of groundwater resources is more essential than ever, and is achieved through mathematical models which offer us a better understanding of physical phenomena as well as their prediction. Hydrogeological systems are generally complex thus characterized by a highly variable dynamic over time and space. These complexities have attracted the attention of many hydrogeologists and many sophisticated models that can handle these issues and describe these systems accurately were developed. Unfortunately, modeling groundwater recharge is still a challenge in groundwater resource management. Generally, groundwater models are used to simulate aquifers flow without a good estimation of recharge and its spatial-temporal distribution. As groundwater recharge rates show spatial-temporal variability due to climatic conditions, land use, and hydrogeological heterogeneity, these methods have limitations in dealing with these characteristics. To overcome these limitations, a coupled model which simulates flow in the unsaturated zone and recharge as well as groundwater flow was developed. The flow in the unsaturated zone is solved either with resolution of Richards equation or with empirical models while the diffusivity equation governs flow in the saturated zone. Robust numerical methods were used to solve these equations: we apply non-conforming finite element to solve the diffusivity equation and we used an accurate and efficient method for solving the Richards equation. [...]

Une augmentation à long terme des réserves de l'aquifère dans la région de Niamey a été mise en évidence en dépit de la sécheresse qui sévit dans cette région depuis 30 ans. Les profondes modifications environnementales depuis 50 ans sont à l'origine d'un accroissement du ruissellement et indirectement de l'augmentation de la recharge de la nappe libre. Le manque de données significatives sur le ruissellement à large échelle et l'absence de relation simple avérée entre la recharge et la pluviométrie a conduit les travaux de caractérisation de la hausse de la nappe phréatique vers une confrontation des modélisations de surface avec les données piézométriques.<br /> Dans cet objectif, une régionalisation du modèle de fine échelle abc-rwf a été réalisée à partir d'un jeu simulé de 50 bassins endoréiques, pour lesquels des relations entres les descripteurs physiographiques et le ruissellement calculé sont apparues significatives. La formulation d'une fonction représentant fidèlement le ruissellement simulé d'un bassin, peut être étendue à un ensemble de bassins par l'adjonction de fonctions prédictives non linéaires paramétrées sur la base de seulement 5 descripteurs physiographiques des bassins. Le modèle de méso-échelle Zarbhy permet ainsi un calcul simplifié du ruissellement des 377 bassins endoréiques identifiés sur les 5000 km² du bassin tout en limitant l'imprécision supplémentaire par rapport au modèle abc rwf.<br />La réalisation du couplage surface-souterrain soulève une insuffisance de la recharge forcée par le ruissellement calculé sur certaines années pour permettre la simulation de la hausse de la nappe dans des gammes de paramétrisation admissibles. Bien que l'extrême variabilité des processus hydrologiques disparaisse en partie à l'échelle régionale, la variabilité intersaisonnière demeure encore difficile à reproduire.<br />Le schéma de la recharge ponctuelle spatialement distribuée est compatible avec la modélisation hydrodynamique réalisée, qui dans son état actuel ne permet cependant pas de quantifier avec précision l'influence de la recharge indirecte sur la dynamique de hausse de la nappe.<br />Face à l'insuffisance de la recharge calculée, une étude combinant géophysique électrique, géochimie et modélisation hydrologique a été entreprise. Elle suggère le lessivage d'une bande minéralisée entre 5 et 10 m de profondeur à l'aplomb des multiples chenaux sableux de la zone d'épandage de Wankama. Cette expérimentation prometteuse ouvre la voie à de nouvelles investigations vers la caractérisation des processus de recharge. L'existence d'une infiltration profonde (>25 m) au moins épisodique, sous les zones d'épandage sableuses peut être envisagée.