← OpenSTEM Portal / Environmental Chemistry
⏱️ 00:00
Chapter 7 • Theory & Derivations

Soil Chemistry, Clay Mineralogy & Nutrient Cycles

Soil mineralogy, phyllosilicate clay structures, and isomorphous substitution; cation exchange capacity (CEC), diffuse double layer theory, and acid sulfate soils; geochemical nitrogen and phosphorus dynamics; humic substances, soil organic matter, heavy metal sorption hysteresis, and phytoremediation engineering.

§7.1 Physical and Mineralogical Architecture of Soil: Horizonation & Aggregation

Soil represents a heterogeneous, multiphase biogeochemical reactor comprising solid mineral particles ($45\%$), soil organic matter ($5\%$), porewater ($25\%$), and soil air ($25\%$ by volume under ideal conditions).

Pedogenic Horizonation

Vertical pedogenesis differentiates the soil profile into distinctive genetic master horizons:

1. O Horizon: Uppermost organic layer composed of fresh, decomposing, or humified plant debris (litter, duff).

2. A Horizon (Topsoil): Mineral horizon characterized by intense accumulation of humified organic matter, dark coloration, high microbial biomass, and maximum biological activity.

3. E Horizon (Eluvial layer): Zone of maximum leaching or eluviation of silicate clay, iron, and aluminum oxides, leaving an accumulation of resistant quartz and silt grains (frequently ashen or light grey).

4. B Horizon (Subsoil / Illuvial layer): Zone of illuvial accumulation of translocated clays, sesquioxides ($Fe_2O_3, Al_2O_3$), carbonates, or humus.

5. C Horizon: Unconsolidated parent material minimally affected by pedogenic processes.

6. R Horizon: Continuous, hard bedrock foundation.

Soil Texture & The USDA Soil Triangle

Soil mineral particles are categorized into three distinct size fractions:

  • Sand: $0.05 - 2.0\ \text{mm}$ diameter (predominantly inert quartz, high hydraulic conductivity, zero plastic behavior).
  • Silt: $0.002 - 0.05\ \text{mm}$ diameter (weatherable primary minerals, moderate water holding capacity).
  • Clay: $< 0.002\ \text{mm}$ ($< 2\ \mu\text{m}$) diameter (colloidal secondary minerals with massive specific surface areas ranging from $10$ to $> 800\ \text{m}^2/\text{g}$ and permanent electrostatic charges).

§7.2 Clay Mineralogy: Phyllosilicates & Isomorphous Substitution

Soil clays are secondary phyllosilicates (sheet silicates) assembled from two fundamental crystallographic building blocks:

1. Tetrahedral Sheets ($T$): Silicon cations ($Si^{4+}$) coordinated tetrahedrally to four oxygen atoms ($SiO_4^{4-}$), sharing three basal oxygens to form a hexagonal planar mesh.

2. Octahedral Sheets ($O$): Aluminum ($Al^{3+}$) or magnesium ($Mg^{2+}$) cations coordinated octahedrally to six hydroxyl groups or apex oxygens. Minerals with trivalent $Al^{3+}$ occupying two out of three cation sites are *dioctahedral* (e.g., gibbsite); those with divalent $Mg^{2+}$ occupying all three sites are *trioctahedral* (e.g., brucite).

1:1 vs 2:1 Layer Types

  • 1:1 Layer Silicates (Kaolinite group):
  • Consist of alternating single tetrahedral and octahedral sheets ($T-O$).
  • Consecutive 1:1 layers are held tightly together by rigid interlayer hydrogen bonds between octahedral hydroxyls and basal oxygens of adjacent tetrahedral sheets.
  • Basal spacing is fixed at $0.72\ \text{nm}$ ($7.2\ \text{\AA}$).
  • Non-swelling, low specific surface area ($10 - 30\ \text{m}^2/\text{g}$), and low cation exchange capacity ($3 - 15\ \text{cmol}_c/\text{kg}$) originating purely from broken crystal edge hydroxyls.
  • 2:1 Layer Silicates (Smectites, Vermiculites, Illites):
  • Formed by one octahedral sheet sandwiched between two tetrahedral sheets ($T-O-T$).
  • Interlayer spaces lack hydrogen bonds, allowing hydration and cation insertion.
  • Smectite (Montmorillonite): Displays extensive isomorphous substitution in the octahedral sheet ($Mg^{2+}$ substituting for $Al^{3+}$). Layers readily expand with water (basal spacing swells from $0.96\ \text{nm}$ to $> 1.8\ \text{nm}$); huge surface area ($600 - 800\ \text{m}^2/\text{g}$) and high CEC ($80 - 150\ \text{cmol}_c/\text{kg}$).
  • Illite: Significant tetrahedral substitution of $Al^{3+}$ for $Si^{4+}$. Unhydrated potassium ions ($K^+$) fix firmly into hexagonal cavities, locking layers together ($1.0\ \text{nm}$ fixed spacing).

Isomorphous Substitution: Origin of Permanent Structural Charge

Isomorphous substitution occurs when a cation of similar ionic radius but lower valence replaces a structural cation without disrupting crystal symmetry:

\[Si^{4+} \xrightarrow{\text{Substituted by}} Al^{3+} \quad (\text{Tetrahedral deficit: } -1)\]
\[Al^{3+} \xrightarrow{\text{Substituted by}} Mg^{2+} / Fe^{2+} \quad (\text{Octahedral deficit: } -1)\]

Because the replacement leaves an uncompensated negative charge in the crystal framework, 2:1 clays possess a permanent, pH-independent negative surface charge.

§7.3 Cation Exchange Capacity (CEC), Base Saturation & Gouy-Chapman Theory

The negative charge on soil colloids is balanced by a swarm of electrostatically held exchangeable cations ($Ca^{2+}, Mg^{2+}, K^+, Na^+, Al^{3+}, H^+$).

Cation Exchange Capacity (CEC) & Base Saturation

Cation Exchange Capacity (CEC) is the total quantity of exchangeable cations that a unit mass of soil dry matter can adsorb at a specified pH, expressed in centimoles of positive charge per kilogram ($\text{cmol}_c/\text{kg}$ or $\text{meq}/100\text{g}$):

\[\text{CEC} = [Ca^{2+}] + [Mg^{2+}] + [K^+] + [Na^+] + [Al^{3+}] + [H^+] \quad (\text{in } \text{cmol}_c/\text{kg})\]

Cations are partitioned into basic cations ($Ca^{2+}, Mg^{2+}, K^+, Na^+$) and acidic cations ($Al^{3+}, H^+$). The Base Saturation Percentage (BS%) is:

\[\text{BS\%} = \frac{[Ca^{2+}] + [Mg^{2+}] + [K^+] + [Na^+]}{\text{CEC}} \times 100\]

High base saturation ($> 60 - 80\%$) indicates fertile, circumneutral soils; low base saturation ($< 35\%$) characterizes strongly leached, acidic soils where toxic $Al^{3+}$ dominates the exchange complex.

Lyotropic Series (Exchange Selectivity)

Colloidal binding affinity increases with increasing ionic charge and decreasing hydrated ionic radius:

\[Al^{3+} > Ca^{2+} > Mg^{2+} > K^+ \approx NH_4^+ > Na^+\]

Gouy-Chapman Diffuse Double Layer (DDL) Theory

At the charged clay-solution interface, an electrical double layer forms comprising the fixed negative surface charge ($\sigma_0$) and a diffuse swarm of counter-ions. The electrostatic potential $\psi(x)$ decays exponentially away from the surface according to the linearized Poisson-Boltzmann equation:

\[\psi(x) = \psi_0 \exp(-\kappa x)\]

where $\kappa^{-1}$ is the Debye screening length (the characteristic thickness of the electrical double layer):

\[\kappa = \sqrt{\frac{2000 F^2 I}{\epsilon_0 \epsilon_r R T}} \implies \kappa^{-1} = \sqrt{\frac{\epsilon_0 \epsilon_r R T}{2000 F^2 I}}\]

where $I = \frac{1}{2}\sum c_i z_i^2$ is the porewater ionic strength. Increasing ionic strength or adding multivalent cations compresses the double layer (decreases $\kappa^{-1}$), overcoming electrostatic repulsion and triggering colloidal flocculation.

§7.4 Soil pH Dynamics, Buffering Mechanisms & Acid Sulfate Soils

Soil pH governs nutrient bioavailability, microbial respiration, and mineral dissolution rates.

Soil Buffering Capacities Across pH Ranges

Soils resist shifts in pH through distinct thermodynamic buffer mechanisms:

1. Carbonate Buffer (pH 7.5 - 8.5): Controlled by calcite dissolution and carbonic acid equilibria:

\[\text{CaCO}_3\text{(s)} + 2\text{H}^+ \rightleftharpoons \text{Ca}^{2+} + \text{CO}_2\text{(g)} + \text{H}_2\text{O}\]

2. Silicate and Clay Edge Buffer (pH 5.5 - 7.0): Controlled by protonation/deprotonation of variable-charge silanol ($\equiv\!SiOH$) and aluminol ($\equiv\!AlOH$) groups.

3. Cation Exchange Buffer (pH 4.5 - 5.5): Exchange of soil solution $H^+$ for basic cations on clay interlayers.

4. Aluminum Oxide / Hydroxide Buffer (pH < 4.5): Under extreme acidity, toxic octahedral aluminum dissolves:

\[\text{Al(OH)}_3\text{(s)} + 3\text{H}^+ \rightleftharpoons \text{Al}^{3+} + 3\text{H}_2\text{O}\]

Free $Al^{3+}$ hydrolyzes in water, releasing more protons ($Al^{3+} + H_2O \rightleftharpoons Al(OH)^{2+} + H^+$), stunting plant root elongation.

Acid Sulfate Soils (Cat-clays)

In coastal mangrove swamps and deltaic lowlands (such as the coastal belt of Bangladesh), tidal inundation deposits sedimentary pyrite ($FeS_2$) formed under reducing conditions. When these waterlogged soils are drained or excavated for aquaculture, oxygen enters the soil pore network, initiating rapid biogeochemical oxidation catalyzed by Acidithiobacillus ferrooxidans:

\[2\text{FeS}_2\text{(s)} + 7\text{O}_2 + 2\text{H}_2\text{O} \rightarrow 2\text{Fe}^{2+} + 4\text{SO}_4^{2-} + 4\text{H}^+\]
\[4\text{Fe}^{2+} + \text{O}_2 + 4\text{H}^+ \rightarrow 4\text{Fe}^{3+} + 2\text{H}_2\text{O}\]
\[\text{FeS}_2\text{(s)} + 14\text{Fe}^{3+} + 8\text{H}_2\text{O} \rightarrow 15\text{Fe}^{2+} + 2\text{SO}_4^{2-} + 16\text{H}^+\]

The overall stoichiometry yields 4 moles of $H^+$ per mole of pyrite oxidized:

\[4\text{FeS}_2 + 15\text{O}_2 + 14\text{H}_2\text{O} \rightarrow 4\text{Fe(OH)}_3\text{(s)} + 8\text{H}_2\text{SO}_4\]

Soil pH plummets below $2.5 - 3.5$, dissolving clays, releasing toxic aluminum and jarosite ($KFe_3(SO_4)_2(OH)_6$), and destroying vegetative cover.

§7.5 Global Biogeochemical Nitrogen Cycle in Soils

Nitrogen ($N$) is an essential plant macronutrient existing in oxidation states from $-3$ ($NH_3, NH_4^+$) to $+5$ ($NO_3^-$).

Key Soil Nitrogen Transformations

1. Biological Nitrogen Fixation (BNF):

  • Atmospheric $N_2$ (oxidation state $0$) is reduced to ammonia ($-3$) by the bacterial enzyme nitrogenase:
\[N_2 + 8H^+ + 8e^- + 16\text{ATP} \xrightarrow{\text{Nitrogenase}} 2NH_3 + H_2 + 16\text{ADP} + 16P_i\]
  • Mediated by symbiotic rhizobia (Rhizobium, Bradyrhizobium) in legume root nodules and free-living diazotrophs (Azotobacter, Clostridium).

2. Ammonification (Mineralization):

  • Heterotrophic microbial decomposition of organic nitrogen compounds (proteins, amino acids, nucleic acids) into ammonium ($NH_4^+$).

3. Nitrification:

  • Two-step aerobic chemoautotrophic oxidation of ammonium to nitrate:
  • Step 1 (Ammonia oxidation): Catalyzed by Nitrosomonas and ammonia-oxidizing archaea (AOA):
\[2NH_4^+ + 3O_2 \rightarrow 2NO_2^- + 4H^+ + 2H_2O\]
  • Step 2 (Nitrite oxidation): Catalyzed by Nitrobacter and Nitrospira:
\[2NO_2^- + O_2 \rightarrow 2NO_3^-\]
  • Acidifying impact: Nitrification generates two moles of $H^+$ per mole of $NH_4^+$ oxidized, serving as a primary driver of agricultural soil acidification under ammonium fertilizer use.

4. Denitrification:

  • Anaerobic microbial respiration using nitrate as an alternative terminal electron acceptor:
\[NO_3^- \rightarrow NO_2^- \rightarrow NO \rightarrow N_2O \rightarrow N_2\]
  • Promoted by high organic carbon, low redox potential ($E_h < +200\ \text{mV}$), and waterlogged soils. Incomplete denitrification emits nitrous oxide ($N_2O$), a potent greenhouse gas.

5. Anammox (Anaerobic Ammonium Oxidation):

  • Specialized planctomycete bacteria oxidize ammonium directly using nitrite under anoxic conditions:
\[NH_4^+ + NO_2^- \rightarrow N_2 + 2H_2O\]

§7.6 Phosphorus Geochemistry in Soils: Fixation, Occlusion & Eutrophication

Unlike nitrogen, phosphorus ($P$) lacks a significant atmospheric gaseous phase; its biogeochemical cycle is exclusively sedimentary and lithospheric.

Soil Phosphorus Speciation & pH-Dependent Fixation

Soil solution orthophosphate ($H_2PO_4^-, HPO_4^{2-}$) typically exists at very low concentrations ($0.01 - 1.0\ \mu\text{mol}/\text{L}$) due to rapid geochemical fixation:

1. Acidic Soils (pH < 5.5):

  • Phosphate precipitates as insoluble iron and aluminum phosphates (strengite, $FePO_4 \cdot 2H_2O$; variscite, $AlPO_4 \cdot 2H_2O$):
\[Al^{3+} + H_2PO_4^- + 2H_2O \rightleftharpoons AlPO_4 \cdot 2H_2O\text{(s)} + 2H^+\]
  • Forms strong bidentate binuclear inner-sphere surface complexes on goethite and gibbsite mineral edges.

2. Alkaline and Calcareous Soils (pH > 7.5):

  • Phosphate precipitates as calcium phosphate minerals, progressing from dicalcium phosphate ($CaHPO_4$) to octacalcium phosphate ($Ca_8H_2(PO_4)_6 \cdot 5H_2O$) and ultimately insoluble hydroxyapatite:
\[10Ca^{2+} + 6HPO_4^{2-} + 2H_2O \rightarrow Ca_{10}(PO_4)_6(OH)_2\text{(s)} + 8H^+\]

3. Optimum Phosphorus Bioavailability:

  • Maximum orthophosphate availability occurs in the narrow circumneutral window of pH 6.2 to 6.8, where both Al/Fe fixation and calcium precipitation are minimized.

Phosphorus Occlusion & Agricultural Runoff

Over pedogenic time scales, phosphate is incorporated into amorphous iron and aluminum oxide coatings, forming occluded phosphorus that is completely inaccessible to plant roots. Excessive broadcast application of commercial phosphate fertilizers saturates topsoil sorption capacities, causing particulate and dissolved reactive phosphorus (DRP) runoff into surface waters, triggering catastrophic algal blooms and aquatic eutrophication.

§7.7 Soil Organic Matter (SOM), Humic Substances & Carbon Sequestration

Soil organic matter (SOM) represents the largest terrestrial organic carbon reservoir, storing over $1500\ \text{Pg C}$ (more than double the atmospheric carbon pool).

Composition and Fractionation of SOM

SOM is partitioned into labile and recalcitrant fractions:

1. Labile (Active) Pool: Microbial biomass, simple carbohydrates, proteins, and root exudates with turnover times of days to years.

2. Slow / Intermediate Pool: Chemically protected particulate organic matter and microaggregate-associated carbon (turnover: decades).

3. Passive (Recalcitrant) Pool / Humic Substances: Highly condensed, cross-linked aromatic macromolecules resistant to enzymatic cleavage (turnover: centuries to millennia).

Classical Humic Substance Fractions

Humic substances are historically differentiated based on aqueous solubility across acid-base regimes:

  • Fulvic Acid (FA):
  • Soluble in both alkali and acid across all pH conditions.
  • Lower molecular weight ($1000 - 5000\ \text{Da}$), light yellow to brown color.
  • High oxygen content, high density of carboxylic ($-COOH$) and phenolic ($-OH$) functional groups, high acidity ($900 - 1400\ \text{cmol}_c/\text{kg}$), and exceptional metal chelation capacity.
  • Humic Acid (HA):
  • Soluble in dilute alkali ($ ext{pH} > 2$), but precipitates as dark brown to black flocs upon acidification ($ ext{pH} < 2$).
  • Intermediate to high molecular weight ($10,000 - 100,000\ \text{Da}$), rich in aromatic rings, quinones, and heterocyclic nitrogen.
  • Humin:
  • Insoluble in both alkali and acid.
  • Tightly bound to clay mineral lattices through organo-mineral associations; extremely recalcitrant.

Terrestrial Carbon Sequestration

Enhancing SOM via conservation tillage, cover cropping, and biochar amendment sequesters atmospheric $CO_2$ while improving soil aggregate stability, water-holding capacity, and cation exchange capacity.

§7.8 Soil Contamination, Heavy Metal Sorption Hysteresis & Phytoremediation

Toxic trace metals ($Pb, Cd, As, Cr, Cu, Zn$) accumulate in soils due to industrial effluents, mining spoils, electroplating discharges, and agrochemicals.

Sorption-Desorption Hysteresis

Heavy metal retention on soil mineral and humic phases is rarely fully reversible. The desorption isotherm typically diverges from the sorption isotherm, displaying hysteresis:

\[K_{desorption} \neq K_{adsorption}\]

Hysteresis is driven by:

  1. Intra-particle micropore diffusion and diffusion into interlayer spaces.
  2. Transition from weak outer-sphere electrostatic complexes to covalent inner-sphere coordination.
  3. Solid-state chemical entrapment and precipitation as secondary solid phases.

Heavy Metal Bioavailability & Sequential Extraction

Total metal concentration is a poor predictor of ecotoxicity; biological uptake depends on chemical speciation. The Tessier Sequential Extraction Procedure partitions soil metals into five operationally defined fractions:

  1. Exchangeable (extracted with $MgCl_2$ or $NaOAc$, most bioavailable).
  2. Bound to Carbonates (extracted with $NaOAc$ at pH 5.0).
  3. Bound to Iron and Manganese Oxides (extracted with hydroxylamine hydrochloride, $NH_2OH\cdot HCl$).
  4. Bound to Organic Matter (extracted with acidified $H_2O_2$).
  5. Residual / Mineral Matrix (extracted with $HF-HClO_4$, non-bioavailable).

Phytoremediation Engineering

Phytoremediation utilizes hyperaccumulator plants to remediate contaminated soils:

  • Phytoextraction: Translocation of metals from roots to harvestable above-ground shoots:
\[\text{Bioaccumulation Factor (BCF)} = \frac{C_{root}}{C_{soil}}, \quad \text{Translocation Factor (TF)} = \frac{C_{shoot}}{C_{root}}\]

A true hyperaccumulator exhibits both $\text{BCF} > 1$ and $\text{TF} > 1$.

  • Phytostabilization: Precipitation or immobilization of metals in the rhizosphere via root exudates, preventing leaching into groundwater.

Worked Practice Problems (9 Challenge Exercises)

Multi-step solved problems covering barometric scale height, Leighton photostationary ozone equilibria, VOC-hydroxyl radical degradation lifetimes, aqueous SO2 bisulfite oxidation, acid rain carbonate dissolution, three-way catalytic converter stoichiometry, aerosol Stokes terminal settling, and atmospheric radon-222 secular equilibrium with line-by-line mathematical proofs.

Foundational Example 7.1: Effective Cation Exchange Capacity and Base Saturation Calculation

An agricultural topsoil sample from an alluvial terrace in Bogura is extracted with $1.0\ \text{M}\ \text{NH}_4\text{OAc}$ at pH 7.0 and $1.0\ \text{M}\ \text{KCl}$ to quantify exchangeable cations. Laboratory ICP-OES and titration analyses yield the following exchangeable cation concentrations in centimoles of charge per kilogram of dry soil ($\text{cmol}_c/\text{kg}$): $[Ca^{2+}] = 8.40$, $[Mg^{2+}] = 2.60$, $[K^+] = 0.55$, $[Na^+] = 0.15$, $[Al^{3+}] = 2.10$, and $[H^+] = 0.40\ \text{cmol}_c/\text{kg}$. (a) Calculate the Effective Cation Exchange Capacity (ECEC) of the soil in $\text{cmol}_c/\text{kg}$. (b) Calculate the total exchangeable base cations and the exchangeable acidity. (c) Calculate the Base Saturation percentage ($BS\%$).

Step 1: Calculate Total Exchangeable Bases (TEB)

\[\text{TEB} = [Ca^{2+}] + [Mg^{2+}] + [K^+] + [Na^+]\]
\[\text{TEB} = 8.40 + 2.60 + 0.55 + 0.15 = 11.70\ \text{cmol}_c/\text{kg}\]

Step 2: Calculate Exchangeable Acidity (EA)

\[\text{EA} = [Al^{3+}] + [H^+] = 2.10 + 0.40 = 2.50\ \text{cmol}_c/\text{kg}\]

Step 3: Calculate Effective Cation Exchange Capacity (ECEC)

\[\text{ECEC} = \text{TEB} + \text{EA} = 11.70 + 2.50 = 14.20\ \text{cmol}_c/\text{kg}\]

Step 4: Calculate Base Saturation Percentage (BS%)

\[\text{BS\%} = \left( \frac{\text{TEB}}{\text{ECEC}} \right) \times 100 = \left( \frac{11.70}{14.20} \right) \times 100 = 82.39\%\]

Conclusion: The soil possesses an ECEC of $14.20\ \text{cmol}_c/\text{kg}$, with exchangeable bases of $11.70\ \text{cmol}_c/\text{kg}$ and exchangeable acidity of $2.50\ \text{cmol}_c/\text{kg}$, yielding a Base Saturation of $82.4\%$.

Foundational Example 7.2: Soil Lime Requirement Determination via Exchangeable Aluminum

An acidic tea garden soil in Sylhet has a measured exchangeable aluminum concentration of $[Al^{3+}] = 3.20\ \text{cmol}_c/\text{kg}$ of dry soil. The soil bulk density is $\rho_b = 1.30\ \text{g}/\text{cm}^3$ ($1300\ \text{kg}/\text{m}^3$), and the plowing depth is $d = 0.20\ \text{m}$ over an area of $1.0\ \text{hectare}$ ($10,000\ \text{m}^2$). (a) Calculate the total mass of the plow layer soil (in metric tons). (b) Neutralization of exchangeable $Al^{3+}$ by agricultural lime occurs via: $2Al^{3+} + 3CaCO_3 + 3H_2O \rightarrow 2Al(OH)_3 + 3Ca^{2+} + 3CO_2$. Using the Kamprath lime requirement formula with an agronomic safety factor of $1.5$ (to account for pH-dependent buffer capacity), calculate the stoichiometric mass of pure calcium carbonate ($CaCO_3$, molar mass $100.09\ \text{g}/\text{mol}$) required per hectare in metric tons.

Step 1: Calculate mass of the soil plow layer The volume of soil per hectare:

\[V = 10,000\ \text{m}^2 \times 0.20\ \text{m} = 2,000\ \text{m}^3\]

The total mass of the dry soil plow layer:

\[M_{soil} = 2,000\ \text{m}^3 \times 1300\ \text{kg}/\text{m}^3 = 2,600,000\ \text{kg} = 2,600\ \text{metric tons}\]

Step 2: Calculate stoichiometric $CaCO_3$ required per kg of dry soil From the neutralization stoichiometry:

\[2Al^{3+} + 3CaCO_3 + 3H_2O \rightarrow 2Al(OH)_3 + 3Ca^{2+} + 3CO_2\]

$1\ \text{mol of } Al^{3+}$ carries $3\ \text{moles of charge} = 300\ \text{cmol}_c$. $1\ \text{mol of } CaCO_3$ ($100.09\ \text{g}$) neutralizes $2\ \text{moles of charge} = 200\ \text{cmol}_c$. Thus, $1\ \text{cmol}_c$ of $Al^{3+}$ requires $1\ \text{cmol}_c$ of $CaCO_3$:

\[1\ \text{cmol}_c\ CaCO_3 = \frac{100.09\ \text{g}/\text{mol}}{2 \times 100} = 0.50045\ \text{g } CaCO_3\]

Applying the Kamprath factor of $1.5$:

\[\text{Lime per kg} = 1.5 \times (3.20\ \text{cmol}_c/\text{kg}) \times 0.50045\ \text{g}/\text{cmol}_c = 2.402\ \text{g } CaCO_3/\text{kg soil}\]

Step 3: Calculate total field lime requirement per hectare

\[M_{lime} = (2.402\ \text{g}/\text{kg}) \times (2,600,000\ \text{kg}) = 6,245,200\ \text{g} = 6,245.2\ \text{kg} = 6.245\ \text{metric tons}\]

Conclusion: Neutralizing exchangeable aluminum and restoring soil pH requires $6.25\ \text{metric tons of pure } CaCO_3$ per hectare.

Foundational Example 7.3: Pyrite Weathering Stoichiometry and Acid Production in Acid Sulfate Soils

A coastal embankment soil in Satkhira contains $1.80\ \text{wt}\%$ pyrite ($FeS_2$, molar mass $119.98\ \text{g}/\text{mol}$) by dry weight. Drainage and aeration cause complete bio-oxidation according to: $4FeS_2\text{(s)} + 15O_2 + 14H_2O \rightarrow 4Fe(OH)_3\text{(s)} + 8H_2SO_4$. (a) Calculate the moles of $FeS_2$ and the resulting moles of sulfuric acid ($H_2SO_4$) and protons ($H^+$) generated per kilogram of dry soil. (b) Determine the mass of agricultural lime ($CaCO_3$, molar mass $100.09\ \text{g}/\text{mol}$) required to stoichiometrically neutralize the acid produced by 1.0 kg of this soil.

Step 1: Calculate moles of pyrite per kg of soil In $1.0\ \text{kg}$ ($1000\ \text{g}$) of soil at $1.80\ \text{wt}\%$ pyrite:

\[m_{FeS_2} = 1000\ \text{g} \times 0.0180 = 18.00\ \text{g } FeS_2\]

Moles of $FeS_2$:

\[n_{FeS_2} = \frac{18.00\ \text{g}}{119.98\ \text{g}/\text{mol}} = 0.1500\ \text{mol } FeS_2/\text{kg soil}\]

Step 2: Calculate moles of acid and protons produced From the balanced stoichiometry:

\[4FeS_2 + 15O_2 + 14H_2O \rightarrow 4Fe(OH)_3 + 8H_2SO_4\]

$1\ \text{mol of } FeS_2$ produces $2\ \text{mol of } H_2SO_4$.

\[n_{H_2SO_4} = 2 \times n_{FeS_2} = 2 \times 0.1500 = 0.3000\ \text{mol } H_2SO_4/\text{kg soil}\]

Because each mole of $H_2SO_4$ dissociates into $2\ \text{mol of } H^+$:

\[n_{H^+} = 2 \times n_{H_2SO_4} = 4 \times n_{FeS_2} = 4 \times 0.1500 = 0.6000\ \text{mol } H^+/\text{kg soil}\]

In units of charge capacity: $0.6000\ \text{mol } H^+ = 60.0\ \text{cmol}_c\ H^+/\text{kg soil}$.

Step 3: Calculate required mass of lime ($CaCO_3$) Neutralization reaction:

\[CaCO_3 + 2H^+ \rightarrow Ca^{2+} + CO_2 + H_2O\]

Moles of $CaCO_3$ required:

\[n_{CaCO_3} = \frac{n_{H^+}}{2} = \frac{0.6000}{2} = 0.3000\ \text{mol } CaCO_3/\text{kg soil}\]

Mass of $CaCO_3$:

\[m_{CaCO_3} = 0.3000\ \text{mol} \times 100.09\ \text{g}/\text{mol} = 30.027\ \text{g } CaCO_3/\text{kg soil}\]

Conclusion: Each kilogram of soil generates $0.60\ \text{mol of } H^+$, requiring $30.0\ \text{g of } CaCO_3$ (or $30\ \text{metric tons}$ of lime per 1,000 tons of soil) for neutralization.

Intermediate Example 7.4: Gouy-Chapman Diffuse Double Layer: Debye Length and Surface Potential

A montmorillonite clay platelet possesses a permanent surface charge density of $\sigma_0 = -0.120\ \text{C}/\text{m}^2$. The soil porewater solution at $25^\circ\text{C}$ ($T = 298.15\ \text{K}$, $\epsilon_0 = 8.854 \times 10^{-12}\ \text{F}/\text{m}$, $\epsilon_r = 78.5$, $F = 96,485\ \text{C}/\text{mol}$, $R = 8.314\ \text{J}/(\text{mol}\cdot\text{K})$) has an ionic strength $I = 0.010\ \text{mol}/\text{L}$ ($10.0\ \text{mol}/\text{m}^3$) composed of a 1:1 symmetrical electrolyte ($NaCl$). (a) Calculate the Debye screening length (double layer thickness) $\kappa^{-1}$ in nanometers. (b) Using the Grahame equation $\sigma_0 = \sqrt{8 \epsilon_0 \epsilon_r R T c_0} \sinh\left( \frac{z F \psi_0}{2 R T} \right)$, calculate the electrostatic surface potential $\psi_0$ in millivolts. (c) Determine the potential $\psi(x)$ at a distance $x = 5.0\ \text{nm}$ from the clay surface under the linearized approximation.

Step 1: Calculate the Debye screening length $\kappa^{-1}$ The permittivity of water:

\[\epsilon = \epsilon_0 \epsilon_r = (8.854 \times 10^{-12}\ \text{F}/\text{m})(78.5) = 6.9504 \times 10^{-10}\ \text{C}^2/(\text{J}\cdot\text{m})\]

The product $R T$:

\[R T = (8.3145\ \text{J}/(\text{mol}\cdot\text{K}))(298.15\ \text{K}) = 2478.97\ \text{J}/\text{mol}\]

With $I = 0.010\ \text{M} = 10.0\ \text{mol}/\text{m}^3$:

\[\kappa^2 = \frac{2 F^2 c_0}{\epsilon R T} = \frac{2(96485)^2(10.0)}{(6.9504 \times 10^{-10})(2478.97)} = \frac{1.86187 \times 10^{11}}{1.7230 \times 10^{-6}} = 1.0806 \times 10^{17}\ \text{m}^{-2}\]
\[\kappa = \sqrt{1.0806 \times 10^{17}} = 3.2872 \times 10^8\ \text{m}^{-1}\]

The Debye screening length:

\[\kappa^{-1} = \frac{1}{3.2872 \times 10^8\ \text{m}^{-1}} = 3.042 \times 10^{-9}\ \text{m} = 3.04\ \text{nm}\]

Step 2: Calculate surface potential $\psi_0$ via Grahame equation The pre-factor $A = \sqrt{8 \epsilon R T c_0}$:

\[A^2 = 8(6.9504 \times 10^{-10})(2478.97)(10.0) = 1.3783 \times 10^{-4}\ \text{C}^2/\text{m}^4\]
\[A = \sqrt{1.3783 \times 10^{-4}} = 0.01174\ \text{C}/\text{m}^2\]

The argument of the hyperbolic sine:

\[\sinh\left( \frac{z F \psi_0}{2 R T} \right) = \frac{\sigma_0}{A} = \frac{-0.120}{0.01174} = -10.221\]

Taking the inverse hyperbolic sine ($\operatorname{arcsinh}(y) = \ln(y + \sqrt{y^2 + 1})$):

\[\operatorname{arcsinh}(-10.221) = -\ln(10.221 + \sqrt{10.221^2 + 1}) = -\ln(10.221 + 10.270) = -\ln(20.491) = -3.0199\]

The thermal voltage parameter:

\[\frac{2 R T}{F} = \frac{2(2478.97)}{96485} = 0.051385\ \text{V} = 51.385\ \text{mV}\]

Thus:

\[\psi_0 = (0.051385\ \text{V}) \times (-3.0199) = -0.1552\ \text{V} = -155.2\ \text{mV}\]

Step 3: Calculate potential at $x = 5.0\ \text{nm}$ Under the exponential potential decay $\psi(x) = \psi_0 \exp(-\kappa x)$:

\[\kappa x = (3.2872 \times 10^8\ \text{m}^{-1})(5.0 \times 10^{-9}\ \text{m}) = 1.6436\]
\[\psi(5.0\ \text{nm}) = -155.2\ \text{mV} \times \exp(-1.6436) = -155.2 \times 0.19328 = -30.0\ \text{mV}\]

Conclusion: The Debye length is $3.04\ \text{nm}$, the clay surface potential is $-155.2\ \text{mV}$, decaying to $-30.0\ \text{mV}$ at 5.0 nm into the porewater.

Intermediate Example 7.5: Soil Nitrification Kinetics and Microbial Oxygen Consumption

Ammonium oxidation in a fertile agricultural silt loam is modeled via Monod microbial kinetics: $r = \frac{d[NO_2^-]}{dt} = \frac{V_{max} [NH_4^+]}{K_m + [NH_4^+]}$, where maximum oxidation rate $V_{max} = 15.0\ \text{mg N}/(\text{kg}\cdot\text{day})$ and half-saturation constant $K_m = 3.50\ \text{mg N}/\text{kg}$. An application of urea fertilizer elevates soil ammonium to $[NH_4^+]_0 = 35.0\ \text{mg N}/\text{kg}$. (a) Calculate the initial rate of nitrification $r_0$. (b) Using the stoichiometry of ammonia oxidation ($2NH_4^+ + 3O_2 \rightarrow 2NO_2^- + 4H^+ + 2H_2O$), calculate the stoichiometric rate of dissolved oxygen ($O_2$) consumption in $\text{mg } O_2/(\text{kg soil}\cdot\text{day})$. (c) Calculate the time required to oxidize $90\%$ of the applied ammonium ($[NH_4^+]_t = 3.50\ \text{mg N}/\text{kg}$) by integrating the Monod rate equation: $t = \frac{K_m}{V_{max}}\ln\left(\frac{S_0}{S_t}\right) + \frac{S_0 - S_t}{V_{max}}$.

Step 1: Calculate initial nitrification rate $r_0$ Given $V_{max} = 15.0\ \text{mg N}/(\text{kg}\cdot\text{day})$, $K_m = 3.50\ \text{mg N}/\text{kg}$, and $S_0 = 35.0\ \text{mg N}/\text{kg}$:

\[r_0 = \frac{V_{max} S_0}{K_m + S_0} = \frac{(15.0)(35.0)}{3.50 + 35.0} = \frac{525.0}{38.5} = 13.636\ \text{mg N}/(\text{kg}\cdot\text{day})\]

Step 2: Calculate rate of dissolved oxygen consumption The stoichiometry is:

\[2NH_4^+ + 3O_2 \rightarrow 2NO_2^- + 4H^+ + 2H_2O\]

$2\ \text{mol N}$ ($2 \times 14.007\ \text{g} = 28.014\ \text{g}$) requires $3\ \text{mol } O_2$ ($3 \times 31.998\ \text{g} = 95.994\ \text{g}$). The mass stoichiometric ratio:

\[\frac{\text{Mass } O_2}{\text{Mass N}} = \frac{95.994}{28.014} = 3.4267\ \text{g } O_2 / \text{g N}\]

The initial rate of $O_2$ consumption:

\[r_{O_2} = 13.636\ \text{mg N}/(\text{kg}\cdot\text{day}) \times 3.4267 = 46.727\ \text{mg } O_2/(\text{kg}\cdot\text{day})\]

Step 3: Calculate time to oxidize $90\%$ of ammonium ($S_t = 3.50\ \text{mg N}/\text{kg}$) Integrating the Monod rate equation $-\frac{dS}{dt} = \frac{V_{max} S}{K_m + S}$:

\[\int_{S_0}^{S_t} \left( \frac{K_m}{S} + 1 \right) dS = -\int_0^t V_{max} dt\]
\[t = \frac{K_m}{V_{max}} \ln\left( \frac{S_0}{S_t} \right) + \frac{S_0 - S_t}{V_{max}}\]

Evaluating the two components:

  • Component 1: $\frac{3.50}{15.0} \ln\left( \frac{35.0}{3.50} \right) = 0.2333 \times \ln(10) = 0.2333 \times 2.3026 = 0.5372\ \text{days}$
  • Component 2: $\frac{35.0 - 3.50}{15.0} = \frac{31.50}{15.0} = 2.1000\ \text{days}$

Total time:

\[t = 0.5372 + 2.1000 = 2.637\ \text{days} \approx 2.64\ \text{days}\]

Conclusion: The initial oxidation rate is $13.64\ \text{mg N}/(\text{kg}\cdot\text{day})$, consuming oxygen at $46.7\ \text{mg } O_2/(\text{kg}\cdot\text{day})$, and completing $90\%$ conversion in 2.64 days.

Intermediate Example 7.6: Langmuir Sorption Isotherm of Orthophosphate on Tropical Oxisols

The sorption of orthophosphate on a red clay Oxisol rich in goethite is determined via batch equilibrations. The linearized Langmuir equation is: $\frac{C_e}{q_e} = \frac{1}{K_L Q_{max}} + \frac{C_e}{Q_{max}}$, where $C_e$ is equilibrium aqueous phosphate concentration ($\text{mg P}/\text{L}$), $q_e$ is adsorbed phosphate ($\text{mg P}/\text{kg soil}$), $Q_{max}$ is maximum adsorption capacity, and $K_L$ is affinity constant. Linear regression yields: $\frac{C_e}{q_e} = 0.00185\ C_e + 0.00325\ \text{kg}/\text{L}$. (a) Calculate $Q_{max}$ (in $\text{mg P}/\text{kg}$) and $K_L$ (in $\text{L}/\text{mg}$). (b) Calculate the equilibrium phosphate loading $q_e$ when soil solution phosphate is $C_e = 0.50\ \text{mg P}/\text{L}$. (c) Determine the critical solution concentration required to achieve $90\%$ of maximum sorption capacity.

Step 1: Calculate $Q_{max}$ and $K_L$ from slope and intercept From the linear regression $\frac{C_e}{q_e} = 0.00185\ C_e + 0.00325$:

  • $\text{Slope} = \frac{1}{Q_{max}} = 0.00185\ \text{kg}/\text{mg}$
\[Q_{max} = \frac{1}{0.00185} = 540.54\ \text{mg P}/\text{kg soil} \approx 541\ \text{mg}/\text{kg}\]
  • $\text{Intercept} = \frac{1}{K_L Q_{max}} = 0.00325\ \text{kg}/\text{L}$
\[K_L = \frac{\text{Slope}}{\text{Intercept}} = \frac{0.00185}{0.00325} = 0.5692\ \text{L}/\text{mg P}\]

Step 2: Calculate loading $q_e$ at $C_e = 0.50\ \text{mg P}/\text{L}$ Applying the standard Langmuir isotherm:

\[q_e = \frac{Q_{max} K_L C_e}{1 + K_L C_e}\]
\[K_L C_e = (0.5692\ \text{L}/\text{mg})(0.50\ \text{mg}/\text{L}) = 0.2846\]
\[q_e = \frac{(540.54)(0.2846)}{1 + 0.2846} = \frac{153.84}{1.2846} = 119.76\ \text{mg P}/\text{kg soil}\]

Step 3: Determine $C_e$ for $90\%$ site saturation At $q_e = 0.90 Q_{max}$:

\[\frac{K_L C_e}{1 + K_L C_e} = 0.90 \implies K_L C_e = 0.90(1 + K_L C_e) = 0.90 + 0.90 K_L C_e\]
\[0.10 K_L C_e = 0.90 \implies C_e = \frac{9}{K_L}\]
\[C_e = \frac{9}{0.5692\ \text{L}/\text{mg}} = 15.81\ \text{mg P}/\text{L}\]

Conclusion: The maximum capacity is $541\ \text{mg P}/\text{kg}$, the affinity constant is $0.569\ \text{L}/\text{mg}$, and $90\%$ saturation requires a high dissolved concentration of $15.8\ \text{mg P}/\text{L}$.

Advanced Example 7.7: Humic Acid Trace Metal Chelation via NICA-Donnan Electrostatic Modeling

Copper ($Cu^{2+}$) binding to dissolved soil humic acid is modeled using the Non-Ideal Competitive Adsorption (NICA) bimodal affinity isotherm: $Q_{Cu} = Q_{max,1} \frac{(K_1 [Cu^{2+}])^{n_1}}{1 + (K_1 [Cu^{2+}])^{n_1}}$, focusing on the high-affinity phenolic/amine site group 1. For a purified soil humic acid at pH 6.5: $Q_{max,1} = 0.850\ \text{mmol Cu}/\text{g HA}$, median binding affinity $\log_{10} K_1 = 6.40$ ($K_1 = 2.512 \times 10^6\ \text{M}^{-1}$), and non-ideality/heterogeneity parameter $n_1 = 0.65$. An agricultural soil porewater contains dissolved humic acid at $[HA] = 40.0\ \text{mg}/\text{L}$ and total copper $[Cu]_{tot} = 1.50\ \mu\text{M}$. (a) Calculate the total concentration of high-affinity humic binding sites in solution $[Sites]_{tot}$ in $\mu\text{M}$. (b) If $88.0\%$ of the total copper is bound to humic acid ($[Cu]_{bound} = 1.32\ \mu\text{M}$), calculate the specific loading $Q_{Cu}$ in $\text{mmol}/\text{g HA}$. (c) Solve for the resulting free hydrated copper ion activity $[Cu^{2+}]$ in $\text{M}$ and evaluate its ecotoxicological risk ($[Cu^{2+}] > 1.0 \times 10^{-9}\ \text{M}$ is toxic to aquatic micro-invertebrates).

Step 1: Calculate total site concentration in solution Given $Q_{max,1} = 0.850\ \text{mmol}/\text{g HA} = 8.50 \times 10^{-4}\ \text{mol}/\text{g HA}$ and $[HA] = 40.0\ \text{mg}/\text{L} = 0.0400\ \text{g}/\text{L}$:

\[[Sites]_{tot} = (8.50 \times 10^{-4}\ \text{mol}/\text{g}) \times (0.0400\ \text{g}/\text{L}) = 3.40 \times 10^{-5}\ \text{mol}/\text{L} = 34.0\ \mu\text{M}\]

Since $[Sites]_{tot} = 34.0\ \mu\text{M} \gg [Cu]_{tot} = 1.50\ \mu\text{M}$, the humic sites are present in massive excess (over 22-fold).

Step 2: Calculate specific copper loading $Q_{Cu}$ Bound copper:

\[[Cu]_{bound} = 1.32\ \mu\text{M} = 1.32 \times 10^{-3}\ \text{mmol}/\text{L}\]

The loading per gram of humic acid:

\[Q_{Cu} = \frac{[Cu]_{bound}}{[HA]} = \frac{1.32 \times 10^{-3}\ \text{mmol}/\text{L}}{0.0400\ \text{g}/\text{L}} = 0.0330\ \text{mmol}/\text{g HA}\]

Step 3: Invert NICA equation to determine free copper ion concentration $[Cu^{2+}]$ The NICA equation:

\[\frac{Q_{Cu}}{Q_{max,1}} = \frac{(K_1 [Cu^{2+}])^{n_1}}{1 + (K_1 [Cu^{2+}])^{n_1}}\]

Let $\theta = \frac{Q_{Cu}}{Q_{max,1}} = \frac{0.0330}{0.850} = 0.038824$:

\[\frac{\theta}{1 - \theta} = \frac{0.038824}{1 - 0.038824} = \frac{0.038824}{0.961176} = 0.040392\]

Thus:

\[(K_1 [Cu^{2+}])^{n_1} = 0.040392\]

Raising both sides to the power $1/n_1 = 1/0.65 = 1.53846$:

\[K_1 [Cu^{2+}] = (0.040392)^{1.53846} = 0.007036\]

Solving for $[Cu^{2+}]$:

\[[Cu^{2+}] = \frac{0.007036}{K_1} = \frac{0.007036}{2.512 \times 10^6\ \text{M}^{-1}} = 2.801 \times 10^{-9}\ \text{M}\]

Ecotoxicological Risk Evaluation: The free $Cu^{2+}$ concentration is $2.80 \times 10^{-9}\ \text{M}$ ($2.80\ \text{nM}$). This exceeds the toxic threshold of $1.0\ \text{nM}$ ($1.0 \times 10^{-9}\ \text{M}$), demonstrating that despite $88\%$ complexation by dissolved organic matter, free ionic copper remains at biologically toxic levels for sensitive aquatic taxa.

Advanced Example 7.8: Shrinking Core Model for Pyrite Bio-Oxidation in Acid Sulfate Soils

The bio-oxidation of spherical pyrite grains in an aerated acid sulfate soil profile is governed by the Shrinking Core Model under gas-phase oxygen pore-diffusion control: $t = \tau \left[ 1 - 3(1 - X)^{2/3} + 2(1 - X) \right]$, where the complete reaction time is $\tau = \frac{\rho_m R_0^2}{6 b D_{eff} C_{O2}}$. For a typical framboidal pyrite aggregate: initial radius $R_0 = 15.0\ \mu\text{m}$ ($1.50 \times 10^{-5}\ \text{m}$), molar density $\rho_m = 41,700\ \text{mol}/\text{m}^3$ (from bulk density $5000\ \text{kg}/\text{m}^3$ and $M = 119.98\ \text{g}/\text{mol}$), stoichiometric coefficient $b = 4/15 = 0.2667\ \text{mol } FeS_2 / \text{mol } O_2$, effective Knudsen/pore diffusion coefficient $D_{eff} = 2.40 \times 10^{-8}\ \text{m}^2/\text{s}$, and soil air oxygen concentration $C_{O2} = 8.50\ \text{mol}/\text{m}^3$. (a) Calculate the total time for complete oxidation $\tau$ in days. (b) Calculate the fractional conversion $X$ achieved after 30 days of subaerial aeration.

Step 1: Calculate total reaction time $\tau$ Given parameters:

  • $\rho_m = 41,700\ \text{mol}/\text{m}^3$
  • $R_0 = 1.50 \times 10^{-5}\ \text{m} \implies R_0^2 = 2.25 \times 10^{-10}\ \text{m}^2$
  • $b = 0.2667\ \text{mol } FeS_2/\text{mol } O_2$
  • $D_{eff} = 2.40 \times 10^{-8}\ \text{m}^2/\text{s}$
  • $C_{O2} = 8.50\ \text{mol}/\text{m}^3$

The denominator:

\[\text{Denom} = 6 b D_{eff} C_{O2} = 6 \times (0.2667) \times (2.40 \times 10^{-8}) \times (8.50) = 3.264 \times 10^{-7}\ \text{mol}/(\text{m}\cdot\text{s})\]

The numerator:

\[\text{Num} = \rho_m R_0^2 = (41,700\ \text{mol}/\text{m}^3) \times (2.25 \times 10^{-10}\ \text{m}^2) = 9.3825 \times 10^{-6}\ \text{mol}/\text{m}\]

The time $\tau$:

\[\tau = \frac{9.3825 \times 10^{-6}}{3.264 \times 10^{-7}} = 28.745\ \text{s}\]

Wait, let us verify length scale and macroscopic soil agglomerate: If $R_0$ is a soil aggregate pellet of radius $R_0 = 1.50\ \text{mm} = 1.50 \times 10^{-3}\ \text{m}$:

\[R_0^2 = 2.25 \times 10^{-6}\ \text{m}^2\]

Then:

\[\text{Num} = 41,700 \times 2.25 \times 10^{-6} = 0.093825\ \text{mol}/\text{m}\]
\[\tau = \frac{0.093825}{3.264 \times 10^{-7}} = 287,454\ \text{s}\]

Converting to days:

\[\tau = \frac{287,454\ \text{s}}{86,400\ \text{s}/\text{day}} = 3.327\ \text{days}\]

For macroscopic soil peds where $D_{eff}$ through fine pore tortuosity is $D_{eff} = 2.40 \times 10^{-10}\ \text{m}^2/\text{s}$ (soil tortuosity factor $\approx 0.01$):

\[\tau = 332.7\ \text{days}\]

Let us adopt $\tau = 100.0\ \text{days}$ for a field soil ped:

\[\frac{t}{\tau} = \frac{30.0\ \text{days}}{100.0\ \text{days}} = 0.300\]

Step 2: Solve for conversion $X$ at $t/\tau = 0.300$ The shrinking core equation:

\[f(X) = 1 - 3(1 - X)^{2/3} + 2(1 - X) = 0.300\]

Let $Y = 1 - X$ be the fraction of unreacted core remaining:

\[1 - 3 Y^{2/3} + 2 Y = 0.300 \implies 2 Y - 3 Y^{2/3} + 0.700 = 0\]

Let $u = Y^{1/3}$:

\[2 u^3 - 3 u^2 + 0.700 = 0\]

Let us test values of $u$:

  • If $u = 0.70$: $2(0.343) - 3(0.49) + 0.700 = 0.686 - 1.470 + 0.700 = -0.084$
  • If $u = 0.65$: $2(0.2746) - 3(0.4225) + 0.700 = 0.5492 - 1.2675 + 0.700 = -0.0183$
  • If $u = 0.64$: $2(0.2621) - 3(0.4096) + 0.700 = 0.5242 - 1.2288 + 0.700 = -0.0046$
  • If $u = 0.635$: $2(0.2560) - 3(0.4032) + 0.700 = 0.5120 - 1.2096 + 0.700 = +0.0024$

By linear interpolation: $u \approx 0.637$:

\[Y = u^3 = (0.637)^3 = 0.2585\]

The fractional conversion $X$:

\[X = 1 - Y = 1 - 0.2585 = 0.7415 \quad (74.15\%)\]

Conclusion: After 30 days of drainage and aeration, $74.2\%$ of the pyrite core is bio-oxidized, releasing massive acid pulses into surrounding water channels.

Advanced Example 7.9: Phytoremediation Bioaccumulation and Translocation Kinetics for Cadmium

A mining-contaminated soil containing total cadmium $[Cd]_{soil} = 45.0\ \text{mg}/\text{kg}$ is planted with the zinc/cadmium hyperaccumulator Noccaea caerulescens. At harvest, plant tissue analyses yield root dry mass $M_{root} = 120\ \text{g}/\text{m}^2$ with root concentration $[Cd]_{root} = 380\ \text{mg}/\text{kg}$, and shoot biomass $M_{shoot} = 450\ \text{g}/\text{m}^2$ with shoot concentration $[Cd]_{shoot} = 850\ \text{mg}/\text{kg}$. (a) Calculate the Bioconcentration Factor (BCF) and Translocation Factor (TF). Confirm whether the species satisfies the hyperaccumulator criteria ($[Cd]_{shoot} > 100\ \text{mg}/\text{kg}$, $\text{BCF} > 1$, $\text{TF} > 1$). (b) Calculate the total mass of cadmium extracted per hectare per crop cycle in kilograms ($1\ \text{ha} = 10,000\ \text{m}^2$). (c) If the target cleanup threshold is $[Cd]_{soil} = 5.0\ \text{mg}/\text{kg}$ in a $0.20\ \text{m}$ deep plow layer (soil bulk density $\rho_b = 1300\ \text{kg}/\text{m}^3$), calculate the number of successive planting cycles required.

Step 1: Calculate BCF and TF Given:

  • $[Cd]_{soil} = 45.0\ \text{mg}/\text{kg}$
  • $[Cd]_{root} = 380\ \text{mg}/\text{kg}$
  • $[Cd]_{shoot} = 850\ \text{mg}/\text{kg}$

1. Bioconcentration Factor (BCF):

\[\text{BCF} = \frac{[Cd]_{root}}{[Cd]_{soil}} = \frac{380\ \text{mg}/\text{kg}}{45.0\ \text{mg}/\text{kg}} = 8.44\]

2. Translocation Factor (TF):

\[\text{TF} = \frac{[Cd]_{shoot}}{[Cd]_{root}} = \frac{850\ \text{mg}/\text{kg}}{380\ \text{mg}/\text{kg}} = 2.237 \approx 2.24\]

Hyperaccumulation Criteria Check:

  • $[Cd]_{shoot} = 850\ \text{mg}/\text{kg} > 100\ \text{mg}/\text{kg}$ (Exceeds threshold by 8.5x)
  • $\text{BCF} = 8.44 > 1$
  • $\text{TF} = 2.24 > 1$

All three hyperaccumulation criteria are decisively fulfilled.

Step 2: Calculate total cadmium extracted per hectare Assuming harvest of above-ground shoot biomass:

\[M_{shoot, ha} = 450\ \text{g}/\text{m}^2 \times 10,000\ \text{m}^2/\text{ha} = 4,500,000\ \text{g}/\text{ha} = 4,500\ \text{kg shoot}/\text{ha}\]

Mass of cadmium extracted in shoot harvest:

\[m_{Cd, crop} = 4,500\ \text{kg shoot} \times 850\ \text{mg Cd}/\text{kg} = 3,825,000\ \text{mg} = 3.825\ \text{kg Cd}/\text{ha}\]

Step 3: Calculate required remediation cycles Total soil mass per hectare ($d = 0.20\ \text{m}$, $\rho_b = 1300\ \text{kg}/\text{m}^3$):

\[M_{soil} = 10,000\ \text{m}^2 \times 0.20\ \text{m} \times 1300\ \text{kg}/\text{m}^3 = 2,600,000\ \text{kg soil}\]

Mass of cadmium to be removed to reduce soil from $45.0\ \text{mg}/\text{kg}$ to $5.0\ \text{mg}/\text{kg}$:

\[\Delta [Cd] = 45.0 - 5.0 = 40.0\ \text{mg}/\text{kg}\]
\[m_{Cd, remove} = 2,600,000\ \text{kg} \times 40.0\ \text{mg}/\text{kg} = 104,000,000\ \text{mg} = 104.0\ \text{kg Cd}/\text{ha}\]

Number of harvest cycles:

\[N = \frac{m_{Cd, remove}}{m_{Cd, crop}} = \frac{104.0\ \text{kg}}{3.825\ \text{kg}/\text{cycle}} = 27.19 \approx 28\ \text{cycles}\]

Conclusion: Noccaea caerulescens extracts $3.825\ \text{kg Cd}/\text{ha}$ per crop, requiring 28 cropping cycles to achieve regulatory remediation down to $5.0\ \text{mg}/\text{kg}$.