Engineering Hydrology Gradeup
Engineering Hydrology Gradeup
co
1
www.gradeup.co
HYDROLOGY
1 PRECIPITATION
1. INTRODUCTION
Ridge: It is a line which divides one catchment are from its neighbouring catchment
area. It is also called water divide or divide.
2
www.gradeup.co
3
www.gradeup.co
2. PRECIPITATION
Precipitation denotes all form of water that reach the earth from the atmosphere.
Precipitation is one of the three major phases in hydrologic cycle and one of the major items
of inflow in the hydrologic equation. It includes rain, snow, hail, dew and frost.
2.1 Types and Causes
The following are the different types of precipitation.
(i) Rain: It is the most dominant from of precipitation in Indian. It denotes water
droplets of size varying 0.5 to 6 mm. On the basis of which intensity of rainfall is
classified as follows
Intensity (mm/hr)
Note:
• In India, rainfall data is collected every day at 8.30 am and if rainfall on a particular
day is more than 2.5 mm then that day is known as rainy day.
• Average rain fall per year in India is 119 cm ≃ 120 cm.
(ii) Drizzle: It denotes fine droplets of water whose size is less than 0.5 mm and
intensity is less than 1 mm/hr
(iii) Snow: These are ice crystals having the density of 0.1 gm/cc
(iv) Glaze: When droplet of water comes in contact with cold ground surface, then
water is converted into ice, which is called as glaze or freezing rain.
(v) Sleet: These are frozen rain drops of transparent nature which forms when rain
falls through subfreezing temperature.
(vi) Hail: These are lumps of ice whose size is greater than 8 mm.
4
www.gradeup.co
• Frontal Precipitation: Under certain favourable conditions when a warm air mass
and cold air mass meet, the warmer air mass is lifted over the colder one with the
formation of a front. The ascending warmer air cools adiabatically with the
consequent formation of clouds and precipitation.
• Cyclonic Precipitation: This type of precipitation occurs due to the pressure
differences. A cyclone is a large low-pressure region with circular wind motion. The
rainfall will normally be heavy in the entire region occupied by the cyclone.
The amount of rain collected by a rain gauge in the last 24 hrs is called daily rainfall and the
amount collected in one year is called as annual rainfall. The average value of annual rainfall
for a period of 35 years is called average annual rainfall.
5
www.gradeup.co
4. INDEX OF WETNESS
The actual precipitation for a given year may differ from the long-time mean. To bring out
this important relationship. We often express the actual precipitation for a given year as a
ratio to the mean. Such a factor is called the index of wetness. This factor is used to find
deficiency in rainfall and is given as follows:
𝑅𝑎 𝑖𝑛𝑓 𝑎 𝑙𝑙 𝑖𝑛 𝑎 𝑦𝑒𝑎𝑟
𝐼𝑛𝑑𝑒𝑥 𝑜𝑓 𝑤𝑒𝑡𝑛𝑒𝑠𝑠 = × 100
𝐴𝑣𝑒𝑟𝑎𝑔𝑒 𝑎𝑛𝑛𝑢𝑎𝑙 𝑟𝑎 𝑖𝑛𝑓 𝑎 𝑙𝑙
• If index of wetness is 100%, it indicates a normal year. If it is greater than 100%, it is
called as good year and if it is less than 100%, it is called bad year.
5. DROUGHT
6
www.gradeup.co
6. MEASUREMENT OF RAINFALL
Rainfall is expressed in terms of depth to which water would stand on an area if all the
rainfall is collected on it. Rainfall is measured by an instrument called Rain gauge which is
also known as pluviometer, ombrometor, hyetometer.
Following are the requirements which has to be satisfied prior to installation of rain gauge.
(a) Rain gauge must be surrounded by an open fenced area of at least 5.5 m × 5.5 m
(b) Rain gauge must be installed at a distance of at least 30 m or twice the height of building
or obstruction.
(c) Rain gauge must be installed on a level surface which is free from undulation.
A rain gauge basically consists of cylindrical assembly which is kept in open in order to collect
the rain water. There are basically two types of rain gauges:
(i) Non Recording rain gauge
In India, most commonly used non recording rain gauge is Symon’s rain gauge which has a
collecting diameter of 127 mm. Now a days IMD (Indian meteorological department) is using
fibre glass reinforced polyester type rain gauge. Which comes in two variants having
collecting area of 100 cm2.
(ii) Recording type rain gauge
This type of rain gauge automatically records the rate of precipitation and the time of its
occurrence. It also produces a continuous variation of rainfall with time. Graphically which
can be represented as follows.
In India natural siphon or float type is used as the standard recording type rain gauge.
Latest improvements in rain gauge technology has led to the development of following rain
gauges.
(ii) Telemetering Rain Gauge
It is basically a recording type rain gauge which contains electronic equipment to transfer
rainfall data to the base relation. As such it can be used in far remote and inaccessible
locations.
7
www.gradeup.co
7. NETWORK DENSITY
𝐶𝑣 2
𝑁=( )
∈
Here, Cv = coefficient of variation.
𝜎
𝐶𝑣 = × 100
𝑥̄
∑𝑛
𝑖=1 𝑥𝑖 ∑𝑛
𝑖=1(𝑥𝑖 –𝑥̄ )
2
Here, 𝑥̄ = 𝑎𝑛𝑑 𝜎 = √
𝑛 𝑛−1
8
www.gradeup.co
Example: The average Annual rainfall of 5 rain gauge stations are 90, 55, 42, 40 and
50 cm. It the error in the estimation of rainfall should not exceed 10%. How many
Solution: -
The mean rainfall is obtained as
90 + 55 + 42 + 40 + 50
X or Pm =
5
X = 55.4 cm
Now,
(90 − 55.4)2 + (55 − 55.4)2 +
(42 − 55.4)2 + (40 − 55.4)2 + (50 − 55.4)2
2 =
5 −1
P − Pm
2 = i
m−1
σ2 = 410.8
σ = 20.27
The coefficient of variation is calculated as:
CV =
X
20.27
CV = = 0.366
55.4
Now, optimum no. of rain gauges, N
2
C
N= V
2
0.366
N= ( = 10%)
0.10
N = 13.396 ~ 14
Thus, additional rain gauges required = 14 – 5
= ‘9’
Example: A catchment area has seven rain gauge stations. In a year the annual rainfall
recorded by the gauges are as fallows.
Station P Q R S T U V
Rainfall (cm) 120 14 110 130 150 116 45
(a) Determine the standards error in the estimation of mean rainfall in the existing set
of rain gauges.
9
www.gradeup.co
(b) for a 5% error in the estimation of mean rainfall calculate the minimum number of
additional rain gauge stations to be establish in the catchment.
Solution: -
For the given data,
No. of rain gauges, m = 7
Mean annual rainfall
120 + 140 + 110 + 130 + 150 + 116 + 145
P =
7
P = 130.142
Now, standard deviation,
m
(P
1
i −P )2
m −1 =
m −1
(120 − 130.14)2 + (140 − 130.14)2
+ (130 − 13014)2 + (150 − 130.14)2
+ (116 − 130.14)2 + (145 − 130.14)2
m −1 =
7 −1
8. CONSISTENCY OF RAINFALL
Whenever there is change in condition at rain gauge prevailing, there will be inconsistency in
rainfall data available for that period. This inconsistency may be due to the following reasons.
10
www.gradeup.co
If on a particular day of a year rainfall data at some station ‘x’ could not be recorded (due to
defective rain gauge or otherwise) and it is required to find an approximate value of this
missing data, following methods are adopted.
Few stations close to problematic station ‘x’ are selected and rainfall value at these stations is
noted down as P1, O2, ….., Pn – 1’ Pn. Px is the missing rainfall data whose approximate value
has to be obtained.
Let N1, N2, N3, … and Nx be the normal precipitation values for station ‘1 to m’ and ‘x’ Normal
precipitation is average value of rainfall for a day, month or a year on the basis of record of
last 30 years.
Case 1: when N1, N2, … Nm differs from Nx by less than 10% the value of Px is given as
follows.
𝑃1 + 𝑃2 + 𝑃3 + 𝑃3 +. . . +𝑃𝑚
𝑃𝑥 =
𝑚
Case 2: when one or more of N1, N2,…Nm, differs from Nx by 10% or more, the value of Px is
calculated by:
𝑁𝑥 𝑃1 𝑃2 𝑃𝑚
𝑃𝑥 = [ + +. . . . . . . + ]
𝑚 𝑁1 𝑁2 𝑁𝑚
11
www.gradeup.co
Solution: -
As the normal precipitation of other stations, A, B and C are not within 10% of normal
precipitation at station x
As,
When a catchment or any other area contains more than one rain gauge station, then
average rainfall for the area can be obtained using the following method.
10.1. Arithmetical Mean/Average Method
This method is suitable if rainfall is uniformly distributed and area is not very large.
𝑃1 + 𝑃2 + 𝑃3 +. . . . . +𝑃𝑛
𝑃𝑎𝑣𝑔 =
𝑛
This method does not give very good result and hence is not used very frequently. Any
station outside the area of consideration is not taken into account in this method.
10.2. Thiessen polygon/mean Method
In this method the rainfall recorded at each station is given a weightage on the basis of
an area closest to the station, that is why this method is also known as weightage
area method.
The procedure of determining the weighing are is as follows:
Step 1: The catchment area is drawn to scale and the position of station is also drawn.
Step 2: All the stations are joined together to form a network of triangles.
12
www.gradeup.co
Step 3: Perpendicular bisectors for each of the sides of the triangle are drawn.
These bisectors form a polygon around each station. Stations are then assigned the
area by which they are surrounded.
𝑃1 𝐴1 + 𝑃2 𝐴2 +. . . +𝑃𝑛 𝐴𝑛
𝑃𝑎𝑣𝑔 =
𝐴1 + 𝐴3 +. . . +𝐴𝑛
∑𝑛
𝑖=1 𝑃𝑖 𝐴𝑖
⇒𝑃𝑎𝑣𝑔 =
𝐴
The ratio Ai/A is called the weightage factor for each station.
This method of finding average rain fall is suitable when area is large and rainfall is
non-uniformly distributed. This method is superior to arithmetical mean method.
10.3. Isohyetal Method
An isohyet is a line joining all the points having same value of rainfall and isohyetal
maps are the one which shows contours of equal rainfall magnitude.
In isohyetal method, it is assumed that the precipitation in areas between isohyetal
lines is equal to the mean of the precipitation of at isohyetal lines.
Mathematically, various following cases are possible.
Case 1:
𝑃1 + 𝑃2 𝑃 + 𝑃2 𝑃 + 𝑃𝑛
( ) 𝐴1 + ( 2 ) 𝐴2 +. . . + ( 𝑛−1 ) 𝐴𝑛−1
𝑃𝑎𝑣𝑔 = 2 2 2
𝐴1 + 𝐴2 + 𝐴3 +. . . +𝐴𝑛−1
𝑃0 + 𝑃1 𝑃 + 𝑃2 𝑃 + 𝑃3 𝑃 + 𝑃𝑛 𝑃 + 𝑃𝑛+1
( ) 𝐴0 + ( 1 ) 𝐴1 + ( 2 ) 𝐴3 + ( 𝑛−1 ) 𝐴𝑛−1 + ( 𝑛 )𝐴
𝑃𝑎𝑣𝑔 = 2 2 2 2 2
𝐴0 + 𝐴1 + 𝐴2 . . . . . +𝐴𝑛−1 + 𝐴𝑛
13
www.gradeup.co
Case 3:
𝑃1 + 𝑃2 𝑃 + 𝑃3 𝑃 + 𝑃𝑛
( ) 𝐴1 + ( 2 ) 𝐴2 +. . . + ( 𝑛−1 ) 𝐴𝑛−1 + 𝑃𝑛+1 𝐴𝑛
𝑃𝑎𝑣𝑔 = 2 2 2
𝐴1 + 𝐴2 + 𝐴3 +. . . +𝐴𝑛−1 + 𝐴𝑛
(i) This method gives better result only if there are larger number of stations.
(ii) Inter isohyetal area is a function of rainfall value which will change with any change
of rainfall pattern
(iii) This method requires use of technical and skilled man power.
Note: Among all the methods discussed above orographically derived isohyetal method
is the best method of finding average rainfall.
Example: Five rain gauges stations namely A, B, C, D and E are located on a circular
shape basin of diameter 20 km as shown in the fig. given below compute the mean
annual rainfall at stations A, B, C, D and E are 100 cm, 90 cm, 110 cm, 120 cm and 80
cm respectively.
Solution: -
Diameter of basin = 20 km. we know that in Thiessen polygon method. The polygons
are obtained by joining the perpendicular bisectors of triangles formed, when various
rain gauge stations are joined.
14
www.gradeup.co
(20)
2
− 100
= 4 = 53.54 km2
4
Mean areal rainfall:-
PA A A + PB AB + PC AC + PD AD + PE AE
=
A A + AB + A C + AD + AE
= 100 cm
mean area rainfall
= 100 cm
Example: The Isohyets for annual rainfall over a catchment basin were drawn. The
areas of strips between the Isohyets are given determine the Average precipitation over
the basin
75-85 1500
85-95 2500
95-105 2800
105-115 1200
115-135 800
135-155 600
15
www.gradeup.co
Solution:
75 + 85 85 + 95
1500 + 2500
2 2
95 + 105 105 + 115
+ 2800 + 1200
2 2
115 + 135 135 + 155
+ 800 + 600
2 2
Pavg =
1500 + 2500 + 2800 + 1200 + 800 + 600
(1500 80) + (2500 90) + (2800 100)
+ (1200 110) + (800 125) + (600 145)
PAvg =
9400
PAvg = 100.425 cm Ans.
Example: Calculate the average depth of rainfall for the catchment with an area of 200
km2 using (i) Arithmetic method (ii) Thiessen polygon method.
Station (i) Rainfall (mm) Pi Thiessen polygon area (km2) Description inside (I) outside (0)
1 45 7.5 0
2 20 10.0 0
3 30 45.0 I
4 45 30.0 I
5 50 16.0 0
6 35 45.0 I
7 45 40.0 I
8 40 6.5 0
16
www.gradeup.co
Solution:
In this method we will consider rain gauges which are inside the catchment only.
P1 + P2 + P3 + P4 . . . Pn
P =
n
30 + 45 + 35 + 45
P =
4
145
P = = 36.25 mm
4
volume of rainfall = A P
36.25 3
= 200 × 1000 × 1000 × m
100
= 72.5 × 105 m3
****
17
www.gradeup.co
18
www.gradeup.co
1
www.gradeup.co
ENGINEERING HYDROLOGY
1. INTRODUCTION
A stream can be defined as a flow channel into which the surface runoff from a specified basin
drains. The science and practice of water measurement is called ‘hygrometry’. Among different
parts of hydrological cycle, stream flow measurement is the only one that can be done
accurately.
Stage is defined as water surface evaluation measured fixed value. The process of
measurements of stage is called gauging. Following methods may be used in order to find the
value of a G.
2.1. Staff Gauge
The simplest of stage measurements are made by noting the evaluation of the water
surface in contact with a fixed graduated staff.
It is fixed rigidly to a structure such as an abutment, pier, wall etc. The staff may be
vertical or inclined with clearly and accurately graduated permanent markings.
It is a gauge used to measure the water surface elevation from above the surface such
as from a bridge or similar structure. In this a weight is lowered by a reel to touch the
water surface.
A mechanical counter measures the rotation of the wheel which is proportional to the
length of the wire paid out. The operating range of this kind of gauge is about 25 m.
2
www.gradeup.co
In this gauge, compressed air or gas is made to bleed out at a very small rate through
an outlet placed out at a very small rate through an outlet placed at bottom of the river.
A pressure gauge measures the gas pressure which in turn is equal to the water column
above the outlet. A small change in the water surface elevation as felt as a change in
pressure from the present value at the pressure gauge and this in turn is adjusted by a
servo-mechanism to bring the gas to bleed at the original rate under the new head. The
pressure gauge reads the new water depth which is transmitted to a recorder.
pgas
H=
r
3
www.gradeup.co
3. MEASUREMENT OF VELOCITY
The measurement of velocity is an important aspect of many direct stream flow measurement
techniques.
Following methods may be used in order to find the velocity of a river or a stream.
3.1. Float
In this method, a very simple float device is used which flows along the river surface.
L
Mathematically, V =
t1
Measurement of float is affected by wind or other such obstructions and hence this does
not give very reliable results. Some improvements in the conventional flow has led to the
development of consister float and rod float
V = aNs + b
4
www.gradeup.co
In order to reduce the drag force these are stream lined in shape. Value of this sounding
weight is given as
W = 50 v y
Here, w = weight in ‘N’
V0.2y + V0.8y
Vavg =
2
Case 2: For shallow k/a one point formula
Vavg = V0.6y
Vavg = kvs
5
www.gradeup.co
Example In a triangular channel, the top with and depth of flow were 2m and 0.9m
respectively. Velocity measurement on the centre line at 18 cm and 72 cm below the
free water surface indicated velocities of 0.6 m/s and 0.4 m/s respectively. The discharge
in the channel is
Solution.
A = 1/2 × 2 × 0.9 = 0.9m
VAvg = (V0.2y + V0.8y) / 2
= (0.4 +0.6) / 2 = 0.5
∵ given depths,
18cm is 0.2y depth i.e = 0.2 × 0.9 = 0.18 m = 18 cm
For 72 cm is 0.8y depth i.e = 0.8 × 0.9 = 0.72m = 72 cm
Q = A × V = 0.9 × 0.25
= 0.45 m3/s
4. MEASUREMENT OF DISCHARGE
6
www.gradeup.co
n
Q= Qi
i=1
Q1 = A1V1
Q2 = A2 V2
Qn = An Vn
Here, A1 = w1y1
A2 = w2y2
An = wny2
Here,
2
w2
w1 + 2
w1 =
2w1
w3 + w4
w3 =
2
2
wn−1
wn = 2
Similarly, wn−1 =
2wn
This method is also known as the chemical method depends upon the continuity principle
Mathematically,
7
www.gradeup.co
Example. A 300 g/l solution of common salt was discharged into stream at the rate of
35 l/s at a downstream section, the solution was completely mixed and an equilibrium
concentration of 70 ppm Was reached. If the background concentration was 15 ppm. The
C0 = 15 × 10–6 kg/l
C2 = 70 × 10–6 kg/l
on solving,
Q = 191 m3/s
This is essentially an area velocity with the average velocity being measured by using
ultrasonic signals. In this method two transducers are provided on both the sides of river
or stream
l
t1 = → (i)
C + v cos
l
t2 = → (ii)
C + v cos
l 1 1
v= −
2 cos t1 t
Q = AV
8
www.gradeup.co
Such as the ganga by the current meter method is very time consuming even when the
flow is low or moderate. It is particularly suitable for wild streams whose width is more.
It this method, a boat moves from one bank to the other bank of the river in which current
meter is installed.
The boat is also provided with an angle indicator which records the direction of resultant
velocity.
0 + y1 y1
Q1 = A1Vf = Vb t1 vf = vb0vf t1
2 2
y1
= Vr2 sin cos t1
2
Similarly,
y + y2
Q2 = Vr2 sin t2 1
2
y + y3
Q3 = Vr2 sin t3 2
2
This method finding discharge makes use of law of conservation of energy for fluid
(Bernoulli’s equation) In this method two section are selected such that cross section
9
www.gradeup.co
v12 V2
Z1 + y1 + = Z2 + y2 + 2 + h2
2g 2g
V2 − V22
h2 = (z1 + y1 ) − (z2 + y2 ) + 1
2g
hf = frictional
Here, h2 = hf + hc
hc = eddies
V2 − V22
hf = (h1 − h2 ) + 1 − he
2g
V12 − V22
he = K e
2g
hf hf
-S= ⇒ Q=k
L L
Here, Q is a function of hf, whereas the value of hf depends on Q and hence in order to
find discharge following method is adopted.
(i) To start with, assume V1 = V2, which gives hf = h1 – h2
hf
(ii) Find discharge using Q = k
L
10
www.gradeup.co
Q Q
V= 1 and V = 2
A1 A2
1
K for section (1)→ K1 = A1 R12/3
2
1
K for section (2)→ K2 = A2 R 2/3
2
2
K= k1k 2
Example. In a river carrying a discharge of 142 m3/s, the stage at a station A was 3.6m
and the water surface slope was 1 in 3000. If during a flood, the stage at A was 3.6m
and the water surface slope was 1 in 6000, what was the flood discharge?
Solution.
Q s
s1 s1 142 1 / 3000
Q1/Q2 = = =
s2 s2 Q2 1 / 6000
Q2 = 200 m3/s
Example: During a flood flow the depth of water in a 10 m wide rectangular channel was
found to be 3.0 m and 2.9 m at (two section 200 m apart. The drop in the water-surface
elevation was found to be 0.12 m. Assuming Manning’s coefficient to be 0.025, estimate
the flood discharge through the channel. Assume (h e = 0).
Solution:
Using suffixes 1 and 2 to denote the upstream and downstream section respectively, the
cross-sectional properties are calculated as follows:
11
www.gradeup.co
Section 1 Section 2
y1 = 3.0 m y2 = 2.90 m
A1 = 30 m 2
A2 = 29 m2
P1 = 16 m P2 = 15.8 m
R1 = 1.875 m R2 = 1.835 m
K1 = 1/0.025 × 30 × (1.875)2/3 K2 = 1/0.025 × 29 × (1.835)2/3
= 1824.7 = 1738.9
Q = K Sf − 1781.3 Sf
2
V12 Q
= / 2g
2g 30
2
V22 Q
= / 2g
2g 29
V12 V22
hl = (h1 – h2) + −
2g 2g
V12 V22
hf = fall + −
2g 2g
V12 V22
= 0.12 + −
2g 2g
****
12
www.gradeup.co
13
www.gradeup.co
1
www.gradeup.co
ENGINEERING HYDROLOGY
3 EVAPORATION
1. EVAPORATION
It is the process in which liquid changes to gaseous state at the free surface of liquid below the
boiling point through the transfer of heat energy. Evaporation is basically a cooling process in
which latent heat of vaporization must be provided by water body.
1.1. Factors affecting evaporation
(i) Vapour pressure:
Rate of evaporation is directly proportional to the difference between saturated vapour
pressure (ew) and vapour pressure of air (ea).
(ii) Temperature:
Rate of evaporation is directly proportional to prevailing temperature.
(iii) Wind speed:
Rate of evaporation increases with increase in speed of the wind But if the wind velocity
is very high and is sufficient to remove all vapour on the liquid surface, then any further
increase in wind velocity will not affect the rate of evaporation.
(iv) Atmospheric pressure:
Rate of evaporation is inversely proportional to the atmospheric pressure.
(v) Soluble salts:
If water is saline which means it has salt in it then its vapour pressure decreases which
in turn results in a reduction in the rate of evaporation. As such evaporation from sea
water is about 2.3% less as compared to fresh water.
(vi) Depth of water body:
In summers, shallow depth of water permits rapid increase in temperature during the day
time and hence evaporation is greater in a shallow water bodies.
In winters, temperature of surrounding air falls but the temperature of deep water bodies
doesn’t fall easily and thus water remains relatively warmer consequently in winter
evaporation is more from a deep-water body. Above phenomenon is primarily due to heat
storage of water.
(vii) Surface Area:
Rate of evaporation is directly proportional to surface area.
2
www.gradeup.co
In order to convert pan evaporation values in actual evaporation value a pan coefficient
is used mathematically.
𝐿𝑎𝑘𝑒 𝐸𝑣𝑎𝑝𝑜𝑟𝑎𝑡𝑖𝑜𝑛 = 𝐶𝑃 × 𝑃𝑎𝑛 𝑒𝑣𝑎𝑝𝑜𝑟𝑎𝑡𝑖𝑜𝑛
3
www.gradeup.co
B. Empirical Method:
A large number of empirical equations are available to estimate lake evaporation using
commonly available meteorological data. Most formulae are based on the Dalton’s type
equation and can be expressed in the general form. It is called Meyer’s equation.
𝑣𝑎
𝐸 = 𝐾𝑚(𝑒𝑤 − 𝑒𝑎 ) (1 + )
16
Here, E = Rate of evaporation per day.
Km = constant which depends on size of water body.
ew = Saturated vapour pressure in mm of mercury.
ea = vapour pressure of air (ea) in mm of mercury.
Va = Mean monthly wind velocity in km/hr at a height of about 9m from the ground
surface.
C. Analytical methods:
The analytical methods for the determination of lake evaporation can be broadly
classified into three categories as
(i) Water budget method.
(ii) Energy balance method
(iii) Mass transfer method.
Example: A reservoir with surface area of 250 hectares has saturation vapor pressure at
water surface = 17.54 mm of Hg and actual vapor pressure of air = 6 mm. wind velocity
at 1 m above ground surface = 16 km/hr. Estimate the average daily evaporation from
the lake using Meyers formula.
Solution:
ew = 17.59 mm of Hg
ea = 6 mm
Wind speed of 9 m above ground = Ug
Ul C(1)1/7
=
Ug C(9)1/7
Ug = 21.9 km/hr
21.9
L = 0.36 (17.54 − 6 ) 1 +
16
= 9.84 mm/day.
4
www.gradeup.co
5
www.gradeup.co
Example: Estimate the PET of an area for the season January to April. The area in North
India at a latitude of 30°N with mean monthly temperatures and monthly percentage of
annual day time (hr) as follows.
Take,
K = 0.65
Use the Blaney riddle formula.
Solution:
The temperature converted to Fahrenheit and the calculations are performed in the table
Month Tf Ph Ph T f / 100
Ph T f
Now, = 16.07
100
By Blaney criddle formula.
εr = 2.54 × 16.07 × 0.65
= 26.53 cm
****
6
www.gradeup.co
7
www.gradeup.co
1
www.gradeup.co
HYDROLOGY
4 INFILTRATION
1. INTRODUCTION
Infiltration is the flow of water into the ground through the soil surface. Before infiltration of
water starts, certain initial losses are to be satisfied which are as follows.
1.1. Interception loss
This is the loss which is captured by vegetation and other such obstructions and which
subsequently evaporate some of the water which is trapped as interception loss runs
down the branches of the tree goes down the stem and is called stem fall.
1.2. Depression storage
When the precipitation of a storm reaches the ground, it must fill all the depressions
available on the ground before it can flow on the surface. The volume of water trapped
in these depressions is called Depression Storage.
2. INFILTRATION CAPACITY
The maximum rate at which a given soil at a given time can absorb water is defined as the
infiltration capacity. The infiltration capacity of a soil is high at the beginning of a storm and
how an exponential decay as the time elapses. Graphically it can be expressed as follows.
Equation of this curve has been developed in exponential form by Horton and is known as
Horton’s equation.
𝐹𝑡 = 𝐹𝑓 + (𝐹𝑖 − 𝐹𝑓 )𝑒 −𝑘𝑡
Here,
Ft = Infiltration rate or capacity at time ‘t’.
Ff = Final infiltration rate or capacity.
Fi = Initial infiltration rate or capacity.
K = Decay constant (T –1
or /s or / hr)
2
www.gradeup.co
Example: During a 3hr storm event it was observed that all abstractions other than infiltration
were negligible. The rainfall was idealized as three 1hr storms of intensities 10mm/hr,
20mm/hr and 10mm/hr respectively. The infiltration was idealized as Horton’s curve.
f = 6.8 + 8.7e–t where f is in mm/hr and t is in hr. calculate the depth of effective rainfall
Solution.
(6.8 + 8.7e
−t
= 20 × 1 + 10 × 1 – )dt
1
= 30 – [6.8 × 2 + 8.7 (e –1
– e–3)]
= 13.63 mm Ans.
3. INFILTRATION INDICES
3
www.gradeup.co
Example: A 6 hr rainstorm with hourly intensities of 7, 18, 15, 11 and 3mm/hr produced
a runoff of 39 mm. Calculate ϕ index in mm/hr.
Solution
Assume ϕ ≤3 mm/h
(7 – ) + (18 – ) + (25 – ) +(17 – ) + (11 – )+ (3 – ) = 39 mm
81 – 6 = 39
4 = 32
= 8 mm/hr
Hence index = 8mm/hr
3.2. W-index
This is the average infiltration rate during the entire duration of rainfall. In the calculation
of w-index, the initial losses are separated from total abstractions to refine the 𝜙 𝐼𝑛𝑑𝑒𝑥.
Mathematically it is defined as
𝑝 − 𝑟 − 𝑡𝑎
𝑤=
𝑡𝑐
Here, P = total storm precipitation (cm)
R = total run off (cm)
ta = minor loss in total duration
Tc = duration of the rainfall excess i.e. total.
W = defined average rate of infiltration (cm)
Example: For a storm of rainfall 80 cm, total storm runoff was recorded as 25 cm, taking
initial loss as 10 cm and duration of the rainfall excess as 2 hours, the value of w-index
will be
Solution.
P − R − Il
W=
tr
80 − 25 − 10
= cm/hr
2
45
= = 22.5 cm/hr
2
****
4
www.gradeup.co
5
www.gradeup.co
1
www.gradeup.co
HYDROLOGY
5 RUNOFF
1. INTRODUCTION
It is the graphical variation of discharge with time at any given location on the stream.
2
www.gradeup.co
The factors that affect the shape of the hydrograph can be broadly grouped into climatic factors
and physiographic factors.
3
www.gradeup.co
3.2. Size
Small basins be different from the large ones in terms of the relative importance of
various phases of the runoff phenomenon. In small catchment the overland flow phase is
predominant over the channel flow. Hence the land use and intensity of rainfall have
important role on the peak flood. On large basins these effects are suppressed as the
channel flow phase is more predominant.
3.3. Stream Density
It is defined as the total number of streams per unit area of the catchment.
3.4. Drainage density
This is defined as total length of streams per unit area of the catchment. Some examples
of drainage density are as follows.
4
www.gradeup.co
n=
5
www.gradeup.co
or
Here, n is in cm
A is in km2
t is in hr
O is in m3/sec
Example:
A 4-hour storm occurs over an 80 km2 watershed. The details of the catchment are as follows.
15 10 16 48 22 10
25 15 16 42 20 8
35 21 12 40 18 6
5 16 15 42 18 8
6
www.gradeup.co
2470
Total direct runoff = 115 + 2040 + 315 = 2470 km2 mm = 106 = 2.47 Mm3
103
5. UNIT HYDROGRAPH
The problem of predicting the flood hydrograph resulting from a known storm in a catchment
has received considerable attention. A large number of methods are proposed to solve this
problem and most popular and widely used method is the unit hydrograph method. This method
was first suggested by Sharman in 1932.
The unit hydrograph of a drainage basin is defined as a hydrograph of direct runoff resulting
from 1 cm of effective rainfall applied uniformly over the basin area at a uniform rate during a
specified period of time (D-hr).
7
www.gradeup.co
• Thus one can have a 6-hr unit hydrograph, 12-hr unit hydrograph etc.
1
• A 6-hr unit hydrograph will have an effective rainfall intensity of cm/hr.
6
• The effective rainfall intensity means the rainfall which will produce once runoff.
• In the D-hr unit hydrograph, D should not be more than any of the following (a) time of
concentration, (b) lag time or (c) period of rise
• Volume of water contained inside the unit hydrograph (i.e. area of unit hydrograph) is equal
to (1 cm ×Catchment area)
5.1. Assumption made in the theory of unit hydrograph (As proposed by Sherman)
is as follows:
The following assumptions are made while using the unit hydrograph principle :
1. Effective rainfall should be uniformly distributed over the basin, that is, if there are “N’
rain gauges spread uniformly over the basin, then all the gauges should record almost
same amount of rainfall during the specified time.
2. Effective rainfall is constant over the catchment during the unit time, i.e. intensity is
constant.
3. The direct runoff hydrograph for a given effective rainfall for a catchment is always the
same irrespective of when it occurs. (Time invariance) Hence, any previous rainfall event
is not considered to effect the new rainfall.
8
www.gradeup.co
Note: This antecedent precipitation is otherwise important because of its effect on soil-
infiltration rate, depression and detention storage, and hence, on the resultant
hydrograph.
4. The ordinates of the unit hydrograph are directly proportional to the effective rainfall
hyetograph ordinate. Hence, if a 6-h unit hydrograph due to 1 cm rainfall is given, then
a 6-h hydrograph due to 2 cm rainfall would just mean doubling the unit hydrograph
ordinates. Hence, the base of the resulting hydrograph (from the start or rise up to the
time when discharge becomes zero) also remains the same.
Example: A flood hydrograph of river draining a catchment of 189 km2 due to a 6 hour
isolated storm is in the form of a triangle with a base of 66 hour and a peak ordinate of
30 m3/s occurring at 10 hours from the start. Assuming zero base flow, develop the 6-
hour unit hydrograph for this catchment.
Sol. Catchment area = 189 km2
Flood hydrograph due to a 6h storm
∴ Volume of runoff
1
= 66 3600 30
2
= 3.564 × 106 m3
∴ Effective rainfall
3.564 106
= = 0.018857 m
189 106
= 1.8857 cm
For getting ordinates of 6 h unit hydrograph for the catchment, it is required to divide the
ordinates of the flood hydrograph of 1.8857 and the peak will be 15.91 m 3/s at 10 hr
from start. Its base will also be 66 hr.
9
www.gradeup.co
Step 3: Divide this volume by known area of catchment to get rainfall or rainfall excess
in (cm).
Step 4: Divide the ordinates of DRH by the depth of rainfall excess to obtain ordinates of
UH.
Step 5: Plot the ordinates of UH against time to get the UH of the catchment.
Example: The following are the ordinates of the flood hydrograph from a catchment area
of 780 km2 due to 6 hr storm. Derive the 6 hr unit hydrograph of the basin.
Time (hrs) 6 12 18 24 30 36 42 48 54 60 66 72 78
Discharge (m3/sec) 40 64 245 360 405 350 270 205 145 100 70 50 40
1794 6 60 60
6
100 cm = 4.968 cm
780 10
Therefore, the ordinates of UH are obtained by dividing the ordinates of DRH hydrograph
by rain excess 4.968 cm to get ordinates of UH.
Time (hrs) 6 12 18 24 30 36 42 48 54 60 66 72 78
Ordinates of
0 4.83 35.22 64.42 74.47 62.4 46.29 33.21 21.13 12.077 6.04 2.01 0
UH
10
www.gradeup.co
Since the rainfall rate is equal to the runoff rate at the equilibrium state, it follows that
1 A
Qe = A . km2 − cm / h = 2.778 m3 / s
D D
11
www.gradeup.co
Where A is the area of the basin in km2 and D is the duration of unit hydrograph in hours
which is used in the construction of the S-curve.
Consider two D-h S-curves A and B displaced by Th. If the ordinates of B are subtracted
from that of A, the resulting curve is a DRH produced by a rainfall excess of duration T h
1
and magnitude T cm. Hence if the ordinate difference of A and B, i.e. (SA – SB)
D
are divided by T/D, the resulting ordinates denote a hydrograph due to an ER of 1 cm
and of duration T i.e. a T-h unit hydrograph.
Example :
The 6-hour unit hydrograph of a catchment is triangular in shape with a base width of 64
hours and a peak ordinate of 30 m3/s. Calculate the equilibrium discharge of the S6-
curve of the basin.
Sol. Direct runoff
12
www.gradeup.co
1
64 3600 30 = 3.456 106 m3
2
ER = 1 cm
3.456 106
Area = 102
1
= 345.6 × 106 = 345.6 km2
Equilibrium discharge Qs = 2.778 A/D m 3/s
345.6
Qs = 2.778 = 160.0128 m3 / s
6
5.4. Limitations of Unit Hydrographs
(i) Precipitation must be rainfall only. Snow melt runoff can not be satisfactory
represented by unit hydrograph.
(ii) The catchment should not have unusually large storages in terms of tanks, ponds,
large flood bank storages, etc. which affect the linear relationship between storage and
discharge.
(iii) If the precipitation is decidedly non uniform, unit hydrographs can not be expected
to give good results.
Unit hydrographs can be derived if rainfall and runoff records are available for the basin under
consideration. But there are many basins, which are not gauged and for which unit-graphs may
be required. Hence some method of deriving unit hydrographs for ungauged basins is
necessary.
This is usually done by relating the selected basin characteristics to the unit hydrograph shape.
Once such relations are established between the basin parameters and unit hydrograph
parameters for the basins having sufficient data, the same relations are applied to get the unit
hydrograph of ungauged basins in the same hydro meteorologically homogeneous area from
the known basin parameters. The unit hydrograph thus obtained is known as Synthetic unit
hydrograph.
Note : Unit hydrograph developed from rainfall and stream flow data on a catchment applies
only for that catchment and point on stream where the stream flow data was measured.
Synthetic unit hydrograph is used to develop unit hydrograph for the locations on the stream
in the same catchment or for a nearby catchment of a similar character.
Snyder selected three parameters for development of SUH. Theses parameters relate the
catchment characteristics
(i) Basin time width T
(ii) Peak discharge QP
13
www.gradeup.co
(iii) Lag time i.e. basin lag time tp. (Snyder defined lag time as time interval from mid pt. of
rainfall to peak of UH (instead of centroid).
Where tp is in hr
Ct is a coefficient reflecting slope, land use, and associated storage characteristics of basin. Its
value varies between 1.35 to 1.65, average being 1.5
L = basin length measured along the water course from the basin divide to the gauging station
in km.
Lca = Distance of centroid of catchment from the gauging point (in km)
T is in hr
QP is in m3/s
A = Catchment area in km2
14
www.gradeup.co
Synder used the standard duration tr (or D-hr) in hr for unit hydrograph
tp
tr = Dhr =
5.5
If a synthetic unit hydrograph of other duration then D’hr is required, then lag time , tpr’ is given
by
D '− tr
tpr = tp =
4
To plot the smooth synthetic unit hydrograph, US army crops of engineering gave the width of
SUH as
5.87
W50 = 1.08
Qp
A
3.35 W50
W75 = =
Qp
1.08
1.75
A
Where W50 and W75 are the width of synthetic unit hydrograph in hr at 50% and 75% of QP
15
www.gradeup.co
Example: Derive 3-hr synthetic unit hydrograph of basin with following data with a catchment
area of 2500 km2
Length of main stream = 120 km
Distance of centroid from outlet = 80 km
Ct and Cp of the catchment are assumed to be 1.5 and 0.6 respectively. Use Synder’s method.
Sol. Tp = Ct(LLca)0.3
Ct = 1.5, L = 120 km, Lca = 80 km
tp = 23.484 hrs
Standard duration tr = D
tp 23.484
= = = 4.27 hrs
5.5 5.5
Required duration D’ = 3 hrs
Lag time
D '+ tr
tpr = tp + = tp + 0.25 (D '− tr )
4
Or tpr = 23.484 + 0.25 (3 – 4.27) = 23.16 hrs
0.6 2500 3
Qp = 2.78 = 180 m / sec
23.16
T = 72 + 3 × 23.16 = 141.48 hours
With these values Qp, T and tpr a smooth curve may be drawn.
16
www.gradeup.co
To plot the smooth synthetic unit hydrograph (SUH), US army corps Engineers gave the width
of SUH as :
5.87
W50 = 1.08
Qp
A
3.35
W75 = 1.08
Qp
A
Where W50 and W75 are widths in hrs of SUH at 50 % and 75% of QP
5.87
W50 = 1.08
= 100.62 hrs
180
2500
And
3.35
W75 = 1.08
= 57.42 hrs
180
2500
Hence the synthetic unit hydrograph (SUH) may be drawn as shown in Figure below :
17
www.gradeup.co
● With the help of IUH runoff is calculated using the following expression
tT
Run off = Q ( t ) = u ( t − t ') (I ( t ') dt )
0
****
18
www.gradeup.co
19
www.gradeup.co
1
www.gradeup.co
ENGINEERING HYDROLOGY
1. INTRODUCTION
2
www.gradeup.co
2. SATURATED ZONE
In the saturated zone, groundwater fills the pore spaces completely, and water is storedas in
a reservoir, having a hydrostatic pressure variation throughout itsdepth with atmospheric
pressure assumed to exist at the water table.
All earth materials, from soils to rocks have pore spaces. Although these pores are completely
saturated with water below the water table, from the groundwater utilization point of view,
only such material through which water moves easily and hence can be extracted with ease
are significant. On this basis, the saturated formation are classified into four categories:
1. Aquifer,
2. Aquitard
3. Aquiclude and
4. Aquifuge.
2.1. Aquifer
An aquifer is a saturated formation of earth material which not only stores water but
yields it in sufficient quantity relatively easily due to its high permeability. Deposits of
sand and gravel from good aquifers.
2.2. Aquitard
It is a formation through which only seepage is possible and thus the yield is insignificant
compared to an aquifer. A sandy clay unit is an example of aquitard.
2.3. Aquiclude
Formations like clay which is highly porous but not permeable due to very small size of
pores.
2.4. Aquifuge
It is a geological formation which is neither porous nor permeable. Massive compact rock
without any fractures is an aquifuge.
3
www.gradeup.co
3. TYPE OF AQUIFERS
• In unconfined aquifer the water table goes down if water is withdrawn from the aquifer
storage and the water table moves up if water is added into the aquifer storage.
• The water level in a large diameter dug wells tapping unconfined aquifer represents
water table. This aquifer is also known as water table aquifer or phreatic aquifer.
• A special case of unconfined aquifer is known as perched aquifer.
A perched aquifer is formed when the infiltrated rain water is intercepted within the
zone of Aeration by an impermeable layer and a local zone of saturation is formed. The
upper surface of such local zone of saturation is known as perched water table. The
perched aquifer occurs at higher elevation than the regional water table.
4
www.gradeup.co
The ground water within a confined aquifer occurs under pressure (known as confined
pressure or artisan pressure) greater than atmospheric pressure. When such confined
aquifer is pierced by a well, the water rises in the well due to release of pressure within
the confined aquifer. The level up to which water will rise in the well is known as
potentiometric level. This potentiometric level indicates the magnitude of pressure within
the confined aquifer. If the potentiometric level is above the ground surface a flowing
well results.
The area from which the infiltrated water enters the confined aquifer is known as
Recharge area.
Full penetrating well means the well which penetrates upto the bottom of the aquifer so that
flow is more or less radial.
5
www.gradeup.co
At a radial distance r from the well, if h is the piezometric head, the velocity of flow by Darcy’s
law is
dh
Vr = K
dr
The cylindrical surface through which this velocity occurs is 2πrB. Hence
dh
Q = (2rB ) K
dr
Q dr
= dh
2KB r
Integrating between limits r1 and r2 with the corresponding piezometric heads being h 1 and h2
respectively.
Q r
ln 2 = (h2 − h1 )
2KB r1
2KB (h2 − h1 )
Q= r
ln r2
1
This is the equilibrium equation for the steady flow in a confined aquifer. This equation is
popularly known as Theme’s equation.
Further, at the edge of the zone of influence, s =0, r 2 = R and h2= H;at the well wall r1 = rw,
h1 = hw and s1 = sw. Hence
2KBSw
Q=
R
ln
rw
6
www.gradeup.co
dh
Vr = K
dr
dh
Q = (2rh) Vr = 2r Kh
dr
Q dr
= h dh
2K r
Integrating between limits r1 and r2 where the water table depths are h 1 and h2 respectively
and on rearranging
Q=
(
k h22 − h12 )
r2
This is the equilibrium equation for a well in an unconfined aquifer (Thiems formula). As at the
edge of the zone of influence of radius R, H = saturated thickness of the aquifer. Hence
Q=
(
k H2 − h2w )
R
ln
rw
R is normally between 300-500 m.
where hw = depth of water in the pumping well of radius rw.
7
www.gradeup.co
Example: A 30-cm well completely penetrates an unconfined aquifer of saturated depth 40m.
After a long period of pumping at a steady rate of 1500 1pm, the drawdown in two observation
wells 25 and 75 m from the pumping well were found to be 3.5 and 2.0 m respectively.
Determine the transmissivity of the aquifer. What is the drawdown at the pumping well?
Solution:
1500 10−3
(a) Q = = 0.025 m3/s
60
h2 = 40.0 2.0 = 38.0 r2 = 75 m
h1 = 40.0 3.5 = 36.5 r1 = 25 m
From Eq. (9.49),
Q=
(
K h22 − h12 )
r2
ln
r1
K (38 ) − (36.5)
2 2
0.025 =
75
ln
25
K = 7.823 × 10–5 m/s
T = KH = 7.823 × 103 × 40 = 3.13 ×103 m2/s
(b) At the pumping well, rw = 0.15 m
Q=
(
K H12 − h2w )
r1
ln
rw
6. RECUPERATING TEST
Although the pumping test gives accurate value of safe yield, it sometimes becomes very
difficult to adjust the rate of pumping, so as to keep the well water level constant. In such
circumstances, recuperation test is adopted.
In this method, the water is first of all drained from the well at a fast rate so as to cause
sufficient drawdown. The pumping is then stopped. The water level in the well will start rising.
The time taken by the water to comeback to its normal level or some other measured level is
then noted. The discharge can then be worked out as below:
8
www.gradeup.co
AB = Static water level in the well before the pumping was started.
CD = Water level in the well when the pumping was stopped.
s1 = Depression head in the well at the time the pumping was stopped.
EF = Water level in the well at the noted time (say after a time T form when the pumping is
stopped).
s2 = Depression head in the well at time T after the pumping is stopped.
Let X.X be the position of water level at any time ti after the pumping was stopped and let the
corresponding depression head be s. Let ds be the decrease in depression head in a time dt
after the time t. Hence, in a time
t after the pumping is stopped, the water level recuperates by (s1 – s). It again recuperates by
ds in a time dt after this.
∴Volume of water entering the well in the small interval of time (dt)
=dV=A.dt (1)
where A is the cross-sectional area of the well at the bottom.
Also, if Q is the rate of recharge into the well at the time t under a depression heads, then the
volume of water entering the well in this small time interval is
= dV=Q.dt
But Q s
∴ Q=C’.s (2)
where C’ is a constant depending on the soil through which the water enters the well.
∴dV = C’.s.dt (3)
Equating (1) and (3), we get
– A.ds = C’.s dt
(The –v sign indicates that s decreases as t increases)
C.dt ds
or =−
A s
9
www.gradeup.co
s2
= −2.3log10
s1
C 2.3 s
∴ = log10 1
A T s2
C
Knowing the values of s1,s2 and T from the above test, the value of can be calculated. C’ is
A
called the specific capacity of the open well. Knowing the value of the discharge Q for a well
under a constant depression head H can be calculated as follows:
Q = C’.s
C
or Q = A.s
A
2.3 s
or Q = log 10 1 A.s
R s2
A and s ate known, the discharge for any amount of drawdown (s) can be easily worked out.
Example: During a recuperation test, the water level in an open well was depressed by
pumping by 2.5 meters and is recuperated by an amount of 1.6 meters in 70 minutes.
(a) Determine the yield from a well of 3 m diameter under a depression head of 3.5 meters.
(b) Also determine the diameter of the well to yield 10 liters/second under a depression head
of 2.5 meters.
Sol.
C 2.3 s
= log10 1
A T s2
where s1= Initial drawdown = 2.5 m
s2 = Final drawdown
= 2.5 – 1.6= 0.9m
T = Time = 70 minutes
70
= hr = 1.167 hr
60
C 2.3 2.5
∴ = log10
A 1.167 0.9
= 0.875 m3/hr/m2/m of depression head
10
www.gradeup.co
(a) Yield from a well of 3 m diameter, under a depression head of 3.5 m is given by
C
Q = A.s
A
= 0.875 32 3.5
4
= 21.65 m3/hr
= 6.02 liters/sec Ans.
Q = 10 liters/sec,
10 60 60 3
= m hr
1000
= 36 m3/hr
s = 2.5 m
C
Q = A.s
A
∴ 36 = (0.875) × A × 2.5
36
or A = = 16.46 m2
0.875 2.5
2
or d = 16.46
4
or d = 4.58 m; Say 4.6 m
Hence, the diameter of the required well = 4.6 m. Ans.
Example: During the recuperation test of a 4.0 m open well a recuperation of the depression
head from 2.5 m to 125 m was found to take place in 90 minutes. Determine the (i) specific
capacity per unit well area and (ii) yield of the well for a safe drawdown of 2.5 m (iii) what
would be the yield from a well of 5.0 m diameter for a drawdown of 2.25 m?
2
Sol. A= ( 4.0 ) = 12.566m2
4
K0 1 H
= ln 1
A Tr H2
****
11
www.gradeup.co
12