L'étude des écoulements souterrains et du transport de contaminant en conditions de terrain est souvent très difficile, notamment en raison de l'hétérogénéité souvent importante des formations aquifères. Pour mieux comprendre les processus en jeu, les paramètres hydrodynamiques et hydrodispersifs sont souvent estimés en laboratoire dans des conditions entièrement contrôlées. L'objectif de ce mémoire est de vérifier s'il est possible de reproduire l'écoulement et le transport de contaminant dans une maquette hydrogéologique munie d'un simulateur de pluie en laboratoire dont les sables auraient été caractérisés en laboratoire. Cinq sables sont d'abord étudiés par des techniques classiques: des courbes granulométriques ainsi que des essais de Darcy et de traçage au NaCI dans un dispositif de Darcy de 43 cm de longueur. Ces techniques ont permis de déterminer les conductivités hydrauliques et les dispersivités des sables. Des essais d'écoulement et des traçages au NaCI sont ensuite réalisés en laboratoire dans une maquette hydrogéologique représentant un aquifère de dimensions 2 m x 1.32 m x 0.25 m. La maquette est d'abord remplie d'un seul des cinq sables pour représenter un aquifère homogène. Des expériences d'écoulement et de transport sont réalisés dans cet aquifère, avec et sans la recharge imposée au moyen d'un simulateur de pluie. La maquette hydrogéologique est par la suite vidée et remplie à nouveau d'une combinaison des cinq sables disposés en blocs de 4 cm x 4 cm x 25 cm suivant un variogramme exponentiel. Les expériences d'écoulement et de transport sont réalisées avec cet aquifère hétérogène, à nouveau avec et sans recharge. Les charges, les débits et les concentrations mesurés dans toutes ces expériences sont ensuite simulés avec un modèle Modflow et un modèle MT3D. Les paramètres calés pour reproduire les résultats de laboratoire sont comparés à ceux déterminés expérimentalement. Les résultats montrent que les conductivités hydrauliques estimées au moyen des courbes granulométriques sont peu représentatives des valeurs obtenues par l'expérience de Darcy. Les conductivités hydrauliques calées dans Modflow avec et sans recharge, pour les aquifères homogène et hétérogène sont quant à elles très similaires aux valeurs obtenues avec le dispositif de Darcy. Ceci indique que les paramètres hydrodynamiques mesurés en laboratoires sont fiables et représentatives des sables testés. Par ailleurs, les dispersivités mesurées dans le dispositif de Darcy sont jusqu'à 60% supérieures à celles calées dans MT3D pour reproduire les traçages réalisés en laboratoire. Ceci s'explique par la différence d'échelle entre le dispositif de Darcy (43 cm) et les blocs de l'aquifère hétérogène (4 cm). Ces résultats confirment qu'il est beaucoup plus difficile d'estimer en laboratoire les propriétés hydrodispersives d'un matériau aquifère que ses propriétés hydrodynamiques. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Traçage, Simulateur de pluie, Maquette hydrogéologique.

La péninsule gaspésienne est composée de petits bassins versants dont la dynamique fluviale est très sensible aux changements climatiques et environnementaux. Dans la perspective des changements climatiques en cours et anticipés, la variabilité des régimes d'écoulement, notamment ceux à charge sédimentaire liés aux cônes de déjection, devient moins prévisible et peut constituer un risque pour la sécurité publique, risque exacerbé par l'emplacement de ces cônes à l'interface versant-vallée. La présente étude constitue une analyse dendrogéomorphologique d'un petit cône de déjection en milieu forestier dans le but d'en caractériser la dynamique en termes de fréquence et d'intensité. Une démarche rétrospective faisant appel à la dimension historique des événements de forte intensité est privilégiée. La stratigraphie et la sédimentologie des dépôts étudiés indiquent que le développement du cône est associé à l'occurrence d'écoulements hyperconcentrés, processus hydrogéomorphologique encore peu connu au Québec. L'utilisation de méthodes issues de la dendrochronologie retrace la dynamique spatiotemporelle au sein des chenaux d'écoulement avec une résolution annuelle pour la période récente. Ces données fournissent un portrait précis des écoulements de forte intensité pour l'ensemble du XXe siècle et indiquent une récurrence approximative de neuf ans du processus à l'échelle du cône. Les déclencheurs météorologiques probables de ces événements extrêmes incluent des épisodes de pluie torrentielle ou une fonte nivale importante au printemps. Les résultats présentés apportent un éclairage nouveau sur la réponse des hydrosystèmes forestiers en contexte de changements climatiques. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Cône de déjection, dendrogéomorphologie, changements climatiques, écoulement hyperconcentré, risque.

La prévision des débits d'étiage est une question importante dans la gestion des bassins versants. Pendant les périodes de basses eaux, l'écoulement de base peut devenir une composante majeure des débits en rivière. Dans les modèles hydrologiques, la composante de l'écoulement souterrain est souvent très simplifiée. Ceci découle le plus souvent d'une caractérisation insuffisante des aquifères et les temps de calculs prohibitifs des modèles intégrant de manière complète l'écoulement souterrain. Cette thèse développe un modèle représentant les écoulements souterrains peu profonds et profonds qui nécessitent peu de paramètres et dont les calculs sont efficaces et fiables. Le modèle proposé peut être utilisé à la place d'un modèle d'écoulement souterrain en différences finies ou en éléments finis, mais il a été conçu dans le but d'être incorporé à un modèle hydrologique de bassin versant. Le modèle d'écoulement transitoire de versant "hillslope-storage Boussinesq" (hsB) (Troch et al., 2003, Paniconi et al., 2003) est sélectionné pour représenter l'écoulement souterrain peu profond à l'échelle locale. Le modèle hsB est couplé avec le modèle GLOW d'écoulement permanent basé sur la méthode des éléments analytiques (EA) (Strack, 1989; Haitjema, 1995) qui représente l'écoulement régional profond 2D horizontal. L'approche de couplage utilisée nécessite d'intégrer un terme de percolation à la base du modèle hsB. Ce terme de percolation est représenté par un flux vertical de Darcy à travers un aquitard hypothétique séparant le versant local de l'aquifère 2D régional. Afin de mieux comprendre les facteurs contrôlant la percolation, les facteurs tels que la géométrie du versant, l'inclinaison de la base, les propriétés des aquifères et les conditions limites sont évalués à l'aide d'un modèle 3D basé sur l'équation de Richards (Paniconi et Putti, 1994 ; Camporese et al., 2010). Les observations principales de cette analyse sont: i) l'eau peut circuler vers le bas ou vers le haut entre les aquifères de peu profonds et profonds et séparer les versants en trois zones distinctes: une zone de flux descendant à l'amont, une zone de flux descendant à l'aval et une zone de transition entre les deux; ii) l'inclinaison des versants et leur géométrie déterminent la partition des flux échangés. Ces résultats sont utilisés dans la mise en œuvre du couplage entre les modèles hsB et EA, où chaque versant et l'aquifère sous-jacent sont subdivisés en trois zones de percolation constante. Le modèle hsB/EA est testé sur différents types de versants et sur un bassin hypothétique formé de deux versants convergeant vers un cours d'eau. La comparaison avec les résultats d'un modèle numérique 3D basé sur l'équation de Richards et utilisé comme référence démontrent: i) que les charges hydrauliques, les taux de percolation et les débits aux exutoires sont généralement bien simulés; ii) de meilleurs résultats ont été obtenus pour les versants peu ou très peu inclinés aux géométries uniformes et divergentes; iii) les débits cumulés sont simulés de manière convenable pour le bassin hypothétique. Les écarts entre le modèle de référence sont attribués au fait que la zone non-saturée n'est pas représentée dans le modèle hsB et à l'hypothèse de Dupuit-Forchheimer du modèle EA qui néglige l'écoulement vertical dans la nappe profonde. Le fait que le modèle EA ne permette pas de simuler l'écoulement en régime transitoire est une autre limitation du modèle hsB/EA (e.g. Kuhlman et Neuman, 2009). Une application du modèle hsB/EA sur un bassin versant de 30 km2 situé dans la région de Covey Hill au sud du Québec a été réalisée. Dans cette application du modèle hsB/EA, les débits de base sont assez bien reproduits à l'exutoire du bassin versant pendant les périodes d'étiage. Toutefois, des écarts importants sont observés au cours des débits de pointe. Les erreurs sur les charges sont également non négligeables dans les zones où un gradient vertical a été observé. Ces erreurs peuvent être attribuées en parties aux limitations du modèle développé. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : hydrologie des versants, percolation, équation de Boussinesq, équation de Richards, éléments analytiques, eaux souterraines stratifié, modèle couplé

La rivière Noire est un sous-bassin de la rivière Yamaska. Elle est située dans la zone de transition entre deux formations géologiques principales du Québec, soit les Basses-Terres du Saint-Laurent et les montagnes appalachiennes. Dans ce secteur les dépôts meubles sont discontinus et d'étendue limitée; le socle rocheux est très fracturé et hétérogène. Cette région rurale fortement exploitée par l'agriculture intensive et l'industrie porcine présente un fort potentiel de contamination d'origine agricole de l'eau souterraine. Un échantillonnage d'eau souterraine effectué sur l'ensemble du bassin a permis de quantifier les paramètres physico-chimiques de l'eau, les principaux ions, incluant les teneurs en nitrates, ainsi que les contenus en ¹⁸O et deutérium. Ces données sont des indicateurs utiles pour identifier les types d'eau et pour établir des relations entrée-sortie sur l'aquifère. Ces informations, combinées aux données géologiques et à la carte piézométrique, permettent d'identifier des sous-secteurs dans les directions d'écoulement souterrain. Le modèle Modflow a été utilisé pour simuler les scénarios de géologie et de recharge du bassin et ainsi tester la pertinence des directions d'écoulement identifiées précédemment. Ainsi, la rivière Noire et certains de ses tributaires majeurs ont été identifiés comme étant des drains importants de l'aquifère. L'étude du bassin a montré que la direction régionale de l'écoulement de l'eau souterraine est orientée vers le nord-ouest et est principalement régie par la topographie et les propriétés du roc. L'aquifère est vulnérable puisqu'il est de nature libre. D'ailleurs, 31 % des puits échantillonnés présentent des concentrations en nitrates supérieures à 1,0 mg N-NO₃/l probablement d'origine anthropique. La complémentarité des approches utilisées dans le cadre de ce projet a été concluante en ce qui a trait à la caractérisation des écoulements souterrains sur un bassin régional. La modélisation représente adéquatement les caractéristiques du bassin mises en évidence par la caractérisation hydrogéologique et la géochimie. Une meilleure compréhension de la circulation des eaux souterraines sur le bassin versant de la rivière Noire aidera à délimiter les zones les plus vulnérables du bassin, ce qui permettra de réduire le potentiel de contamination de cet aquifère. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Hydrogéologie, Géochimie, Piézométrie, Modélisation, Modflow.

Selon le Groupe International d'Experts sur le Climat (lPCC, 2007), les changements climatiques vont entraîner une intensification du cycle hydrologique à l'échelle globale et un accroissement des précipitations dans certaines régions du monde, notamment celles situées aux latitudes moyennes et élevées. Des changements survenant au niveau de la quantité de précipitation saisonnière ainsi qu'au niveau de l'intensité et de la fréquence des événements extrêmes ont un impact direct sur l'amplitude des écoulements fluviaux saisonniers et sur la période d'occurrence et la fréquence des inondations et des sécheresses. De tels changements auront des impacts significatifs sur les ressources hydriques régionales. Cette étude se concentre sur la validation et l'évaluation des changements projetés au niveau des écoulements fluviaux moyens et au niveau de la période d'occurrence et de la fréquence des écoulements extrêmes, i.e. les écoulements de fort débit (crue) et de faible débit (étiage), pour les bassins canadiens sélectionnés. Cela se fait en utilisant un ensemble de simulations du Modèle régional du climat canadien correspondant au climat actuel (1961-1990) et à un climat futur (2041-2070) basé sur le scénario SRES A2. La validation est effectuée en évaluant les erreurs de performance et celles dues au pilotage, causées respectivement par la dynamique interne et la physique du modèle et par les erreurs associées au pilotage du modèle à ses frontières. Les résultats suggèrent des erreurs de performance positives des écoulements annuels moyens pour les bassins sans régulation situés dans la partie ouest du Canada (toujours supérieur à 30% sauf pour le bassin de l'Athabasca ou la différence n'est que de 4%) en raison d'une surestimation de l'équivalent en eau de la neige (SWE). Les erreurs dues au pilotage sont, en général, plus petites que les erreurs de performance (le coefficient d'habileté S est inférieur à 85% pour 12 des 14 bassins dans le cas des erreurs de performance alors que ce n'est le cas que de 2 bassins pour les erreurs de pilotage) et présentent sauf pour les bassins situés plus au sud, un biais négatif (pouvant aller jusqu'à -25%) La validation des étiages suggère que le modèle a quelques difficultés pour reproduire l'amplitude observée et la période d'occurrence des étiages, tandis qu'au niveau des crues, le modèle reproduit raisonnablement la période d'occurrence, quoique avec quelques différences entre les amplitudes observées et modélisées. En général, les résultats suggèrent une augmentation de l'amplitude de l'écoulement hivernal et un pic de fonte de neige survenant plus tôt (une à deux semaines) pour les bassins situés plus au nord, de même que des changements significatifs quant aux caractéristiques des crues et des étiages. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Changement climatique, crue, étiage, modèle climatique régional, écoulement fluvial.

Le concept de vulnérabilité de l'eau souterraine est généralement défini comme « la capacité du milieu à transmettre un polluant depuis la surface jusqu'à la zone saturée ». Bien qu'indispensable pour la protection et la gestion de la ressource en eau souterraine, ce concept montre des limites importantes puisqu'il n'intègre presque jamais la dynamique des écoulements souterrains et rarement l'agencement lithostratigraphique des dépôts quaternaires. Ce mémoire vise à déterminer la vulnérabilité de l'aquifère au roc de la zone Bécancour. L'étude est réalisée par l'application de l'indice DRASTIC et par le développement d'un indice de vulnérabilité basé sur une caractérisation hydrogéochimique de la zone d'étude. La comparaison des deux indices ainsi que la pertinence des paramètres utilisés est discutée et mise en perspective avec l'occurrence des nitrates dans les puits de la zone d'étude. Enfin, les concepts utilisés pour l'estimation de la vulnérabilité sont mis en perspective en s'appuyant notamment sur les résultats d'un modèle numérique synthétique d'écoulement et de transport en zone saturée. La recharge moyenne pour l'ensemble de la zone étudiée est estimée à 103 mm/an. La recharge est principalement contrôlée par la présence de till sur le territoire où elle atteint des valeurs voisines de 100 mm/an. Des taux de recharge inférieurs à 40 mm/an sont estimés pour la partie aval du bassin, en raison de la présence de dépôts fins de la Mer de Champlain. Des maximums de 250 mm/an sont évalués dans le piémont appalachien aux endroits où les dépôts quaternaires sont minces. Les résultats de l'étude hydrogéochimique confirment la présence de zones de recharge privilégiées dans le piémont appalachien avec la présence d'eaux souterraines peu minéralisées de type Ca-HCO3. Des conditions de confinement croissantes vers le Fleuve Saint-Laurent sont également mises en évidence avec des eaux de plus en plus minéralisées et de plus en plus basiques évoluant successivement vers des eaux de type Na-HCO3 puis Na-Cl. Bien qu'aucun dépassement de la norme de potabilité pour les nitrates (10 mg N-NO3/L) n'ait été détecté, les résultats montrent la présence de nitrates dans 30% des forages échantillonnés (concentrations supérieures à 1 mg N-NO3/L, attribuables à des sources anthropiques). Parmi les forages affectés, 37% sont des puits tubulaires au roc, la plupart situés en zone amont du bassin, dans le piémont appalachien. L'indice de vulnérabilité DRASTIC (Aller et al., 1987) est calculé sur la zone d'étude, et montre une cohérence limitée avec le contexte géomorphologique et avec l'occurrence des nitrates sur le bassin. L'indice hydrogéochimique montre une vulnérabilité plus importante en amont qu'en aval du bassin, mais permet de discerner des zones de vulnérabilité localement élevées en aval du bassin, sans doute exposées à des fenêtres locales de recharge, notamment pour la zone sud-sud-ouest. L'utilisation d'un modèle numérique d'écoulement-transport synthétique en zone saturée a permis de confirmer que le taux de recharge est un paramètre majeur de la vulnérabilité et que son effet est couplé avec le potentiel de dilution à l'échelle régionale. Les zones simultanément associées à de forts taux de recharge et à de faibles taux de dilution représenteraient des zones plus vulnérables à la contamination depuis la surface. Ce constat est plausible avec les résultats du bassin de la zone Bécancour où l'on observe l'occurrence de nitrates dans le piémont appalachien plutôt qu'en aval du bassin, bien que la densité des activités anthropiques et agricoles soit plus denses proche du Fleuve Saint-Laurent. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : vulnérabilité eau souterraine, recharge, hydrogéochimie, nitrates, modélisation numérique, DRASTIC, zone Bécancour

À l'échelle de la planète, les milieux humides représentent entre 4 et 6 % de la superficie des terres émergées. Au Canada, cette proportion atteint près de 15 %. Ces écosystèmes, dont l'hydrologie est étroitement liée à l'hydrogéologie de surface, sont fragiles et essentiels au maintien de la biodiversité. Le présent projet de maîtrise est réalisé dans le cadre d'un projet de caractérisation des eaux souterraines financé par le MDDEP et a pour objectif principal de comprendre le rôle des tourbières dans la dynamique de l'aquifère de la partie basse de la zone Bécancour. Une cartographie de l'ensemble des tourbières du territoire a été réalisée et leur degré de perturbation a été déterminé par photo-interprétation. La même démarche a été entreprise pour les photos aériennes de 1966, afin de connaître l'évolution des superficies dans le temps. Les résultats montrent qu'il y a 10 485 ha de tourbières intactes et 7 015 ha de tourbières légèrement perturbées. Les perturbations les plus fréquentes sont les canaux de drainage et les chemins forestiers. De plus, 5 433 ha de tourbière ont subi des perturbations irréversibles entre 1966 et 2010 sur la partie basse de la zone Bécancour, principalement en lien avec l'exploitation de la canneberge. Le contexte géomorphologique de neuf tourbières a été étudié au moyen de sondages manuels, de levés géophysiques par géoradar et d'analyses granulornétriques. Les résultats montrent que les tourbières de la zone d'étude reposent toutes sur un substrat au moins partiellement perméable. Les substrats des tourbières à l'étude sont composés à 44 % sur un substrat de sable éoliens et de till remanié, à 33 % sur un dépôt de till remanié, à 11 % sur un dépôt deltaïque et de till remanié et à 11 % sur un dépôt marin et de till remanié. Les neuf tourbières étudiées ont une profondeur maximale variant entre 1,63 et 6,40 m et une profondeur moyenne entre 1,10 et 2,94 m. Les flux entrants et sortants des neuf sites ont été quantifiés à l'aide de l'approche segmentée de l'équation de Darcy. Les flux échangés entre les tourbières et l'aquifère se situent entre -1097 et 409 m3/j, cinq sites alimentant l'aquifère, trois sites recevant l'eau de l'aquifère et un n'ayant pas d'échange entre les milieux. Représentant en moyenne entre 4 et 11 % des apports totaux de la tourbière, il a été possible d'extrapoler à l'ensemble des tourbières, et d'estimer les flux transitant entre les deux milieux à l'échelle de la zone d'étude. Le flux transitant de l'aquifère vers la tourbière est estimé entre 0,0033 et 0,0099 rn3/s et le flux de la tourbière vers l'aquifère est estimé entre 0,0056 et 0,0168 m3/s. Ces proportions suggèrent que la tourbière joue un rôle tampon relativement limité. Localement, ces flux permettent de garder les niveaux piézométriques relativement stables en redistribuant l'eau au système du surface ou souterrain. ______________________________________________________________________________

On estime qu'environ 50% des milieux humides mondiaux ont été détruits et convertis en terre agricole et en milieu urbain. Ces milieux humides sont d'importants écosystèmes qui jouent non seulement un rôle fondamental dans la préservation de la biodiversité, mais aussi dans le maintien de la qualité et de la quantité d'eau potable contenue dans les rivières et aquifères environnants. Ce projet de recherche avait pour objectif de quantifier les échanges entre le complexe tourbeux de Lanoraie, les rivières présentes sur le territoire et l'aquifère des sables du delta de Sorel. Une compilation des données existantes et des mesures de terrain (GPS différentiel, relevés de coupes stratigraphiques, niveaux de nappe, débits) ont permis de caractériser la tourbière et l'aquifère régional. À l'aide de ces données, un modèle hydrostratigraphique 3D a été construit. Quatre types d'écoulement ont été identifiés : 1) un écoulement suivant la topographie de l'aquifère vers la zone minérotrophe de la tourbière, 2) un écoulement suivant la topographie de la zone ombrotrophe à la zone minérotrophe, 3) un écoulement alimentant les rivières présentes sur le territoire (Saint-Joseph, Point-du-Jour, Saint-Jean, Saint-Antoine et Bras du sud-ouest) et 4) un écoulement de la tourbière à l'aquifère sableux de l'ancien delta de Sorel. L'ensemble de ces travaux a été utilisé pour construire un modèle numérique en régime permanent à l'aide du logiciel MIKE SHE. Focalisant sur les échanges aquifère-tourbière-rivière et l'impact potentiel des changements climatiques extrêmes sur les niveaux de nappes et les débits de base des différentes rivières, le modèle a permis de montrer qu'il existe des liens hydriques importants entre l'aquifère du delta de Sorel, le complexe des tourbières de Lanoraie et les rivières environnantes. Par exemple, les résultats du modèle numérique montrent que 1) 52% de l'eau qui circule annuellement dans la tourbière vient directement de l'aquifère sableux de l'ancien delta de Sorel, 2) les débits de base des différentes rivières environnantes sont constitués majoritairement (48% à 100%) d'eau ayant préalablement circulée dans le complexe tourbeux de Lanoraie et 3) des diminutions de 20% et 50% de la recharge annuelle causeraient des diminutions de 16% et 41% respectivement, en moyenne, des débits de base de la région, ce qui pourrait avoir des conséquences importantes sur les réserves en eau et l'écologie des rivières. Ce projet de recherche montre que la préservation des réserves hydriques de la région de Lanoraie passe non seulement par la protection de la ressource en eau potable et des rivières, mais aussi celle du complexe tourbeux. Ceci est encore plus vrai dans un contexte de changement climatique et de pressions anthropiques. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Tourbières, aquifère, rivière, Mike SHE, changement climatique, Lanoraie

La région de l'Abitibi renferme d'abondants dépôts fluvioglaciaires de sables et de graviers, notamment des eskers, reconnus comme des aquifères exceptionnels étant donné leurs caractéristiques hydrogéologiques intrinsèques et la faible perméabilité du roc sous-jacent. Une meilleure compréhension des régimes hydrauliques relatifs à ces systèmes est nécessaire afin d'établir une gestion adaptée de la ressource hydrique. Ainsi, une étude ayant recours aux isotopes des gaz rares fut initiée dans les eskers de la région d'Amos. Douze échantillons d'eaux furent prélevés dans l'esker Saint-Mathieu-Berry (SMB), la moraine d'Harricana, l'esker de Barraute et dans la plaine argileuse située entre ces formations. Spécifiquement, des analyses sur les concentrations en gaz rares ont été effectuées. Les résultats des rapports 20Ne/4He contre 3He/4He indiquent un mélange entre une composante radiogénique enrichie en 4He, tritiogénique enrichie en 3He et atmosphérique. Le rapport mesuré 3He/4He normalisé à celui de l'atmosphère (R/Ra), où Ra=1,386*10-6, varie entre 0.224±0.012 pour l'eau d'un puits privé prélevée dans le roc (P1) et 1.849±0.036 pour le puits d'Amos. Les âges 3H-3He, obtenus grâce aux contenus en tritium et en 3He mesurés, dévoilent la présence d'eaux modernes (≤32 ans) situées principalement dans les formations fluvioglaciaires. La moraine d'Harricana, marquée par une signature atmosphérique, présente les eaux les plus jeunes variant entre 6.44±1.30 et 9.80±1.29 ans. Les âges 3H-3He montrent une relation linéaire avec la profondeur des puits pour l'esker Saint-Mathieu-Berry (SMB) et la moraine d'Harricana. Cette relation, appuyant l'utilisation du modèle d'écoulement piston pour la plupart des puits, permet de calculer des gradients verticaux d'infiltration variant entre 0.9 et 7 m/an. Elle permet d'identifier un phénomène de dispersion/diffusion dans les puits du secteur sud de l'esker de SMB appuyé par la comparaison entre les concentrations en tritium et la courbe des précipitations à Ottawa. De ce fait, la datation 3H-3He ne semble pas valide pour les puits du secteur sud. La calibration du taux de production en 4He permet l'estimation des temps de séjour (U-Th)/4He des eaux dans les aquifères. Ces derniers montrent l'occurrence de masses d'eaux intermédiaires (50 à 6000 ans) et fossiles (≥6000 ans). Quatre échantillons de l'esker de SMB montrent des eaux possiblement d'âge intermédiaire (50 à 738 ans) selon les analyses (U-Th)/4He, mais modernes selon les analyses 3H/3He. Ces puits, en grande partie situés dans une zone de faille de type indéterminé, montre une signature d'un mélange entre une eau ancienne circulant dans les zones fragilisées du socle granitiques sous-jacent l'esker et de l'eau jeune. Les eaux de l'esker ne présentent pas de signature typique des eaux profondes circulant dans la plaine argileuse. D'autre part, l'âge apparent (U-Th)/4He le plus élevé de la zone d'étude est d'environ 8000 ans. Cet âge correspond à celui d'un puits situé dans le socle fracturé de la plaine argileuse. La salinité et le rapport 3He/4He mis en fonction de la profondeur des puits suggèrent un mélange entre de l'eau minéralisée riche en 4He radiogénique et de l'eau douce moderne plus superficielle dans l'aquifère du socle rocheux de la plaine. Cela peut également indiquer une augmentation graduelle de la salinité engendrée lors du transport de l'eau puisque le temps de séjour augmente aussi avec la distance longitudinale par rapport à la moraine d'Harricana. Ainsi, la moraine semble représenter une zone de recharge importante à la fois pour l'esker de Barraute et pour le socle de la plaine argileuse. Notons que les eaux les plus minéralisées et riches en 4He radiogénique sont celles provenant du socle fracturé Archéen de la plaine argileuse. Les températures de recharge obtenues varient entre -0.10±0.61°C et 11.64±0.73°C. Pour la plupart des puits, elles diminuent avec l'augmentation de l'élévation estimée pour la recharge. Les échantillons situés dans la plaine montrent les températures de recharge les plus basses (entre -0.10 et 2.91°C), alors que celles de l'esker SMB se situent entre 6.19 et 11.6°C. La paléotempérature pour la moraine d'Harricana est plus faible en surface (3.87°C) qu'en profondeur (9.52°C). ______________________________________________________________________________

Depuis les dernières décennies, l'hydroélectricité est devenue une source majeure d'électricité pour le Québec, tout en représentant une source importante d'exportation d'énergie pour cette province canadienne. Malgré le fait qu'il est de plus en plus évident que toutes les formes d'électricité ont des répercussions environnementales, les impacts cumulatifs des grands barrages hydroélectriques du nord ne sont pas bien étudiés. Il y a encore de grandes questions à aborder, particulièrement concernant les effets cumulatifs et à long terme sur les écosystèmes estuariens. Les rivières Romaine et Sainte-Marguerite se situent sur la côte-nord du Québec, et tandis que la rivière Sainte-Marguerite est déjà harnachée, des barrages sur la rivière Romaine sont présentement en construction. Les deux rivières coulent directement dans le golfe du Saint-Laurent, un écosystème marin ayant une diversité animale et végétale importante. En utilisant deux approches différentes, les études de laboratoires et les images satellites, l'objectif de cette étude était de déterminer les impacts des barrages hydroélectriques sur les écosystèmes marins. Au laboratoire, notre but était de déterminer comment le taux de croissance et la photosynthèse de deux espèces de phytoplancton d'eau douce retrouvées dans l'estuaire de la rivière Romaine sont affectés par les changements de salinité et de température. Avec l'aide d'images satellites provenant du Landsat-TM5 bande 6 et les images SeaWiFS, il a été possible d'approfondir nos connaissances quant aux changements entraînés sur l'estuaire par le barrage de la rivière Sainte-Marguerite, avant et après la mise en eau en 1998. Le taux de croissance et la photosynthèse d'Ankistrodesmus falcatus et Fragilaria crotonensis ont été affectés par une augmentation de salinité. Le taux de croissance a diminué linéairement avec l'augmentation de la salinité et une diminution significative du rendement photochimique (Φ’M) a été observé en fonction de ce dernier. La température de l'eau a aussi affecté le taux de croissance ainsi que la photosynthèse, et une interaction a été notée entre la salinité et la température. Ces résultats ont été corroborés avec les images satellites. Dans le cas de la rivière Sainte-Marguerite, nous avons déterminé que la superficie du panache diminue linéairement en fonction des débits de la rivière. En plus, il y avait une différence significative entre la biomasse de phytoplancton dans l'estuaire avant et après la construction du barrage. Les résultats suggèrent que les répercussions environnementales peuvent être importantes parce que les débits d'eau douce qui arrivent dans le golf du Saint-Laurent ont déjà diminué à cause des barrages hydroélectriques. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : photosynthèse, salinité, température, fluorescence, estuaires, A. falcatus, F. crotonensis, rivière Sainte-Marguerite, rivière Romaine, hydroélectricité, Landsat-TM5, SeaWiFS, production primaire