Crop phenology is an important parameter for crop growth monitoring, yield prediction, and growth simulation. The dynamic threshold method is widely used to retrieve vegetation phenology from remotely sensed vegetation index time series. However, crop growth is not only driven by natural conditions, but also modified through field management activities. Complicated planting patterns, such as multiple cropping, makes the vegetation index dynamics less symmetrical. These impacts are not considered in current approaches for crop phenology retrieval based on the dynamic threshold method. Thus, this paper aimed to (1) investigate the optimal thresholds for retrieving the start of the season (SOS) and the end of the season (EOS) of different crops, and (2) compare the performances of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) in retrieving crop phenology with a modified version of the dynamic threshold method. The reference data included SOS and EOS ground observations of three major crop types in 2015 and 2016, which includes rice, wheat, and maize. Results show that (1) the modification of the original method ensures a 100% retrieval rate, which was not guaranteed using the original method. The modified dynamic threshold method is more suitable to retrieve crop SOS/EOS because it considers the asymmetry of crop vegetation index time series. (2) It is inappropriate to retrieve SOS and EOS with the same threshold for all crops, and the commonly used 20% or 50% thresholds are not the optimal thresholds for all crops. (3) For single and late rice, the accuracies of the SOS estimations based on EVI are generally higher compared to those based on NDVI. However, for spring maize and summer maize, results based on NDVI give higher accuracies. In terms of EOS, for early rice and summer maize, estimates based on EVI result in higher accuracies, but, for late rice and winter wheat, results based on NDVI are closer to the ground records.

The leaf area index (LAI) is not only an important parameter used to describe the geometry of vegetation canopy but also a key input variable for ecological models. One of the most commonly used methods for LAI estimation is to establish an empirical relationship between the LAI and the vegetation index (VI). However, the LAI-VI relationships had high seasonal variability, and they differed among phenophases and VIs. In this study, the LAI-VI relationships in different phenophases and for different VIs (i.e., the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI) and near-infrared reflectance of vegetation (NIRv)) were investigated based on 82 site-years of LAI observed data and the Moderate Resolution Imaging Spectroradiometer (MODIS) VI products. Significant LAI-VI relationships were observed during the vegetation growing and declining periods. There were weak LAI-VI relationships (p > 0.05) during the flourishing period. The accuracies for the LAIs estimated with the piecewise LAI-VI relationships based on different phenophases were significantly higher than those estimated based on a single LAI-VI relationship for the entire vegetation active period. The average root mean square error (RMSE) ± standard deviation (SD) value for the LAIs estimated with the piecewise LAI-VI relationships was 0.38 ± 0.13 (based on the NDVI), 0.41 ± 0.13 (based on the EVI) and 0.41 ± 0.14 (based on the NIRv), respectively. In comparison, it was 0.46 ± 0.13 (based on the NDVI), 0.55 ± 0.15 (based on the EVI) and 0.55 ± 0.15 (based on the NIRv) for those estimated with a single LAI-VI relationship. The performance of the three VIs in estimating the LAI also varied among phenophases. During the growing period, the mean RMSE ± SD value for the estimated LAIs was 0.30 ± 0.11 (LAI-NDVI relationships), 0.37 ± 0.11 (LAI-EVI relationships) and 0.36 ± 0.13 (LAI-NIRv relationships), respectively, indicating the NDVI produced significantly better LAI estimations than those from the other two VIs. In contrast, the EVI produced slightly better LAI estimations than those from the other two VIs during the declining period (p > 0.05), and the mean RMSE ± SD value for the estimated LAIs was 0.45 ± 0.16 (LAI-NDVI relationships), 0.43 ± 0.23 (LAI-EVI relationships) and 0.45 ± 0.25 (LAI-NIRv relationships), respectively. Hence, the piecewise LAI-VI relationships based on different phenophases were recommended for the estimations of the LAI instead of a single LAI-VI relationship for the entire vegetation active period. Furthermore, the optimal VI in each phenophase should be selected for the estimations of the LAI according to the characteristics of vegetation growth.

Agricultural land use and cropping patterns are closely related to food production, soil degradation, water resource management, greenhouse gas emission, and regional climate alterations. Methods for reliable and cost-efficient mapping of cropping pattern, as well as their changes over space and time, are therefore urgently needed. To cope with this need, we developed a phenology-based method to map cropping patterns based on time-series of vegetation index data. The proposed method builds on the well-known ‘threshold model’ to retrieve phenological metrics. Values of four phenological parameters are used to identify crop seasons. Using a set of rules, the crop season information is translated into cropping pattern. To illustrate the method, cropping patterns were determined for three consecutive years (2008–2010) in the Henan province of China, where reliable validation data was available. Cropping patterns were derived using eight-day composite MODIS Enhanced Vegetation Index (EVI) data. Results show that the proposed method can achieve a satisfactory overall accuracy (~84%) in extracting cropping patterns. Interestingly, the accuracy obtained with our method based on MODIS EVI data was comparable with that from Landsat-5 TM image classification. We conclude that the proposed method for cropland and cropping pattern identification based on MODIS data offers a simple, yet reliable way to derive important land use information over large areas.

The long-term Normalized Difference Vegetation Index (NDVI) time-series data set generated from the Advanced Very High Resolution Radiometers (AVHRR) has been widely used to monitor vegetation activity change. The third version of NDVI (NDVI3g) produced by the Global Inventory Modeling and Mapping Studies (GIMMS) group was released recently. The comparisons between the new and old versions should be conducted for linking existing studies with future applications of NDVI3g in monitoring vegetation activity change. Based on simple and piecewise linear regression methods, this study made a comparative analysis between NDVIg and NDVI3g for monitoring vegetation activity change and its responses to climate change in the middle and high latitudes of the Northern Hemisphere during 1982–2008. Our results indicated that there were large differences between NDVIg and NDVI3g in the spatial patterns for both the overall changing trends and the timing of Turning Points (TP) in NDVI time series, which spread over almost the entire study region. The average NDVI trend from NDVI3g was almost twice as great as that from NDVIg and the detected average timing of TP from NDVI3g was about one year later. Although the general spatial patterns were consistent between two data sets for detecting the responses of growing-season NDVI to temperature and precipitation changes, there were large differences in the response magnitude, with a higher response magnitude to temperature in NDVI3g and an opposite response to precipitation change for the two data sets. These results demonstrated that the NDVIg data set may underestimate the vegetation activity change trend and its response to climate change in the middle and high latitudes of the Northern Hemisphere during the past three decades.

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.

Accurate data about the spatial distribution and planting area of maize is important for policy making, economic development, environmental protection and food security under climate change. This paper proposes a new identification method for spring maize based on spectral and phenological features derived from the moderate resolution imaging spectroradiometer (MODIS) land surface reflectance time-series data. The method focused on the spectral differences of different land cover types in the specific phenological phases of spring maize by testing the selections and combinations of classification metrics, feature extraction methods and classifiers. Taking Liaoning province, a representative planting region of spring maize in Northeast China, as the study area, the results indicated that the combined multiple metrics, including the red reflectance, near-infrared reflectance and normalized difference vegetation index (NDVI), were conducive to the maize identification and were better than any single metric. With regard to the feature extraction and selection, maize identification based on different phenological features selected with prior knowledge was more efficient than that based on statistical features derived from the principal component analysis. Compared with the maximum likelihood classification method, the decision tree classification based on expert knowledge was more suitable for phenological features selected from some prior knowledge. In summary, discriminant rules were defined with those phenological features from multiple metrics, and the decision tree classification was used to identify maize in the study area. The producer’s accuracy of maize identification was 98.57%, and the user’s accuracy was 81.18%. This method can be potentially applied to an operational identification of maize at large scales based on remote sensing time-series data.