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.     

TOPMODEL在珠江流域布柳河流域的应用及其与新安江模型的比较

TOPMODEL是一种以地形为基础的半分布式流域水文模型.对珠江流域布柳河流域的DEM信息进行处理,提取流域的水系、子流域边界、地形指数及水流路径距离的分布,将TOPMODEL应用于该流域中.另外将新安江模型也应用于该流域进行比较.此外,分析了两种模型结构差异所带来的模拟功能差异.两种模型模拟结果精度差异不大,而TOPMODEL实现了空间产流面积分布的可视化.     

TOPMODEL: A critique

TOPMODEL (a TOPography based hydrological MODEL) is now 20 years old and has been the subject of numerous applications to a wide variety of catchments. This paper represents a critical review of some of the issues involved in application of the TOPMODEL concepts, including the basic assumptions involved; the derivation of topographic index distributions from digital terrain data; additional model components; meaning and calibration of the model parameters; and issues involved in model validation and predictive uncertainty. The aim is to provoke a thoughtful approach to hydrological modelling and the interaction of modelling and field work. Some recommendations are made for future modelling practice. © 1997 John Wiley & Sons, Ltd.     

The in(a/tan/β) index: How to calculate it and how to use it within the topmodel framework

Topographic indices may be used to attempt to approximate the likely distribution of variable source areas within a catchment. One such index has been applied widely using the distribution function catchment model, TOPMODEL, of Beven and Kirkby (1979). Validation of the spatial predictions of TOPMODEL may be affected by the algorithm used to calculate the model's topographic index. A number of digital terrain analysis (DTA) methods are therefore described for use in calculating the TOPMODEL topographic index, In(a/tanβ) (a = upslope contributing area per unit contour; tanβ = local slope angle). The spatial pattern and statistical distribution of the index is shown to be substantially different for different calculation procedures and differing pixel resolutions. It is shown that an interaction between hillslope contributing area accumulation and the analytical definition of the channel network has a major influence on calculated In(a/tanβ) index patterns. A number of DTA tests were performed to explore this interaction. The tests suggested that an ‘optimum’ channel initiation threshold (CIT) may be identified for positioning river headwaters in a raster digital terrain model (DTM). This threshold was found to be dependent on DTM grid resolution. Grid resolution is also suggested to have implications for the validation of spatial model predictions, implying that ‘optimum’ TOPMODEL parameter sets may be unique to the grid scale used in their derivation. Combining existing DTA procedures with an identified CIT, a procedure is described to vary the directional diffusion of contributing area accumulation with distance from the channel network.     

Investigation on the properties of the relationship between rare and extreme rainfall and flood volumes, under some distributional restrictions

The fact that rainfall data are usually more abundant and more readily regionalized than streamflow data has motivated hydrologists to conceive methods that incorporate the hydrometeorologial information into flood frequency analyses. Some of them, particularly those derived from the French GRADEX method, involve assumptions concerning the relationship between extreme rainfall and flood volumes, under some distributional restrictions. In particular, for rainfall probability distributions exhibiting exponential-like upper tails, it is possible to derive the shape and scale of the probability distribution of flood volumes by hypothesizing the basic properties of such a relationship, under rare and/or extreme conditions. This paper focuses on a parametric mathematical model for the relationship between rare and extreme rainfall and flood volumes under exponentially-tailed distributions. The model is analyzed and fitted to rare and extreme events derived from hydrological simulation of long stochastically-generated synthetic series of rainfall and evaporation for the Indaiá River basin, located in south-central Brazil. The paper also provides a sensitivity analysis of the model parameters in order to better understand flood events under rare and extreme conditions. By working with hydrologically plausible hypothetical events, the modeling approach proved to be a useful way to explore extraordinary rainfall and flood events. The results from this exploratory analysis provide grounds to derive some conclusions regarding the relative positions of the upper tails of the probability distributions of rainfall and flood volumes.     

Wavelet‐based regularization of the extracted topographic index from high‐resolution topography for hydro‐geomorphic applications

High‐resolution topography, e.g. 1‐m digital elevation model (DEM) from light detection and ranging (LiDAR), offers opportunity for accurate identification of topographic features of relevance for hydrologic and geomorphologic modelling. Yet, the computation of some derived topographic properties, such as the topographic index (TI), is characterized by daunting challenges that hamper the full exploration of topography‐based models. Particular problems, for example, arise when a distributed (or semi‐distributed) rainfall–runoff model is applied to high‐resolution DEMs. Indeed, the characteristic dependency between landscape resolution and the computed TI distribution results in the formation of un‐physical, unconnected saturated zones, which in turn cause unrealistic representations of rainfall–runoff dynamics. In this study, we present a methodology based on a multi‐resolution wavelet transformation that, by means of a soft‐thresholding scheme on the wavelet coefficients, filters the noise of high‐resolution topography to construct regularized sets of locally smoother topography on which the TI is computed. While the methodology needs a somewhat arbitrary definition of the wavelet coefficients threshold, our study shows that when the information content (entropy) of the TI distribution is used as a filtering efficiency metric, a critical threshold automatically emerges in the landscape reconstruction. The methodology is demonstrated using 1‐m LiDAR data for the Elder Creek River basin in California. While the proposed case study uses a TOPMODEL approach, the methodology can be extended to different topography‐based models and is not limited to hydrological applications. Copyright © 2012 John Wiley & Sons, Ltd.     

Semi‐distributed flood runoff model at the subcontinental scale for southwestern Iran

Rainfall–runoff modelling, as a surface hydrological process, on large‐scale data‐poor basins is currently a major topic of investigation that requires the model parameters be identified by using basin physical characteristics rather than calibration. This paper describes the application of the TOPMODEL framework accompanied by a kinematic wave model to the Karun River sub‐basins in southwestern Iran with just one conceptual parameter for calibration. ISLSCP1, HYDRO1K and Reynolds data sets are presented in a geographical information system and used as data sources for meteorological information, hydrological features and soil characteristics of the study area respectively. The results show that although the model developed can adequately predict flood runoff in the catchment with only one calibrated parameter, it is suggested that the effect of surface reservoirs be considered in the proposed model. Copyright © 2007 John Wiley & Sons, Ltd.     

Capturing the essential spatial variability in distributed hydrological modelling: Infiltration parameters

Selecting the correct resolution in distributed hydrological modelling at the watershed scale is essential in reducing scale-related errors. The work presented herein uses information content (entropy) to identify the resolution which captures the essential variability, at the watershed scale, of the infiltration parameters in the Green and Ampt infiltration equation. A soil map of the Little Washita watershed in south-west Oklahoma, USA was used to investigate the effects of grid cell resolution on the distributed modelling of infiltration. Soil-derived parameters and infiltration exhibit decreased entropy as resolutions become coarser. This is reflected in a decrease in the maximum entropy value for the reclassified/derived parameters vis a vis the original data. Moreover, the entropy curve, when plotted against resolution, shows two distinct segments: a constant section where no entropy was lost with decreasing resolution and another part which is characterized by a sharp decrease in entropy after a critical resolution of 1209 m is reached. This methodology offers a technique for assessing the largest cell size that captures the spatial variability of infiltration parameters for a particular basin. A geographical information system (GIS) based rainfall-runoff model is used to simulate storm hydrographs using infiltration parameter maps at different resolutions as inputs. Model results up to the critical resolution are reproducible and errors are small. However, at resolutions beyond the critical resolution the results are erratic with large errors. A major finding of this study is that a large resolution (1209 m for this basin) yields reproducible model results. When modelling a river basin using a distributed model, the resolution (grid cell size) can drastically affect the model results and calibration. The error structure attributable to grid cell resolution using entropy as a spatial variability measure is shown.     

A network‐index‐based version of TOPMODEL for use with high‐resolution digital topographic data

This paper describes the preliminary development of a network‐index approach to modify and to extend the classic TOPMODEL. Application of the basic Beven and Kirkby form of TOPMODEL to high‐resolution (2·0 m) laser altimetric data (based upon the UK Environment Agency's light detection and ranging (LIDAR) system) to a 13·8 km² catchment in an upland environment identified many saturated areas that remained unconnected from the drainage network even during an extreme flood event. This is shown to be a particular problem with using high‐resolution topographic data, especially over large appreciable areas. To deal with the hydrological consequences of disconnected areas, we present a simple network index modification in which saturated areas are only considered to contribute when the topographic index indicates continuous saturation through the length of a flow path to the point where the path becomes a stream. This is combined with an enhanced method for dealing with the problem of pits and hollows, which is shown to become more acute with higher resolution topographic data. The paper concludes by noting the implications of the research as presented for both methodological and substantive research that is currently under way. Copyright © 2004 John Wiley & Sons, Ltd.     

A MATLAB implementation of TOPMODEL

The MATLAB SIMULINK programming language is applied to the TOPMODEL rainfall–runoff model. SIMULINK requires a good recognition of model dynamics, which has been achieved here in a version based on the first TOPMODEL (Beven and Kirkby, 1979). Introducing the topographic index distribution in a vector form allows the generalization and simplification of the SIMULINK structure. The SIMULINK version of TOPMODEL has a very easy to understand graphical representation, which shows, in a straightforward way, all the physical interactions that take place in the model. Moreover, owing to its modular structure it is easy to add new and/or develop old submodels, depending on the available data and the goal of the modelling. In the example given here TOPMODEL was extended by two submodels representing the soil moisture and evaporation distribution in the catchment. Preparation of the data and presentation of the results is done in MATLAB. Discharge predictions and spatial patterns of hydrological response are demonstrated for a separate validation period. © 1997 John Wiley & Sons, Ltd.     

Flood simulation using different sources of rainfall in the Huong River,Vietnam / Simulation d'inondation à l'aide de différentes sources d'information pluviométrique dans le bassin de la Rivière Huong,Vietnam

Abstract

The possibility of simulating flooding in the Huong River basin, Vietnam, was examined using quantitative precipitation forecasts at regional and global scales. Raingauge and satellite products were used for observed rainfall. To make maximum use of the spatial heterogeneity of the different types of rainfall data, a distributed hydrological model was set up to represent the hydrological processes. In this way, streamflow simulated using the rainfall data was compared with that observed in situ. The forecast on a global scale showed better performance during normal flow peak simulations than during extreme events. In contrast, it was found that during an extreme flood peak, the use of regional forecasts and satellite data gives results that are in close agreement with results using raingauge data. Using the simulated overflow volumes recorded at the control point downstream, inundation areas were then estimated using topographic characteristics. This study is the first step in developing a future efficient early warning system and evacuation strategy.     

Longitudinal distributions of river flood power: the combined automated flood,elevation and stream power (CAFES) methodology

Stream power can be an extremely useful index of fluvial sediment transport, channel pattern, river channel erosion and riparian habitat development. However, most previous studies of downstream changes in stream power have relied on field measurements at selected cross‐sections, which are time consuming, and typically based on limited data, which cannot fully represent important spatial variations in stream power. We present here, therefore, a novel methodology we call CAFES (combined automated flood, elevation and stream power), to quantify downstream change in river flood power, based on integrating in a GIS framework Flood Estimation Handbook systems with the 5 m grid NEXTMap Britain digital elevation model derived from IFSAR (interferometric synthetic aperture radar). This provides a useful modelling platform to quantify at unprecedented resolution longitudinal distributions of flood discharge, elevation, floodplain slope and flood power at reach and basin scales. Values can be resolved to a 50 m grid. CAFES approaches have distinct advantages over current methodologies for reach‐ and basin‐scale stream power assessments and therefore for the interpretation and prediction of fluvial processes. The methodology has significant international applicability for understanding basin‐scale hydraulics, sediment transport, erosion and sedimentation processes and river basin management. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling hydrological processes influenced by soil,rock and vegetation in a small karst basin of southwest China

Hydrological processes in karst basins are controlled by permeable multimedia, consisting of soil pores, epikarst fractures, and underground conduits. Distributed modelling of hydrological dynamics in such heterogeneous hydrogeological conditions is a challenging task. Basing on the multilayer structure of the distributed hydrology‐soil‐vegetation model (DHSVM), a distributed hydrological model for a karst basin was developed by integrating mathematical routings of porous Darcy flow, fissure flow and underground channel flow. Specifically, infiltration and saturated flow movement within epikarst fractures are expressed by the ‘cubic law’ equation which is associated with fractural width, direction, and spacing. A small karst basin located in Guizhou province of southwest China was selected for this hydrological simulation. The model parameters were determined on the basis of field measurement and calibrated against the observed soil moisture contents, vegetation interception, surface runoff, and underground flow discharges from the basin outlet. The results show that due to high permeability of the epikarst zone, a significant amount of surface runoff is only generated after heavy rainfall events during the wet season. Rock exposure and the epikarst zone significantly increase flood discharge and decrease evapotranspiration (ET) loss; the peak flood discharge is directly proportional to the size of the aperture. Distribution of soil moisture content (SMC) primarily depends on topographic variations just after a heavy rainfall, while SMC and actual ET are dominated by land cover after a period of consecutive non‐rainfall days. The new model was able to capture the sharp increase and decrease of the underground streamflow hydrograph, and as such can be used to investigate hydrological effects in such rock features and land covers. Copyright © 2011 John Wiley & Sons, Ltd.     

On some aspects of peaks-over-threshold modeling of floods under nonstationarity using climate covariates

This paper discusses some aspects of flood frequency analysis using the peaks-over-threshold model with Poisson arrivals and generalized Pareto (GP) distributed peak magnitudes under nonstationarity, using climate covariates. The discussion topics were motivated by a case study on the influence of El Niño–Southern Oscillation on the flood regime in the Itajaí river basin, in Southern Brazil. The Niño3.4 (DJF) index is used as a covariate in nonstationary estimates of the Poisson and GP distributions scale parameters. Prior to the positing of parametric dependence functions, a preliminary data-driven analysis was carried out using nonparametric regression models to estimate the dependence of the parameters on the covariate. Model fits were evaluated using asymptotic likelihood ratio tests, AIC, and Q–Q plots. Results show statistically significant and complex dependence relationships with the covariate on both nonstationary parameters. The nonstationary flood hazard measure design life level (DLL) was used to compare the relative performances of stationary and nonstationary models in quantifying flood hazard over the period of records. Uncertainty analyses were carried out in every step of the application using the delta method.     

Challenges of hydrological analysis for water resource development in semi-arid mountainous regions: case study in Iran

Abstract

This study modified the BTOPMC (Block-wise TOPMODEL with the Muskingum-Cunge routing method) distributed hydrological model to make it applicable to semi-arid regions by introducing an adjustment coefficient for infiltration capacity of the soil surface, and then applied it to two catchments above the dams in the Karun River basin, located in semi-arid mountain ranges in Iran. The application results indicated that the introduced modification improved the model performance for simulating flood peaks generated by infiltration excess overland runoff at a daily time scale. The modified BTOPMC was found to fulfil the need to reproduce important signatures of basin hydrology for water resource development, such as annual runoff, seasonal runoff, low flows and flood flows. However, it was also very clear that effective model use was significantly constrained by the scarcity of ground-gauged precipitation data. Considerable efforts to improve the precipitation data acquisition should precede water resource development planning.

Editor D. Koutsoyiannis     

Modelling the hydrological response of mediterranean catchments,Prades, Catalonia. The use of distributed models as aids to hypothesis formulation

This paper attempts to extend the physical arguments underlying the distributed TOPMODEL concepts in an application to the strongly seasonal contributing area responses in two adjacent small mediterranean catchments in the Prades region of Catalonia, Spain. A perceptual model of hydrological response in these catchments is used to suggest possible modifications of the model in a hypothesis testing framework, including an attempt to modify the topographic index approach to reflect the expansion of the effective area of subsurface flow during the wetting-up sequence. It is found that slight improvements in modelling efficiency are possible but that different model parameter distributions are appropriate for different parts of the record. The model was much more successful for the catchment producing the higher runoff volumes. © 1997 John Wiley & Sons, Ltd.     

A bivariate gamma distribution for use in multivariate flood frequency analysis

A gamma distribution is one of the most frequently selected distribution types for hydrological frequency analysis. The bivariate gamma distribution with gamma marginals may be useful for analysing multivariate hydrological events. This study investigates the applicability of a bivariate gamma model with five parameters for describing the joint probability behavior of multivariate flood events. The parameters are proposed to be estimated from the marginal distributions by the method of moments. The joint distribution, the conditional distribution, and the associated return periods are derived from marginals. The usefulness of the model is demonstrated by representing the joint probabilistic behaviour between correlated flood peak and flood volume and between correlated flood volume and flood duration in the Madawask River basin in the province of Quebec, Canada. Copyright © 2001 John Wiley & Sons, Ltd.     

The bivariate lognormal distribution to model a multivariate flood episode

Complex hydrological events such as floods always appear to be multivariate events that are characterized by a few correlated variables. A complete understanding of these events needs to investigate joint probabilistic behaviours of these correlated variables. The lognormal distribution is one of frequently selected candidates for flood‐frequency analysis. The multivariate lognormal distribution will serve as an important tool for analysing a multivariate flood episode. This article presents a procedure for using the bivariate lognormal distribution to describe the joint distributions of correlated flood peaks and volumes, and correlated flood volumes and durations. Joint distributions, conditional distributions, and the associated return periods of these random variables can be readily derived from their marginal distributions. The approach is verified using observed streamflow data from the Nord river basin, located in the Province of Quebec, Canada. The theoretical distributions show a good fit to observed ones. Copyright © 2000 John Wiley & Sons, Ltd.     

变化环境下珠江流域洪水频率变化特征、成因及影响(1951-2010年)

气候变化和人类活动导致珠江流域水文变化,变化前后洪水频率分布显著不同.运用滑动秩和(Mann-Whitney U test)结合Brown-Forsythe、滑动T、有序聚类和Mann-Kendall检验法,并用累积距平曲线法获取年最大流量序列详细信息,综合确定样本最佳变化节点,并对水文变化成因做了系统分析.在此基础上,对整体序列、变化前后序列用线性矩法推求广义极值分布参数以及不同重现期设计流量.结果表明:(1)西江大部以及北江流域最佳变化节点在1991年左右;东江流域最佳变化节点与该流域内3大控制性水库建成时间基本吻合;(2)变化后,西江、北江年最大流量持续增加,洪峰强度增大,尤其是西江干流年最大流量显著增加;东江流域年最大流量显著减小,洪峰强度降低;(3)变化后,西江与北江洪水风险增加,尤其是下游珠三角地区本身受人类活动显著影响,加之西江与北江持续增加的洪水强度,珠三角地区发生洪水的强度及频次加剧,而东江洪水风险减小.此研究对于珠江流域在变化环境下的洪水风险评估与防洪抗灾具有重要意义.