Soil frost effects on soil water and runoff dynamics along a boreal forest transect: 1. Field investigations

To determine how soil frost changes flowpaths of runoff water along a hillslope, a transect consisting of four soil profiles directed towards a small stream in a mature forest stand was investigated at Svartberget, near Vindeln in northern Sweden. Soil temperature, unfrozen water content, groundwater level and snow depth were investigated along the transect, which started at the riparian peat, and extended 30 m upslope into mineral soils. The two, more organic‐rich profiles closest to the stream had higher water retention and wetter autumn conditions than the sandy mineral soils further upslope. The organic content of the soil influenced the variation in frost along the transect. The first winter (1995–96) had abnormally low snow precipitation, which gave a deep frost down to 40–80 cm, whereas the two following winters had frost depths of 5–20 cm. During winter 1995–96, the two organic profiles close to the stream had a shallower frost depth than the mineral soil profile higher upslope, but a considerably larger amount of frozen water. The fraction of water that did not freeze despite several minus degrees in the soil was 5–7 vol.% in the mineral soil and 10–15 vol.% in the organic soil. From the measurements there were no signs of perched water tables during any of the three snowmelt periods, which would have been strong evidence for changed water flowpaths due to soil frost. When shallow soil layers became saturated during snowmelt, especially in 1997 and 1998, it was because of rising groundwater levels. Several rain on frozen ground events during spring 1996 resulted in little runoff, since most of the rain either froze in the soil or filled up the soil water storage. Copyright © 2001 John Wiley & Sons, Ltd.     

Effects of frozen soil and snow cover on cold‐season soil water dynamics in Tokachi,Japan

Despite the potential impact of winter soil water movements in cold regions, relatively few field studies have investigated cold‐season hydrological processes that occur before spring‐onset of snowmelt infiltration. The contribution of soil water fluxes in winter to the annual water balance was evaluated over 5 years of field observations at an agricultural field in Tokachi, Hokkaido, Japan. In two of the winters, soil frost reached a maximum depth of 0·2 m (‘frozen’ winters), whereas soil frost was mostly absent during the remaining three winters (‘unfrozen’ winters). Significant infiltration of winter snowmelt water, to a depth exceeding 1·0 m, occurred during both frozen and unfrozen winters. Such infiltration ranged between 126 and 255 mm, representing 28–51% of total annual soil water fluxes. During frozen winters, a substantial quantity of water (ca 40 mm) was drawn from deeper layers into the 0–0·2 m topsoil layer when this froze. Under such conditions, the progression and regression of the freezing front, regulated by the thickness of snow cover, controlled the quantity of soil water flux below the frozen layer. During unfrozen winters, 13–62 mm of water infiltrated to a depth of 0·2 m, before the spring snowmelt. These results indicate the importance of correctly evaluating winter soil water movement in cold regions. Copyright © 2010 John Wiley & Sons, Ltd.     

Soil frost effects on soil water and runoff dynamics along a boreal transect: 2. Simulations

A physically based SVAT‐model was tested with soil and snow physical measurements, as well as runoff data from an 8600 m² catchment in northern Sweden in order to quantify the influence of soil frost on spring snowmelt runoff in a moderately sloped, boreal forest. The model was run as an array of connected profiles cascading to the brook. For three winter seasons (1995–98) it was able to predict the onset and total accumulation of the runoff with satisfactory accuracy. Surface runoff was identified as only a minor fraction of the total runoff occurring during short periods in connection with ice blocking of the water‐conducting pores. Little surface runoff, though, does not mean that soil frost is unimportant for spring runoff. Simulations without frost routines systematically underestimated the total accumulated runoff. The possibility of major frost effects appearing in response to specific combinations of weather conditions were also tested. Different scenarios of critical initial conditions for the winter, e.g. high water saturation and delayed snow accumulation leading to an increased frost penetration, were tested. These showed that under special circumstances there is potential for increased spring runoff due to soil frost. Copyright © 2001 John Wiley & Sons, Ltd.     

Integrating observations and models to determine the effect of seasonally frozen ground on hydrologic partitioning in alpine hillslopes in the Colorado Rocky Mountains,USA

This study integrated spatially distributed field observations and soil thermal models to constrain the impact of frozen ground on snowmelt partitioning and streamflow generation in an alpine catchment within the Niwot Ridge Long-Term Ecological Research site, Colorado, USA. The study area was comprised of two contrasting hillslopes with notable differences in topography, snow depth and plant community composition. Time-lapse electrical resistivity surveys and soil thermal models enabled extension of discrete soil moisture and temperature measurements to incorporate landscape variability at scales and depths not possible with point measurements alone. Specifically, heterogenous snowpack thickness (~0–4 m) and soil volumetric water content between hillslopes (~0.1–0.45) strongly influenced the depths of seasonal frost, and the antecedent soil moisture available to form pore ice prior to freezing. Variable frost depths and antecedent soil moisture conditions were expected to create a patchwork of differing snowmelt infiltration rates and flowpaths. However, spikes in soil temperature and volumetric water content, as well as decreases in subsurface electrical resistivity revealed snowmelt infiltration across both hillslopes that coincided with initial decreases in snow water equivalent and early increases in streamflow. Soil temperature, soil moisture and electrical resistivity data from both wet and dry hillslopes showed that initial increases in streamflow occurred prior to deep soil water flux. Temporal lags between snowmelt infiltration and deeper percolation suggested that the lateral movement of water through the unsaturated zone was an important driver of early streamflow generation. These findings provide the type of process-based information needed to bridge gaps in scale and populate physically based cryohydrologic models to investigate subsurface hydrology and biogeochemical transport in soils that freeze seasonally.     

Mechanism of the influence of soil freezing depth on winter runoff

The depth of soil freezing in river watersheds is a factor governing winter runoff formation. The freezing depth regulates the redistribution of stored soil moisture between thawed and frozen soil layers. The moisture stored in the thawed soil layer is spent for winter runoff, while that in the frozen zone forms snowmelt runoff. The depth of soil freezing has considerably decreased over the period of climate changes, resulting in an increase in winter runoff and greater losses of snowmelt runoff.     

Evaluating the hydrological and hydrochemical responses of a High Arctic catchment during an exceptionally warm summer

The Arctic has experienced substantial warming during the past century with models projecting continued warming accompanied by increases in summer precipitation for most regions. A key impact of increasing air surface temperatures is the deepening of the active layer, which is expected to alter hydrological processes and pathways. The aim of this study was to determine how one of the warmest and wettest summers in the past decade at a High Arctic watershed impacted water infiltration and storage in deeply thawed soil and solute concentrations in stream runoff during the thaw period. In June and July 2012 at the Cape Bounty Watershed Observatory, we combined active layer measurements with major ion concentrations and stable isotopes in surface waters to characterize the movement of different runoff sources: snowmelt, rainfall, and soil water. Results indicate that deep ground thaw enhanced the storage of infiltrated water following rainfall. Soil water from infiltrated rainfall flowed through the thawed transient layer and upper permafrost, which likely solubilized ions previously stored at depth. Subsequent rainfall events acted as a hydrological flushing mechanism, mobilizing solutes from the subsurface to the surface. This solute flushing substantially increased ion concentrations in stream runoff throughout mid to late July. Results further suggest the importance of rainfall and soil water as sources of runoff in a High Arctic catchment during mid to late summer as infiltrated snowmelt is drained from soil following baseflow. Although there was some evaporation of surface water, our study indicates that flushing from solute stores in the transient layer was the primary driver of increased ion concentrations in stream runoff and not evaporative concentration of surface water. With warmer and wetter summers projected for the Arctic, ion concentrations in runoff (especially in the late thaw season), will likely increase due to the deep storage and subsurface flow of infiltrated water and subsequent flushing of previously frozen solutes to the surface.     

Snow and frost: implications for spatiotemporal infiltration patterns – a review

Vast regions of the northern hemisphere are exposed to snowfall and seasonal frost. This has large effects on spatiotemporal distribution of infiltration and groundwater recharge processes as well as on the fate of pollutants. Therefore, snow and frost need to be central inherent elements of risk assessment and management schemes. However, snow and frost are often neglected or treated summarily or in a simplistic way by groundwater modellers. Snow deposition is uneven, and the snow is likely to sublimate, be redistributed and partly melt during the winter influencing the mass and spatial distribution of snow storage available for infiltration, the presence of ice layers within and under the snowpack and, therefore, also the spatial distribution of depths and permeability of the soil frost. In steep terrain, snowmelt may travel downhill tens of metres in hours along snow layers. The permeability of frozen soil is mainly influenced by soil type, its water and organic matter content, and the timing of the first snow in relation to the timing of sub‐zero temperatures. The aim with this paper is to review the literature on snow and frost processes, modelling approaches with the purpose to visualize and emphasize the need to include these processes when modelling, managing and predicting groundwater recharge for areas exposed to seasonal snow and frost. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling Canadian prairie wetland hydrology using a semi‐distributed streamflow model

A hydrological model (SLURP) that was designed for simulating hydrological processes taking place in large river basins was, with minimal modification, used successfully to simulate water level variations over a 28‐year period (1969–1996) for a 3‐ha prairie wetland in Saskatchewan. The model calculates a water balance based on precipitation, snowmelt, evaporation, surface runoff and subsurface flow on a daily time‐step. The model was first calibrated for two periods (1969–1973 for cropland and 1987–1990 for grassland), then it was applied to records outside the calibration periods. The model reproduced the wetland water level variations during a 28‐year period with good accuracy. The wetland water levels were most sensitive to the infiltration coefficient of surface soil under frozen conditions and to maximum soil moisture storage. The applicability of the model and the calibrated parameters to a smaller wetland, with an area of 0·24 ha, was examined. This simulation indicated that scale effects are important, probably largely in relation to snow redistribution by wind. Copyright © 2000 John Wiley & Sons, Ltd.     

Subsurface lateral flow generation in aspen and conifer‐dominated hillslopes of a first order catchment in northern Utah

Mountain headwater catchments in the semi‐arid Intermountain West are important sources of surface water because these high elevations receive more precipitation than neighboring lowlands. This study examined subsurface runoff in two hillslopes, one aspen dominated, the other conifer dominated, adjacent to a first order stream in snow‐driven northern Utah. Snow accumulation, soil moisture, trenchflow and streamflow were examined in hillslopes and their adjacent stream. Snow water equivalents (SWEs) were greater under aspen stands compared to conifer, the difference increasing with higher annual precipitation. Semi‐variograms of shallow spatial soil moisture patterns and transects of continuous soil moisture showed no increase in soil moisture downslope, suggesting the absence of subsurface flow in shallow (～12 cm) soil layers of either vegetation type. However, a clear threshold relationship between soil moisture and streamflow indicated hillslope–stream connectivity, deeper within the soil profile. Subsurface flow was detected at ～50 cm depth, which was sustained for longer in the conifer hillslope. Soil profiles under the two vegetation types varied, with deep aspen soils having greater water storage capacity than shallow rocky conifer soils. Though SWEs were less under the conifers, the soil profile had less water storage capacity and produced more subsurface lateral flow during the spring snowmelt. Copyright © 2010 John Wiley & Sons, Ltd.     

Controls on old and new water contributions to stream flow at some nested catchments in Vermont,USA

Factors controlling the partitioning of old and new water contributions to stream flow were investigated for three events in four catchments (three of which were nested) at Sleepers River Research Watershed in Danville, Vermont. In the 1993 snowmelt period, two‐component isotopic hydrograph separations showed that new water (meltwater) inputs to the stream ranged widely from 41 to 74%, and increased with catchment size (41 to 11 125 ha) (with one exception) and with open land cover (0–73%). Peak dissolved organic carbon concentrations and relative alkalinity dilution in stream water ranked in the same order among catchments as the new water fractions, suggesting that new water followed shallow flow paths. During the 1994 snowmelt, despite similar timing and magnitude of melt inputs, the new‐water contribution to stream flow ranged only from 30 to 36% in the four catchments. We conclude that the uncommonly high and variable new water fractions in streamwater during the 1993 melt were caused by direct runoff of meltwater over frozen ground, which was prevalent in open land areas during the 1993 winter. In a high‐intensity summer rainstorm in 1993, new water fractions were smaller relative to the 1993 snowmelt, ranging from 28 to 46%, but they ranked in the identical catchment order. Reconciliation of the contrasting patterns of new–old water partitioning in the three events appears to require an explanation that invokes multiple processes and effects, including:

• 1. topographically controlled increase in surface‐saturated area with increasing catchment size;
• 2. direct runoff over frozen ground;
• 3. low infiltration in agriculturally compacted soils;
• 4. differences in soil transmissivity, which may be more relevant under dry antecedent conditions.

These data highlight some of the difficulties faced by catchment hydrologists in formulating a theory of runoff generation at varying basin scales. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of snow and soil frost melting on the concentrations of suspended solids and phosphorus in two rural watersheds in Western Finland

Hydrological processes are known to have a considerable effect on nutrient transport from agricultural land to watercourses. In cold temperate regions peak discharges are caused not only by storm conditions but also by melting of snow and frost. The objective of this work was to investigate the effects of snow and frost melt on concentrations of phosphorus and suspended solids. The samples were taken using flowweighted automatic sampling techniques from two agricultural drainage basins. During the beginning of the snowmelt period the concentration of suspended solids was rather low by comparison with the total phosphorus concentration and the discharge. The different behaviour compared with the relationships found during storm conditions was probably caused by continuous extraction of the soil surface by low ionic concentration melt water, and to some extent by leaching from plant residues. The gradual decrease of total phosphorus concentration during the frost thawing period could be attributed to the gradually increasing proportion of the subsurface and ground water discharge in the runoff.     

Controls on evapotranspiration from jack pine forests in the Boreal Plains Ecozone

The exchanges of water, energy and carbon between the land surface and the atmosphere are tightly coupled, so that errors in simulating evapotranspiration lead to errors in simulating both the water and carbon balances. Areas with seasonally frozen soils present a particular challenge due to the snowmelt-dominated hydrology and the impact of soil freezing on the soil hydraulic properties and plant root water uptake. Land surface schemes that have been applied in high latitudes often have reported problems with simulating the snowpack and runoff. Models applied at the Boreal Ecosystem Research and Monitoring Sites in central Saskatchewan have consistently over-predicted evapotranspiration as compared with flux tower estimates. We assessed the performance of two Canadian land surface schemes (CLASS and CLASS-CTEM) for simulating point-scale evapotranspiration at an instrumented jack pine sandy upland site in the southern edge of the boreal forest in Saskatchewan, Canada. Consistent with past reported results, these models over-predicted evapotranspiration, as compared with flux tower observations, but only in the spring period. Looking systematically at soil properties and vegetation characteristics, we found that the dominant control on evapotranspiration within these models was the canopy conductance. However, the problem of excessive spring ET could not be solved satisfactorily by changing the soil or vegetation parameters. The model overestimation of spring ET coincided with the overestimation of spring soil liquid water content. Improved algorithms for the infiltration of snowmelt into frozen soils and plant-water uptake during the snowmelt and soil thaw periods may be key to addressing the biases in spring ET.     

Rain on snow induced urban floods in cold maritime climate: Risk,indicators and trends

Urban floods pose a societal and economical risk. This study evaluated the risk and hydro-meteorological conditions that cause pluvial flooding in coastal cities in a cold climate. Twenty years of insurance claims data and up to 97 years of meteorological data were analysed for Reykjavík, Iceland (64.15°N; <100 m above sea level). One third of the city's wastewater collection system is combined, and pipe grades vary from 0.5% to 10%. Results highlight semi-intensive rain (<7 mm/h; ≤3 year return period) in conjunction with snow and frozen ground as the main cause for urban flood risk in a climate which undergoes frequent snow and frost cycles (avg. 13 and 19 per season, respectively). Floods in winter were more common, more severe and affected a greater number of neighbourhoods than during summer. High runoff volumes together with debris remobilized with high winds challenged the capacity of wastewater systems regardless of their age or type (combined vs. separate). The two key determinants for the number of insurance claims were antecedent frost depth and total precipitation volume per event. Two pluvial regimes were particularly problematic: long duration (13–25 h), late peaking rain on snow (RoS), where snowmelt enhanced the runoff intensity, elongated and connected independent rainfall into a singular, more voluminous (20–76 mm) event; shorter duration (7–9 h), more intensive precipitation that evolved from snow to rain. Closely timed RoS and cooling were believed to trigger frost formation. A positive trend was detected in the average seasonal snow depth and volume of rain and snowmelt during RoS events. More emphasis, therefore, needs to be placed on designing and operating urban drainage infrastructure with regard to RoS co-acting with frozen ground. Furthermore, more detailed, routine monitoring of snow and soil conditions is important to predict RoS flood events.     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

Coupling snow accumulation and melt rate modules of monthly water balance models with the Jazim monthly water balance model

Several models for simulation of water balance processes in semi-arid mountainous basins were developed by coupling different modules of existing water balance models (WBM). Snow accumulation and snowmelt rate relationships extracted from the McCabe-Markstrom, Guo, Rao-Al Wagdany and WASMOD-M WBMs, originally developed for basins with humid climate, were coupled with the Jazim WBM, primarily developed for arid basins. Karaj Basin, central Iran, with snowy autumn–winter and dry summer periods, was selected to assess model performance. The model parameters were optimized using a genetic algorithm (GA). All coupled models performed better than the non-modified (original) WBMs in the study basin. The coupled Jazim–McCabe-Markstrom model provided the best performance in simulating low and high monthly flows. It estimated the snowmelt runoff values more accurately than other proposed coupled models because the linear relationships used in the snow module of the McCabe-Markstrom model are more compatible with snow variations in the Karaj Basin.     

Integrated Surface and Subsurface Hydrological Modeling with Snowmelt and Pore Water Freeze–Thaw

For the simulation of winter hydrological processes a gap in the availability of flow models existed: one either had the choice between (1) physically-based and fully-integrated, but computationally very intensive, or (2) simplified and compartamentalized, but computationally less expensive, simulators. To bridge this gap, we here present the integration of a computationally efficient representation of winter hydrological processes (snowfall, snow accumulation, snowmelt, pore water freeze–thaw) in a fully-integrated surface water-groundwater flow model. This allows the efficient simulation of catchment-scale hydrological processes in locations significantly influenced by winter processes. Snow accumulation and snowmelt are based on the degree-day method and pore water freeze–thaw is calculated with a vertical heat conduction approach. This representation of winter hydrological processes is integrated into the fully-coupled surface water-groundwater flow model HydroGeoSphere. A benchmark for pore water freeze–thaw as well as two illustrative examples are provided.     

Synergistic effect of vegetation and air temperature changes on soil water content in alpine frost meadow soil in the permafrost region of Qinghai‐Tibet

Seasonal changes over 2 years (2004–2006) in soil moisture content (θ_v) of frozen alpine frost meadow soils of the Qinghai‐Tibet plateau permafrost region under three different levels of vegetation cover were investigated. Vegetation cover and air temperature changes had significant effects (synergistic effect) on θ_v and its distribution in the soil profile. During periods of soil freezing or thawing, the less the vegetation cover, the quicker the temperature drop or rise of soil water, and the shorter the duration of the soil water freeze–thaw response in the active soil layer. Under 30% and 65% vegetation cover the amplitude of variation in θ_v during the freezing period was 20–26% greater than that under 93% cover, while during the thawing period, it was 1·5‐ to 40·5‐fold greater. The freezing temperature of the surface soil layer, ^fT_s, was 1·6 °C lower under 30% vegetation cover than under 93% vegetation cover. Changes in vegetation cover of the alpine frost meadow affected θ_v and its distribution, as well as the relationship between θ_v and soil temperature (T_s). As vegetation cover decreased, soil water circulation in the active layer increased, and the response to temperature of the water distribution across the soil profile was heightened. The quantity of transitional soil phase water at different depths significantly increased as vegetation cover decreased. The influence of vegetation cover and soil temperature distribution led to a relatively dry soil layer in the middle of the profile (0·70–0·80 m) under high vegetation cover. Alpine meadow θ_v and its pattern of distribution in the permafrost region were the result of the synergistic effect of air temperature and vegetation cover. Copyright © 2007 John Wiley & Sons, Ltd.     

Coupling of a simultaneous heat and water model with a distributed hydrological model and evaluation of the combined model in a cold region watershed

Snow and frozen soil prevail in cold regions worldwide, and the integration of these processes is crucial in hydrological models. In this study, a combined model was developed by fully coupling a simultaneous heat and water model with a geomorphologically based distributed hydrological model. The combined model simulates vertical and lateral water transfer as well as vertical heat fluxes and is capable of representing the effects of frozen soil and snowmelt on hydrological processes in cold regions. This model was evaluated by using in situ observations in the Binggou watershed, an experimental watershed for cold region hydrology of the Watershed Allied Telemetry Experimental Research Project. Results showed that the model was able to predict soil freezing and thawing, unfrozen soil water content, and snow depth reasonably well. The simulated hydrograph was in good agreement with the in situ observation. The Nash–Sutcliffe coefficient of daily discharge was 0.744 for the entire simulation period, 0.472 from April to June, and 0.711 from June to November. This model can improve our understanding of hydrological processes in cold regions and assess the impacts of global warming on hydrological cycles and water resources. Copyright © 2012 John Wiley & Sons, Ltd.