USGS Science⌕ Search

SEARCH · USGS Science

Results for “Surveying and Land Information Systems”

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 757 records · Page 42Linked to original sources
USGS2012

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Aynak mineral district in Afghanistan: Chapter E 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 Aynak mineral district, which has copper 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,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. 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 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 Aynak) and the WGS84 datum. The final image mosaics were subdivided into four overlapping tiles or quadrants because of the large size of the target area. The four image tiles (or quadrants) for the Aynak 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 Aynak study area, five subareas were designated for detailed field investigations (that is, the Bakhel-Charwaz, Kelaghey-Kakhay, Kharuti-Dawrankhel, Logar Valley, and Yagh-Darra/Gul-Darra subareas); these subareas were extracted from the area’s image mosaic and are provided as separate embedded geotiff images.

Kabul;Logar;Paktya;Wardak↗
USGS2012

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↗
USGS2013

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Kunduz mineral district in Afghanistan: Chapter S 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 Kunduz mineral district, which has celestite 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,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 Kunduz) 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 Kunduz 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.

Kunduz Mineral District↗
USGS2012

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↗
USGS2013

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Dudkash mineral district in Afghanistan: Chapter R 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 Dudkash mineral district, which has industrial mineral 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,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 Dudkash) 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 Dudkash 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.

Dudkash Mineral District↗
USGS1985

Environmental systems and management activities on the Kennedy Space Center, Merritt Island, Florida: results of a modeling workshop

In the early 1960's, the National Aeronautics and Space Administration (NASA) began purchasing 140,000 acres on Merritt Island, Florida, in order to develop a center for space exploration. Most of this land was acquired to provide a safety and security buffer around NASA facilities. NASA, as the managing agency for the Kennedy Space Center (KSC), is responsible for preventing or controlling environmental pollution from the Federal facilities and activities at the Space Center and is committed to use all practicable means to protect and enhance the quality of the surrounding environment. The Merritt Island National Wildlife Refuge was established in 1963 when management authority for undeveloped lands at KSC was transferred to the U.S. Fish and Wildlife Service. In addition to manage for 11 Federally-listed threatened and endangered species and other resident and migratory fish and wildlife populations, the Refuge has comanagement responsibility for 19,000 acres of mosquito control impoundments and 2,500 acres of citrus groves. The Canaveral National Seashore was developed in 1975 when management of a portion of the coastal lands was transferred from NASA to the National Park Service. This multiagency jurisdiction on Merritt Island has resulted in a complex management environment. The modeling workshop described in this report was conducted May 21-25, 1984, at the Kennedy Space Center to: (1) enhance communication among the agencies with management responsibilities on Merritt Island; (2) integrate available information concerning the development, management, and ecology of Merritt Island; and (3) identify key research and monitoring needs associated with the management and use of the island's resources. The workshop was structured around the formulation of a model that would simulate primary management and use activities on Merritt Island and their effects on upland, impoundment, and estuarine vegetation and associated wildlife. The simulation model is composed of four connected submodels. The Uplands submodel calculates changes in acres and structural components of vegetation communities resulting from succession, fire, facilities development, and shuttle launch depositions, as well as the quantity and quality of surface runoff and aquifer input to an impoundment and an estuary. The Impoundment submodel next determines water quality and quantity and changes in vegetation resulting from water level manipulation and prescribed burning. The Estuary submodel than determines water quality parameters and acres of seagrass beds. Finally, the Wildlife submodel calculates habitat suitability indices for key species of interest, based on vegetation conditions in the uplands and impoundments and on several hydrologic parameters. The model represents a hypothetical management unit with 2,500 acres of uplands, a 600-acre impoundment, and a 1,500-acre section of estuary. Two management scenarios were run to analyze model behavior. The scenarios differ in the frequency of shuttle launches and prescribed burning, the extent of facilities development, the amount of land disposed waste material applied, and the nature and timing of impoundment water level control. Early in a model development project, the process of building the model is usually of greater benefit than the model itself. The model building process stimulates interaction among agencies, assists in integrating existing information, and helps identify research needs. These benefits usually accrue even in the absence of real predictive power in the resulting model. Open communication occurs among the Federal, State, and local agencies involved with activities on Merritt Island and the agencies have a cooperative working relationship. The workshop provided an opportunity for all of these agencies to meet at one time and have focused discussions on the key environmental and multiagency resource management issues. The workshop framework helped to integrate information and assumptions from a number of disciplines and agencies. This integration occurred in the computer simulation model and among workshop participants as submodel linkages were developed and scenario results discussed. A number of research needs were identified at the workshop during the model building and testing exercises and associated discussions. These needs were based on the informed judgement of researchers and managers familiar with Merritt Island or similar areas, rather than on a comprehensive literature review of sensitivity analysis of the preliminary model developed at the workshop. Some of the needs can be addressed by interpreting the results of completed studies from similar geographic areas as they relate to Merritt Island, while other will require additional research studies on Merritt Island. Major research needs associated with the Upland submodel include behavior of the near-surface aquifer, factors limiting slash pine regeneration, frequency and effects of natural fire on various cover types, cumulative effects of shuttle launches, and fate in upland soils of nitrogen and phosphorous from land applied waste material. Key Impoundment submodel needs include documentation of vegetation changes in response to altered water depth, salinity, and nutrient concentrations and better specification of the functional characteristics of impoundments as chemical filters. Important information gaps identified in the Estuary submodel include a more complete analysis of factors contributing to phytoplankton abundance, evaluation of sources of turbidity other than phytoplankton, and identification and quantification of factors limiting seagrass distribution. Primary research needs associated with the Wildlife submodel include a survey of breeding habitat, production data, and harvest data for mottled ducks; data on the emigration and immigration of juvenile mullet (and other transient fish) in the impoundment; the contribution of various seagrasses to habitat requirements of sea trout; and the effects of dissolved oxygen on survival of juvenile sea trout. Ideally, the modeling workshop process is iterative in nature. Periods between workshops are used for research, data collection, and model refinement. Each workshop integrates information collected since the last workshop and produces a more credible model that is more useful in evaluating management alternatives. Participants felt that continued application of this process would help provide ongoing integration and communication among agencies and would allow each agency's planning and management activities to be viewed within the context of an overall assessment.

Report↗
USGS2013

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Herat mineral district in Afghanistan: Chapter T 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 Herat mineral district, which has barium and limestone 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,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 1,000-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 Herat) 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 Herat 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 Herat study area, one subarea was designated for detailed field investigations (that is, the Barium-Limestone subarea); this subarea was extracted from the area's image mosaic and is provided as separate embedded geotiff images.

Herat↗
USGS2012

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↗
USGS2012

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Kharnak-Kanjar mineral district in Afghanistan: Chapter K 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 Kharnak-Kanjar mineral district, which has mercury 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. 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 1,000-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 Kharnak-Kanjar) 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 Kharnak-Kanjar 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 Kharnak-Kanjar study area, three subareas were designated for detailed field investigations (that is, the Koh-e-Katif Passaband, Panjshah-Mullayan, and Sahebdad-Khanjar subareas); these subareas were extracted from the area's image mosaic and are provided as separate embedded geotiff images.

Farah;Ghor;Daykundi↗
USGS2013

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↗
USGS2012

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Haji-Gak mineral district in Afghanistan: Chapter C 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 Haji-Gak mineral district, which has iron ore 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,2006,2007), 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 co-registered 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 Haji-Gak) and the WGS84 datum. The final image mosaics were subdivided into three overlapping tiles or quadrants because of the large size of the target area. The three image tiles (or quadrants) for the Haji-Gak 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 Haji-Gak study area, three subareas were designated for detailed field investigations (that is, the Haji-Gak Prospect, Farenjal, and NE Haji-Gak subareas); these subareas were extracted from the area's image mosaic and are provided as separate embedded geotiff images.

Bamyan;Parwan;Wardak↗
USGS2012

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 DS consists of the locally enhanced ALOS image mosaics for each of the 24 mineral project areas (referred to herein as areas of interest), whose locality names, locations, and main mineral occurrences are shown on the index map of Afghanistan (fig. 1). 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, 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. PRISM image orthorectification for one-half of the target areas 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). 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 SPARKLE logic, which is described in Davis (2006). Each of the four-band images within each 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 specified radius that was usually 500 m. 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 (either 41 or 42) and the WGS84 datum. Most final image mosaics were subdivided into overlapping tiles or quadrants because of the large size of the target areas. The image tiles (or quadrants) for each area of interest 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. Approximately one-half of the study areas have at least one subarea designated for detailed field investigations; the subareas were extracted from the area's image mosaic and are provided as separate embedded geotiff images.

Data Series↗
USGS2014

Aquifers of Arkansas: protection, management, and hydrologic and geochemical characteristics of groundwater resources in Arkansas

Sixteen aquifers in Arkansas that currently serve or have served as sources of water supply are described with respect to existing groundwater protection and management programs, geology, hydrologic characteristics, water use, water levels, deductive analysis, projections of hydrologic conditions, and water quality. State and Federal protection and management programs are described according to regulatory oversight, management strategies, and ambient groundwater-monitoring programs that currently (2013) are in place for assessing and protecting groundwater resources throughout the State. Physical attributes, groundwater geochemistry, and groundwater quality are described for each of the 16 aquifers of the State. Information in regard to the hydrology and geochemistry of each of the aquifers is summarized from about 550 historical and recent publications. Additionally, more than 8,000 sites with groundwater-quality data were obtained from the U.S. Geological Survey National Water Information System and the Arkansas Department of Environmental Quality databases and entered into a spatial database to investigate distribution and trends in chemical constituents for each of the aquifers. The 16 aquifers of the State were divided into two major physiographic regions of the State: the Coastal Plain Province (referred to as Coastal Plain) of eastern and southern Arkansas, which includes 11 of the 16 aquifers, and the Interior Highlands Division (referred to as Interior Highlands) of western Arkansas, which includes the remaining 5 aquifers. The 11 aquifers in the Coastal Plain consist of various geologic units that are Cenozoic in age and consist primarily of Cretaceous, Tertiary, and Quaternary sands, gravels, silts, and clays. Groundwater in the Coastal Plain represents one of the most valuable natural resources in the State, driving the economic engines of agriculture, while also supplying abundant water for commercial, industrial, and public-supply use. In terms of age from youngest to oldest, the aquifers of the Coastal Plain include Quaternary alluvial aquifers, including the Mississippi River Valley alluvial aquifer (the most important aquifer in Arkansas in terms of volume of use and economic benefits), the Jackson Group (a regional confining unit that served for decades as an important source of domestic supply), and the Cockfield, Sparta, Cane River, Carrizo, Wilcox, Nacatoch, Ozan, Tokio, and Trinity aquifers. The Mississippi River Valley alluvial aquifer accounts for approximately 94 percent of all groundwater used in the State, and the aquifer is used primarily for irrigation purposes. The Sparta aquifer is the second most important aquifer in terms of use, and the aquifer was used in the past dominantly as a source of public and industrial supply, although increasing irrigation use is occurring because of critically declining water levels in the Mississippi River Valley alluvial aquifer. Other aquifers of the Coastal Plain generally are used as important local sources of domestic, industrial, and public supply, in addition to other minor uses. Water quality generally is good for all aquifers of the Coastal Plain, except for elevated iron concentrations and localized areas of high salinity. The high salinity results from intrusion from underlying formations, evapotranspiration processes in areas of low recharge, and inadequate flushing in downgradient areas of residual salinity from deposition in marine environments. Trends in the spatial distribution of individual chemical constituents are related to position along the flow path for most aquifers of the Coastal Plain. These trends include elevated iron and nitrate concentrations with lower pH values and dissolved solids in groundwater from the outcrop areas, transitioning to lower iron and nitrate (related to changes in redox) and higher pH and dissolved solids (dominantly from the dissolution of carbonate minerals) in groundwater downgradient from outcrop areas. Groundwater generally trended from a calcium- to a sodium-bicarbonate water type with increasing cation exchange along the flow path. The Interior Highlands of western Arkansas has less reported groundwater use than other areas of the State, reflecting a combination of factors. These factors include prevalent and increasing use of surface water, less intensive agricultural uses, lower population and industry densities, lesser potential yield of the resource, and lack of detailed reporting. The overall low yields of aquifers of the Interior Highlands result in domestic supply as the dominant use, with minor industrial, public, and commercial-supply use. Where greater volumes are required for growth of population and industry, surface water is the greatest supplier of water needs in the Interior Highlands. The various aquifers of the Interior Highlands generally occur in shallow, fractured, well-indurated, structurally modified bedrock of this mountainous region of the State, as compared to the relatively flat-lying, unconsolidated sediments of the Coastal Plain. In terms of age from youngest to oldest, the aquifers of the Interior Highlands include: the Arkansas River Valley alluvial aquifer, the Ouachita Mountains aquifer, the Western Interior Plains confining system, the Springfield Plateau aquifer, and the Ozark aquifer. Spatial trends in groundwater geochemistry in the Interior Highlands differ greatly from trends noted for aquifers of the Coastal Plain. In the Coastal Plain, the prevalence of long regional flow paths results in regionally predictable and mappable geochemical changes along the flow paths. In the Interior Highlands, short, topographically controlled flow paths (from hilltops to valleys) within small watersheds represent the predominant groundwater-flow system. As such, dense data coverage from numerous wells would be required to effectively characterize these groundwater basins and define small-scale geochemical changes along any given flow path for aquifers of the Interior Highlands. Changes in geochemistry generally were related to rock type and residence time along individual flow paths. Dominant changes in geochemistry for the Ouachita Mountains aquifer and the Western Interior Plains confining system are attributed to rock/water interaction and changes in redox zonation along the flow path. In these areas, groundwater evolves along flow paths from a calcium- to a sodium-bicarbonate water type with increasing reducing conditions resulting in denitrification, elevated iron and manganese concentrations, and production of methane in the more geochemically evolved and strongest reducing conditions. In the Ozark and Springfield Plateau aquifers, rapid influx of surface-derived contaminants, especially nitrogen, coupled with few to no attenuation processes was attributed to the karst landscape developed on Mississippian- and Ordovician-age carbonate rocks of the Ozark Plateaus. Increasing nitrate concentrations are related to increasing agricultural land use, and areas of mature karst development result in higher nitrate concentrations than areas with less karst features.

Arkasas↗
USGS2021

Estimates of public-supply, domestic, and irrigation water withdrawal, use, and trends in the Upper Rio Grande Basin, 1985 to 2015

The Rio Grande flows approximately 670 miles from its headwaters in the San Juan Mountains of south-central Colorado to Fort Quitman, Texas, draining the Upper Rio Grande Basin (URGB) study area of 32,000 square miles that includes parts of Colorado, New Mexico, and Texas. Parts of the basin extend into the United Mexican States (hereafter “Mexico”), where the Rio Grande forms the international boundary between Texas and the State of Chihuahua, Mexico. The URGB was chosen as a focus area study (FAS) for the U.S. Geological Survey (USGS) National Water Census (NWC) as part of the WaterSMART initiative. The objective of the USGS NWC under WaterSMART is to focus on the technical aspects of providing information and tools to stakeholders so that they can make informed decisions on water availability. This report contains water-use withdrawal estimates of groundwater and surface water for public-supply, self-supplied domestic, and irrigation water use for years 1985–2015 at 5-year intervals for the 22 drainage basins at the subbasin 8-digit hydrologic unit code (HUC-8) level. Data for additional categories of self-supplied industrial, mining, livestock, aquaculture, thermoelectric, and hydroelectric water use are provided in the accompanying data release to illustrate total withdrawals for the URGB. The additional category data are provided in this report only for the year 2015. Deliveries of water from public-supply systems to domestic users are included and are the only water-delivery data presented in this report. Consumptive use for irrigation is reported for all HUC-8 subbasins for 2015 and for select HUC-8s in the other years beginning in 1985 (the irrigation category includes irrigation for both crop and golf). Water transported outside of the URGB (interbasin transfers) is not included as part of the withdrawals and are not accounted for in any category of use within the URGB. Estimated total withdrawals for all the water-use categories (including hydroelectric) in 2015 as reported in the USGS compilations in the URGB were 3,152.10 million gallons per day (Mgal/d). Surface water was the dominant source of water used in the URGB, providing about 71 percent of total withdrawals. Nearly all withdrawals were from freshwater sources; there was a small amount of saline groundwater that was used for public supply and self-supplied industrial, which were all reported in Texas. The proportions of total 2015 withdrawals from States in the URGB are 46 percent each in Colorado and New Mexico and 8 percent in Texas. A comparison of 2015 water withdrawals for the URGB—for the categories of public supply, self-supplied domestic, self-supplied industrial, thermoelectric, irrigation, livestock, mining, aquaculture, and hydroelectric—showed that irrigation is the dominant water use, at 74 percent of total withdrawals. Other water-use categories in the URGB that use about 1 percent or greater of the total water use by volume are public supply (9 percent) and self-supplied domestic and aquaculture (each about 1 percent). This report focuses on the higher volume, consumptively used categories of public supply, self-supplied domestic, and irrigation. A discussion on basin population provides context for the categories of public-supply and self-supplied domestic water use. The population in the part of the basin in the United States grew from 1.36 to 2.26 million people between 1985 and 2015. With the city of Ciudad Juarez, Chihuahua, Mexico, included, the total population of the URGB grew from an estimated 2.01 to 3.66 million people between 1985 and 2015. The largest concentrations of population are in New Mexico, Texas, and Chihuahua, with 98 percent of the total number of people in the basin in 1985 and 99 percent of the total in 2015 residing in these states. Albuquerque, El Paso, and Ciudad Juarez are the largest cities in the basin. Total withdrawals for public supply in the URGB averaged 277 Mgal/d from 1985 to 2015. About 60 percent of the URGB total public-supply withdrawals occurred in New Mexico, which averaged 170 Mgal/d. Groundwater provided 92 and 70 percent of the total withdrawals for public supply in 1985 and 2015, respectively. Deliveries to domestic users from public suppliers are reported for all drainage basins and years, and account for part of the total public-supply withdrawals. In the URGB, domestic deliveries from public suppliers increased from 1985 to 1995; since 2005, domestic deliveries from public supply have declined. The total populations served by public suppliers in the URGB increased by 90 percent from 1985 to 2015. In the URGB, more people were served by public-supply systems than were self-supplied, and the percentage of people on public-supply systems ranged from 81 to 92 percent from 1985 to 2015. Total domestic withdrawals in the URGB (deliveries plus self-supply withdrawals) ranged from 177.49 to 234.83 Mgal/d and peaked in 2005. Domestic use decreased from 2005 to 2010 by 17 percent and remained less than 200 Mgal/d in 2015. The per-capita daily use for the entire URGB fluctuated between the reporting years, but overall, domestic per-capita use across the basin has declined 46 percent from 145 gallons per capita daily (gpcd) in 1985 to 79 gpcd in 2015. Total irrigation withdrawals in the URGB had a mean value of 2,767.66 Mgal/d from 1985 to 2015 and withdrawals peaked in 1995 at 3,416.84 Mgal/d. Over the 30-year period, irrigation source water in the URGB has ranged from 69 to 84 percent surface water, and the most surface water diverted in the basin for irrigation was in 1995. Groundwater withdrawals for irrigation have fluctuated but overall decreased by 13 percent between 2005 and 2015. Slightly more than one-half of total irrigation withdrawals within the URGB occurred in Colorado, with a mean of 57 percent from 1985 to 2015. From the peak of water withdrawals in 1995 to the conclusion of this study in 2015, total irrigation withdrawals across the study area decreased by 32 percent. The total number of irrigated lands in the URGB from 1985 to 2015 had a mean of about 800 thousand acres, and more irrigated lands were consistently located in the headwaters of the URGB in the San Luis Valley, Colorado than the remainder of the study basin. In the 30-year period, Colorado had a mean of 68 percent of total irrigated lands whereas irrigated acres in New Mexico had a mean of 26 percent and the remaining 7 percent were in Texas. Since 2000, the number of irrigated acres in the URGB has fluctuated, but overall has decreased by 12 percent. More land was irrigated with surface systems (surface irrigation includes flood, furrow, and gated pipe systems, hereafter collectively termed “surface”) in the URGB than with other irrigation system types. Across the 30-year period, from 62 to 88 percent of total irrigated lands had surface-system irrigation, and surface systems covered a mean of 69 percent of the URGB’s acres. Microirrigation systems, predominantly in New Mexico and Texas, compose 0.2 percent or less of the irrigated acres in the basin and were first reported in 1995. From 1985 to 2015, the surface systems decreased in the basin by about 38 percent, and the number of acres of sprinkler and microirrigation systems increased. Acres irrigated by sprinkler systems (predominately center pivot systems) have increased 179 percent from about 99 thousand acres in 1985 to 275 thousand acres in 2015. In this dataset, the number of sprinkler acres surpassed the number of surface irrigated acres in 2000. Within the San Luis Valley in Colorado, the acreage of surface irrigation has decreased, and sprinkler irrigation has increased over the 30-year period. In the New Mexico part of the URGB, surface irrigation is reported as the dominant system type, where irrigation by surface systems accounts for 97–98 percent of how water is provided to crops. As in New Mexico, crops in Texas are irrigated primarily by surface systems. The mean of the mean simulated actual evapotranspiration (ETa) for crops in 2015 across the basin was highest for durum wheat at an estimated 36.00 inches per year (in/yr), and lowest for vegetables at an estimated 19.48 in/yr. Alfalfa and irrigated grass pastures mean ETa had a mean of 31.4 and 27.58 in/yr, respectively, for the basin. Pecans and peppers, both signature crops in the Rio Grande Basin, each had a mean ETa of 30.67 and 30.38 in/yr of mean. In general, mean ETa values for crops were lower in the HUCs of the Colorado San Luis Valley (13010001, 13010002, 13010003 and 13010004) which are more northerly and at higher elevations. The mean ETa for crops increased in the HUCs that are more southerly and at lower elevations in the basin.

New Mexico↗
USGS2013

Identification of metrics to monitor salt marsh integrity on National Wildlife Refuges in relation to conservation and management objectives

Executive Summary Most salt marshes in the US have been degraded by human activities, and threats from physical alterations, surrounding land-use, species invasions, and global climate change persist. Salt marshes are unique and highly productive ecosystems with high intrinsic value to wildlife, and many National Wildlife Refuges (NWRs) have been established in coastal areas to protect large tracts of salt marsh and wetland-dependent species. Various management practices are employed routinely on coastal NWRs to restore and enhance marsh integrity and ensure ecosystem sustainability. Prioritizing NWR salt marshes for application of management actions and choosing among multiple management options requires scientifically-based methods for assessing marsh condition. Monitoring is integral to structured decision-making (SDM), a formal process for decomposing a decision into its essential elements. Within a natural resource context, SDM involves identifying management objectives, alternative management actions, and expected management outcomes. The core of SDM is a set of criteria for measuring system performance and evaluating management responses. Therefore, use of SDM to frame natural resource decisions leads to logical selection of monitoring attributes that are linked explicitly to management needs. We used SDM to guide selection of variables for monitoring the ecological integrity of salt marshes within the National Wildlife Refuge System (NWRS). Our objectives were to identify indicators of salt marsh integrity that are effective across large geographic regions, responsive to a wide range of threats, and feasible to implement within funding and staffing constraints of the NWRS. In April, 2008, we engaged interdisciplinary experts in a week-long rapid prototyping SDM workshop to define the essential elements of salt marsh management decisions on refuges throughout the northeastern, southwestern, and northwestern US, corresponding to respective Regions 5, 2, and 1 of the US Fish and Wildlife Service (FWS). Through this process we identified measurable attributes for monitoring salt marsh ecosystems that are integrated into conservation practice and target management objectives. The following salt marsh attributes were identified through the SDM process either for describing state condition to determine management needs or for evaluating the achievement of management objectives: historical condition and geomorphic setting; ditch density; surrounding land use; ratio of open water area to vegetation area; rate of pesticide application; environmental contaminant concentration; change in marsh surface elevation relative to sea level rise; tidal range and groundwater level; surface topography; salinity; and species composition and abundance of vegetation, invasive species, invertebrates, nekton, and breeding and wintering birds. The identified attributes were too broadly defined to serve as operational monitoring variables. Therefore, we tested specific metrics for quantifying most of these attributes in summers of 2008 and 2009. The first four attributes in the above list can be characterized by office-based analysis of existing GIS data layers. The remaining attributes require field-based methods for assessment. We were forced to exclude a small number of attributes from field tests due to inconsistent data (pesticide application rate, environmental contaminant concentrations) or requirements that exceeded the scope of this project (change in marsh surface elevation; surface topography; benthic invertebrates; wintering birds). We evaluated potential metrics for evaluating all remaining field attributes. In partnership with NWRS biologists, we tested rapid versus intensive metrics for monitoring field attributes (tidal range and groundwater level; marsh surface elevation; salinity; and species composition and abundance of vegetation, invasive species, nekton, and breeding birds) at coastal refuges throughout FWS Region 5. Seven refuges participated in metric testing in 2008: Rachel Carson (ME), Parker River (MA), Wertheim (NY), E. B. Forsythe (NJ), Bombay Hook (DE), Prime Hook (DE), and Eastern Shore of Virginia Complex (VA). These seven and two additional refuges participated in metric testing in 2009: Rhode Island Complex (RI) and Stewart B. McKinney (CT). We based all field metrics on existing protocols for salt marsh assessment. Sampling locations were determined randomly within delineated marsh study units (MSUs) at each refuge. Detailed field methods are provided in appendices to this report. Measurements for individual metrics were averaged across samples within MSUs during each year of sampling. Each year, correlation or regression analysis was conducted on average measurements across MSUs within each attribute set to identify redundant metrics. Statistical redundancy between a pair of metrics within an attribute set (i.e., correlation or regression slopes with p-values < 0.05) was considered justification for eliminating one of the pair from the regional set of monitoring metrics. Decisions regarding metric elimination versus retention were based on feasibility of monitoring, considering such factors as sampling time, resources required, and potential for regional standardization in implementation. The result of these tests is a reduced suite of monitoring metrics that targets NWRS management decisions and is practicable for implementing on a regional scale. Based on these tests, we recommend the following list of metrics for monitoring integrity of NWRS salt marshes (marsh attribute category is in parentheses): (historical condition and geomorphic setting) position of marsh in the landscape, marsh shape, degree of fill and/or fragmentation, degree of tidal flushing, amount of aquatic edge; (ditch density) ranking of ditch density from none to severe; (surrounding land use) relative proportion of agricultural land in a 150-m buffer around the marsh, relative proportion of natural land in a 150-m buffer around the marsh, relative proportion of natural land in a 1-km buffer around the marsh; (ratio of open water area to vegetation area) ratio of open water to emergent herbaceous wetlands within the marsh; (marsh surface elevation) elevation referenced to NAVD88 in a representative area of the marsh; (tidal range and groundwater level) percent of time the marsh surface is flooded during deployment of a continuous water-level monitor at a representative marsh location, mean depth of surface flooding as measured by a continuous water-level monitor at a representative location; (salinity) salinity measured in surface water; (vegetation community) vegetation species richness using the point-intercept method in 100-m diameter survey plots, percent cover of various marsh community types within 100-m diameter survey plots; (invasive species abundance) percent cover of invasive plant species measured using the point-intercept method in 100-m diameter survey plots; (nekton community) nekton density, nekton species richness, length of Fundulus heteroclitus ; (breeding bird community) abundance of Willets counted per point during standard call-broadcast surveys, summed abundance of tidal marsh obligate species (Clapper Rail, Willet, Saltmarsh Sparrow, Seaside Sparrow) counted per point during standard call-broadcast surveys. Metrics describing the historical condition, geomorphic setting, and broad landscape features can be assessed using existing GIS databases. Our results support use of rapid methods to assess the majority of field metrics; only those used to describe the nekton community must be measured using intensive methods (throw traps or ditch nets, dependant on habitat configuration). Implementation of these metrics for quantitative assessment of NWRS salt marsh integrity in FWS Region 5 requires developing sampling designs for each refuge. Additionally, it is important to determine how the monitoring information will be used within a management context. SDM should be used to complete the analysis of salt marsh management decisions. The next steps would involve 1) prioritizing and weighting the management objectives; 2) predicting responses to individual management actions in terms of objectives and metrics; 3) using multiattribute utility theory to convert all measurable attributes to a common utility scale; 4) determining the total management benefit of each action by summing utilities across objectives; and 5) maximizing the total management benefits within cost constraints for each refuge. This process would allow the optimum management decisions for NWRS salt marshes to be selected and implemented based directly on monitoring data and current understanding of marsh responses to management actions. Monitoring the outcome of management actions would then allow new monitoring data to be incorporated into subsequent decisions.

Connecticut, Delaware, Maine, Massachusetts, New J↗
USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗
USGS2020

Fate and behavior tools related to inland spill response—Workshop on the U.S. Geological Survey’s role in Federal science support

Executive Summary There is a growing body of tools available for science support for determining the fate and behavior of industrial and agricultural chemicals that are rapidly injected (“spilled”) into aquatic environments. A 2-day roundtable-style workshop was held by the U.S. Geological Survey (USGS) in Middleton, Wisconsin, in December 2017 to describe and explore existing Federal science support for spill fate and behavior tools used for inland spills, ongoing and new fate and behavior studies, and science gaps in planning and response tools as part of the USGS Midcontinent Region’s efforts to include spill response as part of its strategic plans. A total of 28 attendees representing a variety of Federal, State, and regional entities presented on programs and tools used in various aspects of spill response. Most programs and tools discussed were for spills in riverine environments but tools and applications for spills in lakes, on land surfaces, in urban storm sewer networks, and groundwater also were discussed. A primary workshop focus was to facilitate communication and increase potential for future collaboration among agencies for inland spill science support. The role and need for more USGS science support within the inland spill community was discussed. Enhanced communication is needed within the USGS and the U.S. Department of the Interior science programs, as well as within and among other agencies that do emergency planning and response. A main conclusion of the workshop was that there are untapped resources of the USGS outlined in the agency’s science strategy that could strengthen science support for fate and behavior tools in inland areas, especially in the Upper Mississippi River, Ohio River, and Great Lakes Basins where large freshwater resources overlap with dense corridors of oil and hazardous substances, with transportation networks, and with large populations centers. Fate and behavior tools are being developed quickly for inland spill response by multiple Federal agencies in partnership with local and regional entities. Applicability of these tools ranges from planning and preparedness, to the early stages of spill response for protection of human life and property, and to the application of monitoring and models to assess the long-term consequences of spills. Key findings from the workshop, with an emphasis on potential further development of USGS science support, include the following: •The national and regional response to spills occurs within an established system that must be respected by all parties involved in spill response. The USGS’s role is to support spill responders who are physically working at a spill scene, deploying booms and using other efforts to contain and recover spilled materials. •The USGS has tools that have been used throughout spill response operations, from early response to recovery and restoration. Developing a more formal role for the USGS to participate in science support for inland spills on a consistent basis is a desired outcome. This will require the USGS to improve internal and external communication and would be best accomplished by assigning one or more coordinator positions within the agency to plan and oversee USGS spill-response efforts. More involvement of the USGS on National and Regional Response Teams, especially in the realm of the Science and Technology Subcommittees, will gofar in increasing external communication and integration of fate and behavior tools. •Rapid response to spills requires modeling and mapping of plumes and associated time-of-travel estimation for a range of stream sizes across the United States. Many existing models use USGS streamgage data and the USGS National Hydrography Dataset. Nearly all existing models would benefit from updated linkages to USGS StreamStats and its soon-to-be released time-of-travel estimates,real-time velocity, stream morphology, and slope data. Integrating USGS tools with those from other agencies could be done to better serve the larger spill response community. • A problem is that existing models to rapidly predict plume extent, as well as more followup/longer-term fate and transport models, can be unknown or unavailable to spill responders. Thus, creating and strengthening linkages among USGS scientists skilled at using these tools is needed to support spill response with the on-scene responders. • Research for inland spill fate and behavior done outside of an immediate spill response can assist with spill planning and preparedness by (1) revealing sites likely to experience spills in the future (high-risk sites) and (2) understanding how a spilled substance might behave under a range of environmental conditions. However, USGS research on this topic has been scarce and subject to funding availability. Examples include the 2010 Line 6B Spill release into the Kalamazoo River in Michigan, where the USGS provided science support for a variety of fate and behavior tools for stream and impoundment environments. A long-term research site in Bemidji, Minnesota, provides important insights into transformations and longevity of spilled oil in groundwater and groundwater-surface water interactions. • Linking stream models to other components of this inland environment, including groundwater, overland flow, and karst, is needed. Stream network data can be linked to underground conduits such as storm sewers and karst groundwater systems. Stream models can also be linked with geospatial data such as that contained in U.S. Environmental Protection Agency’s interactive mapping tools. • The USGS is uniquely qualified to collect water-quality data during spills in the United States because of its many geographically dispersed water science centers, its knowledge and preparedness for flood measurement and documentation, and its cadre of skilled water-quality employees. Rapid-deployment gages, used for floods, could also be used for spills if they included spill-specific sensors. Coordinated expertise at USGS water and environmental science centers can be used for monitoring spill effects and for assessing risk to water quality and ecological communities. • Scientists at the USGS have proven capable of providing science coordination and technical assistance within the Incident Command Structure at the request of the lead on-scene coordinator. This external coordination, as well as internal communication within USGS Water, Hazards, and Ecosystems Mission Areas, could be improved by establishing and naming a USGS spills coordinator. Scott Morlock, Jo Ellen Hinck, and Faith Fitzpatrick are currently (2017) serving in informal coordination roles in addition to their traditional duties.

Open-File Report↗
USGS2022

Refining sources of polychlorinated biphenyls in the Back River watershed, Baltimore, Maryland, 2018–2020

Older urban landscapes present unique and complex stressors to urban streams and their habitats through the introduction of legacy and emerging toxic contaminants. Contaminant sources are often associated with various developed land uses such as older residential areas, active and former industrial sites, contaminated sites, and effluents from municipal wastewater treatment plant discharges. These landscapes have a history of legacy contaminant use such as polychlorinated biphenyls (PCBs) resulting in impacts to sediment and water in these complex environments. Despite the ban of PCBs in new commercial use in 1979, PCB contamination is still widespread in the environment, with many fish consumption advisories throughout the Chesapeake Bay region based on elevated PCBs. Several watersheds in the Baltimore region have mandated reductions in PCBs per total maximum daily loads in tidal waters of the watersheds in order to promote compliance with water quality standards. Some of these mandated reductions (for example, regulated watershed runoff) specified in the total maximum daily loads are the responsibility of the local jurisdictions as part of their phase 1 National Pollutant Discharge Elimination System municipal separate storm sewer system permit. In cooperation with the Baltimore City Department of Public Works and Maryland Department of the Environment, the U.S. Geological Survey and University of Maryland, Baltimore County conducted a study from 2018 to 2020 to refine the sources of PCBs from the City of Baltimore into Back River and to use the results to improve the conceptual site model of PCBs in the Back River watershed. PCB concentrations in the water column of the nontidal streams in Back River watershed are relatively consistent throughout both tributaries, with greater concentrations detected in samples collected from Moores Run but greater loads estimated in samples collected from Herring Run. PCB concentrations measured in the bed sediments and analysis of the flux between sediment porewater (hereafter porewater) and surface water within the tributaries suggest that there are no stationary legacy sources within the stream channels. The bulk of PCB mass entering the system from these nontidal tributaries appears to be introduced primarily during storm events. While only one storm event was sampled and concentrations were quantified only in Herring Run, solids captured during the storm were characterized by increases in PCB mass and overall suspended solids concentrations. Although the bioavailability of the PCB-associated sediment is unknown, this mechanism appears to warrant additional attention to better understand how concentrations vary under different storm conditions and temporally. The importance of contaminated stormwater in loading to Herring Run is further supported by the PCB concentrations in storm drain sediments collected near the tributary, which were present in higher concentrations and were characterized by different homolog signatures compared to that in bed sediments. The observations in the tributaries differed from PCB concentrations and sediment characteristics downstream from the City of Baltimore boundary, in the upper tidal area of the main stem of Back River, particularly at the passive sampler locations BRT–1 and BRT–3. This depositional environment is characterized by higher organic content in sediments and higher concentrations of PCBs in porewater, which result in the possible flux of contaminants from sediment to the water column. This flux is generally opposite of that observed in the nontidal tributaries and the farthest upstream tidal site (BRT–2) and may be a result of the possible settling of sediment particles introduced via suspended solids in stormwater. Despite an observed considerable reduction in overall PCB mass loading to and from the Back River Wastewater Treatment Plant (BRWWTP) (and similar reductions observed in biosolids) compared to the estimates previously reported from 2015, effluent from the BRWWTP continues to be a primary source of PCBs to Back River. The current study confirmed the likeliness of fat, oil, and grease deposits within the miles of sewer pipe as a source of PCBs to the BRWWTP influent. The differences between PCB concentrations in fat, oil, and grease deposits found in pipes (during replacement) compared to that of the BRWWTP suggest that legacy deposits may contain higher PCB concentrations and may act as a source of PCBs to passing sewage, eventually entering the BRWWTP. Variation in freely dissolved concentrations in the sewer system was apparent through the analysis of PCBs in the primary pump stations using passive sampling, with the largest contribution to the influent attributed to a single pump station and associated piping. The contribution of PCBs to Herring Run and Moores Run via sanitary sewer overflows compared to the BRWWTP effluent is negligible, similar to reports from another large urban wastewater treatment plant. Therefore, decreased occurrence of sanitary sewer overflows is not expected to largely decrease PCB loads. Results of this study suggest that targeted, sediment-capture best management practices in Back River watershed could be an effective way to reduce PCB mass loading assuming that deposited contaminated sediments are effectively isolated. Recent studies of some common urban best management practices such as bioretention have shown removal of PCBs within the stormwater control structures. In addition, appropriately timed street sweeping practices with appropriate collection equipment may be an effective way to reduce contaminants such as PCBs from road runoff sources. Reductions in concentrations and mass loading within the sewer system measured in this study compared to that estimated 5 years prior reflect the possible success of ongoing gray infrastructure management actions. Reductions may be attributable to enhanced nutrient reduction upgrades to the BRWWTP and extensive capital improvements and maintenance to the sewer system. This study employed a combined sampling approach and a variety of sampling methods to include low-density polyethylene passive samplers, high-volume water samples, and grab samples of both water and sediment to characterize the PCB inputs to Herring Run, Moores Run, and Back River. Incorporating the passive samplers provided a time-weighted average of the freely dissolved concentration in the surface water, porewater, WWTP influent and effluent, and pump station influent over the deployment period with picogram per liter detection limits. A similar monitoring approach from this study could be implemented within other subwatersheds or municipal separate storm sewer system jurisdictions to assist in refining primary sources of PCBs in order to inform appropriate mitigation approaches.

Maryland↗