Jump to main content.


Center for Subsurface Modeling Support

Home

Image: GIS Research
Important Information

Potential Limitations of Four Domenico-Based Fate and Transport Models

  • BIOCHLOR
  • BIOSCREEN
  • FOOTPRINT
  • REMChlor
  • REMFuel

The Center for Subsurface Modeling Support (CSMoS) provides public domain groundwater and vadose zone modeling software and services to public agencies and private companies throughout the United States. CSMoS is located at the Robert S. Kerr Center for Environmental Research in Ada, Oklahoma.

The primary goals of CSMoS are to provide direct technical support to EPA and state decision makers in subsurface model applications, and to manage and support the groundwater models and databases resulting from the research at EPA’s National Risk Management Research Laboratory (NRMRL).

The research includes the:

  • Fate and transport of contaminants in the subsurface
  • Development of methodologies for protection and restoration of groundwater quality
  • Evaluation of subsurface remedial technologies

As a result, a major focus of CSMoS is the coordinating the use of models for risk assessment, site characterization, remedial activities, wellhead protection, and the application of Geographic Information Systems. In these ways, CSMoS performs an active role in protecting, restoring, and preserving our nation's groundwater resources.

Modeling Services

In its effort to apply models for a better understanding and better resolutions of groundwater problems, CSMoS integrates the expertise of numerous individuals and organizations in all aspects of the environmental field. The center is supported by NRMRL scientists and engineers whose specialties include hydrogeology, chemistry, soil science, biology, environmental engineering, and computer programming. CSMoS provides assistance in the modeling areas of conceptualization, development, application, distribution, and training and education.

Top of page

Technical Assistance

CSMoS is an integral part of the NRMRL's Technical Support Center. CSMoS distributes and services all models and databases developed by NRMRL and provides general support on model application to groundwater and vadose zone problems. Technical assistance activities include developing educational documents, providing training courses, and distributing update notices and other pertinent information for all software developed at NRMRL, as well as software developed under laboratory grants and contracts.

CSMoS provides direct technical assistance for a broad spectrum of modeling applications. Models and databases are available for:

  • Site Characterization – chemical and soil characterization databases provide values for required model input data
  • Geostatistics – models produce statistics and semi-variograms for kriging of site data
  • Saturated Zone Ground Water Flow and Fate/Transport – two- and three-dimensional, single-phase and multiphase flow models provide fate and transport for unconfined and confined aquifers
  • Unsaturated Zone Ground Water Flow and Fate/Transport – one-dimensional single and multiphase flow with fate and transport for the vadose zone
  • Combined Unsaturated/Saturated Flow and Transport – two- and three-dimensional, single-phase and multiphase flow models with fate and transport for the vadose and saturated zones
  • Natural Attenuation – one-dimensional screening models to assess whether natural attenuation of fuel hydrocarbons or chlorinated solvents is occurring at a site and a two-dimensional model for predicting the natural attenuation of fuel hydrocarbons at a site
  • Wellhead Protection Assessment – models that enable steady-state and time-related capture zones to be determined for numerous pumping well scenarios in various site settings

These models can be used to assist in:

  • Performing site characterization
  • Conducting groundwater flow and transport simulations
  • Selecting groundwater remediation options at Resource Conservation and Recovery Act sites, and for Superfund studies
  • Determining wellhead protection areas

Top of page

Download Software Image: Disclaimer

Important Information

Potential Limitations of Four Domenico-Based Fate and Transport Models

  • BIOCHLOR
  • BIOSCREEN
  • FOOTPRINT
  • REMChlor

Free Public Domain Ground Water and Vadose Zone Models

The Center for Subsurface Modeling Support (CSMoS) distributes various public domain groundwater and vadose zone models developed by EPA’s Ground Water and Ecosystem Restoration Division. There are three ways to get models and manuals:

  • Click the links in the table below to view model specifications, descriptions, minimum system requirements, and files available for downloading.
  • Send blank CDs or ZIP disks to CSMoS. CSMoS will copy the requested models or manuals to disks and return them to you.
  • Contact Rob Earle, 580-436-8656, for other arrangements

CSMoS provides free technical assistance for many of the models listed in the table below.

Click on the "Mailing List" tab above to join the CSMoS e-mail list. Receive information about new releases of public domain groundwater models, version updates, bugs, fixes, and other important items.

Want more? Click or the "Search Database" tab above to search our database of commercially available groundwater modeling software.

CSMoS Models
Model Version Release Date Platform Additional Information
2DFATMIC 1.0 August 1997 DOS 2-D subsurface flow/transport
3DFATMIC 1.0 August 1997 DOS 3-D subsurface flow/transport
BIOCHLOR 2.2 June 2002 Windows 95/98/NT, Excel 1-D Domenico screening model
BIOPLUME II 1.1 October 1989 DOS 2-D USGS MOC transport
BIOPLUME III 1.0 September 1997 Windows 95/98 2-D USGS MOC transport with Windows GUI
BIOSCREEN 1.4 July 1997 Windows 95/98/NT, Excel 3-D Domenico transport
CHEMFLO 1.3 August 1990 DOS 1-D vadose zone numerical transport
CHEMFLO-2000   April 2003 Windows
95/98/NT/2000/XP, Linux,
Solaris, Mac OS X
Updated CHEMFLO
CZAEM 1.1 August 1996 DOS Analytic element capture zone model
EvalMNA   December 2011 Windows 95/98/NT/2000/XP Supplemental information to EPA 600/R-11/204
FOOTPRINT 1.0 June 2008 Windows 98/XP 2-D transport of BTEX and ethanol
GEOEAS 1.2.1 April 1991 DOS, UNIX Geostatistical analysis
GEOPACK 1.0.e January 1990 DOS Geostatistical analysis
HSSM-DOS 1.1 April 1994 DOS Multiphase LNAPL flow/transport
HSSM-Windows 1.2.e September 1997 Windows 95/98 Multiphase LNAPL flow/transport
HSSM en Español 1.2.e September 1997 Windows 95/98 Multiphase LNAPL flow/transport
Infiltration Models     Mathcad Unsaturated/saturated zone transport
MODFLOW Manual   February 1993 DOS MODFLOW practice problems
MOFAT 2.0.a May 1991 DOS 2-D multiphase transport
MT3D 1.11 January 1992 DOS 3-D numerical transport
NAPL Simulator 2.0 November 2010 DOS Simulation of NAPL spills
OWL 1.2 March 2004 Windows 95/98/NT/2000/XP Monitoring well locator
PESTAN 4.0   Windows Simulate leaching of pesticides
Regression MNA   December 2011 Windows 95/98/NT/2000/XP Supplemental information to EPA 600/R-11/204
REMFuel 1.0 February 2012
Windows 98 or greater Simulating the effects for fuel hydrocarbons
REMChlor 1.0 December 2007 .Net Simulate transient plum remediation
RETC 1.1 November 1994 DOS Estimate soil model parameters
RITZ 2.12 January 1998 DOS Simulate vadose zone transport
UTCHEM
Exit EPA Disclaimer
9.0 July 2000 Windows 3-D multiphase flow/transport
Virulo 1.0 August 2002 Windows, UNIX, Mac Probabilistic virus leaching model
VLEACH 2.2.a May 2007 Windows 1-D vadose zone leaching model
WhaEM        
WhAEM2000 3.2 June 2005 Windows 95/98/NT/2000/XP Analytical element capture zone model
WHPA 2.2 September 1993 DOS Finite-difference capture zone model

Model Descriptions

2DFATMIC NAPL Simulator
3DFATMIC OWL
BIOCHLOR PESTAN
BIOPLUME II Regression MNA
BIOPLUME III REMFuel
BIOSCREEN REMchlor
CHEMFLO RETC
EvalMNA RITZ
FOOTPRINT UTCHEM
GEOEAS Virulo
GEOPACK VLEACH
HSSM WhAEM
MODFLOW Manual WhAEM2000
MOFAT WHPA
MT3D  

2DFATMIC

Application
This model simulates subsurface flow, fate, and transport of contaminants that are undergoing chemical or biological transformations. The model is applicable to transient conditions in both saturated and unsaturated zones.

Processes
The flow module is a Galerkin finite element solution of Richard's equation. The transport module is a hybrid Lagrangian-Eulerian approach, with an adapted zooming and peak capturing algorithm.

Miscellaneous
This model can almost eliminate spurious oscillation, numerical dispersion, and peak clipping due to advective transport.

Contact
Professor G.T. Yeh
814-863-2931

Top of page

3DFATMIC

Application
This model simulates subsurface flow, fate, and transport of contaminants that are undergoing chemical or biological transformations. The model is applicable to transient conditions in both saturated and unsaturated zones.

Processes
The flow module is a Galerkin finite element solution of Richard's equation. The transport module is a hybrid Lagrangian-Eulerian approach, with an adapted zooming and peak capturing algorithm.

Miscellaneous
This model can almost eliminate spurious oscillation, numerical dispersion, and peak clipping due to advective transport.

Contact
Professor G.T. Yeh
814-863-2931

Top of page

BIOCHLOR

Application
This is a screening model that simulates remediation by natural attenuation of dissolved solvents at chlorinated solvent release sites. BIOCHLOR can be used to simulate solute transport without decay, and solute transport with biodegradation modeled as a sequential first-order process within one or two different reaction zones.

Processes
The program is based on the one-dimensional Domenico analytical solute transport model and includes one-dimensional advection, three-dimensional dispersion, linear adsorption, and biotransformation via reductive dechlorination. Reductive dechlorination is assumed to occur under anaerobic conditions and solvent degradation is assumed to follow a sequential first-order decay process.

Miscellaneous
The program is written as an Excel spreadsheet and, therefore, requires Microsoft Excel to run. The program includes a natural attenuation screening protocol scoring system to help the user determine the potential for reductive dechlorination from site data.

Top of page

BIOPLUME II

Application
This is a two-dimensional contaminant transport model used under the influence of oxygen limited biodegradation.

Processes
Processs include advection, dispersion, sorption, biodegradation (aerobic and anaerobic), and re-aeration.

Miscellaneous
This model is based on the two-dimensional solute transport code USGS-MOC.

Top of page

BIOPLUME III

Application
This is a two-dimensional contaminant transport model for use under the influence of oxygen, nitrate, iron, sulfate, and methanogenic biodegradation.

Processes
Processes include advection, dispersion, sorption, first order decay, biodegradation through instantaneous, first order, zero order, or Monod kinetics.

Miscellaneous
This model is based on the two-dimensional solute transport code USGS-MOC, but has an integrated, Windows-based graphical user interface for seamless preprocessing and output review. Hydrocarbon source and each active electron acceptor are simulated as separate plumes.

Top of page

BIOSCREEN

Application
This is a three-dimensional contaminant transport model for dissolved-phase hydrocarbons in a saturated zone under the influences of oxygen, nitrate, iron, sulfate, and methane limited biodegradation.

Processes
Processes include advection, dispersion, adsorption, first-order decay, and instantaneous reactions under aerobic and anaerobic conditions.

Miscellaneous
This is an easy-to-use screening tool that is based on the Domenico analytical solute transport model.

Top of page

CHEMFLO

Application
This model simulates one-dimensional water and chemical movement in the vadose zone.

Processes
Processes include advection, dispersion, first-order decay, and linear sorption.

Miscellaneous
This is a numerical solution model that handles a variety of boundary conditions.

Top of page

EvalMNA

Application
Comparison of Initial Year of Review Cycle to Final Year of Review Cycle (.xls) (2 pp, 60 KB)

Top of page

FOOTPRINT

Application
This model simulates the length and surface area of BTEX plumes in groundwater, when the plumes are produced from a gasoline spill that contains ethanol.

Processes
This model is based on a modified version of the 1982 Domenico model. This model was modified to allow zero-order decay.

Top of page

GEOEAS

Application
This model performs geostatistical analysis of spatially correlated data.

Components
Components include basic statistics, scatter plots, and linear and nonlinear estimation (kriging).

Top of page

GEOPACK

Application
This model performs geostatistical analysis of spatially correlated data.

Components
Components include basic statistics, variography, and linear and nonlinear estimation (kriging).

Top of page

HSSM

Application
This model simulates flow of light nonaqueous-phase liquid (LNAPL) and transport of a chemical constituent of the LNAPL from the surface to the water table, radial spreading of the LNAPL at the water table, and dissolution and aquifer transport of the chemical constituent.

Processes
Processes include one-dimensional in the vadose zone, radial in the capillary fringe, and two-dimensional vertically averaged analytical solution of the advectiondispersion equation in the saturated zone.

Miscellaneous
This model is based on the KOPT, OILENS, and TSGPLUME models. It is also available in Spanish.

Top of page

MODFLOW Manual

Application
This is an instruction manual for the U.S.Geological Survey’s MODFLOW model.

Components
The manual contains a series of 20 problem sets that illustrate by example the use of MODFLOW, including modeling principles, input/output specifics, available options, rules of thumb, and common modeling mistakes.

Top of page

MOFAT

Application
This is a two-dimensional flow and transport model of three fluid phases: water, nonaqueous-phase liquid, and gas.

Processes
Processes include advection, dispersion, diffusion, sorption, decay, and mass transfer.

Miscellaneous
This model is finite-element; numerous parameters are required.

Top of page

MT3D

Application
This model is a three-dimensional contaminant transport in the saturated zone.

Processes
Processes include advection, dispersion, nonlinear sorption, first-order irreversible decay, and biodegradation.

Miscellaneous
This is a numerical solution model that uses a hybrid method of characteristics and particle tracking. It handles a variety of discretization schemes and boundary conditions and includes the MODFLOW program.

Top of page

NAPL Simulator

Application
This model simulates the contamination of soils and aquifers that results from the release of organic liquids, commonly referred to as nonaqueous-phase liquids (NAPLs).

Processes
The simulator is applicable to three interrelated zones: a vadose zone, which is in contact with the atmosphere; a capillary zone; and a water-table aquifer zone. Three mobile phases are accommodated: water, NAPL, and gas. The three-phase k-S-P submodel accommodates capillary and fluid entrapment hysteresis. NAPL dissolution and volatilization are accounted for through rate-limited mass transfer submodels.

Miscellaneous
The numerical solution is based on a Hermite collocation finite element discretization. The simulator provides an accurate solution of a coupled set of nonlinear partial differential equations that are generated by combining fundamental balance equations with constitutive thermodynamic relationships.

Contact
Joseph Guarnaccia
732-914-2516

Top of page

OWL

Application
This model evaluates existing monitoring well networks and optimizes the selection new monitoring well locations.

Processes
Processes include calculation of plume migration, based on variations in groundwater flow magnitude and direction. Optimal locations are selected, based on simple analysis of plume locations and the existing monitoring well network.

Top of page

PESTAN

Application
This program performs vadose zone modeling of the transport of organic (pesticide) contaminants.

Processes
Processes include advection, dispersion, first-order decay, and linear sorption.

Miscellaneous
This is a screening model; few parameters are required.

Top of page

RegressionMNA

Application
Regression MNA (.xls) (3 pp, 56 KB)

Top of page

REMFuel

Application
This model simulates the transient effects of groundwater source and plume remediation for fuel hydrocarbons

Processes
This is a contaminant source model based on a power function relationship between source mass and source discharge for multiple fuel costituents, and it can consider partial source remediation at any time after the initial release. The source model serves as a time-dependent mass flux boundary condition to the analytical plume model, where flow is assumed to be one-dimensional. The plume model for each fuel component simulates first order sequential decay and production of one daughter species. REMFuel can also simulate zero order or Monod's kinetics for decay of fuel components in the plume. The decay rates and other reaction coefficients are variable functions of time and distance in the plume. This approach allows for flexible simulation of enhanced plume remediation that may be temporary in time, limited in space, and which may have different effects on different contaminant species in the plume.

Top of page

REMchlor

Application
This model simulates the transient effects of groundwater source and plume remediation.

Processes
This is a contaminant source model based on a power-function relationship between source, mass, and discharge. It can consider partial source remediation at any time after release. The source model is a time-dependent, mass-flux boundary condition to the analytical plume model (one-dimensional flow), which simulates first-order sequential decay and production of several species. This model also calculates cancer risks posed by carcinogenic species.

Top of page

RETC

Application
This model estimates the soil-water retention curve, unsaturated hydraulic conductivity, and soil model parameters.

Processes
This model uses the parametric equations of Brooks-Corey and van Genuchten, and the pore-size distribution models of Maulem and Burdine.

Miscellaneous
This is an analytical model; it requires relatively few parameters.

Top of page

RITZ

Application
This program performs vadose zone modeling of the transport of contaminants associated with oily wastes.

Processes
Processes include water movement, volatilization, degradation, sorption, and leaching.

Miscellaneous
This is an analytical model that requires relatively few parameters.

Top of page

UTCHEM Exit EPA Disclaimer

Application
Originally a three-dimensional finite difference model for multiphase flow, transport and chemical flooding, the UTCHEM code has been modified to transform the model into a general-purpose NAPL simulator. Appropriate physical, chemical, and biological process models have been incorporated into the simulator to create a threedimensional multiphase, multicomponent model that is capable of simulating the fate and transport of NAPLs in the saturated and unsaturated zones of aquifers. The model can be used to simulate the actual field operation of remediation activities, such as surfactant remediation and bioremediation, as well as laboratory experiments with large-scale aquifer models.

Processes
UTCHEM is capable of modeling transient and steady-state, three-dimensional flow and mass transport in the groundwater (saturated) and vadose (unsaturated) zones of aquifers. Physical, chemical, and biological process models important in describing the fate and transport of NAPLs in contaminated aquifers have been incorporated into the simulator. These include multiple organic NAPL phases, the dissolution or mobilization of NAPLs by nondilute remedial fluids, chemical and microbiological transformations, and changes in fluid properties as a site is remediated. The model allows for nonequilibrium interphase mass transfer, sorption, geochemical reactions, and the temperature dependence of pertinent chemical and physical properties. It can simulate the flow and transport of remedial fluids whose density, temperature, and viscosity are variable, including surfactants, co-solvents, and other enhancement agents. The biodegradation model includes inhibition, sequential use of electron acceptors, and co-metabolism, and can be used to model a very general class of bioremediation processes.

Miscellaneous
Biodegradation capabilities have been added to describe the transformation of organic contaminants from NAPL sources. This model can accommodate multiple substrates, electron acceptors, and biological species. A new multiphase capillarypressure and relative permeability function has been added to allow the use of either Brooks-Corey or Van Genuchten capillary pressure functions. New organic and tracer components have been added, as well as additional water tracer components and gas-phase tracers. The number of oil/water tracers has been expanded to allow any number of tracer components. The geochemical option has been extended to allow the modeling of any solid or aqueous species. UTCHEM uses a solution scheme analogous to the Implicit Pressure Explicit Saturation (IMPES) routine, where the pressure is solved for implicitly but concentrations instead of saturations are solved in a flash routine. An energy balance equation includes heat flow between the reservoir and the over- and under-burden rocks.

Contact
Dr. Gary Pope
University of Texas, Department of Petroleum Engineering
Austin, Texas

Top of page

Virulo

Application
This is a probabilistic screening model for predicting leaching of viruses in unsaturated soils.

Processes
Monte Carlo is employed to generate ensemble simulations of virus attenuation. The probability of failure is generated to achieve a user-chosen degree of attenuation.

Top of page

VLEACH

Application
This model simulates one-dimensional water and chemical movement in vadose zone.

Processes
Processes include advection, sorption, vapor-phase diffusion, and three-phase equilibration.

Miscellaneous
This is a numerical solution and screening model.

Top of page

WhAEM

Application
This program delineates capture zones and isochrones of groundwater residence time for the purpose of wellhead protection.

Processes
Processes include steady state, homogeneous, isotropic, advection, and dispersion.

Miscellaneous
This model consists of two independent executables: Geographic Analytic Element Preprocessor (GAEP) and Capture Zone Analytic Element Model (CZAEM).

Top of page

WhAEM2000

Application
This model delineates capture zones and isochrones of groundwater residence time for the purpose of source water area delineation and wellhead protection.

Processes
Processes include steady state, homogeneous, isotropic, advection.

Miscellaneous
This model includes a full 32-bit graphical user interface for Windows (95, 98, NT), on-line help, and a tutorial. Working with the WhAEM 2000 document demonstrates modeling practice for a Vincennes, Indiana, case study.

Contact
Stephen R. Kraemer
Phone: 706-355-8340
Fax: 706-355-8302

Top of page

WHPA

Application
This model simulates capture zones for pumping wells.

Processes
Processes include steady state and horizontal flow.

Miscellaneous
This model performs particle tracking with analytical, semi-analytical, and numerical modules.

Top of page

Mailing List

Occasionally, CSMoS sends e-mails to registered users announcing noteworthy information, such as bugs, fixes, and new releases of groundwater models. If you would like to join the mailing list, please complete the form below.

Your Name:
Company Name:
Your Phone Number:
Street Address:
City:
State/Province:
Country:
Zip Code:
E-mail:

Type of Organization:


Select the items for which you’d like to be on the COSMos mailing list:

2D-3D FATMIC
CZAEM
MOFAT
RITZ
BIOCHLOR FOOTPRINT MT3D UTCHEM
BIOPLUME II GEOEAS NAPL Simulator Virulo
BIOPLUME III GEOPACK OWL VLEACH
BIOSCREEN HSSM PESTAN WhAEM
CHEMFLO Infiltration Models REMchlor WhAEM 2000
CHEMFLO 2000 MODFLOW RETC WHPA
Training Courses

Type the number from the image below, to help stop spam (required).

Captcha Image: type 962501 in the box.

Office of Research & Development | National Risk Management Research Laboratory


Local Navigation


Jump to main content.