Search USGS⌕ Search

SEARCH · Search USGS

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 721 records · Page 40Linked to original sources

Hydrology of Cache Valley, Cache County, Utah, and adjacent part of Idaho, with emphasis on simulation of ground-water flow

A hydrologic investigation of Cache Valley was done to better understand the ground-water system in unconsolidated basin-fill deposits and the interaction between ground water and surface water. Ground-water recharge occurs by infiltration of precipitation and unconsumed irrigation water, seepage from canals and streams, and subsurface inflow from adjacent consolidated rock and adjacent unconsolidated basin-fill deposit ground-water systems. Ground-water discharge occurs as seepage to streams and reservoirs, spring discharge, evapotranspiration, and withdrawal from wells. Water levels declined during 1984-90. Less-than-average precipitation during 1987-90 and increased pumping from irrigation and public-supply wells contributed to the declines. A ground-water-flow model was used to simulate flow in the unconsolidated basin-fill deposits. Data primarily from 1969 were used to calibrate the model to steady-state conditions. Transient-state calibration was done by simulating ground-water conditions on a yearly basis for 1982-90. A hypothetical simulation in which the dry conditions of 1990 were continued for 5 years projected an average lO-foot water-level decline between Richmond and Hyrum. When increased pumpage was simulated by adding three well fields, each pumping 10 cubic feet per second, in the Logan, Smithfield, and College Ward areas, water-level declines greater than 10 feet were projected in most of the southeastern part of the valley and discharge from springs and seepage to streams and reservoirs decreased.

Idaho, Utah↗

Assessment of channel changes, model of historical floods, and effects of backwater on flood stage, and flood mitigation alternatives for the Wichita River at Wichita Falls, Texas

In cooperation with the City of Wichita Falls, the U.S. Geological Survey assessed channel changes on the Wichita River at Wichita Falls, Texas, and modeled historical floods to investigate possible causes and potential mitigation alternatives to higher flood stages in recent (2007 and 2008) floods. Extreme flooding occurred on the Wichita River on June 30, 2007, inundating 167 homes in Wichita Falls. Although a record flood stage was reached in June 2007, the peak discharge was much less than some historical floods at Wichita Falls. Streamflow and stage data from two gages on the Wichita River and one on Holliday Creek were used to assess the interaction of the two streams. Changes in the Wichita River channel were evaluated using historical aerial and ground photography, comparison of recent and historical cross sections, and comparison of channel roughness coefficients with those from earlier studies. The floods of 2007 and 2008 were modeled using a one-dimensional step-backwater model. Calibrated channel roughness was larger for the 2007 flood compared to the 2008 flood, and the 2007 flood peaked about 4 feet higher than the 2008 flood. Calibration of the 1941 flood yielded a channel roughness coefficient (Manning's n ) of 0.030, which represents a fairly clean natural channel. The step-backwater model was also used to evaluate the following potential mitigation alternatives: (1) increasing the capacity of the bypass channel near River Road in Wichita Falls, Texas; (2) removal of obstructions near the Scott Avenue and Martin Luther King Junior Boulevard bridges in Wichita Falls, Texas; (3) widening of aggraded channel banks in the reach between Martin Luther King Junior Boulevard and River Road; and (4) reducing channel bank and overbank roughness. Reductions in water-surface elevations ranged from 0.1 foot to as much as 3.0 feet for the different mitigation alternatives. The effects of implementing a combination of different flood-mitigation alternatives were not investigated.

Texas↗

Hydrogeology and simulated effects of ground-water withdrawals from the Floridan aquifer system in Lake County and in the Ocala National Forest and vicinity, north-central Florida

The hydrogeology of Lake County and the Ocala National Forest in north-central Florida was evaluated (1995-2000), and a ground-water flow model was developed and calibrated to simulate the effects of both present day and future ground-water withdrawals in these areas and the surrounding vicinity. A predictive model simulation was performed to determine the effects of projected 2020 ground-water withdrawals on the water levels and flows in the surficial and Floridan aquifer systems. The principal water-bearing units in Lake County and the Ocala National Forest are the surficial and Floridan aquifer systems. The two aquifer systems generally are separated by the intermediate confining unit, which contains beds of lower permeability sediments that confine the water in the Florida aquifer system. The Floridan aquifer system has two major water-bearing zones (the Upper Floridan aquifer and the Lower Floridan aquifer), which generally are separated by one or two less-permeable confining units. The Floridan aquifer system is the major source of ground water in the study area. In 1998, ground-water withdrawals totaled about 115 million gallons per day in Lake County and 5.7 million gallons per day in the Ocala National Forest. Of the total ground water pumped in Lake County in 1998, nearly 50 percent was used for agricultural purposes, more than 40 percent for municipal, domestic, and recreation supplies, and less than 10 percent for commercial and industrial purposes. Fluctuations of lake stages, surficial and Floridan aquifer system water levels, and Upper Floridan aquifer springflows in the study area are highly related to cycles and distribution of rainfall. Long-term hydrographs for 9 lakes, 8 surficial aquifer system and Upper Floridan aquifer wells, and 23 Upper Floridan aquifer springs show the most significant increases in water levels and springflows following consecutive years with above-average rainfall, and significant decreases following consecutive years with below-average rainfall. Long-term (1940-2000) hydrographs of lake and ground-water levels and springflow show a slight downward trend; however, after the early 1960's, this downward trend generally is more pronounced, which corresponds with accumulating rainfall deficits and increased development. The U.S. Geological Survey three-dimensional ground-water flow model MODFLOW-2000 was used to simulate ground-water flow in the surficial and Floridan aquifer systems in Lake County, the Ocala National Forest, and adjacent areas. A steady-state calibration to average 1998 conditions was facilitated by using the inverse modeling capabilities of MODFLOW-2000. Values of hydrologic properties from the calibrated model were in reasonably close agreement with independently estimated values and results from previous modeling studies. The calibrated model generally produced simulated water levels and flows in reasonably close agreement with measured values and was used to simulate the hydrologic effects of projected 2020 conditions. Ground-water withdrawals in the model area have been projected to increase from 470 million gallons per day in 1998 to 704 million gallons per day in 2020. Significant drawdowns were simulated in Lake County from average 1998 to projected 2020 conditions: the average and maximum drawdowns, respectively, were 0.5 and 5.7 feet in the surficial aquifer system, 1.1 and 7.6 feet in the Upper Floridan aquifer, and 1.4 and 4.3 feet in the Lower Floridan aquifer. The largest drawdowns in Lake County were simulated in the southeastern corner of the County and in the vicinities of Clermont and Mount Dora. Closed-basin lakes and wetlands are more likely to be affected by future pumping in these large drawdown areas, as opposed to other areas of Lake County. However, within the Ocala National Forest, drawdowns were relatively small: the average and maximum drawdowns, respectively, were 0.1 and 1.0 feet in the surficial aquifer system, 0.2 and

Water-Resources Investigations Report↗

Development of a digital model of ground-water flow in deeply weathered crystalline rock, Chester County, Pennsylvania

The model developed in this study simulates .recharge to, flow through, and discharge from the water-table aquifer in the upper Pickering Creek basin, a 5.98-square-mile basin representative of most of Chester County, Pennsylvania. The two-dimensional finite-difference model of Trescott, Pinder, and Larson was used with slight modification. The way ground-water evapotranspiration varies with depth was modified, and a minimum transmissivity was used to prevent the model from "going dry" during water-table simulations. Tests of the model showed good correlation with field data, but also showed that certain additional types of information would be useful to make the model more reliable. Particularly helpful would be more water-level and evapotranspiration data, synoptic base-flow measurements in the subbasins, and better information on the variation of permeability and specific yield with depth. Predictive simulations made by the model show the effects of imposing various changes on the hydrologic system. Residential development was simulated in 1 square mile of the basin by combinations of the following conditions: 1-acre, 1/2-acre, and 1/4-acre lot size; wastewater disposal by septic systems or by public sewers; domestic wells or public-supply wells for water supply; and 0, 1, and 5 percent reduction in recharge in the developed area. Where water is supplied locally by wells, decreased lot size and reduced recharge result in increased drawdowns and reduced base flow. When wastewater is exported by sewers, these effects are magnified. The simulations by the model show the relative severity of the effects of various development schemes.

Pennslyvania↗

Documentation of a digital spatial data base for hydrologic investigations, Broward County, Florida

Geographic information systems have become an important tool in planning for the protection and development of natural resources, including ground water and surface water. A digital spatial data base consisting of 18 data layers that can be accessed by a geographic information system was developed for Broward County, Florida. Five computer programs, including one that can be used to create documentation files for each data layer and four that can be used to create data layers from data files not already in geographic information system format, were also developed. Four types of data layers have been developed. Data layers for manmade features include major roads, municipal boundaries, the public land-survey section grid, land use, and underground storage tank facilities. The data layer for topographic features consists of surveyed point land-surface elevations. Data layers for hydrologic features include surface-water and rainfall data-collection stations, surface-water bodies, water-control district boundaries, and water-management basins. Data layers for hydrogeologic features include soil associations, transmissivity polygons, hydrogeologic unit depths, and a finite-difference model grid for south-central Broward County. Each data layer is documented as to the extent of the features, number of features, scale, data sources, and a description of the attribute tables where applicable.

Florida↗

Hydrology and geochemistry of a slag-affected aquifer and chemical characteristics of slag-affected ground water, northwestern Indiana and northeastern Illinois

Slag is a by-product of steel manufacturing and a ubiquitous fill material in northwestern Indiana. Ground water associated with slag deposits generally is characterized by high pH and elevated concentrations of many inorganic water-quality constituents. The U.S. Geological Survey, in cooperation with the Indiana Department of Environmental Management, conducted a study in northwestern Indiana from June 1995 to September 1996 to improve understanding of the effects of slag deposits on the water quality of a glacial-outwash aquifer. The Bairstow Landfill, a slag-fill deposit overlying the Calumet aquifer near Hammond, Indiana, was studied to represent conditions in slag-deposit settings that are common in northwestern Indiana. Ground water from 10 observation wells, located in four nests at the site, and surface water from the adjacent Lake George were analyzed for values of field-measured parameters and concentrations of major ions, nutrients, trace elements, and bulk properties. Solid-phase samples of slag and aquifer sediment collected during drilling were examined with X-ray diffraction and geochemical digestion and analysis. Concentrations of calcium, potassium, sodium, and sulfate were highest in wells screened partly or fully in slag. Potassium concentrations in ground water ranged from 2.9 to 120 milligrams per liter (mg/L), were highest in water from slag deposits, and decreased with depth. The highest concentrations for aluminum, barium, molybdenum, nickel, and selenium were in water from the slag. Silica concentrations were highest in wells screened directly beneath the slag-aquifer interface, and magnesium concentrations were highest in intermediate and deep aquifer wells. Silica concentrations in shallow and intermediate aquifer wells ranged from 27 to 41 mg/L and were about 10 times greater than those in water from slag deposits. The highest concentrations for chromium, lead, and zinc were in ground water from immediately below the slag-aquifer interface. The solid-phase analyses indicated that calcite, dolomite, and quartz generally were present throughout the slag-aquifer system; barian celestite, cristobalite, manganese-bearing calcite, and minrecordite were present in fewer samples. Trace elements that are liberated from the slag may be incorporated as impurities during precipitation of major minerals, sorbed onto clays and other grainsize fractions not analyzed as part of this study, or present in low-abundance minerals that were not detected by the X-ray analysis. Mass-balance and speciation programs were used to identify geochemical processes that may be occurring as water infiltrates through the slag, flows into the aquifer, and discharges into Lake George. The geochemical models indicate that precipitation of calcite may be occurring where slag-affected water enters the aquifer. Models also indicate that dolomite precipitation and clay-mineral dissolution may be occurring at the slag-aquifer interface; however, dolomite precipitation is generally believed to require geologically long time periods. Silica may be dissolving where slag-affected ground water enters the aquifer and may be precipitating where slag-affected ground water discharges to the lakebed of Lake George. In addition to the site-specific study, a statistical analysis of regional water quality was done to compare ground water in wells affected and unaffected by slag. When com-pared to wells in background locations in the Calumet aquifer, wells screened in slag across northwestern Indiana and northeastern Illinois generally had relatively higher pH and specific-conductance values and relatively higher concentrations of alkalinity, dissolved solids, suspended solids, total organic carbon, calcium, potassium, sodium, chloride, aluminum, barium, and possibly magnesium, sulfate, chromium, cobalt, copper, cyanide, manganese, mercury, nickel, and vanadium. When compared to wells in slag and wells in background locations, ground water from immediately beneath or immediately downgradient from slag had relatively high concentrations of arsenic and silica. Water-quality characteristics in ground water at the Bairstow Landfill were similar to water-quality characteristics in slag-contact and slag-affected wells throughout northwestern Indiana.

Illinois, Indiana↗

Hydrogeology and simulation of source areas of water to production wells in a colluvium-mantled carbonate-bedrock aquifer near Shippensburg, Cumberland and Franklin Counties, Pennsylvania

This report presents the results of a study by the U.S. Geological Survey in cooperation with the Shippensburg Borough Authority to evaluate the source areas of water to production wells in a colluvium-mantled carbonate-bedrock aquifer in Cumberland and Franklin Counties, Pa. The areal extent of the zone of contribution was simulated for three production wells near Shippensburg, Pa. by use of a ground-water-flow model. A 111-square-mile area was selected as the model area and includes areas of the South Mountain Section and the Great Valley Section of the Valley and Ridge Physiographic Province. Within the model area, the geologic units in the South Mountain area are predominantly metamorphic rocks and the geologic units in the Great Valley are predominantly carbonate rocks. Hydrologic and geologic information were compiled to establish a conceptual model of ground-water flow. Characteristics of aquifer materials were determined, and streamflow and water levels were measured. Streamflow measurements in November 2003 showed all streams lost water as they flowed from South Mountain over the colluvium-mantled carbonate aquifer into the Great Valley. Some streams lost more than 1 cubic foot per second to the aquifer in this area. The Shippensburg Borough Authority owns three production wells in the model area. Two wells, Cu 969 and Fr 823, are currently (2004) used as production wells and produce 500,000 and 800,000 gallons per day, respectively. Well Cu 970 is intended to be brought on line as a production well in the future. Water levels were measured in 43 wells to use for model calibration. Water-level fluctuations and geophysical logs indicated confined conditions in well Cu 970. Ground-water flow was simulated with a model that consisted of two vertical layers, with five zones in each layer. The units were hydrostratigraphic units that initially were based on geologic formations, but boundaries were adjusted during model calibration. Model calibration resulted in a root mean square error of 9.8 feet. A parameter-estimation package was used during model calibration to estimate three parameters. The parameter estimation resulted in a value of 233 feet per day for horizontal hydraulic conductivity of the highly fractured carbonate rocks and sandy colluvium in layer 1; 3.97 feet per day for horizontal hydraulic conductivity of the ridge-forming unit in layer 1; and a value of 1.73 for horizontal anisotropy in both layers. The calibrated model was used to delineate the areal extent of the zone of contribution for wells Cu 969 and Fr 823. Although well Cu 970 is not currently (2004) being used, the areal extent of its zone of contribution also was simulated without additional model calibration. The shape of the areal extent of the zone of contribution was similar for each well and included an area that extended from the well southwest along the Tomstown Formation, and then extended southeast into the metamorphic rocks of South Mountain. The contributing areas from the watersheds of losing streams were also delineated because losing stream reaches bisect the areal extent of the zones of contribution. Spatial uncertainty of the areal extent of the zone of contribution was illustrated using a Monte-Carlo analysis. The model was run 1,000 times using randomly generated parameter sets that were normally distributed within the confidence interval around the optimal values for the three estimated parameters. The model converged and had a reasonable water budget for 980 of the model runs. For each of those 980 model runs, the recharge area was determined, and the results for all runs were compiled and contoured. The results of the Monte-Carlo analysis were compared to the results of the deterministic model, illustrating that the deterministic model has the greatest certainty in the area closest to each well in the Tomstown Formation. The areas farther from the well, upgradient, and in the metamorphic rocks have a higher degree

Pennsylvania↗

Simulated interactions between the proposed Narrows Reservoir and the water-table aquifer along the South Platte River, Morgan County, Colorado

A computer model, including a ground-water-flow component and a mass balance reservoir-operations component, was developed to simulate the proposed Narrows Reservoir and the adjacent alluvial aquifer of the South Platte River, Morgan County, Colo. This model, using a weekly time step, simulated the transient interactions of these two systems for an initial-fill condition and general-operational condition. A sensitivity analysis was made to test the effects of possible errors in the description of aquifer characteristics on the model results. The initial-fill simulation indicated that to fill the reservoir when hydraulic connection between the surface-water system and the aquifer is simulated would take 2 additional years than if no connection was assumed. Simulated ground-water return flow to the river downstream of the proposed dam was about 85 percent of the estimated maximum values, computed under steady-state conditions in an earlier study. The general operation simulation indicated that during the period of lowest reservoir contents the aquifer provided about 80,000 acre-feet of recoverable storage to the reservoir 's capacity. Average return flow was only about 70 percent of the estimated maximum values computed in the earlier, steady-state analysis. Monthly ground-water outflow from the reservoir to the aquifer ranged from 55,200 to -10,400 acre-feet. Parameters tested for sensitivity, in order of decreasing sensitivity, were: Hydraulic conductivity, specific yield, hydraulic connection between reservoir and aquifer, local recharge, boundary conditions, and dam permeability. The probable error in the parameters tested would not seem to warrant significant new data collection. (USGS)

Colorado↗

Simulation of ground-water flow and areas contributing recharge to extraction wells at the Drake Chemical Superfund Site, City of Lock Haven and Castanea Township, Clinton County, Pennsylvania

Extensive remediation of the Drake Chemical Superfund Site has been ongoing since 1983. Contaminated soils were excavated and incinerated on site between 1996 and 1999. After 1999, remedial efforts focused on contaminated ground water. A ground-water remediation system was started in November 2000. The source area of the contaminated ground water was assumed to be the zone 1 area on the Drake Chemical site. The remedial system was designed to capture ground water migrating from zone 1. Also, the remediation system was designed to pump and treat the water in an anoxic environment and re-infiltrate the treated water underground through an infiltration gallery that is hydrologically downgradient of the extraction wells. A numerical ground-water flow model of the surrounding region was constructed to simulate the areas contributing recharge to remedial extraction wells installed on the Drake Chemical site. The three-dimensional numerical flow model was calibrated using the parameter-estimation process in MODFLOW-2000. The model included three layers that represented three poorly sorted alluvial sediment units that were characterized from geologic well and boring logs. Steady-state ground-water flow was simulated to estimate the areas contributing recharge to three extraction wells for three different pumping scenarios--all wells pumping at 2 gallons per minute, at approximately 5 gallons per minute, and at 8 gallons per minute. Simulation results showed the contributing areas to the three extraction wells encompassed 92 percent of zone 1 at a pumping rate of approximately 5 gallons per minute. The contributing areas did not include a very small area in the southwestern part of zone 1 when the three extraction wells were pumped at approximately 5 gallons per minute. Pumping from a fourth extraction well in that area was discontinued early in the operation of the remediation system because the ground water in that area met performance standards. The areas contributing recharge to the three extraction wells did encompass zone 1 at a pumping rate of 8 gallons per minute. At pumping rates of 2 gallons per minute, the contributing areas for the three extraction wells did not encompass zone 1.

Pennsylvania↗

Initial Results from a Study of Climatic Changes and the Effect on Wild Sheep Habitat in Selected Study Areas of Alaska

Climate change theorists have projected striking changes in local weather on earth due to increases in temperature. These predicted changes may cause melting glaciers and ice caps, rising sea levels, increasing desertification and other environmental changes which seem likely to affect presumed indicator species as harbingers of more significant changes. Wild sheep, even though they are one of the more successful mammalian taxa since Pleistocene times, exhibit a suite of adaptations to glacier driven environments which may be presumed to render them sensitive to environmental changes. The authors began investigation with these assumptions by comparing changes, as determined by satellite imagery, in glacier extent in our study areas in Denali National Park, Alaska, during the last 30 years. Our findings showed the extent of glacial retreat in Alaska during this time period was approximately 40-50 percent as measured by ablation zone and retreat of terminal moraines. During the first half of this 30-year period, Dall sheep (Ovis dalli dalli) populations were stable at historically recorded highs. In the early to mid-1990s, Dall sheep populations in Alaska declined from an historical estimated high of 75,000 sheep to the presently estimated 40-50,000. The declines seemed to be weather related, on the basis of the presumption that lamb survival rates are primarily weather-mediated in Alaska. Changes in local weather appear, at this point, to be correlated with oscillation in the Pacific Current in the Northern Pacific ocean. Of course, changes in local weather affect forage abundance and quality seasonally. In investigating a possible linkage of weather to seasonal forage abundance and quality, we also investigated changes in snow and ice extent and distribution, as well as increased water runoff associated with permafrost and depleted glaciers. Databases were assembled from a wide variety of remotely sensed satellite data, ground-based observations, and historical data bases relating to Dall sheep habitats in selected study areas. Alaska's sheep habitats are typified by long, narrow bands of mountainous uplifts generally arrayed west-to-east, and perpendicular to prevailing south-to-north weather-front movements. Classic Dall sheep habitat occurs on snow-shadowed slopes within these narrow mountainous habitats. On the basis of these data, we offer an explanatory hypothesis relating Dall sheep welfare to weather and climate-influenced nutrition and a monitoring scheme, which should produce data sufficient to test the robustness of this hypothesis. If correlated with population changes, the methods used in our comparative observations may provide long-term monitoring tools for wildlife managers and be applicable in other widely-dispersed wild sheep habitats. If no significant correlations emerge from our modeling exercises, the notion that wild sheep are a sufficiently sensitive species to be seen as an indicator species will have to be reexamined.

Open-File Report↗

Hydrology of Laguna Joyuda, Puerto Rico

A study was conducted by the U.S. Geological Survey to define the hydraulic and hydrologic characteristics of the Laguna Joyuda system (in southwestern Puerto Rico) and to determine the water budget of the lagoon. This shallow-water lagoon is connected to the sea by a single canal. Rainfall and evaporation, surface-water, groundwater, and tidal-flow data were collected from December 1, 1985, to April 30, 1988. A conceptual hydrologic model of the lagoon was developed and discharge measurements and modeling were undertaken to quantify the different flow components. The water balance during the 29-month study period was determined by measuring and estimating the different hydrologic components: 4.14 million cubic meters rainfall; 5.38 million cubic meters evaporation; 1.1 8 million cubic meters surface water; and 0.34 million cubic meters ground water. A total of 18.9 million cubic meters ebb flow (tidal outflow) was discharged from the lagoon and 14.4 million cubic meters flood flow (tidal inflow) entered through the canal during the study. Seawater inflow accounted for 71 percent of the water into the lagoon. The storage volume of the lagoon was about 1.55 million cubic meters. The lagoon's hydrologic-budget residual was 4.22 million cubic meters, whereas the sum of the estimated errors for the different hydrologic components amounted to 4.51 million cubic meters. Average flushing rate for the lagoon was estimated at 72 days. During the study, the specific conductance of the lagoon water ranged from 32,000 to 52,000 microsiemens per centimeter at 25 degrees Celsius, whereas the specific conductance of local seawater is about 45,000 to 55,000 microsiemens.

Water-Resources Investigations Report↗

Site Selection for a Deep Monitor Well, Kualapuu, Molokai, Hawaii

Management of the ground-water resources near Kualapuu on the island of Molokai, Hawaii, is hindered by the uncertainty in the vertical salinity structure in the aquifer. In the State of Hawaii, vertical profiles of ground-water salinity are commonly obtained from deep monitor wells, and these profiles are used to estimate the thicknesses of the freshwater part of the ground-water flow system and the freshwater-saltwater transition zone. Information from a deep monitor well would improve the understanding of the ground-water flow system and the ability to effectively manage the ground-water resources near Kualapuu; however, as of mid-1999 no deep monitor wells had been drilled on the island of Molokai. Selection of an appropriate site for drilling a deep monitor well in the Kualapuu area depends partly on where future ground-water development may occur. Simulations using an areally two-dimensional, steady-state, sharp-interface ground-water flow model previously developed for the island of Molokai, Hawaii, indicate that the southeastern part of the Kualapuu area is a possible area of future ground-water development because (1) withdrawals from this area have a small effect on water levels at existing wells in the Kualapuu area (relative to effects from withdrawals in other parts of the Kualapuu area that are outside of the dike complex), and (2) model-calculated water levels in this part of the Kualapuu area are high relative to water levels in other parts of the Kualapuu area that are outside of the dike complex. Additional site-selection criteria include (1) ground-water level, (2) ground-surface altitude, (3) land classification, ownership, and accessibility, (4) geology, (5) locations of existing production wells, and (6) historical ground-water quality information. A deep monitor well in the Kualapuu area will likely be most useful for management purposes if it is located (1) in the vicinity of future ground-water development, (2) in an area where water levels are between about 8 and 12 feet above sea level, (3) at a ground-surface altitude that is between about 1,000 and 1,100 feet, (4) on government-owned land, (5) outside of the dike complex and as far from known volcanic vents as possible, (6) at least about 1,000 feet from, but within the same hydrogeologic setting as, existing or proposed production wells, and (7) east of well 0902-01. A viable area for drilling a deep monitor well is about a half mile southeast of existing wells 0801-01 to -03 and a half mile north of a known volcanic vent, Puu Luahine.

Water-Resources Investigations Report↗

Revision of the documentation for a model for calculating effects of liquid waste disposal in deep saline aquifers

The model developed under this contract is a modified version of the deep well disposal model developed by INTERCOMP Resource Development and Engineering, Inc., for the U.S. Geological Survey (A model for calculating effects of liquid waste disposal in deep saline aquifers). The model is a finite-difference numerical solution of the partial differential equations describing (1) single phase fluid flow in aquifer, (2) energy transport by convection and conduction, and (3) contaminant transport dissolved in the fluid by convection and dispersion. Both the energy and the contaminant transports include molecular diffusion and hydrodynamic dispersion. The objective of the contract was to modify the original version of the model for more general applications. Some of the major features are as follows: (1) Fluid density is permitted to be a function of fluid pressure, temperature and contaminant concentration, the viscosity can be described as a function of pressure and temperature. (2) Aquifer heterogeneities in the hydrological properties can be described on a numerical grid block basis. (3) Free water surface is permitted to exist in the aquifer for shallow ground water applications. (4) Contaminant may absorb on rock surface or decay according to a first order reaction. The absorption is described by a linear adsorption isotherm function of rock type in the aquifer. (5) The model is extremely flexible in providing a wide choice of boundary conditions. These include natural flow in the aquifer, aquifer influence functions around the perimeter of the grid in recognition that the gridded region does not have no-flow boundaries, heat losses into the overlying impermeable strata, and the wellbore heat and pressure drop calculations coupled to the aquifer flow equations. (6) The model offers the option of selecting as iterative or direct solution technique, and selecting central or backward finite-difference approximations in both time and space. (7) Virtually any aquifer can be modeled by proper grid block description in three dimensions. In addition, the model is fully transient. The major limitation of the model arises using second-order correct (central-difference) finite-difference approximation in space. To avoid numerical oscillations in the solution, the user must restrict grid block and time step sizes depending upon the magnitude of the dispersivity.

Water-Resources Investigations Report↗

Hydrogeologic facies characterization of an alluvial fan near Fresno, California, using geophysical techniques

DBCP (1,2-dibromo-3-chloropropane) contamination in the sole source aquifer near Fresno, California, has significantly affected drinking-water supplies. Borehole and surface geophysical data were integrated with borehole textural data to characterize the Kings River alluvial fan sediments and to provide a framework for computer modeling of pesticide transport in ground water. Primary hydrogeologic facies units, such as gravel, coarse sand or gravel, fine sand, and silt and clay, were identified in cores collected from three borings located on a 4.6-kilometer transect of multilevel monitoring wells. Borehole geophysical logs collected from seven wells and surface geophysical surveys were used to extrapolate hydrogeologic facies to depths of about 82meters and to correlate the facies units with neighboring drilling sites. Thickness ranged from 0.3to 13 meters for sand and gravel units, and from 0.3 to 17 meters for silt and clay. The lateral extent of distinct silt and clay layers was mapped using shallow seismic reflection and ground-penetrating radar techniques. About 3.6 kilometers of seismic reflection data were collected; at least three distinct fine-grained layers were mapped. The depth of investigation of the seismic survey ranged from 34 to 107 meters below land surface, and vertical resolution was about 3.5 meters. The ground-penetrating radar survey covered 3.6kilometers and imaged a 1.5-meters thick, continuous fine-grained layer located at a depth of about 8 meters. Integrated results from the borehole sediment descriptions and geophysical surveys provided a detailed characterization over a larger areal extent than traditional hydrogeologic methods alone.

Open-File Report↗

Hydrogeologic framework refinement, ground-water flow and storage, water-chemistry analyses, and water-budget components of the Yuma area, southwestern Arizona and southeastern California

The ground-water and surface-water system in the Yuma area in southwestern Arizona and southeastern California is managed intensely to meet water-delivery requirements of customers in the United States, to manage high ground-water levels in the valleys, and to maintain treaty-mandated water-quality and quantity requirements of Mexico. The following components in this report, which were identified to be useful in the development of a ground-water management model, are: (1) refinement of the hydrogeologic framework; (2) updated water-level maps, general ground-water flow patterns, and an estimate of the amount of ground water stored in the mound under Yuma Mesa; (3) review and documentation of the ground-water budget calculated by the Bureau of Reclamation, U.S. Department of the Interior (Reclamation); and (4) water-chemistry characterization to identify the spatial distribution of water quality, information on sources and ages of ground water, and information about the productive-interval depths of the aquifer. A refined three-dimensional digital hydrogeologic framework model includes the following hydrogeologic units from bottom to top: (1) the effective hydrologic basement of the basin aquifer, which includes the Pliocene Bouse Formation, Tertiary volcanic and sedimentary rocks, and pre-Tertiary metamorphic and plutonic rocks; (2) undifferentiated lower units to represent the Pliocene transition zone and wedge zone; (3) coarse-gravel unit; (4) lower, middle, and upper basin fill to represent the upper, fine-grained zone between the top of the coarse-gravel unit and the land surface; and (5) clay A and clay B. Data for the refined model includes digital elevation models, borehole lithology data, geophysical data, and structural data to represent the geometry of the hydrogeologic units. The top surface of the coarse-gravel unit, defined by using borehole and geophysical data, varies similarly to terraces resulting from the down cutting of the Colorado River. Clay A is nearly the same as the previous conceptual hydrogeologic model definition (Olmsted and others, 1973), except for a minor westward extension from the city of Yuma. Clay B is extended to the southerly international boundary and increased in areal extent by about two-thirds of the original extent (Olmsted and others, 1973). The other hydrogeologic units generally are the same as in the previous conceptual hydrogeologic model. Before development, the Colorado and Gila Rivers were the sources of nearly all the ground water in the Yuma area through direct infiltration of water from river channels and annual overbank flooding. After construction of upstream reservoirs and clearing and irrigation of the floodplains, the rivers now act as drains for the ground water. Ground-water levels in most of the Yuma area are higher now than they were in predevelopment time. A general gradient of ground-water flow toward the natural discharge area south of the Yuma area still exists, but many other changes in flow are evident. Ground water in Yuma Valley once flowed away from the Colorado River, but now has a component of flow towards the river and Mexicali Valley. A ground-water mound has formed under Yuma Mesa from long-term surface-water irrigation; about 600,000 to 800,000 acre-ft of water are stored in the mound. Ground-water withdrawals adjacent to the southerly international boundary have resulted in water-level declines in that area. The reviewed and documented water budget includes the following components: (1) recharge in irrigated areas, (2) evapotranspiration by irrigated crops and phreatophytes, (3) ground-water return flow to the Colorado River, and (4) ground-water withdrawals (including those in Mexicali Valley). Recharge components were calculated by subtracting the amount of water used by crops from the amount of water delivered. Evapotranspiration rates were calculated on the basis of established methods, thus were appropriate for input to the ground-water flow model developed by the Bureau of Reclamation (William Greer, hydrologist, Bureau of Reclamation, written commun., 2005). Evapotranspiration by crops and phreatophytes were calculated by using crop coefficient methods and meteorological data. Other methods of calculating evapotranspiration rates by using combinations of satellite imagery and ground-based data could be used for higher spatial and temporal resolution. Ground-water return flow during years of low flow on the Colorado River (1972–82, 1987–92, and 1994–96) averaged 79,000 acre-ft per year. Ground-water withdrawal data for 1970–99 were similar to other estimates made by the U.S. Geological Survey for the Yuma area. New water-chemistry data were collected in 12 wells and 8 canals/drains to characterize spatial patterns in chemical constituents, determine isotopic ages of water, infer possible sources of ground water, and locate the vertical intervals of the aquifer that contribute most water to wells. Depth-dependent samples were collected at one of the wells (YM-10). A large quantity of water-quality data were compiled from Bureau of Reclamation and U.S. Geological Survey records and merged into the U.S. Geological Survey National Water Information System database. New samples were analyzed for major ions, nutrients, stable isotopes of oxygen and hydrogen, tritium ( 3 H), and carbon-14 ( 14 C) (along with C 13 /C 12 ratios). Light values of oxygen-18 ( 18 O) and deuterium ( 2 H, D) in well 242-2 indicate recharge from the Colorado River. Heavy water samples from wells 242-22, CADC, and Mesa del Sol indicate local recharge sources. Tritium data indicate there is young water in wells in the valleys and near the edge of Yuma Mesa, while older water is found far from the Colorado River. 14 C data indicate that water from wells near the southerly international boundary is at least several thousand years old.

Arizona, California↗

Preliminary evaluation of ground-water flow in Bear Creek Valley, the Oak Ridge Reservation, Tennessee

Bear Creek Valley, Tennessee contains hazardous waste disposal sites where contaminants leach into ground and surface water. Groundwater flow and the potential migration of contaminants is poorly understood. The Valley is underlain by calcareous shale that contains limestone units. Ridges to the north and south are underlain by interbedded sandstones, siltstone and shale, and by massive, siliceous dolomite, respectively. The bedrock, which dips about 45 degrees southeast, is overlain by regolith to a maximum thickness of 80 ft. Observed hydraulic conductivities for the regolith range from 0.01 to 13 ft/day, and for the bedrock, from 0.001 to 11 ft/day. Groundwater flow is probably toward streams and is preferential along strike because of an areal anisotropy in hydraulic conductivity. A cross sectional groundwater flow model was used to test the conceptualized flow system and to help identify areas where additional data are needed. The preliminary model shows a pattern of recharge at both ridges, flow toward the valley, and upward flow that discharges into Bear Creek. Final model values of hydraulic conductivity in the bedrock range from 0.01 to 0.1 ft/day and reflect an areal anisotropy ratio of 1:5. Simulated recharge was 10 inches/year. (USGS)

Tennessee↗

Preliminary digital model of the Arikaree aquifer in the Sweetwater River basin, central Wyoming

In central Wyoming, Potentially large supplies of ground water are available in the Sweetwater River basin from the Arikaree aquifer, which consists of the upper part of the White River, the Arikaree, and the Ogallala Formations. A preliminary digital model was developed for the Arikaree aquifer using a small amount of poorly distributed data, an estimated distribution of recharge, and a conceptual model of the Arikaree aquifer flow system. Calibration of the model was based on reproduction of the potentiometric surface and the base flow of the Sweetwater River in November 1975. Calculated steady-state hydraulic heads were within 50 feet of the observed heads in about 98 percent of the nodes. The calculated leakage from the Arikaree aquifer to he Sweetwater River in the western area was within about 12 percent of the leadage determined by gain and loss studies. In order to develop a comprehensive digital model that would respond to hydraulic stress in nearly the same manner as the actual aquifer flow system, measured responses of the aquifer to stress are needed. Also needed are additional data on aquifer characteristics, recharge to the aquifer, and stream-aquifer relationships. (Woodard-USGS)

Water-Resources Investigations Report↗

Description and comparison of selected models for hydrologic analysis of ground-water flow, St. Joseph River basin, Indiana

The Indiana Department of Natural Resources (IDNR) is developing water-management policies designed to assess the effects of irrigation and other water uses on water supply in the basin. In support of this effort, the USGS, in cooperation with IDNR, began a study to evaluate appropriate methods for analyzing the effects of pumping on ground-water levels and streamflow in the basin 's glacial aquifer systems. Four analytical models describe drawdown for a nonleaky, confined aquifer and fully penetrating well; a leaky, confined aquifer and fully penetrating well; a leaky, confined aquifer and partially penetrating well; and an unconfined aquifer and partially penetrating well. Analytical equations, simplifying assumptions, and methods of application are described for each model. In addition to these four models, several other analytical models were used to predict the effects of ground-water pumping on water levels in the aquifer and on streamflow in local areas with up to two pumping wells. Analytical models for a variety of other hydrogeologic conditions are cited. A digital ground-water flow model was used to describe how a numerical model can be applied to a glacial aquifer system. The numerical model was used to predict the effects of six pumping plans in 46.5 sq mi area with as many as 150 wells. Water budgets for the six pumping plans were used to estimate the effect of pumping on streamflow reduction. Results of the analytical and numerical models indicate that, in general, the glacial aquifers in the basin are highly permeable. Radial hydraulic conductivity calculated by the analytical models ranged from 280 to 600 ft/day, compared to 210 and 360 ft/day used in the numerical model. Maximum seasonal pumping for irrigation produced maximum calculated drawdown of only one-fourth of available drawdown and reduced streamflow by as much as 21%. Analytical models are useful in estimating aquifer properties and predicting local effects of pumping in areas with simple lithology and boundary conditions and with few pumping wells. Numerical models are useful in regional areas with complex hydrogeology with many pumping wells and provide detailed water budgets useful for estimating the sources of water in pumping simulations. Numerical models are useful in constructing flow nets. The choice of which type of model to use is also based on the nature and scope of questions to be answered and on the degree of accuracy required.

Indiana, Michigan↗