USGS ScienceSearch

USGS · 70192089

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product

Abstract

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Zhuosen Wang, Crystal B. Schaaf, Qingson Sun, JiHyun Kim, Angela M. Erb, Feng Gao, Miguel O. Roman, Yun Yang, Shelley Petroy, Jeffrey Taylor, Jeffrey G. Masek, Jeffrey T. Morisette, Xiaoyang Zhang, Shirley A. Papuga. 2017. Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product. https://doi.org/10.1016/j.jag.2017.03.008

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

KEEP EXPLORING

Related USGS reports

Satellite embeddings for crop type classification: A comparative examination

Embedding datasets encode complex relationships among multiple sources of Earth observation data into a compact format. Here, we evaluated the utility of a 10-m global Satellite Embedding product (SE) for classifying crop types in central California for the year 2020. We compared the classification accuracy of a random forest model based exclusively on the SE layer to an existing random forest model with multiple imagery inputs. Our results showed the SE-based classification had higher agreement with the reference dataset (California Department of Water Resources crop map) than the classification based on Landsat and National Agricultural Imagery Program inputs (94.7% versus. 91.9% overall accuracy, respectively). The performance of individual crop types was consistent across models, ranging from high agreement for rice (98.4% versus 98% accuracy) to lower agreement for pasture, grain, and fallow/young perennial classes (< 65% accuracy in both models). The SE-based workflow used three times less cloud-based computational resources and represented substantial savings of predictor development time. Geospatial embedding products can aid classification efforts by reducing predictor development time and processing demands while maintaining classification accuracy.

California

The impacts of cover crop biomass on satellite-based detectability of cover crops in Maryland

Cover crop adoption in the U.S. has increased over the past decades, increasing the need to quantify their performance and environmental benefits. While remote sensing (RS)-based approaches for detecting cover crop presence have been developed, there has been limited research on how cover crops with varied biomass and management practices influence detectability. Using unique field-level cover crop presence and biomass datasets in Maryland, U.S., we investigated how RS-based detectability changes for cover crops with varied aboveground biomass, planting, and termination dates. Specifically, we proposed a time-integrated satellite-based greenness feature from Harmonized Landsat-8 and Sentinel-2 (HLS) time series from 2017 to 2021 to estimate biomass of cover crops and evaluated their detectability using a phenology-based cover crop detection framework. The impacts of cover crop planting and termination dates on cover crop biomass were also analyzed. Our results demonstrated that Normalized Difference Vegetation Index (NDVI) estimated biomass with higher accuracy compared to other vegetation indices, and the time-integrated model estimated biomass with higher accuracy than the single-date “snapshot” linear model (R2 from 0.53 to 0.66 and RMSE from 922 kg/ha to 747 kg/ha). While the snapshot models were species sensitive, the time-integrated models showed strong robustness across different species. Detectability increased with cover crop biomass, as detected cover crops averaged 963.3 ±719.5 kg/ha compared to 297.2 ± 209.0 kg/ha for non-detected cover crops. Detection accuracy reached 96.1% for fields exceeding 500 kg/ha, compared with an overall accuracy of 62.7%. Earlier planting and later termination increased biomass and detectability, with biomass rising by 4.14 kg/ha/day (p < 0.01). This study demonstrates how management practices affect cover crop biomass and detectability via satellite time series and provides insights that can inform management of cover crops and monitoring of their effects on agroecosystems.

Maryland

Aligning legacy NLCD land cover maps based on Landsat Collection 1 to Collection 2

The transition from Landsat Collection 1 to Collection 2 introduced significant improvements in radiometric and geometric accuracy. However, the improvements cause location misalignment between the existing Landsat-derived land cover products and the new collection. The legacy National Land Cover Database (NLCD) has been used as a cornerstone land cover source for a variety of research. Therefore, a method aligning the legacy NLCD product to Collection 2 is required to ensure its continuity and consistency of service. We developed a strategy to not only align legacy NLCD to match new Collection 2 geometric locations but also improve land cover labeling in the region that was affected by the geometric shifts. The method identifies boundary pixels of homogeneous land cover patches as potential problem areas that are likely impacted by geometric shifts and generates candidate labels from 3 × 3 window with the target pixel at the center and segmentation-derived majority label. Standard phenology patterns of each candidate land cover type are established based on the random samples except boundary pixels within a 1000-pixels × 1000-pixels processing window region. The phenological distance to each standard land cover type pattern is calculated through a penalty dynamic time warping (DTW) method for each target pixel in the boundary region. Finally, the method determines the most suitable label based on the phenological distance from the candidate labels. Both visual and accuracy assessment results demonstrate that the alignment preserves the overall land cover patterns in the original legacy NLCD product while reducing the spatial discrepancies between the Landsat Collection 2 and land cover. In addition, it enhances the accuracy of land cover labeling of boundary pixels. The overall accuracy (OA) was increased by 7% in the land cover boundary regions after alignment. The quality and confusion matrix comparison between the alignment results and the original legacy NLCD confirm the reliability of the method. Our alignment method has the potential to serve as a framework for aligning other Landsat-derived land cover products to future collections.

conterminous United States