Groundwater, Soils, and Sediments · Study · FE Environmental · FE → PE Prep
Groundwater, Soils, and Sediments
8% of exam
Aquifer properties and hydrogeology, Darcy's law and seepage velocity, well drawdown (Theis, Jacob, Thiem, Dupuit), and soil, sediment, and groundwater remediation.
4 concepts
A. Basic hydrogeology
Hydrogeology and Aquifer Properties
Confined, unconfined, and perched aquifers, aquitards, the porosity/specific-yield/specific-retention split, the water table vs the potentiometric surface, and recharge and discharge.
Unlocks with an access pass — one-time payment, no auto-renew. View passes
Pass holders
B. Groundwater flow
Darcy’s Law and Seepage Velocity
Darcy’s law, hydraulic conductivity, transmissivity and storativity, and the crucial gap between the Darcy flux and the seepage velocity a contaminant actually travels at.
Almost every groundwater question on the FE Environmental exam stands on a single 1856 experiment: water moving through saturated soil discharges in direct proportion to the hydraulic gradient driving it. Get Darcy's law right and you can compute well yields, plume arrival times, capture zones, and seepage under dams; get the velocity definitions backwards and you will be off by a factor of three to five on every travel-time problem. The single most common — and most expensive — error in this whole topic is reporting the Darcy flux q when the question asks how fast the water (or a contaminant) actually moves. This concept builds the fluency to keep the discharge, the specific flux, and the seepage velocity cleanly separated. The governing relations are in the FE Reference Handbook — Civil Engineering under 'Darcy's Law'.
Darcy’s law and hydraulic conductivity
Darcy found that volumetric discharge Q through a porous medium is proportional to the cross-sectional area A and to the hydraulic gradient −dh/dx, with the proportionality constant being the hydraulic conductivity K. The minus sign encodes the physics: groundwater flows from high head to low head, so positive discharge accompanies a falling head. K (units of velocity, ft/s or m/s) bundles together both the medium — grain size, sorting, pore connectivity — and the fluid's density and viscosity. Clean gravel may carry K≈102m/d while a clay aquitard sits near 10−4m/d, a span of six orders of magnitude that explains why one layer is an aquifer and the one beneath it confines.
Q=−KAdxdh
Specific discharge — the Darcy flux
Divide discharge by the full cross-sectional area and you get the specific discharge q, also called the Darcy velocity or superficial velocity. It has units of velocity but it is not a real velocity — it is the discharge that would occur if water moved through the entire face of the medium, solids included. q is exactly what you need for flux and mass-loading calculations (volume of water per unit total area per time), and it is what Darcy's law delivers directly. Treating q as the speed of a water particle is the canonical mistake, because water cannot flow through the grains, only around them.
q=AQ=−Kdxdh
Seepage velocity — how fast water really moves
Water threads only through the connected pore space, so its actual average velocity is larger than q. Dividing the Darcy flux by the effective porosity ne — the fraction of the bulk volume made of interconnected, flow-carrying pores — gives the average seepage velocity (also called the linear or pore velocity). Because ne is typically 0.1
Effective vs total porosity
Total porosity n counts all void space; effective porosity ne counts only the interconnected voids that actually conduct flow. Dead-end pores and water bound to grain surfaces hold water but do not transmit it, so ne≤n always — and the gap can be large in fine-grained or poorly connected media. Clay can have a total porosity near
Transmissivity and storativity
When you treat a whole aquifer rather than a point, two integrated properties take over. Transmissivity T=Kb is the conductivity times the saturated thickness b (units L2/T, e.g. m2
Gradient and flow direction
The hydraulic gradient is the change in head per unit distance, i=−dh/dx, and flow follows the steepest descent of head, perpendicular to the equipotential (equal-head) contours. With three monitoring wells you can triangulate both the magnitude and the compass direction of the gradient: contour the heads, then drop a perpendicular from high to low. Field problems often give you head differences across a known distance — keep the gradient dimensionless by carrying head and distance in the same length units, and remember the gradient points downhill in head, which is the direction the plume will go.
Exam strategy
Read the verb: 'flux', 'specific discharge', or 'Darcy velocity' means q=Ki; 'how fast', 'travel time', 'pore' or 'seepage velocity' means v=q/ne. Write down which one is asked before you compute. Build the gradient as a clean dimensionless ratio (equal length units top and bottom), and watch K
Key equations
Darcy’s law (discharge)Q=−KAdxdh
Q = volumetric discharge (m³/s, ft³/s); K = hydraulic conductivity (m/s, ft/s); A = total cross-sectional area; dh/dx = hydraulic gradient. Minus sign: flow goes down-gradient.
Specific discharge (Darcy flux)q=AQ=−Kdxdh=Ki
Seepage (pore) velocityv=neq=neKi
Hydraulic gradienti=−dxdh=Lh1−h2
TransmissivityT=Kb
Conductivity times saturated thickness b; units L2/T (m²/d, ft²/s). Aquifer-scale measure of water-transmitting capacity.
Aquifer discharge through a widthQ=Tiw=Kbiw
Discharge through a vertical strip of aquifer of width w and thickness b; folds thickness into T
Problem. A confined sand aquifer has hydraulic conductivity K=15m/d and effective porosity ne=0.25. Two wells 500m apart along the flow path show a head difference of 2.0m. Find (a) the Darcy flux, (b) the seepage velocity, and (c) the time for a conservative tracer to travel between the wells.
Problem. An unconfined aquifer is 12m thick with K=25m/d. Estimate its transmissivity, then find the groundwater underflow passing through a 200m-wide section where the gradient is 0.003.
Solution.
Common pitfalls
•Reporting the Darcy flux q when the question asks for velocity or travel time — divide by ne to get the seepage velocity v. This factor of 3–10 is the single most common error in the topic.
•Using total porosity n instead of effective porosity ne in v=q/ne; bound water and dead-end pores do not transmit flow, so ne≤n.
•Forgetting the time-unit mismatch in K: handbook conductivities are often ft/s or m/s, but problems want ft/day or m/day — convert before computing (×86,400 s/day).
•Confusing confined and unconfined storativity: confined S is tiny (∼10−4, elastic storage), while unconfined S≈Sy∼0.2
•Treating transmissivity as a velocity or a conductivity; T=Kb has units L2/T and already includes the saturated thickness, so do not multiply by b again.
•Dropping the sign/direction of the gradient: flow goes toward LOW head. Building i with head and distance in different length units (ft head over m distance) silently corrupts every result.
Davis & Cornwell, Introduction to Environmental Engineering — groundwater flow fundamentals
C. Drawdown
Well Drawdown: Theis, Jacob, Thiem, and Dupuit
Steady drawdown in confined (Thiem) and unconfined (Dupuit) aquifers, the cone of depression and radius of influence, transient drawdown (Theis/Cooper-Jacob), and specific capacity.
Unlocks with an access pass — one-time payment, no auto-renew. View passes
Pass holders
D. Remediation of soil, sediment, and/or groundwater
Contaminant Transport and Remediation
Advection-dispersion, the retardation factor that slows a sorbing plume, partitioning from foc and Koc, and the matching of remediation technologies to contaminant behavior.
Unlocks with an access pass — one-time payment, no auto-renew. View passes
Pass holders
to
0.35
, the seepage velocity is roughly three to ten times the Darcy flux. This is the velocity that governs how fast a conservative tracer or the leading edge of a plume travels, and therefore the velocity behind every contaminant travel-time problem.
v=neq=neKdxdh
0.5
yet an effective porosity an order of magnitude smaller. For seepage velocity you must use
ne
; mixing in the total porosity inflates the denominator and underestimates how fast water moves.
/d); it measures how much water the full thickness of the aquifer can transmit and is the natural parameter in well equations. Storativity (storage coefficient)
S
is the volume of water released from or taken into storage per unit aquifer area per unit drop in head — dimensionless. In a confined aquifer the water comes from elastic expansion of water and compression of the skeleton, so
S
is tiny, roughly
10−5
to
10−3
. In an unconfined aquifer storage is mostly gravity drainage of pores, so
S
approaches the specific yield, typically
0.05
to
0.30
— hundreds of times larger.
T=KbS=Ssb(confined),S≈Sy(unconfined)
's time unit — handbook
K
is often ft/s while the problem wants ft/day, a factor of
86,400
. For travel time across a distance
L
, use
t=L/v=Lne/(Ki)
, never
L/q
. For aquifer-scale discharge through a strip of width
w
, combine with transmissivity as
Q=Tiw
, which already folds in the thickness.
Discharge per unit total area; units of velocity but NOT a particle speed. Use for flux/mass-loading. i=−dh/dx is the gradient (dimensionless).
Actual average water velocity through the interconnected pores; ne = effective porosity. This is the velocity for travel-time and plume problems.
Head drop per unit distance along the flow path, dimensionless. Use consistent length units for head and distance.
.
Volume released per unit area per unit head drop, dimensionless. Confined S∼10−5–10−3; unconfined S≈Sy∼0.05–0.30.
Time for water (or a non-sorbing tracer) to travel distance L at the seepage velocity. Uses ne, not q.
.
(b) Seepage velocity:
v=neq=0.250.060=0.240m/d
.
(c) Travel time:
t=vL=0.240500=2.08×103d≈5.70yr
.
Sanity check:
v
is
4×
the flux because
1/ne=1/0.25=4
; using
q
instead of
v
would have wrongly given
8,330
days. Units:
(m/d)/(dimensionless)=m/d
, and
m÷(m/d)=d
.
Transmissivity: T=Kb=25×12=300m2/d.
Underflow through the strip: Q=Tiw=300×0.003×200=180m3/d.
Final: Q=1.80×102m3/d.
Sanity check via Darcy's law directly: A=bw=12×200=2400m2, q=Ki=25×0.003=0.075m/d, so Q=qA=0.075×2400=180m3/d. Agreement confirms T=Kb simply repackages thickness. Units: (m2/d)()(m)=m3/d.