Data time-step dependency of conceptual rainfall—streamflow model parameters: an empirical study with implications for regionalisation

Abstract

Employing 1-, 2-, 4-, 6-, 12- and 24-hourly data sets for two catchments (10.6 and 298 km²) in Wales, the calibrated parameters of a unit hydrograph-based model are shown to change substantially over that range of data time steps. For the smaller basin, each model parameter reaches, or approaches, a stable value as the data time step decreases, providing a straightforward method of estimating time-step independent model parameter values. For the larger basin, the model parameters also reach, or approach, stable values using hourly data, but, for reasons given in the paper, interpretation of the results is more difficult. Model parameter sensitivity analyses are presented that give insights into the relative precision on the parameters for both catchments. The paper discusses the importance of accounting for model parameter data time-step dependency in pursuit of a reduction in the uncertainty associated with estimates of flow in ungauged basins, and suggests that further work along these lines be undertaken using different catchments and models.     

A comparison of three approaches to spatial generalization of rainfall–runoff models

The generalization of the parameters of rainfall–runoff models, to enable application at ungauged sites, is an important and ongoing area of research. This paper compares the performance of three alternative methods of generalization, for two parameter‐sparse conceptual models (PDM and TATE), specifically for use in flood frequency estimation using continuous simulation. Two of the methods are based on fitting regression relationships between catchment properties and calibrated parameter values, using weighted or sequential regression (with weights based on estimates of calibration uncertainty), and the third is based on the use of pooling groups, defined through measures of site‐similarity based on catchment properties. The study uses a relatively large sample of catchments in Britain. For the PDM, the site‐similarity method performs best, but not greatly better than either regression method, so there may be cases where the use of regression would be preferable. For the TATE model, weighted regression performs best (with a very similar level of performance to that of the PDM with site‐similarity), whereas site‐similarity performs worst (due to poor performance for catchments with higher baseflow), indicating that the choice of model and generalization method should not be separated. The use of sequential regression, which was developed to try to allow for parameter interdependence, shows no clear advantage for either model. Other than the poor performance of the TATE model with site‐similarity for catchments with a higher baseflow index, there are no clear relationships between performance of any model/method and catchment type. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessment of spatial transferability of process‐based hydrological model parameters in two neighbouring catchments in the Himalayan Region

Estimating the hydrological regime of ungauged catchments in the Himalayan region is challenging due to a lack of sufficient monitoring stations. In this paper, the spatial transferability of the model parameters of the process‐oriented J2000 hydrological model was investigated in 2 glaciated subcatchments of the Koshi river basin in eastern Nepal. The catchments have a high degree of similarity with respect to their static landscape features. The model was first calibrated (1986–1991) and validated (1992–1997) in the Dudh Koshi subcatchment. The calibrated and validated model parameters were then transferred to the nearby Tamor catchment (2001–2009). Sensitivity and uncertainty analyses were carried out for both subcatchments to discover the sensitivity range of the parameters in the two catchments. The model represented the overall hydrograph well in both subcatchments, including baseflow, rising and falling limbs; however, the peak flows were underestimated. The efficiency results according to both Nash–Sutcliffe (ENS) and the coefficient of determination (r²) were above 0.84 in both catchments (1986–1997 in Dudh Koshi and 2001–2009 in Tamor). The ranking of the parameters in respect to their sensitivity matched well for both catchments while taking ENS and log Nash–Sutcliffe (LNS) efficiencies into account. However, there were some differences in sensitivity to ENS and LNS for moderately and less‐sensitive parameters, although the majority (13 out of 16 for ENS and 16 out of 16 for LNS) had a sensitivity response in a similar range. The generalized uncertainty likelihood estimation results suggest that the parameter uncertainty are most of the time within the range and the ensemble mean matches very good (ENS: 0.84) with observed discharge. The results indicate that transfer of the J2000 parameters to a neighbouring catchment in the Himalayan region with similar physiographic landscape characteristics is viable. This indicates the possibility of applying a calibrated process‐based J2000 model to other ungauged catchments in the Himalayan region, which could provide important insights into the hydrological system dynamics and provide much needed information to support water resources planning and management.     

Uncertainty in index flood modelling due to calibration data sizes

Regional frequency analysis is an important tool in estimating design flood for ungauged catchments. Index flood is an important component in regionalized flood formulas. In the past, many formulas have been developed based on various numbers of calibration catchments (e.g. from less than 20 to several hundred). However, there is a lack of systematic research on the model uncertainties caused by the number of calibration catchments (i.e. what is the minimum number of calibration catchment? and how should we choose the calibration catchments?). This study uses the statistical resampling technique to explore the impact of calibration catchment numbers on the index flood estimation. The study is based on 182 catchments in England and an index flood formula has been developed using the input variable selection technique in the data mining field. The formula has been used to explore the model uncertainty due to a range of calibration catchment numbers (from 15 to 130). It is found that (1) as expected, the more catchments are used in the calibration, the more reliable of the models developed are (i.e. with a narrower band of uncertainty); (2) however, poor models are still possible with a large number of calibration catchments (e.g. 130). In contrast, good models with a small number of calibration catchments are also achievable (with as low as 15 calibration catchments). This indicates that the number of calibration catchments is only one of the factors influencing the model performance. The hydrological community should explore why a smaller calibration data set could produce a better model than a large calibration data set. It is clear from this study that the information content in the calibration data set is equally if not more important than the number of calibration data. Copyright © 2011 John Wiley & Sons, Ltd.     

Can climate variability information constrain a hydrological model for an ungauged Costa Rican catchment?

Long‐term hydrological data are key to understanding catchment behaviour and for decision making within water management and planning. Given the lack of observed data in many regions worldwide, such as Central America, hydrological models are an alternative for reproducing historical streamflow series. Additional types of information—to locally observed discharge—can be used to constrain model parameter uncertainty for ungauged catchments. Given the strong influence that climatic large‐scale processes exert on streamflow variability in the Central American region, we explored the use of climate variability knowledge as process constraints to constrain the simulated discharge uncertainty for a Costa Rican catchment, assumed to be ungauged. To reduce model uncertainty, we first rejected parameter relationships that disagreed with our understanding of the system. Then, based on this reduced parameter space, we applied the climate‐based process constraints at long‐term, inter‐annual, and intra‐annual timescales. In the first step, we reduced the initial number of parameters by 52%, and then, we further reduced the number of parameters by 3% with the climate constraints. Finally, we compared the climate‐based constraints with a constraint based on global maps of low‐flow statistics. This latter constraint proved to be more restrictive than those based on climate variability (further reducing the number of parameters by 66% compared with 3%). Even so, the climate‐based constraints rejected inconsistent model simulations that were not rejected by the low‐flow statistics constraint. When taken all together, the constraints produced constrained simulation uncertainty bands, and the median simulated discharge followed the observed time series to a similar level as an optimized model. All the constraints were found useful in constraining model uncertainty for an—assumed to be—ungauged basin. This shows that our method is promising for modelling long‐term flow data for ungauged catchments on the Pacific side of Central America and that similar methods can be developed for ungauged basins in other regions where climate variability exerts a strong control on streamflow variability.     

Uncertainty based modeling of rainfall-runoff: Combined differential evolution adaptive Metropolis (DREAM) and K-means clustering

Simulation of rainfall-runoff process in urban areas is of great importance considering the consequences and damages of extreme runoff events and floods. The first issue in flood hazard analysis is rainfall simulation. Large scale climate signals have been proved to be effective in rainfall simulation and prediction. In this study, an integrated scheme is developed for rainfall-runoff modeling considering different sources of uncertainty. This scheme includes three main steps of rainfall forecasting, rainfall-runoff simulation and future runoff prediction. In the first step, data driven models are developed and used to forecast rainfall using large scale climate signals as rainfall predictors. Due to high effect of different sources of uncertainty on the output of hydrologic models, in the second step uncertainty associated with input data, model parameters and model structure is incorporated in rainfall-runoff modeling and simulation. Three rainfall-runoff simulation models are developed for consideration of model conceptual (structural) uncertainty in real time runoff forecasting. To analyze the uncertainty of the model structure, streamflows generated by alternative rainfall-runoff models are combined, through developing a weighting method based on K-means clustering. Model parameters and input uncertainty are investigated using an adaptive Markov Chain Monte Carlo method. Finally, calibrated rainfall-runoff models are driven using the forecasted rainfall to predict future runoff for the watershed. The proposed scheme is employed in the case study of the Bronx River watershed, New York City. Results of uncertainty analysis of rainfall-runoff modeling reveal that simultaneous estimation of model parameters and input uncertainty significantly changes the probability distribution of the model parameters. It is also observed that by combining the outputs of the hydrological models using the proposed clustering scheme, the accuracy of runoff simulation in the watershed is remarkably improved up to 50% in comparison to the simulations by the individual models. Results indicate that the developed methodology not only provides reliable tools for rainfall and runoff modeling, but also adequate time for incorporating required mitigation measures in dealing with potentially extreme runoff events and flood hazard. Results of this study can be used in identification of the main factors affecting flood hazard analysis.     

Application of L‐moments and Bayesian inference for low‐flow regionalization in Sefidroud basin,Iran

Reliable estimation of low flows at ungauged catchments is one of the major challenges in water‐resources planning and management. This study aims at providing at‐site and ungauged sites low‐flow frequency analysis using regionalization approach. A two‐stage delineating homogeneous region is proposed in this study. Clustering sites with similar low‐flow L‐moment ratios is initially conducted, and L‐moment‐based discordancy and heterogeneity measures are then used to detect unusual sites. Based on the goodness‐of‐fit test statistic, the best‐fit regional model is identified in each hydrologically homogeneous region. The relationship between mean annual 7‐day minimum flow and hydro‐geomorphic characteristics is also constructed in each homogeneous region associated with the derived regional model for estimating various low‐flow quantiles at ungauged sites. Uncertainty analysis of model parameters and low‐flow estimations is carried out using the Bayesian inference. Applied in Sefidroud basin located in northwestern Iran, two hydrologically homogeneous regions are identified, i.e. the east and west regions. The best‐fit regional model for the east and west regions are generalized logistic and Pearson type III distributions, respectively. The results show that the proposed approach provides reasonably good accuracy for at‐site as well as ungauged‐site frequency analysis. Besides, interval estimations for model parameters and low flows provide uncertainty information, and the results indicate that Bayesian confidence intervals are significantly reduced when comparing with the outcomes of conventional t‐distribution method. Copyright © 2013 John Wiley & Sons, Ltd.     

Using micro-catchment experiments for multi-local scale modelling of nature-based solutions

The Q-natural flood management project has co-developed with the Environment Agency 18 monitored micro-catchments (~1 km²) in Cumbria, UK installing calibrated flumes aimed at quantifying the potential shift in observed flows resulting from a range of nature-based-solutions installed by local organizations. The small-scale reduces the influence of variability characterizing larger catchments that would otherwise mask any such shifts, which we attempt to relate to a shift in model parameters. This paper demonstrates an approach to applying donor-parameter-shifts obtained from modelling two of the paired micro-catchments to a much larger scale, in order to understand the potential for improved distributed modelling of nature-based solutions in the form of additional tree-planting. The models include a rainfall-runoff model, Dynamic Topmodel, and a 2D hydrodynamic model, JFlow, permitting analysis of changes in hillslope processes and channel hydrodynamics resulting from a range of distributed measures designed to emulate natural hydrological processes that evaporate, store or infiltrate flows. We report on attempts to detect shift in hydrological response using one of the paired-micro-catchment moorland versus forestry sites in Lorton using Dynamic Topmodel. A donor-parameter-shift approach is used in a hypothetical experiment to represent new woodland in a much larger catchment, although testing all combinations of spatial planting strategies, responses to multiple-extremes, failure-modes and changes to synchronization becomes intractable to support good decision making. We argue that the problem can be re-framed to use donor-parameter-shifts at multi-local-scale catchments above communities known to be at risk, commensurate with most of the evidence of NbS impacts being effective at the small scale (ca. 10 km²). This might lead to more effective modelling to help catchment managers prioritize those communities-at-risk where there is more evidence that NbS might be effective.     

Machine learning based identification of dominant controls on runoff dynamics

Hydrological models used for flood prediction in ungauged catchments are commonly fitted to regionally transferred data. The key issue of this procedure is to identify hydrologically similar catchments. Therefore, the dominant controls for the process of interest have to be known. In this study, we applied a new machine learning based approach to identify the catchment characteristics that can be used to identify the active processes controlling runoff dynamics. A random forest (RF) regressor has been trained to estimate the drainage velocity parameters of a geomorphologic instantaneous unit hydrograph (GIUH) in ungauged catchments, based on regionally available data. We analyzed the learning procedure of the algorithm and identified preferred donor catchments for each ungauged catchment. Based on the obtained machine learning results from catchment grouping, a classification scheme for drainage network characteristics has been derived. This classification scheme has been applied in a flood forecasting case study. The results demonstrate that the RF could be trained properly with the selected donor catchments to successfully estimate the required GIUH parameters. Moreover, our results showed that drainage network characteristics can be used to identify the influence of geomorphological dispersion on the dynamics of catchment response.     

G) Personalia

Abstract

This paper proposes a framework for identifying the parameters of a lumped routing model in small to medium sized catchments where lateral inflows can be large but poorly defined. In a first step, a priori estimates of the parameters are made based on topography, aerial photographs, flood marks and field surveys. In a second step, runoff data are analysed of reservoir release events and convective events where no rainfall in the direct catchments occurred. In a third step the routing model is calibrated to the results of hydrodynamic models for scenarios of different magnitudes. In a fourth step, these pieces of information are combined, allowing for soft expert judgement to be incorporated. In a fifth step, the routing parameters are fine tuned to observed flood events where lateral inflows are estimated by a rainfall—runoff model. The framework is illustrated by the Kamp flood forecasting system in Austria that has been in operational use since 2006.     

Application of distributed hydrological models for predictions in ungauged basins: a method to quantify predictive uncertainty

Stream flow predictions in ungauged basins are one of the most challenging tasks in surface water hydrology because of nonavailability of data and system heterogeneity. This study proposes a method to quantify stream flow predictive uncertainty of distributed hydrologic models for ungauged basins. The method is based on the concepts of deriving probability distribution of model's sensitive parameters by using measured data from a gauged basin and transferring the distribution to hydrologically similar ungauged basins for stream flow predictions. A Monte Carlo simulation of the hydrologic model using sampled parameter sets with assumed probability distribution is conducted. The posterior probability distributions of the sensitive parameters are then computed using a Bayesian approach. In addition, preselected threshold values of likelihood measure of simulations are employed for sizing the parameter range, which helps reduce the predictive uncertainty. The proposed method is illustrated through two case studies using two hydrologically independent sub‐basins in the Cedar Creek watershed located in Texas, USA, using the Soil and Water Assessment Tool (SWAT) model. The probability distribution of the SWAT parameters is derived from the data from one of the sub‐basins and is applied for simulation in the other sub‐basin considered as pseudo‐ungauged. In order to assess the robustness of the method, the numerical exercise is repeated by reversing the gauged and pseudo‐ungauged basins. The results are subsequently compared with the measured stream flow from the sub‐basins. It is observed that the measured stream flow in the pseudo‐ungauged basin lies well within the estimated confidence band of predicted stream flow. Copyright © 2013 John Wiley & Sons, Ltd.     

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

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

Derivation and verification of empirical catchment response time equations for medium to large catchments in South Africa

Despite uncertainties and errors in measurement, observed peak discharges are the best estimate of the true peak discharge from a catchment. However, in ungauged catchments, the catchment response time is a fundamental input to all methods of estimating peak discharges; hence, errors in estimated catchment response time directly impact on estimated peak discharges. In South Africa, this is particularly the case in ungauged medium to large catchments where practitioners are limited to use empirical methods that were calibrated on small catchments not located in South Africa. The time to peak (T_P), time of concentration (T_C) and lag time (T_L) are internationally the most frequently used catchment response time parameters and are normally estimated using either hydraulic or empirical methods. Almost 95% of all the time parameter estimation methods developed internationally are empirically based. This paper presents the derivation and verification of empirical T_P equations in a pilot scale study using 74 catchments located in four climatologically different regions of South Africa, with catchment areas ranging from 20 km² to 35 000 km². The objective is to develop unique relationships between observed T_P values and key climatological and geomorphological catchment predictor variables in order to estimate catchment T_P values at ungauged catchments. The results show that the derived empirical T_P equation(s) meet the requirement of consistency and ease of application. Independent verification tests confirmed the consistency, while the statistically significant independent predictor variables included in the regressions provide a good estimation of catchment response times and are also easy to determine by practitioners when required for future applications in ungauged catchments. It is recommended that the methodology used in this study should be expanded to other catchments to enable the development of a regional approach to improve estimation of time parameters on a national‐scale. However, such a national‐scale application would not only increase the confidence in using the suggested methodology and equation(s) in South Africa, but also highlights that a similar approach could be adopted internationally. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing the relative importance of parameter and forcing uncertainty and their interactions in conceptual hydrological model simulations

Predictions of river flow dynamics provide vital information for many aspects of water management including water resource planning, climate adaptation, and flood and drought assessments. Many of the subjective choices that modellers make including model and criteria selection can have a significant impact on the magnitude and distribution of the output uncertainty. Hydrological modellers are tasked with understanding and minimising the uncertainty surrounding streamflow predictions before communicating the overall uncertainty to decision makers. Parameter uncertainty in conceptual rainfall-runoff models has been widely investigated, and model structural uncertainty and forcing data have been receiving increasing attention. This study aimed to assess uncertainties in streamflow predictions due to forcing data and the identification of behavioural parameter sets in 31 Irish catchments. By combining stochastic rainfall ensembles and multiple parameter sets for three conceptual rainfall-runoff models, an analysis of variance model was used to decompose the total uncertainty in streamflow simulations into contributions from (i) forcing data, (ii) identification of model parameters and (iii) interactions between the two. The analysis illustrates that, for our subjective choices, hydrological model selection had a greater contribution to overall uncertainty, while performance criteria selection influenced the relative intra-annual uncertainties in streamflow predictions. Uncertainties in streamflow predictions due to the method of determining parameters were relatively lower for wetter catchments, and more evenly distributed throughout the year when the Nash-Sutcliffe Efficiency of logarithmic values of flow (lnNSE) was the evaluation criterion.     

新安江模型在资料匮乏的长江中下游山区中小流域洪水预报应用

水文资料匮乏流域的洪水预报(PUBs)是水文科学与工程中一个尚未解决的重大挑战.中国湿润山区中小流域大多是水文资料匮乏的流域,在此地区进行洪水预报的重要手段之一就是水文模型参数的估计.对基于参数物理意义的估算方法(以下简称物理估算法)及两种区域化方法进行了研究,将其用于新安江模型参数的估算及移植.皖南山区的29个中小流...     

Regional calibration of a watershed model

Abstract

As watershed models become increasingly sophisticated and useful, there is a need to extend their applicability to locations where they cannot be calibrated or validated. A new methodology for the regionalization of a watershed model is introduced and evaluated. The approach involves calibration of a watershed model to many sites in a region, concurrently. Previous research that has sought to relate the parameters of monthly water balance models to physical drainage basin characteristics in a region has met with limited success. Previous studies have taken the two-step approach: (a) estimation of watershed model parameters at each site, followed by (b) attempts to relate model parameters to drainage basin characteristics. Instead of treating these two steps as independent, both steps are implemented concurrently. All watershed models in a region are calibrated simultaneously, with the dual objective of reproducing the behaviour of observed monthly streamflows and, additionally, to obtain good relationships between watershed model parameters and basin characteristics. The approach is evaluated using 33 basins in the southeastern region of the United States by comparing simulations using the regional models for three catchments which were not used to develop the regional regression equations. Although the regional calibration approach led to nearly perfect regional relationships between watershed model parameters and basin characteristics, these “improved” regional relationships did not result in improvements in the ability to model streamflow at ungauged sites. This experiment reveals that improvements in regional relationships between watershed model parameters and basin characteristics will not necessarily lead to improvements in the ability to calibrate a watershed model at an ungauged site.     

Regionalization of constraints on expected watershed response behavior for improved predictions in ungauged basins

Approaches to modeling the continuous hydrologic response of ungauged basins use observable physical characteristics of watersheds to either directly infer values for the parameters of hydrologic models, or to establish regression relationships between watershed structure and model parameters. Both these approaches still have widely discussed limitations, including impacts of model structural uncertainty. In this paper we introduce an alternative, model independent, approach to streamflow prediction in ungauged basins based on empirical evidence of relationships between watershed structure, climate and watershed response behavior. Instead of directly estimating values for model parameters, different hydrologic response behaviors of the watershed, quantified through model independent streamflow indices, are estimated and subsequently regionalized in an uncertainty framework. This results in expected ranges of streamflow indices in ungauged watersheds. A pilot study using 30 UK watersheds shows how this regionalized information can be used to constrain ensemble predictions of any model at ungauged sites. Dominant controlling characteristics were found to be climate (wetness index), watershed topography (slope), and hydrogeology. Main streamflow indices were high pulse count, runoff ratio, and the slope of the flow duration curve. This new approach provided sharp and reliable predictions of continuous streamflow at the ungauged sites tested.     

基于卫星遥感的太湖蓝藻水华时空分布规律认识

由于大尺度水文模型和无资料区水文研究是当前国际水文研究的重点和难点,通过参数区域化方法来估计大尺度区域和无资料区的模型参数值成为了研究的热点之一将HBV模型应用于东江流域及其子流域,采用代理流域法和全局乎均法来估计该区域内无资料流域的模型参数研究表明：HBV模型能较好得用于东江流域径流模拟;交叉检验中,较小的序和ME值对应的参数,其转移效果不一定比较大的R^2和ME值对应的参数转移效果差;全局平均法中,面积权重平均值和泰森多边形插值后平均并不能明显改进子流域算术平均值估计无资料流域的模型参数的模拟结果;两者都能有效用于东江流域无资料流域的参数估计,且效果相差不大。