LWC Lecture Notes
LWC Lecture Notes
ABE 4237
SOIL AND WATER CONSERVATION ENGINEERING - Agriculture accounts for between 60-80% of total water use
wherein about 40% of the water for irrigation was not
- The application of engineering principles to the solution of soil available to crops
and water management. - In areas where water is limited, soil-water should be
- Conservation practices can be beneficial; however, it can also conserved by modified tillage and crop management
be detrimental to the environment techniques, level terracing, contouring, reservoirs, water
- Soil and water resources is a worldwide problem (Napier et al., harvesting, and other means to reduce evaporative losses.
2000) especially in developing countries where people are
more concerned with survival than conservation (Troeh, F.R,
et al., 2004)
Two ways in Increasing Crop Production
Importance of Soil and Water Conservation Engineering
- To develop new land
Soil – commonly referred to as earth or dirt, a mixture of - To improve the productivity of present cropland
organic matter, mineral, gases, liquids, and other organisms
that support life.
Water – primary component of the biosphere, that is vital for Roles of Agricultural Engineers in Soil & Water Conservation
agriculture, domestic use, recreation, wildlife, and others.
- To develop and execute a conservation plan:
*Soil and water is vital for sustaining life and crucial as the soil – plant – water – environment interactions
population is increasing. - Must be familiar on the social economic backgrounds
relating to soil and water conservation
- Have a full understanding and updated on the various
Engineering Problems government structures and mechanisms that have been
developed in implementing soil & water conservation
Erosion Control programs
- Primary caused by human activities by the removal of - Should be familiar with the principles of mapping and
natural vegetation (urban-rural areas) classifying land use (GSIS), and weather records and
- Important in the maintenance of crop yields (in the USA predictions
the amount of soil eroded annually is estimated to be 1M –
200,000 ha of farmland) Soil and Water Conservation around the World
- Caused by water and wind (climate, soil type, vegetative
- Many countries have developed a conservation energy like
cover, topography, conservation practice)
the Natural Resources Conservation Services in the United
*Erosion Control Practices States within a national department or ministry of
agriculture. There have been good results however, there are
1. Planting, tillage, and harvesting operations on or early still some inconsistency and little information to be done to
on the contour reduce soil deterioration and water losses
2. Planting of intertilled crops in alternate strips - Transfer of conservation technology problems include
3. Constructing of cross-slope channels geographic and cultural differences
4. Trees barriers - Reduce the chlorofluorocarbons (CFCs) and greenhouse
5. Mulching gases – the Paris Agreement
6. Vegetative waterways - Shifting cultivation & Conservation – slash-and-burn
7. Stabilizing gullies
*Irrigation
1. Collect daily rainfall data for a period of ten years or more. For 1. Shallow well areas – 5.1M ha.
each year, get the weekly total rainfall and arrange as WK1
(Jan 1-7), WK2 (Jan 8-14) and so on until WK 52 (Dec 24-31). 2. Deep well areas – 12.3M ha.
2. Decide what a Dry (D) or Wet (W) week is by selecting a
Water Resources
threshold value of rainfall for a dry/wet week. A week is
considered a dry week when rainfall is less than 20mm, and/or Four largest river basins
less than 0.1mm rainfall in a day
3. For each year; 1. Cayagan River Basin (Region 2) – damage area = 25,649 km2
a. Identify the D and W weeks
b. Count the frequency or the number of dry weeks F(D) 2. Mindanao River Basin (Region 12) – DA =23,169 km2
and wet weeks F(W).
c. Count the frequency or the number of times a Dry week is 3. Agusan River (region 10) - DA = 10,921 km2
preceded by another dry week F(D/D). 4. Pampanga River Basin (Region 3) = DA = 9245 km2
d. Count the frequency or the number of times a Wet week
is preceded by another Wet week F(W/W Smallest
4. Calculate and tabulate the following probabilities for each
week of each year: 1. Buayan-Malungon River Basin – DA=1434 km2
a. Initial probabilities:
i. Probability of a dry week, P(D) P(D) = F(D)/n Major Lakes
ii. Probability of a wet week, P(W) P(W) = F(W/n) 1. Laguna de Bay and Taal Lake in Southern Luzon
Dams (storage capacity exceeds 50M m3 and structural height is 4. Small Water Impounding Projects (SWIPs)
more than 30m).
The Philippines is not big and has never been any small that its
country’s land and water resources has to support a population
growing at an average annual rate of 2.3%.
Example:
1. Temporary Structures
2. Permanent Structures
Temporary Structures
Permanent Structures
Specific Energy
The total energy of water with respect to the channel bottom. The
discharge of water in an open channel may flow at two depths
known as alternate depth with the same total energy head. Specific
energy is sometimes used than total energy (Henderson, 1966).
Where:
Capacity
Outlet Protection
1. Suitable wing walls Used as a culvert has the simple function of providing for
passage of water under an embankment.
2. Side walls
Serves to lower water through a considerable change in
3. Head wall extension elevation and to decrease the energy of the falling water
when combined with a riser or drop-inlet.
4. Toe walls
Drop-inlet is used as gully control structures
Three Types of Permanent Structures Withstand hydraulic pressures and require special
attention to structural design (inverted siphon)
1. Drop Spillway
Its capacity structures are relatively low
2. Drop inlet (Pipe spillway) Only desirable when the discharge from the structure is to
be restricted.
3. Chute Spillway
Culverts
Is a weir structure which flow passes through the weir opening
and falls on an apron and then passes into the downstream It may be controlled by inlet section or conduit, depends
channel. on the type of flow.
Installed in channels to establish a permanent control elevation Headwater may be above or below the top of the inlet
below which an eroding stream cannot lower the channel. section
It controls the stream grade from the spillway crest through the Pipe flow will occur when the slope is less than neutral
entire ponded reach stream. slope and entrance capacity is not limiting
Drop limited at 3m. Flumes or drop-inlet pipe spillways are If the inlet is submerged and outlet is not, the slope of the
used for greater drops. conduit is greater than neutral slope, then orifice flow will
prevail.
The drop spillway capacity is given by the weir formula:
Erosion
Discharge capacity for Orifice flow
A geomorphic process that detaches and removes materials from its
q= AC √ 2 gH primary location by natural erosive agent.
Where: Mechanics of Soil Erosion
q = discharge capacity The soil erosion process undergoes three main processes
(mechanisms) such as the detachment, transport, and deposition.
A = conduit cross-sectional area
Three main mechanisms of Soil Erosion
H = head causing flow, 0.6
1. Detachment (without it, erosion will not even start)
C = orifice discharge coefficient, 0.6 for a sharp-edged
2. Transport (without it, erosion will be much limited
3. Deposition (there will be no deposited soil particle/materials)
Note: if an orifice plate is installed in a pipe the differential head 1. Idle bare land areas are prone to wind and rain erosion
measurements are used to calculate the flow rate, however, if an 2. Growing of more crops accelerate soil erosion
orifice is configured as an end-cap, a head measurement taken 3. Cultivation along slope of land
about 0.6 m upstream can be used. 4. Faulty and inefficient methods of irrigation
5. Removal of organic matter and plant nutrients by inappropriate
SAMPLE PROBLEM cropping pattern
A 457 mm diameter that is 20 m long galvanized iron culvert with a 6. Deforestation
square-edged entrance will be using in constructing a control 7. Forest fire
structures. This have an inlet elevation of 150 m, outlet invert 8. Overgrazing (pasture-patches)
elevation of 149.10 m, headwater elevation of 152.30 m, and the Types of Soil Erosion
tailwater elevation is 148.50 m. Determine the capacity of a full
pipe flow and compared the value of actual slope to the neutral 1. Geological Erosion
slope. 2. Accelerated Erosion
Geological Erosion
Accelerated Erosion
Sheet Erosion Soil erosion by water may be grouped in those that are due to:
o Uniform removal of soil in thin layers resulting from sheet or 1. The action of rainfall (splash erosion)
overland flow. 2. Detaching and transporting actions of runoff (overland flow
o During overland flow, detached soil particles in the flowing scour, rill, gully, and channel erosion)
water will scour other soil particles and transport them further Agronomical Measures to control water erosion
down slope.
o uniform removal of soil in thin layers resulting from sheet or Soil conservation is a preservation technique in which the soil
overland flow. disturbance and losses are minimized by applying conservation
o Its soil detaching and transporting capabilities are small techniques as well as improvement of the soil (or by using it within
its capabilities). Agronomical measures is one of the more m1 and m2 = exponents with estimated values of 1.4 and 0.6,
economical, long lasting, and effective technique. This type of respectively
conservation measures function by reducing the impact of
raindrops, thus reducing soil erosion as well as the surface runoff Ellison (1945 and 1947) – advocated the development of
and enhance infiltration rates. erosion research on splash and erosion. He suggested that
splash erosion is a function of soil and rainfall properties
Widely used agronomical measures for water erosion control
4.33 1.07 0.65
1. Contour cropping – used in sloppy condition to control soil SS= A V D I
losses due to water erosion) Where:
2. Strip cropping (contour strip cropping) – alternate
cropping, best potential for erosion control SS = soil splash in grams in 30 minutes
3. Mulching – minimize rain splash, decrease evaporation, A = soil constant
control weeds, reduce temperature. I = intensity of rainfall in inches per hour, and
D = drop diameter in mm.
Wind Erosion
1. Surface creep (7-25%) – larger diameter soil particles (0.2- Musgrave equation – known as the slope-practice equation
2mm) moves very close to the ground surface and roll along. which was formulated by a committee which was formed at
(collision) about the same time to integrate the results of previous studies.
2. Saltation – (55-72%) – fine textured soil particles (0.1 –
m1 m2 m3
0.5mm) bouncing over the surface causing abrasion and E=0.0053 IRS L P30
attrition or an avalance effect. (breaking of particles into
smaller particles). Where:
3. Suspension (3-10%) – transport of finer particles (0.02-0.1
E = soil loss in mm per year
mm) in the atmosphere forming dust storm due to turbulent
action of the wind. I = inherent erodibility of a soil at 10% slope and 22m slope length
Three phases of wind erosion R = vegetal cover factor
1. Initiation of movement and air velocity and turbulence S = slope in percent
2. Transportation and particle size, wind velocity and range
of sizes L = slope length, m
3. Deposition occurs when weight>others and lower wind
P30 = maximum rainfall in mm over a 30-minute period m1, m2,
velocity, barrier etc
and m3 = constant exponents
Factors affecting Erosion by Wind
Smith and Whitt (1947 and 1948) – developed an expression
1. Climate – precipitation, temperature, wind, humidity, and for the claypan soils of Missouri, expressed sheet erosion as:
the air density
𝐸 = 𝐶𝑆𝐿𝐾𝑃
2. Soil – texture, structure, density of particles, density of
soil mass, organic matter, moisture content. Where:
3. Vegetation – height, density of cover, type of vegetation,
seasonal distribution. E = average annual sheet erosion loss
4. Topography – degree of slope, slope length, size and
C = average annual rotation soil loss from plots
shape of watershed
S,L,K and P = factors to adjust the plot soil loss C for slope
Soil Loss Estimation
steepness, slope length, soil group and conversation practices,
Soil losses vary considerably with the type of erosion. respectively.
There are several investigation and studies that have attempted
to relate the physical characteristics of rainfall (e.g. raindrop
size, and kinetic energy) and soil characteristics to soil splash. Note: The above equations are of empirical nature and were found
Others include the effect of cover, winds, slope, surface to be too localized or regional in applicability. As a result, the
configuration and all other sub-processes as identified by erosion research efforts in the US were consolidated starting in
Ellison (1945). 1954 (more erosion data were collected and all available
Several methods have also been proposed for estimating information were compiled).
erosion by water, however, there’s no suitable method found
for predicting wind erosion losses.
Most of the soil erosion equations were developed in US, The Universal Soil Loss Equation
basically, splash and sheet erosion estimation.
First equation proposed by Zingg (1940) who expressed soil Developed in the period from 1945 until 1965 to estimate
loss from cropped areas as follows: losses based on statistical analyses of field plot data from
small plots.
m1 m2 Gained wide acceptance in the United States and has been
E=C S L
adopted by soil conservationists in other countries
Where: Its popularity is due to its simplicity and the absence of
alternative equations allowing a wider range of variations of
E = soil loss per unit area from a land slope of unit width
the parameters that are contributing to soil erosion
C = constant reflecting the combined effect of rainfall, soil cover
and cover management practices The USLE is very simple to use
S = slope
L = horizontal length of slope and
It estimates sheet erosion as the product of a series of terms for Important elements of it are the quantity (duration and
rainfall, soil, slope gradient, slope length, crop and cover distribution) and intensity of precipitation
management and conservation factors
Energy of Rainfall
This permits the extensive tabulations of individual factors,
incorporating vast amounts of research data. KE of falling raindrops is 0.5mV2
In fact, extensive tables and nomographs recommending the Erosive energy of raindrop is reached after 5m fall
range of values for the different field conditions for the various
Maximum energy is reached after a 20m
parameters in the equation are now available (including from
the Philippines) Quantity of Rainfall
The USLE was developed from standard agricultural plots
which has recommended dimension of 72.6 feet long and 6 The higher the rainfall, the higher the erosion hazard
feet wide, with 9% slope bare and cultivated fallow, and
plowed up and down the slope.
The validity of the USLE depends on how close the conditions Intensity of Rainfall
from which it is being applied at are to the range of conditions
30-360 mm/hr intensities ~ 10% of rain is erosive
from which it was developed.
(approx.)
However, applying the equation in areas with conditions very
>100 mm/hr, all the rain is erosive
much different from those where it was developed or in areas
5% is erosive in temperature; 40% in tropical
with somewhat similar conditions but with factor values not
yet determined is somewhat dangerous. Calculating Rainfall Erosivity
In areas where the USLE applies, it serves as a useful tool
for soil conservation. Among these are: rainfall erosivity is related to the kinetic energy of
1. Predict soil erosion losses rainfall.
2. Guide the selection of agricultural practices such as two methods are widely used for computing the erosivity
cropping and management systems of rainfall: EI30 Index Method and KE>25 Index Method.
3. Guide in determining the on-site effects of land use and 1. Calculating Rainfall Erosivity
crop management changes. EI30 Index (Wischmeier et. al., 1958) found that rain KE is
4. Provide baseline data for conservation planning the factor most closely related to erosion after analyzing 8250
plot-year records from 35 conservation stations in the US.
• Based on the fact that the product of KE of the storm
and the 30-min max rainfall intensity gives the best soil
USLE was first developed by by Wischmeier and Smith (1965).
loss estimation
A=RKLSCP • Rainfall erosivity factor value can be computed using:
Where: E I 30=KE x I 30
A = average annual soil loss, t/ha/yr Where:
R = rainfall erosivity factor, mm KE = rainfall kinetic energy
K = soil erodibility factor, t/ha/hr/in
I 30= maximum rainfall intensity for a 30-minute period.
L = slope length factor, dimensionless
*Kinetic energy can be computed using:
S = slope steepness factor, dimensionless
KE = 0.119 + 0.0873log I
C = cover management factor. dimensionless
Where:
P = conservation practice factor, dimensionless
KE = kinetic energy, MJ/ha-mm
USLE was modified by William in 1975 to MUSLE by replacing
the rainfall energy factor (R) with another factor called as runoff I = intensity of rainfall, mm/h
factor. USLE was also modified by David (1986), and David and
Collado (1987) to suite locally available information and prevailing Note: I 30 index method was developed under American
environmental conditions in the absence of any applicable method condition and is not suitable for tropical and sub-tropical
for estimating soil erosion rates. zones for estimating the erosivity.
E=RKSCP
( )
m
Pi = daily rainfall total greater than 25mm x
L=
A = 0.002 (to be compatible with R using USLE 22.13
m=2
Soil Erodibility Factor (K) Where:
It is the rate of soil loss on a standard plot of 9% slope, 72.6 ft L = length factor
long. x = slope length, m
Few K values were obtained from studies on the soil m = an exponent = 0.5 for S ≥ 5
conservation experiment from follows plots and from row- = 0.4 for 5 < S > 3
crops, corrected for vegetative cover while many were = 0.3 for 3 ≤ S ≥ 1
obtained from small plots using a rainfall simulator. = 0.2 for S < 1
Soil that are high in silt, low in clay and low in organic matter S = percent slope
content are the most erodible
Soil becomes less erodible with decrease in its silt fraction,
regardless of whether sand or clay increases Combining the effects of slope and slope length yields:
OM improves soil aggregation and the infiltration capacity of
soils
( )
m
x 2
L= (0.065+ 0.045 S +0.0065 S )
22.13
This equation was normalized to a standard plot 22.12 m
in length with 9% slope
As USLE gains wide acceptance and applications, so does
the LS factor.
m
LS=a+ bS
Where:
LS = slope-slope length factor
S = land slope, %
a, b = constant equal to 0.1 and 0.21, respectively
m = exponent equal to 4/3
Effect of Cover
1. Interception of rainfall
2. Decrease runoff velocity and the eroding action of water
3. Root effects in increasing granulation and porosity
4. Vegetative biological activities – influence on soil
porosity
5. Transpiration of water – subsequent drying out of the soil
Vegetated Waterways
Are constructed to protect the soil from erosive resistant
of runoff from sloping lands in order to prevent the
formation of rill to gully channel.
May be used route storm runoff from structures including
building roofs, highways. (may be used to lessen the
impact of storm to runoff)
Filtered the water passing through the soil surface.
More complex design as compared to the concreted
channel line.
Design consideration
1. Design flow – dependent on the type of channel,
condition and density of vegetation, and erosive factor.
2. Soil properties – texture, structures, porosity.
3. Topography
4. Channel capacity – based on the estimated runoff – 10
year storm record
5. Channel shape
A. Parabolic – usually selected for channels in
natural waterways
B. Trapezoidal is most easily constructed on
artificial waterways that is in a terrace outlet
along a fence line.
C. Triangular reduce sedimentation but with a
high velocity can damage the bottom of the
waterway.
Runoff Estimation
CIA
Q p=
360
Where:
Qp = peak runoff rate, m3/sec
C = coefficient
I = intensity of rainfall, mm/hr
SEDIMENT TRANSPORT AND SEDIMENTATION
A = area, hectare
Example Problem:
1. A concrete lined canal will be conveying a water from the well to
a tomato farm with a distance of 150 meters. It has a design of 0.5
m drop elevation over the 150 m canal length and roughness
coefficient of 0.018. Determine the velocity of water given the
bottom width of a rectangular canal of 0.5 m with a depth of 0.25
Load
m.
one of the most important sediment of stream flow that
moves along the channel bed and in almost contact with
2. The peak runoff rate for a 0.5 cm/hr rain intensity at 267 ha
the stream bed.
watershed was observed as 2.5 m3/sec. Estimate the runoff
coefficient. Represents the armor that resists erosion
Suspended Load
is carried in the fluid away from the bed.
Materials that are carried by the turbulent forces of the
flowing water
Coarse sediment: Particles diameter > 0.20 mm
Medium sediment: Particles diameter between 0.20 to
0.075mm
Fine sediment: Particle diameter < 0.075mm
Sediment Load of a Stream as per Grain Size Soil Loss, t/ha for 18 River Basins in the Philippines
Wash load
The fine fraction which is mostly in suspension
Bed Material Load
coarse fraction that normally found in the bed
Estimation of Sediment Discharges
Actual monitoring of the sediment load through sediment
sampling
Using sediment discharge formulas
Pumping samplers
Use of photocell to measure the attenuation of light beam
by the suspended particles
Use of x-ray probe
Terrace Layout
1. Select suitable outlet – such as permanent
features, densely wooded areas
2. Locate terraces – usually lay out top terrace first
a. Drainage area above top terrace should be
less than 3 acres
b. If the top of the hill comes to a point, the VI
may be increase to 1.5 times the regular VI.
c. If short abrupt changes in slope occur, the top
terrace should be placed just above the break
a. Stake out top terrace using 0.3% for entire line of terrace
b. Determine spacing for each of the remaining terraces
• Mark each terrace with single stakes
c. Select key terrace
• criteria for selection: other terraces may be made parallel
to it with minimum correction
d. Stake out key terrace (0.3% grade)