Monitoring, field experiments, and geochemical modeling of Fe(II) oxidation kinetics in a stream dominated by net-alkaline coal-mine drainage, Pennsylvania, USA

Abstract

Watershed-scale monitoring, field aeration experiments, and geochemical equilibrium and kinetic modeling were conducted to evaluate interdependent changes in pH, dissolved CO2, O2, and Fe(II) concentrations that typically take place downstream of net-alkaline, circumneutral coal-mine drainage (CMD) outfalls and during aerobic treatment of such CMD. The kinetic modeling approach, using PHREEQC, accurately simulates observed variations in pH, Fe(II) oxidation, alkalinity consumption, and associated dissolved gas concentrations during transport downstream of the CMD outfalls (natural attenuation) and during 6-h batch aeration tests on the CMD using bubble diffusers (enhanced attenuation). The batch aeration experiments demonstrated that aeration promoted CO2 outgassing, thereby increasing pH and the rate of Fe(II) oxidation. The rate of Fe(II) oxidation was accurately estimated by the abiotic homogeneous oxidation rate law −d[Fe(II)]/dt=k1·[O2]·[H+]−2·[Fe(II)] that indicates an increase in pH by 1 unit at pH 5–8 and at constant dissolved O2 (DO) concentration results in a 100-fold increase in the rate of Fe(II) oxidation. Adjusting for sample temperature, a narrow range of values for the apparent homogeneous Fe(II) oxidation rate constant (k1′) of 0.5–1.7 times the reference value of k1=3×10−12mol/L/min (for pH 5–8 and 20°C), reported by Stumm and Morgan (1996), was indicated by the calibrated models for the 5-km stream reach below the CMD outfalls and the aerated CMD. The rates of CO2 outgassing and O2 ingassing in the model were estimated with first-order asymptotic functions, whereby the driving force is the gradient of the dissolved gas concentration relative to equilibrium with the ambient atmosphere. Although the progressive increase in DO concentration to saturation could be accurately modeled as a kinetic function for the conditions evaluated, the simulation of DO as an instantaneous equilibrium process did not affect the model results for Fe(II) or pH. In contrast, the model results for pH and Fe(II) were sensitive to the CO2 mass transfer rate constant (kL,CO2a). The value of kL,CO2a estimated for the stream (0.010min−1) was within the range for the batch aeration experiments (0–0.033min−1). These results indicate that the abiotic homogeneous Fe(II) oxidation rate law, with adjustments for variations in temperature and CO2 outgassing rate, may be applied to predict changes in aqueous iron and pH for net-alkaline, ferruginous waters within a stream (natural conditions) or a CMD treatment system (engineered conditions).