Influence of rainfall and soil properties spatial aggregation on extreme flash flood response modelling: An evaluation based on the Sesia river basin,North Western Italy

High resolution radar rainfall fields and a distributed hydrologic model are used to evaluate the sensitivity of flood and flash flood simulations to spatial aggregation of rainfall and soil properties at catchment scales ranging from 75 to 983 km². Hydrologic modeling is based on a Hortonian infiltration model and a network-based representation of hillslope and channel flow. The investigation focuses on three extreme flood and flash flood events occurred on the Sesia river basin, North Western Italy, which are analysed by using four aggregation lengths ranging from 1 to 16 km. The influence of rainfall spatial aggregation is examined by using the flow distance as a spatial coordinate, hence emphasising the role of river network in the averaging of space–time rainfall. The effects of reduced and distorted rainfall spatial variability on peak discharge have been found particularly severe for the flash flood events, with peak errors up to 35% for rainfall aggregation of 16 km and at 983 km² catchment size. Effects are particularly remarkable when significant structured rainfall variability combines with relatively important infiltration volumes due to dry initial conditions, as this emphasises the non-linear character of the rainfall–runoff relationship. In general, these results confirm that the correct estimate of rainfall volume is not enough for the accurate reproduction of flash flood events characterised by large and structured rainfall spatial variability, even at catchment scales around 250 km². However, accurate rainfall volume estimation may suffice for less spatially variable flood events. Increasing the soil properties aggregation length exerts similar effects on peak discharge errors as increasing the rainfall aggregation length, for the cases considered here and after rescaling to preserve the rainfall volume. Moreover, peak discharge errors are roughly proportional to runoff volume errors, which indicates that the shape of the flood wave is influenced in a limited way by modifying the detail of the soil property spatial representation. Conversely, rainfall aggregation may exert a pronounced influence on the discharge peak by reshaping the spatial organisation of the runoff volumes and without a comparable impact on the runoff volumes.     

Coupling the Xinanjiang model to a kinematic flow model based on digital drainage networks for flood forecasting

The Xinanjiang model, which is a conceptual rainfall‐runoff model and has been successfully and widely applied in humid and semi‐humid regions in China, is coupled by the physically based kinematic wave method based on a digital drainage network. The kinematic wave Xinanjiang model (KWXAJ) uses topography and land use data to simulate runoff and overland flow routing. For the modelling, the catchment is subdivided into numerous hillslopes and consists of a raster grid of flow vectors that define the water flow directions. The Xinanjiang model simulates the runoff yield in each grid cell, and the kinematic wave approach is then applied to a ranked raster network. The grid‐based rainfall‐runoff model was applied to simulate basin‐scale water discharge from an 805‐km² catchment of the Huaihe River, China. Rainfall and discharge records were available for the years 1984, 1985, 1987, 1998 and 1999. Eight flood events were used to calibrate the model's parameters and three other flood events were used to validate the grid‐based rainfall‐runoff model. A Manning's roughness via a linear flood depth relationship was suggested in this paper for improving flood forecasting. The calibration and validation results show that this model works well. A sensitivity analysis was further performed to evaluate the variation of topography (hillslopes) and land use parameters on catchment discharge. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydrological processes in macrocatchment water harvesting in the arid region of Tunisia: the traditional system of tabias/Processus hydrologiques au sein d’un aménagement de collecte des eaux dans la région aride tunisienne: le système traditionnel des tabias

Abstract

Abstract In arid Tunisia, a tabia system is a traditional macrocatchment water harvesting system. It consists of a runoff area, which occupies two thirds of the slope and is traditionally used for grazing; and one to five cropped plots within U-shaped soil banks arranged in a cascade in the third downstream area. These ?run-on? areas accumulate and store the occasional runoff. Each soil bank is constructed with a discharge weir that allows modification of the flooded area and discharge of excess water towards downstream plots. Such a harvesting system, located in an area with 140 mm annual rainfall, was instrumented during four hydrological years (1995–1999) and 45 rainfall events were recorded. Eleven of these events gave a measurable inflow to at least one of the four plots. The observations showed that the traditional tabia system reduced total surface runoff from the catchment to essentially zero. The harvesting system significantly reduced peaks of surface runoff within the catchment, which also reduced erosion hazards. The cultivated area of about 5% of the total catchment could be supplied by a harvested water amount corresponding to about seven times the amount of each rainfall event larger than 20 mm.     

Radar rainfall estimation for flash flood forecasting in small urban watersheds

Radar rainfall estimation for flash flood forecasting in small, urban catchments is examined through analyses of radar, rain gage and discharge observations from the 14.3 km² Dead Run drainage basin in Baltimore County, Maryland. The flash flood forecasting problem pushes the envelope of rainfall estimation to time and space scales that are commensurate with the scales at which the fundamental governing laws of land surface processes are derived. Analyses of radar rainfall estimates are based on volume scan WSR-88D reflectivity observations for 36 storms during the period 2003–2005. Gage-radar analyses show large spatial variability of storm total rainfall over the 14.3 km² basin for flash flood producing storms. The ability to capture the detailed spatial variation of rainfall for flash flood producing storms by WSR-88D rainfall estimates varies markedly from event to event. As spatial scale decreases from the 14.3 km² scale of the Dead Run watershed to 1 km² (and the characteristic time scale of flash flood producing rainfall decreases from 1 h to 15 min) the predictability of flash flood response from WSR-88D rainfall estimates decreases sharply. Storm to storm variability of multiplicative bias in storm total rainfall estimates is a dominant element of the error structure of radar rainfall estimates, and it varies systematically over the warm season and with flood magnitude. Analyses of the 7 July 2004 and 28 June 2005 storms illustrate microphysical and dynamical controls on radar estimation error for extreme flash flood producing storms.     

Flooding and geomorphic impacts in a mountain torrent: Raise Beck,central Lake District,England

Raise Beck is a mountain torrent located in the central Lake District fells, northern England (drainage area of 1·27 km²). The torrent shows evidence of several major flood events, the most recent of which was in January 1995. This event caused a major channel avulsion at the fan apex diverting the main flood flow to the south, blocking the A591 trunk road and causing local flooding. The meteorological conditions associated with this event are described using local rainfall records and climatic data. Records show 164 mm of rainfall in the 24 hours preceding the flood. The peak flood discharge is reconstructed using palaeohydrological and rainfall–runoff methods, which provide discharge values of 27–74 m³ s^?1, and 4–6 m³ s^?1, respectively. The flood transported boulders with b‐axes up to 1400 mm. These results raise some important general questions about flood estimation in steep mountain catchments. The geomorphological impact of the event is evaluated by comparing aerial photographs from before and after the flood, along with direct field observations. Over the historical timescale the impact and occurrence of flooding is investigated using lichenometry, long‐term rainfall data, and documentary records. Two major historical floods events are identified in the middle of the nineteenth century. The deposits of the recent and historical flood events dominate the sedimentological evidence of flooding at Raise Beck, therefore the catchment is sensitive to high magnitude, low frequency events. Following the 1995 flood much of the lower catchment was channelized using rip‐rap bank protection, re‐establishing flow north towards Thirlmere. The likely success of this management strategy in containing future floods is considered, based on an analysis of channel capacities. It is concluded that the channelization scheme is only a short‐term solution, which would fail to contain the discharge of an event equivalent to the January 1995 flood. Copyright © 2002 John Wiley & Sons, Ltd.     

Hydrological Sciences Journal Contents 2002

Abstract

Flood frequency estimation is crucial in both engineering practice and hydrological research. Regional analysis of flood peak discharges is used for more accurate estimates of flood quantiles in ungauged or poorly gauged catchments. This is based on the identification of homogeneous zones, where the probability distribution of annual maximum peak flows is invariant, except for a scale factor represented by an index flood. The numerous applications of this method have highlighted obtaining accurate estimates of index flood as a critical step, especially in ungauged or poorly gauged sections, where direct estimation by sample mean of annual flood series (AFS) is not possible, or inaccurate. Therein indirect methods have to be used. Most indirect methods are based upon empirical relationships that link index flood to hydrological, climatological and morphological catchment characteristics, developed by means of multi-regression analysis, or simplified lumped representation of rainfall–runoff processes. The limits of these approaches are increasingly evident as the size and spatial variability of the catchment increases. In these cases, the use of a spatially-distributed, physically-based hydrological model, and time continuous simulation of discharge can improve estimation of the index flood. This work presents an application of the FEST-WB model for the reconstruction of 29 years of hourly streamflows for an Alpine snow-fed catchment in northern Italy, to be used for index flood estimation. To extend the length of the simulated discharge time series, meteorological forcings given by daily precipitation and temperature at ground automatic weather stations are disaggregated hourly, and then fed to FEST-WB. The accuracy of the method in estimating index flood depending upon length of the simulated series is discussed, and suggestions for use of the methodology provided.

Editor D. Koutsoyiannis     

Development of a thresholds approach for real‐time flash flood prediction in complex geomorphological river basins

Flash floods represent one of the deadliest and costliest natural disasters worldwide. The hydrological analysis of a flash flood event contributes in the understanding of the runoff creation process. This study presents the analysis of some flash flood events that took place in a complex geomorphological Mediterranean River basin. The objective of the present work is to develop the thresholds for a real‐time flash flood forecasting model in a complex geomorphological watershed, based on high‐frequency data from strategically located hydrological and meteorological telemetric stations. These stations provide hourly real‐time data which were used to determine hydrological and meteorological parameters. The main characteristics of various hydrographs specified in this study were the runoff coefficients, lag time, time to peak, and the maximum potential retention. The estimation of these hydrometeorological parameters provides the necessary information in order to successfully manage flash floods events. Especially, the time to peak is the most significant hydrological parameter that affects the response time of an oncoming flash flood event. A study of the above parameters is essential for the specification of thresholds which are related to the geomorphological characteristics of the river basin, the rainfall accumulation of an event, the rainfall intensity, the threshold runoff through karstic area, the season during which the rainfall takes place and the time intervals between the rainstorms that affect the soil moisture conditions. All these factors are combined into a real‐time‐threshold flash flood prediction model. Historical flash flood events at the downstream are also used for the validation of the model. An application of the proposed model is presented for the Koiliaris River basin in Crete, Greece. Copyright © 2011 John Wiley & Sons, Ltd.     

A statistically based runoff‐yield model coupling infiltration excess and saturation excess mechanisms

A statistically based runoff‐yield model is proposed in this paper. The model considers spatial heterogeneities of rainfall, soil infiltration capacity and soil water storage capacity that are main factors controlling runoff‐yield process. It assumes that the spatial variation of rainfall intensity at each time step can be characterized by a probability density function, which is estimated by matching the hyetograph through goodness‐of‐fit measure, whereas the spatial heterogeneities of soil infiltration capacity and soil water storage capacity are described by parabola‐type functions. Surface runoff is calculated according to infiltration excess mechanism; the statistical distribution of surface runoff rate can be deduced with the joint distribution of rainfall intensity and soil infiltration rate, thus obtaining a quasi‐analytical solution for surface runoff. Based on saturation excess mechanism, the groundwater flow (flows below the ground are collectively referred to as groundwater flow) is calculated by infiltration and the probability distribution of soil water storage capacity. Consequently, the total runoff is composed of infiltration excess and saturation excess runoff components. As an example, this model is applied to flood event simulation in Dongwan catchment, a semi‐humid region and a tributary of Yellow River in China. It indicates that the proposed runoff‐yield model could achieve acceptable accuracy. Copyright © 2011 John Wiley & Sons, Ltd.     

Assessment of overland flow variation and blue water production in a farmed semi-arid water harvesting catchment

Upgrading agriculture in semi-arid areas and ensuring its sustainability require an optimal management of rainfall partition between blue and green waters in the farmed water harvesting catchment. The main objective of this study is to analyze the influence of heterogeneous land use on the spatial and temporal variation of rainfall partitioning and blue water production within a typical farmed catchment located in north-eastern Tunisia. The catchment has an area of 2.6 km² and comprises at its outlet a dam, which retains the runoff water in a reservoir. Overland flow and soil water balance components were monitored during two cropping seasons (2000/2001 and 2001/2002) on a network of eleven plots of 2 m² each with different land use and soil characteristics. The hydrological balances of both the catchment and reservoir have been monitored since 1994.Observed data showed a very large temporal and spatial variability of overland flow within the catchment reflecting the great importance of total rainfall as well as land use. During the 2001/2002 season the results showed a large variation of the number of observed runoff events, from 27 to 39, and of the annual overland flow depths, from 8 mm (under vineyard on calcaric cambisols) up to 43 mm (under shrubs-pasture on haplic regosols), between the plots. The annual runoff amounts were moderate; they always corresponded to less than 15% of the annual rainfall amount whatever the observation scale. It was also observed that changes in land use in years with similar rainfall could lead to significant differences in blue water flow. An attempt for predicting the overland flow by the general linear regression approach showed an r² of 31%, the predictors used are the class of soil infiltration capacity, the initial moisture saturation ratio of the soil surface layer and the total rainfall amounts.These experimental results indicate that the variation in land use in a semi-arid catchment is a main factor of variation in soil surface conditions and explain the major role played by the former on hydrological behavior of the upstream area and on rainfall partition between overland flow and infiltration. Therefore, to predict the water harvesting capacities in terms of blue water production of a farmed catchment in semi-arid areas it seems essential to consider precisely its land use and its temporal evolution related to management practices.     

考虑降雨空间变化的随机产汇流模型

引入两个负指数型差值函数,估计降雨量的概率分布,以此描述流域降雨空间变异性问题.将降雨量空间统计分布与垂向混合产流模型耦合进行产流量计算,即对地表径流,采用超渗产流模式,根据降雨与土壤下渗能力的联合分布推求其空间分布;对地面以下径流,采用蓄满产流模式,以地表渗入量的均值作为输入,进行简化处理以提高其实用性;最终推导出总产流量概率分布函数计算公式.将流域概化成一个线性水库,并根据随机微分方程理论,推导任一计算时段洪水流量的概率分布,从而构建了一个完整的随机产汇流模型.以淮河支流黄泥庄流域为例进行应用研究,结果表明,该模型可提供洪水过程的概率预报,可用于防洪风险分析,若以概率分布的期望值作为确定性预报,亦具有较高精度.     

Artificial neural networks for daily rainfall—runoff modelling

Abstract

The application of artificial neural network (ANN) methodology for modelling daily flows during monsoon flood events for a large size catchment of the Narmada River in Madhya Pradesh (India) is presented. The spatial variation of rainfall is accounted for by subdividing the catchment and treating the average rainfall of each subcatchment as a parallel and separate lumped input to the model. A linear multiple-input single-output (MISO) model coupled with the ANN is shown to provide a better representation of the rainfall-runoff relationship in such large size catchments compared with linear and nonlinear MISO models. The present model provides a systematic approach for runoff estimation and represents improvement in prediction accuracy over the other models studied herein.     

Changes in flood regime by use of the modified curve number method

Abstract

This paper describes the use of a simple two stage rainfall-runoff model in which a curve number (CN) principle is used to calculate the soil water content and, subsequently, the rainfall contribution to direct runoff and groundwater flow. The maximum soil water retention, S, is used to express various characteristics of a catchment (infiltration rate, soil cover and land use, as in the CN method) relevant to flood formation. Using historical flood events, the model is calibrated, and the statistical distribution parameters of peak flows determined. With the same historical input data scenarios (rainfall), sets of flood hydrographs are simulated for various values of the parameter S, and corresponding distribution parameters of peak flows are determined. This procedure is used to demonstrate possible changes in flood regime to be expected due to changes of the catchment soil properties and its vegetation cover. A case study is presented for the River Hron catchment, area 582 km², in the mountainous region of central Slovakia.     

Estimating flash flood discharge in an ungauged mountain catchment with 2D hydraulic models and dendrogeomorphic palaeostage indicators

There is still wide uncertainty about past flash‐flood processes in mountain regions owing to the lack of systematic databases on former events. This paper presents a methodology to reconstruct peak discharge of flash floods and illustrates a case in an ungauged catchment in the Spanish Central System. The use of dendrogeomorphic evidence (i.e. scars on trees) together with the combined use of a two‐dimensional (2D) numerical hydraulic model and a terrestrial laser scan (TLS) has allowed estimation of peak discharge of a recent flash flood. The size and height distribution of scars observed in the field have been used to define three hypothetical scenarios (S_min or minimum scenario; S_med or medium scenario; and S_max or maximum scenario), thus illustrating the uncertainty involved in peak‐discharge estimation of flash floods in ungauged torrents. All scars analysed with dendrogeomorphic techniques stem from a large flash flood which took place on 17 December 1997. On the basis of the scenarios, peak discharge is estimated to 79 ± 14 m³ s^?1. The average deviation obtained between flood stage and expected scar height was ? 0·09 ± 0·53 m. From the data, it becomes obvious that the geomorphic position of trees is the main factor controlling deviation rate. In this sense, scars with minimum deviation were located on trees growing in exposed locations, especially on unruffled bedrock where the model predicts higher specific kinetic energy. The approach used in this study demonstrates the potential of tree‐ring analysis in palaeohydrology and for flood‐risk assessment in catchments with vulnerable goods and infrastructure. Copyright © 2010 John Wiley & Sons, Ltd.     

Summer flows in experimental catchments with different forest covers, Chile

Runoff and peak flows in four experimental catchments with different land uses are analyzed for summer periods. The catchments have a rainy temperate climate with annual precipitations between 2000 and 2500 mm, 70% of which is concentrated in the winter period between May and August. The final harvest of the forest plantation in one of these catchments generated increases in summer runoff. Also, differences between the maximum instantaneous discharge and the flow at the beginning of the storm then almost duplicated those registered in rainfall events of similar magnitude when the catchment was fully forested. Runoff analysis in this catchment is difficult because the two post-harvesting summer periods are much wetter than the two pre-harvesting ones but a double mass analysis shows the effect of harvesting clearly. In a paired catchment study, low cover in one of the two neighbour catchments explains higher direct runoff and base flows although lower maximum instantaneous specific discharge occurred in the less vegetated but larger catchment. Low vegetation cover explains increases in summer flows, although the size, topography, rainfall conditions, road density, extent of affected area and runoff generation processes play an important role in the hydrological effects of different land uses.     

Rainfall intensity affects runoff responses in a semi-arid catchment

The source and hydrochemical makeup of a stream reflects the connectivity between rainfall, groundwater, the stream, and is reflected to water quantity and quality of the catchment. However, in a semi-arid, thick, loess covered catchment, temporal variation of stream source and event associated behaviours are lesser known. Thus, the isotopic and chemical hydrographs in a widely distributed, deep loess, semi-arid catchment of the northern Chinese Loess Plateau were characterized to determine the source and hydrochemical behaviours of the stream during intra-rainfall events. Rainfall and streamflow were sampled during six hydrologic events coupled with measurements of stream baseflow and groundwater. The deuterium isotope (²H), major ions (Cl⁻, SO₄²⁻, NO₃⁻, Ca²⁺, K⁺, Mg²⁺, and Na⁺) were evaluated in water samples obtained during rainfall events. Temporal variation of ²H and Cl⁻ measured in the groundwater and stream baseflow prior to rainfall was similar; however, the isotope compositions of the streamflow fluctuated significantly and responded quickly to rainfall events, likely due to an infiltration excess, overland dominated surface runoff during torrential rainfall events. Time source separation using ²H demonstrated greater than 72% on average, the stream composition was event water during torrential rainfall events, with the proportion increasing with rainfall intensity. Solutes concentrations in the stream had loglinear relationships with stream discharge, with an outling anomaly with an example of an intra-rainfall event on Oct. 24, 2015. Stream Cl⁻ behaved nonconservative during rainfall events, temporal variation of Cl⁻ indicated a flush and washout at the onset of small rainfall events, a dilution but still high concentration pattern in high discharge and old water dominated in regression flow period. This study indicates rainfall intensity affects runoff responses in a semi-arid catchment, and the stored water in the thick, loess covered areas was less connected with stream runoff. Solute transport may threaten water quality in the area, requiring further analysis of the performance of the eco-restoration project.     

Hydrological impact and potential flooding of convective rain cells in a semi-arid environment

Abstract

t Spatio-temporal storm properties have a large impact on catchment hydrological response. The sensitivity of simulated flash floods to convective rain-cell characteristics is examined for an extreme storm event over a 94 km² semi-arid catchment in southern Israel. High space–time resolution weather radar data were used to derive and model convective rain cells that then served as input into a hydrological model. Based on alterations of location, direction and speed of a major rain cell, identified as the flooding cell for this case, the impacts on catchment rainfall and generated flood were examined. Global sensitivity analysis was applied to identify the most important factors affecting the flash flood peak discharge at the catchment outlet. We found that the flood peak discharge could be increased three-fold by relatively small changes in rain-cell characteristics. We assessed that the maximum flash flood magnitude that this single rain cell can produce is 175 m³/s, and, taking into account the rest of the rain cells, the flash flood peak discharge can reach 260 m³/s.

Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Morin, E. and Yakir, H., 2013. Hydrological impact and potential flooding of convective rain cells in a semi-arid environment. Hydrological Sciences Journal, 59 (7), 1275–1284. http://dx.doi.org/10.1080/02626667.2013.841315     

Process controls on the statistical flood moments ‐ a data based analysis

In this paper, the controls of different indicators on the statistical moments (i.e. mean annual flood (MAF), coefficient of variation (CV) and skewness (CS)) of the maximum annual flood records of 459 Austrian catchments are analysed. The process controls are analysed in terms of the correlation of the flood moments within five hydrologically homogeneous regions to two different types of indicators. Indicators of the first type are static catchment attributes, which are associated with long‐term observations such as mean annual precipitation, the base flow index, and the percentage of catchment area covered by a geological unit or soil type. Indicators of the second type are dynamic catchment attributes that are associated with the event scale. Indicators of this type used in the study are event runoff coefficients and antecedent rainfall. The results indicate that MAF and CV are strongly correlated with indicators characterising the hydro‐climatic conditions of the catchments, such as mean annual precipitation, long‐term evaporation and the base flow index. For the catchments analysed, the flood moments are not significantly correlated with static catchment attributes representing runoff generation, such as geology, soil types, land use and the SCS curve number. Indicators of runoff generation that do have significant predictive power for flood moments are dynamic catchment attributes such as the mean event runoff coefficients and mean antecedent rainfall. The correlation analysis indicates that flood runoff is, on average, more strongly controlled by the catchment moisture state than by event rainfall. Copyright © 2008 John Wiley & Sons, Ltd.     

Estimating Flood Quantiles on the Basis of Multi-Event Rainfall Simulation — Case Study

This paper presents an approach to estimating the probability distribution of annual discharges Q based on rainfall-runoff modelling using multiple rainfall events. The approach is based on the prior knowledge about the probability distribution of annual maximum daily totals of rainfall P in a natural catchment, random disaggregation of the totals into hourly values, and rainfall-runoff modelling. The presented Multi-Event Simulation of Extreme Flood method (MESEF) combines design event method based on single-rainfall event modelling, and continuous simulation method used for estimating the maximum discharges of a given exceedance probability using rainfall-runoff models. In the paper, the flood quantiles were estimated using the MESEF method, and then compared to the flood quantiles estimated using classical statistical method based on observed data.     

Sediment transport in Swiss torrents

Sediment loads have been measured in six Swiss mountain torrents over several decades. Most of these torrent catchments are situated in the prealpine belt. They have catchment areas of between 0·5 and 1·7 km². Bedslopes at the measuring sites vary between 5 and 17 per cent, and peak discharges up to 12 m³ s⁻¹ have been recorded. Geophone sensors installed in the Erlenbach stream allow bedload transport activity to be monitored and sediment volumes associated with each flood event to be determined. A detailed analysis of the measurements in this stream results in an empirical equation in which the sediment load per flood event is expressed as a function of the effective runoff volume (discharges above the threshold for bedload motion) and of the normalized peak discharge. For the total of 143 investigated flood events in the Erlenbach stream, the deviation of the predicted from the measured value is within a factor of two for more than two-thirds of all events. A distinction can be made between summer and winter events in analysing the bedload transport events. The summer events, mainly caused by thunderstorms, transport comparatively larger sediment loads than the winter events. For the other investigated streams, the periods of the deposited sediment volume surveys cover in general several flood events. An analysis is performed analogous to that for the Erlenbach stream. The sediment loads show a similar dependency on the two factors effective runoff volume and normalized peak discharge. However, the exponents of these factors in the power law expressions differ from stream to stream. A comparison of the investigated stream shows that some of the variation can be explained by considering the bedslope above the measuring site. The inclusion of a bedslope factor is in agreement with laboratory investigations on bedload transport. © 1997 John Wiley & Sons, Ltd.     

Rainfall–runoff response and event-based runoff coefficients in a humid area (northwest Spain)

Abstract

The hydrological response of a small agroforestry catchment in northwest Spain (Corbeira catchment, 16 km²) is analysed, with particular focus on rainfall events. Fifty-four rainfall–runoff events, from December 2004 to September 2007, were used to analyse the principal hydrological patterns and show which factors best explain the hydrological response. The nonlinearity between rainfall and runoff showed that the variability in the hydrological response of the catchment was linked to the seasonal dynamics of the rainfall and, to a lesser extent, to evapotranspiration. The runoff coefficient, estimated as the ratio between direct runoff and rainfall volume, on an event basis, was analysed as a function of rainfall characteristics (amount and intensity) and the initial catchment state conditions prior to an event, such as pre-event baseflow and antecedent rainfall index. The results revealed that the hydrological response depends both on the soil humidity conditions at the start of the event and on rainfall amount, whereas rainfall intensity presented only a significant correlation with discharge increment. The antecedent conditions seem to be a key point in runoff production, and they explain much of the response. The hydrographs are characterized by a steep rising limb, a relatively narrow peak discharge and slow recession limb. These data and the observations suggest that the subsurface flow is the dominant runoff process.

Editor Z.W. Kundzewicz; Associate editor T. Wagener

Citation Rodríguez-Blanco, M.L., Taboada-Castro, M.M. and Taboada-Castro, M.T., 2012. Rainfall–runoff response and event-based runoff coefficients in a humid area (northwest Spain). Hydrological Sciences Journal, 57 (3), 445–459.