Search USGSSearch

SEARCH · Search USGS

Results for “Surface Water Techniques”

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 109 records · Page 6Linked to original sources

Guidelines for surveying bankfull channel geometry and developing regional hydraulic-geometry relations for streams of New York State

Many disturbed streams within New York State are being restored in an effort to provide bank and bed stability and thereby decrease sedimentation and erosion. Efforts to identify and provide accurate indicators for stable-channel characteristics for ungaged streams have been hampered by the lack of regional equations or relations that relate drainage area to bankfull discharge and to channel depth, width, and cross-sectional area (bankfull hydraulic-geometry relations). Regional equations are needed to confirm bankfull hydraulic-geometry, assess stream stability, evaluate restoration needs, and verify restoration design for ungaged streams that lack stage-to-discharge ratings or historic peak-flow records. This report presents guidelines for surveying bankfull channel geometry at USGS stream gages and developing regional hydraulic-geometry relations (equations) for wadeable streams in New York. It summarizes methods to (1) compile and assess existing hydrologic, geometric, photographic, and topographic data, (2) conduct stream-reconnaissance inspections, (3) identify channel-bankfull characteristics, (4) conduct longitudinal and cross-section surveys, (5) measure stream discharge, (6) develop and refine bankfull hydraulic-geometry equations, and (7) analyze and assure data completeness and quality. The techniques primarily address wadeable streams with either active or discontinued surface-water and crest-stage gages. The relations can be applied to ungaged or actively gaged streams that are wadeable, and may be extended to non-wadeable streams (with some limitations) if they have drainage areas comparable to those used to develop the relations.

Open-File Report

Integrating surface and borehole geophysics in ground water studies: An example using electromagnetic soundings in south Florida

Time domain surface electromagnetic soundings, borehole induction logs, and other borehole logging techniques are used to construct a realistic model for the shallow subsurface hydraulic properties of unconsolidated sediments in south Florida. Induction logs are used to calibrate surface induction soundings in units of pore water salinity by correlating water sample specific electrical conductivity with the electrical conductivity of the formation over the sampled interval for a two‐layered aquifer model. Geophysical logs are also used to show that a constant conductivity layer model is appropriate for the south Florida study. Several physically independent log measurements are used to quantify the dependence of formation electrical conductivity on such parameters as salinity, permeability, and clay mineral fraction. The combined interpretation of electromagnetic soundings and induction logs was verified by logging three validation boreholes, confirming quantitative estimates of formation conductivity and thickness in the upper model layer, and qualitative estimates of conductivity in the lower model layer.

Journal of Environmental & Engineering Geophysics

Preliminary Investigation of Paleochannels and Groundwater Specific Conductance using Direct-Current Resistivity and Surface-Wave Seismic Geophysical Surveys at the Standard Chlorine of Delaware, Inc., Superfund Site, Delaware City, Delaware, 2008

The U.S. Geological Survey (USGS), in cooperation with Region III of the U.S. Environmental Protection Agency (USEPA) and the State of Delaware, is conducting an ongoing study of the water-quality and hydrogeologic properties of the Columbia and Potomac aquifers and the extent of cross-aquifer contamination with benzene; chlorobenzene; 1,2-dichlorobenzene; 1,4-dichlorobenzene; and hydrogen chloride (hydrochloric acid when dissolved in water) in the vicinity of the Standard Chlorine of Delaware, Inc. (SCD), Superfund Site, Delaware City, Delaware. Surface geophysical surveys and well data were used to identify and correlate low-permeability units (clays) across the site and to search for sand and gravel filled paleochannels that are potential conduits and receptors of contaminated groundwater and (or) Dense Non-Aqueous Phase Liquid (DNAPL) contaminants. The combined surveys and well data were also used to characterize areas of the site that have groundwater with elevated (greater than 1,000 microsiemens per centimeter) specific conductance (SC) as a result of contamination. The most electrically conductive features measured with direct-current (DC) resistivity at the SCD site are relatively impermeable clays and permeable sediment that are associated with elevated SC in groundwater. Many of the resistive features include paleochannel deposits consisting of coarse-grained sediments that are unsaturated, have low (less than 200 microsiemens per centimeter) SC pore water, or are cemented. Groundwater in uncontaminated parts of the Columbia aquifer and of the Potomac aquifer has a low SC. Specific-conductance data from monitoring wells at the site were used to corroborate the DC-resistivity survey results. For comparison with DC-resistivity surveys, multi-channel analysis of surface wave (MASW) surveys were used and were able to penetrate deep enough to measure the Columbia aquifer, which is known to have elevated SC in some places. MASW survey results respond to solid material stiffness; clays and cemented sediments will have a higher velocity than silts, sands, and gravels (in order of increasing hydraulic conductivity). Geophysical surveys detected elevated SC associated with contamination of the surficial Columbia aquifer. Groundwater with elevated SC over ambient (by an order of magnitude) produced a decrease in measured resistivity at the SCD site. Where SC data are not available from wells, it is not known if a low resistivity value measured with DC resistivity alone results from the geologic material (clay) or elevated SC in groundwater (in sand or gravel). Seismic surface waves used as part of the MASW technique are not affected by water content or quality and are used herein to distinguish between sand and clay when SC is high. Through concurrent interpretation of MASW and DC-resistivity surveys, information was gained about water quality and lithology over large areas at the SCD site.

Open-File Report

Ground-water hydrology, historical water use, and simulated ground-water flow in Cretaceous-age Coastal Plain aquifers near Charleston and Florence, South Carolina

A quasi-three-dimensional, transient, digital, ground-water flow model representing the Coastal Plain aquifers of South Carolina, has been constructed to assist in defining the ground- water-flow system of Cretaceous aquifers near Charleston and Florence, S.C. Both cities are near the centers of large (greater than 150 feet) potentiometric declines in the Middendorf aquifer. In 1989, the diameter of the depressions was approximately 40 miles at Charleston and 15 miles at Florence. The potentiometric decline occurred between predevelopment (1926) and 1982 near Florence, and between predevelopment (1879) and 1989 near Charleston. The city of Charleston does not withdraw water from these aquifers; however, some of the small communities in the area use these aquifers for a potable water supply. The model simulates flow in and between four aquifer systems. The model has a variable-cell-size grid, and spans the Coastal Plain from the Savannah River in the southwest to the Cape Fear Arch in the northeast, and from the Fall Line in the northwest to approximately 30 miles offshore to the southeast. Model-grid cell size is 1 by 1 mile in a 48 by 48 mile area centered in Charleston, and in a 36 by 48 mile area centered in Florence. The model cell size gradually increases to a maximum of 4 by 4 miles outside the two study areas. The entire grid consists of 115 by 127 cells and covers an area of 39,936 square miles. The model was calibrated to historical water-level data. The calibration relied on three techniques: (1) matching simulated and observed potentiometric map surfaces, (2) statistical comparison of observed and simulated heads, and (3) comparison of observed and simulated well hydrographs. Systematic changes in model parameters showed that simulated heads are most sensitive to changes in aquifer transmissivity. Eight predictive ground-water-use scenarios were simulated for the Mount Pleasant area, which presently (1993) uses the Middendorf aquifer as a sole-source of potable water. These simulations use various combinations of spatial distribution, and injection of treated wastewater effluent for existing and future Middendorf aquifer wells.

South Carolina

Simulating the effects of ground-water withdrawals on streamflow in a precipitation-runoff model

Precipitation-runoff models are used to assess the effects of water use and management alternatives on streamflow. Often, ground-water withdrawals are a major water-use component that affect streamflow, but the ability of surface-water models to simulate ground-water withdrawals is limited. As part of a Hydrologic Simulation Program-FORTRAN (HSPF) precipitation-runoff model developed to analyze the effect of ground-water and surface-water withdrawals on streamflow in the Ipswich River in northeastern Massachusetts, an analytical technique (STRMDEPL) was developed for calculating the effects of pumped wells on streamflow. STRMDEPL is a FORTRAN program based on two analytical solutions that solve equations for ground-water flow to a well completed in a semi-infinite, homogeneous, and isotropic aquifer in direct hydraulic connection to a fully penetrating stream. One analytical method calculates unimpeded flow at the stream-aquifer boundary and the other method calculates the resistance to flow caused by semipervious streambed and streambank material. The principle of superposition is used with these analytical equations to calculate time-varying streamflow depletions due to daily pumping. The HSPF model can readily incorporate streamflow depletions caused by a well or surface-water withdrawal, or by multiple wells or surface-water withdrawals, or both, as a combined time-varying outflow demand from affected channel reaches. These demands are stored as a time series in the Watershed Data Management (WDM) file. This time-series data is read into the model as an external source used to specify flow from the first outflow gate in the reach where these withdrawals are located. Although the STRMDEPL program can be run independently of the HSPF model, an extension was developed to run this program within GenScn, a scenario generator and graphical user interface developed for use with the HSPF model. This extension requires that actual pumping rates for each well be stored in a unique WDM dataset identified by an attribute that associates each well with the model reach from which water is withdrawn. Other attributes identify the type and characteristics of the data. The interface allows users to easily add new pumping wells, delete exiting pumping wells, or change properties of the simulated aquifer or well. Development of this application enhanced the ability of the HSPF model to simulate complex water-use conditions in the Ipswich River Basin. The STRMDEPL program and the GenScn extension provide a valuable tool for water managers to evaluate the effects of pumped wells on streamflow and to test alternative water-use scenarios. Copyright ASCE 2004.

Conference Paper

Hydrogeology of, and simulation of ground-water flow In, the Pohatcong Valley, Warren County, New Jersey

A numerical ground-water-flow model was constructed to simulate ground-water flow in the Pohatcong Valley, including the area within the U.S. Environmental Protection Agency Pohatcong Valley Ground Water Contamination Site. The area is underlain by glacial till, alluvial sediments, and weathered and competent carbonate bedrock. The northwestern and southeastern valley boundaries are regional-scale thrust faults and ridges underlain by crystalline rocks. The unconsolidated sediments and weathered bedrock form a minor surficial aquifer and the carbonate rocks form a highly transmissive fractured-rock aquifer. Ground-water flow in the carbonate rocks is primarily downvalley towards the Delaware River, but the water discharges through the surficial aquifer to Pohatcong Creek under typical conditions. The hydraulic characteristics of the carbonate-rock aquifer are highly heterogeneous. Horizontal hydraulic conductivities span nearly five orders of magnitude, from 0.5 feet per day (ft/d) to 1,800 ft/d. The maximum transmissivity calculated is 37,000 feet squared per day. The horizontal hydraulic conductivities calculated from aquifer tests using public supply wells open to the Leithsville Formation and Allentown Dolomite are 34 ft/d (effective hydraulic conductivity) and 85 to 190 ft/d (minimum and maximum hydraulic conductivity, respectively, yielding a horizontal anisotropy ratio of 0.46). Stream base-flow data were used to estimate the net gain (or loss) for selected reaches on Brass Castle Creek, Shabbecong Creek, three smaller tributaries to Pohatcong Creek, and for five reaches on Pohatcong Creek. Estimated mean annual base flows for Brass Castle Creek, Pohatcong Creek at New Village, and Pohatcong Creek at Carpentersville (from correlations of partial- and continuous-record stations) are 2.4, 25, and 45 cubic feet per second (ft3/s) (10, 10, and 11 inches per year (in/yr)), respectively. Ground-water ages estimated using sulfur hexafluoride (SF6), chlorofluorocarbon (CFC), and tritium-helium age-dating techniques range from 0 to 27 years, with a median age of 6 years. Land-surface and ground-water water budgets were calculated, yielding an estimated rate of direct recharge tothe surficial aquifer of about 23 in/yr, and an estimated net recharge to the ground-water system within the area underlain by carbonate rock (11.4 mi2) of 29 in/yr (10 in/yr over the entire 33.3 mi2 basin). A finite-difference, numerical model was developed to simulate ground-water flow in the Pohatcong Valley. The four-layer model encompasses the entire carbonate-rock part of the valley. The carbonate-rock aquifer was modeled as horizontally anisotropic, with the direction of maximum transmissivity aligned with the longitudinal axis of the valley. All lateral boundaries are no-flow boundaries. Recharge was applied uniformly to the topmost active layer with additional recharge added near the lateral boundaries to represent infiltration of runoff from adjacent crystalline-rock areas. The model was calibrated to June 2001 water levels in wells completed in the carbonate-rock aquifer, August 2000 stream base-flow measurements, and the approximate ground-water age. The ground-water-flow model was constructed in part to test possible site contamination remediation alternatives. Four previously determined ground-water remediation alternatives (GW1, GW2, GW3, and GW4) were simulated. For GW1, the no-action alternative, simulated pathlines originating in the tetrachloroethene (PCE) and trichloroethene (TCE) source areas within the Ground-Water Contamination Site end at Pohatcong Creek near the confluence with Shabbecong Creek, although some particles went deeper in the aquifer system and ultimately discharge to Pohatcong Creek about 10 miles downvalley in Pohatcong Township. Remediation alternatives GW2, GW3, and GW4 include ground-water withdrawal, treatment, and reinjection. The design for GW2 includes wells in the TCE and PCE source areas that wit

New Jersey

A non-intrusive approach for efficient stochastic emulation and optimization of model-based nitrate-loading management decision support

Use of physically-motivated numerical models like groundwater flow-and-transport models for probabilistic impact assessments and optimization under uncertainty (OUU) typically incurs such a computational burdensome that these tools cannot be used during decision making. The computational challenges associated with these models can be addressed through emulation. In the land-use/water-quality context, the linear relation between nitrate loading and surface-water/groundwater nitrate concentrations presents an opportunity for employing an efficient model emulator through the application of impulse-response matrices. When paired with first-order second-moment techniques, the emulation strategy gives rise to the “stochastic impulse-response emulator” (SIRE). SIRE is shown to facilitate non-intrusive, near-real time, and risk-based evaluation of nitrate-loading change scenarios, as well as nitrate-loading OUU subject to surface-water/groundwater concentration constraints in high decision variable and parameter dimensions. Two case studies are used to demonstrate SIRE in the nitrate-loading context.

Environmental Modelling and Software

Ground-water recharge in the arid and semiarid southwestern United States

Ground-water recharge in the arid and semiarid southwestern United States results from the complex interplay of climate, geology, and vegetation across widely ranging spatial and temporal scales. Present-day recharge tends to be narrowly focused in time and space. Widespread water-table declines accompanied agricultural development during the twentieth century, demonstrating that sustainable ground-water supplies are not guaranteed when part of the extracted resource represents paleorecharge. Climatic controls on ground-water recharge range from seasonal cycles of summer monsoonal and winter frontal storms to multimillennial cycles of glacial and interglacial periods. Precipitation patterns reflect global-scale interactions among the oceans, atmosphere, and continents. Large-scale climatic influences associated with El Niño and Pacific Decadal Oscillations strongly, but irregularly, control weather in the study area, so that year-to-year variations in precipitation and ground-water recharge are large and difficult to predict. Proxy data indicate geologically recent periods of naturally occurring multidecadal droughts unlike any in the modern instrumental record. Any anthropogenically induced climate change will likely reduce ground-water recharge through diminished snowpack at higher elevations. Future changes in El Niño and monsoonal patterns, both crucial to precipitation in the study area, are highly uncertain in current models. Current land-use modifications influence ground-water recharge through vegetation, irrigation, and impermeable area. High mountain ranges bounding the study area—the San Bernadino Mountains and Sierra Nevada to the west, and the Wasatch and southern Colorado Rocky Mountains to the east—provide external geologic controls on ground-water recharge. Internal geologic controls stem from tectonic processes that led to numerous, variably connected alluvial-filled basins, exposure of extensive Paleozoic aquifers in mountainous recharge areas, and distinct modes of recharge in the Colorado Plateau and Basin and Range subregions. The chapters in this professional paper present (first) an overview of climatic and hydrogeologic framework (chapter A), followed by a regional analysis of ground-water recharge across the entire study area (chapter B). These are followed by an overview of site-specific case studies representing different subareas of the geographically diverse arid and semiarid southwestern United States (chapter C); the case studies themselves follow in chapters D–K. The regional analysis includes detailed hydrologic modeling within the framework of a high-resolution geographic-information system (GIS). Results from the regional analysis are used to explore both the distribution of ground-water recharge for mean climatic conditions as well as the influence of two climatic patterns—the El Niño-Southern Oscillation and Pacific Decadal Oscillation—that impart a high degree of variability to the hydrologic cycle. Individual case studies employ a variety of geophysical and geochemical techniques to investigate recharge processes and relate the processes to local geologic and climatic conditions. All of the case studies made use of naturally occurring tracers to quantify recharge. Thermal and geophysical techniques that were developed in the course of the studies are presented in appendices. The quantification of ground-water recharge in arid settings is inherently difficult due to the generally low amount of recharge, its spatially and temporally spotty nature, and the absence of techniques for directly measuring fluxes entering the saturated zone from the unsaturated zone. Deep water tables in arid alluvial basins correspond to thick unsaturated zones that produce up to millennial time lags between changes in hydrologic conditions at the land surface and subsequent changes in recharge to underlying ground water. Recent advances in physical, chemical, isotopic, and modeling techniques have fostered new types of recharge assessments. Chemical and isotopic techniques include an increasing variety of environmental tracers that are useful and robust. Physically based techniques include the use of heat as a tracer and computationally intensive geophysical imaging tools for characterizing hydrologic conditions in the unsaturated zone. Modeling-based techniques include spatially distributed water-budget computations using high-resolution remotely sensed and ground-based geographic data. Application of these techniques to arid and semiarid settings in the southwestern United States reveals distinct patterns of recharge corresponding to geologic setting, climatic and vegetative history, and land use. Analysis of recharge patterns shows that large expanses of alluvial basin floors are drying out under current climatic conditions, with little to no recharge to underlying ground water. Ground-water recharge occurs mainly beneath upland catchments in which thin soils overlie permeable bedrock, ephemeral channels in which flow may average only several hours per year, and active agricultural areas. The chapters in this professional paper represent a coordinated attempt to develop a better understanding of one of the Nation's most critical yet difficult-to-quantify renewable resources.

Professional Paper

Evaluating the sources of water to wells: Three techniques for metamodeling of a groundwater flow model

For decision support, the insights and predictive power of numerical process models can be hampered by insufficient expertise and computational resources required to evaluate system response to new stresses. An alternative is to emulate the process model with a statistical “metamodel.” Built on a dataset of collocated numerical model input and output, a groundwater flow model was emulated using a Bayesian Network, an Artificial neural network, and a Gradient Boosted Regression Tree. The response of interest was surface water depletion expressed as the source of water-to-wells. The results have application for managing allocation of groundwater. Each technique was tuned using cross validation and further evaluated using a held-out dataset. A numerical MODFLOW-USG model of the Lake Michigan Basin, USA, was used for the evaluation. The performance and interpretability of each technique was compared pointing to advantages of each technique. The metamodel can extend to unmodeled areas.

Illinois, Indiana, Michigan, Ohio, Wisconsin

Comparison of gas chromatography/mass spectrometry and immunoassay techniques on concentrations of atrazine in storm runoff

Gas chromatography/mass spectrometry (GC/MS) and enzyme-linked immunosorbent assay (ELISA) techniques were used to measure concentrations of dissolved atrazine in 149 surface-water samples. Samples were collected during May 1992–September 1993 near the mouth of the White River (Indiana) and in two small tributaries of the river. GC/MS was performed on a Hewlett-Packard 5971 A, with electron impact ionization and selected ion monitoring of filtered water samples extracted by C-18 solid phase extraction; ELISA was performed with a magnetic-particle-based assay with photometric analysis. ELISA results compared reasonably well to GC/MS measurements at concentrations below the Maximum Contaminant Level for drinking water set by the U.S. Environmental Protection Agency (3.0 μg/L), but a systematic negative bias was observed at higher concentrations. When higher concentration samples were diluted into the linear range of calibration, the relation improved. A slight positive bias was seen in all of the ELISA data compared to the GC/MS results, and the bias could be partially explained by correcting the ELISA data for cross reactivity with other triazine herbicides. The highest concentrations of atrazine were found during the first major runoff event after the atrazine was applied. Concentrations decreased throughout the rest of the sampling period even though large runoff events occurred during this time, indicating that most atrazine loading to surface waters in the study area occurs within a few weeks after application. Use of brand names is for identification purposes only and does not constitute endorsement by the U.S. Geological Survey, the Uniroyal Chemical Company, or Wichita State University.

Indiana

Section 3. The SPARROW Surface Water-Quality Model—Theory, application and user documentation

SPARROW (SPAtially Referenced Regressions On Watershed attributes) is a watershed modeling technique for relating water-quality measurements made at a network of monitoring stations to attributes of the watersheds containing the stations. The core of the model consists of a nonlinear regression equation describing the non-conservative transport of contaminants from point and diffuse sources on land to rivers and through the stream and river network. The model predicts contaminant flux, concentration, and yield in streams and has been used to evaluate alternative hypotheses about the important contaminant sources and watershed properties that control transport over large spatial scales. This report provides documentation for the SPARROW modeling technique and computer software to guide users in constructing and applying basic SPARROW models. The documentation gives details of the SPARROW software, including the input data and installation requirements, and guidance in the specification, calibration, and application of basic SPARROW models, as well as descriptions of the model output and its interpretation. The documentation is intended for both researchers and water-resource managers with interest in using the results of existing models and developing and applying new SPARROW models. The documentation of the model is presented in two parts. Part 1 provides a theoretical and practical introduction to SPARROW modeling techniques, which includes a discussion of the objectives, conceptual attributes, and model infrastructure of SPARROW. Part 1 also includes background on the commonly used model specifications and the methods for estimating and evaluating parameters, evaluating model fit, and generating water-quality predictions and measures of uncertainty. Part 2 provides a user's guide to SPARROW, which includes a discussion of the software architecture and details of the model input requirements and output files, graphs, and maps. The text documentation and computer software are available on the Web at http://usgs.er.gov/sparrow/sparrow-mod.html .

Techniques and Methods

Guidelines for field-measured water-quality properties

The “National Field Manual for the Collection of Water-Quality Data” (NFM) provides guidelines and procedures for U.S. Geological Survey (USGS) personnel who collect data used to assess the quality of the Nation’s surface-water and groundwater resources. This chapter, NFM A6.0, provides guidance and protocols for the measurement of field parameters on site, which include the selection of sites and methods for measurement in groundwater and surface water and procedures for measurement and reporting. It updates and supersedes USGS Techniques of Water-Resources Investigations, book 9, chapter A6.0, version 2.0, by Franceska D. Wilde. Field parameters are routinely measured when water samples are collected, are often measured continually at USGS streamgages, and are regularly measured during laboratory and field experiments. The field methods for measuring field parameters described in this chapter are applicable to most natural waters. Before 2017, the NFM chapters were released in the USGS Techniques of Water-Resources Investigations series. Effective in 2018, new and revised NFM chapters are being released in the USGS Techniques and Methods series; this series change does not affect the content and format of the NFM. More information is in the general introduction to the NFM (USGS Techniques and Methods, book 9, chapter A0) at https://doi.org/10.3133/tm9A0 . The authoritative current versions of NFM chapters are available in the USGS Publications Warehouse at https://pubs.er.usgs.gov/ . Comments, questions, and suggestions related to the NFM can be addressed to nfm@usgs.gov .

Techniques and Methods

Specific electrical conductance

Electrical conductance is a measure of the capacity of a substance to conduct an electrical current. The specific electrical conductance (conductivity) of water is a function of the types and quantities of dissolved substances it contains, normalized to a unit length and unit cross section at a specified temperature. This section of the National Field Manual (NFM) describes U.S. Geological Survey (USGS) guidance and protocols for measurement of conductivity in ground and surface waters.

Techniques of Water-Resources Investigations

Hydrogeology of the area near the J4 test cell, Arnold Air Force Base, Tennessee

The U.S. Air Force operates a major aerospace systems testing facility at Arnold Engineering Development Center (AEDC) in Coffee County, Tennessee. Dewatering operations at one of the test facilities, the J4 test cell, has affected the local ground-water hydrology. The J4 test cell is approximately 100 feet in diameter, extends approximately 250 feet below land surface, and penetrates several aquifers. Ground water is pumped continuously from around the test cell to keep the cell structurally intact. Because of the test cell's depth, dewatering has depressed water levels in the aquifers surrounding the site. The depressions that have developed exhibit anisotropy that is controlled by zones of high permeability in the aquifers. Additionally, contaminants - predominately volatile organic compounds - are present in the ground-water discharge from the test cell and in ground water at several other Installation Restoration Program (IRP) sites within the AEDC facility. The dewatering activities at J4 are drawing these contaminants from the nearby sites. The effects of dewatering at the J4 test cell were investigated by studying the lithologic and hydraulic characteristics of the aquifers, investigating the anisotropy and zones of secondary permeability using geophysical techniques, mapping the potentiometric surfaces of the underlying aquifers, and developing a conceptual model of the ground-water-flow system local to the test cell. Contour maps of the potentiometric surfaces in the shallow, Manchester, and Fort Payne aquifers (collectively, part of the Highland Rim aquifer system) show anisotropic water-level depressions centered on the J4 test cell. This anisotropy is the result of features of high permeability such as chert-gravel zones in the regolith and fractures, joints, and bedding planes in the bedrock. The presence of these features of high permeability in the Manchester aquifer results in complex flow patterns in the Highland Rim aquifers near the J4 test cell. The occurrence, distribution, and orientation of these features has a great effect on ground-water flow to the J4 test cell. The depression caused by dewatering extends out horizontally through the aquifers along the most permeable pathways. Since the aquifers above the Chattanooga Shale are not separated by distinct confining units, areas in adjacent aquifers above and below these zones of high permeability in the Manchester aquifer are also dewatered. Conditions in all Highland Rim aquifers approximate steady-state equilibrium because ground-water withdrawal at the test cell has been continuous since the late 1960's. The average ground-water discharge from the dewatering system at the J4 test cell was 105 gallons per minute, for 1992-95. The ground-water capture areas in each aquifer extend into all or parts of landfill #2 and leaching pit #2 (IRP site 1), the main testing area (IRP site 7), and the old fire training area (IRP site 10). IRP sites 8 and 12 are outside the ground-water capture areas. Of the 35 sampled wells in the J4 area, 10 produced water samples containing chlorinated organic compounds such as 1,2-dichloroethane (1,2-DCA), 1,1-dichloroethylene (1,1-DCE), and trichloroethylene (TCE) in concentrations which exceeded the Tennessee Department of Environment and Conservation Maximum Contaminant Levels (MCL's) for public water-supply systems. The highest concentrations were detected in samples from well AEDC-274 with 45 micrograms per liter (mg/L) 1,2-DCA, 320 mg/L 1,1-DCE, and 1,200 mg/L TCE. These compounds are synthetic and do not occur naturally in the environment. A sample of the ground-water discharge from the J4 test cell also contained concentrations of these compounds that exceed MCL's. Chlorinated organic compounds, including 1,2-DCA; 1,1-DCE; and TCE also have been detected at IRP sites 1, 7, 8, nd 10. The six dewatering wells surrounding the J4 test cell penetrate the Chattanooga Shale and are open to the Highland Rim aquifer system, there

Water-Resources Investigations Report

Hydrogeology of the regional aquifer near Flagstaff, Arizona, 1994-97

Sandstones, siltstones, and limestones that are Pennsylvanian to Permian in age underlie the southern part of the Colorado Plateau near Flagstaff, Arizona, and contain a complex regional aquifer that has become increasingly important as a source of water for domestic, municipal, and recreational uses. Ground-water flow in the regional aquifer is poorly understood in this area because (1) depth of the aquifer limits exploratory drilling and testing and (2) the geologic structure increases the complexity of the aquifer characteristics and the ground-water flow system. Four methods were used to improve the understanding of the hydrogeology of the regional aquifer near Flagstaff. • Remote-sensing techniques and geologic mapping provided data to identify many structural features that indicate a more complex structural environment and history than previously realized. • Data from surface-geophysical techniques that included ground-penetrating radar, seismic reflection and seismic refraction, and square-array resistivity, verified that some of the geologic structures expressed at land surface propagate deep into the subsurface and through the principal water-bearing zones of the regional aquifer at near-vertical angles. • A well and spring inventory, borehole-geophysical methods, and well and aquifer tests provided additional information relating aquifer and ground-water flow characteristics to geologic structure. • Water-chemistry data, which included major ion, nutrient, trace-element, and radioactive and stable-isotope analyses, provided an independent means of verifying the hydrogeologic characteristics of the aquifer and were used to determine recharge and discharge areas, groundwater movement, and ground-water age. Ground-water recharge occurs throughout the area but is greatest at higher altitudes where precipitation is greater and in areas where heavily fractured rock units of the aquifer are exposed. The estimated annual average recharge to the regional aquifer in the study area is about 290,000 acre-feet. Ground water flows laterally and vertically through pore spaces in the rock and along faults and other fractures from high-altitude areas in the southern part of the study area to regional drains north of the study area along the Little Colorado and Colorado Rivers, and to drains south of the study area along Oak Creek and the Verde Valley. Ground-water discharge in these areas—about 400,000 acre-feet per year—exceeds the annual recharge to the aquifer in the Flagstaff area, but ground water from areas outside the study area contributes to this discharge as well. The saturated thickness of the regional aquifer averages about 1,200 feet, and the amount of water in storage could be as much as 4,800,000 acrefeet, or about 10 percent of the total volume of the aquifer. The quality of water in the regional aquifer in terms of dissolved-solids concentrations is good for most uses throughout the area. Dissolvedsolids concentrations generally are less than 500 milligrams per liter. Water in the regional aquifer is primarily a calcium magnesium bicarbonate type. In some areas near the Rio de Flag, the water has significant nitrate and chloride components, which indicate direct recharge in these areas from the Rio de Flag. Oxygen and deuterium data indicate a common recharge source for water in the aquifer and that some sites receive recharge from surface waters where evaporation has occurred. Estimated carbon-14 ages and tritium activities indicate ground-water ages from less than 200 years in the Lake Mary area to more than 5,000 years in the Wupatki area. The regional aquifer is heterogeneous and anisotrophic and has a complex ground-water flow system. The most productive water-bearing material tends to be fine- to medium-grained sandstones, and ground-water flow and potential well yields are related to geologic structure. Fracturing associated with structural deformation increases recharge locally and also increases the potential for high well yields. Surface-geophysical techniques provided information on the orientation of high-angle, deep-seated structure in the saturated zone. Borehole-geophysical data identified horizontal to near-horizontal fractures as significant components of the fracture-flow system not apparent in the surface-geophysical data. Structural features that strike northwest appear to be areas that have the greatest potential for high well yields. A north-northeastward-striking structure may be just as promising, but additional data are needed to verify this relation.

Arizona

Groundwater simulation and management models for the upper Klamath Basin, Oregon and California

The upper Klamath Basin encompasses about 8,000 square miles, extending from the Cascade Range east to the Basin and Range geologic province in south-central Oregon and northern California. The geography of the basin is dominated by forested volcanic uplands separated by broad interior basins. Most of the interior basins once held broad shallow lakes and extensive wetlands, but most of these areas have been drained or otherwise modified and are now cultivated. Major parts of the interior basins are managed as wildlife refuges, primarily for migratory waterfowl. The permeable volcanic bedrock of the upper Klamath Basin hosts a substantial regional groundwater system that provides much of the flow to major streams and lakes that, in turn, provide water for wildlife habitat and are the principal source of irrigation water for the basin's agricultural economy. Increased allocation of surface water for endangered species in the past decade has resulted in increased groundwater pumping and growing interest in the use of groundwater for irrigation. The potential effects of increased groundwater pumping on groundwater levels and discharge to springs and streams has caused concern among groundwater users, wildlife and Tribal interests, and State and Federal resource managers. To provide information on the potential impacts of increased groundwater development and to aid in the development of a groundwater management strategy, the U.S. Geological Survey, in collaboration with the Oregon Water Resources Department and the Bureau of Reclamation, has developed a groundwater model that can simulate the response of the hydrologic system to these new stresses. The groundwater model was developed using the U.S. Geological Survey MODFLOW finite-difference modeling code and calibrated using inverse methods to transient conditions from 1989 through 2004 with quarterly stress periods. Groundwater recharge and agricultural and municipal pumping are specified for each stress period. All major streams and most major tributaries for which a substantial part of the flow comes from groundwater discharge are included in the model. Groundwater discharge to agricultural drains, evapotranspiration from aquifers in areas of shallow groundwater, and groundwater flow to and from adjacent basins also are simulated in key areas. The model has the capability to calculate the effects of pumping and other external stresses on groundwater levels, discharge to streams, and other boundary fluxes, such as discharge to drains. Historical data indicate that the groundwater system in the upper Klamath Basin fluctuates in response to decadal climate cycles, with groundwater levels and spring flows rising and declining in response to wet and dry periods. Data also show that groundwater levels fluctuate seasonally and interannually in response to groundwater pumping. The most prominent response is to the marked increase in groundwater pumping starting in 2001. The calibrated model is able to simulate observed decadal-scale climate-driven fluctuations in the groundwater system as well as observed shorter-term pumping-related fluctuations. Example model simulations show that the timing and location of the effects of groundwater pumping vary markedly depending on the pumping location. Pumping from wells close (within a few miles) to groundwater discharge features, such as springs, drains, and certain streams, can affect those features within weeks or months of the onset of pumping, and the impacts can be essentially fully manifested in several years. Simulations indicate that seasonal variations in pumping rates are buffered by the groundwater system, and peak impacts are closer to mean annual pumping rates than to instantaneous rates. Thus, pumping effects are, to a large degree, spread out over the entire year. When pumping locations are distant (more than several miles) from discharge features, the effects take many years or decades to fully impact those features, and much of the pumped water comes from groundwater storage over a broad geographic area even after two decades. Moreover, because the effects are spread out over a broad area, the impacts to individual features are much smaller than in the case of nearby pumping. Simulations show that the discharge features most affected by pumping in the area of the Bureau of Reclamation's Klamath Irrigation Project are agricultural drains, and impacts to other surface-water features are small in comparison. A groundwater management model was developed that uses techniques of constrained optimization along with the groundwater flow model to identify the optimal strategy to meet water user needs while not violating defined constraints on impacts to groundwater levels and streamflows. The coupled groundwater simulation-optimization models were formulated to help identify strategies to meet water demand in the upper Klamath Basin. The models maximize groundwater pumping while simultaneously keeping the detrimental impacts of pumping on groundwater levels and groundwater discharge within prescribed limits. Total groundwater withdrawals were calculated under alternative constraints for drawdown, reductions in groundwater discharge to surface water, and water demand to understand the potential benefits and limitations for groundwater development in the upper Klamath Basin. The simulation-optimization model for the upper Klamath Basin provides an improved understanding of how the groundwater and surface-water system responds to sustained groundwater pumping within the Bureau of Reclamation's Klamath Project. Optimization model results demonstrate that a certain amount of supplemental groundwater pumping can occur without exceeding defined limits on drawdown and stream capture. The results of the different applications of the model demonstrate the importance of identifying constraint limits in order to better define the amount and distribution of groundwater withdrawal that is sustainable.

Oregon;California

Numerical-simulation and conjunctive-management models of the Hunt-Annaquatucket-Pettaquamscutt stream-aquifer system, Rhode Island

Numerical-simulation and optimization techniques were used to evaluate alternatives for the conjunctive management of ground- and surface-water resources of the Hunt­-Annaquatucket­-Pettaquamscutt stream-aquifer system in central Rhode Island. Ground-water withdrawals from the Hunt-Annaquatucket-Pettaquamscutt aquifer exceeded 8 million gallons per day during months of peak water use during 1993­98, and additional withdrawals have been proposed to meet growing demands from within and outside of the system boundary. The system is defined by the Hunt­-Annaquatucket­-Pettaquamscutt aquifer, which is composed of glacial stratified deposits, and the network of rivers, brooks, and ponds that overlie and are in hydraulic connection with the aquifer. Nearly all of the water withdrawn, however, is derived from depletions of flow in the rivers, brooks, and ponds that overlie the aquifer. Streamflow depletions are of concern to environmental agencies because of the adverse effects that reductions in streamflow can have on aquatic and riparian ecosystems. A conjunctive-management model of the stream-aquifer system was developed to simultaneously address the water-demand and streamflow-depletion issues. The objective of the model was to maximize total ground-water withdrawal from the aquifer during July, August, and September. These three months are generally the time of year when water-supply demands are largest and streamflows are simultaneously lowest. Total withdrawal from the aquifer was limited by a set of constraints specified in the model. These constraints were (1) maximum rates of streamflow depletion in the Hunt, Annaquatucket, and Pettaquamscutt Rivers; (2) minimum monthly water demands of each of three water-supply systems that withdraw water from the aquifer; and (3) minimum and maximum withdrawal rates at each supply well. The conjunctive-management model was formulated mathematically as a linear program. The model was solved by a response-matrix technique that incorporates the results of transient, numerical simulation of the stream-aquifer system into the constraint set of the linear program. The basis of the technique was the assumption that streamflow-depletion rates in each river were a linear function of ground-water-withdrawal rates at each well. This assumption was shown to be valid for the conditions evaluated in this study, primarily because of the very high transmissivity of the aquifer near many of the wells pumped for water supply. A transient, numerical model of the system was developed to simulate an average annual cycle of monthly withdrawal and hydrologic conditions representative of the 56-year period 1941­96. The transient model was used to generate characteristic streamflow-depletion responses in each river to simulated withdrawals at each well; these characteristic responses, or response coefficients, were then incorporated directly into the streamflow-depletion constraints of the linear program. Four sets of applications of the conjunctive-management model were made to determine whether total ground-water withdrawal from the aquifer during July, August, and September could be increased over the current total withdrawal for alternative definitions of the maximum rates of streamflow depletion allowed in the Hunt, Annaquatucket, and Pettaquamscutt Rivers. Current conditions were defined as the average monthly withdrawal rates at each supply well, water demands of each of the three water-supply systems, and estimated streamflow-depletion rates during the 6-year period 1993­98. Total withdrawal from all wells in the system from July through September during 1993­98 was 506.5 million gallons. Estimated streamflow-depletion rates for 1993­98 were calculated by use of the transient model, with the 1993­98 average monthly withdrawal rates specified at each supply well. Streamflow-depletion rates calculated for July, August, and September averaged 25 percent of the model-calculated pre-withdrawal streamflow rates for the Hunt River, 19 percent for the Annaquatucket River, and 7 percent for the Pettaquamscutt River. The first set of applications of the model were made with the current estimated rates of streamflow depletion in the Hunt, Annaquatucket, and Pettaquamscutt Rivers. Results of these applications indicated that total withdrawal from the aquifer during July, August, and September could be increased from about 8 to 18 percent (from 546.0 to 596.3 million gallons) over the current total withdrawal. The increased withdrawal would require modifications to the current annual withdrawal schedule of each supply well and, for the 18-percent increase, a modified network of supply wells that would include two new wells in the Annaquatucket River Basin. A second set of model applications then was made to determine if current estimated rates of streamflow depletion in the Hunt River could be reduced without increasing current estimated rates of streamflow depletion in the Annaquatucket or Pettaquamscutt Rivers. Decreases in the current rates of streamflow depletion in the Hunt River would result in increased streamflow in the river during these three months. Results showed that current rates of streamflow depletion in the Hunt River during July, August, and September could be decreased from 5 to 15 percent, depending on whether the existing or modified well network was used. Subsequent model applications indicated that substantial increases in total ground-water withdrawal from the aquifer are possible, but would require increased rates of streamflow depletion in the Annaquatucket and Pettaquamscutt Rivers. Maximum increases in the July through September withdrawal from the aquifer of about 39 to 50 percent (from 705.1 to 760.3 million gallons) over the current total withdrawal were calculated when streamflow-depletion rates in the Annaquatucket and Pettaquamscutt Rivers were allowed to increase from current estimated rates to a maximum of 25 percent of the model-calculated pre-withdrawal streamflow for each river during July, August, and September. Alternatively, it was shown that current estimated rates of streamflow depletion in the Hunt River during July, August, and September could be reduced by as much as 35 percent for the maximum allowed increases in streamflow depletion in the Annaquatucket and Pettaquamscutt Rivers; maximum increased withdrawal from the aquifer, however, would range from 8 to 18 percent over the current total withdrawal for the 35-percent reduction in streamflow-depletion rates in the Hunt River. Results of the different applications of the model demonstrate the usefulness of coupling numerical-simulation and optimization techniques for regional-scale evaluation of water-resource management alternatives. The results of the evaluation must be viewed, however, within the limitations of the quality of data available for the Hunt­-Annaquatucket­-Pettaquamscutt stream-aquifer system and representation of the system by a simulation model. An additional limitation of the analysis was the use of an average annual cycle of monthly withdrawal and hydrologic conditions. Ground-water withdrawal strategies may need to be modified to meet streamflow-depletion constraints during extreme hydrologic events, such as droughts. Contributing areas and sources of water to the supply wells also were delineated by use of a steady-state model of the stream-aquifer system. The model was developed to simulate long-term-average ground-water flow and ground-water/ surface-water interactions in the system during the 56-year period 1941­96. Sources of water to the wells consisted of precipitation and wastewater recharge to the aquifer, streamflow leakage from natural stream-channel losses, streamflow leakage caused by induced infiltration, and lateral ground-water inflow from till and bedrock upland areas.

Rhode Island

Effects of surface-water and groundwater inflows and outflows on the hydrology of the Tsala Apopka Lake Basin in Citrus County, Florida

The U.S. Geological Survey, in cooperation with the Southwest Florida Water Management District, initiated a study to quantify the inflows and outflows in the Floral City, Inverness, and Hernando pools of the Tsala Apopka Lake Basin in Citrus County, Florida. This study assesses hydrologic changes in pool stages, groundwater levels, spring flows, and streamflows caused by the diversion of streamflow from the Withlacoochee River to the Tsala Apopka Lake Basin through water-control structures. A surface-water/groundwater flow model was developed using hydraulic parameters for lakes, streams, the unsaturated zone, and the underlying surficial and Upper Floridan aquifers estimated using an inverse modeling calibration technique. After calibration, the model was used to assess the relation between inflows and outflows in the Tsala Apopka Lake Basin and changes in pool stages. Simulation results using the calibrated surface-water/groundwater flow model showed that leakage rates from the pools to the Upper Floridan aquifer were largest at the deep lake cells and that these leakage rates to the Upper Floridan aquifer were the highest in the model area. Downward leakage to the Upper Floridan aquifer occurred beneath most of the extent of the Floral City, Inverness, and Hernando pools. These leakage rates depended on the lakebed leakance and the difference between lake stages and heads in the Upper Floridan aquifer. Leakage rates were higher for the Floral City pool than for the Inverness pool, and higher for the Inverness pool than for the Hernando pool. Lakebed leakance was higher for the Floral City pool than for the Hernando pool, and higher for the Hernando pool than for the Inverness pool. Simulation results showed that the average recharge rate to the surficial aquifer was 10.3 inches per year for the 2004 to 2012 simulation period. Areas that recharge the surficial aquifer covered about 86 percent of the model area. Simulations identified areas along segments of the Withlacoochee River and within land-surface depressions that receive water from the surficial aquifer. Recharge rates were largest in physiographic regions having a deep water table. Simulated heads in the Upper Floridan aquifer indicated the general flow directions in the active flow model area were from the northeast toward the southwest and then westward toward the coast, and from the southeast toward the northwest and then westward toward the coast, consistent with flow directions inferred from the estimated potentiometric surface map for May 2010. The largest inflow in the water budget of the Upper Floridan aquifer was downward leakage from the overlying hydrogeologic unit. The largest outflow in the water budget of the Upper Floridan aquifer was spring flow. The calibrated surface-water and groundwater flow model was used to simulate hydrologic scenarios that included changes in rainfall rates, projected increases in groundwater pumping rates for 2025 and 2035, no flow for the 2004–12 period through the eight water-control structures in the Tsala Apopka Lake Basin, and the removal of the Inglis Dam and the Inglis Bypass Spillway on Lake Rousseau. Scenario simulation results were compared to annual average calibrated water levels and flows from 2004 to 2012. Simulated declines in the Tsala Apopka Lake pool stages under the 10-percent lower rainfall scenario were about 0.8, 0.3, and 1.3 feet (ft) for the Floral City, Inverness, and Hernando pools, respectively. Simulated groundwater levels under the same scenario declined up to 5.4 ft in the surficial aquifer and up to 2.9 ft in the Upper Floridan aquifer. Under the projected increases in groundwater pumping rates for 2035 that represented an increase of 36 percent from average 2004 to 2012 pumping rates, the simulated declines in the Floral City, Inverness, and Hernando pool stages were, in downstream order, 0.02, 0.06, and 0.04 ft. The largest drawdown under the projected increases in groundwater pumping rates for 2035 was 2.1 ft in the surficial aquifer and about 1.8 ft in the Upper Floridan aquifer. A scenario of decreased rainfall by 10 percent caused greater declines in water levels and pool stages than projected increases in groundwater pumping rates. The simulation with no flow through the eight Tsala Apopka Lake water-control structures resulted in simulated declines in average pool stage of 1.8, 1.9, and 0.5 ft in the Floral City, Inverness, and Hernando pools, respectively. The simulated removal of the two water-control structures in Lake Rousseau caused flow to increase at Rainbow Springs by 28 cubic feet per second, an increase of 4.7 percent from the average calibrated flow for 2004 to 2012.

Florida