0% found this document useful (0 votes)
31 views19 pages

EECQ - 4242 - Stream Flow Hydrograph Analysis - Part 2

The unit hydrograph is defined as the direct runoff hydrograph resulting from 1 cm of uniformly distributed excess rainfall over the drainage area for an effective duration. It is based on assumptions of linearity, time-invariance, and spatial uniformity. The unit hydrograph can be derived from rainfall and runoff data using the principle of superposition and solving the discrete convolution equation. An example shows how the half-hour unit hydrograph is obtained by solving the set of equations relating excess rainfall pulses to runoff pulses. The derived unit hydrograph ordinates may then be adjusted proportionally to ensure the total runoff depth equals 1 cm.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
31 views19 pages

EECQ - 4242 - Stream Flow Hydrograph Analysis - Part 2

The unit hydrograph is defined as the direct runoff hydrograph resulting from 1 cm of uniformly distributed excess rainfall over the drainage area for an effective duration. It is based on assumptions of linearity, time-invariance, and spatial uniformity. The unit hydrograph can be derived from rainfall and runoff data using the principle of superposition and solving the discrete convolution equation. An example shows how the half-hour unit hydrograph is obtained by solving the set of equations relating excess rainfall pulses to runoff pulses. The derived unit hydrograph ordinates may then be adjusted proportionally to ensure the total runoff depth equals 1 cm.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 19

EECQ 4242: Eng.

Hydrology 1B

Unit hydrograph
1.2. The Unit Hydrograph
The unit hydrograph is the unit pulse response function of a linear hydrologic system. The unit
hydrograph of a watershed is defined as:
a direct runoff hydrograph (DRH) resulting from 1 cm (in SI units) of excess rainfall generated
uniformly over the drainage area at a constant rate for an effective duration.
The unit hydrograph is a simple linear model that can be used to derive the hydrograph resulting from
any amount of excess rainfall. The following basic assumptions are inherent in this model:
1. The excess rainfall has a constant intensity within the effective duration.
2. The excess rainfall is uniformly distributed throughout the whole drainage area.
3. The base time of the DRH (the duration of direct runoff) resulting from an excess rainfall of a
given duration is constant.
4. The ordinates of all DRH's of a common base time are directly proportional to the total amount
of direct runoff represented by each hydrograph.
5. For a given watershed, the hydrograph resulting from a given excess rainfall reflects the
unchanging characteristics of the watershed.
Under natural conditions, the above assumptions cannot be perfectly satisfied. However, when the
hydrologic data to be used are carefully selected so that they come close to meeting the above
assumptions, the results obtained by the unit hydrograph model are generally acceptable for practical
purposes
Assumption 1: the storms selected for analysis should be of short duration, since these will most likely
produce an intense and nearly constant excess rainfall rate, yielding a well-defined single-peaked
hydrograph of short time base.
Assumption 2: the unit hydrograph may become inapplicable when the drainage area is too large to
be covered by a nearly uniform distribution of rainfall. In such cases, the area has to be divided and
each subarea analyzed for storms covering the whole subarea.
Assumption 3: the base time of the direct runoff hydrograph (DRH) is generally uncertain but depends
on the method of baseflow separation. The base time is usually short if the direct runoff is considered
to include the surface runoff only; it is long if the direct runoff also includes subsurface runoff.
Assumption 4: the principles of superposition and proportionality are assumed so that the ordinates
𝑄𝑛 of the DRH may be computed by equation 22. However, actual hydrologic data are not truly linear
and when you apply equation 22, the resuting hydrograph is only an approximation. It is however
satisfactory in many practical cases.
Assumption 5:, the unit hydrograph is considered unique for a given watershed and invariable with
respect to time. This is the principle of time invariance, which, together with the principles of
superposition and proportionality, is fundamental to the unit hydrograph model. Unit hydrographs

1
EECQ 4242: Eng. Hydrology 1B

are applicable only when channel conditions remain unchanged, and watersheds do not have
appreciable storage. This condition is violated when the drainage area contains many reservoirs, or
when the flood overflows into the flood plain, thereby producing considerable storage.
The principles of linear system analysis form the basis of the unit hydrograph method. Table 1 shows
a comparison of linear system concepts with the corresponding unit hydrograph concepts. In
hydrology, the step response function is commonly called the S-hydrograph, and the impulse
response function is called the instantaneous unit hydrograph which is the hypothetical response to
a unit depth of excess rainfall deposited instantaneously on the watershed surface.

1.3.1. Derivation of the Unit Hydrograph


The discrete convolution equation (22) allows the computation of direct runoff 𝑄𝑛 given excess rainfall
𝑃𝑚 and the unit hydrograph 𝑈𝑛−𝑚+1
𝑛≤𝑀

𝑄𝑛 = ∑ 𝑃𝑚 𝑈𝑛−𝑚+1 (1)
𝑚=1

The reverse process, called deconvolution, is needed to derive a unit hydrograph given data on 𝑃𝑚 and
𝑄𝑛 .
Suppose that there are M pulses of excess rainfall and N pulses of direct runoff in the storm considered;
then N equations can be written for 𝑄𝑛 , 𝑛 = 1, 2, … , 𝑁 in terms of 𝑁 − 𝑀 + 1 unknown values of the
unit hydrograph as shown in Table 2
Table 2. The set of equations for discrete time convolution
𝑛≤𝑀

𝑄𝑛 = ∑ 𝑃𝑚 𝑈𝑛−𝑚+1 ; 𝑛 = 1, 2, … , 𝑁
𝑚=1

If 𝑄𝑛 and 𝑃𝑚 are given and 𝑈𝑛−𝑚+1 is required, the set of equations in Table 2 is overdetermined,
because there are more equations (N) than unknowns (𝑁 − 𝑀 + 1).

2
EECQ 4242: Eng. Hydrology 1B

Table 1: Comparison of linear system and unit hydrograph concepts

3
EECQ 4242: Eng. Hydrology 1B

Example 1.3.1. Find the half-hour unit hydrograph using the excess rainfall hyetograph and direct
runoff hydrograph given in Table 3 below.
Solution. The ERH and DRH in Table 3 have 𝑀 = 3 𝑎𝑛𝑑 𝑁 = 11 pulses respectively.
The number of pulses in the unit hydrograph is 𝑁 − 𝑀 + 1 = 11 − 3 + 1 = 9
Substituting the ordinates of the ERH and DRH into the equations in Table 2 yields a set of 11
simultaneous equations.
These equations may be solved by Gauss elimination to give the unit hydrograph ordinates. Gauss
elimination involves isolating the unknown variables one by one and successively solving for them. In
this case, the equations can be solved from top to bottom, working with just the equations involving
the first pulse 𝑃1, starting with
𝑄1 428
𝑈1 = = = 404 𝑐𝑓𝑠⁄𝑖𝑛
𝑃1 1.06
Table 3. Excess rainfall hyetograph and direct runoff hydrograph
Time (1/2 h) Excess rainfal Direct runoff
(in) (cfs)
1 1.06 428
2 1.93 1923
3 1.81 5297
4 9131
5 10625
6 7834
7 3921
8 1846
9 1402
10 830
11 313

𝑄2 − 𝑃2 𝑈1 1923 − 1.93 × 404


𝑈2 = = = 1079 𝑐𝑓𝑠⁄𝑖𝑛
𝑃1 1.06
𝑄3 − 𝑃3 𝑈1 − 𝑃2 𝑈2 5297 − 1.81 × 404 − 1.93 × 1079
𝑈3 = = = 2343 𝑐𝑓𝑠⁄𝑖𝑛
𝑃1 1.06
And similarly, for the remaining ordinates:
9131 − 1.81 × 1079 − 1.93 × 2343
𝑈4 = = 2506 𝑐𝑓𝑠⁄𝑖𝑛
1.06
10625 − 1.81 × 2343 − 1.93 × 2506
𝑈5 = = 1460 𝑐𝑓𝑠⁄𝑖𝑛
1.06
7834 − 1.81 × 2506 − 1.93 × 1460
𝑈6 = = 453 𝑐𝑓𝑠⁄𝑖𝑛
1.06
3921 − 1.81 × 1460 − 1.93 × 453
𝑈7 = = 381 𝑐𝑓𝑠⁄𝑖𝑛
1.06
4
EECQ 4242: Eng. Hydrology 1B

1846 − 1.81 × 453 − 1.93 × 381


𝑈8 = = 274 𝑐𝑓𝑠⁄𝑖𝑛
1.06
1402−1.81×381−1.93×274
𝑈9 = = 173 𝑐𝑓𝑠⁄𝑖𝑛
1.06

Table 4. the ½ hour Unit hydrograph derived


n 1 2 3 4 5 6 7 8 9
𝑈𝑛 (𝑐𝑓𝑠⁄𝑖𝑛) 404 1079 2343 2506 1460 453 381 274 173

The derived unit hydrograph is given in table 4. Solutions may be similarly obtained by focusing on
other rainfall pulses. The depth of direct runoff in the unit hydrograph can be checked and found to
equal to 1.00 inch as required
In cases where the derived unit hydrograph does not meet this requirement, the ordinates are
adjusted by proportion so that the depth of direct runoff is 1inch (or cm)

In general, the unit hydrographs obtained by solutions of the set of equations in Table 2 for different
rainfall pulses are not identical. To obtain a unique solution a method of successive approximation can
be used. It involves four steps:
1. assume a unit hydrograph, and apply it to all excess rainfall blocks of the hyetograph except
the largest
2. subtract the resulting hydrograph from the actual DRH, and reduce the residual to unit
hydrograph terms
3. compute a weighted average of the assumed unit hydrograph and the residual unit hydrograph,
and use it as the revised approximation for the next trial
4. repeat the previous three steps until the residual unit hydrograph does not differ by more than
a permissible amount from the assumed hydrograph.
The resulting unit hydrograph may show erratic variations and even have negative values. If this occurs,
a smooth curve may be fitted to the ordinates to produce an approximation of the unit hydrograph.
Erratic variation in the unit hydrograph may be due to nonlinearity in the effective rainfall-direct runoff
relationship in the watershed, and even if this relationship is truly linear, the observed data may not
adequately reflect this. Also, actual storms are not always uniform in time and space, as required by
theory, even when the excess rainfall hyetograph is broken into pulses of short duration.

1.3.2. Unit hydrograph application


The unit hydrograph may be applied to find the direct runoff and streamflow hydrographs. The
following steps are applied:
a. A rainfall hyetograph is selected
b. The abstractions are estimated
c. The excess rainfall hyetograph is calculated

5
EECQ 4242: Eng. Hydrology 1B

The time interval used in defining the excess rainfall hyetograph ordinates must be the same
as that for which the unit hydrograph was specified.
d. The discrete convolution equation may then be used to yield the direct runoff hydrograph.
e. By adding an estimated baseflow to the direct runoff hydrograph, the streamflow hydrograph
is obtained.

Example. Calculate the streamflow hydrograph for a storm of 6 in excess rainfall, with 2 in in the first
half-hour, 3 in in the second half-hour and 1 in in the third half-hour. Use the half-hour unit hydrograph
computed in the previous Example and assume the baseflow is constant at 500 cfs throughout the
flood. Check that the total depth of direct runoff is equal to the total excess precipitation
(𝑤𝑎𝑡𝑒𝑟𝑠ℎ𝑒𝑑 𝑎𝑟𝑒𝑎 = 7.03 𝑚𝑖 2 ).
Solution. The calculation of the direct runoff hydrograph by convolution is shown in Table 5. The unit
hydrograph ordinates from Table 3 are laid out along the top of the table and the excess precipitation
depths down the left side.
The time interval is in ∆𝑡 = 0.5 ℎ intervals.
For the first time interval, 𝑛 = 1 in Equation (1), and
𝑄1 = 𝑃1 𝑈1 = 2.00 × 404 = 808 𝑐𝑓𝑠
For the second time interval,
𝑄2 = 𝑃2 𝑈1 + 𝑃1 𝑈2 = 3.00 × 404 + 2.00 × 1079 = 1212 + 2158 = 3370 𝑐𝑓𝑠
𝑄1 = 𝑃1 𝑈1 = 2.00 × 404 = 808 𝑐𝑓𝑠
For the third time interval,
𝑄3 = 𝑃3 𝑈1 + 𝑃2 𝑈2 + 𝑃2 𝑈2 = 1.00 × 404 + 3.00 × 1079 + 2.00 × 2343
= 404 + 3237 + 4686 = 8327 𝑐𝑓𝑠
Etc.

The total direct runoff volume is


𝑁

𝑉𝑑 = ∑ 𝑄𝑛 ∆𝑡 = 54438 × 0.5 𝑐𝑓𝑠. ℎ


𝑛=1

𝑓𝑡 3 . ℎ 3600 𝑠
= 54438 × 0.5 × = 9.80 × 107 𝑓𝑡 3
𝑠 1ℎ
Check with the corresponding depth of direct runoff:
The corresponding depth of direct runoff is found by dividing by the runoff volume by the watershed
area
𝑊𝑎𝑡𝑒𝑟𝑠ℎ𝑒𝑑 𝑎𝑟𝑒𝑎, 𝐴 = 7.03 𝑚𝑖 2 = 7.03 𝑥 52802 𝑓𝑡 2 = 1.96 × 108 𝑓𝑡 2

6
EECQ 4242: Eng. Hydrology 1B
Table 5. Calculation of the direct runoff hydrograph and streamflow hydrograph for the example
*baseflow = 500 cfs
Time Excess Unit hydrograph ordinates (cfs/in) Direct Streamflow
(1/2 h) rainfall runoff (cfs)
1 2 3 4 5 6 7 8 9
(in) (cfs)
404 1079 2343 2506 1460 453 381 274 173
n=1 2.00 808 808 1308
2 3.00 1212 2158 3370 3870
3 1.00 404 3237 4686 8327 8827
4 1079 7029 5012 13120 13620
5 2343 7518 2920 12781 13281
6 2506 4380 906 7792 8292
7 1460 1359 762 3581 4081
8 453 1143 548 2144 2644
9 381 822 346 1549 2049
10 274 519 793 1293
11 173 173 673
Total 54438

7
EECQ 4242: Eng. Hydrology 1B

𝑉𝑑 9.80 × 107
𝑟𝑑 = = = 0.500 𝑓𝑡 = 6.00 𝑖𝑛
𝐴 1.96 × 108
which is equal to the total depth of excess precipitation as required.
The streamflow hydrograph is found by adding the 500 cfs baseflow to the direct runoff hydrograph,
as shown on the right-hand side of Table 5 and graphically in Figure 1.

Figure 1. Streamflow hydrograph from a storm with excess rainfall pulses of duration 0.5 h and amount
2 in, 3 in, and 1 in, respectively. Total streamflow = baseflow + direct runoff (Previous Example).

1.3.3. Other methods of deriving the unit hydrograph:


1. Matrix calculation
2. Linear regression
3. Linear programming

1.3.4. Synthetic unit hydrograph


The unit hydrograph developed from rainfall and streamflow data on a watershed applies only for that
watershed and for the point on the stream where the streamflow data were measured.
Synthetic unit hydrograph procedures are used to develop unit hydrographs for other locations on the
stream in the same watershed or for nearby watersheds of a similar character. There are three types
of synthetic unit hydrographs:
1. Those relating hydrograph characteristics (peak flow rate, base time, etc.) to watershed
characteristics
2. Those based on a dimensionless unit hydrograph
3. Those based on models of watershed storage. (this will not be discussed in this course)

8
EECQ 4242: Eng. Hydrology 1B

Snyder’s synthetic unit hydrograph


This is based on an empirical study of Snyder in 1938, and later by U.S Army Corps of Engineers in 1938
on watersheds varying in size from 30 – 30000 km2.
They found synthetic relations for some characteristics of a standard unit hydrograph [Fig.2].

Figure 2. Snyder’s synthetic unit hydrograph (a) standard unit hydrograph (𝑡𝑝 = 5.5 𝑡𝑟 ) (b) required
unit hydrograph (𝑡𝑝 𝑟 ≠ 5.5 𝑡𝑟 )
From the relations, five characteristics of a required unit hydrograph [Fig. 2] for a given excess rainfall
duration may be calculated:
Snyder defined a standard unit hydrograph as one whose rainfall duration 𝑡 𝑟 is related to the basin lag
𝑡𝑝 by:
𝑡𝑝 = 5.5 𝑡𝑟 (2)
1. The basin lag is
𝑡𝑝 = 𝐶1 𝐶𝑡 (𝐿𝐿𝑐 )0.3 (3)
where 𝑡𝑝 is in hours, L is the length of the main stream in kilometers (or miles) from the outlet to the
upstream divide, 𝐿𝑐 is the distance in kilometers (miles) from the outlet to a point on the stream
nearest the centroid of the watershed area, 𝐶1 = 0.75, and 𝐶𝑡 is a coefficient derived from gaged
watersheds in the same region.
2. The peak discharge per unit drainage area in 𝑚3 ⁄𝑠 . 𝑘𝑚2 of the standard unit hydrograph is
𝐶2 𝐶𝑝
𝑞𝑝 = (4)
𝑡𝑝
where 𝐶2 = 2.75 and 𝐶𝑝 is a coefficient derived from gaged watersheds in the same region.
9
EECQ 4242: Eng. Hydrology 1B

To compute 𝐶𝑡 and 𝐶𝑝 for a gaged watershed, the values of L and 𝐿𝑐 are measured from the basin map.
From a derived unit hydrograph of the watershed are obtained values of its effective duration 𝑡𝑟 in
hours, its basin lag 𝑡𝑝 𝑟 in hours, and its peak discharge per unit drainage area, 𝑞𝑝 𝑟 𝑖𝑛 𝑚3 ⁄𝑠 . 𝑘𝑚2 . 𝑐𝑚

If (𝑡𝑝 𝑟 = 5.5 𝑡𝑟 ), 𝑡ℎ𝑒𝑛 𝑡𝑟 = 𝑡𝑝 𝑎𝑛𝑑 𝑞𝑝 𝑟 = 𝑞𝑝 ,


and 𝐶𝑡 and 𝐶𝑝 are computed by Eqs. (3) and (4).
If 𝑡𝑝 𝑟 is quite different from 5.5 𝑡𝑟 , the standard basin lag is
𝑡𝑟 − 𝑡𝑅
𝑡𝑝 = 𝑡𝑝 𝑟 + (5)
4
and Eqs. (2) and (5) are solved simultaneously for 𝑡𝑟 𝑎𝑛𝑑 𝑡𝑝 . The
values of 𝐶𝑡 and 𝐶𝑝 are then computed from (3) and (4) with 𝑞𝑝 𝑟 = 𝑞𝑝 and 𝑡𝑝 𝑟 = 𝑡𝑝 .
When an ungaged watershed appears to be like a gaged watershed, the coefficients 𝐶𝑡 and 𝐶𝑝 for the
gaged watershed can be used in the above equations to derive the required synthetic unit hydrograph
for the ungaged watershed.
3. The relationship between 𝑞𝑝 and the peak discharge per unit drainage area 𝑞𝑝 𝑟 of the required
unit hydrograph is
𝑞𝑝 𝑡𝑝
𝑞𝑝 𝑅 = (6)
𝑡𝑝 𝑟
4. The base time 𝑡𝑏 in hours of the unit hydrograph can be determined using the fact that the area
under the unit hydrograph is equivalent to a direct runoff of 1 cm. Assuming a triangular shape
for the unit hydrograph, the base time may be estimated by
𝐶3
𝑡𝑏 = ; 𝐶3 = 5.56 (7)
𝑞𝑝 𝑅
5. The width in hours of a unit hydrograph at a discharge equal to a certain percent of the peak
discharge 𝑞𝑝 𝑅 is given by
𝑊 = 𝐶𝑤 𝑞𝑝 𝑅−1.08 (8)
where 𝐶𝑤 = 1.22 for the 75-percent width, and 2.14 for the 50-percent width. Usually one-third of
this width is distributed before the unit hydrograph peak time and two-thirds after the peak.
Example 1. From the basin map of a given watershed, the following quantities are measured:
L = 150 km, L c = 75 km, Drainage area = 3500 km2
From the unit hydrograph derived for the watershed, the following are determined:
𝑡 𝑅 = 12 ℎ , 𝑡 𝑝 𝑅 = 34 ℎ, 𝑝𝑒𝑎𝑘 𝑑𝑖𝑠𝑐ℎ𝑎𝑟𝑔𝑒 = 157.5 𝑚3 ⁄𝑠
Determine the coefficients𝐶𝑡 and 𝐶𝑝 for the synthetic unit hydrograph of the watershed.
Solution:

10
EECQ 4242: Eng. Hydrology 1B

𝐹𝑟𝑜𝑚 𝑡ℎ𝑒 𝑔𝑖𝑣𝑒𝑛 𝑑𝑎𝑡𝑎, 5.5𝑟𝑅 = 66 ℎ, 𝑤ℎ𝑖𝑐ℎ 𝑖𝑠 𝑞𝑢𝑖𝑡𝑒 𝑑𝑖𝑓𝑓𝑒𝑟𝑒𝑛𝑡 𝑓𝑟𝑜𝑚 𝑡 𝑝 𝑅 = 34 ℎ.


𝑡𝑟 − 𝑡𝑅 𝑡𝑟 − 12
𝑡𝑝 = 𝑡𝑝 𝑟 + 𝑖. 𝑒. 𝑡𝑝 = 34 + (𝑖)
4 4
𝑅𝑒𝑐𝑎𝑙𝑙 𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛 2: 𝑡𝑝 = 5.5 𝑡𝑟 (𝑖𝑖)
Solving (i) and (ii) simultaneously:
𝑡𝑟 = 5.9 ℎ 𝑎𝑛𝑑 𝑡𝑝 = 32.5 ℎ
Use equation 3 to calculate 𝐶𝑡
𝑡𝑝 = 𝐶1 𝐶𝑡 (𝐿𝐿𝑐 )0.3 𝑖. 𝑒. 32.5 = 0.75 × 𝐶𝑡 × (150 × 75)0.3
𝐶𝑡 = 2.65
The peak discharge per unit area is:
157.5
𝑞𝑝 𝑅 = = 0.045 𝑚3 ⁄𝑠 . 𝑐𝑚
3500
The coefficient 𝐶𝑝 is calculated by Eq. (4) 𝑤𝑖𝑡ℎ 𝑞𝑝 = 𝑞𝑝 𝑅 𝑎𝑛𝑑 𝑡𝑝 = 𝑡𝑝 𝑅
𝐶2 𝐶𝑝 2.75 × 𝐶𝑝
𝑞𝑝 𝑅 = ; 0.045 = ; 𝐶𝑝 = 0.56
𝑡𝑝 𝑅 34.0
Example 2. Compute the six-hour synthetic unit hydrograph of a watershed given the following:
Drainage area = 2500 km2 ; L = 100 km, L c = 50 km,
This watershed is a sub-drainage area of the watershed in Example 1.
Solution
The values of Ct and Cp determined from exaplme 1 will be used for this watershed:
𝐶𝑡 = 2.65 ; 𝐶𝑝 = 0.56
Calculate 𝑡𝑝 from equation 3:

𝑡𝑝 = 𝐶1 𝐶𝑡 (𝐿𝐿𝑐 )0.3 𝑖. 𝑒. 𝑡𝑝 = 0.75 × 2.65 × (100 × 50)0.3 = 25.5 ℎ


Calculate 𝑡𝑟 from equation 2:
𝑡𝑝 25.5
𝑡𝑝 = 5.5 𝑡𝑟 𝑖. 𝑒. 𝑡𝑟 = = = 4.64 ℎ
5.5 5.5
𝐹𝑜𝑟 𝑎 6 − ℎ𝑜𝑢𝑟 𝑢𝑛𝑖𝑡 ℎ𝑦𝑑𝑟𝑜𝑔𝑟𝑎𝑝ℎ, 𝑡𝑅 = 6 ℎ;
From Equation (9)
𝑡𝑟 − 𝑡𝑅 𝑡𝑟 − 𝑡𝑅 4.64 − 6
𝑡𝑝 = 𝑡𝑝 𝑅 + ∴ 𝑡𝑝 𝑅 = 𝑡𝑝 − = 25.5 − = 25.8 ℎ
4 4 4
𝐶2 𝐶𝑝 2.75 × 0.56
From Equation (4): 𝑞𝑝 = = = 0.0604 𝑚3 ⁄𝑠 . 𝑘𝑚2
𝑡𝑝 25.5
𝑞𝑝 𝑡𝑝 0.0604 × 25.5
𝐹𝑟𝑜𝑚 𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛 (6): 𝑞𝑝 𝑅 = = = 0.597 𝑚3 ⁄𝑠 . 𝑘𝑚2 . 𝑐𝑚
𝑡𝑝 𝑟 25.8
Or
11
EECQ 4242: Eng. Hydrology 1B

𝑇ℎ𝑒 𝑝𝑒𝑎𝑘 𝑑𝑖𝑠𝑐ℎ𝑎𝑟𝑔𝑒 = 0.0597 × 2500 = 149.2 𝑚3 ⁄𝑠 . 𝑐𝑚


The widths of the unit hydrograph are given by equation (8): W = Cwq p R−1.08
At 75% of peak discharge, 𝑊 = 1.22 × 𝑞𝑝 𝑅−1.08 = 1.22 × 0.0597−1.08 = 25.6 ℎ
At 50% of peak discharge, 𝑊 = 2.14 × 𝑞𝑝 𝑅−1.08 = 2.14 × 0.0597−1.08 = 44.9 ℎ
𝐶3
𝑇ℎ𝑒 𝑏𝑎𝑠𝑒 𝑡𝑖𝑚𝑒 𝑖𝑠 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑡𝑒𝑑 𝑢𝑠𝑖𝑛𝑔 𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛 (7): 𝑡𝑏 = ; 𝐶3 = 5.56
𝑞𝑝 𝑅
5.56
𝑡𝑏 = = 93 ℎ
0.0597
The hydrograph from this calculation is shown in the figure 3 below: A check can be done to ensure
that it represents a depth of direct runoff of 1 cm.

Figure 3. Synthetic unit hydrograph calculated by Snyder’s method in example 2.

Regionalization of UH parameters
A further innovation in the use of Snyder's method has been the regionalization of unit hydrograph
parameters. Espey, Altman and Graves (1977) developed a set of generalized equations for the
construction of 10-minute unit hydrographs using a study of 41 watersheds ranging in size from 0.014
to 15 mi2, and in impervious percentage from 2 to 100 percent. Of the 41 watersheds, 16 are located
in Texas, 9 in North Carolina, 6 in Kentucky, 4 in Indiana, 2 each in Colorado and Mississippi, and 1 each
in Tennessee and Pennsylvania. The equations are:
𝑇𝑝 = 3.1𝐿0.23 𝑆 −0.25 𝐼 −0.18 𝜑1.57
𝑄𝑝 = 31.62 × 103 𝐴0.96𝑇 −1.07

𝑇𝐵 = 125.89 × 103 𝐴𝑄𝑝 −0.95

12
EECQ 4242: Eng. Hydrology 1B

𝑊50 = 16.22 × 103 𝐴0.93 𝑄𝑝 −0.92 (9)

Figure 4. Watershed conveyance factor 𝜑 as a function of channel roughness and watershed


imperviousness. (Adapted from Espey, Altman, and Graves, 1977.)
where
L = the total distance (in feet) along the main channel from the point being considered to the upstream
watershed boundary
S = the main channel slope (in feet per foot), defined by H/0.8L, where H is the difference in elevation
between A and B. A is the point on the channel bottom at a distance of 0.2L downstream from the
upstream watershed boundary; B is a point on the channel bottom at the downstream point
being considered
I = the impervious area within the watershed (in percent), assumed equal to 5 percent for an
undeveloped watershed
O = the dimensionless watershed conveyance factor, which is a function of percent impervious and
roughness (Fig. 4)
A = the watershed drainage area (in square miles)
Tp = the time of rise to the peak of the unit hydrograph from the beginning of runoff (in minutes)
Qp = the peak flow of the unit hydrograph (in cfs/in)
TB = the time base of the unit hydrograph (in minutes)
W50 = the width of the hydrograph at 50 percent of Qp (in minutes)
W75 = the width of at 75 percent of Qp (in minutes)

SCS dimensionless unit hydrograph


The SCS dimensionless hydrograph is a synthetic unit hydrograph in which the discharge is expressed
by the ratio of discharge q to peak discharge qp and the time by the ratio of time t to the time of rise
of the unit hydrograph, Tp.

13
EECQ 4242: Eng. Hydrology 1B

Given the peak discharge and lag time for the duration of excess rainfall, the unit hydrograph can be
estimated from the synthetic dimensionless hydrograph for the given basin.
Figure 5 (a) shows such a dimensionless hydrograph, prepared from the unit hydrographs of a variety
of watersheds. The values of qp and Tp may be estimated using a simplified model of a triangular unit
hydrograph as shown in Figure 7.5 (b) where the time is in hours and the discharge in 𝑚3 ⁄𝑠 (Soil
Conservation Service, 1972)
The Soil Conservation Service suggests the time of recession may be approximated as 1.67 Tp. As the
area under the unit hydrograph should be equal to a direct runoff of 1 cm
, it can be shown that:
𝐶𝐴
𝑞𝑝 = (10)
𝑇𝑝
where C = 2.08
A is the drainage area in (𝑘𝑚2 )
Further, a study of unit hydrographs of many large and small rural watersheds indicates that the basin
lag:
𝑏𝑎𝑠𝑖𝑛 𝑙𝑎𝑔 𝑡 𝑝 ≅ 0.6𝑇𝐶 ,
where TC is the time of concentration of teh watershed
As shown in Figure 7.5 (b), time of rise 𝑇 𝑝 can be expressed in terms of lag time 𝑡 𝑝 and the duration
of effective rainfall 𝑡 𝑟 :
𝑡𝑟
𝑇𝑝 = + 𝑡𝑝 (11)
2

Figure 5. Soil Conservation Service synthetic unit hydrographs (a) Dimensionless hydrograph and (b)
triangular unit hydrograph. (Source: Soil Conservation Service, 1972.)

14
EECQ 4242: Eng. Hydrology 1B

Example 7.7.3. Construct a 10-minute SCS unit hydrograph for a basin of area 3.0 𝑘𝑚2 and time of
concentration 1.25 h.
Solution.
𝑇ℎ𝑒 𝑑𝑢𝑟𝑎𝑡𝑖𝑜𝑛 𝑡 𝑟 = 10 𝑚𝑖𝑛 = 0.166 ℎ,
𝐿𝑎𝑔 𝑡𝑖𝑚𝑒 𝑡𝑝 = 0.6𝑇𝑐 = 0.6 𝑥 1.25 = 0.75 ℎ,
0.166
𝑅𝑖𝑠𝑒 𝑡𝑖𝑚𝑒 𝑇𝑝 = 𝑡 𝑟 𝑙2 + 𝑡𝑝 = + 0.75 = 0.833 ℎ.
2
𝐶𝐴 2.08 𝑥 3.0
𝐹𝑟𝑜𝑚 𝐸𝑞. (10), 𝑞𝑝 = = = 7.49 𝑚3 ⁄𝑠 . 𝑐𝑚
𝑇𝑝 0.833
The dimensionless hydrograph in Fig. 5 may be converted to the required dimensions by multiplying
the values on the horizontal axis by 𝑇𝑝 and those on the vertical axis by 𝑞𝑝 . Alternatively, the triangular
unit hydrograph can be drawn with 𝑡𝑏 = 2.61𝑇𝑝 = 2.22 h. The depth of direct runoff is checked to equal
1 cm.

1.3.5 Unit hydrographs for different rainfall durations


When a unit hydrograph of a given excess-rainfall duration is available, the unit hydrographs of other
durations can be derived.
If other durations are integral multiples of the given duration, the new unit hydrograph can be easily
computed by application of the principles of superposition and proportionality.
A general method of derivation applicable to unit hydrographs of any required duration may be used
on the basis of the principle of superposition. This is the S-hydrograph method.
The curve assumes a deformed S shape and its ordinates ultimately approach the rate of excess rainfall
at a time of equilibrium. This step response function g(t) can be derived from the unit pulse response
function h(t) of the unit hydrograph, as follows:
Recall, the response at time t to a unit pulse of duration 1.3.5 beginning at time 0 is:
1
ℎ (𝑡 ) = [𝑔(𝑡) − 𝑔(𝑡 − ∆𝑡) (12)
∆𝑡
Similarly, the response at time t to a unit pulse beginning at time 1.3.5 is equal to
ℎ(𝑡 − ∆𝑡), i. e. ℎ (𝑡) lagged by ∆t time units:
1
ℎ(𝑡 − ∆𝑡) = [𝑔(𝑡 − ∆𝑡) − 𝑔(𝑡 − 2∆𝑡) (13)
∆𝑡
and the response at time t to a third unit pulse beginning at time 2∆𝑡 is
1
ℎ(𝑡 − 2∆𝑡) = [𝑔(𝑡 − 2∆𝑡) − 𝑔(𝑡 − 3∆𝑡) (14)
∆𝑡
Continuing this process indefinitely, summing the resulting equations, and rearranging, yields the unit
step response function, or S-hydrograph, as shown in

15
EECQ 4242: Eng. Hydrology 1B

Fig. 7.6.(a):
𝑔(𝑡) = ∆𝑡[ℎ(𝑡) + ℎ(𝑡 − ∆𝑡) + ℎ(𝑡 − 2∆𝑡) + ⋯ ] (15)

Figure 7.6 Using the S-hydrograph to find a unit hydrograph of duration ∆𝑡 ' from a unit hydrograph of
duration ∆𝑡
where the summation is multiplied by ∆𝑡 so that g(t) will correspond to an input rate of 1, rather than
1⁄∆𝑡 as used for each of the unit pulses.
Theoretically, the S-hydrograph derived should be a smooth curve because the input excess rainfall is
assumed to be at a constant, continuous rate.
The summation process will result in an undulatory form if:
1. there are errors in the rainfall abstractions or baseflow separation,

16
EECQ 4242: Eng. Hydrology 1B

2. the actual duration of excess rainfall is not the derived duration for the unit hydrograph.
3. there is nonuniform temporal and areal distribution of rainfall;
4. when the natural data are not linear, the resulting unstable system oscillations may produce
negative ordinates. In such cases, an optimization technique may be used to obtain a smoother
unit hydrograph.
After the S-hydrograph is constructed, the unit hydrograph of a given duration can be derived as
follows:
Step 1: Advance, or offset, the position of the S-hydrograph by a period equal to the desired duration
∆𝑡 and call this S-hydrograph
an offset S-hydrograph, g'(t) [Fig. 7.6(b)], defined by:
𝑔′ (𝑡) = (𝑔(𝑡 − ∆𝑡 ′) (16)

g'(t) = g(t-At') (7.8.5)


Step 2: The difference between the ordinates of the original S-hydrograph and the offset S-hydrograph,
divided by ∆𝑡 ′, gives the desired unit hydrograph [Fig. 7.6(c)]:
1
ℎ ′ (𝑡 ) = [𝑔(𝑡) − 𝑔(𝑡 − ∆𝑡 ′ )] (17)
∆𝑡 ′
Example 7.8.1. Use the 0.5-hour unit hydrograph in Table 3 (from Example 1.3.1) to produce the S-
hydrograph and the 1.5-h unit hydrograph for this watershed.
Solution.
The 0.5-h unit hydrograph is shown in column 2 of Table 6.
The S-hydrograph is found using equation (15) with ∆𝑡 = 0.5 h.
For t = 0.5 h,
𝑔(𝑡) = ∆𝑡[ℎ(𝑡)] 𝑖. 𝑒. 𝑔(𝑡) = 0.5 × 404 = 202 𝑐𝑓𝑠
For t = 1 h,
𝑔(𝑡) = ∆𝑡[ℎ(𝑡) + ℎ(𝑡 − 0.5)] 𝑖. 𝑒. 𝑔(𝑡) = 0.5 × (1079 + 404) = 742 𝑐𝑓𝑠
For t = 1.5 h,
𝑔(𝑡) = ∆𝑡[ℎ(𝑡) + ℎ(𝑡 − 0.5) + ℎ(𝑡 − 1.0)]
𝑖. 𝑒. 𝑔(𝑡) = 0.5 × (2343 + 1079 + 404) = 1913 𝑐𝑓𝑠
Etc
This is recorded in column 3 of Table 6
The S-hydrograph is offset by ∆𝑡 ′ = 1.5 h (column 4) to give 𝑔(𝑡 − ∆𝑡 ′), and the difference divided
by∆𝑡 ′ to give the 1.5-h unit hydrograph ℎ′ (𝑡)) (column 5). For example:

17
EECQ 4242: Eng. Hydrology 1B

For t = 2.0
1
ℎ (𝑡 ) = [3166 − 202] = 1976 𝑐𝑓𝑠⁄𝑖𝑛
1.5
Table 6. Calculation of a 1.5-h unit hydrograph by the S-hydrograph method (from previous example)

1.3.6. Instantaneous unit hydrograph


If the excess rainfall is of unit amount and its duration is infinitesimally small, the resulting hydrograph
is an impulse response function called the instantaneous unit hydrograph (IUH).
For an IUH, the excess rainfall is applied to the drainage area in zero time. Of course, this is only a
theoretical concept and cannot be realized in actual watersheds, but it is useful because the IUH
characterizes the watershed's response to rainfall without reference to the rainfall duration.
Therefore, the IUH can be related to watershed geomorphology
The convolution integral is
𝑡

𝑄(𝑡) = ∫ 𝑢(𝑡 − 𝜏)𝐼 (𝜏)𝑑𝜏 (1.3.6.1)


0

If the quantities 𝐼(𝜏) 𝑎𝑛𝑑 𝑄(𝑡) have the same dimensions, the ordinate of the IUH must have
dimensions [𝑇 −1 ].
With l = t − τ, the IUH has the following properties:
0 ≤ u(l) ≤ some positive peak value for l > 0
u(l) = 0some positive peak value for l ≤ 0
u(l) → 0 as l → 0
(1.3.6.2)

18
EECQ 4242: Eng. Hydrology 1B

∞ ∞

∫ 𝑢(𝑙 )𝑑𝑙 = 1 𝑎𝑛𝑑 ∫ 𝑢(𝑙 ) 𝑙 𝑑𝑙 = 𝑡𝐿


0 0

The quantity 𝑡𝐿 is the lag time of the IUH. It can be shown that 𝑡𝐿 gives the time interval between the
centroid of an excess rainfall hyetograph and that of the corresponding direct runoff hydrograph. Note
the difference between 𝑡𝐿 and the variable 𝑡𝑝 used for synthetic unit hydrograph lag time − 𝑡𝑝
measures the time from the centroid of the excess rainfall to the peak, not the centroid, of the direct
runoff hydrograph.
The ideal shape of an IUH that of a single-peaked direct-runoff hydrograph, however, an IUH can have
negative and undulating ordinates.
There are several methods to determine an IUH from a given ERH and DRH.
A. For an approximation, the IUH ordinate at time t is simply set equal to the slope at time t of an
S-hydrograph constructed for an excess rainfall intensity of unit depth per unit time. This
procedure is because the S-hydrograph is an integral curve of the IUH i.e., its ordinate at time
t is equal to the integral of the area under the IUH from 0 to t. The IUH so obtained is in general
only an approximation because the slope of an S-hydrograph is difficult to measure accurately.
B. The IUH can be determined by various methods of mathematical inversion, using, for example,
a. orthogonal functions such as Fourier series or Laguerre functions.
b. integral transforms such as the, Laplace transform, the Fourier transform, and the Z
transform
c. mathematical modeling related to watershed geomorphology

19

You might also like