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.     

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.     

Further testing of the Integrated Hydrology Model (InHM): event‐based simulations for a small rangeland catchment located near Chickasha,Oklahoma

In the paper that is the foundation for this study, VanderKwaak and Loague (2001. Water Resources Research 37 : 999–1013) reported a demonstration of a fully coupled comprehensive physics‐based hydrologic‐response model, InHM (Integrated Hydrology Model), for two rainfall‐runoff events from the small rangeland catchment known as R‐5. The InHM simulations reported herein address (in three phases) limitations in the VanderKwaak and Loague (2001. Water Resources Research 37 : 999–1013) simulations. In Phase I, a new finite‐element mesh was selected to represent R‐5. In Phase II, with the new mesh in place, evaporation was considered for the R‐5 events. In Phase III, with the new mesh in place and evaporation considered, the geology of R‐5 was approximated. Each phase, compared with the results reported by VanderKwaak and Loague (2001. Water Resources Research 37 : 999–1013), shows a change in the simulated near‐surface response. The performance of InHM for 15 R‐5 events is also reported herein. The results from two stages of model calibration are presented. The uncertainty in initial soil‐water content estimates for event‐based simulation is shown to be a major limitation for physics‐based models. The performance of InHM, relative to past event‐based simulation efforts with a quasi‐physically based rainfall‐runoff model, is better for both peak stormflow and the time to peak stormflow, but worse for stormflow depth. The InHM simulations reported here set the stage for continuous simulation of near‐surface response for the R‐5 catchment with InHM. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

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.     

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.     

A model for distributed GIUH‐based flow routing on natural and anthropogenic hillslopes

Attempts to reduce the number of parameters in distributed rainfall–runoff models have not yet resulted in a model that is accurate for both natural and anthropogenic hillslopes. We take on the challenge by proposing a distributed model for overland flow and channel flow based on a combination of a linear response time distribution and the hillslope geomorphologic instantaneous unit hydrograph (GIUH), which can be calculated with only a digital elevation model and a map with field boundaries and channel network as input. The spatial domain is subdivided into representative elementary hillslopes (REHs) for each of which we define geometric and flow velocity parameters and compute the GIUH. The catchment GIUH is given by the sum of all REH responses. While most distributed models only perform well on natural hillslopes, the advantage of our approach is that it can also be applied to modified hillslopes with for example a rectangular drainage network and terrace cultivation. Tests show that the REH‐GIUH approach performs better than classical routing functions (exponential and gamma). Simulations of four virtual hillslopes suggest that peak flow at the catchment outlet is directly related to drainage density. By combining the distributed flow routing model with a lumped‐parameter infiltration model, we were also able to demonstrate that terrace cultivation delays the response time and reduces peak flow in comparison to the same hillslope, but with a natural stream network. The REH‐GIUH approach is a first step in the process of coupling distributed hydrological models to erosion and water quality models at the REH (associated with agricultural management) and at the catchment scale (associated with the evaluation of the environmental impact of human activities). It furthermore provides a basis for the development of models for large catchments and urban or peri‐urban catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

APPLICATION OF UNIT HYDROGRAPH TECHNIQUES TO SOLUTE TRANSPORT IN CATCHMENTS

Current models of solute movement in catchments are based on rainfall–runoff models and are consequently biased towards processes which determine the magnitude and timing of water flux. It is shown here that the instantaneous unit hydrograph (IUH), or runoff response function, obtained from a hydrograph is fundamentally different from the residence time distribution which governs the response to solutes/tracers. Using hydrometric and tracer data obtained from a small (25 ha) catchment in the humid tropics a modification of the IUH technique is demonstrated which also allows approximate modelling of the tracer data. New features of the modified conceptual model are identified with known hillslope processes.     

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.     

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.     

Application of the geomorphological instantaneous unit hydrograph theory to development of forecasting models in Poland

ABSTRACT

Geomorphological instantaneous unit hydrograph (GIUH) theory has been applied for the estimation of the parameters of two conceptual models: a linear cascade model and a Laurenson-type model. Conceptual models, especially the linear cascade model, are more convenient for operational forecasting than the original GIUH model. A single linear reservoir model is suggested, with limited storage to represent the subsurface flow component. Subsurface flow is significant in Polish mountainous river catchments. Preliminary results of applying the model to operational flood forecasting are described.     

Correlation of the velocity parameter of three geomorphological instantaneous unit hydrograph models

Three geomorphological instantaneous unit hydrograph (GIUH) models are investigated. These GIUHs were derived as a function of watershed geomorphological characteristics. The geomorphological input parameters of the models were determined for 10 basins in Indiana, USA. The three GIUH models were used to simulate 187 observed runoff hydrographs of these basins. The non-measurable velocity parameters of the GIUH models were optimized and the relationships between them were investigated. The results of the investigation show that the velocity parameters of the three models are correlated. © 1998 John Wiley & Sons, Ltd.     

The influence of elevation error on the morphometrics of channel networks extracted from DEMs and the implications for hydrological modelling

Stream network morphometrics have been used frequently in environmental applications and are embedded in several hydrological models. This is because channel network geometry partly controls the runoff response of a basin. Network indices are often measured from channels that are mapped from digital elevation models (DEMs) using automated procedures. Simulations were used in this paper to study the influence of elevation error on the reliability of estimates of several common morphometrics, including stream order, the bifurcation, length, area and slope ratios, stream magnitude, network diameter, the flood magnitude and timing parameters of the geomorphological instantaneous unit hydrograph (GIUH) and the network width function. DEMs of three UK basins, ranging from high to low relief, were used for the analyses. The findings showed that moderate elevation error (RMSE of 1·8 m) can result in significant uncertainty in DEM‐mapped network morphometrics and that this uncertainty can be expressed in complex ways. For example, estimates of the bifurcation, length and area ratios and the flood magnitude and timing parameters of the GIUH each displayed multimodal frequency distributions, i.e. two or more estimated values were highly likely. Furthermore, these preferential estimates were wide ranging relative to the ranges typically observed for these indices. The wide‐ranging estimates of the two GIUH parameters represented significant uncertainty in the shape of the unit hydrograph. Stream magnitude, network diameter and the network width function were found to be highly sensitive to elevation error because of the difficulty in mapping low‐magnitude links. Uncertainties in the width function were found to increase with distance from outlet, implying that hydrological models that use network width contain greater uncertainty in the shape of the falling limb of the hydrograph. In light of these findings, care should be exercised when interpreting the results of analyses based on DEM‐mapped stream networks. Copyright © 2007 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.     

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.     

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.     

A semi‐distributed parallel‐type linear reservoir rainfall‐runoff model and its application in Taiwan

The main purpose of this paper is to introduce a semi‐distributed parallel surface rainfall‐runoff conceptual model. In this paper, a general solution of the instantaneous unit hydrograph (IUH) has been derived successfully for N linearly connected reservoirs, each having a different storage constant. The solution is a function of geomorphologic parameters, meteorologic factors and roughness coefficients. The model also takes into account the hydrologic response which is influenced by outflow downstream of a reservoir. For calibration, the shuffled complex evolution (SCE) algorithm is used to search for the global optimal parameters of the model. Because of the parallel structure, the mean roughness parameter of the channel becomes a “conceptual” parameter without a real physical meaning. To evaluate the adaptability of the model adopted, three watersheds around the city of Taipei in Taiwan were chosen to test the effectiveness of the model. The study provides an appropriate rainfall‐runoff model for planning flood mitigation in Taiwan. Copyright © 1999 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.     

The application of electrical conductivity as a tracer for hydrograph separation in urban catchments

Two‐component hydrograph separation was performed on 19 low‐to‐moderate intensity rainfall events in a 4·1‐km² urban watershed to infer the relative and absolute contribution of surface runoff (e.g. new water) to stormflow generation between 2001 and 2003. The electrical conductivity (EC) of water was used as a continuous and inexpensive tracer, with order of magnitude differences in precipitation (12–46 µS/cm) and pre‐event streamwater EC values (520–1297 µS/cm). While new water accounted for most of the increased discharge during storms (61–117%), the contribution of new water to total discharge during events was typically lower (18–78%) and negatively correlated with antecedent stream discharge (r² = 0·55, p < 0·01). The amount of new water was positively correlated with total rainfall (r² = 0·77), but hydrograph separation results suggest that less than half (9–46%) of the total rainfall on impervious surfaces is rapidly routed to the stream channel as new water. Comparison of hydrograph separation results using non‐conservative tracers (EC and Si) and a conservative isotopic tracer (δD) for two events showed similar results and highlighted the potential application of EC as an inexpensive, high frequency tracer for hydrograph separation studies in urban catchments. The use of a simple tracer‐based approach may help hydrologists and watershed managers to better understand impervious surface runoff, stormflow generation and non‐point‐source pollutant loading to urban streams. Copyright © 2007 John Wiley & Sons, Ltd.