Search

At high latitudes, lake-atmosphere interactions are disrupted for several months of the year by the presence of an ice cover. By isolating the water column from the atmosphere, ice, typically topped by snow, drastically alters albedo, surface roughness, and heat exchanges relative to the open water period, with major climatic, ecological, and hydrological implications. Lake models used to simulate the appearance and disappearance of the ice cover have rarely been validated with detailed in situ observations of snow and ice. In this study, we investigate the ability of the physically-based 1D Canadian Small Lake Model (CSLM) to simulate the freeze-up, ice-cover growth, and breakup of a small boreal lake. The model, driven offline by local weather observations, is run on Lake Piché, 0.15 km2 and 4 m deep (47.32°N; 71.15°W) from 25 October 2019 to 20 July 2021, and compared to observations of the temperature profile and ice and snow cover properties. Our results show that the CSLM is able to reproduce the total ice thickness (average error of 15 cm) but not the ice type-specific thickness, underestimating clear ice and overestimating snow ice. CSLM manages to reproduce snow depth (errors less than 10 cm). However, it has an average cold bias of 2°C and an underestimation of average snow density of 34 kg m−3 . Observed and model freeze-up and break-up dates are very similar, as the model is able to predict the longevity of the ice cover to within 2 weeks. CSLM successfully reproduces seasonal stratification, the mixed layer depth, and surface water temperatures, while it shows discrepancies in simulating bottom waters especially during the open water period. [ABSTRACT FROM AUTHOR]

In the boreal forest of eastern Canada, winter temperatures are projected to increase substantially by 2100. This region is also expected to receive less solid precipitation, resulting in a reduction in snow cover thickness and duration. These changes are likely to affect hydrological processes such as snowmelt, the soil thermal regime, and snow metamorphism. The exact impact of future changes is difficult to pinpoint in the boreal forest, due to its complex structure and the fact that snow dynamics under the canopy are very different from those in the gaps. In this study, we assess the influence of a low-snow and warm winter on snowmelt dynamics, soil freezing, snowpack properties, and spring streamflow in a humid and discontinuous boreal catchment of eastern Canada (47.29° N, 71.17° W; ≈ 850 m a.m.s.l.) based on observations and SNOWPACK simulations. We monitored the soil and snow thermal regimes and sampled physical properties of the snowpack under the canopy and in two forest gaps during an exceptionally low-snow and warm winter, projected to occur more frequently in the future, and during a winter with conditions close to normal. We observe that snowmelt was earlier but slower, top soil layers were cooler, and gradient metamorphism was enhanced during the low-snow and warm winter. However, we observe that snowmelt duration increased in forest gaps, that soil freezing was enhanced only under the canopy, and that snow permeability increased more strongly under the canopy than in either gap. Our results highlight that snow accumulation and melt dynamics are controlled by meteorological conditions, soil freezing is controlled by forest structure, and snow properties are controlled by both weather forcing and canopy discontinuity. Overall, observations and simulations suggest that the exceptionally low spring streamflow in the winter of 2020–2120 was mainly driven by low snow accumulation, slow snowmelt, and low precipitation in April and May rather than enhanced percolation through the snowpack and soil freezing. [ABSTRACT FROM AUTHOR]

Rain-on-snow events can cause severe flooding in snow-dominated regions. These are expected to become more frequent in the future as climate change shifts the precipitation from snowfall to rainfall. However, little is known about how winter rainfall interacts with an evergreen canopy and affects the underlying snowpack. In this study, we document 5 years of rain-on-snow events and snowpack observations at two boreal forested sites of eastern Canada. Our observations show that rain-on-snow events over a boreal canopy lead to the formation of melt–freeze layers as rainwater refreezes at the surface of the sub-canopy snowpack. They also generate frozen percolation channels, suggesting that preferential flow is favoured in the sub-canopy snowpack during rain-on-snow events. We then used the multi-layer snow model SNOWPACK to simulate the sub-canopy snowpack at both sites. Although SNOWPACK performs reasonably well in reproducing snow height (RMSE = 17.3 cm), snow surface temperature (RMSE = 1.0 °C), and density profiles (agreement score = 0.79), its performance declines when it comes to simulating snowpack stratigraphy, as it fails to reproduce many of the observed melt–freeze layers. To correct for this, we implemented a densification function of the intercepted snow in the canopy module of SNOWPACK. This new feature allows the model to reproduce 33 % more of the observed melt–freeze layers that are induced by rain-on-snow events. This new model development also delays and reduces the snowpack runoff. In fact, it triggers the unloading of dense snow layers with small rounded grains, which in turn produces fine-over-coarse transitions that limit percolation and favour refreezing. Our results suggest that the boreal vegetation modulates the sub-canopy snowpack structure and runoff from rain-on-snow events. Overall, this study highlights the need for canopy snow property measurements to improve hydrological models in forested snow-covered regions. [ABSTRACT FROM AUTHOR]

Accurately modeling the interactions between inland water bodies and the atmosphere in meteorological and climate models is crucial, given the marked differences with surrounding landmasses. Modeling surface heat fluxes remains a challenge because direct observations available for validation are rare, especially at high latitudes. This study presents a detailed evaluation of the Canadian Small Lake Model (CSLM), a one-dimensional mixed-layer dynamic lake model, in reproducing the surface energy budget and the thermal stratification of a subarctic reservoir in eastern Canada. The analysis is supported by multiyear direct observations of turbulent heat fluxes collected on and around the 85-km2 Romaine-2 hydropower reservoir (50.7°N, 63.2°W) by two flux towers: one operating year-round on the shore and one on a raft during ice-free conditions. The CSLM, which simulates the thermal regime of the water body including ice formation and snow physics, is run in offline mode and forced by local weather observations from 25 June 2018 to 8 June 2021. Comparisons between observations and simulations confirm that CSLM can reasonably reproduce the turbulent heat fluxes and the temperature behavior of the reservoir, despite the one-dimensional nature of the model that cannot account for energy inputs and outputs associated with reservoir operations. The best performance is achieved during the first few months after the ice break-up (mean error = −0.3 and −2.7 W m−2 for latent and sensible heat fluxes, respectively). The model overreacts to strong wind events, leading to subsequent poor estimates of water temperature and eventually to an early freeze-up. The model overestimated the measured annual evaporation corrected for the lack of energy balance closure by 5% and 16% in 2019 and 2020. Significance Statement: Freshwater bodies impact the regional climate through energy and water exchanges with the atmosphere. It is challenging to model surface energy fluxes over a northern lake due to the succession of stratification and mixing periods over a year. This study focuses on the interactions between the atmosphere of an irregular shaped northern hydropower reservoir. Direct measurements of turbulent fluxes using an eddy covariance system allowed the model assessment. Turbulent fluxes were successfully predicted during the open water period. Comparison between observed and modeled time series showed a good agreement; however, the model overreacted to high wind episodes. Biases mostly occur during freeze-up and breakup, stressing the importance of a good representation of the ice cover processes. [ABSTRACT FROM AUTHOR]

This study presents a probabilistic model that partitions the precipitation phase based on hourly measurements from a network of radar-based disdrometers in eastern Canada. The network consists of 27 meteorological stations located in a boreal climate for the years 2020–2023. Precipitation phase observations showed a 2-m air temperature interval between 0–4 °C where probabilities of occurrence of solid, liquid, or mixed precipitation significantly overlapped. Single-phase precipitation was also found to occur more frequently than mixed-phase precipitation. Probabilistic phase-guided partitioning (PGP) models of increasing complexity using random forest algorithms were developed. The PGP models classified the precipitation phase and partitioned the precipitation accordingly into solid and liquid amounts. PGP_basic is based on 2-m air temperature and site elevation, while PGP_hydromet integrates relative humidity. PGP_full includes all the above data plus atmospheric reanalysis data. The PGP models were compared to benchmark precipitation phase partitioning methods. These included a single temperature threshold model set at 1.5 °C, a linear transition model with dual temperature thresholds of –0.38 and 5 °C, and a psychrometric balance model. Among the benchmark models, the single temperature threshold had the best classification performance (F1 score of 0.74) due to a low count of mixed-phase events. The other benchmark models tended to over-predict mixed-phase precipitation in order to decrease partitioning error. All PGP models showed significant phase classification improvement by reproducing the observed overlapping precipitation phases based on 2-m air temperature. PGP_hydromet and PGP_full displayed the best classification performance (F1 score of 0.84). In terms of partitioning error, PGP_full had the lowest RMSE (0.27 mm) and the least variability in performance. The RMSE of the single temperature threshold model was the highest (0.40 mm) and showed the greatest performance variability. An input variable importance analysis revealed that the additional data used in the more complex PGP models mainly improved mixed-phase precipitation prediction. The improvement of mixed-phase prediction remains a challenge. Relative humidity was deemed the least important input variable used, due to consistent near water vapor saturation conditions. Additionally, the reanalysis atmospheric data proved to be an important factor to increase the robustness of the partitioning process. This study establishes a basis for integrating automated phase observations into a hydrometeorological observation network and developing probabilistic precipitation phase models. [ABSTRACT FROM AUTHOR]

As the vegetation in the Arctic changes, tundra ecosystems along the southern border of the Arctic are becoming greener and gradually giving way to boreal ecosystems. This change is affecting local populations, wildlife, energy exchange processes between environmental compartments, and the carbon cycle. To understand the progression and the implications of this change in vegetation, satellite measurements and surface models can be employed. However, in situ observational data are required to validate these measurements and models. This paper presents observational data from two nearby sites in the forest–tundra ecotone in the Tasiapik Valley near Umiujaq in Northern Quebec, Canada. One site is on a mixture of lichen and shrub tundra. The associated data set comprises 9 years of meteorological, soil and snow data as well as 3 years of eddy covariance data. The other site, 850 m away, features vegetation consisting mostly of tall shrubs and black spruce. For that location, 6 years of meteorological, soil and snow data are available. In addition to the data from the automated stations, profiles of snow density and specific surface area were collected during field campaigns. The data are available at 10.1594/PANGAEA.964743 (Domine et al., 2024). [ABSTRACT FROM AUTHOR]

Rain-on-snow events can cause severe flooding in snow–dominated regions. These are expected to become more frequent in the future as climate change shifts the precipitation from snowfall to rainfall. However, little is known about how winter rainfall interacts with an evergreen canopy and affects the underlying snowpack. In this study, we document 5 years of rain-on-snow events and snowpack observations at two boreal forested sites of eastern Canada. Our observations show that rain-on-snow events over a boreal canopy leads to the formation of melt–freeze layers as rainwater refreezes at the surface of the sub–canopy snowpack. They also generate frozen percolation channels, suggesting that preferential flow is favored in the sub–canopy snowpack during rain-on-snow events. We then used the multi–layer snow model SNOWPACK to simulate the sub–canopy snowpack at both sites. Although SNOWPACK performs reasonably well in reproducing snow height (RMSE = 17.3 cm), snow surface temperature (RMSE = 1.0 °C), and density profiles (agreement score = 0.79), its performance declines when it comes to simulating snowpack stratigraphy, as it fails to reproduce many of the observed melt–freeze layers. To correct for this, we implemented a densification function of the intercepted snow in the canopy module of SNOWPACK. This new feature allows 27 of the 32 observed melt–freeze layers induced by rain-on-snow events to be formed by the model, instead of only 18 with the original canopy module. This new model development also delays and reduces the snowpack runoff. Indeed, it triggers the unloading of dense unloaded snow layers with small rounded grains, which in turn produces fine–over–coarse transitions that limit percolation and favor refreezing. Our results show that the boreal vegetation modulates the sub–canopy snowpack structure and runoff from rain-on-snow events. Overall, this study highlights the need for canopy snow properties measurements to improve hydrological models in forested snow–covered regions. [ABSTRACT FROM AUTHOR]

Peatlands occupy around 13% of the land cover of Canada, and thus they play a key role in the water balance at high latitudes. They are well known for having substantial water loss due to evapotranspiration. Since measurements of evapotranspiration are scarce over these environments, hydrologists generally rely on models of varying complexity to evaluate these water exchanges in the global watershed balance. This study quantifies the water budget of a small boreal peatland-dominated watershed. We assess the performance of three evapotranspiration models in comparison with in situ observations and the impact of using these models in the hydrological modeling of the watershed. The study site (~1[km.sup.2]) is located in the eastern James Bay lowlands, Quebec, Canada. During summer 2012, an eddy flux tower measured evapotranspiration continuously, while a trapezoidal flume monitored streamflow at the watershed outlet. We estimated evapotranspiration with a combinational model (Penman), a radiation-based model (Priestley-Taylor), and a temperature-based model (Hydro-Quebec), and performed the hydrological modeling of the watershed with HYDROTEL, a physically based semi-distributed model. Our results show that the Penman and Priestley-Taylor models reproduce the observations with the highest precision, while a substantial drop in performance occurs with the Hydro-Québec model. However, these discrepancies did not appear to reduce the hydrological model efficiency, at least from what can be concluded from a 3-month modeling period. HYDROTEL appears sensitive to evapotranspiration inputs, but calibration of model parameters can compensate for the differences. These findings still need to be confirmed with longer modeling periods. Résumé : Les tourbières occupent environ 13 % de la couverture terrestre au Canada et jouent donc un rôle clé dans l'équilibre hydrique à hautes latitudes. Il est bien établi que les tourbières subissent d'importantes pertes en eau en raison de l'évapotranspiration. Comme des mesures de cette dernière sont rares dans ces milieux, les hydrologues ont souvent recours à des modèles de complexité variable pour évaluer ces échanges hydriques dans le bilan hydrique global d'un bassin versant. L'étude quantifie le bilan hydrique d'un petit bassin versant boreal dominé par des tourbieres. Nous évaluons la performance de trois modèles d'évapotranspiration au vu d'observations in situ et l'incidence de l'utilisation de ces modeles dans la modélisation hydrologique du bassin versant. Le site à l'étude (~1[km.sup.2]) est dans lapartie est des basses terres de la baie James (Québec, Canada). Durant l'été 2011, une tour de mesure des flux turbulents a été utilisée pour mesurer l'évapotranspiration en continu, alors que la surveillance des debits à la sortie du bassin versant était assurée par un canal jaugeur trapézoïdal. Nous avons estimé l'évapotranspiration à l'aide d'un modele combinatoire (Penman), d'un modele base sur le rayonnement (Priestley-Taylor) et d'un modèle basé sur la temperature (Hydro Quebec), et avons effectué la modélisation hydrologique du bassin versant à l'aide d'HYDROTEL, un modèle semi-distribué à base physique. Nos résultats montrent que les modèles de Penman et de Priestley-Taylor reproduisent les observations avec la plus grande précision, alors qu'une baisse de performance considérable est observée pour le modèle d'Hydro Québec. Ces discordances ne semblent toutefois pas réduire l'efficacité du modele hydrologique, du moins à la lumière des conclusions obtenues pour une période de modélisation de trois mois. HYDROTEL semble être sensible aux entrées relatives à l'évapotranspiration, mais le calage des paramètres du modèle peut compenser ces différences. Ces constatations restent à confirmer à la lumière de périodes de modélisation plus longues. [Traduit par la Rédaction], Introduction Wetlands are estimated to cover between 5 and 6 million [km.sup.2], equivalent to roughly 4% of the emerged Earth surface (Matthews and Fung 1987; Aselmann and Crutzen 1989). In [...]

In the boreal forest, winter temperatures are projected to increase substantially by 2100, resulting in a reduction in snow cover thickness and duration. These changes are likely to affect hydrological processes such as snowmelt, the soil thermal regime, and snow metamorphism. The exact impact of future changes is difficult to pinpoint in the boreal forest, due to its complex structure, and the fact that snow dynamics under the canopy are very different from those in the gaps. Although the effects of warmer winters on snow-related processes are well documented, their interactions to influence the spring runoff in evergreen forest remain poorly understood. In this observational study, we assess the influence of a low-snow and warm winter on snowmelt dynamics, soil freezing, snowpack properties, and spring streamflow in a humid and discontinuous boreal catchment of eastern Canada (47.29° N, 71.17° W, » 850 m ASL). We monitored the soil and snow thermal regimes and sampled physical properties of the snowpack under the canopy and in two forest gaps during an exceptionally low-snow and warm winter, plausibly representative of future winters, and during a winter with conditions close to normal. We observe that snowmelt was earlier but slower, top soil layers were cooler, and gradient metamorphism was enhanced during the low-snow and warm winter. However, we observe that snowmelt duration increased in forest gaps, that soil freezing was enhanced only under the canopy, and that snow permeability increased more strongly under the canopy than in either gap. Overall, we observe that the spring streamflow discharge was significantly reduced in the warmest year due to a slower melt and low precipitation in April and May. Our results, based solely on field observations, highlight the complex effects of warmer winters on snow hydrology in discontinuous boreal forests. [ABSTRACT FROM AUTHOR]

As the vegetation in the Arctic changes, tundra ecosystems along the southern border of the Arctic are becoming greener and gradually giving way to boreal ecosystems. This change is affecting local populations, wildlife, energy exchange processes between environmental compartments, and the carbon cycle. To understand the progression and the implications of this change in vegetation, satellite measurements and surface models can be employed. However, in situ observational data are required to validate these measurements and models. This paper presents observational data from two nearby sites in the forest-tundra ecotone in the Tasiapik Valley near Umiujaq in Northern Quebec, Canada. One site is on a mixture of lichen and shrub tundra. The associated data set comprises 9 years of meteorological, soil and snow data as well as 3 years of eddy covariance data. The other site, 850 m away, features vegetation consisting mostly of tall shrubs and black spruce. For that location, 6 years of meteorological, soil and snow data are available. In addition to the data from the automated stations, profiles for snow density and specific surface area were collected during field campaigns. The data are available at https://doi.org/10.1594/PANGAEA.946538 (Lackner et al., 2022b). [ABSTRACT FROM AUTHOR]