Geology ReportsSearch

SEARCH · Geology Reports

Results for “Ground Water Modeling Investigation Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20Linked to original sources

Contributing recharge areas to water-supply wells at Wright-Patterson Air Force Base, Ohio

Wright-Patterson Air Force Base, in southwestern Ohio, has operated three well fields--Area B, Skeel Road, and the East Well Fields--to supply potable water for consumption and use for base activities. To protect these well fields from contamination and to comply with the Ohio Wellhead Protection Plan, the Base is developing a wellhead-protection program for the well fields. A three-dimensional, steady-state ground-water-flow model was developed in 1993 to simulate heads in (1) the buried-valley aquifer system that is tapped by the two active well fields, and in (2) an upland bedrock aquifer that may supply water to the wells. An advective particle-tracking algorithm that requires estimated porosities and simulated heads was used to estimate ground-water-flow pathlines and traveltimes to the active well fields. Contributing recharge areas (CRA's)--areas on the water table that contribute water to a well or well field--were generated for 1-, 5-, and 10-year traveltimes. Results from the simulation and subsequent particle tracking indicate that the CRA's for the Skeel Road Well Fields are oval and extend north- ward, toward the Mad River, as pumping at the well field increases. The sizes of the 1-, 5-, and 10-year CRA's of Skeel Road Well Field, under maximum pumping conditions, are approximately 0.5, 1.5 and 3.2 square miles, respectively. The CRA's for the Area B Well Field extend to the north, up the Mad River Valley; as pumping increases at the well field, the CRA's extend up the Mad River Valley under Huffman Dam. The sizes of the 1-, 5-, and 10-year CRA's of Area B Well Field, under maximum pumping conditions, are approximately 0.1, 0.5, and 0.9 square miles, respectively. The CRA's for the East Well Field are affected by nearby streams under average pumping conditions. The sizes of the 1-, 5-, and 10-year CRA's of the East Well Field, under maximum pumping conditions, are approximately 0.2, 1.2, and 2.4 square miles, respectively. However, as pumping increases at the East Well Field, the ground-water-flow model develops numerical instabilities which limit the usefulness of the CRA's. Sensitivity analyses show that variation of horizontal hydraulic conductivity and porosity in the upland bedrock does not affect the CRA's of the Skeel Road Well Field but does have a slight affect on the CRA's of the Area B Well Field. Uncertainties in horizontal hydraulic conductivity and porosity of the valley-train deposits have the largest affect on the size and shape of the CRA's of the Skeel Road Well Field. The position and size of the CRA's of Area B are probably also controlled by induced infiltration from the nearby Mad River and by pumping at the Rohrer's Island Well Field. However, uncertainty in riverbed conductance, which affects induced infiltration, does not significantly affect the size and shape of these CRA's. Pumping centers not included in the ground-water-flow model do not appreciably affect the CRA's of the Area B and Skeel Road Well Fields under normal pumping. The pumping centers, located near Huffman Dam, will probably limit the northern extent of teh CRA's of Area B Well Field under greater than normal pumping conditions. The CRA's of the East Well Field will propagate farther to the northeast and southwest as a result of the increased pumping-related stress to the aquifer system.

Water-Resources Investigations Report

Phosphorus and nitrogen concentrations and loads at Illinois River south of Siloam Springs, Arkansas, 1997-1999

Water-quality sampling consisting of every other month (bimonthly) routine sampling and storm event sampling (six storms annually) is used to estimate annual phosphorus and nitrogen loads at Illinois River south of Siloam Springs, Arkansas. Hydrograph separation allowed assessment of base-flow and surfacerunoff nutrient relations and yield. Discharge and nutrient relations indicate that water quality at Illinois River south of Siloam Springs, Arkansas, is affected by both point and nonpoint sources of contamination. Base-flow phosphorus concentrations decreased with increasing base-flow discharge indicating the dilution of phosphorus in water from point sources. Nitrogen concentrations increased with increasing base-flow discharge, indicating a predominant ground-water source. Nitrogen concentrations at higher base-flow discharges often were greater than median concentrations reported for ground water (from wells and springs) in the Springfield Plateau aquifer. Total estimated phosphorus and nitrogen annual loads for calendar year 1997-1999 using the regression techniques presented in this paper (35 samples) were similar to estimated loads derived from integration techniques (1,033 samples). Flow-weighted nutrient concentrations and nutrient yields at the Illinois River site were about 10 to 100 times greater than national averages for undeveloped basins and at North Sylamore Creek and Cossatot River (considered to be undeveloped basins in Arkansas). Total phosphorus and soluble reactive phosphorus were greater than 10 times and total nitrogen and dissolved nitrite plus nitrate were greater than 10 to 100 times the national and regional averages for undeveloped basins. These results demonstrate the utility of a strategy whereby samples are collected every other month and during selected storm events annually, with use of regression models to estimate nutrient loads. Annual loads of phosphorus and nitrogen estimated using regression techniques could provide similar results to estimates using integration techniques, with much less investment.

Water-Resources Investigations Report

Hydrogeology and simulation of the effects of reclaimed-water application in west Orange and southeast Lake counties, Florida

Wastewater reclamation and reuse has become increasingly popular as water agencies search for alternative water-supply and wastewater-disposal options. Several governmental agencies in central Florida currently use the land-based application of reclaimed water (wastewater that has been treated beyond secondary treatment) as a management alternative to surface-water disposal of wastewater. Water Conserv II, a water reuse project developed jointly by Orange County and the City of Orlando, began operation in December 1986. In 1995, the Water Conserv II facility distributed approximately 28 Mgal/d of reclaimed water for discharge to rapid-infiltration basins (RIBs) and for use as agricultural irrigation. The Reedy Creek Improvement District (RCID) began operation of RIBs in September 1990, and in 1995 these RIBs received approximately 6.7 Mgal/d of reclaimed water. Analyses of existing data and data collected during the course of this study were combined with ground-water flow modeling and particle-tracking analyses to develop a process-oriented evaluation of the regional effects of reclaimed water applied by Water Conserv II and the RCID RIBs on the hydrology of west Orange and southeast Lake Counties. The ground-water flow system beneath the study area is a multi-aquifer system that consists of a thick sequence of highly permeable carbonate rocks overlain by unconsolidated sediments. The hydrogeologic units are the unconfined surficial aquifer system, the intermediate confining unit, and the confined Floridan aquifer system, which consists of two major permeable zones, the Upper and Lower Floridan aquifers, separated by the less permeable middle semiconfining unit. Flow in the surficial aquifer system is dominated regionally by diffuse downward leakage to the Floridan aquifer system and is affected locally by lateral flow systems produced by streams, lakes, and spatial variations in recharge. Ground water generally flows laterally through the Upper Floridan aquifer aquifer to the north and east. Many of the lakes in the study area are landlocked because the mantled karst environment precludes a well developed network of surface-water drainage. The USGS three-dimensional ground-water flow model MODFLOW was used to simulate ground-water flow in the surficial and Floridan aquifer systems. A steady-state calibration to average 1995 conditions was performed by using a parameter estimation program to vary values of surficial aquifer system hydraulic conductivity, intermediate confining unit leakance, and Upper Floridan aquifer transmissivity. The calibrated model generally produced simulated water levels in close agreement with measured water levels and was used to simulate the hydrologic effects of reclaimed-water application under current (1995) and proposed future conditions. In 1995, increases of up to about 40 ft in the water table and less than 5 ft in the Upper Floridan aquifer potentiometric surface had occurred as a result of reclaimed-water application. The largest increases were under RIB sites. An average traveltime of 10 years at Water Conserv II and 7 years at the RCID RIBs was required for reclaimed water to move from the water table to the top of the Upper Floridan aquifer. Approximately 67 percent of the reclaimed water applied at the RCID RIB site recharged the Floridan aquifer system, whereas 33 percent discharged from the surficial aquifer system to surface-water features; 99 percent of the reclaimed water applied at Water Conserv II recharged the Floridan aquifer system, whereas only 1 percent discharged from the surficial aquifer system to surface-water features. The majority of reclaimed water applied at both facilities probably will ultimately discharge from the Floridan aquifer system outside the model boundaries. Proposed future conditions were assumed to consist of an additional 11.7 Mgal/d of reclaimed water distributed by the Water Conserv II and RCID facilities. Increases of up to about 20 ft in the water table and 2 ft in the potentiometric surface of the Upper Floridan aquifer were simulated. The directions of reclaimed water movement through the ground-water system generally were similar to those under 1995 conditions. However, the greater reclaimed-water application rate at the RCID RIBs caused approximately half of the RCID reclaimed water to discharge to surface-water features and half to recharge the Floridan aquifer system.

Florida

Estimation of shallow ground-water recharge in the Great Lakes basin

This report presents the results of the first known integrated study of long-term average ground-water recharge to shallow aquifers (generally less than 100 feet deep) in the United States and Canada for the Great Lakes, upper St. Lawrence, and Ottawa River Basins. The approach used was consistent throughout the study area and allows direct comparison of recharge rates in disparate parts of the study area. Estimates of recharge are based on base-flow estimates for streams throughout the Great Lakes Basin and the assumption that base flow in a given stream is equal to the amount of shallow ground-water recharge to the surrounding watershed, minus losses to evapotranspiration. Base-flow estimates were developed throughout the study area using a single model based on an empirical relation between measured base-flow characteristics at streamflow-gaging stations and the surficial-geologic materials, which consist of bedrock, coarse-textured deposits, fine-textured deposits, till, and organic matter, in the surrounding surface-water watershed. Model calibration was performed using base-flow index (BFI) estimates for 959 stations in the U.S. and Canada using a combined 28,784 years of daily streamflow record determined using the hydrograph-separation software program PART. Results are presented for watersheds represented by 8-digit hydrologic unit code (HUC, U.S.) and tertiary (Canada) watersheds. Recharge values were lowest (1.6-4.0 inches/year) in the eastern Lower Peninsula of Michigan; southwest of Green Bay, Wisconsin; in northwestern Ohio; and immediately south of the St. Lawrence River northeast of Lake Ontario. Recharge values were highest (12-16.8 inches/year) in snow shadow areas east and southeast of each Great Lake. Further studies of deep aquifer recharge and the temporal variability of recharge would be needed to gain a more complete understanding of ground-water recharge in the Great Lakes Basin.

Great Lakes Basin

Stream-aquifer interactions in the Straight River area, Becker and Hubbard counties, Minnesota

The Straight River, in north-central Minnesota, is a trout stream having cold, clear water. The 75-square-mile Straight River watershed contributes flow to the stream. The watershed is underlain by highly transmissive surficial and confined-drift aquifers. Ground-water discharge from these aquifers sustains flow in the Straight River, and the cold water supports a population of trout. Water withdrawals from these aquifers are increasing in response to changes in land use from dry-land to irrigated fanning. Degradation of the stream's habitat for trout could result from the following: a decrease in ground-water discharge to the stream caused by ground-water withdrawals for irrigation, an increase in ground-water temperature resulting from percolation of irrigated water to the ground-water system, and introduction of agricultural chemicals to the stream through ground-water flow or runoff. Physical data indicate a hydraulic connection between the stream and the surficial aquifer. Discharge of the Straight River increases from about 25 cubic feet per second at the outfall from a reservoir near the headwaters to about 51 cubic feet per second near the mouth. The rate of streamflow gain during summer decreases downstream, possibly as a result of ground-water withdrawal for irrigation. The water table and potentiometric surface of the uppermost confined-drift aquifer generally slope to the southeast and locally toward rivers and lakes; gradients decline to about 5 feet per mile from spring to summer. Daily fluctuations of stream temperature are as great as 15 degrees Celsius during the summer, primarily in response to changes in air temperature. Ground-water discharge to the Straight River decreases stream temperature during the summer. Results of simulations from a stream-temperature model indicate that daily changes in stream temperature are strongly influenced by solar radiation, wind speed, stream depth, and ground-water inflow. Results of simulations from ground-water-flow and stream-temperature models developed for the investigation indicate a significant decrease in ground-water flow could result from ground-water withdrawal at rates similar to those measured during 1988. This reduction in discharge to the stream could result in an increase in stream temperature of 0.5 to 1.5 degrees Celsius. Nitrate concentrations in shallow wells screened at the water table, in some areas, are locally greater than the limit set by the Minnesota Pollution Control Agency. Nitrate concentrations in water from deeper wells and in the stream are low, generally less than 1.0 milligram per liter.

Minnesota

Potential hydrologic effects of peat mining in the Red Lake Peatlands, north-central Minnesota— A project plan

Peat is being considered for fuel in Minnesota. This study will investigate the potential effects of large-scale surface mining of peat on the hydrology and water quality of Upper Red Lake and the Tamarac River. The major aspects of the study are the characterization of the surface-water and groundwater hydrology and water quality, including the trace-metal content of the peat. Data will be collected to construct two- and three-dimensional digital models to simulate the movement of ground water and its relation to surface water in the peatlands, streams, and lakes. After the model is calibrated with field data, it will be used to evaluate the effect of mining peat on the hydrology and water quality of the Upper Red Lake and Tamarac River.

Minnesota

Simulated effects of proposed ground-water pumping in 17 basins of east-central and southern Nevada

The Las Vegas Valley Water District filed 146 applications in 1989 to pump about 180,800 acre- ft/yr in 17 basins for use in Las Vegas Valley. A previously constructed, two-layer computer model of the carbonate-rock province area was configured to simulate transient conditions and used to develop first approximations of the possible effects of these withdrawals. Simulations were made using the phased pumping schedule proposed by the water district. Ground-water-level declines of several hundred feet could ultimately develop in the basins scheduled to supply most of the pumped ground water. Simulated declines in the carbonate-rock aquifer were somewhat larger than simulated declines in the overlying basin-fill deposits. Decreases in simulated regional spring flow were shown in several cells including those representing the Muddy River Springs, Hiko-Crystal-Ash spring area, and the Ash Meadows spring area. Model simulations show flow decreases of about 11 percent, 14 percent, and 2 percent, respectively, at these springs after almost 100 years of pumping. Simulated evapotranspiration also decreased in many basins, with the largest decreases occurring in the basins where ground-water withdrawals were greatest. These basins include Railroad, Spring, and Snake Valleys. The largest decrease in simulated evapotranspiration occurred in Railroad Valley, 64 percent after almost 100 years of pumpage. Model-sensitivity tests indicate that long-term results were relatively insensitive to variations in values used for aquifer storage. The adequacy of the model to simulate the effects of this proposed pumping will remain untested until actual pumping stresses have been in place long enough to cause measurable effects within the system.

Water-Resources Investigations Report

Evaluation of drawdown and sources of water in the Mississippi River alluvium caused by hypothetical pumping, Muscatine, Iowa

A study was conducted to evaluate drawdown and volumetric changes in sources of water in the Mississippi River alluvium caused by hypothetical pumping. A steady-state, ground-water flow model was constructed for a previous study to simulate February 1993 hydrologic conditions, which were assumed to be an acceptable estimate of the ground-water system at equilibrium. The flow model was modified for this study to simulate six hypothetical pumping scenarios: five pumping scenarios to simulate hypothetical pumping at five pumping scenario sites, and a total pumping scenario to simulate cumulative hypothetical pumping from the five pumping scenario sites. The evaluation of drawdown for the six hypothetical pumping scenarios indicates that hypothetical pumping causes simulated drawdown that varies from about 10 ft to greater than 50 ft relative to February 1993 conditions at the five hypothetical pumping scenario sites. The simulated drawdown is less than half of the estimated saturated thickness of the alluvium during February 1993 at these sites. The primary sources of water (inflows) to the alluvium needed to balance the increased ground-water withdrawals (outflows) caused by hypothetical pumping are a combination of increased river leakage and decreased leakage to Muscatine Slough. Compared to February 1993 conditions, larger inflow rates occur as river leakage from the Mississippi River for the six hypothetical pumping scenarios. However, smaller outflow rates to Muscatine Slough compared to February 1993 conditions indicate that an important source of water for hypothetical pumping is ground-water discharge that would have become streamflow in the slough. Increased pumping at the hypothetical pumping scenario sites could affect long-term water quality and hydrology in the study area. The greater amounts of river leakage might affect overall ground-water quality in the alluvium. The lesser amounts of ground water being discharged to streamflow could have a long-term impact on the hydrology of the slough and adjacent wetland areas. The simplified steady-state flow model does not account for dynamic (transient) conditions (natural or development-related). The steady-state model does not indicate time needed to reach new equilibrium conditions. Attaining equilibrium might take many years and is complicated by varying climatic and hydrologic conditions; noncontinuous pumping and pumping that is cycled among well fields; and changing and seasonally varying irrigation pumpage.

Illinois, Iowa

Simulation of nutrient and sediment concentrations and loads in the Delaware inland bays watershed: Extension of the hydrologic and water-quality model to ungaged segments

Rapid population increases, agriculture, and industrial practices have been identified as important sources of excessive nutrients and sediments in the Delaware Inland Bays watershed. The amount and effect of excessive nutrients and sediments in the Inland Bays watershed have been well documented by the Delaware Geological Survey, the Delaware Department of Natural Resources and Environmental Control, the U.S. Environmental Protection Agency’s National Estuary Program, the Delaware Center for Inland Bays, the University of Delaware, and other agencies. This documentation and data previously were used to develop a hydrologic and water-quality model of the Delaware Inland Bays watershed to simulate nutrients and sediment concentrations and loads, and to calibrate the model by comparing concentrations and streamflow data at six stations in the watershed over a limited period of time (October 1998 through April 2000). Although the model predictions of nutrient and sediment concentrations for the calibrated segments were fairly accurate, the predictions for the 28 ungaged segments located near tidal areas, where stream data were not available, were above the range of values measured in the area. The cooperative study established in 2000 by the Delaware Department of Natural Resources and Environmental Control, the Delaware Geological Survey, and the U.S. Geological Survey was extended to evaluate the model predictions in ungaged segments and to ensure that the model, developed as a planning and management tool, could accurately predict nutrient and sediment concentrations within the measured range of values in the area. The evaluation of the predictions was limited to the period of calibration (1999) of the 2003 model. To develop estimates on ungaged watersheds, parameter values from calibrated segments are transferred to the ungaged segments; however, accurate predictions are unlikely where parameter transference is subject to error. The unexpected nutrient and sediment concentrations simulated with the 2003 model were likely the result of inappropriate criteria for the transference of parameter values. From a model-simulation perspective, it is a common practice to transfer parameter values based on the similarity of soils or the similarity of land-use proportions between segments. For the Inland Bays model, the similarity of soils between segments was used as the basis to transfer parameter values. An alternative approach, which is documented in this report, is based on the similarity of the spatial distribution of the land use between segments and the similarity of land-use proportions, as these can be important factors for the transference of parameter values in lumped models. Previous work determined that the difference in the variation of runoff due to various spatial distributions of land use within a watershed can cause substantialloss of accuracy in the model predictions. The incorporation of the spatial distribution of land use to transfer parameter values from calibrated to uncalibrated segments provided more consistent and rational predictions of flow, especially during the summer, and consequently, predictions of lower nutrient concentrations during the same period. For the segments where the similarity of spatial distribution of land use was not clearly established with a calibrated segment, the similarity of the location of the most impervious areas was also used as a criterion for the transference of parameter values. The model predictions from the 28 ungaged segments were verified through comparison with measured in-stream concentrations from local and nearby streams provided by the Delaware Department of Natural Resources and Environmental Control. Model results indicated that the predicted edge-of-stream total suspended solids loads in the Inland Bays watershed were low in comparison to loads reported for the Eastern Shore of Maryland from the Chesapeake Bay watershed model. The flatness of the terrain and the low annual surface runoff are important factors in determining the amount of detached sediment from the land that is delivered to streams. The highest predicted total suspended solids loads were found in the southern part of the watershed, where the values are associated with high total streamflow and a high surface-runoff component, and related to soil and aquifer permeability and land use. Nutrient loads from model segments in the southern part of the Inland Bays watershed were also higher than those measured in the northern part of the basin, due to relatively high runoff and the substantial amount of available organic fertilizer (animal waste) that results in over-application of organic fertilizer to crops. Time series of simulated hourly concentrations indicated a seasonal pattern in the simulated base flow for total nitrogen, with the lowest values occurring during the summer and the highest values during the winter months. Total phosphorus and total-suspended-solids concentrations were less seasonal and were more storm-dependent; in general, base-flow concentrations of total phosphorus and total suspended solids were low. During storm events, the total nitrogen concentrations tended to be diluted and total phosphorus concentrations tended to rise sharply. Nitrogen was transported mainly in the aqueous phase and largely through ground water, whereas phosphorus was strongly associated with sediment, which washes off during rainfall events.

Delaware

Generalized estimates from streamflow data of annual and seasonal ground-water-recharge rates for drainage basins in New Hampshire

This report presents regression equations to estimate generalized annual and seasonal ground-water-recharge rates in drainage basins in New Hampshire. The ultimate source of water for a ground-water withdrawal is aquifer recharge from a combination of precipitation on the aquifer, ground-water flow from upland basin areas, and infiltration from streambeds to the aquifer. An assessment of ground-water availability in a basin requires that recharge rates be estimated under `normal' conditions and under assumed drought conditions. Recharge equations were developed by analyzing streamflow, basin characteristics, and precipitation at 55 unregulated continuous record stream-gaging stations in New Hampshire and in adjacent states. In the initial step, streamflow records were analyzed to estimate a series of annual and seasonal ground-water-recharge components of streamflow in each drainage basin evaluated in this study. Regression equations were then developed relating the series of annual and seasonal ground-water-recharge values to the corresponding series of annual and seasonal precipitation values as determined at the centroid of each drainage basin. This resulted in one equation for each of the 55 basins for each of the four seasonal periods and the annual period, or a total of 275 regression equations. Average annual and seasonal precipitation data for 1961-90 were then used to compute a set of normalized ground-water-recharge values that reflected the long-term average annual and seasonal variations (normalized) and mean recharge characteristics of each drainage basin. Ordinary-least-squares regression was applied in the process of selecting 10 out of 93 possible basin and climatic characteristics for further testing in the development of the equations for computing the generalized estimate of annual and seasonal ground-water recharge based on the set of normalized recharge values. Generalized-least-squares regression was used for the final parameter estimation and error evaluation. The following basin and climatic characteristics were found to be statistically significant predictors for at least one of the dependent variables: average annual, summer, and spring precipitation as determined at U.S. Geological Survey stream-gaging stations; average annual basin-centroid precipitation; average mean annual basin temperature; average minimum winter basin temperature; percent coniferous forest in a basin; percent mixed coniferous and deciduous forest in a basin; average fall basin-centroid precipitation; and average annual snowcover. These 10 basin and climatic characteristics were selected because they were statistically significant based on several statistical parameters that evaluated which combination of characteristics contributed the most to the predictive accuracy of the regression-equation models. A geographic information system is required to measure the values of the predictor variables for the equations developed in the study. The average annual normalized ground-water recharge was 21.0 in. This value was determined by generalized-least-squares (GLS) regression analysis for all of the basins used in the normalized ground-water recharge analysis for precipitation from 1961-90. The average winter (January 1-March 15) ground-water recharge was 4.3 in., average spring (March 16-May 31) ground-water recharge was 9.0 in., average summer (June 1-October 31) ground-water recharge was 4.0 in., and average fall (November 1-December 31) ground-water recharge was 3.6 in. Normalized ground-water recharge ranged annually from 12.3 to 31.8 in., for winter from 2.30 to 7.82 in., for spring from 5.16 to 13.7 in., for summer from 1.45 to 10.2 in., and for fall from 2.21 to 6.06 in.

New Hampshire

Digital model of ground-water flow in the Piceance Basin, Rio Blanco and Garfield Counties, Colorado

The digital model used to simulate ground-water flow in the aquifer system in the basin drained by Piceance and Yellow Creeks in northwestern Colorado is described in detail. The model is quasi three-dimensional in that it simulates ground-water flow in a multiaquifer system by assuming horizontal flow in the aquifers and vertical flow through the confining layers separating the aquifers. The model uses the iterative alternating-direction implicit procedure to solve the finite-difference flow equations. The digital model is documented by a program listing and flow charts. Data used in the model and sample output are presented to document the simulation of steady-state flow in the aquifer system. The variables used in the computer program and program options are discussed in detail. (Woodard-USGS)

Water-Resources Investigations Report

A coupled surface-water and ground-water flow model (MODBRANCH) for simulation of stream-aquifer interaction

Ground-water and surface-water flow models traditionally have been developed separately, with interaction between subsurface flow and streamflow either not simulated at all or accounted for by simple formulations. In areas with dynamic and hydraulically well-connected ground-water and surface-water systems, stream-aquifer interaction should be simulated using deterministic responses of both systems coupled at the stream-aquifer interface. Accordingly, a new coupled ground-water and surface-water model was developed by combining the U.S. Geological Survey models MODFLOW and BRANCH; the interfacing code is referred to as MODBRANCH. MODFLOW is the widely used modular three-dimensional, finite-difference ground-water model, and BRANCH is a one-dimensional numerical model commonly used to simulate unsteady flow in open- channel networks. MODFLOW was originally written with the River package, which calculates leakage between the aquifer and stream, assuming that the stream's stage remains constant during one model stress period. A simple streamflow routing model has been added to MODFLOW, but is limited to steady flow in rectangular, prismatic channels. To overcome these limitations, the BRANCH model, which simulates unsteady, nonuniform flow by solving the St. Venant equations, was restructured and incorporated into MODFLOW. Terms that describe leakage between stream and aquifer as a function of streambed conductance and differences in aquifer and stream stage were added to the continuity equation in BRANCH. Thus, leakage between the aquifer and stream can be calculated separately in each model, or leakages calculated in BRANCH can be used in MODFLOW. Total mass in the coupled models is accounted for and conserved. The BRANCH model calculates new stream stages for each time interval in a transient simulation based on upstream boundary conditions, stream properties, and initial estimates of aquifer heads. Next, aquifer heads are calculated in MODFLOW based on stream stages calculated by BRANCH, aquifer properties, and stresses. This process is repeated until convergence criteria are met for head and stage. Because time steps used in ground-water modeling can be much longer than time intervals used in surface- water simulations, provision has been made for handling multiple BRANCH time intervals within one MODFLOW time step. An option was also added to BRANCH to allow the simulation of channel drying and rewetting. Testing of the coupled model was verified by using data from previous studies; by comparing results with output from a simpler, four-point implicit, open-channel flow model linked with MODFLOW; and by comparison to field studies of L-31N canal in southern Florida.

Techniques of Water-Resources Investigations

Delineation of the Troy Bedrock Valley and particle-tracking analysis of ground-water flow underlying Belvidere, Illinois

The U.S. Geological Survey, U.S. Environmental Protection Agency, and Illinois Environmental Protection Agency began a study of the hydrogeology, flow system, and distribution of contaminants in the aquifers underlying Belvidere, Ill., and vicinity in 1992. As part of the study, the ancestral Troy Bedrock Valley, located about 1.5 miles west of Belvidere, was identified as an important part of the ground-water-flow system. In the deepest parts of the valley, the basal Glenwood confining unit may be absent; thick deposits of sand and gravel that infill part of the valley may directly overlie the sandstone St. Peter aquifer, a regionally important source of water for public supply. With few deep wells open to the St. Peter aquifer present in the valley to provide necessary geologic information, tritium and other water-chemistry data were collected from eight wells to possibly delineate areas where the confining unit may be absent; the data also provide baseline water-quality information for an area expecting changes in land use and increases in water withdrawal. Also as part of the study, particle-tracking analysis was done using an available flow model to (1) identify possible discharge locations of ground water and contaminants and (2) delineate areas contributing recharge to the Belvidere municipal wells. This report presents and interprets water-chemistry data collected during December 2000 and presents results of particle-tracking analysis. Ground water in samples from two of four wells open to the St. Peter aquifer appears to have recharged after 1954, suggesting that the Glenwood confining unit may be absent near the wells. Other hydrogeologic and water-chemistry data, however, were inconclusive or contradictory. Concentrations of iron, manganese, and lead exceeded maximum contaminant levels in five or less samples, but materials associated with the water-distribution systems appear to contribute to the elevated concentrations above natural levels. Particle-tracking analysis indicates that most ground-water flow beneath possible contaminant-source areas discharges from the glacial drift aquifer to the Kishwaukee River. Most of the source areas are in or near Belvidere and are within 1,500 feet of the river. The analysis also indicates the possibility that in parts of the study area, some ground water does not discharge to the river, but flows beneath the Kishwaukee River in the underlying carbonate Galena-Platteville aquifer. Ground water that discharges to the one municipal well open to the glacial drift aquifer is estimated to travel over 1 mile in less than 25 years. Simulated residence (travel) times of ground water from the base of the glacial drift aquifer to the six municipal wells open, in part, to the Galena-Platteville aquifer, are estimated at less than about 40 years. Because fractures in this aquifer are unaccounted for in the flow model, actual areas contributing recharge are likely larger and travel times faster than those simulated for most of the municipal wells. Tritium data indicate that, in general, travel times from the land surface to the deepest parts of the Galena-Platteville aquifer are less than 46 years. Methyl tertiary-butyl ether data indicate that travel times to the upper part of the aquifer may be less than 16 years. The water-quality-based estimates of travel time generally are consistent with the estimates from particle-tracking analysis.

Illinois

Effects of Withdrawals on Ground-Water Levels in Southern Maryland and the Adjacent Eastern Shore, 1980-2005

Ground water is the primary source of water supply in most areas of Maryland?s Atlantic Coastal Plain, including Southern Maryland. The counties in this area are experiencing some of the most rapid growth and development in the State, resulting in an increased demand for ground-water production. The cooperative, basic water-data program of the U.S. Geological Survey and the Maryland Geological Survey has collected long-term observations of ground-water levels in Southern Maryland and parts of the Eastern Shore for many decades. Additional water-level observations were made by both agencies beginning in the 1970s, under the Power Plant Research Program of the Maryland Department of Natural Resources. These long-term water levels commonly show significant declines over several decades, which are attributed to ground-water withdrawals. Ground-water-level trends since 1980 in major Coastal Plain aquifers such as the Piney Point-Nanjemoy, Aquia, Magothy, upper Patapsco, lower Patapsco, and Patuxent were compared to water use and withdrawal data. Potentiometric surface maps show that most of the declines in ground-water levels can be directly related to effects from major pumping centers. There is also evidence that deep drawdowns in some pumped aquifers may be causing declines in adjacent, unpumped aquifers. Water-level hydrographs of many wells in Southern Maryland show linear declines in levels year after year, instead of the gradual leveling-off that would be expected as the aquifers equilibrate with pumping. A continual increase in the volumes of water being withdrawn from the aquifers is one explanation for why they are not reaching equilibrium. Although reported ground-water production in Southern Maryland has increased somewhat over the past several decades, the reported increases are often not large enough to account for the observed water-level declines. Numerical modeling simulations indicate that a steady, annual increase in the number of small wells could account for the observed aquifer behavior. Such wells, being pumped at rates below the minimum legal reporting threshold of 10,000 gallons per day, might be the source of the additional withdrawals. More detailed water-use data, especially from domestic wells, central-pivot irrigation wells, and other small users not currently reporting withdrawals to the State, may help to determine the cause of the aquifer declines.

Scientific Investigations Report

Water quality and the effects of changes in phosphorus loading, Red Cedar Lakes, Barron and Washburn Counties, Wisconsin

The Red Cedar Lakes consist of three mainstem lakes (Balsam, Hemlock and Red Cedar) on the Red Cedar River in Barron and Washburn Counties, Wisconsin. These lakes are productive because of high phosphorus loading, and classified as mesotrophic to eutrophic. Because of concerns that the water quality of these lakes was degrading, three cooperative studies were conducted by the U.S. Geological Survey between 1993 and 2003. As part of these studies, water quality in the lakes was documented in 1993?94, 1996?97, and 2000?01, and water and phosphorus budgets were constructed for water year 2001. Historical water-quality data indicated that the lakes have changed little since the late 1980s. A detailed phosphorus budget indicated that most of the 14,100 pounds of phosphorus input to the lakes during 2001 came from the upstream lakes, streams draining relatively undeveloped land upstream of Hemlock Lake, and ground water. Simulation results from two water-quality models (BATHTUB and WiLMS) indicated that about a 50-percent reduction in phosphorus loading from that measured in 2001 is required for all three lakes to be classified as mesotrophic; therefore, appreciable improvements in the water quality would require improvements in the water quality of the upstream lakes. Although the water quality of the lakes has not changed appreciably in recent years and major improvements in water quality are unlikely without major improvements to upstream lakes, continued efforts to protect the susceptible watershed are necessary to maintain the current level of water quality.

Wisconsin

Hydrogeology and ground-water/surface water interactions in the Des Moines River valley, southwestern Minnesota, 1997-2001

Increased water demand in and around Windom led the U.S. Geological Survey, in cooperation with the Minnesota Department of Natural Resources, local water suppliers, and Cottonwood County, to study the hydrology of aquifers in the Des Moines River Valley near Windom. The study area is the watershed of a 30-kilometer (19-mile) reach of the Des Moines River upstream from Windom. Based on stratigraphic analysis, two hydrologically and genetically separate surficial aquifers underlie the study area. The Windom aquifer has a saturated thickness of 34 meters (111 feet), and the Des Moines aquifer has a saturated thickness of 33 meters (108 ft). The surficial aquifers are relatively isolated from deeper aquifers by till, but some leakage probably occurs. Recharge to the aquifers is from areal recharge, from Cottonwood Lake, and from edge recharge. Pumping at the Windom well field induces substantial amounts of Cottonwood Lake water into the aquifer. During this study, the water level in a well located between two Red Rock wells and the river was lower than the river level during two periods. During those periods, water in the Des Moines River had the potential to recharge the aquifer. Discharge from the aquifers is primarily to municipal wells, the Des Moines River, and other surface waters. Most of the ground-water samples collected in the study area consisted of calcium-magnesium bicarbonate waters. Corn and soybean herbicides and their degradates were detected at low concentrations in 14 of 27 ground-water samples and in all 3 river samples. Metolachlor ethane sulfonic acid was the most commonly detected compound and also was detected at the highest concentrations. Nutrient concentrations in ground-water samples were skewed low with high outliers, and nutrient concentrations in river samples generally were less than analytical reporting limits. Nearly all recharge to the aquifer in the ground-water simulation was from edge recharge (80 percent). Calibrated net areal recharge ranged from 17 to 30 percent of the average annual precipitation. Isotopic composition of ground water and Cottonwood Lake water indicated about one-half of the water withdrawn from the Windom aquifer is from Cottonwood Lake. Scenarios tested with the calibrated model involved increased ground-water withdrawals and changes in recharge to simulate drier or wetter weather conditions. Doubling the withdrawals from all wells in the model had a small effect except in the Windom well-field area. Maximum head declines in the Red Rock well field and the Jeffers city well were less than 40 centimeters (15 inches). In the Windom well field, the maximum head decline was 11 meters (36 feet). The Windom well field does not induce recharge from the Des Moines River. The addition of a new well that pumped 2,000 cubic meters per day (0.44 million gallons per day) in the Augusta Lake Valley area caused a 0.83-meter-deep (2.72-foot-deep) cone of depression that extended to the valley walls. The drought scenario and the high-precipitation scenario resulted in head changes in the northern part of the Augusta Lake Valley area, in the southwestern part of the Red Rock area, and near the valley edges. Long-term withdrawals of water for public supplies may cause a net decrease in ground-water discharge to surface water. Water that does not evaporate, or that is not exported, is discharged to the Des Moines River but with changed water quality. Because ground-water and surface-water qualities in the study area are similar, the ground-water discharge probably has little effect on river water quality.

Minnesota

Hydrologic conditions and assessment of water resources in the Turkey Creek watershed, Jefferson County, Colorado, 1998-2001

The 47.2-square-mile Turkey Creek watershed, in Jefferson County southwest of Denver, Colorado, is relatively steep with about 4,000 feet of relief and is in an area of fractured crystalline rocks of Precambrian age. Water needs for about 4,900 households in the watershed are served by domestic wells and individual sewage-disposal systems. Hydrologic conditions are described on the basis of contemporary hydrologic and geologic data collected in the watershed from early spring 1998 through September 2001. The water resources are assessed using discrete fracture-network modeling to estimate porosity and a physically based, distributed-parameter watershed runoff model to develop estimates of water-balance terms. A variety of climatologic and hydrologic data were collected. Direct measurements of evapotranspiration indicate that a large amount (3 calendar-year mean of 82.9 percent) of precipitation is returned to the atmosphere. Surface-water records from January 1, 1999, through September 30, 2001, indicate that about 9 percent of precipitation leaves the watershed as streamflow in a seasonal pattern, with highest streamflows generally occurring in spring related to snowmelt and precipitation. Although conditions vary considerably within the watershed, overall watershed streamflow, based on several records collected during the 1940's, 1950's, 1980', and 1990's near the downstream part of watershed, can be as high as about 200 cubic feet per second on a daily basis during spring. Streamflow typically recedes to about 1 cubic foot per second or less during rainless periods and is rarely zero. Ground-water level data indicate a seasonal pattern similar to that of surface water in which water levels are highest, rising tens of feet in some locations, in the spring and then receding during rainless periods at relatively constant rates until recharged. Synoptic measurements of water levels in 131 mostly domestic wells in fall of 2001 indicate a water-table surface that conforms to topography. Analyses of reported well-construction records indicate a median reported well yield of 4 gallons per minute and a spatial distribution for reported well yield that has relatively uniform conditions of small-scale variability. Results from quarterly samples collected in water year 1999 at about 112 wells and 22 streams indicate relatively concentrated calcium-bicarbonate to calcium-chloride type water that has a higher concentration of chloride than would be expected on the basis of chloride content in precipitation and evapotranspiration rates. Comparison of the 1999 data to similar data collected in the 1970's indicates that concentrations for many constituents appear to have increased. Reconnaissance sampling in the fall of 2000 indicates that most ground water in the watershed was recharged recently, although some ground water was recharged more than 50 years ago. Additional reconnaissance sampling in the spring and fall of 2001 identified some compounds indicative of human wastewater in ground water and surface water. Outcrop fracture measurements were used to estimate potential porosities in three rock groups (metamorphic, intrusive, and fault zone) that have distinct fracture characteristics. The characterization, assuming a uniform aperture size of 100 microns, indicates very low potential fracture porosities, on the order of hundredths of a percent for metamorphic and intrusive rocks and up to about 2 percent for fault-zone rocks. A fourth rock group, Pikes Peak Granite, was defined on the basis of weathering characteristics. Short-term continuous and synoptic measurements of streamflow were used to describe base-flow characteristics in areas of the watershed underlain by each of the four rock groups and are the basis for characterization of base flow in a physically based, distributed-parameter watershed model. The watershed model, the Precipitation-Runoff Modeling System (PRMS), was used to characterize hydrologic conditions

Water-Resources Investigations Report

Water quality and the effects of changes in phosphorus loading to Muskellunge Lake, Vilas County, Wisconsin

Muskellunge Lake is a productive, eutrophic lake because of high nutrient loading. Historical data indicate that water quality has only slightly degraded since the early 1970s, possibly because of phosphorus input from effluent from septic systems. A detailed phosphorus budget for the lake indicated that most of the phosphorus comes from natural sources?ground water and surface water flowing through relatively undeveloped areas surrounding the lake. Modeling results indicated that the natural input of phosphorus was sufficient to maintain the lake's eutrophic condition. Analysis of sediment cores confirmed that only small changes in nutrient and algal concentrations have occurred over the past 100 years; however, the analysis indicated that the macrophyte community has increased over this time period. The aeration system, installed to alleviate winter anoxia, maintains aerobic conditions throughout the main bays of the lake.

Wisconsin