Large Uncertainty on Forest Area Change in the Early 21st Century among Widely Used Global Land Cover Datasets

: Forests play an important role in the Earth’s system. Understanding the states and changes in global forests is vital for ecological assessments and forest policy guidance. However, there is no consensus on how global forests have changed based on current datasets. In this study, ﬁve global land cover datasets and Global Forest Resources Assessments (FRA) were assessed to reveal uncertainties in the global forest changes in the early 21 st century. These datasets displayed substantial divergences in total area, spatial distribution, latitudinal proﬁle, and annual area change from 2001 to 2012. These datasets also display completely divergent conclusions on forest area changes for di ﬀ erent countries. Among the datasets, total forest area changes range from an increase of 1.7 × 10 6 km 2 to a decrease of 1.6 × 10 6 km 2 . All the datasets show deforestation in the tropics. The accuracies of the datasets in detecting forest cover changes were evaluated by a global land cover validation dataset. The spatial patterns of accuracies are inconsistent among the datasets. This study calls for the development of a more accurate database to support forest policies and to contribute to global actions against climate change.


Introduction
As one of the most widely distributed land cover types, forests have a total area of 40.6 × 10 6 km 2 , which accounts for 31% of the world's total land area [1]. Forests play an important role in maintaining the balance of the global ecosystem by acting as a carbon sink [2,3], prompting water conservation [4], providing a habitat for species [5], improving landscape functions [6], and regulating the climate [7]. Forests have experienced great changes, including deforestation and afforestation. Understanding states of global forest and its changes help to provide guidelines for forest conservation, protection, and management [8,9].
According to the Food and Agriculture Organization (FAO) of the United Nations (UN), more than 1.29 × 10 6 km 2 of forest have been lost globally since 1990 [1]. Fires, insects, severe weather events, and ecological evolution are important natural factors that cause changes in forest distribution [10,11]. Human activities have significant influence on forests. This includes sourcing wood for industrial and were compared. The global land cover validation dataset from the United States Geological Survey (USGS) was utilized to evaluate the accuracy of the land cover datasets.

VCF
The Making Earth System Data Records for Use in Research Environments (the MEaSURES) VCF products use a non-parametric trend analysis in each pixel to measure land cover changes [27]. VCF provides the annual global fraction of tree canopy (TC) cover, short vegetation (SV) cover, and bare ground (BG) cover from 1982 to 2016. This dataset is based on the AVHRR and other supplementary data, including Landsat Enhanced Thematic Mapper Plus (ETM+), the Moderate Resolution Imaging Spectroradiometer (MODIS), and other very high-resolution satellite images [27]. This dataset can be obtained from Global Land Analysis and Discovery [35]. In this study, TC was selected to calculate forest area changes from 2001 to 2012 (Table 1).

Hansen
Hansen is the high-resolution global maps of twenty-first-century forest cover, including tree canopy data for the year 2000, global forest cover gain between 2000 and 2012, year of gross forest cover loss event, and data mask at 30 m × 30 m resolution [28,36]. Forest cover gain is displayed as a binary value that provides information on whether forest gain occurred in each pixel from 2000 to 2012.
Hansen also provides information of the years when net forest loss occurred in each pixel. Hansen was produced based on 654,178 growing season Landsat 7 ETM+ images. We used the forest cover gain for 2001-2012, forest cover loss for 2001-2017, and land mask (whether the specific pixel represents land or not) to calculate the global forest area change from 2001 to 2012 (Table 1) [28].

MCD12C1
Version 6 of MCD12C1 product provides land cover data using three land cover classification schemes, i.e., International Geosphere-Biosphere Programme (IGBP), University of Maryland (UMD), and Leaf Area Index (LAI). For each land cover classification scheme, MCD12C1 provides global land cover fraction at a 0.05 × 0.05 • spatial resolution from 2001 to 2018. In this study, the IGBP classification scheme was utilized to calculate global forest area changes. Forest cover is defined by the pixels of evergreen broadleaf forest, deciduous needle leaf forest, or deciduous broadleaf forest layers [37,38].

CCI-LC
CCI-LC product is produced by the Land Cover (LC) project of the Climate Change Initiative (CCI), which is led by the Europe Space Agency (ESA) [39]. Three original data types were used to develop CCI-LC, including AVHRR data from 1992 to 1999, Systeme Probatoire d'Observation de la Terre-Vegetation (SPOT-VGT) data from 1999 to 2013, and Project for On-Board Autonomy Vegetation (PROBA-V) data for 2013, 2014, and 2015. CCI-LC is at a 300 m × 300 m resolution, based on the Glob Cover unsupervised classification chain, which utilizes a machine learning algorithm [40]. CCI-LC uses the hierarchical classification system of the United Nations Land Cover Classification System (UN-LCCS), which was developed by the FAO.

LUH2
LUH2 is the new generation of global land-use forcing datasets that provide a standard format of historical land use and future projections for climate models. The dataset provides historical reconstructions of land-use, which can be used with future projections via Earth system models. 2.1.6. FRA Reports FRA reports are the most comprehensive forest assessment datasets. The dataset provides the global forest resources and their changes every five year. FRA reports are based on country level government inventories and remote sensing data. FRA reports provided every five-year assessment about the global forest resources and their changes. FRA reports were widely used for forest conditions popularizations, policy guidance, and land cover data accuracy validations. FRA report 2000, FRA report 2010, FRA report 2015, and FRA report 2020 were selected to compare in the forest area change with the global land cover datasets at country level.

USGS Global Land Cover Validation Data
Global land cover validation data were developed by USGS based on high resolution commercial satellite imagery. Over 1100 individual scenes of commercial data imagery were utilized to produce the validation data from several satellites, including QuickBird, WorldView-1, WorldView-2, IKONOS, OrbView-3, and GeoEye-1. The USGS validation dataset, includes 500 5 km × 5 km samples at a global scale, with samples locations that were randomly selected based on the modified Köppen climate zones and population [43,44]. The validation dataset provides the dominant land cover type in each pixel at a 2 m spatial resolution, spanning 2001-2014.

Data Processing
To ensure that the global land cover datasets at different temporal and spatial scales (VCF, Hansen, MCD12C1, CCI-LC, and LUH2) are comparable, all of the land cover datasets were aggregated to the same temporal period, spatial resolution, and spatial coverage. As Hansen exhibits the least latitude coverage (80 • N-60 • S) and shortest temporal period (2001-2012) among the five land cover datasets, all of the datasets were adjusted to the same spatiotemporal extent (latitude: 80 • N-60 • S; time period: 2001-2012). We also processed all the five datasets at a 0.05 • resolution following the procedures outlined below.
Hansen provides information on when forest loss occurred from 2001 to 2017 and whether forest gain occurred between 2001 and 2012 in each 30 m × 30 m pixel. The number of pixels of forest loss, forest gain, and the land pixel in every 200 × 200 Hansen block (matched to a 0.05 × 0.05 • pixel) were calculated. The forest change ratio in a 0.05 × 0.05 • pixel was calculated by dividing the number of net change pixels (the number of forest gain pixels minus the number of forest loss pixels) by the total number of land pixels: where i and j refers to the row number and column number of the 0.05 × 0.05 • pixel; num gain refers to the number of forest gain in the pixel [i,j]; num loss refers to the number of forest loss event range from 1 to 12 in the pixel; num land refers to the number of land surface in the pixel. CCI-LC provides the land cover types in each 250 m × 250 m pixel. In this study, regions with the class code 30, 40, 50, 60, 70, 80, 90, and 100 were defined as forest areas based on the legend system. We counted the pixel number of the above layers as forest cover in every 18 × 18 CCI-LC block (matched to a 0.05 × 0.05 • pixel), representing the forest cover ratio in each 0.05 • pixel.
Forest cover ratio of VCF was obtained from the Percent Tree Cover layer. For the MCD12C1, evergreen needle leaf forest, evergreen broadleaf forest, deciduous needle leaf forest, and deciduous broadleaf forest land cover types were selected as forest cover under Land_Cover_Type_1_Percent for the IGBP scheme. For LUH2, the variables of primf (forested primary land) and secdf (potentially forested secondary land) in the states layer were selected to obtain the percentage of forest cover in each pixel. For these three datasets, forest area change in each pixel was calculated by multiplying the forest cover percentage by the area of the pixel.
To assess the accuracy of the land cover datasets, the forest cover percentage of each sample was calculated. First, we transformed the projection of the global land cover validation data from Universal Transverse Mercator (UTM) Projection to latitude and longitude. Second, the corresponding row number and column number for each sample were calculated based on the latitude and longitude. Then the numbers of forest pixel in each sample were calculated. The forest cover percentage in the sample is equal to the pixel numbers of forest divided by the total pixel numbers in the sample. We only found date information for 300 samples out of the 500 samples, such that the accuracy validation was based on the 300 samples with date information.
We also compared the forest area changes at country level. We calculated the forest area for the years 2001 and 2012 for each country; forest area change for each country was calculated by subtracting forest area of 2001 from that of 2012.

Statistical Indicators
To understand the spatial distribution of forest area changes for the five land cover datasets, we calculated the standard deviation and deviation among the five land cover datasets to describe the magnitude of the difference. For pixel [i, j] , the standard deviation and deviation were defined as follows: where i and j refers to the row number and column number, respectively; x i,j,k refers to the forest area change for pixel [i, j] of dataset k; x i,j refers to the average of forest area change in the five land cover datasets for pixel [i, j] .

Total Forest Area Change from 2001 to 2012
There was substantial inconsistency in total global forest area change among the five global land cover datasets. Figure 1 displays the total forest area change from 2001 to 2012 in 80 • N-60 • S for the five land cover datasets and FRA reports. The VCF was the only dataset that displayed an increase in forest area (1.7 × 10 6 km 2 ). Hansen exhibited the largest forest area decrease, up to 1.6 × 10 6 km 2 , followed by the forest losses of 0.4 × 10 6 km 2 for MCD12C1, 0.2 × 10 6 km 2 for CCI-LC, and 0.1 × 10 6 km 2 for the LUH2. The gap in the forest area change between VCF and Hansen was 3.3 × 10 6 km 2 , accounting for 8.8% of the total forest area worldwide (4.1 × 10 7 km 2 ) based on FRA report 2015 [31].
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 20 in forest area (1.7 × 10 6 km 2 ). Hansen exhibited the largest forest area decrease, up to 1.6 × 10 6 km 2 , followed by the forest losses of 0.4 × 10 6 km 2 for MCD12C1, 0.2 × 10 6 km 2 for CCI-LC, and 0.1 × 10 6 km 2 for the LUH2. The gap in the forest area change between VCF and Hansen was 3.3 × 10 6 km 2 , accounting for 8.8% of the total forest area worldwide (4.1 × 10 7 km 2 ) based on FRA report 2015 [31]. Among them, VCF indicated the largest forest gain, followed by FRA. In the sub tropics, VCF indicated the largest forest gain (1.9 × 10 5 km 2 ), and Hansen indicated the largest forest loss (1.5 × 10 5 km 2 ). As for the tropics (23.5°S-23.5°N), all the datasets indicated forest loss. FRA reported the largest forest loss area (10.0 × 10 5 km 2 ), which approximately seven times that of VCF (1.2 × 10 5 km 2 ). Hansen also indicated a large forest loss (8.5 × 10 5 km 2 ). The tropics was the only climatic zone where all the datasets indicated a same trend (decrease) in forest change. This conclusion was consistent with previous studies: Hansen et al. [45] quantified the gross forest cover loss in the tropics as 47.6 × 10 4 km 2 from 2000 to 2005, and Zeng et al. [24] found that there was a rapid forest loss in Southeast Asia in the early 21st century (about 29.3 × 10 4 km 2 ). In general, VCF displayed increases in forest area in the boreal, temperate, and subtropical zones; it also indicated the least forest gain in the tropics. Hansen reported the largest forest loss among these datasets for all climatic zones except for the tropics.
We also compared the latitudinal profiles of forest area changes based on the five datasets. In Figure 2e, VCF exhibited forest area gains in the northern hemisphere and the southern temperate zone but identified forest losses in the tropics and the southern hemisphere. VCF also had the largest latitudinal average in forest area change. Hansen showed gross forest loss along the latitudes. The latitudinal profiles of forest area changes for MCD12C1, CCI-LC, and LUH2 fluctuated with latitude, but all the three datasets indicated overall forest losses. Most of the inconsistencies occurred in the northern hemisphere, where all the datasets demonstrated decreases in the 10°S-10°N latitudinal range, which corresponds to the Amazon and Southeast Asia. The difference in the latitudinal profile of forest area change between VCF and Hansen was the largest, which was consistent with the results displayed in Figure 1. , which approximately seven times that of VCF (1.2 × 10 5 km 2 ). Hansen also indicated a large forest loss (8.5 × 10 5 km 2 ). The tropics was the only climatic zone where all the datasets indicated a same trend (decrease) in forest change. This conclusion was consistent with previous studies: Hansen et al. [45] quantified the gross forest cover loss in the tropics as 47.6 × 10 4 km 2 from 2000 to 2005, and Zeng et al. [24] found that there was a rapid forest loss in Southeast Asia in the early 21st century (about 29.3 × 10 4 km 2 ). In general, VCF displayed increases in forest area in the boreal, temperate, and subtropical zones; it also indicated the least forest gain in the tropics. Hansen reported the largest forest loss among these datasets for all climatic zones except for the tropics. We also compared the latitudinal profiles of forest area changes based on the five datasets. In Figure 2e, VCF exhibited forest area gains in the northern hemisphere and the southern temperate zone but identified forest losses in the tropics and the southern hemisphere. VCF also had the largest latitudinal average in forest area change. Hansen showed gross forest loss along the latitudes. The latitudinal profiles of forest area changes for MCD12C1, CCI-LC, and LUH2 fluctuated with latitude, but all the three datasets indicated overall forest losses. Most of the inconsistencies occurred Remote Sens. 2020, 12, 3502 7 of 18 in the northern hemisphere, where all the datasets demonstrated decreases in the 10 • S-10 • N latitudinal range, which corresponds to the Amazon and Southeast Asia. The difference in the latitudinal profile of forest area change between VCF and Hansen was the largest, which was consistent with the results displayed in Figure 1.
Forest area changes in different continents were displayed in Figure 3. VCF demonstrated the largest forest gains in Asia, Northern America, Europe, and Oceania. Hansen demonstrated the largest forest losses in Asia, Northern America, Europe, and Oceania. In Africa and South America, regions with severe tropical rainforests, most datasets displayed large forest loss. For other continents, there were inconsistencies in the sign and magnitude of the forest area changes among the datasets. Forest area changes in different continents were displayed in Figure 3. VCF demonstrated the largest forest gains in Asia, Northern America, Europe, and Oceania. Hansen demonstrated the largest forest losses in Asia, Northern America, Europe, and Oceania. In Africa and South America, regions with severe tropical rainforests, most datasets displayed large forest loss. For other continents, there were inconsistencies in the sign and magnitude of the forest area changes among the datasets.    Figure 4b, the most apparent difference was that Hansen indicated forest losses of in more pixels than that in VCF, especially in Southeast Asia, Congo, and the Amazon. Combining with Figure 4c-e, the magnitudes of forest area changes in MCD12C1, CCI-LC, and LUH2 were much lower than those of VCF and Hansen. Figure 4f displays the ensemble mean of forest area change in the five land cover datasets. The Congo Rainforest, the Amazon Rainforest, and Southeast Asia displayed severe deforestation, while Eastern China, Siberia, and Eastern USA displayed forest gain. apparent difference was that Hansen indicated forest losses of in more pixels than that in VCF, especially in Southeast Asia, Congo, and the Amazon. Combining with Figure 4c-e, the magnitudes of forest area changes in MCD12C1, CCI-LC, and LUH2 were much lower than those of VCF and Hansen. Figure 4f displays the ensemble mean of forest area change in the five land cover datasets. The Congo Rainforest, the Amazon Rainforest, and Southeast Asia displayed severe deforestation, while Eastern China, Siberia, and Eastern USA displayed forest gain.  To assess the inconsistencies of forest area change between VCF and Hansen, which displayed the largest inconsistency, the differences in forest change trends between the two datasets were displayed in Figure 6. Both VCF and Hansen demonstrated forest loss in the Amazon, South Africa, and the majority of Southeast Asia; most of these areas were located in the tropics. In contrast, both VCF and Hansen demonstrated forest gain in Europe, Eastern USA, and Western Siberia. VCF and Hansen demonstrated different change trends in forest area in other areas. Several hot spots, where heavy forest changes or large differences occurred in the base map, were also selected to compare differences among the five datasets. The hot spots include North America (37.0°N-50.0°N, 95.0°W- To assess the inconsistencies of forest area change between VCF and Hansen, which displayed the largest inconsistency, the differences in forest change trends between the two datasets were displayed in Figure 6. Both VCF and Hansen demonstrated forest loss in the Amazon, South Africa, and the majority of Southeast Asia; most of these areas were located in the tropics. In contrast, both VCF and Hansen demonstrated forest gain in Europe, Eastern USA, and Western Siberia. VCF and Hansen demonstrated different change trends in forest area in other areas. Several hot spots, where heavy forest changes or large differences occurred in the base map, were also selected to compare differences among the five datasets. The hot spots include North America (

Annual Forest Area Change in Five Datasets from 2001 to 2012
Hansen displayed a relatively stable annual forest area change, while other datasets exhibited irregular changes (Figure 7). Hansen was the only dataset that displayed net forest loss for all years during the research period, and its annual forest loss areas were greater than 5.0 × 10 4 km 2 . The subfigure displayed the magnitude of annual forest area change for these datasets, which represented the range of annual forest area change. VCF had the largest magnitude in forest area change. Forest gain exhibited by VCF was larger than 2.0 × 10 6 km 2

Annual Forest Area Change in Five Datasets from 2001 to 2012
Hansen displayed a relatively stable annual forest area change, while other datasets exhibited irregular changes (Figure 7). Hansen was the only dataset that displayed net forest loss for all years during the research period, and its annual forest loss areas were greater than 5.0 × 10 4 km 2 . The subfigure displayed the magnitude of annual forest area change for these datasets, which represented the range of annual forest area change. VCF had the largest magnitude in forest area change. Forest gain exhibited by VCF was larger than 2.0 × 10 6 km 2

Comparison of Forest Area Change at Country Level
Forest area changes for these datasets at the country level were calculated. Table 2 lists the forest area changes for countries with land areas over 10.0 × 10 7 km 2 . According to Table 2

Comparison of Forest Area Change at Country Level
Forest area changes for these datasets at the country level were calculated. Table 2 lists the forest area changes for countries with land areas over 10.0 × 10 7 km 2 . According to Table 2, the datasets exhibited substantial inconsistencies at the country level. Hansen demonstrated decreases in forest area for all the selected countries. Hansen and FRA indicated very close forest area changes in Brazil, Australia, Kazakhstan, India, and Argentina, but there were huge inconsistencies in other countries. The increased forest area indicated by VCF was the largest among the five datasets. Forest loss was indicated in Brazil and Kazakhstan for VCF. From a country perspective, Hansen demonstrated reductions in Russia and India, while the other datasets exhibited increases in these two countries. All the datasets indicated forest loss in Brazil and Argentina, while they exhibited different patterns of forest changes in other countries.

Accuracy Assessment Using Global Land Cover Validation
In this study, five global land cover datasets and FRA reports were evaluated to reveal the inconsistencies in global forest changes from the perspectives of total area, spatial distribution, latitudinal profile, and annual area change from 2001 to 2012. There were huge inconsistencies in global forest changes for the above aspects among the selected datasets. Thus, assessments of data accuracy were essential to understanding the reliabilities of the selected datasets. However, the selected land cover datasets in this study have highly different spatial resolutions (Table 1). Certain very high resolution satellite images (e.g., Google Earth images and QuickBird images) were effective tools to validate several high-resolution datasets (e.g., Hansen data, 30 m spatial resolution; CCI-LC data, 250 m spatial resolution) based on visual interpretation [24] but not suitable to validate datasets with lower resolution (e.g., LUH2 data, with a spatial resolution of 0.25 × 0.25 • ).
Here, the global land cover validation data from USGS was applied to evaluate the accuracies of the selected land cover datasets. Forest cover percentage in each sample was calculated and compared with the forest cover percentages of the land cover datasets (Figure 8). Spatially, underestimation of forest cover percentage commonly occurred in the Amazon, Western USA, and Eastern Europe; on the contrary, overestimation occurred in Eastern Asia, Western Europe, and Northern USA. The correlations between the land cover datasets and the validation data were not very strong: the correlation coefficients were 0.63, 0.77, 0.59, 0.68, and 0.37 for VCF, Hansen, MCD12C1, CCI-LC, and LUH2, respectively. Hansen had the highest correlation with the validation data, while LUH2 had the lowest correlation. Zeng et al. [24] also proved that Hansen had a pretty good accuracy in Southeast Asia across Google Earth imagery, an overall accuracy of 98.4% can be achieved. Therefore, Hansen is recommended to be utilized in forest and forest changes estimates.
We note that Hansen provided the forest gain information as a binary value from 2000 to 2012, so that years when forest gain occurred were not available. To calculate forest cover percentage of Hansen, the pixel number of forest gain in the validation pixel (0.05 × 0.05 • ) for each year was assumed to be equal in this study. The accuracy of Hansen forest cover percentage based on this hypothesis was lower than the actual accuracy.

Possible issues in VCF
Previous studies have widely reported that the global forest area has decreased over the past several decades [23,28,31]. According to FRA reports, the total global forest area in the world was 41.3 × 10 6 km 2 , 40.6 × 10 6 km 2 , 40.3 × 10 6 km 2 , 40.2 × 10 6 km 2 , and 40.0 × 10 6 km 2 in 1990, 2000, 2005, 2010, and 2015, respectively [31]. It is also indicated that the gross forest cover losses in the boreal zone and the temperate zone were 3.5 × 10 5 km 2 and 1.7 × 10 5 km 2 from 2000 to 2005 [45]. However, the results based on VCF were contradictory to the general understanding of global forest changes, which indicated an increase in global forest area [27]. VCF revealed an increase of 1.7 × 10 6 km 2 in total forest area, whereas Hansen, FRA, MCD12C1, CCI_LC, and LUH2 demonstrated decreases of 1.6 × 10 6 km 2 , 0.6 × 10 6 km 2 , 0.4 × 10 6 km 2 , 0.2 × 10 6 km 2 , and 0.1 × 10 6 km 2 , respectively (Figure 1). VCF also exhibited an order of magnitude inter annual forest area change than that in the other land cover datasets. According to these results, inconsistencies among the datasets primarily occurred between VCF and the other datasets. Another implausible behavior of VCF was that its forest area change fluctuated erratically in 2004, 2005, and 2009, indicating the occurrence of rapid forest gains or forest losses (Figure 7). For trees, defined as vegetation taller than 5 m in VCF, it is unreasonable to expect such a rapid change within a short period. Remote Sens. 2020, 12, x FOR PEER REVIEW 14 of 20 Figure 8. Accuracy assessment of the land cover datasets using USGS's global land cover validation data. Spatial pattern of differences in forest cover percentage between the land cover datasets and global land cover validation data (a1-a5) and correlation coefficient between the land cover datasets and global land cover validation data (b1-b5) for VCF (1), Hansen (2), MCD12C1 (3), CCI-LC (4), and LUH2 (5). Size of dots in (a1-a5) represent the difference of forest cover percentage between the land cover datasets and validation data. 1.6 × 10 km , 0.6 × 10 km , 0.4 × 10 km , 0.2 × 10 km , and 0.1 × 10 km , respectively (Figure 1). VCF also exhibited an order of magnitude inter annual forest area change than that in the other land cover datasets. According to these results, inconsistencies among the datasets primarily occurred between VCF and the other datasets. Another implausible behavior of VCF was that its forest area change fluctuated erratically in 2004, 2005, and 2009, indicating the occurrence of rapid forest gains or forest losses (Figure 7). For trees, defined as vegetation taller than 5 m in VCF, it is unreasonable to expect such a rapid change within a short period.
To evaluate these inconsistencies between VCF and the other land cover datasets, we compared the spatial distributions of forest area changes from VCF with the land cover from MCD12C1 in 2004 and 2005, during which VCF reported drastic fluctuations in forest area. VCF indicated that the forest area increased by 2.4 × 10 5 (Figure 10b). According to the Google Earth images, short vegetation appeared in 2005, but the tree cover did not change significantly (far less than VCF indicated) (Figure 10c-e). We also validated pixels in Southeast Asia and the US with similar results. The validation indicated that VCF cannot distinguish short vegetation (e.g., crops) from forest cover sufficiently. The seasonal growth of short vegetation, such as crops may be one reason for the inconsistencies between VCF and the other datasets.

Impacts of Tree Cover Definition and Land Cover Classification System
The differences in concepts and definitions of "forest" between different land cover datasets lead to different results in forest detection, forest ranges, and forest transitions [46]. The definitions of "forest" in the selected land cover datasets were different from each other. VCF describes the "tree canopy" as consisting of tall vegetation (≥5 m in height) [27]. Hansen defines the tree as all vegetation taller than 5 m in height [28]. The IGBP classification system in MCD12C1 describes "forest" as land dominated by trees with a percent cover >60% and height exceeding 2 m [38]. CCI-LC defines forest as a parcel or unit of land of at least 0.5 hectares in size that is covered by 10% or more trees that are 5 m or taller [39]. FRA defines forest as an area over 0.5 ha with a minimum tree cover of 10%. However, "forest" in FRA reports represents a kind of land use rather than physical trees compared to the satellite-based land cover datasets [1]. An area will be classified as forest if it is registered as "forest" land use even if there is no tree. Thus, forest area measurement in different land cover datasets were not consistent. For example, some short shrubs, which are taller than 2 m, but shorter than 5 m, will be counted as forest in MCD12C1 but will not be counted in VCF and Hansen. Besides, different land cover classification systems in land cover datasets also lead to inconsistencies in forest area, and we cannot compare forest area change directly. In addition, for spatial resolution of the land cover datasets, algorithms to identify forest can also lead to inconsistencies in forest change.
Remote Sens. 2020, 12, x FOR PEER REVIEW 16 of 20 To investigate the rapid forest gain or loss in 2004, 2005, and 2009 indicated by VCF, we utilized the Google Earth images as reference data to validate the reliability of VCF. The fine resolution Google Earth images displayed land surface conditions for a visual inspection of land surface changes. We selected a target pixel in China (110.525°E, 21.775°N) as an example (the black dot in Figure 10a,b), where VCF reported a large forest increase (+19.2%) in 2004 ( Figure 10a) and a substantial forest decrease (−16.6%) in 2005 (Figure 10b). According to the Google Earth images, short vegetation appeared in 2005, but the tree cover did not change significantly (far less than VCF indicated) (Figure 10c-e)). We also validated pixels in Southeast Asia and the US with similar results. The validation indicated that VCF cannot distinguish short vegetation (e.g., crops) from forest cover sufficiently. The seasonal growth of short vegetation, such as crops may be one reason for the inconsistencies between VCF and the other datasets.

Impacts of Tree Cover Definition and Land Cover Classification System
The differences in concepts and definitions of "forest" between different land cover datasets lead to different results in forest detection, forest ranges, and forest transitions [46]. The definitions of "forest" in the selected land cover datasets were different from each other. VCF describes the "tree canopy" as consisting of tall vegetation (≥5 m in height) [27]. Hansen defines the tree as all vegetation The inconsistent definition of "forest" was not the major factor for the inconsistencies in the overall global forest area change. VCF and Hansen have the same definition of "forest", but the inconsistencies between them were the largest. VCF displayed the largest forest gain, while Hansen displayed the largest forest loss. So the differences in definition of "forest" were not the decisive factor for the inconsistencies. Another problem for the satellite-based land cover datasets is that they do not assess what happened to the land after deforestation. Regrowth after deforestation is a slow, incremental process compared to deforestation; the latter is easy to spot when comparing satellite images. In addition, it is more difficult to identify deforestation from year to year. Identifying regrowth and low shrubs is also challenging. Therefore, there is no clear answer on how global forest area has changed in the early 21st century. We recommend interpreting these datasets with caution and performing detailed local-scale spatial validation before use. We also highly advocate more accurate forest cover data with higher resolution and a better classification system.

Conclusions
This study analyzed global forest area changes from 2001 to 2012 using five global land cover datasets and FRA reports. To evaluate the inconsistencies in global forest area changes among these datasets, global forest total area change, spatial distribution, latitudinal profile, and annual area change were compared. We found large inconsistencies on forest changes among these datasets. The majority of the datasets indicated decreases in total global forest area, but VCF demonstrated a net increase trend, which was conflicted with the other datasets and previous studies. The inconsistencies in forest changes primarily occurred in the northern hemisphere, while the southern hemisphere showed more consistent results. These datasets demonstrated forest loss in the tropics, but inconsistent forest changes in other zones. Global land cover validation data from USGS was applied to evaluate the accuracies of the land cover datasets. The results indicated that all of the land cover datasets displayed mediocre accuracies in forest cover. Hansen had the highest correlation with validation data, while LUH2 had the lowest correlation. Therefore, Hansen is recommended for forest cover and forest cover change estimations.
The differences in definition of "forest", classification system, spatial resolution, and algorithm can affect the detection of forest and its changes. Due to these differences and the large uncertainties among these datasets, we need to interpret the results based on these land cover datasets with caution. Validations and assessments of forest and its changes using ground-level field data or very high-resolution imagery (e.g., Google Earth or USGS' validation dataset) are needed. In addition, some methods including data fusion and cross-dataset calibration are useful to reconcile the inconsistencies of the differences in forest definition [47,48].