Runoff estimation for an ungauged catchment using geomorphological instantaneous unit hydrograph (GIUH) models

A geomorphological instantaneous unit hydrograph (GIUH) is derived from the geomorphological characteristics of a catchment and it is related to the parameters of the Clark instantaneous unit hydrograph (IUH) model as well as the Nash IUH model for deriving its complete shape. The developed GIUH based Clark and Nash models are applied for simulation of the direct surface run‐off (DSRO) hydrographs for ten rainfall‐runoff events of the Ajay catchment up to the Sarath gauging site of eastern India. The geomorphological characteristics of the Ajay catchment are evaluated using the GIS package, Integrated Land and Water Information System (ILWIS). The performances of the GIUH based Clark and Nash models in simulating the DSRO hydrographs are compared with the Clark IUH model option of HEC‐1 package and the Nash IUH model, using some commonly used objective functions. The DSRO hydrographs are computed with reasonable accuracy by the GIUH based Clark and Nash models, which simulate the DSRO hydrographs of the catchment considering it to be ungauged. Inter comparison of the performances of the GIUH based Clark and Nash models shows that the DSRO hydrographs are estimated with comparable accuracy by both the models. Copyright © 2007 John Wiley & Sons, Ltd.     

Comparative analysis of a geomorphology‐based instantaneous unit hydrograph in small mountainous watersheds

A conceptual insytnataneous unit hydrograph (IUH) based on geomorphologival association of linear reservoirs (GR) previously developed by the authors has been compared with other IUH models: a distributed GR variation (GR(v)), the Nash IUH, the Chutha and Dooge IUH, and the Troutman and Karlinger IUH for the analysis of direct runoff hydrographs recorded in three experimental watershed of the north of Spain. The comparison was made through a calibration‐validation process in which a leave‐one‐out cross‐validation method was applied. The results indicate the satisfactory performance of all the models, with the advantage for the GR model of the dependence on only one parameter, which can be identified from the watershed and event characteristics. This property makes its use easier than that of other models. Copyright © 2011 John Wiley & Sons, Ltd.     

Application of the Kalman filter to the Nash model

The Nash model was used for application of the Kalman filter. The state vector of the rainfall–runoff system was constituted by the IUH (instantaneous unit hydrograph) estimated by the Nash model and the runoff estimated by the Nash model using the Kalman filter. The initial values of the state vector were assumed as the average of 10% of the IUH peak values and the initial runoff estimated from the average IUH. The Nash model using the Kalman filter with a recursive algorithm accurately predicted runoff from a basin in Korea. The filter allowed the IUH to vary in time, increased the accuracy of the Nash model and reduced physical uncertainty of the rainfall–runoff process in the river basin. © 1998 John Wiley & Sons, Ltd.     

Estimation of Nash's IUH parameters using stochastic search algorithms

In rainfall–runoff studies, it is often necessary to change the duration of a given unit hydrograph. Nash's Instantaneous Unit Hydrograph (IUH) is an ideal method that eliminates the hydrograph duration. This paper presents the results of the application of search algorithms, namely a genetic algorithm and hill climbing, to develop the IUH that minimizes the error between the observed and generated hydrographs. Also the performance of these methods has been compared with that of the classical method used for estimation of IUH, namely the method of moments. The genetic algorithm is a popular search procedure for function optimization that applies the mechanics of natural genetics and natural selection to explore a given search space. Hill climbing is an optimization technique that belongs to the family of local search and algorithms can be used to solve problems that have many solutions, with some solutions better than others. The results obtained from both the genetic algorithm and hill climbing algorithm for estimation of Nash's IUH parameters were compared with the results obtained by the method of moments for storms from two river basins that are located in different climatic regions. It was found that both the genetic algorithm and hill climbing provided improved and consistent results. Copyright © 2008 John Wiley & Sons, Ltd.     

A self‐organization algorithm for real‐time flood forecast

The group method of data handling (GMDH) algorithm presented by A. C. Ivakhnenko and colleagues is an heuristic self‐organization method. It establishes the input–output relationship of a complex system using a multilayered perception‐type structure that is similar to a feed‐forward multilayer neural network. This study provides a step towards understanding and evaluating a role for GMDH in the investigation of the complex rainfall–runoff processes in a heterogeneous watershed in Taiwan. Two versions of the revised GMDH model are implemented: a stepwise regression procedure and a recursive formula. Eleven typhoon events in the Shen‐cei Creek watershed, Taiwan, are used to build the model and verify its usefulness. The prediction results of the revised GMDH models and the instantaneous unit hydrograph (IUH) model are compared. Based on the criteria of forecasting precision and the rate and time of peak error, a much better performance is obtained with the revised GMDH models. Copyright © 1999 John Wiley & Sons, Ltd.     

Exponential data fitting applied to infiltration, hydrograph separation, and variogram fitting

Most lumped rainfall-runoff models separate the interflow and groundwater components from the measured runoff hydrograph in an attempt to model these as hydrologic reservoir units. Similarly, rainfall losses due to infiltration as well as other abstractions are separated from the measured rainfall hyetograph, which are then used as inputs to the various hydrologic reservoir units. This data pre-processing is necessary in order to use the linear unit hydrograph theory, as well as for maintaining a hydrologic budget between the surface and subsurface flow processes. Since infiltration determines the shape of the runoff hydrograph, it must be estimated as accurately as possible. When measured infiltration data is available, Horton’s exponential infiltration model is preferable due to its simplicity. However, estimating the parameters from Horton’s model constitutes a nonlinear least squares fitting problem. Hence, an iterative procedure that requires initialization is subject to convergence. In a similar context, the separation of direct runoff, interflow, and baseflow from the total hydrograph is typically done in an ad hoc manner. However, many practitioners use exponential models in a rather “layer peeling” fashion to perform this separation. In essence, this also constitutes an exponential data fitting problem. Likewise, certain variogram functions can be fitted using exponential data fitting techniques. In this paper we show that fitting a Hortonian model to experimental data, as well as performing hydrograph separation, and total hydrograph and variogram fitting can all be formulated as a system identification problem using Hankel-based realization algorithms. The main advantage is that the parameters can be estimated in a noniterative fashion, using robust numerical linear algebra techniques. As such, the system identification algorithms overcome the problem of convergence inherent in iterative techniques. In addition, the algorithms are robust to noise in the data since they optimally separate the signal and noise subspaces from the observed noisy data. The algorithms are tested with real data from field experiments performed in Surinam, as well as with real hydrograph data from a watershed in Louisiana. The system identification techniques presented herein can also be used with any other type of exponential data such as exponential decays from nuclear experiments, tracer studies, and compartmental analysis studies.     

Analysis of geomorphologic and hydrological characteristics in watershed saturated areas using topographic‐index threshold and geomorphology‐based runoff model

The objective of this paper is to investigate the variation of geomorphology and runoff characteristics in saturated areas under different partial contributing area (PCA) conditions. Geomorphologic information and hydrologic records from two mid‐size watersheds in northern Taiwan were selected for analysis. The PCA ratio in the watershed during a storm was assumed equal to the ratio of the surface‐flow volume to the direct runoff volume from measured hydrologic data. The extents of PCA regions were then determined by using a topographic‐index threshold. Consequently, the geomorphologic factors in saturated and unsaturated areas could be calculated using a digital elevation model, and these factors could then be linked to a geomorphology‐based IUH model for runoff simulation, which can consider both the surface‐ and subsurface‐flow processes in saturated and unsaturated areas, respectively. The results show that geomorphologic characteristics in the saturated areas vary significantly with different PCA ratios especially for higher order streams. A large PCA ratio results in a sharp hydrograph because the quick surface flow dominates the runoff process, whereas the hydrologic response in a low PCA case is dominated by the delayed subsurface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

Use of Nash's IUH and DEMs to identify the parameters of an unequal‐reservoir cascade IUH model

Under the assumption that hydrograph generation was affected by n linear reservoirs with the same value of storage coefficient k, Nash proposed the formulation of the Instantaneous Unit Hydrograph (IUH), which has been widely used in rainfall–runoff simulation and flood forecasting. However, the assumption of the parameter k having the same value in all reservoirs is obviously unphysical as it results in the estimated value of n not being integral. In this study, for parameter n integral, the different k value for each reservoir was derived using the Laplace transform and developing a general rule for the equation of the IUH of any order. The relationship between parameter k and the slope of the river channel estimated using digital elevation model (DEM) data is established, the parameter estimation procedures are given. As in most unit hydrograph studies, only isolated storm events are considered here. Seventeen flood events in three catchments were selected for the case studies. Application results show that the proposed method is slightly better than Nash's IUH with higher model efficiency and smaller absolute relative errors. This work provides a new methodology for the formulation of the IUH. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of a regionalized Clark IUH model with limited hydrologic data availability

To aid prediction of the flow hydrograph in a basin with limited data, a practical approach to determining a regionalized Clark instantaneous unit hydrograph (IUH) model is presented. The proposed model is described in terms of the synthetic time–area concentration curve, the concentration time, and a special regional similarity value that is valid in the whole basin. The latter was estimated from a Monte Carlo testing procedure based on the normal probability distribution of transformed regional similarity values composed of the time of concentration and the storage coefficient in gauged basins. The time–area concentration curve and the concentration time were calculated from a rational equation as in conventional methods. The method of transformation adopted was the Box–Cox power transformation, which is known to make non‐normal values resemble normal data. By introducing the regional similarity value into a Clark IUH, a statistically best estimate of IUH for given data conditions and its quantified degree of uncertainty were realized. The Wi River basin in Korea was used to test the applicability of the regionalized Clark IUH. The performance of the suggested methodology was evaluated by assuming an ungauged sub‐basin at the site. The results showed that the IUH model developed in this work was an effective tool, predicting a reliable hydrograph within the study area even though only limited data were available. Copyright © 2008 John Wiley & Sons, Ltd.     

Book Review

Abstract

The instantaneous unit hydrograph (IUH) of a watershed is the result of one instantaneous unit of rainfall excess distributed uniformly over the watershed. Although the geomorphological characteristics of the basin remain relatively constant, the variable characteristics of storms cause variations in the shape of the resulting hydrographs. It is, therefore, inadequate to use one typical IUH to represent the hydrological response generated from any specific storm. In this study, a variable IUH was derived that directly reflects the time-varying rainfall intensity during storms. The rainfall intensity used to generate the variable IUH at time t is the mean rainfall intensity occurring from the time t—T _c to t in which T _c is the watershed time of concentration. Hydrological records from three watersheds in Taiwan were used to demonstrate the applicability of the proposed model. The results show that better simulations can be obtained by using the proposed model than by using the conventional unit hydrograph method, especially for concentrated rainstorm cases.     

A technique for deriving an average event unit hydrograph from streamflow—only data for ephemeral quick-flow-dominant catchments

A technique is presented for empirical and parametric estimation of an average event unit hydrograph response curve solely from measured streamflow data for use in rainfall-runoff models with a focus on modelling flow for natural resource management. As the technique does not require rainfall data, the unit hydrograph can be derived solely from streamflow data at a temporal resolution appropriate for the response of the catchment. The response curves derived for a number of stations show a power law relationship of decay in flow after peak that can be described using a three-parameter function. The approach is best suited to ephemeral streamflow regimes dominated by surface and near-surface runoff (quick-flow component) where flow events are well-separated and largely independent of each other. Analysis of the derived unit hydrograph for 28 quick-flow-dominated catchments illustrates the range of parameter values obtained in fitting the power law.     

The meaning of oscillations in unit hydrograph S-curves

ABSTRACT

The Duhamel superposition integral is used to obtain some exact solutions for unit hydrograph applications. These equations and numerical examples are used to show that oscillations will occur in an S-curve when the time step is less than the excess rainfall duration if the measured hydrograph differs from a hydrograph that would be obtained by solving a linear differential equation with time-independent coefficients. The implications of this result with regard to the calculation of the instantaneous unit hydrograph (IUH) are discussed.     

A practical estimation of Clark IUH parameters using root selection and linear programming

The values of the parameters of the Clark instantaneous unit hydrograph (IUH) are often relying on the subjective decision of the researcher, which leads to large variations of their values. Therefore, an objective method minimizing the subjective judgement in the IUH modelling procedure while providing a reduced range of acceptable values is proposed. The proposed method uses a basin average IUH to mitigate the robustness problem of the Clark IUH parameters. Using linear system theory, the z‐transform is applied to the average IUH and then the IUH polynomial is factored into the recession and time‐area curve (TAC) components based on a convolution relation between the Clark IUH parameters. During this calculation, the root selection method was adopted to verify the storage coefficient R from the recession component and a linear programming technique was applied for determining the TAC for the basin of interest. The Wi River basin was used to test the applicability of the proposed method. The results showed that the components of a single reservoir and the TAC for Clark IUH were separated effectively, and acceptable values for the parameters were obtained. Copyright © 2011 John Wiley & Sons, Ltd.     

Reservoir attenuation of floods from ungauged basins

Abstract

In a typical reservoir routing problem, the givens are the inflow hydrograph and reservoir characteristic functions. Flood attenuation investigations can be easily accomplished using a hydrological or hydraulic routing of the inflow hydrograph to obtain the reservoir outflow hydrograph, unless the inflow hydrograph is unavailable. Although attempts for runoff simulation have been made in ungauged basins, there is only a limited degree of success in special cases. Those approaches are, in general, not suitable for basins with a reservoir. The objective of this study is to propose a procedure for flood attenuation estimation in ungauged reservoir basins. In this study, a kinematic-wave based geomorphic IUH model was adopted. The reservoir inflow hydrograph was generated through convolution integration using the rainfall excess and basin geomorphic information. Consequently, a fourth-order Runge-Kutta method was used to route the inflow hydrograph to obtain the reservoir outflow hydrograph without the aid of recorded flow data. Flood attenuation was estimated through the analysis of the inflow and outflow hydrographs of the reservoir. An ungauged reservoir basin in southern Taiwan is presented as an example to show the applicability of the proposed analytical procedure. The analytical results provide valuable information for downstream flood control work for different return periods.     

Two‐dimensional continuous simulation of spatiotemporally varied hydrological processes using the time‐area method

Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrological effects of within-catchment heterogeneity of drainage density

Local drainage density (dd) has been traditionally defined as the inverse of twice the distance one has to walk before encountering a channel. This formalization easily allows to derive raster-based maps of dd extracted straight off from digital elevation model data. Maps of local dd, which are continuous in space, are able to reveal the appearance of strong heterogeneities in the geological and geomorphological properties of natural landscapes across different scales. In this work we employ the information provided by these spatial maps to study the potential effects of the within-catchment variability of dd on the hydrologic response. A simple power law relationship between runoff yield at the local scale and the value of dd has been adopted; the hypothesis is supported by a large number of past empirical observations and modeling. The novel framework proposed (ddRWF) embeds this spatially variable runoff weight in the well-known Rescaled Width Function (RWF) framework, based on the more general geomorphological theory of the hydrologic response. The model is applied to four sub-basins in the Cascade Range Region (Oregon, USA) where strong contrasts in dissection patterns due the underlain geology have been broadly addressed in previous literature. The ddRWF approach is compared with the classic RWF in terms of shape, moments and peak of the simulated hydrograph response. Results hint that the variability of runoff yield due to the heterogeneity of dd (i.e. of hillslope lengths) determines a more rapid concentration of runoff, which implies shorter lag times, larger skewness and higher peak floods, especially in the case hillslope velocity is much smaller than channel velocity. The potential of the proposed framework relies on accounting for spatially variable losses related to geomorphologic heterogeneity in lumped rainfall–runoff models, still keeping the simple and robust structure of the IUH approach.     

Analytical model for a geomorphological instantaneous unit hydrograph

The rainfall‐runoff modelling of a river basin can be divided into two processes: the production function and the transfer function. The production function determines the proportion of gross rainfall actually involved in the runoff. The transfer function spreads the net rainfall over time and space in the river basin. Such a transfer function can be modelled using the approach of the geomorphological instantaneous unit hydrograph (GIUH). The effectiveness of geomorphological models is actually revealed in rainfall‐runoff modelling, where hydrologic data are desperately lacking, just as in ungauged basins. These models make it possible to forecast the hydrograph shape and runoff variation versus time at the basin outlet. This article is an introduction to a new GIUH model that proves to be simple and analytical. Its geomorphological parameters are easily available on a map or from a digital elevation model. This model is based on general hypotheses on symmetry that provide it with multiscale versatile characteristics. After having validated the model in river basins of very different nature and size, we present an application of this model for rainfall‐runoff modelling. Since parameters are determined relying on real geomorphological data, no calibration is necessary, and it is then possible to carry out rainfall‐runoff simulations in ungauged river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

An investigation of lag times for rainfall–runoff–sediment yield events in small river basins / Une analyse des temps de réponse d'événements pluie–débit–transport solide au sein de petits bassins versants

Abstract

River basin lag time (LAG), defined as the elapsed time between the occurrence of the centroids of the effective rainfall intensity pattern and the storm runoff hydrograph, is an important factor in determining the time to peak and the peak value of the instantaneous unit hydrograph, IUH. In the procedure of predicting a sedimentgraph (suspended sediment load as a function of time), the equivalent parameter is the lag time for the sedimentgraph (LAG_s ), which is defined as the elapsed time between the occurrence of the centroids of sediment production during a storm event and the observed sedimentgraph at the gauging station. Results of analyses of rainfall, runoff and suspended sediment concentration event data collected from five small Carpathian basins in Poland and from a 2.31-ha agricultural basin, in central Illinois, USA have shown that LAG_s was, in the majority of cases, smaller than LAG, and that a significant linear relationship exists between LAG_s and LAG.     

Towards more systematic perceptual model development: a case study using 3 Luxembourgish catchments

The synthesis of experimental understanding of catchment behaviour and its translation into qualitative perceptual models is an important objective of hydrological sciences. We explore this challenge by examining the cumulative understanding of the hydrology of three experimental catchments and how it evolves through the application of different investigation techniques. The case study considers the Huewelerbach, Weierbach and Wollefsbach headwater catchments of the Attert basin in Luxembourg. Subsurface investigations including bore holes and pits, analysis of soil samples and Electrical Resistivity Tomography measurements are presented and discussed. Streamflow and tracer data are used to gain further insights into the streamflow dynamics of the catchments, using end‐member mixing analysis and hydrograph separation based on dissolved silica and electrical conductivity. We show that the streamflow generating processes in all three catchments are controlled primarily by the subsolum and underlying bedrock. In the Huewelerbach, the permeable sandstone formation supports a stable groundwater component with little seasonality, which reaches the stream through a series of sources at the contact zone with the impermeable marls formation. In the Weierbach, the schist formation is relatively impermeable and supports a ‘fill and spill’‐type of flow mechanism; during wet conditions, it produces a delayed response dominated by pre‐event water. In the Wollefsbach, the impermeable marls formation is responsible for a saturation‐excess runoff generating process, producing a fast and highly seasonal response dominated by event water. The distinct streamflow generating processes of the three catchments are represented qualitatively using perceptual models. The perceptual models are in turn translated into quantitative conceptual models, which simulate the hydrological processes using networks of connected reservoirs and transfer functions. More generally, the paper illustrates the evolution of perceptual models based on experimental fieldwork data, the translation of perceptual models into conceptual models and the value of different types of data for processes understanding and model representation. Copyright © 2014 John Wiley & Sons, Ltd.     

Does the incorporation of process conceptualization and tracer data improve the structure and performance of a simple rainfall‐runoff model in a Scottish mesoscale catchment?

A geomorphological instantaneous unit hydrograph (GIUH) rainfall‐runoff model was applied in a 31 km² montane catchment in Scotland. Modelling was based on flow path length distributions derived from a digital terrain model (DTM). The model was applied in two ways; a single landscape unit response based on the DTM alone, and a two‐landscape unit response, which incorporated the distribution of saturated areas derived from field‐validated geographic information system (GIS) analysis based on a DTM and soil maps. This was to test the hypothesis that incorporation of process‐information would enhance the model performance. The model was applied with limited multiple event calibration to produce parameter sets which could be applied to a spectrum of events with contrasting characteristics and antecedent conditions. Gran alkalinity was used as a tracer to provide an additional objective measure for assessing model performance. The models captured the hydrological response dynamics of the catchment reasonably well. In general, the single landscape unit approach produced the best individual model performance statistics, though the two‐landscape unit approach provided a range of models, which bracketed the storm hydrograph response more realistically. There was a tendency to over‐predict the rising limb of the hydrograph, underestimate large storm event peaks and anticipate the hydrograph recession too rapidly. Most of these limitations could be explained by the simplistic assumptions embedded within the GIUH approach. The modelling also gave feasible predictions of stream water chemistry, though these could not be used as a basis for model rejection. Nevertheless, the study suggested that the approach has potential for prediction of hydrological response in ungauged montane headwater basins. Copyright © 2007 John Wiley & Sons, Ltd.