Geology topics

Qiang Zhou

Publications and source records attributed to Qiang Zhou.

17 recordsLinked to original sources

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↗

USGS2025

Towards seamless global 30-meter terrestrial monitoring: Evaluating 2022 cloud free coverage of harmonized Landsat and Sentinel-2 (HLS) V2.0

Global observations at 30-m ground sampling distance (GSD) are now possible at a cadence of 1-3 days by combining Landsat 8 and 9 with Sentinel-2A and -2B satellites. Previous studies characterizing pixel-level Landsat-class measurement frequency used data from different sources but offered little information on observation availability after rigorous quality screening. This study examined the coverage frequency of HLS V2.0 data for 2022, the first year all four satellites data were available. These data have had quality control filtering and harmonization, and therefore reflect the spatial-temporal distribution of usable observations. On average, HLS data provide observations every 1.6 days at the global scale, and 2.2 days in the data-scarce tropical regions, regardless of cloud cover. The global mean and median cloud-free observations were 69 and 64, respectively. The frequency of good-quality observations varies geographically and seasonally due to changes in satellite swath overlap, cloud frequency, and solar illumination. High latitudes (>~75°N) exhibit the highest number of cloud-free observations between March and September. However, data are unavailable during winter months due to low solar elevation angles and boreal regions have a lower number of clear observations in the summer months. The tropical regions have the lowest number of clear observations. More frequent HLS observations could improve terrestrial monitoring. We mapped the monthly and weekly number of clear observations globally to show where HLS data could support monthly or sub-weekly time series applications.

IEEE Geoscience and Remote Sensing Letters↗

USGS2025

Land change, fire, and climate weaken carbon sink in the conterminous U.S.

The land carbon sink of the conterminous United States was evaluated using a bottom-up modeling framework and 30-meter land change data from 1985 to 2020. This cross-scale, cross-landscape, and cross-system approach tracked fractional land cover changes and applied regional model calibration. Results show average terrestrial and aquatic carbon sinks of +110 ± 37 and +19 ± 0.5 teragrams of carbon per year, respectively. The terrestrial carbon sink, showing no clear trend, peaked in the 1990s, with more years as a carbon source since 2000, contradicting recent national and global studies. Land change had the largest impact (−70 ± 5.5 teragrams of carbon per year), exceeding impacts of climate (−33 ± 48 teragrams of carbon per year), wildfire (−7.7 ± 2.4 teragrams of carbon per year), and erosion transport (−1.9 ± 0.13 teragrams of carbon per year). The positive CO2 fertilization effect (+69 ± 12 teragrams of carbon per year) was insufficient to maintain the carbon sink strength. Our framework reveals key paths of carbon loss, with implications for carbon budget and energy policies in the United States and beyond.

conterminous United States↗

USGS2022

Tree regrowth duration map from LCMAP collection 1.0 land cover products in the conterminous United States, 1985–2017

Forest covers about one-third of the land area of the conterminous United States (CONUS) and plays an important role in offsetting carbon emissions and supporting local economies. Growing interest in forests as relatively cost-effective nature-based climate solutions, particularly restoration and reforestation activities, has increased the demand for information on forest regrowth and recovery following natural and anthropogenic disturbances (e.g., fire, harvest, or thinning). However, a wall-to-wall mapping of the CONUS tree regrowth duration at an annual time interval and 30-m resolution is still challenging. In this study, we utilized the annual land cover products to develop a dataset to quantify forest regrowth duration for CONUS over 1985–2017. The land cover data used to derive the tree regrowth duration map is from the primary land cover product in the U.S. Geological Survey’s Land Change Monitoring, Assessment, and Projection (LCMAP) collection. The LCMAP product used all available Landsat images to detect disturbances over forest and classify Grass/Shrub to Tree Cover transitions on an annual basis. The average regrowth duration was then calculated for each pixel. The regrowth duration map was validated using human interpreted annual reference data that were collected independently. The validation results show one-year of underestimation and 6-year standard deviation of error between the reference data and regrowth duration map. In southeastern CONUS, where major tree regrowth activities have been observed, our map showed higher accuracy with less than one-year bias and 3.6 years standard deviation of error. Forest in the southeast took around 5 years to recover, which was faster than other regions of CONUS. Many pixels had multiple disturbances during the 33-year study period in the region. The spatial pattern of the tree regrowth indicated intense harvesting activities in this region. The Pacific Northwest coast region was the second main area of tree regrowth, but this region often took multiple decades to recover. Given increasing interest in forests as nature-based climate solutions, the tree regrowth duration map can be used to assess reforestation activities as well as forest recovery following natural disturbance and harvesting.

Conterminous United States↗

USGS2022

Implementation of the CCDC algorithm to produce the LCMAP Collection 1.0 annual land surface change product

The increasing availability of high-quality remote sensing data and advanced technologies have spurred land cover mapping to characterize land change from local to global scales. However, most land change datasets either span multiple decades at a local scale or cover limited time over a larger geographic extent. Here, we present a new land cover and land surface change dataset created by the Land Change Monitoring, Assessment, and Projection (LCMAP) program over the conterminous United States (CONUS). The LCMAP land cover change dataset consists of annual land cover and land cover change products over the period 1985-2017 at 30-meter resolution using Landsat and other ancillary data via the Continuous Change Detection and Classification (CCDC) algorithm. In this paper, we describe our novel approach to implement the CCDC algorithm to produce the LCMAP product suite composed of five land cover and five land surface change related products. The LCMAP land cover products were validated using a collection of ~ 25,000 reference samples collected independently across CONUS. The overall agreement for all years of the LCMAP primary land cover product reached 82.5%. The LCMAP products are produced through the LCMAP Information Warehouse and Data Store (IW+DS) and Shared Mesos Cluster systems that can process, store, and deliver all datasets for public access. To our knowledge, this is the first set of published 30m annual land cover and land cover change datasets that span from the 1980s to the present for the United States. The LCMAP product suite provides useful information for land resource management and facilitates studies to improve the understanding of terrestrial ecosystems and the complex dynamics of the Earth system. The LCMAP system could be implemented to produce global land change products in the future.

Earth System Science Data↗

USGS2022

A novel regression method for harmonic analysis of time series

Harmonic analysis of time series is an important technique in remote sensing to reveal seasonal land surface dynamics. However, frequency selection in the harmonic analysis is often difficult because high-frequency components are useful for delineating seasonal dynamics but sensitive to noise and gaps in time series. On the other hand, it is challenging to obtain temporally continuous satellite data with high quality because of atmospheric contamination. We developed a novel regression method named Harmonic Adaptive Penalty Operator (HAPO) for harmonic analysis of unevenly distributed time series. We introduced a new penalty function to minimize unexpected fluctuations in the model, which can substantially reduce the overfitting issue of regression in time series with temporal gaps. Specifically, the new penalty function minimizes the length of the model curve and the value range difference between the model and the time series observations. We compared HAPO with three widely used regression methods (OLS: Ordinary Least Squares; LASSO: Least Absolute Shrinkage and Selection Operator; and Ridge) in different scenarios using Landsat time series data across the United States. First, we evaluated methods using the Landsat surface reflectance time series within a single year. HAPO showed low and consistent monthly Root Mean Square Deviation (RMSD) values, in which most of the time RMSD of predicted reflectance were less than 0.04. More importantly, HAPO showed consistent and less bias given varying density and irregularity of time series. Second, we evaluated methods using multi-year time series. HAPO was a better predictor of relatively short time series (< 4 years) with steady low RMSD values. When a longer time series ( 4 years) was used, all four methods showed similar RMSD values, but HAPO outperformed the other methods if there were temporal gaps. Therefore, for places with large seasonal observation gaps or for time series that are relatively short (less than 4 years), HAPO can provide more consistent and accurate results in harmonic analysis of time series.

ISPRS Journal of Photogrammetry and Remote Sensing↗

USGS2022

Development of the LCMAP annual land cover product across Hawai'i

Following the completion of land cover and change (LCC) products for the conterminous United States (CONUS), the U.S. Geological Survey's (USGS’s) Land Change Monitoring, Assessment, and Projection initiative has broadened the capability of characterizing continuous historical land change across the full Landsat records for Hawaiʻi at 30-meter resolution. One of the challenges of implementing the LCMAP framework to process annual land cover maps in Hawaiʻi is to collect sufficient high-quality training data. Although multiple datasets depicting land cover information are available in Hawaiʻi, they covered limited time frames and were produced from various remote sensing sources with different, classification categories, spatial resolution, and mapping accuracies. No solo product is suitable to provide LCMAP training data labels on its own. In this paper, we focused on enhancing the LCMAP training datasets to generate land cover products from 2000 to 2019 in Hawaiʻi. A total of 200 independent reference data plots were generated and manually interpreted for validating the mapping results produced by the training datasets. The results revealed that using the appropriate filter of multiple products as training data pools improved the classification model performance. The effect of training datasets (e.g., spatial coverage, quality) on accuracies for different land cover types were summarized. The LCMAP land surface change products for Hawaiʻi are available at https://doi.org/10.5066/P91E8M23 .

Hawaii↗

USGS2022

Monitoring and characterizing multi-decadal variations of urban thermal condition using time-series thermal remote sensing and dynamic land cover data

Urban development and associated land cover and land use change alter the thermal, hydrological, and physical properties of the land surface. Assessments of surface urban heat island (UHI) usually focused on using remote sensing and land cover data to quantify UHI intensity and spatial distribution within a certain period. However, the mechanisms and complex interactions in landscape dynamics and land surface thermal features are still being assessed. In this study, we developed and implemented a novel approach to characterize landscape thermal conditions by focusing on UHI intensity and its spatiotemporal variation using the recently available time series of Landsat land surface temperature and land cover change products. We analyzed land surface temperature changes in urban and surrounding non-urban lands to quantify the UHI intensity and landscape thermal conditions in the Atlanta and Minneapolis metropolitan areas of the United States. Our results revealed that UHI intensities had averages of 3.4 °C and 3.3 °C in the Atlanta and Minneapolis metropolitan areas, respectively. The dominant land cover type in rural areas and urban imperviousness cover determines the UHI intensity. Increasing trends of 0.04 °C/year and 0.01 °C/year in UHI intensity between 1985 and 2018 were found in Atlanta and Minneapolis, respectively. The UHI intensity variations in 1985 and 2018 suggest that the magnitudes and temporal variations of UHI intensity averaged from all urban land cover classes are close to the UHI intensity estimated from the low intensity urban area only while the UHI intensities are more than 2 °C larger in medium to high and high intensity urban areas. The UHI intensities estimated from the maximum temperature that have statistically significant increasing trends suggest that the maximum temperature is a good element for measuring UHI effect. Urban land cover dynamics play an important role in controlling temporal variation of UHI and the UHI hotspots. Our findings support the scientific value of implementing the prototype approach as an objective framework to quantify and monitor UHI intensity at a large geographic extent.

Remote Sensing of Environment↗

USGS2021

Urban heat island and its regional impacts using remotely sensed thermal data – A review of recent developments and methodology

Many novel research algorithms have been developed to analyze urban heat island (UHI) and UHI regional impacts (UHIRIP) with remotely sensed thermal data tables. We present a comprehensive review of some important aspects of UHI and UHIRIP studies that use remotely sensed thermal data, including concepts, datasets, methodologies, and applications. We focus on reviewing progress on multi-sensor image selection, preprocessing, computing, gap filling, image fusion, deep learning, and developing new metrics. This literature review shows that new satellite sensors and valuable methods have been developed for calculating land surface temperature (LST) and UHI intensity, and for assessing UHIRIP. Additionally, some of the limitations of using remotely sensed data to analyze the LST, UHI, and UHI intensity are discussed. Finally, we review a variety of applications in UHI and UHIRIP analyses. The assimilation of time-series remotely sensed data with the application of data fusion, gap filling models, and deep learning using the Google Cloud platform and Google Earth Engine platform also has the potential to improve the estimation accuracy of change patterns of UHI and UHIRIP over long time periods.

Land↗

USGS2021

The effects of urban land cover dynamics on urban heat Island intensity and temporal trends

Assessments of surface urban heat island (UHI) have focused on using remote sensing and land cover data to quantify UHI intensity and spatial distribution within a certain time period by including land cover information. In this study, we implemented a prototype approach to characterize the spatiotemporal variations of UHI using time series of Landsat land surface temperature products and annual land change information. We analyzed UHI distribution and change in Sioux Falls, South Dakota, in the north-central United States and found that the mean UHI intensity in the region was as large as 2.2°C during the period 1986–2017 with an increasing trend of 0.02°C per year within the area with a 5-km non-urban extent. The UHI intensity associated with high intensity urban land cover usually is stronger than with low intensity urban land cover. We evaluated the impact of different non-urban reference extents on UHI variation using different non-urban buffers. The result also suggests that the overall temporal trends of UHI intensity are almost the same when using a 5-km or 10-km non-urban buffer surrounding the urban core. The prototype approach provides a framework to consistently quantify UHI and monitor its change to a large geographic extent.

South Dakota↗

USGS2021

Modeling watershed carbon dynamics as affected by land cover change and soil erosion

Process-based ecosystem carbon cycle models typically incorporate vegetation growth, vegetation mortality, and soil respiration as well as the biotic and environmental drivers that influence these variables. However, few spatially explicit process models can efficiently incorporate the influence of land cover change and carbon lateral movement at regional scales or high spatial resolution. This study uses the Land Use and Carbon Scenario Simulator (LUCAS) to demonstrate the development of a fast ecosystem model that not only considers the basic carbon cycle but also incorporates the impact of land cover change, soil erosion, and soil deposition. As input to the LUCAS modeling framework, we used the integrated biosphere simulator (IBIS) to simulate a non-spatial reference carbon cycling scenario without considering land cover change for the Nisqually River watershed in the northwestern United States. We then used the Land Change Monitoring, Assessment, and Projection (LCMAP) remotely sensed 30-m sequential land cover data to generate annual land change history for the Nisqually River area from 1985 to 2017 and used the Unit Stream Powered Erosion and Deposition model (USPED) to estimate annual soil carbon lateral movement. Finally, we combined the annual carbon outputs from IBIS, the land change history from LCMAP, and the soil erosion and deposition from USPED within the LUCAS simulation framework. Results showed that from 1985 to 2017, along with the dynamic land cover changes, total ecosystem biomass carbon increased from 11.4 to 18.6 TgC, mainly due to forest growth. Total ecosystem soil carbon declined from 31.7 to 29.7 TgC, but the overall loss in soil carbon was not uniform across land cover types. Forestland (forest sector) and grassland lost carbon, while wetland, developed land and agricultural land gained carbon. Forest, grassland, and developed land lost 0.0553 TgC during the study period (1.73 Gg of C per year; 1 Gg = 0.001 Tg) from erosion, while wetland gained 0.0071 TgC (0.22 Gg C per year) from deposition. Agricultural land was neutral in terms of soil erosion.

Washington↗

USGS2021

A novel automatic phenology learning (APL) method of training sample selection using multiple datasets for time-series land cover mapping

The long record of Landsat imagery, which is the cornerstone of Earth observation, provides an opportunity to monitor land use and land cover (LULC) change and understand the interactions between the climate and earth system through time. A few change detection algorithms such as Continuous Change Detection and Classification (CCDC) have been developed to utilize all available Landsat images for change detection and characterization at local or global scales. However, the reliable, rapid, and reproducible collection of training samples have become a challenge for time series land cover classification at a large scale. To meet the challenge, we proposed an automatic phenology learning (APL) method with the assumption that the temporal profiles of samples within the same land cover type are the same or similar at a local scale to generate evenly distributed training samples automatically. We designed the method to build land cover patterns for each category based on consensus samples derived from multiple existing scientific datasets including LANDFIRE's (LF) Existing Vegetation Type (EVT), USGS National Land Cover Database (NLCD), National Agricultural Statistics Service (NASS) Cropland Data Layer (CDL), and National Wetlands Inventory (NWI). Then we calculated the Time-Weighted Dynamic Time Warping (twDTW) distance between any undefined samples and land cover patterns in the same geographical region as prior knowledge. Finally, we selected the optimal land cover category for each undefined sample from the land cover products based on the designed criteria iteratively using the twDTW distance as an indicator. The method was applied in the footprint of 10 selected Landsat Analysis Ready Data (ARD) tiles in the eastern and western conterminous United States (CONUS) to produce annual land cover maps from 1985 to 2017. The accuracy assessment and visual comparison revealed that the APL method can generate reliable training samples without any manual interpretation, producing better land cover results especially for the grass/shrub and wetland land cover classes. Applying the APL method, the overall accuracy of the annual land cover maps was improved by 2% over the accuracy of Land Change Monitoring, Assessment, and Projection (LCMAP) Collection 1.0 Science Products in the research regions. Our results also indicate that the APL method provides an approach for best use of different land cover products and meets the requirement of intensive sampling for training data collection.

Remote Sensing of Environment↗

USGS2020

Methods for rapid quality assessment for national-scale land surface change monitoring

Providing rapid access to land surface change data and information is a goal of the U.S. Geological Survey. Through the Land Change Monitoring, Assessment, and Projection (LCMAP) initiative, we have initiated a monitoring capability that involves generating a suite of ten annual land cover and land surface change datasets across the United States at a 30-m spatial resolution. During the LCMAP automated production on a tile-by-tile basis, erroneous data can occasionally be generated due to hardware or software failure. While crucial to assure the quality of the data, rapid evaluation of results at the pixel level during production is a substantial challenge because of the massive data volumes. Traditionally, product quality relies on the validation after production, which is inefficient to reproduce the whole product when an error occurs. This paper presents a method for automatically evaluating LCMAP results during the production phase based on fourteen indices to quickly find and flag erroneous tiles in the LCMAP products. The methods involved two types of comparisons: comparing LCMAP values across the temporal record to measure internal consistency and calculating agreement with multiple intervals of the National Land Cover Database (NLCD) data to measure the consistency with existing products. We developed indices on a tile-by-tile basis in order to quickly find and flag potential erroneous tiles by comparing with surrounding tiles using local outlier factor analysis. The analysis integrates all indices into a local outlier score (LOS) to detect erroneous tiles distinct from neighbor tiles. Our analysis showed that the methods were sensitive to partially erroneous tiles in the simulated data with a LOS higher than 2. The rapid quality assessment methods also successfully identified erroneous tiles during the LCMAP production, in which land surface change results were not properly saved to the products. The LOS map and indices for rapid quality assessment also point to directions for further investigations. A map of all LOS values by tile for the published LCMAP shows all LOS values are below 2. We also investigated tiles with high LOS to ensure the distinction with neighboring tiles was reasonable. An index in this study shows the overall agreement between LCMAP and NLCD on a tile basis is above 71.5% and has an average at 89.1% across the 422 tiles in the conterminous U.S. The workflow is suitable for other studies with a large volume of image products.

Remote Sensing↗

USGS2020

Gap fill of Land surface temperature and reflectance products in Analysis Ready Data

The recently released Landsat Analysis Ready Data (ARD) over the United States provides the opportunity to investigate landscape dynamics using dense time series observations at 30-m resolution. However, the dataset often contains data gaps (or missing data) because of cloud contamination or data acquisition strategy. We present a new algorithm that focuses on data gap filling using clear observations from orbit overlap regions. Multiple linear regression models were established for each pixel time series to estimate stable predictions and uncertainties. The model's training data came from stratified random samples based on the time series similarity between the pixel and data from the overlap regions. The algorithm was evaluated using four tiles (5,000 × 5,000 30-m pixels for each tile) from 2018 land surface temperature data (LST) in Atlanta, Georgia. The accuracy was assessed using 1,000 randomly masked pixels and daily air temperature from eight ground stations. Both assessments showed the r2 value above 0.75, except two stations with mixed Landsat pixels. We also compared our results with the eMODIS LST product in terms of annual mean temperature. The two maps showed a similar spatial pattern at the region level, but our results showed more spatial detail in the urban area that matched the pattern of impervious surface. We also applied the method on ARD surface reflectance bands at Fairbanks, Alaska, to illustrate its improvements in surface reflectance products and in land change modeling. This approach can also be applied to other datasets, vegetation indexes, or spectral reflectance bands of other sensors.

Remote Sensing↗

USGS2020

Training data selection for annual land cover classification for the LCMAP initiative

The U.S. Geological Survey’s Land Change Monitoring, Assessment, and Projection (LCMAP) initiative characterizes changes in land cover, use, and condition with the goal of producing land change information that improves understanding of the earth system and provides insight into the impacts of land change on society. For LCMAP, all available high-quality data from the Landsat archive is used in a time series approach to detect land surface change. Annual thematic land cover maps are produced by classifying time series models. In this paper, we describe optimization of the classification method used to derive the thematic land cover product. We investigated the influences of auxiliary data, sample size, and training from different sources such as the U.S. Geological Survey’s Land Cover Trends project and National Land Cover Database (NLCD 2001 and NLCD 2011). Results were evaluated and validated based on independent data from the training dataset. We found that refining auxiliary data effectively reduced artifacts in the thematic land cover map that are related to data availability (i.e., SLC-off). The classification accuracy and stability were improved considerably by using a total of 20 million training pixels with a minimum of 600,000 and a maximum of 8 million training pixels per class. Finally, the NLCD 2001 training data delivered the best classification accuracy. Comparing to the original LCMAP classification strategy (Trends training data, 20,000 samples), the optimized classification strategy considerably improved the annual land cover map accuracy.

Remote Sensing↗

USGS2020

Lessons learned implementing an operational continuous United States national land change monitoring capability: The Land Change Monitoring, Assessment, and Projection (LCMAP) approach

Growing demands for temporally specific information on land surface change are fueling a new generation of maps and statistics that can contribute to understanding geographic and temporal patterns of change across large regions, provide input into a wide range of environmental modeling studies, clarify the drivers of change, and provide more timely information for land managers. To meet these needs, the U.S. Geological Survey has implemented a capability to monitor land surface change called the Land Change Monitoring, Assessment, and Projection (LCMAP) initiative. This paper describes the methodological foundations and lessons learned during development and testing of the LCMAP approach. Testing and evaluation of a suite of 10 annual land cover and land surface change data sets over six diverse study areas across the United States revealed good agreement with other published maps (overall agreement ranged from 73% to 87%) as well as several challenges that needed to be addressed to meet the goals of robust, repeatable, and geographically consistent monitoring results from the Continuous Change Detection and Classification (CCDC) algorithm. First, the high spatial and temporal variability of observational frequency led to differences in the number of changes identified, so CCDC was modified such that change detection is dependent on observational frequency. Second, the CCDC classification methodology was modified to improve its ability to characterize gradual land surface changes. Third, modifications were made to the classification element of CCDC to improve the representativeness of training data, which necessitated replacing the random forest algorithm with a boosted decision tree. Following these modifications, assessment of prototype Version 1 LCMAP results showed improvements in overall agreement (ranging from 85% to 90%).

Remote Sensing of Environment↗

USGS2019

Monitoring landscape dynamics in central U.S. grasslands with harmonized Landsat-8 and Sentinel-2 time series data

Remotely monitoring changes in central U.S. grasslands is challenging because these landscapes tend to respond quickly to disturbances and changes in weather. Such dynamic responses influence nutrient cycling, greenhouse gas contributions, habitat availability for wildlife, and other ecosystem processes and services. Traditionally, coarse-resolution satellite data acquired at daily intervals have been used for monitoring. Recently, the harmonized Landsat-8 and Sentinel-2 (HLS) data increased the temporal frequency of the data. Here we investigated if the increased data frequency provided adequate observations to characterize highly dynamic grassland processes. We evaluated HLS data available for 2016 to (1) determine if data from Sentinel-2 contributed to an improvement in characterizing landscape processes over Landsat-8 data alone, and (2) quantify how observation frequency impacted results. Specifically, we investigated into estimating annual vegetation phenology, detecting burn scars from fire, and modeling within-season wetland hydroperiod and growth of aquatic vegetation. We observed increased sensitivity to the start of the growing season (SOST) with the HLS data. Our estimates of the grassland SOST compared well with ground estimates collected at a phenological camera site. We used the Continuous Change Detection and Classification (CCDC) algorithm to assess if the HLS data improved our detection of burn scars following grassland fires and found that detection was considerably influenced by the seasonal timing of the fires. The grassland burned in early spring recovered too quickly to be detected as change events by CCDC; instead, the spectral characteristics following these fires were incorporated as part of the ongoing time-series models. In contrast, the spectral effects from late-season fires were detected both by Landsat-8 data and HLS data. For wetland-rich areas, we used a modified version of the CCDC algorithm to track within-season dynamics of water and aquatic vegetation. The addition of Sentinel-2 data provided the potential to build full time series models to better distinguish different wetland types, suggesting that the temporal density of data was sufficient for within-season characterization of wetland dynamics. Although the different data frequency, in both the spatial and temporal dimensions, could cause inconsistent model estimation or sensitivity sometimes; overall, the temporal frequency of the HLS data improved our ability to track within-season grassland dynamics and improved results for areas prone to cloud contamination. The results suggest a greater frequency of observations, such as from harmonizing data across all comparable Landsat and Sentinel sensors, is still needed. For our study areas, at least a 3-day revisit interval during the early growing season (weeks 14–17) is required to provide a >50% probability of obtaining weekly clear observations.

Minnesota↗