Groundwater Arsenic Geochemistry, Oceanic Pollution & Ecotoxicology

The Bengal Basin represents the world's most catastrophic natural environmental catastrophe involving geogenic groundwater contamination. Understanding the distribution of arsenic across Bangladesh and West Bengal demands rigorous hydrostratigraphic and sedimentological contextualization.

Hydrostratigraphy and Geomorphic Architecture

The Bengal Basin comprises vast deltaic and alluvial floodplains deposited by the Ganges, Brahmaputra, and Meghna (GBM) river systems during the Late Pleistocene to Holocene epochs:

1. Holocene Alluvial Sediments (< 10,000 years BP):

• Typically grey, fine-to-medium sand, silt, and clay capping layers.
• Characterized by high concentrations of labile solid-phase organic matter (buried peat, plant debris, and marsh sediments) and low redox potential ($E_h < -50\ \text{mV}$).
• Host shallow aquifers (typically 10 to 50 meters below ground surface) displaying alarming aqueous arsenic levels exceeding $50\ \mu\text{g}/\text{L}$ (the national standard of Bangladesh) and frequently surpassing $500\ \mu\text{g}/\text{L}$, well above the World Health Organization guideline of $10\ \mu\text{g}/\text{L}$.

2. Pleistocene Sediments (Pleistocene Dupi Tila and Madhupur Formations):

• Comprise brown, red, or yellowish weathered sands, clays, and gravels.
• Deposited under oxidizing conditions or subjected to prolonged subaerial weathering and leaching during sea-level lowstands (e.g., Last Glacial Maximum ~20,000 BP).
• Iron minerals exist as crystalline goethite ($\alpha\text{-FeOOH}$) and hematite ($\alpha\text{-Fe}_2\text{O}_3$), which firmly sequester arsenic, resulting in arsenic-safe drinking water ($[\text{As}] < 10\ \mu\text{g}/\text{L}$).

Lithological Controls and Aquifer Confinement

The shallow alluvial aquifers are frequently capped by impermeable or semi-permeable surficial silty-clay aquitards (the 'upper clay layer'). This capping layer restricts the atmospheric diffusion of dissolved oxygen ($O_2$), generating closed-system anoxic hydrochemical environments that promote biogeochemical reducing reactions.

Arsenic is a toxic metalloid (Group 15, atomic number 33) occurring in natural aquatic systems predominantly in two oxidation states: trivalent arsenite [$\text{As(III)}$] and pentavalent arsenate [$\text{As(V)}$].

Inorganic Arsenate [$\text{As(V)}$] Speciation

Arsenate exists as arsenic acid ($H_3AsO_4$) and its conjugate base oxyanions, exhibiting three stepwise thermodynamic dissociation constants at $298.15\ \text{K}$:

In circumneutral groundwater (pH $6.5 - 8.5$), $\text{As(V)}$ exists primarily as the negatively charged divalent and monovalent oxyanions $H_2AsO_4^-$ and $HAsO_4^{2-}$. Because of their electrostatic charges, these species bind strongly to positively charged mineral surfaces via inner-sphere complexation.

Inorganic Arsenite [$\text{As(III)}$] Speciation

Arsenite exists as arsenous acid ($H_3AsO_3$), an uncharged, neutral pyramidal molecule under normal aquifer pH conditions:

Because $pK_{a1} = 9.22$, at typical groundwater pH ($6.5 - 8.0$), more than $95\%$ of $\text{As(III)}$ exists as the uncharged species $H_3AsO_3^0$. Neutral $H_3AsO_3^0$ does not experience significant Coulombic attraction toward charged mineral edges, rendering arsenite significantly more mobile and difficult to remove by conventional coagulants than arsenate.

Extensive multidisciplinary geochemical research has established that the microbially mediated reductive dissolution of iron(III) oxyhydroxides is the primary driver of high dissolved arsenic concentrations in Holocene aquifers of the Bengal Basin.

The Source of Geogenic Arsenic

Arsenic is not introduced via anthropogenic industrial discharge or pesticide application; rather, it originates in Himalayan metamorphic and ophiolitic headwaters. Fluvial erosion carries arsenic-bearing iron oxyhydroxide colloids downstream, depositing them as coatings on mineral sand grains in the delta plain.

Biogeochemical Reaction Mechanism

In buried Holocene sediments, buried reactive sedimentary organic matter ($CH_2O$) serves as the primary electron donor for indigenous dissimilatory metal-reducing bacteria (DMRB), notably species from the genera Geobacter (e.g., G. metallireducens), Shewanella (e.g., S. putrefaciens), and Desulfuromonas:

As the poorly crystalline hydrous ferric oxide (ferrihydrite, $5\text{Fe}_2\text{O}_3 \cdot 9\text{H}_2\text{O}$) dissolves into soluble ferrous iron ($\text{Fe}^{2+}$), the arsenic adsorbed to its mineral lattice is released into the interstitial porewater:

Role of Dissimilatory Arsenate-Reducing Microorganisms (DARM)

Simultaneously, specialized microorganisms equipped with periplasmic arsenate reductase enzymes ($ArrAB$) utilize $\text{As(V)}$ as a terminal respiratory electron acceptor:

This enzymatic reduction converts sorbed or dissolved arsenate into the far more mobile uncharged arsenite.

Arsenic mobility in aquifer matrices is heavily governed by surface complexation equilibria on mineral surfaces (iron, aluminum, and manganese oxides, and clay edges). Oxyanions naturally present in groundwater compete directly for available surface coordination sites.

Competitive Oxyanions: Phosphate ($PO_4^{3-}$), Silicate ($H_4SiO_4$), and Bicarbonate ($HCO_3^-$)

1. Phosphate ($PO_4^{3-}$):

• Possesses identical tetrahedral coordination geometry and similar thermochemical properties to arsenate ($AsO_4^{3-}$).
• Forms bidentate binuclear inner-sphere complexes ($\equiv\!\text{Fe}_2\text{O}_2\text{PO}_2$) with binding constants comparable to or exceeding those of arsenate:

• Elevated dissolved phosphate concentrations (arising from fertilizer runoff or organic matter remineralization) displace arsenate from sediment binding sites, driving arsenic into solution.

2. Silicate ($H_4SiO_4 / H_3SiO_4^-$):

• Present at high baseline concentrations in alluvial groundwater ($0.2 - 0.7\ \text{mmol}/\text{L}$) due to silicate mineral weathering.
• Competes moderately for binding sites, decreasing arsenic adsorption capacity across ferrihydrite and goethite surfaces.

3. Bicarbonate ($HCO_3^-$):

• Generated abundantly during organic matter oxidation (reaching concentrations of $300 - 800\ \text{mg}/\text{L}$).
• Promotes carbonate complexation and forms secondary iron carbonate minerals (siderite, $\text{FeCO}_3$), altering the surface charge and available specific surface area.

Thermodynamic stability fields of arsenic species in aquatic systems are quantified using electrochemical potential ($E_h$) and pH relationships, represented visually in Pourbaix diagrams.

Nernstian Formulation of the As(V)/As(III) Couple

The redox half-reaction connecting the predominant aqueous species at circumneutral pH is:

The standard reduction potential at standard conditions ($T = 298.15\ \text{K}$, $P = 1\ \text{bar}$, activities $a = 1$) is $E^\circ = +0.658\ \text{V}$.

Applying the Nernst equation:

Substituting constants $R = 8.314\ \text{J}/(\text{mol}\cdot\text{K})$, $T = 298.15\ \text{K}$, $F = 96485\ \text{C}/\text{mol}$, and $n = 2$:

At equal activities ($a_{H_3AsO_3} = a_{H_2AsO_4^-}$), the phase boundary simplifies to:

At $\text{pH} = 7.0$:

If the measured groundwater redox potential drops below $+37\ \text{mV}$, thermodynamically stable arsenic shifts entirely from $\text{As(V)}$ to mobile $\text{As(III)}$.

Geochemical Speciation Software (PHREEQC)

Modern hydrogeochemistry computes multi-component equilibrium using chemical speciation software (such as USGS PHREEQC), simultaneously solving mass balance, charge balance, activity coefficient corrections via the Davies or B-dot equations, and surface complexation models (Diffuse Double Layer or Generalized Two-Layer Models).

Engineering interventions to eliminate arsenic from drinking water supplies operate across household, community, and in-situ municipal scales.

Coagulation-Coprecipitation

The most widespread chemical treatment involves ferric chloride ($\text{FeCl}_3$) or alum ($\text{Al}_2(\text{SO}_4)_3 \cdot 18\text{H}_2\text{O}$) dosing:

1. Ferric ions rapidly hydrolyze to form an amorphous ferric hydroxide precipitate:

1. $\text{As(V)}$ oxyanions adsorb strongly onto the floc surface via ligand exchange:

1. Pre-oxidation Requirement: Because uncharged $\text{As(III)}$ binds poorly to ferric flocs, a pre-oxidation step (using hypochlorite, ozone, or solar Fenton processes) is mandatory to oxidize $\text{As(III)}$ to $\text{As(V)}$ prior to coagulation.

Fixed-Bed Adsorption Columns

Fixed-bed contactors utilize granular ferric hydroxide (GFH), activated alumina ($\text{Al}_2\text{O}_3$), or iron-coated sand. Breakthrough curves are governed by the bed volumes treated ($BV$) and the Thomas or Bohart-Adams kinetic transport models:

Subterranean In-situ Arsenic Remediation (SAR)

In-situ remediation avoids hazardous chemical sludge disposal at the surface. Oxygenated, aerated water is periodically injected down the tubewell into the aquifer surrounding the intake screen. Dissolved oxygen oxidizes naturally present $\text{Fe}^{2+}$ into insoluble $\text{Fe(OH)}_3$ coatings on sediment grains:

Subsequent groundwater extraction draws water through this newly formed reactive iron-mineral filter, sequestering arsenic in situ within the subterranean formation.

The oceans cover $71\%$ of the Earth's surface and represent the ultimate sink for chemical weathering fluxes. Marine chemistry is governed by conservative vs non-conservative chemical dynamics and density-driven stratification.

Chemical Composition of Seawater & Marcet's Principle

Seawater has an average practical salinity ($S$) of $35\ \text{g}/\text{kg}$ ($35\ \text{PSU}$). Marcet's Principle (the Law of Constant Proportions) establishes that the ratios of major conservative ions ($Cl^-$, $Na^+$, $SO_4^{2-}$, $Mg^{2+}$, $Ca^{2+}$, $K^+$) remain constant throughout the world's open oceans regardless of absolute salinity:

Ocean Water Column Stratification

The oceanic water column is physically partitioned into three distinct vertical zones:

1. Surface Mixed Layer (0 - 200 m): Wind-stirred, photic zone, saturated with $O_2$, depleted in macronutrients ($N, P, Si$) due to phytoplankton uptake.

2. Pycnocline (Thermocline & Halocline, 200 - 1000 m): Zone of rapid density increase, accompanied by a sharp temperature drop (thermocline) and salinity variation (halocline).

3. Deep Abyssal Layer (> 1000 m): Cold ($0 - 4^\circ\text{C}$), high-pressure, nutrient-rich water originating from polar deep-water formation (Thermohaline Circulation / Global Conveyor Belt).

Estuarine Mixing Dynamics

In estuaries, fresh river water ($S \approx 0$) mixes with saline seawater ($S \approx 35$). Chemical species are categorized by their behavior across the salinity gradient:

• Conservative Mixing: Linear relationship between concentration $C$ and salinity $S$:

• Non-conservative Mixing: Non-linear profiles exhibiting positive deviation (addition, due to desorption from riverine colloids) or negative deviation (removal, due to flocculation and sedimentation, common for humic acids, iron, and trace metals).

Ecotoxicology evaluates the transport, transformation, biological uptake, and toxicological outcomes of anthropogenic contaminants across marine and aquatic ecosystems.

Major Marine Contaminant Classes

1. Trace Heavy Metals ($Hg, Cd, Pb, Cu$):

• Methylmercury ($CH_3Hg^+$) represents the most bioavailable and neurotoxic marine organometal, forming stable covalent complexes with protein sulfhydryl ($-SH$) groups.

2. Petroleum Hydrocarbons and PAHs:

• Polycyclic aromatic hydrocarbons (benzo[a]pyrene, chrysene) induce DNA adduct formation and carcinogenesis via cytochrome P450 monooxygenase enzymatic bioactivation.

3. Marine Plastics and Microplastics (< 5 mm):

• Act as physical stressors and vector carriers that concentrate hydrophobic toxic chemicals (POPs, phthalates, bisphenol A) by several orders of magnitude relative to surrounding seawater.

Molecular and Cellular Biomarkers

Biomarkers provide early warning diagnostics of toxic exposure before organismal mortality occurs:

• Metallothioneins (MTs): Low molecular weight, cysteine-rich proteins synthesized in liver/hepatopancreas tissues that chelate toxic divalent metal cations ($Cd^{2+}, Hg^{2+}, Pb^{2+}$).
• Ethoxyresorufin-O-deethylase (EROD) Induction: Measures the catalytic activity of the CYP1A enzyme family, serving as a specific biomarker for exposure to coplanar PCBs, dioxins, and PAHs.
• Acetylcholinesterase (AChE) Inhibition: Biomarker indicating organophosphate and carbamate neurotoxin contamination.