"Problem solving in engineering hydrology" is primarily proposed as an addition and a supplementary guide to fundamentals of engineering hydrology. Nevertheless, it can be sourced as a standalone problem solving text in engineering hydrology. The book targets university students and candidates taking first degree courses in any relevant engineering field or related area. The document is valued to have esteemed benefits to postgraduate students and professional engineers and hydrologists. Likewise, it is expected that the book will stimulate problem solving learning and quicken self-teaching. By writing such a script it is hoped that the included worked examples and problems will guarantee that the booklet is a precious asset to student-centered learning. To achieve such objectives immense care was paid to offer solutions toselected problems in a well-defined, clear and discrete layout exercising step-by-step procedure and clarification of the related solution employing vital procedures, methods, approaches, equations, data, figures and calculations.
Discover the world's research
- 20+ million members
- 135+ million publications
- 700k+ research projects
Join for free
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
1
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
2
Problem Solving in
Engineering Hydrology
By
Dr. Eng. Faris Gorashi Faris
&
Prof. Dr. Eng. Isam Mohammed Abdel-Magid Ahmed
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
3
Publisher: University of Dammam Press
© Authors, Dammam, 2013
Revised edition, 2013 (Electronically).
Dr. Faris Gorashi Faris: Deputy Dean, Faculty of Engineering and Technology
Infrastructure, (Infrastructure University Kuala Lumpur, (IUKL), Jalan Ikram –
Uniten, 43000 Kajang , Selangor MALAYSIA, Fax : +603-89256361, Tel: +603-
8738 3388 Ext 312 (off), +6019 640 7203, E.mail: dr.faris@gmail.com,
faris@iukl.edu.my
Prof. Dr. Eng. Isam Mohammed Abdel-Magid Ahmed: Professor of water
resources and environmental engineering, Building 800, Room 240 Environmental
Engineering Department, College of Engineering, University of Dammam, Box
1982, Dammam 31451, KSA, Fax: , Phone: , E-
mail: iahmed@ud.edu.sa, isam_abdelmagid@yahoo.com, Web site:
http://www/sites.google.com/site/isamabdelmagid
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
4
Preface
Engineering hydrology books, manuscripts and text material are abundant addressing
perceptions, fundamentals, methods, models, designs and associated scientific concepts.
However, rare are the books that deal with problem solving and practical complications in the
field of engineering hydrology.
This book has been written to tackle material described in the 2-credit hour course of
"Engineering Hydrology, ENVEN 431" given during the first semester for Junior Engineering
students of the department of environmental engineering of the College Of Engineering at
Dammam University as well as a 3-credit course Hydraulics and Hydrology, BEC 208 offered in
the second year for Bachelor of Civil Engineering students at the Infrastructure University Kuala
Lumpur. The named courses contents encompassed: Hydrologic cycle, precipitation and water
losses. Catchment's characteristics and runoff processes. Flood estimation and control.
Hydrograph analysis. Flood & reservoir routing. Groundwater occurrence, distribution,
movement, exploration and recharge, well hydraulics and design, interaction of ground and
surface water. Differential equations of groundwater flow. Darcy Law, solutions of the steady
and unsteady flow, differential equations for confined and unconfined flows. Pumping test
design. Groundwater models, leaky aquifers. Saltwater intrusion.
The book also covered the course contents of the subject "Hydrology and water gauging" CIVL
561 offered for the fifth year students of the department of civil engineering of Sudan University
for Science and Technology. This particular course covered: Elements of hydrology.
Precipitation. Surface runoff. Evaporation and evapotranspiration. Infiltration and groundwater
flow. Hydrological measurement networks. Hydrograph. Flood routing.
Likewise, the book covered "hydrology" Course contents offered at the civil engineering
department of the college of engineering at United Arab Emirates University.
The book also covered material subject dealt with while teaching the course in "hydrology" at
School of Engineering and Technology Infrastructure, Kuala Lumpur Infrastructure University
College.
Objectives of the book are meant to fulfill the main learning outcomes for students registered in
named courses, which covered the following:
• Solving problems in hydrology and making decisions about hydrologic issues that
involve uncertainty in data, scant/incomplete data, and the variability of natural materials.
• Designing a field experiment to address a hydrologic question.
• Evaluating data collection practices in terms of ethics.
• Interpret basic hydrological processes such as groundwater flow, water quality issues,
water balance and budget at a specific site at local and regional scales based on available
geological maps and data sets.
• Conceptualizing hydrogeology of a particular area in three dimensions and be able to
predict the effects on a system when changes are imposed on it.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
5
Learning outcomes are expected to include the following:
• Overview of essential concepts encountered in hydrological systems.
• Developing a sound understanding of concepts as well as a strong foundation for their
application to real-world, in-the-field problem solving.
• Acquisition of knowledge by learning new concepts, and properties and characteristics of
water.
• Cognitive skills through thinking, problem solving and use of experimental work and
inferences
• Numerical skills through application of knowledge in basic mathematics and supply
issues.
• Student becomes responsible for their own learning through solution of assignments,
laboratory exercises and report writing.
Knowledge to be acquired from the course is expected to incorporate:
• Use real world data to develop a water budget for unfamiliar basins.
• Interpret groundwater resources data and correlate the geology with the groundwater
regime.
• Evaluate ways in which water influences various geologic processes
• Identify interconnections in hydrological systems and predict changes.
• Predications of hydrological terms influencing the hydrological cycle.
• Collection of data, analysis and interpretation
Cognitive skills to be developed are expected to incorporate the following:
• Capturing ability of reasonable scientific judgment and concepts of appropriate decision
making.
• Students will be able to apply the knowledge of hydrology that they have learnt in this
course in practical environmental engineering domain.
• Students should be able to design and apply necessary procedures and precautions to
produce durable hydro logical systems.
The majority of incorporated problems within the book represent examination problems that
were given by the author to students at University of United Arab Emirates (UAE), Sudan
University for Science and Technology (SUST), University of Dammam (UoD) and
Infrastructure University Kuala Lumpur (IUKL). This is to expose students to different and
design constraints and prevailing settings.
"Problem solving in engineering hydrology" is primarily proposed as an addition and a
supplementary guide to fundamentals of engineering hydrology. Nevertheless, it can be sourced
as a standalone problem solving text in engineering hydrology. The book targets university
students and candidates taking first degree courses in any relevant engineering field or related
area. The document is valued to have esteemed benefits to postgraduate students and
professional engineers and hydrologists. Likewise, it is expected that the book will stimulate
problem solving learning and quicken self-teaching. By writing such a script it is hoped that the
included worked examples and problems will guarantee that the booklet is a precious asset to
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
6
student-centered learning. To achieve such objectives immense care was paid to offer solutions to
selected problems in a well-defined, clear and discrete layout exercising step-by-step procedure
and clarification of the related solution employing vital procedures, methods, approaches,
equations, data, figures and calculations.
The authors acknowledge support, inspiration and encouragement from many students,
colleagues, friends, institutions and publishers. The authors salute the motivation and stimulus
help offered to them by Dean Dr. Abdul-Rahman ben Salih Hariri, Dean College of Engineering
of University of Dammam.. Thanks are also extended to Miss Azlinda binti Saadon from the
Faculty of Civil Engineering at The Infrastructure University Kuala Lumpur for her cooperation
and support. Special and sincere vote of thanks would go to Mr. Mugbil ben Abdullah Al-Ruwais
the director of Dammam University Press and his supreme staff for their patience, dedication and
orderly typing of the book
Dr. Eng. Faris Gorashi Faris
Prof. Dr. Eng. Isam Mohammed Abdel-Magid Ahmed
Malaysia, KSA, 2013
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
7
Table of contents
Page
Preface 4
Table of contents 7
List of tables 8
List of figures
List of Appendices
Abbreviations, notations, symbols and terminology used in the book
Chapter (1) Meteorological data (Humidity, temperature, radiation and wind)
1.1 Humidity
1.2 Relative humidity
1.3 Dew point
1.4 Water Budget
1.5 Wind
1.6 Theoretical Exercises
1.7 Problem solving in meteorological data
16
16
16
17
18
19
21
Chapter (2) Precipitation
2.1 Precipitation
2.2 Interpolation of rainfall records
2.3 Methods used to find the height of rainfall
2.4 Theoretical Exercises
2.5 Problem solving in precipitation
25
25
26
31
32
Chapter (3) Evaporation and transpiration
3.1 Evaporation
3.2 Theoretical Exercises
3.3
Problem solving in evaporation
39
41
42
Chapter (4) Infiltration and percolation
4.1 Introduction
4.2 Antecedent precipitation index
4.3 Theoretical Exercises
4.4
Problem solving in Infiltration and percolation
43
47
48
49
Chapter (5) Groundwater flow
5.1 Groundwater
5.2 Theoretical Exercises
5.3
Problem solving in groundwater
52
62
63
Chapter (6) Surface water runoff
6.1 Runoff
6.2 Flow mass curve (Ripple diagram, S-curve)
6.3 Hydrograph
6.4 Theoretical Exercises
6.5
Problem solving in surface runoff
66
67
71
78
79
Chapter (7) General exercises
7.1 Complete missing titles
7.2 True/false sentences
7.3 Underline the best answer
7.4 Rearrange groups
7. Complete missing titles
7. Match the words or phrases 7.
83
83
84
85
86
86
86
References and bibliography for further reading
97
Useful Formulae 99
Appendices 101
About the authors 109
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
8
List of Tables
Table number Description Page
2.1 Forms of precipitation. 25
5.1 Expected well yield. 52
List of Figures
Figure number
Description Page
6.1 Mass curve or Rippl diagram 68
6.2 Divisions of hydrograph 71
List of Appendices
Appendix Description Page
A Nomograph for determining evaporation from free water
surface according to the formula of Penman 102
B Agot's values of short-wave radiation flux RA at the outer
limit of the atmosphere in g cal/cm2/day as a function of the
month of the year & latitude.
103
C Soil descriptions with permeability coefficient. 104
D Saturation vapour pressure as a function of temperature. 105
E Some Physical Properties of Water at Various Temperatures. 106
F Conversion Table 107
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
9
Abbreviations, notations, symbols and
terminology used in the book
A = Area, Area of groundwater basin normal to direction of water flow, Drainage area, Area of
catchment basin (m
2
, km
2
, ha).
A = Cross-sectional area
A
i
= Area of polygon surrounding station number i, located in its middle, m
2
.
a = An exponent varying with surface roughness and stability of the atmosphere.
a = Index of surface connected porosity.
a, b, c = Hydraulic monitoring (gauging) stations to measure rainfall.
a = Rain slope constant.
a = Artificial recharge or loss of irrigation.
a and b = Constants.
a = Height (m) between zero on the gauge and the elevation of zero flow.
B = Heat exchange between soil and surface.
B = Stream roughness factor.
b = An empirical constant.
C = Coefficient representing the ratio of runoff to rainfall.
C = An empirical constant, Constant.
C = Coefficient controlling rate of decrease of loss-rate function.
C = Rational coefficient of surface flow.
C = Weight of salt solution passing sampling point per second.
C = Chezy roughness coefficient
C
2
= Weight of salt passing sampling point per second.
C
b
= Concentration of tracer element in the river at start of injection.
C
i
= Concentration of injected tracer element (within the stream).
C
m
= Concentration of tracer component at measurement point (at equilibrium).
C
p
= Specific heat of air at constant pressure.
C
t
= Coefficient depending on units of drainage basin characteristics.
C
1
= Weight of concentrated solution added per second.
c, n = Locality constants.
c
p
= Constant.
D = Distance from ocean, m.
D = Discharge from stage, m.
D = Rate on outflow, m
3
/s.
D
x
= Time period selected, m/s.
d = Surface-layer depth.
dh/dt = Change in stage during measurement, m/s.
dP = Pressure change with temperature, Pa.
dS/dt = Time rate of change of storage.
E = Amount of evaporation.
E = Evaporation from surface of river basin.
E = Evapotranspiration.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
10
E = Amount of evaporation during a month, mm, inches.
E = Actual vapor pressure, mmHg.
E = Exponent, constant.
E = Actual vapor pressure of air at temp t.
EV = Total evaporation from soil and plants for specified time period.
EV = Evaporation, cm.
E
a
= Evaporation from the earth, and from storage in low-lying areas.
E
a
= Open water evaporation per unit time, mm/day.
E
a
= Evaporation, mm/day.
E
a
= Aerodynamic term.
E
b
= Net energy lost by body of water through exchange of long wave radiation between
atmosphere and body of water.
E
e
= Energy utilized for evaporation.
E
h
= Energy conducted from body of water to atmosphere as sensible heat.
E
o
= Evaporation, mm/day.
E
o
= Evaporation of lake, mm/day.
E
q
= Increase in stored energy in the body of water.
E
r
= Reflected solar radiation.
E
s
= Solar radiation incident to water surface.
E
s
= Saturated vapor pressure.
E
T
= Evaporation from open surface of water (or equivalent in heat energy).
E
T
= Evaporation rate.
E
v
= Net energy adverted into the body of water.
E
w
= Saturation vapor pressure at temp t
w
of surface water of lake, mmHg.
e = Median value for vapor pressure, mbar.
e = Vapor pressure.
e = Base of natural logarithms.
e = Actual vapor pressure.
e = Actual vapor pressure at a defined height above surface.
e = Vapor pressure of air (monthly average), inches Hg.
e = Natural algorithm base.
e
s
= Median value for saturation pressure, mbar.
e
a
= Vapor pressure of air, mbar.
e
s
= Saturation vapor pressure.
e
s
= Saturated vapor pressure (mbar) when temperature is T
w
e
s
= Saturated vapor pressure at surface temperature.
e
s
= Saturation flexibility at surface, mm Hg.
e
s
- e = d, = Lack of saturation Millbar.
e
w
= Partial pressure of the gas for wet bulb temperature.
e
s
= Saturation vapor pressure of air when temperature (t) ° C, mbar, mm Hg.
e
s
= Saturation vapor pressure (monthly average), inches Hg.
e
2
= Vapor flexibility at a height of 2 meters.
F = Leakage.
F = Rate of mass infiltration at time t.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
11
F = Total infiltration.
f (u) = A function in terms of wind speed at a some standard height.
f = Infiltration capacity (percolation) for time t.
f
c
= Constant rate of infiltration after long wetting.
f
o
= Initial infiltration capacity.
f
c
= Final infiltration capacity = apparent saturated conductivity.
f = Natural recharge (rainfall – transpiration, surface runoff & infiltration).
F = Loss rate , mm per hour.
f = Possible head loss.
GI = Growth index of crop, percent of maturity.
H = Ponding depth.
H = Equivalent evaporation of total radiation on surface of plants.
H = Depth of confined aquifer
H = Depth of groundwater basin (m)
H = Saturated thickness of aquifer, m
H = Difference in elevation between two points after subtracting projections, feet.
h = Relative humidity, %
h = Depth below original water level.
h = Stage (gauge height).
I = Surface inflow.
I = Input (volume /time).
I = Part obstructed or trapped from rainfall.
I = Average intensity of rainfall, mm/hr.
I = Rainfall intensity, cm/hour.
I = Inflow to reach.
I = Rate of inflow, m
3
/s.
It = Index value at t days later, mm.
Io = Initial value to the index, mm.
i = Rainfall intensity, in or mm per hour.
i = Groundwater flow through area under consideration.
i = Hydraulic gradient.
K = Dimensionless constant.
K= Coefficient, storage constant (s), slope of relationship of storage-weighted discharge relation.
K
o
= Loss coefficient at start of storm.
K
s
= Effective hydraulic conductivity.
k = A recession constant.
k = Empirical constant for rate of decrease in infiltration capacity.
k = Proportionality factor, coefficient of permeability, or hydraulic conductivity (has dimensions
of velocity).
k = Coefficient of permeability of the aquifer, m/s.
k = Constant of proportionality = reciprocal of slope of storage curve.
L = Latent heat of evaporation.
L = Latent heat per mole of changing state of water.
L = Specific heat of the vapor.
L = Latent heat of evaporation of water.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
12
L = Accumulated loss during the storm, in or mm.
l = Distance in direction of stream line.
L = Depth to wetting front.
L = Maximum distance to entry point.
L = Length of stream, km.
L = Length from remotest point in basin to the outlet.
L = Main stream distance from outlet to divide.
l
ca
= Stream distance from outlet to point opposite to basin centroid.
M
o
= Average annual flow, m
3
/s.
mcs = Water content at saturation that equals porosity (water content at residual air saturation).
mc = Water content at any instant.
N = Normal annual rainfall, mm.
N = Occurrence frequency, usually estimated once every N years = 10/n.
N = Net water infiltration rate resulting from rainfall.
N = Days of the month.
n = Number of days of the month (days of freezing not counted).
n = Number of occurrences in 10 years.
n = Retain coefficient equivalent to coefficient of friction.
n = Number of stations.
n/D = Cloudness ratio.
ne = Effective porosity, dimensionless.
O = Output, volume/time.
o = Groundwater outside borders of region.
P = Total pressure of humid air.
P '= Pressure of dry air.
P = Precipitation during a specific time period, year water.
P = Amount of rainfall, mm.
P = Amount of precipitation (rainfall).
P = Rainfall.
P = Total storm rainfall.
P = Atmospheric pressure, mbar.
P
e
= Net rain representing portion of rainfall that reaches water courses for direct surface runoff.
P '= Decrease in rain.
P
x
= Missing measurement or inaccurate record from station x, mm.
P
av
= Arithmetic mean of rain, mm.
P
i
= Amount of rainfall in station I, mm.
P
mean
= Average rainfall in the region, mm.
P
i
= Record of rainfall in station I, mm.
P ¯ = Average rainfall depth over the area.
P = Point rainfall depth measured at centre of the area.
Q = Total runoff (surface and ground) in basin during allotted time.
Q = Flow rate, flow
Q = Maximum rate of runoff, m3/s.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
13
O = Surface outflow.
Q = Water flow rate.
Q = Total storm (surface) runoff.
Q = Steady state discharge from well
Q = Discharge as a function of time t, which changes from time t
1
to time t
2
Q
i
= Average mean daily flow during a month i.
Q = Average annual flow.
Q = Peak discharge estimation expected to occur after heavy rains in the catchment basin, L/s.
Q = River flow.
Qp = Discharge by probability of P , m
3
/s.
Qa = Actual discharge (measured), m
3
/s
Q = Steady state discharge (discharge from rating curve).
Q
t
= Discharge at end of time t.
Q
a
= Discharge at start of period.
Q = Discharge.
q = Flow rate of injected tracer to stream.
q = Ground water production from wells and disposal channels.
q = Flow in aquifer per unit width of basin, m
3
/s/m.
q
p
= Peak flow.
R = Bowen ratio.
R = Gas constant, L×atmosphere/K.
R = Radiation balance
R
A
= Agot's value of solar radiation arriving at atmosphere.
R
1
and R
2
= Working values, representing indices of storage.
r = Addition and increase of underground water due to egress of surface water.
r = Distance from well.
r
a
= Aerodynamics resistance.
r
H
= Hydraulic radius, m.
r
s
= Net physiological resistance.
S = Slope of energy line.
S = Storage coefficient.
S = Change in storage = Total precipitation over collecting surface.
S = Effective surface retention.
S = Volume of water, storage.
S = Stored volume of surface and ground for specified time period.
S = Change in reserved moisture in river basin.
S = Available storage in surface layer.
S = Absolute slope, m/m.
S = Stream slope, %
S = Steady stage energy gradient at the time of measurement, m/m.
S = Storage, m
3
S
d
= Storage in low-lying areas.
S
f
= Suction (capillary) head at the wetting front.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
14
s = Increase in storage.
s = Drawdown.
T = Absolute temperature, K.
T = Transmissivity, m
2
/day.
T = Number of seconds in a year (86400 = the number of seconds per day).
T = Base time (time from beginning and end of the flood).
T = Tme base for unit hydrograph, day.
T
r
=Residence time.
Tw = Water surface temperature,
o
C.
Ta = Air temperature,
o
C.
t = Average dry temperature of the month, ° C.
t = Dry temperature.
t = Rainfall duration, minutes.
t = Time, day.
t
p
= Basin lag, hours.
t
p
= Time lag in the basin, hours.
t
c
= Time of entry, minutes.
t
r
= Unit rain period (duration of unit hydrograph).
t
pR
= Basin lag for a storm of duration t
R
, hours.
t
w
= Temperature of wet bulb temperature (humid temperature).
(t*) = Gamma inverse function for storm time.
tf = Total time during which rainfall intensity is greater than W.
t
o
= Temperature of evaporation surface.
t
2
= Air temperature at a height of 2 meters.
V = Amount of moisture leaking for the period.
V = Volume of runoff.
V
f
= Final molar volume for states of liquidity and gaseous.
V
i
= Primary molar volume for states of liquidity and gaseous.
v = Flow velocity, m/s.
v = Wind speed, m/w.
v = Velocity of water flow ( = specific velocity).
v' = Average pore velocity (actual or real velocity), m/s.
v = Relative (specific) velocity, m/s.
v = Specific velocity in the horizontal x direction, m.
V
s
= Surface velocity.
V
av
= Average velocity in a section of flow.
U = Groundwater inflows and outflows.
U – O = Quantity of flow in river (surface and underground).
U = Velocity of flood wave (flood wave celerity), m/s
u = Wind speed, miles/hour.
u = Wind speed in the region at height z
o
, m /s.
u
o
= Wind speed at least height (at anemometer) z
o
, , m/s.
u
1
= Wind speed at a height of 1 meter, m/s.
u
2
= Median value of wind speed at a height of 2 m above water surface, m/s.
u
2
= Wind speed at a height of 2 meters, m/s.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
15
u
2
= Wind speed at height of 2 m, m/s.
u
6
= Wind velocity at a height of 6 meters above the surface, m/s.
u (T, t) = Ordinate of the T-hour unit hydrograph derived from those of the t hour unit
hydrograph.
W = W–index = Average infiltration rate during the time rainfall intensity exceeds the capacity
rate.
W = Upper limit of the value of soil humidity that do not move by capillary action from the soil
to the surface, %
W
v
= Volumetric soil moisture, %
X = Dimensionless constant for a certain river reach, dimensionless.
x = Distance along flow line, m.
x, k = Constants derived from observed portion of curve.
y
o
, y
n
= Flows at beginning and end of flood, respectively.
Y
1
, y
3
= Odd flows.
Y
4
, y
6
= Even flows.
z = Height, m.
z
o
= Height of anemometer (minimum height), m
α = Constant.
α = Coefficient of aquifer
γ = Psychrometer constant.
γ = Proportionality factor.
γ = Psychrometer constant fixed device to measure humidity.
∆h = Change in humidity for the period.
∆t = Difference between evaporation surface and air temperature.
ρ = Density of water
θ = Psychometric difference, ° C.
∆ =-Slope of vapor pressure curve at t = tan α
∆s/ ∆ t = Rate of change in reach storage with respect to time.
∆t = Routing Period.
ρ
a
= Air density.
ρ
w
= Density of water.
φ = Head loss over an appropriate base level.
φ = Potential head loss, m.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
16
Chapter One
Meteorological data (Humidity, temperature,
radiation and wind)
1.1 Humidity
Maximum value of water vapor that can exit in any one space is a function of temperature, and is
practically independent of the coexistence of other gases. When maximum amount of water
vapor (for a given temperature) is contained in a given space, the space is said to be saturated.
The pressure exerted by a vapor in a saturated space is called "saturation vapor pressure. "
1.2 Relative humidity
Relative humidity is the percentage of actual vapor pressure to saturation vapor pressure. It is the
ratio of amount of moisture in a given space to the amount a space could contain if saturated.
1.3 Dew point
Dew point denotes the temperature at which space becomes saturated when air is cooled under
constant pressure and with constant water vapor pressure. Also dew point is the temperature
having a saturation vapor pressure, e
s
equal to the existing vapor pressure, e.
Example 1.1
An air mass is at a temperature of 20
o
C with relative humidity of 75%. Using the following
equation find: saturation vapor pressure, actual vapor pressure, the deficit
in saturation and dew point. (B.Sc., DU, 2011)
Solution
1) Data given: T = 20
o
C, h = 75 %
2) Find from tables (see annex) the value of saturated vapor pressure at a temperature of 20
o
C,
Value of saturated vapor pressure e
s
= 17.53 mm Hg
3) Substitute data in humidity equation:
Where:
h = Relative humidity, %, which describes the ability of air to absorb additional moisture
at a given temperature.
e = Actual vapor pressure
e
s
= Saturation vapor pressure. (Numerical value of e
s
are to be found from tables).
h = 100 * e/e
s
, find real vapor pressure:
75 = 100 × e ÷ 17.53
Real vapor pressure e = 13.1475 mm Hg.
es
e
xh 100 =
es
e
xh 100 =
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
17
4) Find saturation deficit as follows:
Saturation deficit: e - e
s
= 17.53 - 13.1475 = 4.3825 mm Hg.
5) Find dew point as the temperature at which the values of e
s
and e are equal.
Since: e
s
= 13.1475 mm Hg, then one can find the temperature, and then the dew point
can be found from tables and for the value of e
s
= 13.1475, dew point = 15.4
o
C.
1.4 Water Budget
Hydrological cycle is a very complex series of processes. Nonetheless, under certain well-
defined conditions, the response of a watershed to rainfall, infiltration, and evaporation can be
calculated if simple assumptions can be made.
Example 1.2
For a given month, a 121 ha lake has 0.43 m
3
/s of inflow, 0.37 m
3
/s of outflow, and the total
storage increase of 1.97 ha-m. A USGS gauge next to the lake recorded a total of 3.3 cm
precipitation for the lake for the month. Assuming that infiltration loss is insignificant for the
lake, determine the evaporation loss, in cm over the lake for the month.
Solution
Given, P = 3.3 cm, I = 0.43 m
3
/s, Q = 0.37 m
3
/s, S = 1.97 ha-m, Area, A = 121 ha
Solving the water balance for Inflow, I and Outflow, O in a lake gives, for evaporation, E ;
Inflow, I =
( )
×
ha m
ha
hour
s
day
hour
month
day
month
s
m
1
00010
121
1
3600
1
24
1
30
143.0
2
3
Inflow, I = 0.92 m = 92 cm
Outflow, Q =
( )
×
ha m
ha
hour
s
day
hour
month
day
month
s
m
1
00010
121
1
3600
1
24
1
30
137.0
2
3
Outflow, Q = 0.79 m = 79 cm
Additional data :
Direct Precipitation, P = 3.3 cm
Storage, S = ha mha
121 .97.1 = 0.0163 m = 1.63cm
From water budget equation,
STEGRP
In this case,
I + P – Q – E = !S
Therefore,
Evaporation, E:
E = I + P – Q - !S
E = 92 + 3.3 – 79 – 1.63
E = 14.67 cm
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
18
EXAMPLE 1.3
At a particular time, the storage in a river reach is 60 x 10
3
m
3
. At that time, the inflow into the
reach is 10 m
3
/s and the outflow is 16 m
3
/s. after two hours, the inflow and the outflow are 18
m
3
/s and 20 m
3
/s respectively. Determine the change in storage during two hours period and the
storage volume after two hours.
Solution
Given, I
1
= 10 m
3
/s
I
2
= 18 m
3
/s
O
1
= 16 m
3
/s
O
2
= 20 m
3
/s
S
1
= 60 x 10
3
m
3
!t = 2 hours x 60 min x 60 s = 7200 s
1) Change in storage, !S during two hours :
I - O = !S
S
OOII
=
−
)
()
(
2121
S
=
−
)
2016
()
1810
(
Therefore,
m
S3
4 −=
∆
s
m
S72004 3 ×−=∆
3
80028 mS −=∆
2) Storage Volume after two hours, S
2
:
! S = S
2
- S
1
Rearrange,
S
2
= !S + S
1
S
2
= - 28 800 + 60 x 10
3
S
2
= 31 200 m
3
1.5 Wind
Wind denotes air flowing nearly horizontally. Winds are mainly the result of horizontal
differences in pressure. In absence of other factors tending to influence wind, it should be
expected that its direction would be from areas of high pressure towards areas of low pressure
and that its speed would vary with the pressure gradient.
Example 1.4
At a given site, a long-term wind speed record is available for measurements at heights of 10m
and 15m above the ground. For certain calculations of evaporation the speed at 2m is required,
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
19
so it is desired to extend the long-term record to the 2m level. For one set of data the speeds at
10 and 15 m were 7.68 and 8.11 m/s respectively (B.Sc., DU, 2013).
i. What is the value of the exponent relating the two speeds and elevations?
ii. What speed would you predict for the 2 m level?
Solution
• Data given: z
1
= 10 m, z
2
= 15, u
10
= 7.68 m/s, u
15
= 8.22 m/s. Required: a, u
2
• Use equation to determine the constant a:
Where:
u = Wind speed in the region at height zo m /s
uo = Wind speed at least height (at anemometer) zo . m/s
zo = Height of anemometer (minimum height), m
a = An exponent varying with surface roughness and stability of the atmosphere, (usually
ranging between 0.1 to 0.6 in the surface boundary layer).
Or: a*Log (10/15) = Log (7.68/8.11) which yields a = 0.134
• Use equation
Which yields: u
2
= 6.19 m/s
1.6 Theoretical Exercises
1) Write briefly about three of the following: (B.Sc., UAE, 1989).
a) The hydrologic cycle and the principal influences that prevent its phases from repeating
themselves in an identical fashion from year to year.
b) Factors that influence evaporation. Explain briefly how each factor affects the rate and
time of occurrence of evaporation.
c) Types of precipitation. On what category would you place rain rainfall at El-Ain.
d) Infiltration indexes and their importance.
2) Show how the following relationship can be used to find the wind speed from the height
above the ground. (B.Sc., DU, 2011)
3) State Buys Ballot law (B.Sc., DU, 2013).
Solution
"Low-pressure zone in the northern part of the globe lies north of the viewer who stands
giving his back to the wind, and it is located on his right side in the southern part of the
globe, as a result of the impact of the earth's rotation".
4) It is necessary to specify height above sea level when doing any measurement of wind due
to surface friction factors and water surfaces through which wind is blowing. The
a
z z
u
u
oo
=
a
=15
10
11.8 68.7
a
z
z
u
u
=
2
10
2
10
134.0
2
1068.7
2
=
u
a
z z
u
u
oo
=
a
z
z
u
u
oo
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
20
relationship between wind speed and height can be found from the following relationship.
Define the shown terms. (B.Sc., DU, 2012)
5) Define and explain: Dew point, Ballton's law, Orographic precipitation and an isohyetal
map. (B.Sc., UAE, 1989, B.Sc., DU, 2012)).
Solution
• Is the temperature at which e
s
= e.
• If an observer stands with his back to the wind, the lower pressure is on his left in
the northern hemisphere and his right in the southern hemisphere.
• Results from mechanical lifting over mountain barriers. Most is deposited on the
windward slopes.
• Map showing lines of equal rainfall amount.
6) Define: relative humidity, dew point and wind (B.Sc., DU, 2012).
Solution
Relative humidity is the percentage of actual vapor pressure to saturation vapor pressure or
is the ratio of amount of moisture in a given space to the amount a space could contain if
saturated.
Dew point: Denotes the temperature at which space becomes saturated when air is cooled
under constant pressure and with constant water vapor pressure. Also it is the temperature
having a saturation vapor pressure, e
s
equal to the existing vapor pressure, e.
Wind: denotes air flowing nearly horizontally.
7) What is hydrology? What are the benefits of this science in practical life?
8) Indicate general factors affecting climate, and methods of measurement. Illustrate your
answer with appropriate sketches.
9) Mention advantages of measuring temperature.
10) What is meant by "areas of low and high pressure"? How are they detected?.
11) Write a detailed report on each of the following: moisture, relative humidity, water vapor
pressure, and dew point.
12) Outline marked differences between methods of measuring humidity.
13) What is the benefit of measuring solar radiation?
14) Explain ways of cloud formation.
15) Explain water vapor condensation in air. What are the benefits of the process?
16) What factors affect evaporation?
17) Write a detailed report on instruments for measuring evaporation.
18) How does both surface water and groundwater affect evaporation?
19) Mention different types of precipitation and their presence in practice.
20) How precipitation is measured? Compare its measurement devices.
21) What are elements of a hydrological Meteorological station?
a
o
z z
o
u
u
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
21
1.7 Problem solving in meteorological data
Pressure
1) A hill is of height of 4000 meters and of rising dry air temperature of 14.8 degrees Celsius.
Condensation occurs at a height of 2000 meters. Calculate the pressure and saturation drop.
Find air temperature when it reaches the top of the mountain and its temperature when it
reaches top of hill and its temperature when it returns to the bottom.
Vapor pressure
2) An air mass is at a temperature of 18
o
C with relative humidity of 80 %. Using the following
equation find: (B.Sc., DU, 2012)
i. Saturation vapor pressure.
ii. Actual vapor pressure in mm Hg, m bar, and Pa
iii. The deficit in saturation
iv. Dew point.
Solution
1) data: T = 18
o
C, h = 80 %
2) Find from tables the value of saturated vapor pressure at a temperature of 18
o
C,
value of saturated vapor pressure e
s
= 15.46 mm Hg
3) Substitute data in h = 100 * e / e
s
find real vapor pressure:
80 = 100 × e ÷ 15.46
Actual or real vapor pressure e = 12.386 mm Hg = 12.386*1.33 = 16.45 m bar =
16.45*10
2
= 1.6*10
3
N/m
2
= 1.6 kPa
4) Find saturation deficit as follows:
Saturation deficit: e - e
s
= 15.46 - 12.386 = 3.07 mm Hg.
5) Find dew point as the temperature at which the values of e
s
and e are equal.
Since: e
s
= 12.386 mm Hg, then one can find the temperature, and then the dew point
can be found from tables and for the value of e
s
= 12.386 dew point = 13.2
o
C.
3) An air mass is at a temperature of 24.5 ºC with relative humidity of 58%. Determine: (B.Sc.,
UAE, 1989).
Saturation vapour pressure.
Saturation deficit.
Actual vapour pressure (in mbar, mmHg, and Pa).
Dew point.
Solution
For a temperature of 24.5 ºC from table or figure
i] Saturation pressure e
s
= 23.05 mmHg
ii] Relative humidity, h =
Saturation deficit = e
s
– e = 23.05 – 13.37 = 9.68 mmHg
iii] Actual vapour pressure = 13.37 mmHg = 13.37×1.33 = 1778 mbar
= 1778×10
2
= 1.78 kN/m
2
= 1.78 kPa.
iv] Dew point is the temperature at which e = e
s
∴ e
s
= 13.37 mmHg and from table:
e
e
xh 100 =
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
22
Dew point = 15.7 ºC
4) A mass of air is at a temperature of 22°C and relative humidity of 81%. Find: saturated
vapor pressure, real vapor pressure, saturation deficit, and dew point (Ans. "#, 16.05,
" mmHg, "°C).
5) An air mass is at a temperature of 28
o
C with relative humidity of 70%. Find:
a) Saturation vapor pressure
b) Actual vapor pressure in mbar and mm Hg
c) Saturation deficit
d) Dew point
e) Wet-bulb temperature.
6) To a mass of air at a temperature of 20
o
C and relative humidity of 80%, find saturated vapor
pressure, real vapor pressure, deficit in saturation, and dew point. What is the difference in
your answer if the temperature changed to 25
o
C and relative humidity to 75%.
7) An air mass is at a temperature of 20
o
C with relative humidity of 75%. Using the following
equation find:
i. Saturation vapor pressure
ii. Actual vapor pressure
iii. The deficit in saturation
iv. Dew point
(ans. 22.53 mm Hg, 16.8975 mm Hg, 5.63 mm Hg, 19.4
o
C).
Wind
8) Anemometers at 2.5 m and 40 m on a tower record wind speeds of 2 and 5 m/s respectively.
Compute wind speeds at 5 and 30 m. (B.Sc., DU, 2012)
Solution
b)
9) At a given site, a long-term wind speed record is available for measurements at heights of 5
m and 10 m above the ground. For certain calculations of evaporation the speed at 2 m is
required, so it is desired to extend the long-term record to the 2 m level. For one set of data
the speeds at 5 and 10 m were 5.51 and 6.11 m/s respectively (B.Sc., DU, 2012).
i. What is the value of the exponent relating the two speeds and elevations?
ii. What wind speed would you predict for the 2 m level?
Solution
1. Data given: z
1
= 5 m, z
2
= 10, u
5
= 5.51 m/s, u
10
= 6.11 m/s r.
2. Required: a, u
2
e
e
xh 100 =
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
23
3. Use equation to determine the constant a:
Or a Log 10/15 = Log 5.51/6.11 which yields a = 0.149
4. Use equation
Which yields u
2
= 4.8 m/s
10) At a given site, a long-term wind speed record is available for measurements at heights of 5
m and 10 m above the ground. For certain calculations of evaporation the speed at 2 m is
required, so it is desired to extend the long-term record to the 2 m level. For one set of data
the speeds at 5 and 10 m were 5.11 and 5.61 m/s respectively. (B.Sc., DU, 2011)
iii. What is the value of the exponent relating the two speeds and elevations?
iv. What speed would you predict for the 2 m level?
Solution
1. Data given: z
1
= 5 m, z
2
= 10, u
10
= 5.11 m/s, u
15
= 5.61 m/s r.
2. Required: a, u
2
3. Use equation to determine the constant a:
U
10
/u
15
= (5/10)
a
= 5.11/5.61
Or a Log 5/10 = Log 5.11/5.61 which yields a = 0.134677
4. Use equation u
10
/u
2
= 5.11/u
2
= (5/2)
0.134667
,
Which yields U
2
= 4.51 m/s
11) Anemometers at 2m and 50m on a tower record wind speeds of 2 and 5 m/s respectively.
Compute wind speeds at 5 and 30 m. (B.Sc., UAE, 1989).
Solution
12) Anemometers at 2.5 m and 40 m on a tower record wind speeds of 2 and 5 m/s respectively.
Compute wind speeds at 5 and 30 m. (ans. 2.51, 4.54 m/s).
13) Anemometers at 2m and 50m on a tower record wind speeds of 2 and 5 m/s respectively.
Compute wind speeds at 5 and 30 m. (ans. 2.6, 4.6 m/s).
14) Theoretical Wind speed is measured for two heights 3 and 4 meters and the following values
were found 2.5 and 3 m/s, respectively. Find wind speed for a height of two meters.
15) Measured wind speed in a monitoring station for two heights 1 and 2 meters has shown the
values of 1.6 and 2.2 m/s, respectively. Find wind speed at a height of two meters and ten
meters.
a
z
z
u
u
=
10
5
10
5
a
z
u
u
=2
2
55 149.0
2
5
2
51.5
=
u
a
=10
5
11.6 51.5
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
24
Water budget
16) At a water elevation of 6391 ft, Kenyir Lake has a volume of 2, 939, 000 ac-ft, and a surface
area, A of 48, 100 ac. Annual inputs to the lake include 8.0 in. of direct precipitation, runoff
from gauged streams of 150, 000 ac-ft per year, and ungauged runoff and groundwater
inflow of 37, 000 ac-ft per year. Evaporation is 45 inch. per year.
Make a water budget showing inputs, in ac-ft per year.
(*Note : 1 ac = 43, 560 ft
2
, 1 cm = 0.394 in, 1 in = 0.083 ft)
Solution
Given, Water Elevation, EL = 6391 ft
Volume of lake, V = 2, 939, 000 ac-ft
Area, A = 48, 100 ac
1) Annual input, I to the lake :
a) Direct precipitation, P = 8 in / year
=
( )
ac
in ft
in 100,48
1
083.0
)8(
= 31, 938.40 ac-ft / year
b) Runoff from gauged streams = 150, 000 ac-ft / year
c) Runoff from ungauged = 37, 000 ac-ft / year
Therefore,
Total annual input, I = 31, 938.40 + 150, 000 + 37, 000
= 218, 938.40 ac-ft / year.
2) Outflow, O :
a) Evaporation, E = 45 in / year
E =
( )
ac
in ft
in 100,48
1
083.0
)45(
E = 179, 653.50 ac-ft / year
3) Water Budget :
! S = I - O
! S = 218, 938.40 - 179, 653.50
! S = 39, 284.90 ac-ft / year (positive)
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
25
Chapter Two
Precipitation
2.1 Precipitation
Precipitation denotes all kinds of rainfall on surface of earth from vapor in the atmosphere (most
of humidity at a height of 8 km from the Earth's surface).
Table 2.1 Forms of precipitation [1,2, 10]
2.2 Interpolation of rainfall records
Sometimes records may be lost from the measuring or monitoring station for a specific day or
several days because of the absence of station operator (observer) or because of instrumental
failure or malfunction or damage in the recording devices, for any other reason. In order not to
lose information, it is best to use an appropriate way to estimate the amount of rain in these days
in calculating monthly and annual totals. Procedure for these estimates are based depending on
simultaneous records for three stations close to and as evenly spaced around the station with
missing records as possible. This station should be equi-distant from the three stations and the
following conditions should be achieved:
• If the normal annual precipitation at each of the these stations is within ten percent of that
of the station with missing records, a simple arithmetic average of the precipitation at the
three stations is used for estimating missing record of the station.
• If the normal annual precipitation at any one of the three stations differs from that of the
station with missing records by more than ten percent, the normal ratio method is used
[15].
Example 2.1
1) Rain gauge X was out of operation for a month during which there was a storm. The rainfall
amounts at three adjacent stations A, B, and C were 37, 42 and 49 mm. The average annual
precipitation amounts for the gauges are X = 694, A = 726, B = 752 and C = 760 mm.
Using the Arithmetic method, estimate the amount of rainfall for gauge X.
(mm)
(mm)
Rate of fall (mm/hr) Form of precipitation 2.5
2.8 – 7.6
> 7.6
Light
moderate
Heavy
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
26
Solution
If
N
x
= 694
Then 10% from 694 =
Therefore, precipitation allowed = (694 – 69.4) ~ (694 + 69.4) mm
= 624.6 mm ~ 763.4 mm
Since all annual precipitations (726, 752 and 760) mm are within the ranges, Arithmetic
Method can be applied :
{ }
mmP
x
7.42494237
1=++=
Example 2.2
One of four monthly-read rain gauges on a catchment area develops a fault in a month when the
other three gauges record 48, 58 and 69 mm respectively. If the average annual precipitation
amounts of these three gauges are 741, 769 and 855 mm respectively and of the broken gauge
707 mm, estimate the missing monthly precipitation at the latter (B.Sc., DU, 2013).
Solution
• Data: P
a
= 48 mm, P
b
= 58 mm, P
c
= 69 mm, N
a
= 741 mm, N
b
= 769 mm, N
c
= 855 mm,
N
x
= 707 mm.
• Find the value of rainfall during the storm at station (a) by using the following equation:
Example 2.3
The records of precipitation of hydraulic monitoring stations (x) in a rainy day are missing. The
data indicate that the estimates of rainfall at three stations (b), (c) and (d) adjacent to the station
(x) are equal to: 80, 70 and 60 mm, respectively. If the average annual rainfall at stations (a) and
(b) and (c) and (d) is: 650, 240, 320 and 140 mm, respectively, find the value of rainfall during
the rain storm in station (x) (B.Sc., DU, 2012) .
Solution
1. Data: Pb = 80 mm, Pc = 70 mm, Pd = 60 mm, Nx = 650 mm, Nb = 240 mm, Nc = 320
mm, Nd = 140 mm.
2. Find the value of rainfall during the storm at station (a) by using the following equation:
Px = 212.5 mm
2.3) Methods used to find the height of rainfall
Methods used to find the height of rainfall Include: Arithmetic mean method, Thiessen polygon
method and isohyets method.
NP
N
NP
N
NP
N
Pc
cx
b
bx
a
ax
x
++
=
mm
xxx
NP
N
NP
N
NP
N
P
c
cx
b
bx
a
ax
x
52
3855
69707
769
58707
741
48707
3
69
===
++
++
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
27
Example 2.4
Using the data given below, estimate the average precipitation using Thiessen Polygon
method.
Stations Area (km
) Precipitation (mm) Area x Precipitation (km
2
mm)
A 72 90 6 480
B 34 110 3 740
C 76 105 7 980
D 40 150 6 000
E 76 160 12 160
F 92 140 12 880
G 46 130 5 980
H 40 135 5 400
I 86 95 8 170
J 6 70 420
TOTAL 568 1 185 69 210
Solution
AreaPolygonTotal necipitatioStationeachforAreaPolygon
ionPercipitatAverage ∑
=Pr
69210
= ionPercipitatAverage
mmionPercipitatAverage 8.121
Example 2.5
The precipitation on a catchment in Dubai of area 95 km
2
is sampled in table (1). Determine the
precipitation recorded by station number (7) if the mean precipitation, as computed by Thiessen
method, amounts to 98mm. (B.Sc., UAE, 1989).
Rain gauge Recorded rainfall,
Feb. 1968 (mm) Thiessen polygon on area (km
)
1 84 4.0
2 90 4.0
3 120 10.0
4 86 5.1
5 87 15.1
6 76 30.6
7 X 6.2
8 131 20.0
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
28
Solution
a)
Rain gauge Recorded rainfall,
Feb. 1968 (mm) Thiessen polygon on
area (km
2
) Product
1 84 4.0 336
2 90 4.0 360
3 120 10.0 1200
4 86 5.1 438.6
5 87 15.1 1313.7
6 76 30.6 2325.6
7 X 6.2 6.2x
8 131 20.0 2620
Example 2.6
Find area precipitation by isohyetal method for a certain catchment area given the following
data (B.Sc., DU, 2012).
Isohyet, in
Area enclosed within basin boundary, sq. mile
6.8
6 20
5 97
4 213
3 410
2 602
1.5 633
Solution
Isohyet,
in
Col. (1)
Area enclosed within
basin boundary, sq. mile
Col. (2)
Net area, sq. mile
Col. (3)
(b
n
-b
n-1
)
Average
precipitation
Col. (4)
(a1+a2)/2
Precipitation
volume
Col. (3)*Col.
(4)
6.8
6 20 20 (6.8+6)/2 = 6.4 128
5 97 97-20 = 77 5.5 423.5
4 213 213-97 = 116 4.5 522
3 410 410-213 = 197 3.5 689.5
2 602 602-410 =192 2.5 480
1.5 633 633-602 = 31 1.75 54.25
Total area = 633 Total volume
= 2297.25
P = 2297.25/633 = 3.63 in
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
29
Example 2.7
The following table shows rainfall observations taken for a certain area. Estimate the mean
basin precipitation (B.Sc., DU, 2012) .
Table: Rainfall observations.
(in)
within basin boundary
(ha)
precipitation between
isohyets
> 6 50 6.3
> 5 120 5.5
> 4 250 4.7
> 3 450 3.6
> 2 780 2.7
> 1 999 1.4
< 1 1020 0.8
Solution
Determine corresponding net area as in table below:
Isohyets
(in) Total area enclosed
within basin
boundary (ha)
Corresponding
net area
Estimated mean
precipitation between
isohyets
> 6 50 50 6.3
> 5 120 70 5.5
> 4 250 130 4.7
> 3 450 200 3.6
> 2 780 330 2.7
> 1 999 219 1.4
< 1 1020 21 0.8
= 11.06 "
Example 2.8
The following rainfall observations were taken for a certain area: (B.Sc., UAE, 1989).
Isohyets (in) Total area enclosed within
basin boundary (ha) Estimated mean precipitation
between isohyets
> 6 40 6.2
> 5 110 5.5
> 4 236 4.5
> 3 430 3.5
> 2 772 2.5
> 1 990 1.5
< 1 1010 0.9
Estimate the mean basin precipitation
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
30
Solution
Determine corresponding net area as presented in the following table:
Isohyets
(in) Total area enclosed
within basin
boundary (ha)
Corresponding
net area
Estimated mean
precipitation between
isohyets
> 6 40 40 6.2
> 5 110 70 5.5
> 4 236 126 4.5
> 3 430 194 3.5
> 2 772 342 2.5
> 1 990 218 1.5
< 1 1010 20 0.9
Example 2.9
Use the isohyetal method to determine the average precipitation depth within the basin for the
storm.
0
P1 =20mm
P2 =30mm P 3 =40mm
P4 =50mm
P5 =60mm
P6 =70mm
A1
A3
A4
A5
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
31
Solution
Isohyetal
interval Average
Precipitation
(mm)
Area
(km
2.
)
Area x
Precipitation
(km
2.
mm)
< 10 10 0 0
10 – 20 15 84 1 260
20 – 30 25 75 1 875
30 – 40 35 68 2 380
40 – 50 45 60 2 700
50 – 60 55 55 3 025
60 – 70 65 86 5 590
TOTAL 428 16 830
Area necipitatioArea
ionPercipitatAverage ∑
=Pr
16830
= ionPercipitatAverage
mmionPercipitatAverage 3.39
2.4 Theoretical Exercises
1) Distinguish between the different forms of precipitation.
Solution
Snow: Complex ice crystals. A snowflake consists of ice crystals.
Hail: Balls of ice that are about 5 to over 125 mm in diameter. Hailstones have the potential
for agricultural & other property damage.
Sleet: Results from freezing of raindrops and is usually a combination of snow and rain.
Rain: Consist of liquid water drops of a size 0.5 mm to about 7 mm in diameter.
Drizzle: Very small, numerous and uniformly dispersed water drops that appear to float
while following air currents. Drizzle drops are considered to be less than 0.5 mm diameter. It
is also known as warm precipitation.
2) Outline sources of error when recording readings and record-keeping by precipitation gauge
(B.Sc., DU, 2011)
Solution
• Mistakes in reading scales of gauge,
• Displacement of some collected water during recording of reading (stick creep),
• Rainfall splash from collector,
• Loss of water needed to moisten (wetting) funnel and inside surface,
• Any change in the receiving area of precipitation, and inclination of gauge (e.g. dents in
collector rim)
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
32
• Deficiency of instruments due to wind or otherwise, or the absence of regular
maintenance.
3) List five natural and man-made factors that may influence the distribution of rainfall in a
given area. (B.Sc., UAE, 1989).
4) Indicate how the following equation illustrates the method of weighting, by the ratio of the
normal annual precipitation values to estimate missing records for a neighboring
monitoring station. (B.Sc., DU, 2011)
5) Define the terms used in the equation that illustrates the method of weighting to estimate
missing records for a neighboring monitoring station (B.Sc., DU, 2012).
6) What are the expected sources of error when recording precipitation readings in
hydrological meteorological stations?
7) Explain methods of rain calculation? Indicate advantages and shortcomings of each method.
8) How do you assess rainfall data lost in a station as compared to data allocated to the
neighboring stations?
9) What is the relationship between intensity, duration and frequency of rainfall?
2.5 Problem solving in precipitation
Arithmetic average rainfall
1) Find average rainfall for the rainfall data observed in the hydrological measuring stations
described in the following table
2) Average rainfall in a specific region is equal to 310 mm according to the data shown in the
table below for estimated rainfall in four hydrological monitoring stations: What is the
average rainfall in station number (4) using arithmetic average method for estimating
rainfall? Amount of rain, mm Station number 1 2 # 3 $ 4
Thiesen polygon
3) Find average precipitation in watershed according to the information and data that are
described in the following table using average arithmetic method and Thiesen polygon .
(Ans. ", ).
3NP
N
NP
N
NP
N
Pc
c
x
b
b
x
a
a
x
x
++
=
3NP
N
NP
N
NP
N
Pc
c
x
b
b
x
a
a
x
x
++
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
33
Area, % Rain gauge,
cm
Station
8
11
20
18
33
10
62
74
87
112
118
107
1
2
3
4
5
6
4) List five natural and man-made factors that may influence the distribution of rainfall in a
given area. The precipitation on a catchment in Dubai of area 95 km
2
is sampled in table (1).
Determine the precipitation recorded by station number (7) if the mean precipitation, as
computed by Thiessen method, amounts to 98mm. (ans. 115.5 mm).
Rain gauge Recorded rainfall,
Feb. 1968 (mm) Thiessen polygon on area
(km
2
)
1 84 4.0
2 90 4.0
3 120 10.0
4 86 5.1
5 87 15.1
6 76 30.6
7 X 6.2
8 131 20.0
5) Find precipitation by Thiessen method for the observed precipitation records shown in table
(4) (B.Sc., DU, 2013).
Table: Precipitation records.
Observed precipitation, in Area of corresponding polygon within
basin boundary, sq. mile
1.59 12
2.36 112
2.84 109
3.61 122
2.46 19
3.9 84
6.11 94
5.61 68
∑
=
∑
=
=n
iA i
n
iP i
Ai
P
mean
1
1
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
34
Solution Observed
precipitation,
in
(1)
Area of corresponding
polygon within basin
boundary, sq. mile
(2)
Percent total
area
(3)
Weighted
precipitation,
in
Col. (1) x Col.
(3)
1.59 12 1.9 0.03
2.36 112 18.1 0.43
2.84 109 17.6 0.5
3.61 122 19.7 0.71
2.46 19 3.1 0.08
3.9 84 13.5 0.53
6.11 94 15.1 0.93
5.61 68 11 0.62
Total = 620 100 3.83
Average P = 3.83 in
Method of isohyets
6) Calculate the average rainfall of the data set in the following table using isohyet method
rainfall.
Area between isohyets, km
Isohyet, cm
10
#
#
#
7) Find average rainfall for the following data using isohytal method,
Area between isohyets, km
Isohyet, cm
# #
8) Find average precipitation in catchment area according to the information and data that are
described in the following table using the method of isohyets .
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
35
Estimated average
rain between
isohyets, inch
total area bounded
by boundaries of the
basin (ha)
Isohyets, inch
6.5
5.5
4.5
3.5
2.5
1.5
0.9
#
#
6 <
5 <
4 <
3 <
2 <
1 <
1 >
9) Table shows rainfall observations taken for a certain area. Estimate the mean basin
precipitation. (ans. 11.06").
Isohyets (in) Total area enclosed within
basin boundary (ha) Estimated mean precipitation
between isohyets
> 6 50 6.3
> 5 120 5.5
> 4 250 4.7
> 3 450 3.6
> 2 780 2.7
> 1 999 1.4
< 1 1020 0.8
10) The following rainfall observations were taken for a certain area:
Isohyets (in) Total area enclosed within
basin boundary (ha) Estimated mean
precipitation between
isohyets
> 6 40 6.2
> 5 110 5.5
> 4 236 4.5
> 3 430 3.5
> 2 772 2.5
> 1 990 1.5
< 1 1010 0.9
Estimate the mean basin precipitation (ans. 3').
Rainfall intensity
11) Find average rainfall intensity over an area of 6 square kilometers during a one-hour storm
of frequency of once every ten years (assuming point rain of 30 mm, and inverse gamma
5.6) (Ans. 26 mm).
12) What is the average rainfall intensity over an area of 5 square kilometers during a 60 minute
storm with a frequency of once in 10 years (take p = 25 mm).
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
36
13) From the precipitation data given, estimate cumulative rainfall and rainfall intensity.
Time (min) 0
10
20
30
40
50
60
70
80
90
Rainfall (cm) 0
0.18
0.21
0.26
0.32
0.37
0.43
0.64
1.14
3.18
Time (min) 100
110
120
130
140
150
160
170
180
Rainfall (cm) 1.65
0.81
0.52
0.26
0.42
0.36
0.28
0.19
0.17
Solution
Time
(min) Rainfall
(cm) Cumulative
Rainfall
(cm)
Rainfall
Intensity
(cm/hr)
0 0 0
10 0.18 0.18 1.08
20 0.21 0.39 1.26
30 0.26 0.65 1.56
40 0.32 0.97 1.92
50 0.37 1.34 2.22
60 0.43 1.77 2.58
70 0.64 2.41 3.84
80 1.14 3.55 6.84
90 3.18 6.73 19.08
100 1.65 8.38 9.90
110 0.81 9.19 4.86
120 0.52 9.71 3.12
130 0.42 10.13 2.52
140 0.36 10.49 2.16
150 0.28 10.77 1.68
160 0.24 11.01 1.44
170 0.19 11.20 1.14
180 0.17 11.37 1.02
API
14) The value of the API for a certain station reached 42 mm in first of August, and an amount
of 46 mm of rain fell on the fifth of August, also rain of 28 mm fell in the seventh of august
and 34 mm on the eighth day of August. Find value of the API guide on 12th August noting
that k is equal to 0.92. Calculate the amount of the index assuming no rain in the same
period. Draw curve of change of index with time. (Ans. 85.3, 1.7 mm).
15) The value of the API for a station reached 55 mm in the first of August, an amount of 57
mm of rain fell on the sixth of August, 31 mm of rain fell in the eighth of August and 20 mm
fell in the ninth day of August. Find value of API index for the 15th day of August, noting
that k is equal to 0.92. Calculate the amount of index assuming no rain during the same
period. Draw curve of change of index with time.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
37
16) The antecedent precipitation index (API) is given by: I
e
= I
o
k
t
, Define the parameters shown
in the equation. The API for a station was 60 mm on the first of December, 64 mm rain fell
on the seventh of December, 40 mm rain fell on the ninth of December, and 30 mm rain fell
on the eleventh of December. Compute the API for 20
th
December if k = 0.85. Determine the
API if no rain fell. (ans. 24.1, 2.7 mm).
Missing record
17) A precipitation station (X) was inoperative for some time during which a storm occurred.
The storm totals at three stations (A), (B) and (C) surrounding X, were respectively 6.60,
4.80 and 3.30 cm. The normal annual precipitation amounts at stations X, A, B and C are
respectively 65.6, 72.6, 51.8 and 38.2 cm. Estimate storm precipitation for station (X).
18) Rainfall data taken from the records of hydrological monitoring stations (a), (b), and (c), and
(d) gave values of 17, 18, 24 and 18 mm respectively. Thiesen polygon method was chosen
to calculate the average rainfall. The area of each of the polygons surrounding each station is
as follows: 24, 18, 32 and 35 square kilometers for stations (a), (b), and (c), (d), respectively.
Find average amount of rainfall in the region.
19) Rain record is missed from a hydraulic station (m) in a windy day. However, the estimates
of rainfall at three stations (m1), (m2) and (m3) surrounding station (m) is equal to 25, 50
and 35 mm respectively. Note that the average annual precipitation at the stations (m), (m1),
(m 2) and (m3) equals 500 and 700 and 590 and 440 mm, respectively. Find the value of
rainfall during the storm at station (m).
20) Four hydrological monitoring stations are of similar conditions registered the data shown in
the following table. Estimate missing recod.
Station 1 2 3 4
Mean annual precipitation,
mm 350 1200 1100 1000
Present precipitation, mm 750 ? 760 810
21) Outline sources of error when recording readings and record-keeping by precipitation gauge.
Indicate how the following equation illustrates the method of weighting, by the ratio of the
normal annual precipitation values to estimate missing records for a neighboring monitoring
station.
The records of precipitation of hydraulic monitoring stations (x) in a rainy day are
missing. The data indicate that the estimates of rainfall at three stations (b), (c) and
(d) adjacent to the station (x) are equal to: 70, 608 and 50 mm, respectively. If the
average annual rainfall at stations (a) and (b) and (c) and (d) is: 600, 250, 510 and
120 mm, respectively, find the value of rainfall during the rain storm in station (x).
(ans. 162.9 mm).
22) The average annual precipitation amounts for the gauges A, B, C and D are 1120, 935, 1200
and 978 mm. In year 1975, station D was out of operation. Station A, B, and C recorded
rainfall amounts of 107, 89 and 122 mm respectively. Estimate the amount of precipitation
for station D in year 1975.
3NP
N
NP
N
NP
N
Pc
c
x
b
b
x
a
a
x
x
++
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
38
Stations
Normal Annual
Precipitation (mm) Amounts of Precipitation
Year 1975 (mm)
A 1120 107
B 935 89
C 1200 122
D 978 X
Solution
If
N
x
= 978
Then 10% from 694 =
Therefore, precipitation allowed = (978 – 97.8) ~ (978 + 97.8) mm
= 880.2 mm ~ 1075.8 mm
Since the average annual precipitations amounts for the gauges A and C exceeded 1075.8 mm,
therefore Normal Ratio Method is used :-
mmP x3.95
1200
122
935
89
1120
107
3
978 =
++=
23) The records of precipitation of hydraulic monitoring stations (x) in a rainy day are missing.
The data indicate that the estimates of rainfall at three stations (b), (c) and (d) adjacent to the
station (x) are equal to: 70, 608 and 50 mm, respectively. If the average annual rainfall at
stations (a) and (b) and (c) and (d) is: 600, 250, 510 and 120 mm, respectively, find the value
of rainfall during the rain storm in station (x). (B.Sc., DU, 2011)
Solution
1. Data: Pb = 70 mm, Pc = 60 mm, Pd = 50 mm, Nx = 600 mm, Nb = 250 mm, Nc = 510
mm, Nd = 120 mm.
2. Find the value of rainfall during the storm at station (a) by using the following equation:
Pa = (1 ÷ 3) × {(600 × 70 ÷ 250) + (600 × 60 ÷ 510) + (600 × 50 ÷ 120)} = 162.9 mm
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
39
Chapter Three
Evaporation and transpiration
3.1 Evaporation
Evaporation is the process by which water is changed from the liquid or solid state into the
gaseous state through the transfer of heat energy [15]. Evapotranspiration: the process by which
water is evaporated from wet surfaces and transpired by plants [15].
Methods of estimating evaporation or transpiration include: Evaporation pans, empirical
formulas, methods of water budget, methods of mass transfer and methods of energy budget
[13].
Example 3.1
The daily potential evapotranspiration from a field crop at latitude 20° N in November has been
predicted through use of the nomogram for the solution of Penman's equation to be equal to 2.8
mm. Estimate the mean wind speed at height of 2m, u
2
, m/s. Governing conditions are as
outlined in the following table (B.Sc., DU, 2013).
Table: Evaporation data.
Item Value
Mean air temperature, T,
C 22
Sky cover 70 % cloud
mean humidity, h, % 60
Ratio of potential evapotranspiration to potential
evaporation (E
) 0.8
Where:
E
T
= Evaporation from open surface of water (or equivalent in heat energy)
∆ =-Slope of vapor pressure curve at t = tan α
H = Equivalent evaporation of the total radiation on the surface of plants (the final amount - net-
of energy finally remaining at a free water surface)
γ = psychrometer constant fixed device to measure humidity (= 0.66 if the temperature is
measured ic centigrade and e in mbar)
E = Actual vapor pressure of air at temp t
E
a
= Aerodynamic term (the term ventilation) depends on the air and low pressure steam
(evaporation for the hypothetical case of equal temperatures of air and water)
t = tan α = temperature
n/D = Cloudness ratio = (actual hours of sunshine) ÷ (possible hours of sunshine)
R
A
= Agot's value of solar radiation arriving at the atmosphere (assuming no clouds and a
perfectly transparent atmosphere).
h = Relative humidity
u
2
= Wind speed at height of 2 m (m/s)
hutE
D
n
htE
A
R
D
n
tE
D
n
tE
T
E,
2
,
4
,,
3
,,
2
,
1+
+
+
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
40
Solution
• Data: t = 22 °C, h = 0.6, sky Cover = 0.7, n ÷ D = 1 - sky cover Drag = 1 - 0.7 = 0.3,
evaporation = 2.8 m/day, evapotranspiration = 0.8 E
o
.
• From nomogram for temperature t = 22 ° C & the value of n/D = 0.3 then, E
1
= - 2.2
mm/day
• for Temperature t = 22 °C & tables for altitude 20 ° North In November, Agot's factor
RA = 666 gm cal/cm
2
/day & for values of n ÷ D = 0.3, E
2
= 2.6 mm/day
• For temperature t = 22 °C, n ÷ D = 0.3, h = 0.6 then: E
3
= 1.2 mm/day
• Potential evapotranspiration is given as = 2.8 mm/day
Since evapotranspiration = 0.8*evaporation, then evaporation predicted by Pennman is:
0.8xE
o
= 2.8 mm/day
Or E
o
= 3.5 mm/day
E
o
3.5 = E
1
+ E
2
+ E
3
+ E
4
= -2.2 + 2.6 + 1.2 + E
4
This gives E
4
= 1.9 mm/day
• For a temperature t = 22 °C, E
4
= 1.9 mm/day humidity & h = 0.6, determine the wind
speed at height of 2m, u
2
= 3.6 m/s then:
Example 3.2
Use the nomogram for the solution of Penman's equation to predict the daily potential
evapotranspiration from a field crop at latitude 20°N in December, under the following
conditions: (B.Sc., DU, 2012)
• Mean air temperature = 20
o
C
• Mean h = 70 %
• Sky cover = 60 %cloud
• Mean u
2
= 2.5 m/s
• Ratio of potential evapotranspiration to potential evaporation = 0.8
Solution
1. Data: latitude 20°N in December, t = 20° C, h = 0.7, sky Cover = 0.6, n ÷ D = 1 - sky cover
Drag = 1 - 0.6 = 0.4, u
2
= 2.5 m/s, ET = 0.8E
2. From nomogram for temperature t = 20° C and the value of n/D = 0.4 then, E
1
= - 2.4
mm/day
3. for Temperature t = 20° C and tables for altitude 20° North In December, Agot's factor R
A
=
599 gm cal/cm
2
/day and for values of n ÷ D = 0.4 , E
2
= 2.6 mm/day
4. For temperature t = 20° C, h = 0.7, n ÷ D = 0.4, then: E
3
= 1.4 mm/day
5. For Temperature t = 20° C, speed u
2
= 2.5 m/s and humidity h = 0.7, then: E
4
= 1.1 mm/day
6. Find E
o
from the equation:
E
o
= E
1
+ E
2
+ E
3
+ E
4
= - 2.4 + 2.6 + 1.4 + 1.1 = 2.7 mm/day
Evapotranspiration = 0.7 E
o
= 0.8 x 2.7 = 2.16 mm/day.
Example 3.3
Use the nomogram for the solution of Penman's equation to predict the daily potential
evapotranspiration from a field crop at latitude 20°N in November, under the following
conditions: (B.Sc., DU, 2012)
• Mean air temperature = 20
o
C
• Mean h = 80 %
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
41
• Sky cover = 70 %cloud
• Mean u
2
= 2.5 m/s
• Ratio of potential evapotranspiration to potential evaporation = 0.7
Solution
1) Data: sky Cover = 0.7, n ÷ D = 1 - sky cover Drag = 1 - 0.7 = 0.3, t = 20%C, h =
0.8, u
2
= 2.5 m/s, h = 0.8, Ratio of evapotranspiration to evaporation = 0.7.
2) From nomogram for temperature t = 20%C and the value of n/D = 0.3 then, E
1
= -
2.1 mm/day
3) For altitude 20% North In November, from table find Agot's factor R
A
= 666 gm
cal/ cm
2
/day. Then, for Temperature t = 20%C and for values of n ÷ D = 0.3 and for
R
A
= 666 the value of E
2
= 2.3 mm/day
4) For temperature t = 20%C, h = 0.8 and n ÷ D = 0.3, then: E
3
= 1.2 mm/day
5) For Temperature t = 20%C, and the speed u
2
= 2.5 m / s and humidity h = 0.8, then:
E
4
= 0.55 mm/day
6) Find E
o
from the equation:
E
o
= E
1
+ E
2
+ E
3
+ E
4
= - 2.1 + 2.3 + 1.2 + 0.55 = 1.85 mm/day
Evapotranspiration = 0.7 Eo = 0.7 x 1.85 = 1.295 mm/day.
3.2 Theoretical Exercises
10) Outline importance of the evaporation process (B.Sc., DU, 2012).
Solution
• Affects study of water resources due to its impact on yield of river basin&
• Affects the necessary capacity of reservoirs&
• Determines size of pumping plant&
• Affects water consumptive use of water by crops&
• Affects yield of groundwater supplies
11) Penman equation of evaporation from free water surfaces is based on two requirements that
need to be met for continuous evaporation to occur: there must be a supply of energy to
provide latent heat of vaporization, and there must be some mechanism for removing the
vapor, once produced. Define the parameters shown in the equation (B.Sc., DU, 2012)
12) What the main differences between evaporation and transpiration?
13) What are Penman's assumptions to estimate amount of evaporation?
14) Explain the different factors affecting the evaporation Process.
a- Sun Radiation
Process of conversion from water into vapor continue if the energy is there (sun light as
the heat energy), and cloud will be the barrier of evaporation.
b- Wind
If water evaporates to the atmosphere, the level between land surface and air will be
saturated by vapor, then the process of evaporation stops. In order for the process to
continue, the saturation level must be changed with dry air. This changing is possible if
there is wind which blows the component of vapor.
hutEh
D
n
tE
D
n
A
RtE
D
n
tE
T
E,
2
,
4
,,
3
,,
2
,
1+
+
+
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
42
c- Relative Humidity
If the relative humidity is up, so then the potential of air to absorb water is lower so
then the evaporation is lower as well.
d- Temperature
If the temperature of air and soil is high, therefore the process of evaporation will be
faster.
3.3 Problem solving in evaporation
1) Use Penman nomogram to solve its equation to predict daily evapotranspiration expected
from field plants at latitude 40 degrees north in April under the following conditions:
intermediate air temperature 20 degrees Celsius, average humidity 70%, coverage of sky
60% clouds, relative speed at a height of two meters 2.5 m/s, and expected rate of
transpiration to expected evaporation 0.7. What is the difference in the result at latitude of
40 degrees south? (Ans. 2.5 mm/day).
2) Use the Penman nomogram to solve its equation to predict expected daily evapotranspiration
from field plants at latitude 60 degrees north in March and June, for an expected rate of
evapotranspiration 65 percent of the expected evaporation, under the following conditions:
June
March
5.5
2.5
Intermediate air temperature (° C)
Medium relative humidity (%)
Sky coverage (% clouds)
Relative velocity of wind (m/s)
3) Use the nomogram for the solution of Penman's equation to predict the daily potential
evapotranspiration from a field crop at latitude 20°N in July, under the following conditions:
Mean air temperature = 20°C
Mean h = 80 %
Sky cover = 70 %cloud
Mean u2 = 2.5 m/s
Ratio of potential evapotranspiration to potential evaporation = 0.8
4) Find amount of water evaporating from the surface of a lake during the month of March if
the median value of the maximum temperature of the water surface of the lake during this
period is 18 degrees Celsius, the temperature of dry air at a height of 2 meters above the lake
surface is 20 degrees Celsius, humidity tw = 17 degree C, and the mean wind speed is 1.5
m/s at the same height noting that the lake did not freeze during the monitoring period.
5) Use the nomogram for the solution of Penman's equation to predict the daily potential
evapotranspiration from a field crop at latitude 20°N in November, under the following
conditions:
• Mean air temperature = 20
o
C
• Mean h = 80 %
• Sky cover = 70 %cloud
• Mean u
2
= 2.5 m/s
• Ratio of potential evapotranspiration to potential evaporation = 0.7
(Ans. 1.295 mm/day).
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
43
Chapter Four
Infiltration and percolation
4.1) Infiltration
Infiltration is the movement of water through the soil surface into the soil. Percolation is the
movement of water through the soil.
Factors affecting infiltration include: surface entry, transmission through soil, depletion of
available storage capacity in soil and ccharacteristics of permeable medium.
Important indicators used to estimate infiltration (as an average rate throughout the period of
storm and rain): Average infiltration method, φ-index, W-index and ratio to surface flow
method.
Φ - index is that rate of rainfall above which the rainfall volume equals the runoff volume. The
index reflects the average infiltration rate, and is found from a time-rainfall intensity curve.
Example 4.1
a) For a total rainfall distributed as shown in the table (1), the '-index of the catchment area
for a certain surface runoff is found to be 9 mm/hr (B.Sc., DU, 2012).
Table (1) Rainfall intensity versus time.
Determine:
i. The value of surface runoff producing a φ−index of 9 mm/hr.
ii. The total rainfall in the catchment area.
b) If the same rainfall had been distributed as shown in table (2), compute the value of
φ−index that would yield the same surface runoff. Comment on your answers.
Table (2) Rainfall intensity versus time.
Time Rainfall intensity (mm/hr)
0 0
1 5
2 15
3 20
4 20
5 14
6 1
Time Rainfall intensity (mm/hr)
0 0
1 9
2 16
3 27
4 10
5 8
6 5
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
44
Solution
a) φ-index is the average rainfall intensity above which the volume of rainfall equals the volume
of runoff.
b) Surface runoff = shaded area = (15-9)*1 + (20-9)*1 + (20-9)*1 + (14-9)*1 = 6*1 + 11* 1 +
11* 1 + 5*1 = 33 mm
Total rainfall = 5*1 + 15* 1 +20*1 + 20*1 + 14* 1 + 1* 1 =
75 mm
c)
Rate of rainfall above ' runoff volume
∴ (9 - ')×1 + (16 - ')×1 + (27 - ')×1 + (10 - ')×1 + (8 - ')×1 = 33 mm
∴ 70 - 5' = 33
∴ ' = 7.4 mm/hr
c) The results show that one determination of the f-index is of limited value and that many such
determinations should be made, and averaged, before the index is used.
40 m
20 m
5 m
Time
= 9
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
45
Example 4.2
For a total rainfall of 44 mm, distributed as shown in the table (3), establish the ' index of the
catchment area for a surface runoff of 19.5 mm. (B.Sc., DU, 2012)
Table (3) Rainfall intensity versus time.
Solution
Rate of rainfall above ' runoff volume
∴ (21 - ')×1 + (9 - ')×1 + (6 - ')×1 = 19.5 mm
∴ 36 - 3' = 19.5
∴ ' = 5.5 mm/hr
Example 4.3
A storm with 10.0 cm precipitation produced a direct runoff of 5.8 cm. Given the time distribution of the
storm as below, estimate the φ-index of the storm.
Incremental rainfall in each hour (cm)
Solution
Total Infiltration = 10.0 cm - 5.8 cm = 4.2 cm
Assume time of rainfall excess, te = 8 hr (for the first trial)
Then,
Time Rainfall intensity (mm/hr)
0 0
1 4
2 21
3 9
4 6
5 4
mm/hr
Intensity
mm/hr
9
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
46
hrcmIndex /525.0
8.50.10
,=
=
φ
But this value of φ makes the rainfalls of the first hour and eight hour ineffective as their magnitude is
less than 0.525 cm/hr.
The value of te is therefore modified.
Assume time of rainfall excess, te = 6 hr (for the second trial)
In this period, Infiltration = 10.0 - 0.4 - 0.5 - 5.8 = 3.3 cm
hrcmIndex /55.0
3.3
,==
φ
This value of φ is satisfactory as it gives,
te = 6 hr, and by calculating rainfall excesses.
Total Rainfall excess = 5.8 cm = total runoff
Example 4.4
Following data was obtained from a catchment of an area of 300 km2 .
each hour (cm)
Solution
Estimate the φ index of the storm if the volume of surface runoff was 16.5 x 106 m3.
Total Precipitation, P
= 0.5+1.1+1.3+2.1+2.0+1.8+1.0+0.4
= 10.2 cm
Runoff, R
= (16.5 x 106 m3) / (300 x 106 m2)
= 0.055 m = 5.5 cm
Total Infiltration, I = P - R
Assume time of rainfall excess, te = 8 hr
(for the first trial)
Then,
=
=
e
tRP
Index
φ
,
(10.2 – 5.5) / 8 = 0.5875 cm / hr
But this value of φ makes the rainfalls of the
first hour and eight hour ineffective as their
magnitude is less than 0.5875 cm/hr.
- The value of te is therefore modified.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
47
- Assume time of rainfall excess, te = 6 hr
(for the second trial)
In this period, Infiltration =
In this period, Infiltration =
=
=
e
tRP
Index
φ
,
This value of φ is satisfactory as it gives,
te = 6 hr, and by calculating rainfall excesses.
Time from
start (h) 1 2 3 4 5 6 7 8
Rainfall
Excess
(cm)
0 0.47 0.67 1.47 1.37 1.17 0.37 0
Total Rainfall excess = 5.5 cm = total runoff
4.2 Antecedent precipitation index, API
Example 4.5
The API for a station was 39 mm on first December, 61 mm rain fell on fifth December, 42 mm
on 7 December & 19 mm on 8 December. Compute the anticedent precipitation index for 12
December, if k = 0.89, and for same data assuming no rain fell (B.Sc., DU, 2013).
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
48
Where:
I
t
= Index value at t days later (mm)
I
o
= Initial value to the index (mm)
k = A recession constant, of magnitude between 0.85 to 0.98 (usually taken for a value of about
0.92).
t = Time (day).
Solution
• Data: amount of precipitation in different days of the month of December.
• Since I
t
= I
o
k
t
& assuming t = 1 produces: It = I
o
k. Which means that: the index for any
day is equal to that value in the previous day multiplied by the coefficient k. then if rain
fell in the day, the amount of rain is added to the index.
• For the given data index can be estimated as follows:
On the first day of December, I
1
= 39 mm
Just on the fifth day of December I
5
= 39*(0.89)
4
, I
5
= 24.5 mm
• On the fifth day of December rain should be added: I
5
= 24.5+ 61= 85.5mm
• Just on the seventh day of December I
7
= 85.5* (0.89)
2
, I
7
= 67.7 mm
• On the seventh day of December rain should be added : I
7
= 67.7 + 42 = 109.7mm
• Just on the eighth day of December I
8
= 109.7*(0.89)
1
, I
8
= 97.6 mm
• On the eighth day of December rain should be added: I
8
= 97.6 + 19= 116.6 mm
• In the twelfth day of December a I
12
= 116.6*(0.89)
4
, I
12
= 73.2 mm (with rain in
different days)
• In the case of lack of rain, the value is equal to the directory: I
12
= 39*(0.89)
11
, I
12
= 10.8
mm (with no rain)
4.3) Theoretical Exercises
1) The antecedent precipitation index (API) is given by: I
e
= I
o
k
t
, Define the parameters
shown in the equation (B.Sc., UAE, 1989, B.Sc., DU, 2012).
Solution
I
t
= index value at t days later
I
o
= initial value of index
k = a recession constant.
2) Define φ -index. (B.Sc., DU, 2012)
3) Illustrate major factors that may affect infiltration of water within a certain catchment
area. (B.Sc., UAE, 1989).
4) What are the factors affecting infiltration of water into the soil. How is it measured?
5) What is the difference between infiltration and percolation? Indicate how both are
measured.
6) Elaborate on indicators used to estimate infiltration. Which one is the best to be used?
Why?
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
49
4.4) Problem solving in infiltration
φ
φφ
φ-index
1) Find value of φ -index in a catchment area to make a surface runoff of 19.5 mm, knowing
that the total value of the rain for an area equal to 44 mm as indicated in the following
table: (Ans. 5.5 mm/hr). Rainfall intensity,
mm/hr
Time, hr
0
#
0
#
2) For a total rainfall of 150mm, distributed as shown in table (2), establish the ' index of
the catchment area for a surface runoff of 72mm. (B.Sc., UAE, 1989).
Time Rainfall intensity (mm/hr)
0 0
1 16
2 29
3 35
4 16
5 8
6 19
7 27
Compute the change in ' index for the same catchment in (b) above if a second storm
gave rise to a runoff equivalent to 68mm while the hourly rainfall of the storm indicated
the following distribution: (B.Sc., UAE, 1989). (ans. 11.7 mm/hr,).
Time Rainfall intensity (mm/hr)
0 0
1 15
2 43
3 61
4 24
5 12
3) Determine the ' index for the catchment area as affected by the two storms indicated in
(a) and (b) above. (B.Sc., UAE, 1989).
Solution
b) Rate of rainfall above ' runoff volume
∴ (16 - ')×1 + (29 - ')×1 + (35 - ')×1 + (16 - ')×1 + (19 - ')×1 + (27 - ')×1
∴ 142 - 6' = 72
∴ ' =
a) (43 – '
2
)×1 + (61 – '
2
)×1 + (24 – '
2
)×1 = 68
128 – 3'
2
= 68
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
50
60 =3'
2
'
2
= 20mm/hr
i.e. ' increased by 71%
API
4) The API for a station was 60 mm on the first of December, 64 mm rain fell on the
seventh of December, 40 mm rain fell on the ninth of December, and 30 mm rain fell on
the eleventh of December. Compute the API for 20
th
December if k = 0.85. Determine the
API if no rain fell (B.Sc., DU, 2012) .
Solution
Since I
t
= I
o
k
t
, letting t = 1, k
i
= kI
o
, I
o
= 60 mm
Just on 7
th
, index I
6
== 60x(0.85)
6
= 22.6 mm
After rainfall I
7
= 22.6 + 64 = 86.6 mm
Just on 9
th
, index = 86.6x(0.85)
2
= 62.6 mm
After rainfall I
9
= 62.6 + 40 = 102.6 mm
Just on 11
th
, index = 102.x(0.85)
2
= 74.1 mm
After rainfall I
11
= 74.1 + 30 = 104.1 mm
On 20
th
day, I
20
= 104.1x(0.85)
9
= 24.1 mm
If no rain fell
I
20
= 60x(0.85)
19
= 2.7 mm
5) The API for a station was 55mm on the first of August, 57mm rain fell 0n the sixth of
August, 31mm rain fell on the eighth of August, and 20mm rain fell on the 9
th
August.
Compute the API for 15
th
August if k = 0.92. Determine the APU if no rain fell. (B.Sc.,
UAE, 1989).
Solution
I
t
= index value at t days later
I
o
= initial value of index
k = a recession constant.
Since I
t
= I
o
k
t
, letting t = 1, k
i
= kI
o
, I
o
= 55mm
Just on 6
th
, index I
6
== 55x(0.92)
5
= 36.3 mm
After rainfall I
6
= 36.3 + 57 = 93.3 mm
Just on 8
th
, index = 93.3x(0.92)
2
= 79 mm
After rainfall I
8
= 79 + 31 = 110mm
Just on 9
th
, index = 110x(0.92)
1
= 101.2 mm
After rainfall I
9
= 101.2 + 20 = 121.2 mm
0
1
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
51
On 15
th
day, I
15
= 121.2x(0.92)
6
73.5 mm
If no rain fell
I
15
= 55x(0.92)
14
= 17.1 mm
6) The value of API for rainfall monitoring stations reached 42 mm in the first day of
October, and an amount of 46 mm of rain fell in the fifth day of October, and an another
amount of rain of 28 mm fell on the seventh day and 34 mm fell on the eighth day of
October. Find the value of the API on 12 October assuming that k is equal to 0.92.
Calculate the amount of index assuming no rain fell in the same period.
7) The API for a station was 55mm on the first of August, 57mm rain fell 0n the sixth of
August, 31mm rain fell on the eighth of August, and 20mm rain fell on the 9
th
August.
Compute the API for 15
th
August if k = 0.92. Determine the API if no rain fell. (ans. 17.1
mm).
8) The value of the API for a certain station reached 42 mm in first of August, and an
amount of 46 mm of rain fell on the fifth of August, also rain of 28 mm fell in the seventh
of august and 34 mm on the eighth day of August. Find value of the API guide on 12th
August noting that k is equal to 0.92. Calculate the amount of the index assuming no rain
in the same period. Draw curve of change of index with time. (Ans. 85.3, 1.7 mm).
9) The value of the API for a station reached 55 mm in the first of August, an amount of 57
mm of rain fell on the sixth of August, 31 mm of rain fell in the eighth of August and 20
mm fell in the ninth day of August. Find value of API index for the 15th day of August,
noting that k is equal to 0.92. Calculate the amount of index assuming no rain during the
same period. Draw curve of change of index with time.
10) The antecedent precipitation index (API) is given by: Ie = Iokt, Define the parameters
shown in the equation. The API for a station was 60 mm on the first of December, 64 mm
rain fell on the seventh of December, 40 mm rain fell on the ninth of December, and 30
mm rain fell on the eleventh of December. Compute the API for 20
th
December if k =
0.85. Determine the API if no rain fell. (ans. 24.1, 2.7 mm).
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
52
Chapter Five
Groundwater Flow
5.1 Groundwater
Groundwater represents that part of reserved water in a porous aquifer, resulting from infiltration
of rainfall through the soil and its penetration to the underlying strata.
Example 5.1
Outline mo st important factors affecting productivity of a well (B.Sc., DU, 2013).
Solution
Most important factors affecting productivity of a well:
• Lowering of groundwater within aquifer (drawdown aspects).
• Dimensions of aquifer & its lateral extent.
• Ground water storage.
• Transmissivity & specific yield or storage coefficient of aquifer.
• Conditions of flow (steady or unsteady).
• Depth of well.
• Establishment of well & methods of construction, properties & condition.
Table (5.1) Expected well yield.
Expected well yield, m
/d Well diameter, cm
< 500 15 400 to 1000 20 800 to 2000 25 2000 to 3500 30 3000 to 5000 35 4500 to 7000 40 6500 to 10000 50 8500 to 17000 60
Example 5.1
A well is 30 cm diameter and penetrates 50 m below the static water table. After 36 hr of
pumping at 4.0 m
3
/minute the water level on a test well 200 m distance is lowered by 1.2m and
in a well 40 m away the drawdown is 2.7 m (B.Sc., DU, 2013).
i) Determine radius of zero draw-down.
ii) Find coefficient of permeability.
iii) Compute draw-down in the pumped well.
iv) What is the transmissibility of the aquifer?
Using expected well yield estimates as presented in table (5.1), comment about yield of the well
as related to its diameter, & suggest a more suitable well diameter. Explain and validate your
answer (B.Sc., DU, 2013).
r
R
Ln
h
2
H (k
Q
2
−
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
53
Solution
• Given: D = 0.3 m, H = 50 m, r
1
= 40 m, x
1
= 2.7 m, r
2
= 200 m, x
2
= 1.2 m, Q
o
= 4000 L/min.
Find h
1
= h - x
1
= 50 – 1.2 = 47.3 m, and h
2
= h - x
2
= 50 – 2.7 = 48.8 m.
Use the following equation for both observation wells:
Where:
Q = Flow through well at a distance r.
h = Depth below original water level.
R = radius of zero draw down.
k = coefficient of permeability.
H = depth of aquifer.
By substituting given values into the previous equation, then:
• This yields R = 751.5 m.
• Find the permeability coefficient by using the data of one of the wells.
r
R
Ln
h
2
H (k
Q
2
−
=
well2nd
2
wellist
2
r
R
Ln
h
2
H (k
r
R
Ln
h
2
H (k
−
=
−
∴
well2nd
2
wellist
2
200
R
Ln
48.8
2
05
40
R
Ln
47.3
2
05
−
=
−
∴
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
54
Thus, for h = 50 m, ho = 47.3 m, r = 40 m, R = 751.5 m, Q = 4000*10
-3
*60*24 = 5760
m
3
/day,
• Depth of the water in the pumped well may be found as:
Or
This yields, h = 41.7 m
Determine the draw-down at the well as:
r
o
= h - h
1
= 50 – 41.7 = 8.3 m.
• From table for a well diameter of 30 cm, yield is 2000 – 3500 m
3
/d.
The amount of water abstracted is Q
o
= 4000 L/min = 4000*60*24/1000 = 5760 m
3
/d. A
better design would be selected a well of diameter of 40 cm (giving a yield of 5867 m
3
/d for
computed drawdown sat well (O.k. between 4500 - 7000).
Example 5.2
A fully-penetrating well, with an outside diameter of 50 cm, discharges a constant 3 m
3
/min
from an aquifer whose coefficient of transmissibility is 0.03 m
2
/s. The aquifer is in contact with
a lake 2 km away & has no other source of supply (B.Sc., DU, 2012).
i. Estimate the drawdown at the well surface. (Take R
0
as twice distance between
aquifer & lake).
ii. Using expected well yield estimates as presented in table (1), comment about yield of
the well as related to its diameter, & suggest a more suitable well diameter. Explain
your answer.
iii. Give a suggestion for type of aquifer soil with assuming an aquifer thickness of 20
m.
Solution a) Data: Q
0
= 3 m
3
/min = 3*60*24 = 5760 m
3
/day, R
o
= 2xL = 2x2000 = 4000m, D = 0.5
m, r
o
= 0.5/2 = 0.25 m, T = 0.03 m
3
/s = 1.8 m
2
/min
b) Using equation
c)
Hence, drawdown at well face is 2.57 m
• From table of well yield this diameter (= 0.5m) gives a yield in range of 6500 – 10000
m
3
/d. Nonetheless, needed constant flow ought to be 4320 m
3
/d. This suggests using a well
of diameter of 35 cm to have a yield between 3000 to 5000 m
3
/d, thus satisfying need.
r o
R
T
o
Q
r o
R
kH
o
Q
o
Sln
2
ln
2
ππ
==
25.0
4000
ln
8.12 3
×
=
π
o
S
m/d5.02
47.3
2
05 (
40
751.5
5760Ln
h
2
H (
1
r
R
QLn
k2
1
2
1
=
−
=
−
=
o
r
R
Ln
h
2
H (k
o
Q
2
o
−
=0.15
751.5
Ln
(x20.5
5760
2
05
o
r
R
Ln
(k
Q
2
H
2
0
h−=−=
daypermetercubic8675
)2/4.0( 5.751
2
7.41
2
505.20 =
−×
=
Ln
o
Q
π
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
55
• The actual well yield according to aquifer properties and computed drawdown for a well
diameter of 35 cm would then be:
• Determine permeability k = T/H = (1.8*60*24)/20 = 129.6 m/d
• From figure of range of permeability in natural soils, this value of k suggests fine
gravel soil description.
Example 5.3
A fully-penetrating well, with an outside diameter of 0.3 m, discharges a constant 4 m
3
/min
from an aquifer whose coefficient of transmissibility is 1.4 m
2
/min. The aquifer is in contact
with a lake 1.5 km away & has no other source of supply(B.Sc., DU, 2012).
i. Estimate the drawdown at the well surface. (Take R
0
as twice distance between aquifer &
lake).
ii. Using expected well yield estimates as presented in table (1), comment about yield of the
well as related to its diameter, & suggest a more suitable well diameter. Explain your
answer.
iii. Give a suggestion for type of aquifer soil with assuming an aquifer thickness of 25
m.
Solution
d) Theim's equation assumptions
• Aquifer homogeneous, isotropic & extended to infinity.
• Well penetrates thickness of aquifer carrying water & diverting water from it.
• Transmissivility (hydraulic conductivity) is constant in each place & does not depend on
time.
• Pumping is conducted at a steady rate for a period so that it may assume a stable
condition.
• Stream lines are radial (horizontal).
• Flow is laminar.
e) Data: Q
0
= 4 m
3
/min = 4*60*24 = 5760 m
3
/day, R
o
= 2xL = 2x1500 = 3000m, D = 0.3
m, r
o
= 0.3/2 = 0.15 m, T = 1.4 m
2
/min
f) Using equation
g)
Hence, drawdown at well face is 4.5 m
•
• From table of well yield this dia (= 0.3m) gives a yield in range of 2000 – 3500 m
3
/d.
Nonetheless, needed constant flow ought to be 5760 m
3
/d. This suggests using a well of
diameter of 40 cm to have a yield between 4500 to 7000 m
3
/d, thus satisfying need.
• The actual well yield according to aquifer properties and computed drawdown for a well
diameter of 40 cm would then be:
..50003000
3
4167
175.0
3000 57.2*24*60*8.122 KOandbetweendm
Ln
o
S
o
r
R
Ln T
o
Q=
×
===
ππ
r o
R
T
o
Q
r o
R
kH
o
Q
o
Sln
2
ln
2
ππ
== 15.0
3000
ln
4.12 4
×
=
π
o
S
..5760
3
5929
2.0
3000 5.4*24*60*4.122 KOdm
Ln
o
S
o
r
R
Ln T
o
Q>=
×
===
ππ
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
56
• Determine permeability k = T/H = (1.4*60*24)/25 = 80.6 m/d
• From figure of range of permeability in natural soils, this value of k suggests coarse
sand soil description.
Example 5.4
The static level of water table in an unconfined aquifer was 30 m above the underlying
impermeable stratum. A 150 mm diameter well, penetrating the aquifer to its full depth, was
pumped at the rate of 20 litres per second. After several weeks of pumping, the drawdown in
observation wells 20 m and 50 m from the well were 3.5 m and 2 m respectively, and the
observed drawdowns were increasing very slowly (B.Sc., DU, 2012).
a) Assuming equilibrium conditions, estimate the hydraulic conductivity and transmissivity
of the aquifer.
b) Estimate the drawdown just outside the pumped well.
c) What will be the yield of a 300 mm diameter well which will produce the same
drawdowns just outside the well and at the 50 m distance observation well in (2)? What
would be the drawdown at a nearer observation well?
Solution
Determine Q = 20x10
-3
x60x60x24 = 1728 m
3
/d
∴ R = 183.5 m
= 6.16 m/d = 7.13×10
-5
m/s
T = kH = 6.16×30 = 184.9 m
2
/d
b) Depth of water in well, = 203.8
H = 14.3 m
Drawdown just outside the = 30 – 14.3 = 15.7 m.
c) Yield
h = 28 m.
Drawdown at a nearer observation well = 30 – 28 = 2 m
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
57
Example 5.5
A well of diameter 30 cm penetrates an aquifer, water depth in it 15 meters before pumping.
When pumping is being done at a rate of 3000 liters per minute, the drawdowns in two
observation wells 40 and 90 meters away from the well are found to be 1.3 and 0.7 meters,
respectively (B.Sc., DU, 2011).
i. Find radius of zero draw down.
ii. Determine the coefficient of permeability, and
iii. Compute drawdown in the well.
iv. Give a suggestion for type of aquifer soil for the estimated permeability.
v. Comment on yield of well compared to its diameter. Give suggestions for
improvement?
Solution
1) Data: Data on wells.
2) Use following equation for each well
r = 40 m, s = 1.3 m, h = 15 – 1.3 = 13.7 m
r = 90 m, s= 0.7m, h = 15 – 0.7 = 14.3 m
i) Then: R can be found = 242.2 m
ii) Find coefficient permeability from the equation to the data of the well about 40 m
away: H = 15 m, s = 1.3 m, h = 15 - 1.3 = 13.7 m, R = 254.9 m, Q = 3000 L/min
= 3 m
3
/min = 3*60*24 = 4320 m
3
/day
iii) To find the depth of the water in the well, use equation with r
0
= 0.3/2 = 0.15 and k
= 66.36 m/d:
From which: h = 8.7 m
Then drop in level = drawdown at well face = s
o
= 15 – 8.7 = 6.3 m
iv) Suggest type of aquifer natural soil for a permeability of 66.36 m/day from given
diagram to be "coarse sand".
v) yield of well for a diameter of 30 cm ought to be between 2000 to 3500 m3/d.
Pumping done at a rate of 4320 m3/d which suggests a suitable diameter of 35
cm. Taking this diameter would yield a flow of:
r
R
Ln
hHk
o
Q
−π
=
2
0
2
wellndwellist
o
r
R
Ln
hHk
o
r
R
Ln
hHk
2
2
0
2
0
22
−π
=
−π
∴
2
2
1
2
90
3.14
2
15
40
7.13
2
15
−
=
−
R
Ln
R
Ln
dm
π
Ln
hH π
o
r
R
QLn
k/36.66=
7.13
2
15
40 29.242
4320
=
2
=
22
0
09.75
15.0 9.254
9.24 8.1576
2
15
22
0=
π
−=
π
−= Ln
xr
R
Ln
k
Q
Hh
dm
Ln
x π
o
r
R
Ln
hHk π
o
Q/34322=
15.0 2.242
.488
2
1536.66
=
2
=
22
0
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
58
Thus a well diameter of 35 cm is suggested for improvement to yield a flow rate of
4322 m3/d (within yield of 3000 to 5000, O.K.).
Example 5.6
A catchment area is undergoing a prolonged rainless period. The discharge of the stream
draining it is 140 m
3
/s after 15 days without rain, and 70 m
3
/s after 45 days without rain. Derive
the equation of the depletion curve in the form Q
t
= Q
o
*e
-at
and estimate the discharge 90 days
without rain (B.Sc., DU, 2012).
Solution 1) Data: values of flow rate Q
a
and Q
t
, after one month and eight days, Q
15
=
140, Q
45
= 70 m
3
/s. Q
90
?
2) Find base flow hydrograph equation using the equation: Q
t
= Q
o
* e
-at
Insert given values in the equation as shown below in equations 1 and 2.
140 = Q
o
× e
-15a
(1)
70 = Q
o
× e
-45a
(2)
Dividing equations 1 and 2 to find values of coefficient : a = 0.0231 /day
3) Substitute in one of the equations 1 or 2 to find the value of the primary
discharge Q
o
as follows:
Q
o
= 140 ÷ e-15x0.0231 = 198 m
3
/s
Basal flow hydrograph formula becomes: Q
t
= 198 * e
-0.0231t
4) Find rate of the stream after a period of 90 days in the watercourse
compensation in the basal flow hydrograph equation obtained in step 4
above: Q
t
= 198 * e
-0.0231 * t
, Q
90
= 198 * e
-90 × 0.0231
= 24.7 m
3
/s
Example 5.7
The static level of water table in an unconfined aquifer was 33.5m above the underlying
impermeable stratum. A 150mm diameter well, penetrating the aquifer to its full depth, was
pumped at the rate of 25 litres per second. After several weeks of pumping, the drawdown in
observation wells 20m and 50m from the well were 3.55m and 2.27m respectively, and the
observed drawdowns were increasing very slowly. (B.Sc., UAE, 1989).
a. Assuming equilibrium conditions, estimate the hydraulic conductivity and
transmissivity of the aquifer.
b. Estimate the drawdown just outside the pumped well.
c. What will be the yield of a 300mm diameter well which will produce the same
drawdowns just outside the well and at the 50m distance observation well in (ii)?
What would be the drawdown at a nearer observation well?
Solution
a)
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
59
225.2475 = 146.9371
= 0.6523
∴ R = 279 m
= 7.45×10
-5
m/s = 6.44m/d
T = kH = 6.44×33.5 = 215.7 m
2
/d
ii] Depth of water in well,
= 419.7
H = 20.49 m
Drawdown = 33.5 – 20.49 = 13.01 m.
iii] a)
b)
h = 31.2 m.
Drawdown = 33.5 – 31.2 =2.3m
Example 5.8
A well of diameter 0.3 m contains water to a depth of 50 m before pumping commences. After
completion of pumping the draw-down in a well 20 m away is found to be 5 m, while the draw-
down in another well 40 m further away reached 3 m. For a pumping rate of 2500 L/minute,
determine: • radius of zero draw-down.
• coefficient of permeability, and
• draw-down in the pumped well (B.Sc., DU, 2013).
Using expected well yield estimates as presented in table (5.1), comment about yield of the well
as related to its diameter, & suggest a more suitable well diameter. Explain and validate your
answer(B.Sc., DU, 2013) .
r
R
Ln
h
2
H (k
Q
2
−
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
60
Solution
• Given: D= 0.3 m, H = 50 m, r
1
= 20 m, x
1
= 5 m, r
2
= 40 m, x
2
= 3 m, Q
o
= 2500 L/min.
Find h
1
= h - x
1
= 50 - 5 = 45 m, and h
2
= h - x
2
= 50 - 3 = 47 m.
Use the following equation for both observation wells:
By substituting given values into the previous equation, then:
• This yields R = 119.7 m.
• Find the permeability coefficient by using the data of one of the wells.
Thus, for h = 50 m, ho = 45 m, r = 20 m, R = 119.7 m, Q = 2500*10
-3
*60*24 = 3600
m
3
/day,
• Depth of the water in the pumped well may be found as:
Or
This yields, h = 27 m
r
R
Ln
h
2
H (k
Q
2
−
=
well2nd
2
wellist
2
r
R
Ln
h
2
H (k
r
R
Ln
h
2
H (k
−
=
−
∴
well2nd
2
wellist
2
40
R
Ln
74
2
05
20
R
Ln
54
2
05
−
=
−
∴
m/d4.32
54
2
05 (
20
119.7
3600Ln
h
2
H (
1
r
R
QLn
k
2
1
2
1
=
−
=
−
=
o
r
R
Ln
h
2
H (k
o
Q
2
o
−
=0.15
119.7
Ln
(x4.32
3600
2
05
o
r
R
Ln
(k
Q
2
H
2
0
h−=−=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
61
Determine the draw-down at the well as:
r
o
= h - h
1
= 50 - 27 = 23 m.
• From table for a well diameter of 30 cm, yield is 2000 – 3500 m
3
/d.
The amount of water abstracted is Q
o
= 2500 L/min = 2500*60*24/1000 = 3600 m
3
/d. A
better design would be selected a well of diameter of 35 cm (giving a yield of 3683 m
3
/d for
computed drawdown sat well (O.k. between 3000 - 5000).
Example 5.9
a) A well penetrates into an unconfined aquifer having a saturated depth of 100 m. The
discharge is 250 litres per minute at 12 m drawdown. Assuming equilibrium flow conditions
and a homogeneous aquifer, estimate the discharge at 18 m drawdown. The distance from
the well where the drawdown influence are not appreciable may be taken to be equal for
both cases. (UAE, 1990).
Solution:
a) – Pollutants. - Difficulty.
b)
daypermetercubic3683
)2/35.0( 7.119
2
27
2
5032.4 =
−×
=
Ln
o
Q
π
100 m
Q = 250 L/min
H = 100 m
Q = ?
Q = 250 L/min
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
62
Example 5.10
A 0.3 m diameter well penetrates vertically through an aquifer to an impervious strata which is
located 18 m below the static water table. After a long period of pumping at a rate of 1 m
3
/min.,
the drawdown in test holes 14 and 40 m from the pumped well is found to be 2.62 and 1.5 m,
respectively.
i) Determine coefficient of permeability of the aquifer.
ii) What is the transmissibility of the aquifer?
iii) Compute the specific capacity of the pumped well. (UAE, 1990).
Solution:
a) r, Q, H, T, specific yield storage, lateral extent of aquifer, well
b)
i]
ii] T = kH = 13.5x18 = 243 m
3
/m.d
iii] depth in well h
2
= H
2
-
h = 9.1 m
drawdown (r) = H –h = 18-9.1 = 8.9m
iv] specific capacity = m
2
/s
5.2 Theoretical Exercises
1) Outline main differences between confined and unconfined aquifer (B.Sc., DU, 2012).
2) Define the term " Transmissibility of aquifer".
3) Define parameters that may influence yield of wells.
4) Comment on yield of well compared to its diameter. Give suggestions for improvement?
5) Define factor that may influence yield of wells. (B.Sc., UAE 1989,, UAE 1990, B.Sc., DU, 2012).
Solution
Drawdown, Flow conditions, Depth of penetration, Transmissivity, Specific yield, Storage
coefficient, Lateral extent of aquifer-, Construction and condition of well.
6) Outline most important factors affecting productivity of a well (B.Sc., DU, 2013).
Solution
Most important factors affecting productivity of a well:
• Lowering of groundwater within aquifer (drawdown aspects).
• Dimensions of aquifer & its lateral extent.
• Ground water storage.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
63
• Transmissivity & specific yield or storage coefficient of aquifer.
• Conditions of flow (steady or unsteady).
• Depth of well.
• Establishment of well & methods of construction, properties & condition.
7) In your opinion, what is the most significant source of groundwater pollution in this country?
Explain why is groundwater contamination so difficult to detect and clean up? (UAE, 1990).
8) Define parameters shown in Theims equation: (B.Sc., DU, 2012)
9) Define terms used in the equation used to estimate rate of constant pumping from a well penetrating
an unconfined aquifer (B.Sc., DU, 2011)
Solution
Q
o
= Steady state discharge from well, m
3
/s.
k = Coefficient of permeability, m/s.
H = Saturated thickness of aquifer, m.
h
o
= Depth below original water level at surface of well, m.
R = Radius of zero draw down (effective radius of drawdown), m.
r
o
= Radius of well, m.
10) Comment about Theim's equation (equilibrium equation) assumptions (B.Sc., DU, 2012).
11) Write briefly about groundwater in the kingdom of Saudi Arabia with reference to the following:
Occurrence, Distribution, Quantity & volumes of aquifers, Quality, Potential hazards & pollution
indicators.
12) What are the factors affecting groundwater flow within a basin?
13) What are Dubois's assumptions? How to use them?
14) What is the difference between stable and unstable flow?
15) What is the difference between a water table well and an artesian well? (Draw sketches).
5.3 Problem solving in groundwater
Well yield, drawdown
1) A 30 cm well is pumped at the rate of Q m
3
/minute. At observation wells 15 m and 30 m
away the drawdowns noted are 75 and 60 cm, respectively. The average thickness of the
aquifer at the observation wells is 60 m and the coefficient of its permeability amounts to 26
m/day.
i) Find the coefficient of transmissibility of the aquifer?
ii) Determine the rate of pumping Q.
iii) Compute specific capacity of the well.
iv) What is the drawdown in the pumped well? (UAE, 1989). (Ans. 65 m
2
/hr, 1.5
m
3
/min, 0.73 m
2
/min, 2 m).
2) A 30 cm well penetrates 45 m below the static water table. After a long period of pumping
at a rate of 1200 Lpm, the drawdown in the wells 20 and 45 m from the pumped well is
found to be 3.8 and 2.4 m respectively.
a) Determine the transmissibility of the aquifer.
b) What is the drawdown in the pumped well? (UAE, 1989). (Ans. 7.1 m
2
/hr, 13.5 m).
3) A 30 cm well serves a community of 1100 capita with a water consumption rate of 250
l/day. Under steady-state conditions the drawdown in the well was 1.83 m. The well
r o
R
T
o
Q
r o
R
kH
o
Q
o
Sln
2
ln
2π
=
π
=
r
R
Ln
o
hHk
o
Q
−
=
22
π
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
64
penetrates to an impermeable stratum 32 m below the water table of a homogeneous-
isotropic unconfined aquifer. Compute the change in the discharge in l/s for a well
drawdown of 1.83 m if the diameter of the well was: a) 20.0 cm; and b) 50 cm. Assume that
the radius of influence in all cases is 760 m. (SQU, 1991).
4) A community of 15,000 capita and fire demand of 35 L/s for six hours is to be served by a
50 cm diameter well. The well is constructed in a confined aquifer with a uniform thickness
of 15 m and hydraulic conductivity of 100 m/d. Two observation wells are installed at radial
distance of 50 m and 150 m. The drawdowns in the wells are 1.7 m and 1.3 m respectively.
Find:
a) the discharge of the well.
b) water consumption (l/c/d).
c) the power needed to lift the water to the ground surface if the original
piezometer level is 20 m below the ground surface. (SQU, 1991).
5) Gravity well is 50 cm in diameter. The depth of water in the well is 30 meters before the
start of pumping. Upon pumping at a rate of 2,100 liters per minute, Water level into a well
located about 10 meters away fell by 3 meters, and in another well 20 meters away it
reached 1.5 meters. Find: the zero draw down, permeability coefficient, and the drawdown
in the well. Draw cone of depression attributed to pumping of this well. (Ans. "m, #
m
3
/d, 5.2m).
6) The static level of the surface of the water in an unconfined aquifer is 33.5 meters above the
underneath impervious layer. Pumping is conducted at a rate of 20 liters per second in a 150
mm diameter well penetrating through all the depth of the aquifer. After a few weeks of
pumping in observation wells at 20 meters and 50 meters distance the water table decreased
by 3.55 meters and 2.27 meters, respectively and drawdowns increased gradually and
slowly. Assuming equilibrium conditions, find: hydraulic permeability coefficient,
drawdown directly outside the pumped well. Draw cone of depression attributed to pumping
of this well. How much is the yield from a 300 mm diameter well that can achieve the same
reduction in level outside it and at 50 meters from the observation well. Determine the
amount of drawdown in the neighboring observation wells?
7) A well of diameter of 0.3 meters was drilled to impervious base in the center of a circular
island of radius of one kilometer in a large lake. The well penetrates deeply in a sandstone
aquifer which have a thickness of 20 meters which lies below a layer of impermeable clay.
Permeability of sandstone equals 20 meters per day. Find steady flow on the assumption
that drawdown for the pezometic surface is equal to 2.5 meters in the well. (Assume R
0
=
twice distance between aquifer and lake). Comment about yield of well as related to its
diameter, and suggest a more suitable well diameter. Explain your answer.
8) A fully penetrating well, with an outside diameter of 0.2 m discharges water at a constant
rate of 6 cubic meters per minute from an aquifer whose coefficient of transmissivility is
20m
10m
3m
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
65
120 square meters per hour. The aquifer is in contact with a lake 1.6 kilometer away and has
no other source of supply. What is drawdown at the well surface? (R
0
can be estimated to
equal twice the distance between the aquifer and the lake).
9) Tube well 30 cm diameter penetrates an unconfined aquifer. Find the discharge from the
tube well on the assumption that the drawdown does not exceed 2 meters, the effective
height of the screen under conditions of the above mentioned drawdown is equal to 10
meters, the coefficient of permeability of the aquifer of 0.05 cm/s, and the radius of zero
drawdown Ro (effective radius of drawdown) is equivalent to 300 meters.
10) A gravity well of diameter 30 cm, water depth in it 20 meters before pumping is started.
When pumping is being done at a rate of 1800 liters per minute, the drawdown in a well 15
meters away is 2.5 m; and in another well 30 meters away is 1.5 m. Find:
• Radius of zero draw down.
• The coefficient of permeability, and
• Drawdown in the well.
• Give a suggestion for type of aquifer soil with compute permeability.
11) Find amount of water from a well of a diameter of 15 cm penetrating to the end of an
unconfined groundwater aquifer at a distance of 30 m from the ground surface, note that its
permeability is 3 meters/hour, and the level of drawdown in the well is 2.5 m, and radius at
which the level of groundwater vanishes is 400 meters.
12) A well of diameter 30 cm penetrates into groundwater aquifer. Water depth in it reaches 15
meters. The level of fall of water in two observation wells at distances 30 and 60 meters
from the well was found to be 1.2 and 0.5 meters, respectively, when the water pumping
rate is 4000 liters per minute. Find the permeability coefficient of the basin and the level of
drawdown in the well after pumping.
Aquifer recharge
13) A well was drilled an impervious base in center of a circular island of diameter 1.5 km in a
large lake. The well penetrates totally in the sandstone aquifer which is 18 meters thick
located under a layer of impervious clay. Permeability of sandstone equals 15 meters per
day. Find stable flow if the decline in pesimoetric surface should not exceed 2 meters in a
well of diameter 0.3 meters.
14) Suppose there are two canals, at different levels, separated by a strip of land 900 meter
wide, of permeability that reaches 0.3 m/hr. If one of the canals is 1.3 m higher than the
other and the depth of the aquifer is 15 m below the lower canal to an impermeable base,
find inflow into, or abstraction from, each canal per unit width. take annual rainfall as 2.1
m/annum and assume 85 percent infiltration.
15) A rain gauge in a particular area showed that the average annual rainfall was 610 mm. 75
percent of this rain seeps into the ground. There are two channels separated by a piece of
land of width 1 kilometer and permeability 0.35 m/h, one of the canals is higher than the
other by about 1.2 meters. Geological surveys in the region pointed out to the presence of a
groundwater basin under the piece of land with an average depth of 10 meters. Find the rate
of outflow from each channel to the aquifer.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
66
Chapter Six
Surface runoff
6.1 Runoff
Runoff is that part of the precipitation, as well as any other flow contributions, which appears in
surface stream of either perennial or intermittent form. It is the flow collected from a drainage
basin or water shed, and it appears at an outlet of the basin. Specifically, it is the flow, which is
the stream flow unaffected by artificial diversions, storage or any works of man in or on the
stream channel, or in the drainage basin or water shed.
Runoff divisions according to source of flow [1,2,15] include: surface flow (runoff), sub surface
flow (interflow, sub surface storm flow, storm seepage) and groundwater flow (groundwater
runoff).
Example 6.1
An unregulated river has monthly mean flows (in m
3
/s) as presented in the table (B.Sc., DU,
2013). Table: Monthly records of stream flow.
Month Monthly mean flow, m
/s
January 5.4
February 8.3
March 9.1
April 8.8
May 6.3
June 6.9
July 10.2
August 13.7
September 19.4
October 16.7
November 11.0
December 21.9
i) Allowing compensation water of 4.0 m
3
/s and reservoir losses of 0.5 m
3
/s, what storage
capacity of reservoir is required to ensure that, on average, no water is spilled? Assume
30-day months.
ii) What would the average net yield of the reservoir then be? (Hint: average net yield of
reservoir = demand – compensation – losses).
Solution
• Data: regular consumption of 230 m
3
/min, & data of monthly water flow.
• Find cumulative total flow as shown in the following table:
• Draw mass curve for data by plotting cumulative flow values as a variable with time.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
67
Month
(1) Flow, m
/s
(2) Flow, m
/day
(3)
= col (2)*60*60*24*30
Accumulated flow
in month, m
3
(4)
Jan 5.4 13996800 13996800
Feb 8.3 21513600 35510400
Mar 9.1 23587200 59097600
April 8.8 22809600 81907200
May 6.3 16329600 98236800
June 6.9 17884800 1.16E+08
July 10.2 26438400 1.43E+08
Aug 13.7 35510400 1.78E+08
Sept 19.4 50284800 2.28E+08
Oct 16.7 43286400 2.72E+08
Nov 11 28512000 3E+08
Dec 21.9 56764800 3.57E+08
• Determine monthly flow (column 3) by multiplying each value in column (2) by
60*60*24*30 to obtain flow in m
3
per month, taking a 30day month as assumed.
• Find accumulated flow as shown in column (4).
• Draw mass curve of reservoir.
• Draw demand line for no water spillage on curve.
• Find storage for no spillage as maximum ordinate between mass curve and demand =
66*10
6
m
3
= 66 Mm
3
.
• Determine demand = 356.9*10
6
/12 = 29.74*10
6
m
3
/month = 29.74*10
6
/(30*24*60*60) =
11.5 m
3
/s
• Compute average net yield of reservoir = demand – compensation – losses = 11.5 – 4 –
0.5 = 7 m
3
/s.
6.2 Flow mass curve (Ripple diagram, S-curve)
This is a graph of the cumulative values of hydrological quantity (such as: runoff) plotted against
time or data. Mass curve represents studying the effects of storage on the regime of a stream and
of determining the regulated flow by means of a flow summation curve. Mass curve is a curve on
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
68
which the ordinate of any point represents the total amount of water that has flowed past a given
station on a stream during the length of time represented by the magnitude of the abscissa to the
same point on the curve.
accumulative inflow
e
accumulation of storage draft line
d
Volume, V depletion of storage reservoir full
required storage, S
reservoir full c
Accumulated discharge b
point of critical period reservoir empty
slope representing constant outflow rate (demand)
أ t1 t2
time, t
Fig. 7.7 Mass curve or Rippl diagram
Example 6.2
1) Outline benefits of a mass curve (flow duration curve or Rippl diagram).
2) A water reservoir is designed to collect water from the adjacent catchment basin & to
regulate water use across an average regular flow of 230 cubic meters per minute. The table
below shows the monthly records of the stream flow. Find amount of storage needed to
keep up with regular consumption assuming no loss of water (B.Sc., DU, 2013).
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
69
Monthly records of stream flow.
Month Volume of water, million cubic meter
January 8
February 65
March 45
April 35
May 25
June 12
July 2
August 3
September 9
October 45
November 67
December 77
Solution
• Data: regular consumption of 230 m
3
/min, & data of monthly water flow.
• Find cumulative total flow as shown in the following table:
• Draw mass curve for data by plotting cumulative flow values as a variable with time.
Month Volume of
water (million
cubic meter)
Cumulative Volume
of water (million cubic
meter)
1 8 8
2 65 73
3 45 118
4 35 153
5 25 178
6 12 190
7 2 192
8 3 195
9 9 204
10 45 249
11 67 316
12 77 393
• Find the value of the annual use rate (for the month of December) =
230 (m
3
/min) × 60 (minutes/hour) × 24 (hours/day) × 365 (day/year) = 120.89 × 10
6
m
3
/year = 121 million m
3
/year.
• Draw draft line of uniform use from the point of origin to the point (a) on the mass curve.
• Draw a line parallel to the draft line from the point where the reservoir is full (b), & then
find the value of minimum required storage for the reservoir to keep pace with
consumption = 20 × 10
6
m
3
.
Example 6.3
The peak water consumption on the day of maximum water usage as follows: (B.Sc., UAE,
1989).
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
70
Time L/s Time L/s
Midnight 220 13 640
1 210 14 630
2 180 15 640
3 140 16 640
4 130 17 670
5 120 18 740
6 200 19 920
7 350 20 840
8 500 21 500
9 600 22 320
10 640 23 280
11 700 Midnight 220
Noon 660
i. Calculate hourly cumulative consumption values.
ii. Plot (i) to a mass diagram curve.
iii. What is the constant 24 hour pumping rate.
iv. Compute required storage capacity to equalize demand over the 24 hour period.
Solution
a) A graph of cumulative values of hydrologic quantities was plotted against time or
data: Reservoir condition, Demand, Spillage
b)
Time Hourly consumption Comulative consumption
L×106
L/s L×10
Midnight 220 792 0.792
1 210 756 1.548
2 180 648 2.196
3 140 504 2.7
4 130 468 3.168
5 120 432 3.6
6 200 729 4.32
7 350 1260 5.58
8 500 1800 7.38
9 600 2160 9.54
10 640 2304 11.844
11 700 2520 14.364
12 660 2376 16.74
13 640 2304 19.044
14 630 2268 21.312
15 640 2304 23.616
16 640 2304 25.92
17 670 2412 28.32
18 740 2664 30.996
19 920 3312 34.308
20 840 3024 37.332
21 500 1800 39.132
22 320 1152 40.284
23 280 1008 41.292
Midnight 220 792 42.084
11590/24= 487
Average =
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
71
6.3 Hydrograph
A hydrograph is a graph showing change of runoff volume (or stage, flow or discharge, velocity,
or any other properties of water flow) w.r.t. time.
Example 6.4
The ordinates of discharge hydrograph measured on a 1-hour interval at Lahat Bridge. Pari river, Ipoh.
(a) Concentration curve
b) Peak of hydrograph
c) Recession flood
d) Depletion curve
Groundwater recharge curve
Q
Fig. 6.2 Divisions of hydrograph
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
72
Constant Discharge Method
For t < ts Qb = Q
For ts ≤ t ≤ te Qb = qs
For t > te Qb = Q
Direct Runoff = Total Runoff - Base Flow
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
73
Equivalent Runoff Depth
Equivalent runoff depth = Total Direct Runoff X Time (s) + Area
Equivalent runoff depth
= (26.95 x 60 x 60) m
3
/ (271 x 10
6
) m
2
= 0.36 x 10
-3
m
= 0.36 x 10
-3
m x 100 cm/m
= 0.04 cm
Example 6.5
Solve the previous problem using the Concave Base Flow Method
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
74
Important data taken from the hydrograph
r
r
Notes: Qb at tp is taken as qm
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
75
Example of calculation at 11:00 pm
Qb = 13.22 + ( 11 – 10) [ (13.22 – 16.15) / (10 – 9) ] = 10.29 m3/s
Example of calculation at 1:00 am
Qb = = 7.36 + (1am – 12 am) [(16.55-7.36)/(1am – 12 am)] = 16.55 m3/s
Equivalent Runoff Depth
Equivalent runoff depth = Total Direct Runoff X Time (s) + Area
Equivalent runoff depth
= (26.95 x 60 x 60) m3 / (271 x 106) m2
= 0.36 x 10-3 m
= 0.36 x 10-3 m x 100 cm/m
= 0.04 cm
The general shape of the hydrograph of any ri
Example 6.6
A rain of 35mm fell on a drainage basin during a one-hour period beginning 01-00 hr. The
drainage area is 78.2 km
2
. The Following table lists the measured discharge from the area
with respect to time: (B.Sc., UAE, 1989).
Hour Discharge (m3/s) Hour Discharge (m3/s)
0 0 10 2.83
1 1 11 2.38
2 1.2 12 2.02
3 2.97 13 1.76
4 24.1 14 1.64
5 51.5 15 1.47
6 32.3 16 1.39
7 11.2 17 1.27
8 5.95 18 0
9 3.77
i. Develop a one hour unit hydrograph.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
76
ii. Compute and plot the hydrograph of surface runoff for two periods of heavy
rain occurring:
First storm: of 100mm rain between midnight and 02-00 hr; and 40mm rain
between 02-00 and 03-00 hr.
Second storm: of 89mm rain between 04-00 and 05-00 hr. Assume a constant loss
rate of 18mm and a constant base flow of 8m
3
/s
Comment about your graph.
Solution
a) A hydrograph of direct runoff resulting from 1mm of effective rainfall generated
evenly over the basin area at a uniform rate.
Uniform effective rainfall.
Uniform distribution over area.
Constant time duration
Ordinate ) Q
Reflects physical characteristics of basin.
b)
Total runoff = Q
a
.t = 8.264×18 (hr)×3600 = 535500 m
3
Measured discharge is to be divided by 6.85 to obtain UH
Time m3/s UH,
m3 /s/m Time Effective rainfall×UH Base
flow,
m3 /s
Total
rainfall,
m3 /s
82×UH 22×UH 71×UH
0 0 0 24 0 8 8
1 1 0.15 1 12.3 8 20.3
2 1.2 0.18 2 14.76 0 8 22.76
3 2.97 0.43 3 35.26 3.3 8 43.56
4 24.1 3.52 4 288.64 3.96 0 8 300.6
5 51.5 7.5 5 615 9.46 10.65 8 643.11
6 32.3 4.7 6 385.4 77.44 12.78 8 483.95
7 11.2 1.6 7 131.2 165 30.53 8 334.73
8 5.95 0.87 8 71.34 103.4 249.92 8 432.66
9 3.77 0.55 9 45.1 35.2 532.5 8 620.8
10 2.83 0.41 10 33.62 19.41 333.7 8 394.44
11 2.38 0.35 11 28.7 12.1 113.6 8 162.4
12 2.02 0.29 12 23.78 9.02 61.77 8 102.57
13 1.76 0.26 13 21.32 7.7 39.05 8 76.07
14 1.64 0.24 14 19.68 6.38 29.11 8 63.17
15 1.47 0.21 15 17.22 5.72 24.85 8 55.79
16 1.39 0.2 16 16.4 5.28 20.59 8 50.27
17 1.27 0.19 17 15.58 4.62 18.46 8 46.66
18 0 0 18 0 4.4 17.04 8 29.44
19 4.18 14.91 8 27.09
20 0 14.2 8 22.2
21 13.49 8 21.49
22 0 8 8
Example 6.7
Draw the Unit Hydrograph for a river which flows in a midsize catchment of 120 km2. After a 3-hour
storm, the following flow hydrograph was recorded at the catchment outlet flow station. Assume a
constant base flow of 10 m3 /s
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
77
Hydrograph
(Q) (m3/s)
of DRH
(m3 /s)
/ ER
of 3-hr UH
(m3 /s)
∑ DRH = 200 m3/s
Volume of Direct Runoff Hydrograph (DRH) = 60 X 60 X 3 X 200 = 2.16 Mm3
Drainage Area = 120 km2 = 120 X 106 m2
Runoff Depth = Effective Rainfall (ER) = rainfall excess =
= Volume of DRH / Area = ( 2.16 X106 ) / ( 120 X106 ) 0.018 m = 1.8 cm
Approximately 2cm
Check Area under UH = 1 cm over catchment area
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
78
of 3-hr UH
(m3 /s)
∑ DRH = 100 m3/s
Volume = 60 x 60 x 3 x 100 km2 = 1.08 x 106 m3
Drainage area = = 120 km2 = 120 x 106 m2
Runoff Depth = Effective Rainfall (ER) = rainfall excess =
= Volume of DRH / Area = (1.08 x 106) / ( 120 X106 ) = 0.01m = 1.0 cm (OK)
6.4 Theoretical Exercises
1) Write at length on each of the following: surface flow, subsurface and base flow,
catchment area and use of radioactive materials.
2) Indicate ways used to measure surface flow? Which one would be preferred?
3) What are the assumptions included in the rational formula to estimate surface runoff?
4) Define mass curve "Rippl diagram". What information does it provide? (B.Sc., UAE,
1989).
5) Outline benefits of a mass curve (flow duration curve or Rippl diagram) (B.Sc., DU,
2013).
6) What is meant by "Unit Hydrograph"? Indicate assumptions involved in the hydrograph
theory. (B.Sc., UAE, 1989).
7) Write briefly about each of: calibration curve, Rippl curve, and unit hydrograph.
8) How do you calculate peak flows upon availability and unavailability of measurements?
9) Write at length on each of the following: peak flows, and hydrograph.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
79
10) What are the sections of a hydrograph? Indicate methods of its analysis?
11) What are the ways used to separate base flow in a hydrograph?
12) Write a detailed report on each of: unit hydrograph, instantaneous unit hydrograph,
artificial water curve.
6.5 Problem solving in surface runoff
1) Suppose that the gauge shows a rise at the rate of 0.2 m/hr during a discharge
measurement of 100 m
3
/s and the channel is such that this rate of rise may be assumed to
apply to 1000 m reach of river between the measurement site and the reach control. let
the average width of the channel in the reach be 100 m, find the rate of change of
storage, and discharge measurement to be plotted on the rating curve.
2) A river discharge measurement made during a flood indicated Qa = 2500 cubic meters
per second. During the measurement, which took two hours, the gauge height increased
from 30.2 meters to 30.4 meters. Level readings taken at water surface 340 meters
upstream and 260 meters downstream of the observation site differed by 100 mm. the
river was 350 meters wide with an average depth of 3 meters at the time of
measurement. At what coordinate should the measurements be plotted on the rating
curve.
3) A water reservoir is designed to collect water from the adjacent catchment basin and to
regulate water use across an average regular flow of 8000 cubic meters per hour. The
following table shows the monthly records of the stream flow. Find amount of storage
needed to keep up with regular consumption assuming no loss of water
cubic meter)
(million cubic
meter)
4) A water reservoir is designed to collect water from the adjacent catchment basin and to
regulate water use across an average regular flow of 250 cubic meters per minute. The
following table shows the monthly records of the stream flow. Find amount of storage
needed to keep up with regular consumption assuming no loss of water.
Volume of water
(million cubic meter)
Month Volume of water
(million cubic meter)
Month
10
February
5
January
10
April
10
March
20
June
10
May
40
August
30
July
80
October
60
September
5
December
10
November
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
80
1) A catchment area is undergoing a prolonged rainless period. The discharge of the stream
draining it is 140 m
3
/s after 15 days without rain, and 70 m
3
/s after 45 days without rain.
Derive the equation of the depletion curve in the form Q
t
= Q
o
*e
-at
and estimate the
discharge 90 days without rain. (ans. Q
t
= 198 * e
-0.0231t
,24.7 m
3
/s).
2) A catchment area is undergoing a prolonged rainless period. The discharge of the stream
draining it is 100 m
3
/s after 10 days without rain, and 50 m
3
/s after 40 days without rain.
Derive the equation of the depletion curve and estimate the discharge 120 days without
rain. (ans. α = 0.0231 /day, Q
t
= 126.e
-0.0231t
, Q
t
= 7.9 m
3
/s)
3) A catchment area is undergoing a prolonged rainless period. The discharge of the stream
draining it is 4100 m
3
/min after ten days without rain, and 1200 m
3
/min after one month
without rain. Find:
• The equation of depletion curve.
• Estimate the discharge after a period of four months, and a period of six months
without rain.
• Draw depletion curve to scale.
4) A reservoir of water is built to collect the amount of rain water falling in the neighboring
basin and to organize a regular supply at a flow rate of 9000 cubic meters an hour.
Records of stream flow indicated the monthly data listed in the following table. Find
amount of storage required for organizing consumption (assuming no loss of water).
Volume of water, million m
Month January # February March April May June July August September October # November December
5) Define mass curve "Rippl diagram". What information does it provide? The peak water
consumption on the day of maximum water usage as follows:
Time L/s Time L/s
Midnight 220 13 640
1 210 14 630
2 180 15 640
3 140 16 640
4 130 17 670
5 120 18 740
6 200 19 920
7 350 20 840
8 500 21 500
9 600 22 320
10 640 23 280
11 700 Midnight 220
Noon 660
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
81
v. Calculate hourly cumulative consumption values.
vi. Plot (i) to a mass diagram curve.
vii. What is the constant 24 hour pumping rate.
viii. Compute required storage capacity to equalize demand over the 24 hour
period. (ans. 487 L/s).
6) A reservoir is designed to collect water precipitated in the adjacent watershed basin adjacent
to regulate supply with an average regular flow equal to 284 cubic meters per minute. The
following table shows monthly capacity records of the stream in millions cubic meters. Find
amount of storage needed to cope with regular consumption assuming no loss of water.
volume of water
(million cubic
meters*
month volume of water
(million cubic
meters)
Month
49.4 February 7.4 January 34.6 April 31.1 March 8.6 June 29.6 May 1.2 August 2.5 July 58 October 18.5 September 82.7 December 76.5 November
13) Measuring the flow of a river during the flood has shown that Qa = 2700 cubic meters per
second. During the measurement, which lasted for two hours, the level increased from
40.48 meters to 40.36 meters. Readings on the surface of the water level at distances of
390 meters and 310 meters upstream and downstream from the observation site differed
by about 100 mm. the river width is 400 meters, and its average depth is 3.5 meters
during the instant of measurement. Find coordinates of the point that should be taken for
measurements in the calibration curve. (Ans. 40.42, 2639).
14) Flow rate in a stream discharging water shed runoff is 50 cubic meters per second after
five days without rain. The flow rate in the stream reaches half that amount after 20 days
without rain. Determine the equation of basal flow hydrograph. Find the amount of flow
rate after a period of six months in the watercourse.
15) The rate of flow in a water course discharging an water shed is 3950 m
3
/min after ten
days of no rain, and the flow rate is 1190 m
3
/min after one month with no rain. Draw the
time hydrograph for base flow baseband; and find the amount of flow rate after three
months and a five-month period in the watercourse. (Ans. 0.54, 0.015 m
3
/s)
16) Calculate and draw one-hour unit hydrograph for a certain area and noting that the
discharge area is 70 square kilometers and the runoff from one peak for a precipitation of
the 20 mm as represented in the following table:
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
82
Measured runoff, m
/s Time, hr 0 3.8 11.2 # 12.3 11.1 8.2 5.2 3.6 3.2 2.8 2.6 2.3 2.2 # 2.1 2.0 1.9 1.8 1.7 1.6 1.5 1.4 #
#
17) Calculate and draw one-hour unit hydrograph for to a certain area noting that 35 mm of
rain fell in a discharge area of 78.2 square kilometers over a period of one hour started at
01:00. The following table shows the measured runoff through the area in the unit of
time: Measured flow, m
/s Time, hr Measured flow, m
/s Time, hr #" #" #"# # "# # " #" " #" " " " #" "# "# " "
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
83
Chapter Seven
General exercises
7.1 Complete missing titles
a) Complete missing titles by using the following words and phrases ( Cyclonic
precipitation. Relative humidity. Residence time. Unconfined aquifer. Depression
storage. Speed. Detention storage) (B.Sc., DU, 2012)
1) When rainfall is more than infiltration, evaporation and evapotranspiration, ponds gather
in low-lying areas representing depression storage.
2) Filling of ponds and depressions Leads water to flood on the surface of the earth
representing detention storage.
3) Residence time is the average duration for a water molecule to pass through a subsystem
of the hydrologic cycle.
4) Relative humidity describes the ability of air to absorb additional moisture at a given
temperature.
5) In absence of other factors tending to influence wind, it should be expected that its
direction would be from areas of high pressure towards areas of low pressure and that its
speed would vary with the pressure gradient.
6) Cyclonic precipitation is linked to passage over low temperature areas or altitude,
resulting in lifting of hot air masses over cold masses.
7) Dupuit assumptions are used to get an approximate solution for uni-dimensional flow in
an unconfined aquifer.
7.2 True/false sentences
b) Indicate whether the following sentences are true (T) or false (F ): (B.Sc., DU,
2012)
1) The intensity and frequency of the hydrological cycle is independent of geography and
climate. (F)
2) The role of applied hydrology is to provide guidance for planning and management of
water resource. (T)
3) Earth's rotation affects movement of wind and rotation of seasons. (T)
4) Trees reduce wind speeds and cool the air as they lose moisture and reflect heat upwards
from their leaves. (T)
5) Frost or ice is dew falling from the sky and freezes on the ground. (T)
6) Rate of evaporation decreases as the specific gravity increases. (T)
7) Infiltration devices estimate quality of infiltration rather than the quantity. (T)
8) A confined or artesian aquifer is one that is separated from the surface by an aquiclude or
aquitard. (T)
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
84
c) Indicate whether the following sentences are true (T) or false ( F): (B.Sc., DU,
2011)
1) Hydrology concerns quantity and quality of moving water accumulating on land, soil &
rock adjacent to the surface. (T)
2) At night, treetops act as radiating surfaces, and the soil beneath is protected from
excessive heat losses. (T)
3) In a large city, the amount of heat which is produced annually roughly equals the solar
radiation reaching an equivalent area. (F)
4) Winds are mainly the result of horizontal differences in pressure. (T)
5) Thiessen method is suitable to find number of medium height to a particular area served
by a network of fixed stations in number and positions. (T)
6) Method of isohyets maps give results close to those obtained by the Thiessen network
method. (T)
7) Intensity of rainfall indicates quantity of rain falling in a given time. (T)
8) Clouds slow up the process of evaporation. (T)
9) Rate of evaporation is more for salt water than for fresh water. (F)
10) The role of applied hydrology is to provide guidance for planning and management of
water resource. (T)
7.3 Underline the best answer
d) Underline the best word or phrase (between brackets) to offer a useful sentence.
(B.Sc., DU, 2011)
1) (Residence time /Time of concentration/Rainfall duration) Average travel time for water
to pass through a subsystem of the hydrologic cycle.
2) (Relative humidity /Dew point/Humidity) is the ratio of amount of moisture in a given
space to the amount a space could contain if saturated.
3) A (psychrometer/barometer/desiccant) is an instrument used to determine the value of
humidity..
4) Actinometer and Radiometer are general names for instruments used to measure
(intensity /duration/frequency) of radiant energy.
5) Winds are mainly the result of (vertical/horizontal /direct) differences in pressure..
6) An instrument for measuring the speed or force of the wind is called (wind
vane/anemometer/pyranometer).
7) One of the most important reasons leading to condensation of vapor is (heating by air
masses mix/ Dynamic or adiabatic heating/Contact and radiational cooling ).
8) Moist air moving up the side of a mountain facing the prevailing wind causes
precipitation to fall in a process known as(conventional/orographic/Cyclonic) lifting.
9) One of the factors that influence the measurement of rainfall, especially solid part, is
(wind /pressure/latitude).
10) Of appropriate recommendations to determine the number of monthly rain gauge stations
those proposed by (Bleasdale/Buys Ballot/Thiessen).
e) Underline the best word or phrase (between brackets) to offer a useful sentence.
(B.Sc., DU, 2012)
11) (Orographic precipitation / Conventional precipitation / Cyclonic precipitation)
occurs due to obstruction of topographical barriers (mountains, natural hills) for
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
85
winds laden with moisture and its mechanical lifting to top layers and then its
expansion & cooling resulting in rainfall.
12) (Horton / Penman / Agot) equation shows a mathematical method to estimate the
infiltration capacity curve.
13) Φ
ΦΦ
Φ - index / W-index / Antecedent Precipitation Index) is that rate of rainfall above
which the rainfall volume equals the runoff volume.
14) (Water-bearing layers / Man-made reservoirs / Cones of depression) are
geological features with permeability and they have components that allow
significant movement of water through them.
15) (Meteoric water / Juvenile water /Magmatic water) denotes water in or recently
from the atmosphere.
7.4 Rearrange groups
f) Rearrange group (I) with the corresponding relative ones of group (II) in the area
allocated for the answer in table (1). (B.Sc., DU, 2012)
Table (1) Matching of relevant words or phrases..
Group (I) Rearranged group (II) Group (II)
Hydrology Water science Radiometer
Holton Infiltration Baric
Psychrometer Humidity Penman
Connate water Fossil water Unconfined
Condensation Contact cooling Water science
Perched aquifer Unconfined Limestone
Actinometer Radiometer Fossil water
Evapotraspiration Penman Infiltration
Karst Limestone Contact cooling
Buys-Ballot Baric Humidity
g) Rearrange group (I) with the corresponding relative ones of group (II) in the area
allocated for the answer. (B.Sc., DU, 2011)
Group (I) Rearranged group (II) Group (II)
Humidity Psychrometer Rainfall depth
Contact cooling Dew Convective precipitation
Acinometer Radiant energy Anemometer
Convectional
precipitation Convective precipitation Hygrometer
Wind speed Anemometer Frontal precipitation
Wind direction Wind vane Dew
Orographic precipitation Windward slopes Radiant energy
Relative humidity Hygrometer Wind vane
Cyclonic precipitation Frontal precipitation Psychrometer
Rain gauge Rainfall depth Windward slopes
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
86
7.5 Complete missing titles
h) Complete missing titles by using the following words and phrases (Dew point, Storage,
Speed, Humidity, Temperature, Runoff, Actinometer, Water, Science, Wind) (B.Sc.,
DU, 2011)
1) Hydrology is derived from the Greek words: hydro (meaning water) and logos (meaning
dealing with science).
2) When rainfall is more than infiltration, evaporation and evapotranspiration, ponds gather
in low-lying areas representing Depression storage.
3) The water cycle consists of precipitation, evaporation, evapotranspiration and runoff.
4) Earth's rotation affects movement of wind and rotation of seasons.
5) Relative humidity is the percentage of actual vapor pressure to saturation vapor pressure.
6) Dew point denotes the temperature at which space becomes saturated when air is cooled
under constant pressure & with constant water vapor pressure.
7) Shading effect of trees tends to keep daily maximum temperature somewhat lower.
8) Actinometer and Radiometer are general names for instruments used to measure
intensity of radiant energy..
9) Speed would vary with the pressure gradient.
7.6 Match the words or phrases
Adiabatic Snow Orographic Precipitation Convective
Rain gauge Cyclone Dew Winds Baric
i) Match the words or phrases above with the definitions below. (B.Sc., DU, 2011)
1. Are mainly the result of horizontal differences in pressure.. ( Winds )
2. Baric wind law is also known as Buys-Ballot's law. ( Baric )
3. In this type of cooling heat is not added from outside sources. ( Adiabatic )
4. It forms when cloud droplets (or ice particles) in clouds grow and combine to
become so large that the updrafts in the clouds can no longer support them, and
they fall to the ground. (Precipitation)
5. Represents frozen water pouring from clouds in small pieces. ( Snow )
6. Denotes water vapor that condenses in the cold layers of the atmosphere during
the night and falls to the ground in small droplets. (Dew)
7. In this type of precipitation most rain is deposited on the windward slopes.
Other parts are located under the rain shadow and remain dry. (Orographic)
8. This type of precipitation results from lifting of air converging into a low
pressure area. (Cyclone)
9. This type of precipitation is local, and its intensity varies from light rain
showers to dangerous and destructive thunderstorms. (Convective)
10. It is a hydrological instrument used in measurements to calculate rainfall for the
region. (Rain gauge )
* +,-../012.,34.,,-+5.,06"7/89:;<8.=>?@ABCDE
#"*
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
87
#* FGH-4I.F@E## ° ,/F.J12K,L1+M,N8,OK,LPQ,-"R
STUB1+SVVWO+,-../012.,34.,,-+5.,06"XQ,/,YVZ1+MN
#"* W
*[\P]Y8 T ^## ° +2 h ^R #\:H-QMN8OKL8_`1QB@Q-## ° 2 MN8OKL8_ es ^"#aA=
\ST]Y8b. h = 100*e/es PSVVWOKLQ- 81 = 100×e ÷ 19.82 SVVWOKL e ^ 16.05 "aA= \PS8EMNSTUBQ- PMNSTUB e - es ^"# – "^" "aA=
\_>Q/X1S:W-4FZ?9XQ/YVZQ-@E= es 1 e. P3?81 es ^"
XQ,/,YVZ4,Bc@,08=,;@1+:W-44Bc@08dTaA=
8V1`1QB@ es ^"+
^XQ/YVZ" ° 2
*
2QOe
38/4SZ nomogram ,VH],GZ@,M,_.8S.f/@F0d4W
gK9Y,e.I.F:H-4Ph1ijWG?FkST`8k-4
#Ye.8.J1+.l-4 I8,YG1+R9,,/,1+DW,eR
@mnG#"M,_.8,OM,_.8f,/,Z1+o" K,9,,B/S,Tp,n,"
g$[./--4
+-../012.34.-+5.06#"*
W
^I8,YGP]Y8*q"+DW,e t ^# ° +2 h ^" + u2
^#"+o2 n r D ^ – "
^"+
h ^" • :H-Qh-.8/@ t ^# ° 8_12 n r D ^"Q-1? E1
^\#"2.= • Q1:,,H,,- t ^# ° g,,K,,O`1Q,,B@,,12 ° ,,Q,,-1?+,,?S,,T1`8,,k
P]./_ RA
^ 2.=es.E=-
=V1 n r D ^"4Bc@08 E2
^" 2.= • :H-Q t ^# ° +2 n r D ^"+ h ^" PQ-1? E3
^"2.= • :H-Q t ^# ° .J1+2 h ^" e1+ u2
^#"PQ-o2 E4
^"2.= • Q- Eo
P48@
Eo = E1 + E2 + E3 + E4 = -2.3 + 3.7 + 1.3 + 1.1 = 3.6 mm/day E trans = 0.7 E o = 0.7 x 3.6 =2.5 mm/day. • gKOE[./--4
* tK_,,,1+].,,,3.Z,,,_1+XQ,,,/,,,YVZPf,,,L11h,,,",,,Y8s1,,,G,,,Y1+s,,,u
+-../012.34.-+5.06
2002 "*
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
88
*
FGH-4I.F@E#" ° /F.J12KL1+MN8OKLPQ-"R
,-../012.,34.,,-+5.06"XQ/YVZ1+MNSTUB1+SVVWO+ 2002 "*
* 38/4SZ2QOe nomogram ,VH],GZ@M_.8S.f/@F0d4W
gK9M,_.8,OM,_.8f,/,/+.Z.15sFkST`8k-4
P,,h1,,ij,,WG+M,,_.8,,O,,8
34.,,,,-+5.,,06+,,-../012.,, 2002 "*
5 .Z. ,,,,-46Y,,,,e.I.,,,,F:,,,,H,,,,-4
*.l
*R6Ye.8/.J
*DWeR6I8YG
*o6v/
5.5
2.5
* 34.,-+5.06"VY/STK_5_pJ@ABCDE+,-../012.,
#"*
* 2QOeS`1QBSTWL.8]Z1].8wVJDEeg.HSTKe.8K_Q-
+-../012.34.-+5.06"@;VJ1SWKe.8SVJ#"*
YW8*=e6Y85V*R6H8 1 2 3 4 5 6
62 74 87 112 118 107
8 11 20 18 33 10
!
\ ]YW8STYx]QVGP]Y8 . #\ P48@SWKe.8VJ2QOeVY/8JFYxKe.Q- Pav = % (Pi /n) Pav ^
6 62 " 74 87 12 18 107 r*^""=e
\YxKe.Q-P@GVJ2QOeVY/8JF
Pmean=(A1/A)P1+(A2/A)P2 +..+ (An /A)Pn
6^K,,,,,e.8K_,,,,,#y y y##yy yr*^
""=e
* HKe.8Y8:QkQ-,N,E:,44,,eFGQ1A`zM.E
]./,,e,,YV/,,YgT,,6,,0,,-1+=,,"2.,,34.,,,,-+5.,,06"*
+-../01#"*
W
P]Y8 A ^+M.E t ^+V_4 t* ^"+ P ^"=
P48@Ke.8Y8{8eZQ-1?
^#:QH1eST= "
mm26
6.5 63.0
130
*
t
A3.0
1PP =
−=
−=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
89
* K_48_ API |YW8#QVKVe1+}Y~?@`1xST= 46 Y8@=
jY>8E+}Y~?@}OST 28 1d/MSTYx@= 34 @,@,•2.,ST=
K_48_Q-"}Y~? API 2.#3€w8}Y~? k s1G"#QVDH"Q
,-+5.,06"@U,M,QG9/W/=e":n}nZSTY?`.Y>2QgT
+-../012.34.#"* W
" 3?8 It =I0 k t
gn1 t ^P•/ It =I0 k {,GQ,,8_s1,G2.,sxd,Z?S,/s‚"
ST1ƒa2.ST8V,
k „‚,>,8_h,ƒG2.,s?S,T,Y?j,Y>…c=,;@,1"
"QYx
#" P:Y8]Z@
}Y~?@`1x2.ST I1
^#= • }Y~?@}O2.ST:k I5 =42(0.92)4= 30.1 mm • Y>s‚Y8TLcS/}Y~?@}O2.ST I5 ^ 30.1 46 ^ 76.1 = • }Y~?@M2.ST:k I7 =76.1(0.92)2= 64.4 mm • Y>s‚Y8TLcS/}Y~?@M2.ST I7 ^ 64.4 28 ^#"= • }Y~?@@•2.ST:k I8 =92.4(0.92)1 = 85 mm • Y>s‚Y8TLcS/}Y~?@@•2.ST I8 ^ 34 ^= • }Y~?@NSZ•2.ST I12 =119(0.92)4 = 85.3 mm
" PY?`.Y>2QHST I12 =42(0.92)11 = 1.7 mm
" 9/W/8=e
* „QVSWYe[ZoQHtBVY/8XT4QVQ-",Y8,8_3?=,…c+=,
s1GM_.8S034.,,-+5.,06PS,`1Q,BS,T@,8=,VDH=
+-../012.#"*
*e6@U*e=6Yx:Qk n†
#
n†
#
W
\ ^SWYe[Z‡`Q1+@UMYx:QkGP]Y8"= . #\ x:Qk=e?"ˆ`•8Hˆ0k9@8E@UMY \ P482QOeXT4Q-1?
SWY[Z‡`Q^XT@9Z4xYx`Q^XT@9xYx`Q "^6#\ φ *6\ φ *6\ φ * Ps1XT44Bc@08=;@1 φ ^"e=
#* DE$…,81$,ƒT?,VJs?",,VY/S,TK_,5,_p,J@ABC ,-+5.,06
+-../012.34. 2002 "*
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
90
* 2QOeS`1QBSTWL.8]Z1].8wVJDEeg.HSTKe.8K_Q-
.-+5.06"s18Y8‰.YVJ+-../012.34 2002 "*
Y8‰.Y
*†.6s18
@:.uW80H8
*0>6g.W41QH
@QV8Ke.8Y8
*†.6‰.YO 6 <
5 <
4 <
3 <
2 <
1 <
1 >
#
#
6.5
5.5
4.5
3.5
2.5
1.5
0.9 #$% K_48_ API W8|YQVKVe1+}Y~?@`1xST= Y8@=
jY>8E+}Y~?@54ST1d/@•STYx@= 20 @,M,e2.,S,T=
K_48_Q-"}Y~? API 2.3€w8}Y~? k s1G"#QQVDH"
,-+5.,06"@U,M,QG9/W/=e":n}nZSTY?`.Y>2QgT
+-../012.34. 2002 "* &'() * Q3€w81xŠ48@*sY_aT46A.BWST.-3USTV8aTQ4Q-aT Q [Zte1+%g.W…nZ1+ k sYV„BG‡1+ r [./8gnOZ1+ S +
3UOa81 H .V/81+ T +-../012.34.-+5.06"#"*
#* >Y_ 4 :U- UE ST 08 :QV 9c l ]nH"":E :W ST .E
Gd08e|s‚SBW3USTEl"0@JVJjWGs‚1w
s1GSBW…nZfY[./8gnOZ3E…cV8aTQQ-"2.ST
A1B3?DB‡s1U#>Y_SlST"-+5.06"34.
+-../012.#"*
*
-+5.06$A.BW~‹1A.WB ST.-3U[Z‡@]_1n=>?
+-../012.34.#"*
* +-../012.34.-+5.06"Œm.Z?@ABCDE#"*
* >Y_…-l"=elSTI8a8`Q8•ƒI-cQ/1"•ƒQ_
#F/ QGlST [./8bnOZ +V_QST QV8l ST|1 Ž
QGX?#S.HF/bnOZw"snu[./8gnOZY_•uZPQ-" zero draw down / 1 +gnOZ ‰1O =e "l ST [./8 gnOZ1 +…n
+-../012.34.-+5.06"l„‚F•ƒI-@[./8#"*
!
482QOe*q
φ
5
10
20
15
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
91
o
r
R
Ln
h
2
Hk
o
Q
2
0
−π
=
r = 10 m, h = 0 – 3 = 27 m r = 20 m, h = 0 - 1.5 = 28.5 m
wellnd2
2
0
wellist
2
0
o
r
R
Ln
h
2
Hk
o
r
R
Ln
h
2
Hk
−π
=
−π
∴
l
22
20
R
Ln
5.28
2
30
10
R
Ln
27
2
30
−
=
−
PF/1 R ^"
Q,,GS,,,,l],,Z,,48@,,,,…n/,,Q,,-,, H ^+2 h ^# +2 r ^+2 R ^ 41.5 +2 Q ^#^V_4#"2
^V_4#2
2.
d/m8
27
2
30
105.41
Ln3024
h
2
H
o
r
R
QLn
k
22
0
=
−π
=
−π
=
8a84BtP482QOelSTI
615
25.0 5.41
Ln
8x
30242
40
r
R
Ln
k
Q
2
H
2
0
h=
π
−=
π
−=
PF/1 h ^#"
^[./8STgnOZ‡
– #"^"#
* $ST.Bl-Zc9:;<8.=>? n +-../012.34.-+5.06
## "*
* ST.-3UST I8fY@EX.8A.BW~""dWG08VYp.T
`Q8•ƒX-?#>Y_lSTZ•STQ1"3UOa8EG=
Q9_‹ST •ƒ@Me?Mƒ#1QV8[./8bnOZF/
"1#"# ‡QUG1DG9w3A.Gh1•gT"IK1wBQG]LnOZ
"I8 F/ •ƒ S l q :k [./8 gnOZ1 +01QF …n/ PQ-
>Y_l-ZQV"l„‚F•ƒI-@[./8gnOZ‰1O=e=
[./8STgnOZ‡}nZaVWGF/08Q91F-:k1"_8l@
+-../012.34.-+5.06$:1B8_8lST[./8gnOZQV
2002 "*
*!+,-.
* +,,-../012.,,34.,,,,-+5.06$dGQO,,e=,,>?,,1$[,,Z‡9,,/W/Q,,.T,,
#"*
#* 4Q,,t=,,i/d,,1,,B8DE,,g.,,WS,,T,,JF„,,8],,8EM,,8BS,,m4.,,=8,,†
`4,=i,/aT4Ke.8#X,B8a,T4]zB,eS,`1Q,B@,1",V_QS,Tw,0w,
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
92
,,e‡,,E.82Az,,3U,,OQ,,VQ,,-"D,,08,,3.8,,:Q,,VF,,NS,,8=i,,/8'zF
+-../012.34.-+5.06"I83QVT2QgT#"* FN3.6I8=BH
*D0
FN3.6I8=BH
*D0 / 7.4 T 49.4
5 31.1 c 34.6
. 29.6 .Z. 8.6
.. 2.5 }Y~? 1.2
8e 18.5 .E? 58
8T.Z 76.5 84 82.7
!
\ =i/8'zFe‡P]Y8#2
"FNI8aT4]Z1+V_4 #\ PS`1QBSTfL.8ES0S8EaTQQ- FN3.6I8=BH
*D0
3.6S8EaTQ
*D0
/ 7.4 7.4
T 49.4 56.8
5 31.1 87.9
c 34.6 122.5
. 29.6 152.1
.Z. 8.6 160.7
.. 2.5 163.2
}Y~? 1.2 175.2
8e 18.5 193.7
.E? 58 251.7
8T.Z 76.5 328.2
84 82.7 410.9
\ aTQ=_=e]Z*KYO60aT49/W/=e@,8E@U/8ES8E
"`•8H=e9 \ ^*84FN6s./=i/8'zFe‡`Q8_Q- #6 2
*V_4y*eV_46y#*2.e6yIIII^*/e2.6y
2
"/e \ 9*?6YV/9c†xYVZ@=i/8DWK=e0aT49/W/ . \ @U,OG,_?,8_Q,-=,;+*[6w,l83U,OFT3.0SYV/@DWKOwA.wY=e
III^'zFe‡E.8m4.8[.Yy
2
"
* +-../012.34.-+5.06$SXB8SWY[Z‡'QV"•E#"*
* _3?3?3ƒn`zFZ3-5 Qa
^# 8es‚5VI/;?1"Z•STD0
@5V8mnG4A@"9c"9I8fYe9X.8]I_"
T1.,Wj,n,_8M,_.@„Z4?1„BG9?S",8
F/gdV8Ke.1+"S/SYV/];QHcQ-"]eViWI/;?
+-../012.34.-+5.06":89/W/ST]eVFL1#"*
! ^gP]Y8^a8+2"+2 dt ^#e
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
93
dh = 0. – 0. = 0.12 m
a8^H8y^g"y^2
#
I8Ye.8 v¯ ^ Qa r A ^# 6ry"^*"o2
P48@bn-.eQ- s/m51.293.1x3.1
A
a
Q
3.1U ===
^8"r ^"y
–
22
48@`Q83BQ-P
s/
3
m2639
4
10x43.1x53.2
5
10x67.1
1
2700
dt
dh
US
A
1
a
Q
orrected
Q=
−
−
+
=
+
=
6^5V8Ke.8mnG‡Q-""r*#^"#
6:89/W/9];QHt=;@1"#+#*
*
DEeg.H•u2.VSXBST[Z‡`Q|2
:,N9,ƒQ,,V_4
‡`Q|18E+:Y~2?dT[Z2
"Y?3zY>~@Fk9ƒQV_4
@F,k:,T9,ƒQ,[,Z‡`Q,Q,VQ-1ŽsQV[ZzS/US89/W/8=e
+-../012.34.-+5.06"S8XB8Fk?8:T1#"* ! \ `Q=_P]Y8[Z‡ Qt 1 Qa
+F,k12,?,Z8;9,ƒQ Q10 ^+
Q30 ^
2
V_4
#\ P482QOesQV[Z‡h-1Q>4Q- Qt = Qo *e-% t
@48ST„Z4?fL.8E48ST:Y8=Vg.1#"
^ Qo
y e
\y%
6*
^
Q o
y e
\y%
6#*
@488V1#g.W8_4Bc@08 % X1 :
a ^"
\ @48XQHcSTg.1?#SQ‡hu8_4Bt Qo PSGŒE Qo 6^rr*
e
\y"
^#
2
o
\ PsQV[Z‡h-1Q>4fuG Qt = 120*e-0.06*t
2
3
4
40
م 1.5 م
م
s/m
5
10x67.1
s3600x2 m12.0
dt
dh −
==
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
94
\ h-1Q,>,4S,Tb.,S,8XB8,F,k?,;z;:,T9,ƒQaTQ`QQVQ-
:.YOSTF`.uW=GSsQV[Z‡P„z? Qt = 120*e -0.06*t
Q
^
#
r
e
\y"
^"
2
o
\ :.YO09ƒQaTQ`Q4Bc@08„z?PS8XB8Fk?8:T
d1 Q
^
#
r
e
\y"
^"
2
"o
* FnuGH3€w8d8e1/VY/8:QH.h-1Q>DH "M.E
K_,:Q,H1,8_@3B1
20 34.,,-+5.,06PS,`1Q,BS,T,•8,8E=,
+-../012.#"*
@U
*e6
5V83B
26
*o
0
3.8
# 11.2
12.3
11.1
8.2
5.2
3.6
3.2
2.8
2.6
2.3
# 2.2
2.1
2.0
1.9
1.8
1.7
1.6
1.5
# 1.4
#
!
_:kh-F/*4"unSee?3-Q-.‡=;@1+Y81UQ1
5V83BKe. Qav
P:Y8]Z@
s/
3
m93.3
21 5.82
aV
Q==
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
95
5V83BKe.^S03By@U
3
m297000
h
s
3600
h
21
s
3
m
93.3t
av
Q
T
Q=××=⋅=
S03B^3B{8erH8
mm24.4m
3
1024.4
6
1070
297000
depthrunoff =
−
×=
×
=
QV8`1QBST5V83B=V=;@1"#":QH.h-1Q>9`.uW
PKe.8aTQ
Qav
^#"r#^"2
o
@U
*e6
5V8aTQ
26
*o
:QH.h-1Q>
26
*o
#
#
#
#
0 3.8 11.2 12.3 11.1 8.2 5.2 3.6 3.2 2.8 2.6 2.3 2.2 2.1 2.0 1.9 1.8 1.7 1.6 1.5 1.4 0
0 0.9 2.64 2.9 2.62 1.93 1.23 0.85 0.75 0.66 0.61 0.54 0.52 0.5 0.47 0.45 0.42 0.4 0.4 0.35 0.33 0
5.82Q =
∑
* iZS,,,T,,,-Q8],,,LT‡f,,,L11":Q,,,H.h-1Q,,,F4.,,,uV8,,,"h-1Q,,,F,,,
+-../012.34.-+5.06 2002 "*
* K_,,Y83€w8d8e1/VY/8:QH.h-1Q>DHH,S,T=,
FH•uG "#,]?Q,:Q,H1e:TM.EPS,`1Q,B1"
3B@+,-../012.34.-+5.06P@U:QH1STH85V8 2002 "*
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
96
@U
*e6
5V83B
26
*o
*e6@U5V83B
26
*o
#"
#"
#"# # #"#
#" "
#" "
" "
#" "
"# "#
"
"
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
97
References and bibliography for further
reading
1) Abdel-Magid, I. M., & Ibrahim A. A., Hydrology, Sudan University for Science and
Technology Publishing House, SUSTPH, Khartoum, Sudan, 2002 (In Arabic).
2) AlHassan, E. E., Adam, O. M. A., Ahmed, M. G. and Abdel-Magid, I.M., Hydrology and
hydraulics, submitted for publishing Sudan University for Science and Technology
Publishing House, SUSTPH (In Arabic).
3) Bedient, P. B., Huber, W. C. and Vieux, B. E., Hydrology and flood plain analysis,
Prentice Hall Inc., Upper Saddle River, NJ, 2008.
4) Bruce, J. P. and Clark, R. H., Introduction to hydrometeorology, Pergamon, Oxford,
1966.
5) Brutsaert, W., Hydrology: An Introduction, Cambridge University Press, NewYork,
2005.
6) Gupta, R. S., Hydrology and hydraulic systems, Waveland Press, Inc., 2001.
7) Hammer, M. J., and MacKichan, K. A., Hydrology and quality of water resources, John
Wiley and Sons, NewYork, 1981.
8) Hiscock, K., Hydrogeology: Principles and Practice, Wiley-Blackwell.
9) LaMoreaux P. E., Soliman, M. M, Memon, B. A., LaMoreaux, J. W. and Assaad, F. A.,
Environmental hydrology, CRC Press Publishing, Taylor & Francis Group, LLC, Boca
Raton, FL, 2009.
10) Linsley, R. K. Kohler, M. A. and Paulhus, J. L. H., Hydrology for Engineers, McGraw-
Hill Series in Water Resources and Environmental Engineering, 1982.
11) Linsely, R. K.; Kohler, M. A. and Paulhus, J. L. H., Applied Hydrology, Tata McGraw-
Hill Pub. Co., New Delhi, 1983
12) Maidment, D., Handbook of hydrology, McGraw-Hill Professional Publisher; 1993.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
98
13) Raudkivi, A. K., Hydrology – An advanced introduction to hydrological processes and
modeling, Pergamon Press, Oxford, 1979.
14) Shaw, E.M., Hydrology in practice, Van Nostrand Reinhold Co., Berkshire, 1983.
15) Subramanya, K. Engineering Hydrology, 3rd Edition, International Edition, McGraw-
Hill, London, 2009.
16) Ven Te Chow, Ed., Handbook of Applied Hydrology: A Compendium of Water
Resources Technology, McGraw Hill Book Co., New York, 1964.
17) Viessman, W. and Lewis, G. L., Introduction to hydrology, Fifth Edition, Prentice Hall
Publisher; 2012.
18) Ward, A. D. and Trimble, S. W., Environmental Hydrology, CRC Press Lewis
Publishers, Boca Raton, FL, 2004.
19) Wilson, E. M., Engineering Hydrology, Macmillan Education, London, 1990.
20) Wisler, C. O. and Brater, E. F., Hydrology, John Wiley and Sons, New York, 1959.
21) WMO, Guide to Hydrological Practices, 5th ed., WMO-No. 168. Chapters 20, 21, 1994.
22) WMO, Terakawa, A., Hydrological Data Management: Present State and Trends, WMO-
No. 964. 2003, Available through http://www.wmo.int/.
23) WMO, Guide to Meteorological Instruments and Methods of Observation, 2008,
Available through http://www.wmo.int/.
24) Faris, G.F. and Abdel-Magid, I.M., Fundamentals of engineering hydrology, under
preparation.
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
99
Useful Formulae
1) e = P - P'
2) ew - e = γ (t - t w)
3)
x
x a
a
x b
b
x c
c
P
N P
NN P
NN P
N
=
+ +
3
4) av i
i
n
PP
n
=∑
=1
5)
mean
P
i
Ai
P
i
n
i
A
i
n
==
∑
=
∑
1
1
6) i
t b
=+
7) i
tn
=
8) iP
t
=60
9) P Nt =× ×
−
0 282
5
1214 10
600 254
.
..
10)
11)
( )
it Nt =
× −
60 202 3 2 54
0 282.
. .
12) EV = b(es - e)
13) EV = f(u). (es – e)
14) T a
EH E
=
+
∆
γ
15)
16) f = fc + (fo – fc )e-kt
17)
( )
F fdt f t o
fc
f
f
ke
c
tk t
f
= = + −
∫−
01
18)
c
Fo c
f
f f
k
=
−
19) It = Io *kt
20) v = k*i
21) n
ve
=
22)
d
k
v
−=
23) q kH d
dx
= − .
24) k
−=φ
∑
=
÷
∑
=
=+++= n
1i i
A
n
1i i
P
i
A
n
P
A
n
A
......
2
P
A
2
A
1
P
A
1
A
mean
P
*
t
A3.0
1
P
P−=
h,
2
u,t
4
Eh,
D
n
,t
3
E
D
n
,
A
R,t
2
E
D
n
,t
1
E
T
E+
+
+
=
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
100
25)
q kh dh
dx
= −
26)
qk h h
x
o
=−
2 2
2
27) 22
hNx
kax b = − + +
28)
π
=− 1
r 2
r
Ln
kH2 o
Q
2
S
1
S
29)
−π
=
1
r 2
r
Ln
2
S
1
SH2
Q
k
30)
1
r 2
r
Ln
2
S
1
S2
Q
kHT
−π
==
31) r o
R
ln
T2 o
Q
r o
R
ln
kH2 o
Q
Sπ
=
π
=
32)
o
r
R
Ln
2
o
h
2
Hk
o
Q
−π
=
33)
−−
π
=− 2
r
2
R
k2
N
r
R
Ln
k
Q
2
h
2
Ho
34) Q = π R2 N
35) Q = 27.78 C I A
36) Q = 2.78 C I A
37)
( )
ic T
t d
m
n
=+
38)
+= dt
dh
US
1
1
Q
a
Q
39) A
a
Q
3.1U =
40)
+= dt
dh
S
a
Q3.1 A
1
Q
a
Q
41) Q = k ( h – a ) x
42) S
H
rACQ ⋅=
43)
44) N = b.A0.2
45) 0.5 (I1 + I2 )∆ t – 0.5 (D1 + D2 )∆ t = S2 – S1
46) S = K[XI + (1 - X)D] 47)
t
S
DI ∆
∆
=−
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
101
Appendices
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
102
Problem Solving in Engineering Hydrology Faris Gorashi Faris & Isam Mohammed Abdel-Magid
103
Appendix (B) Agot's values of short-wave radiation flux R
A
at the outer limit of the
atmosphere in g cal/cm
2
/day as a function of the month of the year & latitude.
Year Dec. Nov. Oct. Sept. august July June May April Mar Feb. Jan Latitude
(Degrees)
N 90
# # # # # # # # # # # # # Equator # # # # # # # #