Assessing uncertainties in a conceptual water balance model using Bayesian methodology / Estimation bayésienne des incertitudes au sein d’une modélisation conceptuelle de bilan hydrologique

Abstract

Abstract The aim of this study was to estimate the uncertainties in the streamflow simulated by a rainfall–runoff model. Two sources of uncertainties in hydrological modelling were considered: the uncertainties in model parameters and those in model structure. The uncertainties were calculated by Bayesian statistics, and the Metropolis-Hastings algorithm was used to simulate the posterior parameter distribution. The parameter uncertainty calculated by the Metropolis-Hastings algorithm was compared to maximum likelihood estimates which assume that both the parameters and model residuals are normally distributed. The study was performed using the model WASMOD on 25 basins in central Sweden. Confidence intervals in the simulated discharge due to the parameter uncertainty and the total uncertainty were calculated. The results indicate that (a) the Metropolis-Hastings algorithm and the maximum likelihood method give almost identical estimates concerning the parameter uncertainty, and (b) the uncertainties in the simulated streamflow due to the parameter uncertainty are less important than uncertainties originating from other sources for this simple model with fewer parameters.     

The comparative study of multi‐site uncertainty evaluation method based on SWAT model

With the recent development of distributed hydrological models, the use of multi‐site observed data to evaluate model performance is becoming more common. Distributed hydrological model have many advantages, and at the same time, it also faces the challenge to calibrate over‐do parameters. As a typical distributed hydrological model, problems also exist in Soil and Water Assessment Tool (SWAT) parameter calibration. In the paper, four different uncertainty approaches – Particle Swarm Optimization (PSO) techniques, Generalized Likelihood Uncertainty Estimation (GLUE), Sequential Uncertainty Fitting algorithm (SUFI‐2) and Parameter Solution (PARASOL) – are taken to a comparative study with the SWAT model applied in Peace River Basin, central Florida. In our study, the observed river discharge data used in SWAT model calibration were collected from the three gauging stations at the main tributary of the Peace River. Behind these approaches, there is a shared philosophy; all methods seek out many parameter set to fit the uncertainties due to the non‐uniqueness in model parameter evaluation. On the basis of the statistical results of four uncertainty methods, difficulty level of each method, the number of runs and theoretical basis, the reasons that affected the accuracy of simulation were analysed and compared. Furthermore, for the four uncertainty method with SWAT model in the study area, the pairwise correlation between parameters and the distributions of model fit summary statistics computed from the sampling over the behavioural parameter and the entire model calibration parameter feasible spaces were identified and examined. It provided additional insight into the relative identifiability of the four uncertainty methods Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Non‐stationary hydrological model parameters: a framework based on SOM‐B

The application of stationary parameters in conceptual hydrological models, even under changing boundary conditions, is a common yet unproven practice. This study investigates the impact of non‐stationary model parameters on model performance for different flow indices and time scales. Therefore, a Self‐Organizing Map based optimization approach, which links non‐stationary model parameters with climate indices, is presented and tested on seven meso‐scale catchments in northern Germany. The algorithm automatically groups sub‐periods with similar climate characteristics and allocates them to similar model parameter sets. The climate indices used for the classification of sub‐periods are based on (a) yearly means and (b) a moving average over the previous 61 days. Classification b supports the estimation of continuous non‐stationary parameters. The results show that (i) non‐stationary model parameters can improve the performance of hydrological models with an acceptable growth in parameter uncertainty; (ii) some model parameters are highly correlated to some climate indices; (iii) the model performance improves more for monthly means than yearly means; and (iv) in general low to medium flows improve more than high flows. It was further shown how the gained knowledge can be used to identify insufficiencies in the model structure. Copyright © 2015 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.     

Diagnostic evaluation of multiple hypotheses of hydrological behaviour in a limits‐of‐acceptability framework for 24 UK catchments

Testing competing conceptual model hypotheses in hydrology is complicated by uncertainties from a wide range of sources, which result in multiple simulations that explain catchment behaviour. In this study, the limits of acceptability uncertainty analysis approach used to discriminate between 78 competing hypotheses in the Framework for Understanding Structural Errors for 24 catchments in the UK. During model evaluation, we test the model's ability to represent observed catchment dynamics and processes by defining key hydrologic signatures and time step‐based metrics from the observed discharge time series. We explicitly account for uncertainty in the evaluation data by constructing uncertainty bounds from errors in the stage‐discharge rating curve relationship. Our study revealed large differences in model performance both between catchments and depending on the type of diagnostic used to constrain the simulations. Model performance varied with catchment characteristics and was best in wet catchments with a simple rainfall‐runoff relationship. The analysis showed that the value of different diagnostics in constraining catchment response and discriminating between competing conceptual hypotheses varies according to catchment characteristics. The information content held within water balance signatures was found to better capture catchment dynamics in chalk catchments, where catchment behaviour is predominantly controlled by seasonal and annual changes in rainfall, whereas the information content in the flow‐duration curve and time‐step performance metrics was able to better capture the dynamics of rainfall‐driven catchments. We also investigate the effect of model structure on model performance and demonstrate its (in)significance in reproducing catchment dynamics for different catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessment of the SWMM model uncertainties within the generalized likelihood uncertainty estimation (GLUE) framework for a high‐resolution urban sewershed

This research incorporates the generalized likelihood uncertainty estimation (GLUE) methodology in a high‐resolution Environmental Protection Agency Storm Water Management Model (SWMM), which we developed for a highly urbanized sewershed in Syracuse, NY, to assess SWMM modelling uncertainties and estimate parameters. We addressed two issues that have long been suggested having a great impact on the GLUE uncertainty estimation: the observations used to construct the likelihood measure and the sampling approach to obtain the posterior samples of the input parameters and prediction bounds of the model output. First, on the basis of the Bayes' theorem, we compared the prediction bounds generated from the same Gaussian distribution likelihood measure conditioned on flow observations of varying magnitude. Second, we employed two sampling techniques, the sampling importance resampling (SIR) and the threshold sampling methods, to generate posterior parameter distributions and prediction bounds, based on which the sampling efficiency was compared. In addition, for a better understanding of the hydrological responses of different pervious land covers in urban areas, we developed new parameter sets in SWMM representing the hydrological properties of trees and lawns, which were estimated through the GLUE procedure. The results showed that SIR was a more effective alternative to the conventional threshold sampling method. The combined total flow and peak flow data were an efficient alternative to the intensive 5‐min flow data for reducing SWMM parameter and output uncertainties. Several runoff control parameters were found to have a great effect on peak flows, including the newly introduced parameters for trees. Copyright © 2013 John Wiley & Sons, Ltd.     

Uncertainty quantization of meteorological input and model parameters for hydrological modelling using a Bayesian-based integrated approach

The traditional treatment of uncertainty in hydrological modelling primarily attributes it to model parameters, but rarely systematically considers meteorological input errors, especially in quantifying the impact of meteorological input errors on parameter uncertainty. This study developed a Bayesian-based integrated approach to quantitatively investigate uncertainties in meteorological inputs (precipitation and temperature) and model parameters as well as the variation in parameter uncertainty due to meteorological input errors. Additionally, we analysed the propagation from these uncertainties to runoff response in snowmelt and non-snowmelt periods. The applicability and advantages of this approach were presented by applying of the Soil and Water Assessment Tool to the Shitoukoumen Reservoir Catchment. Differential Evolution Adaptive Metropolis-Markov Chain Monte Carlo was applied for the straightforward Bayesian inference the uncertainties of meteorological inputs and model parameters. On this basis, multilevel factorial analysis technology was used to quantitatively investigate the specific impact of the model parameters' individual and interactive effects due to meteorological input errors. Finally, the impact of meteorological input errors and model parameter uncertainty on the model performance were analysed and quantified systematically. The results showed that the meteorological input errors could affect the random characteristics of multiple model parameters. Moreover, meteorological input errors could further affect the model parameters' effects on annual average runoff. Overall, the above results have significant implications in enhancing hydrological model to simulate/predict runoff and understanding hydrological processes during different periods.     

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.     

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.     

Predictive uncertainty of a nearshore bed evolution model

Predictions of nearshore depth evolution using process-based numerical simulation models contain inherent uncertainties owing to model structural deficiencies, measurement errors, and parameter uncertainty. This paper quantifies the parameter-induced predictive uncertainty of the cross-shore depth evolution model Unibest-TC by applying the Bayesian Generalised Likelihood Uncertainty Estimation methodology to modelling depth evolution at Egmond aan Zee (Netherlands). This methodology works with multiple sets of parameter values sampled uniformly in feasible parameter space and assigns a likelihood value to each parameter set. Acceptable simulations (i.e., based on parameter sets with a nonzero likelihood) were found for a wide range of parameter values owing to parameter interdependence and insensitivity. The 95% uncertainty prediction interval of bed levels after the 33 days prediction period was largest (0.5–1 m) near the sandbar crests that characterize the Egmond depth profile, reducing to near-zero values in the sandbar troughs and the offshore area. The prediction interval built up during storms (when sediment transport rates are largest) and remained the same or even reduced slightly during less-energetic conditions. The prediction uncertainty ranges bracket the observations near the inner-bar crest, its seaward flank, and at the seaward flank of the outer bar, suggesting that elsewhere model structural errors (and, potentially, measurement errors) dominate over parameter errors. The interdependence and the non-Gaussian marginal posterior distribution functions of the free model parameters cast doubt on the ability of commonly applied multivariate normal distribution functions to estimate parameter uncertainty.     

Estimation of catchment yield and associated uncertainties due to climate change in a mountainous catchment in Australia

This paper examines the impacts of climate change on future water yield with associated uncertainties in a mountainous catchment in Australia using a multi‐model approach based on four global climate models (GCMs), 200 realisations (50 realisations from each GCM) of downscaled rainfalls, 2 hydrological models and 6 sets of model parameters. The ensemble projections by the GCMs showed that the mean annual rainfall is likely to reduce in the future decades by 2–5% in comparison with the current climate (1987–2012). The results of ensemble runoff projections indicated that the mean annual runoff would reduce in future decades by 35%. However, considerable uncertainty in the runoff estimates was found as the ensemble results project changes of the 5^th (dry scenario) and 95^th (wet scenario) percentiles by ?73% to +27%, ?73% to +12%, ?77% to +21% and ?80% to +24% in the decades of 2021–2030, 2031–2040, 2061–2070 and 2071–2080, respectively. Results of uncertainty estimation demonstrated that the choice of GCMs dominates overall uncertainty. Realisation uncertainty (arising from repetitive simulations for a given time step during downscaling of the GCM data to catchment scale) of the downscaled rainfall data was also found to be remarkably high. Uncertainty linked to the choice of hydrological models was found to be quite small in comparison with the GCM and realisation uncertainty. The hydrological model parameter uncertainty was found to be lowest among the sources of uncertainties considered in this study. Copyright © 2015 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.     

On the significance of crop‐type information for the simulation of catchment hydrology

Spatially distributed hydrological models require information on the land cover pattern, as it directly influences the generation of run‐off. However, it is not clear which detail of land cover information is suitable for simulating the catchment hydrology. A better understanding of the relationship between the land cover detail and the hydrological processes is therefore required as it would enhance a successful application of the hydrological model. This study investigates the relevance of accurate information about the crop types in the catchment for the simulation of run‐off, baseflow and evapotranspiration. Results reveal that adding knowledge about the crop type in the model simulation is redundant when area‐averaged water flux predictions are intended. On the contrary, when a spatial distribution of water flux predictions is desired, it is meaningful to increase the land cover detail because an increased heterogeneity in the land cover map induces an increased heterogeneity in the hydrological response and locally reduces the uncertainty in the model prediction up to 50%. To assess the land cover uncertainty and to locate the regions for which the reduction in prediction uncertainty is the highest, the thematic uncertainty map (constructed through fuzzy logic) is shown to be a useful tool. Copyright © 2014 John Wiley & Sons, Ltd.     

Post‐processing of medium‐range probabilistic hydrological forecasting: impact of forcing,initial conditions and model errors

The impact of errors in the forcing, errors in the model structure and parameters, and errors in the initial conditions is investigated in a simple hydrological ensemble prediction system. The hydrological model is based on an input nonlinearity connected with a linear transfer function and forced by precipitation forecasts from the European Centre for Medium‐Range Weather Forecast (ECMWF) Ensemble Prediction System (EPS). The post‐processing of the precipitation and/or the streamflow using information from the reforecasts performed by ECMWF is tested. For this purpose, hydrological reforecasts are obtained by forcing the hydrological model with the precipitation from the reforecast data. In the present case study, it is found that the post‐processing of the hydrological ensembles with a statistical model fitted on the hydrological reforecasts improves the verification scores better than the use of post‐processed precipitation ensembles. In the case of large biases in the precipitation, combining the post‐processing of both precipitation and streamflow allows for further improvements. During the winter, errors in the initial conditions have a larger impact on the scores than errors in the model structure as designed in the experiments. Errors in the parameter values are largely corrected with the post‐processing. Copyright © 2014 John Wiley & Sons, Ltd.     

Discrete principal‐monotonicity inference for hydro‐system analysis under irregular nonlinearities,data uncertainties,and multivariate dependencies. Part I: methodology development

Hydrological system analyses are challenged by complexities of irregular nonlinearities, data uncertainties, and multivariate dependencies. Among them, the irregular nonlinearities mainly represent inexistence of regular functions for robustly simulating highly complicated relationships between variables. Few existing studies can enable reliable simulation of hydrological processes under these complexities. This may lead to decreased robustness of the constructed models, unfeasibility of suggestions for human activities, and damages to socio‐economy and eco‐environment. In the first of two companion papers, a discrete principal‐monotonicity inference (DPMI) method is proposed for hydrological systems analysis under these complexities. Normalization of non‐normally distributed samples and invertible restoration of modelling results are enabled through a discrete distribution transformation approach. To mitigate data uncertainties, statistical inference is employed to assess the significance of differences among samples. The irregular nonlinearity between the influencing factors (i.e. predictors) and the hydrological variable of interest (i.e. the predictand) is interpreted as piecewise monotonicity. Monotonicity is further represented as principal monotonicity under multivariate dependencies. Based on stepwise classification and cluster analyses, all paired samples representing the responsive relationship between the predictors and the predictand are discretized as a series of end nodes. A prediction approach is advanced for estimating the predictand value given any combination of predictors. The DPMI method can reveal evolvement rules of hydrological systems under these complexities. Reliance of existing hydro‐system analysis methods on predefined functional forms is removed, avoiding artificial disturbances, e.g. empiricism in selecting model functions under irregular nonlinearities, on the modelling process. Both local and global significances of predictors in driving the evolution of hydrological variables are identified. An analysis of interactions among these complexities is also achieved. The understanding obtained from the DPMI process and associated results can facilitate hydrological prediction, guide water resources management, improve hydro‐system analysis methods, or support hydrological systems analysis in other cases. The effectiveness and advantages of DPMI will be demonstrated through a case study of streamflow simulation in Xingshan Watershed, China, in another paper. Copyright © 2016 John Wiley & Sons, Ltd.     

Integrating glaciers raster‐based modelling in large catchments hydrological balance: the Rhone case study

A raster‐based glacier sub‐model was successfully introduced in the distributed hydrological model FEST‐WB to simulate the water balance and surface runoff of large Alpine catchments. The glacier model is based on temperature‐index approach for melt, on linear reservoir for melt water propagation into the ice and on mass balance for accumulation; the initialization of the volume of ice on the basin was based on a formulation depending on surface topography. The model was first tested on a sub‐basin of the Rhone basin (Switzerland), which is for 62% glaciated; the calibration and validation were based on comparison between simulated and observed discharge from 1999 to 2008. The model proved to be suitable to simulate the typical discharge seasonality of a heavily glaciated basin. The performance of the model was also tested by simulating discharge in the whole Swiss Rhone basin, in which glaciers contribution is not negligible, in fact, in summer, about the 40% of the discharge is due to glacier melt. The model allowed to take into account the volume of water coming from glaciers melt and its simple structure is suitable for analysis of the effects of climate change on hydrological regime of high mountain basins, with available meteorological forcing from current RCM. Copyright © 2012 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     

Multi-period calibration of a semi-distributed hydrological model based on hydroclimatic clustering

Changing climatic conditions contribute to a time varying nature of hydrological responses over different temporal scales. The temporal dynamics of hydrological systems bring uncertainties into hydrological simulation which are different to uncertainties from spatial heterogeneity of soil and land use. This study develops a new approach to improve the calibration of hydrological based on hydroclimatic similarities. Six climatic indexes are integrated using Principal Component Analysis and Fuzzy C-mean Clustering methods to transform hydrological years into hydroclimatic periods. Parameter sets of SWAT model are calibrated independently for each period and used together to generate continuous simulation for a prairie watershed in southern Canada. Results indicate that the multi-period model exhibits comprehensive advantages over the traditional single-period model under various flow conditions. The simulation ability of the model is improved through using period-specific parameter sets in fitting the observations to compensate for deficiencies in the model structure or input data.