Drainage Basin Hydrology Advanced

Spec mapping: AQA 7037, Paper 1 (Physical), §3.1.1 — the drainage basin hydrological cycle as an open system (inputs, outputs, stores, flows), the water balance, and the factors influencing change in the cycle over time and space. This depth lesson assumes the basic open-system model and pushes into the physics of infiltration, the Thornthwaite soil-moisture budget, antecedent conditions, baseflow separation and recession analysis. AOs exercised: AO1 (precise process knowledge — Horton infiltration, field capacity, the baseflow index), AO2 (explaining how catchment properties translate inputs into a runoff response) and AO3 (constructing a soil-moisture budget, separating baseflow, solving the water-balance equation). Synoptic links run to Hazards (flood generation) and the storm-hydrograph lesson that follows.

A drainage basin is the area drained by a river and its tributaries, bounded by the watershed. Unlike the global cycle, it is an open system: water crosses the boundary as input (precipitation) and output (river discharge, evapotranspiration, deep seepage). The art of advanced basin hydrology is to treat the basin as a transfer function — a machine that converts a precipitation input into a discharge output, with the conversion governed by stores (interception, soil moisture, groundwater) that delay, store and release water. Mastery of the quantitative behaviour of that machine is what distinguishes top responses.

---

The Drainage Basin as an Open System

flowchart TB
  P[Precipitation INPUT] --> INT[Interception store]
  P --> DIR[Direct channel
precipitation]
  INT -->|throughfall + stemflow| SURF[Surface]
  SURF -->|infiltration| SM[Soil moisture store]
  SURF -->|overland flow| CH[Channel store]
  SM -->|throughflow| CH
  SM -->|percolation| GW[Groundwater store]
  GW -->|baseflow| CH
  CH --> Q[River discharge OUTPUT]
  SM -->|evapotranspiration| ET[Atmosphere OUTPUT]
  INT -->|evaporation| ET
  GW -->|deep seepage| OUT[Beyond watershed OUTPUT]

Inputs, outputs, stores and transfers

• Inputs: precipitation (rain, snow, sleet, hail); solar energy driving evapotranspiration and snowmelt.
• Outputs: river discharge $Q$Q (measured at the gauging station in cumecs, m³/s); evapotranspiration $E$E; deep groundwater seepage across the boundary (often small in short-term analysis but not always).
• Stores: interception (~0.5–3 mm per storm in deciduous woodland, much more in conifers), surface storage (puddles, ponds, wetlands), soil moisture, groundwater, channel storage, and seasonal snow/ice.
• Transfers (flows): infiltration (surface into soil), percolation (soil/rock down to the water table), throughflow (lateral, often along permeability contrasts), groundwater/baseflow (slow lateral flow through saturated rock to the channel), overland flow (surface runoff), stemflow and throughfall (intercepted water reaching the ground), and channel flow.

The teaching power of the open-system model is that every store imposes a delay. Water routed through the groundwater store may take months to reach the channel; water that runs off the surface arrives in hours. The partitioning of input between fast and slow pathways therefore sets the entire hydrological character of the basin — and that partitioning is decided first at the soil surface, by infiltration.

---

Infiltration and Infiltration Capacity

Infiltration is the entry of water into the soil; the maximum rate at which a soil can absorb water is its infiltration capacity (mm/hr).

Horton's infiltration model

Robert Horton (1933) showed that infiltration capacity is highest at the start of a storm and decays exponentially as pores fill, towards a constant final rate:

$f(t) = f_{c} + (f_{0} - f_{c})\,e^{-kt}$f(t)=fc​+(f0​−fc​)e−kt

where $f(t)$f(t) is the infiltration capacity at time $t$t, $f_{0}$f0​ is the initial capacity, $f_{c}$fc​ is the final (equilibrium) capacity — essentially the saturated hydraulic conductivity — and $k$k is a soil-specific decay constant.

• $f_{0}$f0​ ranges from >200 mm/hr for well-structured sandy soils to <5 mm/hr for compacted clay.
• $f_{c}$fc​ is the long-run rate once the soil is wetted up; it is what matters in prolonged rainfall.

The mechanism of decay matters: as the wetting front advances, the suction (matric potential) gradient that pulls water in weakens, and pores progressively fill, so gravity alone drives the residual flux at $f_{c}$fc​.

Factors affecting infiltration capacity

Factor	Effect
Soil texture	Large pores in sand → high rates; tiny pores in clay → low rates
Soil structure	Macropores (worm burrows, root channels, cracks) allow rapid bypass flow
Antecedent moisture	Wetter soil = less available pore space = lower capacity
Vegetation cover	Roots maintain structure; litter shields the surface from raindrop impact
Land use	Urban sealing and agricultural compaction sharply reduce capacity
Surface crusting	Raindrop impact on bare soil seals the surface
Slope angle	Steeper slopes give water less residence time to infiltrate
Frozen ground	Concrete frost can render even sandy soil near-impermeable

Key distinction. Infiltration-excess (Hortonian) overland flow occurs when rainfall intensity exceeds infiltration capacity — typical of intense convective storms on crusted or compacted surfaces. Saturation-excess overland flow occurs when the soil is saturated from below, so any further rain runs off regardless of intensity — typical of valley bottoms, footslopes and areas of shallow water table. Which mechanism dominates controls where in the basin runoff is generated (the variable source area concept: the contributing area expands during a storm as saturated zones spread up the valley).

---

Soil Moisture Dynamics

Key thresholds

• Saturation: all pores water-filled.
• Field capacity: moisture remaining after gravitational (free) water has drained, usually 1–2 days after rain. The soil is moist but not waterlogged; water is held at ~ −0.033 MPa.
• Wilting point: the moisture content (~ −1.5 MPa) below which plant roots can no longer extract water; the remainder is held too tightly in micropores.
• Available water capacity (AWC): field capacity minus wilting point — the water actually available to plants. AWC is high in loams (good balance of pore sizes), low in both sands (drains too freely) and clays (holds water too tightly).

The Thornthwaite soil-moisture budget

C. W. Thornthwaite's (1948) budget is the standard A-Level framework for tracking soil moisture through the year by comparing precipitation ($P$P) with potential evapotranspiration ($PE$PE, the atmospheric demand if water were unlimited):

• When $P > PE$P>PE, the surplus first recharges the soil to field capacity, then generates runoff/percolation (soil-moisture surplus).
• When $PE > P$PE>P, plants draw on stored soil water; once that is exhausted, actual evapotranspiration (AE) falls below PE and a soil-moisture deficit (SMD) opens up. AE < PE is the hallmark of a supply-limited period.
• Utilisation is the drawdown of stored water; recharge is its replenishment in autumn.

A schematic annual budget for a temperate-maritime (UK) site:

Month (UK)	$P$P vs $PE$PE	Soil-moisture state	Dominant process
Dec–Feb	$P > PE$P>PE	At/above field capacity	Surplus → runoff, recharge
Mar–Apr	$P \approx PE$P≈PE	Near field capacity	Transition; recharge ends
May–Jul	$PE > P$PE>P	Falling; deficit develops	Utilisation; AE < PE
Aug	$PE > P$PE>P	Maximum deficit	Peak SMD; plant water stress
Sep–Oct	$P > PE$P>PE	Deficit shrinking	Recharge begins
Nov	$P > PE$P>PE	Back to field capacity	Recharge complete; surplus resumes

Soil-moisture deficit (SMD) is the depth of water needed to restore the soil to field capacity. In south-east England SMD can reach 100–120 mm in a dry summer. SMD has direct consequences for irrigation demand, crop yield, and — counter-intuitively — flood risk: a severely dried soil can become hydrophobic (water-repellent), especially in peaty or organic soils, so the first autumn storms may run off rather than infiltrate.

---

Antecedent Moisture and Its Significance

Antecedent moisture is the water already present in the soil before a storm, and it is among the most powerful controls on hydrological response.

High antecedent moisture (e.g. a saturated catchment after a wet fortnight):

• reduced infiltration capacity (pores already occupied);
• rapid generation of saturation-excess overland flow;
• higher and faster peak discharge; shorter lag time;
• a large proportion of rainfall becomes quickflow.

Low antecedent moisture (e.g. after a dry summer):

• high initial infiltration capacity;
• possible hydrophobicity in very dry organic soils, paradoxically increasing early runoff;
• lower peak discharge; longer lag;
• much of the rainfall absorbed into soil-moisture recharge rather than runoff.

This is why two meteorologically identical storms can produce utterly different floods: the same 80 mm of rain falling on a saturated November catchment versus a parched August catchment generates very different hydrographs. Antecedent conditions are the single most examinable reason for that contrast (and a recurring theme in the case studies of the next lesson).

---

Baseflow and Quickflow Separation

River discharge is conventionally separated into two components.

Baseflow: groundwater discharge to the channel — sustained, slow, relatively constant; it maintains flow in dry spells and dominates in permeable catchments (chalk, limestone, sandstone).

Quickflow (stormflow): the rapidly arriving component — overland flow, rapid throughflow and direct channel precipitation; it forms the hydrograph peak and dominates in impermeable, steep or saturated catchments.

The Baseflow Index (BFI)

The BFI is the long-term proportion of total runoff supplied by baseflow:

$\text{BFI} = \frac{\text{baseflow volume}}{\text{total runoff volume}}$BFI=total runoff volumebaseflow volume​

• Chalk catchments (e.g. River Itchen, Hampshire): BFI ≈ 0.90–0.95 — almost all flow is baseflow, so the river is flashy-resistant and reliable in drought.
• Clay catchments (e.g. upper Severn headwaters): BFI ≈ 0.30–0.45 — large quickflow contribution, flashy response.
• Urban catchments: BFI can fall below 0.20.

The BFI is a single number that captures a catchment's "personality" and is widely used by the UK National River Flow Archive. It links geology directly to flood and drought behaviour — a permeable, high-BFI basin floods rarely but sustains summer flow; an impermeable, low-BFI basin floods readily but can run very low between storms.

---

Recession Curves and Baseflow Separation as a Skill

A recession curve (the falling limb after rain ceases) describes the draining of the catchment's stores. It is approximately exponential:

$Q_{t} = Q_{0}\,e^{-\alpha t}$Qt​=Q0​e−αt

where $Q_{t}$Qt​ is discharge at time $t$t, $Q_{0}$Q0​ is discharge at the start of recession, and $\alpha$α is the recession constant. A steep recession (large $\alpha$α) indicates rapid drainage from an impermeable or urban catchment; a gentle recession (small $\alpha$α) reflects the slow release of groundwater from a permeable catchment — a chalk stream's recession can extend for weeks.

AO3 worked exemplar: separating baseflow and reading recession

Suppose a gauging station records the following daily mean discharge through a storm:

Day	Discharge $Q$Q (m³/s)	Phase
0	4	Pre-storm baseflow
1	9	Rising limb
2	22	Approaching peak
3	28	Peak discharge
4	19	Falling limb
5	12	Falling limb
6	8	Late recession
7	6	Approaching baseflow

Describe. Discharge rose from a baseflow of 4 m³/s to a peak of 28 m³/s on day 3, then fell back towards ~6 m³/s by day 7.

Manipulate. The storm's quickflow contribution can be approximated by subtracting an assumed straight baseflow line (rising gently from 4 to ~6 m³/s) from the total. The peak quickflow on day 3 is roughly $28 - 5 = 23$28−5=23 m³/s. As a percentage increase above pre-storm flow:

$\frac{28 - 4}{4} \times 100 = 600\%.$428−4​×100=600%.

Discharge rose six-fold. The recession constant between days 4 and 6 can be estimated from $Q_{4}=19$Q4​=19 and $Q_{6}=8$Q6​=8:

$\alpha = \frac{\ln(Q_{4}/Q_{6})}{t} = \frac{\ln(19/8)}{2} = \frac{0.865}{2} \approx 0.43\ \text{day}^{-1}.$α=tln(Q4​/Q6​)​=2ln(19/8)​=20.865​≈0.43 day−1.

Explain. The six-fold, two-day rise and the fairly steep recession ($\alpha \approx 0.43$α≈0.43) together indicate a flashy, low-BFI catchment — most of the storm water travelled by quickflow, with little groundwater buffering.

Evaluate. Straight-line baseflow separation is an approximation; in reality baseflow may rise during the event (bank storage) and the "true" separation is ambiguous. Daily means also hide the true instantaneous peak. The analysis is indicative, not exact — and saying so is an AO3 evaluation mark.

---

The Water Balance Equation

The master equation of basin hydrology:

$P = Q + E \pm \Delta S$P=Q+E±ΔS

with $P$P precipitation, $Q$Q runoff, $E$E evapotranspiration and $\Delta S$ΔS the change in storage (positive = filling, negative = depleting), all in mm or m³.

Worked example

A basin receives 1,200 mm of precipitation in a year; evapotranspiration is 520 mm and discharge at the outlet totals 640 mm. Then:

$\Delta S = P - Q - E = 1{,}200 - 640 - 520 = +40\ \text{mm}.$ΔS=P−Q−E=1,200−640−520=+40 mm.

The basin gained 40 mm of storage — groundwater levels rose or soil moisture increased. Over a long enough record $\Delta S \to 0$ΔS→0, so a persistent positive or negative residual signals either measurement error or genuine long-term storage change (e.g. groundwater depletion).

Limitations

• $E$E is hard to measure directly; $PE$PE and $AE$AE diverge whenever soil moisture is limiting, so substituting one for the other introduces error.
• Deep seepage across the watershed (especially in karst, where surface and groundwater divides differ) is difficult to quantify.
• Short records may miss the long-term trend that $\Delta S$ΔS is trying to capture.

---

Located contrast: a chalk catchment versus a clay catchment

The clearest way to see how physical characteristics control basin behaviour is to contrast two real, neighbouring catchment types in southern England.