Soil–Water–Plant Relationships
Where crops get their water. Learn soil texture and structure, bulk density and porosity, gravimetric and volumetric moisture, soil water potential and the retention curve, saturation, field capacity, permanent wilting point, available water and management allowed depletion, infiltration and hydraulic conductivity, root zones, evapotranspiration basics, soil salinity and measuring soil moisture (gravimetric, tensiometer, TDR/capacitance sensors) — with Python for available water and irrigation interval.
Soil is the crop's water reservoir. How much water it can hold between field capacity (after drainage) and wilting point (plants can no longer extract water) decides how often to irrigate.
Understanding soil–water relationships is the basis of irrigation scheduling, drainage design and efficient water use.
An analogy
Soil is like a sponge in a sink: when soaked it drips (gravity drainage) until it reaches field capacity; plants then suck water out until what remains is held so tightly they can't pull it (wilting point). The water between is the plant's 'bank balance'.
Soil physical properties
Texture (sand, silt, clay proportions — USDA textural triangle), structure (aggregation), bulk density (1.1–1.7 g/cm³; compaction raises it and restricts roots), porosity (~40–55%). Clay soils hold more water but drain slowly; sandy soils drain fast and hold little.
Soil water constants
Key moisture levels:
- Saturation — all pores filled
- Field capacity (~ −33 kPa) — after free drainage (1–3 days)
- Permanent wilting point (~ −1500 kPa) — plants cannot recover
- Available water = FC − PWP (sand ~60 mm/m, loam ~150 mm/m, clay ~200 mm/m)
- Management allowed depletion p ≈ 0.3–0.6 depending on crop sensitivity
Infiltration and movement
Infiltration rate starts high and decreases to a steady basic rate (sand > 30 mm/h, clay < 5 mm/h). Water moves by Darcy's law with unsaturated hydraulic conductivity K(θ) that drops sharply as soil dries (Richards' equation). Infiltration governs choice of irrigation method and application rate to avoid runoff.
Measurement and salinity
Measure moisture by gravimetric oven drying (reference), tensiometers (suction), TDR/FDR/capacitance probes, neutron probes, and satellite soil moisture (SMAP, Sentinel-1). Salinity (ECe) reduces water availability via osmotic potential; leaching requirement LR = EC_w/(5EC_e − EC_w) guides extra irrigation. Sodic soils need gypsum and drainage.
Formulas
- Volumetric moisture:
θ = w·ρ_b/ρ_w - Porosity:
n = 1 − ρ_b/ρ_s (ρ_s ≈ 2.65 g/cm³) - Total available water:
TAW = 1000·(θ_FC − θ_PWP)·Z_r (mm) - Readily available water:
RAW = p·TAW - Irrigation interval:
interval = RAW / ET_c - Infiltration (Kostiakov):
I = k·t^a
Worked example: Volumetric moisture
Given: w = 0.18, ρ_b = 1.5
Find: θ.
- 0.18 × 1.5
Answer: 0.27
Worked example: Available water
Given: FC 0.30, PWP 0.12, root depth 0.9 m
Find: TAW.
- 0.18 × 0.9 × 1000
Answer: 162 mm
Worked example: Irrigation depth
Given: θ now 0.22, FC 0.30, root 0.6 m
Find: Net depth to refill.
- 0.08 × 600
Answer: 48 mm
Available water and irrigation interval
# Soil-water-plant: available water, bulk density and moisture
FC, PWP = 0.32, 0.15 # volumetric water content at field capacity and wilting point
root_m = 0.6
TAW = (FC - PWP) * root_m * 1000 # mm total available water
p = 0.5 # allowable depletion fraction
RAW = p * TAW
print(f"total available water {TAW:.0f} mm, readily available {RAW:.0f} mm")
# Gravimetric -> volumetric
w, rho_b = 0.20, 1.40 # gravimetric moisture, bulk density g/cm3
theta = w * rho_b
porosity = 1 - rho_b / 2.65
print(f"volumetric moisture {theta:.2f}, porosity {porosity:.2f}, saturation {theta/porosity:.0%}")
# Days between irrigations if crop uses 6 mm/day
print(f"irrigation interval at ETc 6 mm/day: {RAW/6:.1f} days")
# total available water 102 mm, readily available 51 mm
# volumetric moisture 0.28, porosity 0.47, saturation 59%
# irrigation interval at ETc 6 mm/day: 8.5 daysA loam with 0.6 m roots stores ~102 mm of available water; allowing 50% depletion before stress gives ~51 mm, i.e. irrigate about every 8 days in peak demand. Sandy soils hold less and need more frequent, smaller irrigations.
Common mistakes
- Wrong: Using gravimetric moisture as volumetric. — Right: Multiply by bulk density.
- Wrong: Applying water faster than infiltration. — Right: Match application rate to soil intake.
- Wrong: Ignoring salinity build-up. — Right: Include a leaching fraction and drainage.
Applications
- Irrigation scheduling: When and how much to irrigate.
- Sensor-based automation: Soil-moisture triggered drip.
- Drainage design: Waterlogged and saline lands.
- Soil health: Compaction and organic matter management.
Key takeaways
- θ = w·ρ_b; n = 1 − ρ_b/2.65.
- TAW = (FC − PWP)·Z_r; RAW = p·TAW.
- Interval = RAW/ET_c.
- Match application to infiltration; manage salinity.
Going further
- Richards' equation: Unsaturated flow modelling (HYDRUS).
- Pedotransfer functions: Estimating hydraulic properties from texture.
- Remote-sensing soil moisture: SMAP and Sentinel-1 retrievals.