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     

Contribution to hydrological modelling of small Amazonian catchments: application of rainfall–runoff models to simulate flow duration curves

Abstract

The objective of this study is to analyse three rainfall–runoff hydrological models applied in two small catchments in the Amazon region to simulate flow duration curves (FDCs). The simple linear model (SLM) considers the rainfall–runoff process as an input–output time-invariant system. However, the rainfall–runoff process is nonlinear; thus, a modification is applied to the SLM based on the residual relationship between the simulated and observed discharges, generating the modified linear model (MLM). In the third model (SVM), the nonlinearity due to infiltration and evapotranspiration is incorporated into the system through the sigmoid variable gain factor. The performance criteria adopted were a distance metric (δ) and the Nash-Sutcliffe coefficient (R²) determined between simulated and observed flows. The good results of the models, mainly the MLM and SVM, showed that they could be applied to simulate FDCs in small catchments in the Amazon region.

Editor D. Koutsoyiannis; Associate editor A. Montanari

Citation Blanco, C.J.C., Santos, S.S.M., Quintas, M.C., Vinagre, M.V.A., and Mesquita, A.L.A., 2013. Contribution to hydrological modelling of small Amazonian catchments: application of rainfall–runoff models to simulate flow duration curves. Hydrological Sciences Journal, 58 (7), 1–11.     

A spatially distributed Clark’s unit hydrograph based hybrid hydrologic model (Distributed-Clark)

ABSTRACT

A hybrid hydrologic model (Distributed-Clark), which is a lumped conceptual and distributed feature model, was developed based on the combined concept of Clark’s unit hydrograph and its spatial decomposition methods, incorporating refined spatially variable flow dynamics to implement hydrological simulation for spatially distributed rainfall–runoff flow. In Distributed-Clark, the Soil Conservation Service (SCS) curve number method is utilized to estimate spatially distributed runoff depth and a set of separated unit hydrographs is used for runoff routing to obtain a direct runoff flow hydrograph. Case studies (four watersheds in the central part of the USA) using spatially distributed (Thiessen polygon-based) rainfall data of storm events were used to evaluate the model performance. Results demonstrate relatively good fit to observed streamflow, with a Nash-Sutcliffe efficiency (E_NS) of 0.84 and coefficient of determination (R²) of 0.86, as well as a better fit in comparison with outputs of spatially averaged rainfall data simulations for two models including HEC-HMS.     

Rainfall–runoff simulation using a normalized antecedent precipitation index

Abstract

The normalized antecedent precipitation index (NAPI) model by Heggen for the prediction of runoff yield is analytically derived from the water balance equation. Heggen's model has been simplified further to a rational form and its performance verified with the Soil Conservation Service Curve Number (SCS-CN) model. The simplified model has three coefficients specific to a watershed, and requires two inputs: rainfall and the derived parameter, NAPI. The characteristic behaviour of the NAPI has resonance with the curve number (CN) of the SCS model. The proposed NAPI model was applied to three watersheds in the semi-arid region of India to simulate runoff yield. The model showed improved correlation between the observed and predicted runoff data compared to the SCS-CN model. The F test and paired t test also confirmed the reliability of the model with significance levels of 0.01 and 0.001%, respectively. The proposed model could be used successfully for rainfall–runoff modelling in a watershed.

Citation Ali, S., Ghosh, N. C. & Singh, R. (2010) Rainfall–runoff simulation using a normalized antecedent precipitation index. Hydrol. Sci. J. 55(2), 266–274.     

Study of discontinuities in hydrological data using catastrophe theory

Abstract

Modelling and prediction of hydrological processes (e.g. rainfall–runoff) can be influenced by discontinuities in observed data, and one particular case may arise when the time scale (i.e. resolution) is coarse (e.g. monthly). This study investigates the application of catastrophe theory to examine its suitability to identify possible discontinuities in the rainfall–runoff process. A stochastic cusp catastrophe model is used to study possible discontinuities in the monthly rainfall–runoff process at the Aji River basin in Azerbaijan, Iran. Monthly-averaged rainfall and flow data observed over a period of 20 years (1981–2000) are analysed using the Cuspfit program. In this model, rainfall serves as a control variable and runoff as a behavioural variable. The performance of this model is evaluated using four measures: correlation coefficient, log-likelihood, Akaike information criterion (AIC) and Bayesian information criterion (BIC). The results indicate the presence of discontinuities in the rainfall–runoff process, with a significant sudden jump in flow (cusp signal) when rainfall reaches a threshold value. The performance of the model is also found to be better than that of linear and logistic models. The present results, though preliminary, are promising in the sense that catastrophe theory can play a possible role in the study of hydrological systems and processes, especially when the data are noisy.

Citation Ghorbani, M. A., Khatibi, R., Sivakumar, B. & Cobb, L. (2010) Study of discontinuities in hydrological data using catastrophe theory. Hydrol. Sci. J. 55(7), 1137–1151.     

Linear and nonlinear input/output models for karstic springflow and flood prediction at different time scales

Karstic formations function as three-dimensional (3D) hydrological basins, with both surface and subsurface flows through fissures, natural conduits, underground streams and reservoirs. The main characteristic of karstic formations is their significant 3D physical heterogeneity at all scales, from fine fissuration to large holes and conduits. This leads to dynamic and temporal variability, e.g. highly variable flow rates, due to several concurrent flow regimes with several distinct response times. The temporal hydrologic response of karstic basins is studied here from an input/output, systems analysis viewpoint. The hydraulic behaviour of the basins is approached via the relationship between hydrometeorological inputs and outputs. These processes are represented and modeled as random, self-correlated and cross-correlated, stationary time processes. More precisely, for each site-specific case presented here, the input process is the total rainfall on the basin and the output process is the discharge rate at the outlet of the basin (karstic spring). In the absence of other data, these time processes embody all the available information concerning a given karstic basin. In this paper, we first present a brief discussion of the physical structure of karstic systems. Then, we formulate linear and nonlinear models, i.e. functional relations between rainfall and runoff, and methods for identifying the kernel and coefficients of the functionals (deterministic vs. statistical; error minimisation vs. polynomial projection). These are based mostly on Volterra first order (linear) or second order (nonlinear) convolution. In addition, a new nonlinear threshold model is developed, based on the frequency distribution of interannual mean daily runoff. Finally, the different models and identification methods are applied to two karstic watersheds in the french Pyrénées mountains, using long sequences of rainfall and spring outflow data at two different sampling rates (daily and semi-hourly). The accuracy of nonlinear and linear rainfall-runoff models is tested at three time scales: long interannual scale (20 years of daily data), medium or seasonal scale (3 months of semi-hourly data), and short scale or “flood scale” (2 days of semi-hourly data). The model predictions are analysed in terms of global statistical accuracy and in terms of accuracy with respect to high flow events (floods).     

Sediment yield variability in central Tunisia: a quantitative analysis of its controlling factors

Abstract

This paper refers to the quantification and prediction of the sedimentation rate of 26 hillside-dam reservoirs in Central Tunisia. The objectives of the study are to develop a simple and practical methodology to identify controlling factors of sedimentation, and to propose a regionalization from the study sites. Principal component analysis (PCA) and complementary multi-dimensional statistical methods are used to relate highly variable area-specific sediment yield to hydro-morphometric, lithological, geomorphological and anthropogenic characteristics of catchments. It appears that catchment area is not the main controlling factor of sedimentation in the studied area. The overall slope index, drainage network characteristics and runoff parameters are also important in characterizing sediment yield. Applied to the annual sedimentation rate series, PCA resulted in retaining the first three principal axes, explaining 65% of the total variance. Statistical methods showed that the overall slope index, the total drainage length, the compacity index and the runoff parameters are as important for the sedimentation quantification. This allowed a graphical clustering of the study zone into three distinct groups having similar behaviours: (i) watersheds characterized by high sediment transport rates and high runoff coefficients, (ii) basins distinguished by relatively low values of both flow discharge and sediment transport rates, and (iii) watersheds with an intermediate sediment yield, especially characterized by relatively high relief. In a second step, a multiple regression model including the four characteristic catchment properties was developed, presenting a valuable tool to predict area-specific sediment yield from catchments in central Tunisia. This model shows reasonable efficiency with an absolute prediction error of 81%.

Citation Ayadi, I., Abida, H., Djebbar, Y. & Mahjoub, M. R. (2010) Sediment yield variability in central Tunisia: a quantitative analysis of its controlling factors. Hydrol. Sci. J. 55(3), 446–458.     

Comparison of data-driven modelling techniques for river flow forecasting

Abstract

Accurate forecasting of streamflow is essential for the efficient operation of water resources systems. The streamflow process is complex and highly nonlinear. Therefore, researchers try to devise alterative techniques to forecast streamflow with relative ease and reasonable accuracy, although traditional deterministic and conceptual models are available. The present work uses three data-driven techniques, namely artificial neural networks (ANN), genetic programming (GP) and model trees (MT) to forecast river flow one day in advance at two stations in the Narmada catchment of India, and the results are compared. All the models performed reasonably well as far as accuracy of prediction is concerned. It was found that the ANN and MT techniques performed almost equally well, but GP performed better than both these techniques, although only marginally in terms of prediction accuracy in normal and extreme events.

Citation Londhe, S. & Charhate, S. (2010) Comparison of data-driven modelling techniques for river flow forecasting. Hydrol. Sci. J. 55(7), 1163–1174.     

A combined methodology for estimating the potential natural aquifer recharge in an arid environment

ABSTRACT

An innovative methodology that combines an indirect physiography-based method for determining the runoff coefficient at a sub-basin scale and a water balance model applied on a daily time scale was developed to calculate the natural groundwater recharge in three watersheds within the Oum Zessar arid area, Tunisia. The effective infiltration was calculated as part of the water surplus by considering the average available water content (AWC) of soil and an average runoff coefficient for each sub-basin. The model indicates that the sub-basins covered mainly by the “artificial” soils of tabias and jessour, characterized by average AWC values greater than 150 mm, did not contribute to natural groundwater recharge over the 10-year period (2003–2012) considered. The estimated volume for the Triassic aquifer amounted to about 4.5 hm³ year^?1, which is consistent with previous studies. For the Jurassic and Cretaceous aquifers, the estimated volumes amounted to about 200 dm³ year^?1.     

Prediction of streamflow regimes over large geographical areas: interpolated flow–duration curves for the Danube region

ABSTRACT

Flow–duration curves (FDCs) are essential to support decisions on water resources management, and their regionalization is fundamental for the assessment of ungauged basins. In comparison with calibrated rainfall–runoff models, statistical methods provide data-driven estimates representing a useful benchmark. The objective of this work is the interpolation of FDCs from ~500 discharge gauging stations in the Danube. To this aim we use total negative deviation top-kriging (TNDTK), as multi-regression models are shown to be unsuitable for representing FDCs across all durations and sites. TNDTK shows a high accuracy for the entire Danube region, with overall Nash-Sutcliffe efficiency values computed in a leave-p-out cross-validation scheme (p equal to one site, one-third and half of the sites), all above 0.88. A reliability measure based on kriging variance is attached to each interpolated FDC at ~4000 prediction nodes. The GIS layer of regionalized FDCs is made available for broader use in the region.     

New fuzzy neural network–Markov model and application in mid- to long-term runoff forecast

ABSTRACT

In this paper, a mid- to long-term runoff forecast model is developed using an ideal point fuzzy neural network–Markov (NFNN-MKV) hybrid algorithm to improve the forecasting precision. Combining the advantages of the new fuzzy neural network and the Markov prediction model, this model can solve the problem of stationary or volatile strong random processes. Defined error statistics algorithms are used to evaluate the performance of models. A runoff prediction for the Si Quan Reservoir is made by utilizing the modelling method and the historical runoff data, with a comprehensive consideration of various runoff-impacting factors such as rainfall. Compared with the traditional fuzzy neural networks and Markov prediction models, the results show that the NFNN-MKV hybrid algorithm has good performance in faster convergence, better forecasting accuracy and significant improvement of neural network generalization. The absolute percentage error of the NFNN-MKV hybrid algorithm is less than 7.0%, MSE is less than 3.9, and qualification rate reaches 100%. For further comparison of the proposed model, the NFNN-MKV model is employed to estimate (training and testing for 120-month-ahead prediction) and predict river discharge for 156 months at Weijiabao on the Weihe River in China. Comparisons among the results of the NFNN-MKV model, the WNN model and the SVR model indicate that the NFNN-MKV model is able to significantly increase prediction accuracy.

Editor D. Koutsoyiannis; Associate editor Y. Gyasi-Agyei     

A Modified Neural Network for Improving River Flow Prediction/Un Réseau de Neurones Modifié pour Améliorer la Prévision de L'Écoulement Fluvial

Abstract

Artificial neural network (ANN) models provide huge potential for simulating nonlinear behaviour of hydrological systems. However, the potential of ANN is yet to be fully exploited due to the problems associated with improving the model generalization performance. Generalization refers to the ability of a neural network to correctly process input data that have not been used for calibrating the neural network model. In the hydrological context, better generalization performance implies higher precision of forecasting. The primary objectives of this study are to explore new measures for improving the generalization performance of an ANN-based rainfall–runoff model, and to evaluate the applicability of the new measures. A modified neural network model (entitled goal programming (GP) neural network) for modelling the rainfall–runoff process has been developed, in which three enhancements are made as compared to the widely-used backpropagation (BP) network. The three enhancements are (a) explicit integration of hydrological prior knowledge into the neural network learning; (b) incorporation of a modified training objective function; and (c) reduction of network sensitivity to input errors. Seven watersheds across a range of climatic conditions and watershed areas in China were selected for examining the alternative networks. The results demonstrate that the GP consistently outperformed the BP both in the calibration and verification periods and three proposed measures yielded improvement of performance.     

Hydrological model comparison and assessment: criteria from catchment scales and temporal resolution

ABSTRACT

This study examines the performance of three hydrological models, namely the artificial neural network (ANN) model, the Hydrologiska Byråns Vattenbalansavdelning-D (HBV-D) model, and the Soil and Water Integrated Model (SWIM) over the upper reaches of the Huai River basin. The assessment is done by using databases of different temporal resolution and by further examining the applicability of SWIM for different catchment sizes. The results show that at monthly scale the performance of the ANN model is better than that of HBV-D and SWIM. The ANN model can be applied at any temporal scale as it establishes an artificial precipitation–runoff relationship for various time scales by only using monthly precipitation, temperature and runoff data. However, at daily scale the performance of both HBV-D and SWIM are similar or even better than the ANN model. In addition, the performance of SWIM at a small catchment size (less than 10 000 km²) is much better than at a larger catchment size. In view of climate change modelling, HBV-D and SWIM might be integrated in a dynamical atmosphere-water-cycle modelling rather than the ANN model due to their use of observed physical links instead of artificial relations within a black box.

Editor D. Koutsoyiannis; Associate editor D. Hughes     

Artificial neural networks for daily rainfall—runoff modelling

Abstract

The application of artificial neural network (ANN) methodology for modelling daily flows during monsoon flood events for a large size catchment of the Narmada River in Madhya Pradesh (India) is presented. The spatial variation of rainfall is accounted for by subdividing the catchment and treating the average rainfall of each subcatchment as a parallel and separate lumped input to the model. A linear multiple-input single-output (MISO) model coupled with the ANN is shown to provide a better representation of the rainfall-runoff relationship in such large size catchments compared with linear and nonlinear MISO models. The present model provides a systematic approach for runoff estimation and represents improvement in prediction accuracy over the other models studied herein.     

Modelling spatio-temporal patterns of sediment yield and connectivity in a semi-arid catchment with the WASA-SED model

Abstract

Rainfall–runoff induced soil erosion causes important environmental degradation by reducing soil fertility and impacting on water availability as a consequence of sediment deposition in surface reservoirs used for water supply, particularly in semi-arid areas. However, erosion models developed on experimental plots cannot be directly applied to estimate sediment yield at the catchment scale, since sediment redistribution is also controlled by the transport conditions along the landscape. In particular, representation of landscape connectivity relating to sediment transfer from upslope areas to the river network is required. In this study, the WASA-SED model is used to assess the spatial and temporal patterns of water and sediment connectivity for a semi-arid meso-scale catchment (933 km²) in Brazil. It is shown how spatial and temporal patterns of sediment connectivity within the catchment change as a function of landscape and event characteristics. This explains the nonlinear catchment response in terms of sediment yield at the outlet.

Citation Medeiros, P. H. A., Güntner, A., Francke, T., Mamede, G. L. & de Araújo, J. C. (2010) Modelling spatio-temporal patterns of sediment yield and connectivity in a semi-arid catchment with the WASA-SED model. Hydrol. Sci. J. 55(4), 636–648.     

Improving kriging of groundwater level data using nonlinear normalizing transformations—a field application

Abstract

The present research study investigates the application of nonlinear normalizing data transformations in conjunction with ordinary kriging (OK) for the accurate prediction of groundwater level spatial variability in a sparsely-gauged basin. We investigate three established normalizing methods, Gaussian anamorphosis, trans-Gaussian kriging and the Box-Cox method to improve the estimation accuracy. The first two are applied for the first time to groundwater level data. All three methods improve the mean absolute prediction error compared to the application of OK to the non-transformed data. In addition, a modified Box-Cox transformation is proposed and applied to normalize the hydraulic heads. The modified Box-Cox transformation in conjunction with OK is found to be the optimal spatial model based on leave-one-out cross-validation. The recently established Spartan semivariogram family provides the optimal model fit to the transformed data. Finally, we present maps of the groundwater level and the kriging variance based on the optimal spatial model.

Editor D. Koutsoyiannis; Associate editor A. Montanari

Citation Varouchakis, E.A., Hristopoulos, D.T., and Karatzas, G.P., 2012. Improving kriging of groundwater level data using nonlinear normalizing transformations—a field application. Hydrological Sciences Journal, 57 (7), 1404–1419.     

Performance of the COSERO precipitation–runoff model under non-stationary conditions in basins with different climates

Abstract

This study is a contribution to a model intercomparison experiment initiated during a workshop at the 2013 IAHS conference in Göteborg, Sweden. We present discharge simulations with the conceptual precipitation–runoff model COSERO in 11 basins located under different climates in Europe, Africa and Australia. All of the basins exhibit some form of non-stationary conditions, due, for example, to warming, droughts or land-cover change. The evaluation of the daily discharge simulations focuses on the overall model performance and its decomposition into three components measuring temporal dynamics, mean flow volume and distribution of flows. Calibration performance is similarly high as in previous COSERO applications. However, when looking at evaluation periods independent of the calibration, the model performance drops considerably, mainly due to severely biased discharge simulations in semi-arid basins with strong non-stationarity in rainfall. Simulations are more robust in European basins with humid climates. This highlights the fact that hydrological models frequently fail when simulations are required outside of calibration conditions in basins with non-stationary conditions. As a consequence, calibration periods should be sufficiently long to include both wet and dry periods, which should yield more robust predictions.     

Prediction and modelling of rainfall–runoff during typhoon events using a physically-based and artificial neural network hybrid model

Abstract

The accurate prediction of hourly runoff discharge in a watershed during heavy rainfall events is of critical importance for flood control and management. This study predicts n-h-ahead runoff discharge in the Sandimen basin in southern Taiwan using a novel hybrid approach which combines a physically-based model (HEC-HMS) with an artificial neural network (ANN) model. Hourly runoff discharge data (1200 datasets) from seven heavy rainfall events were collected for the model calibration (training) and validation. Six statistical indicators (i.e. mean absolute error, root mean square error, coefficient of correlation, error of time to peak discharge, error of peak discharge and coefficient of efficiency) were employed to evaluate the performance. In comparison with the HEC-HMS model, the single ANN model, and the time series forecasting (ARMAX) model, the developed hybrid HEC-HMS–ANN model demonstrates improved accuracy in recursive n-h-ahead runoff discharge prediction, especially for peak flow discharge and time.     

A comparative study on a simple two-parameter monthly water balance model and the Kajiyama formula for monthly runoff estimation

ABSTRACT

A two-parameter monthly water balance model to simulate runoff can be used for a water resources planning programme and climate impact studies. However, the model estimates two parameters of transformation of time scale (c) and of the field capacity (SC) by a trial-and-error method. This study suggests a modified methodology to estimate the parameters c and SC using the meteorological and geological conditions. The modified model is compared with the Kajiyama formula to simulate the runoff in the Han River and International Hydrological Programme representative basins in South Korea. We show that the estimated c and SC can be used as the initial or optimal values for the monthly runoff simulation study in the model.

EDITOR M.C. Acreman; ASSOCIATE EDITOR S. Kanae     

Regional runoff prediction through aggregation of first-order hydrological process knowledge: a case study

Abstract

The dominant processes concept was used to develop a regionally applicable rainfall—runoff model. The first-order runoff processes are identified through a combination of field investigations, physico-geographical analysis of the research area, the Alzette River basin in the Grand-Duchy of Luxembourg, and discharge data series analysis. Lithology appeared to be the major source of discrepancy in hydrological behaviour over the total area. As a result, the hydrological behaviour of each lithological substratum was characterized and conceptualized into a parsimonious model structure. The runoff signals were calibrated against the hourly-recorded discharge series of eight sub-basins, with parameter sensitivity and correlation analysis outlining the need for minor corrections to the model structure. Validation against another set of 10 sub-basins showed good results for the regional parameter set, with an average loss in efficiency (R_eff) of 0.04, compared to the reference model, with a mean R_eff of 0.79. Due to an up-scaling effect, inducing variations in the dominance of particular runoff processes, some anomalies were found in the performance of individual runoff characteristics. In this respect, limiting the application of the model to a certain spatial scale gives a high reliability of the prediction of the dynamics of hourly runoff in ungauged basins within the study area.