GeoSysBRNS - A flexible multidimensional reactive transport model for simulating biogeochemical subsurface processes

The description of reactive transport processes in subsurface environments requires a sound understanding of both the biogeochemical complexity of the system and the spatially resolved transport of reactive species. However, most existing reactive transport models, for example in the field of contaminant hydrology, are specialized either in the simulation of the reactive or of the flow and transport processes. In this paper, we present and test the coupling of two highly flexible codes for the simulation of reactive transport processes in the subsurface: the Biogeochemical Reaction Network Simulator (BRNS), which contains a solver for kinetically and thermodynamically constrained biogeochemical reactions, and GeoSys/RockFlow, a multidimensional finite element subsurface flow and transport simulator. The new model, named GeoSysBRNS, maintains the full flexibility of the original models. The coupling is handled using an operator splitting scheme, which allows the reactive solver to be compiled into a problem specific library that is accessed by the transport simulator at runtime. The accuracy of the code coupling within GeoSysBRNS is demonstrated using two benchmark problems from the literature: a laboratory experiment on organic carbon degradation in a sand column via multiple microbial degradation pathways, and a dispersive mixing controlled bioreactive transport problem in aquifers, assuming three different reaction kinetics.

[1]  C. Steefel,et al.  Long-term fluxes of reactive species in macrotidal estuaries: Estimates from a fully transient, multicomponent reaction-transport model , 1997 .

[2]  Young-Jin Kim,et al.  A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions. , 2007, Journal of contaminant hydrology.

[3]  Wilfried Pfingsten,et al.  Efficient modeling of reactive transport phenomena by a multispecies random walk coupled to chemical equilibrium , 1996 .

[4]  Jean-François Gaillard,et al.  Biogeochemical dynamics in aquatic sediments , 1996 .

[5]  Robert A. Berner,et al.  Early Diagenesis: A Theoretical Approach , 1980 .

[6]  A. Chilakapati,et al.  RAFT: A simulator for ReActive Flow and Transport of groundwater contaminants , 1995 .

[7]  P. Regnier,et al.  A knowledge‐based reactive transport approach for the simulation of biogeochemical dynamics in Earth systems , 2005 .

[8]  D. A. Barry,et al.  MODFLOW/MT3DMS‐Based Reactive Multicomponent Transport Modeling , 2003, Ground water.

[9]  I M Young,et al.  Interactions and Self-Organization in the Soil-Microbe Complex , 2004, Science.

[10]  Sebastian Bauer,et al.  A process-oriented approach to computing multi-field problems in porous media , 2004 .

[11]  Michael J. Carr,et al.  Oxygen isotope constraints on the sources of Central American arc lavas , 2005 .

[12]  Olaf Kolditz,et al.  Assessing measurement uncertainty of first‐order degradation rates in heterogeneous aquifers , 2006 .

[13]  Olaf Kolditz,et al.  Object‐oriented finite element analysis of thermo‐hydro‐mechanical (THM) problems in porous media , 2007 .

[14]  Olaf Kolditz,et al.  Numerical analysis of coupled hydrosystems based on an object-oriented compartment approach , 2008 .

[15]  M. Chijimatsu,et al.  Results from an International Simulation Study on Coupled Thermal, Hydrological, and Mechanical Processes near Geological Nuclear Waste Repositories , 2008 .

[16]  Olaf Kolditz,et al.  Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions , 2008 .

[17]  Pierre Regnier,et al.  Quantitative interpretation of pH distributions in aquatic sediments: A reaction-transport modeling approach , 2005 .

[18]  Henning Prommer,et al.  Modelling the fate of oxidisable organic contaminants in groundwater. In C. T. Miller, M. B. Parlange, and S. M. Hassanizadeh (editors), , 2002 .

[19]  Olaf Kolditz,et al.  Results from an International Simulation Study on Coupled Thermal, Hydrological, and Mechanical Processes near Geological Nuclear Waste Repositories , 2008 .

[20]  Pierre Regnier,et al.  Modeling the impact of microbial activity on redox dynamics in porous media , 2005 .

[21]  Wolfgang Kinzelbach,et al.  Simulation of reactive processes related to biodegradation in aquifers. 2. Model application to a column study on organic carbon degradation , 1998 .

[22]  Wolfgang Kinzelbach,et al.  Simulation of reactive processes related to biodegradation in aquifers: 1. Structure of the three-dimensional reactive transport model , 1998 .

[23]  Timothy R. Ginn,et al.  Modeling microbial processes in porous media , 2000 .

[24]  Adam Brun,et al.  Modelling of transport and biogeochemical processes in pollution plumes: literature review and model development , 2002 .

[25]  B. Boudreau.  Diagenetic models and their implementation , 1997 .

[26]  Albert J. Valocchi,et al.  Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady state , 2007 .

[27]  Olaf Kolditz,et al.  Comments on “Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady-state” by O.A. Cirpka and A.J. Valocchi , 2009 .

[28]  David W. Blowes,et al.  Multicomponent reactive transport modeling in variably saturated porous media using a generalized formulation for kinetically controlled reactions , 2002 .

[29]  Urs von Gunten,et al.  Biogeochemical changes in groundwater-infiltration systems: Column studies , 1993 .

[30]  T. P. Clement,et al.  Modeling Multispecies Reactive Transport in Ground Water , 1998 .

[31]  Murugesu Sivapalan,et al.  Influence of Indian Ocean sea surface temperature variability on southwest Western Australian winter rainfall , 2006 .

[32]  Carlo H. R. Heip,et al.  Reactive transport in surface sediments I. Mexity and software quality , 2003 .

[33]  Albert J. Valocchi,et al.  Reply to comments on “Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady state” by H. Shao et al. , 2009 .

[34]  Pierre Regnier,et al.  Modeling Microbially Induced Carbon Degradation in Redox-Stratified Subsurface Environments: Concepts and Open Questions , 2007 .

[35]  Carl I. Steefel,et al.  Modeling complex multi-component reactive-transport systems: towards a simulation environment based on the concept of a Knowledge Base , 2002 .