De par le monde, la gestion des réservoirs hydroélectriques repose largement sur des modèles hydrologiques qui simulent les entrées et les sorties d’eau du système. L’évaporation sur ces réservoirs est typiquement négligée ou calculée de la même façon que l’évapotranspiration sur l’ensemble du bassin versant, bien que les processus régissant ces deux phénomènes diffèrent. Une meilleure quantification de l’évaporation suscite un intérêt grandissant, dû aux changements climatiques qui provoquent une augmentation des pertes en eau par évaporation sur les réservoirs. Cette étude examine divers modèles d’évaporation, et ce, en tenant compte de la réalité dans laquelle les modèles hydrologiques sont opérés de façon générale, c’est-à-dire avec une faible disponibilité de données d’entrée. L’étude se concentre principalement sur deux réservoirs situés dans la zone boréale canadienne, soit Eastmain-1 et Romaine-2. Un jeu de données observées d’évaporation obtenues grâce à la technique de covariance des tourbillons, est utilisé comme référence. De toutes les approches considérées, le modèle de transfert de masse avec un nombre de Dalton constant de 1,2 x 10⁻³ est le plus performant à un pas de temps horaire. Les cumuls journaliers sont également estimés avec une bonne précision autant lors du réchauffement de la masse d’eau que lors de la période de relâchement d’énergie. Des analyses de sensibilité montrent que l’utilisation de ce modèle avec des données limitées est possible, du moins pour les deux réservoirs étudiés. Ces données limitées consistent en la température de l’air et l’humidité relative prises sur la terre ferme à proximité du réservoir, la température de la surface de l’eau mesurée in situ et la vitesse du vent à une hauteur de 10 m prise depuis la berge du réservoir face aux vents dominants. L’application du modèle sur le lac Léman, situé dans une région plus tempérée, donne aussi de bons résultats. / Hydroelectric reservoirs around the world are managed using hydrological models to calculate water inflows and outflows. Evaporation on these reservoirs is typically neglected or calculated in the same way as evapotranspiration over the entire watershed, although the processes governing the two phenomena are different. A better quantification of evaporation is of growing interest in Quebec and internationally due to climate change, which is causing an increase in water losses by evaporation on reservoirs. This study therefore examines various evaporation models for open water bodies, taking into account the practical context in which hydrological models are generally operated, that is with low availability of input data. The study focuses mainly on two reservoirs with contrasting morphometries and located in the Canadian boreal zone, Eastmain-1 and Romaine-2. For this purpose, a set of observed evaporation data obtained using a direct measurement technique called eddy covariance, is used as a reference. Of all the approaches considered, the mass transfer model with a constant Dalton number of 1.2 x 10⁻³ gives the most accurate estimation of evaporation at hourly time steps. Daily totals are also estimated with good accuracy both during the warming of the water body in spring and during the energy release period in fall. Sensitivity analysis show that the use of this model with limited data is possible, at least for the two reservoirs studied. These limited data consist of air temperature and relative humidity taken on land in the vicinity of the reservoir, water surface temperature measured in situ, and wind speed at a height of 10 m taken from the bank of the reservoir facing the prevailing winds. The application of the model on Lake Leman, located in a more temperate region, also gives good results.

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.