Parameter and rating curve uncertainty propagation analysis of the SWAT model for two small Mediterranean catchments

Abstract

The SWAT model was tested to simulate the streamflow of two small Mediterranean catchments (the Vène and the Pallas) in southern France. Model calibration and prediction uncertainty were assessed simultaneously by using three different techniques (SUFI-2, GLUE and ParaSol). Initially, a sensitivity analysis was conducted using the LH-OAT method. Subsequent sensitive parameter calibration and SWAT prediction uncertainty were analysed by considering, firstly, deterministic discharge data (assuming no uncertainty in discharge data) and secondly, uncertainty in discharge data through the development of a methodology that accounts explicitly for error in the rating curve (the stage?discharge relationship). To efficiently compare the different uncertainty methods and the effect of the uncertainty of the rating curve on model prediction uncertainty, common criteria were set for the likelihood function, the threshold value and the number of simulations. The results show that model prediction uncertainty is not only case-study specific, but also depends on the selected uncertainty analysis technique. It was also found that the 95% model prediction uncertainty interval is wider and more successful at encompassing the observations when uncertainty in the discharge data is considered explicitly. The latter source of uncertainty adds additional uncertainty to the total model prediction uncertainty.

Editor D. Koutsoyiannis; Associate editor D. Gerten

Citation Sellami, H., La Jeunesse, I., Benabdallah, S., and Vanclooster, M., 2013. Parameter and rating curve uncertainty propagation analysis of the SWAT model for two small Mediterranean watersheds. Hydrological Sciences Journal, 58 (8), 1635?1657.     

A soft hydrological monitoring approach for comparing runoff on a network of small poorly gauged catchments

Catchments in many parts of the world are either ungauged or poorly gauged, and the dominant processes governing their streamflow response are still poorly understood. The analysis of runoff coefficients provides essential insight into catchment response, particularly if both range of catchments and a range of events are compared. This paper investigates how well the hydrological runoff of 11 small, poorly gauged catchments with ephemeral streams (0·1‐0·6 km²) can be compared using estimated runoff with the associated uncertainty. Data of rainfall and water depth at a catchment's outlet were recorded using automatic logging equipment during 2008‐2009. The hydrological regime is intermittent and the annual precipitation ranged between 569 and 727 mm. Discharge was estimated using Manning's equation and channel cross‐section measurements. Innovative work has been performed under controlled experimental conditions to estimate Manning's coefficient values for the different cover types observed in studied streams: non‐aquatic vegetations (giant reed, bramble and thistle), grass and coarse granular deposits. The results show that estimates derived using roughness coefficients differ from those previously established for larger streams with aquatic vegetation. Catchment runoff was compared at both the event and the annual scale. The results indicate significant variability between the catchment's responses. This variability allows for classification in spite of all the uncertainty associated with runoff estimation. This study highlights the potential of using a network of poorly gauged catch ments. From almost no catchment understanding the proposed methodology allows to compare poorly gauged catchments and highlights similarity/dissimilarity between catchment responses. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial scale effect on seasonal streamflows in permafrost catchments on the Qinghai–Tibet Plateau

The scale issue is of central concern in hydrological processes to understand the potential upscaling or downscaling methodologies, and to develop models for scaling the dominant processes at different scales and in different environments. In this study, a typical permafrost watershed in the Qinghai‐Tibet Plateau was selected. Its hydrological processes were monitored for 4 years from 2004 to 2008, measuring the effects of freezing and thawing depth of active soil layers on runoff processes. To identify the nature and cause of variation in the runoff response in different size catchments, catchments ranging from 1·07 to 112 km² were identified in the watershed. The results indicated that the variation of runoff coefficients showed a ‘V’ shape with increasing catchment size during the spring and autumn seasons, when the active soil was subjected to thawing or freezing processes. A two‐stage method was proposed to create runoff scaling models to indicate the effects of scale on runoff processes. In summer, the scaling transition model followed an exponential function for mean daily discharge, whereas the scaling model for flood flow exhibited a linear function. In autumn, the runoff process transition across multiple scales followed an exponential function with air temperature as the driving factor. These scaling models demonstrate relatively high simulation efficiency and precision, and provide a practical way for upscaling or downscaling runoff processes in a medium‐size permafrost watershed. For permafrost catchments of this scale, the results show that the synergistic effect of scale and vegetation cover is an important driving factor in the runoff response. Copyright © 2011 John Wiley & Sons, Ltd.     

Appropriate temporal resolution of precipitation data for discharge modelling in pre-alpine catchments

Precipitation time series with high temporal resolution are desired for hydrological modelling and flood studies. Yet the choice of an appropriate resolution is not straightforward because the use of too high a temporal resolution increases the data requirements, computational costs and, presumably, associated uncertainty, while performance improvement may be indiscernible. In this study, the effect of averaging hourly precipitation on model performance and associated uncertainty is investigated using two data sources: station network precipitation (SNP) and radar-based precipitation (RBP). From these datasets, time series of different temporal resolutions were generated, and runoff was simulated for 13 pre-alpine catchments with a bucket-type model. Our results revealed that different temporal resolutions were required for an acceptable model performance depending on the catchment size and data source. These were 1–12 h for small (16–59 km²), 3-21 h for medium (60–200 km²), and 24 h for large (200–939 km²) catchments.     

Linking physiography and evaporation using the isotopic composition of river water in 16 Canadian boreal catchments

Highly seasonal boreal catchments are hydrologically complex and generally data poor and, hence, are ripe for investigation using tracer‐aided hydrologic models. The influence of physiography on isotopic metrics was assessed to identify the catchment characteristics dominating evaporative enrichment. A multiyear stable isotope of water dataset was collected at the outlets of 16 boreal catchments in central Canada ranging in area from 12 to 15,282 km². Physiographic characteristics were obtained through raster analysis of freely available land cover images, stream networks, and digital elevation models. Correlation analysis indicated that as the percentage coverage of open water increased, so too did the evaporative effects observed at the catchment outlet. Correlation to wetland metrics indicated that increasing the percentage coverage of wetlands can reduce or increase evaporative effects observed, depending on the isotopic metric used and the corresponding drainage density, catchment slope, and presence of headwater lakes. The slopes of river evaporative‐mixing lines appear to reflect multifaceted relationships, strongest between catchment slope, headwater lakes, and connected wetlands, whereas mean line‐conditioned excess is more directly linked to physiographic variables. Hence, the slopes of river evaporative‐mixing lines and mean line‐conditioned excess are not interchangeable metrics of evaporative enrichment in a catchment. Relationships identified appear to be independent of catchment scale. These results suggest that adequate inclusion of the distribution of open water throughout a catchment, adequate representation of wetland processes, catchment slope, and drainage density are critical characteristics to include in tracer‐aided hydrologic models in boreal environments in order to minimize structural uncertainty.     

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.     

Hydrological modelling of a small watershed using generated rainfall in the soil and water assessment tool model

An adequately tested soil and water assessment tool (SWAT) model was applied to the runoff and sediment yield of a small agricultural watershed in eastern India using generated rainfall. The capability of the model for generating rainfall was evaluated for a period of 18 years (1981–1998). The watershed and subwatershed boundaries, drainage networks, slope, soil series and texture maps were generated using a geographical information system (GIS). A supervised classification method was used for land‐use/cover classification from satellite imageries. Model simulated monthly rainfall for the period of 18 years was compared with observations. Simulated monthly rainfall, runoff and sediment yield values for the monsoon season of 8 years (1991–1998) were also compared with their observed values. In general monthly average rainfall predicted by the model was in close agreement with the observed monthly average values. Also, simulated monthly average values of surface runoff and sediment yield using generated rainfall compared well with observed values during the monsoon season of the years 1991–1998. Results of this study revealed that the SWAT model can generate monthly average rainfall satisfactorily and thereby can produce monthly average values of surface runoff and sediment yield close to the observed values. Therefore, it can be concluded that the SWAT model could be used for developing a multiple year management plan for the critical erosion prone areas of a small watershed. Copyright © 2004 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.     

Hydrological regime of remote catchments with extreme gradients under accelerated change: the Baker basin in Patagonia

Abstract

The Baker basin (27 000 km²) is located in one of the most pristine and remote areas of the planet. Its hydrological regime is poised to undergo dramatic changes in the near future due to hydropower development and climate change. The basin contains the second-largest lake in South America, and part of a major icefield. This study documents the natural baseline of the Baker River basin, discusses the main hydrological modes and analyses the potential for sustainable management. Annual precipitation varies several-fold from the eastern Patagonian steppes to the North Patagonian Icefield. The westernmost sub-basins are strongly governed by glacier melt with a peak discharge in the austral summer (January–March). The easternmost sub-basins have a much more seasonal response governed by quicker snowmelt in spring (November–December), while they exhibit low flows typical for semi-arid regions during summer and autumn. Topography, vegetation and wetlands may also influence streamflow. The strong spatio-temporal gradients and variability highlight the need for further monitoring, particularly in the headwaters, especially given the severe changes these basins are expected to undergo. The great diversity of hydrological controls and climate change pose significant challenges for hydrological prediction and management.

Editor Z.W. Kundzewicz

Citation Dussaillant, J.A., Buytaert, W., Meier, C., and Espinoza, F. 2012. Hydrological regime of remote catchments with extreme gradients under accelerated change: the Baker basin in Patagonia. Hydrological Sciences Journal, 57 (8), 1530–1542.     

Application of connectivity theory to model the impact of terrace failure on runoff in semi‐arid catchments

Terraces are a common feature of Mediterranean landscapes. In many places they are no longer maintained so that the number of intact terraces is in prolonged decline. The aim of this paper is to examine the effect of terrace removal and failure on hydrological connectivity and peak discharge in an agricultural catchment (475 ha) in south‐east Spain. The situation of 2006 is compared to that in 1956 and to a scenario without terraces (S2). The spatial distribution of concentrated flow was mapped after four storms in 2006. The degree of connectivity was quantified by means of connectivity functions and related to storm characteristics, land use and topography. For 1956, 2006 and scenario S2, connectivity functions and peak discharge to the river were determined for a storm with a return period of 8·2 years. The results show that the decrease in intact terraces has led to a strong increase in connectivity and discharge. The contributing area to the river system has increased by a factor 3·2 between 1956 and 2006. If all terraces were to be removed (scenario S2), the contributing area may further increase by a factor 6·0 compared to 2006. The spatial extent of concentrated flow and the degree of connectivity are related to storm magnitude as expressed by the erosivity index (EI₃₀). Although a large part of the concentrated flow (25–50%) occurs on dirt roads, it appears that croplands become a major source of runoff with increasing rainfall. The results suggest that connectivity theory can be used to improve rainfall–runoff models in semi‐arid areas. Copyright © 2009 John Wiley & Sons, Ltd.     

太湖1889洪水年的流域水文模拟

由于受到水文观测资料时间短的限制,目前难以认识百年遇机率的极端洪水.为此,本文根据19世纪末历史文献的洪水灾害记录,利用流域水文模型,对太湖1889洪水年的流域产流、入湖汇流等水文特征和过程进行模拟.本研究设计了三套模拟实验:首先在现代气候控制实验基础上对1988-2002年时间系列和特大洪水年进行水文模拟和模型率定校验;其次,采用长江下游19世纪末的气候观测资料驱动,对极端年份1889年逐日洪水过程模拟;最后,为减少1 a洪水年模拟的不确定性,还采用蒙特卡罗Bootstrap法模拟了15 a的流域气候场,在5475 d样本下进行特征年份的水文模拟.模拟结果表明,1889年洪水期间产流在当年6月底达到最大,1%频率的径流深达8.6 mm/d,95%CI的误差在-2.94～3.26 mm/d之间.汇入太湖径流同期达到最大,1%频率的洪水流量达到1286.9 m3/s,95%CI的误差在-128.3～165.7 m3/s之间.根据洪水Log-Normal概率分布,计算1889洪水年的重现期为149 a.经Bootstrap法对误差置信区的模拟,95%CI检验在70～175 a间的重现期可信.该研究为延长20世纪洪水序列、拓展对百年时间尺度的特大洪水的认识提供了动力学模拟方面的科学依据.     

Runoff production processes in small alpine catchments within the unconsolidated Pleistocene sediments of the Lainbach area (Upper Bavaria)

In the Lainbach catchment, unconsolidated Pleistocene moraine sediments are widely distributed. Because of the great natural risk of floods, together with extreme loads of sediments, investigations of runoff production processes have been conducted in this area. At hillslope scale three test sites with different states of soil development and vegetation cover were instrumented with V‐shaped weirs, precipitation gauges and measurement devices for electrical conductivity (EC) of discharge water. The EC has been used as a geochemical tracer for hydrograph separation, since the statistical relationship between content of dissolved Ca²⁺, Mg²⁺ cations and EC is highly significant for different stages of runoff. This method allows hydrograph separation at high temporal resolution for both the rising and falling limb of the hydrograph. The following results of the investigations can be resumed. If relief conditions are similar, the effectiveness of runoff production decreases with an increasing density of vegetation cover. The runoff delivery ratio decreases as well as the peaks of runoff. In contrast, concentration times of hillslope catchments are equal, even if vegetation cover is of great density and soils are well developed. As a reason for the short reaction times, different runoff production processes have been detected. On bare ground, infiltration excess overland flow intensified by surface sealing processes is the main source for quick runoff. On hillslopes well covered by vegetation, translatory flow processes indicated by soil water with high solute contents force a rapid runoff reaction only a few minutes after rainfall has begun. It is to be assumed that translatory flow is a runoff production process typical for hillslopes covered by vegetation in a steep alpine relief. By means of the areal distribution of the topographic index, concentration of runoff production on a small part of the catchment has been demonstrated for hillslopes densely covered by vegetation. The investigations have shown that there is a lack of studies on runoff production processes in steep alpine relief, as well as a deficit of methods to quantify hydraulic properties of coarse‐grained soils with a wide grain size distribution. Copyright © 2003 John Wiley & Sons, Ltd.     

Runoff evolution in the Swiss Alps: projections for selected high‐alpine catchments based on ENSEMBLES scenarios

The Alps are often referred to as the ‘water tower of Europe’. In Switzerland, many branches of the economy, especially the hydropower industry, are closely linked to and dependent on the availability of water. Assessing the impact of climate change on streamflow runoff is, thus, of great interest. Major efforts have already been made in this respect, but the analyses often focus on individual catchments and are difficult to intercompare. In this article, we analysed nine high‐alpine catchments spread over the Swiss Alps, selected for their relevance to a wide range of morphological characteristics. Runoff projections were carried out until the end of the current century by applying the Glacier Evolution Runoff Model (GERM) and climate scenarios generated in the framework of the ENSEMBLES project. We focused on assessing the uncertainty induced by the unknown climate evolution and provided general, statistically based statements, which should be useful as a ‘rule of thumb’ for analyses addressing questions related to water management. Catchments with a high degree of glacierization will undergo the largest changes. General statements about absolute variations in discharge are unreliable, but an overall pattern, with an initial phase of increased annual discharge, followed by a phase with decreasing discharge, is recognizable for all catchments with a significant degree of glacierization. In these catchments, a transition from glacial and glacio‐nival regime types to nival will occur. The timing of maximal annual runoff is projected to occur before 2050 in all basins. The time of year with maximal daily discharges is expected to occur earlier at a rate of 4·4 ± 1·7 days per decade. Compared to its present level, the contribution of snow‐ and icemelt to annual discharge is projected to drop by 15 to 25% until the year 2100. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrological modelling of the Kangsabati river under changed climate scenario: case study in India

India is a large developing country with nearly two‐thirds of the population depending directly on the climate‐sensitive sectors such as agriculture, fisheries and forests. A very well‐calibrated Soil and Water Assessment Tool (R² = 0·9968, NSE = 0·91) was exercised over the Kangsabati river watershed in Bankura district of West Bengal, India, for a year including monsoon and non‐monsoon period in order to evaluate projected parameters for agricultural activities. Evapotranspiration, transmission losses, potential evapotranspiration and lateral flow to reach are evaluated from the years 2041–2050 in order to generate a picture for sustainable development of the river basin and its inhabitants. The projected climate change under various scenarios is likely to have implications on food production, water supply, biodiversity and livelihoods. India has a significant stake in scientific advancement as well as an international understanding to promote mitigation and adaptation. This requires improved scientific understanding, capacity building, networking and broad consultation processes. This paper is a commitment towards the planning, management and development of the water resources of the Kangsabati river by presenting detailed future scenarios of the Kangsabati river basin over the mentioned time period. The major findings of this paper were that of all the chosen projected parameters, transmission losses, soil water content, potential evapotranspiration, evapotranspiration and lateral flow to reach, display an increasing trend over the time period of years 2041–2050. Copyright © 2009 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.     

Testing the applicability of a kinematic wave-based distributed hydrological model in two climatically contrasting catchments

Abstract

Steep mountainous areas account for 70% of all river catchments in Japan. To predict river discharge for the mountainous catchments, many studies have applied distributed hydrological models based on a kinematic wave approximation with surface and subsurface flow components (DHM-KWSS). These models reproduce observed river discharge of catchments in Japan well; however, the applicability of a DHM-KWSS to catchments with different geographical and climatic conditions has not been sufficiently examined. This research applied a DHM-KWSS to two river basins that have different climatic conditions from basins in Japan to examine the transferability of the DHM-KWSS model structure. Our results show that the DHM-KWSS model structure explained flow regimes for a wet river basin as well as a large flood event in an arid basin; however, it was unable to explain long-term flow regimes for the arid basin case study.     

动力和统计预测信息融合预测方法及对我国夏季降水预测的检验

短期气候预测中如何将气候模式和统计方法的预测结果科学、客观的集成起来,一直是非常重要的问题.本文针对动力模式和统计方法预测结果相结合的问题,引入资料同化中信息融合的思想,采用最优内插同化方法,实现了动力模式和统计季节降水预测结果的融合.检验表明,对1982-2015年我国夏季降水百分率的回报,融合预测结果与观测的平均空间相关系数可达0.44,分别较统计预测和CFSv2模式统计降尺度订正的技巧提高了0.1左右,而均方根误差较两者可以降低5%~20%.可见,该方法可以进一步提升对我国夏季降水的预测技巧,具有显著的业务应用价值.

     

Assessing the effectiveness of imperviousness on stormwater runoff in micro urban catchments by model simulation

Imperviousness, considered as a critical indicator of the hydrologic impacts of urbanization, has gained increasing attention both in the research field and in practice. However, the effectiveness of imperviousness on rainfall–runoff dynamics has not been fully determined in a fine spatiotemporal scale. In this study, 69 drainage subareas <1 ha of a typical residential catchment in Beijing were selected to evaluate the hydrologic impacts of imperviousness, under a typical storm event with a 3‐year return period. Two metrics, total impervious area (TIA) and effective impervious area (EIA), were identified to represent the impervious characteristics of the selected subareas. Three runoff variables, total runoff depth (TR), peak runoff depth (PR), and lag time (LT), were simulated by using a validated hydrologic model. Regression analyses were developed to explore the quantitative associations between imperviousness and runoff variables. Then, three scenarios were established to test the applicability of the results in considering the different infiltration conditions. Our results showed that runoff variables are significantly related to imperviousness. However, the hydrologic performances of TIA and EIA were scale dependent. Specifically, with finer spatial scale and the condition heavy rainfall, TIA rather than EIA was found to contribute more to TR and PR. EIA tended to have a greater impact on LT and showed a negative relationship. Moreover, the relative significance of TIA and EIA was maintained under the different infiltration conditions. These findings may provide potential implications for landscape and drainage design in urban areas, which help to mitigate the runoff risk. 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.