Geology ReportsSearch

SEARCH · Geology Reports

Results for “Series in Remote Sensing”

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 307 records · Page 17Linked to original sources

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Panjsher Valley mineral district in Afghanistan: Chapter M in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts in Afghanistan

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

Panjsher Valley Mineral District

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 South Bamyan mineral district in Afghanistan

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

South Bamyan Mineral District

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Khanneshin mineral district in Afghanistan: Chapter A in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts in Afghanistan

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

Helm;Nimroz

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

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

Ghazni1 Mineral District

System characterization report on the Satellogic NewSat multispectral sensor

Executive Summary This report addresses system characterization of Satellogic’s NewSat satellite (also known as ÑuSat) and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the specific sensing system being evaluated; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions. Since 2016, Satellogic has launched 17 NewSat satellites. All NewSat satellites have four-band imagers with a 1-meter (m) ground sample distance, and values in pixels are identical to values in meters. All NewSats have been launched into Sun-synchronous orbits of about 475 kilometers, with inclinations of about 97.5 degrees. The satellites have expected lifetimes of about 3 years. More information on the Satellogic satellites and sensors is available in the “2020 Joint Agency Commercial Imagery Evaluation—Remote Sensing Satellite Compendium” and from the manufacturer at https://satellogic.com/ . The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior), radiometric, and spatial performances. Results of these analyses indicate that the NewSat satellites have an interior geometric performance in the range of −0.119 (−0.119 pixel) to 0.020 m (0.020 pixel) in easting and −0.148 (−0.148 pixel) to 0.014 m (0.014 pixel) in northing in band-to-band registration, an exterior geometric performance of −9.04 (−9.04 pixels) to −5.84 m (−5.84 pixels) in easting and 1.25 (1.25 pixels) to 3.11 m (3.11 pixels) in northing offset in comparison to Sentinel-2, an exterior geometric performance using ground control points of a 6.5-m circular error (95 percent), a radiometric performance in the range of 0.034 to 0.081 in offset and 0.652 to 0.808 in slope, and a spatial performance in the range of 1.61 to 1.76 pixels for full width at half maximum, with a modulation transfer function at a Nyquist frequency in the range of 0.081 to 0.138.

Open-File Report

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

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

Farah Mineral District

Geostatistical estimation of signal-to-noise ratios for spectral vegetation indices

In the past 40 years, many spectral vegetation indices have been developed to quantify vegetation biophysical parameters. An ideal vegetation index should contain the maximum level of signal related to specific biophysical characteristics and the minimum level of noise such as background soil influences and atmospheric effects. However, accurate quantification of signal and noise in a vegetation index remains a challenge, because it requires a large number of field measurements or laboratory experiments. In this study, we applied a geostatistical method to estimate signal-to-noise ratio (S/N) for spectral vegetation indices. Based on the sample semivariogram of vegetation index images, we used the standardized noise to quantify the noise component of vegetation indices. In a case study in the grasslands and shrublands of the western United States, we demonstrated the geostatistical method for evaluating S/N for a series of soil-adjusted vegetation indices derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. The soil-adjusted vegetation indices were found to have higher S/N values than the traditional normalized difference vegetation index (NDVI) and simple ratio (SR) in the sparsely vegetated areas. This study shows that the proposed geostatistical analysis can constitute an efficient technique for estimating signal and noise components in vegetation indices.

ISPRS Journal of Photogrammetry and Remote Sensing

A cross comparison of spatiotemporally enhanced springtime phenological measurements from satellites and ground in a northern U.S. mixed forest

Cross comparison of satellite-derived land surface phenology (LSP) and ground measurements is useful to ensure the relevance of detected seasonal vegetation change to the underlying biophysical processes. While standard 16-day and 250-m Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation index (VI)-based springtime LSP has been evaluated in previous studies, it remains unclear whether LSP with enhanced temporal and spatial resolutions can capture additional details of ground phenology. In this paper, we compared LSP derived from 500-m daily MODIS and 30-m MODIS-Landsat fused VI data with landscape phenology (LP) in a northern U.S. mixed forest. LP was previously developed from intensively observed deciduous and coniferous tree phenology using an upscaling approach. Results showed that daily MODIS-based LSP consistently estimated greenup onset dates at the study area (625 m × 625 m) level with 4.48 days of mean absolute error (MAE), slightly better than that of using 16-day standard VI (4.63 days MAE). For the observed study areas, the time series with increased number of observations confirmed that post-bud burst deciduous tree phenology contributes the most to vegetation reflectance change. Moreover, fused VI time series demonstrated closer correspondences with LP at the community level (0.1-20 ha) than using MODIS alone at the study area level (390 ha). The fused LSP captured greenup onset dates for respective forest communities of varied sizes and compositions with four days of the overall MAE. This study supports further use of spatiotemporally enhanced LSP for more precise phenological monitoring.

Wisconsin

Groundwater levels, geochemistry, and water budget of the Tsala Apopka Lake system, west-central Florida, 2004–12

Tsala Apopka Lake is a complex system of lakes and wetlands, with intervening uplands, located in Citrus County in west-central Florida. It is located within the 2,100 square mile watershed of the Withlacoochee River, which drains north and northwest towards the Gulf of Mexico. The lake system is managed by the Southwest Florida Water Management District as three distinct “pools,” which from upstream to downstream are referred to as the Floral City Pool, Inverness Pool, and Hernando Pool. Each pool contains a mixture of deep-water lakes that remain wet year round, ephemeral (seasonal) ponds and wetlands, and dry uplands. Many of the major deep-water lakes are interconnected by canals. Flow from the Withlacoochee River, when conditions allow, can be diverted into the lake system. Flow thorough the canals can be used to control the distribution of water between the three pools. Flow in the canals is controlled using structures, such as gates and weirs. Hydrogeologic units in the study area include a surficial aquifer consisting of Quaternary-age sediments, a discontinuous intermediate confining unit consisting of Miocene- and Pliocene-age sediments, and the underlying Upper Floridan aquifer, which consists of Eocene- and Oligocene-age carbonates. The fine-grained quartz sands that constitute the surficial aquifer are generally thin, typically less than 25 feet thick, within the vicinity of Tsala Apopka Lake. A thin, discontinuous, sandy clay layer forms the intermediate confining unit. The Upper Floridan aquifer is generally unconfined in the vicinity of Tsala Apopka Lake because the intermediate confining unit is discontinuous and breached by numerous karst features. In the study area, the Upper Floridan aquifer includes the upper Avon Park Formation and Ocala Limestone. The Ocala Limestone is the primary source of drinking water and spring flow in the area. The objectives of this study are to document the interaction of Tsala Apopka Lake, the surficial aquifer, and the Upper Floridan aquifer; and to estimate an annual water budget for each pool and for the entire lake system for 2004–12. The hydrologic interactions were evaluated using hydraulic head and geochemical data. Geochemical data, including major ion, isotope, and age-tracer data, were used to evaluate sources of water and to distinguish flow paths. Hydrologic connection of the surficial environment (lakes, ponds, wetlands, and the surficial aquifer) was quantified on the basis of a conceptualized annual water-budget model. The model included the change in surface water and groundwater storage, precipitation, evapotranspiration, surface-water inflow and outflow, and net groundwater exchange with the underlying Upper Floridan aquifer. The control volume for each pool extended to the base of the surficial aquifer and covered an area defined to exceed the maximum inundated area for each pool during 2004–12 by 0.5 foot. Net groundwater flow was computed as a lumped value and was either positive or negative, with a negative value indicating downward or lateral leakage from the control volume and a positive value indicating upward leakage to the control volume. The annual water budget for Tsala Apopka Lake was calculated using a combination of field observations and remotely sensed data for each of three pools and for the composite three pool area. A digital elevation model at a 5-foot grid spacing and bathymetric survey data were used to define the land-surface elevation and volume of each pool and to calculate the changes in inundated area with change in lake stage. Continuous lake-stage and groundwater-level data were used to define the change in storage for each pool. The rainfall data used in the water-budget calculations were based on daily radar reflectance data and measured rainfall from weather stations. Evapotranspiration was computed as a function of reference evapotranspiration, adjusted to actual evapotranspiration using a monthly land-cover coefficient (based on evapotranspiration measurements at stations located in representative landscapes). Surface-water inflows and outflows were determined using stage data collected at a series of streamgages installed primarily at the water-control structures. Discharge was measured under varying flow regimes and ratings were developed for the water-control structures. The discharge data collected during the study period were used to calibrate a surface-water flow model for 2004–12. Flows predicted by the model were used in the water-budget analysis. Net groundwater flow was determined as the residual term in the water-budget equation. The results of the water-budget analysis indicate that rainfall was the largest input of water to Tsala Apopka Lake, whereas evapotranspiration was the largest output. For the 2004–12 analysis period, surface-water inflow accounted for 11 percent of the inputs, net groundwater inflow accounted for 1 percent of inputs (annual periods with positive net groundwater flow were included as inputs, while annual periods with negative net groundwater flow were counted as outputs), and rainfall accounted for the remaining 88 percent. For the same period, the outputs consisted of 2 percent surface-water outflow, 12 percent net groundwater outflow, and 86 percent evapotranspiration. Net groundwater inflows and surface-water/groundwater storage were negligible during the water-budget period but could be important components of the budget in individual years. The net groundwater flow was negative (downward) for 8 out of the 9 years modeled (2004–12), indicating that the Tsala Apopka Lake study area was primarily a recharge area for the underlying Upper Floridan aquifer during this time period. Groundwater-level elevation in paired wells (adjacent wells completed in the surficial aquifer and Upper Floridan aquifer) typically was higher in the surficial aquifer than the Upper Floridan aquifer. However, hydraulic head data indicate that the surficial aquifer often has discharge potential to the surface-water system, especially in the low lying areas near the major lakes. Surficial-aquifer water levels were often higher than lake stages, especially during wet periods, which is likely an indication of aquifer-to-lake seepage in these areas. East of the major lakes, hydraulic head data were nearly equal in the surficial aquifer and Upper Floridan aquifer, which is an indication that the Upper Floridan aquifer is unconfined. Based on deuterium and oxygen stable isotope data collected in December 2011 and December 2012, there was no evidence of recharge to the Upper Floridan aquifer from the wetlands east of the major lakes; aquifer isotopic ratios did not indicate an enriched source, which is typical of lake and wetland sources. West of the major lakes, there was evidence of enriched isotopic ratios in water samples from the Upper Floridan aquifer. Differences in hydraulic head at paired wells in the surficial aquifer and Upper Floridan aquifer indicated that the surficial aquifer has the potential to recharge the Upper Floridan aquifer in the western part of the pools and west of the major lakes.

Florida

High-resolution digital elevation model of lower Cowlitz and Toutle Rivers, adjacent to Mount St. Helens, Washington, based on an airborne lidar survey of October 2007

The lateral blast, debris avalanche, and lahars of the May 18th, 1980, eruption of Mount St. Helens, Washington, dramatically altered the surrounding landscape. Lava domes were extruded during the subsequent eruptive periods of 1980–1986 and 2004–2008. More than three decades after the emplacement of the 1980 debris avalanche, high sediment production persists in the Toutle River basin, which drains the northern and western flanks of the volcano. Because this sediment increases the risk of flooding to downstream communities on the Toutle and lower Cowlitz Rivers, the U.S. Army Corps of Engineers (USACE), under the direction of Congress to maintain an authorized level of flood protection, continues to monitor and mitigate excess sediment in North and South Fork Toutle River basins to help reduce this risk and to prevent sediment from clogging the shipping channel of the Columbia River. From October 22–27, 2007, Watershed Sciences, Inc., under contract to USACE, collected high-precision airborne lidar (light detection and ranging) data that cover 273 square kilometers (105 square miles) of lower Cowlitz and Toutle River tributaries from the Columbia River at Kelso, Washington, to upper North Fork Toutle River (below the volcano's edifice), including lower South Fork Toutle River. These data provide a digital dataset of the ground surface, including beneath forest cover. Such remotely sensed data can be used to develop sediment budgets and models of sediment erosion, transport, and deposition. The U.S. Geological Survey (USGS) used these lidar data to develop digital elevation models (DEMs) of the study area. DEMs are fundamental to monitoring natural hazards and studying volcanic landforms, fluvial and glacial geomorphology, and surface geology. Watershed Sciences, Inc., provided files in the LASer (LAS) format containing laser returns that had been filtered, classified, and georeferenced. The USGS produced a hydro-flattened DEM from ground-classified points at Castle and Coldwater Lakes. Final results averaged about two laser last-return points per square meter. As reported by Watershed Sciences, Inc., vertical accuracy is 10 centimeters (cm) at the 95-percent confidence interval on bare road surfaces; however, over natural terrain, USGS found vertical accuracy to be 10–50 cm. This USGS data series contains the bare-earth lidar data as 1- and 10-meter (m) resolution Esri grid files. Digital-elevation data can be downloaded (1m_DEM.zip and 10m_DEM.zip), as well as a 1-m resolution hillshade image with pyramids (1m_hillshade.zip). These geospatial data files require geographic information system (GIS) software for viewing.

Washington

High-resolution digital elevation model of Mount St. Helens crater and upper North Fork Toutle River basin, Washington, based on an airborne lidar survey of September 2009

The lateral blast, debris avalanche, and lahars of the May 18th, 1980, eruption of Mount St. Helens, Washington, dramatically altered the surrounding landscape. Lava domes were extruded during the subsequent eruptive periods of 1980–1986 and 2004–2008. More than three decades after the emplacement of the 1980 debris avalanche, high sediment production persists in the North Fork Toutle River basin, which drains the northern flank of the volcano. Because this sediment increases the risk of flooding to downstream communities on the Toutle and Cowlitz Rivers, the U.S. Army Corps of Engineers (USACE), under the direction of Congress to maintain an authorized level of flood protection, built a sediment retention structure on the North Fork Toutle River in 1989 to help reduce this risk and to prevent sediment from clogging the shipping channel of the Columbia River. From September 16–20, 2009, Watershed Sciences, Inc., under contract to USACE, collected high-precision airborne lidar (light detection and ranging) data that cover 214 square kilometers (83 square miles) of Mount St. Helens and the upper North Fork Toutle River basin from the sediment retention structure to the volcano's crater. These data provide a digital dataset of the ground surface, including beneath forest cover. Such remotely sensed data can be used to develop sediment budgets and models of sediment erosion, transport, and deposition. The U.S. Geological Survey (USGS) used these lidar data to develop digital elevation models (DEMs) of the study area. DEMs are fundamental to monitoring natural hazards and studying volcanic landforms, fluvial and glacial geomorphology, and surface geology. Watershed Sciences, Inc., provided files in the LASer (LAS) format containing laser returns that had been filtered, classified, and georeferenced. The USGS produced a hydro-flattened DEM from ground-classified points at Castle, Coldwater, and Spirit Lakes. Final results averaged about five laser last-return points per square meter. As reported by Watershed Sciences, Inc., vertical accuracy is 10 centimeters (cm) at the 95-percent confidence interval on bare road surfaces; however, over natural terrain, USGS found vertical accuracy to be 10–50 cm. This USGS data series contains the bare-earth lidar data as 1- and 10-meter (m) resolution Esri grid files. Digital-elevation data can be downloaded (1m_DEM.zip and 10m_DEM.zip), as well as a 1-m resolution hillshade image with pyramids (1m_hillshade.zip). These geospatial data files require geographic information system (GIS) software for viewing.

Washington

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

Modeling PSInSAR time series without phase unwrapping

In this paper, we propose a least-squares-based method for multitemporal synthetic aperture radar interferometry that allows one to estimate deformations without the need of phase unwrapping. The method utilizes a series of multimaster wrapped differential interferograms with short baselines and focuses on arcs at which there are no phase ambiguities. An outlier detector is used to identify and remove the arcs with phase ambiguities, and a pseudoinverse of the variance-covariance matrix is used as the weight matrix of the correlated observations. The deformation rates at coherent points are estimated with a least squares model constrained by reference points. The proposed approach is verified with a set of simulated data.

IEEE Transactions on Geoscience and Remote Sensing

System characterization report on the German Aerospace Center (DLR) Earth Sensing Imaging Spectrometer (DESIS)

Executive Summary This report addresses system characterization of the German Aerospace Center (DLR) Earth Sensing Imaging Spectrometer (DESIS) and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present the methodology and procedures for characterization and the technical and operational information about the specific sensing system being evaluated. These reports also provide a description of data measurements, data retention practices, and data analysis results and provide system characterization conclusions. In partnership with Teledyne Brown Engineering, DLR built the DESIS hyperspectral instrument, which Teledyne Brown Engineering then integrated onto its International Space Station-based imaging platform, the Multi-User System for Earth Sensing. DLR developed the processing software and, together with Innovative Imaging and Research, completes the validation and calibration of the data products. DESIS was launched in 2018, and the data are used for scientific research in atmospheric physics and Earth sciences. The DESIS sensor contributes to the scientific and commercial utilization of the International Space Station and helps to further hyperspectral remote sensing technologies for future satellites. More information on DLR satellites and sensors is included within the “2020 Joint Agency Commercial Imagery Evaluation—Remote Sensing Satellite Compendium” and at https://www.dlr.de/DE/Home/home_node.html . The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior), radiometric, and spatial performances. Results of these analyses indicate that DESIS has an interior geometric performance of less than a 3.30-meter (less than 0.11 pixel) root mean square error in band-to-band registration, an exterior geometric performance in the range of a 2.40- (0.08 pixel) to 17.40-meter (0.58 pixel) offset in comparison to the Landsat 8 Operational Land Imager, a radiometric performance in the range of −0.013 to 1.011 (offset and slope), and a spatial performance for band 130 of 1.5 pixels at full width at half maximum, with a modulation transfer function at a Nyquist frequency of 0.167.

Open-File Report

System characterization report on the BlackSky Global multispectral sensor

Executive Summary This report addresses system characterization of the BlackSky Global satellites and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the specific sensing system being evaluated; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions. The BlackSky Global satellites are three-band multispectral imagers (red, green, and blue multispectral bands plus a panchromatic band) with a 0.8- to 0.9-meter (m) pixel ground sample distance for the assessed satellites. BlackSky Global satellites 9 and 12–17 were launched in March and December 2021, respectively, into a Sun-synchronous orbit of 430–450 kilometers with an inclination of 42–53 degrees and a swath width of 6 kilometers at nadir. Each Global satellite has an expected lifetime of about 3 years. More information on the BlackSky Global satellites is available in the “Land Remote Sensing Satellites Online Compendium” ( https://calval.cr.usgs.gov/apps/compendium ) and from BlackSky at Real-Time Space-Based Intelligence (blacksky.com) The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior) and spatial performances. Results of these analyses indicate that the assessed BlackSky Global satellites have an interior geometric performance in the range of −0.011 m (−0.012 pixel) to 0.007 m (0.008 pixel) in easting and −0.018 m (−0.020 pixel) to 0.012 m (0.013 pixel) in northing in band-to-band registration; an exterior geometric performance using ground control points of 8.0-m circular error (95-percent certainty) for orthorectified products and 10.7- to 17.4-m circular error (95-percent certainty) for nonorthorectified products, depending on the geolocation metadata used; and a spatial performance in the range of 1.70 to 2.43 pixels for full width at half maximum, with a modulation transfer function at a Nyquist frequency in the range of 0.032 to 0.084.

Open-File Report

Phenological classification of the United States: A geographic framework for extending multi-sensor time-series data

This study introduces a new geographic framework, phenological classification, for the conterminous United States based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) time-series data and a digital elevation model. The resulting pheno-class map is comprised of 40 pheno-classes, each having unique phenological and topographic characteristics. Cross-comparison of the pheno-classes with the 2001 National Land Cover Database indicates that the new map contains additional phenological and climate information. The pheno-class framework may be a suitable basis for the development of an Advanced Very High Resolution Radiometer (AVHRR)-MODIS NDVI translation algorithm and for various biogeographic studies.

Remote Sensing

Use of ASTER and MODIS thermal infrared data to quantify heat flow and hydrothermal change at Yellowstone National Park

The overarching aim of this study was to use satellite thermal infrared (TIR) remote sensing to monitor geothermal activity within the Yellowstone geothermal area to meet the missions of both the U.S. Geological Survey and the Yellowstone National Park Geology Program. Specific goals were to: 1) address the challenges of monitoring the surface thermal characteristics of the > 10,000 spatially and temporally dynamic thermal features in the Park (including hot springs, pools, geysers, fumaroles, and mud pots) that are spread out over ~ 5000 km 2 , by using satellite TIR remote sensing tools (e.g., ASTER and MODIS), 2) to estimate the radiant geothermal heat flux (GHF) for Yellowstone's thermal areas, and 3) to identify normal, background thermal changes so that significant, abnormal changes can be recognized, should they ever occur (e.g., changes related to tectonic, hydrothermal, impending volcanic processes, or human activities, such as nearby geothermal development). ASTER TIR data (90-m pixels) were used to estimate the radiant GHF from all of Yellowstone's thermal features and update maps of thermal areas. MODIS TIR data (1-km pixels) were used to record background thermal radiance variations from March 2000 through December 2010 and establish thermal change detection limits. A lower limit for the radiant GHF estimated from ASTER TIR temperature data was established at ~ 2.0 GW, which is ~ 30–45% of the heat flux estimated through geochemical thermometry. Also, about 5 km 2 of thermal areas was added to the geodatabase of mapped thermal areas. A decade-long time-series of MODIS TIR radiance data was dominated by seasonal cycles. A background subtraction technique was used in an attempt to isolate variations due to geothermal changes. Several statistically significant perturbations were noted in the time-series from Norris Geyser Basin, however many of these did not correspond to documented thermal disturbances. This study provides concrete examples of the strengths and limitations of current satellite TIR monitoring of geothermal areas, highlighting some specific areas that can be improved. This work provides a framework for future satellite-based thermal monitoring at Yellowstone and other volcanic and geothermal systems.

Idaho, Montana, Oregon, Wyoming