Geology Reports⌕ Search

Geology topics

Gail L. Schmidt

Publications and source records attributed to Gail L. Schmidt.

13 recordsLinked to original sources
USGS2025

Leveraging airborne imaging spectroscopy and multispectral satellite imagery to map glacial sediment plumes in Kachemak Bay, Alaska

Study Region Kachemak Bay is a fjord-type estuary in the northern Gulf of Alaska. Water quality and habitat characteristics are strongly influenced by freshwater and sediment input from multiple glacierized catchments. Study Focus We present a new method combining imaging spectroscopy from an airborne survey with Landsat and Sentinel-2 imagery to map water surface turbidity originating from glacial runoff based on spectral abundance. We compare the spectral characteristics of turbid glacial water to clear water and generate a high resolution reference map of glacial turbidity in Kachemak Bay. This informs the subsequent analysis of a homogenized, Rayleigh corrected time series of Landsat and Sentinel-2 images and seasonal patterns of turbidity. New Hydrological Insights for the Region Our results provide the most comprehensive data set on water surface turbidity in Kachemak Bay to date and improve understanding of spatial and seasonal variability of glacial turbidity in a data sparse region. July and August have the largest plumes with median sizes around 150 km 2 "> 2 , or around a quarter of Kachemak Bay. Plume sizes typically decrease with decreasing glacier runoff in September and October. We show that imaging spectroscopy aids assessments of turbid water in glacial marine catchments across scales. Leveraging high resolution spectral information allows for water color analyses that are customized to local conditions and catchment characteristics as well as scalable to wider regions.

Alaska↗
USGS2025

The Harmonized Landsat and Sentinel-2 version 2.0 surface reflectance dataset

Frequent multispectral observations of sufficient spatial detail from well-calibrated spaceborne sensors are needed for large-scale terrestrial monitoring. To meet this demand, the NASA Harmonized Landsat and Sentinel-2 (HLS) project was initiated in early 2010s to produce comparable 30-m surface reflectance from the US Landsat 8 Operational Land Imager (OLI) and the European Copernicus Sentinel-2A MultiSpectral Instrument (MSI), and currently from two OLI and two MSI sensors, by applying atmospheric correction to top-of-atmosphere (TOA) reflectance, masking out clouds and cloud shadows, normalizing bi-directional reflectance view angle effects, adjusting for sensor bandpass differences with the OLI as the reference, and providing the harmonized data in a common grid. Several versions of HLS dataset have been produced in the last few years. The recent improvements on almost all the harmonization algorithms had prompted a production of a new HLS dataset, tagged Version 2.0, which was completed in the summer of 2023 and for the first time takes on a global coverage (except for Antarctica). The HLS V2.0 data record starts in April 2013, two months after Landsat 8 launch. For 2022, the first whole year two Landsat and two Sentinel-2 satellites were available, HLS provides a global median of 66 cloud-free observations over land, substantially more than from a single sensor. This paper describes the HLS algorithm improvements and assesses the harmonization efficacy by examining how the reflectance difference between contemporaneous Landsat and Sentinel-2 observations was successively reduced by each harmonization step. The assessment was conducted on 545 pairs of globally distributed same-day Landsat/Sentinel-2 images from 2021 to 2022. Compared to the TOA data, the HLS atmospheric correction slightly increased the reflectance relative difference between Landsat and Sentinel-2 for most of the spectral bands, especially for the two blue bands and the green bands. The subsequent bi-directional reflectance view angle effect normalization effectively reduced the between-sensor reflectance difference present in the atmospherically corrected data for all the spectral bands, and notably to a level below the TOA differences for the red, near-infrared (NIR), and the two shortwave infrared (SWIR) bands. The bandpass adjustment only had a modest effect on reducing the between-sensor reflectance difference. In the final HLS products, the same-day reflectance difference between Landsat and Sentinel-2 was below 4.2% for the red, NIR, and the two SWIR bands, all smaller than the difference in the TOA data. However, the between-sensor differences for the two blue and the green bands remain slightly higher than in TOA data, and this reflects the difficulty in accurately correcting for atmospheric effects in the shorter wavelength visible bands. The data consistency evaluation on a suite of commonly used vegetation indices (VI) calculated from the HLS V2.0 reflectance data showed that the between-sensor VI difference is below 4.5% for most of the indices. These results suggest that the harmonization is robust and the HLS V2.0 data are adequate for quantitative terrestrial applications.

Remote Sensing of Environment↗
USGS2023

The 50-year Landsat collection 2 archive

The Landsat global consolidated data archive now exceeds 50 years. In recognition of the need for consistently processed data across the Landsat satellite series, the U.S. Geological Survey (USGS) initiated collection-based processing of the entire archive that was processed as Collection 1 in 2016. In preparation for the data from the now successfully launched Landsat 9, the USGS reprocessed the Landsat archive as Collection 2 in 2020. This paper describes the rationale for, and the contents and advancements provided by Collection 2, and highlights the differences between the Collection 1 and Collection 2 products. Notably, the Collection 2 products have improved geolocation and, for the first time, the USGS provides a global inventory of Level 2 surface reflectance and surface temperature products. Also for the first time, the USGS used a commercial cloud computing architecture to efficiently process the archive and enable direct cloud access of the Landsat products. The paper concludes with discussion of likely improvements expected in Collection 3 in preparation for the Landsat Next mission that is planned for launch in the early 2030s.

Science of Remote Sensing↗
USGS2020

The Landsat Burned Area algorithm and products for the conterminous United States

Complete and accurate burned area map data are needed to document spatial and temporal patterns of fires, to quantify their drivers, and to assess the impacts on human and natural systems. In this study, we developed the Landsat Burned Area (BA) algorithm, an update from the Landsat Burned Area Essential Climate Variable (BAECV) algorithm. Here, we present the BA algorithm and products, changes relative to the BAECV algorithm and products, and updated validation metrics. We also present spatial and temporal patterns of burned area across the conterminous U.S., how burned area varies in relation to the number of operational Landsat sensors, and a comparison with other burned area datasets, including the BAECV, Monitoring Trends in Burn Severity (MTBS), GeoMAC, and Moderate Resolution Imaging Spectroradiometer (MODIS) MCD64A1.006 data. The BA algorithm identifies burned areas in analysis ready data (ARD) time-series of Landsat imagery from 1984 through 2018 using machine learning, thresholding, and image segmentation. Validation with reference data from high-resolution commercial satellite imagery resulted in omission and commission error rates averaging 19% and 41%, respectively. In comparison, validation with Landsat reference data had omission and commission error rates averaging 40% and 28%, respectively when burned areas in cultivated crops and pasture/hay land-cover types were excluded. Both validation tests documented lower commission error rates relative to the BAECV products. The amount of burned area detected varies not only in response to climate but also with the number of operational sensors and scenes collected. The combined amount of burned area detected by multiple sensors was larger than from any individual sensor, but there was no significant difference between individual sensors. Therefore, we used BA products from individual sensors to assess trends over time and all available sensors to compare with other existing BA products. From 1984 through 2018, annual burned area averaged 30,000 km2, ranged between 14,000 km2 in 1991 and 46,500 km2 in 2012, and increased over time at a rate of 356 km2/year. Compared to existing burned area products, the new Landsat BA products identified 29% more burned area than the BAECV products (1984–2015), 183% more than the MTBS/GeoMAC products (1984–2018), and 56% more than the MCD64A1.006 products (2003–2018). The products had similar patterns of year-to-year variability; the R2 values of linear regressions between annual burned area were >0.70 with the BAECV products and the MTBS/GeoMAC products, but somewhat lower for the MCD64A1.006 product (R2 = 0.66). The BA products are routinely produced as new Landsat data are collected and provide a unique data source to monitor and assess the spatial and temporal patterns and the impacts of fire.

Remote Sensing of Environment↗
USGS2020

New operational national satellite burned area product

Introduction Lack of consistent spatial and temporal fire information with relevant spatial resolution hinders land management and broad-scale assessments of fire activity, especially in the eastern United States and the Great Plains where fi re is important ecologically and culturally. Remote sensing can be used to monitor fi re activity, augment existing fi re data, and fill information gaps. In particular, Landsat offers one of the most complete time series of remote sensing data sets as the Landsat satellites with spectral bands useful for mapping fi res and burn severity have been operational since 1984. Furthermore, Landsat satellite imagery collect data at a resolution useful for on-the-ground comparisons and management decisions. Methods A gradient-boosting regression model algorithm was used to predict burn probabilities (BP), indicating the likelihood that a pixel had burned in a fire. Then the algorithm translated the burn probability images to burn classification (BC) images using thresholding and region growing. Burned areas smaller than 5 acres were removed to reduce noise. The BP and BC products were generated for Landsat scenes collected from 1984 to present with <= 80% cloud cover and are available for download from USGSs EarthExplorer website (www.earthexplorer. usgs.gov). Annual summaries were also produced from the scene-level BC and BP products. These include: 1) the maximum burn probability value from all scenes in a year, (2) the burn classification count, or the number of times a pixel was classified as burned in a year, (3) the burn date, or the day of year of the first Landsat scene a pixel was classified as burned, and (4) a burn classification which removes burned areas still visible from previous years with a filter. Results The Landsat Burned Area Products identified 183% more burned area than the Monitoring Trends in Burn Severity and 56% more burned area than the MODIS MCD64A1.006 products2. From 1984 through 2018, annual burned area mapped by the Landsat Burned Area Products averaged 30,000 km2, ranged between 14,000 km2 in 1991 and 46,500 km2 in 2012, and increased over time at a rate of 356 km2/year. The Landsat Burned Area products consistently map burned areas 5 acres or larger over time while capturing unburned islands within fi res. Such characteristics are critical for understanding long-term impacts of fi res on human and natural systems. The products are being routinely produced as new Landsat data are collected and demonstrate the utility remote sensing offers for monitoring landscape change.

Report↗
USGS2017

Cloud detection algorithm comparison and validation for operational Landsat data products

Clouds are a pervasive and unavoidable issue in satellite-borne optical imagery. Accurate, well-documented, and automated cloud detection algorithms are necessary to effectively leverage large collections of remotely sensed data. The Landsat project is uniquely suited for comparative validation of cloud assessment algorithms because the modular architecture of the Landsat ground system allows for quick evaluation of new code, and because Landsat has the most comprehensive manual truth masks of any current satellite data archive. Currently, the Landsat Level-1 Product Generation System (LPGS) uses separate algorithms for determining clouds, cirrus clouds, and snow and/or ice probability on a per-pixel basis. With more bands onboard the Landsat 8 Operational Land Imager (OLI)/Thermal Infrared Sensor (TIRS) satellite, and a greater number of cloud masking algorithms, the U.S. Geological Survey (USGS) is replacing the current cloud masking workflow with a more robust algorithm that is capable of working across multiple Landsat sensors with minimal modification. Because of the inherent error from stray light and intermittent data availability of TIRS, these algorithms need to operate both with and without thermal data. In this study, we created a workflow to evaluate cloud and cloud shadow masking algorithms using cloud validation masks manually derived from both Landsat 7 Enhanced Thematic Mapper Plus (ETM +) and Landsat 8 OLI/TIRS data. We created a new validation dataset consisting of 96 Landsat 8 scenes, representing different biomes and proportions of cloud cover. We evaluated algorithm performance by overall accuracy, omission error, and commission error for both cloud and cloud shadow. We found that CFMask, C code based on the Function of Mask (Fmask) algorithm, and its confidence bands have the best overall accuracy among the many algorithms tested using our validation data. The Artificial Thermal-Automated Cloud Cover Algorithm (AT-ACCA) is the most accurate nonthermal-based algorithm. We give preference to CFMask for operational cloud and cloud shadow detection, as it is derived from a priori knowledge of physical phenomena and is operable without geographic restriction, making it useful for current and future land imaging missions without having to be retrained in a machine-learning environment.

Remote Sensing of Environment↗
USGS2017

Mapping burned areas using dense time-series of Landsat data

Complete and accurate burned area data are needed to document patterns of fires, to quantify relationships between the patterns and drivers of fire occurrence, and to assess the impacts of fires on human and natural systems. Unfortunately, in many areas existing fire occurrence datasets are known to be incomplete. Consequently, the need to systematically collect burned area information has been recognized by the United Nations Framework Convention on Climate Change and the Intergovernmental Panel on Climate Change, which have both called for the production of essential climate variables (ECVs), including information about burned area. In this paper, we present an algorithm that identifies burned areas in dense time-series of Landsat data to produce the Landsat Burned Area Essential Climate Variable (BAECV) products. The algorithm uses gradient boosted regression models to generate burn probability surfaces using band values and spectral indices from individual Landsat scenes, lagged reference conditions, and change metrics between the scene and reference predictors. Burn classifications are generated from the burn probability surfaces using pixel-level thresholding in combination with a region growing process. The algorithm can be applied anywhere Landsat and training data are available. For this study, BAECV products were generated for the conterminous United States from 1984 through 2015. These products consist of pixel-level burn probabilities for each Landsat scene, in addition to, annual composites including: the maximum burn probability and a burn classification. We compared the BAECV burn classification products to the existing Global Fire Emissions Database (GFED; 1997–2015) and Monitoring Trends in Burn Severity (MTBS; 1984–2013) data. We found that the BAECV products mapped 36% more burned area than the GFED and 116% more burned area than MTBS. Differences between the BAECV products and the GFED were especially high in the West and East where the BAECV products mapped 32% and 88% more burned area, respectively. However, the BAECV products found less burned area than the GFED in regions with frequent agricultural fires. Compared to the MTBS data, the BAECV products identified 31% more burned area in the West, 312% more in the Great Plains, and 233% more in the East. Most pixels in the MTBS data were detected by the BAECV, regardless of burn severity. The BAECV products document patterns of fire similar to those in the GFED but also showed patterns of fire that are not well characterized by the existing MTBS data. We anticipate the BAECV products will be useful to studies that seek to understand past patterns of fire occurrence, the drivers that created them, and the impacts fires have on natural and human systems.

continental United States↗
USGS2014

Development of a generic auto-calibration package for regional ecological modeling and application in the Central Plains of the United States

Process-oriented ecological models are frequently used for predicting potential impacts of global changes such as climate and land-cover changes, which can be useful for policy making. It is critical but challenging to automatically derive optimal parameter values at different scales, especially at regional scale, and validate the model performance. In this study, we developed an automatic calibration (auto-calibration) function for a well-established biogeochemical model—the General Ensemble Biogeochemical Modeling System (GEMS)-Erosion Deposition Carbon Model (EDCM)—using data assimilation technique: the Shuffled Complex Evolution algorithm and a model-inversion R package—Flexible Modeling Environment (FME). The new functionality can support multi-parameter and multi-objective auto-calibration of EDCM at the both pixel and regional levels. We also developed a post-processing procedure for GEMS to provide options to save the pixel-based or aggregated county-land cover specific parameter values for subsequent simulations. In our case study, we successfully applied the updated model (EDCM-Auto) for a single crop pixel with a corn–wheat rotation and a large ecological region (Level II)—Central USA Plains. The evaluation results indicate that EDCM-Auto is applicable at multiple scales and is capable to handle land cover changes (e.g., crop rotations). The model also performs well in capturing the spatial pattern of grain yield production for crops and net primary production (NPP) for other ecosystems across the region, which is a good example for implementing calibration and validation of ecological models with readily available survey data (grain yield) and remote sensing data (NPP) at regional and national levels. The developed platform for auto-calibration can be readily expanded to incorporate other model inversion algorithms and potential R packages, and also be applied to other ecological models.

Ecological Informatics↗
USGS2007

A multi-scale segmentation approach to filling gaps in Landsat ETM+ SLC-off images

On 31 May 2003, the Landsat Enhanced Thematic Plus (ETM+) Scan Line Corrector (SLC) failed, causing the scanning pattern to exhibit wedge-shaped scan-to-scan gaps. We developed a method that uses coincident spectral data to fill the image gaps. This method uses a multi-scale segment model, derived from a previous Landsat SLC-on image (image acquired prior to the SLC failure), to guide the spectral interpolation across the gaps in SLC-off images (images acquired after the SLC failure). This paper describes the process used to generate the segment model, provides details of the gap-fill algorithm used in deriving the segment-based gap-fill product, and presents the results of the gap-fill process applied to grassland, cropland, and forest landscapes. Our results indicate this product will be useful for a wide variety of applications, including regional-scale studies, general land cover mapping (e.g. forest, urban, and grass), crop-specific mapping and monitoring, and visual assessments. Applications that need to be cautious when using pixels in the gap areas include any applications that require per-pixel accuracy, such as urban characterization or impervious surface mapping, applications that use texture to characterize landscape features, and applications that require accurate measurements of small or narrow landscape features such as roads, farmsteads, and riparian areas.

International Journal of Remote Sensing↗
USGS2006

The MODIS reprojection tool

The MODIS Reprojection Tool (MRT) is designed to help individuals work with MODIS Level-2G, Level-3, and Level-4 land data products. These products are referenced to a global tiling scheme in which each tile is approximately 10° latitude by 10° longitude and non-overlapping (Fig. 9.1). If desired, the user may reproject only selected portions of the product (spatial or parameter subsetting). The software may also be used to convert MODIS products to file formats (generic binary and GeoTIFF) that are more readily compatible with existing software packages. The MODIS land products distributed by the Land Processes Distributed Active Archive Center (LP DAAC) are in the Hierarchical Data Format - Earth Observing System (HDF-EOS), developed by the National Center for Supercomputing Applications at the University of Illinois at Urbana Champaign for the NASA EOS Program. Each HDF-EOS file is comprised of one or more science data sets (SDSs) corresponding to geophysical or biophysical parameters. Metadata are embedded in the HDF file as well as contained in a .met file that is associated with each HDF-EOS file. The MRT supports 8-bit, 16-bit, and 32-bit integer data (both signed and unsigned), as well as 32-bit float data. The data type of the output is the same as the data type of each corresponding input SDS.

Book chapter↗
USGS2005

Landsat 7 scan line corrector-off gap-filled product development

The Landsat 7 Enhanced Thematic Mapper Plus (ETM+) scan line corrector (SLC) failed on May 31, 2003, causing the scanning pattern to exhibit wedge-shaped scan-to-scan gaps. The ETM+ has continued to acquire data with the SLC powered off, leading to images that are missing approximately 22 percent of the normal scene area. To improve the utility of the SLC-off data, the U.S. Geological Survey (USGS) developed new products that use the data from multiple ETM+ scenes to provide complete ground coverage. These gap-filled products were developed and deployed in two phases. The gaps in the Phase I products are filled with data from imagery collected previously with a functional SLC (SLC-on). A single SLC-on scene provides complete coverage of the scan gaps, making the gap-filling procedure straightforward. Several radiometric adjustment techniques for matching the SLC-on fill scene to the SLC-off primary scene were evaluated for performance, processing speed, and ease of implementation. A simple local histogram matching method was adopted as a result of this evaluation. The Phase II products use data from multiple SLC-off scenes to fill the scan gaps with more recent data. Because the locations of the scan gaps are different for each SLC-off scene, the gap-filling process must account for scan gap interactions. The Phase II product development included a more comprehensive study of candidate radiometric adjustment techniques. This study showed that the histogram matching method used in Phase I, with minor refinements, provided the best overall performance and was adopted for Phase II as well.

Conference Paper↗