Date of this Version
WATER RESOURCES RESEARCH, VOL. 34, NO. 10, PAGES 2501-2514, OCTOBER 1998
Multisolute reactive transport was investigated in chemically heterogeneous systems to determine the influence of the shapes of the reactive heterogeneities (iron oxide inclusions) by comparison of two-dimensional heterogeneous experiments with spatially averaged models that had differing inclusion characterization. Eleven reactions were considered in this system, starting with absorption of the initial solute (CoIIEDTA) to Fe oxides, followed by two competing surface reactions: oxidation forming CoIIIEDTA and Fe dissolution forming FeIIIEDTA and Co2+. Spatial moments of the eight mobile species were compared between data and models. One spatially averaged model (homogenous equivalent), which incorporated inclusion mass only, significantly under predicted oxidation (up to 74%), the influence of reaction kinetics, and species retardation. In contrast, the ensemble average model (incorporating inclusion mass and length) well predicted speciation, retardation, and skewness. This large difference in prediction between two spatially averaging models was caused by the lack of incorporation of contact time of solutes with iron oxides in the homogenous equivalent model and the importance of the contact time with the differing timescales of reactions.