Soil Water, Darcy Permeability and Seepage Flow Nets

Outcomes

After completing this unit the learner should be able to define the governing state variables, calculate representative quantities with correct units, explain when the idealized relation breaks down, interpret the corresponding laboratory models, and defend a field investigation or verification decision.

Core lessons

Lesson 2.1 — Hydraulic head and Darcy flow

Central question. Groundwater movement follows hydraulic head differences; seepage direction depends on total head, not simply depth.

Governing relationships and definitions. Hydraulic head h=z+u/γw; hydraulic gradient i=Δh/L; for saturated laminar 1D flow Q=k i A. In stratified soil, vertical and horizontal equivalent conductivities differ.

Worked calculation / observation. A 0.35-m head loss over 7 m gives i=0.05. With k=2×10⁻⁵ m/s and A=3 m², Q=3×10⁻⁶ m³/s. If the head difference doubles, ideal Darcy discharge doubles.

Why interpretation matters. The sign of hydraulic gradient is a bookkeeping choice: choose a positive axis and state whether head declines along it. The Darcy proportionality should be applied only after this convention is fixed. Students often compute a correct magnitude but assign flow in the wrong direction after comparing depths instead of total hydraulic heads.

Field case: A dewatering wellpoint system near an adjacent basement. In a dewatering wellpoint system near an adjacent basement, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a dewatering wellpoint system near an adjacent basement that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a dewatering wellpoint system near an adjacent basement. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. The law is not valid for every high-gradient or very coarse flow regime. Also avoid replacing field hydraulic conductivity with a laboratory value without scale assessment. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Lesson 2.2 — Layered permeability and seepage force

Central question. Layering creates anisotropy at the engineering scale even when each layer is treated as isotropic.

Governing relationships and definitions. For horizontal flow parallel to layers, kh=(Σki Hi)/(ΣHi); for vertical flow, kv=(ΣHi)/(ΣHi/ki). Seepage force per unit volume has magnitude iγw.

Worked calculation / observation. For two equal layers, k1=10⁻⁴ and k2=10⁻⁶ m/s, kh=5.05×10⁻⁵ m/s but kv≈1.98×10⁻⁶ m/s. This order-of-magnitude difference changes dewatering predictions.

Why interpretation matters. Hydraulic conductivity is a soil–fluid system parameter. Temperature changes fluid viscosity, and structure changes the connected pore pathways; the parameter should never be described as a constant intrinsic to grain type alone. In field work, a single laboratory specimen may miss higher-conductivity seams that dominate the whole seepage problem.

Field case: A temporary cofferdam with upward seepage at the toe. In a temporary cofferdam with upward seepage at the toe, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a temporary cofferdam with upward seepage at the toe that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a temporary cofferdam with upward seepage at the toe. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. Harmonic averaging assumes steady one-dimensional flow and layers spanning the analyzed cross-section. It does not represent preferential cracks. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Lesson 2.3 — Total and effective vertical stress

Central question. Soil skeleton stress, not total overburden alone, controls consolidation and frictional resistance.

Governing relationships and definitions. At depth z with a hydrostatic water table at zw, total stress σv is the integral of bulk unit weight; pore water pressure u=γw max(z−zw,0); effective vertical stress σ′v=σv−u.

Worked calculation / observation. Take 2 m moist soil at 18 kN/m³ above water and 3 m saturated soil at 20 kN/m³ below: σv=96 kPa; u=29.43 kPa; σ′v=66.57 kPa at 5 m.

Why interpretation matters. Effective stress is a force-balance concept: the mineral skeleton carries the difference between total normal traction and pore-water pressure under appropriate assumptions. It is not a separate field instrument reading. Designers infer it by combining overburden, groundwater observations and possible excess pore pressures induced by construction or weather.

Field case: A bridge abutment with monitored groundwater levels. In a bridge abutment with monitored groundwater levels, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a bridge abutment with monitored groundwater levels that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a bridge abutment with monitored groundwater levels. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. A perched aquifer, artesian head or upward seepage breaks the hydrostatic assumption. A field piezometer measures water head, not simply the water table inferred from color. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Lesson 2.4 — Capillarity and hydraulic gradients

Central question. Water may rise above a free water surface because of surface tension and pore geometry, changing suction and near-surface state.

Governing relationships and definitions. The capillary tube relation hc=2T cosθ/(ρw g r) is an idealized capillary analogy. Soil pore-size distributions require a retention curve; suction is not universally equivalent to a single capillary tube.

Worked calculation / observation. For water surface tension 0.072 N/m, contact angle 0°, radius 0.10 mm, hc≈0.147 m; reducing radius by ten increases the tube-predicted rise tenfold.

Why interpretation matters. An upward gradient can reduce effective stress, potentially toward a condition of boiling or piping under special circumstances. The ideal critical gradient i_cr=(Gs−1)/(1+e) is a useful vertical, cohesionless, saturated-soil benchmark; it is not a comprehensive filter or erosion design criterion.

Field case: A clay–sand sequence with differing seepage routes. In a clay–sand sequence with differing seepage routes, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a clay–sand sequence with differing seepage routes that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a clay–sand sequence with differing seepage routes. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. Do not use a single pore radius to predict unsaturated hydraulic conductivity or real ground suction. Capillary rise may affect sampling, subgrades and salt transport. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Laboratory activities

Laboratory 3: Layered Darcy Seepage. Head loss and equivalent conductivities of horizontal layers. First predict the direction of change for the output when one control changes. Hold all other controls fixed; record a baseline, two variations, and a short discussion of why the governing assumptions are or are not credible. Navigate to ../labs/lab-03.html. Important: Steady, saturated, laminar flow in two continuous horizontal strata.

Laboratory 4: Effective Stress with Water Table. Piecewise total, pore and effective vertical stresses. First predict the direction of change for the output when one control changes. Hold all other controls fixed; record a baseline, two variations, and a short discussion of why the governing assumptions are or are not credible. Navigate to ../labs/lab-04.html. Important: Hydrostatic pressure with a single horizontal groundwater surface.

Integrated design conversation

Imagine that the site data are incomplete and construction will begin in phases. Construct two defensible interpretations of the subsurface conditions. For each, state the best estimate of the controlling material parameter and a credible low and high value. Explain how a field observation or laboratory test could reject one interpretation. The deliverable is not only a number but a documented chain from field evidence to calculation, checking, and the resulting construction choice.

Applied dossier: seepage at an excavation toe

A temporary excavation intersects a permeable sand layer above a lower-permeability silty layer. The groundwater head outside the excavation is 2.8 m above the interior control level, and the approximate vertical seepage path through a silty base is 6 m. A first screening gradient is i=2.8/6≈0.467. This ratio is dimensionless. If a representative vertical hydraulic conductivity is 3×10⁻⁷ m/s, the Darcy velocity magnitude is about 1.4×10⁻⁷ m/s, equivalent to 12.1 mm/day. Darcy velocity is discharge divided by gross cross-sectional area, not the mean pore-water velocity. The latter would require division by the flow-effective porosity and a more complete transport model.

Suppose mineral specific gravity is 2.66 and void ratio 0.72. The elementary saturated cohesionless critical-gradient estimate is (Gs−1)/(1+e)≈0.965. A gradient of 0.467 is below this benchmark, but that fact does not certify the excavation safe from piping, heave or erosion. The gradient may be concentrated around sheet-pile tips; the actual flow path is not straight and uniform, and soil erosion can initiate below the ideal global critical gradient. Real assessments must examine filters, exit gradients, geometry, stratification and staged pumping.

The stress calculation at the base provides a second perspective. If a soil slice carries 100 kPa total vertical overburden and an unusually high upward pressure of 75 kPa at the critical point, the effective vertical stress is only 25 kPa. The skeleton therefore has much less normal confinement than a total-stress chart suggests. It is the spatial distribution and evolution of this pore pressure that matters, not a single water-table mark placed on a geological section.

An investigation proposal should include two piezometers screened at different depths and a record of pumping response. If the deeper instrument reads a higher hydraulic head, upward seepage is plausible; if equal heads are observed, the assumed pressure field must be revised. Report a range of conductivities and investigate coarse seams because an order-of-magnitude variation in a thin layer can dominate the flow. The model supports a defensible measurement plan, not a finished dewatering design.

Technical methods notebook: groundwater energy and layered media

Head as energy per unit weight. Groundwater flow is driven by spatial variation in hydraulic head. In a saturated pore-water domain the piezometric head is the sum of elevation head z and pressure head u/γw. If a piezometer screen is 3 m below a chosen datum but its water column stands 5 m above the same datum, the measured total head is 5 m; it is not −3 m. Explicit datums prevent false claims about upflow and downflow. Compute the difference between two screened heads and divide by path length to obtain a hydraulic gradient. Flow direction follows decreasing head, although the coordinate sign must be defined.

Darcy discharge and velocity. The magnitude Q=k i A uses the hydraulic conductivity in m/s, dimensionless gradient and total area normal to flow in m². The result is m³/s. If a gross cross-section is 2 m², k=10⁻⁵ m/s and i=0.2, then Q=4×10⁻⁶ m³/s. A travel-time discussion needs the effective pore velocity, which is approximately specific discharge divided by effective porosity under an appropriate flow model. The two velocities should not be labeled interchangeably. Neither relation predicts a contaminant breakthrough curve without dispersion, sorption and boundary conditions.

Layer equivalence. Parallel layers carry flow under a common head gradient; their discharges add. The relevant equivalent horizontal conductivity is thickness-weighted arithmetic mean, kh=ΣkiHi/ΣHi. Perpendicular layers carry the same steady discharge but divide the head loss; the effective vertical conductivity is the thickness-weighted harmonic mean, kv=ΣHi/Σ(Hi/ki). Take a 1 m high-conductivity layer at 10⁻⁴ m/s and a 3 m silty layer at 10⁻⁶ m/s. Then kh=(0.0001×1+0.000001×3)/4=2.575×10⁻⁵ m/s. The vertical equivalent is 4/(1/0.0001+3/0.000001)≈1.329×10⁻⁶ m/s. Treating both paths with the arithmetic mean would drastically overestimate vertical flow.

Stress consequences. Integrate total overburden from unit-weight measurements and subtract the measured or modeled pore pressure. Under simple static groundwater, u=γw(z−zw) below the water level, and zero gauge pressure is often assigned above it. Real partially saturated soils may have suction and perched water. When water heads are not hydrostatic, pressure at a given point must be obtained from its piezometric head rather than from the distance to a visual water line. That correction may alter predicted excavation stability or settlement.

Design boundary. A critical hydraulic gradient is not a comprehensive erosion threshold. Backward erosion piping can evolve along preferential flow paths and filter gradations matter. A simplistic classroom safety ratio using an average gradient should never be used to certify an embankment, excavated base or dam exit. A site-specific seepage geometry, monitoring plan and appropriate erosion assessment are required.

Measurement task. Sketch two piezometers at different elevations and invent physically consistent head readings for upward, downward and zero vertical flow. For each sketch, calculate the average vertical gradient, specify the sign and show what changes if the screen locations are swapped. Explain which observation—the screen depth or the water level in the instrument—determines hydraulic head.

Unit 2: Soil water and seepage

Permeability, seepage, effective stress and capillary effects.

Outcomes

After completing this unit the learner should be able to define the governing state variables, calculate representative quantities with correct units, explain when the idealized relation breaks down, interpret the corresponding laboratory models, and defend a field investigation or verification decision.

Core lessons

Lesson 2.1 — Hydraulic head and Darcy flow

Central question. Groundwater movement follows hydraulic head differences; seepage direction depends on total head, not simply depth.

Governing relationships and definitions. Hydraulic head h=z+u/γw; hydraulic gradient i=Δh/L; for saturated laminar 1D flow Q=k i A. In stratified soil, vertical and horizontal equivalent conductivities differ.

Worked calculation / observation. A 0.35-m head loss over 7 m gives i=0.05. With k=2×10⁻⁵ m/s and A=3 m², Q=3×10⁻⁶ m³/s. If the head difference doubles, ideal Darcy discharge doubles.

Why interpretation matters. The sign of hydraulic gradient is a bookkeeping choice: choose a positive axis and state whether head declines along it. The Darcy proportionality should be applied only after this convention is fixed. Students often compute a correct magnitude but assign flow in the wrong direction after comparing depths instead of total hydraulic heads.

Field case: A dewatering wellpoint system near an adjacent basement. In a dewatering wellpoint system near an adjacent basement, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a dewatering wellpoint system near an adjacent basement that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a dewatering wellpoint system near an adjacent basement. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. The law is not valid for every high-gradient or very coarse flow regime. Also avoid replacing field hydraulic conductivity with a laboratory value without scale assessment. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Lesson 2.2 — Layered permeability and seepage force

Central question. Layering creates anisotropy at the engineering scale even when each layer is treated as isotropic.

Governing relationships and definitions. For horizontal flow parallel to layers, kh=(Σki Hi)/(ΣHi); for vertical flow, kv=(ΣHi)/(ΣHi/ki). Seepage force per unit volume has magnitude iγw.

Worked calculation / observation. For two equal layers, k1=10⁻⁴ and k2=10⁻⁶ m/s, kh=5.05×10⁻⁵ m/s but kv≈1.98×10⁻⁶ m/s. This order-of-magnitude difference changes dewatering predictions.

Why interpretation matters. Hydraulic conductivity is a soil–fluid system parameter. Temperature changes fluid viscosity, and structure changes the connected pore pathways; the parameter should never be described as a constant intrinsic to grain type alone. In field work, a single laboratory specimen may miss higher-conductivity seams that dominate the whole seepage problem.

Field case: A temporary cofferdam with upward seepage at the toe. In a temporary cofferdam with upward seepage at the toe, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a temporary cofferdam with upward seepage at the toe that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a temporary cofferdam with upward seepage at the toe. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. Harmonic averaging assumes steady one-dimensional flow and layers spanning the analyzed cross-section. It does not represent preferential cracks. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Lesson 2.3 — Total and effective vertical stress

Central question. Soil skeleton stress, not total overburden alone, controls consolidation and frictional resistance.

Governing relationships and definitions. At depth z with a hydrostatic water table at zw, total stress σv is the integral of bulk unit weight; pore water pressure u=γw max(z−zw,0); effective vertical stress σ′v=σv−u.

Worked calculation / observation. Take 2 m moist soil at 18 kN/m³ above water and 3 m saturated soil at 20 kN/m³ below: σv=96 kPa; u=29.43 kPa; σ′v=66.57 kPa at 5 m.

Why interpretation matters. Effective stress is a force-balance concept: the mineral skeleton carries the difference between total normal traction and pore-water pressure under appropriate assumptions. It is not a separate field instrument reading. Designers infer it by combining overburden, groundwater observations and possible excess pore pressures induced by construction or weather.

Field case: A bridge abutment with monitored groundwater levels. In a bridge abutment with monitored groundwater levels, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a bridge abutment with monitored groundwater levels that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a bridge abutment with monitored groundwater levels. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. A perched aquifer, artesian head or upward seepage breaks the hydrostatic assumption. A field piezometer measures water head, not simply the water table inferred from color. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Lesson 2.4 — Capillarity and hydraulic gradients

Central question. Water may rise above a free water surface because of surface tension and pore geometry, changing suction and near-surface state.

Governing relationships and definitions. The capillary tube relation hc=2T cosθ/(ρw g r) is an idealized capillary analogy. Soil pore-size distributions require a retention curve; suction is not universally equivalent to a single capillary tube.

Worked calculation / observation. For water surface tension 0.072 N/m, contact angle 0°, radius 0.10 mm, hc≈0.147 m; reducing radius by ten increases the tube-predicted rise tenfold.

Why interpretation matters. An upward gradient can reduce effective stress, potentially toward a condition of boiling or piping under special circumstances. The ideal critical gradient i_cr=(Gs−1)/(1+e) is a useful vertical, cohesionless, saturated-soil benchmark; it is not a comprehensive filter or erosion design criterion.

Field case: A clay–sand sequence with differing seepage routes. In a clay–sand sequence with differing seepage routes, the engineer should first sketch the geometry, state the loading or groundwater history, and identify which of the measured quantities are direct observations. Next identify the controlling stress, volume, hydraulic or strength variable from this lesson. A prediction is only as defensible as the link between the laboratory variable and its field counterpart. For example, the stress relevant to settlement may be an effective stress after excess water pressure dissipates, whereas the load relevant to a temporary wall brace may be a short-term total-stress action. Mixing those questions can produce a plausible-looking but physically inconsistent answer.

From observation to a defensible model. A practical investigation has three stages. Observation: identify the features of a clay–sand sequence with differing seepage routes that a field crew can actually record, such as layer thickness, water-table history, moisture state or deformation. Interpretation: select the idealized law described above, explain the chosen signs and units, and identify which soil parameters need independent tests. Prediction: evaluate a baseline and at least two credible alternatives. A single baseline is rarely an adequate basis for choosing a foundation treatment, drainage provision or construction method.

Verification practice. A useful checking sequence is dimensional analysis, limiting behavior and independent arithmetic. Insert the stated units before using a formula; a conductivity given in cm/s cannot be entered as m/s unchanged. Ask whether a computed stress or discharge tends toward zero when its driving load is removed. Then compare with a second hand calculation, a plotted numerical response or a bounding case. A graph from an interactive model provides feedback, not evidence that the assumed governing equation is correct for the site.

Reasoning exercise. Draw the physical control volume or stress path for a clay–sand sequence with differing seepage routes. Label the measured variables, the inferred variables and the quantities that could change after construction. Which variable would be the cheapest to measure more accurately? Describe how two alternative assumptions could change the outcome in opposite directions. Finally explain which decision would be postponed until a more representative specimen, groundwater observation or construction-stage measurement becomes available.

Model boundary. Do not use a single pore radius to predict unsaturated hydraulic conductivity or real ground suction. Capillary rise may affect sampling, subgrades and salt transport. Present the result with an explicit number of significant digits justified by the inputs; do not confuse that numerical precision with geological certainty. Professional design requires the governing standard, site-specific investigation and independent review.

Laboratory activities

Laboratory 3: Layered Darcy Seepage. Head loss and equivalent conductivities of horizontal layers. First predict the direction of change for the output when one control changes. Hold all other controls fixed; record a baseline, two variations, and a short discussion of why the governing assumptions are or are not credible. Navigate to ../labs/lab-03.html. Important: Steady, saturated, laminar flow in two continuous horizontal strata.

Laboratory 4: Effective Stress with Water Table. Piecewise total, pore and effective vertical stresses. First predict the direction of change for the output when one control changes. Hold all other controls fixed; record a baseline, two variations, and a short discussion of why the governing assumptions are or are not credible. Navigate to ../labs/lab-04.html. Important: Hydrostatic pressure with a single horizontal groundwater surface.

Integrated design conversation

Imagine that the site data are incomplete and construction will begin in phases. Construct two defensible interpretations of the subsurface conditions. For each, state the best estimate of the controlling material parameter and a credible low and high value. Explain how a field observation or laboratory test could reject one interpretation. The deliverable is not only a number but a documented chain from field evidence to calculation, checking, and the resulting construction choice.

Applied dossier: seepage at an excavation toe

A temporary excavation intersects a permeable sand layer above a lower-permeability silty layer. The groundwater head outside the excavation is 2.8 m above the interior control level, and the approximate vertical seepage path through a silty base is 6 m. A first screening gradient is i=2.8/6≈0.467. This ratio is dimensionless. If a representative vertical hydraulic conductivity is 3×10⁻⁷ m/s, the Darcy velocity magnitude is about 1.4×10⁻⁷ m/s, equivalent to 12.1 mm/day. Darcy velocity is discharge divided by gross cross-sectional area, not the mean pore-water velocity. The latter would require division by the flow-effective porosity and a more complete transport model.

Suppose mineral specific gravity is 2.66 and void ratio 0.72. The elementary saturated cohesionless critical-gradient estimate is (Gs−1)/(1+e)≈0.965. A gradient of 0.467 is below this benchmark, but that fact does not certify the excavation safe from piping, heave or erosion. The gradient may be concentrated around sheet-pile tips; the actual flow path is not straight and uniform, and soil erosion can initiate below the ideal global critical gradient. Real assessments must examine filters, exit gradients, geometry, stratification and staged pumping.

The stress calculation at the base provides a second perspective. If a soil slice carries 100 kPa total vertical overburden and an unusually high upward pressure of 75 kPa at the critical point, the effective vertical stress is only 25 kPa. The skeleton therefore has much less normal confinement than a total-stress chart suggests. It is the spatial distribution and evolution of this pore pressure that matters, not a single water-table mark placed on a geological section.

An investigation proposal should include two piezometers screened at different depths and a record of pumping response. If the deeper instrument reads a higher hydraulic head, upward seepage is plausible; if equal heads are observed, the assumed pressure field must be revised. Report a range of conductivities and investigate coarse seams because an order-of-magnitude variation in a thin layer can dominate the flow. The model supports a defensible measurement plan, not a finished dewatering design.

Technical methods notebook: groundwater energy and layered media

Head as energy per unit weight. Groundwater flow is driven by spatial variation in hydraulic head. In a saturated pore-water domain the piezometric head is the sum of elevation head z and pressure head u/γw. If a piezometer screen is 3 m below a chosen datum but its water column stands 5 m above the same datum, the measured total head is 5 m; it is not −3 m. Explicit datums prevent false claims about upflow and downflow. Compute the difference between two screened heads and divide by path length to obtain a hydraulic gradient. Flow direction follows decreasing head, although the coordinate sign must be defined.

Darcy discharge and velocity. The magnitude Q=k i A uses the hydraulic conductivity in m/s, dimensionless gradient and total area normal to flow in m². The result is m³/s. If a gross cross-section is 2 m², k=10⁻⁵ m/s and i=0.2, then Q=4×10⁻⁶ m³/s. A travel-time discussion needs the effective pore velocity, which is approximately specific discharge divided by effective porosity under an appropriate flow model. The two velocities should not be labeled interchangeably. Neither relation predicts a contaminant breakthrough curve without dispersion, sorption and boundary conditions.

Layer equivalence. Parallel layers carry flow under a common head gradient; their discharges add. The relevant equivalent horizontal conductivity is thickness-weighted arithmetic mean, kh=ΣkiHi/ΣHi. Perpendicular layers carry the same steady discharge but divide the head loss; the effective vertical conductivity is the thickness-weighted harmonic mean, kv=ΣHi/Σ(Hi/ki). Take a 1 m high-conductivity layer at 10⁻⁴ m/s and a 3 m silty layer at 10⁻⁶ m/s. Then kh=(0.0001×1+0.000001×3)/4=2.575×10⁻⁵ m/s. The vertical equivalent is 4/(1/0.0001+3/0.000001)≈1.329×10⁻⁶ m/s. Treating both paths with the arithmetic mean would drastically overestimate vertical flow.

Stress consequences. Integrate total overburden from unit-weight measurements and subtract the measured or modeled pore pressure. Under simple static groundwater, u=γw(z−zw) below the water level, and zero gauge pressure is often assigned above it. Real partially saturated soils may have suction and perched water. When water heads are not hydrostatic, pressure at a given point must be obtained from its piezometric head rather than from the distance to a visual water line. That correction may alter predicted excavation stability or settlement.

Design boundary. A critical hydraulic gradient is not a comprehensive erosion threshold. Backward erosion piping can evolve along preferential flow paths and filter gradations matter. A simplistic classroom safety ratio using an average gradient should never be used to certify an embankment, excavated base or dam exit. A site-specific seepage geometry, monitoring plan and appropriate erosion assessment are required.

Measurement task. Sketch two piezometers at different elevations and invent physically consistent head readings for upward, downward and zero vertical flow. For each sketch, calculate the average vertical gradient, specify the sign and show what changes if the screen locations are swapped. Explain which observation—the screen depth or the water level in the instrument—determines hydraulic head.

---