www.nature.com/scientificreports OPEN Historical grassland desertification changes in the Horqin Sandy Land, Northern China (1985–2013) Received: 19 August 2016 1,2 2 2 2 3 2 4 Jinya Li , Bin Xu , Xiuchun Yang , Zhihao Qin , Lina Zhao , Yunxiang Jin , Fen Zhao & Accepted: 26 April 2017 Jian Guo Published: xx xx xxxx Since rural reforms in the 1980s, both the state and local governments of China have devoted great efforts to combating desertification through a number of eco-environmental restoration campaigns, resulting in burgeoning contention at all levels of government and sparking public concern. Monitoring and accurately assessing the statuses and trends of grassland desertification are important for developing effective restoration strategies. The Horqin Sandy Land (HSL), a very typical desertified grassland (DG) with better hydrothermal conditions among sandy lands in north China, was recently selected (1985–2013) to assess the spatiotemporal dynamic performances of grassland desertification before and after implementing restoration projects. Landsat images (TM/ETM +/OLI), field investigations and expert review were integrated to form a classification scheme for the HSL. Then, spectral mixture analysis and the decision-tree method were used to extract bare-sand ratios and vegetation cover fraction dynamics. A favourable phenomenon of DG was seen to be reversed in an accelerated pace during 2001–2013, despite challenge from both climatic and anthropogenic factors. However, overexploitation of grassland (especially for farming) and ground water for irrigation has led to remarkable decreases in the ground water level in recent decades, which should be highly concerning regarding the formulation of restoration campaigns in the sandy land. Desertification is a type of land degradation that occurs in arid, semi-arid, and sub-humid regions due to the synthetic action of climatic variations and human activities . Governments around the world have publicly regard desertification as a serious threat (approximately 32% of the world’s population, 67% of countries, and 40% of the land area worldwide are directly or indirectly ae ff cted by desertification) and pledged to reduce the rate of land degradation since the implementation of the Plan of Action to Combat Desertification (PACD) . However, little progress has been made, and the Parties to the United Nations Convention to Combat Desertification (UNCCD) 3, 4 still have no accurate estimates of the extent of desertification . Desertification monitoring and assessment (M & A) must be conducted if the Parties to the UNCCD are to have reliable information that supports monitoring the effectiveness of their activities . Remote sensing for monitoring land degradation was selected as one of three topics to be discussed at the fourth Drylands Deserts and Desertification “DDD” conference held at Ben Gurion University in November 4, 6 2014 . New ideas and scientic fi support for combating desertic fi ation are urgently needed because the challenges involved in combating desertification are still considerable . In many studies, desertification monitoring has been conducted based on changes in vegetation, mainly detected by the NDVI (normalized differential vegetation index), which is derived from red and near-infrared wavebands . Vegetation indices (VIs) are robust, provide an empirical method for monitoring the dynamics of vegetation, and are easily obtained. However, establishing a direct relationship between VIs and grassland desertification is difficult because VIs are vulnerable to soil back- 8–11 ground die ff rences and other factors, such as precipitation fluctuations . Additionally, in China, commonly used methods include visual interpretation and field investigation, which are expensive, labour intensive, and State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, 100085, P.R. China. Key Laboratory of Agri-informatics of the Ministry of Agriculture, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, 3 4 100081, P.R. China. Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, P.R. China. Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese academy of Science, Beijing, 100086, China. Correspondence and requests for materials should be addressed to B.X. (email: email@example.com) Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 1 www.nature.com/scientificreports/ Figure 1. Status and spatial distribution of grassland desertification over four periods in HSL. The maps were generated by ArcGIS software (Version 10.0, ESRI, Redlands, CA, USA, http://www.esri.com/). 12, 13 oen limi ft ted at spatial-temporal scales . Given that pixels in remote sensing images are usually mixed, espe- cially in arid and semiarid environments, extracting information at the sub-pixel scale is important. Spectral mixture analysis (SMA), which is a subpixel classification method, has been promoted as an efficient method for 14–16 deriving coverage information from multispectral/hyperspectral images . However, the application of SMA in desertification research requires further investigation. Grasslands serve important and global roles because they cover approximately one fifth of the terrestrial land surface and balance the global carbon (C) budget; however, grasslands are especially vulnerable to changes . Given a population of more than 1.3 billion people and limited farmland area, efficient and sustainable grassland management is a survival issue for China because China contains the second largest area of grassland in the world (approximately 400 million ha) . However, 90% of China’s available natural grasslands exhibit varying degrees of degradation, with most grasslands being distributed in 498 counties in 18 provinces in northwest, north, and northeast China, all of which are poverty-stricken regions with fragile ecosystems. e HS Th L is the largest sandy land and has the best hydrothermal conditions among the sandy lands in China but faces high anthropogenic pressure under rapid economic development, especially animal husbandry and reclamation. It is very helpful to study the dynamic changes in grassland desertification under the interactions between climate change and anthropogenic factors for understanding the man-earth relationship and for devel- oping effective solutions. e Th specific objectives of this study were (1) to design a grassland desertification grading system applicable to HSL; (2) to investigate the suitability of Landsat images and SMA methods for grassland desertification monitoring; and (3) to monitor and analyse the dynamics of grassland desertification during the previous 30 years in the sandy land. Results Status and spatial distribution of grassland desertification during four periods. The status and spatial distribution of grassland desertification during four periods in the HSL are shown in Fig. 1, and the sta- tistical results are shown in Table 1. The results show that the landcover in the study area is mainly grassland, which accounts for 70% of the entire territory. Non-, slightly, moderately, and severely DG account for 55%, 9%, Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 2 www.nature.com/scientificreports/ 1985 1992 2001 2013 Land Type Area Proportion Area Proportion Area Proportion Area Proportion SlDG 10,431 7.55 14,618 10.59 14,923 10.81 10,410 7.55 MDG 5,374 3.89 6,749 4.89 6,794 4.92 3,265 2.37 SeDG 3,087 2.24 5,342 3.87 3,765 2.73 2,883 2.09 Total DG 18,892 13.68 26,709 19.34 25,482 18.45 16,559 12.00 Non-DG 90,894 65.81 72,659 52.62 71,855 52.03 70,171 50.87 Total Grassland 109,786 79.49 99,368 71.96 97,337 70.48 86,730 62.87 Farmland 20,932 15.16 31,677 22.94 34,236 24.79 44,098 31.97 Forestland 2,306 1.67 2,949 2.14 3,122 2.26 3,520 2.55 Water 2,579 1.87 2,115 1.53 891 0.65 815 0.59 Salinized Land 1,370 0.99 780 0.56 928 0.67 246 0.18 Human Settlement 453 0.33 683 0.49 1,115 0.81 1,993 1.44 Bare Land 679 0.49 522 0.38 467 0.34 548 0.40 Total Non-Grassland 28,320 20.51 38,727 28.04 40,759 29.52 51,221 37.13 Table 1. Areas and proportions of each class over four periods in the Horqin Sandy Land (Area: km , Proportion: % of Horqin Sandy Land). SlDG: Slightly desertified grassland, MDG: Moderately desertified grassland, SeDG: Severely desertified grassland, DG: Desertified grassland. 4%, and 3%, respectively, of the total grassland area. The total DG area covered approximately 15% of the HSL 2 2 and exhibited an increase of 7,817 km from 1985 to 1992 before gradually decreasing to 16,559 km in 2013. The distribution of DG in the HSL represents obvious zonality. Moderately desertified grassland (MDG) and severely desertified grassland (SeDG) are mainly distributed southwest of HSL, where Ongniud, Naiman, Hure, and Aohan are involved in the administrative division. In addition, SeDGs are mainly located in the triangle region that is formed by Wufendi-Wudan-Wuduntaohai in Ongniud and the Xar Moron and Jiaolai Rivers. Dynamic spatio-temporal changes in grassland desertification over four periods. Figure 2 was generated through raster calculation in ArcGIS software (Version 10.0, ESRI, Redlands, CA, USA, http://www. esri.com/) and change detection and decision tree in ENVI software (Version 4.7, Exelis Visual Information Solutions Inc., Boulder, Colorado, USA, http://www.exelisvis.com/) to represent the dynamic changes over four periods. The statistical results are shown in Table 2. Through further processing of DG data and other types of land data from various periods, an annual rate was obtained (Table 3). From 1985 to 1992, the DGs with ‘Reversed’ and ‘Significantly Reversed’ degrees were mainly distributed in the central part of the Bairin Left Banner and the northern part of the Bairin Right Banner. However, the desertification ‘Deteriorated’ and ‘Seriously Deteriorated’ grasslands were mainly distributed in the northeast regions of Ongniud, the northwest regions of Naiman and the west regions of Kailu County. During this period, the annual changes of all three DG levels were positive, especially for the SeDGs (10.44%). From 1992 to 2001, the ‘Reversed’ and ‘Significantly Reversed’ degrees were mainly interspersed among the Naiman Banner, the northern part of the Aohan Banner and the eastern part of the Ar Horqin Banner, while the ‘Deteriorated’ and ‘Seriously Deteriorated’ degrees were mainly distributed in the mid-eastern part of the Ongniud Banner and the northeast part of the Hure Banner. During this period, slightly desertified grassland (SlDG) and MDG continued to expand, but the range was reduced while SeDG and Total-DG began to shrink, suggesting an overall reversion in grassland desertification. From 2001 to 2013, the grassland desertification in HSL was mainly ‘Reversed’ and ‘Significantly Reversed’ (with annual reduction rates of 2.52%, 4.33% and 1.95% in SlDG, MDG and SeDG, respectively), except for some ‘Deteriorated’ and ‘Seriously Deteriorated’ DGs in the Horqin Left Back Banner, the central part of the Hure Banner and the Bairin Right Banner. During all three periods, the rapidly expanding farmland and human settlement areas and continuously decreasing open-water surface area demonstrated the growing pressure of land use and the potential risk of another grassland desertification. Discussion e HS Th L is located in the agriculture and pasture interlaced zone of northern China. The harsh natural environ- ment and backward economic conditions make the eco-environment here more sensitive and vulnerable. The landscape, vegetation, soil, and particularly the land use in this zone express some transitivity and volatility, which are sensitive to human activities. Mismanagement and climate variations may lead to unsustainable development 19, 20 and soil degradation, especially in arid and semi-arid regions such as Horqin . In arid and semi-arid regions, precipitation is always the most limiting factor. The average temperature increased before 2004 and then decreased (Fig. 3b). The annual precipitation varied following an M shape trend, and low amounts of rainfall occurred aer 1998. Th ft e upward trend in the annual temperature and the opposite trend in annual precipitation made the environment unsuitable for vegetation growth and the recovery of deser- tification. Furthermore, the uneven distribution of precipitation throughout the year also resulted in limitations. Overall, 70.87% of the yearly precipitation occurred during the three summer months, while only 3.24% and 2.64% of the total yearly precipitation occurred during the spring and winter months (Fig. 3a). The temperature changes followed the same distribution as precipitation. The average annual wind speed generally decreased, exhibiting an opposite distribution to precipitation. During the grass-withering period, which includes low Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 3 www.nature.com/scientificreports/ Figure 2. Dynamic spatio-temporal changes in grassland desertification over four periods. The maps were generated by ArcGIS software (Version 10.0, ESRI, Redlands, CA, USA, http://www.esri.com/). 1985–1992 1992–2001 2001–2013 1985–2013 Change Direction Area Proportion Area Proportion Area Proportion Area Proportion Deteriorated 11,291 8.17 8,890 6.43 5,456 3.95 6,963 5.04 Seriously Deteriorated 4,019 2.91 2,308 1.67 1,142 0.83 2,488 1.80 Total Deteriorated 15,311 11.08 11,198 8.11 6,598 4.78 9,451 6.84 Reversed 5,510 3.99 11,400 8.25 11,256 8.15 6,612 4.79 Significantly Reversed 1,165 0.84 3,442 2.49 4,017 2.91 3,814 2.76 Total Reversed 6,675 4.83 14,842 10.74 15,273 11.06 10,425 7.55 Grassland to Non-Grassland 17,527 12.69 11,071 8.01 18,199 13.17 32,174 23.29 Non-Grassland to Grassland 7,121 5.15 9,039 6.54 7,723 5.59 9,258 6.70 Stable 91,516 66.24 92,002 66.60 90,357 65.40 76,842 55.62 Table 2. Spatio-temporal dynamic changes over different periods (area: km ; proportion: %). Total Non- Salinized Human Bare Annual rate SlDG MDG SeDG DG DG Farmland Forestland Water Land settlement Land 1985–1992 5.73 3.66 10.44 5.91 −2.87 7.33 3.98 −2.57 −6.15 7.25 −3.30 1992–2001 0.23 0.07 −3.28 −0.51 −0.12 0.90 0.65 −6.43 2.11 7.03 −1.17 2001–2013 −2.52 −4.33 −1.95 −2.92 −0.20 2.40 1.06 −0.71 −6.12 6.56 1.45 Table 3. Annual change rates for each major category in the Horqin Sandy Land over different periods (units: %). Positive values represent an increase, and negative values represent a decrease. Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 4 www.nature.com/scientificreports/ Figure 3. Trends of climatic and human-related factors in the study area. The meteorological data are averaged values from the weather stations in Ongniud, Naiman, Jarud, Barin left and Tongliao, and the socioeconomic data are the averaged values of the statistic almanac data from these banners. (a) is the monthly averaged value and its standard deviation from 1980 to 2014. (b) and (c) are the annual average values. vegetation coverage but strong winds, the bare surfaces are vulnerable to wind erosion, particularly during the spring as temperature and evapotranspiration increase. The long drought and windy season, uneven precipitation distribution and global warming cause this eco-environment in the HSL to be fragile and sensitive . Grassland desertification is a result of the interactions between climate change and anthropogenic factors. Climate change serves as a background to the occurrence and development of desertification and guides the evolution of desertification over long temporal and spatial scales. If the harsh environment provides the back- ground for desertification, human activities provide disturbance factors. With socioeconomic development and the burgeoning demands of people for a better life, land is facing more anthropogenic pressure as an energy source. However, at specific times and locations, human activities can also influence the process of desertifi- cation and even inhibit or reverse desertification processes through particular behaviours or measures based on understanding the mechanisms of desertification. Facing the continued worsening of grassland desertifica- tion, state and local governments have made great efforts to reduce grassland desertification and have created a series of national ecological projects, such as the Grain for Green Project (1999–present), the Beijing and Tianjin Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 5 www.nature.com/scientificreports/ Sandstorm Source Treatment Project (2001–2010), and the Three-North Shelterbelt Project (1978–present). In particular, the grazing ban (2003–present), which prevents livestock from trampling or eating the grass in grass- land areas, has effectively reduced grassland desertification . In the Naiman banner, which contains the main DG with the Ongniud banner in HSL, by 2012, a total of RMB 400 million yuan from these projects had been 2 2 2 invested to cover 3497 km , including 303 km for the Three-North Shelterbelt Project, 711 km for the Grain for 2 2 the Green Project, 1933 km for the Restore Pasturage to Natural Grassland Project and 550 km for the Watershed Project . Figure 3 shows the continuous population growth and fluctuating increase in the year-end number of livestock, which significantly increased following 2001 and should be a negative factor for reversing desertifica- tion. However, as shown in this paper and in similar studies , desertification has been reversed in recent decades. e eff Th ects of these eco-environmental projects are considered to play important roles in vegetation restoration 22, 24–28 and desertification reversion . However, the total area of farmland is fluctuating and generally increasing. Large areas of rich grassland have been reclaimed to obtain economic value as grazing has become restricted. In the absence of enough fertile land, some lowland areas among sandy dunes have been expropriated, where the soil is thin with poor structure. Although grazing has been restricted and degraded grasslands have low productivity, herdsmen have compen- 29–31 sated by increasing for the number of livestock (Fig. 3c). The groundwater table also decreased linearly , set- ting up the area for the next grassland degradation. To maintain sustainable eco-environmental restoration and the effectiveness of treatment projects, on no account can we ignore the interests of the residents in the ae ff cted region . e r Th elationships between man and earth form the foundation for the survival and progress of humans and the environment. To address this problem, the Reward-compensation Mechanism of Grassland Ecology Protection (2011–present) was created to secure the futures of these eco-environmental projects. RMB 77.36 billion yuan th th was invested during the “12 Five-year Plan” period (2011–2015), and investment increased during the “13 Five-year Plan” period (2016–2020) according to the Ministry of Agriculture of P.R. China and the Ministry 33, 34 of Finance of P.R. China . A certain effect has been obtained, e.g., decreasing farmland and the inhibition of livestock growth at the end of the year. However, further efforts are still needed to achieve some interactions and synergy between these projects. Conclusion Despite the challenging climate conditions, the monitoring results presented in this paper show that grassland desertification has been in a state of reversion in the last two decades, including rapid reductions in DGs (with annual reduction rates of 2.52%, 4.33%, 1.95% and 2.92% in SlDG, MDG, SeDG and Total-DG, respectively) aer ft 2001. This indicates that the eco-environmental restoration projects are effective and successful. Nonetheless, the overexploitation of grasslands should be considered, especially the falling water tables. As always, there is a trade-off between environment and development; more eo ff rt and wisdom are needed for further restoration campaigns in the sandy land. To improve or develop effective and timely management strategies, rapid and continuous monitoring and assessment of grassland desertification are necessary. In this paper, satellite images were shown to be valuable for historical and large-scale monitoring. However, despite the objective information derived by SMA and its comparative advantages over vegetation cover based methods, several limitations in the accuracy still exist that are based on low spatial and temporal resolution images. Further efforts are still needed to overcome the instan- taneous nature of satellite images and to develop more effective methods for information extraction. In addition, the appropriate indicators and grading system are indispensable for desertification monitoring. However, the for - mulation of a universal grassland desertification monitoring indicator system is a challenging and ambitious task due to the diverse regional desertification characteristics, technological levels, and cognition levels. In this paper, the bare sand ratio was used as a single index to evaluate the grassland desertification. Although it can reduce the impact of precipitation fluctuations and its comparative advantages over vegetation cover-based methods, the treatment method is simple. Whether a single factor can comprehensively reflect the desertification status and how to obtain other direct response desertification information indicators through remote sensing data require further exploration. Materials and Methods Study area. The study area is located in the western region of Northeast China and covers a total area of 4 2 13.81 × 10 km , from longitude 117°50′E to 124°05′E and latitude 41°40′N to 46°00′N and including 18 admin- istrative counties and cities (Fig. 4). The study area is located in a transition area from the Northeast Plain to the Inner Mongolia Plateau, which extends from the hilly area of the southwestern portion of the Greater Hinggan Mountains to the West Liaohe River Plain. In addition, the study area is in a typical agro-pastoral transition zone and has a temperate continental semi-arid monsoonal climate with an annual mean precipitation of 340–450 mm and an annual mean potential evaporation of 1500–2500 mm. The mean annual temperature was 5.8–6.4 °C (greater than 7 °C in the south and 3–4 °C in the north). The precipitation throughout the year was extremely une- ven, primarily occurring during a few summer months (May to September), and the peak precipitation occurred in July or August, which corresponded to the period with the greatest amount of sunlight. This period is called the “hot rainy season”, with 85% of the yearly precipitation and 45% of the yearly sunshine days. The average annual −1 wind speed is 3.4–4.5 m s , with most windy days and windstorms occurring between March and May, which is before the rainy season. The landscape is characterized by farmland alternating with natural grassland, and the soils in the area are primarily eolian sandy soils (cover the 84.3% of the total area) and chestnut soils. The native vegetation consists of sparse-tree grasslands, which have been degraded to typical steppe and temperate steppe deserts, and the grassland soils consist of degraded sandy chestnut soils and eolian sandy soils, which are coarse in texture, poor Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 6 www.nature.com/scientificreports/ Figure 4. Location of the study area. The maps were generated by ArcGIS software (Version 10.0, ESRI, Redlands, CA, USA, http://www.esri.com/). in structure, low in soil nutrient content and weak in water and nutrient preserving capacity, producing an envi- ronment that is less favourable for vegetation restoration and desertification reversion. Data acquisition and pre-processing. Landsat data are of moderate spatial resolution and oer t ff he long- est continuous global record of space-based surface observations, creating a historical archive with unmatched quality, detail, coverage and length and offering a unique opportunity to observe anthropogenic and natural 35, 36 changes at local to global scales . The successful launch of Landsat 8 extended the 40-year Landsat record by at least 5 years, which is very important for global change research . Furthermore, Landsat images are free and widely applied, which reduces costs and allows for large scale use, especially in areas with low economic value, such as grassland desertification areas . Based on the data quality (cloud coverage of less than 10% or free of clouds), the precipitation data and our goal of analysing the influences of implementing a series of ecological protection and recovery policies for grassland desertification before and after 2002, we selected Landsat Thematic Mapper (TM5) images that were acquired in 1985, 1992 and 2001, a set of Landsat Enhanced Thematic Mapper plus (ETM + 7) images that were acquired in 2001 and Landsat Operational Land Imager (OLI) images that were acquired in 2013 (13 scenes for each period with 30-m resolution). Given the large size of the study area and the frequent occurrence of clouds over the study area, obtaining all images within a given year was difficult; thus, some images were selected that were obtained around the year. es Th e images were acquired during the region’s typical growing season (July and early September, see Supplementary Table S1 for details) and were downloaded from the United States Geological Survey (USGS) web site (http://glovis.usgs.gov) for free. A geometric correction was applied to the Landsat images of 2001 by selecting ground control points (GCPs) from 1:100,000 topographic maps to maintain the rectification error to less than one pixel. Next, these corrected images were used as base images to warp the images from 1985, 1992 and 2013 using the image-to-image regis- 18, 39 tration method . The nearest neighbour assignment was applied, yielding a RMS error of less than 0.5 pixels. Radiometric calibration was performed to calibrate these images to radiance data via the ENVI software. Then, the radiance images were converted into either band-interleaved-by-pixel (BIP) or band-interleaved-by-line (BIL) interleaves to prepare them for atmospheric correction by using the FLAASH (Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercube) module, which is also commercially available as an add-on to the ENVI soft- ware package sold by ITT Exelis Visual Information Solutions . Radiometric calibration was performed to eliminate the errors caused by the sensor, and the FLAASH module was used to eliminate the errors caused by atmospheric absorption and scattering and was used to convert the radiance data of the Landsat images to land surface albedo . Unlike many other atmospheric correction algorithms in which atmospheric parameters rely on synchronous measurements with remote sensing images, FLAASH incorporates the MODTRAN (Moderate Resolution Atmospheric Transmission) radiation transfer code and computes the atmospheric parameters based only on the image spectrum feature. A unique MODTRAN model atmosphere and a selectable aerosol type are Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 7 www.nature.com/scientificreports/ Vegetation community characteristics Grassland desertification Vegetation Vegetation Bare-sand intensity classification composition coverage (%) ratio (%) Geomorphologic features Psammophytes become Relatively moderate sand, fixed SlDG the main accompanying 45–60 15–30 sand dunes. species. Psammophytes become Moderate sand, small blowout 25–45 30–50 MDG the dominant species. pits or semi-fixed sand dunes. Vegetation is very Medium and large sand dunes, sparse, only a few SeDG <25 >50 large blowout pits, semi-mobile psammophytes or mobile sand dunes. remaining. Table 4. Remote sensing grading system of grassland desertification in Horqin. combined for dier ff ent images to ensure acceptable precision . Considering the locations and acquisition times of the images used in this study, we choose the Mid-Latitude Summer (MLS) as the MODTRAN model atmosphere. Other data from the study area, including meteorological data (e.g., rainfall, temperature) obtained from the China Meteorological Data Sharing Service System (http://cdc.cma.gov.cn/) of the National Science and Technology Infrastructure, socioeconomic data (e.g., population, cattle stocks) obtained from the national and provincial Statistic Almanac, annual land-use maps (Data Sharing Infrastructure of Earth System Science, http:// www.geodata.cn/), and grassland type and soil type data, were collected and collated. To develop a grassland desertification classification system for the study area and interpretation symbols for visual interpretation, a two-week field survey of the study area was performed in August 2013, during which extensive ground data were obtained. These data primarily consisted of spatial-location information, geomorphologic features, sample plot and quadrat information, and grazing and desertification conditions around the sample plot. Grading system for Horqin sandy desertification. To monitor the distribution and dynamic changes of grassland desertification and assess policies, a dynamic monitoring system is required with the contributions of advancing technologies and knowledge . However, the formulation of a universal grassland desertification moni- toring indicator system is a challenging and ambitious task due to the diverse regional desertification characteris- 43–45 tics, technological levels, and cognition levels . Dispersed and patchy vegetation cover and sand sheets are the main landscape characteristics of areas undergoing grassland desertic fi ation and serve as good visual indicators of environmental changes and the severity of desertification . Based on relevant studies of grassland desertification 47, 48 criteria and the national standards of Classification Indexes of Degradation, Desertification and Salinization of Natural Grasslands (GB19377-2003), we used the vegetation cover and bare-sand ratio as primary criteria for determining the desertification status. The threshold for the rank of each indicator was set according to the sta- tistical results from the field sample plots and quadrat data, which contain information regarding desertification conditions (classified based on geomorphologic features; the DG and Non-desertified grassland (Non-DG) were distinguished by selecting the apical, stable, and closest regions to the primary grassland community in the Da Qing-Gou Valley Natural Protective Region as the baseline) and the corresponding vegetation coverage data and bare-sand ratios (Table 4). The bare sand ratio was the preferred option to identify desertification because the area with low vegetation coverage caused by non-desertification factors may be overestimated in desertification degree 6, 49 via vegetation coverage . For the sake of contrastive analysis for other studies, the vegetation grading values in Table 4 were retained because they are commonly used in other studies. Furthermore, we identified the following ten land cover classes in the study area (note that shrubland was treated as grassland) combined with land use mapping, grassland type and desertification characters: Farmland, Forestland, Human settlement land, Salinized land, Water, Non-DG land, SlDG, MDG, SeDG and Bare lands . Extraction method for indicators. Pixels in remote sensing images of arid and semi-arid environments, 50, 51 particularly in areas undergoing desertification, usually have multiple components . The linear spectral mixture model (LSMM), a SMA that is widely used due to its simplicity, reasonable level of effectiveness and 52–54 interpretability , was selected to calculate the proportions of each target surface feature within a pixel. In the LSMM, the spectral signature of a given pixel is assumed to be a linear combination of each end-member spec- trum and its relative abundance, which is represented by Eq. (1): ρ= () λρ F () λλ + () i ∑ j ii j=1 (1) where j = 1, 2, … m is the pixel component (end-member) and i = 1, 2, … n is the spectral band. In addition, the 55, 56 m n number of end-members was constrained by the dimensionality of the satellite images ( ≤ + 1) ; ρ() λ is the composite reflectance of the mixed pixels in band i; is the reflectance of the end-member j in band i; ρ () λ F i j refers to the abundance of the end-member j in the pixel (a parameter to be estimated); and () λ is the difference between the actual and modelled reflectance. In addition, represents the best-fit coefficient that minimizes the RMS error, which is obtained using Eq. (2): Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 8 www.nature.com/scientificreports/ Figure 5. Averaged reflectances of endmembers. Figure 6. Scatter plot correction between the field-measured and LSMM-estimated vegetation fractions in (( ελ )) i=1 i RMS = (2) where n is the number of bands and () λ is the residual term in band i. Appropriate selection of the end-members representing the surface components is crucial for the success of LSMM . LSMM with inadequate end-members easily results in the misdivision of the un-modelled land cover components into modelled components (end-members), while many end-members become sensitive to noise . The endmember should be representative and an effective component of most of the images within the pixel to produce physically realistic proportions, minimize the RMS error and force the derived fractions of the 55, 59 end-members constrained to non-negative values with a sum of 1 . In this study, we conducted a compre- hensive field investigation to identify the optimal or most representative surface types and identify preferable 54, 60 locations for the development of a spectral library from the Landsat images . Based on the field investigations, “pure” spectral end-members of vegetation, bare sand and bare soil were extracted using the minimum noise frac- tion (MNF) and the Pixel Purity Index algorithm (PPI) by iteratively testing different endmember combinations 61–64 (Fig. 5), which is a common approach for selecting end-members and is available through ENVI . To focus on desertification information and simplify the process of endmember selection, the normalized difference infrared index (NDII( =− ρρ )/() ρρ + , where ρ and ρ are the reflectances in the 819 1649 819 1649 819 1649 electromagnetic spectrum at 819 nm and 1649 nm, respectively) , was used to mask farmland; the Tasseled Cap (TC), band5 and NDVI were used to mask forestland; the modified normalized differential water index (MNDWI) was used to mask water areas; and the TC (the first band of the TC) and MNDWI were used to mask salinized land. Next, the LSMM algorithm was applied to the masked image, and a grading desertification map was produced from the bare sand fraction map using the decision-tree method based on the grading system (see Supplementary Figure S2 for details). Spatio-temporal change detection. To reveal the mechanisms responsible for grassland desertification spatio-temporal dynamic changes, raster calculation available in the ArcGIS and change detection and decision tree available in ENVI were performed. Before this perform, the farmland, forestland, human settlement land, water, salinized land, and bare land classes were merged as non-grassland. Then, the following criteria were used: an area with an increase in the degree of desertification is designated as ‘Deteriorated’ (e.g., a change from SlDG to MDG); an area with a cross-level increase in the degree of desertification is designated as ‘Seriously Deteriorated’ (e.g., a change from SlDG to SeDG). Similarly, the ‘Reversed’ and ‘Significantly Reversed’ classifications were Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 9 www.nature.com/scientificreports/ defined in the same way. An area without any changes in the degree of desertification is referred to as ‘stable’, except where a change from grassland to non-grassland or vice versa occurs. Field survey. To develop interpretation symbols, a grading system consistent with conditions in the HSL, and field truth data as references for validation, a field survey was conducted on 5–20 August, 2013. The dates of the field survey nearly coincided with of the dates of the Landsat images, which ensured the reliability of the interpre- tation symbols, grading system and result validation. In this article, 125 plots were selected that contained herbs and low shrubs (measuring 1 × 1 m) or shrubs (measuring 10 × 10 m). Plots were georeferenced using GPS, and the vegetation canopy was orthogonally projected to the ground to estimate the vegetation coverage. Aer ex ft clud- ing the plots used for end-member selection and for interpretation symbols development, 81 plots remained for validation. The accuracy was estimated from a scatter plot correction to compare the field vegetation cover in each plot with the pixel values extracted from the vegetation fraction image (derived from LSMM) by using the extract-value-to-points function via the ArcGIS software (Fig. 6). References 1. UNCCD. United Nations convention to combat desertic fi ation in those countries experiencing serious drought and/or desertic fi ation, particularly in Africa (Geneva, 1994). 2. Reynolds, J. F. et al. Global desertification: building a science for dryland development. Science 316, 847–851 (2007). 3. Reynolds, J. F. et al. Scientific concepts for an integrated analysis of desertification. Land Degrad. Dev. 22, 166–183 (2011). 4. Grainger, A. Is land degradation neutrality feasible in dry areas? J. Arid Environ. 112, 14–24 (2015). 5. Badreldin, N. & Goossens, R. A satellite-based disturbance index algorithm for monitoring mitigation strategies effects on desertification change in an arid environment. Mitig. Adapt. Strateg. Glob. Change 20, 263–276 (2015). 6. Gnacadja, L. New challenges in science and policies to combat desertification. J. Arid Environ. 112, 1–4 (2015). 7. Karnieli, A., Qin, Z., Wu, B., Panov, N. & Yan, F. Spatio-temporal dynamics of land-use and land-cover in the Mu Us sandy land, China, using the change vector analysis technique. Remote Sensing 6, 9316–9339 (2014). 8. Dawelbait, M. & Morari, F. Limits and potentialities of studying dryland vegetation using the optical remote sensing. Ital. J. Agron. 3, 97–106 (2010). 9. Wessels, K. J., van den Bergh, F. & Scholes, R. J. Limits to detectability of land degradation by trend analysis of vegetation index data. Remote Sens. Environ. 125, 10–22 (2012). 10. Higginbottom, T. P. & Symeonakis, E. Assessing land degradation and desertification using vegetation index data: current frameworks and future directions. Remote Sensing 6, 9552–9575 (2014). 11. Paz-Kagan, T., Panov, N., Shachak, M., Zaady, E. & Karnieli, A. Structural changes of desertified and managed shrubland landscapes in response to drought: spectral, spatial and temporal analyses. Remote Sensing 6, 8134–8164 (2014). 12. Zhang, X., Shang, K., Cen, Y., Shuai, T. & Sun, Y. Estimating ecological indicators of karst rocky desertification by linear spectral unmixing method. Int. J. Appl. Earth Obs. 31, 86–94 (2014). 13. Wang, T. et al. Spatial-temporal changes of sandy desertified land during last 5 decades in northern China. Acta Geogr. Sin. 59, 203–212 (2004). 14. Sweeney, S., Ruseva, T., Estes, L. & Evans, T. Mapping cropland in smallholder-dominated savannas: integrating remote sensing techniques and probabilistic modeling. Remote Sensing 7, 15295–15317 (2015). 15. Somers, B. & Asner, G. P. Tree species mapping in tropical forests using multi-temporal imaging spectroscopy: wavelength adaptive spectral mixture analysis. Int. J. Appl. Earth Obs. 31, 57–66 (2014). 16. Thorp, K. R., French, A. N. & Rango, A. Effect of image spatial and spectral characteristics on mapping semi-arid rangeland vegetation using multiple endmember spectral mixture analysis (MESMA). Remote Sens. Environ. 132, 120–130 (2013). 17. Allard, V. et al. e r Th ole of grazing management for the net biome productivity and greenhouse gas budget (CO , N O and CH ) of 2 2 4 semi-natural grassland. Agri., Ecosyst. Environ. 121, 47–58 (2007). 18. Li, J. Y. et al. Spatiotemporal variations in grassland desertification based on landsat images and spectral mixture analysis in Yanchi County of Ningxia, China. IEEE J. Sel. Top. Appl. Earth Observations Remote Sensing 7, 4393–4402 (2014). 19. Verstraete, M. M., Brink, A. B., Scholes, R. J., Beniston, M. & Stao ff rd Smith, M. Climate change and desertification: where do we stand, where should we go? Glob. Planet. Change 64, 105–110 (2008). 20. Han, Z. et al. Change trends for desertified lands in the Horqin sandy land at the beginning of the twenty-first century. Environ. Earth Sci. 59, 1749–1757 (2010). 21. Zhao, H. L. Desertification processes and its restoration mechanisms in the Horqin sand land. (Ocean Press, 2003). 22. Miao, R. et al. Effectiveness of shrub planting and grazing exclusion on degraded sandy grassland restoration in Horqin sandy land in Inner Mongolia. Ecol. Eng. 74, 164–173 (2015). 23. Shen, Z. Naiman Banner has adopted a variety of mechanisms and models to combat desertification. J. Inner Mongolia Forest 1, 10 (2014). 24. Li, X., Wang, H., Zhou, S., Sun, B. & Gao, Z. Did ecological engineering projects have a significant effect on large-scale vegetation restoration in Beijing-Tianjin sand source region, China? A remote sensing approach. Chin. Geogr. Sci. 26, 216–228 (2016). 25. Ge, X. et al. Quantifying the synergistic effect of the precipitation and land use on sandy desertification at county level: a case study in Naiman Banner, northern China. J. Environ. Manage. 123, 34–41 (2013). 26. Tan, M. & Li, X. Does the green Great Wall effectively decrease dust storm intensity in China? A study based on NOAA NDVI and weather station data. Land Use Policy 43, 42–47 (2015). 27. Wang, T., Xue, X., Zhou, L. & Guo, J. Combating aeolian desertification in northern China. Land Degrad. Develop. 26, 118–132 (2015). 28. Yan, Q., Zhu, J., Zheng, X. & Jin, C. Causal effects of shelter forests and water factors on desertification control during 2000–2010 at the Horqin sandy land region, China. J. For. Res. 26, 33–45 (2015). 29. Zheng, X., Zhu, J. J., Yan, Q. L. & Song, L. N. Effects of land use changes on the groundwater table and the decline of Pinus sylvestris var. Mongolica plantations in southern Horqin sandy land, Northeast China. Agr. Water Manage. 109, 94–106 (2012). 30. Wang, H., Li, W., Yang, Y. & Jiang, D. Groundwater dynamics in plain area in Tongliao, Inner Mongolia. Inner Mongolia Water Resources 4, 111–113 (2011). 31. Tao, S. et al. Rapid loss of lakes on the Mongolian Plateau. Proc. Natl Acad. Sci. USA 112, 2281–2286 (2015). 32. Lu, G. & Liu, Y. The implementation effect of grassland ecological premium policy. J. Arid Land Res. Environ. 29–32 (2015). 33. Ministry of Finance of the P. R. China. A new round of subsidy and incentives system of grassland ecology protection. http://nys. mof.gov.cn. (Accessed: 23 February 2016) (2016). 34. Ministry of Agriculture of the P. R. China. Implementation guidance for a new round of subsidy and incentives system of grassland ecology protection (2016–2020). http://www.moa.gov.cn. (Accessed: 1 March 2016) (2016). 35. Alatorre, L. C. & Beguería, S. Identification of eroded areas using remote sensing in a badlands landscape on marls in the central Spanish Pyrenees. Catena 76, 182–190 (2009). Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 10 www.nature.com/scientificreports/ 36. Lanfredi, M. et al. Early identification of land degradation hotspots in complex bio-geographic regions. Remote Sensing 7, 8154–8179 (2015). 37. Roy, D. P. et al. Landsat-8: science and product vision for terrestrial global change research. Remote Sens. Environ. 145, 154–172 (2014). 38. Mohamed, I. N. L. & Verstraeten, G. Analyzing dune dynamics at the dune-field scale based on multi-temporal analysis of landsat- TM images. Remote Sens. Environ. 119, 105–117 (2012). 39. Wang, T., Yan, C. Z., Song, X. & Xie, J. L. Monitoring recent trends in the area of aeolian desertified land using landsat images in China’s Xinjiang region. ISPRS J. Photogramm. Remote Sens 68, 184–190 (2012). 40. Perkins, T. Speed and accuracy improvements in FLAASH atmospheric correction of hyperspectral imagery. Opt. Eng. 51, 111707 (2012). 41. Song, C., Woodcock, C. E., Seto, K. C., Lenney, M. P. & Macomber, S. A. Classification and change detection using landsat TM data? Remote Sens. Environ. 75, 230–244 (2001). 42. Matthew, M. W. et al. Status of atmospheric correction using a MODTRAN4-based algorithm in Algorithms for multispectral, hyperspectral, and ultraspectral imagery VI (eds Shen, S. S. & Descour, M. R.) 199–207 (SPIE Proceedings, 2000). 43. Wu, B. & Ci, L. J. Landscape change and desertification development in the Mu Us Sandland, Northern China. J. Arid Environ. 50, 429–444 (2002). 44. Sommer, S. et al. Application of indicator systems for monitoring and assessment of desertification from national to global scales. Land Degrad. Dev. 22, 184–197 (2011). 45. Mbow, C., Brandt, M., Ouedraogo, I., de Leeuw, J. & Marshall, M. What four decades of Earth observation tell us about land degradation in the Sahel? Remote Sensing 7, 4048–4067 (2015). 46. Diouf, A. & Lambin, E. F. Monitoring land-cover changes in semi-arid regions: remote sensing data and field observations in the Ferlo, Senegal. J. Arid Environ. 48, 129–148 (2001). 47. Gao, S. W., Wang, B. F., Zhu, L. Y., Wang, J. H. & Zhang, Y. G. Monitoring and evaluation indicator system on sandy desertification of China. Scientia Silvae Sinicae 34, 1–10 (1998). 48. Wang, J. H. & Sun, S. H. Classification of desertification types and its quantification evaluation system. Arid Environ. Monit. 10, 129–137 (1996). 49. Elmore, A. J., Mustard, J. F., Manning, S. J. & Lobell, D. B. Quantifying vegetation change in semiarid environments. Remote Sens. Environ. 73, 87–102 (2000). 50. Dawelbait, M. & Morari, F. Monitoring desertification in a Savannah region in Sudan using landsat images and spectral mixture analysis. J. Arid Environ. 80, 45–55 (2012). 51. Chen, X., Vierling, L., Rowell, E. & DeFelice, T. Using lidar and effective LAI data to evaluate IKONOS and landsat 7 ETM+ vegetation cover estimates in a ponderosa pine forest. Remote Sens. Environ. 91, 14–26 (2004). 52. Collado, A. D., Chuvieco, E. & Camarasa, A. Satellite remote sensing analysis to monitor desertification processes in the crop- rangeland boundary of Argentina. J. Arid Environ. 52, 121–133 (2002). 53. Yang, J., Weisberg, P. J. & Bristow, N. A. Landsat remote sensing approaches for monitoring long-term tree cover dynamics in semi- arid woodlands: comparison of vegetation indices and spectral mixture analysis. Remote Sens. Environ. 119, 62–71 (2012). 54. Li, J. et al. Monitoring and analysis of grassland desertification dynamics using landsat images in Ningxia, China. Remote Sens. Environ. 138, 19–26 (2013). 55. e Th seira, M. A., Thomas, G. & Sannier, C. A. D. An evaluation of spectral mixture modelling applied to a semi-arid environment. Int. J. Remote Sens. 23, 687–700 (2002). 56. Fan, F. & Deng, Y. Enhancing endmember selection in multiple endmember spectral mixture analysis (MESMA) for urban impervious surface area mapping using spectral angle and spectral distance parameters. Int. J. Appl. Earth Obs. 33, 290–301 (2014). 57. Tompkins, S., Mustard, J. F., Pieters, C. M. & Forsyth, D. W. Optimization of endmembers for spectral mixture analysis. Remote Sens. Environ. 59, 472–489 (1997). 58. Sabol, D. E., Adams, J. B. & Smith, M. O. Quantitative subpixel spectral detection of targets in multispectral images. J. Geophys. Res. Planets 97, 2659–2672 (1992). 59. Van de Voorde, T., De Roeck, T. & Canters, F. A comparison of two spectral mixture modelling approaches for impervious surface mapping in urban areas. Int. J. Remote Sens. 30, 4785–4806 (2009). 60. Fernández-Manso, A., Quintano, C. & Roberts, D. Evaluation of potential of multiple endmember spectral mixture analysis (MESMA) for surface coal mining ae ff cted area mapping in die ff rent world forest ecosystems. Remote Sens. Environ. 127, 181–193 (2012). 61. Okin, G. S. Relative spectral mixture analysis — a multitemporal index of total vegetation cover. Remote Sens. Environ. 106, 467–479 (2007). 62. Jiapaer, G., Chen, X. & Bao, A. A comparison of methods for estimating fractional vegetation cover in arid regions. Agr. Forest Meteorol. 151, 1698–1710 (2011). 63. Ghrefat, H. A. & Goodell, P. C. Land cover mapping at Alkali Flat and Lake Lucero, White Sands, New Mexico, USA using multi- temporal and multi-spectral remote sensing data. Int. J. Appl. Earth Obs. 13, 616–625 (2011). 64. Pan, J. & Li, T. Extracting desertification from landsat TM imagery based on spectral mixture analysis and albedo-vegetation feature space. Nat. Haz. 68, 915–927 (2013). 65. Zucca, C., Wu, W., Dessena, L. & Mulas, M. Assessing the effectiveness of land restoration interventions in Dry lands by Multitemporal remote sensing – a case study in ouled DLIM (Marrakech, Morocco). Land Degrad. Develop. 26, 80–91 (2015). Acknowledgements This work was jointly supported by the National Natural Science Foundation of China [NSFC, 31372354] and [NSFC, 41571105]; the Grassland Monitoring and Supervision Centre, Ministry of Agriculture of P.R. China [grant number 425-1]; and the Center of Supervision Management of Grassland, Ministry of Agriculture of P.R. China. Author Contributions J.L. and B.X. designed the study, analysed the data and wrote the manuscript. J.L.,Y.J., F.Z. and J.G. interpreted the images. J.L. and L.Z. validated the results. B.X., X.Y. and Z.Q. reviewed the manuscript and gave comments and suggestions to improve the manuscript. All authors performed the field investigations. Additional Information Supplementary information accompanies this paper at doi:10.1038/s41598-017-03267-x Competing Interests: The authors declare that they have no competing interests. Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 11 www.nature.com/scientificreports/ Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 Scientific Repo R ts | 7: 3009 | DOI:10.1038/s41598-017-03267-x 12
Scientific Reports – Springer Journals
Published: Jun 7, 2017
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
All the latest content is available, no embargo periods.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud