Numerical modeling of an abiotic hyporheic mixing-dependent reaction: Chemical evolution of mixing and reactant production zones.

Mixing-dependent reactions occur where groundwater and surface water mix in shallow sediments (hyporheic zone) and can attenuate contaminants along upwelling flowpaths, thus reducing transport to surface water. Here we used MODFLOW/SEAM3D to numerically simulate prior laboratory observations of a mixing-dependent reaction between sodium sulfite (NaSO) and dissolved oxygen (DO) to produce sodium sulfate (NaSO). This reaction is not common in nature but is used as a surrogate for mixing-dependent DO consuming reactions of environmental significance. We evaluated how location and thickness of mixing zones and reaction product production zones dynamically respond to variations in hydraulic and chemical boundary conditions and reaction kinetic rate. Sensitivity analysis showed that location and thickness of mixing zones and reactant production zones were most sensitive to changes in the balance of hydrologic inflow from groundwater and surface water (inflow ratio). Mixing zone thickness for reactive DO calibrated to experimental data was thinner than that for the "DO tracer" (identical source location and concentration as DO but conservative tracer), indicating that as DO is consumed its mixing zone narrows. The SO production zone was consistently thicker than the DO mixing zone. Small changes in mixing/production zone thicknesses were linked to large changes in mass consumed and produced, indicating the potential for simpler field metrics like thickness to act as surrogates for more challenging measurements such as contaminant flux or consumption in monitoring natural attenuation. This study improves understanding of the evolution of hyporheic mixing-dependent reaction zones that occur even under steady state hydraulics, emphasizing their complex controls.