Next Article in Journal
Improving Colored Dissolved Organic Matter (CDOM) Retrievals by Sentinel2-MSI Data through a Total Suspended Matter (TSM)-Driven Classification: The Case of Pertusillo Lake (Southern Italy)
Next Article in Special Issue
Detection and Attribution of Vegetation Dynamics in the Yellow River Basin Based on Long-Term Kernel NDVI Data
Previous Article in Journal
MSSFF: Advancing Hyperspectral Classification through Higher-Accuracy Multistage Spectral–Spatial Feature Fusion
Previous Article in Special Issue
Combination of UAV Photogrammetry and Field Inventories Enables Description of Height–Diameter Relationship within Semi-Arid Silvopastoral Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A 20-Year Analysis of the Dynamics and Driving Factors of Grassland Desertification in Xilingol, China

1
State Key Laboratory of Remote Sensing Science, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
2
University of Chinese Academy of Sciences, Beijing 101408, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(24), 5716; https://doi.org/10.3390/rs15245716
Submission received: 1 November 2023 / Revised: 6 December 2023 / Accepted: 10 December 2023 / Published: 13 December 2023
(This article belongs to the Special Issue Remote Sensing of Arid/Semiarid Lands II)

Abstract

:
Grassland desertification stands as an ecological concern globally. It is crucial for desertification prevention and control to comprehend the variation in area and severity of desertified grassland (DGL), clarify the intensities of conversion among DGLs of different desertification levels, and explore the spatial and temporal driving factors of desertification. In this study, a Desertification Difference Index (DDI) model was constructed based on albedo-EVI to extract desertification information. Subsequently, intensity analysis, the Geo-detector model, and correlation analysis were applied to analyze the dynamics and driving factors of desertification. The results showed the following: (1) Spatially, the DGL in Xilingol exhibited a zonal distribution. Temporally, the degree of DGL decreased, with the proportion of severely and moderately desertified areas decreasing from 51.77% in 2000 to 37.23% in 2020, while the proportion of nondesertified and healthy areas increased from 17.85% in 2000 to 37.40% in 2020; (2) Transition intensities among different desertification levels were more intense during 2000–2012, stabilizing during 2012–2020; (3) Meteorological factors and soil conditions primarily drive the spatial distribution of DDI, with evapotranspiration exhibiting the most significant influence (q-value of 0.83), while human activities dominate interannual DDI variations. This study provides insights into the conversion patterns among different desertification levels and the divergent driving forces shaping desertification in both spatial and temporal dimensions in Xilingol.

Graphical Abstract

1. Introduction

Desertification denotes the degradation of land in arid, semiarid, and semihumid regions, resulting from a combination of factors including climate change and human activities [1]. With social development and climate change, desertification has become an increasingly severe issue and is now one of the most serious challenges faced by humanity [2]. Statistics showed that area of desertification was approximately 3.6 × 107 km2 in 2009 globally, constituting 25% of the Earth’s total land area. Currently, land desertification impacts two-thirds of the world’s nations and nearly 20% of the global population [3]. In the 21st century, land desertification has emerged as a major environmental challenge, significantly hindering the prospects of sustainable development [4]. China stands as one of the countries most significantly affected by desertification on a global scale, experiencing the largest land area impacted, the largest population affected, and the most severe wind and sand hazards. Presently, the extent of desertified land in China encompasses 1.72 × 106 km2, representing 17.93% of the entire national land area. Desertification is distributed across 920 counties in 30 provinces (autonomous regions or municipalities) [5]. The top five provincial-level administrative regions with the largest desertification areas are Xinjiang, Inner Mongolia, Tibet, Qinghai, and Gansu. Among these regions, Inner Mongolia’s desertification area comprises 23.7% of the national desertified land, positioning it as one of the areas most heavily impacted by desertification [6]. Furthermore, Inner Mongolia’s grassland, as China’s largest grassland, plays a crucial role as a significant livestock production base. The management and protection of grasslands are of paramount importance for maintaining China’s ecological security [7]. Consequently, the control of desertified grassland (DGL) has garnered widespread attention from the government and various sectors of society.
Confronted with the escalating challenge of grassland desertification, precise, timely, and efficacious monitoring of DGL assumes paramount significance in comprehending the dynamics of the desertification processes and devising strategies for its prevention and control [8]. Field surveys were once the main method for monitoring vegetation changes but were challenging to conduct over long time periods and large spatial scales. Additionally, they required substantial human, financial, and material resources [9]. The application of remote sensing has provided a new possibility for continuous monitoring of DGL over extended periods and large spatial extents [10,11].
At present, the extraction of remote sensing information for desertification primarily relies on visual interpretation or automatic classification using machine learning methods [12]. This process is integral to the monitoring, mapping, analysis, and assessment of desertification. Traditional visual interpretation methods are not only susceptible to human factors but also entail high work intensity and low efficiency. The application of automatic classification methods in large-scale regions is constrained by limitations on classification accuracy. Simultaneously, conventional remote sensing information extraction methods exhibit a low level of utilization, thereby impeding the optimal contribution of abundant remote sensing information to desertification monitoring. Zeng et al. introduced the Desertification Difference Index (DDI) model [13], which utilizes the albedo-NDVI feature space to streamline and improve the monitoring of desertification on a large scale. This model, characterized by easy data acquisition and good temporal and spatial continuity, resolved the problem of data accessibility in most desertification index systems. Therefore, it has been widely adopted by researchers [14,15,16,17,18,19].
However, due to nonlinear stretching, NDVI suppresses high-value areas, making it prone to saturation [20]. Furthermore, NDVI only considers the influence of near-infrared and red bands, failing to effectively control the impact of atmospheric and soil factors on vegetation reflectance spectra. The Enhanced Vegetation Index (EVI) has improved and optimized NDVI, further mitigating the effects of atmospheric and soil factors [21,22]. Moreover, existing studies predominantly explore the drivers of desertification from either a temporal or spatial dimension alone, with limited comprehensive analyses of desertification drivers from both temporal and spatial perspectives. Thus, this study employed the EVI to construct the DDI and comprehensively investigated the drivers of desertification in the Xilingol grassland from both spatial distribution and temporal change perspectives. This approach aims to provide the public with a clearer and more comprehensive understanding.
The factors influencing desertification are complex and diverse. Studies have confirmed the substantial influence of meteorological factors on land desertification, including temperature, precipitation, wind speed, and evapotranspiration [14,23,24,25,26,27,28]. Regions prone to desertification generally experience large temperature differences, scarce precipitation, and high wind speeds, which directly or indirectly lead to sparse vegetation and fragile ecological environments, hindering effective desertification prevention [29]. Apart from meteorological factors, various environmental factors such as different topographies, vegetation cover types, soil types, and surface soil moisture levels can also influence the occurrence of land desertification [30,31,32]. Moreover, human activities, including overgrazing, intensive cultivation, and extensive industrial land development, significantly contribute to land desertification [33,34,35].
This study employed intensity analysis and a geographical detector model to investigate the dynamics of desertification and the spatial distribution drivers of desertification. In previous research, land-use dynamic models and transition matrices have commonly been used to elucidate desertification dynamics [36,37,38]. However, these methods require the separate determination of mutual conversions between different land categories using transition matrices, which hinders a systematic and quantitative evaluation of desertification dynamics. In contrast, intensity analysis demonstrates significant advantages in terms of systematic and quantitative assessments, making it widely utilized by researchers [39,40,41,42,43]. The Geo-detector model, a statistical method for exploring spatial heterogeneity and its driving factors, not only investigates the individual impact of a single factor on desertification but also detects the interaction effects between two variables influencing desertification. This model is widely utilized in research pertaining to economics, population dynamics, land use, and other fields [26,39,40,41,43].
This study aimed to do the following: (1) quantify the variations in DGL area and analyze the spatial patterns of DGL, (2) elucidate the processes of change and conversion patterns in the DGLs of different levels, and (3) explore and compare the driving factors of the DGL in Xilingol from 2000 to 2020 from both temporal and spatial dimensions. This study holds significant value in enhancing our understanding of the spatial distribution characteristics, driving factors, and interactions related to grassland desertification. It provides valuable insights for subsequent research and desertification prevention efforts.

2. Materials and Methods

2.1. Study Area

Xilingol League is located in central Inner Mongolia, spanning from 41.4°N to 46.6°N and from 111.1°E to 119.7°E (Figure 1), with a total area of 206,000 km2. The region displays a temperate continental climate characterized by arid and semiarid features, including an average annual temperature of 2.2 °C, approximately 280 mm of average annual precipitation, and elevations ranging from 760 m to 1925 m above sea level. The terrain is predominantly plateau, sloping from south to north, featuring low hills and basins in the eastern and southern sectors. A distinct zonal pattern is observed in soil and vegetation coverage, with soils transitioning from black calcareous soil to light and dark chestnut soils from southeast to northwest. This transition aligns with noticeable variations in vegetation coverage, progressing from forests to meadow steppe, typical steppe, and desert steppe. Grasslands, covering around 87% of the entire region, dominate the landscape [43,44,45].
Xilingol League comprises twelve county-level administrative regions, with a total resident population of 1.1071 million. The twelve banners and counties include Xilinhot City (XC), Erenhot City (EC), East Ujimqin Banner (EUB), West Ujimqin Banner (WUB), Abag Banner (AB), Sunite Left Banner (SLB), Sunite Right Banner (SRB), Xianghuang Banner (XHB), Zhengxiangbai Banner (ZXBB), Zhenglan Banner (ZLB), Taibus Banner (TB), and Duolun County (DLC) [43].

2.2. Data and Processing

The data used in this study can be categorized into six types, including vegetation index, reanalysis meteorological data, soil data, topographical data, human activities data, and statistics data. The specific data sources and details are provided in Table 1. Regarding the EVI data, the study utilizes the annual maximum values, while precipitation is represented by the annual cumulative values. Surface temperature, wind speed, evapotranspiration, surface albedo, and soil moisture are all represented by annual averages. For population density, the average values for the years 2000–2020 were calculated. Spatial distribution data for livestock density were available only for the years 2020 and 2015; hence, the livestock density data used in this study represent the average for these two years. Additionally, population, livestock numbers (end of June) and GDP statistics for Xilingol from 2000 to 2020 were obtained by referring to the Xilingol Statistical Yearbook. These data have been well applied in relevant studies conducted in the Inner Mongolia region [14,26,40,43,44,46].
To ensure consistency in spatial resolution for input data, all geographical spatial data in this study were resampled to 250 m. In terms of statistical data, a normalization process was applied to the acquired data.

2.3. Methods

The workflow of this study comprised five main components (Figure 2). Firstly, data input, which encompassed topographic data, vegetation index, geography-based human activities data, soil data, data from statistical yearbooks, and reanalysis meteorological data. Secondly, data preprocessing was conducted, involving reprojection, resampling, cropping, calculation, and standardization. Thirdly, the Enhanced Vegetation Index (EVI) and surface albedo were utilized to construct the DDI, with subsequent threshold delineation of the DDI results. Finally, based on the DDI and threshold delineation outcomes, the study employed intensity analysis, the geographical detector model, and correlation analysis to analyze the spatiotemporal dynamics and drivers of desertification in the Xilingol grassland.

2.3.1. Desertification Difference Index (DDI) Model Construction

According to the DDI model proposed by Zeng et al. [13], albedo-NDVI exhibit a significant negative linear correlation in a one-dimensional feature space [13,15,16,17,18]. In this study, the model was enhanced by replacing NDVI with EVI (Figure 3), effectively reducing the impact of atmosphere and soil. A higher DDI value indicates a lower degree of desertification. The basic equations for this modified model are presented in Equations (1)–(3):
DDI = α·EVI-Albedo
Albedo = k·EVI + b
α k = 1
In the equations, k is the slope, b is the constant term, and α is the reciprocal of k.
Researchers tend to employ the natural breakpoint method to classify DDI values after calculation of DDI [14,18]. While this approach is feasible when analyzing data for a specific year, it lacks meaningful implications for long-term data comparisons or trend analyses. Therefore, in this study, a uniform threshold for desertification classification was determined based on human–machine interaction, as shown in Table 2. Subsequently, 200 random points within the study area were selected, and the DDI results for the years of 2006, 2012, and 2020 were validated using Google Earth historical images. The accuracy rates were 82%, 85.5%, and 80% respectively, meeting the experimental requirements.

2.3.2. Data Normalization and Weighting

To explore the driving forces of temporal variations in DDI, factors were categorized into two groups, and normalization was performed on each group. Assuming that the impact weights of each factor on the composite value in both groups are equal, in this study, we employed the method of average weighting to calculate the composite values for meteorological and human activity factors, respectively. The calculation is formulated as shown in Equations (4) and (5):
x n = x i x m i n x m a x x m i n ( i = 1 , 2 , , n )
X = 1 n ( x 1 + x 2 + + x n )
In the equations, x n is the normalized value of each factor, x i is the original value, x m i n is the minimum value, x m a x is the maximum value, and X is the composite value of each group.

2.3.3. Spearman Correlation Coefficient

The study employed the Spearman correlation coefficient to analyze the correlation between variables. DDI averages for each two-year period from 2000 to 2020 were computed, with values for intervening years derived as the mean of adjacent two-year DDI values. To assess the predominant influence of meteorological factors versus human activities on DDI changes, separate weighted averages were calculated for both sets of factors. Subsequent testing revealed nonnormal distributions for both sets of values, leading to the selection of the Spearman correlation coefficient for correlation analysis [47,48]. These were calculated according to Equations (6) and (7):
ρ = 1 6   d i 2 n ( n 2 1 )
d i = x i y i
In these equations, ρ is the Spearman correlation coefficient, n is the sample size, and x i and y i are the variables.

2.3.4. Intensity Analysis

The intensity analysis method is primarily employed for quantitatively analyzing variations in changes among diverse land categories [42,49]. Intensity analysis involves three levels: interval level, category level, and transition level [42].
The mathematical notations used in intensity analysis include the following: i, representing a category at the initial time point within a specific time period; j, representing a category at the final time point within a specific time period; J, denoting the total number of categories; m, representing the losing category in a specified transition; n, representing the gaining category in a specified transition; T, denoting the number of time point; and t, representing the initial time point within [Yt, Yt+1], with t ranging from 1 to T − 1.
  • Interval level
Two indicators are introduced at this level: St, representing the annual percentage of actual changed area over the total area within each time interval, denoting the average annual change rate over T − 1 intervals; and U, representing the annual percentage of theoretical changed area over the total study area per time interval, spanning the entire study period. The terms “slow” and “fast” were employed to characterize the relative rates of change within these time intervals. These were calculated according to Equations (8) and (9):
S t = { j = 1 J [ ( i = 1 J c tij ) C t i i ] } / [ j = 1 J ( i = 1 J C t i j ) ] Y t + 1 Y t × 100 %
U = t = 1 T 1 { j = 1 J [ ( i = 1 J c tij ) C t i i ] } / [ j = 1 J ( i = 1 J C t i j ) ] Y T Y 1 × 100 %
2.
Category level
At this level, two indicators are employed. Firstly, the reduction intensity L t i signifies the annual percentage of land class i area reduced from the initial total area of class i during the study period [ Y t , Y t + 1 ] . Secondly, the increase intensity G t j of each land class j represents the annual percentage of area added to class j from the final total area of class j within the study period [ Y t , Y t + 1 ] , indicating changes from non-j land use types to land class j. L t i and G t j are calculated according to Equations (10) and (11).
L t i = [ (   j = 1 J c t i j ) C t i i ] / ( Y t + 1 Y t )   j = 1 J C t i j × 100 %
G t j = [ (   i = 1 J C tij ) C t j j ] / ( Y t + 1 Y t )   i = 1 J C t i j × 100 %
3.
Transition level
At this level, four change indicators are derived from two perspectives. Firstly, an increase in land class n is assessed using observed transition intensities R t i n from land class i to land class n, where J − 1 transfer intensities exist from other land classes to land class n within each time interval. The average increase intensity W t n of land class n is calculated, indicating the average transfer intensity from other land classes to land class n, weighted by their respective areas. Secondly, a decrease in land class m is evaluated through observed transition intensities Q t m j from land class m to land class j. The average decrease intensity V t m of land class m is computed. R t i n , W t n , Q t m j , and V t m are calculated according to Equations (12)–(15):
R t i n = C t i n ( Y t + 1 Y t ) j = 1 J C t i j × 100 %
W t n = [ ( i = 1 J C t i n ) C t n n ] ( Y t + 1 Y t ) j = 1 J [ ( i = 1 J C t i j ) C t n j ] × 100 %
Q t m j = C t m j ( Y t + 1 Y t ) i = 1 J C t i j × 100 %
V t m = [ ( j = 1 J C t m j ) C t m m ] ( Y t + 1 Y t ) i = 1 J [ ( j = 1 J C t i j ) C t i m ] × 100 %

2.3.5. Geographical Detector

The geographical detector is a powerful tool for conducting driver and factor analyses [50,51]. The model comprises four detectors: risk detector, factor detector, ecological detector, and interaction detector. In this study, the factor detector, ecological detector, and interaction detector were employed [42].
  • Factor Detector
The q-value (ranging from 0 to 1) quantifies how well factor x explains the spatial variation of attribute y. A higher q-value indicates stronger explanatory power, determined through a significant F-test.
q = 1 h = 1 L N h σ h 2 N σ 2
Within the equation, “h = 1 … L” signifies distinct strata or layers for both the independent variable x and the dependent variable y, essentially representing categories or partitions. Additionally, N h stands for the quantity of units in layer h, while N denotes the total number of units across the entire region. Furthermore, σ h 2 represents the variance of y values within layer h, and σ 2 denotes the overall variance of y values across the entire region.
2.
Ecological Detector
Ecological detector is employed to determine if notable distinctions exist in the spatial distribution impacts on attribute y between two factors. This evaluation is quantified through the F-statistic. These are calculated according to Equations (17) and (18):
F = N x 1 ( N x 2 1 ) SSW x 1 N x 2 ( N x 1 1 ) SSW x 2
SSW x 1 = h = 1 L 1 N h σ h 2 N σ 2 , SSW x 2 = h = 1 L 2 N h σ h 2 N σ 2
In these equations, N x 1   and   N x 2 are the number of factors for x1 and x2, respectively; S S W x 1   and   S S W x 2 are the accumulation of squares within each grouping established by the two factors; and L 1   and   L 2 signify the count of classes or categories for x1 and x2, respectively.
3.
Interaction Detector
Interaction detector assesses interactions between risk factors x, determining whether x1 and x2 jointly influence the explanatory power of dependent variable y, or if their effects on y are independent. It calculates q-values for x1 and x2 (q(x1) and q(x2)) and the interaction q-value (q(x1x2)) formed by the interplay of x1 and x2. Comparing q(x1x2) with q(x1) and q(x2) reveals the interaction type. Further details on interaction detector can be found in Table 3.

3. Results

3.1. Analysis of DDI Distribution and Spatiotemporal Trends

3.1.1. Spatiotemporal Distribution Characteristics of DDI

In this study, we calculated the DDI value for Xilingol every two years, and the values were categorized into five groups (Figure 4). Spatially, desertification levels exhibited a systematic east-to-west intensification, aligning with the distribution of meadow steppe, typical steppe, and desert steppe in Xilingol (Figure 5). Temporally, proportions of DGL levels were calculated for each period (Figure S1). It was observed that the area of level I showed a fluctuating decreasing trend, declining from 22.03% in 2000 to 15.15% in 2020, and the proportion of level II areas decreased from 29.74% in 2000 to 22.08% in 2020. In contrast, the areas of level IV and level V increased annually. The combined proportion of these two levels rose from 17.85% in 2000 to 37.40% in 2020, and the combined area of the two was more than double the initial total. These results indicate an amelioration in the degree of grassland desertification over the course of 20 years.

3.1.2. Spatiotemporal Variation of DGL

A statistical analysis was conducted on areas exhibiting various degrees of desertification (Table S1). From 2000 to 2010, level I increased from 44,518.60 km2 to 56,903.04 km2, indicating a worsening trend, particularly intensifying from 2004 to 2010. During this period, levels I, II, and III expanded from 28,395.93 km2, 44,344.86 km2, and 39,155.07 km2 in 2004 to 56,903.04 km2, 48,520.86 km2, and 55,841.82 km2 in 2010, respectively.
Contrastingly, from 2010 to 2020, desertification showed a decline. The most substantial decrease occurred in level I, decreasing from 56,903.04 km2 to 30,606.31 km2. Simultaneously, levels IV and V significantly increased from 29,872.40 km2 and 10,912.78 km2 in 2010 to 54,934.55 km2 and 20,633.46 km2 in 2020, respectively.
Across 2000–2020, levels I, II, and III areas decreased, while levels IV and V increased significantly. Notably, level IV expanded from 27,780.69 km2 in 2000 to 53,273.86 km2 in 2020, a 91.78% increase. Level V grew from 8281.21 km2 to 22,428.92 km2 in 2020, nearly tripling its original size.
Linear fitting was employed to observe changes in different degrees of grassland desertification (DGL) areas (Figure 6). Notably, levels I, II, and III decreased, while Levels IV and V increased. Examining the absolute values of the fitted lines’ slopes revealed change rates, which ranked as follows: level V (1115.49) > level I (1109.25) > level IV (832.93) > level II (620.98) > level III (217.66). Consequently, it can be inferred that from 2000 to 2020, level V increased the fastest, while level I decreased the fastest, indicating an improvement in the desertification situation of Xilingol.

3.2. Intensity Analysis

In this study, based on the years of DDI calculation and to accurately capture the dynamics of desertification changes, the period from 2000 to 2020 was divided into five cycles for intensity analysis. Intensity analysis comprises three levels: interval level, category level, and transition level.

3.2.1. Interval Level

At the interval level (Figure 7), the sections to the left and right of the zero value depict the proportion of the total area changed over the cycle and the proportion of average annual change area during the cycle, respectively. This illustrates the overall conversion intensity among different desertification levels within a specific period. The red dashed line serves as a uniform intensity indicator, representing the overall change intensity from 2000 to 2020; values surpassing this threshold indicate fast conversion, while values below it suggests slow conversion.
It is evident that the intensity during the periods 2000–2004, 2004–2008, and 2008–2012 significantly exceeded the uniform value. Consequently, during these three periods, there were drastic changes among different desertification levels, with the period from 2000 to 2004 exhibiting the highest intensity, while the intensity was relatively mild from 2016 to 2020.

3.2.2. Category Level

At the category level (Figure 8), the values to the left of zero represent the increase in area from other levels of DGL converted to the specific level during the period and the decrease in area from the specific level of DGL converted to other levels during the period. The values to the right of zero represent the average annual gain and loss intensity for the specific level of DGL during the period. The red dashed line denotes the uniform intensity, determined by the average annual change intensity S from the interval level. If a bar surpasses the uniform line, it indicates a relatively active change in the DGL category during the specified time interval; otherwise, the change is relatively steady.
It can be observed that during the 2000–2004 period (Figure 8A), the loss intensities for levels III and IV exceeded the average, and the gain intensities for levels III, IV, and V surpassed the average. This suggests intense conversions among levels III, IV, and V during this period. In the 2004–2008 period (Figure 8B), the gain and loss intensities for almost all desertification levels were close to or greater than the average, making it one of the most active periods for conversions among DGL categories. In the 2008–2012 period (Figure 8C), the gain intensities for levels IV and V were significantly greater than their respective loss intensities, while the loss intensity for level I was significantly higher than the gain intensity. This indicates a notable improvement in the desertification situation during this period. In the 2012–2016 period (Figure 8D), the gain intensities of levels I, II, and III were greater than the loss intensities, while the loss intensities for levels IV and V were greater than the gain intensities, suggesting an exacerbation of desertification during this period. In the 2016–2020 period (Figure 8E), the loss intensities for levels I and III were greater than the gain intensity, and only the gain intensity for level IV and the loss intensity for level III significantly exceeded the average, indicating a stable and positive trend during this period.
In summary, the DGL categories exhibited acute changes in intensity during the 2000–2004, 2004–2008, and 2008–2012 periods. Desertification trends were positive during the 2000–2004 and 2008–2012 periods, but deteriorated during the 2004–2008 and 2012–2016 periods.

3.2.3. Transition Level

At the transition level, the study analyzed statistical data pertaining to the changes in grassland areas during a specific period (Figure 9). With Figure 9A taken as an example, values to the left of zero denote the proportion of the level I area that originated from levels II, III, IV, and V, while values to the right represent the portion of the level I area that transformed into levels II, III, IV, and V.
Figure 9B indicates that during the 2004–2008 and 2012–2016 periods, among all areas transitioning out of level II, the majority shifted to level I. Simultaneously, in all areas transitioning into level II, the largest portion came from level III, which converted into level II. Similar trends can be observed in Figure 9C,D during the same time periods, signifying a consistent pattern of land transitioning from lower to higher desertification levels. This suggests an exacerbation of desertification trends during 2004–2008 and 2012–2016.
In terms of transition intensity (Figure 10), with Figure 10A taken as an example, values to the left of zero signify the intensity of transition from level I to other DGL levels, whereas values to the right indicate the intensity of transition from other levels to level I. The red dashed line to the left illustrates the average transition intensity from level I to other levels, and the corresponding line to the right represents the average transfer intensity from other levels to level I. Figure 10 presents a specific range of values, and values outside this range are annotated.
From Figure 10A, it can be observed that the intensity values of transitions from level I to level II and from level II to level I were the maximum values among all transition intensities involving level I. Specifically, during the 2000–2004, 2008–2012, and 2016–2020 periods, the intensity of transitions from level I to level II was greater than the intensity of transitions from level II to level I, indicating a decrease in the area of level I and an increase in the area of level II during these three periods. A similar situation is evident in Figure 10C, where the intensity of transitions from level III to level IV and level V is greater than the intensity of transitions from level IV and level V to level III, suggesting an increase in the areas of level IV and level V.
In summary, there was an overall improvement in desertification trends during 2000–2004, 2008–2012, and 2016–2020. However, there was a certain degree of exacerbation of desertification during 2004–2008 and 2012–2016, with level III exhibiting the highest intensity of gain and loss. Notably, the transition intensities between adjacent levels consistently surpassed the average transfer intensity. This pattern underscores that grassland desertification is a gradual process, characterized by shifts between neighboring levels.

3.3. Driving Factors of Grassland Desertification

3.3.1. Driving Factors of DDI Spatial Distribution

In the model inputs, the dependent variable (y) was represented by the calculated average value of the multiyear DDI. The independent variables (x) encompassed four categories of factors constituting a total of 14 factors: reanalysis meteorological data, soil data, topographic data, and human activity data (Figure 11). Notably, land surface temperature, precipitation, wind velocity, evapotranspiration, soil moisture, and population density were represented by average values spanning 2000 to 2020. Conversely, soil type, soil erosion type, soil erosion intensity, elevation, slope, and aspect were treated as single values, as these factors were not expected to undergo significant changes during the study period. Geographic distribution data for livestock load employed the average values of currently available data, while land use and cover changes (LUCC) were based on the latest classification results from the IGBP (International Geosphere-Biosphere Programme) in the year 2020.
Within the study area, 3000 random points were generated, and numerical values for each factor and DDI were meticulously extracted at these points. Following thorough data validation, these values were input into the Geo-detector model for rigorous calculations.
In factor detector, it was observed that meteorological factors and soil factors, on the whole, exhibited significant explanatory power concerning the spatial distribution of DDI (Figure 12). Among these factors, evapotranspiration displayed the highest q-value, reaching 0.83, followed by precipitation, soil type, and land surface temperature, which also played crucial roles in driving DDI distribution. In contrast, aspect exhibited the weakest explanatory power.
Ecological detector is used to determine whether there are significant differences in the spatial distribution impacts of each pair of independent variables on the dependent variable. The results indicated that (Figure 13) concerning surface temperature, notable distinctions were present in its impact on the spatial distribution of the DDI when compared to precipitation, evapotranspiration, and soil type. Additionally, significant variations were observed in the influence of precipitation and evapotranspiration on the spatial distribution of DDI. In the context of wind velocity, there were significant differences in its impact on the spatial distribution of DDI in comparison to evapotranspiration, soil type, soil moisture, soil erosion type, elevation, and livestock. The effects on the spatial distribution of DDI also exhibited significant disparities for the soil erosion intensity–DEM, slope–livestock, and aspect–livestock factor pairs. For other factor pairs, there was no significant difference in their impact on the spatial distribution of DDI.
Interaction detector (Figure 14) indicated that the combined effect of any two factors in this study enhanced the explanatory power of the spatial distribution of DDI. Among them, the joint effect of evapotranspiration and soil type exhibited the strongest explanatory power on the spatial distribution of DDI, reaching 0.87, followed by the combined effect of evapotranspiration and land surface temperature, which reached 0.85. Overall, the combination of water, heat, and soil exerted a dominant driving role in the spatiotemporal distribution of desertification.

3.3.2. Driving Factors of DDI Temporal Variation

To explore the driving factors of DDI temporal variation, seven factors were selected and categorized into meteorological factors and human activities. Meteorological factors included annual average land surface temperature (LST), average annual precipitation, average wind velocity, and average evapotranspiration in Xilingol each year. Human activities included population, GDP, and livestock headcount (end of June). The study analyzed the numerical changes and trends of these factors from 2000 to 2020 (Figure 15). As depicted in Figure 15, the DDI mean value showed an increasing trend from 2000 to 2020, indicating a positive direction in desertification to some extent. Both population and GDP exhibited noticeable upward trends, while livestock headcount exhibited significant fluctuations with an overall decreasing trend, possibly due to local grazing control policies. The trends for LST and wind velocity changes were unclear, while precipitation and evapotranspiration showed significant increasing trends.
These two sets of data were then separately correlated with DDI to explore which set exerted a stronger driving force on the temporal variation in DDI. The results, as shown in Table 4, indicated that human activities exerted a stronger driving force on the temporal variation in DDI, with a significance level of 0.012, which is less than 0.05.

4. Discussion

4.1. Dynamics of Grassland Desertification in Xilingol

This study utilized the albedo-EVI to construct a DDI model and conducted a spatiotemporal analysis of desertification dynamics on the Xilingol grassland. Furthermore, quantitative analysis of desertification intensity was performed using intensity analysis. These findings are of significant importance for understanding the trends in DGL variations in the Xilingol and for implementing effective desertification prevention measures.
Presently, research employing the albedo-EVI model for DDI construction is limited. Existing findings indicate that NDVI exhibits saturation in regions with high vegetation [52], impeding the model’s ability to accurately discern actual vegetation information in these areas and consequently diminishing its capacity to precisely evaluate desertification degrees. With an increase in vegetation cover, particularly in moderately desertified areas, NDVI is subject to soil factors, resulting in the oversight of some newly added vegetation information [53]. This oversight leads to misclassification, wherein the model erroneously identifies moderately desertified areas as severely desertified. Consequently, in comparison to the albedo-NDVI model, the albedo-EVI model demonstrates enhanced accuracy in extracting desertification information.
From Figure 6 and Figure S1 and Table S1, it can be seen that desertification degree on the Xilingol grassland improved between 2000 and 2020. This finding aligns with the results reported by Batunacun et al. and Li et al. [54,55]. However, due to differences in the selected indices and threshold delineation used to characterize the degree of desertification, the classification results exhibit variations. For instance, some studies utilize NDVI values or vegetation coverage to represent land degradation or the degree of desertification [26,40,43,44]. It is essential to clarify that the primary focus of these studies is not on achieving precise classification of desertified land but rather on reflecting the spatiotemporal trends of desertification, and this study shares a similar focus. Furthermore, there are currently no established national or industry standards for the classification of lands with different levels of desertification. As a result, this study did not delve into an in-depth investigation of desertification level classification.
The intensity analysis results show that at the interval level, dynamic changes in different levels of DGL have gradually stabilized. At the category level, conversions among different levels of DGL were particularly intense during 2004–2008, indicating an increasing trend in grassland desertification. Transitions at the desertification level mainly occurred between adjacent levels, with substantial year-to-year variations, possibly attributed to meteorological conditions. In the context of global climate change and frequent extreme weather events, local microclimates in Xilingol are significantly influenced, resulting in unstable annual precipitation [56]. The grassland’s responsiveness to water availability induces rapid improvements in grassland quality following rainfall, consequently reducing the degree of desertification.

4.2. Driving Factors of Grassland Desertification in Xilingol

Geodetector and correlation analyses were used to comprehensively explore the driving factors influencing spatial and temporal variations in desertification on the Xilingol grassland. Unlike some prior studies that focused on specific factors, this research systematically classified driving factors into meteorological, soil, topographical, and human activity categories. Surface soil moisture and soil erosion factors, often overlooked in previous studies, were found to impact the spatial distribution of desertification. Additionally, livestock data included goats, sheep, cattle, and horses, recognizing their significant influence on grassland quality, thereby enhancing the scientific rigor of the analysis.
Spatially, the factor detector results showed that evapotranspiration was the most significant driving factor for the spatial distribution in desertification on the Xilingol grassland, which aligns with the results of Han et al. [28]. Since most studies did not consider ET as one of the influencing factors in the Geo-detector model, there are limited experimental results available for comparison and validation. However, evapotranspiration is strongly correlated with temperature and precipitation. In this study, precipitation, soil type, and land surface temperature were also major driving factors for DDI spatial distribution, which is consistent with the results of Wang et al., Xu et al., and Zhu et al. [25,40,44].
Precipitation, soil type, and land surface temperature emerged as pivotal factors influencing desertification. Insufficient precipitation induces soil aridity, impeding vegetation growth, escalating soil wind erosion and erosion, and intensifying desertification [57,58]. Different soil types exhibit varying water retention and wind erosion resistance capacities, influencing the rate of desertification [58,59]. Elevated land surface temperatures accelerate soil moisture evaporation, adversely affecting vegetation growth, and potentially triggering meteorological and climatic changes that exacerbate desertification trends [59,60]. The interplay of these three factors collectively contributes to vegetation degradation and soil impoverishment, accelerating the progression of desertification. A comprehensive consideration of these factors is crucial for scientifically preventing and mitigating desertification.
Meteorological factors, particularly evapotranspiration and precipitation, were identified as the most significant driving factors for spatial distribution differences in grassland desertification. Figure 11B illustrates minimal precipitation in the western part of Xilingol, indicating more severe drought conditions compared to the eastern areas. Due to water scarcity, the western region experiences poorer vegetation growth, resulting in lower vegetation cover and evapotranspiration, thereby leading to more severe desertification.
In terms of temporal variation, Spearman correlation analysis results indicated that human activities were the main driving factors for temporal changes in desertification on the Xilingol grassland. Similar conclusions were drawn by Li and Xie [36], emphasizing the crucial role of human activities in temporal desertification variation. The impact of human activities could be positive or negative. Figure 15 shows that from 2000 to 2020, both the population and GDP in the Xilingol region exhibited an overall upward trend, suggesting intensified land use and potential land desertification. However, the urbanization rate in Xilingol increased from 35.60% in 2000 to 73.88% in 2020. This indicates that during this period, a significant number of herders moved to cities, potentially reducing human interference with the grassland and contributing to the alleviation of grassland desertification. This aspect requires further investigation.
In recent years, China has actively promoted various ecological projects, including the project of “Integrated management of mountains, rivers, forests, farmlands, lakes, grasslands, and deserts”. The importance of grassland management has been emphasized in governmental documents, in which grasslands are described as the “skin of the earth” [61]. Xilingol is a vital ecological conservation area in northern China. Ecological condition of the grassland has gradually improved with the efforts of the official and public. This improvement in grassland quality aligns with the results obtained from intensity analysis.

4.3. Limitations and Future Work

This study has certain limitations that should be acknowledged. In terms of data selection, despite various driving factors being considered, comprehensiveness remains an issue. Notably, policy factors, such as rotational grazing, forbidden grazing, or grazing suspension, which can significantly impact grassland desertification in Xilingol, were not included due to challenges in quantification and incorporation into spatial distribution data suitable for the Geo-detector analysis.
Data utilization, particularly regarding livestock data, poses limitations, as the majority of such data were sourced from statistical yearbooks lacking rasterized geographic data products. Given the spatially varying distributions required for geographic detectors, this study relied on mean values from only two obtainable years, potentially influencing the results.
In the analysis of spatial changes and driving forces, the human activity factors derived from statistical yearbook data lacked uniformity. Although justified for certain reasons, such as limitations in data format and accuracy concerns, these uncertainties may impact the results and represent an area for improvement in future research.
Additionally, uncertainty exists in the method of threshold division. The absence of industry standards for defining desertification levels leads to variations in results as scholars define thresholds based on their own research. Addressing this issue is essential in future work.
The validation of desertification classification in this study relied on visual comparison using Google historical imagery, a method subject to subjectivity. Future research should consider validation through more objective approaches such as field investigations.
In future work, addressing the limitations of the current study requires comprehensive efforts. First, there is a need to incorporate policy factors, such as rotational grazing and grazing suspension, into spatial distribution data to better understand their impact on grassland desertification. Second, efforts should be directed toward obtaining or developing rasterized geographic data products to enhance the accuracy and reliability of livestock data. Third, exploring alternative sources or methodologies for deriving human activity factors can address the lack of uniformity in data. Establishing industry standards for desertification level definition is crucial for ensuring consistency across studies. Lastly, improving the validation process of desertification classification through more objective methods, such as field investigations or remote sensing techniques, will contribute to the overall robustness of future research in this domain.

4.4. Contributions and Significance

Through quantifying variations in desertified grassland (DGL) area, analyzing spatial patterns, and elucidating conversion processes among different desertification levels, this study provides crucial insights into the dynamics of desertification in both temporal and spatial dimensions. These findings offer valuable information for targeted interventions in environmental protection, aiding in the prevention and control of desertification in Xilingol by identifying areas with varying levels of desertification and understanding the intensity and patterns of conversion. The study’s emphasis on meteorological factors, soil conditions, and human activities as driving forces underscores the importance of adopting sustainable land management practices. Regarding climate change, the identification of evapotranspiration as a significant driver highlights the interconnectedness of land-use changes and climate dynamics [62]. This research contributes to the broader discourse on environmental sustainability and climate resilience, providing a foundation for evidence-based policies and practices to address grassland desertification amid global environmental challenges.

5. Conclusions

This study systematically examined the spatiotemporal dynamics of the Desertification Difference Index (DDI) in Xilingol over the 2000–2020 period. The investigation employed intensity assessments on grassland exhibiting diverse levels of desertification, considering intervals, categories, and transition levels. Geodetector and correlation analyses were applied to elucidate the driving forces behind the spatial distribution and temporal variation of DDI. The results indicated that the overall degree of desertification in Xilingol showed an improvement from 2000 to 2020.
Analysis revealed a discernible zonal pattern in the spatial distribution of DDI, exhibiting a gradual decrease from east to west, aligning with the distinct distribution of grassland types in Xilingol—specifically, meadow steppe, typical steppe, and desert steppe. Over the study duration, DDI demonstrated a consistent increasing trend, corresponding to a reduction in the extent of level I and level II areas, juxtaposed with an expansion in areas characterized as level IV and level V. Notably, the proportion of level I and level II decreases from 51.77% in 2000 to 37.23% in 2020, while the proportion of level IV and level V increased from 17.85% in 2000 to 37.40% in 2020. However, it is important to mention that the area corresponding to level V underwent a reduction during the 2000–2020 period.
The intensity analysis further discerned noteworthy trends, indicating rapid changes in DDI from 2000 to 2012, followed by a deceleration from 2012 to 2020. At the category level, the conversion intensity for level III, level IV, and level V reached its pinnacle from 2008 to 2012. In the transition level analysis, level III exhibited the highest intensity of both gain and loss.
The study emphasizes the influential role of meteorological factors and soil conditions in shaping the spatial distribution of DDI. Notably, evapotranspiration emerged as the dominant factor, with the highest q-value at 0.83, followed by precipitation, soil type, and land surface temperature. In contrast, temporal variations were predominantly attributable to human activities.
In conclusion, this research significantly contributes to advancing our understanding of the spatiotemporal dynamics of desertification in Xilingol. It offers valuable insights into the conversion patterns among different levels of desertification within Degraded Grassland (DGL), elucidating the distinct driving factors influencing the spatial distribution and temporal variation of desertification. The implications of these findings extend significantly to ecological research and conservation efforts in the region.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/rs15245716/s1, Figure S1: Proportions of DGL at different desertification levels; Table S1: Statistics of areas with different desertification levels (unit: km2).

Author Contributions

Conceptualization, J.L. and M.X.; data curation, X.Y.; formal analysis, J.L.; funding acquisition, M.X.; investigation, M.X., K.W. and Y.Y.; methodology, J.L. and X.Y.; project administration, C.C.; resources, X.G.; software, H.G.; supervision, C.C.; validation, X.G. and K.W.; visualization, H.G. and Y.Y.; writing—original draft, J.L.; writing—review and editing, C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the project of the National Key R&D Program of China (grant number 2021YFB3901104) and the National Natural Science Foundation of China (grant number 41971394).

Data Availability Statement

Data are contained within the article and supplementary materials.

Acknowledgments

We would like to thank the Resource and Environmental Science and Data Center, Chinese Academy of Sciences, for data support. We also extend our sincere gratitude to Cao and Xu for their guidance on the structure and language of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. UNCCD. United Nations Convention to Combat Desertification in Those Countries Experiencing Serious Drought and/or Desertification, Particularly in Africa; UN: Paris, France, 1994. [Google Scholar]
  2. Grainger, A. The Threatening Desert: Controlling Desertification; Routledge: London, UK, 2013. [Google Scholar]
  3. Feng, Q.; Tian, Y.; Yu, T.; Yin, Z.; Cao, S. Combating desertification through economic development in northwestern China. Land Degrad. Dev. 2019, 30, 910–917. [Google Scholar] [CrossRef]
  4. Liang, X.; Li, P.; Wang, J.; Shun Chan, F.K.; Togtokh, C.; Ochir, A.; Davaasuren, D. Research Progress of Desertification and Its Prevention in Mongolia. Sustainability 2021, 13, 6861. [Google Scholar] [CrossRef]
  5. National Forestry and Grassland Administration. Combating Desertification. Available online: https://www.forestry.gov.cn/main/6193/20220302/153916629457788.html (accessed on 6 October 2023).
  6. Tu, Z.; Li, M.; Sun, T. The Status and Trend Analysis of Desertification and Sandification. For. Resour. Manag. 2016, 1, 1–5+13. [Google Scholar] [CrossRef]
  7. Wang, Y. Building a Beautiful China Should Pay Attention to Grassland Ecological Governance and Protection. Available online: https://paper.cntheory.com/html/2020-10/07/nw.D110000xxsb_20201007_1-A7.htm (accessed on 7 October 2023).
  8. Reynolds, J.F.; Smith, D.M.S.; Lambin, E.F.; Turner, B.; Mortimore, M.; Batterbury, S.P.; Downing, T.E.; Dowlatabadi, H.; Fernández, R.J.; Herrick, J.E. Global desertification: Building a science for dryland development. Science 2007, 316, 847–851. [Google Scholar] [CrossRef] [PubMed]
  9. Kapfer, J.; Hédl, R.; Jurasinski, G.; Kopecký, M.; Schei, F.H.; Grytnes, J.A. Resurveying historical vegetation data–opportunities and challenges. Appl. Veg. Sci. 2017, 20, 164–171. [Google Scholar] [CrossRef]
  10. Wulder, M.A.; White, J.C.; Goward, S.N.; Masek, J.G.; Irons, J.R.; Herold, M.; Cohen, W.B.; Loveland, T.R.; Woodcock, C.E. Landsat continuity: Issues and opportunities for land cover monitoring. Remote Sens. Environ. 2008, 112, 955–969. [Google Scholar] [CrossRef]
  11. Rogan, J.; Chen, D. Remote sensing technology for mapping and monitoring land-cover and land-use change. Prog. Plan. 2004, 61, 301–325. [Google Scholar] [CrossRef]
  12. Feng, K.; Wang, T.; Liu, S.; Kang, W.; Chen, X.; Guo, Z.; Zhi, Y. Monitoring desertification using machine-learning techniques with multiple indicators derived from MODIS images in Mu Us Sandy Land, China. Remote Sens. 2022, 14, 2663. [Google Scholar] [CrossRef]
  13. Zeng, Y.; Xiang, N.; Feng, Z.; Xu, H. Albedo-NDVI Space and Remote Sensing Synthesis Index Models for Desertification Monitoring. Sci. Geogr. Sin. 2006, 26, 75–81. [Google Scholar] [CrossRef]
  14. Xie, Y.H.; Li, Z.H.; Zhu, L.; Zhou, Y.J.; Liu, H.Y.; Yu, T. Dynamic monitoring of desertification in response to climatic factors: A case study from the Gelintan Oasis on the southeastern edge of the Tengger Desert, China. Geocarto Int. 2022, 37, 9453–9473. [Google Scholar] [CrossRef]
  15. Wei, H.; Wang, J.; Han, B. Desertification Information Extraction Along the China–Mongolia Railway Supported by Multisource Feature Space and Geographical Zoning Modeling. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 392–402. [Google Scholar] [CrossRef]
  16. Wei, H.; Wang, J.; Cheng, K.; Li, G.; Ochir, A.; Davaasuren, D.; Chonokhuu, S. Desertification Information Extraction Based on Feature Space Combinations on the Mongolian Plateau. Remote Sens. 2018, 10, 1614. [Google Scholar] [CrossRef]
  17. Li, Y.S.; Li, H.K.; Xu, F. Spatiotemporal changes in desertified land in rare earth mining areas under different disturbance conditions. Environ. Sci. Pollut. Res. 2021, 28, 30323–30334. [Google Scholar] [CrossRef] [PubMed]
  18. Feng, Y.Y.; Wang, S.H.; Zhao, M.S.; Zhou, L.M. Monitoring of Land Desertification Changes in Urat Front Banner from 2010 to 2020 Based on Remote Sensing Data. Water 2022, 14, 1777. [Google Scholar] [CrossRef]
  19. Becerril-Pina, R.; Diaz-Delgado, C.; Mastachi-Loza, C.A.; Gonzalez-Sosa, E. Integration of remote sensing techniques for monitoring desertification in Mexico. Hum. Ecol. Risk Assess. 2016, 22, 1323–1340. [Google Scholar] [CrossRef]
  20. Wang, Z.; Liu, C.; Alfredo, H. From AVHRR-NDVI to MODIS-EVI: Advances in vegetation index research. Acta Ecol. Sin. 2003, 23, 979–987. [Google Scholar] [CrossRef]
  21. Shang, X.; He, Z.; Zhang, T. Spatial and temporal variations and the driving mechanism of the enhanced vegetation index. J. For. Environ. 2020, 40, 478–485. [Google Scholar] [CrossRef]
  22. Li, Y.; Tao, C.; Zhao, Q. A Remote Sensing Classification Method of Land Desertification Based on EVI-Albedo Space. Remote Sens. Inf. 2022, 37, 8–14. [Google Scholar]
  23. Zhao, X.; Zhou, D.; Fang, J. Satellite-based studies on large-scale vegetation changes in China. J. Integr. Plant Biol. 2012, 54, 713–728. [Google Scholar] [CrossRef]
  24. Feng, Q.; Ma, H.; Jiang, X.; Wang, X.; Cao, S. What Has Caused Desertification in China? Sci. Rep. 2015, 5, 15998. [Google Scholar] [CrossRef]
  25. Zhu, L.; Meng, J.; Zhu, L. Applying Geodetector to disentangle the contributions of natural and anthropogenic factors to NDVI variations in the middle reaches of the Heihe River Basin. Ecol. Indic. 2020, 117, 106545. [Google Scholar] [CrossRef]
  26. Na, R.; Du, H.; Na, L.; Shan, Y.; He, H.S.; Wu, Z.; Zong, S.; Yang, Y.; Huang, L. Spatiotemporal changes in the Aeolian desertification of Hulunbuir Grassland and its driving factors in China during 1980–2015. Catena 2019, 182, 104123. [Google Scholar] [CrossRef]
  27. Zhou, W.; Yang, H.; Huang, L.; Chen, C.; Lin, X.; Hu, Z.; Li, J. Grassland degradation remote sensing monitoring and driving factors quantitative assessment in China from 1982 to 2010. Ecol. Indic. 2017, 83, 303–313. [Google Scholar] [CrossRef]
  28. Han, J.; Wang, J.; Chen, L.; Xiang, J.; Ling, Z.; Li, Q.; Wang, E. Driving factors of desertification in Qaidam Basin, China: An 18-year analysis using the geographic detector model. Ecol. Indic. 2021, 124, 107404. [Google Scholar] [CrossRef]
  29. Wang, F.; Gou, R.; Qi, G.; Wang, S.; Li, B.; Guo, X. Effect of Climate Change on Desertification in Qaidamu Basin. Arid Meteorol. 2007, 25, 28–33. [Google Scholar] [CrossRef]
  30. Han, J.; Dai, H.; Gu, Z. Sandstorms and desertification in Mongolia, an example of future climate events: A review. Environ. Chem. Lett. 2021, 19, 4063–4073. [Google Scholar] [CrossRef] [PubMed]
  31. Jiang, L.; Jiapaer, G.; Bao, A.; Kurban, A.; Guo, H.; Zheng, G.; De Maeyer, P. Monitoring the long-term desertification process and assessing the relative roles of its drivers in Central Asia. Ecol. Indic. 2019, 104, 195–208. [Google Scholar] [CrossRef]
  32. Wang, X.; Chen, F.; Hasi, E.; Li, J. Desertification in China: An assessment. Earth-Sci. Rev. 2008, 88, 188–206. [Google Scholar] [CrossRef]
  33. Ibáñez, J.; Martínez, J.; Schnabel, S. Desertification due to overgrazing in a dynamic commercial livestock–grass–soil system. Ecol. Model. 2007, 205, 277–288. [Google Scholar] [CrossRef]
  34. Hilker, T.; Natsagdorj, E.; Waring, R.H.; Lyapustin, A.; Wang, Y. Satellite observed widespread decline in Mongolian grasslands largely due to overgrazing. Glob. Chang. Biol. 2014, 20, 418–428. [Google Scholar] [CrossRef]
  35. D’Odorico, P.; Bhattachan, A.; Davis, K.F.; Ravi, S.; Runyan, C.W. Global desertification: Drivers and feedbacks. Adv. Water Resour. 2013, 51, 326–344. [Google Scholar] [CrossRef]
  36. Freier, K.P.; Schneider, U.A.; Finckh, M. Dynamic interactions between vegetation and land use in semi-arid Morocco: Using a Markov process for modeling rangelands under climate change. Agric. Ecosyst. Environ. 2011, 140, 462–472. [Google Scholar] [CrossRef]
  37. Vieira, R.M.D.; Tomasella, J.; Barbosa, A.A.; Martins, M.A.; Rodriguez, D.A.; Rezende, F.S.; Carriello, F.; Santana, M.D. Desertification risk assessment in Northeast Brazil: Current trends and future scenarios. Land Degrad. Dev. 2021, 32, 224–240. [Google Scholar] [CrossRef]
  38. Xu, D.; Song, A.; Tong, H.; Ren, H.; Hu, Y.; Shao, Q. A spatial system dynamic model for regional desertification simulation–A case study of Ordos, China. Environ. Model. Softw. 2016, 83, 179–192. [Google Scholar] [CrossRef]
  39. Yang, Z.W.; Gao, X.; Lei, J.Q.; Meng, X.Y.; Zhou, N. Analysis of spatiotemporal changes and driving factors of desertification in the Africa Sahel. Catena 2022, 213, 106213. [Google Scholar] [CrossRef]
  40. Wang, Y.; Guo, E.; Kang, Y.; Ma, H. Assessment of Land Desertification and Its Drivers on the Mongolian Plateau Using Intensity Analysis and the Geographical Detector Technique. Remote Sens. 2022, 14, 6365. [Google Scholar] [CrossRef]
  41. Qi, K.; Zhu, J.J.; Zheng, X.; Wang, G.G.; Li, M.C. Impacts of the world’s largest afforestation program (Three-North Afforestation Program) on desertification control in sandy land of China. GISci. Remote Sens. 2023, 60, 2167574. [Google Scholar] [CrossRef]
  42. Aldwaik, S.Z.; Pontius, R.G. Intensity analysis to unify measurements of size and stationarity of land changes by interval, category, and transition. Landsc. Urban Plan. 2012, 106, 103–114. [Google Scholar] [CrossRef]
  43. Qian, K.; Yuan, Q.Z.; Han, J.C.; Leng, R.; Wang, Y.S.; Zhu, K.H.; Lin, S.; Ren, P. A remote sensing monitoring method for alpine grasslands desertification in the eastern Qinghai-Tibetan Plateau. J. Mt. Sci. 2020, 17, 1423–1437. [Google Scholar] [CrossRef]
  44. Xu, M.R.; Zhang, J.; Li, Z.H.; Mo, Y. Attribution analysis and multi-scenario prediction of NDVI drivers in the Xilin Gol grassland, China. J. Arid Land 2022, 14, 941–961. [Google Scholar] [CrossRef]
  45. Li, H.; Yang, X. Temperate dryland vegetation changes under a warming climate and strong human intervention—With a particular reference to the district Xilin Gol, Inner Mongolia, China. Catena 2014, 119, 9–20. [Google Scholar] [CrossRef]
  46. Yang, X.; Xu, B.; Zhu, X.; Jin, Y.; Li, J.; Zhao, F.; Chen, S.; Guo, J.; Ma, H.; Yu, H. A monitoring indicator system for remote sensing of grassland vegetation growth and suitability evaluation—A case study of the Xilingol Grassland in Inner Mongolia, China. Int. J. Remote Sens. 2015, 36, 5105–5122. [Google Scholar] [CrossRef]
  47. Nemati Nasab, M.R.; Sattari Sarbangholi, H.; Pakdelfard, M.R.; Jamali, S. Structural equation modeling (SEM) the relationship between environmental quality and social cohesion components by explaining the mediating role of social resilience in urban cultural spaces. Int. J. Nonlinear Anal. Appl. 2022, 14, 153–167. [Google Scholar]
  48. Wang, K.; Delbaere, K.; Brodie, M.A.; Lovell, N.H.; Kark, L.; Lord, S.R.; Redmond, S.J. Differences between gait on stairs and flat surfaces in relation to fall risk and future falls. IEEE J. Biomed. Health Inform. 2017, 21, 1479–1486. [Google Scholar] [CrossRef] [PubMed]
  49. Huang, B.; Huang, J.; Pontius, R.G., Jr.; Tu, Z. Comparison of Intensity Analysis and the land use dynamic degrees to measure land changes outside versus inside the coastal zone of Longhai, China. Ecol. Indic. 2018, 89, 336–347. [Google Scholar] [CrossRef]
  50. Du, Z.; Yu, L.; Chen, X.; Gao, B.; Yang, J.; Fu, H.; Gong, P. Land use/cover and land degradation across the Eurasian steppe: Dynamics, patterns and driving factors. Sci. Total Environ. 2023, 909, 168593. [Google Scholar] [CrossRef]
  51. André, M. The ARTEMIS European driving cycles for measuring car pollutant emissions. Sci. Total Environ. 2004, 334, 73–84. [Google Scholar] [CrossRef]
  52. Liu, G.; Liu, H.; Yin, Y. Global patterns of NDVI-indicated vegetation extremes and their sensitivity to climate extremes. Environ. Res. Lett. 2013, 8, 025009. [Google Scholar] [CrossRef]
  53. Piedallu, C.; Cheret, V.; Denux, J.-P.; Perez, V.; Azcona, J.S.; Seynave, I.; Gégout, J.-C. Soil and climate differently impact NDVI patterns according to the season and the stand type. Sci. Total Environ. 2019, 651, 2874–2885. [Google Scholar] [CrossRef]
  54. Batunacun; Nendel, C.; Hu, Y.; Lakes, T. Land-use change and land degradation on the Mongolian Plateau from 1975 to 2015-A case study from Xilingol, China. Land Degrad. Dev. 2018, 29, 1595–1606. [Google Scholar] [CrossRef]
  55. Li, S.; Xie, Y. Investigating Coupled Impacts of Climate Change and Socioeconomic Transformation on Desertification by Using Multitemporal Landsat Images: A Case Study in Central Xilingol, China. IEEE Geosci. Remote Sens. Lett. 2013, 10, 1244–1248. [Google Scholar] [CrossRef]
  56. Zhang, Y.; Ya, X.; Wang, R.; Zou, Y.; Dong, X. Assessing and mapping human well-being for sustainable development amid drought and flood hazards: Dadu River Basin of China. Environ. Sci. Pollut. Res. 2022, 29, 90719–90737. [Google Scholar] [CrossRef] [PubMed]
  57. Saad, A.M.A.; Shariff, N.M.; Gairola, S. Nature and causes of land degradation and desertification in Libya: Need for sustainable land management. Afr. J. Biotechnol. 2011, 10, 13680–13687. [Google Scholar]
  58. Huang, J.; Li, Y.; Fu, C.; Chen, F.; Fu, Q.; Dai, A.; Shinoda, M.; Ma, Z.; Guo, W.; Li, Z. Dryland climate change: Recent progress and challenges. Rev. Geophys. 2017, 55, 719–778. [Google Scholar] [CrossRef]
  59. Fu, G.; Qiu, X.; Xu, X.; Zhang, W.; Zang, F.; Zhao, C. The role of biochar particle size and application rate in promoting the hydraulic and physical properties of sandy desert soil. Catena 2021, 207, 105607. [Google Scholar] [CrossRef]
  60. Kalma, J.D.; McVicar, T.R.; McCabe, M.F. Estimating land surface evaporation: A review of methods using remotely sensed surface temperature data. Surv. Geophys. 2008, 29, 421–469. [Google Scholar] [CrossRef]
  61. National Forestry and Grassland Administration. Improve the Grassland Governance System and Promote the Construction of Ecological Civilization. Available online: http://www.forestry.gov.cn/cys/3/20210816/163246652498542.html (accessed on 6 October 2023).
  62. Sarker, S. Fundamentals of Climatology for Engineers: Lecture Note. Eng 2022, 3, 573–595. [Google Scholar] [CrossRef]
Figure 1. Location of study area. The upper-left map displays the location of Xilingol within China, while the main map illustrates the topography and administrative subdivisions of Xilingol.
Figure 1. Location of study area. The upper-left map displays the location of Xilingol within China, while the main map illustrates the topography and administrative subdivisions of Xilingol.
Remotesensing 15 05716 g001
Figure 2. Workflow of the study. The study consisted of four parts: input data, data preprocessing, model construction, and analysis.
Figure 2. Workflow of the study. The study consisted of four parts: input data, data preprocessing, model construction, and analysis.
Remotesensing 15 05716 g002
Figure 3. Linear fitting of albedo-EVI (the figure shows the result from 2020).
Figure 3. Linear fitting of albedo-EVI (the figure shows the result from 2020).
Remotesensing 15 05716 g003
Figure 4. Spatiotemporal distribution of grassland desertification degree in Xilingol.
Figure 4. Spatiotemporal distribution of grassland desertification degree in Xilingol.
Remotesensing 15 05716 g004
Figure 5. Spatial distribution of vegetation in Xilingol.
Figure 5. Spatial distribution of vegetation in Xilingol.
Remotesensing 15 05716 g005
Figure 6. Linear fitting for the area variation of different levels of DGL.
Figure 6. Linear fitting for the area variation of different levels of DGL.
Remotesensing 15 05716 g006
Figure 7. Interval level of intensity analysis.
Figure 7. Interval level of intensity analysis.
Remotesensing 15 05716 g007
Figure 8. Category level of intensity analysis for (A) 2000–2004, (B) 2004–2008, (C) 2008–2012, (D) 2012–2016, and (E) 2016–2020).
Figure 8. Category level of intensity analysis for (A) 2000–2004, (B) 2004–2008, (C) 2008–2012, (D) 2012–2016, and (E) 2016–2020).
Remotesensing 15 05716 g008
Figure 9. The sources of increased area and destinations of decreased area for different levels of DGL during all periods ((A) level I, (B) level II, (C) level III, (D) level IV (E), and level V).
Figure 9. The sources of increased area and destinations of decreased area for different levels of DGL during all periods ((A) level I, (B) level II, (C) level III, (D) level IV (E), and level V).
Remotesensing 15 05716 g009
Figure 10. Transition level of intensity analysis ((A) level I, (B) level II, (C) level III, (D) level IV, and (E) level V).
Figure 10. Transition level of intensity analysis ((A) level I, (B) level II, (C) level III, (D) level IV, and (E) level V).
Remotesensing 15 05716 g010
Figure 11. Input factors of the Geo-detector model. In the figure, (A): LST (Land Surface Temperature), (B): PRE (Precipitation), (C): WV (Wind Velocity), (D): ET (Evapotranspiration), (E): POP (Population), (F): SET (Soil Erosion Type), (G): SEI (Soil Erosion Intensity), (H): SM (Soil Moisture), (I): DEM (Elevation), (J): SLP (Slope), (K): ASP (Aspect), (L): LIV (Livestock), (M): ST (Soil Type), (N): LUCC (Land Use and Cover Change).
Figure 11. Input factors of the Geo-detector model. In the figure, (A): LST (Land Surface Temperature), (B): PRE (Precipitation), (C): WV (Wind Velocity), (D): ET (Evapotranspiration), (E): POP (Population), (F): SET (Soil Erosion Type), (G): SEI (Soil Erosion Intensity), (H): SM (Soil Moisture), (I): DEM (Elevation), (J): SLP (Slope), (K): ASP (Aspect), (L): LIV (Livestock), (M): ST (Soil Type), (N): LUCC (Land Use and Cover Change).
Remotesensing 15 05716 g011
Figure 12. Results of factor detector.
Figure 12. Results of factor detector.
Remotesensing 15 05716 g012
Figure 13. Results of ecological detector.
Figure 13. Results of ecological detector.
Remotesensing 15 05716 g013
Figure 14. Results of interaction detector.
Figure 14. Results of interaction detector.
Remotesensing 15 05716 g014
Figure 15. Temporal variation of factors from 2000 to 2020 (the black dashed line is the result of linear fitting).
Figure 15. Temporal variation of factors from 2000 to 2020 (the black dashed line is the result of linear fitting).
Remotesensing 15 05716 g015
Table 1. Data description and sources.
Table 1. Data description and sources.
CategoryNameCodeSourceSpatial ResolutionTemporal ResolutionUnit
Vegetation indexEnhanced Vegetation IndexEVIhttps://lpdaac.usgs.gov/ (accessed on 4 September 2023)250 m16 day/
Meteorological factorsLand surface temperatureLSThttps://lpdaac.usgs.gov (accessed on 4 September 2023)1000 mPer day°C
PrecipitationPREhttps://chc.ucsb.edu/data0.05°Per daymm
Wind velocityWVhttps://www.climatologylab.org/terraclimate.html (accessed on 4 September 2023)4638.3 mPer monthm/s
EvapotranspirationEThttps://lpdaac.usgs.gov/ (accessed on 4 September 2023)500 m8 daykg/m2
Albedo/https://lpdaac.usgs.gov/ (accessed on 4 September 2023)500 mPer day/
Soil factorsSoil typeSThttp://www.resdc.cn (accessed on 5 September 2023)shapefile//
Soil erosion typeSEThttp://www.resdc.cn (accessed on 6 September 2023)shapefile//
Soil erosion intensitySEIhttp://www.resdc.cn (accessed on 5 September 2023)shapefile//
Soil moistureSMhttps://disc.gsfc.nasa.gov (accessed on 5 September 2023)27,830 m/kg/m2
Topographical factorsDigital elevation modelDEMhttps://cmr.earthdata.nasa.gov/ (accessed on 4 September 2023)30 m/m
SlopeSLPCalculated from DEM30 m/°
AspectASPCalculated from DEM30 m//
Human activitiesPopulationPOPhttps://sedac.ciesin.columbia.edu/ (accessed on 10 September 2023)1000 mPer yearpersons/km2
LivestockLIVhttps://www.fao.org/ (accessed on 10 September 2023)10,000 m/heads/km2
Land use/cover changeLUCChttps://lpdaac.usgs.gov/ (accessed on 10 September 2023)500 mPer year/
StatisticsPopulation/http://tjj.xlgl.gov.cn/ (accessed on 11 September 2023)/Per yeark
GDP/http://tjj.xlgl.gov.cn/ (accessed on 11 September 2023)/Per yearbillion RMB
Livestock/http://tjj.xlgl.gov.cn/ (accessed on 12 September 2023)/Per yearmillion
Table 2. Classification standards for desertification level.
Table 2. Classification standards for desertification level.
CategoryValue of DDI
I (Severely Desertified)DDI < −0.45
II (Moderately Desertified)−0.45 ≤ DDI < −0.25
III (Slightly Desertified)−0.25 ≤ DDI < −0.05
IV (Nondesertified)−0.05 ≤ DDI < 0.15
V (Healthy)DDI ≥ 0.15
Table 3. Definition of interaction detector.
Table 3. Definition of interaction detector.
DescriptionInteraction
P(x1x2) < min(P(x1), P(x2))Nonlinear weaken
min(P(x1), P(x2)) < P(x1x2) < max(P(x1), P(x2))Uni-weaken
P(x1x2) > max(P(x1), P(x2)) and P(x1x2) < P(x1) + P(x2)Bi-enhance
P(x1x2) > P(x1) + P(x2)Nonlinearly enhance
P(x1x2) = P(x1) + P(x2)Independent
Table 4. Results of correlation analysis.
Table 4. Results of correlation analysis.
Human ActivitiesMeteorological Factors
ρ 0.5360.111
Significance0.0120.632
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, J.; Cao, C.; Xu, M.; Yang, X.; Gao, X.; Wang, K.; Guo, H.; Yang, Y. A 20-Year Analysis of the Dynamics and Driving Factors of Grassland Desertification in Xilingol, China. Remote Sens. 2023, 15, 5716. https://doi.org/10.3390/rs15245716

AMA Style

Li J, Cao C, Xu M, Yang X, Gao X, Wang K, Guo H, Yang Y. A 20-Year Analysis of the Dynamics and Driving Factors of Grassland Desertification in Xilingol, China. Remote Sensing. 2023; 15(24):5716. https://doi.org/10.3390/rs15245716

Chicago/Turabian Style

Li, Jingbo, Chunxiang Cao, Min Xu, Xinwei Yang, Xiaotong Gao, Kaimin Wang, Heyi Guo, and Yujie Yang. 2023. "A 20-Year Analysis of the Dynamics and Driving Factors of Grassland Desertification in Xilingol, China" Remote Sensing 15, no. 24: 5716. https://doi.org/10.3390/rs15245716

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