We implemented the widely used CXTFIT code in Excel to provide flexibility and added sensitivity and uncertainty analysis functions to improve transport parameter estimation and to facilitate model discrimination for multi-tracer experiments on structured soils. Analytical solutions for one-dimensional equilibrium and nonequilibrium convection dispersion equations were coded as VBA functions so that they could be used as ordinary math functions in Excel for forward predictions. Macros with user-friendly interfaces were developed for optimization, sensitivity analysis, uncertainty analysis, error propagation, response surface calculation, and Monte Carlo analysis. As a result, any parameter with transformations (e.g., dimensionless, log-transformed, species-dependent reactions, etc.) could be estimated with uncertainty and sensitivity quantification for multiple tracer data at multiple locations and times. Prior information and observation errors could be incorporated into the weighted nonlinear least squares method with a penalty function. Users are able to change selected parameter values and view the results via embedded graphics, resulting in a flexible tool applicable to modeling transport processes and to teaching students about parameter estimation. The code was verified by comparing to a number of benchmarks with CXTFIT 2.0. It was applied to improve parameter estimation for four typical tracer experiment data sets in the literature using multi-model evaluation and comparison. Additional examples were included to illustrate the flexibilities and advantages of CXTFIT/Excel. The VBA macros were designed for general purpose and could be used for any parameter estimation/model calibration when the forward solution is implemented in Excel. A step-by-step tutorial, example Excel files and the code are provided as supplemental material.
[1]
Bill Jelen,et al.
VBA and Macros for Microsoft Excel
,
2004
.
[2]
William H. Press,et al.
Numerical Recipes: FORTRAN
,
1988
.
[3]
Jack C. Parker,et al.
Improving parameter estimation for column experiments by multi-model evaluation and comparison.
,
2009
.
[4]
William H. Press,et al.
Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition
,
1987
.
[5]
Eungyu Park,et al.
Dissolved plume attenuation with DNAPL source remediation, aqueous decay and volatilization--analytical solution, model calibration and prediction uncertainty.
,
2008,
Journal of contaminant hydrology.
[6]
J. Neter,et al.
Applied Linear Regression Models
,
1983
.
[7]
W. Press,et al.
Numerical Recipes in Fortran: The Art of Scientific Computing.@@@Numerical Recipes in C: The Art of Scientific Computing.
,
1994
.
[8]
Dani Or,et al.
Nonlinear Parameter Estimation Using Spreadsheet Software
,
1998
.
[9]
M. V. Genuchten,et al.
Flux-Averaged and Volume-Averaged Concentrations in Continuum Approaches to Solute Transport
,
1984
.
[10]
C. Tiedeman,et al.
Effective Groundwater Model Calibration: With Analysis of Data, Sensitivities, Predictions, and Uncertainty
,
2007
.
[11]
C. Tiedeman,et al.
Effective Groundwater Model Calibration
,
2007
.
[12]
N. Bohr.
MONTE CARLO METHODS IN GEOPHYSICAL INVERSE PROBLEMS
,
2002
.
[13]
W. Gander,et al.
Adaptive Quadrature—Revisited
,
2000
.
[14]
Jack C. Parker,et al.
Determining transport parameters from laboratory and field tracer experiments
,
1984
.