USGS ScienceSearch

USGS · 70221549

Detecting subtle change from dense landsat time series: Case studies of mountain pine beetle and spruce beetle disturbance

Abstract

In contrast to abrupt changes caused by land cover conversion, subtle changes driven by a shift in the condition, structure, or other biological attributes of land often lead to minimal and slower alterations of the terrestrial surface. Accurate mapping and monitoring of subtle change are crucial for an early warning of long-term gradual change that may eventually result in land cover conversion. Freely accessible moderate-resolution datasets such as the Landsat archive have great potential to characterize subtle change by capturing low-magnitude spectral changes in long-term observations. However, past studies have reported limited success in accurately extracting subtle changes from satellite-based time series analysis. In this study, we introduce a supervised framework named ‘PIDS’ to detect subtle forest disturbance from a comprehensive Landsat data archive by leveraging disturbance-based calibration sites. PIDS consists of four components: (1) Parameter optimization; (2) Index selection; (3) Dynamic stratified monitoring; and (4) Spatial consideration. PIDS was applied to map the early stage of bark beetle infestations (i.e., a lower per-pixel fraction of trees cover that show visual signs of infestation), which are a typical example of subtle change in conifer forests. Landsat Analysis Ready Data were used as the time series inputs for mapping mountain pine beetle and spruce beetle disturbance between 2001 and 2019 in Colorado, USA. PIDS-detection map assessment showed that the overall performance of PIDS (namely ‘F1 score’) was 0.86 for mountain pine beetle and 0.73 for spruce beetle, making a substantial improvement (> 0.3) compared to other approaches/products including COntinuous monitoring of Land Disturbance, LandTrendr, and the National Land Cover Database forest disturbance product. A sub-pixel analysis of tree canopy mortality percentage was performed by linking classified high-resolution (0.3- and 1-m) aerial imagery and 30-m PIDS-detection maps. Results show that PIDS typically detects mountain pine beetle infestation when ≥56% of a Landsat pixel is occupied by red-stage canopy mortality (one year after initial infestation), and spruce beetle infestation when ≥55% is occupied by gray-stage mortality (two years after initial infestation). This study addresses an important methodological goal pertinent to the utility of event-based reference samples for detecting subtle forest change, which could be potentially applied to other types of subtle land change.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 37.020098201368114° to 40.91351257612758° latitude; -109.0283203125° to -106.14990234375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Su Ye, John Rogan, Zhe Zhu, Todd Hawbaker, Sarah J. Hart, Robert A. Andrus, Arjan J.H. Meddens, Jeffery A. Hicke, J. Ronald Eastman, Dominik Kulakowski. 2021. Detecting subtle change from dense landsat time series: Case studies of mountain pine beetle and spruce beetle disturbance. https://doi.org/10.1016/j.rse.2021.112560

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

On-demand global Landsat evapotranspiration product: Development, evaluation, and dissemination

Global actual evapotranspiration (ET) is one of the essential climate variables needed to understand and manage the relationships among food, energy, and water resources. The U.S. Geological Survey Earth Resources Observation and Science (EROS) Center launched a provisional ET product in 2020, offering on-demand, field-scale global coverage derived from Landsat data through the EROS Science Processing Architecture (ESPA) platform. The ESPA interface provides ET data for cloud-free Landsat overpasses starting in 1982 with Landsat 4 through the current Landsat 9. The ET data are delivered as a Provisional Level-3 Science product created using the Operational Simplified Surface Energy Balance (SSEBop) model. Landsat surface temperature and reference ET are the main model drivers along with vegetation index and net radiation for model parameterization. A large volume of Landsat-based ET orders (e.g., over 1,200,000 images from June 2020 through December 2025) around the world indicate increasing awareness and application of the ET data. The ESPA platform enables land and water resource managers and researchers to access a first-order ET product without requiring advanced knowledge of remote sensing technology or evapotranspiration modeling. We present the methodology and workflow of the on-demand Landsat ET product and its performance evaluations over diverse hydro-climatic settings. The product can help estimate field-scale consumptive water use and thus quickly and consistently assess historical water use, allocation, and budget to inform water management under changing environments. Future ET data aggregated to monthly and seasonal time scales are expected to enhance integration with decision-making tools and procedures.

Remote Sensing of Environment

A framework for integrating spatiotemporal deep learning methods with landsat for annual land cover and impervious surface mapping

Land cover information is essential for understanding Earth’s surface dynamics and how vegetation, water, soil, climate, and terrain interact. The National Land Cover Database (NLCD) has been the authoritative source for consistent U.S. land cover mapping. To extend NLCD’s temporal resolution and reduce production latency, we developed the Land Cover Artificial Mapping System (LCAMS)—a prototype spatiotemporal deep learning framework piloted as the foundation for the new Annual NLCD. LCAMS builds on concepts from legacy NLCD and the U.S. Geological Survey Land Change Monitoring, Assessment, and Projection (LCMAP) initiatives. It employs a loosely coupled two-stage architecture consisting of independent but functionally interdependent spatial and temporal models. Spatial models extract per-year information from Landsat data, while the temporal models refine the spatial outputs to enforce inter-annual consistency—critical for reliable land change monitoring. LCAMS produces annual 30 m resolution land cover and impervious surface outputs, with region-specific fine-tuning to generalize across diverse landscapes and temporal dynamics. Validation was conducted using an independent dataset of 1925 randomly sampled plots from five U.S. Landsat Analysis Ready Data (ARD) tiles spanning 1985-2021, selected for spatial and temporal variability. This dataset was used consistently to evaluate LCAMS, Legacy NLCD, and LCMAP. Using the NLCD legend, LCAMS achieved 72.1 ± 1.60% overall agreement, compared to 71.1 ± 1.7% agreement for Legacy NLCD. Using the LCMAP legend, LCAMS achieved 83.4 ± 1.22% agreement, compared to 84.6 ± 1.11% agreement for LCMAP. Overall, LCAMS delivers comparable accuracy while offering higher thematic resolution, longer temporal coverage, and automated production of annual 30 m CONUS land cover.

Remote Sensing of Environment

Towards global mapping of dynamic surface water extents using Sentinel-1 SAR data

We introduce a fully automated and scalable method for mapping surface water extents from single-acquisition Sentinel-1 synthetic aperture radar (SAR) imagery. This approach integrates adaptive thresholding of radiometric terrain-corrected SAR backscatter data, fuzzy-logic classification, region growing, dark land estimation, and a bimodality test to minimize false positives in low-backscattering areas and false negatives in high-backscattering areas. By combining these steps, the algorithm achieves classification accuracies exceeding 85% in detecting surface water extents across diverse environmental conditions. Accuracy was first assessed at meter scale using 52 PlanetScope scenes acquired worldwide in September–October 2019; the algorithm achieved 93% overall accuracy, 86% user's accuracy, and 94% producer's accuracy. Global robustness was then evaluated by processing every Sentinel-1 acquisition from 1 to 12 November 2023 and cross-comparing the resulting maps with 6561 temporally matched observational products for end-users from remote sensing analysis (OPERA) dynamic surface water extent from Harmonized Landsat and Sentinel-2 (DSWx-HLS) products. This large-scale test yielded 90% user's and 94% producer's accuracies, confirming reliable performance at continental extent. Additional case studies demonstrate the algorithm's ability to handle surface water extent in sand-dominated deserts, to track seasonal amplitude in Folsom Lake (California), drought-induced loss in Cerro Prieto Reservoir (Mexico), and rapid filling of the Grand Ethiopian Renaissance Dam. These results show that the method scales across local to global domains and maintains high accuracy, providing a practical tool for near-real-time monitoring of floods, droughts, and water-resource management. Because the approach is sensor-agnostic, it can be ported to forthcoming L- and S-band missions such as NASA-ISRO synthetic aperture radar (NISAR), broadening its applicability to future hydrologic observations.

Remote Sensing of Environment