4 5985348647112411029
4 5985348647112411029
4 5985348647112411029
HYDROLOGY
This page is intentionally left blank
ELEMENTARY ENGINEERING
HYDROLOGY
M. J. DEODHAR
Emeritus Fellow
AICTE, New Delhi
India
This book is sold subject to the condition that it shall not, by way of trade or otherwise, be lent, resold,
hired out, or otherwise circulated without the publisher’s prior written consent in any form of binding or
cover other than that in which it is published and without a similar condition including this condition
being imposed on the subsequent purchaser and without limiting the rights under copyright reserved
above, no part of this publication may be reproduced, stored in or introduced into a retrieval system, or
transmitted in any form or by any means (electronic, mechanical, photocopying, recording or otherwise),
without the prior written permission of both the copyright owner and the above-mentioned publisher of
this book.
ISBN 978-81-317-0805-7
First Impression
Published by Dorling Kindersley (India) Pvt. Ltd, licencees of Pearson Education in South Asia.
Preface ix
About the Author x
1. Introduction 1
2. Hydrometeorology 16
3. Evaporation and Transpiration 25
4. Precipitation 51
5. Infiltration 98
6. Hydrograph 125
7. Runoff 177
8. Floods 209
9. Discharge Measurement 236
10. Flood Routing 284
11. Groundwater 315
Appendices 346
Bibliography 387
Index 389
CONTENTS
2. Hydrometeorology 16 5. Infiltration 98
2.1 Definition 16 5.1 Definition 98
2.2 Vapour Pressure 18 5.2 Process of Infiltration 99
2.3 Relative uHmidity 18 5.3 Factors Affecting nfi I ltration 99
2.4 Wind 19 5.4 Measurement of Infiltration 101
2.5 Temperature 20 5.5 Expression of Infiltration 114
2.6 Air Mass 21 eRview uQestions 121
eRview uQestions 22 uNmerical uQestions 122
Multiple Choice Questions 23 Multiple Choice Questions 123
Professor M. J. Deodhar graduated in 1957 and completed his postgraduate studies in Hydraulics
and Dam Engineering in 1959 from the University of Pune. He started his career as Assistant
Research Officer in Maharashtra Engineering Research Institute, Nasik, and later became Professor
and Head of the Civil Engineering Department in the Government College of Engineering at Karad,
Amaravati and Pune. He retired as Principal of K. K. Wagh College of Engineering, Nasik.
He has been Visiting Professor and also Emeritus Fellow, AICTE, New Delhi. He has worked as
technical expert for UPSC, MPSC, UP PSC and Himachal Pradesh PSC, and as Chairman, Board of
Studies in Civil Engineering. He has also been a member of Faculty of Engineering and Technology
and the Academic Council of universities in Maharashtra.
He has published technical articles and also presented papers at the international conferences,
the prominent being the International Association of Hydraulic Engineering and Research held at
San Francisco and Beijing. He has also conducted and attended short-term courses on behalf of the
Indian Society for Technical Education, New Delhi, and UNESCO.
With over 50 years of research and teaching experience behind him, Professor Deodhar
was awarded the ‘best teacher’ award and ‘outstanding engineer’ award by the Institution of
Engineers, Nasik Local Centre. Presently, he is working as Technical Advisor for Coastal Power
Consultants, Pune.
Introduction
1
Chapter Outline
g 1.1 DEFINITION
Hydrology today denotes a field of science that covers many more branches of science than it did
earlier. Till 1906, hydrology referred to the study of underground water only, as distinct from the
present usage that covers both groundwater as well as surface water. In a much broader sense,
hydrology refers to the study of water.
The term ‘hydrology’ is derived from two Greek words hydor and logas meaning ‘water’ and
‘science’, respectively. So, in simple terms, hydrology is a science related to water.
The Merriam-Webster dictionary describes hydrology as a science dealing with the properties,
distribution and circulation of water on the surface of land, in the soil, in the underlying rocks and in
the atmosphere, particularly with respect to evaporation and precipitation.
2 ELEMENTARY ENGINEERING HYDROLOGY
C. O. Wisler and E. F. Brater define it as The science that deals with the processes governing the
depletion and replenishment of the water resources of the land areas of the earth.
The ad-hoc panel on hydrology of the Federal Council for Science and Technology (established
by the president of the United States of America in 1959) recommended the following definition
for hydrology:
Hydrology is the science that treats waters on the earth, their occurrence, circulation,
distribution, their chemical and physical properties and their environment including their
relation to living things.
In short, what happens to the rain is the basis of the definition of the science of hydrology. It
should not be confused with hydraulics, which deals with the mechanics of water. Hydrology has
got a variety of practical applications. Therefore, it should not be treated as a pure science.
The total quantity of water available on the earth is estimated as 1348.25305 106 km3, and
this may cover the earth to an average depth of 2.73 km—assuming earth a uniform sphere of
12,800 km in diameter.
Man understood some hydrological phenomena such as precipitation, runoff, and so on. Studies
were done related to these in the prehistoric times also, and were based on logical thinking, but
as a science, hydrology was developed only recently. Those studies and investigation of nature’s
processes were initially the tasks of philosophers, and were explained by them based on the assump-
tions and deductions therefrom.
4 ELEMENTARY ENGINEERING HYDROLOGY
Serial no. Water occurrence Estimated water (km3) Estimated water (%)
A. Salt water
1 Oceans 1307.410 106
2 Salt water lakes and island seas 0.100 106
Total sea water 1307.510 106 96.9781
B. Freshwater
3 Glaciers, Polar ice caps 30.4300 106
4 Atmospheric moisture 0.0140 106
5 Hydrated earth minerals 0.0040 106
6 Water content in plants and animals 0.0011 106
B1. Non-utilizable freshwater 30.4491 106 2.2584
7 Freshwater lakes 0.1246 106
8 Rivers 0.00115 106
B2. Surface resources of freshwater 0.12575 106 0.0093
9 Soil moisture 0.0375 106
10 Groundwater up to 800-m depth 4.998 106
11 Groundwater below 800-m depth 5.6309 106
B3. Groundwater resources 10.1682 106 0.7542
B4. Total freshwater resources 40.74305 106 3.0219
Total water resources of the earth 1348.25305 106 100
Men, in earlier days, noticed that rainfall annually occurred occasionally for a short span of
time. They also saw there were some perennial rivers and springs. The problem for the early think-
ers was ‘Wherefrom these rivers and springs received their supply?’ The solution of these problems
was based on the following two assumptions:
1. Rainfall was inadequate to account for all the surface water in rivers and springs.
2. The earth was impervious below a certain depth from the surface.
Homer (1000 BC), Thales (650 BC) and Plato (400 BC) put forth the idea of subterranean flow
from sea to land. By subterranean water, they meant that there is a continuous underground
supply from sea to land, and during the flow the sea water loses its salinity. They also believed
that the subterranean water flow supplied water to springs. Aristotle (350 BC) also assumed that
there was an atmosphere below the ground. According to him, water vapour came from the seas
and condensation took place, and this was the major source for the springs and rivers during
the non-rainy days. Lucretius (90 BC), Pliny (30 BC) and Seneca ( 65 AD) accepted this theory.
Vitruvius, who lived during the time of Christ, had mentioned that the source of rivers and
springs was rain or ice only. But nobody believed this, and the theories proposed by the early
thinkers were believed.
Leonardo da Vinci (1442–1519) and Bernard Palissy (1509–1589) correctly understood the
hydrological process, and proposed the infiltration theory. Vitruvius, Vinci and Palissy are, there-
fore, considered as the pioneers in advocating correct hydrological principles.
INTRODUCTION 5
5. Period of modernization (1800–1900) The experiments were modernized, and the founda-
tion of modern science of hydrology was laid.
6. Period of empiricism (1900–1930) In this period, empiricism in hydrology became more
evident. A number of empirical formulae were suggested.
7. Period of rationalization (1930–1950) In this period, rational analysis of hydrological prob-
lems was used instead of empiricism.
8. Period of theorization (1950–till date) In this period, theoretical approaches have been used
extensively in hydrological problems. Many hydrological principles have been subjected to math-
ematical analysis. Use of sophisticated instruments and computer techniques has been developed.
These phases indicate that hydrology is a young science with many important problems only
imperfectly understood, and research in this field is still going on.
Water consists of two hydrogen and one oxygen molecules strongly bonded together. It is an inor-
ganic liquid, and the only substance occurring naturally as solid, liquid and gas. Water has the
highest surface tension compared to other liquids, except mercury. Some chemical properties of
water are mentioned below:
• Water is the universal solvent and many compounds and salts dissolve in it.
• Pure water boils at 100 °C and freezes at 0 °C.
• Water by itself is a bad conductor of electricity; however, dissolved ions like Ca, Mg, Na, Cl, and so
on, have the ability to conduct electricity.
• Water has the greatest thermal conductivity except mercury.
• The physical properties of water are as follows:
Mass density 1000 kg/m3, weight density 9.81 kN/m3.
INTRODUCTION 7
Atmosphere
Evapotranspiration
Precipitation
Sublimation
Precipitation
Precipitation
Precipitation
Evaporation
Evaporation
Vegetation
Precipitation
Snowpack Through fall
and
icecaps
Melt Surface
runoff
Land surface
Vapor Stroms,
Infiltration Exfiltration
diffusion lakes and
rivers
Soil Interflow
Surface
Precolation capillary rise
runoff
Groundwater
Aquifers flow
Ocean
The cyclic movement of water from sea and back to the sea is known as hydrologic cycle,
hydrological cycle or water cycle. This cyclic order of events may be explained as follows:
The water from the sea evaporates due to solar radiation. It becomes lighter than air and hence
moves up in the atmosphere. Here, under favourable conditions, clouds are formed. These clouds
move upwards and then over the land. Precipitation from the clouds occurs over the land area.
Water, thus, precipitated gets collected in streams and rivers and then flows back to the sea. All
these events are repeated continuously. Thus, the hydrological cycle is a case of circular infinity.
The hydrologic cycle has no beginning. Nobody knows for sure when it started and from which
stage it started. But it is surmised that it started some 3500 million years ago. It is the most distinc-
tive and important process on this planet. It is driven by the energy of the sun, influenced by the
gravity of the earth and the physical, chemical and biological properties of water.
The cyclic order of events does occur in nature, but there are some short cuts also at different
stages. Figure 1.1 shows the engineering representation of the hydrological cycle.
Figure 1.2 shows the descriptive representation of the hydrological cycle.
The circulation of water penetrates the three phases of the earth system:
1. Hydrosphere: Bodies of water that cover the surface of the earth
2. Atmosphere: The gaseous envelop above the hydrosphere
3. Lithosphere: Rocks below the hydrosphere
8 ELEMENTARY ENGINEERING HYDROLOGY
Infiltration
Ground water Interflow Ocean
Table 1.2 The quantitative analysis of the hydrological cycle on the earth
A. Precipitation
1 Precipitation on ocean PO 360 106 346,000 961
2 Precipitation on land PL 150 10 6
99,000 660
Total precipitation P 510 106 445,000 872
B. Evaporation
3 Evaporation from oceans EO 360 106 383,000 1063
4 Evaporation from land EL 150 106 62,000 413
6
Total evaporation E 510 10 445,000 872
C. Runoff R 150 106 37,000 247
P PO PL
E EO EL
R PL EL
The average global precipitation is 872 mm, which is equivalent to 445,000 km3 of water. The
average atmospheric moisture is 14,000 km3. This means that the atmospheric moisture is replaced
32 times in a year, or the residence time of atmosphere moisture is 10 days. Figure 1.3 shows the
disposition of global annual average precipitation.
INTRODUCTION 9
Atmosphere
660 mm 961 mm
413 mm
1063 mm
Land
Ocean 6 2
Ocean
247 mm Ocean
6 2
9. Surface runoff The precipitated water after meeting all the requirements flows in a stream.
This is known as surface runoff or overland flow.
10. Interflow When water infiltrates, it starts moving laterally towards a stream and appears on
the surface. This is known as interflow. It is above the groundwater table. The velocity of flow is
very low as compared to the surface flow. It is also known as sub-surface storm flow, sub-surface
runoff, storm seepage and secondary base flow.
11. Groundwater flow The infiltrated water may reach the saturated zone of water below ground
and may get stored between its pores and voids between particles. This water stored eventually may
flow towards a stream. This happens when the groundwater table is above the stream water level.
However, when the stream water level is above the groundwater table, there may be a flow from
stream to groundwater. This movement of water below the ground is known as groundwater flow.
Any phenomenon that undergoes continuous changes is generally termed as process. All hydro-
logical phenomena change with time and thus they are termed as hydrological processes. These
hydrological processes are complex and involve many variables that are interdependent and can be
classified as follows:
1. Deterministic
2. Stochastic or probabilistic
INTRODUCTION 11
In the study of hydrology, the following modern developments and techniques are used:
• Remote sensing
• Artificial neural network
• Geographical information system
• Genetic algorithm
• Fuzzyl ogic
• No gradient information is required that avoids the mathematical complexity of conventional non-
liner optimization methods.
• GA reduces the computational time during evolution process due to the possibility of parallel
computation.
Following are the key steps in GA:
1. Generating an initial population
2. Describing a coding scheme for all the variables
3. Running the flow simulation finite element model for all the sets of variables
4. Computing the fitness function from objective functions
5. Performing the operations with genetic operators
6. Terminationc ondition
REVIEW QUESTIONS
1. Define hydrology. Also discuss its importance in modern times.
2. Is hydrology a pure science? Discuss.
3. What are the sub-branches of hydrology?
4. State and discuss the supplementary sciences for the study of hydrology.
5. What are the overlapping areas of hydrology?
6. Explain the practical use of hydrology in other fields.
7. State and explain the influence of hydrology in other sciences.
8. Write a note on water budget of the world.
9. Write a detailed note on the history of hydrology.
10. Write a note on hydrology in ancient India.
11. Explain the different phases of development of hydrology as suggested by Prof. Ven Te Chow.
12. Explain with a neat sketch hydrological cycle.
13. State and explain the different processes in a hydrological cycle.
14. Analyse the hydrological cycle quantitatively.
15. What are the water-holding elements in a hydrological cycle?
16. What are the liquid-transport phases of a hydrological cycle?
17. What are the vapour-transport phases of a hydrological cycle?
14 ELEMENTARY ENGINEERING HYDROLOGY
4. Assuming earth to be a uniform sphere, the total quantity of water may cover to an average
depth of
(a) 1.73km (b) 2.73km
(c) 3.73km (d) 4.73km
2
Chapter Outline
g 2.1 DEFINITION
The word atmosphere is derived from two Greek words, atmos and spheria. Atmos means vapour and
spheria means sphere or spherical earth as its extended meaning. Atmosphere, therefore, means the
gaseous envelope surrounding the earth.
Hydrometeorology is the study of that atmospheric process, which affects the hydrologic characters
of a region.
is termed as humidity. Besides water, ozone (O3) and carbon dioxide (CO2), exist in the atmosphere
in varying proportions.
140
Temperature
120 Exposphere
100 Ionosphere
Height in KM
Mesopause
80
Pressure Mesosphere/
60
Chemosphere
Stratopause
40
Stratosphere
20 Tropopause
Troposphere
0
⫺80 ⫺60 ⫺40 ⫺20 0 20 40 60 80 100 120 140 160
Temperature in 0 c
0 1 2 3 4 5 6 7 8 9 10
Pressure in 1000 mb
3. Mesosphere: This is a warm layer and is above the stratosphere and the thickness is roughly 30
km. It is also known as chemosphere. The region at the top of the mesosphere is called mesopause.
4. Ionosphere: This layer is above the mesosphere; the pressure in this layer is very low. It consists
of a number of layers with different electrical characteristics.
5. Exosphere: The ionosphere gradually merges into the outermost shell called exosphere.
Water exists in the form of vapour in the atmosphere. Vapour is lighter than air.
The partial pressure exerted by this vapour is known as vapour pressure. The more the
vapour content, the more will be the vapour pressure. The vapour pressure is normally denoted
by ‘e’. For any temperature, there is a limit for the water vapour that can be held in the atmo-
sphere. For a given temperature, when the air holds the vapour to its maximum capacity, the
vapour pressure is said to be saturated vapour pressure at that temperature and is generally
denoted by ‘es’.
Saturated vapour pressure increases with temperature.
Air with vapour is lighter than dry air.
Vapour pressure is normally expressed in millibars (mb)
1 mb ⫽ 102 N/m2
and 1.33 mb ⫽ 1 mm of mercury.
Relative humidity is the ratio of actual vapour pressure and saturated vapour pressure for that tem-
perature, and is generally denoted by ‘h’ and in terms of percentage.
Naturally, it is always less than 100.
2.3.5 HYGROGRAPH
It is noticed that human or animal hair reacts to the change in the humidity. It expands as humidity
increases and contracts with a decrease in humidity. This variation is very uniform. This prop-
erty is used to measure humidity and is automatically recorded on a graph. This instrument is
known as hygrograph. Sometimes hygrograph is also equipped with a thermograph and is called a
thermohygrograph.
g 2.4 WIND
Air in motion is called wind. The horizontal component parallel to the earth surface is called wind
and the vertical component is referred to as air current.
ooo
g 2.5 TEMPERATURE
The temperature indicates the heat energy in soil, water and air mass. It affects the hydrologic
processes like evaporation, precipitation, infiltration, etc.
Air mass is a vast and deep body of air. The temperature, humidity, etc. are homogeneous in the
air mass. The air mass may move because of the pressure gradient. It can be categorized into the
following two types:
• When the air mass moves over a warmer surface, it is called cold air mass.
• When the air mass moves over a cooler surface, it is called warm air mass.
The cyclone may cover a circular area of 500 km2 and its centre is called eye of the storm. A
cyclone in the northern hemisphere revolves always anticlockwise, while in the southern hemi-
sphere it revolves clockwise. The wind speed may be to the extent of 300 km/h and may pick up
moisture from the land and ocean surface and may cause heavy precipitation.
2.6.5 ANTICYCLONE
When high pressure is developed at a location, the air mass moves away spirally outwards and
then it is called anticyclone. The spiral movement is anticlockwise in the northern hemisphere and
clockwise in the southern hemisphere.
REVIEW QUESTIONS
1. Discuss the permanent and variable constituents of the atmosphere.
2. Describe the vertical structure of the atmosphere.
3. Discuss the factors that affect temperature at a location.
4. Explain the working of an anemometer.
5. What is humidity? How is it measured?
6. Discuss the effect of wind on various hydrological processes.
7. Explain the effect of temperature on various hydrological processes.
8. Write short notes on
a. Vapour pressure b. Saturation deficit
c. Relative humidity d. Measurement of temperature
e. Eye of a cyclone f. Lapse rate
g. Variation of wind with the height h. Dew-point temperature
above ground
i. Wet-bulb temperature j. Measurement of humidity
9. Differentiate between
a. Warm front and cold front b. Average temperature and normal temperature
c. Cold air mass and warm air mass d. Cyclone and anticyclone
e. Tropical cyclone and extratropical f. Dry bulb thermometer and wet-bulb
cyclone thermometer
g. Wind and air current h. Isobar and isotherm
i. Hygrograph and thermograph
HYDROMETEOROLOGY 23
2. A pressure of 1 mb is equal to
(a) 100 N/m2 (b) 10 N/m2
(c) 1 N/m2 (d) 0.1 N/m2
*** ***
Evaporation
and Transpiration
3
Chapter Outline
g 3.1 DEFINITION
Evaporation is an important process in the hydrologic cycle preceding precipitation. It is the pro-
cess by which water in the liquid form transforms into vapour through the transfer of energy. The
exact laws governing this phenomenon are not completely known.
The transformation from solid to the vapour state without passing through the usual intermedi-
ate liquid state is known as sublimation. In the atmosphere, evaporation occurs from free water
surfaces like (1) seas, (2) lakes and rivers, (3) moisture in solid surfaces of (4) land and soil,
(5) vegetation, (6) snowfields, (7) glaciers and (8) even from the falling rain drops.
26 ELEMENTARY ENGINEERING HYDROLOGY
Large amounts are spent in constructing a dam and creating a reservoir to store water. Canals are also
constructed to carry water. Some water stored in these reservoirs and flowing in canals is lost due to
evaporation. The amount of the loss of water due to evaporation can be substantial (about 30%).
It is, therefore, necessary to estimate the loss due to evaporation to assess the water available
for use in different seasons.
The factors affecting evaporation are (1) temperature, (2) solar radiation, (3) wind, (4) barometric
pressure and altitude, (5) dissolved solids, (6) turbidity, depth of water, (7) shape of surface, (8) extent,
i.e., total area of water surface, (9) colour of water, (10) velocity of water, and (11) waves at water
surface. These are discussed below in detail.
3.3.1 TEMPERATURE
Evaporation is directly affected by the temperature of water as well as that of the air above. The
kinetic energy of the molecules of water increases with temperature. Similarly, the saturated vapour
pressure of the air film near the water surface and the vapour pressure of the air above increases.
Since evaporation depends upon the difference of these two vapour pressures, change in tem-
perature may not affect evaporation in the same proportion and, as such, effect of temperature on
evaporation is somewhat complicated.
Evaporation observed in a lake for one complete year indicates that for the same mean temperature
observed in different months, the evaporation figures differ. However, evaporation increases with tem-
perature, even though a close relation between evaporation and temperature cannot be established.
3.3.3 WIND
The air above the water surface receives the evaporated water molecules, and the rate of evapo-
ration depends upon the moisture content in the air. If the water molecules in the air above are
removed quickly, then the rate of evaporation increases. Removal of water molecules in the air
above water surface is done by the wind and naturally depends on the velocity of wind. The more
the velocity of wind, the faster will be the process of removal of water vapour particles and more
will be the evaporation rate.
However, there is a limiting velocity of wind up to which the rate of evaporation increases with
the increase in the wind velocity. Further increase has no effect on the evaporation rate. This higher
limit of wind velocity depends on other factors also.
28 ELEMENTARY ENGINEERING HYDROLOGY
The velocity of wind referred above is naturally near the water surface, because the variation in the
velocity of wind in the vertical direction is not a straight-line relation. Wind speed observed by an ane-
mometer at a higher level will have to be modified considering the variation in the vertical direction.
So also, the velocity of wind will affect the evaporation rate if the air movement is turbulent. The
water vapour will be removed quickly by turbulent diffusion. If the airflow is laminar, it will have no
significant effect on the removal of water vapour and indirectly on the evaporation rate.
125 A
Evaporat on n mm
100
B
75
50
25
0
Jul Aug Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul
Period
The following methods are generally adopted to evaluate the rate of evaporation from a reservoir.
• Water budget method
• Energy budget method
• Mass transfer method
• Actualobs ervations
• Evaporationf ormulae
Example 3.1
The catchment area of an irrigation tank is 70 km2. The constant water spread during October 2006
was 2 km2. During that month, the uniform precipitation over the catchment was recorded to be 100
mm. 50% of the precipitation reaches the tank. The irrigation canal discharges at a uniform rate of
1.00 m3/s in the month of October.
Assuming seepage losses to be 50% of the evaporation losses, find out the daily rate of evapora-
tion for October 2006.
Solution:
100
Total inflow 70 106 _____
1000 0.5 3.50 10 m
6 3
Example 3.2
The surface area of a reservoir in m2 given by A 100 y2, where y is the depth of water in metres
in the reservoir. In one week, the water depth in the reservoir has reduced from 10 to 9 m. Find the
average hourly rate of evaporation. Assume seepage loss to be 40% of the evaporation loss.
EVAPORATION AND TRANSPIRATION 31
Solution:
Therefore, A1 Area at depth 10 m
100 102 10,000 m2
A2 100 92 8100 m2 ________
1 冤A A (A A ) 冥
Average surface area __3 1 2 √ 1 _____________
2
__
1
3 冢10,000 8100 √ 10,000 8100 冣
9033.33 m2
Loss of water 9033.33 1
9033.33 m3
Since seepage loss 0.4 evaporation loss
Therefore, evaporation loss seepage loss 1.4 Evaporation loss
9033.33 103
Therefore, rate of evaporation ____________________
1.4 9033.33 7 24
4.25 mm/h/m2
Example 3.3
The catchment area of a reservoir is 10.0 km2. A uniform precipitation of 0.5 cm/h for 2.0 h was
observed on 7th of July. 50% of the runoff reached the reservoir. A canal carrying a discharge of
1.25 m3/s is taken from the reservoir. The rate of evaporation observed was 0.7 mm/m2/h. The seep-
age loss was observed to be 50% of the evaporation loss. Find the change in the reservoir level on
4th July from 8:00 a.m. to 6:00 p.m. if the water spread of the reservoir was 0.476 km2.
Solution:
0.5
Inflow into the reservoir 10.0 106 ____
100 2.0 0.5 5.0 10 m
4 3
Total outflow from the reservoir 4.5 104 0.5 104 5.0 104m 3
i nflow.
Hence, there will not be any change in the reservoir water level.
Example 3.4
A trapezoidal channel of bed width 4.0 m and side slopes 1:1 carries water at a depth of 2.0 m. The rate
of evaporation observed was 0.35 mm/m2/h. Find the daily loss due to evaporation from the canal in a
length of 10 km in ha m.
Solution:
Top width of water level 4.0 2.0 2.0 8.0 m
Total water spread in the canal in 10.0 km 8.0 10.0 1000 8 104m 2
32 ELEMENTARY ENGINEERING HYDROLOGY
Example 3.5
A reservoir having a water spread of 0.5 km2 and water storage of 1 million m3 receives the Sun’s radi-
ation at the average rate of 0.025 Langley per second. 5% of this radiation is reflected back into the
atmosphere. The rise in the average temperature of the water is 2 °C during 9:00 a.m. to 5:00 p.m.
Assuming no exchange of heat between water and ground as well as the atmosphere, find the
average rate of evaporation.
Solution:
1 Langley 1 cal/cm2
Latent heat of evaporation: 539.55 cal/g
Specific heat of water: 1 cal/g/C
(A) Energy received from the Sun 0.5 106 104 0.025 8 3600
3.6 1012c al
(B) Energy received by water 0.95 3.6 1012
3. 42 1012c al
(C) Energy utilized by water 2 1 106 106c al
2.0 1012c al
EVAPORATION AND TRANSPIRATION 33
Porous sphere 1
2
3
4
Waterproof system 5
6
7
8
9
Air vent 1
Rubber stopper 2
3
Rubber stopper 4
5
6
7
8
Glass tube 9
1
2
3
4
5
6
7
8
9
Bottle
3.4.5.1 Pan
A pan is a metal container, square or circular, and of uniform cross-section. Normally it is circular.
The diameter may range from 300 to 1500 mm.
This pan is filled with water and the loss of water from this pan in a specified period is mea-
sured. The rate of evaporation observed here is correlated to the evaporation from a reservoir.
EVAPORATION AND TRANSPIRATION 35
Example 3.6
The drop in water level in the case of a 1.5-m diameter evaporation pan during 24 h is 10 mm. The
precipitation recorded during this period was 15 mm. Find the rate of evaporation from the pan.
Solution:
( 1.5)2 1.767 m2
Area of pan __
4
EP
Evaporation in 24 h ______
Area
10 15
Therefore, Rate of evaporation __________
1.767 24 0.59 mm/h/m
2
Example 3.7
The amount of water added to a 1.0-m diameter circular land pan over a period of 10 h to maintain
the water level constant is 9.56 lit. Find the hourly rate of evaporation from the pan.
Solution:
D2 __
Area of pan __
( 1)2 0.785 m2
4 4
103 103
________________
9.56
Rate of evaporation per hour 0.785 10
Rate of evaporation 1.21 mm/h/m2
Example 3.8
Observations taken on a 1.0-m diameter circular land pan on 7 July 2006 from 8:00 a.m. to 6:00 p.m.
were as follows:
1. Quantity of water taken out of the pan 5.0 lit
2. Precipitation during this time period 20 mm
If there was no change in the water level in the pan, find the rate of evaporation.
36 ELEMENTARY ENGINEERING HYDROLOGY
Solution:
1. Area of pan __
4 1 0.785 m
2 2
Example 3.9
In the case of a 1.5-m diameter circular pan, following observations were taken from 8:00 a.m. to
6:00 p.m.
1. Quantity of water added to keep the water level in the pan constant is 5.0 lit
2. Precipitation during the observation period is 10 mm
3. Leakage from the pan is 2.0 lit
Find the rate of evaporation from the pan.
Solution:
Area of the pan __
4 1.5 1.767 m
2 2
Example 3.10
The rate of evaporation from a 1.25-m pan having a pan coefficient of 0.8, was 1 mm/m2/h. Find
the total evaporation from a reservoir in a week, having a water spread of 2.0 ha.
Solution:
1
Evaporation from the reservoir during a week _____
1000 0.8 2 10 7 24
4
2688 m3
EVAPORATION AND TRANSPIRATION 37
Example 3.11
The monthly evaporation (in mm) observed on a pan from January to December 2006 is as follows:
111, 126, 127, 132, 141, 146, 148, 143, 138, 126, 118, 114
The water spreads of reservoir in January 2006 and December 2006 were 2.56 and 2.69 km2,
respectively. Assuming a pan coefficient of 0.8, find the total loss of water due to evaporation in
2006 from the reservoir. Neglect all other losses.
Solution:
_______
Mean spread of reservoir __1 A A A A
(
3 1 2 2)
√ 1 __________
1 (2.56 2.69 2.56 2.69 )
__ √
3
2.624 km 2
Annual loss due to evaporation 111 126 127 132 141 146 148 143 138
126 118 1 14 1570 mm
1570
Therefore, evaporation loss in 2006 2.624 106 _____
1000 0.8
3.295 10 m 3 3.295 million m3
6
0.5
0.6
Pan coefficient
0.7
0.8
0.9
1.
0 1 2 3 4 5 6 7 8
Diameter of pan in ‘m’
The standard Floating Pan of 900-mm2 area and 450-mm depth is as shown in Fig. 3.6. The pan
coefficient varies from 0.7 to 0.82.
122 cm
25.4 cm
20 cm
Wooden support
15 cm
G.L
7.5 cm
25.5 cm
Steel stiffener
121 cm
The standard type is made of ungalvanized unpainted iron sheet with an area of 900 mm2 and 450
mm in depth. The rim is kept 100 mm above water level and the water level same as that of ground,
as shown in Fig. 3.7. The pan coefficient varies between 0.75 and 0.86.
Stilling well
10 cm
61 cm
183 cm
The various variables used in the formulae in estimating evaporation from a reservoir are as
follows:
E Evaporation in mm/day
es Saturated vapour pressure in mm Hg
ea Actual vapour pressure in mm Hg
V Wind velocity at the water surface in km/h
Pa Mean barometric pressure in mm Hg
Ta Average air temperature in °C
Tw Average water temperature in °C
The formulae are as follows:
1. Fitzgeraldf ormula
E (0.4 0.124 V )( es ea)
Example 3.12
For a reservoir the following observations were taken:
1. Saturation vapour pressure of air 30 mm Hg
2. Relative humidity 0.5
3. Velocity of wind at the rate of 0.5 m above ground 25 km/h
4. Water spread of reservoir 10 km2
5. Atmospheric pressure 10.4 m of water
6. Average air temperature 29 °C
7. Average water temperature 27 °C
Evaluate monthly evaporation by using different formulae.
EVAPORATION AND TRANSPIRATION 41
Solution:
The monthly evaporation rate was calculated by four empirical formulae.
1. Fitzgerald’ s equation
E (0.4 0.124 V )( es – ea)
Actual vapour pressure of air
__________________________
Relative humidity
Structure vapour pressure of air
ea ___
ea
0.5 __
es 30 since es 30 mm Hg
Therefore, ea 15 mm Hg
Assuming average velocity of wind
V Velocity of wind at the rate of 0.5 m above ground
E (0.4 0.124 V )( es ea)
(0.4 0.124 25) (30 – 15)
(0.4 3.1) (15)
E 52.5 mm/day
Therefore, monthly evaporation loss E 30 water spread of reservoir
52.5
_____
1000 30 10 10
6
( ) 10 0 15
V10 V0 5 ___
0.5
V10 V0 5 1.567
V10 25 1.567
V10 39.18 km/h
E C(es – ea) (1 0.0621 V10)
15 (30 – 15) (1 0.0 621 39.18)
15 15 3.433
772.42 mm/month
Therefore, monthly evaporation loss E water spread of reservoir
772.42
______
1000 10 10
6
7.14 106 m3
7.14 million m3
Evaporation from a reservoir can be reduced by as much as 30%. The different methods to control
evaporation from a reservoir are as follows.
g 3.6 TRANSPIRATION
Plants absorb water from the soil through their roots. Mineral salts are also absorbed in dilute solu-
tion using water as its vehicle. This solution is transported through roots and stems to the leaves
where plant food is produced from this sap, CO2 from atmosphere and using energy from the Sun
through chlorophyll. The plant food thus produced is again distributed using water as a vehicle for
the cell growth and tissue building.
Most of the water thus absorbed by the plants is discharged back into the atmosphere in the form
of vapour. This process is known as transpiration. It is observed that only 1% of water sucked in
by the roots is retained by the plants.
3.6.4 EVAPOTRANSPIRATION
In studying water balance from an area, it is difficult to separate evaporation and transpiration and
hence these two processes are combined and treated as a single and termed as evapotranspiration.
Lysimeter: It consists of a circular tank, its diameter ranging from 600 mm to 3000 mm. It is
buried in the ground so that its top is flush with the ground. It is filled with the soil similar to the
field conditions, and the crop or the tree of which AET is to be found is grown in the lysimeter. A
typical lysimeter is shown in Fig. 3.8.
Water is added to the lysimeter and an account of water added, etc. is kept. Then AET is calcu-
lated as follows:
P W O AET S
where, P Precipitation during the period of observations
W Waters upplied
O Waterdr ained
AET Actuale vapotranspiration
S Change in soil moisture in lysimeter
Lysimeters are expensive and the process of observation is time-consuming.
1.5 to 2 m
A
AREA = A = 1 m 2
Inflow–outflow measurements: This method is applicable to large areas, such as up to 100 ha.
An account of the inflow of water is maintained for a specific period. The inflow may be precipita-
tion or surface inflow or even groundwater inflow. During the same period, an account of outflow
is maintained. The outflow may be the surface outflow or the groundwater outf low. Naturally, the
difference between these two will be the loss due to evapotranspiration.
If the observation period is small, then it is assumed that the groundwater inflow and outflow
are equal. In this method, it is very difficult to measure inflow, outflow and groundwater storage to
the desired accuracy.
Immediately after the rain stops, evaporation from the soil starts. Evaporation from a saturated
soil surface is approximately the same as that from a water surface at that temperature. As the soil
begins to dry, the evaporation goes on reducing. Finally it virtually stops as there is no supply of
water. If the groundwater table is within 1.0 m from the ground, evaporation from soil is noticed
since there is some supply of water.
Evaporation opportunity is defined as ‘the ratio of evaporation from soil and evaporation from
equivalent water surface’. It is dimensionless and is generally expressed in terms of percentage.
Thus,
Evaporation from soil
Evaporation opportunity ___________________________________ 100
Evaporation from equivalent water surface
Evaporation from soil is measured by using a lysimeter of size 1 m 1 m 1 m. The soil
is filled flush with the tank edge and the groundwater table is maintained at a specific level. The
evaporation from the soil can be determined by weighing the tank at known intervals and from the
quantity of water added to maintain the groundwater table.
REVIEW QUESTIONS
1. Define evaporation. Explain the process of evaporation.
2. Explain Dalton’s law of evaporation.
3. Discuss the factors that affect evaporation.
4. How is evaporation from a reservoir estimated?
5. Explain the water budget method of estimation of evaporation from a reservoir.
6. Explain the energy budget method of estimation of evaporation from a reservoir.
7. Explain the mass transfer method of estimation of evaporation from a reservoir.
8. What is an atmometer? Discuss the ones normally used.
9. What is a pan? How are observations taken on a pan?
10. What is a pan coefficient? Discuss the factors that affect pan coefficient.
11. Explain land pan with the help of a neat sketch.
12. Explain sunken pan with the help of a neat sketch.
13. Explain floating pan with the help of a neat sketch.
14. Discuss the merits and demerits of the different types of pans.
15. Discuss the various parameters involved in the different formulae used to estimate evapora-
tion from a reservoir.
16. State and explain the different formulae normally used to estimate evaporation from a reservoir.
EVAPORATION AND TRANSPIRATION 47
17. Can evaporation from a reservoir be reduced? Explain the different methods.
18. Write a detailed note on the use of surface films to reduce evaporation from a reservoir.
19. Explain the process of transpiration. Discuss the different methods to reduce transpiration.
20. Discuss the different methods to measure transpiration.
21. What is evapotranspiration? Explain the different methods to estimate it from a catchment.
22. Discusss oile vaporation.
23. Write a note on the pattern of evaporation over India.
24. Write short notes on
a. Potentiale vapotranspiration. b. Factorsa ffectingt ranspiration.
c. Livingstonea tmometer. d. Pitcha tmometer.
e. Transpirationr atio. f. Phytometer.
g. Potometer. h. Lysimeter.
i. Panc oefficient. j. Factorsa ffectingt ranspiration.
k. Annuale vaporationi nI ndia. l. Atmometer.
m. Evaporationoppor tunity.
25. Differentiatebe tween
a. Evaporation,s ublimation b. Phytometera ndpot ometer.
and transpiration.
c. Land pan, sunken pan d. Evapotranspiration and potential
and floating pan. evapotranspiration.
NUMERICAL QUESTIONS
1. The uniform precipitation over a catchment area covering 80 km2 received a total precipitation
of 80 mm in a month. 60% of the precipitation reached the reservoir. The spread of reservoir
during that month was 2.5 km2. The irrigation canal discharges at a constant rate of 1.1 m3/s. If
the seepage losses were 60% of the evaporation losses, find the daily rate of evaporation.
Ans: 8.24 mm/h/m2
2. In the case of a reservoir, the surface area in m2 is given by A 90 y2, y is the average depth of
water in the reservoir in metres. In a week, the depth of water reduced from 9.0 m to 8.5 m. If
the evaporation loss is double the seepage loss, find the rate of hourly evaporation loss.
Ans: 2 mm/h/m2
3. The bed width of a trapezoidal channel is 5.0 m and the side slopes are 1:1. It carries water
at a depth of 1.8 m. The rate of evaporation observed on a 1.5-m diameter floating pan was
0.25 mm/h/m2. Find the loss of evaporation in a day from the canal in a length of 5 km. Assume
the pan coefficient to be 0.8.
Ans: 0.0206 ha m
4. During a test, the drop in water level in the case of a land pan of diameter 1.0 m in 24 h was
observed to be 15 mm. The precipitation recorded during this period was 10 mm. Find the rate
of evaporation from the pan.
Ans: 1.32 mm/h/m2
48 ELEMENTARY ENGINEERING HYDROLOGY
5. A test was conducted on a 0.5-diameter sunken pan. The quantity of water added to maintain the
water level constant over a period of 10 h was 950 cc. Find the rate of evaporation.
Ans: 0.5 mm/h/m2
6. The following observations were taken on a 1.0-m diameter land pan from 8 a.m. to 5 p.m.
(i) Precipitation during the observation period 15 mm
(ii) Leakage from the pan 1.5 l
(iii) Quantity of water added to keep the level constant 4.0 l
Find the rate of evaporation.
Ans: 2.02 mm/h/m2
7. The following observations were taken from a reservoir:
(i) Velocity of wind at 0.5 m above ground 24 m/s
(ii) Atmospheric pressure 10.4 m of water
(iii) Average air temperature 30 °C
(iv) Average water temperature 27 °C
(v) Relative humidity 0.5
(vi) Water spread of reservoir 15 km2
(vii) Saturation vapour pressure of air 30 mm Hg
Find the monthly evaporation using empirical formulae.
Ans: (i) 22.50 106m 3
(ii) 11.25 106m 3
(iii) 4.87 106m 3
(iv) 10.36 106m 3
8. Transpiration is confined to
(a) Dayt imehour s (b) Nighthour s
(c) All the day (d) Any of the three
10. Out of the water sucked by the roots, the percentage retained by plants is
(a) 1 (b) 5
(c) 10 (d) 15
15. Evaporationi sa
(a) Coolingpr ocess (b) Heating process
(c) Combined process (d) None of the three
50 ELEMENTARY ENGINEERING HYDROLOGY
*** ***
Precipitation
4
Chapter Outline
g 4.1 DEFINITION
All water flowing or stored on the land surface or subsurface is derived directly or indirectly from
precipitation, i.e. rainfall including snowfall. It is one of the important basic processes in the hydro-
logical cycle.
Precipitation is the general term for all the moisture emanating from the clouds and falling
on the ground. Precipitation data are of utmost importance to hydrologists as they form the basis
of all hydrological studies. Variation of rainfall distribution over time and space creates extreme
problems like droughts and floods. Study of precipitation, therefore, requires great attention.
52 ELEMENTARY ENGINEERING HYDROLOGY
Four conditions are necessary for the formation of precipitation. These are:
1. Mechanism to cool the air
2. Mechanism to produce condensation
3. Mechanism to produce growth of droplets
4. Mechanism to produce accumulation of moisture of sufficient quantity
2. Ice crystal theory In the high atmosphere where condensation takes place, the temperature
is 0 oC or even lower. The droplets thus formed are in the form of ice crystals or even water
droplets. There is a difference in saturation pressure over the ice crystals and the water droplets.
Thus evaporation of water droplets and condensation over the ice crystals occurs.
This process is comparatively slow and occurs in very cool (super cooled) clouds, i.e. having
temperatures lower than 0 oC.
4.4.6 SUNSPOTS
There are spots on the Sun and scientists have observed that these sunspots have a cycle of variation
of 11 years, i.e. the nature of the sunspots are repeated after every 11 years. Some scientists believe
that there is a correlation between the hydrological cycle and the sunspots cycle, but have failed to
prove any exact correlation.
The movement of such whirling air results in cyclone formation. Some features of a cyclone are
mentioned below:
• A cyclone may be defined as a whirling mass of air at the centre of which the barometric pres-
sure is low.
• It is also known as hurricane, blizzard, typhoon and tornado according to the whirling speeds. An
anticline similarly whirls round a high-pressure centre.
• A cyclone comprises of a very large mass of air spread over an area from 10 km2 to as high as
1500 km2 moving with a velocity of 100 km/h or even more.
• The whirling air mass carries water vapour with it. The central portion of the cyclone acts as a
chimney through which the air gets lifted, expanded and cooled, and the water vapour gets con-
densedc ausingpr ecipitation.
Thus, precipitation caused due to a cyclone is known as cyclonic precipitation and may result
into a drizzle or into heavy precipitation covering a large area.
The instrument used to measure rainfall is known as rain gauge. All forms of precipitation
are measured on the basis of the vertical depth of water that would accumulate on a plain
surface, if the precipitation remains where it falls. It is measured in millimetres or tenths of
millimetres.
In olden days, a rain gauge was known as hyetometer, ombrometer or pluviometer.
the basic data on which all the designs are based and hence utmost accuracy has to be observed
during its measurement. The precipitation per unit time is the intensity of precipitation.
It is measured as mm/h or cm/day.
4.6.3 HYETOGRAPH
A bar chart of time versus precipitation is known as hyetograph. The ordinate graph presents the
rainfall in a year drawn to some scale at the corresponding year.
Specific and ordinate graphs are shown in Figs. 4.2 and 4.3, respectively.
200
160
120
80
40
0 1 2 3 4 5 6 7 8 9 10 11 12
Example 4.1
The hourly precipitation data during a storm are as follows:
Ti e (h)
m 0 1 2 3 4 5 6 7 8 9 10
Precipitation (mm) 0 30 25 50 5 0 10 15 25 20 0
Solution:
The mass curve coordinates will be as follows:
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12
Precipitation (mm) 0 30 55 105 110 110 120 135 160 180 180 180 180
The hyetograph and the mass curve are shown in Fig. 4.2.
The rainfall observed at a station by a rain gauge is known as point rainfall or station rainfall. The
rainfall observed by the non-recording or recording rain gauges is point rainfall and hence these
gauges are known as point precipitation gauges. The third one, i.e. by the radar, gives an instanta-
neous picture of the rainfall over an area and not the point precipitation.
Sharp edged
rim
125 mm
Funnel
Receiver
300 mm
Body
100 mm 50 mm
600 mm
600 mm
10
8
6
4
2
0
8 10 12 14 16 18 20 22 24 2 4 6 8
Time in hours
The vertical line indicates that the siphon has worked and a horizontal line indicates no precipi-
tation. The slope of the chart indicates the intensity of rainfall.
The siphoning action requires 15 seconds and record of rainfall during this period is lost. So also the
rainfall record may be lost during the time taken in changing the paper, if there is rainfall during that time.
Example 4.2
The chart fixed to an automatic float-type rain gauge is shown in Fig 4.5
Find 1. Hourlypr ecipitation
2. Dailypr ecipitation
3. Time when the siphon was operated
4. Period of no precipitation
5. Maximum intensity of precipitation
Solution:
1. The hourly precipitation as read from the chart will be as follows:
Time 8 9 10 11 12 13 14 15 16 17 18 19
Cumulative ppt (mm) 0 0.25 0.5 1 2.5 4 4 4 4.9 6 7.4 8.5
Precipitation (mm) 0 0.25 0.25 0.5 1.5 1.5 0 0 0.9 1.1 1.4 1.1
Time 20 21 22 23 24 1 2 3 4 5 6 7 8
Cumulative ppt (mm) 10 0.5 1 1.8 2.8 4.2 6 6 6 6.5 7 7.5 8
Precipitation (mm) 1.5 0.5 0.5 0.8 1 1.4 1.8 0 0 0.5 0.5 0.5 0.5
30 cm Receiver
Funnel
Stand
Measuring jar
bucket. However, for low intensity of rainfall, the results may not be very accurate since no-rainfall
period is not recorded in this rain gauge.
Receiver
Funnel
Case
Bucket
Platform
Spring
fixed on a drum. The drum is driven mechanically by a spring clock and makes one revolution a
day at a uniform speed. Thus, the trace of the pen on the paper records the weight of rainwater and
indirectly the rainfall.
After a certain amount of precipitation, the mechanism reverses the travel of the pen. However,
the movement of the drum remains unaltered. Figure 4.8 shows a typical rainfall chart of a weigh-
ing bucket-type rain gauge.
The horizontal line indicates the no-rainfall period. Reversal of line has to be taken into
account. With this, the total precipitation as well as the intensity of precipitation can be calcu-
lated accurately.
This type of rain gauge can be used for the measurement of snow, rainfall, hail, and so on.
Example 4.3
The chart fixed to an automatic weighing bucket-type rain gauge reads as in Fig. 4.8:
Find 1. Hourlypr ecipitation
2. Dailypr ecipitation
3. Time when the pointer reverted
4. Period of no precipitation
5. Maximum intensity of precipitation
Solution:
1. The hourly precipitation as read from the chart will be as follows:
Time 8 9 10 11 12 13 14 15 16 17 18 19
Cumulative ppt (mm) 0 0 0 5 18 30 37 40 44 48 56 66
Precipitation (mm) 0 0 0 5 13 12 7 3 4 4 8 10
Time 20 21 22 23 24 1 2 3 4 5 6 7 8
Cumulative ppt (mm) 70 72 78 83 84 87 89 96 100 70 65 55 52
Precipitation (mm) 4 2 6 5 1 3 2 7 4 30 5 10 3
PRECIPITATION 65
Sometimes there might be breaks or gaps in the rainfall data of some stations. These gaps may be
due to the following reasons:
• Absence of an observer
• Instrumentsf ailure
• Unapproachable circumstances such as flooding or washing out of the approach road, etc., result-
ing in non-recording of the rainfall
The missing rainfall data may be supplemented by the following methods:
• Arithmetic average method
• Normal ratio method
• Weighted average method
Example 4.4
In a catchment area, daily precipitation was observed by 11 rain gauge stations. On 2 August 2005,
the observations indicated that one rain gauge was out of order. The observations taken by the 10
rain gauge stations are as follows:
Station A B C D E F G H I J K
Precipitation (mm) 21 23 19 20 23 24 19 ? 21 22 18
Solution:
Since there is not much variation in the precipitation data, a simple arithmetic average of the
precipitation observed at the 10 remaining stations was taken as under.
21 23 19 20 23 24 19 21 22 18
_____________________________________________
Precipitation at H 10
____
210
10 21 mm
[
NA
__ ___
B
NA
___
PA (missing) n1 N PB N PC N Pn
C
NA
___
n
]
where, NA Average of A excluding the missing period
NB NC Average of B and C excluding the missing period
PB PC Precipitation at B and C during the missing period
n number of adjacent rain gauge stations considered
For this method, a minimum of three adjacent rain gauge stations are considered.
Example 4.5
The average annual precipitation at five rain gauge stations in a catchment is as follows:
Station P Q R S T
Average annual precipitation (mm) 2400 2332 2431 2207 2231
However, the precipitation at station P was not available for the year 1996 because the rain
gauge was out of order. The precipitation observed at the other stations for 1996 was as follows:
Station P Q R S T
Precipitation (mm) ? 2113 2200 2028 2095
Solution:
__
1 _____ _____ _____ _____
(
2400 2400 2400 2400
Precipitation at P (in mm) 4 2332 2113 2431 2200 2207 2028 2231 2095 )
__
1
4 (2174 2172 2205 2253)
2201 mm
where, r1, r2 and r3 the distances from A of the adjacent stations B, C and D
PB, PCa nd PD the rainfall at the adjacent stations B, C and D for the missing period.
This method is also known as United States National Weather Service (USNWS) method.
Example 4.6
The location coordinates in km of the five rain gauge stations w.r.t. X are as follows:
Station X A B C D
x and y coordinates (km) w.r.t. X 0,0 20,25 40,15 30,20 25,15
The annual precipitation at X for the year 2005 is missing. The annual precipitation at the remain-
ing four stations for 2005 is as follows:
Station A B C D
Precipitation (mm) 2735 2805 2680 2560
Solution:
The relative positions of the stations A, B, C and D w.r.t. X are shown in Fig. 4.9.
AX √(202 25 2) 32.01 km
BX √(402 15 2) 42.72 km
CX √(302 20 2) 36.05 km
DX √(252 15 2) 29.15 km
PRECIPITATION 69
20 km A
B 40 km
25 km
15 km
15 km
20 km
25 km D
C 30 km
It may happen that the rainfall recorded by a rain gauge station is doubtful. It then becomes neces-
sary to verify the rainfall record of this station. This is known as verification of consistency of a rain
gauge. This may by due to the following reasons:
• Change in the location of the rain gauge
• Change in the surroundings, namely, growth of trees, buildings, and so on.
• Change in the instrument
• Fault developed in the rain gauge
The verification can be done by the double mass curve method.
4500
1941
4000
B B´
1945
Cumulative annual rainfall at station x
3500
3480
α´
3000
1950
α
1953
2500
1955
C
2000
1960
1500
1965
tan α´ = 1.17
tan α = 0.85
1000
1970
500
1975
0 A
500 1000 1500 2000 2500 3000 3500 4000 4500
Cumulative average annual rainfall at base stations
Normally, one should get a straight line through origin if the record at A is correct. If there
is inconsistency at A from a particular year, the slope of the straight line may change from that
year. It may, therefore, be concluded that the records of A are incorrect from that year and need
modification.
The slope of the straight line is maintained and extended and the record of A is corrected accord-
ingly. This procedure cannot be applied for studies for storm rainfall or daily rainfall.
A double mass curve plotted for a specific study is shown in Fig. 4.10.
Example 4.7
The average annual precipitation data of six rain gauge stations in a catchment area during
1991–2000 are as follows:
Station A B C D E F
Year Precipitation (mm)
The data observed at station C were doubtful because of some topographical changes there. Verify
whether the data observed at C is consistent and correct it, if necessary.
Solution:
The mass curve coordinates of precipitation at C and also the combined mass curve coordinates of
precipitation at A, B, D, E and F will be as follows:
Mass curve Corrected mass for Corrected Corrected C
Serial no. for C A, B, D, E and F curve for C precipitation at C (mm) Remarks
1 378 3200 — —
2 780 6689 — —
3 1168 9796 — —
4 1601 13501 — —
5 1975 16693 — —
6 2421 20495 — —
7 2823 23895 — — ←Graph
8 3540 28092 3302 479 changes
9 4240 31890 3756 454 slope
10 4920 35496 4173 417
The precipitation over a catchment area is never uniform. This becomes quite clear from the figures of
the average depth of precipitation of the various rain gauge stations in the catchment area. One of the
72 ELEMENTARY ENGINEERING HYDROLOGY
basic requirements in the study of a catchment area is the average depth of precipitation over the
entire catchment.
This is also known as equivalent uniform depth of rainfall. The average depth of precipitation
can be calculated by the following methods:
• Arithmetic mean method
• Thiessen polygon method
• Isohyetalm ethod
PA PB PC …
_________________
Pav n
where, Pav Average precipitation over the catchment area
PA, PB, PC Average precipitation at various rain gauge stations A, B, C, …
This method may be adopted when the catchment area is flat and there is not much variation
in the average values of the precipitation of all the rain gauge stations. This method is very simple
for calculations. However, no weightage is given neither to the influence area of the individual rain
gauge nor the topography. This method is also known as unweighted mean method.
PA AA PB AB PC AC …
_________________________
Pav
A
where, Pav Average precipitation over the catchment area
PA, PB, PC Average precipitation at various rain gauge stations
PRECIPITATION 73
AA, AB, AC Area of the polygon surrounding the rain gauge stations A, B and C
A Catchment area AA AB AC…
Pav can also be mentioned as follows:
AA AB
___ AC
___
Pav ___ PA A PB A PC …
A
AA ___
___ AB AC
___
and the factors A , A a nd A are called Thiessen’s weights.
In this method, weightage to each rain gauge is given depending upon its influence area. Rain
gauges having hydrologic homogeneity outside the area, if any, are also considered for more accu-
racy. However, no consideration is given to the topography.
Example 4.8
Figure 4.12 shows a typical layout of a catchment area ABCDF. Six rain gauge stations are estab-
lished at A, B, C, D, E and F, as shown in the figure. The precipitation observed at these six stations
in July 2004 is as follows:
Station A B C D E F
Precipitation (mm) 100 120 130 180 125 150
Find the average precipitation over the catchment during July 2004 by Thiessen’s polygon method.
A 4 km B
150° 90°
m
4k
4 km
F 60°
E
4k
m
150° 90°
D 4 km C
A P B
S •E Q
F W
D R C
Solution:
Join the locations of the adjacent rain gauge stations, AB, BC, CD, DA, AE, BE, CE, DE, DF
and FA, and draw their perpendicular bisectors. The perpendicular bisectors will meet at P, Q,
R, S, T, V and W.
Thus, each rain gauge station will be surrounded by a polygon and its area will be as follows:
1 A APSWV 1 2 2 __
__ 1 4 √__
1 __ 3 3.15
2 3 2
2 B BQP 1 2 2 2.00
__________
2
3 C CRQ 1 2 2 2.00
__________
2
4 D DTWSR 1 2 2 __
__ 1 4√__
1 __ 3 3.15
2 3 2
5 E PQRS 4 4 8.00
______
2
__
6 F FVWT 1 _________
__ 1 4 √3 1.15
3 2
Total 19.45 km2
5. The average precipitation of the area between two isohyets may be calculated as follows:
i(2a b)
________ i(a 2b)
________
P1 B or A
3(a b) 3(a b)
Example 4.9
The daily precipitation data observed at four rain gauge stations located inside a catchment area on
2 August 2005 are as follows:
Station B D F H
Precipitation (mm) 30 35 50 45
So also the daily precipitation data observed at four other rain gauge stations that are meteorologi-
cally similar but outside the catchment area on the same day are as follows:
Station A C E G
Precipitation (mm) 40 25 45 55
PRECIPITATION 77
Fig. 4.14 Rain gauge surrounded by a polygon inside the catchment area
Figure 4.14 shows the catchment area and the locations of the rain gauge stations. If the catchment
area is 54 km2, find the daily average precipitation over the catchment.
Solution:
The average daily precipitation over the catchment will be calculated by three different
methods as stated below.
1. Arithmetical mean method
2. Thiessen’s polygon method
3. Isohyetalm ethod
1. Arithmetical mean method The average daily precipitation of the catchment area will be
the average of the four rain gauge stations located inside the catchment area.
30 35 50 45
_________________
Average daily precipitation of the catchment area 4 mm
____
160
4 40 mm
2. Thiessen’s polygon method Join the locations of the adjacent rain gauge stations by
straight lines as, AB, AH, BH, BC, BD, CD, DE, DF, EF, FG, GH and HF.
Draw perpendicular bisectors of these lines to meet
• Inside the catchment area at M, N, P, Q, R, S
• The catchment area boundary at a, b, c, d, e, f, g, h, i, j, k, l, m, n, y
78 ELEMENTARY ENGINEERING HYDROLOGY
Each rain gauge will be surrounded by a polygon inside the catchment area and its area will be as
follows:
Therefore, the average daily precipitation over the catchment area on 2 August 2005 will be
(3.540) (11.630) (1.525) (11.935) (3.245) (9.750) (2.055) (10.645)
__________________________________________________________________________________
3.5 11.6 1.5 11.9 3.2 9.7 2.0 10.6
_____
2158
54 39.96 mm
3. Isohyetal method Taking into consideration the location of the rain gauge stations inside
(B, D, F, H) as well as outside the catchment area (A, C, E, G) and the precipitation recorded
by them, isohyets of 30, 35, 40, 45, 50 mm are drawn and area between them was measured
accurately (Fig. 4.15).
Therefore, the average daily precipitation over the catchment area on 2 August 2005
3010 32.511 37.511.3 42.59.3 47.57.6 504.8
_________________________________________________________
54
300.0 357.5 423.75 395.25 361.0 240
__________________________________________
54
38.47 mm
Rainfall is the most fundamental data used in the hydrological studies and hence a well distributed
network of rain gauge stations is essential. The average area of influence of the rain gauge stations
is indicated as rain gauge density or network density. The density of rain gauge stations in an area
may be decided taking into consideration the following points:
• Variation in the rainfall in the area. If the area is plain and if there is not much variation in rainfall,
the number of rain gauge stations may be small.
• The nature of study for which rainfall data is required
• Cost involved in establishing and maintaining the rain gauge stations
(Pi Pav)2
_________
Standard deviation sx 2 n–1 (i 1, 2, 3, … n)
sx 100
Cv Coeficient of variance _________
Pav
N Number of optimum rain gauge stations ( Cv x)2
x Percentage permissible error in the estimation of average rainfall.
Additional rain gauge stations N – n
It may be noted that both Cv and x be mentioned as percentages.
The additional rain gauge stations may be located in addition to the existing rain gauge stations,
so that they all are evenly distributed over the entire catchment area.
Example 4.10
In a catchment area covering 100 km2, the average annual precipitation observed at five rain gauge
stations is as under.
Station 1 2 3 4 5
Precipitation (mm) 750 1000 900 650 500
Find the number of additional rain gauge stations and also the rain gauge density if the permissible
error is 10%.
Solution:
750 1000 900 650 500
___________________________
Average precipitation 5
760 mm
(750 760)2 (1000 760)2 (900 760)2 (650 760)2 (500 760)2
__________________________________________________________________
sx2 51
_______
157000
4 39250
______
Therefore, sx √ 39250 198.1
_____
198.1
Cv 760 100 26.06%
Since permissible error is 10%,
Number of rain gauge stations required ( Cv/10)2 ( 26.06/10)2
6.79
Additional rain gauge stations required 6.79 5 1.79 艐2
Rain gauge density 100/7 14.30km 2/rainga uge.
Example 4.11
In a catchment area covering 5600 km2, the zone-wise existing rain gauge stations were as follows:
Zone A B C D E Total
Additional nine rain gauge stations are to be installed. Indicate the zone-wise distribution of these
nine additional rain gauge stations.
Solution:
The total number rain gauge stations in the catchment N will be 7 9 16.
The additional nine rain gauge stations will be distributed proportional to the area of the zones
and also taking into consideration the number of existing rain gauge stations, as follows.
Zone A B C D E Total
P Area in decimal of total area 0.143 0.089 0.286 0.447 0.035 1.000
N P rounded 2 1 5 7 1 16
Number of existing rain gauge stations 1 1 2 3 0 7
Additional rain gauge 1 0 3 4 1 9
Example 4.12
The average annual precipitation in millimetres observed over a catchment area from 1980 to 1995
is as follows:
1149, 1260, 1425, 1680, 1200, 1400, 1645, 1500,
1155, 945, 875, 1080, 1480, 1625, 1500, 1470
Construct a 3-year moving average curve and plot it along with the original data.
82 ELEMENTARY ENGINEERING HYDROLOGY
For an area, the maximum depth of precipitation that may occur for a specific duration is known
as possible maximum precipitation or probable maximum precipitation. The usual short form used
is PMP. This data is required for estimating maximum possible flood from a catchment area. The
PMP for a known duration can be correlated as follows:
PMP Pa K sx
where, Pa Average precipitation over the area
sx Standard deviation of rainfall series
K A constant for the area and is in the neighbourhood of 15.
The data of world’s greatest rainfall was analysed and it was observed that it follows a curve
PMP 42 .16 D0 475
where, PMP Precipitation in centimetres
D Duration in hours
Probability
Probability generally denoted by p is reciprocal of the return period, i.e. p m/n 1/ Tr.
Frequency
Probability expressed in terms of percentage is frequency, i.e. frequency p 100 m/n 100
84 ELEMENTARY ENGINEERING HYDROLOGY
Example 4.13
The precipitation in millimetres observed at a rain gauge station for the last 32 years is as follows:
988, 966, 935, 1007, 992, 1050, 975, 920, 1035, 990, 1095, 1015, 986, 927, 1003, 1055, 1135,
955, 1001, 1045, 1090, 997, 1040, 1100, 948, 972, 1012, 950, 1070, 982, 929, 960
Find 1. The return period and frequency of the precipitation of 997 mm
2. The precipitation of return period of 1.33 and its frequency
Solution:
The precipitation figures in millimetre were arranged in a descending order as follows:
Serial no. 1 2 3 4 5 6 7 8 9 10
1135 1100 1095 1090 1070 1055 1050 1045 1040 1035
Serial no. 11 12 13 14 15 16 17 18 19 20 21
1015 1012 1007 1003 1001 997 992 990 988 986 982
Serial no. 22 23 24 25 26 27 28 29 30 31 32
975 972 966 960 955 950 948 935 929 927 920
It is observed that most intense storms last for a short time. As the intensity reduces, the duration
of the storm increases and vice versa. The study of intensity and its duration is known as intensity
duration analysis.
_______
C
I
(t a)b
Example 4.14
A storm occurred over a catchment area as under:
Time (min) 0 10 20 30 40 50 60 70 80 90
Precipitation (mm) 0 19 22 7 20 23 33 28 8 6
Solution:
Time (min) Cumulative precipitation in mm at time interval
10 20 30 40 50 60 70 80 90
0 0
10 19
20 22 41
30 7 29 48
40 20 27 49 68
50 23 43 50 72 91
60 33 56 76 83 105 124
70 28 61 84 104 111 133 152
80 8 36 69 92 112 119 141 160
90 6 14 42 75 98 118 125 147 166
Maximum intensity
198 183 168 156 134 133 130 120 110
(mm/h)
Rainfall over India is very erratic in terms of both time and space. The coefficient of variation of
the annual rainfall varies between 15 and 36, and hence it is said that India’s prosperity is a gamble
in the monsoon’s rains.
The average annual precipitation over India is 1140 mm. Both the maximum and the minimum
precipitations in the world are observed in India. The important features affecting precipitation over
India are the following:
1. The orographic features
2. The wind currents
The seasons in the Indian subcontinent can be divided into four major seasons.
1. South-west monsoon (June to September)
2. Post-monsoon (October to December)
3. Winter season (January to February)
4. Pre-Monsoon (March to May)
The chief characteristics of these seasons are as follows.
4.13.2 POST-MONSOON
As the southwest monsoon retreats, low-pressure areas develop in the Bay of Bengal and a north-
easterly flow of air that picks up moisture in the Bay of Bengal is formed. This air mass strikes
the east coast of the southern peninsula, particularly Tamil Nadu, and causes precipitation in those
areas. In addition, during this period several tropical cyclones form in the Bay of Bengal and to a
lesser extent in the Arabian Sea. These strike the coastal areas and cause intense precipitation.
4.13.4 PRE-MONSOON
There is hardly any precipitation during this season. Convective cells cause some thunderstorms,
mainly in Kerala, West Bengal and Assam. Some cyclone activity, primarily in the eastern coast,
also occurs. Figure 4.23 shows the average annual precipitation over India.
REVIEW QUESTIONS
1. Define precipitation. Explain its importance in the study of hydrology.
2. What are the different forms of precipitation?
3. Discuss the factors affecting precipitation.
4. Explain the different types of precipitation.
5. Explain artificial precipitation. Why is it not followed on a large scale in India?
6. Discuss the considerations for selecting the site for a rain gauge.
7. Explain the different types of rain gauges. Discuss their merits and demerits.
PRECIPITATION 91
8. Explain with the help of a neat sketch, the working of a non-recording rain gauge. Discuss its
merits and demerits.
9. Explain with the help of a neat sketch, the working of a siphon bucket-type rain gauge.
10. Explain with the help of a neat sketch, the working of a tipping bucket-type rain gauge.
11. Explain with a neat sketch, the working of a weighing bucket-type rain gauge.
12. Discuss the difficulties experienced during the measurement of precipitation.
13. Why supplementing of rainfall data is required? State the different methods in use.
14. Explain verification of consistency of a rain gauge.
15. Explain the double mass curve method.
16. How the average precipitation over a catchment is calculated? Discuss the different methods
with their merits and demerits.
17. Explain the ISI standards for rain gauge density in a catchment. Also, discuss the procedure
to calculate the optimum number of rain gauge stations in a catchment.
18. What is a moving average curve? What are its advantages? Explain the procedure to draw it.
19. Explain the depth–area–duration analysis.
20. What is an intensity–duration curve?
21. Write a detailed note on precipitation over India.
22. Write a note on frequency of a storm.
23. Write short notes on the following:
a . Terminalve locity b. Orographicpr ecipitation
c. Cloud seeding d. History of measurement of precipitation
e . Use of radar for measurement of precipitation f. Index of wetness
g. Classification of the intensity of rainfall h. Rain gauge density
i. Probable maximum precipitation j. Station year method
k. Recurrence interval of a storm l. Isopluvial map
m. Transposition of a storm
24. Differentiatebe tween
a. Record graph sheets of siphon bucket-type b. Forms of precipitation and types
rain gauge and a weighing bucket-type of precipitation.
rain gauge.
c. Cyclonic precipitation and orographic d. Non-recording-type rain gauge and
precipitation. recording-type rain gauge.
e. Siphon-typer ainga ugea ndw eighing-type
rain gauge.
NUMERICAL QUESTIONS
1. During a storm, the hourly precipitation data observed was as follows:
Time (h) 0 1 2 3 4 5 6 7 8 9
Precipitation (mm) 0 25 20 45 30 10 0 15 10 0
2. The precipitation coordinates in mm observed on an automatic float-type rain gauge are as follows:
Time (h) 8 9 10 11 12 13 14 15 16 17 18 19 20
Precipitation (mm) 0 3 11 21 35 40 40 50 65 75 84 88 92
Time (h) 21 22 23 24 1 2 3 4 5 6 7
Precipitation (mm) 100 4 8 11 21 21 21 28 32 40 42
Time (h) 21 22 23 24 1 2 3 4 5 6 7
Precipitation (mm) 100 97 85 80 75 75 75 72 68 65 62
The precipitation at B was missing for the year 1998 since the rain gauge was out of order.
The precipitation observed at other stations for 1998 are as follows:
Station A B C D E
Precipitation (mm) 2465 ? 2492 2285 2511
The annual average precipitation at the Station P was missing for the year 2004.
The average annual precipitation at other stations for the year 2004 is as follows:
Station P Q R S T
Precipitation (mm) ? 2735 2800 2831 2655
Ans: The precipitation data at Station A is wrong from the year 1996 and the corrected
precipitation figures from 1990 in mm are 1388, 1412, 1398, 1443, 1384, 1456,
1412, 1489, 1464, 1427, respectively.
7. A catchment area is in the form of an equilateral triangle ABC of side 10 km. Four rain gauge
stations are located at A, B, C and D. Station D is the centroid of the triangle. The average annual
precipitation observed at these stations in mm are 1145, 1252, 1184 and 1056, respectively. Find
the average annual precipitation of the catchment area by all the three methods.
Ans: 1101 mm
8. The average annual precipitations of the five rain gauge stations in a catchment are as follows:
Station A B C D E
Precipitation (mm) 1100 890 1000 1370 1805
94 ELEMENTARY ENGINEERING HYDROLOGY
If the catchment covers an area of 110 km2, find the number of additional rain gauge
stations required. Acceptable permissible error is 10%.
Ans: 4
9. The average annual precipitation in mm observed over a catchment area for the last 16 years is
as follows:
1154, 1265, 1430, 1683, 1205, 1390, 1650, 1497, 1160, 940, 870, 1080, 1490, 1630, 1505, 1475
Construct a 3-year moving average curve and compare with the original data.
10. The average annual precipitation in mm observed at a rain gauge station for the last 32 years
is as follows:
966, 992, 920, 1045, 1012, 1070, 960, 935, 1050, 1035, 1015, 948, 972, 950, 929, 1007, 975,
990, 1095, 986, 988, 927, 1003, 1055, 982, 1135, 955, 1001, 1040, 1090, 997, 1100
Find: (i) The return period and the frequency of the precipitation of 1045 mm
(ii) The precipitation of the return period of 2 and its frequency
Ans: (i) 4,25%
(ii) 997, 50%
11. The storm occurred over a catchment area is a under:
Time (s) 0 60 120 180 240 300 360 420 480 540 100
Precipitation (mm) 0 20 23 8 21 24 35 27 9 5 0
5
Chapter Outline
g 5.1 DEFINITION
Infiltration may be defined as entry and movement of water through the land surface into the sub-
strata below.
5.1.2 ABSTRACTIONS
When there is precipitation, it may or may not result in overland flow into a stream depending upon
its intensity and duration. The part of precipitation that is not available as surface runoff is referred
INFILTRATION 99
to as precipitation loss or abstraction. The abstractions include (1) interception, (2) depression
storage, (3) evaporation and (4) infiltration.
Generally, interception and depression storages are together termed as ‘surface retention’. Knowledge
of these abstractions, their rate, and so on, is necessary for the determination of the surface flow.
When rainwater falls on the ground, there is some resistance offered by the soil surface for the
entry of rainwater and also to the flow of water through the soil. There are cracks, vertical as well
as lateral in the soil, so also there are some voids between the soil particles, which are ordinarily
occupied by air or water. Water flows through these cracks and gaps until it reaches the saturated
zone below. Naturally, the passage of water experiences some resistance.
The rate at which water enters the ground surface and then flows downwards is known as
infiltration rate. This rate is high in the beginning because it has to meet the requirements of the dry
soil. However, it attains a steady constant lower value after a passage of time.
The unit for the rate of infiltration is mm/h.
• Fine soil particles are carried down due to rain water resulting in choking of the pore spaces in
the soil, and consequently resulting in reducing the infiltration rate.
Uniform precipitation
Moist soil
Time
The rate of infiltration is initially high. It goes on reducing with time and after some time it becomes
steady. A usual graph of the rate of infiltration is shown in Fig. 5.1. The rate of infiltration for a soil
is measured in the field as well as in the laboratory. These are known as infiltrometers.
The most common types are the following:
• Flooding-type infiltrometers
• Sprinkling-type infiltrometers or rain simulators
30 cm
10 cm G. L .
50 cm
60 cm
30 cm
10 cm
G. L .
15 cm
Example 5.1
The quantity of water added to a double-ring infiltrometer of 1.00-m diameter at 30 min interval to
keep the water level constant is as follows:
Find: 1. Rate of infiltration for every 30 min and plot the graph
2. Average rate of infiltration
104 ELEMENTARY ENGINEERING HYDROLOGY
27
26
22
21
20
19
0 1 2 3
Time in hours
Solution:
12 0.786 m2
Area of ring __
4
10 103
________ ___
60
Infiltration rate for 1st 30 min 0.786 103 30 25.44 mm/h
9.2103 ___
60
Infiltration rate for next 30 min _________
0.786 10 30 23.40 mm/h
3
3
________
8.610 ___
60
Infiltration rate for next 30 min 0.786 10 30 21.88 mm/h
3
3
________
8.210 ___
60
Infiltration rate for next 30 min 0.786 10 30 20.86 mm/h
3
8103
_______ ___
60
Infiltration rate for next 30 min 0.786 103 30 20.35 mm/h
The graph of time in hours versus infiltration rate is as shown in Fig. 5.4.
Total quantity of water added in 150 min till a steady state was achieved 10.0 9.2 8.6 8.2
8.0 44.0 lit
44 103 103 60
___________________
Therefore, average rate of infiltration 150 0.786 22.39 mm/h/m2
Catchment area
Ridge line
Discharge
measurement site
The surface runoff is measured very accurately at the catchment outlet and thus knowing all the
abstractions such as ‘evaporation, interception, depression storage, surface detention’ and surface
runoff, the infiltration rate during the time of experiment can be worked out.
Consider the following particulars of a specific experiment.
A small catchment was selected for the experiment as shown in Fig. 5.5.
The experiment started at 8:00 a.m. by sprinkling water over the entire catchment area at a uni-
form rate, which was stopped at 12:00 noon. There was no surface runoff noticed up to 8:30 a.m.
The discharge was noticed at 8:30 a.m. and it went on increasing up to 10:00 a.m. and remained
constant up to 12:00 noon. When the sprinkling of water, i.e. rainfall, was stopped, the discharge
went on reducing and finally stopped at 2:00 p.m. as shown in Fig. 5.6.
Recession curve
8 8-30 9 10 11 12 1 2
Time in hours
8 9 10 11 12 1 2
Time in hours
mm/hour
Mass curve of ‘Surface retention’
8 8.30 9 10 11 12
Time in hours
in this abstraction after 8:30 a.m. Surface runoff started only after the requirement of the surface
retention was met with. The surface retention and its mass curve is shown in Fig. 5.8.
The discharge was observed from 8:30 a.m. onwards. It went on increasing up to 10:00 a.m. It
was constant from 10:00 a.m. to 12:00 noon and then reduced from 12:00 noon to 2:00 p.m. The
discharge as well as its mass curve is shown in Fig. 5.9.
The SD is in mm/s, i.e. water stored in the channel is a function of depth in a channel. More
the discharge in the channel, more will be the depth in the channel and more will be the storage in
the channel. Thus, indirectly SDwill be a function of discharge. Here the discharge was started at
8:30 a.m. Similarly, SD also started at 8:30 a.m. It went on increasing up to 10:00 a.m. It remained
constant from 10:00 a.m. to 12:00 noon and then was reduced from 12:00 noon to 2:00 p.m.
Discharge in mm/hour
Discharge
8 8.30 9 10 11 12 1 2
Time in hours
Q3
Recession curve
Discharge in mm/hour
Q2
Q1
8 8.30 9 10 11 12 1 2
Time in hours
A3 A2 A1
Since initially and also after the experiment the channel carried no water, the water stored in the
channel from 8:30 a.m. to 12:00 a.m. was drained completely from 12:00 noon to 2:00 p.m.
When the precipitation stopped at 12:00 a.m., there was still flow in the stream. This was due
to draining the water stored in the channel, i.e. SD. Immediately after the precipitation stopped, the
discharge was maximum. SD was also maximum. The total volume of water stored in the channel
was the area under the recession curve. As the discharge reduced, so also did the storage in the
channel. The storage in the channel will be the area under the recession curve. The recession curve
is shown in Fig. 5.10.
The area under the recession curve is the total channel storage.
A A1 A2 A3
When the discharge is Q1, the channel storage SD is A1.
When the discharge is Q2, the channel storage SD is A1 A2.
When the discharge is Q3, the channel storage SD is A1 A2 A3, which is the area under the
recession curve.
Thus, the relation between Q and SD can be established as shown in Fig. 5.11.
SD can be calculated from this graph and then correlated with time. SD thus calculated will be
in terms of m3 and will have to be worked out in terms of mm over the entire catchment area for
further calculations. The surface detention and its mass curve are shown in Fig. 5.12.
A combined mass curve of (SR Q SD) and mass curve of P is plotted on the graph paper and
the difference between these two mass curves is the infiltration rate w.r.t. time as shown in Fig. 5.13.
From this figure, the infiltration rate w.r.t. time can be plotted. The infiltration rate is more at
the beginning and goes on reducing, and after some time it becomes steady as shown in Fig. 5.14.
INFILTRATION 109
Surface detention in mm
Discharge in mm/hour
A1 = A2
A1
8 8.30 9 10 11 12 2
Time in hours
A2
mm/h
Infiltration rate
Combined mass
curve of Sr+Q+Sd
8 8.30 9 10 11 12 1 2
Time in hours
Example 5.2
A field experiment to assess the infiltration capacity of an area was conducted by following the
Sprinkling type infiltrometer technique’. The discharge was measured accurately over a notch at
10 minutes interval.
1 2 3 4
Sr. no. Time Discharge in c. c. /s Discharge in cm / h
1 8–00 0.00 0.00
2 –10 0.00 0.00
3 –20 1.88 0.15
4 –30 5.00 0.40
5 –40 7.38 0.59
6 –50 -10.00 0.80
7 9–00 12.88 1.03
8 –10 16.00 1.28
9 –20 19.25 1.54
10 –30 22.50 1.80
11 40 25.00 2.00
12 –50 27.25 2.18
13 10–00 28.75 2.30
14 –10 29.75 2.38
15 –20 30.75 2.46
16 –30 31.25 2.50
17 –40 31.75 2.50
18 –50 31.75 2.50
19 11–00 17.00 1.36
20 –10 12.25 0.98
21 –20 9.50 0.76
22 –30 7.88 0.63
23 –40 6.63 0.53
24 –50 5.75 0.46
25 12–00 4.63 0.37
26 –10 4.25 0.34
27 –20 3.88 0.31
28 –30 2.50 0.20
29 –40 1.88 0.15
30 –50 0.00 0.00
The discharge measured on the notch in c. c. /s, was converted for simplicity, into depth in cm /h
uniformly distributed over the entire catchment area of 4.5 m2, for these calculations. 12.5 c.c. /s is
equivalent to 1.00 cm /h as follows
1 4.5 106
____________
100 3600 12.5 c.c./s
112 ELEMENTARY ENGINEERING HYDROLOGY
1 2 3 4 5 6 7 8567 938 10
Mass Q mass
curve Mass Mass* Mass** SR mass Mass
rainfall curve curve SD curve SD mass curve F Infiltration
Sr. no. Time cm Q cm Q cm SR cm
cm cm cm rate cm
1 8–00 4 0.00 0.00 0.00 0.00 0 4.00 4.00
2 –10 8 0.00 0.00 0.00 1.15 1.15 6.85 2.85
3 –20 12 0.15 0.15 0.15 2.35 2.65 9.35 2.50
4 –30 16 0.40 0.55 1.50 2.35 4.40 11.60 2.25
5 –40 20 0.59 1.14 1.83 2.35 6.32 13.68 2.08
6 –50 24 0.80 1.94 4.08 2.35 8.37 15.63 1.95
7 9–00 28 1.03 2.97 5.20 2.35 10.52 17.48 1.85
8 –10 32 1.28 4.25 6.14 2.35 12.74 19.26 1.78
9 –20 36 1.54 5.79 6.90 2.35 15.04 20.96 1.70
10 –30 40 1.80 7.59 7.47 2.35 17.41 22.59 1.63
11 –40 44 2.00 9.59 7.88 2.35 19.82 24.18 1.59
12 –50 48 2.18 11.77 8.14 2.35 22.26 25.74 1.56
13 10–00 52 2.30 14.07 8.30 2.35 24.72 27.28 1.54
14 –10 56 2.38 16.45 8.40 2.35 27.20 28.80 1.52
15 –20 60 2.46 18.91 8.43 2.35 29.69 30.31 1.51
16 –30 64 2.50 21.41 8.43 2.35 32.19 31.81 1.50
17 –40 68 2.50 23.91 8.43 2.35 34.69 33.31 1.50
18 –50 72 2.50 26.41 8.43 2.35 37.19 34.81 1.50
19 11–00 72 1.36 27.77 6.37 2.35 36.49 35.51 0.70
20 –10 72 0.98 28.75 4.94 2.35 36.04 35.96 0.45
21 –20 72 0.76 29.51 3.88 2.35 35.74 36.26 0.30
22 –30 72 0.64 30.15 3.04 2.35 35.54 36.46 0.20
23 –40 72 0.53 30.68 2.36 2.35 35.39 36.61 0.15
24 –50 72 0.46 31.14 1.80 2.35 35.29 36.71 0.10
25 12–00 72 0.37 31.51 1.38 2.35 35.24 36.76 0.05
26 –10 72 0.34 31.85 1.04 2.35 35.24 36.76 0.00
27 –20 72 0.31 32.16 0.73 2.35 35.24 36.76 0.00
28 –30 72 0.20 32.36 0.53 2.35 35.24 36.76 0.00
29 –40 72 0.15 32.51 0.15 2.35 35.24 36.74 0.00
30 –50 72 0.00 32.51 0.00 2.35 35.24 36.74 0.00
The experiment was started at 8:00 a.m. by starting sprinkling water at the rate of 4.00 cm / h
over the entire catchment at a uniform rate. There was no flow of water upto 8:10 a.m. since the
requirement of SR [Interception and depression storage] was to be met with.
INFILTRATION 113
It was assumed plot SR is equivalent to 2.35 cm. The discharge started flowing at 8:20 a.m.
but was very low because the requirement of SR was not fully met with. It was assumed that the
requirement of SR of 2.35 cm. was completely met with upto 8:20 a.m.
The discharge as it started at 8:20 went on increasing upto 10:30 a.m. and reached 2.50 cm and
remained constant from 10:30 a.m. upto 10:50 a.m.. This indicated that the infiltration rate has
reached a steady state.
From 10:50 a.m., sprinkling of water was stopped.
Naturally the discharge started reducing from 10:50 a.m. and was completely stopped at 12:50 p.m.
Fig. 5.15. shows the discharge observed.
The storage capacity SD of the channel w.r.t. discharge flowing in the channel is calculated from
the falling graph of discharge from 10:50 a.m. onwards and is shown in Fig. 5.16.
4.0
Discharge in cm / h
3.0
2.0
1.0
0.0
8.00 9.00 10.00 11.00 12.00
Time in hours
9.0
8.0
7.0
6.0
Storage in cm
5.0
4.0
3.0
2.0
1.0
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0
Discharge in cm / h
4.0
Infiltration rate in cm / h
3.0
2.0
1.0
0.0
8.00 9.00 10.00 11.00 12.00 1.0
Time in hours
5.5.1 EQUATION
Infiltration rate is higher at the beginning. It goes on reducing w.r.t. time, and finally it attains a
steady rate as shown in Fig. 5.14.
Horton has suggested the following equation for the curve
f fc ( fo fc ) ekt
where, f Infiltration rate at time t
fo Initial infiltration rate
INFILTRATION 115
Example 5.3
The value of k in the Horton’s equation for infiltration is 2. The maximum and the minimum rates
of infiltration are 2 cm/h and 0.5 cm/h. Plot the infiltration rate curve.
Solution:
The Horton’s equation will be:
f 0.5 (2.0 0.5) e2t 0.5 1.52t and the hourly rate of infiltration will be as under:
2
lInfiltration rate in mm/hour
f0.5(2.0 0.5)e2t
1
0 1 2 3
Time in hours
5.5.2.1 index
The most common way of expressing infiltration rate is index.
Now, P E SR SD Q I
or, I P (E SR SD) Q
where, P Precipitation
E Evaporation
SR Surface retention (interception depression storage)
SD Surface detention (finally it flows as surface runoff)
Q Surface runoff
I Infiltration rate
Neglecting E, SR and SD, one can say that
index P Q
Indirectly, index includes all abstractions, e.g. evaporation, surface retention, and so on.
Figure 5.19 shows definition sketch of index.
60
Excess rainfall
50
Precipitation in mm/ hour
40
30
-index
20
10
0
1 2 3 4 5 6 7 8 9 10 11 12 13
Time in Hours
5.5.2.2 W index
W index is the rate of infiltration when all other abstractions are accounted, viz,
W index P (E SR SD Q)
Therefore, W index index
5.5.2.3 Wmin index
It is the minimum rate of infiltration when a uniform stage after stabilization is attained.
Wmin index after it is stabilized
k ndex,
i where k will be always less than one
5.5.2.4 fAV
This approach is a slight modification of index. In a storm, when there is no precipitation or
when the precipitation is very low, infiltration is still going on because of the previous high pre-
cipitation. A provision is, thus, made for this infiltration when precipitation is very low from the
previous precipitation when it is high. Consider a storm as shown in Fig. 5.20.
Rainfall exces
fav
No rain period
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13
Time in hours
A provision for infiltration when there is no precipitation, i.e. between 4 to 5 h, is made from the
storm precipitation between 3 to 4 h and hence is deducted from the excess rainfall from 3 to 4 h. fAV
is calculated and then the excess rainfall is worked out. It is slightly less than index.
Example 5.4
The average precipitation during a storm over a catchment area of 10 km2 is as follow:
1. 40 mm/h for 1 h
2. 60 mm/h for 1 h
3. 30 mm/h for 1 h
118 ELEMENTARY ENGINEERING HYDROLOGY
The resulting hydrograph was plotted on a graph paper with the following scale.
1 mc 1 h on x-axis
1 mc 10 m3/s on y-axis
If the area of the hydrograph was measured as 30 cm2, find the index of infiltration.
Solution:
Total runoff observed 30 10 3600
1.08 106 m3
Runoff from the storm assuming rate of infiltration to be
Example 5.5
The storm over a catchment of 50 km2 was having the following intensity:
1. 40 mm/h for 1 h
2. 70 mm/h for 2 h
3. 30 mm/h for 1 h
The catchment area had infiltration rate as follows:
a. 20% area 10 mm/h
b. 60% area 15 mm/h
c. Balance impervious
Find the runoff due to the storm.
Solution:
A) The runoff from the area having 10 mm/h will be
(7010) 2 ___________
(40 10) 1 ___________ (3010) 1
( ____________
1000 1000 1000 ) 0.2 50 10 6
30 120 20
_____________
1000 0.2 50 106 1.7 106 m3
INFILTRATION 119
25 110 15
_____________
1000 30 106 4.5 106 m3
70 2 ______
40 1 ______ 30 1
( ______
1000 1000 1000 ) 0.2 50 10
6
_____
210
1000 10 106 2.1 106 m3
Example 5.6
A catchment area having an average rate of infiltration of 15 mm/h experienced the storm of the
following intensities.
1. 50 mm/h for 2 h
__
1
2. 30 mm/h for 2 h
The resulting runoff was 10 106 m3.
Find the catchment area.
Solution:
(30 15) 0.5
(50 15) 2 ______________
[ ____________
1000 1000 ] A 10 10 10
6 6
70 7.5
Therefore, ( ________
1000 ) A 10
10 1000
_________
A 77.5 129.03 km2
Example 5.7
For a catchment area of 12 km2, a 7-h storm was as follows:
Ti (h) m e 1 2 3 4 5 6 7
Precipitation (mm) 20 40 0 30 50 40 5
120 ELEMENTARY ENGINEERING HYDROLOGY
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Discharge (m3/s) 0 7 17 32 66 55 45 36 28 21 15 10 6 3 1 0
Assume evaporation loss to be 0.6 mm/h/m2and seepage loss equal to 50% of the evaporation loss.
Find the index fAV index and W index.
Solution:
Total precipitation (20 40 0 30 50 40 5) 185 mm
∑ Discharge observed
(0 7 17 32 66 55 45 36 28 21 15 10 6 3 1 0) 342 m3/s
index
Assume index x mm/h
Assuming that the value of index will be more than 5 mm/s, the last rainfall slab of 5 mm/s
is neglected.
(20 x) 1 (40 x) 1 (30 x) 1 (50 x) 1 (40 x) 1
{ ________________________________________________________________
1000 } 12 10 6
342 60 60
3.6
________
342
⬖ 180 5x 12
342 3.6
_________
5x 180 12 180 102.6
__________
180102.6
Therefore, x 5 15.48 mm/h
fAV
Assume fAV y mm/h
Even though there is no precipitation between 2 and 3 h, there will be infiltration. This is taken
into consideration. So also the precipitation between 6 and 7 h is far less than the fAV (expected).
Infiltration during this period is also taken into consideration.
W index
Assume W index z mm/h
Evaporation loss seepage loss 1.5 5 0.6 4.5 mm
(180 4.5 5z)
_______________
1000 12 106 342 60 60
INFILTRATION 121
Example 5.8
The rate of precipitation in mm/h observed over a catchment of 30 km2 for successive 30 min is as
follows:
16, 20, 24, 36, 28, 12, 4
If the value of 22 mm/h, find
(1) Total precipitation and (2) Runoff in ha. m
Solution:
Since the precipitation observed is for a duration of 1/2 hour, the total precipitation will be as follows:
__
1 ____
140
2 (16 20 24 36 28 12 4) 2 70 mm
Since the rate of infiltration is higher than the first, second, sixth and seventh slab, all the precipi-
tation will be lost in the infiltration and the runoff will be only from third, fourth and fifth slab. It
will be as follows:
__
1
Runoff 2 [ (24 22)(36 22)(28 22)]
__
1
2 ( 2 14 6) 22 11 mm
11 30 106
_____________
1000 33 104 m3 33 ha m
REVIEW QUESTIONS
1. Define infiltration. Explain the process of infiltration.
2. Discuss the factors affecting infiltration.
122 ELEMENTARY ENGINEERING HYDROLOGY
3. Explain with the help of neat sketches the flooding-type infiltrometers. What are their advan-
tages and disadvantages?
4. Explain sprinkling-type infiltrometer.
5. Discuss the infiltration indices.
6. Write short notes on the following:
a. Dominant abstraction b. Disadvantages of tube infiltrometers
c. Horton’s equation of infiltration rate d. Estimation of index
e. Excess rainfall f. Limitations of sprinkling-type infiltrometer
7. Differentiate between the following:
a. Infiltration and percolation b. Single-tube infiltrometer and double-tube
infiltrometer
c. index and fAV index d. W index and Wmin index
NUMERICAL QUESTIONS
1. The quantity of water added to a double-ring infiltrometer with 1.0-m inside diameter at
30-min interval to keep the water level constant is as under:
Time (min) 0 30 60 90 120 150 180
Quantity of water added (l) 0 9.5 9.0 8.4 8.0 7.9 7.9
Find, (i) Rate of infiltration for every 30 min and plot the graph
(ii) Average rate of infiltration
Ans: (ii) 11.49 mm/h/m2
2. The value of K in the Horton’s equation for infiltration is 2.2. And the maximum and mini-
mum rates of infiltration are 2.1 cm/h and 0.5 cm/h. Plot the infiltration rate curve.
3. A catchment area of 13 km2 experiences a storm of average precipitation as: (1) 45 mm/h for
1 h, (2) 50 mm/h for 1.5 h, (3) 25 mm/h for 1 h and (4) 8 mm/h for 1 h.
The resulting hydrograph was plotted on a graph with the following scales:
1 mc 1 h on x-axis
1 mc 10 m3/s on y-axis
The area under the hydrograph thus plotted was measured and was found to be 40 cm2. Find
the index of infiltration.
Ans: 10 mm/h
4. The storm over a catchment of 60 km2 was having the intensity as: (1) 45 mm/h for 1 h,
(2) 60 mm/h for 2 h and (3) 35 mm/h for 1 h.
The catchment area had infiltration rates as under:
(a) for 25% area 10 mm/h
(b) for 50% area 15 mm/h
(c) Balance impervious
Find the runoff due to this storm.
Ans: 8.53 106 m3
INFILTRATION 123
5. For a catchment area of 12 km2, a storm of 6-h duration was observed as under:
Time (h) 1 2 3 4 5 6
Precipitation (mm) 22 41 0 32 55 6
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
3 0 8 19 34 68 58 48 40 25 19 15 11 6 3 0
Discharge (m /s)
Assume the evaporation loss to be 1 mm/day/m2 and the seepage loss equal to 50% of the
evaporation loss. Find, index, fAV and W index.
Ans: index 10.95 mm/h
fAV index 7.30 mm/h
W index 9.45 mm/h
*** ***
Hydrograph
6
Chapter Outline
g 6.1 DEFINITION
A hydrograph may be defined as a graphical representation of time versus discharge (Fig. 6.1).
6.1.1 UNITS
The unit for time on the x-axis may be:
1. Hours
2. Days
3. Months
126 ELEMENTARY ENGINEERING HYDROLOGY
This unit is selected depending upon the purpose and the nature of the study. If the study is
floods or flood routing, then the scale may be in hours. If the study is for estimation of runoff for a
year, then it may be in days or in months.
On the y-axis, the unit for discharge is m3/s. Sometimes, the discharge is expressed in cm/s,
i.e. the depth of water per unit area of the catchment area per second expressed in cm/s/m2.
Example 6.1
The scales for the hydrograph are
1 cm 1 h on x-axis
1 cm 10 m3/s on y-axis
If the area under the hydrograph is 45 cm2, (vide Fig. 6.2) find the runoff.
Solution:
Volume of water 45 (1 3600) (1 10) m3
1620 103 m3 1.62 106 m3 1.62 million m3 162 ha m
6.1.3 ISOLATED STORM
When subsequent storm does not occur before the runoff of the previous storm ceases, that storm
is called isolated storm (Fig. 6.3).
due to draining of this water stored in the river channel, i.e. draining of surface retention, and hence
the nature of the recession curve is always the same for a specific catchment. It depends only on the
basin characteristics of the catchment area and is, therefore, independent of the precipitation.
The curve from D to E is called recession curve.
An equation in the following form is derived for the recession curve:
Qt+1 Qt K (tt0)
where, Qt Discharge in m3/s at time t
Qt+1 Discharge in m3/s at time t0
K A constant known as recession constant
By observing recession curves of a number of storm runoffs of a specific catchment, a
master recession curve may be prepared with the help of a tracing paper. This is also known as
composite recession curve, master depletion curve or type curve. Point E denotes the end of sur-
face runoff due to the storm precipitation. The discharge in the stream after E is again the base
flow. It is rather difficult to locate the point E. (Theoretically, the recession curve is tangential
to the base flow at E.)
Thus, the storm hydrograph is from A to E and the time interval from A to E is known as base
period. It may be noted that the period for the rising curve is much less than that of the falling curve.
The flood hydrograph observed during the storm is a combination of the flood runoff and the base
flow. It is necessary to study the excess rainfall and its corresponding runoff from a catchment. Hence,
it is necessary to separate the base flow from the combined hydrograph to arrive at flood runoff.
When the flood occurs, the water level in the river rises substantially. The groundwater inflow,
if any, depends mainly on the difference between the level of groundwater table and water surface
level in the channel. It is a major portion of the base flow. When the water level in the channel rises,
this inflow may reduce during the floods and may be less than what it was at the beginning.
A typical hydrograph with base flow is shown in Fig. 6.6.
E since the recession curve and the base flow curve have different equations and hence will
have different slopes.
Join AE by a straight line.
The flow below AE is the base flow and above AE it is the flood runoff.
Example 6.2
Analyse the flood hydrograph (Fig. 6.14).
Solution:
The hydrograph is during an isolated storm. The base flow is separated from the flood hydro-
graph. The base flow goes on reducing and is negative for some period. It then increases and
becomes positive and reaches the point E.
134 ELEMENTARY ENGINEERING HYDROLOGY
The negative flow indicates that prior to the storm there is flow from groundwater to the river
channel. However, during floods the water level in the river channel is raised abruptly, and hence
water flows from river channel to groundwater, i.e. the base flow is negative.
The flood hydrograph will be as shown in Fig. 6.15.
When the base flow is separated from the combined hydrograph, the storm runoff will be arrived
at. The excess rainfall causing this storm runoff can be calculated as follows:
The hydrograph of the storm runoff is plotted and the volume of the storm runoff is calculated
by measuring the area under the storm hydrograph and also taking into consideration the scales
adopted for plotting the hydrograph for time and discharge on x- and y-axes as discussed in Sub-
section 6.1.2.
Assuming a value for index, the storm runoff is also calculated from the storm precipitation in terms
of and equating it to the storm runoff already calculated. The value of can thus be calculated.
Example 6.3
Calculate the excess rainfall when
P1, P2, P3 The storm rainfall in cm/h
3 cm/h; 5 cm/h; 4 cm/h
A Catchment area 50 km2
AH Area of storm hydrograph 100 cm2
136 ELEMENTARY ENGINEERING HYDROLOGY
Solution:
Assume index of infiltration x cm/h.
Thus,
(3 x) 1 (5 x) 1 (4 x) 1 50 106
______________________________________________ 100 3600 10
100
100 3600 100 10 ____
_____________________ 360
12 3x 50 7.2
50 106
Therefore, 3x 12 7.2 4.8
Therefore, x 1.6 cm/h
And excess rainfall will be 1.4 cm/h; 3.4 cm/h; 2.4 cm/h.
It may be noted that the excess rainfall and the flood hydrograph start at the same time
(Fig 6.16).
For identical precipitation, hydrographs observed for different catchments will have different
shapes. This is because the shape of the hydrograph of every catchment depends on the character-
istics of the catchment.
An attempt is made to correlate the hydrograph of each catchment with the precipitation.
This is done by the unit hydrograph theory.The theory of unit hydrograph was first presented
by L. K. Sherman in Engineering News Record in April 1932. Originally, this theory was known
HYDROGRAPH 137
as unit graph. However, this title unit graph was misinterpreted, and hence it was modified as
unit hydrograph.
Subsequent to the introduction of this theory by Sherman, it underwent a number of modifica-
tions, but the basic principle as presented by Sherman remained the same.
the unit precipitation is 1 cm and the period of precipitation is 1 h. (Naturally, the intensity of
precipitation is 1 cm/h.)
The surface runoff, in case of a 1 cm–1 h unit hydrograph from a catchment area A km2, will
be as below:
The surface runoff (1 cm/100) A 106m 3
( A/100) 106m 3
A 104m 3.
Example 6.4
The unit hydrograph coordinates of a 1 cm–1 h unit hydrograph are as follows:
Time (h) 0 1 2 3 4 5 6 7 8 9 10
3 0 6 13 22 16 11 7 4 2 1 0
Discharge (m /s)
Solution:
(A) For the required flood hydrograph, the time of precipitation is the same as that of unit
hydrograph, but the intensity of precipitation is double that of the unit hydrograph.
The hourly unit hydrograph coordinates are for a precipitation of 1 cm/h for 1 h. The
flood hydrograph coordinates for a precipitation of 2 cm/h for 1 h will be double the coor-
dinates of unit hydrographs, since the intensity of precipitation of the flood hydrograph
required is double the intensity of precipitation of the unit hydrograph.
The flood hydrograph coordinates will be as follows:
Time (h) 0 1 2 3 4 5 6 7 8 9 10
Discharge(m3/s) 0 12 26 44 32 22 14 8 4 2 0
The unit hydrograph and the flood hydrograph will be as shown in Fig. 6.18.
HYDROGRAPH 139
(B) By definition of the unit hydrograph, the precipitation over the entire catchment is 1.0 cm.
Therefore, the total expected runoff from catchment A would be A 106 1/100 m3.
A 104m 3
The runoff calculated from the unit hydrograph would be (hourly ordinates) 3600
These two figures should be equal
Therefore, A 104 82 3600
A 29.52 km2.
Example 6.5
Find the flood hydrograph from a catchment for excess rainfall of 2 cm/h for 1 h followed by
4 cm/h for 1 h followed by 3 cm/h for 1 h.
Given 1 h–1 cm unit hydrograph coordinates
Ti (h) m e 0 1 2 3 4 5 6 7 8 9 10
3 0 6 13 22 16 11 7 4 2 1 0
Discharge (m /s)
Solution:
The flood hydrograph of the excess precipitation mentioned will thus be a summation of three
flood hydrographs, viz.
140 ELEMENTARY ENGINEERING HYDROLOGY
i. The coordinates of the flood hydrograph due to 2 cm/h for 1 h will be as follows:
Time (h)
Serial no. 0 1 2 3 4 5 6 7 8 9 10
1 Unit hydrograph 0 6 13 22 16 11 7 4 2 1 0
coordinates (m3/s)
2 Flood hydrograph 0 12 26 44 32 22 14 8 4 2 0
coordinates (m3/s)
ii. The flood hydrograph coordinates for 4 cm/h for 1 h excess rainfall would be as follows:
Time (h)
Serial no. 0 1 2 3 4 5 6 7 8 9 10
1 Unit hydrograph 0 6 13 22 16 11 7 4 2 1 0
coordinates (m3/s)
2 Flood hydrograph 0 24 52 88 64 44 28 16 8 4 0
coordinates (m3/s)
Time (h)
Serial no. 0 1 2 3 4 5 6 7 8 9 10
1 Unit hydrograph 0 6 13 22 16 11 7 4 2 1 0
coordinates (m3/s)
2 Flood hydrograph 0 18 39 66 48 33 21 12 6 3 0
coordinates (m3/s)
This hydrograph will start 2 h later than the hydrograph due to 2 cm/h and 1 h later than the
hydrograph due to 4 cm/h.
Multiplications of unit hydrograph coordinates by the intensity of respective precipita-
tion are because of the assumed principle of linearity.
The combined hydrograph can be calculated by simple addition of these three hydrographs,
considering the proper time lag, as follows:
T x t 10 3 1 12 h
Example 6.6
Derive a flood hydrograph from a 1 cm–1 h unit hydrograph for the storm of excess rainfall of
2 cm/h for 1 h followed by 4 cm/h for 2 h then there is a gap of 1 h then 1 cm/h for 1 h.
The 1 cm–1 h unit hydrograph coordinates are as follows:
Time (h) 0 1 2 3 4 5 6 7 8 9 10
3
Discharge (m /s) 0 6 13 22 16 11 7 4 2 1 0
Solution:
The calculations for the flood hydrograph are as per the table below. It may be noted that the
precipitation of 4 cm/h for 2 h is split into two parts each of 4 cm/h for 1 h and their hydrographs
worked separately. The hydrograph due to 1 cm/h precipitation is taken after a gap of 2 h since
there is a gap of 1 h after the 4-cm/h precipitation.
142 ELEMENTARY ENGINEERING HYDROLOGY
1 2 3 4 5 6 7
3 3 3 3 3
Time (h) A (m /s) B (m /s) C (m /s) D (m /s) E (m /s) F (m3/s)
0 0 0 – – – 0
1 6 12 0 – – 12
2 13 26 24 0 – 50
3 22 44 52 24 – 120
4 16 32 88 52 0 172
5 11 22 64 88 6 180
6 7 14 44 64 13 135
7 4 8 28 44 22 102
8 2 4 16 28 16 64
9 1 2 8 16 11 37
10 0 0 4 8 7 19
11 – – 0 4 4 8
12 – – – 0 2 2
13 – – – – 1 1
14 – – – – 0 0
6.4.6 S HYDROGRAPH
Consider a 1 cm–1 h unit hydrograph as follows: (Fig. 6.21).
Time (h) 0 1 2 3 4 5 6 7 8 9 10
Discharge (m3/s) 0 6 13 22 16 11 7 4 2 1 0
Now imagine that the unit excess rainfall of the same intensity occurs for a period that tends to
infinity means the excess rainfall of that intensity occurs continuously up to infinite time.
This excess precipitation tending to infinity can be divided into time intervals equal to unit
hydrograph excess precipitation of 1 cm/h for 1 h one after another, tending to infinity.
1 2 3 4 5 6
1 cm/h 1 cm/h 1 cm/h 1 cm/h 1 cm/h 1 cm/h
→∞
144 ELEMENTARY ENGINEERING HYDROLOGY
→ⴥ
Naturally, each slab of excess precipitation will have one hydrograph, which is equivalent to
unit hydrograph one after another with a time lag of 1 h as follows:
The resulting hydrograph coordinates will be as shown in Table 6.1.
The S hydrograph coordinates will be as follows:
3. The unit hydrograph theory is not applicable for surface runoff originated from snow and ice.
4. The unit hydrograph theory is applicable to in-bank floods only. If the flood water overtops the
bank, this theory will not be applicable.
5. The unit hydrograph theory is applicable for catchments less then 5000 km2.
6. The theory is not applicable to narrow elongated catchments because it is not possible to have a
uniform precipitation over the entire catchment.
7. The theory is not applicable if there are storages on the channel or on its tributaries in the catch-
ment upstream of the gauging station.
At times, the time of excess precipitation in case of a unit hydrograph, keeping the unit precipita-
tion the same, is required to be changed depending upon the storm period t.
This time of excess precipitation may be increased or decreased.
A 2 h–1 cm unit hydrograph is required. (Intensity of precipitation will be 1/2 cm/h for 2 h.) This
can be derived by four methods as follows:
6.5.1.1 Method 1
The step-by-step procedure is as follows:
1. Plot the given 1 h–1 cm unit hydrograph.
2. Plot the same 1 h–1 cm unit hydrograph by a time lag of 1 h.
3. Add together these two flood hydrograph coordinates. Thus, resultant hydrograph after addition
will be a 2 h–2 cm hydrograph (Fig. 6.24).
4. Divide the resultant hydrograph by 2 to arrive at the required 2 h–1 cm hydrograph (The intensity of
precipitation will be 1/2 cm per hour for 2 h.) This procedure will be clear from Table 6.2.
The same procedure can also be followed graphically.
This new flood hydrograph will be 2 h–1 cm unit hydrograph with an intensity of 1/2 cm/h for
2 h, resulting into unit precipitation of 1 cm.
It should be noted that,
1. Maximum discharge has reduced.
2. The base period has increased by 1 h (2 1 1 h).
3. The total of the new unit hydrograph coordinates (2 h–1 cm) remains the same (82 m3/s).
HYDROGRAPH 147
Unit Hydrograph
coordinates Required unit
Unit Hydrograph after a lag of 1 h hydrograph
Serial no. Time (h) coordinates (m3/s) (m3/s) Total (m3/s) (m3/s)
1 0 0 0 0 0
2 1 6 0 6 3
3 2 13 6 19 9.5
4 3 22 13 35 17.5
5 4 16 22 38 19.0
6 5 11 16 27 13.5
7 6 7 11 18 9.0
8 7 4 7 11 5.5
9 8 2 4 6 3.0
10 9 1 2 3 1.5
11 10 0 1 1 0.5
12 11 0 0 0 0
3 3 3
Total 82 m /s 82 m /s 164 m /s 82 m3/s
6.5.1.2 Method 2
The same procedure can be followed by a different step-by-step procedure as follows:
1. Plot the given 1 h–1 cm unit hydrograph.
2. Divide these unit hydrograph coordinates by 2, so that a 1 h–1/2 cm hydrograph is obtained.
3. Plot this 1 h–1/2 cm hydrograph by a time lag of 1 h.
4. Add these two 1 h–1/2 cm hydrographs to arrive at the required 2 h–1 cm unit hydrograph.
148 ELEMENTARY ENGINEERING HYDROLOGY
1 0 0 0 0 0
2 1 6 3 0 3
3 2 13 6.5 3 9.5
4 3 22 11.0 6.5 17.5
5 4 16 8.0 11.0 19.0
6 5 11 5.5 8.0 13.5
7 6 7 3.5 5.5 9.0
8 7 4 2.0 3.5 5.5
9 8 2 1.0 2.0 3.0
10 9 1 0.5 1.0 1.5
11 10 0 0 0.5 0.5
12 11 0 0 0 0
Total 82 41 41 82
Use of S hydrograph
S hydrograph can also be used to increase the time of precipitation.
6.5.1.3 Method 3
The step-by-step procedure will be as follows:
1. Prepare an S hydrograph of 1 h–1 cm unit hydrograph and plot the S hydrograph.
2. Plot the same S hydrograph by a time lag of 2 h, say S1.
3. Deduct the second, i.e. S1 hydrograph from the first one S.
4. After deducting S1 hydrograph from S hydrograph, divide the resultant hydrograph by 2 to arrive
__
1
at the required 2 h–1/2 cm unit hydrograph. The intensity of precipitation will be 2 c m/h.
The procedure will be clear from Table 6.4.
The same procedure can also be followed graphically as shown in Fig. 6.26.
6.5.1.4 Method 4
The same procedure as followed in Section 6.5.1.3. in a different step-by-step procedure can be
followed as follows:
1. Prepare the S hydrograph of the given 1 h–1 cm unit hydrograph and plot it.
2. Divide these S hydrograph coordinates by 2 so that it will be a 1 h–1/2 cm S1 hydrograph. Plot
this S1 hydrograph with a time lag of 2 h.
3. Deduct this S1 hydrograph from the original S hydrograph to arrive at the required 2 h–1/2 cm
__
1
unit hydrograph. (The intensity of precipitation will be 2 cm for 2 h.)
The procedure will be clear from Table 6.5.
The same procedure can be followed graphically as shown in Fig. 6.27.
Example 6.7
Derive a 1 cm–5 h unit hydrograph from a 1 cm–3 h unit hydrograph.
The 1 cm–3 h unit hydrograph coordinates are as under:
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12
Discharge 0 2.00 6.34 13.66 17.00 16.32 11.34 7.34 4.34 2.32 1.00 0.33 0
(m3/s)
Solution:
The given unit hydrograph period 3 h and the required unit hydrograph period 5 h are not exact mul-
tiples of each other. Therefore, the required unit hydrograph can be derived only from the S curve.
The step-by-step solution will be as follows: vide Table 6.6
1. Write down the given unit hydrograph coordinates. (Since the given unit hydrograph is 1 cm–3 h,
__
1
the intensity of the precipitation is 3 cm/h.) (column A)
2. Repeat these coordinates each after 3 h to obtain the S hydrograph. (Even though the given unit
hydrograph is of 3 h, the unit hydrograph coordinates are at 1-h interval.) (columns B, C, D, E, F)
6.5.2.1 Method 1
The step-by-step procedure will be as follows:
1. Prepare an S hydrograph from the given 2 h–1 cm unit hydrograph and plot it.
2. Plot this S hydrograph with a time lag of 1 h, say S1 hydrograph and plot it.
3. Deduct this S1 hydrograph from S hydrograph. After deduction, multiply it by 2 to arrive at the
required 1 h–1 cm unit hydrograph (Fig. 6.28).
0 0 0 0 0
1 3.0 3 0 6
2 9.5 9.5 3 13
3 17.5 20.5 9.5 22
4 19.0 28.5 20.5 16
5 13.5 34.0 28.5 11
6 9.0 37.5 34.0 7
7 5.5 39.5 37.5 4
8 3.0 40.5 39.5 2
9 1.5 41 40.5 1
10 0.5 41 41 0
Total 82 82
6.5.2.2 Method 2
The same procedure can be followed in a different step-by-step procedure as below:
1. Prepare the S hydrograph from the given 2 h–1 cm unit hydrograph and plot it.
2. Multiply the coordinates of this S hydrograph by 2 so that a 2 h–2 cm S1 hydrograph is obtained.
3. Enter this S1 hydrograph with a time lag of 1 h.
4. Deduct the S1 hydrograph from the S hydrograph to arrive at the required 1 h–1 cm unit hydro-
graph (Fig. 6.29).
Example 6.8
Derive a 1 cm–3 h unit hydrograph from a 1 cm–4 h unit hydrograph.
The 1 cm–4 h unit hydrograph coordinates are as under:
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12 13
Discharge 0 1.50 4.75 10.25 14.25 15.50 14.00 9.50 6.00 3.50 1.75 0.75 0.25 0
3
(m /s)
Solution:
The given unit hydrograph period 4 h and the required unit hydrograph period 3 h are not exact
multiples of each other. Therefore, the required unit hydrograph can be derived only from the
S curve.
The step-by-step solution will be as follows: vide Table 6.9
1. Write down the given unit hydrograph coordinates. (Since the given unit hydrograph is
__
1
1 cm–4 h, the intensity of the precipitation is 4 cm/h.) (column A)
2. Repeat these coordinates each after 4 h to obtain the S hydrograph. (Even though the given unit
hydrograph is of 4 h, the unit hydrograph coordinates are at 1-h interval.) (columns B, C, D)
HYDROGRAPH 155
Unit hydrograph may be derived from the observed rainfall and its resultant hydrograph. For this,
the available data will have to be scanned before it is used.
The step-by-step procedure is as follows:
1. From the observed data, select an isolated storm and its resultant hydrograph.
2. Check the precipitation data for the accuracy of any rain gauge station. If required, evalu-
ate the missing data. Find the average hourly precipitation over the catchment and plot the
hyetograph.
If the unit hydrograph required is for t hours, then the hyetograph may be divided into
slabs of t hours. Assuming a index for infiltration, evaluate the average hourly excess
rainfall hyetograph.
3. Check the resultant hydrograph for its accuracy. Separate the base flow. Then the correct
observed flood hydrograph will be arrived at.
Then the observed corrected data of excess precipitation and its resultant storm hydrograph
can be used for further analysis.
156 ELEMENTARY ENGINEERING HYDROLOGY
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12
3 0 12 50 114 159 152 106 69 41 22 10 3 0
Discharge (m /s)
Q0 0
Q1 2U1 12
Q2 2U2 4 U1 50
Q3 2U3 4 U2 3 U1 114
Q4 2U4 4 U3 3 U2 159
Q5 2U5 4 U4 3 U3 152
Q6 2U6 4 U5 3 U4 106
Q7 2U7 4 U6 3 U5 69
Q8 2U8 4 U7 3 U6 41
Q9 2U9 4 U8 3 U7 22
Q10 2U104 U9 3 U8 10
Q11 2U114 U10 3U9 3
Q12 2 U124 U11 3U10 0
Q0 0 U0 U0 0
Q1 12 2U1 U1 6
Q2 50 2U2 4 U1 2U2 4 6 U2 13
Q3 114 2U3 4 U2 3 U1 2U3 4 13 3 6 U3 22
Q4 159 2U4 4 U3 3 U2 2U4 4 22 3 13 U4 16
Q5 152 2U5 4 U4 3 U3 2U5 4 16 322 U5 11
Q6 106 2U6 4 U5 3 U4 2U6 4 11 3 16 U6 7
Q7 69 2U7 4 U6 3 U5 2U7 4 7 3 11 U7 4
Q8 41 2U8 4 U7 3 U6 2U8 4 4 3 7 U8 2
Q9 22 2U9 4 U8 3 U7 2U9 4 2 3 4 U9 1
Q10 10 2U10 4U9 3U8 2U10 41 3 2 U10 0
Q11 3 2U11 4U10 3U9 31 U11 0
From the first equation, we get Q0 0. Substitute its value in the next equation and we get
Q1 6. Thus, by substituting the unit hydrograph coordinate, calculated in the previous step, in the
next step, all the unit hydrograph coordinates can be evaluated.
158 ELEMENTARY ENGINEERING HYDROLOGY
The procedure was first suggested by Prof. Dr Ven Te Chow, hence this method is known as Ven
Te Chow’s forward substitution method.
This backward substitution method was also suggested by Prof. Dr Ven Te Chow, and hence it
is known as Ven Te Chow’s backward substitution method.
For both the methods, the unit hydrograph is checked so that the area of the hydrograph is one
unit for unit excess precipitation uniformly distributed over the basin. The unit hydrograph thus
derived by both the methods is shown in Fig. 6.31.
There are a number of assumptions made in the unit hydrograph derivation. So also there might
be some errors in the observed precipitation and the discharge data, and hence unit hydrographs
derived from the different observed data for the same catchment may not tally with each other.
Thus, an average of all the unit hydrographs derived is worked out for further studies.
The procedure followed is as stated below:
Plot all the unit hydrographs on a simple graph paper
t1t2t3…tn
______________
tp= n
where, t1, t2, t3 Time in hours of the maximum discharge from the beginning of the differ-
ent unit hydrographs derived
tp Time in hours of the maximum discharge from the start of the final average
unit hydrograph
Q1, Q2, Q3 Maximum flood discharge of the various unit hydrographs in m3/s
Q Maximum flood discharge of the final average unit hydrograph
T1T2T3…Tn
________________
T n
Example 6.9
The unit hydrographs derived from different storms for a catchment are as follows:
Derive the average unit hydrograph.
324
Average tp ________ 3 h
3
8 10 9
__________
Average T 3 9h
23 22 21
Qm ____________
3 22 m3/s
Example 6.10
The peak of the flood hydrograph resulting from a 3-h 2.5-cm/hour storm is 320 m3/s. If the
value of infiltration is 5 mm/s, find the peak value of the 3 cm–3 h unit hydrograph. Assume the
base flow to be 20 m3/s.
Solution:
After deducting the base flow, the peak discharge of the flood hydrograph resulting due to the
storm only will be 320 20 300 m3/s.
The excess uniform precipitation after deducting the infiltration rate will be 2.5 0.5 2 cm/h
for 3 h. The excess precipitation for a 3 cm–3 h unit hydrograph will have to be 1 cm/h for 3 h.
__
1
Therefore, the peak discharge for a 3 cm–3 h unit hydrograph will be 2 300 150 m3/s.
When sufficient observed data are not available for the derivation of a unit hydrograph, then the
unit hydrograph is framed based on the catchment characteristics.
tp Time between the CG of effective rainfall to the peak discharge in hours (basin lag in
hours)
Qp Peak discharge in m3/s
T Base period in hours
W75 Width of unit hydrograph for 75% of Qp
W50 Width of unit hydrograph for 50% of Qp
Ct A constant depending upon the slope and SD of the catchment area. (Value is normally
between 1.35 to 1.65.)
Cp A coefficient ranging from 0.56 to 0.69
tp Ct (L·Lc)0 3
Qp 2.778 (Cp /tp) A
T 5.455 tp
W50 2.14 (Qp /A)1 08
W75 1.22 (Qp/A)1 08
t 2 tp/11
The widths of unit hydrographs W50 and W75, respectively, may be divided into two parts such
that 1/3 lies in the rising curve and 2/3 in the recession curve.
The unit hydrograph can be framed, based on the above set of equations. As a check, the area of
the unit hydrograph may be checked and corrected for a value of unit precipitation (Fig. 6.35).
Example 6.11
Derive a unit hydrograph for the catchment area shown in Fig 6.36
Solution:
Here A 250 km2, L 20 km, Lc 5.0 km
To suit the characteristics of the catchment, the values of the of constants assumed are Cp 0.6,
Ct 1.50.
HYDROGRAPH 163
Example 6.12
Derive a synthetic unit hydrograph for a catchment area having a shape of an isosceles triangle of
base width 6 km and height 10 km with the main stream along the medium bisecting the base.
Assume Ct 1.5 and Cp 0.6
Solution:
__
1
Here, i. Area of catchment 2 10 6 = 30 km2
ii. Length of the main stream 10 km
HYDROGRAPH 165
iii. Length of the main stream from a point nearest to the CG up to outlet 6.67 km
From these characteristics, the following parameters can be derived:
tp Ct (L Le)0 3 1.5 (10 6.67)0 3 5.28h 5.3 h
Qp 2.778 (Cp/tp) A 2.778 (0.6/5.0) 30 10.31 m3/s
T 5.455 tp 5.455 5.0 27.275 27 h
t 2 tp/11 2 5/11 0.909 1.0 h
W50 2.14 (Qp/A)1 08 2.14 (10.31/30)1 08 6.78 h
W75 1.22 (Qp/A)1 08 1.22 (10.31/30)1 08 3.86 h
The shape of the catchment area and the unit hydrograph thus derived will be as shown is Fig. 6.38.
The area under the unit hydrograph measured 42 cm2
Volume of unit hydrograph 42 2 3600 1 30.24 104m 3
This should be the volume due to unit uniform precipitation over the catchment.
x A 10 6
___________
Therefore, 100 30.24 104
x 1.008 1.0 cm
Therefore, the unit precipitation for the unit hydrograph is 1 cm.
Qp Maximum discharge
t Time when Qt is observed from the beginning of the unit hydrograph
tp Time when the maximum discharge is observed
The shape of the unit hydrograph thus derived is in agreement of the hydrograph expected for
the catchment. For averaging the unit hydrographs derived for various catchments, this method is
more effective. The unit hydrograph derived for the neighbouring catchments can be converted into
a dimensionless form and transposed to an ungauged catchment.
The dimensionless unit hydrograph derived for a basin is shown in Fig. 6.39.
t/tp 0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
Q/Qp 0 0.43 1.0 0.66 0.32 0.16 0.075 0.036 0.018 0.009 0
Example 6.13
The unit hydrograph coordinates of a basin are as under. Prepare the distribution graph.
Time (h) 0 1 2 3 4 5 6 7 8 9 10
Discharge (m3/s) 0 6 13 22 16 11 7 4 2 1 0
Solution:
The percentage runoff during each hour is calculated as under. The table is self-explanatory.
1 2 3 4 5
Unit hydrograph
Serial no. Time (h) Average runoff (m3/s) Runoff (%)
(m3/s)
1 0 0 – –
2 1 6 3 3.66
3 2 13 9.5 11.59
4 3 22 17.5 21.34
5 4 16 19.0 23.17
6 5 11 13.5 16.46
7 6 7 9.0 10.98
8 7 4 5.5 6.70
9 8 2 3.0 3.66
10 9 1 1.5 1.83
11 10 0 0.5 0.61
Total 82 82 100
Example 6.14
The distribution graph coordinates as stated below are to be applied to a neighbouring catchment cov-
ering an area of 360 km2. Evaluate the flood hydrograph when the excess rainfall is 2.0 cm/h for 2 h.
Time (h) 0 1 2 3 4 5 6 7 8 9 10
Distribution graph (%) 0 4 11 21 23 16 11 7 4 2 1
Solution:
The base period of the unit hydrograph 10 h.
If precipitation of p cm/h occurs uniformly distributed over a catchment of A km2, then the
uniform discharge may be Q.
Therefore,
p A 106 2.77 p A
Q____ _______
100 3600
____
x
Therefore, discharge with x% distribution = 100 2.77 p A
0.0277 x p A
The flood hydrograph due to the precipitation of 2 cm/h for 2 h is tabulated as under.
For simplicity, sometimes, a unit hydrograph is expressed in the form of a triangle with its
base parallel to the x-axis. In this case, the base of the triangle represents the base period and
the height of the triangle represents the peak discharge of the unit hydrograph as shown in
Fig. 6.42.
Example 6.15
The 1 h–1 cm unit hydrograph of a catchment is in the form of a triangle of base width of 15 h and
peak discharge of 80 m3/s occurring after 5 h from the start. Find the catchment area.
Solution:
The unit hydrograph diagram will be as shown in Fig. 6.43.
Let the catchment area Akm 2
i. The total runoff will be the area under the graph.
Area under the graph __12 15 80 600 m3/s 1 h 600 m3/s 3600 s
216 104m 3.
i. The runoff from the catchment area A km2 from a precipitation of 1 cm/h for1h
(A 106) 1/100 A 104m 3
170 ELEMENTARY ENGINEERING HYDROLOGY
The characteristics of a unit hydrograph such as base time, maximum discharge, etc. depend upon
the following factors:
• Catchmentc haracteristics
• Intensity of the unit precipitation
• Time of unit precipitation
In order to simplify the unit hydrograph, it is assumed that the unit precipitation occurs uni-
formly over the entire catchment instantaneously, i.e. t 0 or t → 0.
It is a conceptual fictitious unit hydrograph.
This is a hypothetical case. It is incorrect to assume that the unit precipitation occurs instantly.
However, further analysis of unit hydrograph theory is simplified. With this assumption, one factor
(time of unit precipitation) influencing unit hydrograph is eliminated. Thus, a unit hydrograph of a
catchment area, resulting from a unit precipitation occurring instantly, uniformly distributed over
the entire catchment, is known as instantaneous unit hydrograph as shown in Fig. 6.44.
It is normally denoted as IUH and its shape is similar to a single peak hydrograph.
It may be noted that the basic principles of linearity and time invariance, assumed in the unit
hydrograph theory, hold good for IUH also.
Indirectly, a unit hydrograph of time of precipitation of t → 0 is an IUH.
REVIEW QUESTIONS
1. Define hydrograph. Discuss the normal units used to draw a hydrograph.
2. What is a complex storm? Sketch a complex storm and the resulting hydrograph.
3. Sketch a hydrograph resulting from an isolated storm and explain its features.
4. What is base flow? Why is it necessary to separate it from the hydrograph? Explain the dif-
ferent methods to separate it.
5. Explain with a neat sketch as to how the base flow is separated from a hydrograph resulting
from a complex hydrograph.
6. Define unit hydrograph. Discuss the different definitions involved in this theory.
7. Explain the assumptions made in the unit hydrograph theory.
8. Explain how the unit hydrograph is used to find the flood hydrograph.
9. Discuss the limitations of unit hydrograph theory.
10. What is an S hydrograph? How is it derived from a unit hydrograph?
11. How is an S hydrograph used to change the period of the unit hydrograph?
12. Discuss the different methods followed to increase the time of the unit hydrograph.
13. Discuss the different methods to derive a unit hydrograph from observed data.
14. What is averaging of unit hydrograph? Why is it necessary?
15. What is a synthetic unit hydrograph? Why is it necessary?
16. Explain Snyder’s method of framing a synthetic unit hydrograph.
17. Explain distribution graph.
18. Explain dimensionless unit hydrograph.
19. Write a detailed note on instantaneous unit hydrograph.
20. Write short notes on:
a. Excessr ainfall b. Masterde pletionc urve
c. Specifications of unit hydrograph d. Collin’s method of derivation of unit
hydrograph
e. Ven Te Chow’s forward substitution f. Ven Te Chow’s backward substitution
method method
g. Dimensionlessuni thydr ograph h. Distributiongr aph
i. Principle of linearity j. Principle of time invariance
k. Time of concentration
172 ELEMENTARY ENGINEERING HYDROLOGY
NUMERICAL QUESTIONS
6.1. The flood hydrograph co-ordinates resulting from an isolated storm and the base flow
co-ordinates after separation are as under. Study and discuss the hydrograph.
Time in Hours 0 1 2 3 4 5 6 7 8 9 10
Discharge in m3/s 10 14 18 24 16 13 1 10 9 10 10
6.2. The unit hydrograph co-ordinates of a 1 hour 1 cm unit hydrograph are as follows
Time in Hours 0 1 2 3 4 5 6 7 8 9 10 11 12
Discharge in m3/s 0 5 13 25 40 29 20 13 8 3 2 1 0
Time
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
in Hours
Discharge
5 25 92 331 586 594 450 348 225 130 70 42 18 9 7 5
in m3/s
Time in Hours 0 1 2 3 4 5 6 7 8 9 10 11 12 13
Discharge in m3/s 0 2.5 9.0 19.0 32.5 34.5 24.5 16.5 10.5 5.5 2.5 1.5 0.5 0
Derive a 3 hours, 1cm, unit hydrograph. Calculate also the catchment area.
Ans: The 3 hours 1 cm, unit hydrograph co-ordinates will be as follows.
Time
0 1 2 3 4 5 6 7 8 9 10 11 12 13
in Hours
Discharge
0 1.67 6.00 14.34 26.00 31.34 29.67 20.67 13.67 8.00 4.34 2.00 1.00 0.34
in m3/s
Time in Hours 0 1 2 3 4 5 6 7 8 9 10 11 12 13
3
Discharge in m /s 0 2.5 9.0 19.0 32.5 34.5 24.5 16.5 10.5 5.5 2.5 1.5 0.5 0
Time in Hours 0 1 2 3 4 5 6 7 8 9 10 11 12
Discharge in m3/s 0 5 13 25 40 29 20 13 8 3 2 1 0
Time
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
in Hours
Discharge
5 19 79 191 352 472 430 300 198 122 64 34 19 9 5 5
in m3/s
Time in Hours 0 1 2 3 4 5 6 7 8 9 10 11 12
Discharge in m3/s 0 5 13 25 40 29 20 13 8 3 2 1 0
6.6. Derive a 1 hour, 1cm unit hydrograph for a catchment with the following details.
i) Catchment area ____________________________ 1 00 km 2
ii) Length of main stream _______________________ 18 km
iii) Distance along the main stream from the outlet upto a point nearest to the C. G. of the
catchment area. ______________________________ 5 km.
Assume Cp and Ct in Snyder’s method to be ‘0.5 and 1.4 resp.
Ans: The required 1 hour, 1 cm unit hydrograph details will be as follows
(1) Maximum discharge :— Qp 2 42 m3/s,
(2) Time of maximum discharge :— tp 5 4 hours.
(3) Base period :— T 5 22 hours.
(4) Time of unit precipitation :— 1 hour.
(5) Q50 5 21 m3/s, W50 5 5.46 hours.
(6) Q75 32 m3/s, W75 3.11 hours.
6.7. A 1 hour, 1cm unit hydrograph co-ordinates are as under
Time in Hours 0 1 2 3 4 5 6 7 8 9 10 11 12
Discharge in m3/s 0 5 13 25 40 29 20 13 8 3 2 1 0
174 ELEMENTARY ENGINEERING HYDROLOGY
Time
0 1 2 3 4 5 6 7 8 9 10 11 12
in Hours
Run-off
0 1.57 5.66 11.94 20.44 21.69 15.40 10.30 6.60 3.46 1.57 0.94 0.31
in Percentage
6.8. The 1 hour, 1 cm unit hydrograph of a catchment area is in a triangular shape. The details as
follows
(i) Catchment a rea ________________________________ 25 km 2
(ii) Base period ___________________________________ 10 hours
(iii) Maximum discharge ____________________________ 21 m3/s
(iv) Time of max discharge ___________________________ 3 hours
Calculate the flood hydrograph for the following storm.
4 cm / hour for 1 hour followed by 7 cm / hour for 2 hours followed by 2 cm/hour for 1 hour.
Assume a uniform base flow of 5 m3/s
Ans: The flood hydrograph will be as under
Time in Hours 0 1 2 3 4 5 6 7 8 9 10 11 12 13
Discharge in m3/s 5 33 110 236 336 366 326 266 206 146 86 38 11 5
17. In a dimensionless unit hydrograph, the base period ratio is generally taken to be
(a) 3 (b) 4
(c) 5 (d) 6
176 ELEMENTARY ENGINEERING HYDROLOGY
21. The distribution graph for a catchment is practically the same irrespective of
(a) Intensity of precipitation (b) Duration of precipitation
(c) Catchmenta rea (d) Characteristicsof c atchment
22. In the instantaneous unit hydrograph, the effect of the following parameter is
eliminated.
(a) Intensity of precipitation (b) Duration of precipitation
(c) Catchmenta rea (d) Characteristicsof c atchment
23. Principle of linearity means that the coordinates of the flood hydrograph are directly proportional to
(a) Intensity of precipitation (b) Duration of precipitation
(c) Catchment area (d) All of the above
25. In Snyder’s method, the widths of the unit hydrograph for 75% and 50% of maximum dis-
charge are divided in the rising and falling curves in the ratio of
(a) 1:1 (b) 1:2
(c) 1:3 (d) 1:4
*** ***
Runoff
7
Chapter Outline
g 7.1 DEFINITION
‘Runoff’ is that part of precipitation that appears in a drainage channel as surface flow in a peren-
nial or an intermittent form. It is that part of water, which can be used for engineering purposes and
hence is also known as yield of catchment.
The yield from a catchment is generally expressed in terms of volume, in a season or a year.
178 ELEMENTARY ENGINEERING HYDROLOGY
Catchment area may be defined as the area from which the surface runoff is derived. It is also
known as the watershed area, drainage area, drainage basin or simply basin or catchment.
RUNOFF 179
The unit for the catchment area is km2. If the catchment area is less than 25 km2, it is mentioned in
terms of hectares. The catchment areas formed by the divide lines at A and B are shown in Fig. 7.3.
The smallest drainage channels near their origin are grouped under Order 1. When two channels
of Order 1 meet, a channel of Order 2 is formed. Similarly, when two channels of Order 2 meet, a
channel of Order 3 is formed and so on. This ordering of stream is shown in Fig. 7.4. Thus, the order
of the main stream indicates the extent of branching in the catchment area. It is dimensionless.
Example 7.1
The details of a catchment collected from a toposheet are as follows.
Here the scale of the map is 1:50,000.
The required data were collected from the map and are given as under.
1. Area of the catchment: 28 cm2 (on the map) 7.0 km2 (actual)
2. Length of the main stream: 12 cm (on the map) 6 km (actual)
3. Length of the ridge line: 27 cm (on the map) 13.5 km (actual)
4. Level of the highest point on the ridge: 325 m
5. Level of the lowest point, i.e. at the outlet: 230 m
6. The order of all the streams was marked on the map and their lengths were measured as follows:
Solution:
Taking into consideration the scale of the map, 1:50,000, the corresponding prototype values
are also mentioned.
The diameter of a circle whose perimeter is equal to the length of the ridge line
D 13.5
____
13.5
Therefore, D 4.29 km
D2
____
The diameter of a circle whose area is equal to the catchment area 4 7.0
Therefore, D ___ ( )
28 1/2 2.98 km
The characteristics of the catchments will be as follows:
S S S ...S
S[ ]
1/2 1/2 1/2 1/2 2
______________________
1 2 3 n
Thus, n
where, S Slope of the catchment area
S1, S2, S3, ... Slope of the main stream in various segments
n Number of segments
Example 7.2
The L section along the main stream of a channel is shown in Fig. 7.6. Find the average slope of
the channel.
RUNOFF 183
Solution:
Divide the L section into eight equal horizontal segments. Each will be of 10 km. The bed slope
of each segment is as follows:
(a) 1 in 10,000 (b) 3 in 10,000 (c) 5 in 10,000 (d) 6 in 10,000 (e) 8 in 10,000 (f) 13 in 10,000
(g) 21 in 10,000 (h) 38 in 10,000
The average slope S of the channel will be
S11/2S21/2S31/2...
[ ]
2
S _________________
N
]
2
8
1/2 1/2 1/2 1/2 1/2
The catchment areas are classified depending upon the shape as follows (Fig. 7.7):
1. Fan shape
2. Fern leaf-type or elongated
Fan shape: The shape of a catchment may be similar to a fan, i.e. circular. This will affect the
flood intensity from the catchment.
Fern leaf-type: The shape of a catchment may be similar to a fern leaf, i.e. elongated. This will
affect the flood intensity from the catchment.
Karst topography
Regions underlain by soluble rock formations like limestone have characteristic undulating surfaces
with conical round hillocks and circular sinks. Such regions are said to have Karst topography
and are known as Karstic regions. The runoff from such areas is excluded from the surface runoff
RUNOFF 185
estimation, since the precipitation over such areas may be lost into the underground passages to far
distant streams, lakes or even sea.
It k I t1
The value of k ranges from 0.85 to 0.98. The API is taken into consideration while assessing the
runoff or flood from a catchment area.
Whenever there is precipitation, the requirement of interception and depression storage is first met
with. The infiltration also starts immediately. If the rate of precipitation is more than these abstrac-
tions, overland flow results. Otherwise, all precipitation water is lost in these abstractions. This
overland flow is collected in small streams leading to a channel or a river.
Before and after this storm, there may be some discharge in the river. This may be due to inter-
flow, delayed interflow or flow from the groundwater appearing into the channel. The surface flow
observed in the river immediately after the precipitation is the combined flow due to the precipita-
tion and the flow in the river before the precipitation.
RUNOFF 187
The combined effect of climate, soil structure and geology of the catchment area is noticed in the
network of channel or stream pattern. The usual patterns observed are as follows (Fig. 7.10):
• Dendrite type or tree like
• Radial pattern
• Trellis pattern
• Pinnate type
• Anastomising pattern
• Parallel
• Rectangular
These stream patterns have a significant effect on the draining period and the flow intensity of
the surface runoff and indirectly on infiltration and total runoff.
7.8.1 PRECIPITATION
The runoff is clearly a function of precipitation, its intensity, its duration and its coverage. More
the intensity, more will be the runoff. The infiltration rate reduces after some time; hence more the
duration, proportionately more will be the runoff. Similarly, more the area covered by the storm,
more will be the runoff.
Direction of movement of a storm over the catchment area has a definite effect on the
runoff. If the storm moves in the direction of the flow, the base period of hydrograph will
be less and more peak flow may be expected. On the other hand, if the storm moves against the flow direc-
tion, then the base period will be comparatively more and less peak flow may be expected (Fig. 7.11).
In case of a fan-shaped catchment area, the base period of the resulting hydrograph will be less
and thus more peak flow may be expected. In case of an elongated catchment, the base period of
the resulting hydrograph will be comparatively more and thus more will be the infiltration losses
and less will be the runoff (Fig. 7.12).
The average value of (1) Manning’s coefficient, (2) the shape factor, (3) stream order, (4) drainage
density, (5) stream density, (6) circularity ratio, (7) elongation ratio, (8) compactness coefficient, (9) the
slope of the channel and so on of the catchment thus definitely influence the runoff.
190 ELEMENTARY ENGINEERING HYDROLOGY
The sources and components of runoff from a catchment area are as follows:
7.9.3 INTERFLOW
All the water, which infiltrates into the ground, may not reach up to the groundwater table
if it meets a local impervious layer. This water may move laterally as per the slope of this
impervious local layer and may find an outlet into the stream. This component of runoff is
known as interflow.
192 ELEMENTARY ENGINEERING HYDROLOGY
If this inflow is observed in a short period, it is called prompt interflow; otherwise it is called
delayed interflow. This inflow moves slowly as compared to the overland flow and its contribution
depends upon the orientation of the impervious layer as well as the soil characteristics. This is also
known as subsurface runoff, subsurface storm flow, storm seepage, secondary base flow, through
flow or quick return flow.
7.9.5 REGENERATION
When water is supplied to the fields for irrigation, a small portion of it is utilized for the growth of
the crops. Remaining portion infiltrates. A part of the infiltrated water finds its way into the adjoin-
ing stream and is called regeneration. The remaining flows to the groundwater.
For the regenerated water, the lag time is small. Normally, regeneration is considered to be 20%
of the water applied for irrigation.
The observations covered areas that were not having heavy rainfall but were having a precipita-
tion up to 1100 mm. Hence, these percentages are not applicable for high-rainfall areas.
Example 7.3
The average annual precipitation over a catchment of 120 km2 and its percentage is as follows:
Serial no. % of area Average annual precipitation (cm)
A1 20 50
A2 15 75
A3 40 90
A4 25 110
Solution:
The Binnie’s coefficients corresponding to the average annual precipitation for each area and
the corresponding runoff will be as follows:
1 Urban 0.05--0.5
2 Forest 0.2--0.5
3 Commercial and industrial 0.9
4 Parks and farms 0.3--0.5
5 Concrete pavements 0.85
Example 7.4
A catchment area of 150 km2 has the following type of distribution:
Area Type %
A1 Urban 20
A2 Forest 35
A3 Commercial and industrial 30
A4 Concrete pavement 15
Find the annual runoff from the catchment when the average annual precipitation is 90 cm.
Solution:
The runoff coefficient K for the respective type of catchment and the annual runoff from each
type of catchment will be as follows:
A1 150 0.20 30.0 0.25 30.0 106 0.9 0.25 6.75 106
A2 150 0.35 52.5 0.45 52.5 106 0.9 0.45 21.26 106
A3 150 0.30 45.0 0.90 45.0 106 0.9 0.90 36.45 106
A4 150 0.15 22.5 0.85 22.5 106 0.9 0.85 17.21 106
Total 81.67 106
1 6.25 – – 8
2 12.50 – 6 12
3 25.00 3 11 18
4 37.50 6 16 25
5 50.00 10 22 34
6 75.00 20 37 55
7 100.00 30 50 70
Example 7.5
The daily rainfall over a catchment of 140 km2 is 5, 30 and 80 mm, respectively. The area is
classified as follows:
1. Dry—30%
2. Damp—45%
3. Wet—25%
Find the runoff from this storm using Strange’s coefficient. The area is classified as good catchment.
Solution:
The rainfall and the corresponding Strange’s coefficient for the classification of the catchment
will be as follows:
The runoff, in mm, for each day for the three classifications will be as follows:
Runoff (mm)
Serial no. Day Rainfall (mm) Dry Damp Wet
1 1 5 -- -- 0.08 5 0.40
2 2 30 0.035 30 1.05 0.13 30 3.9 0.20 30 6.00
3 3 80 0.230 80 18.40 0.40 80 32.0 0.62 80 49.60
Total 19.45 35.9 56.00
196 ELEMENTARY ENGINEERING HYDROLOGY
Since the area is classified as ‘good’, the runoff from the catchment will be increased by 25%
and the runoff from the catchment will be 1.25 5.039 106
6.30 10 6 m3
Table 7.4 Percentage coefficient of small catchment areas depending on the nature of the catchment
Serial no. Class Nature of catchment Runoff (%)
1 A Flat, cultivated black cotton soil area 10
2 B Flat partly cultivated consisting of different soils 15
3 C Average land 20
4 D Hills and plains with little cultivation 35
5 E Very hilly and steep slopes with hardly any cultivation 45
Example 7.6
A catchment covering an area of 200 km2 has the following classification:
A1. Flat cultivated land 20%
A2. Hills and plain with little cultivation 40%
A3. Very hilly and steep slopes with no cultivation 40%
RUNOFF 197
The average precipitation over the catchment during a storm for 3 days was as under
1. Light rain 10 mm
2. Average rain 50 mm
3. Heavy rain 150 mm
Find the runoff from the catchment.
Solution:
The runoff is in metre for the classification of the catchment, as per Barlow’s table, the coefficient
being as follows:
Nature
Serial Nature of Runoff of Rainfall Rainfall
no. Day catchment coefficient Rainfall coefficient (m) Runoff (m)
1 1 A1 0.10 Light 0.7 0.01 0.1 0.7 0.01 0.0007
2 2 A1 0.10 Average 1.0 0.05 0.1 1.0 0.05 0.005
3 3 A1 0.10 Heavy 1.5 0.15 0.10 1.5 0.15 0.0225
0.0282 m
Total
4 1 A2 0.35 Light 0.8 0.01 0.35 0.8 0.01 0.0028
5 2 A2 0.35 Average 1.0 0.05 0.35 1.0 0.05 0.0175
6 3 A2 0.35 Heavy 1.7 0.15 0.35 1.7 0.15 0.0893
Total 0.1096 m
7 1 A3 0.45 Light 0.8 0.01 0.45 0.8 0.01 0.0036
8 2 A3 0.45 Average 1.0 0.05 0.45 1.0 0.05 0.0225
9 3 A3 0.45 Heavy 1.8 0.15 0.45 1.8 0.15 0.122
Total 0.1481 m
The runoff in m3 from the three classifications of the catchment will be as under.
Nature of
Serial no. catchment Area (m2) Runoff (m) Runoff (m3)
1 A1 200 106 0.2 0.282 20 106 0.0282 0.564 106
2 A2 200 10 0.4
6 0.1096 40 106 0.1096 4.384 106
3 A3 200 16 0.4
6 0.1481 40 106 0.1481 5.924 106
Total 10.872 106 m3
• Nature of catchment
• Average temperature
The formulae normally followed are as follows:
Example 7.7
A catchment area of 200 km2 has the following classification:
1. Ghat-fed—40%
2. Plain—60%
Find the runoff, if the annual average precipitation is 750 mm.
Solution:
The runoff is calculated by using the Inglis formula.
1. Ghat-fed area:
____
750
R 0.85 P 30.5 0.85 10 30.5 33.25 cm 0.34 m
2. Plain area:
P (P 17.8) __________________
_____________ 750/10(750/1017.8) 75 57.2
_________
R 254 254 254
16.88 cm
0.168 m
Total runoff 200 0.4 106 0.34 200 0.6 106 0.168
27.2 106 20.16 106
47.36 106 m3
RUNOFF 199
Example 7.8
The mean monthly temperature and monthly precipitation figures observed for the year 2004 for a
catchment area, covering an area of 160 km2, are as under.
Month Jan. Feb. Mar. Apr. May Jun. Jul. Aug. Sept. Oct. Nov. Dec.
Precipitation (cm) 5 5 2 0 3 13 33 30 16 2 1 1
Temperature (oC) 13 17 21 28 32 35 32 30 29 27 20 14
Solution:
Using Khosla’s formula, the monthly precipitation, monthly temperature, monthly loss and
monthly runoff, for the year 2004, are tabulated as follows:
____
33.5
Therefore, total annual runoff 160 106 100 53.6106 m3
200 ELEMENTARY ENGINEERING HYDROLOGY
Example 7.9
The classification of a catchment covering an area of 110 km2 is as follows:
1. Rocky, impermeable—10%
2. Cultivated—60%
3. Forest—30%
Using the rational formula and the coefficients suggested by Richard, find the annual runoff
from the catchment when the average annual precipitation is 800 mm.
Solution:
The catchment area, its classification, Richard’s coefficients and the annual runoff is tabulated
below.
Serial no. Catchment area (m2) Richard’s coefficient Annual runoff (m2)
1 0.1 110 106 11 106 0.9 0.9 11 106 0.8 7.92 106
2 0.6 110 106 66 106 0.5 0.5 66 106 0.8 26.40 106
3 0.3 110 106 33 106 0.15 0.15 33 106 0.8 3.96 106
Total 38.28 106
RUNOFF 201
m A constant
p Rainfall in mm
x Interception of the straight-line relationship on the x-axis.
The relation thus derived will be applicable for a specific catchment and is normally not used for
other catchment areas. There are several changes carried out by humans every year and the relation
has to be modified every year to account for them. However, it gives a rough estimate.
Example 7.10
The annual runoff and the annual precipitation observed over a catchment are as follows:
Serial no. 1 2 3 4 5 6 7 8 9 10 11 12
Precipitation (mm) 1000 1700 1200 1500 1360 1620 1100 1370 1630 1420 1400 1260
Runoff (million m3) 553 1105 651 903 750 992 605 781 1005 875 854 751
202 ELEMENTARY ENGINEERING HYDROLOGY
Derive the rainfall–runoff relation for the catchment and also find the possible runoff when the
precipitation will be 1160 mm.
Solution:
The rainfall–runoff for the catchment is plotted on a simple graph paper as shown in Fig. 7.14. From
the graph the interception on the x-axis is 320 mm and the slope is 0.76. Thus the relation works to be
as follow:
R 0.76 (P – 320)
where, R Runoff in million m3
P Annual average precipitation in mm
From this relation, the runoff for a precipitation of 1160 mm will be 638 106 m3.
In the case when the annual runoff figures are observed and are available for n years, these
figures are arranged in a descending order. The serial number m of a specific value of runoff
is known as ranking. It can then be said that this specific value of runoff or more than that is
available for m years out of n years and the percentage dependability of this value will be as
follows:
m
__
% dependability n 100
For instance, if the data of annual runoff are available for 40 years and if these are arranged in
a descending order, then the percentage dependability of the runoff value at Serial no. 30 will be
30/40 100 75%.
Similarly, the percentage dependability of 50% will be the runoff figure at Serial no. 20.
RUNOFF 203
Example 7.11
The annual yield in million m3 from a catchment for the last 20 years is as follows: 240, 350, 680,
400, 290, 160, 110, 105, 300, 600, 800, 540, 280, 130, 320, 810, 740, 220, 190 and 480.
Find the following:
(1) 80% dependable yield
(2) 45% dependable yield
(3) The % dependability of the yield of 540 million m3
(4) Plot the flow–duration curve
Solution:
The annual yield in million m3 for the 20 years is arranged in descending order. The yield for
(1) 80% dependability will be at 200.80 16 (serial order) and (2) 45% dependability will be
at 20 0.45 9 (serial order).
204 ELEMENTARY ENGINEERING HYDROLOGY
All the water after precipitation will not flow to the river. There are some abstractions. After meet-
ing the requirements of these abstractions, the balance water will flow to the river as runoff. Thus
runoff will be less than the precipitation.
The runoff coefficient is the ratio of runoff and precipitation. It will be always be less than 1.
way into the river channel as interflow. Part of water precipitated returns to the atmosphere through
evapo-transpiration.
This process, which is part of the hydrological cycle, is known as runoff cycle.
REVIEW QUESTIONS
1. Define runoff. Discuss the various components of runoff.
2. How is the average slope of a catchment area calculated?
3. Explain watershed leakage, with the help of a neat sketch. How does it affect the runoff from
the catchment?
4. How are the catchment areas classified based on the shape of the catchment? Discuss their
effect on the runoff.
5. Explain bank storage, with the help of a neat sketch.
6. What do you understand by Karst topography?
7. How are the streams classified based on the flow consideration?
8. Discuss the factors affecting runoff.
9. What are the sources of runoff?
10. How is runoff from a catchment estimated?
11. State and explain the different formulae used to estimate runoff from a catchment. Discuss
their limitations also.
12. Explain Binnie’s table to estimate the runoff from a catchment. What are its limitations?
13. What are Strange’s coefficients? Discuss their limitations.
14. What is a dependable yield? How is it calculated for a given dependability?
15. What is a flow-duration curve? How is it constructed?
16. Write short notes on the following:
a. Ridge line b. Valley storage
c. Groundwater divide line d. Antecedent precipitation index
e. Drainage area f. Process of runoff
g. Stream order h. Stream pattern
i. Drainage density j. Rainfall runoff correlation
k. Steam density l. Runoff cycle
m. Length of stream n. Ephemeral stream
o. Form factor p. Compactness coefficient
q. Circulatory ratio r. Relief of catchment area
s. Elongation ratio t. Shape factor
17. Is there any relation between area of the catchment and its length? Discuss.
18. Distinguish between the following:
a. Topographic divide and groundwater b. Shape of the catchment and pattern of the
divide catchment
c. Influent stream and effluent stream d. Intermittent stream and perennial stream
e. Drainage density and stream density f. Form factor and shape factor
g. Circulatory ratio and elongation ratio h. Fan-shaped catchment and elongated
catchment
i. Overland flow and interflow j. Binnie’s table and Barlow’s table for the
estimation of runoff
206 ELEMENTARY ENGINEERING HYDROLOGY
NUMERICAL QUESTIONS
7.1. The levels taken along the cross-section of a channel at a regular interval of 1000 m are as
follows.
Section 0 1 2 3 4 5 6 7
Precipitation in mm 541 581 739 942 838 1022 703 903 781 603
Run-off in mm 263 341 438 705 581 743 442 603 526 408
Derive a relation between precipitation and run-off, for the catchment. Also calculate the
possible run-off when the precipitation is 720 mm
Ans:
(1) The relation between run-off and precipitation is R 0.909 ( P 203 )
(2) The run-off when the precipitation is 720 mm 470 mm
7.4. The annual yield in million m3 from a catchment for the last 20 years is as follows
1410 1310 1100 935 1370 1020 1203 930 1440 1290
980 1400 1176 1425 1340 1261 1130 1300 1480 1280
Calculate
(1) 75 % dependable yield
(2) 50 % dependable yield
(3) the % dependability for 1300 106 m3
Ans:
(1) 75 % dependable yield ______________ 1130 106 m3
(2) 50 % dependable yield ______________ 1290 106 m3
(3) the % dependability for 1300 106 m3 ______________ 45 %
*** ***
Flood
8
Chapter Outline
g 8.1 DEFINITION
Any flow of water that is relatively high and that may overtop the natural or the artificial banks in
a reach of a channel may be called as flood.
Water may overtop the bank at some location where it may be treated that the channel is in
floods. However, at some other location of the same channel, water may be contained within the
banks of the channel. Thus, the same channel may not be considered in floods at that location.
In large rivers, an arbitrary level is decided keeping in mind the habitation or the farmland and
other properties so that when water level rises above this level, the river is designated to be in flood.
210 ELEMENTARY ENGINEERING HYDROLOGY
The factors that may affect the flood intensity are discussed below:
8.2.1 PRECIPITATION
It is obvious that precipitation will affect the flood. More the intensity, more the duration and more
the coverage of the storm, more will be the flood.
and the same storm, the time of concentration will be more. Thus, the storm hydrograph base period
will be more and the peak discharge will be comparatively less.
Figure 8.1 shows the two flood hydrographs in the case of a fan-shaped catchment and an
elongated catchment.
• Flash flood: A flood of short duration and abrupt rise with a relatively high-peak rate of flow,
usually resulting from a high intensity of rainfall, is known as ‘Flash flood’.
• N-year flood: A flood that has a probability of being equalled or exceeded once in N years is
known as an ‘N-year flood’.
Example 8.1
In a village, the maximum water level reached in the last 32 years was recorded on the
compound wall of the Town hall as 105.00 m corresponding to the bed of the river at
100.00 m.
A survey was conducted at this cross section and levels were taken. These are as follows:
Level (m) 106.00 104.00 101.00 100.00 101.00 103.00 105.00 106.0
*
River bed
The river bed slope at this location was observed to be 1 in 2500. The Manning’s coefficient
corresponding to the topography of the river was 0.030.
Find the flood discharge corresponding to the maximum water level reached.
214 ELEMENTARY ENGINEERING HYDROLOGY
Solution:
The cross section of the river is plotted to scale as shown in Fig. 8.3.
The area of flow corresponding to this HFL calculated by a plannimeter was 112.5 m2.
The wetted perimeter corresponding to this HFL was 35.6 m.
_____
112.5
Therefore, hydraulic mean depth R 35.6 3.16 m
Using Manning’s formula
3.162/3 ( 1/2500)1/2
R2/3 S1/2 112.5 __________________
Now, Q A V A __________
n 0.03
162.0 m /s
3
The flood discharge from this curve for the required basin may be calculated. Following equation
for such an envelope curve is derived.
__________
3010 A
Q (8.1)
(277 A)0 78
where, Q Flood discharge in m3/s
A Catchment area in km2
The flood discharge used for plotting the curve are all observed ones and hence the flood
discharge calculated from such an envelope curve has to be increased suitably depending upon
the importance of the structure.
Example 8.2
The catchment area of a proposed bandhara is 100 km2. The maximum average annual flood
discharges in the adjoining area with hydrologically homogeneous catchments are as follows:
Catchment A B C D E F G
2 32.0 120.0 60.0 90 155 80.0 140.0
Catchment area (km )
Maximum flood
1600 3300 2100 2900 3700 2800 3500
discharge (m3/s)
Find the maximum annual flood discharge at the site of the bandhara.
Solution:
On a simple graph paper, catchment area vs discharge, is plotted as shown in Fig. 8.5.
A curve by free hand covering all these points is drawn (envelope curve).
Example 8.3
For a catchment area of 100 km2, find the annual average flood discharge by using the envelope
curve equation.
216 ELEMENTARY ENGINEERING HYDROLOGY
Solution:
The envelope curve equation is as follows:
___________
3010 A
Q
(277 A)0 78
Substituting the value of A,
3010 100
_____________ 3010 100
Q __________
(277 100)0 78 102.22
Q 2944 m3/s.
Example 8.4
Find the maximum flood discharge by using the empirical formulae for a catchment area of 80 km2
having a maximum precipitation intensity of 4 cm/h.
Assume
C in Dicken’s formula = 22
C in Ali Nawaz Jung Bahadur formula = 50
C in rational formula = 0.35
Solution:
The flood discharge Q by using the various empirical formulae will be as follows:
1. Dicken’s formula
Q CA3/4 22 803/4 22 26.74
588 m3/s
2. Inglis formula
124 A
___________ 124 80
____________ 124 80
________
Q _________ __________ 9.51
√ (A 10. 4) √ (80 10.4)
1043 m3/s
3. Ali Nawaz Jung Bahadur formula
___
1
0 925 14 log (0 3906 A)
Q C(0.39006 A) 10
___
1
0 925 14 log10 (0 3906 80)
50 (0.39006 80)
50 31.20 819
837 m3/s
4. Rational formula
Q 2.778 CAI 2.778 0.35 80 4
311 m3/s
8.4.4 UNIT HYDROGRAPH PROCEDURE
The unit hydrograph theory and procedure can be used to estimate the flood. For this purpose, the
following data are required:
1. Possible maximum precipitation: The possible maximum precipitation (PMP) that can be
expected over the basin can be worked out from the available data. India Meteorological
Department (IMD) has divided India in various zones for this purpose, PMP in each zone has
been indicated. This data can also be used.
218 ELEMENTARY ENGINEERING HYDROLOGY
2. Unit hydrograph: The unit hydrograph derived for the catchment area up to the basin outlet is
necessary. If such a unit hydrograph is not available, then synthetic unit hydrograph up to the
basin outlet may be derived.
By using the PMP and the unit hydrograph, the flood hydrograph up to the basin outlet can be
estimated.
The peak discharge in the flood hydrograph thus derived may be considered as the design
flood. With the use of this unit hydrograph procedure, not only the flood discharge is estimated
but the flood hydrograph thus derived can also be used for the design of a spillway of a dam for
flood routing purposes.
Example 8.5
For a catchment area of 50 km2, the 1 h–1 cm unit hydrograph coordinates are as follows:
Ti (h) m e 0 1 2 3 4 5 6 7 8 9 10
Discharge
0 6 28 35 27 20 14 9 5 2 0
(m3/s)
The worst possible storm covering the entire catchment area is 2 cm/h for 3 h. Find the maximum
flood discharge.
Solution:
Since the unit hydrograph is for 1 h, the storm will be divided into three parts each one of 1-h duration.
Discharge due to these three parts of the storm with a time lag of 1 h will be as follows:
Unit
hydrograph Discharge
Serial no. Time (h) coordinates due Total
First part Second part Third part
1 0 0 0 — — 0
2 1 6 12 0 — 12
3 2 28 56 12 0 68
4 3 35 70 56 12 138
5 4 27 54 70 56 180
6 5 20 40 54 70 164
7 6 14 28 40 54 122
8 7 9 18 28 40 86
9 8 5 10 18 28 56
10 9 2 4 10 18 32
11 10 0 0 4 10 14
12 11 — — 0 4 4
13 12 — — — — 0
2. Hazen’s formula,
________
n
Tr m 0.50
3. Weibull’s formula,
n1
_____
Tr m
4. Gumbell’s formula,
__________
n
Tr m C 1
Here, C is called gumbell’s correction. Its value depends on m and n as given in the following
table.
m/n 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.04
C 1 0.95 0.88 0.845 0.78 0.73 0.68 0.59 0.52 0.4 0.28
Once recurrence interval Tr is calculated, the probability p can be calculated, since Tr = 1/p.
Weibul’s formula is generally used.
Example 8.6
The maximum annual flood discharges in m3/s observed at a proposed dam site from 1973 to 2004
(32 years) are as follows:
395, 619, 766, 422, 282, 887, 705, 528, 520, 436, 697, 624, 496, 489, 598, 359, 696, 726, 527, 310,
408, 721, 814, 459, 440, 632, 343, 634, 464, 373, 289, 371
Find the return period for all the discharges by using the following formulae:
(1) California formula, (2) Hazen’s formula, (3) Weibul’s formula and (4) Gumbell’s formula.
Solution:
The four formulae referred to are as follows.
__
n
1. California formula, Tr m
________
n
2. Hazen’s formula, Tr m 0.50
n1
_____
3. Weibul’s formula, Tr m
__________
n
4. Gumbell’s formula, Tr m C 1
The series of 32 discharges was arranged in the descending order. The return period
calculated by the first three formulae are tabulated in Table A. Since the procedure for the
fourth method is slightly different, the return period worked out by this method is tabulated
in Table B.
FLOOD 221
Table A
Discharge in
descending Tr
Serial no. Discharge(m3/s) order (m3/s)
California Hazen’s Weibull’s
formula formula formula
1 395 887 32 64 33
2 619 814 16 21.33 16.5
3 766 766 10.67 12.80 11
4 422 726 8 9.14 8.25
5 282 721 6.40 7.11 6.6
6 887 705 5.33 5.81 5.5
7 528 697 4.57 4.92 4.71
8 705 686 4 4.27 4.12
9 520 634 3.55 3.76 3.67
10 436 632 3.2 3.36 3.3
11 697 624 2.9 3.04 3
12 624 619 2.67 2.78 2.75
13 496 598 2.46 2.56 2.54
14 589 589 2.28 2.37 2.36
15 598 528 2.13 2.2 2.2
16 359 527 2 2.06 2.06
17 686 520 1.88 1.94 1.94
18 726 496 1.77 1.82 1.83
19 527 464 1.68 1.73 1.74
20 310 459 1.6 1.64 1.65
21 408 440 1.52 1.56 1.57
22 721 436 1.45 1.49 1.5
23 814 422 1.39 1.42 1.43
24 459 408 1.33 1.36 1.38
25 440 395 1.28 1.31 1.32
26 632 373 1.23 1.25 1.27
27 343 371 1.18 1.21 1.22
28 634 359 1.14 1.16 1.18
29 464 343 1.1 1.12 1.13
30 373 310 1.06 1.08 1.1
31 289 289 1.03 1.05 1.06
32 371 282 1 1.02 1.03
222 ELEMENTARY ENGINEERING HYDROLOGY
Table B
Discharge in
Discharge descending
m__ Tr
Serial no. (m3/s) order (m3/s) n C mC1
1 395 887 0.030 0.27 0.27 118.5
2 619 814 0.0620 0.29 1.29 24.8
3 766 766 0.09 0.38 2.38 13.44
4 422 726 0.125 0.42 3.42 9.35
5 282 721 0.156 0.46 4.46 7.17
6 887 705 0.187 0.51 5.51 5.8
7 528 697 0.218 0.52 6.52 4.9
8 705 686 0.25 0.53 7.53 4.25
9 520 634 0.289 0.595 8.55 3.74
10 436 632 0.312 0.58 9.58 3.34
11 697 624 0.344 0.59 10.59 3.02
12 624 619 0.375 0.64 11.64 2.75
13 496 598 0.406 0.65 12.65 2.53
14 589 589 0.437 0.69 13.69 2.33
15 598 528 0.469 0.71 14.71 2.17
16 359 527 0.5 0.73 15.73 2.03
17 686 520 0.531 0.75 16.75 1.91
18 726 496 0.562 0.77 17.77 1.86
19 527 464 0.593 0.78 18.78 1.71
20 310 459 0.625 0.79 19.79 1.62
21 408 440 0.656 0.81 20.81 1.54
22 721 436 0.689 0.84 21.84 1.46
23 814 422 0.718 0.85 22.85 1.4
24 459 408 0.75 0.86 23.75 1.35
25 440 395 0.781 0.87 24.87 1.29
26 632 373 0.812 0.89 25.89 1.24
27 343 371 0.843 0.91 26.91 1.19
28 634 359 0.875 0.93 27.93 1.14
29 464 343 0.906 0.95 28.95 1.11
30 373 310 0.937 0.967 29.97 1.07
31 289 289 0.969 0.99 30.99 1.03
32 377 282 1.0 1.0 31.00 1.02
All these formulae will give the recurrence interval and probability for any value covered in the
series. But if the recurrence interval and the probability of any value more than the values in the
series are required, then this can be obtained by extrapolation by using the following graphs:
FLOOD 223
1. Natural scale
2. Semi-log
3. Log-log
The results obtained by using log-log plot are considered more accurate.
Example 8.7
For the data in Example 8.6, find the discharge with a return period of 50.
Solution:
The return periods for all the 32 discharges have been worked out by four different formulae in
Example 8.6. These results will be extrapolated on graph paper. The following three types of graph
papers will be used.
1. Natural scale graph paper (Fig. 8.6)
2. Semi-log scale graph paper (Fig. 8.7)
3. Log-log scale graph paper (Fig. 8.8)
Fig. 8.6 A simple graph of the discharge for the return period of 50
Fig. 8.8 A log-log graph of the discharge for the return period of 50
The discharges calculated from the three graph papers for a return period of 50 years are
tabulated below.
pr (1 p)n
FLOOD 225
Example 8.8
For a dam, the designed discharge of 1000 m3/s has a return period of 50 years. Find the
probability that
(a) This discharge occurs twice in 20 years.
(b) This discharge occurs once in 15 years.
(c) This discharge will not occur at all in 25 years.
Solution:
1
p __
50 0.02
Note: The procedure discussed in this sub-chapter is applicable to all hydrological parameters,
e.g. precipitation, runoff, etc.
For the estimation of floods of higher recurrence intervals and probabilities, Gumbell’s distribution
is widely used. Gumbell considered the annual flood series.
Let Q1, Q2, Q3, … Qn be the discharge series of n figures arranged in descending order. As per
this distribution, the probability of occurrence of a discharge figure is given by
y
p 1 e e (8.6)
where, p Probability,
y Reduced variate that is a dimensionless number
The value of y a(Q Qf ) (8.7)
Here, a and Qf are the parameters of the distribution.
a n/sx
__
and Qf Q (yn / n) sx
__
where, Q and sx are the mean and standard deviation, respectively, of the discharge series and
n and yn depend on n, i.e. the number of discharge series. These values are given in following
table.
226 ELEMENTARY ENGINEERING HYDROLOGY
n n yn n n yn n n yn
Yn reduced mean
n reduced standard deviation
y Tr 1
______
Therefore, e e Tr
Taking logarithms of both sides,
(Tr 1)
_______
ey loge Tr
(Tr 1)
_______ (Tr 1)
_______ Tr
_______
ey loge 2.303 log10 2.303 log10
Tr Tr (Tr 1)
Again, taking logarithms of both sides,
Tr
(______
y loge 2.303 loge log10 T 1
r
)
⬖ [ Tr
______
y loge 2.303 loge log10 T 1
r
( )]
Therefore, [ Tr
______
y 0.834 2.303 log10 log10 T 1
r
( )] (8.8)
yT a (QT Qf)
__
n yn
___
[ ___
{
( sx ) QT Q n sx ( ) }]
__ yn
___
sx ___ ( )
yT ( n ) QT Q n sx
__ yn
___
sx ___( )
Therefore, QT Q yT ( n ) n sx
__ y yn
QT Q sx [ ( ___
n ) ( n ) ]
T ___
__ y yn
Q ( _______
T
n ) sx (8.9)
__
Q KT sx.
KT is known as frequency factor (yT yn)/n.
The step-by-step procedure to follow Gumbell’s distribution is as follows:
1. Arrange the discharge series for
__ n years in the descending order.
2. Find out the mean discharge Q and standard deviation sx of the discharge series.
3. From Table 8.2, find the values of yn and n corresponding to n.
4. From Eq. 8.8, find the value of yT for the desired value of Tr.
5. From Eq. 8.9. find the value of QT .
Thus, QT will be the discharge having a return period of Tr .
Gumbell’s distribution procedure can also be used by following a graphical solution. A special
graph Gumbell probability paper, generally known as Powell’s probability paper, can be used. This
graph paper is shown in Fig. 8.9.
On the x-axis, the return period Tr is plotted and on the y-axis, the discharge is plotted. This
graph is a straight line. From this graph, the discharge for any return period can be calculated.
This type of distribution is used for anticipating flood of higher return period. It is commonly used
in the USA.
Let the observed flood discharge series be Q1, Q2, Q3, …, Qn of n figures, arranged in the
descending order. This discharge series is converted into a new log series, viz. log Q1, log Q2,
log Q3, …. Of course, the number of terms, i.e. n will remain unchanged. The base for the loga-
rithms should be 10.
This new log series will now be considered further. The value of the variate yT for a return
period Tr is given as follows:
_
yT = y K sx (8.10)
where, yT = Reduced variate
_
y = Mean of the log series
K = Frequency factor
sx = Standard deviation of the log series
The coefficient of skew g will be as follows:
_
n∑ (y y)3
_______________
g=
(n 1) (n 2)sx3
The frequency factor K is a function of Tr and g.
The following table gives the values of K for various values of Tr and g.
From the value of g calculated for the log series and Tr, the value of yT can be calculated from
Eq. 8.10.
From yT , the value of QT can be calculated as follows:
QT = Antilog ( yT) (8.11)
The step-by-step procedure to follow the Log-Pearson type III distribution is as follows:
1. Arrange the flood discharge series of n terms in descending order.
2. Convert this series into logarithmic series,
_
taking the base as 10.
3. For this log series, find out the mean, y ; standard deviation, sx and the skew coefficient, g.
4. For the required return period Tr and g, find out the value of K from Table 8.3.
5. Find yT from Eq. 8.10.
6. Obtain the value of QT from equation 8.11.
QT is the discharge with the required return period Tr.
Example 8.9
For the data given in Example 8.6, find the discharge with a return period of 50 years using
1. Gumbell’s distribution 2. Log-Pearson Type III distribution.
Solution:
1. Gumbell’s distribution
From the data, the following were computed:
(a) Number of terms in the series n
__ 32
(b) Mean discharge Q 535 m3/s
(c) Standard deviation sx 162.93
(d) From Table 8.2, for n 32 n 1.1193
and yn 0.5380
Now, { ( Tr
_____
yT 0.834 2.303 log10 log10 T 1
r
)}
Substituting the value of Tr 50,
We get,
{ ______
(50
yT 0.834 2.303 log10 log10 50 1 )}
3.9028
(yT yn)
__
Now, QT Q+ sx {
________
n }
Substituting the values of sx, yT , yn and n,
(3.9028 0.5380)
We get, QT 535162.93 { ________________
1.1193 } 535 490
1025 m3/s.
230 ELEMENTARY ENGINEERING HYDROLOGY
Table C
3
Serial no. Q (m /s) Q (m3/s) log10Q (m3/s)
1 395 887 2.9479
2 619 814 2.9106
3 766 766 2.8842
4 422 726 2.8609
5 282 721 2.8579
6 887 705 2.8481
7 528 697 2.8432
8 785 686 2.8363
9 520 634 2.8020
10 436 632 2.8007
11 697 624 2.7951
12 624 619 2.7916
13 496 598 2.7767
14 589 589 2.7701
15 598 528 2.7226
16 359 527 2.7218
17 686 520 2.7160
18 726 496 2.6954
19 527 464 2.6665
20 310 459 2.6618
21 408 440 2.6434
22 721 438 2.6394
23 814 422 2.6253
24 459 408 2.6106
25 440 395 2.5965
26 632 373 2.5717
27 343 371 2.5693
28 634 359 2.5550
29 464 343 2.5352
30 373 310 2.4913
31 289 289 2.4608
32 377 282 2.4502
0.195
For Tr 50 and g 0.195, K 1.946
Now, _
yT y+ K sx 2.7080 + (1.946 0.1377) 2.9759
QT Antilog 2.9759 946 m3/s.
g 8.7 RISK
Any hydraulic structure is designed for a discharge having a probability of p. So also this hydraulic
structure is expected to have a useful life of N years. There is a probability that this discharge or
more than that may occur in the lifetime of the structure. Thus, there is some sort of risk involved
in this.
This risk R is given as follows:
R 1 (1 p)N
__ N
or,
1
R1 1T ( r
)
On percentage basis, R [1 (1 p)N] 100
or, [ ( __
1
R 1 1T
r
) ] 100
N
Naturally, in practice the acceptable risk is decided by the economic and policy considerations.
1 N
__
And reliability will be 1 R (1 p)N 1 T
r
( )
Example 8.10
A cofferdam is designed for a return period of 30 years. It will require 3 years to complete the dam.
Find 1. The risk that the cofferdam may get washed away.
2. If the risk is expected to be only 5%, what will be the return period for which the
cofferdam should be designed?
Solution:
__ N 3
( 1
1. R 1 1 T
r
) 1 ( 1 ___301 ) 8.74%.
232 ELEMENTARY ENGINEERING HYDROLOGY
__ 3
( 1
2. 0.05 1 1 T
r
)
Therefore, Tr 59 years
REVIEW QUESTIONS
1. What is flood? Discuss the various features of flood.
2. Discuss the various causes of floods.
3. What are the factors affecting flood?
4. How floods are classified?
5. When there are no observed data available, how the flood is estimated at a river outlet?
6. What is an envelop curve? How is it used to estimate for estimating flood at a site?
7. Which are the different empirical formulae used to estimate the flood discharge?
Write a note on the rational formula used to estimate the flood at a site.
9. Explain how the unit hydrograph theory is used to estimate flood discharge from a
catchment area.
10. Differentiate between annual flood series and partial duration flood series.
11. Explain:
a. Recurrence interval
b. Probability
c. Frequency of floods
12. Explain the procedure to evaluate the probability of a specific discharge in an observed
discharge series. Discuss the various formulae used to asses these.
13. The probability of a specific discharge, which is more than the observed discharge series,
is to be calculated. Discuss the various methods.
14. Describe the method of estimating a Tr-year flood using Gumbell’s distribution.
15. Describe the procedure of estimating a Tr-year flood using Log-Pearson type
distribution.
16. How is the risk in the design that a structure may fail in its lifetime calculated?
17. Write short notes on:
a. Frequency factor
b. Rank of observation
c. N-year flood
d. Effect of the shape of catchment on the flood
e. Gumbell probability paper
NUMERICAL QUESTIONS
1. In a village, the maximum water level reached on the wall of the school temple is 105 m. The
cross section of the river is in the form of a symmetrical triangle with vertical angle 90 and
apex at a level of 100 m.
Assuming the slope of the river to be 1 in 3600 and the Manning’s coefficient to be 0.035, find
the discharge corresponding to this maximum water reached.
Ans: 17.38 m3/s
FLOOD 233
2. The catchment area of a proposed weir is 80 km2. The maximum flood discharges in the
adjoining area with hydrological homogeneous catchments are as follows:
Catchment A B C D E F
Catchment area (km2) 103 176 48 141 72 220
3. Using the empirical formulae, find the maximum flood discharge for a catchment of 100 km2,
having the maximum intensity of precipitation of 4.0 cm/s.
Ans: Inglis formula 1180 m3/s; Dicken’s formula 790 m3/s;
Rational formula 388 m3/s; Ali Nawaz Jung Bahadur formula 979 m3/s
4. For a weir, the design discharge of 500 m3/s has a return period of 30 years.
Find: (i) This discharge occurs twice in 25 years
(ii) This discharge occurs once in 25 years
(iii) This discharge will not occur at all in 20 years
Ans: (i) 12.7%, (ii) 56%, (iii) 51%
5. The maximum flood discharges in m3/s observed at a proposed dam site for 1985–2004
(20 years) are as follows:
825, 480, 713, 671, 451, 870, 462, 648, 601, 425, 890, 620, 540, 580, 705, 465, 790, 680,
733, 540
Find the discharge with a return period of 25 years using
(i) Gumbell’s distribution
(ii) Log-Pearson type III distribution
Ans: (i) 1020 m3/s, (ii) 931 m3/s
6. A cofferdam is a design period of 25 years. It will require four rears to complete the dam.
Find: (i) The risk that the cofferdam may get washed away
(ii) If the risk is expected to be only 5%, the return period for which the cofferdam
should be designed
Ans: (i) 15% (ii) Tr 76 years
8. A structure having a life of N years is designed for a return period of Tr years. The probability
that it will not fail during its life period is
__
1
(a) 1 __1 (b) 1 T
Tr r
__
1 N
(c) (1 ) __
1 N
(d) (1 T )
Tr r
9. Dicken’s formula for calculating the maximum flood discharge from a catchment is given by
(a) Q CA4/3 (b) Q CA3/4
(c) Q CA2/3 (d) Q CA3/2
10. Inglis formula for calculating the maximum flood discharge from a catchment is given by
124 (A __ 10.4)
_____________ 124 (A __ 10.4)
_____________
(a) Q (b) Q
√A √A
___________
124 A ___________
124 A
(c) Q _________ (d) Q _________
√ (A 10.4) √ (A 10.4)
FLOOD 235
11. As per Gumbell’s distribution, the probability of occurrence of a discharge figure is given by
y y
(a) p = 1 e e (b) p = 1 e e
y y
(c) p = 1 e e (d) p = 1 e e
12. For a hydraulic structure designed for a return period of Tr years and having a useful life of
N years, the risk that this flood may occur during the useful life is
[ ( )]
(a) R = 1 1 T
__
1 N
r
100 [ ( )]
(b) R = 1 1 T
__
1 N
r
100
r r
*** ***
Discharge
Measurement
9
Chapter Outline
g 9.1 DEFINITION
Discharge is defined as the rate of flow, i.e. volume flowing per unit time.
following some empirical formulae, but this will not be accurate. Even for the derivation of an
empirical formula, observed discharge data are required.
The surface water resource can be evaluated by measuring discharge in a stream continuously
over a period. The surface water resource available at a site can be evaluated by scanning the
observed data. This will give more realistic estimate and may involve less error.
(11,000)
________
3600 (2510 ) 6.9 m /s
6 3
The depth of water flowing in a river is known as The Stage of the River at that time. The stage
of the channel is measured because from this observation, the area of flow can be calculated. The
stage of the river is expressed in metres with datum as mean sea level. However, sometimes it is
also mentioned with datum as lowest bed level.
Inclined gauges
Gauges along the slope of the cross section of the stream may be fixed as shown in Fig. 9.1. While
marking the depth of flow on the gauge, the side slope of the stream has to be taken into consider-
ation. Such gauges will not cause any obstruction to the flow.
In both the cases, the water level will have to be observed from the banks of the stream. Know-
ing the water level and the bed level of the stream, the depth of flow can be evaluated.
By following this procedure, the stage of the river can be found out at that moment. However, the
maximum and the minimum stages reached cannot be found out unless observed at those instances.
a week as the case may be. The maximum and minimum water levels reached and the time when reached
are automatically recorded.
In recent years, the water levels thus recorded are transferred by a wireless arrangement to a central
recording station, where all meteorological data in that region are collected.
Crest gauge
When it is not possible to install an automatic water level recorder due to any reason, a crest gauge
may be installed in the observation well. This crest gauge will record the maximum water level
reached. This can be done by two ways:
1. It consists of a small float that rises as the water level rises in the well but will not move down as the
water level in the well reduces. Thus, the maximum water level reached in the well is recorded.
2. One of the sides of the well or a bridge pier may be painted by a water-soluble paint, of course,
well protected from rain. As the water level rises, the portion of the paint coming in contact with
the water will get washed away and thus the maximum water level reached can be obtained.
Example 9.1
After doing precise survey, the cross section of the stream at a gauging site is as shown in Fig. 9.5.
Prepare the two graphs: (1) depth vs area and (2) depth vs hydraulic mean depth.
Solution:
The cross section was divided into five segments each of 1-m height and the area of each segment
was measured accurately. So also the length of the cross section coming in contact with the water
on either side was measured and tabulated.
Serial no. Depth Segment Area (m2) ∑ Area PL (m) PR(m) ∑ (PL PR) (m) HMD(m)
(m) (m2)
1 1.0 A1 5.0 5.0 — 9.9 9.9 0.505
2 2.0 A2 9.0 14.0 2.2 2.0 14.1 0.992
3 3.0 A3 13.5 27.5 1.8 1.7 17.6 1.560
4 4.0 A4 20.5 48.0 1.9 1.7 21.2 2.264
5 5.0 A5 24.0 72.0 1.8 1.6 24.6 2.926
Area
_________
Here, HMD
(PLPR)
Two graphs are drawn as (1) Depth vs Hydraulic Mean Deapth (HMD) and (2) Depth vs area,
on simple graph papers as shown in Fig. 9.6.
These pulses are reflected back by the bed of the channel. The time of dispatch and receipt of pulses
is recorded. The depth of flow then can be worked out from this time interval and the speed of the
sonic pulse in water. Hence, the name of the equipment is ‘Echo sounder’.
Such observations are taken at different locations on the cross section preferably at equal dis-
tance and the cross section of the stream is plotted.
Alternately, the boat may be moved along the cross section at right angles to the flow at a
constant speed. For this, two theodolites on either bank may be installed to guide the boat properly,
or any other ranging equipment or procedure may be adopted. Observations of depth by the echo
sounder may be continuously taken and plotted on a graph paper. With a continuous plotter, the
cross section of the stream at a specific water level may be obtained directly on the graph paper.
Echo sounder can also be used for the study of silting of reservoirs.
The basic principle followed in this method is that the area of flow is measured and the velocity
of flow is evaluated by using the Manning’s formula. For this purpose, the water surface slope is
observed and the Manning’s ‘n’ is assumed. Then the continuity equation, Q A V is used.
where, Q Discharge in m3/s
A Area of cross section in m2 at right angles to the velocity of flow
V Velocity of flow in m/s
The procedure to be followed is as follows:
1. Select two cross sections of the stream at a distance of about 150 m. For the selection of this
reach, the criteria as explained in Section 9.3.1 are followed.
It is assumed that steady uniform flow exists between these two cross sections.
2. Observe the depths of flow at these two cross sections as well as the water surface level accu-
rately between these two cross sections.
3. Find the area of flow and also the wetted perimeter at these two cross sections from the pre-
determined stage vs area and stage vs wetted perimeter curves.
DISCHARGE MEASUREMENT 245
4. Find the mean area of flow Am and the mean wetted perimeter Wp.
5. From the water levels at the two cross sections, find accurately S, i.e. the slope of the water
surface. (This will be the slope of energy also, since ‘steady uniform flow’ is assumed between
these two cross sections.)
6. Calculate the mean hydraulic radius Rm as
Am
Rm ___
Wp
7. Assume a suitable value of the Manning’s n for the reach of the stream between these two cross
sections.
8. Calculate the average velocity of flow using the Manning’s formula as
(Rm S )
2/3
___________
1/2
Vm n
where, Vm Average velocity of flow in m/s
Rm Mean hydraulic radius thus calculated in m
S Energy slope
n Assumed value of Manning’s coefficient
9. Calculate the discharge as
Q Vm Am
The discharge is calculated by using the area and the slope of the water surface and hence it is
known as Area–slope Method. To achieve accuracy, the procedure may be repeated and the average
value of discharge of all the trials may be finalized.
Example 9.2
While taking discharge observations on a stream by area–slope method, the following observations
were taken:
1. Distance between two observation sections of the stream 250 m
2. Depth of flow at the u/s cross section A 3.50 m
246 ELEMENTARY ENGINEERING HYDROLOGY
Solution:
From the graphs of depth vs area of cross section and depth vs wetted perimeter, the area of cross
section and the wetted perimeter at both A and B were evaluated as follows:
110 108
_________
Therefore, Mean area 2 109 m2
76.2 75.8
__________
Mean wetted perimeter 2 76.0 m
____
109
HMD Rm 76 1.434 m
In this method also the basic continuity equation, Q A V, is used to estimate the discharge
flowing in a stream. The values of A and V are assessed separately and then the discharge is
calculated.
( V0 2D V0 8D )
____________
2
( V0 2D 2V0 6D V0 8D )
____________________
4
*
0.6D from the water surface, i.e. 0.4D from the bed of the channel.
Example 9.3
Observations for velocity of flow taken over a vertical of a stream at various depths are as follows.
Depth from water surface (m) 0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
Velocity of flow (m/s) 1.47 1.49 1.52 1.48 1.41 1.34 1.24 1.14 1.00 0.82 0.00
Find the average velocity. Plot the velocity profile, and derive an equation for the velocity profile.
248 ELEMENTARY ENGINEERING HYDROLOGY
Solution:
The velocity profile will be as shown in Fig. 9.8.
1.47 1.49 1.52 1.48 1.41 1.34 1.24 1.14 1.00 0.82 0.00
_______________________________________________________________
11
_____
12.91
11 1.17 m/s
V0 2D 2 V0 6D V0 8D )
(______________________ 1.52(21.24)1.00 ____
___________________ 5.00
V m
4 4 4 1.25 m/s
Vm V0 6D 1.24 m/s
DISCHARGE MEASUREMENT 249
logV
____
or, n
log d
The value of n for three observations will be as follows.
• As the float crosses the cross sections P-P and Q-Q, the observers at these cross sections should
stop their respective stopwatches. Thus, the time of travel from the first cross section P-P to the
next Q-Q, i.e. to travel a distance L can be calculated.
Let the time taken by the float to travel each of the segment be T1, T2, T3,…. Then the surface
velocity in each of the segment will be
L
__ L
__ L
__
Vs1 T ; Vs2 T ; Vs3 T
1 2 3
It may be noted that T1, T2, T3,… will not be equal, since the surface velocity over a cross section
is maximum at the centre and goes on decreasing towards the banks.
Normally, the average velocity over a depth is less than the surface velocity, i.e.
Va 0.85 Vs.
DISCHARGE MEASUREMENT 251
• The average velocity in each segment, Va1, Va2, Va3,… is calculated and then the discharge in the
stream is calculated as:
Q Q1 Q2 Q3 …
Example 9.4
Observations for discharge measurement taken on a stream by float method are as follows:
1. Number of segments 5
2. Depth of flow observed at the centre of each segments (m) 0.75, 1.20, 1.75, 1.35, 0.60
3. Width of the stream at the water surface 19.70 m
4. Width of each segment (m) 3.1, 4.5, 5.0, 4.3, 2.8
5. Distance between the two cross sections 100 m
6. Time taken by the surface floats in each segment is as follows:
Segment Surface velocity (m/s) Average velocity Area of flow (m2) Discharge
(m/s) (m3/s)
1 100/131 0.763 0.649 0.753.1 2.325 1.508
Total 18.707
DISCHARGE MEASUREMENT 253
(ii) The velocity observations may be taken at 0.2D and 0.8D also. The average velocity may
be worked out as
( V0 2D V0 8D )
____________
Va 2
This is known as Two-point Method.
(iii) The velocity observation may be taken at 0.2D, 0.6D and 0.8D from the water surface and
the average velocity can be worked out as
( V0 2D 2V0 6D V0 8D )
____________________
Va 4
This is known as Three-point Method and normally followed for deep flow streams.
(iv) Observations may be taken by lowering the current meter in each of the segment with
a uniform speed from the water surface till it reaches the bottom and lifted gradually
upwards with the same uniform speed up to the water surface.
The total time for the movement of the current meter as well as its speed is recorded. The
average velocity in each segment is calculated as follows:
N
__
Va ab T ( )
where, Va Average velocity in each segment in m/s
N Number of rotations in time T required for downward and upward movement of
the current meter
a, b Calibration constants
This method is known as Integrated Method of Velocity Observation.
Air-line correction
• The air-line correction is the correction for the portion of the chain from the point of suspension
up to the water surface.
• This correction depends on the vertical angle of drift and the vertical length from the point of
suspension to the water surface.
The correction will be Ka ( sec 1) 100
Thus, after applying the air-line correction, the length of the chain from the suspension point up to the
water surface level will be
d1 Y1 (10.01 Ka)
DISCHARGE MEASUREMENT 257
Wet-line correction
• The chain under the water surface will not be having a straight profile but a curved one, depending
upon the velocity of water, unit weight of water, unit weight of the chain material and also the
weight of current meter and fish weight.
• The vertical depth of water from the water surface up to the current meter will be as follows:
d2 Y2 (1 0.01 Kw)
where, Y2 Vertical depth of water from the water surface up to the current meter
d2 Inclined length of the cable from the water surface up to the current meter
Kw Wet-line correction. It is dimensionless.
The two correction factors, viz Ka and Kw will mainly depend on angle . The standard values in
percentages are given in Table 9.1.
Table 9.1 Standard values of K a and K w for different values of
Serial no. (degree) Ka (%) Kw (%) Remarks
1 4 0.24 0.06 —
2 8 0.98 0.32 —
3 12 2.23 0.72 —
4 16 4.03 1.28 —
5 20 6.42 2.04 —
6 24 9.46 2.96 —
7 28 13.26 4.08 —
Example 9.5
Observations are to be taken for the velocity of flow by lowering a current meter from bridge with
the following details:
1. Depth of flow 5.00 m
2. Distance between the index point and the water surface 4.00 m
3. Angle made by the cable with the vertical 24
Find the length of the cable to reach the bed of the channel.
Solution:
9.46*
_____
Air-line correction 4 100 0.378 m
2.96*
_____
Wet-line correction 5 100 0.148 m
*
Taken from the Table 9.1 for an angle of 24
Therefore, length of the cable from the index point up to the bed of the channel
4 .378 5 .148
9.526 m.
Example 9.6
Observations taken on a stream for the measurement of discharge by a current meter are as
follows:
1. Width of the stream at the water surface 19.70 m
2. Number of segments 5
3. Width of each segment (m) 3.1, 4.5, 5.0, 4.3, 2.8
4. Average depth of flow of each segment (m) 0.75, 1.20, 1.75, 1.35, 0.60
The observed velocity of flow at various depths from the water surface are as shown in Fig. 9.15.
Find the discharge.
Solution:
The calculations done for the discharge are as under.
V0 2D V0 8D V__________________
___________ 2V0 6D V0 8D
Here, VA V0 6D
0 2D
2 4
2. Based on the velocity observed at various depths in each segment isovels at equal interval are
drawn on the cross section as shown in Fig. 9.16.
An Isovel May be Defined as a Contour of Equal Velocity, i.e., a Line Joining Locations with
Same Flow Velocity.
The area between two consecutive isovels is measured accurately and tabulated. The discharge in
the stream is calculated as
Q Q 1 Q2 Q 3 …
A V A V A V …
1 a1 2 a2 3 a3
Example 9.7
Velocity observations taken over a cross section of a stream at various locations are as shown in
Fig. 9.17. The width of the cross section at the surface and the maximum depth are 23.67 m and
5.5 m, respectively.
Draw the isovels and evaluate the discharge.
Solution:
The isovels were drawn taking into consideration the velocity and the location. The area between
the consecutive isovels was measured accurately. The results are tabulated below.
Serial no. Isovel range (m) Average of isovel (m) Area (m2) Discharge (m3/s)
1 1.2 and above 1.20 1.41 1.692
2 1.2 and 1.1 1.15 6.21 7.142
3 1.1 and 1.0 1.05 8.75 9.188
4 1.0 and 0.9 0.95 5.45 5.177
5 0.9 and below 0.90 1.85 1.665
Total 23.67 24.864
It is normally connected to a differential U-tube mercury manometer and the velocity head is
measured in terms of the depth of the fluid. The velocity of flow at depth d can be equated as follows:
_____
v k 冪(2gh)
where, v Velocity of flow in m/s
h Rise of water column in the Pitot tube in m
k Coefficient of Pitot tube to account for the loss of head. (This is dimensionless and
always less than one.)
n1 q n0 Q n2 (Q q)
it is called Plateau Method. When it is added all at once, it is called Gulp Method. In the case of
gulp method, chemical is mixed all at once. Then the observations of the normality of the stream
water after thorough mixing should be taken at a regular interval, say t, till no trace of chemical
is observed in the stream water.
In this case, the discharge in the stream Q can be calculated by the equation:
n1 V
Q
(∑n2 n0) t
where, V Volume of chemical added
The method can be used for small streams.
It is simple because it does not involve measurements of depth, velocity, etc.
Example 9.8
20 mg/cc of salt at the rate of 10 cc/s was discharged in a stream having a concentration of the same
salt of 0.1 parts per billion (ppb). The concentration after thorough mixing was found to be 5 ppb.
Find discharge in the stream.
Solution:
q 10 cc/s 1105 m3/s
____
20 ___
1
Concentration of salt added in the stream n1 0.02
106 109
___
0.1
Concentration of salt in the flowing water n0 0.1109
109
___
5
Concentration after the thorough mixing n2 9 5 109
10
q(n1n2) ____________________
5 9
________ 110 (0.02510 )
(n2n0) 9
510 0.110 9
40.81 m3/s
(1 mg/1 m3 1 part 1 per billion)
Example 9.9
30 litres of water with a salt concentration of 20 mg/cc was added to a flowing stream in gulp.
Concentration of the stream water after thorough mixing was measured on the d/s every after 3 min
and was noticed to be 0, 5, 10, 20, 15, 5, 0 ppb.
Find the discharge in the stream.
DISCHARGE MEASUREMENT 263
Solution:
Since the normality of the stream water n0 before mixing is not given, it is assumed to be zero.
____
20 30 103
_________
6
109
_____________________________________________________
10
Q
3 60 [(5 0) (10 0) (20 0) (15 0) (5 0)]
____________________________________________________
109
600 103
600 106 ___________
____________ 3 60 55 60.60 m /s
3
3 60 55
___________
109
9.7.1 NOTCHES
Notches are used when the discharge is small. The different notches are as follows:
1. Triangular notch
2. Rectangular notch
3. Trapezoidal notch
4. Cippolette notch
Example 9.10
In a laboratory, following 4 notches were installed in parallel with the crest of all the notches at the
same level:
The water level above the crest was 150 mm. Assume coefficient of discharge in all the cases
to be 0.6.
1. Triangular notch with semi-vertical angle of 45
Solution:
The discharge over the notches will be as follows:
1. Triangular notch
___
Cd 冪 2g t an 2 H5/2 8 ___________________________
8______________________ 0.6 4.43 tan 45o 0.155/2
Q 15 15
0.0123 m3/s.
2. Rectangular notch
___
2 Cd(L0.1nH) 冪2g H3/2
Q ______________________
3
2 0.6(0.20.1 2 0.15) 4.43 0.153/2
_______________________________________
3
0.0175 m3/s
3. Trapezoidal notch
___ ___
Cd 冪2g LH3/2 ____________________
2________________ 8 Cd√ 2g tan H5/2
Q 3 15
2 0.6 4.43 0.15 0.15
___________________________ 3/2
8 0.6 4.43 tan 30o 0.155/2
_____________________________
3 15
0.0154 0.00712
0.02252 m3/s
4. Cippolette notch
___
Cd 冪2g LH3/2
2________________
Q 3
2 0.6 4.43 0.225 0.153/2
_____________________________
3
0.0230 m3/s
Therefore, total discharge 0.0123 0.0175 0.02252 0.0230
0.07532 m3/s 75.32 lit/s
Example 9.11
A 0.5-m-wide SWF and a 0.25-m-wide Parshall flume are used to measure discharge in a
branch canal. The throat level of both the flumes is at the same level. Find the combined
discharge when the head is 0.4 m. Assume Cd equal to 1.706 and 1.72 for SWF and partial
flume, respectively.
Solution:
Q Cd L H3/2 Cd L H1 49 1.706 0.5 0.43/2 1.72 0.25 0.41 49
0.215 0.109
0.324 m3/s
9.7.3 WEIRS
The existing weirs constructed across a stream can be used for the discharge measurement. The
normal types of weirs generally constructed are the following:
(a) broad-crested weir and (b) ogee-shaped weir.
9.7.3.1 Broad-crested weir
Broad-crested weirs are normally of low heights, constructed over the entire length of the channel.
It will act as a broad-crested weir only if B > 2/3 H (Fig. 9.26).
Example 9.12
A broad-crested weir is constructed across a river (Fig. 9.27). The details are as under:
(1) Length of the weir 10 m
(2) Width of the weir 1.0 m
Find the discharge over the weir when the head over the weir is 0.5 m.
DISCHARGE MEASUREMENT 269
Solution:
Q 1.49L H3/2 1.4910.00.53/2 5.26 m3/s.
Example 9.13
A gated ogee spillway has been provided to a dam, the details are as under.
1. Spillway crest level 100.0 m
2. Number and size of gates 4, each 6 m 9 m
270 ELEMENTARY ENGINEERING HYDROLOGY
Solution:
Q Cd (Lnk*H ) H3/2
2.05 (4 9 0.1 8 5) (105 100)3/2
2.05 32 53/2
733.4 m3/s
(k is assumed equal to 0.1.)
A three-dimensional composite model of the stream to a suitable scale may be constructed in the
laboratory. By running the model, a stage-discharge curve may be obtained for the cross section
where the discharge in the prototype is required. This stage-discharge curve may be used for the
prototype by using the discharge scale. This curve may be extrapolated if required.
Before using the stage-discharge curve observed on the model, it should be seen that the model
is proved to prototype conditions.
In this method, a reach of the stream is reduced to a rectangular section. This length of the reach L
should be approximately equal to the bed width of the rectangular channel (Fig. 9.30).
Two transducers are located on each bank at a depth d. The line joining the transducers should
make an angle of 45
with the axis of the stream. These two transducers should emit and also
receive the ultrasonic signals. When an ultrasonic signal is emitted from A to B, it will travel with a
velocity equal to u Vp and the time of travel T1 to travel for this signal from A to B will be
L
T1 _______
(uVp)
DISCHARGE MEASUREMENT 271
Similarly, T2 will be the time of travel for the signal from B to A and can be mentioned as
T2 L/(uVp ).
_____
L
(
___
1 ___
Thus Vd 2 cos T T
1
1
2
)
This method is known as single path gauging method.
Alternatively, a series of transducers may be installed at different depths and thus the vertical
profile of velocity can be obtained to work out the average velocity and consequently the discharge,
more accurately. This is known as Multiple Path Gauging Method.
The method can be used for river width up to 500 m and an accuracy up to 2% can be obtained.
With this method, immediate results are obtained. However, the observations may be affected
due to weed growth, salinity of water, air entrainment, temperature of water, etc.
This method is based on the Faraday’s principle that An Electromagnetic Force is induced in the
conductor (water) when it passes through a normal magnetic field.
A long conductor is buried at the bottom of the stream at right angles to the flow throughout
the cross section. A current I is allowed to flow in this conductor. This will produce a controlled
vertical magnetic field (Fig. 9.31) and this magnetic field will produce an electromagnetic force
(e.m.f.) depending upon the following:
(a) Depth of flow (b) Discharge in the stream (c) Current I
272 ELEMENTARY ENGINEERING HYDROLOGY
Two sensitive electrodes installed on the banks will measure the small voltage produced due to
this e.m.f. The discharge in the stream can be calculated after the calibration as follows:
Ed n
Q K1 ( ______
I K ) 2
When a large river is in spate, it is very difficult to take discharge observations particularly because it is
very difficult to keep the boat steady at a location, to lower the current meter and take the velocity obser-
vations at the same time. A technique called Moving-boat Technique for discharge measurement was
developed by the US Geological Survey in 1969. In this method, the boat is kept moving and velocity
observations are continuously taken. Hence, the method is named moving-boat technique (Fig. 9.32).
The boat is moved at an angle with the cross section P-P to Q-Q, with a uniform velocity Vv.
However, the flow having a velocity Vf will also act on the boat. Thus, the boat will actually move
DISCHARGE MEASUREMENT 273
along the resultant of Vv and Vf , i.e. Vb. Vb must be along a well-defined cross section at right angles
to the flow. Since the velocity of flow Vf goes on increasing from the bank towards the centre of the
stream, the boatman has to change the so that Vb is along the fixed cross section.
Normally, a propeller-type current meter is used, as it is free to move about its vertical axis. The
current meter is maintained at a depth of 0.5 m below the water surface. There is a vertical vane
attached to the current meter so that it aligns itself in the direction of flow. An echo sounder is also
fixed to the boat for recording the depth of flow or the bed profile.
The boatman starts the boat from the bank at a velocity Vv and keeps the boat always along the pre-
determined cross section by changing the angle . To keep the boat along the cross section, he ranges
the boat with the help of the posts P,P and Q,Q located on the banks. During the movement of the boat,
observations of Vv, and the depths of flow at equal intervals of time for various segments are taken.
The discharge in each segment will be dQ
dQ dAVf
(d d )
[ ________
1
2
2
]WV f
(d d )
[ ________ ] V sin cos dt
1 2
v
2
2
Since all the values on the right hand side are observed during the ride of the boat, dQ can be
evaluated.
Total discharge, Q ∑dQ
The operation is repeated for the return journey of the boat. These operations are repeated and the
average value of discharge is calculated. The method is used for streams having more than 3-m depth
of flow. The method is strenuous as it is difficult to move the boat along a fixed cross section in a
high velocity flow. However, there is no alternative but to follow this method for very large rivers.
274 ELEMENTARY ENGINEERING HYDROLOGY
Example 9.14
Observations for discharge measurement were taken on a river by moving-boat technique. The boat
was moved with a constant speed of 1.2 m/s and observations of angle and the depth of flow were
taken every after 30 s. Find the discharge flowing through the river.
Solution:
The discharge calculations are tabulated below.
Serial Time (s) Angle Depth Average Average Sin Cos Discharge
no. (degree) (m) (degree) depth (m)
1 0 0 0 — — — — —
2 30 20 2.6 10 1.3 0.1736 0.9848 9.601
3 30 24 3.2 22 2.9 0.3746 0.9271 43.508
4 30 28 4.8 26 4.0 0.4383 0.8987 68.065
5 30 36 6.2 32 5.5 0.5299 0.8480 106.766
6 30 40 8.0 38 7.1 0.6156 0.7880 148.787
7 30 32 6.2 36 7.1 0.5877 0.8090 145.829
8 30 26 4.4 29 5.3 0.4848 0.8746 97.080
9 30 20 2.0 23 3.2 0.3907 0.9205 49.716
10 30 00 0.0 10 1.0 0.1736 0.9848 7.385
Total 676.837
When a number of discharge observations along with the stage observations are taken at a gauging
site, these observations are plotted on a simple graph, with discharge on the x-axis and stage on
the y-axis. Once a relation of stage and discharge is established, it becomes easy to calculate the
discharge flowing in the channel.
Such a graph is known as Stage Dischar ge Relation or a Rating Curve . Normally, this is
parabolic in nature as shown in Fig. 9.33.
Example 9.15
The rating curve coordinates obtained at a gauging site are as follows:
Stage (m) 4.0 4.9 6.5 7.2 7.7
3
Discharge (m /s) 40 60 100 120 140
Solution:
(i) Let the equation of the rating curve be
Q k dn
Substituting the values of Q 40 m3/s and Q 100 m3/s in the assumed equation
We get, 40 k 4n and 100 k 6.5n
276 ELEMENTARY ENGINEERING HYDROLOGY
The rating curve coordinates were also plotted on a log–log paper. It was noticed that the graph
is a straight line. This straight line was extended beyond and the required details were observed as
follows:
(i) Discharge when the stage is 8.6 m 180 m3/s
(ii) Stage corresponding to a discharge of 200 m3/s 9.0 m
The stream flow data is the basic data. All further designs are based on these data and hence
sufficient number of stream gauging stations should be established. The WMO has suggested the
following criteria.
(a) Arid zones
Minimum 1 station per 5,000–20,000 km2
(b) For temperate Mediterranean and tropical zones
Minimum 1 station per 1,000–2,500 km2
(c) Mountainous areas of temperate Mediterranean and tropical zones
Minimum 1 station per 300–1,000 km2
A stream gauging station should be established at the following locations:
(i) Downstream of a possible dam or a hydraulic structure
(ii) At the inlet or the outlet of a natural lake
(iii) Upstream of a city subject to flood damage
(iv) A point just above the tidal influence
(v) On a tributary just u/s of its confluence with the main river
278 ELEMENTARY ENGINEERING HYDROLOGY
REVIEW QUESTIONS
1. Define discharge. Discuss the various units used for it.
2. What are the main objectives of the measurement of discharge?
3. Discuss the importance of the measurement of discharge in the study of hydrology.
4. What is water potential? Discuss its units.
5. What is a stage? Discuss the different methods followed to measure the stage of a channel.
6. Explain with a neat sketch a water-stage recorder.
7. State the different methods adopted to measure the discharge in a channel and explain
in details any one of them.
8. Discuss the criteria for the selection of a site for discharge measurement.
9. Discuss the procedure followed to evaluate the cross-sectional area and the hydraulic mean
depth at a cross section of a channel.
10. Explain how the Manning’s coefficient is used to find the discharge flowing in a channel.
11. Explain with a neat sketch the velocity profile along a cross section of a channel.
12. Explain why steady uniform flow is assumed in the case of area–slope method of discharge
measurement.
13. Discuss the criteria for the selection of a float for discharge measurement.
14. Describe the procedure followed to evaluate the surface velocity.
15. Discuss the precautions to be taken while taking observations of velocity by floats.
16. Explain the construction of a cup-type current meter. Also explain how it is used to evaluate
the velocity of flow.
17. Explain the construction of a propeller-type current meter. Also explain how it is used to
evaluate the velocity of flow.
18. Discuss the various methods of lowering the current meter to the desired depth for taking the
velocity observations.
19. Discuss the corrections due the flowing water to be applied, when a current meter is used to
observe velocity of flow.
20. What is an isovel? Sketch a pattern of isovels along a cross section of a channel.
21. Explain with a neat sketch the working of a Pitot tube.
22. Explain the salt-titration method of discharge measurement. Also discuss the properties of the
salt/chemical to be used in this method.
23. Explain with neat sketches the different notches that can be used for the measurement of
discharge.
24. Explain with neat sketches the different types of flumes that can be used for discharge
measurement.
25. Explain the hydraulic model method of discharge measurement.
26. Explain with a neat sketch the ultrasonic method of discharge measurement.
27. Explain with a neat sketch the electromagnetic induction method of discharge measurement.
28. Explain with a neat sketch the moving-boat technique for discharge measurement.
29. Discuss the advantages and the disadvantages of ultrasonic method and electromagnetic
method relative to the method by current meter.
30. Write short notes on:
a. Basic equation of discharge measurement
b. Suspension weight gauge method of measuring the depth of flow in a channel
c. Float type water-stage recorder
DISCHARGE MEASUREMENT 279
NUMERICAL QUESTIONS
1. A test was conducted for the measurement of discharge for a stream. The details are as follows:
(i) Distance between two sections 200 m
(ii) Depth of flow at u/s Section A 2.52 m
(iii) Depth of flow at the d/s section 2.61 m
(iv) Water level at Section A 256.45 m
(v) Water level at Section B 256.25 m
(vi) Manning’s coefficient (assumed) 0.035
Assume, Area of flow in m2 8 d2
and wetted perimeter in m 12 d
Here d is in metres.
Calculate the discharge flowing in the stream.
Ans: Q 68.13 m3/s
2. Observations were taken for the velocity of flow in a channel at various depths as follows:
Depth of flow (m) 0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0
Velocity of flow (m/s) 0 0.025 0.12 0.22 0.45 0.73 1.24 1.91 3.35
Find the length of the cable from the index point up to the bed.
Ans: Length of the cable 7.56 m
5. Isovels were drawn on the cross section of a channel covering an area of 25 m2. The percentage
of distribution of the area between the two consecutive isohyets is as follows:
Time (s) 0 30 30 30 30 30 30 30 30
Angle (degree) 0 20 26 34 42 32 24 18 0
Depth (m) 0 2.0 3.0 4.2 5.8 6.6 5.0 3.2 0
V0 2DV0 8D
__________ V0 2DV0 4D
__________
(c) 2 (d) 2
8. The three-point method of the measurement of velocity is the average velocity as
V0 2D2V0 6DV0 8D
________________ V0 2D2V0 4DV0 8D
________________
(a) 4 (b) 4
V0 2D2V0 4DV0 6D
________________ V0 4D2V0 6DV0 8D
________________
(c) 4 (d) 4
9. The unit of water potential is
(a) Million m3 (b) ha m
(c) m per unit area of catchment (d) Any of these three
10. 1 million m3 is equal to
(a) 1 101 ha m (b) 1 102 ha m
(c) 1 103 ha m (d) 1 104 ha m
11. Crest gauge is used to record:
(a) Maximum water level reached (b) Minimum water level reached
(c Average water level (d) Water level at the moment
12. A rating curve of a stream shows the variation of discharge w.r.t.
(a) Depth of flow (b) Velocity of flow
(c) Area of flow (d) All the above
13. The air-line and wet-line corrections depend on
(a) Velocity of flow (b) Unit weight of water
(c) Weight of chain and the current meter (d) All the above
14. The salt titration method is suitable for
(a) Small streams (b) Meandering rivers
(c) Very big streams (d) Rivers in flood
15. The science of discharge measurement is known as
(a) Hydrography (b) Hydrometry
(c) Hydrometeorology (d) None of the above
16. An observation well is necessary when the stage is to be measured by
(a) Vertical gauges (b) Inclined gauges
(c) Float gauges (d) None of the above
17. In the area–slope method, the flow is assumed to be
(a) Unsteady and non-uniform (b) Unsteady and uniform
(c) Steady and non-uniform (d) Steady and uniform
18. The method that is not a direct stream-flow measurement technique is
(a) Area–velocity method (b) Area–slope method
(c) Salt titration method (d) Ultrasonic method
DISCHARGE MEASUREMENT 283
*** ***
Flood Routing
10
Chapter Outline
g 10.1 DEFINITION
Sometimes it becomes essential to know the depth of flow at a specific location of a channel. This
can be done by following the flood-routing procedure.
When a river is in floods, the depth of the flow and the hydrograph at different locations along
the river course are not the same, since the flow in the river is non-uniform and unsteady. In such a
case, the depth of the flow and the hydrograph at each section depend on: (a) the depth of the flow
and the hydrograph at a section u/s of it and (b) the hydraulic characteristics of the river in between
the two sections such as the slope of the river, the length of the reach, the Manning’s coefficient,
the cross-sectional area and so on.
Thus, the depth of the flow and the hydrograph at the lower section can be calculated having
the above information.
Similarly, the outflow from a reservoir will depend on the following factors:
1. The inflow hydrograph
2. The storage capacity of the reservoir
3. The discharging capacity of the spillway weir
FLOOD ROUTING 285
Thus, the outflow hydrograph from a reservoir, whether gated or ungated, can be determined
by knowing (1) the inflow hydrograph, (2) the hydraulic characteristics of the reservoir and (3) the
discharging capacity of the spillway weir.
Flood routing can, therefore, be defined as the procedure for calculating the depth of flow
and the hydrograph at a section knowing the flood hydrograph at the u/s section and the
hydraulic characteristics of the channel in between or the hydraulic characteristics of the
reservoir.
The inflow-hydrograph ordinates at a time interval of t will be I1, I2 and I3, and these will be
known since the inflow hydrograph is known. The outflow ordinates may be Q1, Q2, Q3 and so on.
However, these are not known and are to be calculated. (I I )
_______
Over a time interval t, the inflow volume will be t.
1 2
2
(Q1 Q2)
_________
And, over the same time interval t, the outflow volume will be 2 t.
The difference between the inflow volume and the outflow volume will be the change in the
storage S2 – S1. Here, neither the energy equation nor the momentum equation is applied, but a
relation between inflow–outflow and storage is used.
Thus,
(I1I2)
______ (Q1Q2)
t ________ t S2 S1 (10.2)
2 2
Upper Framework
dQ Scale
I Q K.dQ/dt
dt
Under carriage
Outflow
hydrograph Inflow hydrograph
Origin of
outflow
hydrograph Time
lag
Self-Balancing
Potentiometer P1
R4 R1
Phototube Bulb
dS1
dS1 I
dt
dt
R2
A.C. Line
C1
Transformer C2 R3
P0
dS2 dS1
dt Q
dt
R5
In a river during floods, the flow is non-uniform and unsteady. This type of flow is very difficult
to solve. The hydraulic characteristics vary from stage to stage and also from channel to channel.
There may be a lateral inflow or outflow, also.
All these conditions make the analysis of the problem increasingly difficult. However, neglect-
ing all the above conditions and assuming changes occurring gradually with time, the problem can
be simplified and solved. It is also known as the stream-flow routing.
I
Wedge storage KX (I Q)
A
Q
Prism
storage KQ
Section A B
B O
Assuming that the cross-sectional area of the flow is directly proportional to the discharge, the
prism storage will be dependant only on the outflow and can be equated to K outflow, where K
is called the storage time constant. The dimension of K will be time unit and will depend on the
length and the hydraulic characteristics of the river reach. Normally, it is in hours. It represents the
time of travel wave through the channel reach.
The wedge storage will be a fraction of the volume of the prism (I – Q) x (I – Q), where
x 1. In the rising flood, this part will be positive, since I Q and in the falling flood, this
part will be negative, since I Q. Now, the prism storage K Q and the wedge storage
x K (I – Q). Here, x is called the weighing factor and will be dimensionless having a range of
0 x 0.3.
C0, C1 and C2 are called the routing coefficients. It may be checked that C0 C1 C2 1.
C0, C1 and C2 can be evaluated, if t, K and x are known and then from Eq. (10.5), Q2 can be
calculated. S2 can be calculated from Eq. (10.4), once Q2 is known.
Then, knowing C0, C1, C2, Q2 and S2, we can calculate Q3 and S3.
Knowing C0, C1, C2, Q3 and S3, we can calculate Q4 and S4. Thus, this forward substitution
procedure can be followed to find the outflow-hydrograph ordinates, Q1, Q2, Q3, … .
This outflow hydrograph will be the inflow hydrograph for the next river reach. The inflow and
the outflow hydrographs observed in a specific case are shown in Fig. 10.4.
Inflow hydrograph
Attenuation
Outflow hydrograph
Discharge in m3/s
TI
TO
QI QO
10.3.2 SELECTION OF Δt
Over a span of time t, the average inflow volume is calculated by assuming a straight-line
variation between I1 and I2 and also the outflow volume, by assuming a straight-line relation
between Q1 and Q2. The accuracy will, therefore, depend on t. If t is small, it will give more
accurate results. t should also be less than the time of travel between the reach of the river, i.e
sections A and B.
Normally, t is taken to be 1/2 or 1/3 time of the travel between the reaches of the river
selected.
Care should be taken so that the peak of the inflow hydrograph is not missed while selecting t.
290 ELEMENTARY ENGINEERING HYDROLOGY
Example 10.1
In case of an open channel, the inflow hydrograph at a location is as follows:
Ti (h) m e 0 1 2 3 4 5 6 7 8 9 10
Discharge (m3/s) 10 16 28 45 37 30 24 19 15 12 10
The values of K and x are 4 h and 0.3, respectively. Route the hydrograph.
Solution:
Assume the time interval t 1 h
Therefore,
(Δt/K) 2x
______________ 1/4 2 0.3
______________ 0.35
______
C0 1.65 0.212
2(1x) (Δt/K) 2(1 0.3) 1/4
(Δt/K) 2x
______________ 1/4 2 0.3
______________ ____
0.85
C1 0.515
2(1x) (Δt/K ) 2(1 0.3) 1/4 1.65
2(1 x)(Δt/K )
________________ 2(1 0.3) (1 4) ____
________________ 1.15
C2 0.696
2(1 x) (Δt/K ) 2(1 0.3) (1/4) 1.65
Now,
Q2 C0 I2 C1 I1 C2 Q1 0.212 I2 0.515 I1 0.696 Q1
The step-by-step method for flood routing was followed and the results are tabulated below.
Observations:
1. The maximum inflow and outflow discharges are 45 m3/s and 30.29 m3/s, respectively.
2. The inflow and the outflow base periods are 10.0 h and 17.0 h, respectively.
The values of K and x will have to be evaluated to find the outflow. These can be calculated from
an observed set of inflow and outflow hydrographs in the reach of the river.
The observed inflow and outflow hydrographs in a specific case are shown in Fig. 10.6.
The values of K and x can be calculated by three methods.
First method
___
dS
Now, IQ
dt
and S K [x I (1 x) Q] (10.3)
Therefore,
___
dS
dt [ __
dI
K x (1 x)
dt
dQ
___
dt ]
IQ K [ x (1 x)
dt ]
__
dI dQ
___
Therefore, (10.6)
dt
When the two hydrographs cross each other, then I Q.
Therefore, from Eq. (10.6) we get
dQ/dt
____________
x
(dQ dtdI dt)
292 ELEMENTARY ENGINEERING HYDROLOGY
INFLOW HYDROGRAPH
WATER ENTERING A1 A2
Discharge in m3/sec
THE STORAGE
A1 OUTFLOW
HYDROGRAPH
A2
WATER
LEAVING THE STORAGE
Time in hours
Time in hours
n
m
m K DAYS
n
Storage in m3
Here,
dQ
___
Slope of the outflow hydrograph at the crossing
dt
__
dI
Slope of the inflow hydrograph at the crossing
dt
Once x is evaluated, K can be evaluated from Equation (10.6) by selecting two points on the two
hydrographs for a specific value of T and substituting the values of I, Q, dI/dt, dQ/dt and x.
This method is not accurate since finding a slope to a curve is approximate.
FLOOD ROUTING 293
Second method
The values of K and x can be evaluated from Eq. (10.6), by solving the two simultaneous equations.
This can be done by selecting two points on each of the inflow and the outflow hydrographs for two
specific times T1 and T2, so that I, Q, dI/dt and dQ/dt can be evaluated at these two points. Then the
only two unknowns are K and x.
Third method
From the observed inflow and outflow hydrographs, it can be seen that initially the inflow will be
more than the outflow. At one point, the inflow will be equal to the outflow and then the outflow
will be more than the inflow.
When the inflow is more than the outflow, the difference in the volume will be temporarily
stored in the river reach and when the outflow will be more than the inflow, this storage will be
drained. The storage in the reach will be the shaded part, as shown in Fig. 10.6.
Then, some values of x are assumed and a graph of S versus [x I (1 x) Q] is plotted for each
value of x.
Normally, the graphs for the various values of x will represent loops. But for a certain value
of x, the graph will be a straight line. This value of x is finalized when the loop is straight or almost
straight. The slope of this line will be K. Thus, x and K can be evaluated.
Example 10.2
The observed values of the inflow and the outflow at two cross sections of a stream are as follows:
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
3 10 16 28 45 37 30 24 19 15 12 10 10 10 10 10 10 10 10
Inflow (m /s)
3 10 9 8 11 23 29 30 29 27 24 21 18 15 14 13 12 11 10
Outflow (m /s)
Calculate the values of K and x, for the Muskingum method of flood routing through channel. The
inflow and the outflow hydrographs are as shown in Fig. 10.7.
Solution:
It can be seen that the two hydrographs cross each other at 5.1 h, when the discharge is
29.3 m3/s. The values of K and x can be evaluated by the three different methods as follows:
First method
dQ/dt
_____________
x
(dQ/dt dI/dt)
The values of dQ/dt and dI/dt, where the two hydrographs cross each other, were calculated and
were 2.250 and 6.071, respectively.
______________
2.250 _____
2.250
Therefore, x 8.321 0.27
2.250(6.071)
Now, IQK x [ __
dI
dt
(1 x)
dQ
___
dt ]
0
Two points on the two hydrographs for the same value of time are selected as follows:
At t 3 h,
__
dI dQ
___
I 45 m3/s, Q 11 m3/s, 0, 8.62 and x 0.27
dt dt
Substituting these values in the equation, we get
___
34
K 6.2 5.4 h
Second method
Two points are selected from each of the two hydrographs, and the slopes at these four points are
also calculated as follows:
__
dI dQ
___
T3h I 45 m3/s Q 11 m 3/s 0.0 8.82
dt dt
__
dI dQ
___
T2h I 28 m3/s Q 8 m3/s 13.38 –0.1
dt dt
Now we have,
__
dI
[
I Q K x (1 x)
dt
dQ
___
dt
0 ]
Substituting these values in this equation, we get the following two equations:
45 – 11 K [x 0 (1 – x) 8.82] (i)
and 28 – 8 K [x 13.38 (1 – x) (–1)] (ii)
FLOOD ROUTING 295
______
179.8
Therefore, x 634.72 0.283
____
34
Therefore, K 6.18 5.5 h
Third method
From the plot of the two hydrographs, the additional storage in the channel is calculated for every hour.
A graph of S vs [x I (1 – x) Q], for some assumed values of x, is plotted. For this, x 0.10, 0.30
and 0.50 are assumed. The calculations are tabulated below.
(A) x 0.10
(B) x 0.30
Serial no. S (m3/h) I (m3/s) Q (m3/s) xI (1 x) Q x I (1 x) Q
1 3.5 16 9 4.80 6.30 11.10
2 17 28 8 8.40 5.60 14.00
3 44 45 11 13.50 7.70 21.20
4 68 37 23 11.50 16.10 27.20
5 76 29.3 29.3 8.79 20.51 29.30
6 73 24 30 7.20 21.00 28.20
7 65 19 29 5.70 20.30 26.00
8 54 15 27 4.50 18.90 23.40
9 42 12 24 3.60 16.80 20.40
10 30.5 10 21 3.00 14.70 17.70
11 20 10 18 3.00 12.60 15.60
12 12.5 10 15 3.00 10.50 13.50
13 8 10 14 3.00 9.80 12.80
14 4.5 10 13 3.00 9.10 12.10
15 2 10 12 3.00 8.40 11.40
16 0.5 10 11 3.00 7.70 10.70
17 0 10 10 3.00 7.00 10.00
(C) x 0.50
Serial no. S (m3/h) I (m3/s) Q (m3/s) xI (1 x) Q x I (1 x) Q
1 3.5 16 9 8.00 4.50 12.50
2 17 28 8 14.00 4.00 18.00
3 44 45 11 22.50 5.50 28.00
4 68 37 23 18.50 11.50 30.00
5 76 29.3 29.3 14.66 14.66 29.32
6 73 24 30 12.00 15.00 27.00
7 65 19 29 9.50 14.50 24.00
8 54 15 27 7.50 13.50 21.00
9 42 12 24 6.00 12.00 18.00
10 30.5 10 21 5.00 10.50 15.50
11 20 10 18 5.00 9.00 14.00
12 12.5 10 15 5.00 7.50 12.50
13 8 10 14 5.00 7.00 12.00
14 4.5 10 13 5.00 6.50 11.50
15 2 10 12 5.00 6.00 11.00
16 0.5 10 11 5.00 5.50 10.50
17 0 10 10 5.00 5.00 10.00
FLOOD ROUTING 297
Three graphs of S vs [x I (1 x) Q] were plotted for x 0.10, 0.30 and 0.50. These are shown
in Fig. 10.8.
It can be seen that for x 0.3, the graph is practically a line. Hence, this value of x is accepted.
The slope of the straight line K 53/13 4.07 h. The values of x and K, evaluated by the three
methods, are tabulated below:
x K (h)
First method 0.270 5.40
Consider an ungated spillway of a reservoir, vide Fig. 10.9. It is full up to its spillway crest weir and
flood water enters into the reservoir.
As the inflow starts into the reservoir, it will be temporarily stored in the reservoir and, hence,
the water level in the reservoir will rise depending upon its storage capacity. When the water level
rises, a head over the spillway weir develops. As the head develops, outflow over the weir starts. All
these phenomena will occur simultaneously. This outflow, therefore, will not be equal to the inflow,
since some volume from the inflow is stored in the reservoir.
The reservoir water level will be assumed to be horizontal. The head developed over the weir
will depend on the water level in the reservoir. Thus, the water level in the reservoir, the storage in
the reservoir and the outflow discharge, all will be dependent on each other.
Since reservoir water level is assumed to be horizontal, the flood routing through reservoir is,
therefore, also called level pool routing.
The outflow from the reservoir will, therefore, depend on the following conditions:
1. Inflow hydrograph
2. The absorption capacity of reservoir
3. The discharging capacity of the spillway weir
The outflow hydrograph can be evaluated by the Puls method.
Crest of weir
Inflow
(I1I2)
______ (Q1Q2)
________
2 Δt 2 Δt S2 S1 (10.2)
∆t should be selected such that it will be approximately equal to tp/5, where tp is the time of
maximum inflow discharge.
Care should also be taken that the peak inflow discharge is not missed.
Thus,
(I1I2)
______
2 Δt Average inflow volume in m in time interval ∆t
3
(Q1Q2)
________
And 2 Δt Average outflow volume in m3 in time interval ∆t
Here I1 and I2 are known, since these are coordinates of the inflow hydrograph. Q1 and S1 are
also known at the beginning of the time.
Equation (10.2) can be rearranged as follows:
(I1 I2)
_______ Δt S Q __
Δt
2 Δt S1Q1 __ 2 2 2
(10.7)
2
In this equation, all the terms in the LHS are known. And in the RHS, S2 and Q2 are not
known. However, a relation between the storage and the outflow can be established and Q2 can be
calculated.
300 ELEMENTARY ENGINEERING HYDROLOGY
Elvation in R.L.
Volume in million m3
The storage in a reservoir and the outflow from a reservoir are related to the water level in the
reservoir. These are known as hydraulic characteristics of the reservoir.
These relations are as follows:
1. The storage capacity curve of the reservoir (Fig. 10.10).
2. The discharging capacity of the spillway weir (Fig. 10.11).
Elvation in R.L.
Discharge in m3/s
Ot
S
2
Outflow in m3/s
Ot
S
2
From these two curves, for a specific value of water level, S and Q can be calculated. This means
that S and Q are related to each other, and their values can be calculated for different water levels.
Thus, a relation of S vs Q can be derived as shown in Fig. 10.12.
On the same graph, in addition, a relation Q vs (S 2O ∆t) (Fig. 10.12) can also be established.
Here, 2O ∆t can be calculated, since ∆t is known and O can be read on the y-axis.
The LHS can be calculated from Eq. (10.7) and then the RHS, i.e. (S2 Q2/2 ∆t), can be known.
Q2 and S2 can be calculated from the graphs in Fig. 10.12. From the value of S2, in the graph in
Fig.10.10, the water level in the reservoir can also be calculated. Thus Q2, S2 and the water level in
the reservoir are known.
The next step will be,
(I2I3)
______ Δt
__ Δt
___
2 Δt Q2 2 S3 Q3 2
From this equation, S3, Q3 and the water level in the reservoir can be calculated.
Then the next step will be,
(I3I4)
______ Δt
__ Δt
___
2 Δt Q3 2 S4 Q4 2
From this equation, S4, Q4 and the water level in the reservoir can be calculated.
Thus, following this step-by-step procedure, Q1, Q2, Q3 and so on, can be evaluated, which is
the outflow hydrograph.
Care should be taken, so that the peak of the inflow hydrograph is not missed while selecting Δt.
The inflow and the outflow hydrographs observed in a specific case are shown in Fig. 10.13.
The following observations can be made:
• The maximum outflow discharge is less than the maximum inflow discharge. This is known as
flood attenuation.
302 ELEMENTARY ENGINEERING HYDROLOGY
Inflow hydrograph
To
Qi Qo
Ti
Time in hours
To Ti
Qi Qo
• The time of maximum outflow discharge is more than that of the maximum inflow discharge.
This difference is known as translation or lag.
• The base period of the outflow hydrograph is more than that of the inflow hydrograph.
• The outflow discharge is the maximum when the outflow hydrograph crosses the inflow
hydrograph.
• The inflow volume is equal to the outflow volume.
Case no. 1
In this case, as flood enters the reservoir, the water level in the reservoir will start rising. As
the water level rises, the water will simultaneously start overflowing the spillway weir. Ini-
tially, the outflow will be less than the inflow and the balance will be absorbed temporarily in
the reservoir, resulting in the rise of water level in the reservoir. This will cause an increase
in the outflow.
As the inflow discharge reduces and is less than the outflow, the water level in the reservoir
will start receding. This will reduce the outflow. Finally, first the inflow and then the outflow will
reduce to zero.
As there is no change in the water level in the reservoir, there will not be any change in the stor-
age. The inflow volume and the outflow volume will be equal.
It should be noted that the outflow will start to recede when it crosses the inflow hydrograph.
The inflow hydrograph, the outflow hydrograph and the water level in the reservoir are shown
in Fig. 10.14.
Definition sketch
Reservoir level
Outflow hydrograph
Inflow hydrograph
Time in hours
Fig. 10.14 Inflow and outflow hydrographs and the water level
304 ELEMENTARY ENGINEERING HYDROLOGY
Example 10.3
An ungated spillway has its crest at 550 m. The elevation, storage and discharging capacity of the
spillway are as follows:
Elevation
(m) 550 551 557 562 565 568 570 572 574 575 576 577
Storage
(ha m) 0 3 8 13.6 20.8 31.0 40.0 52.0 66.0 70.0 76.0 80.0
Discharge
0 4 12 20 32 45 62 92 115 140 160 200
(m3/s)
The inflow into the reservoir, when the water level in the reservoir is at 550 m, is as follows:
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Discharge
0 12 50 160 270 240 205 150 95 68 36 23 14 7 0
(m3/s)
Solution:
For the calculations of the flood routing, Δt was assumed to be 1 h. The curves, reservoir level vs
(a) discharging capacity and (b) storage, were plotted as shown in Fig. 10.15.
From Fig. 10.15, the two relations—discharge vs (a) storage and (b) storage O/2 ∆t—were
established and plotted in Fig. 10.16.
Storage Q/2
Serial no. Discharge (m3/s) Storage (ha m) Q/2 ∆t (ha m) ∆t (ha m)
The step-by-step procedure was followed for the flood routing and is tabulated below.
I I
Time I1 I2
( ______
1
2 )
2
∆t
Q1
1
Q
( S ___
2 )
1
2
Q
∆t ( S ___ ) ∆t
2
2ˆ
Q2 S2
Water
level
(h) (m3/s) (m3/s) ha m (m3/s) ha m ha m (m3/s) (ha m) (m)
The inflow and the outflow hydrographs are as shown in Fig. 10.17.
Flood attenuation 270 – 237 33 m3/s
Maximum water level reached 577.5 m
Case no. 2
In this case, when the flood occurs, the reservoir level is lower than the crest level of the spillway.
The spillway is ungated or the gates are yet to be installed.
As the inflow discharge starts flowing into the reservoir, it will be absorbed in the reservoir. The
water level in the reservoir will go on rising and will continue to rise until it reaches the spillway crest
level. Upto this stage, there will not be any outflow and all the inflow will be absorbed in the reservoir.
Similar to Case no. 1, as the inflow increases further, the water level in the reservoir will increase
above the crest level and will develop a head over the spillway weir and the weir will start func-
tioning. The outflow will start flowing over the weir. The inflow and the outflow will not be equal.
Initially, the outflow will be less than the inflow and the difference will be temporarily absorbed in
the reservoir resulting in an increase in the reservoir water level. This will increase the outflow.
As the inflow reduces and is less than the outflow, the water level in the reservoir will also reduce.
A stage will be reached when the outflow and the inflow will reduce to zero and the water level in the
reservoir will be up to crest level.
In such a case, the inflow volume will be the summation of the outflow volume and the volume
of water absorbed into the reservoir initially from the original water level up to the crest of the
spillway weir.
In this case, also, the outflow will start receding when the inflow and the outflow hydrographs
cross each other.
The inflow hydrograph, the outflow hydrograph, the water absorbed in the reservoir and the water
level in the reservoir are shown in Fig. 10.18.
Case no. 3
In this case, the spillway gates are installed. There is no flood lift and the water level is up to the
FRL or the HFL. (Normally, the water level in the reservoir is kept as high as possible, so that the
storage in the reservoir is the maximum.)
As the inflow occurs, the spillway gates will be operated partially and the outflow will be kept
equal to the inflow, thus maintaining the water level constant at the FRL. As the inflow increases,
the gates will be operated and opened more and more, keeping the water level at FRL. Thus the
outflow will be equal to the inflow. For the maximum discharge, the gates will be fully opened and
there will be no change in the water level. It will be at FRL only.
As the inflow decreases, the gates will be lowered keeping the water level constant at FRL, and
the inflow will be equal to the outflow. When the inflow stops, the gates will be closed completely
and the outflow will also stop. In such a case, since there is no flood lift, the outflow will always be
equal to the inflow and the water level will be constant at FRL.
Thus, there is no flood routing in this case and the inflow and outflow hydrographs will be the
same as shown in Fig. 10.19.
Case no. 4
In this case, there is a flood lift provided in the design. HFL is higher than the FRL by an amount
equal to the flood lift. The water level is at the FRL, with gates completely closed.
As the inflow occurs, the gates will be operated partially to maintain the water level at FRL so
that the inflow is equal to the outflow. As inflow increases, the gates will be operated such that
the water level is maintained constant at the FRL. Thus the outflow will be equal to the inflow.
308 ELEMENTARY ENGINEERING HYDROLOGY
Definition sketch
Reservoir level
Outflow hydrograph
Inflow hydrograph
0
Time in hours
Fig. 10.18 Inflow and outflow hydrographs with the reservoir and crest levels
A stage will be reached when the gates will be completely opened with the water level is at the
FRL, and the outflow will be equal to the inflow. As inflow increases still further, the water level
in the reservoir will increase and the outflow will also increase, but will be less than the inflow.
The water level will rise above the FRL and will reach HFL. when the outflow will be the maxi-
mum. This point will be reached when the inflow and outflow hydrographs cross each other.
As the inflow reduces further, a stage will be reached such that the outflow is equal to the inflow.
From this point, the outflow and the water level will start receding, but the inflow will be less than
the outflow. When the water level reaches the FRL, the gates will have to be operated to maintain the
water level in the reservoir constant at FRL. In this case, the inflow will be equal to the outflow.
Finally, the gates will be closed completely; the inflow and the outflow will stop simultaneously,
with the water level in the reservoir at the FRL.
Since the water level remains unchanged, the inflow volume and the outflow volume will be
equal. The inflow and the outflow hydrographs and the water level in the reservoir are shown in
Fig. 10.20.
Case no. 5
In this case, the spillway is gated with a flood lift provision in the design. The HFL is above the
FRL, by a margin equal to the flood lift. The water level in the reservoir is below the FRL.
FLOOD ROUTING 309
HFL=FRL
Definition sketch
Discharge in m3/s
Inflow outflow
Reservoir level F R L H F L
0 Time in hours
Fig. 10.19 Inflow and outflow hydrographs when there is no flood routing
As the inflow starts into the reservoir, the water level will go on increasing, absorbing the flood in
the reservoir. The gates will remain in the closed position and there will be no outflow. When the water
level reaches the FRL, the spillway gates will be operated partially, maintaining the water level at the
FRL. Here the inflow will be equal to the outflow. A stage will be reached when the gates are fully open
and the water level is at the FRL. This is the limit when the inflow is equal to the outflow.
For the higher inflow, the water level will go on increasing absorbing water in the reser-
voir above the FRL, with a result that both the water level in the reservoir and the outflow will
increase. The inflow after reaching the maximum will start receding. A stage will be reached
when the outflow will be equal to the inflow. From this point, the water level and hence the
outflow will start receding. When the water level reaches the FRL, the gates will be operated
partially so that the water level will not drop below the FRL. In this case, the inflow and the
outflow will be equal.
Finally the inflow will stop. The gates will be closed the outflow will also stop. The inflow
volume will be equal to the summation of the outflow volume and the volume absorbed in the
reservoir from the initial water level up to the FRL.
Figure 10.21 shows the inflow and the outflow hydrographs, the volume absorbed by the
reservoir and the water level in the reservoir.
310 ELEMENTARY ENGINEERING HYDROLOGY
HFL
FRL
Definition sketch
Inflow hydrograph
Discharge in m3s
Outflow hydrograph
Reservoir level
HFL
FRL
0
Time in hours
Fig. 10.20 Inflow and outflow hydrographs and water level in the reservoir
FLOOD ROUTING 311
HFL
FRL
Initial
reservoir level
Definition sketch
Reservoir level
HFL
FRL
Time in hours
Fig. 10.21 Inflow and outflow hydrographs and the water level in the reservoir
312 ELEMENTARY ENGINEERING HYDROLOGY
REVIEW QUESTIONS
1. Define flood routing. Explain its practical use.
2. Which are the different methods used for flood routing?
3. Explain the basic equation used in hydrologic flood routing.
4. Which are the basic equations used in hydraulic flood routing? Discuss them in brief.
5. Explain the Muskingum method of flood routing.
6. How is the value of ∆t in the Muskingum method selected?
7. How are the values of K and x in the Muskingum method determined?
8. Explain the Puls method of flood routing through the reservoir.
9. Discuss the basic data required for the Puls method.
10. Describe how the relation between the storage capacity of the reservoir and the discharging
capacity of the spillway weir established.
11. Point out the specific observations between the inflow hydrograph and the outflow hydrograph
in the case of flood routing through the reservoir.
12. Discuss the possible practical cases of flood routing for a dam having a spillway in the following
situations:
a. Without gates b. With gates without a flood lift
c. With gates with a flood lift
13. Comment on the following fact: In the case of flood routing through the reservoir, the outflow
hydrograph crosses the inflow hydrograph when the outflow discharge is the maximum.
14. Discuss how the gates provided for a spillway affect the outflow hydrograph.
15. Discuss how the flood lift provided for a spillway affects the outflow hydrograph.
16. What is the appropriate time interval to be adopted in reservoir routing?
17. In the Muskingum-storage equation, what do the parameters K and x represent?
18. Derive the expressions for the routing coefficients, C0, C1 and C2 in the Muskingum method.
19. Write short notes on the following:
a. Prism and wedge storage b. Flood routing through channel
c. Flood routing through reservoir d. Storage time constant
e. Flood attenuation f. Weighing factor
g. Level pool routing h. Reservoir lag
20. Differentiate between the following:
a. Hydraulic routing and hydrologic routing b. Channel routing and reservoir routing
c. Prism storage and wedge storage d. Inflow hydrograph and the outflow
hydrograph
NUMERICAL QUESTION
Time (h) 0 1 2 3 4 5 6 7 8 9
Discharge
8.0 14.0 26.0 43.0 35.0 25.0 18.0 15.0 12.0 8.0
(m3/s)
The values of K and x for this specific channel section are 3 h and 0.3, respectively. Find the
hydrograph at the bridge location and the flood attenuation.
FLOOD ROUTING 313
Time (h) 0 1 2 3 4 5 6 7 8 9 10 11
Discharge
8.0 7.06 7.83 12.10 25.13 30.44 29.44 25.53 21.96 18.77 14.65 12.11
(m3/s)
12 13 14 15 16 17 18
10.53 9.56 8.96 8.59 8.36 8.21 8.13
9. The condition satisfied by the three routing coefficients in the Muskingum method is
(a) C0 C1 C2 0 (b) C0 C1 C2 1
(c) C0 C1 C2 10 (d) C0 C1 C2 100
10. Out of the following which is the Muskingum-storage equation?
(a) S [x Q (1 –x) I] (b) S [x Q – (1 –x) I]
(c) S [x I (1 –x ) Q] (d) S [x I – (1 –x) Q]
11 . The Muskingum-routing equation is
(a) Q2 C0 Q1 C1 I2 C2 I1 (b) Q2 C0 I2 C1 I1 C2 Q1
(c) Q2 C0 I1 C1 Q1 C2 I2 (d) Q2 C0 Q1 C1 I1 C2 I2
*** ***
Groundwater
11
Chapter Outline
g 11.1 DEFINITION
Groundwater may be defined as water that occurs below the surface of the earth. It is also known
as sub-surface water.
Water present in the rock at the time of its formation is known as connate water. It is highly
saline. It is also known as fossil water or interstitial water.
Juvenile water is formed chemically within the earth. It is in a very small quantity and is also
known as primitive water.
The level below which the soil is saturated with water is called groundwater table. The groundwater
is divided into two zones (1) unsaturated zone and (2) saturated zone.
These two zones are separated by water table where the pressure is atmospheric, as shown
in Fig. 11.1.
Ground surface
Capillary zone
Water table
Zone of
saturation
Ground water
Bed rock
Groundwater occurs at various locations below the earth surface, depending on the formation of
ground. The various formations are shown in Fig. 11.2.
318 ELEMENTARY ENGINEERING HYDROLOGY
Non-flowing
artesian well
Piezometric Gro
surface und
surf
Flowing artesian we ace
Water ll
table Permeable unconfined aquifer
Water table
Impermeable strata
11.3.1 AQUIFER
Formation of ground that contains water and may transmit water in usable quantity is known as
aquifer. The aquifer may be unconfined or confined.
Unconfined aquifer
An aquifer where the water table is the upper surface limit and extends below till the impermeable
rock strata is called the unconfined aquifer, as shown in Fig. 11.3. It is also known as free aquifer,
phreatic aquifer, water table aquifer and non-artesian aquifer.
The water level in a well is an unconfined aquifer will be upto the watertable.
Ground surface
Water table
Unconfined
aquifer Well
Permeable strata
Confined aquifer
When an aquifer is sandwiched between two impermeable layers, it is known as a confined aquifer.
It will not have a free water table, and the aquifer will be under pressure as shown in Fig. 11.4.
If there is a well in this layer, the water in the well will rise up to the piezometric head. If the
piezometric head is above the ground, then the water from the well in this layer will flow over the
ground. Then the well is called as free-flowing well or flowing well.
If the water level in a well in this layer is above the upper-confining layer level, but below the
ground level, then such a well is called as artesian well.
Groundwater may move to the ground surface at a very small rate through faults, permeable
material in joints, discontinuities, and so on. It is then called a spring.
11.3.2 AQUICLUDE
A geological formation that may contain water because of high porosity but cannot transmit it is
called an aquiclude.
11.3.3 AQUITARD
A geological formation that has poor permeability, but through which seepage is possible, and is
insignificant as compared to an aquifer is known as an aquitard.
11.3.4 AQUIFUGE
A geological formation that neither contains nor transmits water is called an aquifuge. Solid rock
is an aquifuge.
Land Surface
Impervious stratum
General water table
Unconfined Aquifer
11.4.1 POROSITY
Any soil sample contains some pores or voids. The porosity of the sample is the ratio of the volume
of voids and the volume of soil sample. It is a dimensionless number and is normally denoted by n
and is expressed as percentage.
Vv
___
Thus, porosity n V 100
where, Vv Volume of pores or voids
V Volume of sample
Porosity is classified as follows:
20% or more—large, 5–20%—medium and less than 5%—small.
The primary porosity of a material is one that existed when the material was formed. The
secondary porosity results from fractures, joints, solution channels and so on.
11.4.2 PERMEABILITY
Permeability is the property of an aquifer to transmit water through its pores.
The horizontal permeability and the vertical permeability may differ.
11.4.3 DARCY’S LAW
The velocity of the flow of water in an aquifer is very low and hence the flow is laminar. Darcy’s
law is applicable here. The law states, The rate of flow per unit area of an aquifer is proportional
to the gradient of potential flow in the direction of flow.
GROUNDWATER 321
___
dh
i.e., v K i K
dl
where, v Velocity of flow in m/s
K Coefficient of permeability
i Hydraulic gradient
___
dh
dl
The negative sign indicates that the flow is in the falling head.
Coefficient of permeability, K, is also known as hydraulic conductivity and has the dimensions
as those of the velocity of flow. It depends on the fluid property as well as the property of the
aquifer soil sample. K is given by the following identity:
cd 2g
_____
K
Some quantity of water will be retained by the soil sample and will not be drained under gravity.
This is due to molecular attraction and surface tension. Specific retention is defined as the ratio of
the volume of water retained divided by the volume of sample. It is denoted by Sr.
W
Thus, Sr ___f
V
where, Sr Specific retention
Wf Volume of water retained by the soil sample
V Volume of soil sample
Sr is dimensionless. It is also known as field capacity.
Naturally, n Sf Sr
where, n Porosity
Sf Specific yield
Sr Specific retention
11.4.6 TRANSMISSIBILITY
It is the product of the coefficient of permeability and the thickness of the aquifer.
Thus, T K b
It has the dimensions of m2/s. It is also known as transmissivity.
For a unit width
Q area velocity b 1 K ( di
dh)
___
T (
___
di
dh )
K Hydraulic conductivity of the soil
T Transmissibility of the soil
Consider two rivers, A and B, flowing in parallel at a distance L from each other with a water level
difference of h, as shown in Fig. 11.6.
Water will flow from river A to river B in only one direction, i.e. ‘x’ direction, and hence it
is also known as unidirectional flow.
The velocity of flow will depend on the following factors:
1. h Difference in the water levels in the two rivers
2. K Hydraulic conductivity of soil
For the analysis of the flow, the following assumptions are made:
1. The flow is steady
2. The flow is horizontal and uniform in vertical direction
The characteristics of the flow will vary depending on the aquifers in between the two rivers
as mentioned below:
1. Unconfined aquifer
2. Confined aquifer
GROUNDWATER 323
Plan
__
K __
K __
K
q L 2 h21 2 h22 = 2 (h21 h22)
___
K
Therefore, q 2L (h1 h2) (h1 h2)
Example 11.1
Two rivers, A and B, are separated by an unconfined aquifer having a hydraulic conductivity of
10 m/day. The rivers flow at a uniform distance of 1000 m. If the depth of flow in both the rivers is
15 and 10 m, respectively, find the rate of flow from A to B.
5m
1000 m
15 m
K 10 m/day
10 m
Datum
A B
Solution:
The flow between the two rivers is shown in Fig. 11.7.
___
K ________
10
Therefore, q 2L (h21 h22) 2 1000 ( 152 102) 0.625 m3/m /day
h1
h h2
B
L
Datum
(h1 h2)
_________
The velocity of flow v K L
(h1 h2)
________ (h1 h2)
________
Therefore, q B K L T L
since T Transmissibility B K
Example 11.2
Two rivers, P and Q, are separated by a confined aquifer of thickness 2 m. The hydraulic conduc-
tivity of the confined aquifer is 15 m/day. If the depth of flow in both the rivers is 15 and 10 m,
respectively, find the rate of flow from P to Q.
Solution:
The flow between the two rivers is shown in Fig. 11.9.
Impervious layer
15 m
10 m
2m
1000 m
Datum
BK
______ 2 15
________
Q L ( h1 h2) 1000 ( 15 10) 0.15 m /m /day
3
A well is a hole or a shaft, normally vertical, excavated in the ground for bringing water to the
surface. Wells are generally classified as: (a) open wells and (b) tube wells.
Shallow well
An open well that is constructed to tap water from the topmost water-bearing stratum, i.e. from the
unconfined aquifer, is called shallow well. The bottom of such a well will not rest on an impervious
layer as shown in Fig. 11.10. The water level in such a well will be equal to the level of the water
table. These wells are also known as water table wells, unconfined wells or gravity wells.
Deep well
A deep well is constructed to rest on an impervious layer, drawing water from the aquifer below it.
The impervious layer provides a support to the wall of the well, as shown in Fig. 11.11. The yield
from a deep well is comparatively more than a shallow well and is relatively pure. The water level
in such wells is equal to the piezometric head in the water-bearing strata.
This classification of wells is purely technical and has nothing to do with the actual depth of the
wells. A shallow well may be deeper than a deep well. Based on lining, open wells can also be classi-
fied as follows: (1) unlined wells, (2) wells with pervious lining and (3) wells with impervious lining.
Ground surface
Water table
Impervious layer
Aquifer
Strainer well
This is the most common and widely used type. The term ‘tube well’ is generally referred to this
type of well. In this type, a special designed wire mesh with filter pack covering the water-bearing
strata is used to draw water free from silt and sand. The pipe in this portion is perforated. The pipe
is generally plugged at the bottom, to permit settlement of the sand passing through the strainers.
In such type of wells, the inflow is radial. Figure 11.12 shows a strainer well.
Ground surface
Water table
Blind pipe
Impervious layer
Cavity well
This is a special type of tube well in which water is not drawn through the sides but is drawn
through the bottom of the well where a cavity is formed as shown in Fig. 11.13. Normally this type
is used when a clay layer is met with. The zone of the flow in this case is spherical.
Compressed air
Main pipe
Casing pipe
(removed later)
Top soil
Gravel shoulder
Slots
Water table
Water table
Stream
Effluent stream
Stream
Water table
Influent stream
The theory of groundwater movement into a well is commonly known as well hydraulics. The move-
ment of groundwater is in accordance with the Darcy’s law, coupled with the hydraulic principles.
When the water table is higher than the bed of the stream, groundwater will flow to the stream.
Such a stream is known as effluent stream. On the other hand, when the water table is lower than
the bed of the stream, water will flow from the stream to the groundwater. Such a stream is called
influent stream. Both types are shown in Fig. 11.15.
In both the cases, the flow will depend on the soil characteristics. Similar is the case of a well.
When water flows into a well, the well is known as discharging well; when water flows from the
well to the adjacent ground, it is called recharging well.
The groundwater movement into a well is classified as under:
1. Steady radial flow in an unconfined aquifer.
2. Steady radial flow in a confined aquifer.
In both the cases, the following assumptions are made:
• The flow is horizontal and uniformly distributed in a vertical section.
• The velocity of flow is proportional to the tangent of the hydraulic gradient.
• The well fully penetrates the aquifer.
• The discharge taken from the well is constant.
• The aquifer is homogeneous and isotropic.
When water is pumped out from this well at a constant rate, Q, the water level in the well will
naturally go on receding. As water level lowers, water starts flowing in radially from the adjoining
area. The inflow will depend on the difference between the water table and water level in the well.
Since this difference will be initially low, the inflow will be less than the outflow. Thus, the water
level in the well will go on receding as pumping continues. A stage will come when the water level in
the well will be steady, indicating that the inflow is equal to the outflow. This will be a steady state.
The decrease in the water level in the well is known as drawdown or depression head. The water
table in the adjoining area, due to radial flow, will fall and assume a conical shape. This is called
cone of depression. The area up to which the cone of depression extends is called the area of influ-
ence, and the radius at which the cone of depression starts is called the radius of influence.
If there are any other wells in the cone of depression, the water levels in these wells will also fall
depending on their distance from the pumping well. Such wells are known as observation wells.
In the steady state, the depth of the flow at a distance r from the well will be h. Being a radial
flow, the area of flow will be equal to 2rh.
According to Darcy’s law, the velocity of flow K i, i being the hydraulic gradient dh/dr, i.e.,
the slope of the drawdown curve at a distance of h.
___
dh
Therefore, the velocity of flow K i K
dr
___
dh
Now, Q area velocity 2rh K
dr
Q
____ __
dr
or hdh 2K r
The depths of the flow are h1 and h2 in the two observation wells at distances r1 and r2. Integration
of the above equation in these limits yields
r1
h1
Q
____ __
dr
h2
hdh 2K r2 r
Q
H h1
h2 h
hw
r
r1 rw
r2 R
hdh
H R
Q
____ __
dr
hw rw 2K r
Q
____ ___
R
Therefore, (H 2 h2w)/2 2K loge r
w
K(H 2 h2w)
____________
Q
loge (R/rw)
where, H Thickness of the aquifer in m
hw Depth of flow at the pumping well in m
R Radius of influence in m
rw Radius of the well in m
The numerator K (H 2 h2w), in this equation, can be written as K (H hw) (H hw).
Here, (H hw) sw , will be the drawdown in the well, and (H hw)/2 will be the average
thickness of the aquifer. Thus, (H hw)/2 K T.
Therefore, K (H 2 h2w) will be equal to 2T sw , where T is the transmissibility of the
aquifer.
The formula will reduce to
2Tsw
_________
Q
loge( R/rw )
These equations were first derived by Dupuit in 1863, and then modified by Thiem in 1906.
An empirical formula may be used to evaluate R as follows:
__
R 3000 sw √ K
332 ELEMENTARY ENGINEERING HYDROLOGY
Example 11.3
Water is pumped out at the rate of 2400 lit/min, from a well of 0.25 m diameter, penetrating fully
in an aquifer of 32 m thickness. The drawdown observed in two adjoining wells at 15 m and
110 m from the pumping well is 7.0 m and 0.5 m, respectively. Determine the average hydraulic
conductivity.
Solution:
K(h21 h22)
___________
Q
loge (r1/r2)
Example 11.4
Water was pumped out from a well of diameter 1.2 m, from an unconfined aquifer of 50 m thickness.
The well penetrates fully in the aquifer and the depth of water in the well is 40 m. The effect of
pumping was reaching up to a distance of 450 m. Find the discharge if K 10 m/day.
Solution:
K (H 2 hw)
__________
2
Q
loge (R/rw)
GROUNDWATER 333
Ground surface
Q
Draw-down
h1 curve
h2 H
hw
r
r1 rw
Aquifer b r2
R
dh
H R
Q
_____ __
dr
r
hw rw 2bK
Q
____ __
R
H hw 2T l oge r
w
Now, H hw sw
2Tsw
_________
Therefore, Q
loge (r/rw)
Example 11.5
Water was pumped out from a well in a confined aquifer 10 m thick, having a hydraulic conductiv-
ity of 1.5 m/day. The drawdown observed in the two adjoining wells at 10 m and 60 m from the
pumping well was 3 m and 0.05 m, respectively. Find the constant rate of pumping.
Solution:
2T (s2 s1)
___________
Q
loge (r1/r2)
Here, T Kb 1.5 10 15 m2/day 15 86400 1.736 104 m2/s
s2 3 m and s1 0.05 m
r1 60 m and r2 10 m
loge (r1/r2) loge (60/10) 1.79
GROUNDWATER 335
Example 11.6
A well 0.5 m in diameter penetrates fully a confined aquifer of thickness 25 m, having a conductivity
of 1.6/day. It is expected that the drawdown in the well be limited to 2.5 m. Assuming that the
radius of influence is 300 m, find the maximum discharge that can be pumped out of the well.
Solution:
2Tsw
_________
Q
loge (R/rw)
Here, T K b 1.6 25 m2/day 1.6 25 86400 0.00462 m2/s.
R 300 m and rw 0.5 2 0.25 m and sw 2.5 m
loge (R/rw) loge (300/0.25) 7.090
Substituting the above values, we get
0.00462 2.5
_________________
2
Q 7.09
0.01026 m3/s 613.8 lit/min
different media, with or without water, is different. By analysing and properly interpreting the col-
lected data of velocities of the shock waves, the possibility of groundwater can be assessed.
Safe yield from a groundwater basin is the amount of water, which can be annually withdrawn from
a basin without producing any undesirable effect. This safe yield is the amount of water that enters
into the basin annually. Any withdrawal more than the safe yield is called overdraft.
When water more than the safe yield is withdrawn from a basin, the excess comes from the stor-
age within the aquifer. Such permanent depletion of storage within the aquifer is known as mining
of groundwater.
When pumping is started from a well, the water level in the well goes on receding. This drop in
water level in the well is known as drawdown or depression head.
GROUNDWATER 337
When the depression is more than a critical limit, dislodging of soil particles from the adjoining
strata is noticed due to high velocity of water entering into the well. The depression head, when
such dislodging is noticed, is known as critical depression head and corresponding to this head, the
yield is known as maximum yield or critical yield.
Safe depression head, also known as working head, is the head when there is no possibility of
any dislodging of soil particles. Normally, it is taken as one-third of the critical depression head,
and the yield corresponding to this working head is known as maximum safe yield.
The rate of inflow in a well Q is proportional to the depression head H.
__
k
i.e., Q H or Q k H or Q A A H or Q = C A H
where, k is the constant of proportionality having dimensions of m2/s and A is the area of
cross section of the well. Normally, A is taken as 4/3 of the cross-sectional area of the well, since
additional cavity is formed at the bottom of the well.
Here, C k/A, is known as specific yield or specific capacity. Its dimensions are T 1.
Example 11.7
A pumping test was conducted for an open well of diameter 3.5 m. The water was pumped out at a
constant rate of 300 lit/min. Find the specific yield.
Solution:
QCAH
Here, Q 300 lit/min 300 10 60 0.005 m3/s
–3
D 3.5 m
D2 ________
____ 3.52
A 4 4 9.6211 m2
Substituting the above values, we get
0.005 C 9.6211 3.5
∴ C 1.48 10–4/s 0.5345/h.
338 ELEMENTARY ENGINEERING HYDROLOGY
We get
__
K
A ___
dh
h
0 h1
GROUNDWATER 339
___
K T h2 h1
Or A [t]0= loge [h]h1 = loge [h]h2
___
K h
Therefore, A T loge ___1
h2
__
K
Since A C
__
1 h1
__ _____
2.303 h1
__
C T loge = T log10
h2 h2
Example 11.8
A recuperation test was conducted on an open well 5.0 m in diameter. The water levels observed
during the test were as follows:
1. Groundwater table level 250.0 m
2. Water level when the pumping was stopped 243.0 m
3. Water level in the well 2 h after pumping was stopped 245.0 m
Find the safe yield of the well if the working head is 3.0 m.
Solution:
QCAH
Here, h1 250 243 7.0 m and h2 250 245 5.0 m
T 2.0 h
H 3.0 m
340 ELEMENTARY ENGINEERING HYDROLOGY
D 5.0 m
52
D2 ______
_______
A 4 4 19.64 m
2
_____
2.303 h1 _____
__ 2.303 __
7
C T log10 2 log10 5
h2
1.152 0.1461 0.1683/h
Substituting the above values, we get
Q 0.1683 19.64 3.0
9.918 m3/h 2.755 lit/s
When two wells are located close to each other, and when pumping is started in both the wells, a
cone of depression will be developed for each well. If these cones of depression overlap each other,
the yield from each of the wells will naturally be reduced. When two or more cones of depression
of the pumping wells overlap, affecting their yield, the phenomenon is known as well interference
and is shown in Fig. 11.19.
Q1 Q2
Draw-down curve
Draw-down curve Composite for Q2 only
for Q1 only draw-down
curve
h B
b
Confined aquifer
If pumping of fresh groundwater is done near the sea, the sea water being heavier tries to occupy
the lower level space vacated by the fresh water and causes serious water quality problems. This is
known as sea water intrusion. There is a plane that separates the fresh water from the sea water, and
it is known as interface. Normally, it is parabolic in shape as shown in Fig. 11.20.
Sea water is slightly heavier than the fresh water due to the dissolved salts. Specific gravity of
sea water is 1.025 and that of fresh water is 1.000. As the fresh groundwater table increases and
GROUNDWATER 341
Ground surface
Water table
Sea level
Y
Fresh ground water Ocean
Approx. 40Y
Salt ground water
Interface
Natural equilibrium
Pumped well
Cone of depression
Fresh ground
water Ocean
Interface
Effect of pumping
approaches the earth surface by 1.00 m, the sea water interface is pushed back by 40.00 m and
vice versa. When the sea water intrusion is noticed, it can be pushed back by recharging the dug
wells near the seashore by fresh water.
Water may be withdrawn from an aquifer by (1) pumps, (2) discharge to surface water systems or
by (3) evaporation and transpiration (to a minor extent). If this water is not replenished, the aquifer
storage will reduce resulting in depletion of groundwater table.
Normally, this replenishment occurs naturally by the following ways:
• Infiltration from precipitation
• Surface bodies like reservoirs, streams, and so on
Artificial recharge is adhered to meet the deficit of withdrawal and replenishment or even
to increase the aquifer yield. It is also followed to control the saltwater intrusion along the sea
coast.
The following recharge methods are employed, taking into consideration the geological situa-
tion and economic consideration.
342 ELEMENTARY ENGINEERING HYDROLOGY
REVIEW QUESTIONS
1. Discuss the different sources of groundwater.
2. What are the different parameters of an aquifer?
3. Explain the Darcy’s law.
4. Explain with the help of neat sketches the following:
a. Aquifer b. Aquiclude
c. Aquitard d. Aquifuse
5. Derive an expression for the constant discharge pumped from a well in an unconfined aquifer
in terms of the following:
a. Drawdown in two adjacent wells
b. Drawdown of the pumping well
6. Derive an expression for the constant discharge pumped from a well in a confined aquifer in
terms of the following:
a. Drawdown in two adjacent wells
b. Drawdown of the pumping well
7. What is salt water intrusion? How can it be controlled?
8. Discuss the different methods of groundwater exploration.
9. Why is recharging of wells necessary? Discuss the different methods.
10. Explain the pumping test to estimate the safe yield from an open well.
11. Explain the recuperation test to estimate the safe yield of an open well.
12. What are the peculiarities of groundwater?
13. Discuss the assumptions made in the analysis of steady radial flow into a well.
14. Write short notes on the following:
a. Perched water table b. Cone of depression
c. Well interference d. Salt water intrusion
e. Safe yield from a well f. Hydraulic conductivity
g. Storage coefficient h. Porosity of an aquifer
i. Overdraft from a basin j. Mining of groundwater
k. Capillary zone in an aquifer l. Area of influence
m. Critical de pression he ad n. Specific capacity
15. Differentiate between the following:
a. Confined aquifer and unconfined aquifer b. Zone of saturation and zone of aeration
c. Specific yield and specific retention d. Open well and tube well
e. Shallow well and deep well f. Strainer well, cavity well and slotted well
GROUNDWATER 343
g. Pumping test and recuperation test h. Discharging well and charging well
i. Pervious layer and impervious layer j. Groundwater table and piezometric level
k. Primary porosity and secondary porosity l. Artesian well and flowing well
m. Safe depression head and critical n. Cone of depression and depression
depression head head
NUMERICAL QUESTIONS
1. Water was pumped at a constant discharge from a well of diameter 0.5 m located in an 35-m
thick unconfined aquifer having a hydraulic conductivity of 6 m/day. The drawdown observed
in two observation wells located at 16 m and 115 m from the pumping well were 5.0 and 0.5 m
respectively. Find the constant discharge.
Ans: Q 0.032 m3/s
2. Water at the rate of 100 lit/min from a well of diameter 1.0 m, located in a confined aquifer 8.0-m
thick, was pumped. The drawdown observed in two observation wells at 10.0 m and 60.0 m from the
pumping well were 3.5 m and 0.1 m respectively. Find the hydraulic conductivity of the aquifer.
Ans: K 1.5 m/day
3. A recuperation test was conducted on a well of 4.0-m diameter. The water levels observed were
as follows:
(i) Ground w ater t able l evel 300 m
(ii) Water l evel w hen t he pum ping w as s topped 292 m
(iii) Water levelling the well 2 h after the pumping was stopped 295 m
Find the safe yield of the well if the working head is 3.5 m.
Ans: Q 10.23 m3/h
*** ***
Appendix I
INTERNATIONAL SYSTEM OF UNITS [SI]
The 11th General Conference on Weights and Measures (1960) adopted the name Système Interna-
tional d’Unités (International System of Units, international abbreviation SI), for the recommended
practical system of units of measurement. The SI is not static but evolves to match the world’s
increasingly demanding requirements for measurement.
There are seven well-defined units called the base units, which by convention are regarded as
dimensionally independent. These are:
‘The metre, The kilogram, The second, The ampere, The kelvin, The mole, and the Candela.’
The SI list presented here contains only those units that are used in engineering hydrology.
As per the 11th Conference Générale des Poids et Mesures (CGPM), the dimensions of the base
units related to engineering hydrology are as follows:
Metre Kilogram Second
The metre is the length of The kilogram is the unit of The second is the duration of
the path travelled by light in mass; it is equal to the mass of 9 192 631 770 periods of the
vacuum during a time interval the international prototype of radiation corresponding to the
of 1/299 792 458 of a second. the kilogram. transition between the two
hyperfine levels of the ground
state of the caesium 133 atom.
The units formed by combining base units according to the algebraic relations linking the cor-
responding quantities are called the derived units.
The details of the basic and the derived units related to engineering hydrology are as follows:
*** ***
Appendix II
Length
1)1 in 2.54 cmi)1 cm 0.394 in
2)1 ft 0.305 mii)1 m 3.281 ft 1.093 yd
3)1 yd 0.914 miii)1 km 0.621 mile
4)1 mile 1.609 km
Area
1)l in 2
6.452 cm2 i)1 cm 2
0.155 in2
2)1 ft 2
0.093 m2 ii)1 m 2
10.764 ft2 1.195 yd2
3)1 mile 2
2.588 km2 iii)1 km 2
0.386 mile2
4)1 acre 0.405 haiv)1 ha 2.469 acre
Volume
1)1 in 3 16.393 cci)1 cc 0.061 in3
2)1 c ft 0.0283 m3 ii)1 m 3
35.330 ft3 1.306 yd3
3)1 mile 3
4.165 km3 iii)1 km 3
0.240 mile3
4)1 acre-ft 0.123 ha-miv)1 ha-m 8.130 acre-ft
Temperature
1) (F 32) 5
____________ ( C 9)
_______
C i) 32 F
9 5
Logarithms
1)log 10
e 0.4343i)ln e
10 2.303
2)log 10
x 0.4342 ln e xii)ln e
x 2.303 log 10 x
*** ***
Appendix III
ROMAN NUMERALS
The Romans who were active in trade and commerce right from the time they learnt to write
devised a system of symbols to indicate numbers. The difference between the Roman and the
Arabic numericals [The one which is used today] is that the Romans did not have a symbol of zero.
The following chart gives the several basic Roman numerals and explains how to determine the
value of other numbers.
I The easiest way to note down a number is to make marks. Thus, I means 1, II means 2, III means 3.
V However, since four strokes seemed too many, the Romans moved on to the symbol V for 5. To obtain 4, they
decided to place I before V—IV.
Thus, placing any small number before a basic Roman number indicates subtraction and placing a number
after the basic Roman number indicates addition. Following this pattern, VI, VII and VIII indicate 6, 7 and
8 respectively.
X The symbol X indicated the number 10. Placing I before X—IX—gives 9 and placing I after X—XI— gives
11. Numbers in the teens, twenties and thirties follow the same pattern. In these sets, the number of X indicates
the number of tens. So, XXVI is 26 and XXXIV is 34.
L L stands for 50. To obtain 40, which is 10 subtracted from 50, it will be written as XL. And thus 60, 70, and
80 are LX, LXX and LXXX.
C C stands for centum, the Latin word for 100. This is used in words like “century” and “cent”. As per the rule,
90 will be written as XC. Like the X’s and L’s, the C’s are tacked on to the beginning of numbers to indicate
how many hundreds are there. CCCLXIX is 369.
D D stands for 500. As can be guessed, CD means 400. So CDXLVIII is 448.
M M is 1,000. This symbol finds a lot of use when Roman numerals are used to indicate dates. For instance, 2007
will be written as MMVII.
— Large numbers were indicated by putting a horizontal line over them. It means to multiply the number by
1000. Thus V means 5000.
*** ***
Appendix IV
GREEK ALPHABET
The Greek alphabet is used as mathematical symbols in most of the scientific studies. The following
table gives the upper- as well as the lower-case letters.
*** ***
Appendix V
METRIC PREFIXES
In practice, several parameters which are very large or very small are required to be referred to.
These are mentioned as multiples of 10. In SI, a number of metric prefixes are used to mention
their largeness or smallness. These are as follows.
Prefix Symbol
million million [Trillion]10 12
tera-T
thousand million [Billion]10 9
giga-G
million10 6
mega-M
thousand10 3
kilo-k
hundred10 2
hecto-h
ten10deca-da
tenth10 –1
deci-d
hundredth10 –2
centi-c
thousandth10 –3
milli-m
millionth10 –6
micro-
thousand millionth10 –9
nano-n
million millionth10 –12
pico-p
thousand million millionth10 –15
femto-f
million million millionth10 –18
atto-a
Since hecto-, deca-, deci-, and centi- are not multiples of powers of three (10 3), these prefixes are
not approved and should be avoided in formal texts.
*** ***
Appendix VI
KEY TO SURVEY MAPS
The Survey of India was set up as the National Survey and Mapping Organization in 1767 under the
Ministry of Science and Technology, Government of India, with headquarters at Dehradun. This
organization publishes various types of maps.
The maps published by the Survey of India are of two series.
A) India and Adjacent Countries (I A C) series
B) International Map of the World (I M W) series
A) The I. A. C. series
In this series, the maps are numbered. Each numbered sheet covers 4 in latitude and 4 in longi-
tude. The numbering starts from NW corner 40 Latitude North and 44 Longitude and increases
from north to south. Sheets covering sea are not numbered.
The Indian Union is located between 8-4' to 37-6' north latitude and 68-7' to 97-25' east
longitude. Sheets covering India are numbered from 39 to 88.
These sheets are known as ‘I/M or million sheets’. The representative scale of these sheets is
1:1,000,000. The map of India showing this division is as shown in Fig. AVI.1.
Each I/M sheet is divided into 16 equal parts and are designated from A to P. Each part there-
fore covers 1 in latitude and 1 in longitude and hence is called ‘Degree Sheet’. For example, the
I/M sheet 47 is divided into 16 degree sheets as shown in Fig. AVI.2.
A E I M
B F J N
47
47D C G K O 47P
D H L P
Fig. AVI.2
An accurate survey is done and each degree sheet is plotted to a scale of 1 : 250 000. All the
details such as rivers, roads, railways, buildings; and administrative boundaries such as taluka, dis-
trict, state; and also forest area, tanks, reservoirs and other important landmarks located in that area
are shown in different colours. Contours at an interval of 100 m are normally shown.
Each degree sheet is further divided into 16 equal parts numbered from 1 to 16. For instance, the
degree sheet 47 I is divided as shown in Fig. AVI.3.
I 1 5 9 13 I
47 47
2 14
2 6 10 14
47 I
3 7 11 15
4 8 12 16
Fig. AVI.3
Here each part covers 0-15' in latitude and 0-15' in longitude and is known as Topo-Sheet and
is plotted to a scale of 1:50,000. All the details shown on a degree sheet are also shown on a topo-
sheet. The contours are at an interval of 20 m.
A topo-sheet to a scale of 1:50000 is further divided into 4 equal parts each 0-7.5' in latitude
and 0-7.5' in longitude and are plotted to a scale of 1:25000. All possible details as shown on the
degree sheet as well as on the topo-sheet are also shown in this map. Contours at an interval of
10 m are shown. If the ground is relatively plain, then the contour interval is maintained as 5 m.
The topo-sheet ‘47A/4’, for instance, is divided as shown in Fig. AVI.4.
A A
47 NW 47 NE
4 4
NW NE
A
47
4
SW SE
Fig. AVI.4
354 APPENDIX VI
Latitude line
Longitude line
Body of map
Left-hand Right-hand
symbol table symbol table
Restricted note
Ref. to
sheet Scale and footnotes
Index of sheets Adm. index
Fig. AVI.5
All maps are provided with the marginal information as shown in Fig. AVI.5 as follows:
1) Title/central heading.
2) District heading.
3) Edition: Year of publication.
4) Sheet number.
5) Legend. [This table explains the various symbols used.]
6) Footnotes: Any information which cannot be shown on the map is given here.
7) Index to sheets: Sometimes, the adjoining sheet has to be referred to. In such case, the adjoining
sheet can be easily located from the ‘Index to sheets’ shown on each of the sheet as shown in
Fig. AVI.6. The central shaded sheet is the one under reference and the adjoining ones are indi-
cated in the index.
8) Administrative index.
Fig. AVI.6
APPENDIX VI 355
Some of the I/M maps, degree sheets and toposheets are classed as RESTRICTED for security
purposes, and cannot be obtained directly. These can be obtained only after permission from the
‘Ministry of Defense, Government of India’.
This I A C series is normally used in India.
B) The I. M. W series
The system is applicable globally (to whole of the world).
The entire area of the world is divided into sheets each covering an area of 4 latitude and 6
longitude.
The geographical positions are defined by two letters ‘N’ for the northen hemisphere and ‘S’ for
the south hemisphere.
The 90 latitudes in each of the hemisphere are divided into horizontal bands, each of 4 Lati-
tudes and are identified by alphabetical letters from ‘A’ onwards. Thus the first band from 0 to 4
as ‘A’; 4 to 8 as ‘B’; 8 to 12 as ‘C’ and so on.
There are two sets of brands ‘A. B. C. -----’ , on either side of 0 Latitude i.e. the equator.
Similarly, the 360 Longitudes are divided into 60 vertical bands each of 6 Longitudes and are num-
bered from 1 to 60. However the numbering starts from 180 Longitude and goes from west to east.
Thus there will be two vertical bands ‘30 and 31’ on either side of 0 Longitude i.e. Greenwich.
The Mercator projection of the Earth is shown in Fig. AVI.7.
Latitude
180° 174° 168° 162° 24° 18° 12° 6° 0° 6° 12° 18° 24° 162° 168° 174° 188° Longitude
28°
24° G
20° F
16° E
N
12° D
NC33
8° C
4° B
0° A Equator
4° A
8° B
12° C
SD28
16° D
S
20° E
24° F
28° G
Greenwich
1 2 3 27 28 29 30 31 32 33 34 58 59 60
Fig. AVI.7
356 APPENDIX VI
A B C D E F
G H I J K L
M N O P Q R
S T U V W X
Fig. AVI.8
Each part will cover an area of 1 Latitude and 1 Longitude and hence it is also termed as
______
ND43
‘Degree Sheet’. The degree sheet in this International System will be identified as H . This
means that the degree sheet is in the north hemisphere, in the ‘D’ band of Latitudes i.e. in the vertical
band ‘43’ i.e. between 12 to 16 of Longitudes, and ‘H’ sub-part.
These are plotted to a scale of 1:250 000 and are identical to the degree sheet stated in the I A C
system except the identification number.
Thus, the degree sheet ‘47 I’ in the IAC system will have the identification as shown in Fig. AVI.9.
20° N
A E I M A B C D E F
B F J N G H I J K L
47 E
C G K O M N O P Q R
D H L P S T U V W X
16° N
72°E 47 76°E 72°E 43 78°E
47 I NE 43
C
Fig. AVI.9
NE43
The index to sheets will be similar to the I A C system. For the degree sheet _____
A , the Index to
sheets will be as in Fig. AVI.10.
Fig. AVI.10
APPENDIX VI 357
*** ***
Appendix VII
MULTIFACETED WATER
Ballast water: Huge tanks from oil-rich countries are used to transport oil. These tanks
are emptied in different countries. After emptying the oil, during their
return journey, the tanks are filled with the available seawater to maintain
the balance of the tanks. These tanks are emptied in the sea at the loading
oil-rich country to refill them again with oil. This sea water which is
used for balancing purpose and emptied in the sea at the oil-rich country
is called ballast water. It is noticed that this ballast water causes some
environmental problems.
Black water: Waste water from toilets is called black water.
Blue water: Water in rivers, flowing ground water, etc. is called blue water.
Gray water: Waste water from bathrooms, bathtubs, showers, wash basins, laundry,
kitchen sinks, etc. is called gray water.
Green water: Water from soil/land utilized for bio mass generation by way of evapo-
transpiration, etc. is called green water.
Hard water: Water from rivers, springs, wells, etc. contain a certain amount of
dissolved salts such as calcium, magnesium salts. Such water that will not
produce good amount of lather with soap is called hard water.
Heavy water: The oxide of heavy hydrogen (Deuterium) ‘D2O’ is called Heavy water. It
is like normal water in appearance. This was discovered by ‘Uray’ in 1932.
It is used in Atomic energy studies.
Magnetized water: It is noticed that after magnetization, some of the physical and chemical
properties of water change. The extent of magnetization depends on the
magnetic poles used, intensity of magnetic field, duration, amount of water
used, etc.
Some believe that this magnetic water is useful for healing some diseases.
Soft water: Water with a very low content of dissolved salts is called soft water. Some
people claim that their thirst is not satisfied by drinking this type of water.
Virtual water: When any item/commodity is imported by a nation, the water required
for growing or producing that item is not accounted for in water balance
e.g., when one kg of wheat is imported by a nation, then 1.2 m3 of water
required for growing that one kg of wheat is indirectly imported. This 1.2
m3 of water is called virtual water.
*** ***
Appendix VIII
A – INDIAN PERIODICALS
B – FOREIGN PERIODICALS
*** ***
Appendix IX
WATER RESOURCES OF INDIA
The Indian Union lies between 8°-4' to 37°-6' North latitude and 68°-7' to 97°-25' East longitude.
India is also known as ‘Bharat, Hindustan, Bharat-varsh, Arya-varta’.
The Indian Union consists of 28 states, six centrally governed union territories and one national
capital territory.
States: Andhra Pradesh, Arunachal Pradesh, Assam, Bihar, Chhatisgarh, Goa, Gujrat, Haryana,
Himachal Pradesh, Jammu and Kashmir, Jharkhand, Karnataka, Kerala, Madhya Pradesh, Maharashtra,
Manipur, Meghalaya, Mizoram, Nagaland, Orissa, Punjab, Rajasthan, Sikkim, Tamil Nadu, Tripura,
Uttar Pradesh, Uttarakhand, West Bengal.
Centrally governed union territories: Andaman and Nicobar, Chandigarh, Dadra and Nagar
Haveli, Daman and Diu, Lakshadweep, Pondicherry.
National capital territory: Delhi.
Physiographically India is divided into seven divisions 1) The Northern mountains,
2) The Great plains, 3) The Central highlands, 4) The Peninsular plateau, 5) The East coast belt,
6) The West coast belt, 7) The islands.
Total land area : 328 million hectares
3.28 million km2
[2.11% of total world area]
Land frontier : 15,200 km
Coastline : 6,100 km
Average annual precipitation : 1140 mm
Maximum average annual precipitation : 10860 mm [At Mawsyram in Assam
near Cherrapunji]
Minimum average annual precipitation : 30 mm [At Churu dist. Bikaner
Rajasthan]
Annual volume of water received over the country : 3700 km3
370 106 hectare m
Evaporation losses [33%] : 1230 km3
123 106 hectare m
Seepage in subsoil [22%] : 800 km3
80 106 hectare m
Run-off in the rivers [45%] : 1670 km3
167 106 hectare m
APPENDIX IX 363
*** ***
Appendix X
INDIA’S NATIONAL WATER POLICY
India adopted a national water policy (NWP) in September 1987 and revised it in 2002, under the
aegis of National Water Resources Council, established in 2000. It recognized water as one of the
crucial elements in the developmental planning. The policy document, inter alia, lays down that
planning and development of this natural resource, i.e., water, needs to be governed by a national
perspective. Resource planning, in case of water, has to be done for hydrological unit such as drain-
age basin as a whole or for a sub-basin. All individual developmental projects and proposals should
be formulated by states and considered within the framework of such an overall plan for a basin or
a sub-basin so that the best possible combinations of options can be made.
The policy document lays down that water should be made available to water-short areas by
transfer from other areas including from one river basin to another, based on a national perspec-
tive, after taking into account the requirements of areas/basins. Transfer of water from one river
to another, especially if it involved inter-state transfer, has always been a sensitive issue amongst
people and the states.
On the issue of equity and social justice, the policy lays down that water allocation in an irriga-
tion system should be done so that disparities in the availability between the head-reach and tail-
end farms and between large and small farms should be obviated by adoption of a rotational water
distributional system and supply of water on a volumetric basis, subject to certain ceilings.
For implementing the above programme, the policy document lays down that appropriate orga-
nizations should be established for planned development and management of a river basin as a
whole.
The policy recognizes water as a prime natural resource, a basic human need and a precious
national asset. The policy has laid down that, in planning and operation of systems, water allocation
priorities should be as follows.
Drinking water
Irrigation
Hydro-power
Navigation
Industrial and other uses
However, these priorities might be modified, if necessary, in particular regions with reference to
area-specific considerations.
The policy stipulates that adequate drinking water facilities should be provided to the entire
population, both in urban and rural areas. Irrigation and multi-purpose projects should invariably
include a drinking water component, wherever there is no alternative source of drinking water.
Drinking water needs of human beings and animals should be provided for all irrigation projects.
APPENDIX X 365
The policy further provides that economic development and activities including agricultural,
industrial and urban development should be planned with due regard to the constraints imposed
by the configuration of water availability. There should be a water zoning of the country and the
economic activities should be guided and regulated in accordance with such zoning.
The implementation of the policy guidelines is overseen by the National Water Resources
Council in which all chief ministers are members and the prime minister is the chairman. The
National Water Board headed by the Secretary, Union Ministry of Water Resources, and the
representatives from all the states assist the Council. The National Water Board has so far final-
ized a number of policy papers on water-related issues for the consideration of the National Water
Resources Council.
*** ***
Appendix XI
IS Number/DOC
Number Title
IS 1191:2003 Hydrometric determination—Vocubulary
| ISO 772:1988
IS 1192:1981 Velocity area Methods for measurement of flow of water in open channels (first
| ISO 748:1979 revision)
IS Number/DOC
Number Title
IS 6062:1971 Method of measurement of flow of water in open channels using standing wave flume-fall
IS 6063:1971 Method of measurement of flow of water in open channels using standing wave flume
IS 6064:1971 Specification for sounding and suspension
IS 6330:1971 Recommendation for liquid flow measurement in open channels by weirs and flumes-
| ISO 3847 end depth method for estimation of flow in rectangular channels with a free overfall
(approximate method)
IS 8389:2003 Code of practice for installation and use of raingauges, recording (second revision)
IS 9108:1979 Liquid flow measurement in open channels using thin plate weirs
| ISO 1438-1
IS 9116:2002 Specification for water stage recorder (float type) (first revision)
IS 9119:1979 Method for flow estimation by jet characteristics (approximate method)
IS 9163(Part 1):1979 Dilution methods of measurement of steady flow Part 1 constant rate injection method
| ISO 9555-1:1973
IS 9922:1981 Guide for selection of method for measuring flow in open channels
| ISO 8363:1980
IS 12752:1989 Guidelines for the selection of flow gauging
| ISO 8368:1980
IS 13083:1991 Liquid flow measurement in open channels—Flat-V
/ ISO 4377:1990
IS 13084:1991 Liquid flow measurement in open channels—Round nose horizontal crest weirs
/ ISO 4374:1990
IS 13371:1992 Code of Practice for calibration (rating) of rotating element current meters instraight
| ISO 3455:1976 open tank
IS 14359:1996 Echo sounders for water depth measurements
| ISO 4366:1979
IS 14371:1996 Measurement of liquid flow in open channels—Parshall and SANIIRI flumes
/ ISO 9826:1982
IS 14573:1998 Liquid flow measurement in open channels—Velocity area methods—
/ ISO 1088:1985 Collection and processing of data for determination of errors in measurement
IS 14574:1998 Measurement of liquid flow in open channels by weirs and flumes- end depth method for
/ ISO 4371:1984 estimation of flow in non rectangular channels with a free overfall (approximate method)
IS 14615 (Part 1):1999 Measurement of fluid flow by means of pressure differential devices Part 1: Orifice plates,
/ ISO 5167-1:1991 nozzles and venturi tubes inserted in circular cross-section conduits running full
IS 14673:1999 Liquid flow measurement in open channels by weirs and flumes—Triangular profile weirs
| ISO 4360:1984
IS 14869:2000 Liquid flow measurement in open channels—Rectangular, trapezoidal and U-shape
/ ISO 4359:1983 flumes
IS 14973:2001 Measurement of fluid flow in closed conduits—Velocity area method using Pilot Static
/ ISO 3966:1997 Tubes
IS 14974:2001 Liquid flow measurement in open channels by weirs and flumes—Rectangular broad-
/ ISO 3846:1989 crested weirs
IS 14975:2001 Measurement of fluid flow in open channels—Stream lind triangular profile weirs
/ ISO 9827:1994
IS 15117:2002 Hydrometric determination—Cable way system for stream gauging
/ ISO 4375:2000
368 APPENDIX XI
IS Number/DOC
Number Title
IS 15118:2002 Measurement of liquid flow in open channel—Water level measuring devices
/ ISO 4373:1995
IS 15119 (Part 1):2002 Measurement of liquid flow in open channels—Part 1 Establishment and operation of a
/ ISO 1100-1:1996 gauging station (superseding IS 2914:1964)
IS 15119 (Part 2):2002 Measurement of liquid flow in open channels—Part 2 Determination of the stage-discharge
/ ISO 1100-2:1998 relation (superseding IS 2914:1964)
IS 15123:2002 Hydrometric determination—Flow measurement in open channels using structures-
/ ISO 4362:1999 Trapezoidal broad-crested weirs (superseding IS/ISO 4362:1992)
IS 15352:2003 Liquid flow measurement in open channels—Position fixing equipment for hydrometric
/ ISO 6420:1984 boats
IS 15353:2003 Liquid flow measurement in open channels by weirs and flumes—V-shaped broad-
/ ISO 8333:1985 crested weirs
IS 15362:2002 Liquid flow measurement in open channels—Flow measurements in open channels using
| 114139:2000 structures-compound gauging structure (Adoption of ISO 14139:2000)
IS 15454:2004 Liquid flow measurement in open channels—Velocity-area method using a restricted
number of verticals (based on ISO/TR 9823:1990)
IS 15527:2004 Measurement of liquid flow in open channels—Measurement in meandering rivers and in
streams with unstable boundaries (based on ISO/TR 9210:1992)
DOC.WRD 1(338) Measurement of liquid flow in open channels—Field measurement of discharge in large
rivers and floods (based on ISO 9825:1994)
DOC.WRD 1(358) Measurement of total discharge in open channels—Stage-fall-discharge relationship
(adoption of ISO 9123:2001)
DOC.WRD 1(389) Measurement of liquid flow in open channels—Method of specifying performance of
hydrometric equipment’ (adoption of ISO
DOC.WRD 1(333) Liquid flow measurement in open channels velocity-area methods—
investigation of total error (based on ISO/TR 7178:1983)
DOC.WRD 1(335) Measurement of liquid flow in open channel—Determination of the wet-line correction
(based on ISO/TR 9209:1989)
DOC.WRD 1(390) Measurement of liquid flow in open channels—Electromagnetic current meters (adoption
of ISO/TR 11974:1987)
DOC.WRD 1(337) Draft Indian Standard Measurement of fluid flow—Evaluation of uncertainties (based on
ISO/TR 5168:1998)
DOC.WRD 1(361) Hydrometric determinations—Unstable channels and ephemeral streams (based on
ISO/TR
IS 4410 (Part 1):1991 Glossary of terms relating to river valley projects: Part 1 Irrigation practice (first revision)
IS 4410 (Part 11/ Glossary of terms relating to river valley projects Part 11 Hydrology Section 1 General
Sec 1):1972 terms
IS 4410 (Part 11/ Glossary of terms relating to river valley projects: Part 11 Hydrology Section 2
Sec 2):1972 Precipitation and run-
IS 4410 (Part 11/ Glossary of terms relating to river valley projects: Part 11 Hydrology: Section 3:
Sec 3):1973 Infiltration and water
IS 4410 (Part 11/ Glossary of terms relating to river valley projects: Part 11 Hydrology Section 4
Sec 4):1973 Hydrographs
IS 4410 (Part 11/ Glossary of terms relating to river valley projects: Part 11 Hydrology Section 5 Floods
Sec 5): 1977
(continued)
APPENDIX XI 369
IS Number/DOC
Number Title
IS 4410 (Part 11/ Glossary of terms relating to river valley projects: Part 11 Hydrology Section 6 Ground
Sec 6):1994 water
IS 14476 (Part 1):1998 Test pumping of water wells—Code of practice Part 1 General
IS 14476 (Part 2):1998 Test pumping of water wells—Code of practice Part 2 Hydrogeological considerations
IS 14476 (Part 3/ Test pumping of water wells—Code of practice Part 3 Pre-test planning Sec 1 General
Sec 1):1998 aspects
IS 14476 (Part 5):1998 Test pumping of water wells—Code of practice Part 5 Pumping test
IS 14476 (Part 8):1998 Test pumping of water wells—Code of practice Part 8 Water level and discharge
measuring devices
DOC.WRD 3(370) Guidelines for artificial recharge to ground water
DOC.WRD 3(456) Manual method for measurement of a ground water level in a well
IS 6939:1992 Methods for determination of evaporation from reservoirs (first revision)
DOC.WRD 10(354) Determination of volume of water & water level in lakes & reservoirs
DOC.WRD 13(378) Guidelines for adopting coefficient of friction (Rugosity coefficient) for design of canals
*** ***
Appendix XII
1. Assamese: i) (pani)
ii) (jal)
iii) (neer)
2. Bengali: i) (bari)
ii) (neer)
iii) (pani)
4. Gujrati: i) (pani)
ii) (jal)
6. Kannada: i) (niiru)
ii) (udak)
iii) (teertha)
7. Kashmiri: i) (poni)
ii) (aab)
9. Maithili: i) (pani)
21. Sanskrit: There are as many as 70 words for water. Each one being originally Sanskrit and
not derived from any other language. These are as follows.
ok% Vãha lye~ Salam
okfj Vãri ÅtZe~ Urjam
lfyye~ Salilam ?k‘re~ Ghrutam
deye~ Kamalam okte~ Vãjam
i;% Payaha vki% Aapaha
dhykye~ Keelãlam lfjye~ Sarilam
ve‘re~ Amrutam lkEc% Sambaha
Hkqoue~ Bhuvanam vU/ke~ Andham
oue~ Vanam dq’ke~ Kusham
dcU/ke~ Kabandham dcqZje~ Karburam
mnde~ Udakam {khje~ Ksheeram
ikFk% Pãthaha rkeje~ Tãmaram
iq"dje~ Pushkaram gªhosje~ Hriveram
loje~ Savaram l´~pye~ Sanchalam
vEHk% Ambhaha lEoje~ Sambaram
v.kZ% Arnaha tMe~ Jadam
rks;e~ Toyam thoue~ Jeevanam
nde~ Dakam dkseye~ Komalam
dEcye~ Kambalam rkeje~ Tãmaram
L;Unue~ Syandanam ?kuhlkle~ Ghanisãsam
vki% Aapaha d‘Iisje~ Krupperarn
lnue~ Sadanam jsikye~ Repãlam
bjk Irã ’kEcje~ Shambaram
vEcq Ambu ?kujl% Ghanarasaha
de~ Kam tyihFke~ Jalapeetham
tM% Jadaha diU/ke~ Kapandham
lksee~ Somam pUnzksjle~ Chandrorasam
APPENDIX XII 373
1½ *fuEueuqljfr tya*
As the waters (of different rivers) enter the ocean, which though full on all sides, remains undis-
turbed.
Spill way is the guaranty of safety for a reservoir of waters which is a pool having plenty of
water.
The rivers having delicious water at the source became unpotable when they merge into the
ocean.
A huge cloud consumes salty water and turns the same into sweet water.
APPENDIX XIII 375
7½ xrksnds lsrqcU/kks
Clouds collect the water of the ocean keeping away its salinity.
Water given at the roots of trees manifests itself high up at the tips of the leaves.
Drops of water on the leaves of lotus are very transient. Same is the case of human life.
On this earth there are three precious things; Water, food and a witty word.
Great rivers go away from the mountains and collectively flow into the sea.
A tree which is watered at the roots shows the fruits up in the branches.
376 APPENDIX XIII
Continuously striken, even by soft water, regions of rugged mountains give away, hard hills are
eroded.
Water works like a medicine, when the food is not digested, but it is invigorating after the food is
digested.
Hence, for the strengthening of appetite, man should drink a small quantity of water, a number of
times.
For the mountains water works like an old age. (As it decreases human energy.)
The clearing nut, when dropped into turbid water, removes the turbidity of water and also vanishes
itself.
27½ vYirks;ÜoyRdqEHk%
28½ vk;q:nŒpu;U=kksiefefr
Avarta (cloud) is confined to a particular locality. Under Samvarta cloud water is everywhere,
under Puskara cloud water is scarce and under Drona cloud the water is abundant on earth.
Gravity is the cause for falling liquids and solids. It is invisible and is inferred by the falling
motion.
Foundation should either be levelled by the instrument, or the foundation should be filled with
water (and the level should be checked according to water level).
Place copper plates in an earthen pot; cover it with powdered coal and moistened saw dust spread
zinc powder and cover it with mercury. Due to chemical interaction positive and negative electric-
ity is produced. Due to this, water is decomposed into oxygen and hydrogen.
378 APPENDIX XIII
*** ***
Appendix XIV
WISDOM OF WATER
Scientists, philosophers and statesmen, have realized the importance of ‘water’ and have quoted it exten-
sively in their speeches and writings.
A few of them have been quoted here.
If we could ever competitively, at a cheap rate, get fresh water from saltwater, ... (this) would be in
the long-range interests of humanity which could really dwarf any other scientific accomplishments.
—John F. Kennedy
Water and air, the two essential fluids on which all life depends, have become global garbage cans.
—Jacques Cousteau (1910–1997)
There is nothing softer and weaker than water, and yet there is nothing better for attacking hard and
strong things.
For this reason there is no substitute for it.
—Lao-Tzu (C.B.C. 550)
Water that has been begged for does not quench the thirst.
—Geographical Origin Uganda
Water is H2O, hydrogen two parts, oxygen one, but there is also a third thing that makes water and
nobody knows what that is.
—D. H. Lawrence (1885–1930), Pansies, 1929
Water is life’s matter and matrix, mother and medium. There is no life without water.
—Albert Szent-Gyorgyi quotes (Hungarian Biologist,
1937 Nobel Prize for Medicine 1893–1986)
All streams flow into the sea, yet the sea is never full. To the place the streams come from, there they
return again.
—Ecclesiastes 1:7 from New International Version of The Bible
Nothing is softer or more flexible than water, yet nothing can resist it
—Lao Tzu quotes [Chinese Taoist philosopher. Founder of Taoism,
wrote “Tao Te Ching” (also “The Book of the Way”) 600 BC–531 BC]
Water, like religion and ideology, has the power to move millions of people. Since the very birth of
human civilization, people have moved to settle close to it. People move when there is too little of it.
People move when there is too much of it. People journey down it. People write, sing and dance about
it. People fight over it. And all people, everywhere and every day, need it.
—Mikhail Gorbachev, President of Green Cross International
quoted in Peter Swanson’s Water; The Drop of Life, 2001.
Water has no taste, no color, no odor; it cannot be defined, art relished while ever mysterious. Not nec-
essary to life, but rather life itself. It fills us with a gratification that exceeds the delight of the senses.
—Antoine De Saint-Exupery (1900–1944), Wind, Sand, and Stars, 1939
Water is the one substance from which the earth can conceal nothing; it sucks out its innermost secrets
and brings them to our very lips.
—Jean Giraudoux (1882–1944), The Madwomen of Chaillot, 1946
Throughout the history of literature, the guy who poisons the well has been the worst of all
villains.
—Author unknown
Wetlands have a poor public image…. Yet they are among the earth’s greatest natural assets…
mankind’s waterlogged wealth.
—Edward Maltby, Waterlogged Wealth, 1986
Life originated in the sea, and about eighty percent of it is still there.
—Isaac Asimov, Isaac Asimov’s Book of Science and Nature Quotations, 1988
APPENDIX XIV 381
The estuary is the point where man, the sea—his immemorial ally and adversary—and the land meet
and challenge each other.
—U.S. Department of the Interior, National
Estuarine Pollution Study, November 1969
In every glass of water we drink, some of the water has already passed through fishes, trees, bacteria,
worms in the soil, and many other organisms, including people…. Living systems cleanse water and
make it fit, among other things, for human consumption.
—Elliot A. Norse, in R.J. Hoage, ed., Animal Extinctions, 1985
Between earth and earth’s atmosphere, the amount of water remains constant; there is never a drop
more, never a drop less. This is a story of circular infinity, of a planet birthing itself.
—Linda Hogan, Northern Lights, Autumn 1990
The stone in the water knows nothing of the hill which lies parched in the sun.
—African Proverb
The highest good is like water. Water gives life to the ten thousand things and does not strive. It flows in
places men reject and so is like the Tao.
—Excerpt from the Tao Te Ching, chapter 8
It is a fascinating and provocative thought that a body of water deserves to be considered as an organ-
ism in its own right.
—Lyall Watson, Supernature
Water links us to our neighbor in a way more profound and complex than any other.
—John Thorson
We must begin thinking like a river if we are to leave a legacy of beauty and life for future generations.
—David Brower quoted by E-Wire, 7 Apr 2000
382 APPENDIX XIV
Rain is grace; rain is the sky condescending to the earth; without rain, there would be no life.
—John Updike, Self-Consciousness: Memoirs, 1989
Water helped ancient man learn those first lessons about the rights of others and responsibility to a
larger society… . It became part of the moral and mental legacy parents passed on to their children.
—M. Meyer, Water in the Hispanic Southwest
Truths are first clouds; then rain, then harvest and food.
—Henry Ward Beecher
Any river is really the summation of the whole valley. To think of it as nothing but water is to ignore the
greater part.
—Hal Borland, This Hill, This Valley
Water is the basis of life and the blue arteries of the earth! Everything in the non-marine environment
depends on freshwater to survive.
—Sandra Postel, Sandra Postel, Global Water Policy Project,
Grist Magazine 26 Apr 04
Clean water is not an expenditure of Federal funds; clean water is an investment in the future of our
country.
—Bud Shuster, U.S. Representative,
quoted in The Washington, Post, 1 Sep 1987
No one can see their reflection in running water. It is only in still water that we can see.
—Taoist proverb
As water runs towards the sword, money towards the rich man’s hand.
—Proverb
We used to think that energy and water would be the critical issues for the next century. Now we think
water will be the critical issue.
—Mostafa Tolba of Egypt,
former heard of the United Nations Environment Program
To trace the history of a river… is to trace the history of the soul, the history of the mind descending and
arising in the body.
—Gretel Ehrlich
Water is fundamental for life and health. The human right to water is indispensable for leading a healthy
life in human dignity. It is a pre-requisite to the realization of all other human rights.
—The United Nations Committee on Economic, Cultural and
Social Rights, Environment News Service, 27 Nov 02
Rivers are roads which move, and which carry us whither we desire to go
—By Blaise Pascal–famous water sayings
Water is sometimes sharp and sometimes strong, sometimes acid and sometimes bitter, sometimes sweet
and sometimes thick or thin, sometimes it is seen bringing hurt or pestilence, sometime health-giving,
sometimes poisonous. It suffers change into as many natures as are the different places through which it
passes. And as the mirror changes with the colour of its subject, so it alters with the nature of the place,
becoming noisome, laxative, astringent, sulfurous, salty, incarnadined, mournful, raging, angry, red,
yellow, green, black, blue, greasy, fat or slim. Sometimes it starts a conflagration, sometimes it extinguishes
one; is warm and is cold, carries away or sets down, hollows out or builds up, tears or establishes, fills or
empties, raises itself or burrows down, speeds or is still; is the cause at times of life or death, or increase or
privation, nourishes at times and at others does the contrary; at times has a tang, at times is without savor,
sometimes submerging the valleys with great floods. In time and with water, everything changes
—Leonardo da Vinci
384 APPENDIX XIV
Too much or too little of water and its quality has always been responsible for the largest number
of diseases and deaths.
The movement of heavenly bodies despite their distances from earth have presented fewer difficul-
ties to me than the movement of water which is within my reach
—Galileo
Rain is a blessing when it falls gently on parched fields, turning the earth green, causing the birds to
sing.
—Donald Worster, Meeting the Expectation of Land, 1984
The heart of the wise man lies quiet like limpid water .
—Geographical Origin Cameroon
It is really important to solve the problem of rational utilization and distribution of water supplies.
I dare say, the shortage of fresh water is the major ecological problem of this moment.
—Mikhail Gorbachev, President of Green Cross International
quoted in Peter Swanson’s Water: The Drop of Life, 2001
The good rain, like the bad preacher, does not know when to leave off.
—Ralph Waldo Emerson
Water, taken in moderation, cannot hurt anybody.
—Mark Twain
APPENDIX XIV 385
Water links us to our neighbor in a way more profound and complex than any other.
—John Thorson
Aquifer: a mysterious, magical and poorly defined area beneath the surface of the earth that either
yields or withhold vast or lesser quantities of standing/flowing water, the quantity and/or quality of
which is dependent on who is describing it or how much money may be at stake.
—R. Radden, “Watershed Resources”, Jan. 2002
Irrigation of the land with seawater desalinated by fusion power is ancient It’s called rain
—Michael McClary
We forget that the water cycle and the life cycle are one.
—Jacques Cousteau
You should respect each other and refrain from disputes; you should not, like Water and oil, repel each
other, but should, like milk and water, mingle together.
—Buddha quotes, (Hindu Prince Gautams ***
the founder of Buddhism, 563–483 B.C.)
The wise adapt themselves to circumstances, as water moulds itself to the pitcher
—Anonymous (Chinese Proverb) quotes
Like swift water, an active mind never stagnates.
The wise man of Miletus thus declared the first of things is water
—J.S. Blackie, 1877
386 APPENDIX XIV
Nothing on earth is so weak and yielding as water, but for breaking down the firm and strong it has no
equal.
—Lao-Tsze
*** ***
Bibliography
Chitale, M. A. 2000. A Blue Revolution, Mumbai: Bhavans Book University.
Chow, V. T., D. R. Maidment, and L. W. Mays, 1973. Applied Hydrology, New York: McGraw-Hill.
Chow, V. T. 1964. Handbook of Applied Hydrology, New York: McGraw-Hill.
Kulandaiswamy, V. C. 1971. Lecture Notes on Hydrology, Guindy: ISTE Summer School.
Linsley, R. K., M. A. Kohler, and J. L. H. Paulhus, 1988. Hydrology for Engineers, New York:
McGraw-Hill.
Manning, J. C. 1987. Applied Principles of Hydrology, Princeton: Merrill Publishing Company.
Mutreja, K. N. 1982. Applied Hydrology, New Delhi: Tata McGraw-Hill.
National Water Policy. 2002. Ministry of Water Resources, Government of India.
Nemec, J. 1973. Engineering Hydrology, New Delhi: Tata McGraw-Hill.
Patra, K. C. 2002. Hydrology and Water Resources Engineering, Delhi: Narosa Publishing House.
Ponce, V. M. 1989. Engineering Hydrology, Englewood Cliffs: Prentice Hall.
Raghunath, H. M. 1983. Groundwater, New Delhi: Wiley Eastern Limited.
Raghunath, H. M. 1985. Hydrology, New Delhi: Wiley Eastern Limited.
Reddy, P. J. 1992. A Textbook of Hydrology, New Delhi: Laxmi Publications.
Sharma, R. K. 1983. A Textbook of Hydrology and Water Resources, New Delhi: Dhanpat Rai & Sons.
Subramanya, K. 1994. Engineering Hydrology, 2nd ed., New Delhi: Tata McGraw-Hill.
Todd, D. K. 1959. Ground Water Hydrology, Tokyo: Toppan Company Ltd.
Varshney, R. S. 1986. Engineering Hydrology, Roorkee: Nem Chand & Bros.
Viessman, W. Jr., G. L. Lewis, and J. W. Knapp, 1989. Introduction to Hydrology, New York: Harper & Row,
Publishers Inc.
Wanielista, M., R. Kersten, and R. Eaglin, 1997. Hydrology, Hoboken: John Wiley & Sons Inc.
Ward, R. C. and M. Robinson, 1989. Principles of Hydrology, New York: McGraw-Hill.
Water Resources of India. 1987. Central Water Commission, Government of India.
Wisler, C. O. and E. F. Brater, 1959. Hydrology, New York: John Wiley & Sons.
Index
A Artificial neural network, 11
T h2 Artificial precipitation, 57
K dh , 338
A∫ ∫
= − Atmosphere, 7, 16, 17
0 h1
h density and pressure of, 17
∆t, 289 existence of, 17
selection of, 289 its constituents, 16
⌽ index, 116 carbon dioxide, 17
Abstractions, 98 humidity, 17
depression storage, 99 inert gases, 16
evaporation, 99 nitrogen, 16
infiltration, 99 oxygen, 16
interception, 99 ozone, 17
Actual observations, 33 vertical structure of, 17
Aeration zone, 316 Automatic rain gauges, 60
Air, 19, 21, 53 Automatic water level recorder, 240
mechanism to cool, 53 Average daily temperature, 21
current, 19
front, 21 B
mass, 21 Bank storage, 185
Air-line corrections, 255 Barlow’s method, 196
Ali Nawaz Jung Bahadur formula, 216 Barometric pressure, 28
Altitude, 28 Base flow graph, 133
Anemometer, 19 extension of, 133
Annual flood series, 219 Base flow, 128
Antecedent precipitation index, 186, 211 Base period, 129
Anticyclone, 22 Basic flow, See Base flow
Apparent velocity, 321 Basin, 210
Aquiclude, 319 shape of, 210
Aquifer, 318, 320 size of, 210
parameters of, 320 Binnie’s table, 192
Aquifuge, 319 Blizzard, 56
Aquitard, 329 Broad-crested weir, 268
Area of influence, 330
Area–slope method, 244, 245 C
demerits of, 245 California formula, 219
merits of, 245 Capillary fringe, See Capillary zone
Area–velocity method, 246 Capillary zone, 317
area of flow, 246 Catchment, 188
velocity of flow, 246 geographical characteristics of, 188
Arithmetic average method, 66 meteorological characteristics of, 188
Arithmetic mean method, 72 shape of, 188
Artesian well, 319 size of, 188
390 INDEX
Q
Quick return flow, 192 S
S hydrograph, 143
R Safe depression head, 337
Radio-isotope method, 336 Salt, 262
Radius of influence, 330 properties of, 261
Rain gauge, 57–59, 69 Salt titration method, 261
selecting site for, 58 Salt water, 4
INDEX 395