A spatially distributed Clark’s unit hydrograph based hybrid hydrologic model (Distributed-Clark)

ABSTRACT

A hybrid hydrologic model (Distributed-Clark), which is a lumped conceptual and distributed feature model, was developed based on the combined concept of Clark’s unit hydrograph and its spatial decomposition methods, incorporating refined spatially variable flow dynamics to implement hydrological simulation for spatially distributed rainfall–runoff flow. In Distributed-Clark, the Soil Conservation Service (SCS) curve number method is utilized to estimate spatially distributed runoff depth and a set of separated unit hydrographs is used for runoff routing to obtain a direct runoff flow hydrograph. Case studies (four watersheds in the central part of the USA) using spatially distributed (Thiessen polygon-based) rainfall data of storm events were used to evaluate the model performance. Results demonstrate relatively good fit to observed streamflow, with a Nash-Sutcliffe efficiency (E_NS) of 0.84 and coefficient of determination (R²) of 0.86, as well as a better fit in comparison with outputs of spatially averaged rainfall data simulations for two models including HEC-HMS.     

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.     

On-line estimation of unit hydrographs using the wavelet-based LMS algorithm / Estimation en ligne des hydrogrammes unitaires grâce à l'algorithme des moindres carrés moyens à base d'ondelettes

Abstract

This study applies the discrete wavelet transform (DWT) to decompose the unit hydrograph, thereby generating parsimonious reparameterizations of the unit hydrograph. A model compression method is then employed to significantly compress the unit hydrograph requiring that fewer coefficients be estimated. Moreover, a wavelet-based linearly constrained least mean squares (WLCLMS) algorithm is also used to estimate on-line the wavelet coefficients of the unit hydrograph. The updated wavelet coefficients of the unit hydrograph, convoluted with effective rainfall input in the wavelet domain, allow for accurate prediction of one-step-ahead runoff in the time domain. The proposed approach allows the unit hydrographs to vary in time and accurately predicts runoff from a basin in Taiwan, thus making it highly promising for flood forecasting.     

Using recurrent neural networks to reconstruct rainfall‐runoff processes

Many novel techniques for reconstructing rainfall‐runoff processes require hydrometeorologic and geomorphologic information for modelling. However, certain information is not always measurable. In this paper, we employ a special recurrent neural network to reconstruct the rainfall‐runoff process by using collected rainfall data. In addition, we propose an indirect system identification to overcome the drawback of a traditional, time‐consuming trial‐and‐error search. The indirect system identification is an efficient method to recognize the structure of a recurrent neural network. The unit hydrograph can be derived directly from the weights of the network due to a state‐space form embedded in the recurrent neural network. This improves the link between the weights of the network and the physical concepts that most neural networks fail to connect. The case studies of 41 events from 1966 to 1997 have been implemented in Taiwan's Wu‐Tu watershed, where the runoff path‐lines are short and steep. Two recurrent neural networks and one state‐space model are utilized to simulate the rainfall‐runoff processes for comparison. The results are validated by four criteria: coefficient of efficiency; peak discharge error; time to peak arrival error; total discharge volume error. The resulting data from the recurrent neural network reveal that the neural network proposed herein is appropriate for hydrological systems. Copyright © 2005 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 design storm concept in flood control design and planning

The design storm approach, where the subject criterion variable is evaluated by using a synthetic storm pattern composed of identical return frequencies of storm pattern input, is shown to be an effective approximation to a considerably more complex probabilistic model. The single area unit hydrograph technique is shown to be an accurate mathematical model of a highly discretized catchment with linear routing for channel flow approximation, and effective rainfalls in subareas which are linear with respect to effective rainfall output for a selected “loss” function. The use of a simple “loss” function which directly equates to the distribution of rainfall depth-duration statistics (such as a constant fraction of rainfall, or a ?-index model) is shown to allow the pooling of data and thereby provide a higher level of statistical significance (in estimating T-year outputs for a hydrologic criterion variable) than use of an arbitrary “loss” function. The above design storm unit hydrograph approach is shown to provide the T-year estimate of a criterion variable when using rainfall data to estimate runoff.     

A new baseflow separation method based on analytical solutions of the Horton infiltration capacity curve

The separation of baseflow is an important issue in hydrology. The objective of this paper is to develop a new baseflow separation method based on the Horton infiltration capacity curve. For this purpose, the mathematical expressions of three parameters for the Horton infiltration capacity curve were derived in terms of rainfall and runoff data, and the lag time of the effective rainfall in the unsaturated zone and the groundwater flow routing equation are also presented. With these equations, the baseflow hydrographs at the outlet of the basin can be separated. The flow chart of the proposed method for baseflow hydrograph separation is given. Three recent baseflow separation methods, i.e. digital filter, non‐linear reservoir and the Boussinesq equation, were chosen as parallel schemes to compare with the proposed method. Rainfall‐runoff data from four watersheds located in different climatic regions in China were selected and used as case studies. Test and application results indicate that the proposed baseflow hydrograph separation method is in accordance with the hydrological physical process. The proposed method is comparable with current available methods and reduces some of the subjective aspects for the rising limb of the baseflow hydrograph, and it is useful for unit hydrograph analysis and for the study of the rainfall‐runoff relationship. Copyright © 2006 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.     

GEOMORPHOLOGICAL TRANSFER FUNCTION CALCULATED FROM DIGITAL ELEVATION MODELS FOR DISTRIBUTED HYDROLOGICAL MODELLING

Recently, several attempts have been made to relate the hydrological response of a catchment to its morphological and topographical features using different hypotheses to model the effect of the drainage network. Several transfer functions were developed and some of these are based on the theory of a linear model, the geomorphological unit hydrograph. The aim of this paper is to present a methodology to automatically identify the transfer function, using digital elevation models for applications in distributed hydrological modelling. The transfer function proposed herein is based on the Hayami approximation solution of the diffusive wave equation especially adapted for the routing hydrograph through a channel network. The Gardon d’Anduze basin, southern France, was retained for applications. Digital elevation models were used to extract the channel network and divide the basin into subcatchments. Each subcatchment produces, at its own outlet, an impulse response which is routed to the outlet of the whole catchment using the diffusive wave model described by two parameters: celerity and diffusivity functions of geometrical characteristics of the channel network. Firstly, a geomorphological unit hydrograph obtained by routing a homogeneous effective rainfall was compared with the unit hydrograph identified by a lumped model scheme, then the distributed model was applied to take into account the spatial variability of effective rainfall in the catchment. Results show that this new method seems to be adapted for distributed hydrological modelling; it enables identification of a transfer function response for each hydrological unit, here subcatchments, and then simulation of the contribution of each unit to the hydrograph at the outlet. © 1997 by John Wiley & Sons, Ltd.     

Laboratory examination of linearity in rainfall–runoff relationships

ABSTRACT

Certain rainfall–runoff models, e.g. the unit hydrograph, assume linear relationships between the variables. These are proportionality of runoff discharges to (net) rainfall depth and linear summations of discharges resulting from (net) rainfalls during different time intervals or over different sectors of a watershed. This study examines the validity of these assumptions by use of an extensive two-dimensional laboratory experimentation. The results indicate that proportionality would be found under high rainfall intensity through a long duration. Spatial summations would more likely yield correct discharges in cases where rainfall duration is equal to, or is longer than, the time of concentration. Temporal summations may yield correct discharges when rainfall duration is longer than one half of the time of concentration. Here, the time of concentration is determined at the beginning of gradual approach of the discharge towards the equilibrium state.     

Rainfall—runoff response,event-based runoff coefficients and hydrograph separation

Abstract

Event-based runoff coefficients can provide information on watershed response. They are useful for catchment comparison to understand how different landscapes “filter” rainfall into event-based runoff and to explain the observed differences with catchment characteristics and related runoff mechanisms. However, the big drawback of this important parameter is the lack of a standard hydrograph separation method preceding its calculation. Event-based runoff coefficients determined with four well-established separation methods, as well as a newly developed separation method, are compared and are shown to differ considerably. This signifies that runoff coefficients reported in the literature often convey less information than required to allow for catchment classification. The new separation technique (constant-k method) is based on the theory of linear storage. Its advantages are that it is theoretically based in determining the end point of an event and that it can also be applied to events with multiple peaks. Furthermore, it is shown that event-based runoff coefficients in combination with simple statistical models improve our understanding of rainfall—runoff response of catchments with sparse data.     

Effect of the duration and direction of storm movement on planar flow with full and partial areal coverage

Using kinematic wave equations, analytical solutions are derived for flow resulting from storms moving either up or down the plane and covering it fully or partially. By comparing the flow resulting from a moving storm with that from a stationary storm of the same duration and areal coverage, the influence of storm duration, direction and areal coverage is investigated. It is found that the direction, duration and areal coverage of storm movement have a pronounced effect on the discharge hydrograph. The runoff hydrographs resulting from storms moving downstream are quite different from those from storms moving upstream. Likewise, the areal coverage of the storm has a pronounced effect on the runoff hydrograph. Copyright © 2002 John Wiley & Sons, Ltd.     

FLOOD MODELLING WITH GIS-DERIVED DISTRIBUTED UNIT HYDROGRAPHS

The concept of a spatially distributed unit hydrograph is based on the fact that the unit hydrograph can be derived from the time–area curve of a watershed by the S-curve method. The time–area diagram is a graph of cumulative drainage area contributing to discharge at the watershed outlet within a specified time of travel. Accurate determination of the time–area diagram is made possible by using a GIS. The GIS is used to describe the connectivity of the links in the watershed flow network and to calculate distances and travel times to the watershed outlet for various points within the watershed. Overland flow travel times are calculated by the kinematic wave equation for time to equilibrium; channel flow times are based on the Manning and continuity equations. To account for channel storage, travel times for channel reaches are increased by a percentage depending on the channel reach length and geometry. With GIS capability for rainfall mapping, the assumption of a uniform spatial rainfall distribution is no longer necessary; hence the term, spatially distributed unit hydrograph. An example of the application for the Waiparous Creek in the Alberta Foothills is given. IDRISI is used to develop a simple digital elevation model of the 229 km² watershed, using 1 km × 1 km grid cells. A grid of flow directions is developed and used to create an equivalent channel network. Excess rainfall for each 1 km × 1 km cell is individually computed by the Soil Conservation Service (SCS) runoff curve method and routed through the equivalent channel network to obtain the time–area curve. The derived unit hydrograph gave excellent results in simulating an observed flood hydrograph. The distributed unit hydrograph is no longer a lumped model, since it accounts for internal distribution of rainfall and runoff. It is derived for a watershed without the need for observed rainfall and discharge data, because it is essentially a geomorphoclimatic approach. As such, it allows the derivation of watershed responses (hydrographs) to inputs of various magnitudes, thus eliminating the assumption of proportionality of input and output if needed. The superposition of outputs is retained in simulating flood hydrographs by convolution, since it has been shown that some non-linear systems satisfy the principle of superposition. The distributed unit hydrograph appears to be a very promising rainfall runoff model based on GIS technology.     

Hydrograph separation in a mountainous catchment — combining hydrochemical and isotopic tracers

Runoff components of the Zastler catchment (18\4 km², southern Black Forest, Germany) were analysed with hydrograph separations using stable oxygen isotopes and dissolved silica. It was shown that event water and components with low silica contributed only small amounts to total runoff. In addition, comparison of the two‐component hydrograph separations showed that the low‐silica components are generated by both event water and pre‐event water fractions, depending on the state of the system. A modified three‐component hydrograph separation method was introduced using dissolved silica and ¹⁸O. During storm events an interaction of three runoff components having distinct silica concentrations could be shown. Based on the geological and geomorphological genesis of the study site, it was appropriate to assign (i) the low silica component to the riparian zones and impermeable areas, (ii) the medium silica component to the periglacial debris cover and (iii) the high silica component to the crystalline detritus and crystalline hard rock. Exact quantification of the runoff components remained difficult. However, runoff components with medium silica concentrations reacted very sensitively and intensely. The contribution of this component to total runoff is comparatively large. This shows the important role of the periglacial debris to runoff generation of the study site and emphasizes the importance of runoff generation processes occurring in this reservoir. Copyright © 2000 John Wiley & Sons, Ltd.     

A stochastic integral equation analog of rainfall-runoff processes for evaluating modeling uncertainty

In this paper a very general rainfall-runoff model structure (described below) is shown to reduce to a unit hydrograph model structure. For the general model, a multi-linear unit hydrograph approach is used to develop subarea runoff, and is coupled to a multi-linear channel flow routing method to develop a link-node rainfall-runoff model network. The spatial and temporal rainfall distribution over the catchment is probabilistically related to a known rainfall data source located in the catchment in order to account for the stochastic nature of rainfall with respect to the rain gauge measured data. The resulting link node model structure is a series of stochastic integral equations, one equation for each subarea. A cumulative stochastic integral equation is developed as a sum of the above series, and includes the complete spatial and temporal variabilities of the rainfall over the catchment. The resulting stochastic integral equation is seen to be an extension of the well-known single area unit hydrograph method, except that the model output of a runoff hydrograph is a distribution of outcomes (or realizations) when applied to problems involving prediction of storm runoff; that is, the model output is a set of probable runoff hydrographs, each outcome being the results of calibration to a known storm event.     

A stochastic integral equation analog of rainfall-runoff processes for evaluating modeling uncertainty

In this paper a very general rainfall-runoff model structure (described below) is shown to reduce to a unit hydrograph model structure. For the general model, a multi-linear unit hydrograph approach is used to develop subarea runoff, and is coupled to a multi-linear channel flow routing method to develop a link-node rainfall-runoff model network. The spatial and temporal rainfall distribution over the catchment is probabilistically related to a known rainfall data source located in the catchment in order to account for the stochastic nature of rainfall with respect to the rain gauge measured data. The resulting link node model structure is a series of stochastic integral equations, one equation for each subarea. A cumulative stochastic integral equation is developed as a sum of the above series, and includes the complete spatial and temporal variabilities of the rainfall over the catchment. The resulting stochastic integral equation is seen to be an extension of the well-known single area unit hydrograph method, except that the model output of a runoff hydrograph is a distribution of outcomes (or realizations) when applied to problems involving prediction of storm runoff; that is, the model output is a set of probable runoff hydrographs, each outcome being the results of calibration to a known storm event.     

Relating hydrograph components to rainfall and streamflow: a case study from northern Taiwan

Abstract

This study investigates the characteristics of hydrograph components from a watershed in Taiwan. Hydrograph components were modelled by using a model of three serial reservoirs with one parallel reservoir. Mean rainfall was calculated by using the block kriging method. The model parameters for 38 events were calibrated by using the shuffled complex evolution optimization algorithm. The model verification was made using 18 events. Based on the study results, the following findings were obtained: (1) for single-peak events, times to peak of hydrograph components are an increasing power function of the peak time of rainfall; (2) peak discharges of hydrograph components are linearly proportional to that of total runoff, and the ratios of quick and slow runoff are approximately 83% and 17% of total runoff, respectively; and (3) the total volume of quick runoff component is 52% of total runoff and that of slow runoff is 27%.

Editor D. Koutsoyiannis

Citation Li, Y.-J., Cheng, S.-J. Pao, T.-L. and Bi, Y.-J., 2012. Relating hydrograph components to rainfall and streamflow: a case study from northern Taiwan. Hydrological Sciences Journal, 57 (5), 861–877.     

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.     

Runoff and suspended sediment yields from an unpaved road segment,St John,US Virgin Islands

Unpaved roads are believed to be the primary source of terrigenous sediments being delivered to marine ecosystems around the island of St John in the eastern Caribbean. The objectives of this study were to: (1) measure runoff and suspended sediment yields from a road segment; (2) develop and test two event‐based runoff and sediment prediction models; and (3) compare the predicted sediment yields against measured values from an empirical road erosion model and from a sediment trap. The runoff models use the Green–Ampt infiltration equation to predict excess precipitation and then use either an empirically derived unit hydrograph or a kinematic wave to generate runoff hydrographs. Precipitation, runoff, and suspended sediment data were collected from a 230 m long, mostly unpaved road segment over an 8‐month period. Only 3–5 mm of rainfall was sufficient to initiate runoff from the road surface. Both models simulated similar hydrographs. Model performance was poor for storms with less than 1 cm of rainfall, but improved for larger events. The largest source of error was the inability to predict initial infiltration rates. The two runoff models were coupled with empirical sediment rating curves, and the predicted sediment yields were approximately 0·11 kg per square meter of road surface per centimetre of precipitation. The sediment trap data indicated a road erosion rate of 0·27 kg m^?2 cm^?1. The difference in sediment production between these two methods can be attributed to the fact that the suspended sediment samples were predominantly sand and silt, whereas the sediment trap yielded mostly sand and gravel. The combination of these data sets yields a road surface erosion rate of 0·31 kg m^?2 cm^?1, or approximately 36 kg m^?2 year^?1. This is four orders of magnitude higher than the measured erosion rate from undisturbed hillslopes. The results confirm the importance of unpaved roads in altering runoff and erosion rates in a tropical setting, provide insights into the controlling processes, and provide guidance for predicting runoff and sediment yields at the road‐segment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Storm‐event rainfall–runoff modelling approach for ungauged sites in Taiwan

This work develops a top‐down modelling approach for storm‐event rainfall–runoff model calibration at unmeasured sites in Taiwan. Twenty‐six storm events occurring in seven sub‐catchments in the Kao‐Ping River provided the analytical data set. Regional formulas for three important features of a streamflow hydrograph, i.e. time to peak, peak flow, and total runoff volume, were developed via the characteristics of storm event and catchment using multivariate regression analysis. Validation of the regional formulas demonstrates that they reasonably predict the three features of a streamflow hydrograph at ungauged sites. All of the sub‐catchments in the study area were then adopted as ungauged areas, and the three streamflow hydrograph features were calculated by the regional formulas and substituted into the fuzzy multi‐objective function for rainfall–runoff model calibration. Calibration results show that the proposed approach can effectively simulate the streamflow hydrographs at the ungauged sites. The simulated hydrographs more closely resemble observed hydrographs than hydrographs synthesized using the Soil Conservation Service (SCS) dimensionless unit hydrograph method, a conventional method for hydrograph estimation at ungauged sites in Taiwan. Copyright © 2008 John Wiley & Sons, Ltd.