Multi-temporal and multi-source cartography of the glacial cover of Nevado Coropuna (Arequipa,Peru) between 1955 and 2003

This study reports on the glacial cover evolution of the Nevado Coropuna between 1955 and 2003, based on Peruvian topographic maps and satellite images taken from the Landsat 2 and 5 multispectral scanner (MSS), Landsat 5 Thematic Mapper (TM) and Landsat 7 (ETM+). The normalized difference snow index has been applied to these images to estimate the glacierized area of Coropuna. The satellite-based results show that the glacier area was 105 ± 16 km² in 1975 (Landsat 2 MSS), which then reduced to 96 ± 15 km² in 1985 (Landsat 5 MSS), 64 ± 8 km² in 1996 (Landsat 5 TM) and 56 ± 6 km² in 2003 (Landsat 5 TM). Altogether, between 1955 and 2003, Coropuna lost 66 km² of its glacial cover, which represents a mean retreat of 1.4 km² year^?1, that is, a loss of 54% in 48 years (11% loss per decade). The maximum rate of retreat occurred during the 1980s and 1990s, a phenomenon probably linked with the pluviometric deficit of El Niño events of 1983 and 1992.     

Frontal changes in the Manimahesh and Tal Glaciers in the Ravi basin,Himachal Pradesh,northwestern Himalaya (India), between 1971 and 2013

The Manimahesh and Tal Glaciers are located in the Budhil fifth-order sub-basin of the Ravi, Himachal Himalaya, Northwestern Himalaya (India). These glaciers were analysed using high- (Corona KH-4A) to medium- (Landsat TM/ETM+/OLI, ASTER) spatial resolution satellite data between 1971 and 2013, along with extensive field measurements (2011–2014) of frontal changes. The results show that the Manimahesh and Tal Glaciers retreated by 157 ± 34 m (4 ± 1 m year^–1) and 45 ± 34 m (1 ± 1 m year^–1), respectively, whereas, the total area lost is estimated at 0.21 ± 0.01 km² (0.005 km² year^–1) and 0.010 ± 0.003 km² (0.0002 km² year^–1), respectively, between 1971 and 2013. The rate of retreat is significantly lower than that previously reported. Our field measurements (2011–2014) also suggest a retreating trend and validate the measured glacier changes using remotely sensed temporal data.     

Production of Landsat ETM+ reference imagery of burned areas within Southern African savannahs: comparison of methods and application to MODIS

Accurate production of regional burned area maps are necessary to reduce uncertainty in emission estimates from African savannah fires. Numerous methods have been developed that map burned and unburned surfaces. These methods are typically applied to coarse spatial resolution (1 km) data to produce regional estimates of the area burned, while higher spatial resolution (<30 m) data are used to assess their accuracy with little regard to the accuracy of the higher spatial resolution reference data. In this study we aimed to investigate whether Landsat Enhanced Thematic Mapper (ETM+)‐derived reference imagery can be more accurately produced using such spectrally informed methods. The efficacy of several spectral index methods to discriminate between burned and unburned surfaces over a series of spatial scales (ground, IKONOS, Landsat ETM+ and data from the MOderate Resolution Imaging Spectrometer, MODIS) were evaluated. The optimal Landsat ETM+ reference image of burned area was achieved using a charcoal fraction map derived by linear spectral unmixing (k = 1.00, a = 99.5%), where pixels were defined as burnt if the charcoal fraction per pixel exceeded 50%. Comparison of coincident Landsat ETM+ and IKONOS burned area maps of a neighbouring region in Mongu (Zambia) indicated that the charcoal fraction map method overestimated the area burned by 1.6%. This method was, however, unstable, with the optimal fixed threshold occurring at >65% at the MODIS scale, presumably because of the decrease in signal‐to‐noise ratio as compared to the Landsat scale. At the MODIS scale the Mid‐Infrared Bispectral Index (MIRBI) using a fixed threshold of >1.75 was determined to be the optimal regional burned area mapping index (slope = 0.99, r ² = 0.95, SE = 61.40, y = Landsat burned area, x = MODIS burned area). Application of MIRBI to the entire MODIS temporal series measured the burned area as 10 267 km² during the 2001 fire season. The char fraction map and the MIRBI methodologies, which both produced reasonable burned area maps within southern African savannah environments, should also be evaluated in woodland and forested environments.     

Changes in snow and glacier cover in an arid watershed of the western Kunlun Mountains using multisource remote-sensing data

Snow and glaciers in the mountain watersheds of the Tarim River basin in western China provide the primary water resources to cover the needs of downstream oases. Remote sensing provides a practical approach to monitoring the change in snow and glacier cover in those mountain watersheds. This study investigated the change in snow and glacier cover in one such mountain watershed using multisource remote-sensing data, including the Moderate Resolution Imaging Spectroradiometer (MODIS), Landsat (Multispectral Scanner (MSS), Thematic Mapper (TM), and Enhanced Thematic Mapper Plus (ETM+)), Corona, and Google Earth^TM imagery. With 10 years’ daily MODIS snow-cover data from 2002 to 2012, we used two de-cloud methods before calculating daily snow-cover percentage (SCP), annual snow-cover frequency (SCF), and annual minimum snow-cover percentage (AMSCP) for the watershed. Mann–Kendall analysis showed no significant trend in any of those snow-cover characterizations. With a total of 22 Landsat images from 1967 to 2011, we used band ratio and supervised classification methods for snow classification for Landsat TM/ETM+ images and MSS images, respectively. The Landsat snow-cover data were divided into two periods (1976–2002 and 2004–2011). Statistical tests indicated no significant difference in either the variance or mean of SCPs between the two periods. Three glaciers were identified from Landsat images of 1998 and 2011, and their total area increased by 12.6%. In addition, three rock glaciers were also identified on both the Corona image of 1968 and the Google high-resolution image of 2007, and their area increased by 2.5%. Overall, based on multisource remote-sensing data sets, our study found no evidence of significant changes in snow and glacier cover in the watershed.     

Mapping land-use and land-cover change along Bolivia's Corredor Bioceánico with CBERS and the Landsat series: 1975–2008

This study uses a combination of Landsat series data (Multispectral Scanner or MSS, Thematic Mapper or TM and Enhanced Thematic Mapper or ETM+) to map land-use and land-cover change (LULCC) from 1975 to 2001. It extends the land change record to 2008 using Chinese–Brazil Earth Resources Satellite (CBERS)-2 and CBERS-2B data on a multi-scene level. It also establishes a methodology to correct for systematic distortion inherent in CBERS imagery without the loss of information present in Landsat 7 ETM+ imagery post-2003. Image analysis focuses on a 63 000 km² strip of land along a main highway and railroad in southeastern Bolivia named the Corredor Bioceánico. This strip of land is one of the most important agriculturally driven deforestation hotspots in Latin America. It is also located in one of the most poorly understood forest biomes in the world in terms of LULCC – Southern Hemisphere seasonally dry tropical forests – which have very high conservation values globally. Over the 33-year study period, approximately 12 000 km² of forest was lost among the three sub-regions – which is an area nearly the size of Connecticut. Evidence suggests that agriculture-driven deforestation is pushing into sensitive areas threatening globally important ecosystems such as those in the Chaco, Chiquitano and Pantanal as well as noteworthy protected areas. The results also show that imagery of CBERS-2 and CBERS-2B can help to fill the imagery gap created by Landsat ETM's Scan Line Corrector (SLC) failure in 2003. They can help to extend the land change record forward in time.     

Change analysis of land surface temperature based on robust statistics in the estuarine area of Pearl River (China) from 1990 to 2000 by Landsat TM/ETM+ data

The estuarine area of Pearl River that has taken great changes in land cover since 1990 is a typical area for studying the change of land surface temperature (LST). The LST of the years 1990 and 2000 in this area was estimated from the data of Landsat TM/ETM+ band 6, respectively, and three scales, corresponding to high, normal and low temperature ranges, were divided by a robust statistical method. The results show that the area of high temperature range in 2000 has an increase of 250 km² compared with the year 1990. The urban‐used land and the bare land are the main land cover types constituting the high temperature range area.     

Lake-area mapping in the Tibetan Plateau: an evaluation of data and methods

Lake area derived from remote-sensing data is a primary data source, because changes in lake number and area are sensitive indicators of climate change. These indicators are especially useful when the climate change is not convoluted with a signal from direct anthropogenic activities. The data used for lake-area mapping is important, to avoid introducing unnecessary uncertainty into long-term trends of lake-area estimates. The methods for identifying waterbodies from satellite data are closely linked to the quality and efficiency of surface-water differentiation. However, few studies have comprehensively considered the factors affecting the selection of data and methods for mapping lake area in the Tibetan Plateau (TP), nor of evaluating their consequences. This study tests the dominant data sets (Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS) data) and the methods for automated waterbody mapping on 14 large lakes (>500 km²) distributed across different climate zones of the TP. Seasonal changes in lake area and data availability from Landsat imagery are evaluated. Data obtained in October is optimal because in this month the lake area is relatively stable. The data window can be extended to September and November if insufficient data is available in October. Grouping data into three-year bins decreases the effects of year-to-year seasonal variability and provides a long-term trend that is suitable for time series analysis. The Landsat data (Multispectral Scanner, MSS; Thematic Mapper, TM; Enhanced Thematic Mapper Plus, ETM+; and Operational Land Imager, OLI) and MODIS data (MOD09A1) showed good performance for lake-area mapping. The Otsu method is used to determine the optimal threshold for distinguishing water from non-water features. Several water extraction indices, namely NDWI_McFeeters, NDWI_Xu, and AWEI_non-shadow, yielded high overall classification accuracy (92%), kappa coefficient (0.83), and user’s accuracy (~90%) for lake-water classification using Landsat data. The MODIS data using NDWI_McFeeters and NDWI_Xu showed consistent lake area (r² = 0.99) compared with Landsat data on the corresponding date with root mean square error (RMSE) values of 86.87 and 103.33 km² and mean absolute error (MAE) values of 25.7 and 29.04 km², respectively. The MODIS data is suitable for great lake mapping, which is the case for the large lakes in the TP. Although automated water extraction indices exhibited high accuracy in separating water from non-water, visual examination and manual editing are still necessary. Combined with recent Chinese high-resolution satellites, these remotely sensed imageries will provide a wealth of data for studies of lake dynamics and long-term lake evolution in the TP.     

Changes in the glaciers of Chandra–Bhaga basin,Himachal Himalaya,India, between 1980 and 2010 measured using remote sensing

This study reports the glacier changes of Chandra–Bhaga basin, northwest Himalaya, India, from 1980 to 2010. Satellite remote-sensing data from the Landsat Multispectral Scanner (MSS) and Thematic Mapper (TM), the Linear Imaging Self Scanning Sensor (LISS) and Advanced Wide Field Sensor (AWiFS) of the Indian Remote Sensing (IRS) series, and the Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) were used to study the changes in glacier parameters such as glacier area, length, snout elevation, and the impact of glacier topographical parameters (glacier slope, aspect, and altitude range) on the glacier changes. It was found that the total glaciated area had shrunk to 368.2 km² in 2010 from 377.6 km² in 1980, a loss of 2.5%. The average position of glacier terminuses retreated by 465.5 ± 169.1 m from 1980 to 2010 with an average rate of 15.5 ± 5.6 m year^?1. The decadal scale analysis showed that the average rate of retreat had increased the most in the recent decade. A moraine-dammed lake located in the study region was found to have expanded in area from (0.65 ± 0.01) km² in 1980 to (1.26 ± 0.03) km² in 2010. Glaciers with steep slope and less altitude range have lost more area than the glaciers having gentle slope and greater altitude range.     

Surveying sandy deserts and desertified lands in north‐western China by remote sensing

Sandy deserts and desertified lands (SDDL) cover most of north‐western China, and desertification is a severe environmental issue in this region. In this study, we classify SDDL into mobile, semi‐mobile, semi‐anchored, and anchored classes, and divide the corresponding severity of desertification into four levels: extremely severe, severe, moderate, and slight. Using SDDL databases derived from Landsat TM and ETM images in 1986 and 2000, we discuss the evolution of SDDL and the processes and causes of desertification. The desertified area increased by 16 101 km² between 1986 and 2000, and most of the desertified land was classified as severely or extremely severely desertified. However, 3437 km² of SDDL were also rehabilitated during this period. The area of SDDL increased in the Qinghai and Xinjiang regions, and in the western part of Inner Mongolia, but decreased in the Shaanxi, Gansu, and Ningxia regions. The causes of desertification differed among regions, and different measures to control desertification are proposed for the different regions.     

Estimation of chlorophyll‐a concentration and trophic states for inland lakes in Northeast China from Landsat TM data and field spectral measurements

Landsat TM data and field spectral measurements were used to evaluate chlorophyll‐a (Chl‐a) concentration levels and trophic states for three inland lakes in Northeast China. Chl‐a levels were estimated applying regression analysis in the study. The results obtained from the field reflectance spectra indicate that the ratio between the reflectance peak at 700 nm and the reflectance minimum at 670 nm provides a relatively stable correlation with Chl‐a concentration. Their determination of coefficients R ² is 0.69 for three lakes in the area. From Landsat TM data, the results show that the most successful Chl‐a was estimated from TM3/TM2 with R ² = 0.63 for the two lakes on 26 July 2004, from TM4/TM3 with R ² = 0.89 for the two lakes on 14 October 2004, and from the average of TM2, TM3 and TM4 with R ² = 0.72 for the three lakes tested on 13 July 2005. These results are applicable to estimate Chl‐a from satellite‐based observations in the area. We also evaluate the trophic states of the three lakes in the region by employing Shu's modified trophic state index (TSI_M) for the Chinese lakes' eutrophication assessment. Our study presents the TSI_M from different TM data with R ² more than 0.73. The study shows that satellite observations are effectively applied to estimate Chl‐a levels and trophic states for inland lakes in the area.     

Long‐term monitoring of land cover changes based on Landsat imagery to improve hydrological modelling in West Africa

The spatial and temporal variability of land cover changes is a fundamental parameter to integrate when modelling water resources in order to reproduce the relations between rainfall and surface flow more precisely. This is particularly important in West Africa, where the land cover has been changing for more than 40 years under the combined impact of climatic effects and human activities. In this study, we evaluated the potential of Landsat imagery to monitor the vegetation cover in the upper Niger watershed (120 000 km²) using archive images from MSS, TM and ETM+ sensors covering three periods of time around 1975, 1985, and 2000. Because of the heterogeneity of the acquisition dates, the spatial and spectral resolution of the images, and the scale of analysis, we chose a simple system of classification. Pretreatments were applied to reduce variations between the images. Vegetation indices (NDVI) were then calculated and subsequently thresholded using the same land‐cover classification system. The thresholds were then optimized by automated recursive calculations of confusion matrices and control parcels. Our results revealed that although the accuracy was not perfect, it was nevertheless possible to estimate changes using an unconventional spatio‐temporal scale. The resulting changes were characterized by a moderate trend to deforestation with a corresponding increase in bare soils, soils with sparse vegetation, and shrublands. The spatial layers produced were then combined with a soil map to incorporate changes in surface conditions in the hydrological modelling of the Niger River.     

Change detection technique based on the radiometric rotation controlled by no‐change axis,applied on a semi‐arid landscape

We demonstrate a new digital technique for change detection that integrates field survey data in a methodology. The field survey was specifically designed to control the operation of this new technique for change detection. The Radiometric Rotation Controlled by No‐change Axis (RCNA) technique was applied with spectral bands in the visible spectral range without atmospheric correction. The 190.37 km² study area is located in the centre of the semi‐arid region of northeast Brazil. This region is occupied by a natural shrub vegetation and low arboreal type steppe. This is known locally as ‘caatinga’, and in this area the predominant use is grazing by cattle and goat. Two TM/Landsat images in digital format of the dry seasons of 1984 and 1996 were used together with two ETM+ /Landsat images for the dry seasons of 1999 and 2001. Three change maps were generated, all with accuracy index values of around 0.70, which is acceptable for areas with complex patterns like semi‐arid lands. In conclusion, there is no need to apply atmospheric correction methods to a pair of multi‐temporal images to obtain an acceptable accuracy of change detection in semi‐arid regions. Although this technique can combine images from different sensors, this requires that both have similar spectral ranges.     

Measurement of the left‐lateral displacement of Ms 8.1 Kunlun earthquake on 14 November 2001 using Landsat‐7 ETM+ imagery

An imageodesy study has been carried out, using pre‐ and post‐event Landsat‐7 Enhanced Thematic Mapper Plus (ETM+) images, to reveal regional co‐seismic displacement caused by the Ms 8.1 Kunlun earthquake in November 2001. The two Landsat scenes, Kusai Lake and Buka Daban, cover an area of some 57 600 km² (320 km W–E and about 180 km N–S), which includes most of the fault rupture zone. The co‐seismic displacement measured in the Kusai Lake scene shows that the average left‐lateral shift along the Kunlun fault is 4.8 m (ranging from 1.5 to 8.1 m) and the maximum shift appears west of the Kusai Lake. The splayed nature of the fault to the west of Buka Daban, where the fault splits into three branches, causes the displacement pattern to become complicated. Here the average left‐lateral shift, between the south side of the southern branch and the north side of the northern branch, is 4.6 m (ranging from 1.0 to 8.2 m). Our results also illustrate that the south side of the fault is the ‘active’ block, moving significantly in an east–south‐easterly direction, relative to the largely ‘stable’ northern block.     

Analysis of land use/cover changes and urban expansion of Nairobi city using remote sensing and GIS

We used three Landsat images together with socio‐economic data in a post‐classification analysis to map the spatial dynamics of land use/cover changes and identify the urbanization process in Nairobi city. Land use/cover statistics, extracted from Landsat Multi‐spectral Scanner (MSS), Thematic Mapper (TM) and Enhanced Thematic Mapper plus (ETM+) images for 1976, 1988 and 2000 respectively, revealed that the built‐up area has expanded by about 47?km². The road network has influenced the spatial patterns and structure of urban development, so that the expansion of the built‐up areas has assumed an accretive as well as linear growth along the major roads. The urban expansion has been accompanied by loss of forests and urban sprawl. Integration of demographic and socio‐economic data with land use/cover change revealed that economic growth and proximity to transportation routes have been the major factors promoting urban expansion. Topography, geology and soils were also analysed as possible factors influencing expansion. The integration of remote sensing and Geographical Information System (GIS) was found to be effective in monitoring land use/cover changes and providing valuable information necessary for planning and research. A better understanding of the spatial and temporal dynamics of the city's growth, provided by this study, forms a basis for better planning and effective spatial organization of urban activities for future development of Nairobi city.     

Defining the Landsat TM/ETM+ and CHRIS/PROBA spectral regions in which turbidity can be retrieved in inland waterbodies using field spectroscopy

This article describes the results obtained by an existing campaign in which in situ spectroradiometric measurements using a GER1500 field spectroradiometer, Secchi disk depth, and turbidity measurements (using a portable turbidity meter) were acquired at Asprokremmos Reservoir in Paphos District, Cyprus. Field spectroradiometric and water quality data span 18 sampling campaigns during the period May 2010–October 2010. By applying several regression analyses between ‘In-Band’ mean reflectance values against turbidity values for all spectral bands corresponding to Landsat TM/ETM+ (Bands 1 to 4) and CHRIS/PROBA (Bands A1 to A62), the highest correlation was found for Landsat TM/ETM+ Band 3 (R² = 0.85) and for CHRIS/PROBA Bands A30 to A32 (R² = 0.90).     

Mapping burned areas from Landsat TM/ETM+ data with a two-phase algorithm: Balancing omission and commission errors

Maps of burned area have been obtained from an automatic algorithm applied to a multitemporal series of Landsat TM/ETM+ images in two Mediterranean sites. The proposed algorithm is based on two phases: the first one intends to detect the more severely burned areas and minimize commission errors. The second phase improves burned patches delimitation using a hybrid contextual algorithm based on logistic regression analysis, and tries to minimize omission errors. The algorithm was calibrated using six study sites and it was validated for the whole territory of Portugal (89,000 km²) and for Southern California (70,000 km²). In the validation exercise, 65 TM/ETM+ scenes for Portugal and 35 for California were used, all from the 2003 fire season. A good agreement with the official burned area perimeters was shown, with kappa values close to 0.85 and low omission and commission errors (< 16.5%). The proposed algorithm could be operationally used for historical mapping of burned areas from Landsat images, as well as from future medium resolution sensors, providing they acquire images in two bands of the Short Wave Infrared (1.5-2.2 μm).     

Technical Note: Mapping and forecasting of North Indian winter fog: an application of spatial technologies

In India, the Indo‐Gangetic plain (part of Northern India) is invariably affected by dense fog in the winter months every year due to typical meteorological, environmental and prevailing terrain conditions. Pollution also plays an important role in the formation of fog (smoke+fog = smog) in India. Using National Oceanic and Space Administration‐advanced very high resolution radiometer data the fog‐affected regions in Northern India were delineated and the spatial extent of fog for the winter months of the years 2002–03, 2003–04 and 2004–05 (December–February) were studied and mapped. Forecast for future fog based on the analysis of satellite and meteorological (air temperature, relative humidity and wind speed) data was also done. The fog‐affected areas were classified into maximum‐fog‐affected area, moderately fog‐affected area and least fog‐affected area. It has been found that in the winter months of the years 2002–03, 2003–04 and 2004–05, the fog‐affected area in Northern India was about 867 000 km², 625 000 km² and 706 800 km² respectively. The maximum fog‐affected area was found to be 606 400 km², the moderately fog‐affected area was found to be 230 400 km² and the least fog‐affected area was found to be 404 500 km². Further, based on meteorological parameters, such as temperature, humidity and wind speed along with elevation data was used to derive an approach for future fog prediction in this region.     

An operational radiometric calibration procedure for the Landsat sensors based on pseudo-invariant target sites

The Statewide Landcover and Trees Study (SLATS) use both Landsat-7 ETM+ and Landsat-5 TM imagery to monitor short-term woody vegetation changes throughout Queensland, Australia. In order to analyse more subtle long-term vegetation change, time-based trends resulting from artefacts introduced by the sensor system must be removed. In this study, a reflectance-based vicarious calibration approach using high-reflectance, pseudo-invariant targets in western Queensland was developed. This calibration procedure was used to test the existing calibration models for ETM+ and TM, and develop a consistent operational calibration procedure which provides calibration information for the MSS sensors. Ground based data, sensor spectral response functions and atmospheric variables were used as input to MODTRAN radiative transfer code to estimate top-of-atmosphere radiance. The estimated gains for Landsat-7 ETM+ (1999-2003), -5 TM (1987-2004), -5 MSS (1984-1993) and -2 MSS (1979-1982) are presented. Results confirm the stability and accuracy of the ETM+ calibration, and the suitability of this data as a radiometric standard for cross-calibration with TM. Vicarious data support the use of the existing TM calibration model for the red and two shortwave-infrared bands. However, alternative models for blue, green and near-infrared bands are presented. The models proposed differ most noticeably at dates prior to 1995, with differences in estimated gains of up to 9.7%, 10.8% and 6.9% for the blue, green and near-infrared bands respectively. Vicarious gains for Landsat-2 MSS and Landsat-5 MSS are presented and are compared with those applied by the on-board calibration system. Updated calibration coefficients to scale MSS data to the SLATS vicarious measurements are given. The removal of time based calibration trends in the SLATS data archive will enable the measurement of vegetation changes over the 26 year period covered by Landsat -2, -5 and -7.     

Characterizing changes in grassland desertification based on Landsat images of the Ongniud and Naiman Banners,Inner Mongolia

Desertification is treated as an important and critical environmental hazard. In the face of increasingly serious grassland desertification, China has made great efforts to combat desertification and a series of key national ecological projects has been launched. However, accurate, timely, and effective monitoring and assessment are required to determine whether these projects work well. The Horqin sandy land represents the largest area of sandy land in China. In this article, the Naiman and Ongniud Banners were studied because these contain the main desertified grassland in Horqin. Next, a desertification classification and grading system was designed for the Horqin sandy land after conducting fieldwork. Based on spectral mixture analysis (SMA) and decision-tree methods, we interpreted Landsat Thematic Mapper/Enhanced Thematic Mapper Plus/Operational Land Imager (TM/ETM+/OLI) images of the study area from four years: 1985, 1992, 2001, and 2013. Overall, the following results were obtained. The total area of desertified grassland in the Naiman and Ongniud Banners increased from 5979 km² in 1985 to 9195 km² in 1992 (an increase of 53.79%) and then decreased to 7828 km² in 2001 and to 6023 km² in 2013. The changes in the areas of desertified grassland with various degrees of desertification displayed the same trends as that of the total desertified grassland area. The severely desertified grassland expanded from 1872 km² in 1985 to 3723 km² in 1992 before shrinking to 2189 km² in 2013. The annual rates of expansion of desertified grassland during the periods 1985–1992, 1992–2001, and 2001–2013 were 7.68%, ?1.65%, and ?1.92%, respectively, and the corresponding expansion rates of severely desertified grassland were 14.12%, ?3.48%, and ?1.19%, respectively. Both the desertified grasslands and the areas with various degrees of desertification displayed significant expansion during the period 1985–1992. Since 1992, this trend has reduced. During the study period, the changes in temperature and precipitation did not benefit the reversal of grassland sandy desertification. Furthermore, the growing population and expansion of livestock production and farming inhibited such reversal. However, the results presented in this article suggest that a reversal in grassland sandy desertification has been occurring since 1992. The results indicate that ecological engineering measures have helped reverse desertification and have promoted the restoration of grassland vegetation.