USGS ScienceSearch

SEARCH · USGS Science

Results for “Mineral Investigations Resource Map”

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.

318 records · Page 18Linked to original sources

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

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

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

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

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

Hydrogeologic characterization, groundwater chemistry, and vulnerability assessment, Ute Mountain Ute Reservation, Colorado and Utah

The U.S. Geological Survey, in cooperation with the Ute Mountain Ute Tribe (UMUT), initiated a study in 2016 to increase understanding of the hydrogeology and chemistry of groundwater within select areas of the Ute Mountain Ute Reservation (UMUR) in Colorado and Utah, identify vulnerabilities to the system and other natural resources, and outline information needs to aid in the understanding and protection of groundwater resources. The results presented for this study can be used to support the UMUT’s goal of protecting their vital groundwater resources on the UMUR. Hydrogeologic conditions were characterized for the surficial aquifer contained in Quaternary-age unconsolidated surficial deposits and the Dakota aquifer contained in the Cretaceous-age Dakota Sandstone. In the surficial aquifer, median depth to water ranges from about 5.4 to 17.2 feet below land surface in the Farm and Ranch Enterprise area and 11 to 34 feet below land surface in the Towaoc area, and the water table slopes generally southwest or south. A map of depth to the top of the Dakota Sandstone was constructed from existing well data. Depths range from zero in outcrop areas to more than 3,000 feet below land surface on mesas in the southeastern part of the UMUR. Groundwater-chemistry data were collected by the UMUT from 13 springs and 31 wells from 1996 through 2017. Specific conductance was much lower for samples from springs than from wells; median values were 512 and 6,024 microsiemens per centimeter at 25 degrees Celsius, respectively. Spring samples were well oxygenated. A few well samples were anoxic (dissolved oxygen concentrations less than 0.5 milligrams per liter [mg/L]), indicating reducing conditions in the aquifer. About 75 percent of spring samples had fresh water (total dissolved solids concentrations less than 1,000 mg/L), and about 85 percent of well samples had brackish or highly saline water (total dissolved solids concentrations greater than 1,000 mg/L). Water type for springs on the Ute Mountains was calcium bicarbonate. Lower-altitude springs had a calcium-sulfate water type. Most well samples had sodium as the dominant cation, and sulfate, bicarbonate, and chloride as the dominant anions. Fluoride concentrations in about 45 percent of well samples were greater than an agricultural-use standard of 2 mg/L. Nitrate plus nitrite concentrations in most spring and well samples were less than about 1.6 mg/L per liter. Concentrations in samples from wells in the irrigated agricultural area were elevated; the maximum concentration was 78.5 mg/L. About one-half of the trace-element samples had concentrations that were less than laboratory reporting limits. Only aluminum, arsenic, and selenium in spring samples, and boron and selenium in well samples, were detected at concentrations greater than surface-water standards or water-quality standards for agricultural use of groundwater. Only three organic compounds, the pesticides alachlor and atrazine and the volatile organic compound di(2-ethylhexyl) phthalate, were detected in well samples. The Escherichia coli bacteria was detected in 47 and 23 percent of samples from wells and springs, respectively. The E. coli detections included samples from three culturally significant springs, which did not meet the UMUT cultural-use standard of total absence of E. coli. Tritium and carbon-14 were the primary environmental tracers used for interpreting groundwater ages for Lopez 2 Spring and five wells (AP–1, 5000 Block, Cottonwood Spring, Goodknight, and SE Toe). Water from the AP–1 well contained a mixture of pre- and post-1950s recharge. Tritium and carbon-14 recharge ages for Lopez 2 Spring (post-1950s in age), Goodknight and SE Toe wells (pre-1950s in age), and Cottonwood Spring well (primarily pre-1950s in age) are supported by helium-4 data. The helium-4 data for the 5000 Block well are inconsistent with the tritium and carbon-14 age of pre-1950s recharge because of interference caused by high methane concentrations in the water. Springs and surficial deposits are more vulnerable to contamination from anthropogenic chemicals than deeper bedrock wells. Bedrock aquifers are vulnerable in areas where the geologic formations containing the aquifers are exposed at the land surface. Groundwater in deep bedrock aquifers is likely thousands of years old and is not currently affected by present-day land uses. Both shallow and deep groundwater are vulnerable to naturally occurring salts and minerals, such as of total dissolved solids, major ions, nitrate, and trace elements. Effects of a changing climate on water resources and other ecological characteristics of the UMUR could include changes in evapotranspiration, a decrease in snowpack, decreased aquifer recharge and flow of springs, a decrease in soil moisture, and increased occurrence of wildfires and forest mortality. Of particular interest for the UMUT are possible effects of a changing climate on medicinal and culturally important plants and springs Several information needs were identified during this study that would aid in the understanding and protection of groundwater resources on the UMUR. These include well-completion information for bedrock wells, the collection of environmental tracer data at additional wells, the addition of methane and hydrocarbon analysis to well sampling plans, and the resampling of springs and wells that were last sampled in 2002 or earlier.

Colorado

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

Borehole-geophysical and hydraulic investigation of the fractured-rock aquifer near the University of Connecticut Landfill, Storrs, Connecticut, 2000 to 2001

An integrated borehole-geophysical and hydraulic investigation was conducted at the former landfill area near the University of Connecticut in Storrs, Connecticut, where solvents and landfill leachate have contaminated a fractured-bedrock aquifer. Borehole-geophysical techniques and hydraulic methods were used to characterize the site bedrock lithology and structure, fractures, and hydraulic properties. The geophysical and hydraulic methods included conventional logs, borehole imaging, borehole radar, flowmeter under ambient- and stressed hydraulic conditions, and discrete-zone hydraulic testing, sampling, and monitoring. The conventional geophysical-logging methods included caliper, deviation, electromagnetic induction, gamma, specific conductance, and fluid temperature. The advanced methods included optical and acoustic imaging of the borehole wall, heat-pulse flowmeter, and directional radar reflection. Borehole-geophysical methods were used to further define conductive features identified with surface-geophysical methods in the first phase of the investigation. The results of the surface- and borehole-geophysical logging were evaluated in an iterative and integrated manner to develop a conceptual model of ground-water flow at the site. The rock type, foliation, and fractures at the site were characterized from high-resolution optical televiewer (OTV) images of rocks penetrated by the boreholes and were compared to drilling logs and conventional geophysical logs. The rocks are interpreted as fine- to mediumgrained quartz-feldspar-biotite-garnet gneiss and schist with local intrusions of quartz diorite and pegmatite and minor concentrations of sulfide mineralization similar to rocks described as the Bigelow Brook Formation on regional geologic maps. Layers containing high concentrations of sulfide minerals appear as high electrical conductivity zones on electromagnetic-induction and borehole-radar logs. Foliation in the rocks generally strikes to the southwest and northeast, and dips to the northwest and southeast consistent with previous investigations in this area. The orientation of foliation, however, varies locally and with depth in some of the boreholes. These results are consistent with geologic mapping that has identified small-scale folding. The orientations of the transmissive fractures identified in the six boreholes logged for this investigation are similar to the fracture orientations mapped in a previous investigation. Many of these fractures are oriented with a north-northwest strike and have a shallow dip to the west. Other transmissive fractures have a southwest strike and dip at shallow angles to the northwest, and some strike roughly east-west and dip to the north and south. Flowmeter logging was used to identify transmissive fractures and to estimate the hydraulic properties in the boreholes. Ambient down flow was measured in one borehole, and ambient up flow and down flow were measured in another borehole. The other four bedrock boreholes did not have measurable vertical flow. Under low-rate pumping conditions (0.25 to 0.5 gallons per minute), one to three inflow zones were identified in each well. Commonly, fractures that are active under ambient conditions contribute to the well under pumping conditions. The ambient conditions were incorporated into the determination of the relative proportions of transmissivity. Specific capacity and transmissivity were determined for these open-hole low-rate pumping tests. Quasi-steady-state water levels were reached in four of the boreholes, including MW201R, MW204R, MW302R, and W202-NE. When pumped at low-rate conditions for 0.5 to 4 hours, the specific capacity ranged from 0.03 to 0.18 gallons per minute per foot. The open-hole transmissivity estimates ranged from 4.9 to 30 feet squared per day (ft2/d). Open-hole transmissivity was determined for boreholes that did not reach quasi-steady-state conditions under low-rate pumping conditions. Transmissivity was estimated for MW201R, MW202R, and MW203R using non-equilibrium methods, pumping rate, and the transient drawdown data to estimate the open-hole transmissivity. Transmissivity in these boreholes ranged from 0.98 to 3.2 ft 2 /d. The transmissivity and head of individual fractures or zones of fractures were estimated from heat-pulse flowmeter data acquired under ambient and stressed conditions. In the absence of ambient flow, data from two profiles of heat-pulse flowmeter data under two different stressed conditions were used to estimate the transmissivity and head of individual fracture zones. Only two boreholes, MW302R and W202-NE, had sufficient data for these analyses. The estimated transmissivity of individual transmissive zones ranged from 1.2 to 9.2 ft 2 /d. The transmissivity values determined by this numerical simulation method were less than the open-hole estimations, which were 15 and 30 ft 2 /d. Transmissivity also was measured directly over discrete intervals of the borehole using a straddle-packer apparatus and constant-rate pumping tests. Pumping rates were less than or equal to 0.25 gallons per minute. These discretezone single-hole pumping tests were conducted over a short period of time, usually about 30 minutes to 1 hour in duration. Pumping continued until the test zone reached a steady-state water level or until it was determined that the zone could not yield water at the pumped rate. The estimated transmissivity of individual transmissive zones ranged from about 0.21 to 11 ft 2 /d. The zone at a depth of 197 feet in W202-NE was the only zone that had discrete-interval testing with a straddle packer and sufficient heat-pulse flowmeter data for modeling the flow and estimating transmissivity and head. The two methods produced similar results. The straddle-packer method estimated a transmissivity of 4.7 ft 2 /d, and the heat-pulse flowmeter modeling results estimated a transmissivity of 6.9 ft 2 /d. A comparison of the transmissivity estimates indicate estimates typically are within an order of magnitude. The heat-pulse flowmeter methods used in this investigation to determine transmissivity of the boreholes and the individual fractures measure only the upper two or three orders of magnitude of transmissivity. Hence, other fractures in these boreholes permit the movement of water; their transmissivities, however, are lower than the detection limits of the methods that were used for this investigation and very small compared to the transmissive fractures that were studied. The data collected in this investigation were used to design discrete-zone monitoring systems for four of the boreholes used for monitoring. The results of the investigation are useful for refining the conceptual site model of ground-water flow, and for providing critical information for interpreting the results of water-quality sampling.

Connecticut

Ground-water supplies of the Ypsilanti area, Michigan

As of the date of this report (August 1945), the major water users in the Ypsilanti area are: (1) the city of Ypsilanti, (2) the Willow Run bomber plant, built by the Federal Government and operated by the Ford Motor Co., and (3) the war housing project of the Federal Public Housing Authority, designated in this report the Willow Run Townsite. The city, bomber plant, and townsite have required large quantities of water for domestic and industrial uses, and the necessary water supplies have been developed from wells. The Federal Works Agency had the responsibility of deciding whether the existing water facilities were adequate to meet the expected demands and determining the character of any additional public water-supply facilities that might be constructed with Federal assistance. In order to appraise the ground-water resources of the area the Federal Works Agency requested the Geological Survey to investigate the adequacy of the existing supplies and the availability of additional water. The present report is the result of the investigation, which was made in cooperation with the Michigan Geological Survey Division. The water supplies of the three major users are obtained from wells penetrating glacial and associated sands and gravels. Supplies for the city of Ypsilanti and the Willow Run bomber plant are obtained from wells in the valley of the Huron River; the supply for the Willow Run Townsite is obtained from wells penetrating glacial gravels underlying the upland northeast of the valley. The bedrock formations of the area either yield little water to wells or yield water that is too highly mineralized for most uses. The water supply for the bomber plant is obtained from three closely spaced, highly productive wells at the northern edge of the Huron River, a little more than 3 miles southeast of Ypsilanti. The water receives complete treatment in a modern treatment plant. River water also can be treated and has been used occasionally in the winter and spring. The average daily pumpage during periods of maximum production at the bomber plant has been 4.5 to 4.75 million gallons. On June 30, 1945, production of bombers was suspended, and the plant went on a. maintenance basis. The water supply of the bomber-plant well field is replenished by recharge from precipitation and from the Huron River. The evidence shows that recharge from the river is one of the principal sources of water and gives assurance both of the adequacy of the present supply and of the availability of additional water if needed. The safe yield of the three existing wells is estimated to be not less than 6 million gallons per day. The Ypsilanti public water supply is obtained from three tubular wells drilled in 1943, which replaced a number of suction-pumped tubular wells and one large dug well. All the wells penetrate sand and gravel in the bend of the Huron River in the southeastern part of Ypsilanti. The water is treated in a modern treatment plant completed in 1939. The average daily pumpage in million gallons was about 1.68 in 1942, 1.70 in 1943, and 1.66 in 1944. Considerable water was furnished to the Willow Run bomber plant from the Ypsilanti public-supply system during the period from August 1941 through March 1943. The available information indicates that the water pumped from the Ypsilanti well field is replenished by ground-water flow from adjacent stretches of the Huron Valley and from the upland areas outside the valley, from precipitation on the valley in the vicinity of the well field, and possibly from the Huron River. It is believed that sufficient water can be obtained from the well field to meet the expected demand for a considerable time. The safe yield of the present wells is estimated to be not less than 3 million gallons per day, and detailed pumping tests might show that still larger supplies could be safely developed. The water supply of the Willow Run Townsite is obtained from four wells in two well fields about 2 miles apart, one well in the northwest or Prospect and Geddes Road field, and three wells in the southeast or Wiard Road field. The pumpage was originally expected to be 2 to 3 million gallons per day, but it averaged only about 450,000 gallons per day from March 1943 through June 1945. The evidence afforded by logs of wells and by pumping tests indicates that the water-bearing gravels at the townsite are covered by relatively impermeable materials and thus that the rate of recharge is low. However, only relatively small declines in water level have occurred during more than 2 years' operation of the wells, indicating that recharge may occur nearby. The safe yield of the present wells is estimated to be 1.0 to 1.5 million gallons per day, and detailed pumping tests might show it to be somewhat greater. The water supply of the Willow Run Townsite has the smallest potential capacity of the three major ground-water supplies in the area; however, the demand has been relatively small, and no difficulty should arise unless this demand increases greatly. The investigation involved the drilling of 13 test wells to locate additional ground-water supplies, on the assumption that the wartime demand for water in the Ypsilanti area might increase beyond the capacity of the present sources. All 13 wells were drilled at sites selected by the Geological Survey. Two wells on the Willow Run Townsite were drilled by the Federal Public Housing Authority and 11 at other sites by the Federal Works Agency. Records of the wells are given in the report and discussed with respect to the availability of water at the different sites; similarly, the results of the controlled pumping tests made on four of the test wells drilled by the Federal Works Agency are analyzed and discussed with respect to the availability of additional water. The combined results of the test drilling and the pumping tests show that emergency supplies of several million gallons per day can be developed at the sites of the test wells. The best site shown by the test drilling is on the south bank of the Huron River opposite the bomber-plant well field, where a Well with a capacity of several thousand gallons per minute could be constructed if necessary. Three wells in the outwash-filled valley now occupied by Fleming Creek, on two of which pumping tests were made, show that emergency supplies ranging from a few hundred thousand to a million gallons per day could be developed at these test-well sites if necessary. Smaller Supplies ranging from a gallon or two per minute to perhaps 100,000 or 200,000 gallons per day could be developed at the sites of the remaining test Wells. Pumping tests made on the supply wells of the three major water users by the consulting-engineering firms who designed the systems are analyzed and discussed briefly by the authors. The report includes maps and graphs showing the surficial geology of the area, the location of wells, and fluctuations of water level in selected wells. Alsp included are tables giving data on water levels and pumpage, chemical analyses of representative ground waters, and records of wells.

Michigan

Groundwater and surface-water interaction, water quality, and processes affecting loads of dissolved solids, selenium, and uranium in Fountain Creek, near Pueblo, Colorado, 2012–2014

In 2012, the U.S. Geological Survey, in cooperation with the Arkansas River Basin Regional Resource Planning Group, initiated a study of groundwater and surface-water interaction, water quality, and loading of dissolved solids, selenium, and uranium to Fountain Creek near Pueblo, Colorado, to improve understanding of sources and processes affecting loading of these constituents to streams in the Arkansas River Basin. Fourteen monitoring wells were installed in a series of three transects across Fountain Creek near Pueblo, and temporary streamgages were established at each transect to facilitate data collection for the study. Groundwater and surface-water interaction was characterized by using hydrogeologic mapping, groundwater and stream-surface levels, groundwater and stream temperatures, vertical hydraulic-head gradients and ratios of oxygen and hydrogen isotopes in the hyporheic zone, and streamflow mass-balance measurements. Water quality was characterized by collecting periodic samples from groundwater, surface water, and the hyporheic zone for analysis of dissolved solids, selenium, uranium, and other selected constituents and by evaluating the oxidation-reduction condition for each groundwater sample under different hydrologic conditions throughout the study period. Groundwater loads to Fountain Creek and in-stream loads were computed for the study area, and processes affecting loads of dissolved solids, selenium, and uranium were evaluated on the basis of geology, geochemical conditions, land and water use, and evapoconcentration. During the study period, the groundwater-flow system generally contributed flow to Fountain Creek and its hyporheic zone (as a single system) except for the reach between the north and middle transects. However, the direction of flow between the stream, the hyporheic zone, and the near-stream aquifer was variable in response to streamflow and stage. During periods of low streamflow, Fountain Creek generally gained flow from groundwater. However, during periods of high streamflow, the hydraulic gradient between groundwater and the stream temporarily reversed, causing the stream to lose flow to groundwater. Concentrations of dissolved solids, selenium, and uranium in groundwater generally had greater spatial variability than surface water or hyporheic-zone samples, and constituent concentrations in groundwater generally were greater than in surface water. Constituent concentrations in the hyporheic zone typically were similar to or intermediate between concentrations in groundwater and surface water. Concentrations of dissolved solids, selenium, uranium, and other constituents in groundwater samples collected from wells located on the east side of the north monitoring well transect were substantially greater than for other groundwater, surface-water, and hyporheic-zone samples. With one exception, groundwater samples collected from wells on the east side of the north transect exhibited oxic to mixed (oxic-anoxic) conditions, whereas most other groundwater samples exhibited anoxic to suboxic conditions. Concentrations of dissolved solids, selenium, and uranium in surface water generally increased in a downstream direction along Fountain Creek from the north transect to the south transect and exhibited an inverse relation to streamflow with highest concentration occurring during periods of low streamflow and lowest concentrations occurring during periods of high streamflow. Groundwater loads of dissolved solids, selenium, and uranium to Fountain Creek were small because of the small amount of groundwater flowing to the stream under typical low-streamflow conditions. In-stream loads of dissolved solids, selenium, and uranium in Fountain Creek varied by date, primarily in relation to streamflow at each transect and were much larger than computed constituent loads from groundwater. In-stream loads generally decreased with decreases in streamflow and increased as streamflow increased. In-stream loads of dissolved solids and selenium increased between the north and middle transects but generally decreased between the middle and south transects. By contrast, uranium loads generally decreased between the north and middle transects but increased between the middle and south transects. In-stream load differences between transects appear primarily to be related to differences in streamflow. However, because groundwater typically flows to Fountain Creek under low-flow conditions, and groundwater has greater concentrations of dissolved solids, selenium, and uranium than surface water in Fountain Creek, increases in loads between transects likely are affected by inflow of groundwater to the stream, which can account for a substantial proportion of the in-stream load difference between transects. When loads decreased between transects, the primary cause likely was decreased streamflow as a result of losses to groundwater and flow through the hyporheic zone. However, localized groundwater inflow likely attenuated the magnitude by which the in-stream loads decreased. The combination of localized soluble geologic sources and oxic conditions likely is the primary reason for the occurrence of high concentrations of dissolved solids, selenium, and uranium in groundwater on the east side of the north monitoring well transect. To evaluate conditions potentially responsible for differences in water quality and redox conditions, physical characteristics such as depth to water, saturated thickness, screen depth below the water table, screen height above bedrock, and aquifer hydraulic conductivity were compared by using Wilcoxon rank-sum tests. Results indicated no significant difference between depth to water, screen height above bedrock, and hydraulic conductivity for groundwater samples collected from wells on the east side of the north transect and groundwater samples from all other wells. However, saturated thickness and screen depth below the water table both were significantly smaller for groundwater samples collected from wells on the east side of the north transect than for groundwater samples from other wells, indicating that these characteristics might be related to the elevated constituent concentrations found at that location. Similarly, saturated thickness and screen depth below the water table were significantly smaller for groundwater samples under oxic or mixed (oxic-anoxic) conditions than for those under anoxic to suboxic conditions. The greater constituent concentrations at wells on the east side of the north transect also could, in part, be related to groundwater discharge from an unnamed alluvial drainage located directly upgradient from that location. Although the quantity and quality of water discharging from the drainage is not known, the drainage appears to collect water from a residential area located upgradient to the east of the wells, and groundwater could become concentrated in nitrate and other dissolved constituents before flowing through the drainage. High levels of nitrate, whether from anthropogenic or natural geologic sources, could promote more soluble forms of selenium and other constituents by affecting the redox condition of groundwater. Whether oxic conditions at wells on the east side of the north transect are the result of physical characteristics or of groundwater inflow from the alluvial drainage, the oxic conditions appear to cause increased dissolution of minerals from the shallow shale bedrock at that location. Because ratios of hydrogen and oxygen isotopes indicate evaporation likely has not had a substantial effect on groundwater, constituent concentrations at that location likely are not the result of evapoconcentration.

Colorado

Potential effects of energy development on environmental resources of the Williston Basin in Montana, North Dakota, and South Dakota—Water resources

The Williston Basin has been a leading oil and gas producing area for more than 50 years. While oil production initially peaked within the Williston Basin in the mid-1980s, production rapidly increased in the mid-2000s, largely because of improved horizontal (directional) drilling and hydraulic fracturing methods. In 2012, energy development associated with the Bakken Formation was identified as a priority requiring collaboration toward improved timeliness of issuing permits for new wells combined with reasonable measures to maintain environmental quality. Shortly thereafter, the Bakken Federal Executive Group was created to address common challenges associated with energy development. The Bakken Federal Executive Group partner agencies identified a gap in current understanding of the cumulative environmental challenges attributed to energy development throughout the area, resulting in an effort to aggregate scientific data and identify additional research and information needs related to natural resources within areas of energy development in the Williston Basin. As part of this effort, water resources in the area (including groundwater; streams and rivers; and lakes, reservoirs, and wetlands) were characterized and described in terms of physical occurrence, flow characteristics, recharge, water quality, and water use. Similarly, waters produced during energy-development activities also were characterized even though these waters are not considered usable resources within the area. Groundwater resources were characterized by the major hydrogeologic units, or aquifers, identifying the units that supply most groundwater used for domestic, stock, agricultural, and industrial purposes. The groundwater characterization included other deeper hydrogeologic units in the Williston Basin that may be a useable source of water with treatment, have utility as a reservoir for reinjection of produced waters, or be a source of minerals and energy resources. A generalized groundwater budget and flow system identifying the sources of recharge (stream infiltration, precipitation, and movement [leakage] from other aquifers) and the general groundwater flow direction is included for each of the major hydrogeologic units. Rivers and streams within the Williston Basin with 10 or more years of continuous streamflow data were identified. For a subset of these sites, streamflow characteristics, including the monthly and annual mean flow, were generated to identify seasonal and interannual changes in streamflow and thus provide information on the drivers and reliability of streamflow at the seasonal or multiyear scale. Daily streamflow and annual extreme flows (peak and low flow) also were estimated for the subset of sites. The daily streamflow and annual extreme flow values provide information on short-term or extreme events that are relevant to infrastructure design and evaluating spills, leaks, or accidental discharges of water or petroleum products. Surface-water features (lakes, ponds, and wetlands) were classified using the Cowardin system and identified on the National Wetlands Inventory maps generated by the U.S. Fish and Wildlife Service. The spatial distribution of the surface-water features was analyzed by State, county, and specifically in comparison to the Prairie Pothole Region. The proximity of the surface-water features to energy development infrastructure (specifically oil or gas well pads) was evaluated. It was determined that, although oil or gas wells are often near a surface-water feature, most surface-water features do not have wells nearby, with the exception of wells in the Prairie Pothole Region. Water-quality data were aggregated from two data sources: (1) the Water-Quality Portal, sponsored by the U.S. Geological Survey (USGS), U.S. Environmental Protection Agency (EPA), and National Water Quality Monitoring Council; and (2) a data compilation completed as part of the USGS National Water-Quality Assessment project. The Water-Quality Portal integrates publicly available water-quality data from databases maintained by the USGS, EPA, and U.S. Department of Agriculture, including water-quality data from Tribal, State, and local databases. Water-quality data for 15 commonly measured water-quality constituents were aggregated for groundwater, rivers and streams, and lakes and reservoirs. For each aggregated dataset (groundwater, rivers and streams, and lakes and reservoirs), analyses of the water-quality data included summary statistics, maps of spatial distribution of constituent values, boxplots of constituent values by timeframe or hydrogeologic unit, spatial comparisons of site locations and constituent values to petroleum well density, and comparisons of the constituent values measured to EPA drinking-water standards/guidelines. Produced water includes all fluids brought to the surface along with the targeted hydrocarbons as part of the oil and gas exploration and extraction processes. These fluids may include formation water (waters that co-exist with rock/oil/gas), hydraulic fracturing fluids, and other combinations of water and chemicals used during oil and gas well drilling, development, treatments, recompletions, and workovers. Produced water datasets were aggregated from two sources: the USGS National Produced Waters Geochemical database (ver. 2.1) and a series of projects focused specifically on sampling produced water in the Williston Basin from 2010 to 2014. The National Produced Waters Geochemical database was useful for a general understanding of produced-water chemistry. Produced waters are characterized by extreme salinity and contain elevated concentrations of other constituents (including arsenic, barium, cadmium, lead, zinc, radium-226/radium-228, and ammonium) that could negatively affect water and aquatic resources if released. Produced waters also have a generally unique chemical (isotopic) signature that may be useful in tracking water from different geologic units; for example, the oxygen/deuterium and strontium ratio values measured in brine waters from the Bakken Formation are distinct from brines collected from other geologic units in the Williston Basin. Water-use information related to energy production in the area also was aggregated and summarized. The summary of water use is not limited to oil and gas production but includes water used to produce all types of energy resources in the Williston Basin, including coal/lignite, thermoelectric power, oil and gas, hydropower, biomass and biofuels, wind, geothermal, and solar. Each State has its own methods for regulating and reporting water usage within its jurisdiction. These methods can introduce problems when examining water use from sources, such as the Missouri River or Fox Hills aquifer, that are shared across political boundaries. Without the one-to-one match for usage types and amounts used from a water source, it is difficult to develop a comprehensive water budget for the water source being evaluated. A large amount of freshwater is required to prepare a well for oil and gas well production; in some cases, 3 to 7 million gallons of water are needed per well. The EPA estimates that hydraulic fracturing in the Williston Basin uses between 70 to 140 billion gallons per year. Water also is used for myriad other purposes related to ancillary oil and gas extraction. In addition to water used for immediate energy development, the expanded human workforce migrating into the area and other support staff who have moved into the area during the development also use water. Research and information needs were identified that could be relevant in the evaluation of the effects of energy development on water resources. Information needs related to the evaluation of groundwater resources include the following: improved potentiometric-surface maps for glacial units; availability of a uniform stream network digital geographic coverage that spans the international boundary with Canada; enhanced surface-water use information with regards to the gain and loss of streamflow to shallow groundwater, which would increase understanding groundwater and surface-water interactions; and expanded geophysical assessments. Gaps in the availability of streamflow data include the lack of information on ice-jam flooding despite potential for effects to infrastructure (pipelines, roads, and facilities) and an understanding of the cumulative effects of largely undocumented stock and diversion dams. Although this study resulted in the aggregation of a large quantity of water-quality data, the availability of consistently collected, systematically processed and reported data over large parts of the Williston Basin is sparse. Few samples have been analyzed for constituents that may indicate the effect of energy development on water resources. Constituents that could be considered include boron, chloride, bromide, iodine, fluoride, manganese, lithium, radium, strontium isotopes, volatile organic compounds, and isotopes of inorganic ions (such as hydrogen and carbon). Collaboration between Tribal, Federal, State, and local entities to identify a common study design, common monitoring constituents, and consistent sampling locations would generate datasets with broad utility and would likely result in overall cost savings for monitoring over time. Similarly, there is a need for standardized sample collection, processing, laboratory analytical methods, and the collection of ancillary data for produced waters sampling. Additional characterization of the range of chemical, microbial, and isotopic compositions and quantities of “end-member” produced waters, and the collection of time-series datasets to document the changes in produced waters during and after well development also were needs identified during this study. Water-use estimates would be improved through the implementation of comprehensive studies of water use from groundwater and surface-water sources using consistent methodologies across the Williston Basin. The submission of chemical and water data related to hydraulic fracturing collected by the oil and gas industry would add to the quantity of available data. Consistent implementation of regulations and monitoring controls across political boundaries (State, county, and international) would further improve the consistency of data available for the estimates of water use.

Montana, North Dakota, South Dakota

Swatara Creek basin of southeastern Pennsylvania: An evaluation of its hydrologic system

Local concentrations of population in the Swatara Creek basin of Pennsylvania find it necessary to store, transport, and treat water because local supplies are either deficient or have been contaminated by disposal of wastes in upstream areas. Water in the basin is available for the deficient areas and for dilution of the coal-mine drainage in the northern parts and the sewage wastes in the southern parts. Swatara Creek drains 576 square miles just east of Harrisburg, Pa., and is the largest tributary to the Susquehanna River from the north side below Harrisburg. It rises in the southern Pocono Mountains and flows southwestward across the Lebanon Plateau. On an average day Swatara Creek discharges more than 630 million gallons into the Susquehanna River at Middletown, Pa. In a year this amounts to about 23 inches of water over the entire basin and is the residual from an average annual precipitation of 45.5 inches. During an average year the flow in Swatara Creek from the upper third of the basin above Harper Tavern is always greater than 1,300 mgd (million gallons per day) for at least 15 days and is always greater than 25 mgd for at least 350 days. The daily streamflow from the basin averages 1.1 mgd per sq mi, but yields from different areas range from 0.97 to 1.22 mgd per sq mi. These variations are caused chiefly by differences in precipitation and land cover. The area of lowest yield is in the valleys west of Tremont, and the highest yields are in the Upper and Lower Little Swatara Creek subbasins. At high and medium stages the chemical character of the water in the streams is suitable for public and private supplies. At lower stages, defending on the areas and the amounts of contamination by coal-mine drainage and sewage pollution, the natural flow may require some treatment. At low stages the chemical characteristics of the natural flow not affected by man is almost identical with that of the ground water in the area drained by the stream. In general, the total dissolved solids range from about 25 to 400 parts per million and the hardness is as much as about 300 parts per million. The ground-water increment to the base flow of Swatara Creek averages about 240 mgd, or about 8.8 inches annually, for the basin. Generally, ground-water supplies in amounts of less than 0.5 mgd can be developed south of Blue Mountain. Supplies of several million gallons per day have been developed for industrial use from the permeable limestones in the south-central part of the basin. More intensive investigation in other parts of the basin would indicate areas where supplies of more than 0.5 mgd could be developed from properly spaced wells. The chemical character of water from wells depends largely on the host rock. In highly soluble rocks water contains large amount of dissolved solids; in more resistant rocks concentrations are lower. The chemical character of unpolluted ground water generally reflects the composition of the more readily soluble minerals in the local geologic environment. Areas contaminated by septic- tank effluent may have above normal amounts of nitrate and detergent products. Except where polluted, most ground water is suitable for public and industrial uses without extensive treatment. Sites for storage of surface water exist in the part of the basin lying in the valley and ridge area. As much as 30 to 40 percent of the annual flow could be impounded for release as low-flow augmentation for dilution of mine drainage and other wastes in the basin. Low sediment yields of supplying drainage areas would ensure a long life expectancy of reservoirs at these sites. Overbank flooding of the main stem of the Swatara Creek and its tributaries has occurred many times in the past. However, it has not been a hazard because urban development has not encroached on the flood plain. An inundation map of the August 1933 flood provides a basis that urban planners may use to avoid future damage. As water in the Swatara Creek moves downstream to the Susquehanna River, the flow is influenced consecutively by a large annual rainfall on the northern valley and ridge area, the wastes of surface and subsurface coal-mining activities, and less annual rainfall on the part of the basin lying in the Lebanon Plateau area; the flow is supplemented and further influenced by many tributaries and by the industrial and domestic wastes that are carried by these secondary streams. The annual precipitation ranges from 52 inches at the east edge and 49 inches at the west edge of the mountainous part of the basin to about 41 inches at the southwestern part at Middletown. The rainfall generally is adequate during the growing season to mature the crops. The mean annual temperature at Lebanon is about 52°F, and the growing season is about 180 days. In this report the basin has been divided into eight hydrologic zones, leased on runoff, natural use of water, and chemical character of water. Four zones lie in the valley and ridge area, three lie in the Lebanon Plateau area, and one lies in the highland along the southeastern basin boundary. In each of the zones the hydrologic characteristics are virtually the same, but they may be completely different from those in adjacent zones. The boundaries of the zones generally coincide with boundaries between geologic formations, and the areas in each zone include rocks of similar influence on water. Streams in zone 4 at the northeast edge of the plateau have the highest average surface runoff from 1.2 to 1.1 mgd per sq mi whereas those in zone 2 at the northwest edge of the valley and ridge area have the lowest, about 1.0 mgd. Streams in zone 8, along the southeast edge of the basin, have the largest sustained low-flow yield, about 0.26 to 0.19 mgd per sq mi; those in zone 5 overlying the Martinsburg Shale east of Harrisburg have the smallest sustained low-flow yields, 0.03 to 0.01 mgd. Streams in the limestone area of zone 7 have the greatest range in low-flow yields in any one zone from 0.60 to 0 mgd per sq mi. Low-flow yields in zones 1 through 4 range from 0.13 to 0.03 mgd per sq mi. Surface flows from zones 1 and 2 are generally acidic and contain high concentrations of sulfate, iron, and total dissolved solids especially where contaminated with mine wastes. Surface flows from zones 3 and 4 are dilute, slightly alkaline, and suitable for public water supplies. Surface flows from zones 5, 6, and 7 are alkaline and contain moderate concentrations of dissolved solids with waters of highest hardness occurring in zone 7. Surface flows from zone 8 are dilute to moderately mineralized and are relatively high in silica concentration. Nitrate concentrations are high in surf Fee flows below sewage outfalls and in ground water contaminated by septic tank effluent and industrial wastes. Average annual sediment yields of 550 to 650 tons per square mile are characteristic of zones 1 and 2 where strip mining has destroyed the forest cover and coal culm is carried into the streams. From agricultural lands on the Martinsburg Shale in zones 5 and 6, annual sediment yields range from 300 to 350 tons per square mile; but from agricultural lands on the siliceous rocks in zone 8 and zones 3 and 4 in the valley and ridge area, the sediment yield ranges from 200 to 250 tons annually per square mile. Lowest annual sediment yields in the basin are in the forested areas of siliceous rocks in zones 2, 3, 4, and 5, and in the sinkhole topography of the limestones in zone 7 where the yield ranges from 30 to 35 tons and 50 to 60 tons per square mile, respectively. The amount of ground water that can be developed in the basin is dependent on the ability of the underlying rocks to yield water to wells. More than 300 gpm (gallons per minute) can be obtained from wells in alluvial materials in the valley bottoms and in some of the limestones where large solution channels and fractures are penetrated by the wells. From 50 to 300 gpm can be obtained from wells in loosely cemented sandstones and in fractured limestones. From 10 to 50 gpm can be developed from wells in the shales and harder sandstones. The most dense rocks will yield from 1 to 10 gpm from fractures and crevices. Most wells yield water from the upper 350 feet of the formation, for this part contains the most fractures or solution channels. Studies show that the velocity at which a contaminant will move downstream in the basin is related to the discharge of the stream at the time. At a stream discharge of about 400 mgd at Pine Grove, a contaminant in Swatara Creek would require about 40 hours to move from Pine Grove to Middletown. As a result of dispersion and dilution, the maximum concentration of the contaminant at Middletown would be less than 20 percent the concentration at Pine Grove under these conditions. An evaluation of the availability of water in the basin indicates that about I,239 mgd enters as precipitation, 630 mgd leaves as streamflow, 580 mgd is evaporated and transpired, and 56 mgd is diverted for use by man. Not all the diversions for man's use are lost to the basin, as about 27 mgd is returned as sewage for reuse. About one-fourth of the waste water is returned to the ground and the remainder to stream drainageways. Of that diverted by man, 11.6 mgd is used for public supply and 44.4 mgd for industrial and private supplies. Diversions of streamflow furnish 86 percent of the public supply and 27 percent of the industrial supply, and ground-water sources yield the remainder. Municipal and private sewage treatment plants are upgrading the waste water in many places, but no provisions are being made for treatment other than natural dilution and assimilation for the 15 mgd of coal-nine drainage in the northern part of the basin. Technology for economic treatment of mine water is not available at this time, although research in this field is being done. Urbanization eastward from Harrisburg and around Lebanon has increased the population density of the basin. Densities of 500 people per square mile and water use exceeding 2.0 mgd per sq mi can be expected in the future. By the year 2000 the population of the basin may increase 60 percent; and if the per capita rate of use increases 0.5 percent per year the domestic requirements for water will be about two times the present use, or 23 mgd. Similarly, if the present 1:4 ratio of domestic use to industrial use of water continues, at least 89 mgd will be needed for industry in the future. Although an increase to twice the present use of water can be foreseen, or 112 mgd, water for the dilution and assimilation of wastes from treatment systems are not included. Providing water for dilution of wastes from treatment plants has not been a problem, but in the future the amounts needed for this purpose will be greater as the population increases. As water becomes more valuable, treatment of sewage wastes to reduce the biochemical-oxygen-demand load by at least 80 to 90 percent will be necessary to conserve water for more productive uses. As much as 100 mgd may be needed for waste dilution in the basin by year 2000. The present trends in suburban and light industrial development will probably persist in the basin. Problems arising through changes in economic value of water, conflicts in use, and alternatives in development are typical of those confronting the manager of a water-resource system.

Pennsylvania