ModelGIS

DESCRIPTION

Introduction to ModelGIS
  ModelGIS Design
  Comparison to Traditional Methods
  GIS Interface
  MODFLOW Interface
ModelGIS Technical Paper
  ModelGIS Introduction
  ModelGIS Case Study
  ModelGIS Interface
  Other Utilities and Future Modifications
  Conclusion
ModelGIS Requirements

Introduction to ModelGIS

ModelGIS is a comprehensive, menu-driven interface for ground-water modeling within the ARC/INFO GIS software on UNIX workstations. ModelGIS allows model grids, model input, and model output to be created, stored, and displayed within ARC/INFO. A complete preprocessor for the MODFLOW ground-water code is included, and required components for other finite-difference models (model grid, lists, and matrices) can be created. This powerful interface increases modeling accuracy and efficiency, particularly when model regridding is required.

ModelGIS Design
There are two major components to ModelGIS. The first component is GWArc, a menu-driven interface to ARC/INFO for designing variable-spaced model grids, creating list and matrix data that form the basis of model input, resolving unwanted feature overlaps (such as a well in an inactive cell), and displaying model input and output with ARC/INFO's sophisticated map-making capabilities. The second component is GWPre, a menu-driven interface that compiles lists, matrices, and all other required parameters into a model input file. GWPre contains a complete interface for the MODFLOW ground-water flow code. Complete interfaces for other flow and transport codes are planned for GWPre in the future.

Comparison to Traditional Methods
The construction of a ground-water flow model can be a labor-intensive process. Many types of data that vary in space and/or time must be assigned by model layer, row and column. Historically, this has been accomplished via light tables and text editors, and more recently, via CAD-based preprocessors. Use of GIS software allows more sophisticated assignment of model features. Data are stored in the GIS in the form of arcs, points, polygons and surfaces. For instance, rivers are stored as arcs, with each reach assigned parameters such as width, bottom elevation, sediment thickness, etc. ModelGIS not only determines the proper grid cells to assign each river by using ARC/INFO overlay functions, but automatically determines the length of each river segment within each cell and calculates conductance. It also assigns bottom elevation and/or stage based on interpolation from known points. These functions are rapid, accurate, and repeatable. If the model grid is altered, the same overlay operations can be quickly performed again.

GIS Interface
The most difficult aspect of model input is to generate the numerous matrices and lists. Matrix data are assigned at each grid cell such as hydraulic conductivity or aquifer bottom. List data are assigned at select grid cells such as rivers or wells. ModelGIS allows both types of data to be created by overlaying ARC/INFO data with the model grid.

MODFLOW Interface
Each MODFLOW package has its own form menu for entering parameters and joining the matrices and lists created with the GIS interface. ModelGIS then compiles all the input file for the various packages into a comprehensive model input data set.

ModelGIS, An Interface for the USGS Three-Dimensional Finite-Difference Model MODFLOW Technical Paper

ModelGIS was originally developed to facilitate the creation of a local groundwater model from a regional model. The code evolved from GWZOOM, which provided the technology to accommodate the desire to create submodels from regional models. ModelGIS is an interface used to link ARC/INFO to groundwater models. The USGS code MODFLOW is a modular finite-difference computer simulation code capable of simulating groundwater flow in three dimensions. It is an excellent tool used for the prediction of flow in complex groundwater systems. MODFLOW's modular structure allows for flexibility in its use, however the development of a complex model is limited by the lack of an efficient data preprocessor. The hydrogeologic complexity and detail of some sites facilitates the need for an interface to a batch-oriented groundwater simulator. The use of a Geographic Information System (GIS) provides a powerful and efficient means of data preparation and visualization of simulation results. Arc Macro Language (AML) and FORTRAN 77 are used to create ModelGIS which generates model grids, model layer elevations, aquifer properties, surface water data, and model output. ModelGIS may be used to compute regions of transmissivity and vertical leakance from hydrostratigraphic zones containing discrete properties of sand, silt, and clay. A case study is presented demonstrating the capabilities and flexibility of ModelGIS. This new approach allows for the development of a more accurate and complex numerical model that previously could not have been discretized to such a high level of detail.

ModelGIS Introduction
Groundwater models are powerful tools used to analyze many groundwater problems. For example, they may be used for reconnaissance studies prior to field investigations, data interpretation after field programs, predictive studies involving estimation of future groundwater system behavior, and, in combination with contaminant transport models, for evaluation of plume migration and remediation. An understanding of the hydrogeologic conditions of an aquifer system is necessary for the conceptualization and development of groundwater flow models. In order for groundwater flow models to be able to simulate actual aquifer behavior, they must have comprehensive data structures that allow for the utilization of various types of data describing the hydrogeologic conditions of the system. When, for example, parameter values describing the variability within the groundwater flow system show temporal trends in hydrologic stresses as well as past and future trends in water levels, there is a need to combine the spacial characteristics of the aquifer into a numerical representation. The complexities of the data make a Geographic Information System (GIS) a valuable tool for use in the development of models because of its ability to create, store, analyze, and present relational data.

ModelGIS Case Study
In order to demonstrate the functionality of ModelGIS, a groundwater flow model was constructed utilizing hydrogeologic data at the Davis Well Field. The purpose of this model is to provide possible explanations for water-quality degradation in two wells in the Davis Well Field. Located in the southwestern corner of Tennessee, the Davis Well Field is 1 of 10 producing fields operated by MLGW. Production began in 1971, with an estimated 13 million gallons per day (MGD) currently being withdrawn from 14 municipal wells. Production is from the Memphis aquifer. The hydrogeology in the vicinity of the site is quite complex, consisting of consolidated and semi-consolidated sediments of Cretaceous and Tertiary age dipping westward into the north-central portion of the Mississippi embayment, and southward along the axis of the trough of a broad syncline that approximates the Mississippi River. Loess covered bluffs serve both as a physiographic and hydrostratigraphic boundary between the fluvial deposits of the Gulf Coastal Plain to the west and the alluvial deposits of the Mississippi Alluvial Plain to the east (see figure below).

ModelGIS Interface
Converting data from ARC/INFO into MODFLOW data sets is accomplished using a series of ARC/INFO commands, AML functions, and FORTRAN 77 programming. The ModelGIS interface is written using the standards that ESRI developed for ArcTools (ESRI's menu interface to ARC/INFO). The figure below shows the main ModelGIS menu. GWArc converts coverages to model-based data while GWPre assembles the data into MODFLOW data packages. An option to execute the MODFLOW program is also available. Site features are defined in the GIS as points (wells), lines (rivers and streams), polygons (lakes and hydraulic zones), and three-dimensional surfaces (groundwater elevations and hydrostratigraphic surfaces).

Once the GWArc module is invoked, a menu appears containing several options as shown in the figure below. The construction of the finite-difference model grid is usually the first activity associated with the development of a model. The user is presented with the choice of creating a grid interactively or from a MODFLOW row and column file (found in the BCF package). ModelGIS has the flexibility to create a grid using uniform or variably-spaced cells as well as a standard or rotated grid. The MODFLOW row and column file was used to construct the model grid.

After creation of the grid, the user may construct model-based coverages by selecting the Model Data option in the GWArc menu. The creation of matrix-oriented data is accomplished by intersecting a polygon coverage or a three-dimensional surface with the model grid. The Polygon Intersection menu in the figure below allows the user to create matrix-oriented data from polygon coverages. An option is also presented to assign each grid block a value based on the grid center rule or a maximum area rule. Other types of data stored as polygons such as recharge zones, wetlands, water bodies or hydraulic conductivity zones may also be intersected with the grid using the same procedure.

Three-dimensional surfaces most commonly created are potentiometric surfaces and hydrostratigraphic surfaces. The creation of matrix-based coverages is accomplished by the intersection of polygon coverages or three-dimensional surfaces with the model grid.

List-oriented data is created by intersecting a point, line, or polygon coverage with the model grid. Wells, rivers, water bodies, or hydrogeologic zones may be conceptually modeled by using MODFLOW's list-oriented well, river, or drain packages. For MODFLOW's river package, layer, row, column, stage, and conductance are required for each cell defined as a river. In reality, river stage gage measurements may only exist at a few locations along a river. ModelGIS has the ability to linearly interpolate along the river between data points in order to provide data points for all intersected grid cells (see figure below). The calculation of the conductance term in the river package is a function of the length of the river reach; therefore, ModelGIS can easily determine the value of conductance. Also, if the river is defined as a polygon, linear interpolation along a river's center line and extrapolation to adjacent grid cells may be performed. Linear interpolation between stages and calculation of conductances is one of the more powerful features of ModelGIS. These options permit the determination of new interpolated values based on revised stages as well as accurate recreation of the river package based on a new model discretization.

Another unique feature of ModelGIS allows for the generation of three-dimensional surfaces from existing hydrogeologic data. A five-layer model was developed for the Davis Well Field from the interpretation of geophysical logs and two-dimensional profiles. The water table aquifer, a sand and gravel unit and the Jackson-upper Claiborne confining unit, predominantly comprised of clay, were modeled as layers one and two, respectively. The Memphis Sand aquifer, modeled as three layers, is approximately 800-feet thick. Layers three and four each have a thickness of 100 feet while layer five is approximately 600-feet thick. Within the Memphis Sand aquifer, well screens are variably placed due to the heterogeneity of this aquifer unit. 

The figure below illustrates the process used to generate three-dimensional surfaces from two-dimensional profiles using ModelGIS. To determine layer elevations, data is extracted from each profile using a sampling interval, usually the same as the smallest grid spacing. Then, the elevation data is related to the same interval along the profile line (plan view) and stored as a point attribute. The resulting point coverage is used to generate a three-dimensional surface.

Once the structural aspects of the aquifer are included in the model, parameters describing the flow through the model layers must be assigned. The creation of a hydraulic conductivity surface based on a harmonic mean of the percentage of stratigraphic types is performed using a ModelGIS utility. A model profile containing the different hydrostratigraphic zones is intersected with a data sample interval which is then linked to the location along the profile. This procedure is similar to the previous method of converting profile data into a surface. A three-dimensional surface containing the harmonic mean data is then intersected with the model grid and multiplied by the hydraulic conductivity coverage. The use of the profile utility allows for the creation of data sets using data from all available sources thus resulting in a model most accurately depicting the actual conditions.

After all model data sets are created, the user can invoke the GWPre module. This module allows for the preparation of matrix and list-oriented data into MODFLOW format (see figure below). The construction of the basic package, flow package, and other packages are performed in this module. The basic module allows the selection of other packages, determination of the grid size (based on grid dimensions in the GIS workspace), linking of boundary, and starting head arrays as well as input of other non-matrix and non-list based model data. The BCF module has features similar to that of the basic module as well as built-in logic to assist in the creation of the BCF package. Depending on the layer type, ModelGIS automatically determines the files required for the package. Once all packages have been assembled, the user returns to the main module and executes MODFLOW.

Once the model is executed, ModelGIS can also be used to evaluate model results and subsequently modify existing data. These changes may be made to individual model grid cells or an entire matrix based on model results. For example, a model run may generate dry nodes or cells where the water level is below the model layer in discrete locations due to the oscillation of the solution package. The resultant water levels can be read by ModelGIS. Based on the stratigraphic elevation of the layer, the water level at the locations of the dry nodes may be adjusted to ensure grid cell saturation. This new matrix of water level data can then be used for the next model run resulting in fewer dry nodes and a more rapid convergence.

Other Utilities and Future Modifications
Several modules currently exist in ModelGIS that were not used in this case study but can be used in other model applications. Two specific modules include a coupled surface-water and groundwater flow model (MODBRANCH), and a parameter estimation module (PEST). MODBRANCH consists of two programs, BRANCH a one-dimensional numerical model used to simulate unsteady flow in open channel networks developed by the USGS and using the groundwater flow model MODFLOW. PEST, a nonlinear parameter estimation module, adjusts model parameters in order to minimize the discrepancies between observed and measured data. Modules currently under development include interfaces for the MODPATH particle pathline and endpoint processor and a new version of MODFLOW that simulates three-dimensional contaminant transport (MODFLOWT).

Conclusion
ModelGIS is a powerful tool that can be used to develop data sets for most finite-difference models. It provides a dynamic connection between data stored in a GIS and the generation of complex data sets. Using simple menus, calculations of model parameters, construction of model layers, and development of other model features may be performed. From the interpretation of model results, ModelGIS can be used to modify individual grid cells or an entire matrix. Based on this case study, the development of a sophisticated numerical groundwater model can be easily developed using ModelGIS. The use of ModelGIS expedites the conversion of the hydrogeologic framework from GIS point, line, and polygon coverages into model-based data. The modular construction of ModelGIS allows for the evolution of additional utilities such as the contaminant transport code MODFLOWT and the parameter estimation code PEST, that assist in environmental studies.

ModelGIS Requirements: ModelGIS is currently available for SUN Workstations and requires ARC/INFO and the TIN module for a variety of operations. It does not utilize any other "add-on" modules for ARC/INFO (such as GRID or COGO). A version is being developed for Windows NT.