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.     

Development of probability distributions for urban hydrologic model parameters and a Monte Carlo analysis of model sensitivity

This paper proposes an approach to estimating the uncertainty related to EPA Storm Water Management Model model parameters, percentage routed (P_R) and saturated hydraulic conductivity (K_sat), which are used to calculate stormwater runoff volumes. The methodology proposed in this paper addresses uncertainty through the development of probability distributions for urban hydrologic parameters through extensive calibration to observed flow data in the Philadelphia collection system. The established probability distributions are then applied to the Philadelphia Southeast district model through a Monte Carlo approach to estimate the uncertainty in prediction of combined sewer overflow volumes as related to hydrologic model parameter estimation. Understanding urban hydrology is critical to defining urban water resource problems. A variety of land use types within Philadelphia coupled with a history of cut and fill have resulted in a patchwork of urban fill and native soils. The complexity of urban hydrology can make model parameter estimation and defining model uncertainty a difficult task. The development of probability distributions for hydrologic parameters applied through Monte Carlo simulations provided a significant improvement in estimating model uncertainty over traditional model sensitivity analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Automated calibration applied to watershed‐scale flow simulations

Complexity in simulating the hydrological response in large watersheds over long times has prompted a significant need for procedures for automatic calibration. Such a procedure is implemented in the basin‐scale hydrological model (BSHM), a physically based distributed parameter watershed model. BSHM simulates the most important basin‐scale hydrological processes, such as overland flow, groundwater flow and stream–aquifer interaction in watersheds. Here, the emphasis is on estimating the groundwater parameters with water levels in wells and groundwater baseflows selected as the calibration targets. The best set of parameters is selected from within plausible ranges of parameters by adjusting the values of hydraulic conductivity, storativity, groundwater recharge and stream bed permeability. The baseflow is determined from stream flow hydrographs by using an empirical scheme validated using a chemical approach to hydrograph separation. Field studies determined that the specific conductance for components of the composite hydrograph were sufficiently unique to make the chemical approach feasible. The method was applied to the Big Darby Creek Watershed, Ohio. The parameter set selected for the groundwater system provides a good fit with the estimated baseflow and observed water well data. Copyright © 1999 John Wiley & Sons, Ltd.     

Parameter estimation and uncertainty analysis of SWAT model in upper reaches of the Heihe river basin

Heihe river basin, the second largest inland river basin in China, has attracted more attention in China due to the ever increasing water resources and eco‐environmental problems. In this article, SWAT (Soil and Water Assessment Tool; http://www.brc.tamus.edu/swat/ ) model was applied to upper reaches of the basin for better understanding of the hydrological process over the watershed. Parameter uncertainty and its contribution on model simulation are the main foci. In model calibration, the aggregate parameters instead of the original parameters in SWAT model were used to reduce the computing effort. The Bayesian approach was employed for parameter estimation and uncertainty analysis because its posterior distribution provides not only parameter estimation but also uncertainty analysis without normality assumption. The results indicated that: (1) SWAT model performs satisfactorily in this watershed as a whole, although some low and high flows were under‐ or overestimated, particularly in dry (e.g. 1991) and wet (e.g. 1996) years; (2) all calibrated parameters were not normally distributed (essentially positively or negatively skewed) and the parameter uncertainties were relatively small; and (3) the contributions of parameter uncertainty on model simulation uncertainty were relatively small. Copyright © 2009 John Wiley & Sons, Ltd.     

Climate impacts on river flow: projections for the Medway catchment,UK, with UKCP09 and CATCHMOD

The potential impact of climate change on areas of strategic importance for water resources remains a concern. Here, river flow projections for the River Medway, above Teston in southeast England are presented, which is just such an area of strategic importance. The river flow projections use climate inputs from the Hadley Centre Regional Climate Model (HadRM3) for the time period 1960–2080 (a subset of the early release UKCP09 projections). River flow predictions are calculated using CATCHMOD, the main river flow prediction tool of the Environment Agency (EA) of England and Wales. In order to use this tool in the best way for climate change predictions, model setup and performance are analysed using sensitivity and uncertainty analysis. The model's representation of hydrological processes is discussed and the direct percolation and first linear storage constant parameters are found to strongly affect model results in a complex way, with the former more important for low flows and the latter for high flows. The uncertainty in predictions resulting from the hydrological model parameters is demonstrated and the projections of river flow under future climate are analysed. A clear climate change impact signal is evident in the results with a persistent lowering of mean daily river flows for all months and for all projection time slices. Results indicate that a projection of lower flows under future climate is valid even taking into account the uncertainties considered in this modelling chain exercise. The model parameter uncertainty becomes more significant under future climate as the river flows become lower. This has significant implications for those making policy decisions based on such modelling results. Copyright © 2010 John Wiley & Sons, Ltd.     

Towards hydrological model calibration and validation: simulation of stable water isotopes using the isoWATFLOOD model

Calibration and validation of hydrological models is a challenge, particularly in remote regions that are minimally gauged. This paper develops a novel methodology for large‐scale (>1000 km²) hydrological model calibration and validation using stable water isotopes founded on the rigorous constraints imposed by the need to conserve both water mass and stable isotopes simultaneously. The isoWATFLOOD model is applied to five basins within the Fort Simpson, Northwest Territories region of northern Canada to simulate stream discharge and oxygen‐18 signals over a 3‐year period. The isotopic variation of river discharge, runoff components, and evaporative fractionation are successfully simulated on both a seasonal and continual basis over the watershed domain to demonstrate the application of isotope tracers to regional hydrologic calibration. The intended application of this research is to remote, large‐scale basins, showing promise for improving predictions in minimally gauged basins and climate change research where traditional, rigorous approaches to constraining parameter uncertainty may be impractical. This coupled isotope‐hydrological (i.e. iso‐hydrological) approach to modelling reduces the number of possible parameterizations, resulting in potentially more physically‐based hydrological predictions. isoWATFLOOD provides a tool for water resource managers and utilities to use operationally for water use, allocation, and runoff generation estimations. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimation of uncertainty sources in the projections of Lithuanian river runoff

Particular attention is given to the reliability of hydrological modelling results. The accuracy of river runoff projection depends on the selected set of hydrological model parameters, emission scenario and global climate model. The aim of this article is to estimate the uncertainty of hydrological model parameters, to perform sensitivity analysis of the runoff projections, as well as the contribution analysis of uncertainty sources (model parameters, emission scenarios and global climate models) in forecasting Lithuanian river runoff. The impact of model parameters on the runoff modelling results was estimated using a sensitivity analysis for the selected hydrological periods (spring flood, winter and autumn flash floods, and low water). During spring flood the results of runoff modelling depended on the calibration parameters that describe snowmelt and soil moisture storage, while during the low water period—the parameter that determines river underground feeding was the most important. The estimation of climate change impact on hydrological processes in the Merkys and Neris river basins was accomplished through the combination of results from A1B, A2 and B1 emission scenarios and global climate models (ECHAM5 and HadCM3). The runoff projections of the thirty-year periods (2011–2040, 2041–2070, 2071–2100) were conducted applying the HBV software. The uncertainties introduced by hydrological model parameters, emission scenarios and global climate models were presented according to the magnitude of the expected changes in Lithuanian rivers runoff. The emission scenarios had much greater influence on the runoff projection than the global climate models. The hydrological model parameters had less impact on the reliability of the modelling results.     

Sensitivity analysis and uncertainty assessment in water budgets simulated by the variable infiltration capacity model for Canadian subarctic watersheds

In this study, we evaluate uncertainties propagated through different climate data sets in seasonal and annual hydrological simulations over 10 subarctic watersheds of northern Manitoba, Canada, using the variable infiltration capacity (VIC) model. Further, we perform a comprehensive sensitivity and uncertainty analysis of the VIC model using a robust and state-of-the-art approach. The VIC model simulations utilize the recently developed variogram analysis of response surfaces (VARS) technique that requires in this application more than 6,000 model simulations for a 30-year (1981–2010) study period. The method seeks parameter sensitivity, identifies influential parameters, and showcases streamflow sensitivity to parameter uncertainty at seasonal and annual timescales. Results suggest that the Ensemble VIC simulations match observed streamflow closest, whereas global reanalysis products yield high flows (0.5–3.0 mm day⁻¹) against observations and an overestimation (10–60%) in seasonal and annual water balance terms. VIC parameters exhibit seasonal importance in VARS, and the choice of input data and performance metrics substantially affect sensitivity analysis. Uncertainty propagation due to input forcing selection in each water balance term (i.e., total runoff, soil moisture, and evapotranspiration) is examined separately to show both time and space dimensionality in available forcing data at seasonal and annual timescales. Reliable input forcing, the most influential model parameters, and the uncertainty envelope in streamflow prediction are presented for the VIC model. These results, along with some specific recommendations, are expected to assist the broader VIC modelling community and other users of VARS and land surface schemes, to enhance their modelling applications.     

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.     

Understanding the time‐varying importance of different uncertainty sources in hydrological modelling using global sensitivity analysis

Simulations from hydrological models are affected by potentially large uncertainties stemming from various sources, including model parameters and observational uncertainty in the input/output data. Understanding the relative importance of such sources of uncertainty is essential to support model calibration, validation and diagnostic evaluation and to prioritize efforts for uncertainty reduction. It can also support the identification of ‘disinformative data’ whose values are the consequence of measurement errors or inadequate observations. Sensitivity analysis (SA) provides the theoretical framework and the numerical tools to quantify the relative contribution of different sources of uncertainty to the variability of the model outputs. In traditional applications of global SA (GSA), model outputs are aggregations of the full set of a simulated variable. For example, many GSA applications use a performance metric (e.g. the root mean squared error) as model output that aggregates the distances of a simulated time series to available observations. This aggregation of propagated uncertainties prior to GSA may lead to a significant loss of information and may cover up local behaviour that could be of great interest. Time‐varying sensitivity analysis (TVSA), where the aggregation and SA are repeated at different time steps, is a viable option to reduce this loss of information. In this work, we use TVSA to address two questions: (1) Can we distinguish between the relative importance of parameter uncertainty versus data uncertainty in time? (2) Do these influences change in catchments with different characteristics? To our knowledge, the results present one of the first quantitative investigations on the relative importance of parameter and data uncertainty across time. We find that the approach is capable of separating influential periods across data and parameter uncertainties, while also highlighting significant differences between the catchments analysed. Copyright © 2016 The Authors. Hydrological Processes. Published by John Wiley & Sons Ltd.     

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.     

Parameter Uncertainty Propagation in a Rainfall–Runoff Model; Case Study: Karoon-III River Basin

Conceptual hydrological models are popular tools for simulating land phase of hydrological cycle. Uncertainty arises from a variety of sources such as input error, calibration and parameters. Hydrologic modeling researches indicate that parametric uncertainty has been considered as one of the most important source. The objective of this study was to evaluate parameter uncertainty and its propagation in rainfall-runoff modeling. This study tried to model daily flows and calculate uncertainty bounds for Karoon-III basin, Southwest of Iran, using HEC-HMS (SMA). The parameters were represented by probability distribution functions (PDF), and the effect on simulated runoff was investigated using Latin Hypercube Sampling (LHS) on Monte Carlo (MC). Three chosen parameters, based on sensitivity analysis, were saturated-hydraulic-conductivity (Ks), Clark storage coefficient (R) and time of concentration (t _c ). Uncertainty associated with parameters were accounted for, by representing each with a probability distribution. Uncertainty bounds was calculated, using parameter sets captured from LHS on parameters PDF of sub-basins and propagating to the model. Results showed that maximum reliability (11%) resulted from Ks propagating. For three parameters, underestimation was more than overestimation. Maximum sharpness and standard deviation (STD) was resulted from propagating Ks. Cumulative Distribution Function (CDF) of flow and uncertainty bounds showed that as flow increased, the width of uncertainty bounds increased for all parameters.     

Observational and predictive uncertainties for multiple variables in a spatially distributed hydrological model

In this study, uncertainty in model input data (precipitation) and parameters is propagated through a physically based, spatially distributed hydrological model based on the MIKE SHE code. Precipitation uncertainty is accounted for using an ensemble of daily rainfall fields that incorporate four different sources of uncertainty, whereas parameter uncertainty is considered using Latin hypercube sampling. Model predictive uncertainty is assessed for multiple simulated hydrological variables (discharge, groundwater head, evapotranspiration, and soil moisture). Utilizing an extensive set of observational data, effective observational uncertainties for each hydrological variable are assessed. Considering not only model predictive uncertainty but also effective observational uncertainty leads to a notable increase in the number of instances, for which model simulation and observations are in good agreement (e.g., 47% vs. 91% for discharge and 0% vs. 98% for soil moisture). Effective observational uncertainty is in several cases larger than model predictive uncertainty. We conclude that the use of precipitation uncertainty with a realistic spatio‐temporal correlation structure, analyses of multiple variables with different spatial support, and the consideration of observational uncertainty are crucial for adequately evaluating the performance of physically based, spatially distributed hydrological models.     

Bayesian uncertainty analysis in hydrological modeling associated with watershed subdivision level: a case study of SLURP model applied to the Xiangxi River watershed,China

Uncertainty analysis in hydrological modeling would help to better implement decision-making related to water resources management, which relies heavily on hydrologic simulations. However, an important concern will be raised over the uncertainty associated with watershed subdivision broadly applied in distributed/semi-distributed hydrological models since scale issues would significantly affect model performance, and thus, lead to dramatic variations in simulations. To fully understand the uncertainty associated with watershed subdivision level, however, is still a tough work confronting researchers because of complex modeling processes and high computation requirements. In this study, we analyzed this uncertainty within a formal Bayesian framework using a Markov Chain Monte Carlo method based on Metropolis–Hastings algorithm. In a case study using the semi-distributed land use-based runoff processes hydrologic model in the Xiangxi River watershed, results showed that the variation in the simulated discharges due to parameter uncertainty was much smaller than that due to parameter and model uncertainty under different watershed subdivision levels defined using aggregated simulation areas (ASAs). However, the posterior probability distribution of model parameters varied in response to subdivision levels, and four parameters (i.e. maximum infiltration rate, retention constant for slow store, maximum capacity for slow store, and retention constant for fast store) were identified with smaller uncertainty. Although the uncertainty in the simulated discharge due to parameter and model uncertainty varied little across subdivisions, the simulation uncertainty only due to parameter uncertainty was found to be reduced through increasing the subdivisions. In addition, the coarsest subdivision level (7 ASAs) was not sufficient for obtaining satisfying simulations in the Xiangxi River watershed, but inappreciable improvement was achieved through increasing the level among finer subdivisions. Moreover, it was demonstrated that increasing subdivision level would have no advantage of improving the reliability of hydrological simulations beyond the threshold (45 ASAs). The findings of this research may shed light on the design of operational hydrological forecasting in the Three Gorges Reservoir region with profound socio-economic implications.     

Uncertainty assessment of a process-based integrated catchment model of phosphorus

Despite the many models developed for phosphorus concentration prediction at differing spatial and temporal scales, there has been little effort to quantify uncertainty in their predictions. Model prediction uncertainty quantification is desirable, for informed decision-making in river-systems management. An uncertainty analysis of the process-based model, integrated catchment model of phosphorus (INCA-P), within the generalised likelihood uncertainty estimation (GLUE) framework is presented. The framework is applied to the Lugg catchment (1,077 km²), a River Wye tributary, on the England–Wales border. Daily discharge and monthly phosphorus (total reactive and total), for a limited number of reaches, are used to initially assess uncertainty and sensitivity of 44 model parameters, identified as being most important for discharge and phosphorus predictions. This study demonstrates that parameter homogeneity assumptions (spatial heterogeneity is treated as land use type fractional areas) can achieve higher model fits, than a previous expertly calibrated parameter set. The model is capable of reproducing the hydrology, but a threshold Nash-Sutcliffe co-efficient of determination (E or R ²) of 0.3 is not achieved when simulating observed total phosphorus (TP) data in the upland reaches or total reactive phosphorus (TRP) in any reach. Despite this, the model reproduces the general dynamics of TP and TRP, in point source dominated lower reaches. This paper discusses why this application of INCA-P fails to find any parameter sets, which simultaneously describe all observed data acceptably. The discussion focuses on uncertainty of readily available input data, and whether such process-based models should be used when there isn’t sufficient data to support the many parameters.     

Modelling spatial variations in dissolved oxygen in fine‐grained streams under uncertainty

This paper presents a novel triple‐layer model, called VART DO‐3L, for simulation of spatial variations in dissolved oxygen (DO) in fine‐grained streams, characterized by a fluid mud (fluff or flocculent) layer (an advection‐dominated storage zone) as the interface between overlying stream water and relatively consolidated streambed sediment (a diffusion‐dominated storage zone). A global sensitivity analysis is conducted to investigate the sensitivity of VART DO‐3L model input parameters. Results of the sensitivity analysis indicate that the most sensitive parameter is the relative size of the advection‐dominated storage zones (A_s/A), followed by a lumped reaction term (R) for the flocculent layer, biological reaction rate (μ_o) in diffusive layer and biochemical oxygen demand concentration (L) in water column. In order to address uncertainty in model input parameters, Monte Carlo simulations are performed to sample parameter values and to produce various parameter combinations or cases. The VART DO‐3L model is applied to the Lower Amite River in Louisiana, USA, to simulate vertical and longitudinal variations in DO under the cases. In terms of longitudinal variation, the DO level decreases from 7.9 mg l at the Denham Springs station to about 2.89 mg l^?1 at the Port Vincent station. In terms of vertical variation, the DO level drops rapidly from the overlying water column to the advection‐dominated storage zone and further to the diffusive layer. The DO level (C_F) in the advective layer (flocculent layer) can reach as high as 40% of DO concentration (C) in the water column. The VART DO‐3L model may be applied to similar rivers for simulation of spatial variations in DO level. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling the hydrologic effects of dynamic land‐use change using a distributed hydrologic model and a spatial land‐use allocation model

This study develops a novel approach for modelling and examining the impacts of time–space land‐use changes on hydrological components. The approach uses an empirical land‐use change allocation model (CLUE‐s) and a distributed hydrological model (DHSVM) to examine various land‐use change scenarios in the Wu‐Tu watershed in northern Taiwan. The study also uses a generalized likelihood uncertainty estimation approach to quantify the parameter uncertainty of the distributed hydrological model. The results indicate that various land‐use policies—such as no change, dynamic change and simultaneous change—have different levels of impact on simulating the spatial distributions of hydrological components in the watershed study. Peak flow rates under simultaneous and dynamic land‐use changes are 5·71% and 2·77%, respectively, greater than the rate under the no land‐use change scenario. Using dynamic land‐use changes to assess the effect of land‐use changes on hydrological components is more practical and feasible than using simultaneous land‐use change and no land‐use change scenarios. Furthermore, land‐use change is a spatial dynamic process that can lead to significant changes in the distributions of ground water and soil moisture. The spatial distributions of land‐use changes influence hydrological processes, such as the ground water level of whole areas, particularly in the downstream watershed. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial uncertainty in bias corrected climate change projections and hydrogeological impacts

The question of which climate model bias correction methods and spatial scales for correction are optimal for both projecting future hydrological changes as well as removing initial model bias has so far received little attention. For 11 climate models (CMs), or GCM/RCM – Global/Regional Climate Model pairing, this paper analyses the relationship between complexity and robustness of three distribution‐based scaling (DBS) bias correction methods applied to daily precipitation at various spatial scales. Hydrological simulations are forced by CM inputs to assess the spatial uncertainty of groundwater head and stream discharge given the various DBS methods. A unique metric is devised, which allows for comparison of spatial variability in climate model bias and projected change in precipitation. It is found that the spatial variability in climate model bias is larger than in the climate change signals. The magnitude of spatial bias seen in precipitation inputs does not necessarily correspond to the magnitude of biases seen in hydrological outputs. Variables that integrate basin responses over time and space are more sensitive to mean spatial biases and less so on extremes. Hydrological simulations forced by the least parameterized DBS approach show the highest error in mean and maximum groundwater heads; however, the most highly parameterised DBS approach shows less robustness in future periods compared with the reference period it was trained in. For hydrological impacts studies, choice of bias correction method should depend on the spatial scale at which hydrological impacts variables are required and whether CM initial bias is spatially uniform or spatially varying. Copyright © 2015 John Wiley & Sons, Ltd.