Assessment of the Effectiveness of Sand-Control and Desertiﬁcation in the Mu Us Desert, China

: The ﬁrst successful sand-control was achieved in the Mu Us Desert by local people in the 1950–1960s, and their experience and approach have been extended to the whole Ordos and Northern China since then. The objective of this paper is to assess comprehensively the effectiveness of sand-control in 15 counties in and around Mu Us using multitemporal satellite images and socioeconomic data. After atmospheric correction, Landsat TM and OLI images were harnessed for land cover classiﬁcation based on the ground-truth data and for derivation of the GDVI (generalized difference vegetation index) to extract the biophysical changes of the managed desert and desertiﬁcation. Climatic, socioeconomic, environmental and spatial factors were selected for coupling analysis by multiple linear and logistic regression models to reveal the driving forces of desertiﬁcation and their spatial determinants. The results show that from 1991 to 2020, 8712 km 2 or 63% of the desert has been converted into pastures and shrublands with a greenness increase of 0.3509 in GDVI; the effectiveness of sand-control is favored by the rational agropastoral activities and policies; though desertiﬁcation occurs locally, it is associated with both climatic and socioeconomic factors, such as wind water availability, distance to animal husbandry.


Introduction
Desertification is a phenomenon of land degradation resulting from by human activities and climate variations in arid, semi-arid and partially sub-humid areas [1,2]. With more than about four million km 2 of dryland, China has been suffering from desertification since centuries ago. Though sand-control through afforestation and plantation of shrubs and herbaceous vegetation has been widely effectuated in the dryland areas since the 1950s, there is still large area of desert in China, reaching 739,200 km 2 [3]. This poses a serious threat to the sustainable socioeconomic development in the arid region.
Remote sensing (RS) technology has become an important data source and technical means for desertification monitoring due to its wide detection range and observation periodicity since the 1970s [2,[4][5][6][7][8][9][10]. It has been successfully applied to distinguish desert encroachment or climate events of the south border movement of the Saharan Desert [5], desertification in Sahel [4], in western Rajasthan in India in 1990s [11], and in Mu Us [2,8]. Remote sensing approaches include, but are not limited to, differencing-and-thresholding technique and supervised classification [7,8,[12][13][14], in which the quantitative evaluation of desertification based on remote sensing indicators is made possible [15,16]. Commonlyused indicators are the normalized difference vegetation index (NDVI), fractional vege- The average annual rainfall is about 317 mm, mainly concentrated in July-September [2,8], while annual evaporation reaches 2500 mm, much stronger than precipitation. Winds from the northwest blow on average 230 days per year, and those exceeding the speed of 17 m/s (Beaufort Scale-8) may occur during more than 40 days [8]. The mean annual temperature is about 6.0-8.5 °C, but there is a significant difference between summer and winter and between day and night.
The vegetation type is mainly shrubs and herbs in the sandy land, and rice and maize are the main cultivation in the Yinchuan Plain with which the greenness reaches its maximum in July and August. The special terrain and climate conditions, as well as the overgrazing in the early periods, have made the study area largely desertified or susceptible to desertification.
There are rich mineral resources including coal, oil and gas in the study area, and that is why the latter was designated as one of the energy bases of China by the central government in 1999 [8]. The average annual rainfall is about 317 mm, mainly concentrated in July-September [2,8], while annual evaporation reaches 2500 mm, much stronger than precipitation. Winds from the northwest blow on average 230 days per year, and those exceeding the speed of 17 m/s (Beaufort Scale-8) may occur during more than 40 days [8]. The mean annual temperature is about 6.0-8.5 • C, but there is a significant difference between summer and winter and between day and night.
The vegetation type is mainly shrubs and herbs in the sandy land, and rice and maize are the main cultivation in the Yinchuan Plain with which the greenness reaches its maximum in July and August. The special terrain and climate conditions, as well as the overgrazing in the early periods, have made the study area largely desertified or susceptible to desertification.
There are rich mineral resources including coal, oil and gas in the study area, and that is why the latter was designated as one of the energy bases of China by the central government in 1999 [8].
Since 1959, local people started to conduct sand-control activities and many largescale combating desertification campaigns have been successfully undertaken by individual/household companies and national teams since then.

Data and Processing
To achieve our purposes, a combined method of post-classification differencing and GDVI differencing was proposed to quantitatively evaluate the spatial distribution and degree of sand-control. A set of socioeconomic, climatic and environmental factors were selected to quantitatively analyze the driving forces of the spatial heterogeneity of desert governance in different periods using multiple linear regression and logistic regression models.

Data Satellite Data
Landsat TM, ETM+ and OLI images from 1987 to 2020 with path/row numbers 128/33-34 and 129/33-34 and 16 tiles of digital elevation model (DEM) data (ASTERV003, 30 m resolution) of the study area were obtained from the USGS data server (glovis.usgs.gov, accessed on 16 June 2020) and NASA (earthdata.nasa.gov, accessed on 16 June 2020). Landsat images were mostly acquired in summer from July to September with cloud cover less than 1% ( Table 1) in terms of the phenological features of crops and vegetation in the continental climate zone. Due to frequent cloudy weather in summer, it is difficult to acquire cloud-free images for the adjacent scenes of different paths in the same year. Hence, selection of images depended on the availability of cloud-free images and there was difference in acquisition dates of the adjacent scenes of images. The obtained Landsat images were utilized for the following supervised classification and derivation of different vegetation index, GDVI, and land surface temperature (LST). High resolution images available on Google Earth were harnessed to derive road vector data for 1991,1999, 2010 and 2020.

Socioeconomic Data
It is difficult to monitor and quantify human activity, and the best way is to use annual socioeconomic data, the finally quantified expression of human activity during a one-year period, for research [2]. For the study area, county-level socioeconomic data are available from the Statistic Yearbooks published by the government. The data include the total sown area, meat product, total number of livestock at the end of year, sheep stock, per capita net income of farmers and herdsmen, gross output of agriculture, forestry, animal husbandry and fishery, gross domestic product (GDP) and total population. For this reason, the Statistic Yearbooks of Inner Mongolia, Ningxia and Shaanxi from 1990 to 2020 were collected. Due to the change of administrative divisions, we combined the data of Xixia, Xingqing and Jinfeng districts into Yinchuan City; Wuda, Hainan and Haibowan districts were merged into Wuhai City; and Dawukou and Huinong districts were incorporated into Shizuishan City. Finally, socioeconomic data of 15 counties and cities were prepared for analysis.

Meteorological Data
The monthly average temperature, average precipitation, average wind speed, maximum wind speed and sunshine duration of four stations, e.g., Yulin, Etog, Otogqian and Uxin, were obtained from the website of National  (4) sandcontrol achievement by local and national companies in the whole study area. In total, more than 2000 observation points were collected with GPS location. One surprising thing observed is that sandy land has greatly decreased from 2000 to 2021 and converted mostly into shrublands and forests in dunes and grasslands in interdunes thanks to the sand-control activity.

Processing Procedures Satellite Data Processing
Atmospheric correction of Landsat images was conducted using the COST model developed by Chavez (1996) [40], where the band minimum was employed to remove the haze effect and at-satellite radiance of all bands was converted into surface reflectance.
Land surface temperature (LST) was derived from the thermal band of Landsat TM and OLI data [41,42].
The dryland-tailored vegetation index, GDVI, is shown in Equation (1). When n = 1, GDVI = NDVI; when n = 2, it is suitable for characterizing dryland biomes including shrubland, woodland and forest; when n = 3, it can be used for monitoring the degradation and desertification in sparsely vegetated area [22,23]. Compared with NDVI and other vegetation indices, GDVI has higher sensitivity and dynamic range, and is able to identify more effectively subtle differences in vegetation greenness in low vegetated areas [22]. Therefore, GDVI of n = 3 was produced for the successive analysis.
where, ρ NIR is the spectral reflectance of the near-infrared band, and ρ R is that of the red band. Slope and aspect were produced from the DEM data and used for land cover classification [24] and desertification modeling as a part of the environmental factors.

Land Cover Classification
A dataset of eleven bands, including elevation, slope, aspect, LST, GDVI and six spectral bands (blue, green, red, near infrared, shortwave infrared 1 (SWIR 1) and SWIR 2) through layer-stacking function was composed for each observation year [24].
Based on field survey and knowledge, we defined 30-32 initial classes of groundtruth samples (regions of interest, ROIs) for different scenes and observation years, and half of the samples (about 5-6% of the total scene) were used for training, i.e., training samples (TS), and the remaining half for validation, or rather, validation samples (VS). Sand-control and desertified land (or deserts) were the important classes to be classified out. After classification, all initial classes were merged into 12 classes in terms of the similarity, Remote Sens. 2022, 14, 837 6 of 24 namely sandy land, shrublands, grasslands, waters, artificial (built up) areas, bare soil, loess land, woodlands, coal mines, saline land, croplands and fallows. Forests were sporadic and less than several pixels, and merged into woodlands. Finally, the sandy land was extracted as a mask, and all the remaining classes were merged into non-sand land for further analysis.
Supervised classification using maximum likelihood approach was conducted, and the results were validated with VS.

Processing of the Driving Factors
The distance from the random points in the managed areas and desertified land to roads, water bodies, residential areas and farmland was obtained by proximity analysis. The change of water bodies and mining area were obtained by the statistics of classified images. Population density was calculated by Equation (2).
where, Y is the population density, X is total population, and A is the area.

Post-Classification Differencing
Post-classification differencing is a method for identification of land cover change proposed by Wu [12]. The same type of land cover of the different observation time, deserts and controlled sandy land in this case, was extracted from the classified images and a differencing was applied to obtain the change in space. Equation (3) was used to calculate the difference. Taking desertification as an example, when ∆S = 0, it means no change in sandy land during this period; when ∆S = 1, the desert or sandy land has expanded; when ∆S = −1, the sandy land is reduced and controlled [8,12]. This method can be used to monitor the overall change of sandy land and identify accurately the area and spatial distribution of sandy land, and at the same time, quantify which land cover types have been converted into sandy land or vice versa. This will allow to produce the gain/loss change matrix of land cover of the study area in different observation periods.

∆GDVI
Based on the calculation of GDVI of different time images, a differencing was effectuated under the masks of deserts or sand-control extracted from the classified images: When ∆GDVI < 0, desertification has occurred as vegetation vigor has decreased; if ∆GDVI > 0, indicating a vegetation increase, or rather, sandy land has been controlled and converted into vegetated area such as grasslands and shrublands.

Multiple Linear Regression Analysis
Multiple linear stepwise regression is an approach to analyze the relationship between a dependent variable and several independent variables, and to find which independent variable(s) contribute(s) most to the dependent or which factor(s) has(have) played the most important role in the event occurrence (Equation (5)) [2,6,43].
where, Y i is the dependent variable, X 1i . . . X ki are the independent variables, and β 0 , β 1 . . . , β k are the regression coefficients of their corresponding variables, and U i is the random error. The dependent variable in this case is the desertified land or sand-control area at county-level, and the independent variables are county-level socioeconomic and climatic factors.

Logistic Regression Analysis
Logistic regression model (LRM) is capable of dealing with both categorical and continuous variables and can effectively achieve spatially explicit analysis to reveal the probability of certain change in space (i.e., with probability of 1 indicating event occurred, and 0 indicating it has not occurred) associated with multiple independent variables, i.e., spatial determinants. It was first used in the field of disease diagnosis and has been extended to land change research to reveal the spatial determinants of land cover change in environmental geography and geological disaster prediction in recent decades [2,[44][45][46]. However, it has been rarely used in the quantitative analysis of desertification driving force because of the particularity of its dependent variable. The formula (6) was used to establish the LRM.
where, P is the probability, X 1 . . . , X i are the independent variables, and β 0 , β 1 . . . , β i are the regression coefficients of their corresponding variables. Before logistic regression analysis, it is necessary to undertake the collinearity diagnosis to check whether there is collinearity among the independent variables using either tolerance (TOL) or variable inflation factor (VIF) as an indicator. If VIF < 10 or TOL > 0.1, there is no evident collinearity and the result of the logistic regression analysis shall be reliable.
In order to examine the importance of factors influencing sand-control and desertification, the following preprocessing is required. First, rasterization of the independent variables, e.g., GDP, per capita rural net income, meet product and so on are converted into county-level raster; then, creation of the random points in the controlled and uncontrolled areas of sandy lands, the desertified areas and the unchanged areas in line with the ratio of 1:1. The points in the desertified areas and sand-control areas were assigned a value of 1, and in non-desertification and non-control areas were assigned a value of 0, which were used as dependent variables. At last, the values of independent and dependent variables were both extracted and input into SPSS, a package for statistic analysis, for logistic regression modeling.
The model was tested in two aspects: (1) the Hosmer-Lemeshow (H-L) test was used to evaluate the degree of fit of the model; if H-L > 0.05, the model was considered to have a good fit; and (2) using significance (Sig), e.g., Sig < 0.05, to test whether the independent variable in the model has a significant effect on the dependent variable. Finally, the absolute values of the regression coefficients or odds ratios (OR) were used for judging the influence degree of the independent variables on the dependent variable.

Dynamic Situation of the Desert
As shown in Table 2, the overall accuracy (OA) of the land cover classification vs. the validation samples is about 90.3-92.8% with a kappa coefficient of 0.87-0.90. The desert or sandy land is mainly distributed in the banners Otog, Otogqian and Uxin, among which the last has the largest area (Table 3 and Figure 2). In the past 30 years, sandy land has decreased by 8712.23 km 2 , accounting for 63.05% of the total sandy land though there were slight fluctuation in some counties in different periods, e.g., from 1991 to 1999, sandy land increased by 291.6 km 2 in Yanchi, and from 1999 to 2010, sandy land in Yinchuan, Otogqian and Wuzhong had an increase of 32.22 km 2 , 65.13 km 2 and 3.59 km 2 , respectively, while a decrease appeared in other areas in the whole period.

Effectiveness of the Combating Desertification
Between 1991-2020, a total of 9140.44 km 2 of sandy land has been managed and converted into pastures, shrublands, and even forests ( Table 4). The main sand-control areas are located in the banners (counties) Uxin, Otog and Otogqian, with an area of 3022.76 km 2 , 2090.48 km 2 and 1637.36 km 2 respectively. As shown in Figures 3 and 4, the area of sand-control was much larger than the desertified land.
Between 1991-2020, a total of 9140.44 km 2 of sandy land has been managed and converted into pastures, shrublands, and even forests ( Table 4). The main sand-control areas are located in the banners (counties) Uxin, Otog and Otogqian, with an area of 3022.76 km 2 , 2090.48 km 2 and 1637.36 km 2 respectively. As shown in Figures 3 and 4, the area of sand-control was much larger than the desertified land.  Most of the sand-control was effectuated in the western margin of the desert, i.e., the west part of the banners Otog and Otogqian, followed by the southern area in Dingbian, Lingwu and Yanchi (Figures 3 and 4). The desertification occurred mainly in Yanchi and Lingwu in 1991-1999 (Figure 4a), but it was concentrated in the west of Otogqian and Otog in 1999-2010 (Figure 4b), and it appeared as small dots in the east of the study area    As seen in Table 4, the vegetation greenness represented by ΔGDVI in the con  As seen in Table 4, the vegetation greenness represented by ΔGDVI in the contro areas in the past 30 years has gained an increase of 0.3518 with a minimum of 0.145 2010-2020, while the greenness had decreased by 0.0691 in the desertified area, and Most of the sand-control was effectuated in the western margin of the desert, i.e., the west part of the banners Otog and Otogqian, followed by the southern area in Dingbian, Lingwu and Yanchi (Figures 3 and 4). The desertification occurred mainly in Yanchi and Lingwu in 1991-1999 (Figure 4a), but it was concentrated in the west of Otogqian and Otog in 1999-2010 (Figure 4b), and it appeared as small dots in the east of the study area in 2010-2020, and a large piece of desertification was not observed in the last period ( Figure 4c).
As seen in Table 4, the vegetation greenness represented by ∆GDVI in the controlled areas in the past 30 years has gained an increase of 0.3518 with a minimum of 0.1456 in 2010-2020, while the greenness had decreased by 0.0691 in the desertified area, and the biggest decrease (−0.1237) appeared in 1999-2010.
Presented in Tables 5-8, sandy land was mostly transferred into shrublands and grasslands with an area of 6359.79 km 2 , accounting for 69.58%, from 1991 to 2020. This conversion took up 4943.51 km 2 , 4088.33 km 2 and 2606.92 km 2 , accounting for 89.27%, 82.27% and 76.22%, respectively, in the three observed periods 1991-2000, 2000-2010 and 2010-2020. It is obvious that sand-control by planting shrubs such as Salix cheilophila, Tamarix ramosissima, Hedysarum mongolicum, Caragana korshinskii Kom, Hippophae rhamnoides L. and Lycium chinense Miller and herbaceous vegetation including Medicago sativa L., Artemisia desertorum, Artemisia sieversiana, Astragalus adsurgens Pall. has achieved a remarkable output.  show that the extent of sand-control was much larger than that of desertification, whereas it was the opposite in Yanchi during 1991-1999. From 1991 to 2020, sand-control seems to have had the best effectiveness in Uxin with an increase of 0.0897 in GDVI, and the lowest effectiveness in Shizuishan, with an increase of 0.0002. Meanwhile, ∆GDVI of the desertified area was close to 0 ( Figure 5).  5 show that the extent of sand-control was much larger than that of desertification, whereas it was the opposite in Yanchi during 1991-1999. From 1991 to 2020, sand-control seems to have had the best effectiveness in Uxin with an increase of 0.0897 in GDVI, and the lowest effectiveness in Shizuishan, with an increase of 0.0002. Meanwhile, ΔGDVI of the desertified area was close to 0 ( Figure 5).

Determinants of Sand-Control and Desertification
As shown in Table 9, 15 factors from three aspects were used for stepwise linear regression modeling with the area of sand-control and desertification as the dependent variable. Taking R 2 > 0.60 as the model testing standard, seven models were finally obtained (Table 10).

Determinants of Sand-Control and Desertification
As shown in Table 9, 15 factors from three aspects were used for stepwise linear regression modeling with the area of sand-control and desertification as the dependent variable. Taking R 2 > 0.60 as the model testing standard, seven models were finally obtained (Table 10).
Note: Here Y 1 , Y 2 , Y 3 and Y 4 are the areas of sand-control in each county in 1991-1999, 1999-2010, 2010-2020 and 1991-2020, respectively, while Y 5 , Y 6 , and Y 7 are the areas of desertification in each county in 1991-1999, 2010-2020 and 1991-2020, respectively. Of the corresponding period, ∆X 10 is the difference in annual mean temperature, ∆X 15 the difference in coal mine area, ∆X 2 the difference in meat product, ∆X 13 the difference in sunshine duration, ∆X 5 the difference in per capita net income of farmers and herdsmen, ∆X 3 the difference in sheep number, ∆X 9 the difference in annual precipitation, ∆X 12 the difference in average wind speed, and ∆X 14 the difference in area of water body.
In terms of the multiple linear regression models in Table 10, it is seen that sandcontrol area is closely associated with the difference in temperature, mining area, meat product, sunshine hours, per capita net income of farmers and herdsmen, sheep number, and precipitation between the two observation years. More concretely, the increases in mining area, temperature and sunshine hours are negatively correlated with the sandcontrolled area, while the increase in meat (pork, beef and mutton) product, per capita net income of farmers and herdsmen, sheep number and precipitation are positively correlated with the sand-controlled area.
The desertified area is related to the average wind speed and water. More exactly, the increase in average wind speed is the main driving force of desertification, while the increase in water areas are the favorable factors to reverse desertification, that is, lead to a decrease in desertified area.

Spatial Variability of Sand-Control
Different from those used for multiple linear regression analysis, here 19 socioenvironmental factors shown in Table 11 were used for logistic regression modeling to understand the different spatial distribution of sand-control and desertification in different periods, and the obtained models that have passed the Hosmer-Lemeshow (HL) test are presented in Table 12. It is worth mentioning that the spatial factors are the different distances at the start of the observed period, and the socioeconomic and climate factors are the increments of the same period which have been rastorized to county-level as depicted before. The optimal discretization method was used to discretize the meteorological and spatial factors, and Z-Score was applied to normalize the socioeconomic factors, which were all taken as independent variables. Temperature, meat product, per capita net income of farmers and herdsmen and sheep number, elevation, distance from road and distance from cropland are the important factors driving spatial differentiation of sand-control. In general, low elevation, increases in temperature, and increase in per capita net income of farmers and herdsmen and proximity to roads created a favorable condition for sand-control activity. It is worth mentioning that prior to 1999, increases in meat product and sheep number constrained sand-control, or rather, led to desertification, but after 1999 they have become favorable factors as controlled deserts have served for cropping and forage production. Prior to 2010, the closer to cropland, the higher the possibility of sandy land to be controlled; while it was the opposite after 2010, as combating desertification activity gradually moved into the heart areas of deserts.
Based on the regression coefficients and odds ratio (OR), it is possible to distinguish the importance of the socio-environmental factors (independent variables) in the sand-control or desertification event. As seen in Tables 10 and 12-14 the important factors causing the spatial variability of desertification include temperature, average wind speed, precipitation, per capita net income of farmers and herdsmen, GDP, total sown area and year-end number of big livestock. The areas far from cities and roads with low elevation and slope seem more susceptible to desertification than those close to cities and roads because the latter is easier to be managed and controlled. The increase in GDP, net income of farmers and herdsmen, the year-end number of big livestock and the total sown area for grain production are the driving forces of desertification. Additionally, the decrease in temperatures and reduction in water availability and increase in wind speeds increase the probability of desertification.  (14) 1999-2010 16 ) 0.567 (15) 2010-2020 0.723 (16) 1991-2020 Desertification 1991-1999 (18) 1999-2010 (19) 2010-2020 (20) 1991-2020 Note: P 1 , P 2 , P 3 and P 4 are the probability of sand-control occurrence respectively in 1991-1999, 1999-2010, 2010-2020 and 1991-2020, while P 5 , P 6 , P 7 , and P 8 are the probability of desertification appearance respectively in 1991-1999, 1999-2010, 2010-2020, and 1991-2020. Of the corresponding period, ∆X 1 is the difference in total sown area, ∆X 2 the difference in meat product, ∆X 3 the difference in sheep number, ∆X 4 the difference of year-end big livestock, ∆X 5 the difference of per capita net income of farmers and herdsmen, ∆X 6 the difference of gross output of farming, forestry, animal husbandry and fishery, ∆X 7 the difference of gross domestic product (GDP), ∆X 8 the difference of precipitation, ∆X 9 the difference of temperature, ∆X 10 the difference of maximum wind speed, ∆X 11 the difference of average wind speed, ∆X 12 the difference of sunshine duration, X 13 the distance from road, X 14 the distance from city, X 15 the distance from water resource, X 16 the distance from cropland, X 17 the elevation, X 18 the aspect, and X 19 the slope. Note: (1) β denotes the regression coefficient and SE standard error; Wals means Wald test, and Sig significance; OR is odds ratio. (2) Factors or variables whose variable inflation factor (VIF) was largely less than 2.5 were used for logistic regression modeling, and those with VIF > 10 were removed prior to modeling. Note: (1) β denotes the regression coefficient and SE standard error; Wals means Wald test, and Sig significance; OR is odds ratio. (2) Factors or variables whose VIF was largely less than 2.5 were used for logistic regression modeling, and those with VIF > 10 were removed prior to modeling.

Spatiotemporal Variability of the Sand-Control Effectiveness
As demonstrated above, the sand-control exhibits obvious spatiotemporal variability. This is probably because of the difference in regional implementation of the ecological restoration projects or environmental protection policies [8,33,47,48]. In 1999, the Grain for Green Program (GGP), or the Conversion of Cropland to Forest and Grassland Program, was promulgated by the Central Government and first implemented in Shaanxi, Ningxia and Inner Mongolia while the third phase of the Three-North Shelter Forest Program (TNSFP) was reaching its end. In addition, the local governments, e.g., the banners Otogqian and Otog, issued the "Grazing ban and rotation policy with a subsidy system" in 2001 [8]. At the same time, a large number of families or household companies started to devote themselves to the sand-control activities. The fourth and fifth phase of the TNSFP, which were focused on desertification prevention and control, was achieved respectively in 2010 and 2020. Moreover, the Central Government invested another 13.6 billion yuan in 2011 to eight main pastoral areas in Northern China for implementation of a protection policy with subsidy and incentive mechanisms, including subsidies for the grazing ban, incentives for balancing forage productivity and livestock, and for good performance of herdsmen (http://www.moa.gov.cn, accessed on 8 June 2021) [49]. These programs, policies and spontaneous sand-control activities have played a positive role in desertification reversal since the 1990s, but with a spatial variability in space as there was a difference in implementation progress of these initiatives. For this reason, the proportion of sand-control has increased significantly, despite of its spatial variability. Actually, a large piece of the degraded grassland in the banners Uxin, Otog and Otogqian has been restored.

Driving Forces for Desertification and Sand-Control
Some authors consider that climate is the main driving force of desertification [50], while others believe that human activities exert more impact on the latter [18]. Our research indicates that both human activities and climate condition including animal husbandry (i.e., increase in big livestock and, in particular, sheep number), GDP growth, variation in wind speed and temperature (Tables 10, 12 and 14) have played an important role in this process. More concretely, linear regression analysis reveals that desertification is positively correlated with the maximum and average wind speed and negatively correlated with the area of water (Table 10).
These analyses demonstrate that variation in strong wind speed, coal mining, grazing and reduction in water availability may have driven the occurrence of desertification.
As demonstrated in Table 14, logistic regression modeling uncovers that in different periods, the probability of desertification in space is different and the roles of the socioenvironmental factors are also different. However, in the past 30 years, the contribution of ∆X 8 , local reduction in precipitation, ∆X 11 , increase in wind speed, and X 13 , distance to roads, have played a positive role, while X 19 , slope, takes a negative part in the desertification process, meaning that steep slope does not favor desertification. Therefore, desertification is associated with both climate condition and human activity, and this finding confirms the correctness of the definition of "desertification" by UNCCD in 1994.
Increase in precipitation is conducive to vegetation growth, thus promoting sandcontrol and vegetation recovery. It is worth mentioning that previous studies have shown that over-reclamation and over-stock will inevitably exacerbate desertification [2,8,30], whereas our study reveals that since 1999, the increase in meat product including pork, beef and mutton is a favorable sand-control indicator, and since 2010, the increase of per capita income of farmers and herdsmen has been also conducive for combating desertification. This is because in the frame of the TNSFP, GGP (1999), the Ecological Grassland Protection Project (2011) and other national ecological projects, sand-control has not only promoted vegetation recovery to reverse desertification, but also produced economic value as forages such as Medicago sativa and Astragalus adsurgens pall for dry-lot feeding [8], and plantations of Elaeagnus angustifolia Linn., Lycium barbarum L., Ephedra sinica Stapf [2], etc., have brought a significant improvement in net rural income. Hence, sand-control favors both agricultural and pastoral activities.

Spatially Explicit Probability of Sand-Control and Desertification
Logistic regression modeling discovers that the area with flat terrain and low elevation is convenient for transportation and irrigation of planted trees, shrubs and herbage, and thus shows high probability of sand-control events. This is consistent with our field survey ( Figure 6). Owing to the developed transportation, high accessibility and convenience for management, the areas close to cities, roads and farmland have a high probability of being managed. Meanwhile, those far from cities, roads and farmland are prone to desertification due to less human management and costly transport; this is similar to the results of Feng et al. (2021) [33]. Water source is essential for vegetation growth, and the closer to water, the lower the cost for irrigation and the higher the possibility of being controlled. The areas with irrational agricultural activity and animal husbandry had a high probability of desertification before 1999 because profit-seeking did not go with environmental protection. As Wu et al. (2013) revealed [8], reclamation for cropping and abandonment after 2-3 years of cultivation, land was left with erosion by both wind and water, and herdsmen preferred to breed the maximum possible number of sheep and big livestock to compete for the public grassland resource to maximize their economic profit, which had been discussed by Hardin in his article "The Tragedy of the Commons" [51]. However, since the large-scale implementation of policies, such as rotational grazing and the grazing-ban policy with a subsidy mechanism in 2001, the probability of successful sand-control in these areas has been significantly increased, in particular after 2010. This is in agreement with our field survey in July 2021 (Figure 7), which shows the difference between the rotational grazing areas and the grazing-ban area.
promoted vegetation recovery to reverse desertification, but also produced economi value as forages such as Medicago sativa and Astragalus adsurgens pall for dry-lot feedin [8], and plantations of Elaeagnus angustifolia Linn., Lycium barbarum L., Ephedra sinica Stap [2], etc., have brought a significant improvement in net rural income. Hence, sand-contro favors both agricultural and pastoral activities.

Spatially Explicit Probability of Sand-Control and Desertification
Logistic regression modeling discovers that the area with flat terrain and low eleva tion is convenient for transportation and irrigation of planted trees, shrubs and herbage and thus shows high probability of sand-control events. This is consistent with our field survey ( Figure 6). Owing to the developed transportation, high accessibility and conven ience for management, the areas close to cities, roads and farmland have a high probabilit of being managed. Meanwhile, those far from cities, roads and farmland are prone to des ertification due to less human management and costly transport; this is similar to the re sults of Feng et al. (2021) [33]. Water source is essential for vegetation growth, and th closer to water, the lower the cost for irrigation and the higher the possibility of bein controlled. The areas with irrational agricultural activity and animal husbandry had high probability of desertification before 1999 because profit-seeking did not go with en vironmental protection. As Wu et al. (2013) revealed [8], reclamation for cropping and abandonment after 2-3 years of cultivation, land was left with erosion by both wind and water, and herdsmen preferred to breed the maximum possible number of sheep and bi livestock to compete for the public grassland resource to maximize their economic profi which had been discussed by Hardin in his article "The Tragedy of the Commons" [51 However, since the large-scale implementation of policies, such as rotational grazing and the grazing-ban policy with a subsidy mechanism in 2001, the probability of successfu sand-control in these areas has been significantly increased, in particular after 2010. Thi is in agreement with our field survey in July 2021 (Figure 7), which shows the differenc between the rotational grazing areas and the grazing-ban area.

Conclusions
This paper presents an integrated remote sensing-based study on the effectiveness of sand-control and desertification in the Mu Us Desert, and spatially explicit models aiding to explain spatial determinants of sand-control and desertification were developed. We found that a total of 8712.23 km 2 of sandy land (63.05% of the study area) has been controlled and converted into shrublands and grasslands with different greenness, where GDVI has increased by 0.3509 on average from 1991 to 2020. This activity has brought not only recovery of vegetation cover, but also growth in rural income and meat products, allowing an improvement of the livelihood of herdsmen, and thus created a win-win situation. This is a great success of combating desertification by local people and the Chinese governments, and this experience deserves to be extended to other desert control in Northern China and even to other dryland countries across the world. However, despite the decrease in desert area, desertification is still taking place locally and more than 5000 km 2 of deserts remain unmanaged in the study area.
Multiple linear regression models illustrate that the rural per capita income, wind speed and water source play a role in desertification process, and logistic regression models reveal that the local reduction in precipitation, increase in wind speed and the distance to roads are the key factors influencing the desertification in the past three decades. Anyway, desertification is associated with both climate and human activity, such as wind speed, precipitation, water availability, distance to roads, rural per capita income and animal husbandry. The climatic part, such as variation in wind speed and precipitation. can be hardly managed, but it is possible to control animal husbandry (e.g., rational sheep number and livestock), and to make relatively remote areas accessible by road construction and available water. In this way, existing desertification may be reversed and future desertification could be avoided. The findings of this study may provide advice for local decision-makers in taking measures to conduct sand-control of the next step in Mu Us or elsewhere.
Different from the existing MEDALUS and DesertWatch projects, which were aimed at desertification vulnerability assessment and scenario modeling, this study attempted to reveal the desertification mechanism by two kinds of regression modeling, and achieved a quantitative desertification assessment with a deepened understanding on its driving forces. We can say that this study is an extension or a complement of the two mentioned projects in the continental climate dryland environment. Probably, a combination of these two approaches may lead to a more comprehensive research on desertification. This will be the topic of our next research in the Ordos region.

Conclusions
This paper presents an integrated remote sensing-based study on the effectiveness of sand-control and desertification in the Mu Us Desert, and spatially explicit models aiding to explain spatial determinants of sand-control and desertification were developed. We found that a total of 8712.23 km 2 of sandy land (63.05% of the study area) has been controlled and converted into shrublands and grasslands with different greenness, where GDVI has increased by 0.3509 on average from 1991 to 2020. This activity has brought not only recovery of vegetation cover, but also growth in rural income and meat products, allowing an improvement of the livelihood of herdsmen, and thus created a win-win situation. This is a great success of combating desertification by local people and the Chinese governments, and this experience deserves to be extended to other desert control in Northern China and even to other dryland countries across the world. However, despite the decrease in desert area, desertification is still taking place locally and more than 5000 km 2 of deserts remain unmanaged in the study area.
Multiple linear regression models illustrate that the rural per capita income, wind speed and water source play a role in desertification process, and logistic regression models reveal that the local reduction in precipitation, increase in wind speed and the distance to roads are the key factors influencing the desertification in the past three decades. Anyway, desertification is associated with both climate and human activity, such as wind speed, precipitation, water availability, distance to roads, rural per capita income and animal husbandry. The climatic part, such as variation in wind speed and precipitation. can be hardly managed, but it is possible to control animal husbandry (e.g., rational sheep number and livestock), and to make relatively remote areas accessible by road construction and available water. In this way, existing desertification may be reversed and future desertification could be avoided. The findings of this study may provide advice for local decision-makers in taking measures to conduct sand-control of the next step in Mu Us or elsewhere.
Different from the existing MEDALUS and DesertWatch projects, which were aimed at desertification vulnerability assessment and scenario modeling, this study attempted to reveal the desertification mechanism by two kinds of regression modeling, and achieved a quantitative desertification assessment with a deepened understanding on its driving forces. We can say that this study is an extension or a complement of the two mentioned projects in the continental climate dryland environment. Probably, a combination of these two approaches may lead to a more comprehensive research on desertification. This will be the topic of our next research in the Ordos region. Data Availability Statement: Data used for this research will be available upon reasonable request to the corresponding author.