Heterogeneity in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream

Spatial and temporal variability in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream was investigated. Four locations in a 150‐m reach of the stream were studied using hydrometric and hydrochemical tracing techniques. A high degree of hydrological connectivity between the riparian hillslope and the stream channel was indicated at two locations, where hydrochemical changes and hydraulic gradients indicated that the hyporheic zone was dominated by upwelling ground water. The chemistry of ground water reflected relatively long residence times and reducing conditions with high levels of alkalinity and conductivity, low dissolved oxygen (DO) and nitrate. At the other locations, connectivity was less evident and, at most times, the hyporheic zone was dominated by downwelling stream water characterized by high DO, low alkalinity and conductivity. Substantial variability in hyporheic chemistry was evident at fine (<10 m) spatial scales and changed rapidly over the course of hydrological events. The nature of the hydrochemical response varied among locations depending on the strength of local ground water influence. It is suggested that greater emphasis on spatial and temporal heterogeneity in ground water–surface water interactions in the hyporheic zone is necessary for a consideration of hydrochemical effects on many aspects of stream ecology. For example, the survival of salmonid eggs in hyporheic gravels varied considerably among the locations studied and was shown to be associated with variation in interstitial chemistry. River restoration schemes and watershed management strategies based only on the surface expression of catchment characteristics risk excluding consideration of potentially critical subsurface processes. Copyright © 2002 John Wiley & Sons, Ltd.     

Temporal and spatial response of hyporheic zone geochemistry to a storm event

Although there has been recent focus on understanding spatial variability in hyporheic zone geochemistry across different morphological units under baseflow conditions, less attention has been paid to temporal responses of hyporheic zone geochemistry to non‐steady‐state conditions. We documented spatial and temporal variability of hyporheic zone geochemistry in response to a large‐scale storm event, Tropical Storm Irene (August 2011), across a pool–riffle–pool sequence along Chittenango Creek in Chittenango, NY, USA. We sampled stream water as well as pore water at 15 cm depth in the streambed at 14 locations across a 30 m reach. Sampling occurred seven times at daily intervals: once during baseflow conditions, once during the rising limb of the storm hydrograph, and five times during the receding limb. Principal component analysis was used to interpret temporal and spatial changes and dominant drivers in stream and pore water geochemistry (n = 111). Results show the majority of spatial variance in hyporheic geochemistry (62%) is driven by differential mixing of stream and ground water in the hyporheic zone. The second largest driver (17%) of hyporheic geochemistry was temporal dilution and enrichment of infiltrating stream water during the storm. Hyporheic sites minimally influenced by discharging groundwater (‘connected’ sites) showed temporal changes in water chemistry in response to the storm event. Connected sites within and upstream of the riffle reflected stream geochemistry throughout the storm, whereas downstream sites showed temporally lagged responses in some conservative and biogeochemically reactive solutes. This suggests temporal changes in hyporheic geochemistry at these locations reflect a combination of changes in infiltrating stream chemistry and hyporheic flowpath length and residence time. The portion of the study area strongly influenced by groundwater discharge increased in size throughout the storm, producing elevated Ca²⁺ and concentrations in the streambed, suggesting zones of localized groundwater inputs expand in response to storms. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatial variability in specific discharge and streamwater chemistry during low flows: Results from snapshot sampling campaigns in eleven Swiss catchments

Catchments consist of distinct landforms that affect the storage and release of subsurface water. Certain landforms may be the main contributors to streamflow during extended dry periods, and these may vary for different catchments in a given region. We present a unique dataset from snapshot field campaigns during low‐flow conditions in 11 catchments across Switzerland to illustrate this. The catchments differed in size (10 to 110 km²), varied from predominantly agricultural lowlands to Alpine areas, and covered a range of physical characteristics. During each snapshot campaign, we jointly measured streamflow and collected water samples for the analysis of major ions and stable water isotopes. For every sampling location (basin), we determined several landscape characteristics from national geo‐datasets, including drainage area, elevation, slope, flowpath length, dominant land use, and geological and geomorphological characteristics, such as the lithology and fraction of quaternary deposits. The results demonstrate very large spatial variability in specific low‐flow discharge and water chemistry: Neighboring sampling locations could differ significantly in their specific discharge, isotopic composition, and ion concentrations, indicating that different sources contribute to streamflow during extended dry periods. However, none of the landscape characteristics that we analysed could explain the spatial variability in specific discharge or streamwater chemistry in multiple catchments. This suggests that local features determine the spatial differences in discharge and water chemistry during low‐flow conditions and that this variability cannot be assessed a priori from available geodata and statistical relations to landscape characteristics. The results furthermore suggest that measurements at the catchment outlet during low‐flow conditions do not reflect the heterogeneity of the different source areas in the catchment that contribute to streamflow.     

Heat, chloride, and specific conductance as ground water tracers near streams

Commonly measured water quality parameters were compared to heat as tracers of stream water exchange with ground water. Temperature, specific conductance, and chloride were sampled at various frequencies in the stream and adjacent wells over a 2-year period. Strong seasonal variations in stream water were observed for temperature and specific conductance. In observation wells where the temperature response correlated to stream water, chloride and specific conductance values were similar to stream water values as well, indicating significant stream water exchange with ground water. At sites where ground water temperature fluctuations were negligible, chloride and/or specific conductance values did not correlate to stream water values, indicating that ground water was not significantly influenced by exchange with stream water. Best-fit simulation modeling was performed at two sites to derive temperature-based estimates of hydraulic conductivities of the alluvial sediments between the stream and wells. These estimates were used in solute transport simulations for a comparison of measured and simulated values for chloride and specific conductance. Simulation results showed that hydraulic conductivities vary seasonally and annually. This variability was a result of seasonal changes in temperature-dependent hydraulic conductivity and scouring or clogging of the streambed. Specific conductance fits were good, while chloride data were difficult to fit due to the infrequent (quarterly) stream water chloride measurements during the study period. Combined analyses of temperature, chloride, and specific conductance led to improved quantification of the spatial and temporal variability of stream water exchange with shallow ground water in an alluvial system.     

The influence of climate on water chemistry states and dynamics in rivers across Australia

For effective water quality management and policy development, spatial variability in the mean concentrations and dynamics of riverine water quality needs to be understood. Using water chemistry (calcium, electrical conductivity, nitrate-nitrite, soluble reactive phosphorus, total nitrogen, total phosphorus and total suspended solids) data for up to 578 locations across the Australian continent, we assessed the impact of climate zones (arid, Mediterranean, temperate, subtropical, tropical) on (i) inter-annual mean concentration and (ii) water chemistry dynamics as represented by constituent export regimes (ratio of the coefficients of variation of concentration and discharge) and export patterns (slope of the concentration-discharge relationship). We found that inter-annual mean concentrations vary significantly by climate zones and that spatial variability in water chemistry generally exceeds temporal variability. However, export regimes and patterns are generally consistent across climate zones. This suggests that intrinsic properties of individual constituents rather than catchment properties determine export regimes and patterns. The spatially consistent water chemistry dynamics highlights the potential to predict riverine water quality across the Australian continent, which can support national riverine water quality management and policy development.     

A Permanent Multilevel Monitoring and Sampling System in the Coastal Groundwater Mixing Zones

To study the spatial and temporal variability of water dynamics and chemical reactions within the coastal groundwater mixing zones (CGMZs), high‐resolution periodical and spatial groundwater sampling within CGMZs is needed. However, current samplers and sampling systems may require heavy driving machines to install. There is also possible contamination from the metal materials for current samplers and sampling systems. Here, a permanent multilevel sampling system is designed to sample coastal groundwater within CGMZs. This cost‐effective system consists of metal‐free materials and can be installed easily. The system is tested in Po Sam Pai and Tingkok, Tolo Harbor and Hong Kong. Major ions, nutrients, stable isotopes and radium and radon isotopes were analyzed and the data provided scientific information to study the fresh‐saltwater interface fluctuations, and temporal variations and spatial heterogeneity of geochemical processes occurred within CGMZs. The reliable spatial and temporal data from the sampling system demonstrate that the system functions well and can provide scientific data for coastal aquifer studies.     

Influence of water level fluctuation on the spatial dynamic of the water chemistry at Lake Fertõ/Neusiedler See

Lake Fertõ/Neusiedler See is a steppe lake situated on the Hungarian-Austrian border with a surface area of 309 km². It is an extremely shallow lake with regulated outflow. 54% of the whole lake and 85% of the Hungarian part (75 km²) is covered by reed. There are numerous reedless areas (inner ponds) of variable size within the reed belt, which is enmeshed with canals connecting the inner ponds with the open water areas. Fundamental differences were recorded in the water chemistry of the three water type areas, open water, canals and inner ponds. The degree of these differences is basically determined by the distance from the open water areas and the water level of the actual years. An extremely high, an average and an extremely low water level year was examined, in order to present the spatial and temporal differences in the water chemistry of the investigated water bodies. Differences due to the water level fluctuation between years in the three different water types are discussed on the basis of electrical conductivity, cation, and anion concentration changes of the water. Our results support the concept of high instability of Lake Fertõ and describe habitats with remarkable temporal and spatial variability. The high influence of water level on habitat conditions at different sites of the lake recommend a carefull water level regulation strategy.     

Stratigraphic control of flow and transport characteristics

Ground water flow and travel time are dependent on stratigraphic architecture, which is governed by competing processes that control the spatial and temporal distribution of accommodation and sediment supply. Accommodation is the amount of space in which sediment may accumulate as defined by the difference between the energy gradient and the topographic surface. The temporal and spatial distribution of accommodation is affected by processes that change the distribution of energy (e.g., sea level or subsidence). Fluvial stratigraphic units, generated by FLUVSIM (a stratigraphic simulator based on accommodation and sediment supply), with varying magnitudes and causes of accommodation, were incorporated into a hydraulic regime using MODFLOW (a ground water flow simulator), and particles were tracked using MODPATH (a particle-tracking algorithm). These experiments illustrate that the dominant type of accommodation process influences the degree of continuity of stratigraphic units and thus affects ground water flow and transport. When the hydraulic gradient is parallel to the axis of the fluvial system in the depositional environment, shorter travel times occur in low-total accommodation environments and longer travel times in high-total accommodation environments. Given the same total accommodation, travel times are longer when sea-level change is the dominant process than those in systems dominated by subsidence.     

Stratigraphic control of flow and transport characteristics

Ground water flow and travel time are dependent on stratigraphic architecture, which is governed by competing processes that control the spatial and temporal distribution of accommodation and sediment supply. Accommodation is the amount of space in which sediment may accumulate as defined by the difference between the energy gradient and the topographic surface. The temporal and spatial distribution of accommodation is affected by processes that change the distribution of energy (e.g., sea level or subsidence). Fluvial stratigraphic units, generated by FLUVSIM (a stratigraphic simulator based on accommodation and sediment supply), with varying magnitudes and causes of accommodation, were incorporated into a hydraulic regime using MODFLOW (a ground water flow simulator), and particles were tracked using MODPATH (a particle-tracking algorithm). These experiments illustrate that the dominant type of accommodation process influences the degree of continuity of stratigraphic units and thus affects ground water flow and transport. When the hydraulic gradient is parallel to the axis of the fluvial system in the depositional environment, shorter travel times occur in low-total accommodation environments and longer travel times in high-total accommodation environments. Given the same total accommodation, travel times are longer when sea-level change is the dominant process than those in systems dominated by subsidence.     

Sensitivity to acidification of subalpine ponds and lakes in north‐western Colorado

Although acidifying deposition in western North America is lower than in many parts of the world, many high‐elevation ecosystems there are extremely sensitive to acidification. Previous studies determined that the Mount Zirkel Wilderness Area (MZWA) has the most acidic snowpack and aquatic ecosystems that are among the most sensitive in the region. In this study, spatial and temporal variability of ponds and lakes in and near the MZWA were examined to determine their sensitivity to acidification and the effects of acidic deposition during and after snowmelt. Within the areas identified as sensitive to acidification based on bedrock types, there was substantial variability in acid‐neutralizing capacity (ANC), which was related to differences in hydrological flowpaths that control delivery of weathering products to surface waters. Geological and topographic maps were of limited use in predicting acid sensitivity because their spatial resolution was not fine enough to capture the variability of these attributes for lakes and ponds with small catchment areas. Many of the lakes are sensitive to acidification (summer and autumn ANC < 100 µeq L^?1), but none of them appeared to be threatened immediately by episodic or chronic acidification. In contrast, 22 ponds had minimum ANC < 30 µeq L^?1, indicating that they are extremely sensitive to acidic deposition and could be damaged by episodic acidification, although net acidity (ANC < 0) was not measured in any of the ponds during the study. The lowest measured pH value was 5·4, and pH generally remained less than 6·0 throughout early summer in the most sensitive ponds, indicating that biological effects of acidification are possible at levels of atmospheric deposition that occurred during the study. The aquatic chemistry of lakes was dominated by atmospheric deposition and biogeochemical processes in soils and shallow ground water, whereas the aquatic chemistry of ponds was also affected by organic acids and biogeochemical processes in the water column and at the sediment–water interface. These results indicate that conceptual and mechanistic acidification models that have been developed for lakes and streams may be inadequate for predicting acidification in less‐understood systems such as ponds. Copyright © 2004 John Wiley & Sons, Ltd.     

A geographic data model for representing ground water systems

The Arc Hydro ground water data model is a geographic data model for representing spatial and temporal ground water information within a geographic information system (GIS). The data model is a standardized representation of ground water systems within a spatial database that provides a public domain template for GIS users to store, document, and analyze commonly used spatial and temporal ground water data sets. This paper describes the data model framework, a simplified version of the complete ground water data model that includes two-dimensional and three-dimensional (3D) object classes for representing aquifers, wells, and borehole data, and the 3D geospatial context in which these data exist. The framework data model also includes tabular objects for representing temporal information such as water levels and water quality samples that are related with spatial features.     

The spatial and temporal variability of groundwater recharge in a forested basin in northern Wisconsin

Recharge varies spatially and temporally as it depends on a wide variety of factors (e.g. vegetation, precipitation, climate, topography, geology, and soil type), making it one of the most difficult, complex, and uncertain hydrologic parameters to quantify. Despite its inherent variability, groundwater modellers, planners, and policy makers often ignore recharge variability and assume a single average recharge value for an entire watershed. Relatively few attempts have been made to quantify or incorporate spatial and temporal recharge variability into water resource planning or groundwater modelling efforts. In this study, a simple, daily soil–water balance model was developed and used to estimate the spatial and temporal distribution of groundwater recharge of the Trout Lake basin of northern Wisconsin for 1996–2000 as a means to quantify recharge variability. For the 5 years of study, annual recharge varied spatially by as much as 18 cm across the basin; vegetation was the predominant control on this variability. Recharge also varied temporally with a threefold annual difference over the 5‐year period. Intra‐annually, recharge was limited to a few isolated events each year and exhibited a distinct seasonal pattern. The results suggest that ignoring recharge variability may not only be inappropriate, but also, depending on the application, may invalidate model results and predictions for regional and local water budget calculations, water resource management, nutrient cycling, and contaminant transport studies. Recharge is spatially and temporally variable, and should be modelled as such. Copyright © 2009 John Wiley & Sons, Ltd.     

Deducing the spatial variability of exchange within a longitudinal channel water balance

Developing an appropriate data collection scheme to infer stream–subsurface interactions is not trivial due to the spatial and temporal variability of exchange flowpaths. Within the context of a case study, this paper presents the results from a number of common data collection techniques ranging from point to reach scales used in combination to better understand the spatial complexity of subsurface exchanges, infer the hydrologic conditions where individual influences of hyporheic and groundwater exchange components on stream water can be characterized, and determine where gaps in information arise. We start with a tracer‐based, longitudinal channel water balance to quantify hydrologic gains and losses at a sub‐reach scale nested within two consecutive reaches. Next, we look at groundwater and stream water surface levels, shallow streambed vertical head gradients, streambed and aquifer hydraulic conductivities, water chemistry, and vertical flux rates estimated from streambed temperatures to provide more spatially explicit information. As a result, a clearer spatial understanding of gains and losses was provided, but some limitations in interpreting results were identified even when combining information collected over various scales. Due to spatial variability of exchanges and areas of mixing, each technique frequently captured a combination of groundwater and hyporheic exchange components. Ultimately, this study provides information regarding technique selection, emphasizes that care must be taken when interpreting results, and identifies the need to apply or develop more advanced methods for understanding subsurface exchanges. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatio‐temporal variability of the isotopic input signal in a partly forested catchment: Implications for hydrograph separation

The isotopic composition of precipitation (D and ¹⁸O) has been widely used as an input signal in water tracer studies. Whereas much recent effort has been put into developing methodologies to improve our understanding and modelling of hydrological processes (e.g., transit‐time distributions or young water fractions), less attention has been paid to the spatio‐temporal variability of the isotopic composition of precipitation, used as input signal in these studies. Here, we investigated the uncertainty in isotope‐based hydrograph separation due to the spatio‐temporal variability of the isotopic composition of precipitation. The study was carried out in a Mediterranean headwater catchment (0.56 km²). Rainfall and throughfall samples were collected at three locations across this relatively small catchment, and stream water samples were collected at the outlet. Results showed that throughout an event, the spatial variability of the input signal had a higher impact on hydrograph separation results than its temporal variability. However, differences in isotope‐based hydrograph separation determined preevent water due to the spatio‐temporal variability were different between events and ranged between 1 and 14%. Based on catchment‐scale isoscapes, the most representative sampling location could also be identified. This study confirms that even in small headwater catchments, spatio‐temporal variability can be significant. Therefore, it is important to characterize this variability and identify the best sampling strategy to reduce the uncertainty in our understanding of catchment hydrological processes.     

Spatial and temporal variability of pore‐water pressures in residual soil slopes in a tropical climate

It is critical to understand and quantify the temporal and spatial variability in hillslope hydrological data in order to advance hillslope hydrological studies, evaluate distributed parameter hydrological models, analyse variability in hydrological response of slopes and design efficient field data sampling networks. The spatial and temporal variability of field‐measured pore‐water pressures in three residual soil slopes in Singapore was investigated using geostatistical methods. Parameters of the semivariograms, namely the range, sill and nugget effect, revealed interesting insights into the spatial structure of the temporal situation of pore‐water pressures in the slopes. While informative, mean estimates have been shown to be inadequate for modelling purposes, indicator semivariograms together with mean prediction by kriging provide a better form of model input. Results also indicate that significant temporal and spatial variability in pore‐water pressures exists in the slope profile and thereby induces variability in hydrological response of the slope. Spatial and temporal variability in pore‐water pressure decreases with increasing soil depth. The variability decreases during wet conditions as the slope approaches near saturation and the variability increases with high matric suction development following rainfall periods. Variability in pore‐water pressures is greatest at shallow depths and near the slope crest and is strongly influenced by the combined action of microclimate, vegetation and soil properties. Copyright © 2002 John Wiley & Sons, Ltd.     

Analysis of recharge-induced geochemical change in a contaminated aquifer

Recharge events that deliver electron acceptors such as O2, NO3, SO4, and Fe3+ to anaerobic, contaminated aquifers are likely important for natural attenuation processes. However, the specific influence of recharge on (bio)geochemical processes in ground water systems is not well understood. The impact of a moderate-sized recharge event on ground water chemistry was evaluated at a shallow, sandy aquifer contaminated with waste fuels and chlorinated solvents. Multivariate statistical analyses coupled with three-dimensional visualization were used to analyze ground water chemistry data (including redox indicators, major ions, and physical parameters) to reveal associations between chemical parameters and to infer processes within the ground water plume. Factor analysis indicated that dominant chemical associations and their interpreted processes (anaerobic and aerobic microbial processes, mineral precipitation/dissolution, and temperature effects) did not change significantly after the spring recharge event of 2000. However, the relative importance of each of these processes within the plume changed. After the recharge event, the overall importance of aerobic processes increased from the fourth to the second most important factor, representing the variability within the data set. The anaerobic signatures became more complex, suggesting that zones with multiple terminal electron-accepting processes (TEAPs) likely occur in the same water mass. Three-dimensional visualization of well clusters showed that water samples with similar chemical associations occurred in distinct water masses within the aquifer. Water mass distinctions were not based on dominant TEAPs, suggesting that the recharge effects on TEAPs occurred primarily at the interface between infiltrating recharge water and the aquifer.     

Mapping multi‐year ground water salinity patterns in irrigation areas using time‐series analysis of ground water salinity maps

In irrigation areas, ground water salinity (GWS) levels may vary depending on the amount and quality of irrigation water applied and on the activity of the drainage system. GWS plays a vital role in irrigation systems by influencing whether the environment is suitable for plant growth. Hence, it is necessary to monitor changes in GWS both temporally and spatially. Maps are generally used to visualize this information. However, evaluation of temporal and spatial variations of GWS can be difficult because of the necessity of assessing many maps together to understand both temporal and spatial changes. In this study, a data assessment method that can be used for multi‐year ground water salinity evaluations is presented. The method looks at the spatial and temporal relationships between the main salinity classes present in the study area, their typical locations (i.e. areas where the salinity classes are most frequently located), and the alternate salinity classes in those locations in any of the years of the time series. As a case study, the method was applied to multi‐year (1990–2000) GWS observations in the Mustafakemalpasa irrigation project (19 370 ha) in the Marmara region of north‐western Turkey. Copyright © 2007 John Wiley & Sons, Ltd.     

Use of temperature profiles and stable isotopes to trace flow lines: Nagaoka area, Japan

In this study, we use borehole temperature data and stable isotopes to delineate the flow system and estimate the effect of urbanization in the Nagaoka area of Japan. Temperature profiles were measured four times in observation wells during the period 2000-2001 and compared with those measured in the same wells during the period 1977-1983 (Taniguchi 1986). Water was sampled in both observation and pumping wells during the same period. The temporal and spatial variability in temperature indicate clearly the effect of urban warming and heavy pumping on the ground water system. Urban warming caused higher temperatures recently as compared to the older values, and pumping caused induced recharge from the river to the ground water. The stable isotope data show the ground water flow system is divided into shallow, intermediate, and deep systems, and that land use and infiltration rate are affecting the shallow flow system.     

Synoptic sampling reveals spatial stability in flow and chemistry patterns despite large seasonal variability in a subarctic mountain catchment

Seasonality plays a critical role in cold mountain regions as variation in air temperature, ground thermal status, and precipitation phase alter the rate, timing and magnitude of hydrological and chemical transport. Additionally, cold mountain catchments can have highly variable topography, geology, permafrost, and landcover, which intrinsically add to this irregularity. Understanding how external and internal variability act to control mass fluxes requires sampling at a high spatial resolution over time, which rarely occurs in cold remote regions. In this work, we conduct five snapshot sampling surveys across 34 subcatchments during the ice-free period in Wolf Creek Research Basin (a mesoscale montane subarctic catchment) and two additional winter surveys across a subset of sites to assess the drivers of variability in stream chemistry and discharge. We sampled for specific conductance (SpC), major ions, and dissolved organic carbon (DOC) and used statistical metrics and Bayesian mixing analysis to quantify patterns of flow and chemistry across space and time. Our results indicate patterns in both flow and chemistry remain largely consistent across seasons for all solutes. However, there was weaker correlation of chemistry between sites, suggesting asynchronous behaviour within the catchment. There was evidence of increasing production of ions and DOC along the stream network during high spring flows but not during low flows. Although concentrations and flows exhibit high seasonality in subarctic mountains, this seasonal variability does not alter spatial patterns that arise from highly variable catchment characteristics.     

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.