Steven C. Chapra | Geology Reports

USGS1996

Reactive solute transport in streams: 1. Development of an equilibrium- based model

An equilibrium-based solute transport model is developed for the simulation of trace metal fate and transport in streams. The model is formed by coupling a solute transport model with a chemical equilibrium submodel based on MINTEQ. The solute transport model considers the physical processes of advection, dispersion, lateral inflow, and transient storage, while the equilibrium submodel considers the speciation and complexation of aqueous species, precipitation/dissolution and sorption. Within the model, reactions in the water column may result in the formation of solid phases (precipitates and sorbed species) that are subject to downstream transport and settling processes. Solid phases on the streambed may also interact with the water column through dissolution and sorption/desorption reactions. Consideration of both mobile (water-borne) and immobile (streambed) solid phases requires a unique set of governing differential equations and solution techniques that are developed herein. The partial differential equations describing physical transport and the algebraic equations describing chemical equilibria are coupled using the sequential iteration approach.

Water Resources Research↗

USGS1996

Reactive solute transport in streams: 2. Simulation of a pH modification experiment

We present an application of an equilibrium-based solute transport model to a p H-modification experiment conducted on the Snake River, an acidic, metal-rich stream located in the Rocky Mountains of Colorado. During the experiment, instream p H decreased from 4.2 to 3.2, causing a marked increase in dissolved iron concentrations. Model application requires specification of several parameters that are estimated using tracer techniques, mass balance calculations, and geochemical data. Two basic questions are addressed through model application: (1) What are the processes responsible for the observed increase in dissolved iron concentrations? (2) Can the identified processes be represented within the equilibrium-based transport model? Simulation results indicate that the increase in iron was due to the dissolution of hydrous iron oxides and the photoreduction of ferric iron. Dissolution from the streambed is represented by considering a trace compartment consisting of freshly precipitated hydrous iron oxide and an abundant compartment consisting of aged precipitates that are less soluble. Spatial variability in the solubility of hydrous iron oxide is attributed to heterogeneity in the streambed sediments, temperature effects, and/or variability in the effects of photoreduction. Solubility products estimated via simulation fall within a narrow range ( pK sp from 40.2 to 40.8) relative to the 6 order of magnitude variation reported for laboratory experiments ( pK sp from 37.3 to 43.3). Results also support the use of an equilibrium-based transport model as the predominate features of the iron and p H profiles are reproduced. The model provides a valuable tool for quantifying the nature and extent of p H-dependent processes within the context of hydrologic transport.

Water Resources Research↗