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

This study obtains hydrologic, transport and geochemical parameters from an experimental dataset involving transient unsaturated water flow and solute transport, using an automatic inverse optimization algorithm. The data come from controlled laboratory experiments on the transport of major cations during water absorption into horizontal soil columns. Experimental data consisted of the depth profiles of water contents, Cl concentrations, and 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. 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.