Multicomponent reactive transport modeling in variably saturated porous media using a generalized formulation for kinetically controlled reactions

[1] A generalized formulation for kinetically controlled reactions has been developed and incorporated into a multicomponent reactive transport model to facilitate the investigation of a large variety of problems involving inorganic and organic chemicals in variably saturated media. The general kinetic formulation includes intra-aqueous and dissolution-precipitation reactions in addition to geochemical equilibrium expressions for hydrolysis, aqueous complexation, oxidation-reduction, ion exchange, surface complexation, and gas dissolution-exsolution reactions. The generalized approach allows consideration of fractional order terms with respect to any dissolved species in terms of species activities or in terms of total concentrations, which facilitates the incorporation of a variety of experimentally derived rate expressions. Monod and inhibition terms can be used to describe microbially mediated reactions or to limit the reaction progress of inorganic reactions. Dissolution-precipitation reactions can be described as surface-controlled or transport-controlled reactions. The formulation also facilitates the consideration of any number of parallel reaction pathways, and reactions can be treated as irreversible or reversible processes. Two groundwater contamination scenarios, both set in variably saturated media but with significantly different geochemical reaction networks, are investigated and demonstrate the advantage of the generalized approach. The first problem focuses on a hypothetical case study of the natural attenuation of organic contaminants undergoing dissolution, volatilization, and biodegradation in an unconfined aquifer overlaid by unsaturated sediments. The second problem addresses the generation of acid mine drainage in the unsaturated zone of a tailings impoundment at the Nickel Rim Mine Site near Sudbury, Ontario, and subsequent reactive transport in the saturated portion of the tailings.

[1]  James E. Szecsody,et al.  Groundwater flow, multicomponent transport and biogeochemistry: development and application of a coupled process model , 2000 .

[2]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[3]  Peter C. Lichtner,et al.  Continuum formulation of multicomponent-multiphase reactive transport , 1996 .

[4]  E A Sudicky,et al.  Multicomponent simulation of wastewater-derived nitrogen and carbon in shallow unconfined aquifers. I. Model formulation and performance. , 2001, Journal of contaminant hydrology.

[5]  J. Bahr,et al.  Direct comparison of kinetic and local equilibrium formulations for solute transport affected by surface reactions , 1987 .

[6]  Keith Loague,et al.  A compartmentalized solute transport model for redox zones in contaminated aquifers: 2. Field‐scale simulations , 2000 .

[7]  E. Frind,et al.  Modelling of sulfide oxidation with reactive transport at a mine drainage site , 2000 .

[8]  F. Schwartz,et al.  Simulating the in situ oxidative treatment of chlorinated ethylenes by potassium permanganate , 2000 .

[9]  G. Heron,et al.  Anaerobic microbial redox processes in a landfill leachate contaminated aquifer (Grindsted, Denmark) , 1998 .

[10]  R. Jakobsen,et al.  Redox zoning, rates of sulfate reduction and interactions with Fe-reduction and methanogenesis in a shallow sandy aquifer, Rømø, Denmark , 1999 .

[11]  J. Rubin,et al.  Transport of reacting solutes in porous media: Relation between mathematical nature of problem formulation and chemical nature of reactions , 1983 .

[12]  Lei Chou,et al.  Steady-state kinetics and dissolution mechanisms of albite , 1985 .

[13]  S. Brantley,et al.  Chemical weathering rates of pyroxenes and amphiboles , 1995 .

[14]  R. J. Millington,et al.  Gas Diffusion in Porous Media , 1959, Science.

[15]  P. Huyakorn,et al.  Techniques for Making Finite Elements Competitve in Modeling Flow in Variably Saturated Porous Media , 1984 .

[16]  Stephen Jacob Kowall.  The DOE Vadose Zone Science And Technology Roadmap: A National Program to Address Characterization, Monitoring and Simulation of Subsurface Contaminant Fate and Transport , 2001 .

[17]  Bruce E. Rittmann,et al.  Natural attenuation strategy for groundwater cleanup focuses on demonstrating cause and effect , 2001 .

[18]  A. Lasaga.  Kinetic theory in the earth sciences , 1998 .

[19]  S. Brantley,et al.  Chapter 4. CHEMICAL WEATHERING RATES OF PYROXENES AND AMPHIBOLES , 1995 .

[20]  D. Blowes,et al.  Microbial populations associated with the generation and treatment of acid mine drainage , 2000 .

[21]  LorettaY Li,et al.  Evaluation of vadose zone biodegradation of BTX vapours , 2000 .

[22]  D. L. Parkhurst,et al.  User's guide to PHREEQC (Version 2)-a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 1999 .

[23]  J. Schnoor.  Kinetics of chemical weathering : a comparison of laboratory and field weathering rates , 1990 .

[24]  Hans-Jørgen Albrechtsen,et al.  Characterization of redox conditions in groundwater contaminant plumes , 2000 .

[25]  P. Lichtner,et al.  Surface reaction versus diffusion control of mineral dissolution and growth rates in geochemical processes , 1989 .

[26]  Philippe Van Cappellen,et al.  Kinetic modeling of microbially-driven redox chemistry of subsurface environments : coupling transport, microbial metabolism and geochemistry , 1998 .

[27]  J. Allison,et al.  MINTEQA2/PRODEFA2, a geochemical assessment model for environmental systems: Version 3. 0 user's manual , 1991 .

[28]  K. Loague,et al.  A compartmentalized solute transport model for redox zones in contaminated aquifers: 1. Theory and development , 2000 .

[29]  Robert W. Gillham,et al.  Pyrite oxidation in carbonate-buffered solution: 2. Rate control by oxide coatings , 1990 .

[30]  R. Berner.  Rate control of mineral dissolution under Earth surface conditions , 1978 .

[31]  René Therrien,et al.  Improved three-dimensional finite-element techniques for field simulation of variably saturated flow and transport , 1993 .

[32]  J. Gould.  The kinetics of hexavalent chromium reduction by metallic iron , 1982 .

[33]  M. V. Genuchten,et al.  Using Texture and Other Soil Properties to Predict the Unsaturated Soil Hydraulic Functions , 1988 .

[34]  R. Knapp.  Spatial and temporal scales of local equilibrium in dynamic fluid-rock systems , 1989 .

[35]  R. Baker,et al.  Estimation of hydrocarbon biodegradation rates in gasoline-contaminated sediment from measured respiration rates , 2000 .

[36]  A. Blum,et al.  Feldspar dissolution kinetics , 1995 .

[37]  E. Sudicky,et al.  Simulation of biodegradable organic contaminants in groundwater. 1. Numerical formulation in principal directions. , 1990 .

[38]  Donald I. Siegel,et al.  Crude oil in a shallow sand and gravel aquifer-III , 1993 .

[39]  D. W. Pollock,et al.  Gas transport in unsaturated porous media: The adequacy of Fick's law , 1989 .

[40]  D. Blowes,et al.  The hydrogeochemistry of the Nickel Rim mine tailings impoundment, Sudbury, Ontario , 2000 .

[41]  Robert L. Street,et al.  A Groundwater Mass Transport and Equilibrium Chemistry Model for Multicomponent Systems , 1985 .

[42]  Kathryn L. Nagy,et al.  Dissolution and precipitation kinetics of sheet silicates , 1995 .

[43]  Larry W. Lake,et al.  Effect of partial local equilibrium on the propagation of precipitation/dissolution waves , 1993 .

[44]  David W. Blowes,et al.  Reactive transport modeling of an in situ reactive barrier for the treatment of hexavalent chromium and trichloroethylene in groundwater , 2001 .

[45]  V. S. Tripathi,et al.  A critical evaluation of recent developments in hydrogeochemical transport models of reactive multichemical components , 1989 .

[46]  K. Pruess,et al.  MODELING MULTIPHASE NON-ISOTHERMAL FLUID FLOW AND REACTIVE GEOCHEMICAL TRANSPORT IN VARIABLY SATURATED FRACTURED ROCKS: 1. METHODOLOGY , 2001 .

[47]  Donald L. Suarez,et al.  Two-dimensional transport model for variably saturated porous media with major ion chemistry , 1994 .

[48]  P. Singer,et al.  Acidic Mine Drainage: The Rate-Determining Step , 1970, Science.

[49]  R. Garrels,et al.  Comparative study of the kinetics and mechanisms of dissolution of carbonate minerals , 1989 .

[50]  Gregory B. Davis,et al.  A model of oxidation in pyritic mine wastes: part 1 equations and approximate solution , 1986 .

[51]  A. Blum,et al.  Chapter 7. FELDSPAR DISSOLUTION KINETICS , 1995 .

[52]  Carl I. Steefel,et al.  Chapter 2. APPROACHES TO MODELING OF REACTIVE TRANSPORT IN POROUS MEDIA , 1996 .

[53]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[54]  C. Miller,et al.  Simulation of Solute Transport in a Chemically Reactive Heterogeneous System' Model Development and Application , 1983 .

[55]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[56]  R. Schwarzenbach,et al.  Environmental Organic Chemistry , 1993 .

[57]  H. Sverdrup,et al.  Weathering of primary silicate minerals in the natural soil environment in relation to a chemical weathering model , 1988, Water, Air, and Soil Pollution.

[58]  Jeanne M. VanBriesen,et al.  Multicomponent transport with coupled geochemical and microbiological reactions: model description and example simulations , 1998 .

[59]  A. Bourg,et al.  Speciation of Fe(II) and Fe(III) in Contaminated Aquifer Sediments Using Chemical Extraction Techniques. , 1994, Environmental science & technology.

[60]  C. Steefel,et al.  Approaches to modeling of reactive transport in porous media , 1996 .

[61]  A. Jennings,et al.  Multicomponent equilibrium chemistry in groundwater quality models , 1982 .

[62]  S. P. Neuman,et al.  SATURATED-UNSATURATED SEEPAGE BY FINITE ELEMENTS , 1973 .

[63]  Peter C. Lichtner,et al.  Continuum model for simultaneous chemical reactions and mass transport in hydrothermal systems , 1985 .

[64]  C. Appelo,et al.  Geochemistry, groundwater and pollution , 1993 .

[65]  Howard W. Reeves,et al.  Multicomponent mass transport with homogeneous and heterogeneous chemical reactions: Effect of the chemistry on the choice of numerical algorithm: 1. Theory , 1988 .

[66]  E. Frind,et al.  Sulfide mineral oxidation and subsequent reactive transport of oxidation products in mine tailings impoundments: A numerical model , 1996 .

[67]  A. Walter,et al.  Modeling of multicomponent reactive transport in groundwater-2 , 1994 .

[68]  K. L. Nagy,et al.  Chapter 5. DISSOLUTION AND PRECIPITATION KINETICS OF SHEET SILICATES , 1995 .

[69]  S F Thornton,et al.  Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer. , 2001, Journal of contaminant hydrology.

[70]  E. Frind,et al.  Simulation of contaminant transport in three dimensions: 2. Dimensionality effects , 1987 .

[71]  P. Bedient,et al.  Transport of dissolved hydrocarbons influenced by oxygen‐limited biodegradation: 1. Theoretical development , 1986 .

[72]  E. Frind,et al.  Advective-dispersive transport of dense organic vapors in the unsaturated zone 1 , 1990 .

[73]  E. Frind,et al.  Advective‐dispersive transport of dense organic vapors in the unsaturated zone: 2. Sensitivity analysis , 1990 .

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

[75]  Peter A. Forsyth,et al.  Mechanisms Controlling Vacuum Extraction Coupled With Air Sparging for Remediation of Heterogeneous Formations Contaminated by Dense Nonaqueous Phase Liquids , 1995 .

[76]  R Kahnt,et al.  Modelling the closure-related geochemical evolution of groundwater at a former uranium mine. , 2001, Journal of contaminant hydrology.

[77]  O. Levenspiel.  Chemical Reaction Engineering , 1972 .

[78]  C. Steefel,et al.  A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution rea , 1994 .

[79]  P. C. Lichtner.  Scaling properties of the time-space kinetic mass transport equations and the local equilibrium limit , 1993 .