Assessing the Effects of Time Interpolation of NDVI Composites on Phenology Trend Estimation

Abstract

The accurate evaluation of shifts in vegetation phenology is essential for understanding of vegetation responses to climate change. Remote-sensing vegetation index (VI) products with multi-day scales have been widely used for phenology trend estimation. VI composites should be interpolated into a daily scale for extracting phenological metrics, which may not fully capture daily vegetation growth, and how this process affects phenology trend estimation remains unclear. In this study, we chose 120 sites over four vegetation types in the mid-high latitudes of the northern hemisphere, and then a Moderate Resolution Imaging Spectroradiometer (MODIS) MCD43A4 daily surface reflectance data was used to generate a daily normalized difference vegetation index (NDVI) dataset in addition to an 8-day and a 16-day NDVI composite datasets from 2001 to 2019. Five different time interpolation methods (piecewise logistic function, asymmetric Gaussian function, polynomial curve function, linear interpolation, and spline interpolation) and three phenology extraction methods were applied to extract data from the start of the growing season and the end of the growing season. We compared the trends estimated from daily NDVI data with those from NDVI composites among (1) different interpolation methods; (2) different vegetation types; and (3) different combinations of time interpolation methods and phenology extraction methods. We also analyzed the differences between the trends estimated from the 8-day and 16-day composite datasets. Our results indicated that none of the interpolation methods had significant effects on trend estimation over all sites, but the discrepancies caused by time interpolation could not be ignored. Among vegetation types with apparent seasonal changes such as deciduous broadleaf forest, time interpolation had significant effects on phenology trend estimation but almost had no significant effects among vegetation types with weak seasonal changes such as evergreen needleleaf forests. In addition, trends that were estimated based on the same interpolation method but different extraction methods were not consistent in showing significant (insignificant) differences, implying that the selection of extraction methods also affected trend estimation. Compared with other vegetation types, there were generally fewer discrepancies between trends estimated from the 8-day and 16-day dataset in evergreen needleleaf forest and open shrubland, which indicated that the dataset with a lower temporal resolution (16-day) can be applied. These findings could be conducive for analyzing the uncertainties of monitoring vegetation phenology changes.

A Novel Remote Sensing-Based Modeling Approach for Maize Light Extinction Coefficient Determination

“…A cubic convolution gap-filling approach was utilized to provide temporal interpolated LAI data for the days without measurements during the period studied, similar to the interpolation methods implemented by [92,93]. The interpolation method for the measured LAI was explicitly applied to each data collection time series per measurement station within the frequently irrigated fields at the LIRF and IIC.…”

Assessing multi‐sensor hourly maize evapotranspiration estimation using a one‐source surface energy balance approach

“…The h c measurements were considered from the ground surface to the upper leaves. A cubic convolution gap-filling approach provided temporally interpolated LAI and h c data for the days without measurements, similar to the interpolation methods implemented by Vorobiova and Chernov (2017) and Li et al (2021). The interpolation method for measured LAI and h c was explicitly applied to each data collection time series per measurement station within the frequently irrigated fields at LIRF and IIC.…”

Multi-year monitoring of wheat phenology and effect of climate change in the south Asian region using Sentinel-2 NDVI time series analysis

“…Tillering, jointing and maturity are further detected from the NDVI trend but for a single year. Different interpolation methods 13 including, piece-wise, linear, polynomial, gaussian and spline are compared for 8day and 16day data to detect phenology using maximum change, curve change and dynamic threshold method based on NDVI values, smoothed using SG filter and maximum value composite. The authors highlight that the type interpolation method has no significant effect on results however, using a more frequent data of less than 8days can accurately detect the phenological stages using rate of change methods and threshold methods with threshold set at 20 to 30 %.…”