Macropore flow in frozen soils plays a critical role in partitioning snowmelt at the land surface and modulating snowmelt-driven hydrological processes. Previous descriptions of macropore flow processes in frozen soil do not explicitly represent the physics of water and heat transfer between macropores and the soil matrix, and there is a need to adapt recent conceptual and numerical models of unfrozen macropore flow to account for frozen ground. Macropores remain air filled under partially saturated conditions, allowing preferential flow and meltwater infiltration prior to ground thaw. Nonequilibrium gravity-driven flow can rapidly transport snowmelt to depths below the frost zone or, alternatively, infiltrated water may refreeze in macropores and restrict preferential flow. As with unfrozen soils, models of water movement in frozen soil that rely solely on diffuse flow concepts cannot adequately represent unsaturated macropore hydraulics. Dual-domain descriptions of unsaturated flow that explicitly define macropore hydraulic characteristics have been successful under unfrozen conditions but need refinement for frozen soils. In particular, because pore connectivity and hydraulic conductivity are influenced by ice content, modeling schemes specifying macropore-matrix interactions and refreezing of infiltrating water are critical. This review discusses the need for research on the interacting effects of macropore flow and soil freezethaw and the integration of these concepts into a framework of coupled heat and water transfer. As a result, it proposes a conceptual model of unsaturated flow in frozen macroporous soils that assumes two interacting domains (macropore and matrix) with distinct water and heat transfer regimes. Abbreviations: SFC, soil freezing characteristic; SMC, soil moisture characteristic. Citation: Mohammed, A.A., B.L. Kurylyk, E.E. Cey, and M. Hayashi. 2018. Snowmelt infiltration and macropore flow in frozen soils: Overview, knowledge gaps, and a conceptual framework. Vadose Zone J. 17:180084.

Snowmelt is a major source of groundwater recharge in cold regions. Throughout many landscapes snowmelt occurs when the ground is still frozen; thus frozen soil processes play an important role in snowmelt routing, and, by extension, the timing and magnitude of recharge. This study investigated the vadose zone dynamics governing snowmelt infiltration and groundwater recharge at three grassland sites in the Canadian Prairies over the winter and spring of 2017. The region is characterized by numerous topographic depressions where the ponding of snowmelt runoff results in focused infiltration and recharge. Water balance estimates showed infiltration was the dominant sink (35 %-85 %) of snowmelt under uplands (i.e. areas outside of depressions), even when the ground was frozen, with soil moisture responses indicating flow through the frozen layer. The refreezing of infiltrated meltwater during winter melt events enhanced runoff generation in subsequent melt events. At one site, time lags of up to 3 d between snow cover depletion on uplands and ponding in depressions demonstrated

Infiltration into frozen soil plays an important role in soil freeze-thaw and snowmeltdriven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near‐surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three‐dimensional model of a permafrost plateau that includes the supra‐permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug‐in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin‐up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady‐state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.

The infiltrability of frozen soils modulates the partitioning of snowmelt between infiltration and runoff in cold regions. Preferential flow in macropores may enhance infiltration, but flow dynamics in frozen soil are complicated by soil heat transfer processes. We developed a dual-permeability model that considers the interacting effects of freeze-thaw and preferential flow on infiltration and runoff generation in structured soils. This formulation was incorporated into the fully integrated groundwatersurface water model HydroGeoSphere, to represent water-ice phase change in macropores such that porewater freezing is governed by macropore-matrix heat exchange. Model performance was evaluated against laboratory experiments and synthetic test cases designed to examine the effects of preferential flow on snowmelt partitioning between infiltration, runoff, and drainage. Simulations were able to reproduce experimental observations of rapid infiltration and drainage behavior due to macropores very well, and approximated soil thaw to an acceptable degree. Simulation of measured data highlighted the importance of macropore hydraulic conductivity, as well as macropore-matrix heat and water transfer, on controlling preferential flow dynamics. Test cases replicated a range of snowmelt partitioning behavior commonly observed in frozen soils, including subsurface conditions that produce rapid infiltration and deeper drainage, the contrast between limited vs. unlimited infiltration responses to snowmelt, and the temporal evolution of runoff generation. This study demonstrates the important influence that water freezing along preferential flowpaths can have on infiltrability and runoff characteristics in frozen soils and provides a physically based description of this mechanism that links infiltration behavior to hydraulic and thermal properties of structured soils.