Two-dimensional transport model for variably saturated porous media with major ion chemistry

We present the development and demonstrate the use of the two-dimensional finite element code UNSATCHEM-2D for modeling major ion equilibrium and kinetic nonequilibrium chemistry in variably saturated porous media. The model is intended for prediction of major ion chemistry and water and solute fluxes for soils under transient conditions. Since the solution chemistry in the unsaturated zone is significantly influenced by variations in water content, temperature, and CO2 concentrations in the soil gas, all these variables are also calculated by the model. The major variables of the chemical system are Ca, Mg, Na, K, SO4, Cl, NO3, alkalinity, and CO2. The model accounts for equilibrium chemical reactions between these components such as complexation, cation exchange, and precipitation-dissolution. For the precipitation-dissolution of calcite and dissolution of dolomite, either equilibrium or multicomponent kinetic expressions are used which include both forward and back reactions. Other dissolution-precipitation reactions considered include gypsum, hydromagnesite, and nesquehonite. Since the ionic strength of soil solutions can often reach high values, both modified Debye-Huckel and Pitzer expressions were incorporated into the model to calculate single ion activities. The need for an iterative coupling procedure between the solute transport and chemical modules is demonstrated with an example which considers root water uptake and irrigation using moderately saline water. The utility of the model is further illustrated with two-dimensional simulations with surface and subsurface irrigation from a line source.

[1]  van Bavel,et al.  a Soil Aeration Theory Based on Diffusion , 1951 .

[2]  J. Quirk,et al.  Permeability of porous solids , 1961 .

[3]  E. Jong,et al.  CALCULATION OF SOIL RESPIRATION AND ACTIVITY FROM CO2 PROFILES IN THE SOIL , 1972 .

[4]  B. Jones,et al.  WATEQ: A COMPUTER PROGRAM FOR CALCULATING CHEMICAL EQUILIBRIA OF NATURAL WATERS , 1973 .

[5]  E. Bresler,et al.  Simultaneous transport of solutes and water under transient unsaturated flow conditions , 1973 .

[6]  Kenneth S. Pitzer,et al.  Thermodynamics of electrolytes. I. Theoretical basis and general equations , 1973 .

[7]  S. P. Neuman,et al.  Finite element simulation of flow in saturated - unsaturated soils considering water uptake by plants , 1974 .

[8]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[9]  D. Hillel,et al.  Simulation of Profile Water Storage as Related to Soil Hydraulic Properties1 , 1976 .

[10]  W. G. Gray,et al.  Finite Element Simulation in Surface and Subsurface Hydrology , 1977 .

[11]  D. Suarez.  Ion Activity Products of Calcium Carbonate in Waters Below the Root Zone1 , 1977 .

[12]  J. D. Rhoades,et al.  Reducing water quality degradation through minimized leaching management , 1977 .

[13]  R. Feddes,et al.  Simulation of field water use and crop yield , 1978 .

[14]  David L. Parkhurst,et al.  The kinetics of calcite dissolution in CO 2 -water systems at 5 degrees to 60 degrees C and 0.0 to 1.0 atm CO 2 , 1978 .

[15]  A. Lasaga,et al.  The treatment of multi-component diffusion and ion pairs in diagenetic fluxes , 1979 .

[16]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[17]  C. W. Robbins,et al.  A Combined Salt Transport-Chemical Equilibrium Model for Calcareous and Gypsiferous Soils , 1980 .

[18]  C. W. Robbins,et al.  Calculating cation exchange in a salt transport model. , 1980 .

[19]  Robert L. Street,et al.  Transport of ion‐exchanging solutes in groundwater: Chromatographic theory and field simulation , 1981 .

[20]  A. Gureghian.  A two-dimensional finite-element solution for the simultaneous transport of water and multi solutes through a nonhomogeneous aquifer under transient saturated-unsaturated flow conditions , 1981 .

[21]  Description of soil chemistry during transient solute transport , 1981 .

[22]  E. Busenberg,et al.  The kinetics of dissolution of dolomite in CO 2 -H 2 O systems at 1.5 to 65 degrees C and O to 1 atm PCO 2 , 1982 .

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

[24]  J. Rhoades,et al.  The Apparent Solubility of Calcium Carbonate in Soils1 , 1982 .

[25]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

[26]  N. Møller,et al.  The prediction of mineral solubilities in natural waters: The Na-K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O system to high ionic strengths at 25°C , 1984 .

[27]  Larry W. Lake,et al.  Precipitation and dissolution of solids attending flow through porous media , 1984 .

[28]  T. N. Narasimhan,et al.  Groundwater contamination from an inactive uranium mill tailings pile: 2. Application of a dynamic mixing model , 1984 .

[29]  Aaron A. Jennings,et al.  Multisolute mass transport with chemical interaction kinetics , 1985 .

[30]  D. S. Ward,et al.  Finite Element Matrix and Mass Balance Computational Schemes for Transport in Variably Saturated Porous Media , 1985 .

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

[32]  G. Yeh,et al.  FEMA: a Finite Element Model of Material Transport through Aquifers , 1985 .

[33]  Paul R. Bloom,et al.  An evaluation of rate equations for calcite precipitation kinetics at pCO2 less than 0.01 atm and pH greater than 8 , 1985 .

[34]  R. S. Schechter,et al.  Interactions of precipitation/dissolution waves and ion exchange in flow through permeable media , 1986 .

[35]  David Russo,et al.  Simulation of Leaching of A Gypsiferous‐Sodic Desert Soil , 1986 .

[36]  John H. Weare,et al.  The prediction of borate mineral equilibria in natural waters: Application to Searles Lake, California , 1986 .

[37]  W. Inskeep,et al.  Kinetics of Calcite Precipitation in the Presence of Water-soluble Organic Ligands1 , 1986 .

[38]  F. Morel,et al.  MINEQL-EIR. Chemical Equilibrium Composition of Aqueous Systems , 1986 .

[39]  R. Förster.  A multicomponent transport model , 1986 .

[40]  Linda M. Abriola,et al.  Modeling Contaminant Transport in the Subsurface: An Interdisciplinary Challenge (Paper 7R0031) , 1987 .

[41]  D. Solomon,et al.  The annual carbon dioxide cycle in a montane soil: Observations, modeling, and implications for weathering , 1987 .

[42]  John L. Nieber,et al.  Oxygen, Carbon Dioxide, and Water Transfer in Soils: Mechanisms and Crop Response , 1988 .

[43]  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 .

[44]  Tomas Vogel.  SWMII-Numerical model of two-dimensional flow in a variably saturated porous medium. , 1988 .

[45]  P.S. Huyakorn,et al.  VAM2D: Variably saturated analysis model in two dimensions , 1989 .

[46]  T. N. Narasimhan,et al.  Redox-Controlled Multiple-Species Reactive Chemical Transport 2. Verification and Application , 1989 .

[47]  T. N. Narasimhan,et al.  Redox‐controlled multiple‐species reactive chemical transport: 1. Model development , 1989 .

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

[49]  V. S. Tripathi,et al.  HYDROGEOCHEM: A coupled model of HYDROlogic transport and GEOCHEMical equilibria in reactive multicomponent systems , 1990 .

[50]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[51]  Jacob Rubin.  Solute transport with multisegment, equilibrium‐controlled reactions: A feed forward simulation method , 1990 .

[52]  Gour-Tsyh Yeh,et al.  A Model for Simulating Transport of Reactive Multispecies Components: Model Development and Demonstration , 1991 .

[53]  J. T. McCord.  Application of Second‐Type Boundaries in Unsaturated Flow Modeling , 1991 .

[54]  Chin-Fu Tsang,et al.  A summary of subsurface hydrological and hydrochemical models , 1991 .

[55]  Donald L. Suarez,et al.  Modeling of carbon dioxide transport and production in soil: 2. Parameter selection, sensitivity analysis, and comparison of model predictions to field data , 1993 .

[56]  J. Simunek,et al.  UNSATCHEM-2D code for simulating two-dimensional variably saturated water flow, heat transport, carbon dioxide production and transport, and multicomponent solute transport with major ion equilibrium and kinetic chemistry , 1993 .

[57]  Donald L. Suarez,et al.  Modeling of Carbon Dioxide Transport and Production in Soil 1. Model Development , 1993 .