ABSTRACTMethods used to quantify evapotranspiration (ET) from Sphagnum‐dominated peatlands often assume that soil moisture is not a limiting factor; actual ET (AET) equals potential ET (PET). However, soil moisture can become limiting as peatlands dry, lowering AET below PET and necessitating the use of a surface resistance term in AET estimations. Quantifying and calculating surface resistance is a challenge for the non‐vascular plant surfaces such as those dominated by Sphagnum moss. This paper explores and quantifies the ecohydrological processes that drive Sphagnum resistance to ET. It is hypothesized that a relationship exists between the Sphagnum moss resistance and the ratio of unsaturated to saturated hydraulic conductivity (K‐ratio) for boreal peatlands, where the K‐ratio is a proxy for the hydrophysical properties of the porous medium. An empirical relationship between Sphagnum moss resistance and the K‐ratio was developed from data collected from a boreal peatland and implemented in the cold regions hydrological model. Empirically modelled resistance values (0–800 s m−1) did not match well with estimates from inverting observations and the Penman–Monteith (PM) algorithm (0–5000 s m−1). Difficulties in validating resistance values were possibly due to lack of moisture limiting conditions although this is seemingly contradicted by the alpha value being less than 1. Priestley–Taylor (PT) and PM algorithms in CRHM were used to estimate AET and compared with each other and with observations from an onsite eddy covariance (EC) system. The PT algorithm, using a site‐specific alpha value (0.75) performed the best with a mean difference of 9.4% (±12.0%) when compared to EC measurements of AET. The PM algorithm consistently overestimated EC measurements with a mean difference of 68.4% (±50.0%), even with a moss resistance incorporated into its use. The performance of PM algorithm is impeded by the uncertainty in quantifying Sphagnum resistance. Reducing this uncertainty should be a focus of future studies, as it does not require the use of a site‐specific alpha value.

AbstractEnergy absorption and flow through a nest is an important aspect of embryonic development in many reptile species including turtles. To date, few studies have explicitly attempted to quantify the energy flow through turtle nests, opting instead for the simplified approach offered by temperature index models. However, the quantification of the energy can provide an explicit abiotic link that can link biological models to biometeorological and ecohydrological processes and models. We investigated the energy flow through turtle nests occupying different bedrock morphologies within a Canadian Shield Rock Barren landscape, in Ontario, Canada. The taxons studied were Spotted Turtle (Clemmys guttata), Midland Painted Turtle (Chrysemys picta marginata), and Blanding's Turtle (Emydoidea blandingii). Nest temperature and soil moisture were measured in 2018 and 2019 using sensors placed in the soil adjacent to 12 turtle nest cavities. Three main rock morphologies were identified for each nest location, Crevice, Ledge, and Flat types, that are in order of decreasing bedrock percentage contact with the nest site. Ground heat flux and change in heat storage were determined using the calorimetric method for each nest, while the direction of energy flux between the atmosphere and the underlying rock was also determined. The Crevice nest morphology experienced the lowest ground heat flux on average (1.56 × 10−1 W m−2) and lowest cumulative heat storage (230 MJ) compared to the Flat (440 MJ) and Ledge (331 MJ) nests. However, over the diurnal cycle, large heat gains by Flat nests were mostly balanced out by nighttime heat losses. While Crevice nests saw the lowest daily heat storage gains, they experienced much lower heat losses over the evening period compared to the other nest types. Furthermore, we found that 59% of the energy is directed from the underlying bedrock into the Crevice nest, highlighting the importance of the bedrock in controlling thermal dynamics in the turtle nesting habitat. The lower variability in energy parameters for Crevice nest types can be attributed to higher amounts of nest‐to‐bedrock contact, compared to the flat nest types. Our results indicate that Crevice morphology may be ideal for turtles nesting at their northern limits because minimal heat loss during the evening can result in a more stable thermal incubation environment. Future conservation and habitat restoration efforts should consider the importance of bedrock morphology and prioritize the protection of Crevice nest sites. Furthermore, this work highlights important opportunities for potential interdisciplinary work between ecologists, climatologists, biologists, and hydrologists, specifically the integration of ecohydrological and biological models. This work also underscores the potential uncertainty of climate change impacts on turtle egg hatching success and nest sex ratios.

Bare soil evaporation has been studied extensively, but less is certain regarding how site‐specific features, especially the overstory tree canopy and ground covers, mediate evaporation processes. Inspired by recent advances on modelling bare soil evaporative efficiency (SEE), this study explored SEE over a range of soil substrates and ground cover types, with and without the presence of an overstory canopy in three mesic ecosystems in Canadian Rocky Mountains. A significant relationship was found between the critical soil water content and ground cover mass fractions across various ground cover types, both with and without the presence of an overstory canopy. This relationship is expected to be prevalent across various ecosystems. Moreover, a simple approach for modelling SEE of vegetated surfaces and a correction method to account for below‐canopy SEE is also proposed. The model yields satisfactory simulations, and the approach is expected to be widely applicable, given the strength that its parameters are easily acquired, and its formulations are simple and straightforward. While the model may be particularly suited to mesic ecosystems, the underlying mechanism of SEE suggests that this model can also be applied in dryer conditions. This approach will greatly improve ET parameterization in land‐surface models (LSMs) and increase our knowledge of the global water cycle and ecosystem responses under climate change impacts.

The modeling of seasonal ground ice (SGI) freeze/thaw a common feature in boreal peatlands, has often been completed using a unidirectional approach, where melting is driven by energy inputs from the surface. However, bi-directional melt is known to occur, and can potentially increase the spring melt rate. Accurate modelling of the timing of ice-free conditions in peatlands is important because SGI can impede spring infiltration and lead to substantial spring snowmelt runoff from peatlands. However, when modelling melt only from above, erroneous results in the model estimation of ice-free conditions can occur, which can lead to knock on-effects for modelling peatland hydrological function. Furthermore, as the climate warms, it is unclear how this role of SGI may change in the future. This study used the Stefan Equation to model unidirectional and bi-directional melt to assess which performed better in modelling the timing of ice-free conditions compared to observed values (BI: 3.9 +/- 2.1 days, UNI: 9.0 +/- 4.7). Including bi-directional melt improved model performance by reducing this difference by approximately 5 days. Model performance for SGI freeze/thaw cycles were similar, with BI being slightly more accurate in freezing (RMSE:2.7 cm versus 3.3 cm) and melting (RMSE: 2.6 cm vs 3.7 cm) compared to the unidirectional approach. While the model improvement in the timing of ice-free conditions was substantial, careful consideration is needed in determining when a peatland is functionally ice free in future modelling studies. The Stefan Equation was found to be most sensitive to changes in soil moisture, compared to ground surface temperature and peat porosity, likely due to the relationship between thermal conductivity and frozen and liquid water content. Comparisons with future climate change projections suggest that the timing of ice-free conditions 'could shift by as much as 2 weeks earlier in the 2050's and by almost a month earlier in the 2080's. However, the timing of snowfall, and rain on snow events continues to be a source of model uncertainty. Future studies should work to investigate the potential positive feedbacks this could create. In conclusion, the Stefan Equation presents a relatively easy path for incorporating bi-directional melt into peatland models. This process should be included in peatland ecohydrological models in order to properly model the timing of melt and ice-free conditions.

Surface energy budgets are important to the ecohydrology of complex terrain, where land surfaces cycle in and out of shadows creating distinct microclimates. Shading in such environments can help regulate downstream flow over the course of a growing season, but our knowledge on how shadows impact the energy budget and consequently ecohydrology in montane ecosystems is very limited. We investigated the influence of horizon shade on the surface energy fluxes of a subalpine headwater wetland in the Canadian Rocky Mountains during the growing season. During the study, surface insolation decreased by 60% (32% due to evolving horizon shade and 28% from seasonality). The influence of shade on the energy budget varied between two distinct periods: (1) Stable Shade, when horizon shade was constant and reduced sunlight by 2 h per day; and (2) Dynamic Shade, when shade increased and reduced sunlight by 0.18 h more each day, equivalent to a 13% reduction in incoming shortwave radiation and 16% in net radiation. Latent heat flux, the dominant energy flux at our site, varied temporally because of changes in incoming radiation, atmospheric demand, soil moisture and shade. Horizon shade controlled soil moisture at our site by prolonging snowmelt and reducing evapotranspiration in the late growing season, resulting in increased water storage capacity compared to other mountain wetlands. With the mounting risk of climate-change-driven severe spring flooding and late season droughts downstream of mountain headwaters, shaded subalpine wetlands provide important ecohydrological and mitigation services that are worthy of further study and mapping. This will help us better understand and protect mountain and prairie water resources.

A new flow for Canadian young hydrologists: Key scientific challenges addressed by research cultural shiftsCaroline Aubry-Wake1, Lauren D. Somers2,3, Hayley Alcock4, Aspen M. Anderson5, Amin Azarkhish6, Samuel Bansah7, Nicole M. Bell8, Kelly Biagi9, Mariana Castaneda-Gonzalez10, Olivier Champagne9, Anna Chesnokova10, Devin Coone6, Tasha-Leigh J. Gauthier11, Uttam Ghimire6, Nathan Glas6, Dylan M. Hrach11, Oi Yin Lai14, Pierrick Lamontagne-Halle3, Nicolas R. Leroux1, Laura Lyon3, Sohom Mandal12, Bouchra R. Nasri13, Nataša Popović11, Tracy. E. Rankin14, Kabir Rasouli15, Alexis Robinson16, Palash Sanyal17, Nadine J. Shatilla9, 18, Brandon Van Huizen11, Sophie Wilkinson9, Jessica Williamson11, Majid Zaremehrjardy191 Centre for Hydrology, University of Saskatchewan, Saskatoon, SK, Canada2 Civil and Environmental Engineering, Massachusetts Institute of Technology, MA, USA3 Department of Earth and Planetary Sciences, McGill University, Montreal QC4 Department of Natural Resource Science, McGill University, Montreal, QC, Canada5 Department of Earth Sciences, Simon Fraser University, Burnaby, BC, Canada6 School of Engineering, University of Guelph, Ontario, ON, Canada7 Department of Geological Sciences, University of Manitoba, Winnipeg, Canada8 Centre for Water Resources Studies, Department of Civil & Resource Engineering, Dalhousie University, Halifax, NS, Canada9 School of Geography and Earth Sciences, McMaster University, Hamilton, ON, Canada.10 Department of Construction Engineering, École de technologie supérieure, Montreal, QC, Canada11 Department of Geography & Environmental Management, University of Waterloo, Waterloo, ON, Canada12 Department of Geography and Environmental Studies, Ryerson University, Toronto, ON, Canada13 Department of Mathematics and Statistics, McGill University, Montréal, Qc, Canada14 Geography Department, McGill University, Montreal, QC, Canada15 Meteorological Service of Canada, Environment and Climate Change Canada, Dorval, QC, Canada16 Department of Geography and Planning, University of Toronto, Toronto, ON17 Global Institute for Water Security, University of Saskatchewan.18 Lorax Environmental Services Ltd, Vancouver, BC, Canada.19 Department of Earth and Atmospheric Sciences, University of Alberta, Edmonton, AB, Canada

Peatlands in the Western Boreal Plains act as important water sources in the landscape. Their persistence, despite potential evapotranspiration (PET) often exceeding annual precipitation, is attributed to various water storage mechanisms. One storage element that has been understudied is seasonal ground ice (SGI). This study characterized spring SGI conditions and explored its impacts on available energy, actual evapotranspiration, water table, and near surface soil moisture in a western boreal plains peatland. The majority of SGI melt took place over May 2017. Microtopography had limited impact on melt rates due to wet conditions. SGI melt released 139mm in ice water equivalent (IWE) within the top 30cm of the peat, and weak significant relationships with water table and surface moisture suggest that SGI could be important for maintaining vegetation transpiration during dry springs. Melting SGI decreased available energy causing small reductions in PET (<10mm over the melt period) and appeared to reduce actual evapotranspiration variability but not mean rates, likely due to slow melt rates. This suggests that melting SGI supplies water, allowing evapotranspiration to occur at near potential rates, but reduces the overall rate at which evapotranspiration could occur (PET). The role of SGI may help peatlands in headwater catchments act as a conveyor of water to downstream landscapes during the spring while acting as a supply of water for the peatland. Future work should investigate SGI influences on evapotranspiration under differing peatland types, wet and dry spring conditions, and if the spatial variability of SGI melt leads to spatial variability in evapotranspiration.

Melting seasonal ground ice (SGI) in western Boreal Plains (WBP) peatlands can reduce the available energy at the surface by reducing potential evapotranspiration (PET). PET often exceeds annual precipitation in the WBP. Including this effect in hydrological models may be important in assessing water deficits. However, SGI melt and the timing of ice-free conditions vary spatially, which suggests PET spatial variability could be influenced by SGI. Understanding this potential linkage can help improve site scale PET in peatland hydrological models. The objectives of this paper were (a) to quantify the effect of ice thickness and melt rate on peatland PET; (b) quantify the spatial variability of SGI thickness and melt rate across spatial scales; and (c) assess how/if spatial variability in SGI thickness/melt rate affects site scale PET. Results from the sensitivity analysis indicated that SGI thickness had a bigger impact on reducing PET compared with the melt rate. Two SGI thickness values were used that were observed on site: 0.32 m, which was measured in a more treed area, and 0.18 m, which was in a more open area. The 0.32 m had an average PET reduction of 14 mm (+/- 0.7), over the month of May, compared with 9 mm (+/- 1 mm) when there was 0.18 m of SGI, which are 13.7 and 8.8% reductions, respectively. SGI thickness and melt rate, both exhibited large- and small-scale spatial variability. At the large scale, spatial patterns in SGI thickness appeared to be influenced by extensive shading from the adjacent hillslopes. Small scale, SGI thickness may be a function of tree proximity and the snowpack. Finally, net radiation, rather than SGI, appeared to be the main driver behind PET spatial variability. This work enhances our conceptual understanding of the role of SGI in WBP peatlands. Future work can use the findings to better inform peatland hydrological models, allowing for better representation of peatlands in regional-scale models.

Warming conditions across Canada's subarctic and arctic regions are causing permafrost landforms to thaw, resulting in rapid land cover change, including conversion of peat plateaus to wetland and thermokarst. These changes have important implications for northern ecosystems, including shifting controls on carbon uptake and release functions, as well as altering evapotranspiration (ET) rates, which form feedbacks with climatic change. Four landforms (peat plateau, sedge lawn, channel fen, and a thermokarst shoreline collapse scar) in the Hudson Bay Lowlands, northern Manitoba, were instrumented for weekly chamber measurements of carbon dioxide (CO2) and water vapor flux over a summer season (May to September 2014). Relative to other landforms, thermokarst CO2 exchange was characterized by high respiration rates early in the season, which decreased and were offset later in the season by CO2 uptake driven by sedge productivity. For all landforms, ET peaked post-snowmelt during rapid active layer thaw, and decreased throughout the growing season, controlled primarily by atmospheric vapor deficits. This work shows distinct differences in CO2 exchange and ET between intact and thawing permafrost features. While representative of small-scale processes in a single study region over one growing season, the results presented in this study have important implications for our understanding of ecohydrological and biogeochemical functioning of subarctic landscapes under future climates.