Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Computers and Geology”

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 631 records · Page 35Linked to original sources

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↗

Methods for estimating selected low-flow frequency and mean annual flow statistics at gaged and ungaged locations on streams in Georgia, North Carolina, and South Carolina

The U.S. Geological Survey, in cooperation with the Georgia Department of Natural Resources (Environmental Protection Division), North Carolina Department of Environmental Quality (Division of Water Resources), North Carolina Department of Public Safety (Office of Recovery and Resiliency), and South Carolina Department of Environmental Services, updated low-flow frequency, mean annual flow, and flow-duration statistics at 843 streamgages in and near Georgia, North Carolina, and South Carolina. The low-flow frequency statistics are annual minimum 1-day average flow with a 10-year recurrence interval (1Q10), annual minimum 7-day average flow for 2- and 10-year recurrence intervals (7Q2 and 7Q10, respectively), and annual minimum 30-day average flow with 2- and 3-year recurrence intervals (30Q2 and 30Q3, respectively). Monthly 1Q10 and 7Q10, and W7Q10 flow statistics for the winter period (November–March) also are presented. By using data from 604 of the streamgages on streams with streamflows that are not substantially affected by regulation or diversion and are not tidally influenced, regional regression equations were developed to predict flow statistics with prediction intervals at ungaged locations on streams with those same criteria. The regional regression analysis included data from 132 streamgages from adjacent States Alabama, Florida, Tennessee, and Virginia. The final regional regression equations include variables such as drainage area, streamflow variability, precipitation, percentage of impervious area, and percentage of the basin in various ecoregions. The low-flow statistics for the streamgages analyzed and the regional regression equations will be integrated into the U.S. Geological Survey StreamStats application ( https://www.usgs.gov/streamstats ) for Georgia, North Carolina, and South Carolina. StreamStats generates basin characteristics needed to compute low-flow frequency statistics for ungaged locations. A trend analysis of annual minimum 7-day average flows was done for 78 streamgages with at least 30 years of continuous record. Trends were evaluated for 30-, 50‑, 70-, and 90-year periods, ending in climate year 2021, and independence and short- and long-term persistence assumptions were considered. For all trend analysis assumptions, most streamgages did not exhibit significant trends in annual minimum 7-day average flows. Trends in annual precipitation and air temperature were similarly evaluated for the period 1895–2021 to assess the variability of climate for Georgia, North Carolina, and South Carolina.

Georgia, North Carolina, South Carolina↗

User's manual for the National Water-Quality Assessment Program Invertebrate Data Analysis System (IDAS) software, version 5

The Invertebrate Data Analysis System (IDAS) software was developed to provide an accurate, consistent, and efficient mechanism for analyzing invertebrate data collected as part of the U.S. Geological Survey National Water-Quality Assessment (NAWQA) Program. The IDAS software is a stand-alone program for personal computers that run Microsoft Windows(Registered). It allows users to read data downloaded from the NAWQA Program Biological Transactional Database (Bio-TDB) or to import data from other sources either as Microsoft Excel(Registered) or Microsoft Access(Registered) files. The program consists of five modules: Edit Data, Data Preparation, Calculate Community Metrics, Calculate Diversities and Similarities, and Data Export. The Edit Data module allows the user to subset data on the basis of taxonomy or sample type, extract a random subsample of data, combine or delete data, summarize distributions, resolve ambiguous taxa (see glossary) and conditional/provisional taxa, import non-NAWQA data, and maintain and create files of invertebrate attributes that are used in the calculation of invertebrate metrics. The Data Preparation module allows the user to select the type(s) of sample(s) to process, calculate densities, delete taxa on the basis of laboratory processing notes, delete pupae or terrestrial adults, combine lifestages or keep them separate, select a lowest taxonomic level for analysis, delete rare taxa on the basis of the number of sites where a taxon occurs and (or) the abundance of a taxon in a sample, and resolve taxonomic ambiguities by one of four methods. The Calculate Community Metrics module allows the user to calculate 184 community metrics, including metrics based on organism tolerances, functional feeding groups, and behavior. The Calculate Diversities and Similarities module allows the user to calculate nine diversity and eight similarity indices. The Data Export module allows the user to export data to other software packages (CANOCO, Primer, PC-ORD, MVSP) and produce tables of community data that can be imported into spreadsheet, database, graphics, statistics, and word-processing programs. The IDAS program facilitates the documentation of analyses by keeping a log of the data that are processed, the files that are generated, and the program settings used to process the data. Though the IDAS program was developed to process NAWQA Program invertebrate data downloaded from Bio-TDB, the Edit Data module includes tools that can be used to convert non-NAWQA data into Bio-TDB format. Consequently, the data manipulation, analysis, and export procedures provided by the IDAS program can be used to process data generated outside of the NAWQA Program.

Techniques and Methods↗

Incorporating uncertainty in susceptibility criteria into probabilistic liquefaction hazard analysis

Most conventional approaches for assessing liquefaction triggering hazards generally rely on simplified procedures that involve identifying liquefaction susceptible layers and calculating a factor of safety against liquefaction (FSL) in each layer. Such procedures utilize deterministic semi-empirical models for standard penetration test (SPT), cone penetrometer test (CPT), or shear wave velocity (Vs)-based subsurface data. This general approach largely neglects considerable uncertainties in ground shaking, as well as aleatory variabilities and epistemic uncertainties inherent to liquefaction susceptibility and triggering prediction. A more robust methodology introduced by Kramer and Mayfield (2007) is known as probabilistic liquefaction hazard analysis (PLHA), which integrates the full ground motion hazard space with probabilistic forms of liquefaction triggering models (e.g., Boulanger and Idriss 2014), resulting in the computation of FSL profiles with consistent return periods. Multiple PLHA computational platforms have been developed over the years, with the computational framework from Makdisi (2021) serving as the basis for a new Liquefaction Hazard Tool under development at the U.S. Geologic Survey (USGS). Despite significant improvements in recent years to the availability of seismic hazard data and probabilistic triggering and effects models, the issue of incorporating uncertainty in characterizing liquefaction susceptibility remains a challenge. Most compositional susceptibility criteria (i.e., whether or not the soil exhibits sand-like behavior) currently in use are presented as deterministic bounds based on in-situ or laboratory test data; similarly, determination of soil saturation is often based on a single groundwater level from in-situ testing. As a result, the same types of binary decisions must be made in PLHA as in more conventional methods. With the expansion and availability of field and laboratory data pertaining to liquefaction through resources such as the Next Generation Liquefaction (NGL) project, there exists the potential for an improved set of susceptibility models for CPT, SPT, and Vs-based applications. Presented here is a brief discussion on how probabilistic susceptibility modeling can be accommodated in PLHA calculations, as well as how the use of multiple models can be leveraged within a logic tree to improve the representation of epistemic uncertainty in liquefaction hazard analysis.

Conference Paper↗

Dissolved solids and sodium in water from the High Plains aquifer in parts of Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming

In 1978, the U.S. Geological Survey initiated a 5-year study of the High Plains regional aquifer to provide: (1) Hydrologic information needed to evaluate the effects of continued ground-water development; and (2) computer models to predict aquifer response to changes in ground-water development. The plan of study for the High Plains Regional Aquifer-System Analysis was described by Weeks (1978). A description of the High Plains aquifer and a map of the 1978 water table were presented by Gutentag and Weeks (1980). Maps of the bedrock geology, altitude of aquifer base, and saturated thickness of the High Plains aquifer were published by Weeks and Gutentag (1981). Water-level and saturated-thickness changes, from predevelopment to 1980, were mapped by Luckey, Gutentag, and Weeks (1981). This report describes the areal distribution of dissolved solids and sodium in the water of the High Plains aquifer. Data used in this study were provided by the U.S. Geological Survey and State agencies in each of the eight States in the High Plains. Their contribution is an integral part of this investigation.

Hydrologic Atlas↗

Tectonic and metallogenic model for northeast Asia

This document describes the digital files in this report that contains a tectonic and metallogenic model for Northeast Asia. The report also contains background materials. This tectonic and metallogenic model and other materials on this report are derived from (1) an extensive USGS Professional Paper, 1765, on the metallogenesis and tectonics of Northeast Asia that is available on the Internet at http://pubs.usgs.gov/pp/1765/; and (2) the Russian Far East parts of an extensive USGS Professional Paper, 1697, on the metallogenesis and tectonics of the Russian Far East, Alaska, and the Canadian Cordillera that is available on the Internet at http://pubs.usgs.gov/pp/pp1697/. The major purpose of the tectonic and metallogenic model is to provide, in movie format, a colorful summary of the complex geology, tectonics, and metallogenesis of the region. To accomplish this goal four steps were taken: (1) 13 time-stage diagrams, from the late Neoproterozoic (850 Ma) through the present (0 Ma), were adapted, generalized, and transformed into color static time-stage diagrams; (2) the 13 time-stage diagrams were placed in a computer morphing program to produce the model; (3) the model was examined and each diagram was successively adapted to preceding and subsequent diagrams to match the size and surface expression of major geologic units; and (4) the final version of the model was produced in successive iterations of steps 2 and 3. The tectonic and metallogenic model and associated materials in this report are derived from a project on the major mineral deposits, metallogenesis, and tectonics of the Northeast Asia and from a preceding project on the metallogenesis and tectonics of the Russian Far East, Alaska, and the Canadian Cordillera. Both projects provide critical information on bedrock geology and geophysics, tectonics, major metalliferous mineral resources, metallogenic patterns, and crustal origin and evolution of mineralizing systems for this region. The major scientific goals and benefits of the projects are to: (1) provide a comprehensive international data base on the mineral resources of the region that is the first extensive knowledge available in English; (2) provide major new interpretations of the origin and crustal evolution of mineralizing systems and their host rocks, thereby enabling enhanced, broad-scale tectonic reconstructions and interpretations; and (3) promote trade and scientific and technical exchanges between North America and eastern Asia.

Summit, Colorado↗

A computer program (MODFLOWP) for estimating parameters of a transient, three-dimensional ground-water flow model using nonlinear regression

This report documents a new version of the U.S. Geological Survey modular, three-dimensional, finite-difference, ground-water flow model (MODFLOW) which, with the new Parameter-Estimation Package that also is documented in this report, can be used to estimate parameters by nonlinear regression. The new version of MODFLOW is called MODFLOWP (pronounced MOD-FLOW*P), and functions nearly identically to MODFLOW when the ParameterEstimation Package is not used. Parameters are estimated by minimizing a weighted least-squares objective function by the modified Gauss-Newton method or by a conjugate-direction method. Parameters used to calculate the following MODFLOW model inputs can be estimated: Transmissivity and storage coefficient of confined layers; hydraulic conductivity and specific yield of unconfined layers; vertical leakance; vertical anisotropy (used to calculate vertical leakance); horizontal anisotropy; hydraulic conductance of the River, Streamflow-Routing, General-Head Boundary, and Drain Packages; areal recharge rates; maximum evapotranspiration; pumpage rates; and the hydraulic head at constant-head boundaries. Any spatial variation in parameters can be defined by the user. Data used to estimate parameters can include existing independent estimates of parameter values, observed hydraulic heads or temporal changes in hydraulic heads, and observed gains and losses along head-dependent boundaries (such as streams). Model output includes statistics for analyzing the parameter estimates and the model; these statistics can be used to quantify the reliability of the resulting model, to suggest changes in model construction, and to compare results of models constructed in different ways.

Open-File Report↗

Low-flow characteristics of streams in the Puget Sound region, Washington

Periods of low streamflow are usually the most critical factor in relation to most water uses. The purpose of this report is to present data on low-flow characteristics of streams in the Puget Sound region, Washington, and to briefly explain some of the factors that influence low flow in the various basins. Presented are data on low-flow frequencies of streams in the Puget Sound region, as gathered at 150 gaging stations. Four indexes were computed from the flow-flow-frequency curves and were used as a basis to compare the low-flow characteristics of the streams. The indexes are the (1) low-flow-yield index, expressed in unit runoff per square mile; (2) base-flow index, or the ratio of the median 7-day low flow to the average discharge; (3) slope index, or slope of annual 7-day low-flow-frequency curve; and (4) spacing index, or spread between the 7-day and 183-day low-flow-frequency curves. The indexes showed a wide variation between streams due to the complex interrelation between climate, topography, and geology. The largest low-flow-yield indexes determined--greater than 1.5 cfs (cubic feet per second) per square mile--were for streams that head at high altitudes in the Cascade and Olympic Mountains and have their sources at glaciers. The smallest low-flow-yield indexes--less than 0.5 cfs per square mile--were for the small streams that drain the lowlands adjacent to Puget Sound. Indexes between the two extremes were for nonglacial streams that head at fairly high altitudes in areas of abundant precipitation. The base-flow index has variations that can be attributed to a basin's hydrogeology, with very little influence from climate. The largest base-flow indexes were obtained for streams draining permeable unconsolidated glacial and alluvial sediments in parts of the lowlands adjacent to Puget Sound. Large volume of ground water in these materials sustain flows during late summer. The smallest indexes were computed for streams draining areas underlain by relatively impermeable igneous, sedimentary, and metamorphic rocks or by relatively impermeable glacial till. Melt water from snow and ice influences the index for streams which originate at glaciers, and result in fairly large indexes--0.25 or greater. The slope index is influenced principally by the character of the geologic materials that underlie the basin. The largest slope indexes were computed for small streams that drain areas underlain by compact glacial till or consolidated sedimentary rocks. In contrast, lowland streams that flow through areas underlain by unconsolidated alluvia and glacial deposits have the smallest indexes. Small slope indexes also are characteristic of glacial streams and show the moderating effect of the snow and ice storage in the high mountain basins. The spacing indexes are similar to the slope indexes in that they are affected by the character of the geologic materials underlying a basin. The largest spacing indexes are characteristic of small streams whose basins are underlain by glacial till or by consolidated sedimentary rocks. The smallest indexes were computed for some lowland streams draining areas underlain by permeable glacial and alluvial sediments. The indexes do not appear to have a definite relation to each other. The low-flow-yield indexes are not related to either the slope or spacing indexes because snow and ice storage has a great influence on the low-flow-yield index, while the character of the geologic materials influences the slope and spacing indexes. A relation exists between the slope and spacing indexes but many anomalies occur that cannot be explained by the geology of the basins.

Open-File Report↗

Nutrient loads in the Lost River and Klamath River Basins, south-central Oregon and northern California, March 2012–March 2015

The U.S. Geological Survey and Bureau of Reclamation collected water-quality data from March 2012 to March 2015 at locations in the Lost River and Klamath River Basins, Oregon, in an effort to characterize water quality and compute a nutrient budget for the Bureau of Reclamation Klamath Reclamation Project. The study described in this report resulted in the following significant findings: Total phosphorus (TP), total nitrogen (TN), 5-day biochemical oxygen demand (BOD5), and 5-day carbonaceous biochemical oxygen demand (CBOD5) loads, calculated using the U.S. Geological Survey LOADEST software package at the upper and lower boundaries of the Klamath Reclamation Project, indicated higher loads at the upper boundary on the southern end of Upper Klamath Lake upstream of the Bureau of Reclamation A Canal diversion compared to the lower boundary on the Klamath River downstream of Keno Dam. Accounting for the diversion of loads down A Canal, BOD5 and CBOD5 loads decreased between these two sites during irrigation season, indicating that the Klamath Reclamation Project is not a large source of oxygen-demanding material and that much of the oxygen demand at study site FMT, the northern boundary of the study area, has been expressed by the time the same water passes through site KRK, the southern boundary of the study area. An evaluation of the nutrient balance along the Klamath River flowpath from sites FMT to KRK indicated that, during irrigation season in the 3 years of the study period (March 2012–March 2015), more loads of TP, TN, BOD5, and CBOD5 were being diverted from the Klamath River than were being added to the Klamath River from the combination of Klamath Straits Drain, regulated point sources along the Klamath River, and internal loading from the bottom sediments in the river. By contrast, during non-irrigation seasons, more loads were added to the Klamath River than were diverted through Ady and North Canals, and this difference primarily was due to additional loads to the river from the Lost River Diversion Channel. At the Lost River Diversion Channel, BOD5 loads were higher during irrigation season than non-irrigation season in all three study years owing to the high concentrations of oxygen-demanding cyanobacterial biomass from the seasonal blooms of Aphanizomenon flos-aquae in the Klamath River and Upper Klamath Lake. The difference between the two seasons was particularly large in years 2 and 3, when the low flows of these two drought years resulted in smaller nonirrigation period loads than in year 1. CBOD5 loads also were higher during irrigation season in years 2 and 3 than during non-irrigation season, indicating that the largest oxygen demand was coming from senescence of Aphanizomenon flos-aquae cells that are present in the Klamath River during the summer. However, during irrigation season in year 1, CBOD5 loads were lower than in the non-irrigation season, which may indicate that at times high concentrations of ammonia or cellular organic nitrogen leaving Upper Klamath Lake contribute a large nitrogenous oxygen demand as well. The smallest loads were computed for the farthest upstream sites in the Lost River Basin, suggesting that the upper Lost River Basin does not contribute substantial loads of TP, TN, BOD5, and CBOD5 to the Klamath Reclamation Project. Median concentrations of BOD5 and CBOD5 were lowest among the upper Lost River Basin sites and highest at site PPD (however, this comparison is based on only four samples collected at site PPD over the 3-year study). Median concentrations of BOD5 and CBOD5 also were elevated at sites KSDH (6.60 and 4.70 milligrams per liter [mg/L], respectively) and KSD97 (4.47 and 3.45 mg/L, respectively). The highest maximum BOD5 and CBOD5 concentrations were reported at the Lost River Diversion Channel (39.0 and 26.5 mg/L, respectively) when water was flowing from the Klamath River toward the Klamath Reclamation Project, and site FMT (25.0 and 23.9 mg/L, respectively), the study site at the southern end of Upper Klamath Lake. Carbonaceous oxygen demand, as represented by CBOD5, typically dominated the composition of the samples at all sites. The highest concentrations of dissolved organic carbon were present at sites KSDH (the headworks of Klamath Straits Drain) and KSD97 (Klamath Straits drain before it enters the Klamath River), and PPD (outlet of Tule Lake). Median concentrations of TN and TP at the upper Lost River Basin sites in years 1 and 2 were variable, but site MCRV showed a smaller range of values in those years compared to the other upper Lost River Basins sites, and an overall lower median concentration during irrigation seasons in years 1 and 2, suggesting that Gerber Reservoir does not contribute high concentrations of nutrients to the Lost River during irrigation season. Total Maximum Daily Load (TMDL) load allocations for TP and TN in Klamath Straits Drain were exceeded in all three study years. BOD5 load allocations were exceeded in years 1 and 2, but not year 3. TMDL load allocations for TP were exceeded in the Lost River Diversion Channel for all 3 years. Load allocations for TN were exceeded in year 1, but not in years 2 and 3. BOD5 loads were less than the TMDL load allocation for all three study years. The dearth of samples collected at the Klamath Straits Drain just downstream of the Lower Klamath National Wildlife Refuge did not allow for direct assessment of the Klamath Straits Drain acting as a nutrient source or sink. TP, TN, BOD5, and CBOD5 loads estimated during the study period likely were smaller than long-term average conditions because of persistent drought conditions in the Upper Klamath Basin. The study results, therefore, fail to characterize loads from the Klamath Reclamation Project to the Klamath River that could be present in typical years, and suggest the need for load assessments during average or aboveaverage streamflow years.

California, Oregon↗

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↗

Geochemical orientation for mineral exploration in the Hashemite Kingdom of Jordan

This report is a supplement to previous accounts of geochemical exploration conducted in the Hashemite Kingdom of Jordan by the Natural Resources Authority of the Royal Government of Jordan and the U.S. Geological Survey. The field work on which this report is based was sponsored by the U.S. Agency for International Development, U.S. Department of State. Procedures used in collecting various kinds of rocks, ores, slags, eluvial and alluvial sediments, heavy-mineral concentrates, and organic materials for use as geochemical sample media are summarized, as are the laboratory procedures followed for the analysis of these sample materials by semiquantitative spectrographic, atomic absorption, fluorometric, and X-ray diffraction methods. Geochemical evaluations of the possibilities for economic mineral deposits in certain areas are presented. The results of these preliminary investigations open concepts for further use in geochemical exploration in the search for metallic mineral deposits in Jordan. Perhaps the most desirable new activity would be hydrogeochemical exploration for uranium and base metals, accompanied by interpretation of such remote-sensing data as results of airborne radiometric surveys and computer-enhanced LANDSAT imagery. For more conventional approaches to geochemical exploration, however, several fundamental problems regarding proper choice of geochemical sample media for different geologic and geographic parts of the Country must be solved before effective surveys can be made. The present results also show that such common geochemical exploration techniques as the determination of the trace-element contents of soils, plant ash, and slags have direct application also toward the resolution of several archaeological problems in Jordan. These include the relation of trace-elements chemistry of local soils to the composition of botanic remains, the trace-elements composition of slags to the technological development of the extractive metallurgy of copper and iron in the region, and the use of charcoal from slags for the C-14 dating of periods of archaeometallurgical activity. Less directly, interpretations based on the distribution in time and space of the archaeometallurgical activities of the region might add to the knowledge of early climatic conditions and vegetative cover of the area.

Open-File Report↗

Audio-magnetotelluric methods in reconnaissance geothermal exploration

An audio-magnetotelluric (AMT) system has been developed by the U.S. Geological Survey for low-cost reconnaissance exploration of geothermal regions. This is an electromagnetic sounding technique in which the scalar or Cagniard resistivity is computed at 12 frequencies logarithmically spaced from 7.5 to 18 600 Hz. Our system uses natural source fields except at the two upper frequencies of 10 200 and 18 600 Hz where artificial VLF sources are available. As a reconnaissance technique we use AMT surveys in conjunction with regional gravity, magnetic, and telluric surveys. The exploration depth is a function of the resistivities of the lithologic section, but typically ranges from the surface to 0.2 km in low-resistivity areas and to greater than 2 km in high-resistivity regions. Results of the initial reconnaissance AMT surveys provide a rational basis for deciding on the extent of costlier follow-up surveys. As part of the U.S. Geological Survey geothermal program, surveys were conducted in Long Valley and Surprise Valley, California; the Vale, Ore-Weiser, Idaho region; and Bruneau-Grand View, Raft River, and Island Park regions of Idaho. AMT surveys in five additional known geothermal resource areas (KGRA's) have been scheduled for completion by May 1975. In the Raft River and Bruneau-Grand View regions and Long Valley, follow-up electrical surveys substantiated the effectiveness of the AMT technique for reconnaissance surveying.

Conference Paper↗

Determination of dilution factors for discharge of aluminum-containing wastes by public water-supply treatment facilities into lakes and reservoirs in Massachusetts

Dilution of aluminum discharged to reservoirs in filter-backwash effluents at water-treatment facilities in Massachusetts was investigated by a field study and computer simulation. Determination of dilution is needed so that permits for discharge ensure compliance with water-quality standards for aquatic life. The U.S. Environmental Protection Agency chronic standard for aluminum, 87 micrograms per liter (μg/L), rather than the acute standard, 750 μg/L, was used in this investigation because the time scales of chronic exposure (days) more nearly match rates of change in reservoir concentrations than do the time scales of acute exposure (hours). Whereas dilution factors are routinely computed for effluents discharged to streams solely on the basis of flow of the effluent and flow of the receiving stream, dilution determination for effluents discharged to reservoirs is more complex because (1), compared to streams, additional water is available for dilution in reservoirs during low flows as a result of reservoir flushing and storage during higher flows, and (2) aluminum removal in reservoirs occurs by aluminum sedimentation during the residence time of water in the reservoir. Possible resuspension of settled aluminum was not considered in this investigation. An additional concern for setting discharge standards is the substantial concentration of aluminum that can be naturally present in ambient surface waters, usually in association with dissolved organic carbon (DOC), which can bind aluminum and keep it in solution. A method for dilution determination was developed using a mass-balance equation for aluminum and considering sources of aluminum from groundwater, surface water, and filter-backwash effluents and losses caused by sedimentation, water withdrawal, and spill discharge from the reservoir. The method was applied to 13 reservoirs. Data on aluminum and DOC concentrations in reservoirs and influent water were collected during the fall of 2009. Complete reservoir volume was determined to be available for mixing on the basis of vertical and horizontal aluminum-concentration profiling. Losses caused by settling of aluminum were assumed to be proportional to aluminum concentration and reservoir area. The constant of proportionality, as a function of DOC concentration, was established by simulations in each of five reservoirs that differed in DOC concentration. In addition to computing dilution factors, the project determined dilution factors that would be protective with the same statistical basis (frequency of exceedance of the chronic standard) as dilutions computed for streams at the 7-day-average 10-year-recurrence annual low flow (the 7Q10). Low-flow dilutions are used for permitting so that receiving waters are protected even at the worst-case flow levels. The low-flow dilution factors that give the same statistical protection are the lowest annual 7-day-average dilution factors with a recurrence of 10 years, termed 7DF10s. Determination of 7DF10 values for reservoirs required that long periods of record be simulated so that dilution statistics could be determined. Dilution statistics were simulated for 13 reservoirs from 1960 to 2004 using U.S. Geological Survey Firm-Yield Estimator software to model reservoir inputs and outputs and present-day values of filter-effluent discharge and aluminum concentration. Computed settling velocities ranged from 0 centimeters per day (cm/d) at DOC concentrations of 15.5 milligrams per liter (mg/L) to 21.5 cm/d at DOC concentrations of 2.7 mg/L. The 7DF10 values were a function of aluminum effluent discharged. At current (2009) effluent discharge rates, the 7DF10 values varied from 1.8 to 115 among the 13 reservoirs. In most cases, the present-day (2009) discharge resulted in receiving water concentrations that did not exceed the standard at the 7DF10. Exceptions were one reservoir with a very small area and three reservoirs with high concentrations of DOC. Maximum permissible discharges were determined for water-treatment plants by adjusting discharges upward in simulations until the 7DF10 resulted in reservoir concentrations that just met the standard. In terms of aluminum flux, these discharges ranged from 0 to 28 kilograms of aluminum per day.

Massachusetts↗

Description of the National Hydrologic Model for use with the Precipitation-Runoff Modeling System (PRMS)

This report documents several components of the U.S. Geological Survey National Hydrologic Model of the conterminous United States for use with the Precipitation-Runoff Modeling System (PRMS). It provides descriptions of the (1) National Hydrologic Model, (2) Geospatial Fabric for National Hydrologic Modeling, (3) PRMS hydrologic simulation code, (4) parameters and estimation methods used to compute spatially and temporally distributed default values as required by PRMS, (5) National Hydrologic Model Parameter Database, and (6) model extraction tool named Bandit. The National Hydrologic Model Parameter Database contains values for all PRMS parameters used in the National Hydrologic Model. The methods and national datasets used to estimate all the PRMS parameters are described. Some parameter values are derived from characteristics of topography, land cover, soils, geology, and hydrography using traditional Geographic Information System methods. Other parameters are set to long-established default values and computation of initial values. Additionally, methods (statistical, sensitivity, calibration, and algebraic) were developed to compute parameter values on the basis of a variety of nationally-consistent datasets. Values in the National Hydrologic Model Parameter Database can periodically be updated on the basis of new parameter estimation methods and as additional national datasets become available. A companion ScienceBase resource provides a set of static parameter values as well as images of spatially-distributed parameters associated with PRMS states and fluxes for each Hydrologic Response Unit across the conterminuous United States.

Techniques and Methods↗

Results of pumping tests on artesian wells in the Milwaukee - Waukesha area, Wisconsin

As a result of a bill passed by the Wisconsin State Legislature in 1945, ground-water investigations in Wisconsin have been under way since February 1946 under the terms of an agreement between the U. S. Geological Survey and the University of Wisconsin. Pumping tests on wells that yield water from the Ordovician and Cambrian sandstones underlying the Milwaukee-Waukesha area have been made as a part a larger regional investigation. The purpose of the tests has been to determine the water-bearing characteristics of the aquifer. These characteristics, the coefficients of transmissibility and storage, are used to determine the effect on water levels caused by changes in the rate of withdrawal from the aquifer. Average coefficients of transmissibility and storage determined from the results of 47 pumping tests at five different locations are 23,800 gallons a day per foot, and 0.00039, respectively. The amount of drawdown in the water level at any point caused by pumping a well for a given length of time may be computed by the nonequilibrium formula, using the coefficients and correcting for the effects of boundaries and of any changes in the character of the aquifer. Further study of the geology is needed to determine the location of the recharge area, the location of possible boundaries, and changes in the character of the aquifer. Collection of water-level and pumpage data is continuing and will serve as a check of the computations using the coefficients of transmissibility and storage.

Wisconsin↗

Computer Program for Point Location And Calculation of ERror (PLACER)

A program designed for point location and calculation of error (PLACER) was developed as part of the Quality Assurance Program of the Federal Highway Administration/U.S. Geological Survey (USGS) National Data and Methodology Synthesis (NDAMS) review process. The program provides a standard method to derive study-site locations from site maps in highwayrunoff, urban-runoff, and other research reports. This report provides a guide for using PLACER, documents methods used to estimate study-site locations, documents the NDAMS Study-Site Locator Form, and documents the FORTRAN code used to implement the method. PLACER is a simple program that calculates the latitude and longitude coordinates of one or more study sites plotted on a published map and estimates the uncertainty of these calculated coordinates. PLACER calculates the latitude and longitude of each study site by interpolating between the coordinates of known features and the locations of study sites using any consistent, linear, user-defined coordinate system. This program will read data entered from the computer keyboard and(or) from a formatted text file, and will write the results to the computer screen and to a text file. PLACER is readily transferable to different computers and operating systems with few (if any) modifications because it is written in standard FORTRAN. PLACER can be used to calculate study site locations in latitude and longitude, using known map coordinates or features that are identifiable in geographic information data bases such as USGS Geographic Names Information System, which is available on the World Wide Web.

Open-File Report↗

An evaluation of selected extraordinary floods in the United States reported by the U.S. Geological Survey and implications for future advancement of flood science

Thirty flood peak discharges determine the envelope curve of maximum floods documented in the United States by the U.S. Geological Survey. These floods occurred from 1927 to 1978 and are extraordinary not just in their magnitude, but in their hydraulic and geomorphic characteristics. The reliability of the computed discharge of these extraordinary floods was reviewed and evaluated using current (2007) best practices. Of the 30 flood peak discharges investigated, only 7 were measured at daily streamflow-gaging stations that existed when the flood occurred, and 23 were measured at miscellaneous (ungaged) sites. Methods used to measure these 30 extraordinary flood peak discharges consisted of 21 slope-area measurements, 2 direct current-meter measurements, 1 culvert measurement, 1 rating-curve extension, and 1 interpolation and rating-curve extension. The remaining four peak discharges were measured using combinations of culvert, slope-area, flow-over-road, and contracted-opening measurements. The method of peak discharge determination for one flood is unknown. Changes to peak discharge or rating are recommended for 20 of the 30 flood peak discharges that were evaluated. Nine floods retained published peak discharges, but their ratings were downgraded. For two floods, both peak discharge and rating were corrected and revised. Peak discharges for five floods that are subject to significant uncertainty due to complex field and hydraulic conditions, were re-rated as estimates. This study resulted in 5 of the 30 peak discharges having revised values greater than about 10 percent different from the original published values. Peak discharges were smaller for three floods (North Fork Hubbard Creek, Texas; El Rancho Arroyo, New Mexico; South Fork Wailua River, Hawaii), and two peak discharges were revised upward (Lahontan Reservoir tributary, Nevada; Bronco Creek, Arizona). Two peak discharges were indeterminate because they were concluded to have been debris flows with peak discharges that were estimated by an inappropriate method (slope-area) (Big Creek near Waynesville, North Carolina; Day Creek near Etiwanda, California). Original field notes and records could not be found for three of the floods, however, some data (copies of original materials, records of reviews) were available for two of these floods. A rating was assigned to each of seven peak discharges that had no rating. Errors identified in the reviews include misidentified flow processes, incorrect drainage areas for very small basins, incorrect latitude and longitude, improper field methods, arithmetic mistakes in hand calculations, omission of measured high flows when developing rating curves, and typographical errors. Common problems include use of two-section slope-area measurements, poor site selection, uncertainties in Manning’s n -values, inadequate review, lost data files, and insufficient and inadequately described high-water marks. These floods also highlight the extreme difficulty in making indirect discharge measurements following extraordinary floods. Significantly, none of the indirect measurements are rated better than fair, which indicates the need to improve methodology to estimate peak discharge. Highly unsteady flow and resulting transient hydraulic phenomena, two-dimensional flow patterns, debris flows at streamflow-gaging stations, and the possibility of disconnected flow surfaces are examples of unresolved problems not well handled by current indirect discharge methodology. On the basis of a comprehensive review of 50,000 annual peak discharges and miscellaneous floods in California, problems with individual flood peak discharges would be expected to require a revision of discharge or rating curves at a rate no greater than about 0.10 percent of all floods. Many extraordinary floods create complex flow patterns and processes that cannot be adequately documented with quasi-steady, uniform one-dimensional analyses. These floods are most accurately described by multidimensional flow analysis. Within the U.S. Geological Survey, new approaches are needed to collect more accurate data for floods, particularly extraordinary floods. In recent years, significant progress has been made in instrumentation for making direct discharge measurements. During this same period, very little has been accomplished in advancing methods to improve indirect discharge measurements. Greater use of paleoflood hydrology could fill many shortcomings of U.S. Geological Survey flood science today, such as enhanced knowledge of flood frequency. Additional links among flood runoff, storm structure, and storm motion would provide more insight to flood hazards. Significant improvement in understanding flood processes and characteristics could be gained from linking radar rainfall estimation and hydrologic modeling. Additionally, more could be done to provide real-time flood-hazard warnings with linked rainfall/runoff and flow models. Several important recommendations are made to improve the flood-documentation capability of the U.S. Geological Survey. When very large discharges are measured by current meter or hydroacoustics, water-surface slope should be measured as well. This measurement would allow validation of roughness values that can significantly extend the discharge range of verified Manning’s n for 1-dimensional and 2-dimensional flow analyses. At least two of the floods investigated may have had flow so unstable that large waves affected the interpretation of high-water marks. Instability criteria should be considered for hydraulic analysis of large flows in high-gradient, smooth channels. The U.S. Geological Survey needs to modernize its toolbox of field and office practices for making future indirect discharge measurements. These practices could include, first and foremost, a new peak-flow file database that allows greater description and interpretation of flow events, such as stability criteria in high-gradient, smooth channels, debris flow documentation, and details of flood genesis (hurricane, snowmelt, rain-on-snow, dam failure, and the like). Other modernized practices could include (a) establishment of calibrated stream reaches in chronic flash flood basins to expedite indirect computation of flow; (b) development of process-based theoretical rating curves for streamflow-gaging stations; (c) adoption of step-backwater models as the standard surface-water modeling tool for U.S. Geological Survey field offices; (d) development and support for multidimensional flow models capable of describing flood characteristics in complex terrain and high-gradient channels; (e) greater use of the critical-depth method in appropriate locations; (f) deployment of non-contact instruments to directly measure large floods, rather than attempting to reconstruct them; (g) increased use of paleoflood hydrology; and (h) assurance that future collection of hydro-climatic data meets the needs of more robust watershed models.

Scientific Investigations Report↗

Regional Aquifer-System Analysis— Appalachian Valley and Piedmont

The Regional Aquifer-System Analysis Program, RASA, represents a systematic effort to study a number of the Nation’s most important aquifer systems, which, in aggregate, underlie much of the country and which represent an important component of the Nation’s total water supply. In general, the boundaries of these studies are identified by the hydrologic extent of each system and, accordingly, transcend the political subdivisions to which investigations have often arbitrarily been limited in the past. The broad objective for each study is to assemble geologic, hydrologic, and geochemical information, to analyze and develop an understanding of the system, and to develop predictive capabilities that will contribute to the effective management of the system. The use of computer simulation is an important element of the RASA studies to develop an understanding of the natural, undisturbed hydrologic system and the changes brought about in it by human activities and to provide a means of predicting the regional effects of future pumping or other stresses. The final interpretive results of the RASA Program are presented in a series of U.S. Geological Survey Professional Papers that describe the geology, hydrology, and geochemistry of each regional aquifer system. Each study within the RASA Program is assigned a single Professional Paper number beginning with Professional Paper 1400. This paper, Professional Paper 1422, represents the Regional Aquifer-System Analysis— Appalachian Valley and Piedmont. It is published as several individual volumes over several years.

Professional Paper↗