One-Dimensional Transport with Equilibrium Chemistry (OTEQ): A Reactive Transport Model for Streams and Rivers Completed

Within OTEQ, 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 (waterborne) 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. The model's ability to simulate pH, precipitation/dissolution, and pH-dependent sorption provides a means of evaluating the complex interactions between instream chemistry and hydrologic transport at the field scale.

OTEQ is generally applicable to solutes which undergo reactions that are sufficiently fast relative to hydrologic processes ("Local Equilibrium"). Although the definition of "sufficiently fast" is highly solute and application dependent, many reactions involving inorganic solutes quickly reach a state of chemical equilibrium. Given a state of chemical equilibrium, inorganic solutes may be modeled using OTEQ's equilibrium approach. This equilibrium approach is facilitated through the use of an existing database that describes chemical equilibria for a wide range of inorganic solutes. In addition, solute reactions not included in the existing database may be added by defining the appropriate mass-action equations and the associated equilibrium constants. As such, OTEQ provides a general framework for the modeling of solutes under the assumption of chemical equilibrium. Despite this generality, most OTEQ applications to date have focused on the transport of metals in streams and small rivers. The OTEQ documentation is therefore focused on metal transport. Potential model users should note, however, that additional applications are possible.