Using an inverse modelling approach with equifinality control to investigate the dominant controls on snowmelt nutrient export

There is great interest in modelling the export of nitrogen (N) and phosphorus (P) from agricultural fields because of ongoing challenges of eutrophication. However, the use of existing hydrochemistry models can be problematic in cold regions because models frequently employ incomplete or conceptually incorrect representations of the dominant cold regions hydrological processes and are overparameterized, often with insufficient data for validation. Here, a process‐based N model, WINTRA, which is coupled to a physically based cold regions hydrological model, was expanded to simulate P and account for overwinter soil nutrient biochemical cycling. An inverse modelling approach, using this model with consideration of parameter equifinality, was applied to an intensively monitored agricultural basin in Manitoba, Canada, to help identify the main climate, soil, and anthropogenic controls on nutrient export. Consistent with observations, the model results suggest that snow water equivalent, melt rate, snow cover depletion rate, and contributing area for run‐off generation determine the opportunity time and surface area for run‐off–soil interaction. These physical controls have not been addressed in existing models. Results also show that the time lag between the start of snowmelt and the arrival of peak nutrient concentration in run‐off increased with decreasing antecedent soil moisture content, highlighting potential implications of frozen soils on run‐off processes and hydrochemistry. The simulations showed TDP concentration peaks generally arriving earlier than NO₃ but also decreasing faster afterwards, which suggests a significant contribution of plant residue Total dissolved Phosphorus (TDP) to early snowmelt run‐off. Antecedent fall tillage and fertilizer application increased TDP concentrations in spring snowmelt run‐off but did not consistently affect NO₃ run‐off. In this case, the antecedent soil moisture content seemed to have had a dominant effect on overwinter soil N biogeochemical processes such as mineralization, which are often ignored in models. This work demonstrates both the need for better representation of cold regions processes in hydrochemical models and the model improvements that are possible if these are included.     

Hydrological responses to climatic variability in a cold agricultural region

Extended severe dry and wet periods are frequently observed in the northern continental climate of the Canadian Prairies. Prairie streamflow is mainly driven by spring snowmelt of the winter snowpack, whilst summer rainfall is an important control on evapotranspiration and thus seasonality affects the hydrological response to drought and wet periods in complex ways. A field‐tested physically based model was used to investigate the influences of climatic variability on hydrological processes in this region. The model was set up to resolve agricultural fields and to include key cold regions processes. It was parameterized from local and regional measurements without calibration and run for the South Tobacco Creek basin in southern Manitoba, Canada. The model was tested against snow depth and streamflow observations at multiple scales and performed well enough to explore the impacts of wet and dry periods on hydrological processes governing the basin scale hydrological response. Four hydro‐climatic patterns with distinctive climatic seasonality and runoff responses were identified from differing combinations of wet/dry winter and summer seasons. Water balance analyses of these patterns identified substantive multiyear subsurface soil moisture storage depletion during drought (2001–2005) and recharge during a subsequent wet period (2009–2011). The fractional percentage of heavy rainfall days was a useful metric to explain the contrasting runoff volumes between dry and wet summers. Finally, a comparison of modeling approaches highlights the importance of antecedent fall soil moisture, ice lens formation during the snowmelt period, and peak snow water equivalent in simulating snowmelt runoff.     

Observations of urban snow properties in Calgary,Canada

Urban winter hydrology has garnered very little attention owing to the general notion that high‐intensity rainfalls are the major flood‐generating events in urban areas. As a result, few efforts have been made to research urban snow and its melt characteristics. This study investigates the characteristics of urban snow that differentiate it from rural snow, and makes recommendations for incorporating these characteristics into an urban snowmelt model. A field study was conducted from the fall of 2001 to the spring of 2002 in the city of Calgary, Canada. Snow depths and densities, soil moisture, soil temperature, snow albedo, net radiation, snow evaporation, and surface temperature were measured at several locations throughout the winter period. The combination of urban snow removal practices and the physical elements that exist in urban areas were found to influence the energy balance of the snowpack profoundly. Shortwave radiation was found to be the main source of energy for urban snow; as a consequence, the albedo of urban snow is a very important factor in urban snowmelt modelling. General observations lead to the classification of snow as one of four types: snow piles, snow on road shoulders, snow on sidewalk edges, and snow in open areas. This resulted in the development of four separate functions for the changing snow albedo values. A study of the frozen ground conditions revealed that antecedent soil moisture conditions had very little impact on frozen ground, and thus frozen ground very nearly always acts as a near impervious area. Improved flood forecasting for urban catchments in cold regions can only be achieved with accurate modelling of urban winter runoff that involves the energy balance method, incorporating snow redistribution and urban snow‐cover characteristics, and using small time steps. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydrological and hydrochemical response of a small canadian shield catchment to late winter rain-on-snow events

A study was undertaken during the winter of 1990–1991 in a small (3.7 ha) Canadian Shield catchment to examine the hydrological and hydrochemical response during rain-on-snow events. The results are presented of two large (37.9 and 34.6 mm) rain-on-snow events occurring in early and late March 1991. Peak and total runoff and the groundwater response from the two events are significantly different. Hydrological data indicate that these differences can be attributed to a combination of meteorological (temperature) and physical conditions (antecedent snowpack ripeness, soil moisture and groundwater levels). An immature snowpack (low temperature and density) combined with low antecedent soil moisture conditions significantly reduced the magnitude of the net hydrological input and runoff from the catchment during the early March event, whereas a more mature snowpack and high antecedent soil moisture conditions led to a large runoff event during late March. During both rain-on-snow events a significant portion of the pre-event snowpack chemical load was lost. Based on the maximum snowpack chemical load measured before the events, the two large rain-on-snow events and a brief mid-March warm period during which there were two much smaller rain-on-snow events removed 78% of the hydrogen ion and 63% of the sulphate and nitrate load from the snowpack, while only reducing snowpack water equivalence by 7%. A two-component (rain and snowmelt) isotopic (δ¹⁸O SMOW %₀) separation of snowmelt lysimeter water during the two events indicated that snowmelt was an important (50 and 65%, respectively) water source available for infiltration and runoff at the snow-soil interface. Considering the high hydrogen ion loadings to the catchment during these two events (3.3 and 3.0 mequiv.m^?2, respectively) streamflow pH was not significantly reduced due to an increase in the discharge of well-buffered groundwater. A two-component isotopic hydrograph separation of peak stream discharge during the 2–3 March event indicated that 75% of the total flow was groundwater. In mid-latitude acid-sensitive catchments, winter rain-on-snow events are an important hydrological occurrence due to their ability to elute much of the chemical load (H⁺, SO₄, NO₃) from the snowpack before the onset of spring melt when the maximum annual hydrological input typically occurs.     

One‐year‐long runoff forecast by a single snowpack evaluation

A degree‐day‐based model is presented for a 1 year ahead runoff forecast, with 1 day time steps. The input information is a single snowpack evaluation collected at the beginning of the snowmelt season. The snow‐cover dynamics, the key information for long‐term snowmelt forecast, are described by the snow‐line dynamics, i.e. by the movements of the downhill snowpack limit. The snowmelt volume, estimated by the snow‐line dynamics, is the exogenous input of an autoregressive transformation model. The model is calibrated by a least‐squares procedure on the basis of observed daily runoff data and the corresponding measurements of the snowpack volume (one measurement per year). A real‐world case study on the Alto Tunuyan River basin (2380 km², Argentinean Andes) is presented. The 1 year ahead Alto Tunuyan River runoff patterns, computed for both calibration and validation periods, reveal high agreement with observed streamflows. Copyright © 2004 John Wiley & Sons, Ltd.     

Testing a spatially distributed tracer‐aided runoff model in a snow‐influenced catchment: Effects of multicriteria calibration on streamwater ages

Integrating stable isotope tracers into rainfall‐runoff models allows investigation of water partitioning and direct estimation of travel times and water ages. Tracer data have valuable information content that can be used to constrain models and, in integration with hydrometric observations, test the conceptualization of catchment processes in model structure and parameterization. There is great potential in using tracer‐aided modelling in snow‐influenced catchments to improve understanding of these catchments' dynamics and sensitivity to environmental change. We used the spatially distributed tracer‐aided rainfall‐runoff (STARR) model to simulate the interactions between water storage, flux, and isotope dynamics in a snow‐influenced, long‐term monitored catchment in Ontario, Canada. Multiple realizations of the model were achieved using a combination of single and multiple objectives as calibration targets. Although good simulations of hydrometric targets such as discharge and snow water equivalent could be achieved by local calibration alone, adequate capture of the stream isotope dynamics was predicated on the inclusion of isotope data in the calibration. Parameter sensitivity was highest, and most local, for single calibration targets. With multiple calibration targets, key sensitive parameters were still identifiable in snow and runoff generation routines. Water ages derived from flux tracking subroutines in the model indicated a catchment where runoff is dominated by younger waters, particularly during spring snowmelt. However, resulting water ages were most sensitive to the partitioning of runoff sources from soil and groundwater sources, which was most realistically achieved when isotopes were included in the calibration. Given the paucity of studies where hydrological models explicitly incorporate tracers in snow‐influenced regions, this study using STARR is an important contribution to satisfactorily simulating snowpack dynamics and runoff generation processes, while simultaneously capturing stable isotope variability in snow‐influenced catchments.     

Seasonal snowpack dynamics and runoff in a cool temperate forest: lysimeter experiment in Niigata,Japan

Seasonal snowpack dynamics are described through field measurements under contrasting canopy conditions for a mountainous catchment in the Japan Sea region. Microclimatic data, snow accumulation, albedo and lysimeter runoff are given through the complete winter season 2002–03 in (1) a mature cedar stand, (2) a larch stand, and (3) a regenerating cedar stand or opening. The accumulation and melt of seasonal snowpack strongly influences streamflow runoff during December to May, including winter baseflow, mid‐winter melt, rain on snow, and diurnal peaks driven by radiation melt in spring. Lysimeter runoff at all sites is characterized by constant ground melt of 0·8–1·0 mm day⁻¹. Rapid response to mid‐winter melt or rainfall shows that the snowpack remains in a ripe or near‐ripe condition throughout the snow‐cover season. Hourly and daily lysimeter discharge was greatest during rain on snow (e.g. 7 mm h⁻¹ and 53 mm day⁻¹ on 17 December) with the majority of runoff due to rainfall passing through the snowpack as opposed to snowmelt. For both rain‐on‐snow and radiation melt events lysimeter discharge was generally greatest at the open site, although there were exceptions such as during interception melt events. During radiation melt instantaneous discharge was up to 4·0 times greater in the opening compared with the mature cedar, and 48 h discharge was up to 2·5 times greater. Perhaps characteristic of maritime climates, forest interception melt is shown to be important in addition to sublimation in reducing snow accumulation beneath dense canopies. While sublimation represents a loss from the catchment water balance, interception melt percolates through the snowpack and contributes to soil moisture during the winter season. Strong differences in microclimate and snowpack albedo persisted between cedar, larch and open sites, and it is suggested further work is needed to account for this in hydrological simulation models. Copyright © 2005 John Wiley & Sons, Ltd.     

Snowmelt-generated runoff and soil erosion in Fife,Scotland

The hydrology and contrasting erosional responses of two snowmelt events on arable farmland in Fife, Scotland, are compared. Snowmelt-generated runoff in January 1993 caused widespread soil erosion across eastern Scotland. Gullying was exemplified by three sites in Fife, where thaw of a drifted snowpack was augmented by rainfall to produce a larger erosive response than meteorological data alone would have predicted. Up to 127 m³ of soil was lost from individual gullies in fields sown to winter cereals. In February 1996 snowfall of comparable depth again covered the field area, but a more uniform snowpack, slower thaw, greater crop cover and lower rainfall during the thaw phase combined to lessen the impact of erosion. These case studies demonstrate the complexity of the erosion/runoff relationship for rain on snow events, in which erosional severity depends not just on snow depth but on snow distribution, thaw rate and the amount and timing of rainfall during the thaw phase. © 1998 John Wiley & Sons, Ltd.     

Quantifying the role of the snowpack in generating water available for run‐off during rain‐on‐snow events from snow pillow records

Rain‐on‐snow events have generated major floods around the world, particularly in coastal, mountainous regions. Most previous studies focused on a limited number of major rain‐on‐snow events or were based primarily on model results, largely due to a lack of long‐term records from lysimeters or other instrumentation for quantifying event water balances. In this analysis, we used records from five automated snow pillow sites in south coastal British Columbia, Canada, to reconstruct event water balances for 286 rain‐on‐snow events over a 10‐year period. For large rain‐on‐snow events (event rainfall >40 mm), snowmelt enhanced the production of water available for run‐off (WAR) by approximately 25% over rainfall alone. For smaller events, a range of antecedent and meteorological factors influenced WAR generation, particularly the antecedent liquid water content of the snowpack. Most large events were associated with atmospheric rivers. Rainfall dominated WAR generation during autumn and winter events, whereas snowmelt dominated during spring and summer events. In the majority of events, the sensible heat of rain contributed less than 10% of the total energy consumed by snowmelt. This analysis illustrated the importance of understanding the amount of rainfall occurring at high elevations during rain‐on‐snow events in mountainous regions.     

Hydroclimate controls over seasonal sediment yield in two adjacent High Arctic watersheds

Interannual variations in seasonal sediment transfer in two High Arctic non‐glacial watersheds were evaluated through three summers of field observations (2003–2005). Total seasonal discharge, controlled by initial watershed snow water equivalence (SWE) was the most important factor in total seasonal suspended sediment transfer. Secondary factors included melt energy, snow distribution and sediment supply. The largest pre‐melt SWE of the three years studied (2004) generated the largest seasonal runoff and disproportionately greater suspended sediment yield than the other years. In contrast, 2003 and 2005 had similar SWE and total runoff, but reduced runoff intensity resulted in lower suspended sediment concentrations and lower total suspended sediment yield in 2005. Lower air temperatures at the beginning of the snowmelt period in 2003 prolonged the melt period and increased meltwater storage within the snowpack. Subsequently, peak discharge and instantaneous suspended sediment concentrations were more intense than in the otherwise warmer 2005 season. The results for this study will aid in model development for sediment yield estimation from cold regions and will contribute to the interpretation of paleoenvironmental records obtained from sedimentary deposits in lakes. Copyright © 2007 John Wiley & Sons, Ltd.     

Influence of melt‐freeze‐cycles on the radionuclide transport in homogeneous laboratory snowpack

Radionuclides released to the environment and deposited with or onto snow can be stored over long time periods if ambient temperature stays low, particularly in glaciated areas or high alpine sites. The radionuclides will be accumulated in the snowpack during the winter unless meltwater runoff at the snow base occurs. They will be released to surface waters within short time during snowmelt in spring. In two experiments under controlled melting conditions of snow in the laboratory, radionuclide migration and runoff during melt‐freeze‐cycles were examined. The distribution of Cs‐134 and Sr‐85 tracers in homogeneous snow columns and their fractionation and potential preferential elution in the first meltwater portions were determined. Transport was associated with the percolation of meltwater at ambient temperatures above 0 °C after the snowpack became ripe. Mean migration velocities in the pack were examined for both nuclides to about 0.5 cm hr^?1 after one diurnal melt‐freeze‐cycle at ambient temperatures of ?2 to 4 °C. Meltwater fluxes were calculated with a median of 1.68 cm hr^?1. Highly contaminated portions of meltwater with concentration factors between 5 and 10 against initial bulk concentrations in the snowpack were released as ionic pulse with the first meltwater. Neither for caesium nor strontium preferential elution was observed. After recurrent simulated day‐night‐cycles (?2 to 4 °C), 80% of both radionuclides was released with the first 20% of snowmelt within 4 days. 50% of Cs‐134 and Sr‐85 were already set free after 24 hr. Snowmelt contained highest specific activities when the melt rate was lowest during the freeze‐cycles due to concentration processes in remaining liquids, enhanced by the melt‐freeze‐cycling. This implies for natural snowpack after significant radionuclide releases, that long‐time accumulation of radionuclides in the snow during frost periods, followed by an onset of steady meltwater runoff at low melt rates, will cause the most pronounced removal of the contaminants from the snow cover. This scenario represents the worst case of impact on water quality and radiation exposure in aquatic environments.     

Spatio‐temporal controls of snowmelt and runoff generation during rain‐on‐snow events in a mid‐latitude mountain catchment

A network of 30 standalone snow monitoring stations was used to investigate the snow cover distribution, snowmelt dynamics, and runoff generation during two rain‐on‐snow (ROS) events in a 40 km² montane catchment in the Black Forest region of southwestern Germany. A multiple linear regression analysis using elevation, aspect, and land cover as predictors for the snow water equivalent (SWE) distribution within the catchment was applied on an hourly basis for two significant ROS flood events that occurred in December 2012. The available snowmelt water, liquid precipitation, as well as the total retention storage of the snow cover were considered in order to estimate the amount of water potentially available for the runoff generation. The study provides a spatially and temporally distributed picture of how the two observed ROS floods developed in the catchment. It became evident that the retention capacity of the snow cover is a crucial mechanism during ROS. It took several hours before water was released from the snowpack during the first ROS event, while retention storage was exceeded within 1 h from the start of the second event. Elevation was the most important terrain feature. South‐facing terrain contributed more water for runoff than north‐facing slopes, and only slightly more runoff was generated at open compared to forested areas. The results highlight the importance of snowmelt together with liquid precipitation for the generation of flood runoff during ROS and the large temporal and spatial variability of the relevant processes. Copyright © 2015 John Wiley & Sons, Ltd.     

A temperature‐precipitation‐based model of thirty‐year mean snowpack accumulation and melt in Oregon,USA

High‐resolution, spatially extensive climate grids can be useful in regional hydrologic applications. However, in regions where precipitation is dominated by snow, snowmelt models are often used to account for timing and magnitude of water delivery. We developed an empirical, nonlinear model to estimate 30‐year means of monthly snowpack and snowmelt throughout Oregon. Precipitation and temperature for the period 1971–2000, derived from 400‐m resolution PRISM data, and potential evapotranspiration (estimated from temperature and day length) drive the model. The model was calibrated using mean monthly data from 45 SNOTEL sites and accurately estimated snowpack at 25 validation sites: R² = 0·76, Nash‐Sutcliffe Efficiency (NSE) = 0·80. Calibrating it with data from all 70 SNOTEL sites gave somewhat better results (R² = 0·84, NSE = 0·85). We separately applied the model to SNOTEL stations located < 200 and ≥ 200 km from the Oregon coast, since they have different climatic conditions. The model performed equally well for both areas. We used the model to modify moisture surplus (precipitation minus potential evapotranspiration) to account for snowpack accumulation and snowmelt. The resulting values accurately reflect the shape and magnitude of runoff at a snow‐dominated basin, with low winter values and a June peak. Our findings suggest that the model is robust with respect to different climatic conditions, and that it can be used to estimate potential runoff in snow‐dominated basins. The model may allow high‐resolution, regional hydrologic comparisons to be made across basins that are differentially affected by snowpack, and may prove useful for investigating regional hydrologic response to climate change. Published in 2011 by John Wiley & Sons, Ltd.     

Incorporating accumulated temperature and algorithm of snow cover calculation into the snowmelt runoff model

An accurate simulation of snowmelt runoff is of much importance in arid alpine regions. Data availability is usually an obstacle to use energy‐based snowmelt models for the snowmelt runoff simulation, and temperature‐based snowmelt models are more appealing in these regions. The snow runoff model is very popular nowadays, especially in the data sparse regions, because only temperature, precipitation and snow cover data are required for inputs to the model. However, this model uses average temperature as index, which cannot reflect the snowmelt simulation in the high altitude band. In this study, the snow runoff model is modified on the basis of accumulated active temperature. Snow cover calculation algorithm is added and is no longer needed as input but output. This makes the model able to simulate long‐time runoff and long‐time snow cover variation in every band. An examination of the improved model in the Manas River basin showed that the model is effective. It can reproduce the behaviour of the hydrology and can reflect the actual snow cover fluctuation. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrologic flowpaths during snowmelt in forested headwater catchments under differing winter climatic and soil frost regimes

Changes in hydrologic flowpaths have important impacts on the timing, magnitude and hydrochemistry of run‐off during snowmelt in forested catchments, but how flowpaths are affected by variation in winter climate and the irregular presence of soil frost remains poorly understood. The depth and extent of soil frost may be expected to increase as snowpack decreases or develops later because of climate change. In this study, we used end‐member mixing analysis to determine daily contributions of snow, forest floor soil water and groundwater to stream run‐off during snowmelt under different soil frost regimes resulting from interannual and elevational variation at the Hubbard Brook Experimental Forest in New Hampshire, USA. We observed greater routing of run‐off through forest floor flowpaths during early snowmelt in 2011, when the snowpack was deep and soil frost was minimal, compared with the early snowmelt in 2012 under conditions of deep and extensive soil frost. The results indicate that widespread soil frost that penetrated the depth of the forest floor decreased the flow signal through the shallowest subsurface flowpaths, but did not reduce overall infiltration of melt waters, as the contribution from the snow‐precipitation end‐member was similar under both conditions. These results are consistent with development of granular soil frost which permits vertical infiltration of melt waters, but either reduces lateral flow in the forest floor or prevents the solute exchange that would produce the typical chemical signature of shallow subsurface flowpaths in streamwater. Copyright © 2016 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.     

Development and application of a spatially‐distributed Arctic hydrological and thermal process model (ARHYTHM)

A process‐based, spatially distributed hydrological model was developed to quantitatively simulate the energy and mass transfer processes and their interactions within arctic regions (arctic hydrological and thermal model, ARHYTHM). The model first determines the flow direction in each element, the channel drainage network and the drainage area based upon the digital elevation data. Then it simulates various physical processes: including snow ablation, subsurface flow, overland flow and channel flow routing, soil thawing and evapotranspiration. The kinematic wave method is used for conducting overland flow and channel flow routing. The subsurface flow is simulated using the Darcian approach. The energy balance scheme was the primary approach used in energy‐related process simulations (snowmelt and evapotranspiration), although there are options to model snowmelt by the degree‐day method and evapotranspiration by the Priestley–Taylor equation. This hydrological model simulates the dynamic interactions of each of these processes and can predict spatially distributed snowmelt, soil moisture and evapotranspiration over a watershed at each time step as well as discharge in any specified channel(s). The model was applied to Imnavait watershed (about 2·2 km²) and the Upper Kuparuk River basin (about 146 km²) in northern Alaska. Simulated results of spatially distributed soil moisture content, discharge at gauging stations, snowpack ablations curves and other results yield reasonable agreement, both spatially and temporally, with available data sets such as SAR imagery‐generated soil moisture data and field measurements of snowpack ablation, and discharge data at selected points. The initial timing of simulated discharge does not compare well with the measured data during snowmelt periods mainly because the effect of snow damming on runoff was not considered in the model. Results from the application of this model demonstrate that spatially distributed models have the potential for improving our understanding of hydrology for certain settings. Finally, a critical component that led to the performance of this modelling is the coupling of the mass and energy processes. Copyright © 2000 John Wiley & Sons, Ltd.     

Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

Transformations of runoff chemistry in the Arctic tundra,Northwest Territories,Canada

The transformation of snowmelt water chemical composition during melt, elution and runoff in an Arctic tundra basin is investigated. The chemistry of the water flowing along pathways from the surface of melting snow to the 95·5 ha basin outlet is related to relevant hydrological processes. In so doing, this paper offers physically based explanations for the transformation of major ion concentrations and loads of runoff water associated with snowmelt and rainfall along hydrological pathways to the stream outlet. Late‐lying snowdrifts were found to influence the ion chemistry in adjacent reaches of the stream channel greatly. As the initial pulse of ion‐rich melt water drained from the snowdrift and was conveyed through hillslope flowpaths, the concentrations of most ions increased, and the duration of the peak ionic pulse lengthened. Over the first 3 m of overland flow, the concentrations of all ions except for NO increased by one to two orders of magnitude, with the largest increase for K⁺, Ca²⁺ and Mg²⁺. This was roughly equivalent to the concentration increase that resulted from percolation of relatively dilute water through 0·25 m of unsaturated soil. The Na⁺ and Cl^? were the dominant ions in snowmelt water, whereas Ca²⁺ and Mg²⁺ dominated the hillslope runoff. On slopes below a large melting snowdrift, ion concentrations of melt water flowing in the saturated layer of the soil were very similar to the relatively dilute concentrations found in surface runoff. However, once the snowdrift ablated, ion concentrations of subsurface flow increased above parent melt‐water concentrations. Three seasonally characteristic hydrochemical regimes were identified in a stream reach adjacent to late‐lying snowdrifts. In the first two stages, the water chemistry in the stream channel strongly resembled the hillslope drainage water. In the third stage, in‐stream geochemical processes, including the weathering/ion exchange of Ca²⁺ and Mg²⁺, were the main control of streamwater chemistry. Copyright © 2006 John Wiley & Sons, Ltd.     

Seasonal dynamics in shallow freshwater pond‐peatland hydrochemical interactions in a subarctic permafrost environment

Terrestrial and aquatic ecological productivity are often nutrient limited in subarctic permafrost environments. High latitude regions are experiencing significant climatic change, including rapid warming and changing precipitation patterns, which may result in changes in nutrient dynamics within terrestrial and aquatic systems and hydrochemical transport between them. The objective of this research was to characterize changes in runoff quantity and quality within, and between peatlands and ponds throughout the snow‐free summer season. Two ponds and their catchments were monitored over the snow‐free season to measure changes in hydrologic storage, and to determine how water chemistry changed with the evolution of the frost table depth. Thresholds in hydrologic storage combined with frost table position (which inhibited infiltration and storage) produced nonlinear responses for runoff generation through highly conductive shallow peat layers while deeper, less conductive layers retarded flow. Greater inputs were required to exceed hydrologic storage (fill and spill) as a deepening frost table increased the hydrologically active portion of the soil, leading to seasonal variability in runoff pathways between peatlands and ponds. Runoff contributions to ponds were an integral component of the snow‐free water balance during the study period, contributing up to 60% of all snow‐free inputs. Groundwater chemistry (and pond chemistry following runoff events when ponds were connected with peatlands) reflected the different depths of peat and mineral soil accessed throughout the season. This work has improved scientific understanding of the combined controls of hydrologic inputs and ground frost on runoff and nutrient transport between peatlands and ponds, and sheds insight into how nutrient dynamics in cold regions may evolve under a changing climate.