MODFLOW

From Wikipedia, the free encyclopedia

MODFLOW simulation
Enlarge
MODFLOW simulation

MODFLOW is the U.S. Geological Survey modular finite-difference flow model, which is a computer code that solves the groundwater flow equation. The program is used by hydrogeologists to simulate the flow of groundwater through aquifers. The code is free software, written primarily in Fortran, and can compile and run on DOS, Windows or Unix-like operating systems.

3-dimensional grid
Enlarge
3-dimensional grid

Since its original development in the early 1980s, the USGS have released five major releases, and is now considered to be the de facto standard code for aquifer simulation. Currently, there are at least four actively developed commercial and non-commercial graphical user interfaces for MODFLOW.

Contents

[edit] Groundwater flow equation

The governing partial differential equation used in MODFLOW is:

\frac{\partial}{\partial x} \left[ K_{xx} \frac{\partial h}{\partial x} \right] + \frac{\partial}{\partial y} \left[ K_{yy} \frac{\partial h}{\partial y} \right] + \frac{\partial}{\partial z} \left[ K_{zz} \frac{\partial h}{\partial z} \right] + W = S_{S} \frac{\partial h}{\partial t}

where

Kxx, Kyy and Kzz are the values of hydraulic conductivity along the x, y, and z coordinate axes (L/T)
h is the potentiometric head (L)
W is a volumetric flux per unit volume representing sources and/or sinks of water, where negative values are extractions, and positive values are injections (T-1)
SS is the specific storage of the porous material (L-1); and
t is time (T)

[edit] Limitations

\mathbf{K} = \begin{bmatrix} K_{xx} & 0 & 0 \\ 0 & K_{yy} & 0 \\ 0 & 0 & K_{zz}\end{bmatrix} \

[edit] Versions

Cover image from McDonald & Harbaugh (1984), which illustrates a computer surrounded by modules and arrays used by MODFLOW
Enlarge
Cover image from McDonald & Harbaugh (1984), which illustrates a computer surrounded by modules and arrays used by MODFLOW

[edit] MODFLOW (original version)

The first version was released on December 28, 1983 (McDonald & Harbaugh, 1984), and was coded entirely in FORTRAN 66. The source code for this version is not available on the internet.

[edit] MODFLOW-88

This version of MODFLOW was rewritten in FORTRAN 77, and was originally released in 1987 July 24. The current version of MODFLOW-88 is 2.6, released in 1996 September 20.

MODPATH, was initially developed in 1989 to post-process the steady-state MODFLOW-88 data to determine three-dimensional pathlines of particles. This innovation has been indispensable for the fields of contaminant hydrogeology. It is still used as a post-processor in recent versions of MODFLOW.

A separate program, MODFLOWP, was developed in 1992 to estimate various parameters used in MODFLOW. This program was eventually built into MODFLOW-2000.

[edit] MODFLOW-96

MODFLOW-96 (version 3.0) was originally released in 1996 December 3, and is a cleaned-up and revised continuation of MODFLOW-88. There are three final releases of MODFLOW-96:

Several graphical interfaces were first developed using the MODFLOW-96 code.

[edit] MODFLOW-2000

MODFLOW-2000 (version 1.0; version numbering was reset) was released on 2000 July 20, which merged MODFLOWP and HYDMOD codes into the main program. Many new packages and enhancements were also included, including new solvers, stream and saturated flow packages. The internal design concepts also changed from previous versions, such that packages, processes and modules are distinct. This version is coded in a mixture of FORTRAN 77, Fortran 90, and C. MODFLOW-2000 can also be compiled for parallel computing, which can allow multiple processors to be used to increase model complexity and/or reduce simulation time.

There are four current version of MODFLOW-2000:

[edit] MODFLOW-2005

This is the present development version of MODFLOW, and does not yet have a major commercial graphical interface. It is written primarily in Fortran 90 and C. The current versions include:

[edit] Packages

Name Long name Version introduced
Required packages
BAS Basic original
OC Output Control original
Groundwater flow packages
BCF Block-Centered Flow original
LPF Layer-Property Flow MODFLOW-2000 (1.0)
HUF Hydrogeologic Unit Flow MODFLOW-2000 (1.1)
Solvers
SIP Strongly Implicit Procedure original
SOR Slice Successive Over-Relaxation original
DE4 Direct Solver MODFLOW-88 (2.5)
PCG Preconditioned Conjugate-Gradient MODFLOW-88
LMG Link-AMG MODFLOW-2000 (1.4)
GMG Geometric Multigrid Solver MODFLOW-2000 (1.15.00)
Head-dependent flux boundary packages
GHB General-Head Boundary original
DRN Drain original
DRT Drain Return MODFLOW-2000 (1.1)
RIV River original
EVT Evapotranspiration original
ETS Evapotranspiration Segments MODFLOW-2000 (1.1)
RES Reservoir MODFLOW-88 (2.6)
LAK Lake MODFLOW-2000 (1.1)
STR Stream MODFLOW-88
SFR Streamflow-Routing MODFLOW-2000 (1.14.00)
MNW Multi-Node, Drawdown-Limited Well MODFLOW-2000 (1.11)
DAF DAFLOW MODFLOW-2000 (1.11)
Other stress packages
WEL Well original
CHD Constant-Head Boundary MODFLOW-88
FHB Flow and Head Boundary MODFLOW-96 (3.2)
RCH Recharge original
HFB Horizontal Flow Barrier MODFLOW-88
Miscellaneous packages
DIS Discretization MODFLOW-2000 (1.0)
GAG Gage MODFLOW-2000
HYD HYDMOD MODFLOW-2000 (1.1)
IBS Interbed-Storage MODFLOW-88
KDEP Hydraulic-Conductivity Depth-Dependence Capability MODFLOW-2000 (1.12)
LMT Link-MT3DMS MODFLOW-2000 (1.5)
LVDA Model-Layer Variable-Direction Horizontal Anisotropy Capability MODFLOW-2000 (1.12)
SUB Subsidence and Aquifer-System Compaction MODFLOW-2000 (1.12)
UTL Utility original
Observation process input files
OBS Input File For All Observations MODFLOW-2000
HOB Head-Observation MODFLOW-2000
DROB Drain Observation MODFLOW-2000
DTOB Drain Return Observation MODFLOW-2000
RVOB River Observation MODFLOW-2000
GBOB General-Head-Boundary Observation MODFLOW-2000
CHOB Constant-Head Flow Observation MODFLOW-2000
ADV Advective-Transport Observation MODFLOW-2000 (1.0)
STOB Stream Observation MODFLOW-2000
Obsolete packages
GFD General Finite-Difference MODFLOW-88 to 96
TLK Transient Leakage MODFLOW-88 to 96
  •   due to licensing restrictions, the USGS is no longer able to publicly distribute the Algebraic Multi-Grid

[edit] Commercial graphical interfaces

There are several graphical commercial versions of MODFLOW, which include the compiled MODFLOW code with modifications. These programs aid the input of data for creating MODFLOW models, however there is limited compatibility between the programs.

All current versions of these products run only on Microsoft Windows, however previous versions of GMS (up to Version 3.1) were compiled for several Unix platforms.

The USGS has a free MODFLOW-GUI, which is updated often, however it requires a non-free version of Argus ONE.

[edit] References

  • Harbaugh, A.W., Banta, E.R., Hill, M.C., and McDonald, M.G., 2000, MODFLOW-2000, the U.S. Geological Survey modular ground-water model — User guide to modularization concepts and the Ground-Water Flow Process: U.S. Geological Survey Open-File Report 00-92, 121 p. [1]
  • Harbaugh, A.W., and McDonald, M.G., 1996a, User's documentation for MODFLOW-96, an update to the U.S. Geological Survey modular finite-difference ground-water flow model: U.S. Geological Survey Open-File Report 96-485, 56 p. PDF, PS.
  • Harbaugh, A.W., and McDonald, M.G., 1996, Programmer's documentation for MODFLOW-96, an update to the U.S. Geological Survey modular finite-difference ground-water flow model: U.S. Geological Survey Open-File Report 96-486, 220 p. PDF, PS.
  • McDonald, M.G., and Harbaugh, A.W., 1984, A modular three-dimensional finite-difference ground-water flow model: U.S. Geological Survey Open-File Report 83-875, 528 p.
  • McDonald, M.G., and Harbaugh, A.W., 1988, A modular three-dimensional finite-difference ground-water flow model: U.S. Geological Survey Techniques of Water-Resources Investigations, book 6, chap. A1, 586 p. English, Chinese.

[edit] External links

In other languages