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.

Ce travail de thèse s’est concentré sur la modélisation de l’effet d’un couvert de neige saisonnier sur le comportement thermique de grands ouvrages de génie civil construits dans les régions nordiques, ainsi qu’à son application sur des ouvrages construits avec des matériaux à grande porométrie. À l’heure actuelle, la présence de neige est souvent négligée dans les simulations numériques de régimes thérmiques, sauf pour les études portant sur les avalanches ou le régime nival qui utilisent des modèles conceptuels ou des modèles très détaillés, basés sur la physique. L’utilisation des modèles très détaillés dans l’industrie n’est pas avantageuse pour modéliser de grands ouvrages de génie civil en raison du grand nombre d'équations différentielles partielles et de paramètres parfois difficiles à estimer. De plus, l’introduction de la neige comme un milieu dans le domaine de calcul impose de grands défis pour les simulations en continu sur plusieurs années, car, comme la neige est saisonnière, le domaine qui représente la neige doit être éliminé après sa fonte. À cet égard, l’objectif principal de cette recherche est de développer un nouvel outil pour modéliser, en continu sur plusieurs années, l’effet thermique d’un couvert de neige sur des ouvrages de génie civil, qui puisse être aisément appliqué dans des cas pratiques. Pour cela, les transferts de chaleur et de masse dans la neige ont été étudiés afin d'identifier les modes de transfert qui ont une influence significative sur la température du sol. Ensuite, un outil numérique dans le logiciel FlexPDE a été établi pour modéliser le transfert de chaleur par conduction et convection, entre le sol, la neige et l'atmosphère en continu sur toute l'année. Cet outil considère également l’effet de la neige (frontière fermée ou ouverte) sur la convection dans des matériaux granulaires grossiers, et l’effet de la pluie et du changement de phase de l’eau sur le bilan d’énergie du couvert de neige. Cet outil emploie la méthode du bilan d’énergie de surface qui est considéré comme étant une condition aux limites de Neumann pour la température. L’outil a également servi à réaliser une analyse thermique d’un barrage en remblai et à démontrer la présence des cellules de convection dans l’enrochement et l’influence de la convection d’air sur l’extraction de chaleur de la fondation. Un autre objectif de cette recherche est d'établir un modèle simple et précis de la conductivité thermique pour tous les types de neige, y compris la neige artificiellement manipulée et compactée. La plupart des modèles existants de la conductivité thermique de la neige sont développés grâce à des techniques de régression qui ont l'inconvénient de ne pas respecter les limites physiques de la neige. Pour intégrer ces limites physiques dans un modèle simple, le concept de conductivité thermique relative est utilisé dans cette étude. Ce modèle est examiné avec des données publiées et réévaluées et avec les résultats de tests effectués au laboratoire de l'Université Laval. Le modèle proposé permet d’estimer la conductivité thermique de tous les types de neige, et ce, avec une grande fiabilité. / This work focused on the modeling of the effect of a seasonal snow cover on the thermal behavior of large engineering structures built in the northern regions, and its application for structures built with materials of large porometry. Nowadays, the presence of snow is often neglected in thermal numerical simulations, except for investigating the phenomenon of avalanches or production of water after snowmelt in the mountainous regions, while using conceptual models or very detailed models based on physics. In the industry, modeling a large structure by detailed models have two principals drawback. First, because of large numbers of partial differential equations and parameters, which are sometimes difficult to estimate. The introduction of snow as a medium in the computational domain also imposes great challenges for the continuous simulations for consecutive years, because seasonal snow is present only a few months a year and the domain representing the snow must be removed after melting. In this regard, the main objective of this research is to develop a new tool for the modeling of snow thermal effect on geotechnical structure continuously for several years that can also be applied in simple practical cases. For this, transfers of the heat and mass in snow are first studied to identify the most important transfer modes that significantly affect soil temperature. Then, a numerical tool using the FlexPDE software has been established to model the heat transfer by conduction and convection between soil, snow and atmosphere continuously over the entire year. This tool also considers the effect of snow (close or open boundary) on the air convection in the coarse granular materials, the rain effect and the water phase change in the energy balance of the snow cover. This tool uses the surface energy balance as the Neumann boundary condition for temperature. The tool is also served for thermal analysis of an embankment dam and demonstrates the presence of convection cells and the influence of the air convection on the heat extraction of foundation. Another purposes of this research is to establish a simple and accurate model of thermal conductivity for all types of snow, including artificially manipulated and compacted snow. Most of the existing model of the snow’s thermal conductivity are developed through regression techniques which have the drawback of not respecting the physical limits of snow. To integrate these physical limits in a simple model, the relative thermal conductivity concept is used in this study. This model is verified with published data and further validated with the results of the tests performed in the laboratory of Laval University. The proposed model estimates the thermal conductivity of all types of snow with great reliability.

Depuis plusieurs années, les changements climatiques ont une grande influence sur la dégradation du pergélisol, laquelle impacte directement la stabilité des infrastructures. Le Nunavik est la partie nord de la province de Québec où il n'y a pas de routes ou de voies ferrées reliant les villages. Le transport aérien et maritime est donc essentiel pour assurer la communication entre les villages et le reste de la province. La construction de remblais de transport affecte inévitablement le régime thermique du sol et peut entraîner une dégradation thermique du pergélisol sous-jacent, ce qui entraîne une perte importante des capacités structurelles et fonctionnelles de l'infrastructure. Pour réduire les impacts du dégel du pergélisol sur les infrastructures de transport, plusieurs techniques de protection ont été développées. Ces différentes techniques sont basées sur la réduction de l’apport de chaleur sous l’infrastructure et sur l’augmentation de l’extraction de chaleur du remblai. Ce document présente l’’analyse des résultats d'un suivi à long terme de trois sites adaptés en utilisant des données climatologiques et des données sur la température du sol présentées graphiquement avec des analyses détaillées et la validation de la stabilité thermique avec des abaques de conception. À Tasiujaq les trois méthodes de mitigation qui ont été installées dans la pente du remblai de la piste d’atterrissage sont le remblai à pente douce, le remblai à convection d'air et le drain thermique, et les résultats de l’analyse démontrent comment la pente douce a la meilleure performance, dans les conditions de climat à Tasiujaq. À la piste d'atterrissage de Puvirnituq, le projet étudie un remblai à convection d'air, et conclue que la hauteur actuelle du remblai permet de garantir la stabilité thermique en laissant une marge de sécurité pour les futurs changements climatiques. Sur la route menant à l’aéroport de Salluit, la méthode de mitigation du pergélisol à l'étude est un drain thermique selon les analyses effectuées, la stabilisation thermique n'était pas nécessaire, mais elle a été bénéfique en offrant une grande marge de sécurité à cette importante infrastructure. / For several years now, climate change has had a major influence on the degradation of permafrost, which has a direct impact on infrastructure stability. Nunavik is the northern part of the province of Quebec where there are no roads or railways connecting the villages. Air and marine transportation are therefore essential to ensure communication between the villages and the rest of the province. The construction of transportation embankments inevitably affects the thermal regime of the ground and can lead to thermal degradation of the underlying permafrost, resulting in a significant loss of structural and functional capacity of the infrastructure. To reduce the impacts of permafrost thawing on transportation infrastructure, several protection techniques have been developed. These different techniques are based on reducing the heat input under the infrastructure and increasing the heat extraction from the embankment. The analysis of the results of a long-term monitoring of three test sites is presented in this document. In Tasiujaq, the three mitigation methods that have been installed in the slope of the airstrip embankment are a gentle slope embankment, an air convection embankment and a heat drain, and the results of the analysis demonstrate how the gentle slope performs best under Tasiujaq's climate conditions. At the Puvirnituq airstrip, the project is studying an air convection embankment, finding that the current height of the embankment provides thermal stability by leaving a one-metre safety margin for future climate change. On the road to the Salluit airport, the permafrost mitigation method under study is a heat drain, which, according to the analyses performed, was not necessary in terms of heat balance, but was beneficial in generating a large safety margin for this important infrastructure.

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.

Les infrastructures de transport jouent un rôle majeur dans le développement socio-économique des régions nordiques. La construction de remblais combinés aux changements climatiques engendrent des impacts négatifs sur le pergélisol sous-jacent, causant la dégradation des infrastructures. Des techniques de mitigation ont été proposées et testées pour limiter la dégradation du pergélisol. Toutefois, il y a peu d'informations disponibles sur les procédures de conception ou sur des lignes directrices. Le but de cette recherche est de développer des outils d'ingénierie améliorés pour les techniques de stabilisation convectives, se concentrant sur le remblai à convection d'air (ACE) et le drain thermique. Plus spécifiquement, l'approche du bilan thermique est proposée pour déterminer la condition thermique de remblais conventionnels et pour permettre la sélection de la technique de mitigation la plus appropriée pour extraire la chaleur en excès transmise au sol, si le système est actuellement instable, ou pour donner un facteur de sécurité considérant une instabilité future estimée. Quatre modèles thermiques ont été développés et calibrés à partir de données mesurées sur des sites expérimentaux. Un abaque de bilan thermique pour les remblais conventionnels ainsi que plusieurs abaques de capacité d'extraction de chaleur pour l'ACE et le drain thermique ont été développés et validés utilisant, respectivement, des simulations numériques et les données de sites d'essais du Yukon et du Nord du Québec, Canada. Ces abaques permettent aux concepteurs et ingénieurs de concevoir des remblais à convection d'air et des drains thermiques optimisés pour limiter, ou même éviter, le dégel du pergélisol à des sites spécifiques. / Transportation infrastructure plays a vital role in the social and economic development of northern regions. The construction of embankments and climate change can lead to negative impacts on the underlying permafrost, causing degradation of the infrastructure. Mitigation techniques have been proposed and tested to limit permafrost degradation. However, there is limited information on the design procedures or guidelines. The purpose of this research is to develop improved engineering tools for convective stabilization techniques, focussing on air convection embankment (ACE) and heat drain. More specifically, the heat balance approach is proposed to determine the thermal condition of conventional embankments and to allow the selection of suitable mitigation techniques to extract the amount of extra heat flowing into the foundation, if the system is currently unstable, or give a safety factor considering estimated future instability. Four thermal models have been built and calibrated with field data from experimental sites. One heat balance chart for conventional embankments and, several heat extraction capacity charts for ACE and heat drain have been developed using numerical simulations and have been validated using data from test sites in Yukon and Northern Quebec, Canada. These charts allow designers and engineers to design optimized ACE and heat drain to limit or even avoid the thawing of permafrost at specific sites.

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.

L’avancée des espèces arbustives en région subarctique influence l’environnement abiotique des écosystèmes. Cependant, l’impact de la structure verticale des arbustes est encore méconnu. Notre objectif était de quantifier le développement de la structure verticale du bouleau glanduleux afin de vérifier son effet sur le couvert nival et la température au sol. Notre échantillonnage stratifié verticalement nous a permis de déterminer que la hauteur des peuplements est positivement associée à leur âge et à leur croissance. Alors que la biomasse ligneuse est maximale à la base des individus, la biomasse foliaire suit le développement en hauteur des individus. La hauteur des peuplements est la variable la plus plausible pour expliquer l’épaisseur du couvert nival et la température au sol en hiver (relation positive). Durant l’été, cette température est toutefois mieux expliquée par la biomasse ligneuse entre 30 et 40 cm de hauteur (relation négative). L’avancée des espèces arbustives aura des conséquences sur l’environnement abiotique des écosystèmes subarctiques. / In subarctic ecosystems, shrubification has great impacts on the abiotic environment. However, impacts of the vertical shrub structure are not well known. Our objective was to characterise the vertical structure of Betula glandulosa Michx., to increase our understanding of the impact of the shrub growth form on snow depth and soil temperature. Our stratified sampling allowed us to determine that shrub height was positively associated with the age of the dominant stems but also with the vertical growth rate. Moreover, woody biomass was greater in the lower stratum while foliar biomass followed vertical growth. Shrub height was the most plausible variable explaining snow thickness and soil temperature during winter. However, the woody biomass between 30 and 40 cm best explained summer temperature. According to our results, shrub expansion would have important consequences on the abiotic environment of subarctic ecosystems.

L’apport de chaleur induit par les constructions routières est susceptible d’affecter le régime thermique du pergélisol sous-jacent et de causer des dommages importants. Sur l’Alaska Highway au Yukon (Canada), des problèmes d’instabilités sont fréquemment observés et certains sont localisés aux abords des ponceaux. En autorisant la libre circulation de l’air et de l’eau à travers une conduite, des perturbations thermiques additionnelles sont infligées au sol environnant. De plus, une mauvaise gestion du drainage de l’eau aux abords des ouvrages aggrave la sévérité des instabilités. Afin de documenter ces perturbations, deux ponceaux ont été instrumentés sur l’Alaska Highway au Yukon. Le premier est situé au site expérimental de Beaver Creek; le deuxième est proche de la frontière avec l’Alaska. Les températures du sol autour des ponceaux et la température de l’air dans les ponceaux ont été enregistrées pendant une année; tandis que la température de l’eau et le débit ont été enregistrés pendant deux printemps et étés. La collecte de ces données a permis de concevoir un modèle 2D sur le logiciel TEMP/W représentant le remblai et le ponceau de Beaver Creek. En parallèle, les données ont servi à valider un modèle mathématique établi basé sur des principes physiques. Ce modèle permet de calculer le flux de chaleur transmis au sol en fonction du débit et de la température de l’eau. Par la suite, des simulations de variation de débit et de température de l’eau ont été appliquées dans le modèle mathématique. Les différents flux de chaleur correspondants ont été appliqués sur le modèle thermique 2D afin de déterminer l’influence du débit et de la température de l’eau sur les températures du sol sous le ponceau. Des simulations similaires ont été conduites en faisant varier des paramètres de conception comme la présence ou non de l’isolant. Ce mémoire présente donc trois objectifs principaux; 1) l’amélioration des connaissances relatives à l’effet de la circulation de l’eau à travers un ponceau sur le régime thermique du sol sous-jacent; 2) la documentation sur le régime thermique autour de ponceaux routiers; et 3) l’étude de l’influence du débit et de la température de l’eau sur la stabilité thermique du sol. Finalement, une méthode de conception de système de drainage, permettant de limiter l’impact d’un ponceau sur la profondeur de dégel du sol sous-jacent, est présentée dans le présent mémoire. / Heat sources induced by road construction are likely to affect the thermal regime of underlying permafrost and cause significant damage. On the Alaska Highway in Yukon (Canada), instability problems have been observed and some are located close to culverts. Culverts allow free water and air circulation through the buried pipe. Thus, additional thermal disturbances occur at that location. In addition, poor water management around the embankment exacerbates instabilities. To document these disturbances, two culverts were monitored; one is located at the Beaver Creek test site and the second is close to the U.S Border. Air temperature within and soil temperatures around the culverts were recorded for one year, while water temperature and flow rate were recorded during two thawing seasons. The data collection allowed the design of a 2D model in TEMP/W representing thermal conditions around the Beaver Creek culvert. Then, these data were used to validate an established mathematical model based on the physical principles of heat transfer. The mathematical model calculates the heat flux transmitted to the ground as a function of water temperature and flow rate. Subsequently, variations of flow rate and water temperature were applied to the mathematical model. The corresponding heat fluxes were applied to the 2D thermal model in TEMP/W to determine the influence of flow rate and water temperature on soil temperatures below the culvert. Similar simulations were conducted with design parameters, such as the presence or absence of the insulation. This thesis 1) enhanced knowledge of the effects of flow rate and water temperature through a culvert on the soil thermal regime; 2) documented the thermal regime around culverts; and 3) studied the influence of flow rate and water temperature on the thermal regime. Finally, a drainage system design method is presented which aims to minimize the impact of a culvert on the thaw depth.