L'écotone forêt-toundra (EFT) est un environnement où la forêt boréale amorce une transition vers la toundra arctique. Longue de 13 400 km, cette interface forme probablement la plus grande zone de transition sur Terre et couvre de grandes parties du nord de l'Eurasie et de l'Amérique du Nord. Sa taille en fait d'ailleurs un facteur d'importance grandissante avec le réchauffement global du climat sur Terre. En effet, les propriétés fort différentes de la forêt et de la toundra se font sentir sur le bilan énergétique de surface, qui décrit les échanges d'énergie et de masse à l'interface sol-atmosphère. En somme, ce bilan reflète la répartition du rayonnement net en flux de chaleur sensible et latente, ainsi qu'en flux de chaleur dans le sol. Jusqu'à présent, ce bilan a été peu étudié dans l'EFT, bien que son rôle crucial de couplage entre l'atmosphère et le sol soit essentiel pour le régime thermique et hydrologique de la surface terrestre. L'objectif de cette étude est d'analyser le bilan d'énergie de surface d'un site du Bas-Arctique dans l'EFT, sur la côte est de la baie d'Hudson à l'est du Canada, en été et en hiver. Pour ce faire, nous utilisons les données recueillies par une tour micrométéorologique utilisant l'approche de la covariance des tourbillons. Le site choisi est la vallée Tasiapik, d'une longueur de 4.5 km, où les parties hautes de la vallée sont couvertes par une toundra arbustive qui se transforme en forêt boréale vers les parties basses de la vallée. De plus, nous comparons les observations recueillies avec des simulations produites à l'aide des modèles de surface ISBA et SVS pour l'été et du modèle de neige Crocus en hiver. La comparaison avec ces modèles de surface est particulièrement importante comme ils sont utilisés avec des modèles atmosphériques pour générer des prévisions météorologiques et des projections climatiques. En été, nous avons constaté que 23% du rayonnement net était converti en flux de chaleur latente sur notre site, 35% en flux de chaleur sensible et environ 15% en flux de chaleur du sol. Ces résultats contrastent avec ceux de six sites FLUXNET répartis dans l'Arctique, où la majeure partie du rayonnement net sert à alimenter le flux de chaleur latente, et ce, même s'ils ont tous une précipitation annuelle bien inférieure à celle du site étudié. Nous attribuons ce comportement à la conductivité hydraulique élevée du sol (présence de sédiments littoraux et intertidaux), typique de ce que l'on trouve dans les régions côtières de l'est de l'Arctique canadien. Les modèles de surface terrestre ISBA et SVS surestiment la teneur en eau de surface de ces sols, mais parviennent à simuler avec précision les flux de chaleur turbulents, notamment le flux de chaleur sensible et, dans une moindre mesure, le flux de chaleur latente. En hiver, la couverture neigeuse modifie complètement le bilan énergétique de surface. Les pertes par rayonnement infrarouge sont en grande partie compensées par le flux de chaleur sensible, tandis que le flux de chaleur latente est minime. À la surface du couvert, le flux de chaleur dans la neige est similaire en amplitude au flux de chaleur sensible. Comme la couverture neigeuse stocke très peu de chaleur, l'amplitude du flux de chaleur dans la neige est comparable au flux de chaleur dans le sol. Dans l'ensemble, Crocus est en mesure de reproduire le bilan énergétique observé, mais démontre quelques défaillances lors de la simulation des flux de chaleur turbulents à un pas de temps horaire en conditions atmosphériques stables. Comme les deux types de couvert végétal dans l'EFT, la toundra et la forêt, ont un effet contrasté sur la couverture neigeuse au sol, nous avons analysé les propriétés du manteau neigeux à chacun des deux environnements. Nous avons à nouveau eu recours au modèle de neige Crocus pour simuler les propriétés de la neige sur les deux sites. D'abord, nos observations montrent que la hauteur et la densité de la neige diffèrent considérablement d'un site à l'autre. Sur le site forestier, la neige est environ deux fois plus épaisse qu'au site de toundra, et la densité diminue du sol vers la surface de la neige, alors que le contraire est observé sur le site de toundra. Crocus n'est pas en mesure de reproduire ces profils de densité dans sa configuration standard. En ajustant le paramétrage de la densité de la neige fraîche, de la compaction en présence de végétation et la neige soufflée, nous arrivons à simuler des profils comparables aux observations. Nous émettons l'hypothèse que le transport de la vapeur d'eau est le mécanisme dominant qui façonne le profil de densité sur le site de la toundra, alors qu'une combinaison de la compaction due au poids des couches superficielles et du transport de la vapeur d'eau est responsable du profil de densité sur le site forestier. Ce processus n'est pas inclus dans Crocus et les ajustements compensent partiellement ce fait et sont, dans une certaine mesure, spécifiques au site, ce qui complique l'application à plus grande échelle des modifications mises en œuvre ici. / The forest-tundra ecotone (FTE) is an environment where the boreal forest begins a transition to the arctic tundra. With an extent of about 13,400 km, this interface is probably the largest transition zone on Earth and covers large parts of northern Eurasia and North America. Its size makes it a factor of increasing importance for the Earth's climate with global warming. Indeed, the very different properties of forest and tundra are reflected in the surface energy balance, which describes the energy and mass exchanges at the soil-atmosphere interface. In short, this balance reflects the distribution of net radiation into sensible and latent heat fluxes, as well as heat fluxes into the soil. To date, this balance has been little studied in the FTE, despite its crucial role in coupling the atmosphere and the soil which is essential for the thermal and hydrological regime of the land surface. The objective of this study is to analyze the surface energy balance at a Low-Arctic site in the FTE, on the east coast of Hudson Bay in eastern Canada, in summer and winter. To do so, we use data collected by a micrometeorological tower using the eddy covariance approach. The study site is the Tasiapik valley, 4.5 km long, where the upper parts of the valley are covered by shrub tundra that transforms into a boreal forest towards the lower parts of the valley. In addition, we compare the collected observations with simulations produced using the ISBA and SVS surface models in summer and the Crocus snow model in winter. The comparison with surface models is particularly important as they are used with atmospheric models to generate weather forecasts and climate projections. In summer, we found that 23% of the net radiation was converted to latent heat flux at our site, 35% to sensible heat flux, and about 15% to ground heat flux. This contrasts with six FLUXNET sites across the Arctic, where most of the net radiation is used to drive the latent heat flux, even though they all have much lower annual precipitation than our study site. We attribute this behavior to the high hydraulic conductivity of the soil (presence of littoral and intertidal sediments), typical of the coastal regions of the eastern Canadian Arctic. The ISBA and SVS land surface models overestimate the surface water content of these soils, but are able to accurately simulate turbulent heat fluxes, including sensible heat flux and, to a lesser extent, latent heat flux. In winter, the snow cover completely changes the surface energy balance. Energy losses due to longwave radiation are largely offset by the sensible heat flux, while the latent heat flux is minimal. At the surface of the snow cover, the heat flux in the snow is similar in magnitude to the sensible heat flux. As the snow cover stores very little heat, the magnitude of the heat flux in the snow is comparable to the heat flux in the ground. Overall, Crocus is able to reproduce the observed energy balance, but shows some shortcomings when simulating turbulent heat fluxes at an hourly time step under stable atmospheric conditions. Since the two vegetation types in the FTE, tundra and forest, have a contrasting effect on the snow cover on the ground, we analyzed the snowpack properties at each of the two environments. Again, we used the Crocus snow model to simulate snow properties at both sites. First, our observations show that snow height and density differ significantly between the two sites. At the forest site, the snow is about twice as thick as at the tundra site, and the density decreases from the ground to the snow surface, while the opposite is observed at the tundra site. Crocus is not able to reproduce these density profiles in its standard configuration. By adjusting the density settings for fresh snow, compaction in the presence of vegetation and blown snow, we are able to simulate profiles comparable to the observations. We hypothesize that water vapor transport is the dominant mechanism shaping the density profile at the tundra site, while a combination of compaction due to the overburden weight and water vapor transport is responsible for the density profile at the forest site. This process is not included in Crocus and the adjustments partially compensate for this and are to some extent site-specific, making it difficult to apply the modifications implemented here on a larger scale.

Aucun résumé disponible ______________________________________________________________________________

Le schéma de surface canadien (version 3.5) a été testé sous un sol agricole avec des observations provenant d’une station météorologique située dans un champ de pommes de terre. Les observations couvrent une saison complète en 2007, et la période de tubérisation du plant pour l’année 2010. Les performances du modèle sur les composantes des bilans énergétique et hydrique ont pu être évaluées. Sur les deux saisons à l’étude, le modèle a surestimé les flux latents et produit un flux de chaleur dans le sol de trop grande amplitude, comparativement à celui observé. L’évapotranspiration et le drainage produits par CLASS ont par la suite été comparés à ceux simulés par le modèle d’écoulement vertical BV3C, un sous-modèle du modèle hydrologique HYDROTEL. BV3C a semblé produire une évapotranspiration surestimée dont la dynamique ne suit pas celle observée. Finalement, une simulation dont les intrants météo ont été remplacés par le produit régional du modèle global multi-échelle, pour la même localisation, a été testée. Les composantes des bilans énergétique et hydrique simulés par CLASS et ISBA, un autre schéma de surface utilisé au Canada, ont pu être comparés. Lors de ce test, ISBA a montré un plus grand ruissellement que CLASS, ainsi que de faibles quantités d’eau drainée, en désaccord avec les observations disponibles. / The performance of the Canadian Land Surface Scheme (CLASS) was assessed using meteorological observations from stations located in a potato field during a whole season in 2007, and from tuber initiation to tuber bulking in 2010. The scheme’s performance has been evaluated for both the energetic and hydric balance. Over the two seasons, CLASS overestimated the latent flux and generated a ground heat flux of higher amplitude compared to the observed one. The evapotranspiration and the drainage generated by CLASS were then compared to those of BV3C, a submodel of the hydrological model HYDROTEL. BVC3 seemed to have produced an overestimated evapotranspiration, which the dynamics did not match with those of the observed one. Finally, both CLASS and ISBA, an other land surface scheme in use in Canada, were run with data derived from the regional product of the global environmental multiscale model, for the same location. During this test, ISBA showed a higher flow than CLASS. ISBA also showed little quantity of drained water in disagreement with the available observations.

Cette étude a pour but de comparer et d'évaluer les bilans hydrique et énergétique de surface des versions opérationnelle et en développement du Modèle Régional Canadien du Climat (MRCC) avec les ré-analyses atmosphériques et les observations de surface disponibles. Les versions opérationnelle et en développement du MRCC se distinguent par l'emploi du module de paramétrage des processus physiques du Modèle de Circulation Générale Canadien (MCGC) de deuxième et troisième génération respectivement, développé par le Centre Canadien de modélisation et de l'analyse Climatique (CCmaC). Parmi les améliorations apportées à la physique du MCGC III, soulignons l'implantation d'un nouveau schéma de surface appelé CLASS (Canadian LAnd Surface Scheme), un modèle à trois couches de sol avec traitement explicite des couverts végétal et nival; il remplace le modèle hydrologique Beautified Bucket et le régime thermique de force de rappel à une couche de sol utilisés dans la version opérationnelle du MRCC. Les simulations sont issues de la phase 1c du projet PIRCS couvrant l'ensemble des États-Unis entre 1987 et 1994. L'analyse est réalisée sur les bassins-versants des fleuves Mississippi et Columbia. Certains flux de surface et champs diagnostics d'intérêts sont comparés avec les ré-analyses du NCEP/NCAR sur 50 ans et du ECMWF sur 40 ans. Les taux de précipitation, la température à l'abri et l'écart diurne de température à l'abri sont également comparés avec des climatologies d'observations. Les estimés mensuels de ruissellement sur les deux bassins-versants proviennent de mesures de débits fluviaux et d'aires de drainage avec quelques modifications. Enfin, des analyses de la profondeur de neige et des flux radiatifs de surface sont également utilisées. Les résultats montrent que la version en développement du MRCC tend à simuler une meilleure climatologie que la version opérationnelle, particulièrement en ce qui concerne l'évapotranspiration, l'écart diurne de température à l'abri et la précipitation estivale. Cependant, un biais froid persistant de température à l'abri sur les deux bassins est associé à un couvert de neige excessif suivi d'une importante crête de ruissellement au printemps. De par sa construction et contrairement au schéma CLASS, le modèle hydrologique Beautified Bucket favorise la glace dans le sol au dépend de la couverture de neige. Des suggestions de modifications simples pour chacune des versions du MRCC sont également proposées. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Bilan hydrique de surface, Bilan énergétique de surface, MRCC opérationnel, MRCC en développement, Bassin-versant du Mississippi, Bassin-versant du Columbia.

Titre de l'écran-titre (visionné le 25 mars 2024) / Au cours des dernières décennies, la télédétection a grandement facilité la documentation des tendances associées au raccourcissement de la période d'englacement des lacs en réponse au changement climatique. Bien que ces investigations soient diversifiées, certaines lacunes persistent, notamment l'exclusion des petits lacs et des réservoirs des études globales, et ce malgré leur abondance dans l'hémisphère nord. Ainsi, l'objectif de ce projet est de caractériser l'englacement des lacs et réservoirs à travers l'hémisphère nord grâce à l'imagerie satellitaire Sentinel-2A. Il vise également à examiner les tendances spatiales des dates de phénologie et à évaluer l'influence des caractéristiques géographiques, morphométriques et climatiques associées aux plans d'eau. Ces travaux reposent sur le développement d'une vaste base de données regroupant les dates de phénologie de la glace de 3702 lacs et de 1028 réservoirs à travers l'hémisphère nord de 2019 à 2023. En outre, cette base de données a permis la calibration d'un modèle prédictif de l'englacement des plans d'eau. Ce modèle, de type forêt aléatoire, a démontré une bonne performance avec une erreur absolue moyenne de 11 et 13 jours pour la prise et la fonte de la glace, respectivement. Cette marge d'erreur s'explique principalement par les causes d'incertitudes associées à la détection de glace par imagerie optique, notamment la gestion de la couverture nuageuse, mais aussi la détection de glace noire en climat aride. Parmi les variables d'entrée du modèle, la latitude et l'accumulation maximum des degrés-jours de gel constituent les paramètres déterminants dans la prédiction des dates de gel et de fonte de la glace. Bien que l'utilisation des capteurs optiques présente des défis spécifiques, le potentiel de la télédétection reste significatif dans le cadre d'études globales sur la glace, offrant également la possibilité d'améliorer les modèles climatiques actuels. L'étude de cas du réservoir de la Romaine-2 illustre un exemple concret de la synergie entre la télédétection et la modélisation. / In recent decades, remote sensing has considerably facilitated the documentation of trends related to the shortening of lake ice duration in response to climate change. Although these studies are diverse, certain gaps remain, in particular the exclusion of small lakes and reservoirs from global studies, despite their abundance in the Northern Hemisphere. Using Sentinel-2 satellite imagery, this project aims to characterize lake and reservoir freeze-up period across the Northern Hemisphere, and to assess the influence of their geographic, morphometric and climatic characteristics. This work is based on the development of an extensive database of ice phenology dates for 3702 lakes and 1028 reservoirs across the Northern Hemisphere from 2019 to 2023. This database was also used to calibrate a predictive model of lake and reservoir ice conditions. This random forest model performed well, with mean absolute errors of 11 and 13 days for freeze-up and break-up, respectively. This error is mainly due to uncertainties associated with ice detection from optical imagery, particularly cloud cover management, but also black ice detection in arid climates. Among the model input variables, latitude and maximum accumulation of freezing degree-days were the determining parameters in predicting freeze-up and break-up dates. Although the use of optical sensors presents specific challenges, the potential of remote sensing in the context of global ice studies remains significant and offers the possibility of improving current climate models. The case study of the Romaine-2 reservoir provides a concrete example of the synergy between remote sensing and modeling.

La forêt boréale recouvre 30% de la surface du Canada et 14% de la surface émergée de la terre. Les changements climatiques vont fortement l’affecter, et ces écosystèmes vont significativement impacter la climatologie et l’hydrologie mondiale avec des échanges importants d’eau, d’énergie et de carbone entre le sol et l’atmosphère. Il est maintenant crucial de bien comprendre le bilan d’énergie de surface sur ce biome pour efficacement prévoir son comportement et son évolution dans un climat changeant. Beaucoup d’études ont analysé le bilan d’énergie de la forêt boréale, mais des lacunes importantes subsistent : il y a peu d’études en terrain non plat, ou dans des zones recevant d’importantes précipitations, ou encore avec des mesures à des échelles spatiales variées, et encore moins des combinaisons de ces trois possibilités. Le principal but de cette thèse est d’attaquer ces lacunes de front en offrant l’analyse du bilan d’énergie et de l’évapotranspiration d’une forêt boréale humide couvrant une topographie prononcée, et ce à plusieurs échelles spatiales (point : ~m²; locale : ~ha; régionale : ~km²). Les résultats sont principalement basés sur une campagne de mesures ayant cours à la forêt Montmorency de l’Université Laval, à 80 km au nord de Québec, Canada. La forêt est une sapinière à bouleau blanc avec des arbres à divers degrés de maturité. Là, deux tours à flux mesurant tous les termes du bilan d’énergie sont installés depuis l’automne 2015. Trois objectifs spécifiques sont associés à trois échelles spatiales de mesure ou de modélisation dans un gradient allant de l’échelle ponctuelle vers l’échelle régionale. Dans un premier objectif, l’hétérogénéité spatiale du couvert forestier est caractérisée par des mesures de rayonnement solaire sous-canopée. Ensuite, l’évaluation de la densité de végétation permet de paramétrer un schéma de surface pour obtenir la variabilité de l’évapotranspiration et de ses composantes. Les résultats montrent que même si la transmission du rayonnement est très variable de point en point (entre 7% et 69% sur toute la saison), une moyenne spatiale à l’échelle locale représente bien la zone. Les résultats de modélisation indiquent qu’une forêt plus dense cause légèrement plus d’évapotranspiration totale, car elle évapore plus de précipitations interceptées et transpire plus. Une forêt plus dense évapore toutefois moins d’eau au sol, ce qui peut mener à une humidité du sol accrue dans des conditions de sécheresse momentanée. / Dans le deuxième objectif, l’impact de fortes précipitations sur le bilan d’énergie à l’échelle locale et l’évapotranspiration en forêt boréale est évalué. Pour ce faire, le site principal de la forêt Montmorency est d’abord comparé avec 13 sites en forêt boréale dans le monde sur la base du bilan d’énergie et de l’évapotranspiration. La forêt Montmorency est le site qui reçoit le plus de précipitations avec ~1600 mm a⁻¹. Pour tous les sites, la précipitation reçue est liée à l’évapotranspiration annuelle, et c’est aussi le site principal qui a la plus forte évapotranspiration, avec ~550 mm a⁻¹. Grâce à des mesures précises des débits sortants du bassin versant de 3.5 km² contenant les sites de mesure, le bilan hydrique est clairement établi : l’eau excédentaire provenant des précipitations est principalement évacuée par des débits sortants dans les cours d’eau du bassin versant, à une hauteur de ~1050 mm a⁻¹. Pour le troisième objectif, la méthode de la scintillométrie à deux longueurs d’onde est testée sur le site d’étude et ses mesures du bilan d’énergie à l’échelle régionale sont comparées à celles à l’échelle locale. Les scintillomètres sont installés à travers une vallée où est également présente une des deux tours à flux. Les faisceaux électromagnétiques parcourent 1347 m à une hauteur variant entre 5 et 100 m et une hauteur efficace de 88 m. Les résultats montrent que les deux montages expérimentaux ont une faible concordance au niveau des paramètres de structure météorologiques, mais une concordance plus qu’acceptable au niveau des flux turbulents. Pour ces derniers, la corrélation entre les scintillomètres et la tour à flux est optimale lorsque les faisceaux électromagnétiques sont entièrement inclus dans la couche de surface atmosphérique. Cependant, comme la hauteur des faisceaux est hautement variable, ceux-ci se retrouvent la plupart du temps en partie dans la couche de surface atmosphérique, ce qui mène à une corrélation tout de même acceptable dans ces circonstances. Néanmoins, les mesures des scintillomètres sont souvent de qualité douteuse lors des périodes nocturnes et lorsque l’atmosphère est stable. En bref, la forêt boréale étudiée exhibe un bilan d’énergie et une évapotranspiration significativement différents d’autres sites dans des biomes semblables recensés dans la littérature. Cette thèse apporte des précisions importantes sur ce type d’environnement. De plus, la thèse offre des outils méthodologiques rigoureux pour évaluer le bilan d’énergie à diverses échelles spatiales et élabore sur le passage entre ces échelles, une contribution à ne pas négliger pour les modélisateurs hydrologiques et du climat au Canada et dans le monde. / For the third objective, the two-wavelength scintillometry method is evaluated at the study site and its regional energy balance measurements are compared to those at the local scale. The scintillometers are installed across a valley where one of the two flux towers is localized. The scintillometers’ electromagnetic beams travel 1347 m at a height varying between 5 and 100 m and an effective height of 88 m. The results show that the two experimental systems have a low agreement in terms of the meteorological structure parameters, but a more than acceptable agreement for the turbulent fluxes. For the latter, the correlation between scintillometers and flux tower is optimal when the electromagnetic beams are entirely included in the atmospheric surface layer. However, since the beam height is highly variable, they are more often than not partially present in the atmospheric surface layer anyway, which leads to a correlation that is still acceptable in these circumstances. However, measurements of scintillometers are often unrealistic during nocturnal periods and when the atmosphere is stable. In short, the studied boreal forest exhibits an energy balance and evapotranspiration significantly different from other sites in similar biomes referenced in the literature. This thesis provides important details on this type of environment. In addition, the thesis offers rigorous methodological tools to assess the energy balance at various spatial scales and elaborates on the possibility of upscaling and/or downscaling results, a contribution not to be overlooked for hydrological and climate modelers in Canada and around the world. / The boreal forest covers 30% of Canada's surface and 14% of the earth's land surface. Climate change will severely affect it, and these ecosystems will in turn impact climate and global hydrology with significant exchanges of water, energy and carbon between the soil and the atmosphere. It is now crucial to understand the surface energy balance of this biome to effectively predict its behavior and evolution in a changing climate. Many studies have analyzed the energy balance of the boreal forest, but significant gaps remain: there are little studies in non-flat terrain, or in areas receiving significant rainfall, or with measurements at various spatial scales, let alone combinations of these three possibilities. The main aim of this thesis is to fill these gaps with a rigorous analysis of the energy balance and evapotranspiration of a boreal forest covering a pronounced topography, and this at several spatial scales (point: ~m², local: ~ha, regional: ~km²). The results are mainly based on a measurement campaign taking place at the Montmorency Forest of Université Laval, 80 km north of Québec, Canada. The forest is a balsam fir – white birch forest with trees of varying degrees of maturity. There, two flux towers are measuring all the energy balance terms since autumn 2015. Three specific objectives are associated with three spatial scales of measurement or modeling in a gradient from the point scale to the regional scale. In a first objective, the spatial heterogeneity of the forest cover is characterized by sub-canopy solar radiation measurements. Then, the vegetation density evaluation makes it possible to parameterize a land-surface scheme to obtain the variability of the evapotranspiration and its components. The results show that even though the transmission of radiation is highly variable from point to point (seasonal average between 7% and 69%), a spatial average at the local scale represents the area quite well. Modeling results indicate that a denser forest causes slightly more total evapotranspiration because it evaporates more intercepted precipitation and generates more transpiration. A denser forest, however, evaporates less water on the ground, which can lead to increased soil moisture under conditions of momentary drought. In the second objective, the impact of heavy rainfall on the local energy balance and evapotranspiration in the boreal forest is evaluated. To do this, the main site of Montmorency Forest is first compared with 13 boreal forest sites around the world on the basis of energy balance and evapotranspiration. The Montmorency Forest is the site receiving the most rainfall with ~1600 mm y⁻¹. For all sites, the precipitation received is positively related to annual evapotranspiration, which means the main site has the highest evapotranspiration rates, with ~550 mm y⁻¹. With accurate measurements of the outflow from the 3.5 km² watershed containing the Montmorency Forest measurement sites, the water balance is clearly established: excess water from precipitation is mainly discharged through outflows of the watershed, to an extent of ~1050 mm y⁻¹.

Les réservoirs constituent des ouvrages de retenue d'eau qui permettent de supporter une vaste gamme d'activités humaines telles que la production d'électricité et d'eau potable, l'irrigation, la navigation, etc. L'Est de la région boréale canadienne compte une part importante de réservoirs hydroélectriques. Le contexte actuel des changements climatiques affecte fortement ces écosystèmes et modifie ainsi l'hydrologie et la climatologie régionale via les échanges d'eau, d'énergie et de gaz à effet de serre (dioxyde de carbone, méthane, vapeur d'eau) avec l'atmosphère. Parmi ces échanges, l'évaporation, une composante essentielle des modèles climatique et hydrologique, demeure encore difficile à estimer à ce jour. Le bilan hydrique d'un réservoir reflète l'équilibre entre les flux massiques entrant et sortant, et permet d'anticiper l'évolution des volumes d'eau disponibles pour supporter les différents usages anthropiques. Il se compose des flux entrant et sortant, de la précipitation, de l'évaporation, qui mis ensemble peuvent faire varier le volume d'eau stocké. Via l'évaporation, le bilan d'énergie d'un réservoir est couplé au bilan hydrique, ce qui en motive son étude. Le bilan d'énergie d'un réservoir compare les flux de chaleurs advectifs, turbulents (sensible et latent) et de rayonnement net. Même si beaucoup d'études ont analysé et quantifié le bilan énergétique des plans d'eau, des lacunes demeurent. En effet, peu d'entre elles ont été effectuées sur des réservoirs hydroélectriques, de surcroît en zone boréale. De plus, lorsque disponibles, les observations ont généralement des portées spatiale et temporelle limitées. L'objectif principal de cette thèse est de pallier ces insuffisances en réalisant une analyse des échanges eau-atmosphère d'un réservoir hydroélectrique profond à plusieurs échelles spatiales (locale ~ ha; régionale ~ km²) et temporelles (journalière, mensuelle et annuelle), puis en quantifiant les bilans de masse et d'énergie en incluant les échanges advectifs liés au turbinage du réservoir. Notre démarche expérimentale s'appuie sur une campagne de mesures réalisée sur le réservoir Romaine-2 (50.68°N, 63.25°O), exploité par Hydro-Québec depuis 2015, situé à 243 m d'altitude au Québec, Canada. Le réservoir présente des profondeurs moyenne de 44 m et maximale de 101 m, une superficie maximale de 85.6 km² avec un marnage annuel maximal de 17 m. Deux tours à flux mesurant le bilan d'énergie thermique ont été déployées de juin 2018 à juin 2022, l'une sur la berge et la seconde sur un quai flottant ancré au fond du réservoir et déployé en période d'eau libre chaque année. L'objectif principal de la thèse est décliné en trois objectifs spécifiques. Le premier objectif consiste à évaluer l'applicabilité sur un plan d'eau d'une méthode récente de mesure des flux turbulents à l'échelle régionale (~ km²), soit la scintillométrie à deux longueurs d'onde. Les résultats sont comparés avec ceux de la méthode de référence à l'échelle locale (~ ha), la covariance des tourbillons, réalisée à partir d'un quai flottant. La méthode de scintillométrie repose sur deux couples d'émetteurs/récepteurs installés de part et d'autre du réservoir Romaine-2 et émettant deux faisceaux situés pour l'un dans l'infrarouge et pour l'autre dans le domaine des micro-ondes, sur une distance de 1745 m et à une hauteur approximative de 10 m au-dessus de la surface du plan d'eau. Les résultats révèlent une concordance acceptable des flux de chaleur sensible, mais moins probante quant aux flux de chaleur latente qui sont surestimés par rapport à la méthode locale de covariance des tourbillons. L'empreinte de mesure plus large des scintillomètres peut expliquer ces différences en captant une plus grande hétérogénéité dans les flux. Enfin, la différence de température eau-air se révèle être un bon indicateur du régime de stabilité et par conséquent de la direction (i.e., signe) du flux de chaleur sensible, initialement mal attribué par la scintillométrie. Pour le deuxième objectif, les variabilités journalière, mensuelle et interannuelle de l'évaporation mesurée à l'échelle locale au-dessus du réservoir Romaine-2 sont quantifiées et analysées. Les résultats montrent un déphasage de 12 heures entre les flux de chaleurs sensible et latente pendant la période d'eau libre. Le flux de chaleur sensible répond avant tout à la différence de température eau-air qui est maximale la nuit et minimale l'après-midi, tandis que le flux de chaleur latente est corrélé à l'énergie apportée par le rayonnement solaire qui est maximale l'après-midi et minimale la nuit. Annuellement, l'évaporation atteint 590 mm en moyenne (minimum de 555 mm, et maximum de 656 mm), ce qui représente environ 51 % de la précipitation annuelle. 84% de l'eau est évaporée entre août et décembre, période pendant laquelle le réservoir largue intensément sa chaleur stockée dans une atmosphère plus froide. L'évaporation annuelle cumulée a dû être corrigée à la hausse par application de la fraction d'énergie perdue sur une année énergétique (aucun stockage net). Pour le troisième objectif, l'évolution temporelle du régime thermique du réservoir est caractérisée par des mesures de profil thermique de la colonne d'eau, à l'aide de deux chaînes de thermistors, et ce entre juin 2018 et juin 2022. Les données de turbinage offrent la possibilité de quantifier le bilan hydrique du réservoir ainsi que le bilan d'énergie complet associé. Les résultats montrent que les couches supérieures affichent des décalages d'amplitudes thermique et temporel avec les couches plus profondes. La variabilité interannuelle de la thermocline reste faible, mais sa profondeur et son gradient thermique varient en fonction du niveau d'eau et du turbinage. Le bilan hydrique est dominé en entrée et en sortie par les débits turbinés (61.4% du réservoir amont et 88.0% via la centrale en aval), tandis que le bilan d'énergie est principalement gouverné en entrée par le rayonnement net (62.3%) et en sortie dans les mêmes proportions par les flux de chaleurs sensible et latente (41.2%) et le débit sortant du réservoir (37.4%). Cette thèse offre donc une analyse méthodique et structurée de résultats obtenus sur une longue période de campagne en milieu éloigné, sur la base de méthodes de mesure peu usitées jusqu'alors dans une région climatique parfois hostile d'un réservoir hydroélectrique en milieu boréal côtier, la basse Côte Nord du Québec. / Reservoirs are water retention structures that support a wide range of human activities such as power generation, drinking water, irrigation, navigation, etc. The eastern part of the Canadian boreal region has a significant amount of hydroelectric reservoirs. The current context of climate change strongly affects these ecosystems and thus modifies regional hydrology and climatology through the exchange of water, energy and greenhouse gases (carbon dioxide, methane, water vapour) with the atmosphere. Among these exchanges, evaporation, an essential component of climate and hydrological models, remains difficult to estimate to this day. The water balance of a reservoir represents the balance between incoming and outgoing mass flows, and allows to anticipate the evolution of the volumes of water available to support the different anthropic uses. It is composed of lateral and upstream inflows and outlet flows (i.e. spillway and turbine), precipitation and evaporation, which together can vary the volume of water stored. Through evaporation, the energy balance of a reservoir is coupled to the water balance, which motivates its study. The energy balance of a reservoir compares advective, turbulent (sensible and latent) and net radiation heat fluxes. Although many studies have analysed and quantified the energy balance of water bodies, there are still gaps. Indeed, few of them have been carried out on hydroelectric reservoirs, especially in boreal zones. Moreover, when available, the observations generally have limited spatial and temporal scope. The main objective of this thesis is to overcome these shortcomings by analyzing the water-atmosphere exchanges of a deep hydroelectric reservoir at several spatial (local ~ ha; regional ~ km²) and temporal (daily, monthly and annual) scales, and then quantifying the mass and energy balances by including the advective exchanges related to the reservoir turbining. Our experimental approach is based on a measurement campaign carried out on the Romaine-2 reservoir (50.68°N, 63.25°W), operated by Hydro-Québec since 2015, located at 243 m altitude in Quebec, Canada. The reservoir has an average depth of 44 m and a maximum depth of 101 m, a maximum surface area of 85.6 km² with a maximum tidal range of 17 m. Two flux towers measuring the thermal energy balance were deployed from June 2018 to June 2022, one on the shore and the second on a floating raft anchored to the bottom of the reservoir and deployed during the open water period each year. The main objective of the thesis is broken down divided into three specific objectives. The first objective is to evaluate the applicability on a water body of a recent method for measuring turbulent flows on a regional scale (~ km²), i.e. two-wavelength scintillometry. The scintillometry method is based on two transmitter/receiver installed on either side of the Romaine-2 reservoir and emitting two beams, one in the infrared and the other in the microwave bands, over a distance of 1745 m and at a height of approximately 10 m above the water surface. The results show acceptable agreement for sensible heat fluxes, but less agreement for latent heat fluxes which are overestimated compared to the local eddy covariance method. The larger footprint of the scintillometers may explain these differences by capturing greater heterogeneity in the fluxes. Finally, the water-air temperature difference turns out to be a good indicator of the stability regime and therefore of the direction (i.e., sign) of the sensible heat flux, initially poorly assigned by scintillometry. For the second objective, the daily, monthly and interannual variabilities of evaporation measured at the local scale above the Romaine-2 reservoir are quantified and analysed. Results show a 12-hour phase shift between sensible and latent heat fluxes during the open water period. The sensible heat flux responds primarily to the water–air temperature difference which is maximum at night and minimum in the afternoon, while the latent heat flux is related to the energy provided by solar radiation which is maximum in the afternoon and minimum at night. Annually, evaporation averages 590 mm (minimum and maximum 555 mm and 656 mm respectively), which represents about 51% of the annual precipitation. 84% of the water is evaporated between August and December, when the reservoir intensively releases its stored heat into a colder atmosphere. The annual cumulative evaporation is then corrected upwards by considering the ratio of the energy budget over an energy year (no net storage). As a third objective, the temporal trends of the reservoir thermal regime is characterized by thermal profile measurements of the water column using two thermistors chains, between June 2018 and June 2022. The turbining data provides the opportunity to quantify the water budget of the reservoir and the associated full energy budget. Results show that there are temporal thermal amplitude lags between the surface and the deeper layers. The interannual variability of the thermocline remains small, but the depth and thermal gradient vary with water level and turbining. The water budget is dominated at the inlet and outlet by the turbined flows (61.4% from the upstream reservoir and 88.0% via the downstream power station), while the energy budget is mainly governed at the inlet by net radiation (62.3%) and at the outlet by both the sensible and latent heat fluxes (41.2%) and reservoir flow (37.4%). Ultimately, this thesis provides a methodical and structured analysis of results obtained over a long period of fieldwork in a remote environment. It is based on measurement technics that have not been widely used up to now in a hostile climatic region of a coastal boreal hydroelectric reservoir, the lower Côte Nord of Québec.