Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Optimal Spectral Domain Selection for Maximizing Archaeological Signatures: Italy Case Studies
Previous Article in Journal
CMOS Image Sensor with a Built-in Lane Detector
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Assessment of Desertification Using Landsat Data on a Regional Scale – A Case Study in the Ordos Plateau, China

1
Collage of Resources and Environmental Sciences, Nanjing Agricultural University, Nanjing 210095, Jiangsu Province, P.R. China
2
Institute of Scientific and Technical Information of China, Beijing 100038, P.R. China
3
State Key Lab of Resources and Environmental Information System, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, P.R. China
*
Author to whom correspondence should be addressed.
Sensors 2009, 9(3), 1738-1753; https://doi.org/10.3390/s90301738
Submission received: 2 December 2008 / Revised: 4 January 2009 / Accepted: 6 February 2009 / Published: 12 March 2009
(This article belongs to the Section Remote Sensors)

Abstract

:
Desertification is a serious threat to the ecological environment and social economy in our world and there is a pressing need to develop a reasonable and reproducible method to assess it at different scales. In this paper, the Ordos Plateau in China was selected as the research region and a quantitative method for desertification assessment was developed by using Landsat MSS and TM/ETM+ data on a regional scale. In this method, NDVI, MSDI and land surface albedo were selected as assessment indicators of desertification to represent land surface conditions from vegetation biomass, landscape pattern and micrometeorology. Based on considering the effects of vegetation type and time of images acquired on assessment indictors, assessing rule sets were built and a decision tree approach was used to assess desertification of Ordos Plateau in 1980, 1990 and 2000. The average overall accuracy of three periods was higher than 90%. The results showed that although some local places of Ordos Plateau experienced an expanding trend of desertification, the trend of desertification of Ordos Plateau was an overall decrease in from 1980 to 2000. By analyzing the causes of desertification processes, it was found that climate change could benefit for the reversion of desertification from 1980 to 1990 at a regional scale and human activities might explain the expansion of desertification in this period; however human conservation activities were the main driving factor that induced the reversion of desertification from 1990 to 2000.

1. Introduction

Desertification is recognized as one of the most serious social-economic-environmental issues in arid, semi-arid and dry sub-humid areas, and in more than one hundred countries, about one billion of the six billion world population are affected by desertification [1,2]. Since the United Nations Conference on Desertification held in 1977, the definition and assessing methodology for desertification have been controversial. Although the definition of desertification put forward by United Nations Convention to Combat Desertification [3] in 1994, that “desertification is land degradation in arid, semiarid and dry sub-humid areas resulting from various factors, including climatic variations and human activities” has been widely accepted, there is no consensus concerning assessment systems and methodology on the proper way to assess it and this may result in desertification being measured by different researchers using different assessment methods that cannot be compared for the same region [4]. Meanwhile, there is a potential problem that qualitative indicators and assessment methods may produce human-induced errors when assessed by different people, which can also reduce accuracy for desertification measurement. Therefore, there is a pressing need to develop a quantitative methodology, especially from a remote sensing point of view, to accurately measure desertification on both regional and global scales.
Climatic variations including rainfall, temperature, sunlight, wind etc. and human activities like intensive land use, overgrazing, over-cutting etc. can lead to changes in land surface conditions during desertification processes [58]. As a consequence of climatic variations and human activities that affect desertification, vegetation conditions in arid, semiarid and dry sub-humid environment have greatly changed, and these changes not only include the reduction of vegetation cover, density and biomass, but also are characterized by structural configurations of vegetation types and landscape patterns [5,9,10]. Vegetation dynamics can also induce the change of soil properties including physical and chemical properties of soil [11,12], and micrometeorological conditions of land surface like temperature, albedo etc. [13], which usually have been used to characterize desertification. These changes of land surface conditions make the spectral characteristics of desertification land vary greatly to different degrees, which can be captured by satellite sensors and this may be fundamental for quantitatively assessing desertification by means of the indices derived from satellite images [14]. Because the accuracy of soil properties acquired from satellite images is very variable [15,16], much previous research has focused on constructing image-based indices to retrieve vegetation and micrometeorological conditions of land surface to monitor desertification at different scales [1719].
Image-based vegetation indices such as Normalized Difference Vegetation Index (NDVI), Soil Adjusted Vegetation Index (SAVI), and Modified SAVI (MSAVI) are virtually always used to assess vegetation conditions in desertification monitoring [18,20]. Meanwhile, some indices that can reflect vegetation patterns and spatial heterogeneity, such as the landscape connectivity index [21], Moving Standard Deviation Index (MSDI) [22,23] etc. are adopted for desertification assessment. There are classical methods for retrieving the micrometeorological conditions of land surfaces such as land surface temperature (LST) and albedo from satellite image by using thermal infrared band or other band combinations to accurately get micrometeorological parameters of land surfaces [20,24]. Based on these indices retrieved from images, various assessment models such as decision tree classification [20], unsupervised classification [25], Radial Basis Function Network (RBFN) [26] etc. were used to determine desertification status. However, many researchers have focused on assessment indicators or methods themselves, but little research has been conducted on the effect of environment heterogeneity and seasonal variation on the values of assessment indicators. Liu divided North China into four sub-regions according to climate variation and built different indices system for each region to assess desertification [20], but to the authors’ knowledge few researchers have made similar attempts on a regional scale.
Landsat series data can provide long-term, high quality multi-band images for environmental monitoring and assessment at regional, national and global scales [27,28]. In this paper, we selected the Ordos Plateau in China as the research region and attempted to develop a quantitative method for desertification assessment by developing indicators from Landsat data that consider the effect of environment heterogeneity and seasonal variation on values determined using the indicator system.

2. Study Area

The Ordos Plateau is located in the southwest of the Inner Mongolia Autonomous Region of China, ranging from 37°41′∼40°51′ N and 106°42′∼111°31′ E (Figure 1). Ordos Plateau (hereinafter shortened to “Ordos”) has an area of approximately 86,752 km2 and population about 1,395,000. It contains an ecotone with air circulation changing from temperate continental climate, high-pressure ridge (Mongolia-Siberian) in the northwest to the influence of monsoonal trough in the southeast. This air circulation makes Ordos have the climate with an annual sunshine duration of 2,716–3,194 hours; an average annual temperature of 5.3–8.7°C, and an annual evaporation of 2,000–3,000mm. The average annual rainfall of Ordos is 170–450 mm, which varies from west to east with most of the annual rainfall falling from July to September. The climate transition in Ordos make its vegetation diverse and mainly includes five vegetation or landscape types: irrigated farmland along Yellow River in the north, where dominated vegetation type is crop such as corn and wheat; temperate steppe in the east, where the plant community coverage can reach 40–50% and the dominant grass species are Stipa bungeana, Aster altaicus and S. breviflora; temperate deciduous scrubs in the southeast, where covers a majority of Mu Us sandy land and the dominant scrub species are Artemisia giraldii and Caragana korshinskii; steppe shrub desert in the southwest, where shrub and semi-shrub species are usually lower and dominant specie is Caragana tibetica; desert in the northwest, where covered by drift dune and the dominant species are Agriophyllum squarrosum, Corispermum spp., salsola spp.. Long-term intensive land use and dry climate condition make desertification become the most serious environmental problem in Ordos, which has Mu Us sandy land (in the southeast) and Qubqi desert (in the northwest) covering approximately 43,000 km2, and accounting for 48% of the whole region [29,30].

3. Materials and Methodology

3.1. Data Source

Landsat data including seven scenes (September 1st, 1979, path/row 137/32; August 19, 1978, path/row 137/33; August 16, 1980, path/row 138/32; October 9, 1980, path/row 138/33; August 20, 1978, path/row 138/34; September 22, 1980, path/row 139/32; October 9, 1979, path/row 139/33) of Landsat MSS (Multi Spectral Scanner, MSS) in 1980 (here 1980 means the period before and after 1980 and the same meaning as that in 1990 and 2000), seven scenes (August 13, 1990, path/row 127/32; September 11, 1989, path/row 127/33; August 7, 1991, path/row 128/32; September 15, 1990, path/row 128/33; August 23, 1991, path/row 128/34; August 24, 1989, path/row 129/32; August 30, 1991, path/row 129/33) of Landsat TM (Thematic Mapper, TM) in 1990 and seven scenes (June 13, 2000, path/row 127/32; August 6, 2000, path/row 127/33; November 11, 2000, path/row 128/32; November 11, 2000, path/row 128/33; September 22, 1999, path/row 128/34; June 11, 2000, path/row 129/32; August 12, 1999, path/row 129/33) of Landsat ETM (Enhanced Thematic Mapper plus, ETM+) in 2000 covering all of Ordos were used for monitoring desertification from 1980 to 2000. All spectral data were collected in summer and autumn, mainly from June to November. Radiometric calibration, geometric correction, and cloud removal were carried out for all of these images, which were all geo-referenced to the WGS_1984_UTM Projected Coordinate System with a geometric precision of 0.5 pixel. A 1:1,000,000 vegetation map of Ordos was acquired from the 1:1,000,000 vegetation map of China, which was used for vegetation and landscape classification in Ordos. Meteorological data including monthly averaged temperature, rainfall, etc. of thirty-four meteorological stations in and around Ordos were collected from China Central Meteorological Bureau. Social-economic census data of Ordos were collected from local statistical yearbooks.

3.2. Indicator Selection and Acquisition

The changes of land surface conditions induced by desertification are mainly characterized by the change of vegetation biomass or cover, landscape pattern and micrometeorological conditions. Therefore, the indicators that reflect these changes can be selected to assess desertification. Considering the strong correlation among indicators in assessing one aspect of land surface condition, such as among NDVI, SAVI and vegetation cover in assessing vegetation biomass [31], we only selected one indicator to represent one aspect of land surface conditions. In this study, NDVI, MSDI and land surface albedo were used to characterize vegetation biomass, landscape pattern and micrometeorological conditions of land surface respectively.
In this study, all of the three indicators were acquired through analysis of Landsat images for desertification assessment. To quantify vegetation biomass for desertification, NDVI was selected and calculated as the ratio of the difference between near infrared (NIR) and red (RED) bands to the sum of them (Formula 1):
NDVI = NIR RED NIR + RED
Desertification is characterized not only by the change of vegetation biomass, but also by the change of landscape pattern or heterogeneity. In this study, landscape pattern or heterogeneity in desertification was assessed by using MSDI, which had been proven to be a robust indicator to assess landscape heterogeneity for land degradation. The MSDI was calculated by passing a 3×3 filter across red wavelengths band (MSS band2 and TM/ETM+ band3) and calculating the standard deviation for every nine-pixel window. Then the standard deviation was placed onto a new map at the same location as the central pixel of each nine-pixel window (Formula 2) [22,23]. Where N represent pixel number of filter, and here N =9; DNi is the digital number of pixel i in each nine-pixel window and DN ¯ is the average digital number value of each nine-pixel window:
MSDI = i = 1 N ( DN i DN ¯ ) 2 N
Land surface albedo is an important indicator that determines energy budget and the change of micrometeorological conditions like temperature, aridity/humidity etc. of land affected by desertification [13,32]. Some studies showed the increasing of land surface albedo implies a degradation of land quality [33]. In this study, we adopted broad-band albedo that determined by the combination of narrow-band albedo to assess the micrometeorological conditions of land surface. Narrow-band albedo for each band (except band 6 for Landsat TM/ETM+) of Landsat images was determined by using the dark object subtraction method technique [34,35]. Broad-band albedo was then calculated according to its relationship with each narrow-band [36,37]. The broad-band albedo calculation for Landsat MSS and Landsat TM/ETM+ are expressed by Formulas (3) and (4):
Albedo MSS = 0.228 α 1 + 0.217 α 2 + 0.182 α 3 + 0.373 α 4
Albedo TM / ETM + = 0.356 α 1 + 0.130 α 3 + 0.373 α 4 + 0.085 α 5 + 0.072 α 7 0.0018
where AlbedoMSS and AlbedoTM / ETM+ represent the broad-band albedo of Landsat MSS and Landsat TM/ETM+; α1–7 are narrow-band albedo for each band of Landsat MSS (band 1–4) and Landsat TM/ETM+.

3.3. Assessment Method

In this study, the land desertification of Ordos was classified into five grades: non-desertification, low desertification, medium desertification, high desertification and severe desertification according to a frequently-used grading system in China [38,39]. Our goal was to separate each desertification grade by using the differences among the indicators and its combinations of desertification land. Decision Tree (DT) is a suitable approach for classification because it can employ tree structured rules to recursively partition the data set feature spaces into increasingly homogenous classes based on a splitting criterion. DT classification has been extensively used for vegetation mapping, ecological modeling, soil mapping, and in remote sensing studies [40,41]. In this study, we used DT classification to assess the status of desertification of Ordos.
Information from ancillary data sources may help to increase assessment accuracy, especially in large area land-cover mappings that involve cover types characterized by high levels of within-class variability [42]. Land surface vegetation conditions in Ordos have great heterogeneity due to the rainfall and temperature gradients in this region [30], which may produce major classification errors when using the same rules for the whole region. Therefore, vegetation type was used as ancillary information for desertification assessment. Meanwhile, the date of image may also be served as ancillary data due to variations of indicators for a desertification grade at different times in a year.
According to the desertification classification system used in the National Key Basic Project “Process and Controlling of Aeolian Desertification in North China” [38,39], ground information has been used to define the training and accuracy checking sets, which would be used to make rule sets for assessment and accuracy checking. Fifty training sets and ten checking sets were selected for each grade of desertification in each vegetation sub-region (Figure 1) according to the date of the images. In order to represent more information, these sets must be equally distributed in a specific landscape. In this study, the temporal unit of images used to make rule sets was month and mainly between June and November. The surface information of desertification land used to define the training sets was extracted from the three indicator bands for all images. According to frequencies distribution of the five desertification grade for each vegetation sub-region in a certain period, the rules for each indicator can be established and the DT can be built. Images were classified in ERDAS by using the knowledge engineer module with total agreement of rules in assessment. The single assessed map of whole Ordos for each of the three periods was created by using mosaic images tool in ERDSA to merge all classified results of each period with selecting maximum function in overlap areas. Then accuracy was assessed and the training sets and rule sets might be adjusted according to visual examination and accuracy assessment until better assessment accuracy was achieved.

4. Results and Analysis

4.1. Desertification Assessment and Accuracy Checking

The rule sets of indicators for each grade of desertification in every vegetation sub-region are listed in Tables 1 and 2 (considering the huge size of the total rule sets, only rule sets for August, October and November for irrigated farmland, temperate deciduous scrubs and desert are listed). Great differences exist for defining a desertification grade from the three indictors at a regional scale due to various vegetation sub-regions, the date of images and different satellite sensors. According to the rule sets, the change of NDVI and albedo had a linear relationship with the change of desertification grade, which showed that the reversion of desertification was characterized by the increasing of NDVI or the decreasing of albedo. In this study, water body was classified into non-desertification and can be accurately classified by setting the rule of that region with negative value of NDVI. MSDI was used in this study to characterize desertification from the difference of landscape pattern or heterogeneity in different desertification grades. However, the change of landscape pattern or heterogeneity did not have a linear relationship with the change of desertification. The change of landscape pattern or heterogeneity in desertification land always had a trend like a parabola, which indicated that the MSDI of desertification land initially increased and then declined with the expansion of desertification. In this study, because severe and non-desertification can be accurately identified by using only NDVI and albedo, MSDI was mainly used to distinguish the low, medium and high desertification. The rule sets for different desertification grades varied greatly due to seasonal variation, especially for NDVI and albedo, which had a higher NDVI value and lower albedo value in summer (August) than that in autumn (October and November) respectively. According to the rule sets, Decision Tree was built to assess the status of desertification of Ordos. The Decision Tree model for desertification assessment in temperate deciduous scrubs sub-region in August by using Landsat TM/ETM+ was shown in Figure 2, where landscape 3 is the temperate deciduous scrubs sub-region. Models for desertification assessment in other sub-regions were developed in a similar manner as in Figure 2 and are not presented here.
Accuracy checking of assessing results for 1980, 1990 and 2000 had been carried out by using confusion matrix. A confusion matrix for accuracy checking of desertification assessment results in the three periods was listed in Table 3. The satisfactory classification results have been obtained with the average overall accuracy of three periods higher than 0.9 and the highest overall accuracy and kappa statistic in the assessment for 1990, which reached 0.926 and 0.9075 respectively. Due to the relative low quality of images in 2000, the overall accuracy and kappa statistic were the lowest, which reached 0.9 and 0.875 respectively. The overall accuracy and kappa statistic in the assessment for 1980 were little higher than that of 2000. For the accuracy of each grade, the severe desertification grade was total correctly assessed and this might attribute to that severe desertification land covered by drifting sand dunes, which was characterized by the lowest NDVI and highest albedo that made these regions easily and accurately classified. The land with a medium desertification grade always had the most complex land surface condition, which made some parts of these regions classified into low and high desertification grade and the accuracy of medium desertification having the lowest, albeit still good, accuracy.

4.2. The Processes and Causes of Desertification from 1980 to 2000

The assessment results of desertification in 1980, 1990 and 2000 are shown in Figure 3. Owing to the differences of spatial resolution between Landsat MSS (57 m) and TM/ETM+ (28.5 m), the assessment result of desertification in 1980 was resampled to the same resolution as that in 1990 and 2000 for comparison. According to areas of different desertification grades in three periods (Figure 4), the desertification of Ordos displayed a reversing trend, although some local places experienced a desertification expansion from 1980 to 2000. From 1980 to 2000, the area of non-desertification had a significant increase and low desertification had changed little. Based on analyzing the area and percentage of regions experienced desertification reversion and expansion in each sub-region (Table 4), all sub-regions had a reversing trend except steppe shrub sub-region from 1980 to 1990 and the proportion of reversing desertification in irrigated farmland was higher than that of other sub-regions. From 1990 to 2000, desertification in irrigated farmland and desert sub-regions had an expanding trend; however the other three sub-regions experienced a reversing trend with the highest reversing percentage found in the temperate steppes.
The change of desertification is caused by both climate change and human activities. In this study, climate data from thirty-four meteorological stations in and around Ordos and some social-economic indicator were used to analyze the impact of climate change and human activities on desertification from 1980 to 2000. Aridity/humidity conditions of land surface can reflect the integrated impact of various climate factors [43], which can be used to explain the effect of climate change on vegetation growth and the processes of desertification. To investigate the impacts of climate change on desertification, aridity/humidity index that calculated as the ratio of potential evapotranspiration to rainfall [44] was used to reflect the effects of climate change on aridity/humidity condition of land surface and the processes of desertification. In this study, climate data from thirty-four meteorological stations in and around Ordos were interpolated to grid data with 1 km resolution and a linear trend of aridity/humidity index in every pixel was calculated to analyze the impacts of climate change in Ordos from a spatial point of view [Figures 5(a) and (b)].
Meanwhile, the trend of regional averaged rainfall and temperature were also illustrated [Figures 5(c) and (d)]. Due to the social-economic data are always collected based on administrative unit, so the sum of social-economic indicator data (livestock number and area of afforestation) of all counties in Ordos was used in temporal serial analysis to reflect the impacts of human activities on desertification (Figure 6). From 1980 to 1990, although the temperature had increased somewhat, the rainfall had a significant increase and the aridity/humidity index of the whole region showed a decreasing trend, which indicated the climate of this period was more humid and would benefit vegetation growth and the reversion of desertification. The change of human activities in this period showed that although livestock number decreased somewhat, the area of afforestation decreased significant especially in the late of 1980s; meanwhile, related vegetation protection policies were not extensively conducted in this period and land-use was still intensive. Based on integrated analysis of the trend and characteristics of climate change and human activities, human activities might explain the expansion of desertification in Ordos from 1980 to 1990 at a regional scale. From 1990 to 2000, the temperature showed an increasing trend and the rainfall had a decreasing trend, which made the aridity/humidity index of the whole region show an increasing trend and it indicated the climate of this period was more arid and would speed up the expansion of desertification. Meanwhile, livestock numbers in Ordos decreased gradually and the area of afforestation had an obviously increase in this period, which might attribute to ecological projects such as Three-North Forest Shelterbelt Program conducted in Ordos since 1990s [45,46]. Based on integrated analysis of the trend and characteristics of climate change and human activities, human activities might explain the reversion of desertification in Ordos from 1980 to1990 at a regional scale.

5. Conclusions

This paper focused on assessing desertification from a quantitative viewpoint by means of remote sensing. In this study, assessing indicators were selected according to three different aspects from the responses of land surface conditions to desertification. Based on considering the differences of indicators for assessing desertification in different region and time, vegetation or landscape types, and dates of image acquisition were employed for building rule sets in a Decision Tree model at a regional scale, which improved the accuracy of assessment. In this study, average overall accuracy of assessments in three periods was higher than 90%.
Although some local places of Ordos experienced an expansion of desertification, the trend of desertification of Ordos was an overall decrease in from 1980 to 2000. For the change of desertification in sub-regions, temperate steppe and temperate deciduous scrubs sub-regions showed a stable reversion of desertification from 1980 to 2000; however, the other sub-regions did not show a stable reversion trend in this 20-year period. The reversing trend of desertification from 1980 to 1990 was mainly induced by climate change; however human activity was the factor that dominated the reversion of desertification from 1990 to 2000 at a regional scale.
The concepts and methods put forward in this study not only can quantitatively assess desertification at a regional scale, but also can be used in other land degradation assessment at different temporal and spatial scales. Improvement of the indicator system and use of more ancillary data may be helpful to improve the assessing accuracy.

Acknowledgments

This research was supported and funded by the Knowledge Innovation Project of The Chinese Academy of Sciences (Grant NO.KSCX2-YW-G-036), a special program of China Earthquake Administration (Grant NO.200708007-03) and the Chinese National Academic Institutional Commonweal Research Programs (Grant NO.2005DIA3J006). We are also very grateful to Wenya Li from Nanjing Normal University for correcting the English of this manuscript.

References and Notes

  1. UNSO Office to Combat Desertification and Drought. Aridity zones and dryland populations: An assessment of population levels in the world’s drylands; UNSO/UNDP: New York, 1997. [Google Scholar]
  2. Thomas, D.S.G. Science and desertification debate. J. Arid. Environ 1997, 37, 599–608. [Google Scholar]
  3. UNCCD. United Nations convention to combat desertification in countries experiencing serious drought and/or desertification, particularly in Africa; A/AC.247/27,; Pairs, 1994; p. 4. [Google Scholar]
  4. Veron, S.R.; Paruelo, J.M.; Oesterheld, M. Assessing desertification. J. Arid. Environ 2006, 66, 751–763. [Google Scholar]
  5. Sivakumar, M.V.K. Interactions between climate and desertification. Agric. For. Meteorol 2007, 142, 143–155. [Google Scholar]
  6. Oba, G.; Post, E.; Stenseth, N. Sub-saharan desertification and productivity are linked to hemispheric climate variability. Glob. Change. Biol 2001, 7, 241–246. [Google Scholar]
  7. Javier, J.; Martinez, J.; Schnabel, S. Desertification due to overgrazing in a dynamic commercial livestock-grass-soil system. Ecol. Model 2007, 205, 277–288. [Google Scholar]
  8. Zeidler, J.; Hanrahan, S.; Scholes, M. Land-use intensity affects range condition in arid to semi-arid Namibia. J. Arid. Environ 2002, 52, 389–403. [Google Scholar]
  9. Li, X.R.; Jia, X.H.; Dong, G.R. Influence of desertification on vegetation pattern variations in the cold semi-arid grasslands of Qinghai-Tibet Plateau, North-west China. J. Arid. Environ 2006, 64, 505–522. [Google Scholar]
  10. Wu, B.; Ci, L.J. Landscape change and desertification development in the Mu Us Sandland, Northern China. J. Arid. Environ 2002, 50, 429–444. [Google Scholar]
  11. Zhao, H.L.; Zhou, R.L.; Zhang, T.H.; Zhao, X.Y. Effects of desertification on soil and crop growth properties in Horqin sandy cropland of Inner Mongolia, north China. Soil. Till. Res 2006, 87, 175–185. [Google Scholar]
  12. Su, Y.Z.; Zhao, W.Z.; Su, P.X.; Zhang, Z.H.; Wang, T.; Ram, R. Ecological effects of desertification control and desertified land reclamation in an oasis-desert ecotone in an arid region: A case study in Hexi Corridor, northwest China. Ecol. Eng 2007, 29, 117–124. [Google Scholar]
  13. Li, S.G.; Harazono, Y.; Oikawa, T.; Zhao, H.L.; He, Z.Y.; Chang, X.L. Grassland desertification by grazing and the resulting micrometeorological changes in Inner Mongolia. Agr. Forest. Meteorol 2000, 102, 125–137. [Google Scholar]
  14. Pinet, P.C.; Kaufmann, C.; Hill, J. Imaging spectroscopy of changing Earth’s surface: a major step toward the quantitative monitoring of land degradation and desertification. C. R. Geosci 2006, 338, 1042–1048. [Google Scholar]
  15. Udelhoven, T.; Emmerling, C.; Jarmer, T. Quantitative analysis of soil chemical properties with diffuse reflectance spectrometry and partial least-square regression: A feasibility study. Plant Soil 2003, 251, 319–329. [Google Scholar]
  16. Lagacherie, P.; Baret, F.; Feret, J.B.; Netto, J.M.; Robbez-Masson, J.M. Estimation of soil clay and calcium carbonate using laboratory, field and airborne hyperspectral measurements. Remote Sens. Envrion 2008, 112, 825–835. [Google Scholar]
  17. Runnstrom, M.C. Rangeland development of the Mu Us sandy land in semiarid China: An analysis using Landsat and NOAA remote sensing data. Land Degrad. Dev 2003, 14, 189–202. [Google Scholar]
  18. Geerken, R.; Ilaiwi, M. Assessment of rangeland degradation and development of a strategy for rehabilitation. Remote Sens. Envrion 2004, 90, 490–504. [Google Scholar]
  19. Sun, D.F.; Dawson, R.; Li, H.; Li, B.G. Modeling desertification change in Minqin county, China. Environ. Monit. Assess 2005, 108, 169–188. [Google Scholar]
  20. Liu, A.X.; Wang, J.; Liu, Z.J.; Wang, J. Monitoring desertification in arid and semi-arid areas of China with NOAA-AVHRR and MODIS data. Geoscience and Remote Sensing Symposium, 2005. IGARSS apos; 05. Proceedings. 2005 IEEE international 2005, 4, 2362–2364. [Google Scholar]
  21. Sun, D.F.; Dawson, R.; Li, H.; Wei, R.; Li, B.G. A landscape connectivity index for assessing desertification: a case study of Minqin County, China. Landscape Ecol 2007, 2, 531–543. [Google Scholar]
  22. Tanser, F.C.; Palmer, A.R. The application of a remotely-sensed diversity index to monitor degradation patterns in a semi-arid, heterogeneous, South African landscape. J. Arid. Environ 1999, 43, 477–484. [Google Scholar]
  23. Jafari, R.; Lewis, M.M.; Ostendorf, B. An image-based diversity index for assessing land degradation in an arid environment in South Australia. J. Arid. Environ 2008, 72, 1282–1293. [Google Scholar]
  24. Julien, Y.; Sobrino, J.A.; Verhoef, W. Changes in land surface temperatures and NDVI values over Europe between 1982 and 1999. Remote Sens. Envrion 2006, 103, 43–55. [Google Scholar]
  25. Lira, J. A model of desertification process in a semi-arid environment employing multi-spectral images. In Progress in Pattern Recognition, Image Analysis and Applicaton, 1st Ed; Sanfeliu, A., Martinez-Trinidad, J.F., Carrasco-Ochoa, J.A., Eds.; Springer: Berlin Heidelberg, Germany, 2004; Volume 3287, pp. 249–258. [Google Scholar]
  26. Huang, Q.H.; Cai, Y.L. Assessment of karst rocky desertification using radial basis function network model and GIS technique: a case study of Guizhou Province, China. Environ. Geol 2006, 49, 1173–1179. [Google Scholar]
  27. Qi, S.Z.; Cai, Y.M. Mapping and assessment of degraded land in the Heihe River basin, arid northwestern China. Sensors 2007, 7, 2565–2578. [Google Scholar]
  28. Wulder, M.A.; White, J.C.; Goward, S.N.; Masek, J.G.; Irons, J.R.; Herold, M.; Cohen, W.B.; Loveland, T.R.; Woodcock, C.E. Landsat continuity: Issues and opportunities for landcover monitoring. Remote Sens. Envrion 2008, 112, 955–969. [Google Scholar]
  29. Runnstrom, M.C. Is northern China winning the battle against desertification? Satellite remote sensing as a tool to study biomass trends on the Ordos Plateau in Semiarid China. Ambio 2000, 29, 468–476. [Google Scholar]
  30. Zheng, Y.R.; Xie, Z.X.; Jiang, L.H.; Shimizu, H.; Rimmington, G.M.; Zhou, G.S. Vegetation responses along environmental gradients on the Ordos plateau, China. Ecol. Res 2006, 21, 396–404. [Google Scholar]
  31. Bechtel, A.; Puttmann, W.; Carlson, T.N.; Ripley, D.A. On the relation between NDVI, Fractional Vegetation Cover and Leaf Area Index. Remote Sens. Envrion 1997, 62, 241–252. [Google Scholar]
  32. Jackson, R.D.; Idso, S.B.; Otterman, J. Surface albedo and desertification. Science 1975, 189, 1012–1015. [Google Scholar]
  33. Riobinove, C.J.; Chavez, P.S.; Gehring, D.; Holmgren, R. Arid land monitoring using Landsat albedo difference images. Remote Sens. Envrion 1981, 11, 133–156. [Google Scholar]
  34. Chavez, P.S., Jr. An improved dark-object subtraction technique for atmospheric scattering correction of multispectral data. Remote Sens. Envrion 1988, 24, 459–479. [Google Scholar]
  35. Chavez, P.S., Jr. Image-based atmospheric corrections–Revisited and improved. Photogramm. Eng. Remote Sens 1996, 62, 1025–1036. [Google Scholar]
  36. Musick, H.B. Temporal change of Landsat MSS albedo estimates in arid rangeland. Remote Sens. Envrion 1986, 20, 107–120. [Google Scholar]
  37. Liang, S.L. Narrowband to broadband conversions of land surface albedo I Algorithms. Remote Sens. Envrion 2000, 76, 213–238. [Google Scholar]
  38. Wang, T.; Zhu, Z.D.; Wu, W. Sandy desertification in the north of China. Sci. China Ser. D 2002, 45, 23–34. [Google Scholar]
  39. Wang, T.; Wu, W.; Xue, X.; Sun, Q.W.; Chen, G.T. Study of spatial distribution of sandy desertification in North China in recent 10 years. Sci. China Ser. D 2004, 47, 78–88. [Google Scholar]
  40. Friedl, M.A.; Brodley, C. Decision tree classification of land cover from remotely sensed data. Remote Sens. Envrion 1997, 61, 399–409. [Google Scholar]
  41. Rogan, J.; Franklin, J.; Roberts, D.A. A comparison of methods for monitoring multitemporal vegetation change using Thematic Mapper imagery. Remote Sens. Envrion 2002, 80, 143–156. [Google Scholar]
  42. Mclver, D.K.; Friedl, M.A. Using prior probabilities in decision-tree classification of remotely sensed data. Remote Sens. Envrion 2002, 81, 253–261. [Google Scholar]
  43. Wu, S.H.; Yin, Y.H.; Zheng, D.; Yang, Q.Y. Aridity/humidity status of land surface in China during the last three decades. Sci. China Se. r D 2005, 48, 1510–1518. [Google Scholar]
  44. Yin, Y.H.; Wu, S.H.; Zheng, D.; Yang, Q.Y. Regional difference of aridity/humidity conditions change over China during the last thirty years. Chin. Sci. Bull 2005, 50, 2226–2233. [Google Scholar]
  45. Wu, W. Study on process of desertification in Mu Us sandy land for last 50 years, China. J. Desert Res 2001, 21, 164–169. [Google Scholar]
  46. Wang, B.; Ding, G.D.; Gu, X.H.; Ma, S.L. Initial study on vegetation recovering effects in Maowusu Hinterland – take Wushen Banner of Inner Mongolia as an example. Res Soil Water Conserv 2007, 14, 237–242. [Google Scholar]
Figure 1. General location of study area: Ordos Plateau.
Figure 1. General location of study area: Ordos Plateau.
Sensors 09 01738f1
Figure 2. Decision Tree for desertification assessment of temperate deciduous scrubs sub-region (landscape = 3) in August, 1990. “non”, “low”, “medium”, “high” and “severe” are used to represent desertification grades for short. NDVI, MSDI and albedo are the indicators used to build the rules.
Figure 2. Decision Tree for desertification assessment of temperate deciduous scrubs sub-region (landscape = 3) in August, 1990. “non”, “low”, “medium”, “high” and “severe” are used to represent desertification grades for short. NDVI, MSDI and albedo are the indicators used to build the rules.
Sensors 09 01738f2
Figure 3. Desertification maps of Ordos in 1980 (a), 1990 (b) and 2000 (c).
Figure 3. Desertification maps of Ordos in 1980 (a), 1990 (b) and 2000 (c).
Sensors 09 01738f3
Figure 4. Statistics of the desertification area of Ordos in 1980, 1990 and 2000. “non”, “low”, “medium”, “high” and “severe” are used to represent desertification grades for short.
Figure 4. Statistics of the desertification area of Ordos in 1980, 1990 and 2000. “non”, “low”, “medium”, “high” and “severe” are used to represent desertification grades for short.
Sensors 09 01738f4
Figure 5. The linear trend of aridity/humidity index in Ordos from 1980 to 1990 (a), from 1990 to 2000 (b); the change of rainfall (c) and temperature (d) in Ordos from 1980 to 2000.
Figure 5. The linear trend of aridity/humidity index in Ordos from 1980 to 1990 (a), from 1990 to 2000 (b); the change of rainfall (c) and temperature (d) in Ordos from 1980 to 2000.
Sensors 09 01738f5
Figure 6. The change of livestock number (a) and the area of afforestation (b) in Ordos from 1980 to 2000.
Figure 6. The change of livestock number (a) and the area of afforestation (b) in Ordos from 1980 to 2000.
Sensors 09 01738f6
Table 1. The rule sets for desertification assessment by using Landsat MSS images.
Table 1. The rule sets for desertification assessment by using Landsat MSS images.
MonthDesertification GradeIrrigated farmlandTemperate deciduous scrubsDesert

NDVIMSDIAlbedoNDVIMSDIAlbedoNDVIMSDIAlbedo

Augustnon<0 or >0.45>10–0.22<0 or >0.45>10–0.19<0.2 or >0.4>20–0.24
low0.35–0.450–40.22–0.240.3–0.450–40.19–0.210.28–0.40–80.24–0.26
medium0.23–0.350–60.24–0.260.2–0.30–70.21–0.230.18–0.280–60.26–0.28
high0.1–0.230–50.26–0.280.1–0.20–50.23–0.260.06–0.180–40.28–0.30
severe0–0.10–3>0.280–0.1>1>0.260–0.060–3>0.30
Octobernon<0 or >0.3>10–0.26<0 or >0.35>10–0.24<0.15 or >0.3>10–0.26
low0.2–0.30–50.26–0.290.22–0.350–30.24–0.270.2–0.30–50.26–0.29
medium0.145–0.200–80.29–0.330.16–0.220–50.27–0.30.14–0.20–40.29–0.31
high0.085–0.1450–60.33–0.360.08–0.160–40.3–0.330.04–0.140–30.31–0.34
severe0–0.0850–4>0.360–0.08>1>0.330–0.040–2>0.34
“non”, “low”, “medium”, “high” and “severe” are used to represent desertification grades for short
Table 2. The rule sets for desertification assessment by using Landsat TM/ETM+ images.
Table 2. The rule sets for desertification assessment by using Landsat TM/ETM+ images.
MonthDesertification GradeIrrigated farmlandTemperate deciduous scrubsDesert

NDVIMSDIAlbedoNDVIMSDIAlbedoNDVIMSDIAlbedo

Augustnon<0 or >0.5>10–0.16<0 or >0.5>10–0.175<0.28 or >0.4>10–0.2
low0.4–0.50–30.16–0.180.4–0.50–30.175–0.190.32–0.40–50.16–0.18
medium0.32–0.40–60.18–0.200.32–0.40–50.19–0.2050.26–0.320–40.18–0.20
high0.24–0.320–40.20–0.220.25–0.320–40.205–0.220.22–0.260–30.20–0.22
severe0–0.240–3>0.220–0.25>1>0.220–0.220–2>0.22
Novembernon<0 or >0.27>10–0.35<0 or >0.25>10–0.35<0.18 or >0.24>30–0.37
low0.22–0.270–80.35–0.40.22–0.250–80.35–0.40.20–0.240–80.37–0.42
medium0.18–0.220–60.4–0.450.19–0.220–60.4–0.450.17–0.200–60.42–0.48
high0.15–0.180–50.45–0.50.16–0.190–40.45–0.50.15–0.170–50.48–0.52
severe0–0.150–4>0.50–0.16>1>0.50–0.150–4>0.52
“non”, “low”, “medium”, “high” and “severe” are used to represent desertification grades for short
Table 3. The confusion matrix and accuracy assessment for desertification in 1980, 1990 and 2000.
Table 3. The confusion matrix and accuracy assessment for desertification in 1980, 1990 and 2000.
YearDesertification GradenonlowmediumhighseveretotalProducers accuracyUsers accuracy
1980non94420010094.00%97.92%
low48853010088.00%88.00%
medium08864010086.00%85.15%
high04687310087.00%88.78%
severe0000100100100.00%95.24%
total9610010198105500

Overall accuracy: 0.91; kappa statistic: 0.8875

1990nonlowmediumhighseveretotal
non94510010094.00%97.92%
low29080010090.00%90.00%
medium05887010088.00%87.13%
high00491510091.00%92.86%
severe0000100100100.00%95.24%
total9610010198105500

Overall accuracy: 0.926; kappa statistic: 0.9075

2000nonlowmediumhighseveretotal
non92341010092.00%95.83%
low28666010086.00%86.00%
medium07858010085.00%84.16%
high03787310087.00%88.78%
severe0000100100100.00%95.24%
total9610010198105500

Overall accuracy: 0.9; kappa statistic: 0.875
Average overall accuracy of three periods: 91.2%; average kappa statistic of three periods: 0.89; “non”, “low”, “medium”, “high” and “severe” are used to represent desertification grades for short
Table 4. Area and percentage of desertification reversion and expansion in sub-regions from 1980 to 2000.
Table 4. Area and percentage of desertification reversion and expansion in sub-regions from 1980 to 2000.
sub-regiondesertification change from 1980 to 1990desertification change from 1990 to 2000

reversedexpandedreversedexpanded

Area (km2)Percentage (%)Area (km2)Percentage (%)Area (km2)Percentage (%)Area (km2)Percentage (%)

irrigated farmland1555.337.8%1100.126.8%1354.232.9%1359.733.1%
temperate steppe3240.234.6%3091.133.0%3708.539.6%2604.927.8%
temperate deciduous scrubs15601.528.0%14594.126.2%20057.536.0%13132.323.6%
steppe shrub1887.920.6%2667.229.1%3066.533.4%2547.327.8%
desert2037.723.7%1816.021.2%1965.822.9%2086.924.3%
total2432.623268.530152.521731.1

Share and Cite

MDPI and ACS Style

Xu, D.; Kang, X.; Qiu, D.; Zhuang, D.; Pan, J. Quantitative Assessment of Desertification Using Landsat Data on a Regional Scale – A Case Study in the Ordos Plateau, China. Sensors 2009, 9, 1738-1753. https://doi.org/10.3390/s90301738

AMA Style

Xu D, Kang X, Qiu D, Zhuang D, Pan J. Quantitative Assessment of Desertification Using Landsat Data on a Regional Scale – A Case Study in the Ordos Plateau, China. Sensors. 2009; 9(3):1738-1753. https://doi.org/10.3390/s90301738

Chicago/Turabian Style

Xu, Duanyang, Xiangwu Kang, Dongsheng Qiu, Dafang Zhuang, and Jianjun Pan. 2009. "Quantitative Assessment of Desertification Using Landsat Data on a Regional Scale – A Case Study in the Ordos Plateau, China" Sensors 9, no. 3: 1738-1753. https://doi.org/10.3390/s90301738

APA Style

Xu, D., Kang, X., Qiu, D., Zhuang, D., & Pan, J. (2009). Quantitative Assessment of Desertification Using Landsat Data on a Regional Scale – A Case Study in the Ordos Plateau, China. Sensors, 9(3), 1738-1753. https://doi.org/10.3390/s90301738

Article Metrics

Back to TopTop