Chapter 3 - Part 1
Chapter 3 - Part 1
Chapter 3 - Part 1
Hydrometeorology
Prepared By:
NATHANIEL R. ALIBUYOG
August 2020
Module 3
Precipitation
Introduction
Precipitation replenishes surface water bodies, renews soil moisture for plants, and recharges aquifers. Its
principal forms are rain and snow. The relative importance of these forms is determined by the climate of
the area under consideration. Hydrologic modeling and water resources assessments depend upon a
knowledge of the form and amount of precipitation occurring in a region of concern over a time period of
interest.
Learning Outcomes
Learning
Input
The conditions for precipitation to take place may be summarized stepwise as follows:
1. Supply of moisture
2. Cooling to below point condensation
3. Condensation
4. Growth of particles
The supply of moisture is obtained through evaporation from wet surfaces, transpiration from
vegetation or transport from elsewhere. The cooling of moist air may be through contact with a cold
earth surface causing dew, white frost, mist or fog, and loss of heat through long wave radiation (fog
patches). However, much more important is the lifting of air masses under adiabatic conditions
(dynamic cooling) causing a fall of temperature to near its dew point.
1. Convection, due to vertical instability of the air. The air is said to be unstable if the temperature
gradient is larger than the adiabatic lapse rate. Consequently, a parcel moving up obtains a
temperature higher than its immediate surroundings. Since the pressure on both is the same
the density of the parcel becomes less than the environment and buoyancy causes the parcel
to ascend rapidly. Instability of the atmosphere usually results from the heating of the lower air
layers by a hot earth surface and the cooling of the upper layers by outgoing radiation.
Convective rainfall is common in tropical regions and it usually appears as a thunderstorm in
temperate climates during the summer period. Rainfall intensities of convective storms can be
very high locally; the duration, however is generally short.
2. Orographic lifting. When air passes over a mountain it is forced to rise which may cause
rainfall of the windward slope. As a result of orographic lifting rainfall amounts are usually
highest in the mountainous part of the river basin.
3. Frontal Lifting. The existence of an area with low pressure causes surrounding air to move
into the depression, displacing low pressure air upwards, which may then be cooled to dew
point. If cold air is replaced by warm air (warm front) the frontal zone is usually large and the
rainfall of low intensity and long duration. A cold shows a much steeper slope of the interface
of warm and cold air usually resulting in rainfall of shorter duration and higher intensity.
4. Cyclones, tropical depressions or hurricanes. These are active depressions, which gain
energy while moving over warm ocean water, and which dissipate energy while moving over
land or cold water. They may cause torrential rains and heavy storms. Typical characteristics
of these tropical depressions are high intensity rainfall of long duration (several days).
5. Convergence. The Inter-tropical Convergence Zone (ITCZ) is the tropical region where the
air masse originating from the Tropic of Cancer and Capricon converge and lift. In the tropics,
the position of the ITCZ governs the occurrence of wet and dry seasons. In July, the ITCZ lies
to the north of the equator and in January it lies to the south. In the tropics the position of the
ITCZ determines the main rain-bringing mechanism which is also called monsoon. Hence, the
ITCZ is also called the Monsoon Trough.
Estimate of the amount of precipitation that might occur over a given region with favorable conditions
are often useful. These may be obtained by calculating the amount of water contained in a column of
atmosphere extending up from the earth’s surface. This quantity is known as the precipitable water,
W, although it cannot all be removed from the atmosphere by natural processes.
The amount of precipitable water in the atmosphere is equal to the water mass contained in a column,
having a square base of 1 cm, between elevation zero and some height z. In equation form, it is given
as
"
𝑊 = ∫# 𝜌! 𝑑𝑧
3. Cyclonic Precipitation. 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 annual heating of the earth’s surface. Cyclonic precipitation may be classified as frontal
and nonfrontal.
a. Nonfrontal precipitation - is produce when air is lifted through horizontal convergence
of the inflow into a low-pressure area.
b. Frontal precipitation - results from the lifting of warm air over cold air at the contact zone
between air masses having different characteristics.
c. Warm front – if the air masses are moving so that warm air replaces colder air
d. Cold front - if the cold air displaces warm air
e. Stationary front - if the front is not in motion.
4. Thunderstorms. Thunderstorm cells develop from vertical air movements associated with
intense surface heating or orographic effects.
2. Snow - Precipitation in the form of ice crystals which usually combine to form flakes, with an
intensity density of 0.1 g/cm3
3. Drizzle - Rain droplets of size less than 0.5 mm and rain intensity of less than 1 mm/h is known
as drizzle
4. Glaze - When rain or drizzle touches ground at 0oC, glaze or freezing rain is formed
5. Sleet - It is frozen raindrops of transparent grains which form when rain falls through air at
subfreezing temperature
6. Hail - It is a showery precipitation in the form of irregular pellets or lumps of ice of size more
than 8 mm
The instrument used to collect and measure the precipitation is called rain gauge.
Hyetograph – plot of rainfall intensity against time. The hyetograph is derived from the mass curve and is
usually represented by a bar chart. It is a very convenient way of representing the characteristics of a storm
and is particularly important in the development of design storms to predict extreme floods.
Mass Curve of Rainfall -plot of accumulated precipitation against time, plotted in chronological order.
Mass curves of rainfall are very useful in extracting the information on the duration and magnitude of a
storm. Also, intensities at various time intervals in a storm can be obtained by the slope of the curve.
Point Rainfall – It is also known as station rainfall. It refers to the rainfall data of a station.
Moving average. Moving average is a technique for smoothening out the high frequency fluctuations of a
time series and to enable the trend, if any, to be noticed. The basic principle is that a window of time range
m years is selected. Starting from the first set of m years of data, the average of the data of m years is
calculated and placed in the middle year of the range m. The window is next moved sequentially one time
unit (year) at a time and he mean of m terms in the window is determined at each window location. The
value of m can be 3 or more years; usually ad odd value.
Hyetograph Mass Curve of Rainfall
For most hydrologic analyses, it is important to know the areal distribution of precipitation. The reliability
of rainfall measured at one gauge in representing the average depth over a surrounding area is a function
of:
a. The distance from the gauge to the center of the representative area
b. The size of the area
c. Topography
d. The nature of the rainfall of concern
e. Local storm pattern characteristics
The following methods are used to measure the average precipitation over an area:
1. Arithmetic mean method
2. Theissen polygon method
3. Isohyetal method
4. Inverse distance weighting
Arithmetic Mean Method – Simplest method for determining areal average. In this method, only stations
within the watershed boundary are considered.
!
𝑃" = " ∑"
#$! 𝑃#
Thiessen Polygon Method - In this method the rainfall recorded at each station is given a weightage on
the basis of an area closes to the station. The procedure of determining the weighing area is as follows:
Consider a catchment area as shown below containing three raingauge stations. There are three stations
outside the catchment but in its neighborhood. The catchment area is drawn to scale and the positions of
the six stations marked on it. Stations 1 to 6 are joined to form a network of triangles. Perpendicular
bisectors for each of the sides of the triangle are drawn. These bisectors form a polygon around each station.
The boundary of the catchment, if it cuts the bisector is taken as the outer limit of the polygon. Thus, for
station 1, the bounding polygon is abcd. For station 2, kade is taken as the bounding polygon. These
bounding polygons are called Thiessen polygons. The areas of the polygons are determined either with a
planimeter or by using an overlay grid.
∑ "
&! '!
𝑃" = !#$
∑'
The Thiessen polygon method is superior to the arithmetic average method as some weightage is given to
the various stations on a rational basis. Further, the raingauge stations outside the catchment are also used
effectively.
Isohyetal Method. An isohyets 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 𝑃" is to be determined are then marked on the plot at appropriate stations. Neighboring stations
outside the catchment are also considered. The isohyets of various values are then drawn by considering
point rainfalls as guide and interpolating between them by the eye. The procedure is similar to the drawing
of elevation contours based on spot levels.
The area between two isohyets are then determined with a planimeter. If the isohyets go out of the
catchment, the catchment boundary is used as the bounding line. The average value of the rainfall indicated
by two isohyets is assumed to be acting over the inter-isohyet area. Thus, P1, P2, …, Pn are the values of
isohyets and if A1, A2, …, An-1 are the inter-isohyets areas respectively, then the mean precipitation over
the catchment of area A is given by
% &% % &% % &%
'$ ( $ ' )*'' ( ' ()*⋯*')*$ ( )*$ ))
𝑃" = ' '
'
'
The Isohyet method is superior to other two methods especially when stations are large in number.
Inverse DIstance Weighing (IDW) Method. Prediction at a point is more influenced by nearby
measurements than that by distant measurements. The prediction at an ungauged point is inversely
proportional to the distance to the measurement points.
Steps:
a. Computed distance (di) from ungauged point to all measurement points
𝑃
∑-#$! . #, /
𝑑#
𝑃" =
1
∑-#$! 0 2
𝑑# ,
Station 1 2 3 4 5
Coordinates (30,80) (70,100) (100,140) (130,100) (100,70)
Precipitation 85.0 135.2 95.3 146.4 102.2
(cm)
Thiessen polygon
Solution:
The catchment area is drawn to scale and the stations are marked on it. The stations are joined to form a set
of triangles and the perpendicular bisector of each side is then drawn. The Thiessen-polygon area enclosing
each station is then identified. It may be noted that Station 1 in this problem does not have any area of
influence in the catchment. The areas of various Thiessen polygons are determined either by a planimeter
or by placing an overlay grid.
Solution:
For the first area consisting of a station surrounded by a closed isohyet, a precipitation value of 12.0 cm is
taken. For all other areas, the mean of two bounding isohyets are taken.
Before using the rainfall records of a station, it is necessary to first check the data for continuity and
consistency. The continuity of a record may be broken with missing data due to many reasons such as
damage or fault in a raingauge during a period. The missing data can be estimated by using the data of the
neighboring stations. In these calculations the Normal Rainfall is used as standard 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 are recomputed every decade. Thus, the term Normal annual precipitation at Station
A means the average annual precipitation at A based on a specified 30-year of record.
Given the annual precipitation values, P1, P2, P3, …, Pm at neighboring 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.
If the normal annual precipitations at various stations are within about 10% of the normal annual
precipitation at station X, then a simple arithmetic average procedure is followed to estimate Px. Thus
!
𝑃. = / [𝑃! + 𝑃, + ⋯ + 𝑃0 ]
If the normal precipitation vary considerably, then Px is estimated by weighing the precipitation at the
various stations by the ratios of normal annual precipitations. This method, known as the Normal Ratio
Method, gives Px as
Example 3. The normal annual rainfall at station A, B, C, and D in a basin are 80.97, 67.59, 76.28, and
92.01 cm respectively. In the year 2019, the station D was inoperative and the stations A, B, and C recorded
annual precipitation of 91.11, 72.23, and 79.89 cm, respectively. Estimate the annual rainfall at station D
in that year.
Solution:
As the normal rainfall values vary more than 10% , the normal ration method is adopted. Hence,
𝑃. = 99.48 𝑐𝑚
If the conditions relevant to the recording of a raingauge station have undergone a significant change during
the period of record, inconsistency would arise in the rainfall data of that station. This inconsistency would
be felt from the time the significant change took place.
The checking for inconsistency of a record is done by the Double Mass Curve Technique. This technique
is based on the principle that when each recorded data comes from the same parent population, they are
consistent.
Steps:
𝑀<
𝑃<. = 𝑃. > @
𝑀9
Example 4. Annual rainfall data for station M as well as the average annual rainfall values for a
group of ten neighboring stations located in a meteorologically homogenous region are given
below.
Test the consistency of the annual rainfall data of station M and correct the record if there is any
discrepancy. Estimate the mean annual precipitation at station M.
Solution:
a. The data is sorted in descending order of the year, starting from the latest year 1979.
b. Cumulative value of station M rainfall (∑ 𝑃0 ) and the ten station average rainfall values
(∑ 𝑃9:; ) are calculated as shown in Table below.
c. The data is then plotted with ∑ 𝑃0 on the y-axis and ∑ 𝑃9:; on the x-axis to obtain a double
mass curve plot as shown below
d. The value of the year corresponding to the plotted points is also noted on the plot.
e. It is that the data plots as two straight lines with a break of grade at the year 1969. This
represent a change in the regime of the station M after the year 1968.
f. The slope of the best straight line for the period 1979 – 1968 is Mc = 1.029 and the slope of
the best straight line for the period 1968 – 1950 is Ma = 0.8779
g. The correction ratio to bring the old record (1950 – 1968) to the current (post 1968) regimes
is = Mc/Ma = 1.029/0.8779 = 1.173
h. Each of the pre 1960 annual rainfall value is multiplied by the correction ration of 1.173 to get
the adjusted value.
i. The adjusted values at station M are shown in Col. 5 of Table. The finalized values of Pm
(rounded off to nearest mm) for all the 30 years of record are shown in Col. 7
j. The mean annual precipitation at station M (based on the corrected time series) = (19004/30)
= 633.5 mm
1 2 3 4 5 6 7 8
YEA Pm Pave SPm SPave Adjusted Finalized Adjusted
R Values of Pm values of Cumulative
Pm value of Pm
1979 612 588 612 588 612 612
1978 426 410 1038 998 426 1038
1977 825 787 1863 1785 825 1863
1976 685 653 2548 2438 685 2548
1975 356 377 2904 2815 356 2904
1974 568 570 3472 3385 568 3472
1973 438 390 3910 3775 438 3910
1972 386 400 4296 4175 386 4296
1971 497 490 4793 4665 497 4793
1970 635 590 5428 5255 635 5428
1969 375 350 5803 5605 375 5803
1968 596 646 6399 6251 698.92 699 6502
1967 573 650 6972 6901 671.95 672 7174
1966 999 1140 7971 8041 1171.51 1172 8345
1965 1244 1400 9215 9441 1458.82 1459 9804
1964 679 770 9894 10211 796.25 796 10600
1963 828 950 10722 11161 970.98 971 11571
1962 504 580 11226 11741 591.03 591 12162
1961 531 600 11757 12341 622.70 623 12785
1960 415 480 12172 12821 486.66 487 13272
1959 503 575 12675 13396 589.86 590 13862
1958 493 560 13168 13956 578.13 578 14440
1957 431 490 13599 14446 505.43 505 14945
1956 479 540 14078 14986 561.72 562 15507
1955 699 800 14777 15786 819.71 820 16327
1954 472 540 15249 16326 553.51 554 16880
1953 462 520 15711 16846 541.78 542 17422
1952 95 110 15806 16956 111.41 111 17533
1951 578 660 16384 17616 677.81 678 18211
1950 676 780 17060 18396 792.73 793 19004
Total of Pm = 19004 mm
Mean of Pm = 633.5 mm
30000
Cumulative annual rainfall at station M (mm)
Old
25000
Adjusted
20000 New
Linear (Old)
15000
Linear (New)
10000 y = 0.8779x + 917.93
5000
y = 1.0295x + 12.48
0
0 5000 10000 15000 20000
Cumulative ten station average (mm)