The use of agent based modelling techniques in hydrology: determining the spatial and temporal origin of channel flow in semi‐arid catchments

Information on the spatial and temporal origin of runoff entering the channel during a storm event would be valuable in understanding the physical dynamics of catchment hydrology; this knowledge could be used to help design flood defences and diffuse pollution mitigation strategies. The majority of distributed hydrological models give information on the amount of flow leaving a catchment and the pattern of fluxes within the catchment. However, these models do not give any precise information on the origin of runoff within the catchment. The spatial and temporal distribution of runoff sources is particularly intricate in semi‐arid catchments, where there are complex interactions between runoff generation, transmission and re‐infiltration over short temporal scales. Agents are software components that are capable of moving through and responding to their local environment. In this application, the agents trace the path taken by water through the catchment. They have information on their local environment and on the basis of this information make decisions on where to move. Within a given model iteration, the agents are able to stay in the current cell, infiltrate into the soil or flow into a neighbouring cell. The information on the current state of the hydrological environment is provided by the environment generator. In this application, the Connectivity of Runoff Model (CRUM) has been used to generate the environment. CRUM is a physically based, distributed, dynamic hydrology model, which considers the hydrological processes relevant for a semi‐arid environment at the temporal scale of a single storm event. During the storm event, agents are introduced into the model across the catchment; they trace the flows of water and store information on the flow pathways. Therefore, this modelling approach is capable of giving a novel picture of the temporal and spatial dynamics of flow generation and transmission during a storm event. This is possible by extracting the pathways taken by the agents at different time slices during the storm. The agent based modelling approach has been applied to two small catchments in South East Spain. The modelling approach showed that the two catchments responded differently to the same rainfall event due to the differences in the runoff generation and overland flow connectivity between the two catchments. The model also showed that the time of travel to the nearest flow concentration is extremely important for determining the connectivity of a point in the landscape with the catchment outflow. Copyright © 2007 John Wiley & Sons, Ltd.     

Contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments

We examined the contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments by conducting detailed hydrochemical observations in five catchments that ranged from zero to second order. End‐member mixing analysis (EMMA) was performed to identify the geographical sources of stream water. Throughfall, hillslope groundwater, shallow bedrock groundwater, and deep bedrock groundwater were identified as end members. The contribution of each end member to storm runoff differed among the catchments because of the differing quantities of riparian groundwater, which was recharged by the bedrock groundwater prior to rainfall events. Among the five catchments, the contribution of throughfall was highest during both baseflow and storm flow in a zero‐order catchment with little contribution from the bedrock groundwater to the riparian reservoir. In zero‐order catchments with some contribution from bedrock groundwater, stream water was dominated by shallow bedrock groundwater during baseflow, but it was significantly influenced by hillslope groundwater during storms. In the first‐order catchment, stream water was dominated by shallow bedrock groundwater during storms as well as baseflow periods. In the second‐order catchment, deeper bedrock groundwater than that found in the zero‐order and first‐order catchments contributed to stream water in all periods, except during large storm events. These results suggest that bedrock groundwater influences the upscaling of storm‐runoff generation processes by affecting the linkages of geomorphic units such as hillslopes, riparian zones, and stream channels. Our results highlight the need for a three‐dimensional approach that considers bedrock groundwater flow when studying the upscaling of storm‐runoff generation processes. Copyright © 2014 John Wiley & Sons, Ltd.     

Stormflow generation in two headwater catchments in eastern Amazonia,Brazil

Throughout the tropics, and the Amazon region in particular, only a few experimental studies have identified the main hydrological pathways and response to storm events. This study identifies the hydrological response patterns and quantifies the main runoff generating processes for two headwater catchments in eastern Amazonia, an area of low relief. Over an 18 month study period, 245 and 55 rainfall–runoff events at the respective headwater catchments were analysed. The rainfall‐runoff regression lines for both catchments revealed a remarkably strong linear correlation between event rain total and runoff volume. The area contributing to stormflow was proven to be constant in extent at approximately 0·6% of the catchment and to coincide with the exact extent of the riparian wetland zone. The soils of the surrounding hillslopes were found to be highly permeable oxisols. Indications of secondary permeability due to the deep root system of the secondary vegetation were also found. Copyright © 2007 John Wiley & Sons, Ltd.     

Concentration–discharge relationships describe solute and sediment mobilization,reaction, and transport at event and longer timescales

Concentration–discharge (C‐Q) relationships reflect material sources, storage, reaction, proximity, and transport in catchments. Differences in hydrologic pathways and connectivity influence observed C‐Q patterns at the catchment outlet. We examined solute and sediment C‐Q relationships at event and interannual timescales in a small mid‐Atlantic (USA) catchment. We found systematic differences in the C‐Q behaviour of geogenic/exogenous solutes (e.g., calcium and nitrate), biologically associated solutes (e.g., dissolved organic carbon), and particulate materials (e.g., total suspended solids). Negative log(C)–log(Q) regression slopes, indicating dilution, were common for geogenic solutes whereas positive slopes, indicating concentration increase, were common for biologically associated solutes. Biologically associated solutes often exhibited counterclockwise hysteresis during events whereas geogenic solutes exhibited clockwise hysteresis. Across event and interannual timescales, solute C‐Q patterns are linked to the spatial distribution of hydrologic sources and the timing and sequence of hydro‐biogeochemical source contributions to the stream. Groundwater is the primary source of stormflow during the earliest and latest stages of events, whereas precipitation and soil water become increasingly connected to the stream near peakflow. This sequence and timing of flowpath connectivity results in dilution and clockwise hysteresis for geogenic/exogenous solutes and concentration increase and counterclockwise hysteresis for biologically associated solutes. Particulate materials demonstrated positive C‐Q slopes over the long‐term and clockwise hysteresis during individual events. Drivers of particulate and solute C‐Q relationships differ, based on longitudinal and lateral expansion of active channels and changing shear stresses with increasing flows. Although important distinctions exist between the drivers of solute and sediment C‐Q relationships, overall solute and sediment C‐Q patterns at event and interannual timescales reflect consistent catchment hydro‐biogeochemical processes.     

Runoff generation in zero-order precambrian shield catchments: The stormflow response of a heterogeneous landscape

Hydrometric and isotopic (oxygen-18) observations were used to delineate the runoff processes operating in several headwater catchments on the Precambrian Shield of Canada. The catchments comprise patches of conifer forest situated on thin soils among areas of lichen-covered granitic bedrock. Horton overland flow occurred from the lichen-bedrock areas in all precipitation events that exceeded 4–6 mm. Runoff from the forest stands occurred mainly as subsurface stormflow, but in some instances saturation overland flow was observed. The occurrence of saturation overland flow was controlled by the topography of the bedrock beneath the forest soils. The area contributing runoff and the pathway by which water was conveyed to the catchment outflow switched from the open lichen-bedrock areas producing overland flow on the rising limb of the storm hydrograph to the forest stands contributing subsurface stormflow on the recession limb of the hydrograph. The areal extent and position of the landscape units in the basin were important to the rate and magnitude of stormflow production. Runoff was generated from the catchments only during and immediately after snowmelt and/or rainfall events. The catchments were dry and/or frozen for about 70% of the year.     

Temporal patterns and controls on runoff magnitude and solution chemistry of urban catchments in the semiarid southwestern United States

Urban expansion and the scarcity of water supplies in arid and semiarid regions have increased the importance of urban runoff to localized water resources. However, urban catchment responses to precipitation are poorly understood in semiarid regions where intense rainfall often results in large runoff events during the short summer monsoon season. To evaluate how urban runoff quantity and quality respond to rainfall magnitude and timing, we collected stream stage data and runoff samples throughout the 2007 and 2008 summer monsoons from four ephemeral drainages in Tucson, Arizona. Antecedent rainfall explained 20% to 30% of discharge (mm) and runoff ratio in the least impervious (22%) catchment but was not statistically related to hydrologic responses at more impervious sites. Regression models indicated that rainfall depth, imperviousness and their combined effect control discharge and runoff ratios (p < 0.01, r² = 0.91 and 0.75, respectively). In contrast, runoff quality did not vary with imperviousness or catchment size. Rainfall depth and duration, time since antecedent rainfall and event and cumulative discharge controlled runoff hydrochemistry and resulted in five specific solute response patterns: (i) strong event and seasonal solute mobilization (solute flush), (ii) event chemostasis and strong seasonal flush, (iii) event chemostasis and weak seasonal flush, (iv) event and seasonal chemostasis and (v) late seasonal flush. Our results indicate that hydrologic responses of semiarid catchments are controlled by rainfall partitioning at the event scale, whereas wetting magnitude, frequency and timing alter solute stores readily available for transport and control temporal runoff quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of hillslope topography on runoff response in a small catchment in the Fudoji Experimental Watershed,central Japan

We investigated the role of different hillslope units with different topographic characteristics on runoff generation processes based on field observations at two types of hillslopes (0·1 ha): a valley‐head (a convergent hillslope) and a side slope (a planar hillslope), as well as at three small catchments having two types of slopes with different drainage areas ranging from 1·9 to 49·7 ha in the Tanakami Mountains, central Japan. We found that the contribution of the hillslope unit type to small catchment runoff varied with the magnitude of rainfall. When the total amount of rainfall for a single storm event was < 35 mm, runoff in the small catchment was predominantly generated from the side slope. As the amount of rainfall increased (>35 mm), the valley‐head also began to contribute to the catchment runoff, adding to runoff from the side slope. Although the direct runoff from the valley‐head was greater than that from the side slope, the contribution from the side slope was quantitatively greater than that from the valley‐head due to the proportionally larger area occupied by the side slope in the small catchment. The storm runoff responses of the small catchments reflected the change in the runoff components of each hillslope unit as the amount of rainfall increased and rainfall patterns changed. However, similar runoff responses were found for the small catchments with different areas. The similarity of the runoff responses is attributable to overlay effects of different hillslope units and the similar composition ratios of the valley‐head and side slope in the catchments. This study suggests that the relative roles of the valley‐head and side slope are important in runoff generation and solute transport as the catchment size increases from a hillslope/headwater to a small catchment. Copyright © 2011 John Wiley & Sons, Ltd.     

HYDROCHEMICAL DECIPHERING OF STREAMFLOW GENERATION IN SEMI-ARID EAST AFRICA

There is a critical lack of knowledge regarding the dynamics of streamflow generation in the semi-arid tropics, particularly in Africa. In this project runoff mechanisms in forested and non-forested degraded catchments in northern Tanzania were studied using combined hydrometrical and hydrochemical methods. Following the hydrochemical identification of several flowpaths contributing to runoff, hydrograph separation by an end-member model based on K and Ca was undertaken. Results from the forested catchment indicate that stormflow was dominated by event water (about 75%), via overland flow and throughflow. The proportion of pre-event water (groundwater) displaced into the stream by a suggested riparian groundwater ridge mechanism varied, depending on the rainfall characteristics. In the non-forested, degraded catchment, nearly all stormflow was event water, and groundwater discharge was unaffected by rainfall. It is suggested that macropore flow is pivotal to the transmission of rainfall to runoff via throughflow, particularly in semi-arid tropical areas.     

Sources and sinks of nutrients in a New Zealand hill pasture catchment I. Stormflow generation

The processes of stormflow generation were studied in a hill pasture catchment near Hamilton, New Zealand. Although rainfall was relatively evenly distributed throughout the year, stormflow was highly seasonal and over 65 per cent occurred during the winter. Three main processes contributing to stormflow were identified which could be related to soil type and physiographic position. On gleyed soils derived from rhyolitic colluvium, saturation overland flow was the dominant process. Hydrographs from ‘Whipkey’ throughflow troughs also indicated that there was a subsurface response (saturated wedge) from this soil type. On steeper convex slopes, more permeable soils were derived from weathered greywacke. The presence of ephemeral springs on the hillslopes and direct observation during storm events indicated that storm runoff was generated as return flow from this soil. It was noted that nitrate concentrations from subsurface sources were 5–10 times higher than surface runoff. This difference in concentration was utilized in a chemical mixing equation which partitioned stormflow sources. This was compared with the stormflow predicted from rain falling on to saturated areas. There was good agreement between the two models for winter-spring events with respect to the volumes of surface runoff predicted, however the saturated areas model underestimated total stormflow. The results of the study are briefly discussed in terms of the potential for water quality management.     

Spatial delineation of groundwater–surface water interactions through intensive in‐stream profiling

The dominance of ‘old’ pre‐event water in headwater storm runoff has been recorded in numerous upland catchment studies; however, the mechanisms by which this pre‐event water enters the stream channel are poorly understood. Understanding these processes is fundamental to determining the controls on surface water quality and associated impacts on stream ecology. Previous studies in the upland forested catchment of the Afon Hafren (River Severn) at Plynlimon, mid‐Wales, identified an active bedrock groundwater system that was discharging into the stream channel during storm response. Detailed analysis showed that these discharges were small and could not account for the majority of pre‐event storm water response identified at this site; pre‐event storm runoff had to be sourced predominantly from further upstream. An intensive stream survey was used to determine the spatial nature of groundwater–surface water (GW–SW) interactions in the Hafren Catchment. Detailed physico‐chemical in‐stream profiling identified a marked change in water quality indicating a significant discrete point of bedrock groundwater discharge upstream of the Hafren Transect study site. The in‐stream profiling showed the importance of high spatial resolution sampling as a key to understanding processes of GW–SW interaction and how quick and cost‐effective measurements of specific electrical conductance of stream waters could be used to highlight in‐stream heterogeneity. This approach is recommended for use in headwater catchments for initial characterisation of the stream channel in order to better locate instrumentation and to determine more effective targeted sampling protocols in upland catchment research. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatial and temporal variability of the hydrological response in a small Mediterranean research catchment (Vallcebre,Eastern Pyrenees)

This paper analyses the spatial and temporal variability of the hydrological response in a small Mediterranean catchment (Cal Rodó). The first part of the analysis focuses on the rainfall–runoff relationship at seasonal and monthly scale, using an 8‐year data set. Then, using storm‐flow volume and coefficient, the temporal variability of the rainfall–runoff relationship and its relationship with several hydrological variables are analysed at the event scale from hydrographs observed over a 3‐year period. Finally, the spatial non‐linearity of the hydrological response is examined by comparing the Cal Rodó hydrological response with the Can Vila sub‐catchment response at the event scale. Results show that, on a seasonal and monthly scale, there is no simple relationship between rainfall and runoff depths, and that evapotranspiration is a factor that introduced some non‐linearity in the rainfall–runoff relationship. The analysis of monthly values also reveals the existence of a threshold in the relationship between rainfall and runoff depths, denoting a more contrasted hydrological response than the one usually observed in humid catchments. At the event scale, the storm‐flow coefficient has a clear seasonal pattern with an alternance between a wet period, when the catchment is hydrologically responsive, and a dry summer period, when the catchment is much less reactive to any rainfall. The relationship between the storm‐flow coefficient and rainfall depth, rainfall maximum intensity and base‐flow shows that observed correlations are the same as those observed for humid conditions, even if correlation coefficients are notably lower. Comparison with the Can Vila sub‐catchment highlights the spatial heterogeneity of the rainfall‐runoff relationship at the small catchment scale. Although interpretation in terms of runoff processes remains delicate, heterogeneities between the two catchments seem to be related to changes in the ratio between infiltration excess and saturation processes in runoff formation. Copyright © 2007 John Wiley & Sons, Ltd.     

A case study in catchment hydrochemistry: Conflicting interpretations from hydrological and chemical observations

Soil water, stream water, groundwater and rain water were sampled through a storm event in a moorland catchment. Samples were analysed for major ions and deuterium. Chloride and deuterium are used as tracers to enable separation of the stream runoff hydrograph into three components: rain, soil and groundwater. The results indicate that rain water arrives in the stream quickly during the event and contributes a significant volume to the runoff peak. The chemical signal in the rain water is, however, significantly damped, apparently due to mixing with soil water held in the catchment before the event. This is further modified before reaching the stream, apparently through mixing with a deeper groundwater component. Interpretation of tracer, chemistry and hydrological data to present an integrated picture of catchment hydrochemical response is difficult due to problems in the chemical and conceptual definition of the flow components.     

Seasonal changes in runoff generation in a small forested mountain catchment

This study aimed to investigate the seasonal variability of runoff generation processes, the sources of stream water, and the controls on the contribution of event water to streamflow for a small forested catchment in the Italian pre‐Alps. Hydrometric, isotopic, and electrical conductivity data collected between August 2012 and August 2013 revealed a marked seasonal variability in runoff responses. Noticeable differences in runoff coefficients and hydrological dynamics between summer and fall/spring rainfall events were related to antecedent moisture conditions and event size. Two‐component and three‐component hydrograph separation and end‐member mixing analysis showed an increase in event water contributions to streamflow with event size and average rainfall intensity. Event water fractions were larger during dry conditions in the summer, suggesting that stormflow generation in the summer consisted predominantly of direct channel precipitation and some saturated overland flow from the riparian zone. On the contrary, groundwater and hillslope soil water contributions dominated the streamflow response during wet conditions in fall. Seasonal differences were also noted between event water fractions computed based on isotopic and electrical conductivity data, likely because of the dilution effect during the wetter months. Copyright © 2014 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.     

Stormflow generation and flowpath characteristics in an Amazonian rainforest catchment

The Amazon basin covers an area of roughly 7 × 10⁶ km² and encompasses diverse soil – landscape types with potentially differing hydrological behaviour. This study was conducted in the Ultisol landscape of the western Amazon basin in Peru. Processes of stormflow generation were investigated on an event basis in a first‐order rainforest catchment to establish a causal link between soil physical and precipitation characteristics, hillslope flowpaths and stormflow hydrograph attributes. A sharp decrease in soil hydraulic conductivity with depth and high rainfall intensity and frequency favour rapid near‐surface flowpaths, mainly in the form of saturation‐excess overland flow and return flow. The latter results in an almost random occurrence of overland flow, with no obvious topographic control. Hillslope flowpaths do not vary much with respect to the hydrograph attributes time of rise, response time, lag time and centroid lag time. They have the same response time as streamflow, but a somewhat lower time of rise and significantly shorter lag times. The recession constant for hillslope hydrographs is about 10 min, in contrast to the streamflow recession constants of 28, 75 and 149 min. Stormflow generation in this Ultisol rainforest catchment differs strongly from that reported for Oxisol rainforest catchments. These two soilscapes may define a spectrum of possible catchment hydrological behaviour in the Amazon basin. Copyright © 2000 John Wiley & Sons, Ltd.     

Effect of bedrock flow on catchment rainfall‐runoff characteristics and the water balance in forested catchments in Tanzawa Mountains,Japan

In this paper, we examined the role of bedrock groundwater discharge and recharge on the water balance and runoff characteristics in forested headwater catchments. Using rigorous observations of catchment precipitation, discharge and streamwater chemistry, we quantified net bedrock flow rates and contributions to streamwater runoff and the water balance in three forested catchments (second‐order to third‐order catchments) underlain by uniform bedrock in Japan. We found that annual rainfall in 2010 was 3130 mm. In the same period, annual discharge in the three catchments varied from 1800 to 3900 mm/year. Annual net bedrock flow rates estimated by the chloride mass balance method at each catchment ranged from ?1600 to 700 mm/year. The net bedrock flow rates were substantially different in the second‐order and third‐order catchments. During baseflow, discharge from the three catchments was significantly different; conversely, peak flows during large storm events and direct runoff ratios were not significantly different. These results suggest that differences in baseflow discharge rates, which are affected by bedrock flow and intercatchment groundwater transfer, result in the differences in water balance among the catchments. This study also suggests that in these second‐order to third‐order catchments, the drainage area during baseflow varies because of differences between the bedrock drainage area and surface drainage area, but that the effective drainage area during storm flow approaches the surface drainage area. Copyright © 2012 John Wiley & Sons, Ltd.     

Identifying and linking source areas of flow and P transport in dairy‐grazed headwater catchments,North Island,New Zealand

We examined the applicability of the critical‐source area (CSA) concept to the dairy‐grazed 192‐ha Upper Toenepi catchment and its 8·7‐ha Kiwitahi sub‐catchment, New Zealand. We evaluated if phosphorus (P) transport from land into stream is dominated by saturation‐excess (SE) and infiltration‐excess (IE) runoff during stormflow and by sub‐surface (<1·5 m depth) flows during baseflow. We measured stream flow and shallow groundwater levels, collected monthly stream, tile drain (TDA) and groundwater samples, and flow‐proportional stream samples from the Kiwitahi sub‐catchment, and determined their dissolved reactive phosphorus (DRP) and total phosphorus (TP) concentrations. In the Kiwitahi sub‐catchment, during storm events, IE contributions were significant. Contributions from SE appeared significant in the Upper Toenepi catchment. However, in both catchments, sub‐surface contributions dominated stormflow and baseflow periods. Absence of water table at the surface and the water table gradient towards the stream indicated that P transport during events was not limited to surface runoff. The dynamics of the groundwater table and the occurrence of SE areas were influenced by proximity to the stream and hillslope positions. Baseflow accounted for 42% of the annual flow in the Kiwitahi sub‐catchment, and contributed 37 and 52% to the DRP and TP loads, respectively. The P transport during baseflow appeared equally important as P losses from CSAs during stormflow. The close resemblance in P levels between groundwater and stream samples during baseflow demonstrates the importance of shallow groundwater for stream flow. In the Upper Toenepi catchment, contributions from effluent ponds (EFFs) dominated P loads. Management strategies should focus on controlling P release from EFFs, and on decreasing Olsen P concentrations in soil to minimize leaching of P via sub‐surface flow to streams. Research is needed to quantify the role of sub‐surface flow as well as to expand management strategies to minimize P transfers during stormflow and baseflow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

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.