Introduction To Hydrology
Introduction To Hydrology
Hydrology is an earth science. It is the science that deals with the occurrence,
circulation and distribution of water of the earth and the earth’s atmosphere. As a
branch of earth science, it is concerned with the water in streams and lakes, rainfall
and snowfall, snow and ice on the land and water occurring below the earth’s surface
in the pores of the soil and rocks. In a general sense, hydrology is a very broad subject
of an inter-disciplinary nature drawing support from allied sciences, such as
meteorology, geology, statistics, chemistry, physics and fluid mechanics.
1
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
1) Determination of Maximum Probable Flood
2) Determination of Water yield from a basin
3) Study the groundwater development
4) Determination of maximum intensity of the storm
1.2.3 The main applications of engineering hydrology are:
1. Hydrology provides guidance for undergoing proper planning and
management of water resources.
2. Calculates rainfall, surface runoff, and precipitation.
3. It determines the water balance for a particular region.
4. It mitigates and predicts flood, landslide and drought risk in the region.
5. It estimates the water resource potential of the river basins
6. Enables real-time flood forecasting and flood warning..
7. Hydrology analyses the variations observed in the catchments by bringing a
relationship between the surface water and groundwater resources of the
catchment.
8. Hydrology studies the required reservoir capacity that is necessary for
irrigation and municipal water supply purpose during drought conditions.
9. It is used in the design and operation of hydraulic structures
10. It is used for hydropower generation.
11. Brings measures to control erosion and sediment
1.3 Hydrologic Cycle
Water occurs on the earth in all its three states, which are: liquid, solid and gaseous,
and in various degrees of motion. Evaporation of water from water bodies such as
oceans and lakes, formation and movement of clouds, rain and snowfall, streamflow
and groundwater movement are some examples of the dynamic aspects of water. The
various aspects of water related to the earth can be explained in terms of a cycle
known as the hydrologic cycle.
A convenient starting point to describe the cycle is in the oceans. Water in the oceans
evaporates due to the heat energy provided by solar radiation. The water vapor moves
upwards and forms clouds. While much of the clouds condense and fall back to the
oceans as rain, a part of the clouds is driven to the land areas by winds. There they
condense and precipitate onto the landmass as rain, snow, hail, sleet, etc. A part of the
precipitation may evaporate back to the atmosphere even while falling. Another part
2
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
may be intercepted by vegetation, structures and other such surface modifications
from which it may be either evaporated back to the atmosphere or move down to the
ground surface. A portion of the water that reaches the ground enters the earth’s
surface through
infiltration, enhances the moisture content of the soil and reaches the groundwater
body. Vegetation sends a portion of the water from under the ground surface back to
the atmosphere through the process of transpiration. The precipitation reaching the
ground surface after meeting the needs of infiltration and evaporation moves down the
natural slope over the surface and through a network of gullies, streams and rivers to
reach the ocean. The groundwater may come to the surface through springs and other
outlets after spending a considerably longer time than the surface flow. The portion of
the precipitation which by a variety of paths above and below the surface of the earth
reaches the stream channel is called runoff. Once it enters a stream channel, runoff
becomes stream flow.
The sequence of events as above is a simplistic picture of a very complex cycle that
has been taking place since the formation of the earth. It is seen that the hydrologic
cycle is a very vast and complicated cycle in which there are a large number of paths
of varying time scales. Further, it is a continuous recirculating cycle in the sense that
there is neither a beginning nor an end or a pause. Each path of the hydrologic cycle
involves one or more of the following aspects: (i) transportation of water, (ii)
temporary storage, and (iii) change of state. For example, (a) the process of rainfall
has the change of state and transportation, and (b) the groundwater path has storage
and transportation aspects.
Of the many processes involved in the hydrologic cycle, the most important are:
1. Precipitation
2. Runoff
3. Evaporation
4. Transpiration
3
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
5. Infiltration
6. Storage on the land surface (depression storage, ponds, Lakes, Ponds,
reservoir storage, etc.)
7. Soil moisture storage
8. Ground water storage
4
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
If the density of the inflow, outflow and storage volumes are the same,
ΔI- ΔO = ΔS
1.1
ΔI = inflow volume of water into the problem area during the time period, ΔO =
outflow volume of water from the problem area during the time period, and ΔS =
change in the storage of the water volume over and under the given area during the
given period. In applying this continuity equation [Eq. (1.)] to the paths of the
hydrologic cycle involving change of state, the volumes considered are the equivalent
volumes of water at a reference temperature. In hydrologic calculations, the volumes
are often expressed as average depths over the catchment area.
While realizing that all the terms in a hydrological water budget may not be known to
the same degree of accuracy, an expression for the water budget of a catchment for a
time interval Δt is written as:
P – R – G – E – T = ΔS
1.2
In this, P = precipitation, R = surface runoff, G = net groundwater flow out of the
catchment, E = evaporation, T = transpiration and ΔS = change in storage.
The storage S consists of three components as
S = Ss + Ssm + Sg 1.3
where Ss = surface water storage,
Ssm = water in storage as soil moisture, and
Sg= water in storage as groundwater.
1.5 World Water Balance
The total quantity of water in the world is estimated to be about 1386 million cubic
kilometers’ (M Km3). About 96.5% of this water is contained in the oceans as saline
water. Some of the water on the land amounting to about 1% of the total water is also
saline. Thus, only about 35 M Km 3 of fresh of fresh water is available. Out of this
about 10.6M km3 is both liquid and fresh and the remaining 24.4 M km3 is contained
in frozen state as ice in the polar regions and on mountain tops and glaciers. An
estimated distribution of water on the earth is given in Table below.
6
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Total:
(a) All kinds of
510 1386 100
water (b) Fresh
2.5 100
water 148.8 35
Table from WORLD WATER BALANCE AND WATER RESOURCES OF THE
EARTH, © UNESCO, 1975. Reproduced by the permission of UNESCO.
7
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
2 PRECIPITATION
2.1 Introduction
The term precipitation denotes all forms of water that reach the earth from the
atmosphere. The usual forms are rainfall, snowfall, hail, frost and dew. Of all these,
only the first two contribute significant amounts of water. Rainfall being the
predominant form of precipitation causing stream flow, especially the flood flow in a
majority of rivers. The magnitude of precipitation varies with time and space.
Differences in the magnitude of rainfall in various parts of a country at a given time
and variations of rainfall at a place in various seasons of the year are obvious and
need no elaboration. It is this variation that is responsible for many hydrological
problems, such as floods and droughts.
Condensation of water vapors in the air mass causes precipitation. Due to adiabatic
cooling, the rising air mass containing a significant number of water vapors becomes
saturated. Two processes are considered to be capable of supporting the growth of
droplets of sufficient mass(droplets from about500 to 4000 µm in diameter) to
overcome air resistance and consequently fall to the earth as precipitation. These are
known as the ice crystal process (Bergeron Process) and the coalescence process
(Collision – Coalescence Process)
Condensation of water vapors causes clouds to develop. In every cloud, there are
updrafts and downdrafts. The updraft determines cloud growth and height. The
stronger the updraft, the taller the cloud. As the liquid water accumulates, the updraft
weakens and the downdraft strengthens. Thereafter, Precipitation results. While all
clouds contain water, only some generate precipitation. Precipitated moisture from
clouds sometimes evaporates before reaching the earth's surface. Precipitation
happens when cloud droplets or ice crystals grow large enough to overcome
atmospheric updrafts. It signifies that specific processes are at work in a cloud that
produces rain.
8
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
These mechanisms can be explained by two processes:
1. Bergeron Process
The clouds contain ice crystals and super-cooled water droplets. When an ice crystal
collides with super-cooled water, it freezes the water. This method uses two properties
of water.
First property- A cloud's water droplets do not freeze at 0°C but remain liquid until -
40°C termed as Super-cooled water if disturbed, super-cooled water will freeze. It
needs nuclei to freeze on. But freezing nuclei are rare in the atmosphere. Hence, with
an increase in temperature, some water droplets will turn into ice. Encouraging the
formation of all-ice clouds, a single ice crystal is added to a supercooled water droplet
cloud.
Second property- Over ice crystals, the saturation vapor pressure (es) is lower than
over water. Between water and ice crystals, a vapor pressure gradient is created. Ice
crystals are formed at the expense of supercooled water. When these ice crystals get
large enough, they begin to tumble out of the cloud. These ice crystals melt and fall as
rain before reaching the ground.
This method is suited to clouds whose base does not exceed the freezing point. Warm
clouds are these. These clouds have several various-sized cloud droplets. Large drops
expand at the expense of smaller ones. So, they crash with the smaller droplets and
form part of it. Updrafts and downdrafts repeatedly lift and lower cloud droplets in a
huge cloud. So, these drops grow swiftly. Raindrops must have a diameter of at least
100µm.Cloud droplets collide to generate the size of a drizzle droplet. More collisions
produce larger drops and rain. Small droplets of uniform size may be present in
clouds that do not produce precipitation. A condition like this could lead to colloidal
stability in the clouds. Because of the small size of the droplets, a collision may not
occur and these cloud droplets can fall slowly and uniformly without colliding. As a
result, any clouds that lack the required size of cloud droplets may fail to produce
precipitation.
For precipitation to form: (i) the atmosphere must have moisture, (ii) there must be
sufficient nuclei present to aid condensation, (iii) weather conditions must be good for
9
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
condensation of water vapor to take place, and (iv) the products of condensation must
reach the earth.
Some of the common forms of precipitation are rain, snow, drizzle, glaze, sleet and
hail.
1. Rain:
10
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Rain is a form of precipitation that is in the form of water drops of a size that is larger
than 0.5mm. The maximum raindrop size is about 6mm. Drops of larger size break up
into smaller drops as it falls down on the Earth’s surface. Rainfall is the predominant
form of precipitation and therefore, the term precipitation is used synonymously with
rainfall.
2. Snow
Snow is another important form of precipitation. Snow consists of ice crystals which
usually combine to form flakes. When fresh, snow has an initial density varying from
0.06 to 0.15 g/cm3 and it is usual to assume an average density of 0.1 g/cm3.
3. Drizzle
A fine sprinkle of numerous water droplets of size less than 0.5 mm and intensity less
than 1 mm/h is known as drizzle. In this, the drops are so small that they appear to
float in the air.
4. Glaze
When rain or drizzle comes in contact with cold ground at around 0ºC, the water
drops freeze to form an ice coating called glaze or freezing rain.
5. Sleet
It is frozen raindrops of transparent grains which form when rain falls through air at
subfreezing temperature. In Britain, sleet denotes precipitation of snow and rain
simultaneously.
6. Hail
11
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
It is a showery precipitation in the form of irregular pellets or lumps of ice of size
more than 8 mm. Hails occur in violent thunderstorms in which vertical currents are
very strong.
1. Convective Precipitation:
2. Orographic Precipitation
Orographic precipitation results from the mechanical lifting of moist horizontal air
currents over natural barriers such as mountain ranges. Factors that are important in
this process include land elevation, local slope, orientation of land slope, and distance
from the moisture source.
12
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Cyclonic precipitation is associated with the movement of air masses from high
pressure regions to low-pressure regions. These pressure differences are created by
the unequal heating of the earth's surface. Cyclonic precipitation maybe classified as
frontal or non-frontal. Any barometric low can produce non frontal precipitation as air
is lifted through horizontal convergence of the inflow into a low-pressure area.
Frontal precipitation results from the lifting of warm air over cold air at the contact
zone between air masses having different characteristics. If the air masses are moving
so that warm air replaces colder air, the front is known as a warm front; if , on the
other hand, cold air displaces warm air, the front is said to be cold. If the front is not
in motion, it is said to be a stationary front.
13
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
An air mass moving from ocean to land results in increase fraction over land. While
moving towards land, these air masses rise gradually resulting in condensation and
rainfall over a limited area. Sometime these storm result in highly intense rainfall.
The characteristics of rainfall are the amount, the intensity, the duration, the frequency
or return period, and the seasonal distribution.
The amount is of course important to the overall hydrologic cycle and replenishment
of the soil water, and the amount is an accumulation or product of the intensity times
the duration. For example, the amount may be the same for a high intensity short
duration rainfall as it is for a low intensity long duration rainfall.
However, the intensity and duration can have a large influence on whether the rainfall
infiltrates or becomes surface runoff. Higher rainfall intensity produces larger size
raindrops which have more impact energy, and thus higher intensity storms can
damage delicate vegetation and bare soil. High intensity storms can literally displace
soil particles, causing soil crusting or starting the soil erosion process. High intensity
14
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
storms may also overwhelm the soils’ ability to infiltrate the rainfall at the same rate,
causing infiltration-excess runoff.
The duration refers to the length of time rainfall occurs. A high intensity rainfall for a
short duration may affect tender seedlings, but it will not likely have much effect on
soil erosion and runoff. Rainfall of longer duration can significantly affect infiltration,
runoff, and soil erosion processes.
The frequency, or more specifically, the return period refers to how often rainfall
occurs at a particular amount or intensity and duration. For example, rainfall return
periods are referred to as 100 year-1 hour rainfall or 100 year-24-hour rainfall to
define the probability that a given amount will fall within a given time period.
The seasonal distribution of rainfall refers to the time of year when various rainfall
amounts occur. The seasonal distribution determines when surface runoff or deep
percolation are most likely to occur or if irrigation is needed. Since the seasonal
distribution of rainfall varies in different parts of the country, practices used or
recommended in one part of the country may not necessarily be appropriate in
another.
The rainfall data is a prerequisite for the study of water resources and design of water
harvesting structures. Therefore, its correct measurement is essential. Since most
estimates are based on rainfall data, information regarding the depth and intensity of
rainfall is of great importance.
Precipitation is expressed in terms of the depth to which rainfall water would stand on
an area if all the rain were collected on it. Thus, 1 cm of rainfall over a catchment area
of l km2 represents a volume of water equal to 104 m3. In the case of snowfall, an
equivalent depth of water is used as the depth of precipitation. The precipitation is
collected and measured in a rain gauge. Terms such as pluviometer, ombrometer and
hyetometer are also sometimes used to designate a rain gauge.
15
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
For setting up a rain gauge, the following considerations are important:
The ground must be level and in the open and the instrument must present a
horizontal catch surface.
The gauge must be set as near the ground as possible to reduce wind effects
but it must be sufficiently high to prevent splashing, flooding, etc.
The instrument must be surrounded by an open fenced area of at least 5.5 m *
5.5 m. No object should be nearer to the instrument than 30 m or twice the
height of the obstruction.
Rain gauges can be broadly classified into two categories as (i) nonrecording rain
gauges, and (ii) recording gauges.
1. Nonrecording Gauges
The nonrecording gauge extensively used is the Symons’ gauge. It essentially consists
of a circular collecting area of 12.7 cm (5.0 inch) diameter connected to a funnel. The
rim of the collector is set in a horizontal plane at a height of 30.5 cm above the ground
level. The funnel discharges the rainfall catch into a receiving vessel. The funnel and
receiving vessel are housed in a metallic container. Figure below shows the details of
the installation. Water contained in the receiving vessel is measured by a suitably
graduated measuring glass, with an accuracy up to 0.1 mm.
16
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Figure 7: Non recording Raingauge
2. Recording Gauges
Recording gauges produce a continuous plot of rainfall against time and provide
valuable data of intensity and duration of rainfall for hydrological analysis of storms.
The following are some of the commonly used recording rain gauges.
This is a 30.5 cm size rain gauge adopted for use by the US Weather Bureau. The
catch from the funnel falls onto one of a pair of small buckets. These buckets are so
balanced that when 0.25 mm of rainfall collects in one bucket, it tips and brings the
other one in position. The water from the tipped bucket is collected in a storage can.
The tipping actuates an electrically driven pen to trace a record on the clockwork-
driven chart. The water collected in the storage can is measured at regular intervals to
provide the total rainfall and also serve as a check. It may be noted that the record
from the tipping bucket gives data on the intensity of rainfall. Further, the instrument
is ideally suited for digitalizing of the output signal.
In this rain gauge, the catch from the funnel empties into a bucket mounted on a
weighing scale. The weight of the bucket and its contents are recorded on a
clockwork-driven chart. The clockwork mechanism has the capacity to run for as long
as one week. This instrument gives a plot of the accumulated rainfall against the
elapsed time, i.e., the mass curve of rainfall. In some instruments of this type, the
17
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
recording unit is so constructed that the pen reverses its direction at every preset
value, say 7.5 cm (3 in.) so that a continuous plot of storm is obtained.
This type of recording rain gauge is also known as float-type gauge. Here, the rainfall
collected by a funnel-shaped collector is led into a float chamber causing a float to
rise. As the float rises, a pen attached to the float through a lever system records the
elevation of the float on a rotating drum driven by a clockwork mechanism. A syphon
arrangement empties the float chamber when the float has reached a preset maximum
level.
A typical chart from this type of rain gauge is shown in Fig. below. This chart shows a
rainfall of 53.8 mm in 30 h. The vertical lines in the pen-trace correspond to the
sudden emptying of the float chamber by syphon action which resets the pen to zero
level. It is obvious that the natural syphon-type recording rain gauge gives a plot of
the mass curve of rainfall.
18
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
These rain gauges are of the recording type and contain electronic units to transmit the
data on rainfall to a base station both at regular intervals and on interrogation. The
tipping-bucket type rain gauge, being ideally suited, is usually adopted for this
purpose. Any of the other types of recording rain gauges can also be used equally
effectively. Telemetering gauges are of utmost use in gathering rainfall data from
mountainous and generally inaccessible places.
The meteorological radar is a powerful instrument for measuring the areal extent,
location and movement of rain storms. Further, the amounts of rainfall over large
areas can be determined through the radar with a good degree of accuracy. The radar
emits a regular succession of pulses of electromagnetic radiation in a narrow beam.
Since the catching area of a rain gauge is very small compared to the areal extent of a
storm, it is obvious that to get a representative picture of a storm over a catchment, the
number of rain gauges should be as large as possible, i.e., the catchment area per
gauge should be small. On the other hand, economic considerations to a large extent
19
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
and other considerations, such as topography, accessibility, etc. to some extent, restrict
the number of gauges to be maintained. Hence, one aims at an optimum density of
gauges from which reasonably accurate information about the storms can be obtained.
Towards this, the World Meteorological Organization (WMO) recommends the
following densities.
In arid and polar zones: 1 station for 1500–10,000 km2 depending on the
feasibility.
Ten per cent of rain gauge stations should be equipped with self-recording gauges to
know the intensities of rainfall.
If there are already some rain gauge stations in a catchment, the optimal number of
stations that should exist to have an assigned percentage of error in the estimation of
mean rainfall is obtained by statistical analysis as:
2
Cv
N=( ) 2.1
έ
Cv =( 100*σm -1)/ P
2.2
σ m−1=√ ¿ ¿
20
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Pi = precipitation magnitude in the ith station
Consider the existing m raingauge. They have a mean rainfall of P and a coefficient
of variation of Cv. To know the percentage of error (έx) of estimation of mean in the
existing system of m raingauge, Eq. (2.3) can be rewritten (by changing N by m) as
2
Cv
m=( )
εex
2.4
εex =Cv/√ m
2.5
In the above Eq. (2.5), the term εex represents the expected error (in percentage) in
the estimation of the mean, P. It is a measure of the accuracy of estimation of mean
precipitation in the existing system and is called standard error in the estimation of the
mean. In calculating the number of rain gauges N for a given level of error by using
Eq. (2.1) it is usual to take εex as 10%. It is seen that if the value of εex is small, the
number of rain gauges required will be more. According to WMO recommendations,
at least 10% of total rain gauges should be of self-recording type.
Before using the rainfall data in application, it is necessary to first check the data for
continuity and consistency. The continuity of record may be broken with missing data
due to many reasons such as damage or fault in the rain gauge during a period. If the
conditions of installation of a rain gauge undergo a change, the records of the period
before and after the change would be inconsistent. Methods of pre-cleaning the
records for such problems of continuity and consistency are indicated in the following
sections.
The point observation from a precipitation gage may have a short break in the record
because of instrument failure or absence of the observer. Thus, it is often necessary to
estimate the missing record using data from the neighboring station. The following
methods are most commonly used for estimating the missing records.
21
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
1. Simple Arithmetic Method
2. Normal Ratio Method
3. Modified normal ratio method
4. Inverse distance method
In the estimation of missing data from a rain gauge station, performance of a group of
neighbouring stations including the one with missing data are considered. A
comparison of the recordings of these stations are made by using their normal rainfall
as standard of comparison. The normal rainfall is the average value of rainfall at a
particular date, month or year over a specified 30-year period. The 30-year normal is
recomputed every decade. Thus, the term normal annual precipitation at a station A
means the average annual precipitation at station A based on a specified 30-year of
record. Insertion of missing data to a station record must be done sparingly. If too
many data are estimated, the quality of the total data set may be diluted due to
interpolation. sometimes, if too many gaps exist in a record, it may be worthwhile to
neglect that station than to have a station record with too much of interpolated data. A
WMO guideline states that not more than 5 to 10% of a record should have
interpolated data. A commonly used procedure for estimating missing data of a station
is given below.
1
Px= [P 1+ P 2+…+ Pm ]
M
2.6
(b) Normal Precipitation: Given the annual precipitation values P1, P2, P3,
…..Pm at neighbouring M stations 1, 2, 3, …, M respectively, it is required to
find the missing annual precipitation Px at a station X not included in the
above M stations. Further, the normal annual precipitations N1, N2, …Ni,…at
each of the above (M+1) stations, including station X, are known.
22
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Nx P1 P 2 Pm
Px= [ + +…+ ]
m N1 N 2 Nm
2.7
∑ D1i /b
i=1
2.8
Where r is normal rainfall, D is the distance between the index station i and
the gauge station with missing data or ungagged station, n is the number of
index stations and b is the constant by which the distance is weighted
(normally 1.5-2.0) commonly used D0.5.
(d) Inverse Distance Method: The inverse distance method has been expected
to be the most accurate method as compare to other two methods discussed
above.
Amount of rainfall to be estimated at a location is a function of;
I. rainfall measured at the surrounding index stations
II. distance to each index station from the ungauged location
Rainfall rx at station x is given by;
n
r
∑ ( Dib )
i=1 i
r x= n
∑ (D1b )
i=1 i
2.9
b = 2 is commonly used.
23
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
change took place. Some of the common causes for inconsistency of record are
the station X (i.e., ∑Px) and the accumulated values of the average of the group of
bases
stations (i.e., ∑Pav) are calculated starting from the latest record. Values of ∑Px are
plotted against ∑Pav for various consecutive time periods. A decided break in the
slope of the resulting plot indicates a change in the precipitation regime of station X.
The precipitation values at station X beyond the period of c6
Mc
Pcx =P x
Ma
2.10
In this way, the older records are brought to the new regime of the station. It is
apparent that the more homogeneous the base station records are, the more accurate
will be the corrected values at station X. A change in the slope is normally taken as
24
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
significant only where it persists for more than five years. The double-mass curve is
also helpful in checking systematic arithmetical errors in transferring rainfall data
from one record to another.
There are different homogeneity tests used to test the homogeneity of the rainfall data.
Some of them are Standard normal homogeneity test (SNHT), Buishand range (BR)
test, Pettitt test, and von Neumann ratio (VNR) test are selected.
Over a large area like watershed (or catchment) of a stream, there will be several such
stations, and the average depth of rainfall over the entire area can be estimated by one
of the following methods:
The most direct approach is to use the arithmetic average of gauged quantities. This
procedure is satisfactory if gauges are uniformly distributed and the topography is
flat.
Thus,
25
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
P 1+ P2+ P 3+…+ Pn 1
P=
n n
∑ Pi
2.11
In this method, the rainfall recorded at each station is given a weightage on the basis
of an area closest to the station. The entire area is divided into number of triangular
areas by joining adjacent rain gauge stations with straight lines
n
The ratio Ai/A is called the weightage factor for each station.
3. Isohyetal Method:
26
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
An isohyet is a line joining points of equal rainfall magnitude. In the Isohyetal
method, the catchment area is drawn to scale and the raingauge stations are marked.
The recorded values for which areal average P is to be determined are then marked
on the plot at appropriate stations. Neighbouring stations outside the catchment are
also considered. The isohyets of various values are then drawn by considering point
rain-falls as guides and interpolating between them by the eye (Fig. 11). The
procedure is similar to the drawing of elevation contours based on spot levels.
27
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
The isohyet method is superior to the other two methods especially when the stations
are large in number.
The rainfall data measured at a place using a measuring device is known as point (or
station) rainfall data. For small areas of less than 50 km 2 , point rainfall may be taken
as the average depth over the area. In large areas, there will be a network of rain-
gauge stations. Depending upon the need, data can be listed as daily, weekly, monthly,
seasonal or annual values for various periods. The variation of rainfall with respect to
time may be shown graphically by:
a) Hyetographs
A mass curve of rainfall is a plot of cumulative depth of rainfall against time. From
the mass curve, the total depth of rainfall and intensity of rainfall at any instant of
time can be found.
The amount of rainfall for any increment of time is the difference between the
ordinates at the beginning and end of the time increments,
The intensity of rainfall at any time is the slope of the mass curve (i.e., i=
ΔP/Δt) at that time.
28
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
A mass curve of rainfall is always a rising curve and may have some horizontal
sections which indicates periods of no rainfall. The mass curve for the design storm is
generally obtained by maximizing the mass curves of the severe storms in the basin
can determine magnitude and duration of storm. Slope of curve give intensity at
various time.
A flow-duration curve (FDC) represents the relationship between the magnitude and
frequency of daily, weekly, monthly (or some other time interval of) streamflow for a
particular river basin, providing an estimate of the percentage of time a given
streamflow was equaled or exceeded over a historical period. An FDC provides a
simple, yet comprehensive, graphical view of the overall historical variability
associated with streamflow in a river basin.
29
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Simply, it is a graphical representation of runoff. It is denoted by FDC. It relates flow
rate with duration.
Discharge and time are calculated at first. FDC is then plotted which is used to learn
or study floods, evaluate many dependable stream data, design drainage works, etc.
In the FDC:
30
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
7. It used in the calculation of the sediment load and dissolve solids load of
stream.
8. It used for comparing stemflow data.
9. It is used for evaluation many dependable flows in water engineering projects.
Evaporation is the process by which water is transferred from the land and water
masses of the earth to the atmosphere. Transpiration is the counterpart for plants. It is
the processes by which soil moisture taken up by vegetation is eventually evaporated
as it exists at plant pores. Evaporation and transpiration combined
(evapotranspiration) generally constitute the largest components of losses in rainfall
runoff sequences. Accordingly, good estimates of evapotranspiration are request for
hydrologic modelling.
3.2 Evaporation
Evaporation is the process during which a liquid change into a gas. The process of
evaporation of water in nature is one of the fundamental components of the
hydrological cycle by which is one of the vapors through absorption of heat energy.
This is the only form of moisture transfer from land and oceans into the atmosphere.
31
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
3.2.1 Factors Affecting the Rate of Evaporation
There are four primary factors that affect the rate of evaporation; temperature, surface
area, wind, and humidity.
The wet clothes dry faster on a hot sunny day because of the high temperature. When
the temperature is high, the kinetic energy of the water particles increases, leading to
the quicker conversion of liquid water into its gaseous state. This explains why hot
water evaporates more quickly than cold water.
When the temperature of the liquid (water, in this case) is increased, it results in an
increase in the kinetic energy of the individual molecules that constitute the liquid.
This increase in energy makes it easier for the liquid molecules to overcome the
intermolecular forces of attraction (that holds the liquid together) and escape into the
atmosphere as a gas. It can be noted that the temperature of the surrounding can also
contribute towards faster evaporation since hot environments can transfer their heat to
the liquid.
Suppose there is water in a wide shallow bowl, and the same amount of water is also
kept in a narrow glass. In which one of the two situations will the water evaporate
quickly? The water in the bowl will evaporate much faster than in the glass because
the larger surface area gives more space for the water particles to evaporate. You
might have also noticed that the wet clothes dry quickly when they are spread out
properly on the clothesline compared to the clothes that are not spread out properly.
This is again due to the larger surface area of the spread-out wet clothes, leading to
faster water evaporation.
Evaporation is a surface phenomenon, which implies that the process only occurs at
the surface of the liquid. During evaporation, the molecules of the liquid that are
present at its surface overcome the intermolecular forces of attraction to break away
from the liquid and escape into the atmosphere as a gas (or vapor).
Liquids are known to not have any definite shape – they assume the shape of their
container. Therefore, the greater the surface area of the container (or the greater the
surface occupied by the liquid), the greater the number of liquid molecules present at
32
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
the surface. The larger the number of liquid molecules is at the surface, the larger the
number of molecules that will break away from the liquid and become a gas at a given
point of time, which contributes to an increase in the rate of evaporation.
When the clothes are spread out on a clothesline on a windy day, they will dry up
quicker than on a still day. This is because the water vapor in the surrounding air will
move away from the clothes by the wind, leading to faster water evaporation from wet
clothes.
Evaporation increases the humidity of the atmosphere that immediately surrounds the
liquid. This humid air takes some time to dissipate into the rest of the atmosphere. The
presence of a breeze, a powerful wind, or some other form of air circulation can speed
up this process and make the environment of the liquid less humid. Therefore, by
decreasing the humidity of the liquid’s surrounding, a powerful breeze or wind can
increase the rate at which the liquid evaporates.
Why is it that high humidity makes you feel even hotter on a hot day? Humidity is a
measure of the amount of water vapor in the air. We already sweat on a hot day;
however, if it is windy, we feel cooler because of the evaporation of the sweat.
However, sweat evaporates more slowly in high humidity because there is already a
lot of water vapor in the surrounding air. The slower evaporation rate reduces the
potential for evaporative cooling. This is why the wet clothes will not dry up quickly
when there is high humidity in the surrounding air.
There is a limit to how much water vapor the atmosphere can hold. This limit is
proportional to the temperature (which implies that hot air can hold more water vapor
than cold air). Evaporation can be viewed as an equilibrium process through which
the amount of water vapor in the atmosphere reaches an equilibrium with the amount
of water in the surface of the liquid. Therefore, the greater the amount of water vapor
in the atmosphere over the liquid, the slower the rate of evaporation. Also, it can be
noted that an increase in temperature with constant humidity will contribute to an
increase in the rate of evaporation since hot air can hold a greater amount of water
vapor.
33
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
3.3 Transpiration
Transpiration is the process by which water vapor leaver the living plant body and
enters the atmosphere. It involves continuous movement of water from the soil into
roots, through the stem and cut through the leaves to the atmosphere. The process
includes cuticular transpiration or direct evaporation in to the atmosphere from moist
membranes through the cuticle, and stomatal transpiration or outward diffusion into
the atmosphere through the stomata and lenticels vapor previously evaporated from
imbibed membranes, into intracellular space within the plant.
Plants put down roots into the soil to draw up water and nutrients into its stems and
leaves. Some of this water is returned to the air by transpiration. Transpiration rates
vary widely depending on weather and other conditions, such as
Type of plant: Plants transpire water at different rates. Some plants which grow in
arid regions, such as cacti and succulents, conserve precious water by transpiring less
water than other plants.
Soil type and saturation: Clay particles are small (smaller than 0.002 mm), holding
onto water whereas sand particles which are large (0.05-2 mm) release water readily
(think of how water disappears into the sand quickly at the beach). When moisture is
lacking, plants can begin to senesce (premature aging, which can result in leaf loss)
and transpire less water.
Humidity: As the relative humidity of the air surrounding the plant rises the
transpiration rate falls. It is easier for water to evaporate into dryer air than into more
saturated air.
34
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Wind & air movement: Increased movement of the air around a plant will result in a
higher transpiration rate. Wind will move the air around, with the result that the more
saturated air close to the leaf is replaced by drier air.
Surface Area of the Leaves: A leaf having more surface area will show more
transpiration rate than the leaf with a lesser surface area.
The evapotranspiration rate is normally expressed in millimeters (mm) per unit time.
The rate expresses the amount of water lost from a cropped surface in units of water
depth. The time unit can be an hour, day, decade, month or even an entire growing
period or year.
35
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Evapotranspiration is defined as the total losses of water from vegetation – both as
evaporation from the soil and transpiration from the plants. In the cropped area it is
difficult to determine evaporation separately from soil surface and transpiration from
the plants and usually both the processes are combined in a single term called
evapotranspiration.
Availability of water is also very important. When water is adequately available at the
site of evaporation i.e., soil or plant surface, the rate of evaporation is primarily
controlled by meteorological factors like solar radiation, wind, temperature, vapor
pressure deficit etc. constituting evaporative demand of the atmosphere.
Radiation
Radiation is the main source of energy for the ET process. It depends on the global
solar radiation flux density and vegetation albedo. A darker vegetation absorbs more
incident solar radiation and evapotranspiration more. A thorough understanding of the
factors controlling the energy balance of a cropped soil enables making accurate
estimates or predictions of ET and irrigation water requirements. It also facilitates
more effective irrigation water management.
Temperature
Over the course of a day, an increase of the air temperature causes an increase on the
saturation deficit triggering a higher evaporative demand in the air, and leading to
high ET rates
36
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Air relative humidity acts in conjunction with temperature. The higher relative
humidity, the lesser the evaporative demand and, therefore, the lower.
Wind
Advection represents the horizontal transport of energy from a drier area to another
more humid, and such additional energy is utilized in the ET process. Wind also helps
remove water vapor near the plants to other regions.
Crop factors such as crop species, radiation reflection coefficient, leaf area index
(LAI) in different growth stages of the plant, plant height and root depth (depth of the
radicular system) are influenced on the crop ET, which are described details as
follows:
Crop species
This factor is related to the foliar architecture (spatial distribution of the leaves),
internal resistance of the plant to water transport, and other morphological aspects
(number, size, and distribution of stomata, etc.), which exert a direct influence on ET.
Radiation reflection influences directly net radiation availability for the ET process.
The darker the vegetation, the lower the reflection coefficient and the higher net
radiation.
LAI in different growth stages of a plant is directly related to the size of the
transpiring foliar surface, for the larger leaf area the larger the transpiring surface, and
the higher the potential for water use.
Plant height
Plant height also influences the ET. Taller and rougher plants interact more efficiently
with the atmosphere in motion, extracting more energy from the air and, therefore,
increasing ET.
37
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Rooting depth of plant is directly related to the volume of soil explored by the roots,
aiming at meeting the atmospheric hydric demand. A superficial radicular system, for
exploring a smaller soil volume, keeps the crop more susceptible to drying periods.
Crop orientations
Soil texture and soil structure are both unique properties of the soil that will have a
profound effect on the behavior of soils. Both the properties influence the crop ET
through influencing the water holding capacity of the soils. Clay soils have higher
water holding capacity than sandy soils and are proficient of preserving a more
persistent crop ET rate for longer.
Chemical/physical impediments
Inhibitions system limit the growth of the radicular system of a plant, affecting the
root system of plants to explore a wider volume of soil both in dry and rainy seasons.
In the rainy season, soil with any physical obstructions gets soaking wet suffocating
the roots. While in the dry season, the volume of available water is reduced from the
roots of a plant; as a results root system of a plant is deepening into the soil for
searching available water.
38
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
3.5 evaporation and transpiration process
39
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
40
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
The pan is placed on a wooden platform of height 15cm above ground
level to allow free air circulation below the pan
Evaporation is measured by measuring the depth of water in a stilling well
with a hook gauge
41
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
(3) Colorado Sunken Pan
920mm square pan made of unpainted GI sheet, 460mm deep, and buried into the
ground within 100mm of the top
Main advantage of this pan – its aerodynamic and radiation characteristics are
similar to that of a lake
Disadvantages – difficult to detect leaks, expensive to install, extra care is needed
to keep the surrounding area free from tall grass, dust etc
Evaporation pans are not exact models of large reservoirs. Their major drawbacks are
the following: – They differ from reservoirs in the heat storage capacity and heat
transfer characteristics from the sides and the bottom (sunken and floating pans aim to
minimize this problem). Hence evaporation from a pan depends to some extent on its
size (Evaporation from a pan of about 3m dia is almost the same as that from a large
lake whereas that from a pan of about 1m dia is about 20% in excess of this). – The
height of the rim in an evaporation pan affects wind action over the water surface in
the pan. Also, it casts a shadow of varying size on the water surface. – The heat
transfer characteristics of the pan material is different from that of a reservoir. Hence
evaporation measured from a pan has to be corrected to get the evaporation from a
42
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
large lake under identical climatic and exposure conditions. Lake Evaporation = Pan
Coefficient (Cp) *Pan Evaporation.
Accounts for the effect of pressure in addition to the wind speed effect
Wind Velocity
In the lower part of the atmosphere, up to a height of about 500m above the
ground level, wind velocity follows the one-seventh power law as
44
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
45
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
3.6.2 Measurements of Transpiration
46
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Owing to the difficulty in obtaining accurate direct measurement of pan evaporation
under field conditions, evaporation is often predicted on the basis of climatological
data. The approaches followed are to relate the magnitude and variation of
evapotranspiration to one or more climatic factors (temperature, day length, humidity,
wind, sunshine etc.). The more commonly used empirical formulae in estimating
evapotranspiration are:
a) Blaney-Criddle Method
b) Thornthwaite Method
c) Hargreaves’ Method
The Blaney-Criddle method is recommended for periods of one month or longer.
a. Blaney-Criddle Method
This method requires the use of only two factors, temperature and information of day
light hours which is a factor based purely on the latitude of the place. Using Blaney-
Criddle approach, potential evapotranspiration can be expressed as follows, in metric
unit:
(3.1)
Where, PET = potential evapotranspiration, mm of water per day (mean value over
the month)
P= monthly percent of total day time hours of the year
T= mean monthly temp. In °C (Average of daily max and minvalues)
b. Thornthwaite Method
Thornthwaite method is based on the assumption of an exponential relationship
between mean monthly temperature and mean monthly consumptive use.
3.2
Where, e = unadjusted PET (cm per month)
t = mean air temperature (°C)
I = annual or seasonal heat index, the summation of 12 values of monthly
heat indices (i) when, i = ( t / 5 )1.514
3.3
c. Hargreaves’ Method
Hargreaves based on his work on data from grass lysimeter, proposed the following
relationship to estimate ET,
3.4
Where, PET = reference crop potential consumptive use
47
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
t = mean daily temperature (°C)
Rs= incident solar radiation in lang lay/day, it can be calculated using the
following relationship,
(13.7)
3.5
Where, S is the percent possible sunshine hour and R so is the clear daysolar radiation
in langlay/day.
d. FAO Penman-Monteith Method
The FAO Penman-Monteith method is used to estimate reference evapotranspiration.
The equation is:
3.6
Where, ET0 = reference evapotranspiration [mm day-1]
Rn = net radiation at the crop surface [MJ m-2 day-1]
T = mean daily air temperature at 2 m height [°C]
u2 = wind speed at 2 m height [m/s]
es = saturation vapour pressure [kPa]
ea = actual vapour pressure [kPa]
es - ea = saturation vapour pressure deficit [kPa]
Δ = slope of vapour pressure curve [kPa°C-1]
Υ = psychometric constant [kPa°C-1]
The reference evapotranspiration, ET0, provides a standard to which:
a) Evapotranspiration at different periods of the year or in another region can be
compared.
b) Evapotranspiration of other crops can be related.
4 INFILTRATION
4.1 Introduction
Infiltration is that process by which precipitation moves down ward through the
surface of the earth and replenishes soil moisture, recharges aquifers, and ultimately
supports streamflow during dry periods. Along with interception, depression storage,
and storm period evaporation, it determines the availability, if any, of the precipitation
input for generating over land flows. Furthermore, infiltration rates influence the
48
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
timing of over land flow inputs to channelized systems. Accordingly, infiltration is an
important component of any hydrologic model.
During the rain infiltration loss occurs quickly almost exclusively from the water that
has reached the ground surface. The water infiltrating into the soil moves downward
through larger soil pores under the force of gravity. The smaller surface pores take in
water by capillarity. The downward moving water is also sucked in by capillary pores.
The gravitational water moves towards the ground water following the path of least
resistance. When the capillary pores at the surface are filled and intake capacity
reduced infiltration rate decreases. As a trend the rate of infiltration is high in the
beginning. It decreases rapidly in the initial stages and then slowly till it approaches a
nearly constant rate in about 30 to 90 minutes depending upon the type of soil.
The typical infiltration curves for soils with different land uses viz. bare soil, soil with
grass cover and soil covered with grain crop are shown in fig below.
This trend of infiltration process is observed not only because of filling up surface
capillary pores with water but also due to changes in the soil such as dispersion of
aggregates, puddling of the surface layer, impact of the rain drops, swelling of
colloids, closing of soil cracks, type of vegetal cover etc.
49
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
The process of infiltration is affected by many different factors. Important factors are
the following:
It is already made clear that the water cannot continue to enter soil more rapidly than
it is transmitted downward. The conditions at the surface, therefore, cannot increase
infiltration unless the transmission capacity of the soil profile is adequate.
The continuity of non-capillary or large pores provides easy paths for percolating
water. If the subsoil formation has coarse texture the water may infiltrate into the soil
so quickly that no water will be left for runoff even if rainfall is quite heavy. On the
contrary clayey soils after soaking some water in the initial stages of the rainfall may
swell considerably. It makes the soil almost watertight and infiltration may get
reduced to practically negligible extent.
Even if the subsoil has excellent under drainage but at the surface soil pores are sealed
due to turbid water or by in wash of fine soil particles it may prevent entry of water
into the soil and infiltration rate will be low.
When the soil is fairly dry the rate of infiltration into the soil is quite high. The
infiltration rate diminishes as the soil-moisture storage capacity is exhausted. After
this infiltration rate equals transmission rate. The rate of infiltration in early phases of
a rainfall will be less if the soil pores are still filled from previous rain storm.
If saturated soil mass gets frozen due to severe low temperature it becomes nearly
impermeable. It affects the infiltration.
50
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
(vi) Human Activities on Soil Surface:
51
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
52
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Also, a dry soil can absorb more water than one whose pores are already full (Fig.
above). The land use has a significant influence on fp. For example, a forest soil rich
in organic matter will have a much higher value of f p under identical conditions than
the same soil in an urban area where it is subjected to compaction.
Surface of Entry
At the soil surface, the impact of raindrops causes the fines in the soil to be displaced
and these in turn can clog the pore spaces in the upper layers of the soil. This is an
important factor affecting the infiltration capacity. Thus, a surface covered with grass
54
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
and other vegetation which can reduce this process has a pronounced influence on the
value of fp.
Fluid Characteristics
Water infiltrating into the soil will have many impurities, both in solution and in
suspension. The turbidity of the water, especially the clay and colloid content, is an
important factor and such suspended particles block the fine pores in the soil and
reduce its infiltration capacity. The temperature of the water is a factor in the sense
that it affects the viscosity of the water which in turn affects the infiltration rate.
Contamination of the water by dissolved salts can affect the soil structure and in turn
affect the infiltration rate.
4.3.2.1 Modelling Infiltration capacity
4.1
Thus, the curve Fp(t) vs time in Fig. above is the mass curve of infiltration. It may be
df p (t)
f p (t)=
dt
4.2
55
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Many equations have been proposed to express the curves f p(t) or Fp(t) for use in
hydrological analysis. In this section four such equations will be described.
1. Horton’s Equation
Horton defined infiltration capacity as “the maximum rate at which a given soil can
absorb rainfall when the soil is in a specified condition. Based on observations of
infiltration in the field, he deduced that the decrease in infiltration capacity over time
during an infiltration event was due to flow-restricting changes occurring in a thin
layer at the soil surface. He believed the key changes included: 1) packing of the soil
surface by raindrop impact, 2) swelling of the soil, and 3) plugging of surface pores
with fine materials.
Horton expressed the decay of infiltration capacity with time as an exponential decay
given by
−k h t
f p=f c +(f 0−f c )e for 0 ≥ t ≤ tc 4.3
where fp = infiltration capacity at any time t from the start of the rainfall
Kh= Horton’s decay coefficient which depends upon soil characteristics and
vegetation cover
The difficulty of determining the variation of the three parameters f 0, fc and kh with
soil characteristics and antecedent moisture conditions preclude the general use of Eq
4.3.
56
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
infiltration capacity could be expressed as
−1
1
f p= st 2 + K 4.5
2
where a and b are local parameters with a > 0 and 0 < b < 1. The infiltration capacity
would now be expressed by Eq. (3.6) as
( b −1 )
f p=( ab ) t
4.7
Green and Ampt proposed a model for infiltration capacity based on Darcy’s law as
η SC
f p=K (1+ )
Fp
4.8
n
f p=m+
Fp
4.9
Infiltration indices are the parameters for measuring the various infiltration losses that
occur through various means. The process of water entering the soil is called
infiltration. The infiltration occurs on the ground surface. Below the surface, the
penetration further is called percolation. Whatever rainfall occurs on the surface of the
earth, some quantities infiltrate.
57
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Types of Infiltration Indices
The average infiltration rate is called infiltration index and two types of indices are in
commonly used:
1. ф Index:
The ф index is that portion of average rate of rainfall during any storm which gets
lost by the processes of interception, depression storage and infiltration taken
together. It can, therefore, be defined as that rate of average rainfall during any
storm beyond which the volume of remaining rainfall equals the volume of direct
surface runoff. The index can be calculated from a hyetograph (time versus
intensity of rainfall graph) of the storm in such a way that the rainfall volume in
excess of this rate will equal the volume of the storm runoff Fig. 24.
If the rainfall intensity throughout the storm remains equal to or more than ф
index then the ф index represents basin recharge because ф index represents sum
total of infiltration, interception and depression storage.
Ф =1/t (P-R)
4.10
Where
P = Precipitation (rainfall).
R = Surface runoff
2. W Index:
58
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
This index gives the average rate of infiltration for that time period of the storm
rainfall during which rainfall intensity is greater than W. Thus, it can be said to be
refinement over ф index which apart from infiltration also includes interception
and depression storage.
W = (P-Q-S)/t
59
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
5 RUNOFF
5.1 Introduction
Runoff means the draining or flowing off of precipitation from a catchment area
through a surface channel. It thus represents the output from the catchment in a given
unit of time.
Consider a catchment area receiving precipitation. For a given precipitation, the
evapotranspiration, initial loss, infiltration and detention storage requirements will
have to be first satisfied before the commencement of runoff. When these are
satisfied, the excess precipitation moves over the land surfaces to reach smaller
channels. This portion of the runoff is called overland flow and involves building up
of a storage over the surface and draining off of the same. Usually, the lengths and
depths of overland flow are small and the flow is in the laminar regime. Flows from
several small channels join bigger channels and flows from these in turn combine to
form a larger stream, and
so on, till the flow reaches the catchment outlet. The flow in this mode, where it
travels all the time over the surface as overland flow and through the channels as
open-channel flow and reaches the catchment outlet is called surface runoff. A part of
the precipitation that infilters moves laterally through upper crusts of the soil and
returns to the surface at some location away from the point of entry into the soil.
5.2 Components of runoff
60
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
The amount of interflow depends on the geological conditions of the catchment. A
fairly pervious soil overlying a hard impermeable surface is conducive to large
interflows. Depending upon the time delay between the infiltration and the outflow,
the interflow is sometimes classified into prompt interflow, i.e., the interflow with the
least time lag and delayed interflow.
Another route for the infiltered water is to undergo deep percolation and reach the
groundwater storage in the soil. The groundwater follows a complicated and long path
of travel and ultimately reaches the surface. The time lag, i.e., the difference in time
between the entry into the soil and outflows from it is very large, being of the order of
months and years. This part of runoff is called groundwater runoff or groundwater
flow. Groundwater flow provides the dry-weather flow in perennial streams.
Based on the time delay between the precipitation and the runoff, the runoff is
classified into two categories; as
61
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
is expressed by water balance equation as
RN = (Ro – Vr) + Vd + E + EX + ΔS
5.1
where RN = natural flow volume in time Δt
Ro = observed flow volume in time Δt at the terminal site
Vr= volume of return flow from irrigation, domestic water supply and
industrial use
Vd = volume diverted out of the stream for irrigation, domestic water supply
and industrial use
E = net evaporation losses from reservoirs on the stream
EX = net export of water from the basin
ΔS = change in the storage volumes of water storage bodies on the stream
Factors affecting runoff
The main factors affecting the runoff from a catchment area are:
i. Precipitation characteristics
ii. Shape and size of catchment
iii. Topography
iv. Geologic characteristics
v. Meteorological characteristics
vi. Storage characteristics of a catchment
Precipitation Characteristics
Precipitation is the most important factor, which affects runoff. The important
characteristics of precipitation are duration, intensity and areal distribution.
Duration: Total runoff depends on the duration of rainstorm. For a given rainfall
intensity and other conditions, a longer duration rainfall event will result in more
runoff.
Intensity: Rainfall intensity influences both rate and volume of runoff. The runoff
volume and also runoff rate will be greater for an intense rainfall event than for less
intense event.
Areal Distribution: It also influences both the rate and volume of runoff. Generally,
the maximum rate and volume of runoff occurs when the entire watershed contributes.
Shape and Size of Catchment
The runoff from a catchment depends upon the size, shape and location of the
catchment. The following are the general observations:
a) More intense rainfall events are generally distributed over a relatively smaller area,
i.e., larger the area lower will be the intensity of rainfall.
62
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
b) The peak normally decreases as the area of the basin increase. (Peak flow per unit
area)
c) larger basins give a more constant minimum flow than the smaller ones. (Effect of
local rains and greater capacity of the ground-water reservoir)
d) Fan shaped catchments give greater runoff because tributaries are nearly of same
size and hence time of concentration of runoff is nearly same. On the contrary,
discharges over fern leaf arrangement of tributaries are distributed over long period
because of the different lengths of tributaries.
63
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
A study of the annual hydrographs of streams enables one to classify the streams in
to three classes as (i) perennial, (ii) intermittent, and (iii) ephemeral.
A perennial stream is one that always carries some flow in it (Fig. 24). There is
considerable amount of groundwater flow throughout the year. Even during the dry
seasons, the water table will be above the bed of the stream
A catchment, also known as drainage basin, drainage area and watershed, is the basic
hydrologic unit in the analysis of runoff phenomenon. Physical characteristics of the
catchment such as area, shape, slope, and drainage channel pattern in the catchment
are some of the major static characteristics that affect the volume of the surface runoff
and shape of the runoff hydrograph from a catchment due to a storm. Study of the
catchment characteristics forms an important component of the subject of
geomorphology and is studied under the topic head of quantitative geomorphology.
Some of the important geometrical parameters of the catchment as related to runoff
process are briefly described below.
1. Area
2. Stream Order
It is a classification reflecting the pattern of branches that unite to form the trunk
stream leaving the catchment. The smallest stream at the start of the network is
designated as of ‘order 1’. Two channels of order 1 when joined produce a stream of
‘order 2’.Two channels of order 2 when joined produce a stream of order 3, and so on.
Note that when a lower order stream (say order 2) meets a higher order stream (say
order 3) the order of the resulting stream is still the higher order stream entering the
65
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
confluence with, order 3 in this case. The trunk stream discharging out of the drainage
basin has the highest order. Fig. 27 shows schematically a drainage basin with a
stream system of order 4.
The ratio of the number of streams (Ns) of all orders to the area of the basin (A) is
Ns
known as stream density (Sd) of the basin. Thus, Sd = .
A
5.2
Stream density represents the number of streams per unit area and is indicative of the
pattern of drainage channels
in the catchment.
4. Drainage Density, Dd
Drainage density is the ratio of the total length of streams of all orders within a basin
66
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
The drainage density is a measure of how well the drainage basin is drained by the
stream network.
The drainage density is called coarse, medium or fine depending on the value of Dd
Dd > 5 km/km2.
Drainage density depends upon the climate and physical characteristics and geology
of the catchment. Generally, impermeable soil material in the catchment, thin land
cover and steep slope cause more channels to move the runoff rapidly. Hence, a high
drainage density may indicate one or more of the following:
Maximum basin relief is the elevation difference (in meters) between the catchment
outlet and the highest point on the basin perimeter.
6. Slope
Usually, a slope profile prepared along the main stream is used to characterize the
slope of the catchment. For a given stretch of the stream, the ratio of horizontal
distance between the two end parts of the stretch to the difference in elevation
between these two points gives the slope of the stretch. As a rule, the catchment slope
is the highest at the beginning of the stream and gradually decreases as one moves
along the stream to the basin outlet. Slope is an important parameter in many
watershed simulation models.
7. Length, L
The length of the catchment is defined as the length of the main stream measured
from the basin outlet to the remotest point on the basin boundary. The main stream is
identified by starting from the basin outlet and moving up the catchment. At any
branching point the largest order branch is taken. If there is a branch of two streams of
the same order, the one with the largest catchment area is taken as the main steam.
67
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
8. Shape
There are a number of ways of designating and quantifying the shape of a catchment.
Some of the commonly used shape parameters are shown in Table 5 along with their
definitions. In this table A = catchment area, L = catchment length and P = perimeter
of the catchment.
9. Hypsometric Curve
68
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
The relationship between rainfall in a period and the corresponding runoff is quite
complex and is influenced by a host of factors relating to the catchment and climate.
Further, there is the problem of paucity of data which forces one to adopt simple
correlations for adequate estimation of runoff. One of the most common methods is to
correlate seasonal or annual measured runoff values (R) with corresponding rainfall
(P) values. A commonly adopted method is to fit a linear regression line between R
and P and to accept the result if the correlation coefficient is nearer unity. The
equation of the straight-line regression between runoff R and rainfall P is
R = Ap+b
5.6
a=N ¿ ¿
5.7
be calculated as
r =N ¿ ¿ 5.9
69
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
The value of r lies between 0 and 1 as R can have only positive correlation with P.
The value of 0.6 < r < 1.0 indicates good correlation. Further, it should be noted that R
≥ 0. The term r2 is known as the coefficient of determination.
70
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
AI' A2, A3 ••· An are areas of n sub-parts of the watershed having runoff coefficients
Cl' C2 ' C3 ..Cn,
II. Curve Number Method
The curve number is an empirical parameter used in hydrology for predicting direct
runoff or infiltration from rainfall excess. The curve number method was developed
by the U.S. Soil Conservation Service or SCS. The number is still popularly known
as a "SCS runoff curve number" .The SCS Runoff Curve Number (CN) method is
described as follows, the SCS runoff equation is:
2
(P−0.2 S )
Q=
P+0.8 S
5.12
2
(P−Ia)
Q=
( P−Ia ) +S
5.13
Where,
Q = Runoff, mm;
P = Rainfall, mm;
S = Potential maximum retention after runoff begins, mm; and
Ia = Initial abstraction, mm.
Ia = O.2S
S is related to the soil and cover conditions of the watershed through the CN by the
following relation:
25400
S= −254
CN
5.14
The value of CN ranges from 0 for most permeable surface and 100 for the
impervious (concrete) surface. For paved surface, S will be zero inducing runoff equal
to rainfall. The major factors that determine CN are the hydrologic soil group, cover
type, treatment, hydrologic condition, and antecedent runoff condition. The runoff
depth as a function of rainfall depth and CN is shown in Fig. 29.
71
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
72
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
73
INTRODUCTION TO HYDROLOGY HANDOUT
Department of Water Resources and Irrigation Engineering by Habtamu M.
Time-area method is based on the concept of time-area histogram, i.e., a
histogram of contributing catchment subareas (segmentations of the catchment).
To develop a time-area histogram, the catchment's time of concentration is divided
into a number of time intervals based on the hydrology time step. Cumulative time
at the end of each time interval is used to divide the catchment into segments
delimited by isochrones lines, i.e., the loci of points of equal travel time to the
catchment outlet. Based on the end times of each time interval and the area of the
corresponding segment the time-area histogram table can be built and used in the
hydrologic routing, as shown in the following figure.
As stated above, the method assumes that the hydrology computation has time
intervals (time segments) within the time of concentration of the catchment, and
based on these time segments the total area is also divided into corresponding sub-
areas (area segments). In general, a minimum of 2 time/area segments are needed
for good results, therefore the hydrology time step should be less than half of the
time of concentration. If user specifies a larger time step the model will give a
warning and results are less reliable.
74