11.1 Chemical Fate and Transport Mechanisms (Partitioning, Darcy's Law, Advection)

Chemical Fate and Transport Mechanisms (Partitioning, Darcy's Law, Advection)

Understanding the physical, chemical, and hydrogeological mechanisms governing contaminant migration through environmental media is fundamental to environmental management, site remediation, and regulatory compliance. When hazardous chemicals are released into the subsurface, their ultimate fate, velocity, and spatial distribution are dictated by intermedia phase partitioning, groundwater hydraulics, and solute transport processes.

A Certified Hazardous Materials Manager (CHMM) must master the quantitative formulas and conceptual frameworks used to predict contaminant mobility, calculate groundwater flow velocities, evaluate soil retardation, and delineate Light and Dense Non-Aqueous Phase Liquid (LNAPL/DNAPL) plumes.

---

1. Chemical Partitioning & Intermedia Transport Parameters

When a chemical enters an environmental system, it distributes (partitions) dynamically among four primary environmental phases:

1. Gaseous Phase (soil vapor, atmospheric air)
2. Aqueous Phase (pore water, groundwater, surface water)
3. Solid Phase (soil minerals, natural organic matter)
4. Non-Aqueous Phase Liquid (NAPL) (free-phase pure product)

+-----------------------------------------------------------------------------+
|                   FOUR-PHASE INTERMEDIA PARTITIONING SYSTEM                 |
|                                                                             |
|                      [ GASEOUS PHASE (Soil Vapor) ]                         |
|                               ^         ^
|                Henry's Law   /           \  Vaporization /                  |
|                Constant (Kh)/             \ Volatilization                  |
|                            v               v                                |
|    [ AQUEOUS PHASE (Groundwater) ] <-----> [ NAPL (Free-Phase Product) ]    |
|                            ^      Dissolution                               |
|             Distribution   /                                                |
|             Coeff (Kd)    /                                                 |
|                          v                                                  |
|                 [ SOLID PHASE (Soil Matrix / Organic Carbon) ]              |
+-----------------------------------------------------------------------------+

Vapor Pressure ($P_v$) and Henry's Law Constant ($K_H$)

• Vapor Pressure ($P_v$): The pressure exerted by a chemical's vapor in equilibrium with its pure liquid or solid phase at a given temperature (typically measured in $\text{mm Hg}$, $\text{atm}$, or $\text{Pa}$). High vapor pressure ($> 10,\text{mm Hg}$ at $20^\circ\text{C}$) denotes high volatility (e.g., vinyl chloride, methylene chloride).
• Henry's Law Constant ($K_H$ or $H$): Quantifies the equilibrium partitioning of a chemical between the aqueous phase and the gaseous phase.
  • Dimensional Form: $H = \frac{P_i}{C_w}$ where $P_i$ is the partial pressure in the gas phase ($\text{atm}$) and $C_w$ is the aqueous concentration ($\text{mol/m}^3$ or $\text{mol/L}$). Typical units are $\text{atm}\cdot\text{m}^3/\text{mol}$.
  • Dimensionless Form ($K_H'$ or $H'$): $K_H' = \frac{C_{\text{air}}}{C_{\text{water}}} = \frac{H}{R_u \cdot T}$ where $R_u$ is the universal gas constant ($8.2057 \times 10^{-5},\text{atm}\cdot\text{m}^3/(\text{mol}\cdot\text{K})$) and $T$ is temperature in Kelvin ($K$).

Water Solubility ($S_w$)

Water solubility ($S_w$, expressed in $\text{mg/L}$ or $\text{g/L}$) indicates the maximum mass of a compound that will dissolve in pure water at a specified temperature ($20^\circ\text{C}$ or $25^\circ\text{C}$):

• Hydrophilic Compounds ($S_w > 1,000,\text{mg/L}$): Highly mobile in groundwater, resist soil sorption (e.g., MTBE $S_w \approx 43,000,\text{mg/L}$, Acetone $S_w = \text{miscible}$, Ethanol $S_w = \text{miscible}$). Produce elongated dissolved plumes.
• Hydrophobic Compounds ($S_w < 10,\text{mg/L}$): Low aqueous solubility, tend to partition into soil organic carbon or accumulate as pure NAPL (e.g., Benzo[a]pyrene $S_w = 0.0016,\text{mg/L}$, PCBs $S_w < 0.1,\text{mg/L}$). Slow plume propagation; high sorption.

Octanol-Water Partition Coefficient ($K_{ow}$ and $\log K_{ow}$)

The Octanol-Water Partition Coefficient is the ratio of chemical concentration in $n$-octanol (representing lipid/organic phases) to its concentration in water at equilibrium: $K_{ow} = \frac{C_{\text{octanol}}}{C_{\text{water}}}$ Because values span many orders of magnitude, it is expressed logarithmically as $\mathbf{\log K_{ow}}$:

• $\mathbf{\log K_{ow} < 1.0}$: Hydrophilic, highly water-soluble, minimal sorption to soils, low bioaccumulation potential.
• $\mathbf{\log K_{ow} = 1.0\text{ to }3.0}$: Moderately hydrophobic, mobile in sandy aquifers with low organic carbon.
• $\mathbf{\log K_{ow} > 3.0}$: Strongly lipophilic, hydrophobic, high sorption to soil organic matter, high bioaccumulation factor in aquatic organisms and human adipose tissue (e.g., DDT $\log K_{ow} \approx 6.9$, PCBs $\log K_{ow} \approx 5.5\text{--}7.5$, Dioxins $\log K_{ow} \approx 6.8$).

Soil Organic Carbon Partition Coefficient ($K_{oc}$) and Distribution Coefficient ($K_d$)

• $K_{oc}$ ($\text{mL/g}$ or $\text{L/kg}$): The chemical-specific partition coefficient normalized to organic carbon content, quantifying a chemical's intrinsic affinity to bind to organic carbon in soil or aquifer sediment.
• $K_d$ (Soil-Water Distribution Coefficient, $\text{mL/g}$ or $\text{L/kg}$): The linear adsorption coefficient representing the ratio of contaminant mass sorbed to solid soil ($C_s$, $\text{mg/kg}$) to dissolved aqueous concentration ($C_w$, $\text{mg/L}$): $K_d = \frac{C_s}{C_w} = K_{oc} \cdot f_{oc}$ where $\mathbf{f_{oc}}$ is the mass fraction of natural organic carbon in the soil matrix (dimensionless, e.g., $1%\text{ organic carbon} \implies f_{oc} = 0.01$).

[!NOTE] In mineral soils with very low organic carbon ($f_{oc} < 0.001$ or $0.1%$), sorption is governed primarily by direct mineral surface interactions (clay minerals, iron/aluminum oxyhydroxides) rather than organic carbon partitioning.

Contaminant Retardation Factor ($R$) and Solute Velocity ($v_c$)

As groundwater transports dissolved contaminants through porous media, reversible sorption onto soil organic matter slows contaminant migration relative to the bulk groundwater flow. This process is quantified by the Retardation Factor ($R$):

$R = 1 + \frac{\rho_b \cdot K_d}{\theta} = 1 + \frac{\rho_b \cdot (K_{oc} \cdot f_{oc})}{\theta}$

Where:

• $R$ = Retardation factor (dimensionless; $\mathbf{R \ge 1.0}$)
• $\rho_b$ = Dry soil bulk density (typically $1.4\text{ to }1.8,\text{g/cm}^3$ or $\text{kg/L}$)
• $K_d$ = Soil-water distribution coefficient ($\text{cm}^3/\text{g}$ or $\text{L/kg}$)
• $\theta$ (or $n$) = Total or effective porosity (dimensionless, typically $0.20\text{ to }0.40$)

The linear contaminant plume migration velocity ($v_c$) is calculated by dividing the average linear groundwater velocity ($v_w$ or $v_x$) by the Retardation Factor:

$v_c = \frac{v_w}{R}$

• When $R = 1.0$, the contaminant does not sorb (conservative tracer, e.g., Chloride, Bromide, Tritium) and migrates at the exact velocity of groundwater ($v_c = v_w$).
• When $R = 10.0$, the contaminant plume advances at only $\frac{1}{10}\text{th}$ ($10%$) of the groundwater seepage velocity.

+-----------------------------------------------------------------------------+
|                 WORKED EXAMPLE: RETARDATION & PLUME VELOCITY                |
|                                                                             |
|   Given:                                                                    |
|   - Contaminant: Trichloroethene (TCE) with Koc = 126 mL/g                 |
|   - Aquifer Soil Organic Carbon fraction: foc = 0.005 (0.5%)                |
|   - Dry Soil Bulk Density: rho_b = 1.60 g/cm^3                              |
|   - Effective Porosity: theta = 0.25 (25%)                                  |
|   - Average Linear Groundwater Velocity: vw = 120 ft/year                   |
|                                                                             |
|   Step 1: Calculate Distribution Coefficient (Kd)                           |
|     Kd = Koc * foc = (126 mL/g) * 0.005 = 0.63 mL/g (or cm^3/g)            |
|                                                                             |
|   Step 2: Calculate Retardation Factor (R)                                  |
|     R = 1 + (rho_b * Kd) / theta                                            |
|     R = 1 + (1.60 g/cm^3 * 0.63 cm^3/g) / 0.25                              |
|     R = 1 + (1.008) / 0.25 = 1 + 4.032 = 5.032                              |
|                                                                             |
|   Step 3: Calculate Contaminant Plume Velocity (vc)                         |
|     vc = vw / R = (120 ft/year) / 5.032 = 23.85 ft/year                     |
|                                                                             |
|   Conclusion: While groundwater travels 120 ft/year, the TCE plume         |
|   advances at only ~23.85 ft/year due to soil sorption.                     |
+-----------------------------------------------------------------------------+

---

2. Groundwater Hydrogeology & Transport Equations

Groundwater flow in saturated porous media is governed by fundamental hydrogeologic laws formulated by Henry Darcy in 1856.

+-----------------------------------------------------------------------------+
|                          DARCY'S LAW SCHEMATIC                              |
|                                                                             |
|     Piezometer 1                                 Piezometer 2               |
|        |                                            |                       |
|      [ h1 ]                                         |                       |
|        |     \                                      |                       |
|        |      \  Hydraulic Gradient (i = dh / L)    |                       |
|        |       \                                  [ h2 ]                    |
|        |========\===================================|                       |
|        |                                            |                       |
|   +----+--------------------------------------------+----+                  |
|   |    |-----> Flow Q                               |    |                  |
|   |    |       Porous Medium Cross-Section Area (A) |    | Aquifer          |
|   |    |       Hydraulic Conductivity (K)           |    | Thickness (b)    |
|   +----+--------------------------------------------+----+                  |
|        |<------------------- L -------------------->|                       |
+-----------------------------------------------------------------------------+

Darcy's Law

Darcy's Law states that the volumetric discharge rate ($Q$) of water through a porous medium is directly proportional to the hydraulic conductivity ($K$), cross-sectional area ($A$), and hydraulic gradient ($i$):

$Q = -K \cdot A \cdot \frac{dh}{dl} = -K \cdot A \cdot i$

Where:

• $Q$ = Volumetric flow rate / discharge ($\text{m}^3/\text{day}$ or $\text{ft}^3/\text{day}$)
• $K$ = Hydraulic conductivity ($\text{m/day}$, $\text{cm/s}$, or $\text{ft/day}$)
• $A$ = Cross-sectional area perpendicular to groundwater flow ($\text{m}^2$ or $\text{ft}^2$; where $A = W \cdot b$, width $\times$ aquifer thickness)
• $i = \frac{dh}{dl} = \frac{h_1 - h_2}{L}$ = Hydraulic gradient (dimensionless ratio of head loss over distance $L$)

Darcy Velocity (Specific Discharge, $q$)

The Darcy velocity ($q$, also called specific discharge or Darcy flux) represents the apparent volumetric discharge per unit total cross-sectional area:

$q = \frac{Q}{A} = -K \cdot i$

[!WARNING] Critical Exam Distinction: Darcy velocity ($q$) is a fictitious velocity because it assumes water flows uniformly across the entire cross-sectional area (both solid grains and pores). In reality, water flows only through connected, open pore channels.

Average Linear Groundwater Seepage Velocity ($v_x$ or $v_w$)

To calculate the true physical velocity at which water and non-sorbing dissolved solutes travel through pore spaces, Darcy velocity must be divided by the effective porosity ($n_e$):

$v_x = \frac{q}{n_e} = \frac{K \cdot i}{n_e} = \frac{K}{n_e} \cdot \frac{dh}{dl}$

Where:

• $v_x$ = Average linear seepage velocity ($\text{m/day}$, $\text{cm/s}$, or $\text{ft/day}$)
• $n_e$ = Effective porosity (dimensionless, representing the fraction of interconnected, drainable pore space; $n_e < n_{\text{total}}$)

+-----------------------------------------------------------------------------+
|                  WORKED EXAMPLE: SEEPAGE VELOCITY & TRAVEL TIME             |
|                                                                             |
|   Given:                                                                    |
|   - Hydraulic conductivity: K = 2.0 * 10^-3 cm/s                            |
|   - Upgradient head h1 = 105.0 ft, Downgradient head h2 = 99.0 ft          |
|   - Distance between wells: L = 300 ft                                      |
|   - Effective porosity: ne = 0.20 (20%)                                     |
|   - Distance to nearest municipal drinking water well: D = 1,200 ft         |
|                                                                             |
|   Step 1: Calculate Hydraulic Gradient (i)                                  |
|     i = (h1 - h2) / L = (105.0 ft - 99.0 ft) / 300 ft = 6.0 / 300 = 0.020   |
|                                                                             |
|   Step 2: Convert K to ft/day                                               |
|     1 cm/s = 2,834.65 ft/day                                                |
|     K = (2.0 * 10^-3 cm/s) * 2,834.65 (ft/day)/(cm/s) = 5.67 ft/day         |
|                                                                             |
|   Step 3: Calculate Average Linear Seepage Velocity (vx)                    |
|     vx = (K * i) / ne = (5.67 ft/day * 0.020) / 0.20                        |
|     vx = 0.1134 / 0.20 = 0.567 ft/day (approx 207 ft/year)                  |
|                                                                             |
|   Step 4: Calculate Groundwater Travel Time (t) to Municipal Well           |
|     t = D / vx = 1,200 ft / (0.567 ft/day) = 2,116 days (approx 5.8 years) |
+-----------------------------------------------------------------------------+

---

3. Solute Fate and Transport Processes in Porous Media

The 3D transport of dissolved chemical species through saturated groundwater systems is mathematically modeled by the Advection-Dispersion Equation (ADE):

$\frac{\partial C}{\partial t} = \nabla \cdot (D \nabla C) - \nabla \cdot (v_x C) - \frac{\rho_b}{\theta}\frac{\partial C_s}{\partial t} \pm R_{\text{rxn}}$

This equation couples five distinct physical-chemical mechanisms:

+-----------------------------------------------------------------------------+
|                   FIVE CONTAMINANT TRANSPORT MECHANISMS                     |
|                                                                             |
|   1. ADVECTION        : Transport via bulk fluid flow along hydraulic       |
|                         gradient at velocity vx.                            |
|   2. MECHANICAL       : Spreading due to velocity variations within         |
|      DISPERSION         tortuous pore channels (longitudinal & transverse). |
|   3. MOLECULAR        : Spreading along concentration gradients according   |
|      DIFFUSION          to Fick's Laws (dominant in low-permeability clays).|
|   4. ADSORPTION /     : Partitioning of solute onto solid soil particles,   |
|      DESORPTION         causing plume retardation (Kd, R).                  |
|   5. CHEMICAL / BIO   : Abiotic reactions (hydrolysis, redox) and microbial |
|      TRANSFORMATION     biodegradation (first-order decay).                 |
+-----------------------------------------------------------------------------+

1. Advection

The physical movement of dissolved solute mass along with the bulk flowing groundwater at average linear velocity $v_x$. Advection carries the center of mass of the contaminant plume.

2. Mechanical Dispersion

Mixing and spreading of the solute front caused by microscopic velocity variations in the porous matrix. Dispersion occurs because:

• Fluid travels faster in the center of pore throats than along friction-bound grain walls.
• Pore channels vary in diameter.
• Fluid particles travel along tortuous, divergent pathways.
• Longitudinal dispersivity ($\alpha_L$) produces spreading along the flow direction; Transverse dispersivity ($\alpha_T$) produces lateral and vertical plume spreading.

3. Molecular Diffusion (Fick's Laws)

The spontaneous movement of chemical molecules from regions of higher concentration to lower concentration due to random thermal Brownian motion, independent of bulk fluid velocity:

• Fick's First Law (Steady-State Flux): $J = -D_e \cdot \frac{dC}{dx}$ where $J$ is mass flux ($\text{mg}/(\text{m}^2\cdot\text{s})$), $D_e$ is the effective diffusion coefficient in porous media ($D_e = \omega \cdot D_0$, accounting for tortuosity $\omega$), and $\frac{dC}{dx}$ is concentration gradient.
• Fick's Second Law (Transient Diffusion): $\frac{\partial C}{\partial t} = D_e \cdot \frac{\partial^2 C}{\partial x^2}$

[!NOTE] In high-permeability sands and gravels, mechanical dispersion completely dominates molecular diffusion. In low-permeability clays, shales, and aquitards where groundwater seepage velocity approaches zero ($v_x \approx 0$), molecular diffusion is the sole mechanism of contaminant transport and matrix back-diffusion.

4. Hydrodynamic Dispersion ($D_L$)

In field hydrogeology, mechanical dispersion and molecular diffusion are combined into the coefficient of hydrodynamic dispersion ($D_L$): $D_L = \alpha_L \cdot v_x + D^*$

5. Biological & Chemical Transformation (Decay Kinetics)

Many organic contaminants degrade via microbial metabolism (biodegradation) or abiotic chemical reactions (hydrolysis, oxidation-reduction). Most environmental transformations follow first-order decay kinetics:

$C(t) = C_0 \cdot e^{-k \cdot t}$

Where:

• $C(t)$ = Concentration at time $t$
• $C_0$ = Initial concentration at $t = 0$
• $k$ = First-order degradation rate constant ($\text{time}^{-1}$, e.g., $\text{day}^{-1}$ or $\text{year}^{-1}$)
• Half-life ($t_{1/2}$): The time required for $50%$ of the contaminant mass to degrade: $t_{1/2} = \frac{\ln(2)}{k} = \frac{0.693}{k}$

---

4. Non-Aqueous Phase Liquids (NAPLs): LNAPL vs. DNAPL Dynamics

When immiscible liquid hydrocarbons or chlorinated solvents are released, they do not initially dissolve; they migrate as separate-phase Non-Aqueous Phase Liquids (NAPLs). Their subsurface fate is governed by their Specific Gravity (SG), interfacial tension, viscosity, and capillary entry pressure.

+-----------------------------------------------------------------------------------------+
|                         LNAPL VS DNAPL SUBSURFACE PROFILES                              |
|                                                                                         |
|               [ LNAPL RELEASE (Gasoline/BTEX) ]    [ DNAPL RELEASE (TCE/PCE) ]          |
|                             |                                    |                      |
|  VADOSE ZONE                v                                    v                      |
|  (Unsaturated)       Residual Droplets                    Residual Ganglia              |
|                             |                                    |                      |
|  ===========================v====================================v==================== |
|  WATER TABLE          [ LNAPL POOL ]                             |                      |
|  (Capillary Fringe)  Floats & Smears                             |                      |
|                     (Smear Zone forms)                           v                      |
|                                                        Penetrates Saturated             |
|  SATURATED ZONE                                        Zone via Gravity                 |
|  (Groundwater Flow)  Dissolved Plume ======>                     |                      |
|                                                                  v                      |
|                                                        Dissolved Plume ======>          |
|                                                                  |                      |
|  ================================================================v==================== |
|  BASAL AQUITARD                                            [ DNAPL POOL ]               |
|  (Impermeable Clay)                                     Accumulates in Troughs          |
+-----------------------------------------------------------------------------------------+

Light Non-Aqueous Phase Liquids (LNAPLs)

• Physical Criterion: Specific gravity less than water ($\mathbf{\text{SG} < 1.00}$, density $< 1.0,\text{g/cm}^3$).
• Representative Compounds: Gasoline, diesel fuel, jet fuel (JP-4, JP-8), motor oil, benzene ($\text{SG} = 0.88$), toluene ($\text{SG} = 0.87$), ethylbenzene ($\text{SG} = 0.87$), xylenes ($\text{SG} = 0.86\text{--}0.88$) (BTEX).
• Migration Behavior:
  • Released at surface, migrates downward through the unsaturated vadose zone under gravity, leaving behind residual globules.
  • Upon reaching the capillary fringe and saturated water table, the buoyancy of the lighter liquid prevents significant penetration into the water table.
  • Spreads radially and laterally atop the water table, forming an LNAPL "pancake" or mobile free-phase lens.
  • As seasonal water tables rise and fall, the free-phase LNAPL smears vertically across soil pores, creating a thick, highly contaminated smear zone that acts as a continuous long-term secondary source of dissolved groundwater contamination.

Dense Non-Aqueous Phase Liquids (DNAPLs)

• Physical Criterion: Specific gravity greater than water ($\mathbf{\text{SG} > 1.00}$, density $> 1.0,\text{g/cm}^3$).
• Representative Compounds: Tetrachloroethene / PCE ($\text{SG} = 1.62$), Trichloroethene / TCE ($\text{SG} = 1.46$), 1,1,1-Trichloroethane ($\text{SG} = 1.34$), Methylene chloride ($\text{SG} = 1.33$), Carbon tetrachloride ($\text{SG} = 1.59$), Chloroform ($\text{SG} = 1.49$), Creosote / Coal tar ($\text{SG} = 1.05\text{--}1.15$), Polychlorinated Biphenyls / PCBs ($\text{SG} = 1.20\text{--}1.55$).
• Migration Behavior:
  • High density and low kinematic viscosity enable DNAPLs to overcome the capillary entry pressure of water-saturated pores.
  • Sinks vertically downward directly through the water table and through the entire saturated aquifer thickness, leaving a trail of disconnected residual droplets (ganglia) held by capillary forces.
  • Downward vertical migration continues until the DNAPL encounters an impermeable confining layer (clay aquitard, competent bedrock surface).
  • Critical Structural Control: DNAPL pools in bedrock or clay troughs and migrates along the down-dip slope of the confining geologic unit, which is often completely opposite to the direction of groundwater hydraulic gradient.