USGS Science⌕ Search

SEARCH · USGS Science

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

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↗

Conflict minerals from the Democratic Republic of the Congo—Gold supply chain

The U.S. Geological Survey (USGS) analyzes mineral and metal supply chains to identify and describe major components of material flows from ore extraction, through intermediate forms, to a final product. Supply chain analyses may be used to identify risks to the United States associated with the supply of critical and strategic minerals and metals and to provide greater supply chain transparency so that policymakers have the fact-based information needed to formulate public policy. This fact sheet focuses on the gold supply chain. The USGS National Minerals Information Center (NMIC) has been asked by governmental and non-governmental organizations to provide information about tantalum, tin, tungsten, and gold (collectively known as &ldquo;3TG minerals&rdquo; ) processing facilities worldwide in response to U.S. legislation aimed at identifying and removing the supply chain links associated with the trade of these metals and minerals among armed groups in the Democratic Republic of the Congo (DRC) and adjacent countries. Post-beneficiation processing plants (generally called smelters and refineries) for tantalum, tin, and tungsten (3T) mineral ores and concentrates were identified by company and industry association representatives as being the link in the 3T mineral supply chain through which these minerals can be traced to their source of origin (mine). Tungsten processing plants were the subject of the first fact sheet in a series of USGS reports about 3TG minerals, which was published by the NMIC in August 2014 (Berm&uacute;dez-Lugo, 2014). Background information about historical conditions and the voluntary due diligence of multinational stakeholders for minerals from conflict-affected and high-risk areas is presented in the tungsten fact sheet. The current fact sheet, the fourth and last in the series about 3TG minerals, focuses on the gold supply chain. Processing of the 3T mineral concentrates requires substantial infrastructure and capital and generally is done at relatively few specialized facilities that are not located at the mine site; primary and secondary processors typically are at separate locations. Gold, however, can easily be processed into semi-refined products at or near the mine site and has a high unit value in any form, which allows it to be readily exported through undocumented channels, making it more difficult to track to the mine or region of origin. To put this in perspective, 30 kilograms (66 pounds) of 83 percent pure gold (20 carat) would form a cube measuring 12 centimeters per side (about the size of a small tissue box) and, at a price of $1,200 per ounce, would be worth nearly $1 million. By contrast, the equivalent value of tungsten concentrates would weigh about 45 metric tons (t) (100,000 pounds). Once conflict sourced gold has been combined with gold from other mines and scrap at a refiner, there is no feasible way to distinguish the source of the gold. Thus, once the gold leaves the immediate area of production, it is nearly indistinguishable from gold products mined in other areas.

Fact Sheet↗

Determination of (4-methylcyclohexyl)methanol isomers by heated purge-and-trap GC/MS in water samples from the 2014 Elk River, West Virginia, chemical spill

A heated purge-and-trap gas chromatography/mass spectrometry method was used to determine the cis - and trans -isomers of (4-methylcyclohexyl)methanol (4-MCHM), the reported major component of the Crude MCHM/Dowanol™ PPh glycol ether material spilled into the Elk River upriver from Charleston, West Virginia, on January 9, 2014. The trans -isomer eluted first and method detection limits were 0.16-μg L −1 trans -, 0.28-μg L −1 cis -, and 0.4-μg L −1 Total (total response of isomers) 4-MCHM. Estimated concentrations in the spill source material were 491-g L −1 trans - and 277-g L −1 cis -4-MCHM, the sum constituting 84% of the source material assuming its density equaled 4-MCHM. Elk River samples collected ⩽ 3.2 km downriver from the spill on January 15 had low (⩽2.9 μg L −1 Total) 4-MCHM concentrations, whereas the isomers were not detected in samples collected 2 d earlier at the same sites. Similar 4-MCHM concentrations (range 4.2–5.5 μg L −1 Total) occurred for samples of the Ohio River at Louisville, Kentucky, on January 17, ∼630 km downriver from the spill. Total 4-MCHM concentrations in Charleston, WV, office tap water decreased from 129 μg L −1 on January 27 to 2.2 μg L −1 on February 3, but remained detectable in tap samples through final collection on February 25 indicating some persistence of 4-MCHM within the water distribution system. One isomer of methyl 4-methylcyclohexanecarboxylate was detected in all Ohio River and tap water samples, and both isomers were detected in the source material spilled.

West Virginia↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Zarkashan mineral district in Afghanistan: Chapter G 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 Zarkashan mineral district, which has copper 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 (©JAXA,2006,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. 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 Zarkashan) 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 Zarkashan 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 Zarkashan study area, three subareas were designated for detailed field investigations (that is, the Mine Area, Bolo Gold Prospect, and Luman-Tamaki Gold Prospect subareas); these subareas were extracted from the area’s image mosaic and are provided as separate embedded geotiff images.

Zarkashan Mineral District↗

Hurricane Sandy impacts on coastal wetland resilience

The goal of this research was to evaluate the impacts of Hurricane Sandy on surface elevation trends in estuarine marshes located across the northeast region of the United States from Virginia to Maine using data from an opportunistic (in other words, not strategic) and collaborative network (from here on, an opportunistic network) of surface elevation table-marker horizon (SET-MH) stations. First, we built a data-base of metadata for 965 individual stations from 96 unique geographical locations that included the location, geomorphic setting, and wetland type for each SET-MH station. The dominant estuarine settings included in the analyses were back-barrier lagoonal marshes and emergent marshes along embayments and tidal tributaries. We then calculated prestorm elevation trends to compare to poststorm elevation measurements to determine the storm impact on each station trend. We hypothesized that the effect of Hurricane Sandy on marsh elevation trends would differ by position relative to landfall (right or left) and distance from landfall in southern New Jersey, as both of these variables influence the presence or absence of storm surge as a result of the physical characteristics of tropical cyclones (in other words, strongest winds typically occur to the right of landfall). Storm surge was spatially less extensive and less deep (~1 meter [m]) in marshes located to the left (in other words, south) of landfall compared to marshes located to the right (in other words, north) of landfall where storm surge covered a larger area and was deeper (3–4 m). About 63 percent of 223 eligible stations had a poststorm trend that was similar to the prestorm trend (in other words, less than ±5 millimeters [mm]), indicating little storm impact on elevation trends at those sites. The remaining 37 percent of stations exhibited significant poststorm deviations from the prestorm trend (in other words, greater than ±5 mm). Of these, stations located to the left of landfall had a significant and greater deviation in their elevation trend, and the deviation was more likely to be positive (elevation gain) compared to marshes located to the right of landfall, which had a significant deviation in their elevation trend that was more likely to be negative (elevation loss). This finding is directly related to storm surge impacts on marsh sediment deposition, where deep storm surge (3–4 m) results in sediment deposition in habitats inland of coastal marshes but less so in the marshes themselves. Substrate compaction by the storm surge over-burden may have contributed to elevation loss, but this was not measured because sufficient marker horizon data were not available for analysis. In contrast, to the left of landfall the wind-driven flooding of sediment laden water pushed into the headwaters of rivers and small bays with an ~1 m surge, and resulted in more prevalent sediment deposition on the marsh surfaces and elevation gain. In general, the findings support previous research showing that the physical characteristics of the storm (for example, wind speed, storm surge height, impact angle of landfall) combined with the local wetland conditions (for example, marsh productivity, groundwater level, tide height) are important factors determining a storm’s impact on soil elevation, and that the soil elevation response can vary widely among multiple wetland sites impacted by the same storm and among different storms for the same wetland site. The final objective of this project was to create a framework using metadata from the opportunistic network of SET-MH stations that could be used to develop a strategic monitoring network designed to address specific climate change impacts and related phenomena identified by land managers and stakeholders. We evaluated the spatial distribution and density of SET-MH stations in relation to geographic coverage, marsh setting, availability of public land, and historical storm surge footprints and hurricane return intervals in order to identify gaps in our understanding of risk and our ability to assess it. Analyses revealed that the general geographic coverage of SET-MH stations is limited given the low percentage of marsh patches with stations, low density of stations, the clumped distribution of stations, and the often limited and uneven distribution of stations in wetlands with a high historical frequency of hurricane strikes and storm surge impacts. These findings can be used by managers and planners to inform the creation of a strategic monitoring network that can, in turn, inform management and adaptation plans for coastal resources in the region. Final plan designs will need to consider financial and infrastructural support required for station maintenance, as well as data collection and management over the long term.

Open-File Report↗

Documentation of programs used to determine a wetlands hydroperiod from model-simulated water-surface elevations

A technique has been developed to determine a wetlands hydroperiod by comparing simulated water levels from a ground-water flow model and land- surface elevation data through a geographic information system. The simulated water levels are compared with the land-surface elevation data to determine the height of the water surface above or below land surface for the area of interest. Finally, the hydroperiod is determined for established time periods using criteria specified by the user. The program application requires the use of geographic information system software (ARC/INFO), including the TIN and GRID subsystems of the software. The application consists of an ANSI compatible C program to translate ground- water data output from the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model (MODFLOW) into a format that can be used as input for the geographic information system programs (AML's). The application uses ARC/INFO AML programs and ARC/INFO menu interface programs to create digital spatial data layers of the land surface and water surface and to determine the hydroperiod. The technique can be used to evaluate and manage wetlands hydrology.

Open-File Report↗

Variations on a method for evaluating decadal-scale changes in the groundwater quality of two GAMA coastal study units 2004–14, California GAMA Priority Basin Project

Decadal changes in groundwater quality in two study units on the north-central California coast were evaluated by the Priority Basin Project (PBP) of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. Groundwater samples collected from wells during 2004–05 were compared on a pair-wise basis to samples collected from the same wells during 2014. The data set consisted of paired-samples from 50 public supply wells in 2 GAMA-PBP study units (25 wells each in the North San Francisco Bay and the Monterey Bay and Salinas Valley Basins), with analytical results for 160 water-quality constituents. Statistical analysis was done on grouped results for the 59 constituents that were detected in at least 10 percent of the samples during either sampling period to evaluate decadal-scale change by a step-trend analysis. The data for both of the sampling periods were processed three different ways, resulting in three variations of the paired results to be submitted for statistical analyses. The first evaluation method variation processed data only to facilitate comparison of data when one or both of the results was a non-detection. The second and third variations applied the additional requirement that differences between initial and decadal-sample results exceed a defined threshold to prevent small differences from supporting the conclusion of a step trend. One method for setting the difference threshold between initial-sampling and resampling results is based on criteria used by the GAMA-PBP to determine whether or not replicate results are acceptable. The other difference-threshold-setting method uses a calculated confidence interval around each result based on demonstrated analytical variability for the constituents during each sampling period. Finally, constituents for which decadal-scale changes were statistically significant were identified using the Wilcoxon-Pratt signed-rank test on each of the three evaluation method variations. Step trends were identified by at least 1 of the 3 method variations for 33 constituents. After considering other factors, however, such as water-quality context and the results of quality-control samples, it was concluded that decadal changes were meaningful for 14 constituents in at least 1 of the 2 study units. Constituents for which step trends indicated meaningful increases were dissolved oxygen, total nitrogen, nitrate, orthophosphate, calcium, chloride, sulfate, iron, and lithium. Constituents for which step trends indicated meaningful decreases were temperature, arsenic, lead, the isotopic ratio of carbon-13, and tritium.

California↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Dusar-Shaida mineral district in Afghanistan: Chapter I 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 Dusar-Shaida mineral district, which has copper and tin 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,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’ 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’ local zone (41 for Dusar-Shaida) and the WGS84 datum. The final image mosaics were subdivided into eight overlapping tiles or quadrants because of the large size of the target area. The eight image tiles (or quadrants) for the Dusar-Shaida 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 Dusar-Shaida study area, three subareas were designated for detailed field investigations (that is, the Dahana-Misgaran, Kaftar VMS, and Shaida subareas); these subareas were extracted from the area’ image mosaic and are provided as separate embedded geotiff images.

Farah;Herat↗

Improvement in precipitation-runoff model simulations by recalibration with basin-specific data, and subsequent model applications, Onondaga Lake Basin, Onondaga County, New York

Water-resource managers in Onondaga County, New York, are faced with the challenge of improving the water quality of Onondaga Lake, which has the distinction of being one of the most contaminated lakes in the United States. To assist in this endeavor, during 2003-07 the U.S. Geological Survey (USGS), in cooperation with the Onondaga Lake Partnership, developed a precipitation-runoff model of the 285-square-mile Onondaga Lake Basin with the computer program Hydrological Simulation Program-Fortran (HSPF). The model was intended to provide a tool whereby the processes responsible for the generation of loads of sediment and nutrients that are transported to Onondaga Lake could be better understood. This objective was only partly attained because data for calibration of the model were available from monitoring sites only at or near the mouths of the major tributaries to Onondaga Lake; no calibration data from headwater subbasins, where the loads originated, were available. To address this limitation and thereby decrease the uncertainty in the simulated results that were associated with headwater processes, the USGS conducted a 3-year (2005-08) basinwide study to assess the quality of surface water in the Onondaga Lake Basin. The study quantified the relative contributions of nonpoint sources associated with the major land uses and land covers in the basin and also monitored known sources and presumed sinks of sediment and nutrient loads, which previously had not been evaluated. The use of the newly acquired data to recalibrate the HSPF model resulted in improvements in the simulation of processes in the headwater subbasins, including suspended-sediment, orthophosphate, and phosphorus generation and transport. Simulation of streamflows in small subbasins was improved by adjusting model parameter values to match base flows, storm peaks, and storm recessions more precisely than had been done with the original model. Simulated recessional and low flows were either increased or decreased as appropriate for a given stream, and simulated peak flows generally were lowered in the revised model. The use of suspended-sediment concentrations rather than concentrations of the surrogate constituent, total suspended solids, resulted in increases in the simulated low-flow sediment concentrations and, in most cases, decreases in the simulated peak-flow sediment concentrations. Simulated orthophosphate concentrations in base flows generally increased but decreased for peak flows in selected headwater subbasins in the revised model. Compared with the original model, phosphorus concentrations simulated by the revised model were comparable in forested subbasins, generally decreased in developed and wetland-dominated subbasins, and increased in agricultural subbasins. A final revision to the model was made by the addition of the simulation of chloride (salt) concentrations in the Onondaga Creek Basin to help water-resource managers better understand the relative contributions of salt from multiple sources in this particular tributary. The calibrated revised model was used to (1) compute loading rates for the various land types that were simulated in the model, (2) conduct a watershed-management analysis that estimated the portion of the total load that was likely to be transported to Onondaga Lake from each of the modeled subbasins, (3) compute and assess chloride loads to Onondaga Lake from the Onondaga Creek Basin, and (4) simulate precolonization (forested) conditions in the basin to estimate the probable minimum phosphorus loads to the lake.

Onondaga, New York↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Ghazni1 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 Ghazni1 mineral district, which has spectral reflectance anomalies indicative of clay, aluminum, gold, silver, mercury, and sulfur 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, 2008, 2009), 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 (42 for Ghazni1) and the WGS84 datum. The images for the Ghazni1 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.

Ghazni1 Mineral District↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Ghazni2 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 Ghazni2 mineral district, which has spectral reflectance anomalies indicative of gold, mercury, and sulfur 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, 2008, 2009), 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 (42 for Ghazni2) and the WGS84 datum. The images for the Ghazni2 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.

Ghazni2 Mineral District↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Ghunday-Achin mineral district in Afghanistan, in Davis, P.A, compiler, 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 Ghunday-Achin mineral district, which has magnesite and talc 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,2008,2009), 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 Ghunday-Achin) 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 Ghunday-Achin 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 Ghunday-Achin study area, two subareas were designated for detailed field investigations (that is, the Achin-Magnesite and Ghunday-Mamahel subareas); these subareas were extracted from the area’s image mosaic and are provided as separate embedded geotiff images.

Data Series↗

The next big earthquake in the Bay Area may come sooner than you think—Are you prepared?

Many of us breathed a little easier after October 17, 1989. The Lorna Prieta earthquake, 7.1 on the Richter scale, meant that the big one, talked about for decades, had finally happened. And, bad as it was, we had survived. There are two things wrong with that. First, Lorna Prieta was not the big one. It was a moderately big one, certainly destructive to some parts of the Bay Area, but nowhere near the size of the great San Francisco earthquake of 1906. Second, having an earthquake like Lorna Prieta has little to do with the likelihood of having another one on a different fault, somewhere else in the area. The inevitability of a damaging earthquake still confronts everybody in the Bay Area, and we still risk substantial damage. A new study, released in July 1990 by the United States Geological Survey, says that there is a 67 percent chance of another earthquake the size of Lorna Prieta during the next 30 years and that the quake could strike at any time, including today. In other words, scientists think that a magnitude 7 or larger earthquake is now twice as likely to happen as not to happen. This is a substantial increase, since in 1988, scientists thought the chance for such an earthquake was 50 percent Gust as likely to occur as not to occur) within 30 years. The new report also says that the next one will most likely strike farther north than Lorna Prieta, somewhere between San Jose and Santa Rosa on either side of the Bay. The epicenter of the October 1989 quake was in a sparsely populated area. The next one, according to the study, will likely be centered in a more populated area. During the Lorna Prieta earthquake, shaking was so severe in the Santa Cruz Mountains that a van overturned, treetops snapped off, and many people were thrown to the ground. Because the next one is expected to strike closer to an urban area, it will cause much more damage. Fortunately, there is something we can do about it. By taking actions, such as those described in this booklet, we can drastically reduce the losses and we can make the Bay Area a safer place to live. Earthquake damage is particularly great in certain locations and in certain buildings. Most locations and most modem buildings are relatively safe. By identifying the greatest hazards, we can set priorities for using our limited resources most effectively to reduce them. The choice is ours.

General Interest Publication↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Kundalyan mineral district in Afghanistan: Chapter H 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 Kundalyan mineral district, which has porphyry copper 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 (©JAXA,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 Kundalyan) and the WGS84 datum. The final image mosaics were subdivided into five overlapping tiles or quadrants because of the large size of the target area. The five image tiles (or quadrants) for the Kundalyan 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 Kundalyan study area, three subareas were designated for detailed field investigations (that is, the Baghawan-Garangh, Charsu-Ghumbad, and Kunag Skarn subareas); these subareas were extracted from the area’s image mosaic and are provided as separate embedded geotiff images.

Kundalyan Mineral District↗