Ce travail s'inscrit dans le cadre de recherche sur les maladies vectorielles de Lyme et Virus du Nil au sein de l'Agence de Santé Publique du Canada (ASPC) ayant pour finalité d'évaluer et de cartographier les risques sanitaires associés à ces maladies infectieuses liées au climat aux échelles municipales, provinciales et fédérale. Dans ce contexte, cette recherche vise à démontrer la faisabilité, la pertinence et les enjeux de recourir aux méthodes de réduction d'échelle pour obtenir à une haute résolution spatio-temporelle (100/30 m et 1 jour) avec au plus des marges d'erreur de 2 unités, des déterminants climatiques et microclimatiques (DCMC) en milieu hétérogène du Canada. Un cadre méthodologique d'application des méthodes de réduction d'échelle, Random Forest Regression (RFR), Thermal sharpening (TsHARP), Pixel block intensity modulation (PBIM), a été proposé pour estimer la température de surface (LST) de MODIS 1000 m à 100/30 m. Des expérimentations basées sur cette approche ont été effectuées sur trois sites au Québec à différentes époques. Les résultats, spatialement représentatifs, ont été validés avec les températures de l'air et celles prises par de Landsat 08 avec des marges d'erreur autour de 2°C. L'analyse des résultats démontre la capacité effective des méthodes de réduction d'échelle à discriminer la LST dans l'espace. Toutefois, dans le contexte du projet de l'ASPC, ces résultats sont non concluants à 100/30 m en l'absence d'une plus-value significative au plan spatial de LST. Cette analyse a conduit à discuter des enjeux temporels, spatiaux, méthodologiques et de gestion de gros volumes de données en lien avec la réduction d'échelle dans le contexte du projet. / This research is part of the Public Health Agency of Canada's (PHAC) research on Lyme and West Nile Virus vector-borne diseases, which aims to assess and map the health risks associated with these climate-related infectious diseases at the municipal, provincial and federal levels. In this context, this research aims to demonstrate the feasibility, relevance and challenges of using downscaling methods to obtain high spatial and temporal resolution (100/30 m and 1 day), with margins of error of no more than 2 units, of climatic and microclimatic determinants (CMDs) in a heterogeneous Canadian environment. A methodological framework for the application of downscaling methods, Random Forest Regression (RFR), Thermal sharpening (TsHARP), Pixel block intensity modulation (PBIM), has been proposed to estimate the surface temperature (LST) from MODIS 1000 m to 100/30 m. Experiments with our approach were carried out at three sites in Quebec at different times. The spatially representative results were validated with air and Landsat 08 temperatures with error margins around 2°C. The analysis of our results demonstrates the effective capacity of downscaling methods to discriminate LST in space. However, in the context of the ASPC project, these results are inconclusive at 100/30 m in the absence of a significant, expected increase in the spatial accuracy of LST. This analysis led to a discussion of the temporal, spatial, methodological and large data volume management issues related to downscaling in the context of the project.

Titre de l'écran-titre (visionné le 29 mai 2023) / L'amplification polaire du changement climatique, qui conduit les hautes latitudes à se réchauffer près de deux fois plus vite que les régions tempérées, entraîne la diminution progressive de certaines zones englacées du globe comme celle du pergélisol en Arctique. En profondeur, le dégel du pergélisol modifie le régime hydrologique du bassin-versant associé, impactant les ressources en eau et les écosystèmes aquatiques qui en dépendent. En surface, il se traduit par une augmentation de la couche active et une fragilisation de la stabilité du sol mettant en danger la durabilité des infrastructures locales. Dans ce contexte, il parait indispensable d'établir une carte de suivi long-terme de l'évolution du pergélisol à l'échelle de l'arctique afin d'anticiper les dangers associés à son dégel. Cependant, l'accès difficile et les conditions extrêmes de l'Arctique complexifient grandement l'observation directe, continue et à grande échelle du pergélisol. Si certaines cartes du pergélisol existent déjà, leur résolution est faible et elles restent statiques dans le temps. Cette thèse présente des méthodes indirectes de suivi annuel du dégel du pergélisol, à l'échelle de l'arctique et avec une résolution de l'ordre du bassin-versant puis du pixel. Ces méthodes sont l'analyse de récession de bassin-versant et l'interférométrie d'images radar de surface (InSAR), appliquées à des régions pergélisolées. D'une part, l'analyse de récession reflète, indirectement et à l'échelle du bassin-versant, la dynamique du pergélisol souterrain. En effet, selon Brutsaert (2005), le temps de récession d'un bassin-versant pergélisolé serait inversement proportionnel à la profondeur de dégel du pergélisol. Après analyse de 336 bassins arctiques, cette relation de proportionnalité se révèle complexifiée par l'étendu du pergélisol, les conditions climatiques du bassin, sa topographie et ses propriétés hydrauliques. Ces complexifications hydrologiques sont simulées avec Hydrogeosphere (HGS) à travers un model conceptuel 3D de bassin-versant arctique. Les simulations montrent que l'ouverture de tàlik sous le cours d'eau ainsi que la pente des berges du cours d'eau participent à l'augmentation de la verticalité et de la connectivité des écoulements souterrains ainsi qu'à la non-linéarité du réservoir. Ces résultats limitent l'applicabilité de la méthode d'analyse de récession basée sur le modèle de Brutsaert (2005) aux bassins-versants plats et sous-couverts de pergélisol continu uniquement. Pour les autres bassins-versants arctiques, d'autres modèles d'analyse de récession doivent être développés pour pouvoir en déduire la dynamique du pergélisol. Par ailleurs, la méthode InSAR permet de traduire la vitesse déformation du terrain arctique en vitesse de dégel du pergélisol souterrain. Cette étude s'appuie sur la méthode d'analyse de Daout et al. (2017), appliquée à plus large échelle, sur une région de 390,000 km² dans les basses terres de la baie d'Hudson. Pour cette région, l'étude estime une épaisseur de couche active de 64 cm et un taux d'affaissement de terrain de 2.5 mm/an sur 2015-2021. Si ces résultats nécessitent 15 ans de données satellites pour pouvoir être traduits en taux de dégel du pergélisol, la méthode InSAR a l'avantage d'être applicable pour n'importe quel terrain arctique et ce quel que soit l'étendu du pergélisol souterrain ou sa topographie. Chacune des méthodes offre des solutions prometteuses de cartographie haute résolution de la dynamique du pergélisol sur une large couverture spatiale et temporelle. Cette thèse donne les clefs de traduction d'un événement de récession ou d'une déformation de surface en un taux de dégel du pergélisol associé. Finalement, la transdisciplinarité de cette étude permet de confronter les résultats issus des différentes approches pour mieux les valider. Si la flexibilité et le réalisme du model 3D HGS nous ont permis de comparer les résultats issus des simulations avec ceux issus des analyses de récession des 336 bassins arctiques, il serait intéressant d'en faire de même pour les résultats issus de l'analyse InSAR. / Due to polar amplification of climate change, high latitudes are warming up twice as fast as the rest of the world. This warming leads to permafrost thawing, which modifies the subsurface hydrologic regime of draining watersheds, therefore affecting water resources and related aquatic ecosystems. At the surface, permafrost thawing increases the thickness of the overlying active layer, decreases soil stability, and jeopardizes local infrastructure sustainability. In this context, building a map of long-term permafrost evolution in Arctic regions is essential to anticipate such hazards. However, direct, continuous, and large-scale monitoring of permafrost in the Arctic is very challenging because of difficult access and unfavorable environmental conditions. Even though maps of permafrost exist, their resolution is low and they are static in time. This research presents indirect methodologies to annually update maps of Arctic permafrost at catchment or pixel scales. These methods are the recession analysis of permafrost-dominated catchments and the interferometry of Arctic-surface SAR images (InSAR). On the one hand, recession analysis allows to reflect permafrost thawing dynamics since recession slopes of Arctic catchment hydrographs are linearly related to permafrost thawing depth under Brutsaert (2005) assumptions. The recession analysis of 336 Arctic catchments reveals that the linear relationship may be complexified by the extent of the permafrost itself, the climatic conditions, the landscape topography as well as the hydraulic properties of the catchment. Hydrological complexifications and mechanisms involved when permafrost extent decreases and landscape hillslope increases are simulated with a Hydrogeosphere (HGS) 3D model of a conceptual Arctic catchment. Simulations show that river-tàlik opening and valley incision both clearly increase the verticality and the connections of subsurface flow paths as well as the non-linearity of the reservoir. These results limit the applicability of the recession analysis to flat catchments underlain by continuous permafrost only. For catchments with steeper slopes and with discontinuous permafrost, other models than Brutsaert (2005) should be used to interpret permafrost dynamics from recession analysis. On the other hand, InSAR analysis translates ground surface displacement into permafrost thawing rate. The proposed method is based on previous local studies and tested on the 390,000 km² region of the Hudson Bay Lowlands. For this region, the study reveals a median active layer thickness of 64 cm and a median rate of subsidence of 2.5 mm/yr over 2015-2021. While this method requires at least 15 years to translate active layer thickness and subsidence rate into permafrost thawing rate, the method has the advantage to be applicable to any type of permafrost extent or topography. Both methods provide great opportunities to map permafrost dynamics with large temporal and spatial coverage at a fine resolution. This research provides the keys to relate Arctic river hydrograph dynamics and Arctic surface deformation to permafrost thawing rate. Finally, the research transdisciplinarity allows for method comparison and validation. The flexibility and realism of the HGS model allow to compare permafrost-related outputs derived from both simulations and real-catchment hydrograph analysis. In the future, these outputs could be compared to InSAR results.