Next Article in Journal
Estimation of Vegetation Leaf-Area-Index Dynamics from Multiple Satellite Products through Deep-Learning Method
Previous Article in Journal
Increasing Negative Impacts of Climatic Change and Anthropogenic Activities on Vegetation Variation on the Qinghai–Tibet Plateau during 1982–2019
Previous Article in Special Issue
Time Series Analysis of Landsat Data for Investigating the Relationship between Land Surface Temperature and Forest Changes in Paphos Forest, Cyprus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrating Remote Sensing and Spatiotemporal Analysis to Characterize Artificial Vegetation Restoration Suitability in Desert Areas: A Case Study of Mu Us Sandy Land

1
School of Geography and Environment, Jiangxi Normal University, Nanchang 330022, China
2
Graduate School of Global Food Resources, Hokkaido University, Sapporo 060-0809, Japan
3
Graduate School of Environmental Science, Hokkaido University, Sapporo 060-0810, Japan
4
National Engineering Research Center for Geographic Information System, China University of Geosciences, Wuhan 430074, China
5
College of Surveying and Geo-Informatics, Tongji University, Shanghai 200092, China
6
Frontiers Science Center for Intelligent Autonomous Systems, Tongji University, Shanghai 200092, China
7
College of Computer and Data Science, Fuzhou University, Fuzhou 350108, China
8
Center for Geographic Analysis, Harvard University, Cambridge, MA 02138, USA
9
Department of Geomatics, Istanbul Technical University, Istanbul 36626, Turkey
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(19), 4736; https://0-doi-org.brum.beds.ac.uk/10.3390/rs14194736
Submission received: 27 July 2022 / Revised: 25 August 2022 / Accepted: 6 September 2022 / Published: 22 September 2022
(This article belongs to the Special Issue Big Earth Observation Data Analysis for Environment Monitoring)

Abstract

:
One of the major barriers to hindering the sustainable development of the terrestrial environment is the desertification process, and revegetation is one of the most significant duties in anti-desertification. Desertification deteriorates land ecosystems through species decline, and remote sensing is becoming the most effective way to monitor desertification. Mu Us Sandy Land is the fifth largest desert and the representative area under manmade vegetation restorations in China. Therefore, it is essential to understand the spatiotemporal characteristics of artificial desert transformation for seeking the optimal revegetation location for future restoration planning. However, there are no previous studies focusing on exploring regular patterns between the spatial distribution of vegetation restoration and human-related geographical features. In this study, we use Landsat satellite data from 1986 to 2020 to achieve annual monitoring of vegetation change by a threshold segmentation method, and then use spatiotemporal analysis with Open Street Map (OSM) data to explore the spatiotemporal distribution pattern between vegetation occurrence and human-related features. We construct an artificial vegetation restoration suitability index (AVRSI) by considering human-related features and topographical factors, and we assess artificial suitability for vegetation restoration by mapping methods based on that index and the vegetation distribution pattern. The AVRSI can be commonly used for evaluating restoration suitability in Sandy areas and it is tested acceptable in Mu Us Sandy Land. Our results show during this period, the segmentation threshold and vegetation area of Mu Us Sandy Land increased at rates of 0.005/year and 264.11 km2/year, respectively. Typically, we found the artificial restoration vegetation suitability in Mu Us area spatially declines from southeast to northwest, but eventually increases in the most northwest region. This study reveals the revegetation process in Mu Us Sandy Land by figuring out its spatiotemporal vegetation change with human-related features and maps the artificial revegetation suitability.

1. Introduction

Drylands occupy 41% of the land surface on earth and their existence influences the life of 38% of the world’s people and imposes a great threat to global land sustainability [1]. Ecosystems in worldwide arid/semiarid areas appear to suffer from various levels of degradation, commonly described as “desertification” [2]. Desertification causes deteriorated socioeconomic and environmental issues, such as severe shortage of water and food and loss of biodiversity [3]. In Sustainable Development Goals (SDGs), deforestation and desertification are viewed as two of the major barriers to achieving sustainable development, especially under concurrent backgrounds of global warming and population surging in developing regions [4]. Combating desertification is a worldwide duty to conserve natural resources and sustainable development, while this may be challenged by fragile ecological environments and excessive human activities [5].
Figuring out the evolution of main deserts in China is significant to understanding land degradation problems in the central-east area. China is one of the countries that are seriously suffering from desertification. Based on a survey from 1994 to 1996, the total desert area in China was 2.62 million km2, accounting for about 27.3% of China’s land area [4,6]. One spatial characteristic of desertification in China is that it mainly dominates in an agro-pastoral transitional zone, where local people make overuse inland water and grass resources. Worse, unreasonable patterns of irrigation and grazing exacerbate soil salinization, which indirectly exacerbates the irreversibility of desertification [6]. In the past few decades, China has set more and more environmental protection programs to save northwest lands from desertification, including the Nature Reserve Development Program, Forest Eco-Compensation Program, Natural Forest Conservation Program, etc. [5].
Mu Us Sandy Land is the fifth biggest desert in China and has been reported as the first desert being eliminated by mankind. However, in the last century, related research indicates the desert area in Mu Us desert increased from 12,900 km2 in the 1950s to 41,110 km2 in the middle of the 1970s, which was probably caused by the surging population, over-grazing, over-reclaiming, and over-cutting [7]. This increasing trend of desertification has been reversed since the end of the 1970s, mainly due to the implementation and protection policies launched by both the central government and local government. Residents are encouraged in green-reserving measures to plant trees and grass (such as Tree Planting Campaign), build artificial banks to block moving sands, and applied actions in hydrologic solutions [8,9]. Based on previous desert dynamic monitoring, the desert land area of the Mu Us Sandy Land decreased from 32,568 km2 in 1986 (accounting for 67.5% of the total monitoring area) to 30,650 km2 in 1993 (accounting for 63.5% of the total monitoring area), which is mainly contributed by the conduction of the “Three North” Shelterbelt Project and the positive efforts of local officials and people to prevent desertification [10]. More compelling studies also revealed that the net area of desert land averagely decreased and the average vegetation coverage on desert land increased in recent years [9]. On the other hand, the reduction of sandy area and rising vegetation area in Mu Us Sandy Land showed that human construction activities played a positive role in the desert transformation and ecological restoration, and this greening trend is supposed to be continued in the future. Therefore, to better adapt future greening policies in Mu Us area, an optimal plan for artificial restoration development on a spatial scale is needed.
Desertification always occurs at a large global or regional scale [11]. Therefore, to save time and labor costs, remote sensing technologies are commonly applied in monitoring and assessing both desertification and anti-desertification processes. For example, in Africa, researchers acquire remote sensing archives to extract continental-scale indicators for degradation susceptibility in the Saharan desert, which effectively reflects the spatiotemporal trend of local desertification [12]. The vegetation index is a crucial value for identifying and evaluating the desertification process, so previous studies are conducted based on a comparison of the Normalized Difference Vegetation Index (NDVI) value and other related impact factors, such as precipitation [13] and Topographic Wetness Index (TWI) [14]. Applications for detecting land cover change develop with the advanced earth observation technologies. Preceding data were interpreted artificially by pictures taken from aircraft and complementary field-based observation [15,16,17]. Nowadays, aerospace databases showing land cover change are widely used for research purposes on land degradation and desertification [4,9]. Public optical satellite data with low resolution, such as Climate Change Initiative Land Cover (CCI-LC) and Moderate Resolution Imaging Spectroradiometer (MODIS), have been commonly used in monitoring the dynamics of the land surface [18,19,20,21]. Although these data products offer quality information in data-limited areas, the medium accuracy still cannot adapt to a smaller regional extent. Moreover, its insufficient time span brings a barrier to long-term earth surface monitoring.
In the past, there have mainly been two types of research on land surface monitoring by remote sensing in the Mu Us Sandy Land: (1) to directly classify the land cover classes, and then use these multi-temporal classification results for comparative analysis [4]; and (2) to compare spectral indices such as NDVI and Desertification Vulnerability Index (DVI) to reflect the gradual surface change. Both methods can effectively reflect the land cover status of the Mu Us Sandy Land to analyze and evaluate the local vegetation restoration [9]. However, vegetation and desertification indexes can only be used to reflect changes in vegetation growth or desertification level but its “non-vegetation” definition makes it difficult to use for spatiotemporal analysis or aggregation. In addition, previous studies have shown that the human and natural environmental system of the Mu Us Sandy Land has been undergoing frequent changes in recent years, which complicates the situation of restoring vegetation [22], but most existing remote sensing work on its land cover change can only cover several periods in past 30 to 40 years. These studies can generally reflect its overall changes in a long-term specific span, but for the short-term, such as on an annual or monthly basis, it can only provide a limited effect. Therefore, there is still a need to monitor the land surface on a higher time resolution for further understanding of the dynamic vegetation area in Mu Us Sandy Land.
In recent years, researchers have explored the causes of vegetation restoration in the Mu Us Sandy Land at multiple levels. Liu has used panel data to analyze and evaluate the relative role of climate, socio-economic conditions, and policy effects at a province level, showing the crucial positive effect of policy implementations on vegetation and ecological restoration [23,24]. Xiu and his group quantitatively analyzed the effect of the ecological engineering projects on the vegetation restoration of Mu Us and emphasized the importance of human activities to the growth of regional vegetation under the background of climate change [25]. Sun used linear regression methods to compare the vegetation recovery effects of climate factors and afforestation and found afforestation activities are more effective [26]. Most views believe that human interference imposed a significant effect on vegetation restoration, but there is no spatially explicit evidence to demonstrate the sensitivity of vegetation restoration to human activities, because most studies have made the “restoration process” a confirmed assumption and related it with vegetation changes [6,24,27]. Therefore, it is particularly important to explore the relationship between spatiotemporal vegetation changes and explicit geographical features of human activities.
A series of previous studies comprehensively assessed the vegetation suitability in Mu Us Sandy Land and worldwide dry areas by modeling and remote sensing methods. For example, Wang and his team integrated hydrological and ecological models to evaluate the sustainable vegetation restoration Loess Plateau of China, based on natural factors such as rainfall and soil [28]. McVicar and his group simulated vegetation suitability in Loess Plateau by considering land uses, topographical characteristics, and water use [29]. Peng and his research team exploited hidden vegetation suitability at the vegetation scale in terrestrial ecosystem models to apply to revegetation programs under diverse climate change scenarios [30]. In Africa, remote sensing-based vegetation indices were proven to be effective indicators to assess the reactions of terrestrial ecoregion drought [31]. However, to the best of our knowledge, there are no studies focusing on mapping vegetation suitability based on different revegetation convenience levels by considering the spatiotemporal distribution of human footprint features. The revegetation convenience refers to the different difficulties of construction in revegetation activities, such as different slopes causing different difficulties in restoration measurements. To bridge this knowledge gap, this study combines methods of satellite remote sensing and spatial and temporal analysis to explore the impact of human-related features on vegetation change.
There are three objectives on which our research is focused: (1) to achieve annual vegetation and non-vegetation classification in Mu Us Sandy Land; (2) to explore the relationship between spatial vegetation distribution patterns and human-related footprint features; (3) to construct a method to assess the artificial vegetation restoration suitability in desert areas and apply it in Mu Us case. Therefore, we designed our study contents based on the above three purposes, firstly, we exploit an Otsu-based threshold segmentation method to automatically extract annual vegetation pixels based on an improved Enhanced Normalized Difference Vegetation Index (ENDVI), for providing a data base of image classification and following spatiotemporal analysis. Secondly, we selected the geographic information features related to human activities and used the spatiotemporal analysis methods to perform vegetation distribution in these features to explore the regular pattern of vegetation change in Mu Us Sandy Land. Thirdly, we construct a novel AVRSI and assess the vegetation suitability of manmade restoration in the Mu Us Sandy Land by analyzing the influence of the dominant characteristics of human activities on the distribution of vegetation. Our research provides a scientific reference for making more sustainable vegetation restoration policies or projects in the Mu Us Sandy Land and other deserts.

2. Materials and Methods

2.1. Study Area

Mu Us Sandy Land, lying on the Ordos Plateau in China, administratively belonged to the Inner Mongolia Autonomous Region, Shaanxi Province, and the Ningxia Hui Autonomous Region. Its district is mainly located between 107°21′E and 111°30′E, and 37°27′N and 39°22′N, with an elevation from 0 m to 1614 m a.s.l (Figure 1) and an area of about 42 thousand km2 [6]. Continental semi-arid climate is its representative climate characteristic. The annual temperature ranges from 6 °C to 8.5 °C, with monthly mean temperatures of 22 °C in July and −11 °C in January. Annual precipitation varies from 440 mm in the southeast to 250 mm in the northwest, which is up to 80%, particularly from June to August.
The main type of land surface with vegetation is sandy grassland and the main plantation for cultivation are corn, buckwheat, potato, etc. [4]. There are mainly four landscape types in Mu Us Sandy Land: sandstone hills, meadow steppe, and active and fixed sand dunes [32]. Several rivers lie across the southeast of the region and gather into the Yellow River. Lakes are mainly distributed in the internal land of the area, but most of them are sodic lakes with chloride, and only a few are freshwater lakes [9].

2.2. Data Preparations

The process of vegetation classification is conducted through the cloud-based platform Google Earth Engine (GEE) and Pixel Information Expert Engine (PIE) (https://engine.piesat.cn/ (accessed on 9 September 2021)). Our study period lasts for 35 years; therefore, the historical images for monitoring vegetation areas are derived from Landsat Thematic Mapper (TM), Landsat Enhanced Thematic Mapper Plus (ETM+), and Landsat Operational Land Imager (OLI) with 30 m spatial resolution. All Landsat Tier1 level data are created using the best available processing level for each scene. The Landsat images with rare clouds at the Tier1 level were collected based on time and space conditions. According to the processing results, there are no data available covering the study area before 1986, so the images selected in the 35 years of data from 1986 to 2020 were counted. Considering climate factors and period consistency, images synthesized every summer are used as a study database. Totally 2785 scenes are used in our image synthesis process. The satellite data sources, period of selected scenes, and the number of scenes for each year are described in Table S1 (Supplementary File).
Human information features include roads, railways, human settlements, service stations, and waterways, which are collected from OSM and updated to the 2020 version. We also selected slope as a topographical feature for mapping the artificial suitability of vegetation restoration, which is derived from the Shuttle Radar Topography Mission (SRTM) at a scale of 30 m.

2.3. Methods

This research consists of two main sections: annual vegetation detection by remote sensing and spatiotemporal analysis with human-related factors. The final aim of this study is to map artificial vegetation restoration suitability based on the two above methods. Figure 2 indicates the data input, classification algorithm, related spectral index, and specific steps of each section.

2.3.1. Annual Vegetation Detection by Remote Sensing

In this section, we first construct ENDVI by original NDVI, NDSI, and MNDWI and then use the Otsu method to classify vegetation and non-vegetation areas in Mu Us Sandy Land. Last, we made statistics on the analysis results. This analysis will provide a panel reference for past vegetation change in Mu Us Sandy Land and this result will be used for the following analysis of the relationship between vegetation change and human features in the next section.
After image pre-processing, such as radiation correction, image mosaic, and band synthesis, we calculate three spectral indices NDVI, Modified Normalized Difference Water Index (MNDWI), and Normalized Difference Soil Index (NDSI) that, respectively, characterize vegetation, water, and bare land. Then, we unified these three indexes as values from −1 to 1 and mask the area of identified water bodies. To expand the background value difference between vegetation and bare land, we subtract NDSI from the value of NDVI, constructing a new ENDVI threshold under the background of weakened bare soil. This Otsu-based segmentation method with spectral index is commonly used in earth surface detection, but a proper index for differentiating vegetation and non-vegetation is still under exploration in desert areas. Therefore, after the above processing of coupled indices, the difference of ENDVI between vegetation and non-vegetation will become larger, which makes it more adaptive in arid areas. Lastly, we normalize the ENDVI value from 0 to 1 before processing it into vegetation classification. The related spectral index equations and derivation process of the ENDVI formula are shown in the following:
E N D V I = N D V I N D S I  
N D V I = N I R R e d N I R + R e d
N D S I =   S W I R N I R     S W I R + N I R  
Among them, Red, NIR, and SWIR refer to the red band, near-infrared red band, and shortwave infrared band of satellite images, respectively.
The Otsu method is a nonparametric and unsupervised method of automatic threshold selection for picture segmentation, which is used to maximize the difference between vegetation and non-vegetation. Assuming that image pixels can be divided into two parts based on distinct gray levels: background and objects. When the variance of the gray value between the objects and the background reaches the maximum, the threshold can be the optimal segmentation threshold. Currently, the difference between the target and the background is the largest, and the segmentation is the most effective [33]. It can be briefly introduced that:
Assuming that the gray level of an image is L, which presents as a level range [0,1,2...L − 1]. The total number of pixels in the image is N, n i represents the number of pixels with gray level i, namely:
N = i = 0 L 1 n i
Pi represents the probability of the expected point with gray level i, namely:
p i = n i N
If the pixels are divided into two classes, object α1 and background α2, by the threshold t, α1 is composed of pixels whose gray value is between [0, t], and α2 is composed of pixels whose gray value is between [t + 1, L − 1]. Then the probabilities of these two classes w1 and w2 are:
w 1 = i = 0 t P i
w 2 = i = t + 1 L 1 P i = 1 w 1
The means of these two classes can be computed as:
μ 1 = i = 0 t i P i w 1 = μ n w 1
μ 2 = i = t + 1 L 1 i P i w 2 = μ N μ n 1 w 1
Among them:
μ N = i = 0 L 1 i P i
μ n = i = 0 t i P i
If σ b 2 is used to represent the between-class variance between the two classes, then it is:
σ b 2 = w 1 μ 1 μ N 2 + w 2 μ 2 μ N 2 = w 1 μ 1 2 + μ N 2 + μ N 2 w 1 + w 2 2 w 1 μ 1 + w 2 μ 2 = w 1 w 2 μ 1 μ 2 2
Let the t values be taken in sequence from the range of [0, L − 1], when σ b 2 becomes largest, the corresponding t value is the optimal threshold obtained by the Otsu algorithm [34,35].
The Otsu algorithm is to process the pixels of the entire area, and then calculate the optimal threshold, which would be ineffective if an imbalanced proportion of the object and background pixels occurs. In this study, to realize a more equal ratio of the object and background pixels, canny edge detection was performed, and buffer zones were identified based on the contour line before employing Otsu segmentation. NDVI value is a general indicator to demonstrate the greenness of vegetation visually and numerically on an observed land surface. In general, NDVI values below 0.1 signify a low level of vegetation cover, namely barren areas such as rock, sand, and snow. Moderate values from 0.2 to 0.3 correspond to shrub and grassland, and relatively high levels from 0.6 to 0.8 always indicate more bushy forests [36]. In this study, vegetation and non-vegetation area are divided by the Otsu ENDVI threshold. Therefore, to detect the distribution of vegetation as much as possible, we assume vegetation areas possess two attributes: the ENDVI value is greater than the Otsu ENDVI threshold and the original NDVI is over 0.1.

2.3.2. Spatiotemporal Analysis for Vegetation Suitability Mapping

By reviewing works of literature, we found that the exact actions of ecological restorations mainly include constructing windscreens to secure sand and planting shrubs and grasses [6], and the effective implementation of these measures is strongly related to the existence of transportation and settlement services and water resources; thus, we selected five human-related features to reflect the ecological restoration footprint features, including roads, railways, human settlements, service stations, and waterways. Euclidean distance refers to the shortest straight-line distance between pixels in the geographical study, which is commonly used for referring distance convenience between two locations [37,38]. Due to the limitations of the availability and format of analysis data, we adopted the Euclidean distance to conduct long time-series spatial statistics to explore the rules of spatial accessibility between each pixel and human-related features in the study area.
To explore the relationship between vegetation distribution and the above geographical features of human activities, we use the Euclidean distance method to analyze the average distance from each vegetation pixel to these features. First, we use the int function in ArcGIS Pro to convert human-related features in vector form into the same raster form. Then, we calculated the Euclidean distance from each grid to these human-related features. The selected human-related features for time-series analysis are only 2020-based data; therefore, if the analyzed period for the past is too long, these features might show a big difference from 2020-based data. Therefore, we chose the period of ten years only from 2011 to 2020 for analysis to avoid excessive feature changes from affecting the analysis.

2.3.3. Constructing an AVRSI

Typically, land suitability is defined in terms of specific activities or land uses [39,40,41]. Due to worldwide changing economic development and social construction, the theories, methods, and practices of land suitability assessment have made rapid progress, and the application ranges of land suitability evaluation are also constantly improving [41,42]. Based on previous studies, the main fields of land suitability evaluation application can be basically divided into four aspects: agricultural land, forestry and animal husbandry land, urban construction land, tourism development land, and reclamation and consolidation land [43,44,45,46,47].
We compute the AVRSI and generate a suitability map based on the following three principles by referring distance analysis of spatiotemporal vegetation change: (1) accessibility and rationality of human-related feature data—we select the geographical data which can be acquired and closely related to restoration activities; (2) correlation between distances from human-related features and vegetation suitability can be scaled—we assume that vegetation area is expanding in the study area and its growth will linearly decline with the distance from human-related features; and (3) the effect weights of each human-related features on vegetation restoration can be quantified.
In this study, we focus on combining the distribution of geographical features of ecological restoration activities and the law of land development and utilization to construct an AVRSI and map the suitability of artificial vegetation restoration, to achieve labor-consuming artificial production of vegetation in desert areas. The defined expression of AVRSI is described below:
AVRSI = T i = 1 n x i α i
where n is the total number of selected human-related features; α i is the effect weights of each human-related feature; x i indicates the distance convenience for artificial restoration; T is the coefficient ratio of topographical factors on the targeted grid.

2.3.4. Mapping Artificial Vegetation Restoration Suitability

Based on the above AVRSI computing principles for mapping the suitability of anthropogenic vegetation restoration, we use the functions of GIS to spatially express the suitability of the Mu Us Sandy Land: (1) we compute the Euclidean distance of each grid for the five human-related features and use the Int function to define each grid value as the Euclidean distance value (five human-related features correspond to five sets of Euclidean distance values), and we use the natural breaks (Jenks) method to classify the distance convenience, which is arranged as very high, high, moderate, low, and very low (quantitatively defined as 9:7:5:3:1 in descending order); (2) we determine the weights of each characteristic factor through the buffer analysis of vegetation distribution, and then we define the ratio of their weights as the ratio of their rounded regression equation slopes, namely, railway: roads: waterways: human settlements: service station is 1:2:3:2:1. In addition, we define the slope grades into 5 levels based on the Manual of Detailed Geomorphological Mapping [48]: 0°–2°, 2°–4°, 5°–14°, 15°–35°, and 36°–83°, and then we define the suitability coefficient ratio of these five grades as 1:0.9:0.7:0.5:0.1 based on our understanding of the influence of slope on human construction activities; and (3) we first calculate the weighted distance from each grid to these five human-related features through the raster calculator function, and then multiply the data of the modified weighted distance by the slope coefficient ratio of the corresponding grid, as shown in Equation (13).

3. Results

3.1. ENDVI

The ENDVI values distribution and the corresponding number of pixels are shown in Figure S1 (Supplementary File) and Figure 3. From the two figures, we can find the overall ENDVI value in the Mu Us Sandy Land was relatively distributed between 0.4 and 0.6, and the number of pixels with ENDVI below 0.4 became less and the number of pixels with ENDVI over 0.6 became more. In addition, the ENDVI value points of which the count of pixels peaks every year are also increasing.
The threshold segmentation values of ENDVI for distinguishing vegetation from non-vegetation in Mu Us Sandy Land from 1986 to 2020 are shown in Table S2 (Supplementary File). During this period, the overall fluctuation of the segmentation threshold was only around 0.02 (Figure 3), which indicates a relatively small difference and supports the acceptable effect on vegetation classification based on the Otsu segmentation method. The ENDVI threshold of Mu Us Sandy Land has increased significantly (R2 = 0.48, p < 0.01) at a change rate of 0.005/year from 1986 to 2020 (Figure 4a).

3.2. Annual Vegetation Change

Figure 4b and Figure S2 (Supplementary File) show the results of the annual vegetation change segmented by the Otsu method, of which general change is similar to that of the “greenness degree”. Both results show a fluctuant trend before 2000 but an obvious uptrend after 2000. The vegetation area of Mu Us Sandy Land has obviously increased (R2 = 0.35, p < 0.01) at a rate of 264.11 km2/year from 1986–2020 (Figure 4b). Regardless of the temporal change, the spatial pattern of vegetation distribution shows a general increase trend from southeast to northwest.

3.3. Spatiotemporal Analysis between Time-Series Vegetation Dynamic and Human-Related Factors

We computed the Euclidean distance from each grid to these human-related features, as shown in Figure 5, which we could use as spatial accessibility between vegetation and human-related features for restoration suitability. Mean Euclidean distances from the vegetation grid to human-related features are calculated annually as Figure 6.
From the results of the line chart in Figure 6, we found that the average Euclidean distance from the vegetation grid to these human-related features increased over the targeted ten years, showing that the overall boundaries of vegetation growth were expanding toward areas where it was not before, and the start points of these movements are human-related features. Therefore, the spatiotemporal characteristics of human ecological restoration measures are to restore from the areas closer to the human-related features to the area far from these features, which indicates that the difficulty of vegetation restoration increases with the distance between the targeted restoration location and the human-related features.
In addition, we analyze the relationship between vegetation distribution and human-related features in 2020 by establishing a buffer zone (10 km in total) with an interval of 1 km. The result is displayed in Figure 7. We found that as the distance from the geographic elements of human activities increased the proportion of vegetation grids in the total buffer grid decreased.
We assume that the distance from any locations in the Mu Us Sandy Land to these five geographical elements of human activities is linearly negatively correlated with the suitability of vegetation restoration. This result shows that human activities played a positive role in the process of vegetation restoration in Mu Us from the perspective of distance analysis. Based on the results shown in Figure 7, in 2020, we use the ratio of the slope of each trend of these features to represent the weight ratio of each element in mapping the suitability of vegetation restoration.
As one of the dominant topographical features, the slope reflects the steepness of the land surface, which has a significant impact on human-driven restoration activities [4,49,50]. As the slope in Mu Us sandy land shows nearly unchanged over the last 40 years, we count the slope values of each vegetation grid from 1986 to 2020 and calculate their averages for each year to explore the impact on vegetation distribution at Mu Us Sandy Land. The result is shown in Figure 8.
The spatial range of natural vegetation growth should be from the area with a low slope to the area with a higher slope. However, the slope of the vegetation grid increases over years from 1986 to 2020, following the changes in ENDVI threshold and vegetation area over the targeted years, which indicated that human-related restoration promoted the vegetation expansion. From this perspective, we assume that the suitability of anthropogenic vegetation restoration decreases with the rising slope. This summary provides guidance for our subsequent vegetation restoration suitability.

3.4. Vegetation Restoration Suitability Mapping

The suitability results of anthropogenic vegetation restoration are shown in Figure 9. We can find that the overall spatial distribution of vegetation restoration suitability is decreasing from the southeast to the northwest, and after reaching the lowest point of suitability in the central–western area, it shows an upward trend in the northwest eventually. The probable causes of the restoration differences are attributed to the distribution of more dense waterways and flatter plains in the southeast, while the terrain in the central and western regions shows more large undulating slopes and is far from the water source. However, more dense road distribution and gentler slope lead to upward suitability in the most northwest region.
Integrating the annual changes in the vegetation distributions in the Mu Us Sandy Land with the regular patterns of artificial transformation found in result 3.3 that is, vegetation tends to be restored closer to waterways, roads, and railways, this suitability result is well-matched with the actual expectations, and it also shows that the vegetation suitability index designed in Section 2.3.3 is applicable to the progress of artificial vegetation restoration in the Mu Us Sandy Land over recent years.

4. Discussion

4.1. Anti-Desertification Process in Mu Us Sandy Land

Before the 21st century, human activities, such as agricultural and animal husbandry activities, and mineral resources development, strongly disturbed the natural environment in the Mu Us region, leading to the instability of the “greenness degree” and vegetation area [22,51]. Since the beginning of this century, when human activities have become ecological restoration activities, desertification showed the characteristics of reversal in this desert. Under the governmental implementations of the multipolicy, such as the “Grain for Green Project”, “The Three Northern Regions (northeastern, northwestern, and northern China) Shelter”, and “mine ecological restoration”, the vegetation area and the “greenness degree” have significantly increased in Mu Us Sandy Land [24,51,52]. This environmental benefit can be indicated from the result of Figure 4a and Figure S2 by our remote sensing work, the upward trend of the segmentation threshold and vegetation area indicates that the “greenness degree” and vegetation area in the Mu Us Sandy Land are continuously developed, which is basically consistent with the previous studies by many researchers [25,53,54].
The driving factors for the vegetation restoration in Mu Us Sandy Land mainly include natural driving factors and anthropogenic driving factors. The phenomenon of desertification existed throughout the quaternary period; thus, the impact of human activities on the landscape change before the modern era only accounted for a small effect of the entire desertification process [55]. Arid and windy extreme events deteriorate the stability of ecosystems, exacerbating soil barrenness, wind erosion, and vegetation degeneration [56], but the previous study clarified the climate factor is one of the minimally significant factors influencing the ecological restoration process in the Mu Us desert areas [24], which is the main reason our study exclude climate change factor as consideration for suitability mapping. On the other hand, along with the increasing frequency of human activities, the improvement of productivity levels, and the expansion of the scope and content of agricultural activities in arid areas, anthropogenic factors have also become crucial factors in the changes of desert vegetation [57]. These causes of vegetation degradation in arid areas include excessive grazing, deforestation, unreasonable farming, and overuse of water resources caused by the increasing population and inordinate utilization of natural resources. However, the impact of anthropogenic factors on desert vegetation can be bidirectional. We found that the changes in the “greenness degree” and vegetation area in Mu Us show a time-series fluctuating increase. The fluctuation shows that the unreasonable socioeconomic activities that still exist in the Mu Us area, mainly occurring the unreasonable agricultural and animal husbandry development, indiscriminate logging, deforestation, intensive mining activities, and possible sand dune movement, while the results show an overall increasing trend indicates that ecological restoration activities dominate the change of Mu Us vegetation [51].

4.2. Future Restoration

Combining the suitability mapping results in Figure 10 with the vegetation cover in 2020 (the most recent targeted year in this study), we can obtain the distribution of the difficulty level of implementation of anthropogenic vegetation restoration, as the results shown in Figure 10.
If the vegetation restoration measures in the Mu Us Sandy Land continue to be promoted in the future, the non-vegetation areas in the southeast and most northwest areas show better adaptability, indicating that the obstacles to vegetation restoration are relatively small and it would be more efficient to achieve significant results there. Therefore, these areas with higher suitability should be restored to vegetation first, rather than those areas with lower suitability.
Our study outcomes are aiming to save labor and economical costs during future ecological; however, due to the lack of data for specific costs for restoration processes, such as transportation of plants and irrigation, we cannot give a quantitative analysis for evaluation of its economic feasibility. Therefore, making a more comprehensive assessment of artificial restoration suitability by collecting panel data for specific restoration costs is significant in our future work.

4.3. Limitations and Future Work

This study uses the, x i , α i , and T (Equation (13)) as influencing factor coefficients to a novel vegetation suitability index for the anthropogenically transforming desert, which acts as a novel method to achieve mapping vegetation suitability in the Mu Us region based on results of remote sensing interpretation with an originally designed ENDVI and spatiotemporal analysis.
However, there exist the following limitations in this study: (1) the spatial resolution of the remote sensing data used in this study is 30 m, which indicates that there are still some details that need to be explained more accurately; (2) due to the availability of data, the selected anthropogenic geographical features have relatively limited influence on the distribution of vegetation restoration and also probably change in future, causing the impracticality of suitability mapping; (3) meanwhile, the overlapping effects of multiple factors cannot be quantitatively valued; (4) although the Mu Us Sandy Land is a representative study area due to its special agropastoral location and on-going vegetation restoration activities, it is inevitable that regional differences still exist. Typically, since a previous study has explained the climate factor accounts for limited weight in Mu Us vegetation change [24], this study is a special case for transforming desert which is less influenced by climate conditions; (5) since the restoration measures and financial costs of different vegetation species are different, the lack of subclassification of vegetation classes will lead to a certain difference between the modeled suitability and the actual suitability.
We are supposed to solve the above limitations from the following aspects: (1) to use satellite images with a higher spatial resolution for monitoring; (2) to investigate more data sources and update them for optimizing the completeness of selected human-related features and topographical factors, and if deems necessary, consider more additional other related factors; (3) to use proper weight analysis methods to rank vegetation change driven factors; (4) when applying the vegetation suitability mapping methods to tropical and subtropical regions, the climatic and geomorphological characteristics, such as surface temperature, precipitation, soil moisture content, soil texture, and wind force, are supposed to be additionally considered; and (5) to subdivide vegetation types by higher resolution remote sensing images or ground observation, and use this more detailed vegetation classification results to comprehensively make artificial restoration suitability for different vegetations.

5. Conclusions

This study uses a remote sensing method based on an automatic threshold segmentation algorithm to monitor annual vegetation change in Mu Us Sandy Land from 1986 to 2020 by an introduced novel enhanced normalized vegetation index. Based on analysis of the threshold segmentation and vegetation classification results, we have found, albeit with some fluctuations, that both the “greenness degree” and vegetation area show an obvious overall upward trend, especially after 2000. Then, we use distance analysis methods to investigate regular patterns between time-series vegetation distribution and human activity footprint features. These features include roads, railways, human settlements, service stations, waterways, and topographic factors. The results show that the vegetation restoration is expanding from the area with a lower slope near human-related features to the area with a higher slope far from these features, and this pattern effectively explains the vegetation restoration did cause by human activities from the perspective of spatiotemporal analysis. Moreover, we exploit the buffer zone function to analyze the relationship between vegetation distribution and human-related features, of which results indicate that the area portion of vegetation occurrence decreases with the distance from human-related features to the target location, and waterways, roads, and human settlements act as more influential geospatial features. Last, we construct an AVRSI for manmade desert transformation to map vegetation suitability based on distance analysis of human-related features and topographical factors. In the case of the Mu Us Sandy Land, we find artificial restoration vegetation suitability spatially decreases from southeast to northwest, but it shows a considerably increasing trend in the most northwest region. This study offers an effective means for quantitatively evaluating the resistance strength and expected effect of manmade desert transformation and provides data support and reference suggestions for relevant governments and organizations to formulate future policies and plans for artificial desert transformation.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/rs14194736/s1, Table S1. Data sources, period of selected scenes, and the number of scenes for each year from 1986 to 2020; Table S2. Annual segmentation thresholds for Otsu method; Figure S1. Relationship between ENDVI value and number of corresponding pixels; Figure S2. Annual detected vegetation in the Mu Us region from 1986 to 2020.

Author Contributions

Z.C.: Writing—original draft, Data curation, Formal analysis, Investigation, Methodology, Visualization. M.H.: Writing—review and editing, Investigation, Methodology, Conceptualization, Software, Funding acquisition, Project administration, Supervision. C.X.: Writing—review and editing. S.Q.: Writing—review and editing. W.D.: Writing—review and editing. D.Z.: Writing—review and editing. O.A.: Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Nature Science Foundation of China (NSFC) program (No. 42201438, 42001372, 41867012, 41971351), the JST SPRING (Grant Number JPMJSP2119), the fundamental research funds for the central universities (No. 2042020kf0011), Open Fund of National Engineering Research Center for Geographic Information System, China University of Geosciences, Wuhan 430074, China (Grant No. 2021KFJJ07), the China Postdoctoral Science Foundation (2021T140513, 2019M661621), the National Key R&D Program (2021YFF0704400), the Special Fund of Hubei Luojia Laboratory (No. 220100034), and the talent startup fund of Jiangxi Normal University.

Data Availability Statement

The data presented in this study are cited within the article.

Acknowledgments

We appreciate the constructive suggestions and comments from the editor and anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reynolds, J.F.; Stafford Smith, D.M.S.; Lambin, E.F.; Turner, B.L.; Mortimore, M.; Batterbury, S.P.; Downing, T.E.; Dowlatabadi, H.; Fernández, R.J.; Herrick, J.E.; et al. Global Desertification: Building a Science for Dryland Development. Science 2007, 316, 847–851. [Google Scholar] [CrossRef]
  2. Hillel, D.; Rosenzweig, C. DESERTIFICATION. In Encyclopedia of Soils in the Environment; Elsevier: Amsterdam, The Netherlands, 2005; pp. 382–389. [Google Scholar] [CrossRef]
  3. Jiang, L.; Bao, A.; Jiapaer, G.; Guo, H.; Zheng, G.; Gafforov, K.; Kurban, A.; De Maeyer, P. Monitoring land sensitivity to desertification in Central Asia: Convergence or divergence? Sci. Total Environ. 2019, 658, 669–683. [Google Scholar] [CrossRef]
  4. Li, Y.; Cao, Z.; Long, H.; Liu, Y.; Li, W. Dynamic analysis of ecological environment combined with land cover and NDVI changes and implications for sustainable urban–rural development: The case of Mu Us Sandy Land, China. J. Clean. Prod. 2017, 142, 697–715. [Google Scholar] [CrossRef]
  5. Chen, F.; Li, G.; Zhao, H.; Jin, M.; Chen, X.; Fan, Y.; Liu, X.; Wu, D.; Madsen, D. Landscape evolution of the Ulan Buh Desert in northern China during the late Quaternary. Quat. Res. 2014, 81, 476–487. [Google Scholar] [CrossRef]
  6. 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] [CrossRef]
  7. Salvati, L.; Kosmas, C.; Kairis, O.; Karavitis, C.; Acikalin, S.; Belgacem, A.; Solé-Benet, A.; Chaker, M.; Fassouli, V.; Gokceoglu, C.; et al. Assessing the effectiveness of sustainable land management policies for combating desertification: A data mining approach. J. Environ. Manag. 2016, 183, 754–762. [Google Scholar] [CrossRef]
  8. Gao, Z.; Li, C.; Yang, X.; Ci, L. Engineering and Technological Measures for Combating Desertification. In Desertification and Its Control in China; Springer: Berlin/Heidelberg, Germany, 2010; pp. 427–463. [Google Scholar] [CrossRef]
  9. 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 Sens. 2014, 6, 9316–9339. [Google Scholar] [CrossRef]
  10. Wu, W. Study on process of desertification in Mu Us sandy land for last 50 years. China J. Desert Res. 2001, 21, 164. Available online: http://www.desert.ac.cn/EN/Y2001/V21/I2/164 (accessed on 6 July 2021).
  11. Prince, S.D. Where Does Desertification Occur? Mapping Dryland Degradation at Regional to Global Scales. In The End of Desertification? Disputing Environmental Change in the Drylands; Springer: Berlin/Heidelberg, Germany, 2016; pp. 225–263. [Google Scholar] [CrossRef]
  12. Symeonakis, E.; Drake, N. Monitoring desertification and land degradation over sub-Saharan Africa. Int. J. Remote Sens. 2004, 25, 573–592. [Google Scholar] [CrossRef]
  13. Hountondji, Y.; Sokpon, N.; Ozer, P. Analysis of the vegetation trends using low resolution remote sensing data in Burkina Faso (1982–1999) for the monitoring of desertification. Int. J. Remote Sens. 2006, 27, 871–884. [Google Scholar] [CrossRef]
  14. Lamchin, M.; Lee, W.-K.; Jeon, S.W.; Lee, J.-Y.; Song, C.; Piao, D.; Lim, C.H.; Khaulenbek, A.; Navaandorj, I. Correlation between Desertification and Environmental Variables Using Remote Sensing Techniques in Hogno Khaan, Mongolia. Sustainability 2017, 9, 581. [Google Scholar] [CrossRef]
  15. Chen, Z.; Huang, M.; Zhu, D.; Altan, O. Integrating Remote Sensing and a Markov-FLUS Model to Simulate Future Land Use Changes in Hokkaido, Japan. Remote Sens. 2021, 13, 2621. [Google Scholar] [CrossRef]
  16. Huang, M.; Chen, N.; Du, W.; Wen, M.; Zhu, D.; Gong, J. An on-demand scheme driven by the knowledge of geospatial distribution for large-scale high-resolution impervious surface mapping. GIScience Remote Sens. 2021, 58, 562–586. [Google Scholar] [CrossRef]
  17. Liu, Y.Y.; Van Dijk, A.; De Jeu, R.A.M.; Canadell, J.; McCabe, M.; Evans, J.; Wang, G. Recent reversal in loss of global terrestrial biomass. Nat. Clim. Chang. 2015, 5, 470–474. [Google Scholar] [CrossRef]
  18. He, Y.; Lee, E.; Warner, T.A. A time series of annual land use and land cover maps of China from 1982 to 2013 generated using AVHRR GIMMS NDVI3g data. Remote Sens. Environ. 2017, 199, 201–217. [Google Scholar] [CrossRef]
  19. Hua, T.; Zhao, W.; Liu, Y.; Wang, S.; Yang, S. Spatial Consistency Assessments for Global Land-Cover Datasets: A Comparison among GLC2000, CCI LC, MCD12, GLOBCOVER and GLCNMO. Remote Sens. 2018, 10, 1846. [Google Scholar] [CrossRef]
  20. Sulla-Menashe, D.; Gray, J.M.; Abercrombie, S.P.; Friedl, M.A. Hierarchical mapping of annual global land cover 2001 to present: The MODIS Collection 6 Land Cover product. Remote Sens. Environ. 2019, 222, 183–194. [Google Scholar] [CrossRef]
  21. Zhang, H.K.; Roy, D.P. Using the 500 m MODIS land cover product to derive a consistent continental scale 30 m Landsat land cover classification. Remote Sens. Environ. 2017, 197, 15–34. [Google Scholar] [CrossRef]
  22. Li, Y.; Wang, J.; Liu, Y.; Long, H. Problem regions and regional problems of socioeconomic development in China: A perspective from the coordinated development of industrialization, informatization, urbanization and agricultural modernization. J. Geogr. Sci. 2014, 24, 1115–1130. [Google Scholar] [CrossRef]
  23. Liu, Q.; Zhao, Y.; Zhang, X.; Buyantuev, A.; Niu, J.; Wang, X. Spatiotemporal Patterns of Desertification Dynamics and Desertification Effects on Ecosystem Services in the Mu Us Desert in China. Sustainability 2018, 10, 589. [Google Scholar] [CrossRef]
  24. Liu, Q.; Zhang, Q.; Yan, Y.; Zhang, X.; Niu, J.; Svenning, J.-C. Ecological restoration is the dominant driver of the recent reversal of desertification in the Mu Us Desert (China). J. Clean. Prod. 2020, 268, 122241. [Google Scholar] [CrossRef]
  25. Xiu, L.; Yan, C.; Li, X.; Qian, D.; Feng, K. Monitoring the response of vegetation dynamics to ecological engineering in the Mu Us Sandy Land of China from 1982 to 2014. Environ. Monit. Assess. 2018, 190, 543. [Google Scholar] [CrossRef]
  26. Sun, Z.; Mao, Z.; Yang, L.; Liu, Z.; Han, J.; Wanag, H.; He, W. Impacts of climate change and afforestation on vegetation dynamic in the Mu Us Desert, China. Ecol. Indic. 2021, 129, 108020. [Google Scholar] [CrossRef]
  27. Huang, Y.; Wang, N.-A.; He, T.; Chen, H.; Zhao, L. Historical desertification of the Mu Us Desert, Northern China: A multidisciplinary study. Geomorphology 2009, 110, 108–117. [Google Scholar] [CrossRef]
  28. Wang, C.; Wang, S.; Fu, B.; Lü, Y.; Liu, Y.; Wu, X. Integrating vegetation suitability in sustainable revegetation for the Loess Plateau, China. Sci. Total Environ. 2020, 759, 143572. [Google Scholar] [CrossRef]
  29. McVicar, T.R.; Van Niel, T.; Li, L.; Wen, Z.; Yang, Q.; Li, R.; Jiao, F. Parsimoniously modelling perennial vegetation suitability and identifying priority areas to support China’s re-vegetation program in the Loess Plateau: Matching model complexity to data availability. For. Ecol. Manag. 2010, 259, 1277–1290. [Google Scholar] [CrossRef]
  30. Peng, S.; Yu, K.; Li, Z.; Wen, Z.; Zhang, C. Integrating potential natural vegetation and habitat suitability into revegetation programs for sustainable ecosystems under future climate change. Agric. For. Meteorol. 2019, 269–270, 270–284. [Google Scholar] [CrossRef]
  31. Robinson, T.; Rogers, D.; Williams, B. Mapping tsetse habitat suitability in the common fly belt of Southern Africa using multivariate analysis of climate and remotely sensed vegetation data. Med. Veter. Ѐntomol. 1997, 11, 235–245. [Google Scholar] [CrossRef]
  32. Kobayashi, T.; Nakayama, S.; Wang, L.-M.; Li, G.-Q.; Yang, J. Socio-ecological analysis of desertification in the Mu-Us Sandy Land with satellite remote sensing. Landsc. Ecol. Eng. 2005, 1, 17–24. [Google Scholar] [CrossRef]
  33. Otsu, N. A threshold selection method from gray-level histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef]
  34. Liu, D.; Yu, J. Otsu Method and K-means. In Proceedings of the 2009 9th International Conference on Hybrid Intelligent Systems, Shenyang, China, 12–14 August 2009; Volume 1, pp. 344–349. [Google Scholar] [CrossRef]
  35. Yuan, X.-C.; Wu, L.-S.; Peng, Q. An improved Otsu method using the weighted object variance for defect detection. Appl. Surf. Sci. 2015, 349, 472–484. [Google Scholar] [CrossRef]
  36. Gandhi, G.M.; Parthiban, S.; Thummalu, N.; Christy, A. Ndvi: Vegetation Change Detection Using Remote Sensing and Gis—A Case Study of Vellore District. Procedia Comput. Sci. 2015, 57, 1199–1210. [Google Scholar] [CrossRef]
  37. Gou, L.; Zhang, J.; Li, N.; Wang, Z.; Chen, J.; Qi, L. Weighted assignment fusion algorithm of evidence conflict based on Euclidean distance and weighting strategy, and application in the wind turbine system. PLoS ONE 2022, 17, e0262883. [Google Scholar] [CrossRef]
  38. Yang, Z.; Liao, W.; Wang, H.; Bak, C.L.; Chen, Z. Improved Euclidean Distance Based Pilot Protection for Lines with Renewable Energy Sources. IEEE Trans. Ind. Inform. 2022, 3203. [Google Scholar] [CrossRef]
  39. Bozdağ, A.; Yavuz, F.; Günay, A.S. AHP and GIS based land suitability analysis for Cihanbeyli (Turkey) County. Environ. Earth Sci. 2016, 75, 813. [Google Scholar] [CrossRef]
  40. Rossiter, D.G. A theoretical framework for land evaluation. Geoderma 1996, 72, 165–190. [Google Scholar] [CrossRef]
  41. Akpoti, K.; Higginbottom, T.P.; Foster, T.; Adhikari, R.; Zwart, S.J. Mapping land suitability for informal, small-scale irrigation development using spatial modelling and machine learning in the Upper East Region, Ghana. Sci. Total Environ. 2022, 803, 149959. [Google Scholar] [CrossRef]
  42. Jiao, S.; Zhang, X.; Xu, Y. A review of Chinese land suitability assessment from the rainfall-waterlogging perspective: Evidence from the Su Yu Yuan area. J. Clean. Prod. 2017, 144, 100–106. [Google Scholar] [CrossRef]
  43. Akpoti, K.; Kabo-Bah, A.T.; Zwart, S.J. Agricultural land suitability analysis: State-of-the-art and outlooks for integration of climate change analysis. Agric. Syst. 2019, 173, 172–208. [Google Scholar] [CrossRef]
  44. Ebrahimi, M.; Nejadsoleymani, H.; Daneshvar, M.R.M. Land suitability map and ecological carrying capacity for the recognition of touristic zones in the Kalat region, Iran: A multi-criteria analysis based on AHP and GIS. Asia-Pac. J. Reg. Sci. 2019, 3, 697–718. [Google Scholar] [CrossRef]
  45. Nurda, N.; Noguchi, R.; Ahamed, T. Change Detection and Land Suitability Analysis for Extension of Potential Forest Areas in Indonesia Using Satellite Remote Sensing and GIS. Forests 2020, 11, 398. [Google Scholar] [CrossRef]
  46. Parry, J.A.; Ganaie, S.A.; Bhat, M.S. GIS based land suitability analysis using AHP model for urban services planning in Srinagar and Jammu urban centers of J&K, India. J. Urban Manag. 2018, 7, 46–56. [Google Scholar] [CrossRef]
  47. Zhou, J.; Cao, X. What is the policy improvement of China’s land consolidation? Evidence from completed land consolidation projects in Shaanxi Province. Land Use Policy 2020, 99, 104847. [Google Scholar] [CrossRef]
  48. Verstappen, H.T. Old and New Trends in Geomorphological and Landform Mapping. Dev. Earth Surf. Processes 2011, 15, 13–38. [Google Scholar] [CrossRef]
  49. Proulx, R.; Fahrig, L. Detecting human-driven deviations from trajectories in landscape composition and configuration. Landsc. Ecol. 2010, 25, 1479–1487. [Google Scholar] [CrossRef]
  50. Zhang, Z.; Huang, X.; Zhou, Y. Factors influencing the evolution of human-driven rocky desertification in karst areas. Land Degrad. Dev. 2020, 32, 817–829. [Google Scholar] [CrossRef]
  51. Li, N.; Yan, C.Z.; Xie, J.L. Remote sensing monitoring recent rapid increase of coal mining activity of an important energy base in northern China, a case study of Mu Us Sandy Land. Resour. Conserv. Recycl. 2015, 94, 129–135. [Google Scholar] [CrossRef]
  52. Gao, W.; Zheng, C.; Liu, X.; Lu, Y.; Chen, Y.; Wei, Y.; Ma, Y. NDVI-based vegetation dynamics and their responses to climate change and human activities from 1982 to 2020: A case study in the Mu Us Sandy Land, China. Ecol. Indic. 2022, 137, 108745. [Google Scholar] [CrossRef]
  53. Guo, Z.-C.; Wang, T.; Liu, S.-L.; Kang, W.-P.; Chen, X.; Feng, K.; Zhang, X.-Q.; Zhi, Y. Biomass and vegetation coverage survey in the Mu Us sandy land-based on unmanned aerial vehicle RGB images. Int. J. Appl. Earth Obs. Geoinf. ITC J. 2021, 94, 102239. [Google Scholar] [CrossRef]
  54. Ju, C.; Cai, T.; Yang, X. Topography-based modeling to estimate percent vegetation cover in semi-arid Mu Us sandy land, China. Comput. Electron. Agric. 2008, 64, 133–139. [Google Scholar] [CrossRef]
  55. Wang, L.; Wu, Q.; Fu, Z.; Jiang, G.; Liu, Y.; Xu, K. Aeolian sand accumulation and land desertification over the past 1,500 years as revealed by two aeolian dunes in the Beiluhe Basin on interior Qinghai-Tibet Plateau, SW China. Quat. Int. 2022, 613, 101–117. [Google Scholar] [CrossRef]
  56. Cao, S. Why Large-Scale Afforestation Efforts in China Have Failed To Solve the Desertification Problem. Environ. Sci. Technol. 2008, 42, 1826–1831. [Google Scholar] [CrossRef] [PubMed]
  57. Le Houérou, H.N. Climate change, drought and desertification. J. Arid Environ. 1996, 34, 133–185. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The elevation map of Mu Us Sandy land (Source: Shuttle Radar Topography Mission (SRTM)).
Figure 1. The elevation map of Mu Us Sandy land (Source: Shuttle Radar Topography Mission (SRTM)).
Remotesensing 14 04736 g001
Figure 2. The flowchart of the methodology.
Figure 2. The flowchart of the methodology.
Remotesensing 14 04736 g002
Figure 3. Annual change of segmentation threshold from 1986 to 2020.
Figure 3. Annual change of segmentation threshold from 1986 to 2020.
Remotesensing 14 04736 g003
Figure 4. Linear regression results of ENDVI threshold and vegetation area change from 1986 to 2020. (a) ENDVI threshold; (b) vegetation area changes.
Figure 4. Linear regression results of ENDVI threshold and vegetation area change from 1986 to 2020. (a) ENDVI threshold; (b) vegetation area changes.
Remotesensing 14 04736 g004
Figure 5. Spatial accessibility between each pixel and human-related features including (a) roads, (b) railways, (c) human settlements, (d) service stations, and (e) waterways.
Figure 5. Spatial accessibility between each pixel and human-related features including (a) roads, (b) railways, (c) human settlements, (d) service stations, and (e) waterways.
Remotesensing 14 04736 g005aRemotesensing 14 04736 g005b
Figure 6. Mean Euclidean distance from vegetation grid to human-related features between 2011 and 2020.
Figure 6. Mean Euclidean distance from vegetation grid to human-related features between 2011 and 2020.
Remotesensing 14 04736 g006
Figure 7. Percentages of vegetation pixels in each 1 km buffer area in 2020.
Figure 7. Percentages of vegetation pixels in each 1 km buffer area in 2020.
Remotesensing 14 04736 g007
Figure 8. Annual average slope in each vegetation pixel from 1986 to 2020.
Figure 8. Annual average slope in each vegetation pixel from 1986 to 2020.
Remotesensing 14 04736 g008
Figure 9. Suitability map for artificial vegetation restoration.
Figure 9. Suitability map for artificial vegetation restoration.
Remotesensing 14 04736 g009
Figure 10. Example of how to utilize the artificial vegetation restoration suitability map for future planning (empty region indicates the area has been covered by vegetation in 2020).
Figure 10. Example of how to utilize the artificial vegetation restoration suitability map for future planning (empty region indicates the area has been covered by vegetation in 2020).
Remotesensing 14 04736 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Z.; Huang, M.; Xiao, C.; Qi, S.; Du, W.; Zhu, D.; Altan, O. Integrating Remote Sensing and Spatiotemporal Analysis to Characterize Artificial Vegetation Restoration Suitability in Desert Areas: A Case Study of Mu Us Sandy Land. Remote Sens. 2022, 14, 4736. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14194736

AMA Style

Chen Z, Huang M, Xiao C, Qi S, Du W, Zhu D, Altan O. Integrating Remote Sensing and Spatiotemporal Analysis to Characterize Artificial Vegetation Restoration Suitability in Desert Areas: A Case Study of Mu Us Sandy Land. Remote Sensing. 2022; 14(19):4736. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14194736

Chicago/Turabian Style

Chen, Zhanzhuo, Min Huang, Changjiang Xiao, Shuhua Qi, Wenying Du, Daoye Zhu, and Orhan Altan. 2022. "Integrating Remote Sensing and Spatiotemporal Analysis to Characterize Artificial Vegetation Restoration Suitability in Desert Areas: A Case Study of Mu Us Sandy Land" Remote Sensing 14, no. 19: 4736. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14194736

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

Article Metrics

Back to TopTop