How does landscape structure influence catchment transit time across different geomorphic provinces?

Despite an increasing number of empirical investigations of catchment transit times (TTs), virtually all are based on individual catchments and there are few attempts to synthesize understanding across different geographical regions. Uniquely, this paper examines data from 55 catchments in five geomorphic provinces in northern temperate regions (Scotland, United States of America and Sweden). The objective is to understand how the role of catchment topography as a control on the TTs differs in contrasting geographical settings. Catchment inverse transit time proxies (ITTPs) were inferred by a simple metric of isotopic tracer damping, using the ratio of standard deviation of δ¹⁸O in streamwater to the standard deviation of δ¹⁸O in precipitation. Quantitative landscape analysis was undertaken to characterize the catchments according to hydrologically relevant topographic indices that could be readily determined from a digital terrain model (DTM). The nature of topographic controls on transit times varied markedly in different geomorphic regions. In steeper montane regions, there are stronger gravitational influences on hydraulic gradients and TTs tend to be lower in the steepest catchments. In provinces where terrain is more subdued, direct topographic control weakened; in particular, where flatter areas with less permeable soils give rise to overland flow and lower TTs. The steeper slopes within this flatter terrain appear to have a greater coverage of freely draining soils, which increase sub‐surface flow, therefore increasing TTs. Quantitative landscape analysis proved a useful tool for inter‐catchment comparison. However, the critical influence of sub‐surface permeability and connectivity may limit the transferability of predictive tools of hydrological function based on topographic parameters alone. Copyright © 2009 John Wiley & Sons, Ltd.     

Conceptualization of runoff processes using a geographical information system and tracers in a nested mesoscale catchment

Tracer investigations were combined with a geographical information system (GIS) analysis of the 31 km² Girnock catchment (Cairngorm Mountains, Scotland) in order to understand hydrological functioning by identifying dominant runoff sources and estimating mean residence times. The catchment has a complex geology, soil cover and topography. Gran alkalinity was used to demonstrate that catchment geology has a dominant influence on baseflow chemistry, but flow paths originating in acidic horizons in the upper soil profiles controlled stormflow alkalinity. Chemically based hydrograph separations at the catchment scale indicated that ～30% of annual runoff was derived from groundwater sources. Similar contributions (23–36%) were estimated for virtually all major sub‐basins. δ¹⁸O of precipitation (mean: ? 9·4‰; range: ? 16·1 to ? 5·0‰) and stream waters (mean: ? 9·1‰; range: ? 11·6 to ? 7·4‰) were used to assess mean catchment and sub‐basin residence times, which were in the order ～4–6 months. GIS analysis showed that these tracer‐based diagnostic features of catchment functioning were consistent with the landscape organization of the catchment. Soil and HOST (Hydrology of Soil Type) maps indicated that the catchment and individual sub‐basins were dominated by hydrologically responsive soils, such as peats (Histosol), peaty gleys (Histic Gleysols) and rankers (Umbric Leptosols and Histosols). Soil cover (in combination with a topographic index) predicted extensive areas of saturation that probably expand during hydrological events, thus providing a high degree of hydrological connectivity between catchment hillslopes and stream channel network. This was validated by aerial photographic interpretation and groundtruthing. These characteristics of hydrological functioning (i.e. dominance of responsive hydrological pathways and short residence times) dictate that the catchment is sensitive to land use change impacts on the quality and quantity of streamflows. It is suggested that such conceptualization of hydrological functioning using tracer‐validated GIS analysis can play an important role in the sustainable management of river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

Is groundwater response timing in a pre‐alpine catchment controlled more by topography or by rainfall?

Groundwater levels in steep headwater catchments typically respond quickly to rainfall, but the timing of the response may vary spatially across the catchment. In this study, we investigated the topographic controls and the effects of rainfall and antecedent conditions on the groundwater response timing for 51 groundwater monitoring sites in a 20‐ha pre‐alpine catchment with low permeability soils. The median time to rise and median duration of recession for the 133 rainfall events were highly correlated to the topographic characteristics of the site and its upslope contributing area. The median time to rise depended more on the topographic characteristics than on the rainfall characteristics or antecedent soil wetness conditions. The median time to rise decreased with Topographic Wetness Index (TWI) for sites with TWI < 6 and was almost constant for sites with a higher TWI. The slope of this relation was a function of rainfall intensity. The rainfall threshold for groundwater initiation was also a function of TWI and allowed extrapolation of point measurements to the catchment scale. The median lag time between the rainfall centroid and the groundwater peak was 75 min. The groundwater level peaked before peak streamflow at the catchment outlet for half of the groundwater monitoring sites, but only by 15 to 25 min. The stronger correlations between topographic indices and groundwater response timing in this study compared to previous studies suggest that surface topography affects the groundwater response timing in catchments with low permeability soils more than in catchments with more transmissive soils. Copyright © 2015 John Wiley & Sons, Ltd.     

Catchment transit times and landscape controls—does scale matter?

Mean transit times (MTTs) can give useful insights into the internal processes of hydrological systems. However, our understanding of how they vary and scale remains unclear. We used MTT estimates obtained from δ¹⁸O data from 20, mostly nested, contrasting catchments in North East Scotland, ranging from 1 to 1700 km². The estimated MTTs ranged between 270 and 1170 days and were used to test a previously developed multiple linear regression (MLR) model for MTT prediction based on metrics of soil cover, landscape organization and climate. We show that the controls on MTT identified by the MLR model hold with the independent data from these 20 sites and that the MLR can be used to predict MTT in ungauged montane catchments. The dominant controls also remain unchanged over four orders of magnitude of catchment size, suggesting no major change of dominant flow paths and mixing processes at larger scales. This is consistent with the fact that only the variance of MTT, rather than MTT, showed a scaling relationship. MTTs appeared to converge with increasing catchment scale, apparently due to the integration of heterogeneous headwater responses in larger downstream catchments. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimating runoff from a glacierized catchment using natural tracers in the semi‐arid Andes cordillera

This paper presents a methodology for hydrograph separation in mountain watersheds, which aims at identifying flow sources among ungauged headwater sub‐catchments through a combination of observed streamflow and data on natural tracers including isotope and dissolved solids. Daily summer and bi‐daily spring season water samples obtained at the outlet of the Juncal River Basin in the Andes of Central Chile were analysed for all major ions as well as stable water isotopes, δ¹⁸O and δD. Additionally, various samples from rain, snow, surface streams and exfiltrating subsurface water (springs) were sampled throughout the catchment. A principal component analysis was performed in order to address cross‐correlation in the tracer dataset, reduce the dimensionality of the problem and uncover patterns of variability. Potential sources were identified in a two‐component U‐space that explains 94% of the observed tracer variability at the catchment outlet. Hydrograph separation was performed through an Informative‐Bayesian model. Our results indicate that the Juncal Norte Glacier headwater sub‐catchment contributed at least 50% of summer flows at the Juncal River Basin outlet during the 2011–2012 water year (a hydrologically dry period in the Region), even though it accounts for only 27% of the basin area. Our study confirms the value of combining solute and isotope information for estimating source contributions in complex hydrologic systems, and provides insights regarding experimental design in high‐elevation semi‐arid catchments. The findings of this study can be useful for evaluating modelling studies of the hydrological consequences of the rapid decrease in glacier cover observed in this region, by providing insights into the origin of river water in basins with little hydrometeorological information. Copyright © 2016 John Wiley & Sons, Ltd.     

Stormflow generation and flowpath characteristics in an Amazonian rainforest catchment

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

Isotopic and geochemical tracers reveal similarities in transit times in contrasting mesoscale catchments

There is a need for more isotopic tracer studies at the mesoscale to extend our understanding of catchment transit times and their associated controls beyond smaller experimental sites (typically < 10 km²). This paper, therefore, examines the isotope hydrology of six mesoscale (10¹–10² km²) sub‐catchments of the 2000 km² basin of the River Dee in northern Scotland. All the catchments were upland in character (mean altitude > 400 m) with similar suites of soil coverage (predominantly regosols, gleys, peats and podzols), although the relative distribution varied, as did the presence of other landscape features such as aquifers in Quaternary drifts and lakes. Input–output relationships of δ¹⁸O in precipitation and runoff revealed contrasting responses and differential damping which were broadly consistent with catchment characteristics. The mean transit times (MTTs) were estimated using a convolution integral with a Gamma distribution as the transfer function. These varied from 528 days in the most responsive catchments to > 800 days in catchments where the tracer signature was most damped. Shorter MTTs were found in sub‐catchments with a higher percentage cover of responsive soils (i.e. regosols, gleys and peats), whilst sub‐catchments with longest MTTs had a higher coverage of free‐draining podzolic and alluvial soils, as well as significant amount of stored water either in fluvio‐glacial aquifers or large lakes. The MTT of all six catchments had the same order of magnitude; this contrasts with studies in the Scottish Highlands with smaller (<10 km²) catchments where MTT has been shown to vary between 60 and 1200 days. Copyright © 2010 John Wiley & Sons, Ltd.     

Topographic,pedologic and climatic interactions influencing streamflow generation at multiple catchment scales

Dominant flow pathways (DFPs) in mesoscale watersheds are poorly characterized and understood. Here, we make use of a conservative tracer (Gran alkalinity) and detailed information about climatic conditions and physical properties to examine how temporally and spatially variable factors interact to determine DFPs in 12 catchments draining areas from 3.4 to 1829.5 km² (Cairngorms, Scotland). After end‐member mixing was applied to discriminate between near surface and deep groundwater flow pathways, variation partitioning, canonical redundancy analyses and regression models were used to resolve: (i) What is the temporal variability of DFPs in each catchment?; (ii) How do DFPs change across spatial scales and what factors control the differences in hydrological responses?; and (iii) Can a conceptual model be developed to explain the spatiotemporal variability of DFPs as a function of climatic, topographic and soil characteristics? Overall, catchment characteristics were only useful to explain the temporal variability of DFPs but not their spatial variation across scale. The temporal variability of DFPs was influenced most by prevailing hydroclimatic conditions and secondarily soil drainability. The predictability of active DFPs was better in catchments with soils supporting fast runoff generation on the basis of factors such as the cumulative precipitation from the seven previous days, mean daily air temperature and the fractional area covered by Rankers. The best regression model R² was 0.54, thus suggesting that the catchments’ internal complexity was not fully captured by the factors included in the analysis. Nevertheless, this study highlights the utility of combining tracer studies with digital landscape analysis and multivariate statistical techniques to gain insights into the temporal (climatic) and spatial (topographic and pedologic) controls on DFPs. Copyright © 2011 John Wiley & Sons, Ltd.     

Projecting climate change impacts on stream flow regimes with tracer‐aided runoff models ‐ preliminary assessment of heterogeneity at the mesoscale

The northern mid‐high latitudes form a region that is sensitive to climate change, and many areas already have seen – or are projected to see – marked changes in hydroclimatic drivers on catchment hydrological function. In this paper, we use tracer‐aided conceptual runoff models to investigate such impacts in a mesoscale (749 km²) catchment in northern Scotland. The catchment encompasses both sub‐arctic montane sub‐catchments with high precipitation and significant snow influence and drier, warmer lowland sub‐catchments. We used downscaled HadCM3 General Circulation Model outputs through the UKCP09 stochastic weather generator to project the future climate. This was based on synthetic precipitation and temperature time series generated from three climate change scenarios under low, medium and high greenhouse gas emissions. Within an uncertainty framework, we examined the impact of climate change at the monthly, seasonal and annual scales and projected impacts on flow regimes in upland and lowland sub‐catchments using hydrological models with appropriate process conceptualization for each landscape unit. The results reveal landscape‐specific sensitivity to climate change. In the uplands, higher temperatures result in diminishing snow influence which increases winter flows, with a concomitant decline in spring flows as melt reduces. In the lowlands, increases in air temperatures and re‐distribution of precipitation towards autumn and winter lead to strongly reduced summer flows despite increasing annual precipitation. The integration at the catchment outlet moderates these seasonal extremes expected in the headwaters. This highlights the intimate connection between hydrological dynamics and catchment characteristics which reflect landscape evolution. It also indicates that spatial variability of changes in climatic forcing combined with differential landscape sensitivity in large heterogeneous catchments can lead to higher resilience of the integrated runoff response. Copyright © 2012 John Wiley & Sons, Ltd.     

On the relationships between catchment scale and streamwater mean residence time

The relationship between streamwater mean residence time (MRT) and landscape characteristics is poorly understood. We used tritium (³H) to define our MRT. We tested the hypothesis that baseflow water MRT increases with increasing absolute catchment size at the Maimai catchments. These catchments are simple hydrologic systems relative to many catchments around the world, with uniformly wet climatic conditions, little seasonality, uniform and nearly impermeable bedrock, steep short hillslopes, shallow soils, and well‐characterized hillslope and catchment hydrology. As a result, this is a relatively simple system and an ideal location for new MRT‐related hypothesis testing. Whilst hydrologists have used ³H to estimate water age since the 1960s nuclear testing spike, atmospheric ³H levels have now approached near background levels and are often complicated by contamination from the nuclear industry. We present results for ³H sampled from our set of nested catchments in nuclear‐industry‐free New Zealand. Because of high precision analysis, near‐natural atmospheric ³H levels, and well‐characterized rainfall ³H inputs, we were able to estimate the age of young (i.e. less than 3 years old) waters. Our results showed no correlation between MRT and catchment size. However, MRT was correlated to the median sub‐catchment size of the sampled watersheds, as shown by landscape analysis of catchment area accumulation patterns. These preliminary findings suggest that landscape organization, rather than total area, is a first‐order control on MRT and points the way forward for more detailed analysis of how landscape organization affects catchment runoff characteristics. Copyright © 2003 John Wiley & Sons, Ltd.     

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

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

Testing of the SIBERIA landscape evolution model using the Tin Camp Creek,Northern Territory,Australia, field catchment

The SIBERIA landscape evolution model was used to simulate the geomorphic development of the Tin Camp Creek natural catchment over geological time. Measured hydrology, erosion and geomorphic data were used to calibrate the SIBERIA model, which was then used to make independent predictions of the landform geomorphology of the study site. The catchment, located in the Northern Territory, Australia is relatively untouched by Europeans so the hydrological and erosion processes that shaped the area can be assumed to be the same today as they have been in the past, subject to the caveats regarding long‐term climate fluctuation. A qualitative, or visual comparison between the natural and simulated catchments indicates that SIBERIA can match hillslope length and hillslope profile of the natural catchments. A comparison of geomorphic and hydrological statistics such as the hypsometric curve, width function, cumulative area distribution and area–slope relationship indicates that SIBERIA can model the geomorphology of the selected Tin Camp Creek catchments. Copyright 2002 © Environmental Research Institute of the Supervising Scientist, Commonwealth of Australia.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

The importance of aspect for modelling the hydrological response in a glacier catchment in Central Asia

Understanding how explicit consideration of topographic information influences hydrological model performance and upscaling in glacier dominated catchments remains underexplored. In this study, the Urumqi glacier no. 1 catchment in northwest China, with 52% of the area covered by glaciers, was selected as study site. A conceptual glacier‐hydrological model was developed and tested to systematically, simultaneously, and robustly reproduce the hydrograph, separate the discharge into contributions from glacier and nonglacier parts of the catchment, and establish estimates of the annual glacier mass balance, the annual equilibrium line altitude, and the daily catchment snow water equivalent. This was done by extending and adapting a recently proposed landscape‐based semidistributed conceptual hydrological model (FLEX‐Topo) to represent glacier and snowmelt processes. The adapted model, FLEX^G, allows to explicitly account for the influence of topography, that is, elevation and aspect, on the distribution of temperature and precipitation and thus on melt dynamics. It is shown that the model can not only reproduce long‐term runoff observations but also variations in glacier and snow cover. Furthermore, FLEX^G was successfully transferred and up‐scaled to a larger catchment exclusively by adjusting the areal proportions of elevation and aspect without the need for further calibration. This underlines the value of topographic information to meaningfully represent the dominant hydrological processes in the region and is further exacerbated by comparing the model to a model formulation that does not account for differences in aspect (FLEX^G,nA) and which, in spite of satisfactorily reproducing the observed hydrograph, does not capture the influence of spatial variability of snow and ice, which as a consequence reduces model transferability. This highlights the importance of accounting for topography and landscape heterogeneity in conceptual hydrological models in mountainous and snow‐, and glacier‐dominated regions.     

Use of a parsimonious rainfall–run‐off model for predicting hydrological response in ungauged basins

A new parameter parsimonious rainfall–run‐off model, the Distance Distribution Dynamics (DDD) model, is used to simulate hydrological time series at ungauged sites in the Lygne basin in Norway. The model parameters were estimated as functions of catchment characteristics determined by geographical information system. The multiple regression equations relating catchment characteristics and model parameters were trained from 84 calibrated catchments located all over Norway, and all model parameters showed significant correlations with catchment characteristics. The significant correlation coefficients (with p‐value < 0.05) ranged from 0.22 to 0.55. The suitability of DDD for predictions in ungauged basins was tested for 17 catchments not used to estimate the multiple regression equations. For ten of the 17 catchments, deviations in Nash–Sutcliffe efficiency (NSE) criteria between the calibrated and regionalised model were less than 0.1, and for two calibrated catchments within the Lygne basin, the deviations were less than 0.08. The median NSE for the regionalized DDD for the 17 catchments for two time series was 0.66 and 0.72. Deviations in NSE between calibrated and regionalised models are well explained by the deviations between calibrated and regressed parameters describing spatial snow distribution and snowmelt respectively. The quality of the simulated run‐off series for the ungauged sites in the Lygne basin was assessed by comparing flow indices describing high, medium and low flow estimated from observed run‐off at the 17 catchments and for the simulated run‐off series. The indices estimated for the simulated series were generally well within the ranges defined by the 17 observed series. Copyright © 2014 John Wiley & Sons, Ltd.     

A broad‐scale assessment of the effect of improved soil management on catchment baseflow index

Although the importance of sustainable soil management is recognized, there are many threats to soils including widespread soil structural degradation. This reduces infiltration through the soil surface and/or the percolation of water through the soil profile, with important consequences for crop yields, nutrient cycling and the hydrological response of catchments. This article describes a broad‐scale modelling approach to assess the potential effect that improved agricultural soil management, through reduced soil structural degradation, may have on the baseflow index (BFI) of catchments across England and Wales. A daily soil–water balance model was used to simulate the indicative BFI of 45 696 thirty‐year model runs for different combinations of soil type, soil/field condition, land cover class and climate which encapsulate the variability across England and Wales. The indicative BFI of catchments was then calculated by upscaling the results by spatial weighting. WaSim model outputs of indicative BFI were within the 95% confidence intervals of the national‐average BFI values given for the Hydrology of Soil Type (HOST ? ) classes for 26 of the 28 classes. At the catchment scale, the concordance correlation coefficient between the BFI from the WaSim model outputs and those derived from HOST was 0·83. Plausible improvements in agricultural soil/field condition produced modest simulated increases of up to 10% in the indicative BFI in most catchments across England and Wales, although for much of southern and northern England the increases were less than 5%. The results suggest that improved soil management might partially mitigate the expected adverse effects of climate change on baseflow to rivers. Healthy, well‐functioning soils produce many additional benefits such as better agricultural yields and reduced pollutant movement, so improved soil management should provide win‐win opportunities for society, agricultural systems and the environment and provide resilience to some of the expected environmental impacts of climate change. Copyright © 2011 John Wiley & Sons, Ltd.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Nonparametric catchment clustering using the data depth function

ABSTRACT

The clustering of catchments is important for prediction in ungauged basins, model parameterization and watershed development and management. The aim of this study is to explore a new measure of similarity among catchments, using a data depth function and comparing it with catchment clustering indices based on flow and physical characteristics. A cluster analysis was performed for each similarity measure using the affinity propagation clustering algorithm. We evaluated the similarity measure based on depth–depth plots (DD-plots) as a basis for transferring parameter sets of a hydrological model between catchments. A case study was developed with 21 catchments in a diverse New Zealand region. Results show that clustering based on the depth–depth measure is dissimilar to clustering on catchment characteristics, flow, or flow indices. A hydrological model was calibrated for the 21 catchments and the transferability of model parameters among similar catchments was tested within and between clusters defined by each clustering method. The mean model performance for parameters transferred within a group always outperformed those from outside the group. The DD-plot based method was found to produce the best in-group performance and second-highest difference between in-group and out-group performance.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR A. Viglione     

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.     

Impact of grid resolution on the integrated and distributed response of a coupled surface–subsurface hydrological model for the des Anglais catchment,Quebec

Digital elevation models (DEMs) at different resolutions (180, 360, and 720 m) are used to examine the impact of different levels of landscape representation on the hydrological response of a 690‐km² catchment in southern Quebec. Frequency distributions of local slope, plan curvature, and drainage area are calculated for each grid size resolution. This landscape analysis reveals that DEM grid size significantly affects computed topographic attributes, which in turn explains some of the differences in the hydrological simulations. The simulations that are then carried out, using a coupled, process‐based model of surface and subsurface flow, examine the effects of grid size on both the integrated response of the catchment (discharge at the main outlet and at two internal points) and the distributed response (water table depth, surface saturation, and soil water storage). The results indicate that discharge volumes increase as the DEM is coarsened, and that coarser DEMs are also wetter overall in terms of water table depth and soil water storage. The reasons for these trends include an increase in the total drainage area of the catchment for larger DEM cell sizes, due to aggregation effects at the boundary cells of the catchment, and to a decrease in local slope and plan curvature variations, which in turn limits the capacity of the watershed to transmit water downslope and laterally. The results obtained also show that grid resolution effects are less pronounced during dry periods when soil moisture dynamics are mostly controlled by vertical fluxes of evaporation and percolation. Copyright © 2010 John Wiley & Sons, Ltd.