Probabilistic calibration of a distributed hydrological model for flood forecasting

Abstract

The complexity of distributed hydrological models has led to improvements in calibration methodologies in recent years. There are various manual, automatic and hybrid methods of calibration. Most use a single objective function to calculate estimation errors. The use of multi-objective calibration improves results, since different aspects of the hydrograph may be considered simultaneously. However, the uncertainty of estimates from a hydrological model can only be taken into account by using a probabilistic approach. This paper presents a calibration method of probabilistic nature, based on the determination of probability functions that best characterize different parameters of the model. The method was applied to the Real-time Interactive Basin Simulator (RIBS) distributed hydrological model using the Manzanares River basin in Spain as a case study. The proposed method allows us to consider the uncertainty in the model estimates by obtaining the probability distributions of flows in the flood hydrograph.

Citation Mediero, L., Garrote, L. & Martín-Carrasco, F. J. (2011) Probabilistic calibration of a distributed hydrological model for flood forecasting. Hydrol. Sci. J. 56(7), 1129–1149.     

Development of a distributed hydrological model to facilitate watershed management

To facilitate precise and cost-effective watershed management, a simple yet spatially and temporally distributed hydrological model (DHM-WM) was developed. The DHM-WM is based on the Mishra-Singh version of the curve number method, with several modifications: The spatial distribution of soil moisture was considered in moisture updating; the travel time of surface runoff was calculated on a grid cell basis for routing; a simple tile flow module was included as an option. The DHM-WM was tested on a tile-drained agricultural watershed in Indiana, USA. The model with the tile flow module performed well in the study area, providing a balanced water budget and reasonable flow partitioning. The daily coefficient of determination and Nash-Sutcliffe coefficient were 0.58 and 0.56, for the calibration period, and 0.63 and 0.62 for the validation period. The DHM-WM also provides detailed information about the source areas of flow components, the travel time and pathways of surface runoff.

EDITOR A. Castellarin; ASSOCIATE EDITOR F.-J. Chang     

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.     

Development of distributed time-variant gain model for nonlinear hydrological systems

Hydrological science is a branch of the earth sci-ences.To study the complexities of hydrological pro-cesses and the associated environmental problems,a systematic approach is desired.For instance,Dooge published Linear Theory of Hydrological System[1]in1973.Singh(1988)published Hydrological Systems[2],with its Chinese version[3]being translated by the Yel-low River Conservancy Commission in2000.Ge(1999)carried out systematic studies on the hydro-logical linear system theory,and produced …     

Development and application of a distributed catchment-scale hydrological model for the River Ythan,NE Scotland

Mathematical models are being used to develop a decision support system for integrated management of the Ythan catchment in NE Scotland. One component of this has involved the development of a distributed catchment-scale hydrological model. The model is based on subsurface flow routing and calculates the contribution to stream flow from each 50 m×50 m cell in the 548 km² catchment. It uses two topographic parameters, slope and distance to stream following the main line of flow, and five physical parameters. The topographic analysis and distributed flow accumulation are performed by linking the single cell model with a geographic information system. Preliminary results from a three-year simulation of daily flows indicate that the model successfully predicts the main characteristics of the catchment flow. © 1998 John Wiley & Sons, Ltd.     

A BTOP model to extend TOPMODEL for distributed hydrological simulation of large basins

Topography is a dominant factor in hillslope hydrology. TOPMODEL, which uses a topographical index derived from a simplified steady state assumption of mass balance and empirical equations of motion over a hillslope, has many advantages in this respect. Its use has been demonstrated in many small basins (catchment areas of the order of 2–500 km²) but not in large basins (catchment areas of the order of 10 000–100 000 km²). The objective of this paper is to introduce the Block‐wise TOPMODEL (BTOP) as an extension of the TOPMODEL concept in a grid based framework for distributed hydrological simulation of large river basins. This extension was made by redefining the topographical index by using an effective contributing area af(a) (0?f(a)?1) per unit grid cell area instead of the upstream catchment area per unit contour length and introducing a concept of mean groundwater travel distance. Further the transmissivity parameter T₀ was replaced by a groundwater dischargeability D which can provide a link between hill slope hydrology and macro hydrology. The BTOP model uses all the original TOPMODEL equations in their basic form. The BTOP model has been used as the core hydrological module of an integrated distributed hydrological model YHyM with advanced modules of precipitation, evapotranspiration, flow routing etc. Although the model has been successfully applied to many catchments around the world since 1999, there has not been a comprehensive theoretical basis presented in such applications. In this paper, an attempt is made to address this issue highlighted with an example application using the Mekong basin. Copyright © 2007 John Wiley & Sons, Ltd.     

Lumped parameter sensitivity analysis of a distributed hydrological model within tropical and temperate catchments

Parameter sensitivity of the Distributed Hydrology‐Soil‐Vegetation Model (DHSVM) was studied in two contrasting environments: (1) Pang Khum Experimental Watershed (PKEW) in tropical northern Thailand; and (2) Cedar River basin (CRB) in Washington State of the temperate US Pacific Northwest. The analysis shows that for both basins, the most sensitive soil parameters were porosity, lateral saturated hydraulic conductivity, and the exponential decrease rate of lateral saturated hydraulic conductivity with soil depth. The most sensitive vegetation parameters were leaf area index, vegetation height, vapour pressure deficit, minimum stomatal resistance (for both grassland and forest scenarios), hemisphere fractional coverage, overstory fractional coverage, and trunk space (for the forest scenario only). Parameter sensitivity was basin‐specific, with the humid, temperate CRB being more influenced by vegetation parameters, while tropical PKEW was more influenced by soil properties. Increases and decreases in parameter values resulted in opposite and unequal changes in bias and root mean square error (RMSE), indicating the non‐linearity of physical process represented in the hydrological model. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial scale effect on the upper soil effective parameters of a distributed hydrological model

Nonlinear dynamics and spatial variability in hydrological systems make the formulation of scaling theories difficult. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter H_u and the saturated hydraulic conductivity parameter k_s from microscale (sub‐grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of H_u and k_s was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre‐established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the nonlinearity of the system from microscale to the mesoscale via non‐stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.     

A distributed monthly hydrological model for integrating spatial variations of basin topography and rainfall

Hydrological models at a monthly time‐scale are important tools for hydrological analysis, such as in impact assessment of climate change and regional water resources planning. Traditionally, monthly models adopt a conceptual, lumped‐parameter approach and cannot account for spatial variations of basin characteristics and climatic inputs. A large requirement for data often severely limits the utility of physically based, distributed‐parameter models. Based on the variable‐source‐area concept, we considered basin topography and rainfall to be two major factors whose spatial variations play a dominant role in runoff generation and developed a monthly model that is able to account for their influences in the spatial and temporal dynamics of water balance. As a hybrid of the Xinanjiang model and TOPMODEL, the new model is constructed by innovatively making use of the highly acclaimed simulation techniques in the two existing models. A major contribution of this model development study is to adopt the technique of implicit representation of soil moisture characteristics in the Xinanjiang model and use the TOPMODEL concept to integrate terrain variations into runoff simulation. Specifically, the TOPMODEL topographic index ln(a/tanβ) is converted into an index of relative difficulty in runoff generation (IRDG) and then the cumulative frequency distribution of IRDG is used to substitute the parabolic curve, which represents the spatial variation of soil storage capacity in the Xinanjiang model. Digital elevation model data play a key role in the modelling procedures on a geographical information system platform, including basin segmentation, estimation of rainfall for each sub‐basin and computation of terrain characteristics. Other monthly data for model calibration and validation are rainfall, pan evaporation and runoff. The new model has only three parameters to be estimated, i.e. watershed‐average field capacity WM, pan coefficient η and runoff generation coefficient α. Sensitivity analysis demonstrates that runoff is least sensitive to WM and, therefore, it can be determined by a prior estimation based on the climate and soil properties of the study basin. The other two parameters can be determined using optimization methods. Model testing was carried out in a number of nested sub‐basins of two watersheds (Yuanjiang River and Dongjiang River) in the humid region in central and southern China. Simulation results show that the model is capable of describing spatial and temporal variations of water balance components, including soil moisture content, evapotranspiration and runoff, over the watershed. With a minimal requirement for input data and parameterization, this terrain‐based distributed model is a valuable contribution to the ever‐advancing technology of hydrological modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessment on the function of reservoirs for flood control during typhoon seasons based on a distributed hydrological model

A grid‐based distributed hydrological model, PDTank model, is used to simulate hydrological processes in the upper Tone River catchment. The Tone River catchment often suffers from heavy rainfall events during the typhoon seasons. The reservoirs located in the catchment play an important role in flood regulation. Through the coupling of the PDTank model and a reservoir module that combines the storage function and operation function, the PDTank model is used for flood forecasting in this study. By comparing the hydrographs simulated using gauging and radar rainfall data, it is found that the spatial variability of rainfall is an important factor for flood simulation and the accuracy of the hydrographs simulated using radar rainfall data is slightly improved. The simulation of the typhoon flood event numbered No. 9 shows that the reservoirs in the catchment attenuate the peak flood discharge by 423·3 m³/s and validates the potential applicability of the distributed hydrological model on the assessment of function of reservoirs for flood control during typhoon seasons. Copyright © 2011 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.     

Coupling of a simultaneous heat and water model with a distributed hydrological model and evaluation of the combined model in a cold region watershed

Snow and frozen soil prevail in cold regions worldwide, and the integration of these processes is crucial in hydrological models. In this study, a combined model was developed by fully coupling a simultaneous heat and water model with a geomorphologically based distributed hydrological model. The combined model simulates vertical and lateral water transfer as well as vertical heat fluxes and is capable of representing the effects of frozen soil and snowmelt on hydrological processes in cold regions. This model was evaluated by using in situ observations in the Binggou watershed, an experimental watershed for cold region hydrology of the Watershed Allied Telemetry Experimental Research Project. Results showed that the model was able to predict soil freezing and thawing, unfrozen soil water content, and snow depth reasonably well. The simulated hydrograph was in good agreement with the in situ observation. The Nash–Sutcliffe coefficient of daily discharge was 0.744 for the entire simulation period, 0.472 from April to June, and 0.711 from June to November. This model can improve our understanding of hydrological processes in cold regions and assess the impacts of global warming on hydrological cycles and water resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Alternative climate data sources for distributed hydrological modelling on a daily time step

Two major criteria in choosing climate data for use in hydrological modelling are the period of record of the data set and the proximity of the collection platform(s) to the basin under study. Conventional data sets are derived from weather stations; however, in many cases there are no weather stations sufficiently close to a basin to be representative of climate conditions in that basin. In addition, it is often the case either that the period of record for the weather station(s) does not cover the period of the proposed simulation or that there are gaps in the data. Therefore, the objectives of this study are to investigate alternative climate data sources for use in hydrological modelling and to develop a protocol for creating hydrological data sets that are spatially and temporally harmonized. The methods we used for constructing daily, spatially distributed, climatic data sets of precipitation, maximum and minimum temperature, wind speed, solar radiation, potential evapotranspiration, and relative humidity are described. The model used in this study was the Soil and Water Assessment Tool implemented on the Mimbres River Basin located in southwestern New Mexico, USA, for the period 2003–2006. Our hydrological simulations showed that two events in January and February 2005 were missed, while an event in August 2006 was well simulated. We have also investigated the usefulness of several other precipitation data sets and compared the simulation results. Copyright © 2010 John Wiley & Sons, Ltd.     

Grid parameterization of a conceptual distributed hydrological model through integration of a sub-grid topographic index: necessity and practicability

Abstract

Grid-based distributed models have become popular for describing spatial hydrological processes. However, the influence of non-homogeneity within a grid on streamflow simulation was not adequately addressed in the literature. In this study, we investigated how the statistical characteristics of soil moisture storage within a grid impacts on streamflow simulations. The spatial variation of the topographic index, TI, within a grid was used to determine parameter B of the statistical curve of soil moisture storage in the Xinanjiang model. For comparison of influences of the non-homogeneity within a grid on streamflow simulation, two parameterization schemes of soil moisture storage capacity were developed: a grid-parameterization scheme for a distributed model and a catchment-averaged scheme for a semi-distributed model. The practicability and usefulness of the grid-parameterization method were evaluated through model comparisons. The two models were applied in Jiangwan experimental catchment Zhejiang Province, China. Streamflow discharge data at the catchment outlet from 1971 to 1986 at different temporal resolutions, e.g. 15 min and daily time step, were used for model calibration and validation. Statistical results for different grid scales demonstrated that the mean and variation of TI and B decline significantly as the grid scale increases. The simulated streamflow discharges of the two models were similar and the semi-distributed model outperformed the distributed model slightly when the streamflow at the outlet of the catchment was used as the only basis for comparison. In addition, a relatively larger bias in the predicted discharges between these two models was observed along with an abrupt increase of soil moisture saturation ratio. A further analysis of the simulated soil moisture content distribution revealed that the distributed model can provide a reasonable representation of the variable source area concept, which was justified to some extent by the field experiment data.

Editor D. Koutsoyiannis

Citation Liu, J.T., Chen, X., Wu, J.C., Zhang, X.N., Feng, D.Z. and Xu, C.-Y., 2012. Grid parameterization of a conceptual, distributed hydrological model through integration of a sub-grid topographic index: necessity and practicability. Hydrological Sciences Journal, 57 (2), 282–297.     

Simulating saturation‐excess surface run‐off in a semi‐distributed hydrological model

The level of complexity, and the number of parameters, to include in a hydrological model is a relatively contentious issue in hydrological modelling. However, it can be argued that explicitly representing important run‐off generation processes can improve the practical value of a model's outputs. This paper explores the benefits of including a new function into an existing semi‐distributed hydrological model (the Pitman model) that is widely used in the sub‐Saharan Africa region. The new function was designed to represent saturation‐excess surface run‐off processes at subcatchment scales and was motivated by the evidence of dambo (low topography riparian areas) type features in many sub‐Saharan river basins. The results for uncertainty versions of the model, with and without the new function, were compared for 25 catchments, which were divided up into those where evidence of dambos exists and those where there is no such evidence. The results suggest that the new function certainly improves the model results for the catchments where dambos exist, but not in situations where saturation‐excess surface run‐off is not expected to occur. The overall conclusion is therefore that the addition of the new function is justified.     

Development and application of a spatially‐distributed Arctic hydrological and thermal process model (ARHYTHM)

A process‐based, spatially distributed hydrological model was developed to quantitatively simulate the energy and mass transfer processes and their interactions within arctic regions (arctic hydrological and thermal model, ARHYTHM). The model first determines the flow direction in each element, the channel drainage network and the drainage area based upon the digital elevation data. Then it simulates various physical processes: including snow ablation, subsurface flow, overland flow and channel flow routing, soil thawing and evapotranspiration. The kinematic wave method is used for conducting overland flow and channel flow routing. The subsurface flow is simulated using the Darcian approach. The energy balance scheme was the primary approach used in energy‐related process simulations (snowmelt and evapotranspiration), although there are options to model snowmelt by the degree‐day method and evapotranspiration by the Priestley–Taylor equation. This hydrological model simulates the dynamic interactions of each of these processes and can predict spatially distributed snowmelt, soil moisture and evapotranspiration over a watershed at each time step as well as discharge in any specified channel(s). The model was applied to Imnavait watershed (about 2·2 km²) and the Upper Kuparuk River basin (about 146 km²) in northern Alaska. Simulated results of spatially distributed soil moisture content, discharge at gauging stations, snowpack ablations curves and other results yield reasonable agreement, both spatially and temporally, with available data sets such as SAR imagery‐generated soil moisture data and field measurements of snowpack ablation, and discharge data at selected points. The initial timing of simulated discharge does not compare well with the measured data during snowmelt periods mainly because the effect of snow damming on runoff was not considered in the model. Results from the application of this model demonstrate that spatially distributed models have the potential for improving our understanding of hydrology for certain settings. Finally, a critical component that led to the performance of this modelling is the coupling of the mass and energy processes. Copyright © 2000 John Wiley & Sons, Ltd.     

The value of environmental modelling languages for building distributed hydrological models

An evaluation is made of the suitability of programming languages for hydrological modellers to create distributed, process‐based hydrological models. Both system programming languages and high‐level environmental modelling languages are evaluated based on a list of requirements for the optimal programming language for such models. This is illustrated with a case study, implemented using the PCRaster environmental modelling language to create a distributed, process‐based hydrological model based on the concepts of KINEROS‐EUROSEM. The main conclusion is that system programming languages are not ideal for hydrologists who are not computer programmers because the level of thinking of these languages is too strongly related to specialized computer science. A higher level environmental modelling language is better in the sense that it operates at the conceptual level of the hydrologist. This is because it contains operators that identify hydrological processes that operate on hydrological entities, such as two‐dimensional maps, three‐dimensional blocks and time‐series. The case study illustrates the advantages of using an environmental modelling language as compared with system programming languages in fulfilling requirements on the level of thinking applied in the language, the reusability of the program code, the lack of technical details in the program, a short model development time and learnability. The study shows that environmental modelling languages are equally good as system programming languages in minimizing programming errors, but are worse in generic application and performance. It is expected that environmental modelling languages will be used in future mainly for development of new models that can be tailored to modelling aims and the field data available. Copyright © 2002 John Wiley & Sons, Ltd.