Land Cover Based Landscape Pattern Dynamics of Anhui Province Using GlobCover and MCD 12 Q 1 Global Land Cover Products

The development and free distribution of global land cover (GLC) products have greatly assisted in the evolution and analysis of relationships between land cover and landscape pattern. In this study, GlobCover and MCD12Q1 GLC datasets of 2005 and 2009 were comparatively used to analyze the variation of land cover in Anhui Province, China at both the class and landscape scale. The land cover classification schemes of both datasets were firstly reclassified to six types of forestland, grassland, wetland, cropland, artificial area, and others, and then FRAGSTATS was used to calculate the landscape indices. The results showed that from 2005 to 2009, the area density of ‘cropland’ landscape decreased, and it increased for ‘wetland’ and ‘artificial area’. The landscape fragmentation of ‘forestland’ and ‘grassland’ were larger. Moreover, over the same period, the class edge (CE) of ‘cropland’ was diminished; while the CE of ‘wetland’ was enhanced and the aggregation became larger. Conversely, the aggregation and shape complexity of ‘artificial area’ remained the same. The clumpiness index (CLUMPY) of ‘cropland’ varied from 0.8995 to 0.9050, indicating a higher aggregation and more concentrated distribution. The heterogeneity index (HT) value of MCD12Q1 and GlobCover datasets varied, respectively, from 0.9642 to 0.9053 and from 0.8867 to 0.8751, demonstrating that the landscape heterogeneity of Anhui Province was reduced from 2005 to 2009. Driving force analysis (DFA) was just performed for ‘artificial area’, ‘cropland’, and ‘wetland’ according to the 2005–2009 statistical yearbook data, because they were apt to be affected by human activities over a relatively short period of time.


Introduction
Landscape ecology is largely founded on the notion that environmental patterns strongly influence ecological processes [1].It is a rapidly growing science of quantifying the ways in which ecosystems interact, of establishing a link between activities in one region and repercussions in another region [2].As a part and a focus of landscape ecology, landscape pattern has been paid more attention.Landscape pattern mainly refers to the shape, ratio, and spatial features of landscape elements, and the basic characteristics is landscape heterogeneity [3].It can be quantified in a variety of ways depending on the type of data collected, the manner in which it is collected, and the objectives of the investigation [4][5][6].At a large spatial scale, the development and integration of remote sensing (RS), geographic information system (GIS), and global positioning system (GPS) have greatly facilitated the evolution and quantification of landscape pattern analysis [7][8][9].Land cover is generally considered to be the important ecological indicator for investigating the Earth s resources and many studies have been performed concerning such an issue [10].Consequently, analysis of spatial variations in landscape patterns of land cover can be very useful for evaluating the regional ecological system qualitatively and quantitatively.
In recent years, many studies have focused on analyzing the characteristics and dynamic changes of landscape patterns in urban-suburban areas, agriculture, forest, dryland, wetland, etc., through various remote sensing imagery.Vogelmann et al. [11] investigated the effects of Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper Plus (ETM+) radiometric and geometric calibrations and corrections on landscape characterization.Fichera et al. [12] characterized the dynamics of land cover pattern and its changes during a fifty-year period , using the aerial photos (1954), and Landsat scenes (Multispectral Scanner (MSS) 1975, TM 1985and 1993, ETM+ 2004).Li et al. [13] investigated how landscape composition and configuration would affect urban heat island (UHI) in the Shanghai metropolitan region of China, based on the analysis of land surface temperature (LST) in relation to normalized difference vegetation index (NDVI), vegetation fraction (Fv), and percent impervious surface area (ISA).Qian et al. [14] quantified spatiotemporal pattern of urban greenspace using the most commonly used Landsat TM data with 30 m resolution and 2.5 m high spatial resolution imagery.del Castillo et al. [15] analyzed the spatiotemporal changes in forest cover in Moncayo Natural Park (Spain) from 1987 to 2010 using RS techniques, GIS and quantitative indices of landscape ecology.It is obvious that most studies have been performed on the evaluation of landscape pattern and changes therein with respect to a certain land cover type.Conversely, it is just highly crucial to evaluate the landscape pattern of primary land cover types at provincial, national, continental, and even global scales.The development of remote sensing has facilitated the identification of land cover and the assessment of landscape pattern [16][17][18][19].
Identification of land cover types is always the first requisite to monitor and evaluate the landscape pattern [20].The continuous development of global land cover (GLC) products, depending on remote sensing technology, has provided the land cover classification datasets at regional and global scales.Some typical GLC products have been widely used to investigate landscape pattern, such as the International Geosphere-Biosphere Program Data and Information Systems (IGBP-DIS) [21], University of Maryland (UMD) [22], Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5.1 Land Cover type (hereafter referred to as MCD12Q1) [23], GLC2000 [24], GlobCover [25], and GlobeLand30 [26].GLC products are usually produced using different remotely sensed imagery (sensors), methodologies, and validation techniques.It has become increasingly significant to evaluate the landscape pattern dynamics of land cover by harmonizing multi-source and multitemporal GLC products.
To measure or quantify the landscape pattern, it is necessary to select and calculate appropriate indicators or 'metrics' through readily available data and software [27,28].According to the characteristics of remotely sensed imagery, some studies have been carried out to characterize the landscape pattern and its dynamics of land use and land cover (LULC) using multi-stage remote sensing images and built-up patch density metrics.Ji et al. [29] characterized long-term trends and patterns of urban sprawl using both the remotely sensed data and landscape metrics.Gillanders et al. [30] identified potential and the limitations of landscape pattern indices for land cover spatial pattern analysis using three or more image dates.Zhang et al. [31] investigated the relation between soil erosion and landscape patterns using 12 landscape indices (patch index, patch cohesion index, modified Simpson's evenness index, and the aggregation index) in FRAGSTATS.In summarizing, the scale effect that the landscape response of LULC can be better revealed within appropriate spatial units.Selection of landscape pattern indices greatly depends on the spatial distribution of study objectives and landscape effect.In addition, FRAGSTATS has been widely used in exploring the landscape pattern of remote sensed based LULC, due to its popularity and good compatibility with ArcGIS.
In our study, Anhui Province, China is used as the study area and two typical GLC products of MCD12Q1 2005/2009 and GlobCover 2005/2009 are selected.When harmonizing the land cover types due to different land cover classification schemes (LCCSs), the changes in landscape patterns of primary land cover types are investigated qualitatively and quantitatively at both the class and landscape scale using the landscape structure analysis software FRAGSTATS 4.2.Additionally, the driving force analysis (DFA) of some susceptible land cover types is also performed to assist in finding out the primary influence factors.

Study Area
Anhui Province is located in the mid-latitude zone of east China at 114  1a).It is about 450 km wide from east to west, with a north-south length of 570 km, and a total area of 139 thousand and 600 km 2 , which lies in the hinterland of the Yangtze River Delta.The province is situated in the transition zone from alternating subtropical to temperate zone, with a mild and humid climate characterized by four distinct seasons.There are three main geomorphic features: plains (Huaihe & Yanjiang); hills (Jianghuai, southern Anhui & western hills); and mountains (Western mountainous region).It is a relatively ideal pilot area to assess the availability and possibility of investigating the land cover based landscape patterns using the GLC products.As a big agricultural province, cropland is the top land use type in the study area.The province is geographically divided by the Yangtze and Huaihe rivers into three natural areas of Wanzhong, Wanbei, and Wannan (Figure 1b).types due to different land cover classification schemes (LCCSs), the changes in landscape patterns of primary land cover types are investigated qualitatively and quantitatively at both the class and landscape scale using the landscape structure analysis software FRAGSTATS 4.2.Additionally, the driving force analysis (DFA) of some susceptible land cover types is also performed to assist in finding out the primary influence factors.

Study Area
Anhui Province is located in the mid-latitude zone of east China at 114°54′ ~ 119°37′ E longitude and 29°41′ ~ 34°38′ N latitude (Figure 1a).It is about 450 km wide from east to west, with a north-south length of 570 km, and a total area of 139 thousand and 600 km 2 , which lies in the hinterland of the Yangtze River Delta.The province is situated in the transition zone from alternating subtropical to temperate zone, with a mild and humid climate characterized by four distinct seasons.There are three main geomorphic features: plains (Huaihe & Yanjiang); hills (Jianghuai, southern Anhui & western hills); and mountains (Western mountainous region).It is a relatively ideal pilot area to assess the availability and possibility of investigating the land cover based landscape patterns using the GLC products.As a big agricultural province, cropland is the top land use type in the study area.The province is geographically divided by the Yangtze and Huaihe rivers into three natural areas of Wanzhong, Wanbei, and Wannan (Figure 1b).

Technical Route
Several steps are required to finish the dynamic monitoring of land cover based landscape patterns on a provincial scale (Figure 2).

Technical Route
Several steps are required to finish the dynamic monitoring of land cover based landscape patterns on a provincial scale (Figure 2).

Data Sources and Preprocessing
Considering the spatial-temporal availability of current GLC maps, two kinds of GLC products were selected including GlobCover and MCD12Q1 (Table 1).GlobCover  MCD12Q1 provides data characterizing five global land cover classification systems.They describe land cover properties derived from observations spanning a year's input of observation data from the Terra and Aqua satellites applied to depict land cover types.Specifically, the five land cover classification systems are respectively IGBP global vegetation classification scheme (Land Cover Type 1), UMD scheme (Land Cover Type 2), MODIS-derived leaf area index (LAI) and fraction of photosynthetically active radiation (LAI/fPAR) scheme (Land Cover Type 3), MODIS-derived net primary production (NPP) scheme (Land Cover Type 4), and plant functional type (PFT) scheme (Land Cover Type 5).Land Cover Type 1 is considered the optimum scheme to study the land cover in Anhui Province, China [32].

Data Sources and Preprocessing
Considering the spatial-temporal availability of current GLC maps, two kinds of GLC products were selected including GlobCover and MCD12Q1 (Table 1).GlobCover  MCD12Q1 provides data characterizing five global land cover classification systems.They describe land cover properties derived from observations spanning a year's input of observation data from the Terra and Aqua satellites applied to depict land cover types.Specifically, the five land cover classification systems are respectively IGBP global vegetation classification scheme (Land Cover Type 1), UMD scheme (Land Cover Type 2), MODIS-derived leaf area index (LAI) and fraction of photosynthetically active radiation (LAI/fPAR) scheme (Land Cover Type 3), MODIS-derived net primary production (NPP) scheme (Land Cover Type 4), and plant functional type (PFT) scheme (Land Cover Type 5).Land Cover Type 1 is considered the optimum scheme to study the land cover in Anhui Province, China [32].
To cover the whole study area, four scenes are required with the track numbers of h27v05, h27v06, h28v05, and h28v06.The original MCD12Q1 products are stored in hierarchical data format (HDF) with the sinusoidal (SIN) projection.It is highly necessary to perform some preprocessing for matching the GlobCover products.MRT was employed to finish the image mosaicking, format conversion, reprojection, and resampling.Here, the MODIS HDF was converted into Geotiff, while the projection was converted from SIN to World Geodetic System 1984 (WGS84)/Universal Transverse Mercator (UTM).In addition, the image subsetting and reclassification were also completed in ENVI (ENvironment for Visualizing Images).

GlobCover Products and Pre-Treatment
GlobCover 2005 (covering December 2004-June 2006), the first 300 m GLC product, was released by the European Space Agency (ESA) in 2008.GlobCover 2009 (covering January-December 2009) was released on 21 December 2010.These data are derived from the 300 m medium resolution imaging spectrometer (MERIS) sensor, on board the ENVISAT satellite mission based on a multi-dimensional iterative clustering algorithm.MERIS was built by the ESA at the Cannes Mandelieu Space Center located in both the towns of Cannes and Mandelieu in France.The GLC product includes 22 land cover types, according to the LCCS and the overall classification accuracy is 73%.The available GlobCover product has been produced using the WGS84 datum and is freely distributed in the Geotiff format.It was firstly projected to the UTM coordinate system, and then was just preprocessed by subsetting the image and reassigning digital number (DN) values to the reclassified land cover types.To match the minimum spatial resolution of 500 m of MCD12Q1, it was resampled to 500 m from the original resolution of 300 m using the nearest-neighbor resampling method.

Harmonization of LCCSs
To generate comparable GLC maps, it is highly necessary to reclassify the land cover categories due to different classification schemes.A total of six land cover types were acquired in ENVI by harmonizing the various LCCSs for both GlobCover and MCD12Q1 (Table 2).The original land cover types were reclassified into 'forestland, 'grassland', 'cropland', 'wetland', 'artificial area', and 'others'.In our study, the decision tree was constructed to achieve our goal in ENVI [33].To effectively realize the evolution analysis of regional landscape pattern, it is necessary to select the appropriate landscape indices [34][35][36][37].The dynamics of landscape pattern on land cover in Anhui Province was characterized, using the landscape structure analysis software FRAGSTATS 4.2, defining patches using a four-neighbor rule.FRAGSTATS is a computer software program designed to compute a wide variety of landscape metrics for categorical map patterns [38].Metrics are grouped to six types according to the aspect of landscape pattern measured.They are area and edge metrics, shape metrics, core area metrics, contrast metrics, aggregation metrics, and diversity metrics.Within each of these groups, metrics are further grouped into patch, class, and landscape metrics.In this way, we selected the landscape indices at both the class and landscape scale by comprehensively considering the structure and composition of landscape and regional scale effect (Table 3).1)) is fundamental measures of landscape composition; specifically, how much of the landscape is comprised of a particular patch type [39].NP (Equation ( 2)) simply measures the extent of subdivision or fragmentation of the patch type.PD (Equation (3)) is a limited, but fundamental, aspect of landscape pattern.The corresponding formulae are given below where a ij is the area of patch ij; A is the total landscape area; and n i is the number of patches in the landscape of patch type (class) i.

ED, LSI, and PAFRAC
ED (Equation ( 4)) measures the edge length over the unit area.LSI (Equation ( 5)) measures the perimeter-to-area ratio for the landscape as a whole, which is identical to the habitat diversity index [40].PAFRAC (Equation ( 6)) is appealing because it reflects shape complexity across a range of spatial scales (patch sizes) [41].The relevant equations are (10, 000) (4) where e ik is the total length of edge in landscape between patch types (classes) i and k; e i is the total length of the patch edge of the i type; a ij is the ij type patch area; and p ij is the perimeter (m) of patch ij.

Aggregation Index
CLUMPY (Equations ( 7) and ( 8)) is a class-level only metric computed such that it ranges from −1 when the patch type is maximally disaggregated to 1 when the patch type is maximally clumped [42].It is calculated from the adjacency matrix, which shows the frequency with which different pairs of patch types appear side-by-side on the map.IJI (Equation ( 9)) is based on patch adjacencies, not cell adjacencies like the contagion index [43].The corresponding equations are where g ii is the number of like adjacencies (joins) between pixels of patch type (class) i based on the double-count method.;g ik is the number of adjacencies (joins) between pixels of patch type (class) i and k based on the double-count method; P i is the proportion of the landscape occupied by patch type (class) i; e ik is the sum of the edge lengths between the i and k type patches, and m is the number of patch types (classes) present in the landscape.

Connectivity Index
COHESION (Equation ( 10)) measures the physical connectedness of the corresponding patch type as [44] where p ij * is the perimeter of patch ij in terms of number of cell surfaces; a ij * is the area of patch ij in terms of number of cells; and Z is the total number of cells in the landscape.

Metric Selection of Spatial Landscape Pattern
In this section, F (Equation ( 11)), HT (Equation ( 12)), and MPFD (Equation ( 13)) were selected to study the spatial pattern of landscape pattern in Anhui Province.Here, F is a quantitative index to determine the degree of landscape fragmentation, which characterizes the degree of disturbance of human activities and is related to biodiversity.HT characterizes the complexity of the landscape, by analyzing its composition and function.MPFD is another measure of shape complexity.Mean fractal dimension approaches one for shapes with simple perimeters and approaches two when shapes are more complex [45].The corresponding formulae are where P is the number of patches of the landscape; Q is the area average of the landscape type; P i is the i type of patch area to whole landscape area proportion; m is the type of landscape; P ij is the patch perimeter of type ij, in units of m; a ij is the patch area of type ij; and N is the number of patches of a certain type.

Comparison of CA, NP, and PD
In ArcGIS, the image format of the four land cover images from 2005 to 2009 was converted into a grid of 500 × 500 m grid cells.In addition, the CA, NP and PD were calculated by FRAGSTATS 4.2 on the scale pattern of landscape patches, as shown in Table 4.For the 'forestland', the CA, NP, and PD of MCD12Q1 are significantly increasing, while the CA of GlobCover is decreasing and the NP and PD are increasing.The two datasets also show different trends for the 'grassland'.The CA of MCD12Q1 is greatly increasing, whereas the NP and PD are markedly reducing.Conversely, the CA, NP, and PD of GlobCover have a slightly increasing trend.For the 'cropland', the trend of the two datasets is the same, with the values increasing and the NP and PD slightly reducing.For the 'wetland', the CA of both datasets demonstrate a decreasing tendency, in contrast to the increasing NP and PD.These changes are particularly apparent in the GlobCover.For the 'artificial area', the

Comparison of CA, NP, and PD
In ArcGIS, the image format of the four land cover images from 2005 to 2009 was converted into a grid of 500 × 500 m grid cells.In addition, the CA, NP and PD were calculated by FRAGSTATS 4.2 on the scale pattern of landscape patches, as shown in Table 4.For the 'forestland', the CA, NP, and PD of MCD12Q1 are significantly increasing, while the CA of GlobCover is decreasing and the NP and PD are increasing.The two datasets also show different trends for the 'grassland'.The CA of MCD12Q1 is greatly increasing, whereas the NP and PD are markedly reducing.Conversely, the CA, NP, and PD of GlobCover have a slightly increasing trend.For the 'cropland', the trend of the two datasets is the same, with the values increasing and the NP and PD slightly reducing.For the 'wetland', the CA of both datasets demonstrate a decreasing tendency, in contrast to the increasing NP and PD.
These changes are particularly apparent in the GlobCover.For the 'artificial area', the two indices of the two datasets are minimal.For the 'others', the increasing value of MCD12Q1 is trending opposite to the NP and PD.The three indices of GlobCover are increasing, and the CA, NP, and PD of MCD12Q1 for both periods are comparatively larger.In summarizing, the increase and decrease of CA for 'cropland' and 'wetland', respectively, are opposite to the NP and PD.In 2005, the order of the values based on MCD12Q1 was 'cropland' > 'forestland' > 'grassland' > 'wetland' > 'artificial area' > 'others', and in 2009, the order was 'cropland' > forestland' > 'wetland' > 'artificial area' > 'others', indicating 'cropland' and 'grassland' patch type area is expanding, whereas, 'forestland' patch type area is greatly reduced.In 2005, the order of NP and PD values based on MCD12Q1 was 'grassland' > 'forestland' > 'cropland' > 'wetland' > 'artificial area' > 'others', and in 2009, the order was 'forestland' > 'grassland' > 'cropland' > 'wetland' > 'artificial area' > 'others', indicating the number of patches is increasing in 'forestland', while, in contrast, 'grassland' patch number reduced.However, based on the GlobCover, the CA, NP, and PD of the sequence are unchanged, and the specific gravity is low, and the two data results are different.
In order to further measure the degree of landscape fragmentation, the landscape fragmentation index F i is introduced.The formula is where P i is the number of patches of the ith type, and Q is the average of the area of all landscape types.As shown in Table 5, a comprehensive analysis by combining the landscape index and landscape fragmentation, can be found that the F i value changes in the 'forestland', 'grassland', and 'others' during the four years is larger, the F i change of 'cropland' and 'wetland' is very small, and the F i value of 'artificial area' almost remained unchanged.The results show that the corresponding index change trend of the landscape type, with a smaller change of the value of F i is more consistent with the change of the NP and PD of the two datasets.The change of the F i is larger, the inconsistency of the corresponding index change trend is bigger.Although the spatial resolution of the MCD12Q1 and GlobCover data are resampled to 500 m, the scale effect, due to the significantly different LCCSs, can still affect the analysis of landscape fragmentation when the landscape pattern evolution analysis is carried out in an area with a large landscape fragmentation.Therefore, the two datasets show a significant difference in the analysis results.As shown in Table 6, through the horizontal comparison of ED and LSI of Anhui Province in 2005-2009, we obtain a similar trend with the NP and PD.For the 'forestland', the ED of MCD12Q1 is decreasing, while the LSI and PAFRAC are increasing, indicating that the landscape edge of 'forestland' is decreasing, and the dispersion and shape complexity increasing.The ED, LSI, and PAFRAC of GlobCover are reduced, indicating that the edge of the landscape, dispersion and shape complexity are all reduced.For the 'grassland', the same three values of MCD12Q1 are reduced, while the ED and LSI of GlobCover are increased, and the PAFRAC decreases, indicating that the shape complexity of 'grassland' is diminished.For the 'cropland', the trend of the two datasets are the same, and the three values are smaller.For the 'wetland', the three values of the two datasets become larger.For the 'artificial area', the two values of the two datasets are minimal.For the 'others' type, the ED and PAFRAC of MCD12Q1 become larger, and the LSI becomes smaller.The three values of GlobCover become larger, and the LSI and PAFRAC of MCD12Q1 for both periods (2005 and 2009) are larger than those of the corresponding GlobCover.At the same time, the ED, LSI, and PAFRAC of the different land cover types of the same data in Table 6 were compared vertically.As shown in Figure 4, we can get the order change of the three indices from large to small in four years.The proportion of the three exponential values for each land cover type based on GlobCover varies little, while the ED and LSI based on MCD12Q1 vary widely.
The order based on the size of the two datasets are 'cropland' > 'forestland' > 'grassland' > 'wetland' > 'artificial area' > 'others'.However, within the four years, according to MCD12Q1, 'cropland' significantly increased and 'grassland' significantly reduced, whereas comparatively smaller changes in these land cover types were apparent, when viewing the GlobCover.This indicates that the edge length of 'cropland' in MCD12Q1 is larger in the unit area, and that of 'grassland' between the heterogeneous landscape elements in the unit area is smaller.At the same time, compared with GlobCover, 'grassland' in MCD12Q1 was significantly smaller over the four years, and 'cropland' became larger, indicating that the shape of 'grassland' tends to be regular and the distribution is more concentrated, the shape becomes irregular, and the distribution becomes more dispersed.While the two datasets were not significantly changed, that is, the complexity of patch shape did not change significantly.

Comparison of Aggregation
As shown in Table 7, the CLUMPY of 'forestland' in MCD12Q1 are reduced from 0.7971, in 2005, to 0.2585 in 2009.Also, patch types tend to be randomly distributed, while the degree of aggregation decreases and the IJI increases from 0.6032, in 2005, to 0.6471 in 2009, indicating that the number of adjacent patch types has increased.The CLUMPY and IJI of 'forestland' in GlobCover increase over the four-year period, indicating that in 2009, both the degree of aggregation of the patch types and the number of adjacent patch types increased compared to 2005.MCD12Q1 and GlobCover show a different trend for CLUMPY, which is due to measuring the clustering index using the node matrix.The scale and resolution of the image granularity will affect the number of The order based on the size of the two datasets are 'cropland' > 'forestland' > 'grassland' > 'wetland' > 'artificial area' > 'others'.However, within the four years, according to MCD12Q1, 'cropland' significantly increased and 'grassland' significantly reduced, whereas comparatively smaller changes in these land cover types were apparent, when viewing the GlobCover.This indicates that the edge length of 'cropland' in MCD12Q1 is larger in the unit area, and that of 'grassland' between the heterogeneous landscape elements in the unit area is smaller.At the same time, compared with GlobCover, 'grassland' in MCD12Q1 was significantly smaller over the four years, and 'cropland' became larger, indicating that the shape of 'grassland' tends to be regular and the distribution is more concentrated, the shape becomes irregular, and the distribution becomes more dispersed.While the two datasets were not significantly changed, that is, the complexity of patch shape did not change significantly.

Comparison of Aggregation
As shown in Table 7, the CLUMPY of 'forestland' in MCD12Q1 are reduced from 0.7971, in 2005, to 0.2585 in 2009.Also, patch types tend to be randomly distributed, while the degree of aggregation decreases and the IJI increases from 0.6032, in 2005, to 0.6471 in 2009, indicating that the number of adjacent patch types has increased.The CLUMPY and IJI of 'forestland' in GlobCover increase over the four-year period, indicating that in 2009, both the degree of aggregation of the patch types and the number of adjacent patch types increased compared to 2005.MCD12Q1 and GlobCover show a different trend for CLUMPY, which is due to measuring the clustering index using the node matrix.The scale and resolution of the image granularity will affect the number of nodes and, hence, the CLUMPY may vary, according to the scale used.Likewise, the CLUMPY of the two datasets also show different trends in the region of the larger area of landscape fragmentation, such as 'grassland' and 'others', where the IJI becomes larger and adjacent.The number of patch types increased, while 'cropland', 'wetland', and 'artificial area', with a smaller landscape fragmentation, showed the same cluster index.By comparing the CLUMPY and IJI of the two datasets in Table 7, the differences in the CLUMPY of the 'forestland' and 'grassland' between the two datasets are evident, and for 'cropland', the IJI of the difference has also very significant changes (Figure 5).nodes and, hence, the CLUMPY may vary, according to the scale used.Likewise, the CLUMPY of the two datasets also show different trends in the region of the larger area of landscape fragmentation, such as 'grassland' and 'others', where the IJI becomes larger and adjacent.The number of patch types increased, while 'cropland', 'wetland', and 'artificial area', with a smaller landscape fragmentation, showed the same cluster index.By comparing the CLUMPY and IJI of the two datasets in Table 7, the differences in the CLUMPY of the 'forestland' and 'grassland' between the two datasets are evident, and for 'cropland', the IJI of the difference has also very significant changes (Figure 5).

Connectivity Comparison
This section chooses COHESION to study the connectivity of landscape patches.As shown in Table 8, the COHESION of the 'forestland' in MCD12Q1 dropped from 99.35 in 2005 to 62.12 in 2009, indicating the distribution became more and more broken.Conversely, the 'forestland' of COHESION in GlobCover slightly increased from 99.14 to 99.22.The results of two datasets are different, due to the spatial resolution and the LCCSs.The same reason also results in different results for the 'grassland' and 'others' of the two datasets.The COHESION of 'cropland', 'wetland', and 'artificial area' show negligible change, and have a higher patch cohesion index, and the natural connectivity is better and the distribution is greater.In section, the landscape pattern of Anhui Province is analyzed by using MCD12Q1 and GlobCover, and the landscape index of the above landscape-scale is analyzed.Table 9 lists

Driving Force Analysis
The current integrity of the planet is being stressed beyond its biological capacity, and it is more essential now to understand the interaction between human activities and natural landscapes than ever [46].DFA is considered to be an important tool for investigating the landscape pattern and changes [47].The changes in the natural environment, human and social activities will impact on the structure of land cover.These changes can also cause the spatial changes of landscape pattern.In comparison with natural driving factors, human activities produce effects on landscape pattern dynamics for land cover types, especially for a relatively short period of time [48][49][50].
Landscape change associated with exponential population growth poses major challenges to coupled human and natural systems [51].'Artificial area', 'wetland', and 'cropland' are the primary land covers that are apt to be affected by human activities.We just consider the affecting factors of

Conclusions
Both MCD12Q1 and GlobCover of 2005 and 2009 are comparatively used to identify the landscape pattern dynamics on land cover at a provincial scale.To form unified land cover types, it is highly necessary to harmonize the original LCCSs and generate the same resolution maps for the two datasets.It is hardly inevitable that some obvious differences can be found for some land cover types (e.g., forestland, cropland), due to various LCCSs, remotely sensed imagery and validation techniques.It will be more convincing to track the dynamics of landscape pattern by selecting certain types with the minimum difference among different GLC products.In general, we can find that there are slight differences for the reclassified land cover types between the two datasets over a relatively short time period.Nevertheless, there are still obvious differences for some land cover types.For example, there are significant differences for 'forestland' and 'cropland' in comparison with other types.They can be also used to reflect the dynamics of landscape pattern to a certain degree.By contrast, it will have more significance to compare the landscape pattern dynamics of land cover using a certain GLC product with longer time series (e.g., MCD12Q1).In addition, DFA is also an important tool to investigate the changes of landscape patterns.The impact of human activates on landscape pattern (e.g., cropland, artificial area, wetland) is more significant compared with the natural factors during a relatively short period of time.
Author Contributions: Jinling Zhao has gathered the experimental data and processed the time series global land cover products; Jie Wang wrote the introduction; Yu Jin unified the land cover classification schemes of GlobCover and MCD12Q1; Lingling Fan finished some of the results part for landscape-scale analysis; Chao Xu selected landscape indices; Dong Liang performed the driving force analysis of landscape pattern change; Linsheng Huang designed the general technical workflow and wrote the conclusion.All of the authors have read and approved the final manuscript.

Figure 1 .
Figure 1.Geographic location (a) and land cover types and subdivisions (b) of Anhui Province, China.
The MODIS Reprojection Tool (MRT) is firstly used to finish the preprocessing of MCD12Q1 including imaging mosaicking, projection transformation, data format conversion.Two datasets of MCD12Q1 2005/2009 and GlobCover 2005/2009 are jointly processed to form a unified basis by image subsetting, projection transformation, and resampling.Subsequently, landscape indices at both the class and landscape scale are selected to dynamically analyze the land cover based landscape patterns, according to the landscape features of the study area.Finally, DFA is performed to explain the landscape pattern dynamics of 'artificial area', 'cropland', and 'wetland'.

Figure 1 .
Figure 1.Geographic location (a) and land cover types and subdivisions (b) of Anhui Province, China.
The MODIS Reprojection Tool (MRT) is firstly used to finish the preprocessing of MCD12Q1 including imaging mosaicking, projection transformation, data format conversion.Two datasets of MCD12Q1 2005/2009 and GlobCover 2005/2009 are jointly processed to form a unified basis by image subsetting, projection transformation, and resampling.Subsequently, landscape indices at both the class and landscape scale are selected to dynamically analyze the land cover based landscape patterns, according to the landscape features of the study area.Finally, DFA is performed to explain the landscape pattern dynamics of 'artificial area', 'cropland', and 'wetland'.Sustainability 2018, 10, x FOR PEER REVIEW 4 cover classification schemes from Land Cover Classification System (LCCS) for GlobCover and IGBP global vegetation classification scheme for MCD12Q1 products have just two periods of 2005 and 2009, while MCD 12Q1 has successive temporal coverage (V051) during 2001-2013.Consequently, we have to select the years of 2005 and 2009.In addition, both of the datasets have relatively nearest spatial resolution in comparison with other GLC products.

Figure 2 .
Figure 2. General workflow chart in this study.
products have just two periods of 2005 and 2009, while MCD 12Q1 has successive temporal coverage (V051) during 2001-2013.Consequently, we have to select the years of 2005 and 2009.In addition, both of the datasets have relatively nearest spatial resolution in comparison with other GLC products.

3. 2 . 5 .
Analysis of Spatial Landscape Pattern the changes of F, HT, and MPFD over the period 2005-2009.The F value of MCD12Q1 was reduced from 0.0668 in 2005 to 0.0600 in 2009, while that of GlobCover increased from 0.0524 in 2005 to 0.0563 in 2009.It is shown that the fragmentation of landscape in Anhui Province is smaller, the GlobCover data becomes larger in MCD12Q1, and the fragmentation index of MCD12Q1 is larger than that of GlobCover and has a different scale effect.At the same time, the HT of the two datasets were decreased, indicating that the landscape heterogeneity of Anhui Province was reduced.The diversity also reduced from 2005 to 2009, and MCD12Q1 was more sensitive to landscape heterogeneity and decreased more.While the MPFD of the two datasets are larger, and the values in 2005-2009 are slightly lower.Landscape patch type 'forestland', 'grassland', and 'cropland' area changes significantly during 2005 to 2009, indicating that human activities on the land reform frequently.
the three types according to the 2005-2009 statistical yearbook data.The transportation length of railway, road, and inner river is just used to reflect the changes of 'artificial area'.It was 80,747 km in 2005, according to Anhui Statistical Yearbook, but increased to 157,316 km in 2009.It is obvious that human activities have greatly affected the landscape pattern.The study area is divided into five major landforms in the Huaihe River Plain, Jianghuai Hilly and Hilly Areas, Wanxi Mountain Hills, Wuliang Plain, and Wannan Hilly Mountains.The main nature reserves are mainly distributed in the plains along the Yangtze River and the hilly areas of southern Anhui.During the years 2005-2009, the number of nature reserves in the study area increased from 31 to 38, and the total area increased from 0.35 km 2 to 0.44 km 2 .According to the landscape index analysis, the NP of the 'wetland' increased by 125 and 62, respectively, while the PD increased by 0.05 and 0.02.Also, ED, LSI and PAFRAC changed from 181.98, 40.6286, and 1.5147 in 2005 to 164.39, 36.3044, and 1.5116, in 2009, respectively.According to the 2005-2009 statistical yearbook data, population density increased from 466.76 people/km 2 to 486.75 people/km 2 , and agricultural population also increased annually from 51.48 million in 2005 to 52.78 million in 2009.Over the four-year period, the cropland area increased from 40,924.51 to 41,712.22 km 2 .The land use mode changed from extensive to intensive, corresponding to the F decreasing from 0.0668 in 2005 to 0.0600 in 2009.Conversely, the CLUMPY of 'cropland' increased from 0.8995 to 0.9095.

ENVI MRT Driving force analysis for artificial area, cropland and wetland due to human activities
FRAGSTATSFigure 2. General workflow chart in this study.

Table 1 .
Summarization and comparison of MCD12Q1 and GlobCover.

Table 1 .
Summarization and comparison of MCD12Q1 and GlobCover.

Table 3 .
Description of the selected landscape metrics used in the study.

Table 8 .
Analysis of landscape connection based on MCD12Q1 and GlobCover in Anhui Province.

Table 9 .
Analysis of landscape indices based on MCD12Q1 and GlobCover data in Anhui Province of 2005-2009, at landscape-scale.