Next Article in Journal
Remote Marker-Based Tracking for UAV Landing Using Visible-Light Camera Sensor
Next Article in Special Issue
Village Building Identification Based on Ensemble Convolutional Neural Networks
Previous Article in Journal
GPS Satellite Orbit Prediction at User End for Real-Time PPP System
Previous Article in Special Issue
A Synergy Cropland of China by Fusing Multiple Existing Maps and Statistics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Differences in Phenology Extracted from the Enhanced Vegetation Index and the Leaf Area Index

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Authors to whom correspondence should be addressed.
Sensors 2017, 17(9), 1982; https://doi.org/10.3390/s17091982
Submission received: 21 June 2017 / Revised: 24 August 2017 / Accepted: 28 August 2017 / Published: 30 August 2017
(This article belongs to the Special Issue Sensors and Smart Sensing of Agricultural Land Systems)

Abstract

:
Remote-sensing phenology detection can compensate for deficiencies in field observations and has the advantage of capturing the continuous expression of phenology on a large scale. However, there is some variability in the results of remote-sensing phenology detection derived from different vegetation parameters in satellite time-series data. Since the enhanced vegetation index (EVI) and the leaf area index (LAI) are the most widely used vegetation parameters for remote-sensing phenology extraction, this paper aims to assess the differences in phenological information extracted from EVI and LAI time series and to explore whether either index performs well for all vegetation types on a large scale. To this end, a GLASS (Global Land Surface Satellite Product)-LAI-based phenology product (GLP) was generated using the same algorithm as the MODIS (Moderate Resolution Imaging Spectroradiometer)-EVI phenology product (MLCD) over China from 2001 to 2012. The two phenology products were compared in China for different vegetation types and evaluated using ground observations. The results show that the ratio of missing data is 8.3% for the GLP, which is less than the 22.8% for the MLCD. The differences between the GLP and the MLCD become stronger as the latitude decreases, which also vary among different vegetation types. The start of the growing season (SOS) of the GLP is earlier than that of the MLCD in most vegetation types, and the end of the growing season (EOS) of the GLP is generally later than that of the MLCD. Based on ground observations, it can be suggested that the GLP performs better than the MLCD in evergreen needleleaved forests and croplands, while the MLCD performs better than the GLP in shrublands and grasslands.

1. Introduction

Phenology is the study of the timing of recurring biological events, the causes of their timing with regard to biotic and abiotic forces, and the interactions among phases of the same or different species [1]. Sprouting and flowering of plants in spring, color changes of leaves in fall, bird migration and nesting, insect hatching, and animal hibernation are all examples of phenological events [2]. Phenology is an important indicator of global change and the carbon cycle because of its direct effects on vegetation photosynthesis, carbon sequestration and land–atmosphere water and energy exchange [3]. Vegetation phenology indicates the responses of plants to seasonal and interannual variations of climate, hydrology, soil and anthropogenic factors [4]. Vegetation phenology also affects the climate system by influencing the seasonality of albedo; surface roughness length; canopy conductance; and fluxes of water, energy, CO2 and biogenic volatile organic compounds [5]. The extraction of vegetation phenology is important for research on global climate change and terrestrial ecosystems [6]. The phenology products retrieved from remote-sensing data have been widely applied in many fields, such as climatic change [7,8,9,10], biomass monitoring [11,12,13,14] and farm management [15,16,17,18].
Remote-sensing vegetation phenological metrics are generally derived from satellite time-series data of vegetation parameters. Vegetation indexes (VIs) are one of the most commonly used parameters, including the normalized difference vegetation index (NDVI) [4,19,20,21] and enhanced vegetation index (EVI) [22,23,24]. Generally, the NDVI tends to lose sensitivity over dense canopies because of saturation [25], while the EVI has a larger dynamic range than the NDVI and is more resistant to atmospheric and soil background effects [26]. Many studies have explored other VIs to indicate the growing season transitions, such as the soil adjusted vegetation index (SAVI) [27], MERIS terrestrial chlorophyll index (MTCI) [28], perpendicular vegetation index (PVI) [9] and wide dynamic range vegetation index (WDRVI) [25]. While VI is an integrated indicator of the vegetation, the soil background, and the observation-illumination geometry, some vegetation physiological parameters have also been used to detect vegetation phenology, such as the leaf area index (LAI) [3], the fraction of photosynthetically active radiation (fPAR) [29] and albedo [30]. Compared with the VIs, these parameters provide more explicit information on the biophysical characteristics of vegetation changes. Among these parameters, the LAI is most widely used because of the availability of multiple remote-sensing products and its direct indication of physical and biological processes related to vegetation dynamics at global and regional scales, such as energy exchange and water and carbon cycling [31,32].
In addition, many methods have been developed to monitor vegetation phenology using satellite time-series data. The methods can be grouped into six classes: threshold-based, delayed moving average, curve fitting, curve derivative, phenological cumulative frequency and principal component analysis. The threshold-based method defines the start of season (SOS) using a predefined or relative reference value [4,33,34]. The delayed moving-average method determines the SOS based on an autoregressive moving average (CARMA) model [35]. The curve-fitting method fits the time-series remote-sensing data with a smooth model and then determines vegetation phenology using the fitted function [22,36,37,38,39,40,41]. The curve-derivative method defines the time point of the time-series curve with the maximal second derivative value as the SOS [42]. The phenological cumulative-frequency method uses field phenology observations and the phenological frequency distribution pattern to determine the threshold of plant phenology events [43]. The principal component analysis method extracts vegetation phenology based on the characteristics of the first principal component of the datasets employing empirical orthogonal function (EOF) decomposition with NDVI [44].
Compared with the strong demand for vegetation phenology mapping for global and regional research and applications, the development of remote-sensing vegetation phenology products has been slower than other remote-sensing parameters, such as the LAI and VIs. Few regional-to-global remote-sensing vegetation phenology products have been generated despite the value of vegetation phenology for Earth-system monitoring and modeling [45,46]. The MODIS (Moderate Resolution Imaging Spectroradiometer) Global Land Cover Dynamics product (MLCD) is the only standard global product and is generated by NASA from MODIS-EVI time-series data [22,47].
The MLCD is based on the piecewise logistic model, which is the most widely used model in remote-sensing phenology extraction as it can effectively restrain the noise effect and requires no threshold or empirical constraints. Scholars have performed some validations for the MLCD. Comparisons with ground data suggest that the MLCD performs well, but the end-of-season metrics associated with vegetation senescence and dormancy have higher uncertainties than the start-of-season metrics [47]. In evergreen forests, the MLCD is unable to provide consistent phenologic patterns, and considerable data are missing, especially in tropical rainforests [48]. Against this background, Verger et al. [49] characterized the baseline phenology of the vegetation at the global scale from the GEOCLIM climatology of the LAI estimated from the 1-km SPOT-VEGETATION time series, where the start- and end-of-season were identified using 30% and 40% thresholds of the LAI amplitude values, respectively. Although the spatial patterns of the LAI-derived phenology agreed well with those from the MODIS-EVI and -NDVI, the timing of the start, end, and length of season differed by approximately one month at the global scale. The higher uncertainties appear in areas with limited seasonality expressed in the satellite signal and systematic biases due to the differences in the methodologies and datasets. Some other studies have also shown significant discrepancies in vegetation phenological metrics derived from time series of different vegetation parameters [50,51,52]. However, it is not clear which data sources and vegetation parameters are best suited for extracting vegetation phenological metrics accurately.
While EVI and LAI are the most widely used vegetation parameters that can be used for remote-sensing phenological extraction, this paper aims to assess the differences in phenological information extracted using the EVI and LAI time series and to determine whether either EVI or LAI time series performs well for all vegetation types over a large scale. The objective is to compare the phenology detection results in China from the MODIS-EVI and the GLASS-LAI. To achieve this objective, the GLASS-LAI phenology product (GLP) was generated using the same algorithm as the MODIS-EVI phenology product (MLCD) in China from 2001 to 2012. The two phenology products were compared in different vegetation types in China and then evaluated using ground observations from forests, grasslands and croplands. Our work lays the foundation for uniting multisource data and for improving remote-sensing phenology products in the future.

2. Data and Methods

2.1. Phenology Extraction Based on the EVI and LAI

The MLCD product was used as the vegetation phenology extracted from the EVI for analysis in the paper. The MLCD provides yearly 500-m resolution global estimates of vegetation phenology [47] and is available for the public at https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mcd12q2. The MLCD identifies phenophase transition dates based on the logistic functions fit-to-time series of the EVI, which were calculated from composite 8-day normalized BRDF (Bidirectional Reflectance Distribution Function)-adjusted reflectance data [26]. Specifically, a time series of the EVI is assembled for each pixel, the data undergo a gap-filling and smoothing process, periods of sustained EVI increase or decrease are identified, logistic models are fit to the time series, and transition dates are identified as local maxima and minima in the rate of curvature change of the fitted logistic function [22]. The MLCD has four layers corresponding to key phenological stages: onset greenness increase, onset greenness maximum, onset greenness decrease, and onset greenness minimum. MLCD images over China from 2001 to 2012 were acquired for the study. Due to its different spatial resolution from the GLP, the MLCD was upscaled to 1000 m by averaging the value from every four pixels for statistics of different values between GLP and MLCD. But in other analyses, MLCD was used with its original spatial resolution.
The vegetation phenology extracted from the LAI was generated from the GLASS-LAI (http://glass-product.bnu.edu.cn/) that provides yearly 1000-m resolution global estimates of vegetation phenology, which is called the GLASS-LAI-based phenology product (GLP). GLASS-LAI was estimated from MODIS and CYCLOPES LAI time series using general regression neural networks, with 8-day temporal resolution and 1-km spatial resolution. Based on multitemporal data, GLASS-LAI are spatially and temporally continuous with no gaps in time series and no high-frequency temporal variability. As a result of cloud, atmospheric and snow contamination, Savitzky–Golay filtering was used to smooth and gap-fill the CYCLOPES LAI. The MODIS reflectance was reprocessed to remove remaining effects of cloud contamination and other factors [53]. The GLP applied the same algorithm as MLCD (the piecewise logistic model) to fit the LAI time series, which is a flexible, repeatable, and realistic means of monitoring seasonal and interannual vegetation dynamics from remote-sensing data at large scales. The piecewise logistic equation is represented as Equation (1):
y ( t ) = c 1 + e a + b t + d
where t is time in days, y(t) is the LAI value at time t, a and b are fitting parameters, c + d is the maximum LAI value, and d is the initial background LAI value. GLP products with 1-km resolution were generated over China from 2001 to 2012.
For both the MLCD and GLP products, the vegetation dynamics can be divided into four stages: greenup, maturity, senescence and dormancy [22], as shown in the Figure 1. The key transition dates define the main phases of vegetation dynamics, including the start of the growing season (SOS), the start of the peak (PS), the end of the peak (PE) and the end of the growing season (EOS): (1) SOS is the date of onset of photosynthetic activity; (2) PS is the date when the greenness reaches its maximum; (3) PE represents the date when the photosynthetic activity and the greenness begin to decrease rapidly; and (4) EOS is the date when the physiological activity approaches zero [14,54].

2.2. Remote-Sensing Auxiliary Datasets

The nadir BRDF-adjusted reflectance (NBAR) data from MODIS (MCD43A4) were used to compute the EVI, as shown in Equation (2). The MODIS-EVI products and the GLASS-LAI products were used to compare the difference of vegetation growing states reflected by the LAI and EVI time series. Besides, the land cover type (MOD12Q1) was downloaded from NASA’s Earth Observing System Data Center website. The MOD12Q1 was used to identify the land cover type over pixels, as shown in Figure 2.
EVI = 2.5 × ρ N I R ρ R E D ρ N I R + 6.0 ρ R E D 7.5 ρ B L U E + 1
where ρ N I R , ρ R E D and ρ B L U E are the reflectances of near infrared, red and blue bands respectively.

2.3. Ground Phenological Observations

The ground phenological observations were taken as the true values to evaluate the accuracy of the phenology products. However, the ground-based phenology was observed in a number of individual plants, while the remote-sensing phenology represented the integrated phenological characteristics of a plant community in one pixel. Ground validation of remote-sensing measurements with coarse resolution entails considerable difficulties. To improve the reliability of the statistical analysis based on ground observations, ground phenological observations meeting the following requirements were selected from the China Meteorological Administration (CMA) and the Chinese Ecosystem Research Network (CERN):
(1)
Spatial representation. The poor relationship between ground and satellite phenology due to data-scale issues is a drawback of satellite phenology because of the small chance of a single-point ground observation being representative of an entire area at the remote-sensing scale (typically ≥1 km in remote-sensing phenological studies) [4,35]. Consequently, the phenological homogeneity and subdued topography of field sites must be ensured in comparison with the remote-sensing data [41]. The phenological homogeneity requires the phenophases of dominant species at one site be similar (less than 30 days), and the topography is checked by Google Earth to avoid sites in mountains as far as possible. The dominant species were selected after considering the distribution and quantity of the community based on the instruction files from CERN and CMA.
(2)
Data integrity. The selected ground sites should have phenological phase continuity and few missing records.
The ground phenological observations were selected according to the above requirements, as shown in Table 1. Each site corresponds to a typical vegetation type, including evergreen needleleaved forest (ENF), evergreen broadleaved forest (EBF), deciduous needleleaved forest (DNF), deciduous broadleaved forest (DBF), shrub, grass and crop. Rice was chosen as a representative crop with two sites, because the SY (Shenyang) site from CERN can provide harvest and JT (Jiutai) site from CMA can provide grain-filling. The mean phenological dates of the dominant species were calculated and treated representative of this region. The distribution of the ground phenological observation sites is shown in Figure 2.

2.4. Methods for Evaluating the Phenology Products

To analyze the differences in the GLP and the MLCD, the mean phenological metrics for each pixel in each year were calculated as the spatial pattern of China. Furthermore, the mean and standard deviation of the difference values between the two products were calculated based on land cover type data. In addition, the mean bias (MB) and mean absolute bias (MAB) were used to evaluate the remote-sensing phenological metrics based on ground-station observations [55], as shown in Equations (3) and (4).
M B = i = 1 N ( r s o b s ) N
M A B = i = 1 N | r s o b s | N
where rs is the remote-sensing phenological date; obs is the ground observation; and N is the number of rs or obs.

3. Analysis and Results

3.1. Comparison between the GLP and the MLCD over China

Figure 3a–d presents the mean phenological metrics for China for the 2001–2012 period derived from the GLP and MLCD data. In general, the GLP and the MLCD have consistent patterns in the north of China, but inconsistent patterns are exhibited in the south. For both the GLP and the MLCD, the missing data exist inevitably and vary in different regions.

3.1.1. Missing Data

In general, the ratio of missing data is 8.3% for the GLP, which is obviously less than the 22.8% for the MLCD, as shown in Figure 4. The missing ratios of the MLCD and the GLP differ obviously over different vegetation types. For the MLCD product, the missing ratio is highest in EBF (53.6%), followed by ENF (48.73%), which is much higher than that of GLP (6.2% and 4.0% respectively). Nevertheless, the missing ratio of shrub is the highest for GLP (55.3%), which is higher than that of MLCD (14.9%). In addition, the missing ratio of MLCD in crop is relatively high, especially for the double-cycle crop in the North China Plain (Figure 2); invalid results are also treated as missing data.
As shown in Figure 5, the patterns of the missing ratios depend strongly on the latitude between 22° N and 53° N in China. In general, the missing ratios of the MLCD are higher than those of the GLP along each latitude zone except from 37° N to 43° N. The missing ratio of the MLCD gradually decreases as the latitude decreases, especially in midlatitude zones (30° N–53° N), and remained inordinately high between 22° N and 30° N (approximately 50%), where much data are missing in the Sichuan Basin (as shown in Figure 2). However, the missing ratio of GLP largely depends on desert vegetation. For the GLP, the latitudinal distribution of missing data is consistent with the distribution of grass or shrub, as shown in Figure 5 and Figure 6.
There are several factors causing the missing data. The first is incomplete of time-series source data caused by clouds or snow. Cloud cover severely compromises the quality of NBAR EVI time series in extensive areas of tropics and subtropics [56], as a consequence, many phenological metrics are missing in these areas. In addition, the seasonal change characteristics of evergreen forest are not obvious, so it is hard to extract the phenology information of EBF and ENF. While GLASS-LAI exhibits spatial-temporal continuity [53], there is no gaps between GLASS-LAI time series. This is the main reason why the missing ratio of GLP is less than that of MLCD generally, but the missing ratio of GLP in desert vegetation (e.g., shrub) is relatively high because the range of the seasonal variation in the LAI time series is small. Therefore, fitting by the logistics method fails easily.

3.1.2. Difference Comparison between the GLP and the MLCD

Specifically, the difference between the GLP and the MLCD varies among different vegetation types, as shown in Figure 7. Overall, the SOS of the GLP is earlier than that retrieved by the MLCD in most vegetation types, while the EOS of the GLP is later, including shrub. Among the vegetation types, the difference between the GLP and the MLCD in deciduous broadleaved forest is the smallest, among which the difference in SOS is smaller than that of EOS (the means and standard deviations are −7.2 ± 17.7 and 5.5 ± 36.5, respectively). The greatest difference occurs in evergreen broadleaved forest, among which the difference in SOS is larger than that of EOS (the means and standard deviations are −84.3 ± 76.5 and 38.1 ± 73.3, respectively). In addition, the difference in SOS or EOS between the GLP and the MLCD in croplands is relatively large (the means and standard deviations are −14.5 ± 53.6 and 38.8 ± 61.9, respectively), as double-cycle crops are poorly detected in MLCD.
Obviously, the differences between the GLP and the MLCD become stronger as the latitude decreases, as shown in Figure 2 and Figure 8. The difference between the GLP and the MLCD largely depends on different vegetation types. The greatest difference between the GLP and the MLCD is found below 26° N, where EBF is mainly distributed. As the ratio of EBF decreases, the difference decreases. Because DBF is mainly distributed above 40° N, the difference between the GLP and the MLCD become small above 40° N. The difference in the SOS below 40° N is larger than that in the EOS, and the differences in the SOS and the EOS are similar above 40° N.

3.2. Accuracy Assessment Based on Field Data

To analyze the variances and relationships between ground observations and the phenological information extracted from the MODIS-EVI and GLASS-LAI products, the EVI and LAI time series at representative sites were chosen as reference to compare the phenological results, as shown in Figure 9. The information from ground phenology observation sites is presented in Table 1.
For evergreen forests, significant discrepancies occur in the LAI and EVI time series (Figure 9a,b). Obviously, irregular fluctuations exist in the EVI of ENF and EBF, which are more likely due to a pattern of noise caused by clouds in the EVI data structure rather than to vegetation growing changes. Thus, the phenological metrics extracted by MODIS-EVI may be discontinuous, or even missing, especially in the tropical climate zone. Apart from noise, in many areas with evergreen vegetation, the seasonal change characteristics in EVI is too subtle to retrieve phenology.
For deciduous forests, the LAI time series agree with the EVI time series overall, but the LAI time series have fewer fluctuations than the EVI time series in the peak (Figure 8c,d). The differences between LAI and EVI time series in DNF is larger than that in DBF. In DNF, although the LAI start to increase earlier than the EVI, the EVI is more sensitive than the LAI during the greenup. Therefore, the EVI generally reaches maturity faster than the LAI. And also, the EVI start to decrease later than the LAI, but decreases rapidly to the minimum.
The LAI and the EVI of shrublands are relatively lower. Although the range of seasonal variation of the EVI time series is smaller than that of the LAI, the precision of the current LAI products is 0.1, which is easier for presenting the “ladder” time series than the EVI product, whose precision is 0.01. Thus, the EVI shows the subtle dynamic changes of grasslands more clearly. As shown in Figure 9e, the inflection point of the EVI time series is close to bud burst and leaf defoliation of shrub. When using the LAI, it is difficult to show the subtle dynamic changes, causing a time offset in shrub phenology. Especially for shrubs in desert, the ranges of the seasonal variations in EVI or LAI are so small that the logistics model is difficult to fit to extract phenological information.
For croplands, although the LAI start to increase earlier than the EVI, the EVI is more sensitive to vegetation growth than the LAI during the greenup. Therefore, he EVI generally reaches maturity faster than the LAI, but the peak of the LAI time series is closer to the ground observations, as shown in Figure 9g,h.
In order to take the ground phenology observations as the truth values to evaluate the accuracy of the phenology products, the phenophases observed from ground observations corresponding to remote-sensing phenological metrics are determined, as shown in Table 2. All the observed ground phenology stages are the beginning of the phenophases in both CERN and CMA, so the phenophases observed from ground observations that were closest to the remote-sensing phenological metrics were chosen based on Figure 9. According to previous studies, the key remote-sensing phenological metrics might correspond to the four growing stages for crop [57,58]. For herbs and deciduous broadleaf trees, no phenophases correspond to PS and PE.
Figure 10 compares the MB and MAB values between the remote-sensing phenological metrics and the ground data at typical vegetation sites.
For evergreen forests, the phenological metrics of MLCD are missing at the HS and DHS sites; only the phenological metrics of GLP are presented, as shown in Figure 10a,b. Specifically, LAI can extract the phenology information of ENF especially for SOS with high accuracy. The MAB of the GLP for the SOS is 4.3 days (MB = 1.0 day), and the MAB of the GLP for the EOS is 13.6 days (MB = −1.0 day). However, GLP performs poorly in EBF, the MAB of the GLP for the SOS is 39.2 days (MB = 39.2 days), and the MAB of the GLP for the EOS is 34.0 days (MB = 34.0 days).
For deciduous forests, the accuracy of GLP or MLCD for SOS is higher than that for EOS, as shown in Figure 10c,d. The SOS of MLCD in DBF is very close to the 1:1 line, suggesting that the SOS that extracted by EVI is in good agreement with ground observations. Specifically, the MAB of the MLCD for the SOS is 4.1 days (MB = −1.9 days), and the MAB of the GLP for the SOS is 9.0 days (MB = −8.5 days). The EOS of GLP is later than that of MLCD, and their MABs are 22.9 days and 17.3 days respectively.
For shrublands, because LAI is hard to extract the phenology information in desert at SPT (Shapotou) site, only a few phenological metrics of GLP are presented, as shown in Figure 10e. Specifically, the MABs of the MLCD for the SOS and EOS are 17.9 days (MB = −4.4 days) and 20.6 days (MB = −11.4 days). Compared to shrublands, the MABs of the MLCD and the GLP are smaller. The MAB of the MLCD for the SOS is 15.1 days (MB = −8.6 days), and the MAB of the GLP for the SOS is 14.1 days (MB = 8.1 days). The EOS of GLP is later than that of MLCD, and their MABs are 24.0 days and 16.4 days respectively.
For croplands, the distributions of MLCD for SOS and EOS are similar to that of GLP. The MABs of the MLCD and the GLP for the SOS are 11.0 days (MB = 10.5 days) and 10.4 days (MB = 10.4 days) respectively. The MABs of the MLCD and the GLP for the EOS are 8.1 days (MB = 7.9 days) and 7.2 days (MB = 7.0 days) respectively. However, the accuracy of GLP for the PS is much higher than that of MLCD, since the MAB of the GLP and the MLCD for the PS is 6.1 days (MB = −5.1 days) and 18.2 days (MB = 18.2 days) respectively. Compared to the PS, the MABs of both phenology products are larger for PE.

4. Discussion and Conclusions

This paper assessed the differences in phenological information extracted from the EVI and LAI time series and explored whether the EVI or the LAI time series performs well for all vegetation types over a large scale in China. The GLP and the MLCD were compared in different regions and evaluated based on ground observations.
The missing ratio of the MLCD is more than that of the GLP overall. The main reason is incompletion of time-series source data caused by clouds or snow for the MLCD product. While GLP is not affected by clouds or snow, because the GLASS-LAI is provided with spatial-temporal continuity, there are no gaps in time series and no high-frequency temporal variability. However, the missing ratio of GLP in desert vegetation (e.g., shrub) is relatively high. Because the range of seasonal variation of LAI time series is small, fitting by the logistics method easily fails. Besides, in many areas with evergreen vegetation, the annual variation in EVI is too subtle to retrieve phenology.
Accuracy assessment is always an important concern in any remote-sensing-based analysis, especially one performed at a large scale. To improve the reliability of the statistical analysis based on ground observations, we selected sites based on the requirements of spatial representation and data integrity. The accuracy assessment based on representative sites in this paper show that both phenology products can provide the best accuracy with ground observations for the SOS of DBF and the EOS of rice. Because the leaf defoliation period for DBF lasts for one or two months, and the EOS extracted from LAI may be close to the end of leaf defoliation, a lag can occur between the EOS and the start of leaf defoliation from ground observations. Besides, previous research has shown that the end-of-season metrics of the MLCD associated with vegetation senescence and dormancy are highly uncertain [47,56]. In this paper, it is shown that the EOS in croplands has less uncertainties than that in forests. Thus, it can be inferred that the uncertainties of EOS in forest may come from its complex structure and composition, as the composition of cropland is relatively homogeneous. For grasslands or shrublands, the range of seasonal variation is relatively small, the logistics model couldn’t fit well with EVI or LAI time series; as a consequence, the accuracy of grasslands or shrublands is worse than that of deciduous forests or croplands. LAI is easier for presenting the “ladder” time series than the EVI and even failing in fitting when the maximum of LAI is less than 0.3. In addition, the growth of grass or shrub in arid regions or desert is sensitive to precipitation, which may result in irregular fluctuations in EVI time series. Our analysis also shows that the GLP performs better than the MLCD in ENF and croplands, especially for the PS and PE of crop. Although GLP has less missing data in EBF, the disagreement exist between vegetation phenological metrics and ground observations. However, the MLCD performs better than the GLP in shrublands and grasslands. The precision of the current LAI product is 0.1, indicating that it does not show subtle dynamic changes as well as the EVI product, whose precision is 0.01.
It is noted that MLCD is derived from 500-m, 8-day composite EVI data; the GLP is derived from 1000-m, 8-day composite LAI data; while ground observations are daily point-based observations. The temporal and spatial scales of ground observations make them difficult to compare directly with remote-sensing retrievals. Although studied sites are the best representation for each station and the surrounding area, mismatches between phenology products and ground-based observations at both temporal and spatial scales are unavoidable [56]. In addition, the monitoring method may lead to a mismatch with the ground observations; therefore, we used a flexible, repeatable, and realistic means (the piecewise logistic models) to monitor seasonal and inter-annual dynamics in vegetation from remote-sensing data at large scales. However, the method may also lead to an offset in some situations when the remote-sensing time series do not match the piecewise logistic models.
Besides, even though based on the same algorithm and mathematical curvature changes for the same spatial/vegetation type, phenological metrics derived from EVI and LAI time series may have different biophysical meanings. It seems that the SOS and EOS that are extracted from LAI are close to bud burst and leaf defoliation respectively, while the SOS and EOS that are extracted from EVI is close to leaf unfolding and leaf coloring respectively at deciduous broadleaved forest site. Thus, the SOS that is extracted from LAI is earlier than that from EVI, while the EOS that is extracted from LAI is later than that from EVI overall according to the analysis results.
In conclusion, both the LAI and the EVI have advantages in different vegetation types. It should be considered that uniting multisource data is an effective way to improve the accuracy and validity of remote-sensing phenology products. However, for evergreen forests and sparse vegetation, extracting phenological metrics is still challenging using the LAI or EVI time series based on piecewise logistic models. New methods or data sources must be adopted for these land cover types. In addition, due to the coarse spatial resolution of current remote-sensing phenology products, other plants are inevitably present within one square kilometer, which may lead to variability in phenological developmental stages and a weak relationship between remote-sensing data and ground observations. These findings are of great value for improving the spatial resolution of remote-sensing phenology products to promote their application and development in the future.

Acknowledgments

This work was supported by the National Basic Research Program of China (No. 2013CB733401), and the National Natural Science Foundation of China (No. 41671374).

Author Contributions

C.W. designed the research, performed the analyses and wrote the paper. J.L. and Q.L. conceived the research, provided the theoretical guidance for the subject study and reviewed and edited this paper. B.Z. and S.W. collected and preprocessed the remote sensing data. C.X. collected the field data and wrote the code of the model.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Helmut, L. Phenology and seasonality modeling. Soil Sci. 1975, 120, 461. [Google Scholar]
  2. Dubé, P.A.; Perry, L.P.; Vittum, M.T. Instructions for phenological observations: Lilac and honeysuckle. In Vermont Agricultural Experiment Station Bulletin 692; University of Vermont: Burlington, VT, USA, 1984. [Google Scholar]
  3. Peñuelas, J.; Rutishauser, T.; Filella, I. Phenology feedbacks on climate change. Science 2009, 324, 887–888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. White, M.A.; Thomton, P.E.; Running, S.W. A continental phenology model for monitoringvegetation responses to interannual climatic variability. Glob. Biogeochem. Cycles 1997, 11, 217–234. [Google Scholar] [CrossRef]
  5. Richardson, A.D.; Keenan, T.F.; Migliavacca, M.; Ryu, Y.; Sonnentag, O.; Toomey, M. Climate change, phenology, and phenological control of vegetation feedbacks to the climate system. Agric. For. Meteorol. 2013, 169, 156–173. [Google Scholar] [CrossRef]
  6. Chen, X.; Hu, B.; Yu, R. Spatial and temporal variation of phenological growing season and climate change impacts in temperate eastern china. Glob. Chang. Biol. 2005, 11, 1118–1130. [Google Scholar] [CrossRef]
  7. Suzuki, R.; Nomaki, T.; Yasunari, T. West-east contrast of phenology and climate in northern asia revealed using a remotely sensed vegetation index. Int. J. Biometeorol. 2003, 47, 126–138. [Google Scholar] [PubMed]
  8. Tateishi, R.; Ebata, M. Analysis of phenological change patterns using 1982–2000 advanced very high resolution radiometer (avhrr) data. Int. J. Remote Sens. 2004, 25, 2287–2300. [Google Scholar] [CrossRef]
  9. Guyon, D.; Guillot, M.; Vitasse, Y.; Cardot, H.; Hagolle, O.; Delzon, S.; Wigneron, J.-P. Monitoring elevation variations in leaf phenology of deciduous broadleaf forests from spot/vegetation time-series. Remote Sens. Environ. 2011, 115, 615–627. [Google Scholar] [CrossRef]
  10. Stöckli, R.; Vidale, P.L. European plant phenology and climate as seen in a 20-year avhrr land-surface parameter dataset. Int. J. Remote Sens. 2004, 25, 3303–3330. [Google Scholar] [CrossRef]
  11. Xiao, X.; Hollinger, D.; Aber, J.; Goltz, M.; Davidson, E.A.; Zhang, Q.; Moore, B., III. Satellite-based modeling of gross primary production in an evergreen needleleaf forest. Remote Sens. Environ. 2004, 89, 519–534. [Google Scholar] [CrossRef]
  12. Kimball, J.; Zhao, M.; McDonald, K.; Running, S. Satellite remote sensing of terrestrial net primary production for the pan-arctic basin and Alaska. Mitig. Adapt. Strateg. Glob. Chang. 2006, 11, 783–804. [Google Scholar] [CrossRef]
  13. Picard, G.; Quegan, S.; Delbart, N.; Lomas, M.R.; Toan, T.; Woodward, F. Bud-burst modelling in siberia and its impact on quantifying the carbon budget. Glob. Chang. Biol. 2005, 11, 2164–2176. [Google Scholar] [CrossRef]
  14. Galina, C.; David, S.; Braswell, B.H.; Xiao, X. Spatial analysis of growing season length control over net ecosystem exchange. Glob. Chang. Biol. 2005, 11, 1777–1787. [Google Scholar]
  15. Duchemin, B.; Hadria, R.; Rodriguez, J.C.; Lahrouni, A.; Khabba, S.; Boulet, G.; Mougenot, B.; Maisongrande, P.; Watts, C. Spatialisation of a crop model using phenology derived from remote sensing data. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS ’03), Toulouse, France, 21–25 July 2003; pp. 2200–2202. [Google Scholar]
  16. Lobell, D.B.; Asner, G.P.; Ortiz-Monasterio, J.I.; Benning, T.L. Remote sensing of regional crop production in the yaqui valley, mexico: Estimates and uncertainties. Agric. Ecosyst. Environ. 2003, 94, 205–220. [Google Scholar] [CrossRef]
  17. Fan, J.-L.; Wu, B.-F. A methodology for retrieving cropping index from NDVI profile. J. Remote Sens. 2004, 8, 628–636. [Google Scholar]
  18. De Beurs, K.M.; Henebry, G.M. Land surface phenology, climatic variation, and institutional change: Analyzing agricultural land cover change in Kazakhstan. Remote Sens. Environ. 2004, 89, 497–509. [Google Scholar] [CrossRef]
  19. Schwartz, M.D.; Reed, B.C.; White, M.A. Assessing satellite derived start of season measures in the conterminous USA. Int. J. Climatol. 2002, 22, 1793–1805. [Google Scholar] [CrossRef]
  20. Goward, S.N.; Tucker, C.J.; Dye, D.G. North American vegetation patterns observed with the NOAA-7 advanced very high resolution radiometer. Vegetatio 1985, 64, 3–14. [Google Scholar] [CrossRef]
  21. White, M.A.; Schwartz, M.D.; Running, S.W. Young students, satellites aid understanding of climate-biosphere link. EOS 2000, 81, 1–5. [Google Scholar] [CrossRef]
  22. Zhang, X.; Friedl, M.A.; Schaaf, C.B.; Strahler, A.H.; Hodges, J.C.F.; Gao, F.; Reed, B.C.; Huete, A. Monitoring vegetation phenology using modis. Remote Sens. Environ. 2003, 84, 471–475. [Google Scholar] [CrossRef]
  23. Zhang, X.; Goldberg, M.D. Monitoring fall foliage coloration dynamics using time-series satellite data. Remote Sens. Environ. 2011, 115, 382–391. [Google Scholar] [CrossRef]
  24. Ma, X.; Huete, A.; Yu, Q.; Coupe, N.R.; Davies, K.; Broich, M.; Ratana, P.; Beringer, J.; Hutley, L.B.; Cleverly, J.; et al. Spatial patterns and temporal dynamics in savanna vegetation phenology across the north Australian tropical transect. Remote Sens. Environ. 2013, 139, 97–115. [Google Scholar] [CrossRef]
  25. Gitelson, A.A. Wide dynamic range vegetation index for remote quantification of biophysical characteristics of vegetation. J. Plant Physiol. 2004, 161, 165–173. [Google Scholar] [CrossRef] [PubMed]
  26. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  27. Zhao, B.; Yan, Y.; Guo, H.; He, M.; Gu, Y.; Li, B. Monitoring rapid vegetation succession in estuarine wetland using time series MODIS-based indicators: An application in the Yangtze river delta area. Ecol. Indic. 2009, 9, 346–356. [Google Scholar] [CrossRef]
  28. Atkinson, P.M.; Jeganathan, C.; Dash, J.; Atzberger, C. Inter-comparison of four models for smoothing satellite sensor time-series data to estimate vegetation phenology. Remote Sens. Environ. 2012, 123, 400–417. [Google Scholar] [CrossRef]
  29. Garrity, S.R.; Bohrer, G.; Maurer, K.D.; Mueller, K.L.; Vogel, C.S.; Curtis, P.S. A comparison of multiple phenology data sources for estimating seasonal transitions in deciduous forest carbon exchange. Agric. For. Meteorol. 2011, 151, 1741–1752. [Google Scholar] [CrossRef]
  30. Wang, Z.; Schaaf, C.B.; Sun, Q.; Kim, J.H.; Erb, A.M.; Gao, F.; Román, M.O.; Yang, Y.; Petroy, S.; Taylor, J.R.; et al. Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from landsat and the MODIS BRDF/NBAR/ALBEDO product. Int. J. Appl. Earth Obs. Geoinf. 2017, 59, 104–117. [Google Scholar] [CrossRef]
  31. Chen, J.M.; Pavlic, G.; Brown, L.; Cihlar, J.; Leblanc, S.G.; White, H.P.; Hall, R.J.; Peddle, D.R.; King, D.J.; Trofymow, J.A.; et al. Derivation and validation of Canada-wide coarse-resolution leaf area index maps using high-resolution satellite imagery and ground measurements. Remote Sens. Environ. 2002, 80, 165–184. [Google Scholar] [CrossRef]
  32. Fassnacht, K.S.; Gower, S.T.; Norman, J.M.; McMurtric, R.E. A comparison of optical and direct methods for estimating foliage surface area index in forests. Agric. For. Meteorol. 1994, 71, 183–207. [Google Scholar] [CrossRef]
  33. Justice, C.O.; Townshend, J.R.G.; Holben, B.N.; Tucker, C.J. Analysis of the phenology of global vegetation using meteorological satellite data. Int. J. Remote Sens. 1985, 6, 1271–1318. [Google Scholar] [CrossRef]
  34. Dall’Olmo, G.; Karnieli, A. Monitoring phenological cycles of desert ecosystems using NDVI and LST data derived from NOAA-AVHRR imagery. Int. J. Remote Sens. 2002, 23, 4055–4071. [Google Scholar] [CrossRef]
  35. Reed, B.C.; Brown, J.F.; VanderZee, D.; Loveland, T.R.; Merchant, J.W.; Ohlen, D.O. Measuring phenological variability from satellite imagery. J. Veg. Sci. 1994, 5, 703–714. [Google Scholar] [CrossRef]
  36. Jonsson, P.; Eklundh, L. Seasonality extraction by function fitting to time-series of satellite sensor data. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1824–1832. [Google Scholar] [CrossRef]
  37. Jönsson, P.; Eklundh, L. Timesat—A program for analyzing time-series of satellite sensor data. Comput. Geosci. 2004, 30, 833–845. [Google Scholar] [CrossRef]
  38. Fisher, J.I.; Richardson, A.D.; Mustard, J.F. Phenology model from surface meteorology does not capture satellite-based greenup estimations. Glob. Chang. Biol. 2007, 13, 707–721. [Google Scholar] [CrossRef]
  39. Song, C.; Ke, L.; You, S. Comparison of three NDVI time-series fitting methods based on timesat: Taking the grassland in northern Tibet as case. Remote Sens. Technol. Appl. 2011, 26, 147–155. [Google Scholar]
  40. Lu, L.; Guo, H. Extraction method of winter wheat phenology from time series of spot/vegetation data. Trans. Chin. Soc. Agric. Eng. 2009, 25, 174–179. [Google Scholar]
  41. Fisher, J.I.; Mustard, J.F.; Vadeboncoeur, M.A. Green leaf phenology at Landsat resolution: Scaling from the field to the satellite. Remote Sens. Environ. 2006, 100, 265–279. [Google Scholar] [CrossRef]
  42. Yu, F.; Price, K.P.; Ellis, J.; Shi, P. Response of seasonal vegetation development to climatic variations in eastern central Asia. Remote Sens. Environ. 2003, 87, 42–54. [Google Scholar] [CrossRef]
  43. Chen, X.; Tan, Z.; Schwartz, M.D.; Xu, C. Determining the growing season of land vegetation on the basis of plant phenology and satellite data in northern China. Int. J. Biometeorol. 2000, 44, 97–101. [Google Scholar] [CrossRef] [PubMed]
  44. Wen, G.; Fu, C.B. Large scale features of the seasonal phenological responses to the monsoon climate in east China: Multi-year average results. Sci. Atmos. Sin. 2000, 24, 676–682. [Google Scholar]
  45. Richardson, A.D.; Anderson, R.S.; Arain, M.A.; Barr, A.G.; Bohrer, G.; Chen, G.; Chen, J.M.; Ciais, P.; Davis, K.J.; Desai, A.R.; et al. Terrestrial biosphere models need better representation of vegetation phenology: Results from the north American carbon program site synthesis. Glob. Chang. Biol. 2012, 18, 566–584. [Google Scholar] [CrossRef] [Green Version]
  46. Morison, J.I.; Morecroft, M.D. Plant Growth and Climate Change; John Wiley & Sons: New York, NY, USA, 2008. [Google Scholar]
  47. Ganguly, S.; Friedl, M.A.; Tan, B.; Zhang, X.; Verma, M. Land surface phenology from MODIS: Characterization of the collection 5 global land cover dynamics product. Remote Sens. Environ. 2010, 114, 1805–1816. [Google Scholar] [CrossRef]
  48. Hmimina, G.; Dufrêne, E.; Pontailler, J.Y.; Delpierre, N.; Aubinet, M.; Caquet, B.; de Grandcourt, A.; Burban, B.; Flechard, C.; Granier, A.; et al. Evaluation of the potential of modis satellite data to predict vegetation phenology in different biomes: An investigation using ground-based NDVI measurements. Remote Sens. Environ. 2013, 132, 145–158. [Google Scholar] [CrossRef]
  49. Verger, A.; Filella, I.; Baret, F.; Peñtilde, J. Vegetation baseline phenology from kilometric global LAI satellite products. Remote Sens. Environ. 2016, 178, 1–14. [Google Scholar] [CrossRef]
  50. White, M.A.; De Beurs, K.M.; Didan, K.; Inouye, D.W.; Richardson, A.D.; Jensen, O.P.; O’Keefe, J.; Zhang, G.; Nemani, R.R.; Van Leeuwen, W.J.D.; et al. Intercomparison, interpretation, and assessment of spring phenology in north America estimated from remote sensing for 1982–2006. Glob. Chang. Biol. 2009, 15, 2335–2359. [Google Scholar] [CrossRef]
  51. Schwartz, M.D.; Hanes, J.M. Intercomparing multiple measures of the onset of spring in eastern north America. Int. J. Climatol. 2010, 30, 1614–1626. [Google Scholar] [CrossRef]
  52. White, K.; Pontius, J.; Schaberg, P. Remote sensing of spring phenology in northeastern forests: A comparison of methods, field metrics and sources of uncertainty. Remote Sens. Environ. 2014, 148, 97–107. [Google Scholar] [CrossRef]
  53. Xiao, Z.; Liang, S.; Wang, J.; Chen, P.; Yin, X.; Zhang, L.; Song, J. Use of general regression neural networks for generating the glass leaf area index product from time-series MODIS surface reflectance. IEEE Trans. Geosci. Remote Sens. 2014, 52, 209–223. [Google Scholar] [CrossRef]
  54. Melaas, E.K.; Richardson, A.D.; Friedl, M.A.; Dragoni, D.; Gough, C.M.; Herbst, M.; Montagnani, L.; Moors, E. Using fluxnet data to improve models of springtime vegetation activity onset in forest ecosystems. Agric. For. Meteorol. 2013, 171, 46–56. [Google Scholar] [CrossRef]
  55. Xiao, W.; Sun, Z.; Wang, Q.; Yang, Y. Evaluating MODIS phenology product for rotating croplands through ground observations. J. Appl. Remote Sens. 2013, 7, 073562. [Google Scholar] [CrossRef]
  56. Zhang, X.; Friedl, M.A.; Schaaf, C.B. Global vegetation phenology from moderate resolution imaging spectroradiometer (MODIS): Evaluation of global patterns and comparison with in situ measurements. J. Geophys. Res. Biogeosci. 2006, 111, G04017. [Google Scholar] [CrossRef]
  57. Xin, J.; Yu, Z.; van Leeuwen, L.; Driessen, P.M. Mapping crop key phenological stages in the north China plain using NOAA time series images. Int. J. Appl. Earth Obs. Geoinf. 2002, 4, 109–117. [Google Scholar] [CrossRef]
  58. Xiao, D.; Moiwo, J.; Tao, F.; Yang, Y.; Shen, Y.; Xu, Q.; Liu, J.; Zhang, H.; Liu, F. Spatiotemporal variability of winter wheat phenology in response to weather and climate variability in china. Mitig. Adapt. Strateg. Glob. Chang. 2015, 20, 1191–1202. [Google Scholar] [CrossRef]
Figure 1. The phenological metrics of vegetation dynamics.
Figure 1. The phenological metrics of vegetation dynamics.
Sensors 17 01982 g001
Figure 2. Land cover map and phenological observation sites distribution in China.
Figure 2. Land cover map and phenological observation sites distribution in China.
Sensors 17 01982 g002
Figure 3. Patterns of mean phenological metrics of the GLP and the MLCD. (a) SOS of GLP; (b) EOS of GLP; (c) SOS of MLCD; (d) EOS of MLCD.
Figure 3. Patterns of mean phenological metrics of the GLP and the MLCD. (a) SOS of GLP; (b) EOS of GLP; (c) SOS of MLCD; (d) EOS of MLCD.
Sensors 17 01982 g003
Figure 4. Missing ratios in different vegetation types.
Figure 4. Missing ratios in different vegetation types.
Sensors 17 01982 g004
Figure 5. Missing ratios in different latitudinal bands.
Figure 5. Missing ratios in different latitudinal bands.
Sensors 17 01982 g005
Figure 6. Proportion of vegetation in different latitudinal bands.
Figure 6. Proportion of vegetation in different latitudinal bands.
Sensors 17 01982 g006
Figure 7. The means and standard deviations of the difference values between the GLP and the MLCD in different vegetation types.
Figure 7. The means and standard deviations of the difference values between the GLP and the MLCD in different vegetation types.
Sensors 17 01982 g007
Figure 8. The difference values between the GLP and the MLCD in different latitudinal bands.
Figure 8. The difference values between the GLP and the MLCD in different latitudinal bands.
Sensors 17 01982 g008
Figure 9. Comparison between the LAI and the EVI time series for different vegetation types based on the ground observations. (a) HS (ENF); (b) DHS (EBF); (c) BJF (DNF); (d) CBS (DBF); (e) SPT (shrub); (f) NM (grass); (g) SY (rice); (h) JT (rice).
Figure 9. Comparison between the LAI and the EVI time series for different vegetation types based on the ground observations. (a) HS (ENF); (b) DHS (EBF); (c) BJF (DNF); (d) CBS (DBF); (e) SPT (shrub); (f) NM (grass); (g) SY (rice); (h) JT (rice).
Sensors 17 01982 g009aSensors 17 01982 g009b
Figure 10. Comparison between the GLP and the MLCD for different vegetation types based on the ground observations. (a) HS (ENF); (b) DHS (EBF); (c) BJF (DNF); (d) CBS (DBF); (e) SPT (Shrub); (f) NM (Grass); (g) SY + JT (rice).
Figure 10. Comparison between the GLP and the MLCD for different vegetation types based on the ground observations. (a) HS (ENF); (b) DHS (EBF); (c) BJF (DNF); (d) CBS (DBF); (e) SPT (Shrub); (f) NM (Grass); (g) SY + JT (rice).
Sensors 17 01982 g010
Table 1. Information for ground phenological observation sites in China.
Table 1. Information for ground phenological observation sites in China.
Station NameCodeVegetation TypeDominant SpeciesLonLatSourceYears
ShenyangSYcroprice123.36041.520CERN2004–2009
JiutaiJTcroprice125.80044.170CMA2003–2010
NaimanNMgrasshorsetail116.67643.550CERN2005–2010
ShapotouSPTshrubherbage105.00337.470CERN2002–2012
HeshanHSENFMasson’s pine, cedar112.90022.681CERN2004–2009
DinghushanDHSEBFCastanea henryi, Schima superba, Aporosa yunnanensis, Cryptocarya chinensis, Acmena acuminatissima112.53942.144CERN2004–2009
BeijingBJFDNFChinese pine, larch115.42539.958CERN2003–2011
ChangbaishanCBSDBFMeng gu oak128.10941.403CERN2003–2010
Table 2. One-to-one correspondence between remote-sensing phenological metrics and phenophases observed from ground observations.
Table 2. One-to-one correspondence between remote-sensing phenological metrics and phenophases observed from ground observations.
Vegetation TypeSOSPSPEEOS
evergreen treebud burst--leaf coloring
deciduous treebud burst--leaf defoliation
shrubbud burst--leaf defoliation
herbemergence--withering
riceregreeningheadinggrain-fillingharvest

Share and Cite

MDPI and ACS Style

Wang, C.; Li, J.; Liu, Q.; Zhong, B.; Wu, S.; Xia, C. Analysis of Differences in Phenology Extracted from the Enhanced Vegetation Index and the Leaf Area Index. Sensors 2017, 17, 1982. https://doi.org/10.3390/s17091982

AMA Style

Wang C, Li J, Liu Q, Zhong B, Wu S, Xia C. Analysis of Differences in Phenology Extracted from the Enhanced Vegetation Index and the Leaf Area Index. Sensors. 2017; 17(9):1982. https://doi.org/10.3390/s17091982

Chicago/Turabian Style

Wang, Cong, Jing Li, Qinhuo Liu, Bo Zhong, Shanlong Wu, and Chuanfu Xia. 2017. "Analysis of Differences in Phenology Extracted from the Enhanced Vegetation Index and the Leaf Area Index" Sensors 17, no. 9: 1982. https://doi.org/10.3390/s17091982

APA Style

Wang, C., Li, J., Liu, Q., Zhong, B., Wu, S., & Xia, C. (2017). Analysis of Differences in Phenology Extracted from the Enhanced Vegetation Index and the Leaf Area Index. Sensors, 17(9), 1982. https://doi.org/10.3390/s17091982

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop