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.     

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.     

Gridded climate data products are an alternative to instrumental measurements as inputs to rainfall–runoff models

Rainfall–runoff models are widely used to predict flows using observed (instrumental) time series of air temperature and precipitation as inputs. Poor model performance is often associated with difficulties in estimating catchment‐scale meteorological variables from point observations. Readily available gridded climate products are an underutilized source of temperature and precipitation time series for rainfall–runoff modelling, which may overcome some of the performance issues associated with poor‐quality instrumental data in small headwater monitoring catchments. Here we compare the performance of instrumental measured and E‐OBS gridded temperature and precipitation time series as inputs in the rainfall–runoff models “PERSiST” and “HBV” for flow prediction in six small Swedish catchments. For both models and most catchments, the gridded data produced statistically better simulations than did those obtained using instrumental measurements. Despite the high correspondence between instrumental and gridded temperature, both temperature and precipitation were responsible for the difference. We conclude that (a) gridded climate products such as the E‐OBS dataset could be more widely used as alternative input to rainfall–runoff models, even when instrumental measurements are available, and (b) the processing applied to gridded climate products appears to provide a more realistic approximation of small catchment‐scale temperature and precipitation patterns needed for flow simulations. Further research on this issue is needed and encouraged.     

Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

MODELLING RUNOFF GENERATION ON SMALL AGRICULTURAL CATCHMENTS: CAN REAL WORLD RUNOFF RESPONSES BE CAPTURED?

This paper focuses on the problem of quantifying real world catchment response using a distributed model and discusses the ability of the model to capture that response. The rainfall–runoff responses of seven small agricultural catchments in the eastern wheatbelt region of south-western Australia are examined. The variability in runoff generation and the factors that contribute to that variability (i.e. rainfall intensity, soil properties and topography) are investigated to determine if their influence can be captured in a mathematical model. The spatially distributed rainfall–runoff model used in this study is based on the TOPMODEL concepts of Beven and Kirkby (1979), and simulates runoff generation by both the infiltration excess and saturation excess mechanisms. Simulations with the model revealed the highly complex nature of catchment response to rainfall events. Runoff generation was highly heterogeneous in both space and time, with the runoff response being governed by the spatial variability of soil properties and topography, and by the temporal variation in rainfall intensity. Although the model proved capable of simulating catchment response for many events, the investigation has demonstrated that not all aspects of the variability associated with agricultural catchments (particularly the effects of land management) can be captured using this relatively simple model. © 1997 by John Wiley & Sons, Ltd     

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

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

Simulated impacts of climate change and land‐clearing on runoff from a small Sahelian catchment

In the Sahel, there are few long‐term data series available to estimate the climatic and anthropogenic impacts on runoff in small catchments. Since 1950, land clearing has enhanced runoff. The question is whether and by how much this anthropogenic effect offsets the current drought. To answer this question, a physically based distributed hydrological model was used to simulate runoff in a small Sahelian catchment in Niger, from the 1950–1998 rain‐series. The simulation was carried out for three soil surface states of the catchment (1950, 1975 and 1992). The catchment is characterized by an increase in cultivated land, with associated fallow, from 6% in 1950 to 56% in 1992, together with an increase in the extent of eroded land (from 7 to 16%), at the expense of the savanna. Effects of climate and land use are first analysed separately: irrespective of the land cover state, the simulated mean annual runoff decreases by about 40% from the wet period (1950–1969) to the dry period (1970–1998); calculated on the 1950–1998 rainfall‐series, the changes that occurred in land cover between 1950 and 1992 multiplies the mean annual runoff by a factor close to three. The analysis of a joint climatic and anthropogenic change shows that the transition from a wet period under a ‘natural’ land cover (1950) to a dry period under a cultivated land cover (1992) results in an increase in runoff of the order of 30 to 70%. At the scale of a small Sahelian catchment, the anthropogenic impact on runoff is probably more important than that of drought. This figure for relative increase in runoff contributions to ponds, preferential sites of seepage to groundwater, is less than that currently estimated for aquifer recharge, which has been causing a significant continuous water table rise over the same period. Copyright © 2004 John Wiley & Sons, Ltd.     

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

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

Advancing understanding of runoff and sediment transfers in agricultural catchments through simultaneous observations across scales

Our understanding of the effect of scale on runoff and sediment transfers within catchments is currently limited by a lack of available data. A multi‐scale dataset of 17 rainfall events collected simultaneously at four spatial scales within a small agricultural catchment in 2005–2006 is presented. Analysis using exploratory techniques and a two‐step, zero‐inflated lognormal mixed‐effects regression model, has demonstrated that event responses, and event response characteristics representing runoff and sediment peaks and area‐normalized yields, are scale dependent, and hence cannot be transferred directly between scales. Runoff and sediment yields increase as scale increases, and it is proposed that this effect, which differs from that observed in the few other studies of scale effects undertaken, is due to increasing connectivity within the catchment, and the dominance of preferential flow pathways including through macropores and field drains. The processes contributing to scale dependence in the data, and the possibility that certain processes dominate at particular scales, are discussed. The data presented here help to improve our spatial understanding of runoff and sediment transport in small agricultural catchments, and provide examples of the type of spatial dataset and the type of analysis that are essential if we are to develop models which are able to predict runoff and soil erosion accurately, and allow us to manage runoff and sediment transport effectively across scales. Copyright © 2011 John Wiley & Sons, Ltd.     

On the estimation of antecedent wetness conditions in rainfall–runoff modelling

The analysis of the physical processes involved in a conceptual model of soil water content balance is addressed with the objective of its application as a component of rainfall–runoff modelling. The model uses routinely measured meteorological variables (rainfall and air temperature) and incorporates a limited number of significant parameters. Its performance in estimating the soil moisture temporal pattern was tested through local measurements of volumetric water content carried out continuously on an experimental plot located in central Italy. The analysis was carried out for different periods in order to test both the representation of infiltration at the short time‐scale and drainage and evapotranspiration processes at the long time‐scale. A robust conceptual model was identified that incorporated the Green–Ampt approach for infiltration and a gravity‐driven approximation for drainage. A sensitivity analysis was performed for the selected model to assess the model robustness and to identify the more significant parameters involved in the principal processes that control the soil moisture temporal pattern. The usefulness of the selected model was tested for the estimation of the initial wetness conditions for rainfall–runoff modelling at the catchment scale. Specifically, the runoff characteristics (runoff depth and peak discharge) were found to be dependent on the pre‐event surface soil moisture. Both observed values and those estimated by the model gave good results. On the contrary, with the antecedent wetness conditions furnished by two versions of the antecedent precipitation index (API), large errors were obtained. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of road‐generated storm runoff on a small catchment response

The impact of road‐generated runoff on the hydrological response of a zero‐order basin was monitored for a sequence of 24 storm events. The study was conducted in a zero‐order basin (C1; 0·5ha) with an unpaved mountain road; an adjacent unroaded zero‐order basin (C2; 0·2 ha) with similar topography and lithology was used to evaluate the hydrological behaviour of the affected zero‐order basin prior to construction of the road. The impact of the road at the zero‐order basin scale was highly dependent on the antecedent soil‐moisture conditions, total storm precipitation, and to some extent rainfall intensity. At the beginning of the monitoring period, during dry antecedent conditions, road runoff contributed 50% of the total runoff and 70% of the peak flow from the affected catchment (C1). The response from the unroaded catchment was almost insignificant during dry antecedent conditions. As soil moisture increased, the road exerted less influence on the total runoff from the roaded catchment. For very wet conditions, the influence of road‐generated runoff on total outflow from the roaded catchment diminished to only 5·4%. Both catchments, roaded and unroaded, produced equivalent amount of outflow during very wet antecedent conditions on a unit area basis. The lag time between the rainfall and runoff peaks observed in the unroaded catchment during the monitoring period ranged from 0 to 4 h depending on the amount of precipitation and antecedent conditions, owing mainly to much slower subsurface flow pathways in the unroaded zero‐order basin. In contrast, the lag time in the roaded zero‐order basin was virtually nil during all storms. Copyright © 2009 John Wiley & Sons, Ltd.     

Storm‐event rainfall–runoff modelling approach for ungauged sites in Taiwan

This work develops a top‐down modelling approach for storm‐event rainfall–runoff model calibration at unmeasured sites in Taiwan. Twenty‐six storm events occurring in seven sub‐catchments in the Kao‐Ping River provided the analytical data set. Regional formulas for three important features of a streamflow hydrograph, i.e. time to peak, peak flow, and total runoff volume, were developed via the characteristics of storm event and catchment using multivariate regression analysis. Validation of the regional formulas demonstrates that they reasonably predict the three features of a streamflow hydrograph at ungauged sites. All of the sub‐catchments in the study area were then adopted as ungauged areas, and the three streamflow hydrograph features were calculated by the regional formulas and substituted into the fuzzy multi‐objective function for rainfall–runoff model calibration. Calibration results show that the proposed approach can effectively simulate the streamflow hydrographs at the ungauged sites. The simulated hydrographs more closely resemble observed hydrographs than hydrographs synthesized using the Soil Conservation Service (SCS) dimensionless unit hydrograph method, a conventional method for hydrograph estimation at ungauged sites in Taiwan. Copyright © 2008 John Wiley & Sons, Ltd.     

Rainfall–runoff modelling of ephemeral streams in the Valencia region (eastern Spain)

This paper presents preliminary results from the application of a transfer‐function rainfall–runoff model to ephemeral streams in Mediterranean Spain. Flow simulations have been conducted for two small catchments (Carraixet and Poyo basins), located in close proximity to one another yet with significantly different geological characteristics. Analysis of flow simulations for a number of high‐flow events has revealed the dominant influence of the rainfall on the catchment response, particularly for high‐rainfall events. Particular success has been attained modelling the highest magnitude events in both catchments and for all events in the faster responding (Poyo) catchment. In order to investigate the viability of the model for forecasting floods in ungauged catchments, additional investigations have been conducted by calibrating the model for one catchment (donor catchment) and then applying it to another (receptor catchment). The results indicate that this can be successful when either the donor catchment is a fast response catchment or when the model is calibrated using a high‐magnitude event in the donor catchment, providing that the modelled receptor catchment event is of a lower magnitude. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of catchment properties on runoff coefficient in a karst area of southwest China

Frequent human activities and climate change in the karst region of southwest China since the 1950s have led to the investigation of response of runoff to climate and catchment properties. Runoff coefficient (Rc) as an expression variable of the catchment response to rainfall is important to describe runoff dynamics and to estimate available streamflow for utilization. In this study, the equations of Rc associated with its attributors of climate condition and catchment property were derived using the Budyko framework. The equations were used to estimate relationship between the Rc and the attributors in the karst catchments in Guizhou province of southwest China. Analysis in the selected 23 karst catchments demonstrates that the spatial distribution of Rc is dominated by the catchment properties, such as the catchment properties of geology, slope and land use and land cover, rather than climate condition of drought index. Correlation analysis indicates that the catchment with a large slope usually has a high value of Rc, and a large proportion of carbonate rock in a catchment reduces Rc in the study area. Temporal increasing trend of Rc during 1961–2000 was found for most catchments in the study area. This increasing trend was primarily resulted from changes of catchment properties, e.g. deforestation in large areas of Guizhou province during the 1950s–1980s. Copyright © 2013 John Wiley & Sons, Ltd.     

Distribution of surface imperviousness in small urban catchments predicts runoff peak flows and stream flashiness

Urban growth is a global phenomenon, and the associated impacts on hydrology from land development are expected to increase, especially in peri‐urban catchments. It is well understood that greater peak flows and higher stream flashiness are associated with increased surface imperviousness and storm location. However, the effect of the distribution of impervious areas on runoff peak flow response and stream flashiness of peri‐urban catchments has not been well studied. In this study, a new geometric index, Relative Nearness of Imperviousness to the Catchment Outlet (RNICO), is defined to correlate imperviousness distribution of peri‐urban catchments with runoff peak flows and stream flashiness. Study sites include 21 suburban catchments in New York representing a range of drainage area from 5 to 189 km² and average imperviousness from 10% to 48%. On the basis of RNICO, all development patterns are divided into 3 classes: upstream, centralized, and downstream. Results showed an obvious increase in runoff peak flows and decrease in time to peak when moving from upstream to centralized and downstream urbanization classes. This indicates that RNICO is an effective tool for classifying urban development patterns and for macroscale understanding of the hydrologic behavior of small peri‐urban catchments, despite the complexity of urban drainage systems. We also found that the impact of impervious distribution on runoff peak flows and stream flashiness decreases with catchment scale. For small catchments (A < 40 km²), RNICO was strongly correlated with the average (R² = .95) and maximum (R² = .91) gaged peak flows due to the relatively efficient subsurface routing through stormwater and sewer networks. Furthermore, the Richards–Baker stream flashiness index in small catchments was positively correlated with fractional impervious area (R² = .84) and RNICO (R² = .87). For large catchments (A > 40 km²), the impact of impervious surface distribution on peak flows and stream flashiness was negligible due to the complex drainage network and great variability in travel times. This study emphasizes the need for greater monitoring of discharge in small peri‐urban catchments to support flood prediction at the local scale.     

Conceptual rainfall–runoff model with a two‐parameter,infinite characteristic time transfer function

A two‐parameter transfer function with an infinite characteristic time is proposed for conceptual rainfall–runoff models. The large time behaviour of the unit response is an inverse power function of time. The infinite characteristic time allows long‐term memory effects to be accounted for. Such effects are observed in mountainous and karst catchments. The governing equation of the model is a fractional differential equation in the limit of long times. Although linear, the proposed transfer function yields discharge signals that can usually be obtained only using non‐linear models. The model is applied successfully to two catchments, the Dud Koshi mountainous catchment in the Himalayas and the Durzon karst catchment in France. It compares favourably to the linear, non‐linear single reservoir models and to the GR4J model. With a single reservoir and a single transfer function, the model is capable of reproducing hysteretic behaviours identified as typical of long‐term memory effects. Computational efficiency is enhanced by approximating the infinite characteristic time transfer function with a sum of simpler, exponential transfer functions. This amounts to partitioning the reservoir into several linear sub‐reservoirs, the output discharges of which are easy to compute. An efficient partitioning strategy is presented to facilitate the practical implementation of the model. Copyright © 2015 John Wiley & Sons, Ltd.     

Interaction between runoff – bedrock groundwater in a steep headwater catchment underlain by sedimentary bedrock fractured by gravitational deformation

Recent studies have suggested the importance of the bedrock groundwater (BG) contribution in storm runoff in headwater catchments. However, few such studies have been conducted, and the study of different types of bedrock conditions is still ongoing. The role of BG in storm runoff is still poorly understood, particularly in headwater catchments underlain by relatively deep fractured bedrock. This study aims to clarify this role using hydrometric and hydrochemical observations of BG via boreholes and catchment discharge. The responses of the BG to rainfall are demonstrated to be fast and independent of the sediment cover. The BG exhibits different responses and flow paths that are largely controlled by the bedrock fracture system. The storm runoff in the studied catchment is characterized by rapid discharge response generally followed by a delayed discharge response. The peak of the delayed discharge is much faster than that observed in previous studies, and it is well correlated with the BG levels. A hydrograph separation was performed for two storm events using three end members: rainfall, shallow BG and deep BG. The results demonstrate that the delayed discharge is primarily composed of deep BG. Moreover, a significant contribution of shallow BG is observed during large precipitation events. Although we observed no physical evidence of direct contributions of BG in the catchment, the calculations presented in this study demonstrate that the BG controls the hydrological and hydrogeological response of the catchment to rainfall events. Copyright © 2015 John Wiley & Sons, Ltd.     

A soft hydrological monitoring approach for comparing runoff on a network of small poorly gauged catchments

Catchments in many parts of the world are either ungauged or poorly gauged, and the dominant processes governing their streamflow response are still poorly understood. The analysis of runoff coefficients provides essential insight into catchment response, particularly if both range of catchments and a range of events are compared. This paper investigates how well the hydrological runoff of 11 small, poorly gauged catchments with ephemeral streams (0·1‐0·6 km²) can be compared using estimated runoff with the associated uncertainty. Data of rainfall and water depth at a catchment's outlet were recorded using automatic logging equipment during 2008‐2009. The hydrological regime is intermittent and the annual precipitation ranged between 569 and 727 mm. Discharge was estimated using Manning's equation and channel cross‐section measurements. Innovative work has been performed under controlled experimental conditions to estimate Manning's coefficient values for the different cover types observed in studied streams: non‐aquatic vegetations (giant reed, bramble and thistle), grass and coarse granular deposits. The results show that estimates derived using roughness coefficients differ from those previously established for larger streams with aquatic vegetation. Catchment runoff was compared at both the event and the annual scale. The results indicate significant variability between the catchment's responses. This variability allows for classification in spite of all the uncertainty associated with runoff estimation. This study highlights the potential of using a network of poorly gauged catch ments. From almost no catchment understanding the proposed methodology allows to compare poorly gauged catchments and highlights similarity/dissimilarity between catchment responses. Copyright © 2011 John Wiley & Sons, Ltd.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.