robustly detect change with different magnitudes (>0.1 NDVI) within time series with different noise levels(0.01 –0.07σ) and seasonal amplitudes (0.1 –0.5 NDVI). Additionally, BFAST was applied to 16-day NDVI