CE3141 Module 1 - UNIT 3 and 4
CE3141 Module 1 - UNIT 3 and 4
Infiltration can be schematically modelled in two situations. One is under low intensity
rainfall, and the other under high intensity rainfall.
Given these,
Try to understand the infiltration that occurs in soil due to light rainfall and compare this
to when rain is continuous like during the monsoon season. Is there a difference? In the
models given above, can you differentiate what it means by actual and potential
infiltration?
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 34
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
Infiltration is the flow of water into the ground through the soil surface. Infiltration volume
is subtracted from a precipitation event in order to determine the net volume of rainfall,
or rainfall excess, which is equivalent to the direct runoff from a watershed area.
The distribution of soil moisture within the soil profile during the downward movement of
water is illustrated by the figure below.
𝒇 = 𝒇𝒄 + (𝒇𝟎 − 𝒇𝒄 )𝒆 𝒌𝒕
Horton’s Equation
Where:
𝑓= infiltration capacity 𝑓 = initial infiltration capacity
𝑓 = final/ equilibrium infiltration capacity 𝑘=empirical constant
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 35
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
Illustrative example 1. The initial infiltration capacity of a watershed is estimated as 1.5
in./hr., and the time constant is taken to be 0.35/ hr. The equilibrium capacity is 0.2
in./hr. Use Horton’s equation to find (a) the values of 𝑓 at 𝑡 = 10 min, 30 min, 1 hr., 2 hrs.,
and 6 hrs., and (b) the total volume of infiltration over the 6-hr period.
(1.30)
𝐹 = 0.20 𝑖𝑛 ℎ𝑟 (6ℎ𝑟) + 1 − 𝑒−(0.35/ℎ𝑟)(6ℎ𝑟) 𝑭 = 𝟒. 𝟒𝟓𝟗𝟒 𝒊𝒏
0.35
To check,
𝐹 = 0.20 + (1.30)𝑒 .
𝑑𝑡 = 0.2𝑡 + (1.30 − 0.35)𝑒 . ]6
0
𝑭 = 𝟒. 𝟒𝟓𝟗𝟒 𝒊𝒏
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 36
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
Philip (1957), using analytical solutions for unsaturated flow, developed two
equations of the form,
Considering the vertical column of soil of unit horizontal cross-sectional area between
the surface and depth L, the cumulative infiltration, the Green and Ampt equation for
cumulative infiltration is given by:
𝐹(𝑡)
𝐹(𝑡) = 𝐾𝑡 + 𝜓∆𝜃 ln 1 +
𝜓∆𝜃
where ∆𝜃 = (1 − 𝑠 )𝜃 ,
𝑠 = effective saturation and, 𝜃 = 𝜂 − 𝜃 is effective porosity
In 1983, Rawls, Brakensiek, and Miller used the Brooks-Corey equation to analyze
approximately 5000 soil horizons across the United States and determined average
values of the Green-Ampt parameters for different soil classes.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 37
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
It can be seen in that as the soil particles get finer, from sands to clays, the saturated
hydraulic conductivity Ks decreases, the average wetting front suction c increases
(negatively), and porosity us is variable.
Illustrative example 2. Compute the infiltration rate 𝑓(𝑡) and cumulative infiltration 𝐹(𝑡)
after one hour of infiltration into a silt loam soil that initially had an effective saturation of
30 percent. Assume water is ponded to a small but negligible depth on the surface.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 38
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
The cumulative infiltration at 𝑡 = 1 ℎ𝑟 is determined from the equation,
𝐹(𝑡)
𝐹(𝑡) = 𝐾𝑡 + 𝜓∆𝜃 ln 1 +
𝜓∆𝜃
𝑐𝑚 𝐹(𝑡)
𝐹(𝑡) = 0.65 𝑐𝑚 ℎ𝑟 (1ℎ𝑟) + 5.6745 ln 1 + 𝑐𝑚
ℎ𝑟 5.6745
ℎ𝑟
𝑭(𝒕) = 𝟑. 𝟏𝟔𝟓𝟔 𝒄𝒎
𝑐𝑚
𝜓∆𝜃 5.6745
𝑓(𝑡) = 𝐾 + 1 = 0.65 𝑐𝑚 ℎ𝑟 1 + ℎ𝑟 = 1.8152 𝑐𝑚/ℎ𝑟
𝐹(𝑡) 3.1656 𝑐𝑚
The infiltration methods discussed above uses the assumption that water is ponded to a
small depth on the soil surface so all the water that the soil can infiltrate is available at
the surface. However, this ponding only results when the rainfall intensity is greater than
the infiltration capacity of the soil.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 39
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
Infiltration rate 𝑓(𝑡) and cumulative infiltration
𝐹(𝑡) are related by the Green-Ampt equation by,
𝜓∆𝜃
𝑓(𝑡) = 𝐾 +1
𝐹(𝑡)
The figure shows that cumulative infiltration at
ponding time 𝐹 is equal to
𝐹 = 𝑖𝑡 and 𝑓 = 𝑖,
𝑲𝝍∆𝜽
𝒕𝒑 =
𝒊(𝒊 − 𝑲)
Illustrative example 3. Compute the ponding time and the depth of water infiltrated at
ponding for a silt loam soil of 30 percent initial effective saturation, subject to rainfall
intensity of 5 cm/h.
Since the parameters given are the same as in illustrative example 3, we may use the
values for
𝐾 = 0.65 𝑐𝑚 ℎ𝑟 and 𝜓∆𝜃 = 5.6745
In solving for the ponding time,
𝑐𝑚
𝐾𝜓∆𝜃 0.65 𝑐𝑚 ℎ𝑟 5.6745
𝑡 = = ℎ𝑟 = 0.1696 ℎ𝑟 (≈ 10.1760 𝑚𝑖𝑛𝑠)
𝑖(𝑖 − 𝐾) 𝑐𝑚 𝑐𝑚 𝑐𝑚
5 5 − 0.65
ℎ𝑟 ℎ𝑟 ℎ𝑟
This means that it takes 10.1760 minutes to saturate the soil.
In this case, the infiltration rate 𝑓(𝑡) is equal to rainfall intensity 𝑖 at ponding.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 40
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
However, after ponding 𝑡 > 𝑡 , the infiltration depth increases exponentially
according to this equation,
𝜓∆𝜃 + 𝐹
𝐹 − 𝐹 − 𝜓∆𝜃 ln =𝐾 𝑡−𝑡
𝜓∆𝜃 + 𝐹
where infiltration depth varies with time. Only the cumulative infiltration is unknown at at
any time after ponding. The infiltration rate is later obtained using the equation,
𝜓∆𝜃
𝑓(𝑡) = 𝐾 +1
𝐹(𝑡)
Illustrative example 4. Calculate the cumulative infiltration 𝐹 and the infiltration rate 𝑓
after one hour (𝑡 = 1 ℎ𝑟) of rainfall of intensity 5 cm/hr on a silt loam soil with an initial
effective saturation of 30 percent.
Since the parameters given are the same as in illustrative examples 3 and 4, we may use
the values for
Substituting,
𝜓∆𝜃 + 𝐹
𝐹 − 𝐹 − 𝜓∆𝜃 ln =𝐾 𝑡−𝑡
𝜓∆𝜃 + 𝐹
𝑐𝑚
𝑐𝑚 5.6745 +𝐹
𝐹 − 0.8479 𝑐𝑚 − 5.6745 ln ℎ𝑟 = 0.65 𝑐𝑚 ℎ𝑟 (1 − 0.1696 ℎ𝑟)
ℎ𝑟 𝑐𝑚
5.6745 + 0.8479 𝑐𝑚
ℎ𝑟
𝑭 = 𝟑. 𝟎𝟏𝟔𝟓 𝒄𝒎
This is the cumulative infiltration (depth) after an hour of rainfall of intensity of 5 cm/hrl.
Comparing the results in example number 2 and 4, the infiltration depth obtained under
continuous ponding is more than the infiltrated depth after an hour of 5cm/hr of rainfall. This is
because of the time it took for ponding to occur in the first 0.1696 hr and during this period, the
infiltration rate is also less than its potential value.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 41
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
Another infiltration method is the use of
the ∅-index. It assumes that infiltration loss
is uniformly distributed across the rainfall
pattern where the area above the ∅-
index sums to the volume of direct runoff
and the area below sums to the volume
of infiltration loss.
Illustrative example 5. Use the rainfall data below to determine the ∅-index for a
watershed that is 0.875 square miles, where the runoff volume is 228.7 ac-ft.
Infiltration when taken away from precipitation is equal to the direct runoff from a watershed
area, 𝑃 − 𝐼 = 𝑅. This may be written as,
Given the runoff volume of 228.7 ac-ft, its equivalent depth can be computed by
dividing it the runoff volume with the watershed area,
228.7 𝑎𝑐𝑟𝑒 𝑓𝑡
𝑟𝑢𝑛𝑜𝑓𝑓 𝑑𝑒𝑝𝑡ℎ = = 0.4084 𝑓𝑡 ≈ 4.9007 𝑖𝑛
640 𝑎𝑐𝑟𝑒𝑠
0.875𝑚𝑖 ×
𝑚𝑖
Substituting this in our first equation,
(𝑖 − ∅)∆𝑡 = 4.9007 𝑖𝑛
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 42
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
We can now determine the ∅-index by trial and error. In choosing a starting point, it is
logical to choose any of the rainfall intensities.
The first trial yields a ∅-index of 0.9599 in/hr. But since it is not equal to the original
assumption of 1.1 in/hr., another iteration is needed.
We can use the outcome of the first iteration as the starting point of the next iteration
such that we remove the 4th and 5th terms (since they correspond to intensities less
than 0.9599 in/hr.)
(1.4 − ∅)2ℎ𝑟 + (2.3 − ∅)3ℎ𝑟 + (1.1 − ∅)2ℎ𝑟 + (0.7 − ∅)3ℎ𝑟 + (0.3 − ∅)2ℎ𝑟 = 4.9007 𝑖𝑛
The second trial yields a ∅-index of 0.9999 in/hr ≠ 0.9599 in/hr (from the 1st iteration)
With a ∅-index = 0.9999 in/hr., the equation remains the same which means that
(1.4 − ∅)2ℎ𝑟 + (2.3 − ∅)3ℎ𝑟 + (1.1 − ∅)2ℎ𝑟 + (0.7 − ∅)3ℎ𝑟 + (0.3 − ∅)2ℎ𝑟 = 4.9007 𝑖𝑛
It gives us the same ∅-index of 0.999 in/hr. This is the closest approximation of the
infiltration index ∅.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 43
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
The volume of infiltration losses can easily be estimated by solving the shaded area
under the plot of the ∅-index, and multiplying by the watershed area.
𝑓𝑡 640 𝑎𝑐𝑟𝑒𝑠
𝑉𝑜𝑙𝑢𝑚𝑒 𝑜𝑓 𝑖𝑛𝑓𝑖𝑙𝑡𝑟𝑎𝑡𝑖𝑜𝑛 = 9.70 𝑖𝑛 × × 0.875𝑚𝑖 × = 452.6667 𝑎𝑐. 𝑓𝑡
12𝑖𝑛 𝑚𝑖
To check, the runoff volume should be equal to the shaded area above the ∅-index
line multiplied by the watershed area.
𝐴𝑟𝑒𝑎 𝑎𝑏𝑜𝑣𝑒 ∅ 𝑖𝑛𝑑𝑒𝑥 𝑙𝑖𝑛𝑒 = 1.4 𝑖𝑛 ℎ𝑟 − 1 𝑖𝑛 ℎ𝑟 2ℎ𝑟 + 2.3 𝑖𝑛 ℎ𝑟 − 1 𝑖𝑛 ℎ𝑟 3ℎ𝑟 + 1.1 𝑖𝑛 ℎ𝑟 − 1 𝑖𝑛 ℎ𝑟 2ℎ𝑟 = 4.9 𝑖𝑛
𝑓𝑡 640 𝑎𝑐𝑟𝑒𝑠
𝑉𝑜𝑙𝑢𝑚𝑒 𝑜𝑓 𝑟𝑢𝑜𝑛𝑜𝑓𝑓 = 4.90 𝑖𝑛 × × 0.875𝑚𝑖 × = 228.6667 𝑎𝑐. 𝑓𝑡
12𝑖𝑛 𝑚𝑖
EFFECT OF INFILTRATION
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 44
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
Methods of measuring Infiltration
Flooding type infiltrometers - used to measure the rate of water infiltration into soil
or other porous media
Check out the video on the use of infiltrometer from this link
https://www.youtube.com/watch?v=PYvfTxQhbOQ
Rainfall simulator
Hydrograph analysis
1. A 5-hr storm over a 15-ac basin produces a 5-in. rainfall: 1.2 in./hr for the first hour,
2.1 in./hr for the second hour, 0.9 in./hr for the third hour, and 0.4 in./hr for the last
2 hr. Determine the infiltration that would result from the Horton model with 𝑘 =
1.1/hr, 𝑓 = 0.2 in./hr, and 𝑓 = 0.9 in./hr. Plot the overland flow for this condition in
in./hr vs. t.
2. Parameters in Philip's equation for a clay soil are S = 45 cm-h-1/2 and K = 10 cm/h.
Determine the cumulative infiltration and the infiltration rate at 0.5-hour
increments for a 3-hour period. Plot both as functions of time. Plot the infiltration
rate as a function of the cumulative infiltration. Assume continuously ponded
conditions.
3. For a sandy loam soil, calculate the infiltration rate (cm/h) and depth of
infiltration (cm) after one hour if the effective saturation is initially 40 percent,
using the Green-Ampt method. Assume continuously ponded conditions.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 45
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
4. Compute the ponding time and cumulative infiltration at ponding for a clay
loam soil with a 25 percent initial effective saturation subject to a rainfall intensity
of 1 cm/hr.
Compare the infiltration methods discussed in class, which did you appreciate most?
Share your thoughts in the forum.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 46
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
UNIT 4: EVAPORATION
Pre-assessment Quiz. For the following, write T if the statement is true and F is false.
____ 1. Water lost through evaporation and transpiration is in the form of water vapor.
____ 2. Transpiration occurs only in the daytime.
____ 3. Evaporation continues all day and all night at a constant rate.
____ 4. Evaporation stops when air is fully saturated.
____ 5. Water intercepted by plant leaves that turn to water vapor still comprise
transpiration.
Evaporation is the process by which liquid water from water bodies (rivers, ponds, lakes)
or soil is converted to vapor. It is commonly expressed in mm/h or in/hr.
The rate at which water is converted into vapor is called the rate of vaporization. While
the rate at which water vapor is converted to liquid water, is called the rate of
condensation.
When the vapor pressure is low, vaporization takes place, and so does evaporation. If
the rate of vaporization equals condensation, the air is saturated and evaporation
stops. The difference between saturated vapor pressure (eS) and vapor pressure of the
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 47
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
air (ea) determines the rate of evaporation. As this difference increases, evaporation
also increases, and vice versa.
Factors that affect evaporation in open water include solar radiation, temperature of
the water and air, difference in vapor pressure between water and the overlying air,
and wind speed across the lake.
The amount of water evaporated from a water surface is estimated by the following
methods: (a) using evaporimeter data, (b) empirical formulas and (c) analytical
formulas
(a) Evaporimeters are water-containing pans exposed to the atmosphere wherein loss
of water is monitored at regular intervals. Examples of these are
Each pan has a coefficient given by Cp, so that lake evaporation is estimated as,
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 48
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.
(b) Empirical Equations commonly used meteorological data. Some formulas are
enumerated below:
(c) Analytical Methods include the water budget, energy budget and mass transfer
methods
Between the three approaches, analytical methods can provide good results.
However, they involve parameters that are difficult to assess. Empirical equations can at
best give approximate values of the correct order of magnitude. In view of the above,
pan measurements find wide acceptance in practice.
Accomplish the Summative Quiz for Module 1 and the Midterm Examination.
Property of and for the exclusive use of SLU. Reproduction, storing in a retrieval system, distributing, uploading or posting online, or transmitting 49
in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise of any part of this document, without the prior written
permission of SLU, is strictly prohibited.