Incorporating Vegetation Type Transformation with NDVI Time-Series to Study the Vegetation Dynamics in Xinjiang

Time-series normalized difference vegetation index (NDVI) is commonly used to conduct vegetation dynamics, which is an important research topic. However, few studies have focused on the relationship between vegetation type and NDVI changes. We investigated changes in vegetation in Xinjiang using linear regression of time-series MOD13Q1 NDVI data from 2001 to 2020. MCD12Q1 vegetation type data from 2001 to 2019 were used to analyze transformations among different vegetation types, and the relationship between the transformation of vegetation type and NDVI was analyzed. Approximately 63.29% of the vegetation showed no significant changes. In the vegetation-changed area, approximately 93.88% and 6.12% of the vegetation showed a significant increase and decrease in NDVI, respectively. Approximately 43,382.82 km2 of sparse vegetation and 25,915.44 km2 of grassland were transformed into grassland and cropland, respectively. Moreover, 17.4% of the area with transformed vegetation showed a significant increase in NDVI, whereas 14.61% showed a decrease in NDVI. Furthermore, in areas with NDVI increased, the mean NDVI slopes of pixels in which sparse vegetation transferred to cropland, sparse vegetation transferred to grassland, and grassland transferred to cropland were 9.8 and 3.2 times that of sparse vegetation, and 1.97 times that of grassland, respectively. In areas with decreased NDVI, the mean NDVI slopes of pixels in which cropland transferred to sparse vegetation, grassland transferred to sparse vegetation were 1.75 and 1.36 times that of sparse vegetation, respectively. The combination of vegetation type transformation NDVI time-series can assist in comprehensively understanding the vegetation change characteristics.

Estimating Actual Evapotranspiration over Croplands Using Vegetation Index Methods and Dynamic Harvested Area

Advances in estimating actual evapotranspiration (ETa) with remote sensing (RS) have contributed to improving hydrological, agricultural, and climatological studies. In this study, we evaluated the applicability of Vegetation-Index (VI) -based ETa (ET-VI) for mapping and monitoring drought in arid agricultural systems in a region where a lack of ground data hampers ETa work. To map ETa (2000–2019), ET-VIs were translated and localized using Landsat-derived 3- and 2-band Enhanced Vegetation Indices (EVI and EVI2) over croplands in the Zayandehrud River Basin (ZRB) in Iran. Since EVI and EVI2 were optimized for the MODerate Imaging Spectroradiometer (MODIS), using these VIs with Landsat sensors required a cross-sensor transformation to allow for their use in the ET-VI algorithm. The before- and after- impact of applying these empirical translation methods on the ETa estimations was examined. We also compared the effect of cropping patterns’ interannual change on the annual ETa rate using the maximum Normalized Difference Vegetation Index (NDVI) time series. The performance of the different ET-VIs products was then evaluated. Our results show that ETa estimates agreed well with each other and are all suitable to monitor ETa in the ZRB. Compared to ETc values, ETa estimations from MODIS-based continuity corrected Landsat-EVI (EVI2) (EVIMccL and EVI2MccL) performed slightly better across croplands than those of Landsat-EVI (EVI2) without transformation. The analysis of harvested areas and ET-VIs anomalies revealed a decline in the extent of cultivated areas and a loss of corresponding water resources downstream. The findings show the importance of continuity correction across sensors when using empirical algorithms designed and optimized for specific sensors. Our comprehensive ETa estimation of agricultural water use at 30 m spatial resolution provides an inexpensive monitoring tool for cropping areas and their water consumption.

A Rapid Model (COV_PSDI) for Winter Wheat Mapping in Fallow Rotation Area Using MODIS NDVI Time-Series Satellite Observations: The Case of the Heilonggang Region

Rapid and accurate monitoring of spatial distribution patterns of winter wheat over a long period is of great significance for crop yield prediction and farmland water consumption estimation. However, weather conditions and relatively long revisit cycles often result in an insufficient number of continuous medium-high resolution images over large areas for many years. In addition, the cropland pattern changes frequently in the fallow rotation area. A novel rapid mapping model for winter wheat based on the normalized difference vegetation index (NDVI) time-series coefficient of variation (NDVI_COVfp) and peak-slope difference index (PSDI) is proposed in this study. NDVI_COVfp uses the time-series index volatility to distinguish cultivated land from background land-cover types. PSDI combines the key growth stages of winter wheat phenology and special bimodal characteristics, substantially reducing the impact of abandoned land and other crops. Taking the Heilonggang as an example, this study carried out a rapid mapping of winter wheat for four consecutive years (2014–2017), and compared the proposed COV_PSDI with two state-of-the-art methods and traditional methods (the Spectral Angle Mapping (SAM) and the Iterative Self-Organizing Data Analysis Technique Algorithm (ISODATA)). The verification results revealed that the COV_PSDI model improved the overall accuracy (94.10%) by 4% compared with the two state-of-art methods (90.80%, 89.00%) and two traditional methods (90.70%, 87.70%). User accuracy was the highest, which was 93.74%. Compared with the other four methods, the percentage error (PE) of COV_PSDI for four years was the lowest in the same year, with the minimum variation range of PE being 1.6–3.6%. The other methods resulted in serious overestimation. This demonstrated the effectiveness and stability of the method proposed in the rapid and accurate extraction of winter wheat in a large area of fallow crop rotation region. Our study provides insight for remote sensing monitoring of spatiotemporal patterns of winter wheat and evaluation of “fallow rotation” policy implementation.

Temporal Information Extraction for Afforestation in the Middle Section of the Yarlung Zangbo River Using Time-Series Landsat Images Based on Google Earth Engine

Afforestation is one of the most efficient ways to control land desertification in the middle section of the Yarlung Zangbo River (YZR) valley. However, the lack of a quantitative way to record the planting time of artificial forest (AF) constrains further management for these forests. The long-term archived Landsat images (including the Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), and Operational Land Imager (OLI)) provide a good opportunity to capture the temporal change information about AF plantations. Under the condition that there would be an abrupt increasing trend in the normalized difference vegetation index (NDVI) time-series curve after afforestation, and this characteristic can be thought of as the indicator of the AF planting time. To extract the indicator, an algorithm based on the Google Earth Engine (GEE) for detecting this trend change point (TCP) on the maximum NDVI time series within the growing season (May to September) was proposed. In this algorithm, the time-series NDVI was initially smoothed and segmented into two subspaces. Then, a trend change indicator Sdiff was calculated with the difference between the fitting slopes of the subspaces before and after each target point. A self-adaptive method was applied to the NDVI series to find the right year with the maximum TCP, which is recorded as the AF planting time. Based on the proposed method, the AF planting time of the middle section of the YZR valley from 1988 to 2020 was derived. The detected afforestation temporal information was validated by 222 samples collected from the field survey, with a Pearson correlation coefficient of 0.93 and a root mean squared error (RMSE) of 2.95 years. Meanwhile, the area distribution of the AF planted each year has good temporal consistency with the implementation of the eco-reconstruction project. Overall, the study provides a good way to map AF planting times that is not only helpful for sustainable management of AF areas but also provides a basis for further research on the impact of afforestation on desertification control.

Land Cover and Crop Classification Based on Red Edge Indices Features of GF-6 WFV Time Series Data

Time series of vegetation indices can be utilized to capture crop phenology information, and have been widely used in land cover and crop classification, phenological feature extraction, and planting structure monitoring. This is of great significance for guiding agricultural production and formulating agricultural policies. According to the characteristics of the GF-6 satellite’s newly-added red edge bands, wide field view and high-frequency imaging, the time series of vegetation indices about multi-temporal GF-6 WFV data are used for the study of land cover and crop classification. In this study, eight time steps of GF-6 WFV data were selected from March to October 2019 in Hengshui City. The normalized difference vegetation index (NDVI) time series and 10 different red edge spectral indices time series were constructed. Then, based on principal component analysis (PCA), using two feature selection and evaluation methods, stepwise discriminant analysis (SDA) and random forest (RF), the red edge vegetation index of normalized difference red edge (NDRE) was selected. Seven different lengths of NDVI, NDRE and NDVI&NDRE time series were reconstructed by the Savizky-Golay (S-G) smoothing algorithm. Finally, an RF classification algorithm was used to analyze the influence of time series length and red edge indices features on land cover and crop classification, and the planting structure and distribution of crops in the study area were obtained. The results show that: (1) Compared with the NDRE red edge time series, the NDVI time series is more conducive to the improvement of the overall classification accuracy of crops, and NDRE can assist NDVI in improving the crop classification accuracy; (2) With the shortening of NDVI and NDRE time series, the accuracy of crop classification is gradually decreased, and the decline is gradually accelerated; and (3) Through the combination of the NDVI and NDRE time series, the accuracy of crop classification with different time series lengths can be improved compared with the single NDVI time series, which is conducive to improving the classification accuracy and timeliness of crops. This study has fully tapped the application potential of the new red edge bands of GF-6 WFV time series data, which can provide references for crop identification and classification of time series data such as NDVI and red edge vegetation index of different lengths. At the same time, it promotes the application of optical satellite data with red edge bands in the field of agricultural remote sensing.

Estimation of the Conifer-Broadleaf Ratio in Mixed Forests Based on Time-Series Data

Most natural forests are mixed forests, a mixed broadleaf-conifer forest is essentially a heterogeneously mixed pixel in remote sensing images. Satellite missions rely on modeling to acquire regional or global vegetation parameter products. However, these retrieval models often assume homogeneous conditions at the pixel level, resulting in a decrease in the inversion accuracy, which is an issue for heterogeneous forests. Therefore, information on the canopy composition of a mixed forest is the basis for accurately retrieving vegetation parameters using remote sensing. Medium and high spatial resolution multispectral time-series data are important sources for canopy conifer-broadleaf ratio estimation because these data have a high frequency and wide coverage. This paper highlights a successful method for estimating the conifer-broadleaf ratio in a mixed forest with diverse tree species and complex canopy structures. Experiments were conducted in the Purple Mountain, Nanjing, Jiangsu Province of China, where we collected leaf area index (LAI) time-series and forest sample plot inventory data. Based on the Invertible Forest Reflectance Model (INFORM), we simulated the normalized difference vegetation index (NDVI) time-series of different conifer-broadleaf ratios. A time-series similarity analysis was performed to determine the typical separable conifer-broadleaf ratios. Fifteen Gaofen-1 (GF-1) satellite images of 2015 were acquired. The conifer-broadleaf ratio estimation was based on the GF-1 NDVI time-series and semi-supervised k-means cluster method, which yielded a high overall accuracy of 83.75%. This study demonstrates the feasibility of accurately estimating separable conifer-broadleaf ratios using field measurement data and GF-1 time series in mixed broadleaf-conifer forests.

Optimal Estimate of Global Biome—Specific Parameter Settings to Reconstruct NDVI Time Series with the Harmonic ANalysis of Time Series (HANTS) Method

Terrestrial remote sensing data products retrieved from radiometric measurements in the optical and thermal infrared spectrum such as vegetation spectral indices can be heavily contaminated by atmospheric conditions, including cloud and aerosol layers. This contamination results in gaps or noisy observations. The harmonic analysis of time series (HANTS) has been widely used for time series reconstruction of remote sensing imagery in recent decades. To use HANTS model, a series of parameters, such as number of frequencies (NF), fitting error tolerance (FET), degree of over-determinedness (DoD), and regularization factor (Delta), need to be defined by users. These parameters provide flexibilities, but also make it difficult for non-expert users to determine appropriate settings for specific applications. This study systematically evaluated the reconstruction performance of the model under different parameter setting scenarios by simulating pixel-wise reference and noisy NDVI time series. The results of these numerical experiments were further used to identify optimal settings and improve global NDVI reconstruction performance. The results suggested optimal settings for different areas (local optimization). If a user opts to use unique settings for global reconstruction, the setting NF = 4, FET = 0.05, DoD = 5, and Delta = 0.5 can produce the best performance across all setting scenarios (global optimization). In addition, several internal improvements, such as dynamic weighting scheme, polynomial and inter-annual harmonic components, and ancillary attributes of input data can be used to further improve the performance of reconstruction. With these results, future non-expert users can easily determine appropriate settings of HANTS for specific applications in different regions.