Water Resources Engineering PDF
Water Resources Engineering PDF
HYDROLOGIC CYCLE
1.1 Introduction
1.2 Components of Hydrologic Cycle
1.3 Scales for study of hydrologic cycle
1.4 Mathematical Representation of the Hydrologic Cycle
1.5 Global water balance
1.6 Influence of Human Activities and Land use Changes on Hydrologic Cycle
1.7 Impact of climate changes on the hydrologic cycle
1.8 Closure
Keywords: Hydrologic cycle, precipitation, water balance, components, water balance, human
impact, climate change,
1.1 Introduction
Water can occur in three physical phases: solid, liquid, and gas and is found in nature in all these
phases in large quantities. Depending upon the environment of the place of occurrence, water
can quickly change its phase.
A number of cycles are operating in nature, such as the carbon cycle, the nitrogen cycle,
and several biogeochemical cycles. The Hydrologic Cycle, also known as the water cycle, is one
such cycle which forms the fundamental concept in hydrology. Hydrologic cycle was defined by
the National Research Council (NRC, 1982) the as “the pathway of water as it moves in its
various phases to the atmosphere, to the earth, over and through the land, to the ocean and back
to the atmosphere”. This cycle has no beginning or end and water is present in all the three states
(solid, liquid, and gas). A pictorial view of the hydrological cycle is given in Fig. 1.1. The
science of hydrology primarily deals with the land portion of the hydrologic cycle; interactions
with the oceans and atmosphere are also studied. NRC (1991) called the hydrologic cycle as the
integrating process for the fluxes of water, energy, and the chemical elements.
The hydrologic cycle can be visualized as a series of storages and a set of activities that
move water among these storages. Among these, oceans are the largest reservoirs, holding about
97% of the earth’s water. Of the remaining 3% freshwater, about 78% is stored in ice in
Antarctica and Greenland. About 21% of freshwater on the earth is groundwater, stored in
sediments and rocks below the surface of the earth. Rivers, streams, and lakes together contain
less than 1% of the freshwater on the earth and less than 0.1% of all the water on the earth.
Hydrologic
H cycle
c consid
ders the mo otion, loss, and rechargge of the eearth's wateers. It
connects the atmosphere and tw wo storages ofo the earth system: the oceans, andd the land spphere
(lithosph
here and ped dosphere). The
T water ev vaporated frrom the eartth and the ooceans enterrs the
atmospheere. Water leeaves the atmosphere th hrough preciipitation. Thhe oceans receive water from
the atmoosphere by means
m of preecipitation and
a from thhe land throuugh rivers aand ground w water
flow. Water
W goes out
o of ocean ns only thro ough evapooration. The water leavves land thrrough
evapotrannspiration, streamflow,
s and ground water flow. Evaporationn and precippitation proccesses
take placce in the veertical plane while streaamflow and ground watter flow occcur mostly iin the
horizontaal plane.
The
T exchangee of water amonga the oceans, land,, and the atm mosphere waas termed ass ‘the
turnover’’ by Shikhlo omanov (199 99). This turnnover affectts the global patterns of the movemeent of
ocean waaters and gases in the attmosphere, thereby
t greaatly influenccing climate. Since wateer is a
very goood solvent, chemistry is an integral part
p of the hhydrologic ccycle. Usuallly, rain and snow
are consiidered as thee purest form
m of water allthough thesse may also bbe mixed wiith pollutants that
are preseent in the atm
mosphere. Duuring the jouurney on eartth, many cheemical comppounds are m mixed
with watter and conssequently thee water quallity undergooes a changee. One can aalso visualizze the
hydrologgic cycle as a perpetual distillation
d and
a pumpingg system in which the gglaciers and snow
packs aree replenishedd and rivers get water off good qualitty.
We need to study the hydrologic cycle since water is essential for survival of life and is
an important input in many economic activities. From the use point of view, the land phase of the
hydrologic cycle is the most important.
In view of the complexities and extensive coverage, the study of the complete hydrologic
cycle is truly interdisciplinary. For instance, the atmospheric part is studied by meteorologists,
the pedospheric part by soil scientists, the lithosphere part by geologists, and the part pertaining
to oceans falls in the domain of oceanographers. A host of other professionals study hydrologic
cycle: agricultural engineers, energy managers, ecologists and environmentalists, public health
officers, industrialists, chemists, and inland navigation managers.
The major components of the hydrologic cycle are precipitation (rainfall, snowfall, hale,
sleet, fog, dew, drizzle, etc.), interception, depression storage, evaporation, transpiration,
infiltration, percolation, moisture storage in the unsaturated zone, and runoff (surface runoff,
interflow, and baseflow).
Evaporation of water takes place from the oceans and the land surface mainly due to solar
energy. The moisture moves in the atmosphere in the form of water vapour which precipitates
on land surface or oceans in the form of rain, snow, hail, sleet, etc. A part of this precipitation is
intercepted by vegetation or buildings. Of the amount reaching the land surface, a part infiltrates
into the soil and the remaining water runs off the land surface to join streams. These streams
finally discharge into the ocean. Some of the infiltrated water percolates deep to join
groundwater and some comes back to the streams or appears on the surface as springs.
This immense movement of water is mainly driven by solar energy: the excess of
incoming radiation over the outgoing radiation. Therefore, sun is the prime mover of the
hydrologic cycle. The energy for evaporation of water from streams, lakes, ponds and oceans and
other open water bodies comes from sun. A substantial quantity of moisture is added to the
atmosphere by transpiration of water from vegetation. Living beings also supply water vapor to
the atmosphere through perspiration. Gravity has an important role in the movement of water on
the earth’s surface and anthroprogenic activities also have an increasingly important influence on
the water movement.
An interesting feature of the hydrologic cycle is that at some point in each phase, there
usually occur: (a) transportation of water, (b) temporary storage, and (c) change of state. For
example, in the atmospheric phase, there occurs vapor flow, vapor storage in the atmosphere and
condensation or formation of precipitation created by a change from vapor to either the liquid or
solid state. Moreover, in the atmosphere, water is present in the vapor form while it is mostly
(saline) liquid in the oceans.
Global scale
From a global perspective, the hydrologic cycle can be considered to be comprised of three
major systems; the oceans, the atmosphere, and the landsphere. Precipitation, runoff and
evaporation are the principal processes that transmit water from one system to the other. This
illustration depicts a global geophysical view of the hydrologic cycle and shows the interactions
between the earth (lithosphere), the oceans (hydrosphere), and the atmosphere. The study at the
global scale is necessary to understand the global fluxes and global circulation patterns. The
results of these studies form important inputs to water resources planning for a national, regional
water resources assessment, weather forecasting, and study of climate changes. These results
may also form the boundary conditions of small-scale models/applications.
Catchment Scale
While studying the hydrologic cycle on a catchment scale, the spatial coverage can range from a
few square km to thousands of square km. The time scale could be a storm lasting for a few
hours to a study spanning many years. When the water movement of the earth system is
considered, three systems can be recognized: the land (surface) system, the subsurface system,
and the aquifer (or geologic) system. When the attention is focused on the hydrologic cycle of
the land system, the dominant processes are precipitation, evapotranspiration, infiltration, and
surface runoff. The land system itself comprises of three subsystems: vegetation subsystem,
structural subsystem and soil subsystem. These subsystems subtract water from precipitation
through interception, depression and detention storage. This water is either lost to the
atmospheric system or enters subsurface system. The exchange of water among these subsystems
takes place through the processes of infiltration, exfiltration, percolation, and capillary rise.
Fig. 1.2 shows the schematic of the hydrologic cycle at global scale, in the earth system, and
micro-scale view of the cycle in the land system. Fig. 1.3 gives a schematic presentation of the
hydrologic cycle of the earth system. Detailed schematic of the hydrologic cycle in the land
system is shown in Fig. 1.4.
Precipitation
Evaporation
Precipitation
ET
Surface
Land system runoff
Exfiltration Infiltration
Structural Soil
Vegetation system
system Depression
system Inter- Surface
ception storage
runoff
Throughflow
Infiltration
Moisture supply
Fig. 1.4 A detailed schematic of the hydrologic cycle in the land system.
The range of spatial and temporal dimensions of many processes related to the hydrologic cycle
is shown in Fig. 1.5.
1.4 Mathematical Representation of the Hydrologic Cycle
The quantities of water going through the various components of the hydrologic cycle can be
evaluated by the so-called hydrologic equation, which is a simple spatially-lumped continuity or
water budget equation:
I - Q = ΔS (1.1)
where I = inflow of water to a given area during any given time period, Q = outflow of water
from the area during the selected time period, and ΔS = change in storage of water in the given
area during the time period. If I and Q vary continuously with time, then equation (1.1) can be
written as
GLOBAL
GLOBAL GLOBAL
WEATHER C02
SYSTEMS VARIATIONS
10000 km
DEVELOPMENT
OF MAJOR
RIVER BASINS
SOIL
1000 km FORMATION
RUNOFF
CYCLE
SPACE
SHALLOW
10 km GROUND
WATER
NUTRIENT
CIRCULATION
CYCLES
THUNDER
1 km
STORM
TIME
t t
or S (t ) S (0) I (t )dt Q(t )dt
0 0
where P is precipitation, QSI is surface inflow, QGI is ground water inflow, E is evaporation from
the watershed, QSO is surface water outflow, QGO is ground water outflow, S is change in the
storage of water in the watershed, and is a discrepancy term. For large watersheds, QGI and QGO
are usually negligible. The discrepancy term is included in equation (1.4) because the sum of all
other terms may not be zero due to measurement errors and/or simplifying assumptions.
However, a small value of discrepancy term does not necessarily means that all other terms have
been correctly measured/estimated.
Depending on the specific problem, the terms of equation (1.1) may be further
subdivided. For example, when applying the hydrologic equation for short time intervals, the
change in the total water storage (S) may be subdivided into several parts: changes of moisture
storage in the soil (M), in aquifers (G), in lakes and reservoirs (L), in river channels (SC),
in glaciers (SG), and in snow cover (SS). Thus, S can be expressed as:
The hydrologic equation may be applied for any time interval; the computation of the mean
annual water balance for a basin being the simplest, since it is possible to disregard changes in
water storages in the basin (S), which are difficult to measure and compute. In general, the
shorter the time interval, the more stringent are the requirements for measurement or
computation of the components and the more subdivided are the terms of equation (1.1). This
results in a complex equation which is difficult to close with acceptable errors.
The hydrologic equation may be applied for areas of any size, but the complexity of
computation greatly depends on the extent of the area under study. The smaller is the area, the
more complicated is its water balance because it is difficult to estimate components of the
equation. Finally, the components of the hydrologic equation may be expressed in terms of the
mean depth of water (mm), or as a volume of water (m3), or in the form of flow rates (m3/s or
mm/s).
Watershed changes can be categorised as point changes and non-point changes. Structural
changes, such as dam construction, channel improvement, and detention storage are examples of
point changes and affect watershed response in terms of evaporation, seepage, residence time,
etc. Forestry, agriculture, mining, and urbanization are non-point land use changes. A qualitative
discussion of the hydrologic consequences due to watershed changes is given in the following.
Natural system
Oceans
Decreasing channel roughness increases flow velocity and peak discharge for the same
channel size. The removal of vegetation, lining of the channel, and proper maintenance can
greatly reduce roughness. The other alterations, such as straightening the channel, maintenance
of bands, or increasing slope, significantly affect travel time and flow velocity. Depending upon
the bed material, infiltration through the bed and banks also modifies flow characteristics.
1.8 Closure
Water is central to the environment and is considered to be the driver of nature. All studies
dealing with water are concerned with some part or the other of the hydrologic cycle. The
occurrence and availability of water is closely linked with human development. As a result of
many human activities, the features/characteristics of many components of the hydrologic cycle
are undergoing changes which are likely to result in long-term changes in the climate. It is
necessary to predict the magnitude and extent of these changes so that the required ameliorative
measures can be initiated well in time.
References
NRC (1982). Scientific Basis of Water Resource Management. National Research Council,
National Academy Press, Washington D.C.
NRC (1991). Opportunities in the Hydrologic Sciences. Committee on ‘Opportunities in the
Hydrologic Sciences’ of National Research Council, National Academy Press, Washington
D.C.
Seckler, David, Upali Amarasinghe, Molden David, Radhika de Silva, and Randolph Barker.
(1998). World water demand and supply, 1990 to 2025: Scenarios and issues. Research
Report 19. International Water Management Institute. Colombo, Sri Lanka.
Shiklomanov, I. A. (1999). World water resources: modern assessment and outlook for the 21st
century. (Prepared in the framework of IHP, Unesco). State Hydrology Institute, St.
Petersburg.
Singh,V.P. (1992). Elementary Hydrology. Prentice Hall, Engelwood Cliffs, New Jersey.
Unesco (1978). World water balance and water resources of the Earth. United Nations
Educational, Scientific, and Cultural Organization, Paris.
WMO (2008). Guide to Hydrological Practices, Volume I: Hydrology – From Measurement to
Hydrological Information. WMO-No. 168, World Meteorological Organization, Geneva.
Module 2
Measurement and Processing of Hydrologic Data
2.1 Introduction
2.1.1 Methods of Collection of Hydrologic Data
2.2 Classification of Hydrologic Data
2.2.1 Time-Oriented Data
2.2.2 Space-Oriented Data
2.2.3 Relation-Oriented Data
2.3 Design of Hydrometeorological Data Networks
2.3.1 Classification of Observation Networks
2.4 Precipitation Networks
2.5 Stream Gauging Networks
2.5.1 Network Design Process
2.5.2 Criteria for Location of Stations
2.5.3 Evaluation and Adequacy of Networks
2.5.4 Site Selection Surveys
2.5.5 General site selection guidelines
2.5.6 Criteria for Water Level Gauging Sites
2.5.7 Criteria for Streamflow Measurement Sites
2.5.8 Criteria for Natural Control Sites
2.5.9 Criteria for Artificial Control Sites
2.5.10 Bureau of Indian Standards (BIS) criteria for selection of river gauging sites
2.5.11 World Meteorological Organisation (WMO) criteria for selection of site
2.5.12 International Standards for Hydrometry
2.6 Errors in Hydrological Observations
2.6.1 Sources of Errors
2.6.2 Secondary Errors of Measurement
2.7 Validation of Hydrologic Data
2.7.1 Levels of Validation
2.7.2 Primary Validation
2.7.3 Secondary Validation
2.7.4 Hydrological Validation
2.7.5 Validation of Climatic Data
2.7.6 Single series tests of homogeneity
2.7.7 Multiple stations validation
2.A Definitions of terms related to measurement errors
2.11 References
2.1 Introduction
Data are the foundations on which any analysis rests. The practice of hydrological measurements
is very old. Kautilya initiated systematic precipitation measurements in India in the fourth
century BC. Streamflow was probably first monitored by Hero of Alexandria in the first century
AD. Equipment and techniques of hydrologic data collection have evolved with growth in
technology and water sciences.
For a water resources study, one needs data of a number of variables in the vertical as
well as horizontal planes. The data needed for water resources development come from a vast
swath of disciplines: hydro-meteorologic, geomorphologic, agricultural, pedologic, geologic,
hydrologic, social, economic, ecological and environmental sciences, etc. Hydrometeorologic
data include rainfall, snowfall, temperature, radiation, humidity, vapor pressure, sunshine hours,
wind velocity, and pan evaporation. Agricultural data include crop cover, irrigation application,
and fertilizer application. Pedologic data include soil type and texture; soil particle size; porosity;
moisture content; steady-state infiltration, and saturated hydraulic conductivity. Geologic data
include stratigraphy, lithology, and structural controls. Frequently, data on the type, depth and
areal extent of aquifers are needed. Ecological and environmental data includes water quality
variables, aquatic plants and animals and their habitats. Each data set is examined with respect to
homogeneity, completeness, and accuracy.
Data processing also includes aggregation of data observed at a given time interval (e.g.
hourly) to a different interval (e.g., daily and daily to monthly) or disaggregation, i.e., conversion
from a long to short (say daily to hourly) time step is also carried out.
Typical stages in hydrological data processing are:
Scrutiny of raw data;
data entry to computer, validation, and correction; and
data archival and dissemination.
Direct Measurement
This is the most common way to measure hydrometeorological variables, such as precipitation
and streamflow. A gauging site is established and is equipped with the devices that can measure
the variable(s) of interest. In case of manual observations, an observer visits the site, measures
the values of the concerned variables, and records or transmits them to the controlling office for
processing and storage. On the other hand, at an automated hydrologic or weather station the
seasons can measure a number of hydrometeorological variables and store/transmit the data to
the controlling office without any human intervention. The equipment may be programmed to
transmit the data at selected time interval or it can be interrogated as per the needs to get the data.
With improvement in communication technology, it is possible to get the desired data from the
stations widely spread over an area at a central place in real-time.
Remote Sensing
In this technique, the data about an object are obtained without coming in physical contact with
the object. This technique is now very commonly used to provide spatial data of terrain features.
Similarly, weather radars are being increasingly used for measurement of precipitation.
OBSERVATION INPUT
STORE
FEEDBACK RETRIEVE
DATA
PROCESSING
OPERATIONS
UTILISE
Observation and input, processing and storage, Retrieval and use, Feedback,
Hydrologic data can also be classified as time varying or time non-varying data. The time
non-varying or static data includes most space-oriented data which do not change (or very-very
slowly change) with time, for example catchment topographic map, soil map, etc. Some features,
such as river network and land use in a catchment, might gradually change with time and can be
considered as semi-static. A brief description of each type of data is presented next.
Depending on the frequency of observations, the time-series data can also be classified as:
Equidistant time-series data which are the measurements made at regular intervals
(hourly, daily); the reported values may be instantaneous, accumulated or averaged.
Cyclic time-series data are the data measured at irregular intervals of time but the time
sequence is repeated regularly. For example, the observation of at many places river stage
is measured daily at 08:30 and 17:30 hrs.
Values of non-equidistant data series are observed when some specified event takes
place. For example, in a tipping bucket rain gauge, the bucket tips after a certain depth of
rain has fallen and the value is recorded.
Since the requirement of water resources data depends on their end use, it is difficult to
formulate general rules on network design. While designing hydrologic networks, the decisions
to be taken are:
Areal networks are established to get spatial characteristic of the variables over an area
while the linear networks are created for rivers, canals, etc. On the basis of purpose, the
networks can be classified in three categories: basic (to get the fundamental characteristics of the
variables of interest), specific (to gather data for some specific purpose, e.g., a reservoir project),
and temporary (which are in operation for a short period of time).
The optimum density of a precipitation gauge network depends on the purpose for which the data
are to be used. For example, measurements of precipitation for flood forecasting require denser
networks as compared to that for rainfall-runoff modeling. WMO (2008) has recommended the
following (Table 2.1) as minimum network densities for precipitation stations.
Table 2.1 Recommended minimum densities of stations (area in km2 per station) [Source: WMO
(2008)].
Physiographic unit Precipitation Evaporation
Non-recording Recording
Coastal 900 9000 50000
Mountains 250 2500 50000
Interior plains 575 5750 5000
Hilly/undulating 575 5750 50000
Small islands 25 250 50000
Urban areas - 10–20 -
Polar/arid 10000 100000 100000
The optimum network should make it possible to derive required information with desired
accuracy. The optimum number of rain gauge stations (N) in a network is given by (BIS 1968):
N = [Cv /p]2 (2.1)
where Cv = the coefficient of variation of the precipitation values of the existing rain gauge
stations, and p = the allowable maximum percentage error in the estimate of basin mean rainfall.
A typical value of p is 10 percent. Here, Cv is computed by
Cv = 100*s/Pm (2.2)
In which s is the standard deviation and Pm is the mean rainfall of the existing stations.
Obviously, a decrease in the percentage error would mean an increase in the number of gauges
required. Mukherjee and Kaur (1987) have proposed a modified form of eq. (2.1) by including
the mean correlation (r) of precipitation over the area
WMO recommends that the precipitation (amount and form) should be measured with an
accuracy of 3–7% and rainfall intensity with 1 mm/hr at the 95 per cent confidence interval.
Snow depth below 20 cm should be accuracy of less than 1 cm and depth above 20 cm should
not have more than 10% error. The recommended accuracy for evaporation range 2–5% and for
wind speed 0.5 m/sec.
Example 2.1: A catchment has 6 rain gauges and the annual rainfall at these has been measured
as 750, 540, 465, 493, 421, and 780 mm. Find out the optimum number of rain gauges for the
basin if the error of estimation is required to be kept below 10%.
Solution: For the data given, mean = 574.83 mm and standard deviation = 152.59 mm. Thus
Cv = 100*152.59/574.83 = 26.54
Hence, using eq. (2.1), the optimum number of rain gauges for the basin (N) is
N = (26.54/10)2 = 7.04
This means that 7 rain gauge stations are required in the basin and the existing network of 6 rain
gauges is slightly inadequate. It needs to be strengthened by adding one new gauge so that the
estimate of rainfall depth has stipulated accuracy.
Every major stream should be gauged near its mouth and its major tributaries should also
be gauged as feasible. Naturally, gauging depends on the existing and likely development in the
basin. According to WMO, the first gauging station is selected at the most upstream location
where the drainage area is about 1300 km2. The second station is located at a point in the
downstream direction where the drainage area is approximately doubled. The WMO
recommendations for a minimum density of stream gauging stations are given in Table 2.2.
Table 2.2 Recommended minimum densities of stations (area in km2 per station) [Source: WMO
(2008)].
Physiographic unit Streamflow Sediments Water quality
Coastal 2750 18300 55000
Mountains 1000 6700 20000
Interior plains 1875 12500 37500
Hilly/undulating 1875 12500 47500
Small islands 300 2000 6000
Urban areas - - -
Polar/arid 20000 200000 200000
Stations are also established in the basin at the sites where significant changes in the
volume of flow are noticed, for example downstream of the confluence of a major tributary or at
the outflow point of a lake etc. In case a suitable location is not available downstream of the
confluence, the sites can be located upstream of the confluence, preferably on the tributary.
While establishing sites at the downstream of confluence, it should be ensured that no other
small stream joins the main river before the station so that correct assessment of the contribution
of the tributary to the main river is obtained.
The distance between two consequent stations on the same river may vary from about 50
km to several hundred kilometers, depending on many factors. The drainage area computed from
origin up to two consecutive observation sites on a large river should preferably differ by more
than ten percent so that the difference in quantities of flow at the two stations is significant.
Sometimes stations are set up due to hydrological significance, say, to determine the flow
contribution from a typical catchment.
A different approach is recommended for small independent rivers which flow directly
into the sea (for example, the rivers in Western Ghats). In such cases, the first hydrological
observation station is to be established on a stream that is typical of the region and further
stations are added to the network to cover the area and obtain information about the variability.
Stream in the area whose flows are low should not be avoided from the network. Absence of
stations from low flow streams may lead to wrong assessment of the surface water potential of
the area if it has been evaluated just on the basis of the data from the high flow streams. Thus,
great care is to be exercised to ensure that all distinct hydrological features are adequately
covered by the gauging network.
An existing gauging network needs periodic review. The developments that take place in
the basin like construction of new water resources development projects may warrant addition or
closure of the sites. Often the rivers are polluted by the discharge of affluents from industries. A
need may also arise to establish stations to monitor the quality of water in the river.
Regarding the accuracy desired in measuring river water depth and discharge, WMO
recommends that the water depth measurement should have accuracy of about 2%, velocity of
flow 2–5%, and discharge about 5%. Suspended sediment concentration should be estimated
with accuracy of 10%.
These steps and related topics are further elaborated in the following sections.
Based on the answers, stations which can be added, deleted or relocated are identified.
Financial Aspects
In addition to technical financial considerations are also important in network design because the
stations cannot be established without adequate money, equipment cannot be purchased and
operated, and staff cannot be hired. Hence, after the preliminary design of the network has been
completed, the expenditure to establish stations and the cost of operating them should be
estimated. These monetary requirements should match with the budget so that the proposed
network is sustainable. In case of deficit in the budget, the network should be re-aligned or
additional budget should be arranged. Stations in the network may be prioritized to best attain the
objectives, given the constraints. Table 2.3 gives a possible categorization of stations based on
their relative importance.
Prior to removing/shifting an existing station in the network, the main users and beneficiaries of
the data should be consulted.
By now, the target location for the gauging station will have already been identified on a
1:250,000 or similar map in earlier steps. However, 1:250,000 is too small a scale for final site
selection purposes. Large-scale topographic maps (1:50,000) should be checked to identify
possible sites within the target zone. Reconnaissance surveys should be undertaken by an
experienced hydrologist along with a person familiar with the area. As the hydraulic conditions
and river characteristics vary considerably from non-monsoon to monsoon season,
reconnaissance survey in both the seasons would facilitate correct decision on the suitability of
the site. When the establishment of site cannot wait that long, the suitability of the site for
hydrological observation could be decided after single inspection assisted by toposheet studies
and other relevant field investigations including measurements of width and depth.
At sites of interest, ownership of the land and approach should be ascertained. The site
shall be accessible in all seasons and all weather. It is important to use updated maps since most
surveys were completed several decades ago and things may have undergone large changes.
Recent situation can be obtained from remote sensing images or internet sites such as Google
Earth.
Information on the historical high flood level should also be collected (by local enquiry
and / or by examining the available landmarks) during the inspection. An all weather accessible
site located in a straight uniform reach free from weeds, rock outcrop, pools and back water
effect with stable non-overflowing banks with flow confined to single channel normal to the
selected cross-section of measurement would be an ideal site for stream gauging.
On completion of the reconnaissance surveys, one or more locations are shortlisted for
further consideration. After this, field surveys are carried out and the cross-section of the
proposed site is surveyed. If artificial controls (e.g. a weir) are planned, it will be necessary to
survey the river for some distance upstream and downstream to ascertain the impact on flows
and water levels. It will also be important to understand what type of control exists and to make
sure that the location will not be impacted by variable backwater effect of any structure.
1) Generally the use of artificial controls should be limited to small but important rivers (<
100 m wide) and for special investigations in artificial channels.
2) Existing structures may be adapted for the purpose of flow measurement, wherever
feasible.
3) The sensitivity of upstream area to increased levels should be assessed. For example, will
the installation of the structure cause a potential, increased risk of flooding.
4) A minimum length of straight approach channel of five times the maximum width of the
water surface is recommended for most structures, except for thin plate weirs where ten
times the maximum channel width is recommended. However, research has shown that
for triangular profile weirs accurate results can be obtained even if the weir is only twice
the channel width from an upstream bend.
5) Thin plate weirs are particularly sensitive to upstream velocity distribution.
6) Like all controls, it is essential that the structure creates a sensitive stage-discharge
relationship. In wider rivers, this can be a problem at low flows. Structures such as the
triangular profile flat “v” weirs provide such sensitivity.
7) The discharge coefficients of many structures vary when the velocity head in the
upstream approach channel becomes large in relation to the depth of flow. A
dimensionless number which describes this is the Froude number (Fr). To prevent water
surface instability in the approach channel the Froude number should generally not
exceed 0.5.
8) The design of the structure should be such as to minimise upstream sediment deposition
or downstream scouring. In rivers with high bed loads the use of structures which
significantly reduce the stream velocity is not recommended.
9) On rivers which are navigable or those which are important fish migration routes the use
of flow measurement structures should be avoided and some other form of flow
measurement considered.
2.5.10 Bureau of Indian Standards (BIS) criteria for selection of river gauging sites
The ideal requirements for a good gauging site as enunciated in the standard IS 1192-1981
"Velocity - Area methods for measurement of flow of water in open channels" are given below.
The accuracy of measurement of discharge by velocity area method is increased if the site is
selected considering these aspects.
The site selected should comply, as far as possible, with the following essential requirements:
a) The reach of the open channel at the gauging site shall be straight and of uniform cross
section and slope, as far as possible, in order to avoid irregularities in velocity
distribution. The length of the reach need not be more than 1600 m and should not be less
than 400m. When the length of the straight channel is restricted, it is recommended for
current meter measurements and the straight length upstream of the measuring cross
section should be twice that on the downstream.
(Note: In case of artificial channel, the minimum length of straight reach should
preferably be such as to give a drop in water level of 0.06 m. or the minimum length
should be equal to four times the width of the channel, whichever is larger.)
b) The depth of water in the selected reach shall be sufficient to provide for the effective
immersion of either the current meters or floats, whichever are to be used.
c) When near a confluence, the site, if located on a tributary shall be sufficiently upstream
preferably beyond the backwater effect; and if located on the main stream, upstream or
downstream of the confluence it shall be beyond the disturbances due to the tributary.
d) The site should be easily accessible at all times of the year.
In addition to the above requirements, the following points shall be taken into
consideration as desirable requirements in the selection of the gauging site.
a. The flow should be confined in a single channel and there should be no overflow as far as
possible. Where this is not possible, the site in which minimum number of channels exist
and the flood plain has minimum width should be preferred.
b. Where these requirements cannot be met (for instance- when in alluvial rivers the river
bed is changing during the period of measurement, or when, under flood conditions, the
river is not confined to a single channel in embankments), a gauging site shall be chosen
such that the bed change and/or overflow is minimum. Floodplain, if cannot be avoided,
shall be of minimum width, as smooth as possible, with a distinct channel, and clear of
bushes and trees. The flow in the over bank or floodplain section (s) shall be measured
separately and added, treating the whole as a composite section.
c. The site shall be remote from any bend or natural or artificial obstruction if disturbances
of the flow are likely to be caused thereby.
d. The orientation of the reach should be such that the direction of flow is as close as
possible normal to that of the prevailing wind.
e. Sites at which there is a tendency for vortex formation should be avoided.
f. The site should, as far as possible, be free from trees and obstructions which may
interfere with flow and clear vision during observation.
g. The site shall be free from aquatic growth which is likely to interfere with the
measurement of depth and the current meter reading.
h. The site shall be away from the back water zone caused by any structure on the river.
i. The site should be sufficiently away from the disturbance caused by rapids and falls. etc.
It needs to be stressed that no statistical analysis can replace correct observations because
spurious and systematic errors cannot be eliminated by such analysis. Only random errors can be
characterized by statistical means.
Errors
E in hydrologic meeasurements can be claassified in tthree categoories: system matic,
random, and spurious. Fig. 2.3 giives a graph
hical depictioon of errors. These are ddiscussed in ddetail
in the folllowing.
Regarding
R prrecipitation, WMO (19882) listed thee following errors for wwhich adjusttment
needs to be made to get a near accurate
a estiimate of preecipitation frrom a measuured precipittation
report.
The first six errors listed above are systematic and are listed in order of general importance. The
net error due to blowing and drifting snow and due to splash in and out of water can be either
negative or positive while net systematic errors due to the wind field and other factors are
negative. Since for liquid precipitation the errors listed at (e) and (f) above are near zero, the
general model for adjusting the data from most gauges takes the form
The data needed to make the adjustments include wind speed, drop size, precipitation intensity,
air temperature, humidity and other characteristic of the gauge site.
Random Errors
Random errors vary in an unpredictable manner, both in magnitude and sign, when repeated
measurements of the same variable are made under the same conditions. Random errors cannot
be eliminated, but their impacts can be reduced by repeated measurements of the variable. These
are equally distributed about the mean or ‘true’ value. The errors of individual readings may be
large or small, e.g., the errors in a staff gauge reading where the water surface is subject to wave
action. Usually, they compensate with time or are minimized by taking a sufficient number of
measurements. The uncertainty of the arithmetic mean computed from n independent
measurements is several times smaller than the uncertainty of a single measurement. The
distribution of random errors can usually be assumed to be normal (Gaussian). For certain cases,
normal distribution can or should be replaced by other statistical distributions. These errors can
be identified by a statistical-outlier test that gives a rejection criterion.
In measuring rainfall, random errors could arise due to spilling of the water when
transferring it to the measuring jar, leakage into or out of the receiver, observational error etc.
The others random errors which could be due to observer include:
i. misreading and transposing digits,
ii. misrecording because of faulty memory,
iii. recording the data at the wrong place on the recording sheet,
iv. making readings at improper interval,
v. incorrect dating of the report,
vi. incorrectly reading or communicating the data to a reporting centre, etc.
It appears, therefore, that computerization and automation may be solution to reduce the
error. However, even without human intervention chances of erroneous reading in case of
precipitation may be possible because of
i. evaporation from gauge,
ii. overflowing gauge,
iii. mechanical or electrical mal-functions.
Spurious errors
These arise due to human mistakes or instrument malfunction or some abnormal external cause.
Reported data appear to be clearly in error. Sometimes the errors become obvious, for example,
wrong placement of decimal and the data can be easily corrected in such cases. In other cases,
the concerned measurements may have to be discarded. For example, an animal may drink water
from the evaporation pan and introduce errors in the data. Sometimes, such errors may be readily
detected but it may not be easy to correct them.
2 2 2
Q Q Q
(eq )
2
e x e y ez
x y z
(2.5)
where ∂Q/∂x, ∂Q/∂y and ∂Q/∂z are the partial differentials of the function expressing explicitly
the relationship of the dependent variable with the independent variables.
Data validation is the means by which data are checked to ensure that the corrected
values are the best possible representation of the true values of the variable. Data validation
procedure includes primary and secondary data validation. Validation of hydrologic data must
never be considered as a purely statistical or mathematical exercise. Staff involved in it must
have a background in hydrology and must understand the field practices. To understand the
source of errors, one must understand the method of measurement or observation in the field and
the typical errors of given instruments and techniques. Knowledge of the method of
measurement or observation influences our view of why the data are suspect.
The input variables in an analysis may be directly measured (e.g., rainfall) or they may be
derived using a relationship with one or more variables (e.g., discharge that has been obtained
from a rating curve). In the latter case the error in the variable (discharge) depends both on field
measurements and the error in the relationship. An error may also be introduced if the
relationship is no longer valid or the values are extrapolated outside the applicable range.
Validation involves different types of comparisons of data and includes the following:
Multiple stations/data
between two or more measurements at nearby stations, e.g. flow at two sites along a
river; and
between measurements of different but related variables, e.g., rainfall and river flow.
Note that none of the procedures of data validation are absolutely objective and there is
no guarantee that all the errors will be captured and removed. They are basically tools to screen
out suspect data which are to be further examined by other tests and corroborative facts. When it
is ascertained that a particular value is incorrect, an alternative value that is likely to be closer to
the true value is substituted. Since each hydrological variable has distinct characteristics, it is
necessary that specific validation techniques be designed for each variable. Further validation a
pure statistical exercise; the properties and behavior of the variable under consideration should
always be kept in mind.
When hydrological variables have a high auto-correlation, such as ground water levels, or
the data has high correlation with neighboring stations, the validation and data correction can be
carried out with a higher level of confidence. However, processes such as convective rainfall
show a great temporal and spatial variability. It is difficult to ascertain the behavior of such
processes with the desired degree of confidence. Based on available information and statistical
properties, if it is not possible to conclude whether the suspected value is erroneous or not, such
value is not changed but is flagged as doubtful. All doubtful data are further validated on the
basis of additional information.
Ideally all the hydrological data should be subjected to hydrological validation. For
historical data to which no (or few) checks have been applied, hydrological validation provides
an effective check on the quality and reliability of records. Thorough hydrological validation
requires a high level of professional expertise and can be time consuming. Required man power
and time may not be available always. Therefore, this validation may be applied selectively.
Finally, the validation may be able to identify a particular section of record/ data item that is
unreliable, but it may not always be possible to correct the values.
Trend generally does not become evident for a number of years and so the tests must be
carried out on long data series, often aggregated into annual series. Trend may result from a wide
variety of factors including:
Change of instrumentation
Change of observation practice or observer
Local shift in the site of the station
Growth of vegetation or nearby new buildings affecting exposure of the station
Effects of new irrigation in the vicinity of the station (affecting humidity, temperature
and pan evaporation)
Effects of the urban heat island with growing urbanisation
Global climatic change
The presence of trend does not necessarily mean that part of the data are erroneous but
that the environmental conditions have changed. Unless there is reason to believe that the trend is
due to instrumentation or observation practices or observer, the data should not generally be
altered but the existence of trend noted in the station record.
While examining the current data, the plot should include the past time series of
sufficient length to ensure that there are no discontinuities between one batch of data received
from the station and the next. This will ensure that the data are being entered against the correct
station and correct date/time.
For climatic variables, which have strong spatial correlation, such as temperature, the
graphs will generally run along closely, with nearly the same variation with the mean separation
representing some location factor such as altitude. Abrupt or progressive straying from this
pattern will be evident from the comparative plot, which would not necessarily have been
perceived at primary validation from the inspection of the single station. An example might be
the use of a faulty thermometer, in which there might be an abrupt change in the plot in relation
to other stations. An evaporation pan affected by leakage may show a progressive shift as the
leak develops. This would permit the data processor to delimit the period over which suspect
values should be corrected.
Comparison of series may also help in accepting of values which might have been
suspected in primary validation because they fell outside the warning range. Where two or more
stations display the same behavior, there is strong evidence to suggest that the values are correct.
An example might be the occurrence of an anomalous atmospheric pressure in the vicinity of a
tropical cyclone.
Comparison plots provide a simple means of identifying anomalies but not of correcting
them. This may be done through correlation or regression analysis, spatial homogeneity testing
(nearest neighbor analysis) or double mass analysis.
Residual series
An alternative method of displaying comparative time series is to plot the differences. This
procedure is often applied to river flows along a channel to detect anomalies in the water balance
but it may equally be applied to climate variable to detect anomalies and to flag suspect values or
sequences.
In a plot of the relationship, the suspect values will generally show up as outliers but such
plots provide no indication of the time sequencing of the suspect values. One will not be able to
say whether the outliers were scattered or contained in one block. For seasonal data, the
relationship should be derived for a period within the same season as the suspect values. (The
relationship may change between seasons). The identified suspected values should be removed
before deriving the relationship, which may be applied to compute corrected values to replace
the suspect ones.
Closure
The science of hydrology deals with immense volumes of data of a number of variables. Since
data collection is expensive process, it is necessary that the data collection campaigns are
carefully planned and executed. Errors may creep in the data due to various causes. Therefore,
before the data can be put to use, it needs to be screened and validated to remove these errors.
After the data has undergone quality control checks, it is ready for use in planning, design, and
operation.
Annexure A
2.A Definitions of terms related to measurement errors
Definitions of some important terms related to accuracy have been sourced from WMO (2008)
and are given below.
Accuracy: The word accuracy is generally used to indicate the closeness or the agreement
between an experimentally determined value of a quantity and its true value. It is the extent to
which a measurement agrees with the true value. An accurate result closely agrees with the
actual value for that quantity. In other words, accuracy tells us how close a measurement is to an
accepted standard. Precision describes how well repeated measurements agree with each other. It
tells us how close two or more measurements agree. It is worth mentioning here that precision
does not necessarily indicate anything about the accuracy of the measurements. An experiment is
considered good when it is both precise and accurate.
An experiment is said to have high precision, if it has small random error. It is said to
have high accuracy, if it has small systematic error. There may be four possibilities for
characterizing the obtained experimental data, as shown in Fig. 2.A.1: (B) precise and accurate,
(A) precise and inaccurate, (D) imprecise and accurate, and (C) imprecise and inaccurate. In
hydrological observation, the objective is to reduce both systematic and random errors as much
as possible. However, for economy of effort, one must try to strike a balance between these two
sources of error, giving greater weight to the larger of the two.
B
A
Precision
C D
Accuracy
Fig. 2.A1 Measurement of rainfall by four rain gauges. Gauge A is precise, inaccurate; gauge B
is precise, accurate; gauge C is imprecise, inaccurate; gauge D is imprecise, accurate. The
innermost circle indicates the true value.
Confidence interval: The interval which includes the true value with a prescribed probability and
is estimated by using the sample data.
Confidence level: The probability that the confidence interval includes the true value.
Error: The difference between the result of a measurement and the true value of the quantity
measured. This term is also used for the difference between the result of a measurement and
the best approximation of the true value, rather than the true value itself. The best
approximation may be a mean of many measurements.
Expected value: The best approximation of the true value, which may be a mean of many
measurements.
Measurement: An action intended to assign a number as the value of a physical quantity in stated
units. The result of a measurement is complete if it includes an estimate of the probable
magnitude of the uncertainty.
Precision: The closeness of agreement between independent measurements of a single quantity
obtained by applying a stated measurement procedure several times under prescribed
conditions.
Accuracy has to do with closeness to the truth, precision has to do only with closeness together.
Precision of observation or of reading is the smallest unit of division on a scale of measurement
to which a reading is possible either directly or by estimation.
Random error: That part of the error that varies in an unpredictable manner, in magnitude and in
sign, when measurements of the same variable are made under the same conditions.
Range: The interval between the minimum and maximum values of the quantity to be measured,
for which the instrument has been constructed, adjusted or set. It can be expressed as a ratio of
maximum and minimum measurable values.
2.11 REFERENCES
BIS (1968). Recommendations for establishing network of raingauge stations. Bureau of Indian
Standards, New Delhi.
Herschey, R. W. (Editor). (1978). Hydrometry Principles and Practices. John Wiley and Sons,
New York.
Herschey, R. W. (Editor). (1995). Streamflow Measurement. Chapman & Hall, E & F N SPON,
London.
ISO (1982). Liquid flow measurement in open channels – part 2: Determination of stage-
discharge relation, ISO 1100/2-1982(E). International Standards Organization, Switzerland.
ISO (1983). Measurement of Liquid Flow in Open Channels. ISO Standards Handbook 16.
International Organization for Standardization, Switzerland.
ISO (1996). Hydrometry — Measurement of liquid flow in open channels — Part 1: Guidelines
for selection, establishment and operation of a gauging station. ISO 1100-1:1996,
International Organization for Standardization, Geneva.
WMO(1988). Manual on water quality monitoring -- planning and implementation of sampling
and field testing. Operational Hydrology Report No. 27, WMO No. 680. World
Meteorological Organization, Geneva.
WMO(1989). Management of ground water observation programmes. Operational Hydrology
Report No. 31, WMO No. 705. World Meteorological Organization, Geneva.
WMO (2008). Guide to Hydrological Practices. Volume I: Hydrology – From Measurement to
Hydrological Information. WMO No. 168. World Meteorological Organization, Geneva.
Module 3
Measurement and Processing of Precipitation Data
Keywords: Measurement, Precipitation, Rain Gauges, Tipping Bucket, Weather Radar, Snow,
Processing, Spatial Interpolation, Thiessen Polygon, Isohyetal, Filling-in, Normal Ratio,
Disaggregation
Since the size, shape and exposure affect the precipitation caught by a gauge, standard
gauges are used so that the observations can be compared and there is consistent use of data. The
Symon's rain gauge has been adopted as the standard rain gauge In India. To install the gauge, a
masonry or concrete platform is sunk into the ground and the gauge is placed on it such that the
rim of the gauge is perfectly levelled and is approximately 30 cm above the ground. In a non-
recording instrument, the precipitation measured at 8:30 am on any particular date is entered
against that date and it is understood that this precipitation has been received in 24 hours
preceding the time and day of observation.
Standard or ordinary rain gauges (ORG) are manually read, commonly once a day.
Self-recording rain gauges (SRRG) record the rainfall depth in the form of a continuous plot.
Automatic rain gauges with data logger. Data are stored in digital form either as rainfall at
fixed interval or as timings for each event of rainfall of fixed depth.
An ordinary rain gauge or ORG (see Fig. 3.1) is a simple device consisting of a collector
and a funnel. A standard rain gauge consists of a circular collector funnel with a brass or gun
metal rim and the rainfall collection area of either 200 cm2 (diameter 159.5 mm) or 100 cm2
(diameter 112.8 mm). The funnel leads to a base unit, partly embedded in the ground and
containing, a polythene or glass collector bottle. The gauge is read once (usually) or twice daily
and any rain gathered in the collector is poured into a measuring graduated measuring glass
cylinder to determine rainfall depth in millimetres.
Since rainfall received by the ORG is measured and recorded manually, errors may be
introduced by the observer in many ways. Observer may read the measuring cylinder incorrectly
or may incorrectly note
n the amo ount in the data
d sheet. O
Observer maay also read the gauge aat the
wrong tim me or may noten the am
mount againstt the wrong date. Obserrver may alsso use the w wrong
2 2 2 2
measurinng glass (e.g g., 200 cm glass
g for a 100
1 cm gauuge or 100 ccm glass forr 200 cm gauge.
Hence prrecipitation will
w either bee half or twiice the true vvalue). Besiides human ccause, errorss may
also arisee due to insttrument faultt: if the gaugge rim is dam
amaged, the collection arrea is changged; if
the funneel is blocked, water willl not reach the collectiion bottle – it may oveerflow the riim or
(partly) lost
l by evap poration. Colllector bottlee may be brroken and gaauge might be leaking. Note
that errors due to most
m of these causes wiill be difficuult to detect, particularrly if only a few
records are
a present, unless a veery large errror has beenn committedd in reading or reportinng the
values. Errors
E in the observationss at a stationn can be easi ly detected iif there is a cconcurrent reecord
from a reecording rainn gauge at th he same or a nearby statioon.
Three
T types of
o recordingg rain gauges are in genneral use: weighing typee, float typee, and
tipping bucket
b type. In weighing type instru uments, all the precipittation fallingg is continuously
recorded and it can measure
m all ty
ypes of preccipitation (liqquid or solidd).
3.2.2 Autographic
A c rain gaugee
Natural syphon
s rain gauge is frrequently useed to measuure short-perriod rainfalll. It consistss of a
circular collector
c funnnel and rain nfall recordiing mechaniism. The coollector has a gun metall rim,
2
200 mm diameter (314 cm area)) and is fixed d at 750 mm m above grouund level. Thhe funnel leaads to
a chambeer in which a float is lo ocated. Thiss float rises when rain w water enterss the chamber. A
syphon chamber
c is atttached to th
he float cham
mber. After 110 mm of raiin has fallen, the float risses to
a certainn level and thet syphon action is in nitiated. A ppen is attachhed to the ffloat and records
rainfall on
o a chart placed on a drrum in the formf of rise of the float. After waterr is siphonedd out,
the float returns to th
he original po osition.
The
T drum mo oves with the help of a mechanical
m cclock such tthat it complletes one rottation
in 24 houurs. The horiizontal axis of the chart is marked w with hours annd the vertical axis repreesents
depth of rainfall. Thee chart is chaanged every day at the sset time and the observerr reads the hhourly
rainfalls from the ch hart and reco ords in data sheets. If noo rain has faallen on a giiven day, thee pen
will tracees a horizonttal line on th
he chart reprresenting zerro rainfall. DDuring the peeriods of raiinfall,
the pen produces
p a slloping line; higher is thee intensity of rainfall, steeeper will be the slope oof the
curve pro oduced by th he pen. Rainn water drain ned by the siiphon actionn may be colllected in a vvessel
and meassured when the t chart is changed
c and
d this will be a check for the total rainnfall.
Fig. 3.2 show ws a self recoording rain gauge assembbly.
Fig. 3.2 A self recorrding rain gaauge with chart mountedd on a drum ((top) and raiinfall recordeed by
the instru
ument on a chart.
c
Measurement
M t errors can arise due too instrumentaal or observver faults. Thhe funnel maay be
(partially
y) blocked soo that waterr enters in th
he float cham mber at a raate different from the raainfall
rate, lead
dings to erroors. If the flloat is wrongly adjustedd, syphon acction is initiiated at a raainfall
depth diffferent than 10 mm. In case of very intense
i rainffall, the rise of float andd siphon mayy take
place so frequently thatt individu
ual pen tracces cannot bbe distinguisshed. Furtheer, the clockk may
malfuncttion – it may y stop, or maay be either slow or fastt and the tim mings will bee incorrect. IIf the
float getss sticky with h the walls of the cham mber, rainfalll will not bbe recorded correctly. OOn his
part, the observer maay incorrectly y read the daata from the chart.
Commonly
C an n ORG is allso installed at the site oof recording station and the daily raainfall
recorded at the two stations maay be differrent due to different exxposure connditions, diffferent
levels off the rims, ettc. Traditionaally, the dailly ORG dataa is considerred to be moore correct wwhere
there is a discrepancy y between thhe two.
3.2.3 Tipping
T Bucket Rain Ga auge
A tipping g bucket rainn gauge con nsists of a circular collecctor funnel w
with a brass or gunmetaal rim
that direccts the rain into
i a pair of tipping buckets whichh sit on a kniife edge (Figg. 3.3). Afterr rain
water hass filled one side of the bucket
b with a small amoount (say, 0.0025cm) of raain, the asseembly
tips. An electrical pu ulse is generrated on eacch tilt and iss recorded tto provide ddata of rain ffallen
with timee. At the saame time, raain water beegins to fill tthe second sside of buckket, and so oon. A
data loggger records thet occurren nce of each tip along w with the timee of tip and the data maay be
download ded as per convenience
c and need. Since
S manuall input is noot used in daata observation or
recordingg, there are little
l chances of errors. The
T instrum ment can be eeasily calibraated and cheecked
in field by
b pouring a known amount of waterr in the colleector funnel.
Measurement
M t errors may occur when
n the funnel iis (partially blocked) so that water eenters
in the tip
pping bucketts at a rates which are different
d thaan the rate oof rainfall. IIf the bucketts are
damaged d or out of balance, theey may tip after non-sstandard rainnfall or the tipping maay be
incorrecttly recorded. A faculty reed switch may
m fail to rregister tips or may doubble register tthem.
The gaug d to some fault or lighhtning strike etc. Differennces may arise in
ge electroniccs may fail due
the rainfaall recorded at the daily rain
r gauge (ORG) due too of differennt exposure cconditions.
3.2.4 Precipitation Measurement by Weather Radar
Although rain gauges are in use for a very long time, they suffer from a limitation: a rain gauge
gives the point value of rainfall but in almost all hydrologic studies, spatially averaged values are
required. Weather radars have been developed to overcome this limitation.
The word radar stands for radio detection and ranging. Weather radar sends
electromagnetic waves in all directions. When these waves collide with raindrops, an echo of the
waves is produced which is caught back by the radar. The strength of echo or returning signal
depends on the intensity of the falling rain and the time taken by the echo to reach the radar
depends upon the distance of the rainfall field from the radar.
A radar is calibrated using the data of rain gauges and relationship between the radar
reflectivity Z (mm6/m3 or the sixth power of the diameter of raindrops in mm per cubic meter of
the atmosphere) and rainfall intensity R (mm/h), known as Z-R relationship is established. A
commonly used empirical equation is:
The main advantage of the weather radar is that it can give the estimates of rainfall over
wide areas along with location and movement of storms. Weather radars can also measure
localized storms which may be missed by rain gauges. The range of radar depends on the
hardware; radars with range up to hundreds of km are common these days. The factors that
control the measurement of precipitation are its type, size of raindrops, width of the radar beam,
refraction of beam and atmospheric attenuation. Combined use of rain gauges and radars
provides much more improved rainfall estimates than what any single of them could provide.
A radar can measure precipitation over areas that are difficult to access and gives areal
coverage of rainfall distribution, it is possible to detect the movement of rainfall and most
systems are capable of automatic data processing. The spatial distribution of rainfall is
simultaneously displayed on a monitor. Due to these reasons, the use of weather radars in
precipitation measurement and warning is increasing. These have been found to be particularly
useful in flood forecasting and warning. On the negative side, the measurements from weather
radar are affected by echoes from the ground or ground clutter. In some cases, these radars have
been found to underestimate light rainfall and overestimate heavy rainfall.
Satellites are also being used to estimate precipitation over large areas and in near real-
time. Images from geostationary and polar orbiting satellites along with cloud top temperature,
shape, texture, and cloud history are used for estimation of precipitation. Another possibility is to
combine satellite images with radar data to obtain improved estimates.
3.3 Measurements Related to Snow
Many big rivers of the world are snowfed and their catchments receive a part of precipitation in
solid form or snow. Three variables related to snow are important in water resources studies:
snow depth, area of snow cover, and snow water equivalent. Such data are used to estimate water
content of the snowpack and are of immense use in planning of irrigation and hydropower
projects, reservoir regulation, disaster management, etc.
While measuring snowfall, care should be taken to measure only fresh snow. For this
purpose, a suitable patch on the ground is cleared beforehand or the top of the existing snow
surface is covered with some material (wood, tarpaulin, etc.) so that old and fresh snow can be
demarcated and the depth to this surface is measured. Note that if there is a layer of old snow, it
would not be right to compute the depth of a fresh snowfall from the difference between two
consecutive measurements of the total depth of (fresh and old) snow since old snow keeps on
settling continuously.
A snow pillow is a flat circular container that is filled with non-freezing liquid and is
used to measure the water equivalent of snow. The weight of the snow on the pillow can be
measured by hydrostatic techniques or pressure transducer. The depth of snow may also be
measured in a standard container of uniform cross-section after levelling the snow (without
compressing it). Snow pillows of various dimensions and materials are used to measure the
weight of snow that accumulates. The most common pillows are flat circular container of 3.7 m
diameter made of rubberized material filled with a non-freezing liquid. The pillow is installed on
the surface of the ground, flush with the ground, or buried under a thin layer of soil or sand. To
prevent damage to the equipment and to preserve the snow cover in its natural condition, it is
recommended that the site be fenced. Under normal conditions, snow pillows can be used for 10
years or more.
Hydrostatic pressure inside the pillow is a measure of the weight of the snow on the
pillow. Hydrostatic pressure can be measured by means of a float operated water-level recorder
or a pressure transducer. Data from snow pillows are most reliable when the snow cover does not
contain ice layers, which can cause bridging above the pillows.
In regions receiving heavy snow, commonly snow stakes are installed at different
(accessible) representative locations and are observed to find the depth of fallen snow. Snow
stakes should be installed such that they can be inspected easily from a distance. It needs to be
ensured that the site is the representative of the area and is free from human disturbance. Snow
stakes should be painted in white to minimize increase in melt of surrounding snow. The stakes
are clearly graduated in metres and centimetres and in inaccessible areas, these are provided with
crossbars so that the reading can be taken from a distance.
The water equivalent of snow is the depth of water that would be obtained when the snow
is melted. It can be determined by melting the snow collected in a gauge and measuring the melt
water in the same way as rainfall is measured. It is better to take several representative samples.
Alternatively, cylindrical samples of fresh snow are collected with a suitable snow sampler and
are either weighed or melted. Weighing-type recording gauge are used to determine the water
content of snowfall. During periods of snowfall, funnels of the gauges should be removed so that
the precipitation can fall directly into the receiver. Isotope techniques are also used to estimate
snow water equivalent. However, the instruments used are expensive and adequate safety
measures need to be observed.
In important catchments which receive considerable snow fall, snow surveys are made on
permanent snow courses. A snow course is a permanently marked line where snow surveys are
made. The snow courses should be carefully selected and should be representative of the
catchment area. These should be located in accessible areas where snow falls to the ground
without being intercepted by vegetation, and the site should be protected from strong winds.
During surveys, snow samples are taken and depths are measured at various places at a number
of points along the course.
Processing of data is carried out with two major objectives: to examine the data for errors
and remove them, and to organize the data in a form that is appropriate for subsequent analysis.
The rainfall data may have following errors:
Entries on the wrong day - shifted entries,
Entries made as accumulations,
Missed entries, and
Rainfall measurement missed on days of low rainfall.
Example 3.1: The daily precipitation reported from a station is 358.6 mm and the statistics of
the reporting station are:
The reported daily rainfall value (358.6 mm) is more than the normal monthly rainfall of
the corresponding month and is, therefore, doubtful. This value is more than the mean maximum
1-day ( x ) and ( x +) which are 261.6 mm and 261.6 mm, respectively. The reported daily value
is compared with 1-day Probable Maximum Precipitation (PMP) value. Since it is less than the
PMP (370.8 mm), it is considered possible. At this stage, it is not advisable to reject this value. It
should be flagged and further checked by spatial consistency.
The internal consistency or self consistency checks are applied by using statistical
information based on historical data of the station and current data in case of short duration
rainfall. An example of checking the data by the internal consistency is given here.
The hourly rainfall reported during 3rd hour is suspected. Check its reasonability.
Solution: Although the hourly rainfall reported during 3rd hour is suspected, it could not be ruled
out. We examine the total rainfall for 1-3 hours to check the value in the 3rd hour. The 3-hourly
total rainfall was reported as 54.1 mm and this indicates that the rainfall in the 3rd hour could be
35.3 mm (assuming that the rainfalls in the 1st and 2nd hours are correct). When the 6 hour total
rainfall is reported as 117.4 mm, the value of 35.3 mm is confirmed for the 3rd hour. Further
checking for the erroneous value is carried out similarly.
This example shows that it is a good practice to report the hourly values as well as the
sum of block of hours so that possible errors in conveying the data can be detected and corrected.
where, X 1 and X 2 are the arithmetic mean of the two samples of size n1 and n2 respectively; S
is the unknown population standard deviation estimated from the samples variances s1 and s2 as:
(n1 1) s1 (n 2 1) s 2
S
n1 n 2 2 (3.3)
If the value of the statistic t is less than the tabulated value of Student’s distribution at some
chosen significance level α and (n1+n2-1) degrees of freedom then the hypothesis that “the means
of both the samples are not significantly different” may be accepted at the chosen significance
level.
If the computed F is less than the tabulated value of F distribution at some chosen significance
level α, and n1-1 and n2-1 degrees of freedom then the hypothesis that “the variances of both the
samples are not significantly different” may be accepted at the chosen significance level.
Statistical tests such as t-test and F-tests are discussed in Module 10 in detail.
Rain gauge station failure: Depending on the nature and type of failure of rain gauge, data of
self recording rain gauge and/or the ordinary rain gauge can be utilised to correct and complete
the missing data.
Accumulated rainfall data: If the data of ordinary rain gauge is read after several days, the
recorded total represents accumulated rainfall over a period of several days. Such accumulated
totals can be distributed over the period of accumulation based on the rainfall recorded at the
neighbouring stations over the same period.
Long-term shift in rainfall data: Consistency of a rain gauge station data is examined to
determine if the conditions relevant to recording rainfall at a station have undergone a significant
change during the period of record. In case there is significant change inconsistency would arise
in the rainfall data of that station. This inconsistency would be noticed from the time the
significant changes have taken place. Some of the common causes for inconsistency of record
are:
The double mass curve technique is used to check the consistency of a rainfall series.
This technique is based on the principle that the data are consistent when each recorded data
comes from the same parent population.
To carry out the test, a group of stations in neighborhood of the index station (for which
the consistency test is required) is selected. The data of the annual rainfall of the index station
(X) and also the average rainfall of the group of base stations covering a long period is arranged
in the reverse chronological order, i.e. the latest record as the first entry and the oldest record as
the last entry in the list. The accumulated precipitation of the station X and the accumulated
values of the average of the group of base stations are calculated starting from the latest record.
Individual cumulative rainfall values of the index station are plotted against mean cumulative
values for various consecutive time periods. A break in the slope of the resulting plot indicates a
change in the precipitation regime of the index station. The precipitation values at the index
station prior to the change of regime are corrected by using the initial slope of the line as follows:
PX (i ) m C A (i ) C X (i 1) (3.6)
where PX(i) is the corrected rainfall at the index station X during year i; m is the initial slope of
the line; CA(i) is the cumulative mean rainfall up to the year i; CX(i-1) is the cumulative rainfall
at the index station X till the year i-1.
As example, consider the double mass plots of Bhopal and Basoda (Figures 3.5 and 3.6)
which have been plotted using the annual rainfall data of eight other neighbouring stations. It can
be observed from Fig. 3.5 that the change in slope of the mass curve for Bhopal is insignificant,
whereas, the same is significant in the case of Basoda (Fig. 3.6). Hence, corrections are needed
for Basoda data.
80000
70000
60000
Cum. at Bhopal, mm
50000
40000
30000
20000
10000
0
0 10000 20000 30000 40000 50000 60000 70000 80000
Cum. areal average, mm
Figure 3.5: Double mass curve for Bhopal. X-axis is the accumulated mean precipitation at group
of stations and on the y-axis is accumulated precipitation at Bhopal.
100000
Original Corrected
Cum. at Basoda, mm
80000
60000
40000
20000
0
0 20000 40000 60000 80000 100000
Cum. Average, mm
Fig. 3.6 Double mass curve for Basoda. X-axis has the accumulated mean precipitation at group
of stations and on the y-axis is accumulated precipitation at Basoda. Both original and corrected
curves have been plotted.
3.4.4 Spatial Consistency Check
Rainfall data exhibit some spatial consistency and this forms the basis of investigating the
observed rainfall values. An estimate of the interpolated rainfall value at a station is obtained on
the basis of the weighted average of rainfall observed at the surrounding stations. If the
difference between the observed and the estimated values exceed the expected limiting value,
such values are considered as suspect and are flagged for further investigation and ascertaining
the possible causes of departures.
Spatial consistency checks for rainfall data are carried out by relating the observations
from surrounding stations for the same duration with the rainfall observed at the station. This is
achieved by interpolating the rainfall at the station under question with rainfall data of
neighboring stations. The station being considered is called the test station. The interpolated
value is estimated by computing the weighted average of the rainfall observed at neighboring
stations. Ideally, the stations selected as neighbors should be physically representative of the area
in which the station under scrutiny is situated. The following criteria are used to select the
neighboring stations:
(a) The distance between the test and the neighboring station must be less than a specified
maximum correlation distance;
(b) too many neighboring stations should not be considered for interpolation; and
(c) to reduce the spatial bias in selection, it is advisable to consider an equal number of
stations in each quadrant.
Example 3.3: Rainfall reported at a group of five stations (see Fig. 3.7) is as follows.
Chando Kundla
Virpur Idar
Sabli
During the quality control process, the data at Idar is identified as doubtful. Check this
data for spatial consistency.
Solution: The rainfall at Idar is estimated using the distance power method and compared with
the observed value. From the four quadrants around Idar (Fig. 3.6), the station nearest from each
quadrant is selected for estimation of rainfall at Idar. Using the reference coordinate system, the
distance of each of the estimator stations from Idar is determined and the rainfall at Idar is
estimated.
Since the observed value is very much different from the estimated value, it is rejected
and replaced by the estimated value. Note that there is a possibility that the decimal point was
wrongly placed while recording the data at Idar.
Numerous methods of computing areal rainfall from point measurements have been
developed. While using precipitation data, one often comes across missing data situations. Data
for the period of missing rainfall could be filled using various techniques. Due to the spatial
structure of precipitation data, some type of interpolation making use of the data of nearby
stations is commonly adopted.
Let the precipitation data be available at n stations, spread over an area and Pi be the
observed depth of precipitation at the ith station. Using a linear interpolation technique, an
estimate of precipitation over the area can be expressed by
n
P * = Pi W i (3.7)
i=1
where Wi is the weight of the ith station. The spatial averaging techniques differ in the method of
evaluation of these weights. Weights of an optimal interpolation technique are decided such that
the variance of error in estimation is the minimum.
The most commonly used methods are for Spatial Averaging of Precipitation Data:
The choice of the method is depends on the quality and nature of data, importance of use and
required precision, availability of time and computer. Some of the commonly used methods are
described below.
P W
i=1
i i
P= (3.8)
n
where P is the average catchment precipitation from the data of n stations, Pi is the precipitation
at station i, and Wi is the weight of ith station. If the gauges are relatively uniformly distributed
over the catchment and the rainfall values do not have a wide variation, this technique yields
good results.
in which Wi = Ai/A, where Ai is the area represented by the station i and A is the total catchment
area. Clearly, the weights will sum to unity.
An advantage of this method is that the data of stations outside the catchment may also be
used if these are believed to help in capturing the variation of rainfall in the catchment. The
method works well with non-uniform spacing of stations.
A major drawback of this method is the assumption that precipitation between two
stations varies linearly and the method does not make allowance for variation due to orography.
In this method, the precipitation depth changes abruptly at the boundary of polygons. Also,
whenever a set of stations are added to or removed from the network, a new set of polygons have
to be drawn. The method fails to give any idea as to the accuracy of the results. If a few
observations are missing, it may be more convenient to estimate the missing data than to
construct the new set of polygons.
Example 3.4: For a catchment, the rainfall data at six stations for July month along with their
weights are as given in Table 3.4. Find the weighted average rainfall for the catchment by using
the Thiessen polygon method.
Solution: Using the observed rainfall and station weight, weighted rainfall at each station is
computed. Summation of these values gives the weighted average rainfall for the catchment. The
computations are shown in Table 3.1.
Table 3.1 Estimation of the mean areal rainfall by the Thiessen polygon method.
Berasia
Vidisha
Gairatganj
Raisen
Bhopal
Sehore
Fig. 3.8 The Thiessen polygon method for computing the mean areal rainfall.
Computers are being used to draw isohyetal maps these days, by using special software.
As an example, the isohyetal map for an area is shown in Fig. 3.8. The total depth of
precipitation is computed by measuring the area between successive isohyets, multiplying this
area by the average rainfall of the two ishohyets, and totaling. The average depth of precipitation
is obtained by dividing this sum by the total area. The average depth of precipitation (Pi) over
this area is obtained by:
n
P
i=1
i Ai
P= n
(3.10)
A
i=1
i
where Ai is the area between successive isohyets and Pi is the average rainfall between the two
isohyets.
Fig. 3.8 The
T isohyetaal method for computing
g the mean arreal rainfall..
Iso
ohyetal rangee Averrage value Areaa Volume
(mm) (mm)
( (km2 ) (105m3)
52
21.0 - 500.0 510.5
5 70 357.4
50
00.0 – 300.0 400.0
4 530 2120.0
33
38.0 – 300.0 319.0
3 100 319.0
20
00.0 – 300.0 250.0
2 20800 5200.0
15
58.0 – 200.0 179.0
1 28200 5047.8
56000 13044.2
13044.2
Averaage catchmen
nt rainfall = 2232.9mm
5600
Solution
n: For each issohyet, the average
a value is worked out (the maxximum obseerved rainfalll was
108 cm and
a the minimum 38 cm m). This, multiplied by thhe area encloosed by that isohyet givees the
volume of
o rainfall foor that isohy
yet. Now thhe volumes ffor differentt isohyetals are summedd and
divided by the area of the catchment to get average catchment rainfall. The computations are
shown in Table 3.3.
NR
NR A
* Pi
i=1 i (3.11)
PA=
n
where Pi is the rainfall at surrounding stations, NRA is the normal monthly or seasonal rainfall at
station A, NRi is the normal monthly or seasonal rainfall at station i, and n is the number of
surrounding stations whose data are used for estimation.
Example 3.6: A catchment has four rain gauge stations A, B, C & D. Normal monthly rainfall at
these stations is known. The observed rainfall at the stations B, C & D for a storm event is
known and is given in the following table. Find the missing rainfall at station A.
Solution: The ratio of normal rainfall at station A to normal rainfall at station i or NRA/NRi has
been calculated and is given in table below.
Station A B C D
Normal Rainfall (mm) 331.3 290.8 325.9 360.5
Event Rainfall (mm) ? 98.9 120.5 110.0
NRA/NRi 1 1.14 1.02 0.92
P /D
i=1
i
2
i
(3.12)
PA = n
1/ D
i=1
2
i
Note that the weights go on reducing with distance and approach zero at large distances. A major
shortcoming of this method is that the orographic features and spatial distribution of the variables
are not considered. The extra information, if stations are close to each other, is not properly used.
The procedure for estimating the rainfall data by this technique is illustrated through an example.
If A, B, C, D are the location of stations discussed in the example of the normal ratio method, the
distance of each estimator station (B, C, and D) from station (A) whose data is to be estimated is
computed with the help of the coordinates using the formula:
where x and y are the coordinates of the station whose data is estimated and xi and yi are the co-
ordinates of stations whose data are used in estimation.
Example 3.7: Using the data of Example 3.6, estimate rainfall at station A using the distance
power method.
Solution: Since the coordinates of the stations are known, their distances from station A can be
calculated. The weights 1/Di2 are then computed for each station and the rainfall at station A is
estimated as follows:
570.7 10 3
Rainfall at station A = 113.9 mm
5.0110 3
A mass curve is a graphical plot of accumulated rainfall at a station versus time. Mass
curves of accumulated rainfall at (non-recording) daily stations and recording stations can be
prepared by plotting the accumulated rainfall values against time for the storm duration under
analysis. Now the mass curves of the recording rain gauge stations are compared with those of
the non-recording stations to decide which recording rain gauge is representative of which of the
non-recording rain gauge for the purpose of distributing daily rainfall into hourly rainfall.
Assume that the daily rainfall is observed at 0800 hours. For converting the daily rainfall
into hourly rainfall, the hourly rainfall from 0800 hr to 0800 hr for consecutive days is
accumulated and the rainfall during each hour is expressed as a ratio of the total rainfall during
24 hours (0800 to 0800). These ratios for the identified station are used to distribute the daily
rainfall for the corresponding duration at the non-recording rain gauge stations.
The procedure for distribution of daily rainfall at non-recoding rain gauge stations into
hourly rainfall is explained with the help of an example.
Example 3.8: Daily rainfall data (mm) of four stations for the period 28th August to 30th August
1973 is given below:
Station name 28/08/1973 29/08/1973 30/08/1973
Pendra Road 65.3 80.8 1.1
Dindori 23.2 58.4 1.0
Ghansore 42.0 98.0 20.6
Niwas 30.4 203.0 33.4
Hourly rainfall data of two self-recording rain gauge (SRRG) stations (Mandla and Jabalpur) is
given in Table 3.5 for the period 27-29 August 1973. Disaggregate the daily rainfall data into
hourly values.
300 JABALPUR
Cumulative Rainfall
200
100
MANDLA
(mm)
0
27.8.73 28.8.73 29.8.73
Fig. 3.9 Mass curve of hourly rainfall.
300
NIWAS
GHANSOR
200
PENDRA ROAD
Rainfall (mm)
Cumulative
DINDORI
100
0
27.8.73 28.8.73 29.8.73 30.8.73
Solution: The hourly rainfall data at each of the two SRRG stations is plotted on a graph to
prepare the mass curve of hourly rainfall as shown in Figure 3.9. The daily rainfall data at each
of the four stations is also cumulated and plotted on a graph in Figure 3.10. The points are joined
to form the mass curve of daily rainfall. It may be seen that the mass curve of hourly rainfall has
many kinks as this has been prepared using hourly data. The mass curves of daily rainfall are
compared with those of hourly rainfall to determine which of the daily rain-gauge stations are
represented by which of the SRRG stations. Based on the match of curves, it was decided that
the daily rainfall stations at Pendra Road and Dindori are represented by Mandla while Jabalpur
represents Niwas and Ghansore stations.
To convert the daily rainfall into hourly rainfall, the hourly rainfall from 0900 hr on a day
to 0800 hr for consecutive days is cumulated and the rainfall during each hour is expressed as a
ratio of the total rainfall during 24 hours (0900 to 0800). These ratios are used to distribute the
daily rainfall for the corresponding duration at those rain gauge stations, which are represented
by the SRRG. The daily rainfall distributed for 1 day, i.e. corresponding to 28-29 August 1973 is
given in Table 3.6.
Remarks
It is very well recognized that the areal precipitation estimates derived on the basis of point
rainfall observations from a network of sparsely and unevenly distributed precipitation gauges
could only be regarded as an index of rainfall.
Radar sensed echo intensity reflected by precipitation could provide useful areal estimates of
precipitation over areas with few or no precipitation gauges. Precipitation observed by a dense
network of precipitation gauge specially set up for the purpose on a temporary basis could be
related to echo intensities as sensed by radar to develop reasonable relationships for the area of
interest.
Table 3.6: Daily rainfall distributed into hourly rainfall (28-29 August 1973)
Estimation of mean areal precipitation in mountainous areas by the conventional methods does
not yield the desired results because of the influence of orography, aspect, storm orientation, etc
on the precipitation in areas with pronounced orography. A better way of estimating areal rainfall
in mountainous regions is through the use of isopercental technique.
While designing a dam, it is necessary that the outlet capacity is large enough to safely pass a
flood of certain magnitude. This critical flood is known as the design flood for the structure. The
type of the hydraulic structure is the main criterion to decide the design flood hydrograph. For
this purpose, the structures are classified as:
If long-term runoff data are not available, rainfall data which are generally available for a
longer period are used to estimate the design storm. This design storm is used with a suitable
rainfall-runoff model to obtain a design flood. The rainstorm analysis is the first step in the
design storm estimation procedure. The design storm (rainfall) is a magnitude of rainfall and its
distribution which is developed for the design of specific types of structures. It has three
components, namely, the rainfall amount, the areal distribution of rainfall, and the time
distribution of rainfall.
The Probable Maximum Precipitation (PMP) is defined as “theoretically the greatest
depth of precipitation for a given duration that is physically possible over a given size storm area
at a particular geographical location at a certain time of year”. The Standard Project Storm is
defined as that rain storm which is reasonably capable of occurring in the region of problem
basin. It is generally the most severe rain storm which has occurred in the region of the basin
during the period of available records.
The data on volumes of precipitation during severe storms is important for examining and
studying storms suitable for design purposes. Such information is generally presented in the
form of tables of the maximum average depth of storm precipitation over various standard area
sizes, such as 100 km2, 500 km2 etc. These data are known as the Depth-Area-Duration data and
they could be presented in tabular as well graphical form. WMO (2008) has described
procedures to estimate the probable maximum flood and standard project flood.
Closure
Undoubtedly, precipitation remains the most important and frequently used hydro-
meteorological in studies dealing with water resources. Precipitation is also the variable whose
longest measuring data series are available at a large number of stations. This module gives a
detailed coverage of acquisition and processing of precipitation data.
References
Jain, S.K., and Singh, V.P. (2003). Water Resources Systems Planning and Management.
Elsevier, Amsterdam.
Subramanya, K., Engineering Hydrology, Tata McGraw Hill Education Pvt. Ltd., New Delhi.
WMO(1982). Concepts and techniques in hydrological network design. Operational Hydrology
Report No.19, WMO No. 580, World Meteorological Organization, Geneva.
WMO(1982). Methods of correction for systematic error in point precipitation measurements for
operational use. Operational Hydrology Report No. 21. WMO No. 589. World
Meteorological Organization, Geneva.
WMO(1983). Guide to Meteorological Instruments and Method of Observation. WMO No. 8.
World Meteorological Organization, Geneva.
WMO (2008). Guide to Hydrological Practices. Volume I: Hydrology – From Measurement to
Hydrological Information. WMO No. 168. World Meteorological Organization, Geneva.
Module 4
Measurement and Processing of Streamflow Data
Keywords: Streamflow Data, Station Control, River Stage, Depth Sounding, Water Level
Recorder, Velocity-Area, Float, Current Meter, Slope-Area, Moving Boat Method, Weirs,
Flumes, Dilution, Electromagnetic Method.
For effective and efficient water resource assessment and management and for proper
river basin planning, development of flood forecasting system, etc. reliable, accurate, processed
and easily accessible data base containing data of hydrologic variables is a pre-requisite.
The most important hydrological data for management of surface water resources
pertains to streamflow. Streamflow has served as the lifeline for mankind and continues to do so.
Its importance is also relatively more, since this source is visible in contrast to ground water
which is hidden. Streamflow records are primarily continuous records of flow passing through a
particular location on the stream. Streamflow data are analysed to determine the magnitude and
variability of surface waters. These records are input in planning, design, and operation of
surface water projects and are also used in design of bridges and culverts, flood forecasting
systems, and flood plain delineation.
When the channel cross-section is regular and of a known shape, an estimate of the
discharge can be obtained by simply measuring the averaged flow either across the channel at
some depth or over depth in the center of the channel. Fig. 4.1 shows the classification of various
methods of discharge measurement.
The sections where river measurements are carried out are known as stream gauging
stations. A network of these stations is established to collect data about river flows of a region.
The location of gauging sites primarily depends upon the purpose of data collection. If the site is
needed for a specific project, the general location will be in the vicinity of the project. However,
if a network of gauging stations is to be established to study the general hydrology of a region
and for planning and design of various projects than careful planning is required to identify
locations so that optimum information is obtained for the resources deployed in the data
collection.
River gauging stations are of three types: basic data stations, operational stations, and
special stations. The basic stations are established to collect data for a variety of uses, including
planning and design of projects, and to understand the hydrological characteristics of the region.
The operational stations collect data to run projects and issue forecasts. The objective of special
stations is to meet specific data needs that may arise in cases, such as research, project
investigation, special studies, legal cases, etc. Their operation is terminated when the specific
need is fulfilled. Sometimes, auxiliary stations are set up to augment the database from the
network or region. An auxiliary station may, for example, record only the peak discharge which
occurs at that site during a certain (usually limited) period.
Conventional Advanced
Methods Techniques
Direct Indirect
The number of gauging sites depends on the cost of installation and operation, the value
of the data, watershed size, degree of development, objective of data collection, accuracy,
hydrologic characteristics, etc. Some of these factors are interrelated. The streamflow data that
are of immense use in water resources are river water level or gauge and discharge. A continuous
observation of the river water level or stage may be made with comparative ease and economy.
At important stations, the stage is measured at short intervals and discharge is measured once or
twice each day. At less important stations, only stage measurements are made regularly.
A continuous measurement of discharge in a natural channel is comparatively difficult,
time consuming, expensive and requires special skills. Therefore, the discharge at a site is
measured less often and is estimated by indirect methods. Fortunately, there exists a relation
between stage and discharge at a section. This relation is termed as stage-discharge relationship
or rating curve. This relationship is used to transform the observed stages into discharges. Note
that the reliability of such discharge records is dependent on the reliability of stage data and the
correctness of the stage-discharge relation.
At many sites, the discharge is not a unique function of stage; variables other than stage
must also be simultaneously measured at such sites to correctly estimate discharge. For example,
if variable backwater occurs at a site, the data on water surface slope is required. The slope can
be measured by installing an auxiliary gauge downstream. The rate of change of stage is an
important variable where the flow is unsteady and channel bottom slopes are flat.
When any change in the physical characteristics of the channel downstream to the control
has no effect on the flow at the gauging section, such control is termed as section control.
Natural or artificial local narrowing of the cross-section (waterfalls, rock bar, gravel bar) creating
a zone of acceleration (Fig. 4.2) are some examples of section controls. The section control
necessarily has a critical flow section at a short distance downstream.
Fig. 4.2 A section co
ontrol
A cross sectiion where non acceleratiion of flow occurs or w where the accceleration iis not
sufficientt to prevent passage of disturbances from the ddownstream m to the upstrream directiion is
called as a channel control.
c The length of th he downstreaam reach off the river afffecting the rrating
curve deppends on thee normal or equilibrium
e depth he andd on the eneergy slope S [L ∞ he/S, w
where
5/3 0.5
5 0.5 3/5
he follow
ws from Man nning’s Q =KmBhe S (wide rectaangular channnel) so he = [(Q/KmS ) ].
The length of channeel effective as
a a control increases
i witth dischargee.
An
A artificial control
c is sp
pecifically co
onstructed too stabilize thhe relationshhip between stage
and dischharge. Thesee include weeirs and flum mes, dischargging under ffree flow coonditions. Naatural
section controls
c incllude a ledgee of rock accross a channnel, the brinnk of a watterfall, or a local
constrictiion in width (including bridge
b openiings).
Natural
N contrrols can have a range off geometry aand stabilityy. Some natuural controlss may
have a siingle featuree such as a rock
r ledge across
a the chhannel at thee crest of a w
waterfall, theereby
forming a complete control.
c Thuus complete control
c goveerns the stagge-discharge relation oveer the
whole rannge of stagee experienced d. However, in many casses, station ccontrols are a combinatiion of
section control
c at low
w stages andd a channel control at hiigh stages annd are thus called comppound
or compllex controls.
Where
W the geeometry of a section and d the resultiing stage-disscharge relattionship doees not
change with
w time, it is described d as a stable or permanennt control. SShifting conttrols changee with
time thesse may be seection contro
ols such as boulder, gravvel or sand riiffles which undergo perriodic
or near continuous
c sccour and depposition, or they may bee channel coontrols with erodible bedd and
banks.
The amount of gauging effort and maintenance cost to obtain a record of adequate quality
is much greater for shifting controls than for permanent controls. Since stage discharge
observations require significant effort and money, it is always preferred to select a gauging site
with a section or structure control. A complete control has many advantages: permanence, easy
installation and running of the gauging site, and favourable conditions for current meter
measurements. However, a complete control is not practicable in many cases and one may have
to be content with either channel control or a compound control.
The International Standard Organisation (ISO) has brought out a large number of
standards dealing with measurement of liquid flow in open channels. Technical committee
TC113 of ISO deals with this topic. The publication ISO (1983) is a useful collection of
standards dealing with various aspects of streamflow measurement. WMO has also brought out
many publications related to streamflow measurement. Of course, the individual countries have
their own standards. The topic is covered in detail in Herschey (1986 and 1995) and Boiten
(2008). Site requirements for measurement of discharge using current meters or floats are given
in ISO 748.
River stage
Gauge datum
Figuree 4.3 Cross section of a river
r with gaauge datum aand river stagge.
Usually,
U it iss difficult annd expensivve to measurre the dischharge of a rriver directlyy and
continuouusly. Therefore, one frequently
fr measures
m thee water stagge which iss easy and then
converts the stage vaalues to disch harges. Recoords of stagee are used w
with a stage-ddischarge rellation
to obtain
n the record ds of stream m discharge. Clearly, thhe reliabilityy of the disscharge recoord is
dependennt on the reliiability of th
he stage recorrd and the sttage-discharrge relation.
The
T water lev vel data is measured
m using a varietyy of equipmeent: staff gauuges, autographic
water lev
vel (chart) reecorders, an
nd digital typ
pe water levvel recorderss. These cann be classifiied in
two broad categories: (1) Non-reecording, maanual gaugess, (2) Recordding gauges.
4.2.1 Non-recordin
N ng river sta
age gauges
The simp plest way to measure riv ver stage is by
b the use of a staff gauuge which iss basically a scale
installed such that a portion of itt is always immersed
i inn the water. TThis gauge m may be a veertical
scale attaached to a co
olumn, pillaar, bridge pieer, or other sstructure thaat extends innto the low-w
water
part of thhe channel. Such staff (manual) gaauges are sim mple and innexpensive bbut must bee read
frequentlly.
Where
W the raange of wateer levels ex xceeds the caapacity of a single verttical gauge other
gauges may
m be instaalled in the linel of a cro
oss-section nnormal to thhe direction of flow (seee Fig.
4.2). Thee scales on suuch a series of stepped staff
s gauges should overllap by not leess than 15 ccm. A
ramp gau uge consistss of a scale marked on or securelyy attached too a suitable inclined surrface,
which clo osely follow
ws the contou ur of the rivverbank. Thee ramp gaugge may lie onn one continnuous
slope orr on more slopes. It should lie
on the line
l of a crross-section normall to
the
directioon of
flow.
A
staff or
ramp gauge is an inexpensive, simple, and reliable method of measuring water level. By using it,
water level can be measured by relatively unskilled staff. Ramp gauges amplify surges and
ripples but provides the opportunity to of a higher resolution. A staff gauge can only be used for
spot measurements. It is difficult to obtain readings in the field with a true resolution higher than
±5 mm.
Figure 4.4: Staff gauge – pictorial view (left) and markings (right)
Staff gauges can be either vertical or inclined. Vertical staff gauges are normally
porcelain enameled iron sections with a scale (with graduations of 5 or 10 mm) marked on or
securely attached to a suitable vertical surface graduated every 10 mm. The vertical staff gauge
is used as an inside reference gauge (if installed in a well), or as an outside gauge if installed in
the stream. Where the water level of the river varies over a large range, observations using a
single staff gauge might be difficult. In such cases, the gauge consists of stepped sections
installed at different locations in a line normal to the flow. Each of these stepped gages should
refer to the common datum and they should overlap by not less than 15 cm to ensure continuity
of readings and confirm their consistency with each other. Staff gauges are manually read,
generally each day in the morning in lean season and at (multi) hourly intervals during high
flows.
Manual gauges are to be read at fixed times whereas the recording gauge provides a
continuous data of the variation of stage. In manual observation, commonly the stage is read
daily in fair weather and (multi)hourly in monsoon season. In automatic recording stations, a
continuous record of stage is obtained by utilizing water level sensors interfaced with a analogue
recorder (chart) or a digital recorder (logger or telemetry).
The advantages of the non-recording gauge are low initial cost and ease of installation.
The disadvantages are the need for an observer (and even then data will observed at limited
times) and less accuracy. For a long-term operation, the advantages of a recording gauge
outweigh those of a non-recording gauge. Sometimes, an automatic and a non-recording gauge
are maintained together because the electro-mechanical recording gauge equipments are liable to
breakdowns. The values of stage may be required as a single instantaneous measurement, as a
short series of instantaneous measurements or as a continuous or practically continuous record of
the fluctuations of stage.
Even when a gauge recorder is used, the observer should visit the station from time to
time to ensure satisfactory performance of the sensor and recorder. Observer should note the
time and date of such checks along with the staff gauge and recorder values. It is essential that
the staff gauge itself is maintained such that it is safely accessible and legible to the observer.
Further, the observer should always visit the site following a major high flow event to ensure
continued measurement and recording of data.
For smaller depths and velocities, a wading-rod is used; for greater depths, a sounding
line is used. A sounding rod should be as lightweight as possible but sufficiently strong to
withstand the force exerted by flowing water without undergoing significant deflection or
vibration. It should remain straight and vertical during use; should not cause significant heading
up of water and should not penetrate into the channel bed.
To make observations in flowing water with the help of suspension equipment, the
measuring equipment must be placed at the point of measurement in such a way that it does not
cause appreciable disturbances, irrespective of the depth of water and velocity of flow.
The water level recorders are generally of shaft-angular-input type, and the angular
rotation of the shaft is recorded. The depth of water surface is sensed for automatic recording by
a float in a stilling well which follows the rise and fall of the water level. A gas-purge system
that transmits the pressure head of water in a stream to a manometer is known as a bubble gauge.
A water level recorder gives a continuous record of the water level on a chart from which
values are manually extracted at desired intervals. The data from a digital water level recorder
are stored in an electronic memory and these are downloaded to a computer. The data can either
be recorded at equal intervals of time, usually at (fraction of) an hour interval, or at only those
instants when there is a change in water level by more than a pre-set amount.
A float gauge can provide a direct record of river stage and no external energy source or
battery is required. It provides almost uniform resolution throughout the range and good
accuracy. But being a mechanical device, it is subject to errors from hysteresis and friction.
Construction and maintenance of stilling well is expensive, particularly in places where rivers
carry large amount of sediments.
A chart recorder produces a continuous hard copy record. Quality of the recording
mechanism affects the accuracy of the chart record. Subsequently, the data is manually converted
to digital format; this process is labour insensitive, time consuming, and may introduce errors.
Figure 4..5 Float typee recorder
In the direct determination methods either discharge itself is measured or some variable
on which discharge depends is measured. The commonly used methods are: velocity-area
methods, dilution techniques, electromagnetic method, and ultrasonic method. The first two are
described here.
Q=v*A (4.1)
Since the velocity of flow at a cross-section varies laterally and with depth, it is not
enough to measure the velocity at a single point. Depending on the accuracy required, the width
of the stream is divided into a number of vertical portions (Fig. 4.7). In each of these portions,
the velocity is measured at one or more points along the depth to get a representative velocity in
that portion. The area of the individual portion can be easily calculated if the bed profile and
stage are known. The velocity may be measured by a conventional method (for example, float or
current meter) or by an advanced procedure, for example, the moving boat technique.
Fig. 4.7 Schematic sketch
s for a velocity-area
v a station.
Water
W particlles in a riverr cross-section travel at different veelocities. Thhey are subjeect to
friction as
a they com me into contact with th he sides andd bottom off the channel. Due to these
frictionall effects, waater flows faastest near thhe surface annd center off the channeel (away from m the
immediatte frictionall influences)). A typicall velocity pprofile is shhown in Figg. 4.8 where the
velocity varies as a parabola fro om zero at the channel bottom to a maximum m at (or nearr) the
surface. Also
A shown are the contours of equaal velocity inn typical riveer cross sectiion.
Fraction of Depth (%)
Velocity (m/s)
Figure 4..8 Typical veelocity profile (top) and velocity disttribution (isoovels) in a cchannel (botttom).
4.4.1 Velocity
V Meaasurement by b Floats
A float iss a distinguisshable article that floats on the wateer surface, suuch as a woooden log, a pplastic
bottle paartly filled with
w water, or o branch off a tree. Surrface or neaar-surface flooats are useed for
streamfloow measurem ment. Norm mally, these are
a wooden cylindrical rods of nearrly 0.5 m leength.
They aree shaped succh that they float nearly y vertically with one thhird of the llength protruuding
above thee water surfaface. The floaats are paintted in bright colours for easy identiffication in m
muddy
or turbullent water. At
A times, a floating
f article such as a lump of vvegetation m may be usedd as a
float.
The velocity of the float is equal to the distance between the two cross-sections divided
by the time taken by the float to cover this distance. The mean velocity in the vertical is
determined as the float velocity multiplied by a coefficient which varies between 0.80 and 0.85.
This coefficient depends on the shape of the velocity profile of the river and the depth of
immersion of the float. Discharge for the segment is calculated as the segment area multiplied by
the segment velocity. Total discharge is obtained as the sum of segment discharges and mean
velocity of flow is computed as the total discharge divided by total area.
Floats can rarely be positioned at the desired location in the cross section. If these are
thrown manually, it is difficult to throw them very far and hence only the velocity of water near
the banks can be measured. Float measurements are not very reliable and this method is normally
restricted to an emergency or to measure high discharges when current meter is not available or
can’t be used. Hence, floats should be used only when it is not possible to use any other better
method.
In
n the propelller-type currrent meters (Figure 4.99), a propeller rotates abbout a horizzontal
axis due to force appplied by flowwing water. The revoluttions per timme interval aare recordedd. The
relation between
b revoolutions per second N of
o the currennt meter and the water veelocity v is ggiven
by follow
wing equation
v = a + b*N (4.2)
where b = constant ofo proportionnality and a = starting vvelocity or velocity requuired to overcome
mechaniccal friction. These consstants differ from one ccurrent meteer to the othher as a resuult of
manufactturing variattions as well as changee with time due to weaar and tear. Therefore, each
current meter
m should be recalibraated periodiccally.
The
T velocity measuremen nts made byy using curreent meters arre usually cllassified in tterms
of the meeans used to o cross the stream
s durin
ng measurem ments, such aas wading, ccableway, brridge,
or boat. Wading
W is possible
p in sm
mall streamss of shallow
w depth onlyy; the currentt meter is heeld at
the requiisite depth below
b the su
urface by ann observer wwho stands inn the water.. In narrow well-
defined channels,
c a cableway is stretched from bank to bank weell above thhe flood levvel. A
carriage moving
m overr the cablewway serves ass the observaation platforrm. Bridges are advantaggeous
from the viewpoint of o accessibillity and tran
nsportation, aalthough theese are not tthe best locaations
from hyddraulic pointt of view. Th he velocity measuremen
m nt is perform
med on the doownstream oof the
bridge to
o minimize th he instrumennt damage due
d to drift aand knock aggainst bridgee piers. Boatts are
most satiisfactory for measuremen nts in wide rivers.
r
The
T section liine at the gaauging site iss marked by permanent survey markkings. The ccross-
section along
a this seection line is
i determineed by surveyying with thhe help of ssounding rods or
soundingg weights. When
W the deppth of water is
i more or iff quick and aaccurate deppth measurem
ments
are needeed, an echo sounder
s is ussed.
A discharge measuremen
m nt requires determination
d n of sufficieent point vellocities in a river
cross secction to perm
mit computattion of an average veloccity in the sttream. The sstream is divvided
into a nu
umber of verrtical section
ns. The sum of cross-secctional area oof each sectiion multiplieed by
the averaage velocity of each sectiion gives thee total dischaarge:
n m
Q total Q i (A i v i ) (4.3)
i 1 i 1
where Qtotal
t = total discharge,
d mean
m = number of sections, Ai = cross area of secttion I, vi = m
velocity of
o section i, Qi = discharrge in sectio
on i.
A current meeter measuress the velocitty at a point.. However, the mean veelocity in eaach of
the selectted vertical segment
s is required
r to estimate discharge. The mean velociity in a vertiical is
determinned from velo ocity observvations at on
ne or more p oints in thatt vertical. C
Current meterrs are
held dowwn and posittioned at thee required lo ocation in fl owing waterr by soundinng weights. The
weights are
a connecteed to the currrent meter by y a hanger aand pin assem
mbly.
The
T number of velocity determinatio ons is limiteed to those w
which can bbe made witthin a
reasonable time. If the river stage
s is chaanging rapiidly, one shhould quickkly completee the
measuremment with a minimum chhange in watter stage.
Figure
F 4.10: Procedure for
f a currentt meter meassurement
Usually there are national guidelines detailing how the vertical sections should be chosen
and in which depth the velocity measurements should be made to determine the mean velocity in
the vertical profile. These guidelines aim to ensure homogenous data sets. Two methods that are
frequently used to determine the mean velocity over the vertical profile are:
Two point measurement: The variation for most channels is such that the average of the
velocities at 0.2d and 0.8d below the surface equals the mean velocity in the vertical.
One point measurement: The velocity at 0.6d depth below the surface closely approximates the
mean in the vertical. The adequacy of these assumptions for a particular stream can be tested
by making a detailed vertical velocity measurement.
A measuring reach of the stream is chosen for which the mean cross section area and the
slope of water are determined. The mean velocity is computed by using the open channel flow
formula such as the Manning's formula. Appropriate value of the rugosity coefficient is chosen
depending on the physical conditions of the channel. Discharge is computed as the product of the
mean velocity and the average cross-sectional area of the reach.
For the measurement reach, three things should be known: (i) The cross-sectional
geometry and properties at its ends, (ii) the value of Manning's roughness coefficient n, and (iii)
water-surface elevations at the end sections. The river reach should be straight (to the extent
possible) having stable bed and banks and uniform cross-section over a length of at least five
times the width of the channel. In any case, the length should not be less than about 300m. As far
as possible, the river reach should be long enough so that the difference between water levels at
the upstream and downstream gauges is not less than ten times the uncertainty in the difference.
The flow in the reach should be free from significant disturbances, draw-down or back-water
effect of any structure or tributary joining upstream or downstream and it should be contained
within its banks for all stages at which this method is used. If no uniform reach is available, the
reach should preferably be converging rather than diverging. The measurement site should be
easily accessible at all times.
A minimum of three cross-sections are generally desirable in the selected reach (two
should be at the ends and others in-between, one at the centre). Water surface slope is computed
from the gauge observations at either end of the reach; intermediate gauge(s) are used to confirm
that the slope is uniform throughout the reach. Cross-section area and wetted perimeter are
estimated at each cross-section. Manning’s n should be estimated only by a trained person.
Measurement of slope
The upstream and downstream water level gauges shall be read simultaneously. If there is only
one observer and this is not possible, the upstream gauge shall be read first, followed by the
downstream gauge, followed by a second reading of the upstream gauge. If the value of the
upstream gauge has changed by more than 10 mm between the two sets of readings then the
average of the two upstream readings should be used. Effort should be made to gauge readings to
at least the nearest cm or better.
Slope of water surface is computed from the average of gauge observations at either ends
of the reach. The intermediate gauge(s) are used to confirm that the slope is uniform throughout
the reach. When accurate gauges do not exist or have been damaged, flood marks on the banks
may be used to estimate the slope. In such cases, an effort should be made to locate and use as
many flood marks as possible in the reach with least possible delay after the flood.
If intermediate gauges are installed, the readings at these gauges should also be taken and
the slopes between the upstream and intermediate gauges, and intermediate and downstream
gauges estimated and compared with the slope obtained using the upstream gauges. If there is a
significant difference in the three slope estimates (say greater than 10%) then the validity of the
estimates should be ascertained.
Since slope-area method is frequently used to estimate peak discharges following a flood
event, observation of high-water marks after the passage of a flood is important. This involves
identifying debris marks and then topographic level surveys to find their levels. Several types of
high-water marks may be found, such as drift on banks, mud lines on trees and banks, and drift
in bushes or trees. Each high water mark should be rated as excellent, good, fair or poor; this
could help when interpreting the high-water profile. Sufficient wrack marks need to be identified
on both banks to fully define the water profile over the measuring reach. A graphical profile of
the high-water marks should be made so that irregularities in the profile can be removed.
Sufficient number of points along the measuring reach need to be used to define the surface
water profile.
Estimation of Discharge
Fig. 4.11 shows the longitudinal section of a river between two sections, 1 and 2. The head at a
cross-section consists of water surface elevation and the velocity head. The head loss is made up
of two parts: frictional loss and energy loss due to expansion or contraction.
Fig. 4.11 Channel reeach for the slope-area
s method.
m
where L is the reach length, k is the coefficieent for energgy loss; its vvalue is 1 forr contractionns and
0.5 for ex
xpansions. Estimate
E the hydraulic raadius for botth the upstreeam downstrream sectionns (R1
& R2) using
u the esstimated crooss-sectionaal areas (A)) and wetteed perimeterrs (P): R= A/P.
Accordin ng to Manninng's formulaa, the mean velocity
v in reeach 1-2 is calculated as
v1-2 = (1/n)R2/3
2
S1/2 (4.5)
where R is the hydraaulic mean depth and n iss Manning's roughness ccoefficient. If A is the ccross-
section area, then thee discharge Q is
Q = (1/n)AR2/3
2
S1/2 = K S1/2 (4.6)
Using
U the esttimated hydrraulic radius, conveyancce (K) of thee channel at the upstream
m and
downstreeam sectionss is computeed. Conveyan nce dependss on channell characterisstics. As the flow
in the reach may not be truly uniform, average conveyance of the reach is computed as the
geometric mean of the conveyances of the two end sections (K1 and K2):
1
K1 A1 R12 / 3 (4.7)
n
1
K2 A2 R22 / 3 (4.8)
n
Using the values of K1, K2 and the surface water slope, the following equation is applied
to provide a first approximation of discharge:
Q K K S
1 2 f (4.9)
Using the value of Q, first estimates of the velocity heads at sections 1 & 2 are obtained
by:
2
V2 Q 1
VH1 1 1 (4.10)
2 g A1 2 g
2
V2 Q 1
VH 2 2 2 (4.11)
2 g A2 2 g
If VH1 < VH2, i.e., the velocity head difference is negative, then the reach is contracting
and the following equation can be used to obtain a second estimate of the slope of the energy
line:
S 2g 2g
(4.12)
L
If VH1 > VH2, i.e. the velocity head difference is positive, then the reach is expanding which
should be avoided. However, if an estimate is necessary, the following equation is used:
Using the new estimate of the energy line obtain a second approximation of discharge (Q2). If
(Q2 - Q1)/Q1 > ± 1%, repeat the iteration procedure until the current approximation of Q is within
1 % of the previous estimate.
The slope-area method can be used with some degree of accuracy in open channels with
stable boundaries, or in channels with relatively coarse bed material. This method may also be
used in other cases, such as alluvial channels including channels with over-bank flow or non-
uniform channel cross-sections, subject to the acceptance of large uncertainties involved in the
selection of the value of the rugosity coefficient, such as Manning's roughness coefficient n.
If the reach is substantially uniform and there are insignificant differences in the cross-
sectional areas, A1, A2 .... Am at the chosen sections, the mean area of cross- section for the reach
of the stream may be taken as
Mean wetted perimeter for the reach may also be estimated similarly
Table 4.1 Typical values of roughness coefficient for channels with relatively coarse bed
material
S. N. Type of bed material Size of bed material equivalent Rugosity coefficient 'n’
diameter in mm
1 Gravel 4 to 8, 8 to 20, 20 to 60 0.019 to 0.020, 0.020 to
0.022, 0.022 to 0.027
2 Cobbles and Shingle 60 to 110, 110 to 250 0.027 to 0.030, 0.030 to
0.035
The major limitation of the slope-area method is the difficulty of a correct assessment of
the rugosity (roughness) coefficient n in the Manning's formula. This method can be used with
some degree of accuracy in channels with stable bed and banks having relatively coarse bed
material. This method may also be used in other cases, such as alluvial channels, subject to the
acceptance of larger errors involved in the selection of the value of the rugosity coefficient n.
This method is not suitable for large channels, rivers with very flat slopes and high sediment
concentration, and channels with significant curvature.
The
T principaal differencee between a conventionnal measureement and the movingg-boat
measuremment is in thet method of computing water vvelocity. In the case of a conventtional
techniquee the mean velocity in the segmentts of a crosss-section of the stream is determineed by
point velocities meaasured at seeveral depth h is in the vertical. T
The moving-- boat technnique
measuress the velocityy over the width
w of a segment by suuspending thhe current m
meter at a connstant
depth durring the trav
verse of the boat
b across the
t stream. T The measureed velocity aand the depthh data
are used to determinee the discharrge.
The
T velocity observed at a typical ob bservation pooint in the crross-section (Fig. 4.13), vv, is
the veloccity of waterr past the currrent meter resulting
r froom both streaam flow andd boat movem ment.
It is the vector
v sum of
o the velocitty of water with
w respect to the stream m bed (v) annd the velocity of
the boat with respectt to the streaam bed (vb). The velocityy of water ccan be obtainned by meassuring
the anglee between n the path off the boat an nd a verticaal vane whicch aligns itself in a direection
parallel to the movemment of the water
w past it.
Referring
R to Fig.
F 4.13, Th ocity v, perpeendicular to the boat patth (true courrse) at
he flow velo
each obseervation poin
nt 1, 2, 3,…, can be deteermined from
m the relationship
v = vv sin (4.15)
here v is that compoonent of the stream velocity which iis perpendicuular to the ttrue course oof the
boat evenn though thee direction off flow may not
n be perpenndicular.
ation of Disccharge
Computa
After thee cross-sectiion has beenn selected, the
t width off the stream m is divided into an adequate
number ofo segments such that th here is not much
m variattion betweenn two adjaceent verticalss. Of
course, if
i previous measuremen nts have sh hown uniform mity of botth the crosss-section andd the
velocity distribution then fewer verticals may be takenn. It is betterr if no segmment carries more
than 5 to
o 10 percent of the total discharge. Fig.
F 4.14 shoows the crosss section off a river in wwhich
(n-1) verrticals are drawn.
d The velocity aveeraged overr the verticaal at each seection is knnown.
Considerring the totall area to be divided into
o (n-1) segm
ments, the tottal dischargee is calculateed by
the methood of mid-seection as:
n
Q =
i= 1
( vi a i ) (4.16)
where Q is the total discharge, ai is an individual partiaal cross-secttion area, annd vi is the m
mean
velocity in that areaa. The areaa extends laaterally from
m half the ddistance from m the preceeding
observatiion vertical to half the distance
d to the
t next andd vertically ffrom the watter surface tto the
sounded depth.
Path
h of the
boat
Flow
w vv
vb
Observatioon
verticals v
The sampling data recorded at each observation point provide the necessary information
to determine the velocity of the stream. There are two methods to obtain this velocity, referred to
as method 1 and method 2. Method 1 consists of measuring the angle a between the selected path
of the boat and a vertical vane which aligns itself in a direction parallel to the movement of the
water. An angle indicator attached to the vane assembly indicates angle a.
The reading from the rate indicator unit in pulses per second is used in conjunction with a
calibration table to obtain the vector magnitude vv. Normally data are collected at 30 to 40
observation points in the cross-section for each run. Where practicable, automatic and
simultaneous readings of all required parameters may be recorded
v = vv sin a (4.17)
The solution of equation (1) yields an answer which represents that component of the stream
velocity which is perpendicular to the true course even though the direction of the flow may not
be perpendicular.
Example 4.1: The data from a stream gauging station during a measurement campaign are given
below:
Solution: First, the average width and cross-section area of flow for each segment are computed.
The average width for the first segment (refer to Fig. 4.14), W1 , is
(1 2 / 2) 2
W1 2.0m
2 x1
(1 2 / 2) 2
W N 1 2.0m
2 x1
(2 2)
Wi 2.0m
2
ai d iWi
q i a i vi
4.7.1 Weirs
Weirs are used to control upstream water level (to divert water in another channel) or for
measuring discharge or for both. They produce a critical relationship between stage and
discharge by obstructing channel flow; the head over the crest of the weir can be related to
discharge. Weirs have a defined cross-section (either rectangular or triangular) and hence the
computation of discharge is simple. Flow conditions are more complex in flumes; a section of
critical, high-velocity flow is produced by a lateral smooth-walled constriction of the flow. In
both weirs and flumes, the influence of the water level downstream of the structure on flow
conditions at the upstream point of measurement is prevented. Figure 4.15 shows a typical weir
and flume.
Use of these structures in field conditions is limited by the ranges of head, debris or
sediment load of the stream, and the backwater effects produced by these. Since the size and cost
of the structures increase as the size (width and discharge) of the river increases, such structures
are usually built only in upper and middle reaches of rivers. Site requirements consist of a
reasonably straight (at least for a distance of five times the width) approach channel which
should be free of excessive sedimentation, weeds and other aquatic growth. The structure should
be rigid, water-tight, normal to the flow direction, and should be capable of passing high flows
without any damage to its body. The stage-discharge relation at the site depends on the
geometrical characteristics of such a structure.
Problems arise during the use of such structures on large rivers, particularly those that
carry appreciable quantity of sediments because the sediment may be deposited upstream of the
structure and change the discharge characteristics. Periodic inspections of the structure should be
carried out and any damage should be promptly attended.
Figure 4.15 Typical
T weir and flume
Rating
R of the weir and flumes is carried c out on the bassis of laboraatory calibraation.
Howeverr, the rating should be checked
c afteer constructioon. The hyddraulic princciple that govverns
the use of
o weirs or flumes is thatt these structtures produc e a unique ccontrol sectioon in the floww. At
these stru
uctures, the river
r discharrge is a functtion of the hhead of waterr:
Q = f (H) (4.188)
ge (m3/s), an
where Q is discharg nd H is the head of waater (m) at thhe structuree. For weirss, the
equation becomes:
Q = K Hn (4.19)
Q = Cd g b H1.5 (4.20)
where Cd is the discharge coefficient, g is the acceleration due to gravity, and b is the crest width
(m). The flow that is affected by tailwater conditions is known as drowned or submerged flow.
Discharge under drowned conditions is obtained by applying a reduction factor to the free flow
discharge.
4.7.1 Flumes
A flume is a flow measuring structure formed by a constriction in a channel. The constriction can
be either a narrowing section of the channel or a narrowing section in combination with a hump
in the invert. Critical conditions are produced providing sufficient amount of constriction. In
such cases, a unique stage-discharge relationship exists independent of the downstream
conditions.
Critical-depth flumes have an edge over a weir in situations when sediment material is
being transported by the flow. This sediment gets deposited upstream of the weir and affects the
discharge relation. The critical-depth flumes consisting only of horizontal contraction would
easily carry the material through the flume. Critical-depth flumes can be grouped into two main
categories, viz., long-throated flumes and short-throated flumes.
Constriction of long-throated flumes is sufficiently long (the length of the throat should
be at least twice the maximum head of water that will occur upstream of the flume) so that it
produces small curvatures in the water surface and the flow in the throat is virtually parallel to
the invert of the flume. This condition results in nearly hydrostatic pressure distribution at the
control section (where critical depth occurs) which, in turn, allows analytical derivation of the
stage-discharge relation. This gives the designer the freedom to vary the dimensions of the flume
to meet specific requirements. Such flumes are usually of rectangular, trapezoidal, triangular or
U-shaped cross-section. For a rectangular flume, the discharge of an ideal fluid is expressed as
3/ 2
2
Q g bH 3 / 2
3 (4.21)
here, H represents the upstream energy and b is the typical width dimension for the particular
cross-section shape of the flume. By introducing suitable coefficients, this equation can be
generalized to the following form
3/ 2
2
Q g Cv Cs Cd bh3 / 2
3 (4.22)
where Cv = coefficient to take in to account the velocity head in the approach channels, Cs =
coefficient to take account of the cross-section shape of the flume, Cd = coefficient for energy
loss, and h = depth of water, upstream of the flume, measured relative to the invert level of the
throat (i.e., gauged head).
4.7.2 Short-Throated Flumes
In these flumes, the curvature of the water surface is large and the flow in the throat is not
parallel to the invert of the flume. The principle of operation of these flumes is the same as that
of long-throated flumes, i.e., the critical condition is created at the throat. However, analytical
derivation of the discharge equation is not possible since large curvature of flow produces non-
hydrostatic pressure distribution. Further, energy loss cannot be assessed. Therefore, it is
necessary to rely on direct calibration either in the field or in the laboratory to establish and
verify the discharge equation. The designer has limited freedom in choosing the dimensions of
the flume and has to select one of the standard designs to meet his requirements. Such flumes,
however, require lesser length and, hence, are economical compared to long-throated flumes.
One of the most commonly used short-throated flumes is the Parshall flume.
Parshall of USA designed a short-throated flume with a depressed bottom known as the
Parshall flume after his name. Parshall flumes have gives satisfactory service at water treatment
plants and irrigation projects. These flumes consist of a short parallel throat preceded by a
uniformly converging section and followed by a uniformly expanding section. The floor is
horizontal in the converging section; it slopes downwards in the throat and in the expanding
section, it is inclined upwards. The control section at which the depth is critical occurs near the
downstream end of the contraction.
4.8 Advanced Discharge Measurement Techniques
Conventional methods of flow measurement by the use of open channel sections or calibrated
structures are not always feasible because of channel size, unstable stage-discharge relationships,
or other physical, hydraulic, or environmental considerations (such as aesthetic, navigation and
fisheries impacts) and costs. Three relatively new methods of flow measurement in open
channels that have gained importance in recent times are:
Dilution methods,
Electromagnetic method, and
Ultrasonic method
A tracer is an ion or compound which is introduced into the flow to follow its behavior. A
known quantity of a tracer is introduced in the flow at an upstream section. At a downstream
cross-section of the reach, the concentration of the tracer is measured at regular intervals of time.
The reach selected for measurement should be such that there is no loss or gain of water
and the reach length should be sufficient to achieve complete mixing. The tracer can be
introduced in two ways: (1) constant rate injection, and (2) gulp injection. In the constant rate
injection method, the duration of injection should be such that a steady regime of concentration
is achieved for an adequate duration (about 10 to 15 minutes) in the sampling section.
In the constant injection rate, a tracer of concentration c1 is injected at the upstream end
of the stream reach at rate q. No flow should join the stream in the reach. Water samples are
taken at a downstream point and after an equilibrium concentration c2 is reached, the discharge Q
can be calculated as:
c1 q + c0 Q = c2 (q +Q) (4.23)
c 2 c1
Q q (4.24)
c0 c2
Ideally the tracer used should not be absorbed by the sediments channel bank material
and vegetation. It should not chemically react with any of these, should not be lost by
evaporation and should be non-toxic. The tracer should not be very expensive and it should be
easy to detect it in small quantities.
Section 2
C2
C0 Background concentration
Time
Fig. 4.16 The concept of dilution method of discharge measurement with constant rate injection.
Figure 4.17 shows a sketch of an electromagnetic gauging station where the coil is
placed in the bed and the magnetic field is generated in the x direction, the EMF is in the y
direction, and the streamflow is in the z direction. Faraday’s law of electromagnetic induction
relates the length of a conductor moving in a magnetic field to the EMF generated. The
measures of induced probe voltage, coil current and depth are used to calculate flow as follows:
In practice, most river beds have significant electrical conductivity that will allow
electric currents to flow in the bed. From practical considerations, the induced field will be
spatially limited and electric currents flowing in the area outside the field will have the effect of
reducing the output potential. Both of the above factors have the effect of reducing the signal
and hence the voltage recorded. At an electromagnetic gauging station, it is necessary to
measure the conductivity of bed and water.
Figure 4.17
7 Basic systeem of the eleectromagnettic method (S
Source: ISO 9213).
Earth's
E magnnetic field, battery
b effeccts of the prrobes in the water, and external eleectric
noise cau uses problem
ms in measurrements. To overcome thhese, the maagnetic fieldd produced bby the
coil is rev
versed at reg
gular intervaals. Measurem ments are taaken with thee field switchhed first onee way
and then the other.
A limitation of
o the electrromagnetic method
m is thhat it has to be calibratedd by referennce to
current meter
m gauginng or some other
o form off flow measuurement. Whhen the gaugge is construucted,
its perforrmance need ds to be veriified by com mparison witth the resultss of an alterrnative method of
flow gau uging. The electromagne
e etic river floow gauge esstimates thee mean velocity in the eentire
cross-secction but due to the assumptions an nd limitationns, the instrruments needd calibrationn and
adjustmeent.
Special site
s related consideratio
c ons
The folloowing are the
t special consideratio ons when seelecting the site for ann electromaggnetic
gauge:
a) Thee channel att the measurring section should be sstraight for at least threee times the river
widdth at the measuring
m seection, it sho
ould be sym mmetrical about the riveer centre linee and
preeferably of (n
nearly) trapeezoidal or recctangular shhape;
b) Thee upper limmit of channeel width useed is approxximately 30 m for sitess with below w-bed
meembrane. Th his restriction is based on construcction cost; thhe cost of constructionn will
inccrease signifiicantly with the size of thhe river;
c) Theere should be no majo or source ofo electricall interferencce near the site. Electrrified
raillways, electrrical power grid
g lines an nd radio transsmitters are the most likely sources;
d) Prooximity to heavily
h reinfforced concrrete and sheeet piling maay distort thhe magnetic field
andd should be avoided
a if poossible;
e) It is important to ensure that the membrane is particularly well protected against damage
by vessels moving in the river;
f) The conductivity of the watercourse being measured needs to be uniform throughout.
An electromagnetic gauge will only measure flows passing through the vertical magnetic
field generated by the coil. The velocity of any water passing outside the measuring reach will
register as negative, leading to errors in overall flow calculation. If full range flow calculation is
desired, careful selection of the site is required to ensure it is not by-passed. Training
embankments may be needed to contain the full range of flows to be measured.
Electromagnetic gauges are more suited to rivers where the flow range is relatively small
e.g., high base flows and not very high flows in monsoon. These gauges are particularly suited to
rivers in which considerable weed growth can occur; the gauge will also tolerate high levels of
suspended solids. Since this is a non-intrusive method of flow measurement, it is thus
particularly suited to measuring flows in channels with considerable fish movement.
Electromagnetic gauge are not suited to rivers where base flows can be very low relative
to the size of the channel. The uncertainties can be large when the depth of flow is very low. This
is because the offset, which is introduced in the calibration to allow for the difference between
bed level and the position of the coil, becomes highly significant relative to the depth of flow.
Thus, this method may not work well for many large Indian rivers or non-perennial rivers.
The accuracy of the electromagnetic method depends on the signal processing equipment
detecting and measuring small potentials sensed at the voltage probes. The electromagnetic
gauging station requires on-site calibration by current meter or other means and a relation
established between discharge and output.
Weed growth in a watercourse can change the hydraulic flow conditions in the measuring
reach, and hence the stage-discharge relationship. Seasonal vegetation development may make it
necessary to correct computed discharge data to follow the variation in the stage-discharge
relationship. This can be achieved by carrying out sufficient measurements to compute the
deviation. In some cases, the weed growth patterns may change in such a regular and repeatable
pattern that a seasonal stage-discharge relation may be defined.
In the rivers that carry large volume of sediments, care has to be taken to prevent or
remedy significant siltation and choking of the inlet tube and/ or the stilling well itself. The risk
of siltation is very much dependent on local conditions. The stage-discharge relations of stations
operated in alluvial channels are subject to shifting-control conditions and typically more
frequent measurements are required to define a family of shift curves.
1. These regions mostly have inadequate infrastructure which makes movement and
communication difficult, time-consuming, unreliable, and expensive.
2. The harsh environment, dust, high soil erosion by strong winds, and the sediment carried
by rivers create problems for the conventional equipment. The end result is frequent
malfunctioning of equipment and loss of records.
3. The short duration and rapidity of onset of floods, coupled with (1), imply that the team
intending to measure such flows may miss them unless they are present at the site. The
duration of the peak flow may be so short that it is almost impossible to carry out
measurements.
4. Soft erodible beds make sounding operations difficult when scouring may occur beneath
the sounding weight as it touches the river bed.
5. The channel may shift laterally, isolating gauge and recorder well from the flow.
6. The bed level may vary due to scour and fill, particularly during the passage of a flood.
Such bed variations are difficult to measure and lead to errors in the application of cross-
sectional area rendering the resultant discharge value erroneous.
A detailed and careful planning of logistics, local knowledge, and experience are the key
factors in successfully tackling the river gauging problems of such regions. A combination of
conventional wisdom with new and appropriate technology by employing automation, telemetry,
and remote sensing is necessary to overcome such difficulties. Many equipment with minor
improvements could work well under arid and semi-arid conditions. Data logging systems, with
retrieval and telemetry have proved reliable under harsh environments.
A sizeable number of streamflow-gaging stations in the cold regions are affected by ice
during the winter. The procedure used by the U.S. Geological Survey for measuring discharge
under an ice cover involves averaging streamflow velocity at 0.2D and 0.8D, where D is the
effective depth of flow (measured vertically from the streambed to the underside of the ice
cover). Another approach is to use a coefficient to adjust point velocities at either 0.5D or 0.6D
to give the correct depth-averaged mean. A previous study to determine the streamflow
adjustment coefficients had found that there was considerable variability from station to station.
Hence, it is important that the coefficients for individual stations be determined.
On the basis of bias and root-mean-squared error, the first-visit complete-profile method
with velocity measured at 0.5D was recommended by Walker (1994) for implementation. In this
method, a complete set of vertical profiles for the first measurement of the winter season is used
to establish the adjustment coefficient for an individual site. Thereafter, the same adjustment
coefficient is used for subsequent measurements at that particular station and it was found to
result in nearly unbiased discharge measurements across the streamflow-gaging network.
4.10 Closure
The main objective of river gauging is to get as close as possible to the existing realities of the
field. The choice of measuring methods and the means of implementing them is a function of the
configuration of site, available resources, both human and material, and the degree of accuracy
expected. In difficult conditions, such as equipment breakdowns or dangerous situations,
estimations of speed, photographs taken, observing floodwater marks on permanent structures,
etc., are actions that can mitigate the absence of a true measurement and establish coherent limits
on further evaluation of a site. The selection and adaptation of gauging sites, competence of
personnel, equipment to be used for a job should all serve to advance a project toward optimum
measurement precision. The most important elements for good measurements are motivated,
well trained and properly equipped personnel. Not only the must the measuring equipment itself
be in good condition, but also the annex equipment must be in good working order.
REFERENCES
Boiten, W. (2008). Hydrometry, CRC Press, London.
Herschy, R.W. (2009). Streamflow Measurement, Routledge, London.
ISO 4369:1979. Measurement of liquid flow in open channels - Moving boat method.
International Organization for Standardization, Geneva.
ISO 6418:1985. Liquid flow measurement in open channels - Ultrasonic (Acoustic) velocity
meters. International Organization for Standardization, Geneva.
ISO 772:2011. Liquid flow measurement in open channels - Vocabulary and symbols.
International Organization for Standardization, Geneva.
ISO 748:2007. Liquid flow measurement in open channels - Velocity-area-method. International
Organization for Standardization, Geneva.
ISO 9213:2004. Measurement of total discharge in open channels -- Electromagnetic method
using a full-channel-width coil. International Organization for Standardization, Geneva.
Subramanya, K., Engineering Hydrology, Tata McGraw Hill Education Pvt. Ltd., New Delhi.
Walker, John F. (1994). Methods for Measuring Discharge under Ice Cover. Journal of
Hydraulic Engineering, ASCE, 120(11), 1327-1336.
WMO (2008). Guide to Hydrological Practices. Volume I: Hydrology – From Measurement to
Hydrological Information. WMO No. 168. World Meteorological Organization, Geneva.
Annexure A
Preliminary survey
A detailed examination of a large scale map of the area is required in the first instance; supported
by an aerial survey if necessary. This may be necessary if the area is not easily accessible by
other means. Aerial surveys and or satellite imagery can be used as a basis to select potential
sites. Enquiries should be made to determine if any plans exist for the modification of the river
reach which would modify the stream bed regime and impact the proposed gauging station.
Selection criteria
A list of potential sites shall be established with their advantages and disadvantages identified.
The establishment of the gauge zero shall be chosen so as to avoid negative readings. Thus it
shall be set well below the level of the control feature. This zero point shall be correlated with a
national datum through a station benchmark.
Stage recording
Customer requirements shall dictate the method of recording stage. It may be that a single record
of stage taken daily and read manually will suffice. It is more common to provide a continuous
record of stage utilising water level sensors, such as floats, pressure transducers, and echo
sounders, interfaced with a digital recorder (logger or telemetry) or analogue recorder (chart).
When a recorder is used, visits by the observer should be made from time to time to ensure
satisfactory performance of the sensor and recorder.
Discharge measurements
Discharge measurements using the above techniques shall be related to a stage reading taken at
the beginning and end of the discharge measurement and during the measurement if the stage is
changing rapidly or inconsistently. When sufficient numbers of discharge measurements have
been taken, a stage discharge relationship can be computed (see ISO 1100-2). Subsequent to the
formulation of this stage-discharge relationship, only occasional discharge measurements need
be taken at flows in the normal range to confirm the robustness of the relationship unless the site
is subject to shifting control conditions.
Discharge measurements made using the velocity-area methods can be performed using
rotating element current meters, electromagnetic current meters, acoustic Doppler velocimeters,
or acoustic Doppler current profilers. A variation of velocity area methods is the slope-area
method, which is typically used to compute flood discharge indirectly (see ISO 1070).
Module 5
Measurement and Processing of Meteorological Data
The meteorological factors controlling evaporation are air and water temperature, wind
speed, atmospheric pressure, incoming solar radiation, humidity, and saturation vapor pressure
deficit. In addition, water quality, depth, soil type and nature also influence evaporation.
Meteorological factors that affect transpiration are essentially the same as for evaporation. In
addition, vapour pressure gradient, available soil moisture and plant properties also affect
transpiration.
When installing evaporation pans it is important to ensure that the site of the pan is level
and free of obstruction. To the extent feasible, the ground cover should be maintained as close as
possible to the natural cover common in the area. Obstructions such as trees, buildings, shrubs or
instrument shelters should not be closer than four times the height of the object. The pan should
never be placed on a concrete slab, asphalt floor or gravel. The pan should be installed in such a
way that no object casts any shadow over the pan at any time. The recommended minimum size
of the plot where pan is installed should be 15 m x 20 m. The plot should be fenced to protect the
instruments and to prevent birds and animals from drinking the pan water. The fence should be
constructed so that it does not affect the wind structure over the pan. A barbed wire or iron mesh
fence is preferred. It may be necessary to prevent birds and small animals from drinking water by
using chemical repellants (these should not pollute water) in addition to a wire mesh. To estimate
the error introduced by the wire-mesh screen on the wind field and thermal characteristics of the
pan, readings from the protected pan should be compared with those of a standard pan at a
nearby site.
On days without rain, at daily (or twice-daily) reading time, water is poured into the pan
using a graduated cylinder to bring the level precisely to the top of the pointer gauge. The
volume of water added is recorded and represents the depth of evaporation. On the days when
rain has fallen since the last observation, the rainfall may exceed evaporation and water may
have to be removed from the pan to bring the level to the hook level. The adjacent rain gauge is
used to measure the rainfall. It there is a forecast of heavy rainfall on a given day, a measured
amount of water may be removed from the pan in advance so as to avoid pan overflow.
The water level in the pan must be measured accurately before and after water is added.
The water level may be determined by means of a hook gauge consisting of a movable scale and
vernier fitted with a hook, enclosed in a stilling well in the pan. An alternative arrangement is to
use a float. A calibrated container is used to add or remove water at each observation so as to
maintain the water level to a pre-specified depth.
where P is the depth of precipitation during the period between the two measurements, and d is
the depth of water added (+) to or removed (–) from the pan.
Pan evaporation provides an estimate of open water evaporation. Evaporation data from
pans are frequently used to estimate evaporation from water bodies, such as lakes and reservoirs
and evapotranspiration from an area. But a pan is a really small in size and its material is
different than a natural body. Hence its heat storage characteristics and air dynamics will be quite
different than a large water body. Therefore, evaporation from a pan will be higher than a large
open water body. An estimate of lake or reservoir evaporation (ER) is obtained by multiplying
the pan evaporation by a coefficient which is called the pan coefficient:
ER = KP Epan (5.2)
where KP is the pan coefficient; ER is the depth of evaporation from the reservoir and Epan is the
pan evaporation, both in mm/day. The value of pan coefficient depends on climate, geographical
location, season, size, and depth of the water body. This coefficient generally varies from 0.6 to
0.8. Lower values are typical of dry seasons and arid climates while higher values are
appropriate for humid climates. In the absence of better estimates, a value of 0.7 for the pan
coefficient may be used.
Some automatic evaporation pans are also in use. In these pans, water level is kept
constant by automatically releasing water into the pan from a storage tank or by removing water
from the pan in the case of precipitation. A record of the amount of water added to or removed
from the pan is maintained.
The maximum value of evaporation in the region is used as the first check to screen
doubtful measurements.
The observer may over-or under-fill the pan on a particular day. Such values may be
compensated for on the next day but the value for the current day will be wrong.
Among the instrument errors, leakage is the most serious problem and it occurs usually at
the joint between the pan base and side wall. Small leaks are difficult to detect but may
have a significant systematic effect on evaporation totals over extended period such as a
mouth. Where leakage has been detected, the pan evaporation records for a period of a
few months preceding the detection of leakage should be inspected carefully and marked
as suspect.
Animals and birds may be able to drink water from the pan especially if the covering wire
mesh is damaged.
Algae and dirt in the water will reduce the rate of evaporation.
Errors arise in periods of high rainfall when the depth caught by the rain gauge is
different from the rain fallen in the pan as a result of splash or wind eddies.
To detect error in pan evaporation data which may arise due to the above reasons,
warning and maximum limits on evaporation are fixed. Pan evaporation measurements are less
accurate during the periods of heavy rains but during these periods the rate of evaporation is also
less. Excessively high evaporation values during days should be flagged as doubtful and should
be subject to further validation.
A weighing type lysimeter uses mechanical balance to determine the change in water
content of the control volume. The hydraulics-based equipment employs hydrostatic principles of
weighing, and in the volumetric based ones, ET is measured by the amount of water added or
removed from the control volume to keep constant water content. Note that lysimeters are
difficult and costly to install and maintain. In view of the difficulties in direct measurement of
ET, indirect methods are generally used to estimate it. This requires measurements of
meteorological variables which influence evaporation. Commonly, pan evaporation is multiplied
by a coefficient to get the crop ET.
Evapotranspiration can be estimated by the water budget or heat-budget methods; many
empirical formulae have been developed which are based on meteorological data. Food and
Agricultural Organization (FAO) of the United Nations have adopted the Penman-Monteith
(PM) equation as the standard technique to compute reference ET (Allen et al. 1998).
Researchers are experimenting with satellite data (see, Bastiaanssen and others, 1998) to
estimate regional actual evapotranspiration. Remote-sensing of several important parameters
used to estimate ET is made by measuring the electromagnetic radiation in a particular waveband
reflected or emitted from the Earth’s surface. The soil moisture may be estimated using the
measurement of microwave properties of the soil (microwave emission and reflection or
backscatter from soil). However, there are uncertainties in such soil moisture estimates due to
previously mentioned factors such as surface roughness and vegetative cover.
Type of Purpose
thermometer
Dry bulb To measure ambient air temperature.
thermometer
Wet bulb To measure temperature that will be attained by a volume of air if cooled
thermometer adiabatically to saturation by evaporation of water into it; all latent heat
has been supplied by the volume of air. This temperature is used to
calculate relative humidity.
Maximum To indicate the highest temperature reached since the last setting.
thermometer
Minimum To indicate the lowest temperature reached since the last setting.
thermometer
Two units are commonly followed in measuring and reporting temperature: Centigrade
and Fahrenheit. In India, centigrade is used more commonly. In the case of the dry bulb, wet
bulb and maximum thermometers, observations are of the position of the end of the mercury
column but in the case of the minimum thermometer, the reading is taken of the position of the
end of th
he dumb-belll shaped inndex farthestt from the bbulb. After tthe maximuum and miniimum
thermom
meters have been
b read, theey are reset using
u a standdard proceduure.
The
T wet bulb b is so called
d because th he bulb of thhe thermomeeter is alwayys kept wet. Wet-
bulb temmperature is measured
m by
y using a theermometer w whose bulb is wrapped in cloth (callled a
sock). Th he sock is kept
k wet with h water via wicking or capillary acction. One ennd of the soock is
immersed d in a small cup of wateer and the otther surroun ds the bulb. It is to be eensured that there
is alwayss enough waater in the cu up so that thee sock is keppt wet. Wet-bbulb thermoometers shouuld be
shielded from radiattion as welll as adequatte ventilatioon must be maintained. An accuraccy of
0.3C is enough in most
m temperaature related applicationss.
Regarding
R th
he use of th hermometers, the naturee of intendeed use will determine w which
thermom meter is most suitable. Daata of the dry y bulb thermmometer is thhe most frequuently used input
in hydrologic studies and obserrvations of ambient
a air temperaturee are best caarried out w
with a
mercury--in-glass theermometer. Continuouss records oof temperatuure may bee obtained with
resistance or therm mocouple ellements. Fo or measurinng and reccording waater temperaature,
thermogrraphs, which h produce a continuou us record oof temperatuure are usedd. Usually these
comprisee of a mercu ury-in-steel sensing
s elem
ment immerssed in the waater, which is connectedd to a
recorder (WMO 200 08). Care sho ould be takeen in the insstallation of thermograpphs to ensuree that
measurem ments taken are represen ntative of th
he water tem mperature. Inn case temperature of thhe air
above a water body is needed, it i should bee measured aat 2 m abovve the waterr surface neaar the
centre off the water boody.
Temperature
T is also meaasured contiinuously usiing a thermoograph in w which changges in
temperatuure are recorrded throughh the use of a bi-metallicc strip. The ttemperature is registeredd on a
chart on a clock-driv ven revolvinng drum and d the measurrement (chaart) period m may be eitheer one
day or onne week. Th he observer extracts temmperatures at a selected interval froom the chartt. The
manuallyy observed reading on the dry bu ulb thermom meter is meeasured and recorded aat the
beginning and end of o the chart period and if these diff ffer from thee chart valuee, a correctiion is
applied to
o the chart readings at th
he selected in
nterval.
To measure water surface temperature, several types of thermometers, such as mercury-
in-glass or mercury-in-steel (including maximum and minimum and reversing thermometer),
platinum resistance or thermistor elements with electronic circuit and meter or recorder and
thermocouple thermometers, are used. In general, the precision required for the measurement of
water temperature is ±0.1°C. However, in many circumstances precision of observation of
±0.5°C is adequate.
Thermometer faults usually result in persistent systematic errors in data. Most of the
common errors can be detected by setting up appropriate maximum, minimum, and warning
limits for a particular station and these may vary seasonally. For example, summer maximum
temperature at a place can be expected not to exceed 50°C and the winter maximum temperature
not to exceed 35°C. Of course, these limits will change with location. Other possible checks
include:
Dry bulb temperature should be greater than the wet bulb temperature. On rare occasion
they may be equal.
Maximum temperature should be several degrees more than minimum temperature.
Maximum temperature measured using the maximum thermometer should be greater than
or equal to the maximum temperature recorded by the dry bulb during the interval. The
value of the maximum temperature will be set to the observed maximum on the dry bulb
if this is greater.
Minimum temperature measured using the minimum thermometer should be less than or
equal to the minimum temperature recorded by the dry bulb during the interval. The
value of the minimum will be set to the observed minimum on the dry bulb if this is
lower.
where the mixing ratio r = mv/md. Here, mv is the mass of water vapor and md is the mass of dry
air. The saturation mixing ratio rs is the ratio of the mass of water vapor in a volume of the
saturated air to the mass of dry air. The relative humidity (%) can be defined as the ratio of the
actual vapor pressure of the air to the saturation vapor pressure at the same pressure and
temperature. The relative humidity does not vary rapidly with time. Places close to sea have
higher RH (due to availability of moisture) and a smaller daily variation than inland locations.
RH is most commonly assessed by the joint measurement of dry bulb and wet bulb
temperatures (see section 5.2). From these measurements, the dew point temperature, and actual
and saturated vapor pressures may also be calculated. While the actual vapor pressure may vary
little during the day (except with the incursion of a new air mass), RH has a regular diurnal
pattern with a minimum normally coinciding with the highest temperature (when the saturation
vapor pressure is at its highest). It also shows a regular seasonal variation. RH is calculated from
the wet bulb depression (difference between dry and wet bulb readings) using a set of tables.
Measurement errors for the dry and wet bulb thermometers are the same as for
temperature. If the wet bulb is not properly maintained, it will show a high value of the
temperature from the wet bulb thermometer and consequently RH will be estimated higher.
Errors in the hygrograph may arise from poor calibration.
Wind
W speed is measured by using a cup typee anemomeeter (Fig. 5..3). Cups oof the
anemomeeter rotate due
d to the prressure exertted by windd movement.. The numbeer of rotatioons of
the anemmometer overr a time inteerval is displlayed by a ccounter and oor recorded by a data loogger.
This nummber indicatees the averag ge wind speeed over the time intervaal. Standard practice in India
is to meeasure the wind speed d over a th hree-minute period andd this is coonsidered ass the
instantan
neous wind speed at th hat time. Wind
W speed ffor the day is calculatted from coounter
readings on successive days. Observations
O are made daily in thee morning oor twice-daiily in
morning and evenin ng. Wind sp peed measurrements may ay be instanntaneous; theese may alsso be
accumulaative. By noting the cou unter reading
g at the beginnning and ennd of a periood, average wind
speed forr the time intterval can bee determined
d by dividingg this value bby the time iinterval.
The
T direction n of wind is reported
r as 16
1 points of the compasss; it could bee shown eithher as
numericaal or as alph
habetic code. Note that wind
w directiion is not a variable in formulas ussed to
compute evapotransp piration. Thee wind directtion may inffluence evapporation at a place if the wind
is coming
g from a place whose hu umidity is sig
gnificantly ddifferent.
Typical measureme
m ent errors
Errors in
n wind speed
d might arisee as the obseerver may wrrongly read and report tthe reading oof the
counter. Fault may develop
d in th
he instrumennt due to weear and tearr or poor maaintenance. W Wind
speed maay have largee variation with
w time and d distances aand hence itt is difficult tto set up lim
mits to
screen su
uspected valuues. Wind sp peeds in exccess of 100 kkm per hourr should be cchecked careefully
with otheer meteorolo
ogical variables.
The
T instrumeent commonlly employed d in India forr observatioon of the sunnshine duratiion is
the Camp pbell Stokess sunshine reecorder (Fig. 5.4). It is a glass spherre mounted on a sectionn of a
sphericall bowl. The sphere focu uses sun’s rays on a caard held in tthe grooves of the bow wl and
graduated d in hours. The sun’s rays burn th he card dayy when the sun is shinning. The caard is
changed daily after sunset. Note that, the sunnshine recordder uses the movement oof the sun innstead
of a clocck to form th
he time basiis of the reccord. Differeent grooves in the bowl must be ussed in
winter suummer and the
t equinoxees, taking diffferent card types. The llengths of buurnt traces oon the
sunshine card indicaate the sunsh hine duratioon. Sunshinee duration ddata at requiired resolutiion is
tabulated
d from the caard.
The
T use of th he sun as a timing
t devicce avoids tim ming errors. Errors mayy arise in thee data
from the wrong placcement of ch hart. The errrors may alsso creep in while the obbserver notees the
duration of sunshine from the chaart. The valuues below zeero or greateer than the m
maximum posssible
sunshine hours for th
he location shhould be lag
gged for furthher checkingg.
Typical measuremen
m nt errors
The sunsshine recordeer is simple instrument in
i principle aand the use oof the sun raather than a clock
as a timee base avoid
ds timing errrors. Potential errors maay arise from
m the use off the wrong chart
which may
m result in the bum reaaching the edgee of the chart, beyonnd which it is not registtered.
Mistake may be made by the observer while extraction of information from the chart.
5.6 Closure
This module discusses the methods of observation and processing of data pertaining to several
other variables which are used in hydrological studies. Although data of these variables are
important, unfortunately the network to measure these variables is very weak in our country.
Thus there is a need to strengthen the network to observe these variables and create a strong
database.
With a view to Standardize the observation of meteorological data, WMO (1988) have
issued a document listing recommended practices for meteorological observations.
References
Allen, R.G., Pereira, L.S., Raes, D., and Smith, M. (1998). Crop Evapotranspiration, Irrigation
and Drainage Paper No. 56, Food and Agriculture Organization, Rome, Italy.
Bastiaanssen, W.G.M. (1998). Remote sensing in water resources management: The state of the
art. International Water Management Institute, Colombo, Sri Lanka.
Engman, E.T. and R.J. Gurney (1991). Remote Sensing in Hydrology, Chapman and Hall, London.
Schultz, G.A. (1993). Application of GIS and remote sensing in hydrology. HydroGIS 93:
Application of Geographic Information Systems in Hydrology and Water Resources. IAHS
Pub. No. 211, 127-140.
WMO (1988). General Meteorological Standards and Recommended Practices. WMO - No. 49,
World Meteorological Organization, Geneva, Switzerland.
WMO (2008). Guide to Hydrological Practices. Volume I: Hydrology – From Measurement to
Hydrological Information. WMO No. 168. World Meteorological Organization, Geneva.
Module 6
Measurement and Processing of Water Quality Data
Pure water rarely occurs in nature. In fact, water commonly contains a variety of
constituents. Quality of water changes due to natural and human actions. Although rain water is
pure, it may interact with the matter present in atmosphere and become polluted. After falling on
the ground, water dissolves and carries impurities present on the ground. Infiltrating water comes
in contact with soil and rocks and some of these are dissolved and transported by water. Natural
waters often interact with organic materials. Human activities produce municipal, agricultural,
mining, and industrial waste. Activities such as recreation, social and religious rituals also cause
pollution.
By the time the precipitation water reaches the surface of the earth, it has already
collected a number of substances. Gases may have been absorbed or dissolved in it and smoke
and dust particles may have been picked up. If the atmospheric water has picked up radioactive
or acidic material, it may not be clean and suitable for some uses.
Precipitation water is subject to further changes of quality after hitting the earth surface
and when it travels on or below the ground, depending upon the place where precipitation falls.
Water reacts with carbon dioxide released by bacterial action on organic matter and acquires
acidic properties. In the acidic conditions, many soil and rock constituents are dissolved in water.
Anthropogenic actions influence the quality of water in significant (and mostly harmful)
ways. Dumping of municipal and industrial waste after little or no treatment in natural waters
results in drastic deterioration of the water quality of natural sources. Agricultural related
activities such as irrigation, use of fertilizer, pesticides, herbicides, etc., may lead to diffuse
pollution of both surface waters and ground water. Irrigation return flows contain chemicals and
also tend to increase salt contents in the receiving water. Construction schemes, such as those
connected with river training, flood control, low flow augmentation, etc., considerably influence
the quality regime. Unscientific mining activities often degrade water quality.
A big range of parameters are used to characterize waters quality. Largely the water
quality management objectives and the previous history of the water body will determine
selection of parameters. Some parameters, however, are of special importance and deserve more
attention. The commonly used parameters are discussed in the following sections. An integral
part of water resources management is that water quality standards are formulated, and
implemented, by doing measurements and monitoring.
Common measurements to estimate the presence of physical impurities and their quantity
in water are turbidity, solids, electrical conductance, color, odor, and temperature. Temperature
measurements help understand the value of density, viscosity, vapor pressure, oxygen saturation
value and rates of biochemical degradation. The test for residue is important in the management
of sewage treatment to indicate the physical state of the principal constituent. The solids present
in the dissolved form impact electrical conductivity of the water.
Monitoring objectives determine which variables are to be measured. Different variables
are measured at different frequencies which depend on the rate of changes in the variable as well
as its intended use. Typically, studies are carried out to collect and define base line or
background information, suitability of water for a designated use, or environmental impact
studies.
Generally water quality variables are not measured continuously because of logistical
reasons. Taking water samples, storing and transporting them to the laboratories, and testing
them are expensive. Therefore, sampling is done at regular (but large) time intervals and
sometimes at irregular intervals. Important water quality variables or those which change rapidly
with time are measured more frequently than those which change slowly. For very few water
quality variables, samples are analysed at the observation station. Most analyses require
sophisticated equipment and are, therefore, carried in a laboratory. A wide range of measurement
techniques have been developed for water quality variables. Here, the discussion is categorized
according to physical, chemical, and biological variables.
After the water quality data have been assembled in a data storage system, the next step is
to interpret the data with respect to specific questions and inputs in management decisions. The
most frequently asked questions are:
Solids can be classified by their state, chemical characteristics, and size and distribution.
Solids can be either dissolved (including colloidal and small suspended particles) or suspended
(including settleable). Particles passing a membrane filter with a pore size of about 1.2 micron
are considered dissolved and those retained on the filter are considered suspended.
The total solids content is the sum of dissolved and suspended solids. The amount of total
dissolved solids (TDS) in water is an important indicator of its quality for drinking, irrigation and
industrial use. TDS affects dissolved oxygen concentration and also influences the ability of a
water body to assimilate wastes. Dissolved solids affect ionic strength of water and thereby
impact mobility and transformation of metals. TDS also affects the growth and decay of aquatic
life.
To determine TDS, a sample of known volume is dried in an oven and the weight of the
residue divided by the volume of the sample, gives TDS which is normally expressed in mg/L.
Wide variations in the TDS concentration are observed in natural waters. TDS concentration in
rain water is below 10 mg/L; in river water it may be of the order of hundreds of mg/L.
Temperature of river water generally follows the long term average ambient temperature
of the area. Temperature of water in a lake varies with depth. Cold-water species reside in deep
waters while warm-water species are found in shallow regions near the edges.
Color: The color of water depends on factors such as the type of dissolved material and the
suspended particles. Colors associated with water are many times the result of colloidal
suspension, e.g., tea. Most colors in natural waters result from dissolved material coming from
decaying plants. Many industrial wastes are colored.
Odor: Pure water has no odor. Odor indicates some toxic pollution of water. Water may smell
due to the presence of decaying organic matter. The decaying organic matter may accumulate in
bottom where conditions are suitable for the anaerobic bacteria. Sources of the organics include
plant debris washed into streams, dead animals, microorganisms, and wastewater.
Hardness: Water is classified as soft when the carbonate concentration is below 50 mg/L and is
called as very hard when it exceeds 180 mg/L. According to the guidelines of the World Health
Organization, the hardness of water should not exceed 500 mg/L.
Silica (SiO2) is the main non-ionic mineral found in natural surface water. Its presence
creates problems especially in industrial applications, where it causes scaling. A variety of heavy
metals may also be found. Important among these are arsenic, barium, cadmium, chromium,
lead, mercury, selenium, silver, zinc, and cyanide. These constituents are of concern primarily
because of their toxicity to micro-organisms, plants, and animals. Typically these constituents
come from the discharge of improperly processed industrial wastes, and high concentrations are
often found in wastewater sludges.
Toxic metals and other elements may exist naturally in water. These metals include
Copper (Cu), Chromium (Cr), Mercury (Hg), Lead (Pb), Nickel (Ni), Cadmium (Cd), and
Arsenic (As). In many instances, concentrations of these have increased due to human activity.
Arsenic has been detected in many water sources in eastern India.
Bacteria are single cell protists. Most bacteria can be grouped into four general categories:
spheroid, rod, curved rod or spiral, and filamentous. Fungi are aerobic, multicellur,
nonphotosynthetic, heterophic, eucaryotic protists. Most fungi obtain food from the dead organic
matter. Fungi play an important role in the breakdown of organic materials in terrestrial and
aquatic environments. Consequent to decomposition of organic materials, carbon dioxide is
released to the atmosphere and nitrogen to the terrestrial environment.
6.3 Types of Pollutants
Pollutants can be categorized in five classes: Pathogenic micro-organisms, Organic matter,
Nutrients, Dissolved solids, and Trace compounds.
Organic matter depletes dissolved oxygen (DO) through microbial decomposition. Deficit
of oxygen or anaerobic conditions change ecological regime: fishes may die, water emits foul
odours and unsightly conditions develop. BOD test measures aggregate organic matter in terms
of oxygen equivalence.
Organic matter is classified as that biochemically decomposable and that which cannot be
decomposed. Organic matter is characterized by BOD, COD and volatile solids. BOD test
measures biochemically decomposable or degradable organic matter, the other two measure total
organics. It may be pointed out that the three parameters reflect the aggregate organic
constituents. The BOD test measures oxygen equivalence of organic matter. The equation for
BOD decomposition reaction may be written as:
Bacteria utilize various types of waste organic matter as their food to obtain energy. They
decompose the organic matter to carbon dioxide and water in presence of oxygen. The amount of
oxygen consumed is taken as a measure of the quantity of organic matter. The test is carried out
in a laboratory representation of the aquatic environment; a water sample is added to oxygenated
water and the loss of oxygen measured over a period of time.
Dissolved Oxygen: The DO amount in a surface water body is an important indicator of its
‘health’. The presence of oxygen in water is necessary to maintain favourable conditions for
growth and reproduction of a normal population of fish and other aquatic life. The absence of a
low level of DO in surface waters indicates pollution by organic matter. Under such conditions,
organic matter is decomposed by anaerobic bacteria and obnoxious end-products, such as
hydrogen sulphide, ammonia, etc. are produced
Generally, the dissolved oxygen concentration is near saturation in pristine rivers.
However, concentrations of gases change throughout the day. Usual concentration of CO 2 in
river water could be up to 10 ppm. BOD is generally 1 to 2 ppm for unpolluted river and. pH
values are generally around 7 (neutral). Turbidity is very high in flood season. Quality of ground
water is superior to that of surface water in some respects. For example ground water is free from
suspended solids and objectionable colour. But at the same time, ground water generally contains
higher dissolved solids.
Nutrients are those chemicals which contain elements, such as nitrogen (N), phosphorus
(P), carbon (C), sulphur (S), calcium (Ca), potassium (K), iron (Fe), manganese (Mn), boron (B),
and cobalt (Co). These elements are essential for the growth of living beings. Some of these are
required only in very small quantities and are called micro-nutrients. Plants require relatively
large amounts of C, N, and P. They obtain C from carbon dioxide and N and P from soil or
water. N and P may be limiting in the aquatic or land environment; that is the concentration of
one or other of these species may dictate the biomass of plant species which can survive in a
particular water body.
Nutrients such as nitrogen and phosphorus are essential elements for growth of plants and
animals. Both organic and inorganic forms of these are important. In water quality studies, the
nitrogen-containing compounds of interest are organic nitrogen, ammonia, nitrite and nitrate.
Phosphorus is of importance in aquatic environments. Phosphorous enters into water bodies
through fertilizer, industrial waste and rocks. Phosphorous is vital for all organisms and in many
cases it is the nutrient that limits productivity. Municipal and some industrial wastewaters
contain N and P. When such wastes are added to water bodies, it may result in accelerated
growth of algae (algal blooms or eutrophication). When the algal cells die, it increases the
decomposable organic matter in the water body resulting in consumption of oxygen and the
quality of water degrades.
Nitrogen accounts for about 80% of the gases present in the atmosphere and maintains
equilibrium of nitrogen concentration in open water bodies. Human activities influence nitrogen
in surface water in several ways: A) Nitrogen may be present in wastes discharged into surface
water bodies; B) Runoff from agriculture areas contains nitrogen; C) Nitrogen is also present in
the exhaust of automobiles and industries. The nitrogen solubility in water is very less; of the
order of 15 mg/L. Nitrogen is frequently present as nitrate and ammonia. The usual concentration
of nitrate in streams is 0.5 to 3 mg N/L and of ammonia about 3mg N/L. Runoff from
agricultural areas may have considerably higher concentration of nitrates. Domestic sewage
contains about 15-100 mg/L of the total nitrogen. Nitrification, which is oxidation of ammonia
and nitrite to nitrate, consumes dissolved oxygen in water. Excessive presence of nitrogen in
water can lead to eutrophication.
A large variety of substances exist in water in their ionic form. The major cations are
(positively charged ions) comprise calcium (Ca++), magnesium (Mg++), sodium (Na++) and
potassium (K+) and the associated anions include sulphate (SO4--), bicarbonate (HCO3-) and
chloride (Cl-). The cations having two positive charges are responsible for the hardness of water.
The aggregate salts are measured as total dissolved solids (TDS). Waters having less than 1500
mg/L TDS can be considered fresh waters.
A number of toxic metals and organic compounds may be added to water through
anthropogenic activity. Some metals which are toxic even in small concentrations are cadmium
(Cd), copper (Cu), chromium (Cr), mercury (Hg), lead (Pb), nickel (Ni) and zinc (Zn).
When water from a particular depth is to be collected an arrangement like Meyer’s sample bottle
may be used, Figure 6.2. After the bottle is lowered to the desired depth a slight tug removes the
stopper.
Alternatively, if the DO sampler is lowered quickly to the desired depth, it will collect sample
which would be representative of the water at that depth. Small amount of water which gets in
the bottle from other depths while the sampler is being lowered is likely to flow out when water
overflows into the sampler.
As an example, Table 6.1 gives DO values in BOD bottles, containing identical samples
when incubated at 20 and 27 ˚C over a period of 5 days. It also gives the cumulative oxygen
uptake values. Note that the BOD3 27 ˚C and BOD 5 20 ˚C values are nearly the same.
Table 6.1 DO values in BOD bottles, containing identical samples when incubated at 20 and 27
˚C over a period of 5 days.
Day DO, mg/L 20 Cumulative DO DO, mg/L 27 ˚C Cumulative DO
˚C consumed, mg/L, 20 ˚C consumed, mg/L, 27 ˚C
0 8.1 0 8.1 0
1 6.5 1.6 5.5 2.6
2 5.0 3.1 3.8 4.3
3 3.9 4.2 2.5 5.5
4 3.3 4.8 1.7 6.4
5 2.5 5.6 1.4 6.7
Figure 6.3 BOD curves for two samples: Sample 1 – easily degradable organics, Sample 2 –
difficult to degrade organics.
After a sample of water containing organic material is taken from a water body, the
oxidation reactions in the sample will continue thereby altering its BOD. To avoid this, the
samples should be analyzed as soon as possible after collection. If a delay in the analysis is
expected, samples should be stored around 4 to 5 ˚C which greatly reduces the rate of oxidation
reactions and the BOD does not change significantly.
6.5.3 Inferences Concerning BOD
The BOD test relies on growth of the micro-organisms degrading the organic matter in the
sample. Hence, right conditions to allow the microorganisms to grow without undue
environmental stress should be created. If the pH of the sample is too low (below 6.5) or too high
(above 7.5), the BOD analysis may be affected by the ability of bacteria to grow. Therefore, if
the pH of the samples is below 6.5 or above 7.5, it should be modified by the addition of acid or
alkali so that a pH is within the range 6.5 to 7.5.
Toxic material present in the BOD sample also inhibits the analysis and hence the sample
must be ‘seeded’ with bacteria so that the analysis proceeds normally. Organic compounds are
not the only materials which, when discharged to the aquatic environment, have an oxygen
demand. Ammonia, either free or when released from nitrogen containing organic compounds, is
also oxidised in rivers resulting in depletion of dissolved oxygen. The oxidation of nitrogen
compounds, carried out by nitrifying bacteria (a process known as nitrification) in the BOD
sample can be suppressed by the addition of an inhibitory chemical so that only the BOD
resulting from the oxidation of carbon compounds is determined. If such a chemical is not added,
the resulting BOD may be a combination of the oxygen demand caused by both carbonaceous
matter and ammonia based material in the sample.
Algae have the ability to produce oxygen by photosynthesis. Hence, if the sample
contains algae, they must be stored in the dark.
Example 6.1: Calculate the oxygen demand exerted by a sample of industrial waste whose
ultimate BOD is 256 mg/L, in 1, 2, 3 and 4 days. Assume that 25% of remaining demand is
exerted each day.
Solution: Table below shows the computations of BOD. First column shows the day number and
the second the BOD remaining in the sample. Column 3 gives the BOD exerted for the day and
the last column is the cumulative BOD. The BOD at the end of day t is the BOD at the beginning
of a day t less the BOD exerted on day t.
Time Remaining BOD mg/L BOD exerted for the day Cumulative BOD exerted
day mg/L mg/L
0 256 0 0
1 192 64 64
2 144 48 112
3 108 36 148
4 81 27 175
6.5.4 Factors Affecting the Oxygen Demand Rate
A number of factors control the rate of oxidization of organic matter in the aquatic environment.
These are the composition of the material, water temperature and the concentration of micro-
organisms. Generally the rate of the BOD reaction increases with temperature.
The concentration of micro-organisms present in the water body also affects the rate of
the BOD reaction. Normally there are sufficient numbers of suitable micro-organisms in a
natural water body to allow the BOD reaction to proceed. In some cases, for example, when
water contains chemicals which are toxic to bacteria, a few micro-organisms may be there to
carry out the oxygen demand reaction. If the BOD of such a sample were to be determined, it
would give a value much lower than that given the concentration of organic matter in the sample.
To get the correct results, the sample must be ‘seeded’ with suitable bacteria. Usually a small
amount of settled raw sewage is used as a source of seed. In case the waste is toxic, an
acclimated seed must be developed.
After a sample of water containing organic matter has been collected, it is possible that
natural biochemical reactions within the sample will act to reduce the COD. Therefore, samples
should be tested without delay. If delay in analysing the sample is unavoidable, the sample
should be preserved by acidification to pH = 2 or below.
Since the COD test measures both the biologically degradable and biologically non-
degradable organic matter, the COD value of a sample is always higher than its BOD value.
Further, COD test measures the total oxygen equivalents of the organic matter while the BOD
test measures oxygen equivalents of organic matter which is degraded in 3 days. This 3 day
value is about 70 % of the ultimate BOD. For domestic wastewaters and many decomposable
industrial wastes, COD is about 2.5 times the BOD. For unpolluted surface and ground waters
this ratio is even higher.
The Total Organic Carbon (TOC) content of a water sample is a measure of all the
carbon, which is bound up in the organic matter in the water sample. TOC is related to BOD and
COD as most of the oxygen demand measured during these analyses is due to organic carbon.
Since the BOD of a sample can be related to both the COD and the TOC, it is possible to
estimate the BOD from either the COD or the TOC. For this purpose, it is necessary to establish
a relationship between these variables for a particular sampling point. This is best done by
carrying out BOD and COD or TOC analyses on the sampling point under various conditions to
establish a reliable relationship.
Sediment being carried away by river flow originates from various sources, including the
watershed, agriculture field river bed and banks. The sediment data collected includes the
particle size distribution for the sediment and the sediment transport as bed and suspended load.
Sediment transport rates depend on many factors, and it may less than the sediment transport
capacity because of sediment availability. The total sediment transport by a river is the sum of
the suspended load and bed load. The sediment might originate as the bed material load and the
wash load. Fig. 6.4 shows the classification of sediment as suggested by ISO (ISO 4363: 1993).
Along Bed
Bed the bed load
material
Total
In load
Total load suspension Suspended (t
(origin)
load
Wash load in
suspension
Fig. 6.4 Definition of sediment load and transport (Source ISO 4363: 1993).
Sediment samplers are used to determine the concentration of sediment. The type of
sampler depends upon the need and the purpose. The suspended sediment samplers collect
samples of water and sediment mixture in a river. The sampler is lowered to the desired depth
and is filled with a sample. The sampler should be big enough such that an adequate amount of
sample, say at least 0.5 liter, is filled. A special class of samplers, known as the depth integrating
samplers, is used to obtain the vertical profile. The samplers are filled as they are lowered from
surface to bed and then raised up. The rate of lowering or raising should be constant in both
directions.
In flowing channel, samples are collected using scoop material sampler in shallow rivers
and low velocities. Dredge type bed material samplers are useful for higher depths and velocities.
In dry beds the samples are collected by scrappers after digging a pit of about 1 foot depth in
river bed. The samples collected are dried. Normally 3 to 5 samples are collected during pre-
monsoons monsoon and post-monsoon period. These samples are brought to laboratory for
analysis.
When doing the sampling, the width of the cross-section is divided into sub-areas,
depending on the variability of sediment concentration in the lateral direction. Specially designed
samplers are used to measure bed load sediments. Bed load samplers are lowered to the bottom
and are allowed to be filled up for 5-10 minutes. It is necessary to do sampling a wide range of
discharges so that a representative rating curve relating sediment and water discharge can be
established.
The concentration of sediments is commonly expressed in g/cm3 or kg/m3. The
international standards ISO 3716 and 4363 provide details of sediment sampling. The samples
are usually analyzed in laboratories. After allowing the sediments to settle down, water is
carefully removed from the container and the remaining sediments are oven-dried. The particle
size analysis is first carried out by sieving and then the finer sediments are analyzed using
hydrometer. The sediment transported as suspended load is classified in three categories,
depending on the particle size:
The coarse fraction (particles above 0.2 mm diameter),
The medium fraction (particles between 0.075 and 0.2 mm diameter), and
The fine fraction (particles below 0.075 mm diameter).
Acoustic Doppler Current Profilers (ADCPs) are being increasingly used for streamflow
and suspended sediment measurements. These instruments measure either the attenuation of an
acoustic pulse due to suspended particles or the backscatter of the pulse by the particles. Remote
sensing data are also now being used to estimate suspended sediments in water. The reflectance
properties of water change depend on the concentration and the properties of suspended
sediments and these form the basis to estimate suspended sediments.
where C is the suspended sediment concentration (mg/l), Q is the discharge (m3/s), and a and b
are constants. A typical sediment rating curve is shown in Fig. 6.5. Conventionally, discharge
and sediment concentration are plotted on a log-log graph paper and a straight line is drawn.
Usually, the power equation is log transformed, and linear regression is applied to estimate the
parameters. A least squares method can be used to obtain the values of parameters a and b.
Typically, exponent b in eq. (6.2) lies in the range between 2 and 3.
In a plot between, discharge and sediment concentration there will be a large scatter in
points. One reason behind this scatter is that soil erosion rates in a watershed are not the same
during different seasons of the year. If the scatter is large, it might be necessary to develop
separate rating curves for different seasons or according to streamflow generation mechanisms,
such as rainfall, snowmelt, etc. The data pertaining to rising and falling limbs of the hydrograph
may also be separated to improve the relationship.
A sediment rating curve is mainly applied to obtain the value of sediment concentration
for a given discharge. Along with the flow duration curve at a given location, the sediment rating
curve can also be used to estimate the amount of sediment transport over a period of time, say a
year. Another important use of sediment rating curve is in estimation of the impact of land use
changes and watershed management on sediment yield.
y = 0.1236x - 264.83
R² = 0.5833
Sediment Conc.
Discharge
Table 6.2 Primary water quality criteria for various uses of fresh water
Designated best use Class Criteria
Drinking water source A 1. Total coliform organisms MPN/100mL shall be 50 or less.
without conventional 2. pH between 6.5 and 8.5
treatment but after 3. Dissolved oxygen 6 mg/L or more
disinfection 4. Biochemical oxygen demand 2 mg/L or less
Outdoor bathing B 1. Total coliform organisms MPN/100mL shall be 500 or less
(organised) 2. pH between 6.5 and 8.5
3. Dissolved oxygen 5 mg/L or more
4. Biochemical oxygen demand 3 mg/L or less
Drinking water source C 1. Total coliform organisms MPN/ 100mL shall be 5000 or
with conventional less
treatment followed by 2. pH between 6 and 9
disinfection 3. Dissolved oxygen 4 mg/L or more
4. Biochemical oxygen demand 3 mg/L or less
Propagation of wild life, D 1. pH between 6.5 and 8.5
fisheries 2. Dissolved oxygen 4 mg/L or more
3. Free ammonia (as N) 1.2 mg/L or less
Irrigation, industrial E 1. pH between 6.0 and 8.5
cooling, controlled waste 2. Electrical conductivity less than 2250 micro mhos/cm
disposal 3. Sodium absorption ratio less than 26
4. Boron less than 2mg/L
Source: Central Pollution Control Board, Govt. of India.
6.9 Closure
Monitoring, modeling, and management of water quality of various sources of water has become
very important for India because many of our rivers and lakes are highly polluted. In fact, some
rivers are so polluted that they resemble a dirty water drain. A disturbing recent trend is that the
aquifers are also being polluted by seepage of polluted water from surface and rivers; some
industries are pumping polluted water in the aquifers to avoid treatment and escape detection by
enforcement agencies. Urgent actions are needed to restore water quality of the sources so that
the ecosystems remain in healthy conditions and provide sustainable benefits to the society.
6.10 REFERENCES
McCutcheon, S.C., Martin, J.L, and Barnwell, T.O. (1993). Water Quality. In Handbook of
Hydrology, Edited by D.R. Maidment. McGraw-Hill Inc., New York.
Reichel, G. (1998). Suspended sediment monitoring: Use of Acoustic Doppler Current Profiler.
In Encyclopedia of Hydrology and Water Resources, edited by Herchy, R.W., and Fairbridge,
R.W. Kluwer Academic Publishers, Dordrecht.
Swamee, Prabhata K. and Aditya Tyagi (2000). Describing water quality with aggregate index.
Journal of Environmental Engineering, ASCE, 126(5), 451-455.
WMO(1988). Manual on water quality monitoring -- planning and implementation of sampling
and field testing. Operational Hydrology Report No. 27, WMO No. 680. World
Meteorological Organization, Geneva.
WMO(1989). Management of ground water observation programmes. Operational Hydrology
Report No. 31, WMO No. 705. World Meteorological Organization, Geneva.
WMO (2008). Guide to Hydrological Practices. Volume I: Hydrology – From Measurement to
Hydrological Information. WMO No. 168. World Meteorological Organization, Geneva.
Module 7
Acquisition and Processing of Ground Water and Other Data
Keywords: Ground Water, Data, Time Invariant, Water Level, Pumping Tests, Validation,
Reservoir, Lake, Demographic, Socio-economic, Agriculture, Land Use, Water Use, Demand.
Geological material which can hold water consists of unconsolidated deposits or consolidated
rock. Water exists in the openings in this material. Water will flow depending upon how the
voids are interconnected. Void space may be saturated (totally filled with water) or unsaturated
(not completely filled with water). A geological formation that contains and yields sufficient
quantity of water is known as aquifer. An aquiclude is a formation through which inappreciable
quantities of water can flow.
An important hydraulic property that controls the movement of groundwater is the
hydraulic conductivity. Hydraulic conductivity is defined as the volume of water that will move
in a unit time under a unit hydraulic gradient through a unit area, which results in units of
velocity (distance per unit time). A related term, transmissivity, is defined as the hydraulic
conductivity multiplied by the aquifer thickness. It is important to note that the hydraulic
conductivity is a unit property, whereas transmissivity pertains to the entire aquifer. To
understand the ground water conditions in an area, it is necessary to have correct
hydrogeological information. The greater the volume and accuracy of this information, the better
will be the understanding of the aquifer response, water level changes, hydraulic gradients, etc.
Fence diagrams are used to represent stratigraphic data in three dimensions. They are
similar to cross sections, but rather than interpolating subsurface geology from a map, the
geology between stratigraphic sections or cores drilled into the subsurface is interpolated. Fence
diagrams (see Fig. 7.1) are very effective at demonstrating changes in geological formations and
stratigraphic relationships in a region. To construct a fence diagram, the location of each section
is marked on a paper in a same way as if it was a map, and a vertical scale chosen. A vertical line
representing the length of the section is drawn and the stratigraphic boundaries are marked along
the line. The next step is to choose pairs of sections between which to draw the “fence” based on
relative locations of sections and the lithologic and stratigraphic variations. Generally, those
sections are selected which will present the panel in the most useful orientation and will show the
widest variation in lithologic and stratigraphic relationships. Most sections are connected to two
other sections with panels; some of them may be connected to three and the sections on the edges
may be connected
c to
o only one seection. In so
ome cases, oone of the paanels will bee partially hiidden
behind annother one. Once
O all of the
t useful paanels are commpleted, the fence diagrram will show w the
three dim
mensional geometry of thhe various stratigraphic uunits.
7.1.2 Time
T Variannt Ground Water
W Data
The mosst common time variantt ground waater data arre the time series of w water levels in an
observatiion well and d withdrawaals from a well.
w An exiisting well m may be useed to observve the
water tab ble if the depth of th he well exteends sufficiiently below w the rangee of water level
fluctuatio
ons. Existing pumped wellsw can also
a be usedd for this purpose if thhe annular sspace
between the outer casing
c of th
he well and the pump column alloows smoothh passage oof the
measurinng tape or caable for meaasuring the water
w level. IIf an existinng in use weell is taken uup for
observatiions, measurrement shou uld be carried out after tthe cessationn of pumpinng and passaage of
sufficientt time so thaat the water level in thee well fully rrecovers. Liikewise, withhdrawal of w water
from any y other welll in the vicinity of th he observatiion well should also bbe stopped for a
sufficienttly long timee.
A time seriees of groun nd water daata typicallyy consists oof depth off water table or
piezomettric head meeasured at an
n observation
n well or pieezometer at ddifferent tim
mes. Time intterval
between the successsive data maay be uniforrm or non-uuniform. Annalyses of suuch a time sseries
provides valuable information on trend of the data, ground water potential and its use, and
estimates of aquifer parameters. Joint analysis of time series of water level and rainfall can
provide useful insight about recharge to ground water.
The water level in a well is most often observed manually by dipping a graduated tape
thereby producing a time series (at non-uniform) time intervals. In recent past, automatic water
level recorders are increasingly being used. In India, there are wide variations in frequency of
observations and ranging from three times a year to weekly observations.
Water withdrawal from the pumping wells are estimated by gauges, based on duration of
pumping, or from the quantity of electricity consumed. Of course, specification of the pump and
the pumping head should be known. The discharge of springs can be easily measured by
constructing a weir or V-notch.
The most common and simple manual method is by suspending a weighted line (for
example, a graduated flexible steel or plastic-coated tape) in the well from a fixed point at the
well head. The depth of the ground water level is obtained by subtracting from the total length of
the tape suspended in the well that length which has been below water or wetted. This wetted
part can be identified, for example, by putting chalk on the lower part of the tape before
measurement. Usually approximate depth to water level is known before hand and this serves a
check on the measurement. As depth-to-water level increases and consequently the length of tape
to be used increases, this method becomes cumbersome. Depths-to-water surface of up to 50 m
can be measured without much difficulty. At these greater depths, steel tapes of narrower widths
or lightweight plastic-coated tapes can be used. Depths to water level can be measured with
accuracy of a few mm by most methods and the accuracy usually depends on the depth.
A useful approach to organize and present water level measurements from a network of
observation wells is to produce an accurate map of well locations and then draw contours of
water-level at each well. Two types of maps can be prepared; depth-to-water table measured in a
well from the land surface or the elevation of the water level in the wells relative to an
established datum, such as the mean sea level. Generally, these maps are produced for a single
aquifer.
Depth-to-water Maps
Another type of maps depicts the depth-to-water level in a well relative to land surface. These
maps are referred to as a depth-to water map. These types of maps provide an indication of the
depth to drill to draw water which can be useful in planning utilization of resources. A map
based on the difference in depth to water between two different times would show, for example,
the areal variation of seasonal fluctuations. Of course such maps cannot be used to establish the
direction of groundwater flow because it also depends on the variation of topographic elevation.
Water Table Maps
A map drawn by using the elevation of the water table in a well with reference to a common
datum, such as the mean sea level, is known as a the water table map. Accurate elevation of
measuring point at each observation well is required to produce such maps and this is a difficult
task. Each depth-to-water measurement collected must be subtracted from the elevation of the
measuring point relative to the datum to produce the necessary data. This type of map has a
significant benefit that it can be used to infer the gradient of groundwater flow in many cases.
The accuracy of the map is depends on the accuracy with which the elevations of measuring
point are determined. Most commonly, the elevations are determined from the Survey of India
topographic maps or by carrying out leveling from a nearby benchmark. Alternatively, a GPS
may be used but the elevation data may have considerable errors.
7.3 Well Discharge Measurements
Measurement of a pumping well’s discharge is important in ground water modeling,
management and regulation. The discharge of a pumping well will vary with rise and fall in
groundwater level. This requires repeated measurements to keep track of the rate of change of
discharge. When a pump is turned on, the water level in the well drops quickly, causing the
discharge to vary. Pumping rate becomes stable within a few hours. Water-level changes that
could affect pumping rate can also occur as a result of recharge from precipitation or changes in
pumping of nearby wells.
The simplest method to determine the discharge from a pumping well is by measuring the
time required by discharge to fill a known volume. Dividing that volume by the time yields the
pumping rate. The accuracy of the measurement depends on the accuracy of measuring time.
However, at relatively high discharge rates, a measurement of this type may require logistical
planning and creating a large enough temporary storage space. Another way is to construct a V-
notch in the channel which carries the pumped out water.
In an aquifer test, pumping is controlled and generally held constant for the test duration
and the ground water levels in the pumped well and a nearby observation wells are measured. A
cone of depression is formed in the aquifer due to pumping. Depending on the time required to
achieve steady state conditions, pumping is carried out for an extended period of time.
Spring discharge may vary from zero or a trickle to a large flow which varies over time.
In a monitoring programme, it is essential to prepare an inventory of major and medium springs
and measure their discharges at least once in a month. In the case of small or medium size
springs, a channel is made and a V-notch is constructed for flow measurements.
It is a good practice to record the results of the validation in the database. For this
purpose, a special field may be created and a code, a flag, or a label is inserted after validation
task is completed and the data are considered to be correct. General concepts of data validation
have been discussed in Module 2. This section discusses some aspects specific to ground water
data.
For observation wells, the location and the elevation of the measuring point should be
correct because the measured depth to the groundwater level depends on these. Any shift in the
measuring point should be promptly corrected into the database. The depth of the well should
also be entered in the database and is important in data validation.
After entry, the data entered must be carefully compared with the field records to ensure
that there are no errors. A common method to check the data entry errors is to get the data
entered by two different persons and then compare the data entered by them. The values which
are same are considered to be correct and where there are differences, the corresponding data are
checked form the original records and corrected. It is assumed that the probability of two
different persons making the same mistake is nil.
The geographical locations of the wells or other structures may be checked by plotting
the locations on a map of the area. It is important that the attributes of the map layers are based
on the same master information as the database. Otherwise, correction of the data may give
unwanted results. For example, in the case of administrative names, frequently problems arise
with spellings. Different persons use different spellings and this causes difficulties in search.
Hence, it is important to use the same spelling in all the places.
Before validation of ground water level data, the well data should be validated including
checks on the details of wells. Some of these checks included in primary validations the well
location on a map, the elevation of the measuring point, etc.
Errors in ground water level data may be detected by carrying out simple checks which
flags:
• all groundwater levels which are below the bottom of the well,
• all groundwater levels which are above the top of the well plus 0.25 meter.
Visual inspection of plots of ground water levels usually identifies suspect groundwater
levels by comparison with the levels measured in a nearby well at the same time. An example is
given in Figure 7.2, where two ‘doubtful’ measurements are indicated in the well in upper
panel.
The correction of flagged doubtful measurements is not straightforward. A typical
problem is: What was the correct value? The observed value should only be corrected in case the
‘right’ value can be determined without any doubt. In all other cases where it is not possible to
alive at a ‘correct’ value, the value should be flagged as ‘doubtful’ value.
Water Level (m)
Well No 1
Water Level (m)
Well No 2
The anion and cation are plotted on their individual triangles and then their lines are
extended up onto the diamond diagram. Where the two lines intersect, the point for that water
sample is plotted. This point tells about the water sample in terms of its chemistry and the
controlling factors. If multiple data points from different sites form a line, it can be concluded
that all the water samples are derived from the same source. Ca2+, and Mg2+ is first plotted on the
cation triangle. The relative abundance of Cl-, SO42- and HCO3-+ CO32- is then plotted on the
anion triangle. The two data points on the cation and anion triangles are then combined into the
quadrilateral field that shows the overall chemical property of the water sample (Figure 7.3). The
Piper diaagram not on nly shows graphically
g th
he nature off a given waater sample, but also dicctates
the relationship to otther samples. For examp ple, by classiifying samplles on the Piiper diagram
m, one
can idenntify geologiic units with h chemically y similar wwater, and deefine the evvolution in w water
chemistry y along the flow
f path (Figure 7.3).
Diagrams
D succh as the Pipper diagram will bring oout deviatingg values. In a Piper diaggram,
outlying points indiicate extrem me water qu uality condiitions. However such rresults may also
indicate erroneous
e vaalues for the constituting
g componentts, e.g., the cconcentrationn of a compoonent
may havee been determmined by clo osing the ion
n balance wiith a residuall value.
7.6 Reservoir
R an
nd Lake Datta
Data perrtaining to reservoirs have both time invarriant/slowly--variant andd time invaariant
componeents. The tim me invariant/ slowly varyying data inccludes the eleevation-areaa-storage cappacity
table, an
nd operating g policies. Various
V colu
umns of thee reservoir working tabble, viz., innflow,
outflow through can nals, power house and spillway, eevaporation release to m meet the vaarious
demands, hydropow wer generated d (if any) etc.
e are the time-seriess data. Depending uponn the
purpose ofo the reserv
voir, these daata may be observed
o andd stored at m
monthly, dailly, or hourlyy time
interval.
7.7 Demographi
D ic and Socioo-economic Data
D
The dataa about popu
ulation and economic
e acttivities are aalso needed for planningg water resoources
projects. Census or demographic data are periodically collected in most countries. In India, an
elaborate exercise is carried out every 10 years to gather extensive demographic data which
forms an important input in national planning, including the water resources. The census data
also form the basis to forecast the population growth for use in planning activities. Census 2011
has been completed recently and a large volume of high quality demographic data has been
collected. In this survey, some new information was also gathered and preliminary results have
also become available.
It is essential to know the human resources available, since all development proceeds
through humans, for humans and on behalf of these. The survey should determine the size of the
active and potential labor force for agriculture, together with that for the tertiary sector
contributing to agricultural development. The population features as well as future projections
are needed to estimate the present consumption pattern and future demand of agricultural
products. The level of employment and the income of the entire agricultural and non-agricultural
active population must be found. It is a settled fact that assured irrigation is key in growth of
both mechanized large farms as well as manual labor-oriented small farms.
Normally, the socioeconomic and agriculture data are compiled for two conditions: ‘with’
and ‘without’ project. It is improper to base the analysis on a comparison of conditions before
and after project construction because a decline or improvement of present conditions might
occur even in the absence of the proposed project. This factor should be recognized while
determining the impacts attributable to the project. Consistent assumptions should be applied to
future conditions with and without the project so that comparability is assured. The net
incremental benefit stream should be an accurate reflection of the project's income generating
capacity, or its net contribution to real national income.
On the basis of the nine-fold classification, it is possible to arrive at the area according to
concepts like 'arable land' or 'potential land' available for crop. The 'arable land' would comprise
the 'net area sown' plus the 'current fallows' & 'other fallow land'. Similarly, the 'potential land'
available for cultivation would include, besides the 'arable land', the land under 'culturable
waste', 'permanent pastures & grazing land', & 'miscellaneous tree crops and groves, not included
in the net area sown'. However, such potential land cannot be ascertained without surveys of the
soil types & land-use & the economics of bringing them under cultivation.
The utility of the concepts of 'arable' and 'potential' land can be enhanced if some
indicators of arable and potential land are computed:
Percentage Potential land exploited = Net area sown/Potential land*100
Percentage arable land exploited = Net area sown/arable land*100
To standardize the concepts and definitions, standard definitions of the various categories
of land use have been adopted in India and these are given in Table 7.2.
Identification of areas earmarked for annual crops, orchards, pastures, forests and non-
cultivated along with the knowledge of the system of land tenure (owner occupancy, tenancy,
share cropping), size of farms, and fragmentation, enables preparation of preliminary blueprints
on the development possibilities. These details would also show the constraints on the
availability of land for cultivation. Soil properties would determine the spectrum of crops which
can be grown in the area.
It is also essential to know the present cropping pattern, crop rotations, varieties grown,
yields for each crop, and market price of each crop. These along with soil data are used to
optimize the cropping pattern and estimate irrigation water requirement for a proposed project.
Table 7.2 Standard definition of various categories of land use adopted in India
SN Category Description
1 Forest Includes all lands classed as forests under any legal
enactment dealing with forests or administered as forests.
2 Area under Non- Includes all lands occupied by buildings, roads and railways
agricultural Uses or under water, e.g. river, and canals and other lands used for
non-agriculture purpose.
3 Barren and un-cultivable Includes all barren and un-cultivable land like mountains,
land desert etc.
4 Permanent pastures and Includes all grazing lands where they are permanent pastures
other grazing lands and meadows or not. Village common grazing land is
included under this head.
5 Land under miscellaneous This includes all cultivable land, which is not included in
tree crops and groves etc. ‘Net Area Sown’ but is put to some agricultural uses. Lands
under Casuarina trees, thatching grasses, bamboo bushes, and
other groves for fuel, etc which are not included under
‘Orchards’ are classified under this category.
6 Culturable Wasteland This includes lands available for cultivation. Such lands may
be either fallow or covered with shrubs or jungles, which are
not put to any use. Land once cultivated but not cultivated for
five years in succession should be include in this category at
the end of the five years.
7 Fallow lands other than This includes all lands, which were taken up for cultivation
current fallows but are temporarily out of cultivation for a period of not less
than one year and not more than five years.
8 Current Fallows This represents cropped area, which are kept fallow during
the current year. For example, if any seeding area is not
cropped in the same year again, it may be treated as current
fallows.
9 Net Area Sown This represents the total area sown with crops and orchards.
Area sown more than once in the same year is counted only
once.
Source: Ministry of Statistics and Programme Implementation, Govt. of India.
The data about municipal can be obtained from municipal water works department. To
compute featuring water demand, the population is projected for some future year (say 2050).
Per capital water norms for urban and rural areas are used to compute the total water demand.
Data about industrial water demand are difficult to obtain. In absence of any better estimate, it is
assumed that the industrial water demand is the same as the municipal water demand. The
procedures to estimate of water requirements for agriculture, hydropower, etc. and their future
projections have been described in many text books.
7.11 Closure
Monitoring of ground water regime involves obtaining information on ground water levels and
quality through sampling. This chapter has also discussed some other data which are of immense
value in water resources planning and management. These include reservoir and lake data,
agriculture Data and land use data. Besides hydrology, data from social sector is also needed in
planning and management tasks. To that end, discussion includes demographic and socio-
economic data and water use and demand data.
7.12 REFERENCES
Jain, S.K., and Singh, V.P. (2003). Water Resources Systems Planning and Management.
Elsevier, Amsterdam.
WMO (1989). Management of ground water observation programmes. Operational Hydrology
Report No. 31, WMO No. 705. World Meteorological Organization, Geneva.
WMO (2008). Guide to Hydrological Practices. Volume I: Hydrology – From Measurement to
Hydrological Information. WMO No. 168. World Meteorological Organization, Geneva.
Module 8
ACQUISITION AND MANAGEMENT OF SPATIAL DATA
Keywords: Spatial Data, Point, Linear, Topography, Remote Sensing, Geographical Information
Systems
Large infrastructure is needed to capture, manage, and analyze spatial data. Spatial data
infrastructure consists of technology, hardware and software, policies and procedure/rules, and
people that are necessary to acquire, store and manage, analyze, and share the spatial data among
the users. It provides a structure of practices and relationships among data producers and users
that facilitates data sharing and use. A rigorous national data framework facilitates information
exchange and significantly reduces duplication of effort.
In GIS, maps are called spatial data. Information shown on paper maps can be input is
GIS as spatial data. Example of spatial data are stream network, well locations, villages, cities,
topographic contours, spot elevations, roads, land use, soil, geology, hydrological investigation
locations, hydrological response units etc. Spatial data are classified in to three types, namely
area, line and point. Areas are spatial data that are represented as closed figures, e.g. forests,
lakes, Thiessen polygons etc. Lines are spatial data that are represented as curves e.g. forest
boundaries, lake boundaries, contours, stream network, roads etc. Point spatial data are
represented as points on maps e.g. well locations, rain gauge stations, villages, etc. The data is
scale dependent in some instances, e.g., on small-scale maps a city will be represented as point
data, whereas on large-scale map, it will be represented as area data. A closed line data can be
converted easily to area data in GIS. For example forest boundary data can be converted to forest
land-use map. Point data cannot be converted in to area data. But reverse is true, i.e. area data
can be converted as point data.
The first two types are linked to areal aspects, which will be examined later. The thalweg is itself
to be considered not only as represented in horizontal projection and longitudinal profile, but
also by the way in which it combines with other thalwegs to form a drainage network, which has
its own physiographic characteristics. Some drainage network characteristics are linear, for
example, the bifurcation ratio, while others are areal in nature, such as the drainage density.
The first-order streams are the channels without any tributary; they receive all the flow
from surface overland route. Two first order channels join to form a second-order channel. A
second-order channel receives flow from the first-order channels which form it and the overland
flow. Evidently, a second-order channel carries more flow than a first-order channel. A third-
order channel is formed by the junction of two second-order channels; it receives flow from the
two second-order channels that form it, from direct overland flow, and from first- or second-
order channel(s) that might join it. Thus, a stream of any order has two or more tributaries of the
next lower order and so on. This scheme of ordering of streams is referred to as Horton-Strahler
ordering scheme.
Distance (km)
Fig. 8.1 Stream profiles of a typical Himalyan river and its tributary.
8.4 Area Data
The basin or catchment or watershed area at a point is defined as the area that receives
precipitation and other inputs and, after hydrological processes contributes to runoff at that point.
The watershed boundary or the ridge, directs any precipitation falling within its area towards the
outlet, whereas any precipitation falling outside the boundary drains to a different outlet. The
watershed is usually defined by using contour maps small scale.
The basin perimeter is measured in a GIS or with a curvometer. On a paper map, it is
determined by planimeter. The measured perimeter is a function of the scale and accuracy of the
maps or photographs, the quality of the curvometer, and the care taken in its use Many GIS have
routines to automatically delineate catchment areas by using digital elevation models.
A basin’s physical characteristics are the soil types, the land use and land cover (for
example, crops lakes, swamps, or glaciers), and the type of land use (for example, rural or
urbanized areas, lakes, or swamps). These physical features may be compiled as layers within a
GIS.
Ordinary maps do not give bathymetric data on lakes and reservoirs. Bathymetric survey
is carried out by using a boat and sounding by using echo sounders or tapes. The depths are
referenced to the national or an arbitrary datum. Volume is computed by employing the
trapezoidal, prismoidal, or any other suitable formula.
8.6 Topographic Maps
A topographic map is a detailed representation of the objects present on the surface of the earth
at a particular scale. Conventionally, topographic maps show the land features such as roads,
railway lines, power and other utilities, rivers and lakes, habited areas, etc. Elevations and
changes in it are shown by means of contour lines. Contours are the curves that connect points of
the same elevations.
Topographic maps provide the information about a terrain and thus are very useful for
water resources studies. The Survey of India has carried out extensive survey for almost the
whole India and prepared topographic maps at various scales. To identify a map of a particular
area, a map numbering system has been adopted by Survey of India. The system of identification
is as described here.
For the purpose of nomenclature or reference of the maps, an International Map Series (in
the Latitude range 4° N to 40° N and Longitude range 44° E to 124° E) at the scale of
1:1,000,000 is considered as the base. This base is divided into grids of 4° x 4° and these are
numbered from 1 (at the extreme north-west) to 136. Only land areas are covered and if any 4°
square falls completely in the sea, the area is not covered.
To refer to any Indian Topographic map, each 4° x 4° region is further divided into 16
grids, each covering 1° latitude x 1° longitude. The 1° x 1° grid in the North-West corner is
assigned the letter A and the grids are assigned consecutive letters, column-wise. Thus the letter
for the last grid on the South-East corner will be P. Due to geographic coverage, these maps are
also called degree sheets and are on a scale 1:250,000. Contours are drawn in these maps at
intervals of 100m. To refer to a particular map, the number of the base map grid and the letter
code is used, for example 39N (Fig. 8.2). Approximate area covered by each such map is 11140
sq. km. Maps in this scale are highly useful for planning studies from different sectors, viz.,
water, transportation, administrative.
Fig. 8.2 Nomenclatu
ure scheme of
o 1:250,000 topographicc maps.
The
T 1:250,000 scale maps or degree sheetss are fuurther sub-divvided in twoo ways:
1) Each
E sheet iss divided intto four partss (2 rows byy 2 columnss), each of 330' latitude x 30'
lo
ongitude (1:100,000) deesignated by y cardinal ddirections N NW, NE, SW W, and SE. Such
sh
heets are ideentified, for example
e as 39
3 P/NE (Figg. 8.3).
Approximate
A area covereed by each mapm is 700 ssq. km. Theese are generral purpose maps
an
nd are used by adminisstrators, plan nners, and eengineers. Inn fact, it is tthe most poopular
map
m series for all activitiees of the Gov
vt.
Each
E 1:50,00
00 scale sheeet contains four (2 row
ws by 2 coluumns) 1:25,0000 sheet (77' 1/2
latitude x 7' 1/2 longitude) which
h are numberred NW, NE E, SW, and SSE. Such sheeets are idenntified
as 53 O/114/NE. Apprroximate areea covered by y each such map is 175 sq. km. Thiss is more dettailed
map and is used for planning
p andd developmeent purposes..
Most
M topograaphic maps are fairly old these maay have to bbe updated aas required for a
particularr project. Fo
or large scalee maps, furth
her surveyingg needs to be carried outt.
CARTOSAT – 1 is the first Indian Remote Sensing Satellite capable of providing in-orbit
stereo images. The images can be used for cartographic applications. Cameras of this satellite
have a resolution of 2.5m (http://www.isro.org/satellites/earthobservationsatellites.aspx). The
Cartosat – 1 provides stereo pairs required for generating Digital Elevation Models, Ortho Image
products, and Value added products for various applications of Geographical Information System
(GIS).
An institute under the Indian Council of Agricultural Research (ICAR), the National
Bureau of Soil Survey and Land Use Planning (NBSS&LUP) situated at Nagpur is the premier
institute for soil related activities such as classification, mapping, uniform nomenclature, and
interpretation. Publications of NBSS&LUP include State soil maps (1:250,000 scale), District
soil resource maps (1:50,000 scale), and District land use atlas. More details can be obtained at
http://www.nbsslup.in/.
Maps of land use are prepared and maintained by Agriculture Department and revenue
authorities. These maps are extensively used in planning activities as well as in many other
sectors including water.
We will now briefly describe two proven technologies that are being increasingly used to
collect and manage spatial data. These are geographical information systems and remote sensing.
8.8 Remote
R Senssing
The termm Remote Sensing
S (RS)) denotes th he acquisitioon of data aabout an obbject withoutt any
physical contact betw ween the object and the sensing
s deviice. These daays, the termm RS techniqque is
chiefly used
u to denoote the acquiisition and analysis
a of ssatellite dataa. The main advantage oof the
RS technnology is thaat it providees a broad perspective
p oover a large area. One ccan "see" beeyond
visible electro-magn
e netic (EM) radiation baand, and daata of inacccessible areaas can be eeasily
obtained. Fig. 8.5 sh hows the EM M spectrum. Remote sennsing techniqques use the EM spectruum to
almost itts entire rangge. The interraction of EME radiationn with an objject can reveal a tremenndous
amount of
o informatio on about the object: Whaat is it? Wheere is it? Whhat are its phyysical properrties?
What arre its spatiaal relationsh hips with th he surrounddings? It is also possiible to infer the
characterristics of a to
op thin layerr of the earth
h’s surface.
Fig. 8.5
8 The elecctromagneticc spectrum.
Remote
R sensiing works on n the basic principle
p thatt depending upon the phhysical propeerties,
each object on the Earth
E surfacee reflects, emits
e and abbsorbs varyinng intensitiees of radiatiion at
different EM waveleengths. Usin ng informatiion from theese wavelenngth ranges, it is possibble to
nate between
discrimin n different types
t of grouund objects (e.g., waterr, dry soil, vvegetation, rrocks,
etc.) and
d map theirr distribution n on the grround. Inforrmation aboout the objeect concernned is
obtained by a senso or system lo ocated on a satellite, w which receives electrom magnetic radiiation
either em
mitted by the object or the radiation in
nteracted wiith the objecct.
The
T source of all sign nals collectted by mosst remote sensing insstruments iss the
electromagnetic radiiation (EM).. Sensors th hat rely on tthe Sun to generate alll the EM ennergy
needed too image terrrestrial surfaaces are callled passive ssensors. Sennsors that geenerate theirr own
energy arre called acttive sensors.
Remote sensing is a powerful technique to explore, map and manage the earth resources.
The images collected are processed and interpreted to produce useful data for application in a
wide variety of fields. The main advantage of remote sensing is that it gives a
broad perspective or bird’s eye view of the area. Depending upon the resolution, small to
very large area can be covered in the image. Further, one can "see" the objects in
wavelengths beyond the visible bands. Using remote sensing, images of
inaccessible regions can be easily obtained.
Airborne optical devices (Lidar) are now being used to determine topography more
accurately and at lower cost than conventional aerial photography.
Reflection occurs when radiation is neither absorbed nor transmitted. The reflection of the
energy depends on the properties of the object and surface roughness relative to the wavelength
of the incident radiation. Differences in surface properties allow to distinguish one object from
another.
where incident energy is the amount of incoming radiant energy and reflected energy is the
amount of energy bouncing off the object.
8.8.4 Re
eflectance
e Characterristics of
f Different
t Objects
Reflecta
ance is the perce entage of energy reflected to the total en nergy
incident
t on a body.
b It is the ratio
r betw
ween radiation off f the surrface
(outgoin
ng) to the
e total raadiation incident
i on
n the surf oming). We see
face (inco
an objeect green because only gre een energyy is refl lected inn the vissible
spectrum
m. Pure black
b or pure
p white
e are not colours. Other col lours such
h as
magenta,
, violet, yellow, cyan
c etc. are the mi
mixing of t ic colours RBG
three basi
in diff
ferent propportions. Reflectancce charact
teristics of some oof the fea
ature
are show
wn in Figu
ure 8.7.
Fig. 8.6
6 A typic
cal remote sensing data
d collec
ction and analysis s
system.
Itt is helpful too be familiaar with the baasic spectrall characteristtics of greenn vegetation,, soil,
and wateer as these arre of immensse use in watter sector.
Spectral Reflectancee of Soil. Soiil reflectancce typically iincreases wiith wavelenggth in the visible
portion of
o the spectruum and thenn stays relatiively constannt in the neaar-IR and shhortwave IR,, with
some loccal dips due to water abssorption at 1.4
1 and 1.9 µ µm and due to clay absoorption at 1.44 and
2.2 µm.
To acquire remote sensing data, polar orbiting or sun-synchronous satellites are used. A
satellite orbit in which the orbital plane is near polar (pass through N and S poles) and the
altitude is such that the satellite passes over all places on earth having the same latitude twice in
each orbit at the same local sun-time. Through these satellites, the entire globe is covered on
regular basis and it gives the repetitive coverage on periodic basis. Satellites such as Landsat
series, SPOT series, IRS series, may be grouped in this category.
An image can be analyzed using either the visual or the digital techniques. Visual interpretation
has many limitations and digital image processing (DIP) techniques are almost universally
applied these days. DIP involves the processing of data of an image to improve the image qualities
or to enhance the features of interest with the aid of a computer. The process helps in maximising
clarity, sharpness, and details of objects of interest and leads to better information extraction. It
improves the image’s interpretability. Image processing operations are carried out to remove noise
from the data and enhance certain features based on their spectral response. DIP is a broad subject
and may involve procedures that can be simple as well as quite complex. Basically, each pixel of
an image is mathematically manipulated and the operation may involve more than one image. The
results of computations for each pixel are stored and form a new image. The new digital image
may be subject to further analysis, may be stored or printed.
IRS-1C and IRS-1D were launched in 1995 and 1997 and were the first of the second generation,
operational, multi-sensor satellite missions with better resolution. These were placed in near-polar,
sun-synchronous orbit and had resolution of 5.2-5.8 m and 188 m, respectively. The repetitive time
of the satellites was 22 days, and on-board sensors had with Panochromatic (PAN) and LISS-III
cameras.
IRS-P3 was launched in 1996 on PSLV-D3, had two imaging sensing sensors and one non-
imaging sensors, viz. WiFs, with a resolution of 188 m and swath of 810 km.
IRS- P4, also called as Oceansat, primarily caters to oceanographic applications. The indigenous
Polar Satellite Launch Vehicle (PSLV) launched this satellite in 1999 into polar Sun-synchronous
orbit. The satellite has a high receptivity of 2 days. The payload includes an Ocean Colour Monitor
(OCM), a Multi frequency Scanning Microwave Radiometer (MSMR) and solid state memory for
recording data outside the visibility of a ground station.
IRS- P6 or RESOURCESAT-1 was launched into a polar sun synchronous orbit on October 17,
2003. The satellite carried LISS-3, AWiFS (Advanced Wide Field Sensor), and LISS-4 cameras.
LISS-3 camera had improved spatial resolution in shortwave infrared (SWIR). LISS- 4 was the
new high spatial resolution sensor. The sensor operates in three spectral bands with 5.8 metre
spatial resolution and steerable up to ± 26˚ across track to obtain stereoscopic imagery and
achieve five day revisit capability. LISS-3 operates in three spectral bands in VNIR and one
SWIR band; and an AWiFS operating in three spectral bands in VNIR and one band in SWIR
with 23.5 m spatial resolution. The payloads of this satellite have greatly aided crop/vegetation
and integrated land and water resources related applications.
IRS-P5 (CARTOSAT-1) was launched on May 05, 2005. It has an improved sensor system that
provides 2.5 m resolution with fore-aft stereo capability. This mission caters to the needs of
cartographers and terrain modelling applications. The satellite provides cadastral level information
up to 1:5000 scale and will be useful for making 2.5 m contour maps.
Researchers, practitioners and students in India are extensively using remote sensing
products to obtain spatial data required for water resources studies these days.
8.9 Geographical Information Systems
A Geographic Information System (GIS) is a computer-based system that is used in input,
storage, analysis manipulation, retrieval, and output, of spatial data. These systems consist of
computer hardware and software. GIS are increasingly being used in applications in natural
resources management. These days GIS are also being integrated with remote sensing and GPS.
GIS has origin in manual overlay operations done as early as in 1920s. In 1958, computer
based cartography initiated in USA which culminated in development of first general purpose
mapping software in 1960s. Canada GIS is also cited as first GIS and was developed around the
same time. Present day popular commercial/ open source GIS are: ArcGIS, GRASS, MapInfo,
ERDAS, IDRISI, ILWIS etc. GRASS (Geographic Resource Analysis Support System) is high-
end open source software. Arc GIS is modular high-end commercial software. ERDAS and
ILWIS packages have image processing and GIS capabilities.
GIS systems are now available for standard computers in practical, low-cost formats. The
main cost factor now resides in database compilation, and training and updating of technical
staff. Data capture or input in GIS is costly but it is one-time affair. Commercially available
paper maps may be cheaper than GIS layer. Use of GIS requires investment in computer
hardware, software, and training. GIS handling requires trained manpower. GIS software should
have proper functionality as desired in an application. For example, for hydrological modeling,
DEM analysis should be available in the package.
In raster data, points and lines are represented with finite area and finite width and thus is
not a natural representation. Lines have stepped or zagged appearance. In vector model, points
and lines have infinitesimal area and width respectively. Lines are smooth curves. Raster data
require large storage space. Vector data require small storage space. Thematic maps prepared
from remotely sensed data are available is raster form and are often processed as such. Many
hydrological models use both the representations. For example, thematic maps of catchment
variables and hydrometeorolgical measurements are prepared in raster form. Stream network is
processed in vector form etc. In raster form, value of many catchment variables is scale
dependent. For example, average slope of catchment reduces with increase in raster grid size. In
most GIS, the representations coexist. For example, it is better to capture spatial data from
conventional thematic maps, through visual interpretation of remotely sensed data in vector
form. Thematic maps from digital processing of satellite data may be obtained in raster form.
Topology: Method of representing vector data is called its topology. A line consists of two nodes
and one or more vertices. Nodes are end points of the line. Lines also have directions. Thus,
nodes are referred as ‘from node’ and ‘to node’ depending on direction of the line. Areas are
represented by ‘left area’ and ‘right area’ of each line.
TIN
In TIN model, elevations at the vertices of triangles are used to compute elevation at interior
points of the triangles. Using elevation of the vertices of a triangle, a planner or higher order
surface can be fitted. The surface can be used to derive elevation at points inside the triangle.
TIN model requires Delaunay triangulation. In this, constituent triangles are as equilateral as
possible. Circum-circles of the triangles include no other point of the triangulation. Triangulation
is performed first by constructing Voronoi diagram (Thiessen polygons). Points included in
adjacent polygons are joined to create Delaunay triangulation. Voronoi diagram is drawn using
proximity analysis.
Interpolation
Interpolation is a technique of determining unknown value of a variable at location from known
values at other locations. Interpolation can be used for any spatial variable, e.g. topographic
elevation, pH, SAR, pollutant concentration, groundwater depth and level, population etc.
Known values can be at point, line of area locations. Point data can be spot heights, pH, pollutant
concentration etc. Line data can be topographic contours, etc. Area data can be population
density in regions, etc.
Surface fitting: Here, n-degree polynomial surface is fit between selected known values. The
surface can be used, among the other application, to interpolate values.
Krigging: Variation of spatial variables can be partitioned in three components, namely drift or
structure, small variations and random noise. First component depicts general trend of the data.
Second component represents small variations from the general trend. These variations are
random but spatially autocorrelated. Third component depicts random values that are not
spatially autocorrelated. Kriging is best suited for interpolation of pollutant concentration,
geological and mining variables, e.g. grade of ores, etc. For these data, single smooth
mathematical equations are not suitable. The technique is based on assumption that values in
neighbourhood have generally higher correlated. Apart from the estimate of values, error
estimates are also provided in kriging technique. In presence of large random noise in data, good
semivariogram is not obtained and this results in deterioration in interpolation quality.
A map which is to be geo-referenced is called the source map and the reference map is
the map which has known coordinates. Points whose reference coordinates are known and which
are clearly identifiable on both the source and reference maps are known as control points. For
coordinates of control points in two maps, coefficients of a polynomial transformation equation
are estimated.
Storage: Geographic data are stored in GIS the native format of GIS. For one spatial data, many
computer files are created which contain different information. Attribute data are stored in Data
Base Management System and are linked to geographic objects. Storage of data in the form of
layers looks very attractive from water resources data. Different types of data, such as soil and
land cover, are stored in different map layers. GIS permits analysis of single or multiple layers
and various layers can be overlaid, one on top of another. From a water resources point of view,
spatial variation of data is important, e.g., the variation of soil hydraulic properties.
Statistics: for example, count, length, area, perimeter, shape, centroid, etc. for geographic objects
can be derived in GIS. For continuous surfaces, average, standard deviation, maximum,
minimum, etc. are derived. Summary operation produces zonal statistics for a map. For example,
land use statistics for watershed in a basin can be generated.
Logical operations: Logical operations, namely or, and, not, xor can be performed on maps.
Figure 8..8 shows log gical operatiions. For exaample, landuuse= agricullture and pH
H >= 8 will rresult
in salt aff
ffected agricu
ulture area.
Figure 8.8 Lo
ogical operaations commonly used inn GIS.
Classification: It con
nverts valuees into interv
val. A contiinuous surfaace is input and area m map is
output fo
or the operattion. In the output area map, isolinees, i.e., line of equal vaalues, enclosse the
area. Exaamples of vaarious isolin
nes are contoours, isohyette, isotherm,, isobar etc.,, which reprresent
topographic elevatiion, rainfalll, temperatture, pressuure, respecctively. In reclassificaation;
informatiion of geogrraphic objectt is changed.. For examplle soil seriess map may bbe changed too soil
pH map.
Search/ buffer:
b The operation iss similar to distance, exxcept that at a specified distance ann area
geographhic object is created.
Output: GIS output may be obtained as paper maps or as a picture which may be shared or
printed or inserted in a report or presentation. The output maps may contain various cartographic
elements namely title, legend, graticules or grids, north arrow, scale, annotations, notes, etc. In
one output more than one GIS layer may be included apart from cartographic elements. When
design is saved, it only contains reference to the layers. Thus, if a layer is modified and designed
output map is opened at a later time, the changes are reflected in the output.
The topography of a river basin may be represented in two different ways: as a digital
elevation model or as a triangulated irregular network (TIN). The digital elevation model is a
grid of elevation values that has regular spacing while TIN is a series of points linked into
triangular surfaces that approximate the surface. The spacing of points in TIN are non-uniform,
which allows points to be located on critical terrain features, roads or river banks. The accuracy
of such digital terrain models depends on the source of the data, the point density and
distribution, and other related data used in their development. Conventional contour maps may
be prepared from a digital elevation model or TIN. Orthophotos are images of the landscape
from which features can be referenced to one another. They are digital images produced by
processing aerial photography to geodetic control elevation data to remove all sources of
distortion. The image has the properties of scale and accuracy associated with a map. Such
images can be derived using airborne or satellite sensors.
8.9.6 Applications
Several hydrological software now have a GIS interface. GIS interfaces have been developed for
hydrological models such as HEC models, SWAT, etc. ArcInfo extensions such as Spatial
Analyst, 3D analyst, are useful for hydrological tasks. Script languages are also available to write
interfaces in GIS, for example, Arc Macro Language is the script language in ArcInfo.
GIS techniques facilitate input of spatial data to hydrological model. GIS are being
incorporated in hydrological models to extract and format distributed watershed data. Use of
DEM permits complete physiographic and hydrological depiction of basins. The efficiency of
handling large volumes of data means that more comprehensive and detailed maps, isolines and
themes can be prepared. This is a significant improvement in water resources management as
map preparation is often time-consuming and expensive.
GIS are becoming very common in the field of water resources assessment and
management. Many tasks of data collection, compilation, and interpretation can be facilitated by
means of GIS. In network planning and design, the ability to map quickly and display surface
water and related stations enables a more effective integration to take place. Network maps,
showing basins or stations selected according to record quality, watershed, or operational
characteristics, can be used for both short-term and long-term planning.
Groundwater potential and quality can be studied in GIS environment. Various layers
namely slope, geology, distances to drainage channel, tanks and lineaments, depth to water table,
depth of weathered zone can be overlaid and integrated on GIS environment to obtain
groundwater potential map. Similarly, layers of water quality variables may be created to obtain
quality map.
8.10 Closure
RS and GIS are powerful tools for spatial data collection and management. Spatial information
and its attributes can be stored, analysed and output efficiently through GIS. It has many
applications in natural resources and infrastructure management. GIS also has application in
varied industries namely power, transportation etc. Whereas in many applications, it can be
independently used, it also helpful in hydrological modeling in terms of data preparation. For
such applications, varieties of interfaces, stand-alone programs, embedded programs exist. These
increase productivity and reproducibility in modeling.
Conformal
Lambert conformal, one standard parallel Conic, tangent
Lambert conformal, two standard parallel Conic, secant
Stereographic Azimuthal
Mercator (e.g. Universal Transverse Mercator or UTM) Cylindrical
References
Engman, E.T. and R.J. Gurney (1991). Remote Sensing in Hydrology, Chapman and Hall, London.
Lillesand, T.M. and Kiefer, R.W. (1994). Remote Sensing and Image Interpretation. Wiley, New
York.
Sabins, F.F. (1987). Remote Sensing Principle and Interpretation, W.H. Freeman & Co., San
Francisco.
Module 9
Hydrological Information System
The primary role of a HIS (see Fig. 9.1) is to provide reliable data sets for planning,
design and management of water resource and for research activities. The system should function
in such a manner that it provides the information to users in time and in proper form. Sometimes,
the scope of HIS is extended to provide data to users on a real-time basis for short-term
forecasting or operational purposes.
The data collected for different hydrometeorological variables is called the raw or
observed data. The raw data have to be processed to ensure the reliability of the resulting
information. Both raw and processed data sets have to be properly stored -- processed data for
dissemination and raw data to permit inspection and revalidation in response to queries from
users. Note that the users have a central role in a HIS.
DATA FROM
OTHER SECTORS
HIS
USERS STORAGE
PROCESSED
DATA
Policy Makers and PROCESSING
Planners
Decision Makers
Designers
VALIDATION
Researchers
Other Users RECORDING
REAL
WORLD
Inadequate data may result in costly construction with an unnecessary high margin of
safety or with the equally costly consequences of failure. Ideally there must be a strong linkage
between data/information producers and users, who can be from within the same organization
and external users from other agencies, consultants, research scholars, etc.
The activities under HIS can be broadly classified in the following categories:
i. Assessing the user needs;
ii. Establishing an observational network and operating it;
iii. Data collection, validation, processing, and reporting;
iv. Management of historical data;
v. Data transmission, storage, and dissemination; and
vi. Institutional and human resource development.
In India, many Central and State Government organizations are the major users of the
hydrologic data. There are scores of other governmental, nongovernmental and private agencies
also who use this information. Usually there is a direct link between the objectives of water
resource management and the type of data needed from the Hydrological Information System.
Due to various reasons, hydrological data needs of the users change with time. Therefore, it is
very important to identify the potential data users and regularly update their data needs.
Normally, it is expected that the hydrological information service agencies will satisfy most of
the genuine data needs of the potential users. For ensuring an optimal use of the public resources
spent for maintaining a Hydrological Information Service, it is essential to have a proper balance
between the data needs of different users and mandate of various services supporting HIS.
Often, it is assumed that the information being provided is the same as required by the
potential users. User’s needs may change with time and the HIS would fulfill its obligations only
if there is a continual review of the changing needs of the users.
After finalizing the objectives of the HIS, the observational network has to be planned
and established. This aspect has been discussed in a previous module. It is important to ensure
that the observation networks of different agencies are properly integrated so that duplication of
efforts is avoided.
1. Manual with the observer at a station sending data to the central office.
2. Manual/semi-automatic system where the central office manually interrogates the
automatic field stations through telephone, radio, etc. and receives the data.
3. Pre-program and time system where automatic equipment initiates the transmission of
observations.
4. Automatic event indicator and the station automatically transmitting the specified change
of variable to a central location.
5. Automatic system with station transmitting and central office recording data
continuously.
a) Telephone lines are used wherever feasible. With improvements in information technology,
high rates of data transmission are possible.
b) High frequency radio links are used when land lines are not available or topography is
difficult. The installation cost can be high.
c) A significant development of last few decades has been the use of mobile telephony or
satellites for data transmission. A satellite-based system consists of Data Collection
Platforms (DCP) that are installed at hydrometeorological stations. DCPs are (rechargeable)
battery-operated devices that collect, encode, and communicate the data of the station to a
central location through a satellite link at pre-determined time/frequency. Solar panels are
frequently employed to generate energy. This system is very useful for remote and difficult-
to-access locations.
In some organisations, data collected by the whole network may be kept at a central
place. It may be noted that the above is a general descriptions and there may be case specific
differences.
It is clear from the above that data processing activities are carried out at more than one
level within each agency and this makes it essential to have adequate /communication links
between them with increasing penetrations internet, data transmission has become quite simple
and efficient. The requirement for communication is to be based on a low frequency and high
volume of communication. Information is exchanged between various agencies for data
validation as surface and groundwater networks are operated by different state and central
agencies.
9.4 Data Storage and Retrieval
The vast amounts of data which are observed by incurring huge efforts and resources should be
stored in such that they are easily obtainable and secure against losses. Also, it is much more
useful if the basic data is processed and provided to the end users in standard formats. Common
formats include ASCII files [(data are given as comma separated variables (CSV)] or as
spreadsheet files. This can save a lot of money and effort of the user agencies and they would be
encouraged to use the data to solve different water resources problems. Archival of historic data
is important in any field. When done in a logical manner, it enables the end users to make the
efficient use of data.
Vast quantities of water resources data are collected all over the world every day. In this
computer era, the archival of data may be accomplished in a very efficient and economic way.
The basic and other processed information may be stored on computer media and the same may
then be quickly made available to the users. Hard copies in the form of data year-books may also
be brought out for use by practising engineers, planners and managers. These water year-books
are frequently made available to the users on computer. Besides the processed data, it is
advisable to store the raw data as well which might be needed for research purposes and it might
be realized later that the data validation procedures had missed some aspects.
A typical setup for water resources data management is depicted in Fig. 9.2. The main
components are: a) data entry module to input data from various sources in the database, b) user
interface for data editing, display, and management, and c) applications that can retrieve data or
write to the data base.
Due to large volumes of water resources data, it is necessary that the data are stored such
that minimum space is needed. The current trend is to archive the data on optical/magnetic disks
which can hold large amount of data and do not require stringent environmental conditions for
storage. Now-a-days high capacity hard disks and pen-drives are available and the cost of
hardware has dropped drastically. Of course, the volume of data that are being generated each
day is also progressively increasing.
Data entry
Utility Application
programs 100110100111010110
programs
011011010001001001
Data compression is applied to reduce the requirement of storage space. WMO (1994)
has described such procedures. A number of data compression algorithms have been developed.
For example, daily rainfall is measured to an accuracy of 0.1mm. Rather than storing it as a real
number which requires 4 bytes, the values can be multiplied by 10 and stored as an integer
requiring only 2 bytes. The daily data of one month when rainfall was zero need not be stored as
30 zeros. An efficient way is to use notation ‘30*0’. The database files are not normal ASCII
files that require larger space; these are special types of files. The suitability of a particular
compression technique will depend on the characteristics of the data.
Efficient data search and retrieval module is also necessary so that the requisite data are
quickly fetched from the database. A good retrieval system should provide the user a
combination of options to select the data using criteria, such as by variable, basin, station, time
period or range of values. The user may also be given the option to choose the format of the
output, i.e. tabular, plot or ASCII data files that can be directly input to another software. If
graphs are displayed on the monitor, the user may be given option to print them or store on hard
disk for later use. Adequate security measures should also be built in the retrieval system so that
only authorized users have access to the database. Among the authorized users also, there could
be various categories. Most users are given read-only access and they cannot modify to the
database. A limited group of users are given all privileges, i.e., they can read, modify and delete
data from the database. It is useful to have a log of all users who have accessed the database and
the operations that they have performed so that the source can be identified in case of any
mishap.
Many organizations routinely publish basin-wise data-year books. A typical water year
book consists of description of the basin, its topography, soils, land use etc., major rivers, and
salient features of various water resources projects. Maps are included to illustrate all these
features. The data section contains typically precipitation, streamflow, evaporation, and ground
water data. Periodically, special reports may be published giving long-term statistics of stations
or highlighting special or unusual events, such as floods or droughts.
Of late, many organizations have begun to curtail hard copy publishing due to high costs
and handling problems. With almost all data now available in digital form in a HIS, the hard-
copy publication is not an efficient means of data dissemination. A water-year book can be
conveniently published on line as a PDF file. A browser may also be supplied to handle data
search, display and print. This trend is likely to accelerate further. An important thing to
remember is that the format and content of publications should depend on the need of users. The
contents should be so designed that the need of most data users are answered and the efforts to
handle data requests are reduced. The contents may also depend on the frequency of publication.
Finally, in this age of computers and Internet, it is appropriate that many international
organizations have established databases that can be accessed through Internet. The Global
Runoff Data Centre (GRDC) at the Federal Institute of Hydrology, Germany, has a large archive
of surface water data. The United States Geological Survey (USGS) maintains one of the largest
online hydrologic data base (www.usgs.gov).
Although an integrated HIS that contains data of all the hydrological variables is
desirable, currently in India separate databases are being maintained for surface water data,
ground water data, and hydro-meteorological data. Further, currently organizations belonging to
the Central government and State governments have separate data bases. An integrated data base
is desirable and will have many advantages:
1. An integrated data base will foster co-operation among the various Central and State
Organisations which will lead to improvement in the quality of the data,
2. It will be of great benefits to the users by providing all requisite information from a single
source,
3. It will improve the sustainability of the system and enhance the chances of a uniform
approach in the long term, and
4. It will lead to savings on the expenditure on the hardware, software and manpower.
9.5 Institutional and Human Resources development
Since HIS is a vast system, institutional and human resource development aspects are very
important. The institutions implementing of HIS must be developed in such a manner that the
system is sustainable in a long run. Trained staff is required to carry out different activities under
HIS. In India, observations are taken manually or automatically depending upon the type of
instrument at the station. Persons of various categories (e.g., Junior Engineers, Supervisors,
Technicians, Observers, Helpers, etc.) are engaged for carrying out data collection and for
regular maintenance. Very importantly, staff must have skills required to carry out the desired
tasks. Only those staff should be posted to these centers who have interest and temperament to
work in a data center. Uninterested or incapable people should never be posted in such places.
Regular trainings have to be organized to update the existing staff as well as to train the new
recruits.
Many countries are using the power and opportunities provided by computerized data
processing system and Internet. However, the potential has not been fully exploited in India in
the field of hydrology so far.
9.7 Closure
Two basic tools for integrated management of the environment are modeling and
environmental data. Both tools were available and valid in the past; however, the recent
requirements for integrated environmental management have also led to a significant evolution
of both modeling procedures and data management systems. Current literature provides vast
amounts of studies on modeling of different environmental processes. However, issues related
to data management systems are barely touched in a comprehensive framework. Data
requirements and data availability are mentioned merely as subtopics in most environmental
studies although it is well recognized that data constitute the basis for all environmental
management activities. Most developed countries have well-established databases which can
be accessed easily by the users. Developing countries, on the other hand, do not have extensive
data banks and these many not be easily accessible by the users. There is a need for
harmonization or standardization in development of databases so that data exchange can be
facilitated on regional and global levels.
9.8 References
Hydrology Project, “Basic Surface Water Data Processing, HIS (SW) Operation Manual Vol. 8 –
Data Processing”, Ministry of Water Resources, India, 2003.
Latkovich, V.J. and Leavesley, G.H. (1993). Automated Data Acquisition and Transmission, in
Handbook of Hydrology, edited by D.R. Maidment. McGraw-Hill Book Company, New
York.
WMO (2008). Guide to Hydrological Practices. WMO No. 168. World Meteorological
Organization, Geneva.
Module 10
Statistical Analysis of Hydrological Data
Before commencing the analysis, it is always useful to plot the sample data. The
following are the main steps of statistical analysis of data:
If all the values of a random variable and the corresponding probabilities are known or
found, the relation between these values and probabilities is described by a probability
distribution. Knowing this distribution, the probability of any value of the random variable can
be determined.
10.1.1 Probability
Based on the daily experiences, most people have an intuitive appreciation of the concept of
probability or chance. When a coin is thrown, there are two equally likely possible outcomes:
head or tail. Let the coin be tossed n times and head occurs s times. The ratio s/n is the
probability of occurrence of heads.
Sample space is a collection of all possible random events that might arise from a random
experiment. Sample space S and two sets of events A and B are shown in Fig. 10.1. If two events
A and B do not contain any common sample points, they are said to be mutually exclusive.
If two events A and B are not mutually exclusive, the common set is called their
intersection, denoted as A∩B. The union of two events A and B gives the event which is the
collection of all sample points occurring at least once in either A or B and is denoted as AB.
Fig. 10.2 shows this concept through the Venn diagram. The complement Ac of an event A
consists of all sample points in the sample space of the experiment not included in the event A.
A B
S
A A
B B
A∩B
AB
Fig. 10.2 Venn diagram showing intersection and union of two events.
The notation P[A] is used to denote the probability of a random event A. Now we discuss the
axioms of probability.
Axiom 1: The probability of an event A is a number greater than or equal to zero but less than or
equal to unity:
0 P[A] 1 (10.1)
where A is the event associated with all sample points in the sample space.
Axiom 3: The probability of an event which is the union of two events is:
where AB denotes the union of events A and B which means that either event A occurs or event
B occurs, and A B denotes the intersection of event A and event B. Eq. (4) can be extended to
the union of n events. If A and B are two mutually exclusive (disjointed) events, the probability
of A and B, P[A B], will be zero and eq. (5) becomes:
Statistical independence implies that the occurrence of event E1 has no influence on the
occurrence of event E2.
The probability distribution function (PDF) describes the relative likelihood that a
continuous random variable X takes on different values, and is the derivative of the CDF:
Let us now state some important properties of f(x) and F(x) for continuous and discrete
random variables.
f ( x) 0, x
(10.8)
2. The sum of probabilities of all possible outcomes is equal to 1, i.e., the area under the
PDF is unity.
f ( x ) dx 1
(10.9)
(a)
F(x)
F(x)
CDF
|
x X
f(x)
x
(b)
n
---
Nx
|
0 x X
Fig. 10.3 PDF and CDF of a random variable.
f (x ) 1
i
i (10.10)
where f(xi) represents the probability of X = xi in the sample space if the sample contains finite
and f(xi) observations can be replaced by p(xi). The probability of “X assumes a value x” is
equal to the area under PDF curve to the left of x:
x
3.
P ( X x) F ( X x) F ( x) f ( x)dx
(10.11)
Moreover
P ( a X b) P ( X b) P ( X a )
b a
f ( x)dx f ( x)dx
b
f ( x)dx, for a b (10.12)
a
k
P ( X x k ) p ( xi ) (10.14)
i 1
10.2 Key Statistical Measures of Data
Four principal features which characterize a set of observations on a random variable are:
(i) the central tendency or the value around which all other values are bunched,
(ii) the spread of the sample data around mean,
(iii) the asymmetry or skewness of the spread of data, and
(iv) the peakedness of the data.
These characteristics are expressed in terms of statistical properties which are estimated
from the sample data.
(i) Arithmetic Mean: If x1, x2 … xn represent a series of observations, the mean of this series is:
1 n
x = xi
n i=1 (10.15)
Where x represents the sample mean; the mean of population is generally denoted by .
(ii) Mode: It is the value which occurs most frequently. It is the peak value of the PDF. A data
set may have more than one peak.
(iii) Median: It is the middle value of the ranked observations for a data set. The median divides
the distribution in two equal parts.
(i) Variance: It represents the scatter of the data are about the mean. Variance is computed by:
1 n
s = (xi - x )
2 2
(10.16)
n i =1
A small value of variance implies that values are bunching close to the mean.
(ii) Standard Deviation (SD): The unbiased estimate of population standard deviation (s) is given
computed as the square root of the variance:
1 n
s= [
n i=1
( xi - x )2 ] 0.5 (10.17)
when n < 30, the unbiased estimate of s is found by replacing n by n-1 in the denominator. Greek
letter σ is used to denote the standard deviation of population.
(iii) Coefficient of Variation (CV) is a dimensionless parameter and is obtained by dividing the
standard deviation by the mean:
CV = s / x (10.18)
When the mean of the data is zero, Cv is not defined. This coefficient is useful to compare
different populations. Given two samples of data, the one with larger Cv will have more spread of
the values around the mean.
Example 10.1: Average annual flows (in cumec) at a river gauging site are given in the table
below. Compute the mean, variance, standard deviation, and the coefficient of variation of the
flows.
Year 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980
Flow (cumec) 195.5 145.4 148.1 324.7 205.6 302.9 210.3 194.4 71.2 126.8 216.0
Year 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991
Flow (cumec) 136.4 403.9 145.3 161.6 112.5 110.0 90.2 129.8 80.5 136.3 243.3
Solution: We have a total of 22 values. The mean of the flows can be computed as
The third moment of the data about the mean is used in indicating symmetry and is given
by:
n
1
M3 =
n
(x
i =1
i - x )3
(10.19)
It is easy to see that this moment is zero if the data are symmetrical. Otherwise, M3 will have
certain value, a positive or negative.
Note that because the third central moment has dimensions equal to the cube of the data,
it is not useful while comparing different data sets. Being non-dimensional the coefficient of
skewness does not have this disadvantage and is preferred.
Symmetrical frequency distributions have very small or negligible value for skewness
coefficient Cs, while asymmetrical frequency distributions have either positive or negative
coefficients. When Cs has a small value, it indicates that the probability distribution may be
approximated by the normal distribution since Cs = 0 for this distribution. The symmetrical and
skewed distributions are shown in Fig. 10.4.
f(x1)
f(x2)
f(x3)
x --
Fig. 10.4 Symmetrical and asymmetrical (+ve and –ve) skewed distributions.
Normal distribution has the kurtosis 3. If a data set has a relatively greater concentration near the
mean than the normal distribution, the kurtosis will be greater than 3. Conversely, if the data
have a relatively smaller concentration near the mean than the normal distribution, the kurtosis
will be less than 3.
Example 10.2: Compute the coefficient of skewness and the coefficient kurtosis of the data of
example 10.1.
where w is the size of class interval. In general the number of classes varies between and 25.
(vii) Plot the relative frequencies as well as cumulative relative frequencies with group
interval as abscissa and the relative frequencies or cumulative relative frequencies as
ordinate.
Example 10.3: The annual flow of Sabarmati River at Dharoi is plotted in Fig. 10.5 for the
period 18
868-1965. Pllot the histog
gram and thee cumulativee histogram.
The cum
mulative histoogram of thee annual flow
w of Sabarm
mati River att Dharoi for the period 11868-
1965 (98 years) is plo
otted in Fig. 10.6.
Fig. 10.6
6 Histogram
m of the annual river flow
ws.
A probability distribution gives important information about the data, how the values are
changing, whether they are bunched together or spread out, and whether they are symmetrically
disposed on the X-axis or not. Distribution also tells the relative frequency or proportion of
various X values in the population in the same way that a histogram gives information about a
sample. We now describe the distributions that are commonly used in addressing water resources
problems.
Commonly used distributions in hydrology are the Normal, Log Normal, Extreme Value
type-1 (Gumbel or EV1), Gamma, Pearson Type - III, and Log Pearson Type - III distributions.
A brief description of these distributions is given below.
1
x
(x )2
F(x) = 2 2 du
2 -
exp (10.25)
The origin of the normal organization lies in the central limit theorem which states that if
a sequence of random variables xi, i = 1, 2, …, n are independently and identically distributed
n
with mean μ and standard deviation σ then the distribution of n such random variables Y= x
i 1
i
tends to the normal distribution with mean nμ and standard deviation n , as n becomes large.
This theorem holds good irrespective of the probability distribution of x.
The reduced variate of the normal distribution is defined as Z = (x - )/. The properties
of the reduced variate are mean = 0, standard deviation z = 1, and coefficient of skewness = 0.
Fig. 10.8 shows the normal distribution and the area under the curve for three values of the
reduced variate. As shown, the area under the curve within ± is 68.27%, within ± 2 is
95.45 and it is 99.73 within ± 3.
±
Area = 68.27%
f(x)
± 2, area =
95 45%
± 3, area = 99.73%
| | | | | | |
x
Fig. 10.8 The normal distribution and the area for three values of the standard variate.
Among the probability distributions used in hydrology, the normal distribution is the
most widely. It is also employed in the analysis of variance, estimation of random errors of
hydrologic measurements, hypothesis testing, synthetic generation of random numbers, etc. A
random variable that is made up of the sum of many small independent effects is expected to
follow a normal distribution. Many hydrologic variables are not normally distributed, but
transformations can, in many cases, make them approximately normally distributed. When there
is increase in the time interval over which a hydrologic variable is measured, the variable
approximately follows a normal distribution because the number of causative effects increases.
Example 10.4: Assuming that the data of Sabarmati River follows the normal distribution, find
the parameters of the distribution and plot it.
n: For the datta of Sabarm
Solution mati river, thee mean and S
SD are:
Mean
M of the data
d x = 6655.37 million cubic m.
Standard deviiation = 34
46.9 million cubic m.
Coefficient
C off variation CV = 346.9/665.37 = 0.5221.
Coefficient
C off skewness Cs = 0.76 (poositively skeewed).
Kurtosis
K Ck = 3.65.
Fig. 10.9
9 shows the plot
p of the prrobability distribution off Sabarmati R
River data.
9 Probability
Fig. 10.9 y distribution
n of Sabarmaati River datta.
For the Sabarmati data, the coefficieent of skewn wness Cs is 00.76 or the ddata is posittively
skewed. This is easilly verified by
y Fig. 10.9. Further, kurrtosis Ck forr the data is 33.65 (kurtossis for
the norm
mal distributio
on is 3). Agaain, this can also be veriffied from Fiig. 10.9.
The log-normal distribution has two parameters y and y which can be estimated by
transforming all xi’s to yi’s by
yi = ln xi (10.27)
f(x) = exp {-(x - ) – exp[-(x - )]} - < x < ; - < < ; > 0
where , and are scale and location parameters. The estimates of parameters using the method
of moments are:
1.283
ˆ ; ˆ x - 0.45s (10.29)
s
It has been shown that the EVI distribution is a special case of a distribution known as the
Generalized Extreme Value (GEV) distribution. The CDF of the GEV distribution is given by
x u
1/ k
F ( x) exp 1 k (10.30)
Where k, u, and α are the parameters of the distribution. When k = 0, we get the EV1
distribution. For k<0, the distribution known as EV2 and it is known as EV3 distribution when k
> 0.
According to Gumbel, the probability that an event with magnitude larger than x0 occurs
is (Subramanya 2008):
XT = x + kσx (10.35)
Equations (10.35) and (10.36) assume that a large data series are available to compute the
various statistics. However in practice, the record length is finite. In such cases, the following
equation may be used:
xT x k n n 1 (10.37)
where σn-1 is the standard deviation of the sample of size n. Frequency factor for use with sample
of size n is given as
k n yT y n / s n (10.38)
where y n and sn are reduced mean and reduced standard deviations which are functions of n.
Values of these can be obtained from standard tables that are widely available (see, for example,
Subramanya 2008). Note that as n➝∞, y n ➝ 0.577 and sn ➝ 1.2825.
Example 10.5: From the flow data of a river, the mean and standard deviation were computed
and these turned out to be 660 million cubic m and 330 million cubic m, respectively. Find the
value of parameters of EV1 distribution.
Solution: The mean and standard deviation of the data are 665.37 million cubic m and 346.9
million cubic m, respectively. Therefore, the estimates by the method of moment are:
= 1.2825/330 = 0.0039.
and = 660 – 0.45*330 = 511.5.
Example 10.6: Annual maximum flood discharge data of a river was available for 30 years.
Mean and standard deviation were 5250 m3/s and 1650 m3/s. Compute the flood discharge with a
return period of 100 years by using the Gumbel Extreme Value 1 distribution.
where a, b, and c are scale, shape, and location parameters, respectively, and (b) is a gamma
function. If c = 0, this distribution becomes a two-parameter gamma distribution. Parameters a,
b, and c are related to mean, standard deviation, and coefficient of skewness as (method of
moment estimates)
a = /b (10.39a)
2
b = (2/Cs) (10.39b)
c = - b (10.39c)
To determine flood for a return period T by using the LP3 distribution, the procedure
described below is followed.
Where k = Cs/6. To complete z for a given return period T, exceedance probability p is obtained
as p = 1/T. Now, complete a variable w as
w ln(1 / p 2 ) 0 p 0 .5
when p > 0.5, p in eq. (10.41) is replaced by (1-p) and the negative sign is put before z computed
by eq. (10.42). Now, by following the frequency factor method, the flood for the return period T
years is computed by:
yT y K T s y (10.42)
Example 10.7: For the data of Example 10.1, find the parameters of the Pearson Type III
distribution.
Example 10.8: Logarithms of the annual flood peak data of a river were taken and the mean was
4.146, SD was 0.403 and Cs = -0.07. Find 50 year return period flood by using the LP3
distribution.
r n- r n! r n- r
P r,n = n C r P q = p q (10.43)
r! (n - r)!
where Pr,n is the probability of a random event of a given magnitude and exceedance probability
P occurring r times in n successive years. The probability of the event not occurring at all in n
successive years is:
P1 = 1 - qn = 1 - (1 - p)n (10.45)
Example 10.9: An analysis of data on the maximum one-day rainfall depth at a station indicated
that a depth of 280 mm had a return period of 50 years. Determine the probability of a one-day
rainfall depth equal to or greater than 280 mm occurring (a) once in 20 successive years, and (b)
two times in 15 successive years.
Example 10.10: What is the probability that a 5-year flood will not occur at all in a 10-year
period?
x exp( )
P( X x) , 0, x 0, 1, 2,... (10.46)
x!
The conditions for application of Poisson distribution are: a) the number of events is
discrete, b) two events cannot coincide, c) the mean number of events per unit time is constant,
and d) events are independent. Thus, it can be applied to following situations with p relatively
small and n relatively large to determine the probability of:
Let X be a continuous variable and f(x) its function satisfying some necessary conditions.
The r moment of f(x) about an arbitrary point ‘a’ is denoted as M ra (f). The rth moment of the
th
Fig. 10.10 shows the definition of various terms used in the above equation.
Thus, the zero-order moment is the area under the curve defined by f(x) subject to - <x < .
For probability distribution, this area is unity. If r = 1, eq. (10.48) yields
M 1a ( x a )1 f ( x) dx a
(10.50)
where is the mean. If the moment is taken around the origin, then a = 0, and the first moment
gives the mean. When a = , the rth moment about the mean is expressed by
M r ( x ) r f ( x) dx
(10.51)
f(x)
a x-a
- 0 d x
x
For simplicity of notation, we drop the superscript if the moment is taken about the origin
0 and the familiar terminology of the moments can be written as follows:
M0 =Area
M1 = Mean
M2 = Variance,
M3 = Measurement of skewness of the function
M4 = Kurtosis,
These terms have already been defined earlier.
Example 10.11: The frequency table of annual flows of Sabarmati River is given Table 10.2.
Find the mean and variance of the data by using the method of moments.
f(x
f -4) f(x-3) f(x-2) f(x-1) f(x
f 1) f(x2) f(x3) f(x4)
… -3 -2 -1 1 2 3 …
-
x-2 x2 x
x-3 x3
Disch
harge range Frequenccy Dischaarge range Frequenccy
10
00-200 6 2000-300 9
30
00-400 11 4000-500 9
50
00-600 9 6000-700 9
70
00-800 19 8000-900 6
900-1000 6 10000-1100 1
110
00-1200 5 12000-1300 2
130
00-1400 3 14000-1500 0
150
00-1600 2 16000-1700 0
170
00-1800 1
Solutionn: The sum ofo all frequenncies = (6+1 11+9+19+ ……) = 98. For the first raange of dischharge
(100-200 0), the mean value is 150
0; for the seccond range (2200-300), it is 250 and sso on.
Hence, thhe first mom
ment of the daata = (150*6 6 + 250*9 + 350*11 + … …+ 1750*1)//98
= 664.2 29 cumec.
This is th
he mean (xm) of the annuual flows.
The seco
ond moment about the mean
m gives th
he variance oof the data.
N
( xi x m ) 2 f ( xi )
Second M
2
i 1 N moment
Let f(x; a1, a2,… am) be a PDF of the random variable X with parameters ai, i=1, 2, …, m,
to be estimated. For a random sample of data x1, x2, …xn, drawn from this probability density
function, the joint PDF is defined as
n
f ( xi , x 2 ,....x n ; a1 , a 2 ,....a m ) f ( xi ; a1 , a 2 ,....a m )
i 1 (10.53)
By maximizing the likelihood that the sample under consideration is the one that would
be obtained if n random observations were selected from f(x; a1, a2, … am), the unknown
parameters are determined, and hence the name of the method. The values of parameters so
obtained are known as MLE estimators. Since the logarithm of L attains its maximum for the
same values of ai, i = 1, 2, … m, as does L, the MLE function can also be expressed as
n n
ln L L* ln f ( x1 ; a1 , a2 , .... am ) ln f ( x1 ; a1 , a2 , .... am )
i 1
i 1
(10.55)
Frequently ln[L] is maximized, for it is many times easier to find the maximum of the
logarithm of the maximum likelihood function than that of the normal L.
The procedure for estimating parameters or determining the point where the MLE
function achieves its maximum involves differentiating L or ln L partially with respect to each
parameter and equating each differential to zero. This results in as many equations as the number
of unknown parameters. For m unknown parameters, we get
L(a1 , a2 , ... am )
0
a1
L(a1 , a2 , ... am )
(10.56)
0
a2
L(a1 , a2 , ... am )
0
am
These m equations in m unknowns are then solved for the m unknown parameters.
Example 10.12: Using the method of maximum likelihood, find the parameter of the
exponential distribution for the data of the Sabarmati River in India, given in Example 10.1.
Solution: The probability density function of the one-parameter exponential distribution is given
by
fX(x) = exp(-x) (10.57)
n
ln L( ) n ln( ) xi
i 1 (10.59)
In Example 10.1, the mean of the data was found to be 664.29 cumec. Hence, the estimate of
is:
= 1/664.29 = 1.51×10-3 cumec-1.
Let there be a function Y = f(X; a1, a2,… am), where ai, i = 1, 2, … m, are parameters to
be estimated. The method of least squares (MOLS) involves estimating parameters by
minimizing the sum of squares of all deviations between observed and computed values of Y.
Mathematically, this sum D can be expressed as
n n n
D d12 y0 i yc i y0 i f x; a1 , a2 ,...am
2 2
i 1 i 1 i 1 (10.61)
where y0(i) is the ith observed value of Y, yc(i) is the ith computed value of Y, and n > m is the
number of observations. The minimum of D in eq. (71) can be obtained by differentiating D
partially with respect to each parameter and equating each differential to zero, e.g.,
n
y 0 i f xi ; a1 , a 2 ,...a m
2
i 1
0 (10.62)
a1
The resulting m equations, usually called the normal equations, are then solved for estimation of
m parameters. This method is frequently used to estimate parameters of linear regression model.
Let us consider a sample of size n and arrange the data such that X(i|n) is the ith largest
observation; clearly i = n will be the largest value. Then, for any distribution, the second L-
moment, λ2, is a description of scale based upon the expected difference between two randomly
selected observations:
To compute L-moment measures of skewness three randomly selected observations are used and
for kurtosis, we use four randomly selected observations.
Sample L-moment estimates are often computed using (PWMs). The rth PWM is defined
(Loucks and Beek, 2005) as:
1 n
j 1 n 1
br
j r 1 r
X j
n r
(10.69)
1 n j 1 n
X j
r 1 j r 1 r
r 1
for r = 1, … , n -1.
As with traditional product moments, measures of the coefficient of variation, skewness and
kurtosis of a distribution can be computed with L-moments. Following L-moment ratios are
important:
Example 10.13: Table 10.3 gives annual discharge data of a river for 36 years. Compute sample
L-moments and L-moment ratios, L-CV, L-sk, and L-ku.
Table 10.3 Annual discharge data of a river for 36 years
Year Discharge Year Discharge Year Discharge Year Discharge
1950 400 1959 1390 1968 2291 1977 1499
1951 1100 1960 3300 1969 1340 1978 2598
1952 900 1961 2190 1970 3200 1979 3487
1953 440 1962 935 1971 2200 1980 1234
1954 3000 1963 785 1972 1014 1981 819
1955 2500 1964 501 1973 1790 1982 1210
1956 760 1965 1123 1974 1140 1983 1510
1957 1250 1966 1581 1975 764 1984 1780
1958 1340 1967 959 1976 783 1985 1398
Solution: Equation (10.68) yields estimates of the first three Probability Weighted Moments:
b0 = 1514.19
b1 = 889.16
b2 = 655.38
b3 = 518.64 (10.71)
The sample L-moments can be calculated using the probability weighted moments to obtain:
1 = b0 = 1514.19
2 = 2b1 - b0 = 264.12
3 = 6b2 - 6b1 + b0 = 111.53 (10.72)
4 = 20b3 - 30b2 + 12b1 - b0 = -132.82
Thus, the sample estimates of the L-coefficient of variation, t2, and L-coefficient of skewness, t3,
are:
t2 = 264.12/1514.19 = 0.174
t3 = 111.53/264.12 = 0.422 (10.73)
t4 = -132.82/264.12 = -0.502
10.6 PROBLEMS OF PARAMETER ESTIMATION
The parameters of a distribution function are estimated from the available sample data. But while
doing so, errors may arise due to many reasons. The sample data may contain errors, the
assumption underlying a particular method of parameter estimation may not hold good, and there
may be truncation and round-off errors. All these may result in errors in estimates of parameter.
Each estimate of a parameter is a function of sample parameter data which are observations of a
random variable. Thus, the estimate value of the parameter itself is a random variable with
certain distribution. An estimate obtained from a given set of values can be regarded as an
observed value of the random variable. Thus, the goodness of an estimate can be judged from its
distribution.
Several questions arise in parameter estimation. How should we employ the available
data to obtain the best estimate? What is the best estimate? Are these estimates unique? A
number of statistical properties are available by which to address the above questions. These are
discussed below.
Bias
Let the estimate of parameter a be ac denoted by. Estimate ac will be called an unbiased estimate
of a if the expected value of a, denoted E (ac) = a. In general, an estimate will have a certain bias
b(a) depending on a so that
An unbiased estimate mean b(a) = 0. Note that an individual ac may not be equal to or
close to a even if b(a) = 0. Unbiasedness simply implies that the average of many independent
estimates of a will be equal to a.
The bias in a given quantity is usually measured in dimensionless terms and is often
referred to as standardized bias (or BIAS). Thus, BIAS is defined as
E aˆ a
BIAS (10.75)
a
where n is the number of samples generated or the number of values of the a estimate. The value
of a in eq. (10.75) is the true value of a or the value of parameter a of the population.
Efficiency
An estimate ac of a is said to be efficient if it is unbiased and its variance is at least as small as
that of any other unbiased estimate of a. If there are two estimates of a, say a1 and a2, then the
relative efficiency of a1 with respect to a2 is defined as
E (a1 a ) 2
e 1 (10.77)
E (a2 a) 2
Standard Error
Another dimensionless performance measure frequently used in hydrology is the standard error
(SE), defined as
SE aˆ / a (10.78)
12
1 n 2
aˆ aˆ i E aˆ i
n 1 i 1 (10.79)
where the summations are over n estimates â of a. In Monte Carlo experiments, referred to as
above, for each sample size, a value of SE is obtained. Thus, this measure is similar to the
coefficient of variation.
where E[. ] is the expectation of [. ]. It can be shown that RMSE is related to BIAS and SE as
12
n 1 2
RMSE SE BIAS 2
n (10.81)
RME 0
(10.82)
N i 1 Q0
in which N is the sample size, Q is the observed quantity of a given probability and Qc is the
computed quantity of the same probability. Also, used sometimes is the relative absolute error
defined as
1 N
Q0 Qc
RAE
N
i 1 Qc (10.83)
10.7 HYPOTHESIS TESTING
Many times, while analyzing water resources data, questions arise such as: does the river
discharge at the given site follow normal distribution? Is the river water quality violating the
relevant standards? Are the two given variables significantly correlated? In some cases, a
conclusion may be reached without using a statistical test but sometimes the situation could be
such that well-articulated tests are needed to arrive at a conclusion which is not subjective.
Statistical procedures known as hypothesis testing are employed in these situations.
There are two broad categories of hypothesis tests: parametric tests and non-parametric
tests. In parametric tests, the distribution of the data has to be known and is to be specified. The
examples are the t-test, and the F-tests. In the non-parametric or distribution free tests, it is not
required to specify the distribution of data or the data need not follow a specified distribution.
The Kendall’s Tau test and the Kruskal-Wallis test are the examples of these tests.
The parametric tests begin with the formulation of a hypothesis which is termed as null
hypothesis. It is denoted by H0. As the name suggests this is the hypothesis of no change or no
difference. For example, the null hypothesis could be that the flows at two stations follow the
same distribution or there is no correlation between two given variables. This null hypothesis
may be mathematically written as:
H0: 1 = 0 (10.84)
This statement indicates that the mean 1 of a data set is not different from the mean 0 of
another data.
After the null hypothesis has been chosen, an alternative hypothesis is formulated. The
alternate hypothesis is a statement of some expected departure from the null hypothesis. The
null and the alternate hypothesis are mutually exclusive and all inclusive. The alternate
hypothesis can be stated in the following forms
Ha : 1 > 0 (10.85 a)
Ha : 1 0 (10.85 b)
Ha : 1 < 0 (10.85 c)
Equations (10.85) state that the mean of the population from which the sample was drawn
is greater than, not equal to, or less than the specified population mean, respectively. Some
typical ex
xamples of null
n hypotheesis are: the distributionns of flow att two stationns are differeent or
the two variables
v und
der examinattion are relatted to each oother.
β α
Since the rejeection of nulll hypothesiss implies thaat the speciffic statementt is not suppported
by the data,
d the altternate hypo
othesis mustt be sufficieently generaal. If the nuull hypothessis is
rejected, the truth liesl somewhhere in thee wide set of possibiliities stated in the alteernate
hypothessis. Outcomee of hypothesis testing iss stated as reejection of thhe null hypoothesis against the
failure to
o reject the null
n hypothessis.
After
A the hyppothesis has been formu ulated, the teest statistic is computed.. Test statisttic on
which the decision too accept or reject
r null hy
ypothesis is bbased is com mputed from m the sample data.
A rejectiion region iss the set of all test statiistic values for which H0 will be rejected. Thee null
hypothessis will be rejjected if and
d only if the test statistic value falls iin the rejectiion region.
The
T null hyp pothesis may y be correctt or incorrecct and it miight be acceepted or rejeected.
This prodduces four possible comb binations wh hich are indiicated in the Table 10.4.
The
T characterristics of the data that is to be examinned is quanttified by the test statisticc. The
test checks whether the t behaviorr of the test statistic is siimilar to whhat is expectted (leaving aside
the possibility of thiss to happen due
d to chancce alone) if H0 were true..
Table 10
0.4 Hypothesis testing: possible
p outccomes and thheir probabillities.
Hypothesis
H is correct Hypotheesis is incorrrect
Hypothesis is Correct decision. This Type II error with a
accepted outcome has a probability probability
1-
Hypothesis is Type I error with a Correct decision with
rejected probability of probability 1-
In hydrology, commonly the significance level of 0.05 (1 in 20) or 0.01 (1 in 100) are
adopted. A level of 0.05 implies that the decision of statistical test may be in error one time out
of 20. In terms of normal distribution properties, it corresponds to 5% of the area under the
curve. This concept is illustrated for a two-sided test in Fig. 10.13 in which the test statistic under
null hypothesis is normal and each shaded area near the tails contains 2.5% of the total area. This
shaded area is termed as the area of rejection or the critical region. For the inequality type,
alternative hypothesis in eq. (10.85b) the null hypothesis is rejected if the test statistic falls in the
critical region either because it is too high or too low. The choice of significance level depends
upon the risks associated with a wrong decision.
Values of t-distribution for selected degrees of freedom and significance level are given
in Table 10.5 for ready reference. The value of the t-distribution can be read from the row
corresponding to and the column corresponding to the significance level . For example, for
= 10 the t value is 1.81 for 5% significance level. This implies that 95% of the area of the curve
lies to the left of value 1.81. Since the t-distribution is symmetric, 5% of the area in the left tail is
to the left of t value of –1.81 for = 10. For the case where one is interested in 95% of the area
but with 2.5% in each tail, the critical t value for = 2.5 is 2.23 for = 10.
Table 10.5 Values of t-distribution for selected degrees of freedom and significance level
If the test is conducted when population mean 2 is not known, the value of sample standard
deviation is used in the formula in step c).
Examplee 10.14: Asssume that the annual run noff data at a site follow
ws normal diistribution annd 21
observatiions are avaailable. For these observ
vations, the mean xm iss 15.0 mm and the stanndard
deviation
n sd is 5.0 mm.
m Now, whether,
w at 5% significannce level, thhe mean annuual runoff caan be
considereed to be draw
wn from a poopulation whhose mean iss 17.0 mm?
The
T test statisstic is
t = [n (xm – )/sd]0.5 = [21
1*(15.0 – 17.0)/5.0]0.5 = -1.83.
Since this is a two-tailed test, the null hypothesis will be rejected if the test statistic is
either too high or too low and hence value of t is needed for /2 and = n - 1. From Table 10.5,
t/2,n-1 = t2.5,20 = 2.09. Since |t| = 1.83 which is less than 2.09, the statistic does not fall in the
region of rejection and the null hypothesis is accepted.
Z = (X - )/ (10.86)
If these standardised values are squared and added they follow a new statistic:
n
Y Z i2
i 1 (10.87)
The variable Y follows a chi-square (2) distribution with n degrees of freedom. This
distribution is a special case of the gamma distribution and has a single parameter. However, this
distribution is not symmetric and is always positive. The chi-square tests are single-tailed and
the region of rejection is near the right tail. Table 10.6 lists chi-square values for selected degrees
of freedom. For example, for 10 degrees of freedom, 5% of the area in the right tail (region of
rejection) from 2 values is from 18.31 to .
Table 10.6 Values of Chi-square distribution for selected degrees of freedom and significance
level
k
c2 (Oi Ei ) 2 / Ei
(10.88)
i 1
where Oi and Ei are the observed and expected number of observations in the ith segment and k is
the total number of segments. If p parameters are estimated from data, 2c follows a chi-square
distribution with (k – p – 1) degrees of freedom. If the difference between the actual and
expected observations in the segments is large, it implies that samples were not drawn from the
assumed distribution. Therefore, the null hypothesis that the observations follow the assumed
distribution is rejected if 2c > 21-,k-p-1.
Example 10.15: In a goodness of fit test, the data were divided in 11 classes and the value of 2c
came out to be 10.44. If two parameters of the distribution were computed, test whether the
chosen distribution is appropriate for the data at a significance level of 0.1?
Solution: The degree of freedom is 11 – 2 – 1 = 8. From the table of chi-square values, 210,8 =
12.02. Since this value is greater than 10.44, the null hypothesis cannot be rejected. It is,
therefore, concluded that the chosen distribution properly describes the behavior of data.
10.8 Closure
After the measured data of some variable(s) becomes available, the question arises as to what are
key characteristics of the population, what useful information can be derived from the data,
which probability distribution best represents the data. Further, one may also like to statistically
test the validity of a hypothesis. This module deals with these topics.
REFERENCES
Chow, V.T., Maidment, D.R., and Mays, L.W. (1988). Applied Hydrology, McGraw-Hill Book
Company, New York.
Haan, C.T. (2002). Statistical Methods in Hydrology. Iowa State Press, Ames, U.S.A.
Jain, S.K., and Singh, V.P. (2003). Water Resources Systems Planning and Management.
Elsevier, Amsterdam.
Loucks, D.P., and Beek, E.V. (2005). Water Resources Systems Planning and Management- An
Introduction to Methods, Models, and Applications. Studies and Reports in Hydrology,
Unesco, Paris.
McCuen, R.H. (1989). Hydrologic Analysis and Design. Prentice Hall, New Jersey.
Yevjevich, V. (1972). Probability and Statistics in Hydrology. Water Resources Publications,
Fort Collins, Colorado.