Exploring effects of rainfall intensity and duration on soil erosion at the catchment scale using openLISEM: Prado catchment,SE Spain

In semi‐arid areas, high‐intensity rainfall events are often held responsible for the main part of soil erosion. Long‐term landscape evolution models usually use average annual rainfall as input, making the evaluation of single events impossible. Event‐based soil erosion models are better suited for this purpose but cannot be used to simulate longer timescales and are usually applied to plots or small catchments. In this study, the openLISEM event‐based erosion model was applied to the medium‐sized (～50 km²) Prado catchment in SE Spain. Our aim was to (i) test the model's performance for medium‐sized catchments, (ii) test the ability to simulate four selected typical Mediterranean rainfall events of different magnitude and (iii) explore the relative contribution of these different storms to soil erosion using scenarios of future climate variability. Results show that because of large differences in the hydrologic response between storms of different magnitudes, each event needed to be calibrated separately. The relation between rainfall event characteristics and the calibration factors might help in determining optimal calibration values if event characteristics are known. Calibration of the model features some drawbacks for large catchments due to spatial variability in K_sat values. Scenario calculations show that although ～50% of soil erosion occurs as a result of high frequency, low‐intensity rainfall events, large‐magnitude, low‐frequency events potentially contribute significantly to total soil erosion. The results illustrate the need to incorporate temporal variability in rainfall magnitude–frequency distributions in landscape evolution models. Copyright © 2011 John Wiley & Sons, Ltd.     

An approach to catchment‐scale groundwater nitrate risk assessment from diffuse agricultural sources: a case study in the Upper Bann,Northern Ireland

DRASTIC has drawbacks in groundwater risk assessment that are important in guiding activities to prevention agricultural diffuse groundwater pollution. This paper presents an improved and GIS‐based D‐DRASTIC approach for groundwater nitrate risk assessment from diffuse agricultural sources based on DRASTIC. D‐DRASTIC considers the risk concept, nitrate loading, pollutant transport with runoff, depth to water, net recharge, aquifer media, soil media, topography, impact of the vadose zone media, and the hydraulic conductivity of the aquifer. D‐DRASTIC was developed within an ArcGIS environment and applied to the Upper Bann Catchment, Northern Ireland as a case study. D‐DRASTIC shows that ‘very high’ and ‘high’ zones of groundwater nitrate risk occupy 5% and 11% of the case study area, respectively. When considering groundwater pollution sources and pathways, the results using D‐DRASTIC are helpful in guiding the activities of groundwater pollution prevention at the catchment scale in the context of better implementation of the EU Water Framework Directive. Copyright © 2008 John Wiley & Sons, Ltd.     

Assessing the impacts of climate change on water quantity and quality modelling in small Slovenian Mediterranean catchment – lesson for policy and decision makers

The study aims to address the long‐term impacts of six different downscaled Regional Climate Models (RCM) climate models on the quantity (river flow) and quality (sediment load, total nitrogen load and total phosphorus load) state of surface waters in the river Reka catchment, in the northern Mediterranean. Mediterranean areas are – due to high population density, favourable natural conditions for agriculture, limited water resources, diverse ecosystems biodiversity and expected climate change impacts – a global hotspot in climate research. Additionally, the study area lies on the border with the alpine climate zone, with a strong orographic effect on weather patterns. The location, and a wide range of studied parameters, provides an interesting insight into how various emerging climate change models may impact the status of surface waters and procedures for the governance of water resources. The study contributes to the knowledge and understanding of the climate change impact on the local catchment level, using the ensemble of the RCMs. It opens discussion about the impact of RCM selection on modelling climate changes with catchment models like Soil and Water Assessment Tool. This article also questions the usability of the results for the policy and decision makers in relation to the implementation of the results into short or long‐term water strategies or water/river management plans. Copyright © 2015 John Wiley & Sons, Ltd.     

Sources and transport of nitrate constrained by the isotopic technique in a karst catchment: an example from Southwest China

The concentration and isotopic composition of nitrate were analyzed to improve an understanding of nitrate sources and transformation processes in a typical karstic agricultural field in the Houzhai catchment, Guizhou Province, Southwest China. The results revealed that no distinct spatial pattern of content and isotopic composition of nitrate exists in this karst catchment. Nitrate in surface stream (SFS) had slightly lighter isotopic composition and lower concentration compared with nitrate in subterranean stream (STS) during the dry season. Concentrations of SFS nitrate increased to concentrations similar to those of STS during the wet season. The isotopic values indicated that nitrate were mainly impacted by manure sources during the dry season and influenced by a mix of chemical fertilizer and manure during the wet season. The denitrification rates were roughly estimated based on the isotopic compositions of nitrate after considering volatilization and ignoring assimilation. The calculated result showed that approximately one fifth of nitrate load was removed by denitrification in the catchment. Annual nitrate flux from the outlets accounted for 14.2% of applied total fertilizers used in the catchment, approximately 85% of total transported flux from the catchment in the wet season. Furthermore, chemical weathering processes were enhanced by using nitrogen fertilizer because liberated protons and enhanced HCO₃^? flux were produced through by nitrification. Copyright © 2014 John Wiley & Sons, Ltd.     

Nitrate removal and young stream water fractions at the catchment scale

Despite extensive research on nitrate export and removal, nutrient contamination remains a major threat to water bodies worldwide. At the local scale, nitrate removal is governed by biogeochemical conditions that vary in space and time, making integration to entire landscapes critical. Water transit times have often been used to describe solute transport, but the relation between water age and nitrate removal at the catchment scale is still poorly understood. We test the hypothesis that nitrate removal peaks when the fraction of young water in discharge is at its minimum, because nitrate removal occurs mostly under dry conditions where deeper, older groundwater dominates streamflow. We tested this hypothesis by exploring a detailed water quality record from the Kervidy–Naizin catchment (FR) and comparing the dynamics of nitrate to those of a conservative solute (chloride). We find that estimates of nitrate removal are consistent with previous estimates at the site and they show a good (inverse) correlation with the fraction of streamflow that is younger than 2.5 months. However, this young water fraction cannot be used to predict nitrate removal in the winter–spring period, when no removal is observed regardless of streamflow age. While this leads us to reject our hypothesis during the winter period, it also suggests that water age distributions and their correlation with nitrate removal can possibly reveal distinct sources of stream water at different hydrologic regimes and relevant biogeochemical reactions.     

From precipitation to runoff: stable isotopic fractionation effect of glacier melting on a catchment scale

Stable isotope variability and fractionation associated with transformation of precipitation/accumulation to firn to glacial river water is critical in a variety of climatic, hydrological and paleoenvironmental studies. This paper documents the modification of stable isotopes in water from precipitation to glacier runoff in an alpine catchment located in the central Tibetan Plateau. Isotopic changes are observed by sampling firnpack profiles, glacier surface snow/ice, meltwater on the glacier surface and catchment river water at different times during a melt season. Results show the isotopic fractionation effects associated with glacier melt processes. The slope of the δD‐δ¹⁸O regression line and the deuterium excess values decreased from the initial precipitation to the melt‐impacted firnpack (slope from 9.3 to 8.5 and average d‐excess from 13.4‰ to 7.4‰). The slope of the δD‐δ¹⁸O line further decreased to 7.6 for the glacier runoff water. The glacier surface snow/ice from different locations, which produces the main runoff, had the same δD‐δ¹⁸O line slope but lower deuterium excess (by 3.9‰) compared to values observed in the firnpack profile during the melt season. The δD‐δ¹⁸O regression line for the river water exhibited a lower slope compared to the surface snow/ice samples, although they were closely located on the δD‐δ¹⁸O plot. Isotope values for the river and glacier surface meltwater showed little scatter around the δD‐δ¹⁸O regression line, although the samples were from different glaciers and were collected on different days. Results indicate a high consistency of isotopic fractionation in the δD‐δ¹⁸O relationships, as well as a general consistency and temporal covariation of meltwater isotope values at the catchment scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Catchment‐scale sensitivity and uncertainty in water quality modelling

In this paper, we assess the performance of the catchment model SIMulated CATchment model (SIMCAT), to predict nitrate and soluble reactive phosphorus concentrations against four monitoring regimes with different spatial and temporal sampling frequencies. The Generalised Likelihood Uncertainty Estimation (GLUE) uncertainty framework is used, along with a general sensitivity analysis to understand relative parameter sensitivity. Improvements to model calibration are explored by introducing more detailed process representation using the Integrated Catchments model (INCA) water quality model, driven by the European hydrological predictions for the environment model. The results show how targeted sampling of headwater watercourses upstream of point discharges is essential for calibrating diffuse loads and can exert a strong influence on the whole‐catchment model performance. Further downstream, if the point discharges and loads are accurately represented, then the improvement in the catchment‐scale model performance is relatively small as more calibration points are added or frequency is increased. The higher‐order, dynamic model integrated catchments model of phosphorus dynamics, which incorporates sediment and biotic interaction, resulted in improved whole‐catchment performance over SIMCAT, although there are still large epistemic uncertainties from land‐phase export coefficients and runoff. However, the very large sampling errors in routine monitoring make it difficult to invest confidence in the modelling, especially because we know phosphorous transport to be very episodic and driven by high flow conditions for which there are few samples. The environmental modelling community seems to have been stuck in this position for some time, and whilst it is useful to use an uncertainty framework to highlight these issues, it has not widely been adopted, perhaps because there is no clear mechanism to allow uncertainties to influence investment decisions. This raises the question as to whether it might better place a cost on uncertainty and use this to drive more data collection or improved models, before making investment decisions concerning, for example, mitigation strategies. Copyright © 2016 John Wiley & Sons, Ltd.     

Development and application of a distributed catchment-scale hydrological model for the River Ythan,NE Scotland

Mathematical models are being used to develop a decision support system for integrated management of the Ythan catchment in NE Scotland. One component of this has involved the development of a distributed catchment-scale hydrological model. The model is based on subsurface flow routing and calculates the contribution to stream flow from each 50 m×50 m cell in the 548 km² catchment. It uses two topographic parameters, slope and distance to stream following the main line of flow, and five physical parameters. The topographic analysis and distributed flow accumulation are performed by linking the single cell model with a geographic information system. Preliminary results from a three-year simulation of daily flows indicate that the model successfully predicts the main characteristics of the catchment flow. © 1998 John Wiley & Sons, Ltd.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 3. The large catchment model

This paper presents an application of a long-term, large catchment-scale, water balance model developed to predict the effects of forest clearing in the south-west of Western Australia. The conceptual model simulates the basic daily water balance fluxes in forested catchments before and after clearing. The large catchment is divided into a number of sub-catchments (1–5 km² in area), which are taken as the fundamental building blocks of the large catchment model. The responses of the individual subcatchments to rainfall and pan evaporation are conceptualized in terms of three inter-dependent subsurface stores A, B and F, which are considered to represent the moisture states of the subcatchments. Details of the subcatchment-scale water balance model have been presented earlier in Part 1 of this series of papers. The response of any subcatchment is a function of its local moisture state, as measured by the local values of the stores. The variations of the initial values of the stores among the subcatchments are described in the large catchment model through simple, linear equations involving a number of similarity indices representing topography, mean annual rainfall and level of forest clearing. The model is applied to the Conjurunup catchment, a medium-sized (39·6 km²) catchment in the south-west of Western Australia. The catchment has been heterogeneously (in space and time) cleared for bauxite mining and subsequently rehabilitated. For this application, the catchment is divided into 11 subcatchments. The model parameters are estimated by calibration, by comparing observed and predicted runoff values, over a 18 year period, for the large catchment and two of the subcatchments. Excellent fits are obtained.     

Hydrological behaviour and modelling of a volcanic tropical cultivated catchment

The hydrological behaviour of the cultivated Féfé catchment (17·8 ha) on the tropical volcanic island of Guadeloupe was studied to identify flow paths, to quantify water fluxes, and finally, to build a lumped model to simulate discharge and piezometer levels. The approach combined two steps, an experimental step and a modelling step, which covered two time scales, the annual and the storm event scale. The hydrological measurements were conducted over 2 years. The Féfé catchment is characterized by heavy rainfall (4229 mm year^?1) on permeable Andosols; the results showed that underground flow paths involved two overlapping aquifers, and that the annual water balance in 2003 was shared among outflows of the deep aquifer (42%), evapotranspiration (31%), and streamflow (27%). On the event scale, the surface runoff coefficient ranges between 6·2% and 24·4% depending on antecedent dry or wet moisture conditions. Hortonian overland flow predominated over subsurface and saturation overland flow processes. Recharge of the shallow aquifer is mainly governed by a constant infiltration capacity of the Andosols with depth in the vadose zone. Outflows of this shallow aquifer were the baseflow of the main stream and the recharge of the deep aquifer. Volcanic deposits at Féfé promoted the underground flow path, and cultivated areas seemed to explain the high stormflow values relative to other tropical small catchments under rain forest. A conceptual lumped model integrating runoff, infiltration, evapotranspiration, and fluctuations of the two overlapping aquifers was developed. The model has six parameters and was calibrated and validated on the hydrograph at the outlet and on the two piezometers of the shallow and the deep aquifers. The results show fair to good agreement between measured and simulated variables, and consequently, the model was consistent with the main hydrological processes observed from experimental results in wet conditions. Copyright © 2008 John Wiley & Sons, Ltd.     

Improving bank erosion modelling at catchment scale by incorporating temporal and spatial variability

Bank erosion can contribute a significant portion of the sediment budget within temperate catchments, yet few catchment scale models include an explicit representation of bank erosion processes. Furthermore, representation is often simplistic resulting in an inability to capture realistic spatial and temporal variability in simulated bank erosion. In this study, the sediment component of the catchment scale model SHETRAN is developed to incorporate key factors influencing the spatio‐temporal rate of bank erosion, due to the effects of channel sinuosity and channel bank vegetation. The model is applied to the Eden catchment, north‐west England, and validated using data derived from a GIS methodology. The developed model simulates magnitudes of total catchment annual bank erosion (617–4063 t y^‐1) within the range of observed values (211–4426 t yr^‐1). In addition, the model provides both greater inter‐annual and spatial variability of bank eroded sediment generation when compared with the basic model, and indicates a potential 61% increase of bank eroded sediment as a result of temporal flood clustering. The approach developed within this study can be used within a number of distributed hydrologic models and has general applicability to temperate catchments, yet further development of model representation of bank erosion processes is required. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

TOPCAT‐NP: a minimum information requirement model for simulation of flow and nutrient transport from agricultural systems

Future catchment planning requires a good understanding of the impacts of land use and management, especially with regard to nutrient pollution. A range of readily usable tools, including models, can play a critical role in underpinning robust decision‐making. Modelling tools must articulate our process understanding, make links to a range of catchment characteristics and scales and have the capability to reflect future land‐use management changes. Hence, the model application can play an important part in giving confidence to policy makers that positive outcomes will arise from any proposed land‐use changes. Here, a minimum information requirement (MIR) modelling approach is presented that creates simple, parsimonious models based on more complex physically based models, which makes the model more appropriate to catchment‐scale applications. This paper shows three separate MIR models that represent flow, nitrate losses and phosphorus losses. These models are integrated into a single catchment model (TOPCAT‐NP), which has the advantage that certain model components (such as soil type and flow paths) are shared by all three MIR models. The integrated model can simulate a number of land‐use activities that relate to typical land‐use management practices. The modelling process also gives insight into the seasonal and event nature of nutrient losses exhibited at a range of catchment scales. Three case studies are presented to reflect the range of applicability of the model. The three studies show how different runoff and nutrient loss regimes in different soil/geological and global locations can be simulated using the same model. The first case study models intense agricultural land uses in Denmark (Gjern, 114 km²), the second is an intense agricultural area dominated by high superphosphate applications in Australia (Ellen Brook, 66 km²) and the third is a small research‐scale catchment in the UK (Bollington Hall, 2 km²). Copyright © 2007 John Wiley & Sons, Ltd.     

Quantifying hydrological budgets and pathways in a small upland catchment using a combined modelling and tracer approach

Hydrological budgets and flow pathways have been quantified for a small upland catchment (1.76 km²) in the northeast of Scotland. Water balance calculations for four subcatchments identified spatial variability within the catchment, with an estimated runoff enhancement of up to 25% for the upper western area, compared with the rest of the catchment. Data from spatial hydrochemical sampling, over a range of flow conditions, were used to identify the principal hillslope runoff mechanisms within the catchment. A hydrochemical mixing analysis revealed that runoff emerging from springs in various locations of the hillslope accounted for a significant proportion of flow in the streams, even during storm events. A hydrological model of the catchment was calibrated using the calculated stream flows for four locations, together with results from the mixing analysis for different time points. The calibrated model was used to predict the temporal variability in contributions to stream flow from the hillslope springs and soil water flows. The overall split ranged from 57%:43% spring water:soil water in the upper eastern subcatchment, to 76%:24% in the upper western subcatchment. Copyright © 2006 John Wiley & Sons, Ltd.     

Attributes and uncertainty of dissolved zinc data from a mined catchment

This study evaluated the attributes and uncertainty of non‐point source pollution data derived from synoptic surveys in a catchment affected by inactive metal mines in order to help to identify and select appropriate methods for data analysis/reporting and information use. Dissolved zinc data from the Upper Animas River Basin, Colorado, USA, were the focus of the study. Zinc was evaluated because concentrations were highest relative to national water quality criteria for brown trout, and zinc had the greatest frequency of criteria exceedances compared with other metals. Data attributes evaluated included measurement and model error, sample size, non‐normality, seasonality and uncertainty. The average measurement errors for discharges, concentrations and loadings were 0·15, 0·1 and 0·18, respectively. The 90 and 95% coefficients of confidence intervals for mean concentrations based on a sample size of four were 0·48 and 0·65, respectively, and ranged between 0·15 and 0·23 for sample sizes greater than 40. Aggregation of data from multiple stations decreased the confidence intervals significantly, but additional aggregation of all data increased them as a result of increasing spatial variability. Unit area loading data were approximately log‐normal. Concentration data were right‐skewed but not log‐normal. Differences in median concentrations were appreciable between snowmelt and both storm flow and baseflow, but not between storm flow and baseflow. Differences in unit area loadings between all flow events were large. It was determined that the average concentration and unit area loading values should be estimated for each flow event because of significant seasonality. Time weighted values generally should be computed if annual information is required. The confidence in average concentrations and unit area loadings is dependent on the computation method used. Both concentrations and loadings can be significantly underestimated on an annual basis when using data from synoptic surveys if the first flush of contaminants during the initial snowmelt runoff period is not sampled. The ambient standard for dissolved zinc for all events was estimated as 1600 μg l⁻¹ using the 85th percentile of observed concentration data, with a 90% confidence interval width of 200 μg l⁻¹. Copyright © 1999 John Wiley & Sons, Ltd.     

Modelling the impact of stormwater source control infiltration techniques on catchment baseflow

Stormwater management increasingly recognises the need to emulate, to the maximum extent possible, the flow regime of receiving waters in their pre‐development state. Hydrological models play a central role in assessing the catchment‐scale impacts of alternative stormwater management strategies. However, because of the complexity of physical processes involved in urban hydrology, particularly subsurface flows, the predictive performance of such models is often low. We investigated how the structure of hydrological models influenced the prediction of urbanisation and stormwater management impacts on baseflow. We calibrated three conceptual models of the same reference catchment and compared the modelled flow regime from different stormwater management scenarios, using each of the three model structures. Scenarios were assessed using six metrics, characterising the whole streamflow regime and in particular baseflow. Although the three models of the reference catchment represented the observed hydrograph well, the most complex structure developed using a thorough diagnostic of the catchment behaviour better captured the change in hydrological regime during dry years. Predictions of baseflow changes due to urbanisation varied significantly according to the model structure. Similarly, the models showed distinct responses to the stormwater management scenarios applied, especially for scenarios involving infiltration of stormwater at source. Our results confirm the importance of predicting the consequences of land use changes with conceptual models that are consistent with the hydrological behaviour of the study catchment. Future work should help to quantify the uncertainties due to model structure and thus provide practical guidance to the use of catchment models for assessing stormwater management strategies. Copyright © 2013 John Wiley & Sons, Ltd.     

A catchment scale evaluation of the SIBERIA and CAESAR landscape evolution models

Landscape evolution models provide a way to determine erosion rates and landscape stability over times scales from tens to thousands of years. The SIBERIA and CAESAR landscape evolution models both have the capability to simulate catchment–wide erosion and deposition over these time scales. They are both cellular, operate over a digital elevation model of the landscape, and represent fluvial and slope processes. However, they were initially developed to solve research questions at different time and space scales and subsequently the perspective, detail and process representation vary considerably between the models. Notably, CAESAR simulates individual events with a greater emphasis on fluvial processes whereas SIBERIA averages erosion rates across annual time scales. This paper describes how both models are applied to Tin Camp Creek, Northern Territory, Australia, where soil erosion rates have been closely monitored over the last 10 years. Results simulating 10 000 years of erosion are similar, yet also pick up subtle differences that indicate the relative strengths and weaknesses of the two models. The results from both the SIBERIA and CAESAR models compare well with independent field data determined for the site over different time scales. Representative hillslope cross‐sections are very similar between the models. Geomorphologically there was little difference between the modelled catchments after 1000 years but significant differences were revealed at longer simulation times. Importantly, both models show that they are sensitive to input parameters and that hydrology and erosion parameter derivation has long‐term implications for sediment transport prediction. Therefore selection of input parameters is critical. This study also provides a good example of how different models may be better suited to different applications or research questions. Copyright © 2010 John Wiley & Sons, Ltd and Commonwealth of Australia     

Opportunities for manipulating catchment water balance by changing vegetation type on a topographic sequence: a simulation study

This simulation study explores opportunities to reduce catchment deep drainage through better matching land use with soil and topography, including the ‘harvesting’ (evapotranspiration) of excess water running on to lower land units. A farming system simulator was coupled with a catchment hydrological framework to enable analysis of climate variability and 11 different land‐use options as they impact the catchment water balance. These land‐use options were arranged in different configurations down a sequence of three hydrologically interconnected slope units (uphill, mid‐slope and valley floor land units) in a subcatchment of Simmons Creek, southern New South Wales, Australia. With annual crops, the valley floor land units were predicted to receive 187 mm year^?1 of run‐on water in addition to annual rainfall in 1 in 10 years, and in excess of 94 mm year^?1 in 1 in 4 years. In this valley floor position, predicted drainage averaged approximately 110 mm year^?1 under annual crops and pastures, whereas permanent tree cover or perennial lucerne was predicted to reduce drainage by up to 99%. The planting of trees or lucerne on the valley floor units could ‘harvest’ run‐on water, reducing drainage for the whole subcatchment with proportionately small reduction in land areas cropped. Upslope land units, even though often having shallower soil, will not necessarily be the most effective locations to plant perennial vegetation for the purposes of recharge reduction. Water harvesting opportunities are site specific, dependent on the amounts and frequency of flows of water to lower landscape units, the amounts and frequency of deep drainage on the different land units, the relative areas of the different land units, and interactions with land use in the different slope positions. Copyright © 2007 John Wiley & Sons, Ltd.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 1. Small catchment water balance model

A long-term water balance model has been developed to predict the hydrological effects of land-use change (especially forest clearing) in small experimental catchments in the south-west of Western Australia. This small catchment model has been used as the building block for the development of a large catchment-scale model, and has also formed the basis for a coupled water and salt balance model, developed to predict the changes in stream salinity resulting from land-use and climate change. The application of the coupled salt and water balance model to predict stream salinities in two small experimental catchments, and the application of the large catchment-scale model to predict changes in water yield in a medium-sized catchment that is being mined for bauxite, are presented in Parts 2 and 3, respectively, of this series of papers. The small catchment model has been designed as a simple, robust, conceptually based model of the basic daily water balance fluxes in forested catchments. The responses of the catchment to rainfall and pan evaporation are conceptualized in terms of three interdependent subsurface stores A, B and F. Store A depicts a near-stream perched aquifer system; B represents a deeper, permanent groundwater system; and F is an intermediate, unsaturated infiltration store. The responses of these stores are characterized by a set of constitutive relations which involves a number of conceptual parameters. These parameters are estimated by calibration by comparing observed and predicted runoff. The model has performed very well in simulations carried out on Salmon and Wights, two small experimental catchments in the Collie River basin in south-west Western Australia. The results from the application of the model to these small catchments are presented in this paper.     

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.