Dynamics of stormflow generation—A hillslope‐scale field study in east‐central Pennsylvania,USA

A 40 m × 20 m mowed, grass hillslope adjacent to a headwater stream within a 26‐ha watershed in east‐central Pennsylvania, USA, was instrumented to identify and map the extent and dynamics of surface saturation (areas with the water table at the surface) and surface runoff source areas. Rainfall, stream flow and surface runoff from the hillslope were recorded at 5‐min intervals from 11 August to 22 November 1998, and 13 April to 12 November 1999. The dynamics of the water table (0 to 45 cm depth from the soil surface) and the occurrence of surface runoff source areas across the hillslope were recorded using specially designed subsurface saturation and surface runoff sensors, respectively. Detailed data analyses for two rainfall events that occurred in August (57·7 mm in 150 min) and September (83·6 mm in 1265 min) 1999, illustrated the spatial and temporal dynamics of surface saturation and surface runoff source areas. Temporal data analyses showed the necessity to measure the hillslope dynamics at time intervals comparable to that of rainfall measurements. Both infiltration excess surface runoff (runoff caused when rainfall intensity exceeds soil infiltration capacity) and saturation excess surface runoff (runoff caused when soil moisture storage capacity is exceeded) source areas were recorded during these rainfall events. The August rainfall event was primarily an infiltration excess surface runoff event, whereas the September rainfall event produced both infiltration excess and saturation excess surface runoff. Occurrence and disappearance of infiltration excess surface runoff source areas during the rainfall events appeared scattered across the hillslope. Analysis of surface saturation and surface runoff data showed that not all surface saturation areas produced surface runoff that reached the stream. Emergence of subsurface flow to the surface during the post‐rainfall periods appeared to be a major flow process dominating the hillslope after the August rainfall event. Copyright © 2002 John Wiley & Sons, Ltd.     

Delineating runoff processes and critical runoff source areas in a pasture hillslope of the Ozark Highlands

The identification of runoff contributing areas would provide the ideal focal points for water quality monitoring and Best Management Practice (BMP) implementation. The objective of this study was to use a field‐scale approach to delineate critical runoff source areas and to determine the runoff mechanisms in a pasture hillslope of the Ozark Highlands in the USA. Three adjacent hillslope plots located at the Savoy Experimental Watershed, north‐west Arkansas, were bermed to isolate runoff. Each plot was equipped with paired subsurface saturation and surface runoff sensors, shallow groundwater wells, H‐flumes and rain gauges to quantify runoff mechanisms and rainfall characteristics at continuous 5‐minute intervals. The spatial extent of runoff source areas was determined by incorporating sensor data into a geographic information‐based system and performing geostatistical computations (inverse distance weighting method). Results indicate that both infiltration excess runoff and saturation excess runoff mechanisms occur to varying extents (0–58% for infiltration excess and 0–26% for saturation excess) across the plots. Rainfall events that occurred 1–5 January 2005 are used to illustrate the spatial and temporal dynamics of the critical runoff source areas. The methodology presented can serve as a framework upon which critical runoff source areas can be identified and managed for water quality protection in other watersheds. Copyright © 2008 John Wiley & Sons, Ltd.     

The hydrology of inland valleys in the sub‐humid zone of West Africa: rainfall‐runoff processes in the M'bé experimental watershed

Inland valleys with wet lowlands are an important water source for farming communities in the sub‐humid zone of West Africa. An inland valley and surrounding contributing watershed area located in the sub‐humid zone near M'bé in central Côte d'Ivoire was instrumented to study surface runoff and base flow mechanisms. Four flumes at different distances down the main stream and more than 100 piezometers were installed. Measurements were taken during two rainfall seasons in 1998 and 1999. Under initial wet conditions, a typical single‐peak hydrograph was observed. Under low antecedent moisture conditions, however, runoff was characterized by a double‐peaked hydrograph. The first peak, which occurred during the storm, was caused by rain falling on the saturated valley bottom. The second peak was delayed by minutes to hours from the first peak and consisted of rain flowing via the subsurface of the hydromorphic zone that surrounds the valley bottom. The duration of the delay was a function of the water table depth in the hydromorphic zone before the storm. The volume of the second peak constituted the largest portion of the stream flow. Copyright © 2003 John Wiley & Sons, Ltd.     

Infiltration,runoff and sediment production in blanket peat catchments: implications of field rainfall simulation experiments

Blanket peat covers the headwaters of many major European rivers. Runoff production in upland blanket peat catchments is flashy with large flood peaks and short lag times; there is minimal baseflow. Little is known about the exact processes of infiltration and runoff generation within these upland headwaters. This paper presents results from a set of rainfall simulation experiments performed on the blanket peat moorland of the North Pennines, UK. Rainfall was simulated at low intensities (3–12 mm h^?1), typical of natural rainfall, on bare and vegetated peat surfaces. Runoff response shows that infiltration rate increases with rainfall intensity; the use of low‐intensity rainfall therefore allows a more realistic evaluation of infiltration rates and flow processes than previous studies. Overland flow is shown to be common on both vegetated and bare peat surfaces although surface cover does exert some control. Most runoff is produced within the top few centimetres of the peat and runoff response decreases rapidly with depth. Little vertical percolation takes place to depths greater than 10 cm owing to the saturation of the peat mass. This study provides evidence that the quickflow response of upland blanket peat catchments is a result of saturation‐excess overland flow generation. Rainfall–runoff response from small plots varies with season. Following warm, dry weather, rainfall tends to infiltrate more readily into blanket peat, not just initially but to the extent that steady‐state surface runoff rates are reduced and more flow takes place within the peat, albeit at shallow depth. Sediment erosion from bare peat plots tends to be supply limited. Seasonal weather conditions may affect this in that after a warm, dry spell, surface desiccation allows sediment erosion to become transport limited. Copyright © 2002 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

One‐ and two‐dimensional modelling of overland flow in semiarid shrubland,Jornada basin,New Mexico

Two distributed parameter models, a one‐dimensional (1D) model and a two‐dimensional (2D) model, are developed to simulate overland flow in two small semiarid shrubland watersheds in the Jornada basin, southern New Mexico. The models are event‐based and represent each watershed by an array of 1‐m² cells, in which the cell size is approximately equal to the average area of the shrubs. Each model uses only six parameters, for which values are obtained from field surveys and rainfall simulation experiments. In the 1D model, flow volumes through a fixed network are computed by a simple finite‐difference solution to the 1D kinematic wave equation. In the 2D model, flow directions and volumes are computed by a second‐order predictor–corrector finite‐difference solution to the 2D kinematic wave equation, in which flow routing is implicit and may vary in response to flow conditions. The models are compared in terms of the runoff hydrograph and the spatial distribution of runoff. The simulation results suggest that both the 1D and the 2D models have much to offer as tools for the large‐scale study of overland flow. Because it is based on a fixed flow network, the 1D model is better suited to the study of runoff due to individual rainfall events, whereas the 2D model may, with further development, be used to study both runoff and erosion during multiple rainfall events in which the dynamic nature of the terrain becomes an important consideration. Copyright © 2006 John Wiley & Sons, Ltd.     

A GIS‐based variable source area hydrology model

Effective control of nonpoint source pollution from contaminants transported by runoff requires information about the source areas of surface runoff. Variable source hydrology is widely recognized by hydrologists, yet few methods exist for identifying the saturated areas that generate most runoff in humid regions. The Soil Moisture Routing model is a daily water balance model that simulates the hydrology for watersheds with shallow sloping soils. The model combines elevation, soil, and land use data within the geographic information system GRASS, and predicts the spatial distribution of soil moisture, evapotranspiration, saturation‐excess overland flow (i.e., surface runoff), and interflow throughout a watershed. The model was applied to a 170 hectare watershed in the Catskills region of New York State and observed stream flow hydrographs and soil moisture measurements were compared to model predictions. Stream flow prediction during non‐winter periods generally agreed with measured flow resulting in an average r² of 0·73, a standard error of 0·01 m³/s, and an average Nash‐Sutcliffe efficiency R² of 0·62. Soil moisture predictions showed trends similar to observations with errors on the order of the standard error of measurements. The model results were most accurate for non‐winter conditions. The model is currently used for making management decisions for reducing non‐point source pollution from manure spread fields in the Catskill watersheds which supply New York City's drinking water. Copyright © 1999 John Wiley & Sons, Ltd.     

Scale effects of Hortonian overland flow and rainfall–runoff dynamics in a West African catena landscape

Hortonian runoff was measured from plots with lengths of 1·25 and 12 m, and at watershed level for rainstorms during the 1996 rainy season in cental Côte d'Ivoire, Africa. A clear reduction in runoff coefficients was found with increasing slope lengths, giving order of magnitude differences between runoff measurements at point level (1 m²: 30–50% of total rain) and watershed level (130 ha: 4% of total rain). Runoff reduction from 1·25 and 12 m slopes was reproduced for each major runoff‐producing rainstorm at two different sets of plots, but the reduction was erratic for rainfall events which produced little runoff. In addition, runoff reduction varied wildly from one rainstorm to the next. In the analysis, we show that the spatial variability of runoff parameters causes the erratic behaviour during rainstorms with little runoff. During the more important, larger runoff‐producing events, which give 78% of total runoff, the temporal dynamics of the rainfall–runoff process determine the reduction of runoff coefficients from longer slopes. A simple infiltration/runoff model was used to simulate the field results, thereby confirming the importance of rainfall dynamics as an explanatory factor for measured reduction of runoff coefficients. Copyright © 2000 John Wiley & Sons, Ltd.     

Estimation of temporally averaged sediment delivery ratio using aggradational terraces in headwater catchments of the Waipaoa River,North Island,New Zealand

The sediment delivery ratio was estimated for two periods (28 years and eight years) following reforestation of seven tributary catchments (0·33 to 0·49 km²) in the headwaters of the Waipaoa River basin, North Island, New Zealand. In these catchments, gully erosion, which largely resulted from clearance of the natural forest between 1880 and 1920, is the main source of sediment to streams. Reforestation commenced in the early 1960s in an attempt to stabilize hillslopes and reduce sediment supply. Efforts have been partially successful and channels are now degrading, though gully erosion continues to supply sediment at accelerated rates in parts of the catchment. Data from the area indicate that the sediment delivery ratio (SDR) can be estimated as a function of two variables, ψ (the product of catchment area and channel slope) and A _g (the temporally averaged gully area for the period). Sediment input from gullies was determined from a well defined relationship between sediment yield and gully area. Sediment scoured from channels was estimated from dated terrace remnants and the current channel bed. Terrace remnants represent aggradation during major floods. This technique provides estimates of SDR averaged over periods between large magnitude terrace‐forming events and with the present channel bed. The technique averages out short‐term variability in sediment flux. Comparison of gully area and sediment transport between two periods (1960–1988 and 1988–1996) indicates that the annual rate of sediment yield from gullies for the later period has decreased by 77 per cent, sediment scouring in channels has increased by 124 per cent, and sediment delivered from catchments has decreased by 78 per cent. However, average SDR for the tributaries was found to be not significantly different between these periods. This may reflect the small number of catchments examined. It is also due to the fact that the volume of sediment scoured from channels was very small relative to that produced by gullies. According to the equation for SDR determined for the Waipaoa headwaters, SDR increases with increasing catchment area in the case where A _g and channel slope are fixed. This is because the amount of sediment produced from a channel by scouring increases with increasing catchment area. However, this relationship does not hold for the main stem of the study catchments, because sediment delivered from its tributaries still continues to accumulate in the channel. Higher order channels are, in effect, at a different stage in the aggradation/degradation cycle and it will take some time until a main channel reflects the effects of reforestation and its bed adjusts to net degradation. Results demonstrate significant differences among even low order catchments, and such differences will need to be taken into consideration when using SDR to estimate sediment yields. Copyright © 2001 John Wiley & Sons, Ltd.     

Spatial–temporal variability and hydrologic connectivity of runoff generation areas in a North Alabama pasture—implications for phosphorus transport

This study delineated spatially and temporally variable runoff generation areas in the Sand Mountain region pasture of North Alabama under natural rainfall conditions, and demonstrated that hydrologic connectivity is important for generating hillslope response when infiltration‐excess (IE) runoff mechanism dominates. Data from six rainfall events (13·7–32·3 mm) on an intensively instrumented pasture hillslope (0·12 ha) were analysed. Analysis of data from surface runoff sensors, tipping bucket rain gauge and HS‐flume demonstrated spatial and temporal variability in runoff generation areas. Results showed that the maximum runoff generation area, which contributed to runoff at the outlet of the hillslope, varied between 67 and 100%. Furthermore, because IE was the main runoff generation mechanism on the hillslope, the data showed that as the rainfall intensity changed during a rainfall event, the runoff generation areas expanded or contracted. During rainfall events with high‐intensity short‐ to medium‐duration, 4–8% of total rainfall was converted to runoff at the outlet. Rainfall events with medium‐ to low‐intensity, medium‐duration were found less likely to generate runoff at the outlet. In situ soil hydraulic conductivity (k) was measured across the hillslope, which confirmed its effect on hydrologic connectivity of runoff generation areas. Combined surface runoff sensor and k‐interpolated data clearly showed that during a rainfall event, lower k areas generate runoff first, and then, depending on rainfall intensity, runoff at the outlet is generated by hydrologically connected areas. It was concluded that in IE‐runoff‐dominated areas, rainfall intensity and k can explain hydrologic response. The study demonstrated that only connected areas of low k values generate surface runoff during high‐intensity rainfall events. Identification of these areas would serve as an important foundation for controlling nonpoint source pollutant transport, especially phosphorus. The best management practices can be developed and implemented to reduce transport of phosphorus from these hydrologically connected areas. Copyright © 2009 John Wiley & Sons, Ltd.     

Spatial and temporal variability of groundwater dynamics in a sub‐Mediterranean mountain catchment

The temporal and spatial dynamics of groundwater was investigated in a small catchment in the Spanish Pyrenees, which was extensively used for agriculture in the past. Analysis of the water table fluctuations at five locations over a 6‐year period demonstrated that the groundwater dynamics had a marked seasonal cycle involving a wetting‐up period that commenced with the first autumn rainfall events, a saturation period during winter and spring and a drying‐down period from the end of spring until the end of the summer. The length of the saturation period showed great interannual variability, which was mainly influenced by the rainfall and evapotranspiration characteristics. There was marked spatial variability in the water table, especially during the wetting‐up period, which could be related to differences in slope and drainage area, geomorphology, soil properties and local topography. Areas contributing to runoff generation were identified within the catchment by field mapping of moisture conditions. Areas contributing to infiltration excess runoff were correlated with former cultivated fields affected by severe sheetwash erosion. Areas contributing to saturation excess runoff were characterized by a marked spatial dynamics associated with catchment wetness conditions. The saturation spatial pattern, which was partially related to the topographic index, was very patchy throughout the catchment, suggesting the influence of other factors associated with past agricultural activities, including changes in local topography and soil properties. The relationship between water table levels and stream flow was weak, especially during the wetting‐up period, suggesting little connection between ground water and the hydrological response, at least at some locations. The results suggest that in drier and human‐disturbed environments, such as sub‐Mediterranean mountains, saturation patterns cannot be represented only by the general topography of the catchment. They also suggest that groundwater storage and runoff is not a succession of steady‐state flow conditions, as assumed in most hydrological models. Copyright © 2013 John Wiley & Sons, Ltd.     

Flow behaviour,suspended sediment transport and transmission losses in a small (sub‐bank‐full) flow event in an Australian desert stream

The behaviour of a discrete sub‐bank‐full flow event in a small desert stream in western NSW, Australia, is analysed from direct observation and sediment sampling during the flow event and from later channel surveys. The flow event, the result of an isolated afternoon thunderstorm, had a peak discharge of 9 m³/s at an upstream station. Transmission loss totally consumed the flow over the following 7·6 km. Suspended sediment concentration was highest at the flow front (not the discharge peak) and declined linearly with the log of time since passage of the flow front, regardless of discharge variation. The transmission loss responsible for the waning and eventual cessation of flow occurred at a mean rate of 13.2% per km. This is quite rapid, and is more than twice the corresponding figure for bank‐full flows estimated by Dunkerley (1992) on the same stream system. It is proposed that transmission losses in ephemeral streams of the kind studied may be minimized in flows near bank‐full stage, and be higher in both sub‐bank‐full and overbank flows. Factors contributing to enhanced flow loss in the sub‐bank‐full flow studied included abstractions of flow to pools, scour holes and other low points along the channel, and overflow abstractions into channel filaments that did not rejoin the main flow. On the other hand, losses were curtailed by the shallow depth of banks wetted and by extensive mud drapes that were set down over sand bars and other porous channel materials during the flow. Thus, in contrast with the relatively regular pattern of transmission loss inferred from large floods, losses from low flows exhibit marked spatial variability and depend to a considerable extent on streamwise variations in channel geometry, in addition to the depth and porosity of channel perimeter sediments. Copyright © 1999 John Wiley & Sons, Ltd.     

Solute concentrations of overland flow water in a cultivated field: spatial variations,intra‐ and inter‐storm trends

Major solute concentrations in overland flow water (OFW) were measured in an agricultural field of Brittany (western France). Two storm events were monitored in detail to examine the short time‐scale processes. During one year, samples were taken at different positions on the slope after each storm event to describe the spatial and seasonal variations of OFW chemistry. Although the total dissolved load in OFW is not much higher than in rain water, distinctive features are observed. K⁺, Ca²⁺, NH₄ , Cl⁻ and SOare the major solutes. The main origin of the elements (sea salts, exchangeable soil complex or fertilizers) determined most of the variations observed. Spatial variations along the slope are mainly seen for exchangeable cations, while seasonal variations are predominant for sea salts. Rainfall intensity and suspended sediment load induce strong differences between the two storm events studied in detail. However, the within‐storm variations and the seasonal monitoring show that this relationship is complex. Within‐storm variations suggest that, in addition to desorption processes, mixing with pre‐event water may occur. The lack of a relationship between sediment load and dissolved load is attributed to the high rate of the exchange processes, which has been checked by a simple experiment in vitro. It is concluded that the conditions of the transit of water on the field (velocity, length, status of the surface, crusted or not) may well play a major role in the chemical changes between rain water and OFW. The results suggest that vegetated buffer strips designed to reduce the sediment load only, and not the amount of overland flow, will have little effect on the transfer of dissolved pollutants to the watercourses. Copyright © 1999 John Wiley & Sons, Ltd.     

On the assessment of the management priority of sediment source areas in a debris‐flow catchment

In many Alpine catchments the monitoring and hazard mitigation of debris‐flow events require enormous economic and social resources. To confront these problems, a subjective estimation of the most hazardous zones of the basin could be useful in the best, sustainable planning of protective measures. In this paper, a new methodology is proposed that develops a Management Priority Index (MPI) to rank sediment source areas by their quantitative capability to deliver debris‐flow volumes to a point of interest within the catchment. The MPI sets the intervention priority based on a combination of three sub‐indicators: a susceptibility indicator evaluating the overall catchment predisposition to generate debris flow, a triggering indicator and a volume budget indicator assessing the rate of deliverable volume to a selected outlet. MPI was applied to the basin of the Rio Gadria catchment (Venosta Valley, Bolzano, Italy), an alpine basin with an unlimited sediment supply that is characterized by multiple, very active, shallow landslides and bare soil zones. The proposed ranking method was successfully verified using post‐event surveys and through evidence from consolidation check dams built over many years in the basin. Copyright © 2014 John Wiley & Sons, Ltd.     

Improving the spatial representation of soil properties and hydrology using topographically derived initialization processes in the SWAT model

Topography exerts critical controls on many hydrologic, geomorphologic and biophysical processes. However, many watershed modelling systems use topographic data only to define basin boundaries and stream channels, neglecting opportunities to account for topographic controls on processes such as soil genesis, soil moisture distributions and hydrological response. Here, we demonstrate a method that uses topographic data to adjust spatial soil morphologic and hydrologic attributes: texture, depth to the C‐horizon, saturated conductivity, bulk density, porosity and the water capacities at field (33 kpa) and wilting point (1500 kpa) tensions. As a proof of concept and initial performance test, the values of the topographically adjusted soil parameters and those from the Soil Survey Geographic Database (SSURGO; available at 1 : 20 000 scale) were compared with measured soil pedon pit data in the Grasslands Soil and Water Research Lab watershed in Riesel, TX. The topographically adjusted soils were better correlated with the pit measurements than were the SSURGO values. We then incorporated the topographically adjusted soils into an initialization of the Soil and Water Assessment Tool model for 15 Riesel research watersheds to investigate how changes in soil properties influence modelled hydrological responses at the field scale. The results showed that the topographically adjusted soils produced better runoff predictions in 50% of the fields, with the SSURGO soils performing better in the remainder. In addition, the a priori adjusted soils result in fewer calibrated model parameters. These results indicate that adjusting soil properties based on topography can result in more accurate soil characterization and, in some cases, improve model performance. Copyright © 2016 John Wiley & Sons, Ltd.     

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Relation of sediment transport capacity to stone cover and size in rain‐impacted interrill flow

This study examines how the sediment transport capacity of interrill overland flow varies with stone cover and stone size at two flow intensities. Six series of flume experiments were conducted on two slopes (2° and 10°) with stones of three sizes (28·0, 45·5 and 91·3 mm) serving as roughness elements. Bed sediment size, water discharge and simulated rainfall intensity were the same in all experiments. It was found (1) that transport capacity is positively related to stone size, with the relation becoming stronger as stone cover increases and flow intensity decreases; and (2) that transport capacity is negatively related to stone cover at the high flow intensity and curvilinearly related to stone cover at the low flow intensity. The curvilinear relations are concave‐upward with the lowest transport capacities occurring at stone covers between 0·40 and 0·60. The highest transport capacities are found at stone covers of 0 and 1, with the transport capacity being greater at the former stone cover than at the latter. 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.     

Towards a simple dynamic process conceptualization in rainfall–runoff models using multi‐criteria calibration and tracers in temperate,upland catchments

Empirically based understanding of streamflow generation dynamics in a montane headwater catchment formed the basis for the development of simple, low‐parameterized, rainfall–runoff models. This study was based in the Girnock catchment in the Cairngorm Mountains of Scotland, where runoff generation is dominated by overland flow from peaty soils in valley bottom areas that are characterized by dynamic expansion and contraction of saturation zones. A stepwise procedure was used to select the level of model complexity that could be supported by field data. This facilitated the assessment of the way the dynamic process representation improved model performance. Model performance was evaluated using a multi‐criteria calibration procedure which applied a time series of hydrochemical tracers as an additional objective function. Flow simulations comparing a static against the dynamic saturation area model (SAM) substantially improved several evaluation criteria. Multi‐criteria evaluation using ensembles of performance measures provided a much more comprehensive assessment of the model performance than single efficiency statistics, which alone, could be misleading. Simulation of conservative source area tracers (Gran alkalinity) as part of the calibration procedure showed that a simple two‐storage model is the minimum complexity needed to capture the dominant processes governing catchment response. Additionally, calibration was improved by the integration of tracers into the flow model, which constrained model uncertainty and improved the hydrodynamics of simulations in a way that plausibly captured the contribution of different source areas to streamflow. This approach contributes to the quest for low‐parameter models that can achieve process‐based simulation of hydrological response. Copyright © 2009 John Wiley & Sons, Ltd.     

Influence of soil‐surface types on the overall runoff of the Tabernas badlands (south‐east Spain): field data and model approaches

The Tabernas desert, an extensive badlands area in Almeria province (south‐east Spain), is characterized by a high variability in soil surface cover and soil properties along with important topographical contrasts giving rise to a wide range of hydrological behaviour. A double approach through field monitoring and modelling has been used to ascertain the influence of soil‐surface variability on the overall hydrological response. Small plots were monitored for 3 years to assess runoff from the different surface types. Data provided by the long‐term monitoring of three small catchments formed by different soil surfaces were used to find out the specific contribution of each soil surface to the catchment runoff. A simple spatially distributed model was built to predict runoff generation based on the infiltration rate of each soil‐surface type (defined as terrain units with the same cover, the same soil type and on the same landform). Plot results prove that the soil surface units within the study area behave differently in terms of hydrological response to natural rainfall. These responses are explained by the types of cover, topographical characteristics and soil properties. When runoff events are simple (with one or two runoff peaks), the modelled hydrographs reproduce the hydrographs observed reasonably well, but in complex events (with several runoff peaks) the adjustment is not as good. The model also shows the influence of the spatial distribution of soil surfaces on the overall runoff, aiding exploration of the spatial hydrological relationships among different landscape units. Copyright © 2002 John Wiley & Sons, Ltd.