USGS ScienceSearch

USGS · 70026122

Hydrostratigraphic modeling of a complex, glacial-drift aquifer system for importation into MODFLOW

Abstract

Deposition from at least three episodes of glaciation left a complex glacial-drift aquifer system in central Illinois. The deepest and largest of these aquifers, the Sankoty-Mahomet Aquifer, occupies the lower part of a buried bedrock valley and supplies water to communities throughout central Illinois. Thin, discontinuous aquifers are present within glacial drift overlying the Sankoty-Mahomet Aquifer. This study was commissioned by local governments to identify possible areas where a regional water supply could be obtained from the aquifer with minimal adverse impacts on existing users. Geologic information from more than 2200 existing water well logs was supplemented with new data from 28 test borings, water level measurements in 430 wells, and 35 km of surface geophysical profiles. A three-dimensional (3-D) hydrostratigraphic model was developed using a contouring software package, a geographic information system (GIS), and the 3-D geologic modeling package, EarthVision??. The hydrostratigraphy of the glacial-drift sequence was depicted as seven uneven and discontinuous layers, which could be viewed from an infinite number of horizontal and vertical slices and as solid models of any layer. Several iterations were required before the 3-D model presented a reasonable depiction of the aquifer system. Layers from the resultant hydrostratigraphic model were imported into MODFLOW, where they were modified into continuous layers. This approach of developing a 3-D hydrostratigraphic model can be applied to other areas where complex aquifer systems are to be modeled and is also useful in helping lay audiences visualize aquifer systems.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

B.L. Herzog, D.R. Larson, C.C. Abert, S.D. Wilson, G.S. Roadcap. 2005-12-13. Hydrostratigraphic modeling of a complex, glacial-drift aquifer system for importation into MODFLOW. https://doi.org/10.1111/j.1745-6584.2003.tb02568.x

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Evaluating micrometeorological estimates of groundwater discharge from Great Basin desert playas

Groundwater availability studies in the arid southwestern United States traditionally have assumed that groundwater discharge by evapotranspiration (ETg) from desert playas is a significant component of the groundwater budget. However, desert playa ETg rates are poorly constrained by Bowen Ratio energy budget (BREB) and eddy-covariance (EC) micrometeorological measurement approaches. Best attempts by previous studies to constrain ETg from desert playas have resulted in ETg rates that are within the measurement error of micrometeorological approaches. This study uses numerical models to further constrain desert playa ETg rates that are within the measurement error of BREB and EC approaches, and to evaluate the effect of hydraulic properties and salinity-based groundwater-density contrasts on desert playa ETg rates. Numerical models simulated ETg rates from desert playas in Death Valley, California and Dixie Valley, Nevada. Results indicate that actual ETg rates from desert playas are significantly below the uncertainty thresholds of BREB- and EC-based micrometeorological measurements. Discharge from desert playas likely contributes less than 2 percent of total groundwater discharge from Dixie and Death Valleys, which suggests discharge from desert playas also is negligible in other basins. Simulation results also show that ETg from desert playas primarily is limited by differences in hydraulic properties between alluvial fan and playa sediments and, to a lesser extent, by salinity-based groundwater density contrasts.

Ground Water

A semi-structured MODFLOW-USG model to evaluate local water sources to wells for decision support

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A “semi-structured” approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Illinois, Indiana, Michigan, Wisconsin

Extending the MODPATH algorithm to rectangular unstructured grids

The recent release of MODFLOW-USG, which allows model grids to have irregular, unstructured connections, requires a modification of the particle-tracking algorithm used by MODPATH. This paper describes a modification of the semi-analytical particle-tracking algorithm used by MODPATH that allows it to be extended to rectangular-based unstructured grids by dividing grid cells with multi-cell face connections into sub-cells. The new method will be incorporated in the next version of MODPATH which is currently under development.

Ground Water