Engineering Hydrology 4E by Wilson
Engineering Hydrology 4E by Wilson
Engineering Hydrology 4E by Wilson
Craig
Soil Mechanics
Van Nostrand Reinhold (International)
Smith
Soil Mechanics
Longman
Tomlinson
Foundation Design and Construction
Longman
Macmillan Education
Macmillan Education Ltd
Houndmills, Basingstoke, Hampshire RG21 2XS
1 Introduction 1
1.1 Allied sciences
1.2 The hydrological cycle
1.3 Inventory of Earth's water 3
1.4 Hydrology as applied in engineering 3
2 Meteorological Data S
2.1 Weather and climate 5
2.2 Humidity 6
2.3 Temperature 8
2.4 Radiation 9
2.5 Wind 9
2.6 Precipitation 10
2.7 Forms of precipitation other than rain 15
2.8 The extension and interpretation of data 15
2.9 The meteorological section of the Flood Studies Report (1975) 27
2.10 Probable maximum precipitation (PMP) 34
References 37
Problems 38
5 Groundwater 83
5.1 The occurrence of ground water 83
5.2 Factors of influence 84
5.3 Groundwater flow 88
5.4 The abstraction of groundwater 94
5.5 The yield of wells 96
5.6 Test pumping analysis 105
References 109
Problems 109
Index 343
This fourth edition, like its three predecessors, is written for engineering students
and junior engineers; to introduce them to the principles and practice of engineer-
ing hydrology and to show, through many worked examples, how to approach
the many apparently intractable problems which hydraulic engineers meet.
The last decade has been a time of considerable activity in the subject, follow-
ing the publication of the Flood Studies Report by the Institute of Hydrology.
Examples of this on-going work include: advances in urban hydrology, published
as The Wallingford Procedure; the Low Flow Studies, the Flood Studies Supple-
mentary Reports and the World Flood Study from IOH; and the Manual for
Estimation of Probable Maximum Precipitation from the World Meteorological
Organisation.
Short descriptions of some of these subjects have been included in this
edition to encourage deeper study of the original texts. The opportunity has
been taken of enlarging the lists of problems, re-organising several chapters,
updating references and including relevant new material.
I continue to receive and much appreciate comments, corrections and advice
from correspondents around the world, to whom I am grateful.
Manchester, 1989' E. M. Wilson
ix
Acknowledgements
x
Acknowledgements xi
From the American Geographical Society, New York
Figure 3.2
From the McGraw-Hi1l Book Company
Figure 4.7
From Dr I. G. Littlewood
Figures 6.13 and 6.14
From the American Water Works Association
Figure 2.18
From the Director, Transport and Road Research Laboratory
Tables 10.1 and 10.2
Figure 10.1
From the Director General, U.N. Food and Agriculture Organisation
Tables 3.2, 3.3 and 3.4
Figure 3.4
xii
1 Introduction
The science of hydrology deals with the occurrence and movement of water on
and over the surface of the Earth. It deals with the various forms of moisture
that occur, and the transformation between the liquid, solid and gaseous states
in the atmosphere and in the surface layers of land masses. It is concerned also
with the sea-the source and store of all the water that activates life on this
planet.
Water in the sea evaporates under solar radiation, and clouds of water vapour
move over land areas. Precipitation occurs as snow, hail, rain and condensate in
the form of dew, over land and sea. Snow and ice on land are water in temporary
storage. Rain falling over land surfaces may be intercepted by vegetation and
evaporate back to the atmosphere. Some of it infiltrates into the soil and moves
down or percolates into the saturated ground zone beneath the water-table, or
phreatic surface. The water in this zone flows slowly through aquifers to river
channels or sometimes directly to the sea. The water that infiltrates also feeds
the surface plant life and some gets drawn up into this vegetation where trans-
piration takes place from leafy plant surfaces.
The water remaining on the surface partially evaporates back to vapour, but
the bulk of it coalesces into streamlets and runs as surface runoff to the river
channels. The river and lake surfaces also evaporate, so still more is removed
he-re. Finally, the remaining water that has not infiltrated or evaporated arrives
back at the sea via the river channels. The groundwater, moving much more
Introduction 3
slowly, either emerges into the stream channels or arrives at the coastline and
seeps into the sea, and the whole cycle starts again.
Fresh-water lakes
Rivers 125125}
062
Soil moisture 65
Groundwater 8250
Saline lakes and inland seas 105 0008
Atmosphere 13 0001
Polar ice-caps, glaciers and snow 29200 21
Seas and oceans 1320000 97'25
Of the 06 per cent of total water that is available as fresh water, about half
is below a depth of 800 m and so is not practically available on the surface. This
means that the stock of the Earth's fresh water that is obtainable one way or
another for man's use is about 4 x 106 km 3 and is mainly in the ground. Spread
over the Earth's land surface it would be about 30 m deep.
The four processes with which the hydrologist is mainly concerned are preci-
pitation, evaporation and transpiration, surface runoff or stream flow, and
groundwater flow. He needs to be able to interpret data about these processes
and to predict from his studies the most likely quantities involved in the extreme
cases of flood and drought. He must be able also to express an opinion about the
likely frequency with which such events will occur, since it is on the frequency
of certain values of extreme events that much hydraulic engineering design is
based.
34
I
32 t-
30
/
ry 28 l-
/:i
T
26
E
E 24 I-
<L " 22
'"
~ 20 I-
I
ew,V ~
)1 ~r-
'"C. 18
;; 16
0
"- 14
:
I---+--+-+-t---+--+---J---J,L--t-,, __ ~ / v ~ 'C? ~~ ~r--
V ~ \ I ~
0
>
c: 12 I
.<;> 1 / ~ 1\ I
0
.20
10
I ",I !?. _ _ -<6> __ ! _-+f-'r~~1p-,.,..-+c...-----1
i
8
(/) I ,...- "'":: .3 (e, t)
6
t-~_
--!~-----tH--+-!.
_ i~"'----t
4
2 F-+---+---+--+-____t D --+--+---+-------i--t---i
oL-~~~_L_L_I~~~-~~L_L_~~_L_L~~I~~I~~~~~I~
-10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16 18 2022 24 26 28 30 32 34
Temperature (DC)
ea (mmHg)
t 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
re)
-10 2.15
-9 2.32 2.30 2.29 2.27 2.26 2.24 2.22 2.21 2.19 2.17
-8 2.51 2.49 2.47 2.45 2.43 2.41 2.40 2.38 2.36 2.34
-7 2.71 2.69 2.67 2.65 2.63 2.61 2.59 2.57 2.55 2.53
-6 2.93 2.91 2.89 2.86 2.84 2.82 2.80 2.77 2.75 2.73
-5 3.16 3.14 3.11 3.09 3.06 3.04 3.01 2.99 2.97 2.95
-4 3.41 3.39 3.37 3.34 3.32 3.29 3.27 3.24 3.22 3.18
-3 3.67 3.64 3.62 3.59 3.57 3.54 3.52 3.49 3.46 3.44
-2 3.97 3.94 3.91 3.88 3.85 3.82 3.79 3.76 3.73 3.70
-1 4.26 4.23 4.20 4.17 4.14 4.11 4.08 4.05 4.03 4.00
-0 4.58 4.55 4.52 4.49 4.46 4.43 4.40 4.36 4.33 4.29
0 4.58 4.62 4.65 4.69 4.71 4.75 4.78 4.82 4.86 4.89
1 4.92 4.96 5.00 5.03 5.07 5.11 5.14 5.18 5.21 5.25
2 5.29 5.33 5.37 5.40 5.44 5.48 5.53 5.57 5.60 5.64
3 5.68 5.72 5.76 5.80 5.84 5.89 5.93 5.97 6.01 6.06
4 6.10 6.14 6.18 6.23 6.27 6.31 6.36 6.40 6.45 6.49
5 6.54 6.58 6.54 6.68 6.72 6.77 6.82 6.86 6.91 6.96
6 7.01 7.06 7.11 7.16 7.20 7.25 7.31 7.36 7.41 7.46
7 7.51 7.56 7.61 7.67 7.72 7.77 7.82 7.88 7.93 7.98
8 8.04 8.10 8.15 8.21 8.26 8.32 8.37 8.43 8.48 8.54
9 8.61 8.67 8.73 8.78 8.84 8.90 8.96 9.02 9.08 9.14
10 9.20 9.26 9.33 9.39 9.46 9.52 9.58 9.65 9.71 9.77
11 9.84 9.90 9.97 10.03 10.10 10.17 10.24 10.31 10.38 10.45
12 10.52 10.58 10.66 10.72 10.79 10.86 10.93 11.00 11.08 11.15
13 11.23 11.30 11.38 11.75 11.53 11.60 11.68 11.76 11.83 11.91
14 11.98 12.06 12.14 12.22 12.96 12.38 12.46 12.54 12.62 12.70
15 12.78 12.86 12.95 13.03 13.11 13.20 13.28 13.37 13.45 13.54
16 13.63 13.71 13.80 13.90 13.99 14.08 14.17 14.26 14.35 14.44
17 14.53 14.62 14.71 14.80 14.90 14.99 15.09 15.17 15.27 15.38
18 15.46 15.56 15.66 15.76 15.96 15.96 16.06 16.16 16.26 16.36
19 16.46 16.57 16.68 16.79 16.90 17.00 17.10 17.21 17.32 17.43
20 17.53 17.64 17.75 17.86 17.97 18.08 18.20 18.31 18.43 18.54
21 18.65 18.77 18.88 19.00 19.11 19.23 19.35 19.46 19.58 19.70
22 19.82 19.94 20.06 20.19 20.31 20.43 20.58 20.69 20.80 20.93
23 21.05 21.19 21.32 21.45 21.58 21.71 21.84 21.97 22.10 22.23
24 22.27 22.50 22.63 22.76 22.91 23.05 23.19 23.31 23.45 23.60
25 23.75 23.90 24.03 24.20 24.35 24.49 24.64 24.79 24.94 25.08
26 25.31 25.45 25.60 25.74 25.89 26.03 26.18 26.32 26.46 26.60
27 26.74 26.90 27.05 27.21 27.37 27.53 27.69 27.85 28.00 28.16
28 28.32 28.49 28.66 28.83 29.00 29.17 29.34 29.51 29.68 29.85
29 30.03 30.20 30.38 30.56 30.74 30.92 31.10 31.28 31.46 31.64
30 31.82 32.00 32.19 32.38 32.57 32.76 32.95 33.14 33.33 33.52
8 Engineering Hydrology
Alternatively, if no change were to take place in the humidity of the air
while it was cooled, then P would move horizontally to the left along line ~
until the saturation line was intersected again. At this point P would be satu-
rated, at a new temperature td, the dew-point. Cooling of the air beyond this
point would result in condensation or mist being formed.
If water is allowed to evaporate freely into the air mass, neither of the above
two possibilities occurs. This is because the evaporation requires heat, which is
withdrawn from the air itself. This heat, called the latent heat of evaporation,
hr, is given by the equation
hr =606.5 - 0.695t calJg
So, as the humidity and vapour pressure rise, the temperature of the air falls and
the point P moves diagonally along line a> until saturation vapour pressure is
reached at the point defined by ew and tw.This temperature tw is called the
wet-bulb temperature and is the temperature to which the original air can be
cooled by evaporating water into it. This is the temperature found by a wet-
bulb thermometer.
The relative humidity is now given as
h = e/e s, or as a percentage, h = 100 e/e s per cent
and is a measure of the air's capacity, at its existing temperature, to absorb
further moisture. It is measured by blowing air over two thermometers, one with
its bulb wrapped in wet muslin and one dry. The air flow past the bulb has an
influence on the wet-bulb reading and the two thermometers can either be
whirled around on a string or more conveniently have the air current provided
by a clockwork fan. In this latter case the instrument is called a psychrometer.
The value of e for air temperature t may be obtained from the equation
2.3 Temperature
Air temperature is recorded by thermometers housed in open louvred boxes,
known as Stevenson screens, about 1.25 m above ground. Protection is necessary
from precipitation and the direct rays of the sun.
Many temperature observations are made using maximum and minimum
thermometers. These record, by indices, the maximum and minimum tem-
peratures experienced since the instrument was last set.
Meteorological Data 9
The daily variation in temperature varies from a minimum around sunrise,
-t
to a maximum from to 3 hours after the sun has reached its zenith, after which
there is a continual fall through the night to sunrise again. Accordingly, maximum
and minimum observations are best made in the period from 8 a.m. to 9 a.m.
after the minimum has occurred.
The mean daily temperature is the average of the maximum and minimum
and is normally within a degree of the true average as continuously recorded.
Temperature is measured in degrees Celsius, commonly, though erroneously,
called centigrade. The Fahrenheit scale is still also in common use.
Distribution of temperature. Generally, the nearer to the equator a place is, the
warmer that place is. The effects of the different specific heats of earth and
water, the patterns of oceanic and atmospheric currents, the seasons of the year,
the topography, vegetation and altitude all tend to vary this general rule, and all
need consideration.
2.4 Radiation
Most meteorological recording stations are equipped with radiometers to measure
both incoming short-wave radiation from sun and sky, and net radiation, which
is the algebraic sum of all incoming radiation and the reflected short-wave and
long-wave radiation from the earth's surface. The net radiation is of great import-
ance in evaporation studies, as will be seen in chapter 3.
2.5 Wind
Wind speed and direction are measured by anemometer and wind vane respec-
tively. The conventional anemometer is the cup anemometer formed by a circlet
of three (sometimes four) cups rotating around a vertical axis. The speed of
rotation measures the wind speed and the total revolutions around the axis gives
10 Engineering Hydrology
a measure of wind run, the distance a particular parcel of air travels in a specified
time.
Because of the frictional effects of the ground or water surface over which
the wind is blowing, it is important to specify in any observation of wind, the
height above ground at which it was taken. An empirical relationship between
wind speed and height has been commonly used
u/uo = (z/zo)O.15
where Uo =wind speed at anemometer at height Zo
u =wind speed at some higher level z.
In recent years there has been some effort to standardise observation heights
and in Europe wind speed is usually observed 2 m above the surface.
Figure 2.2 shows an instrument array for making meteorological observations
at regular, short time-intervals. Instruments, which record automatically on
magnetic tape, include net-radiation radiometer, wet and dry bulb thermometers,
wind vane, anemometer and incident solar radiometer at the mast-top.
2.6 Precipitation
The source of almost all our rainfall is the sea. Evaporation takes place from the
oceans and water vapour is absorbed in the air streams moving across the sea's
surface. The moisture-laden air keeps the water vapour absorbed until it cools to
below dew-point temperature when the vapour is precipitated as rain, or if the
temperature is sufficiently low, as hail or snow.
The cause of the fall in temperature of an air mass may be due to convection,
the warm moist air rising and cooling to form cloud and subsequently to precipi-
tate rain. This is called convective precipitation. This is typified by the late
afternoon thunderstorms that develop from day-long heating of moist air, rising
into towering anvil-shaped clouds. Orographic precipitation results from ocean
air streams passing over land and being deflected upward by coastal mountains,
thus cooling below saturation temperature and spilling moisture. Most orographic
rain is deposited on the windward slopes. The third general classification of rain-
fall is cyclonic and frontal precipitation. When low-pressure areas exist, air tends
to move into them from surrounding areas and in so doing displaces low-pressure
air upward, to cool and precipitate rain. Frontal rain is associated with the
boundaries of air masses where one mass is colder than the other and so intrudes
a cool wedge under it, raising the warm air to form clouds and rain. The slope
of these frontal wedges can be quite flat and so rain areas associated with fronts
may be very large.
Figure 2.3 Automatic recording rain gauge set with rim at ground
level and with anti-splash screen
ment of a catchment's rainfall. Bleasdale [4] quotes tables 2.2 and 2.3 as general
guides and comments as follows
The disparity between the two tables is not so great as might appear at first
sight. The first indicates station densities which are reached in important reser-
voired areas and which may well be exceeded in small experimental areas. The
14 Engineering Hydrology
TABLE 2.2 Minimum numbers of rain gauges required in reservoired moorland
areas a
0.8 2 1 2 3
l.6 4 2 4 6
7.8 20 3 7 10
15.6 41 4 11 15
3l.3 81 5 15 20
46.9 122 6 19 25
62.5 162 8 22 30
Area Number of
rain gauges
mile 2 km 2
10 26 2
100 260 6
500 1300 12
1000 2600 15
2000 5200 20
3000 7800 24
second indicates densities which are more appropriate for country-wide net-
works. In the application of the general guidance embodied in this Table [2.3 J it
must be appreciated that any large river basin will almost invariably have within
it a number of sub-basins for which the relatively dense networks would be
recommended. Moreover, the minimum densities suggested would often be
substantially increased in mountainous areas, and would be closely followed
only in areas of low or moderate elevation without complex topography.
Fog. Estimates of amounts of moisture reaching the ground from fog formation
have been made by installing fog collectors over standard rain gauges. Collectors
consist of wire gauge cylinders on which moisture droplets form and run down
into the rain gauge. Comparisons with standard rain-gauge records at the same
locality show differences that are a measure of fog precipitation. The interpre-
tation of such data requires experience and the use of conversion factors, but
can make substantial differences (of the order of 50 to 100 per cent) to precipi-
tation in forest areas.
Dew. Dew collectors have been used in Sweden and Israel to measure dew fall.
They were made as conical steel funnels, plastic coated and with a projected
plan area of about I square metre. Dew ponds are used as a source of water in
some countries. They are simply shallow depressions in the earth lined with
ceramic tiles.
2.8.1 Defmitions. The total annual amount of rain falling at a point is the usual
basic precipitation figure available. For many purposes, however, this is not
adequate and information may be required under any or all of the following
headings.
Paulhus [5] suggests that if rainfall is plotted against duration, both scales
being logarithmic, the world's greatest recorded rainfalls lie on or just under a
straight line whose equation is
R = 16.6 DO.475
where R is rainfall (in.) and D is duration (h), or
R = 422D.475
where R is in mm.
Meteorological Data 17
Apparently British maxima also lie close to a straight line on a similar plot
with values close to one-quarter of the world's values.
British data are presented in table 2.5 and are plotted in figure 2.4. Maximum
recorded rainfall in the British Isles is closely approximated by
R = 106Do.46
where R is rainfall (mm) and D is time period (h).
, t=' h
,.-
r-
-
- "
-,
_._-
~
1--' _.,
I
I' (9).......
10' . ,
of-d
F.i~o
..I ,'.
,I":-t --
~.
1
I
--
V- I I
I i ; f)
I
.t' ._- ~,
,cc}' 1=1=
,(C) ..
Ilj'
(a)
b : I: I ,
I I i
10
f' ! I I
I1
I I I 1
!
I:
I--- ~
r
f
+
I ~
. ----- ----
--r .- I i
i
1 I I i
10 10 10'
Duration (h)
80
70
60
E
.5
~
50
~
'0 40
C
:J
0
E 30
20
10
P 03Y5
-=1- --=1-012=088
P 56
N
o
..._---..... From a 7-year (1956 to 1962) assessment averaged over 20 stations in England
o 0 From a 35- year (1926 10 1960) assessment averaged over a group of 10 stations
From Silhom's original formula
The rainfall depth on each pair of curves falls in the time marked on the dashed time curlleS at the Micated frequenc
~---,---
400
'" \ \ '..;: \ I
300 I
L,!li' \ r'''~
I':"!t\\\ \ ~ '.
200
~
C
., 150
't
c 100
.Q ~j;~utJ
80
2
'" 60
1 ISi-lE 65 m
:::::I
8. '"0 CS.
:::::I
40 ~ -C----==-LC=- ~ ~T ~ <l> CD
c
~
Ji >- CD
::::!.
:::::I
~ 20 -~ ~-! j5 .Q
"a; cc
~ ::t:
'0
" E ~
10 :J ...
Q; 0
i 8 Il: 0-
c ~~~~~~~~~~~ '~~-~-~,~r--~~~~r-
cc
0 6 -<
""2 4
'"
'0 3 ,., \,~ -=1'''b ,- l
a; ~ I ~ ,
.tl
2
5 ~ - f - - - ---+--1 ", I ' ~..< S", 75
z
'I} 100
10 20 30 40 50 60 70 80 100 110 120 130 140 150
Rain intensity (mm/h)
1
150 6
EI25 f---------+-----t------+~"'---~""-----~ 5
E
I
~IOO
-;e
4
~
~u
c
2 ~eQ(" :5 -
f/
050
75
25 +------+- ---------
0 ~------~--------~------~------~--------~~o
500 1000 1500 2000 2500 3000
Average orrual rainfall (mm)
Hence
102
2 51 .<::
.<:: "E
"c E
,., I
f..
25
-.::---
-;;; 08 20
c
c:'"
06 15
OA --- ~c
0
C -e
03 a:
0
0::
02 --
01 -
0-08 -
-
0-06 - -
004
003
002 L------':--~-:-':-..........J....J....w...'--......J..-..I._-'-..L....J---L...U---....J.....J 05
6 2024
Minutes Hours
Durotion
Another method, due to Thiessen [15] , defines the zone of influence of each
station by drawing lines between pairs of gauges, bisecting the lines with perpen-
diculars, and assuming all the area enclosed within the boundary formed by
these intersecting perpendiculars has had rainfall of the same amount as the
enclosed gauge.
A variation of this technique is to draw the perpendiculars to the lines joining
the gauges at points of median altitude, instead of at mid length. This altitude-
co"ected analysis is sometimes held to be a more logical approach but as a rule
produces little difference in result. Either method is more accurate than that of
the simple arithmetic mean but involves much labour. Thiessen polygons are
illustrated in figure 2.12.
A third method is to draw isohyets, or contours of equal rainfall depth. The
areas between successive isohyets are measured and assigned an average value of
Meteorological Data 23
maXimum depth P
t
E
E
:Ec:
'0
.a:
rainfall. The overall average for the area is thus derived from weighted averages.
This method is possibly the best of the three and has the advantage that the
isohyets may be drawn to take account of local effects like prevailing wind and
uneven topography. A typical isohyetal map is shown in figure 2.13, though the
fall recorded is far from typical, this being the heaviest recorded daily fall in the
United Kingdom.
24 Engineering Hydrology
7 .--
V
6
~
V ~I
V I
5 /
t'
y V
V
4 -- '-----
/
V
V
3
2
! 2 4 6 8 10 20 40 60 80 100 200 300
Ti me t (minutes)
2.8.7 Apparent trends in observed data. From several years' records it may
seem that annual rainfall is, say, declining. It is important to know that this
trend is independent of the gauging, and is due to meteorological conditions
Meteorological Data 25
----.... boundory of
)-..... / catchment or
/ . . . <,
I
I ,
orea of interest
f
l- "\ \
I \
/ \
I \
/ \
( \
I
I
I
I
I
only. This may be checked by plotting a double mass curve as shown in figure
2.15.
A sudden divergence from the straight-line correlation, shown by the dashed
line in the figures, indicates that a change has occurred in gauging and that the
meteorology of the region is probably not the cause of the decline. Such a
change might be due to the erection of a building or fence near the gauge, which
changes the wind pattern round the gauge, the planting of trees, the replacement
of one measuring vessel by another, eVen the changing of an observer to one who
uses different procedures.
2.8.8 Trends from progressive averages. Trends can be more clearly discerned
by the use of the simple statistical technique of examining averages over longer
periods, and moving the group averaged one-year at a time.
Suppose the rainfall records at a station over a number of years are as shown
in figure 2.l6. The first five years on the record are averaged and this average
26 Engineering Hydrology
' \.,
...,
i
I in. '~.,
lin -
2in.
Rain gauges
.c-
t
a.
,::~.o~~,,::::~ 0
c::.. /' -
/
/
-------:1-1
"0 '" I ~
I
/ . d
/' .; deme mass F!J
/' ~ cu've for stotlOn B
/
I
/
/
::;.---- ' /
Time ~
t
/
x /
/
c
o
o
/,
C
~
o
:>
c
c
o
abrupt change indicates
o
change of CIrcumstances
2 in gauge or abserver
<V
~
o
:>
E
:>
u
is plotted at the mid point of the group. The next point is obtained by omitting
the first year and averaging years 2 to 6 inclusive, again plotting the average at
the mid point of the group. In this way the wide variations of particular years
are smoothed out and long-term trends may be detected.
The same techniques can be applied to temperature, hours of sunshine, wind
speeds, cloud cover and other data.
2.9.2 Duration. Here again standards are selected. The standard durations are
2-days and 60-minutes. 2-day MS rainfall is mapped in the FSR for the British
Isles using values from 6000 stations: 60-minute MS rainfall is mapped for the
British Isles as a ratio r = 60-minute MS/2-day MS. Maps of both 2-day MS and
r for almost all of Britain and Ireland are reproduced in appendix A.
2.9.3 Areal reduction factor (ARF). This is the factor that, when applied to
point rainfall of specified duration and return period, gives the areal rainfall for
the same duration and return period. ARF does not vary appreciably with
return period and appears to vary only with area and duration. It was found
to be the same for a wide range of geographical locations through variations in
duration from 1 minute to 25 days, and in areas from I to 30000 km 2 . ARF is
tabulated in table 2.8.
Meteorological Data 29
TABLE 2.6 Growth factors MT/MS for England and Wales (13)
Partial duration
series Annual maximum series
MS
(mm) 2M IM M2 MIO M20 MSO MIOO MIOOO MIOOOO
0.5 0.52 0.67 0.76 1.14 1.30 1.51 1.70 2.52 3.75
2 0.49 0.65 0.74 1.16 1.32 1.53 1.74 2.60 3.94
5 0.45 0.62 0.72 1.18 1.35 1.56 1.79 2.75 4.28
10 0.43 0.61 0.70 1.21 1.41 1.65 1.91 3.09 5.01
15 0.46 0.62 0.70 1.23 1.44 1.70 1.99 3.32 5.54
20 0.50 0.64 0.72 1.23 1.45 1.73 2.03 3.43 5.80
25 0.52 0.66 0.73 1.22 1.43 1.72 2.01 3.37 5.67
30 0.54 0.68 0.75 1.21 1.41 1.70 1.97 3.27 5.41
40 0.56 0.70 0.77 1.18 1.37 1.64 1.89 3.03 4.86
50 0.58 0.72 0.79 1.16 1.33 1.58 1.81 2.81 4.36
75 0.63 0.76 0.81 1.13 1.27 1.47 1.64 2.37 3.43
100 0.64 0.78 0.83 1.12 1.24 1.40 1.54 2.12 2.92
150 0.64 0.78 0.84 1.11 1.21 1.33 1.45 1.90 2.50
200 0.64 0.78 0.84 1.10 1.20 1.30 1.40 1.79 2.30
500 0.65 0.79 0.85 1.09 1.15 1.20 1.27 1.52
1000 0.66 0.80 0.86 1.07 l.12 1.18 1.23 1.42
TABLE 2.7 Growth factors MT/M5 for Scotland and Northern Ireland (13)
Partial duration
series Annual maximum series
MS
(mm) 2M IM M2 MIO M20 M50 MIOO MIOOO MIO 000
0.5 0.55 0.68 0.76 1.14 1.30 1.5 1 1.71 2.54 3.78
2 0.55 0.68 0.76 1.15 1.31 1.54 1.75 2.65 4.01
5 0.54 0.67 0.76 1.16 1.34 1.62 1.86 2.94 4.66
10 0.55 0.68 0.75 1.18 1.38 1.69 1.97 3.25 5.36
15 0.55 0.69 0.75 1.18 1.38 1.70 1.98 3.28 5.44
20 0.56 0.70 0.76 1.18 1.37 1.66 1.93 3.14 5.12
25 0.57 0.71 0.77 1.17 1.36 1.64 1.89 3.03 4.85
30 0.58 0.72 0.78 l.l7 1.35 1.61 1.85 2.92 4.60
40 0.59 0.74 0.79 l.16 1.33 1.56 1.77 2.72 4.16
50 0.60 0.75 0.80 1.15 1.30 1.52 1.72 2.57 3.85
75 0.62 0.77 0.82 1.13 1.26 1.45 1.62 2.31 3.30
100 0.63 0.78 0.83 1.12 1.24 1.40 1.54 2.12 2.92
150 0.64 0.79 0.84 l.l0 1.20 1.33 1.45 1.90 2.50
200 0.65 0.80 0.85 1.09 l.l8 1.30 1.40 1.79 2.30
500 0.66 0.80 0.86 1.08 1.14 1.20 1.27 1.52
1000 0.66 0.80 0.86 1.07 l.l2 1.18 1.23 1.42
30 Engineering Hydrology
TABLE 2.8 Areal reduction factor (ARF) [ 13)
Area A (km 2 )
Duration
D 1 5 10 30 100 300 10003000 10000 30000
2.9.4 The use of the method. To obtain point and areal rainfall for any chosen
location, duration and return period, the procedure is as follows.
(a) Identify the point by its National Grid Reference (NGR). The National Grid
is based on a network of 100 km squares. A diagram showing these and the
letters used to designate them is shown in appendix A. A similar diagram for
the Irish Grid is provided for the Irish section of appendix A. Most maps
used by hydrologists will be 1 :25000 or 1 :50 000 series, which use these
letters. FSR maps use only the 100 km square numbers. To give an NGR,
read the western north-south line number of the square where the point
lies, and estimate or measure tenths and hundredths east of it: followed by
the southern east-west line, and estimate or measure tenths and hundredths
north of it.
(b) Determine the corresponding location on the maps of 2-day M5 and rand
extract values.
(c) Using table 2.9 (which is a model for M5 rainfall for various durations), and
interpolating as necessary for r, the values of M5 at the chosen duration,
expressed as a percentage of 2-day MS, can be found.
(d) With this value of M5 (mm), use the appropriate growth factor (table 2.6 or
2.7) to establish the chosen duration MT factor. The application of this
factor to the M5 value gives the point rainfall of chosen duration and
T-years return period.
(e) Now use the ARF (table 2.8) to establish the areal rainfall for the appro-
priate area around the chosen point.
TABLE 2.9 Model for MS rainfall for durations up to 48 hours [13]
r
MS rainfall (amounts as percentages of 2-day MS)
(per
cent) 1 min 2 min 5 min 10 min 15 min 30 min 60 min 2h 4h 6h 12 h 24 h 48 h
w
....
32 Engineering Hydrology
An example will illustrate the use of the method more clearly.
Example 2.2. Estimate the lOO-year return period rainfall of 6-h duration over
a catchment area of 75 km 2 surrounding the NGR point NM950700.
2.95 Storm profIle. The FSR describes a series of storm proftles of varying
probability for both winter and summer. For specific proposals, several, of
varying peakedness, should be used in the application of 'design rain'. However,
for the purposes of this text, and for use in flood estimation in later chapters,
only the 'winter 75 per cent proftle' will be described. This is the storm proftle
that, on average, is more peaked than 75 per cent of all winter storm proftles;
see figure 2.17 and table 2.1l. Rain is assumed to fall at varying intensity, but
symmetrically about the peak intensity.
TABLE 2.10 Factors to give M5 values for hours from rainfall days [13]
Rainfall days 1 2 4 8
Multiplying factor 1.11 1.06 1.03 1.015
Rainfall hours 24 48 96 192
The profIle is usually used in the form of a stepped profIle as shown in figure
7.30 (page 187). When used in this way it is convenient, though not essential, to
divide the rain duration into an uneven number of parts, to permit symmetrical
arrangement about a peak intensity. (See example 7.2 [15] on page 186.)
34 Engineering Hydrology
2.10 Probable maximum precipitation (PMP)
The preceding section describes methods for estimating rainfall for particular
locations with prescribed frequencies. If one imagines lower and lower frequen-
cies then the rainfall amounts apparently would continue to increase. That,
however, is a statistical concept rather than a physical one and begs the question
- is there some maximum possible value, limited by the nature of the earth's
atmosphere and the laws of physics? This uppermost limit, assuming it exists, is
known as probable maximum precipitation (PMP) and is defined [16] as the
theoretically greatest depth of precipitation for a given duration that is physically
possible over a given size of storm area for a particular location and time of year.
This concept of an upper limit to possible precipitation was known previously
(before 1950) as maximum possible precipitation but PMP is now preferred
because, being an estimate, it has a degree of uncertainty.
--
30
500 600
Mean annual maximum rainfall (mm)
It will be noted that this figure has a maximum Km value of 20, but other
investigators have found even higher values in other regions [21]. Km apparently
varies inversely with intensity so the greater the annual rainfall in a particular
climate the less likely that Km will be greater than 15. The more arid the area,
the greater the likelihood.
Detailed deSCriptions of the methods are given in WMO's Operational Hydro-
logy Report No. 1 [17].
2.10.3 General application of data. Many large storms have been analysed,
particularly in the USA, and maximum values of rainfall depth for various
durations and areas have been published [16,21] . Such data are usually presented
as sets of curves, each representing a rainfall depth plotted on rectangular co-
ordinates of storm duration in hours and storm area in square miles. These are a
useful guide to limiting values but must be used with judgement for particular
catchments, since topography and elevation as well as climate may modify the
result appreciably for other regions.
In the light of this section it is interesting to look again at the world maxi-
mum point rainfalls quoted in section 2.8.2 where it is suggested that if rainfall
is plotted against duration, both scales being logarithmic, the world's greatest
observed point rainfalls lie on or just under a straight line whose equation is
R = 16.6Do. 475
where R is a rainfall in inches andD is duration in hours. Reference to figure 2.4
indicates that the extreme point rain falls recorded in the British Isles, extracted
Meteorological Data 37
from more than half a million station years of data, lie close to a line of equation
R = 106Do. 46
where R is in mm and D in hours. The fit is particularly good in the region from
0.2 to 20 hours, which is probably the critical range for application to unit
hydrographs.
References
1. MAIDENS, A. L. New Meteorological Office rain-gauges. Meteorological
Magazine, 94, No. 1114 (May 1965) 142
2. GOODISON, C. E. and BIRD, L. G. Telephone interrogation of rain-gauges.
Meteorological Magazine, 94, No. 114(May 1965) 144
3. GREEN, M. J. Effects of exposure on the catch of rain gauges. Technical
Publication 67, Water Research Association, July 1969
4. BLEASDALE, A. Rain gauge networks development and design with special
reference to the United Kingdom. International Association of Scientific
Hydrology Symposium on Design of Hydrological Networks, Quebec,
1965
5. PAULHUS, J. L. H. Indian ocean and Taiwan rainfalls set new records.
Monthly Weather Rev., 93 (May 1965) 331
6. BILHAM, E. G. The Classification of Heavy Falls of Rain in Short Periods,
H.M.S.O., London, 1962 (republished)
7. A guide for engineers to the design of storm-sewer systems. Road Research
Laboratory, Road Note 35, H.M.S.O., London, 1963
8. HOLLAND, D. J. Rain intensity-frequency relationships in Britain. British
Rainfall 1961 , H.M.S.O., London, 1967
9. YARNALL, D. L. Rainfall intensity-frequency data. U.S. Department of
Agriculture Miscellaneous Publication, 204, Washington D.C., 1935
10. LINSLEY, R. K. and KOHLER, M. A. Variations in storm rainfall over small
areas. Trans. Am, Geophys. Union, 32 (April 1951) 245
11. HOLLAND, D. J. The Cardington rainfall experiment. Meteorological Maga-
zine, 96, No. I 140 (July 1967) 193-202
12. YOUNG, C. P. Estimated rainfall for drainage calculations. LR 595, Road
Research Laboratory, H.M.S.O., London, 1973
13. Natural Environmental Research Council. Flood Studies Report, Vols. I-V,
NERC, 1975
14. The areal reduction factor in rainfall frequency estimation. FSR Suppl.
Report No. 1, Institute of Hydrology, Wallingford, United Kingdom,
1977
15. THIESSEN, A. H. Precipitation for large areas. Monthly Weather Rev., 39,
(J uly 1911) 1082
16. HANSEN, E. M., SCHREINER, L. C. and MILLER, J. F. Application of
probable maximum precipitation estimates - United States east of the
105th meridian. Hydrometeorological Report No. 52, National Weather
Service, National Oceanic and Atmospheric Administration, US Department
of Commerce, Washington D.C., 1982
17. Manual for estimation of probable maximum precipitation. Operational
Hydrology Report No. 1, 2nd edn, World Meteorological Organisation, 1986
18. HERSHFIELD, D. M. Estimating the probable maximum precipitation. Proc.
Am. Soc. Civ. Eng., Journal Hydraulics Division, 87 (1961) 99
38 Engineering Hydrology
19. HERSHFIELD, D. M. Method for estimating probable maximum precipita-
tion. Journal American Waterworks Association, 57 (1965) 965
20. CHOW, V. T. A general formula for hydrologic frequency analysis. Trans.
Am. Geophys. Union, 32 (1961) 231
21. McKAY, G. A. Statistical estimates of precipitation extremes for the Prairie
Provinces, Canada Department of Agriculture, PFRA Engineering Branch,
1965
Further reading
Problems
2.1 An air mass is at a temperature of 28C with relative humidity of 70 per
cent. Determine: (a) saturation vapour pressure, (b) saturation deficit, (c) actual
vapour pressure in mbar and mm Hg, (d) dew-point, and (e) wet-bulb temperature.
2.2 Discuss the relationships between depth, duration and area of rainfall for
particular storms.
2.3 The following are annual rainfall figures for four stations in Derbyshire.
The average values for Cubley and Biggin School have not been established.
Meteorological Data 39
Average 1959 1960
(in.)
(a) Assume departures from normal are the same for all stations. Forecast
the Rodsley 'annual average' from that at Wirksworth over the two years of
record. Compare the result with the established value.
(b) Forecast annual averages for Cubley and Biggin School using both Wirks
worth and Rodsley data.
(c) Comment on the assumption in part (a). Is it reasonable?
2.4 One of four monthly-read rain gauges on a catchment area develops a fault
in a month when the other three gauges record 37, 43 and 51 mm respectively.
If the average annual precipitation amounts of these three gauges are 726, 752
and 840 mm respectively and of the broken gauge 694 mm, estimate the missing
monthly precipitation at the latter.
2.5 Compute the average annual rainfall, in inches depth, on the catchment
area shown
(i) by arithmetic means, (ii) by the Theissen method, and (ill) by plotting iso-
hyets. Comment on the applicability of each method.
2.6 Discuss the setting of rain gauges on the ground and comment on the effect
of wind and rain falling non-vertically on the catch.
2.7 Annual precipitation at rain gauge X and the average annual precipitation
at twenty surrounding rain gauges are listed in the following table
40 Engineering Hydrology
2.8 Plot the data for the mean of the 20 stations in 2.7 as a time series. Then
plot 5-year moving averages and accumulated annual departures from the 36-year
mean. Is there evidence of cyclicity or particular trends?
(a) What is the value of the exponent relating the two speeds and elevations?
(b) What speed would you predict for the 2 m level?
2.12 What is the average rainfall over an area of 8 km 2 during a storm lasting
30 minutes with a frequency of once in 20 years in (a) Oxford, (b) Kumasi.
Does your answer for (b) require qualification?
2.13 The table below lists the annual maximum rainfall over a 6-hour period,
for an observation station, for 20 years. Make a first estimate of the Probable
Maximum Precipitation of 6-hour duration at the station. Is this station in the
UK?
Depths in mm
146 193
203 124
194 176
102 210
173 102
151 131
182 165
115 187
143 121
180 98
3 Evaporation and Transpiration
Wind. As the water vaporises into the atmosphere, the boundary layer between
earth and air, or sea and air, becomes saturated and this layer must be removed
and continually replaced by drier air if evaporation is to proceed. This move-
42
Evaporation and Transpiration 43
ment of the air in the boundary layer depends on wind and so is a function of
wind speed.
Relative humidity. The third factor affecting evaporation is the relative humi-
dity of the air. As the air's humidity rises, its ability to absorb more water
vapour decreases and the rate of evaporation slows. Replacement of the boundary
layer of saturated air by air of equally high humidity will not maintain the
evaporation rate: this will occur only if the incoming air is drier than the air that
is displaced.
3.2 Transpiration
Growing vegetation of all kinds needs water to sustain life, though different
species have very different needs. Only a small fraction of the water needed by a
plant is retained in the plant structure. Most of it passes through the roots to the
stem or trunk and is transpired into the atmosphere through the leafy part of
the plant.
In field conditions it is practically impossible to differentiate between
evaporation and transpiration if the ground is covered with vegetation. The two
processes are commonly linked together and referred to as evapotranspiration.
The amount of moisture that a land area loses by evapotranspiration depends
primarily on the incidence of precipitation, secondly on the climatic factors of
temperature, humidity etc. and thirdly on the type, manner of cultivation and
extent of vegetation. The amount may be increased, for example, by large trees
whose roots penetrate deeply into the soil, bringing up and transpiring water
that would otherwise be far beyond the influence of surface evaporation.
Transpiration proceeds almost entirely by day under the influence of solar
radiation. At night the pores or stomata of plants close up and very little moisture
leaves the plant surfaces. Evaporation, on the other hand, continues so long as a
heat input is available, although it occurs primarily during the day. The other
factor of importance is the availability of a plentiful water supply. If water is
always available in abundance for the plant to use in transpiration, more will be
used than if at times less is available than could be used. Accordingly, a distinc-
tion must be made between potential evapotranspiration and what actually takes
place. Most of the methods of estimation necessarily assume an abundant water
supply and so give the potential figure.
44 Engineering Hydrology
3.3 Methods of estimating evaporation
(i) by a change in the storage within the catchment, either in surface lakes
and depressions or in underground aquifers;
(ii) by a difference in the underground flow into and out of the catchment;
(iii) by evaporation and transpiration.
If the observations are made over a sufficiently long time the significance of S,
which is not cumulative, will decrease and can be ignored if the starting and
finishing points of the study are chosen to coincide as nearly as possible with
the same seasonal conditions. The significance of U cannot be generalised but in
many cases can be assigned second-order importance because of known geological
conditions that preclude large underground flows. In such cases a good esti-
mation of evapotranspiration becomes possible and the method provides a means
of arriving at first approximations.
3.3.2 Energy budget method. This method, like the water budget approach,
involves solving an equation that lists all the sources and sinks of thermal energy
and leaves evaporation as the only unknown. It involves a great deal of instru-
mentation and is still under active development. It cannot be used readily with-
out many data that are not normally available, and so it is a specialist approach.
3.3.3 Empirical formulae. Many attempts have been made to produce satisfac-
tory formulae for the estimation of evaporation. These are usually for evaporation
from an open water surface, as indeed are the more general methods to follow.
The reason for this is simple. Evaporation, if it is to take place, presupposes a
Evaporation and Transpiration 45
supply of water. Whatever the meteorological conditions may be, if there is
no water present then there can be no evaporation. Accordingly, estimating
methods using meteorological data work on the assumption that abundant water
is available; that is, a free water surface exists. The results obtained therefore are
not necessarily a measure of actual but of potential evaporation. Often these two
are the same, as for example, in reservoirs where a free water surface exists.
When evaporation from land surfaces is concerned, the loss of water in this way
clearly depends on availability: rainfall, water-table level, crop or vegetation,
and soil type all have an influence, which can be expressed by applying an
empirical factor, usually less than unity, to the free water surface evaporation.
There are two cases that should be considered:
(i) when the temperature of the water surface is the same as the air tem-
perature;
(ii) when the air and water surface temperatures are different.
In 1948 Penman [1] presented a theory and formula for the estimation of
evaporation from weather data. The theory is based on two requirements, which
must be met if continuous evaporation is to occur. These are: (i) there must be a
supply of energy to provide latent heat of vaporisation; (ii) there must be some
mechanism for removing the vapour, once produced.
The energy supply. During the hours of daylight there is a certain measurable
amount of short-wave radiation arriving at the earth's surface. The amount
depends on latitude, season of the year, time of day and degree of cloudiness.
Assuming there were no clouds and a perfectly transparent atmosphere, the total
radiation to be expected at a point has been given in tabular form by Angot, and
is reproduced in table 3.1 as values of RA.
If Rc = short-wave radiation actually received at the earth from sun and sky
and n/D = ratio of actual/possible hours of sunshine, then Penman gives (for
southern England)
Rc =R A (0.18 + 0.55n/D)
and quotes Kimball (for Virginia, USA)
Rc =RA (0.22 + 0.54n/D)
and Prescott (for Canberra, Australia)
Rc =R A (0.25 + 0.54n/D)
Thus even on days of complete cloud cover (n/D = 0), about 20 per cent of solar
radiation reaches the earth's surface, while on cloudless days about 75 per cent
of radiation gets throUgh.
Part of Rc is reflected as short-wave radiation; the exact amount depends on
the reflectivity of the surface, or, the reflection coefficient r.
TABLE 3.1 Angot's values of short-wave radiation flux RA at the outer limit of the atmosphere in g cal/cm 2 /day
as a function of the month of the year and the latitude a
Latitude Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. Year
(degrees) m
<I
"Cl
N 90 0 0 55 518 903 1077 944 605 136 0 0 0 3540 0
I
80 0 3 143 518 875 1060 930 600 219 17 0 0 3660 .....
60 86 234 424 687 866 983 892 714 494 258 113 55 4850 er
:::I
358 538 663 847 930 1001 941 843 719 528 397 318 6750 I
40 :::I
599 8070 Q.
20 631 795 821 914 912 947 912 887 856 740 666
Equator 844 963 878 876 803 803 792 820 891 866 873 829 8540 ...-t
I
970 1020 832 737 608 580 588 680 820 892 986 978 8070 :::I
20 III
994 1033 6750 "Cl
40 998 963 686 515 358 308 333 453 648 817 ::;-
I
60 947 802 459 240 95 50 77 187 403 648 920 1013 4850 ..
80 981 649 181 9 0 0 0 0 113 459 917 1094 3660 cS"
:::I
S 90 995 656 92 0 0 0 0 0 30 447 932 1110 3540
aFrom: Physical and Dynamical Meteorology by David Brunt, p. 112 (Cambridge University Press, 1944). The SI unit for
RA is joules/m 2 /day. The table in g cal/cm 2 /day is used so that it is compatible with Rijkoort's nomogram. The conversion is
1 g cal/cm 2 = 41.9 kJ /m 2 .
~
48 Engineering Hydrology
If RI = the net amount of radiation absorbed, then (for southern England)
RI=Rc(1-r)=RA(1-r)(0.18+0.55n/D)
In turn, some of RI is re-radiated by the earth as long-wave radiation, particularly
at night when the air is dry and the sky clear. The net outward flow RB may be
expressed empirically as
RB = aTa4(0.47 - 0.077Ye)(0.20 + 0.80n/D)
where a = Lummer and Pringsheim constant = 117.74 x 10-9 g cal/cm 2 /day
Ta = absolute earth temperature = tOe + 273
e = actual vapour pressure of air in mm mercury.
Hence the net amount of energy finally remaining at a free water surface (r = 0.06)
is given by H, where
H = RI -RB
= Rc -rRc -RB
= Rc(1-r)-RB
= R A (0.18 + 0.55n/D)(1 - 0.06) -RB
Therefore
then
H H
E =-- =
o 1 + {3 t~ - t
1 +')'-,--
es - e
Now eliminate t~ - t by substitution, since t~ - t = (e~ - es )/ t:..
where e s = saturation vapour pressure at temperature t
t:.. = slope of vapour pressure curve at t, = tan a (see figure 3.1).
t eo',
.,
~
'"a.
"
0
a.
0
>
Temperature ~
Since
e~ - es = (e~ - e) - (e s - e) (3.6)
50 Engineering Hydrology
and from equation 3.1
Ea = C(es - e)f(u)
while from equation 3.3
Eo = C(e~ - e)f(u)
then
Ea = es - e (3.7)
Eo e~ - e
where Ea = evaporation (in energy terms) for the hypothetical case of equal
temperatures of air and water.
Then by the values of equations 3.6 and 3.7 into equation 3.5
and
from which
E = Ml + rEa
o t. + r
t. has values obtained from the saturation vapour pressure curve, typically as
shown
t. = 0.36
0.61
1.07
1.80
Referring to equations 3.2 and 3.4 for Ea and H respectively, it can be seen that
Eo is now computed from standard meteorological observations of mean air
temperature, relative humidity, wind velocity at a standard height and hours of
sunshine. The formula has been checked in many parts of the world and gives
very good results. Being based on physical principles it is of general application
and gives values that should serve for most project studies until supplemented by
actual evaporation measurements (see section 3.6).
To overcome the computational labour involved in solving the Penman
equation, a nomogram has been designed by P. J. Rijkoort of the Royal Meteoro-
logical Institute, The Netherlands, which enables rapid evaluations to be made.
Evaporation and Transpiration 51
It is reproduced as appendix C at the back of this book by permission of the
designer. The nomogram has been drawn for a slightly different value of Rc
from that used by Penman
Rc = RA(0.20 + 0.48n/D)
instead of Penman's RA (0.18 + 0.55n/D), but any difference will be smaller than
the probable margin of error in cloud cover estimation, so it can be ignored.
Values of RA can be derived for any latitude from table 3.1.
Where actual sunshine records are available, then table 3.2 may be used to
determine n/D accurately. When sunshine records are not available, n/D may be
estimated from assessment of cloud cover in tenths; i.e. ten-tenths cloud is com-
pletely overcast and 0 tenths completely clear. Unity minus this fraction may be
used for n/D. For example, on a day when about seven-tenths of the sky is
covered by cloud, on average, then
n/D=(l-7/10)=0.3
Several observations per day should be made.
TABLE 3.2 Mean daily maximum hours of sunshine for different months and latitudes a
Latitude North Jan. Feb. March April May June July Aug. Sept. Oct. Nov. Dec.
South July Aug. Sept. Oct. Nov. Dec. Jan. Feb. March April May June
50 8.5 10.1 11.8 13.8 15.4 16.3 15.9 14.5 12.7 10.8 9.1 8.1
48 8.8 10.2 11.8 13.6 15.2 16.0 15.6 14.3 12.6 10.9 9.3 8.3
46 9.1 10.4 11.9 13.5 14.9 15.7 15.4 14.2 12.6 10.9 9.5 8.7
44 9.3 10.5 11.9 13.4 14.7 15.4 15.2 14.0 12.6 11.0 9.7 8.9
42 9.4 10.6 11.9 13.4 14.6 15.2 14.9 13.9 12.5 11.1 9.8 9.1
40 9.6 10.7 11.9 13.3 14.4 15.0 14.7 13.7 12.5 11.2 10.0 9.3
35 10.1 11.0 11.9 13.1 14.0 14.5 14.3 13.5 12.4 11.3 10.3 9.8
30 10.4 11.1 12.0 12.9 13.6 14.0 13.9 13.2 12.4 11.5 10.6 10.2
25 10.7 11.3 12.0 12.7 13.3 13.7 13.5 13.0 12.3 11.6 10.9 10.6
20 11.0 11.5 12.0 12.6 13.1 13.3 13.2 12.8 12.3 11.7 11.2 10.9
15 11.3 11.6 12.0 12.5 12.8 13.0 12.9 12.6 12.2 11.8 11.4 11.2
10 11.6 11.8 12.0 12.3 12.6 12.7 12.6 12.4 12.1 11.8 11.6 11.5
5 11.8 11.9 12.0 12.2 12.3 12.4 12.3 12.3 12.1 12.0 11.9 11.8
0 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1
aFrom: Crop Water Requirements. Irrigation and Drainage Paper 24 (United Nations F.A.O., Rome, 1975).
Evaporation and Transpiration 53
potential evapotranspiration from short, close-set vegetation with an adequate
water supply, in the latitudes of the USA.
If tn = average monthly temperature of the consecutive months of the year in
c (where n = 1,2,3, ... , 12) andj = monthly 'heat index', then
j=(~)1.S14 (3.8)
PE x = 16 (JlOt) a
mm per month
where
a = (675 x 10- 9 )]3 - (771 x 1O- 7 )J 2 + (179 x 10- 4 )1 + 0.492 (3.9)
However PE x is a theoretical standard monthly value based on 30 days and
12 hours of sunshine per day. The actual PE for the particular month with
t
average temperature t C) is given by
DT
PE=PE --mm (3.10)
x 360
where D = number of days in the month
T =average number of hours between sunrise and sunset in the month.
The method has been tested by Serra, who suggested that equations (3.8) and
(3.9) may be simplified as follows
f = 0.09tn 3/2
a = 0.016J + 0.5
This method of estimating potential evapotranspiration is empirical and com-
plicated and requires the use of a nomogram for its solution. Thornthwaite
published such a nomogram, which is reproduced in figure 3.2.
The first step is to obtain the heat indexJ. From figure 3.2 obtain the unad-
justed value of potential evapotranspiration by drawing a straight line from the
location's J value through the point of convergence at t = 26.5 cC. (If t is greater
than 26.5 cC, use the table alongside figure 3.2.)PE x for the month can then be
read off, corresponding to its given mean temperature. Twelve values are obtained
for each of the 12 months. These unadjusted values can then be adjusted for day
54 Engineering Hydrology
40
i III I PE.
t (OC)
30 265 1350
pomt of convergence-
-
u
~
25
20 ;:;.
270
275
280
139'5
1437
1478
t=
15
V"'" 28'5 1517
~ 290 155'4
e..
::J
: ;\m V 295
300
1589
1621
--
10 / 30'5 1652
V
Cl.
E 90
80 310 1680
~
-
V
31-5 1707
,.. 70
320 173'1
~ 60
325 1753
C
0 5 0
330 1772
E ~ 335 1790
.. 3 /
c
4 340 1805
0
30 345 1818
:::;: 350 1829
355 1837
2 - 20 360
365
184'3
1847
370 1849
375 1850
380 1850
10 _I I
10 20 30 40 60 80 100 150 20 0
50 70
PE. (mm)
and month length by equation 3.10 and totalled to give annual potential evapo-
trl\llspiration.
It has been found that the method gives reasonably good results whatever the
vegetation cover, though different types of vegetation will affect a particular
locality's true value. The formula is based on temperature, which does not
necessarily correspond to incoming solar radiation immediately, because of the
'heat inertia' of land and water. Transpiration, however, responds directly to
solar radiation. Accordingly, care should be exercised when using the method
to ensure that conditions do not change abruptly in a particular month, though if
figures for many consecutive months are being used, the cumulative differences
are probably negligible.
For project studies the method is a useful complement to the Penman
approach, though it will be found generally to give higher values of potential
evaporation than the latter.
3.7.1 Arable crops. Consumptive use refers to water that is actually used, while
evapotranspiration formulae give potential water use. Reference to figure 3.3
will show that for a particular locality, unless the rain falls when it is needed, a
large water deficit may develop in the growing season despite quite high rainfall.
In a case like this, consumptive use will be less than potential evapotranspiration
and a need for irrigation in the growing period is indicated.
Modern methods of estimation of consumptive use include a development of
the Penman open-water evaporation calculation and a more empirical method
developed from the Blaney-Criddle procedures of the Division of Irrigation and
8~------------------------,
rate of potential
evopotronspiration in
in./rnonth
6~------------~~~-------;
soil moisture
utilisation
-co
preClpilotlon
4~-------+----~---t------~ (in/month)
Q
o'"
soil moisture
accretion
TABLE 3.3 Typical crop coefficient kc for crops at different growth stages and
prevailing climatic conditions a
11
3
10
~ti~midity < 20%
9 20-50%
------>50%
8
1 light daytime wind 0-2 m/s
2 Moderate daytime wind 2-5 m/s
7 3 Strong daytime wind> 5 m/s
>
'"
"0
6
E
E
t-? 5
lLJ
2 3 4 5 6 7 8
Blaney-Criddle factor f ~ P (0.46t (QC) + 8.13) mm/day
Figure 3.4 ETo against f for medium cloudiness (n/D: 0.6-0.8), 3 wind con-
ditions and 3 relative humidities. After Crop Water Requirements: Irr. and
Drainage Paper 24 (United Nations F.A. 0., Rome, 1975)
60 Engineering Hydrology
Table 3.4 gives the mean daily percentage (p) of annual daylight hours for
different latitudes.
The ETo and kc values for both the modified Penman and Blaney-Criddle
methods are compatible, so only one set of kc values is required.
(1) With annual rainfall around 2300 mm, there was a considerably greater
loss from the forested catchment than from the grassland. The forested
catchment losses were about double those of the grassland - about
850 mm against 405 mm on average.
(2) The reason for the heavier forest loss is the interception of rainfall by the
forest canopy and the subsequent evaporation from it.
(3) In drier parts of the country, where the rainfall was about 600 mm
annually, forested catchments were no more likely to lead to losses than
grassland, and in some observations of pine forest were actually less [14].
(4) The Penman formula for evapotranspiration is in good agreement with
the actual figure for grassland catchments but underestimates the figure
considerably for forested catchments in wetter parts of the country and
during wet periods in areas of low rainfall.
(5) Although erosion protection has always been an aim of afforestation, it
is not always achieved. Where open ditching has been adopted to drain
forest areas, the erosion may be greater than for grassland. This was the
case in the Plynlimon catchments, where sediment was measured over
five years.
(6) From observations of lysimeters, it has been found that about 50 to 55
per cent of rainfall at the forest canopy falls through to the forest floor,
while a further 12 to 23 per cent reaches ground level by stem-flow; so
that between 22 and 38 per cent evaporates from the canopy.
(7) Mature forest cover reduces peak flows by about one-third compared
with grassland, but where drainage ditches have been dug for tree planting
there may well be both an increase compared with grassland and shorter
response times, particularly while the trees are immature.
TABLE 3.4 Mean daily percentage (p) of annual daytime hours for different latitudes [10]
Latitude North Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec.
South July Aug. Sept. Oct. Nov. Dec. Jan. Feb. Mar. Apr. May June
60 0.15 0.20 0.26 0.32 0.38 0.41 0.40 0.34 0.28 0.22 0.17 0.13
58 0.16 0.21 0.26 0.32 0.37 0.40 0.39 0.34 0.28 0.23 0.18 0.15
56 0.17 0.21 0.26 0.32 0.36 0.39 0.38 0.33 0.28 0.23 0.18 0.16
m
54 0.18 0.22 0.26 0.31 0.36 0.38 0.37 0.33 0.28 0.23 0.19 0.17 <
11)
"Cl
52 0.19 0.22 0.27 0.31 0.35 0.37 0.36 0.33 0.28 0.24 0.20 0.17 0
...
11)
50 0.19 0.23 0.27 0.31 0.34 0.36 0.35 0.32 0.28 0.24 0.20 0.18 er
48 0.36 0.28 0.24 0.21 0.19 11)
0.20 0.23 0.27 0.31 0.34 0.35 0.32
-=
46 0.20 0.23 0.27 0.30 0.34 0.35 0.34 0.32 0.28 0.24 0.21 0.20 c.
=
44 0.21 0.24 0.27 0.30 0.33 0.35 0.34 0.31 0.28 0.25 0.22 0.20 ...-I
11)
42 0.21 0.24 0.27 0.30 0.33 0.34 0.33 0.31 0.28 0.25 0.22 0.21 en
=
"Cl
:::;"
11)
40 0.22 0.24 0.27 0.30 0.32 0.34 0.33 0.31 0.28 0.25 0.22 0.21
35 0.23 0.25 0.27 0.29 0.31 0.32 0.32 0.30 0.28 0.25 0.23 0.22 o
-=
30 0.24 0.25 0.27 0.29 0.31 0.32 0.31 0.30 0.28 0.26 0.24 0.23
25 0.24 0.26 0.27 0.29 0.30 0.31 0.31 0.29 0.28 0.26 0.25 0.24
20 0.25 0.26 0.27 0.28 0.29 0.30 0.30 0.29 0.28 0.26 0.25 0.25
15 0.26 0.26 0.27 0.28 0.29 0.29 0.29 0.28 0.28 0.27 0.26 0.25
10 0.26 0.27 0.27 0.28 0.28 0.29 0.29 0.28 0.28 0.27 0.26 0.26
5 0.27 0.27 0.27 0.28 0.28 0.28 0.28 0.28 0.28 0.27 0.27 0.27
0 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27
0)
....
62 Engineering Hydrology
Water engineers concerned with land use in reservoir catchments should be aware
that if annual rainfall is higher than average then afforestation will increase
water loss, whatever its other benefits.
References
1. PENMAN, H. L. Natural evaporation from open water, bare soil and grass.
Froc. Roy. Soc., A 193 (April 1948) 120
2. PENMAN, H. L. Evaporation over the British Isles. Quart. J. Roy. Met. Soc.,
76 (1950) 372
3. THORNTHWAITE, C. W. An approach towards a rational classification of
climate. Geographical Rev., 38 (1948) 55
4. British Rainfall 1939 (and subsequent years), H.M.S.O., London
5. LAW, F. The aims of the catchment studies at Stocks Reservoir, Slaidburn,
Yorkshire. Unpublished communication to Pennines Hydrological Group,
Institution of Civil Engineers, September 1970
6. HOUK, I. E. Irrigation Engineering, Vol. 1, Wiley, New York, 1951
7. OLlVIER, H. Irrigation and Climate, Arnold, London, 1961
8. THORNTHWAITE, C. W. The moisture factor in climate. Trans. Am. Soc. Civ.
Eng., 27, No. 1 (February 1946) 41
9. BLANEY, H. F. and CRIDDLE, W. D. Determining water requirements in
irrigated areas from climatological and irrigation data. Div. Irr. Water
Conserv., S.C.S. U.S. Dept. Agr., SCS-TP-96, Washington D.C., 1950
10. DOORENBOS, J. and PRUITT, W. O. Crop water requirements. Irrigation and
Drainage Paper 24. F.A.a., United Nations, Rome, 1975
11. LAW, F. The effect of afforestation upon the yield of water catchment
areas. J. Brit. Waterworks Assoc., 38 (1956) 484
12. Water balance of the headwater catchments of the Wye and Severn 1970-75.
Report No. 33. Institute of Hydrology, Wallingford, United Kingdom,
December 1976
13. NEWSON, M. D. The results of ten years' experimental study on Plynlimon,
Mid-Wales, and their importance for the water industry. J. Inst. Water
Eng. Sci., 33 (1979) 321-33
14. GASH, J. H. C. and STEWART, J. B. The evaporation from Thetford Forest
during 1975. J. Hydrol., 35 (1977) 385-96
Further reading
Problems
3.1 Determine the evaporation from a free water surface using the Penman
equation nomogram for the following cases
3.2 Use the nomogram for the solution of Penman's equation to predict the
daily potential evapotranspiration from a field crop at latitude 400 N in April,
under the following conditions: mean air temperature = 20C; mean h = 70
per cent; sky cover =60 per cent cloud; mean U2 = 2.5 mts; ratio of potential
evapotranspiration to potential evaporation =0.7.
A B Daylight A B Daylight
hours of hours of
year at A year at A
(per cent) (per cent)
Jan. -5 -2 6 July 19 16 11
Feb. 0 2 7 Aug 17 14 10
Mar. 5 3 71 Sept. 13 10 S-!-
Apr. 9 7 S;; Oct. 9 S 7;;
May 13 10 10 Nov. 5 3 7
June 17 15 11 Dec. 0 0 6
(a) at A for April (mean temperature 10C) and November (mean tem-
perature 3C),
(b) at B for June (mean temperature 20C) and October (mean temperature
SoC).
At A the average number of hours between sunrise and sunset is 13 for April
and 9 for November. At B the figures are 14 for June and 10 for October. Use
the Serra simplification for A and the nomogram for B.
3.4 Water has maximum density at 4C; above and below this temperature its
density is less. Consider a deep lake in a region where the air temperature falls
below 4 c in winter.
(a) Describe what will happen in the lake in spring and autumn.
(b) What will be the effect of what happens on (i) the time lag between air
and water temperatures? (ii) the evaporation rate in the various seasons?
(c) Will there be a difference if the winter temperature does not drop below
4 C and if so explain why.
3.5 Using the modified Blaney-Criddle method, determine the crop water
requirement for growing cotton in mid-season in very dry conditions at a location
at latitude 30 0 N in February where the mean daily temperature is 28C, and the
mean wind speed is 4 m/so
3.6 Discuss the advantages and disadvantages of evaporation pans placed above
the ground surface (for example, the U.S. Class A pan) compared to those sunk
in the earth.
3.9 A large reservoir is located in latitude 400 30N. Compute monthly and
annual lake evaporation for the reservoir from the given data using the nomo-
gram of Penman's theory. If the Class A pan evaporation at the reservoir for the
year is 1143 mm, compute the pan coefficient. Assuming that the precipitation
on the lake is as given and that the runoff represents unavoidable spillage of this
precipitation during floods, what is the net annual anticipated loss from the
reservoir per square kilometre of surface in cubic metres per day?
What would the change in evaporation be for the month of July if the reser-
voir was at 40 0 S?
.s:: .s::
"- "-
E E
E E
'----
Low rain intensity
Time ~ Time!~ 8 16 24 32
Slope (per cent)
Soil type fo fe K
(mm/h) (mm/h) (min- l )
surface soil grains are shown, then the governing factor is the head h over the
smallest cross-section of a pore. This continues to increase with rainfall intensity
until a limiting value is reached where runoff prevents any further increase. It
seems unlikely that this limiting condition is often reached in natural conditions.
Previous researchers [3] have found similar results but attributed the increase
in fc at higher rainfall intensities to lack of homogeneity in their experimental
catchment watershed. Others have also emphasised the overriding importance of
the superficia//ayer of a soil [4] .
The infiltration rate of a soil is the sum of percolation and water entering
storage above the groundwater table. Generally the soil is far from saturated and
so storage goes on increasing for very long periods. Accordingly, fc goes on
decreasing under a steady rain intensity for equally long periods.
However the soil will eventually become saturated under persistent rain. All
the reservoirs of interstitial space become filled. When precipitation ceases the
soil gradually loses the 'free' water, to a point where it can sustain the water
content against gravity, i.e. water will no longer drain from it. At this point the
soil is said to be at field capacity.
Exposed soils can be rendered almost impermeable by the compacting impact
of large drops coupled with the tendency to wash very fine particles into the
voids. The surface tends to become 'puddled' and the fc value drops sharply.
Similarly, compaction due to man or animals treading the surface, or to vehicular
traffic, can severely reduce infiltration capacity.
Dense vegetal cover such as grass or forest tends to promote high values of
f c . The dense root systems, all providing ingress to the subsoil, the layer of
organic debris forming a sponge-like surface, burrowing animals and insects
opening up ways into the soil, the cover preventing compaction and the vege-
tation's transpiration removing soil moisture, all tend to help the infiltration
process.
Infiltration and Percolation 69
Other effects that marginally affect the issue are frost heave, leaching out of
soluble salts and drying cracks which increase to, and the entrapping of inter-
stitital air, which decreases to. Temperature has some effect since flow in inter-
stices is laminar and hence viscosity has a direct effect on resistance to flow.
Other things being equal, {o and {e will have higher values during the warmer
seasons of the year.
5 mm over the surface and continually refilled to maintain this depth, the inflow
to the central tube being measured. The purpose of the outer tube is to eliminate
to some extent the edge effect of the surrounding drier soil. Such tests give
useful comparative results but they do not simulate real conditions and have
been largely replaced by sprinkler tests on larger areas. Here the sprinkler simu-
lates rainfall, and runoff from the plot is collected and measured as well as
inflow, the difference being assumed to have infiltrated.
While rain-simulating sprinklers are a good deal more realistic than flooded
rings, there are limitations to the reliability of results thus obtained, which
usually give higher values of {than natural conditions do. For qualitative effects
(for example, comparisons between different conditions of vegetation, soil
types etc) the methods are simple and effective.
Consistent and repeatable results can be obtained by using laboratory catch-
ments with rainfall simulators, where the quantity and thickness of soil is
adequately representative of nature. Nassif and Wilson [1] used 7 tonnes of soil
in a layer 200 mm thick, and measured all inputs, outputs and changes of
storage. Such equipment as theirs gives very good comparative and perhaps
absolute figures for infiltration, but it still does not simulate natural conditions
completely as there is atmospheric pressure at the bottom of the laboratory soil
layer, while this is not so in nature.
70 Engineering Hydrology
4.3.2 Drainage basin rainfall-runoff analysis. Several investigators have attempted
to improve on sprinkled infiltrometers by choosing small 'homogeneous' drain
age basins and carefully measuring precipitation, evaporation and outflow as
surface runoff. By eliminating everything except inflitration, average f values
can be obtained for such basins by techniques presented by Horton [5] and
Sherman [6] .
The difficulty remains of ensuring that there has been neither unrecorded
underground outflow nor variation in underground aquifer storage, so that
although quantitative results are obtained the analysis is intricate and the margin
for error is wide.
nel rain =
quantity of
runoff
storm loss
Time (h) -+
storm is shown plotted on a time base in terms of its average hourly intensity.
The shaded area above the dashed line represents measured runoff, as mm, over
the catchment area. Since the un shaded area below the line is measured rainfall
but did not appear as runoff, it represents all the losses, induding surface deten-
tion and evaporation as well as inflitration. However, inflitration is much the
largest loss in many catchments and although it is a rough and ready method
(since it takes no account of the variation in f with time) it is widely used as a
means of quickly assessing probable runoff from large catchments for particular
rainstorms.
Suppose, however, the same total rainfall had been distributed as shown in
figure 4.5(b). To obtain a runoff of 33 mm above the I>-index line requires the
~ .<:
..... 25 .....
E E
.5 E
z:. ~
.~
c 'i!
J! J!
.5 .5
line to be raised to give an index value of 9 mm/h. It is seen therefore that one
determination of the I>-index is of limited value and that many such deter-
minations should be made, and averaged, before the index is used.
4.3.4 The fay method. This method is a developed version of the I>-index in
that it attempts to allow for depression storage and short rainless periods during
a storm, as well as eliminating all rain periods where the rainfall intensity is less
than the infiltration capacity assumed.
Referring to figure 4.6, the approximate position of the fay line is fixed by
reference to runoff data and the raingraph. The line is then moved vertically
until the various losses are balanced and the runoff values satisfied. The loss
..
.;;; deducled for infillrotion
c depression
.5
Time (h) -+
If t is unity, then any day's value is k times that of the previous day. If precipi-
tation occurs it will increase the value of the index by an amount that is inde-
terminate, since some rain may have left the catchment as surface runoff. The
Infiltration and Percolation 73
amount added to the index should therefore strictly be the basin recharge only,
but the difference in the index by using all the precipitation is usually small.
The progressive daily reduction of the index is due to evapotranspiration,
which alters seasonally, so equation 4.1 is used with a k value that also varies
seasonally. When the index is used to assess the runoff that takes place from a
particular rainstorm, this variation can be incorporated into a graphical coaxial
relationship derived from the analysis of a larger number of observed rainfall and
runoff data on a particular catchment. Unsley et af. [8] give a detailed descrip-
tion of how such relationships can be derived, and figure 4.7 is reproduced to
illustrate the technique.
The diagram is entered at the antecedent precipitation index and a horizontal
line is followed until the particular week-number curve corresponding to the
calendar date is met. From the intersection a line is dropped vertically down to
intersect with the line of appropriate storm duration in hours, and then followed
II- - - i - - - + - - - t - - ' - - + - -- - t- - - I I
~---t---+----t~~-+----+---t----~
{1 .
4.4.4 Catchment wetness index (CWI). The property of catchment wetness has
an important influence on the quantity of net or excess rainfall that provides
runoff and in subsequent chapters is used in the estimation of particular flood
events by the FSR methods. It is convenient, therefore, to define the index here.
A period of 5 days preceding a storm event is held to provide a recent history
of rainfall and a decay index of 0.5 is proposed in the formula for API:
API5 = 0.5 -.!.[
2 Pd-l + 0.5 P d - 2 + (0.5) 2 Pd-3 + (0.5)3 Pd-4 + (0.5)4 Pd-SJJ
Here P denotes daily rainfall and the suffix indicates the relevant day. From this
CWI = 125 + API5 - SMD
For the purposes of this text, it is sufficient to know how to derive average
values of CWI for use in flood-prediction calculations where the 5 days preceding
rainfall is unknown. Such average values have been calculated in the FSR and are
strongly dependent on annual average rainfall. The relationship for the British
Isles is shown in figure 4.9.
4.4.5 Soil classification. A further and perhaps most important property of the
catchment in respect of runoff is the soil cover, its depth, permeability and
slope.
The Soil Surveys of England and Wales, Scotland, Ireland and the Northern
Ireland Ministry of Agriculture have prepared soil maps of the British Isles,
which have been used as a basis for calculating a soil index as a term in many
of the equations recommended in the Flood Studies Report.
Infiltration and Percolation 77
The basis for comparison is 'the winter rain acceptance potential', which is
the reverse of the runoff potential. It is influenced by permeability, position of
the groundwater level and the slope of the terrain. Five classes of winter rain
acceptance are indicated on the soil maps:
Class Winter rain Runoff
acceptance
Very high Very low
2 High Low
3 Moderate Moderate
4 Low High
5 Very low Very high
A soil index is calculated from the proportions of a catchment's area that are.
covered by each class. The precise definition of this index, SOIL, is given in
section 9.4.1.
The guidelines for each classification are as set out in table 4.2, which can be
used for preliminary estimates of ungauged catchments when soil maps are not
available.
Soil classification maps (entitled RP, denoting runoff potential) for the
British Isles are given in appendix A.
TABLE 4.2 Classification of soils by runoff potential (source: Flood Studies Reportt
a Upland and peaty soils are Class 5; urban areas unclassified. Class l, very low runoff; Class 5, very high runoff.
Infiltration and Percolation 79
If a catchment is provided with a number of permanently installed aluminium
tubes at suitable locations, then a single battery-powered instrument can be
carried from site to site to measure the in situ soil moisture at each one, thus
removing much of the guess-work from the evaluation of this parameter.
The neutron probe is normally used for making repeated measurements of
moisture change at the same site and depth and for such a purpose can give a
very high precision. The accuracy of absolute moisture values on the other hand
is heavily dependent on accurate calibration at each site and depth and this is
generally impracticable.
A diagram of the instrument and a photograph of one in use are shown in
figures 4.10 and 4.11 respectively.
Depth Indicator
Cable Clamp- __~*~
LOCk----J... r
.,
0 Q
.,
0
0
..--Access Tube
0
S 0
YN'"t"n So",,,
Figure 4.10 Wallingford soil-moisture probe (courtesy of the Institute
of Hydrology)
80 Engineering Hydrology
References
1. NASSIF, S. H. and WILSON, E. M. The influence of slope and rain intensity
on runoff and infiltration. Bull. Int. Assoc. Sci. Hydrol., 20, No. 4 (1976).
2. HORTON, R. E. The role of infiltration in the hydrologic cycle. Trans. Am.
Geophys. Union, 14, ((1933) 443-60
Infiltration and Percolation 81
3. BOUCHARDEAU, A. and RODIER, J. Nouvelle Methode de determination de
la capacite d'absorption en terrains permeables. La Houille Blanche, No.
A (July/August 1960) 531-6
4. SOR, K. and BERTRAND, A. R. Effects of rainfall energy onthe permeability
of soils. Proc. Am. Soc. Soil Sci., 26, No. 3 (1962)
5. HO,RTON, R. E. Determination of infiltration capacity for large drainage
basins. Trans. Am. Geophys. Union, 18, (1937) 371
6. SHERMAN, L. K. Comparison of F-curves derived by the methods of Sharp
and Holtan and of Sherman and Mayer. Trans. Am. Geophys. Union, 24,
No. 2 (1943) 465
7. BUTLER, S. S. Engineering Hydrology, Prentice-Hall, Englewood Cliffs,
New Jersey, 1957
8. LINSLEY, R. K., KOHLER, M. A. and PAULHUS, J. L. H. Hydrology for
Engineers, McGraw-Hill, New York, 1958, p. 162
9. Estimated soil moisture deficit over Gt. Britain. Explanatory Notes Meteoro-
logical Office, Bracknell (issued twice monthly)
Further reading
BELL, J. P. Neutron probe practice. Report No. 19, Institute of Hydrology,
Wallingford, United Kingdom
GRINDLEY, J. Estimation of soil moisture deficits. Meteorological Magazine,
96, (1967) 97
GRINDLEY, J. Estimation and mapping of evaporation. Int. Assoc. Sci. Hydrol.
Symposium on World Water Balance, Reading. lASH Publication 92, 1970,
pp. 200-213
HORTON, R. E. Analyses of runoff-plot experiments with varying infiltration
capacity. Trans. Am. Geophys. Union, Part IV, (1939) 693
PENMAN, H. L. The dependence of transpiration on weather and soil condi-
tions,l. Soil Sci., 1, (1949) 74
WILM, H. G. Methods for the measurements of infiltration. Trans. Am. Geophys.
Union, Part Ill, (1941) 678
Problems
4.1 Discuss the influence of slope of catchment and rainfall intensity on
infiltration rates under constant rainfall.
4.3 The table below gives the hourly rainfall of three storms that gave rise to
runoff equivalent of 14, 23 and 18.5 mm respectively
82 Engineering Hydrology
2 4 3
2 6 9 8
3 7 15 11
4 10 12 4
5 5 5 12
6 4 3
7 4
8 2
Determine the <P-index for the catchment.
4.4 Why is the method of subtracting infIltration rates from rainfall intensities
to compute hydrographs of runoff not applicable to large natural river basins?
4.6 Use the co-axial relationship of figure 4.7 to determine how the runoff in
this river changes seasonally. Assume that during week number 1 a storm of
5 in. of rain lasting 72 h occurs. Compare what happens with the effects of the
same storm in week number 25, if the antecedent precipitation index in each
case is 1.5 in. Suggest which seasons of the year the weeks are in and explain
why there should be a difference in runoff.
4.7 (a) Describe 'antecedent precipitation index' and 'soil moisture deficit',
indicating how they are derived and used.
(b) List the parameters necessary for the solution of Penman's equation
for open-water evaporation. Describe how you would obtain suitable
values for each of them.
4.8 Calculate the Catchment Wetness Index (CWI) for a point 400 km north and
400 km east of the UK National Grid origin, if the preceding 5 days' rainfall
were (in mm)
12.0,0.5,4.2,0.0,3.5
Compare this with the FSR design values.
5 Groundwater
83
84 Engineering Hydrology
flowing well, and a fracture or flaw in the impermeable overlay will, in such
conditions, result in an artesian spring. Sometimes a small area of impermeable
material exists in a large aquifer. This happens through geological faulting or, for
example, through a lens of clay occurring in otherwise sandy glacial drift. A
small local water-table, called a perched water-table may result and this can
often be a long way above the true phreatic surface.
Some of the modes of occurrence of groundwater described above are illu-
strated in figure 5.1.
/ /
/ / / /
/
/
i rr,
Pe
rrr,
e () b
/ e
on
~ 10000
.;:
o
.c
U
1030
Oensity p (kg 1m 3 )
Viscosity is a measure of the shear strength of a liquid; the lower the viscosity,
the more mobile or penetrating the liquid.
Absolute viscosity, denoted by Il, has units of N s/m2 ; the CGS unit, the poise =
10- 1 N s/m 2 Water at 20C has a viscosity of 1 centipoise (0.01 poise = 10-3
Ns/m 2 ).
Porosity is defined as n =total VOids/total volume and ranges from a few percent
to about 90 per cent. In a granular mass composed of perfect uniform spheres:
in the loosest possible packing, n = 47.6 per cent
in the densest possible packing, n = 26 per cent
Natural soils are, of course, composed of irregular particles of many different
sizes. The more regular the soil, the more porous it tends to be, since in irregular
soils the smaller particles tend to fill the voids in the larger particle packing. It
86 Engineering Hydrology
is therefore, standard procedure in any groundwater survey to analyse the soils
mechanically and plot the particle sizing in a standard way, using a logarithmic
size scale. A typical analysis is shown in figure 5.3. Two soils are plotted: the
more regular soil has the steeper slope and is likely to be more porous.
100
,/"
/'
80
/'
<;
c
.,
.;:
60 /V
0'
0
I I Ilv
.,
C
40
I
.,~
Q.
/,/1 I
/
.--
,// Ij
20
.---I-'
001 0-1 1-0 10
Groin size (mm)
B_ s_----r----::'"'iI':-------,----:-:so=-=n-:;-d----r-----
cloy I---T:-::---r=m:-:-e"i
d ,-um:::-r-=-co:-:a-,s:-:e-t-.f'"n:-:e--'-=m::::ed7-,u"-:m""c-=-oo-,-se-l g r a vel
0-002 0-006 002 0-06 0-2 0-6 2-0mm-
plastiC
,iI' fine
sond
medium coarse
rOve I
When water fills the pores of a soil there is a thin layer, only a few molecules
thick, that coats the particles. This water is not free to move and adheres to the
particles even when the voids have been drained, occupying part of the available
space. This means that the effective porosity, ne, may be less than the true
porosity n. In coarse materials such as gravels there will be no difference between
ne and n but in fine sands the difference may be 5 per cent, or even more in very
fine materials. In most considerations of porosity in the flow of groundwater it
is ne, the effective porosity, that is of importance.
10~
grovels
10
~ fme
103
gravel
~ ~ond
102
t-
>- coarse
0
"0
..... 10'
.5 fine sond
c:
-'C 0
100
C
~ very
r
0-
.'!' fine sond u
u
;;: .,'"
~""
0; 10"
0
u
"
~ 10" g
:0
0
Cl>
~ ,"od, clay
Cl>
a
10-5
pure clay,
10 6
5.3.1 Darcy's Law. Before any mathematical treatment of the flow of ground-
water can be attempted, it is necessary to make certain simplifying assumptions
The fundamental law is Darcy's Law of 1856. This states that the rate of flow
per unit area of an aquifer is proportional to the gradient of the potential head
measured in the direction of flow. Or
The specific velocity is not the true velocity, but is merely Q/A. The actual
velocity of water in the pores is greater than the specific velocity since the path
the water follows through porous media is always longer than a straight line
beween any two points.
If the average real, or pore, velocity is denoted by v then
v = ___d_is_ch_a_r-=.ge___ = ~ = Av =~
area of water passage Ane Ane ne
Hence
pore velocity (average) = specific velocity -;- effective porosity
Since the velocity distribution across a pore is probably parabolic, being
highest in the centre and zero at the edges, the maximum pore velocity vrnax =
2 x v (approximately).
Groundwater 89
So in a typical case where (say) ne =1, then
v=3v and vmax=6v
While these are typical figures only, it is important to remember the order of the
v
velocities, since it is on max that the Reynolds number and the continuance
of laminar flow depends.
5.3.2 Flow in a confined aquifer. Consider now the case of unidirectional flow
in a confined aquifer of permeability k, illustrated in figure 5.5. Groundwater is
--- --
L ~..L.LJ'~~'LL..~
potentia I
piezometers head 4>
4> ImpervIous
I low .. . conllnl ng
H 1.loye?
f 7 /1 ; J } ; , , 1/ ) J ;; i ;; ,,),} J) I ) J J j
flowing from left to right, the energy required to move the water through the
pores is continually using up the available head and so the line of potential head
as indicated by piezometers introduced into the aquifer, is declining. From
Darcy's Law
v = -k drp
x dx
dq =0
dx
and differentiating the equation for q, above
dq d2 rp
-=-kH-
dx dx 2
90 Engineering Hydrology
So
(5.2)
5.3.3 Flow in an aquifer with phreatic surface. Consider now the case of an
aquifer with a phreatic surface, resting on an impermeable base. This case is
illustrated in figure 5.6. Here the first equation from Darcy's Law would be
v = - kd</l
-
S dl
--_qdX
equ"'polenlia:'-
dl
;h=dCP
__
/ - -.......
linps (Ph real ic
sur fnce
assumed impervious
. I ; ~~r!>z~njt~l/b~:e,; , ; ~) ) 77 I
, ; f , 1 ) ; if' 1 f i I I I J J7
dh
q == -kH- (5.3)
dx
Groundwater 91
and
dq _ d 2 (h 2 )*
dx - -tk dx2
-- - ---
net infiltration N
q q+dq
dx
dq = -tk d2 (h 2 ) =N
dx dx 2
and hence
(5.5)
Example 5.1. Suppose there are two canals, at different levels, separated by a
strip of land 1000 m wide, of permeability k =12 m/day as shown in figure 5.8.
If one canal is 2 m higher than the other and the depth of the aquifer is 20 m
below the lower canal to an impermeable base, find the inflow into, or abstrac-
tion from, each canal per metre length. Take annual rainfall as 1.20 m per
annum and assume 60 per cent infiltration.
Assume a reference origin as indicated on figure 5.8. Then the boundary condi-
tions are simply: when x = 0, h = 20; and when x = 1000, h = 22.
N = 60 per cent of 1.2 = 0.72 m/year
= 0.72/365 m/day
------... --t-T-:_::-_-=_~_=-::-
2m
k = 12 m/day
h 20m
IOOOm
h
N 2
h 2 = - -x + CIX + C2 (5.6)
k
When x = 0, h = 20, therefore
400 =0 + 0 + C2
and
C2 =400
and
Cl = 0.084 + 0.164 = 0.248
Now from equation 5.3
q =-kh-
dh
dx
dh =
- 1
-I (N-.2x + 0.248 )
-
dx 2112 k
Atx=O
1
q = -ku'i. - I (0.248) = -6(0.248) = -1.49 m3 JdayJm
I
2112
Atx = 1000
=_ ~ (_ 2000 x 0.72 + 0.248\
q 2 365 x 12 ')
The construction of drilled wells, well screens and gravel packs and the tech-
niques of well development and maintenance are beyond the scope of this book.
Readers requiring information on these topics should consult references [1] and
[2] .
Once the water has entered the well it has to be pumped to the surface. Well
pumps are classified as reciprocating, rotating vertical shaft, jet and air lift
pumps. Rotating vertical shaft are either surface driven or submersible, and can
be centrifugal or rotary positive displacement.
By far the most common application is now the electric submersible centri-
fugal borehole pump, with the electric drive motor directly coupled to the pump
stages in one long pump body that is placed near the bottom of the well. Such
pumps are manufactured in sizes down to 100 mm diameter to supply heads up
to 100 m or more if necessary. Such a 100 m (4 in.) diameter pump would
supply about 4 m 3 /h while a 250 mm (10 in.) diameter pump might supply
30 times as much. A sketch of a typical installation is shown in figure 5.10.
96 Engineering Hydrology
--- - -- --:. - --
=
..
~~ment grout _
.... . ~ ".", \:I ~v " ....
-
- ~ ... . . - :: . "
. !.: ... ' -.:
Aqu i fer
grove l pock
~motor
SI - S2 - . In-
= -Qo r2 (Th" .)
elm s equatIOn (5.7)
2rrkH rl
This is a very important equation and holds true (at least approximately) for all
kinds of underground flow, steady and unsteady, confined and unconfined. The
value of kH is known as the coefficient of transmissibility.
Equation 5.7, indefinitely integrated, yields
S=- ~.lnr+C
2rrkH
and if s = 0 when r =R then
s= ~ln!i (5.8)
2rrkH r
Either of the two equations 5.7 and 5.8 enables the drawndown curve to be
established, provided that the integration constants Qo and R can be determined
from the boundary conditions. Qo is the constant discharge of the pumped well
and so can be measured. R, however, varies from one observation point to
another. In most cases, however, it is the drawdown close to the well that is
important, where R has the value Ro. Some values of Ro for idealised boundaries
are given in figure 5.12, where the shaded areas are land and the blank adjacent
spaces are open water.
This figure illustrates how Ro depends on the distance to open water. As the
distance increases, so does Ro and the drawdown. Indeed finite drawdowns are
possible only when the constant-level open water is a fmite distance away. Also,
98 Engineering Hydrology
I ground surface
. . . . . . . . . . -.... l
... r ~. ~;e'z~~elf:~ -:face
. before pumping
-
H
flow
r r
5.5.2 Steady unconfined flow. When the drawdown is slight compared with the
thickness of the aquifer, the factor kH remains nearly constant and the formula
for steady confined flow may be used. As drawdown increases, the falling water
level reduces the area of transmitting aquifer and the equations now become, in
the notation of figure 5.13: from Darcy's Law
dh
Q = 21Tr . h . k -
dr
and from continuity
Q = Qo = constant
Combining these
h. dh = Qo . dr
21Tk r
/ ground surface
- . .. - - - - . - - -
.....,...,....,.-
" ..
/ / / 'I
r ~----:----. r
-11~
Figure 5.13 Well pumping from an un confined aquifer
100 Engineering Hydrology
Integrating
h 2 = Qo In r + C
1Tk
and if
h =Hat r=R
then
H2 _h 2 = QOln~ (5.9)
1Tk r
The value of R must satisfy the boundary conditions. Then the draw down at the
wellface (H - ho) is deduced by introducing the radius of the well ro since
Example 5.3. A well is drilled to the impermeable base in the centre of a circular
island of 1 mile diameter in a large lake. The well completely penetrates a sand-
stone aquifer 50 ft. thick overlain by impermeable clay. The sandstone has a
permeability of 50 ft./day. What will be the steady discharge if the drawdown of
the piezometric surface is not to exceed 10ft. at the well, which has a diameter
of 1 ft?
For a well in the centre of the circular island, the boundary condition is
s = 0 when r = 2640 ft.
/00 /
piezometric surface
before pumping
--'
K= 50 fUday -
5280'
10 = Qo In 2640
21r x 50 x 50 0.5
Therefore
= 0.212 ft. 3 /s
5.5.3 Steady unconfmed flow with rainfall. When rainfall is present the equa-
tions become: from Darcy's Law
dh
Q = 2nrh . k -
dr
and from continuity
dQ = -2nr . dr. N (where N =net inflltration)
Integrating
Q = -nr 2N + Cl
and Cl may be determined from the condition that when r = ro :.::: 0, Q = Qo.
Therefore
Q = -nr2N + Qo
Substituting this value in the Darcy equation
h dh = Qo . dr _ N . r dr
. 2nk r 2k
and integrating
(5.10)
-----
..... --::-:-:.~.-:- ~
... ' k
. : .2
':":
'
. . . . h
H
r
. .
1
(5.11)
H2 _ h~ = _ N (R2 _ ,2)
2k
l 1
Area within which
1 phreatic surface with
rainfall but no pumping
rainfall supplies \
well output )
x l- . - " - - - phreatic surface with
r--f--H'-+---.... both pumping
and rainfall
-----
H= IOm
k = 20 m/day
Well diameter 06 m
-I--ro m
=
/
03
500 m - - - "
I
Figure 5 .16 Circular island with central well, rainfall and an unconfined
aquifer. Solution by superposition
H2 _ h~ = Qo In ~
1Tk r
2 500
10O - h2 = 21.9 log - (5.13)
r
Inserting r = 0.3 m (at the wellface) in equations 5.12 and 5.13 yields h 1 =
11.18 m and h2 = 5.43 m. Since h =h 1 + h2 - H (by superposition of draw-
downs), it follows that
h = 11.18 + 5.43 - 10
=6.61 m
Therefore
Sw = hl - h
=4.57 m
If a water-divide exists, then all the output of the well is being contributed by
rainfall, since if the sea around the island were contributing, the hydraulic
104 Engineering Hydrology
gradient would be sloping downwards and inwards at every point. So the area
contributing is obtained from
Qo = 1Trx 2 N
25 =1Tr 2 0.004
x--
x 24
Therefore
rx = 218 m (well within the 500 m radius of the land)
Using this value of r, equation 5.12 yields hI = 10.97 and equation 5.13 yields
h2 = 9.60. Then
h=hI+h2-H
= 10.97 + 9.60 - 10.0
= 10.57
Sx = 10.97 - 10.57
=0.40m
Suppose the simple formula of equation 5.8 had been used, thus assuming
constant aquifer thickness. In this case, without rainfall
-
Sw---
Qo 1n-
R -- 600 x 2 .3 1og--
500
21TkH r 21T x 20 x 10 0.3
= 0.477 x 2.3 x 3.223
= 3.54m (compare the value of 4.57 m above)
and
500
Sx = 0.477 x 2.3 x log -
218
= 1.01 x 0.36
= 0.364 (compare the value of 0.40 m above)
It will be realised that the simple formula for the confined case is adequate in
this case except in the immediate neighbourhood of the well. It would, of course,
still be necessary to compute the 'no-pumping' phreatic surface, for the case
with rainfall.
Groundwater 105
5.6 Test pumping analysis
5.6.1 General. Using equation 5.8, the drawdown resulting from groundwater
abstraction from wells can be determined if the boundary conditions and geo-
hydrological constants are known. In particular, the presence of water-bearing
strata that form aquifers, their extent, thickness and permeability are usually
unknown until a number of test borings have been made. Each test boring, as
well as providing information on the underlying geology of an area, may be left
open with a simple porous screen in the bottom of the hole to allow subsequent
observation of water levels.
Although much information can be gleaned from the drilling of these test
bore holes/observation wells, the extent of the yield of an aquifer can finally be
determined only by test pumping from a well. Such a well should be positioned
so that observation wells are placed on either side of it, on a line through the
well, and preferably on two lines at right angles with the pumping well at the
intersection.
Pumping tests can be performed either by pumping from the well at a steady
rate until steady-state conditions are obtained, (that is, there is no appreciable
movement with time in any of the observation wells) and then plotting the data
recorded, or by observation of the rate of change of level in all the wells up to
and including steady-state conditions. Variants include stepped pumping tests
where the discharge is progressively increased at regular intervals [4]. Only the
first of these methods is discussed here. For full treatment of groundwater
recovery, including test pumping refer to Huisman [3], Verruijt [5], and other
specialised texts.
Observation wells should ideally be spaced at increasing intervals from the
pumping well-say at 20, 50, 100, 200 and 500 m-depending on the depth and
expected productivity of the aquifer. Always use the largest capacity pump
available to do the actual pumping. Constant-rate extraction may have to con-
tinue for days, or even weeks and months in some cases before steady-state
conditions are reached. Careful observation of all the wells should, therefore,
be made before starting pumping, at regular intervals throughout the test and
during recovery of the levels after pumping has ceased, until the initial equili-
brium levels are regained.
s=~ln~
2rrkH r
where R depends on both boundary conditions and the point of observation. If
the aquifer is of large extent and r (the radial distance to an observation well) is
not excessive then for r < O.lR
s=~lnRo (5.14)
2rrkH r
106 Engineering Hydrology
where Ro is the integration constant for the well-face, and for all points of
observation has the same value.
If the observed values of drawdown s are now plotted against distances r on
logarithmic paper, then a straight-line relationship is obtained:
kH= ~ 1.15Qo
B 1T
and
A
logRo =-
B
14
12
10 ~
~~
E 08
I"-
I
0-6
04
~
'I '-..
02
~
"t'----
o
10 20 304050 100 200300 500 1000
rIm)
H2 _h 2 = Qo In~
rrk r
The left-hand side of this equation can be written as (H - h) (H + h). But since
(H -- h) = s and hence h = (H - s), then the left-hand side can also be written
as s(2H - s). Therefore
s Qo In R
rrk(2H - s) r
or
s Qo In Ro
(5.15)
2rrk(H - f) r
s' =s (1 __2Hs_) = ~
2rrkH
In Ro
r
and s' can be plotted against r in the manner of figure 5.17.
108 Engineering Hydrology
Example 5.6 Figure 5.18 is a plot of observed drawdown against r for a well
in an unconfined aquifer discharging a constant 0.03 m 3 /s. The depth of the
aquifer below the phreatic surface has been established during drilling as about
20 m. Determine the formation constants.
Obse"ed
drawdown
60
50
~
40
~~
E I " "::::: ::::~
v; 30 Corrected
drowdown
s' ~
20
~
10 ~
t-..r--,....
o
10 20 304050 100 200300 500 1000
rim)
s', the corrected drawdown, is computed for each observed drawdown and from
the straightline relationship with r
s' =A +B logr
As s' = 4.8 when r = 10, hence log r = 1.00
and as s' = 0 when r = 830, hence log r = 2.92
Therefore
A - 2.92B = 0
A - l.OOB = 4.8
Hence
A =7.30 andB =2.50
Also
so
Ro = 830m
Groundwater 109
References
1. CRUSE, K. A review of water-well drilling methods. J. Eng. Geol., 12 (1979)
63
2. STONER, R. F. et al. Economic design of wells. J. Eng. Geol., 12 (1979) 63
3. HUISMAN, L. Groundwater Recovery, Macmillan, London, 1972
4. BRERETON, N. R. Step-drawdown pumping tests for the determination of
aquifer and bore-hole characteristics. Water Resources Council Tech.
Report 103, Washington D.C., January 1979
5. VERRUIT, A. Theory of Groundwater Flow, 2nd edition, Macmillan, London,
1981
Further reading
ARONOVICI, V. S. The mechanical analysis as an index of subsoil permeability.
Proc. A m. Soc. Soil Sci., 11 (1947) 137
CEDERGREEN, H. R. Seepage, Drainage and Flow Nets, John Wiley, New York,
1967
CHILDS, E. C. and COLLIS.GEORGE, N. The permeability of porous materials.
Proc. Roy. Soc., A201 (1950) 392
KIRKHAM, D. Measurement of the hydraulic conductivity of soil in place.
Symposium on Permeability of Soils. American Society for Testing and
Materials, Special Tech. Publication 163, 1955, p. 80
RUSHTON, K. R. and REDSHAW, S. C. Seepage and Groundwater Flow, Wiley,
1979
TODD, D. K. Ground Water Hydrology, John Wiley, New York, 2nd edition,
1980
WENZEL, L. K. Methods for determining permeability of water bearing materials.
U.S. Geological Survey Water Supply Paper 887, Washington D.C., 1942
Problems
5.1 Rainfall of 2.50 m per annum falls on a strip of land 1 km wide lying
between two parallel canals, one of which (canal A) is 3 m higher than the other
(canal B). The infIltration rate is 80 per cent of the rainfall and there is no run-
off. The aquifer that contains the canals is 10 m deep below the level of canal B
and both canals penetrate it fully. It is underlain by a horizontal impermeable
stratum. Compute the discharge per metre length into both canals, assuming
their boundaries are vertical and the aquifer coefficient of permeability
K = 10 m/day.
5.2 A fully -penetrating well of 0.5 m diameter in a confined aquifer has been
test-pumped at a rate of 0.025 m 3 /s until steady-state conditions have been
reached.
Observation wells at various distances from the well s,~ow the following
results:
111
112 Engineering Hydrology
A second group of questions arises in using the curves of runoff frequency
and duration once found; for example
(i) how can the volume of discharge be reduced?
(ii) how can the cost of flood prevention be compared with the damage that
will arise if no measures are taken?
(iii) how valuable is stored flood water in times of drought?
These questions are not directly related and each involves a quite different
and distinct approach, though the same techniques can be used in answering
more than one. In this and the following sections ways in which some of these
questions can be answered will be sought.
(a) Catchment area. The area as defined at the beginning of this section is
usually, but not necessarily, bounded by the topographic water-divide. Figure
6.1 shows a hypothetical cross-section through the topographic water-divide
of a catchment. Because of the underlying geology it is perfectly possible for
areas beyond the divide to contribute to the catchment. The true boundary is
indeterminate, however, because although some of the groundwater on the left
of the divide in the figure may arrive in catchment B, the surface runoff will stay
in catchment A. Here the infiltration capacity of the soil and the intensity of the
rainfall will influence the portion of the rainfall that each catchment will collect.
If runoff is expressed not as a total quantity for a catchment but as a quantity
per unit area (usually m 3 /s per square kilometre or ft. 3 /s per square mile), it is
observed, other things being equal, that peak runoff decreases as the catchment
area increases. This is due to the finite time taken by the water to flow through
the stream channels to the control section (the time of concentration) and also
to the lower average intensity of rainfall as storm size increases (see section
2.8.4.). Similarly, minimum runoff per unit area is increased because of the
greater areal extent of the groundwater aquifers and minor local rainfall.
Surface Runoff 113
topographical
/ rain
/
water divide
catchment catchment
A B
(b) Main stream length (MSL). This is measured in km from the gauging station
or catchment outfall. When measuring from maps it is usual to use a standard
routine to remove subjectivity (for example, set dividers at 0.1 km on a 1:25000
scale map).
(c) Slope of catchment. The more steeply the ground surface is sloping the
more rapidly will surface runoff travel, so that concentration times will be
shorter and flood peaks higher. Infiltration capacities tend to be lower as slopes
get steeper, since vegetation is less dense and soil more easily eroded, thus
accentuating runoff.
Slope can be enumerated by covering a catchment contour map with a
rectilinear grid and evaluating the slope, perpendicular to the contour lines at
each grid point as shown in figure 6.2a. A frequency distribution of these
numbers can then be plotted as in figure 6.2b. Different catchments can be
compared on the same plot, the relatively steep frequency curves indicating
catchments of fast runoff and flat curves the converse.
A simpler method is to describe the slope in m/km between two points on the
main stream. The standard points used in Britain are 10 per cent and 85 per cent
of the main stream length above the point of interest. Slope defined in this way
is written as S1085.
Figure 6.2 (aJ Rectilinear grid to evaluate catchment slope and orientation. (b) Frequency curves
for comparison of catchment steepness
Surface Runoff 115
meridian (say) at each grid point and the subsequent plotting of a circular
frequency diagram like that of figure 6.3, similar to a wind rose.
N
A B le
:;/ ) I ~ VI
I
at
.... 1---
...-
./ ---
t ....-
f..-/ Il
/
~ /J II
."
/
a2 /
VV I
~~V 1-'"
/
o 5 10 15
Time (h)
(f)Annual average rainfall. The standard annual average rainfall (SAAR) for the
period 1941-1970 for the British Isles is given in the ten section maps for Britain
and four for Ireland in appendix A.
(h) Baseflow index (BFl). This is an index calculated as the ratio of the flow
under the separated hydrograph (see section 7.3) to the flow under the total
hydrograph. It is a good indicator of the catchment's underground storage,
which is dependent on the solid geology. Methods of estimation are given in
section 6.6.4 and figure 6.5 illustrates how BFI is calculated from flow data.
Surface Runoff 117
HYDROGRAPH
?;
o
..J
IL
..
>-
..J
o ". .....
.
:z
w
::E /
/
/
/ ' ' .....
//BASEFLOW
/
-~
TIME
BASEFLOW INDEX
(i) Lake and reservoir area. These act as surface water stores and have the
effect of smoothing out the hydrographs of catchments that contain them.
'0
c
er"
Time ~
slope: 1 vertical
in y horizontal
E
E
et
P(mm)----
R = P--x (6.1)
y
so that the annual rainfall can be used to obtain a first approximation to the
annual runoff. Young has analysed the correlation on a world-wide basis [1].
Variations from the straight line may be due to conditions in the preceding
year that gave markedly higher or lower groundwater levels, or to variations in
the seasonal distribution ofrainfall. The method can be used also for wet months
in humid climates when the ground is saturated, but beyond such narrow limits
it is not valid.
Although the application of a relationship like equation 6.1 is restricted, it
can nevertheless be a useful method for estimating total annual runoff or com-
pletely ungauged catchments, if they are in similar climates and of similar size
and character. The relationship has been developed for use with synthetic hydro-
graphs, using catchment wetness index and soil characteristics as parameters
[2] .
Further refinements are possible taking into account the particular period of
the year, the antecedent precipitation index (see section 4.4.2) and the storm
duration as well as depth, so that relationships may be derived for particular
storms on a particular catchment. Coaxial graphs may be produced which take
the various variables into account. A relationship of this form is shown in
figure 4.7.
6.5.1 Definition. A rating curve is a graph that shows the connection between
the water level elevation, or stage of a river channel at a certain cross-section,
and the corresponding discharge at that section. A typical rating curve is shown
in figure 6.8. It can be seen that the curve is drawn through a cluster of points,
each of which represents the results of a river discharge measurement. Such
measurements can be made in a number of ways, of which the most important
are
(l) velOcity-area methods,
(2) flow-measuring structures,
(3) dilution gauging.
6.5.2 Velocity-area methods. These are conventional for medium to large rivers
and involve the use of a cu"ent meter, which is a fluid velocity-measuring
instrument. A small propeller rotates about a horizontal shaft, which is kept
Surface Runoff 121
Discharge 0 --+
Figure 6.9 A modern helix current meter by HUger and Watts Ltd
Surface Runoff 123
from a manned cable-car or directly from the cable, the instrument in this latter
case being positioned by auxiliary cables from the river banks.
Depths should always be measured at the time of velocity observation since a
profile can change appreciably during flood discharges. Observers should also
remember such elementary rules as to observe the stage before and after the
discharge measurement, and to observe the water slope by accurate levelling to
pegs at the river level as far upstream and downstream of the gauging site as is
practicable, up to say 500 m in each direction.
As water velocities increase in high floods the ballasted current meter will be
increasingly swept downstream on an inclined cable. The position of a meter in
these circumstances can be found reasonably accurately if the cable angle is
measured. Ballast can be increased but only within limits. Rods can be used to
suspend the meter but a rigid structure in the boat will then be required to
handle the rods, calling for a stable platform of the catamaran type. Rod vibra-
tion and bending are common in deep rivers unless diameters exceed 50 mm, by
which time the whole apparatus is getting very heavy and unmanageable.
It will be appreciated that since each river is unique, each will require a
careful assessment of its width, depth, likely flood velocities, cable support
facilities, availability of bridges, boats, etc. before a discharge measurement
programme is started.
From many observations on many rivers it has been established that the
variation of velocity integrated over the full depth of a river can be closely
approximated by the mean of two observations made at 0.2 and 0.8 of the depth.
If time and circumstances preclude even two observations at each horizontal
intercept then one reading at 0.6 depth measured from the surface will approxi-
mate the average over the whole depth.
The discharge at the cross-section is best obtained by plotting each velocity
observation on a cross-section of the gauging site with an exaggerated vertical
scale. Isovels or contours of equal velocity are then drawn and the included
areas measured by planimeter. A typical cross-section, so treated, is shown in
figure 6.lOa. Alternatively the river may be subdivided vertically into sections
and the mean velocity of each s~ction applied to its area (figure 6.10b). In this
method, the cross-sectional area of anyone section should not exceed 10 per
cent of the total cross-sectional area.
A check should always be made using the slope-area method of section
6.5.7 (Hi) and a value obtained for Manning's n. In this way a knowledge of the n
values of the river at various stages will be built up, which may prove most
valuable in extending the discharge rating curve subsequently.
To ensure uniformity in the techniques of current-meter gauging the Inter-
national Organisation for Standardisation (ISO) has published various recom-
mendations and, in the United Kingdom, BS 3680 'Measurement of liquid flow
in open channels' refers [3]. In the USA, the U.S. Geological Survey and the
U.S. Bureau of Reclamation have established practice [4-7].
124 Engineering Hydrology
(c)
(b)
The method involves the injection of a chemical into the stream and the
sampling of the water some distance downstream after complete mixing of the
chemical in the water has occurred. The chemical can either be added by
constant-rate injection until the sampling downstream reveals a constant concen-
tration level, or administered in a single dose as quickly as possible, known as
gulp injection. In this case samples over a period of time disclose the concen-
tration-time correlation. In both cases the concentration of chemical in the
samples is used to compute the dilution, and hence the discharge of the stream
can be obtained. Figure 6.12 shows constant-rate injection of sodium dichromate
from a Mariotte bottle (a constant-head device) in a mountain stream.
Analysis of the samples is by an automated colorimetric procedure that
estimates the concentration of very small amounts of the chromium compound
by comparison with a sample of the injection solution. The equipment is expen-
sive and specialised. References [16, 17] give comprehensive guidance.
Another method, developed by Littlewood [18], deserves mention at greater
length as the required equipment is simple and relatively cheap. The method
depends on the electrical conductivity of solutions of common salt (NaCl) in the
stream water and is a version of the relative-dilution-gauging method of Aastad
and S~gnen [19].
126 Engineering Hydrology
1
Conceltt aljoll
Hence
(6.2)
and
V -
Q= - x N (6.3)
Tz
where N is the mean dilution ratio.
It is not necessary to measure the shaded area under the curve of figure 6.13
since, provided the curve is defined by closely spaced ordinates, the average of
these (in the shaded area) is all that is needed. This holds good for any length of
128 Engineering Hydrology
injection time, including a 'gulp', and not just for the constant-rate injection
of figure 6.13 which is simply a special case of a general method.
It is necessary now to move from solution concentration to conductivity but
two problems must be overcome:
(a) the relationship between concentration of salt and conductivity of the
solution is non-linear
(b) natural streams have varying background conductivities.
The first of these is overcome by constructing a graph of conductivity against
concentration which for very weak solutions is virtually linear. This is so for
solvents of initially different conductivities, where the gradients of the lines in
figure 6.14 are practically parallel.
It follows that, for changes of concentration, .:lc, and conductivity, .:lc', in
the linear range of weak solutions
.:lc = K I.:lC'
Kl for NaCl in water is approximately O.5ll1s cm- 1/mg 1-1.
Equation 6.2 may now be written
Q= ~ x Cl
Tz Kl.:lc'
where .:lc' is the average of (Cl - Co). Now from equation 6.3
N= __C-=-l__ = _C_l_
(cz - co) Kl.:lc'
which can be written
,
- Cl
N=Kz= (6.4)
.:lc'
where Kz is a combination of K 1 and a multiplier for converting a particular
injection solution conductivity, c; , to concentration, Cl .
If a small volume, v, of the particular strong solution is added to a larger
volume V* of the streamwater, giving rise to a dilution ratio N*, and the change
in conductivity is recorded as .:lc' *, then we can write
,
N*=K ~ (6.5)
z .:lc' *
whereN* = -V*
v
Now dividing equation 6.4 by equation 6.5 gives
N .:lc' *
N*
Surface Runoff 129
250
200
'E0
(/)
::i.
~
s 100
+'
0
:l
"U
C
.0
u
M
o 20 40 60 BO 100
-1
Concentration of added salt (mg I I
Figure 6.14 Conductivity against concentration of added salt (weak solution) [18]
Q=-
v X -
V*
X
!:le' *
Tz V !:le'
This means that the discharge may be measured by pouring a known volume
(V) of a strong salt-solution into the stream, measuring the change in conduc-
tivity of the water at the downstream end of the mixing length (to find !:le') over
time T z and measuring the change in conductivity of the weak solution obtained
by controlled dilution of the strong solution in a separate operation.
The method is independent of the make and calibration of the conductivity
meter. It is not necessary to know the concentration of the strong solution and
very simple equipment is needed.
Full details of the method, comparative test results and typical quantities and
dimensions are given in reference [18] .
130 Engineering Hydrology
6.5.5 Ultrasonic gauging. Among recent developments in river gauging is the
method of ultrasonics. It can provide continuous unattended discharge measure-
ment or, more conventionally, at some specified regular interval. The method
depends on the transmission of ultrasonic pulses between two sets of trans-
mitter/receivers situated on either bank of the river, and so arranged so that the
line of transmission is at 45, or thereabouts, to the flow direction. Since the
speed of transmission in one direction is greater than in the other, owing to the
component of the water flow velocity along the line of transmission, the pulse
travel times are different. This difference is a function of the mean water velocity
at the pulse level.
The usual arrangement consists of a section of a river with a rectangular lined
cross-section at least as long as its width, with pairs of transmitters at various
levels and a water level recorder. Automatic samples and recording can then
provide a comprehensive discharge record at the station to an accuracy at least as
good as any other method.
For detailed descriptions of installations, methods and analyses of results,
interested readers should consult references [20,21].
The methods of sections 6.5.2, 6.5.3, 6.5.4 and 6.5.5, used singly or in con-
junction, will establish the correlation for any stream or river discharge with
stage.
A rating curve when established in this way, enables a single observation of
stage, made each day at a set time, by an unskilled observer, to be converted irtto
a discharge rate and hence into a finite quantity of water, flowing at the observa-
tion point. The difficulty about rating curves is to obtain enough points at times
of high discharge to enable an accurate correlation to be obtained.
Example 6.1. Suppose that the gauge shows a rise at the rate of 0.2 m/h during
a discharge measurement of 100 m 3 /s and the channel is such that this rate of
rise may be assumed to apply to a 1000 m reach of river between the measure-
ment site and the reach control (the control of a river reach is the section at
which the profile of the river changes).
Let the average width of channel in the reach be 100 m. Then the rate of change
of storage in the reach, dS is given by
Surface Runoff 131
dS = 1000 x 100 x 0.2
= 20000 m 3 /h
= 5.6 m 3 /s.
The discharge measurement should be plotted on the rating curve as 94.4 m 3 /s
(not as 100), since this is the discharge past the control corresponding to the
mean gauge height.
The second reason for the looping of rating curves is the variation in surface
slope that occurs as a flood wave moves along the channel. Figure 6.15 represents
1/-
-----
htdh
The figure represents an advancing flood wave and rising state. The gauge
reading at the commencement of the discharge measurement is h and at its
conclusion, dt, later, h + dh. In this time the wave has advanced Udt, and
AS = dh =dh/dt (6.8)
Udt U
dh/dt being positive for a rising stage, negative for a falling stage. Combining
equations 6.6, 6.7 and 6.8
or
Qa = / ( 1 + dh/dt) (6.9)
Q", US
Surface Runoff 133
If the discharge measurement is taken at a site with two gauges, one upstream
and one downstream, in the same reach, then all the terms in equation 6.9 are
measured except Q and U, Q being the steady-state discharge required and U the
flood wave celerity. There are now several ways to proceed. The first is to take
an empirical value for U. Corbett et al. [4] suggest that in a fairly uniform
channel, in flood conditions, the celerity U is approximately equal to 1 .3 times
the mean water velocity, or
U= 1.3 Qa_
A
where A = cross-sectional area.
From which
J(
Q= Qa (6.10)
l + A . dh/dt)
1.3 QaS
Example 6.2 A river discharge measurement made during a flood indicated
Qa = 3160 m 3 /s. During the measurement, which took 2 h, the gauge height
increased from 50.40 to 50.52 m. Level readings taken at water surface 400 m
upstream and 300 m downstream of the observation site differed by 100 mm.
The river was 500 m wide with an average depth of 4 m at the time of measure-
ment. At what co-ordinate should the measurement be plotted on the rating
curve?
J(l +
3160 3160 = 3074 m3 /s
Qcorr = 1.67 x 10- 5 ) vi(1057)
2.054 x 1.43 x 10-4
and taking the mean gauge height, the corrected co-ordinates are 50.46 m and
3074 m3 Is.
An alternative procedure, due to Boyer [22] is available, where values for
neither U nor S need be obtained. If a sufficient number of observations are
available, including measurements taken during rising and falling stages and in
steady states, then a rating curve can be drawn as a median line through the
uncorrected values. The steady-state discharge Q can now be estimated from the
134 Engineering Hydrology
median curve. Since Qa and dh/dt are measured quantities and therefore known,
equation 6.9 yields the term l/US for each measurement of discharge.
The term l/US is now plotted against stage and an 'average' curve fitted to
the plotted points, as shown in figure 6.17. From the l/US against stage relation-
ship, new values of l/US are obtained and inserted into equation 6.9 to yield the
steady state Q. The new values of Q are then plotted against stage as the corrected,
steady-state curve.
Another method, which uses the observed slope but avoids evaluating U, is
due to Mitchell [23] .
~~
estimated 'norma]' .~.
curve h~ corrected 'normal'
t .r-
.p' curve
~ ..
..f
~ ..
00 -+ 11US
o
Such a curve plots as a straight line on logarithmic paper and so can be easily
extended. At best it is a questionable procedure since there is little theoretical
justification for exponential laws operating and at high flows there may be quite
abrupt changes in cross-section with rising stage.
(ii) Stevens' method [24].
This method is based on Chezy's formula
Q =ACv'(RS)
where A = cross-sectional area
C = Chezy roughness coefficient
R = hydraulic radius
S = slope of the energy line.
If Cv'S is assumed constant and D, the mean depth substituted for R then
Q =kAv'D
Known values of Av'D and Q are plotted and often come close to a straight
line which can be extended. Field values of Av'D above the measured rating
can then be used from the extended line to plot Q against stage points on the
rating curve.
The objection to this method is simply that Cv'S is not a constant. However,
as it takes account of the varying stream geometry it is a more rational proce-
dure than (i).
Typical values of Chezy's C are given in appendix B.
(iii) Slope-area method. This method depends on hydraulic principles and
presupposes that it is practical to drive in pegs or make other temporary elevation
marks at the time of the flow measurement upstream and downstream of the dis-
charge measuring site. These marks can subsequently be used to establish the
water slope. Cross-sectional measurements will yield the area and hydraulic
radius of the section. Then from Manning
2 I
AKi82
Q=--
n
136 Engineering Hydrology
This method is sometimes criticised because of its dependence on the value of
n. Since n for natural streams is about 0.035 an error in n of 0.001 gives an error
in discharge of 3 per cent. This objection may be partially met by plotting n
against stage for all measured discharges, so that the choice of n for high stages
is not arbitrary but is taken from such a plot. If the high flood slope can be
measured then this method is probably the best one.
Typical values of Manning's n are given in appendix B.
It should be emphasised that all methods of rating curve extension are suspect
to some degree and should be resorted to only if flood measurements cannot be
obtained. The latter two methods above are both subject to errors that arise
from alteration of cross-section as a result of flood scour and subsequent low
water deposition; so cross-sections and mean depth measurements should be
taken as near to the time of discharge measurement as possible.
6.6.1 Hydrographs. With an adjusted and well-measured rating curve, the daily
gauge readings can be converted directly to runoff volumes. A typical set of such
daily runoff figures is presented graphically in figure 6.18. Such a presentation is
called a hydrograph. Although figure 6.18 shows a hydrograph with a time base
of many months, hydrographs for smaller catchments can have time bases of
days or even hours.
While floods and droughts are important from many points of view, they
tend, as extremes, to be of comparatively short duration. For many water-
resource investigations it is equally important to know the total volumes of
5400
5200
4800
4400
4000
$- 3600
.
E 3200 J
2'
o
2800
1\ I V\1
1\
A
~ 2400
/\ )U
,
i5 2000 t--
A rt
1600
1200 t--
~\
IV IV-,'
.., -- f3roundwater "~.\
... /
--
800 contribution ..... h
J ..-
400 ~ I ~
o
Feb. March April May June July August September
6.6.2 Mass curves. If the volumes denoted by the product of ordinate and
time interval of a hydrograph are plotted against time by adding each new
volume to the previous total, a cumulative mass curve of runoff is obtained.
Such a curve is shown in figure 6.19.
24
22
V-- -
V
20
18 ~
_ 16
~ 14
/
V
~ 12
::J
~ 10
~.:: B
::J
Cl: 6
l/
V'
4
2 I-
o
Feb. March April May June July August. September
Figure 6.19 Cumulative mass curve of runoff for the river Euphrates at
Hit, February to September 1957
Mass curves are extremely useful in reservoir design studies since they provide
a ready means of determining storage capacity necessary for particular average
rates of runoff and draw off. Suppose for example that the mass curve OA of
figure 6.20 represents the runoff from a catchment that is to be used for base
load hydroelectric development. If the required constant drawoff is plotted on
the same diagram, as line OB, then the required storage capacity to ensure this
rate can be found by drawing the line CD parallel to OB from a point C at the
beginning of the driest period recorded. The storage capacity necessary is
denoted by the maximum ordinate cd. Normally much longer periods, as long
as the record allows, are used for reservoir design and in many instances the
drawoff is not constant nor continuous. In such cases, different techniques based
on the same principles are used [25].
138 Engineering Hydrology
D A
Time (months,yeors,etc.) ~
6.6.3 Runoff accumulation-time curves. Another way of using flow data in the
design of reservoirs is to calculate reservoir filling times from daily flow readings,
Since river flow varies seasonally, the time taken to accumulate particular
volumes of runoff (equivalent to reservoir filling times) will vary with the time
of year that filling commences, as well as with the stored volume required.
Engineers also need to know the probabilities of certain events taking place
and the risks of not meeting particular objectives. These various criteria are met
by runoff accumulation-time frequency diagrams. A typical example is shown in
figure 6.2l.
The derivation of such curves is straightforward but is too long to include
here. A full description, with worked examples of the techniques required, is
given in reference [26].
6.6.4 Flow duration curves. A flow duration curve (FDC) for a particular point
on a river shows the proportion of time during which the discharge there equals
or exceeds certain values. Such a curve is shown in figure 6.22. Flow duration
curves for long periods of runoff are useful for deciding what proportion of flow
should be used for particular purposes, since the area under a curve represents
volume. Storage upstream of the gauging point in the forms of lakes or reservoirs
will modify the FDC of a river that has been previously without such storage.
Surface Runoff 139
Probability that 50% of annual runoff volume
will be accumulated within lime shown
A <-...,.....--,90%
M 70
ru
F 50
ru -30
a. J
E .. 10.
o
.~ 0
c
.~ N
~ 0
~ S
g A
""
.c
i'
"
>-
n
F M A M J A SON D
Date on which runolf accumulation commences
For many rivers the ratio of peak to minimum discharges may be two or more
orders of magnitude and FDCs for points on them are often more conveniently
drawn with the ordinate (Q) to a logarithmic scale and a normal probability scale
used for the frequency axis. On such a graph, if the logarithms of the discharges
are normally distributed, then the FDC plots as a straight line. This is often
approximately the case.
Flow duration curves for different points on a river can be compared when
presented in this more compact way by standardising them. The discharges are
divided firstly by the area of contributing catchment and secondly by weighted
annual rainfall over the catchment. The resulting discharges, in m 3 /s or 1/s, per
unit area per unit annual rainfall, can then be compared directly. Figure 6.23
shows six FDCs for tributaries of the River Severn in England compared in this
way. A composite regional standardised FDC for the upper reaches of the Severn
and its tributaries was subsequently derived from figure 6.23 [27].
Such standardised curves provide a convenient means of synthesising an FDC
for an ungauged point on a river, and for points on other rivers of similar geo-
morphological nature in a similar climate. Obviously, care should be taken that
these two conditions apply.
140 Engineering Hydrology
.
~
for 17 per cent of the time
0
.t:
U
2
'"
i5
Curve tends to be
,, modified in this way
0 20 40 60 80 100
Time (per cent)
Figure 6.22 Flow duration curve for the river Euphrates at Hit,
February to September, 1955 (derived from the data of figure 6.18)
200
t:-.....
<; 100
-E
<11
>-
"- ~~~
N-
E= 50 ~~
"" ~ ~
-=2
c.~
30 ~
<110
~ ~
0
_0 20
'c E
" c
1;;0
a:-=
c 10
~~
"
1;; ~
~
0.
5
~
~ \~
0 3 \'" ~
\
'\ ~ I::---
K
2
101 2 5 10 20 40 60 80 "" ~
90
I'--.
95
10..
9899
Percentage of time equalled or exceeded
The value of Qm (also referred to asADF or average daily flow) may be derived
for a gauged catchment from the daily flow records by summation and averaging.
2.0
, ...~
~
~ ~
1~
~
"" ~~ ~ OB m
=
S
~
0 ~ ~
0.7 ~.
0.1
~ ~ ~ 0.6 ~
r----- I'---~ =
i"'- X
0.5 <
~
"- ~ ~
0
O~ 0
"- ....... ~
'" 0.4 <
Qro
'" I" ~
!
"-......... 0.3
om
0.01
"" '" ~.2
0.1 2 5 10 20 40 60 80 90 95 98 99 99.9
% probability
Figure 6.24 Typical flow duration curves for values of BFI (by Dr A.
Custard, after Low Flow Studies [26]
Surface Runoff 143
If the catchment is ungauged then Qm must be derived by the use of rainfall
records and estimates of evaporation. In the long term, catchment runoff may be
assumed to be the difference between annual average rainfall (SAAR) and actual
evaporation (AE), since infiltration and percolation will appear eventually as
baseflow.
In the UK, SAAR is estimated from the 1 :625,000 annual average rainfall
map for 1941-70 and potential evaporation (PE) from the 1:2 million map of
annual average PE based on the Penman equation, both available from the
Meteorological Office. Elsewhere PE may be calculated from weather station
data as described in chapter 3.
Actual evaporation (AE) then has to be estimated from potential evaporation
by using an adjustment factor I/J (table 6.1). This factor is based on actual data
of rainfall and runoff for stations in the UK with more than 10 years of data,
and refers to the UK. In other parts of the world it may need modification, using
similar methods.
AE=I/Jx PE
When SAAR exceeds 1100 mm, AE is assumed to equal PE, which is related
to a freely transpiring grass surface. Periods when AE is limited by soil moisture
deficit are assumed to be compensated for by AE exceeding PE at other times by
evapotranspiration from other types of vegetation exceeding that of grass.
If local data are available from a similar neighbouring gauged catchment the
runoff rate may be used directly, but if SAAR for the gauged catchment is
different from the one under study, its loss rate in mm should be deducted from
the SAAR for the latter.
References
Further reading
HERSCHY, R. W. New methods of river gauging. Facets of Hydrology (ed. J. C.
Rodda), John Wiley, New York, 1976, Chapter 5
Surface Runoff 145
HORTON, R. E. Erosional development of streams and their drainage basins.
Bull Geol. Soc. Am., 56 (March 1954) 275
Logarithmic plotting of stage-discharge observations. Tech. Note 3. Water
Resources Board, Reading, 1966.
NASH, J. E. and SHAW, B. L. Flood frequency as a function of catchment
characteristics. Inst. Civ. Eng. Symposium on River Flood Hydrology, 1966,
pp. 115-6
RODDA, J. C. The significance of characteristics of basin rainfall and morpho-
metry in a study of floods in the United Kingdom. Int. Assoc. Sci. Hydrol.
Symposium on Floods and their Computation, Leningrad. International
Association of Scientific Hydrology, 85 (1967) 834
SMOOT, G. F. and NOYAK, C. E. Measurement of discharge by the moving
boat method. U.S. Geological Survey Tech. Water Resources Inv., Bk. 3
(1969)All
STRAHLER, A. N. Statistical analysis of geomorphic research. J. Geol., 62, No. 1
(1954)
Problems
6.1 A river gauging gives Q = 4010 m 3 /s. The gauging took 3 h during which the
gauge fell 0.15 m. The slope of the river surface at the gauging site at the time
was 80 mm in 500 m, and the cross-section approximated a shallow rectangle
200 m wide by 11 m deep. What adjusted value of discharge would you use?
What value of n in Manning's formula results?
lOA 50
12.2 65
13.9 77
14.3 80
22.3 150
27.3 180 -0.32
28.1 228 +0.80
30.8 256 +0.525
32.6 225 -0.36
35.2 251 -0.355
38.9 338 + 0.345
40.3 316 -0.22
40.8 352 +0.18
41.5 333 -0.235
42.2 362
Using Boyer's method adjust the figure for slope variation to produce a steady
flow discharge rating curve for the river.
146 Engineering Hydrology
6.3 Explain how observations of river discharge at particular gauge heights can
be corrected so that they fall on a smooth curve, and explain why this is desirable.
A river discharge was measured at Q = 2640 m 3 /s. During the 100 minutes of
the measurement the gauge height rose from 50.40 to 50.52 m. Level readings
upstream and downstream differed by 100 mm in 700 m. The flood wave
celerity was 2.2 m/so Give the corrected rating curve co-ordinates.
0 0 28 0.7 56 0.6
2 2.0 30 0.8 58 1.2
4 3.2 32 0.8 60 1.4
6 2.3 34 0.7 62 1.8
8 2.1 36 0.7 64 2.0
10 1.8 38 0.5 66 2.3
12 2.2 40 0.4 68 3.2
14 0.9 42 0.7 70 3.4
16 0.5 44 0.8 72 3.5
18 0.3 46 0.4 74 3.7
20 0.7 48 0.3 76 2.8
22 0.7 50 0.2 78 2.4
24 0.6 52 0.2 80 2.0
26 1.2 54 0.4
6.5 The average domestic per capita demand for water in an expanding com-
munity is 0.20 m 3 /day. Industrial demand is 30 per cent of total domestic
requirements. The town has 100000 inhabitants now and is expected to double
its population in future.
Water is supplied from a river system with existing storage capacity of 10 7 m 3
and whose mean daily discharges for each month of the year are as follows
(thousands of m 3 )
Surface Runoff 147
Jan. Feb. March April May June July Aug. Sept. Oct. Nov. Dec.
290 250 388 150 64.5 50 64.5 117 283 388 317 385
April May June July Aug. Sept. Oct. Nov. Dec. Jan. Feb. March
220 250 370 670 865 1630 670 530 270 300 280 280
6.7 list eight characteristics of drainage basins affecting their discharge hydro-
graphs and comment on each.
6.8 A river gauging gives Q = 2060 m 3 /s. The gauging took 2 h during which the
gauge fell 0.12 m. The river surface slope at the time was 5 cm in 400 m and the
cross-section of the river at the site was 300 m wide by 4 m deep.
What adjusted discharge would you use? What is Manning's n for this river at
this time, and what does the value indicate about the river's condition?
6.9 The two-year monthly discharges of a river into a reservoir are listed below
(a) If the reservoir is full at the end of February, year 1, what is the permis-
sible draw-off in lis so that the reservoir may be full at the end of June,
year 2?
(b) If the reservoir is full at the beginning of January, year 1, and the draw
off is 60 lis for the first year and 80 lis for the second year, what will be
the state of storage in the reservoir at the end of December, year 2?
(c) What storage capacity will be required in the reservoir in the second case
to ensure the required discharges may be supplied?
6.10 An unregulated river has monthly mean flows (in m 3 /s) as follows:
Jan. Feb. March April May June July Aug. Sept. Oct. Nov. Dec.
5.4 8.3 9.1 8.8 6.3 6.9 10.2 13.7 19.4 16.7 11.0 21.9
Allowing compensation water of 4.0 m 3 /s and reservoir losses of 0.5 m 3 /s, what
storage capacity of reservoir is required to ensure that, on average, no water is
spilled? What would the average net yield of the reservoir then be? Assume
30-day months.
6.11 (a) Draw a flow duration curve from the tabulated data below -the mean
monthly discharge (flows in m 3 /s) of a source of water.
If a hydro-power plant is to be developed at the site to which the data refers and
where the head available is 15 m, what would be a reasonable first estimate of
plant capacity and annual energy production on the basis of the data? Justify
your choice.
(b) What is the mean monthly flow with a return period of once in 10 years?
Q (m 3 Is) 4.5 3.5 2.9 2.5 2.2 1.9 1.6 1.4 1.1 0.5
7 Hydrograph Analysis
.-
~
.;;:; net or effective rain
c
centre of gravity
~~
_E
:25 ;;;::; loss curve
c .If'-"
~-4---+--~---------------------------
"-
.groundwater "
deplelion curve "-
exlrapolated .....--...:
-----
Groundwater contri bution
to stream flow
Time-
An effluent stream on the other hand is fed by the ground water and acts as
a drain for bordering aquifers. The great majority of streams in Britain and
Europe are in this category.
Intermittent streams are those that act as both influent and effluent streams
according to season, tending to dry up in the dry season.
Perennial streams are greatly in the majority, with a low dry-season flow fed
by baseflow, and are mainly effluent streams, though. many perennial rivers
crossing different geological formations of varying permeability and subject
to different climates are both influent and effluent at different parts of their
courses. A good example of this is the river Euphrates in Iraq. Figure 6.18 shows a
part-annual hydrograph of the river Euphrates and the slow seasonal variation of
the baseflow can be observed. This baseflow is derived principally from the
headwaters of the catchment in northern Iraq, Turkey and Syria. At Hit, where
the hydrograph was observed, the river for much of the year is influent.
Bank storage describes the portion of runoff in a rising flood that is absorbed
by the permeable boundaries of a water course above the normal phreatic
surface, it is illustrated in figures 7.3 and 7.4. In the latter figure the direction of
the arrows showing influx of groundwater to the stream will be reversed during
the flood period while the surface level of the stream is above the phreatic
Bank s lorage
v
1
Flood level
normal oulflow
of groundw-ot-er---
surface. As a result the hydrograph of a particular flood might well have a base-
flow contribution as indicated in figure 7.5. Such a separation as is shown there
is virtually impossible to make quantitatively but it is qualitatively correct.
In many natural rivers, depending naturally on bank permeability and the
slope of the phreatic surface, the variation in base flow will be much less than
indicated in figure 7.S and will cause only a slight dip from the extrapolation of
the depletion curve, followed by a gradual rise to a higher-than-initial value as
indicated in figure 7.6.
'"'"o
.<:
<.>
o'" increase in baseflow
after flood
_--
.......... //..-,
"",,
.....
If baseflow
/
\
"'':'</..'"',"- negative baseflow
lime- or bank storage
.,
~
o
.cc
u prob.oble separation
'"
(5 extropolated line
de.pletion curve
-.~.~~.~~.~.~.~.~ ........ .
baseflow
Time~
Time-
(a) ::::t:
~
a
cc
iil
'0
Recession curves ore plaited so that their :::r
lower ends (where the flow is only boseflow) l>
::I
I
___ ~ to":"'.' to. ~...... liM Master depletion curve ~
iij"
replolled linearly
Cl <:)
go .,
...J 2:1~ ~
o
ReSUllill9l~~t!
plot of the ~ D ~II' .:
u
master depletion curve I i
is'"
---
Time units Time units
(b) (c)
Figure 7.7 Derivation of a master depletion curve. (a) Normal hydrograph with recession curves selected. -0&
U'I
(b) Log plot of recession curves. (c) Linear plot of master depletion curve U'I
156 Engineering Hydrology
Time-
I ID I I I I
120 I--+--+-I--+--+-~""'De(ak - 1;)+1'3~ hange Df slope-
.. 1;)......, /
110 1--+-+--+--+-1-t7"I~......
I
-+-
"' ~
"2 ~'....
100 I--+--+--I-+--+-++-rl---+>t/i. 1'1 "',
90 r-+--r-+-~~-+-H~+-'\
I I \. I
I'O~~~~~~~
I I ....; .... '1
6-11 intervals
80
I i\
~ r-+--+----1I-+-+I+--r...,'I+--+-----1I-\-\-+-~ThiS vertical
E 70
I I N' \ I the 6 h time interval
bisects
eo 60 I---l--+-I--hl-+-L'*==*==I=~~ between Oand 0. 6 for
0
.,
.J:.
50
I
surface runoff / \ I the point at change of slope
\1:
u
'0 I
40
)
30 . __ f.--
20
1/ -I-
baseflow
10
2 3 4 5 6 7
Time (days)
critical ratio can be determined and the first point beyond the region of inter-
section on the groundwater side gives a conservative position for N. The subse-
quent hydrograph analysis is not very sensitive to the precise position of Nand
either 0300 h or 0600 h on day 6 would be satisfactory. A straight line is now
drawn to N from the point where the hydrograph started to rise, as before, The
total amount of runoff can now be obtained by measuring the area under the
curve and above the straight line.
The position of N is important also in synthesising hydrographs, as will be
seen in section 7.11, since it partly defines the baselength of the hydrograph.
The baselength is made up of the part before the peak, which depends on the
length of the period of rain and te the period of concentration, and the recession
limb after the peak, which depends primarily on the character of the catchment.
From observations on many natural catchments the position of N can be estab-
lished empirically from table 7.2 in days after the peak of the flood.
250 2
1250 3
5000 4
12500 5
25000 6
"I
Figure 7.10 Proportional principle of the unitgraph
?:-
" ..,
c
.s net rain only
c
~ ~~~~~--~~-----------------------
Time
Having established these principles the concept of unit rain is now introduced.
A unit of rain can be any specified amount, measured as depth on the catch-
ment, usually 1 cm or 1 in. but not exclusively so. The unit rain then must all
appear as runoff in the unit hydrograph. The area under the curve of the hydro-
graph has the dimensions of instantaneous discharge multiplied by time, or
L3
- - x T:: L 3 :: volume of runoff
T
160 Engineering Hydrology
so that although unit rain is spoken of as 1 cm over the whole of the catchment
area the resulting runoff is given in cubic metres, and the quantities involved are
identical. If the unitgraph for a particular catchment, and a particular duration
of rain is known, then from principle 2, the runoff from any other rain of the
same duration can be predicted.
This is a first step towards the complete correlation sought, but if the rainfall
should be of different duration from that of the unitgraph then the unitgraph
must be altered before it can be used.
10
.....
Time (h)
Time(h)-
'.xcE 1/ t , m'77?777.'r77;'V77777."7n'77771'7:"7n777.'r77;'?V7777-"7n'777777-'7777:0
E
278A (km')
o Q.= t,(h)-
o Time-
If the time base of the unitgraph is T h, then steady-state runoff must occur
at T h and so only Tit 1 unitgraphs are necessary to develop constant outflow
and so produce an S-curve. The equilibrium flow, Qe, can easily be obtained
since 1 cm on the catchment is being supplied and removed every t 1 h:
Engineering Hydrology
or Qe = 645A
t1
It will be apparent that the method can be used for altering the unit period
either way, longer or shorter, and that if changing from shorter to longer dura-
tion then t 2 need not be a direct multiple of t l' Although the method has been
described graphically, in practice its application is usually made in tabular form
and example 7.1 illustrates it.
Example 7.1. Given the 4-h unit hydrograph listed in column (2) of table 7.3
derive the 3-h unit hydrograph. The catchment area is 300 km 2
0 0 0 0 0
1 6 6 6 8
2 36 36 36 48
3 66 66 0 66 88
4 91 0 91 6 85 113
5 106 6 112 36 76 101
6 93 36 129 66 63 84
7 79 66 145 91 54 72
8 68 91 159 112 47 63
9 58 112 170 129 41 55
10 49 129 178 145 33 44
11 41 145 186 159 27 36
12 34 159 193 170 23 31
13 27 170 197 178 19 25
14 23 178 201 186 15 20
15 17 186 203 193 10 13.5 a
16 13 193 206 197 9 12a
17 9 197 206 201 5 6.5 a
18 6 201 207 203 4 5.5 a
19 3 203 206 206 0 Oa
20 1.5 206 207 206 1 1.5 a
21 0 206 206 207 - 1
aSlight adjustment is required to the tail of the 3-h unitgraph. This is most easily done
by eye (see figure 7.14).
Hydrograph Analysis 163
120
110 I-
1\
100
/ it,
fi I "
\~:-h unilgraph derived from S-curve
I
90 r-
4-h unilgraph given
80
I/
\~\
-;;; 70 r-
;;,.
'.
'~\
.eo
. 60
<> 50
.c;
u
on
0 40
./
I
'\\
'\~
..
t \\
30 \~,
ff ,t~
20
".~t
10 -#
:j
.~
;~ ...
1 1 ::1"".....L......
00 2 4 6 8 10 12 14 16 18 20 22
Time (h)
It will be noted that Qe = 208 m3 /s, as calculated, agrees very well with the
tabulated S-curve terminal value 207. This is an indication that the 4-h period of
the unitgraph is correctly assessed. Very often with an uneven rainfall distribu-
tion, an attempt has to be made to reduce the net rain to a uniform-intensity
rain of a particular duration. The S-curve can in this way serve as a check on the
chosen value. If the S-curve terminal value had fluctuated wildly and not steadied
to a minor variation it would have indicated an incorrect rainfall-time for the
unitgraph.
Note also that it was not necessary in table 7.3 to set out Tit 1 columns of the
4-h unitgraphs, and add them laterally. The S-curve additions are the S-curve
ordinates shifted in time by 4 h. Since the first 4 h of unitgraph and S-curve are
the same, the S-curve additions and S-curve columns are filled in, in alternate
steps. The effect is the same as setting out rows of unitgraph ordinates succes-
sively staggered by 4 h, since the S-curve additions represent the sum of all
previous unitgraph ordinates.
40
32
C
Cl> 30
u
~
3-
cv
~ 20
0 20
<:
u
17
'"
U
0
; 10
4
05
0
Unit periods
(b)
line are the same and so, in deriving unitgraphs from distribution percentages, a
smooth line must be drawn through the steps to give equal areas.
The distribution graph is therefore less precise than the hydrograph but is
much better suited for iterative processes of derivation, as will be seen in section
7.8.
Hydrograph Analysis 165
7.8 Derivation of the unit hydrograph
The unitgraph for a particular catchment can be derived from the natural hydro-
graph resulting from any storm that covers the catchment and is of reasonably
uniform intensity. If the catchment is very large that is, greater than (say) 5000
km2), it may never be covered by a uniform-intensity storm, since these are
limited in size by meteorolOgical conditions. In such a case the catchment should
be divided up into tributary catchments and the unitgraphs for each of these
determined separately.
The first step is to separate the baseflow from surface runoff (section 7.3)
and plot the runoff and the rain graph on the same time base. The quantity
of net storm rain must then be estimated and its intensity and duration estab-
lished. A check is now made on the quantity of net rain on the catchment and
the amount of runoff under the hydrograph. These should be the same and one
or the other may require adjustment.
The unitgraph can now be obtained by dividing the runoff hydrograph
ordinates by the net rain in cm. The adjusted ordinates represent the unitgraph
for the particular duration established.
It is always advisable to determine several unitgraphs, using separate and
distinct isolated uniform-intensity storms, if available. Natural events like rain-
storms and runoff are affected by a multiplicity of factors and no two are
precisely the same. Frequently the best natural data will be for different rain
durations and the resulting unitgraphs will require to be altered to the same
duration (section 7.6). Once a number of such hydrographs has been obtained
for the same duration, an 'average' or typical unitgraph can be determined as
shown in figure 7.16. The ordinates are not averaged since this would produce
an untypical peak. The peak values of the sel'arate unitgraphs are averaged as are
the values of the time from the beginning of runoff to the peak. These values are
assigned to the average unitgraph which is then sketched in to a median form on
both rising and falling limbs, so that the total area under the curve is equal to
1 cm runoff.
(
~t '2ft3 PI+P2+P3 )
3 3
overage unitgraph
sketched to fit and
to canto in unit runoff
, Time
base flow being assumed previously separated. The process may be illustrated
by reference to figure 7.17.
The first rain period, of duration and intensity t and i I respectively, gives rise
to runoff illustrated by the hypothetical hydrograph bounded by the lower
dashed line. Each ordinate of this hydrograph is ti I times the unitgraph ordinate
U I , U 2 . . . Un. Similarly, the second and third rains of intensity i2 and i3
respectively produce additional runoff whose ordinates in each case are ti 2 and
ti 3 , multiples of the t h unitgraph shifted appropriately in time. If the complex
hydrograph is now defined by ordinates at suitable equal intervals (conveniently
but not essentially fixed as a whole multiple of t h) then the first ordinate of the
unitgraph, UI, is obtained from Ql = tilU1 where Ql is the observed runoff,
hence U I can be found. For the second ordinate, Q2 =ti 1U2 + ti2U1 in which
equation U2 is the only unknown.
The third ordinate is similarly obtained from Q3 = ti 1 U 3 + ti 2U2 + ti 3U1
where U3 is now the only unknown. Proceeding in this way, the t h unitgraph
ordinates can be successively determined.
In the above illustration all the rain periods, although of different intensities,
were assumed to be of the same duration t h. This is a condition for the use of
the method since otherwise other variables U~. U; etc. (the ordinates of the t' h
unitgraph) would be introduced.
Although the method appears simple, since each U ordinate depends on the
preceding ones, and on the assumptions about intensity and duration of rainfall
Hydrograph Analysis 167
..e-
.o
s::.
u
Ci
and deduction of an assumed base flow , errors accumulate and several trials and
restarts may be necessary to find a reasonable unitgraph.
The second method is Que to Collins [5] and is amongst the simplest of
various iterative methods proposed. To illustrate its use a unit hydrograph will
be derived from the data of rainfall and natural discharge for the catchment of
the River Rother at Woodhouse Mill in Yorkshire.
(1) Assemble all rainfall data available for the catchment under considera-
tion and the storm period, including daily observations, recording rain-
gauge records and the synoptic weather maps of the region, if available.
(2) Derive a mean mass curve of rainfall for the catchment, for the period of
rain producing the hydrograph under study. Make provisional separation
of the rain into separate uniform periods.
The catchment of the River Rother at Woodhouse Mill, Yorkshire is
shown in figure 7.18. The continuous rainfall record (mass curve) at
Catchment area =
136 miles 2
14 1 1 J
~
c:'
12
Sulton rain recorder
trace
---/1 ~
/ /
V/
'0
..
.t:.
Ci.
"0
0'8
/1
c
1 rainfall rates assumed
...
.~
W in the analysis
0'6
V
~
1!
"E 0'4
"
u ,,/
0'2
,./
~
~
~
t200 1400 1600 '800 2000 2200 2400 0200 0400 0600 0800
- - - - - 14th May - - - - - -....~tt__--- 15th May--
0'40
0'2
0'05 Vnet rain =065in.
0'1 ~ I index 00()5 in./h
0
t
2200
2000
Lf~
800
I ~/ Observed hydragraph
.f.
EiOO
~
VDerived distribution graph
- I 400
~
SI 200
'"~ 1000
o
s;;
u
800 r
VOlume=206'5l(lcf~- ~
-}
Runoff
i:5'"
r---..
"O65in. on '36 m',z
60 0
400
20C kiC -- - -- - Baseflow
L----f - -
I .1 I I I I I I
o 4 8 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80
Time (h)
(6) Tabulate the relevant data in columns (1)-(9) of a table similar to table
7.4 and provide columns for the number of unit periods in the unit
hydrograph under 'Distribution coefficients'.
The number of columns in table 7.4 under 'Distribution coefficients' is
16. This is the number of unit periods from the beginning of the last rain
period to the end of surface runoff
(7) Assume distribution coefficients of the unitgraph (representing percentage
of total runoff per unit period) and arrange in appropriate columns.
(8) Determine the discharge, which, flowing constantly for one unit period
would just equal 1 in. of net rain on the catchment.
This figure is found in this case to be
(9) The first net rain is multiplied by this discharge and the product is
distributed in percentages across the distribution coefficient columns by
being multiplied in turn by each percentage coefficient. The various
numbers are entered in the columns diagonally.
In this case 0.05 x 21960 = 1098, so in the first column 0.05 x 1098 =
55 and so on. Note that 55 was entered opposite the corresponding rain
and not on the top row, which is ignored.
(l0) The procedure of (9), above, is repeated for all net rains except the largest
for which a dash is entered throughout.
It is purely fortuitous that the largest rain is at the end in this case.
TABLE 7.4 Derivation of unit hydrograph from a multi-period storm
(I) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11 ) (12) (13)
Day Time Period Rain Losses Net Aver- Basf!- Net Distribution CoeJflcients (per <:ent) Resi-
" /J"r
no. (in.) rain (in.) age flow Q 5.0 10.0 14.0 15.0 /3.5 //.0 8.5 6.2 4.0 3.3 2.6 2.1 1.7 1.4 1.0 0.7 duals CCllt
2nd trial coefficients 6.4 10.9 12.6 13.5 13.4 12.5 9.4 6.7 4.1 3. 2.4 1.8 1.3 1.0 0.6 0.4 100.4
3nl trial coefficients 7.2 11.5 12.3 12.8 12.9 12.3 9.5 6.7 4.2 3. 2.5 1.8 1.3 1.0 0.5 0.3
Accepted coefficients 7.0 11.6 12.5 12.9 12.7 12.0 9.6 6.8 4.3 3. 2.5 1.8 1.3 1.0 0.5 0.3
...........
172 Engineering Hydrology
(1 I) The various discharges are now summed horizontally and entered in the
~ column.
(12) The ~ column discharge totals are now subtracted from column (9) and
the remainders entered in the 'Residuals' column. These residuals are
then converted into percentages of the unit distribution graph by dividing
by the discharge of (8) (see page 170) multiplied by the largest rain, and
subsequently multiplying by 100. The sum of the percentage column
should be 100. The percentages that cannot have been influenced by the
largest rainfall are bracketed and redistributed over the other coefficients
so that the total, 100 per cent, remains constant. These percentages
represent the distribution that would be necessary for the largest rain to
make good the net Q of column (9). If they are the same as the assumed
distribution coefficients then the unit distribution graph has been deter-
mined.
The residual 1090 of period 6 for example is converted thus
1090
x 100 == 12.5 per cent
0.40 x 21960
C3 -- -
QC l +PC2
----
Q+P
28 --.
~f-
2 600
2400 ~
~~,
Distribution graph from the
2200 vaccePted coefficients of table 7A
2 I
I
~~,
..
I
\
~ I"""
..
I
~
u I
000 t
i5 I
I
800
I
I
600 I
400
'1.::::......
200
;: ~ I'-c::..-
o 4 8 12 16 20 24 28 32 36 40 44 48 52 56 6064GB
Time (h)
The method is particularly useful when the largest rainfall is very great com-
pared with the others since rapid convergence of the coefficients then takes
place. This was not the case in the example illustrated. Too much relianc~ should
not be placed on unitgraphs derived in this way until they have been used in
practice and/or derived from a series of different storms, since the loss rates
chosen have a critical influence on the resulting unitgraphs.
1-
U(O, t) = d~~t)
that is, the ordinate of the IUH at any time t is given by dSt/dt at time t.
The IUH is a unique demonstration of a particular catchment's response to
rain, independent of duration, just as the unitgraph is its response to rain of a
particular duration. Since it is not time-dependent, the IUH is thus a graphical
expression of the integration of all the catchment parameters of length, shape,
slope condition etc. that control such a response.
The conversion of an IUH to a unitgraph of finite period is simple. The
ordinate of an n-h unitgraph at time t is the average ordinate of the IUH for
n h before t. From figure 7.23, it can be seen that the IUH is divided into n-h
time intervals, and the averages of the ordinates at the beginning and end of each
interval are plotted at the end of the intervals to provide the n-h unitgraph.
The IUH can be used to derive unitgraphs by flood routing as explained in
section 8.6.
Instantaneous unitgroph
which there are no runoff records at all and for which unitgraphs may be
required. In these circumstances hydrographs can be synthesised on the basis of
past experience in other areas and applied as first approximations to the unre-
corded catchment. Such devices are called synthetic unitgraphs.
The original approach is due to Snyder [6] who selected the three parameters
of hydrograph base width, peak discharge and basin lag as being sufficient to
define the unit hydrograph. These are shown in figure 7.24.
Snyder considered the catchment characteristics likely to affect unit hydro-
graph shape as being catchment area, shape of basin, topography, channel slopes,
stream density and channel storage. He eliminated all these parameters except
the first two by including them in a coefficient Ct. He dealt with the size and
shape of catchment by measuring the length of the main stream channel and he
proposed that
tp = Ct (L ca L)o.3
where tp =basin lag in h (that is, the time between mass centre of unit rain of
tr h duration and runoff peak flow.
Lea = distance from gauging station to centroid of catchment area, meas-
ured along the main stream channel to the nearest point, in miles.
L = distance from station to catchment boundary measured along the
main stream channel, in miles.
Ct = a coefficient depending on units and drainage basin characteristics
and varying between 1.8 and 2.2 for the Appalachian Highlands
catchments studied.
176 Engineering Hydrology
-- ..
"0
c: c:
-- G>
0_
~
.~
Ir .S
t,
net storm rain
.. "\
synthetic unitgraph
t
The equation for peak flow (per square mile of area) of the tr unit graph was
given by
640
qp =Cp.~ (7.3)
tp
where qp is in cubic feet per second, and Cp is a coefficient depending on units
and basin characteristics that varies between 0.56 and 0.69 for the Appalachian
catchments and generally approaches its largest value as Ct approaches its lowest
and vice versa.
Therefore the peak flow for the unit graph is
Q =C 640A
p p'
tp
where A = catchment area in square miles.
The duration of surface runoff, or unit hydrograph baselength, T was given
by Snyder by the empirical expression
T=3 +3 (~)
where T is in days and tp in hours. This expression gives a minimum baselength
of 3 days for even small areas, a period much in excess of delay attributable to
channel storage.
Hydrograph Analysis 177
Snyder comments on this as being due to the 'subsurface storm flow', which
has been defined by Hursh [7] as 'that portion of the storm flow which inftl-
trates into the surface soil but moves away from the area through the upper
soil-horizons at a rate much in excess of normal groundwater seepage'. This is
what is referred to in recent literature as interflow and for most practical pur-
poses it is regarded as surface runoff.
t r , the unit rain period, was assumed to equal tp/5.5 in the study, since it was
necessary to choose a single standard in all the catchment for the derivation of
the formulae. This particular value was selected to make the unit of time equal
to a minimum value below which further decrease would have little or no effect
on basin lag or unitgraph peak discharge. If the actual length of the storm is not
equal to tro but is tR, equation 7.3 becomes
where qpR =peak discharge (per square mile) of the tR unitgraph, which allows
for the generally observed reduction in unitgraph peak flows with longer periods
of rain. Snyder proposed subsequently [8] an expression to allow for some
variation in basin lag with variation in effective rainfall duration
tpR = tp + (tR ~ t r)/4
where t pR = basin lag for a storm of duration tR'
linsley subsequently presented data [9] based on a study of Californian
catchments and suggested modifications of Snyder's formulae and gave values of
the various coefficients, as follows.
Basin lag: a new factor tpo was introduced-the basin lag of an instantaneous
storm~and used to derive tpR, which has the same meaning as before.
C' = 0.6/v'(Sst)
Sst = [ 1 I
n 1
J2
(lIS{ + l/Si + ... + l/Sit)
- - - 10 mm intervals
- - - - Intermediate intervals
30 < Minimum
> Maximum
Scale of kilometres
6---------+--------~------ .------+ 0 100 -------I
6 6~
Scale of stnute miles
,
"
, ,- -15.
.
:..... /
'I '
'\
/
-,,
UTM Grid
Zone 30U 4 6
events that equal or exceed a particular level (see also section 9.3.2); the return
periods of floods and their causative rainfall are not necessarily the same, since
there are other factors involved. However, the FSR provides a relationship for a
recommended storm return period to yield a flood peak of required return
period by the design method described, from which figure 7.28 has been derived.
Hydrograph Analysis 181
10 mm intervals
Intermedliue Interval"
12- M.ximum - - -+_--iJ,
>
Scole 01 k,lometre!
oI 100
!
o ~O
Sule or '5t.atutc mires.
11- - - - + -- - -- + - - - - - + - - - - - - - 1
IO----T-----r-----r-~~~-+_---~
Example 7.2. By means of a synthetic unit hydrograph, estimate the flood peak
with a 200-year return period for the River Brathay at Skelwith Force at NGR
NY34 103 2 (figure 7.29).
182 Engineering Hydrology
S<>I. 01 kilome.,etoo
~- r
o
Stale 01
-T
60
"0
o 30h-,-,.;.;.,.;.;rt;;;tt;;;a,m
o
LL
20~~~~~~~HTIrrm
20 304050 100
Rain storm return period (yeors)
t
Kilometres
oL 1
I
Figure 7.29 The catchment area of the River Brathay at Skelwith Force
(1) Outline the catchment area on the 1 :25 000 map and measure the area
(AREA) and length of the main stream MSL:
AREA = 50.4 km 2
MSL = 12.2 km
(2) Channel slope (S1085) is the average slope in m/km between two points
at 10 per cent and 85 per cent of the main stream length measured from
the outlet. Since
elevation at 85 per cent point = 137 m
and
elevation at 10 per cent point = 58 m
therefore
S1085 = 8.6 m/km
(3) The average annual rainfall (SAAR) is obtained by taking a weighted
average over the catchment:
SAAR = 2500 mm
and knowing
ARF (table 2.8) = 0.95
gives I-day MS catchment rainfall = 86.6 mm
Finally deducting s.m.d. (figure 4.8) (= 4 mm) gives 82.6 mm
RSMD = 82.6 mm
or rounding up
RSMD=83 mm
(8) The design storm duration can now beI calculated from equation 7.5:
D =(I + SAAR/1000)Tp
= (1 + 2500/1000)5
= 17.5
It is convenient to make D an odd integer multiple of T.
Hydrograph Analysis 185
Hence
D = 17 h.
(9) It is now necessary to decide the return period of storm that will produce
the appropriate return period peak flow. For a 200-year peak flow the
recommended return period is obtained from figure 7.28: the storm
return period is 240 years.
and
62 per cent of 2-day MS = 0.62 x 13S = 83.7
Hence take
17 h MS as 84 mm.
(11) It is now necessary to convert 17-h MS to 17-h M240. From table 2.6,
for MS = 84 mm, by interpolation M240 requires a growth factor of 1.71,
therefore
17 -h M240 = 84 x 1.71 = 144 mm
(12) This is a point rainfall value and must now be reduced to a catchment
average using ARF (table 2.8). Interpolating
ARF (for 17 h and SO km 2 ) = 0.9S
Therefore
rainfallP= 137 mm
(15) This net rain must now be applied to the unit hydrograph in accordance
with a 75 per cent winter storm profile. This distribution should take
place over the 17-h storm period and so a stepped distribution graph over
17 one-hour periods is required. Each interval will be approximately
6 per cent of duration. The percentage rain is obtained from figure 2.l7.
Each of the increments, except the first (the peak), represents the total of two
equal 1 h periods of rain. This rain is now arranged symmetrically about the
centre line and applied to the synthetic unit hydrograph.
The stepped distribution graph of rain is shown in figure 7.30 and its applica-
tion to the synthetic unitgraph is made in table 7.5.
15
r:
.....
E 10
.5
5
5 10
T.
40
30 Qp
20
TB
10
0 5 10
Time (h)
Figure 7.30 Synthetic 1-h unitgraph and hyetograph for example 7.2
1 0.08 0 0.4 0.7 1.1 1.4 1.8 1.6 1.3 1.1 0.9 0.7 0.4 0.2 0 0 0 0
2 0.17 0 0.8 1.5 2.3 3.0 3.8 3.3 2.8 2.4 1.9 1.4 0.9 0.5 0 0 0
3 0.21 0 0.9 1.9 2.8 3.7 4.7 4.1 3.5 2.9 2.3 1.8 1.2 0.6 0 0
4 0.25 0 1.1 2.2 3.3 4.4 5.6 4.8 4.2 3.5 2.8 2.1 1.4 0.7 0
5 0.42 0 1.9 3.7 5.6 7.5 9.3 8.2 7.0 5.8 4.7 3.5 2.3 1.2
6 0.50 0 2.2 4.4 6.7 8.9 11.1 9.7 8.3 6.9 5.6 4.2 2.8
7 0.75 0 3.3 6.7 10.0 13.3 16.6 14.6 12.5 10.4 8.3 6.2
8 1.14 0 5.1 10.1 15.2 20.2 25.3 22.1 19.0 15.8 12.6
9 1.34 0 6.0 11.9 17.8 23.8 29.7 26.0 22.3 18.6
10 1.14 0 5.1 10.1 15.2 20.2 25.3 22.1 19.0 m
::::J
II 0.75 0 3.3 6.7 10.0 13.3 16.6 14.6 cc
12 0.50 0 2.2 4.4 6.7 8.9 11.1 :i"
13 0.42 0 1.9 3.7 5.6 7.5 CD
CD
14 0.25 0 1.1 2.2 3.3 ...S
15 0.21 0 0.9 1.9 cc
16 0.17 0 0.7
:::J:
17 0.08 0
8.38 0 0.4 1.5 3.5 6.7 11.7 18.3 27.0 39.6 55.9 74.5 92.3 107.6 116.2 116.6 109.9 99.5
-a.a
Base flow
0"
3.3 cc
-<
119.9
The choice of the s.m.d. on which a design storm may fall is significant. It
is not always reasonable to presuppose a saturated catchment, and a hypothetical
inflitration capacity curve, a <I>-index or coaxial correlation graph (for example
figure 4.7) can be used with this in mind. For catchments in the British Isles,
when no other information is available, the mean value of s.m.d. from figure
4.8 should be used.
190 Engineering Hydrology
7.13.2 The l-h unit hydrograph
In recent years the concept of design rainfall has changed to simulate natural
events more precisely. Rain does not fall naturally at uniform intensity for long
periods. With computers now in common use, the simplications that were
previously desirable are no longer necessary since the tedious calculations involved
need no longer be performed manually.
Rainfall events are now specified as total depths (mm) falling with particular
durations and to particular storm profiles, the profile being compared in terms
of their peakedness and seasonality. Thus, the United Kingdom 75 per cent
winter profile is a storm profile that is typical of storm rainfall from October
to March, or the winter season, in Britain, and is described in section 2.9.5.
For design purposes, generally large events are postulated. Rainfall durations
thus tend to be of many hours length. Accordingly, it is frequently convenient
to use a I-h rain increment applied to a I-h unit graph, and to convolute the net
rain profile with the hourly ordinates of the unitgraph, as was done in example
7.2.
References
1. Natural Environment Research Council. Flood Studies Report, Vol. I,
NERC, 1975, Chapter 6
2. SHERMAN, L. K. Stream flow from rainfall by the unitgraph method. Eng.
News Record, 108, (1932) 501
3. BERNARD, M. An approach to determinate stream flow. Trans. Am. Soc.
Civ. Eng., 100,(1935) 347
4. LINSLEY, R. K., KOHLER, M. A. and PAULHUS, J. L. H. Applied Hydrology,
McGraw-Hill, New York, 1949, pp. 448-9
5. COLLINS, W. T. Runoff distribution graphs from precipitation occurring in
more than one time unit. Civ. Eng., 9, No. 9 (September 1939) 559
6. SNYDER, F. F. Synthetic unitgraphs. Trans. Am. Geophys. Union, 19th
Annual Meeting 1938, Part 2, p. 447
7. HURSH, C. R. Discussion on Report of the committee on absorption and
transpiration. Trans. Am. Geophys. Union, 17th Annual Meeting 1936,
p. 296
8. SNYDER, F. F. Discussion in ref. 6
9. LINSLEY, R. K. Application of the synthetic unitgraph in the western moun-
tain States. Trans. Am. Geophys. Union, 24th Annual Meeting 1943,
Part 2, p. 580
10. TAYLOR, A. B.and SCHWARZ, H. E. Unit hydrograph lag and peak flow
related to basin characteristics. Trans. Am. Geophys. Union, 33, (1952)
235
I!. Institution of Civil Engineers. Floods and Reservoir Safety: an engineering
guide, LC.E., London, 1978
12. PAULHUS, J. L. H. and GILMAN, C. S. Evaluation of probable maximum
precipitation. Trans. Am. Geophys. Union, 34, (1953) 701
13. Generalised estimates of probable maximum precipitation over the U.S.
east of the 105th meridian. Hydrometeorological Report No. 23, U.S.
Weather Bureau, Washington D.C., 1947
Hydrograph Analysis 191
14. Generalised estimates of probable maximum precipitation of the United
States west of the 105th meridian for areas to 400 square miles and
durations to 24 hours. Tech. Paper 38, U.S. Weather Bureau, Washington
D.C., 1960
15. Manual for depth duration area analysis of storm precipitation. U.S. Weather
Bureau Co-operative Studies Tech. Paper No. 1, Washington D.C., 1946
16. HERSHFIELD, D. M. Estimating the probable maximum precipitation. Proc.
Am. Soc. Civ. Eng., 87, (September 1961) 99
17. National Environment Research Council. Flood Studies Report, Vol. 11,
NERC, 1975
18. BLEASDALE, A. The distribution of exceptionally heavy daily falls of rain
in the United Kingdom. J. Inst. Water Eng., 17, (February 1963) 45
19. WIESNER, C. J. Analysis of Australian storms for depth, duration, area data.
Rain Seminar, Commonwealth Bureau of Meteorology, Melbourne, 1960
Further reading
BARNES, B. S. Consistency in unitgraphs. Proc. Am. Soc. Civ. Eng., 85, HY8
(August 1959) 39
BUlL, J. A. Unitgraphs for non uniform rainfall distribution. Proc. Am. Soc.
Civ. Eng., 94, HYl (January 1965) 235
MORGAN, P. E. and JOHNSON, S. M. Analysis of synthetic unitgraph methods.
Proc. Am. Soc. Civ. Eng., 88, HY5 (September 1962) 199
MORRIS, W. V. Conversion of storm rainfall to runoff. Proc. Symposium No. 1,
Spillway Design Floods, N.R.C., Ottawa, 1961, p. 172
Problems
7.1 A catchment area is undergoing a prolonged rainless period. The discharge
of the stream draining it is 100m3 /s after 10 days without rain, and 50 m 3 /s
after 40 days without rain. Derive the equation of the depletion curve and
estimate the discharge after 120 days without rain.
7.3 The recession limb of a hydrograph, listed below, is to be divided into run-
off and baseflow. Carry out this separation
(a) by finding the point of discontinuity on the recession limb,
(b) by finding the depletion curve equation and extrapolating back in time.
Comment on your results.
192 Engineering Hydrology
Time Flow Time Flow
(h) (m 3 /s) (h) (m 3 /s)
15 41.1 33 10.0
18 35.8 36 8.3
21 25.0 39 7.0
24 19.2 42 5.8
27 15.1 45 4.9
30 12.2 48 4.1
7.4 Describe how to derive a master depletion curve for a river. What would
you use it for and why?
0 3 24 25
3 3 27 21
6 10 30 17
9 25 33 13.5
12 39 36 10.5
15 43 39 8
18 37 42 5.5
21 30.5 45 4
48 3.9
7.6 Describe how you would proceed to separate baseflow from the hydro-
graph of a stream's discharge.
7.7 Write down the three major principles of unit hydrograph theory illustrating
their application with sketches.
Given below are three unit hydrographs (all values in ft. 3 /s) derived from
separate storms on a small catchment, all of which are believed to have resulted
from 3-h rains. Derive the average unit hydro graph and confirm its validity if
the drainage area is 5.25 square miles.
Hours Storm 1 Storm 2 Storm 3
o o o o
165 37 25
Hydrograph Analysis 193
2 547 187 87
3 750 537 260
4 585 697 505
5 465 608 660
6 352 457 600
7 262 330 427
8 195 255 322
9 143 195 248
10 97 135 183
11 60 90 135
12 33 52 90
13 15 30 53
14 7 12 24
15 0 0 0
7.8 Listed below are three hydrographs derived from three separate uniform-
intensity storms each lasting 3 h. The gross rainfall for storm A was 14 mm, for
storm B was 24 mm, and for storm C was 19 mm. The cI>-index for the catch-
ment is estimated to be It mm/h. Derive the average unit hydrograph for the
catchment and confirm its validity if the drainage area is 13.60 km 2 All values
are in m 3 /s.
0 0 0 0
1 1.84 0.82 0.42
2 6.10 4.16 1.46
3 8.36 11.98 4.35
4 6.52 15.54 8.45
5 5.18 13.56 11.04
6 3.92 10.20 10.04
7 2.92 7.36 7.14
8 2.17 5.68 5.39
9 1.59 4.34 4.16
10 1.08 3.02 3.06
11 0.67 2.00 2.27
12 0.37 1.16 LSO
13 0.17 0.66 0.89
14 0.08 0.26 0.41
15 0 0 0
194 Engineering Hydrology
7.9 The 3-h unit hydrograph derived from a catchment of 14.5 km 2 is given
below
What peak discharge would be expected from a 4-h rainfall at a uniform intensity
of 15 mm/h, followed immediately by a 3-h storm at a uniform intensity of
10 mm/h? Assume a constant storm loss of 3 mm/h and a base flow starting at
1.2 m3 /s at hour 0 and rising at 0.1 m3 /s per hour until after the peak.
7.10 The hydrograph tabulated below was observed for a river draining a 40
square miles catchment, following a storm lasting 3 h
Separate baseflow from runoff and calculate total runoff volume. What was
the net rainfall in inches per hour? Comment on the severity and likely frequency
of such a storm in the United Kingdom.
7.11 The 4-h unit hydrograph for a river-gauging station draining a catchment
area of 5 S4 km 2 , is given below.
Hydrograph Analysis 195
0 0 11 76
1 11 12 62
2 60 13 51
3 120 14 39
4 170 15 31
5 198 16 23
6 184 17 16
7 153 18 11
8 127 19 6
9 107 20 3
10 91 21 0
Make any checks possible on the validity of the unitgraph. Find the probable
peak discharge in the river, at the station from a storm covering the catchment
and consisting of two consecutive 3-h periods of net rain of intensities 12 and
6 mm/h respectively. Assume baseflow rises linearly during the period of runoff
from 30 to 70 m3 /s.
7.12 The 3-h unit hydrograph for a river-gauging station draining an 835 km 2
catchment is listed below.
0 0 12 85
I 22 13 68
2 120 14 55
3 240 IS 42
4 318 16 30
5 298 17 20
6 250 18 12
7 206 19 7
8 174 20 3
9 144 21 0
10 123
11 102
An intermittent storm lasting 7 h covers the catchment, the gross rainfall being:
17 mm/h for 4 h, followed by 12 mm/h for a further 3 h. The I>-index for the
catchment is 7 mm/h. Assuming baseflow is constant at 40 m 3 /s, provide an
estimate of the maximum discharge and its time of occurrence from the start of
the storm.
196 Engineering Hydrology
7.13 Using the data and catchment of 7.11 find the probable peak discharge
in the river, at the station, from a storm covering the catchment and consisting
of three consecutive 2-h periods of rain producing 7, 14 and 12 mm runoff
respectively. Assume baseflow rises from 10m3 /s to 20 m 3 /s during the total
period of runoff.
7.14 The 4-h unit hydrograph for a 550 km 2 catchment is given below
Hours Q Hours Q
(m 3 /s) (m 3 /s)
0 0 11 76
1 11 12 62
2 71 l3 51
3 124 14 40
4 170 15 31
5 198 16 27
6 172 17 17
7 147 18 11
8 127 19 5
9 107 20 3
10 90 21 0
7.15 Using a FSR (1975) synthetic unit hydro graph with Tp = 8 h for a catch-
ment of 350 km 2 , deduce the surface runoff contribution to the discharge of the
river draining the catchment at the end of the fourth hour of continuous rain of
10 mm/h, if the runoff coefficient is 70 per cent. Assume runoff starts at start of
rain.
Define the parameters on which Snyder's original synthetic unit hydrograph
was based and comment on the differences from the FSR synthetic unitgraph
values.
8.1 Introduction
Civilisation has always developed along rivers, whose presence guaranteed access
to and from the sea coast, irrigation for crops, water supplies for urban com-
munities and latterly power development and industrial water supply. The many
advantages have always been counterbalanced by the dangers of floods and, in
the past, levees or flood banks were built along many major rivers to prevent
inundation in the flood season. In more recent times storage reservoirs have
been built as the principles of dam construction became better understood and
other measures like relief channels, storage basins and channel improvements are
continually under construction in many parts of the world. It is important for
such works that estimates can be made of how the measures proposed will
affect the behaviour of flood waves in rivers so that economic solutions can be
found in particular cases. Flood routing is the description applied to this process.
It is a procedure through which the variation of discharge with time at a point
on a stream channel can be determined by consideration of similar data for a
point upstream. In other words it is a process that shows how a flood wave can
be reduced in magnitude and lengthened in time (attenuated) by the use of stor-
age in the reach between the two points.
198
Flood Routing 199
floods. As many data as possible are required for the second method, which is
the one generally used, including flow records at the beginning and end of the
reach and on any tributary streams joining it, and rainfall records over any areas
contributing direct runoff to it.
Storage in the reach of a river is divided into two parts, prism and wedge
storage, which are illustrated in figure 8.1. This is simply because the slope of
the surface is not uniform during floods (see section 6.2).
wedge storage
A 8
slream bell
reach
If the continuity of flow through the reach shown in figure 8.1 is now
considered, it is clear that what enters the reach at point A must emerge at
point H, or temporarily move into storage.
dS
I=D+-
dt
where I = inflow to the reach
D = discharge from the reach
dS/dt = rate of change in reach storage with respect to time.
This equation is approximated, for a time interval t, by
11-+12
- - t - -Dl
--+D2
-t--82 - 8 1 (8.1)
2 2
where subscripts 1 and 2 denote values at the beginning and end respectively
of the time t. The time t is called the routing period and it must be chosen
sufficiently short so that the assumption impliCit in equation 8.1 (that is, the
inflow and outflow hydrographs consist of a series of straight lines), does not
depart too far from actuality. In particular. if t is too long it is possible to miss
the peak of the inflow curve, so the period should be kept shorter than the
travel time of the flood wave crest through the reach. On the other hand, the
shorter the routing period the greater the amount of computation to be done.
200 Engineering Hydrology
8.3 Reservoir routing
If equation 8.1 is now arranged so that all known terms are on one side, the
expression becomes
(8.2)
The routing process consists of inserting the known values to obtain S2 + 1- D2 t
and then deducing the corresponding value of D2 from the relationship connec-
ting storage and discharge. This method was first developed by L. G. PuIs of the
U.S.Army Corps of Engineers.
The simplest case is that of a reservoir receiving inflow at one end and dis-
charging through a spillway at the other. In such a reservoir it is assumed there
is no wedge storage and that the discharge is a function of the surface elevation,
provided that the spillway arrangements are either free-overflow or gated with
fixed gate openings. Reservoirs with sluices can be treated also as simple reser-
voirs if the sluices are opened to defined openings at specified surface-water
levels, so that an elevation-discharge curve can be drawn. The other required
data are the elevation-storage curve of the reservoir and the inflow hydrograph.
72
70 ~ 1/
V V
E68 V
::66 V V
'"
~64
~62
V V
.60
1/
cAse 1/
V /
56
7
544.0 42 4-4 46 48 50 52 o 20 40 60 80 100
Reservoir surfoce orea (km 2) live slorage (m 3 x 10 6 )
(0) (b)
whose centres are at elevation 54.0 and a free overflow spillway 72.5 m long
with crest level at elevation 66.0. The discharge gates are open and the surface
water level is at elevation 63.5 at time t = o. The flood hydrograph of table 8.2,
column (3), is forecast. What will the maximum reservoir level be and when will
it occur?
Flood Routing 201
(1) Assume the ports have a coefficient of discharge, Cd = 0.8, then Q =
2 (CdA v'(2gH)) and at time 0, Q = D (table 8.2, column (5)) = 2(0.8 x
5.7 xv'186.5) = 125 m 3 /s. (Note g = 9.81 m/s2.) Insert this value on first
line of column (5).
(2) Compute the elevation-storage curve of figure 8.2 (b). Remember that
live storage starts at 52.65, the invert level of the discharge ports, and
amounts to 5.5 x lO6 m3 by 54.0. The storage between 54.0 and 56.0 =
mean area of reservoir between these levels x 2 m = 8.33 x 10 6 m 3 .
Successive increments, computed in this way, are plotted cumulatively
with 52.65 as datum.
(4) From figure 8.2(b) and table 8.1, the D-storage curve of figure 8.3 can
now be drawn; that is, the central curve. The abscissa of figure 8.3 is
graduated in 'storage units'. Each storage unit = routing period x 1m 3 Is.
Since the forecast hydrograph of table 8.2, column (3) is given at 6-h
intervals, it is convenient to make this the routing period. Then each
storage unit = 6 x 3600 x I = 21.6 x 10 3 m 3 = "* m 3 /s day. The use of
these units is necessary to keep the dimensions of table 8.2, columns
(4), (6) and (7) compatible.
202 Engineering Hydrology
(5) The routing calculation can now be started in table 8.2. The figures in
heavy type are known.
To begin with, compute column (4) by averaging successive pairs of
inflow values. Now find the value of this parameter corresponding to
D = 125 m 3 /s from the S - tDt curve of figure 8.3. The value is 2095
and this figure is inserted in the first space of column (6).
(6) The figure in column (4) is added to column (6) and the total put in
column (7) (for example, 2095 + 62 = 2157). The left-hand side of
equation 8.2 has now been evaluated. Find the column (7) value on
the S + t Dt curve and read off the corresponding value of D, entering it
in column (5) (for example, from the S + tDt curve find D = 122 corre-
sponding to 2157).
(7) Use this new value of D to find S - t Dt again as in step (5). Note that
this is directly obtainable without using the curve by subtracting the
value of D from the column (7) value in the line above (for example,
2157 - 122 = 2035). The column (4) figure is then added to the new
column (6) to get a new column (7) figure (for example, 127 + 2035 =
2162).
Flood Routing 203
TABLE 8.2 Reservoir routing computation
(1) (2) (3) (4) (5) (6) (7) (8)
S -.l.Dt S +-tDt Surface
Time Routing Inflow (11 +12 )t D (-!- m2j /s (-!- m 3 /s level
(h) period (m 3 /s) 2 (m 3 /s) day) day) (m)
Complete the table and plot the outflow hydrograph (figure 8.4). The
peak outflow should fall on the recession limb of the inflow graph.
The time difference between the peaks of the inflow and discharge
hydrographs is termed reservoir fag and the reduction in peak flows
together with the spreading out of the recession curve is referred to as
attenuation.
(8) The column (8) values of surface water level are derived from the values
of discharge and levels in table 8.1 ..They are most conveniently found by
plotting a graph and reading off the levels corresponding to the values of
column (5).
The maximum water level in the case of this example is 67.5 m occur-
ring about hour 40.
The solution of the storage equation in this case is more complicated than for
the simple reservoir, since wedge storage is involved. Storage is no longer a func-
204 Engineering Hydrology
-.;
;;,-
E
~ 400
2u
~ 300
"-
~ 200
:E
100
0 12 24 ~ ~ w n ~ % ~ 120 ~
Time (h)
12345678910111213141516171819202122
6-h rouling periods
tion of discharge only, as was the case in example 8.1. McCarthy [1], in what
has become known as the Muskingum method, proposed that storage should be
expressed as a function of both inflow and discharge in the form
S=K[xI+(I-x)D] (8.3)
where x =dimensionless constant for a certain river reach
K = storage constant with dimensions of time that must be found from
observed hydrographs of I and D at both stations.
The two constants can be found as follows. Let figure 8.5 represent the simul-
taneous inflow I and outflow D of a river reach. While I> D, water is entering
Time -+
Figure 8.6 Difference diagram for the river reach of figure 8.5
Time -+
Figure 8.7 Mass curve of storage for the river reach of figure 8.5
Now assume a value of X, say x = 0.1 , and compute the value of the expres-
sion (0.11 + 0.9D) for various times and plot these against corresponding S values
taken from figure 8.7. The resulting plot, known as a storage loop is shown in
figure 8.8(a); clearly there is no linear relationship. Take further values of x
(say 0.2, 03 etc.) until a linear relationship is established, as in figure 8.8(c)
when the particular value of x may be adopted. K is now obtained by measuring
the slope of the line.
Care about units is required. It is often helpful to work in somewhat unusual
units, both to save computation and to keep numbers small. For example,
storage S is conveniently expressed in m3 /s day: such a unit is that quantity
obtained from 1 m3 /s flowing for 1 day = 86.4 x 103 m 3 If S is expressed in
m3 /s day and the ordinate of figure 8.8 is in m3 /s then K is in days.
206 Engineering Hydrology
/)
//j
/ /
/ /
/ /
/ / n
/ /
.,
0> 1/
//
o
o ,/
Vi //
// m
/
/ j
/
(b) K= fir days
----------1(0)
- - - - - - - - - 1 ( b)
---------~(c)
xl + (1- x)D (m 3/s)
The following excerpt from Carter and Godfrey [2] concisely sums up the
choice of values for x and K :
The factor x is chosen so that the indicated storage volume is the same
whether the stage is rising or falling. For spill way discharges from a reservoir,
x may be shown to be zero, because the reservoir stage, and hence the storage,
are uniquely defined by the outflow; hence, the rate of inflow has a negligible
influence on the storage in the reservoir at any time. For uniformly progressive
flow, x equals 0.50, and both the inflow and the outflow are equal in weight. In
this wave no change in shape occurs and the peak discharge remains unaffected.
Thus, the value of x will range from 0 to 0.50 with a value of 0.25 as average for
river reaches.
The factor K has the dimension of time and is the slope of the storage-weighted
discharge relation, which in most flood problems approaches a straight line.
Analysis of many flood waves indicates that the time required for the centre of
mass of the flood wave to pass from the upstream end of the reach to the down-
stream end is equal to the factor K. The time between peaks only approximates
the factor K. Ordinarily, the value of K can be determined with much greater
ease and certainty than that of x.
Having obtained values of K and x, the outflow D from the reach can be
obtained, since by combining and simplifying the two equations
11 -+-
- 12t - Dl + D2 t-- S 2 -
S1 (8.1 )
2 2
and
(the latter being equation 8.3 for a discrete time interval) the equation (8.4)
(8.4)
Flood Routing 207
is obtained, where
A worked example illustrating the application of the method is set out below.
Vi 110
"-
'"$ 100 I
"
I
\ .. ---........
... / \
\
.... <Hypothetical routed outflow
....,
., 90 I \ Discharg~,
01
~
I ....
.c 80 I ....
~
'0 70 II '" ,
"0
c 60 I
'"~ I
I
50
.2 I
i: 40
/
30 ,,-:.- __ ",'
20
10
The data are set out in tabular form in table 83. In columns 1 and 2, the
given hydrographs are listed at a routing period interval, taken as 6 h. The
storage units are taken here as (t m 3 /s day) since the routing period is t day.
Colums (4), (S) and (6) are simply tabular statements of the processes illustrated
in figures 8.6 and 8.7.
A value of x is then chosen, in the first instance 0.2, and the value inside the
square brackets of equation 83 is then evaluated in columns (7), (8) and (9).
Columns (6) and (9) are now plotted in figure 8.10 and produce the loop in the
N
40
00
TABLE 8.3 Storage loops calculations
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15)
Hour Inflow Outflow Mean Cumulotive x =0.2 x = 0.25 x =0.3
I D I-D storage storage
(m 3 /s) (m 3 /s) (m 3 /s) (-! m 3 /s day) (t m 3 /s day) 0.21 O.BD total 0.251 0.75D total 0.31 0.7D total
0 22 22 0 0 0 4 17 21 5 16 21 7 15 22
6 23 21 2 1 1 5 17 22 6 16 22 7 15 22
12 35 21 14 8 9 7 17 24 9 16 25 10 15 25
18 71 26 45 29 38 14 21 35 18 19 37 21 18 39 m
~
24 103 34 69 57 95 20 27 47 26 25 51 31 24 55 IQ
30 111 44 67 68 163 22 35 57 28 33 61 33 31 64 :r
CD
CD
36 109 55 54 60 223 22 44 66 27 41 68 33 38 71 ...
42 100 66 34 44 267 20 53 73 25 49 74 30 46 76 :r
IQ
48 86 75 11 22 289 17 60 77 21 56 77 26 52 78 ::J:
-<
Q.
54 71 82 -11 0 289 14 66 80 18 61 79 21 57 78 ...
60 59 85 -26 -18 271 12 68 80 15 64 79 18 59 77 0
0"
66 47 84 -37 -31 240 9 67 76 12 63 75 14 59 73 IQ
-<
72 39 80 -41 -39 201 8 64 72 10 60 70 11 56 67
78 32 73 -41 -41 160 6 58 64 8 55 63 10 51 61
84 28 64 -36 -38 122 6 51 57 7 48 55 8 45 53
90 24 54 -30 -33 89 5 43 48 6 40 46 7 38 45
96 22 44 -22 -26 63 4 35 39 5 33 38 7 31 38
102 21 36 -15 -18 45 4 29 33 5 27 32 6 25 31
108 20 30 -10 -12 33 4 24 28 5 22 27 6 21 27
114 19 25 - 6 - 8 25 4 20 24 5 19 24 6 17 23
120 19 22 - 3 - 4 21 4 18 22 5 16 21 6 15 21
126 18 19 - 1 - 2 19 4 15 19 4 14 18 5 13 18
Flood Routing 209
:300
280
/ J " J
I1 lA 1
2
240
I1 !
2 20
r-
>-
.g 2 00
1/11240-* I
1(1 = mmtsts day
3
II1
~ I 80
60 3 day
"'E I60
-I" 40
20 X=0.2/) X:0.21 X=O.~
If
00
80
~ I 1//
60
40
20
LJ L,1----1
40 m3 ts Ihll
V [/ :. K = ~ = 15 days 7
O~~~~~Wro80 O~~~~~Wro80
o 10 20 ~ 40 50 60 70 80
xl+(I-x)O (m 3/s)
First, compute Co, Cl and C2 from equation 8.5. Use a routing period t = 6,
h =t day as before.
Co = _ (1.5 X 0.25) - (0.5 X 0.25) = _ 0.25 = -0.2
1.5 - (1.5 X 0.25) + (0.5 x 0.25) 1.25
210 Engineering Hydrology
TABLE 8.4 Derivation of discharge
(1) (2) (3) (4) (5) (6)
Hours I -0.212 0.4 11 0.8 DI D
(rn 3 /s) (rn 3 /s) (rn 3 /s) (m 3 Is) (rn'i'fs)
0 31 31.rf1
6 50 -10.0 12.4 24.8 27.2
12 86 -17.2 20.0 21.8 24.6
18 123 -24.6 34.4 19.7 29.5
24 145 -29.0 49.2 23.6 43.8
30 150 -30.0 58.0 35.0 63.0
36 144 -28.8 60.0 50.4 81.6
42 128 -25.6 57.6 65.3 97.3
48 113 -22.6 51.2 77.8 106.4
54 95 -19.0 45.2 85.2 111.4
60 79 ~15.8 38.0 89.1 111.3
66 65 -13.0 31.6 89.0 107.6
72 55 -11.0 26.0 86.1 101.1
78 46 -9.2 22.0 80.9 93.7
84 40 -8.0 18.4 74.9 85.3
90 35 -7.0 16.0 68.3 77.3
96 31 -6.2 14.0 61.8 69.6
102 27 -5.4 12.4 55.7 62.7
108 25 -5.0 10.8 50.2 56.0
114 24 -4.8 10.0 44.8 50.0
120 23 -4.6 9.6 40.0 45.0
126 22 -4.4 9.2 36.0 40.8
aAssumed value.
The similarly calculated values of Cl = 0.4 and C2 = 0.8 check that -0.2 + 0.4 +
0.8 = 1.0. From equation 8.4
D2 = -0.212 +0.41 1 +0.8D 1
11 ,12 etc. are known from the hypothetical flood hydrograph, but DI is unknown.
Assume a value for DI =/ 1 = 31 m3 /s. This will be very nearly correct since the
river is at a low level and will be near to a steady state. Then the first equation to
be solved is
D2 = -0.2(50) + 0.4(31) + 0.8(31)
= -10.0 + 12.4 + 24.8 =27.2
This value of D2 becomes the DI for the next calculation and the values are
tabulated as in table 8.4.
The outflow hydrograph is plotted as a dashed line to a little way beyond the
peak in figure 8.9.
A disadvantage of the Muskingum method is that it does not include effects
such as friction and diffusion, which can vary widely as a river reaches and then
Flood Routing 211
exceeds bank-full conditions. A method that does include these effects, known
as the diffusion method, was devised by Hayami [3]. Cunge [4], however,
showed that the solution to the Muskingum equations can be made to approxi-
mate to the solutions of the diffusion method provided the values of K and x
are chosen appropriately. Essentially Cunge converted the method from one
derived hydrologically to one based on hydraulic principles. He did this by
writing the original Muskingum equations in finite difference form, identifying
the magnitude of the error introduced and allowing this term to simulate the
diffusion of the flood wave into storage areas along the reach.
In this modified method, now termed the Muskingum-Cunge method, the
reach is divided into a series of sub-reaches, typically L/IO long, x is now derived
for each sub-reach from the physical properties of the reach, i.e. average bed
slope, mean channel width including the flood plain into which the flow diffuses,
the average speed of the flood wave w, and the average peak discharge. K is
assumed to be M/w. The inflow hydrograph is then routed as before, through
the sub-reaches in succession, with the values of K and x varying appropriately in
each sub-reach.
For a more detailed description and further information about this method,
the reader should consult references [5,6].
dS =Kd~
dt dt
and
1-
1-
The principles of flood routing can now be used to derive unit hydrographs for a
catchment where almost no rainfall or runoff records exist. The method is not
entirely synthetic since at least one observation of a runoff hydrograph must be
made.
,
, "-. variable K
'"'"
"'-
" ,
" " "-.
" .............. , ......."-
-~
(a) (b)
1- r. and K
Figure 8.13 Graphical routing with variable lag and K
214 Engineering Hydrology
Consider a catchment as a series of sub-areas, each of which, under a sudden
burst of rain, contributes inflow into the system of drainage channels, which
have storage. The instantaneous unit hydro graph is therefore in two parts, the
first representing the inflow of the rain, and the second the gradual withdrawal
of the catchment storage. The dividing line between these two parts can be con-
veniently taken at the point of inflexion on the recession limb as shown in figure
8.14.
The assumption is now made that the catchment discharge Q and the storage
S are directly proportional, so that
S=KQ (8.6)
(that is, equation 83 with x = 0, and Q used instead of D) and
I _Q= dS
dt
where I represents the inflow resulting from the instantaneous rain. Since
dS/dt =K dQ/dt, by differentiating equation 8.6, then
KdQ =I-Q
dt
and using the condition Q = when t = 0, the equation may be solved to
short rainburst
storage
water
leaving
r---------- .--1------
storage
boseftow
Time_
t= T
T+l
QTe-(t-T)/K
Hence
A = [-KQTe-(t-T)/KJ~
= K(QT - QTe- 1/K )
so
The second observation that must be made from the observed hydrograph is
the catchment lag (Td; that is, the maximum travel time through the catchment.
inflexion point
boseflow
/------------1 -
Example 8.3. Given the catchment area of figure 8.l6, of area 250 km 2 , and
the information derived from a short rain hydrograph that TL = 8 hand K = 7.5 h,
derive the 2 -h unit hydrograph
=8h
K =7-5h
catchment area = 250km 2
(2) Measure by planimeter the area of each of the hourly areas. The areas of
the figure are
Hour 2 3 4 5 6 7 8
10 23 39 43 42 40 35 18
(3) Draw the distribution graph of the runoff using the sub-areas as ordinates
and 1-h intervals as abscissa. The result is figure 8.17- the time-area
graph drawn in full lines.
(4) This time-area graph is now treated as the inflow I due to unit net rain
of 1 cm on the catchment of a hypothetical reservoir, situated at the out-
let, with storage equal to that of the catchment. Then
I +1 t _
Q 1 +Q2 . t =S2 -SI
__
1 __
2 (from equation 8.1)
2 2
and
.......
50r
-'
occounl for non-uniform
oreal dlslribulion of rain
o
Cl>
4 r--
20- .
1--'
10-
Time (h)
where
O.St O.5t K - O.St
mo=--- ml=--- m2=----
K + O.5t K+ O.5t K + O.5t
and since a distribution graph is being used and I I = 12 , then
Q2 =m/I + m2QI
where
/ t
m = ----
K + O.St
and in this case
/ 1
m = =-=0.125
7.5+0.5 8
and
(5) Tabulate the data and compute Q2 as in table 8.5. Q2 is the required
synthetic instantaneous unit hydrograph. Compute the conversion con-
stant for column (3).
10 6 X 10- 2
1 cm rain on 1 km 2 in 1 h = = 2.78 m 3 /s
3600
(6) Plot the IUH and 2-h unitgraphs of table 8.5, columns (5) and (6) as
figure 8.18.
0 0 0 0 0 0
1 10 3.5 0 3.5
2 23 8.0 3.1 11.1 5.5
3 39 13.5 9.7 23.2 13.4
4 43 14.9 20.3 35.2 23.1
5 42 14.6 30.8 45.4 34.3
6 40 13.9 39.6 53.5 44.3
7 35 12.1 46.8 58.9 52.2
8 18 6.2 51.4 57.6 55.5
9 0 0 50.5 50.5 54.7
10 0 0 44.1 44.1 50.8
11 0 0 39.6 39.6 45.1
12 0 0 34.6 34.6 39.3
13 0 0 30.2 30.2 34.9
14 0 0 26.4 26.4 30.5
15 0 0 etc. etc. etc.
V IUH
60
/e-...
/ /
,/
\" / " 2-h unitgroph
"'-. r"",
50
I [/ \.,
.",", "
/ I
~ 40
'"S V/
2'"'" 30 / /
u
<J>
Cl
i
I /
/ "-
r~
//
20
10 I
e~
/1
.,?/
o 10 12 14 16
2 4 6 8
Time (h)
Further reading
Problems
8.1 A catchment can be divided into ten sub-areas by isochrones in the manner
shown in the table below, the catchment lag TL being 10 h:
Flood Routing 221
2 345 6 7 8 9 10
14 30 84 107 121 95 70 55 35 20
A single flood recording is available from which the storage coefficient K is
found as 8 h. Derive the 2-h unit hydrograph for the catchment.
8.2 Tabulated below is the inflow I to a river reach where the storage constants
are K = 10 h and x =0:
Time I Time I
(h) (m3js) (h) (m 3js)
0 28.3 40 90.6
5 26.9 45 70.8
10 24.1 50 53.8
15 62.3 55 42.5
20 133.1 60 34.0
25 172.7 65 28.3
30 152.9 70 24.1
35 121.8
Find graphically the outflow peak in time and magnitude. What would be the
effect of making x > 07 Assume outflow at hour 11 is 28.3 m 3 /s and starting to
rise.
8.3 A storm over the catchment shown in the figure generates simultaneously
at A and B the hydrograph listed below:
~
Hours Q (m 3 js) Hours
B
C
Q(m 3js}
-
0 10 24 91
3 35 27 69
6 96 30 54
9 163 33 41
12 204 36 33
15 210 39 27
18 190 42 24
21 129
Use the Muskingum stream flow routing technique to determine the combined
maximum discharge at C. The travel time for the mass centre of the flood
between A and C is 9 h and the factor x = 0.33. Any local inflow is neglected.
222 Engineering Hydrology
8.4 Define the instantaneous unit hydrograph of a catchment area, and describe
how it can be used to derive the n-h unitgraph.
A catchment area is 400 km 2 in total and is made up of the sub-areas bounded
by the isochrones tabulated below:
Sub-area Area
bounded by (km 2 )
isochrone
(h)
15
2 30
3 50
4 75
5 80
6 60
7 45
8 25
9 20
8.5 Listed below is the storm inflow hydrograph for a full reservoir that has an
uncontrolled spillway for releasing flood waters:
3-h 3-h
intervals m 3 /s intervals m 3/s
0 l.5 12 54
1 156 13 45
2 255 14 40
3 212 15 34
4 184 16 28
5 158 17 23
6 136 18 17
7 116 19 11
8 99 20 8.5
9 85 21 5.5
10 74 22 3.0
11 62
Flood Routing 223
Determine the outflow hydrograph for the 48-h period after the start of
the storm. Assume outflow is 1 m3 /s at time O. The storage and outflow
characteristics of the reservoir and spillway are tabulated below:
0 40 18 220
2 52 20 180
4 132 22 150
6 218 24 120
8 294 26 96
10 322 28 80
12 324 30 68
14 308 32 60
16 264 34 54
224 Engineering Hydrology
8.7 The following hydrograph was observed on a river as a result of an hour
of uniform-intensity rain, baseflow having been estimated and subtracted:
The catchment contributing this flow was divided into sub-areas by isochrones,
from which the following time-area correlation was derived:
Hour 1 2 3 4 5 6 7 8 9 10 11
Area (km 2 ) 25 40 80 200 340 300 220 170 110 50 20
Given that Q2 =m'I + m2Ql where m'= tl(K +-!t)and ml = (K --!t)/(K +-tt)
derive the I-h unit hydrograph for the catchment and determine Qp and tp_
Estimate the net rainfall that caused the original hydrograph.
9.1 Introduction
In the previous chapters the various physical processes involved in the hydrologic
cycle have been enumerated and examined in detail. Methods of evaluating each
process have been suggested and often explained, and techniques discussed that
can be used to provide quantitative answers to many questions.
The remaining problem that must now be tackled is how to use this know-
ledge to predict from existing data, however meagre it may be, what will happen
in future. This is a fundamental problem of all engineering design, since the
engineer designs and constructs work to provide for future needs, whether he
be a structural engineer designing an office block, an electrical engineer designing
power systems to meet future electrical demand, or a hydraulic engineer designing
reservoirs to meet future demand for water.
There is one major difference in these three cases. The structural designer is
working with homogeneous materials whose behaviour is known within narrow
limits. His buildings will be used by people whose spacing, dimensions, weight
and behaviour, en masse, can be predicted quite accurately. He has to cope with
natural events only in the form of wind loads, which usually form a small pro
portion of the total load, and earthquakes. For both these eventualities there are
codes of practice and recommendations available to him.
The electrical system designer has to extrapolate the rising demand curve of
recent years, and to examine the trends of industry and personal habits, to
decide how much capacity should be available in future years. While this is a
complex and continuous task, it is almost completely immune from natural
events other than disasters for which he cannot be expected to provide.
The hydraulic engineer on the other hand is dealing, in reservoir design,
almost exclusively with natural events: in the incidence of precipitation, evapora
tion and so on. These events are usually random in nature and may have any or
all the nonnegative values. It is true that if the rainfall at a place is measured
daily for a period of time, a knowledge about what is a probable daily rainfall
225
226 Engineering Hydrology
will be built up, but it will not, however long it goes on, lead to any limiting
possible value of daily rainfall, other than intuitively.
The hydrologist is frequently asked what the maximum possible discharge
of a particular river will be. There is simply no such value. The only answer that
can be given is that from the data available, and making various assumptions, it
would appear that a certain value will not be exceeded on average more than
once in a specific number of years. On such estimates all hydrologic design must
be performed, and this chapter deals with methods whereby some of the uncer-
tainties can be removed or narrowed in range.
9.2.1 Catchment area fonnulae. The particular random variable of river flood
discharge has been of interest to engineers and hydrologists from the earliest
days of hydrology and many formulae have been proposed to define the 'maxi-
mum flood' that could occur for a particular catchment. The formulae are
empirical by nature, derived from observed floods on particular catchments and
usually of the form
Q=CA n
where Q = flood discharge in m 3 /s (or ft. 3 /s)
A = catchment area in km 2 (or mile 2 )
n = an index usually between 0.5 and 1.25
C = a coefficient depending on climate, catchment and units.
An early example of such a formula, due to Dickens, was developed in India
Q = 825ao.75
with Q in ft. 3 /s and a in square miles; but since the formula takes no account of
soil moisture, rainfall, slope, altitude etc., it is clearly of very little value in
general application. This is true of all such formulae although they are frequently
used to obtain a quick first estimate of the order of 'maximum flood' that can
be expected. For such purposes Morgan [1] proposed the formula for a cata-
strophic flood in Scotland and Wales of
Q = 30001l1l.5
where Q is in ft. 3 /s and M is catchment area in square miles, and added the
sophistication of a recurrence period T (in years) by quoting
1
design flood = catastrophic flood x (T/500)3
for cases where the adoption of the catastrophic flood was not justified by
danger to human life or the safety of a dam. A similar formula of the same
type, due to Fuller, has been widely used in the USA:
Qav = CAO.8
Hydrological Forecasting 227
where A is catchment area in square miles
C is a coefficient often taken as 75
Qav is average value of annual flood discharge in ft. 3 /s.
The value of Qav is then substituted in the formula
Qm = Qav(1 + 0.8 log T)
where T is a return period in years and Qm is the 'most probable' annual maxi-
mum flood.
9.2.2 The Rational Method. The introduction of rainfall into a formula might
be expected to improve it, bearing in mind the type of relationship which exists
between rainfall and runoff (figure 6.7 refers).
This kind of direct relationship of runoff to rainfall depths has been used in
the past to determine flood discharges. Mulvaney [2] was the fust to propose
the idea in his work on Irish arterial drainage. It was also the basis of the Uoyd-
Davis method of sewer designs [3] and the Bransby-Williams estimating method
for floods in India [4]. Its use has persisted to the present because of its sim-
plicity.
The formulae are all of the form
Qp = CiA
where = rainfall intensity in a time t
A = catchment area
le = time of concentration, i.e. the time taken for rain falling on the
catchment farthest from the gauging station to arrive there
C = a dimensionless runoff coefficient, whose value depends on the
catchment characteristics
Qp = peak discharge due to the particular rainstorm, and assumed to
occur after time te when the whole catchment area is contributing
(it follows that the rain is assumed to be uniform over the catch-
ment and to last at least for t c).
Bransby-Williams gave a formula for the design rainfall duration D (in hours)
L 5
D=t =--{/-
e d h
2
) (a
where L = greatest distance from the edge of catchment to the outfall
d =is the diameter of a circle of area equal to catchment area
(Lld is, therefore, a dimensionless coefficient of circularity)
a = catchment area in square miles
h = the channel slope (as a percentage) along its greatest length
te =time of concentration in hours.
228 Engineering Hydrology
If i is in in./h, and A is in acres, then Qp is in ft3/S. Similarly if i is in mm/h
and A in km 2 with Qp in m 3 /s, a correcting factor of 0.278 must be introduced
to permit the same C values to be used, Le.
Qp (fe /s) = C i (in./h) A (acres)
i.e. Qp (m 3 /s) = 0.278 C i(mm/h) A (km 2 )
Many values are quoted for C, from 0.1 to 0.9 depending on the nature of the
catchment, but as the runoff depends also on the intensity and duration of the
storm, the catchment wetness etc., the method is of limited value unless data
about a particular catchment are available. Then it may be used to interpolate or
marginally extrapolate the data.
There is, however, one specialised use of the Rational Method where the value
of C is sometimes taken as 1.0 Le. completely impermeable ground. This is in
the design of urban sewerage which is discussed in chapter 10.
A comparison of the Rational Method with those of the Flood Studies Report
[5] indicated that the Rational Method typically yields results about twice the
flood peak of the FSR for lowland catchments and much more as catchments
become smaller and flatter.
9.3.1 Series of events. The next approach is to use the methods of statistics to
extend the available data and hence predict the likely frequency of occurrence
of natural events. Given adequate records, statistical methods will show that
floods of certain magnitudes may, on average, be expected annually, every 10
years, every 100 years and so on. It is important to realise that these extensions
Hydrological Forecasting 229
are only as valid as the data used. It may be queried whether any method of
extrapolation to 100 years is worth a great deal when it is based on (say) 30
years of record. Still more does this apply to the '1000 year flood' and similar
estimates.
Another point for emphasis is the non-cyclical nature of random events (see
also section 9.6). The lOO-year flood (that is, the flood that will occur on
average, once in 100 years) may occur next year, or not for 200 years or may be
exceeded several times in the next 100 years. The accuracy of estimation of the
value of the (say) lOO-year flood depends on how long the record is and, for
flood flows, one is fortunate to have records longer than 30 years. Notwith-
standing these warnings, frequency analysis can be of great value in the inter-
pretation and assessment of events such as flood flows and the risks of their
occurrence in specific time periods.
It is particularly important to define what is meant by an event. For example,
if a river has been gauged every day for 10 years, there will be about 3650
observations. These are not independent random events since the flow on any
one day is dependent to some extent on that of the day before, and so the
observations do not comprise an independent series. The array of these obser-
vations is termed a full series.
Suppose the maximum event was extracted from the 10-year record in each
year. These woulJ constitute an independent series since it is highly unlikely
that the maximum flow of one year is affected by that of a previous year. Even
so, care is necessary, as can be seen from figure 9.1, where water years, measured
between seasons of minimum flow, are marked, as well as calendar years. One
water years
-;- .. : -
Q
*
I 1
calendar y e 0 r s
9.3.2 Probability of the N-year event. The term recurrence interval (also called
the return period), denoted by Tn is the time that, on average, elapses between
two events that equal or exceed a particular level. Putting it another way, the
N-year event, the event that is expected to be equalled or exceeded, on average,
every N years, has a recurrence interval, TT of N years.
As mentioned previously there is no implication that the N-year event occurs
cyclically. It does, however, have a probability of occurrence in any particular
period under consideration.
Let the probability P(X';;;; x) represent the probability that x will not be
equalled or exceeded in a certain period of time.
Then P(X';;;; x)n will represent the probability that x will not be equalled
or exceeded in n such periods.
For an independent series and from the multiple probability rule
P(X';;;;x)n [P(X';;;;x)]n
= [l -P(X~x)]n
Therefore
1 - [1 -P(X~x)]n
Now
T. = 1
r P(X~x)
Hydrological Forecasting 231
therefore
So, for example, the probability of X ~ x, where x is the value of a flood with a
return period of 20 years, occurring in a particular 3-year period is
P(X ~ 20 yr floodh = 1 - [1 - toP
=1 - [0.9SP
= 1 - 0.857
= 0.143 or 14.3 per cent
Table 9.1 shows the probability of the N-year flood occurring in a particular
period.
For example, it can be seen from the table there is a 1 per cent chance of the
200-year flood occurring in the next 2 years and an 8 per cent chance that it
will not occur for the next 500 years.
If the probability P(X :s;;; x)n is defined by a policy ruling, the value of n, the
design period, can be found from
P(X~x)n = 1-(1 - ;J n
1 - P(X ~ x)n = 1-
_1) n
( Tr1)n (T. i;- = - r
n=
Example 9.1. How long may a cofferdam remain in a river, with an even chance
of not being overtopped, if it is designed to be secure against a lO-year flood?
Here, the policy ruling is that there should be an even chance, so P(X ~x)n = 0.50
and Tr = 10 then
9.3.3 Detennining the magnitude of the N -year event by plotting. Having listed
a series of events (for example, maximum floods) they are then each accorded a
ranking m, starting with m = 1 for the highest value, m = 2 for the next highest
and so on in descending order. The recurrence interval Tr can now be computed
from one of a number of formulae, which have been reviewed by Cunnane [7].
That most frequently used in the past is the Weibull formula [8]
T.=n+l (9.1 )
r
m
where m = event ranking and n = number of events, but there are objections to
its use because of the bias it introduces to the largest events in a short series.
Other formulae used are the Californian [9]
n
Tr = -
m
and Hazen's [10]
2n
Tr= - - -
2m -1
about both of which there are reservations. One of the more satisfactory, due to
Gringorten [11] is
_ (n+O.l2)
Tr - - ------- (9.2)
(m - 0.44)
Hydrological Forecasting 233
For a single, simple compromise, Cunnane recommends
n +0.2
Tr = ------
m -0.4
The probability P of an N-year event of return period Tr is
100
P= - per cent (9.3)
Tr
so that once the series has been ranked, its various events can be plotted on
graphs connecting the variable Q and either Tr or P.
It is often assumed that such series are normally distributed, in which case the
plotted points on normal probability paper would lie on a straight line. This
seldom happens for flood series and shallow curves more often result, making
extrapolations more difficult. To overcome this difficulty the variate Q is some-
times plotted logarithmically, which requires logarithmic-normal probability,
or log-normal paper.
Table 9.2 is a listing of the annual maximum mean daily flows of the River
Thames at Teddington Weir for the years 1883-1988. This is a true annual series
with return periods and probabilities calculated from equations 9.2 and 9.3.
Table 9.2 data can be plotted in a variety of ways and figures 9.2 to 9.5 illustrate
the most common.
Return Percentage
Water Q'r ax Rank period probability
year (m Is) m Tr (year) P Log Qmax
-
--
1000
~
~
800 I ~
r------ l---
J..
<5
600
/~ .-<
- l.------
7I
400
200
!
I
,
o
o 20 40 60 80 '00 ,20 '40 ,60 180 200
Return period, Tt (year)
Figure 9.2 Annual maximum mean daily flow of the River Thames at
Teddington, 1883-1988
(b) Q (linear) against Tr (logarithmic) (figure 9.3). This yields a straight line
fitted to all but the lowest values. Although extrapolation is simpler, unless
Tr follows a logarithmic law, extrapolation is not necessarily more accurate
than for figure 9.2.
1200
T
T
1100
1000 1 1
1 1
900 __ f-
800 1 1
700 I l.---"
t.-- I--'
--:.-:r
~ 600 1
oS
c5 500 ..
i-.-.'
V'"......
400
300
200
.?Y
"
lOO 1 11
0
3 4 5
I
6 7 8 910 15 20 30 40 50 60708090100 150 200 300
Return period, T, (years)
Figure 9.3 Annual maximum mean daily flow of the River Thames at
Teddington, 1883-1988 (semi-log)
Hydrological Forecasting 237
(c) Q (linear) against probability (per cent) (figure 9.4). As often happens,
flood series points lie on a shallow curve on probability paper (where a
normal distribution of probability is assumed).
1200
1100
1000
900
800 I V
I I /
~ 700
'1 I I /'
0- 600
500 I
i
I /.
400
/-
.' I I
'.,,,f"- I
300
~I
Figure 9.4 Annual maximum mean daily flow of the River Thames at
Teddington, 1883-1988 (normal probability)
(d) Q (logarithmic) against probability (per cent) (figure 9.5). The curve of
figure 9.4 is now transformed to a straight line.
A variation of the approach in figure 9.4 is to assume that the logarithm
of the variate Q is normally distributed, leading to the use of logarithmic
normal distribution (or log-normal paper)-first used by Whlpple [12].
(e) Other investigators have proposed methods assuming other theoretical
frequency distributions. Gurnbel [13, 14, 15] used extreme-value theory
(EVl) to show that in a series of extreme values Xl, X 2 Xn where the
samples are of equal size and X is an exponentially distributed variable (for
example, the maximum discharge observed in a year's gauge readings), then
the cumulative probability pi that any of the n values will be less than a
particular value X (of return period T) approaches the value
pi == e-e-Y
1000
900 I I
800 I/'
700 V
600 I I I 1.--:
500 I I 1.:/1
1 400 V" I I I
C .l7 I? ~'''I I
300
vi-'"
200 I :;P V I I I
VV
.....-- .......-
lOO
./"
Figure 9.5 Annual maximum mean daily flow of the River Thames at
Teddington, 1883-1988 (log-normal)
T= __1_
1 _.p '
(Note that this argument refers to Gumbel's method. The reader should not
confuse this with the normal usage of Tr = liP where P = probability of
occurrence.)
The event X, of return period Tyears, is now defined as QT, and
QT = Qav + 0(0.78y -- 0.45) (9.4)
where Qav = average of all values of 'annual flood' Qmax
o = standard deviation of the series.
Thus
(9.5)
Powell [16] suggested that if plotting paper is prepared in which the hori-
zontallines are spaced linearly and the vertical lines' spacing is made proportional
to y, then from equation 9.4 QT and T will plot as straight lines. This is the basis
Hvdrological Forecasting 239
TABLE 9.3 y as function of T
T Y T Y
1.01 -1.53 100 4.60
1.58 0.00 200 5.30
2.0 0.37 300 5.70
5.0 1.50 400 5.99
10.0 2.25 500 6.21
20 2.97 1000 6.91
50 3.90 10,000 9.21
1100
I
I
I
lOO 0
I
1
I
90 0
1
I Q~oo _
1 1/
:1
800
V
.' I
VV
70 0
I I
I
1
V
lA:.
I 1
0
I
I
I
500
1
400 1 1
1
I
., .'
/'
v y..
1 1 III 1
w V-
I
c
300 1 11
V I
I
I
200
:' ;,.: I I 11
V- I
/ I
I
:
lOO 11
I
o 11
1.01 1.05 1.1 1.21.31.41.5 4 5 6 7 8 910 15 20 30 40 5060708090 lOO 150 200
Return period, T, (years)
Figure 9.6 Annual maximum mean daily flow of the River Thames at
Teddington. 1883-1988 (Gumbel distribution)
240 Engineering Hydrology
n> 50, when Qav at 2.33 years is included on the line through the points. The
other point Q200, represents the '200-year flood' and is found by inserting the
appropriate values in equation 9.4.
The correspondence between the plotted data and Gumbel's theoretical line is
demonstrated. Gumbel paper should not be used for partial series, which usually
plot better on semi-log paper as used in figure 9.3.
From the plots presented in figures 9.2 to 9.6, it may seem there is little to
choose between the particular plotting papers available. This is very often the
case and investigators should use whichever distribution makes their particular
job of fitting and extrapolation simplest and the line apparently of best fit.
The foregoing is a necessarily brief resume of the methods of plotting flood
events, in current use. For the underlying theories the reader should refer to
original papers and more comprehensive treatments available [5,7,17,18].
Example 9.1. Determine by calculation the mean daily discharge of the River
Thames at Teddington with a lOO-year return period, assuming the annual
series is normally distributed.
From table 9.2: QBV = 329.7 m 3 /s and a = 133.8 m 3 /s. For Tr = 100 years,
Pis 1.0 per cent; and from table 9.4, K = 2.33. Therefore
Probability K Probability K
of exceedance of exceedance
(per cent) (per cent)
9.3.5 Log Pearson Type III distribution. It has already been mentioned that
annual flood series are rarely normally distributed. A histogram of such series
is usually skewed; that is, the mean value does not coincide with the mode (the
value of the variate with largest frequency). Pearson [19] devised a measure of
skewness based on (mean mode)/o and developed a family of curves to des-
cribe degrees of skewness. One of these, the Pearson Type III distribution, when
used together with the logarithm of the variate Q is found to allow many annual
flood series to plot as straight lines. This log Pearson Type III distribution has
been adopted as a standard by several Federal agencies in the USA [20].
N-year events can be calculated in a similar manner to the normal distribution
method, with this time the additional complication of using logarithms of the
variate and a skew coefficient (G), given by
G = fn 2 (~X3l=- 3n(~X)(~~2~2(~X)3J (9.6)
[ n(n - 1)(n - 2)oi
where X = log Q and Ox = standard deviation on n values of X.
Accordingly, to compute a particular return period flood Q from an annual
series the following steps are required.
(a) Transform all (n) values of Q in the series to their logarithms (base 10):
X = log Q
(b) Find the mean of all values of X:
X
av
= ~X
n
242 Engineering Hydrology
(c) Compute the standard deviation of nvalues of X
Ox = J[n : 1 (L~2 -X;v ) J
TABLE 9.5 Values of K' in the Pearson Type III distribution (21J
9.3.6 Estimation of the mean annual flood Q from a partial duration or peaks
over threshold (POT) series [22]. When only a limited amount of data is
available (say 3 to 10 years) a series of peaks over some arbitrary value may be
used instead of an annual maxima series to estimate the value of Q, the mean
annual flood.
The procedure is as follows.
and
Q = qo + SIn X + 0.5772~
The theoretical background and justification for this is given in reference
[5] .
9.4.1 The general equation for Q. This is a method for obtaining estimates of
flood discharge for ungauged catchments through the use of catchment charac
teristics. The characteristic discharge adopted is the mean annual flood, Q. After
consideration of more than 500 catchments in Britain and Ireland and a detailed
examination of the influence of catchment and meteorological characteristics on
Q, the recommended general equation for various hydrometric areas is
Q = C [AREA0.94 STMFRQo.27 SOIU .23 RSMD I 03 SI 085.16 (1 + LAKE)-0.8S]
The notations are as follows. (9.8)
For the Essex, Lee and Thames area only, a different equation is used
12 =0.373 AREAo. 7o STMFROo.s2 (1 + URBAN)2.s (9.9)
where URBAN is the urban fraction of the catchment.
By using equation 9.8 (or 9.9) and the data provided in maps and figures, it is
possible to predict Q for any catchment in the British Isles. Having obtained a
value of Q for an ungauged catchment, an estimate of a flood QT, where T is the
return period of a magnitude chosen for the design problem posed, is then
required.
Figure 9.8 gives dimensionless region curves of Q/Q plotted against the
reduced variate y from the Gumbel EVl distribution. Return period T is also
scaled. By identifying the particular region in which the catchment lies and using
the appropriate region curve, the factor Q/Q is obtained for any particular
return period and hence QT can be obtained. The regions are delineated and the
values of coefficient C of equation 9.8 are given in figure 9.7.
Example 9.4. Estimate the flood frequency curve for the River Wyre at St.
Michaels. This catchment is in Region 10.
/
5-0
1/6/7
4-0
17/
/17/
2
1/4
1//1/;1
1
1.-8
hVA ~ V.
3-0
9
It-IO ./
Q/Q
2-0 ~~ ~ Irish
~~
~
1-0
~
, ......
~~
Figure 9.8 Region curves showing average distribution of Q/Q for each
region (derived from Flood Studies Report, Figure 1.2.14)
From figure 9.8 the Q/Q factors for Region 10 are listed and the flood frequency
curve tabulated:
Tr Q/Q QT
(years) (m 3 /s)
2 0.9 88
5 1.17 115
10 1.36 134
25 1.63 160
100 2.10 206
500 2.62 258
9.4.2 Small catchment e4.uations for (j. Since the publication of the Flood
Studies Report, various authors have considered the possibility of predicting
248 Engineering Hydrology
Q for small catchments of less than 20 km 2 area with an equation of fewer
terms than the six-term general equation (equation 9.8).
Many small catchments are not gauged and engineers designing for small
bridges and culverts often have to resort to empirical equations. However, when
catchments are small, the STMFRQ and S1085 terms are often difficult to
quantify from standard maps.
Poots and Cochrane [23] studied 42 small catchments throughout the British
Isles and concluded that the equation
Q = 0.0136 AREAo.866 RSMD 1A13 SOILl.S21
gave better results statistically than equation 9.8. Poots later concluded [24]
that a slightly different version was a little better:
Q = 0.015 AREAo.S82 RSMD 1A62 SOIL1.904
A similar exercise was carried out by the Institute of Hydrology o~ 47 catch-
ments of similar size [25]. The equations resulting from regressing Q on catch-
ment characteristics were for three terms and four terms respectively:
Q = 0.00066 AREAo.92 SAAR1.22 SOIL2.o (9.10)
and
Q = 0.0288 AREAo.90 RSMD1.23 SOILl. 77 STMFRQo.23 (9.11 )
In this case it was concluded that there was little difference between the predic-
tion accuracy of equations 9.1 0 and 9.11 and the standard six-term equation
9.8, and that the former two should be used only if time was critical.
References
Further reading
Problems
9.2 The annual precipitation data for Edinburgh are given below for the years
1948-1963 inclusive:
Year Precipita tion Year Precipitation
(in.) (in.)
1948 36.37 1956 28.17
Hydrological Forecasting 253
1949 28.01 1957 25.68
1950 28.88 1958 29.51
1951 30.98 1959 18.04
1952 24.41 1960 24.38
1953 23.64 1961 25.29
1954 35.15 1962 25.84
1955 18.08 1963 30.24
(i) estimate the maximum annual rainfall that might be expected in a 20-year
period and a 50-year period;
(ii) define the likelihood of the 20-year maximum being equalled or exceeded
in the 9 years since 1963.
9.3 Discuss the methods in common use for plotting the frequency of flood
discharge in rivers. List the separate steps to be taken to predict, for a particular
river, the flood discharge with the probability of occurrence of 0.005 in any
year. Assume 50 years of stage observations and a number of well-recorded
discharge measurements with simultaneous slope 0 bservations.
9.4 The following table lists in order of magnitude the largest recorded mean
daily discharges of a river with a drainage area of 12 560 km2 .
Year Date Discharge Year Date Discharge
(m 3 /s) (m 3 /s)
1948 29 May 2804 1922 19 May 1716
1948 22 May 2450 1925 20 May 1694
1933 10 June 2305 1924 13 May 1668
1928 26 May 2042 1917 9 June 1609
1932 14 May 2042 1916 19 June 1586
1933 4 June 2016 1912 21 May 1563
1917 17 June 1997 1918 5 May 1495
1947 8 May 1980 1918 10 June 1495
1917 30 May 1974 1929 24 May 1492
1921 20 May 1974 1943 29 May 1478
1927 8 June 1943 1922 26 May 1476
1928 9 May 1861 1919 23 May 1473
1927 17 May 1818 1936 10 April 1433
1917 15 May 1801 1936 5 May 1410
1938 19 April 1796 1923 26 May 1405
1936 15 May 1790 1927 28 April 1314
1922 6 June 1767 1939 4 May 1314
1932 21 May 1762 1934 25 April 1300
1912 20 May 1753 1945 6 May 1257
1938 28 May 1722 1935 24 May 1246
(continued overleaf)
254 Engineering Hydrology
Year Date Discharge Year Date Discharge
(m 3ls) (/Il3 Is)
9.5 The annual rainfall in inches for Woodhead Reservoir for the period of
record 1921-1960 is listed below:
Year Rainfall Year Rainfall
1921 44.48 1941 46.79
1922 56.25 1942 43.38
1923 65.57 1943 45.87
1924 45.72 1944 59.00
1925 44.78 1945 42.74
1926 48.02 1946 55.39
1927 54.48 1947 41.04
1928 51.53 1948 45.03
1929 48.11 1949 44.11
1930 59.03 1950 52.15
1931 60.59 1951 55.43
1932 47.20 1952 48.87
1933 38.16 1953 43.26
1934 45.38 1954 63.13
1935 54.62 1955 36.46
1936 53.03 1956 56.57
i937 40.98 1957 51.79
1938 50.88 1958 51.22
1939 49.95 1959 39.58
1940 46.27 1960 60.66
The mean and standard deviation for the period are 49.69 in. and 7.08 in.
respectively.
Hydrological Forecasting 255
Arrange the data in ranking order. Compute return periods and probability.
Plot the data on probability paper.
(a) What are the 50-year and lOO-year annual rainfalls? How do these com-
pare with prediction made using Gumbel's theory? What qualification needs to
be made if using the latter results?
(b) What is the probability that the 20-year rainfall will be exceeded in a
lO-year, a 20-year and a 40-year period?
(c) A certain waterworks plant is to be designed for a useful life of 50 years.
It can tolerate an occasional rainfall year of 70 in. What is the probability that
this amount may occur during the project's life?
9.6 A river is subject to annual flooding, and from the records of the observed
annual maxima over 30 years, the data plot as a straight line on semi-log paper,
with return period plotted on the logarithmic scale. The maximum recorded
flood is 900 m3 /s and the minimum is ISO m3 /s. If 1000 m 3 /s is selected as a
design discharge, what is the probability of its being exceeded during the next
20 years?
9.7 Below are the annual mean daily inflows (in m3 x 10 6 ) to a reservoir
during 20 consecutive years:
7.31 6.90 6.64 5.08 5.37 5.75 7.30 7.22 6.48 5.20
6.38 5.51 5.65 5.82 5.81 6.30 6.53 6.12 6.06 6.07
The mean and standard deviation of the data are 6.175 and 0.68 (x 106 m3 ) res-
pectively
(a) Plot the data on Gumbe1 probability paper. Fit a straight line to the data
and estimate the lOO-year mean daily inflow. Compare plot with analy-
tical value if y = 4.60 for T = 100 years.
(b) What is the probability of this 100 year value occurring in any period of
10 consecutive years?
9.8 Determine the magnitude of the 200-year flood from a catchment in the
United Kingdom Region 2, given the following information:
catchment area 107 km 2
stream frequency 1.32 junctions/km 2
SI085 24.2 m/km
30 per cent of area is soil class 1,35 per cent soil class 2, 20 per cent soil class
4 and IS per cent soil class 5.
5 per cent of the catchment area drains through lakes
2-day MS rainfall is 100 mm at catchment outlet
256 Engineering Hydrology
r = (60 minute M5/2-day MS) x 100 = 20 per cent
s.m.d. = 4 mm
What spillway capacity would you design for, with an earth dam, at this location-
and why?
9.9 Annual flood discharges in a river for a 24-year period are listed below:
32.6 17.0 59.5
22.7 5.7 56.6
11.3 36.8 20.0
34.0 51.0 11.8
93.4 8.5 85.2
3l.l 73.6 3l.l
19.8 12.8 25.0
25.5 25.4 23.7
Rank the events in magnitude and compute their return periods and prob-
abilities. Plot the data. Estimate the discharge of a lOO-year flood.
What is the probability of such a flood occurring in the next 2 years?
9.10 The Flood Studies Report gives the following equation for a mean annual
flood in the central region of England:
Q= 0.0213 AREAo. 94 STMFRQo.27 SOIL1.23 RSMD1.03 (I + LAKE)-O.85
S1085o. 16
(a) Discuss the significance of each term in the equation.
(b) Describe how you would use the equation to determine the flood with a T-
year return period (QT)'
(c) How would you amend QT for the design of a dam spillway? Indicate what
value T might have in this case and whether you consider QT may be used
directly.
(d) Indicate the relative merits of this method of flood determination and the
unitgraph-rainfall approach. Give examples of situations where each method
might be used.
9.11 Using the data listed in 9.4, determine by use of the log Pearson III dis-
tribution, the maximum annual mean daily discharge, with a return period of
50 years.
MSL = 8 km no urban
S1085 = 5.0 m/km development
RSMD=75 mm
and baseflow rises from 2 m 3 /s at the beginning of rainfall to 5.5 m 3 /s at its end,
thereafter declining at the same rate. Assume response runoff starts simul-
taneously with rainfall.
9.14 The maximum annual discharges for a river are listed below for thirty-three
years of observation. The mean and standard deviation for the period are
131.5 m 3 /s and 83.0 m 3 /s respectively. Rank the data, compute probability and
plot the data on log-normal probability or semi-log paper.
258 Engineering Hydrology
Year Qpeak Year Qpeak
(m 3 /s) (m 3 /s)
(a) Plot the line of best fit to the data and predict the lOO-year and 400-year
flood discharges.
(b) Given Gumbel's equation QT = Qav + a(0.78y - 0.45) and thaty = 4.60
for T = 100 year and y = 5.99 for T = 400 year, compare the results from
Gumbel's theory with your plot.
(c) What probability is there of a flow equalling or exceeding 300 m 3 Is in the
next five years?
9.15 The following table lists annual peak floods for a river gauging station.
Order the events, compute probabilities and plot the results on log-normal
probability paper.
U se T = -n -+0.2
-
r m _ 0.4
9.16 covers material from Chapters 1-9, and is included as an example of the
multiple-choice type.
D
is the method where 'equilibrium flow' is
used as a check on unit hydrograph integrity?
D reduces evaporation?
10.1 Introduction
Figure 10.1 Key plan of sewer system (reproduced from TRRL Road
Note 35 (1976))
system. The longest sewer (measured from the outfall) is numbered 1 and
individual pipe lengths numbered successively 1.0, 1.1, 1.2, etc., starting
at the upstream end. The first branch sewer which joins 1 is numbered 2
and the longest sewer contributing to the branch has its component parts
numbered 2.0, 2.1, 2.2, etc. This process is repeated until every pipe in
the system is numbered.
(b) A table is then prepared, of which table 10.1 is a typical example. The
data of columns 1, 2, 3, 4, 9, 10 and 11 are then entered in it. All of this
information should be available from a site survey and building proposals
and plans. Notice that it is only the directly-contributing areas which are
used. Permeable areas of gardens, woods, etc. are not included.
(c) Column 12 is then completed by adding in branch contributions at
junctions as the table is completed. Column 12 gives the total surface
area contributing to a particular length of sewer.
(d) It is' now necessary to fmd the time of concentration, i.e. the time it takes
for rain falling on the small area of impermeable surface that contributes
to pipe 1.0 to arrive at the manhole at its downstream end. To do this it
is necessary to assume a particular pipe diameter and material for 1.0. An
entry time must be chosen and 2 minutes is recommended for normal
urban areas, increasing to 4 minutes for very large, flattish paved areas.
Knowing the pipe diameter and gradient, the full-bore velocity in the
TABLE 10.1 Rationalformula design sheet,b
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Pipe Difference Time of Time of Rate of lmpermea b1e area ha Rate of Pipe
length in level Length Gradient Velocity flow concentration rainfall flow diameter
number (m) (m) (1 in) (mM (min.) (min.) (mm/h) Buildings, Total (l/s) (mm) Remarks
Roads yards etc. (9 + 10) Cumulative
1.0 1.J0 63.1 57 1.33 0.79 2.79 67.9 0.089 0.053 0.142 0.142 26.8 150
...c:er
I
1.1 1.J2 66.1 59 1.70 0.65 3.44 62.5 0.077 0.109 0.186 0.328 56.9 225
1.2 0.73 84.7 116 1.46 0.97 4.41 57A 0.081 0 0.081 0.409 65.2 300 =
%
2.0 1.40 44.8 32 2.32 0.32 2.32 72.5 0.113 0.081 0.194 0.194 39.1 225 -a
2.1 0.61 49.1 80 1.77 0.46 2.78 67.9 0.045 0.105 0.150 0.344 64.9 300 ...0
3.0 0.98 48.5 49 1.43 0.56 2.56 70.5 0 0.129 0.129 0.129 25.2 150 0"
2.2 1.65 54.3 33 2.74 0.33 3.11 65.5 0.101 0.073 0.174 0.647 117.7 300 Sewer 3 added ca
l.3 1.22 27.7 23 3.29 0.14 4.55 56.9 0.121 0.235 0.356 1.412 223.2 300 Sewers 2 and 3 -<
4.0 0.88 54.9 62 1.66 0.55 2.55 70.5 0.093 0.093 0.186 0.186 36.4 225 added
4.1 0.58 45.7 79 1.48 0.52 3.07 66.3 0.069 0.040 0.109 0.295 54.2 225
1.4 052 22.9 43 2.77 0.14 4.69 55.8 0.069 0 0.069 1.776 275.3 375 Sewer 4 added
m
264 Engineering Hydrology
pipe may be found from published tables [3] using the Colebrook-White
formula and an appropriate k value in that formula, for the pipe material.
This value is entered in column 5 and with the pipe length of column 3,
the time of flow may be calculated and entered in column 6.
(e) The rate of rainfall corresponding to the time of concentration in column
7 (time of entry and time of flow in the case of the first pipe) is then
determined from Meteorological Office tables. These are based on the
Meteorological section of the Flood Studies Report and may be obtained
for any National Grid reference point (approximately in the centre of the
drainage area) from the Meteorological Office, London Road, Bracknell,
Berkshire. Note that it is necessary to choose a design rain frequency.
Suitable frequencies may be between one year for modern separate
surface water systems and 100 years for old combined systems with base-
ment developments. Table 10.2 is an abbreviated typical example.
TABLE 10.2 Rates of rainfall in mmjh for a range of duration and return
period for a specified location in the UK (in Southern England)a
References
Further reading
Since the Flood Studies Report was published, many of the techniques devised
for it have been subsequently developed and their range of application extended.
One of the investigations arising from this development was the World Flood
Study [1] which, in the words of its authors, "was conceived with the aim of
examining and classifying the characteristics of floods in as many countries and
from as wide a range of climates as possible." Much of this work has been
subsequently reported in the literature [2].
Data were collected from individuals and organisations around the world in
70 countries, and were assessed and analysed uniformly to allow comparisons of
results between stations and countries. Both annual maximum flood series and
catchment characteristics were collected.
The results were used to produce growth curves (see figure 9.8) showing the
increasing value of the parameter Q/Q with decreasing frequency of occurrence.
By making the curves non-dimensional in this way, results could be compared by
region, climate, catchment area and annual rai~fall, for homogeneity.
It was, of course, necessary to establish Q, the mean annual flood (termed
MAF in [1, 2]) from the data set available in each region or country, but the
growth curves derived can be used to estimate floods at sites where there are no
gauging stations, always provided that Q may be estimated from catchment
characteristics.
A selection of the growth curves is shown in figures 11.1 to 11.4. It is empha-
sised that these represent averages of broad regions and, for specific projects and
detailed design, more data would always be necessary. The curves presented in
the figures are only some of those available. For example, Thailand on figure
11.4 has five curves representing different altitudes and catchment areas, the one
shown as no. 5 being an average curve for the country, excluding the Malay
peninsula. Those countries marked with an asterisk also have other more rigidly-
defined catchment curves presented in [1] . North and South America feature in
the source publication but are not sampled here.
268
International Flood Frequency Growth Curves 269
Interested readers are recommended to consult the original source publica-
tions for details of the data used, the statistical techniques employed, and the
limitations and qualifications about the use of the curves.
6,--.------"----,---.-----"---.-----.------.---~
1 Romania*
2 Bulgaria
3 Czechoslovakia"
5 4 Poland"
5 Spain
6 Hungary & Yugoslavia"
7 Denmark"
8 Sweden & Norway'"
9 France. Germany.
Netherlands
10 Italy'"
4~~----~~--~----+_----~--~+_--V_--~~~~
6'-~----~,---,---,,----,---~-,-,------,---,
24-~------~---M~--+---~~--~~~+-----+---4
9
10
11
6~-.-------r----.----r-----.----.---.-----.----'
aT
Cl
2~-+------~---H~~A---~~--~----+-----+---~
6,,-r------.----.----.------.--~----._----._--~
1 Eastern Australia
2 Australia"
3 Thailand. Malaysia."
(Malay peninsula)
5 __ 4 Japan"
5 Thailand"
6 Korea
7 Papua New Guinea
8 Indonesia'"
Certain trends are evident from the material shown here. For example, semi-arid
climates have the steepest growth curves since the infrequent heavy rainfall
produces floods relatively much greater than the mean annual flood. Examples
of this in figure 11.2 are the curves for South Africa, Algeria, Morocco and
Tunisia, and in figure 11.3, Saudia Arabia and Jordan. Similarly, the equatorial
climates with high annual rainfall tend to have flatter curves, e.g. the Nigeria,
Ghana, Congo group in figure 11.2 and Indonesia and Papua New Guinea in
figure 11.4.
Clearly, the similarity of climate and topography are not the only factors
influencing the growth curves since Sri Lanka and Kerala State in India are both
subject to the South West monsoon but have appreciably different curves.
Other important factors are area, soils and type of vegetation. The larger the
catchment area, the flatter the curve is a general rule, though the authors of [1]
mention one exception to this in the Malay peninsula.
The Institute of Hydrology in the UK is continuing to expand the database of
flood series and catchment characteristics, and invites interested engineers and
hydrologists throughout the world to send it such information.
References
l. MEIGH, J. and FARQUHARSON, F. A. K. World Flood Study, Phase II.
Institute of Hydrology, Wallingford, UK, Nov. 1985
2. FARQUHARSON, F. A. K. et al. Comparison of flood frequency curves for
many different regions of the World, from Regional Flood Frequency
Analysis (V. P. Singh, ed.), Reidel, Dordrecht, 1987
Further reading
All hydraulic engineering design carries with it an implied structural life during
which the structure is expected to meet its design specification. Some structures
have specified design lives (for example, a cofferdam or a diversion tunnel) but
often such a life is not made explicit. How long, for example, are the design lives
of a canal, a dam, or a breakwater?
The engineering designer has to approach the problem of design of hydraulic
structures from consideration of the damage that would ensue if failure of the
structure occurred and/or if it is possible to exceed the design specification with-
out structural failure. It is not sufficient to consider return period alone; what is
needed is an understanding of the risk of encountering certain conditions during
specific periods of time, and of the consequences of the design flow being
exceeded, including the danger to life and the economic, environmental and
social effects of failure of the structure. Only then can a proper judgement be
formed as to whether or not enough has been done to make the risk acceptable
without economic repercussions that are unacceptable.
It may be difficult to make any quantitative assessment of these factors, but
consideration of a series of steps may at least provide a logical means of com-
parison of the various options open to the designer:
(a) Identification of the events or sequences of events which may lead to
failure, and determination of their probability of occurrence.
(b) Identification of specific features of the structure which might initiate
failure or partial failure (e.g. failure of gates to operate or be operated,
loss of power, etc.) and estimates of their probability.
(c) The likelihood of combinations of events in (a) and (b).
(d) The consequences, including the forecast economic, social and environ-
mental costs of each combination in (c).
The thoroughness with which these studies are made will depend on the
nature of the structure. Dams are the most immediately obvious risks, but canals
274
Design Criteria 275
and their control structures, river bank revetments and sea defence works are all
hydraulic structures liable to need such risk analysis.
In the end, there must be engineering judgement about the acceptable level of
risk. Such judgement can best be exercised using whatever quantification of the
risk can be made, and it should always be made clear to the owner/client con-
cerned what is the basis for a recommended level of installation. A much fuller
discussion of the subject has been published in references [1] and [2].
The discussion of damage, and how it can be described and assessed in this
context, is beyond the scope of this text, but has been studied by Borgman [3],
US Army Corps of Engineers [4] and Penning-Rowsell and Chatterton [5].
12.2 Choice of design return period by consideration of design life and probability
of encountering design flow during this life
In table 9.1, the probability of the N-year flood occurring in a specific period of
years is set out. Table 12.1 is a somewhat expanded version of table 9.1, and a
different terminology is employed to make the concept clear. The return period
of an event (Tr) is related to the arbitrarily chosen design life (L) to give the
probability of encountering the event, or something greater, in that design life
(P). The table is based on the expression
P= 1- (1 - ~JL
TABLE 12.1 Percentage risk of encountering an event within a
particular design life L for different return periods a
1 20 10 5 3 2 1
2 36 19 10 7 4 2 1
3 49 27 14 10 6 3 1
5 67 41 23 16 10 5 2 1
7 79 52 30 21 13 7 3 1 1
10 89 65 40 29 18 10 5 2 1
15 96 79 54 40 26 14 7 3 1
20 99 88 64 49 33 18 IQ 4 2
30 96 78 64 45 26 14 6 3
50 99 92 82 64 39 22 9 5
75 98 92 78 53 31 14 7
100 99 97 87 63 39 18 10
150 99 95 78 53 26 14
200 98 87 63 33 18
300 95 78 45 26
500 99 92 63 39
1000 99 86 63
aWhere no figure appears the risk is either less than 0.5 per cent or more than 99.5 per cent
276 Engineering Hydrology
12.3 Choice of a design value of a rare event
Many of the techniques described in this text were developed to provide the
designer with numerical values of discharge, rainfall etc. that have physical
meaning (that is, they yield velocity, depth, volume and so on).
The various numbers that result from calculations, however, will not all
have equal significance because they have been derived in different ways, and
care needs to be taken before they are used as specifications. Consider the follow-
ing three examples.
References
Further reading
MOSONYI, E. and BUCK, W. Some aspects of risk analysis for improved water
resources planning. Proc. 2nd World Congress Int. Wat. Res. Assoc. New Delhi
1975, Vol. 2, p. 221
Appendix A Rainfall and Soil
Characteristics of the
British Isles
The British and Irish National Grids are illustrated. The rainfall and soil type
maps are at a scale such that there are ten separate sections covering Britain
(with the exception of some Scottish islands) and four covering Ireland.
Each section has national grid lines identified by numbers signifying hundreds
of kilometres from the origin and is itself identified by a large bold number (for
example, southeast England is 3). For convenience, four different maps of each
section are kept together. These are:
Each map has one of these abbreviated titles followed by a bold number that
identifies the section on the key map. Britain is dealt with first and then the
Irish key map and Irish sections follow.
These maps are reproduced from the NERC Flood Studies Report, 1975,
except for RP which is from the 1978 revision of Winter Rain Acceptance
Potential, published as Supplementary Report No. 7.
278
Appendix A 279
1300
HL HM HN HO HP JL
1200
<U.
HQ HR HS HT JQ
I ~u
1100 I-n
HV HW HX HY: HZ JV
: NA NE OA
900L
800
I OF
OL
OQ
600
SOO
400
SK
300
SP
200
SQ
100
SV
0
/,0
"
100 200 lOO 400 700 k~~:::)
SOO 600
FAl SE ORIGI .... OF
NATIONAL GRIO
----------~~~S~A~A~R.:1-------------~
282 Engineering Hydrology
2 3
3+-----------------------+-------------.r--~--~~
2DM5.1
Appendix A 283
r. 1
284 Engineering Hydrology
RP.1
Appendix A 285
3 4
3,
SAAR.2
286 Engineering Hydrology
2DM5.2
Appendix A 287
3 4
3 .".., .........
!O
C~ 40
("'------'-'-":
..... It;
-.. __ .. ;'
" -.....---
~" ...", ../........ i
c)
10
......
". ____ .1
o
\ .......... ~
3I ., ,
r.2
288 Engineering Hydrology
3 .
, 5
I
RP.2
Appendix A 289
,
SAAR.3
290 Engineering Hydrology
5 6
. . "'~
.\
.....4f'; i<
~
J.
FIII"NeE
........
... .... ,
A.- \
G ~-.
...
.
2DM5.3
Appendix A 291
5 6
(
31~~~--+-----~~----~----~----+---------~--~-3
................. .
'.'.
-""
......
..... ---_ .. -
_-4S,,-
t6
-._-,,-
........
....... ...l
t .... 45"
2 2
--_.')
....._-_ ..
'RANCE
ChcaJMl 1.1o.1Id1
UTM
,.3
292 Engineering Hydrology
s
RP.3
Appendix A 293
.L---+-7R.lI}--l 5
SAAR.4
294 Engineering Hydrology
2 3
5,~----------------------~------~-'~----------~
4+-----------------------4------------------------+4
3+-----------------------+-------------~~~~~~
2t-------------------~~1---L-----~--~~~~~~
2
3
2DM5.4
Appendix A 295
2 3
5i+-------------------------1--------,~,,~~--------_t5
4+---------------------------~--------------------------_t,4
3;+-------------------------~--------------~-e~----~~
{:#~
:. .......?
~----------------------------~~~~~~--~--~------~\~~--~~--+Q
2 3
r.4
296 Engineering Hydrology
2 3
RP.4
Appendix A 297
SAAR.5
298 Engineering Hydrology
3 4
~.....
'''''.)
3 4
2DM5.5
Appendix A 299
4 5
0"/0,,'
{ !
\ ......
.'
.... __ .......':
(~......... -".
r.5
300 Engineering Hydrology
RP.5
Appendix A 301
S~~Y--- 5
I:" .. - -
., . .
3
:' \S' ,.;,\ ;"
...
2~~46~ SAAR.6 2
302 Engineering Hydrology
6
5~~>0~~------------------------+---------------~
4-------,.i \ \ ~~~~~-------------+----------------4
, ~
';
.'-...
7
5 6
2DM5.6
Appendix A 303
5 6
5~----.-+'------------------------+-----------------
J(. . . . /?,.
:
~-
41------~+_-r~~--~._------------+_----------------,~
.... ~6-._~
.....
3~~~--~----~~----~'----~~---+--------~.-~---3
.......... _-_. .....
\'"
)
._-45'-'
------":
.,;
.'
5 6
r.6
304 Engineering Hydrology
5 6
( 5
"
2
5 6
RP.6
Appendix A 305
2 3
SAAR.7
306 Engineering Hydrology
5.1---+2 ---.4.::::::I:liL-_ r
, 2DM5.7
Appendix A 307
2 3
""'"", I .15---)
':
.. ,.-
13 ; _, -IS")
"," "'--
! "
,"
\ r
'"
r.7
308 Engineering Hydrology
2 3
RP.7
Appendix A 309
SAAR.8
310 Engineering Hydrology
3 4 5
2DM5.8
Appendix A 311
3 5
81~~------~
....~.-y----~--------------------~
:'~"""'25-""'~i
i'';''
........ .....
:"",_,:
7't---~~~~------~------------------k
.. _.....
'S'"2S.}
t, .......
"! . . . .
\/
~
. ........l>\
3 4 5
r.8
312 Engineering Hydrology
3 4 5
RP.S
Appendix A 313
SAAR.9
314 Engineering Hydrology
2DM5.9
Appendix A 315
2 3
2 3
f.9
316 Engineering Hydrology
2 3
RP.9
Appendix A 317
'7
5
SAAR.10
318 Engineering Hydrology
2DM5.10
Appendix A 319
i
5
91+-------------------------+-------------------------~
,/"'/ ("
f
' .. -"~
..../.
7+---~~--~~~----------4_------------------------~
3 4 5
r.10
320 Engineering Hydrology
I I
S ~'8 4 5
4 5
RP.10
Appendix A 321
A B
4
H
3
N
2
S
1
X Y
o 1 2 3 4
I SAAR.'
324 Engineering Hydrology
- - --~ 12DM5.1
Appendix A 325
2------------------~~~~----_r_r--------~~~~
25
I r.1
326 Engineering Hydrology
2
I RP.1
Appendix A 327
2
()
I SAAR.2
328 Engineering Hydrology
2 3
I~~==~-----,~~~~~~----t----------------l
..-75
2 3
I2DM5.2
Appendix A 329
2
3----~------------~~------~------_+_+------------------3
.. , ....
" .... "
,
,:'
,'",
2~
I r.2
330 Engineering Hydrology
~~----~~------~P/~\T~rt---------------3
I RP.2
Appendix A 331
I SAAR.3
332 Engineering Hydrology
4 - - - - - - - -________-+____________ -.~~~~~
,.
()
ot
.......')
/
(~-~~~---.'
~D iD('
D
2--------------------~~~------~--~~~~~~~--~
I \
2
12DM5.3
Appendix A 333
4------------------~------------~~ __~~--4.
"
2--------------------~~~~~----~~----------~~__+~
2
I r.3
334 Engineering Hydrology
I RP.3
Appendix A 335
()
2---t---L~~~--~--L---~~~~1l-------2
2 3
I SAAR.4
336 Engineering Hydrology
12DM5.4
Appendix A 337
[)
~. !:'------,
..
{.... .3S .... ', ..
I
.
3----r-~=-------~~----~----_+_+---------------3
2 3
I r.4
338 Engineering Hydrology
2 3
3~~-----*--~~--------~~~+--4-4------------------
2 3
I RP.4
Appendix B Typical Values of Manning's
n in Q = (1/n) (AR 213 S1I2)
and Chezy's C* in
V = C-j(RS)
Type of channel n C
(SI units)
Smooth timber 0.011
Cement-asbestos pipes, welded steel 0.012 70-90
Concrete-lined (high-quality formwork) 0.013 60-75
Brickwork well-laid and flush-jointed 0.014
Concrete and cast iron pipes 0.015
Rolled earth: brickwork in poor condition 0.018 40-55
Rough-dressed-stone paved, without sharp bends 0.021 30-45
Natural stream channel, flowing smoothly in clean
conditions 0.030 19-30
Standard natural stream or river in stable
condition 0.035 14-25
River with shallows and meanders and noticeable
aquatic growth 0.Q45
River or stream with rocks and stones, shallow
and weedy 0.060
Slow flowing meandering river with pools, slight
rapids, very weedy and overgrown 0.100
* For a full discussion of Chezy's coefficient C see An Introduction to Engineering
Fluid Mechanics by J. A. Fox, published by The Macmillan Press, London, second
edition, 1977.
339
Answers to Problems
Chapter 2
2.1 (a) 28.32 mm Hg; (b) 8.5 mm Hg; (c) 19.82 mm Hg; (d) 22.0C;
(e) 23.7C.
2.3 (a) 29.99 in.; (b) Cubley 29.76 in.; Biggin School 41.97 in.
2.4 39.1 mm.
2.5 (i) 29.8 in.; (ii) 31.2 in.; (iii) 27.0 in.
2.7 (b) About 1951; (c) assuming earlier period correct, increases it from 279
to 328 mm/year.
2.8 Some evidence of cyclicity, but longer record needed to show it conclu-
sively.
2.9 (a) 0.136; (b) 7.34 m/so
2.10 (a) Once in 4 years on average; (b) 23.8 x 10 3 m3 .
2.11 At X, 75 mm; at Y, 92 mm.
2.12 (a) 20 mm Oxford; (b) 64 mm Kumasi (if equation 2.4 derived for U.K.
applies in Ghana).
2.13 481 mm; no.
Chapter 3
3.1 Amsterdam 4.0 mm/day; Seattle 0.05 mm/day; Khartoum 6.2 mm/day.
3.2 2.5 mm/day.
3.3 (a) April 5.24 cm; November 1.01 cm; (b) June 11.67 cm; October
3.70 cm.
3.5 6.7 mm/day.
3.9 Annual Eo = 966 mm, pan coefficient = 0.85, net annual loss = 1.367 x
103 m3 /km2 /day, evaporation 136 mm less in July at 40 0 S.
3.10 14.076 x 106 /m 3
340
Answers to Problems 341
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9