Improving the calibration strategy of the physically-based model WaSiM-ETH using critical events

Abstract

The use of a physically-based hydrological model for streamflow forecasting is limited by the complexity in the model structure and the data requirements for model calibration. The calibration of such models is a difficult task, and running a complex model for a single simulation can take up to several days, depending on the simulation period and model complexity. The information contained in a time series is not uniformly distributed. Therefore, if we can find the critical events that are important for identification of model parameters, we can facilitate the calibration process. The aim of this study is to test the applicability of the Identification of Critical Events (ICE) algorithm for physically-based models and to test whether ICE algorithm-based calibration depends on any optimization algorithm. The ICE algorithm, which uses the data depth function, was used herein to identify the critical events from a time series. Low depth in multivariate data is an unusual combination and this concept was used to identify the critical events on which the model was then calibrated. The concept is demonstrated by applying the physically-based hydrological model WaSiM-ETH on the Rems catchment, Germany. The model was calibrated on the whole available data, and on critical events selected by the ICE algorithm. In both calibration cases, three different optimization algorithms, shuffled complex evolution (SCE-UA), parameter estimation (PEST) and robust parameter estimation (ROPE), were used. It was found that, for all the optimization algorithms, calibration using only critical events gave very similar performance to that using the whole time series. Hence, the ICE algorithm-based calibration is suitable for physically-based models; it does not depend much on the kind of optimization algorithm. These findings may be useful for calibrating physically-based models on much fewer data.

Editor D. Koutsoyiannis; Associate editor A. Montanari

Citation Singh, S.K., Liang, J.Y., and Bárdossy, A., 2012. Improving calibration strategy of physically-based model WaSiM-ETH using critical events. Hydrological Sciences Journal, 57 (8), 1487–1505.     

Données isotopiques (8786/Sr) et changements hydrologiques depuis 15 000 ans sur l'Altiplano andin

Abstract

In catchments characterized by spatially varying hydrological processes and responses, the optimal parameter values or regions of attraction in parameter space may differ with location-specific characteristics and dominating processes. This paper evaluates the value of semi-distributed calibration parameters for large-scale streamflow simulation using the spatially distributed LISFLOOD model. We employ the Shuffled Complex Evolution Metropolis (SCEM-UA) global optimization algorithm to infer the calibration parameters using daily discharge observations. The resulting posterior parameter distribution reflects the uncertainty about the model parameters and forms the basis for making probabilistic flow predictions. We assess the value of semi-distributing the calibration parameters by comparing three different calibration strategies. In the first calibration strategy uniform values over the entire area of interest are adopted for the unknown parameters, which are calibrated against discharge observations at the downstream outlet of the catchment. In the second calibration strategy the parameters are also uniformly distributed, but they are calibrated against observed discharges at the catchment outlet and at internal stations. In the third strategy a semi-distributed approach is adopted. Starting from upstream, parameters in each subcatchment are calibrated against the observed discharges at the outlet of the subcatchment. In order not to propagate upstream errors in the calibration process, observed discharges at upstream catchment outlets are used as inflow when calibrating downstream subcatchments. As an illustrative example, we demonstrate the methodology for a part of the Morava catchment, covering an area of approximately 10 000 km². The calibration results reveal that the additional value of the internal discharge stations is limited when applying a lumped parameter approach. Moving from a lumped to a semi-distributed parameter approach: (i) improves the accuracy of the flow predictions, especially in the upstream subcatchments; and (ii) results in a more correct representation of flow prediction uncertainty. The results show the clear need to distribute the calibration parameters, especially in large catchments characterized by spatially varying hydrological processes and responses.     

The applicability of a parsimonious model for local and regional prediction of runoff

Abstract

The Hydrological Recursive Model (HRM), a conceptual rainfall-runoff model, was applied for local and regional simulation of hourly discharges in the transnational Alzette River basin (Luxembourg-France-Belgium). The model was calibrated for a range of various sub-basins with a view to analysing its ability to reproduce the variability of basin responses during flood generation. The regionalization of the model parameters was obtained by fitting simultaneously the runoff series of calibration sub-basins after their spatial discretization in lithological contrasting isochronal zones. The runoff simulations of the model agreed well with the recorded runoff series. Significant correlations with some basin characteristics and, noticeably, the permeability of geological formations, could be found for two of the four free model parameters. The goodness of fit for runoff predictions using the derived regional parameter set was generally satisfactory, particularly for the statistical characteristics of streamflow. A more physically-based modelling approach, or at least an explicit treatment of quick surface runoff, is expected to give better results for high peak discharge.     

Performance of temperature and radiation index models for point-scale snow water equivalent (SWE) simulations in the Moroccan High Atlas Mountains

ABSTRACT

Four models of increasing complexity were tested and compared to simulate snow water equivalent at the local scale in the Moroccan High Atlas range. A classical temperature index model (TI) and three enhanced temperature index models that respectively include the potential clear-sky direct radiation (HTI), the incoming solar radiation (ETI-A) and net solar radiation (ETI-B), were subjected to annual and multi-annual calibration and cross-validated over the period 2003–2010. When calibrated yearly, the ETI models could be better transferred to other years, whereas all models, including the simple TI model, showed good transferability when calibrated over a longer period that includes inter-annual climate variability. No strong and recurrent relationships emerged between yearly calibrated model parameters and annual climate conditions. However, strong parameter compensation was observed for the enhanced models, which can be explained partly by the collinearity of air temperature and solar radiation causing equifinality of model parameters.     

A multi-objective calibration approach using in-situ soil moisture data for improved hydrological simulation of the Prairies

ABSTRACT

Traditionally, hydrological models are only calibrated to reproduce streamflow regime without considering other hydrological state variables, such as soil moisture and evapotranspiration. Limited studies have been performed on constraining the model parameters, despite the fact that the presence of a large number of parameters may provide large degree of freedom, resulting in equifinality and poor model performance. In this study, a multi-objective optimization approach is adopted, and both streamflow and soil moisture data are calibrated simultaneously for an experimental study basin in the Saskatchewan Prairies in western Canada. The results of this study show that the multi-objective calibration improves model fidelity compared to the single objective calibration. Moreover, the study demonstrates that single objective calibration performed against only streamflow can fairly mimic the streamflow hydrograph but does not yield realistic estimation of other fluxes such as evapotranspiration and soil moisture (especially in deeper soil layers).     

György Kovacs 1925–1988

Abstract

There has been a trend in recent years towards the development and popularity of physically-based deterministic models. However, the application of such models is not without difficulties. This paper investigates the usefulness of a conceptual single-event model for simulating floods from catchments covering a wide variety of climatic and physiographic areas. The model has been calibrated on a group of catchments and the calibrated parameter values related to physical catchment indices. The resulting quantitative relationships are assessed with respect to their value for estimating the parameter values of the model when calibration is not possible. The results indicate that the technique is likely to provide flood estimations for medium sized catchments (5–150 km²) that are more reliable than several flood estimation methods currently in use in South Africa.     

Simulation of surface runoff and sediment yield using the water erosion prediction project (WEPP) model: a study in Kaneli watershed,Himalaya, India / Simulation de ruissellement de surface et d'érosion à l'aide du modèle WEPP: cas du bassin versant de Kaneli,Himalaya, Inde

Abstract

The Water Erosion Prediction Project (WEPP) model was calibrated and evaluated for estimation of runoff and sediment yield in the data-scarce conditions of the Indian Himalaya. The inputs derived from remote sensing and geographic information system technologies were combined in the WEPP modelling system to simulate surface runoff and sediment yield from the hilly Kaneli watershed. The model parameters were calibrated using measured data on runoff volumes and sediment yield. The calibrated model was validated by producing the monthly runoff and sediment yield simulations and comparing them with data that were not used in calibration. The model was also used to make surface runoff and sediment yield simulations for each of the individual watershed elements, comprising 18 hillslopes and seven channels, and the detailed monthly results for each are presented. Although, no field data on hillslope runoff and sediment yield are currently available for the validation of distributed results produced by the model, the present investigation has demonstrated clearly the applicability of the WEPP model in predicting hydrological variables in a data-scarce situation.     

A simplified climate model with combined atmospheric-hydrological processes

Abstract

Global climate change can be reproduced in detail by using three-dimensional general circulation models (GCMs). However, such complex models require super-computers and extensive hours of computational time for a single attempt at reproducing long term climate change. An alternative approach is to make simplifying assumptions that retain the essential physics for the desired simulation. Energy balance and Radiative-convective models are examples of such models. The model in this study follows the simplified approach using physics-based climate processes as well as interactions between atmospheric and hydrological processes. The vertically and latitudinally averaged mean temperature and mean water vapour content between 30°N-50°N latitudes are considered as atmospheric state variables while soil and sea temperatures and water storage amount are considered for describing the behaviour of the hydrological system. Temperatures in both the atmosphere and ground are calculated by a thermal energy equation that considers the physically-based processes of shortwave radiation, longwave radiation, sensible heat flux, and latent heat flux. Precipitation and evaporation processes transport moisture between the atmosphere and ground. In this study, the radiation parameterization of the simplified climate model is tested in the investigation of the various effects of global warming due to doubling and quadrupling of CO₂. Changes of temperature, soil water content, evaporation rate and precipitation rate are investigated by numerical experiments. The simplified climate model provides acceptable simulation of climate change and holds promise for practical investigations such as the interactions of physical processes in the evolution of drought phenomena.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

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.     

SUR LA PUISSANCE DES CRUES AU JAPON A MADAGASCAR ET A LA REUNION

Abstract

Hydrological modelling has faced the problem of ungauged basins for many years: how does one estimate hydrological characteristics for a river for which there are no data? Whatever the kind of model, it needs at least hydroclimatic input data and discharge data for calibration. However, the Yates model does not need any discharge data for calibration: it is a pre-calibrated model from a vegetation—climate classification map. In the specific context of West and Central Africa, where data are often of poor quality and very scarce, it is interesting to compare the performance of such a model with those of calibrated models, and with observed data. For this study, a platform including different semi-global rainfall—runoff models which allow the estimation of monthly runoff at a spatial resolution of 0.5° × 0.5° was used. The performance of the Yates model is very close to those of calibrated models, so that one can say that this simple model, based simply on a vegetation—climate classification, can be a very useful prediction tool in regions of scarce and unreliable data, such as those of interest to the International Association of Hydrological Sciences (IAHS) initiative on prediction in ungauged basins (PUB). Therefore, this model was applied to a period covering the last 30 years, and to a data set covering the first decades of the 21st century, from a climatic scenario of doubling the CO₂ concentration in the atmosphere. The results show that, in West Africa, where drought conditions have now prevailed for 35 years, water resources should still be decreasing in the future, following the general decreasing trend of rainfall projected by the climatic scenarios.     

Can satellite land surface temperature data be used similarly to river discharge measurements for distributed hydrological model calibration?

Abstract

A new methodology is proposed for the calibration of distributed hydrological models at the basin scale by constraining an internal model variable using satellite data of land surface temperature (LST). The model algorithm solves the system of energy and mass balances in terms of a representative equilibrium temperature that governs the fluxes of energy and mass over the basin domain. This equilibrium surface temperature, which is a critical model state variable, is compared to operational satellite LST, while calibrating soil hydraulic parameters and vegetation variables differently in each pixel, minimizing the errors. This procedure is compared to the traditional calibration using only discharge measurements. The distributed energy water balance model, Flash-flood Event-based Spatially-distributed rainfall–runoff Transformation – Energy Water Balance model (FEST-EWB), is used to test this approach. This methodology is applied to the Upper Yangtze River basin (China) using MODIS LST retrieved from satellite data in the framework of the NRSCC-ESA DRAGON-2 Programme. The calibration procedure based on LST seems to outperform the calibration based on discharge, with lower relative error and higher Nash-Sutcliffe efficiency index on cumulated volume.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

An approach to using snow areal depletion curves inferred from MODIS and its application to land surface modelling in Alaska

Snowcover areal depletion curves inferred from the moderate resolution imaging spectroradiometer (MODIS) are validated and then applied in NASA's catchment‐based land surface model (CLSM) for numerical simulations of hydrometeorological processes in the Kuparuk River basin (KRB) of Alaska. The results demonstrate that the MODIS snowcover fraction f derived from a simple relationship in terms of the normalized difference snow index compares well with Landsat values over the range 20 ≤ f ≤ 100%. For f < 20%, however, MODIS 500 m subpixel data underestimate the amount of snow by up to 13% compared with Landsat at spatial resolutions of 30 m binned to equivalent 500 m pixels. After a bias correction, MODIS snow areal depletion curves during the spring transition period of 2002 for the KRB exhibit similar features to those derived from surface‐based observations. These results are applied in the CLSM subgrid‐scale snow parameterization that includes a deep and a shallow snowcover fraction. Simulations of the evolution of the snowpack and of freshwater discharge rates for the KRB over a period of 11 years are then analysed with the inclusion of this feature. It is shown that persistent snowdrifts on the arctic landscape, associated with a secondary plateau in the snow areal depletion curves, are hydrologically important. An automated method is developed to generate the shallow and deep snowcover fractions from MODIS snow areal depletion curves. This provides the means to apply the CLSM subgrid‐scale snow parameterization in all watersheds subject to seasonal snowcovers. Improved simulations and predictions of the global surface energy and water budgets are expected with the incorporation of the MODIS snow data into the CLSM. Copyright © 2005 John Wiley & Sons, Ltd.     

Impact of model simplifications on soil erosion predictions: application of the GLUE methodology to a distributed event‐based model at the hillslope scale

In this paper, we analyse how the performance and calibration of a distributed event‐based soil erosion model at the hillslope scale is affected by different simplifications on the parameterizations used to compute the production of suspended sediment by rainfall and runoff. Six modelling scenarios of different complexity are used to evaluate the temporal variability of the sedimentograph at the outlet of a 60 m long cultivated hillslope. The six scenarios are calibrated within the generalized likelihood uncertainty estimation framework in order to account for parameter uncertainty, and their performance is evaluated against experimental data registered during five storm events. The Nash–Sutcliffe efficiency, percent bias and coverage performance ratios show that the sedimentary response of the hillslope in terms of mass flux of eroded soil can be efficiently captured by a model structure including only two soil erodibility parameters, which control the rainfall and runoff production of suspended sediment. Increasing the number of parameters makes the calibration process more complex without increasing in a noticeable manner the predictive capability of the model. Copyright © 2015 John Wiley & Sons, Ltd.     

Pilot-point based multi-objective calibration in a surface–subsurface distributed hydrological model

Abstract

In this study, a fully-coupled surface–subsurface, distributed, physics-based hydrological model was calibrated using the pilot-point method. A minimum variance field rule was included in the objective function to regularize the extensive calibration exercise that included 74 parameters (72 associated with pilot points and two spatially-invariant channel parameters). Because the overland and vadose zone systems are not in permanent hydrological connection, the information contained in the observation points may not be accessible by the pilot points at all times, rendering them insensitive to the observations and hindering the calibration process. An analysis of the spatial and temporal variability of parameter sensitivities was done to explore how the information contained in local observations spreads from the observation points to the pilot points, where parameter values are identified. The results show that the channel flow time series is valuable to identify the parameters at all pilot-point locations, indicating that the information in channel flow propagates to the entire basin. However, information in soil moisture measurements is of local extent and thus only valuable to identify the parameters at locations close to the observation point.

Editor D. Koutsoyiannis; Associate editor I. Nalbantis

Citation Maneta, M.P. and Wallender, W.W., 2013. Pilot-point based multi-objective calibration in a surface–subsurface distributed hydrological model. Hydrological Sciences Journal, 58 (2), 390–407.     

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.     

Impact of warming climate on water management for the Ariège River basin (France)

Abstract

The French national project IMAGINE2030 aims to assess future water availability in the Garonne River basin (southwest France) by taking account of changes in both climate and water management in the 2030s. Within this project, two mountainous drainage basins located in the Pyrenees were examined to assess the specific impact of climate change on reservoir management. The Salat River basin at Roquefort, is considered as a proxy (representative of a natural basin), whereas the Ariège River at Foix is influenced by hydropower production in winter and by water releases to sustain low flows in summer. The Cequeau rainfall–runoff model, combined with a simplified model of reservoir management operations, was calibrated on present-day conditions and forced with climate projections derived from the IPCC AR4 report. The results show that a warming climate over the basins induces a decrease in mean annual runoff, a shift to earlier snow melting in mountainous areas and more severe low-flow conditions. The simulations show a decrease in electricity generation. Under two water management scenarios (one “business-as-usual” and the other incorporating an increased downstream water demand in compliance with requirements for increased minimum flow), simulations for the Ariège River basin suggest an earlier filling of the reservoir is necessary in winter to anticipate the increased release from reservoirs in summer to support minimum flow farther downstream.

Editor Z.W. Kundzewicz; Associate editor D. Hughes

Citation Hendrickx, F. and Sauquet, E., 2013. Impact of warming climate on water management for the Ariège River basin (France). Hydrological Sciences Journal, 58 (5), 976–993.     

Simulating the daily discharge of the Mandovi River,west coast of India

Abstract

A hydrological modelling framework was assembled to simulate the daily discharge of the Mandovi River on the Indian west coast. Approximately 90% of the west-coast rainfall, and therefore discharge, occurs during the summer monsoon (June–September), with a peak during July–August. The modelling framework consisted of a digital elevation model (DEM) called GLOBE, a hydrological routing algorithm, the Terrestrial Hydrological Model with Biogeochemistry (THMB), an algorithm to map the rainfall recorded by sparse raingauges to the model grid, and a modified Soil Conservation Service Curve Number (SCS-CN) method. A series of discharge simulations (with and without the SCS method) was carried out. The best simulation was obtained after incorporating spatio-temporal variability in the SCS parameters, which was achieved by an objective division of the season into five regimes: the lean season, monsoon onset, peak monsoon, end-monsoon, and post-monsoon. A novel attempt was made to incorporate objectively the different regimes encountered before, during and after the Indian monsoon, into a hydrological modelling framework. The strength of our method lies in the low demand it makes on hydrological data. Apart from information on the average soil type in a region, the entire parameterization is built on the basis of the rainfall that is used to force the model. That the model does not need to be calibrated separately for each river is important, because most of the Indian west-coast basins are ungauged. Hence, even though the model has been validated only for the Mandovi basin, its potential region of application is considerable. In the context of the Prediction in Ungauged Basins (PUB) framework, the potential of the proposed approach is significant, because the discharge of these (ungauged) rivers into the eastern Arabian Sea is not small, making them an important element of the local climate system.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi

Citation Suprit, K., Shankar, D., Venugopal, V. and Bhatkar, N.V., 2012. Simulating the daily discharge of the Mandovi River, west coast of India. Hydrological Sciences Journal, 57 (4), 686–704.     

Soil hydraulic properties for moisture redistribution in a large-scale heterogeneous landscape

Abstract

Guidelines of effective soil hydraulic parameters were developed to be applicable in simulating average infiltration and subsequent moisture redistribution over a large-scale heterogeneous field. Average large-scale infiltration and redistribution in heterogeneous soils were quantified through multiple simulations of local-scale processes. The effective hydraulic parameters were derived to simulate the average amount of infiltrating water, and to capture the subsequent surface soil moisture redistribution averaged over the large heterogeneous landscape. The results demonstrated that the effective hydraulic parameters typically exhibited a step change from infiltration to redistribution, with the size of the step change being related to the degree of hydraulic parameter heterogeneity and the correlations among the hydraulic parameters. However, the effective hydraulic parameters did not change significantly over time for the moisture redistribution. It was further demonstrated that the size of the step change was smallest for effective saturated hydraulic conductivity.

Editor Z.W. Kundzewicz; Associate editor Y. Guttman

Citation Zhu, J.T. and Sun, D.M., 2012. Soil hydraulic properties for moisture redistribution in a large-scale heterogeneous landscape. Hydrological Sciences Journal, 57 (6), 1196–1206.