Using dual‐domain advective‐transport simulation to reconcile multiple‐tracer ages and estimate dual‐porosity transport parameters

• Published in: Water resources research 2017-06, Vol.53 (6), p.5002-5016
• Main Authors: Sanford, Ward E., Niel Plummer, L., Casile, Gerolamo, Busenberg, Ed, Nelms, David L., Schlosser, Peter
• Format: Article
• Language: English
• Subjects: Advection; Age; Algorithms; Atmospheric models; Carbonate rocks; Carbonates; Chemicals; Computer simulation; Correlation; Dissolved chemicals; dual domain; environmental tracer; Environmental tracers; fractured rock; Fractures; Groundwater; Groundwater age; Karst; Lines; Mathematical models; Modelling; Movement; Parameter estimation; Parameter sensitivity; Parameters; Pollutants; Pollution dispersion; Porosity; River basins; Rivers; Rocks; Simulation; Solute transport; Solutes; Spring; Tracers; Transport; Travel; Void space; Water springs; Watersheds
• ISSN: 0043-1397; 1944-7973
• DOI: 10.1002/2016WR019469

Dual‐domain transport is an alternative conceptual and mathematical paradigm to advection‐dispersion for describing the movement of dissolved constituents in groundwater. Here we test the use of a dual‐domain algorithm combined with advective pathline tracking to help reconcile environmental tracer concentrations measured in springs within the Shenandoah Valley, USA. The approach also allows for the estimation of the three dual‐domain parameters: mobile porosity, immobile porosity, and a domain exchange rate constant. Concentrations of CFC‐113, SF6, 3H, and 3He were measured at 28 springs emanating from carbonate rocks. The different tracers give three different mean composite piston‐flow ages for all the springs that vary from 5 to 18 years. Here we compare four algorithms that interpret the tracer concentrations in terms of groundwater age: piston flow, old‐fraction mixing, advective‐flow path modeling, and dual‐domain modeling. Whereas the second two algorithms made slight improvements over piston flow at reconciling the disparate piston‐flow age estimates, the dual‐domain algorithm gave a very marked improvement. Optimal values for the three transport parameters were also obtained, although the immobile porosity value was not well constrained. Parameter correlation and sensitivities were calculated to help quantify the uncertainty. Although some correlation exists between the three parameters being estimated, a watershed simulation of a pollutant breakthrough to a local stream illustrates that the estimated transport parameters can still substantially help to constrain and predict the nature and timing of solute transport. The combined use of multiple environmental tracers with this dual‐domain approach could be applicable in a wide variety of fractured‐rock settings. Plain Language Summary The movement of dissolved chemicals, such as pollutants, underground has been notoriously difficult to predict in terrains composed of fractured rock because of the complex way in which the chemicals mix between fractures and the host rock. One evidence for this has been that when different atmospheric‐derived tracers are measured in springs in karst terrains they give different apparent times of travel between the land surface and the spring. The authors of this paper have presented a new approach that can both partially characterize the nature of how such chemicals are mixing underground and reconcile the difference in apparent times of travel. The authors took a