Aquatic Chemistry, Carbonate Equilibria & Water Quality Metrics

Water ($\text{H}_2\text{O}$) is the universal solvent of aquatic ecosystems, endowed with anomalous thermodynamic and transport properties governed by its condensed-phase hydrogen-bonding network.

1. Molecular Geometry and Electric Dipole Moment

The isolated gas-phase water molecule possesses $C_{2v}$ point group symmetry with an experimental $\text{H}-\text{O}-\text{H}$ bond angle of $104.5^\circ$ and an $\text{O}-\text{H}$ bond length of $0.0958\text{ nm}$. Electronegativity differences between oxygen ($\chi = 3.44$) and hydrogen ($\chi = 2.20$) induce substantial bond polarization, producing an isolated molecular permanent electric dipole moment:

In condensed liquid water, cooperative mutual polarization enhances the effective molecular dipole moment to:

Liquid water possesses an exceptionally large static dielectric constant (relative permittivity $\varepsilon_r \approx 78.4$ at $298\text{ K}$), dramatically reducing the electrostatic Coulombic interaction energy between dissolved ionic species:

This intense dielectric screening facilitates the spontaneous solvation and electrolytic dissociation of ionic lattices.

2. Tetrahedral Hydrogen-Bonding Topology and Anomalies

Each oxygen atom acts as a dual hydrogen-bond donor (via its two polarized covalent $\text{O}-\text{H}$ bonds) and a dual hydrogen-bond acceptor (via its two $sp^3$-hybridized non-bonding lone pairs). This forms a fluctuating, transient three-dimensional tetrahedral network characterized by:

• An average hydrogen bond energy of $\Delta H_{\text{HB}} \approx 21\text{ kJ/mol}$.
• Density Maximum at $3.98^\circ\text{C}$ ($277.13\text{ K}$): Above $3.98^\circ\text{C}$, normal thermal kinetic expansion dominates; below $3.98^\circ\text{C}$, the progressive formation of an open, low-density tetrahedral cage network (pre-ice $I_h$ structuring) decreases liquid density. This ensures that freshwater lakes freeze from the top down, preserving benthic aquatic ecosystems.
• High specific heat capacity ($c_p = 4.184\text{ J/(g}\cdot\text{K)}$), high enthalpy of vaporization ($\Delta H_{\text{vap}} = 40.66\text{ kJ/mol}$), and high surface tension ($\gamma = 72.8\text{ mN/m}$ at $20^\circ\text{C}$).

The aquatic carbonate system is the principal acid-base buffer regulating natural water bodies.

1. Fundamental Carbonate Equilibria

Dissolved inorganic carbon comprises four interconverting chemical species:

1. Gaseous carbon dioxide dissolution:

where $[\text{CO}_2^*] = [\text{CO}_2(aq)] + [\text{H}_2\text{CO}_3]$.

1. First dissociation:

1. Second dissociation:

The total concentration of dissolved inorganic carbon is:

Defining alpha ionization fractions:

2. Alkalinity and Acid Neutralizing Capacity (ANC)

Total Alkalinity (Alk) represents the equivalent sum of all titratable proton-accepting bases measured relative to a carbonic acid equivalence point ($ ext{pH} \approx 4.5$):

In typical natural waters with $6.0 < \text{pH} < 9.0$:

Alkalinity is a conservative quantity under pressure and temperature changes as well as gas exchange ($ ext{CO}_2$ addition or loss changes $ ext{pH}$ and $C_T$, but leaves Total Alkalinity strictly invariant!).

3. Calcite Saturation and The Langelier Saturation Index

The saturation state with respect to calcite ($ ext{CaCO}_3$) is:

The Langelier Saturation Index (LSI) is:

where $\text{pH}_s$ is the saturation pH:

• $\text{LSI} > 0$: Supersaturated; tendency for scale deposition.
• $\text{LSI} < 0$: Undersaturated; corrosive toward distribution pipes.

Dissolved oxygen (DO) is the master parameter governing aquatic aerobic respiration, biogeochemical redox potential, and the survival of fish and benthic invertebrates.

1. Gas Dissolution Thermodynamics and Henry's Law

The dissolution of molecular oxygen in water is an exothermic process ($\Delta H_{\text{diss}} < 0$):

According to Henry's law:

where $P_{\text{O}_2} = y_{\text{O}_2} (P - P_{\text{H}_2\text{O}})$ is the dry partial pressure of atmospheric oxygen ($y_{\text{O}_2} = 0.2095$). The temperature dependence of Henry's law coefficient is described by the van 't Hoff equation:

Because dissolution is exothermic ($\Delta H_{\text{diss}} \approx -11.7\text{ kJ/mol}$), Henry's solubility coefficient decreases with increasing temperature. At sea level ($1.0\text{ atm}$):

• At $T = 0^\circ\text{C}$ ($273\text{ K}$): $[ ext{O}_2]_{\text{sat}} \approx 14.62\text{ mg/L}$
• At $T = 20^\circ\text{C}$ ($293\text{ K}$): $[ ext{O}_2]_{\text{sat}} \approx 9.09\text{ mg/L}$
• At $T = 35^\circ\text{C}$ ($308\text{ K}$): $[ ext{O}_2]_{\text{sat}} \approx 6.95\text{ mg/L}$

2. Salting-Out Effect in Saline Waters

Dissolved electrolytes reduce the solubility of non-polar gases through electrostrictive hydration of ions (the Setchenow salting-out equation):

where $k_s$ is the Setchenow salting-out coefficient and $I_s$ is the ionic strength of the water. At $20^\circ\text{C}$, full-strength seawater ($S = 35\text{ g/kg}$) holds only $\approx 7.35\text{ mg/L}$ of DO at saturation, compared to $9.09\text{ mg/L}$ in freshwater.

3. Hypoxia and Critical Ecological Thresholds

Aquatic ecology defines three dissolved oxygen regimes:

• Normoxia ($> 6.0\text{ mg/L}$): Healthy conditions supporting diverse teleost fish and macroinvertebrate communities.
• Stress Zone ($2.0 - 5.0\text{ mg/L}$): Physiological distress, reduced swimming speeds, avoidance behavior.
• Hypoxia ($< 2.0\text{ mg/L}$): Acute mortality for fish; benthic fauna mortality.
• Anoxia ($[\text{O}_2] \equiv 0.0\text{ mg/L}$): Complete depletion, initiating anaerobic microbial reduction of $\text{NO}_3^-, \text{Fe}^{3+}, \text{SO}_4^{2-}$, producing toxic hydrogen sulfide ($\text{H}_2\text{S}$).

Biochemical Oxygen Demand (BOD) quantifies the amount of dissolved oxygen consumed by heterotrophic microorganisms during the biological oxidation of organic matter.

1. Biological Oxidation Stoichiometry

The aerobic degradation of carbonaceous organic matter (represented generically as glucose $\text{C}_6\text{H}_{12}\text{O}_6$ or cell biomass $\text{C}_5\text{H}_7\text{NO}_2$) proceeds via:

One mole of glucose ($180.16\text{ g}$) requires $6\text{ moles}$ of $\text{O}_2$ ($192.00\text{ g}$), yielding a theoretical carbonaceous demand of $1.066\text{ g O}_2\text{/g glucose}$.

2. First-Order Kinetics of Carbonaceous BOD Exertion

The rate of deoxygenation is assumed to be proportional to the concentration of remaining biodegradable organic matter $L(t)$ (measured in oxygen equivalents, $ ext{mg/L}$):

Integrating yields the remaining unoxidized substrate:

where $L_0$ is the Ultimate Carbonaceous BOD ($ ext{CBOD}_u$) at $t = 0$, and $k_1$ is the first-order deoxygenation rate constant (base $e$, typically $0.10 - 0.35\text{ day}^{-1}$ at $20^\circ\text{C}$). The cumulative oxygen consumed—the BOD exerted $y(t)$—is:

Using common base 10 rate constants ($K_{10} = k_1 / \ln(10) = k_1 / 2.3026$):

3. The Standard 5-Day BOD Test (BOD5)

Because total oxidation can require 20 to 30 days, regulatory agencies standardize testing at $t = 5\text{ days}$ incubated in darkness at $20^\circ\text{C}$ ($ ext{BOD}_5$):

For a typical municipal wastewater with $k_1 = 0.23\text{ day}^{-1}$:

Thus, $\text{BOD}_5$ captures approximately $68\%$ of the ultimate carbonaceous demand.

4. Nitrogenous BOD (NBOD)

After $\approx 5 - 8\text{ days}$, autotrophic nitrifying bacteria (Nitrosomonas and Nitrobacter) proliferate, oxidizing un-ionized ammonia and ammonium:

Stoichiometrically, oxidizing $1\text{ mole}$ of nitrogen ($14.01\text{ g N}$) consumes $2\text{ moles}$ of $\text{O}_2$ ($64.00\text{ g O}_2$):

In uninhibited BOD bottles, the onset of nitrification produces a secondary surge in oxygen exertion.

While BOD measures only biologically accessible organics over days, Chemical Oxygen Demand (COD) and Total Organic Carbon (TOC) provide rapid, comprehensive chemical metrics.

1. Potassium Dichromate Reflux Chemistry

Chemical Oxygen Demand (COD) quantifies the equivalent oxygen capacity required to chemically oxidize virtually all organic compounds to carbon dioxide and water using boiling hexavalent potassium dichromate ($\text{K}_2\text{Cr}_2\text{O}_7$) in concentrated sulfuric acid:

For a generic organic molecule $\text{C}_n\text{H}_a\text{O}_b\text{N}_c$:

(Nitrogen in amino groups is converted to ammonia, rather than oxidized to nitrate). The theoretical COD per mole of compound is:

2. Analytical Catalysis and Halide Interference

• Silver Sulfate Catalyst ($ ext{Ag}_2 ext{SO}_4$): Added ($~10\text{ g/L}$) to catalyze the oxidation of refractory straight-chain aliphatic hydrocarbons and fatty acids, which otherwise resist dichromate attack.
• Mercuric Sulfate Masking ($ ext{HgSO}_4$): Chloride ions ($\text{Cl}^-$) are readily oxidized by dichromate, creating false-positive COD:

Adding $\text{HgSO}_4$ forms an un-ionized mercuric chloride complex:

masking up to $2000\text{ mg/L}$ of chloride. The unreacted dichromate is back-titrated with standard ferrous ammonium sulfate (FAS, $\text{Fe}(\text{NH}_4)_2(\text{SO}_4)_2$) using ferroin indicator:

3. Total Organic Carbon (TOC) and Biodegradability Index

• TOC Measurement: Sample is acidified to purge Inorganic Carbon ($\text{IC} = \text{CO}_2, \text{HCO}_3^-$), and non-purgeable organic carbon is oxidized catalytically at $680-950^\circ\text{C}$ to $\text{CO}_2$, quantified by a non-dispersive infrared (NDIR) detector.
• Biodegradability Ratio ($\text{BOD}_5 / \text{COD}$):
• $\text{BOD}_5 / \text{COD} > 0.5$: Readily biodegradable (municipal domestic sewage).
• $0.2 < \text{BOD}_5 / \text{COD} < 0.5$: Moderately biodegradable.
• $\text{BOD}_5 / \text{COD} < 0.2$: Refractory, recalcitrant, or toxic (industrial chemical or landfill leachate).

The Streeter-Phelps model (1925) mathematically couples the biochemical consumption of oxygen by biodegradable organic matter with atmospheric reaeration across the river surface.

1. Coupled Differential Equations of the Oxygen Sag

Consider a point-source wastewater discharge into a 1D plug-flow stream of cross-sectional area $A$ and average velocity $u$. Let $x$ be downstream distance and $t = x/u$ be the downstream travel time. Let $D(t) = [ ext{O}_2]_{\text{sat}} - [ ext{O}_2](t)$ be the dissolved oxygen deficit ($ ext{mg/L}$). Two competing first-order rate processes govern the system:

1. Deoxygenation: Oxygen consumed during microbial oxidation of organic matter $L(t)$:

2. Atmospheric Reaeration: Oxygen flux across the gas-liquid interface driven by the deficit:

where $k_2$ is the reaeration rate constant (governed by stream depth $H$ and velocity $u$ via the O'Connor-Dobbins equation: $k_2 = 3.93 u^{0.5} / H^{1.5}$).

The net rate of change of the oxygen deficit is:

2. Analytical Solution of the Classical Streeter-Phelps Equation

This is a first-order linear ordinary differential equation. Multiplying by the integrating factor $\exp(k_2 t)$:

Integrating from $t = 0$ (where $D(0) = D_0$) to $t$:

Dividing through by $\exp(k_2 t)$ yields the classical Streeter-Phelps equation:

3. The Critical Deficit (Dc) and Critical Travel Time (tc)

The point of maximum oxygen depletion occurs where $dD/dt = 0$, meaning deoxygenation exactly balances reaeration ($k_1 L = k_2 D_c$):

Solving for the critical travel time $t_c$:

The critical distance downstream is $x_c = u \times t_c$. The critical maximum oxygen deficit is:

The minimum dissolved oxygen in the river is:

If $[ ext{DO}]_{\text{min}} < 2.0\text{ mg/L}$, a septic, hypoxic fish-kill zone develops.

Water hardness reflects the total concentration of multivalent metallic cations in solution, principally calcium ($ ext{Ca}^{2+}$) and magnesium ($ ext{Mg}^{2+}$).

1. Classification and Geochemical Origin

Hardness originates from the dissolution of limestone ($ ext{CaCO}_3$) and dolomite ($ ext{CaMg}( ext{CO}_3)_2$) by soil carbonic acid:

Hardness is conventionally categorized as:

• Carbonate (Temporary) Hardness: Associated with bicarbonate and carbonate anions; can be precipitated out by boiling:

• Non-Carbonate (Permanent) Hardness: Associated with sulfate ($ ext{SO}_4^{2-}$), chloride ($ ext{Cl}^-$), and nitrate ($ ext{NO}_3^-$) anions; unaffected by boiling.

Total hardness is expressed quantitatively as equivalent concentration of calcium carbonate in $\text{mg/L as CaCO}_3$:

Classification standards:

• Soft: $0 - 60\text{ mg/L as CaCO}_3$
• Moderately Hard: $61 - 120\text{ mg/L}$
• Hard: $121 - 180\text{ mg/L}$
• Very Hard: $> 180\text{ mg/L}$

2. Complexometric Chelation Titration with EDTA

Disodium ethylenediaminetetraacetate ($ ext{Na}_2 ext{H}_2 ext{Y}$) forms exceptionally stable $1:1$ octahedral chelate complexes with divalent metal ions:

Titration is buffered at $\text{pH} = 10.0$ using an $\text{NH}_3/\text{NH}_4\text{Cl}$ buffer. The metallochromic indicator Eriochrome Black T (EBT, $\text{H}_2\text{In}^-$) binds to free magnesium:

Because the EDTA chelate $[ ext{CaY}]^{2-}$ is more stable than $[ ext{MgY}]^{2-}$, EDTA titrates all free $ ext{Ca}^{2+}$ first, then free $ ext{Mg}^{2+}$, and finally displaces EBT from $[ ext{MgIn}]^-$, restoring the brilliant blue color of free $ ext{HIn}^{2-}$ at the sharp endpoint.

The chemical speciation of multivalent elements in water is controlled simultaneously by proton activity ($ ext{pH}$) and electron activity ($ ext{p}arepsilon$).

1. Definition and Thermodynamic Derivation of pe

Analogous to $\text{pH} = -\log_{10} a_{\text{H}^+}$, the dimensionless electron activity $\text{p}\varepsilon$ is defined as:

Consider a generic reduction half-reaction:

The equilibrium constant is:

Taking the logarithm:

Solving for $\text{p}\varepsilon$:

where $\text{p}\varepsilon^\circ = \frac{1}{n} \log_{10} K$. Relating $\text{p}\varepsilon$ to the reduction potential $E_h$ (in Volts):

2. Limits of Water Stability

In natural waters, $\text{p}\varepsilon$ and $\text{pH}$ are bounded by the thermodynamic stability limits of liquid water:

1. Upper Oxidizing Limit (Oxygen Evolution):

2. Lower Reducing Limit (Hydrogen Evolution):

3. The Thermodynamic Redox Ladder

Microorganisms extract metabolic Gibbs free energy by coupling organic carbon oxidation to terminal electron acceptors in order of decreasing $\text{p}\varepsilon$:

1. Aerobic Respiration: $\text{O}_2 / \text{H}_2\text{O}$ ($ ext{p}\varepsilon^\circ(W) = +13.75$ at $ ext{pH}=7$)

2. Denitrification: $\text{NO}_3^- / \text{N}_2$ ($ ext{p}\varepsilon^\circ(W) = +12.65$)

3. Manganese Reduction: $\text{MnO}_2 / \text{Mn}^{2+}$ ($ ext{p}\varepsilon^\circ(W) = +8.9$)

4. Iron Reduction: $\text{FeOOH} / \text{Fe}^{2+}$ ($ ext{p}\varepsilon^\circ(W) = -0.8$)

5. Sulfate Reduction: $\text{SO}_4^{2-} / \text{H}_2\text{S}$ ($ ext{p}\varepsilon^\circ(W) = -3.75$)

6. Methanogenesis: $\text{CO}_2 / \text{CH}_4$ ($ ext{p}\varepsilon^\circ(W) = -4.13$)