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.

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.