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,279 records · Page 71Linked to original sources

Delineation of the hydrogeologic framework of the Big Sioux aquifer near Sioux Falls, South Dakota, using airborne electromagnetic data

The U.S. Geological Survey, in cooperation with the City of Sioux Falls, South Dakota, began developing a groundwater-flow model of the Big Sioux aquifer in 2014 that will enable the City to make more informed water management decisions, such as delineation of areas of the greatest specific yield, which is crucial for locating municipal wells. Innovative tools are being evaluated as part of this study that can improve the delineation of the hydrogeologic framework of the aquifer for use in development of a groundwater-flow model, and the approach could have transfer value for similar hydrogeologic settings. The first step in developing a groundwater-flow model is determining the hydrogeologic framework (vertical and horizontal extents of the aquifer), which typically is determined by interpreting geologic information from drillers’ logs and surficial geology maps. However, well and borehole data only provide hydrogeologic information for a single location; conversely, nearly continuous geophysical data are collected along flight lines using airborne electromagnetic (AEM) surveys. These electromagnetic data are collected every 3 meters along a flight line (on average) and subsequently can be related to hydrogeologic properties. AEM data, coupled with and constrained by well and borehole data, can substantially improve the accuracy of aquifer hydrogeologic framework delineations and result in better groundwater-flow models. AEM data were acquired using the Resolve frequency-domain AEM system to map the Big Sioux aquifer in the region of the city of Sioux Falls. The survey acquired more than 870 line-kilometers of AEM data over a total area of about 145 square kilometers, primarily over the flood plain of the Big Sioux River between the cities of Dell Rapids and Sioux Falls. The U.S. Geological Survey inverted the survey data to generate resistivity-depth sections that were used in two-dimensional maps and in three-dimensional volumetric visualizations of the Earth resistivity distribution. Contact lines were drawn using a geographic information system to delineate interpreted geologic stratigraphy. The contact lines were converted to points and then interpolated into a raster surface. The methods used to develop elevation and depth maps of the hydrogeologic framework of the Big Sioux aquifer are described herein. The final AEM interpreted aquifer thickness ranged from 0 to 31 meters with an average thickness of 12.8 meters. The estimated total volume of the aquifer was 1,060,000,000 cubic meters based on the assumption that the top of the aquifer is the land-surface elevation. A simple calculation of the volume (length times width times height) of a previous delineation of the aquifer estimated the aquifer volume at 378,000,000 cubic meters; thus, the estimation based on AEM data is more than twice the previous estimate. The depth to top of Sioux Quartzite, which ranged in depth from 0 to 90 meters, also was delineated from the AEM data.

South Dakota↗

Pollution of Illinois and Mississippi Rivers by Chicago sewage, a digest of the testimony taken in the case of the State of Missouri v. the State of Illinois and the Sanitary District of Chicago

The testimony taken in the suit of the State of Missouri against the State of Illinois and the sanitary district of Chicago comprises the best symposium on river pollution, its biological and chemical aspects, and its general and special sanitary significance that has ever been assembled. The contentions of both parties to the suit' are supported by the most eminently qualified men in the United States. The evidence presented and the discussions recorded are therefore of unique importance. The final record of testimony occupies 8,000 printed pages, much of which is irrelevant. This digest of testimony is the result of an attempt to recover the valuable material and present it in concise form. A consistent endeavor has been made by the reviewer to eliminate all personal opinions with reference to the issue and to make an impartial presentation of so much of the testimony as in his opinion appears to be relevant and of scientific importance. It will be well to remember in this connection that any digest of so large a volume of testimony must be the result of a final exercise of personal opinion by the reviewer as to those parts which may best be excluded. Naturally opinions will differ on this point; therefore it will be strange if many of those familiar with the case do not find that certain portions of testimony which they consider most important are passed over in this digest without reference. Controversies between counsel, objections to the admission of testimony, legal technicalities and quibbles, badgering cross-examination, and in general all the testimony introduced for purposes of mere corroboration have been disregarded. The object has been to present a faithful statement of the scientific phases of the testimony to the exclusion, if need be, of the legal aspect of the case.

Illinois, Missouri↗

Forecasting stream habitat and Brook Trout responses to climate change in Catoctin Mountain Park

Anticipating and mitigating the effects of climate change is a fundamental challenge for natural resource conservation. In this report, we respond to research needs identified by Catoctin Mountain Park (CATO) for native Brook Trout (Salvelinus fontinalis) conservation and management as part of the US Geological Survey (USGS) Natural Resources Preservation Program in FY15-16. We addressed three overarching research questions: (1) How will anticipated changes in air temperature affect stream habitats? (2) How will changes to stream habitat affect the distribution of Brook Trout? (3) Which stream segments are most and least vulnerable to the effects of climate change? First, we surveyed Brook Trout abundance and fish community composition using electrofishing techniques within three watersheds: Owens Creek, upper Big Hunting Creek, and Blue Blazes Creek (a tributary to Big Hunting Creek). Second, we deployed a network of stream temperature gages to assess spatial variation in stream temperature and groundwater (GW) influence. Third, we used modeling techniques to forecast future stream temperatures that account for GW influences and air temperature scenarios. Fish sampling detected 13 species and 15,345 individual fish, the majority of which were Blacknose Dace (60%), Blue Ridge Sculpin (26%), and Brook Trout (6%). Brook Trout were not observed in Blue Blazes Creek and exhibited higher densities in Owens Creek than upper Big Hunting Creek (average densities = 19 fish/100 m and 4 fish/100 m, respectively). In contrast, Brown Trout were present in Blue Blazes Creek and exhibited greater density in Blue Blazes Creek than either Owens Creek or upper Big Hunting Creek (average densities = 3.0 fish/100 m, 0.3 fish/100 m, and 1.7 fish/100 m, respectively). Brown Trout occurred in sympatry with Brook Trout in Owens Creek and upper Big Hunting Creek, but appeared to have replaced Brook Trout in Blue Blazes Creek. Our fish surveys also revealed important locations for Brook Trout reproduction and young-of-year (YOY) dispersal within the Owens Creek watershed. Our study also revealed surprising differences in the distribution of Blue Ridge Sculpin among CATO streams. This species was abundant in Owens Creek (average density = 83 fish/100 m) but was less common in Blue Blazes Creek (average density = 12 fish/100 m) and was not detected in upper Big Hunting Creek. Histological examination of several specimens from Blue Blazes Creek by V. Blazer at the USGS Leetown Science Center revealed the presence of a novel parasite (Dermosystidium sp.) which has been linked to fish population declines elsewhere (Blazer et al. 2016). The parasite was not detected in Blue Ridge Sculpin samples from Owens Creek, and all trout appeared to be uninfected. Our survey results suggest that Blue Ridge Sculpin have been extirpated from upper Big Hunting Creek and have not recolonized from downstream source populations due to the fish passage barrier of Cunningham Falls. We recommend additional research to (1) evaluate the feasibility of reintroducing Blue Ridge Sculpin into upper Big Hunting Creek and (2) continue monitoring the distribution and potential spread of Dermocystidium in downstream waters. Stream temperatures ranged from 9.6 – 27.6 ºC during baseflow conditions in 2015 and 2016. Sites within upper Big Hunting Creek were consistently warmer than in Owens Creek or Blue Blazes Creek, suggesting an effect of headwater ponds outside CATO on upper Big Hunting Creek temperatures. For instance, in 2016 the maximum observed temperature in upper Big Hunting Creek was 27.6 ºC whereas Owens Creek reached a maximum of 23.7 ºC that year. Stream temperature data also revealed that 2016 was warmer than 2015 throughout the study area but did not exceed thermal tolerance limits for Brook Trout in either year. We estimated the influence of GW on stream temperatures using a statistical modeling approach based on the relationship between daily mean air temperature and stream temperature over time. Results indicated that effects of GW were generally stronger in the Owens Creek watershed than in Blue Blazes or upper Big Hunting Creek. However, we detected substantial spatial variation in GW influence among Owens Creek sites, with stream temperatures at some locations showing relatively little GW influence and others showing very strong influences (and correspondingly small influence of daily mean air temperatures). Although incoming lateral seeps were detected in upper Big Hunting Creek (D. Ferrier, Hood College, personal communication), the strongest effects of GW in the study area were due to GW upwelling within portions of the Owens Creek watershed (i.e., Tributary C in Figure 4) where we also observed high numbers of Brook Trout juveniles. Our results therefore identified potential high-priority areas for Brook Trout conservation in CATO. Finally, we modeled future stream temperatures based on scenarios characterizing GW sensitivity to air temperature and future air temperature increases. Stream temperature forecasts revealed important differences in habitat suitability for Brook Trout within and among watersheds. Big Hunting Creek sites were generally more sensitive to air temperature increases than sites in Owens Creek or Blue Blazes Creek. For instance, an increase in mean annual air temperature of 1.5 ºC (lowest level evaluated) exceeded thermal thresholds for Brook Trout in the majority of sites within that watershed, regardless of GW influence levels. In contrast, an air temperature increase of 1.5 ºC did not exceed thermal thresholds for Brook Trout in Owens Creek. However, modeled air temperature increases of 5 ºC resulted in a loss of Brook Trout thermal suitability throughout the study area. Model results revealed spatially patchy responses to air temperature increases that could provide an early-warning system for trout monitoring designs in CATO.

Maryland↗

Genetic diversity and mate selection in a reintroduced population of gray wolves

The genetic composition of an individual can markedly affect its survival, reproduction, and ultimately fitness. As some wildlife populations become smaller, conserving genetic diversity will be a conservation challenge. Many imperiled species are already supported through population augmentation efforts and we often do not know if or how genetic diversity is maintained in translocated species. As a case study for understanding the maintenance of genetic diversity in augmented populations, I wanted to know if genetic diversity (i.e., observed heterozygosity) remained high in a population of gray wolves in the Rocky Mountains of the U.S. > 20 years after reintroduction. Additionally, I wanted to know if a potential mechanism for such diversity was individuals with below average genetic diversity choosing mates with above average diversity. I also asked whether there was a preference for mating with unrelated individuals. Finally, I hypothesized that mated pairs with above average heterozygosity would have increased survival of young. Ultimately, I found that females with below average heterozygosity did not choose mates with above average heterozygosity and wolves chose mates randomly with respect to genetic relatedness. Pup survival was not higher for mated pairs with above average heterozygosity in my models. The dominant variables predicting pup survival were harvest rate during their first year of life and years pairs were mated. Ultimately, genetic diversity was relatively unchanged > 20 years after reintroduction. The mechanism for maintaining such diversity does not appear related to individuals preferentially choosing more genetically diverse mates. Inbreeding avoidance, however, appears to be at least one mechanism maintaining genetic diversity in this population.

Idaho↗

Continuous borehole strain and pore pressure in the near field of the 28 September 2004 M 6.0 Parkfield, California, earthquake: Implications for nucleation, fault response, earthquake prediction and tremor

Near-field observations of high-precision borehole strain and pore pressure, show no indication of coherent accelerating strain or pore pressure during the weeks to seconds before the 28 September 2004 M 6.0 Parkfield earthquake. Minor changes in strain rate did occur at a few sites during the last 24 hr before the earthquake but these changes are neither significant nor have the form expected for strain during slip coalescence initiating fault failure. Seconds before the event, strain is stable at the 10 −11 level. Final prerupture nucleation slip in the hypocentral region is constrained to have a moment less than 2 × 10 12 N m ( M 2.2) and a source size less than 30 m. Ground displacement data indicate similar constraints. Localized rupture nucleation and runaway precludes useful prediction of damaging earthquakes. Coseismic dynamic strains of about 10 microstrain peak-to-peak were superimposed on volumetric strain offsets of about 0.5 microstrain to the northwest of the epicenter and about 0.2 microstrain to the southeast of the epicenter, consistent with right lateral slip. Observed strain and Global Positioning System ( gps ) offsets can be simply fit with 20 cm of slip between 4 and 10 km on a 20-km segment of the fault north of Gold Hill ( M 0 = 7 × 10 17 N m). Variable slip inversion models using gps data and seismic data indicate similar moments. Observed postseismic strain is 60% to 300% of the coseismic strain, indicating incomplete release of accumulated strain. No measurable change in fault zone compliance preceding or following the earthquake is indicated by stable earth tidal response. No indications of strain change accompany nonvolcanic tremor events reported prior to and following the earthquake.

California↗

Reclamation of mined lands in the western coal region

In 1978, a group of scientists from several Federal agencies examined reclamation work at 22 coal mines in seven western States. The results of these examinations were not used to derive quantitative predictions of the outcome of reclamation work but rather to determine the general requirements for revegetation success. Locally, reclamation efforts are affected by climate, especially precipitation; the landform of the restored surface; the nature of the overburden material; the nature of the surface soil; and the natural ecological system. The goals of reclamation efforts are now broader than ever. Regulations call not only for reducing the steepness of the final surface and establishing a cover of mostly perennial native vegetation, but for restoring the land for specific land uses, achieving diversity both in types of plants and in number of species, and reintroduction of biological and ecological processes. If specific sites are monitored over a long enough period of time, quantitative predictions of success for individual mines may be possible, and such predictions can be included in environmental impact statements to help in the decision-making process. The results of this study indicate that current reclamation objectives can be met when the reclamation effort is designed on the basis of site-specific needs and when existing technology is used.

Circular↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Panjsher Valley mineral district in Afghanistan: Chapter M in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts 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 Panjsher Valley mineral district, which has emerald and silver-iron 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 (©JAXA, 2009, 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 315-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 (42 for Panjsher Valley) and the WGS84 datum. The final image mosaics were subdivided into two overlapping tiles or quadrants because of the large size of the target area. The two image tiles (or quadrants) for the Panjsher Valley 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 Panjsher Valley study area, two subareas were designated for detailed field investigations (that is, the Emerald and Silver-Iron subareas); these subareas were extracted from the area’s image mosaic and are provided as separate embedded geotiff images.

Panjsher Valley Mineral District↗

Air- and stream-water-temperature trends in the Chesapeake Bay region, 1960-2014

Water temperature is a basic, but important, measure of the condition of all aquatic environments, including the flowing waters in the streams that drain our landscape and the receiving waters of those streams. Climatic conditions have a strong influence on water temperature, which is therefore naturally variable both in time and across the landscape. Changes to natural water-temperature regimes, however, can result in a myriad of effects on aquatic organisms, water quality, circulation patterns, recreation, industry, and utility operations. For example, most species of fish, insects, and other organisms, as well as aquatic vegetation, are highly dependent on water temperature. Warming waters can result in shifts in floral and faunal species distributions, including invasive species and pathogens previously unable to inhabit the once cooler streams. Many chemical processes are temperature dependent, with reactions occurring faster in warmer conditions, leading to degraded water quality as contaminants are released into waterways at greater rates. Circulation patterns in receiving waters, such as bays and estuaries, can change as a result of warmer inflows from streams, thereby affecting organisms in those receiving waters. Changes in abundance of some aquatic species and (or) degradation of water quality can reduce the recreational value of water bodies as waters are perceived as less desirable for water-related activities or as sportfish become less available for anglers. Finally, increasing water temperatures can affect industry and utilities as the thermal capacity is reduced, making the water less effective for cooling purposes. Chesapeake Bay is the largest estuary in the United States. Eutrophication, the enrichment of a water body with excess nutrients, has plagued the bay for decades and has led to extensive restoration efforts throughout the bay watershed. The warming of stream water can exacerbate eutrophication through increased release of nutrients from in-stream sediments, so understanding changes in stream-water temperature throughout the bay watershed is critical to resource managers seeking to restore the bay ecosystem. The U.S. Environmental Protection Agency (EPA) uses indicators that “represent the state or trend of certain environmental or societal conditions … to track and better understand the effects of changes in the Earth’s climate” (U.S. Environmental Protection Agency, 2014). Updates to these indicators are published biennially by the EPA. The U.S. Geological Survey (USGS), in cooperation with the EPA, has completed analyses of air- and stream-water-temperature trends in the Chesapeake Bay region to be included as an indicator in a future release of the EPA report.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Detection and measurement of organic lampricide residues

The selective lampricide, 3-trifluoromethyl-4-nitrophenol (TFM), and its synergist, 5,2'-dichloro-4'-nitrosalicylanilide (DCN), are separable from natural waters by anion exchange. The adsorbed compounds can then be recovered from the resin as concentrates by elution with selective solvent mixtures. Measurements of the amounts of lampricides in the final concentrates can be made colorimetrically at 395 mI? for TFM and at 530 mI? for the safranin complex of DCN. TFM has also been separated for quantitative determination from homogenates of whole fish. The fish is first macerated in a blender and then hydrolyzed in hot, 3 N hydrochloric acid. The amount of background color, due to certain components of the fish in the hydrolysate, is reduced by one or a combination of three methods: (1) a series of three extractions with ether, methylene chloride, and benzene; (2) cation exchange followed by methylene chloride extraction; or (3) ether extraction followed by anion exchange and subsequent desorption with amyl acetate-acetic acid.

Technical Report↗

Performance and retention of lightweight satellite radio tags applied to the ears of polar bears ( Ursus maritimus )

Background Satellite telemetry studies provide information that is critical to the conservation and management of species affected by ecological change. Here we report on the performance and retention of two types (SPOT-227 and SPOT-305A) of ear-mounted Argos-linked satellite transmitters (i.e., platform transmitter terminal, or PTT) deployed on free-ranging polar bears in Eastern Greenland, Baffin Bay, Kane Basin, the southern Beaufort Sea, and the Chukchi Sea during 2007–2013. Results Transmissions from 142 out of 145 PTTs deployed on polar bears were received for an average of 69.3 days. The average functional longevity, defined as the number of days they transmitted while still attached to polar bears, for SPOT-227 was 56.8 days and for SPOT-305A was 48.6 days. Thirty-four of the 142 (24%) PTTs showed signs of being detached before they stopped transmitting, indicating that tag loss was an important aspect of tag failure. Furthermore, 10 of 26 (38%) bears that were re-observed following application of a PTT had a split ear pinna, suggesting that some transmitters were detached by force. All six PTTs that were still on bears upon recapture had lost the antenna, which indicates that antenna breakage was a significant contributor to PTT failure. Finally, only nine of the 142 (6%) PTTs—three of which were still attached to bears—had a final voltage reading close to the value indicating battery exhaustion. This suggests that battery exhaustion was not a major factor in tag performance. Conclusions The average functional longevity of approximately 2 months for ear-mounted PTTs (this study) is poor compared to PTT collars fitted to adult female polar bears, which can last for several years. Early failure of the ear-mounted PTTs appeared to be caused primarily by detachment from the ear or antenna breakage. We suggest that much smaller and lighter ear-mounted transmitters are necessary to reduce the risk of tissue irritation, tissue damage, and tag detachment, and with a more robust antenna design. Our results are applicable to other tag types (e.g., iridium and VHF systems) and to research on other large mammals that cannot wear radio collars.

Animal Biotelemetry↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Khanneshin mineral district in Afghanistan: Chapter A in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts 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 Khanneshin mineral district, which has uranium, thorium, rare-earth-element, and apatite 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 (©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 Khanneshin) and the WGS84 datum. The final image mosaics were subdivided into nine overlapping tiles or quadrants because of the large size of the target area. The nine image tiles (or quadrants) for the Khanneshin 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 Khanneshin study area, one subarea was designated for detailed field investigations (that is, the Khanneshin volcano subarea); this subarea was extracted from the area's image mosaic and is provided as separate embedded geotiff images.

Helm;Nimroz↗

Paleoflood investigations to improve peak-streamflow regional-regression equations for natural streamflow in eastern Colorado, 2015

The U.S. Geological Survey (USGS), in cooperation with the Colorado Department of Transportation, developed regional-regression equations for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, 0.2-percent annual exceedance-probability discharge (AEPD) for natural streamflow in eastern Colorado. A total of 188 streamgages, consisting of 6,536 years of record and a mean of approximately 35 years of record per streamgage, were used to develop the peak-streamflow regional-regression equations. The estimated AEPDs for each streamgage were computed using the USGS software program PeakFQ. The AEPDs were determined using systematic data through water year 2013. Based on previous studies conducted in Colorado and neighboring States and on the availability of data, 72 characteristics (57 basin and 15 climatic characteristics) were evaluated as candidate explanatory variables in the regression analysis. Paleoflood and non-exceedance bound ages were established based on reconnaissance-level methods. Multiple lines of evidence were used at each streamgage to arrive at a conclusion (age estimate) to add a higher degree of certainty to reconnaissance-level estimates. Paleoflood or nonexceedance bound evidence was documented at 41 streamgages, and 3 streamgages had previously collected paleoflood data.To determine the peak discharge of a paleoflood or non-exceedanc bound, two different hydraulic models were used. The mean standard error of prediction ( SEP ) for all 8 AEPDs was reduced approximately 25 percent compared to the previous flood-frequency study. For paleoflood data to be effective in reducing the SEP in eastern Colorado, a larger ratio than 44 of 188 (23 percent) streamgages would need paleoflood data and that paleoflood data would need to increase the record length by more than 25 years for the 1-percent AEPD. The greatest reduction in SEP for the peak-streamflow regional-regression equations was observed when additional new basin characteristics were included in the peak-streamflow regional-regression equations and when eastern Colorado was divided into two separate hydrologic regions. To make further reductions in the uncertainties of the peak-streamflow regional-regression equations in the Foothills and Plains hydrologic regions, additional streamgages or crest-stage gages are needed to collect peak-streamflow data on natural streams in eastern Colorado. Generalized-Least Squares regression was used to compute the final peak-streamflow regional-regression equations for peak-streamflow. Dividing eastern Colorado into two new individual regions at –104° longitude resulted in peak-streamflow regional-regression equations with the smallest SEP . The new hydrologic region located between –104° longitude and the Kansas-Nebraska State line will be designated the Plains hydrologic region and the hydrologic region comprising the rest of eastern Colorado located west of the –104° longitude and east of the Rocky Mountains and below 7,500 feet in the South Platte River Basin and below 9,000 feet in the Arkansas River Basin will be designated the Foothills hydrologic region.

Colorado↗

Estimating monthly and annual streamflow statistics at ungaged sites in Idaho

Updated monthly and annual streamflow information for the many ungaged streams throughout Idaho is needed to assist planners and managers with issues regarding fish and wildlife, water rights, and other land and water uses. To provide this information, the U.S. Geological Survey used a multiple-regression analysis to develop equations for estimating daily mean discharge exceeded 80, 50, and 20 percent of the time each month (80-, 50-, and 20-percent monthly exceedances) and mean annual discharge values at ungaged sites. The analysis produced estimating equations that relate specific streamflow statistics to basin characteristics. The standard errors of estimate, in log 10 , of the final estimating equations ranged from a minimum of 0.072 (+18.1 to -15.3 percent) to a maximum of 0.498 (+214.8 to -68.2 percent). The estimating equations might not yield reliable results for sites with basin characteristic values outside of the range of values used to develop the equations. The equations also are not applicable for regulated streams or those that exhibit significant gains and (or) losses as a result of spring flow, seepage through permeable streambeds, or irrigation diversions. The equations generally were more reliable for estimating the high monthly streamflow statistics (20-percent exceedances) than for estimating the low streamflow statistics (80-percent exceedances) in each region and all streamflow statistics in regions where annual streamflow is relatively low.

Idaho;Montana;Nevada;Oregon;Utah;Washington;Wyomin↗

Proceedings of the fourth conference on research in California's National Parks

The papers in this proceedings were selected from the 63 presentations given at the Fourth Biennial Conference on Research in National Parks in California. The overall theme for this meeting was a recurring one: “The Integration of Research into National Park Service Resource Management Decisions.” The conference was held at the University of California, Davis, on 10-12 September 1991 and was sponsored by the National Park Service Cooperative Park Studies Unit and the Institute of Ecology at the University of California, Davis. This proceedings highlights a variety of research and resource management efforts to improve the stewardship of our most treasured landscapes. In the future, it will become increasingly more important for federal and state agencies, university scientists and students, and the public to cooperate fully to improve the quantity and quality of science and resource management programs in units of the National Park System. As many of the papers in this proceedings attest, we must look beyond the political boundaries of protected landscapes to incorporate entire ecosystems. Competing resource uses inside and outside parklands must be reexamined to weave a common thread of biological conservation. As scientists, our studies must bridge the gap from plots to landscapes and from landscapes to regions. Our studies must built on information from species and populations to ecosystems and the processes that influence them. The papers in this proceedings are modest but important contribution to those ideals. Each paper represents original research and has been peer-reviewed. Many agencies, institutions, and individuals contributed in the development of this product. In the planning stages, National Park Service Western Region scientists provided advice and assistance in structuring the conference format. University of California, Davis, graduate students, directed by Sharon Lynch, assisted with logistics at the meeting, and provided general assistance with the paper sessions. We thank the students, technicians, and support staff who labored unselfishly behind the scenes to collect data and transfer information important to the completion of each paper. Financial support for this publication was provided with funds from the Western Region Resource Management and Resource Program and the Washington Office Servicewide Publications Program of the National Park Service. We thank Denny Fenn, Stan Albright, Bruce Kilgore, Dave Cherry, Gene Hester, and Donna O'Leary of the National Park Service and Charles Goldman of the University of California, Davis, Institute of Ecology for their support throughout the development of the proceedings. Finally, we are entering into a new age of research in the U.S. Department of the Interior. Recently, plans were announced for the creation of the National Biological Survey, dedicated to improving the quantity and quality of research necessary for the wise stewardship of our natural resources and wildlands. We welcome our new partners and pledge our full support to both the new and old agencies. The papers presented here--and past research efforts in national parks--provide an important stepping stone to the future. While many of us will be transferred administratively to the National Biological Survey, our emotional ties to the National Park Service mission of preserving resources for future generations will remain with us forever.

California↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Badakhshan mineral district in Afghanistan: Chapter F in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts 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 Badakhshan mineral district, which has 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), 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. For this particular area, PRISM image orthorectification was performed by the Alaska Satellite Facility, applying its photogrammetric software to PRISM stereo images with vertical control points obtained from the digital elevation database produced by the Shuttle Radar Topography Mission (Farr and others, 2007) and horizontal adjustments based on a controlled Landsat image base (Davis, 2006). The 10-m AVNIR multispectral imagery was then coregistered to the orthorectified PRISM images and individual multispectral and panchromatic images were mosaicked into single images of the entire area of interest. The image coregistration was 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 (42 for Badakhshan) and the WGS84 datum. The final image mosaics were subdivided into six overlapping tiles or quadrants because of the large size of the target area. The six image tiles (or quadrants) for the Badakhshan 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 Badakhshan study area, three subareas were designated for detailed field investigations (that is, the Bharak, Fayz-Abad, and Ragh subareas); these subareas were extracted from the area’s image mosaic and are provided as separate embedded geotiff images.

Badakhshan↗

Two-dimensional resistivity investigation of the North Cavalcade Street site, Houston, Texas, August 2003

The North Cavalcade Street site was first developed for wood treating in 1946. By 1955, pentachlorophenol wood preservation services and other support facilities, such as creosote ponds, pentachlorophenol and creosote storage structures, various tanks, lumber sheds, a treatment facility, and other buildings had been added. In 1961, the property was closed. To protect public health and welfare and the environment from release or threatened releases of hazardous substances, the U.S. Environmental Protection Agency added the North Cavalcade Street site to the National Priorities List on October 5, 1984. Between September 1985 and November 1987, the U.S. Environmental Protection Agency conducted a remedial investigation which, through exploratory drilling, determined the locations of two contaminated source areas and a normal fault. During August 2003, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, conducted a two-dimensional (2D) resistivity investigation at the North Cavalcade Street site to provide additional characterization of the dense non-aqueous phase liquids and the lithologies that can influence contaminant migration. The 2D resistivity investigation used a capacitively coupled (CC) resistivity method as a reconnaissance tool to locate geophysical anomalies that could be associated with possible areas of creosote contamination. The inversion results of the CC resistivity survey identified resistive anomalies in the subsurface near the eastern and western contaminated source areas. A direct-current (DC) resistivity survey conducted near the CC resistivity survey confirmed the occurrence of subsurface resistive anomalies. The inversion results of the DC resistivity survey were used to define the subsurface distribution of resistivity at each line. Forward modeling was used as an interpretative tool to relate the subsurface distribution of resistivity from four DC resistivity lines to known, assumed, and hypothetical information on subsurface lithologies. The final forward models were used as an estimate of the true resistivity structure for the field data. The forward models and the inversion results of the forward models show the depth, thickness, and extent of strata as well as the resistive anomalies occurring along the four lines and the displacement of strata resulting from the Pecore Fault along two of the four DC resistivity lines. Ten additional DC resistivity lines show similarly distributed shallow subsurface lithologies of silty sand and clay strata. Eight priority areas of resistive anomalies were identified for evaluation in future studies. The interpreted DC resistivity data allowed subsurface stratigraphy to be extrapolated between existing boreholes resulting in an improved understanding of lithologies that can influence contaminant migration.

Scientific Investigations Report↗

Unraveling the dynamics of magmatic CO2 degassing at Mammoth Mountain, California

The accumulation of magmatic CO 2 beneath low-permeability barriers may lead to the formation of CO 2 -rich gas reservoirs within volcanic systems. Such accumulation is often evidenced by high surface CO 2 emissions that fluctuate over time. The temporal variability in surface degassing is believed in part to reflect a complex interplay between deep magmatic degassing and the permeability of degassing pathways. A better understanding of the dynamics of CO 2 degassing is required to improve monitoring and hazards mitigation in these systems. Owing to the availability of long-term records of CO 2 emissions rates and seismicity, Mammoth Mountain in California constitutes an ideal site towards such predictive understanding. Mammoth Mountain is characterized by intense soil CO 2 degassing (up to ∼1000 t d −1 ) and tree kill areas that resulted from leakage of CO 2 from a CO 2 -rich gas reservoir located in the upper ∼4 km. The release of CO 2 -rich fluids from deeper basaltic intrusions towards the reservoir induces seismicity and potentially reactivates faults connecting the reservoir to the surface. While this conceptual model is well-accepted, there is still a debate whether temporally variable surface CO 2 fluxes directly reflect degassing of intrusions or variations in fault permeability. Here, we report the first large-scale numerical model of fluid and heat transport for Mammoth Mountain. We discuss processes (i) leading to the initial formation of the CO 2 -rich gas reservoir prior to the occurrence of high surface CO 2 degassing rates and (ii) controlling current CO 2 degassing at the surface. Although the modeling settings are site-specific, the key mechanisms discussed in this study are likely at play at other volcanic systems hosting CO 2 -rich gas reservoirs. In particular, our model results illustrate the role of convection in stripping a CO 2 -rich gas phase from a rising hydrothermal fluid and leading to an accumulation of a large mass of CO 2 (∼10 7 –10 8 t) in a shallow gas reservoir. Moreover, we show that both, short-lived (months to years) and long-lived (hundreds of years) events of magmatic fluid injection can lead to critical pressures within the reservoir and potentially trigger fault reactivation. Our sensitivity analysis suggests that observed temporal fluctuations in surface degassing are only indirectly controlled by variations in magmatic degassing and are mainly the result of temporally variable fault permeability. Finally, we suggest that long-term CO 2 emission monitoring, seismic tomography and coupled thermal–hydraulic–mechanical modeling are important for CO 2 -related hazard mitigation.

California↗

Methods for estimating annual exceedance probability discharges for streams in Arkansas, based on data through water year 2013

In 2013, the U.S. Geological Survey initiated a study to update regional skew, annual exceedance probability discharges, and regional regression equations used to estimate annual exceedance probability discharges for ungaged locations on streams in the study area with the use of recent geospatial data, new analytical methods, and available annual peak-discharge data through the 2013 water year. An analysis of regional skew using Bayesian weighted least-squares/Bayesian generalized-least squares regression was performed for Arkansas, Louisiana, and parts of Missouri and Oklahoma. The newly developed constant regional skew of -0.17 was used in the computation of annual exceedance probability discharges for 281 streamgages used in the regional regression analysis. Based on analysis of covariance, four flood regions were identified for use in the generation of regional regression models. Thirty-nine basin characteristics were considered as potential explanatory variables, and ordinary least-squares regression techniques were used to determine the optimum combinations of basin characteristics for each of the four regions. Basin characteristics in candidate models were evaluated based on multicollinearity with other basin characteristics (variance inflation factor < 2.5) and statistical significance at the 95-percent confidence level ( p ≤ 0.05). Generalized least-squares regression was used to develop the final regression models for each flood region. Average standard errors of prediction of the generalized least-squares models ranged from 32.76 to 59.53 percent, with the largest range in flood region D. Pseudo coefficients of determination of the generalized least-squares models ranged from 90.29 to 97.28 percent, with the largest range also in flood region D. The regional regression equations apply only to locations on streams in Arkansas where annual peak discharges are not substantially affected by regulation, diversion, channelization, backwater, or urbanization. The applicability and accuracy of the regional regression equations depend on the basin characteristics measured for an ungaged location on a stream being within range of those used to develop the equations.

Arkansas↗