Dans des conditions de climat humide et lorsque les sols sont peu perméables, les systèmes de drainage souterrains sont généralement utilisés pour contrôler le niveau de la nappe phréatique et améliorer la production agricole. Cependant, les drains souterrains modifient à la fois les voies d'écoulement hydrologique et les taux de transport des nutriments des terres cultivées vers les eaux de surface, pouvant détériorer la qualité des eaux souterraines et de surface. De plus, des macropores sont souvent présents dans les sols composés de till argileux, ce qui génère un flux d'eau rapide et riche en nutriments de la surface du sol vers les drains souterrains. Une approche rentable pour réduire le lessivage des nutriments provenant de l'agriculture consiste à imposer des restrictions uniquement dans les zones vulnérables à la contamination de l'eau. Ces zones peuvent être identifiées à l'aide de modèles hydrologiques distribués. Les résultats obtenus sur de petits bassins versants expérimentaux doivent être simplifiés pour être appliqués à des échelles plus grandes, généralement requises pour l'élaboration de politiques. L'objectif de cette étude était d'examiner les avancées et les limitations de l'inclusion des drains souterrains dans les modèles d'écoulement de surface et souterrain. Les objectifs spécifiques étaient de i) démontrer l'utilisation des estimations de conductivité électrique spécifique (CE), pour améliorer les simulations hydrologiques dans un champ drainé, ii) étudier l'efficacité d'un modèle hydrologique et de transport de soluté tridimensionnel pour simuler un test de traçage de bromure (Br) dans un champ drainé et iii) évaluer différents modèles conceptuels de drains souterrains et d'hétérogénéité du sol pour la simulation numérique du drainage dans un bassin versant agricole au Danemark. Les résultats suggèrent que la simulation de la profondeur de la nappe phréatique peut être améliorée par l'inclusion d'hétérogénéités basées sur des estimations de la CE. L'approche des seepage nodes était appropriée pour simuler les débits de drainage, cependant la précision des simulations était meilleure pour les modèles à l'échelle du terrain. À l'échelle du bassin versant, le fait de ne représenter que les drains principaux est approprié pour pouvoir utiliser des maillages plus grossiers et pour simuler le débit des cours d'eau et les faibles profondeurs des eaux de surface dans les zones drainées. Des résultats similaires ont été obtenus lorsque les seepage nodes ont été appliqués sur l'ensemble des zones agricoles, sans tenir compte de l'emplacement spécifique des drains souterrains. Cette dernière approche peut être appliquée lorsque les drains souterrains ne sont pas cartographiés, ce qui est généralement le cas. Une représentation simplifiée de l'hétérogénéité et de la macroporosité peut expliquer les différences entre ls valeurs observées et simulées des charges hydrauliques, débits de drainage et processus de transport de solutés. Les approches de modélisation étudiées dans cette thèse peuvent améliorer la représentation de la dynamique de l'écoulement souterrain et les simulations du transport de substances agrochimiques lessivées des champs cultivés, telles que le nitrate et phosphate. / Under humid climate conditions and for low-permeability soils, subsurface tile drains are usually employed to lower the water table and enhance agricultural production. However, tile drains alter both the hydrologic flow pathways in agricultural catchments and the rates of nutrient transport from cropland to surface water bodies, potentially impairing the groundwater and surface water quality. Furthermore, macropores are often present in clayey till soils, generating rapid and nutrient-rich water flow from the ground surface to the tile drains. A cost-effective approach to reduce nutrient leaching from agriculture is to impose restrictions only in vulnerable areas to water contamination, which can be identified using distributed hydrological models. Results on small experimental catchments need to be simplified for application on larger scales, usually required for policy-making purposes. The objective of this study was to investigate the outcomes and limitations of including tile drains in surface and subsurface flow models. Specific objectives were to i) demonstrate the use of electrical conductivity (EC) estimates to improve hydrological simulations in a tile-drained field, ii) investigate the efficiency of a three-dimensional hydrological and solute transport model to simulate a bromide (Br) tracer test in a tile-drained field and iii) assess different conceptual models for tile drains and soil heterogeneity for the numerical simulation of tile drainage in an agricultural catchment in Denmark. The results suggest that the simulation of the water table depth can be improved by the inclusion of heterogeneities based on EC estimates. The seepage nodes approach was suitable to simulate drain discharge, however the accuracy of the simulations was better for the field-scale models. At the catchment scale, representing only the main drains was suitable to reduce the mesh refinement and simulate stream flow and low surface water depths in drained areas. Similar results were obtained when seepage nodes were applied all over the agricultural areas, without considering the specific location of tile drains. The later approach can be applied when tile drains are not mapped, which is usually the case. The misrepresentation of heterogeneity and macroporosity may explain the differences between observed and simulated hydraulic heads, drain discharge and solute transport processes. The modeling approaches investigated in this dissertation can improve subsequent simulations of tile drainage and the transport and fate of leached agrochemicals such as nitrate or phosphate.

La pollution des eaux de surface par les pesticides est un problème sérieux. Cette pollution étant en grande partie liée à l’activité agricole, la mise en place de pratiques permettant de limiter les transferts de pesticides vers les eaux de surface représente une solution intéressante. Afin d’apprécier l’efficacité de tels projets il est important de développer des méthodes spécifiques d’évaluation agri-environnementale. Si les méthodes d’évaluation parcellaires sont opérationnelles, celles à l’échelle du bassin versant ne sont pas encore bien adaptées aux besoins des acteurs de terrain. Ce travail a pour objectif de développer un indicateur de risque d’exposition des eaux de surface aux pertes de pesticides par ruissellement à l’échelle du bassin versant. Pour ce faire, nous avons développé deux indicateurs : un indicateur parcellaire évaluant l’impact des pratiques agricoles sur les pertes de pesticides par ruissellement en sortie de parcelle et un indicateur de connectivité estimant la contribution de chaque parcelle au ruissellement global du bassin versant. Le couplage de ces deux indicateurs permet d’obtenir un outil d’évaluation du risque des pratiques phytosanitaires spécifique d’un bassin versant et de comparer des scenarii. Les méthodes de développement des deux indicateurs se distinguent des méthodes classiques par leur fondement sur un couplage indicateur / modèle. L’hypothèse méthodologique principale sur laquelle repose ce travail est qu’à partir des processus décrits par les modèles, il est possible d’explorer une gamme de situations très large et de créer une base de connaissances qui permet de développer nos indicateurs / Surface water pollution by pesticides is a great environmental concern. In order to limit pesticides transfers from agricultural field to surface water, mitigation strategies are implemented. Their effectiveness has then to be assessed by specific tools. Most of pesticide runoff risk assessment tools are simple indicators or simulation models. Both can be performing for field application but most of the available methods are not yet adapted to help catchment managers in their decision process. Simple indicators lack taking into account inherent variability of watersheds structures and mechanistic models demands too many data to be used as routine tools. The aim of this study is then to develop a catchment scale pesticide runoff risk assessment tool to help catchment managers to identify risky situations. We developed two indicators: one to assess agricultural impact on pesticide runoff risk at a field scale and the other field spatial position impact on runoff route from the field to the hydrological network. Aggregating those two indicators lead us to assess catchment fields pesticide runoff risk and allows comparing agronomical or climatic scenarios. We based both indicators on modelling assuming that processes described by model make possible to explore wide factors variation ranges and to build databases on which indicators are developed

La présence de résidus de médicaments a été évaluée dans des eaux superficielles et des effluents de huit stations d'épuration (quatre rejetant dans l'estuaire de Seine dans le cadre du projet Medseine financé par le GIP Seine- Aval, quatre dans le cadre du projet Toxstep financé par l'Anses (ex-Afsset). Une centaine de molécules a été systématiquement recherchée dans les deux types d'eaux en appliquant des méthodologies analytiques impliquant une extraction solide liquide suivie une analyse par LC-MS/MS. Un suivi saisonnier des rejets de médicaments a été effectué également. Des résidus de médicaments ont été systématiquement quantifiés dans les deux types d'eaux. Les concentrations mesurées dans des eaux de surface sont significatives pour certains composés : de 1 à 108 ng/L pour diclofénac, de 2 à 324 ng/L pour ibuprofène. Pour d’autres, les valeurs sont beaucoup plus faibles : de 2 à 13 ng/L pour sulfapyridine, de 1 à 5 ng/Lpour nordiazépam. Dans les effluents de STEP, les anti-inflammatoires, les β-bloquants, les macrolides et les fluoroquinolones sont plus souvent détectés. Leurs concentrations peuvent varier fortement : de 62 ng/L à 10 μg/L pour sotalol, de 39 ng/L à 8 μg/L pour aténolol et plus particulièrement de 32ng/L jusqu’à 127 μg/L pour paracétamol. Par ailleurs, la stabilité des résidus de médicaments vis à vis des rayonnements lumineux a été évaluée. Des irradiations UV (à 254 nm) pouvant potentiellement être utilisés en étape finaleavant rejet dans des zones particulièrement sensibles ont été appliquées (avec des doses allantjusqu'à 1500 J/m2. Par ailleurs, compte tenu du fait que les effluents sont ensuite majoritairement rejetés dans des eaux superficielles, ils sont soumis aux rayonnements solaires et cette exposition a également été évaluée. Le phénomène de photodégradation à 254 nm a été bien observé sur certains composés comme diclofénac, kétoprofène (> 90% dégradé), ciprofloxacine,norfloxacine (jusqu’à 75% dégradé). L’azithromycine, la clarithromycine, la carbamazépine et l’aténolol... sont des composés plus stables avec les rayonnements lumineux. En exposant sous lumière solaire simulée, le diclofénac, la kétoprofène, la ciprofloxacine et la norfloxacine sontaussi des composés sensibles (> 90% dégradé) après 24 h d’exposition (172,8 kJ/m2). La carbamazépine, la sulfaméthoxazole, l’aténolol, bisoprolol et métoprolol sont des composés plus stables observés. / The presence of pharmaceuticals was evaluated in Seine surface waters and in eight wastewater treatment plants (WWTP) effluents (4 WWTP within MEDSEINE - GIP Seine – Aval project and 4 WWTP within TOXSTEP – ANSES project. A hundred molecules was systematically reseached in both types of waters by applying the analytical methodologies (SPE solid phase extraction followed by LC-MS/MS analysis). A seasonal monitoring of pharmaceutical residues was also carried out. The pharmaceutically residues have been systematically quantified in this both types of water. The concentrations in suface of waters are significative for some molecules: 1 – 108 ng/L for diclofenac, 2 – 324 ng/L for ibuprofen and the values are weakly for some molecules: 2 – 13 for sulfapyridine, 1 – 5 for nordiazepam. Non steroidal anti-inflammatory drugs, β-blockers, macrolides and fluoroquinolones are mostly founded in the effluents of WWTP. Their concentrations were measured in the range: 62 ng/L – 10 µg/L for sotalol, 39 ng/L – 8 µg/L for atenolol, 32 ng/L – up to 127 µg/L for paracetamol.Moreover, the stability of pharmaceuticals was evaluated with light radiation . UV irradiation (at 254 nm) potentially used in the final WWTP before discharge in particularly areas have been applied (with a dose up to 1500 J/m2). Sunlight simulations have also been done. The phenomenon of photodegradation at 254 nm has been observed for some molecules such as : diclofenac, ketoprofen (>90% degraded), ciprofloxacine, norfloxacine (up to 75% degraded). Azithromycinem clarithromycine, carbamazepine and atenolol are most stable molecules observed. Exposing under simulated sunlight, diclofenac, ketoprofen, ciprofloxacine and norfloxacine are also sensitive molecules (>90 % degraded) after 24 hours of exposure (172.8 kJ/m2). Carbamazepine, sulfamethoxazole, atenolol, metoprolol and bisoprolol are most stable molecules observed.

Les eaux superficielles sont nécessaires à toute forme de vie en tant que parties intégrantes de tout processus de vie sur Terre. Quantifier les eaux de surface et suivre leurs variations est primordial en raison du lien direct qui existe entre les variables hydrologiques et le changement climatique. La télédétection par satellite, de l’hydrologie continental offre l’opportunité unique d’étudier, depuis l’espace, les processus hydrologiques à différentes échelles (régionale et globale). Dans cette thèse, différentes techniques ont été développées afin d’étudier les variations des eaux superficielles ainsi que d’autres variables hydrologiques, au niveau du bassin inférieur du Mékong (entre le Vietnam et le Cambodge) et ce en utilisant plusieurs estimations satellitaires différentes. Cette thèse s’articule autour de quatre points principaux. Premièrement, l’utilisation d’observations satellitaires dans le visible et dans l’infrarouge (MODIS) est étudiée et comparée afin d’évaluer les eaux de surface au niveau du bassin inférieur du Mékong. Quatre méthodes de classification ont été utilisées afin de différencier les types de surface (inondés ou pas) dans le bassin. Les différentes méthodes ont donné des cartes d’eaux de surface aux résultats semblables en terme de dynamique saisonnière. La classification la plus adaptée aux régions tropicales a été ensuite choisie pour produire une carte des eaux de surface à la résolution de 500 m entre janvier 2001 et aujourd’hui. La comparaison des séries temporelles issues de cette carte et de celles issues du produit de référence MODIS donne une forte corrélation temporelle (> 95%) pour la période 2001-2007. Deuxièmement, l’utilisation des observations issues du satellite SAR Sentinel-1 est examinée à des fins identiques. L’imagerie satellitaire optique est ici remplacée i par les images SAR qui grâce aux longueurs d’ondes utilisées dans le micro-ondes, permettent de « voir » à travers les nuages. Un jeu d’images Landsat-8-sans-nuage est alors utilisé pour entraîner un Réseau de Neurones (RN) afin de restituer des cartes d’eaux de surface par l’utilisation d’un seuillage sur les sorties du modèle RN. Les cartes sont à la résolution spatiale de 30 m et disponibles depuis janvier 2015. Comparées aux cartes de référence Landsat-8-sans-nuage, les sorties de modèles RN montre une très grande corrélation (90%) ainsi qu’une détection "vraie" à 90%. Les cartes restituées d’eaux de surface utilisant la technologie SAR sont enfin comparées aux cartes d’inondation issues de données topographiques. Les résultats montrent une fois encore une très grande consistance entres les deux cartes avec 98% des pixels considérés comme inondés dans cartes SAR se trouvant dans les régions de très grande probabilité d’inondation selon la topographie (>60%). Troisièmement, la variation volumique des eaux de surface est calculée comme le produit de l’étendue de la surface avec la hauteur d’eau. Ces deux variables sont validées à l’aide d’autres produits hydrologiques et montrent de bons résultats. La hauteur d’eau superficielle est linéairement interpolée aux régions non inondées afin de produire des cartes mensuelles à la résolution spatiale de 500 m. La hauteur d’eau est ensuite analysée pour estimer les variations volumiques. Ces résultats montrent une très bonne corrélation avec la variation volumique induite par la mesure du contenu en eau du satellite GRACE (95%) ainsi qu’avec la variation des mesures in situ de débit des rivières. Finalement, deux produits globaux et multi-satellites d’eaux superficielles sont comparés à l’échelle régionale et globale sur la période 1993-2007: GIEMS et SWAMPS. Lorsqu’elles existent, les données auxiliaires sont utilisées afin de renforcer l’analyse. Les deux produits montrent une dynamique similaire, mais 50% des pixels inondés dans SWAMPS se trouvent le long des côtes. / Surface water is essential for all forms of life since it is involved in almost all processes of life on Earth. Quantifying and monitoring surface water and its variations are important because of the strong connections between surface water, other hydrological components (groundwater and soil moisture, for example), and the changing climate system. Satellite remote sensing of land surface hydrology has shown great potential in studying hydrology from space at regional and global scales. In this thesis, different techniques using several types of satellite estimates have been made to study the variation of surface water, as well as other hydrological components in the lower Mekong basin (located in Vietnam and Cambodia) over the last two decades. This thesis focuses on four aspects. First, the use of visible/infrared MODIS/Terra satellite observations to monitor surface water in the lower Mekong basin is investigated. Four different classification methods are applied, and their results of surface water maps show similar seasonality and dynamics. The most suitable classification method, that is specially designed for tropical regions, is chosen to produce regular surface water maps of the region at 500 m spatial resolution, from January 2001 to present time. Compared to reference data, the MODIS-derived surface water time series show the same amplitude, and very high temporal correlation for the 2001-2007 period (> 95%). Second, the use of SAR Sentinel-1 satellite observations for the same objective is studied. Optical satellite data are replaced by SAR satellite data to benefit the ability of their microwave wavelengths to pass through clouds. Free-cloud Landsat-8 satellite imagery are set as targets to train and optimize a Neural Network (NN). Predicted surface water maps (30 m spatial resolution) are built for the studied region from January 2015 to present time, by applying a threshold (0.85) to the output of the NN. Compared to reference free-cloud Landsat-8 surface water maps, results derived from the NN show high spatial correlation (_90%), as well as true positive detection of water pixels (_90%). Predicted SAR surface water maps are also compared to floodability maps derived from topography data, and results show high consistency between the two independent maps with 98% of SAR-derived water pixels located in areas with a high probability of inundation (>60%). Third, the surface water volume variation is calculated as the product of the surface water extent and the surface water height. The two components are validated with other hydrological products, and results show good consistencies. The surface water height are linearly interpolated over inundated areas to build monthly maps at 500 m spatial resolution, then are used to calculate changes in the surface water volume. Results show high correlations when compared to variation of the total land surface water volume derived from GRACE data (95%), and variation of the in situ discharge estimates (96%). Fourth, two monthly global multi-satellite surface water products (GIEMS & SWAMPS) are compared together over the 1993-2007 period at regional and global scales. Ancillary data are used to support the analyses when available. Similar temporal dynamics of global surface water are observed when compared GIEMS and SWAMPS, but _50% of the SWAMPS inundated surfaces are located along the coast line. Over the Amazon and Orinoco basins, GIEMS and SWAMPS have very high water surface time series correlations (95% and 99%, respectively), but SWAMPS maximum water extent is just a half of what observed from GIEMS and SAR estimates. SWAMPS fails to capture surface water dynamics over the Niger basin since its surface water seasonality is out of phase with both GIEMS- and MODIS-derived water extent estimates, as well as with in situ river discharge data.

Thèse doctorat : Génie civil : Lille 1 : 2003. Thèse doctorat : Génie civil : Université Ibnou Zohr d'Agadir : 2003. / Thèse en cotutelle. N° d'ordre (Lille 1) : 3434. Résumé en français et en anglais. Bibliogr. p. 128-139.

Les estrogènes et les bêtabloquants transitent quotidiennement par les réseaux d'assainissement et les stations d'épuration des eaux résiduaires urbaines. Une première partie de ce travail a consisté à développer des méthodes pour analyser ces micropolluants émergents dans les boues de station d'épuration et les matières en suspension. L'efficacité des traitements épuratoires a ensuite été évaluée vis-à-vis de l'élimination de ces molécules. En règle générale, les 5 estrogènes étudiés (estrone, 17α- et 17β-estradiol, estriol, 17α-éthinylestradiol) sont plutôt bien éliminés par les stations d'épuration avec des efficacités d'élimination généralement supérieures à 90%. Pour les 10bêtabloquants analysés (acébutolol, aténolol, bêtaxolol, bisoprolol, métoprolol, nadolol, oxprénolol,propranolol, sotalol et timolol), des comportements différents ont été observés. Alors que l'acébutololet le nadolol sont éliminés (abattements moyens supérieurs à 80%), d'autres molécules comme lesotalol et le propranolol ne sont que peu impactées par les traitements épuratoires (abattements moyens inférieurs à 20%). Pour les autres bêtabloquants étudiés, des abattements intermédiaires (entre 40 et 70%) sont observés. Excepté le propranolol, qui est le plus hydrophobe des bêtabloquants étudiés, les molécules ciblées dans cette étude s'adsorbent peu sur les matières en suspension (en moyenne 90% des molécules se trouvent dans la phase dissoute) et ne sont pas concentrées dans les boues ; les abattements mesurés correspondent donc, sauf pour le propranolol,à une biodégradation et pas à un transfert des micropolluants vers les boues. Les molécules résiduelles non éliminées par le traitement épuratoire rejoignent le milieu aquatique récepteur. Les rejets de station d'épuration représentent en effet une source de contamination pour les milieux aquatiques et l'impact de ces rejets est visible, notamment en termes d'augmentation des concentrations de micropolluants mesurées dans le milieu. Néanmoins, les faibles concentrations mesurées ne semblent pas individuellement représenter de risque environnemental majeur ; seul le propranolol a présenté sur 4 des 15 sites étudiés un quotient de risque supérieur à 1 ce qui indique un risque environnemental probable. Cependant, même si le risque associé à une molécule est faible,chaque molécule présente dans l'environnement contribue au potentiel toxique global des substances présentes dans l'environnement.

Pour une meilleure gestion de l’eau, la Directive Cadre sur l’Eau requiert la prise en compte des relations entre les masses d’eau superficielles et souterraines. Dans ce cadre, la dynamique des échanges entre eaux de surface et eaux souterraines est étudiée, pour un aquifère alluvial du Rhône. Deux approches sont utilisées: une analyse des séries temporelles de niveaux de nappes et de rivières, à l’aide de techniques de traitement du signal, et des modèles d’écoulements numériques, déterministes, à base physique. Ces techniques sont mises en œuvre sur un secteur (Péage-de-Roussillon) à forts enjeux socio-économiques pour l’usage de la ressource en eau. Les résultats sont analysés du point de vue de leurs complémentarités.Une analyse en composantes principales, à partir des signaux piézométriques, a montré que les fluctuations de niveaux du Rhône expliquent la majeure partie des variations de niveau de la nappe. Les analyses corrélatoires et spectrales, ont permis de caractériser la relation existant entre les niveaux du Rhône et de la nappe. Des comportements particuliers de l’hydrosystème ont été identifiés : colmatage du fond des cours d’eau, écoulements transverses aux cours d’eau,…. Ces comportements ont ensuite pu être étudiés, plus en détail, à l’aide d’un modèle hydrodynamique de la nappe qui intègre un module de calcul des écoulements surfaciques. Le modèle permet également de quantifier les flux échangés entre la nappe et le cours d’eau. / For better water management, the Water Framework Directive requires to take into account the relationships between surface water and groundwater bodies.In this frame, the exchanges dynamic between surface water and groundwater is studied, for a Rhône alluvial aquifer. Two sets of tools are employed: a time series analysis of groundwater and rivers levels, using signal processing techniques, and numerical flow models, deterministic and physically based. These techniques are implemented on an area (Péage-de-Roussillon) with high socio-economic stakes regarding water resources. The complementarities among the results are analysed.A principal component analysis, based on piezometric head signals, showed that the fluctuations of the Rhône water level explain most of the groundwater variations. Correlative and spectral analysis were used to characterise the relationship between the Rhône and the groundwater level. Specific behaviours of the hydrosystem were identified: clogging of river beds, transversal flows below river beds,… These behaviours were then studied, in details, using a hydrodynamic model of the aquifer, which incorporates a surface runoff calculation module. The model also allows quantifying the exchange rates between rivers and groundwater.

Les systèmes hydrogéologiques devraient réagir au changement climatique de manière complexe. En région froide, la simulation de l'effet du changement climatique nécessite un modèle hydrologique intégré de pointe. Dans cette recherche, un modèle numérique entièrement couplé en 3D a été développé pour simuler l’écoulement des eaux souterraines et le transport de chaleur dans un bassin versant dans la région d'Umiujaq, dans le nord du Québec, au Canada. Le bassin versant est situé dans une zone de pergélisol discontinue et contient une épaisse couche glaciofluviale à grains grossiers formant un bon aquifère sous une unité gélive de silts marins sensible au gel. Une étude de terrain détaillée a été réalisée pour mesurer les caractéristiques du bassin versant telles que les propriétés hydrauliques et thermiques et la distribution des unités géologiques. Trois méthodes différentes disponibles dans le logiciel PEST sont utilisées pour caler le modèle 3D par rapport aux charges hydrauliques mesurées. Les résultats ont montré que l'utilisation de méthodes de calage simplifiées, telles que la méthode de zonation, n'est pas efficace dans cette zone d'étude, qui est très hétérogène. L’utilisation d’un calage plus détaillé par les méthodes du système PEST de points pilotes a permis de mieux s’adapter aux valeurs observées. Cependant, le temps de calcul était élevé. L'effet de la condition initiale pour la simulation du transport de chaleur est étudié en appliquant une condition initiale différente au modèle. Les résultats montrent que l'inclusion du processus de démarrage dans les simulations produit des températures simulées plus stables. Les zones du modèles à des profondeurs plus élevées, en-dessous de la profondeur de pénétration des variations saisonnières de température, nécessitent des temps de simulation plus longs pour être en équilibre avec les conditions limites appliquées. Les résultats montrent que l'application de la température moyenne de surface en tant que condition limite pour la simulation du transport de chaleur donne un meilleur ajustement aux valeurs observées en été qu'en hiver. En hiver, du fait de l’épaisseur variable de la neige dans le bassin versant, l’utilisation d’une température de surface uniforme diminue la qualité de l’ajustement aux valeurs observées. L'inclusion de l'advection dans la simulation du transport de chaleur accélère le taux d'augmentation de la température. De plus, l'eau chaude qui pénètre dans le sous-sol augmente la température souterraine dans les zones de recharge. Lorsque les eaux souterraines s'écoulent, elles perdent de l'énergie thermique. Par conséquent, le taux d’augmentation de la température dans les zones de décharge est inférieur à celui des zones de recharge. / Groundwater systems are expected to respond to climate change in a complex way. In cold regions, simulating the effect of climate change requires a state-of-the-art integrated hydrologic model. In this research, a fully coupled 3D numerical model has been developed to simulate groundwater-surface water flow and heat transport in a 2-km² catchment in Umiujaq, Nunavik (northern Quebec), Canada. The catchment is located in a discontinuous permafrost zone. It contains a lower aquifer, consisting of a thick coarse-grained glaciofluvial layer, overlain by a frost-susceptible silty marine unit and a perched upper aquifer. Detailed field investigations have been carried out to characterize the catchment, including its hydraulic and thermal properties and the subsurface geology. Three different calibration methods using the inverse calibration code PEST were used to calibrate the 3D flow model against measured hydraulic heads, assuming a fixed distribution of low hydraulic conductivity for discontinuous permafrost blocks. Heat transfer was not considered for this calibration. Results showed that using simplified calibration methods, such as the zoning method, is not efficient in this study area, which is highly heterogeneous. Using a more detailed calibration, such as the pilot-points method of PEST, gave a better fit to observed values. However, the computational time was significantly higher. In subsequent simulations, which included heat transport, different approaches for assigning initial temperatures during model spin-up were investigated. Results show that including the spin-up process in the simulations produces more realistic simulated temperatures. Furthermore, the spin-up improves the model fit to deeper subsurface temperatures because areas of the subsurface below the depth where seasonal surface temperature variations penetrate require longer simulation times to reach equilibrium with the applied boundary conditions. Applying the annual average surface temperature as the boundary condition to the heat transport simulation provided a better fit to observed values in the summer compared to winter. During winter, because of different snow thicknesses throughout the catchment, using a uniform surface temperature results in a poor fit to observed values. v Simulations show that warm water entering the subsurface increases the subsurface temperature in the recharge areas. As groundwater flows through the subsurface, it loses thermal energy. Therefore, discharging water is cooler than recharging water. This causes the rate of temperature rise to be lower in discharge areas than in recharge areas. The modelling results have helped provide insights into 3D simulation of coupled water flowheat transfer processes. Furthermore, it will help users of cryo-hydrogeological models in understanding effective parameters in development and calibration of model to develop their own site-specific models.

Ce projet de maîtrise a comme objectif de caractériser le ruissellement urbain pour un quartier résidentiel de la région de Québec. Un bassin d'orage a été sélectionné dans un quartier résidentiel nouvellement construit et une campagne d'échantillonnage de l'eau de ruissellement a été conduite à l'entrée et la sortie du bassin d'orage existant. L'ouvrage de sortie de ce bassin d'orage a ensuite été modifié pour y augmenter le temps de rétention et en évaluer l'effet sur l'enlèvement de polluants. Il a été trouvé que l'eau de ruissellement présentait des concentrations en polluants, tels les matières en suspension (MeS), métaux lourds et azote ammoniacal, semblables à celles trouvées dans la littérature. Les concentrations plus importantes sont relevées au printemps lors des premiers lessivages de rue après l'hiver. Le bassin d'orage original, ou bassin de rétention, a démontré une efficacité d'enlèvement moyenne de 56% pour les MeS, 20 % pour le zinc et 10 % pour l'azote ammoniacal. Son adaptation a permis de grandement augmenter ces efficacités d'enlèvement. Celles-ci ont atteint 82% pour les MeS, 42% pour le zinc et 84% pour l'azote ammoniacal. L'augmentation du temps de rétention du bassin d'orage permet donc d'améliorer considérablement l'enlèvement des polluants présents dans le ruissellement urbain avant qu'ils n'atteignent le milieu récepteur et permettrait donc d'améliorer la qualité de l'eau des cours d'eau en milieu urbain.

Titre de l'écran-titre (visionné le 27 avril 2023) / Une eau de qualité adéquate est essentielle au maintien de la vie aquatique et à un approvisionnement en eau potable salubre. La plupart des travaux passés portant sur la qualité de l'eau des rivières reposent sur des échantillonnages occasionnels, ce qui peut entraîner une incertitude considérable lors de l'estimation des charges en polluants à l'échelle annuelle. L'objectif de cette étude est d'identifier et quantifier l'impact des conditions hydrométéorologiques (e.g. pluie et fonte de la neige) sur les paramètres physico-chimiques de l'eau de rivière d'un bassin versant mixte agricole-forestier à partir d'observations à haute fréquence. Le site d'étude est le bassin versant de la rivière Beaurivage de 709 km² dans le sud du Québec (46,66°N; 71,29°O), en amont d'une importante prise d'eau municipale, et dont l'occupation du sol est dominée par la forêt (52%) et des champs agricoles (41%). Une campagne de terrain a eu lieu de décembre 2020 à mars 2022, avec des mesures en continu du débit, de la turbidité et de la conductivité électrique dans l'eau, et le suivi des conditions météorologiques. Au cours de 12 événements hydrométéorologiques associés à la dégradation de la qualité de l'eau (10 événements de pluie et 2 crues printanières), des mesures intensives de plusieurs paramètres de la qualité de l'eau ont été effectuées. Nos résultats montrent qu'il est possible d'estimer la concentration en azote total, en matière en suspension et en phosphore total à partir d'une relation empirique avec le débit. Pour chacune de ces relations, un point de rupture à 18,4 m³ s⁻¹ suggère qu'il y a deux phases ou plus dans la dynamique de transport des contaminants. Il est aussi possible d'estimer la concentration en carbone organique total à partir de la conductivité, et le nombre le plus probable de Escherichia coli dans l'eau à partir de la turbidité. Les résultats indiquent aussi que la charge événementielle de l'azote total, du carbone organique total, de la matière en suspension et du phosphore total peut être estimée par une relation empirique avec le 3ᵉ quartile de la distribution des débits enregistrés pendant un événement. Enfin, le bilan de la charge totale annuelle pour l'année 2021 est de 1 [0,5, 1,4] tonne N km⁻² pour l'azote total, de 14 [9, 18] tonnes C km⁻² pour le carbone organique total, de 34 [20, 70] tonnes km⁻² pour la matière en suspension et de 90 [55, 171] kg P km⁻² pour le phosphore total. / An adequate water quality is essential to maintain aquatic life and to supply safe drinking water. Most studies on river water quality are based on occasional sampling, which can bring considerable uncertainty when estimating pollutants concentrations or annual loads. The objective of our research is to identify and quantify the impact of hydrometeorological conditions (e.g. rain and snowmelt) on physico-chemical parameters of a river in a mixt agricultural-forest watershed based on high frequency observations. The study site is the Beaurivage river watershed (709 km² ) in the south of the province of Québec (46,66°N; 71,29°W). This watershed is mostly occupied by forest (52%) and agricultural fields (41%), and it is located upstream of a drinking water facility. A field campaign took place between December 2020 and March 2022 to obtain continuous measurement of the river flow, turbidity and electric conductivity, and the monitoring of meteorological conditions. During 12 hydrometeorological events associated with water quality degradation (10 rain and 2 snowmelt events), measurements of multiple water quality parameters have been carried out. Our results show that it is possible to estimate total nitrogen, suspended sediments and total phosphorus concentration with an empirical relation with the flow. For each of these relations, a fracture point at 18,4 m³ s⁻¹ suggested that there are two or more phases in the dynamic transport of contaminants. Also, it is possible to estimate the total organic carbon concentration with a relation using the conductivity, and the most probable number of E. coli in the river water using the turbidity. The results also indicate that the event load of total nitrogen, total organic carbon, suspended sediments and total phosphorus can be estimated with a empirical relation using the 3ᵗʰ quartile of the event flow distribution. Finally, the total annual load for the year 2021 is 1 tonne N/km² for the total nitrogen, 14 tonnes C/km² for the total organic carbon, 34 tonnes/km² for the suspended sediments and 90 kg P/km² for the total phosphorus.