Uncertainty analysis of flood inundation modelling using GLUE with surrogate models in stochastic sampling

A generalized likelihood uncertainty estimation (GLUE) method incorporating moving least squares (MLS) with entropy for stochastic sampling (denoted as GLUE‐MLS‐E) was proposed for uncertainty analysis of flood inundation modelling. The MLS with entropy (MLS‐E) was established according to the pairs of parameters/likelihoods generated from a limited number of direct model executions. It was then applied to approximate the model evaluation to facilitate the target sample acceptance of GLUE during the Monte‐Carlo‐based stochastic simulation process. The results from a case study showed that the proposed GLUE‐MLS‐E method had a comparable performance as GLUE in terms of posterior parameter estimation and predicted confidence intervals; however, it could significantly reduce the computational cost. A comparison to other surrogate models, including MLS, quadratic response surface and artificial neural networks (ANN), revealed that the MLS‐E outperformed others in light of both the predicted confidence interval and the most likely value of water depths. ANN was shown to be a viable alternative, which performed slightly poorer than MLS‐E. The proposed surrogate method in stochastic sampling is of practical significance in computationally expensive problems like flood risk analysis, real‐time forecasting, and simulation‐based engineering design, and has a general applicability in many other numerical simulation fields that requires extensive efforts in uncertainty assessment. Copyright © 2014 John Wiley & Sons, Ltd.     

Applying ANN emulators in uncertainty assessment of flood inundation modelling: a comparison of two surrogate schemes

Abstract

A generalized likelihood uncertainty estimation (GLUE) framework coupling with artificial neural network (ANN) models in two surrogate schemes (i.e. GAE-S1 and GAE-S2) was proposed to improve the efficiency of uncertainty assessment in flood inundation modelling. The GAE-S1 scheme was to construct an ANN to approximate the relationship between model likelihoods and uncertain parameters for facilitating sample acceptance/rejection instead of running the numerical model directly; thus, it could speed up the Monte Carlo simulation in stochastic sampling. The GAE-S2 scheme was to establish independent ANN models for water depth predictions to emulate the numerical models; it could facilitate efficient uncertainty analysis without additional model runs for locations concerned under various scenarios. The results from a study case showed that both GAE-S1 and GAE-S2 had comparable performances to GLUE in terms of estimation of posterior parameters, prediction intervals of water depth, and probabilistic inundation maps, but with reduced computational requirements. The results also revealed that GAE-S1 possessed a slightly better performance in accuracy (referencing to GLUE) than GAE-S2, but a lower flexibility in application. This study shed some light on how to apply different surrogate schemes in using numerical models for uncertainty assessment, and could help decision makers in choosing cost-effective ways of conducting flood risk analysis.     

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.     

Application of a data‐based mechanistic modelling (DBM) approach for predicting runoff generation in semi‐arid regions

This paper addresses the application of a data‐based mechanistic (DBM) modelling approach using transfer function models (TFMs) with non‐linear rainfall filtering to predict runoff generation from a semi‐arid catchment (795 km²) in Tanzania. With DBM modelling, time series of rainfall and streamflow were allowed to suggest an appropriate model structure compatible with the data available. The model structures were evaluated by looking at how well the model fitted the data, and how well the parameters of the model were estimated. The results indicated that a parallel model structure is appropriate with a proportion of the runoff being routed through a fast flow pathway and the remainder through a slow flow pathway. Finally, the study employed a Generalized Likelihood Uncertainty Estimation (GLUE) methodology to evaluate the parameter sensitivity and predictive uncertainty based on the feasible parameter ranges chosen from the initial analysis of recession curves and calibration of the TFM. Results showed that parameters that control the slow flow pathway are relatively more sensitive than those that control the fast flow pathway of the hydrograph. Within the GLUE framework, it was found that multiple acceptable parameter sets give a range of predictions. This was found to be an advantage, since it allows the possibility of assessing the uncertainty in predictions as conditioned on the calibration data and then using that uncertainty as part of the decision‐making process arising from any rainfall‐runoff modelling project. Copyright © 2001 John Wiley & Sons, Ltd.     

Assessing the uncertainty in distributed model predictions using observed binary pattern information within GLUE

In this paper we extend the generalized likelihood uncertainty estimation (GLUE) technique to estimate spatially distributed uncertainty in models conditioned against binary pattern data contained in flood inundation maps. Untransformed binary pattern data already have been used within GLUE to estimate domain‐averaged (zero‐dimensional) likelihoods, yet the pattern information embedded within such sources has not been used to estimate distributed uncertainty. Where pattern information has been used to map distributed uncertainty it has been transformed into a continuous function prior to use, which may introduce additional errors. To solve this problem we use here ‘raw’ binary pattern data to define a zero‐dimensional global performance measure for each simulation in a Monte Carlo ensemble. Thereafter, for each pixel of the distributed model we evaluate the probability that this pixel was inundated. This probability is then weighted by the measure of global model performance, thus taking into account how well a given parameter set performs overall. The result is a distributed uncertainty measure mapped over real space. The advantage of the approach is that it both captures distributed uncertainty and contains information on global likelihood that can be used to condition predictions of further events for which observed data are not available. The technique is applied to the problem of flood inundation prediction at two test sites representing different hydrodynamic conditions. In both cases, the method reveals the spatial structure in simulation uncertainty and simultaneously enables mapping of flood probability predicted by the model. Spatially distributed uncertainty analysis is shown to contain information over and above that available from global performance measures. Overall, the paper highlights the different types of information that may be obtained from mappings of model uncertainty over real and n‐dimensional parameter spaces. Copyright © 2002 John Wiley & Sons, Ltd.     

Comparison of parameter uncertainty analysis techniques for a TOPMODEL application

Parameter uncertainty in hydrologic modeling is crucial to the flood simulation and forecasting. The Bayesian approach allows one to estimate parameters according to prior expert knowledge as well as observational data about model parameter values. This study assesses the performance of two popular uncertainty analysis (UA) techniques, i.e., generalized likelihood uncertainty estimation (GLUE) and Bayesian method implemented with the Markov chain Monte Carlo sampling algorithm, in evaluating model parameter uncertainty in flood simulations. These two methods were applied to the semi-distributed Topographic hydrologic model (TOPMODEL) that includes five parameters. A case study was carried out for a small humid catchment in the southeastern China. The performance assessment of the GLUE and Bayesian methods were conducted with advanced tools suited for probabilistic simulations of continuous variables such as streamflow. Graphical tools and scalar metrics were used to test several attributes of the simulation quality of selected flood events: deterministic accuracy and the accuracy of 95 % prediction probability uncertainty band (95PPU). Sensitivity analysis was conducted to identify sensitive parameters that largely affect the model output results. Subsequently, the GLUE and Bayesian methods were used to analyze the uncertainty of sensitive parameters and further to produce their posterior distributions. Based on their posterior parameter samples, TOPMODEL’s simulations and the corresponding UA results were conducted. Results show that the form of exponential decline in conductivity and the overland flow routing velocity were sensitive parameters in TOPMODEL in our case. Small changes in these two parameters would lead to large differences in flood simulation results. Results also suggest that, for both UA techniques, most of streamflow observations were bracketed by 95PPU with the containing ratio value larger than 80 %. In comparison, GLUE gave narrower prediction uncertainty bands than the Bayesian method. It was found that the mode estimates of parameter posterior distributions are suitable to result in better performance of deterministic outputs than the 50 % percentiles for both the GLUE and Bayesian analyses. In addition, the simulation results calibrated with Rosenbrock optimization algorithm show a better agreement with the observations than the UA’s 50 % percentiles but slightly worse than the hydrographs from the mode estimates. The results clearly emphasize the importance of using model uncertainty diagnostic approaches in flood simulations.     

Generalized likelihood uncertainty estimation (GLUE) using adaptive Markov Chain Monte Carlo sampling

In the last few decades hydrologists have made tremendous progress in using dynamic simulation models for the analysis and understanding of hydrologic systems. However, predictions with these models are often deterministic and as such they focus on the most probable forecast, without an explicit estimate of the associated uncertainty. This uncertainty arises from incomplete process representation, uncertainty in initial conditions, input, output and parameter error. The generalized likelihood uncertainty estimation (GLUE) framework was one of the first attempts to represent prediction uncertainty within the context of Monte Carlo (MC) analysis coupled with Bayesian estimation and propagation of uncertainty. Because of its flexibility, ease of implementation and its suitability for parallel implementation on distributed computer systems, the GLUE method has been used in a wide variety of applications. However, the MC based sampling strategy of the prior parameter space typically utilized in GLUE is not particularly efficient in finding behavioral simulations. This becomes especially problematic for high-dimensional parameter estimation problems, and in the case of complex simulation models that require significant computational time to run and produce the desired output. In this paper we improve the computational efficiency of GLUE by sampling the prior parameter space using an adaptive Markov Chain Monte Carlo scheme (the Shuffled Complex Evolution Metropolis (SCEM-UA) algorithm). Moreover, we propose an alternative strategy to determine the value of the cutoff threshold based on the appropriate coverage of the resulting uncertainty bounds. We demonstrate the superiority of this revised GLUE method with three different conceptual watershed models of increasing complexity, using both synthetic and real-world streamflow data from two catchments with different hydrologic regimes.     

The role of hydrological modelling uncertainties in climate change impact assessments of Irish river catchments

This study attempts to assess the uncertainty in the hydrological impacts of climate change using a multi-model approach combining multiple emission scenarios, GCMs and conceptual rainfall-runoff models to quantify uncertainty in future impacts at the catchment scale. The uncertainties associated with hydrological models have traditionally been given less attention in impact assessments until relatively recently. In order to examine the role of hydrological model uncertainty (parameter and structural uncertainty) in climate change impact studies a multi-model approach based on the Generalised Likelihood Uncertainty Estimation (GLUE) and Bayesian Model Averaging (BMA) methods is presented. Six sets of regionalised climate scenarios derived from three GCMs, two emission scenarios, and four conceptual hydrological models were used within the GLUE framework to define the uncertainty envelop for future estimates of stream flow, while the GLUE output is also post processed using BMA, where the probability density function from each model at any given time is modelled by a gamma distribution with heteroscedastic variance. The investigation on four Irish catchments shows that the role of hydrological model uncertainty is remarkably high and should therefore be routinely considered in impact studies. Although, the GLUE and BMA approaches used here differ fundamentally in their underlying philosophy and representation of error, both methods show comparable performance in terms of ensemble spread and predictive coverage. Moreover, the median prediction for future stream flow shows progressive increases of winter discharge and progressive decreases in summer discharge over the coming century.     

A Bayesian method for missing rainfall estimation using a conceptual rainfall–runoff model

The estimation of missing rainfall data is an important problem for data analysis and modelling studies in hydrology. This paper develops a Bayesian method to address missing rainfall estimation from runoff measurements based on a pre-calibrated conceptual rainfall–runoff model. The Bayesian method assigns posterior probability of rainfall estimates proportional to the likelihood function of measured runoff flows and prior rainfall information, which is presented by uniform distributions in the absence of rainfall data. The likelihood function of measured runoff can be determined via the test of different residual error models in the calibration phase. The application of this method to a French urban catchment indicates that the proposed Bayesian method is able to assess missing rainfall and its uncertainty based only on runoff measurements, which provides an alternative to the reverse model for missing rainfall estimates.     

Numerical investigation of the impact of uncertainties in satellite rainfall estimation and land surface model parameters on simulation of soil moisture

This study aims at evaluating the uncertainty in the prediction of soil moisture (1D, vertical column) from an offline land surface model (LSM) forced by hydro-meteorological and radiation data. We focus on two types of uncertainty: an input error due to satellite rainfall retrieval uncertainty, and, LSM soil-parametric error. The study is facilitated by in situ and remotely sensed data-driven (precipitation, radiation, soil moisture) simulation experiments comprising a LSM and stochastic models for error characterization. The parametric uncertainty is represented by the generalized likelihood uncertainty estimation (GLUE) technique, which models the parameter non-uniqueness against direct observations. Half-hourly infra-red (IR) sensor retrievals were used as satellite rainfall estimates. The IR rain retrieval uncertainty is characterized on the basis of a satellite rainfall error model (SREM). The combined uncertainty (i.e., SREM + GLUE) is compared with the partial assessment of uncertainty. It is found that precipitation (IR) error alone may explain moderate to low proportion of the soil moisture simulation uncertainty, depending on the level of model accuracy—50–60% for high model accuracy, and 20–30% for low model accuracy. Comparisons on the basis of two different sites also yielded an increase (50–100%) in soil moisture prediction uncertainty for the more vegetated site. This study exemplified the need for detailed investigations of the rainfall retrieval-modeling parameter error interaction within a comprehensive space-time stochastic framework for achieving optimal integration of satellite rain retrievals in land data assimilation systems.     

Assessing the results of scenarios of climate and land use changes on the hydrology of an Italian catchment: modelling study

Hydrological models are recognized as valid scientific tools to study water quantity and quality and provide support for the integrated management and planning of water resources at different scales. In common with many catchments in the Mediterranean, the study catchment has many problems such as the increasing gap between water demand and supply, water quality deterioration, scarcity of available data, lack of measurements and specific information. The application of hydrological models to investigate hydrological processes in this type of catchments is of particular relevance for water planning strategies to address the possible impact of climate and land use changes on water resources. The distributed catchment scale model (DiCaSM) was selected to study the impact of climate and land use changes on the hydrological cycle and the water balance components in the Apulia region, southern Italy, specifically in the Candelaro catchment (1780 km²). The results obtained from this investigation proved the ability of DiCaSM to quantify the different components of the catchment water balance and to successfully simulate the stream flows. In addition, the model was run with the climate change scenarios for southern Italy, i.e. reduced winter rainfall by 5–10%, reduced summer rainfall by 15–20%, winter temperature rise by 1·25–1·5 °C and summer temperature rise by 1·5–1·75 °C. The results indicated that by 2050, groundwater recharge in the Candelaro catchment would decrease by 21–31% and stream flows by 16–23%. The model results also showed that the projected durum wheat yield up to 2050 is likely to decrease between 2·2% and 10·4% due to the future reduction in rainfall and increase in temperature. In the current study, the reliability of the DiCaSM was assessed when applied to the Candelaro catchment; those parameters that may cause uncertainty in model output were investigated using a generalized likelihood uncertainty estimation (GLUE) methodology. The results showed that DiCaSM provided a small level of uncertainty and subsequently, a higher confidence level. Copyright © 2010 John Wiley & Sons, Ltd.     

Quantifying and predicting the accuracy of radar-based quantitative precipitation forecasts

On the basis that hydrological users need to know the forecast uncertainty at the time that the forecast is issued, we computed distributions of radar rainfall forecast uncertainty as a function of forecast lead time, basin size, and forecasted rainfall intensity using data from the US 3-D National Mosaic of radar data. We document how exceptional forecasts such as those of heavy rainfall are generally biased. Since forecast uncertainty is also weather dependent, we tried to find good predictors to help either reduce the forecast uncertainty or better define it. These predictors were based either on characteristics of the current precipitation field or on the performance of the nowcast in the immediate past. The value of some predictors, especially those based on the properties of large-scale rainfall patterns, was significant though modest, the predictors being generally more skillful at characterizing forecast uncertainty than at improving forecast accuracy.     

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.     

Parameter identification and uncertainty analysis for simulating streamflow in a river basin of Eastern India

Uncertainty is inherent in modelling studies. However, the quantification of uncertainties associated with a model is a challenging task, and hence, such studies are somewhat limited. As distributed or semi‐distributed hydrological models are being increasingly used these days to simulate hydrological processes, it is vital that these models should be equipped with robust calibration and uncertainty analysis techniques. The goal of the present study was to calibrate and validate the Soil and Water Assessment Tool (SWAT) model for simulating streamflow in a river basin of Eastern India, and to evaluate the performance of salient optimization techniques in quantifying uncertainties. The SWAT model for the study basin was developed and calibrated using Parameter Solution (ParaSol), Sequential Uncertainty Fitting Algorithm (SUFI‐2) and Generalized Likelihood Uncertainty Estimation (GLUE) optimization techniques. The daily observed streamflow data from 1998 to 2003 were used for model calibration, and those for 2004–2005 were used for model validation. Modelling results indicated that all the three techniques invariably yield better results for the monthly time step than for the daily time step during both calibration and validation. The model performances for the daily streamflow simulation using ParaSol and SUFI‐2 during calibration are reasonably good with a Nash–Sutcliffe efficiency and mean absolute error (MAE) of 0.88 and 9.70 m³/s for ParaSol, and 0.86 and 10.07 m³/s for SUFI‐2, respectively. The simulation results of GLUE revealed that the model simulates daily streamflow during calibration with the highest accuracy in the case of GLUE (R² = 0.88, MAE = 9.56 m³/s and root mean square error = 19.70 m³/s). The results of uncertainty analyses by SUFI‐2 and GLUE were compared in terms of parameter uncertainty. It was found that SUFI‐2 is capable of estimating uncertainties in complex hydrological models like SWAT, but it warrants sound knowledge of the parameters and their effects on the model output. On the other hand, GLUE predicts more reliable uncertainty ranges (R‐factor = 0.52 for daily calibration and 0.48 for validation) compared to SUFI‐2 (R‐factor = 0.59 for daily calibration and 0.55 for validation), though it is computationally demanding. Although both SUFI‐2 and GLUE appear to be promising techniques for the uncertainty analysis of modelling results, more and more studies in this direction are required under varying agro‐climatic conditions for assessing their generic capability. Copyright © 2015 John Wiley & Sons, Ltd.     

Fuzzy computing based rainfall–runoff model for real time flood forecasting

This paper analyses the skills of fuzzy computing based rainfall–runoff model in real time flood forecasting. The potential of fuzzy computing has been demonstrated by developing a model for forecasting the river flow of Narmada basin in India. This work has demonstrated that fuzzy models can take advantage of their capability to simulate the unknown relationships between a set of relevant hydrological data such as rainfall and river flow. Many combinations of input variables were presented to the model with varying structures as a sensitivity study to verify the conclusions about the coherence between precipitation, upstream runoff and total watershed runoff. The most appropriate set of input variables was determined, and the study suggests that the river flow of Narmada behaves more like an autoregressive process. As the precipitation is weighted only a little by the model, the last time‐steps of measured runoff are dominating the forecast. Thus a forecast based on expected rainfall becomes very inaccurate. Although good results for one‐step‐ahead forecasts are received, the accuracy deteriorates as the lead time increases. Using the one‐step‐ahead forecast model recursively to predict flows at higher lead time, however, produces better results as opposed to different independent fuzzy models to forecast flows at various lead times. Copyright © 2004 John Wiley & Sons, Ltd.     

Climate models and hydrological parameter uncertainties in climate change impacts on monthly runoff and daily flow duration curve of a Mediterranean catchment

ABSTRACT

Climate models and hydrological parameter uncertainties were quantified and compared while assessing climate change impacts on monthly runoff and daily flow duration curve (FDC) in a Mediterranean catchment. Simulations of the Soil and Water Assessment Tool (SWAT) model using an ensemble of behavioural parameter sets derived from the Generalized Likelihood Uncertainty Estimation (GLUE) method were approximated by feed-forward artificial neural networks (FF-NN). Then, outputs of climate models were used as inputs to the FF-NN models. Subsequently, projected changes in runoff and FDC were calculated and their associated uncertainty was partitioned into climate model and hydrological parameter uncertainties. Runoff and daily discharge of the Chiba catchment were expected to decrease in response to drier and warmer climatic conditions in the 2050s. For both hydrological indicators, uncertainty magnitude increased when moving from dry to wet periods. The decomposition of uncertainty demonstrated that climate model uncertainty dominated hydrological parameter uncertainty in wet periods, whereas in dry periods hydrological parametric uncertainty became more important.

Editor M.C. Acreman; Associate editor S. Kanae     

Forecast-skill-based simulation of streamflow forecasts

Streamflow forecasts are updated periodically in real time, thereby facilitating forecast evolution. This study proposes a forecast-skill-based model of forecast evolution that is able to simulate dynamically updated streamflow forecasts. The proposed model applies stochastic models that deal with streamflow variability to generate streamflow scenarios, which represent cases without forecast skill of future streamflow. The model then employs a coefficient of prediction to determine forecast skill and to quantify the streamflow variability ratio explained by the forecast. By updating the coefficients of prediction periodically, the model efficiently captures the evolution of streamflow forecast. Simulated forecast uncertainty increases with increasing lead time; and simulated uncertainty during a specific future period decreases over time. We combine the statistical model with an optimization model and design a hypothetical case study of reservoir operation. The results indicate the significance of forecast skill in forecast-based reservoir operation. Shortage index reduces as forecast skill increases and ensemble forecast outperforms deterministic forecast at a similar forecast skill level. Moreover, an effective forecast horizon exists beyond which more forecast information does not contribute to reservoir operation and higher forecast skill results in longer effective forecast horizon. The results illustrate that the statistical model is efficient in simulating forecast evolution and facilitates analysis of forecast-based decision making.     

Updating Real-Time Flood Forecasting Using a Fuzzy Rule-Based Model/Mise à Jour de Prévision de Crue en Temps Réel Grâce à un Modèle à Base de Règles Floues

Abstract

An updating technique is a tool to update the forecasts of mathematical flood forecasting model based on data observed in real time, and is an important element in a flood forecasting model. An error prediction model based on a fuzzy rule-based method was proposed as the updating technique in this work to improve one- to four-hour-ahead flood forecasts by a model that is composed of the grey rainfall model, the grey rainfall—runoff model and the modified Muskingum flow routing model. The coefficient of efficiency with respect to a benchmark is applied to test the applicability of the proposed fuzzy rule-based method. The analysis reveals that the fuzzy rule-based method can improve flood forecasts one to four hours ahead. The proposed updating technique can mitigate the problem of the phase lag in forecast hydrographs, and especially in forecast hydrographs with longer lead times.     

A computationally efficient method for uncertainty analysis of SWAT model simulations

The physically based distributed hydrological models are ideal for hydrological simulations; however most of such models do not use the basic equations pertaining to mass, energy and momentum conservation, to represent the physics of the process. This is plausibly due to the lack of complete understanding of the hydrological process. The soil and water assessment tool (SWAT) is one such widely accepted semi-distributed, conceptual hydrological model used for water resources planning. However, the over-parameterization, difficulty in its calibration process and the uncertainty associated with predictions make its applications skeptical. This study considers assessing the predictive uncertainty associated with distributed hydrological models. The existing methods for uncertainty estimation demand high computational time and therefore make them challenging to apply on complex hydrological models. The proposed approach employs the concepts of generalized likelihood uncertainty estimation (GLUE) in an iterative procedure by starting with an assumed prior probability distribution of parameters, and by using mutual information (MI) index for sampling the behavioral parameter set. The distributions are conditioned on the observed information through successive cycles of simulations. During each cycle of simulation, MI is used in conjunction with Markov Chain Monte Carlo procedure to sample the parameter sets so as to increase the number of behavioral sets, which in turn helps reduce the number of cycles/simulations for the analysis. The method is demonstrated through a case study of SWAT model in Illinois River basin in the USA. A comparison of the proposed method with GLUE indicates that the computational requirement of uncertainty analysis is considerably reduced in the proposed approach. It is also noted that the model prediction band, derived using the proposed method, is more effective compared to that derived using the other methods considered in this study.     

Quantification of parametric uncertainty of ANN models with GLUE method for different streamflow dynamics

This study describes the parametric uncertainty of artificial neural networks (ANNs) by employing the generalized likelihood uncertainty estimation (GLUE) method. The ANNs are used to forecast daily streamflow for three sub-basins of the Rhine Basin (East Alpine, Main, and Mosel) having different hydrological and climatological characteristics. We have obtained prior parameter distributions from 5000 ANNs in the training period to capture the parametric uncertainty and subsequently 125,000 correlated parameter sets were generated. These parameter sets were used to quantify the uncertainty in the forecasted streamflow in the testing period using three uncertainty measures: percentage of coverage, average relative length, and average asymmetry degree. The results indicated that the highest uncertainty was obtained for the Mosel sub-basin and the lowest for the East Alpine sub-basin mainly due to hydro-climatic differences between these basins. The prediction results and uncertainty estimates of the proposed methodology were compared to the direct ensemble and bootstrap methods. The GLUE method successfully captured the observed discharges with the generated prediction intervals, especially the peak flows. It was also illustrated that uncertainty bands are sensitive to the selection of the threshold value for the Nash–Sutcliffe efficiency measure used in the GLUE method by employing the Wilcoxon–Mann–Whitney test.