Next Article in Journal
A Deep Neural Networks Approach for Augmenting Samples of Land Cover Classification
Previous Article in Journal
Unraveling Misunderstandings about Desertification: The Paradoxical Case of the Tabernas-Sorbas Basin in Southeast Spain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatially Explicit Reconstruction of Anthropogenic Grassland Cover Change in China from 1700 to 2000

1
Key Laboratory of Land Surface Pattern and Simulation, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
School of Public Administration, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Land 2020, 9(8), 270; https://doi.org/10.3390/land9080270
Submission received: 7 July 2020 / Revised: 7 August 2020 / Accepted: 10 August 2020 / Published: 12 August 2020
(This article belongs to the Section Land Systems and Global Change)

Abstract

:
Long-term anthropogenic land use and land cover changes (LULCCs) are regarded as an important component of past global change. The past 300 years have witnessed dramatic changes in LULCC in China, and this has resulted in the large-scale conversion of natural vegetation to agricultural landscapes. Studies of past LULCC in China have mainly focused on cropland and forest; however, estimates of grassland cover remain rare due to the scarcity of grassland-related historical documents. Based on a qualitative analysis of trends in grassland cover in China and their driving forces, we devised different reconstruction methods for grassland cover in eastern and western China and then developed a 10 km database of grassland cover in China for the past 300 years. The grassland area in western China decreased from 295.54 × 106 ha in 1700 to 269.78 × 106 ha in 2000 due to the increase in population and cropland, especially in northeastern China (Heilongjiang, Jilin, and Liaoning), Gan-Ning, and Xinjiang. In eastern China, grassland is degraded secondary vegetation characterized by shrub grassland and meadow grassland, which is scattered in the hills and mountains; its area increased from 7.30 × 106 ha in 1700 to 16.43 × 106 ha in 1950 due to the increase in the degraded land caused by deforestation.

1. Introduction

Over the past, land use and land cover change (LULCC) has substantially altered the landscape of the Earth’s terrestrial surface and continuously influenced climate change through a series of biogeophysical and biogeochemical mechanisms [1,2,3,4,5,6,7]. Long-term LULCC, therefore, has important implications for global environmental change. Substantial progress toward reconstructing historical LULCCs at different spatiotemporal scales has been made in recent decades.
There are four well-known global historical land-use scenarios: the History Database of the Global Environment (HYDE) [8]; the anthropogenic land cover change scenario from 8000 years ago to AD 1850 (KK10) [9]; the dataset of cropland and pasture for 1700–2007 (SAGE) [10,11]; and the global agricultural areas and land cover dataset for last millennium (PJ) [12]. Additionally, many regional historical LULCC reconstructions have been developed, including the anthropogenic deforestation dataset for 1000 BC–AD 1850 in Europe [13], land use for 1000 BC–AD 1500 in sub-Saharan Africa [14], cropland expansion and abandonment analysis for 1850–2016 in the continental United States [15], and land use for 1880–2010 in India [16].
Global historical land-use scenarios are dedicated to the reconstruction of long-term LULCC and provide estimates of LULCC at the global and continental scales. However, due to the use of different data sources and methods of reconstruction, the results of these scenarios vary widely [8,9,17], and significant differences have been observed between the global scenarios and regional reconstructions [18,19,20,21]. This has resulted in a poor understanding of past global LULCC and demonstrates the need for promoted historical LULCC reconstructions at the regional level [22].
In China, the reconstruction of past LULCC has mainly focused on cropland [23,24,25] and forest [26,27]. However, most current studies on the history of grassland cover changes in China are qualitative analyses; quantitative reconstructions are rare because of the scarcity of grassland-related historical documents. Statistics on the area of grassland first became available in the middle of the 20th century, and satellite-based grassland monitoring began in the 1970s. Existing long-term grassland reconstructions are limited to a few time periods [28] or sporadic regions [29,30]. Many qualitative studies suggest that historically, large areas of grassland have been cleared as a result of expanded land reclamation in northern China [31,32,33]; however, the temporal and spatial changes of grassland in China remains unclear.
As one of the most widespread forms of land cover on Earth, grassland covers an area of 34.00 × 108 ha, accounting for 24% of all land area worldwide [34]. In China, grasslands are the largest terrestrial ecosystem, and the area of all types of natural grassland accounts for approximately 41% of the national land area [35]. Grassland resources also provide significant economic value and ecological function. Grasslands provide forage for animal production and play an important role in wind-breaking and sand fixation, as well as in water, soil, and biodiversity conservation since Chinese grasslands are mainly distributed in areas characterized by fragile ecological environments and/or significant topographic fluctuations [36]. Changes in grassland cover have caused a series of climatic and ecological environmental problems [37]. Thus, reconstructing long-term grassland cover in China is valuable for understanding past global environmental changes.
The past 300 years have witnessed dramatic LULCCs in China accompanied by rapid population growth and the extensive popularization of drought-tolerant crops from abroad, including corn and potatoes. During this time period, land reclamation has accelerated in the northern borderlands as well as in the hilly and mountainous areas in the south. This has resulted in the large-scale conversion of natural vegetation to agricultural landscapes and the generation of large areas of degraded land. Through decades of research, we have a good understanding of the changes in croplands and forests in China during this period [23,26]. In contrast, we know little about the spatiotemporal changes in grasslands over the past 300 years. Based on a qualitative analysis of the history of grassland changes in China, this study reconstructed historical grassland cover change by synthesizing the potential natural vegetation and data on croplands and forests over the past 300 years. Then, maps of grassland cover for the time period of 1700–2000 were created with a 10 km resolution. The results provide reliable regional data for modeling climate change and improving the global historical LULCC scenarios.

2. Materials and Methods

2.1. Data Sources

China has a long history and abundant historical documents. These historical materials contain considerable LULCC information and play an important role in understanding past LULCC. However, information that can be used to directly estimate the amount of grassland in China is rare. Therefore, in this study, existing historical data on croplands and forests were used to reconstruct historical grassland cover. The datasets used for the reconstruction are introduced as follows: Using Chinese historical tax-cropland area records for 1661–1893 and government inventories for 1913–1996, Ge et al. [38] estimated the provincial cropland area of China over the past 300 years by solving the problem of cropland conversion caused by taxation and the deliberate or inadvertent omission of cropland area during land surveys. Li et al. [23] transformed the provincial cropland area into a spatially explicit dataset—the Chinese historical cropland dataset (CHCD)—based on the suitability of land for cultivation obtained by quantifying factors closely related to the spatial distribution of cropland (including surface slope, altitude, and climatic potential productivity). Based on historical documents (including official history, local chronicles, encyclopedias, travel notes, literator notes, official documents, and collected literati works), contemporary surveys, statistics, and previous studies, He et al. [39] first estimated the Chinese provincial forest area for 1700–2000 through the use of time section correction and multi-source data calibration. Subsequently, the authors devised a forested area allocation model—the Chinese historical forest dataset (CHFD)—and allocated the provincial forested area into grids at a resolution of 10 km based on the principle that the lands with high suitability for cultivation will be deforested first, with deforestation spreading to marginal lands with lower suitability for cultivation [26]. In summary, the CHCD and CHFD depend on trustworthy historical documents and appropriate reconstruction methods; and compared to historical documents and remote sensing-based data, these two datasets can capture historical changes in forest and cropland areas in China well at the provincial and grid scales. See Table 1 for more details.
Remote sensing-derived LULCC data and potential natural vegetation data are important auxiliary data for reconstructing historical grassland cover (Table 1). Using the human–computer interactive interpretation method along with Landsat TM (Thematic Mapper)/ETM (Enhanced Thematic Mapper) digital images covering China, China’s land use/cover datasets (CLUDs) were created [41,42]. The CLUDs contain six land use/cover types (accuracy greater than 94.3%)—cropland, forestland, grassland, waterbody, unused land, and built-up land—along with 25 subclasses (accuracy greater than 91.2%). These CLUDs are available on the Resource and Environment Data Cloud Platform (http://www.resdc.cn).
Potential natural vegetation refers to the vegetation that would most likely exist without anthropogenic activities. At present, the global potential natural vegetation (GPNV) map is widely used to explore historical LULCC [11,12,13,26]. The GPNV data were synthesized from the DISCover land cover dataset and the modeled natural vegetation [11]. The GPNV comprises 15 vegetation types at a 5 min resolution. Reliability evaluation showed that the GPNV spatial pattern for northeastern China is similar to the macro-scale native vegetation pattern [42]. Therefore, we believe that the GPNV map is applicable to China and used the GPNV data to replace the regions dominated by land use.
Topographic and meteorological data were also used in this study. The altitude data were obtained from the digital elevation model (DEM) of the Geospatial Data Cloud (http://www.gscloud.cn). On this basis, the slope data were calculated using ArcGIS 10.2 software. The meteorological data included the gridded data of annual precipitation and annual average temperature from 1980 to 2015 in China. Based on the daily observation data from more than 2400 meteorological stations in China, these data were generated by calculation and spatial interpolation (http://www.resdc.cn).

2.2. Regionalization of Grassland Cover Change in China

Large regional differences in the spatial patterns and changes in grassland cover in China occur due to the limitations of the natural environment (e.g., topography and climate) and the impacts of anthropogenic activity. Statistics show that 79% of grasslands in China occur in the northern temperate zone and the Qinghai-Tibetan plateau alpine zone (I and II in Figure 1) [35]. Historically, these two regions have mainly been inhabited by nomads, and the nomadic way of production has not yet caused serious damage to grassland vegetation over a long period of time [32]. However, in the course of securing the borderlands in western China, successive dynastic governments vigorously developed agriculture in the northern temperate zone, resulting in the extensive conversion of grasslands to croplands. Over the past 2000 years, this region has experienced four large-scale immigration and cultivation events, including during the Qin-Han, Sui-Tang, and the mid-Qing dynasties as well as throughout the mid-20th century [31,33]. Although the cold-dry climate of the Qinghai-Tibetan plateau has seriously hindered the development of agriculture in the past, grassland reclamation has still occurred in a few valleys, including the Yarlung Zangbo River and Huangshui River valleys [29,43]. As a result, the expansion of cropland is the main driver of decreases in the grassland area in the northern temperate zone and the Qinghai-Tibetan plateau alpine zone.
In eastern China, grassland vegetation is dominated by non-zonal secondary growth (III in Figure 1), which accounts for 21% of the total grassland area in China [35]. Compared to western China, grassland in eastern China is characterized by greater diversity and complexity in both distribution and succession [44]. Before agriculture emerged, the landscape in eastern China mainly comprised zonal forest vegetation with a few non-zonal secondary grasslands. Previous studies suggest that large-scale anthropogenic deforestation in the past generated a large area of “degraded land”, apart from parts of the deforested land turning into cropland [45]; secondary grassland then emerged in part of the degraded land in the suitable natural habitat [44,45,46]. Therefore, it is clear that the area of secondary grassland increased as anthropogenic deforestation expanded.
Considering that trends in grassland cover are consistent between the northern temperate zone and the Qinghai-Tibetan plateau alpine zone, these two regions were merged into the western region. Reconstruction methods for historical grassland areas were then devised based on these two divisions (eastern and western China). To facilitate analyses of grassland cover change at the provincial scale, we used the 25 provincial units adopted by the CHCD and CHFD (Figure 1).

2.3. Methodology

The scheme for reconstructing grassland cover in China from 1700 to 2000 is illustrated in Figure 2. In western China, considering that changes in grassland cover have primarily been influenced by land reclamation in the past, we used the proportion of grid-based cropland cover out of the potential extent of natural grassland (PENG) to capture the change in historical grassland cover. For eastern China, changes in secondary grassland cover are primarily the result of interactions between anthropogenic deforestation activities and natural habitats. This study developed a grassland estimation scenario for this region based on these two factors, given that climate variability in China over the past 300 years has been relatively minor [47]. This scenario is based on two assumptions: (1) The habitat range of current eastern secondary grassland among the potential extent of natural forest (PENF) is the historical possible distribution range of secondary grassland; and (2) the parts of non-forest and non-cropland in this habitat range for each period are regarded as secondary grassland.

2.3.1. Determination of the Potential Extent of Natural Forest and Grassland without Anthropogenic Activities

Several studies have demonstrated that most contemporary grassland regions monitored by remote sensing in China are historically grassland regions, and small parts might also be forested [48]. Similarly, most contemporary forest regions detected using remote sensing were also forested regions in the past [26]. Non-forest and non-grassland regions detected by remote sensing, especially croplands and built-up areas, might have also been forests or grasslands in the past [42,49]. Based on this information, we developed a method using CLUDs and GPNV data to determine PENG and PENF in China.
The method used in this study to determine PENG and PENF in the absence of land reclamation involved the following steps: (1) We first obtained data for the grid of contemporary forest, grassland cropland, built-up areas, and potential grassland and forest cover. (2) After potential forest was subtracted from the grid data of contemporary grassland, the remainder was classified as both the past and present areas of grassland distribution. The spatial pattern of the contemporary forest was regarded as the distribution range of forest in both the past and present. (3) Potential forest and grassland cover were superimposed over the current non-forest and non-grassland cover, and the overlapped grid cells were classified as historical forest cover or historical grassland cover in the non-forest and non-grassland regions. (4) Finally, the forest and grassland cover obtained in steps (2) and (3) were combined to enable the determination of PENG and PENF without land reclamation.

2.3.2. Reconstruction Method of Grassland Cover across Western China

Grassland in western China was reconstructed using the method of Yang et al. [42]. In China, both in the past and at present, croplands are mainly distributed around urban built-up areas and rural settlements; the statistics of remote sensing-derived LULCC indicate that the proportion of cropland in western China was greater than or equal to 90% at a resolution of 10 km, the corresponding grassland coverage was mostly less than 1% because the remainder was generally occupied by residential land [42]. Therefore, this study assumed historical grassland coverage to be zero in a grid if the proportion of cropland was greater than 90% in western China.
Based on this assumption, we first established PENG without land use and extracted historical cropland data across western China from CHCD. Grassland cover for 1700–2000 was further reconstructed by overlaying cropland cover over PENG. The equation used to determine the grassland area for western China in grid j (10 km resolution) in year t [GW(j, t)] takes the form:
G W ( j ,   t ) =   ×   [ i = 1 100 G P ( i ,   j ,   t )     C ( j ,   t ) ]
where GP(i, j, t) denotes the ith PENG (1 km) in grid j in year t; C(j, t) refers to the cropland area in grid j in year t; and α denotes the ratio of cropland to grassland. If C(j, t) is greater than 90 km2, the value of α is zero otherwise the value is one.

2.3.3. Reconstruction Method of Grassland Cover across Eastern China

In addition to anthropogenic deforestation, natural habitats have played important roles in changes in secondary grassland cover. Therefore, determining the common habitats suitable for the growth of secondary grassland is essential to reconstruct grassland cover in eastern China. Zhang et al. [50] explored the relationships between habitat and various grassland vegetation types by analyzing current remote sensing-derived grassland cover data along with corresponding topographic data (altitude and slope) and meteorological data (annual precipitation and annual average temperature). The results suggested that grassland vegetation in eastern China mainly comprises shrub grassland and meadow grassland, and their habitats are similar. These two types of grassland are located in the same topography; shrub grassland is mainly distributed in areas characterized by high temperature and precipitation, whereas meadow grassland is primarily distributed in areas with moderate temperature and precipitation. The details of these land cover types are given in Table 2.
The possible extent of secondary grassland vegetation in the PENF was determined based on the above analysis of the habitat of this land cover type. The forest and cropland data for different years across eastern China were then extracted from the CHFD and CHCD, respectively. Finally, we combined the possible extent of secondary grassland and the extracted historical forest and cropland data to reconstruct grassland cover for 1700–2000 in eastern China as follows:
G E ( j ,   t ) = [ i = 1 100 SG P ( i ,   j ,   t )     C ( j ,   t )     F ( j ,   t ) ]
where GE(j, t) denotes the grassland area for eastern China in grid j (10 km) in year t; SGP(i, j, t) refers to the ith possible extent of secondary grassland (1 km) in grid j in year t; and F(j, t) denotes the forested area in grid j in year t.

3. Results

3.1. Potential Extents of Natural Forest and Grassland

PENG and PENF are shown in Figure 3. In the absence of anthropogenic activities, PENG covered an area of 320.85 × 106 ha and was mainly distributed in western and northern China (Figure 3a). Eastern China was dominated by forest cover without land reclamation, and the area of PENF accounted for 397.12 × 106 ha (Figure 3b).

3.2. Total and Regional Grassland Area from 1700 to 2000

Figure 4 shows the total grassland area for China’s mainland from 1700 to 2000. The area of grassland continuously decreased over time because of land reclamation activities. The reconstruction results suggest two phases of change in the grassland area: a slow decrease from 1700 to 1950 followed by a rapid decrease from 1950 to 2000. The grassland area decreased from 302.84 × 106 ha in 1700 to 296.74 × 106 ha in 1950, representing a total loss of 6.10 × 106 ha at an annual loss rate (ALR) of 2.44 × 104 ha. In 2000, the grassland area was 278.80 × 106 ha; this translates to a loss of 17.94 × 106 ha between 1950 and 2000 at an ALR of 35.87 × 104 ha.
At the regional scale, grassland areas in western and eastern China have shown contrasting trends over the last 300 years. Grassland cover in western China exhibited a continuous decrease in this time period with three distinct phases: a rapid decrease between 1700 and 1750, a slight decrease between 1750 and 1850, and a sharp decrease between 1850 and 2000 (Figure 4), with corresponding losses in the grassland at 3.07 × 106, 3.13 × 106, and 19.56 × 106 ha, respectively, and ALRs of 6.13 × 104, 3.13 × 104, and 13.04 × 104 ha, respectively. The grassland area in eastern China exhibited a two-stage trend (Figure 4), a gradual increase in the prophase and a rapid decrease in the anaphase. The area occupied by grasslands in this region increased from 7.30 × 106 ha in 1700 to 16.43 × 106 ha in 1950, representing an increase of 9.12 × 106 ha over this 250-year period. In 2000, the grassland area was 9.03 × 106 ha; thus, 7.40 × 106 ha of grassland was lost between 1950 and 2000.

3.3. Spatial Distribution of Grassland Cover

Figure 5 illustrates the spatial pattern of grassland cover across China from 1700 to 2000. The spatial distribution of grassland over time generally remained constant; western China was dominated by grassland cover because of environmental conditions (e.g., climate, topography, and soil) and limited population levels. However, the reconstruction results indicated large-scale decreases in grassland areas in northeastern China (Heilongjiang, Jilin, and Liaoning), Xinjiang, and the Gan-Ning region. In eastern China, grassland was mainly scattered in the hills and mountains with a gradual increase in area over time.

4. Discussion

4.1. Analysis of Grassland Cover Change Over the Past 300 Years

While it is difficult to quantitatively estimate the area of grassland in China based on historical documents, this study can use these materials to qualitatively assess our reconstruction. The grassland area across western China decreased slowly before the mid-19th century, and the clearing of grassland mainly occurred in Gan-Ning, Inner Mongolia, and northeastern China (Figure 6a). Specifically, the emperors Kangxi and Yongzheng (1661–1735) exempted all taxes except the basic ones to encourage farmers in the Gan-Ning region to open up grasslands [31,51]. This resulted in the large-scale conversion of grasslands into croplands in the central and eastern parts of the Gan-Ning region as enthusiasm for land reclamation was greatly enhanced. For instance, in some local areas, cropland area increased by three times compared to the original cropland area [52]. Throughout this period, refugees entered Inner Mongolia and northeastern China, despite the government’s implementation of the “Prohibit reclamation in northeast China” policy [42]. As a result, the population of northeastern China increased from 1.01 million people in 1776 to 4.19 million people in 1850 [53]. Agriculture developed alongside this population growth, and the increase in cropland was mainly derived from grassland (Figure 6a).
Following the mid-19th century, the resultant limited cropland resources struggled to sustain the level of population, and the northern border areas were successively invaded by Russia and Japan. Against this background of internal revolt and foreign invasion, the Qing government began to encourage immigration into borderlands for the purposes of consolidating the borderlands and relieving conflicts between the rapidly increasing population and limited land resources [54]. Therefore, a large number of immigrants moved into this border area, resulting in the large-scale conversion of grasslands into croplands (Figure 6b) [31].
Data suggest that the grassland area in western China has remained in a stage of continuous decrease since the mid-20th century (Figure 6c). Indeed, Fan et al. [55] suggested that 18.2% of current cropland in China was derived from grassland, and western China experienced three stages of large-scale grassland reclamation throughout this period: the 1950s, the mid-1960s, and the early 1970s.
Secondary grassland increased continuously from the Qing dynasty to the Republic of China (Figure 6a,b). Due to population growth, natural disasters, and chaos caused by wars, some people were driven from their homes, settled in simple tents in the mountains, and lived by land reclamation. These immigrants were thus called shed people. Historically, changes in grassland have been significantly influenced by the shed people. The population increased geometrically since the Qing dynasty. The population of China increased from 160 million in the late 17th century to 400 million by the middle of the 19th century [53]. A large number of people were thus forced into the mountains to make a living. At the same time, crops from overseas that were suitable for planting in the mountains including corn and potatoes became available, leading to intensified land reclamation in these regions [33,56]. The soil fertility in mountainous areas was low; as a result, slope land was frequently abandoned, and new cropland was established. This led to large-scale deforestation in the mountains and the emergence of further secondary grassland on parts of the deforested land. Our reconstructions show that the area of secondary grassland gradually increased in the mountainous areas, including the Yanshan, Funiushan, and Qinling mountains as well as the area surrounding the Sichuan basin, which has been deeply affected by continuous anthropogenic deforestation (Figure 6a,b).
Since the 1950s, sustained reforestation and afforestation efforts have resulted in an increase in forest cover. Large-scale afforestation was implemented by the Chinese government to meet the needs of the economic reconstruction from the 1950s to the 1980s [57]. After the 1980s, timber production was no longer the primary goal of forest management, and the government introduced a series of ecological programs aimed at protecting forest resources [58]. As a result, secondary grassland decreased gradually as forest cover grew (Figure 6c).
The results show that our reconstructed changes in grassland cover over the past 300 years are highly consistent with the information recorded in historical documents, including information on policies, population growth, immigration, and land reclamation.

4.2. Comparisons with CLUDs and the Vegetation Atlas of China

We attempted to verify our model results via comparison with remote sensing-based and field investigation-derived grassland data. Two available contemporary LULCC datasets in China contain data on grassland, CLUDs and the vegetation atlas of China (VAC). Based on Landsat TM/ETM digital images, the CLUDs were developed by using visual interpretation against field investigation data [40,41], and this dataset has strong timeliness. The VAC, which has a scale of 1:1 million, was published in 2001 [59]. This map was compiled by multiple scientific institutions based on a large number of field investigation data throughout the country over the past 50 years and supplemented by remote sensing data. Compared to CLUDs, it has the advantages of high accuracy and less time taken.
Using the CLUDs and VAC, the reliability of the reconstructed grassland data was evaluated at the gross, provincial, and grid scales. The grassland area in the VAC, CLUDs, and our reconstruction in 2000 was 280.00 × 106 ha [36], 290 × 106 ha, and 278.80 × 106 ha, respectively. Thus, our reconstructed area of total grassland was closer to the VAC estimate than the CLUDs value. This is likely because the CHCD and CHFD were the main data sources used to reconstruct historical grassland cover in this study. The cropland and forest data for 1950–2000 in the CHCD and CHFD were derived from national statistics based on field investigations; thus, the data acquisition methods on which national statistics are based are more similar to those of the VAC than those of the CLUDs.
Figure 7 presents the relationships between our reconstruction and contemporary LULCC data at the provincial scale. The linear regression between the CLUDs and VAC grassland areas (Figure 7a) has an r2 value of 0.990 and a slope of 1.034, indicating a strong correlation between them. Strong correlations were also observed between the provincial grassland areas of the CLUDs and this study (r2 = 0.998; slope = 1.014; Figure 7b) and between those of the VAC and this study (r2 = 0.988; slope = 0.961; Figure 7c). The correlations among these datasets with the provincial grassland area less than 10.00 × 106 ha are shown in Figure 7d–f. The values of r2 between the provincial grassland areas of the CLUDs and VAC, the CLUDs and this study, and the VAC and this study were 0.767, 0.869, and 0.759, respectively, while the slopes were 0.790, 0.643, and 0.545, respectively. Overall, the correlations between them were fairly good, indicating that the reconstructed provincial grassland area in this study was close to those estimated in CLUDs and VAC.
The CLUDs grassland distribution was compared to our reconstruction at the grid-scale because data in both are in raster format. Figure 8 shows the overall spatial patterns of grassland cover in these two datasets as well as the differences between them in 2000. The spatial distribution of grassland in our reconstruction generally agreed with the CLUDs distribution; both indicated that western China was dominated by grassland cover (Figure 8a,b). Positive differences (i.e., the reconstructed grassland area is greater than the CLUDs area) are distributed mainly in the western part of northeastern China, Yinshan Mountain in Inner Mongolia, northern Jing-Jin-Ji, Shanxi, Shaanxi, and southern Chuan-Yu (Figure 8c). Negative differences (i.e., the reconstructed grassland area is less than the CLUDs area) are distributed mainly in Taihang Mountain on the border of Shanxi and Hebei, southern Shaanxi, and western Guizhou.
The percentages of grid cells with various ranges of differences between the reconstructed and CLUDs grassland cover are shown in Table 3. Grid cells with differences of −10–0% accounted for 58.55% of all grid cells, while 17.76% of grid cells had differences of 0–10%. Grids with differences exceeding 50% accounted for only 1.12% of all grid cells, while those with differences lower than −50% accounted for only 1.20%. Thus, a larger absolute value of difference corresponded to a smaller percentage of grid cells (Table 3).

4.3. Prospects and Caveats of the Reconstruction Methods

Given the large differences in changes in grassland cover between western and eastern China, this study used zoning to reconstruct historical grassland cover. Based on comparison with historical documents and previous studies, the results of our reconstruction methods well capture the history of grassland changes in China. Due to the complexity of the changes in grassland cover in eastern China, we introduced the concept of habitat and utilized this concept to constrain the possible extent of secondary grassland in eastern China in the past. At present, this approach allows us to distinguish secondary grassland from the large amount of “degraded land” in eastern China.
Notably, China has a long history and a vast territory with disparate natural and human environments among regions, resulting in different land-use practices. The current zoning approach is not sufficient to accurately capture historical grassland cover; and the uncertainty inherent in data sources (including the CHCD, CHFD, and GPNV) can also impact the results. Thus, it should be cautious about applying the dataset in this study at local research. The availability of even more information on changes in grassland cover at the regional scale would contribute to the division of more subregions in the future and the development of more reconstruction methods, which would more accurately reflect historic spatiotemporal changes in Chinese grasslands.
The method for reconstructing grassland cover in western China reflects the fact that changes in grassland area in this region have primarily been forced by land reclamation. However, changes in this land cover have also been influenced by numerous natural and cultural factors, including climate change, desertification, war, and the development of built-up areas. For instance, historical anthropogenic activities coupled with climate change accelerated desertification and greatly affected the grassland area [60]. It remains a challenge to provide quantitative information on these factors with respect to past land cover and land use.
For eastern China, we proposed a method for reconstructing grassland cover based on limited knowledge. As shown in Figure 7d–f, the uncertainties in the grassland area obtained by our method in eastern China are higher than those in western China. However, these uncertainties for eastern China have little effect on the total grassland area and its spatial pattern because grassland cover is mainly distributed in western China. Overall, the history of changes in secondary grassland in eastern China requires further exploration.
Based on some reasonable assumptions, this study determined PENG and PENF by utilizing GPNV and CLUDs. Numerous scholars have also presented quantitative reconstructions of past land cover based on archaeological and paleoecological records [61,62,63,64]. Although these records document the time history of anthropogenic land cover at numerous individual sites, it remains difficult to reconstruct long-term, spatially explicit maps of grassland cover across China at 10 km resolution using archaeological and paleoecological records. Therefore, the synthesis of multiple data sources, including archaeological data, pollen data, historical documents, statistics, and remote sensing images, shows great potential for reconstructing changes in historical land cover.

5. Conclusions

Historically, eastern and western China have exhibited large differences in changes in grassland cover. Based on these differences, this study devised methods for the reconstruction of historical grassland cover in eastern and western China and reconstructed grassland cover for 1700 to 2000. Over the past 300 years, the grassland area in China continuously decreased, and the total decrease in area was 24.03 × 106 ha. Population growth, agricultural expansion, and large scale deforestation were identified as causes of the grassland cover change. Meanwhile, traditional extensive land use over long time periods in China has resulted in a large number of degraded land while transforming natural vegetation into agricultural land. Through comparison with historical document-based qualitative studies, remote sensing-based CLUDs for 2000, and field investigation-derived VAC for 2001, we evaluated the reliability of our reconstructed grassland data at different spatial scales. The results suggest that our method for reconstructing historical grassland cover is reliable, and the reconstructed grassland cover maps accurately reflect the spatiotemporal patterns of Chinese grasslands in the past.

Author Contributions

Conceptualization, F.Y., F.H., and S.L.; methodology, F.Y.; formal analysis, F.Y.; data curation, S.L.; writing—original draft preparation, F.Y.; writing—review and editing, F.Y., F.H., and S.L.; funding acquisition, F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (Grant No. 2017YFA0603304), the National Natural Science Foundation of China (Grant No. 41671149), and the Chinese Academy of Sciences Strategic Priority Research Program (Grant No. XDA19040101).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brovkin, V.; Claussen, M.; Driesschaert, E.; Fichefet, T.; Kicklighter, D.; Loutre, M.F.; Matthews, H.D.; Ramankutty, N.; Schaeffer, M.; Sokolov, A. Biogeophysical effects of historical land cover changes simulated by six Earth system models of intermediate complexity. Clim. Dyn. 2006, 26, 587–600. [Google Scholar] [CrossRef]
  2. Canadell, J.G.; Corinne, L.Q.; Raupach, M.R.; Field, C.B.; Buitenhuis, E.T.; Philippe, C.; Conway, T.J.; Gillett, N.P.; Houghton, R.A.; Gregg, M. Contributions to accelerating atmospheric CO2 growth from economic activity, carbon intensity, and efficiency of natural sinks. Proc. Natl. Acad. Sci. USA 2007, 104, 18866–18870. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Cook, B.I.; Anchukaitis, K.J.; Kaplan, J.O.; Puma, M.J.; Kelley, M.; Gueyffier, D. Pre-Columbian deforestation as an amplifier of drought in Mesoamerica. Geophys. Res. Lett. 2012, 39, L16706. [Google Scholar] [CrossRef] [Green Version]
  4. Ellis, E.C.; Kaplan, J.O.; Fuller, D.Q.; Vavrus, S.; Klein Goldewijk, K.; Verburg, P.H. Used planet: A global history. Proc. Natl. Acad. Sci. USA 2013, 110, 7978–7985. [Google Scholar] [CrossRef] [Green Version]
  5. Houghton, R.A.; House, J.I.; Pongratz, J.; Van der Werf, G.R.; Defries, R.S.; Hansen, M.C.; Quéré, C.L.; Ramankutty, N. Chapter G2 Carbon emissions from land use and land-cover change. Biogeosciences 2012, 9, 5125–5142. [Google Scholar] [CrossRef] [Green Version]
  6. Ellis, E.C.; Beusen, A.H.W.; Goldewijk, K.K. Anthropogenic Biomes: 10,000 BCE to 2015 CE. Land 2020, 9, 129. [Google Scholar] [CrossRef]
  7. Hughes, R.; Weiberg, E.; Bonnier, A.; Finné, M.; Kaplan, J. Quantifying land use in past societies from cultural practice and archaeological data. Land 2018, 7, 9. [Google Scholar] [CrossRef] [Green Version]
  8. Klein Goldewijk, K.; Beusen, A.; Doelman, J.; Stehfest, E. Anthropogenic land use estimates for the Holocene-HYDE 3.2. Earth Syst. Sci. Data 2017, 9, 927–953. [Google Scholar] [CrossRef] [Green Version]
  9. Kaplan, J.O.; Ruddiman, W.F.; Crucifix, M.C.; Oldfield, F.A.; Krumhardt, K.M.; Ellis, E.C.; Ruddiman, W.F.; Lemmen, C.; Klein Goldewijk, K. Holocene carbon emissions as a result of anthropogenic land cover change. Holocene 2011, 21, 775–791. [Google Scholar] [CrossRef] [Green Version]
  10. Ramankutty, N. Global Cropland and Pasture Data from 1700–2007. Available online: http://www.geog.mcgill.ca/nramankutty/Data-sets/Datasets.html (accessed on 15 October 2018).
  11. Ramankutty, N.; Foley, J.A. Estimating historical changes in global land cover: Croplands from 1700 to 1992. Glob. Biogeochem. Cycles 1999, 13, 997–1027. [Google Scholar] [CrossRef]
  12. Pongratz, J.; Reick, C.; Raddatz, T.; Claussen, M. A reconstruction of global agricultural areas and land cover for the last millennium. Glob. Biogeochem. Cycles 2008, 22, GB3018. [Google Scholar] [CrossRef]
  13. Kaplan, J.O.; Krumhardt, K.M.; Zimmermann, N. The prehistoric and preindustrial deforestation of Europe. Quat. Sci. Rev. 2009, 28, 3016–3034. [Google Scholar] [CrossRef]
  14. Kay, A.U.; Kaplan, J.O. Human subsistence and land use in sub-Saharan Africa, 1000BC to AD1500: A review, quantification, and classification. Anthropocene 2015, 9, 14–32. [Google Scholar] [CrossRef]
  15. Yu, Z.; Lu, C. Historical cropland expansion and abandonment in the continental U.S. during 1850 to 2016. Glob. Ecol. Biogeogr. 2017, 27, 322–333. [Google Scholar] [CrossRef]
  16. Tian, H.Q.; Banger, K.; Bo, T.; Dadhwal, V.K. History of land use in India during 1880–2010: Large-scale land transformations reconstructed from satellite data and historical archives. Glob. Planet. Chang. 2014, 121, 78–88. [Google Scholar] [CrossRef] [Green Version]
  17. Gaillard, M.J.; Sugita, S.; Mazier, F.; Trondman, A.K.; Broström, A.; Hickler, T.; Kaplan, J.O.; Kjellström, E.; Kokfelt, U.; Kuneš, P.; et al. Holocene land-cover reconstructions for studies on land cover-climate feedbacks. Clim. Past 2010, 6, 483–499. [Google Scholar] [CrossRef]
  18. He, F.; Li, S.C.; Yang, F.; Li, M.J. Evaluating the accuracy of Chinese pasture data in global historical land use datasets. Sci. China Earth Sci. 2018, 61, 1685–1696. [Google Scholar] [CrossRef]
  19. Kaplan, J.O.; Krumhardt, K.; Gaillard, M.-J.; Sugita, S.; Trondman, A.-K.; Fyfe, R.; Marquer, L.; Mazier, F.; Nielsen, A. Constraining the deforestation history of Europe: Evaluation of historical land use scenarios with pollen-based land cover reconstructions. Land 2017, 6, 91. [Google Scholar] [CrossRef] [Green Version]
  20. Li, S.C.; He, F.N.; Zhang, X.Z.; Zhou, T.Y. Evaluation of global historical land use scenarios based on regional datasets on the Qinghai–Tibet Area. Sci. Total Environ. 2019, 657, 1615–1628. [Google Scholar] [CrossRef]
  21. Yang, F.; He, F.N.; Li, M.J.; Li, S.C. Evaluating the reliability of global historical land use scenarios for forest data in China. J. Geogr. Sci. 2020, 30, 1083–1094. [Google Scholar] [CrossRef]
  22. Gaillard, M.J.; Morrison, K.D.; Madella, M.; Whitehouse, N. Past land-use and land-cover change: The challenge of quantification at the subcontinental to global scales. Pages Mag. 2018, 26, 1–44. [Google Scholar] [CrossRef] [Green Version]
  23. Li, S.C.; He, F.N.; Zhang, X.Z. A spatially explicit reconstruction of cropland cover in China from 1661 to 1996. Reg. Environ. Chang. 2016, 16, 417–428. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, X.; Jin, X.; Du, X.; Xiang, X.; Han, J.; Shan, W.; Fan, Y.; Zhou, Y. Multi-agent model-based historical cropland spatial pattern reconstruction for 1661–1952, Shandong Province, China. Glob. Planet. Chang. 2016, 143, 175–188. [Google Scholar] [CrossRef] [Green Version]
  25. Ye, Y.; Fang, X.Q.; Ren, Y.Y.; Zhang, X.Z.; Chen, L. Cropland cover change in Northeast China during the past 300 years. Sci. China Earth Sci. 2009, 52, 1172–1182. [Google Scholar] [CrossRef]
  26. He, F.N.; Li, S.C.; Zhang, X.Z. A spatially explicit reconstruction of forest cover in China over 1700–2000. Glob. Planet. Chang. 2015, 131, 73–81. [Google Scholar] [CrossRef]
  27. Yang, X.H.; Jin, X.B.; Xiang, X.M.; Fan, Y.T.; Shan, W.; Zhou, Y.K. Reconstructing the spatial pattern of historical forest land in China in the past 300 years. Glob. Planet. Chang. 2018, 165, 173–185. [Google Scholar] [CrossRef] [Green Version]
  28. Ge, Q.S.; Zhao, M.C.; Zheng, J.Y. Land use change of China during the 20th century. Acta Geogr. Sin. 2000, 55, 698–706. [Google Scholar]
  29. Wu, Z.L.; Liu, F.G.; Zhang, Y.L.; Chen, Q.; Zhou, Q.; Yang, D.X. Forest and grassland coverage change in the Yellow River-Huangshui River Valley in northeast Qinghai-Tibet Plateau during the Qing Dynasty. Prog. Geogr. 2016, 35, 768–778. [Google Scholar] [CrossRef]
  30. Ye, Y.; Fang, X.Q. Spatial pattern of land cover changes across Northeast China over the past 300 years. J. Hist. Geogr. 2011, 37, 408–417. [Google Scholar] [CrossRef]
  31. Miao, Y.; Lu, X.S. Primary study on the reasons of grassland reclamation of China during the historical period. Pratacult. Sci. 2008, 25, 124–129. [Google Scholar]
  32. Tao, Y.; Gao, R.P. The changes of grasslands and the animal husbandry vicissitude during the historical periods. Agric. Hist. China 1992, 62–68. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFD9093&filename=ZGNS199203008&v=MjI0NDRmYkt4RjlQTXJJOUZiSVI4ZVgxTHV4WVM3RGgxVDNxVHJXTTFGckNVUjdxZllPWm9GaURuVnIzQlB5ckY= (accessed on 25 October 2018).
  33. Zhao, G. Environmental Changes in China’s History; China Environmental Science Press: Beijing, China, 1996. [Google Scholar]
  34. FAO. FAOSTAT; Food and Agriculture Organization of the United Nations: Rome, Italy, 2015. [Google Scholar]
  35. Su, D.X. The regional distribution and productivity structure of the Chinese grassland resources. Acta Agrestia. Sin. 1994, 2, 71–77. [Google Scholar]
  36. Shen, H.H.; Zhu, Y.K.; Zhao, X.; Geng, X.Q.; Gao, S.Q.; Fang, J.Y. Analysis of current grassland resources in China. Chin. Sci. Bull. 2016, 61, 139–154. [Google Scholar] [CrossRef]
  37. Kang, L.; Han, X.G.; Zhang, Z.B.; Sun, O.J. Grassland ecosystems in China: Review of current knowledge and research advancement. Philos. Trans. R. Soc. B Biol. Sci. 2007, 362, 997–1008. [Google Scholar] [CrossRef] [PubMed]
  38. Ge, Q.S.; Dai, J.H.; He, F.N.; Zheng, J.Y.; Man, Z.M.; Zhao, Y. Spatiotemporal dynamics of reclamation and cultivation and its driving factors in parts of China during the last three centuries. Prog. Nat. Sci. 2004, 14, 605–613. [Google Scholar] [CrossRef]
  39. He, F.N.; Ge, Q.S.; Dai, J.H.; Rao, Y.J. Forest change of China in recent 300 years. J. Geogr. Sci. 2008, 18, 59–72. [Google Scholar] [CrossRef]
  40. Liu, J.Y.; Kuang, W.H.; Zhang, Z.X.; Xu, X.L.; Qin, Y.W.; Ning, J.; Zhou, W.C.; Zhang, S.W.; Dong, L.R.; Zhen, Y.C.; et al. Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s. J. Geogr. Sci. 2014, 24, 195–210. [Google Scholar] [CrossRef]
  41. Ning, J.; Liu, J.Y.; Kuang, W.H.; Xu, X.L.; Zhang, S.W.; Yan, C.Z.; Li, R.D.; Wu, S.X.; Hu, Y.F.; Du, G.M.; et al. Spatiotemporal patterns and characteristics of land-use change in China during 2010–2015. J. Geogr. Sci. 2018, 28, 547–562. [Google Scholar] [CrossRef] [Green Version]
  42. Yang, F.; He, F.N.; Li, S.C.; Li, M.J. Exploring spatiotemporal pattern of grassland cover in western China from 1661 to 1996. Int. J. Environ. Res. Public Health 2019, 16, 3160. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, Y.K.; Tao, J.P.; Liu, F.G.; Zhang, Y.L.; Chen, Q. Reconstruction of cropland spatial pattern in 1830 in the middle reaches of Yarlung Zangbo River Valley. Geogr. Res. 2015, 34, 2355–2367. [Google Scholar]
  44. Fan, J.W. The characteristics of grassland resources in southern mountain areas of China and the ways for utilization. Acta Pratacult. Sin. 1995, 4, 14–18. [Google Scholar]
  45. Ren, J.Z.; Hu, Z.Z.; Zhang, Z.H.; Hou, F.; Chen, Q.G. A preliminary discussion on grassland ecological-economic regions in China. Acta Pratacult. Sin. 1999, 8, 12–22. [Google Scholar]
  46. Zhang, X.S.; Li, B.; Shi, P.J. Development and utilization of grassland resources in southern China. J. Nat. Resour. 1998, 13, 1–7. [Google Scholar]
  47. Ge, Q.S.; Zheng, J.Y.; Hao, Z.X.; Liu, H.L. General characteristics of climate changes during the past 2000 years in China. Sci. China Earth Sci. 2013, 56, 321–329. [Google Scholar] [CrossRef]
  48. Lan, Y. Vegetation succession of the medium and low mountainous areas along the subtropics of the upper Yangtze River in recent 500 years. Geogr. Res. 2010, 29, 1182–1192. [Google Scholar] [CrossRef]
  49. Li, S.C.; He, F.N.; Zhen, Z.X. An approach to spatially explicit reconstruction of historical forest in Northeast China. J. Geogr. Sci. 2014, 24, 1022–1034. [Google Scholar] [CrossRef]
  50. Zhang, Z.X. Remote Sensing Monitoring of Land Cover in China; Star Map Press: Beijing, China, 2010. [Google Scholar]
  51. Meng, F.G. Agricultural Development in Gansu Province in the early Qing Dynasty and Its historical reflection. J. Lanzhou Acad. 2009, 74–77. Available online: http://www.cnki.com.cn/Article/CJFDTotalLZXK200901019.htm (accessed on 25 October 2018).
  52. Chen, S.Z. Lanzhou Fu Chronicle; Wencheng Press: Taipei, Taiwan, 1976. [Google Scholar]
  53. Cao, S.J. Population History of China. Vol. 5, Qing Dynasty Period; Fudan University Press: Shanghai, China, 2001. [Google Scholar]
  54. Wang, X.L. Northwest Mita in Qing Dynasty; Xinjiang People’s Press: Urumqi, China, 2012. [Google Scholar]
  55. Fan, J.W.; Zhong, H.P.; Yun, X.J. Natural grassland reclamation and its ecological impaction in recent 50 years. Grassl. China 2002, 24, 69–72. [Google Scholar]
  56. Chen, B.Q. Forest change for historical period in Jiangxi province, China. Agric. Archaeol. 1985, 2, 204–214. [Google Scholar]
  57. Song, C.H.; Zhang, Y.X. Forest cover in China from 1949 to 2006. In Reforesting Landscapes; Nagendra, H., Southworth, J., Eds.; Springer: Dordrecht, The Netherlands, 2010; pp. 341–356. [Google Scholar]
  58. Zhang, Y.Q. Deforestation and forest transition: Theory and evidence in China. In World Forests From Deforestation to Transition; Palo, M., Vanhanen, H., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2000; pp. 41–65. [Google Scholar]
  59. Hou, X.Y. 1:1,000,000 Vegetation Atlas of China; Science Press: Beijing, China, 2001. [Google Scholar]
  60. Zhang, D.; Deng, H. Historical human activities accelerated climate-driven desertification in China’s Mu Us Desert. Sci. Total Environ. 2020, 708, 134771. [Google Scholar] [CrossRef]
  61. Cao, X.; Tian, F.; Li, F.; Gaillard, M.-J.; Rudaya, N.; Xu, Q.; Herzschuh, U. Pollen-based quantitative land-cover reconstruction for northern Asia covering the last 40 ka cal BP. Clim. Past 2019, 15, 1503–1536. [Google Scholar] [CrossRef] [Green Version]
  62. Li, F.; Gaillard, M.-J.; Cao, X.; Herzschuh, U.; Sugita, S.; Tarasov, P.E.; Wagner, M.; Xu, Q.; Ni, J.; Wang, W.; et al. Towards quantification of Holocene anthropogenic land-cover change in temperate China: A review in the light of pollen-based REVEALS reconstructions of regional plant cover. Earth Sci. Rev. 2020, 203, 103119. [Google Scholar] [CrossRef]
  63. Roberts, N.; Fyfe, R.M.; Woodbridge, J.; Gaillard, M.J.; Davis, B.A.S.; Kaplan, J.O.; Marquer, L.; Mazier, F.; Nielsen, A.B.; Sugita, S.; et al. Europe’s lost forests: A pollen-based synthesis for the last 11,000 years. Sci. Rep. 2018, 8, 716. [Google Scholar] [CrossRef] [PubMed]
  64. Stephens, L.; Fuller, D.; Boivin, N.; Rick, T.; Gauthier, N.; Kay, A.; Marwick, B.; Armstrong, C.G.D.; Barton, C.M.; Denham, T.; et al. Archaeological assessment reveals Earth’s early transformation through land use. Science 2019, 365, 897–902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Land 09 00270 g001
Figure 2. Scheme for reconstructing grassland cover in China from 1700 to 2000.
Figure 2. Scheme for reconstructing grassland cover in China from 1700 to 2000.
Land 09 00270 g002
Figure 3. The potential extents of natural (a) grassland (PENG) and (b) forest (PENF) without land reclamation.
Figure 3. The potential extents of natural (a) grassland (PENG) and (b) forest (PENF) without land reclamation.
Land 09 00270 g003
Figure 4. Changes in the total and regional grassland areas from 1700 to 2000.
Figure 4. Changes in the total and regional grassland areas from 1700 to 2000.
Land 09 00270 g004
Figure 5. Spatial distribution of grassland cover in China from 1700 to 2000.
Figure 5. Spatial distribution of grassland cover in China from 1700 to 2000.
Land 09 00270 g005
Figure 6. Net changes in grassland cover in China from 1700 to 2000.
Figure 6. Net changes in grassland cover in China from 1700 to 2000.
Land 09 00270 g006
Figure 7. Provincial comparisons of grassland areas in 2000 derived from CLUDs (China’s land use/cover datasets), our reconstruction, and VAC (Vegetation Atlas of China). The correlations at the provincial scale: (a) the correlation between the VAC and CLUDs; (b) the correlation between our reconstruction and CLUDs; (c) the correlation between our reconstruction and VAC. The correlations among these datasets with the provincial grassland area less than 10.00 × 106 ha: (d) the correlation between the VAC and CLUDs; (e) the correlation between our reconstruction and CLUDs; (f) the correlation between our reconstruction and VAC.
Figure 7. Provincial comparisons of grassland areas in 2000 derived from CLUDs (China’s land use/cover datasets), our reconstruction, and VAC (Vegetation Atlas of China). The correlations at the provincial scale: (a) the correlation between the VAC and CLUDs; (b) the correlation between our reconstruction and CLUDs; (c) the correlation between our reconstruction and VAC. The correlations among these datasets with the provincial grassland area less than 10.00 × 106 ha: (d) the correlation between the VAC and CLUDs; (e) the correlation between our reconstruction and CLUDs; (f) the correlation between our reconstruction and VAC.
Land 09 00270 g007
Figure 8. Spatial patterns of grassland cover in 2000: our reconstruction (a), CLUDs (b), and the difference between them (c).
Figure 8. Spatial patterns of grassland cover in 2000: our reconstruction (a), CLUDs (b), and the difference between them (c).
Land 09 00270 g008
Table 1. Details of historical cropland (CHCD) and forest (CHFD) data for China, remotely sensed data (China’s land use/cover datasets, CLUDs), and global potential natural vegetation (GPNV) data.
Table 1. Details of historical cropland (CHCD) and forest (CHFD) data for China, remotely sensed data (China’s land use/cover datasets, CLUDs), and global potential natural vegetation (GPNV) data.
DataData SourcesTemporal CoverageSpatial ResolutionReference
CHCDHistorical documents1661–199610 km × 10 kmGe et al. [38]; Li et al. [23]
CHFDHistorical documents1700–200010 km × 10 kmHe et al. [39]; He et al. [26]
CLUDsSatellite images1980–20151 km × 1 kmLiu et al. [40]; Ning et al. [41]
GPNVSatellite images, BIOME3 5′ × 5′Ramankutty and Foley [11]
Table 2. Habitat characteristics of shrub grassland and meadow grassland [50].
Table 2. Habitat characteristics of shrub grassland and meadow grassland [50].
TypeTemperature (°C)Precipitation (mm)Altitude (m)Slope (°)
Shrub grassland8.4–13.2 (90%)>600 (90%)300–2200 (85%)<25 (90.22%)
Meadow grassland2–10 (90%)>550 (90%)300–2200 (90%)<25 (90.22%)
Note: 90%, 85%, and 90.22% refer to the proportions of grassland areas in this grassland vegetation type within the corresponding temperature, precipitation, altitude, and slope range in the statistical analysis of contemporary data (grassland cover, topographic, and meteorological data).
Table 3. Statistical classifications of the differences in grassland cover in 2000 between our reconstruction and CLUDs.
Table 3. Statistical classifications of the differences in grassland cover in 2000 between our reconstruction and CLUDs.
Difference (%)−100 to −50−50 to −40−40 to −30−30 to −20−20 to −10−10 to 0
Percentage (%)1.201.151.823.236.0958.55
Difference (%)0–1010–2020–3030–4040–5050–100
Percentage (%)17.763.942.571.590.981.12

Share and Cite

MDPI and ACS Style

Yang, F.; He, F.; Li, S. Spatially Explicit Reconstruction of Anthropogenic Grassland Cover Change in China from 1700 to 2000. Land 2020, 9, 270. https://doi.org/10.3390/land9080270

AMA Style

Yang F, He F, Li S. Spatially Explicit Reconstruction of Anthropogenic Grassland Cover Change in China from 1700 to 2000. Land. 2020; 9(8):270. https://doi.org/10.3390/land9080270

Chicago/Turabian Style

Yang, Fan, Fanneng He, and Shicheng Li. 2020. "Spatially Explicit Reconstruction of Anthropogenic Grassland Cover Change in China from 1700 to 2000" Land 9, no. 8: 270. https://doi.org/10.3390/land9080270

APA Style

Yang, F., He, F., & Li, S. (2020). Spatially Explicit Reconstruction of Anthropogenic Grassland Cover Change in China from 1700 to 2000. Land, 9(8), 270. https://doi.org/10.3390/land9080270

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

Article Metrics

Back to TopTop