La pollution par les nitrates des eaux de surface et des eaux souterraines a suscité une attention mondiale ces dernières années. La recharge des eaux souterraines via d'infiltration dans les zones cultivées est une source importante de la contamination des eaux de surfaces. Les plaines alluviales supportent une activité agricole intensive et subissent une pollution aux nitrates importante. Il a été démontré que les échanges entre les eaux de nappes et la rivière contribuent à la rétention et/ou transformation de l'azote dans le continuum eaux de surface-eaux souterrains. La compréhension et la quantification des processus d'atténuation des concentrations en nitrates se produisant à l'interface nappe-rivière permettraient d'améliorer la connaissance du cycle de l'azote à l'échelle du bassin versant. Les objectifs de cette thèse sont : 1) quantifier les volumes d'eau échangés entre la rivière et l'aquifère alluvial dans les zones de plaine alluviale à l'échelle du bassin versant et 2) quantifier les taux de dénitrification dans les nappes alluviales à l'échelle du bassin versant et évaluer leurs influences sur les flux de nitrates de la rivière. Un échantillonnage sur le terrain ainsi qu'un travail de modélisation ont été effectués pour atteindre ces objectifs. Les campagnes d'échantillonnages sur le terrain ont eu lieu d'avril 2013 à mars 2014 sur le site d'étude de Monbéqui dans la plaine alluviale de la Garonne. Le modèle Soil and Water Assessment Tool (SWAT) qui est largement utilisé à travers le monde a été choisi pour simuler les processus hydrologiques et le cycle de l'azote. Cependant, la simulation des échanges nappe-rivière qui se produisent dans la plaine alluviale n'est pas prise en compte dans la version originale de SWAT. Premièrement, l'équation de Darcy a été introduite pour simuler les échanges nappe-rivière à partir de la structure Landscape Unit (LU). L'influence des crues débordantes sur ces échanges a également été ajoutée dans le modèle. Le modèle modifié, SWAT-LUD, a d'abord été appliqué à un méandre de la plaine alluviale de la Garonne- Monbéqui (25 km²). Ensuite, un module supplémentaire représentant les processus de dénitrification dans les aquifères peu profonds des nappes alluviales a été développé et ajouté au modèle SWAT-LUD. Les flux de nitrates ainsi que les taux de dénitrification de l'aquifère alluvial de Monbéqui ont alors été quantifiés. Dans un deuxième temps, le modèle SWAT-LUD a été appliqué à l'échelle de la plaine alluviale de la Garonne dans son cours moyen (environ 4600 km²) et l'influence de la dénitrification de l'aquifère alluvial sur les flux de nitrates de la rivière a été quantifiée. Finalement, le modèle a été appliqué sur l'ensemble du bassin versant de la Garonne (environ 51 500 km²) et l'influence des échanges nappes-rivières sur le cycle hydrologique du bassin versant a été évaluée. Les résultats ont montré que le modèle SWAT-LUD pouvait représenter de façon réaliste les échanges nappe-rivière ainsi que les taux de dénitrification dans les aquifères alluviaux à différentes échelles. Les flux échangés vont majoritairement de l'aquifère vers la rivière et contribuent pour 65% du volume total échangé. À l'échelle du bassin versant de la Garonne, le volume annuel échangé entre la nappe et la rivière représente environ 5% du débit annuel de la Garonne. A l'échelle du méandre, le taux annuel de dénitrification dans la zone riparienne a été estimé à environ 130 N-NO3-.ha-1.an-1. Près de 40% des nitrates arrivant dans cette zone ont été dégradés via les processus de dénitrification. Dans le cours moyen de la Garonne, les taux de dénitrification annuels dans ces mêmes zones varient entre 55 et 120 kg N-NO3-. ha-1. an-1. / Nitrate pollution in surface water and groundwater systems has attracted worldwide attention in recent decades. Recharged groundwater infiltrating through cultivated fields is an important source of the nitrate contamination of surface water. As alluvial plains support intensive agricultural activities, they often suffer from groundwater nitrate pollution. The exchanges between surface water and groundwater (SW-GW) were proved contributing to nitrate retention and/or transformation in the land-surface water continuum. The understanding and quantifying of nitrate attenuation processes occurring at the surface-groundwater interface would enhance understanding of nitrogen cycling at the catchment scale. The objectives of this thesis were: 1) quantifying the exchanged SW-GW volume in the floodplain area at the catchment scale and 2) quantifying the shallow aquifer denitrification rate in the floodplain area at the catchment scale and evaluating its influence on the nitrate flux in the river. Field sampling and modelling study were conducted to achieve these objectives. Monthly field work campaigns were carried out from April 2013 to March 2014 at the Monbéqui site in the Garonne river floodplain. The Soil and Water Assessment Tool (SWAT) which has been successfully applied all over the world was chosen to simulate the hydrologic processes and nitrogen cycle. However, the simulation of the water exchange between river and groundwater occurring in the floodplain area was not simulated in the original SWAT model. Firstly, the Darcy's equation was implemented to simulate SW-GW exchanges based on the Landscape Unit (LU) structure in the floodplain area. The influence of flooding on the water exchange was also introduced to the model. The modified model was called SWAT-LUD and was applied to a meander of Garonne floodplain - Monbéqui (around 25 km2). Then, another module representing the denitrification processes in the floodplain shallow aquifer was developed and added to the SWAT-LUD model. The nitrate flux and shallow aquifer denitrification rates in Monbéqui was then quantified. Afterwards, the SWAT-LUD model was applied to the middle floodplain section of the Garonne River (around 4 600 km2) and the influence of shallow aquifer denitrification on the nitrate flux in the river was quantified. Lastly, the model was applied to the entire Garonne catchment (around 51 500 km2) and the hydrologic influence of SW-GW exchanges was evaluated. The results showed that the SWAT-LUD model could satisfactorily represent the SW-GW exchanges and shallow aquifer denitrification rate at different spatial scales. The main water flow direction is from the shallow aquifer to the river, with water flowing in this direction accounted for around 65% of the total exchanged water volume. In the Garonne catchment, the annual total exchanged water volume represented around 5% of the total discharge volume of the Garonne river. For the Monbéqui site, the simulated annual denitrification rate in the riparian zone was around 130 kg N-NO3-ha-1y-1. Around 40% of the nitrate input in this zone was degraded through denitrification. In the middle floodplain section, the annual denitrification rate in the near bank zone ranges from 55 to 120 kg N-NO3-ha-1y-1.

Dans le cadre du programme international AMMA (Analyse Multidisciplinaire de la Mousson Africaine), la haute vallée de l'Ouémé (centre ouest du Bénin), d'une superficie de 10.000 km2, est instrumentée (pluie, débit) depuis 1997. Le sous-bassin de la Donga (586 km2) est le cadre d'observations intensives précisant les processus hydrologiques essentiels et quantifiant les termes du bilan.<br />La recherche entreprise vise (i) à déterminer le fonctionnement hydrologique du bassin de la Donga par confrontation des données de précipitation, d'écoulement de surface, de nappe et d'évapotranspiration potentielle, (ii) à caractériser les écoulements par traçage naturel géochimique (paramètres physiques, éléments majeurs et traces) et isotopique (isotopes stables de l'eau : d18O et d2H) afin de hiérarchiser les processus majeurs mis en jeu. Le croisement des informations hydrodynamiques et hydrogéochimiques permet de reconstituer un schéma cohérent du fonctionnement hydrologique sur le bassin. Le caractère temporaire des eaux de surface, l'asynchronisme entre le tarissement des rivières et la vidange des nappes ainsi que la très faible minéralisation des écoulements traduisent une origine superficielle des débits et une déconnexion apparente de la nappe libre du réseau hydrographique.<br />En accord avec les processus mis en évidence, une décomposition de l'hydrogramme en deux composantes (ruissellement rapide sur surface saturée et flux de subsurface) a été retenue pour 2003 et 2004, deux années à pluviométrie et écoulement contrastés.

Cette thèse s'inscrit dans le développement de la plateforme EauDyssée de modélisation intégrée des hydrosystèmes régionaux, au sein du bassin pilote de la Seine. L'objectif principal est de contribuer à une meilleure simulation des tirants d'eau à l'échelle régionale afin d'améliorer la simulation des interactions nappe-rivière et de mieux quantifier les niveaux piézométriques dans les aquifères. La première partie de la thèse vise à évaluer la sensibilité d'un modèle hydraulique à la précision de la description géomorphologique des lits pour identifier le meilleur compromis entre parcimonie et réalisme et identifier les facteurs morphologiques les plus importants pour obtenir une simulation satisfaisante des tirants d'eau à l'échelle régionale. Cette étude est menée sur le Serein (affluent de l'Yonne), entre les stations limnimétriques de Dissangis et Beaumont, dans un bief bien renseigné (20 sections transversales sur 89 kms). Débits et tirants d'eau sont simulés par le modèle hydraulique HEC-RAS (équations de Saint-Venant 1D), en fonction des apports latéraux simulés par le modèle régional EauDyssée. Les résultats de cette étude montrent qu'un modèle 1D type Saint-Venant n'est pas adapté à la simulation des écoulements à l'échelle régionale. Nous avons donc développé une méthode de changement d'échelle originale, dans laquelle la modélisation fine des processus hydrauliques à haute résolution permet d'améliorer la représentation des profils d'eau en rivière et les interactions nappe-rivière simulées à l'échelle régionale par le modèle intégré EauDyssée. Cette méthodologie de changement d'échelle a été validée dans un sous bassin versant de l'Oise d'une superficie de 4500 km2, pour la période 1990-1995. Nous avons utilisé HEC-RAS pour la modélisation hydraulique d'un tronçon de l'Oise de 188 km, où 420 sections transversales sont disponibles. Le modèle permet d'interpoler des courbes de tarage simulées tous les 200m en moyenne. Ces courbes de tarage sont ensuite projetées sur les mailles rivière du modèle régional EauDyssée (résolution de 1 km), où elles permettent de simuler la fluctuation du niveau d'eau en fonction du débit à l'échelle régionale par EauDyssée. La cote de la surface libre de la rivière définissant sa charge hydraulique, ces fluctuations influencent alors les échanges entre les mailles rivière et les nappes, qui dépendent des gradients de charge verticaux entre rivière et nappe (loi de Darcy). Ce travail montre l'intérêt de l'approche pour mieux évaluer les interactions nappes-rivières à l'échelle régionale avec un faible coût de calcul. Il offre des perspectives intéressantes pour simuler des processus jusque là négligés par le modèle EauDyssée : élimination de nitrate dans les zones humides qui sont souvent situées à la zone de contact entre les nappes souterraines et la rivière, ou l'impact du changement climatique sur le fonctionnement des hydrosystèmes et plus particulièrement sur l'élimination ou le relarguage de polluants par des processus biogéochimiques, ainsi que de mieux estimer les risques d'inondation à l'échelle régionale.

Les interactions entre les eaux souterraines et les eaux de surface sont complexes et jouent un rôle prépondérant dans le fonctionnement des hydrosystèmes, tant en termes quantitatifs que qualitatifs. Ces deux compartiments interagissent à travers différentes interfaces emboitées depuis l’échelle locale, notamment la zone d’interface rivière-nappe, appelée zone hyporhéique (ZH) jusqu’à l’échelle régionale. La ZH se définit comme un buffer entre les eaux souterraines et les eaux de surface. Les variations des flux d’eau et de chaleur au niveau de cette zone modulent les processus biogéochimiques. Le but de cette thèse est d’améliorer la compréhension du régime thermique et hydrologique de la ZH et de sa variation spatiale. Cette problématique est abordée par l’utilisation de stations de mesures haute fréquence des échanges nappe-rivière (MOLONARI) sur le bassin des Avenelles (sous bassin de la Seine (46 km2)) et de la modélisation. Le suivi haute fréquence a permis de réaliser une analyse des données expérimentales présentant le régime hydro-climatique du bassin ainsi que l’évolution spatiale et temporelle du fonctionnement des différents compartiments de l’hydrosystème. Un cas synthétique a été réalisé afin de caractériser les différents facteurs contrôlant le régime thermique de l’interface nappe rivière et la variation de stocks d’énergie au sein de la ZH. La démarche de modélisation suivie pour le cas synthétique fournit un cadre d’analyse des données des cinq stations MOLONARI. Les données acquises sur ces stations ont été utilisées afin de déterminer les propriétés hydrodynamiques et thermiques de la ZH ainsi que des couches géologiques sous-jacentes par inversion à l’aide du modèle METIS couplé à un script de balayage de l’espace des paramètres afin de quantifier les flux d’eau et de chaleur à l’interface nappe-rivière le long d'un corridor hydraulique. Les résultats du modèle mettent en lumière la variabilité spatio-temporelle des échanges de chaleur au niveau des cinq stations MOLONARI. La quantification des flux de chaleur a fourni un nouvel éclairage des interactions entre les eaux de surface et les eaux souterraines à l’interface nappe-rivière du bassin des Avenelles. / The river-aquifer interactions are complex and play a preponderant role in hydrosystems functioning, in both qualitative and quantitative terms. These two compartments interact through various nested interfaces from the local scale in particular the river aquifer interface called the Hyporheic Zone (HZ) to the regional scale. The HZ acts as a buffer between the stream and the groundwater. The water and heat fluxes variation at this zone modulates the biogeochemical processes.The aim of this thesis, is to improve the understanding of the thermal and hydrological regime of the HZ and its spatial variation. This problem is approached by the use of high frequency measurement stations of river-aquifer exchanges (LOMOS) in the Avenelles basin (sub-basin of the Seine basin (46 km2)) and by modeling. The high frequency monitoring allowed to realize an analysis of the experimental data characterizing the hydroclimatic regime of the Avenelles basin as well as the spatial and temporal evolution of the hydrosystem various compartments functioning.A synthetic case was realized to characterize the main factors controlling the thermal regime at the river aquifer interface and the energy storage variation within the ZH. The modeling approach followed by the synthetic case provides an analysis framework of the five LOMOS data.The acquired LOMOS data were used to determine the hydrogeological and thermal properties of the HZ and of the underlying aquifers by inversion by means of a 2D finite element thermo-hydrogeological model (METIS) coupled with a parameters screening script to quantify water and heat fluxes through the stream – aquifer interface over along the stream network.The model results highlight the spatiotemporal variability of the heat exchanges at the five LOMOS. The quantification of heat fluxes provided a new lighting of the stream-aquifer interactions of the Avenelles basin.

Compte tenu de l’évolution démographique et climatique planétaire, la gestion de la ressource en eau constitue un défi majeur auquel la communauté internationale devra faire face au cours du XXIème siècle. A cet effet, l'identification de la continuité hydrique entre les eaux de surface et les eaux souterraines permet l'introduction de la notion de gestion intégrée de la ressource. L'application de ce principe de gestion au bassin de la Seine, à travers l'estimation des échanges nappe-rivière, est rendue possible par la mise en pratique du concept d'interface nappe-rivière emboitées au sein de travaux de modélisation. Pour cela une procédure de modélisation multi-échelles peut être mise en place. Elle vise à intégrer des informations locales au sein de modélisation à une échelle supérieure. Dans ce mémoire, une procédure de modélisation multi-échelles est mise en œuvre. Ce protocole de modélisation est initié par une estimation des flux d'eau régionaux au sein de l'hydrosystème Seine. La cohérence globale de ces flux est garantie par le développement d'une méthodologie de calibration de modèles couplés en deux étapes. Ensuite les informations locales, que sont les hétérogénéités de la plaine alluviale de la Bassée et de la représentation des interfaces nappe-rivière du réseau secondaire, sont intégrées au modèle régional par une procédure de modélisation emboitée et de changement d'échelle des paramètres hydrauliques. La mise en place de cette procédure a finalement permis l'estimation fine des échanges nappe-rivière sur la quasi-totalité (83%) du réseau hydrographique naturel du bassin de la Seine et ainsi de répondre aux recommandations de gestion intégrée de la ressource faites par la directive cadre sur l'eau. / Given the current climate and anthropogenic evolutions, water management becomes one of the greatest challenges of the 21st Century. For that purpose, by identifying hydraulic continuity between surface and subsurface water, the concept of integrated water management can be introduced. In this work this management concept is applied on the Seine basin by quantizing hydrological processes occuring at the nested stream-aquifer interface. The implementatin of the nested interface concept can bedone through multi-scale modeling. This modelling procedure, aimed at embody the local characteristics of the interfaces (such as structural or hydrodynamic heterogeneities) in large scale models. A multi-scale modelling procedures is applied to the regional Seine basin model (70000 km²) in order, to study the hydrodynamic behaviour of the Bassée alluvial plain, and to quantify the stream-aquifer exchanged fluxes at the basin scale. The modelling protocol is initiated with regionals fluxes estimation over Seine hydrosystem. Regional fluxes consistency are assured by a two-step calibration procedure of fully coupled models. Then, the local characteristics of the Bassée alluvial plain, are implemented in the regional model by nested modelling methodology associated with upscaling procedure of hydraulics properties. Finally, the multi-scale modelling procedure lead to quantify distributed stream-aquifer exchanged water fluxes over 83% of the natural river network of the Seine basin, and thus, achieve to answer the integrated water resources management recommandations of the water framework directive.

La prise en compte de l'hétérogénéité spatiale des bassins versants est nécessaire pour répondre aux divers enjeux liés à la gestion des ressources en eaux, tant au niveau quantitatif que qualitatif. Les modèles hydrologiques spatialisés permettent de prendre en compte l'hétérogénéité des surfaces continentales et des données d'entrée dans la modélisation du cycle de l'eau. La revue bibliographique réalisée dans la première partie a permis d'aborder différents aspects de la problématique de la prise en compte de l'hétérogénéité spatiale dans les modèles hydrologiques. Cette revue a débouché sur la proposition d'une méthode de découpage spatiale suivant trois niveaux emboîtés, de manière à adapter la description de l'hétérogénité spatiale des bassins versants à la question posée. Le premier niveau consiste en un découpage en sous-bassins (REWs) organisés autour du réseau hydrographique. Le second niveau de discrétisation (hydro-paysages) permet de décrire l'hétérogénéité spatiale des REWs. Nous avons proposé un cadre méthodologique pour le découpage en hydro-paysages. Le troisième niveau permet d'obtenir les mailles finales de la modélisation après redécoupage éventuel des hydro-paysages en fonction de contraintes numériques. Un exemple sur le Haut-bassin de la Saône permet d'illustrer l'ensemble de la démarche. La méthodologie proposée est associée à une modélisation flexible permettant d'adapter le découpage spatial à l'échelle des données disponibles et à l'échelle des processus hydrologiques à modéliser. Dans la deuxième partie de la thèse, nous avons présenté les outils utilisés pour réaliser des modélisations hydrologiques ‘à la carte', basés sur la plate-forme de modélisation LIQUID qui a servi de support à notre travail. Afin de prendre en compte dans la modélisation, les mailles non-structurées issues du découpage spatial proposé, nous avons développé un modèle de simulation des tranferts latéraux dans la zone saturée. Il résout l'équation de Boussinesq 2D sur un maillage irrégulier mais conforme. Différents tests ont été réalisés pour valider le modèle. Nous avons enfin abordé, sur quelques exemples simples, la problématique du couplage de différents modules en partant de ce module saturé et en étudiant le couplage avec les écoulements dans la rivière et la zone non saturée.

La compréhension des transferts sédimentaires depuis les zones " source " jusqu'aux zones de dépôt des sédiments constitue un enjeu majeur pour la lutte contre l'érosion des sols et la qualité des eaux. Portant sur de petits bassins versants côtiers de l'Ouest du bassin parisien, l'intérêt de ce travail est de présenter une démarche intégrée de la dynamique érosive à l'échelle du bassin versant, depuis l'analyse de la stabilité structurale des sols jusqu'au suivi haute fréquence des transferts de matières en rivière. L'objectif est d'identifier les facteurs contrôlant le ruissellement et l'érosion sur ces sites et de proposer une estimation des apports fluviatiles de matières en suspension (MES) à la Manche orientale. Les résultats obtenus soulignent que malgré des sols limoneux peu structurés générant des taux d'érosion à la parcelle importants, les bilans d'érosion enregistrés en rivière sont parmi les plus faibles du monde (10,1 t.km-2.an-1 pour le bassin de l'Yères, 1,3 t.km-2.an-1 pour le bassin du Dun et 6,6 t.km-2.an-1 pour le bassin de la Ganzeville). La comparaison de ces données avec celles de 214 bassins versants mondiaux montre que ces faibles bilans sont expliqués par la taille des bassins étudiés (< 1000 km2), le relief modéré et le faible débit spécifique. Régionalement, ce dernier facteur est contraint par un fonctionnement hydrologique karstique, qui est à l'origine non seulement du déficit de la balance des écoulements de surface mais aussi d'un stockage intrakarstique de matériaux fins issus de l'érosion des sols. Les résultats mettent en évidence des échanges nappe-rivière très importants sur les sites d'étude : la contribution des eaux souterraines au débit de surface est de l'ordre de 90 %. Concernant les apports fluviatiles de MES à la Manche orientale, l'extrapolation des bilans d'érosion obtenus à l'ensemble des fleuves côtiers du littoral haut-normand conduit à une estimation du flux exporté à la mer de 48 500 t.an-1. Issu de l'érosion des bassins versants côtiers, ce flux de MES apparait très faible tant à l'échelle mondiale que régionale où le flux de MES de la Seine atteint en moyenne 715 000 t.an-1 à Poses en amont de l'estuaire.

Les interactions entre les processus de surface et de sub-surface sont des composantes clés du cycle hydrologique que les modèles hydrologiques doivent représenter pour obtenir des prédictions cohérentes et précises dans un contexte de gestion durable de la ressource en eau. Les modèles hydrologiques intégrés qui décrivent de façon physique les processus et leurs interactions sont de conception récente. La plupart de ces modèles s‘appuient sur l’équation de Richards 3D pour décrire les processus d’écoulement souterrain. Cette approche peut être problématique compte tenu de contraintes importantes sur le maillage et sur la résolution numérique. Ce travail de thèse propose un modèle hydrologique intégré qui s’appuie sur approche innovante à dimension réduite pour simplifier les écoulements de surface et souterrains d'un bassin versant. Les différents compartiments du modèle sont d’abord testés indépendamment puis couplés. Les résultats montrent que l’approche proposée décrit précisément les processus hydrologiques considérés tout en améliorant de façon significative l’efficacité générale du modèle. / Interactions between surface and subsurface flow processes are key components of the hydrological water cycle. Accounting for these interactions in hydrological modelsis mandatory to provide relevant and accurate predictions for water quality and water resources management. Fully-integrated hydrological models that describe with aphysical meaning the hydrological processes and their interactions are recent. Most of these models rely upon the resolution of a 3D Richards equation to describe subsurface flow processes. This approach may become intractable because of the heavy constrains on both meshing and numerical resolution. This PhD proposes a new integrated hydrological model on the idea of dealing with dimensionally reduced flow in both the surface and sub-surface compartments of a watershed. The different compartments of the model are first tested independently and then coupled. The results show that the proposed approach allows for a proper and precise depiction ofthe hydrological processes enclosed in the model while providing significant gain incomputational efficiency.

La modélisation hydraulique avec modèle déterministe est une méthode largement utilisée. Cependant,lamodélisation est un ptocessus complexe, notamrnent pour les aquifères où la quantité et la qualité desdonnées ne sont pas satisfaisantes. Etantune des sources pdncipales de I'eau douce dans la basse valléedu Var, Côte d',{.zur,Frarrce,la nappe libre de lavallée est menacée parla pénurie et la pollution. Maþéle grand nombre d'études effectuées dans cette zone, la dynamique des écoulements souterrains esttoujours patiellement inconnue. Par conséquent, la métropole Nice Côte d'Azur a besoin de développerun système d'aide à la décision (SÂD) à base des modèles numériques afin d'assurer une gestion plusefficace de l'eau souteraine. Un modèle numérique est développé avec FEFLO\ø en tenant compte desprécipitations, de l'é:vapoftanspiration, du pompage de l'eau souterraine, et des échanges nappe-rivière.Le volume d'eau pompée pour l'usage agdcole et le taux d'échange du lit mineur du Var ont été calibrés.Le modèle a été. vabdé. par une simulation de '1.266 jours. Le modèle est utilisé pour simuler les scénatüdes événements d'inondation et de sécheresse, les scénadi de pollution et l'intrusion d'eau de mer dans lavallée. Une première conception du système d'aide à la décision est présentée comme le demier exemplede l'application du modèle. Une interface de couplage est développée en Java et sert à échanger lesdonnées entre le modèle souterrain et le modèle à surface libre élaboré avec MIKE21FM. / Groundwater modeling with deterministic model is a complicated process, especially in complex aquiferswhere the quantity and the quality of the measuted data arc not satisfying. The unconfined alluvialaquifet is the main water resource in the lower valley of Vat river, Ftench Riviera, but it faces a thteat ofshortage and pollution. Despite numerous previous studies, the dynamics of the gtoundwater flow in thealluvial aquifer and the characteristics of the rivet-aquifer exchanges ¿re s':ll partially unknown.Therefore the local u/ater management service requires a decision support system PSS) based onnumerical models to ensure a better groundwater management. A hydraulic model is set up withFEFLO!ø software by considedng ptecþitation, evâpotranspiration, gtoundwatet exftacdon and rivetaquiferexchanges. The non-documented groundwater exttaction fot agticultural use and the transferrates in the dverbed along the river have been calibtated. The model has been validated with asimulation of 7266 days. The model is applied to simulate the scenarios of flood and drought events, thepollution events in the unconfined aquifer in the valley and the seawater inrusion in the estuary of Yarriver. These case studies contributes to increase the knowledge of the aquifet. A fust conception of theDSS tool is presented as the last example of model application. A coupling interface is developed thanksto aJava which enables an automatic exchange of data between the groundwater flow model and thesurface wâter flow model built with MII(E2ltr}'/. More tests should be done to validate the couplinginterface

Des analyses de sensibilité et des estimations de paramètres sont étudiées sur deux études de cas de transfert de masse en milieu poreux. La première partie est consacrée à la sensibilité des écoulements souterrains dans une modélisation des échanges drain-aquifère pour mettre en évidence les différences entre les deux méthodes de discrétisation mises en œuvre. La seconde partie est dédiée à la modélisation de l’écoulement en milieu poreux variablement saturé dans une zone humide artificielle, au calage des paramètres du modèle de van Genuchten-Mualem et à l’évaluation de son efficacité à reproduire des données piézométriques collectées sur le site de l’Ostwaldergraben. La variabilité temporelle des paramètres hydrodynamiques, incluant l’effet d’hystérésis, montre que ceux de la couche active du filtre changent au cours du temps. Ces deux études sont conduites à l’aide de la différenciation automatique. / Sensitivity analyses and parameter estimation are applied to mass transfer in porous media for two remediation facilities. The first part is devoted to the sensitivity analysis of groundwater flows in a modeling of drain-aquifer exchanges to highlight the differences between the two implemented methods of discretization. The second part is dedicated to the modeling of the flow in a variably saturated porous medium in a stormwater constructed wetland, to the calibration of van Genuchten-Mualem parameters and to the evaluation of its efficiency in the reproduction of piezometric data collected on the Ostwaldergraben site. The temporal variability of the hydrodynamic parameters, including the hysteresis effect, shows that the characteristics of the filter layer alters along time. Both studies are carried using automatic differentiation.