Precipitation: Lecture Note Civil Engineering1

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 28

Chapter 2 Precipitation2021

2. PRECIPITATION
2.1 Introduction
Precipitation is that part of atmospheric moisture, which reaches the earth's surface in
different forms. Hydrologists start working when the precipitation reaches the ground.
This connects hydrology with meteorology. There is a great variation of precipitation in
space and time. In India, atmospheric moisture (or weather system) causes good
precipitation during June-October and nearly dry weather during the remaining periods.
Some of the precipitation that might get intercepted while reaching the ground by trees
and buildings and evaporates back is called the initial loss. The other part meets
requirements like depression storage and infiltrates into the ground. The excess rainfall
flows in streams to large water bodies. Factors like soil-type, vegetation, geology and
topography of the area largely determine the quantity of rainfall excess available as
stream flow from the precipitable water.
Nearly one-fourth of the total precipitation that falls on land reaches large water bodies as
direct runoff. The balance three-fourths of water returns back to the atmosphere at
different times as evaporation.
The essential requirements for precipitation to occur are: (i) moisture in the atmosphere,
(ii) presence of nuclei around which condensation of vapour takes place, (iii) dynamic
cooling responsible for condensation of water vapour and (iv) precipitation product must
reach the ground in some form. In the atmosphere, condensation nuclei are available in
plenty. Dust particles and smoke wastes form excellent nuclei. When such a nucleus
collects freezing moisture from the cloud and grows to a size greater than 1 mm, it forms
a rain drop and falls. Years back, cloud seeding was attempted for making artificial rain
in India. Cloud seeding should be attempted only when there is sufficient moisture in the
atmosphere and dynamic cooling exists. In the process, artificial nuclei like portland
cement, salt powder, carbon-dioxide powder, ice-powder, clay or silver iodide powder is
introduced into the cloud by aircraft, balloon, rocket, projectile and other devices. Due to
rapid industrialization, cloud seeding is no more attempted. Owing to the great variability
of meteorological processes, cloud seeding has so far provided inconsistent results.
Location of a region with respect to general circulation, latitude, topography, orographic
features and distance from moisture source are the responsible factors for rainfall to
occur. Orographic barriers force a precipitating system to rain on its windward side and
the leeward side receives less precipitation. Such precipitation distributions are noticed in
western Ghat mountain ranges and at the south of Himalayan ranges of India. West of
Western Ghats and the north-eastern states of India usually receive heavy precipitation
due to the orographic features. Precipitation decreases with altitude and, therefore, we
can say that higher is the altitude lesser is the precipitation. When the prevailing wind
moves over ocean containing warm water, it picks up high moisture. This gives very high
precipitation.

Lecture note Civil Engineering1


Chapter 2 Precipitation2021

During low pressure time, India and the adjoining areas experience such type of heavy
precipitations. While other conditions remain same, a place far away from any coast
receives less precipitation than one closer to sea. Location of a place with respect to the
prevailing wind is very effective in weather systems. Formation of a desert in western
part of India and the floods occurring frequently in north-eastern part are part of this
system. Presence of orographic features and oceans adds to the advantage.
Here precipitation measurement, network design, analysis and presentation of
precipitation data are discussed to the extent a partitioning civil engineer dealing with
water resource projects ought to know.
2.2 Forms of Precipitation
Any product of atmospheric water must reach the surface of earth after condensation.
However, fog and frost are not part of precipitation as they are not falling moisture. Some
common forms of precipitation explained below are: (i) rain, (ii) snow, (iii) drizzle, (iv)
glaze, (v) sleet, (vi) hail and (vii) dew.
(i) Rain: When precipitation reaches the surface of earth in the form of droplets of water,
we call it rain. The size of drops varies from 0.5 mm to 6 mm as drops larger than this
size are found to breakup during their fall in the air. Rain is considered as light if
intensity of rainfall is up to 2.5 mm/h, moderate from 2.5 to 7.5 mmlh and heavy over 7.5
mmlh.
(ii) Snow: It is precipitation in the form of ice-crystals, normally hexagonal in shape.
Snow reaches the earth's surface either separately or combines together to form flakes.
The density of snow is usually 0.10 glcm3, which means that 10 cm of snowfall is
equivalent to 1.0 cm of rainfall.
(iii) Drizzle: Drizzle is defined as water droplets of size less than 0.5 mm. It reaches the
ground with intensity less than 1.00 mm/h. These water droplets are so light that they
appear to be floating in air.
(iv) Glaze: It is the drizzle, which freezes immediately in contact with cold objects of the
earth's surface.
(v) Sleet: Where rain falls through air of subfreezing temperature, the drops freeze to
form grains of ice, called sleet. Sometimes snow and rain precipitates simultaneously.
The rain drops under this circumstance are half frozen.
(vi) Hail: It is the precipitating rain in the form of any irregular form of ice with size
ranging from 5.0 mm to 50 mm or above. Cumulonimbus convective. Clouds with strong
vertical currents are responsible for the formation of hail. The densities of hails are
normally 0.8 gmlcm3• While falling they combine together to form bigger sizes.
(vii) Dew: During nights when surface of the objects on earth cools by radiation, the
moisture present in atmosphere condenses on the surface of these objects forming water
droplets called dew.

Lecture note Civil Engineering2


Chapter 2 Precipitation2021

2.3 Types of Precipitation


Moisture is always present in the atmosphere and there is no shortage of condensation
nuclei in the present days due to rapid industrialization. Adiabatic cooling of the moist air
through lifting is the main cause of condensation. Precipitation, normally classified
according to the factors responsible for lifting the air mass, is of following types: (i)
convective, (ii) orographic, (iii) cyclonic and (iv) thunder storms.
i) Convective
Unequal heating at the surface of earth is the main reason for this type of precipitation. In
summer days air in contact with surface of the earth gets heated up, expands and rises due
to lesser density. Surrounding air rushes to replace it and inturn gets heated up and rises.
In the process, increasing quantities of water vapour are taken up by the air. When the air
mass reaches to a great height, temperature falls below the dew point. At this stage
condensation releases latent heat of 539 caVg of water. This heat is added to the air at the
height of condensation causing further heating to its upper air. This forces the air mass to
move up. As more moist air from the surrounding joins the system, more energy is added.
The vertical air currents develop tremendous velocities and dynamic cooling takes place.
Depending on the moisture content, cooling and other factors, the precipitation intensity
varies from light showers to cloud bursts amounting to 300 mm or more rainfall per hour.
Such type of precipitation covers limited aerial extent and is normally found near
equatorial zone. Sometimes strong upward wind currents exceeding 150 kmph freezes the
rain drops to form hail. Multiple rise and fall of hail due to very strong upward currents
may add to its size. Hail storms in India during summer are very common.

Fig 2.1 Convective Precipitation


ii) Orographic
Orographic or mountain-range barriers cause lifting of the air masses. Dynamic cooling
takes place causing precipitation on the side of the blowing wind as shown in Fig. 3.1.
Precipitation is normally heavier on the windward side and lighter on leeward side.
Lecture note Civil Engineering3
Chapter 2 Precipitation2021

In India, heavy precipitation in Himalayan region and at the western coast are mainly due
to orographic features associated with the south west wind carrying sufficient quantity of
moisture, while passing over Arabian sea. Orographic precipitation gives medium to high
intensity rainfall and continues for longer duration.

Fig. 2.2 Precipitation due to Orographic Features


iii) Cyclonic
A cyclone is a low pressure area surrounded by a larger high pressure area. When low
pressure occurs in an area, especially over large water bodies, air from the surrounding
rushes, causing the air at the low pressure zone to lift. Such a type of cyclone is called
Tropical cyclone or simply cyclone in India, Typhoon in south-east Asia and Hurricane in
America. The cyclone centre is called eye and is a calm area. This zone is surrounded by
a very strong wind zone with wind speed sometimes exceeding 300 km/h. Aerial extent
of the cyclone extends to a few hundred kilometers. The pressure distribution diagram for
this zones in the form of Isobars (lines joining equal pressure points with respect to mean
sea level) are closely spaced. There is steep decrease in pressure towards the eye.
Tropical cyclones originate near equator at 5 to 10 degree latitude and move towards
higher latitudes in a path guided by number of weather factors. Sea water temperature
between 25 and 27°C is the most favorable for formation of cyclones. The system derives
its energy from sea vapour and grows in size. Once the cyclone crosses over to land, the
energy source is cutoff; it becomes week and disappears gradually. The rainfall is
normally heavy in the entire zone travelled by a cyclone. In northern hemisphere,
cyclones move in anti-clockwise direction and in southern hemisphere they move
clockwise. Cyclonic storms move at the rate of about 30-50 km/h and give medium to
high intensity rainfall over a larger area.
A cyclone formed outside the tropical zone near the boundary between warm and cold air
is called extratropical cyclone. Aerial extents of such cyclones are large, while wind and
precipitation characteristics are lesser than the tropical cyclones.

Lecture note Civil Engineering4


Chapter 2 Precipitation2021

An anticyclone is an area of high pressure in which winds tend to blow spirally outward
in clockwise direction in the northern hemisphere and anticlockwise in southern
hemisphere. Weather is normally calm and such anticyclones are not associated with rain.
(iv) Thunder Storms
An air mass which moves from sea to land gets increased friction over land. These air
masses rise gradually as they move inland, giving rise to condensation and precipitation
over a limited area. Winter rainfall in southern part of India and Indonesia are mainly due
to this process. Sometimes thunder storms result in very intense rainfall.
2.4 Measurement of Precipitation
Precipitation is measured as depth of water equivalent from all forms that would
accumulate on a horizontal surface if there are no losses. The vertical depth of water is
expressed in millimeters and tenths in metric system and in inches and hundredths in FPS
system. Precipitation data is a basic input for the study of any water resources system and
should be measured extensively. Due to its great variability in space and time this natural
parameter is recorded continuously. Rainfall is collected and measured in instruments
called rain gauges. If snow is the form of precipitation, then it is collected in snow-
gauges, melted and its water equivalent is recorded. A gauge in its simplest form is a
horizontal circular opening aperture of known cross sectional area in the form of a
cylindrical vessel. The circular opening leads its catch to a collecting and measuring jar.
Three types of instruments generally used for measurement of rainfall are: (i) non-
recording gauge, (ii) recording gauge and (iii) weather radars.
2.3.1. Non-Recording Rain Gauge
Various types of non-recording rain gauges are available, for instance, A Symons type or
a Standard gauge. The two types of gauges are alike but differ in their size proportion and
material. A standard gauge is made up of reinforced fiberglass polyester material of
different combinations of collectors and bottles. Details of installation of a standard type
of rain gauge is shown in Fig.2.3

Fig.2.3 Standard Gauge


Lecture note Civil Engineering5
Chapter 2 Precipitation2021

The circular collector opening has an area of either 100 or 200 cm 2 from which rain
enters into the receiving vessel through a funnel. The collectors are interchangeable. Top
of the circular opening is placed at a standard height of 30 cm above ground level. The
metal container should be fixed to a concrete block of 60 cm x 60 cm x 60 cm as shown
in Fig. 2.3.
The rain catch collected in the bottle is taken out and poured into a graduated measuring
glass jar (chosen accordingly for 100 cm or 200 cm collector), which gives directly the
depth of rainfall for the day. Degree of accuracy of a graduated jar is 0.1mm. If
precipitation on a day is very heavy then more readings should be taken and summed up
to give the rainfall depth for the day along with its final observation taken. Such a system
can also be used for snow measurement but the snow should be melted to arrive at the
equivalent depth of water recorded. This system needs the service of an attending
observer who monitors the gauge at regular intervals, usually daily and hourly during
continuous and heavy precipitations.
Since the above system does not record rain but simply collects, it is called a non-
recording gauge. The ratio of volume of water collected in cm3 divided by the area of
opening of the gauge mouth in cm2 gives the depth of rainfall for the day. If any non-
standard measuring jar is used, then it should be calibrated to give directly the depth of
rainfall in centimeters and millimeters.
2.3.2. Recording Type Rain-Gauges
These rain gauges give a continuous record of rainfall at a place over time. Such gauges
give all the required information ofa storm like the onset and cessation of rain, i.e.,
duration of the storm, intensity and the cumulative rainfall. The recording gauges are
commonly installed along with a non-recording type gauge for the purpose of checking
and calibration. Many types of recording gauges are available in the market. One may use
a tipping bucket type, weighing bucket type or syphon type gauge on the consideration of
their merits and demerits to suit the conditions prevailing over the site. When such a
gauge is fitted with an electronic device to transmit rainfall data to a base station then it is
known as Telemetering gauge.
A. Tipping Bucket Type: A Stevens tipping bucket type of rain gauge consists of a
200 mm collector that directs the rain water through a funnel into a twocompartmental
bucket. The size of each bucket is 0.25 mm of rain. Once rain water fills up a bucket, it
over-balances and the water tips down to the casing

Lecture note Civil Engineering6


Chapter 2 Precipitation2021

Figure 2.4. Tipping bucket type rain gauge


of the container bringing thereby the second bucket to its measuring position beneath the
funnel. A tipping bucket rain gauge is shown in Fig. 2.4. Tipping of the buckets actuates
an electric circuit which records the number of tips during rain. This type of gauge can be
installed at an inaccessible area from where electric pulses generated due to the tipping of
the buckets is recorded at the control room far away from the gauge location.
Disadvantages of such a type of gauge are (i) when tipping of buckets takes place, rainfall
at that instant is not recorded, (ii) very high intensity rainfall gives close signals, which
can make it difficult to record the number of tips and (iii) calibration of tips may change
due to rusting and dirt accumulation.
B. Weighing Bucket Type: This type of gauge can be used for recording rainfall as well
as snow. Rain is collected in a receiver bucket supported on a spring balance. A
mechanical lever arm of the balance is connected with a pen which touches a clock
mounted drum with a graph paper. As it rains, the weight of the bucket gradually
increases. This changes the position of the pan of the balance. With time the pen marks a
line c" the continuously moving graph paper. The record shows the accumulation of
precipitation over time. The recording can be taken after 24 hours or 7 days depending on
the clock and drum size. Such gauges are normally used in USA and are becoming
increasingly popular. A sketch of the weighing bucket type rain gauge is shown in Fig.
2.5.

Lecture note Civil Engineering7


Chapter 2 Precipitation2021

Figure 2.5. Weighing bucket type rain gauge


Disadvantages of such a gauge are (i) when very heavy precipitation occurs, there is good
chance that the bucket will overflow and (ii) such instruments are costly. When the pen
reaches the upper end of the graph paper it reverses its love, thus continuously recording
the precipitation over time.
C. Syphon (Float) Type: Such a rain gauge is shown in Fig. 2.6. Rain entering the
gauge is led to a float chamber through a funnel. With increase in rainwater in the
chamber, the float rises. A pen mounted on the float through a lever system touches a
graph chart warped around the circumference of the drum. The drum is mounted on a
mechanical clock. Clock of the rotating drum is wound either for 24 hours or 7 days after
which the graph chart is to be replaced. By the time the pen reaches the top of the graph
the float also reached the top of the chamber. At this point syphonic action takes place in
the chamber and all the water in the chamber below the float empties. The pen comes
back to its original zero position. If there is no rainfall, the pen moves horizontally over
the graph paper at that level. One syphonic action means 10 mm of rainfall and the time
taken to collect the depth of rain can be noted from the horizontal axis of the graph
paper. A used graph paper for syphon type gauge is shown in Fig. 2.7. The beginning
and the end of the storm, its intensity, duration, distribution of rain and the depth of total
storm precipitation can easily be obtained from the plot of the graph.

Lecture note Civil Engineering8


Chapter 2 Precipitation2021

Figure 2.6.A Sypon (Float) type automatic rain gauge.


Disadvantages of such gauges are: (i) they are costlier than other non-recording type
and (ii) mechanical defects sometimes give erroneous results. There are many
advantages of this type of recording gauge. All types of information about a storm can
be obtained from the cumulative plot of graph paper.

Fig 2.7Used graph paper for siphon type gauge

Rain Gauge Network Design


Lecture note Civil Engineering9
Chapter 2 Precipitation2021

Ideally, a basin should have as many numbers of gauges possible to give a clear
representative picture of the aerial distribution of the precipitation. Factors like economy,
topography, accessibility and rainfall variability govern the number of stations for a
basin. There is no definite rule as to how many gauges are needed for a complete
ungauged basin. To begin with a few pilot gauges can be fixed and after a few years of
data are available, statistical analysis can be carried out to check the adequacy of the
system. World Meteorological Organization (WMO) recommends certain densities of
gauges to be followed for different types of catchments.
Table 1.Israel has the highest network density in the world with one gauge
Regions Minimum area for one Area to be covered
station under ideal under difficult condition
condition in sq. km. per station in sq. km.
1. Flat regions of temperate 600-900 900-3000
Mediterranean and tropical zones.
2. Mountainous regions of temperate 100-250
Mediterranean and tropical zones. 250-1000
3. Small mountainous regions with 25
irregular precipitation.
4. Arid and polar Zones. 1500-10000
Optimum Number of Rain Gauge Stations
Records from all the existing gauges of a basin help to fix the optimum number of
stations. The following statistical analysis helps to obtain number of gauges for a basin
optimally on the basis of an assigned percentage of error in estimating the mean aerial
rainfall.
2
Cv
N=
[ ]Ep
WhereN is the optimal number of stations, Ep is the allowable percentage of error in the
estimation of mean aerial rainfall, Cu is the coefficient of variation of the rainfall from the
existing stations in percentage. Coefficient of variation can be calculated in the following
steps from the data of existing n stations.
(i) Calculate the mean of rainfall 1
from the equation;
Pav =
1
P
σ n−1=
√ (n−1)
∑ ( Pi−P av )2

n∑ i (iii) Compute the coefficient of


(ii) Calculate the standard deviation variation as
as; σ ∗100
C v= n−1
Pav

Lecture note Civil Engineering10


Chapter 2 Precipitation2021

If the allowable percent of error in estimating the mean rainfall is taken higher, then a
basin will require fewer numbers of gauges and vice-versa. The allowable percentage of
error Epis normally taken as 10%. While computing the value of Cuif its value comes 4(ss
than 10%, we can assume the existing stations to be sufficient for the basin. In case N
>n, the additional stations required for the basin can be found as (N - n). Annual rainfall
values are normally used in the above analysis. Additional stations are to be established
at the appropriate locations giving an even distribution over the basin.
Ideal location for rain gauge stations
While setting up any rain Gauge station the following points should be noted.
(i) The site should be on a level ground, i.e., slopping ground, hill tops or hill slopes
are not suitable.
(ii) The site should be an open space
(iii) Horizontal distance between the rain gauge and the nearest objects should be
twice the height of the objects.
(iv) Site should be away from continuous wind forces.
(v) Other meteorological instruments and the fencing of the site should maintain the
step (iii) above.
(vi) The site should be easily accessible. (vii) The gauge should be truly vertical.
(vii) Ten percent of total number of rain gauge stations of any basin should be self-
recording.
(viii) The observer must visit the site regularly to ensure its proper readiness for
measurement
Precipitation measurements are susceptible to the following errors, which can be (a) in
measurement, (b) mistakes in recording, (c) instrumental errors, (d) initial loss in wetting
the gauge (which is equal to 0.25 mm per precipitation or may be 25 mm/year), (e) error
due to rain drop splash, (f) evaporation from the gauge and (g) loss of 1.5% of rainfall if
gauge is inclined 10 degree from its true vertical position.
Example 2.1: A sub-basin has six numbers of rain gauges. Annual rainfall recorded by
the gauges is given below. Considering 10% error in the estimation of mean annual
rainfall, calculate optimum number of gauges required for the sub-basin and check if the
present network is sufficient.
Rain gauge Name A B C D E F
Annual rainfall (cm) 102 77 84 53 66 80

Example 2.2: A sub-basin with area of 1038 sq. km has 7 stations. The normal annual
rainfall depths for all the seven stations are given below. Determine the optimum number
of rain gauge stations to be established in the basin if it is desired to limit the error in the
mean value of rainfall to 10%. Indicate how you are going to distribute the additional
rain gauge stations (if required). Is it possible to have zero percent error in the estimate of
the mean value?

Lecture note Civil Engineering11


Chapter 2 Precipitation2021

2.5 Estimating Missing Data


Failure of any rain gauge or absence of observer from a station causes short break in the
record of rainfall at the station. These gaps are to be estimated first before we use the
rainfall data for any analysis. The surrounding stations located within the basin help to
fill the missing data on the assumption of hydrometeorological similarity of the group of
stations. The general equation of the weightage transmission of the rainfall of the nearby
stations to the missing station (xi) can be represented as:
n
∑ ai Pi
Pxi = i=1n
∑ ai
i =1
Where Pi is the normal rainfall of ithsurrounding station, i =1,2,3, ... , n are the
surrounding gauge numbers which are used for filling the gaps, aithe weighing factor of
the station Pi and Pxiis the data required to be filled up. The methods mostly used in
hydrology for filling the missing data are discussed below.
2.5.1 Arithmetic Mean Method
This method is used when: (i) the normal annual rainfall of the missing station say x is
within 10% of the normal annual rainfall of the surrounding stations, (ii) d~ta of at least
three surrounding stations, called index stations, are available within the basin, (iii) the
index stations should be evenly spaced around the missing station and should be as close
as possible, (iv) The missing rainfall data of station x is computed by simple arithmetic
average of the index stations in the form
1
Px = P +P +. ..+P n )
n( 1 2
In which P1, P2, . . ., Pn are the precipitations of index stations and Pxis that of the missing
station, n is the number of index stations. This is the simplest form of equation where
a1+a2+. . . +an= 1. The word normal means average of 30 years of data, i.e., 30 values of the
latest records. For example, for a station when the last 30 years of June month rainfall is
averaged, we call it as normal rainfall for the month of June for that station.
2.5.2 Normal Ratio Method
This method is used when the normal annual precipitation of the index stations differ by
more than 10% of the missing station. The rainfall of the surrounding index stations are
weighed by the ratio of normal annual rainfalls by using the following equation.
N P1 P 2 P
Px = x
(
n N1 N2
+ +.. .+ n
Nn )
where P1, P2, . . ., Pnare the rainfall data of index stations, N 1, N2, . . ., Nnthe normal annual
rainfall of index stations, Pxand Nxthe corresponding values for the missing station x in
question and n is the number of stations surrounding the station x. Comparing the above
equation with the general equation, the weight coefficient ai for the ithstation is given as;

Lecture note Civil Engineering12


Chapter 2 Precipitation2021

Nx
ai =∑
nN i
2.5.3 Regression Method
A multiple linear regression of the form
Px =a o +a1 P1 +a 2 P 2 +. ..+a n Pn may be established.
The coefficients ao, a1, a2, . . ., an can be calculated by least square method. The equation
can be used to compute rainfall Px of the missing station. Use of this method is
advantageous over the previous two and can be effectively used when a digital computer
is available at a site. A random error component ε ×Syxmay be added to equation above if
large amount of missing data is to be estimated, where EI is the normal random number
with zero mean and unit standard deviation. It can be selected from random number
tables given in any standard mathematical book. Syxis the standard error of estimate. The
standard error of estimate can be estimated from the equation
n
S yx =
√∑ (i=1
Y i−Y ei )2

in which yeiis the estimated value of Yifor given Xi and n is the number of data. Addition
of the last term Et x SyXmaintains the standard deviation of the estimated value of Pxclose
to the observed standard deviation. Use of the above method is not popular because of the
nature of computations involved.
2.5.4 Inverse Distance (US Weather Service) Method
In this method, a rectangular coordinate system is superimposed over the map marked
with rain gauge station in such a way that the origin (0, 0) represents the missing station
as shown in Fig. 2.8. The surrounding index stations lie within the quadrants to the point
for which rainfall is to be estimated. The distance of index stations from the missing
station gives a weightage of the station by which missing rainfall is estimated. The
following relation may be used.
n
∑ W i Pi
Px = i=1n
∑ Wi
i=1
1
W i= 2
Where: D ; D2 =(X2+Y2) is the distance of the station I in X and Y
coordinates taking missing rainfall station at (0,0) position. This is the most acceptable
method and is widely used for determining the missing rainfall for any scientific analysis.
However, the limitation is that it estimates missing rainfall between the highest and
lowest values of the index stations.
Other methods using various numerical interpolations of Isohyetal method, Thiessen
polygon method, Station-year-method, Graphical method and Rational method have their
own limitations for computation of the missing data. Estimation of missing daily rainfall

Lecture note Civil Engineering13


Chapter 2 Precipitation2021

data is less reliable than annual values and such estimates should be subjected to further
statistical analysis. In plain terrain any of the above methods can be suitably used but in
hilly regions, the normal ratio method or linear regression method can be more
effectively used for better results.

Fig. 2.8.
Example 2.4: In a river basin, a station A was inoperative during a storm, while stations
B, C and D surrounding A were in operation, registering 12.3, 14.8 and 11.9 cm of
precipitation. Mean annual precipitation at the four stations A, B, C and D are 1290,
1510, 1680 and 1375 mm respectively. Estimate the missing storm precipitation of station
A by all the methods you know. The coordinates of B, C and D are (6, -4), (8, -6) and (-
4, 4), respectively, y where as the coordinate of A is (0, 0) as shown in Fig. 3.12.
Example 2.5: Normal rainfall of 6 rain gauge stations A1, A2, A3, A4, A5 and A6 are
122, 98, 72, 116, 135 and 110 cms respectively. During a particular storm station A2 was
inoperative due to some mechanical defect. The locations of the stations on a map are
given in Fig. 2.9 and bracket terms give the rainfalls. Calculate the missing rainfall of
station A2 by all the methods and compare the results.

Fig 2.6

Lecture note Civil Engineering14


Chapter 2 Precipitation2021

2.6 Consistency of Rainfall Data


Rainfall data reported from a station may not be consistent always. Over the period of
observation of rainfall record, there could be (i) unreported shifting of the rain gauge site
by as much as 8 km aerially or 30 m in elevation, (ii) significant construction work in the
area might have changed the surroundings, (iii) change in observational procedure
incorporated from a certain period or (iv) a heavy forest fire, earthquake or landslide
might have taken place in that area. Such changes at any station are likely to affect the
consistency of data from a station. Use of double mass-curve checks the consistency of
the record and helps to correct the rain gauge data for the station. In this method, the
accumulated annual rainfall of a particular station is compared with the concurrent
accumulated values of mean rainfall of groups of 5 to 8 surrounding base stations. The
basis of such an exercise is that a group of sample data (for any period) drawn from its
population will be the same. The procedure for double mass-curve analysis is discussed
as follows:
(i) The doubtful station, say A, is marked and the group of stations surrounding it are
identified.
(ii) A table is prepared in which the first column represents the year in decreasing order,
i.e., it starts with the latest year of station A.
(iii) Yearly precipitation values of station A are written in second column.
(iv)In the third column the cumulative rainfall of second column are entered.
(v) Mean yearly precipitation of the group of stations surrounding station A are computed
and entered in the fourth column against the year of col. 1.
(vi)In column five, cumulative precipitation of the group of stations of column four are
computed.
(vii) A graph is plotted taking the cumulative rainfall of the group of station as abscissa
and cumulative rainfall of the station A as ordinate. Consecutive points are joined by
straight line.
(viii) If the consistency of the station A has undergone changes from any year, then
it can be noticed from the slope of the plot. The lines joining the initial points of the
graph A are extended by a dotted line and correction is computed.
(ix)Rainfall of subsequent years from the year of deviation (marked x in the figure) are
corrected by multiplying the correction factor. This exercise helps to bring the older
rainfall data of station A to the new environment. Corrections are to be applied to the
data of station A only when the change of slope of the double mass-curve is observed
for more than five years.

Lecture note Civil Engineering15


Chapter 2 Precipitation2021

Example 2.5: Annual rainfall of station A and the average annual rainfall of five
surrounding stations from 1996 to 1977 are given below. Check the consistency of data
of station A. If data is found inconsistent, then correct the inconsistent data.

Cum. Avr. Ppt Cum Avr. Ppt of


Year Avr.Ppt of A of A Avr. Ppt of group group
1996 1430 1430 1410 1410
1995 1100 2530 1260 2670
1994 1170 3700 1100 3770
1993 1100 4800 1230 5000
1992 1200 6000 1150 6150
1991 1220 7220 1430 7580
1990 1280 8500 1150 8730
1989 750 9250 950 9680
1988 1120 10370 1230 10910
1987 1250 11620 1350 12260
1986 1380 13000 1440 13700
1985 1210 14210 1360 15060
1984 1760 15970 1730 16790
1983 1400 17370 1080 17870
1982 1240 18610 970 18840
1981 1760 20370 1320 20160
1980 1480 21850 1350 21510
1979 1740 23590 1410 22920
1978 1420 25010 1270 24190
1977 1580 26590 1260 25450

Solution
A graph is plotted taking the cumulative annual rainfall of the five surrounding stations in
abscissa and the cumulative annual rainfall of the station A in ordinate. The consecutive
points are joined by a straight line. Deviation of the straight line plot is noticed from the
year 1983. Slope of the line is found as 0.756. Therefore correction is applied to all the
data from 1983 to 1977.

Lecture note Civil Engineering16


Chapter 2 Precipitation2021

Fig. 2.9. Double Mass Curve

2.7 Presentation of Precipitation Data


Rainfall is usually presented in the form of the following graphs. Such graphs are useful
for analysis and design purposes.
2.7.1 Moving Average Curve
In hydrology, rainfall data are plotted chronologically between time in x-axis and
precipitation in y-axis. A rain event is associated with randomness. To overcome the
random component, a simple moving average of order 3 or 5 is used. This helps to isolate
the trend in the rainfall data. If there is any dry or wet cyclic trend associated with
precipitation, then such a trend can be clearly visible from the plot. From the graph, the
wet period mean, overall mean and dry period mean can be identified. Such a method is
applicable to annual series. If Xl, X2, X3, X4, X5 are annual precipitation at a station and
a 5 year moving average is applied to the series then the following moving mean are
computed
( X + X 2+ X 3+ X 4+ X 5 )
X 1= 1
5
( X 2 + X 3+ X 4 + X 5 + X 6 )
X 2=
Similarly 5
( X 3 + X 4+ X 5+ X 6+ X 7 )
X 3=
and 5
2.7.2 Mass Curve
Such a curve can be presented from third year onwards only. For example, if data are
available from 1961 to 1996, then a 5 year moving average can be represented from 1963
to 1994. The first two years data of 1961-1962, and the last two years data of 1995-1996
Lecture note Civil Engineering17
Chapter 2 Precipitation2021

are lost in the moving average process. In hydrology a moving average of more than 5 is
usually not applied as some of the cyclic trends associated with the data are smoothened
out. A plot of moving average is shown in Fig. 2.10

Fig. 2.10: Moving Average Curve


2.7.3 Rainfall Hyetograph
Rainfall intensity means the ratio of rainfall depth with time and is expressed in cm/hour
(cm/h) or mm/hour (mm/h). From the mass curve, a plot between the intensity of rainfall
with time can be obtained. During a storm, intensity always changes with time. On a
mass curve any two points can be marked and the depth of rainfall (y) between these
two points is noted from the y-axis. Time between these two points (M) are recorded
from x-axis. The depth divided by time, i.e., x/t is the intensity of rainfall for the
period under consideration. Usually dtselected is either 1, 2 or 3 hours. For larger
cyclonic storms spreading over a few days t may be selected for 6 or 12 hours.

Fig. 2.11: Rainfall Hyetograph

Lecture note Civil Engineering18


Chapter 2 Precipitation2021

When the plot of rainfall intensity with time is presented in the form of bar graphs then
such a graph is known as hyetograph (Fig. 2.11). The plot is very useful for flood studies
and calculation of rainfall loss indices.
2.7.4 Intensity-Duration-Frequency Curves
An intensity-duration and frequency curve is a three parameter curve in which duration is
taken on x-axis, intensity on y-axis and the return period or frequency as the third
parameter. By fixing the return period of say 10-, 50-, 100-years or any other period, a
particular curve between intensity and duration can be obtained for the area. Through
such a curve, an exponential equation of the following order can be fit.
4 −d
I=CT ( D+b )

Fig. 2.12 Intensity-Duration-Frequency Curves


in which T is the frequency or the return period of the storm of intensity I cm/h and
duration D-hours. Other terms in the equation are constants. For Indian conditions a
varies between 0.15 and 0.7, b from 0.10 to 1.05, C from 3 to 15 and d from 0.75 to 1.25.
Plot of such a curve is shown in Fig. 3.12.
Example 2.6: Rainfall recorded by a self recording rain gauge for the flood event on 28th
August, 1996 at a station is given below. Construct a rainfall mass curve and plot the
hyetograph for the event.
Time (h) 9 9.5 10 10. 11 11.5 12 12. 13 13.25
5 5
Rainfall 0 0.55 2.25 3.5 3.85 4.25 4.5 4.8 5.3 5.4
(cm) 5

2.8. Mean Aerial Rainfall


A rain gauge records rainfall at a geographical point. In most of the hydrologic analysis,
average depth of precipitation over the area under consideration is required to be
computed on hourly, daily, storm period, ten-day, monthly or yearly basis. There are

Lecture note Civil Engineering19


Chapter 2 Precipitation2021

many methods available in literature for computation of average precipitation over the
basin. However, depending on the accuracy and objective of the analysis any of the
followings methods can be used:

(i) Arithmetic average, (ii) Thiessen polygon, (iii) Isohyetal, (iv) Grid point, (v)
Orographic or (vi) Isopercental method.

Other methods available in the literature include the triangular mean weight, trend
surface analysis, reciprocal distance square, modified polygon, analysis of variance and
the double Fourier series. These methods are not extensively used in hydrology. The
rigorous mathematical calculations involved in these methods do not support the
accuracy achieved in arriving at the mean aerial rainfall over a basin.
2.8.1 Arithmetic Mean Method
This method is suitably applied for a basin where the gauges are uniformly distributed
and the individual gauge catches do not vary much from the mean. The basin should be a
reasonably flat area. The assumption made is that all gauges weigh equally. This method
gives fairly good results if the topographic influences on precipitation and aerial
representativeness are considered while selecting the gauge site. It is the simplest form in
which the average depth of precipitation over the basin is obtained by taking simple
arithmetic mean of all the gauged amounts within the basin
P 1 + P2 +. . .+ P n 1 n
Pav = = ∑ Pi
n n i=1
Where P1, P2 . . . Pn are the precipitation recorded by n number of gauges located within
the basin. This method gives a rough estimate of the average precipitation. It does not
account for the topographic and other influences. For use of this method, no gauge station
located outside the boundary of the watershed should be considered.
2.8.2Thiessen Polygon Method
In this method, weightage is given to all the measuring gauges on the basis of their aerial
coverage on the map, thus eliminating the 'discrepancies in their spacing over the basin.
All the stations in and around the basin are considered and a linear variation in the
precipitation between two gauge stations is assumed. Procedure for delineating area for
each station in the basin is described below.
(a) All the gauges in and around the basin are accurately marked on a map drawn to scale.
(b) Consecutive stations are joined by dotted straight lines, forming triangles.
(c) Perpendicular bisectors are drawn to these dotted lines such that the bisectors form a
polygon around each station. For a basin having large number of stations, the
approach should be to start from one end of the map and traverse to the other end
gradually forming polygons. The boundary of the map and the perpendicular bisectors
cutting them form polygons for the periphery stations are shown in Fig. 2.13, while
for inner stations, the bisectors only from the polygons.

Lecture note Civil Engineering20


Chapter 2 Precipitation2021

(d) Each station on the map is thus enclosed by a polygon. A polygon represents an area
for which the station rainfall is the representative.
(e) Area of each polygon is measured by planimetering. Sum of the areas of all the
polygons must be equal to the total area of the basin.
(f) Thiessen weights are computed by dividing the area of each polygon by the total area
of the basin. Thus, sum of Thiessen weights for all stations should be equal to unity.
If there are seven stations in and around the basin then seven thiessen polygons are
drawn. Sum of all the seven Thiessen weights must be equal to unity.
(g) The average precipitation is computed from the relation
A P + A P +.. .+ An P n
Pav = 1 1 2 2
A 1 + A 2 +.. .+ An P n

OR
P =W P +W P +. . .+W P
av 1 1 2 2 n n
Where P1, P2, . . . Pn represents precipitation t stations 1, 2, 3, . . . , n and A1, A2, . . . , An
represents the area of polygons representing the corresponding stations, A is the total area
of basin which is the sum of all the polygons and W1, W2, . . . , Wn= An/A such that
W1+W2+. . . +Wn= 100.

Fig. 2.13: Thiessen Polygon Method


This method suffers from the following limitations.
(i) If the network of a basin is changed, i.e., if there is an addition or removal of a
station from the basin, then a new thiessen diagram is to be drawn.
(ii) Orographic features are not accounted for.
(iii) A linear variation of precipitation between two stations is assumed, whereas the
precipitation is influenced by a large number of meteorological and catchment
characteristics.
(iv) Topographic influences and other barriers are not considered.
However the advantage of this method is that, it is much more accurate than the previous
methods and the procedure of computation becomes simple, once the areas of the
Lecture note Civil Engineering21
Chapter 2 Precipitation2021

polygons are measured. This method is popularly applied to 'l11ost of the field problems.
2.8.3 Isohyetal Method
This method gives more accurate results of the average rainfall of a basin. An .
experienced analyst takes care of the orographic features and storm characteristics while
drawing contours of equal rainfall depths of a basin. The resulting map can represent the
actual rainfall pattern of the storm over the watershed. An Isohyet is a line joining points
of equal rainfall magnitude on a map. Steps in determining the average precipitation of
an area are given below.
(i) A map of the basin is drawn to scale.
(ii) All gauge stations in and around the basin are accurately located on the map.
(iii) Depth of precipitation recorded at each station are marked on the map.
(iv) Isohyets are drawn by eye approximation interpolating the distances between
stations on consideration of orographic, storm characteristics and other factors
that affect the rainfall variability at the place. They follow the principles of
elevation contours drawn on a map. Isohyets do not cross each other. Over a large
area, all isohyets form closed contours. Closely spaced isohyets’ contours indicate
that the precipitation gradient is more. Area bounded by the highest closed
contour represents the eye of a cyclonic storm.
(v) Area between successive isohyets within the basin is measured. (planimetered).
(vi) Average precipitation between two successive isohyets multiplied by the area
bounded by them should be computed for all isohyets covering the area.
(vii) Sum of all such products over the entire basin divided by the total area of the
basin gives average precipitation. This can be computed using the following
equation.
P +P P +P P +P
A 1 1 2 + A 2 2 3 +. ..+ A n n−1 n
2 2 2
Pav =
A 1 + A 2 +. ..+ A n

Lecture note Civil Engineering22


Chapter 2 Precipitation2021

Fig. 2.14: Isohyetal Method


Where P1, P2, ••• , Pnare the isohyetalvalues such that P1and P2bound the area A1, P2and
P3bound the area A2and so on. Sum of all such sub-areas bounded by all isohyets, i.e., Al+
A2+ ... + Anis the total area of the basin as shown in Fig. 2.14. If a linear interpolation
between stations is used while drawing isohyets, then the result of Thiessen-Polygon
method and this method will essentially be the same. Isohyetal method gives more
accurate results than the previous two. To achieve this experienced hydrologist or the
engineer should draw the isohyets with due consideration of the topographical and other
influences. Thorough knowledge of the topography of the basin helps to draw isohyets
more accurately. Demerit of such a system is that for each storm a separate isohyet is to
be drawn. If monthly rainfall is to be computed by isohyetal method then for each month
of each year one set of isohyet is to be drawn thus multiplying the work many fold. Since
this method permits consideration of topographic features, it is invariably used to
calculate the average precipitation
during storms for flood studies of
various projects.
2.5.4 Inverse-distance method
This technique is based on a
rectangular grid covering the
reference area (Figure 2.5). For each
gridpoint the precipitation value is
computed using the neighboring gage
site observations. The quadrant
method is most commonly used (Figure
2.6). A pair of coordinate lines facing
North- South and East-West is placed
upon each gridpoint. The gridpoint precipitation value is computed from the four closest

Lecture note Civil Engineering23


Chapter 2 Precipitation2021

gage stations, one in each quadrant. Areal precipitation may then be determined by
computing the arithmetic mean using all gridpoint precipitation data.

Figure 2.15: Rectangular grid covering the reference area

Figure 2.16: Quadrants for gridpoint R


Precipitation for each gridpoint may be determined by

hP,j[mm] height of precipitation at gridpoint j


Wi,j[-] weight of gridpoint j with respect to station i
hP,i[mm] height of precipitation at station / element i

Lecture note Civil Engineering24


Chapter 2 Precipitation2021

where the weights Wi,j are the relative reciprocal square distances between the gridpoint
and the respective gage station.

Wi,j [-] weight of gridpoint j with respect to station i


di,j [m] distance of station i from gridpoint j
Areal precipitation may be estimated by:

hP [mm] areal precipitation


nG [-] number of gridpoints
hP,j [mm] height of precipitation at gridpoint j

Lecture note Civil Engineering25


Chapter 2 Precipitation2021

2.8.6 Isopercental Method


This method gives most accurate estimate of the average rainfall when the catchment is
affected by orographic and other topographical features. The method is similar to the
isohyetal method of estimation of catchment average rainfall. The procedure to calculate
the average storm of say monthly rainfall Ri over the catchment are outlined below.
(i) Obtain a map of the basin to scale showing the positions of all rain gauge stations.
(ii) Collect the storm or monthly rainfall data Ri and the normal annual rainfall Ni
for the marked stations.
(iii) Convert the storm or monthly rainfall values Ri as the percentage of the
normal annual values Ni using the relation Xi= 100Ri/Ni.
(iv) Draw the isopercental lines, i.e., lines joining the percentage of same rainfall over
the basin as thick lines by selecting a suitable interval of isopercentages. This can
be drawn following the rules of drawing of the contours.
(v) On the same map draw isohyetallines by considering the normal annual rainfalls
over the basin. Let these lines be shown as dotted lines.
(vi) Locate the points where the isopercentallines cut the annual isohyets. (vii)
Calculate the rainfall at the points of intersection of step (6). This can be done by
multiplying the values of (Xnli /1(0) and (N i) in the units of depth of rainfall,
where Xnli is the value of Xi at the paints of intersection
(vii) Now draw the isohyetal patterns of the area on a separate map by considering
these intersection points and their depth of rainfall.
(viii) Obtain the catchment rainfall over the basin by considering the isohyets as per
the procedure outlined in the Section 2.7.3.
The method though gives the most promising alternative of the calculation of average
catchment rainfall in a mountainous regions, the method lacks its popularity due to
repetitive and difficult nature of computations involved. The data requirement for this
method is also more.
Example 2.7: For a catchment of 590 sq. km rainfall for the month of July along with
their areas of influence for the stations are given below. Calculate average precipitation
over the basin by (a) Thiessen-polygon method, (b) Isohyetal method.
Station A B C D E F G H I J K L M
RF (cm) 37.5 1 3 4 5 52.5 56 32.5 4 3 5 3 45
9 5 5 3 0 0 4 2
2
Area of influence (km ) 3 8 9 1 6 79 88 70 6 3 6 5 10
(As per Thiessen-Polygon) 1 5 9 0 6 5

Isohyet <20cm 20-30cm 30-40cm 40-50cm 50-60cm


area (km2) 3 44 82 242 219

Lecture note Civil Engineering26


Chapter 2 Precipitation2021

2.9 Depth –Area Duration (Dad) Curve


Depth of precipitation of a storm is related to the area of its coverage and duration of the
storm. DAD analysis is carried out to obtain a curve relating the depth of precipitation D,
area of its coverage A and duration of occurrence of the storm D. A DAD curve is a
graphical representation of the gradual decrease of depth of precipitation with the
progressive increase of the area of the storm, away from the storm center, for a given
duration taken as the third parameter. This gives a direct relationship between depth, area
and duration of precipitation over the region for which the analysis is carried out. The
purpose of DAD analysis is to determine the maximum precipitating amounts that have
occurred over various sizes of drainage area during the passage of storm periods of say 6
h, 12 h, 24 h or other durations. There are two methods available for carrying out DAD
analysis. The first method, known as mass curve method is not extensively used by
hydrologists. The second method called the incremental-isohyetal method is more
popular. Steps in carrying out DAD analysis of a region by the second method are as
follows:
(i) All the major storms of the area are identified.
(ii) The duration for all the storms are noted (if duration chosen is say I-day, then all
storm's occurring for one day are selected. When a storm has occurred for 3 days,
then the maximum one day precipitation out of the 3 days is noted).
(iii) Isohyetal patterns for all 1 day storms are prepared on maps.
(iv)One day storm is taken up and the area bounded within the highest isohyet is
planimetered. This is called the eye-area of the storm. Next the area bounded
between the largest and the second largest isohyet is planimetered. The depth of
precipitation covering up to the second largest isohyet is obtained by the relation d2=
(Pm1A1+ Pm2A2)/(A1 + A2), where Pml is mean precipitation bounded by the highest
isohyetal area A1and Pm2is the mean precipitation between the first and second
highest isohyets covering area A2, Similarly for the area covering up to the third
isohyet, the progressive depth of precipitation can be obtained by the relation d3=
(Pm1A1+ Pm2 A2+ Pm3)/(A1 +A2+ A3) where Pm3is the mean precipitation between the
second and third highest isohyet covering area A3between them. The procedure is
repeated to cover the remaining isohyets of the area.
(v) All the area-depth precipitations are recorded in a table.
(vi)The steps from (iv) to (v) are repeated for all other 1day storm experiences of the
area.
(vii) A graph is plotted between area as abscissa and maximum average depths of
precipitation as ordinate covering the depth-area data of all 1 day storms of step (vi).
(viii) Such an exercise may also be taken up for 6-h, 12-h, 2- and 3-day storms of the
region. The curves are plotted on the same paper as in step (vii) or separate plots
may be made for each durations.

Lecture note Civil Engineering27


Chapter 2 Precipitation2021

(ix)If a semi-log graph paper is used with area plotted on log scale then the curve will
plot close to a straight line.
The resulting plot is the DAD curve for the region. DAD analysis requires all the
previous storm information and involves tremendous computational efforts. Use of
computers reduces the work load to a great extent. India Meteorological Department on
request can carry out DAD analysis for any portion of the country. Depending on the
availability of data DAD analysis for any region and for any other periods can be
prepared. DAD curves for 12-h, 1 and 2-day storms are shown in Fig. 2.15.
Interesting points to note from the graph are (a) when the area of the storm increases, the
depth of precipitation decreases (b) when the duration of the storm increases for the given
area, the depth of precipitation increases. Use of DAD graph is further discussed in the
flood chapter.
Example 2.8: From the isohyetal map shown in Fig. 2.14and the sub-areas between the
isohyets given in Example 2.7, prepare a DAD curve for the basin if the rain depths are
recorded from a 2-day storm.

Lecture note Civil Engineering28

You might also like