SEARCH · Geology Reports

Results for “Modeling Report”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,369 records · Page 76Linked to original sources

USGS2004

Compilation of Data to Support Development of a Pesticide Management Plan by the Yankton Sioux Tribe, Charles Mix County, South Dakota

The U.S. Environmental Protection Agency is working with the Yankton Sioux Tribe to develop a pesticide management plan to reduce potential for contamination of ground water that may result from the use of registered pesticides. The purpose of this study was to compile technical information to support development of a pesticide management plan by the Yankton Sioux Tribe for the area within the Yankton Sioux Reservation, Charles Mix County, South Dakota. Five pesticides (alachlor, atrazine, cyanazine, metolachlor, and simazine) were selected by the U.S. Environmental Protection Agency for the management plan approach because they had been identified as probable or possible human carcinogens and they often had been associated with ground-water contamination in many areas and at high concentrations. This report provides a compilation of data to support development of a pesticide management plan. Available data sets are summarized in the text of this report, and actual data sets are provided in one Compact Disk?Read-Only Memory that is included with the report. The compact disk contains data sets pertinent to the development of a pesticide management plan. Pesticide use for the study area is described using information from state and national databases. Within South Dakota, pesticides commonly are applied to corn and soybean crops, which are the primary row crops grown in the study area. Water-quality analyses for pesticides are summarized for several surface-water sites. Pesticide concentrations in most samples were found to be below minimum reporting levels. Topographic data are presented in the form of 30-meter digital elevation model grids and delineation of drainage basins. Geohydrologic data are provided for the surficial deposits and the bedrock units. A high-resolution (30-by-30 meters) land-cover and land-use database is provided and summarized in a tabular format. More than 91 percent of the study area is used for row crops, pasture, or hay, and almost 6 percent of the study area is covered by water or wetlands. Average monthly and yearly precipitation data are summarized in a tabular format. Irrigation information associated with permitted and licensed diversion points is provided. A composite of aerial photographs of Charles Mix County is provided. This report also describes and summarizes the data sets and files, and how the data are relevant to development of a pesticide management plan.

Open-File Report↗

USGS1978

Summary appraisals of the nation's ground-water resources – Great Lakes region

The Great Lakes Regions, as a whole, has abundant supplies of water. Nearly 805,000 billion cubic feet of water is contained in the Great Lakes. An additional 35,000 billion cubic feet of potable ground water is available from storage in the region. Estimated ground-water discharge to the streams and lakes of the region is 26 billion gallons per day. Despite this abundance of water, the United States part of the Great Lakes basin is faced with many water-related problems, most of which involve water quality and water supply. Other problems concern periods of low flow in streams, preservation of wetlands, detrimental effects of erosion, and flooding. The significance of ground water in these problems is often overlooked. Ground water can be an alternative to surface water as a source of supply, or it can be used conjunctively with surface water to provide flexibility in water-supply management. Ground water supplied approximately 1,800 million gallons per day of the 39,900 million gallons. per day used in the Great Lakes Region in 1970. The ground-water contribution was only 4.5 percent of the water used. Thus, ground water represents a potential source of supply for much of the region. It also can be used, where conditions permit, to maintain lake levels and flow in streams, to dilute poor quality surface water, and to maintain or create wetlands and ponds. In managing water resources, ground water and surface water should be considered parts of a single system. Management includes not only planning and controlling the development but also monitoring the effects of this development. Recent advances in ground-water hydrology have provided methods to resolve some of the development and management questions that formerly slowed the development of ground water. All of the States in the Great Lakes Region have some regulations to control the development or protect the quality of the ground water. These regulations, however, are not as comprehensive as those governing surface water. Future legislation could be designed to encourage the development of ground water and, at the same time, to protect the resource. Efficient development and management of ground-water resources requires a through knowledge of the system. Reports on ground water are available for about 80 percent of tbe Great Lakes Region. Most of these reports, however, are not sufficiently detailed to be useful in comprehensive planning. As ground-water development continues, quantitative groundwater studies, utilizing models as predictive tools, will enable this development to proceed in an efficient manner.

Illinois, Indiana, Michigan, Minnesota, New York, ↗

USGS1996

Computer input and output files associated with ground-water-flow simulations of the Albuquerque Basin, central New Mexico, 1901-94, with projections to 2020; (supplement one to U.S. Geological Survey Water-resources investigations report 94-4251)

Open-File Report↗

USGS1996

Computer input and output files associated with ground-water-flow simulations of the Albuquerque Basin, central New Mexico, 1901-95, with projections to 2020; (supplement three to U.S. Geological Survey Water-resources investigations report 94-4251)

This report presents the computer input files required to run the three-dimensional ground-water-flow model of the Albuquerque Basin, central New Mexico, documented in Kernodle and others (Kernodle, J.M., McAda, D.P., and Thorn, C.R., 1995, Simulation of ground-water flow in the Albuquerque Basin, central New Mexico, 1901-1994, with projections to 2020: U.S. Geological Survey Water-Resources Investigations Report 94-4251, 114 p.) and revised by Kernodle (Kernodle, J.M., 1998, Simulation of ground-water flow in the Albuquerque Basin, 1901-95, with projections to 2020 (supplement two to U.S. Geological Survey Water-Resources Investigations Report 94-4251): U.S. Geological Survey Open-File Report 96-209, 54 p.). Output files resulting from the computer simulations are included for reference.

Open-File Report↗

USGS2012

Profile measurements and data from the 2011 Optics, Acoustics, and Stress In Situ (OASIS) project at the Martha's Vineyard Coastal Observatory

This report documents data collected by the U.S. Geological Survey (USGS) for the Coastal Model Applications and Field Measurements project under the auspices of the U.S. Navy Office of Naval Research Optics, Acoustics, and Stress In Situ (OASIS) Project. The objective of the measurements was to relate optical and acoustic properties of suspended particles to changes in particle size, concentration, and vertical distribution in the bottom boundary layer near the seafloor caused by wave- and current-induced stresses. This information on the physics of particle resuspension and aggregation and light penetration and water clarity will help improve models of sediment transport, benthic primary productivity, and underwater visibility. There is well-established technology for acoustic profiling, but optical profiles are more difficult to obtain because of the rapid attenuation of light in water. A specially modified tripod with a moving arm was designed to solve this problem by moving instruments vertically in the bottom boundary layer, between the bottom and about 2 meters above the seafloor. The profiling arm was designed, built, and tested during spring and summer 2011 by a team of USGS scientists, engineers, and technicians. To accommodate power requirements and the large data files recorded by some of the optical instruments, the tripod was connected via underwater cable to the Martha's Vineyard Coastal Observatory, operated by the Woods Hole Oceanographic Institution (WHOI). This afforded real-time Internet communication with the embedded computers aboard the tripod. Instruments were mounted on the profiling arm, and additional instruments were mounted elsewhere on the tripod and nearby on the seafloor. The tripod and a small mooring for a profiling current meter were deployed on September 17, 2011, at the Martha's Vineyard Coastal Observatory 12-meter-deep underwater node about 2 kilometers south of Martha's Vineyard, Massachusetts. Divers assisted in the deployment and cleaned the instrument surfaces on the tripod approximately once per week until the tripod and current meter were recovered on October 23, 2011. There was a range of wave and current conditions during the 36-day deployment, including the distant passage of Hurricane Ophelia, several moderate wave events, and a significant local gale that generated wave heights greater than 4 meters at the 12-meter site and knocked over the tripod 3 days before it was recovered. All but one of the instruments functioned well and provided complete datasets. The details of these data and the location of files containing the best basic version of the data are described in this report.

Massachusetts↗

USGS1997

Fate and transport of metam spill in Sacramento River

A mass balance model was developed and applied to the Sacramento River in northern California during the July 1991 Sacramento River metam-sodium spill. The transport and reactions of metam-sodium, a soil fumigant, and the volatile and toxic methyl isothiocyanate (MITC) were simulated during the two-and-a-half days of movement along a 68-km stretch of river. Results from modeling were compared with field data for MITC, which is the only product measured downriver after the spill. Agreement between the simulated and measured values of MITC concentrations were found at Doney Creek (65.9 km downstream). Results illuminated the complexities and unique characteristics associated with the multiple kinetic processes of the chemical plume in the river. In particular, the photolysis of metam-sodium followed zero-order kinetics for high concentrations and first-order kinetics for low concentrations, a unique phenomenon consistent with the finding reported in a laboratory study. Concentrations of metam-sodium for transition from zeroto first-order, obtained by calibration and model sensitivity analyses, were in the same range as those in the reported laboratory results.

Journal of Environmental Engineering↗

USGS2001

MODFLOW-2000, the U.S. Geological Survey modular ground-water model: User guide to the LMT6 package, the linkage with MT3DMS for multi-species mass transport modeling

MODFLOW-2000, the newest version of MODFLOW, is a computer program that numerically solves the three-dimensional ground-water flow equation for a porous medium using a finite-difference method. MT3DMS, the successor to MT3D, is a computer program for modeling multi-species solute transport in three-dimensional ground-water systems using multiple solution techniques, including the finite-difference method, the method of characteristics (MOC), and the total-variation-diminishing (TVD) method. This report documents a new version of the Link-MT3DMS Package, which enables MODFLOW-2000 to produce the information needed by MT3DMS, and also discusses new visualization software for MT3DMS. Unlike the Link-MT3D Packages that coordinated previous versions of MODFLOW and MT3D, the new Link-MT3DMS Package requires an input file that, among other things, provides enhanced support for additional MODFLOW sink/source packages and allows list-directed (free) format for the flow model produced flow-transport link file. The report contains four parts: (a) documentation of the Link-MT3DMS Package Version 6 for MODFLOW-2000; (b) discussion of several issues related to simulation setup and input data preparation for running MT3DMS with MODFLOW-2000; (c) description of two test example problems, with comparison to results obtained using another MODFLOW-based transport program; and (d) overview of post-simulation visualization and animation using the U.S. Geological Survey?s Model Viewer.

Open-File Report↗

USGS2022

What common-garden experiments tell us about climate responses in plants

Common garden experiments are indoor or outdoor plantings of species or populations collected from multiple distinct geographic locations, grown together under shared conditions. These experiments examine a range of questions for theory and application using a variety of methods for analysis. The eight papers of this special feature comprise a cross section of contemporary approaches, summarized and synthesized here by what they tell us about the relationships between climate-related trait spectra and fitness optima. Four of the eight papers are based on field experiments in prairie, desert, Mediterranean and boreal biomes. Representative of many common garden experiments, these experiments reveal consistent evidence of traits varying with population climate provenance, but evidence of a tradeoff between growth and tolerance traits or of consistent fitness optimization at home is scant, in contrast to trait theory. Two synthesis papers highlight dominant patterns of trait divergence, including for an exotic invasive species. One theoretical paper warned that unknown kinship relationships between populations can result in the misidentification of adaptive trait divergence. A third synthesis paper formulated novel and ambitious goals for common-garden studies through including measurement of response variables at multiple levels of biological organization. The featured papers discuss multiple avenues for improving common garden studies. Genomic analysis, together with the quantification of kinship relationships, will continue to reveal the influence of environmental drivers on gene selection. Measuring a more complete set of fitness traits, especially for traits related to regeneration, will permit the development of projection models to explicitly link trait spectra, climate patterns and fitness consequences. More standardized data reporting will additionally improve abilities to synthesize findings across experiments. Testing population performance in competition with other species will produce more robust fitness comparisons between genotypes, especially for slower-growing genotypes in higher-resource environments. Adding gardens in and beyond climatic edge locations will furthermore strengthen the understanding of population failure and species exclusion. Finally, there is unrealized potential in adding ecosystem-level observations to common-garden studies that will enhance integrative analysis across scales of biological organization and scientific domains. Synthesis . With novel, creative designs, data integration and synthesis, common garden experiments will continue to advance the understanding of trait ensembles interacting with climate across scales of biological organization, provide pivotal data for global change models and guide ecological applications such as restoration of habitats for rare and climate sensitive species.

Journal of Ecology↗

USGS1980

Effect of delayed reporting of band recoveries on survival estimates

Brownie et al. (U.S. Fish and Wildl. Serv., Resource Publ. 131, 1978) presented 14 models based on an array of explicit assumptions for the study of survival in avian populations. These methods are replacing the life table methods previously used to estimate survival rates (e.g., Burnham and Anderson, J. Wildl. Manage. , 43: 356-366, 1979). The new methods allow survival or recovery rates, or both, to be constant, time-specific, or time- and age-specific. In studies to estimate survival rates for birds the data are often from recoveries of birds shot or found dead during the hunting season and reported to the Bird Banding Laboratory by sportsmen, conservation agency employees, or the general public. This note examines the bias in estimating annual survival due to a proportion of the recoveries being incorrectly reported a year late. Specifically, a few recoveries each year of, for example, adult male American Widgeon ( Anas americana ) banded in California are reported as being recovered in year i + 1 when in fact they were actually recovered the previous year i. Delayed reporting might typically be caused by people finding a band in their health clothing in the fall of the year and, being embarrassed about their failure to report the band when it was taken, report it a year late not mentioning the actual year of recovery. Heuristically, delayed reporting should bias estimated annual survival rates upwards because it appears from the data that the birds corresponding to the "delayed" recoveries actually lived an additional year.

Journal of Field Ornithology↗

USGS2024

Flood-inundation maps for the Cuyahoga River in and near Independence, Ohio, 2024

Digital flood-inundation maps for a 9.9-mile reach of the Cuyahoga River in and near Independence, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The model was calibrated to the current (2024) stage-streamflow relation (rating curve 43.0) for the USGS streamgage 04208000, Cuyahoga River at Independence, Ohio. The resulting hydraulic model was then used to compute 13 water-surface profiles for water levels (flood stages) ranging from 14.00 to 26.00 feet. The flood stages range from “action stage” to above “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used in combination with a digital elevation model derived from light detection and ranging data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood mitigation strategies and to define flood hazard areas to protect life and property, to coordinate flood response activities such as evacuations and road closures, and to aid postflood recovery efforts.

Ohio↗

USGS2016

Investigating impacts of oil and gas development on greater sage-grouse

The sagebrush ( Artemisia spp.) ecosystem is one of the largest ecosystems in western North America providing habitat for species found nowhere else. Sagebrush habitats have experienced dramatic declines since the 1950s, mostly due to anthropogenic disturbances. The greater sage-grouse ( Centrocercus urophasianus ) is a sagebrush-obligate species that has experienced population declines over the last several decades, which are attributed to a variety of disturbances including the more recent threat of oil and gas development. We developed a hierarchical, Bayesian state-space model to investigate the impacts of 2 measures of oil and gas development, and environmental and habitat conditions, on sage-grouse populations in Wyoming, USA using male lek counts from 1984 to 2008. Lek attendance of male sage-grouse declined by approximately 2.5%/year and was negatively related to oil and gas well density. We found little support for the influence of sagebrush cover and precipitation on changes in lek counts. Our results support those of other studies reporting negative impacts of oil and gas development on sage-grouse populations and our modeling approach allowed us to make inference to a longer time scale and larger spatial extent than in previous studies. In addition to sage-grouse, development may also negatively affect other sagebrush-obligate species, and active management of sagebrush habitats may be necessary to maintain some species.

Wyoming↗

USGS2006

Revised geologic cross sections of parts of the Colorado, White River, and Death Valley regional groundwater flow systems, Nevada, Utah, and Arizona

This report presents revisions to parts of seven of the ten cross sections originally published in U.S. Geological Survey Open-File Report 2006-1040. The revisions were necessary to correct errors in some of the original cross sections, and to show new parts of several sections that were extended and (or) appended to the original section profiles. Revisions were made to cross sections C-C', D-D', E-E', F-F', G-G', I-I', and J-J', and the parts of the sections revised or extended are highlighted below the sections on plate 1 by red brackets and the word "revised," or "extended." Sections not listed above, as well as the interpretive text and figures, are generally unchanged from the original report. Cross section C-C' includes revisions in the east Mormon Mountains in the east part of the section; D-D' includes revisions in the Mormon Mesa area in the east part of the section; E-E' includes revisions in the Muddy Mountains in the east part of the section; F-F' includes revisions from the Muddy Mountains to the south Virgin Mountains in the east part of the section; and J-J' includes some revisions from the east Mormon Mountains to the Virgin Mountains. The east end of G-G' was extended about 16 km from the Black Mountains to the southern Virgin Mountains, and the northern end of I-I' was extended about 45 km from the Muddy Mountains to the Mormon Mountains, and revisions were made in the Muddy Mountains part of the original section. This report contains 10 interpretive cross sections and an integrated text describing the geology of parts of the Colorado, White River, and Death Valley regional groundwater flow systems in Nevada, Utah, and Arizona. The primary purpose of the report is to provide geologic framework data for input into a numerical groundwater model. Therefore, the stratigraphic and structural summaries are written in a hydrogeologic context. The oldest rocks (basement) are Early Proterozoic metamorphic and intrusive crystalline rocks that are considered confining units because of their low permeability. Late Proterozoic to Lower Cambrian clastic units overlie the crystalline rocks and are also considered confining units within the regional flow systems. Above the clastic units are Middle Cambrian to Lower Permian carbonate rocks that are the primary aquifers in the flow systems. The Middle Cambrian to Lower Permian carbonate rocks are overlain by a sequence of mainly clastic rocks of late Paleozoic to Mesozoic age that are mostly considered confining units, but they may be permeable where faulted. Tertiary volcanic and plutonic rocks are exposed in the northern and southern parts of the study area. In the Clover and Delamar Mountains, these rocks are highly deformed by north- and northwest-striking normal and strike-slip faults that are probably important conduits in transmitting groundwater from the basins in the northern Colorado and White River flow systems to basins in the southern part of the flow systems. The youngest rocks in the region are Tertiary to Quaternary basin-fill deposits. These rocks consist of middle to late Tertiary sediments consisting of limestone, conglomerate, sandstone, tuff, and gypsum, and younger Quaternary surficial units consisting of alluvium, colluvium, playa deposits, and eolian deposits. Basin-fill deposits are both aquifers and aquitards. The rocks in the study area were complexly deformed by episodes of Mesozoic compression and Cenozoic extensional tectonism. Some Cretaceous thrust faults and folds of the Sevier orogenic belt form duplex zones and define areas of maximum thickness for the Paleozoic carbonate rocks. Cenozoic faults are important because they are the primary structures that control groundwater flow in the regional flow systems.

Arizona, Nevada, Utah↗

USGS2005

UCODE_2005 and six other computer codes for universal sensitivity analysis, calibration, and uncertainty evaluation constructed using the JUPITER API

This report documents the computer codes UCODE_2005 and six post-processors. Together the codes can be used with existing process models to perform sensitivity analysis, data needs assessment, calibration, prediction, and uncertainty analysis. Any process model or set of models can be used; the only requirements are that models have numerical (ASCII or text only) input and output files, that the numbers in these files have sufficient significant digits, that all required models can be run from a single batch file or script, and that simulated values are continuous functions of the parameter values. Process models can include pre-processors and post-processors as well as one or more models related to the processes of interest (physical, chemical, and so on), making UCODE_2005 extremely powerful. An estimated parameter can be a quantity that appears in the input files of the process model(s), or a quantity used in an equation that produces a value that appears in the input files. In the latter situation, the equation is user-defined. UCODE_2005 can compare observations and simulated equivalents. The simulated equivalents can be any simulated value written in the process-model output files or can be calculated from simulated values with user-defined equations. The quantities can be model results, or dependent variables. For example, for ground-water models they can be heads, flows, concentrations, and so on. Prior, or direct, information on estimated parameters also can be considered. Statistics are calculated to quantify the comparison of observations and simulated equivalents, including a weighted least-squares objective function. In addition, UCODE_2005 can be used fruitfully in model calibration through its sensitivity analysis capabilities and its ability to estimate parameter values that result in the best possible fit to the observations. Parameters are estimated using nonlinear regression: a weighted least-squares objective function is minimized with respect to the parameter values using a modified Gauss-Newton method or a double-dogleg technique. Sensitivities needed for the method can be read from files produced by process models that can calculate sensitivities, such as MODFLOW-2000, or can be calculated by UCODE_2005 using a more general, but less accurate, forward- or central-difference perturbation technique. Problems resulting from inaccurate sensitivities and solutions related to the perturbation techniques are discussed in the report. Statistics are calculated and printed for use in (1) diagnosing inadequate data and identifying parameters that probably cannot be estimated; (2) evaluating estimated parameter values; and (3) evaluating how well the model represents the simulated processes. Results from UCODE_2005 and codes RESIDUAL_ANALYSIS and RESIDUAL_ANALYSIS_ADV can be used to evaluate how accurately the model represents the processes it simulates. Results from LINEAR_UNCERTAINTY can be used to quantify the uncertainty of model simulated values if the model is sufficiently linear. Results from MODEL_LINEARITY and MODEL_LINEARITY_ADV can be used to evaluate model linearity and, thereby, the accuracy of the LINEAR_UNCERTAINTY results. UCODE_2005 can also be used to calculate nonlinear confidence and predictions intervals, which quantify the uncertainty of model simulated values when the model is not linear. CORFAC_PLUS can be used to produce factors that allow intervals to account for model intrinsic nonlinearity and small-scale variations in system characteristics that are not explicitly accounted for in the model or the observation weighting. The six post-processing programs are independent of UCODE_2005 and can use the results of other programs that produce the required data-exchange files. UCODE_2005 and the other six codes are intended for use on any computer operating system. The programs consist of algorithms programmed in Fortran 90/95, which efficiently performs numerical calculations. The model runs required to obtain perturbation sensitivities can be performed using multiple processors. The programs are constructed in a modular fashion using JUPITER API conventions and modules. For example, the data-exchange files and input blocks are JUPITER API conventions and many of those used by UCODE_2005 are read or written by JUPITER API modules. UCODE-2005 includes capabilities likely to be required by many applications (programs) constructed using the JUPITER API, and can be used as a starting point for such programs.

Techniques and Methods↗

USGS2016

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967–2011, was divided into two sub-intervals because of varying water use, 1970–80 and 1981–2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981–2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967–2011; 20,724 acre-ft/yr during 1970–80, and 13,739 acre-ft/yr during 1981–2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967–2011; 19,963 acre-ft/yr during 1970–80; and 12,285 acre-ft/yr during 1981–2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967–2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980–2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980–2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7–279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4–279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981–2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994–2004) during the 1980–2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma↗

USGS1962

Theory of aquifer tests

The development of water supplies from wells was placed on a rational basis with Darcy's development of the law governing the movement of fluids through sands and with Dupuit's application of that law to the problem of radial flow toward a pumped well. As field experience increased, confidence in the applicability of quantitative methods was gained and interest in developing solutions for more complex hydrologic problems was stimulated. An important milestone was Theis' development in 1935 of a solution for the nonsteady flow of ground water, which enabled hydrologists for the first time to predict future changes in ground-water levels resulting from pumping or recharging of wells. In the quarter century since, quantitative ground-water hydrology has been enlarging so rapidly as to discourage the preparation of comprehensive textbooks. This report surveys developments in fluid mechanics that apply to groundwater hydrology. It emphasizes concepts and principles, and the delineation of limits of applicability of mathematical models for analysis of flow systems in the field. It stresses the importance of the geologic variable and its role in governing the flow regimen. The report discusses the origin, occurrence, and motion of underground water in relation to the development of terminology and analytic expressions for selected flow systems. It describes the underlying assumptions necessary for mathematical treatment of these flow systems, with particular reference to the way in which the assumptions limit the validity of the treatment.

Water Supply Paper↗

USGS2011

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗

USGS2012

Ecosystem services valuation to support decisionmaking on public lands—A case study of the San Pedro River watershed, Arizona

This report details the findings of the Bureau of Land Management–U.S. Geological Survey Ecosystem Services Valuation Pilot Study. This project evaluated alternative methods and tools that quantify and value ecosystem services, and it assessed the tools’ readiness for use in the Bureau of Land Management decisionmaking process. We tested these tools on the San Pedro River watershed in northern Sonora, Mexico, and southeast Arizona. The study area includes the San Pedro Riparian National Conservation Area (managed by the Bureau of Land Management), which has been a focal point for conservation activities and scientific research in recent decades. We applied past site-specific primary valuation studies, value transfer, the Wildlife Habitat Benefits Estimation Toolkit, and the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) and Artificial Intelligence for Ecosystem Services (ARIES) models to value locally important ecosystem services for the San Pedro River watershed—water, carbon, biodiversity, and cultural values. We tested these approaches on a series of scenarios to evaluate ecosystem service changes and the ability of the tools to accommodate scenarios. A suite of additional tools were either at too early a stage of development to run, were proprietary, or were place-specific tools inappropriate for application to the San Pedro River watershed. We described the strengths and weaknesses of these additional ecosystem service tools against a series of evaluative criteria related to their usefulness for Bureau of Land Management decisionmaking. Using these tools, we quantified gains or losses of ecosystem services under three categories of scenarios: urban growth, mesquite management, and water augmentation. These results quantify tradeoffs and could be useful for decisionmaking within Bureau of Land Management district or field offices. Results are accompanied by a relatively high level of uncertainty associated with model outputs, valuation methods, and discount rates applied. Further guidance on representing uncertainty and applying uncertain results in decisionmaking would benefit both tool developers and those offices in using ecosystem services to compare management tradeoffs. Decisionmakers and Bureau of Land Management managers at the State-, district-, and field-office level would also benefit from continuing model improvements, training, and guidance on tool use that can be provided by the U.S. Geological Survey, the Bureau of Land Management, and the Department of the Interior. Tradeoffs were identified in the level of effort needed to parameterize and run tools and the amount and quality of information they provide to the decision process. We found the Wildlife Habitat Benefits Estimation Toolkit, Ecosystem Services Review, and United Nations Environment Programme–World Conservation Monitoring Centre Ecosystem Services Toolkit to be immediately feasible for application by the Bureau of Land Management, given proper guidance on their use. It is also feasible for the Bureau of Land Management to use the InVEST model, but in early 2012 the process of parameterizing the model required resources and expertise that are unlikely to be available in most Bureau of Land Management district or field offices. Application of past primary valuation is feasible, but developing new primary-valuation studies is too time consuming for regular application. Value transfer approaches (aside from the Wildlife Habitat Benefits Estimation Toolkit) are best applied carefully on the basis of guidelines described in this report, to reduce transfer error. The ARIES model can provide useful information in regions modeled in the past (Arizona, California, Colorado, and Washington), but it lacks some features that will improve its usability, such as a generalized model that could be applied anywhere in the United States. Eleven other tools described in this report could become useful as the tools more fully develop, in high-profile cases for which additional resources are available for tool application or in case-study regions where place-specific models have already been developed. To improve the value of these tools in decisionmaking, we suggest scientific needs that agencies such as U.S. Geological Survey can help meet—for instance, development and support of data archives. Such archives could greatly reduce resource needs and improve the reliability and consistency of results. Given the rapid state of evolution in the field, periodic follow-up studies on ecosystem services tools would help to ensure that the Bureau of Land Management and other public land management agencies are kept up to date on new tools and features that bring ecosystem services closer to readiness for use in regular decisionmaking.

Arizona;Sonora↗

USGS2001

Hydraulic and mechanical properties affecting ground-water flow and aquifer-system compaction, San Joaquin Valley, California

This report summarizes hydraulic and mechanical properties affecting ground-water flow and aquifer-system compaction in the San Joaquin Valley, a broad alluviated intermontane structural trough that constitutes the southern two-thirds of the Central Valley of California. These values will be used to constrain a coupled ground-water flow and aquifer-system compaction model of the western San Joaquin Valley called WESTSIM. A main objective of the WESTSIM model is to evaluate potential future land subsidence that might occur under conditions in which deliveries of imported surface water for agricultural use are reduced and ground-water pumping is increased. Storage values generally are components of the total aquifer-system storage and include inelastic and elastic skeletal storage values of the aquifers and the aquitards that primarily govern the potential amount of land subsidence. Vertical hydraulic conductivity values generally are for discrete thicknesses of sediments, usually aquitards, that primarily govern the rate of land subsidence. The data were compiled from published sources and include results of aquifer tests, stress-strain analyses of borehole extensometer observations, laboratory consolidation tests, and calibrated models of aquifer-system compaction.

Open-File Report↗