Spatiotemporal Comparison and Validation of Three Global-Scale Fractional Vegetation Cover Products

Fractional vegetation cover (FVC) is an important parameter for many environmental and ecological models. Large-scale and long-term FVC products are critical for various applications. Currently, several global-scale FVC products have been generated with remote sensing data, such as VGT bioGEOphysical product Version 2 (GEOV2), PROBA-V bioGEOphysical product Version 3 (GEOV3) and Global LAnd Surface Satellite (GLASS) FVC products. However, studies comparing and validating these global-scale FVC products are rare. Therefore, in this study, the performances of three global-scale time series FVC products, including the GEOV2, GEOV3, and GLASS FVC products, are investigated to assess their spatial and temporal consistencies. Furthermore, reference FVC data generated from high-spatial-resolution data are used to directly evaluate the accuracy of these FVC products. The results show that these three FVC products achieve general agreements in terms of spatiotemporal consistencies over most regions. In addition, the GLASS and GEOV2 FVC products have reliable spatial and temporal completeness, whereas the GEOV3 FVC product contains much missing data over high-latitude regions, especially during wintertime. Furthermore, the GEOV3 FVC product presents higher FVC values than GEOV2 and GLASS FVC products over the equator. The main differences between the GEOV2 and GLASS FVC products occur over deciduous forests, for which the GLASS product presents slightly higher FVC values than the GEOV2 product during wintertime. Finally, temporal profiles of the GEOV2 and GLASS FVC products show better consistency than the GEOV3 FVC product, and the GLASS FVC product presents more reliable accuracy (R2 = 0.7878, RMSE = 0.1212) compared with the GEOV2 (R2 = 0.5798, RMSE = 0.1921) and GEOV3 (R2 = 0.7744, RMSE = 0.2224) FVC products over these reference FVC data.


Introduction
Fractional vegetation cover (FVC) is defined as the fraction of green vegetation seen from nadir, which can characterize the growth conditions and horizontal density of land surface live vegetation [1][2][3][4][5]. As a significant biophysical parameter involved in surface processes [6], the FVC is widely used for studies of the atmosphere, pedosphere, hydrosphere, ecology, and their interactions [7].
Thus, accurate and stable FVC products with regional and global scales are critical for related studies, such as those on climate change, numerical weather predictions, and land-surface processes [8,9].
Remote sensing technology is a feasible and reliable way for large-scale and long-term FVC generation due to its excellent ability to provide land surface observations. Currently, many FVC estimation algorithms have been developed based on remote sensing data, which can be divided into three major types: empirical methods, pixel unmixing models and machine learning methods [10,11]. Empirical methods build the statistical relationships between the FVC and vegetation indices or specific bands' reflectance through sufficient and reliable sample data. Empirical methods can achieve satisfactory accuracy at the regional scale with specific vegetation types. However, empirical methods become invalid over large-scale regions, in which the various vegetation types and land conditions increase uncertainties in the established relationships [10,12]. Pixel unmixing models assume that each pixel is composed of several components, and the fraction of vegetation composition is the corresponding FVC value of the pixel [6,12,13]. The dimidiate pixel model, as a widely used method of pixel unmixing models, assumes that each pixel can be divided into two parts: vegetation and non-vegetation [12,14]. However, the main limitation of pixel unmixing models is the determination of representative endmembers because of the complex land surface conditions and various spectral characteristics over a large scale [12,15,16]. Recently, machine learning methods have also been widely used to retrieve FVC values because of their computational efficiency and stable performances in nonlinear fitting [7]. Generally, machine learning methods estimate the FVC through training on a representative sample database containing pre-processed reflectance and corresponding simulated land surface parameters data [7]. Several algorithms of machine learning methods are proposed for FVC product generation over regional and global scales with satisfying results [1,10,12].
With these developed FVC estimation algorithms, some FVC products over regional and global scales are generated from remote sensing data, which are summarized in Table 1. Among these FVC products, CNES/POLDER data are generated based on empirical method, which establishes a statistical relationship between the FVC and the simulated reflectance data by means of a kernel-driven model [17,18]. However, the CNES/POLDER FVC product is only available for 1996, 1997, and 2003, which limits time-series designs and long-term earth-monitoring applications. The EUMETSAT/LSA SAF FVC is estimated through a probabilistic spectral mixture analysis method [6,19,20], which is provided with a 3-km spatial resolution and a daily temporal resolution [21]. Nevertheless, this product is only generated over Europe, Africa and South America, which limits its application for global-scale researches. For the Carbon Cycle and Change in Land Observational Products from an Ensemble of Satellites (CYCLOPES) FVC product, neural networks and PROSPECT + SAIL (PROSAIL) canopy radiative transfer model simulations are adopted for FVC estimation [22]. However, significant systematic underestimation was detected in the CYCLOPES FVC product by Verger [23,24]. The VGT bioGEOphysical product Version 1 (GEOV1) FVC algorithm is based on a neural network that treats the top-of-canopy normalized reflectances as input data and the corrected CYCLOPES FVC product as the expectant output values [24,25]. The GEOV1 FVC product is linearly correlated with the CYCLOPES FVC product [26], but the GEOV1 FVC product also has poor spatial continuity with a large amount of missing data over high-latitude regions [12]. In addition, the validation result indicates that the GEOV1 FVC product overestimates FVC values by up to 0.20 in croplands [12,27].
Considering the poor performance of the GEOV1 FVC product, a new version of the global FVC product, known as the VGT bioGEOphysical product Version 2 (GEOV2) FVC product, is proposed based on neural networks [28]. The GEOV2 FVC achieves better performance in terms of the spatiotemporal continuity and accuracy than the GEOV1 FVC product [29]. Additionally, to further improve the spatial resolution and maintain temporal continuity, the GEOV3 FVC product is produced using PROBA-V data with a 300 m spatial resolution, which adopts a neural network algorithm along with data filtering, gap filling and smoothing processes [30,31]. In addition, the Global LAnd Surface Satellite (GLASS) FVC product, under the support of China's National High Technology Research and Development Program, has comparable accuracy and better spatiotemporal continuity compared with the GEOV1 FVC product [7,12,32].
For these existing FVC products, comparisons and validations, especially over global-scale and long-term FVC products, are crucial for their usage and improvement in land surface models [24]. Currently, several studies have been performed to evaluate and compare the performance of these FVC products. For example, Mu et al validated the GEOV1 FVC product over croplands in the Heihe study area, and the results indicated that the GEOV1 FVC product presented systematic overestimation up to approximately 0.2 [27]. Ding et al compared and validated GEOV1 and regional Australian MODIS FVC products over the Australian continent. The spatial comparisons showed that there were robust correlations between these two FVC products. Temporally consistent validation showed that GEOV1 and Australian MODIS FVC products had similar seasonal variation and different magnitudes over several biome types, such as sparse vegetation and open broadleaved deciduous forest [24].Additionally, the SEVIRI and MERIS FVC products were compared over Europe, as well as Africa, by García-Haro et al, and clear differences between the two products were observed over South Africa, where MERIS generally presented lower FVC values than SEVIRI FVC product [33]. As a whole, the aforementioned studies of FVC comparisons and validations are confined to the regional scale or short periods. Studies comparing and validating global and long-term FVC products are rare. However, long time series of global FVC products are crucial in climate and hydrologic modelling, natural hazards monitoring and soil erosion risk assessment [5,9,34]. For various applications, it is important to assess and validate the differences among these FVC products [24]. Therefore, the main objective of this study is to conduct spatiotemporal consistency comparisons and validations for three representative global-scale FVC products, including the GEOV2, GEOV3, and GLASS FVC products. Moreover, the accuracy of these three FVC products is validated through reference FVC data, which are processed by the Implementing Multi-scale Agricultural Indicators Exploiting Sentinels (IMAGINES) project.

GEOV2 FVC Product
The GEOV2 FVC product was derived from SPOT/VEGETATION data from January 1999 to December 2013 and PROBA-V data from January 2014 by the Copernicus Global Land Service. The main purpose of GEOV2 FVC product generation was improving the GEOV1 FVC product in terms of accuracy and spatiotemporal continuity, especially at high-latitude and equatorial areas [28]. For GEOV2 FVC product generation, a neural network algorithm was adopted to obtain instantaneous FVC estimations [29,36,37]. Training samples were collected from the BELMANIP2 sites from VEGETATION data and the corrected CYCLOPES FVC product over the 2003-2007 period [29,38]. After correcting for atmospheric effects by means of the SMAC model [39], the top-of-canopy daily reflectance data with three bands (red, NIR and SWIR) were acquired as input data for FVC estimation [29]. Then, multi-step outlier rejection processes were conducted to remove abnormal estimations, which might be contaminated by atmospheric effects, such as clouds, aerosols, water vapor and ozone [40,41]. Finally, dedicated temporal techniques, combining an adaptive Savitzky-Golay (SG) filter and an innovative climatology fitting method (Consistent Adjustment of the Climatology to Actual Observations, CACAO), were applied to ensure continuity and robustness of the estimated FVC data [28,40,41]. In addition, the GEOV2 FVC product was provided with 10-day temporal resolution and 1-km spatial resolution by the Copernicus Global Land Service (https://land.copernicus.eu/global/ products/fcover). In this study, the GEOV2 FVC product from 2001 to 2016 was adopted for comparison and validation.

GEOV3 FVC Product
To maintain continuity and improve the spatial resolution of the GEOV2 FVC product, a new algorithm was employed to develop the GEOV3 FVC product using PROBA-V observation data under the IMAGINES project [30]. Similar to the GEOV2 FVC algorithm, a neural network algorithm was used for GEOV3 FVC estimation with the blue, red and NIR reflectance as input data. Because there was no available PROBA-V data when this algorithm was proposed, the training process was conducted based on SPOT/VEGETATION data, whose bands were close to PROBA-V data in the VIS-NIR channels [28,30]. In addition, the CYCLOPES FVC product was adopted as the reference FVC values for training samples in GEOV3 FVC algorithm development. Because of the systematic underestimation of the CYCLOPES FVC product, it had a maximum value of 0.6872 (FVC CYCV31 (99%) = 0.6872), which corresponds to full vegetation coverage FVC = 1. To address this problem, a scaling factor was applied to the CYCLOPES FVC product to correct the underestimation (formula 1) [25,28,30]. Compared with GEOV2 FVC product generation, the main discrepancy for the GEOV3 FVC product was that the SWIR reflectance were abandoned as input data. Because SWIR band (700 m) had approximate twice spatial resolution than VIS-NIR bands (300-m) in PROBA-V data [30,42]. The GEOV3 FVC product from 2014 to 2016 was obtained from the Copernicus Global Land Service (https://land.copernicus.eu/global/products/fcover) and was used for comparison and validation in this study.
where FVC CYCV31 is the extracted FVC value from the CYCLOPES V3.1 product, and FVC training is the corresponding corrected FVC value used to train the neural network.

GLASS FVC Product
The GLASS FVC product used in this study was generated by Beijing Normal University, China [12,43,44]. In 2015, Jia et al proposed the GLASS FVC estimation algorithm based on general regression neural networks (GRNNs) for MODIS surface reflectance data with red and NIR bands, and this algorithm achieved comparable accuracy and better performance in terms of the spatial continuity compared with the GEOV1 FVC product [12]. However, the developed GRNN-based FVC estimation algorithm presented unsatisfactory computation efficiency: it usually took over one hour to generate FVC data for one MODIS tile. Therefore, in 2016, Yang et al proposed the multivariate adaptive regression splines (MARS) method for GLASS FVC generation that presented both adequate computational efficiency and comparable accuracy to the GRNN method [7]. Under the direct validation in an agricultural region, the GLASS FVC product was more accurate (R 2 = 0.86, RMSE = 0.087) than the GEOV1 FVC product (R 2 = 0.71, RMSE = 0.193) [27,45]. Moreover, the FVC time series data used for validation were consistent with the whole-crop growing characteristics [45]. The GLASS FVC product had an 8-day temporal resolution, a 500-m spatial resolution, and sinusoidal grid projection, which was released in 'HDF' format by the National Earth System Science Data Sharing Infrastructure, National Science & Technology Infrastructure of China (http://www.geodata.cn/thematicView/GLASS.html). In this study, the GLASS FVC product from 2001 to 2016 was used for comparison and validation.

Field Survey Based Reference Data and MODIS Land Cover Data
To evaluate the performances of different FVC products, in-situ FVC measurements were acquired under the IMAGINES programme, which aimed to continue the innovation and development activities, as well as support the operations of the Copernicus Global Land service (http://fp7-imagines.eu/pages/ services-and-products/ground-data.php). With this purpose, the ground FVC measurement data were collected and processed by the Earth Observation LABoratory (EOLAB) and/or local teams [46,47]. To obtain the reference data, digital hemispherical photographs were taken and processed for the ground FVC data over different vegetation types. Then, an empirical transfer function between the reflectance and ground FVC values was derived based on multiple robust regression for each site. Finally, reference FVC maps were produced using high-spatial-resolution data, such as Landsat-8 and FASat-C data. For kilometric biophysical product assessment, these high-spatial-resolution reference maps were averaged over an area of 3 km×3 km. Because of the accurate generation processes and containing various vegetation types around the world, the reference data had reliable accuracy and were representative to evaluate the performance of the FVC products. Considering the temporal coverage of the adopted FVC products, 29 reference samples were available during 2014-2016 (showed in Table 2) and adopted to perform the evaluation in this study. In addition, the RMSE values of each site were calculated using the ground FVC data and their corresponding estimates from transfer function, which indicated the uncertainty of these reference FVC data. Moreover, the MODIS Land Cover Climate Modelling Grid Product (MCD12C1) was also used as the base map to evaluate the performance of the three FVC products over different vegetation types (https://e4ftl01.cr.usgs.gov/). MCD12C1 was provided at a 0.05 • spatial resolution and an annual temporal resolution, and was aggregated and re-projected from the MODIS Land Cover Type Product (MCD12Q1, 500-m spatial resolution, 8-day temporal resolution) [48,49]. The International Geosphere-Biosphere Programme (IGBP) classification system was adopted as the basic classes for further comparison and validation in this study.

Methodology
In this study, three global scale FVC products, including GEOV2, GEOV3, and GLASS FVC products, were selected for comparison and validation. According to [24] and [50], evaluation and comparison of different land surface parameter products were performed through three main aspects: spatial consistency, temporal consistency and accuracy validation, which drawn reasonable conclusions about the compared products. Thus, the designed methods in this study were mainly referred to these two studies. Because of the multiple spatiotemporal resolutions of these FVC products, it was inconvenient to conduct comparisons and analyses directly. Therefore, all three FVC products were first resampled to a 0.05 • spatial resolution using bilinear interpolation, which has a trade-off between interpolation accuracy and computational efficiency [51]. Additionally, monthly averaged FVC data were also calculated using the resampled FVC products for further comparison.
To conduct the comparison in terms of spatial consistency, the monthly averaged FVC maps of the GEOV2, GEOV3, and GLASS FVC products in January and July of 2015, which showed the representative spatial patterns in winter and summer, were selected to present their spatial patterns and completeness. In addition, the difference maps between different FVC products were calculated to investigate their spatial discrepancies [24,50]. Furthermore, histograms of these pre-processed monthly FVC data from 2014 to 2016 over different vegetation types were generated according to the MODIS land cover product (MCD12C1, IGBP), which aimed to further explain the distribution patterns of the FVC values for different vegetation types [48].
In the temporal consistency comparison, temporal profiles of the mean values for different land cover types were calculated to determine the temporal consistency among these FVC products [24,51]. Moreover, several temporal profiles were extracted from the processed FVC data over the site locations in Table 2. Inter-annual and seasonal variations of these FVC products were analysed through these temporal profiles across different land conditions, as well as vegetation types. In addition, the reference FVC values were also added to the extracted temporal profiles for comparison. Furthermore, the accuracies of these three products were evaluated by direct comparison with the reference FVC data. In this study, these FVC products with original spatiotemporal resolutions were used for accuracy validation. Because the reference data were the averaged FVC values over a 3 km × 3 km region at each site, the mean values of 5 × 5 pixel subset at each site were adopted for the GLASS FVC product, the 3×3 pixel subset for the GEOV2 FVC product and the 9×9 pixel subset for the GEOV3 FVC product. The R-squared (R 2 ) value and root mean square error (RMSE) were used to quantify the uncertainty and accuracy of each FVC product.
where n denotes the number of selected FVC values, which is 29 in this study; FVC Product(i) denotes the ith FVC value extracted from the sites in Table 2; FVC Re f erence(i) denotes the ith FVC value of the reference FVC data in Table 2; and FVC Product denotes the mean value of the extracted FVC data from the corresponding product.

Spatial Consistency
To compare the spatial consistency, the monthly FVC maps of the GEOV2, GEOV3, and GLASS FVC products in January and July of 2015 are shown in Figure 1. The dark grey areas represent missing data. Visually, the GLASS FVC products present excellent spatial completeness and no missing data were found, which is mainly attributed to the pre-processing procedure to obtain continuous and smooth reflectance data before the product generation process [12]. Additionally, the GEOV2 FVC product shows a significant improvement in the spatial completeness over the GEOV1 FVC product, the latter presenting approximately 20% missing data over the BELMANIP2 sites during 1999-2010 [52,53]. However, the GEOV3 FVC product presents poor spatial completeness because of a large fraction of missing values in the high-latitudes are found during winter. Figure 1 also shows a general agreement in the spatial distribution among these three FVC products. In January, high FVC values are mainly concentrated around the equator areas, such as north South America, Central Africa and Indonesia, where tropical and subtropical moist broadleaf forests are the dominant vegetation types [54]. In addition, North America, Asia and Europe have low FVC values. During summer, there is a clear increase in the FVC values in northern and eastern Asia, North America and Europe. Moreover, all three products show low FVC values during both winter and summer over the midwestern Qinghai-Tibet Plateau, the Himalayas and the Arabian Peninsula, which are mainly covered by desert scrublands or montane grasslands [54]. These spatial distributions of these FVC products are in accordance with the global terrestrial ecoregions and seasonal variations.    To quantify the differences in these spatial patterns, the difference maps for these FVC products in January and July of 2015 are generated in Figure 2. Considering the poor spatial completeness of the GEOV3 FVC product, only pixels where both selected FVC products have valid FVC data are used to calculate the difference maps. With the negative values around the equator in Figure 2a-d, both the GEOV2 and GLASS FVC products show lower FVC values than the GEOV3 FVC product. For example, the GEOV3 FVC product presents approximately 0.15 and 0.13 larger than the GEOV2 and GLASS FVC values over the Amazon rain forest region. Generally, GLASS and GEOV2 show clear discrepancies across different hemispheres in January. Specifically, the GLASS FVC values are slightly higher than the GEOV2 FVC values in the Northern Hemisphere but lower in the Southern Hemisphere. In July, negative difference values between the GLASS and GEOV2 FVC products are found in Europe, Southwest Asia and southern North America. In addition, random signs for the difference values are observed over Europe between the GLASS and GEOV3 FVC products in July 2015.   To further evaluate the spatial consistency of these three FVC products, the histograms of the GEOV2, GEOV3, and GLASS FVC values during 2014-2016 with different biome types are generated based on the MODIS land cover data and pre-processed monthly FVC products (Figure 3). Because there are some missing data in the GEOV3 product, only the pixels with available FVC values of all three products in the same period are adopted for statistical analysis. Additionally, several vegetation types, which are mainly distributed in high-latitudes, are also abandoned for further evaluation to avoid the disturbance of the missing data in GEOV3 FVC product (like Figure A1 in Appendix A), and finally the classes, including evergreen broadleaf forest, deciduous needleleaf forest, deciduous broadleaf forests, closed shrubland, grasslands, permanent wetlands, and croplands, as well as barren or sparsely vegetated, are used for further analysis. Generally, histograms of the three FVC products show excellent agreements over the adopted vegetation types, particularly in closed shrubland and croplands, where all three FVC products have similar distribution patterns. The spatial distributions for evergreen broadleaf forest of the GLASS and GEOV2 FVC products are more consistent than that of the GEOV3 FVC product. There were more pixels, which are concentrated between 0.6 and 0.95, in the GLASS and GEOV2 FVC products. In contrast, the GEOV3 FVC product has more pixels with higher values close to 1. For deciduous needleleaf forest, the GEOV2 and GEOV3 FVC products show better consistency with similar frequency distributions. However, the GLASS FVC product has more pixels between 0.65 and 0.9 over deciduous needleleaf forest. For permanent wetlands, the frequencies of the GLASS FVC values between 0.3 and 0.6 are higher than the GEOV2 and GEOV3 values. For closed shrubland and croplands, all three products achieve excellent agreements in the frequency distributions. For grasslands, the GLASS FVC product showed more pixels with low values close to 0. For barren or sparsely vegetated regions, most of the FVC values of these products were less than 0.1.    Figure 4. By visual observation, the varying magnitude of all the FVC products are coincident with the corresponding vegetation types. For evergreen broadleaf forest, the GLASS and GEOV2 FVC values have stable consistency, showing slight seasonality and maintaining high FVC values throughout the year. However, overall, the GEOV3 FVC values are systematically 0.05 higher than the GLASS and GEOV2 FVC values for evergreen broadleaf forests. For deciduous needleleaf and broadleaf forests, the main discrepancies between GLASS and GEOV2 FVC products occur during the seasons with low FVC values, when the GEOV2 presents clearly smaller FVC values than the GLASS product at valley values. The GEOV3 FVC product shows more agreements with the GEOV2 FVC product over deciduous needleleaf forests. Nevertheless, there is better consistency between the GEOV3 and GLASS products in deciduous broadleaf forests. For closed shrubland and croplands, all three FVC products achieve similar temporal variations and seasonal dynamics. For grasslands, the GEOV3 product presents smaller seasonal dynamics and higher FVC values than both the GEOV2 and GLASS FVC products, which present stable temporal consistency. For barren or sparsely vegetated biomes, all of these products have low inter-annual variations, as expected, and the GEOV2 values were slightly higher than the GLASS values.   Figure 5 shows the temporal profiles derived from the GEOV2, GEOV3, and GLASS FVC products over the site locations from Table 2. For better visual presentation and comparison, the reference FVC values of the selected sites are exhibited in Figure 5. In Figure 5a,d, both the GEOV2 and GEOV3 products have similar seasonal dynamics and magnitudes, whereas the GLASS FVC  Figure 5 shows the temporal profiles derived from the GEOV2, GEOV3, and GLASS FVC products over the site locations from Table 2. For better visual presentation and comparison, the reference FVC values of the selected sites are exhibited in Figure 5. In Figure 5a,d, both the GEOV2 and GEOV3 products have similar seasonal dynamics and magnitudes, whereas the GLASS FVC values are smaller than those of GEOV2 and GEOV3 for the rice, tree plantation, tea and coffee vegetation types during the growing season. In Figure 5b,f, the GEOV2 and GLASS FVC products achieve reliable agreement in wheat, barley, corn, sunflower and soybean croplands. However, the GEOV3 values are over 0.15 higher than the GEOV2 and GLASS FVC values during summer. Furthermore, all these FVC products demonstrate satisfactory temporal consistency in Figure 5c,e, which are covered by wheat, corn soybean and grass. Generally, all these time series of different FVC products present consistent seasonal dynamics and magnitude ranges with corresponding vegetation types. Obviously, the GEOV2 and GLASS FVC products show complete temporal continuity, whereas the temporal profiles of the GEOV3 FVC product are discontinuous with poor-quality temporal continuity.

Accuracy Validation
The scatterplots of the three FVC products and the reference FVC values are shown in Figure 6. The RMSE of the reference FVC data is taken as the error bar in Figure 6, which aims to visually indict the uncertainty of reference data. Obviously, the GLASS FVC product has better consistency with the reference data compared with the GEOV2 and GEOV3 FVC data. A regression equation with a 0.7452 slope and a 0.2680 intercept is established between the GEOV2 FVC product and reference FVC values in Figure 6a, which presents high uncertainty and overestimation of small estimations based on the available reference data. In Figure 6b, a clear systematic overestimation (up to approximately 0.2) of the reference data is observed for the GEOV3 FVC product, whose regression equation has a slope of 1.0162 and an intercept of 0.1887 with the reference data. Furthermore, the GLASS FVC product achieves a more reliable accuracy (R 2 = 0.7878, RMSE = 0.1212) than the GEOV2 (R 2 = 0.5798, RMSE = 0.1921) and GEOV3 FVC data (R 2 = 0.7744, RMSE = 0.2224) over the reference FVC values. In addition, considering the uncertainty of the reference FVC values in each field site, there are 8, 3 and 16 sites within the error range of the GEOV2, GEOV3, and GLASS FVC products, respectively. This finding further indicates the better performance of the GLASS FVC product with the available

Accuracy Validation
The scatterplots of the three FVC products and the reference FVC values are shown in Figure 6. The RMSE of the reference FVC data is taken as the error bar in Figure 6, which aims to visually indict the uncertainty of reference data. Obviously, the GLASS FVC product has better consistency with the reference data compared with the GEOV2 and GEOV3 FVC data. A regression equation with a 0.7452 slope and a 0.2680 intercept is established between the GEOV2 FVC product and reference FVC values in Figure 6a, which presents high uncertainty and overestimation of small estimations based on the available reference data. In Figure 6b, a clear systematic overestimation (up to approximately 0.2) of the reference data is observed for the GEOV3 FVC product, whose regression equation has a slope of 1.0162 and an intercept of 0.1887 with the reference data. Furthermore, the GLASS FVC product achieves a more reliable accuracy (R 2 = 0.7878, RMSE = 0.1212) than the GEOV2 (R 2 = 0.5798, RMSE = 0.1921) and GEOV3 FVC data (R 2 = 0.7744, RMSE = 0.2224) over the reference FVC values. In addition, considering the uncertainty of the reference FVC values in each field site, there are 8, 3 and 16 sites within the error range of the GEOV2, GEOV3, and GLASS FVC products, respectively. This finding further indicates the better performance of the GLASS FVC product with the available reference FVC data.
14 of 19 Figure 6. Scatterplots between different FVC products and high-spatial resolution reference FVC data.

Discussion
The comparison and validation for the GEOV2, GEOV3, and GLASS FVC products show that these three global-scale products present satisfying agreements, which have very similar spatiotemporal variations over most areas on a global scale. During the GLASS FVC product generation, a MODIS reflectance re-processing method was developed to obtain the continuous and smooth surface reflectance data by Tang et al [55]. In this method, the contaminated reflectance data were identified through the MODIS quality control data, the temporal characteristics of the spectrum, as well as other auxiliary information, and missing data were filled using an optimum interpolation algorithm [12,55]. With reliable reflectance data, the GLASS FVC product achieved excellent spatial, as well as temporal continuity and completeness. For the GEOV2 FVC product, compositing processes were performed for the FVC estimations instead of the reflectance data, which aimed to reduce the sensitivity of the missing data. Consequently, this strategy tremendously improved the spatial continuity and completeness compared to the GEOV1 FVC product [1,28,53]. For the GEOV3 FVC product, the main reason for the missing data over the high latitudes in winter was the snow coverage and persistent cloudiness, which led to high uncertainty in the PROBA-V reflectance data [56]. Moreover, during the GEOV3 FVC generation, a particular gap-filling method was used over evergreen broadleaf forests because of the large amount of cloud contaminated observations, which might be the main explanation for the overestimation of GEOV3 FVC in equator regions [56,57].
Due to the reprocessing and smoothing processes, both GEOV2 and GLASS FVC products have more agreements and reliable seasonal variations over the most temporal profiles. However, the temporal profiles of the GEOV3 FVC data are slightly noisy in deciduous forests, croplands and grasslands. Because there are no available climatology data, it is difficult to build with limited time series data on a 300-m pixel scale [30]. For accuracy validation, both the GEOV2 and GEOV3 products adopt the corrected CYCLOPES FVC values as training samples. Nevertheless, this procedure may lead to a limitation of accuracy when CYCLOPES FVC do not perform well in some cases [56,57]. Hence, both the GEOV2 and GEOV3 FVC algorithms may need to fine-tune the parameters for highaccuracy estimations. Specifically, real PROBA-V data may be required for neural network training in GEOV3 FVC generation. For the GLASS FVC product, high quality Landsat TM/ETM+ data are used to generate the high-spatial-resolution FVC values as training samples [12]. The sample data are further refined based on the approximate linear relationship between the FVC and NDVI values. With the high-quality sample data, the GLASS FVC product is more reliably accurate in this study.

Conclusion
In this study, the performances of the three representative and available global-scale FVC time series products, including the GEOV2, GEOV3, and GLASS FVC products, were compared and Figure 6. Scatterplots between different FVC products and high-spatial resolution reference FVC data.

Discussion
The comparison and validation for the GEOV2, GEOV3, and GLASS FVC products show that these three global-scale products present satisfying agreements, which have very similar spatiotemporal variations over most areas on a global scale. During the GLASS FVC product generation, a MODIS reflectance re-processing method was developed to obtain the continuous and smooth surface reflectance data by Tang et al [55]. In this method, the contaminated reflectance data were identified through the MODIS quality control data, the temporal characteristics of the spectrum, as well as other auxiliary information, and missing data were filled using an optimum interpolation algorithm [12,55]. With reliable reflectance data, the GLASS FVC product achieved excellent spatial, as well as temporal continuity and completeness. For the GEOV2 FVC product, compositing processes were performed for the FVC estimations instead of the reflectance data, which aimed to reduce the sensitivity of the missing data. Consequently, this strategy tremendously improved the spatial continuity and completeness compared to the GEOV1 FVC product [1,28,53]. For the GEOV3 FVC product, the main reason for the missing data over the high latitudes in winter was the snow coverage and persistent cloudiness, which led to high uncertainty in the PROBA-V reflectance data [56]. Moreover, during the GEOV3 FVC generation, a particular gap-filling method was used over evergreen broadleaf forests because of the large amount of cloud contaminated observations, which might be the main explanation for the overestimation of GEOV3 FVC in equator regions [56,57].
Due to the reprocessing and smoothing processes, both GEOV2 and GLASS FVC products have more agreements and reliable seasonal variations over the most temporal profiles. However, the temporal profiles of the GEOV3 FVC data are slightly noisy in deciduous forests, croplands and grasslands. Because there are no available climatology data, it is difficult to build with limited time series data on a 300-m pixel scale [30]. For accuracy validation, both the GEOV2 and GEOV3 products adopt the corrected CYCLOPES FVC values as training samples. Nevertheless, this procedure may lead to a limitation of accuracy when CYCLOPES FVC do not perform well in some cases [56,57]. Hence, both the GEOV2 and GEOV3 FVC algorithms may need to fine-tune the parameters for high-accuracy estimations. Specifically, real PROBA-V data may be required for neural network training in GEOV3 FVC generation. For the GLASS FVC product, high quality Landsat TM/ETM+ data are used to generate the high-spatial-resolution FVC values as training samples [12]. The sample data are further refined based on the approximate linear relationship between the FVC and NDVI values. With the high-quality sample data, the GLASS FVC product is more reliably accurate in this study.

Conclusions
In this study, the performances of the three representative and available global-scale FVC time series products, including the GEOV2, GEOV3, and GLASS FVC products, were compared and evaluated for spatiotemporal consistency and accuracy. Based on this study, the following conclusions can be drawn: (1) The GLASS FVC product is much more continuous and complete based on the comparison. Additionally, the GLASS FVC product presents more reliable accuracy by direct comparison against the available reference FVC data (R 2 = 0.7878, RMSE = 0.1212). (2) The GEOV2 FVC product shows numerous consistencies with the GLASS FVC product over most regions on a global scale. According to the accuracy validation using the available reference data, high uncertainties are found for low FVC values in the GEOV3 FVC product (R 2 = 0.5798, RMSE = 0.1921). (3) The GEOV3 FVC product shows poor performance in the spatiotemporal completeness and continuity. Moreover, systematic overestimation (up to approximately 0.2) is observed in the GEOV3 FVC product compared with the reference FVC values (R 2 = 0.7744, RMSE = 0.2224). Additionally, the number of the reference data is small in this study and our further work will focus on extensive validation of these FVC products using reference data with more vegetation types and land surface conditions. Author Contributions: K.J. conceptualized this study; D.L and K.J. made the data curation; D.L., K.J. and X.W. provided the methods of this study; K.J offered the resources; K.J. supervised this study; D.L performed the processes of comparison and validation for adopted FVC products; D.L wrote this paper; K.J., X.W., X.Z., M.X., Y.Y., X.Z. and B.W. reviewed and edited this paper.   Figure A1 showed the frequency of missing value over evergreen needle forest in GEOV3 FVC product during 2014 to 2016. Obviously, there were approximately 50% of missing data within half a year over evergreen needle forest in GEOV3 FVC product. Considering the massive missing data of GEOV3 FVC product over evergreen needle forests, it was reasonable to drop this vegetation type for further evaluation.