Inverse optimization of hydraulic, solute transport, and cation exchange parameters using HP1 and UCODE to simulate cation exchange.

Reactive transport modeling is a powerful tool to evaluate systems with complex geochemical relations. However, parameters are not always directly measurable. This study represents one of the first attempts to obtain hydrologic, transport and geochemical parameters from an experimental dataset involving transient unsaturated water flow and solute transport, using an automatic inverse optimization (or calibration) algorithm. The data come from previously published, controlled laboratory experiments on the transport of major cations (Na, K, Mg, Ca) during water absorption into horizontal soil columns that were terminated at different times. Experimental data consisted of the depth profiles of water contents (θ), Cl concentrations, and total aqueous and sorbed concentrations of major cations. The dataset was used to optimize several parameters using the reactive transport model, HP1 and the generic optimization code, UCODE. Although the soil hydraulic and solute transport parameters were also optimized, the study focused mainly on the geochemical parameters because the soil columns were constructed from disturbed soil. The cation exchange capacity and the cation exchange coefficients for two exchange models (Gapon and Rothmund-Kornfeld) were optimized. The results suggest that both calibrated models satisfactorily described the experimental data, although the Rothmund-Kornfeld model fit was slightly better. However, information content and surface response analyses indicated that parameters of the Gapon model are well identifiable, whereas those of the Rothmund-Kornfeld model were strongly correlated. The calibrated geochemical parameters were validated using an independent dataset. In agreement with the identifiability analysis, the Gapon approach was better than the Rothmund-Kornfeld model at calculating the observed concentrations of major cations in the soil solution and on the exchange sites.

[1]  C. J. Smith,et al.  Absorption of artificial piggery effluent by soil: A laboratory study , 2004 .

[2]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[3]  A. N. D. R E A S V O E G E L I N, † V I J A,et al.  Reaction-Based Model Describing Competitive Sorption and Transport of Cd , Zn , and Ni in an Acidic Soil , 2022 .

[4]  R. E. Phillips,et al.  Cation Exchange in Soils , 2018, SSSA Book Series.

[5]  M. Friedel.  Coupled inverse modeling of vadose zone water, heat, and solute transport: calibration constraints, parameter nonuniqueness, and predictive uncertainty , 2005 .

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

[7]  J. Vrugt,et al.  Inverse Modeling of Subsurface Flow and Transport Properties: A Review with New Developments , 2008 .

[8]  W. Bond,et al.  Simulating cation transport during water flow in soil: two approaches , 1993 .

[9]  Florian Kuhnen,et al.  Multicomponent transport of major cations predicted from binary adsorption experiments , 2000 .

[10]  C. Marshall.  Cation Exchange in Soils.Walter P. Kelley , 1950 .

[11]  Diederik Jacques,et al.  MULTICOMPONENT GEOCHEMICAL TRANSPORT MODELING USING HYDRUS‐1D AND HP11 , 2006 .

[12]  W. Bond,et al.  Approximate solutions for cation transport during unsteady, unsaturated soil water flow. , 1990 .

[13]  David R. Anderson,et al.  Multimodel Inference , 2004 .

[14]  C. Tiedeman,et al.  Effective Groundwater Model Calibration: With Analysis of Data, Sensitivities, Predictions, and Uncertainty , 2007 .

[15]  J. W. Biggar,et al.  On the Equivalence of Two Conceptual Models for Describing Ion Exchange During Transport Through an Aggregated Oxisol , 1984 .

[16]  Cation Competition in a Natural Subsurface Material: Modeling of Sorption Equilibria , 2000 .

[17]  Luis S. Pereira,et al.  Field evaluation of a multicomponent solute transport model in soils irrigated with saline waters , 2011 .

[18]  Diederik Jacques,et al.  Modeling Coupled Hydrologic and Chemical Processes: Long‐Term Uranium Transport following Phosphorus Fertilization , 2008 .

[19]  R. Hills,et al.  Parameter estimation of hydraulic properties from one‐step outflow data , 1992 .

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

[21]  Mary C. Hill,et al.  UCODE_2005 and six other computer codes for universal sensitivity analysis, calibration, and uncertainty evaluation constructed using the JUPITER API , 2006 .

[22]  D Mallants,et al.  Operator-splitting errors in coupled reactive transport codes for transient variably saturated flow and contaminant transport in layered soil profiles. , 2006, Journal of contaminant hydrology.

[23]  K. Verburg,et al.  Comparison of Methods for Predicting Ternary Exchange from Binary Isotherms , 1997 .

[24]  M. V. Genuchten,et al.  Selected HYDRUS modules for modeling subsurface flow and contaminant transport as influenced by biological processes at various scales , 2009, Biologia.

[25]  J. J. McKeown.  Non-Linear Parameter Estimation , 1972 .

[26]  Miroslav Šejna,et al.  Development and Applications of the HYDRUS and STANMOD Software Packages and Related Codes , 2008 .

[27]  W. Bond,et al.  On the Rothmund-Kornfeld Description of Cation Exchange , 1995 .

[28]  Chris Smith,et al.  Absorption of gypsum solution by a potassic soil: a data set , 2008 .

[29]  W. Bond,et al.  Extraction Procedure for Determining Solution and Exchangeable Ions on the Same Soil Sample , 1989 .

[30]  R. Mansell,et al.  An algorithm for generating cation exchange isotherms from binary selectivity coefficients , 2001 .

[31]  Diederik Jacques,et al.  Modelling coupled water flow, solute transport and geochemical reactions affecting heavy metal migration in a podzol soil , 2008 .

[32]  Javier Samper,et al.  Inverse problem of multicomponent reactive chemical transport in porous media: Formulation and applications , 2004 .

[33]  G. E. Rayment,et al.  Australian laboratory handbook of soil and water chemical methods. , 1992 .

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

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

[36]  Jirka Šimůnek,et al.  Multicomponent solute transport in soil lysimeters irrigated with waters of different quality , 2006 .

[37]  D. Smiles,et al.  A survey of the cation content of piggery effluents and some consequences of their use to irrigate soils , 2004 .

[38]  M. T. van Genuchten,et al.  Estimating Unsaturated Soil Hydraulic Properties from Tension Disc Infiltrometer Data by Numerical Inversion , 1996 .