Search USGSSearch

Geology topics

R. J. Hunt

Publications and source records attributed to R. J. Hunt.

At least 19 recordsLinked to original sources

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

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

SWB: A modified Thornthwaite-Mather Soil-Water-Balance code for estimating groundwater recharge

A Soil-Water-Balance (SWB) computer code has been developed to calculate spatial and temporal variations in groundwater recharge. The SWB model calculates recharge by use of commonly available geographic information system (GIS) data layers in combination with tabular climatological data. The code is based on a modified Thornthwaite-Mather soil-water-balance approach, with components of the soil-water balance calculated at a daily timestep. Recharge calculations are made on a rectangular grid of computational elements that may be easily imported into a regional groundwater-flow model. Recharge estimates calculated by the code may be output as daily, monthly, or annual values.

Illinois, Michigan, Wisconsin

Relating groundwater to seasonal wetlands in southeastern Wisconsin, USA

Historically, drier types of wetlands have been difficult to characterize and are not well researched. Nonetheless, they are considered to reflect the precipitation history with little, if any, regard for possible relation to groundwater. Two seasonal coastal wetland types (wet prairie, sedge meadow) were investigated during three growing seasons at three sites in the Lake Michigan Basin, Wisconsin, USA. The six seasonal wetlands were characterized using standard soil and vegetation techniques and groundwater measurements from the shallow and deep systems. They all met wetland hydrology criteria (e.g., water within 30 cm of land surface for 5% of the growing season) during the early portion of the growing season despite the lack of appreciable regional groundwater discharge into the wetland root zones. Although root-zone duration analyses did not fit a lognormal distribution previously noted in groundwater-dominated wetlands, they were able to discriminate between the plant communities and showed that wet prairie communities had shorter durations of continuous soil saturation than sedge meadow communities. These results demonstrate that the relative rates of groundwater outflows can be important for wetland hydrology and resulting wetland type. Thus, regional stresses to the shallow groundwater system such as pumping or low Great Lake levels can be expected to affect even drier wetland types. ?? Springer-Verlag 2008.

Hydrogeology Journal

Importance of unsaturated zone flow for simulating recharge in a humid climate

Transient recharge to the water table is often not well understood or quantified. Two approaches for simulating transient recharge in a ground water flow model were investigated using the Trout Lake watershed in north-central Wisconsin: (1) a traditional approach of adding recharge directly to the water table and (2) routing the same volume of water through an unsaturated zone column to the water table. Areas with thin (less than 1 m) unsaturated zones showed little difference in timing of recharge between the two approaches; when water was routed through the unsaturated zone, however, less recharge was delivered to the water table and more discharge occurred to the surface because recharge direction and magnitude changed when the water table rose to the land surface. Areas with a thick (15 to 26 m) unsaturated zone were characterized by multimonth lags between infiltration and recharge, and, in some cases, wetting fronts from precipitation events during the fall overtook and mixed with infiltration from the previous spring snowmelt. Thus, in thicker unsaturated zones, the volume of water infiltrated was properly simulated using the traditional approach, but the timing was different from simulations that included unsaturated zone flow. Routing of rejected recharge and ground water discharge at land surface to surface water features also provided a better simulation of the observed flow regime in a stream at the basin outlet. These results demonstrate that consideration of flow through the unsaturated zone may be important when simulating transient ground water flow in humid climates with shallow water tables.

Ground Water

Comparison of total mercury and methylmercury cycling at five sites using the small watershed approach

The small watershed approach is well-suited but underutilized in mercury research. We applied the small watershed approach to investigate total mercury (THg) and methylmercury (MeHg) dynamics in streamwater at the five diverse forested headwater catchments of the US Geological Survey Water, Energy, and Biogeochemical Budgets (WEBB) program. At all sites, baseflow THg was generally less than 1 ng L −1 and MeHg was less than 0.2 ng L −1 . THg and MeHg concentrations increased with streamflow, so export was primarily episodic. At three sites, THg and MeHg concentration and export were dominated by the particulate fraction in association with POC at high flows, with maximum THg (MeHg) concentrations of 94 (2.56) ng L −1 at Sleepers River, Vermont; 112 (0.75) ng L −1 at Rio Icacos, Puerto Rico; and 55 (0.80) ng L −1 at Panola Mt., Georgia. Filtered (<0.7 μm) THg increased more modestly with flow in association with the hydrophobic acid fraction (HPOA) of DOC, with maximum filtered THg concentrations near 5 ng L −1 at both Sleepers and Icacos. At Andrews Creek, Colorado, THg export was also episodic but was dominated by filtered THg, as POC concentrations were low. MeHg typically tracked THg so that each site had a fairly constant MeHg/THg ratio, which ranged from near zero at Andrews to 15% at the low-relief, groundwater-dominated Allequash Creek, Wisconsin. Allequash was the only site with filtered MeHg consistently above detection, and the filtered fraction dominated both THg and MeHg. Relative to inputs in wet deposition, watershed retention of THg (minus any subsequent volatilization) was 96.6% at Allequash, 60% at Sleepers, and 83% at Andrews. Icacos had a net export of THg, possibly due to historic gold mining or frequent disturbance from landslides. Quantification and interpretation of Hg dynamics was facilitated by the small watershed approach with emphasis on event sampling.

Environmental Pollution

Effects of climate and land management change on streamflow in the driftless area of Wisconsin

Baseflow and precipitation in the Kickapoo River Watershed, located in the Driftless Area of Wisconsin, exhibit a step increase around 1970, similar to minimum and median flows in many other central and eastern USA streams. Potential effects on streamflow due to climatic and land management changes were evaluated by comparing volumetric changes in the hydrologic budget before and after 1970. Increases in precipitation do not fully account for the increase in baseflow, which appears to be offset by a volumetric decrease in stormflow. This suggests that factors that influence the partitioning of precipitation into overland runoff or infiltration have changed. A transition from relatively more intensive to relatively less intensive agricultural land use is generally associated with higher infiltration rates, and likely influences partitioning of flow. Changes in agricultural land management practices in the Driftless Area, which began in the mid-1930s, do not coincide with the abrupt increase in baseflow around 1970. Instead, the timing of hydrologic change appears to coincide with changes in precipitation, whereas the magnitude of the change in baseflow and stormflow was likely amplified by changes in agricultural land management. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Hydrology

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

Using stable isotopes of water and strontium to investigate the hydrology of a natural and a constructed wetland

Wetlands cannot exist without water, but wetland hydrology is difficult to characterize. As a result, compensatory wetland mitigation often only assumes the proper hydrology has been created. In this study, water sources and mass transfer processes in a natural and constructed wetland complex were investigated using isotopes of water and strontium. Water isotope profiles in the saturated zone revealed that the natural wetland and one site in the constructed wetland were primarily fed by ground water; profiles in another constructed wetland site showed recent rain was the predominant source of water in the root zone. Water isotopes in the capillary fringe indicated that the residence time for rain is less in the natural wetland than in the constructed wetland, thus transpiration (an important water sink) was greater in the natural wetland. Strontium isotopes showed a systematic difference between the natural and constructed wetlands that we attribute to the presence or absence of peat. In the peat‐rich natural wetland, δ 87 Sr in the pore water increased along the flowline due to preferential weathering of minerals containing radiogenic Sr in response to elevated Fe concentrations in the water. In the constructed wetland, where peat thickness was thin and Fe concentrations in water were negligible, δ 87 Sr did not increase along the flowline. The source of the peat (on‐site or off‐site derived) applied in the constructed wetland controlled the δ 87 Sr at the top of the profile, but the effects were restricted by strong cation exchange in the underlying fluvial sediments. Based on the results of this study, neither constructed wetland site duplicated the water source and weathering environment of the adjoining natural wetland. Moreover, stable isotopes were shown to be effective tools for investigating wetlands and gaining insight not easily obtained using non‐isotopic techniques. These tools have potential widespread application to wetlands that have distinct isotopic endmember sources.

Ground Water

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

Scale effects of hydrostratigraphy and recharge zonation on base flow

Uncertainty regarding spatial variations of model parameters often results in the simplifying assumption that parameters are spatially uniform. However, spatial variability may be important in resource assessment and model calibration. In this paper, a methodology is presented for estimating a critical basin size, above which base flows appear to be relatively less sensitive to the spatial distribution of recharge and hydraulic conductivity, and below which base flows are relatively more sensitive to this spatial variability. Application of the method is illustrated for a watershed that exhibits distinct infiltration patterns and hydrostratigraphic layering. A ground water flow model (MODFLOW) and a parameter estimation code (UCODE) were used to evaluate the influence of recharge zonation and hydrostratigraphic layering on base flow distribution. Optimization after removing spatial recharge variability from the calibrated model altered base flow simulations up to 53% in watersheds smaller than 40 km2. Merging six hydrostratigraphic units into one unit with average properties increased base flow residuals up to 83% in basins smaller than 50 km2. Base flow residuals changed <5% in watersheds larger than 40 and 50 km2 when recharge and hydrostratigraphy were simplified, respectively; thus, the critical basin size for the example area is ???40 to 50 km2. Once identified for an area, a critical basin size could be used to guide the scale of future investigations. By ensuring that parameter discretization needed to capture base flow distribution is commensurate with the scope of the investigation, uncertainty caused by overextending uniform parameterization or by estimating extra parameter values is reduced. ?? 2006 National Ground Water Association.

Ground Water

Estimating recharge rates with analytic element models and parameter estimation

Quantifying the spatial and temporal distribution of recharge is usually a prerequisite for effective ground water flow modeling. In this study, an analytic element (AE) code (GFLOW) was used with a nonlinear parameter estimation code (UCODE) to quantify the spatial and temporal distribution of recharge using measured base flows as calibration targets. The ease and flexibility of AE model construction and evaluation make this approach well suited for recharge estimation. An AE flow model of an undeveloped watershed in northern Wisconsin was optimized to match median annual base flows at four stream gages for 1996 to 2000 to demonstrate the approach. Initial optimizations that assumed a constant distributed recharge rate provided good matches (within 5%) to most of the annual base flow estimates, but discrepancies of >12% at certain gages suggested that a single value of recharge for the entire watershed is inappropriate. Subsequent optimizations that allowed for spatially distributed recharge zones based on the distribution of vegetation types improved the fit and confirmed that vegetation can influence spatial recharge variability in this watershed. Temporally, the annual recharge values varied >2.5-fold between 1996 and 2000 during which there was an observed 1.7-fold difference in annual precipitation, underscoring the influence of nonclimatic factors on interannual recharge variability for regional flow modeling. The final recharge values compared favorably with more labor-intensive field measurements of recharge and results from studies, supporting the utility of using linked AE-parameter estimation codes for recharge estimation. Copyright ?? 2005 The Author(s).

Ground Water

Investigating surface water-well interaction using stable isotope ratios of water

Because surface water can be a source of undesirable water quality in a drinking water well , an understanding of the amount of surface water and its travel time to the well is needed to assess a well's vulnerability. Stable isotope ratios of oxygen in river water at the City of La Crosse, Wisconsin, show peak-to-peak seasonal variation greater than 4‰ in 2001 and 2002. This seasonal signal was identified in 7 of 13 city municipal wells, indicating that these 7 wells have appreciable surface water contributions and are potentially vulnerable to contaminants in the surface water. When looking at wells with more than 6 sampling events, a larger variation in δ 18 O compositions correlated with a larger fraction of surface water, suggesting that samples collected for oxygen isotopic composition over time may be useful for identifying the vulnerability to surface water influence even if a local meteoric water line is not available. A time series of δ 18 O from one of the municipal wells and from a piezometerlocated between the river and the municipal well showed that the travel time of flood water to the municipal well was approximately 2 months; non-flood arrival times were on the order of 9 months. Four independent methods were also used to assess time of travel. Three methods (groundwater temperature arrival times at the intermediate piezometer, virus-culture results, and particle tracking using a numerical groundwater-flow model) yielded flood and non-flood travel times of less than 1 year for this site. Age dating of one groundwater sample using 3 H– 3 He methods estimated an age longer than 1 year, but was likely confounded by deviations from piston flow as noted by others. Chlorofluorocarbons and SF 6 analyses were not useful at this site due to degradation and contamination, respectively. This work illustrates the utility of stable hydrogen and oxygen isotope ratios of water to determine the contribution and travel time of surface water in groundwater, and demonstrates the importance of using multiple methods to improve estimates for time of travel of 1 year or less.

Journal of Hydrology

Vulnerability of drinking-water wells in La Crosse, Wisconsin, to enteric-virus contamination from surface water contributions

Human enteric viruses can contaminate municipal drinking-water wells, but few studies have examined the routes by which viruses enter these wells. In the present study, the objective was to monitor the municipal wells of La Crosse, Wisconsin, for enteric viruses and determine whether the amount of Mississippi River water infiltrating the wells was related to the frequency of virus detection. From March 2001 to February 2002, one river water site and four wells predicted by hydrogeological modeling to have variable degrees of surface water contributions were sampled monthly for enteric viruses, microbial indicators of sanitary quality, and oxygen and hydrogen isotopes. 18O/ 16O and 2H/1H ratios were used to determine the level of surface water contributions. All samples were collected prior to chlorination at the wellhead. By reverse transcription-PCR (RT-PCR), 24 of 48 municipal well water samples (50%) were positive for enteric viruses, including enteroviruses, rotavirus, hepatitis A virus (HAV), and noroviruses. Of 12 river water samples, 10 (83%) were virus positive by RT-PCR. Viable enteroviruses were not detected by cell culture in the well samples, although three well samples were positive for culturable HAV. Enteroviruses detected in the wells by RT-PCR were identified as several serotypes of echoviruses and group A and group B coxsackieviruses. None of the well water samples was positive for indicators of sanitary quality, namely male-specific and somatic coliphages, total coliform bacteria, Escherichia coli, and fecal enterococci. Contrary to expectations, viruses were found in all wells regardless of the level of surface water contributions. This result suggests that there were other unidentified sources, in addition to surface water, responsible for the contamination.

Applied and Environmental Microbiology

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

Development of Maximum Considered Earthquake Ground Motion Maps

The 1997 NEHRP Recommended Provisions for Seismic Regulations for New Buildings use a design procedure that is based on spectral response acceleration rather than the traditional peak ground acceleration, peak ground velocity, or zone factors. The spectral response accelerations are obtained from maps prepared following the recommendations of the Building Seismic Safety Council's (BSSC) Seismic Design Procedures Group (SDPG). The SDPG-recommended maps, the Maximum Considered Earthquake (MCE) Ground Motion Maps, are based on the U.S. Geological Survey (USGS) probabilistic hazard maps with additional modifications incorporating deterministic ground motions in selected areas and the application of engineering judgement. The MCE ground motion maps included with the 1997 NEHRP Provisions also serve as the basis for the ground motion maps used in the seismic design portions of the 2000 International Building Code and the 2000 International Residential Code . Additionally the design maps prepared for the 1997 NEHRP Provisions , combined with selected USGS probabilistic maps, are used with the 1997 NEHRP Guidelines for the Seismic Rehabilitation of Buildings .

Earthquake Spectra

The importance of ground water in the Great Lakes Region

Ground water is a major natural resource in the Great Lakes Region because it indirectly contributes more than 50 percent of the stream discharge to the Great Lakes. In addition, ground water is the source of drinking water for millions of people in the region, is an important source of supply for agriculture and many industries, and provides a relatively uniform supply of water in some ecologically sensitive areas to sustain plant and animal species. Therefore, to improve our understanding of water-resources issues in the Great Lakes Region, it is important to have a better understanding of the role that ground water plays in the overall hydrologic system of the lakes. The main ground-water resources issues in the Great Lakes Region are related to the amount of ground water, the interaction of ground water and surface water, changes in ground-water quality as development expands, and ecosystem health related to quantity and quality of water. Issues related to the amount of ground water Although the amount of water in the Great Lakes Region is vast, issues related to relatively small quantities of water are being raised more and more often. For example, even though the amount of ground water pumped in the region is small compared to the total amount of water present, ground water is an important source of public-water supply as well as an important source of supply for industrial, agricultural, and domestic needs. Less clearly understood, however, is the relation between the amount of streamflow discharging to the Great Lakes and the large portion of that flow that originates as ground water. The implications of this understanding for water- and land-use practices and, in turn, their effects on water quantity and quality, have not been fully incorporated into a policy framework. To help include information about the implications of the role that ground water plays in addressing regional water issues, a comprehensive analysis of indirect ground-water discharge to the Great Lakes is needed. Direct ground-water discharge to the Great Lakes is not a large factor in water-budget analyses for the Great Lakes. Locally, however, direct ground-water discharge to the Great Lakes may be important, even though the rates and places of discharge are not well known. A long-term evaluation of direct ground-water discharge to the Great Lakes would help place this hydrologic process in proper perspective. Near-shore areas with high rates of direct ground-water discharge may provide valuable habitat for aquatic organisms. Issues related to the interaction of ground water and surface water Withdrawal of ground water removes that water from the watershed when it is consumptively used or when the return flow is discharged to another drainage basin. Under these circumstances, pumping ground water constitutes a diversion of Great Lakes water. Alternatively, ground-water withdrawal could have the opposite effect of diverting ground-water flow into the watershed by altering the ground-water divides. In particular, as withdrawals associated with urban expansion increase, more accurate data on the amount and effects of ground-water use need to be collected. Data on the amounts of ground water pumped both within the watershed and outside, but near the watershed boundaries needs to be collected and evaluated for potential diversion of water to or from the Great Lakes. It is currently thought that both irrigation and ground-water withdrawals near the watershed boundaries constitute relatively small amounts of water; however, both rapidly changing farming practices and rapidly expanding urban communities could alter these amounts in a relatively short timeframe, especially during drought periods. At present, the effects of ground-water withdrawals have been quantified in detail at only a few urban locations. In addition to quantifying the amount of water pumped out of aquifers, it is also important to improve our knowledge of the amount of water that is recharging them. Ground-water recharge rates estimated in earlier studies cover only a small part of the Great Lakes Region. A comprehensive study of ground-water recharge rates for the entire watershed is needed to more completely determine the role of ground water in the hydrologic budget of the Great Lakes. Issues related to changes in ground-water quality as development expands Ground-water quality is as important as quantity for most water uses. As ground-water development proceeds, the possibility of altering the quality of ground water increases. The quality of ground water can be altered when water levels are drawn below the layer that confines the aquifer or by inducing water of lesser quality into an aquifer. Many local studies of these problems have been conducted, but few regional-scale analyses of changes in ground-water quality as a result of ground-water development have been done. Issues related to ecosystem health and quantity and quality of ground water Ground water is essential to maintain wetlands and to provide healthy habitat for other aquatic systems. Wetland hydrology is widely recognized as the primary influence on wetland ecology, development, and persistence, and information about hydrology is essential to understanding and quantifying wetland functions and processes. Studies of the role of ground water in selected wetlands in a range of physiographic settings throughout the Great Lakes watershed are needed to more fully understand the role of wetlands in the Great Lakes Region.

Lake Erie, Lake Huron, Lake Michigan, Lake Ontario

Simulation of the shallow hydrologic system in the vicinity of Middle Genesee Lake, Wisconsin, using analytic elements and parameter estimation

Middle Genesee Lake is a ground-water flow-through lake located in a developing area in southeastern Wisconsin. Because the lake is in good connection with the shallow ground-water system, hydrologic stresses to the shallow ground-water system could adversely affect the lake system. In order to assess the effects of potential stresses on the lake, a study was completed by the U.S. Geological Survey, in cooperation with the Middle Genesee Lake Management District. The objective of the study was to identify areas that contribute ground water to the lake and estimate the hydrologic budget of the lake and hydraulic parameters affecting ground-water flow. A two-dimensional, steady-state analytic element model of the lake and surrounding area was developed using the computer code GFLOW. A parameter estimation model, UCODE, was used to optimize the calibration to measured water levels and streamflow. The calibrated model was used to evaluate the effect of three hypothetical stress scenarios on the stage of Middle Genesee Lake; the simulations were linked to UCODE, which formally incorporated parameter uncertainty into 95-percent confidence intervals around the simulated value. The scenarios included: (1) pumping from upgradient irrigation wells, (2) pumping from Lower Genesee Lake to lower lake levels, and (3) reduction in recharge resulting from development. The results of the simulations demonstrated that lake levels could be affected by hydrologic stresses in the shallow hydrologic system, with effects ranging from a 2.7 feet decline in lake stage resulting from pumping in Lower Genesee Lake to a 0.1 feet decline in lake stage from development in part of the upgradient recharge area. The range of lake stage decline increased when parameter uncertainty was included, from a decline of 3.1 feet for pumping from Lower Genesee Lake to no reduction in lake stage for the development in the recharge area. Whereas these simulated effects are within the natural variation in lake stage, they represent a systematic reduction of ground-water flow to the lake. Therefore, these hypothetical stresses are expected to establish a new, lower, baseline lake stage over which the natural variation due to climatic effects are added and subtracted.

Wisconsin