The compilation of global Landsat data-sets and the ever-lowering costs of computing now make it feasible to monitor the Earth's land cover at Landsat resolutions of 30 m. In this article, we ...describe the methods to create global products of forest cover and cover change at Landsat resolutions. Nevertheless, there are many challenges in ensuring the creation of high-quality products. And we propose various ways in which the challenges can be overcome. Among the challenges are the need for atmospheric correction, incorrect calibration coefficients in some of the data-sets, the different phenologies between compilations, the need for terrain correction, the lack of consistent reference data for training and accuracy assessment, and the need for highly automated characterization and change detection. We propose and evaluate the creation and use of surface reflectance products, improved selection of scenes to reduce phenological differences, terrain illumination correction, automated training selection, and the use of information extraction procedures robust to errors in training data along with several other issues. At several stages we use Moderate Resolution Spectroradiometer data and products to assist our analysis. A global working prototype product of forest cover and forest cover change is included.
This study aims at validating the cloud mask produced by the land surface reflectance code (LaSRC) for Landsat 8 data. To detect clouds in optical satellite imagery, LaSRC uses quality assurance (QA) ...layers, which are produced during the atmospheric correction process. The QA layers include a "cloud mask," which is based on the estimation of a residual metric showing the quality of aerosol inversion, and "high aerosol," which shows the impact of aerosols on the derived surface reflectance. Validation is performed using the "L8 Biome" cloud validation dataset, which is produced by the US Geological Survey, and consists of 96 Landsat 8 scenes distributed globally over 12 different biomes. We show that the LaSRC cloud detection algorithm reliably identifies thick clouds with commission and omission errors less than 4%. Large cloud overdetection errors occur for thin clouds, which is due to the subjectivity of defining and extracting thin clouds in the reference dataset. We conclude this paper with recommendations on using the LaSRC QA layers, and give suggestions on reducing subjectivity, when generating cloud validation datasets.
A combination of Landsat 8 and Sentinel-2 offers a high frequency of observations (3–5 days) at moderate spatial resolution (10–30 m), which is essential for crop yield studies. Existing methods ...traditionally apply vegetation indices (VIs) that incorporate surface reflectances (SRs) in two or more spectral bands into a single variable, and rarely address the incorporation of SRs into empirical regression models of crop yield. In this work, we address these issues by normalizing satellite data (both VIs and SRs) derived from NASA’s Harmonized Landsat Sentinel-2 (HLS) product, through a phenological fitting. We apply a quadratic function to fit VIs or SRs against accumulated growing degree days (AGDDs), which affects the rate of crop development. The derived phenological metrics for VIs and SRs, namely peak, area under curve (AUC), and fitting coefficients from a quadratic function, were used to build empirical regression winter wheat models at a regional scale in Ukraine for three years, 2016–2018. The best results were achieved for the model with near infrared (NIR) and red spectral bands and derived AUC, constant, linear, and quadratic coefficients of the quadratic model. The best model yielded a root mean square error (RMSE) of 0.201 t/ha (5.4%) and coefficient of determination R2 = 0.73 on cross-validation.
Thermal infrared remote sensing observations have been widely used to provide useful information on surface energy and water stress for estimating evapotranspiration (ET). However, the revisit time ...of current high spatial resolution (<100 m) thermal infrared remote sensing systems, sixteen days for Landsat for example, can be insufficient to reliably derive ET information for water resources management. We used in situ ET measurements from multiple Ameriflux sites to (1) evaluate different scaling methods that are commonly used to derive daytime ET estimates from time-of-day observations; and (2) quantify the impact of different revisit times on ET estimates at monthly and seasonal time scales. The scaling method based on a constant evaporative ratio between ET and the top-of-atmosphere solar radiation provided slightly better results than methods using the available energy, the surface solar radiation or the potential ET as scaling reference fluxes. On average, revisit time periods of 2, 4, 8 and 16 days resulted in ET uncertainties of 0.37, 0.55, 0.73 and 0.90 mm per day in summer, which represented 13%, 19%, 23% and 31% of the monthly average ET calculated using the one-day revisit dataset. The capability of a system to capture rapid changes in ET was significantly reduced for return periods higher than eight days. The impact of the revisit on ET depended mainly on the land cover type and seasonal climate, and was higher over areas with high ET. We did not observe significant and systematic differences between the impacts of the revisit on monthly ET estimates that are based on morning or afternoon observations. We found that four-day revisit scenarios provided a significant improvement in temporal sampling to monitor surface ET reducing by around 40% the uncertainty of ET products derived from a 16-day revisit system, such as Landsat for instance.
The Atmospheric Correction Inter-comparison eXercise (ACIX) is an international initiative with the aim to analyse the Surface Reflectance (SR) products of various state-of-the-art atmospheric ...correction (AC) processors. The Aerosol Optical Thickness (AOT) and Water Vapour (WV) are also examined in ACIX as additional outputs of AC processing. In this paper, the general ACIX framework is discussed; special mention is made of the motivation to initiate the experiment, the inter-comparison protocol, and the principal results. ACIX is free and open and every developer was welcome to participate. Eventually, 12 participants applied their approaches to various Landsat-8 and Sentinel-2 image datasets acquired over sites around the world. The current results diverge depending on the sensors, products, and sites, indicating their strengths and weaknesses. Indeed, this first implementation of processor inter-comparison was proven to be a good lesson for the developers to learn the advantages and limitations of their approaches. Various algorithm improvements are expected, if not already implemented, and the enhanced performances are yet to be assessed in future ACIX experiments.
The MODIS instrument provides major advances in moderate resolution earth observation. Improved spatial resolution for land observation at 250 and 500 m and improved spectral band placement provide ...new remote sensing opportunities. NASA has invested in the development of improved algorithms for MODIS, which will provide new data sets for global change research. Surface reflectance is one of the key products from MODIS and is used in developing several higher-order land products. The surface reflectance algorithm builds on the heritage of the Advanced Very High Resolution Radiometer (AVHRR) and SeaWiFS algorithms, taking advantage of the new sensing capabilities of MODIS. Atmospheric correction by the removal of water vapor and aerosol effects provides improvements over previous coarse resolution products and the basis for a new time-series, which will extend through to the NPOESS generation imagers. This paper summarizes the first evaluation of the MODIS surface reflectance product accuracy, in comparison with other data products and in the context of the MODIS instrument performance since launch. The MODIS surface reflectance product will provide an important time-series data set for quantifying global environmental change.
•First long-term global evaluation of the surface reflectance record from AVHRR.•Overall record performance close to proposed specification of 0.071ρ + 0.0071.•No significant variations of surface ...reflectance performance between seasons.•Largest uncertainties found over forest classes.•Strong correlations between band errors over barren and sparsely vegetated areas.
The long-term data record (LTDR) from the Advanced Very High-Resolution Radiometer (AVHRR) provides daily surface reflectance with global coverage from the 1980s to present day, making it a unique source of information for the study of land surface properties and their long-term dynamics. Surface reflectance is a critical input for the generation of products such as vegetation indices, albedo, and land cover. Therefore, it is of utmost importance to quantify its uncertainties to better understand how they might propagate into downstream products. Due to the prolonged length of the surface reflectance LTDR and previous unavailability of a well calibrated reference, no comprehensive evaluation of the complete record has been reported so far. Recently, the United States Geological Survey (USGS) began production of surface reflectance datasets from the Landsat 4–8 satellites, which provide a suitable reference against which the LTDR can be compared to. In this study, we evaluate the LTDRV5 between 1984 and 2011 using surface reflectance data from the Landsat-5 Thematic Mapper (TM5) Collection-1 as a reference. Data from TM5 was obtained from over 740,000 globally distributed scenes which gave a representative set of land surface types and atmospheric conditions. Differences due to observation geometry were accounted for using the Vermote-Justice-Breon (VJB) Bidirectional Reflectance Distribution Function (BRDF) normalization method to adjust the AVHRR surface reflectance to TM5 observation conditions; the spectral response differences were minimized using spectral band adjustment factors (SBAFs) derived from the Earth Observing One (EO-1) Hyperion atmospherically corrected hyperspectral spectra. The performance of the AVHRR record is reported in terms of the accuracy, precision, and uncertainty (APU). Results show that the LTDR performance is close or within the combined uncertainty specification of 0.071ρ + 0.0071, where ρ is the estimated reflectance.
Daily daytime Advanced Very High Resolution Radiometer (AVHRR) 4-km global area coverage data have been processed to produce a Normalized Difference Vegetation Index (NDVI) 8-km equal-area dataset ...from July 1981 through December 2004 for all continents except Antarctica. New features of this dataset include bimonthly composites, NOAA-9 descending node data from August 1994 to January 1995, volcanic stratospheric aerosol correction for 1982-1984 and 1991-1993, NDVI normalization using empirical mode decomposition/reconstruction to minimize varying solar zenith angle effects introduced by orbital drift, inclusion of data from NOAA-16 for 2000-2003 and NOAA-17 for 2003-2004, and a similar dynamic range with the MODIS NDVI. Two NDVI compositing intervals have been produced: a bimonthly global dataset and a 10-day Africa-only dataset. Post-processing review corrected the majority of dropped scan lines, navigation errors, data drop outs, edge-of-orbit composite discontinuities, and other artefacts in the composite NDVI data. All data are available from the University of Maryland Global Land Cover Facility (
http://glcf.umiacs.umd.edu/data/gimms/
).
Differences in the relative spectral response functions of sensors lead to data inconsistencies that should be harmonized before multisensor exploitation. In this paper, we use spectral libraries to ...simulate satellite data and build models to correct them. We then explore and compare different models for coarse and medium spatial resolution optical sensors, including moderate resolution imaging spectroradiometer, advanced very high resolution radiometer (AVHRR), visible infrared imaging radiometer suite, multispectral instrument aboard Sentinel-2, and Operational Land Imager aboard Landsat 8. We found that optimal correction of different bands depends on the model used. For the green and near infrared bands, a multilinear land cover dependent regression improves the accuracy by up to 80.9%. For the red band, a novel exponential dependence of the spectral band adjustment factor with the normalized difference vegetation index (NDVI) provides an accuracy improvement of up to 72.8%. The best way to correct the NDVI value is to use the corrected NIR and red bands using these models. We apply the proposed methods to 445 BELMANIP2 sites using AVHRR data from the long-term data record from 1982-2017. High NDVI pixels result in 30-year trends varying up to 0.06 when comparing uncorrected to spectrally adjusted NDVI. Further application of these methods to NASA's Harmonized Landsat and Sentinel 2 product shows that for the red band and NDVI, our proposed method provides improved accuracy (54.6% and 62.5%) over the linear spectral adjustment currently used.