Hydrogeology: Porosity and Permeability, Aquifer Properties, Darcy’s Law, Groundwater Flow, Well Hydraulics and Exploration
2. Groundwater occurrence, storage and aquifer properties
Below the surface the unsaturated (vadose) zone holds soil water, intermediate water and the capillary fringe; below the water table (pressure equal to atmospheric) the saturated zone holds groundwater. An aquifer yields usable water; an aquitard is a poorly permeable layer that leaks; an aquiclude (clay) transmits negligible water. An unconfined aquifer has a free water table. A confined aquifer lies between aquitards, its water is under pressure and rises in a well to the piezometric surface; where that surface is above ground the well flows (artesian). A perched water table sits on a local aquitard above the main one.
- Specific yield S_y: volume of water an unconfined aquifer yields per unit area per unit fall of the water table, by draining pores; typically 0.01 to 0.3.
- Specific storage S_s: volume released per unit volume of aquifer per unit fall of head by the elastic expansion of water and compaction of the skeleton, S_s = ρg(α + nβ). Storativity S = S_s b is dimensionless: roughly 10⁻⁵ to 10⁻³ for confined aquifers, while for unconfined ones S ≈ S_y.
- Transmissivity T = K b (m²/d) is the flow per unit width under unit gradient through the full saturated thickness b, so Q = T i w across a width w.
Water released is V = S·A·Δh. For an unconfined aquifer of 50 km² with S_y = 0.12, a 2 m fall of the water table releases 50 × 10⁶ × 0.12 × 2 = 12 million m³. A confined aquifer of the same area with S = 2 × 10⁻⁴ releases only 0.2 million m³ for a 20 m head fall: it gives water by pressure relief, not by draining pores, and stays fully saturated. In the water-table fluctuation method the recharge of a season is R = S_y Δh per unit area (plus the draft, if the rise observed is net of pumping).
3. Darcy’s law and its limits
Henry Darcy (1856), working on the sand filters of Dijon, found that discharge is proportional to the area and the head loss and inversely proportional to the length: q = Q/A = −K dh/dl = K i, and Q = K i A. The hydraulic head is h = z + p/(ρg) (elevation plus pressure head; velocity head is negligible). Water moves from high to low total head, not necessarily from high to low pressure. The gradient i = Δh/L is dimensionless.
The Darcy flux q is a volume per unit area per time, not a speed of water. Water moves only through the effective pores, so the seepage (average linear) velocity is v_s = q/n_e, larger than q, and the travel time over a distance L is L/v_s. Example: heads 105.2 m and 104.7 m in wells 250 m apart give i = 0.002; with K = 20 m/d and n_e = 0.25, q = 0.04 m/d, v_s = 0.16 m/d, and 400 m takes 2500 days (about 6.8 years). Three wells are needed to fix the gradient and its direction: flow is down-gradient, perpendicular to the head contours.
- Limits. Darcy’s law holds for laminar flow, up to a Reynolds number (based on grain size) of about 1 to 10. It fails at high velocities in coarse gravel or open fractures, especially near a pumped well, and applies to saturated flow, since K in the unsaturated zone depends on moisture.
- Laboratory K. Constant-head permeameter (sands): K = VL/(A h t). Falling-head permeameter (silts, clays): K = (aL/At) ln(h₁/h₂). Small samples miss fractures and layers, so field pumping tests give the representative K.
- Flow nets and Dupuit. In a flow net, flow lines cross equipotentials at right angles and q = K H N_f/N_d per unit width (N_f flow channels, N_d head drops). For unconfined flow between two water levels h₁ and h₂ a distance L apart, Dupuit gives q = K(h₁² − h₂²)/(2L).
4. Hydraulic conductivity and groundwater flow dynamics
Refraction. At a boundary between media of conductivity K₁ and K₂, flow lines bend such that tan θ₁/tan θ₂ = K₁/K₂ (θ measured from the normal to the boundary): they turn towards the normal in the less conductive layer. Flow concentrates in the high-K layers, which is why lenses of sand and gravel, buried channels and fracture zones carry most of the water.
Governing equations. For steady flow in a homogeneous isotropic medium, head obeys Laplace’s equation, ∇²h = 0. For transient flow in a confined aquifer, ∂²h/∂x² + ∂²h/∂y² = (S/T) ∂h/∂t. The ratio T/S is the hydraulic diffusivity: it is large in confined aquifers (S is small), so a pressure change at one point is felt far away quickly, whereas in an unconfined aquifer the response is slow because S ≈ S_y is large.
Regional flow. Recharge (rain, influent streams, canals, irrigation return) enters at topographic highs and discharges at lows as springs, effluent streams, wetlands and evapotranspiration. The water table is a subdued copy of the topography, and flow systems nest as local, intermediate and regional (Tóth). In a recharge area head falls with depth (downward flow); in a discharge area head rises with depth (upward flow, artesian conditions possible). Natural velocities are small, from centimetres to a few metres per day, so residence times run from years to millennia.
5. Well hydraulics
Pumping lowers the head around a well to form a cone of depression; drawdown s = h₀ − h. For steady radial flow to a fully penetrating well (Thiem, equilibrium) in a confined aquifer, with heads h₁ and h₂ at radii r₁ and r₂, Q = 2πT(h₂ − h₁)/ln(r₂/r₁), that is s₁ − s₂ = [Q/(2πT)] ln(r₂/r₁). For an unconfined aquifer, Q = πK(h₂² − h₁²)/ln(r₂/r₁). Example: Q = 1800 m³/d, T = 600 m²/d, observation wells at 5 m and 50 m give s₁ − s₂ = (1800/(2π × 600)) × ln 10 = 1.10 m.
Theis (non-equilibrium) gives the drawdown at radius r and time t: s = Q W(u)/(4πT), u = r²S/(4Tt). It assumes a confined, homogeneous, isotropic aquifer of uniform thickness and infinite extent, a fully penetrating well of negligible diameter, an initially flat piezometric surface, constant Q and instantaneous release from storage. For u < 0.05 (later times, or small r), the Cooper-Jacob form is a straight line on a semi-log plot: s = [2.303Q/(4πT)] log₁₀(2.25Tt/(r²S)). The drawdown per log cycle Δs gives T = 2.303Q/(4πΔs) ≈ 0.183Q/Δs, and the time t₀ at which the line cuts s = 0 gives S = 2.25Tt₀/r². Example: Q = 1200 m³/d and Δs = 0.5 m give T = 440 m²/d.
Boundaries are handled by image wells. An impermeable boundary (fault, bedrock wall) acts like an image well pumping the same rate, so drawdown rises faster and the semi-log slope doubles. A recharge boundary (perennial river) acts like an image injection well, so drawdown flattens and stabilises. Well performance: the specific capacity is Q/s_w, and the total drawdown in the pumped well is s_w = BQ + CQ², where BQ is the linear aquifer loss and CQ² the well loss from turbulent flow through screen and casing. Drawdowns of neighbouring wells add (superposition). A recovery test, with residual drawdown against t/t′, gives T independently. In an unconfined aquifer the Theis equation is usable only if drawdown is small against the saturated thickness b, after correcting s′ = s − s²/(2b).
6. Groundwater exploration
Exploration proceeds from cheap to costly: a desk study of geology, geomorphology, rainfall and existing wells; remote sensing and GIS (lineaments, drainage density, geomorphic units and land use combined by weighted overlay into potential zones); field hydrogeology (well inventory, water-level contours, springs, quality); geophysics; test drilling with logging; and finally a pumping test. In alluvium and valley fills the targets are thick sand and gravel and buried channels. In karst they are solution-widened joints and conduits. In basalt they are vesicular and jointed zones, weathered tops and inter-flow beds. In crystalline hard rock (granite-gneiss) the water is held in the weathered mantle and the fractures below it, so favourable sites are thick weathering in topographic lows, lineament intersections and fracture zones, while fresh massive rock is dry and a dyke may act as a barrier.
| Method | What it senses | Groundwater use |
|---|---|---|
| Electrical resistivity (VES) | Resistivity layering with depth | Depth to aquifer and bedrock, thickness of weathering, fresh-saline contrast |
| Seismic refraction | P-wave velocity layering | Water table and bedrock in unconsolidated cover |
| Electromagnetic and VLF | Lateral conductivity changes | Fracture and shear zones, clay lenses |
| Gravity and magnetic | Density and susceptibility contrasts | Buried valleys and basins; dykes and faults |
In a vertical electrical sounding the current electrodes are spread symmetrically about a fixed centre. Apparent resistivity is ρ_a = 2πa ΔV/I for the Wenner array (spacing a) and ρ_a = π(L² − l²)/(2l) · ΔV/I for the Schlumberger array (L = AB/2, l = MN/2). Archie’s law links resistivity to the pores of a clean saturated formation: F = ρ₀/ρ_w = a φ^(−m), with a ≈ 1 and m ≈ 2 for clean sands (use the values a question states). For φ = 0.25 and water resistivity 20 Ωm, F = 16 and ρ₀ = 320 Ωm. Resistivity falls with higher porosity, saline water and clay content. Borehole logs complete the picture: resistivity and spontaneous-potential logs identify beds, and the gamma log marks clay.
Key takeaways
- Porosity n = V_v/V, with n = S_y + S_r and n_e ≤ n; K = kρg/μ. Layered ground has K_h = ΣK_i b_i/Σb_i and K_v = Σb_i/Σ(b_i/K_i). High porosity does not mean high permeability.
- Water released is V = S A Δh: unconfined S ≈ S_y (0.01–0.3) drains pores, confined S (about 10⁻⁵–10⁻³) comes from elastic storage. T = Kb, and recharge by water-table fluctuation is S_y Δh plus draft.
- Darcy: q = Ki and Q = KiA; flow follows total head h = z + p/ρg. Seepage velocity v_s = q/n_e. The law holds for laminar flow (Re below about 1–10).
- Flow refracts as tan θ₁/tan θ₂ = K₁/K₂; transient flow follows ∂²h/∂x² + ∂²h/∂y² = (S/T)∂h/∂t. Recharge areas have downward and discharge areas upward head gradients.
- Thiem (confined): s₁ − s₂ = [Q/(2πT)] ln(r₂/r₁). Theis: s = QW(u)/(4πT), u = r²S/(4Tt). Cooper-Jacob (u < 0.05): T = 2.303Q/(4πΔs), S = 2.25Tt₀/r². An impermeable boundary doubles the semi-log slope; a recharge boundary flattens it.
- Explore by geology, geomorphology, lineaments and GIS, then geophysics (VES: ρ_a = 2πaΔV/I for Wenner) and logs. Archie F = ρ₀/ρ_w = aφ^(−m); low resistivity may be clay or saline water, not a good aquifer.
Practice questions (18)
Attempt each one before opening the answer. Every explanation names the tempting wrong option as well as the right one, because that is where marks are lost.
A bed of clay and a bed of well-sorted coarse sand are compared. Which statement is correct?
Show answer
Answer: A — The clay can have the higher total porosity yet a far lower conductivity, because its pores are tiny and hold water tightly
Total porosity of clay (often 40 per cent or more) can exceed that of sand, but the pores are tiny and water is bound to the grains, so K is smaller by many orders of magnitude. Conductivity depends on pore size and connection, not on porosity alone, so the claim of higher K for clay is wrong. Fine grains hold more water against gravity, so clay, not sand, has the larger specific retention.A saturated sand has a porosity of 0.38. After free drainage under gravity it still holds water equal to 9 per cent of its total volume. Its specific yield is
Show answer
Answer: C — 0.29, the porosity minus the specific retention
Porosity splits as n = S_y + S_r. The water that stays after drainage (0.09) is the specific retention, so S_y = 0.38 − 0.09 = 0.29. The value 0.09 is the retention itself, 0.38 would mean that no water is retained at all, and adding the two quantities breaks the relation n = S_y + S_r.When the head in a confined aquifer falls, the water released from storage comes mainly from
Show answer
Answer: D — expansion of the water and compaction of the skeleton as pore pressure falls, the aquifer staying saturated
In a confined aquifer the pores stay full: a head fall lowers the pore pressure, so the water expands slightly and the skeleton compacts, which defines the small storativity S = S_s b. Gravity drainage of pores, with specific yield, belongs to unconfined aquifers. Evaporation from a capillary fringe and drainage from above do not describe release from confined storage.Two piezometers at the same site have their screens at elevations of 90 m and 70 m above a datum. Water stands at 98 m above the datum in the shallow piezometer and at 95 m above the datum in the deep one. The vertical flow between the two screens is
Show answer
Answer: A — downward, because total head is 98 m at the shallow screen and 95 m at the deep screen, so head falls with depth
Total head is the water level above the datum: 98 m at the shallow screen and 95 m at the deep one. Flow goes from high to low total head, so it is downward, a recharge setting. The pressure heads (8 m and 25 m) are a distractor: pressure is higher at depth, yet flow does not follow pressure. Elevation alone does not set the direction either, and different total heads mean there is flow.Which statement about the Darcy flux q and the seepage velocity v_s in a saturated sand is correct?
Show answer
Answer: B — v_s = q/n_e, so pore water moves faster than q suggests, as flow uses only part of the cross-section
The Darcy flux is the discharge divided by the whole cross-sectional area, but water can pass only through the effective pores, so the true average speed is v_s = q/n_e, which exceeds q because n_e < 1. Travel times and tracer speeds therefore use v_s. The flux is not the pore speed, and the relations that multiply by n_e invert the geometry of the flow.Darcy’s linear law is least likely to hold in which of the following situations?
Show answer
Answer: B — Fast flow through coarse gravel or an open fracture near a pumped well, at Reynolds number above 10
Darcy’s law assumes laminar, viscous-dominated flow, valid up to a Reynolds number of roughly 1 to 10. Near a pumped well in gravel or an open fracture, velocities are high enough for inertial and turbulent losses, so the relation between flux and gradient becomes non-linear. Slow flow in fine sand, silt or a regional sandstone aquifer is the standard domain of the law.In a pumping test on a confined aquifer, drawdown plotted against log time first follows one straight line and later a second straight line with twice the slope. This is best explained by
Show answer
Answer: B — a single impermeable boundary, equivalent to an image well pumping at the same rate
An impermeable boundary reflects the cone of depression, which is equivalent to a second well of the same sign pumping at the same rate, so after its effect arrives the drawdown rate doubles and the semi-log slope doubles. A recharge boundary or a leaky aquitard supplies water and flattens the curve, and a delayed-yield effect flattens the middle part rather than doubling the final slope.A fully penetrating well pumps 1500 m³/d from a confined aquifer of transmissivity 300 m²/d. Assume steady radial flow. The difference in drawdown between observation wells at 20 m and 200 m from the well is closest to
Show answer
Answer: C — 1.83 m, from [Q/(2πT)] ln(200/20)
Thiem’s equation gives s₁ − s₂ = [Q/(2πT)] ln(r₂/r₁) = [1500/(2π × 300)] × ln 10 = 0.7958 × 2.3026 = 1.83 m. The value 0.80 uses a base-10 logarithm, which is wrong in this form of the equation, 0.92 uses the Theis denominator 4π, and 11.5 forgets the 2π.A vertical electrical sounding in an alluvial plain shows a layer of low resistivity between 25 m and 60 m depth. Which interpretation is most defensible before drilling?
Show answer
Answer: D — It may be clay or sand saturated with saline water, so logs, lithology and water-quality data are needed to tell them apart
Both clay and saline-water-saturated sand have low resistivity, and only one of them is a useful fresh-water aquifer, so the sounding cannot decide by itself. Fresh-water sand has moderate to high resistivity, not the lowest. Dry sand is highly resistive, not conductive, and fresh unweathered crystalline rock is far more resistive than saturated sediments.In a granite-gneiss terrain, which site is the best first target for a bore well?
Show answer
Answer: D — A topographic low at a lineament intersection, over thick weathering on fractured rock
In crystalline rock the water is held in the weathered mantle and in fractures below it, so the best sites are lows with thick weathering and intersecting lineaments, which concentrate recharge and fracturing. Fresh massive rock, a hill crest or a steep scarp with thin soil offers little storage, little infiltration and few open fractures.Which of the following statements about aquifer storage properties are correct?
Show answer
Answer: A — For the same head decline an unconfined aquifer releases far more water per unit area than a confined one, because S_y is much larger than confined storativity; C — Transmissivity T = Kb has the dimensions of length squared per unit time
S_y (0.01 to 0.3) is orders of magnitude larger than confined S (about 10⁻⁵ to 10⁻³), and T = K × b has units m²/d. Specific retention is the water held against gravity, which is the opposite of what drains freely, and storativity S = S_s b is dimensionless; the quantity with units of inverse length is the specific storage S_s.Which of the following statements about groundwater flow are correct?
Show answer
Answer: B — For flow parallel to layers the effective conductivity is the thickness-weighted arithmetic mean, and for flow across the layers it is the harmonic mean; C — At a boundary between two media, flow lines refract so that tan θ₁/tan θ₂ = K₁/K₂, with the angles measured from the normal
The thickness-weighted arithmetic mean (parallel flow) and harmonic mean (flow across layers) are the standard layered-medium results, and the tangent law of refraction is tan θ₁/tan θ₂ = K₁/K₂. Water moves down total head, so a deeper point may have the higher pressure and still receive no flow from the shallower one. Seepage velocity v_s = q/n_e is larger than the Darcy flux, not smaller.Which of the following statements about well hydraulics are correct?
Show answer
Answer: A — The Theis solution assumes a confined, homogeneous, isotropic aquifer of infinite extent and a fully penetrating well of negligible diameter; D — A recharge boundary such as a perennial river acts like an image injection well and makes the drawdown flatten
The listed Theis assumptions are the standard ones, and a recharge boundary behaves as an image injection well so drawdown flattens. The Cooper-Jacob line is valid for small u (later times), not large u. In s_w = BQ + CQ², the linear term BQ is the aquifer loss and the quadratic term CQ² is the well loss, the reverse of the statement given.Two observation wells 300 m apart in a sand aquifer have heads of 52.4 m and 51.2 m. The hydraulic conductivity is 15 m/d and the effective porosity is 0.20. The time taken by groundwater to travel 600 m along the flow line at the seepage velocity is ____ days.
Numerical answer — type the value.
Show answer
Answer: 2000
Gradient i = (52.4 − 51.2)/300 = 0.004. Darcy flux q = K i = 15 × 0.004 = 0.06 m/d. Seepage velocity v_s = q/n_e = 0.06/0.20 = 0.30 m/d. Time = 600/0.30 = 2000 days. Using q instead of v_s would give 10 000 days.An unconfined aquifer of area 25 km² has specific yield 0.08. Between the pre-monsoon and post-monsoon surveys the water table rose by 2.4 m. This rise is net of a seasonal pumping draft of 1.3 million m³ taken from the aquifer in the same period; other inflows and outflows are negligible. By the water-table fluctuation method, gross recharge (net storage gain plus draft) is ____ million m³, to one decimal place.
Numerical answer — type the value.
Show answer
Answer: 6.1
Net storage gain = S_y × Δh × A = 0.08 × 2.4 × 25 × 10⁶ = 4.8 × 10⁶ m³ = 4.8 million m³. Because that rise is after the draft was withdrawn, gross recharge = 4.8 + 1.3 = 6.1 million m³. Reporting 4.8 would understate recharge by the amount pumped out.Three horizontal layers lie one above another: 4 m thick with K = 12 m/d, 2 m thick with K = 0.3 m/d and 4 m thick with K = 6 m/d. The effective vertical hydraulic conductivity of the sequence, to two decimal places, is ____ m/d.
Numerical answer — type the value.
Show answer
Answer: 1.30
For flow across layers K_v = Σb_i/Σ(b_i/K_i). Total thickness = 10 m. Σ(b_i/K_i) = 4/12 + 2/0.3 + 4/6 = 0.3333 + 6.6667 + 0.6667 = 7.6667 d. K_v = 10/7.6667 = 1.30 m/d. The arithmetic mean, 7.26 m/d, is the horizontal conductivity and is wrong for vertical flow, where the thin 0.3 m/d layer controls.Cooper-Jacob analysis of a pumping test at Q = 2400 m³/d gives a drawdown of Δs = 0.8 m per log₁₀ cycle of time. The straight line cuts zero drawdown at t₀ = 0.002 d in an observation well at r = 30 m. Use T = 2.303Q/(4πΔs) and S = 2.25Tt₀/r². The storativity S, expressed as a multiple of 10⁻³ and to two decimal places, is ____ × 10⁻³.
Numerical answer — type the value.
Show answer
Answer: 2.75
T = 2.303 × 2400/(4π × 0.8) = 5527.2/10.053 = 549.8 m²/d. S = 2.25 × 549.8 × 0.002/30² = 2.474/900 = 2.749 × 10⁻³, i.e. 2.75 × 10⁻³. The slope gives T and only then does the intercept give S; using r² = 30 instead of 900 would give a value 30 times too large.A clean saturated sand has porosity 0.20 and pore-water resistivity 15 Ωm. Use Archie’s law ρ₀ = ρ_w a φ^(−m) with a = 1 and m = 2. The bulk resistivity of the sand is ____ Ωm.
Numerical answer — type the value.
Show answer
Answer: 375
Formation factor F = a φ^(−m) = 0.20^(−2) = 25. Bulk resistivity ρ₀ = F × ρ_w = 25 × 15 = 375 Ωm. Using φ^m instead of φ^(−m) would give 0.6 Ωm, which is far below the pore-water resistivity and therefore impossible for a rock that also contains insulating grains.