Surface and subsurface water contributions to streamflow from a mesoscale watershed in complex mountain terrain

An understanding of surface and subsurface water contributions to streamflow is essential for accurate predictions of water supply from mountain watersheds that often serve as water towers for downstream communities. As such, this study used the end‐member mixing analysis technique to investigate source water contributions and hydrologic flow paths of the 264 km² Boulder Creek Watershed, which drains the Colorado Front Range, USA. Four conservative hydrochemical tracers were used to describe this watershed as a 3 end‐member system, and tracer concentration reconstruction suggested that the application of end‐member mixing analysis was robust. On average from 2009 to 2011, snowmelt and rainwater from the subalpine zone and groundwater sampled from the upper montane zone contributed 54%, 22%, and 24% of the annual streamflow, respectively. These values demonstrate increased rainwater and decreased snow water contributions to streamflow relative to area‐weighted mean values derived from previous work at the headwater scale. Young water (2.3 ± 0.8 months) fractions of streamflow decreased from 18–22% in the alpine catchment to 8–10% in the lower elevation catchments and the watershed outlet with implications for subsurface storage and hydrological connectivity. These results contribute to a process‐based understanding of the seasonal source water composition of a mesoscale watershed that can be used to extrapolate headwater streamflow generation predictions to larger spatial scales.     

Water source dynamics of high Arctic river basins

Arctic river basins are amongst the most vulnerable to climate change. However, there is currently limited knowledge of the hydrological processes that govern flow dynamics in Arctic river basins. We address this research gap using natural hydrochemical and isotopic tracers to identify water sources that contributed to runoff in river basins spanning a gradient of glacierization (0–61%) in Svalbard during summer 2010 and 2011. Spatially distinct hydrological processes operating over diurnal, weekly and seasonal timescales were characterized by river hydrochemistry and isotopic composition. Two conceptual water sources (‘meltwater’ and ‘groundwater’) were identified and used as a basis for end‐member mixing analyses to assess seasonal and year‐to‐year variability in water source dynamics. In glacier‐fed rivers, meltwater dominated flows at all sites (typically >80%) with the highest contributions observed at the beginning of each study period in early July when snow cover was most extensive. Rivers in non‐glacierized basins were sourced initially from snowmelt but became increasingly dependent on groundwater inputs (up to 100% of total flow volume) by late summer. These hydrological changes were attributed to the depletion of snowpacks and enhanced soil water storage capacity as the active layer expanded throughout each melt season. These findings provide insight into the processes that underpin water source dynamics in Arctic river systems and potential future changes in Arctic hydrology that might be expected under a changing climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of climate on alpine stream chemistry and water sources

The resilience of alpine/subalpine watersheds may be viewed as the resistance of streamflow or stream chemistry to change under varying climatic conditions, which is governed by the relative size (volume) and transit time of surface and subsurface water sources. Here, we use end‐member mixing analysis in Andrews Creek, an alpine stream in Rocky Mountain National Park, Colorado, from water year 1994 to 2015, to explore how the partitioning of water sources and associated hydrologic resilience change in response to climate. Our results indicate that four water sources are significant contributors to Andrews Creek, including snow, rain, soil water, and talus groundwater. Seasonal patterns in source‐water contributions reflected the seasonal hydrologic cycle, which is driven by the accumulation and melting of seasonal snowpack. Flushing of soil water had a large effect on stream chemistry during spring snowmelt, despite making only a small contribution to streamflow volume. Snow had a large influence on stream chemistry as well, contributing large amounts of water with low concentrations of weathering products. Interannual patterns in end‐member contributions reflected responses to drought and wet periods. Moderate and significant correlations exist between annual end‐member contributions and regional‐scale climate indices (the Palmer Drought Severity Index, the Palmer Hydrologic Drought Index, and the Modified Palmer Drought Severity Index). From water year 1994 to 2015, the percent contribution from the talus‐groundwater end member to Andrews Creek increased an average of 0.5% per year (p < 0.0001), whereas the percent contributions from snow plus rain decreased by a similar amount (p = 0.001). Our results show how water and solute sources in alpine environments shift in response to climate variability and highlight the role of talus groundwater and soil water in providing hydrologic resilience to the system.     

Quantifying aggregation and change in runoff source in accordance with catchment area increase in a forested headwater catchment

There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km² forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km². The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km², becoming almost constant above 0.1 km². The specific discharge of groundwater showed large variability below 1 km² and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km², it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.     

A comparison of fibre‐optic distributed temperature sensing to traditional methods of evaluating groundwater inflow to streams

There are several methods for determining the spatial distribution and magnitude of groundwater inputs to streams. We compared the results of conventional methods [dye dilution gauging, acoustic Doppler velocimeter (ADV) differential gauging, and geochemical end‐member mixing] to distributed temperature sensing (DTS) using a fibre‐optic cable installed along 900 m of Ninemile Creek in Syracuse, New York, USA, during low‐flow conditions (discharge of 1·4 m³ s^?1). With the exception of differential gauging, all methods identified a focused, contaminated groundwater inflow and produced similar groundwater discharge estimates for that point, with a mean of 66·8 l s^?1 between all methods although the precision of these estimates varied. ADV discharge measurement accuracy was reduced by non‐ideal conditions and failed to identify, much less quantify, the modest groundwater input, which was only 5% of total stream flow. These results indicate ambient tracers, such as heat and geochemical mixing, can yield spatially and quantitatively refined estimates of relatively modest groundwater inflow even in large rivers. DTS heat tracing, in particular, provided the finest spatial characterization of groundwater inflow, and may be more universally applicable than geochemical methods, for which a distinct and consistent groundwater end member may be more difficult to identify. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrologic functioning of the deep critical zone and contributions to streamflow in a high‐elevation catchment: Testing of multiple conceptual models

High‐elevation mountain catchments are often subject to large climatic and topographic gradients. Therefore, high‐density hydrogeochemical observations are needed to understand water sources to streamflow and the temporal and spatial behaviour of flow paths. These sources and flow paths vary seasonally, which dictates short‐term storage and the flux of water in the critical zone (CZ) and affect long‐term CZ evolution. This study utilizes multiyear observations of chemical compositions and water residence times from the Santa Catalina Mountains Critical Zone Observatory, Tucson, Arizona to develop and evaluate competing conceptual models of seasonal streamflow generation. These models were tested using endmember mixing analysis, baseflow recession analysis, and tritium model “ages” of various catchment water sources. A conceptual model involving four endmembers (precipitation, soil water, shallow, and deep groundwater) provided the best match to observations. On average, precipitation contributes 39–69% (55 ± 16%), soil water contributes 25–56% (41 ± 16%), shallow groundwater contributes 1–5% (3 ± 2%), and deep groundwater contributes ~0–3% (1 ± 1%) towards annual streamflow. The mixing space comprised two principal planes formed by (a) precipitation‐soil water‐deep groundwater (dry and summer monsoon season samples) and (b) precipitation‐soil water‐shallow groundwater (winter season samples). Groundwater contribution was most important during the wet winter season. During periods of high dynamic groundwater storage and increased hydrologic connectivity (i.e., spring snowmelt), stream water was more geochemically heterogeneous, that is, geochemical heterogeneity of stream water is storage‐dependent. Endmember mixing analysis and ³H model age results indicate that only 1.4 ± 0.3% of the long‐term annual precipitation becomes deep CZ groundwater flux that influences long‐term deep CZ development through both intercatchment and intracatchment deep groundwater flows.     

Spatio‐temporal tracer variability in the glacier melt end‐member — How does it affect hydrograph separation results?

Geochemical and isotopic tracers were often used in mixing models to estimate glacier melt contributions to streamflow, whereas the spatio‐temporal variability in the glacier melt tracer signature and its influence on tracer‐based hydrograph separation results received less attention. We present novel tracer data from a high‐elevation catchment (17 km², glacierized area: 34%) in the Oetztal Alps (Austria) and investigated the spatial, as well as the subdaily to monthly tracer variability of supraglacial meltwater and the temporal tracer variability of winter baseflow to infer groundwater dynamics. The streamflow tracer variability during winter baseflow conditions was small, and the glacier melt tracer variation was higher, especially at the end of the ablation period. We applied a three‐component mixing model with electrical conductivity and oxygen‐18. Hydrograph separation (groundwater, glacier melt, and rain) was performed for 6 single glacier melt‐induced days (i.e., 6 events) during the ablation period 2016 (July to September). Median fractions (±uncertainty) of groundwater, glacier melt, and rain for the events were estimated at 49±2%, 35±11%, and 16±11%, respectively. Minimum and maximum glacier melt fractions at the subdaily scale ranged between 2±5% and 76±11%, respectively. A sensitivity analysis showed that the intraseasonal glacier melt tracer variability had a marked effect on the estimated glacier melt contribution during events with large glacier melt fractions of streamflow. Intra‐daily and spatial variation of the glacier melt tracer signature played a negligible role in applying the mixing model. The results of this study (a) show the necessity to apply a multiple sampling approach in order to characterize the glacier melt end‐member and (b) reveal the importance of groundwater and rainfall–runoff dynamics in catchments with a glacial flow regime.     

Understanding hydrological processes in glacierized catchments: Evidence and implications of highly variable isotopic and electrical conductivity data

The spatial and temporal characterization of geochemical tracers over Alpine glacierized catchments is particularly difficult, but fundamental to quantify groundwater, glacier melt, and rain water contribution to stream runoff. In this study, we analysed the spatial and temporal variability of δ²H and electrical conductivity (EC) in various water sources during three ablation seasons in an 8.4‐km² glacierized catchment in the Italian Alps, in relation to snow cover and hydro‐meteorological conditions. Variations in the daily streamflow range due to melt‐induced runoff events were controlled by maximum daily air temperature and snow covered area in the catchment. Maximum daily streamflow decreased with increasing snow cover, and a threshold relation was found between maximum daily temperature and daily streamflow range. During melt‐induced runoff events, stream water EC decreased due to the contribution of glacier melt water to stream runoff. In this catchment, EC could be used to distinguish the contribution of subglacial flow (identified as an end member, enriched in EC) from glacier melt water to stream runoff, whereas spring water in the study area could not be considered as an end member. The isotopic composition of snow, glacier ice, and melt water was not significantly correlated with the sampling point elevation, and the spatial variability was more likely affected by postdepositional processes. The high spatial and temporal variability in the tracer signature of the end members (subglacial flow, rain water, glacier melt water, and residual winter snow), together with small daily variability in stream water δ²H dynamics, are problematic for the quantification of the contribution of the identified end members to stream runoff, and call for further research, possibly integrated with other natural or artificial tracers.     

Controls of streamflow generation in small catchments across the snow–rain transition in the Southern Sierra Nevada,California

Processes controlling streamflow generation were determined using geochemical tracers for water years 2004–2007 at eight headwater catchments at the Kings River Experimental Watersheds in southern Sierra Nevada. Four catchments are snow‐dominated, and four receive a mix of rain and snow. Results of diagnostic tools of mixing models indicate that Ca²⁺, Mg²⁺, K⁺ and Cl^? behaved conservatively in the streamflow at all catchments, reflecting mixing of three endmembers. Using endmember mixing analysis, the endmembers were determined to be snowmelt runoff (including rain on snow), subsurface flow and fall storm runoff. In seven of the eight catchments, streamflow was dominated by subsurface flow, with an average relative contribution (% of streamflow discharge) greater than 60%. Snowmelt runoff contributed less than 40%, and fall storm runoff less than 7% on average. Streamflow peaked 2–4 weeks earlier at mixed rain–snow than snow‐dominated catchments, but relative endmember contributions were not significantly different between the two groups of catchments. Both soil water in the unsaturated zone and regional groundwater were not significant contributors to streamflow. The contributions of snowmelt runoff and subsurface flow, when expressed as discharge, were linearly correlated with streamflow discharge (R² of 0.85–0.99). These results suggest that subsurface flow is generated from the soil–bedrock interface through preferential pathways and is not very sensitive to snow–rain proportions. Thus, a declining of the snow–rain ratio under a warming climate should not systematically affect the processes controlling the streamflow generation at these catchments. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of groundwater and topography on stream drying in semi-arid headwater streams

Non-perennial streams comprise over half of the global stream network and impact downstream water quality. Although aridity is a primary driver of stream drying globally, surface flow permanence varies spatially and temporally within many headwater streams, suggesting that these complex drying patterns may be driven by topographic and subsurface factors. Indeed, these factors affect shallow groundwater flows in perennial systems, but there has been only limited characterisation of shallow groundwater residence times and groundwater contributions to intermittent streams. Here, we asked how groundwater residence times, shallow groundwater contributions to streamflow, and topography interact to control stream drying in headwater streams. We evaluated this overarching question in eight semi-arid headwater catchments based on surface flow observations during the low-flow period, coupled with tracer-based groundwater residence times. For one headwater catchment, we analysed stream drying during the seasonal flow recession and rewetting period using a sensor network that was interspersed between groundwater monitoring locations, and linked drying patterns to groundwater inputs and topography. We found a poor relationship between groundwater residence times and flowing network extent (R² < 0.24). Although groundwater residence times indicated that old groundwater was present in all headwater streams, surface drying also occurred in each of them, suggesting old, deep flowpaths are insufficient to sustain surface flows. Indeed, the timing of stream drying at any given point typically coincided with a decrease in the contribution from near-surface sources and an increased relative contribution of groundwater to streamflow at that location, whereas the spatial pattern of drying within the stream network typically correlated with locations where groundwater inputs were most seasonally variable. Topographic metrics only explained ~30% of the variability in seasonal flow permanence, and surprisingly, we found no correlation with seasonal drying and down-valley subsurface storage area. Because we found complex spatial patterns, future studies should pair dense spatial observations of subsurface properties, such as hydraulic conductivity and transmissivity, to observations of seasonal flow permanence.     

Spatial and temporal patterns of snow accumulation and aerial ablation across the McMurdo Dry Valleys,Antarctica

Snow in the McMurdo Dry Valleys is a potential source of moisture for subnivian soils in a cold desert ecosystem. In a water‐limited environment, enhanced soil moisture is expected to provide more favourable conditions for subnivian soil communities. In addition, snow cover insulates the underlying soil from air temperature extremes. Quantifying the spatial and temporal patterns of seasonal snow accumulation and ablation is necessary to understand these dynamics. Repeat high‐resolution imagery acquired for the 2009–2010 austral summer was used to map the seasonal distribution of snow across Taylor and Wright valleys, Southern Victorialand, Antarctica. An edge detection algorithm was used to perform an object‐based classification of snow‐covered area. Coupled with topographic parameters obtained from a 30‐m digital elevation model, unique distribution patterns were characterized for five regions within the neighbouring valleys. Time lapses of snow distribution in each region provide insight into spatially variable aerial ablation rates (change in area of landscape covered by snow) across the region. A strong coastal to interior gradient of decreasing snow‐covered area was evident for both Taylor and Wright valleys. The surrounding regions of Lake Fryxell, Lake Hoare, Lake Bonney, Lake Brownworth, and Lake Vanda exhibited losses of snow‐covered area of 9.61 km² (?93%), 1.63 km² (?72%), 1.07 km² (?97%), 2.60 km² (?82%), and 0.25 km² (?96%), respectively, as measured from peak accumulation in October to mid‐January. Differences in aerial ablation rates within and across local regions suggest that both topographic variation and regional microclimates influence the ablation of seasonal snow cover. Copyright © 2012 John Wiley & Sons, Ltd.     

A hydrochemical approach to quantify the role of return flow in a surface flow‐dominated catchment

Stormflow generation in headwater catchments dominated by subsurface flow has been studied extensively, yet catchments dominated by surface flow have received less attention. We addressed this by testing whether stormflow chemistry is controlled by either (a) the event‐water signature of overland flow, or (b) the pre‐event water signature of return flow. We used a high‐resolution hydrochemical data set of stormflow and end‐members of multiple storms in an end‐member mixing analysis to determine the number of end‐members needed to explain stormflow, characterize and identify potential end‐members, calculate their contributions to stormflow, and develop a conceptual model of stormflow. The arrangement and relative positioning of end‐members in stormflow mixing space suggest that saturation excess overland flow (26–48%) and return flow from two different subsurface storage pools (17–53%) are both similarly important for stormflow. These results suggest that pipes and fractures are important flow paths to rapidly release stored water and highlight the value of within‐event resolution hydrochemical data to assess the full range and dynamics of flow paths.     

Controls on the chemistry of runoff from an upland peat catchment

This study uses 2 years of data from a detailed weekly water sampling programme in a 11·4 km² upland peat catchment in the Northern Pennines, UK. The sampling comprised precipitation, soil‐water samples and a number of streams, including the basin outlet. Samples were analysed for: pH, conductivity, alkalinity, Na, K, Ca, Mg, Fe, Al, Total N, SO₄, Cl and colour. Principal component analysis (PCA) was used to identify end‐members and compositional trends in order to identify controls on the development of water composition. The study showed that the direct use of PCA had several advantages over the use of end‐member mixing analysis (EMMA) as it combines an analysis of mixing and evolving waters without the assumption of having to know the compositional sources of the water. In its application to an upland peat catchment, the study supports the view that shallow throughflow at the catotelm/acrotelm boundary is responsible for storm runoff generation and shows that baseflow is controlled by cation exchange in the catotelm and mixing with a base‐rich groundwater. Copyright © 2003 John Wiley & Sons, Ltd.     

A tracer‐based assessment of hydrological pathways at different spatial scales in a mesoscale Scottish catchment

Geochemically based hydrograph separation techniques were used in a preliminary assessment to infer how runoff processes change with landscape characteristics and spatial scale (1–233 km²) within a mesoscale catchment in upland Scotland. A two‐component end‐member mixing analysis (EMMA) used Gran alkalinity as an assumed conservative tracer. Analysis indicated that, at all scales investigated, acidic overland flow and shallow subsurface storm flows from the peaty soils covering the catchment headwaters dominated storm runoff generation. The estimated groundwater contribution to annual runoff varied from 30% in the smallest (ca 1 km²) peat‐dominated headwater catchment with limited groundwater storage, to >60% in larger catchments (>30 km²) with greater coverage of more freely draining soils and more extensive aquifers in alluvium and other drift. This simple approach offers a useful, integrated conceptualization of the hydrological functioning in a mesoscale catchment, which can be tested and further refined by focused modelling and process‐based research. However, even as it stands, the simple conceptualization of system behaviour will have significant utility as a tool for communicating hydrological issues in a range of planning and management decisions. Copyright © 2002 John Wiley & Sons, Ltd.     

A comparison of high‐resolution specific conductance‐based end‐member mixing analysis and a graphical method for baseflow separation of four streams in hydrologically challenging agricultural watersheds

Quantifying the relative contributions of different sources of water to a stream hydrograph is important for understanding the hydrology and water quality dynamics of a given watershed. To compare the performance of two methods of hydrograph separation, a graphical program [baseflow index (BFI)] and an end‐member mixing analysis that used high‐resolution specific conductance measurements (SC‐EMMA) were used to estimate daily and average long‐term slowflow additions of water to four small, primarily agricultural streams with different dominant sources of water (natural groundwater, overland flow, subsurface drain outflow, and groundwater from irrigation). Because the result of hydrograph separation by SC‐EMMA is strongly related to the choice of slowflow and fastflow end‐member values, a sensitivity analysis was conducted based on the various approaches reported in the literature to inform the selection of end‐members. There were substantial discrepancies among the BFI and SC‐EMMA, and neither method produced reasonable results for all four streams. Streams that had a small difference in the SC of slowflow compared with fastflow or did not have a monotonic relationship between streamflow and stream SC posed a challenge to the SC‐EMMA method. The utility of the graphical BFI program was limited in the stream that had only gradual changes in streamflow. The results of this comparison suggest that the two methods may be quantifying different sources of water. Even though both methods are easy to apply, they should be applied with consideration of the streamflow and/or SC characteristics of a stream, especially where anthropogenic water sources (irrigation and subsurface drainage) are present. Copyright © 2014 John Wiley & Sons, Ltd.     

Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 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.     

Isotopic investigation of runoff generation in a glacierized catchment in northern Sweden

In this study, summer rainfall contributions to streamflow were quantified in the sub‐arctic, 30% glacierized Tarfala (21.7 km²) catchment in northern Sweden for two non‐consecutive summer sampling seasons (2004 and 2011). We used two‐component hydrograph separation along with isotope ratios (δ¹⁸O and δD) of rainwater and daily streamwater samplings to estimate relative fraction and uncertainties (because of laboratory instrumentation, temporal variability and spatial gradients) of source water contributions. We hypothesized that the glacier influence on how rainfall becomes runoff is temporally variable and largely dependent on a combination of the timing of decreasing snow cover on glaciers and the relative moisture storage condition within the catchment. The results indicate that the majority of storm runoff was dominated by pre‐event water. However, the average event water contribution during storm events differed slightly between both years with 11% reached in 2004 and 22% in 2011. Event water contributions to runoff generally increased over 2011 the sampling season in both the main stream of Tarfala catchment and in the two pro‐glacial streams that drain Storglaciären (the largest glacier in Tarfala catchment covering 2.9 km²). We credit both the inter‐annual and intra‐annual differences in event water contributions to large rainfall events late in the summer melt season, low glacier snow cover and elevated soil moisture due to large antecedent precipitation. Together amplification of these two mechanisms under a warming climate might influence the timing and magnitude of floods, the sediment budget and nutrient cycling in glacierized catchments. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of wildfire on source water contributions in Devil Creek,CA: evidence from end‐member mixing analysis

A geochemical and end‐member mixing analysis (EMMA) is undertaken in Devil Canyon catchment, located in southern California, to further understanding of watershed behaviour and source water contributions after an acute and extensive wildfire. Physical and chemical transformations in post‐fire watersheds are known to increase overland flow and decrease infiltration, mainly due to formation of a hydrophobic layer at, or near, the soil surface. However, less is known about subsurface flow response in burned watersheds. The current study incorporates EMMA to evaluate and quantify source water contributions before, and after, a catchment affected by wildfires in southern California during the fall of 2003. Pre‐ and post‐fire stream water data were available at several sampling sites within the catchment, allowing the identification of contributing water sources at varying spatial scales. Proposed end‐member observations (groundwater, overland flow, shallow subsurface flow) were also collected to constrain and develop the catchment mixing model. Post‐fire source water changes are more evident in the smaller and faster responding sub‐basin (interior sampling point). Early post‐fire storm events are dominated by overland flow with no significant soil water or groundwater flow contribution. Inter‐storm streamwater in this smaller basin shows an increase in groundwater and a decrease in soil water. In the larger, baseflow‐dominated system, source water components appear less affected by fire. A slight increase in lateral flow is observed with only a slight decrease in baseflow. Changes in the post‐fire flow regimes affect nutrient loading and chemical response of the basin. Relatively rapid recovery of the chaparral ecosystem is evidenced, with active re‐growth and evapotranspiration evidenced by the fourth post‐fire rainy season. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydrological response to timber harvest in northern Idaho: implications for channel scour and persistence of salmonids

The potential for forest harvest to increase snowmelt rates in maritime snow climates is well recognized. However, questions still exist about the magnitude of peak flow increases in basins larger than 10 km² and the geomorphic and biological consequences of these changes. In this study, we used observations from two nearly adjacent small basins (13 and 30 km²) in the Coeur d'Alene River basin, one with recent, relatively extensive, timber harvest, and the other with little disturbance in the last 50 years to explore changes in peak flows due to timber harvest and their potential effects on fish. Peak discharge was computed for a specific rain‐on‐snow event using a series of physical models that linked predicted values of snowmelt input to a runoff‐routing model. Predictions indicate that timber harvest caused a 25% increase in the peak flow of the modelled event and increased the frequency of events of this magnitude from a 9‐year recurrence interval to a 3·6‐year event. These changes in hydrologic regime, with larger discharges at shorter recurrence intervals, are predicted to increase the depth and frequency of streambed scour, causing up to 15% added mortality of bull trout (Salvelinus confluentus) embryos. Mortality from increased scour, although not catastrophic, may have contributed to the extirpation of this species from the Coeur d'Alene basin, given the widespread timber harvest that occurred in this region. Copyright © 2008 John Wiley & Sons, Ltd.