A rainfall‐runoff model parameterized from GIS and runoff data

A guiding principle in hydrological modelling should be to keep the number of calibration parameters to a minimum. A reduced number of parameters to be calibrated, while maintaining the accuracy and detail required by modern hydrological models, will reduce parameter and model structure uncertainty and improve model diagnostics. In this study, the dynamics of runoff are derived from the distribution of distances from points in the catchments to the nearest stream. This distribution is unique for each catchment and can be determined from a geographical information system. The distribution of distances, will, when a celerity of (subsurface) flow is introduced, provide a distribution of travel times, or a unit hydrograph (UH). For spatially varying levels of saturation deficit, we have different celerities and, hence, different UHs. Runoff is derived from the superposition of the different UHs. This study shows how celerities can be estimated if we assume that recession events represent the combined UHs for different levels of saturation deficit. A new soil moisture routine which estimates saturated and unsaturated volumes of subsurface water and with only one parameter to calibrate is included in the new model. The performance of the new model is compared with that of the Swedish HBV model and is found to perform equally well for eight Norwegian catchments although the number of parameters to be calibrated in the module concerning soil moisture and runoff dynamics is reduced from seven in the HBV model to one in the new model. It is also shown that the new model has a more realistic representation of the subsurface hydrology. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessment of rainfall‐runoff models based upon wavelet analysis

A basic hypothesis is proposed: given that wavelet‐based analysis has been used to interpret runoff time‐series, it may be extended to evaluation of rainfall‐runoff model results. Conventional objective functions make certain assumptions about the data series to which they are applied (e.g. uncorrelated error, homoscedasticity). The difficulty that objective functions have in distinguishing between different realizations of the same model, or different models of the same system, is that they may have contributed in part to the occurrence of model equifinality. Of particular concern is the fact that the error present in a rainfall‐runoff model may be time dependent, requiring some form of time localization in both identification of error and derivation of global objective functions. We explore the use of a complex Gaussian (order 2) wavelet to describe: (1) a measured hydrograph; (2) the same hydrograph with different simulated errors introduced; and (3) model predictions of the same hydrograph based upon a modified form of TOPMODEL. The analysis of results was based upon: (a) differences in wavelet power (the wavelet power error) between the measured hydrograph and both the simulated error and modelled hydrographs; and (b) the wavelet phase. Power difference and wavelet phase were used to develop two objective functions, RMSE(power) and RMS(phase), which were shown to distinguish between simulated errors and model predictions with similar values of the commonly adopted Nash‐Sutcliffe efficiency index. These objective functions suffer because they do not retain time, frequency or time‐frequency localization. Consideration of wavelet power spectra and time‐ and frequency‐integrated power spectra shows that the impacts of different types of simulated error can be seen through retention of some localization, especially in relation to when and the scale over which error was manifest. Theoretical objections to the use of wavelet analysis for this type of application are noted, especially in relation to the dependence of findings upon the wavelet chosen. However, it is argued that the benefits of localization and the qualitatively low sensitivity of wavelet power and phase to wavelet choice are sufficient to warrant further exploration of wavelet‐based approaches to rainfall‐runoff model evaluation. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of the catchment runoff coefficient on the performance of TOPMODEL in rainfall–runoff modelling

Effects of the catchment runoff coefficient on the performance of TOPMODEL in simulating catchment rainfall–runoff relationships are investigated in this paper, with an aim to improve TOPMODEL's simulation efficiency in catchments with a low runoff coefficient. Application of TOPMODEL in the semi‐arid Yihe catchment, with an area of 2623 km² in the Yellow River basin of China, produced a Nash–Sutcliffe model efficiency of about 80%. To investigate how the catchment runoff coefficient affects the performance of TOPMODEL, the whole observed discharge series of the Yihe catchment is multiplied with a larger‐than‐unity scale factor to obtain an amplified discharge series. Then TOPMODEL is used to simulate the amplified discharge series given the original rainfall and evaporation data. For a set of different scale factors, TOPMODEL efficiency is plotted against the corresponding catchment runoff coefficient and it is found that the efficiency of TOPMODEL increases with the increasing catchment runoff coefficient before reaching a peak (e.g. about 90%); after the peak, however, the efficiency of TOPMODEL decreases with the increasing catchment runoff coefficient. Based on this finding, an approach called the discharge amplification method is proposed to enhance the simulation efficiency of TOPMODEL in rainfall–runoff modelling in catchments with a low runoff coefficient. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

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 analysis based on competing linear impulse responses: decomposition of rainfall‐runoff processes

Many recent studies have successfully used neural networks for non‐linear rainfall‐runoff modelling. Due to fundamental limitation of linear structures, approaches employing linear models have been generally considered inferior to the neural network approaches in this area. However, the authors believe that with an appropriate extension, the concept of linear impulse responses can be a viable tool since it enables one to understand underlying dynamics of rainfall‐runoff processes. In this paper, the use of competing impulse responses for rainfall‐runoff analysis is proposed. The proposed method is based on the switch over of competing linear impulse‐responses, each of which satisfies the constraints of non‐negativity and uni‐modality. The computational analyses performed for the rainfall‐runoff data in the Seolma‐Chun experimental basin, Korea showed that the proposed method can yield promising results. Considering the basin characteristics as well as the results from this study, it may be concluded that three impulse responses are enough for rainfall‐runoff analysis. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Multi‐site calibration using a grid‐based event rainfall–runoff model: a case study of the upstream areas of the Nakdong River basin in Korea

Most runoff analyses using a grid‐based distributed model use one parameter group calibrated at the outlet of a watershed, instead of dividing the watershed into subwatersheds. Significant differences between the observed value and the simulation result of the subwatersheds can occur if just one parameter group is used in all subwatersheds that have different hydrological characteristics from each other. Therefore, to improve the simulation results of the subwatersheds within a watershed, a model calibrated at every subwatershed needs to be used to reflect the characteristics of each subwatershed. In this study, different parameter groups were set up for one or two sites using a distributed model, the GRM (Grid based Rainfall‐runoff Model), and the evaluations were based on the results of rainfall–runoff analysis, which uses a multi‐site calibration (MSC) technique to calibrate the model at the outlet of each site. The Hyangseok watershed in Naeseong River, which is a tributary of Nakdong River in Korea, was chosen as the study area. The watershed was divided into five subwatersheds each with a subwatershed outlet that was applied to the calibration sites . The MSC was applied for five cases. When a site was added for calibration in a watershed, the runoff simulation showed better results than the calibration of only one site at the most downstream area of the watershed. The MSC approach could improve the simulation results on the calibrated sites and even on the non‐calibrated sites, and the effect of MSC was improved when the calibrated site was closer to the runoff site. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

An improved AMC‐coupled runoff curve number model

In the Soil Conservation Service Curve Number (SCS‐CN) method, the three levels of antecedent moisture condition (AMC) permit unreasonable sudden jumps in curve numbers, which result into corresponding jumps in the estimated runoff. A few recently developed SCS‐CN‐based models obviate this problem, yet they have several limitations. In this study, such a model incorporating a continuous function for antecedent moisture has been presented. It has several advantages over the other existing SCS‐CN‐based models. Its application to a large dataset from US watersheds showed to perform better than the existing SCS‐CN method and the others based on it. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessing rainfall–runoff models for the management of Lake Verbano

Growing human pressure and potential change in precipitation pattern induced by climate change require a more efficient and sustainable use of water resources. Hydrological models can provide a fundamental contribution to this purpose, especially as increasing availability of meteorological data and forecast allows for more accurate runoff predictions. In this article, two models are presented for describing the flow formation process in a sub‐alpine catchment: a distributed parameter, physically based model, and a lumped parameter, empirical model. The scope is to compare the two modelling approaches and to assess the impact of hydrometeorological information, either observations or forecast, on water resources management. This is carried out by simulating the real‐time management of the regulated lake that drains the catchment, using the inflow predictions provided by the two models. Copyright © 2010 John Wiley & Sons, Ltd.     

A semi‐distributed parallel‐type linear reservoir rainfall‐runoff model and its application in Taiwan

The main purpose of this paper is to introduce a semi‐distributed parallel surface rainfall‐runoff conceptual model. In this paper, a general solution of the instantaneous unit hydrograph (IUH) has been derived successfully for N linearly connected reservoirs, each having a different storage constant. The solution is a function of geomorphologic parameters, meteorologic factors and roughness coefficients. The model also takes into account the hydrologic response which is influenced by outflow downstream of a reservoir. For calibration, the shuffled complex evolution (SCE) algorithm is used to search for the global optimal parameters of the model. Because of the parallel structure, the mean roughness parameter of the channel becomes a “conceptual” parameter without a real physical meaning. To evaluate the adaptability of the model adopted, three watersheds around the city of Taipei in Taiwan were chosen to test the effectiveness of the model. The study provides an appropriate rainfall‐runoff model for planning flood mitigation in Taiwan. Copyright © 1999 John Wiley & Sons, Ltd.     

Error analysis for the evaluation of model performance: rainfall–runoff event time series data

This paper provides a procedure for evaluating model performance where model predictions and observations are given as time series data. The procedure focuses on the analysis of error time series by graphing them, summarizing them, and predicting their variability through available information (recalibration). We analysed two rainfall–runoff events from the R‐5 data set, and evaluated 12 distinct model simulation scenarios for these events, of which 10 were conducted with the quasi‐physically‐based rainfall–runoff model (QPBRRM) and two with the integrated hydrology model (InHM). The QPBRRM simulation scenarios differ in their representation of saturated hydraulic conductivity. Two InHM simulation scenarios differ with respect to the inclusion of the roads at R‐5. The two models, QPBRRM and InHM, differ strongly in the complexity and number of processes included. For all model simulations we found that errors could be predicted fairly well to very well, based on model output, or based on smooth functions of lagged rainfall data. The errors remaining after recalibration are much more alike in terms of variability than those without recalibration. In this paper, recalibration is not meant to fix models, but merely as a diagnostic tool that exhibits the magnitude and direction of model errors and indicates whether these model errors are related to model inputs such as rainfall. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of rainfall uncertainty on the performance of physically based rainfall–runoff models

This paper analyses the effect of rain data uncertainty on the performance of two hydrological models with different spatial structures: a semidistributed and a fully distributed model. The study is performed on a small catchment of 19.6 km² located in the north‐west of Spain, where the arrival of low pressure fronts from the Atlantic Ocean causes highly variable rainfall events. The rainfall fields in this catchment during a series of storm events are estimated using rainfall point measurements. The uncertainty of the estimated fields is quantified using a conditional simulation technique. Discharge and rain data, including the uncertainty of the estimated rainfall fields, are then used to calibrate and validate both hydrological models following the generalized likelihood uncertainty estimation (GLUE) methodology. In the storm events analysed, the two models show similar performance. In all cases, results show that the calibrated distribution of the input parameters narrows when the rain uncertainty is included in the analysis. Otherwise, when rain uncertainty is not considered, the calibration of the input parameters must account for all uncertainty in the rainfall–runoff transformation process. Also, in both models, the uncertainty of the predicted discharges increase in similar magnitude when the uncertainty of rainfall input increase.     

Using recurrent neural networks to reconstruct rainfall‐runoff processes

Many novel techniques for reconstructing rainfall‐runoff processes require hydrometeorologic and geomorphologic information for modelling. However, certain information is not always measurable. In this paper, we employ a special recurrent neural network to reconstruct the rainfall‐runoff process by using collected rainfall data. In addition, we propose an indirect system identification to overcome the drawback of a traditional, time‐consuming trial‐and‐error search. The indirect system identification is an efficient method to recognize the structure of a recurrent neural network. The unit hydrograph can be derived directly from the weights of the network due to a state‐space form embedded in the recurrent neural network. This improves the link between the weights of the network and the physical concepts that most neural networks fail to connect. The case studies of 41 events from 1966 to 1997 have been implemented in Taiwan's Wu‐Tu watershed, where the runoff path‐lines are short and steep. Two recurrent neural networks and one state‐space model are utilized to simulate the rainfall‐runoff processes for comparison. The results are validated by four criteria: coefficient of efficiency; peak discharge error; time to peak arrival error; total discharge volume error. The resulting data from the recurrent neural network reveal that the neural network proposed herein is appropriate for hydrological systems. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

Using high resolution tracer data to constrain water storage,flux and age estimates in a spatially distributed rainfall‐runoff model

Models simulating stream flow and conservative tracers can provide a representation of flow paths, storage distributions and mixing processes that is advantageous for many predictive purposes. Compared with models that only simulate stream flow, tracer data can be used to investigate the internal consistency of model behaviour and to gain insight into model performance. Here, we examine the strengths and weaknesses of a data‐driven, spatially distributed tracer‐aided rainfall‐runoff model. The model structure allowed us to assess the influence of landscape characteristics on the routing and mixing of water and tracers. The model was applied to a site in the Scottish Highlands with a unique tracer data set; ~4 years of daily isotope ratios in stream water and precipitation were available, as well as 2 years of weekly soil and ground water isotopes. The model structure was based on an empirically based, lumped tracer‐aided model previously developed for the catchment. The best model runs were selected from Monte Carlo simulations based on dual calibration criteria using objective functions for both stream isotopes and discharge at the outlet. Model performance for these criteria was reasonable (Nash–Sutcliffe efficiencies for discharge and isotope ratios were ~0.4–0.6). The model could generally reproduce the variable isotope signals in the soils of the steeper hill slopes where storage was low, and damped isotope responses in valley bottom cells with high storage. The model also allowed us to estimate the age distributions of internal stores, water fluxes and stream flow. Average stream water age was ~1.6 years, integrating older groundwater in the valley bottom and dynamic younger soil waters. By tracking water ages and simulating isotopes, the model captured the changes in connectivity driven by distributed storage dynamics. This has substantially improved the representation of spatio‐temporal process dynamics and gives a more robust framework for projecting environmental change impacts. Copyright © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

Testing a spatially distributed tracer‐aided runoff model in a snow‐influenced catchment: Effects of multicriteria calibration on streamwater ages

Integrating stable isotope tracers into rainfall‐runoff models allows investigation of water partitioning and direct estimation of travel times and water ages. Tracer data have valuable information content that can be used to constrain models and, in integration with hydrometric observations, test the conceptualization of catchment processes in model structure and parameterization. There is great potential in using tracer‐aided modelling in snow‐influenced catchments to improve understanding of these catchments' dynamics and sensitivity to environmental change. We used the spatially distributed tracer‐aided rainfall‐runoff (STARR) model to simulate the interactions between water storage, flux, and isotope dynamics in a snow‐influenced, long‐term monitored catchment in Ontario, Canada. Multiple realizations of the model were achieved using a combination of single and multiple objectives as calibration targets. Although good simulations of hydrometric targets such as discharge and snow water equivalent could be achieved by local calibration alone, adequate capture of the stream isotope dynamics was predicated on the inclusion of isotope data in the calibration. Parameter sensitivity was highest, and most local, for single calibration targets. With multiple calibration targets, key sensitive parameters were still identifiable in snow and runoff generation routines. Water ages derived from flux tracking subroutines in the model indicated a catchment where runoff is dominated by younger waters, particularly during spring snowmelt. However, resulting water ages were most sensitive to the partitioning of runoff sources from soil and groundwater sources, which was most realistically achieved when isotopes were included in the calibration. Given the paucity of studies where hydrological models explicitly incorporate tracers in snow‐influenced regions, this study using STARR is an important contribution to satisfactorily simulating snowpack dynamics and runoff generation processes, while simultaneously capturing stable isotope variability in snow‐influenced catchments.     

Comparing theoretically supported rainfall‐runoff erosivity factors at the Sparacia (South Italy) experimental site

Interpreting rainfall‐runoff erosivity by a process‐oriented scheme allows to conjugate the physical approach to soil loss estimate with the empirical one. Including the effect of runoff in the model permits to distinguish between detachment and transport in the soil erosion process. In this paper, at first, a general definition of the rainfall‐runoff erosivity factor REF_e including the power of both event runoff coefficient Q_R and event rainfall erosivity index EI₃₀ of the Universal Soil Loss Equation (USLE) is proposed. The REF_e factor is applicable to all USLE‐based models (USLE, Modified USLE [USLE‐M] and Modified USLE‐M [USLE‐MM]) and it allows to distinguish between purely empirical models (e.g., Modified USLE‐M [USLE‐MM]) and those supported by applying theoretical dimensional analysis and self‐similarity to Wischmeier and Smith scheme. This last model category includes USLE, USLE‐M, and a new model, named USLE‐M based (USLE‐MB), that uses a rainfall‐runoff erosivity factor in which a power of runoff coefficient multiplies EI₃₀. Using the database of Sparacia experimental site, the USLE‐MB is parameterized and a comparison with soil loss data is carried out. The developed analysis shows that USLE‐MB (characterized by a Nash–Sutcliffe Efficiency Index NSEI equal to 0.73 and a root mean square error RMSE = 11.7 Mg ha^?1) has very similar soil loss estimate performances as compared with the USLE‐M (NSEI = 0.72 and RMSE = 12.0 Mg ha^?1). However, the USLE‐MB yields a maximum discrepancy factor between predicted and measured soil loss values (176) that is much lower than that of USLE‐M (291). In conclusion, the USLE‐MB should be preferred in the context of theoretically supported USLE type models.