Snow cover estimation is a principal source of error for spring streamflow simulations in Québec, Canada. Optical and near infrared remote sensing can improve snow cover area (SCA) estimation due to high spatial resolution but is limited by cloud cover and incoming solar radiation. Passive microwave remote sensing is complementary by its near-transparence to cloud cover and independence to incoming solar radiation, but is limited by its coarse spatial resolution. The study aims to create an improved SCA product from blended passive microwave (AMSR-E daily L3 Brightness Temperature) and optical (MODIS Terra and Aqua daily snow cover L3) remote sensing data in order to improve estimation of river streamflow caused by snowmelt with Québec's operational MOHYSE hydrological model through direct-insertion of the blended SCA product in a coupled snowmelt module (SPH-AV). SCA estimated from AMSR-E data is first compared with SCA estimated with MODIS, as well as with in situ snow depth measurements. Results show good agreement (+95%) between AMSR-E-derived and MODIS-derived SCA products in spring but comparisons with Environment Canada ground stations and SCA derived from Advanced Very High Resolution Radiometer (AVHRR) data show lesser agreements (83 % and 74% respectively). Results also show that AMSR-E generally underestimates SCA. Assimilating the blended snow product in SPH-AV coupled with MOHYSE yields significant improvement of simulated streamflow for the aux Écorces et au Saumon rivers overall when compared with simulations with no update during thaw events. These improvements are similar to results driven by biweekly ground data. Assimilation of remotely-sensed passive microwave data was also found to have little positive impact on springflood forecast due to the difficulty in differentiating melting snow from snow-free surfaces. Considering the direct-insertion and Newtonian nudging assimilation methods, the study also shows the latter method to be superior to the former, notably when assimilating noisy data.

Memoria para optar al título de Geóloga / Durante el periodo de vida de los mantos estacionales, parte de la nieve superficial de éste es liberado hacia la atmósfera en forma de flujos de vapor de agua, mediante el proceso de sublimación (Gustafson, 2010), constituyendo una fracción del recurso hídrico almacenado que no se encontrará disponible para ser fundido durante su deshielo. El factor de sublimación es uno de los componentes menos estudiados dentro de los balances hidrológicos de cuencas nivales en Chile. El presente trabajo se realizó en la mina El Teniente ubicada en una cuenca pluvio-nival, en la cual las pérdidas asociadas a la sublimación del manto nival y la evaporación de sus remanentes, constituye el factor de mayor incertidumbre dentro de los balances hidrológicos realizados en la zona. Es por esto, que el presente estudio tiene por objetivo estimar las pérdidas del manto nival, mediante el análisis de sus signaturas isotópicas para la estimación de sublimación y aplicando la metodología del domo para la estimación de evaporación. Gustafson et al., (2010) sugiere que partir de los valores isotópicos del manto de nieve en el momento de máxima acumulación se puede estimar la sublimación durante el invierno Se realizó una primera campaña de terreno entre el 12 de septiembre y 19 de octubre del 2017, en la que se recolectaron 49 muestras de nieve para análisis de isótopos estables y se realizaron mediciones de densidad y altura de nieve. Y una segunda, entre el 07 de noviembre y 26 de diciembre del 2017 en la que se aplicó la metodología del domo para determinar la evaporación. Sumado a lo anterior, se analizaron los datos de las estaciones nivométricas y meteorológicas, para determinar el valor de densidad de nieve fresca y la influencia de los factores meteorológicos en la ocurrencia de los procesos de sublimación y evaporación. A partir de las mediciones físico-químicas realizadas en terreno y asumiendo que entre las nevadas del 30 de septiembre y 04 de octubre el manto no experimentó fusión debido a las bajas temperaturas en la cuenca. Se determina que la densidad de nieve fresca presenta una variabilidad estacional, registrándose durante el invierno promedios mensuales menores que en primavera. El promedio anual de nieve fresca posee un valor de 152 Kg/m^3 y 114 Kg/m^3 para los años 2016 y 2017, respectivamente. Posterior a la depositación de nieve, ésta incrementa su densidad hasta alcanzar valores cercanos a los 500 Kg/m^3. Se estimó un porcentaje de sublimación y evaporación de 25% y 3%, respectivamente, en función de la precipitación anual de nieve. Finalmente, se concluye que el equivalente de nieve en agua de las precipitaciones de nieve han sido subestimadas en los balances hidrológicos realizados anteriormente y que las pérdidas por sublimación del manto de nieve y la evaporación de sus remanentes corresponden a un 28% de la precipitación de nieve anual del año 2017, constituyendo la sublimación el factor preponderante.

Résumé : Cette recherche présente l'analyse spatio-temporelle de l'évolution conjointe de l'augmentation de la présence arbustive et de la dynamique de la fonte de la neige au Nunavik, Québec, Canada. Cette zone est caractérisée par la complexité de l'interaction de multiples changements simultanés de la température, de la couverture de la neige ainsi que de la pousse végétative. La première partie de ce travail consiste à faire l'analyse de l’évolution temporelle de ces multiples changements. Cette analyse a nécessité la mise en place d’une importante base de données climatiques, satellitaires et de couverture de sol à plusieurs échelles, sur une période allant jusqu'à 60 ans, soit de 1950-2012. La deuxième partie du travail consiste à faire l'analyse spatiale à haute résolution de l’influence de la fraction du couvert arbustif sur la fonte de la neige. L'analyse et l'interprétation des résultats obtenus dans la première partie montrent clairement un changement climatique significatif sur la région étudiée, découpée en 3 bandes de latitude correspondant à la toundra, la taïga ouverte et à la taïga forestière, respectivement du Nord au Sud. Ce changement de climat correspond à un réchauffement marqué, entre 0.75°C et 1.57°C par décade entre les zones 1 (toundra) et 3 (taïga forestière) respectivement. On peut noter que la hauteur de neige maximale annuelle a diminué dans les trois zones alors que les précipitations hivernales ont augmenté en zone 1 et 3 sur les 45 dernières années. Les résultats montrent une nette augmentation de la végétation arbustive dans les zones 2 et 3 (LAI plus élevé de 100% dans la zone 3 par rapport à la zone 1). L'impact de la végétation a été analysé à partir de la durée de fonte relative entre le début de la fonte et la disparition de la neige. Il apparait clairement que la végétation active la fonte précocement, allongeant ainsi significativement la durée de fonte (+600%). Cependant, l'impact de la végétation ne retarde pas la date de fin du couvert nival qui est de plus en plus précoce pour les zones 2 et 3. L'analyse spatiale à haute résolution montre que la présence arbustive entraine une date de fin de neige plus précoce par rapport au sol nu. Cette étude démontre clairement que la croissance de la végétation qui résulte du réchauffement climatique impacte la dynamique du couvert nival, aussi affectée par ce réchauffement. Une étude approfondie des processus en causes avec des mesures in situ appuyées par leur modélisation permettrait de mieux comprendre ces phénomènes. // Abstract : This study presents a spatial-temporal analysis of the joint evolution of the increase of shrubiness and the dynamics of snowmelt in Nunavik, Quebec, Canada. This zone is characterized by the complexity of the interaction of multiple changes of temperature, snow cover and vegetation growth. The first part of this study analyzes the temporal evolution of these changes. The analysis required the use of a large database on climate, satellite data and ground cover at multiple scales over a period of up to 60 years, from 1950 to 2012. The second part of the study consists of a spatial high-resolution analysis of the influence of the fraction of shrub cover on snowmelt. The analysis and interpretation of the results clearly show a significant climate change over the study area, divided into three latitudinal transects corresponding to tundra, open taiga and forested taiga. A significant warming of 0.75 ° C and 1.57 ° C per decade was experienced between zones 1 (tundra) and 3 (forested taiga) respectively. The maximum annual snow depth on the ground decreased over the 3 zones studied while winter precipitations increased in zones 1 and 3 over the last 45 years. The results show a significant increase in shrub vegetation in zones 2 and 3. The impact of the vegetation on snow was analyzed with melt duration (from melt onset to complete melt). It appears clearly that the vegetation triggers the melting process earlier and significantly extends the melt duration (+600%). However, the impact of vegetation does not delay the date of the snow cover disappearance. The high-resolution spatial analysis showed that shrubs cause an earlier snow cover disappearance date than bare soil. This study clearly demonstrates that vegetation growth resulting from global warming impacts the snow cover dynamics, which are also affected by global warming. A thorough study of the processes with in-situ measurements supported by models would help gaining a better comprehension of these phenomena.

Le caribou de Peary est l’unité désignable du caribou la plus septentrionale ; sa population a chuté d’environ 70% au cours des trois dernières générations. Le Comité sur la situation des espèces en péril au Canada (COSEPAC) identifie les conditions difficiles d’accès à la nourriture à travers le couvert nival comme le facteur le plus influant contribuant à ce déclin. Cette étude se concentre sur l’établissement d’un outil spatial de caractérisation des conditions nivales pour l’accès à la nourriture du caribou de Peary dans le Nord canadien, utilisant des simulations du couvert nival générées avec le logiciel suisse SNOWPACK à partir des données du Modèle Régional Climatique Canadien. Le cycle de vie du caribou de Peary a été divisé en trois périodes critiques : la période de mise bas et de migration printanière (avril – juin), celle d’abondance de nourriture et de rut (juillet – octobre) et celle de migration automnale et de survie des jeunes caribous (novembre – mars). Les conditions nivales sont analysées et les simulations du couvert nival comparées aux comptes insulaires de caribous de Peary pour identifier un paramètre nival qui agirait comme prédicateur des conditions d’accès à la nourriture et expliquerait les fluctuations des comptes de caribous. Cette analyse conclue que ces comptes sont affectés par des densités de neige au-dessus de 300 kg/m³. Un outil logiciel cartographiant à une échelle régionale (dans l’archipel arctique canadien) les conditions d’accès à la nourriture possiblement favorables et non favorables basées sur la neige est proposé. Des exemples spécifiques de sorties sont données pour montrer l’utilité de l’outil, en cartographiant les pixels selon l’épaisseur cumulée de neige au-dessus de densités de 300 kg/m³, où des épaisseurs cumulées au-dessus de 7000 cm par hiver sont considérées comme non favorables pour le caribou de Peary.

Résumé : Cette recherche présente l'analyse spatio-temporelle de l'évolution conjointe de l'augmentation de la présence arbustive et de la dynamique de la fonte de la neige au Nunavik, Québec, Canada. Cette zone est caractérisée par la complexité de l'interaction de multiples changements simultanés de la température, de la couverture de la neige ainsi que de la pousse végétative. La première partie de ce travail consiste à faire l'analyse de l’évolution temporelle de ces multiples changements. Cette analyse a nécessité la mise en place d’une importante base de données climatiques, satellitaires et de couverture de sol à plusieurs échelles, sur une période allant jusqu'à 60 ans, soit de 1950-2012. La deuxième partie du travail consiste à faire l'analyse spatiale à haute résolution de l’influence de la fraction du couvert arbustif sur la fonte de la neige. L'analyse et l'interprétation des résultats obtenus dans la première partie montrent clairement un changement climatique significatif sur la région étudiée, découpée en 3 bandes de latitude correspondant à la toundra, la taïga ouverte et à la taïga forestière, respectivement du Nord au Sud. Ce changement de climat correspond à un réchauffement marqué, entre 0.75°C et 1.57°C par décade entre les zones 1 (toundra) et 3 (taïga forestière) respectivement. On peut noter que la hauteur de neige maximale annuelle a diminué dans les trois zones alors que les précipitations hivernales ont augmenté en zone 1 et 3 sur les 45 dernières années. Les résultats montrent une nette augmentation de la végétation arbustive dans les zones 2 et 3 (LAI plus élevé de 100% dans la zone 3 par rapport à la zone 1). L'impact de la végétation a été analysé à partir de la durée de fonte relative entre le début de la fonte et la disparition de la neige. Il apparait clairement que la végétation active la fonte précocement, allongeant ainsi significativement la durée de fonte (+600%). Cependant, l'impact de la végétation ne retarde pas la date de fin du couvert nival qui est de plus en plus précoce pour les zones 2 et 3. L'analyse spatiale à haute résolution montre que la présence arbustive entraine une date de fin de neige plus précoce par rapport au sol nu. Cette étude démontre clairement que la croissance de la végétation qui résulte du réchauffement climatique impacte la dynamique du couvert nival, aussi affectée par ce réchauffement. Une étude approfondie des processus en causes avec des mesures in situ appuyées par leur modélisation permettrait de mieux comprendre ces phénomènes. // Abstract : This study presents a spatial-temporal analysis of the joint evolution of the increase of shrubiness and the dynamics of snowmelt in Nunavik, Quebec, Canada. This zone is characterized by the complexity of the interaction of multiple changes of temperature, snow cover and vegetation growth. The first part of this study analyzes the temporal evolution of these changes. The analysis required the use of a large database on climate, satellite data and ground cover at multiple scales over a period of up to 60 years, from 1950 to 2012. The second part of the study consists of a spatial high-resolution analysis of the influence of the fraction of shrub cover on snowmelt. The analysis and interpretation of the results clearly show a significant climate change over the study area, divided into three latitudinal transects corresponding to tundra, open taiga and forested taiga. A significant warming of 0.75 ° C and 1.57 ° C per decade was experienced between zones 1 (tundra) and 3 (forested taiga) respectively. The maximum annual snow depth on the ground decreased over the 3 zones studied while winter precipitations increased in zones 1 and 3 over the last 45 years. The results show a significant increase in shrub vegetation in zones 2 and 3. The impact of the vegetation on snow was analyzed with melt duration (from melt onset to complete melt). It appears clearly that the vegetation triggers the melting process earlier and significantly extends the melt duration (+600%). However, the impact of vegetation does not delay the date of the snow cover disappearance. The high-resolution spatial analysis showed that shrubs cause an earlier snow cover disappearance date than bare soil. This study clearly demonstrates that vegetation growth resulting from global warming impacts the snow cover dynamics, which are also affected by global warming. A thorough study of the processes with in-situ measurements supported by models would help gaining a better comprehension of these phenomena.

No description available.

Long-duration records are necessary to understand and assess the long-term dynamics of natural systems. The purpose of this research was to use dendrochronologic modelling to construct proxy histories of hydroclimatic conditions and Pacific salmon abundance in west central British Columbia. A multi-species regional network of tree ring-width and ring-density measurements was established from new and archived tree-ring chronologies. These chronologies were then used in multivariate linear regression models to construct proxy records of nival river discharge, summer temperature, end-of-winter snow-water equivalent (SWE), the winter Pacific North America pattern (PNA) and Pacific salmon abundance. All proxy hydroclimate records provide information back to 1660 AD. Reconstructions of July-August mean runoff for the Skeena and Atnarko rivers describe below average conditions during the early- to mid-1700s and parts of the early-, mid- and late-1900s. Models describe intervals of above average river discharge during the late-1600s, the early-1700s and 1800s, and parts of the early- and mid-1900s. Fluctuations in proxy reconstructions of July-August mean temperature for Wistaria and Tatlayoko Lake, May 1 SWE at Mount Cronin and Tatlayoko Lake and October-February PNA occurred in near synchrony with the shifts described in runoff records. Episodes of above average runoff were typically associated with periods of enhanced end-of-winter SWE, below average summer temperature and positive winter PNA. A history of Pacific salmon abundance was reconstructed for four species of salmon (chinook, sockeye, chum and pink) that migrate to coastal watersheds of west central British Columbia. Proxy records vary in length and extend from 1400 AD, 1536 AD and 1638 AD to present. Salmon abundance reconstructions varied throughout the past six centuries and described significant collapse in population levels during the early-1400s, the late-1500s, the mid-1600s, the early-1700s, the early-1800s and parts of the 1900s. Wavelet analyses of reconstructed hydroclimate and salmon population records revealed low- and high-frequency cycles in the data. Correlation analyses related reconstructions to atmospheric teleconnection indices describing variability in North Pacific sea surface temperatures and the Aleutian Low pressure centre. To a lesser degree, relationships were also established between reconstructions and the El Niño-Southern Oscillation. Results thus confirm the long-term influence of large-scale ocean and atmospheric circulation patterns on hydroclimate and Pacific salmon abundance in west central British Columbia. The reconstructions introduced in this thesis provide insights about the long-term dynamics of the west central British Columbia environment. Several reconstructions presented in this thesis provide novel contributions to dendrohydroclimatic and paleoecologic research in Pacific North America. Proxy runoff records for the Skeena and Atnarko rivers are the first to be constructed for nival-regime basins in British Columbia. The models of Skeena River runoff and Mount Cronin SWE are additionally the first reconstructions of runoff and snowpack in Pacific North America based on a ring-density chronology, demonstrating the significant contribution that wood density measurements can make to dendrohydroclimate research. The models of Pacific salmon stocks are the first to utilize climate-sensitive tree-ring records to construct a history of regional salmon abundance and thus represent a significant advancement to paleoecological modelling. / Graduate

Situés dans une région clé pour la dynamique du climat global, les glaciers tropicaux sont de très bons indicateurs de la variabilité climatique. Ce travail a permis de spatialiser les modèles de bilan d'énergie CROCUS (pour la neige et la glace) et ISBA (pour le sol non enneigé) sur le bassin versant du glacier Zongo (16°S, 3.3 km2, Bolivie), englacé à 63 % et suivi par le SO/ORE GLACIOCLIM. Ainsi, des jeux de données nivo-météorologiques complets ont été constitués du 01/09/2004 au 31/03/2006 en accordant une attention particulière aux lames d'eau et aux phases des précipitations. Le modèle CROCUS-ISBA a été adapté pour simuler le manteau neigeux éphémère de la moraine et les résultats ont été validés sur plusieurs sites de mesures. Sur le glacier, le modèle CROCUS a permis de simuler correctement le bilan de masse local mesuré par les balises d'ablation. Ensuite, les modèles ont été étendus à l'ensemble du bassin versant et alimentés par les données météorologiques spatialisées. Le bilan net spécifique du glacier simulé est validé par celui évalué par photogrammétrie. Les débits à l'exutoire du bassin versant et leur saisonnalité sont bien restitués par le modèle (8% d'erreur sur les volumes d'eau pour toute la période, critère de Nash quotidien de 0.63). De nombreux résultats sont présentés sur le coefficient d'écoulement de la moraine, la contribution du glacier par tranche d'altitude et la réponse du bassin versant à deux saisons humides contrastées. Cet outil permet ainsi de mieux évaluer et comprendre les processus responsables de la fonte sur le glacier et sur les moraines et de quantifier leurs contributions respectives à la ressource en eau locale.

Chaque jour, des décisions doivent être prises quant à la quantité d'hydroélectricité produite au Québec. Ces décisions reposent sur la prévision des apports en eau dans les bassins versants produite à l'aide de modèles hydrologiques. Ces modèles prennent en compte plusieurs facteurs, dont notamment la présence ou l'absence de neige au sol. Cette information est primordiale durant la fonte printanière pour anticiper les apports à venir, puisqu'entre 30 et 40% du volume de crue peut provenir de la fonte du couvert nival. Il est donc nécessaire pour les prévisionnistes de pouvoir suivre l'évolution du couvert de neige de façon quotidienne afin d'ajuster leurs prévisions selon le phénomène de fonte. Des méthodes pour cartographier la neige au sol sont actuellement utilisées à l'Institut de recherche d'Hydro-Québec (IREQ), mais elles présentent quelques lacunes. Ce mémoire a pour objectif d'utiliser des données de télédétection en micro-ondes passives (le gradient de températures de brillance en position verticale (GTV)) à l'aide d'une approche statistique afin de produire des cartes neige/non-neige et d'en quantifier l'incertitude de classification. Pour ce faire, le GTV a été utilisé afin de calculer une probabilité de neige quotidienne via les mélanges de lois normales selon la statistique bayésienne. Par la suite, ces probabilités ont été modélisées à l'aide de la régression linéaire sur les logits et des cartographies du couvert nival ont été produites. Les résultats des modèles ont été validés qualitativement et quantitativement, puis leur intégration à Hydro-Québec a été discutée. / Every day, decisions must be made about the amount of hydroelectricity produced in Quebec. These decisions are based on the prediction of water inflow in watersheds based on hydrological models. These models take into account several factors, including the presence or absence of snow. This information is critical during the spring melt to anticipate future flows, since between 30 and 40 % of the flood volume may come from the melting of the snow cover. It is therefore necessary for forecasters to be able to monitor on a daily basis the snow cover to adjust their expectations about the melting phenomenon. Some methods to map snow on the ground are currently used at the Institut de recherche d'Hydro-Québec (IREQ), but they have some shortcomings. This master thesis's main goal is to use remote sensing passive microwave data (the vertically polarized brightness temperature gradient ratio (GTV)) with a statistical approach to produce snow maps and to quantify the classification uncertainty. In order to do this, the GTV has been used to calculate a daily probability of snow via a Gaussian mixture model using Bayesian statistics. Subsequently, these probabilities were modeled using linear regression models on logits and snow cover maps were produced. The models results were validated qualitatively and quantitatively, and their integration at Hydro-Québec was discussed.

Le biome boréal, emmagasinant d’importantes quantités de carbone en son sol et recouvrant une majorité du territoire alaskien, fennoscandien et russe, contribue grandement au système climatique. Toutefois, les variabilités climatiques et les propriétés de l’écosystème, notamment en ce qui a trait à la présence ou l’absence de pergélisol, complexifient la quantification de la variabilité des bilans de carbone du biome boréal, au sein duquel se retrouvent des écosystèmes forestiers, lentiques et de zones humides. Ces bilans de carbone sont grandement influencés par le début et la fin de la saison de croissance photosynthétique, étant à leur tour dépendants de plusieurs variables environnementales telles que la température de l’air et du sol, le contenu du sol en eau, les stades de développement de la végétation, etc. Cette recherche vise à quantifier l’impact de ces variabilités environnementales sur la variabilité des moments où se produisent le début et la fin de la saison de croissance photosynthétique, en distinguant les forêts boréales avec et sans pergélisol. La saison de croissance photosynthétique est caractérisée à partir de la productivité primaire brute dérivée de mesures covariance des turbulences provenant de 40 sites-années d’observation à travers la forêt boréale nord-américaine où l’épinette noire est l’espèce d’arbre dominante. Les variables environnementales considérées étaient les températures de l’air et du sol, les stades de développement de la végétation, le couvert nival, le rayonnement photosynthétiquement actif et le contenu du sol en eau. Le cadre statistique choisi incluait le calcul des coefficients de corrélations de Pearson, l’analyse des points communs et la modélisation par équations structurelles. Les résultats de cette étude montrent que la variabilité du début de la saison de croissance dans les sites sans pergélisol est contrôlée directement par la variabilité annuelle des stades de développement de la végétation ainsi que par le moment où survient le dégel du sol. Ce résultat souligne ainsi l’importance de l’accès à l’eau liquide du sol afin que la végétation initie la photosynthèse. Aucune variable environnementale ne pouvait significativement expliquer le contrôle du début de la photosynthèse au sein des sites avec pergélisol. À l’automne, le contenu du sol en eau ainsi que le début du couvert nival influencent directement la variabilité de la fin de la saison de croissance photosynthétique. Il est alors montré que la disponibilité de l’eau peut mener à une cessation plus hâtive de la photosynthèse à l’automne. L’effet de l’apparition du couvert nival est quant à lui opposé dans les sites avec et sans pergélisol. Son retard dans les sites sans pergélisol témoigne d’une température de l’air suffisamment élevée pour que les précipitations tombent sous forme liquide, prolongeant ainsi les activités photosynthétiques. Son retard dans les sites avec pergélisol signifie plutôt des précipitations neigeuses moindres, retardant ainsi l’apparition d’une couche isolante pour le sol, qui aurait pu allonger la saison de croissance photosynthétique. Cette étude contribue à clarifier les processus contrôlant le début et la fin de la saison de croissance photosynthétique et aidera à améliorer la compréhension des effets des changements climatiques sur la force du puits de carbone de la forêt boréale nord-américaine. / The boreal forest, storing large amounts of carbon in its soil and covering a majority of the Alaskan, Canadian, Fennoscandian and Russian territory, is an integral part of the climate system. However, climatic variability and ecosystem properties, particularly with regards to the presence or absence of permafrost, limits our understanding of the carbon balance variability in the boreal biome, which comprises forest, lake and wetland ecosystems. The boreal carbon sink-source strength is greatly influenced by phenological events, including the start and end of the photosynthetic growing season, which are themselves dependent on several environmental variables such as air and soil temperature, soil water content, vegetation development stages, etc. This research aims to provide new insights on the influence of environmental variability on the variability in the timing of the photosynthetic growing season, by broadly distinguishing between boreal forests with and without permafrost. The photosynthetic growing season is characterized using gross primary productivity derived from eddy covariance measurements of net ecosystem carbon dioxide exchange. Data from 40 black spruce- dominated site-years of observation across the North American boreal forest are used. The considered environmental predictors were air and soil temperatures, vegetation development stages, snow cover, photosynthetically active radiation and soil water content. The statistical framework included the calculation of Pearson correlation coefficients, commonality analyses and structural equation modeling. This study shows that the variability in the start of the growing season in permafrost-free sites is directly controlled by the variability in vegetation development stage as well as by the thawing of seasonally frozen ground. This result thus emphasizes the importance of access to liquid soil water for the vegetation to initiate photosynthesis. No environmental variable could significantly explain photosynthesis recovery in sites with permafrost. In fall, the soil water content as well as the start of snow cover directly influence the variability in the end of the photosynthetic growing season. These results suggest that the availability of water can limit photosynthesis in the fall. The effect of snow cover is opposite in sites with and without permafrost. A delay in the appearance of continuous snow cover in sites without permafrost indicates that the air temperature is high enough for precipitation to fall in liquid form and for photosynthesis to continue. In contrast, its delay in sites with permafrost indicates less snowfall, thus delaying the appearance of an insulating layer for the soil, which could have lengthened the photosynthetic growing season. This study sheds light on the controls of the annual variation of the timing of the photosynthetic growing season and will help understanding of the effects of climate change on the strength of the North American boreal forest carbon sink.