Spatiotemporal residual noise in terrestrial earth observation products, often caused by unfavorable atmospheric conditions, impedes their broad applications. Most users prefer to use gap-filled remote sensing products with time series reconstruction (TSR) algorithms. Applying currently available implementations of TSR to large-volume datasets is time-consuming and challenging for non-professional users with limited computation or storage resources. This study introduces a new open-source software package entitled ‘HANTS-GEE’ that implements a well-known and robust TSR algorithm, i.e. Harmonic ANalysis of Time Series (HANTS), on the Google Earth Engine (GEE) platform for scalable reconstruction of terrestrial earth observation data. Reconstruction tasks can be conducted on user-defined spatiotemporal extents when raw datasets are available on GEE. According to site-based and regional-based case evaluation, the new tool can effectively eliminate cloud contamination in the time series of earth observation data. Compared with traditional PC-based HANTS implementation, the HANTS-GEE provides quite consistent reconstruction results for most terrestrial vegetated sites. The HANTS-GEE can provide scalable reconstruction services with accelerated processing speed and reduced internet data transmission volume, promoting algorithm usage by much broader user communities. To our knowledge, the software package is the first tool to support full-stack TSR processing for popular open-access satellite sensors on cloud platforms.

Satellite-based Earth Observation systems archived a variety of vegetation products during the last 50 years, which can reveal regional to global ecosystem dynamics across diverse spatiotemporal scales. The anomaly metrics such as Vegetation Condition Index (VCI) defined by comparing the current vegetation growth condition to historical average status based on long-term EO-based vegetation products were widely used to delineate abnormal vegetation variation exerted by either climatic or anthropogenic factors (e.g., droughts, wildfires). However, currently available long-term vegetation products may differ from each other in terms of sensors (observational platform or spectral bands), bio-physical definitions (e.g., NDVI, EVI, LAI, and VOD), spatiotemporal resolution, as well as the time-spans, which results in inconsistency across these vegetation products. Taking the VCI as an example, this study evaluated the uncertainty of vegetation anomalies detected based on different vegetation products over the middle reach of the Yangtze River by explicitly considering the effect of sensors, biophysical definitions, and time-spans. The preliminary results showed that VCI derived from NDVI products from different sensors (AVHRR vs. MODIS) induced significant inconsistent anomalies over most landscapes. The differences resulting from products with different biophysical definitions (NDVI vs. EVI, LAI, and VOD) are much lower than those from different sensors but still significant over specific areas. As for the time-spans, the 20-year NDVI based VCI presented a considerable reduction in variance over the study area on average compared to VCI calculated based on 5-year NDVI. In summary, caution should be taken when applying EO-based vegetation products for vegetation anomalies mapping, especially for quantitative assessment.