Search USGS⌕ Search

SEARCH · Search USGS

Results for “Final Report”

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

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

At least 1,243 records · Page 69Linked to original sources
USGS2020

Elevation-derived hydrography acquisition specifications

Hydrographic features derived from U.S. Geological Survey (USGS) 3D Elevation Program data, and collected for use by the USGS, must meet the specifications described in this document. The specifications described herein pertain to the final product delivered to the USGS, not to methods used to derive the hydrographic features. The specifications describe the collection area, spatial reference system, attribute table structure, feature codes and values, delineation of hydrographic features, topology, positional assessment, metadata, and delivery formats. A companion document, Elevation-Derived Hydrography—Representation, Extraction, Attribution, and Delineation Rules, defines the fields, domains, and minimum feature collection requirements for hydrography features derived from elevation data. Hydrographic features collected to this specification will be suitable for using as breaklines to hydroflatten digital elevation models, processing for preconflation of features to the National Hydrography Dataset, and using for hydroenforcement of digital elevation models.

Techniques and Methods↗
USGS1985

Hydrogeology, water quality, and ground-water development alternatives in the Beaver-Pasquiset ground-water reservoir, Rhode Island

In a 23 sq mi study area, the Beaver-Pasquiset groundwater reservoir within the Pawcatuck River basin in southern Rhode Island, stratified drift is the only principal geologic unit capable of producing yields > 350 gal/min. Transmissivity of the aquifer ranges from 7,200 to 24,300 sq ft/day. Water table conditions prevail in the aquifer, which is in good hydraulic connection with perennial streams and ponds. A digital model of two-dimensional groundwater flow was used to simulate the interaction between surface water and groundwater, and to evaluate the impact of alternative schemes of groundwater development on groundwater levels, pond levels, and streamflow in the Beaver-Pasquiset groundwater reservoir. Transient simulations of theoretical pumpage were made for a drought period (1963-66) and a wet period (1976-78). The areas most favorable for development of high-capacity wells (350 gal/min or more) are along the Beaver River and near Pasquiset Pond. The water is soft and generally contains < 100 mg/L dissolved solids. Locally, groundwater contains elevated concentrations of iron and manganese (7.5 and 3.7 mg/L, respectively), southeast of Pasquiset Pond, and will require treatment if used for public supply. The groundwater reservoir was simulated with a two-dimensional finite-difference model using a block-centered grid consisting of 33 rows and 75 columns. Differences between measured and simulated water table altitudes for the final steady state run for 21 selected observation wells averaged +0.07 ft. Combined pumping rates for simulation of groundwater development alternatives at eight sites ranged from 3.25 to 7.00 Mgal/d. Pumping rates for individual wells ranged from 0.25 to 1.50 Mgal/d. Transient simulations suggest that the Beaver-Pasquiset groundwater reservoir is capable of sustaining a pumping rate of 4.25 Mgal/d during years of average groundwater recharge with minimal impact on groundwater levels, pond levels, and streamflow. During extreme drought periods (1965 and 1966) it would be necessary to reduce pumpage below 3.25 Mgal/d to maintain flow in both the Beaver River and Pasquiset Brook. (Author 's abstract)

Water-Resources Investigations Report↗
USGS2000

Sources and haloacetic acid/trihalomethane formation potentials of aquatic humic substances in the Wakarusa River and Clinton Lake near Lawrence, Kansas

Gram quantities of aquatic humic substances (AHS) were extracted from the Wakarusa River−Clinton Lake Reservoir system, near Lawrence, KS, to support nuclear magnetic resonance (NMR) experimental studies, report concentrations of dissolved organic carbon (DOC) and AHS, define sources of the AHS, and determine if the AHS yield sufficient quantities of haloacetic acids (HAA5) and trihalomethanes (THM4) that exceed U.S. Environmental Protection Agency (EPA) Maximum Contaminant Levels (MCL) in drinking water. AHS from the Wakarusa River and Clinton Lake originated from riparian forest vegetation, reflected respective effects of soil organic matter and aquatic algal/bacterial sources, and bore evidence of biological degradation and photodegradation. AHS from the Wakarusa River showed the effect of terrestrial sources, whereas Clinton Lake humic acid also reflected aquatic algal/bacterial sources. Greater amounts of carbon attributable to tannin-derived chemical structures may correspond with higher HAA5 and THM4 yields for Clinton Lake fulvic acid. Prior to appreciable leaf-fall from deciduous trees, the combined (humic and fulvic acid) THM4 formation potentials for the Wakarusa River approached the proposed EPA THM4 Stage I MCL of 80 μg/L, and the combined THM4 formation potential for Clinton Lake slightly exceeded the proposed THM4 Stage II MCL of 40 μg/L. Finally, AHS from Clinton Lake could account for most (>70%) of the THM4 concentrations in finished water from the Clinton Lake Water Treatment Plant based on September 23, 1996, THM4 results.

Environmental Science & Technology↗
USGS2002

Cyclic injection, storage, and withdrawal of heated water in a sandstone aquifer at St. Paul, Minnesota--Analysis of thermal data and nonisothermal modeling of short-term test cycles

In May 1980, the University of Minnesota began a project to evaluate the feasibility of storing heated water (150 degrees Celsius) in the Franconia-Ironton Galesville aquifer (183 to 245 meters below land surface) and later recovering it for space heating. The University's steam-generation facilities supplied high-temperature water for injection. The Aquifer Thermal-Energy Storage system is a doublet-well design in which the injection-withdrawal wells are spaced approximately 250 meters apart. Water was pumped from one of the wells through a heat exchanger, where heat was added or removed. This water was then injected back into the aquifer through the other well. Four short-term test cycles were completed. Each cycle consisted of approximately equal durations of injection and withdrawal ranging from 5.25 to 8.01 days. Equal rates of injection and withdrawal, ranging from 17.4 to 18.6 liters per second, were maintained for each short-term test cycle. Average injection temperatures ranged from 88.5 to 117.9 degrees Celsius. Temperature graphs for selected depths at individual observation wells indicate that the Ironton and Galesville Sandstones received and stored more thermal energy than the upper part of the Franconia Formation. Clogging of the Ironton Sandstone was possibly due to precipitation of calcium carbonate or movement of fine-grain material or both. Vertical-profile plots indicate that the effects of buoyancy flow were small within the aquifer. A three-dimensional, anisotropic, nonisothermal, ground-water-flow, and thermal-energy-transport model was constructed to simulate the four short-term test cycles. The model was used to simulate the entire short-term testing period of approximately 400 days. The only model properties varied during model calibration were longitudinal and transverse thermal dispersivities, which, for final calibration, were simulated as 3.3 and 0.33 meters, respectively. The model was calibrated by comparing model-computed results to (1) measured temperatures at selected altitudes in four observation wells, (2) measured temperatures at the production well, and (3) calculated thermal efficiencies of the aquifer. Model-computed withdrawal-water temperatures were within an average of about 3 percent of measured values and model-computed aquifer-thermal efficiencies were within an average of about 5 percent of calculated values for the short-term test cycles. These data indicate that the model accurately simulated thermal-energy storage within the Franconia-Ironton-Galesville aquifer.

Professional Paper↗
USGS1999

External quality-assurance results for the National Atmospheric Deposition Program/National Trends Network, 1995-96

The U.S. Geological Survey operated four external quality-assurance programs for the National Atmospheric Deposition Program/National Trends Network (NADP/NTN) in 1995 and 1996: the intersite-comparison program, the blind-audit program, the interlaboratory- comparison program, and the collocated-sampler program. The intersite-comparison program assessed the precision and bias of pH and specific-conductance determinations made by NADP/NTN site operators. The analytical bias introduced during routine handling, processing, and shipping of wet-deposition samples and precision of analyte values was estimated using a blind-audit program. An interlaboratory-comparison program was used to evaluate differences between analytical results and to estimate the analytical precision of five North American laboratories that routinely analyzed wet deposition. A collocated-sampler program estimated the precision of the overall precipitation collection and analysis system from initial sample collection through final storage of the data. Results of two intersite-comparison studies completed in 1995 indicated 94.6 and 94.4 percent of the onsite pH determinations met the NADP/NTN accuracy goals, whereas 97.2 and 98.3 percent of the specific-conductance determinations were within the established limits. The percentages of onsite determinations that met the accuracy goals in 1996 were slightly less for both pH and specific-conductance than in 1995. In 1996, 93.2 and 87.5 percent of onsite pH determinations met the accuracy goals, whereas the percentage of onsite specific-conductance measurements that met the goals was 93.9 and 94.9 percent.The blind audit program utilizes a paired sample design to evaluate the effects of routine sample handling, processing and shipping on the chemistry of weekly precipitation samples. The portion of the blind audit sample subject to all of the normal onsite handling and processing steps of a regular weekly precipitation sample is referred to as the bucket portion, whereas the portion receiving only minimal handling is referred to as the bottle portion. Throughout the report, the term positive bias in regard to blind-audit results indicates that the bucket portion had a higher concentration than the bottle portion. The paired t-test of 1995 blind-audit data indicated that routine sample handling, processing, and shipping introduced a very small positive bias (a=0.05) for hydrogen ion and specific conductance and a slight negative bias (a =0.05) for ammonium and sodium. In 1995, the median paired differences between the bucket and bottle portions ranged from -0.02 milligram per liter for both ammonium and nitrate to +0.002 milligram per liter for calcium. Although the paired t-test indicated a very small positive bias for hydrogen ion, the median paired difference between the bucket and bottle portions was 0.00 microequivalents per liter, whereas for specific conductance, the median paired difference between the bucket and bottle portions was 0.200 microsiemens per centimeter in 1995. The paired t-test of blind-audit results in 1996 indicated statistically significant bias for 6 of the 10 analytes. Only chloride, nitrate, hydrogen ion, and specific conductance were not biased in 1996. However, the magnitude of the bias in 1996 was very small and only of limited importance from the viewpoint of an analytical chemist or data user. The median paired differences between the bucket and bottle portions ranged from -0.02 milligram per liter for both ammonium and chloride to +0.006 milligram per liter for calcium. For hydrogen ion, the median paired difference between the bucket and bottle portions was -0.357 microequivalent per liter; for specific conductance, the median paired difference between the bucket and bottle portions was 0.00 microsiemens per centimeter in 1996. Surface-chemistry effects due to different amounts of precipitation contacting the sample collection and shipping container surfac

Water-Resources Investigations Report↗
USGS2000

Methods for estimating low-flow statistics for Massachusetts streams

Methods and computer software are described in this report for determining flow duration, low-flow frequency statistics, and August median flows. These low-flow statistics can be estimated for unregulated streams in Massachusetts using different methods depending on whether the location of interest is at a streamgaging station, a low-flow partial-record station, or an ungaged site where no data are available. Low-flow statistics for streamgaging stations can be estimated using standard U.S. Geological Survey methods described in the report. The MOVE.1 mathematical method and a graphical correlation method can be used to estimate low-flow statistics for low-flow partial-record stations. The MOVE.1 method is recommended when the relation between measured flows at a partial-record station and daily mean flows at a nearby, hydrologically similar streamgaging station is linear, and the graphical method is recommended when the relation is curved. Equations are presented for computing the variance and equivalent years of record for estimates of low-flow statistics for low-flow partial-record stations when either a single or multiple index stations are used to determine the estimates. The drainage-area ratio method or regression equations can be used to estimate low-flow statistics for ungaged sites where no data are available. The drainage-area ratio method is generally as accurate as or more accurate than regression estimates when the drainage-area ratio for an ungaged site is between 0.3 and 1.5 times the drainage area of the index data-collection site. Regression equations were developed to estimate the natural, long-term 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, and 50-percent duration flows; the 7-day, 2-year and the 7-day, 10-year low flows; and the August median flow for ungaged sites in Massachusetts. Streamflow statistics and basin characteristics for 87 to 133 streamgaging stations and low-flow partial-record stations were used to develop the equations. The streamgaging stations had from 2 to 81 years of record, with a mean record length of 37 years. The low-flow partial-record stations had from 8 to 36 streamflow measurements, with a median of 14 measurements. All basin characteristics were determined from digital map data. The basin characteristics that were statistically significant in most of the final regression equations were drainage area, the area of stratified-drift deposits per unit of stream length plus 0.1, mean basin slope, and an indicator variable that was 0 in the eastern region and 1 in the western region of Massachusetts. The equations were developed by use of weighted-least-squares regression analyses, with weights assigned proportional to the years of record and inversely proportional to the variances of the streamflow statistics for the stations. Standard errors of prediction ranged from 70.7 to 17.5 percent for the equations to predict the 7-day, 10-year low flow and 50-percent duration flow, respectively. The equations are not applicable for use in the Southeast Coastal region of the State, or where basin characteristics for the selected ungaged site are outside the ranges of those for the stations used in the regression analyses. A World Wide Web application was developed that provides streamflow statistics for data collection stations from a data base and for ungaged sites by measuring the necessary basin characteristics for the site and solving the regression equations. Output provided by the Web application for ungaged sites includes a map of the drainage-basin boundary determined for the site, the measured basin characteristics, the estimated streamflow statistics, and 90-percent prediction intervals for the estimates. An equation is provided for combining regression and correlation estimates to obtain improved estimates of the streamflow statistics for low-flow partial-record stations. An equation is also provided for combining regression and drainage-area ratio estimates to obtain improved e

Massachusetts↗
USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗
USGS2017

Conversing with Pelehonuamea: A workshop combining 1,000+ years of traditional Hawaiian knowledge with 200 years of scientific thought on Kīlauea volcanism

The events surrounding volcanic eruptions and damaging earthquakes in Hawai‘i have often been described in journals, letters, and newspapers articles in the English language; however, the Hawaiian nation was among the most literate of countries in the 19th century, and many Hawaiian-language newspapers were in circulation through all but the earliest decades of the 19th century. Any modern reconstruction of the history of Hawaiian eruptions or earthquakes should take advantage of all available sources, and so we seek to add the Hawaiian-language newspaper articles, journals, stories, and chants to the volcano and earthquake literature. These sources have been used in many recent volcanological studies. Another aspect to the mix of science and traditional Hawaiian values is that many of the volcanic summits in Hawaiʻi are considered sacred to Hawaiians. Hawaiian travelers brought the first Western missionary team to the summit of Kīlauea and advised them of the proper protocols and behaviors while in this sacred area. The missionaries dismissed this advice as native superstition and they began a campaign of aggressively stamping out customs and protocols related to the Hawaiian volcano goddess Pelehonuamea. What has survived as native knowledge of the volcanoes is a few phrases from native guides included in some of the missionaries’ journals, and a few stories. Pualani and Ku ʻ ulei Kanahele provide excellent introductions to the Pelehonuamea chants. In the 21st century, amid a reawakening of Hawaiian culture, modern Hawaiians are demanding protection of their sacred areas, and scientists must be aware of these interests. At the very least, scientists should show respect to Hawaiian values when working in these areas, and should try to minimize disruptions caused by their work. Kaeo Duarte, Peter Mills, and Scott Rowland describe taking this approach in their field work. Traditional knowledge is also contained in place names. It is important not only to preserve old place names and to recover those no longer used, but also to preserve the stories of those places. Bobby Camara talks about the joys and frustrations of getting information on and recovering Hawaiian place names. Finally, we hope that a broader interest in Hawaiian views about locations in Hawaiʻi where physical scientific work is done will be as beneficial to physical scientists as it has been to life scientists investigating Hawaiian lifeforms on land and in the ocean, and that both studies will continue to benefit the native peoples of Hawaiʻi. Note that these proceedings are transcripts of oral presentations illustrated with PowerPoint presentations or charts. Although every effort has been made to assure the accuracy of the oral presentations, there are some gaps where words are not discernible in the voice recordings and are so noted. In other places, bracketed words were added to clarify the speaker’s meaning.

Hawaii↗
USGS2022

Areas contributing recharge to priority wells in valley-fill aquifers in the Neversink River and Rondout Creek drainage basins, New York

In southeastern New York, the villages of Ellenville, Wurtsboro, Woodridge, the hamlet of Mountain Dale, and surrounding communities in the Neversink River and Rondout Creek drainage basins rely on wells that pump groundwater from valley-fill glacial aquifers for public water supply. Glacial aquifers are vulnerable to contamination because they are highly permeable and have a shallow depth to water table. To protect the quality of these water resources, water managers need accurate information about the areas that contribute recharge to production wells that pump from these aquifers. The New York State Department of Environmental Conservation and the New York State Department of Health designated eight priority wells in this region for which water supply protection is of primary concern. The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation and the New York State Department of Health, began an investigation in 2019 with the general objectives of (1) improving understanding of regional groundwater-flow system, (2) delineating areas contributing recharge to eight priority production wells, and (3) quantifying the uncertainty of these contributing areas in a probabilistic way that can be used to inform decision-making related to priority well source-water protection. To complete these objectives, a MODFLOW 6 groundwater model was created encompassing the eight priority wells and the surrounding flow system, which includes parts of the Neversink River and Rondout Creek Basins in Sullivan County and Ulster County, New York. The model was built using Python tools (such as flopy , modflow-setup , and sfrmaker ) that facilitate transparent and repeatable model development using existing datasets. The model parameters were estimated with a stepwise approach using an iterative ensemble smoother implementation of the Parameter ESTimation software PEST++ (version 5.0.0). We evaluated initial “best guess” parameter bounds with a prior Monte Carlo analysis. Results of the first prior Monte Carlo analysis were used to make informed adjustments to model parameter bounds (typically resulting in expanded bounds), and a second prior Monte Carlo analysis was run to identify improved ranges for model parameters during history matching. The history matching effort produced an ensemble of parameter values for the groundwater-flow model that spans the range of values within prior uncertainty bounds. The ensemble is informed by the historical observation data, within a reasonable range of uncertainty on those observations. This history-matched ensemble was used in a particle tracking Monte Carlo analysis to delineate the areas contributing recharge to priority wells. The groundwater-flow and particle tracking (MODPATH7) models were run once for each ensemble member. Deterministic contributing areas computed for each ensemble member were aggregated to produce maps showing the probability that a location contributes recharge to priority wells. Finally, the particle tracking Monte Carlo analysis was repeated for six pumping scenarios, representing a wide range of possible pumping levels, to incorporate uncertainty in future pumping rates related to population growth or other management decisions. Increasing pumping rates generally led to larger contributing recharge areas and larger areas of high probability that a location contributes recharge to priority wells. These maps show the overall uncertainty of the areas contributing recharge to priority wells in the study area and provide a tool for risk-based decision making for protection of well source water.

New York↗
USGS2023

Uses of epistemic uncertainties in the USGS National Seismic Hazard Models

The need for US Geological Survey (USGS) National Seismic Hazard Models (NSHMs) to report estimates of epistemic uncertainties in the hazard (e.g. fractile hazard curves) in all forthcoming releases is increasing. With fractile hazard curves as potential new outputs from the USGS 2023 NSHM, a simultaneous need is to help end-users better understand these epistemic uncertainties and clarify their potential uses. In this article, we address the latter need by (1) characterizing epistemic uncertainties in two updates of the USGS NSHM (2014 for California and 2021 for Hawaii), (2) illustrating a variety of downstream applications of fractile hazard curves in both hazard and risk contexts, and (3) discussing implications from the various types of uncertainties. We found that the epistemic uncertainty in hazard is generally larger for Hawaii than for California, the epistemic uncertainty in hazard can be reasonably approximated with a lognormal distribution for most of the cases considered, and the correlation between epistemic uncertainty in hazard at two different intensity measure levels generally varies with both location and type of intensity measure. Furthermore, we developed models for readily generating approximate fractile hazard curves in California and Hawaii. Finally, given the complexities involved in the hazard modeling process, we developed an open-source interactive tool to enable a broad range of users to independently examine and potentially start using such epistemic uncertainties for their respective applications.

Earthquake Spectra↗
USGS2014

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Farah mineral district in Afghanistan

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Defense Task Force for Business and Stability Operations, prepared databases for mineral-resource target areas in Afghanistan. The purpose of the databases is to (1) provide useful data to ground-survey crews for use in performing detailed assessments of the areas and (2) provide useful information to private investors who are considering investment in a particular area for development of its natural resources. The set of satellite-image mosaics provided in this Data Series (DS) is one such database. Although airborne digital color-infrared imagery was acquired for parts of Afghanistan in 2006, the image data have radiometric variations that preclude their use in creating a consistent image mosaic for geologic analysis. Consequently, image mosaics were created using ALOS (Advanced Land Observation Satellite; renamed Daichi) satellite images, whose radiometry has been well determined (Saunier, 2007a,b). This part of the DS consists of the locally enhanced ALOS image mosaics for the Farah mineral district, which has spectral reflectance anomalies indicative of copper, zinc, lead, silver, and gold deposits. ALOS was launched on January 24, 2006, and provides multispectral images from the AVNIR (Advanced Visible and Near-Infrared Radiometer) sensor in blue (420-500 nanometer, nm), green (520-600 nm), red (610-690 nm), and near-infrared (760-890 nm) wavelength bands with an 8-bit dynamic range and a 10-meter (m) ground resolution. The satellite also provides a panchromatic band image from the PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) sensor (520-770 nm) with the same dynamic range but a 2.5-m ground resolution. The image products in this DS incorporate copyrighted data provided by the Japan Aerospace Exploration Agency ((c)JAXA, 2007, 2008, 2010), but the image processing has altered the original pixel structure and all image values of the JAXA ALOS data, such that original image values cannot be recreated from this DS. As such, the DS products match JAXA criteria for value added products, which are not copyrighted, according to the ALOS end-user license agreement. The selection criteria for the satellite imagery used in our mosaics were images having (1) the highest solar-elevation angles (near summer solstice) and (2) the least cloud, cloud-shadow, and snow cover. The multispectral and panchromatic data were orthorectified with ALOS satellite ephemeris data, a process which is not as accurate as orthorectification using digital elevation models (DEMs); however, the ALOS processing center did not have a precise DEM. As a result, the multispectral and panchromatic image pairs were generally not well registered to the surface and not coregistered well enough to perform resolution enhancement on the multispectral data. Therefore, it was necessary to (1) register the 10-m AVNIR multispectral imagery to a well-controlled Landsat image base, (2) mosaic the individual multispectral images into a single image of the entire area of interest, (3) register each panchromatic image to the registered multispectral image base, and (4) mosaic the individual panchromatic images into a single image of the entire area of interest. The two image-registration steps were facilitated using an automated control-point algorithm developed by the USGS that allows image coregistration to within one picture element. Before rectification, the multispectral and panchromatic images were converted to radiance values and then to relative-reflectance values using the methods described in Davis (2006). Mosaicking the multispectral or panchromatic images started with the image with the highest sun-elevation angle and the least atmospheric scattering, which was treated as the standard image. The band-reflectance values of all other multispectral or panchromatic images within the area were sequentially adjusted to that of the standard image by determining band-reflectance correspondence between overlapping images using linear least-squares analysis. The resolution of the multispectral image mosaic was then increased to that of the panchromatic image mosaic using the SPARKLE logic, which is described in Davis (2006). Each of the four-band images within the resolution-enhanced image mosaic was individually subjected to a local-area histogram stretch algorithm (described in Davis, 2007), which stretches each band's picture element based on the digital values of all picture elements within a 500-m radius. The final databases, which are provided in this DS, are three-band, color-composite images of the local-area-enhanced, natural-color data (the blue, green, and red wavelength bands) and color-infrared data (the green, red, and near-infrared wavelength bands). All image data were initially projected and maintained in Universal Transverse Mercator (UTM) map projection using the target area's local zone (41 for Farah) and the WGS84 datum. The final image mosaics were subdivided into four overlapping tiles or quadrants because of the large size of the target area. The four image tiles (or quadrants) for the Farah area are provided as embedded geotiff images, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image. Within the Farah study area, five subareas were designated for detailed field investigations (that is, the FarahA through FarahE subareas); these subareas were extracted from the area’s image mosaic and are provided as separate embedded geotiff images.

Farah Mineral District↗
USGS1982

Laboratory investigations of the physics of steam flow in a porous medium

Experiments were carried out in the laboratory to test a theory of transient flow of pure steam in a uniform porous medium. This theory is used extensively in modeling pressure-transient behavior in vapor-dominated geothermal systems. Transient, superheated steam-flow experiments were run by bringing a cylinder of porous material to a uniform initial pressure, and then making a step increase in pressure at one end of the sample, while monitoring the pressure-transient breakthrough at the other end. It was found in experiments run at 100?, 125?, and 146?C that the time required for steam-pressure transients to propagate through an unconsolidated material containing sand, silt, and clay was 10 to 25 times longer than predicted by theory. It is hypothesized that the delay in the steam-pressure transient was caused by adsorption of steam in the porous sample. In order to account for steam adsorption, a sink term was included in the conservation of mass equation. In addition, energy transfer in the system has to be considered because latent heat is released when steam adsorption occurs, increasing the sample temperature by as much as 10?C. Finally, it was recognized that the steam pressure was a function of both the temperature and the amount of adsorption in the sample. For simplicity, this function was assumed to be in equilibrium adsorption isotherm, which was determined by experiment. By solving the modified mass and energy equations numerically, subject to the empirical adsorption isotherm relationship, excellent theoretical simulation of the experiments was achieved. The experiments support the hypothesis that adsorption of steam can strongly influence steam pressure-transient behavior in porous media; the results suggest that the modified steam-flow theory, which includes steam adsorption terms, should be used in modeling steam flow in vapor-dominated geothermal systems.

Open-File Report↗
USGS2024

Ecology shapes the genomic and biosynthetic diversification of Streptomyces bacteria from insectivorous bats

Streptomyces are prolific producers of secondary metabolites from which many clinically useful compounds have been derived. They inhabit diverse habitats but have rarely been reported in vertebrates. Here, we aim to determine to what extent the ecological source (bat host species and cave sites) influence the genomic and biosynthetic diversity of Streptomyces bacteria. We analysed draft genomes of 132 Streptomyces isolates sampled from 11 species of insectivorous bats from six cave sites in Arizona and New Mexico, USA. We delineated 55 species based on the genome-wide average nucleotide identity and core genome phylogenetic tree. Streptomyces isolates that colonize the same bat species or inhabit the same site exhibit greater overall genomic similarity than they do with Streptomyces from other bat species or sites. However, when considering biosynthetic gene clusters (BGCs) alone, BGC distribution is not structured by the ecological or geographical source of the Streptomyces that carry them. Each genome carried between 19–65 BGCs (median=42.5) and varied even among members of the same Streptomyces species. Nine major classes of BGCs were detected in ten of the 11 bat species and in all sites: terpene, non-ribosomal peptide synthetase, polyketide synthase, siderophore, RiPP-like, butyrolactone, lanthipeptide, ectoine, melanin. Finally, Streptomyces genomes carry multiple hybrid BGCs consisting of signature domains from two to seven distinct BGC classes. Taken together, our results bring critical insights to understanding Streptomyces -bat ecology and BGC diversity that may contribute to bat health and in augmenting current efforts in natural product discovery, especially from underexplored or overlooked environments.

Microbial Genomics↗
USGS2004

Generalized estimates from streamflow data of annual and seasonal ground-water-recharge rates for drainage basins in New Hampshire

This report presents regression equations to estimate generalized annual and seasonal ground-water-recharge rates in drainage basins in New Hampshire. The ultimate source of water for a ground-water withdrawal is aquifer recharge from a combination of precipitation on the aquifer, ground-water flow from upland basin areas, and infiltration from streambeds to the aquifer. An assessment of ground-water availability in a basin requires that recharge rates be estimated under `normal' conditions and under assumed drought conditions. Recharge equations were developed by analyzing streamflow, basin characteristics, and precipitation at 55 unregulated continuous record stream-gaging stations in New Hampshire and in adjacent states. In the initial step, streamflow records were analyzed to estimate a series of annual and seasonal ground-water-recharge components of streamflow in each drainage basin evaluated in this study. Regression equations were then developed relating the series of annual and seasonal ground-water-recharge values to the corresponding series of annual and seasonal precipitation values as determined at the centroid of each drainage basin. This resulted in one equation for each of the 55 basins for each of the four seasonal periods and the annual period, or a total of 275 regression equations. Average annual and seasonal precipitation data for 1961-90 were then used to compute a set of normalized ground-water-recharge values that reflected the long-term average annual and seasonal variations (normalized) and mean recharge characteristics of each drainage basin. Ordinary-least-squares regression was applied in the process of selecting 10 out of 93 possible basin and climatic characteristics for further testing in the development of the equations for computing the generalized estimate of annual and seasonal ground-water recharge based on the set of normalized recharge values. Generalized-least-squares regression was used for the final parameter estimation and error evaluation. The following basin and climatic characteristics were found to be statistically significant predictors for at least one of the dependent variables: average annual, summer, and spring precipitation as determined at U.S. Geological Survey stream-gaging stations; average annual basin-centroid precipitation; average mean annual basin temperature; average minimum winter basin temperature; percent coniferous forest in a basin; percent mixed coniferous and deciduous forest in a basin; average fall basin-centroid precipitation; and average annual snowcover. These 10 basin and climatic characteristics were selected because they were statistically significant based on several statistical parameters that evaluated which combination of characteristics contributed the most to the predictive accuracy of the regression-equation models. A geographic information system is required to measure the values of the predictor variables for the equations developed in the study. The average annual normalized ground-water recharge was 21.0 in. This value was determined by generalized-least-squares (GLS) regression analysis for all of the basins used in the normalized ground-water recharge analysis for precipitation from 1961-90. The average winter (January 1-March 15) ground-water recharge was 4.3 in., average spring (March 16-May 31) ground-water recharge was 9.0 in., average summer (June 1-October 31) ground-water recharge was 4.0 in., and average fall (November 1-December 31) ground-water recharge was 3.6 in. Normalized ground-water recharge ranged annually from 12.3 to 31.8 in., for winter from 2.30 to 7.82 in., for spring from 5.16 to 13.7 in., for summer from 1.45 to 10.2 in., and for fall from 2.21 to 6.06 in.

New Hampshire↗
USGS2005

Factors influencing mercury in freshwater surface sediments of northeastern North America

We report on an inventory and analysis of sediment mercury (Hg) concentrations from 579 sites across northeastern North America. Sediment Hg concentrations ranged from the limit of detection ca. 0.01-3.7 ??g g -1 (dry weight, d.w.), and the average concentration was 0.19 ??g g-1 (d.w.) Sediment methylmercury concentrations ranged from 0.15 to 21 ng g-1 (d.w.) and the mean concentration was 3.83 ng g -1 (d.w.). Total Hg concentrations (HgT) were greatest in lakes > reservoirs > rivers, although the proportion of Hg as methylmercury showed an inverse pattern. Total Hg was weakly and positively correlated with the sediment organic matter and percent of watershed as forested land, and weakly and negatively correlated with sediment solids content, drainage area, and agricultural land. Sediment methylmercury concentrations were weakly and positively correlated to wetland area, and weakly and negatively correlated to drainage area. Methylmercury, expressed as a percentage of HgT was positively correlated to agricultural land area. For sites with co-located sediment and fish-tissue sampling results, there was no relationship between sediment Hg and fish-tissue Hg. Finally, our data indicate that at least 44% of waters across the region have sediment HgT concentrations in excess of Canadian and United States minimum sediment contaminant guidelines for the protection of aquatic biota. ?? 2005 Springer Science+Business Media, Inc.

Ecotoxicology↗
USGS2024

Direct measurements of sediment geoacoustic properties in the New England Mud Patch and shelf break

This paper reports on an original set of direct sound speed measurements collected with the acoustic coring system in the New England Mud Patch (NEMP) and shelf break area to the south. Cores collected within the NEMP show range-dependence of the mud with slower sound speed and lower attenuation on the west side. In the shelf break region, the highest sound speeds are observed between the 200- and 350-m isobaths. The depth-dependence of the mud layer in the NEMP includes a surficial layer with a negative sound speed gradient of 28 s −1 . The remainder of the mud column has a weak positive sound speed gradient of 6.2 s −1 over an isovelocity layer. Comparison between in situ and ex situ sound speed measurements provides an assessment of the effects of sediment disturbance from gravity coring operations. Small differences in the upper 2.5 m were attributed to the changes in the geoacoustic properties caused by disturbance from the coring process. Below 2.5 m, the average difference is close to zero, suggesting that these sediments were minimally disturbed. Finally, an in situ measurement of shear speed was obtained near the depth of maximum penetration. The shear speed was well correlated with sound speed from approximately the same depth interval.

New England Mud Patch↗
USGS2007

Riparian Land Use/Land Cover Data for Five Study Units in the Nutrient Enrichment Effects Topical Study of the National Water-Quality Assessment Program

This dataset was developed as part of the National Water-Quality Assessment (NAWQA) Program, Nutrient Enrichment Effects Topical (NEET) study for five study units distributed across the United States: Apalachicola-Chattahoochee-Flint River Basin, Central Columbia Plateau-Yakima River Basin, Central Nebraska Basins, Potomac River Basin and Delmarva Peninsula, and White, Great and Little Miami River Basins. One hundred forty-three stream reaches were examined as part of the NEET study conducted 2003-04. Stream segments, with lengths equal to the logarithm of the basin area, were delineated upstream from the downstream ends of the stream reaches with the use of digital orthophoto quarter quadrangles (DOQQ) or selected from the high-resolution National Hydrography Dataset (NHD). Use of the NHD was necessary when the stream was not distinguishable in the DOQQ because of dense tree canopy. The analysis area for each stream segment was defined by a buffer beginning at the segment extending to 250 meters lateral to the stream segment. Delineation of land use/land cover (LULC) map units within stream segment buffers was conducted using on-screen digitizing of riparian LULC classes interpreted from the DOQQ. LULC units were mapped using a classification strategy consisting of nine classes. National Wetlands Inventory (NWI) data were used to aid in wetland classification. Longitudinal transect sampling lines offset from the stream segments were generated and partitioned into the underlying LULC types. These longitudinal samples yielded the relative linear extent and sequence of each LULC type within the riparian zone at the segment scale. The resulting areal and linear LULC data filled in the spatial-scale gap between the 30-meter resolution of the National Land Cover Dataset and the reach-level habitat assessment data collected onsite routinely for NAWQA ecological sampling. The final data consisted of 12 geospatial datasets: LULC within 25 meters of the stream reach (polygon); LULC within 50 meters of the stream reach (polygon); LULC within 50 meters of the stream segment (polygon); LULC within 100 meters of the stream segment (polygon); LULC within 150 meters of the stream segment (polygon); LULC within 250 meters of the stream segment (polygon); frequency of gaps in woody vegetation LULC at the reach scale (arc); stream reaches (arc); longitudinal LULC at the reach scale (arc); frequency of gaps in woody vegetation LULC at the segment scale (arc); stream segments (arc); and longitudinal LULC at the segment scale (arc).

Data Series↗
USGS2002

Historic bluff retreat and stabilization at Flag Harbor, Chesapeake Bay, Maryland

Studies of bluff erosion and slope stability along the western shore of Chesapeake Bay suggest relative evolution from steep, eroding coastal bluffs to stable slopes at angles of repose ca. 35 degrees over decades. Because of the dating methods in those studies, it was impossible to precisely define rates of change. The present study provides historic age control. A pair of small harbor structures were constructed in the early 1950's at Chesapeake Beach, MD to maintain a dredged channel to a small marina occupying a ravine in the Calvert Cliffs. Prior to construction, this section of shoreline was comprised of eroding steep bluffs cut into Miocene-age sediments. Downdrift erosion is now apparent south of the structures as is updrift deposition behind the northern jetty. Since construction the updrift sand body has prograded northward and progressively deposited protective beaches along the toes of the bluffs. Former eroding bluffs nearest the harbor are now stable, vegetated slopes at angles near 35 degrees. Slope angles widen to the north and to the northern limit of the sand body. Beyond this are eroding bluffs standing at angles of 70-80 degrees. The relative time required for eroding bluffs to reach stability is estimated by interpolating the distance and time for the sand body to prograde northward since harbor construction. We measured slope angles at intervals northward from the updrift structure for a distance of 2000 feet. A least squares regression of slope angle vs distance showed progressive decrease in angle from north to south. Actively eroding 70-80 degree bluffs gave way to vegetated, but slumping slopes, and finally to stable 35-degree slopes at the harbor. A relationship between time and distance along the shore allowed us to estimate a stabilization time for this location of 35-40 years. The shortness of this time scale allows us to suggest that attempts to artificially stabilize eroding bluffs along this coast is not a simple task of protecting the toes of slopes from wave action. Once shoreline retreat ends, sloughing of sediment from bluff faces gives way to longer-term landslide processes. The bluff top recedes until a stable 35-degree slope is attained. Thus, simple shoreline protection methods may not preserve property at the bluff edge.

Open-File Report↗