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

## **He Chen 1,2, Zhenzhong Zeng 2,\*, Jie Wu 2, Liqing Peng 3, Venkataraman Lakshmi 4, Hong Yang <sup>5</sup> and Junguo Liu <sup>2</sup>**


Received: 25 September 2020; Accepted: 23 October 2020; Published: 25 October 2020

**Abstract:** 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, five global land cover datasets and Global Forest Resources Assessments (FRA) were assessed to reveal uncertainties in the global forest changes in the early 21st century. These datasets displayed substantial divergences in total area, spatial distribution, latitudinal profile, and annual area change from 2001 to 2012. These datasets also display completely divergent conclusions on forest area changes for different countries. Among the datasets, total forest area changes range from an increase of 1.7 <sup>×</sup> <sup>10</sup><sup>6</sup> km2 to a decrease of 1.6 <sup>×</sup> 106 km2. 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.

**Keywords:** forest area change; data assessment; uncertainty evaluation; inconsistency

#### **1. Introduction**

As one of the most widely distributed land cover types, forests have a total area of 40.6 <sup>×</sup> <sup>10</sup><sup>6</sup> km2, 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 <sup>×</sup> 10<sup>6</sup> km<sup>2</sup> 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 fuel usage [12,13]. Land expansions caused by human activities, e.g., farmland and urban expansion,

are also major causes of deforestation [14]. Such changes can result in significant climate consequences by initiating considerable climate feedbacks in biochemical and biogeophysical processes [7]. For example, the Amazon suffers from the most severe deforestation. As a consequence of the biogeophysical feedback, it has become drier and warmer [15]. Forest changes also affect humans through their ecological functions, such as mitigating the heat island effect [16]. Given that forest changes can greatly affect the Earth systems and human society, it is critical to obtain accurate descriptions of the status and changes of the global forest. However, there is no consensus on how global forests have changed in the past several decades.

There are two main methods to monitor forests at large scales: i.e., remote sensing monitoring and forest censuses. Many global land cover datasets were produced to quantify the global land cover situations and its changes with respect to the advances in remote sensing theory and technology [17]. Large-scale forest monitoring by remote sensing has improved in terms of its availability, accuracy, and spatiotemporal resolution [18,19]. These datasets fostered the development of large-scale forest researches [20,21]. Forest censuses by governments are another commonly used method to assess forest status. The FAO's FRA reports were the most commonly used national scale forest datasets. It can provide a perspective on the impact of policies and other human factors on forests, e.g., how the Three-North Shelter Forest Program in China significantly improved forest coverage in northern China [22].

Based on existing datasets and studies, global forests have substantially changed, including changes in total area [23], spatial distribution [24], and tree biodiversity [25]. However, these studies came to inconsistent and even contradictory conclusions on global forest changes. Hansen et al. [26] detected global forest change with the Advanced Very High Resolution Radiometer (AVHRR) data, revealing a decrease in the tree cover percentage from 1982 to 1999. Among different regions, Latin America and Southeast Asia were the dominant regions for forest loss. This finding was contradictory to that from the FRA. According to Song et al. [27], global tree cover area increased by 2.2 <sup>×</sup> 10<sup>6</sup> km2 from 1982 to 2016; the gross loss and gain of tree cover both increased, but the rate of tree cover gain was greater than that of tree cover loss, especially after the first decades of the 21st century. In contrast, Hansen et al. [28] reported that global forest area decreased by 1.6 <sup>×</sup> 10<sup>6</sup> km<sup>2</sup> from 2001 to 2012 [28]. According to FRA reports, the global forest area decreased by 1.3 <sup>×</sup> 10<sup>6</sup> km<sup>2</sup> between 1990 and 2015, including decreases of 7.3 <sup>×</sup> 10<sup>5</sup> km<sup>2</sup> from 1990 to 2001 and 5.6 <sup>×</sup> 105 km<sup>2</sup> from 2001 to 2015 [29–31]. Current studies reached a consensus: the tropics have undergone rapid forest loss over the past several decades. Keenan et al. [23] concluded the tropical forest area decreased by 1.95 <sup>×</sup> <sup>10</sup><sup>6</sup> km<sup>2</sup> from 1990 to 2015. Ordway et al. [32] found that agricultural expansion increased forest loss in South America and Southeast Asia; Qin et al. [33] also demonstrated a forest loss in the Brazilian Amazon from 2000 to 2017. Current studies also revealed that cropland expansion was an important driver of the forest loss in the tropics [24,34]. However, when it comes to the specific area of forest change, the conclusions are diverse. So far, these studies paid more attention to analyzing the results based on their own datasets, while comparisons between different datasets were often ignored. Comparisons of how forests changed based on the commonly used land cover datasets can help us to understand the states of forest and recognize the inconsistencies among the datasets.

The primary objective of this study was to investigate inconsistencies of global forest changes among five commonly used globally land cover datasets and FRA reports from 2001 to 2012. The selected land cover datasets including Vegetation Continuous Fields (VCF), Global Forest Change (Hansen), Terra and Aqua combined Moderate Resolution Imaging Spectroradiometer Land Cover Climate Modeling Grid (MCD12C1), Land Cover project of the Climate Change Initiative (CCI-LC), and the new generation of Land-Use Harmonization (LUH2). FRA reports also participated in the comparison at country and continental scales. To evaluate the inconsistencies of forest area changes among these datasets, total forest area change, latitudinal profile of forest change, forest area change in different climatic zones, annual area change, and spatial distribution of forest change for the selected datasets

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.

#### **2. Materials and Methods**

#### *2.1. Materials*

#### 2.1.1. 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).

#### 2.1.2. 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].

#### 2.1.3. 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].

#### 2.1.4. 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.

#### 2.1.5. 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. LUH2 provides global land-use states and transitions at a 0.25 × 0.25◦ resolution from 850 to 2100. The classification of LUH2 is more detailed than that of the previous versions (Table 1) [41,42].

**Table 1.** Basic information of Vegetation Continuous Fields (VCF), Global Forest Change (Hansen), Terra and Aqua combined Moderate Resolution Imaging Spectroradiometer Land Cover Climate Modeling Grid (MCD12C1), Land Cover project of the Climate Change Initiative (CCI-LC), and Land-Use Harmonization (LUH2).


#### 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.

#### 2.1.7. 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 ata2m spatial resolution, spanning 2001–2014.

#### *2.2. 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:

$$\text{Change ratio}\_{[i,j]} = \frac{num\_{\text{gain}} - num\_{\text{loss}}}{num\_{\text{load}}},\tag{1}$$

where *i* and *j* refers to the row number and column number of the 0.05 × 0.05◦ pixel; *numgain* refers to the number of forest gain in the pixel [*i*,*j*]; *numloss* refers to the number of forest loss event range from 1 to 12 in the pixel; *numland* 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.

#### *2.3. 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:

$$\text{Standard deviation}\_{[i,j]} = \sqrt{\frac{\sum\_{k=1}^{5} \left(\mathbf{x}\_{i,j,k} - \overline{\mathbf{x}}\_{i,j}\right)^{2}}{5}},\tag{2}$$

$$\text{Deviation}\_{[i,j]} = \frac{\sum\_{k=1}^{5} \left| \mathbf{x}\_{i,j,k} - \overline{\mathbf{x}}\_{i,j} \right|}{5},\tag{3}$$

where *i* and *j* refers to the row number and column number, respectively; *xi*,*j*,*<sup>k</sup>* refers to the forest area change for pixel [*i*, *j*] of dataset *k*; *xi*,*<sup>j</sup>* refers to the average of forest area change in the five land cover datasets for pixel [*i*, *j*] .

#### **3. Results**

#### *3.1. 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 <sup>×</sup> 106 km2). Hansen exhibited the largest forest area decrease, up to 1.6 <sup>×</sup> 106 km2, followed by the forest losses of 0.4 <sup>×</sup> 10<sup>6</sup> km2 for MCD12C1, 0.2 <sup>×</sup> 10<sup>6</sup> km<sup>2</sup> for CCI-LC, and 0.1 <sup>×</sup> 106 km2 for the LUH2. The gap in the forest area change between VCF and Hansen was 3.3 <sup>×</sup> 106 km2, accounting for 8.8% of the total forest area worldwide (4.1 <sup>×</sup> 107 km2) based on FRA report 2015 [31].

**Figure 1.** Total forest area changes in 60◦S–80◦N for VCF, Hansen, MCD12C1, CCI-LC, and LUH2 from 2001 to 2012.

Forest area changes across climatic zones are displayed in Figure 2. In the boreal zone (66.5◦N–80.0◦N in this study), VCF, FRA, and CCI-LC displayed increases in forest area, while the other datasets demonstrated decreases. VCF indicated the largest forest area increase (2.2 <sup>×</sup> 105 km2) (Figure 2a). In the temperate zone (35◦N–66.5◦N and 35◦S–60◦S) (Figure 2b), Hansen indicated forest loss while other datasets indicated forest gain. Among them, VCF indicated the largest forest gain, followed by FRA. In the sub tropics, VCF indicated the largest forest gain (1.9 <sup>×</sup> 105 km2), and Hansen indicated the largest forest loss (1.5 <sup>×</sup> 105 km2). 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 <sup>×</sup> 105 km2), which approximately seven times that of VCF (1.2 <sup>×</sup> 105 km2). Hansen also indicated a large forest loss (8.5 <sup>×</sup> 105 km2). 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 <sup>×</sup> 104 km<sup>2</sup> 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 <sup>×</sup> 104 km2). 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.

**Figure 2.** Forest area change in different climatic zones from 2001 to 2012: (**a**) boreal (66.5◦N–80◦N), (**b**) temperate (35◦N–66.5◦N and 35◦S–60◦S), (**c**) subtropical (23.5◦N–35◦N and 23.5◦S–35◦S), (**d**) tropical (23.5◦S–23.5◦N), and (**e**) the latitudinal profiles of forest area change in 60◦S–80◦N.

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.

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 3.** Forest area change in different continents from 2001 to 2012: (**a**) Asia, (**b**) North America, (**c**) Europe, (**d**) Africa, (**e**) South America, and (**f**) Oceania.

#### *3.2. Spatial Distribution of Forest Change from 2001 to 2012*

Figure 4 displays the spatial distribution of forest area changes in the five land cover datasets. For VCF, forest area increased in Europe, the USA, China, India, and Siberia (Figure 4a). The spatial distribution of Hansen was significantly different from that of VCF. Based on 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.

**Figure 4.** Distributions of forest area changes from 2001 to 2012: (**a**) VCF, (**b**) Hansen, (**c**) MCD12C1, (**d**) CCI-LC, (**e**) LUH2, and (**f**) ensemble mean.

The spatial patterns of standard deviation and deviation of forest area change between the five datasets are demonstrated in Figure 5. Standard deviation and deviation have very similar spatial patterns. The Congo Rainforest, Amazon Rainforest, and Southeast Asia have high standard deviation and deviations, which means discrete distributions of forest change area in these areas. Siberia, India, and Western USA have small standard deviations and deviations in forest area change.

**Figure 5.** Distributions of standard deviation and deviation of forest area changes from 2001 to 2012: (**a**) standard deviation and (**b**) deviation.

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–85.0◦W), Europe (45.0◦N–54.0◦N, 0.0–15.0◦E), Eastern Siberia (55.0◦N–66.5◦N, 90.0◦E–150.0◦E), the Amazon (10.0◦S–10.0◦N, 73.0◦W–40.0◦W), South Africa (5.0◦S–8.0◦N, 10.0◦E–40.0◦E), and Southeast Asia (10.0◦S–28.0◦N, 92.0–140.0◦E). Apart from the Amazon, where all five datasets displayed forest losses, all other regions demonstrated inconsistent trends in forest change (Figure 6d). Forest losses were clear for regions in the southern hemisphere (e.g., the Amazon and the Congo rainforests), but the uncertainties were larger for regions in the northern hemisphere.

**Figure 6.** Difference in forest area change between VCF and Hansen, and comparisons of forest area changes among the five datasets in several hotspots. (**a**) North America (37.0◦N–50.0◦N, 95.0◦W–85.0◦W), (**b**) Europe (45.0◦N–54.0◦N, 0.0–15.0◦E), (**c**) Eastern Siberia (55.0◦N–66.5◦N, 90.0◦E–150.0◦E), (**d**) Amazon (10.0◦S–10.0◦N, 73.0◦W–40.0◦W), (**e**) Southern Africa (5.0◦S–8.0◦N, 10.0◦E–40.0◦E), and (**f**) Southeast Asia (10.0◦S–28.0◦N, 92.0◦E–140.0◦E).

#### *3.3. 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×104 km2. 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 <sup>×</sup> <sup>10</sup><sup>6</sup> km<sup>2</sup> in 2004; forest loss in 2009 was larger than 1.5 <sup>×</sup> <sup>10</sup><sup>6</sup> km2. The number of deforestation years was greater than the number of afforestation years for VCF. Annual forest area

changes in MCD12C1, CCI-LC, and LUH2 fluctuated but exhibited a decreasing trends from 2001 to 2012. LUH2 demonstrated a minor total forest area changes with dramatic annual fluctuations in 2009 and 2010.

**Figure 7.** Annual forest area changes in VCF, Hansen, MCD12C1, CCI-LC, and LUH2 from 2001 to 2012.

#### *3.4. 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 <sup>×</sup> 107 km2. 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.


**Table 2.** Comparison of forest area changes during 2001–2012 for VCF, Hansen, FRA, MCD12C1, CCI-LC, and LUH2 for nine major countries with total land areas over 10.0 <sup>×</sup> <sup>10</sup><sup>7</sup> km2.

#### **4. Discussion**

#### *4.1. 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.

#### *4.2. 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 <sup>×</sup> <sup>10</sup><sup>6</sup> km2, 40.6 <sup>×</sup> <sup>10</sup><sup>6</sup> km2, 40.3 <sup>×</sup> <sup>10</sup><sup>6</sup> km2, 40.2 <sup>×</sup> <sup>10</sup><sup>6</sup> km2, and 40.0 <sup>×</sup> <sup>10</sup><sup>6</sup> km2 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 <sup>×</sup> 105 km<sup>2</sup> and 1.7 <sup>×</sup> 105 km2 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 <sup>×</sup> 106 km<sup>2</sup> in total forest area, whereas Hansen, FRA, MCD12C1, CCI\_LC, and LUH2 demonstrated decreases of 1.6 <sup>×</sup> 106 km2, 0.6 <sup>×</sup> 106 km2, 0.4 <sup>×</sup> 10<sup>6</sup> km2, 0.2 <sup>×</sup> 106 km2, and 0.1 <sup>×</sup> 106 km2, 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.

**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.

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 <sup>×</sup> <sup>10</sup><sup>5</sup> km2 in 2004 (Figure 9a) and decreased by 1.0 <sup>×</sup> <sup>10</sup><sup>5</sup> km<sup>2</sup> in 2005 (Figure 9b), but the forest distribution did not substantially change in 2004 and 2005 (Figure 9c or Figure 9d). MCD12C1 demonstrated that the main land cover types in Siberia were grassland and shrub land. Forest area did not change significantly in this region, while VCF displayed clear forest change in the corresponding period. Minor changes in the Amazon forest were indicated by MCD2C1 data, while VCF indicated a noticeable increase in 2004 and a decrease in 2005 for the same area.

**Figure 9.** VCF forest area change and MCD12C1 land cover in 2004 and 2005, (**a**) forest area change in 2004, (**b**) forest area change in 2005, (**c**) MCD12C1 land cover in 2004, and (**d**) MCD12C1 land cover in 2005.

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.

#### *4.3. 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.

**Figure 10.** Contrast of forest change between VCF and Google Earth images. (**a**) VCF forest area change ratio in 2004, (**b**) VCF forest area change ratio in 2005, (**c**) land surface of target pixel from Google Earth in December 2004, (**d**) land surface of target pixel from Google Earth in December 2005, and (**e**) land surface of target pixel from Google Earth in December 2006.

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.

#### **5. 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].

**Author Contributions:** Conceptualization, Z.Z. and J.L.; methodology, Z.Z.; validation, H.C. and Z.Z.; formal analysis, H.C.; writing—original draft preparation, H.C.; writing—review and editing, Z.Z., J.W., L.P., V.L., H.Y., and J.L.; visualization, H.C. and Z.Z.; supervision, J.L. and Z.Z.; funding acquisition, J.L. and Z.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA20060402), the National Natural Science Foundation of China (NSFC) (Grant No. 41625001; 42071022), and the start-up fund for Dr. Zeng provided by Southern University of Science and Technology (No. 29/Y01296122).

**Acknowledgments:** We thank Bruce Pengra from USGS for providing the global land cover dataset.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
