Distributed rainfall‐runoff modelling for flood frequency estimation and flood forecasting

Nowadays, in the scientific literature many rainfall‐runoff (RR) models are available ranging from simpler ones, with a limited number of parameters, to highly complex ones, with many parameters. Therefore, the selection of the best structure and parameterisation for a model is not straightforward as it is dependent on a number of factors: climatic conditions, catchment characteristics, temporal and spatial resolution, model objectives, etc. In this study, the structure of a continuous semi‐distributed RR model, named MISDc (‘Modello Idrologico Semi‐Distribuito in continuo’) developed for flood simulation in the Upper Tiber River (central Italy) is presented. Most notably, the methodology employed to detect the more relevant processes involved in the modelling of high floods, and hence, to build the model structure and its parameters, is developed. For this purpose, an intense activity of monitoring soil moisture and runoff in experimental catchments was carried out allowing to derive a parsimonious and reliable continuous RR model operating at an hourly (or smaller) time scale. Specifically, in order to determine the catchment hydrological response, the important role of the antecedent wetness conditions is emphasized. The application of MISDc both for design flood estimation and for flood forecasting is reported here demonstrating its reliability and also its computational efficiency, another important factor in hydrological practice. As far as the flood forecasting applications are concerned, only the accuracy of the model in reproducing discharge hydrographs by assuming rainfall correctly known throughout the event is investigated indepth. In particular, the MISDc has been implemented in the framework of Civil Protection activities for the Upper Tiber River basin. Copyright © 2011 John Wiley & Sons, Ltd.     

Flood simulation using the gauge‐adjusted radar rainfall and physics‐based distributed hydrologic model

This paper reports the results of an investigation into flood simulation by areal rainfall estimated from the combination of gauged and radar rainfalls and a rainfall–runoff model on the Anseong‐cheon basin in the southern part of Korea. The spatial and temporal characteristics and behaviour of rainfall are analysed using various approaches combining radar and rain gauges: (1) using kriging of the rain gauge alone; (2) using radar data alone; (3) using mean field bias (MFB) of both radar and rain gauges; and (4) using conditional merging technique (CM) of both radar and rain gauges. To evaluate these methods, statistics and hyetograph for rain gauges and radar rainfalls were compared using hourly radar rainfall data from the Imjin‐river, Gangwha, rainfall radar site, Korea. Then, in order to evaluate the performance of flood estimates using different rainfall estimation methods, rainfall–runoff simulation was conducted using the physics‐based distributed hydrologic model, Vflo^?. The flood runoff hydrograph was used to compare the calculated hydrographs with the observed one. Results show that the rainfall field estimated by CM methods improved flood estimates, because it optimally combines rainfall fields representing actual spatial and temporal characteristics of rainfall. Copyright © 2008 John Wiley & Sons, Ltd.     

Simulation of flood hydrographs based on lumped and semi-distributed models for two tropical catchments in Kenya

Abstract

The design and operation of flood management systems require computation of flood hydrographs for both design floods and flood forecasting purposes, since observed data are usually inadequate for these tasks. This is particularly relevant for most developing countries, i.e. mainly for tropical catchments. One possible way of obtaining information about flood hydrographs is through the use of rainfall—runoff models. Two such models, namely the Bochum model and the Nash Cascade—Diskin Infiltration model, which are semi-distributed and lumped models, respectively, were used in the present study. These models were applied to two catchments in Kenya with drainage areas of 6.71 km² and 26.03 km². A set of 13 selected rainfall—runoff events was used to calibrate and validate the models. The physical parameters required by the models were derived from catchment characteristics using GIS and remote sensing data while the conceptual parameters were obtained by optimization. The flood hydrographs simulated using the parameters so derived indicated that it is possible to use the two models in this tropical environment.     

Sensitivity of hydrological models to uncertainty in rainfall input

Abstract

Different approaches used in hydrological modelling are compared in terms of the way each one takes the rainfall data into account. We examine the errors associated with accounting for rainfall variability, whether in hydrological modelling (distributed vs lumped models) or in computing catchment rainfall, as well as the impact of each approach on the representativeness of the parameters it uses. The database consists of 1859 rainfall events, distributed on 500 basins, located in the southeast of France with areas ranging from 6.2 to 2851 km². The study uses as reference the hydrographs computed by a distributed hydrological model from radar rainfall. This allows us to compare and to test the effects of various simplifications to the process when taking rainfall information (complete rain field vs sampled rainfall) and rainfall–runoff modelling (lumped vs distributed) into account. The results appear to show that, in general, the sampling effect can lead to errors in discharge at the outlet that are as great as, or even greater than, those one would get with a fully lumped approach. We found that small catchments are more sensitive to the uncertainties in catchment rainfall input generated by sampling rainfall data as seen through a raingauge network. Conversely, the larger catchments are more sensitive to uncertainties generated when the spatial variability of rainfall events is not taken into account. These uncertainties can be compensated for relatively easily by recalibrating the parameters of the hydrological model, although such recalibrations cause the parameter in question to completely lose physical meaning.

Citation Arnaud, P., Lavabre, J., Fouchier, C., Diss, S. & Javelle, P. (2011) Sensitivity of hydrological models to uncertainty of rainfall input. Hydrol. Sci. J. 56(3), 397–410.     

The effect of differences between rainfall measurement techniques on groundwater and discharge simulations in a lowland catchment

Several rainfall measurement techniques are available for hydrological applications, each with its own spatial and temporal resolution and errors. When using these rainfall datasets as input for hydrological models, their errors and uncertainties propagate through the hydrological system. The aim of this study is to investigate the effect of differences between rainfall measurement techniques on groundwater and discharge simulations in a lowland catchment, the 6.5‐km² Hupsel Brook experimental catchment. We used five distinct rainfall data sources: two automatic raingauges (one in the catchment and another one 30 km away), operational (real‐time and unadjusted) and gauge‐adjusted ground‐based C‐band weather radar datasets and finally a novel source of rainfall information for hydrological purposes, namely, microwave link data from a cellular telecommunication network. We used these data as input for the, a recently developed rainfall‐runoff model for lowland catchments, and intercompared the five simulated discharges time series and groundwater time series for a heavy rainfall event and a full year. Three types of rainfall errors were found to play an important role in the hydrological simulations, namely: (1) Biases, found in the unadjusted radar dataset, are amplified when propagated through the hydrological system; (2) Timing errors, found in the nearest automatic raingauge outside the catchment, are attenuated when propagated through the hydrological system; (3) Seasonally varying errors, found in the microwave link data, affect the dynamics of the simulated catchment water balance. We conclude that the hydrological potential of novel rainfall observation techniques should be assessed over a long period, preferably a full year or longer, rather than on an event basis, as is often done. Copyright © 2016 The Authors. Hydrological Processes. Published by John Wiley & Sons Ltd.     

Coupling between weather radar rainfall data and a distributed hydrological model for real-time flood forecasting / Couplage de données pluviométriques issues de radars météorologiques et d’un modèle hydrologique distribué pour la prévision de crue en temps réel

Abstract

Abstract After the destructive flood in 1998, the Chinese government planned to build national weather radar networks and to use radar data for real-time flood forecasting. Hence, coupling of weather radar rainfall data and a hydrological (Xinanjiang) model became an important issue. The present study reports on experience in such coupling at the Shiguanhe watershed. After having corrected the radar reflectivity and the attenuation data, the weather radar rainfall was estimated and then corrected in real time using a Kalman filter. In general, the precipitation estimated from weather radar is reasonably accurate in most of the catchment investigated, after corrections as above. Compared to the results simulated by raingauge data, the simulations based on the weather radar data are of similar accuracy. Present research results show that rainfall estimated from the weather radar, the radar data correction method, the method of coupling, and the Xinanjiang model lend themselves well to application in operational real-time flood forecasting.     

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.     

Hydrological appraisal of rainfall estimates from radar,satellite, raingauge and satellite–gauge combination on the Qinhuai River Basin,China

ABSTRACT

Multisource rainfall products can be used to overcome the absence of gauged precipitation data for hydrological applications. This study aims to evaluate rainfall estimates from the Chinese S-band weather radar (CINRAD-SA), operational raingauges, multiple satellites (CMORPH, ERA-Interim, GPM, TRMM-3B42RT) and the merged satellite–gauge rainfall products, CMORPH-GC, as inputs to a calibrated probability distribution model (PDM) on the Qinhuai River Basin in Nanjing, China. The Qinhuai is a middle-sized catchment with an area of 799 km². All sources used in this study are capable of recording rainfall at high spatial and temporal resolution (3 h). The discrepancies between satellite and radar data are analysed by statistical comparison with raingauge data. The streamflow simulation results from three flood events suggest that rainfall estimates using CMORPH-GC, TRMM-3B42RT and S-band radar are more accurate than those using the other rainfall sources. These findings indicate the potential to use satellite and radar data as alternatives to raingauge data in hydrological applications for ungauged or poorly gauged basins.     

Distributed hydrological modelling using weather radar in gauged and ungauged basins

Distributed hydrological modelling using space–time estimates of rainfall from weather radar provides a natural approach to area-wide flood forecasting and warning at any location, whether gauged or ungauged. However, radar estimates of rainfall may lack consistent, quantitative accuracy. Also, the formulation of hydrological models in distributed form may be problematic due to process complexity and scaling issues. Here, the aim is to first explore ways of improving radar rainfall accuracy through combination with raingauge network data via integrated multiquadric methods. When the resulting gridded rainfall estimates are employed as input to hydrological models, the simulated river flows show marked improvements when compared to using radar data alone. Secondly, simple forms of physical–conceptual distributed hydrological model are considered, capable of exploiting spatial datasets on topography and, where necessary, land-cover, soil and geology properties. The simplest Grid-to-Grid model uses only digital terrain data to delineate flow pathways and to control runoff production, the latter by invoking a probability-distributed relation linking terrain slope to soil absorption capacity. Model performance is assessed over nested river basins in northwest England, employing a lumped model as a reference. When the distributed model is used with the gridded radar-based rainfall estimators, it shows particular benefits for forecasting at ungauged locations.     

Using stochastic state‐space modelling to specify and reduce hydrograph error

A Bayesian post‐processor is used to generate a representation of the likely hydrograph forecast flow error distribution using raingauge and radar input to a stochastic catchment model and its deterministic equivalent. A hydrograph ensemble is so constructed. Experiments are analysed using the model applied to the River Croal in north‐west England. It is found that for rainfall input to the model having errors less than 3mm h^?1, corresponding to about a 15% error in peak flow, the stochastic model outperforms the deterministic model. The range of hydrographs associated with the different model simulations and the measured hydrographs are compared. The significant improvement possible using a stochastic approach is demonstrated for a specific case study, although the mean hydrograph derived using the stochastic model has an error range associated with it. Copyright © 2007 John Wiley & Sons, Ltd.     

Estimation of discharges of water flows and debris floods in a small watershed

Floods in small mountainous watersheds cover a wide spectrum of flow. They can range from clear water flows and hyperconcentrated flows to debris floods and debris flows, and calculation of the peak discharge is crucial for predicting and mitigating such hazards. To determine the optimal approach for discharge estimation, this study compared water flow monitoring hydrographs to investigate the performance of five hydrological models that incorporate different runoff yields and influx calculation methods. Two of the models performed well in simulating the peak discharge, peak time, and total flow volume of the water flood. The ratio (γ) of the monitored debris flood discharge (Q_d) to the simulated water flow discharge (Q_w) was investigated. Qualitatively, γ initially increased with Q_w but then decreased when Q_w exceeded a certain threshold, which corresponded to rainfall of 95 and 120 mm in a 6- and 24-h event with a normal distribution of precipitation, respectively. The decrease might be attributable to a threshold of sediment availability being reached, beyond which increased flow rate is not matched by increased sediment input in the large watershed. Uncertainty of hydrological calculation was evaluated by dividing the catchment into sub-basins and adopting different rainfall time steps as input. The efficiency of using a distributed simulation exhibited marginal improvement potential compared with a lumped simulation. Conversely, the rainfall time step input significantly affected the simulation results by delaying the peak time and decreasing the peak discharge. This research demonstrates the applicability of a discharge estimation method that combines a hydrological water flow simulation and an estimation of γ. The results were verified on the basis of monitored flow densities and videos obtained in two watersheds with areas of 2.34 and 32.4 km².     

Effect of growing watershed imperviousness on hydrograph parameters and peak discharge

An increasing impervious area is quickly extending over the Wu‐Tu watershed due to the endless demands of the people. Generally, impervious paving is a major result of urbanization and more recently has had the potential to produce more enormous flood disasters than those of the past. In this study, 40 available rainfall–runoff events were chosen to calibrate the applicable parameters of the models and to determine the relationships between the impervious surfaces and the calibrated parameters. Model inputs came from the outcomes of the block kriging method and the non‐linear programming method. In the optimal process, the shuffled complex evolution method and three criteria were applied to compare the observed and simulated hydrographs. The tendencies of the variations of the parameters with their corresponding imperviousness were established through regression analysis. Ten cases were used to examine the established equations of the parameters and impervious covers. Finally, the design flood routines of various return periods were furnished through use of approaches containing a design storm, block kriging, the SCS model, and a rainfall‐runoff model with established functional relationships. These simulated flood hydrographs were used to compare and understand the past, present, and future hydrological conditions of the watershed studied. In the research results, the time to peak of flood hydrographs for various storms was diminished approximately from 11 h to 6 h in different decrements, whereas peak flow increased respectively from 127 m³ s^?1 to 629 m³ s^?1 for different storm intensities. In addition, this study provides a design diagram for the peak flow ratio to help engineers and designers to construct hydraulic structures efficiently and prevent possible damage to human life and property. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact assessment of climatic and land-use changes on flood runoff in southeast Queensland

Abstract

Over the past century, land-use has changed in southeast Queensland, and when coupled with climatic change, the risk of flooding has increased. This research aims to examine impacts of climate and land-use changes on flood runoff in southeast Queensland, Australia. A rainfall–runoff model, RORB, was calibrated and validated using observed flood hydrographs for one rural and one urbanized catchment, for 1961–1990. The validated model was then used to generate flood hydrographs using projected rainfall based on two climate models: the Geophysical Fluid Dynamics Laboratory Climate Model 2.1 (GFDL CM2.1) and the Conformal-Cubic Atmospheric Model (CCAM), for 2016–2045. Projected daily rainfall for the two contrasting periods was used to derive adjustment factors for a given frequency of occurrence. Two land-use change scenarios were used to evaluate likely impacts. Based on the projected rainfall, the results showed that, in both catchments, future flood magnitudes are unlikely to increase for large flood events. Extreme land-use change would significantly impact flooding in the rural catchment, but not the urbanized catchment.

Editor Z.W. Kundzewicz; Associate editor Y. Gyasi-Agyei     

Integration of gauge and radar rainfall to enable best simulation of hydrological parameters

This study is about use of spatially distributed rain in physically based hydrological models. In recent years, spatially distributed radar rainfall data have become available. The distributed radar rain is used to precisely model hydrologic processes and it is more realistic than the past practice of distribution methods like Thiessen polygons. Radar provides a highly accurate spatial distribution of rainfall and greatly improves the basin average rainfall estimates. However, quantification of the exact amount of rainfall from radar observation is relatively difficult. Thus, the fundamental idea of this study is to apply hourly gauge and radar rainfall data in a distributed hydrological model to simulate hydrological parameters. Hence the comparison is made between the outcomes of the WetSpa model from radar rainfall distribution and gauge rainfall distributed by the Thiessen polygon technique. The comparative plots of the hydrograph and the results of hydrological components such as evapotranspiration, surface runoff, soil moisture, recharge and interflow, reflect the spatially distributed radar input performing well for model outflow simulation.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR F. Pappenberger     

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatial-temporal rainfall modelling for flood risk estimation

Some recent developments in the stochastic modelling of single site and spatial rainfall are summarised. Alternative single site models based on Poisson cluster processes are introduced, fitting methods are discussed, and performance is compared for representative UK hourly data. The representation of sub-hourly rainfall is discussed, and results from a temporal disaggregation scheme are presented. Extension of the Poisson process methods to spatial-temporal rainfall, using radar data, is reported. Current methods assume spatial and temporal stationarity; work in progress seeks to relax these restrictions. Unlike radar data, long sequences of daily raingauge data are commonly available, and the use of generalized linear models (GLMs) (which can represent both temporal and spatial non-stationarity) to represent the spatial structure of daily rainfall based on raingauge data is illustrated for a network in the North of England. For flood simulation, disaggregation of daily rainfall is required. A relatively simple methodology is described, in which a single site Poisson process model provides hourly sequences, conditioned on the observed or GLM-simulated daily data. As a first step, complete spatial dependence is assumed. Results from the River Lee catchment, near London, are promising. A relatively comprehensive set of methodologies is thus provided for hydrological application.     

Flood hydrograph attenuation induced by a reservoir system: analysis with a distributed rainfall‐runoff model

Spatially distributed hydrologic models can be effectively utilized for flood event simulation over basins where a complex system of reservoirs affecting the natural flow regime is present. Flood peak attenuation through mountain reservoirs can, in fact, mitigate the impact of major floods in flood‐prone areas of the lower river valley. Assessment of this effect for a complex reservoir system is performed with a spatially distributed hydrologic model where the surface runoff formation and the hydraulic routing through each reservoir and the river system are performed at a fine spatial and time resolution. The Toce River basin is presented as a case study, because of the presence of 14 active hydroelectric dams that affect the natural flow regime. A recent extreme flood event is simulated using a multi‐realization kriging method for modelling the spatial distribution of rainfall. A sensitivity analysis of the key elements of the distributed hydrologic model is also performed. The flood hydrograph attenuation is assessed. Several possible reservoir storage conditions are used to characterize the initial condition of each reservoir. The results demonstrate how a distributed hydrologic model can contribute to defining strategies for reservoir management in flood mitigation. Copyright © 2004 John Wiley & Sons, Ltd.     

An evaluation of a traditional and a neural net modelling approach to flood forecasting for an upland catchment

This study evaluates two (of the many) modelling approaches to flood forecasting for an upland catchment (the River South Tyne at Haydon Bridge, England). The first modelling approach utilizes ‘traditional’ hydrological models. It consists of a rainfall–runoff model (the probability distributed model, or PDM) for flow simulation in the upper catchment. Those flows are then routed to the lower catchment using two kinematic wave (KW) routing models. When run in forecast‐mode, the PDM and KW models utilize model updating procedures. The second modelling approach uses neural network models, which use a ‘pattern‐matching’ process to produce model forecasts.Following calibration, the models are evaluated in terms of their fit to continuous stage data and flood event magnitudes and timings within a validation period. Forecast times of 1 h, 2 h and 4 h are selected (the catchment has a response time of approximately 4 h). The ‘traditional’ models generally perform adequately at all three forecast times. The neural networks produce reasonable forecasts of small‐ to medium‐sized flood events but have difficulty in forecasting the magnitude of the larger flood events in the validation period. Possible modifications to the latter approach are discussed. © Crown copyright 2002. Reproduced with the permission of Her Majesty's stationery office. Published by John Wiley & Sons, Ltd.     

A probabilistic framework for floodplain mapping using hydrological modeling and unsteady hydraulic modeling

ABSTRACT

Prediction of design hydrographs is key in floodplain mapping using hydraulic models, which are either steady state or unsteady. The former, which require only an input peak, substantially overestimate the volume of water entering the floodplain compared to the more realistic dynamic case simulated by the unsteady models that require the full hydrograph. Past efforts to account for the uncertainty of boundary conditions using unsteady hydraulic modeling have been based largely on a joint flood frequency–shape analysis, with only a very limited number of studies using hydrological modeling to produce the design hydrographs. This study therefore presents a generic probabilistic framework that couples a hydrological model with an unsteady hydraulic model to estimate the uncertainty of flood characteristics. The framework is demonstrated on the Swannanoa River watershed in North Carolina, USA. Given its flexibility, the framework can be applied to study other sources of uncertainty in other hydrological models and watersheds.     

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.