Geology topics

D. T. Feinstein

Publications and source records attributed to D. T. Feinstein.

8 recordsLinked to original sources

USGS2012

Simulation of an urban ground-water-flow system in the Menomonee Valley, Milwaukee, Wisconsin using analytic element modeling

A single-layer, steady-state analytic element model was constructed to simulate shallow ground-water flow in the Menomonee Valley, an old industrial center southwest of downtown Milwaukee, Wisconsin. Project objectives were to develop an understanding of the shallow ground-water flow system and identify primary receptors of recharge to the valley. The analytic element model simulates flow in a 18.3 m (60 ft) thick layer of estuarine and alluvial sediments and man-made fill that comprises the shallow aquifer across the valley. The thin, laterally extensive nature of the shallow aquifer suggests horizontal-flow predominates, thus the system can appropriately be modeled with the Dupuit-Forchheimer approximation in an analytic element model. The model was calibrated to the measured baseflow increase between two USGS gages on the Menomonee River, 90 head measurements taken in and around the valley during December 1999, and vertical gradients measured at five locations under the river and estuary in the valley. Recent construction of the Milwaukee Metropolitan Sewer District Inline Storage System (ISS) in the Silurian dolomite under the Menomonee Valley has locally lowered heads in the dolomite appreciably, below levels caused by historic pumping. The ISS is a regional hydraulic sink which removes water from the bedrock even during dry weather. The potential effect on flow directions in the shallow aquifer of dry-weather infiltration to the ISS was evaluated by adjusting the resistance of the line-sink strings representing the ISS in the model to allow infiltration from 0 to 100% of the reported 9,500 m 3 /d. The best fit to calibration targets was found between 60% (5,700 m 3 /d) and 80% (7,600 m 3 /d) of the reported dry-weather infiltration. At 60% infiltration, 65% of the recharge falling on the valley terminates at the ISS and 35% at the Menomonee River and estuary. At 80% infiltration, 73% of the recharge terminates at the ISS, and 27% at the river and estuary. Model simulations suggest that the ISS has an greater influence on the shallow ground-water flow in the eastern half of valley as compared to the western half. Preliminary three-dimensional simulations using the numerical MODFLOW code show good agreement with the single-layer simulation and supports its use in evaluating the shallow system. Copyright ASCE 2004.

Conference Paper↗

USGS2012

Development and application of a groundwater/surface-water flow model using MODFLOW-NWT for the Upper Fox River Basin, southeastern Wisconsin

The Fox River is a 199-mile-long tributary to the Illinois River within the Mississippi River Basin in the states of Wisconsin and Illinois. For the purposes of this study the Upper Fox River Basin is defined as the topographic basin that extends from the upstream boundary of the Fox River Basin to a large wetland complex in south-central Waukesha County called the Vernon Marsh. The objectives for the study are to (1) develop a baseline study of groundwater conditions and groundwater/surface-water interactions in the shallow aquifer system of the Upper Fox River Basin, (2) develop a tool for evaluating possible alternative water-supply options for communities in Waukesha County, and (3) contribute to the methodology of groundwater-flow modeling by applying the recently published U.S. Geological Survey MODFLOW-NWT computer code, (a Newton formulation of MODFLOW-2005 intended for solving difficulties involving drying and rewetting nonlinearities of the unconfined groundwater-flow equation) to overcome computational problems connected with fine-scaled simulation of shallow aquifer systems by means of thin model layers. To simulate groundwater conditions, a MODFLOW grid is constructed with thin layers and small cell dimensions (125 feet per side). This nonlinear unconfined problem incorporates the streamflow/lake (SFR/LAK) packages to represent groundwater/surface-water interactions, which yields an unstable solution sensitive to initial conditions when solved using the Picard-based preconditioned-gradient (PCG2) solver. A particular problem is the presence of many isolated wet water-table cells over dry cells, causing the simulated water table to assume unrealistically high values. Attempts to work around the problem by converting to confined conditions or converting active to inactive cells introduce unacceptable bias. Application of MODFLOW-NWT overcomes numerical problem by smoothing the transition from wet to dry cells and keeps all cells active. The simulation is insensitive to initial conditions and the water-table trend is smooth across layers. The MODFLOW-NWT code permits rigorous calibration and also robust application of the model to transient scenarios. Runtimes on a 64-bit computer are kept reasonably short by use of updated initial conditions and informed choices of solver parameters. The shallow aquifer system consists of unconsolidated material of varying thickness over Silurian dolomite. The unconsolidated material, largely of glacial origin, contains fine-textured and coarse-textured deposits that vary in permeability over short distances. This study at least partly encompasses the inevitable uncertainty in the hydraulic conductivity zones by developing two models—one favors the continuity of fine-grained deposits and a second favors the continuity of coarse-grained deposits. The separate calibration processes for the fine-favored and coarse-favored models using MODFLOW-NWT and the nonlinear regression algorithms in the parameter estimation (PEST) code produce distinct parameter values for hydraulic conductivity zones, storage parameters, and streambed conductance zones. Both models are applied to a hypothetical scenario involving 27 "riparian" wells completed adjacent to the river channel and open to the shallow aquifer systems along a 10-mile stretch of the Fox River. The results suggest that a riparian well system withdrawing about 9 million gallons per day would induce about one-third to one-half its total discharge from the river, and that this riverbank inducement would appreciably limit drawdown around the hypothetical wells.

Wisconsin↗

USGS2011

Application of the Local Grid Refinement package to an inset model simulating the interaction of lakes, wells, and shallow groundwater, northwestern Waukesha County, Wisconsin

Groundwater use from shallow, high-capacity wells is expected to increase across southeastern Wisconsin in the next decade (2010-2020), owing to residential and business growth and the need for shallow water to be blended with deeper water of lesser quality, containing, for example, excessive levels of radium. However, this increased pumping has the potential to affect surface-water features. A previously developed regional groundwater-flow model for southeastern Wisconsin was used as the starting point for a new model to characterize the hydrology of part of northwestern Waukesha County, with a particular focus on the relation between the shallow aquifer and several area lakes. An inset MODFLOW model was embedded in an updated version of the original regional model. Modifications made within the inset model domain include finer grid resolution; representation of Beaver, Pine, and North Lakes by use of the LAK3 package in MODFLOW; and representation of selected stream reaches with the SFR package. Additionally, the inset model is actively linked to the regional model by use of the recently released Local Grid Refinement package for MODFLOW-2005, which allows changes at the regional scale to propagate to the local scale and vice versa. The calibrated inset model was used to simulate the hydrologic system in the Chenequa area under various weather and pumping conditions. The simulated model results for base conditions show that groundwater is the largest inflow component for Beaver Lake (equal to 59 percent of total inflow). For Pine and North Lakes, it is still an important component (equal, respectively, to 16 and 5 percent of total inflow), but for both lakes it is less than the contribution from precipitation and surface water. Severe drought conditions (simulated in a rough way by reducing both precipitation and recharge rates for 5 years to two-thirds of base values) cause correspondingly severe reductions in lake stage and flows. The addition of a test well south of Chenequa at a pumping rate of 47 gal/min from a horizon approximately 200 feet below land surface has little effect on lake stages or flows even after 5 years of pumping. In these scenarios, the stage and the surface-water outflow from Pine Lake are simulated to decrease by only 0.03 feet and 3 percent, respectively, relative to base conditions. Likely explanations for these limited effects are the modest pumping rate simulated, the depth of the test well, and the large transmissivity of the unconsolidated aquifer, which allows the well to draw water from upstream along the bedrock valley and to capture inflow from the Bark River. However, if the pumping rate of the test well is assumed to increase to 200 gal/min, the decrease in simulated Pine Lake outflow is appreciably larger, dropping by 14 percent relative to base-flow conditions.

Wisconsin↗

USGS2010

Regional groundwater-flow model of the Lake Michigan Basin in support of Great Lakes Basin water availability and use studies

A regional groundwater-flow model of the Lake Michigan Basin and surrounding areas has been developed in support of the Great Lakes Basin Pilot project under the U.S. Geological Survey's National Water Availability and Use Program. The transient 2-million-cell model incorporates multiple aquifers and pumping centers that create water-level drawdown that extends into deep saline waters. The 20-layer model simulates the exchange between a dense surface-water network and heterogeneous glacial deposits overlying stratified bedrock of the Wisconsin/Kankakee Arches and Michigan Basin in the Lower and Upper Peninsulas of Michigan; eastern Wisconsin; northern Indiana; and northeastern Illinois. The model is used to quantify changes in the groundwater system in response to pumping and variations in recharge from 1864 to 2005. Model results quantify the sources of water to major pumping centers, illustrate the dynamics of the groundwater system, and yield measures of water availability useful for water-resources management in the region. This report is a complete description of the methods and datasets used to develop the regional model, the underlying conceptual model, and model inputs, including specified values of material properties and the assignment of external and internal boundary conditions. The report also documents the application of the SEAWAT-2000 program for variable-density flow; it details the approach, advanced methods, and results associated with calibration through nonlinear regression using the PEST program; presents the water-level, drawdown, and groundwater flows for various geographic subregions and aquifer systems; and provides analyses of the effects of pumping from shallow and deep wells on sources of water to wells, the migration of groundwater divides, and direct and indirect groundwater discharge to Lake Michigan. The report considers the role of unconfined conditions at the regional scale as well as the influence of salinity on groundwater flow. Lastly, it describes several categories of limitations and discusses ways of extending the regional model to address issues at the local scale. Results of the simulations portray a regional groundwater-flow system that, over time, has largely maintained its natural predevelopment configuration but that locally has been strongly affected by well withdrawals. The quantity of rainfall in the Lake Michigan Basin and adjacent areas supports a dense surface-water network and recharge rates consistent with generally shallow water tables and predominantly shallow groundwater flow. At the regional scale, pumping has not caused major modifications of the shallow flow system, but it has resulted in decreases in base flow to streams and in direct discharge to Lake Michigan (about 2 percent of the groundwater discharged and about 0.5 cubic foot per second per mile of shoreline). On the other hand, well withdrawals have caused major reversals in regional flow patterns around pumping centers in deep, confined aquifers - most noticeably in the Cambrian-Ordovician aquifer system on the west side of Lake Michigan near the cities of Green Bay and Milwaukee in eastern Wisconsin, and around Chicago in northeastern Illinois, as well as in some shallow bedrock aquifers (for example, in the Marshall aquifer near Lansing, Mich.). The reversals in flow have been accompanied by large drawdowns with consequent local decrease in storage. On the west side of Lake Michigan, groundwater withdrawals have caused appreciable migration of the deep groundwater divides. Before the advent of pumping, the deep Lake Michigan groundwater-basin boundaries extended west of the Lake Michigan surface-water basin boundary, in some places by tens of miles. Over time, the pumping centers have replaced Lake Michigan as the regional sink for the deep flow system. The regional model is intended to support the framework pilot study of water availability and use for the Great Lakes Basin (Reeves, in press).

Illinois, Indiana, Michigan, Wisconsin↗

USGS2006

The importance of diverse data types to calibrate a watershed model of the Trout Lake Basin, Northern Wisconsin, USA

As part of the USGS Water, Energy, and Biogeochemical Budgets project and the NSF Long-Term Ecological Research work, a parameter estimation code was used to calibrate a deterministic groundwater flow model of the Trout Lake Basin in northern Wisconsin. Observations included traditional calibration targets (head, lake stage, and baseflow observations) as well as unconventional targets such as groundwater flows to and from lakes, depth of a lake water plume, and time of travel. The unconventional data types were important for parameter estimation convergence and allowed the development of a more detailed parameterization capable of resolving model objectives with well-constrained parameter values. Independent estimates of groundwater inflow to lakes were most important for constraining lakebed leakance and the depth of the lake water plume was important for determining hydraulic conductivity and conceptual aquifer layering. The most important target overall, however, was a conventional regional baseflow target that led to correct distribution of flow between sub-basins and the regional system during model calibration. The use of an automated parameter estimation code: (1) facilitated the calibration process by providing a quantitative assessment of the model's ability to match disparate observed data types; and (2) allowed assessment of the influence of observed targets on the calibration process. The model calibration required the use of a 'universal' parameter estimation code in order to include all types of observations in the objective function. The methods described in this paper help address issues of watershed complexity and non-uniqueness common to deterministic watershed models. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

USGS2005

Simulating ground water-lake interactions: Approaches and insights

Approaches for modeling lake-ground water interactions have evolved significantly from early simulations that used fixed lake stages specified as constant head to sophisticated LAK packages for MODFLOW. Although model input can be complex, the LAK package capabilities and output are superior to methods that rely on a fixed lake stage and compare well to other simple methods where lake stage can be calculated. Regardless of the approach, guidelines presented here for model grid size, location of three-dimensional flow, and extent of vertical capture can facilitate the construction of appropriately detailed models that simulate important lake-ground water interactions without adding unnecessary complexity. In addition to MODFLOW approaches, lake simulation has been formulated in terms of analytic elements. The analytic element lake package had acceptable agreement with a published LAK1 problem, even though there were differences in the total lake conductance and number of layers used in the two models. The grid size used in the original LAK1 problem, however, violated a grid size guideline presented in this paper. Grid sensitivity analyses demonstrated that an appreciable discrepancy in the distribution of stream and lake flux was related to the large grid size used in the original LAK1 problem. This artifact is expected regardless of MODFLOW LAK package used. When the grid size was reduced, a finite-difference formulation approached the analytic element results. These insights and guidelines can help ensure that the proper lake simulation tool is being selected and applied.

Ground Water↗

USGS2005

The vertical hydraulic conductivity of an aquitard at two spatial scales

Aquitards protect underlying aquifers from contaminants and limit recharge to those aquifers. Understanding the mechanisms and quantity of ground water flow across aquitards to underlying aquifers is essential for ground water planning and assessment. We present results of laboratory testing for shale hydraulic conductivities, a methodology for determining the vertical hydraulic conductivity (Kv) of aquitards at regional scales and demonstrate the importance of discrete flow pathways across aquitards. A regional shale aquitard in southeastern Wisconsin, the Maquoketa Formation, was studied to define the role that an aquitard plays in a regional ground water flow system. Calibration of a regional ground water flow model for southeastern Wisconsin using both predevelopment steady-state and transient targets suggested that the regional Kv of the Maquoketa Formation is 1.8 ?? 10 -11 m/s. The core-scale measurements of the Kv of the Maquoketa Formation range from 1.8 ?? 10-14 to 4.1 ?? 10-12 m/s. Flow through some additional pathways in the shale, potential fractures or open boreholes, can explain the apparent increase of the regional-scale Kv. Based on well logs, erosional windows or high-conductivity zones seem unlikely pathways. Fractures cutting through the entire thickness of the shale spaced 5 km apart with an aperture of 50 microns could provide enough flow across the aquitard to match that provided by an equivalent bulk Kv of 1.8 ?? 10-11 m/s. In a similar fashion, only 50 wells of 0.1 m radius open to aquifers above and below the shale and evenly spaced 10 km apart across southeastern Wisconsin can match the model Kv. Copyright ?? 2005 National Ground Water Association.

Ground Water↗

USGS2004

Simulation of ground-water flow, surface-water flow, and a deep sewer tunnel system in the Menomonee Valley, Milwaukee, Wisconsin

Numerical models were constructed for simulation of ground-water flow in the Menomonee Valley Brownfield, in Milwaukee, Wisconsin. An understanding of ground-water flow is necessary to develop an efficient program to sample ground water for contaminants. Models were constructed in a stepwise fashion, beginning with a regional, single-layer, analytic-element model (GFLOW code) that provided boundary conditions for a local, eight layer, finite-difference model (MODFLOW code) centered on the Menomonee Valley Brownfield. The primary source of ground water to the models is recharge over the model domains; primary sinks for ground water within the models are surface-water features and the Milwaukee Metropolitan Sewerage District Inline Storage System (ISS). Calibration targets were hydraulic heads, surface-water fluxes, vertical gradients, and ground-water infiltration to the ISS. Simulation of ground-water flow by use of the MODFLOW model indicates that about 73 percent of recharge within the MODFLOW domain circulates to the ISS and 27 percent discharges to gaining surface-water bodies. In addition, infiltration to the ISS comes from the following sources: 36 percent from recharge within the model domain, 45 percent from lateral flow into the domain, 15 percent from Lake Michigan, and 4 percent from other surface-water bodies. Particle tracking reveals that the median traveltime from the recharge point to surface-water features is 8 years; the median time to the ISS is 255 years. The traveltimes to the ISS are least over the northern part of the valley, where dolomite is near the land surface. The distribution of traveltimes in the MODFLOW simulation is greatly influenced by the effective porosity values assigned to the various lithologies.

Wisconsin↗