Watershed-scale monitoring, field aeration experiments, and geochemical equilibrium and kinetic modeling were conducted to evaluate interdependent changes in pH, dissolved CO₂, O₂, 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 CO₂ 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 = k₁·[O₂]·[H⁺]⁻²·[Fe(II)] that indicates an increase in pH by 1 unit at pH 5–8 and at constant dissolved O₂ (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 (k₁′) of 0.5–1.7 times the reference value of k₁ = 3 × 10⁻¹² mol/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 CO₂ outgassing and O₂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 CO₂ mass transfer rate constant (k_L,_CO2a). The value of k_L,_CO2a estimated for the stream (0.010 min⁻¹) was within the range for the batch aeration experiments (0–0.033 min⁻¹). These results indicate that the abiotic homogeneous Fe(II) oxidation rate law, with adjustments for variations in temperature and CO₂ 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).