Unit Hydrograph
Unit Hydrograph
Unit Hydrograph
A hydrograph is a graph showing variation of discharge (Q) with time (T),at a particular
point of stream or river.
Discharge is usually expressed in cumecs or hectare meter per day and time is expressed
in hours or days. Discharge is plotted on y-axis and corresponding time is plotted on x-
axis.
COMPONENTS OF HYDROGRAPH:
RISING LIMB
PEAK OF CREST SEGMENT
RECCESSION CURVE
Area
Shape
Slope
Rock type
Soil
Land use
Drainage density
Precipitation/Temp
Tidal conditions
The Unit Hydrograph (UH) is the simplest but at the same time a very powerful tool for
hydrological analysis in general and flood forecasting in particular. The unit hydrograph may
be defined as the direct runoff (outflow) hydrograph resulting from one unit of effective
rainfall, which is uniformly distributed over the basin at a uniform rate during a specified
period of time known as unit time or unit duration. The following paragraphs make this
statement still clearer.
Effective rainfall should be uniformly distributed over the basin, i.e. if there are five
rain gauges in the basin, which represent the areal distribution of rainfall over the basin, all
the five rain gauges should record for almost same amount of rainfall during specified time.
A watershed shown on the right here fully marks this stipulation, while converse is true in
respect of left one.
1
In addition, effective rainfall should be at a uniform rate during the unit duration. If the
average rainfall over a particular basin during 6 hour is 126mm, a unit hydrograph of 6 hours
duration can be derived only if the intensity of rainfall is more or less 21 mm/hour over 6
hours. If the same amount of rainfall is distributed with varied intensity, the unit hydrograph
cannot be precisely estimated by simple method.
The unit quantity of effective rainfall is normally taken as 1mm or 1cm; and the outflow
hydrograph is expressed by discharge in cumec. The unit duration may be of 1hour duration
or more, depending upon the size of the catchment, storm characteristics and operational
facilities. However, the unit duration cannot be more than the time of concentration or basin
lag or period of rise. The concept of time of concentration has been covered in detail later in
the chapter.
2
ASSUMPTIONS IN UNIT HYDROGRAPH THEORY
The following are the basic assumption in the unit hydrograph theory:
The unit hydrograph theory assumes the principle of time invariance. This implies
that the direct runoff hydrograph from a given drainage basin due to a given pattern of
effective rainfall will be always same irrespective of the time, i.e. even if the basin
characteristics change with season etc., the unit hydrograph remains the same.
Unit Hydrograph theory assumes the principle of linearity, superimposition or
proportionality. It means that:
If the ordinates of a unit hydrograph of say 1 hour duration are 0,1,6,4,3,2,1,0 units
respectively, the effective rainfall of 2 units falling in 1 hour will produce a direct
runoff hydrographs having ordinates of 0,2,12,8,6,4,2,0 units.
Secondly, if the effective rainfall of two units occurs in 2 hours, i.e. 1 unit per hour,
the direct runoff hydrograph ordinates will be obtained by summing up the
corresponding ordinates of the two unit hydrographs as shown here.
The unit hydrograph theory cannot be applied to catchment areas greter than 5000km sq.
The unit hydrograph theory cannot be applied to very small catchments with area less than
2km sq..
This theory cannot be applied when the major portion of the storm is in the form of snow.
This theory is not very accurate. the accuracy obtained is +10%.
3
A large number of unit hydrographs are required for a catchment because as the unit
duration changes, the unit hydrograph also changes.
The UH is best derived from the observed hydrograph resulting from a storm which
fulfils the two basic conditions i.e., the rainfall is more or less uniformly distributed
over the basin and has a reasonably uniform intensity. Such a hydrograph will generally form an
isolated peak.
In case, such a hydrograph is not available, the UH has to be derived from the analysis
of an observed multipeaked flood hydrograph resulting from several spells of rainfall
of varying intensities.
When the observed discharge and rainfall data at short interval are not available, the
synthetic UH is derived with the help of basin characteristics.
The steps involved in derivation of UH from the analysis of the flood hydrograph with a
single peak are as follows;
1. Inspect discharge records at watershed outlet and corresponding rainfall events to identify
events exhibiting isolated, well defined and single peak with considerable runoff volume.
Pick up as many sets of such records as available. A plot displaying rainfall and
corresponding rise in flood hydrograph, such as here, can help selection of records.
4
Note: A no break/continuous discharge series, as shown in the plot, is developed by
transforming hourly river stage (also water level) into discharge with the help of rating
equation. Rating equation used for this purpose must be developed for the period to which
flood event belongs to.
A rating equation/curve is an equation that relates discharge with water level observed at a
site, and in mathematically expressed as
Q = c * (G Go)n
Where, c, Go; & n are constants; and G is water level.
While gathering information as listed above, it is recommended that
Storms with rainfall should have been active for duration of around 20 to 30 % of
basin lag. Various studies estimate basin lag as 5075% of Tc, time of concentration.
Later part of this chapter describes ways to estimate Tc
Storms should have generated rainfall excess between 1 cm and 4.5 cm.
2.
A flood hydrograph is a basin (catchment) response driven by occurrence of rainfall event
plus contribution of base flow. Had there been no rainfall over the basin, 'Bulge' (rise) in
flood hydrograph would have not appeared. Secondly, all water that falls over a catchment
does not reach the river/stream because of 'losses'; and only a fraction of it contributes to this
'Bulge'. This bulge is termed as Direct Runoff Hydrograph (DRH). The part that reaches the
stream is called as 'Rainfall excess'. Hydrologist seeks to develop a relationship between
'rainfall excess' and DRH. Apparently, the next step is separation of base flow from flood
hydrograph to compute volume of DRH. Following are couple of methods outlined for
separation of base flow.
5
Fixed base method (ABD)
This method suggests the extension of the base flow line along its general trend before the
rise of the hydrograph up to a point B directly below the runoff hydrograph peak. From B, a
straight line BD is drawn to meet the hydrograph at point D, which is N days away from B in
the time
N (in days) = 0.83 A0.2
Where, A is the area of the drainage basin in square kilometersscale. 'N' is determined by an
empirical relation by Linsley as:
This method requires identification of two additional points on the recession limb of
hydrograph one is inflection point; while the other is point E. At inflection point, curve
changes its concavity. This point also indicates end of surface flow to river. This point
beyond, discharge is a combination of interflow (also called as subsurface flow) and base
flow. After a while, interflow also ceases; and only base flow remains in the river. The 'E'
suggests this stage. Once, these two points are located on the graph, a line from 'E' is drawn
6
backward to meet a vertical line from inflection point. A line ABCE divides the DRH and
base flow. Nevertheless, for flood studies, the base flow component is rather insignificant and
hence does not influence the magnitude of peak runoff substantially. Therefore, inaccuracies
involved in separation of base flow are not crucial in overall flood studies.
3. Computation of direct runoff hydrograph ordinates by deducting base flow ordinates from
that of the corresponding observed flood hydrograph.
4. Scanning and analysis of the rainfall data of all rain gauge stations in and around the basin
with a view to
Obtaining areal rainfall over the catchment by appropriate methods, such as Thiessen
Polygon or Isohyetal technique, and
Estimating phiindex. Volume of DRH equals the product of catchment area and
rainfall excess over the basin. This simple analogy helps us estimate depth of rainfall
excess.
5. As DRH is a consequence of given rainfall excess, say 'x' unit. Estimation of the ordinates
of the UH is obtained by dividing the ordinates of direct runoff hydrograph by 'x' rainfall
excess.
6. This process is repeated for all records picked up for this purpose.
7. It is highly probable that UHs derived for more than one record may differ in duration of
excess rainfall. They need to be converted to an identical duration before attempting step 8. A
discussion on conversion of UH duration has been added toward the end of this module.
7
8. All such UHs are eventually averaged. For this, first peaks, Qp, of all UH are averaged to
give Qp, followed by time to peak, tp and time base, Tb of UHs. All other ordinates are
adjusted in such a way that total runoff volume of UH equals the product of 1cm/mm and
catchment area.
Collins'-Method
Matrix-method
Instantaneous Unit Hydrograph
The discrete convolution equation allows the computation of direct runoff, Qn given excess
rainfall, Pm and the unit hydrograph, Un-m+1.
The reverse process, called deconvolution, can be utilised to derive a unit hydrograph given
data on Pm and Qn. Suppose that there are 'M' pulses or burst of rainfall excess and 'n' pulse
of direct runoff in the storm considered; then N equation can be written for Qn, n = 1,2,….,n,
in terms of (n-m+1) unknown values of the unit hydrograph. If Qn and Pm are given and
8
Un-m+1 is required, the set equations is over determined, because there are more equations
(N) than unknowns (n-m+ 1). The term n <= M in the equation restrains the total nos. of P*U
terms for Qn. In first case, when n is <=M, m = 1, 2, .....n; while in case n is >M, m = 1,2,
….M.
(For more details, reader may refer to Applied Hydrology by Ven Te Chow) Let us derive
term for Q1, assuming total number of rainfall pulse, M = 3 Here, n =1 (M =3), hence, m
=1, therefore,
9
UNIT HYDROGRAPH BY CLARK MODEL
The Clark model uses two parameters, time of concentration, 'Tc', and storage constant, 'K',
and a time area histogram concept. Before, we set out for UH by Clark model, let us first
familiarize ourselves with these new terminologies.
10
Where;
Q1=correspond to the discharge, after separating base flow, at the point of inflection on the
recession limb of flood hydrograph.
Q2=correspond to the discharge, after separating base flow, after time t on the recession limb
of flood hydrograph.’
t=time interval between Q1 and Q2.
Adopting a suitable base flow, the value of K can be computed for different hydrographs and
an average value can be worked out.
The timearea diagram represents the areas that will contribute to the flow at the outlet over
successive periods of time. Once the time of concentration is known, lines of equal time
interval called isochrones can be drawn with assumption that time of travel is directly
proportional to distance from the outlet to isochrones. In the picture shown here, dashed lines
in yellow mark isochrones. The US Army Corps of Engineers recommends following
formula to develop time area table/diagram for estimation of inflow from areas bounded
between successive isochrones. An example at the end of this discussion explains the use of
this equation.
The timearea diagram is considered as the inflow to a hypothetical reservoir (S = KO; this
implies that there is absence of wedge storage) and routed through the reservoir to obtain the
outflow hydrograph which is the required instantaneous UH for the basin. Before routing,
inflow from incremental areas between isochrones is converted into discharge units by
following equation;
Where 'ai' is the area in km2 and 't' is the routing period in hours.
Here, we will take a pause to understand as to how this model produces an outflow with
timearea concept; and inflow is generated by each timearea zone due to instantaneous 1mm
effective rainfall. Let us consider the uppermost part of the catchment. Being uppermost part
of the watershed, it does not receive any outflow (O0=0). Instead, it produces I1 runoff which
takes 1 hr (if tc for the catchment is 6hrs; and isochrones are separated by 1hr each) to reach
at the tip of area just below it with a magnitude of Oi routed by following equation.
Oi = CIi + (1-C) Oi – 1
11
Where C = (t/(K+0.5t)), here, t is routing interval in hr, k is basin storage coefficient; and Ii
and Oi are the inflow and outflow at the end of period ti.
In other words, for uppermost area, outflow is obtained by
Now, let us consider the area downstream of first one. Like above, this part generates inflow,
I2 and receives O1 as outflow from upper area. Using routing equation, this area will produce
O2 as
This process is continued till we reach terminal point of the catchment/study area. The IUH
can be converted to a unit hydrograph of same unit duration as routing interval simply
averaging two instantaneous hydrographs lagged by the selected duration that is,
To obtain unit hydrograph for durations other than routing interval (provided that it is exact
multiple of routing interval t. The following equation is used.
Derivation of IUH
There are various methods for the determination of an IUH from the given effective rainfall
hyetograph and direct runoff hydrograph. But the most common is the model suggested by
Nash in 1957. Nash proposed a conceptual model by considering a drainage basin as 'n'
identical linear reservoirs in series. By routing a unit inflow through the reservoirs a
mathematical equation for IUH can be derived.
The ordinate of the IUH at time t is given by,
Where,
n = no. Of the reservoir; and
K= a reservoir constant, also called as storage coefficient.
The values of K and n in Nash model can be evaluated by the method of moments by using
the following relations.
12
Where,
MDRH1= First moment arm of Direct Runoff Hydrograph (DRH)
MERH1= First moment arm of Effective Rainfall Hyetograph (ERH)
MDRH2= Second moment arm of DRH
MERH2= Second moment arm of ERH
The unit of the ordinates of IUH is per sec (sec1). When the ordinates are multiplied by the
total volume of runoff (in cubic meters) resulting from 1mm of rainfall over the catchment
area, the unit will be cumecs.
For finding the unit hydrograph from IUH, the area under the IUH is plotted with respect to
time at the point. The entire area from the start of IUH at different time interval gives points
of S Curve. If a unit hydrograph of T hour duration is required, the S Curve so arrived at is
shifted by T hour and the difference in the ordinates of the two S Curves is computed and
divided by T. The resulting curve forms the unit hydrograph of T hour duration. To illustrate
derivation of UH by Nash method, an example is presented ahead.
in which tp = the basin lag in hours, from midpoint of effective rainfall duration tr to peak of
a unit graph:
L= the length of the main stream from the outlet to the divide in kms;
Lc =the distance from the outlet to a point on the stream nearest to the centroid of the basin;
and Ct = a coefficient
The location of the center of area may be determined by cutting the basin outline from
cardboard and marking the point of intersection of plumb lines drawn with the map
13
suspended from different corners. The coefficient Ct varies from 1.0 to 2.2 with lower values
associated with basins of steeper slopes.
For the standard duration of effective rainfall tr , Snyder proposed:
For the rains of this duration, he found that synthetic unit hydrograph peak Qp
in cumecs may be obtained from the equation:
For the time base T (in days) of the synthetic unit hydrograph U.S. Army Corps of Engineer
adopted the following expression
Where, N=0.38 and C =1.2 for mountainous drainage areas; 0.72 for foothills;
and 0.35 for valley areas.
If unit hydrographs are available for several areas adjacent to a basin for which a unit
hydrograph is required but for which necessary data are lacking, then transposition of
available unit hydrograph will ordinarily give better results than resorting to a synthetic
procedure. Sherman originally proposed that the ordinates and abscissas of unit hydrograph
for similar basins might be assumed to be proportional to the square roots respective drainage
areas. Further details are available in any textbook on applied hydrology.
14
Based on set of equations recommended by Flood Estimation Reports (FER)
FER for 26 hydrometeorologically homogenous subzones in India are reports jointly brought
by Central Water Commission (CWC); India Meteorological Department (IMD); Research,
Design & Standard Organization (RDSO); and Ministry of Shipping & Transport (MoST). In
each report, a number of mathematical relationships between physiographic parameters and
components of unit hydrograph, derived on multiple regression technique, exist. The
stepbystep procedures as to how to develop UH based on FER are illustrated in paragraphs to
follow.
Step-1
Physiographic parameters
1. Location of catchment area to be identified from Survey of India toposheet and measure
the catchment area (A)
2. Measure the length of the longest stream in Km. (L)
3. Length of the longest stream from a point opposite to C.G. of catchment to the point of
study in Km. (Lc)
4. Compute Equivalent Slope in m/Km. (Seq)
To determine equivalent slope, reader may look at following plot which displays longitudinal
profile and formula used for the purpose. We will use this formula a little later to calculate
Seq.
With advances in information technology in recent years and also with the availability of
Digital Elevation Model (DEM), distillation of physiographic parameters is relatively faster
and accurate. Illustrated example in later part of this chapter demonstrates the application of
GIS technique to deduce these parameters.
Step-2
15
Notations/components used in above table are used to construct an SUH using a plot shown
here. If the volume of UH so defined deviates from the volume of runoff generated by
catchment because of 1cm rainfall excess, falling limb of the hydrograph is suitably modified
without altering the points of synthetic parameters such that the volume of UH equals the
theoretical value. The SUH ordinates at one hour interval after corrections are taken as the
final estimate of SUH. Succeeding example presents stepwise procedure to develop SUH of
1hr duration.
1. Storms should be selected with a simple structure with relatively uniform spatial and
temporal distributions.
2. Watershed sizes should generally fall between 1.0 and 100 mi2 in modern watershed
analysis.
3. Direct runoff should range 0.5 to 2 in.
4. Duration of rainfall excess D should be approximately 25% to 30% of lag time tp
5. A number of storms of similar duration should be analyzed to obtain an average UH for
that duration.
6. Step 5 should be repeated for several rainfall of different durations.
APPLICATIONS
Once a unit hydro graph for a duration try has been developed for a basin ,the storm
hydrograph for that basin for any other storm of different intensity but the same duration
can be easily developed.
To draw a maximum flood hydrograph for the design of spillways.
It can be used for flooding and flood warning based on observed rainfall in the basin.
To determine runoff from a storm of unit duration.
16
CONCLUSION
The results pertaining to each approach are provided in the relevant chapters. There exists
numerous approaches based on which design flood in ungauged catchment are estimated.
Only four approaches have been attempted in this study. In the first approach a synthetic unit
hydrograph was attempted without using any meteorological data except the geomorphologic
parameters. The performance was evaluated incorporating the physical process represented
by the channel routing. It revealed that the unit hydrograph should have a region of its
application. There must be upper and lower limit of the catchment size and slope where such
a unit hydrograph concept is applicable. In the second approach, a rainfall based
methodology was developed. A PDF based model was prepared using the rainfall distribution
of 3, 6, 12 and 24 hours duration. Equating the statistical shape variables of both the recorded
and the fitted rainfall event, design hyetograph of given duration was developed. This design
hyetograph when applied as an input to a rainfall-runoff model or unit hydrograph provides
the design flood hydrograph of same return period. Next, runoff based analysis of the gauged
data of a head catchment in Baitarani Sub-basin was carried out. Both the AMS and the PDS
approaches were employed to determine flood quantile of different return period. In the AMS
analysis only the yearly maximum peaks were considered while in PDS analysis all peaks
above a threshold discharge value are taken into account. While a single probability
distribution is capable of modelling the flood peaks in AMS analysis, PDS analysis requires
two probability distributions; one for the number of peaks in a year and the other for the
magnitude of the peaks above a threshold value. But in both the cases, the analyses are
sensitive to the parameters of the probability distributions involved in the corresponding
analysis. After determining the T-year flood quantile at the gauged location, the information
were transferred to the ungauged site adopting a suitable transfer mechanism as per the
requirement.
REFERENCES
Text book of “IRRIGATION AND WATER POWER ENGINEERING” by
Dr.BC.PUNMIA, Dr.PANDE B.B.LAL, ASHOK KUMAR JAIN, ARUN KUMAR
JAIN.
www.wikipedia.com
17