Question 1 of 5: Confined-Aquifer Storage, Three-Point Hydraulic Gradient, and Density-Driven Flow Across an Aquitard
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Reference texts: Freeze & Cherry, Groundwater (Prentice-Hall, 1979) — Darcy's law, storativity/specific storage, three-point hydraulic-gradient estimation, density-dependent (freshwater-equivalent-head) flow across an aquitard, the Theis and Thiem well equations, Hantush-Jacob leaky-aquifer theory, superposition for wells with variable pumping schedules, and the Dupuit-Forchheimer approximation with areal recharge; Kruseman & de Ridder, Analysis and Evaluation of Pumping Test Data (ILRI, 1994) — Cooper-Jacob straight-line method and the Hvorslev slug-test method; EGBC Geoscience Professional Practice Guidelines for assumption-disclosure conventions on open-book calculations.
Question 1: Confined-Aquifer Storage, Three-Point Hydraulic Gradient, and Density-Driven Flow Across an Aquitard (20 marks)
$1.1\times10^{-6}\ \text{m}^{-1}$, 55 m, 1500 m, 2400 m, 2 m
(b)
Head at A, B (250 m N of A), C (250 m E of A)
180, 196, 188 m.a.s.l.
(c)
$\rho_f$, $\rho_{s,\text{true}}$, aquitard thickness, $\psi_1$ (fresh well column), $\psi_2$ (saline well column)
998, 1150 kg/m³, 22 m, 10 m, 30 m
Find. (a) The volume of water released by a 2 m piezometric decline. (b) The magnitude and direction of the hydraulic gradient from the three water levels. (c) The direction of flow across the aquitard under the two stated density assumptions.
Approach. Part (a) converts specific storativity to a dimensionless storativity ($S=S_sb$) and multiplies by plan area and head decline. Part (b) sets up a local $x$–$y$ (east–north) coordinate system at Well A and solves the two independent one-dimensional gradients (A→C for the east component, A→B for the north component) as a vector. Part (c) converts each well's pressure-head reading to a common freshwater-equivalent total head using $h_f=z+(\rho/\rho_f)\psi$, once with the saline water mis-assumed fresh (case i) and once with its true density (case ii), and compares the two heads to find the flow direction in each case.
Part (a) — storativity and released volume. $S=S_sb=(1.1\times10^{-6})(55)=6.05\times10^{-5}$ (dimensionless). Over the plan area $A=1500\times2400=3.6\times10^6\ \text{m}^2$:
$$V=S\cdot A\cdot\Delta h=(6.05\times10^{-5})(3.6\times10^6)(2)=\boxed{435.6\ \text{m}^3}.$$
Part (b) — gradient components. Placing Well A at the origin, Well B is at $(0,250)$ (north) and Well C at $(250,0)$ (east). Each is a one-dimensional gradient along its own axis:
$$\frac{\partial h}{\partial y}\bigg|_{\text{north}}=\frac{196-180}{250}=0.0640,\qquad \frac{\partial h}{\partial x}\bigg|_{\text{east}}=\frac{188-180}{250}=0.0320.$$
Magnitude and direction. The gradient vector's magnitude and bearing (measured clockwise from north) are:
$$|\nabla h|=\sqrt{0.0640^2+0.0320^2}=\boxed{0.0716\ (7.16\%)},\qquad \theta=\tan^{-1}\!\left(\frac{0.0320}{0.0640}\right)=\boxed{N26.6^{\circ}E}\ \text{(direction of rising head)}.$$
Groundwater flows down-gradient, i.e. exactly opposite this bearing: $\boxed{S26.6^{\circ}W}$, from the Well-B corner of the triangle toward the Well-A/C side.
Part (c) — freshwater-equivalent heads. Take the top of the aquitard as datum ($z_1=0$) so the bottom sits at $z_2=-22\ \text{m}$. The fresh-water well head is unaffected by the density conversion (its own fluid is already the reference fluid):
$$h_{f,1}=z_1+\psi_1=0+10=\boxed{10.00\ \text{m}}.$$
For the saline well, the freshwater-equivalent head is $h_f=z+(\rho/\rho_f)\psi$, applied twice with the two stated densities:
$$\text{(i) }\rho=998:\ h_{f,2}=-22+\left(\frac{998}{998}\right)(30)=\boxed{8.00\ \text{m}}$$
$$\text{(ii) }\rho=1150:\ h_{f,2}=-22+\left(\frac{1150}{998}\right)(30)=\boxed{12.57\ \text{m}}.$$
Direction in each case. Case (i): $h_{f,1}=10.00\ \text{m} > h_{f,2}=8.00\ \text{m}$, so flow is downward, from the fresh aquifer into the saline aquifer. Case (ii): $h_{f,2}=12.57\ \text{m} > h_{f,1}=10.00\ \text{m}$ — the ranking flips, and flow is upward, from the saline aquifer into the fresh aquifer, even though the denser fluid sits below.
Figure 1 — Plan view of the three-well array (Q1b). The head-increase (gradient) direction points N26.6°E toward the Well-B corner; groundwater actually flows the opposite way, down-gradient, S26.6°W.
Figure 2 — The fresh/saline aquifer system of Q1c. With the true saline density (case ii), the freshwater-equivalent head at the base of the aquitard (12.57 m) exceeds the fresh-aquifer head (10.00 m), so flow is upward across the aquitard despite the denser fluid sitting below.
Quantity
Result
(a) Storativity $S$
6.05×10⁻⁵
(a) Volume of water pumped
435.6 m³
(b) Gradient magnitude
0.0716 (7.16%)
(b) Gradient direction (rising head) / flow direction