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 649 records · Page 36Linked to original sources

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 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

Hydrogeology and simulation of flow between the alluvial and bedrock aquifers in the upper Black Squirrel Creek basin, El Paso County, Colorado

Anticipated increases in pumping from the bedrock aquifers in El Paso County potentially could affect the direction and rate of flow between the alluvial and bedrock aquifers and lower water levels in the overlying alluvial aquifer. The alluvial aquifer underlies about 90 square miles in the upper Black Squirrel Creek Basin of eastern El Paso County. The alluvial aquifer consists of unconsolidated alluvial deposits that unconformably overlie siltstones, sandstones, and conglomerate (bedrock aquifers) and claystone, shale, and coal (bedrock confining units) of the Denver Basin. The bedrock aquifers (Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers) are separated by confining units (upper and lower Denver and the Laramie confining units) and overlie a relatively thick and impermeable Pierre confining unit. The Pierre confining unit is assumed to be a no-flow boundary at the base of the alluvial/ bedrock aquifer system. During 1949-90, substantial water-level declines, as large as 50 feet, in the alluvial aquifer resulted from withdrawals from the alluvial aquifer for irrigation and municipal supplies. Average recharge to the alluvial aquifer from infiltration of precipitation and surface water was an estimated 11.97 cubic feet per second and from the underlying bedrock aquifers was an estimated 0.87 cubic foot per second. Water-level data from eight bedrock observation wells and eight nearby alluvial wells indicate that, locally, the alluvial and bedrock aquifers probably are hydraulically connected and that the alluvial aquifer in the upper Black Squirrel Creek Basin receives recharge from the Denver and Arapahoe aquifers but-locally recharges the Laramie-Fox Hills aquifer. Subsurface-temperature profiles were evaluated as a means of estimating specific discharge across the bedrock surface (the base of the alluvial aquifer). However, assumptions of the analytical method were not met by field conditions and, thus, analyses of subsurface-temperature profiles did not reliably estimate specific discharge across the bedrock surface. The vertical hydraulic diffusivity of a siltstone and sandstone in the lower Denver confining unit was estimated, by an aquifer test, to be about 8 x 10'4 square foot per day. Physical and chemical characteristics of water from the bedrock aquifers in the study area generally differ from the physical and chemical characteristics of water from the alluvial aquifer, except for the physical and chemical characteristics of water from one bedrock well, which is completed in the Laramie-Fox Hills aquifer. In the southern part of the study area, physical and chemical characteristics of ground water indicate downward flow of water from the alluvial aquifer to the Laramie-Fox Hills aquifer. A three-dimensional numerical model was used to evaluate flow of water between the alluvial aquifer and underlying bedrock. Simulation of steady-state conditions indicates that flow from the bedrock aquifers to the alluvial aquifer was about 7 percent of recharge to the alluvial aquifer, about 0.87 cubic foot per second. The potential effects of withdrawal from the alluvial and bedrock aquifers at estimated (October 1989 to September 1990) rates and from the bedrock aquifers at two larger hypothetical rates were simulated for a 50-year projection period. The model simulations indicate that water levels in the alluvial aquifer will decline an average of 8.6 feet after 50 years of pumping at estimated October 1989 to September 1990 rates. Increases in withdrawals from the bedrock aquifers in El Paso County were simulated to: (1) Capture flow that currently discharges from the bedrock aquifers to springs and streams in upland areas and to the alluvial aquifer, (2) induce flow downward from the alluvial aquifer, and (3) accelerate the rate of waterlevel decline in the alluvial aquifer.

Colorado

Simulation of net infiltration and potential recharge using a distributed-parameter watershed model of the Death Valley region, Nevada and California

This report presents the development and application of the distributed-parameter watershed model, INFILv3, for estimating the temporal and spatial distribution of net infiltration and potential recharge in the Death Valley region, Nevada and California. The estimates of net infiltration quantify the downward drainage of water across the lower boundary of the root zone and are used to indicate potential recharge under variable climate conditions and drainage basin characteristics. Spatial variability in recharge in the Death Valley region likely is high owing to large differences in precipitation, potential evapotranspiration, bedrock permeability, soil thickness, vegetation characteristics, and contributions to recharge along active stream channels. The quantity and spatial distribution of recharge representing the effects of variable climatic conditions and drainage basin characteristics on recharge are needed to reduce uncertainty in modeling ground-water flow. The U.S. Geological Survey, in cooperation with the Department of Energy, developed a regional saturated-zone ground-water flow model of the Death Valley regional ground-water flow system to help evaluate the current hydrogeologic system and the potential effects of natural or human-induced changes. Although previous estimates of recharge have been made for most areas of the Death Valley region, including the area defined by the boundary of the Death Valley regional ground-water flow system, the uncertainty of these estimates is high, and the spatial and temporal variability of the recharge in these basins has not been quantified. To estimate the magnitude and distribution of potential recharge in response to variable climate and spatially varying drainage basin characteristics, the INFILv3 model uses a daily water-balance model of the root zone with a primarily deterministic representation of the processes controlling net infiltration and potential recharge. The daily water balance includes precipitation (as either rain or snow), snow accumulation, sublimation, snowmelt, infiltration into the root zone, evapotranspiration, drainage, water content change throughout the root-zone profile (represented as a 6-layered system), runoff (defined as excess rainfall and snowmelt) and surface water run-on (defined as runoff that is routed downstream), and net infiltration (simulated as drainage from the bottom root-zone layer). Potential evapotranspiration is simulated using an hourly solar radiation model to simulate daily net radiation, and daily evapotranspiration is simulated as an empirical function of root zone water content and potential evapotranspiration. The model uses daily climate records of precipitation and air temperature from a regionally distributed network of 132 climate stations and a spatially distributed representation of drainage basin characteristics defined by topography, geology, soils, and vegetation to simulate daily net infiltration at all locations, including stream channels with intermittent streamflow in response to runoff from rain and snowmelt. The temporal distribution of daily, monthly, and annual net infiltration can be used to evaluate the potential effect of future climatic conditions on potential recharge. The INFILv3 model inputs representing drainage basin characteristics were developed using a geographic information system (GIS) to define a set of spatially distributed input parameters uniquely assigned to each grid cell of the INFILv3 model grid. The model grid, which was defined by a digital elevation model (DEM) of the Death Valley region, consists of 1,252,418 model grid cells with a uniform grid cell dimension of 278.5 meters in the north-south and east-west directions. The elevation values from the DEM were used with monthly regression models developed from the daily climate data to estimate the spatial distribution of daily precipitation and air temperature. The elevation values were also used to simulate atmosp

Water-Resources Investigations Report

Digital model for simulated effects of ground-water pumping in the Hueco Bolson, El Paso Area, Texas, New Mexico, and Mexico

The Hueco Bolson provides a substantial part of the municipal and industrial water supply of the El Paso area of Texas, New Mexico, and Mexico. Although the supply 6f fresh ground water in the bolson is large, about 10.6 million acre-feet (13,070 hm 3 ) in 1973 in the Texas part of the bolson alone, the supply is being depleted. A two-layer digital model of the Hueco Bolson was developed to duplicate the historic changes in water levels and the predevelopment (1903) water surface, and'to predict the response of the hydrologic system to any plan of water development. The model study showed that a proposed plan of development for 1973-91, in which pumping would be increased by 29 percent in the Texas part of the bolson and by 34 percent in Ciudad Juarez, would cause additional waterlevel declines of as much as 45 feet (13.7 m) in the vicinity of El Paso and 70 feet (21.3 m) in Ciudad Juarez. The study also showed that nearly 60 percent of the water would come from storage in the water-table part of the bolson aquifer (model layer 2) and 28 percent from leakage from the alluvium (model layer 1). By the end of the period, 9.84 million acrefeet (12,133 hm 3 ) of fresh water would be in storage in the Texas part of the bolson, as compared to 10.6 million acre-feet (13,070 hm 3 ) in storage in 1973.

Texas

Potentiometric surface of the Edwards-Trinity aquifer system and contiguous hydraulically connected units, west-central Texas, winter, 1974-75

The potentiometric surface of the Edwards-Trinity aquifer system and contiguous hydraulically connected units (from December 1974 through February 1975) was mapped as part of the Edwards-Trinity Regional Aquifer-System Analysis (RASA) project. A major goal of the Edwards-Trinity RASA project is to understand and describe the regional flow system (Bush, 1986). The development of a digital ground-water flow model of the aquifer system is a key part of the project. This potentiometric map will be used in the calibration of the ground-water flow model and in understanding ground-water movement in the aquifer system. The map depicts the potentiometric surface of the major aquifers of the Edwards-Trinity aquifer system and contiguous units that form a continuous hydraulically connected regional aquifer within the study area in west-central Texas (fig. 1). The potentiometric surface of an aquifer is an imaginary surface defined by contouring locations of equal static head (the altitude to which water will rise in a well). The potentiometric surface map shows the direction of ground-water flow from higher to lower altitude. The study area extends beyond the aquifers of the Edwards-Trinity system to hydrologic divides, including the Colorado River and the Rio Grande (fig. 2). The data used to compile this map were obtained from the Texas Natural Resources Information System on magnetic tape and from Rees and Buckner (1980). The winter of 1974-75 (December 1974 through February 1975) was selected for mapping for two reasons: (1) More water-level data were available throughout the study area for this winter season than for other winter seasons, and (2) during winter there is almost no loss of ground water as a result of evaporation, irrigation withdrawals, and transpiration.

Texas

The Effectiveness of Cattlemans Detention Basin, South Lake Tahoe, California

Lake Tahoe (Nevada-California) has been designated as an 'outstanding national water resource' by the U.S. Environmental Protection Agency, in part, for its exceptional clarity. Water clarity in Lake Tahoe, however, has been declining at a rate of about one foot per year for more than 35 years. To decrease the amount of sediment and nutrients delivered to the lake by way of alpine streams, wetlands and stormwater detention basins have been installed at several locations around the lake. Although an improvement in stormwater and snowmelt runoff quality has been measured, the effectiveness of the detention basins for increasing the clarity of Lake Tahoe needs further study. It is possible that poor ground-water quality conditions exist beneath the detention basins and adjacent wetlands and that the presence of the basins has altered ground-water flow paths to nearby streams. A hydrogeochemical and ground-water flow modeling study was done at Cattlemans detention basin, situated adjacent to Cold Creek, a tributary to Lake Tahoe, to determine whether the focusing of storm and snowmelt runoff into a confined area has (1) modified the ground-water flow system beneath the detention basin and affected transport of sediment and nutrients to nearby streams and (2) provided an increased source of solutes which has changed the distribution of nutrients and affected nutrient transport rates beneath the basin. Results of slug tests and ground-water flow modeling suggest that ground water flows unrestricted northwest across the detention basin through the meadow. The modeling also indicates that seasonal flow patterns and flow direction remain similar from year to year under transient conditions. Model results imply that about 34 percent (0.004 ft3/s) of the total ground water within the model area originates from the detention basin. Of the 0.004 ft3/s, about 45 percent discharges to Cold Creek within the modeled area downstream of the detention basin. The remaining 55 percent of ground water is either consumed by evapotranspiration, is discharged to Cold Creek outside the modeled area downstream of the detention basin, or is discharged directly to Lake Tahoe. Of the 45 percent discharging to Cold Creek, about 9 percent enters directly downstream of the detention basin and 36 percent enters further downstream. Geochemical and microbial data suggest that a seasonal variation of chemical constituents and microbe population size is present at most wells. The geochemical data also indicate that construction of Cattlemans detention basin has not substantially changed the composition of the ground water in the area. High concentrations of ammonia, iron, and dissolved organic carbon, low concentrations of sulfate and nitrate, and large populations of sulfate-reducing microbes imply that the major geochemical process controlling nutrient concentrations beneath the detention basin is sulfate reduction. High concentrations of total nitrogen indicate that oxidation of organic carbon is a second important geochemical process occurring beneath the basin. The influx of surface runoff during spring 2002 apparently provided sufficient oxidized organic carbon to produce iron-reducing conditions and an increase in reduced iron, sulfate, and iron-reducing microorganisms. The increase in recharge of oxygenated water to the ground water system beneath the basin in future intervals of increased recharge may eventually redistribute nutrients and speed up transport of dissolved nutrients from the ground water system to Cold Creek.

Scientific Investigations Report

Estimating Water Fluxes Across the Sediment-Water Interface in the Lower Merced River, California

The lower Merced River Basin was chosen by the U.S. Geological Survey?s (USGS) National Water Quality Assessment Program (NAWQA) to be included in a national study on how hydrological processes and agricultural practices interact to affect the transport and fate of agricultural chemicals. As part of this effort, surface-water?ground-water (sw?gw) interactions were studied in an instrumented 100-m reach on the lower Merced River. This study focused on estimating vertical rates of exchange across the sediment?water interface by direct measurement using seepage meters and by using temperature as a tracer coupled with numerical modeling. Temperature loggers and pressure transducers were placed in monitoring wells within the streambed and in the river to continuously monitor temperature and hydraulic head every 15 minutes from March 2004 to October 2005. One-dimensional modeling of heat and water flow was used to interpret the temperature and head observations and deduce the sw?gw fluxes using the USGS numerical model, VS2DH, which simulates variably saturated water flow and solves the energy transport equation. Results of the modeling effort indicate that the Merced River at the study reach is generally a slightly gaining stream with small head differences (cm) between the surface water and ground water, with flow reversals occurring during high streamflow events. The average vertical flux across the sediment?water interface was 0.4?2.2 cm/day, and the range of hydraulic conductivities was 1?10 m/day. Seepage meters generally failed to provide accurate data in this high-energy system because of slow seepage rates and a moving streambed resulting in scour or burial of the seepage meters. Estimates of streambed hydraulic conductivity were also made using grain-size analysis and slug tests. Estimated hydraulic conductivity for the upstream transect determined using slug tests ranged from 40 to 250 m/day, whereas the downstream transect ranged from 10 to 100 m/day. The range in variability was a result of position along each transect. A relative percent difference was used to describe the variability in estimates of hydraulic conductivity by grain-size analysis and slug test. Variability in applied methods at the upstream transect ranged from 0 to 9 percent, whereas the downstream transect showed greater variability, with a range of 80 to 133 percent.

Scientific Investigations Report

Description and field analysis of a coupled ground-water/surface-water flow model (MODFLOW/BRANCH) with modifications for structures and wetlands in southern Dade County, Florida

A coupled surface-water model (BRANCH) and ground-water model (MODFLOW) model were tested to simulate the interacting wetlands/surface-water/ ground-water system of southern Dade County. Several options created for the MODFLOW ground- ground-water model were used in representing this field situation. The primary option is the MODBRANCH interfacing software, which allows leakage to be accounted for between the MODFLOW ground-water model and the BRANCH dynamic model for simulation of flow in an interconnected network of open channels. A modification to an existing software routine, which is referred to as BCF2, allows cells in MODFLOW to rewet when dry--a requirement in representing the seasonal wetlands in Dade County. A companion to BCF2 is the modified evapotranspiration routine EVT2. The EVT2 routine changes the cells where evapotranspiration occurs, depending on which cells are wet. The Streamlink package represents direct connections between the canals and wetlands at locations where canals open directly into overland flow. Within the BRANCH model, the capability to represent the numerous hydraulic structures, gated spillways, gated culverts, and pumps was added. The application of these modifications to model surface-water/ground-water interactions in southern Dade County demonstrated the usefulness of the coupled MODFLOW/BRANCH model. Ground-water and surface-water flows are both simulated with dynamic models. Flow exchange between models, intermittent wetting and drying, evapotranspiration, and hydraulic structure operations are all represented appropriately. Comparison was made with a simulation using the RIV1 package instead of MODBRANCH to represent the canals. RIV1 represents the canals by user-defined stages, and computes leakage to the aquifer. Greater accuracy in reproducing measured ground- water heads was achieved with MODBRANCH, which also computes dynamic flow conditions in the canals, unlike RIV1. The surface-water integrated flow and transport two-dimensional model (SWIFT2D) was also applied to the southeastern coastal wetlands for comparison with the wetlands flow approximation made in MODFLOW. MODFLOW simulates the wetlands as a highly conductive upper layer of the aquifer, whereas SWIFT2D solves the hydrodynamic equations. Comparison in this limited test demonstrated no specific advantage for either method of representation. However, much additional testing on a wider variety of geometric and hydraulic situations, such as in areas with greater tidal or other dynamic forcing effects, is needed to make definite conclusions. A submodel of the existing southern Dade County model schematization was used to examine water-delivery alternatives proposed by the U.S. Army Corps of Engineers. For this application, the coupled MODFLOW/BRANCH model was used as a design tool. A new canal and several pumps to be tested to maintain lower water levels in a residential area (while water levels in the Everglades are raised) were added to the model schematization. The pumps were assumed to have infinite supply capacity in the model so that their maximum pumping rates during the simulation could be used to determine pump sizes.

Water-Resources Investigations Report

Water-resources-related information for the Oneida Reservation and vicinity, Wisconsin

Water-resources information has been compiled from 82 studies in which data were collected from the Oneida Reservation and vicinity. Forty-seven studies addressed surface-water issues, 33 studies addressed ground-water issues, and 23 studies addressed aquatic-biology issues. Some multidisciplinary studies are included in more than one category. Most of the surface-water studies summarized in this report included both water-quality and flow information. Several surface-water studies provided detailed short-term descriptions of surface- water quality and flow for parts of the Reservation and vicinity. Surface-water and stream-sediment quality data from several data bases have been compiled for this report. Most of the compiled data come from two sites on Duck Creek. Data from Duck Creek were analyzed for trends in concentrations of suspended sediment, dissolved nitrite plus nitrate, and dissolved atrazine. No trends were detected for any of these constituents. Trends in concentration of most constituents in surface-water samples were not calculated because of the short period of data collection at nearly all of the sites. Most of the ground-water reports that were identified included both quality and quantity and flow information. None of the ground-water studies provided a detailed description of ground-water quality for the Reservation as a whole. Several reports provide varied and detailed information for ground-water models that are useful for understanding hydrogeology and ground-water flow for the Reservation and vicinity. Ground-water quality data from 180 wells, compiled from several data bases, provided an incomplete summary of the condition of the drinking- water resources of the Reservation. Only 12 constituents, from a small number of wells, exceeded a USEPA drinking-water limit. Most of the exceedences were for trace metals and organics. No exceedences for pesticides or nitrate were reported; however, pesticide data were collected from only a small number of wells. Most of the aquatic biology studies described in this report include fish data, habitat data, or calculations of biotic index values, most of which comes from Duck Creek. Historical aquatic biology data for the Reservation and vicinity are limited. Most of the 23 studies described here were done since 1992. Most of the biota-quality data compiled for this report come from several sites on Duck Creek and represent a small number of samples. Results of the community survey regarding the water resources of the Oneida Reservation indicate that water usage by Tribal members today has declined when compared to the past. The most common reason given for the decline in usage was pollution. Most of those surveyed perceived Duck Creek as being "polluted," but about 50 percent thought that water quality in the Reservation was improving.

Wisconsin

Simulation of Ground-Water Flow and Optimization of Withdrawals from Aquifers at the Naval Air Station Patuxent River, St. Mary's County, Maryland

Potentiometric surfaces in the Piney Point-Nanjemoy, Aquia, and Upper Patapsco aquifers have declined from 1950 through 2000 throughout southern Maryland. In the vicinity of Lexington Park, Maryland, the potentiometric surface in the Aquia aquifer in 2000 was as much as 170 feet below sea level, approximately 150 feet lower than estimated pre-pumping levels before 1940. At the present rate, the water levels will have declined to the regulatory allowable maximum of 80 percent of available drawdown in the Aquia aquifer by about 2050. The effect of the withdrawals from these aquifers by the Naval Air Station Patuxent River and surrounding users on the declining potentiometric surface has raised concern for future availability of ground water. Growth at Naval Air Station Patuxent River may increase withdrawals, resulting in further drawdown. A ground-water-flow model, combined with optimization modeling, was used to develop withdrawal scenarios that minimize the effects (drawdown) of hypothetical future withdrawals. A three-dimensional finite-difference ground-water-flow model was developed to simulate the ground-water-flow system in the Piney Point-Nanjemoy, Aquia, and Upper Patapsco aquifers beneath the Naval Air Station Patuxent River. Transient and steady-state conditions were simulated to give water-resource managers additional tools to manage the ground-water resources. The transient simulation, representing 1900 through 2002, showed that the magnitude of withdrawal has increased over that time, causing ground-water flow to change direction in some areas. The steady-state simulation was linked to an optimization model to determine optimal solutions to hypothetical water-management scenarios. Two optimization scenarios were evaluated. The first scenario was designed to determine the optimal pumping rates for wells screened in the Aquia aquifer within three supply groups to meet a 25-percent increase in withdrawal demands, while minimizing the drawdown at a control location. The resulting optimal solution showed that pumping six wells above the rate required for maintenance produced the least amount of drawdown in the local potentiometric surface. The second hypothetical scenario was designed to determine the optimal location for an additional well in the Aquia aquifer in the northeastern part of the main air station. The additional well was needed to meet an increase in withdrawal of 43,000 cubic feet per day. The optimization model determined the optimal location for the new well, out of a possible 10 locations, while minimizing drawdown at control nodes located outside the western boundary of the main air station. The optimal location is about 1,500 feet to the east-northeast of the existing well.

Scientific Investigations Report

Geochemical investigation of the Arbuckle-Simpson Aquifer, South-Central Oklahoma, 2004-06

A geochemical reconnaissance investigation of the Arbuckle-Simpson aquifer in south-central Oklahoma was initiated in 2004 to characterize the ground-water quality at an aquifer scale, to describe the chemical evolution of ground water as it flows from recharge areas to discharge in wells and springs, and to determine the residence time of ground water in the aquifer. Thirty-six water samples were collected from 32 wells and springs distributed across the aquifer for chemical analysis of major ions, trace elements, isotopes of oxygen and hydrogen, dissolved gases, and age-dating tracers. In general, the waters from wells and springs in the Arbuckle-Simpson aquifer are chemically suitable for all regulated uses, such as public supplies. Dissolved solids concentrations are low, with a median of 347 milligrams per liter (mg/L). Two domestic wells produced water with nitrate concentrations that exceeded the U.S. Environmental Protection Agency's nitrate maximum contaminant level (MCL) of 10 mg/L. Samples from two wells in the confined part of the aquifer exceeded the secondary maximum contaminant level (SMCL) for chloride of 250 mg/L and the SMCL of 500 mg/L for dissolved solids. Water samples from these two wells are not representative of water samples from the other wells and springs completed in the unconfined part of the aquifer. No other water samples from the Arbuckle-Simpson geochemical reconnaissance exceeded MCLs or SMCLs, although not every chemical constituent for which the U.S. Environmental Protection Agency has established a MCL or SMCL was analyzed as part of the Arbuckle-Simpson geochemical investigation. The major ion chemistry of 34 of the 36 samples indicates the water is a calcium bicarbonate or calcium magnesium bicarbonate water type. Calcium bicarbonate water type is found in the western part of the aquifer, which is predominantly limestone. Calcium magnesium bicarbonate water is found in the eastern part of the aquifer, which is predominantly a dolomite. The major ion chemistry for these 34 samples is consistent with a set of water-rock interactions. Rainfall infiltrates the soil zone, where the host rock, limestone or dolomite, dissolves as a result of uptake of carbon dioxide gas. Some continued dissolution of dolomite and precipitation of calcite occur as the water flows through the saturated zone. The major ion chemistry of the two samples from wells completed in the confined part of the aquifer indicates the water is a sodium chloride type. Geochemical inverse modeling determined that mixing of calcite-saturated recharge water with brine and dissolving calcite, dolomite, and gypsum accounts for the water composition of these two samples. One of the two samples, collected at Vendome Well in Chickasaw National Recreation Area, had a mixing fraction of brine of about 1 percent. The brine component of the sample at Vendome Well is likely to account for the relatively large concentrations of many of the trace elements (potassium, fluoride, bromide, iodide, ammonia, arsenic, boron, lithium, selenium, and strontium) measured in the water sample. Carbon-14, helium-3/tritium, and chlorofluorocarbons were used to calculate ground-water ages, recharge temperatures, and mixtures of ground water in the Arbuckle-Simpson aquifer. Thirty four of 36 water samples recharged the aquifer after 1950, indicating that water is moving quickly from recharge areas to discharge at streams and springs. Two exceptions to this classification were noted in samples 6 and 15 (Vendome Well). Ground-water ages determined for these two samples by using carbon-14 are 34,000 years (site 6) and 10,500 years (site 15). Concentrations of dissolved argon, neon, and xenon in water samples were used to determine the temperature of the water when it recharged the aquifer. The mean annual air temperature at Ada, Oklahoma, is 16 degrees Celsius (C) and the median temperature of the 30 reconnaissance water samples was 18.1 C. The av

Oklahoma

Preconstruction and simulated postconstruction ground-water levels at urban centers in the Red River Navigation Project area, Louisiana

The Red River Valley in Louisiana is 3 to 10 miles wide and is underlain by the Red River alluvial aquifer. This aquifer is in hydraulic connection with the Red River. Precipitation infiltrates the aquifer and water discharges from the aquifer at the Red River and major tributaries. Construction of locks and dams along the Red River will cause new, higher minimum stages for the river for each of the pools. Water levels will rise and outflow from the alluvial aquifer to the river after construction of the locks and dams will be at higher elevations because of the new river stages. Mathematical models of the stream-aquifer system simulate the effects of the lock and dam system for plan B-3, modified. Water-level changes in the aquifer in response to the changes imposed by the system are greatest near the river and are progressively smaller away from the river. The saturated zone in the fine-grained material overlying the aquifer provides the principal potential impact on structures and vegetation in urban areas. As water levels in the aquifer rise, the saturated zone above the aquifer will decrease the thickness. The soil zone in low areas in some communities will be waterlogged. (USGS)

Louisiana

The association of arsenic with redox conditions, depth, and ground-water age in the glacial aquifer system of the northern United States

More than 800 wells in the glacial aquifer system of the Northern United States were sampled for arsenic as part of U.S. Geological Survey National Water-Quality Assessment (NAWQA) studies during 1991-2003. Elevated arsenic concentrations (greater than or equal to 10 micrograms per liter) were detected in 9 percent of samples. Elevated arsenic concentrations were associated with strongly reducing conditions. Of the samples classified as iron reducing or sulfate reducing, arsenic concentrations were elevated in 19 percent. Of the methanogenic samples, arsenic concentrations were elevated in 45 percent. In contrast, concentrations of arsenic were elevated in only 1 percent of oxic samples. Arsenic concentrations were also related to ground-water age. Elevated arsenic concentrations were detected in 34 percent of old waters (recharged before 1953) as compared to 4 percent of young waters (recharged since 1953). For samples classified as both old and methanogenic, elevated arsenic concentrations were detected in 62 percent of samples, as compared to 1 percent for samples classified as young and oxic. Arsenic concentrations were also correlated with well depth and concentrations of several chemical constituents, including (1) constituents linked to redox processes and (2) anions or oxyanions that sorb to iron oxides. Observations from the glacial aquifer system are consistent with the idea that the predominant source of arsenic is iron oxides and the predominant mechanism for releasing arsenic to the ground water is reductive desorption or reductive dissolution. Arsenic is also released from iron oxides under oxic conditions, but on a more limited basis and at lower concentrations. Logistic regression was used to investigate the relative significance of redox, ground-water age, depth, and other water-quality constituents as indicators of elevated arsenic concentrations in the glacial aquifer system. The single variable that explained the greatest amount of variation in the data was redox. Multivariate models that included a redox variable overestimated the percentage of samples with elevated arsenic concentrations because, even though elevated arsenic concentrations were associated with strongly reducing samples, not all strongly reducing samples had elevated arsenic concentrations. Arsenic concentrations and redox conditions differed among four broad areas of the glacial aquifer system. For the East, Central, and West-Central north areas, there was a trend of increasing arsenic concentrations that corresponded to an increase in reducing conditions. For the West-Central south area, arsenic concentrations in oxic samples were higher than for the other areas, possibly because of high concentrations of orthophosphate, which is linked to desorption of arsenic from iron oxides under oxic conditions. The observed differences in arsenic concentrations among broad areas of the glacial aquifer system were generally consistent with a conceptual model developed by Smedley and Kinniburg, who studied or reviewed studies of widespread arsenic contamination in Bangladesh, India, China, Vietnam, Hungary, Argentina, northern Chile and the Southwestern United States.

Scientific Investigations Report

Hydrologic analysis of the High Plains aquifer system in Box Butte County, Nebraska

During the past 40 years, pumpage of ground water for irrigation from the High Plains aquifer system underlying Box Butte County, Nebraska, has resulted in a steady decline of water levels. Consequently, a digital model of the aquifer system was constructed to evaluate various water-management alternatives. The hydraulic conductivity of the aquifer system ranges from 6 to 60 feet per day; the specific yield ranges from 12 to 21 percent; and natural recharge ranges from 0.06 to 4.33 inches annually. Predevelopment saturated thickness (1938) ranged from 190 to 510 feet. Water pumped in 1980 was estimated at 104,000 acre-feet from an estimated recoverable volume of 34.4 million acre-feet in the aquifer system. Results from model simulation predict that the area of water-level declines of 10 feet or more will increase from 336 square miles (1981) to 630 square miles by 1991 if pumpage is increased at the maximum annual rate experienced for the period 1972-81. Maximum water-level declines would increase from 50 feet (1981) to 79 feet (1991). However, pumpage rates held at the 1981 level (no further development) would limit the decline area of 10 feet or more to 530 square miles by 1991 and the maximum decline to 63 feet.

Nebraska

Hydraulic properties and ground-water flow in the St Peter-Prairie du Chien-Jordan aquifer, Rochester area, southeastern Minnesota

The hydraulic properties were updated and their effects on ground-water flow in the St. Peter-Prairie du Chien-Jordan aquifer in the Rochester area in southeastern Minnesota were evaluated, using new information compiled since a study by Delin (1990). Since 1988, new information on the hydrogeology of the ground-water system in the Rochester area has become available from well-drilling and construction activity associated with Rochester's rapid growth. The St. Peter-Prairie du Chien-Jordan aquifer consists of the St. Peter Sandstone, the Prairie du Chien Group (limestones and dolomites), and the Jordan Sandstone. Horizontal hydraulic conductivity and transmissivity were determined from 15 aquifer tests and specific-capacity information compiled for 310 wells. A 140-square-mile area of the aquifer bounded on the west, south, and east by a ground-water divide contributes water to the Rochester, Minnesota, municipal wells. Transmissivities for the St. Peter-Prairie du Chien-Jordan aquifer in the study area range from less than 5,000 square feet per day (ft 2 /d) to greater than 20,000 ft 2 /d. Transmissivities greater than 20,000 ft 2 /d occur in the west-central, northwestern, and east-central parts of the study area. Transmissivities of less than 5,000 ft 2 /d occur in the northern, northeastern, central, and southern parts of the study area. The areas of greatest potential well yield coincide with areas of greatest transmissivity. Delin (1990) developed a ground-water-flow model to simulate flow of ground water in the St. Peter-Prairie du Chien-Jordan aquifer in the Rochester area. The 1988 Rochester model was rerun using revised horizontal hydraulic conductivity arrays in the model, based on the transmissivity distribution determined for this study. The results of the simulations using horizontal hydraulic conductivities based on the transmissivity distribution determined for this study may indicate that transmissivity values derived from specific-capacity information generally are too high. The transmissivity distribution determined for this study, however, is valid as an indicator of the spatial variability of the relative magnitude of transmissivity and potential well yield for the St. Peter-Prairie du Chien-Jordan aquifer in the study area. Water-level changes in wells from January through February 1988 to February through March 1995 ranged from -6.8 to +15.3 feet. Water-level changes in 12 Rochester municipal wells for the same period ranged from -7.4 to +8.0 feet. Water levels in wells generally rose in the northern and eastern parts of the study area and generally declined in the southwestern and western parts. Near Rochester, water levels in wells generally declined near the city boundaries and showed little change or rose in the central part of the city. Water-level changes from 1988 to 1995 near the ground-water divide generally were less than 2 feet, resulting in no appreciable changes in the location of the divide.

Minnesota