Making Landsat time series consistent: Evaluating and improving Landsat analysis ready data

Recently, the United States Geological Survey (USGS) has released a new dataset, called Landsat Analysis Ready Data (ARD), which is designed specifically for facilitating time series analysis. In this study, we evaluated the temporal consistency of this new dataset and recommended several processing streamlines for improving data consistency. Specifically, we examined the impacts of data resampling, cloud/cloud shadow detection, Bidirectional Reflectance Distribution Function (BRDF) correction, and topographic correction on the temporal consistency of the Landsat Time Series (LTS).We have four major observations. First, single-resampled data (ARD) are generally more consistent than double-resampled data (re-projected Collection 1 data), but the difference is very minor. Second, the improved cloud and cloud shadow detection approach (e.g., Fmask 4.0 vs. 3.3) moderately increased data consistency. Third, BRDF correction contributed the most in making LTS consistent. Finally, we corrected the topographic effects by using several widely used algorithms, including Sun-Canopy-Sensor (SCS), a semiempirical SCS (SCS+C), and Illumination Correction (IC) algorithms, however they were found to have very limited or even negative impacts on the consistency of LTS. Therefore, we recommend using Landsat ARD with the improved cloud and cloud shadow detection approach (Fmask 4.0), and with BRDF correction for routine time series analysis.