CXTFIT/Excel-A modular adaptable code for parameter estimation, sensitivity analysis and uncertainty analysis for laboratory or field tracer experiments

Guoping Tang, Melanie A. Mayes, Jack C. Parker, Philip M. Jardine

Research output: Contribution to journalArticlepeer-review

63 Scopus citations

Abstract

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.

Original languageEnglish
Pages (from-to)1200-1209
Number of pages10
JournalComputers and Geosciences
Volume36
Issue number9
DOIs
StatePublished - Sep 2010

Funding

We appreciate the review comments provided by Dr. Philip Meyer and another anonymous reviewer, which significantly helped improve this paper. This research was funded by the U.S. Department of Energy, Office of Science, Office of the Biological and Environmental Research. Oak Ridge National Laboratory is managed by the University of Tennessee-Battelle, LLC , under contract DE-AC05-00OR22725 with the US DOE.

Keywords

  • Equilibrium/nonequilibrium convection dispersion equation
  • Monte Carlo analysis
  • Penalty function
  • Response surface
  • Weighted nonlinear least squares

Fingerprint

Dive into the research topics of 'CXTFIT/Excel-A modular adaptable code for parameter estimation, sensitivity analysis and uncertainty analysis for laboratory or field tracer experiments'. Together they form a unique fingerprint.

Cite this